Ribosome Flow Model on a Ring

Persistence and stability of generalized ribosome flow models with time-varying transition rates

In this paper some important qualitative dynamical properties of generalized ribosome flow models are studied. Ribosome flow models known from the literature are generalized by allowing an arbitrary directed network structure between compartments, and by assuming general time-varying rate functions corresponding to the transitions. Persistence of the dynamics is shown using the chemical reaction network (CRN) representation of the system where the state variables correspond to ribosome density and the amount of free space in the compartments. The L1 contractivity of solutions is also proved in the case of periodic reaction rates having the same period. Further we prove the stability of different compartmental structures including strongly connected ones with entropy-like logarithmic Lyapunov functions through embedding the model into a weakly reversible CRN with time-varying reaction rates in a reduced state space. Moreover, it is shown that different Lyapunov functions may be assigned to the same model depending on the non-unique factorization of the reaction rates. The results are illustrated through several examples with biological meaning including the classical ribosome flow model on a ring.

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.

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.

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.

Optimal Translation Along a Circular mRNA

The ribosome flow model on a ring (RFMR) is a deterministic model for ribosome flow along a circularized mRNA. We derive a new spectral representation for the optimal steady-state production rate and the corresponding optimal steady-state ribosomal density in the RFMR. This representation has several important advantages. First, it provides a simple and numerically stable algorithm for determining the optimal values even in very long rings. Second, it enables efficient computation of the sensitivity of the optimal production rate to small changes in the transition rates along the mRNA. Third, it implies that the optimal steady-state production rate is a strictly concave function of the transition rates. Maximizing the optimal steady-state production rate with respect to the rates under an affine constraint on the rates thus becomes a convex optimization problem that admits a unique solution. This solution can be determined numerically using highly efficient algorithms. This optimization problem is important, for example, when re-engineering heterologous genes in a host organism. We describe the implications of our results to this and other aspects of translation.

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 deterministic model for one-dimensional excluded flow with local interactions

Natural phenomena frequently involve a very large number of interacting molecules moving in confined regions of space. Cellular transport by motor proteins is an example of such collective behavior. We derive a deterministic compartmental model for the unidirectional flow of particles along a one-dimensional lattice of sites with nearest-neighbor interactions between the particles. The flow between consecutive sites is governed by a “soft” simple exclusion principle and by attracting or repelling forces between neighboring particles. Using tools from contraction theory, we prove that the model admits a unique steady-state and that every trajectory converges to this steady-state. Analysis and simulations of the effect of the attracting and repelling forces on this steady-state highlight the crucial role that these forces may play in increasing the steady-state flow, and reveal that this increase stems from the alleviation of traffic jams along the lattice. Our theoretical analysis clarifies microscopic aspects of complex multi-particle dynamic processes.

Optimal Down Regulation of mRNA Translation

Down regulation of mRNA translation is an important problem in various bio-medical domains ranging from developing effective medicines for tumors and for viral diseases to developing attenuated virus strains that can be used for vaccination. Here, we study the problem of down regulation of mRNA translation using a mathematical model called the ribosome flow model (RFM). In the RFM, the mRNA molecule is modeled as a chain of n sites. The flow of ribosomes between consecutive sites is regulated by n + 1 transition rates. Given a set of feasible transition rates, that models the outcome of all possible mutations, we consider the problem of maximally down regulating protein production by altering the rates within this set of feasible rates. Under certain conditions on the feasible set, we show that an optimal solution can be determined efficiently. We also rigorously analyze two special cases of the down regulation optimization problem. Our results suggest that one must focus on the position along the mRNA molecule where the transition rate has the strongest effect on the protein production rate. However, this rate is not necessarily the slowest transition rate along the mRNA molecule. We discuss some of the biological implications of these results.

A model for competition for ribosomes in the cell

A single mammalian cell includes an order of 10(4)-10(5) mRNA molecules and as many as 10(5)-10(6) ribosomes. Large-scale simultaneous mRNA translation induces correlations between the mRNA molecules, as they all compete for the finite pool of available ribosomes. This has important implications for the cell's functioning and evolution. Developing a better understanding of the intricate correlations between these simultaneous processes, rather than focusing on the translation of a single isolated transcript, should help in gaining a better understanding of mRNA translation regulation and the way elongation rates affect organismal fitness. A model of simultaneous translation is specifically important when dealing with highly expressed genes, as these consume more resources. In addition, such a model can lead to more accurate predictions that are needed in the interconnection of translational modules in synthetic biology. We develop and analyse a general dynamical model for large-scale simultaneous mRNA translation and competition for ribosomes. This is based on combining several ribosome flow models (RFMs) interconnected via a pool of free ribosomes. We use this model to explore the interactions between the various mRNA molecules and ribosomes at steady state. We show that the compound system always converges to a steady state and that it always entrains or phase locks to periodically time-varying transition rates in any of the mRNA molecules. We then study the effect of changing the transition rates in one mRNA molecule on the steady-state translation rates of the other mRNAs that results from the competition for ribosomes. We show that increasing any of the codon translation rates in a specific mRNA molecule yields a local effect, an increase in the translation rate of this mRNA, and also a global effect, the translation rates in the other mRNA molecules all increase or all decrease. These results suggest that the effect of codon decoding rates of endogenous and heterologous mRNAs on protein production is more complicated than previously thought. In addition, we show that increasing the length of an mRNA molecule decreases the production rate of all the mRNAs.

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.