1
|
Khazaaleh M, Samarasinghe S, Kulasiri D. A new hierarchical approach to multi-level model abstraction for simplifying ODE models of biological networks and a case study: The G1/S Checkpoint/DNA damage signalling pathways of mammalian cell cycle. Biosystems 2021; 203:104374. [PMID: 33556446 DOI: 10.1016/j.biosystems.2021.104374] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2020] [Revised: 01/19/2021] [Accepted: 01/25/2021] [Indexed: 11/15/2022]
Abstract
Model reduction is an important topic in studies of biological systems. By reducing the complexity of large models through multi-level models while keeping the essence (biological meaning) of the model, model reduction can help answer many important questions about these systems. In this paper, we present a new reduction method based on hierarchical representation and a lumping approach. We used G1/S checkpoint pathway represented in Ordinary Differential Equations (ODE) in Iwamoto et al. (2011) as a case study to present this reduction method. The approach consists of two parts; the first part represents the biological network as a hierarchy (multiple levels) based on protein binding relations, which allowed us to model the biological network at different levels of abstraction. The second part applies different levels (level 1, 2 and 3) of lumping the species together depending on the level of the hierarchy, resulting in a reduced and transformed model for each level. The model at each level is a representation of the whole system and can address questions pertinent to that level. We develop and simulate reduced models for levels-1, 2 and 3 of lumping for the G1/S checkpoint pathway and evaluate the biological plausibility of the proposed method by comparing the results with the original ODE model of Iwamoto et al. (2011). The results for continuous dynamics of the G1/S checkpoint pathway with or without DNA-damage for reduced models of level- 1, 2 and 3 of lumping are in very good agreement and consistent with the original model results and with biological findings. Therefore, the reduced models (level-1, 2 and 3) can be used to study cell cycle progression in G1 and the effects of DNA damage on it. It is suitable for reducing complex ODE biological network models while retaining accurate continuous dynamics of the system. The 3 levels of the reduced models respectively achieved 20%, 26% and 31% reduction of the base model. Moreover, the reduced model is more efficient to run (30%, 44% and 52% time reduction for the three levels) and generate solutions than the original ODE model. Simplification of complex mathematical models is possible and the proposed reduction method has the potential to make an impact across many fields of biomedical research.
Collapse
Affiliation(s)
- Mutaz Khazaaleh
- Complex Systems, Big Data and Informatics Initiative (CSBII), New Zealand
| | - Sandhya Samarasinghe
- Complex Systems, Big Data and Informatics Initiative (CSBII), New Zealand; Centre for Advanced Computational Solutions, Lincoln University, Christchurch, New Zealand.
| | - Don Kulasiri
- Complex Systems, Big Data and Informatics Initiative (CSBII), New Zealand; Centre for Advanced Computational Solutions, Lincoln University, Christchurch, New Zealand
| |
Collapse
|
2
|
Yang Q, Sing-Long CA, Reed EJ. Rapid data-driven model reduction of nonlinear dynamical systems including chemical reaction networks using ℓ 1-regularization. CHAOS (WOODBURY, N.Y.) 2020; 30:053122. [PMID: 32491878 DOI: 10.1063/1.5139463] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/20/2019] [Accepted: 04/17/2020] [Indexed: 05/21/2023]
Abstract
Large-scale nonlinear dynamical systems, such as models of atmospheric hydrodynamics, chemical reaction networks, and electronic circuits, often involve thousands or more interacting components. In order to identify key components in the complex dynamical system as well as to accelerate simulations, model reduction is often desirable. In this work, we develop a new data-driven method utilizing ℓ1-regularization for model reduction of nonlinear dynamical systems, which involves minimal parameterization and has polynomial-time complexity, allowing it to easily handle large-scale systems with as many as thousands of components in a matter of minutes. A primary objective of our model reduction method is interpretability, that is to identify key components of the dynamical system that contribute to behaviors of interest, rather than just finding an efficient projection of the dynamical system onto lower dimensions. Our method produces a family of reduced models that exhibit a trade-off between model complexity and estimation error. We find empirically that our method chooses reduced models with good extrapolation properties, an important consideration in practical applications. The reduction and extrapolation performance of our method are illustrated by applications to the Lorenz model and chemical reaction rate equations, where performance is found to be competitive with or better than state-of-the-art approaches.
Collapse
Affiliation(s)
- Q Yang
- Computer Science and Engineering Department, University of Connecticut, Storrs, Connecticut 06269, USA
| | - C A Sing-Long
- Institute for Mathematical and Computational Engineering, Pontificia Universidad Católica de Chile, Santiago 7820436, Chile
| | - E J Reed
- Department of Materials Science and Engineering, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
3
|
Gupta S, Lee REC, Faeder JR. Parallel Tempering with Lasso for model reduction in systems biology. PLoS Comput Biol 2020; 16:e1007669. [PMID: 32150537 PMCID: PMC7082068 DOI: 10.1371/journal.pcbi.1007669] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2019] [Revised: 03/19/2020] [Accepted: 01/20/2020] [Indexed: 01/08/2023] Open
Abstract
Systems Biology models reveal relationships between signaling inputs and observable molecular or cellular behaviors. The complexity of these models, however, often obscures key elements that regulate emergent properties. We use a Bayesian model reduction approach that combines Parallel Tempering with Lasso regularization to identify minimal subsets of reactions in a signaling network that are sufficient to reproduce experimentally observed data. The Bayesian approach finds distinct reduced models that fit data equivalently. A variant of this approach that uses Lasso to perform selection at the level of reaction modules is applied to the NF-κB signaling network to test the necessity of feedback loops for responses to pulsatile and continuous pathway stimulation. Taken together, our results demonstrate that Bayesian parameter estimation combined with regularization can isolate and reveal core motifs sufficient to explain data from complex signaling systems. Cells respond to diverse environmental cues using complex networks of interacting proteins and other biomolecules. Mathematical and computational models have become invaluable tools to understand these networks and make informed predictions to rationally perturb cell behavior. However, the complexity of detailed models that try to capture all known biochemical elements of signaling networks often makes it difficult to determine the key regulatory elements that are responsible for specific cell behaviors. Here, we present a Bayesian computational approach, PTLasso, to automatically extract minimal subsets of detailed models that are sufficient to explain experimental data. The method simultaneously calibrates and reduces models, and the Bayesian approach samples globally, allowing us to find alternate mechanistic explanations for the data if present. We demonstrate the method on both synthetic and real biological data and show that PTLasso is an effective method to isolate distinct parts of a larger signaling model that are sufficient for specific data.
Collapse
Affiliation(s)
- Sanjana Gupta
- Department of Computational and Systems Biology, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania, United States of America
| | - Robin E C Lee
- Department of Computational and Systems Biology, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania, United States of America
| | - James R Faeder
- Department of Computational and Systems Biology, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania, United States of America
| |
Collapse
|
4
|
|
5
|
Casagranda S, Touzeau S, Ropers D, Gouzé JL. Principal process analysis of biological models. BMC SYSTEMS BIOLOGY 2018; 12:68. [PMID: 29898718 PMCID: PMC6001159 DOI: 10.1186/s12918-018-0586-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/04/2017] [Accepted: 05/15/2018] [Indexed: 11/10/2022]
Abstract
BACKGROUND Understanding the dynamical behaviour of biological systems is challenged by their large number of components and interactions. While efforts have been made in this direction to reduce model complexity, they often prove insufficient to grasp which and when model processes play a crucial role. Answering these questions is fundamental to unravel the functioning of living organisms. RESULTS We design a method for dealing with model complexity, based on the analysis of dynamical models by means of Principal Process Analysis. We apply the method to a well-known model of circadian rhythms in mammals. The knowledge of the system trajectories allows us to decompose the system dynamics into processes that are active or inactive with respect to a certain threshold value. Process activities are graphically represented by Boolean and Dynamical Process Maps. We detect model processes that are always inactive, or inactive on some time interval. Eliminating these processes reduces the complex dynamics of the original model to the much simpler dynamics of the core processes, in a succession of sub-models that are easier to analyse. We quantify by means of global relative errors the extent to which the simplified models reproduce the main features of the original system dynamics and apply global sensitivity analysis to test the influence of model parameters on the errors. CONCLUSION The results obtained prove the robustness of the method. The analysis of the sub-model dynamics allows us to identify the source of circadian oscillations. We find that the negative feedback loop involving proteins PER, CRY, CLOCK-BMAL1 is the main oscillator, in agreement with previous modelling and experimental studies. In conclusion, Principal Process Analysis is a simple-to-use method, which constitutes an additional and useful tool for analysing the complex dynamical behaviour of biological systems.
Collapse
Affiliation(s)
- Stefano Casagranda
- Université Côte d'Azur, Inria, INRA, CNRS, UPMC Univ Paris 06, Biocore team, Sophia Antipolis, France.
| | - Suzanne Touzeau
- Université Côte d'Azur, Inria, INRA, CNRS, UPMC Univ Paris 06, Biocore team, Sophia Antipolis, France.,Université Côte d'Azur, INRA, CNRS, ISA, Sophia Antipolis, France
| | | | - Jean-Luc Gouzé
- Université Côte d'Azur, Inria, INRA, CNRS, UPMC Univ Paris 06, Biocore team, Sophia Antipolis, France
| |
Collapse
|
6
|
Amikiya EA, Banda MK. WITHDRAWN: A stoichiometric method for reducing simulation cost of chemical kinetic models. Comput Chem Eng 2018. [DOI: 10.1016/j.compchemeng.2018.05.026] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
7
|
Borovinskaya ES. New Approach for the Non‐redundant Modeling of Complex Chemical Reactions. CHEM-ING-TECH 2018. [DOI: 10.1002/cite.201700155] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
8
|
Amikiya EA, Banda MK. A stoichiometric method for reducing simulation cost of chemical kinetic models. Comput Chem Eng 2018. [DOI: 10.1016/j.compchemeng.2018.02.020] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
9
|
Yang Q, Sing-Long CA, Reed EJ. Learning reduced kinetic Monte Carlo models of complex chemistry from molecular dynamics. Chem Sci 2017; 8:5781-5796. [PMID: 28989618 PMCID: PMC5625287 DOI: 10.1039/c7sc01052d] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2017] [Accepted: 06/09/2017] [Indexed: 12/16/2022] Open
Abstract
We propose a novel statistical learning framework for automatically and efficiently building reduced kinetic Monte Carlo (KMC) models of large-scale elementary reaction networks from data generated by a single or few molecular dynamics simulations (MD). Existing approaches for identifying species and reactions from molecular dynamics typically use bond length and duration criteria, where bond duration is a fixed parameter motivated by an understanding of bond vibrational frequencies. In contrast, we show that for highly reactive systems, bond duration should be a model parameter that is chosen to maximize the predictive power of the resulting statistical model. We demonstrate our method on a high temperature, high pressure system of reacting liquid methane, and show that the learned KMC model is able to extrapolate more than an order of magnitude in time for key molecules. Additionally, our KMC model of elementary reactions enables us to isolate the most important set of reactions governing the behavior of key molecules found in the MD simulation. We develop a new data-driven algorithm to reduce the chemical reaction network which can be solved either as an integer program or efficiently using L1 regularization, and compare our results with simple count-based reduction. For our liquid methane system, we discover that rare reactions do not play a significant role in the system, and find that less than 7% of the approximately 2000 reactions observed from molecular dynamics are necessary to reproduce the molecular concentration over time of methane. The framework described in this work paves the way towards a genomic approach to studying complex chemical systems, where expensive MD simulation data can be reused to contribute to an increasingly large and accurate genome of elementary reactions and rates.
Collapse
Affiliation(s)
- Qian Yang
- Institute for Computational and Mathematical Engineering , Stanford University , Stanford , 94305 , USA .
| | - Carlos A Sing-Long
- Mathematical and Computational Engineering , School of Engineering , Pontificia Universidad Catolica de Chile , Santiago , Chile .
| | - Evan J Reed
- Institute for Computational and Mathematical Engineering , Stanford University , Stanford , 94305 , USA .
- Department of Materials Science and Engineering , Stanford University , Stanford , 94305 , USA .
| |
Collapse
|
10
|
Maiwald T, Hass H, Steiert B, Vanlier J, Engesser R, Raue A, Kipkeew F, Bock HH, Kaschek D, Kreutz C, Timmer J. Driving the Model to Its Limit: Profile Likelihood Based Model Reduction. PLoS One 2016; 11:e0162366. [PMID: 27588423 PMCID: PMC5010240 DOI: 10.1371/journal.pone.0162366] [Citation(s) in RCA: 57] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2016] [Accepted: 08/22/2016] [Indexed: 01/22/2023] Open
Abstract
In systems biology, one of the major tasks is to tailor model complexity to information content of the data. A useful model should describe the data and produce well-determined parameter estimates and predictions. Too small of a model will not be able to describe the data whereas a model which is too large tends to overfit measurement errors and does not provide precise predictions. Typically, the model is modified and tuned to fit the data, which often results in an oversized model. To restore the balance between model complexity and available measurements, either new data has to be gathered or the model has to be reduced. In this manuscript, we present a data-based method for reducing non-linear models. The profile likelihood is utilised to assess parameter identifiability and designate likely candidates for reduction. Parameter dependencies are analysed along profiles, providing context-dependent suggestions for the type of reduction. We discriminate four distinct scenarios, each associated with a specific model reduction strategy. Iterating the presented procedure eventually results in an identifiable model, which is capable of generating precise and testable predictions. Source code for all toy examples is provided within the freely available, open-source modelling environment Data2Dynamics based on MATLAB available at http://www.data2dynamics.org/, as well as the R packages dMod/cOde available at https://github.com/dkaschek/. Moreover, the concept is generally applicable and can readily be used with any software capable of calculating the profile likelihood.
Collapse
Affiliation(s)
- Tim Maiwald
- Institute of Physics, University of Freiburg, Freiburg im Breisgau, Germany
| | - Helge Hass
- Institute of Physics, University of Freiburg, Freiburg im Breisgau, Germany
| | - Bernhard Steiert
- Institute of Physics, University of Freiburg, Freiburg im Breisgau, Germany
| | - Joep Vanlier
- Institute of Physics, University of Freiburg, Freiburg im Breisgau, Germany
| | - Raphael Engesser
- Institute of Physics, University of Freiburg, Freiburg im Breisgau, Germany
| | - Andreas Raue
- Merrimack Pharmaceuticals, Boston, MA, United States of America
| | - Friederike Kipkeew
- Department of Gastroenterology, Hepatology and Infectiology, University Hospital Duesseldorf, Duesseldorf, Germany
| | - Hans H. Bock
- Department of Gastroenterology, Hepatology and Infectiology, University Hospital Duesseldorf, Duesseldorf, Germany
| | - Daniel Kaschek
- Institute of Physics, University of Freiburg, Freiburg im Breisgau, Germany
| | - Clemens Kreutz
- Institute of Physics, University of Freiburg, Freiburg im Breisgau, Germany
- Center for Biosystems Analysis (ZBSA), University of Freiburg, Freiburg im Breisgau, Germany
| | - Jens Timmer
- Institute of Physics, University of Freiburg, Freiburg im Breisgau, Germany
- Center for Biosystems Analysis (ZBSA), University of Freiburg, Freiburg im Breisgau, Germany
- BIOSS Centre for Biological Signalling Studies, University of Freiburg, Freiburg im Breisgau, Germany
- * E-mail:
| |
Collapse
|
11
|
Tsai KJ, Chang CH. Diagnostics for stochastic genome-scale modeling via model slicing and debugging. PLoS One 2014; 9:e110380. [PMID: 25368989 PMCID: PMC4219680 DOI: 10.1371/journal.pone.0110380] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2013] [Accepted: 09/21/2014] [Indexed: 12/02/2022] Open
Abstract
Modeling of biological behavior has evolved from simple gene expression plots represented by mathematical equations to genome-scale systems biology networks. However, due to obstacles in complexity and scalability of creating genome-scale models, several biological modelers have turned to programming or scripting languages and away from modeling fundamentals. In doing so, they have traded the ability to have exchangeable, standardized model representation formats, while those that remain true to standardized model representation are faced with challenges in model complexity and analysis. We have developed a model diagnostic methodology inspired by program slicing and debugging and demonstrate the effectiveness of the methodology on a genome-scale metabolic network model published in the BioModels database. The computer-aided identification revealed specific points of interest such as reversibility of reactions, initialization of species amounts, and parameter estimation that improved a candidate cell's adenosine triphosphate production. We then compared the advantages of our methodology over other modeling techniques such as model checking and model reduction. A software application that implements the methodology is available at http://gel.ym.edu.tw/gcs/.
Collapse
Affiliation(s)
- Kevin J. Tsai
- Bioinformatics Program, Taiwan International Graduate Program, Academia Sinica, Taipei, Taiwan
- Institute of Biomedical Informatics, National Yang-Ming University, Taipei, Taiwan
| | - Chuan-Hsiung Chang
- Institute of Biomedical Informatics, National Yang-Ming University, Taipei, Taiwan
- Center for Systems and Synthetic Biology, National Yang-Ming University, Taipei, Taiwan
- * E-mail:
| |
Collapse
|
12
|
Rao S, van der Schaft A, van Eunen K, Bakker BM, Jayawardhana B. A model reduction method for biochemical reaction networks. BMC SYSTEMS BIOLOGY 2014; 8:52. [PMID: 24885656 PMCID: PMC4041147 DOI: 10.1186/1752-0509-8-52] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/29/2013] [Accepted: 04/23/2014] [Indexed: 01/01/2023]
Abstract
BACKGROUND In this paper we propose a model reduction method for biochemical reaction networks governed by a variety of reversible and irreversible enzyme kinetic rate laws, including reversible Michaelis-Menten and Hill kinetics. The method proceeds by a stepwise reduction in the number of complexes, defined as the left and right-hand sides of the reactions in the network. It is based on the Kron reduction of the weighted Laplacian matrix, which describes the graph structure of the complexes and reactions in the network. It does not rely on prior knowledge of the dynamic behaviour of the network and hence can be automated, as we demonstrate. The reduced network has fewer complexes, reactions, variables and parameters as compared to the original network, and yet the behaviour of a preselected set of significant metabolites in the reduced network resembles that of the original network. Moreover the reduced network largely retains the structure and kinetics of the original model. RESULTS We apply our method to a yeast glycolysis model and a rat liver fatty acid beta-oxidation model. When the number of state variables in the yeast model is reduced from 12 to 7, the difference between metabolite concentrations in the reduced and the full model, averaged over time and species, is only 8%. Likewise, when the number of state variables in the rat-liver beta-oxidation model is reduced from 42 to 29, the difference between the reduced model and the full model is 7.5%. CONCLUSIONS The method has improved our understanding of the dynamics of the two networks. We found that, contrary to the general disposition, the first few metabolites which were deleted from the network during our stepwise reduction approach, are not those with the shortest convergence times. It shows that our reduction approach performs differently from other approaches that are based on time-scale separation. The method can be used to facilitate fitting of the parameters or to embed a detailed model of interest in a more coarse-grained yet realistic environment.
Collapse
Affiliation(s)
| | | | | | | | - Bayu Jayawardhana
- Systems Biology Center for Energy Metabolism and Ageing, University of Groningen, ERIBA, Antonius Deusinglaan 1 9713 AV Groningen, Netherlands.
| |
Collapse
|
13
|
Perumal TM, Gunawan R. pathPSA: A Dynamical Pathway-Based Parametric Sensitivity Analysis. Ind Eng Chem Res 2014. [DOI: 10.1021/ie403277d] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Thanneer Malai Perumal
- Luxembourg
Centre for Systems Biomedicine, University of Luxembourg, Esch/Alzette 4362, Luxembourg
| | - Rudiyanto Gunawan
- Institute
for Chemical and Bioengineering, ETH Zurich, Zurich 8093, Switzerland
| |
Collapse
|
14
|
Sikalo N, Hasemann O, Schulz C, Kempf A, Wlokas I. A Genetic Algorithm-Based Method for the Automatic Reduction of Reaction Mechanisms. INT J CHEM KINET 2013. [DOI: 10.1002/kin.20826] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- N. Sikalo
- Institute for Combustion and Gas Dynamics-Fluid Dynamics (IVG); University of Duisburg-Essen; 47048 Duisburg Germany
| | - O. Hasemann
- Institute for Combustion and Gas Dynamics-Fluid Dynamics (IVG); University of Duisburg-Essen; 47048 Duisburg Germany
| | - C. Schulz
- Institute for Combustion and Gas Dynamics-Reactive Fluids (IVG); University of Duisburg-Essen; 47048 Duisburg Germany
- Center for Nanointegration Duisburg-Essen (CENIDE); University of Duisburg-Essen; 47048 Duisburg Germany
| | - A. Kempf
- Institute for Combustion and Gas Dynamics-Fluid Dynamics (IVG); University of Duisburg-Essen; 47048 Duisburg Germany
- Center for Nanointegration Duisburg-Essen (CENIDE); University of Duisburg-Essen; 47048 Duisburg Germany
- Center for Computational Sciences and Simulation (CCSS); University of Duisburg-Essen; 47048 Duisburg Germany
| | - I. Wlokas
- Institute for Combustion and Gas Dynamics-Fluid Dynamics (IVG); University of Duisburg-Essen; 47048 Duisburg Germany
| |
Collapse
|
15
|
Perumal TM, Madgula Krishna S, Tallam SS, Gunawan R. Reduction of kinetic models using dynamic sensitivities. Comput Chem Eng 2013. [DOI: 10.1016/j.compchemeng.2013.05.003] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
16
|
Rodriguez-Fernandez M, Rehberg M, Kremling A, Banga JR. Simultaneous model discrimination and parameter estimation in dynamic models of cellular systems. BMC SYSTEMS BIOLOGY 2013; 7:76. [PMID: 23938131 PMCID: PMC3765209 DOI: 10.1186/1752-0509-7-76] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/08/2012] [Accepted: 08/08/2013] [Indexed: 01/06/2023]
Abstract
Background Model development is a key task in systems biology, which typically starts from an initial model candidate and, involving an iterative cycle of hypotheses-driven model modifications, leads to new experimentation and subsequent model identification steps. The final product of this cycle is a satisfactory refined model of the biological phenomena under study. During such iterative model development, researchers frequently propose a set of model candidates from which the best alternative must be selected. Here we consider this problem of model selection and formulate it as a simultaneous model selection and parameter identification problem. More precisely, we consider a general mixed-integer nonlinear programming (MINLP) formulation for model selection and identification, with emphasis on dynamic models consisting of sets of either ODEs (ordinary differential equations) or DAEs (differential algebraic equations). Results We solved the MINLP formulation for model selection and identification using an algorithm based on Scatter Search (SS). We illustrate the capabilities and efficiency of the proposed strategy with a case study considering the KdpD/KdpE system regulating potassium homeostasis in Escherichia coli. The proposed approach resulted in a final model that presents a better fit to the in silico generated experimental data. Conclusions The presented MINLP-based optimization approach for nested-model selection and identification is a powerful methodology for model development in systems biology. This strategy can be used to perform model selection and parameter estimation in one single step, thus greatly reducing the number of experiments and computations of traditional modeling approaches.
Collapse
|
17
|
Apri M, de Gee M, Molenaar J. Complexity reduction preserving dynamical behavior of biochemical networks. J Theor Biol 2012; 304:16-26. [DOI: 10.1016/j.jtbi.2012.03.019] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2011] [Revised: 02/13/2012] [Accepted: 03/15/2012] [Indexed: 12/31/2022]
|
18
|
Yuan Z, Chen B, Sin G, Gani R. State-of-the-art and progress in the optimization-based simultaneous design and control for chemical processes. AIChE J 2012. [DOI: 10.1002/aic.13786] [Citation(s) in RCA: 108] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
19
|
Lou HH, Chen D, Martin CB, Li X, Li K, Vaid H, Singh KD, Gangadharan P. Optimal Reduction of the C1–C3 Combustion Mechanism for the Simulation of Flaring. Ind Eng Chem Res 2012. [DOI: 10.1021/ie2027684] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Helen H. Lou
- Dan F.
Smith Department of Chemical Engineering, ‡Department of Chemistry and Biochemistry, and §Department of Mechanical
Engineering, Lamar University, Beaumont, Texas 77710, United States
| | - Daniel Chen
- Dan F.
Smith Department of Chemical Engineering, ‡Department of Chemistry and Biochemistry, and §Department of Mechanical
Engineering, Lamar University, Beaumont, Texas 77710, United States
| | - Christopher B. Martin
- Dan F.
Smith Department of Chemical Engineering, ‡Department of Chemistry and Biochemistry, and §Department of Mechanical
Engineering, Lamar University, Beaumont, Texas 77710, United States
| | - Xianchang Li
- Dan F.
Smith Department of Chemical Engineering, ‡Department of Chemistry and Biochemistry, and §Department of Mechanical
Engineering, Lamar University, Beaumont, Texas 77710, United States
| | - Kuyen Li
- Dan F.
Smith Department of Chemical Engineering, ‡Department of Chemistry and Biochemistry, and §Department of Mechanical
Engineering, Lamar University, Beaumont, Texas 77710, United States
| | - Hitesh Vaid
- Dan F.
Smith Department of Chemical Engineering, ‡Department of Chemistry and Biochemistry, and §Department of Mechanical
Engineering, Lamar University, Beaumont, Texas 77710, United States
| | - Kanwar Devesh Singh
- Dan F.
Smith Department of Chemical Engineering, ‡Department of Chemistry and Biochemistry, and §Department of Mechanical
Engineering, Lamar University, Beaumont, Texas 77710, United States
| | - Preeti Gangadharan
- Dan F.
Smith Department of Chemical Engineering, ‡Department of Chemistry and Biochemistry, and §Department of Mechanical
Engineering, Lamar University, Beaumont, Texas 77710, United States
| |
Collapse
|
20
|
The Rate-Controlled Constrained-Equilibrium Approach to Far-From-Local-Equilibrium Thermodynamics. ENTROPY 2012. [DOI: 10.3390/e14020092] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
21
|
Liu C, Fu J, Xu Q. Simultaneous mixed-integer dynamic optimization for environmentally benign electroplating. Comput Chem Eng 2011. [DOI: 10.1016/j.compchemeng.2011.07.004] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
22
|
(Jacky) Huang Z, Chu Y, Hahn J. Model simplification procedure for signal transduction pathway models: An application to IL-6 signaling. Chem Eng Sci 2010. [DOI: 10.1016/j.ces.2009.11.035] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
|
23
|
Nagy AM, Mourot G, Marx B, Ragot J, Schutz G. Systematic Multimodeling Methodology Applied to an Activated Sludge Reactor Model. Ind Eng Chem Res 2010. [DOI: 10.1021/ie8017687] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Anca Maria Nagy
- Centre de Recherche en Automatique de Nancy, Nancy Université, 2, Avenue de la Forêt de Haye, 54516 Vandoeuvre-lès-Nancy Cedex France, and Centre de Recherche Public “Henri Tudor”, Laboratoire de Technologies Industrielles et Matériaux, 29, Avenue John F. Kennedy, L-1855 Luxembourg-Kirchberg
| | - Gilles Mourot
- Centre de Recherche en Automatique de Nancy, Nancy Université, 2, Avenue de la Forêt de Haye, 54516 Vandoeuvre-lès-Nancy Cedex France, and Centre de Recherche Public “Henri Tudor”, Laboratoire de Technologies Industrielles et Matériaux, 29, Avenue John F. Kennedy, L-1855 Luxembourg-Kirchberg
| | - Benoît Marx
- Centre de Recherche en Automatique de Nancy, Nancy Université, 2, Avenue de la Forêt de Haye, 54516 Vandoeuvre-lès-Nancy Cedex France, and Centre de Recherche Public “Henri Tudor”, Laboratoire de Technologies Industrielles et Matériaux, 29, Avenue John F. Kennedy, L-1855 Luxembourg-Kirchberg
| | - José Ragot
- Centre de Recherche en Automatique de Nancy, Nancy Université, 2, Avenue de la Forêt de Haye, 54516 Vandoeuvre-lès-Nancy Cedex France, and Centre de Recherche Public “Henri Tudor”, Laboratoire de Technologies Industrielles et Matériaux, 29, Avenue John F. Kennedy, L-1855 Luxembourg-Kirchberg
| | - Georges Schutz
- Centre de Recherche en Automatique de Nancy, Nancy Université, 2, Avenue de la Forêt de Haye, 54516 Vandoeuvre-lès-Nancy Cedex France, and Centre de Recherche Public “Henri Tudor”, Laboratoire de Technologies Industrielles et Matériaux, 29, Avenue John F. Kennedy, L-1855 Luxembourg-Kirchberg
| |
Collapse
|
24
|
He K, Androulakis IP, Ierapetritou MG. On-the-fly reduction of kinetic mechanisms using element flux analysis. Chem Eng Sci 2010. [DOI: 10.1016/j.ces.2009.09.073] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
25
|
Al-Khateeb AN, Powers JM, Paolucci S, Sommese AJ, Diller JA, Hauenstein JD, Mengers JD. One-dimensional slow invariant manifolds for spatially homogenous reactive systems. J Chem Phys 2009; 131:024118. [DOI: 10.1063/1.3171613] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
Hsu SH, Stamatis SD, Caruthers JM, Delgass WN, Venkatasubramanian V, Blau GE, Lasinski M, Orcun S. Bayesian Framework for Building Kinetic Models of Catalytic Systems. Ind Eng Chem Res 2009. [DOI: 10.1021/ie801651y] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
27
|
Maurya MR, Bornheimer SJ, Venkatasubramanian V, Subramaniam S. Mixed-integer nonlinear optimisation approach to coarse-graining biochemical networks. IET Syst Biol 2009; 3:24-39. [PMID: 19154082 DOI: 10.1049/iet-syb:20080098] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Quantitative modelling and analysis of biochemical networks is challenging because of the inherent complexities and nonlinearities of the system and the limited availability of parameter values. Even if a mathematical model of the network can be developed, the lack of large-scale good-quality data makes accurate estimation of a large number of parameters impossible. Hence, coarse-grained models (CGMs) consisting of essential biochemical mechanisms are more suitable for computational analysis and for studying important systemic functions. The central question in constructing a CGM is which mechanisms should be deemed 'essential' and which can be ignored. Also, how should parameter values be defined when data are sparse? A mixed-integer nonlinear-programming (MINLP) based optimisation approach to coarse-graining is presented. Starting with a detailed biochemical model with associated computational details (reaction network and mathematical description) and data on the biochemical system, the structure and the parameters of a CGM can be determined simultaneously. In this optimisation problem, the authors use a genetic algorithm to simultaneously identify parameter values and remove unimportant reactions. The methodology is exemplified by developing two CGMs for the GTPase-cycle module of M1 muscarinic acetylcholine receptor, Gq, and regulator of G protein signalling 4 [RGS4, a GTPase-activating protein (GAP)] starting from a detailed model of 48 reactions. Both the CGMs have only 17 reactions, fit experimental data well and predict, as does the detailed model, four limiting signalling regimes (LSRs) corresponding to the extremes of receptor and GAP concentration. The authors demonstrate that coarse-graining, in addition to resulting in a reduced-order model, also provides insights into the mechanisms in the network. The best CGM obtained for the GTPase cycle also contains an unconventional mechanism and its predictions explain an old problem in pharmacology, the biphasic (bell-shaped) response to certain drugs. The MINLP methodology is broadly applicable to larger and complex (dense) biochemical modules.
Collapse
Affiliation(s)
- M R Maurya
- University of California, San Diego, Department of Bioengineering, La Jolla, CA 92093, USA
| | | | | | | |
Collapse
|
28
|
He K, Ierapetritou MG, Androulakis IP. Integration of on-the-fly kinetic reduction with multidimensional CFD. AIChE J 2009. [DOI: 10.1002/aic.12072] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
29
|
Oscillator model reduction preserving the phase response: application to the circadian clock. Biophys J 2008; 95:1658-73. [PMID: 18487303 DOI: 10.1529/biophysj.107.128678] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Mathematical model reduction is a long-standing technique used both to gain insight into model subprocesses and to reduce the computational costs of simulation and analysis. A reduced model must retain essential features of the full model, which, traditionally, have been the trajectories of certain state variables. For biological clocks, timing, or phase, characteristics must be preserved. A key performance criterion for a clock is the ability to adjust its phase correctly in response to external signals. We present a novel model reduction technique that removes components from a single-oscillator clock model and discover that four feedback loops are redundant with respect to its phase response behavior. Using a coupled multioscillator model of a circadian clock, we demonstrate that by preserving the phase response behavior of a single oscillator, we preserve timing behavior at the multioscillator level.
Collapse
|
30
|
Reinhardt V, Winckler M, Lebiedz D. Approximation of slow attracting manifolds in chemical kinetics by trajectory-based optimization approaches. J Phys Chem A 2008; 112:1712-8. [PMID: 18247506 DOI: 10.1021/jp0739925] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Many common kinetic model reduction approaches are explicitly based on inherent multiple time scales and often assume and directly exploit a clear time scale separation into fast and slow reaction processes. They approximate the system dynamics with a dimension-reduced model after eliminating the fast modes by enslaving them to the slow ones. The corresponding restrictive assumption of full relaxation of fast modes often renders the resulting approximation of slow attracting manifolds inaccurate as a representation of the reduced model and makes the numerical solution of the nonlinear "reduction equations" particularly difficult in many cases where the gap in intrinsic time scales is not large enough. We demonstrate that trajectory optimization approaches can avoid such severe restrictions by computing numerical solutions that correspond to "maximally relaxed" dynamical modes in a suitable sense. We present a framework of trajectory-based optimization for model reduction in chemical kinetics and a general class of reduction criteria characterizing the relaxation of chemical forces along reaction trajectories. These criteria can be motivated geometrically exploiting ideas from differential geometry and fundamental physics and turn out to be highly successful in example applications. Within this framework, we provide results for the computational approximation of slow attracting low-dimensional manifolds in terms of families of optimal trajectories for a six-component hydrogen combustion mechanism.
Collapse
Affiliation(s)
- V Reinhardt
- Interdisciplinary Center for Scientific Computing, Im Neuenheimer Feld 368, D-69120 Heidelberg, Germany
| | | | | |
Collapse
|
31
|
Dong GQ, Jakobowski L, Iafolla MAJ, McMillen DR. Simplification of stochastic chemical reaction models with fast and slow dynamics. J Biol Phys 2007; 33:67-95. [PMID: 19669554 DOI: 10.1007/s10867-007-9043-2] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2006] [Accepted: 05/29/2007] [Indexed: 11/28/2022] Open
Abstract
Biological systems often involve chemical reactions occurring in low-molecule-number regimes, where fluctuations are not negligible and thus stochastic models are required to capture the system behaviour. The resulting models are generally quite large and complex, involving many reactions and species. For clarity and computational tractability, it is important to be able to simplify these systems to equivalent ones involving fewer elements. While many model simplification approaches have been developed for deterministic systems, there has been limited work on applying these approaches to stochastic modelling. Here, we describe a method that reduces the complexity of stochastic biochemical network models, and apply this method to the reduction of a mammalian signalling cascade and a detailed model of the process of bacterial gene expression. Our results indicate that the simplified model gives an accurate representation for not only the average numbers of all species, but also for the associated fluctuations and statistical parameters.
Collapse
Affiliation(s)
- Guang Qiang Dong
- Institute for Optical Sciences and Department of Chemical and Physical Sciences, University of Toronto Mississauga, Mississauga, ON L5L 1C6, Canada
| | | | | | | |
Collapse
|
32
|
Tonse SR, Day MS, Brown NJ. Dynamic reduction of a CH4/air chemical mechanism appropriate for investigating vortex–flame interactions. INT J CHEM KINET 2007. [DOI: 10.1002/kin.20227] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
33
|
Frenklach M, Packard A, Feeley R. Chapter 6 Optimization of Reaction Models with Solution Mapping. MODELING OF CHEMICAL REACTIONS 2007. [DOI: 10.1016/s0069-8040(07)42006-4] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
|
34
|
Maurya M, Bornheimer S, Venkatasubramanian V, Subramaniam S. Reduced-order modelling of biochemical networks: application to the GTPase-cycle signalling module. ACTA ACUST UNITED AC 2006; 152:229-42. [PMID: 16986265 PMCID: PMC3417759 DOI: 10.1049/ip-syb:20050014] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
Abstract
Biochemical systems embed complex networks and hence development and analysis of their detailed models pose a challenge for computation. Coarse-grained biochemical models, called reduced-order models (ROMs), consisting of essential biochemical mechanisms are more useful for computational analysis and for studying important features of a biochemical network. The authors present a novel method to model-reduction by identifying potentially important parameters using multidimensional sensitivity analysis. A ROM is generated for the GTPase-cycle module of m1 muscarinic acetylcholine receptor, Gq, and regulator of G-protein signalling 4 (a GTPase-activating protein or GAP) starting from a detailed model of 48 reactions. The resulting ROM has only 17 reactions. The ROM suggested that complexes of G-protein coupled receptor (GPCR) and GAP--which were proposed in the detailed model as a hypothesis--are required to fit the experimental data. Models previously published in the literature are also simulated and compared with the ROM. Through this comparison, a minimal ROM, that also requires complexes of GPCR and GAP, with just 15 parameters is generated. The proposed reduced-order modelling methodology is scalable to larger networks and provides a general framework for the reduction of models of biochemical systems.
Collapse
Affiliation(s)
- M.R. Maurya
- San Diego Supercomputer Center, 9500 Gilman Drive MC 0505, La Jolla, CA 92093, USA
| | - S.J. Bornheimer
- Departments of Chemistry and Biochemistry and Cellular and Molecular Medicine, University of California, San Diego, 9500 Gilman Drive La Jolla, CA 92093, USA
| | - V. Venkatasubramanian
- Laboratory for Intelligent Process Systems, School of Chemical Engineering, Purdue University, West Lafayette, IN 47907, USA
| | - S. Subramaniam
- San Diego Supercomputer Center, 9500 Gilman Drive MC 0505, La Jolla, CA 92093, USA, the Departments of Chemistry and Biochemistry and Cellular and Molecular Medicine, University of California, San Diego, 9500 Gilman Drive La Jolla, CA 92093, USA and the Department of Bioengineering, University of California, San Diego, 9500 Gilman Drive La Jolla, CA 92093, USA
| |
Collapse
|
35
|
Chachuat B, Singer AB, Barton PI. Global Methods for Dynamic Optimization and Mixed-Integer Dynamic Optimization. Ind Eng Chem Res 2006. [DOI: 10.1021/ie0601605] [Citation(s) in RCA: 82] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Benoît Chachuat
- Automatic Control Laboratory, Swiss Federal Institute of Technology, Station 9, CH-1015 Lausanne, Switzerland, ExxonMobil Upstream Research Company, Houston, Texas, and Process Systems Engineering Laboratory, Department of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Adam B. Singer
- Automatic Control Laboratory, Swiss Federal Institute of Technology, Station 9, CH-1015 Lausanne, Switzerland, ExxonMobil Upstream Research Company, Houston, Texas, and Process Systems Engineering Laboratory, Department of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Paul I. Barton
- Automatic Control Laboratory, Swiss Federal Institute of Technology, Station 9, CH-1015 Lausanne, Switzerland, ExxonMobil Upstream Research Company, Houston, Texas, and Process Systems Engineering Laboratory, Department of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| |
Collapse
|
36
|
Gokulakrishnan P, Lawrence A, McLellan P, Grandmaison E. A functional-PCA approach for analyzing and reducing complex chemical mechanisms. Comput Chem Eng 2006. [DOI: 10.1016/j.compchemeng.2006.02.007] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
37
|
A systematic framework for the design of reduced-order models for signal transduction pathways from a control theoretic perspective. Comput Chem Eng 2006. [DOI: 10.1016/j.compchemeng.2005.10.003] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
38
|
Sirdeshpande AR, Ierapetritou MG, Androulakis IP. Design of flexible reduced kinetic mechanisms. AIChE J 2006. [DOI: 10.1002/aic.690471110] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
39
|
|
40
|
Shaik OS, Kammerer J, Gorecki J, Lebiedz D. Derivation of a quantitative minimal model from a detailed elementary-step mechanism supported by mathematical coupling analysis. J Chem Phys 2005; 123:234103. [PMID: 16392910 DOI: 10.1063/1.2136882] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accurate experimental data increasingly allow the development of detailed elementary-step mechanisms for complex chemical and biochemical reaction systems. Model reduction techniques are widely applied to obtain representations in lower-dimensional phase space which are more suitable for mathematical analysis, efficient numerical simulation, and model-based control tasks. Here, we exploit a recently implemented numerical algorithm for error-controlled computation of the minimum dimension required for a still accurate reduced mechanism based on automatic time scale decomposition and relaxation of fast modes. We determine species contributions to the active (slow) dynamical modes of the reaction system and exploit this information in combination with quasi-steady-state and partial-equilibrium approximations for explicit model reduction of a novel detailed chemical mechanism for the Ru-catalyzed light-sensitive Belousov-Zhabotinsky reaction. The existence of a minimum dimension of seven is demonstrated to be mandatory for the reduced model to show good quantitative consistency with the full model in numerical simulations. We derive such a maximally reduced seven-variable model from the detailed elementary-step mechanism and demonstrate that it reproduces quantitatively accurately the dynamical features of the full model within a given accuracy tolerance.
Collapse
Affiliation(s)
- O S Shaik
- Interdisciplinary Center for Scientific Computing, Im Neuenheimer Feld 368, D-69120 Heidelberg, Germany
| | | | | | | |
Collapse
|
41
|
Lebiedz D, Kammerer J, Brandt-Pollmann U. Automatic network coupling analysis for dynamical systems based on detailed kinetic models. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2005; 72:041911. [PMID: 16383424 DOI: 10.1103/physreve.72.041911] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/09/2005] [Indexed: 05/05/2023]
Abstract
We introduce a numerical complexity reduction method for the automatic identification and analysis of dynamic network decompositions in (bio)chemical kinetics based on error-controlled computation of a minimal model dimension represented by the number of (locally) active dynamical modes. Our algorithm exploits a generalized sensitivity analysis along state trajectories and subsequent singular value decomposition of sensitivity matrices for the identification of these dominant dynamical modes. It allows for a dynamic coupling analysis of (bio)chemical species in kinetic models that can be exploited for the piecewise computation of a minimal model on small time intervals and offers valuable functional insight into highly nonlinear reaction mechanisms and network dynamics. We present results for the identification of network decompositions in a simple oscillatory chemical reaction, time scale separation based model reduction in a Michaelis-Menten enzyme system and network decomposition of a detailed model for the oscillatory peroxidase-oxidase enzyme system.
Collapse
Affiliation(s)
- Dirk Lebiedz
- Interdisciplinary Center for Scientific Computing, Im Neuenheimer Feld 368, D-69120 Heidelberg, Germany.
| | | | | |
Collapse
|
42
|
|
43
|
|
44
|
Androulakis IP. “Store and retrieve” representations of dynamic systems motivated by studies in gas phase chemical kinetics. Comput Chem Eng 2004. [DOI: 10.1016/j.compchemeng.2004.02.038] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
45
|
Lebiedz D. Computing minimal entropy production trajectories: An approach to model reduction in chemical kinetics. J Chem Phys 2004; 120:6890-7. [PMID: 15267587 DOI: 10.1063/1.1652428] [Citation(s) in RCA: 55] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Advanced experimental techniques in chemistry and physics provide increasing access to detailed deterministic mass action models for chemical reaction kinetics. Especially in complex technical or biochemical systems the huge amount of species and reaction pathways involved in a detailed modeling approach call for efficient methods of model reduction. These should be automatic and based on a firm mathematical analysis of the ordinary differential equations underlying the chemical kinetics in deterministic models. A main purpose of model reduction is to enable accurate numerical simulations of even high dimensional and spatially extended reaction systems. The latter include physical transport mechanisms and are modeled by partial differential equations. Their numerical solution for hundreds or thousands of species within a reasonable time will exceed computer capacities available now and in a foreseeable future. The central idea of model reduction is to replace the high dimensional dynamics by a low dimensional approximation with an appropriate degree of accuracy. Here I present a global approach to model reduction based on the concept of minimal entropy production and its numerical implementation. For given values of a single species concentration in a chemical system all other species concentrations are computed under the assumption that the system is as close as possible to its attractor, the thermodynamic equilibrium, in the sense that all modes of thermodynamic forces are maximally relaxed except the one, which drives the remaining system dynamics. This relaxation is expressed in terms of minimal entropy production for single reaction steps along phase space trajectories.
Collapse
Affiliation(s)
- D Lebiedz
- Interdisciplinary Center for Scientific Computing, Im Neuenheimer Feld 368, D-69120 Heidelberg, Germany.
| |
Collapse
|
46
|
Katare S, Caruthers JM, Delgass WN, Venkatasubramanian V. An Intelligent System for Reaction Kinetic Modeling and Catalyst Design. Ind Eng Chem Res 2004. [DOI: 10.1021/ie034067h] [Citation(s) in RCA: 58] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Santhoji Katare
- Center for Catalyst Design, School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907
| | - James M. Caruthers
- Center for Catalyst Design, School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907
| | - W. Nicholas Delgass
- Center for Catalyst Design, School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907
| | - Venkat Venkatasubramanian
- Center for Catalyst Design, School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907
| |
Collapse
|
47
|
Tonse SR, Brown NJ. Dimensionality estimate of the manifold in chemical composition space for a turbulent premixed H2 + air flame. INT J CHEM KINET 2004. [DOI: 10.1002/kin.20002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
48
|
|
49
|
Tonse SR, Moriarty NW, Frenklach M, Brown NJ. Computational economy improvements in PRISM. INT J CHEM KINET 2003. [DOI: 10.1002/kin.10140] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
50
|
Caruthers J, Lauterbach J, Thomson K, Venkatasubramanian V, Snively C, Bhan A, Katare S, Oskarsdottir G. Catalyst design: knowledge extraction from high-throughput experimentation. J Catal 2003. [DOI: 10.1016/s0021-9517(02)00036-2] [Citation(s) in RCA: 63] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|