Algorithms for ribosome traffic engineering and their potential in improving host cells' titer and growth rate

Specific codons control cellular resources and fitness

As cellular engineering progresses from simply overexpressing proteins to imparting complex phenotypes through multigene expression, judicious appropriation of cellular resources is essential. Since codon use is degenerate and biased, codons may control cellular resources at a translational level. We investigate how partitioning transfer RNA (tRNA) resources by incorporating dissimilar codon usage can drastically alter interdependence of expression level and burden on the host. By isolating the effect of individual codons' use during translation elongation while eliminating confounding factors, we show that codon choice can trans-regulate fitness of the host and expression of other heterologous or native genes. We correlate specific codon usage patterns with host fitness and derive a coding scheme for multigene expression called the Codon Health Index (CHI, χ). This empirically derived coding scheme (χ) enables the design of multigene expression systems that avoid catastrophic cellular burden and is robust across several proteins and conditions.

Scaling of nuclear numbers and their spatial arrangement in skeletal muscle cell size regulation

Many cells display considerable functional plasticity and depend on the regulation of numerous organelles and macromolecules for their maintenance. In large cells, organelles also need to be carefully distributed to supply the cell with essential resources and regulate intracellular activities. Having multiple copies of the largest eukaryotic organelle, the nucleus, epitomizes the importance of scaling gene products to large cytoplasmic volumes in skeletal muscle fibers. Scaling of intracellular constituents within mammalian muscle fibers is, however, poorly understood, but according to the myonuclear domain hypothesis, a single nucleus supports a finite amount of cytoplasm and is thus postulated to act autonomously, causing the nuclear number to be commensurate with fiber volume. In addition, the orderly peripheral distribution of myonuclei is a hallmark of normal cell physiology, as nuclear mispositioning is associated with impaired muscle function. Because underlying structures of complex cell behaviors are commonly formalized by scaling laws and thus emphasize emerging principles of size regulation, the work presented herein offers more of a unified conceptual platform based on principles from physics, chemistry, geometry, and biology to explore cell size-dependent correlations of the largest mammalian cell by means of scaling.

Modelling genetic stability in engineered cell populations

Predicting the evolution of engineered cell populations is a highly sought-after goal in biotechnology. While models of evolutionary dynamics are far from new, their application to synthetic systems is scarce where the vast combination of genetic parts and regulatory elements creates a unique challenge. To address this gap, we here-in present a framework that allows one to connect the DNA design of varied genetic devices with mutation spread in a growing cell population. Users can specify the functional parts of their system and the degree of mutation heterogeneity to explore, after which our model generates host-aware transition dynamics between different mutation phenotypes over time. We show how our framework can be used to generate insightful hypotheses across broad applications, from how a device's components can be tweaked to optimise long-term protein yield and genetic shelf life, to generating new design paradigms for gene regulatory networks that improve their functionality.

Modeling mRNA Translation With Ribosome Abortions

We derive a deterministic mathematical model for the flow of ribosomes along a mRNA called the ribosome flow model with extended objects and abortions (RFMEOA). This model incorporates important cellular features such as every ribosome covers several codons and they may detach from various regions along the track due to more realistic biological situations including phenomena of ribosome-ribosome collisions. We prove that the ribosome density profile along the mRNA in the RFMEOA and in particular, the protein production rate converge to a unique steady-state. Simulations of the RFMEOA demonstrate a surprising result that an increase in the initiation rate may sometimes lead to a decrease in the production rate. We believe that this model could be helpful to provide insight into the effects of premature termination on the protein expression and be useful for understanding and re-engineering the translation process.

On the gain of entrainment in the n-dimensional ribosome flow model

The ribosome flow model (RFM) is a phenomenological model for the flow of particles along a one-dimensional chain of n sites. It has been extensively used to study ribosome flow along the mRNA molecule during translation. When the transition rates along the chain are time-varying and jointly T-periodic the RFM entrains, i.e. every trajectory of the RFM converges to a unique T-periodic solution that depends on the transition rates, but not on the initial condition. In general, entrainment to periodic excitations like the 24 h solar day or the 50 Hz frequency of the electric grid is important in numerous natural and artificial systems. An interesting question, called the gain of entrainment (GOE) in the RFM, is whether proper coordination of the periodic translation rates along the mRNA can lead to a larger average protein production rate. Analysing the GOE in the RFM is non-trivial and partial results exist only for the RFM with dimensions n = 1, 2. We use a new approach to derive several results on the GOE in the general n-dimensional RFM. Perhaps surprisingly, we rigorously characterize several cases where there is no GOE, so to maximize the average production rate in these cases, the best choice is to use constant transition rates all along the chain.

Translation in the cell under fierce competition for shared resources: a mathematical model

During translation, mRNAs 'compete' for shared resources. Under stress conditions, during viral infection and also in high-throughput heterologous gene expression, these resources may become scarce, e.g. the pool of free ribosomes is starved, and then the competition may have a dramatic effect on the global dynamics of translation in the cell. We model this scenario using a network that includes m ribosome flow models (RFMs) interconnected via a pool of free ribosomes. Each RFM models ribosome flow along an mRNA molecule, and the pool models the shared resource. We assume that the number of mRNAs is large, so many ribosomes are attached to the mRNAs, and the pool is starved. Our analysis shows that adding an mRNA has an intricate effect on the total protein production. The new mRNA produces new proteins, but the other mRNAs produce less proteins, as the pool that feeds these mRNAs now has a smaller abundance of ribosomes. As the number of mRNAs increases, the marginal utility of adding another mRNA diminishes, and the total protein production rate saturates to a limiting value. We demonstrate our approach using an example of insulin protein production in a cell-free system.

Large-scale mRNA translation and the intricate effects of competition for the finite pool of ribosomes

We present a new theoretical framework for large-scale mRNA translation using a network of models called the ribosome flow model with Langmuir kinetics (RFMLK), interconnected via a pool of free ribosomes. The input to each RFMLK depends on the pool density, and it affects the initiation rate and potentially also the internal ribosome entry rates along each RFMLK. Ribosomes that detach from an RFMLK owing to termination or premature drop-off are fed back into the pool. We prove that the network always converges to a steady state, and study its sensitivity to variations in the parameters. For example, we show that if the drop-off rate at some site in some RFMLK is increased then the pool density increases and consequently the steady-state production rate in all the other RFMLKs increases. Surprisingly, we also show that modifying a parameter of a certain RFMLK can lead to arbitrary effects on the densities along the modified RFMLK, depending on the parameters in the entire network. We conclude that the competition for shared resources generates an indirect and intricate web of mutual effects between the mRNA molecules that must be accounted for in any analysis of translation.

Modeling the ribosomal small subunit dynamic in Saccharomyces cerevisiae based on TCP-seq data

Translation Complex Profile Sequencing (TCP-seq), a protocol that was developed and implemented on Saccharomyces cerevisiae, provides the footprints of the small subunit (SSU) of the ribosome (with additional factors) across the entire transcriptome of the analyzed organism. In this study, based on the TCP-seq data, we developed for the first-time a predictive model of the SSU density and analyzed the effect of transcript features on the dynamics of the SSU scan in the 5′UTR. Among others, our model is based on novel tools for detecting complex statistical relations tailored to TCP-seq. We quantitatively estimated the effect of several important features, including the context of the upstream AUG, the upstream ORF length and the mRNA folding strength. Specifically, we suggest that around 50% of the variance related to the read counts (RC) distribution near a start codon can be attributed to the AUG context score. We provide the first large scale direct quantitative evidence that shows that indeed AUG context affects the small sub-unit movement. In addition, we suggest that strong folding may cause the detachment of the SSU from the mRNA. We also identified a number of novel sequence motifs that can affect the SSU scan; some of these motifs affect transcription factors and RNA binding proteins. The results presented in this study provide a better understanding of the biophysical aspects related to the SSU scan along the 5′UTR and of translation initiation in S. cerevisiae, a fundamental step toward a comprehensive modeling of initiation.

Accelerating Whole-Cell Simulations of mRNA Translation Using a Dedicated Hardware

In recent years, intracellular biophysical simulations have been used with increasing frequency not only for answering basic scientific questions but also in the field of synthetic biology. However, since these models include networks of interaction between millions of components, they are extremely time-consuming and cannot run easily on parallel computers. In this study, we demonstrate for the first time a novel approach addressing this challenge by using a dedicated hardware designed specifically to simulate such processes. As a proof of concept, we specifically focus on mRNA translation, which is the process consuming most of the energy in the cell. We design a hardware that simulates translation in Escherichia coli and Saccharomyces cerevisiae for thousands of mRNAs and ribosomes, which is in orders of magnitude faster than a similar software solution. With the sharp increase in the amount of genomic data available today and the complexity of the corresponding models inferred from them, we believe that the strategy suggested here will become common and can be used among others for simulating entire cells with all gene expression steps.

Determinants of efficient modulation of ribosomal traffic jams

mRNA translation is the process which consumes most of the cellular energy. Thus, this process is under strong evolutionary selection for its optimization and rational optimization or reduction of the translation efficiency can impact the cell growth rate. Algorithms for modulating cell growth rate can have various applications in biotechnology, medicine, and agriculture. In this study, we demonstrate that the analysis of these algorithms can also be used for understanding translation.

We specifically describe and analyze various generic algorithms, based on comprehensive computational models and whole cell simulations of translation, for introducing silent mutations that can either reduce or increase ribosomal traffic jams along the mRNA. As a result, more or less resources are available, for the cell, promoting improved or reduced cells growth-rate, respectively. We then explore the cost of these algorithms' performance, in terms of their computational time, the number of mutations they introduce, the modified genomic region, the effect on local translation rates, and the properties of the modified genes.

Among others, we show that mRNA levels of a gene are much stronger predictors for the effect of its engineering on the ribosomal pool than the ribosomal density of the gene. We also demonstrate that the mutations at the ends of the coding regions have a stronger effect on the ribosomal pool. Furthermore, we report two optimization algorithms that exhibit a tread-off between the number of mutations they introduce and their executing time.

The reported results here are fundamental both for understanding the biophysics and evolution of translation, as well as for developing efficient approaches for its engineering.

A Modelling Framework Linking Resource-Based Stochastic Translation to the Optimal Design of Synthetic Constructs

Simple Summary

In synthetic biology, it is commonplace to design and insert gene expression constructs into cells for the production of useful proteins. In order to maximise production yield, it is useful to predict the performance of these “engineered cells” in advance of conducting experiments. This is typically a complex task, which in recent years has motivated the use of “whole-cell models” (WCMs) that act as computational tools for predicting different aspects of cell growth. Many useful WCMs exist, however a common problem is their over-simplification of ribosome movement on mRNA transcripts during translation. WCMs typically don’t consider that, for constructs with inefficient (“slow”) codons, ribosomes can stall and form “traffic jams”, thereby becoming unavailable for translation of other proteins. To more accurately address these scenarios, we have built a computational framework that combines whole-cell modelling with a detailed account of ribosome movement on mRNA. We show how our framework can be used to link the modular design of a gene expression construct (via its promoter, ribosome binding site and codon composition) to protein yield during continuous cell culture, with a particular focus on how the optimal design can change over time in the presence or absence of “slow” codons.

Abstract

The effect of gene expression burden on engineered cells has motivated the use of “whole-cell models” (WCMs) that use shared cellular resources to predict how unnatural gene expression affects cell growth. A common problem with many WCMs is their inability to capture translation in sufficient detail to consider the impact of ribosomal queue formation on mRNA transcripts. To address this, we have built a “stochastic cell calculator” (StoCellAtor) that combines a modified TASEP with a stochastic implementation of an existing WCM. We show how our framework can be used to link a synthetic construct’s modular design (promoter, ribosome binding site (RBS) and codon composition) to protein yield during continuous culture, with a particular focus on the effects of low-efficiency codons and their impact on ribosomal queues. Through our analysis, we recover design principles previously established in our work on burden-sensing strategies, namely that changing promoter strength is often a more efficient way to increase protein yield than RBS strength. Importantly, however, we show how these design implications can change depending on both the duration of protein expression, and on the presence of ribosomal queues.