The effect of noise on dynamics and the influence of biochemical systems

Understanding a complex system requires integration and collective analysis of data from many levels of organisation. Predictive modelling of biochemical systems is particularly challenging because of the nature of data being plagued by noise operating at each and every level. Inevitably we have to decide whether we can reliably infer the structure and dynamics of biochemical systems from present data. Here we approach this problem from many fronts by analysing the interplay between deterministic and stochastic dynamics in a broad collection of biochemical models. In a classical mathematical model we first illustrate how this interplay can be described in surprisingly simple terms; we furthermore demonstrate the advantages of a statistical point of view also for more complex systems. We then investigate strategies for the integrated analysis of models characterised by different organisational levels, and trace the propagation of noise through such systems. We use this approach to uncover, for the first time, the dynamics of metabolic adaptation of a plant pathogen throughout its life cycle and discuss the ecological implications. Finally, we investigate how reliably we can infer model parameters of biochemical models. We develop a novel sensitivity/inferability analysis framework that is generally applicable to a large fraction of current mathematical models of biochemical systems. By using this framework to quantify the effect of parametric variation on system dynamics, we provide practical guidelines as to when and why certain parameters are easily estimated while others are much harder to infer. We highlight the limitations on parameter inference due to model structure and qualitative dynamical behaviour, and identify candidate elements of control in biochemical pathways most likely of being subjected to regulation

Transient and stochastic dynamics in cellular processes

This Thesis studies different cellular and cell population processes driven by non-linear and stochastic dynamics. The problems addressed here gravitate around the concepts of transient dynamics and relaxation from a perturbed to a steady state. In this regard, in all processes studied, stochastic fluctuations, either intrinsically present in or externally applied to these systems play an important and constructive role, by either driving the systems out of equilibrium, interfering with the underlying deterministic laws, or establishing suitable levels of heterogeneity. The first part of the Thesis is committed the analysis of genetically regulated transient cellular processes. Here, we analyse, from a theoretical standpoint, three genetic circuits with pulsed excitable dynamics. We show that all circuits can work in two different excitable regimes, in contrast to what was previously speculated. We also study how, in the presence of molecular noise, these excitable circuits can generate periodic polymodal pulses due to the combination of two noise induced phenomena: stabilisation of an unstable spiral point and coherence resonance. We also studied an excitable genetic mechanism for the regulation of the transcriptional fluctuations observed in some pluripotency factors in Embryonic Stem cells. In the embryo, pluripotency is a transient cellular state and the exit of cells from it seems to be associated with transcriptional fluctuations. In regard to pluripotency control, we also propose a novel mechanism based on the post-translational regulation of a small set of four pluripotency factors. We have validated the theoretical model, based on the formation of binary complexes among these factors, with quantitative experimental data at the single-cell level. The model suggests that the pluripotency state does not depend on the cellular levels of a single factor, but rather on the equilibrium of correlations between the different proteins. In addition, the model is able to anticipate the phenotype of several mutant cell types and suggests that the regulatory function of the protein interactions is to buffer the transcriptional activity of Oc4, a key pluripotency factor. In the second part of the Thesis we studied the behaviour of a computational cell signalling network of the human fibroblast in the presence of external fluctuations and signals. The results obtained here indicate that the network responds in a nontrivial manner to background chatter, both intrinsically and in the presence of external periodic signals. We show that these responses are consequence of the rerouting of the signal to different network information-transmission paths that emerge as noise is modulated. Finally, we also study the cell population dynamics during the formation of microbial biofilms, wrinkled pellicles of bacteria glued by an extracellular matrix that are one of the simplest cases of self-organised multicellular structures. In this Thesis we develop a spatiotemporal model of cellular growth and death that accounts for the experimentally observed patterns of massive bacterial death that precede wrinkle formation in biofilms. These localised patterns focus mechanical forces during biofilm expansion and trigger the formation of the characteristic ridges. In this sense, the proposed model suggests that the death patterns emerge from the mobility changes in bacteria due to the production of extracellular matrix and the spatially inhomogeneous cellular growth. An important prediction of the model is that matrix productions is crucial for the appearance of the patterns and, therefore for winkle formation. We have also experimentally validated validated this prediction with matrix deficient bacterial strains, which show neither death patterns nor wrinkles.En aquesta Tesi s’estudien diferents processos intracel·lulars i de poblacions cel·lulars regits per dinàmica estocàstica i no lineal. El problemes biològics tractats graviten al voltant el concepte de dinàmica transitòria i de relaxació d’un estat dinàmic pertorbat a l’estat estacionari. En aquest sentit, en tots els processos estudiats, les fluctuacions estocàstiques, presents intrínsecament o aplicades de forma externa, hi tenen un paper constructiu, ja sigui empenyent els sistemes fora de l’equilibri, interferint amb les lleis deterministes subjacents, o establint els nivells d’heterogeneïtat necessaris. La primera part de la Tesi es dedica a l’estudi de processos cel·lulars transitoris regulats genèticament. En ella analitzem des d’un punt de vista teòric tres circuits genètics de control de polsos excitables i, contràriament al que s’havia especulat anteriorment, establim que tots ells poden treballar en dos tipus de règim excitable. Analitzem també com, en presència de soroll molecular, aquests circuits excitables poden generar polsos periòdics i multimodals degut a la combinació de dos fenòmens induïts per soroll: l’estabilització estocàstica d’estats inestables i la ressonància de coherència. D’altra banda, estudiem com un mecanisme genètic excitable pot ser el responsable de regular a nivell transcripcional les fluctuacions que s’observen experimentalment en alguns factors de pluripotència en cèl·lules mare embrionàries. En l’embrió, la pluripotència és un estat cel·lular transitori i la sortida de les cèl·lules d’aquest sembla que està associada a fluctuacions transcripcionals. En relació al control de la pluripotència, presentem també un nou mecanisme basat en la regulació post-traduccional d’un petit conjunt de 4 factors de pluripotència. El model teòric proposat, basat en la formació de complexos entre els diferents factors de pluripotència, l’hem validat mitjançant experiments quantitatius en cèl·lules individuals. El model postula que l’estat de pluripotència no depèn dels nivells cel·lulars d’un únic factor, sinó d’un equilibri de correlacions entre diverses proteïnes. A més, prediu el fenotip de cèl·lules mutants i suggereix que la funció reguladora de les interaccions entre les quatre proteïnes és la d’esmorteir l’activitat transcripcional d’Oct4, un dels principals factors de pluripotència. En el segon apartat de la Tesi estudiem el comportament d’una xarxa computacional de senyalització cel·lular de fibroblast humà en presència de senyals externs fluctuants i cíclics. Els resultats obtinguts mostren que la xarxa respon de forma no trivial a les fluctuacions ambientals, fins i tot en presència d’una senyal externa. Diferents nivells de soroll permeten modular la resposta de la xarxa, mitjançant la selecció de rutes alternatives de transmissió de la informació. Finalment, estudiem la dinàmica de poblacions cel·lulars durant la formació de biofilms, pel·lícules arrugades d’aglomerats de bacteris que conformen un dels exemples més simples d’estructures multicel·lulars autoorganitzades. En aquesta Tesi presentem un model espai-temporal de creixement i mort cel·lular motivat per l’evidència experimental sobre l’aparició de patrons de mort massiva de bacteris previs a la formació de les arrugues dels biofilms. Aquests patrons localitzats concentren les forces mecàniques durant l’expansió del biofilm i inicien la formació de les arrugues característiques. En aquest sentit, el model proposat explica com es formen els patrons de mort a partir dels canvis de mobilitat dels bacteris deguts a la producció de matriu extracel·lular combinats amb un creixement espacialment heterogeni. Una important predicció del model és que la producció de matriu és un procés clau per a l’aparició dels patrons i, per tant de les arrugues. En aquest aspecte, els nostres resultats experimentals en bacteris mutants que no produeixen components essencials de la matriu, confirmen les prediccions

Computational study of resting state network dynamics

Lo scopo di questa tesi è quello di mostrare, attraverso una simulazione con il software The Virtual Brain, le più importanti proprietà della dinamica cerebrale durante il resting state, ovvero quando non si è coinvolti in nessun compito preciso e non si è sottoposti a nessuno stimolo particolare. Si comincia con lo spiegare cos’è il resting state attraverso una breve revisione storica della sua scoperta, quindi si passano in rassegna alcuni metodi sperimentali utilizzati nell’analisi dell’attività cerebrale, per poi evidenziare la differenza tra connettività strutturale e funzionale. In seguito, si riassumono brevemente i concetti dei sistemi dinamici, teoria indispensabile per capire un sistema complesso come il cervello. Nel capitolo successivo, attraverso un approccio ‘bottom-up’, si illustrano sotto il profilo biologico le principali strutture del sistema nervoso, dal neurone alla corteccia cerebrale. Tutto ciò viene spiegato anche dal punto di vista dei sistemi dinamici, illustrando il pionieristico modello di Hodgkin-Huxley e poi il concetto di dinamica di popolazione. Dopo questa prima parte preliminare si entra nel dettaglio della simulazione. Prima di tutto si danno maggiori informazioni sul software The Virtual Brain, si definisce il modello di network del resting state utilizzato nella simulazione e si descrive il ‘connettoma’ adoperato. Successivamente vengono mostrati i risultati dell’analisi svolta sui dati ricavati, dai quali si mostra come la criticità e il rumore svolgano un ruolo chiave nell'emergenza di questa attività di fondo del cervello. Questi risultati vengono poi confrontati con le più importanti e recenti ricerche in questo ambito, le quali confermano i risultati del nostro lavoro. Infine, si riportano brevemente le conseguenze che porterebbe in campo medico e clinico una piena comprensione del fenomeno del resting state e la possibilità di virtualizzare l’attività cerebrale

Multicellular Systems Biology of Development

Embryonic development depends on the precise coordination of cell fate specification, patterning and morphogenesis. Although great strides have been made in the molecular understanding of each of these processes, how their interplay governs the formation of complex tissues remains poorly understood. New techniques for experimental manipulation and image quantification enable the study of development in unprecedented detail, resulting in new hypotheses on the interactions between known components. By expressing these hypotheses in terms of rules and equations, computational modeling and simulation allows one to test their consistency against experimental data. However, new computational methods are required to represent and integrate the network of interactions between gene regulation, signaling and biomechanics that extend over the molecular, cellular and tissue scales. In this thesis, I present a framework that facilitates computational modeling of multiscale multicellular systems and apply it to investigate pancreatic development and the formation of vascular networks. This framework is based on the integration of discrete cell-based models with continuous models for intracellular regulation and intercellular signaling. Specifically, gene regulatory networks are represented by differential equations to analyze cell fate regulation; interactions and distributions of signaling molecules are modeled by reaction-diffusion systems to study pattern formation; and cell-cell interactions are represented in cell-based models to investigate morphogenetic processes. A cell-centered approach is adopted that facilitates the integration of processes across the scales and simultaneously constrains model complexity. The computational methods that are required for this modeling framework have been implemented in the software platform Morpheus. This modeling and simulation environment enables the development, execution and analysis of multi-scale models of multicellular systems. These models are represented in a new domain-specific markup language that separates the biological model from the computational methods and facilitates model storage and exchange. Together with a user-friendly graphical interface, Morpheus enables computational modeling of complex developmental processes without programming and thereby widens its accessibility for biologists. To demonstrate the applicability of the framework to problems in developmental biology, two case studies are presented that address different aspects of the interplay between cell fate specification, patterning and morphogenesis. In the first, I focus on the interplay between cell fate stability and intercellular signaling. Specifically, two studies are presented that investigate how mechanisms of cell-cell communication affect cell fate regulation and spatial patterning in the pancreatic epithelium. Using bifurcation analysis and simulations of spatially coupled differential equations, it is shown that intercellular communication results in a multistability of gene expression states that can explain the scattered spatial distribution and low cell type ratio of nascent islet cells. Moreover, model analysis shows that disruption of intercellular communication induces a transition between gene expression states that can explain observations of in vitro transdifferentiation from adult acinar cells into new islet cells. These results emphasize the role of the multicellular context in cell fate regulation during development and may be used to optimize protocols for cellular reprogramming. The second case study focuses on the feedback between patterning and morphogenesis in the context of the formation of vascular networks. Integrating a cell-based model of endothelial chemotaxis with a reaction-diffusion model representing signaling molecules and extracellular matrix, it is shown that vascular network patterns with realistic morphometry can arise when signaling factors are retained by cell-modified matrix molecules. Through the validation of this model using in vitro assays, quantitative estimates are obtained for kinetic parameters that, when used in quantitative model simulations, confirm the formation of vascular networks under measured biophysical conditions. These results demonstrate the key role of the extracellular matrix in providing spatial guidance cues, a fact that may be exploited to enhance vascularization of engineered tissues. Together, the modeling framework, software platform and case studies presented in this thesis demonstrate how cell-centered computational modeling of multi-scale and multicellular systems provide powerful tools to help disentangle the complex interplay between cell fate specification, patterning and morphogenesis during embryonic development

Extension of Generalized Modeling and Application to Problems from Cell Biology

Mathematical modeling is an important tool in improving the understanding of complex biological processes. However, mathematical models are often faced with challenges that arise due to the limited knowledge of the underlying biological processes and the high number of parameters for which exact values are unknown. The method of generalized modeling is an alternative modeling approach that aims to address these challenges by extracting information about stability and bifurcations of classes of models while making only minimal assumptions on the specific functional forms of the model. This is achieved by a direct parameterization of the Jacobian in the steady state, introducing a set of generalized parameters which have a biological interpretation. In this thesis, the method of generalized modeling is extended and applied to different problems from cell biology. In the first part, we extend the method to include also the higher derivatives at the steady state. This allows an analysis of the normal form of bifurcations and thereby a more specific description of the nearby dynamics. In models of gene-regulatory networks, it is shown that the extended method can be applied to better characterize oscillatory systems and to detect bistable dynamics. In the second part, we investigate mathematical models of bone remodeling, a process that renews the human skeleton constantly. We investigate the connection between structural properties of mathematical models and the stability of steady states in different models. We find that the dynamical system operates from a stable steady state that is situated in the vicinity of bifurcations where stability can be lost, potentially leading to diseases of bone. In the third part of this thesis, models of the MAPK signal transduction pathway are analyzed. Since mathematical models for this system include a high number of parameters, statistical methods are employed to analyze stability and bifurcations. Thereby, the parameters with a strong influence on the stability of steady states are identified. By an analysis of the bifurcation structure of the MAPK cascade, it is found that a combination of multiple layers in a cascade-like way allows for additional types of dynamic behavior such as oscillations and chaos. In summary, this thesis shows that generalized modeling is a fruitful alternative modeling approach for various types of systems in cell biology.Mathematische Modelle stellen ein wichtiges Hilfmittel zur Verbesserung des Verständnisses komplexer biologischer Prozesse dar. Sie stehen jedoch vor Schwierigkeiten, wenn wenig über die zugrundeliegende biologischen Vorgänge bekannt ist und es eine große Anzahl von Parametern gibt, deren exakten Werte unbekannt sind. Die Methode des Verallgemeinerten Modellierens ist ein alternativer Modellierungsansatz mit dem Ziel, diese Schwierigkeiten dadurch anzugehen, dass dynamische Informationen über Stabilität und Bifurkationen aus Klassen von Modellen extrahiert werden, wobei nur minimale Annahmen über die spezifischen funktionalen Formen getätigt werden. Dies wird erreicht durch eine direkte Parametrisierung der Jacobimatrix im Gleichgewichtszustand, bei der neue, verallgemeinerte Parameter eingeführt werden, die eine biologische Interpretation besitzen. In dieser Arbeit wird die Methode des Verallgemeinerten Modellierens erweitert und auf verschiedene zellbiologische Probleme angewandt. Im ersten Teil wird eine Erweiterung der Methode vorgestellt, bei der die Analyse höherer Ableitungen im Gleichgewichtszustand integriert wird. Dies erlaubt die Bestimmung der Normalform von Bifurkationen und hierdurch eine spezifischere Beschreibung der Dynamik in deren Umgebung. In Modellen für genregulatorische Netzwerke wird gezeigt, dass die so erweiterte Methode zu einer besseren Charakterisierung oszillierender Systeme sowie zur Erkennung von Bistabilität verwendet werden kann. Im zweiten Teil werden mathematische Modelle zur Knochenremodellierung untersucht, einem Prozess der das menschliche Skelett kontinuierlich erneuert. Wir untersuchen den Zusammenhang zwischen strukturellen Eigenschaften verschiedener Modelle und der Stabilität von Gleichgewichtszuständen. Wir finden, dass das dynamische System von einem stabilen Zustand operiert, in dessen Nähe Bifurkationen existieren, welche das System destabilisieren und so potentiell Knochenkranheiten verursachen können. Im dritten Teil werden Modelle für den MAPK Signaltransduktionsweg analysiert. Da mathematische Modelle für dieses System eine hohe Anzahl von Parametern beinhalten, werden statistische Methoden angewandt zur Analyse von Stabilität und Bifurkationen. Zunächst werden Parameter mit einem starken Einfluss auf die Stabilität von Gleichgewichtszuständen identifizert. Durch eine Analyse der Bifurkationsstruktur wird gezeigt, dass eine kaskadenartige Kombination mehrerer Ebenen zu zusätzliche Typen von Dynamik wie Oszillationen und Chaos führt. Zusammengefasst zeigt diese Arbeit, dass Verallgemeinertes Modellieren ein fruchtbarer alternativer Modellierungsansatz für verschiedene zellbiologische Probleme ist

Mathematical models of cellular decisions: investigating immune response and apoptosis

The main objective of this thesis is to develop and analyze mathematical models of cellular decisions. This work focuses on understanding the mechanisms involved in specific cellular processes such as immune response in the vascular system, and those involved in apoptosis, or programmed cellular death. A series of simple ordinary differential equation (ODE) models are constructed describing the macrophage response to hemoglobin:haptoglobin (Hb:Hp) complexes that may be present in vascular inflammation. The models proposed a positive feedback loop between the CD163 macrophage receptor and anti-inflammatory cytokine interleukin-10 (IL-10) and bifurcation analysis predicted the existence of a cellular phenotypic switch which was experimentally verified. Moreover, these models are extended to include the intracellular mediator heme oxygenase-1 (HO-1). Analysis of the proposed models find a positive feedback mechanism between IL-10 and HO-1. This model also predicts cellular response of heme and IL-10 stimuli. For the apoptotic (cell suicide) system, a modularized model is constructed encompassing the extrinsic and intrinsic signaling pathways. Model reduction is performed by abstracting the dynamics of complexes (oligomers) at a steady-state. This simplified model is analyzed, revealing different kinetic properties between type I and type II cells, and reduced models verify results. The second model of apoptosis proposes a novel mechanism of apoptosis activation through receptor-ligand clustering, yielding robust bistability and hysteresis. Using techniques from algebraic geometry, a model selection criterion is provided between the proposed and existing model as experimental data becomes available to verify the mechanism. The models developed throughout this thesis reveal important and relevant mechanisms specific to cellular response; specifically, interactions necessary for an organism to maintain homeostasis are identified. This work enables a deeper understanding of the biological interactions and dynamics of vascular inflammation and apoptosis. The results of these models provide predictions which may motivate further experimental work and theoretical study

Compartmental modelling of the Wnt pathway: Elucidating the role of nucleo-cytoplasmic shuttling of beta-catenin and its antagonists

The Wnt pathway plays a critical role in development and disease. The key player of the pathway is b-cat. In the nucleus, the complex formation of b-cat and TCF initiates target gene expression. Its activity is mainly regulated by retention and degradation by its antagonists APC, Axin and GSK3. Based on experimental findings, I develop and investigate compartmental models in order to analyse of the role of nucleo-cytoplasmic shuttling of these proteins in Wnt signalling. I show that the compartmental separation of b-cat and its antagonists yields an increase of the b-cat/TCF concentration