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

Limit cycles in piecewise-affine gene network models with multiple interaction loops

In this paper we consider piecewise affine differential equations modeling gene networks. We work with arbitrary decay rates, and under a local hypothesis expressed as an alignment condition of successive focal points. The interaction graph of the system may be rather complex (multiple intricate loops of any sign, multiple thresholds...). Our main result is an alternative theorem showing that, if a sequence of region is periodically visited by trajectories, then under our hypotheses, there exists either a unique stable periodic solution, or the origin attracts all trajectories in this sequence of regions. This result extends greatly our previous work on a single negative feedback loop. We give several examples and simulations illustrating different cases

Periodic solutions of piecewise affine gene network models: the case of a negative feedback loop

In this paper the existence and unicity of a stable periodic orbit is proven, for a class of piecewise affine differential equations in dimension 3 or more, provided their interaction structure is a negative feedback loop. It is also shown that the same systems converge toward a unique stable equilibrium point in dimension 2. This extends a theorem of Snoussi, which showed the existence of these orbits only. The considered class of equations is usually studied as a model of gene regulatory networks. It is not assumed that all decay rates are identical, which is biologically irrelevant, but has been done in the vast majority of previous studies. Our work relies on classical results about fixed points of monotone, concave operators acting on positive variables. Moreover, the used techniques are very likely to apply in more general contexts, opening directions for future work

A simple rule for the evolution of fast dispersal at the edge of expanding populations

Evolution by natural selection is commonly perceived as a process that favors those that replicate faster to leave more offspring; nature, however, seem to abound with examples where organisms forgo some replicative potential to disperse faster. When does selection favor invasion of the fastest? Motivated by evolution experiments with swarming bacteria we searched for a simple rule. In experiments, a fast hyperswarmer mutant that pays a reproductive cost to make many copies of its flagellum invades a population of mono-flagellated bacteria by reaching the expanding population edge; a two-species mathematical model explains that invasion of the edge occurs only if the invasive species' expansion rate, v₂, which results from the combination of the species growth rate and its dispersal speed (but not its carrying capacity), exceeds the established species', v₁. The simple rule that we derive, v₂ > v₁, appears to be general: less favorable initial conditions, such as smaller initial sizes and longer distances to the population edge, delay but do not entirely prevent invasion. Despite intricacies of the swarming system, experimental tests agree well with model predictions suggesting that the general theory should apply to other expanding populations with trade-offs between growth and dispersal, including non-native invasive species and cancer metastases.First author draf

A model for selection of eyespots on butterfly wings

The development of eyespots on the wing surface of butterflies of the family Nympalidae is one of the most studied examples of biological pattern formation.However, little is known about the mechanism that determines the number and precise locations of eyespots on the wing. Eyespots develop around signaling centers, called foci, that are located equidistant from wing veins along the midline of a wing cell (an area bounded by veins). A fundamental question that remains unsolved is, why a certain wing cell develops an eyespot, while other wing cells do not. We illustrate that the key to understanding focus point selection may be in the venation system of the wing disc. Our main hypothesis is that changes in morphogen concentration along the proximal boundary veins of wing cells govern focus point selection. Based on previous studies, we focus on a spatially two-dimensional reaction-diffusion system model posed in the interior of each wing cell that describes the formation of focus points. Using finite element based numerical simulations, we demonstrate that variation in the proximal boundary condition is sufficient to robustly select whether an eyespot focus point forms in otherwise identical wing cells. We also illustrate that this behavior is robust to small perturbations in the parameters and geometry and moderate levels of noise. Hence, we suggest that an anterior-posterior pattern of morphogen concentration along the proximal vein may be the main determinant of the distribution of focus points on the wing surface. In order to complete our model, we propose a two stage reaction-diffusion system model, in which an one-dimensional surface reaction-diffusion system, posed on the proximal vein, generates the morphogen concentrations that act as non-homogeneous Dirichlet (i.e., fixed) boundary conditions for the two-dimensional reaction-diffusion model posed in the wing cells. The two-stage model appears capable of generating focus point distributions observed in nature. We therefore conclude that changes in the proximal boundary conditions are sufficient to explain the empirically observed distribution of eyespot focus points on the entire wing surface. The model predicts, subject to experimental verification, that the source strength of the activator at the proximal boundary should be lower in wing cells in which focus points form than in those that lack focus points. The model suggests that the number and locations of eyespot foci on the wing disc could be largely controlled by two kinds of gradients along two different directions, that is, the first one is the gradient in spatially varying parameters such as the reaction rate along the anterior-posterior direction on the proximal boundary of the wing cells, and the second one is the gradient in source values of the activator along the veins in the proximal-distal direction of the wing cell

Evolution at the edge of expanding populations

Predicting evolution of expanding populations is critical to control biological threats such as invasive species and cancer metastasis. Expansion is primarily driven by reproduction and dispersal, but nature abounds with examples of evolution where organisms pay a reproductive cost to disperse faster. When does selection favor this 'survival of the fastest?' We searched for a simple rule, motivated by evolution experiments where swarming bacteria evolved into an hyperswarmer mutant which disperses $\sim 100\%$ faster but pays a growth cost of $\sim 10 \%$ to make many copies of its flagellum. We analyzed a two-species model based on the Fisher equation to explain this observation: the population expansion rate ($v$) results from an interplay of growth ($r$) and dispersal ($D$) and is independent of the carrying capacity: $v=2\sqrt{rD}$. A mutant can take over the edge only if its expansion rate ($v_2$) exceeds the expansion rate of the established species ($v_1$); this simple condition ($v_2 > v_1$) determines the maximum cost in slower growth that a faster mutant can pay and still be able to take over. Numerical simulations and time-course experiments where we tracked evolution by imaging bacteria suggest that our findings are general: less favorable conditions delay but do not entirely prevent the success of the fastest. Thus, the expansion rate defines a traveling wave fitness, which could be combined with trade-offs to predict evolution of expanding populations