In silico evolution of signaling networks using rule-based models: bistable response dynamics

One of the ultimate goals in biology is to understand the design principles of biological systems. Such principles, if they exist, can help us better understand complex, natural biological systems and guide the engineering of de novo ones. Towards deciphering design principles, in silico evolution of biological systems with proper abstraction is a promising approach. Here, we demonstrate the application of in silico evolution combined with rule-based modelling for exploring design principles of cellular signaling networks. This application is based on a computational platform, called BioJazz, which allows in silico evolution of signaling networks with unbounded complexity. We provide a detailed introduction to BioJazz architecture and implementation and describe how it can be used to evolve and/or design signaling networks with defined dynamics. For the latter, we evolve signaling networks with switch-like response dynamics and demonstrate how BioJazz can result in new biological insights on network structures that can endow bistable response dynamics. This example also demonstrated both the power of BioJazz in evolving and designing signaling networks and its limitations at the current stage of development.Comment: 24 pages, 7 figure

Synthesis of Biological and Mathematical Methods for Gene Network Control

abstract: Synthetic biology is an emerging field which melds genetics, molecular biology, network theory, and mathematical systems to understand, build, and predict gene network behavior. As an engineering discipline, developing a mathematical understanding of the genetic circuits being studied is of fundamental importance. In this dissertation, mathematical concepts for understanding, predicting, and controlling gene transcriptional networks are presented and applied to two synthetic gene network contexts. First, this engineering approach is used to improve the function of the guide ribonucleic acid (gRNA)-targeted, dCas9-regulated transcriptional cascades through analysis and targeted modification of the RNA transcript. In so doing, a fluorescent guide RNA (fgRNA) is developed to more clearly observe gRNA dynamics and aid design. It is shown that through careful optimization, RNA Polymerase II (Pol II) driven gRNA transcripts can be strong enough to exhibit measurable cascading behavior, previously only shown in RNA Polymerase III (Pol III) circuits. Second, inherent gene expression noise is used to achieve precise fractional differentiation of a population. Mathematical methods are employed to predict and understand the observed behavior, and metrics for analyzing and quantifying similar differentiation kinetics are presented. Through careful mathematical analysis and simulation, coupled with experimental data, two methods for achieving ratio control are presented, with the optimal schema for any application being dependent on the noisiness of the system under study. Together, these studies push the boundaries of gene network control, with potential applications in stem cell differentiation, therapeutics, and bio-production.Dissertation/ThesisDoctoral Dissertation Biomedical Engineering 201

Implications of Rewiring Bacterial Quorum Sensing

Bacteria employ quorum sensing, a form of cell-cell communication, to sense changes in population density and regulate gene expression accordingly. This work investigated the rewiring of one quorum-sensing module, the lux circuit from the marine bacterium Vibrio fischeri. Steady-state experiments demonstrate that rewiring the network architecture of this module can yield graded, threshold, and bistable gene expression as predicted by a mathematical model. The experiments also show that the native lux operon is most consistent with a threshold, as opposed to a bistable, response. Each of the rewired networks yielded functional population sensors at biologically relevant conditions, suggesting that this operon is particularly robust. These findings (i) permit prediction of the behaviors of quorum-sensing operons in bacterial pathogens and (ii) facilitate forward engineering of synthetic gene circuits

Exploring design principles of cellular information processing

As a summary, this work attempts to explore and uncovered design principles of certain dynamics of cellular networks by combining evolution in silico with rule-based modelling approach. Biological systems exhibit complex dynamics, due to the complex interactions in the intra- and inter- cellular biochemical reaction networks. For instance, signalling networks are composed of many enzymes and scaffolding proteins which have combinatorial interactions. These complex systems often generate response dynamics that are essential for correct decision-makings in cells. Especially, these complex interactions are results of long term of evolutionary process. With such evolutionary complexity, systems biologists aim to decipher the structure and dynamics of signalling and regulatory networks underpinning cellular responses; synthetic biologists can use this insight to alter existing networks or engineer de novo ones. Both tasks will benefit from an understanding of which structural and dynamic features of networks can emerge from evolutionary processes, through which intermediary steps these arise, and whether they embody general design principles. As natural evolution at the level of network dynamics is difficult to study, in silico evolution of network models can provide important insights. However, current tools used for in silico evolution of network dynamics are limited to ad hoc computer simulations and models. In my PhD study, with collaborators I construct the BioJazz, an extendable, user-friendly tool for simulating the evolution of dynamic biochemical networks. Unlike previous tools for in silico evolution, BioJazz allows for evolution of cellular networks with theoretically unbounded complexity by combining rule-based modelling with an encoding of networks that is akin to a genome. BioJazz can be used to implement biologically realistic selective pressures, and allows exploration of the space of network architectures and dynamics that implement prescribed physiological functions. It is provided as an open-source tool to facilitate its further development and use. I use this tool to explore the possible biochemical designs for signalling networks displaying ultrasensitive and adaptive response dynamics. By running evolutionary simulations mimicking different biochemical scenarios, we find that enzyme sequestration emerges as a key biochemical mechanism for both dynamics. Detailed analysis of these evolved networks revealed that enzyme sequestration enables both ultrasensitive and adaptive response dynamics. I verified this proposition by designing a generic model of a signalling cycle, featuring two enzymes and a sequestering (scaffold) protein. This simple system is capable of displaying both ultrasensitive and adaptive response dynamics, even more interestingly modulating the system switching between two response dynamics through perturbing the scaffold protein. These results show that enzyme sequestration can be exploited by evolution so to generate diverse response dynamics in signalling networks. From evolutionary simulations towards ultrasensitivity, bistable dynamics emerged as an alternative solution. On one hand, inspired by such results I used the fitness function as an objective function combined with different constraints to design and optimise bistable signalling networks with completely new structure and mechanism. Studying designed bistable signalling network explicates how such bistable network can be experimentally implemented. On the other hand, from studying the evolved bistable networks allosteric enzymes catalysing futile cycles appear to be a new mechanism of bistability in signalling networks. Furthermore, one of the smallest bistable signalling motifs is derived. This motif is composed of one kinase protein with two distinct conformational states and one substrate subject to phosphorylation by the kinase and auto-dephosphorylation reactions. The sufficient and necessary condition on parameters, with which the signalling motif displays bistable response dynamics, is analytically defined. By expanding the systems with more kinases, unlimited multistability emerges with potentials of implementing complex logic gates and cell state transitions. Further exploring the discovered and natural signalling networks implies shared design patterns. Motivated by searching structural boundaries between monostationary and multistationary networks, I performed algorithmic searching of multistationary signalling networks intending to find the sufficient structural conditions for multistationarity in signalling networks

Enzyme sequestration as a tuning point in controlling response dynamics of signalling networks

Signalling networks result from combinatorial interactions among many enzymes and scaffolding proteins. These complex systems generate response dynamics that are often essential for correct decision-making in cells. Uncovering biochemical design principles that underpin such response dynamics is a prerequisite to understand evolved signalling networks and to design synthetic ones. Here, we use in silico evolution to explore the possible biochemical design space for signalling networks displaying ultrasensitive and adaptive response dynamics. By running evolutionary simulations mimicking different biochemical scenarios, we find that enzyme sequestration emerges as a key mechanism for enabling such dynamics. Inspired by these findings, and to test the role of sequestration, we design a generic, minimalist model of a signalling cycle, featuring two enzymes and a single scaffolding protein. We show that this simple system is capable of displaying both ultrasensitive and adaptive response dynamics. Furthermore, we find that tuning the concentration or kinetics of the sequestering protein can shift system dynamics between these two response types. These empirical results suggest that enzyme sequestration through scaffolding proteins is exploited by evolution to generate diverse response dynamics in signalling networks and could provide an engineering point in synthetic biology applications

The Comet Cometh: Evolving Developmental Systems

abstract: In a recent opinion piece, Denis Duboule has claimed that the increasing shift towards systems biology is driving evolutionary and developmental biology apart, and that a true reunification of these two disciplines within the framework of evolutionary developmental biology (EvoDevo) may easily take another 100 years. He identifies methodological, epistemological, and social differences as causes for this supposed separation. Our article provides a contrasting view. We argue that Duboule’s prediction is based on a one-sided understanding of systems biology as a science that is only interested in functional, not evolutionary, aspects of biological processes. Instead, we propose a research program for an evolutionary systems biology, which is based on local exploration of the configuration space in evolving developmental systems. We call this approach—which is based on reverse engineering, simulation, and mathematical analysis—the natural history of configuration space. We discuss a number of illustrative examples that demonstrate the past success of local exploration, as opposed to global mapping, in different biological contexts. We argue that this pragmatic mode of inquiry can be extended and applied to the mathematical analysis of the developmental repertoire and evolutionary potential of evolving developmental mechanisms and that evolutionary systems biology so conceived provides a pragmatic epistemological framework for the EvoDevo synthesis.The electronic version of this article is the complete one and can be found online at: https://link.springer.com/article/10.1007/s13752-015-0203-

Dynamical Models of biological networks

In der Molekularbiologie sind mathematische Modelle von regulatorischen und metabolischen Netzwerken essentiell, um von einer Betrachtung isolierter Komponenten und Interaktionen zu einer systemischen Betrachtungsweise zu kommen. Genregulatorische Systeme eignen sich besonders gut zur Modellierung, da sie experimentell leicht zugänglich und manipulierbar sind. In dieser Arbeit werden verschiedene genregulatorische Netzwerke unter Zuhilfenahme von mathematischen Modellen analysiert. Weiteres wird ein Modell einer in silico Zelle vorgestellt und diskutiert. Zunächst werden zwei zyklische genregulatorische Netzwerke - der klassische Repressilator und ein Repressilator mit zusätzlicher Autoaktivierung – im Detail mit analytischen Methoden untersucht. Um den Einfluß zufällig schwankender Molekülzahlen auf die Dynamik der beiden Systeme zu untersuchen, werden stochastische Modelle erstellt und die beiden oszillierenden Systeme verglichen. Weiteres werden mögliche Auswirkungen von Genduplikationen auf ein einfaches genregulatorisches Netzwerk untersucht. Dazu wird zunächst ein kleines Netzwerk von GATA Transkriptionsfaktoren, das eine zentrale Rolle in der Regulation des Stickstoffmetabolismus in Hefe spielt, modelliert und das Modell mit experimentellen Daten verglichen, um Parameterregionen einschränken zu können. Außerdem werden potentielle Topologien genregulatorischer Netzwerke von GATA Transkriptionsfaktoren in verwandten Fungi mittels sequenzbasierender Methoden gesucht und verglichen. Im letzten Teil der Arbeit wird MiniCellSim vorgestellt, ein Modell einer selbständigen in silico Zelle. Es erlaubt ein dynamisches System, das eine Protozelle mit einem genregulatorischen Netzwerk, einem einfachen Metabolismus und einer Zellmembran beschreibt, aus einer Sequenz abzuleiten. Nachdem alle Parameter, die zur Berechnung des dynamischen Systems benötigt werden, ohne zusätzliche Eingabe nur aus der Sequenzinformation abgeleitet werden, kann das Modell für Studien zur Evolution von genregulatorischen Netzwerken verwendet werden.In this thesis different types of gene regulatory networks are analysed using mathematical models. Further a computational framework of a novel, self-contained in silico cell model is described and discussed. At first the behaviour of two cyclic gene regulatory systems - the classical repressilator and a repressilator with additional auto-activation - are inspected in detail using analytical bifurcation analysis. To examine the behaviour under random fluctuations, stochastic versions of the systems are created. Using the analytical results sustained oscillations in the stochastic versions are obtained, and the two oscillating systems compared. In the second part of the thesis possible implications of gene duplication on a simple gene regulatory system are inspected. A model of a small network formed by GATA-type transcription factors, central in nitrogen catabolite repression in yeast, is created and validated against experimental data to obtain approximate parameter values. Further, topologies of potential gene regulatory networks and modules consisting of GATA-type transcription factors in other fungi are derived using sequence-based approaches and compared. The last part describes MiniCellSim, a model of a self-contained in silico cell. In this framework a dynamical system describing a protocell with a gene regulatory network, a simple metabolism, and a cell membrane is derived from a string representing a genome. All the relevant parameters required to compute the time evolution of the dynamical system are calculated from within the model, allowing the system to be used in studies of evolution of gene regulatory and metabolic networks

Addressing Evolutionary Questions with Synthetic Biology

Synthetic biology emerged as an engineering discipline to design and construct artificial biological systems. Synthetic biological designs aim to achieve specific biological behavior, which can be exploited for biotechnological, medical and industrial purposes. In addition, mimicking natural systems using well-characterized biological parts also provides powerful experimental systems to study evolution at the molecular and systems level. A strength of synthetic biology is to go beyond nature’s toolkit, to test alternative versions and to study a particular biological system and its phenotype in isolation and in a quantitative manner. Here, we review recent work that implemented synthetic systems, ranging from simple regulatory circuits, rewired cellular networks to artificial genomes and viruses, to study fundamental evolutionary concepts. In particular, engineering, perturbing or subjecting these synthetic systems to experimental laboratory evolution provides a mechanistic understanding on important evolutionary questions, such as: Why did particular regulatory networks topologies evolve and not others? What happens if we rewire regulatory networks? Could an expanded genetic code provide an evolutionary advantage? How important is the structure of genome and number of chromosomes? Although the field of evolutionary synthetic biology is still in its teens, further advances in synthetic biology provide exciting technologies and novel systems that promise to yield fundamental insights into evolutionary principles in the near future