The RNA polymerase flow model of gene transcription

Model-based characterization of the equilibrium dynamics of transcription initiation and promoter-proximal pausing in human cells

In metazoans, both transcription initiation and the escape of RNA polymerase (RNAP) from promoter-proximal pausing are key rate-limiting steps in gene expression. These processes play out at physically proximal sites on the DNA template and appear to influence one another through steric interactions. Here, we examine the dynamics of these processes using a combination of statistical modeling, simulation, and analysis of real nascent RNA sequencing data. We develop a simple probabilistic model that jointly describes the kinetics of transcription initiation, pause-escape, and elongation, and the generation of nascent RNA sequencing read counts under steady-state conditions. We then extend this initial model to allow for variability across cells in promoter-proximal pause site locations and steric hindrance of transcription initiation from paused RNAPs. In an extensive series of simulations, we show that this model enables accurate estimation of initiation and pause-escape rates. Furthermore, we show by simulation and analysis of real data that pause-escape is often strongly rate-limiting and that steric hindrance can dramatically reduce initiation rates. Our modeling framework is applicable to a variety of inference problems, and our software for estimation and simulation is freely available.

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.

A theoretical framework to analyse the flow of particles in a dynamical system with stochastic transition rates and site capacities

We study the stochasticity in a dynamical model: ribosome flow model with different site sizes that models the unidirectional movement of particles controlled by transition rates along a lattice having different site sizes. Our work models the parameters as random variables with known distributions and investigates the steady-state flow rate under this notion by using tools from the random matrix theory. Some closed-form theoretical results are derived for the steady-state flow rate under some restrictive assumptions such as random variables being independent and identically distributed. Furthermore, for arbitrary but bounded stochastic transition rates, stochastic site capacities, or both, we establish bounds for the steady-state flow rate. Our analysis can be generalized and applied to study the flow of particles in numerous transport systems in the stochastic environment.

Variability in mRNA translation: a random matrix theory approach

The rate of mRNA translation depends on the initiation, elongation, and termination rates of ribosomes along the mRNA. These rates depend on many "local" factors like the abundance of free ribosomes and tRNA molecules in the vicinity of the mRNA molecule. All these factors are stochastic and their experimental measurements are also noisy. An important question is how protein production in the cell is affected by this considerable variability. We develop a new theoretical framework for addressing this question by modeling the rates as identically and independently distributed random variables and using tools from random matrix theory to analyze the steady-state production rate. The analysis reveals a principle of universality: the average protein production rate depends only on the of the set of possible values that the random variable may attain. This explains how total protein production can be stabilized despite the overwhelming stochasticticity underlying cellular processes.

Controllability Analysis and Control Synthesis for the Ribosome Flow Model

The ribosomal density along different parts of the coding regions of the mRNA molecule affects various fundamental intracellular phenomena including: protein production rates, global ribosome allocation and organismal fitness, ribosomal drop off, co-translational protein folding, mRNA degradation, and more. Thus, regulating translation in order to obtain a desired ribosomal profile along the mRNA molecule is an important biological problem. We study this problem by using a dynamical model for mRNA translation, called the ribosome flow model (RFM). In the RFM, the mRNA molecule is modeled as an ordered chain of $n$ sites. The RFM includes $n$ state-variables describing the ribosomal density profile along the mRNA molecule, and the transition rates from each site to the next are controlled by $n+1$ positive constants. To study the problem of controlling the density profile, we consider some or all of the transition rates as time-varying controls. We consider the following problem: given an initial and a desired ribosomal density profile in the RFM, determine the time-varying values of the transition rates that steer the system to the desired density profile, if they exist. More specifically, we consider two control problems. In the first, all transition rates can be regulated separately, and the goal is to steer the ribosomal density profile and the protein production rate from a given initial value to a desired value. In the second problem, one or more transition rates are jointly regulated by a single scalar control, and the goal is to steer the production rate to a desired value within a certain set of feasible values. In the first case, we show that the system is controllable, i.e., the control is powerful enough to steer the system to any desired value in finite time, and provide simple closed-form expressions for constant positive control functions (or transition rates) that asymptotically steer the system to the desired value. In the second case, we show that the system is controllable, and provide a simple algorithm for determining the constant positive control value that asymptotically steers the system to the desired value. We discuss some of the biological implications of these results.

Enhancing heterologous expression in Chlamydomonas reinhardtii by transcript sequence optimization

Various species of microalgae have recently emerged as promising host-organisms for use in biotechnology industries due to their unique properties. These include efficient conversion of sunlight into organic compounds, the ability to grow in extreme conditions and the occurrence of numerous post-translational modification pathways. However, the inability to obtain high levels of nuclear heterologous gene expression in microalgae hinders the development of the entire field. To overcome this limitation, we analyzed different sequence optimization algorithms while studying the effect of transcript sequence features on heterologous expression in the model microalga Chlamydomonas reinhardtii, whose genome consists of rare features such as a high GC content. Based on the analysis of genomic data, we created eight unique sequences coding for a synthetic ferredoxin-hydrogenase enzyme, used here as a reporter gene. Following in silico design, these synthetic genes were transformed into the C. reinhardtii nucleus, after which gene expression levels were measured. The empirical data, measured in vivo show a discrepancy of up to 65-fold between the different constructs. In this work we demonstrate how the combination of computational methods and our empirical results enable us to learn about the way gene expression is encoded in the C. reinhardtii transcripts. We describe the deleterious effect on overall expression of codons encoding for splicing signals. Subsequently, our analysis shows that utilization of a frequent subset of preferred codons results in elevated transcript levels, and that mRNA folding energy in the vicinity of translation initiation significantly affects gene expression.

A code for transcription elongation speed

The two major steps of gene expression are transcription and translation. While hundreds of studies regarding the effect of sequence features on the translation elongation process have been published, very few connect sequence features to the transcription elongation rate. We suggest, for the first time, that short transcript sub-sequences have a typical effect on RNA polymerase (RNAP) speed: we show that nucleotide 5-mers tend to have typical RNAP speed (or transcription rate), which is consistent along different parts of genes and among different groups of genes with high correlation. We also demonstrate that relative RNAP speed correlates with mRNA levels of endogenous and heterologous genes. Furthermore, we show that the estimated transcription and translation elongation rates correlate in endogenous genes. Finally, we demonstrate that our results are consistent for different high resolution experimental measurements of RNAP densities. These results suggest for the first time that transcription elongation is partly encoded in the transcript, affected by the codon-usage, and optimized by evolution with a significant effect on gene expression and organismal fitness.

Ribosome flow model with extended objects

We study a deterministic mechanistic model for the flow of ribosomes along the mRNA molecule, called the ribosome flow model with extended objects (RFMEO). This model encapsulates many realistic features of translation including non-homogeneous transition rates along mRNA, the fact that every ribosome covers several codons, and the fact that ribosomes cannot overtake one another. The RFMEO is a mean-field approximation of an important model from statistical mechanics called the totally asymmetric simple exclusion process with extended objects (TASEPEO). We demonstrate that the RFMEO describes biophysical aspects of translation better than previous mean-field approximations, and that its predictions correlate well with those of TASEPEO. However, unlike TASEPEO, the RFMEO is amenable to rigorous analysis using tools from systems and control theory. We show that the ribosome density profile along the mRNA in the RFMEO converges to a unique steady-state density that depends on the length of the mRNA, the transition rates along it, and the number of codons covered by every ribosome, but not on the initial density of ribosomes along the mRNA. In particular, the protein production rate also converges to a unique steady state. Furthermore, if the transition rates along the mRNA are periodic with a common period T then the ribosome density along the mRNA and the protein production rate converge to a unique periodic pattern with period T, that is, the model entrains to periodic excitations in the transition rates. Analysis and simulations of the RFMEO demonstrate several counterintuitive results. For example, increasing the ribosome footprint may sometimes lead to an increase in the production rate. Also, for large values of the footprint the steady-state density along the mRNA may be quite complex (e.g. with quasi-periodic patterns) even for relatively simple (and non-periodic) transition rates along the mRNA. This implies that inferring the transition rates from the ribosome density may be non-trivial. We believe that the RFMEO could be useful for modelling, understanding and re-engineering translation as well as other important biological processes.

A deterministic mathematical model for bidirectional excluded flow with Langmuir kinetics

In many important cellular processes, including mRNA translation, gene transcription, phosphotransfer, and intracellular transport, biological "particles" move along some kind of "tracks". The motion of these particles can be modeled as a one-dimensional movement along an ordered sequence of sites. The biological particles (e.g., ribosomes or RNAPs) have volume and cannot surpass one another. In some cases, there is a preferred direction of movement along the track, but in general the movement may be bidirectional, and furthermore the particles may attach or detach from various regions along the tracks. We derive a new deterministic mathematical model for such transport phenomena that may be interpreted as a dynamic mean-field approximation of an important model from mechanical statistics called the asymmetric simple exclusion process (ASEP) with Langmuir kinetics. Using tools from the theory of monotone dynamical systems and contraction theory we show that the model admits a unique steady-state, and that every solution converges to this steady-state. Furthermore, we show that the model entrains (or phase locks) to periodic excitations in any of its forward, backward, attachment, or detachment rates. We demonstrate an application of this phenomenological transport model for analyzing ribosome drop off in mRNA translation.

A resource dependent protein synthesis model for evaluating synthetic circuits

Reliable in silico design of synthetic gene networks necessitates novel approaches to model the process of protein synthesis under the influence of limited resources. We present such a novel protein synthesis model which originates from the Ribosome Flow Model and among other things describes the movement of RNA-polymerase and ribosomes on mRNA and DNA templates, respectively. By analyzing the convergence properties of this model based upon geometric considerations, we present additional insights into the dynamic mechanisms of the process of protein synthesis. Further, we demonstrate how this model can be used to evaluate the performance of synthetic gene circuits under different loading scenarios.

On the Ribosomal Density that Maximizes Protein Translation Rate

During mRNA translation, several ribosomes attach to the same mRNA molecule simultaneously translating it into a protein. This pipelining increases the protein translation rate. A natural and important question is what ribosomal density maximizes the protein translation rate. Using mathematical models of ribosome flow along both a linear and a circular mRNA molecules we prove that typically the steady-state protein translation rate is maximized when the ribosomal density is one half of the maximal possible density. We discuss the implications of our results to endogenous genes under natural cellular conditions and also to synthetic biology.

Predictive biophysical modeling and understanding of the dynamics of mRNA translation and its evolution

mRNA translation is the fundamental process of decoding the information encoded in mRNA molecules by the ribosome for the synthesis of proteins. The centrality of this process in various biomedical disciplines such as cell biology, evolution and biotechnology, encouraged the development of dozens of mathematical and computational models of translation in recent years. These models aimed at capturing various biophysical aspects of the process. The objective of this review is to survey these models, focusing on those based and/or validated on real large-scale genomic data. We consider aspects such as the complexity of the models, the biophysical aspects they regard and the predictions they may provide. Furthermore, we survey the central systems biology discoveries reported on their basis. This review demonstrates the fundamental advantages of employing computational biophysical translation models in general, and discusses the relative advantages of the different approaches and the challenges in the field.

Sensitivity of mRNA Translation

Using the dynamic mean-field approximation of the totally asymmetric simple exclusion process (TASEP), we investigate the effect of small changes in the initiation, elongation, and termination rates along the mRNA strand on the steady-state protein translation rate. We show that the sensitivity of mRNA translation is equal to the sensitivity of the maximal eigenvalue of a symmetric, nonnegative, tridiagonal, and irreducible matrix. This leads to new analytical results as well as efficient numerical schemes that are applicable for large-scale models. Our results show that in the usual endogenous case, when initiation is more rate-limiting than elongation, the sensitivity of the translation rate to small mutations rapidly increases towards the 5′ end of the ORF. When the initiation rate is high, as may be the case for highly expressed and/or heterologous optimized genes, the maximal sensitivity is with respect to the elongation rates at the middle of the mRNA strand. We also show that the maximal possible effect of a small increase/decrease in any of the rates along the mRNA is an increase/decrease of the same magnitude in the translation rate. These results are in agreement with previous molecular evolutionary and synthetic biology experimental studies.

Maximizing protein translation rate in the non-homogeneous ribosome flow model: a convex optimization approach

Translation is an important stage in gene expression. During this stage, macro-molecules called ribosomes travel along the mRNA strand linking amino acids together in a specific order to create a functioning protein. An important question, related to many biomedical disciplines, is how to maximize protein production. Indeed, translation is known to be one of the most energy-consuming processes in the cell, and it is natural to assume that evolution shaped this process so that it maximizes the protein production rate. If this is indeed so then one can estimate various parameters of the translation machinery by solving an appropriate mathematical optimization problem. The same problem also arises in the context of synthetic biology, namely, re-engineer heterologous genes in order to maximize their translation rate in a host organism. We consider the problem of maximizing the protein production rate using a computational model for translation-elongation called the ribosome flow model (RFM). This model describes the flow of the ribosomes along an mRNA chain of length n using a set of n first-order nonlinear ordinary differential equations. It also includes n + 1 positive parameters: the ribosomal initiation rate into the mRNA chain, and n elongation rates along the chain sites. We show that the steady-state translation rate in the RFM is a strictly concave function of its parameters. This means that the problem of maximizing the translation rate under a suitable constraint always admits a unique solution, and that this solution can be determined using highly efficient algorithms for solving convex optimization problems even for large values of n. Furthermore, our analysis shows that the optimal translation rate can be computed based only on the optimal initiation rate and the elongation rate of the codons near the beginning of the ORF. We discuss some applications of the theoretical results to synthetic biology, molecular evolution, and functional genomics.

Maximizing Protein Translation Rate in the Ribosome Flow Model: The Homogeneous Case

Gene translation is the process in which intracellular macro-molecules, called ribosomes, decode genetic information in the mRNA chain into the corresponding proteins. Gene translation includes several steps. During the elongation step, ribosomes move along the mRNA in a sequential manner and link amino-acids together in the corresponding order to produce the proteins. The homogeneous ribosome flow model (HRFM) is a deterministic computational model for translation-elongation under the assumption of constant elongation rates along the mRNA chain. The HRFM is described by a set of n first-order nonlinear ordinary differential equations, where n represents the number of sites along the mRNA chain. The HRFM also includes two positive parameters: ribosomal initiation rate and the (constant) elongation rate. In this paper, we show that the steady-state translation rate in the HRFM is a concave function of its parameters. This means that the problem of determining the parameter values that maximize the translation rate is relatively simple. Our results may contribute to a better understanding of the mechanisms and evolution of translation-elongation. We demonstrate this by using the theoretical results to estimate the initiation rate in M. musculus embryonic stem cell. The underlying assumption is that evolution optimized the translation mechanism. For the infinite-dimensional HRFM, we derive a closed-form solution to the problem of determining the initiation and transition rates that maximize the protein translation rate. We show that these expressions provide good approximations for the optimal values in the n-dimensional HRFM already for relatively small values of n. These results may have applications for synthetic biology where an important problem is to re-engineer genomic systems in order to maximize the protein production rate.