The Timing of Transcriptional Regulation in Synthetic Gene Circuits

Inferring delays in partially observed gene regulation processes

MOTIVATION

Cell function is regulated by gene regulatory networks (GRNs) defined by protein-mediated interaction between constituent genes. Despite advances in experimental techniques, we can still measure only a fraction of the processes that govern GRN dynamics. To infer the properties of GRNs using partial observation, unobserved sequential processes can be replaced with distributed time delays, yielding non-Markovian models. Inference methods based on the resulting model suffer from the curse of dimensionality.

RESULTS

We develop a simulation-based Bayesian MCMC method employing an approximate likelihood for the efficient and accurate inference of GRN parameters when only some of their products are observed. We illustrate our approach using a two-step activation model: an activation signal leads to the accumulation of an unobserved regulatory protein, which triggers the expression of observed fluorescent proteins. With prior information about observed fluorescent protein synthesis, our method successfully infers the dynamics of the unobserved regulatory protein. We can estimate the delay and kinetic parameters characterizing target regulation including transcription, translation, and target searching of an unobserved protein from experimental measurements of the products of its target gene. Our method is scalable and can be used to analyze non-Markovian models with hidden components.

AVAILABILITY AND IMPLEMENTATION

Our code is implemented in R and is freely available with a simple example data at https://github.com/Mathbiomed/SimMCMC.

Systematic inference identifies a major source of heterogeneity in cell signaling dynamics: The rate-limiting step number

Identifying the sources of cell-to-cell variability in signaling dynamics is essential to understand drug response variability and develop effective therapeutics. However, it is challenging because not all signaling intermediate reactions can be experimentally measured simultaneously. This can be overcome by replacing them with a single random time delay, but the resulting process is non-Markovian, making it difficult to infer cell-to-cell heterogeneity in reaction rates and time delays. To address this, we developed an efficient and scalable moment-based Bayesian inference method (MBI) with a user-friendly computational package that infers cell-to-cell heterogeneity in the non-Markovian signaling process. We applied MBI to single-cell expression profiles from promoters responding to antibiotics and discovered a major source of cell-to-cell variability in antibiotic stress response: the number of rate-limiting steps in signaling cascades. This knowledge can help identify effective therapies that destroy all pathogenic or cancer cells, and the approach can be applied to precision medicine.

Mechanobiology of immune cells: Messengers, receivers and followers in leishmaniasis aiding synthetic devices

Cytokines are influential molecules which can direct cells behavior. In this review, cytokines are referred as messengers, immune cells which respond to cytokine stimulus are referred as receivers and the immune cells which gets modulated due to their plasticity induced by infectious pathogen leishmania, are referred as followers. The advantage of plasticity of cells is taken by the parasite to switch them from parasite eliminating form to parasite survival favoring form through a process called as reciprocity which is undergone by cytokines, wherein pro-inflammatory to anti-inflammatory switch occur rendering immune cell population to switch their phenotype. Detailed study of this switch can help in identification of important targets which can help in restoring the phenotype to parasite eliminating form and this can be done through synthetic circuit, finding its wider applicability in therapeutics.

•

Cytokines as messengers for governing reciprocity in infection.

•

Leishmania induces reciprocity modulating the immune cells plasticity.

•

Reciprocity of cytokines identifies important target for therapeutics.

•

Therapeutic targets aiding the design of synthetic devices to combat infection.

Hierarchical Bayesian models of transcriptional and translational regulation processes with delays

MOTIVATION

Simultaneous recordings of gene network dynamics across large populations have revealed that cell characteristics vary considerably even in clonal lines. Inferring the variability of parameters that determine gene dynamics is key to understanding cellular behavior. However, this is complicated by the fact that the outcomes and effects of many reactions are not observable directly. Unobserved reactions can be replaced with time delays to reduce model dimensionality and simplify inference. However, the resulting models are non-Markovian, and require the development of new inference techniques.

RESULTS

We propose a non-Markovian, hierarchical Bayesian inference framework for quantifying the variability of cellular processes within and across cells in a population. We illustrate our approach using a delayed birth-death process. In general, a distributed delay model, rather than a popular fixed delay model, is needed for inference, even if only mean reaction delays are of interest. Using in silico and experimental data we show that the proposed hierarchical framework is robust and leads to improved estimates compared to its non-hierarchical counterpart. We apply our method to data obtained using time-lapse microscopy and infer the parameters that describe the dynamics of protein production at the single cell and population level. The mean delays in protein production are larger than previously reported, have a coefficient of variation of around 0.2 across the population, and are not strongly correlated with protein production or growth rates.

AVAILABILITY AND IMPLEMENTATION

Accompanying code in Python is available at https://github.com/mvcortez/Bayesian-Inference.

CONTACT

kresimir.josic@gmail.com or jaekkim@kaist.ac.kr or cbskust@korea.ac.kr.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Design and Evaluation of Synthetic RNA-Based Incoherent Feed-Forward Loop Circuits

RNA-based regulators are promising tools for building synthetic biological systems that provide a powerful platform for achieving a complex regulation of transcription and translation. Recently, de novo-designed synthetic RNA regulators, such as the small transcriptional activating RNA (STAR), toehold switch (THS), and three-way junction (3WJ) repressor, have been utilized to construct RNA-based synthetic gene circuits in living cells. In this work, we utilized these regulators to construct type 1 incoherent feed-forward loop (IFFL) circuits in vivo and explored their dynamic behaviors. A combination of a STAR and 3WJ repressor was used to construct an RNA-only IFFL circuit. However, due to the fast kinetics of RNA–RNA interactions, there was no significant timescale difference between the direct activation and the indirect inhibition, that no pulse was observed in the experiments. These findings were confirmed with mechanistic modeling and simulation results for a wider range of conditions. To increase delay in the inhibition pathway, we introduced a protein synthesis process to the circuit and designed an RNA–protein hybrid IFFL circuit using THS and TetR protein. Simulation results indicated that pulse generation could be achieved with this RNA–protein hybrid model, and this was further verified with experimental realization in E. coli. Our findings demonstrate that while RNA-based regulators excel in speed as compared to protein-based regulators, the fast reaction kinetics of RNA-based regulators could also undermine the functionality of a circuit (e.g., lack of significant timescale difference). The agreement between experiments and simulations suggests that the mechanistic modeling can help debug issues and validate the hypothesis in designing a new circuit. Moreover, the applicability of the kinetic parameters extracted from the RNA-only circuit to the RNA–protein hybrid circuit also indicates the modularity of RNA-based regulators when used in a different context. We anticipate the findings of this work to guide the future design of gene circuits that rely heavily on the dynamics of RNA-based regulators, in terms of both modeling and experimental realization.

Systematic Quantification of Sequence and Structural Determinants Controlling mRNA stability in Bacterial Operons

mRNA degradation is a central process that affects all gene expression levels, and yet, the determinants that control mRNA decay rates remain poorly characterized. Here, we applied a synthetic biology, learn-by-design approach to elucidate the sequence and structural determinants that control mRNA stability in bacterial operons. We designed, constructed, and characterized 82 operons in Escherichia coli, systematically varying RNase binding site characteristics, translation initiation rates, and transcriptional terminator efficiencies in the 5' untranslated region (UTR), intergenic, and 3' UTR regions, followed by measuring their mRNA levels using reverse transcription quantitative polymerase chain reaction (RT-qPCR) assays during exponential growth. We show that introducing long single-stranded RNA into 5' UTRs reduced mRNA levels by up to 9.4-fold and that lowering translation rates reduced mRNA levels by up to 11.8-fold. We also found that RNase binding sites in intergenic regions had much lower effects on mRNA levels. Surprisingly, changing the transcriptional termination efficiency or introducing long single-stranded RNA into 3' UTRs had no effect on upstream mRNA levels. From these measurements, we developed and validated biophysical models of ribosome protection and RNase activity with excellent quantitative agreement. We also formulated design rules to rationally control a mRNA's stability, facilitating the automated design of engineered genetic systems with desired functionalities.

Stochastic fluctuations in apoptotic threshold of tumour cells can enhance apoptosis and combat fractional killing

Fractional killing, which is a significant impediment to successful chemotherapy, is observed even in a population of genetically identical cancer cells exposed to apoptosis-inducing agents. This phenomenon arises not from genetic mutation but from cell-to-cell variation in the activation timing and level of the proteins that regulates apoptosis. To understand the mechanism behind the phenomenon, we formulate complex fractional killing processes as a first-passage time (FPT) problem with a stochastically fluctuating boundary. Analytical calculations are performed for the FPT distribution in a toy model of stochastic p53 gene expression, where the cancer cell is killed only when the p53 expression level crosses an active apoptotic threshold. Counterintuitively, we find that threshold fluctuations can effectively enhance cellular killing by significantly decreasing the mean time that the p53 protein reaches the threshold level for the first time. Moreover, faster fluctuations lead to the killing of more cells. These qualitative results imply that fluctuations in threshold are a non-negligible stochastic source, and can be taken as a strategy for combating fractional killing of cancer cells.

Control strategies for the timing of intracellular events

While the timing of intracellular events is essential for many cellular processes, gene expression inside a single cell can exhibit substantial cell-to-cell variability, raising the question of how cells ensure precision in event timing despite such stochasticity. We address this question by analyzing a biologically reasonable model of gene expression in the context of first passage time (FPT), focusing on two experimentally measurable statistics: mean FPT (MFPT) and timing variability (TV). We show that (1) transcriptional burst size (BS) and burst frequency (BF) can minimize the TV; (2) translational BS monotonically reduces the MFPT to a nonzero low bound; (3) the timescale of promoter kinetics can minimize both the MFPT and the TV, depending on the ratio of the on-switching rate over the off-switching rate; and (4) positive feedback regulation of any form can all minimize the TV, whereas negative feedback regulation of transcriptional BF or BS always enhances the TV. These control strategies can have broad implications for diverse cellular processes relying on precise temporal triggering of events.

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.

Metabolic perceptrons for neural computing in biological systems

Synthetic biological circuits are promising tools for developing sophisticated systems for medical, industrial, and environmental applications. So far, circuit implementations commonly rely on gene expression regulation for information processing using digital logic. Here, we present a different approach for biological computation through metabolic circuits designed by computer-aided tools, implemented in both whole-cell and cell-free systems. We first combine metabolic transducers to build an analog adder, a device that sums up the concentrations of multiple input metabolites. Next, we build a weighted adder where the contributions of the different metabolites to the sum can be adjusted. Using a computational model fitted on experimental data, we finally implement two four-input perceptrons for desired binary classification of metabolite combinations by applying model-predicted weights to the metabolic perceptron. The perceptron-mediated neural computing introduced here lays the groundwork for more advanced metabolic circuits for rapid and scalable multiplex sensing.

Diversity, versatility and complexity of bacterial gene regulation mechanisms: opportunities and drawbacks for applications in synthetic biology

Gene expression occurs in two essential steps: transcription and translation. In bacteria, the two processes are tightly coupled in time and space, and highly regulated. Tight regulation of gene expression is crucial. It limits wasteful consumption of resources and energy, prevents accumulation of potentially growth inhibiting reaction intermediates, and sustains the fitness and potential virulence of the organism in a fluctuating, competitive and frequently stressful environment. Since the onset of studies on regulation of enzyme synthesis, numerous distinct regulatory mechanisms modulating transcription and/or translation have been discovered. Mostly, various regulatory mechanisms operating at different levels in the flow of genetic information are used in combination to control and modulate the expression of a single gene or operon. Here, we provide an extensive overview of the very diverse and versatile bacterial gene regulatory mechanisms with major emphasis on their combined occurrence, intricate intertwinement and versatility. Furthermore, we discuss the potential of well-characterized basal expression and regulatory elements in synthetic biology applications, where they may ensure orthogonal, predictable and tunable expression of (heterologous) target genes and pathways, aiming at a minimal burden for the host.

Bayesian inference of distributed time delay in transcriptional and translational regulation

MOTIVATION

Advances in experimental and imaging techniques have allowed for unprecedented insights into the dynamical processes within individual cells. However, many facets of intracellular dynamics remain hidden, or can be measured only indirectly. This makes it challenging to reconstruct the regulatory networks that govern the biochemical processes underlying various cell functions. Current estimation techniques for inferring reaction rates frequently rely on marginalization over unobserved processes and states. Even in simple systems this approach can be computationally challenging, and can lead to large uncertainties and lack of robustness in parameter estimates. Therefore we will require alternative approaches to efficiently uncover the interactions in complex biochemical networks.

RESULTS

We propose a Bayesian inference framework based on replacing uninteresting or unobserved reactions with time delays. Although the resulting models are non-Markovian, recent results on stochastic systems with random delays allow us to rigorously obtain expressions for the likelihoods of model parameters. In turn, this allows us to extend MCMC methods to efficiently estimate reaction rates, and delay distribution parameters, from single-cell assays. We illustrate the advantages, and potential pitfalls, of the approach using a birth-death model with both synthetic and experimental data, and show that we can robustly infer model parameters using a relatively small number of measurements. We demonstrate how to do so even when only the relative molecule count within the cell is measured, as in the case of fluorescence microscopy.

AVAILABILITY AND IMPLEMENTATION

Accompanying code in R is available at https://github.com/cbskust/DDE_BD.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Intracellular Noise Level Determines Ratio Control Strategy Confined by Speed-Accuracy Trade-off

Robust and precise ratio control of heterogeneous phenotypes within an isogenic population is an essential task, especially in the development and differentiation of a large number of cells such as bacteria, sensory receptors, and blood cells. However, the mechanisms of such ratio control are poorly understood. Here, we employ experimental and mathematical techniques to understand the combined effects of signal induction and gene expression stochasticity on phenotypic multimodality. We identify two strategies to control phenotypic ratios from an initially homogeneous population, suitable roughly to high-noise and low-noise intracellular environments, and we show that both can be used to generate precise fractional differentiation. In noisy gene expression contexts, such as those found in bacteria, induction within the circuit's bistable region is enough to cause noise-induced bimodality within a feasible time frame. However, in less noisy contexts, such as tightly controlled eukaryotic systems, spontaneous state transitions are rare and hence bimodality needs to be induced with a controlled pulse of induction that falls outside the bistable region. Finally, we show that noise levels, system response time, and ratio tuning accuracy impose trade-offs and limitations on both ratio control strategies, which guide the selection of strategy alternatives.