The origins of time-delay in template biopolymerization processes

Analytic delay distributions for a family of gene transcription models

Models intended to describe the time evolution of a gene network must somehow include transcription, the DNA-templated synthesis of RNA, and translation, the RNA-templated synthesis of proteins. In eukaryotes, the DNA template for transcription can be very long, often consisting of tens of thousands of nucleotides, and lengthy pauses may punctuate this process. Accordingly, transcription can last for many minutes, in some cases hours. There is a long history of introducing delays in gene expression models to take the transcription and translation times into account. Here we study a family of detailed transcription models that includes initiation, elongation, and termination reactions. We establish a framework for computing the distribution of transcription times, and work out these distributions for some typical cases. For elongation, a fixed delay is a good model provided elongation is fast compared to initiation and termination, and there are no sites where long pauses occur. The initiation and termination phases of the model then generate a nontrivial delay distribution, and elongation shifts this distribution by an amount corresponding to the elongation delay. When initiation and termination are relatively fast, the distribution of elongation times can be approximated by a Gaussian. A convolution of this Gaussian with the initiation and termination time distributions gives another analytic approximation to the transcription time distribution. If there are long pauses during elongation, because of the modularity of the family of models considered, the elongation phase can be partitioned into reactions generating a simple delay (elongation through regions where there are no long pauses), and reactions whose distribution of waiting times must be considered explicitly (initiation, termination, and motion through regions where long pauses are likely). In these cases, the distribution of transcription times again involves a nontrivial part and a shift due to fast elongation processes.

A Nonlinear Delay Model for Metabolic Oscillations in Yeast Cells

We introduce two time-delay models of metabolic oscillations in yeast cells. Our model tests a hypothesis that the oscillations occur as multiple pathways share a limited resource which we equate to the number of available ribosomes. We initially explore a single-protein model with a constraint equation governing the total resource available to the cell. The model is then extended to include three proteins that share a resource pool. Three approaches are considered at constant delay to numerically detect oscillations. First, we use a spectral element method to approximate the system as a discrete map and evaluate the stability of the linearized system about its equilibria by examining its eigenvalues. For the second method, we plot amplitudes of the simulation trajectories in 2D projections of the parameter space. We use a history function that is consistent with published experimental results to obtain metabolic oscillations. Finally, the spectral element method is used to convert the system to a boundary value problem whose solutions correspond to approximate periodic solutions of the system. Our results show that certain combinations of total resource available and the time delay, lead to oscillations. We observe that an oscillation region in the parameter space is between regions admitting steady states that correspond to zero and constant production. Similar behavior is found with the three-protein model where all proteins require the same production time. However, a shift in the protein production rates peaks occurs for low available resource suggesting that our model captures the shared resource pool dynamics.

Ribosome Abundance Control in Prokaryotes

Cell growth is an essential phenotype of any unicellular organism and it crucially depends on precise control of protein synthesis. We construct a model of the feedback mechanisms that regulate abundance of ribosomes in E. coli, a prototypical prokaryotic organism. Since ribosomes are needed to produce more ribosomes, the model includes a positive feedback loop central to the control of cell growth. Our analysis of the model shows that there can be only two coexisting equilibrium states across all 23 parameters. This precludes the existence of hysteresis, suggesting that the ribosome abundance changes continuously with parameters. These states are related by a transcritical bifurcation, and we provide an analytic formula for parameters that admit either state.

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.

Trade-offs among transcription elongation rate, number, and duration of ubiquitous pauses on highly transcribed bacterial genes

In this paper, we study the limitations imposed on the transcription process by the presence of short ubiquitous pauses and crowding. These effects are especially pronounced in highly transcribed genes such as ribosomal genes (rrn) in fast growing bacteria. Our model indicates that the quantity and duration of pauses reported for protein-coding genes is incompatible with the average elongation rate observed in rrn genes. When maximal elongation rate is high, pause-induced traffic jams occur, increasing promoter occlusion, thereby lowering the initiation rate. This lowers average transcription rate and increases average transcription time. Increasing maximal elongation rate in the model is insufficient to match the experimentally observed average elongation rate in rrn genes. This suggests that there may be rrn-specific modifications to RNAP, which then experience fewer pauses, or pauses of shorter duration than those in protein-coding genes. We identify model parameter triples (maximal elongation rate, mean pause duration time, number of pauses) which are compatible with experimentally observed elongation rates. Average transcription time and average transcription rate are the model outputs investigated as proxies for cell fitness. These fitness functions are optimized for different parameter choices, opening up a possibility of differential control of these aspects of the elongation process, with potential evolutionary consequences. As an example, a gene's average transcription time may be crucial to fitness when the surrounding medium is prone to abrupt changes. This paper demonstrates that a functional relationship among the model parameters can be estimated using a standard statistical analysis, and this functional relationship describes the various trade-offs that must be made in order for the gene to control the elongation process and achieve a desired average transcription time. It also demonstrates the robustness of the system when a range of maximal elongation rates can be balanced with transcriptional pause data in order to maintain a desired fitness.

Influence of Delayed Conductance on Neuronal Synchronization

In the brain, the excitation-inhibition balance prevents abnormal synchronous behavior. However, known synaptic conductance intensity can be insufficient to account for the undesired synchronization. Due to this fact, we consider time delay in excitatory and inhibitory conductances and study its effect on the neuronal synchronization. In this work, we build a neuronal network composed of adaptive integrate-and-fire neurons coupled by means of delayed conductances. We observe that the time delay in the excitatory and inhibitory conductivities can alter both the state of the collective behavior (synchronous or desynchronous) and its type (spike or burst). For the weak coupling regime, we find that synchronization appears associated with neurons behaving with extremes highest and lowest mean firing frequency, in contrast to when desynchronization is present when neurons do not exhibit extreme values for the firing frequency. Synchronization can also be characterized by neurons presenting either the highest or the lowest levels in the mean synaptic current. For the strong coupling, synchronous burst activities can occur for delays in the inhibitory conductivity. For approximately equal-length delays in the excitatory and inhibitory conductances, desynchronous spikes activities are identified for both weak and strong coupling regimes. Therefore, our results show that not only the conductance intensity, but also short delays in the inhibitory conductance are relevant to avoid abnormal neuronal synchronization.

Kinetic Modeling of Virus Growth in Cells

When a virus infects a host cell, it hijacks the biosynthetic capacity of the cell to produce virus progeny, a process that may take less than an hour or more than a week. The overall time required for a virus to reproduce depends collectively on the rates of multiple steps in the infection process, including initial binding of the virus particle to the surface of the cell, virus internalization and release of the viral genome within the cell, decoding of the genome to make viral proteins, replication of the genome, assembly of progeny virus particles, and release of these particles into the extracellular environment. For a large number of virus types, much has been learned about the molecular mechanisms and rates of the various steps. However, in only relatively few cases during the last 50 years has an attempt been made-using mathematical modeling-to account for how the different steps contribute to the overall timing and productivity of the infection cycle in a cell. Here we review the initial case studies, which include studies of the one-step growth behavior of viruses that infect bacteria (Qβ, T7, and M13), human immunodeficiency virus, influenza A virus, poliovirus, vesicular stomatitis virus, baculovirus, hepatitis B and C viruses, and herpes simplex virus. Further, we consider how such models enable one to explore how cellular resources are utilized and how antiviral strategies might be designed to resist escape. Finally, we highlight challenges and opportunities at the frontiers of cell-level modeling of virus infections.

Formulation of the protein synthesis rate with sequence information

Translation is a central biological process by which proteins are synthesized from genetic information contained within mRNAs. Here, we investigate the kinetics of translation at the molecular level by a stochastic simulation model. The model explicitly includes RNA sequences, ribosome dynamics, the tRNA pool and biochemical reactions involved in the translation elongation. The results show that the translation efficiency is mainly limited by the available ribosome number, translation initiation and the translation elongation time. The elongation time is a log-normal distribution, with the mean and variance determined by the codon saturation and the process of aa-tRNA selection at each codon binding site. Moreover, our simulations show that the translation accuracy exponentially decreases with the sequence length. These results suggest that aa-tRNA competition is crucial for both translation elongation, translation efficiency and the accuracy, which in turn determined the effective protein production rate of correct proteins. Our results improve the dynamical equation of protein production with a delay differential equation that is dependent on sequence information through both the effective production rate and the distribution of elongation time.

Parametric Sensitivity Analysis of Oscillatory Delay Systems with an Application to Gene Regulation

A parametric sensitivity analysis for periodic solutions of delay-differential equations is developed. Because phase shifts cause the sensitivity coefficients of a periodic orbit to diverge, we focus on sensitivities of the extrema, from which amplitude sensitivities are computed, and of the period. Delay-differential equations are often used to model gene expression networks. In these models, the parametric sensitivities of a particular genotype define the local geometry of the evolutionary landscape. Thus, sensitivities can be used to investigate directions of gradual evolutionary change. An oscillatory protein synthesis model whose properties are modulated by RNA interference is used as an example. This model consists of a set of coupled delay-differential equations involving three delays. Sensitivity analyses are carried out at several operating points. Comments on the evolutionary implications of the results are offered.

Elements of biological oscillations in time and space

Oscillations in time and space are ubiquitous in nature and play critical roles in dynamic cellular processes. Although the molecular mechanisms underlying the generation of the dynamics are diverse, several distinct regulatory elements have been recognized as being critical in producing and modulating oscillatory dynamics. These include negative and positive feedback, time delay, nonlinearity in regulation, and random fluctuations ('noise'). Here we discuss the specific roles of these five elements in promoting or attenuating oscillatory dynamics, by drawing on insights from quantitative analyses of natural or synthetic biological networks.

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.

Stochastic modeling of biochemical systems with multistep reactions using state-dependent time delay

To deal with the growing scale of molecular systems, sophisticated modelling techniques have been designed in recent years to reduce the complexity of mathematical models. Among them, a widely used approach is delayed reaction for simplifying multistep reactions. However, recent research results suggest that a delayed reaction with constant time delay is unable to describe multistep reactions accurately. To address this issue, we propose a novel approach using state-dependent time delay to approximate multistep reactions. We first use stochastic simulations to calculate time delay arising from multistep reactions exactly. Then we design algorithms to calculate time delay based on system dynamics precisely. To demonstrate the power of proposed method, two processes of mRNA degradation are used to investigate the function of time delay in determining system dynamics. In addition, a multistep pathway of metabolic synthesis is used to explore the potential of the proposed method to simplify multistep reactions with nonlinear reaction rates. Simulation results suggest that the state-dependent time delay is a promising and accurate approach to reduce model complexity and decrease the number of unknown parameters in the models.

High-order dynamic Bayesian Network learning with hidden common causes for causal gene regulatory network

BACKGROUND

Inferring gene regulatory network (GRN) has been an important topic in Bioinformatics. Many computational methods infer the GRN from high-throughput expression data. Due to the presence of time delays in the regulatory relationships, High-Order Dynamic Bayesian Network (HO-DBN) is a good model of GRN. However, previous GRN inference methods assume causal sufficiency, i.e. no unobserved common cause. This assumption is convenient but unrealistic, because it is possible that relevant factors have not even been conceived of and therefore un-measured. Therefore an inference method that also handles hidden common cause(s) is highly desirable. Also, previous methods for discovering hidden common causes either do not handle multi-step time delays or restrict that the parents of hidden common causes are not observed genes.

RESULTS

We have developed a discrete HO-DBN learning algorithm that can infer also hidden common cause(s) from discrete time series expression data, with some assumptions on the conditional distribution, but is less restrictive than previous methods. We assume that each hidden variable has only observed variables as children and parents, with at least two children and possibly no parents. We also make the simplifying assumption that children of hidden variable(s) are not linked to each other. Moreover, our proposed algorithm can also utilize multiple short time series (not necessarily of the same length), as long time series are difficult to obtain.

CONCLUSIONS

We have performed extensive experiments using synthetic data on GRNs of size up to 100, with up to 10 hidden nodes. Experiment results show that our proposed algorithm can recover the causal GRNs adequately given the incomplete data. Using the limited real expression data and small subnetworks of the YEASTRACT network, we have also demonstrated the potential of our algorithm on real data, though more time series expression data is needed.

Time Delayed Causal Gene Regulatory Network Inference with Hidden Common Causes

Inferring the gene regulatory network (GRN) is crucial to understanding the working of the cell. Many computational methods attempt to infer the GRN from time series expression data, instead of through expensive and time-consuming experiments. However, existing methods make the convenient but unrealistic assumption of causal sufficiency, i.e. all the relevant factors in the causal network have been observed and there are no unobserved common cause. In principle, in the real world, it is impossible to be certain that all relevant factors or common causes have been observed, because some factors may not have been conceived of, and therefore are impossible to measure. In view of this, we have developed a novel algorithm named HCC-CLINDE to infer an GRN from time series data allowing the presence of hidden common cause(s). We assume there is a sparse causal graph (possibly with cycles) of interest, where the variables are continuous and each causal link has a delay (possibly more than one time step). A small but unknown number of variables are not observed. Each unobserved variable has only observed variables as children and parents, with at least two children, and the children are not linked to each other. Since it is difficult to obtain very long time series, our algorithm is also capable of utilizing multiple short time series, which is more realistic. To our knowledge, our algorithm is far less restrictive than previous works. We have performed extensive experiments using synthetic data on GRNs of size up to 100, with up to 10 hidden nodes. The results show that our algorithm can adequately recover the true causal GRN and is robust to slight deviation from Gaussian distribution in the error terms. We have also demonstrated the potential of our algorithm on small YEASTRACT subnetworks using limited real data.

Inferring Time-Delayed Causal Gene Network Using Time-Series Expression Data

Inferring gene regulatory network (GRN) from the microarray expression data is an important problem in Bioinformatics, because knowing the GRN is an essential first step in understanding the inner workings of the cell and the related diseases. Time delays exist in the regulatory effects from one gene to another due to the time needed for transcription, translation, and to accumulate a sufficient number of needed proteins. Also, it is known that the delays are important for oscillatory phenomenon. Therefore, it is crucial to develop a causal gene network model, preferably as a function of time. In this paper, we propose an algorithm CLINDE to infer causal directed links in GRN with time delays and regulatory effects in the links from time-series microarray gene expression data. It is one of the most comprehensive in terms of features compared to the state-of-the-art discrete gene network models. We have tested CLINDE on synthetic data, the in vivo IRMA (On and Off) datasets and the [1] yeast expression data validated using KEGG pathways. Results show that CLINDE can effectively recover the links, the time delays and the regulatory effects in the synthetic data, and outperforms other algorithms in the IRMA in vivo datasets.

Noise analysis of genome-scale protein synthesis using a discrete computational model of translation

Noise in genetic networks has been the subject of extensive experimental and computational studies. However, very few of these studies have considered noise properties using mechanistic models that account for the discrete movement of ribosomes and RNA polymerases along their corresponding templates (messenger RNA (mRNA) and DNA). The large size of these systems, which scales with the number of genes, mRNA copies, codons per mRNA, and ribosomes, is responsible for some of the challenges. Additionally, one should be able to describe the dynamics of ribosome exchange between the free ribosome pool and those bound to mRNAs, as well as how mRNA species compete for ribosomes. We developed an efficient algorithm for stochastic simulations that addresses these issues and used it to study the contribution and trade-offs of noise to translation properties (rates, time delays, and rate-limiting steps). The algorithm scales linearly with the number of mRNA copies, which allowed us to study the importance of genome-scale competition between mRNAs for the same ribosomes. We determined that noise is minimized under conditions maximizing the specific synthesis rate. Moreover, sensitivity analysis of the stochastic system revealed the importance of the elongation rate in the resultant noise, whereas the translation initiation rate constant was more closely related to the average protein synthesis rate. We observed significant differences between our results and the noise properties of the most commonly used translation models. Overall, our studies demonstrate that the use of full mechanistic models is essential for the study of noise in translation and transcription.

Hypoxia augments outgrowth endothelial cell (OEC) sprouting and directed migration in response to sphingosine-1-phosphate (S1P)

Therapeutic angiogenesis provides a promising approach to treat ischemic cardiovascular diseases through the delivery of proangiogenic cells and/or molecules. Outgrowth endothelial cells (OECs) are vascular progenitor cells that are especially suited for therapeutic strategies given their ease of noninvasive isolation from umbilical cord or adult peripheral blood and their potent ability to enhance tissue neovascularization. These cells are recruited to sites of vascular injury or tissue ischemia and directly incorporate within native vascular endothelium to participate in neovessel formation. A better understanding of how OEC activity may be boosted under hypoxia with external stimulation by proangiogenic molecules remains a challenge to improving their therapeutic potential. While vascular endothelial growth factor (VEGF) is widely established as a critical factor for initiating angiogenesis, sphingosine-1-phosphate (S1P), a bioactive lysophospholipid, has recently gained great enthusiasm as a potential mediator in neovascularization strategies. This study tests the hypothesis that hypoxia and the presence of VEGF impact the angiogenic response of OECs to S1P stimulation in vitro. We found that hypoxia altered the dynamically regulated S1P receptor 1 (S1PR1) expression on OECs in the presence of S1P (1.0 μM) and/or VEGF (1.3 nM). The combined stimuli of S1P and VEGF together promoted OEC angiogenic activity as assessed by proliferation, wound healing, 3D sprouting, and directed migration under both normoxia and hypoxia. Hypoxia substantially augmented the response to S1P alone, resulting in ~6.5-fold and ~25-fold increases in sprouting and directed migration, respectively. Overall, this report highlights the importance of establishing hypoxic conditions in vitro when studying ischemia-related angiogenic strategies employing vascular progenitor cells.

The utility of simple mathematical models in understanding gene regulatory dynamics

In this review, we survey work that has been carried out in the attempts of biomathematicians to understand the dynamic behaviour of simple bacterial operons starting with the initial work of the 1960's. We concentrate on the simplest of situations, discussing both repressible and inducible systems and then turning to concrete examples related to the biology of the lactose and tryptophan operons. We conclude with a brief discussion of the role of both extrinsic noise and so-called intrinsic noise in the form of translational and/or transcriptional bursting.

Continuous neural identifier for uncertain nonlinear systems with time delays in the input signal

Time-delay systems have been successfully used to represent the complexity of some dynamic systems. Time-delay is often used for modeling many real systems. Among others, biological and chemical plants have been described using time-delay terms with better results than those models that have not consider them. However, getting those models represented a challenge and sometimes the results were not so satisfactory. Non-parametric modeling offered an alternative to obtain suitable and usable models. Continuous neural networks (CNN) have been considered as a real alternative to provide models over uncertain non-parametric systems. This article introduces the design of a specific class of non-parametric model for uncertain time-delay system based on CNN considering the so-called delayed learning laws analysis. The convergence analysis as well as the learning laws were produced by means of a Lyapunov-Krasovskii functional. Three examples were developed to demonstrate the effectiveness of the modeling process forced by the identifier proposed in this study. The first example was a simple nonlinear model used as benchmark example. The second example regarded the human immunodeficiency virus dynamic behavior is used to show the performance of the suggested non-parametric identifier based on CNN for no fictitious neither academic models. Finally, a third example describing the evolution of hepatitis B virus served to test the identifier presented in this study and was also useful to provide evidence of its superior performance against a non-delayed identifier based on CNN.

DnaA dynamics could be linked with fitness cost in bacteria

Fitness cost is the measure of the metabolic burden of unneeded gene expression. It is defined as the lag in bacterial cells growth harboring unneeded genes relative to unburdened cells. Separate cells can concurrently adapt to the burden, demonstrating a decrease in or even a disappearance of the lag. The precise mechanisms of this adaptation are not clearly understood. One possibility is that an increased amount of free ribosomes "absorb" the unnecessary burden. In this work, the mechanism by which an increased concentration of ribosomes could result in faster growth and mask the unneeded gene expression burden is discussed. The initiation time of chromosome replication by the initiator protein DnaA, for which the accumulation speed depends on the ribosomes amount, is taken into account.

Mechanistically consistent reduced models of synthetic gene networks

Designing genetic networks with desired functionalities requires an accurate mathematical framework that accounts for the essential mechanistic details of the system. Here, we formulate a time-delay model of protein translation and mRNA degradation by systematically reducing a detailed mechanistic model that explicitly accounts for the ribosomal dynamics and the cleaving of mRNA by endonucleases. We exploit various technical and conceptual advantages that our time-delay model offers over the mechanistic model to probe the behavior of a self-repressing gene over wide regions of parameter space. We show that a heuristic time-delay model of protein synthesis of a commonly used form yields a notably different prediction for the parameter region where sustained oscillations occur. This suggests that such heuristics can lead to erroneous results. The functional forms that arise from our systematic reduction can be used for every system that involves transcription and translation and they could replace the commonly used heuristic time-delay models for these processes. The results from our analysis have important implications for the design of synthetic gene networks and stress that such design must be guided by a combination of heuristic models and mechanistic models that include all relevant details of the process.

Stochastic modelling of biochemical systems of multi-step reactions using a simplified two-variable model

Background

A fundamental issue in systems biology is how to design simplified mathematical models for describing the dynamics of complex biochemical reaction systems. Among them, a key question is how to use simplified reactions to describe the chemical events of multi-step reactions that are ubiquitous in biochemistry and biophysics. To address this issue, a widely used approach in literature is to use one-step reaction to represent the multi-step chemical events. In recent years, a number of modelling methods have been designed to improve the accuracy of the one-step reaction method, including the use of reactions with time delay. However, our recent research results suggested that there are still deviations between the dynamics of delayed reactions and that of the multi-step reactions. Therefore, more sophisticated modelling methods are needed to accurately describe the complex biological systems in an efficient way.

Results

This work designs a two-variable model to simplify chemical events of multi-step reactions. In addition to the total molecule number of a species, we first introduce a new concept regarding the location of molecules in the multi-step reactions, which is the second variable to represent the system dynamics. Then we propose a simulation algorithm to compute the probability for the firing of the last step reaction in the multi-step events. This probability function is evaluated using a deterministic model of ordinary differential equations and a stochastic model in the framework of the stochastic simulation algorithm. The efficiency of the proposed two-variable model is demonstrated by the realization of mRNA degradation process based on the experimentally measured data.

Conclusions

Numerical results suggest that the proposed new two-variable model produces predictions that match the multi-step chemical reactions very well. The successful realization of the mRNA degradation dynamics indicates that the proposed method is a promising approach to reduce the complexity of biological systems.

A traffic flow model for bio-polymerization processes

Bio-polymerization processes like transcription and translation are central to proper function of a cell. The speed at which the bio-polymer grows is affected both by the number of pauses of elongation machinery, as well the number of bio-polymers due to crowding effects. In order to quantify these effects in fast transcribing ribosome genes, we rigorously show that a classical traffic flow model is the limit of a mean occupancy ODE model. We compare the simulation of this model to a stochastic model and evaluate the combined effect of the polymerase density and the existence of pauses on the instantaneous transcription rate of ribosomal genes.

Modeling and dynamical analysis of virus-triggered innate immune signaling pathways

The investigation of the dynamics and regulation of virus-triggered innate immune signaling pathways at a system level will enable comprehensive analysis of the complex interactions that maintain the delicate balance between resistance to infection and viral disease. In this study, we developed a delayed mathematical model to describe the virus-induced interferon (IFN) signaling process by considering several key players in the innate immune response. Using dynamic analysis and numerical simulation, we evaluated the following predictions regarding the antiviral responses: (1) When the replication ratio of virus is less than 1, the infectious virus will be eliminated by the immune system’s defenses regardless of how the time delays are changed. (2) The IFN positive feedback regulation enhances the stability of the innate immune response and causes the immune system to present the bistability phenomenon. (3) The appropriate duration of viral replication and IFN feedback processes stabilizes the innate immune response. The predictions from the model were confirmed by monitoring the virus titer and IFN expression in infected cells. The results suggest that the balance between viral replication and IFN-induced feedback regulation coordinates the dynamical behavior of virus-triggered signaling and antiviral responses. This work will help clarify the mechanisms of the virus-induced innate immune response at a system level and provide instruction for further biological experiments.

Mathematical and Computational Modelling of Ribosomal Movement and Protein Synthesis: an overview

Translation or protein synthesis consists of a complex system of chemical reactions, which ultimately result in decoding of the mRNA and the production of a protein. The complexity of this reaction system makes it difficult to quantitatively connect its input parameters (such as translation factor or ribosome concentrations, codon composition of the mRNA, or energy availability) to output parameters (such as protein synthesis rates or ribosome densities on mRNAs). Mathematical and computational models of translation have now been used for nearly five decades to investigate translation, and to shed light on the relationship between the different reactions in the system. This review gives an overview over the principal approaches used in the modelling efforts, and summarises some of the major findings that were made.