1
|
Song YM, Hong H, Kim JK. Universally valid reduction of multiscale stochastic biochemical systems using simple non-elementary propensities. PLoS Comput Biol 2021; 17:e1008952. [PMID: 34662330 PMCID: PMC8562860 DOI: 10.1371/journal.pcbi.1008952] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2021] [Revised: 11/02/2021] [Accepted: 10/04/2021] [Indexed: 11/19/2022] Open
Abstract
Biochemical systems consist of numerous elementary reactions governed by the law of mass action. However, experimentally characterizing all the elementary reactions is nearly impossible. Thus, over a century, their deterministic models that typically contain rapid reversible bindings have been simplified with non-elementary reaction functions (e.g., Michaelis-Menten and Morrison equations). Although the non-elementary reaction functions are derived by applying the quasi-steady-state approximation (QSSA) to deterministic systems, they have also been widely used to derive propensities for stochastic simulations due to computational efficiency and simplicity. However, the validity condition for this heuristic approach has not been identified even for the reversible binding between molecules, such as protein-DNA, enzyme-substrate, and receptor-ligand, which is the basis for living cells. Here, we find that the non-elementary propensities based on the deterministic total QSSA can accurately capture the stochastic dynamics of the reversible binding in general. However, serious errors occur when reactant molecules with similar levels tightly bind, unlike deterministic systems. In that case, the non-elementary propensities distort the stochastic dynamics of a bistable switch in the cell cycle and an oscillator in the circadian clock. Accordingly, we derive alternative non-elementary propensities with the stochastic low-state QSSA, developed in this study. This provides a universally valid framework for simplifying multiscale stochastic biochemical systems with rapid reversible bindings, critical for efficient stochastic simulations of cell signaling and gene regulation. To facilitate the framework, we provide a user-friendly open-source computational package, ASSISTER, that automatically performs the present framework.
Collapse
Affiliation(s)
- Yun Min Song
- Department of Mathematical Sciences, Korea Advanced Institute of Science and Technology, Daejeon, Republic of Korea
- Biomedical Mathematics Group, Institute for Basic Science, Daejeon, Republic of Korea
| | - Hyukpyo Hong
- Department of Mathematical Sciences, Korea Advanced Institute of Science and Technology, Daejeon, Republic of Korea
- Biomedical Mathematics Group, Institute for Basic Science, Daejeon, Republic of Korea
| | - Jae Kyoung Kim
- Department of Mathematical Sciences, Korea Advanced Institute of Science and Technology, Daejeon, Republic of Korea
- Biomedical Mathematics Group, Institute for Basic Science, Daejeon, Republic of Korea
| |
Collapse
|
2
|
Leviyang S, Griva I. Investigating Functional Roles for Positive Feedback and Cellular Heterogeneity in the Type I Interferon Response to Viral Infection. Viruses 2018; 10:v10100517. [PMID: 30241427 PMCID: PMC6213501 DOI: 10.3390/v10100517] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2018] [Revised: 09/16/2018] [Accepted: 09/20/2018] [Indexed: 12/20/2022] Open
Abstract
Secretion of type I interferons (IFN) by infected cells mediates protection against many viruses, but prolonged or excessive type I IFN secretion can lead to immune pathology. A proper type I IFN response must therefore maintain a balance between protection and excessive IFN secretion. It has been widely noted that the type I IFN response is driven by positive feedback and is heterogeneous, with only a fraction of infected cells upregulating IFN expression even in clonal cell lines, but the functional roles of feedback and heterogeneity in balancing protection and excessive IFN secretion are not clear. To investigate the functional roles for feedback and heterogeneity, we constructed a mathematical model coupling IFN and viral dynamics that extends existing mathematical models by accounting for feedback and heterogeneity. We fit our model to five existing datasets, reflecting different experimental systems. Fitting across datasets allowed us to compare the IFN response across the systems and suggested different signatures of feedback and heterogeneity in the different systems. Further, through numerical experiments, we generated hypotheses of functional roles for IFN feedback and heterogeneity consistent with our mathematical model. We hypothesize an inherent tradeoff in the IFN response: a positive feedback loop prevents excessive IFN secretion, but also makes the IFN response vulnerable to viral antagonism. We hypothesize that cellular heterogeneity of the IFN response functions to protect the feedback loop from viral antagonism. Verification of our hypotheses will require further experimental studies. Our work provides a basis for analyzing the type I IFN response across systems.
Collapse
Affiliation(s)
- Sivan Leviyang
- Department of Mathematics and Statistics, Georgetown University, Washington, DC 20057, USA.
| | - Igor Griva
- Department of Mathematical Sciences, George Mason University, Fairfax, VA 22030, USA.
| |
Collapse
|
3
|
Cai C, Zhou J, Sun X, Sun T, Xie W, Cui J. Integrated modeling and analysis of intracellular and intercellular mechanisms in shaping the interferon response to viral infection. PLoS One 2017; 12:e0186105. [PMID: 29020068 PMCID: PMC5636135 DOI: 10.1371/journal.pone.0186105] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2017] [Accepted: 09/25/2017] [Indexed: 12/13/2022] Open
Abstract
The interferons (IFNs) responses to viral infection are heterogeneous, while the underlying mechanisms for variability among cells are still not clear. In this study, we developed a hybrid model to systematically identify the sources of IFN induction heterogeneity. The experiment-integrated simulation demonstrated that the viral dose/type, the diversity in transcriptional factors activation and the intercellular paracrine signaling could strikingly shape the heterogeneity of IFN expression. We further determined that the IFNβ and IFNλ1 induced diverse dynamics of IFN-stimulated genes (ISGs) production. Collectively, our findings revealed the intracellular and intercellular mechanisms contributing to cell-to-cell variation in IFN induction, and further demonstrated the significant effects of IFN heterogeneity on antagonizing viruses.
Collapse
Affiliation(s)
- Chunmei Cai
- Key Laboratory of Gene Engineering of the Ministry of Education, State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-sen University, Guangzhou, P. R. China
- Zhongshan School of Medicine, Sun Yat-sen University, Guangzhou, P. R. China
| | - Jie Zhou
- Key Laboratory of Gene Engineering of the Ministry of Education, State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-sen University, Guangzhou, P. R. China
| | - Xiaoqiang Sun
- Zhongshan School of Medicine, Sun Yat-sen University, Guangzhou, P. R. China
| | | | - Weihong Xie
- Key Laboratory of Gene Engineering of the Ministry of Education, State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-sen University, Guangzhou, P. R. China
| | - Jun Cui
- Key Laboratory of Gene Engineering of the Ministry of Education, State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-sen University, Guangzhou, P. R. China
- Collaborative Innovation Center of Cancer Medicine, Sun Yat-sen University, Guangzhou, P. R. China
- * E-mail:
| |
Collapse
|
4
|
Several Indicators of Critical Transitions for Complex Diseases Based on Stochastic Analysis. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2017; 2017:7560758. [PMID: 28835768 PMCID: PMC5556999 DOI: 10.1155/2017/7560758] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/18/2017] [Revised: 06/19/2017] [Accepted: 06/19/2017] [Indexed: 12/18/2022]
Abstract
Many complex diseases (chronic disease onset, development and differentiation, self-assembly, etc.) are reminiscent of phase transitions in a dynamical system: quantitative changes accumulate largely unnoticed until a critical threshold is reached, which causes abrupt qualitative changes of the system. Understanding such nonlinear behaviors is critical to dissect the multiple genetic/environmental factors that together shape the genetic and physiological landscape underlying basic biological functions and to identify the key driving molecules. Based on stochastic differential equation (SDE) model, we theoretically derive three statistical indicators, that is, coefficient of variation (CV), transformed Pearson's correlation coefficient (TPC), and transformed probability distribution (TPD), to identify critical transitions and detect the early-warning signals of the phase transition in complex diseases. To verify the effectiveness of these early-warning indexes, we use high-throughput data for three complex diseases, including influenza caused by either H3N2 or H1N1 and acute lung injury, to extract the dynamical network biomarkers (DNBs) responsible for catastrophic transition into the disease state from predisease state. The numerical results indicate that the derived indicators provide a data-based quantitative analysis for early-warning signals for critical transitions in complex diseases or other dynamical systems.
Collapse
|
5
|
Kim JK, Sontag ED. Reduction of multiscale stochastic biochemical reaction networks using exact moment derivation. PLoS Comput Biol 2017; 13:e1005571. [PMID: 28582397 PMCID: PMC5481150 DOI: 10.1371/journal.pcbi.1005571] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2017] [Revised: 06/22/2017] [Accepted: 05/12/2017] [Indexed: 12/14/2022] Open
Abstract
Biochemical reaction networks (BRNs) in a cell frequently consist of reactions with disparate timescales. The stochastic simulations of such multiscale BRNs are prohibitively slow due to high computational cost for the simulations of fast reactions. One way to resolve this problem uses the fact that fast species regulated by fast reactions quickly equilibrate to their stationary distribution while slow species are unlikely to be changed. Thus, on a slow timescale, fast species can be replaced by their quasi-steady state (QSS): their stationary conditional expectation values for given slow species. As the QSS are determined solely by the state of slow species, such replacement leads to a reduced model, where fast species are eliminated. However, it is challenging to derive the QSS in the presence of nonlinear reactions. While various approximation schemes for the QSS have been developed, they often lead to considerable errors. Here, we propose two classes of multiscale BRNs which can be reduced by deriving an exact QSS rather than approximations. Specifically, if fast species constitute either a feedforward network or a complex balanced network, the reduced model based on the exact QSS can be derived. Such BRNs are frequently observed in a cell as the feedforward network is one of fundamental motifs of gene or protein regulatory networks. Furthermore, complex balanced networks also include various types of fast reversible bindings such as bindings between transcriptional factors and gene regulatory sites. The reduced models based on exact QSS, which can be calculated by the computational packages provided in this work, accurately approximate the slow scale dynamics of the original full model with much lower computational cost.
Collapse
Affiliation(s)
- Jae Kyoung Kim
- Department of Mathematical Sciences, Korea Advanced Institute of Science and Technology, Daejeon, Korea
- * E-mail: (JKK); , (EDS)
| | - Eduardo D. Sontag
- Department of Mathematics and Center for Quantitative Biology, Rutgers University, New Brunswick, New Jersey, United States of America
- * E-mail: (JKK); , (EDS)
| |
Collapse
|
6
|
Dalet A, Argüello RJ, Combes A, Spinelli L, Jaeger S, Fallet M, Vu Manh TP, Mendes A, Perego J, Reverendo M, Camosseto V, Dalod M, Weil T, Santos MA, Gatti E, Pierre P. Protein synthesis inhibition and GADD34 control IFN-β heterogeneous expression in response to dsRNA. EMBO J 2017; 36:761-782. [PMID: 28100675 DOI: 10.15252/embj.201695000] [Citation(s) in RCA: 53] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2016] [Revised: 12/20/2016] [Accepted: 12/22/2016] [Indexed: 02/01/2023] Open
Abstract
In innate immune responses, induction of type-I interferons (IFNs) prevents virus spreading while viral replication is delayed by protein synthesis inhibition. We asked how cells perform these apparently contradictory activities. Using single fibroblast monitoring by flow cytometry and mathematical modeling, we demonstrate that type-I IFN production is linked to cell's ability to enter dsRNA-activated PKR-dependent translational arrest and then overcome this inhibition by decreasing eIF2α phosphorylation through phosphatase 1c cofactor GADD34 (Ppp1r15a) expression. GADD34 expression, shown here to be dependent on the IRF3 transcription factor, is responsible for a biochemical cycle permitting pulse of IFN synthesis to occur in cells undergoing protein synthesis inhibition. Translation arrest is further demonstrated to be key for anti-viral response by acting synergistically with MAVS activation to amplify TBK1 signaling and IFN-β mRNA transcription, while GADD34-dependent protein synthesis recovery contributes to the heterogeneous expression of IFN observed in dsRNA-activated cells.
Collapse
Affiliation(s)
- Alexandre Dalet
- CNRS, INSERM, CIML, Aix Marseille University, Marseille, France
| | | | - Alexis Combes
- CNRS, INSERM, CIML, Aix Marseille University, Marseille, France
| | - Lionel Spinelli
- CNRS, INSERM, CIML, Aix Marseille University, Marseille, France
| | | | - Mathieu Fallet
- CNRS, INSERM, CIML, Aix Marseille University, Marseille, France
| | | | - Andreia Mendes
- CNRS, INSERM, CIML, Aix Marseille University, Marseille, France
| | - Jessica Perego
- CNRS, INSERM, CIML, Aix Marseille University, Marseille, France
| | | | - Voahirana Camosseto
- CNRS, INSERM, CIML, Aix Marseille University, Marseille, France.,International associated laboratory (LIA) CNRS "Mistra", Marseille, France
| | - Marc Dalod
- CNRS, INSERM, CIML, Aix Marseille University, Marseille, France
| | - Tobias Weil
- Institute for Research in Biomedicine - iBiMED and Aveiro Health Sciences Program, University of Aveiro, Aveiro, Portugal
| | - Manuel A Santos
- International associated laboratory (LIA) CNRS "Mistra", Marseille, France.,Institute for Research in Biomedicine - iBiMED and Aveiro Health Sciences Program, University of Aveiro, Aveiro, Portugal
| | - Evelina Gatti
- CNRS, INSERM, CIML, Aix Marseille University, Marseille, France .,International associated laboratory (LIA) CNRS "Mistra", Marseille, France.,Institute for Research in Biomedicine - iBiMED and Aveiro Health Sciences Program, University of Aveiro, Aveiro, Portugal
| | - Philippe Pierre
- CNRS, INSERM, CIML, Aix Marseille University, Marseille, France .,International associated laboratory (LIA) CNRS "Mistra", Marseille, France.,Institute for Research in Biomedicine - iBiMED and Aveiro Health Sciences Program, University of Aveiro, Aveiro, Portugal
| |
Collapse
|