A mean-field approach for modeling the propagation of perturbations in biochemical reaction networks
Introduction
Systems biology has emerged in recent years as a powerful method for quantitative modeling of biological phenomena. The approach has gained attention because it is generic and it enables the integration of experimental data with mathematical and computational modeling to understand and make predictions about complex biological processes. Using a systems biology approach, one can model virtually any type of interacting system, including protein and gene networks, metabolic networks, cellular interactions, and tissue and organism-level interactions. The effectiveness of the method has contributed to its popularity for addressing questions in a number of different biological settings, including in molecular pathways underlying cancers and other diseases Barillot et al. (2012); Chuang et al. (2010); Ebhardt et al. (2015); Faratian et al. (2009); Greenblum et al. (2012); Hornberg et al. (2006); Kreeger and Lauffenburger (2010); Loscalzo and Barabasi (2011); Mani et al. (2008); Nicholson (2006); Parikshak et al. (2015); Perou and Børresen-Dale (2011); Yarden and Pines (2012), in drug discovery and target identification Arrell and Terzic (2010); Azmi et al. (2010); Berg et al. (2005); Bugrim et al. (2004); Butcher (2005); Butcher et al. (2004); Keskin et al. (2007); Zhu et al. (2008), and in overcoming drug resistance Galluzzi et al. (2014); Wang et al. (2013).
However, one limitation that arises in systems biology approaches, particularly in the context of molecular pathways and biochemical reaction networks, such as gene and protein regulatory networks, is that they often include several kinetic parameters that cannot be determined independently of experimental time series data. The limited availability of quantitative experimental data then leads to a trade-off in the development of the model. On one hand, the goal of a systems biology approach is to develop a mathematical model that is complex enough to capture the observed biological phenomena and to answer questions related to the corresponding protein or gene network. On the other hand, increasing the complexity of the model without sufficient experimental data leads to a model with poor predictive power Ay and Arnosti (2011). In practice, fitting the outputs of the systems biology model to a limited experimental data set can lead to non-unique parameter sets or subsets of parameters that give fits of similar quality. This results in difficulties in interpreting the meanings of the parameters in the underlying molecular pathway or biochemical reaction network.
There are other methods for modeling the temporal evolution of gene regulatory networks, such as boolean models, that require only qualitative experimental data. However, while these types of models are easily extendable to large-scale systems, they tend to produce only qualitative results Ay and Arnosti (2011). On the other hand, several quantitative model reduction techniques have been developed to overcome the challenges posed by large-scale systems biology models. Such methods were recently reviewed in Snowden et al. (2017) and can be roughly categorized into timescale exploitation approaches Briggs and Haldane (1925); Choi et al. (2008); Debussche and Temam (1991); Härdin et al. (2009); Klonowski (1983); Kokotović (1984); Lam (1985); Maas and Pope (1992); Michaelis and Menten (1913); Noel et al. (2012); Prescott and Papachristodoulou (2014); Schneider and Wilhelm (2000); Surovtsova et al. (2009); Tikhonov (1952); Vejchodskỳ et al. (2014); West et al. (2015); Zobeley et al. (2005), reduction via sensitivity analysis Apri et al. (2012); Maurya et al. (2005a); Turanyi et al. (1989); Zhang and Goutsias (2010); Zi (2011), optimisation methods Anderson et al. (2011); Danø et al. (2006); Hangos et al. (2013); Maurya, Bornheimer, Venkatasubramanian, Subramaniam, 2009, Maurya, Scott, Venkatasubramanian, Subramaniam, 2005; Prescott and Papachristodoulou (2012); Taylor et al. (2008), lumping Dokoumetzidis and Aarons (2009); Koschorreck et al. (2007); Kuo and Wei (1969); Li and Rabitz (1990); Li, Rabitz, Tóth, 1994, Li, Tomlin, Rabitz, Tóth, 1994; Sunnåker, Cedersund, Jirstrand, 2011, Sunnåker, Schmidt, Jirstrand, Cedersund, 2010; Tomlin et al. (1994); Wei and Kuo (1969), and singular value decomposition-based approaches Hahn and Edgar (2000); Hardin and van Schuppen (2006); Lall et al. (2002); Liebermeister et al. (2005); Moore (1981); Sootla and Anderson (2014). Other approaches that do not strictly fit into these categories have also been developed for biochemical reaction networks, for example, based on iterative methods Apri et al. (2014); Maiwald et al. (2016); Quaiser et al. (2011); Rao et al. (2013) and differential geometry Transtrum, Qiu, 2014, Transtrum, Qiu, 2016. Each approach has its advantages, limitations, and distinct ranges of applicability Snowden et al. (2017). The optimal technique for a particular system depends on the context of the problem at hand, including the nature of the full-scale model and available experimental data and the goals of the modeller.
Model reduction approaches have also been developed that draw inspiration from mean-field theory in physics and probability theory Chaikin et al. (1995); Curie (1895); Weiss (1907), in which high-dimensional stochastic models are reduced to a simpler approximate model by averaging over the degrees of freedom in the system. This enables the reduction of a many-body problem to a one-body problem by approximating individual interactions by a single averaged effect. In the context of gene regulatory networks, mean-field approaches have been developed to construct simplified models of gene transcription Andrecut and Kauffman (2006); Cao and Grima (2018); Sturrock et al. (2015). However, in the case of Andrecut and Kauffman (2006) and Sturrock et al. (2015), the resulting models take the form of systems of coupled ordinary differential equations with several kinetic parameters and are not particularly suitable for modeling large-scale networks. Indeed, the existing mean-field approaches were applied to the study of small-scale networks with a limited number of feedback loops and are well-suited for modeling gene regulatory networks that only include protein-promoter type interactions.
While several approaches have been developed to construct simplified models of gene regulatory networks, another common issue that arises in the mathematical modeling of genetic pathways is that, in many cases, transcript levels alone are not sufficient to predict protein levels Li and Biggin (2015); Liu et al. (2016). This is because RNA expression often does not correlate with protein expression for a number of reasons, including the availability of amino acids, the formation of higher order protein complexes, and post-translational modifications, such as phosphorylation and proteolytic events. Despite these biological constraints, often it is dubiously assumed, for simplicity, that the rate of protein translation is directly proportional to RNA concentration. Indeed, this is the case for the mean-field approaches discussed above.
In this work, we introduce a simplistic mean-field approach to model the temporal evolution of protein expression levels which partially overcomes the limitations and challenges faced by other modeling techniques. The method relies on the decomposition of a complicated molecular pathway into a simplified uncoupled model of the net effects exerted by the interaction network on a target protein. The method can be easily scaled to large interaction networks and can be generalized to encompass multiple target proteins and arbitrary types of perturbations to the network. Intriguingly, the dominant effects identified by the mean-field model, using only the target protein expression levels, are supported by RNA sequencing data of the regulatory proteins in the pathway, despite the model having no a priori knowledge of this data. Thus this approach offers interpretability to the modeling results and provides a stepping stone toward the integration of protein level data with gene expression data.
The manuscript is organized as follows. In Section 2, we discuss the relevant biological background and protein interactions for the examples considered in this work. In Section 3, we describe the experimental methods and the modeling approach, then present the results and discussion in Section 4. Finally, in Section 5 we give some concluding remarks.
Section snippets
Biology background
This work is motivated by the development of a large-scale mathematical model to study the effects of combination therapy on venetoclax-resistant acute myeloid leukemia (AML) cells Przedborski et al. (2021); Sharon et al. (2019). A genome-wide CRISPR-Cas9 knockout screen revealed that genes involved in mitochondrial translation are potential therapeutic targets for restoring sensitivity to venetoclax in resistant AML cells Sharon et al. (2019). The administration of antibiotics that inhibit
Methods
We developed a mean-field approach to overcome some of the limitations of standard modeling methods in the face of complex molecular pathways and limited experimental data. The approach developed here captures quantitative experimental trends as well as enables the interpretation of the pathway dynamics in response to a network perturbation, which, in this work, is due to the administration of therapeutic agents. The details of the available experimental data and the modeling approach developed
Results
Using the integrative methods described in Section 3, we examined the ability of the mathematical model to capture meaningful trends in the experimental data set for in vitro experiments of molm-13 R2 cells Sharon et al. (2019) treated with venetoclax monotherapy, tedizolid monotherapy, and venetoclax/tedizolid combination therapy. First, we show that the model for cellular drug uptake indeed captures the previously observed in vivo drug uptake curves. Then we illustrate that the long-term in
Discussion
In this work, we constructed an analytical and computational framework for modeling the temporal propagation of perturbations through a gene regulatory network. The framework relies on first identifying the possible regulatory effects that can occur for a protein of interest, referred to in this work as the target protein. These effects are subsequently implemented into a mean-field mathematical model to capture the dominant regulatory behaviour leading to the time evolution of the target
Acknowledgements
M.K. and M.P. acknowledge the financial support from the Canadian Institutes of Health Research (CIHR).
References (172)
- et al.
Model decomposition and reduction tools for large-scale networks in systems biology
Automatica
(2011) - et al.
Complexity reduction preserving dynamical behavior of biochemical networks
J. Theor. Biol.
(2012) - et al.
Endoplasmic reticulum: er stress regulates mitochondrial bioenergetics
The international journal of biochemistry & cell biology
(2012) - et al.
Early prediction of drug metabolism and toxicity: systems biology approach and modeling
Drug Discov. Today
(2004) - et al.
E2f4/5 And p107 as smad cofactors linking the tgf receptor to c-myc repression
Cell
(2002) - et al.
P38mapk: stress responses from molecular mechanisms to therapeutics
Trends Mol Med
(2009) - et al.
The amino terminal and e2f interaction domains are critical for c/ebp-mediated induction of granulopoietic development of hematopoietic cells
Blood
(2003) - et al.
Inertial manifolds and slow manifolds
Appl Math Lett
(1991) - et al.
Reduction of nonlinear models using balancing of empirical gramians and galerkin projections
Proceedings of the 2000 American Control Conference. ACC (IEEE Cat. No. 00CH36334)
(2000) - et al.
The c/ebp homologous protein (chop) transcription factor functions in endoplasmic reticulum stress-induced apoptosis and microbial infection
Front Immunol
(2018)
Mathematical modeling in systems biology: An introduction
Cd93 marks a non-quiescent human leukemia stem cell population and is required for development of mll-rearranged acute myeloid leukemia
Cell Stem Cell
Cyclophilin b is involved in p300-mediated degradation of chop in tumor cell adaptation to hypoxia
Cell Death Differ.
The mitochondrial unfolded protein response, a conserved stress response pathway with implications in health and disease
J. Exp. Biol.
Applications of singular perturbation techniques to control problems
SIAM Rev.
Lumping analysis in monomolecular reaction systems. analysis of approximately lumpable system
Industrial & Engineering chemistry fundamentals
A general analysis of exact nonlinear lumping in chemical kinetics
Chem Eng Sci
A general analysis of approximate nonlinear lumping in chemical kinetics. i. unconstrained lumping
J Chem Phys
Tedizolid for the management of human infections: in vitro characteristics
Clinical infectious diseases
Driving the model to its limit: profile likelihood based model reduction
PLoS ONE
Bcl-2 regulates reactive oxygen species signaling and a redox-sensitive mitochondrial proton leak in mouse pancreatic -cells
Endocrinology
Isolation of hematopoietic stem cells and the effect of cd38 expression during the early erythroid progenitor cell development process
Oncol Lett
An introduction to systems biology: Design principles of biological circuits
Differential expression analysis for sequence count data
Genome Biol.
Mean-field model of genetic regulatory networks
New J Phys
Identifying optimal models to represent biochemical systems
PLoS ONE
Network systems biology for drug discovery
Clinical Pharmacology & Therapeutics
Mathematical modeling of gene expression: a guide for the perplexed biologist
Crit. Rev. Biochem. Mol. Biol.
Proof of concept: network and systems biology approaches aid in the discovery of potent anticancer drug combinations
Mol. Cancer Ther.
Overexpression of the atp binding cassette gene abca1 determines resistance to curcumin in m14 melanoma cells
Mol. Cancer
Gene expression regulation by upstream open reading frames and human disease
PLoS Genet.
Computational systems biology of cancer
Biological complexity and drug discovery: a practical systems biology approach
IEE Proceedings-Systems Biology
A note on the kinetics of enzyme action
Biochem. J
Can cell systems biology rescue drug discovery?
Nat. Rev. Drug Discovery
Systems biology in drug discovery
Nat. Biotechnol.
Linear mapping approximation of gene regulatory networks with stochastic dynamics
Nat Commun
Principles of condensed matter physics
A role for p38 mitogen-activated protein kinase and c-myc in endothelin-dependent rat aortic smooth muscle cell proliferation
Hypertension
Differential regulation of chop translation by phosphorylated eif4e under stress conditions
Nucleic Acids Res.
Differential regulation of chop translation by phosphorylated eif4e under stress conditions
Nucleic Acids Res.
Rapid method to determine intracellular drug concentrations in cellular uptake assays: application to metformin in organic cation transporter 1–transfected human embryonic kidney 293 cells
Drug Metab. Dispos.
New time-scale criteria for model simplification of bio-reaction systems
BMC Bioinformatics
A decade of systems biology
Annu. Rev. Cell Dev. Biol.
An integrated cellular and sub-cellular model of cancer chemotherapy and therapies that target cell survival.
Mathematical biosciences and engineering: MBE
The integrated stress response: from mechanism to disease
Science
Mechanisms and functions of p38 mapk signalling
Biochem. J
Propriétés magnétiques des corps a diverses températures
Phosphorylation of chop (c/ebp homologous protein) by the amp-activated protein kinase alpha 1 in macrophages promotes chop degradation and reduces injury-induced neointimal disruption in vivo
Circ. Res.
Myc, metabolism, cell growth, and tumorigenesis
Cold Spring Harb Perspect Med
Cited by (1)
An integrative systems biology approach to overcome venetoclax resistance in acute myeloid leukemia
2022, PLoS Computational Biology