Generating Effective Models and Parameters for RNA Genetic Circuits

Generative and predictive neural networks for the design of functional RNA molecules

RNA is a remarkably versatile molecule that has been engineered for applications in therapeutics, diagnostics, and in vivo information-processing systems. However, the complex relationship between the sequence and structural properties of an RNA molecule and its ability to perform specific functions often necessitates extensive experimental screening of candidate sequences. Here we present a generalized neural network architecture that utilizes the sequence and structure of RNA molecules (SANDSTORM) to inform functional predictions. We demonstrate that this approach achieves state-of-the-art performance across several distinct RNA prediction tasks, while learning interpretable abstractions of RNA secondary structure. We paired these predictive models with generative adversarial RNA design networks (GARDN), allowing the generative modelling of novel mRNA 5' untranslated regions and toehold switch riboregulators exhibiting a predetermined fitness. This approach enabled the design of novel toehold switches with a 43-fold increase in experimentally characterized dynamic range compared to those designed using classic thermodynamic algorithms. SANDSTORM and GARDN thus represent powerful new predictive and generative tools for the development of diagnostic and therapeutic RNA molecules with improved function.

Layered feedback control overcomes performance trade-off in synthetic biomolecular networks

Layered feedback is an optimization strategy in feedback control designs widely used in engineering. Control theory suggests that layering multiple feedbacks could overcome the robustness-speed performance trade-off limit. In natural biological networks, genes are often regulated in layers to adapt to environmental perturbations. It is hypothesized layering architecture could also overcome the robustness-speed performance trade-off in genetic networks. In this work, we validate this hypothesis with a synthetic biomolecular network in living E. coli cells. We start with system dynamics analysis using models of various complexities to guide the design of a layered control architecture in living cells. Experimentally, we interrogate system dynamics under three groups of perturbations. We consistently observe that the layered control improves system performance in the robustness-speed domain. This work confirms that layered control could be adopted in synthetic biomolecular networks for performance optimization. It also provides insights into understanding genetic feedback control architectures in nature.

Effects of DNA template preparation on variability in cell-free protein production

DNA templates for protein production remain an unexplored source of variability in the performance of cell-free expression (CFE) systems. To characterize this variability, we investigated the effects of two common DNA extraction methodologies, a postprocessing step and manual versus automated preparation on protein production using CFE. We assess the concentration of the DNA template, the quality of the DNA template in terms of physical damage and the quality of the DNA solution in terms of purity resulting from eight DNA preparation workflows. We measure the variance in protein titer and rate of protein production in CFE reactions associated with the biological replicate of the DNA template, the technical replicate DNA solution prepared with the same workflow and the measurement replicate of nominally identical CFE reactions. We offer practical guidance for preparing and characterizing DNA templates to achieve acceptable variability in CFE performance.

BioCRNpyler: Compiling chemical reaction networks from biomolecular parts in diverse contexts

Biochemical interactions in systems and synthetic biology are often modeled with chemical reaction networks (CRNs). CRNs provide a principled modeling environment capable of expressing a huge range of biochemical processes. In this paper, we present a software toolbox, written in Python, that compiles high-level design specifications represented using a modular library of biochemical parts, mechanisms, and contexts to CRN implementations. This compilation process offers four advantages. First, the building of the actual CRN representation is automatic and outputs Systems Biology Markup Language (SBML) models compatible with numerous simulators. Second, a library of modular biochemical components allows for different architectures and implementations of biochemical circuits to be represented succinctly with design choices propagated throughout the underlying CRN automatically. This prevents the often occurring mismatch between high-level designs and model dynamics. Third, high-level design specification can be embedded into diverse biomolecular environments, such as cell-free extracts and in vivo milieus. Finally, our software toolbox has a parameter database, which allows users to rapidly prototype large models using very few parameters which can be customized later. By using BioCRNpyler, users ranging from expert modelers to novice script-writers can easily build, manage, and explore sophisticated biochemical models using diverse biochemical implementations, environments, and modeling assumptions.

This paper describes a new software package BioCRNpyler (pronounced “Biocompiler”) designed to support rapid development and exploration of mathematical models of biochemical networks and circuits by computational biologists, systems biologists, and synthetic biologists. BioCRNpyler allows its users to generate large complex models using very few lines of code in a way that is modular. To do this, BioCRNpyler uses a powerful new representation of biochemical circuits which defines their parts, underlying biochemical mechanisms, and chemical context independently. BioCRNpyler was developed as a Python scripting language designed to be accessible to beginning users as well as easily extendable and customizable for advanced users. Ultimately, we see Biocrnpyler being used to accelerate computer automated design of biochemical circuits and model driven hypothesis generation in biology.

Best Practices for DNA Template Preparation Toward Improved Reproducibility in Cell-Free Protein Production

Performance variability is a common challenge in cell-free protein production and hinders a wider adoption of these systems for both research and biomanufacturing. While the inherent stochasticity and complexity of biology likely contributes to variability, other systematic factors may also play a role, including the source and preparation of the cell extract, the composition of the supplemental reaction buffer, the facility at which experiments are conducted, and the human operator (Cole et al. ACS Synth Biol 8:2080-2091, 2019). Variability in protein production could also arise from differences in the DNA template-specifically the amount of functional DNA added to a cell-free reaction and the quality of the DNA preparation in terms of contaminants and strand breakage. Here, we present protocols and suggest best practices optimized for DNA template preparation and quantitation for cell-free systems toward reducing variability in cell-free protein production.

Functionalizing Cell-Free Systems with CRISPR-Associated Proteins: Application to RNA-Based Circuit Engineering

Cell-free systems have become a compelling choice for the prototyping of synthetic circuits. Many robust protocols for preparing cell-free systems are now available along with toolboxes designed for a variety of applications. Thus far, the production of cell-free extracts has often been decoupled from the production of functionalized proteins. Here, we leveraged a recent protocol for producing an E. coli-based cell-free expression system with two CRISPR-associated proteins, Csy4 and dCas9, expressed prior to harvest. We found that pre-expression did not affect the resulting extract performance, and the final concentrations of the endonucleases matched the level required for synthetic circuit prototyping. We demonstrated the benefits and versatility of dCas9 and Csy4 through the use of RNA circuitry based on a combination of single guide RNAs, small transcriptional activator RNAs, and toehold switches. For instance, we show that Csy4 processing increased 4-fold the dynamic range of a previously published AND-logic gate. Additionally, blending the CRISPR-enhanced extracts enabled us to reduce leakage in a multiple inputs gate, and to extend the type of Boolean functions available for RNA-based circuits, such as NAND-logic. Finally, we reported the use of simultaneous transcriptional and translational reporters in our RNA-based circuits. In particular, the AND-gate mRNA and protein levels were able to be independently monitored in response to transcriptional and translational activators. We hope this work will facilitate the adoption of advanced processing tools for RNA-based circuit prototyping in a cell-free environment.

Effective Biophysical Modeling of Cell Free Transcription and Translation Processes

Transcription and translation are at the heart of metabolism and signal transduction. In this study, we developed an effective biophysical modeling approach to simulate transcription and translation processes. The model, composed of coupled ordinary differential equations, was tested by comparing simulations of two cell free synthetic circuits with experimental measurements generated in this study. First, we considered a simple circuit in which sigma factor 70 induced the expression of green fluorescent protein. This relatively simple case was then followed by a more complex negative feedback circuit in which two control genes were coupled to the expression of a third reporter gene, green fluorescent protein. Many of the model parameters were estimated from previous biophysical studies in the literature, while the remaining unknown model parameters for each circuit were estimated by minimizing the difference between model simulations and messenger RNA (mRNA) and protein measurements generated in this study. In particular, either parameter estimates from published studies were used directly, or characteristic values found in the literature were used to establish feasible ranges for the parameter estimation problem. In order to perform a detailed analysis of the influence of individual model parameters on the expression dynamics of each circuit, global sensitivity analysis was used. Taken together, the effective biophysical modeling approach captured the expression dynamics, including the transcription dynamics, for the two synthetic cell free circuits. While, we considered only two circuits here, this approach could potentially be extended to simulate other genetic circuits in both cell free and whole cell biomolecular applications as the equations governing the regulatory control functions are modular and easily modifiable. The model code, parameters, and analysis scripts are available for download under an MIT software license from the Varnerlab GitHub repository.

An Automated Model Test System for Systematic Development and Improvement of Gene Expression Models

Gene expression models greatly accelerate the engineering of synthetic metabolic pathways and genetic circuits by predicting sequence-function relationships and reducing trial-and-error experimentation. However, developing models with more accurate predictions remains a significant challenge. Here we present a model test system that combines advanced statistics, machine learning, and a database of 9862 characterized genetic systems to automatically quantify model accuracies, accept or reject mechanistic hypotheses, and identify areas for model improvement. We also introduce model capacity, a new information theoretic metric for correct cross-data-set comparisons. We demonstrate the model test system by comparing six models of translation initiation rate, evaluating 100 mechanistic hypotheses, and uncovering new sequence determinants that control protein expression levels. We then applied these results to develop a biophysical model of translation initiation rate with significant improvements in accuracy. Automated model test systems will dramatically accelerate the development of gene expression models, and thereby transition synthetic biology into a mature engineering discipline.

Cell-Free Systems: A Proving Ground for Rational Biodesign

Cell-free gene expression systems present an alternative approach to synthetic biology, where biological gene expression is harnessed inside non-living, in vitro biochemical reactions. Taking advantage of a plethora of recent experimental innovations, they easily overcome certain challenges for computer-aided biological design. For instance, their open nature renders all their components directly accessible, greatly facilitating model construction and validation. At the same time, these systems present their own unique difficulties, such as limited reaction lifetimes and lack of homeostasis. In this Perspective, I propose that cell-free systems are an ideal proving ground to test rational biodesign strategies, as demonstrated by a small but growing number of examples of model-guided, forward engineered cell-free biosystems. It is likely that advances gained from this approach will contribute to our efforts to more reliably and systematically engineer both cell-free as well as living cellular systems for useful applications.

The Evolution of Cell Free Biomanufacturing

Cell-free systems are a widely used research tool in systems and synthetic biology and a promising platform for manufacturing of proteins and chemicals. In the past, cell-free biology was primarily used to better understand fundamental biochemical processes. Notably, E. coli cell-free extracts were used in the 1960s to decipher the sequencing of the genetic code. Since then, the transcription and translation capabilities of cell-free systems have been repeatedly optimized to improve energy efficiency and product yield. Today, cell-free systems, in combination with the rise of synthetic biology, have taken on a new role as a promising technology for just-in-time manufacturing of therapeutically important biologics and high-value small molecules. They have also been implemented at an industrial scale for the production of antibodies and cytokines. In this review, we discuss the evolution of cell-free technologies, in particular advancements in extract preparation, cell-free protein synthesis, and cell-free metabolic engineering applications. We then conclude with a discussion of the mathematical modeling of cell-free systems. Mathematical modeling of cell-free processes could be critical to addressing performance bottlenecks and estimating the costs of cell-free manufactured products.

In vitro implementation of robust gene regulation in a synthetic biomolecular integral controller

Feedback mechanisms play a critical role in the maintenance of cell homeostasis in the presence of disturbances and uncertainties. Motivated by the need to tune the dynamics and improve the robustness of gene circuits, biological engineers have proposed various designs that mimic natural molecular feedback control mechanisms. However, practical and predictable implementations have proved challenging because of the complexity of synthesis and analysis of complex biomolecular networks. Here, we analyze and experimentally validate a synthetic biomolecular controller executed in vitro. The controller ensures that gene expression rate tracks an externally imposed reference level, and achieves this goal even in the presence of certain kinds of disturbances. Our design relies upon an analog of the well-known principle of integral feedback in control theory. We implement the controller in an Escherichia coli cell-free transcription-translation system, which allows rapid prototyping and implementation. Modeling and theory guide experimental implementation with well-defined operational predictability.

Toward a genome scale sequence specific dynamic model of cell-free protein synthesis in Escherichia coli

In this study, we developed a dynamic mathematical model of E. coli cell-free protein synthesis (CFPS). Model parameters were estimated from a dataset consisting of glucose, organic acids, energy species, amino acids, and protein product, chloramphenicol acetyltransferase (CAT) measurements. The model was successfully trained to simulate these measurements, especially those of the central carbon metabolism. We then used the trained model to evaluate the performance, e.g., the yield and rates of protein production. CAT was produced with an energy efficiency of 12%, suggesting that the process could be further optimized. Reaction group knockouts showed that protein productivity was most sensitive to the oxidative phosphorylation and glycolysis/gluconeogenesis pathways. Amino acid biosynthesis was also important for productivity, while overflow metabolism and TCA cycle affected the overall system state. In addition, translation was more important to productivity than transcription. Finally, CAT production was robust to allosteric control, as were most of the predicted metabolite concentrations; the exceptions to this were the concentrations of succinate and malate, and to a lesser extent pyruvate and acetate, which varied from the measured values when allosteric control was removed. This study is the first to use kinetic modeling to predict dynamic protein production in a cell-free E. coli system, and could provide a foundation for genome scale, dynamic modeling of cell-free E. coli protein synthesis.

•

Protein production is biphasic, powered initially by glucose and later by pyruvate.

•

Protein is produced with an energy efficiency of only 12%.

•

Protein productivity is most sensitive to oxidative phosphorylation and glycolysis.

•

Protein production is robust to allosteric control.

Cell-Free Prototyping of AND-Logic Gates Based on Heterogeneous RNA Activators

RNA-based devices controlling gene expression bear great promise for synthetic biology, as they offer many advantages such as short response times and light metabolic burden compared to protein-circuits. However, little work has been done regarding their integration to multilevel regulated circuits. In this work, we combined a variety of small transcriptional activator RNAs (STARs) and toehold switches to build highly effective AND-gates. To characterize the components and their dynamic range, we used an Escherichia coli (E. coli) cell-free transcription-translation (TX-TL) system dispensed via nanoliter droplets. We analyzed a prototype gate in vitro as well as in silico, employing parametrized ordinary differential equations (ODEs), for which parameters were inferred via parallel tempering, a Markov chain Monte Carlo (MCMC) method. On the basis of this analysis, we created nine additional AND-gates and tested them in vitro. The functionality of the gates was found to be highly dependent on the concentration of the activating RNA for either the STAR or the toehold switch. All gates were successfully implemented in vivo, offering a dynamic range comparable to the level of protein circuits. This study shows the potential of a rapid prototyping approach for RNA circuit design, using cell-free systems in combination with a model prediction.

Point-of-care biomarker quantification enabled by sample-specific calibration

Easy-to-perform, relatively inexpensive blood diagnostics have transformed at-home healthcare for some patients, but they require analytical equipment and are not easily adapted to measuring other biomarkers. The requirement for reliable quantification in complex sample types (such as blood) has been a critical roadblock in developing and deploying inexpensive, minimal-equipment diagnostics. Here, we developed a platform for inexpensive, easy-to-use diagnostics that uses cell-free expression to generate colored readouts that are visible to the naked eye, yet quantitative and robust to the interference effects seen in complex samples. We achieved this via a parallelized calibration scheme that uses the patient sample to generate custom reference curves. We used this approach to quantify a clinically relevant micronutrient and to quantify nucleic acids, demonstrating a generalizable platform for low-cost quantitative diagnostics.

Cell-Free Synthetic Biology Platform for Engineering Synthetic Biological Circuits and Systems

Synthetic biology brings engineering disciplines to create novel biological systems for biomedical and technological applications. The substantial growth of the synthetic biology field in the past decade is poised to transform biotechnology and medicine. To streamline design processes and facilitate debugging of complex synthetic circuits, cell-free synthetic biology approaches has reached broad research communities both in academia and industry. By recapitulating gene expression systems in vitro, cell-free expression systems offer flexibility to explore beyond the confines of living cells and allow networking of synthetic and natural systems. Here, we review the capabilities of the current cell-free platforms, focusing on nucleic acid-based molecular programs and circuit construction. We survey the recent developments including cell-free transcription-translation platforms, DNA nanostructures and circuits, and novel classes of riboregulators. The links to mathematical models and the prospects of cell-free synthetic biology platforms will also be discussed.

Evolving methods for rational de novo design of functional RNA molecules

Artificial RNA molecules with novel functionality have many applications in synthetic biology, pharmacy and white biotechnology. The de novo design of such devices using computational methods and prediction tools is a resource-efficient alternative to experimental screening and selection pipelines. In this review, we describe methods common to many such computational approaches, thoroughly dissect these methods and highlight open questions for the individual steps. Initially, it is essential to investigate the biological target system, the regulatory mechanism that will be exploited, as well as the desired components in order to define design objectives. Subsequent computational design is needed to combine the selected components and to obtain novel functionality. This process can usually be split into constrained sequence sampling, the formulation of an optimization problem and an in silico analysis to narrow down the number of candidates with respect to secondary goals. Finally, experimental analysis is important to check whether the defined design objectives are indeed met in the target environment and detailed characterization experiments should be performed to improve the mechanistic models and detect missing design requirements.

Sequence Specific Modeling of E. coli Cell-Free Protein Synthesis

Cell-free protein synthesis (CFPS) is a widely used research tool in systems and synthetic biology. However, if CFPS is to become a mainstream technology for applications such as point of care manufacturing, we must understand the performance limits and costs of these systems. Toward this question, we used sequence specific constraint based modeling to evaluate the performance of E. coli cell-free protein synthesis. A core E. coli metabolic network, describing glycolysis, the pentose phosphate pathway, energy metabolism, amino acid biosynthesis, and degradation was augmented with sequence specific descriptions of transcription and translation and effective models of promoter function. Model parameters were largely taken from literature; thus the constraint based approach coupled the transcription and translation of the protein product, and the regulation of gene expression, with the availability of metabolic resources using only a limited number of adjustable model parameters. We tested this approach by simulating the expression of two model proteins: chloramphenicol acetyltransferase and dual emission green fluorescent protein, for which we have data sets; we then expanded the simulations to a range of additional proteins. Protein expression simulations were consistent with measurements for a variety of cases. The constraint based simulations confirmed that oxidative phosphorylation was active in the CAT cell-free extract, as without it there was no feasible solution within the experimental constraints of the system. We then compared the metabolism of theoretically optimal and experimentally constrained CFPS reactions, and developed parameter free correlations which could be used to estimate productivity as a function of carbon number and promoter type. Lastly, global sensitivity analysis identified the key metabolic processes that controlled CFPS productivity and energy efficiency. In summary, sequence specific constraint based modeling of CFPS offered a novel means to a priori estimate the performance of a cell-free system, using only a limited number of adjustable parameters. While we modeled the production of a single protein in this study, the approach could easily be extended to multiprotein synthetic circuits, RNA circuits, or the cell-free production of small molecule products.

Engineering a Functional Small RNA Negative Autoregulation Network with Model-Guided Design

RNA regulators are powerful components of the synthetic biology toolbox. Here, we expand the repertoire of synthetic gene networks built from these regulators by constructing a transcriptional negative autoregulation (NAR) network out of small RNAs (sRNAs). NAR network motifs are core motifs of natural genetic networks, and are known for reducing network response time and steady state signal. Here we use cell-free transcription-translation (TX-TL) reactions and a computational model to design and prototype sRNA NAR constructs. Using parameter sensitivity analysis, we design a simple set of experiments that allow us to accurately predict NAR function in TX-TL. We transfer successful network designs into Escherichia coli and show that our sRNA transcriptional network reduces both network response time and steady-state gene expression. This work broadens our ability to construct increasingly sophisticated RNA genetic networks with predictable function.

Quantitative Characterization of Translational Riboregulators Using an in Vitro Transcription-Translation System

Riboregulators are short RNA sequences that, upon binding to a ligand, change their secondary structure and influence the expression rate of a downstream gene. They constitute an attractive alternative to transcription factors for building synthetic gene regulatory networks because they can be engineered de novo. However, riboregulators are generally designed in silico and tested in vivo, which provides little quantitative information about their performances, thus hindering the improvement of design algorithms. Here we show that a cell-free transcription-translation (TX-TL) system provides valuable information about the performances of in silico designed riboregulators. We first propose a simple model that provides a quantitative definition of the dynamic range of a riboregulator. We further characterize two types of translational riboregulators composed of a cis-repressed (cr) and a trans-activating (ta) strand. At the DNA level we demonstrate that high concentrations of taDNA poisoned the activator until total shut off, in agreement with our model, and that relative dynamic ranges of riboregulators determined in vitro are in agreement with published in vivo data. At the RNA level, we show that this approach provides a fast and simple way to measure dissociation constants of functional riboregulators, in contrast to standard mobility-shift assays. Our method opens the route for using cell-free TX-TL systems for the quantitative characterization of functional riboregulators in order to improve their design in silico.

Mathematical Modeling of RNA-Based Architectures for Closed Loop Control of Gene Expression

Feedback allows biological systems to control gene expression precisely and reliably, even in the presence of uncertainty, by sensing and processing environmental changes. Taking inspiration from natural architectures, synthetic biologists have engineered feedback loops to tune the dynamics and improve the robustness and predictability of gene expression. However, experimental implementations of biomolecular control systems are still far from satisfying performance specifications typically achieved by electrical or mechanical control systems. To address this gap, we present mathematical models of biomolecular controllers that enable reference tracking, disturbance rejection, and tuning of the temporal response of gene expression. These controllers employ RNA transcriptional regulators to achieve closed loop control where feedback is introduced via molecular sequestration. Sensitivity analysis of the models allows us to identify which parameters influence the transient and steady state response of a target gene expression process, as well as which biologically plausible parameter values enable perfect reference tracking. We quantify performance using typical control theory metrics to characterize response properties and provide clear selection guidelines for practical applications. Our results indicate that RNA regulators are well-suited for building robust and precise feedback controllers for gene expression. Additionally, our approach illustrates several quantitative methods useful for assessing the performance of biomolecular feedback control systems.

Cell-Free Optogenetic Gene Expression System

Optogenetic tools provide a new and efficient way to dynamically program gene expression with unmatched spatiotemporal precision. To date, their vast potential remains untapped in the field of cell-free synthetic biology, largely due to the lack of simple and efficient light-switchable systems. Here, to bridge the gap between cell-free systems and optogenetics, we studied our previously engineered one component-based blue light-inducible Escherichia coli promoter in a cell-free environment through experimental characterization and mathematical modeling. We achieved >10-fold dynamic expression and demonstrated rapid and reversible activation of the target gene to generate oscillatory response. The deterministic model developed was able to recapitulate the system behavior and helped to provide quantitative insights to optimize dynamic response. This in vitro optogenetic approach could be a powerful new high-throughput screening technology for rapid prototyping of complex biological networks in both space and time without the need for chemical induction.

Achieving large dynamic range control of gene expression with a compact RNA transcription-translation regulator

RNA transcriptional regulators are emerging as versatile components for genetic network construction. However, these regulators suffer from incomplete repression in their OFF state, making their dynamic range less than that of their protein counterparts. This incomplete repression causes expression leak, which impedes the construction of larger synthetic regulatory networks as leak propagation can interfere with desired network function. To address this, we demonstrate how naturally derived antisense RNA-mediated transcriptional regulators can be configured to regulate both transcription and translation in a single compact RNA mechanism that functions in Escherichia coli. Using in vivo gene expression assays, we show that a combination of transcriptional termination and ribosome binding site sequestration increases repression from 85% to 98%, or activation from 10-fold to over 900-fold, in response to cognate antisense RNAs. We also show that orthogonal repressive versions of this mechanism can be created through engineering minimal antisense RNAs. Finally, to demonstrate the utility of this mechanism, we use it to reduce network leak in an RNA-only cascade. We anticipate these regulators will find broad use as synthetic biology moves beyond parts engineering to the design and construction of more sophisticated regulatory networks.

Respectful Modeling: Addressing Uncertainty in Dynamic System Models for Molecular Biology

Although there is still some skepticism in the biological community regarding the value and significance of quantitative computational modeling, important steps are continually being taken to enhance its accessibility and predictive power. We view these developments as essential components of an emerging 'respectful modeling' framework which has two key aims: (i) respecting the models themselves and facilitating the reproduction and update of modeling results by other scientists, and (ii) respecting the predictions of the models and rigorously quantifying the confidence associated with the modeling results. This respectful attitude will guide the design of higher-quality models and facilitate the use of models in modern applications such as engineering and manipulating microbial metabolism by synthetic biology.

Population Heterogeneity in the Epithelial to Mesenchymal Transition Is Controlled by NFAT and Phosphorylated Sp1

Epithelial to mesenchymal transition (EMT) is an essential differentiation program during tissue morphogenesis and remodeling. EMT is induced by soluble transforming growth factor β (TGF-β) family members, and restricted by vascular endothelial growth factor family members. While many downstream molecular regulators of EMT have been identified, these have been largely evaluated individually without considering potential crosstalk. In this study, we created an ensemble of dynamic mathematical models describing TGF-β induced EMT to better understand the operational hierarchy of this complex molecular program. We used ordinary differential equations (ODEs) to describe the transcriptional and post-translational regulatory events driving EMT. Model parameters were estimated from multiple data sets using multiobjective optimization, in combination with cross-validation. TGF-β exposure drove the model population toward a mesenchymal phenotype, while an epithelial phenotype was enhanced following vascular endothelial growth factor A (VEGF-A) exposure. Simulations predicted that the transcription factors phosphorylated SP1 and NFAT were master regulators promoting or inhibiting EMT, respectively. Surprisingly, simulations also predicted that a cellular population could exhibit phenotypic heterogeneity (characterized by a significant fraction of the population with both high epithelial and mesenchymal marker expression) if treated simultaneously with TGF-β and VEGF-A. We tested this prediction experimentally in both MCF10A and DLD1 cells and found that upwards of 45% of the cellular population acquired this hybrid state in the presence of both TGF-β and VEGF-A. We experimentally validated the predicted NFAT/Sp1 signaling axis for each phenotype response. Lastly, we found that cells in the hybrid state had significantly different functional behavior when compared to VEGF-A or TGF-β treatment alone. Together, these results establish a predictive mechanistic model of EMT susceptibility, and potentially reveal a novel signaling axis which regulates carcinoma progression through an EMT versus tubulogenesis response.

Using RNA as Molecular Code for Programming Cellular Function

RNA is involved in a wide-range of important molecular processes in the cell, serving diverse functions: regulatory, enzymatic, and structural. Together with its ease and predictability of design, these properties can lead RNA to become a useful handle for biological engineers with which to control the cellular machinery. By modifying the many RNA links in cellular processes, it is possible to reprogram cells toward specific design goals. We propose that RNA can be viewed as a molecular programming language that, together with protein-based execution platforms, can be used to rewrite wide ranging aspects of cellular function. In this review, we catalogue developments in the use of RNA parts, methods, and associated computational models that have contributed to the programmability of biology. We discuss how RNA part repertoires have been combined to build complex genetic circuits, and review recent applications of RNA-based parts and circuitry. We explore the future potential of RNA engineering and posit that RNA programmability is an important resource for firmly establishing an era of rationally designed synthetic biology.