Intrinsic noise profoundly alters the dynamics and steady state of morphogen-controlled bistable genetic switches

During tissue development, patterns of gene expression determine the spatial arrangement of cell types. In many cases, gradients of secreted signaling molecules - morphogens - guide this process. The continuous positional information provided by the gradient is converted into discrete cell types by the downstream transcriptional network that responds to the morphogen. A mechanism commonly used to implement a sharp transition between two adjacent cell fates is the genetic toggle switch, composed of cross-repressing transcriptional determinants. Previous analyses emphasize the steady state output of these mechanisms. Here, we explore the dynamics of the toggle switch and use exact numerical simulations of the kinetic reactions, the Chemical Langevin Equation, and Minimum Action Path theory to establish a framework for studying the effect of gene expression noise on patterning time and boundary position. This provides insight into the time scale, gene expression trajectories and directionality of stochastic switching events between cell states. Taking gene expression noise into account predicts that the final boundary position of a morphogen-induced toggle switch, although robust to changes in the details of the noise, is distinct from that of the deterministic system. Moreover, stochastic switching introduces differences in patterning time along the morphogen gradient that result in a patterning wave propagating away from the morphogen source. The velocity of this wave is influenced by noise; the wave sharpens and slows as it advances and may never reach steady state in a biologically relevant time. This could explain experimentally observed dynamics of pattern formation. Together the analysis reveals the importance of dynamical transients for understanding morphogen-driven transcriptional networks and indicates that gene expression noise can qualitatively alter developmental patterning

The stochastic behavior of a molecular switching circuit with feedback

Background: Using a statistical physics approach, we study the stochastic switching behavior of a model circuit of multisite phosphorylation and dephosphorylation with feedback. The circuit consists of a kinase and phosphatase acting on multiple sites of a substrate that, contingent on its modification state, catalyzes its own phosphorylation and, in a symmetric scenario, dephosphorylation. The symmetric case is viewed as a cartoon of conflicting feedback that could result from antagonistic pathways impinging on the state of a shared component. Results: Multisite phosphorylation is sufficient for bistable behavior under feedback even when catalysis is linear in substrate concentration, which is the case we consider. We compute the phase diagram, fluctuation spectrum and large-deviation properties related to switch memory within a statistical mechanics framework. Bistability occurs as either a first-order or second-order non-equilibrium phase transition, depending on the network symmetries and the ratio of phosphatase to kinase numbers. In the second-order case, the circuit never leaves the bistable regime upon increasing the number of substrate molecules at constant kinase to phosphatase ratio. Conclusions: The number of substrate molecules is a key parameter controlling both the onset of the bistable regime, fluctuation intensity, and the residence time in a switched state. The relevance of the concept of memory depends on the degree of switch symmetry, as memory presupposes information to be remembered, which is highest for equal residence times in the switched states. Reviewers: This article was reviewed by Artem Novozhilov (nominated by Eugene Koonin), Sergei Maslov, and Ned Wingreen.Comment: Version published in Biology Direct including reviewer comments and author responses, 28 pages, 7 figure

Complex and unexpected dynamics in simple genetic regulatory networks

Peer reviewedPublisher PD

Computational Models of the Notch Network Elucidate Mechanisms of Context-dependent Signaling

The Notch signaling pathway controls numerous cell fate decisions during development and adulthood through diverse mechanisms. Thus, whereas it functions as an oscillator during somitogenesis, it can mediate an all-or-none cell fate switch to influence pattern formation in various tissues during development. Furthermore, while in some contexts continuous Notch signaling is required, in others a transient Notch signal is sufficient to influence cell fate decisions. However, the signaling mechanisms that underlie these diverse behaviors in different cellular contexts have not been understood. Notch1 along with two downstream transcription factors hes1 and RBP-Jk forms an intricate network of positive and negative feedback loops, and we have implemented a systems biology approach to computationally study this gene regulation network. Our results indicate that the system exhibits bistability and is capable of switching states at a critical level of Notch signaling initiated by its ligand Delta in a particular range of parameter values. In this mode, transient activation of Delta is also capable of inducing prolonged high expression of Hes1, mimicking the “ON” state depending on the intensity and duration of the signal. Furthermore, this system is highly sensitive to certain model parameters and can transition from functioning as a bistable switch to an oscillator by tuning a single parameter value. This parameter, the transcriptional repression constant of hes1, can thus qualitatively govern the behavior of the signaling network. In addition, we find that the system is able to dampen and reduce the effects of biological noise that arise from stochastic effects in gene expression for systems that respond quickly to Notch signaling

Mammalian Brain As a Network of Networks

Acknowledgements AZ, SG and AL acknowledge support from the Russian Science Foundation (16-12-00077). Authors thank T. Kuznetsova for Fig. 6.Peer reviewedPublisher PD

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