How retroactivity impacts the robustness of genetic networks

This paper studies how retroactivity impacts the robustness of gene transcription networks against parameter perturbations. By employing the linearization technique and the real stability radius, we provide comparisons of the robustness between gene transcription networks with retroactivity and ones without retroactivity. Both numerical and analytical results show that retroactivity tends to decrease such robustness. This finding in turn implies that modular genetic networks tend to be more robust against parameter perturbations.National Science Foundation (U.S.) (NSF-CCF-I058127

Mathematical models and modular composition rules for synthetic genetic circuits

One major challenge in synthetic biology is how to design genetic circuits with predictable behaviors in various biological contexts. There are two limitations to addressing this challenge in mammalian cells. First, models that can predict circuit behaviors accurately in bacteria cells cannot be directly translated to mammalian cells. Second, upon interconnection, the behavior of a module, the building block of a circuit, may be different from its behavior in a standalone setting. In this thesis, I present a bottom-up modeling framework that can be used to predict circuit behaviors in transiently transfected mammalian cells (TTMC). The first part of the framework is based on a novel bin-dependent ODE model that can describe the behavior of modules in TTMC accurately. The second part of the framework rests upon a method of modular composition that allows model-based design of circuits. The efficacies of the bin-dependent model and the method of modular composition are validated via experimental data. The effects of retroactivity, a loading effect that arises from modular composition, on circuit behaviors are also investigated

A modular approach to analyzing biological networks

This monograph addresses the decomposition of biochemical networks into functional modules that preserve their dynamic properties upon interconnection with other modules, which permits the inference of network behavior from the properties of its constituent modules. The modular decomposition method developed here also has the property that any changes in the parameters of a chemical reaction only affect the dynamics of a single module. To illustrate our results, we define and analyze a few key biological modules that arise in gene regulation, enzymatic networks, and signaling pathways. We also provide a collection of examples that demonstrate how the behavior of a biological network can be deduced from the properties of its constituent modules, based on results from control systems theory. We then use this modular decomposition method to analyze the p53 core regulation network, which plays a key role in tumor suppression in many organisms. By decomposing the network into modules, we study the evolution of the p53 core regulation network and conduct a formal analysis of the different network configurations that emerge in the evolutionary path to complexity from putative primordial organisms to more evolved vertebrates. We develop an algorithm to solve the system of equations that describe the network behavior by interconnecting the network modules systematically, as these equations are typically difficult to solve using standard numerical solvers. In the process, we qualitatively compare the distinct types of switching behaviors that each network can exhibit. We demonstrate how our novel model for the core regulation network matches experimentally observed phenomena, and contrast this with the plausible behaviors that primordial network configurations can admit. Specifically, we show that the complexity of the p53 network in humans and evolved vertebrates permits a wide range of behaviors that can bring about distinct cell fate decisions, but that this is not the case for primordial organisms

Modular analysis of signal transduction networks

Modularity, signaling networks, sytems biologyMagdeburg, Univ., Fak. für Verfahrens- und Systemtechnik, Diss., 2007von Julio Sáez RodríguezZsfassung in dt. Sprach

Optimizing enzymatic catalysts for rapid turnover of substrates with low enzyme sequestration

Enzymes are central to both metabolism and information processing in cells. In both cases, an enzyme’s ability to accelerate a reaction without being consumed in the reaction is crucial. Nevertheless, enzymes are transiently sequestered when they bind to their substrates; this sequestration limits activity and potentially compromises information processing and signal transduction. In this article, we analyse the mechanism of enzyme–substrate catalysis from the perspective of minimizing the load on the enzymes through sequestration, while maintaining at least a minimum reaction flux. In particular, we ask: which binding free energies of the enzyme–substrate and enzyme–product reaction intermediates minimize the fraction of enzymes sequestered in complexes, while sustaining a certain minimal flux? Under reasonable biophysical assumptions, we find that the optimal design will saturate the bound on the minimal flux and reflects a basic trade-off in catalytic operation. If both binding free energies are too high, there is low sequestration, but the effective progress of the reaction is hampered. If both binding free energies are too low, there is high sequestration, and the reaction flux may also be suppressed in extreme cases. The optimal binding free energies are therefore neither too high nor too low, but in fact moderate. Moreover, the optimal difference in substrate and product binding free energies, which contributes to the thermodynamic driving force of the reaction, is in general strongly constrained by the intrinsic free-energy difference between products and reactants. Both the strategies of using a negative binding free-energy difference to drive the catalyst-bound reaction forward and of using a positive binding free-energy difference to enhance detachment of the product are limited in their efficacy

Modeling genetic circuit behavior in transiently transfected mammalian cells

Binning cells by plasmid copy number is a common practice for analyzing transient transfection data. In many kinetic models of transfected cells, protein production rates are assumed to be proportional to plasmid copy number. The validity of this assumption in transiently transfected mammalian cells is not clear; models based on this assumption appear unable to reproduce experimental flow cytometry data robustly. We hypothesize that protein saturation at high plasmid copy number is a reason previous models break down and validate our hypothesis by comparing experimental data and a stochastic chemical kinetics model. The model demonstrates that there are multiple distinct physical mechanisms that can cause saturation. On the basis of these observations, we develop a novel minimal bin-dependent ODE model that assumes different parameters for protein production in cells with low versus high numbers of plasmids. Compared to a traditional Hill-function-based model, the bin-dependent model requires only one additional parameter, but fits flow cytometry input-output data for individual modules up to twice as accurately. By composing together models of individually fit modules, we use the bin-dependent model to predict the behavior of six cascades and three feed-forward circuits. The bin-dependent models are shown to provide more accurate predictions on average than corresponding (composed) Hill-function-based models and predictions of comparable accuracy to EQuIP, while still providing a minimal ODE-based model that should be easy to integrate as a subcomponent within larger differential equation circuit models. Our analysis also demonstrates that accounting for batch effects is important in developing accurate composed models.Accepted manuscrip

Reengineering Metaphysics: Modularity, Parthood, and Evolvability in Metabolic Engineering

The premise of biological modularity is an ontological claim that appears to come out of practice. We understand that the biological world is modular because we can manipulate different parts of organisms in ways that would only work if there were discrete parts that were interchangeable. This is the foundation of the BioBrick assembly method widely used in synthetic biology. It is one of a number of methods that allows practitioners to construct and reconstruct biological pathways and devices using DNA libraries of standardized parts with known functions. In this paper, we investigate how the practice of synthetic biology reconfigures biological understanding of the key concepts of modularity and evolvability. We illustrate how this practice approach takes engineering knowledge and uses it to try to understand biological organization by showing how the construction of functional parts and processes can be used in synthetic experimental evolution. We introduce a new approach within synthetic biology that uses the premise of a parts-based ontology together with that of organismal self-organization to optimize orthogonal metabolic pathways in E. coli. We then use this and other examples to help characterize semisynthetic categories of modularity, parthood, and evolvability within the discipline.\ud \ud Part of a special issue, Ontologies of Living Beings, guest-edited by A. M. Ferner and Thomas Prade

Engineered dCas9 with reduced toxicity in bacteria: implications for genetic circuit design

Large synthetic genetic circuits require the simultaneous expression of many regulators. Deactivated Cas9 (dCas9) can serve as a repressor by having a small guide RNA (sgRNA) direct it to bind a promoter. The programmability and specificity of RNA:DNA basepairing simplifies the generation of many orthogonal sgRNAs that, in theory, could serve as a large set of regulators in a circuit. However, dCas9 is toxic in many bacteria, thus limiting how high it can be expressed, and low concentrations are quickly sequestered by multiple sgRNAs. Here, we construct a non-toxic version of dCas9 by eliminating PAM (protospacer adjacent motif) binding with a R1335K mutation (dCas9*) and recovering DNA binding by fusing it to the PhlF repressor (dCas9*_PhlF). Both the 30 bp PhlF operator and 20 bp sgRNA binding site are required to repress a promoter. The larger region required for recognition mitigates toxicity in Escherichia coli, allowing up to 9600 ± 800 molecules of dCas9*_PhlF per cell before growth or morphology are impacted, as compared to 530 ± 40 molecules of dCas9. Further, PhlF multimerization leads to an increase in average cooperativity from n = 0.9 (dCas9) to 1.6 (dCas9*_PhlF). A set of 30 orthogonal sgRNA-promoter pairs are characterized as NOT gates; however, the simultaneous use of multiple sgRNAs leads to a monotonic decline in repression and after 15 are co-expressed the dynamic range is <10-fold. This work introduces a non-toxic variant of dCas9, critical for its use in applications in metabolic engineering and synthetic biology, and exposes a limitation in the number of regulators that can be used in one cell when they rely on a shared resource.United States. Defense Advanced Research Projects Agency (DARPA HR0011-15-C-0084 Living Foundries: 1000 Molecules Program