Complex and unexpected dynamics in simple genetic regulatory networks

Peer reviewedPublisher PD

External Stimuli Mediate Collective Rhythms: Artificial Control Strategies

The artificial intervention of biological rhythms remains an exciting challenge. Here, we proposed artificial control strategies that were developed to mediate the collective rhythms emerging in multicellular structures. Based on noisy repressilators and by injecting a periodic control amount to the extracellular medium, we introduced two typical kinds of control models. In one, there are information exchanges among cells, where signaling molecules receive the injected stimulus that freely diffuses toward/from the intercellular medium. In the other, there is no information exchange among cells, but signaling molecules also receive the stimulus that directionally diffuses into each cell from the common environment. We uncovered physical mechanisms for how the stimulus induces, enhances or ruins collective rhythms. We found that only when the extrinsic period is close to an integer multiplicity of the averaged intrinsic period can the collective behaviors be induced/enhanced; otherwise, the stimulus possibly ruins the achieved collective behaviors. Such entrainment properties of these oscillators to external signals would be exploited by realistic living cells to sense external signals. Our results not only provide a new perspective to the understanding of the interplays between extrinsic stimuli and intrinsic physiological rhythms, but also would lead to the development of medical therapies or devices

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

Robust Design of Biological Circuits: Evolutionary Systems Biology Approach

Artificial gene circuits have been proposed to be embedded into microbial cells that function as switches, timers, oscillators, and the Boolean logic gates. Building more complex systems from these basic gene circuit components is one key advance for biologic circuit design and synthetic biology. However, the behavior of bioengineered gene circuits remains unstable and uncertain. In this study, a nonlinear stochastic system is proposed to model the biological systems with intrinsic parameter fluctuations and environmental molecular noise from the cellular context in the host cell. Based on evolutionary systems biology algorithm, the design parameters of target gene circuits can evolve to specific values in order to robustly track a desired biologic function in spite of intrinsic and environmental noise. The fitness function is selected to be inversely proportional to the tracking error so that the evolutionary biological circuit can achieve the optimal tracking mimicking the evolutionary process of a gene circuit. Finally, several design examples are given in silico with the Monte Carlo simulation to illustrate the design procedure and to confirm the robust performance of the proposed design method. The result shows that the designed gene circuits can robustly track desired behaviors with minimal errors even with nontrivial intrinsic and external noise

A systems engineering approach to model, tune and test synthetic gene circuits

La biología sintética se define como la ingeniería de la biología: el (re)diseño y construcción de nuevas partes, dispositivos y sistemas biológicos para realizar nuevas funciones con fines útiles, que se basan en principios elucidados de la biología y la ingeniería. Para facilitar la construcción rápida, reproducible y predecible de estos sistemas biológicos a partir de conjuntos de componentes es necesario desarrollar nuevos métodos y herramientas. La tesis plantea la optimización multiobjetivo como el marco adecuado para tratar los problemas comunes que surgen en el diseño racional y el ajuste óptimo de los circuitos genéticos sintéticos. Utilizando un enfoque clásico de ingeniería de sistemas, la tesis se centra principalmente en: i) el modelado de circuitos genéticos sintéticos basado en los primeros principios, ii) la estimación de parámetros de modelos a partir de datos experimentales y iii) el ajuste basado en modelos para lograr el desempeño deseado de los circuitos. Se han utilizado dos circuitos genéticos sintéticos de diferente naturaleza y con diferentes objetivos y problemas: un circuito de realimentación de tipo 1 incoherente (I1-FFL) que exhibe la importante propiedad biológica de adaptación, y un circuito de detección de quorum sensing y realimentación (QS/Fb) que comprende dos bucles de realimentación entrelazados -uno intracelular y uno basado en la comunicación de célula a célula- diseñado para regular el nivel medio de expresión de una proteína de interés mientras se minimiza su varianza a través de la población de células. Ambos circuitos han sido analizados in silico e implementados in vivo. En ambos casos, se han desarrollado modelos de estos circuitos basado en primeros principios. Se presta especial atención a ilustrar cómo obtener modelos de orden reducido susceptibles de estimación de parámetros, pero manteniendo el significado biológico. La estimación de los parámetros del modelo a partir de los datos experimentales se considera en diferentes escenarios, tanto utilizando modelos determinísticos como estocásticos. Para el circuito I1-FFL se consideran modelos determinísticos. Aquí, la tesis plantea la utilización de modelos locales utilizando la optimización multiobjetivo para realizar la estimación de parámetros del modelo bajo escenarios con estructura de modelo incompleta. Para el circuito QS/Fb, una estructura controlada por realimentación, el problema tratado es la falta de excitabilidad de las señales. La tesis propone una metodología de estimación en dos etapas utilizando modelos estocásticos. La metodología permite utilizar datos de curso temporal promediados de la población y mediciones de distribución en estado estacionario para una sola célula. El ajuste de circuitos basado en modelos para lograr un desempeño deseado también se aborda mediante la optimización multiobjetivo. Para el circuito QS/Fb se realiza un análisis estocástico completo. La tesis aborda cómo tener en cuenta correctamente tanto el ruido intrínseco como el extrínseco, las dos principales fuentes de ruido en los circuitos genéticos. Se analiza el equilibrio entre ambas fuentes de ruido y el papel que desempeñan en el bucle de realimentación intracelular, y en la realimentación extracelular de toda la población. La principal conclusión es que la compleja interacción entre ambos canales de realimentación obliga al uso de la optimización multiobjetivo para el adecuado ajuste del circuito. En esta tesis además del uso adecuado de herramientas de optimización multiobjetivo, la principal preocupación es cómo derivar directrices para el ajuste in silico de parámetros de circuitos que puedan aplicarse de forma realista in vivo en un laboratorio estándar. Como alternativa al análisis de sensibilidad de parámetros clásico, la tesis propone el uso de técnicas de clustering a lo largo de los frentes de Pareto, relacionando el comprLa biologia sintètica es defineix com l'enginyeria de la biologia: el (re) disseny i construcció de noves parts, dispositius i sistemes biològics per a realitzar noves funcions útils que es basen a principis elucidats de la biologia i l'enginyeria. Per facilitar la construcció ràpida, reproduïble i predictible de aquests sistemes biològics a partir de conjunts de components és necessari desenvolupar nous mètodes i eines. La tesi planteja la optimització multiobjectiu com el marc adequat per a tractar els problemes comuns que apareixen en el disseny racional i l' ajust òptim dels circuits genètics sintètics. Utilitzant un enfocament clàssic d'enginyeria de sistemes, la tesi es centra principalment en: i) el modelatge de circuits genètics sintètics basat en primers principis, ii) l' estimació de paràmetres de models a partir de dades experimentals i iii) l' ajust basat en models per aconseguir el rendiment desitjat dels circuits. S'han utilitzat dos circuits genètics sintètics de diferent naturalesa i amb diferents objectius i problemes: un circuit de prealimentació de tipus 1 incoherent (I1-FFL) que exhibeix la important propietat biològica d'adaptació, i un circuit de quorum sensing i realimentació (QS/Fb) que comprèn dos bucles de realimentació entrellaçats -un intracel·lular i un basat en la comunicació de cèl·lula a cèl·lula- dis-senyat per regular el nivell mitjà d'expressió normal d'una proteïna d'interès mentre es minimitza la seua variació al llarg de la població de cèl·lules. Els dos circuits han estat analitzats in silico i implementats in vivo. En tots dos casos, s'han desenvolupat models basats en primers principis d'aquests circuits. Després es presta especial atenció a delinear com obtenir models d'ordre reduït susceptibles de estimació de paràmetres, però mantenint el significat biològic. L' estimació dels paràmetres del model a partir de les dades experimentals es considera en diferents escenaris, tant utilitzant models determinístics com estocàstics. Per al circuit I1-FFL es consideren models determinístics. La tesi planteja la utilització de models locals utilitzant la optimització multiobjectiu per realitzar l'estimació de parametres del model sota escenaris amb estructura de model incompleta (dinàmica no modelada). Per al circuit de QS/Fb, una estructura controlada per realimentació, el problema tractat és la manca d'excitabilitat dels senyals. La tesi proposa una metodologia de estimació en dues etapes utilitzant models estocàstics. La metodologia permet utilitzar dades de curs temporal promediats de la població i mesures de distribució en estat estacionari d'una sola una cèl·lula. L' ajust de circuits basat en models per aconseguir el rendiment desitjat dels circuits també s' aborda mitjançant la optimització multiobjectiu. Per al circuit QS/Fb, es fa un anàlisi estocàstic complet. La tesi aborda com tenir en compte correctament tant el soroll intrínsec com l' extrínsec, les dues principals fonts de soroll en els circuits genètics sintètics. S' analitza l'equilibri entre dues fonts de soroll i el paper que exerceixen en el bucle de realimentació intracel·lular, les i en la realimentació extracel·lular de tota la població. La principal conclusió es que la complexa interacció entre els dos canals de realimentació fa necessari l' ús de la optimització multiobjectiu per al adequat ajust del circuit. En aquesta tesi, a més de l'ús adequat d'eines d'optimització multiobjectiu, la principal preocupació és com derivar directives per al ajust in silico de paràmetres de circuits que puguin aplicar-se de forma realista en viu en un laboratori estàndard. Així, com a alternativa a l'anàlisi de sensibilitat de paràmetres clàssic, la tesi proposa l'ús de l' tècniques de l'agrupació al llarg dels fronts de Pareto, relacionant el compromís de dessempeny amb les regions en l'espai d'paràmetres.Synthetic biology is defined as the engineering of biology: the deliberate (re)design and construction of novel biological and biologically based parts, devices and systems to perform new functions for useful purposes, that draws on principles elucidated from biology and engineering. Methods and tools are needed to facilitate fast, reproducible and predictable construction of biological systems from sets of biological components. This thesis raises multi-objective optimization as the proper framework to deal with common problems arising in rational design and optimal tuning of synthetic gene circuits. Using a classical systems engineering approach, the thesis mainly addresses: i) synthetic gene circuit modeling based on first principles, ii) model parameters estimation from experimental data and iii) model-based tuning to achieve desired circuit performance. Two gene synthetic circuits of different nature and with different goals and inherent problems have been used throughout the thesis: an Incoherent type 1 feedforward circuit (I1-FFL) that exhibits the important biological property of adaptation, and a Quorum sensing/Feedback circuit (QS/Fb) comprising two intertwined feedback loops -an intracellular one and a cell-to-cell communication-based one-- designed to regulate the mean expression level of a protein of interest while minimizing its variance across the population of cells. Both circuits have been analyzed in silico and implemented in vivo. In both cases, circuit modeling based on first principles has been carried out. Then, special attention is paid to illustrate how to obtain reduced order models amenable for parameters estimation yet keeping biological significance. Model parameters estimation from experimental data is considered in different scenarios, both using deterministic and stochastic models. For the I1-FFL circuit, deterministic models are considered. In this case, the thesis raises ensemble modeling using multi-objective optimization to perform model parameters estimation under scenarios with incomplete model structure (unmodeled dynamics). For the QS/Fb gene circuit, a feedback controlled structure, the lack of excitability of the signals is the problem addressed. The thesis proposes a two-stage estimation methodology using stochastic models. The methodology allows using population averaged time-course data and steady state distribution measurements at the single-cell level. Model-based circuit tuning to achieve desired circuit performance is also addressed using multi-objective optimization. First, for the QS/Fb feedback control circuit, a complete stochastic analysis is performed. Here, the thesis addresses how to correctly take into account both intrinsic and extrinsic noise, the two main sources of noise in gene synthetic circuits. The trade-off between both sources of noise, and the role played by in the intracellular single-cell feedback loop and the extracellular population-wide feedback is analyzed. The main conclusion being that the complex interplay between both feedback channels compel the use of multi-objective optimization for proper tuning of the circuit to achieve desired performance. Thus, the thesis wraps up all the previous results and uses them to address circuit tuning for desired performance. Here, besides the proper use of multi-objective optimization tools, the main concern is how to derive guidelines for circuit parameters tuning in silico that can realistically be applied in vivo in a standard laboratory. Thus, as an alternative to classical parameters sensitivity analysis, the thesis proposes the use of clustering techniques along the optimal Pareto fronts relating the performance trade-offs with regions in the circuits parameters space.This work has been partially supported by the Spanish Government (CICYT DPI2014- 55276-C5-1) and the European Union (FEDER). The author was recipient of the grant Formación de Personal Investigador by the Universitat Politècnica de València, subprogram 1 (FPI/2013-3242). She was also recipient of the competitive grants for pre-doctoral stays Erasmus Student Placement-European Programme 2015, and FPI Mobility program 2016 of the Universitat Politècnica de València. She also received the competitive grant for a pre-doctoral stay Becas de movilidad para Jóvenes Profesores e Investigadores 2016, Programa de Becas Iberoamérica of the Santander Bank.Boada Acosta, YF. (2018). A systems engineering approach to model, tune and test synthetic gene circuits [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/112725TESI

Control Theory for Synthetic Biology: Recent Advances in System Characterization, Control Design, and Controller Implementation for Synthetic Biology

Living organisms are differentiated by their genetic material-millions to billions of DNA bases encoding thousands of genes. These genes are translated into a vast array of proteins, many of which have functions that are still unknown. Previously, it was believed that simply knowing the genetic sequence of an organism would be the key to unlocking all understanding. However, as DNA sequencing technology has become affordable, it has become clear that living cells are governed by complex, multilayered networks of gene regulation that cannot be deduced from sequence alone. Synthetic biology as a field might best be characterized as a learn-by-building approach, in which scientists attempt to engineer molecular pathways that do not exist in nature. In doing so, they test the limits of both natural and engineered organisms

Multisynchronization of Chaotic Oscillators via Nonlinear Observer Approach

The goal of this work is to synchronize a class of chaotic oscillators in a master-slave scheme, under different initial conditions, considering several slaves systems. The Chen oscillator is employed as a benchmark model and a nonlinear observer is proposed to reach synchronicity between the master and the slaves’ oscillators. The proposed observer contains a proportional and integral form of a bounded function of the synchronization error in order to provide asymptotic synchronization with a satisfactory performance. Numerical experiments were carried out to show the operation of the considered methodology

Nonlinearity and stochasticity in biochemical networks

Recent advances in biology have revolutionized our understanding of living systems. For the first time, it is possible to study the behavior of individual cells. This has led to the discovery of many amazing phenomena. For example, cells have developed intelligent mechanisms for foraging, communicating, and responding to environmental changes. These diverse functions in cells are controlled through biochemical networks consisting of many different proteins and signaling molecules. These molecules interact and affect gene expression, which in turn affects protein production. This results in a complex mesh of feedback and feedforward interactions. These complex networks are generally highly nonlinear and stochastic, making them difficult to study quantitatively. Recent studies have shown that biochemical networks are also highly modular, meaning that different parts of the network do not interact strongly with each other. These modules tend to be conserved across species and serve specific biological functions. However, detect- ing modules and identifying their function tends to be a very difficult task. To overcome some of these complexities, I present an alternative modeling approach that builds quantitative models using coarse-grained biological processes. These coarse-grained models are often stochastic (probabilistic) and highly non-linear. In this thesis, I focus on modeling biochemical networks in two distinct biological systems: Dictyostelium discoideum and microRNAs. Chapters 2 and 3 focus on cellular communication in the social amoebae Dictyostelium discoideum. Using universality, I propose a stochastic nonlinear model that describes the behavior of individual cells and cellular populations. In chapter 4 I study the interaction between messenger RNAs and noncoding RNAs, using Langevin equations

Control Theory for Synthetic Biology: Recent Advances in System Characterization, Control Design, and Controller Implementation for Synthetic Biology

Living organisms are differentiated by their genetic material-millions to billions of DNA bases encoding thousands of genes. These genes are translated into a vast array of proteins, many of which have functions that are still unknown. Previously, it was believed that simply knowing the genetic sequence of an organism would be the key to unlocking all understanding. However, as DNA sequencing technology has become affordable, it has become clear that living cells are governed by complex, multilayered networks of gene regulation that cannot be deduced from sequence alone. Synthetic biology as a field might best be characterized as a learn-by-building approach, in which scientists attempt to engineer molecular pathways that do not exist in nature. In doing so, they test the limits of both natural and engineered organisms

Stochastic analysis of nonlinear dynamics and feedback control for gene regulatory networks with applications to synthetic biology

The focus of the thesis is the investigation of the generalized repressilator model (repressing genes ordered in a ring structure). Using nonlinear bifurcation analysis stable and quasi-stable periodic orbits in this genetic network are characterized and a design for a switchable and controllable genetic oscillator is proposed. The oscillator operates around a quasi-stable periodic orbit using the classical engineering idea of read-out based control. Previous genetic oscillators have been designed around stable periodic orbits, however we explore the possibility of quasi-stable periodic orbit expecting better controllability. The ring topology of the generalized repressilator model has spatio-temporal symmetries that can be understood as propagating perturbations in discrete lattices. Network topology is a universal cross-discipline transferable concept and based on it analytical conditions for the emergence of stable and quasi-stable periodic orbits are derived. Also the length and distribution of quasi-stable oscillations are obtained. The findings suggest that long-lived transient dynamics due to feedback loops can dominate gene network dynamics. Taking the stochastic nature of gene expression into account a master equation for the generalized repressilator is derived. The stochasticity is shown to influence the onset of bifurcations and quality of oscillations. Internal noise is shown to have an overall stabilizing effect on the oscillating transients emerging from the quasi-stable periodic orbits. The insights from the read-out based control scheme for the genetic oscillator lead us to the idea to implement an algorithmic controller, which would direct any genetic circuit to a desired state. The algorithm operates model-free, i.e. in principle it is applicable to any genetic network and the input information is a data matrix of measured time series from the network dynamics. The application areas for readout-based control in genetic networks range from classical tissue engineering to stem cells specification, whenever a quantitatively and temporarily targeted intervention is required