1
|
Lu YA, McCann MG, Hu WS, Zhang Q. Multi-cell-line learning for the data-driven construction of mechanistic metabolic models. Biotechnol Bioeng 2024; 121:2833-2847. [PMID: 38831695 DOI: 10.1002/bit.28757] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2023] [Revised: 04/25/2024] [Accepted: 05/19/2024] [Indexed: 06/05/2024]
Abstract
Mammalian cells are commonly used as hosts in cell culture for biologics production in the pharmaceutical industry. Structured mechanistic models of metabolism have been used to capture complex cellular mechanisms that contribute to varying metabolic shifts in different cell lines. However, little research has focused on the impact of temporal changes in enzyme abundance and activity on the modeling of cell metabolism. In this work, we present a framework for constructing mechanistic models of metabolism that integrate growth-signaling control of enzyme activity and transcript dynamics. The proposed approach is applied to build models for three Chinese hamster ovary (CHO) cell lines using fed-batch culture data and time-series transcript profiles. Leveraging information from the transcriptome data, we develop a parameter estimation approach based on multi-cell-line (MCL) learning, which combines data sets from different cell lines and trains the individual cell-line models jointly to improve model accuracy. The computational results demonstrate the important role of growth signaling and transcript variability in metabolic models as well as the virtue of the MCL approach for constructing cell-line models with a limited amount of data. The resulting models exhibit a high level of accuracy in predicting distinct metabolic behaviors in the different cell lines; these models can potentially be used to accelerate the process and cell-line development for the biomanufacturing of new protein therapeutics.
Collapse
Affiliation(s)
- Yen-An Lu
- Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota, USA
| | - Meghan G McCann
- Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota, USA
| | - Wei-Shou Hu
- Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota, USA
| | - Qi Zhang
- Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota, USA
| |
Collapse
|
2
|
Lange E, Kranert L, Krüger J, Benndorf D, Heyer R. Microbiome modeling: a beginner's guide. Front Microbiol 2024; 15:1368377. [PMID: 38962127 PMCID: PMC11220171 DOI: 10.3389/fmicb.2024.1368377] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2024] [Accepted: 05/27/2024] [Indexed: 07/05/2024] Open
Abstract
Microbiomes, comprised of diverse microbial species and viruses, play pivotal roles in human health, environmental processes, and biotechnological applications and interact with each other, their environment, and hosts via ecological interactions. Our understanding of microbiomes is still limited and hampered by their complexity. A concept improving this understanding is systems biology, which focuses on the holistic description of biological systems utilizing experimental and computational methods. An important set of such experimental methods are metaomics methods which analyze microbiomes and output lists of molecular features. These lists of data are integrated, interpreted, and compiled into computational microbiome models, to predict, optimize, and control microbiome behavior. There exists a gap in understanding between microbiologists and modelers/bioinformaticians, stemming from a lack of interdisciplinary knowledge. This knowledge gap hinders the establishment of computational models in microbiome analysis. This review aims to bridge this gap and is tailored for microbiologists, researchers new to microbiome modeling, and bioinformaticians. To achieve this goal, it provides an interdisciplinary overview of microbiome modeling, starting with fundamental knowledge of microbiomes, metaomics methods, common modeling formalisms, and how models facilitate microbiome control. It concludes with guidelines and repositories for modeling. Each section provides entry-level information, example applications, and important references, serving as a valuable resource for comprehending and navigating the complex landscape of microbiome research and modeling.
Collapse
Affiliation(s)
- Emanuel Lange
- Multidimensional Omics Data Analysis, Department for Bioanalytics, Leibniz-Institut für Analytische Wissenschaften - ISAS - e.V., Dortmund, Germany
- Graduate School Digital Infrastructure for the Life Sciences, Bielefeld Institute for Bioinformatics Infrastructure (BIBI), Faculty of Technology, Bielefeld University, Bielefeld, Germany
| | - Lena Kranert
- Institute for Automation Engineering, Otto von Guericke University Magdeburg, Magdeburg, Germany
| | - Jacob Krüger
- Engineering of Software-Intensive Systems, Department of Mathematics and Computer Science, Eindhoven University of Technology, Eindhoven, Netherlands
| | - Dirk Benndorf
- Applied Biosciences and Bioprocess Engineering, Anhalt University of Applied Sciences, Köthen, Germany
| | - Robert Heyer
- Multidimensional Omics Data Analysis, Department for Bioanalytics, Leibniz-Institut für Analytische Wissenschaften - ISAS - e.V., Dortmund, Germany
- Graduate School Digital Infrastructure for the Life Sciences, Bielefeld Institute for Bioinformatics Infrastructure (BIBI), Faculty of Technology, Bielefeld University, Bielefeld, Germany
- Multidimensional Omics Data Analysis, Faculty of Technology, Bielefeld University, Bielefeld, Germany
| |
Collapse
|
3
|
Zhang Y, Toyoda F, Himeno Y, Noma A, Amano A. Cell-specific models of hiPSC-CMs developed by the gradient-based parameter optimization method fitting two different action potential waveforms. Sci Rep 2024; 14:13086. [PMID: 38849433 PMCID: PMC11161598 DOI: 10.1038/s41598-024-63413-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2023] [Accepted: 05/28/2024] [Indexed: 06/09/2024] Open
Abstract
Parameter optimization (PO) methods to determine the ionic current composition of experimental cardiac action potential (AP) waveform have been developed using a computer model of cardiac membrane excitation. However, it was suggested that fitting a single AP record in the PO method was not always successful in providing a unique answer because of a shortage of information. We found that the PO method worked perfectly if the PO method was applied to a pair of a control AP and a model output AP in which a single ionic current out of six current species, such as IKr, ICaL, INa, IKs, IKur or IbNSC was partially blocked in silico. When the target was replaced by a pair of experimental control and IKr-blocked records of APs generated spontaneously in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), the simultaneous fitting of the two waveforms by the PO method was hampered to some extent by the irregular slow fluctuations in the Vm recording and/or sporadic alteration in AP configurations in the hiPSC-CMs. This technical problem was largely removed by selecting stable segments of the records for the PO method. Moreover, the PO method was made fail-proof by running iteratively in identifying the optimized parameter set to reconstruct both the control and the IKr-blocked AP waveforms. In the lead potential analysis, the quantitative ionic mechanisms deduced from the optimized parameter set were totally consistent with the qualitative view of ionic mechanisms of AP so far described in physiological literature.
Collapse
Affiliation(s)
- Yixin Zhang
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Futoshi Toyoda
- Department of Physiology, Shiga University of Medical Science, Otsu, Japan
- Central Research Laboratory, Shiga University of Medical Science, Otsu, Japan
| | - Yukiko Himeno
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan.
- Department of Physiology, Shiga University of Medical Science, Otsu, Japan.
| | - Akinori Noma
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Akira Amano
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| |
Collapse
|
4
|
Iida K, Okada M. Identifying Key Regulatory Genes in Drug Resistance Acquisition: Modeling Pseudotime Trajectories of Breast Cancer Single-Cell Transcriptome. Cancers (Basel) 2024; 16:1884. [PMID: 38791962 PMCID: PMC11119661 DOI: 10.3390/cancers16101884] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2024] [Revised: 05/11/2024] [Accepted: 05/15/2024] [Indexed: 05/26/2024] Open
Abstract
Single-cell RNA-sequencing (scRNA-seq) technology has provided significant insights into cancer drug resistance at the single-cell level. However, understanding dynamic cell transitions at the molecular systems level remains limited, requiring a systems biology approach. We present an approach that combines mathematical modeling with a pseudotime analysis using time-series scRNA-seq data obtained from the breast cancer cell line MCF-7 treated with tamoxifen. Our single-cell analysis identified five distinct subpopulations, including tamoxifen-sensitive and -resistant groups. Using a single-gene mathematical model, we discovered approximately 560-680 genes out of 6000 exhibiting multistable expression states in each subpopulation, including key estrogen-receptor-positive breast cancer cell survival genes, such as RPS6KB1. A bifurcation analysis elucidated their regulatory mechanisms, and we mapped these genes into a molecular network associated with cell survival and metastasis-related pathways. Our modeling approach comprehensively identifies key regulatory genes for drug resistance acquisition, enhancing our understanding of potential drug targets in breast cancer.
Collapse
Affiliation(s)
- Keita Iida
- Institute for Protein Research, Osaka University, Suita 565-0871, Osaka, Japan;
| | | |
Collapse
|
5
|
Lang PF, Penas DR, Banga JR, Weindl D, Novak B. Reusable rule-based cell cycle model explains compartment-resolved dynamics of 16 observables in RPE-1 cells. PLoS Comput Biol 2024; 20:e1011151. [PMID: 38190398 PMCID: PMC10773963 DOI: 10.1371/journal.pcbi.1011151] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2023] [Accepted: 11/24/2023] [Indexed: 01/10/2024] Open
Abstract
The mammalian cell cycle is regulated by a well-studied but complex biochemical reaction system. Computational models provide a particularly systematic and systemic description of the mechanisms governing mammalian cell cycle control. By combining both state-of-the-art multiplexed experimental methods and powerful computational tools, this work aims at improving on these models along four dimensions: model structure, validation data, validation methodology and model reusability. We developed a comprehensive model structure of the full cell cycle that qualitatively explains the behaviour of human retinal pigment epithelial-1 cells. To estimate the model parameters, time courses of eight cell cycle regulators in two compartments were reconstructed from single cell snapshot measurements. After optimisation with a parallel global optimisation metaheuristic we obtained excellent agreements between simulations and measurements. The PEtab specification of the optimisation problem facilitates reuse of model, data and/or optimisation results. Future perturbation experiments will improve parameter identifiability and allow for testing model predictive power. Such a predictive model may aid in drug discovery for cell cycle-related disorders.
Collapse
Affiliation(s)
- Paul F. Lang
- Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| | - David R. Penas
- Computational Biology Lab, MBG-CSIC (Spanish National Research Council), Pontevedra, Spain
| | - Julio R. Banga
- Computational Biology Lab, MBG-CSIC (Spanish National Research Council), Pontevedra, Spain
| | - Daniel Weindl
- Computational Health Center, Helmholtz Zentrum München Deutsches Forschungszentrum für Gesundheit und Umwelt (GmbH), Neuherberg, Germany
| | - Bela Novak
- Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
6
|
Raimúndez E, Fedders M, Hasenauer J. Posterior marginalization accelerates Bayesian inference for dynamical models of biological processes. iScience 2023; 26:108083. [PMID: 37867942 PMCID: PMC10589897 DOI: 10.1016/j.isci.2023.108083] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2023] [Revised: 07/16/2023] [Accepted: 09/25/2023] [Indexed: 10/24/2023] Open
Abstract
Bayesian inference is an important method in the life and natural sciences for learning from data. It provides information about parameter and prediction uncertainties. Yet, generating representative samples from the posterior distribution is often computationally challenging. Here, we present an approach that lowers the computational complexity of sample generation for dynamical models with scaling, offset, and noise parameters. The proposed method is based on the marginalization of the posterior distribution. We provide analytical results for a broad class of problems with conjugate priors and show that the method is suitable for a large number of applications. Subsequently, we demonstrate the benefit of the approach for applications from the field of systems biology. We report an improvement up to 50 times in the effective sample size per unit of time. As the scheme is broadly applicable, it will facilitate Bayesian inference in different research fields.
Collapse
Affiliation(s)
- Elba Raimúndez
- Life and Medical Sciences (LIMES) Institute, University of Bonn, Bonn, Germany
- Technische Universität München, Center for Mathematics, Garching, Germany
| | - Michael Fedders
- Life and Medical Sciences (LIMES) Institute, University of Bonn, Bonn, Germany
| | - Jan Hasenauer
- Life and Medical Sciences (LIMES) Institute, University of Bonn, Bonn, Germany
- Technische Universität München, Center for Mathematics, Garching, Germany
- Helmholtz Zentrum München - German Research Center for Environmental Health, Computational Health Center, Neuherberg, Germany
| |
Collapse
|
7
|
Egert J, Kreutz C. Realistic simulation of time-course measurements in systems biology. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2023; 20:10570-10589. [PMID: 37322949 DOI: 10.3934/mbe.2023467] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
In systems biology, the analysis of complex nonlinear systems faces many methodological challenges. For the evaluation and comparison of the performances of novel and competing computational methods, one major bottleneck is the availability of realistic test problems. We present an approach for performing realistic simulation studies for analyses of time course data as they are typically measured in systems biology. Since the design of experiments in practice depends on the process of interest, our approach considers the size and the dynamics of the mathematical model which is intended to be used for the simulation study. To this end, we used 19 published systems biology models with experimental data and evaluated the relationship between model features (e.g., the size and the dynamics) and features of the measurements such as the number and type of observed quantities, the number and the selection of measurement times, and the magnitude of measurement errors. Based on these typical relationships, our novel approach enables suggestions of realistic simulation study designs in the systems biology context and the realistic generation of simulated data for any dynamic model. The approach is demonstrated on three models in detail and its performance is validated on nine models by comparing ODE integration, parameter optimization, and parameter identifiability. The presented approach enables more realistic and less biased benchmark studies and thereby constitutes an important tool for the development of novel methods for dynamic modeling.
Collapse
Affiliation(s)
- Janine Egert
- Institute of Medical Biometry and Statistics (IMBI), Faculty of Medicine and Medical Center - University of Freiburg, Stefan-Meier-Str. 26, 79104 Freiburg, Germany
- Centre for Integrative Biological Signalling Studies CIBSS, University of Freiburg, 79104 Freiburg, Germany
| | - Clemens Kreutz
- Institute of Medical Biometry and Statistics (IMBI), Faculty of Medicine and Medical Center - University of Freiburg, Stefan-Meier-Str. 26, 79104 Freiburg, Germany
- Centre for Integrative Biological Signalling Studies CIBSS, University of Freiburg, 79104 Freiburg, Germany
- Center for Data Analysis and Modelling (FDM), University of Freiburg, 79104 Freiburg, Germany
| |
Collapse
|
8
|
Geertjens NHJ, de Vink PJ, Wezeman T, Markvoort AJ, Brunsveld L. Straightforward model construction and analysis of multicomponent biomolecular systems in equilibrium. RSC Chem Biol 2023; 4:252-260. [PMID: 37034407 PMCID: PMC10074561 DOI: 10.1039/d2cb00211f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2022] [Accepted: 01/17/2023] [Indexed: 01/19/2023] Open
Abstract
Mathematical modelling of molecular systems can be extremely helpful in elucidating complex phenomena in (bio)chemistry. However, equilibrium conditions in systems consisting of more than two components, such as for molecular glues bound to two proteins, can typically not be analytically determined without assumptions and (semi-)numerical models are not trivial to derive by the non-expert. Here we present a framework for equilibrium models, geared towards molecular glues and other contemporary multicomponent chemical biology challenges. The framework utilizes a general derivation method capable of generating custom mass-balance models for equilibrium conditions of complex molecular systems, based on the simple, reversible biomolecular reactions describing these systems. Several chemical biology concepts are revisited via the framework to demonstrate the simplicity, generality and validity of the approach. The ease of use of the framework and the ability to both analyze systems and gain additional insights in the underlying parameters driving equilibria formation strongly aids the analysis and understanding of biomolecular systems. New directions for research and analysis are brought forward based on the model formation and system and parameter analysis. This conceptual framework severely reduces the time and expertise requirements which currently impede the broad integration of such valuable equilibrium models into molecular glue development and chemical biology research.
Collapse
Affiliation(s)
- Nick H J Geertjens
- Laboratory of Chemical Biology, Department of Biomedical Engineering and Institute for Complex Molecular Systems, Eindhoven University of Technology, Den Dolech 2 Eindhoven 5612 AZ The Netherlands
| | - Pim J de Vink
- Laboratory of Chemical Biology, Department of Biomedical Engineering and Institute for Complex Molecular Systems, Eindhoven University of Technology, Den Dolech 2 Eindhoven 5612 AZ The Netherlands
| | - Tim Wezeman
- Laboratory of Chemical Biology, Department of Biomedical Engineering and Institute for Complex Molecular Systems, Eindhoven University of Technology, Den Dolech 2 Eindhoven 5612 AZ The Netherlands
| | - Albert J Markvoort
- Computational Biology Group, Department of Biomedical Engineering and Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513 Eindhoven 5600 MB The Netherlands
| | - Luc Brunsveld
- Laboratory of Chemical Biology, Department of Biomedical Engineering and Institute for Complex Molecular Systems, Eindhoven University of Technology, Den Dolech 2 Eindhoven 5612 AZ The Netherlands
| |
Collapse
|
9
|
Lao-Martil D, Schmitz JPJ, Teusink B, van Riel NAW. Elucidating yeast glycolytic dynamics at steady state growth and glucose pulses through kinetic metabolic modeling. Metab Eng 2023; 77:128-142. [PMID: 36963461 DOI: 10.1016/j.ymben.2023.03.005] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2022] [Revised: 02/12/2023] [Accepted: 03/05/2023] [Indexed: 03/26/2023]
Abstract
Microbial cell factories face changing environments during industrial fermentations. Kinetic metabolic models enable the simulation of the dynamic metabolic response to these perturbations, but their development is challenging due to model complexity and experimental data requirements. An example of this is the well-established microbial cell factory Saccharomyces cerevisiae, for which no consensus kinetic model of central metabolism has been developed and implemented in industry. Here, we aim to bring the academic and industrial communities closer to this consensus model. We developed a physiology informed kinetic model of yeast glycolysis connected to central carbon metabolism by including the effect of anabolic reactions precursors, mitochondria and the trehalose cycle. To parametrize such a large model, a parameter estimation pipeline was developed, consisting of a divide and conquer approach, supplemented with regularization and global optimization. Additionally, we show how this first mechanistic description of a growing yeast cell captures experimental dynamics at different growth rates and under a strong glucose perturbation, is robust to parametric uncertainty and explains the contribution of the different pathways in the network. Such a comprehensive model could not have been developed without using steady state and glucose perturbation data sets. The resulting metabolic reconstruction and parameter estimation pipeline can be applied in the future to study other industrially-relevant scenarios. We show this by generating a hybrid CFD-metabolic model to explore intracellular glycolytic dynamics for the first time. The model suggests that all intracellular metabolites oscillate within a physiological range, except carbon storage metabolism, which is sensitive to the extracellular environment.
Collapse
Affiliation(s)
- David Lao-Martil
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, Noord-Brabant, 5612AE, the Netherlands
| | - Joep P J Schmitz
- DSM Biotechnology Center, Delft, Zuid-Holland, 2613AX, the Netherlands
| | - Bas Teusink
- Systems Biology Lab, Vrije Universiteit Amsterdam, Amsterdam, Noord-Holland, 1081HZ, the Netherlands
| | - Natal A W van Riel
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, Noord-Brabant, 5612AE, the Netherlands; Amsterdam University Medical Center, University of Amsterdam, Amsterdam, Noord-Holland, 1105AZ, the Netherlands.
| |
Collapse
|
10
|
Martino CF, Jimenez P, Goldfarb M, Abdulla UG. Optimization of parameters in coherent spin dynamics of radical pairs in quantum biology. PLoS One 2023; 18:e0273404. [PMID: 36827400 PMCID: PMC9956872 DOI: 10.1371/journal.pone.0273404] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2021] [Accepted: 08/08/2022] [Indexed: 02/26/2023] Open
Abstract
Identification of the external electromagnetic fields and internal hyperfine parameters which optimize the quantum singlet-triplet yield of simplified radical pairs modeled by Schrödinger system with spin Hamiltonians given by the sum of Zeeman interaction and hyperfine coupling interaction terms are analyzed. A method that combines sensitivity analysis with Tikhonov regularization is implemented. Numerical results demonstrate that the quantum singlet-triplet yield of the radical pair system can be significantly reduced if optimization is pursued simultaneously for both external magnetic fields and internal hyperfine parameters. The results may contribute towards understanding the structure-function relationship of a putative magnetoreceptor to manipulate and enhance quantum coherences at room temperature and leveraging biofidelic function to inspire novel quantum devices.
Collapse
Affiliation(s)
- Carlos F. Martino
- Johns Hopkins University Applied Physics Lab, Laurel, MD, United States of America
- * E-mail:
| | - Pablo Jimenez
- Departamento de Física Médica, Instituto Balseiro, San Carlos de Bariloche, Argentina
| | - Max Goldfarb
- Department of Mathematical Sciences, Florida Institute of Technology, Melbourne, Fl, United States of America
| | - Ugur G. Abdulla
- Department of Mathematical Sciences, Florida Institute of Technology, Melbourne, Fl, United States of America
- Analysis & PDE Unit, Okinawa Institute of Science and Technology, Okinawa, Japan
| |
Collapse
|
11
|
Wieczorek PZ, Starecki K, Gołofit K, Radtke M, Pilarz M. A Tiny Flexible Differential Tension Sensor. SENSORS (BASEL, SWITZERLAND) 2023; 23:1819. [PMID: 36850420 PMCID: PMC9966848 DOI: 10.3390/s23041819] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/09/2023] [Revised: 02/03/2023] [Accepted: 02/03/2023] [Indexed: 06/18/2023]
Abstract
Modern applications of Internet of Things (IoT) devices require cheap and effective methods of measurement of physical quantities. Cheap IoT devices with sensor functionalities can detect a lack or excess of substances in everyday life or industry processes. One possible use of tension sensors in IoT applications is the automated replenishment process of fast moving consumer goods (FMCG) on shop shelves or home retail automation that allows for quick ordering of FMCG, where the IoT system is a part of smart packaging. For those reasons, a growing demand for cheap and tiny tension sensors has arisen. In this article, we propose a solution of a small flexible tension sensor fabricated in an amorphous InGaZnO (a-IGZO) thin-film process that can be integrated with other devices, e.g., near-field communications (NFC) or a barcode radio frequency identification (RFID) tag. The sensor was designed to magnify the slight internal changes in material properties caused by mechanical stress. These changes affect the dynamic electrical properties of specially designed inverters for a pair of ring oscillators, in which the frequencies become stress-dependent. In the article, we discuss and explain the approach to the optimum design of a ring oscillator that manifests the highest sensitivity to mechanical stress.
Collapse
Affiliation(s)
- Piotr Z. Wieczorek
- Faculty of Electronics and Information Technology, Institute of Electronic Systems, Warsaw University of Technology, 00-665 Warsaw, Poland
| | - Krzysztof Starecki
- Faculty of Electronics and Information Technology, Institute of Electronic Systems, Warsaw University of Technology, 00-665 Warsaw, Poland
| | - Krzysztof Gołofit
- Faculty of Electronics and Information Technology, Institute of Electronic Systems, Warsaw University of Technology, 00-665 Warsaw, Poland
| | - Maciej Radtke
- Faculty of Electronics and Information Technology, Institute of Electronic Systems, Warsaw University of Technology, 00-665 Warsaw, Poland
| | | |
Collapse
|
12
|
Lambert B, Lei CL, Robinson M, Clerx M, Creswell R, Ghosh S, Tavener S, Gavaghan DJ. Autocorrelated measurement processes and inference for ordinary differential equation models of biological systems. J R Soc Interface 2023. [PMCID: PMC9943878 DOI: 10.1098/rsif.2022.0725] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2023] Open
Abstract
Ordinary differential equation models are used to describe dynamic processes across biology. To perform likelihood-based parameter inference on these models, it is necessary to specify a statistical process representing the contribution of factors not explicitly included in the mathematical model. For this, independent Gaussian noise is commonly chosen, with its use so widespread that researchers typically provide no explicit justification for this choice. This noise model assumes ‘random’ latent factors affect the system in the ephemeral fashion resulting in unsystematic deviation of observables from their modelled counterparts. However, like the deterministically modelled parts of a system, these latent factors can have persistent effects on observables. Here, we use experimental data from dynamical systems drawn from cardiac physiology and electrochemistry to demonstrate that highly persistent differences between observations and modelled quantities can occur. Considering the case when persistent noise arises owing only to measurement imperfections, we use the Fisher information matrix to quantify how uncertainty in parameter estimates is artificially reduced when erroneously assuming independent noise. We present a workflow to diagnose persistent noise from model fits and describe how to remodel accounting for correlated errors.
Collapse
Affiliation(s)
- Ben Lambert
- Department of Mathematics, University of Exeter, Exeter EX4 4PY, UK
| | - Chon Lok Lei
- Faculty of Health Sciences, Institute of Translational Medicine, University of Macau, Macau, People’s Republic of China
| | - Martin Robinson
- Department of Computer Science, University of Oxford, Oxford OX1 3QG, UK
| | - Michael Clerx
- School of Mathematical Sciences, University of Nottingham, Nottingham NG7 2RD, UK
| | - Richard Creswell
- Department of Computer Science, University of Oxford, Oxford OX1 3QG, UK
| | - Sanmitra Ghosh
- MRC Biostatistics Unit, University of Cambridge, Cambridge CB2 0SR, UK
| | - Simon Tavener
- Department of Mathematics, Colorado State University, Fort Collins, CO 80523, UK
| | - David J. Gavaghan
- Department of Computer Science, University of Oxford, Oxford OX1 3QG, UK
| |
Collapse
|
13
|
Fröhlich F, Gerosa L, Muhlich J, Sorger PK. Mechanistic model of MAPK signaling reveals how allostery and rewiring contribute to drug resistance. Mol Syst Biol 2023; 19:e10988. [PMID: 36700386 PMCID: PMC9912026 DOI: 10.15252/msb.202210988] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2022] [Revised: 11/29/2022] [Accepted: 12/15/2022] [Indexed: 01/27/2023] Open
Abstract
BRAF is prototypical of oncogenes that can be targeted therapeutically and the treatment of BRAFV600E melanomas with RAF and MEK inhibitors results in rapid tumor regression. However, drug-induced rewiring generates a drug adapted state thought to be involved in acquired resistance and disease recurrence. In this article, we study mechanisms of adaptive rewiring in BRAFV600E melanoma cells using an energy-based implementation of ordinary differential equation (ODE) modeling in combination with proteomic, transcriptomic and imaging data. We develop a method for causal tracing of ODE models and identify two parallel MAPK reaction channels that are differentially sensitive to RAF and MEK inhibitors due to differences in protein oligomerization and drug binding. We describe how these channels, and timescale separation between immediate-early signaling and transcriptional feedback, create a state in which the RAS-regulated MAPK channel can be activated by growth factors under conditions in which the BRAFV600E -driven channel is fully inhibited. Further development of the approaches in this article is expected to yield a unified model of adaptive drug resistance in melanoma.
Collapse
Affiliation(s)
- Fabian Fröhlich
- Laboratory of Systems Pharmacology, Department of Systems BiologyHarvard Medical SchoolBostonMAUSA
| | - Luca Gerosa
- Laboratory of Systems Pharmacology, Department of Systems BiologyHarvard Medical SchoolBostonMAUSA,Present address:
Genentech, Inc.South San FranciscoCAUSA
| | - Jeremy Muhlich
- Laboratory of Systems Pharmacology, Department of Systems BiologyHarvard Medical SchoolBostonMAUSA
| | - Peter K Sorger
- Laboratory of Systems Pharmacology, Department of Systems BiologyHarvard Medical SchoolBostonMAUSA
| |
Collapse
|
14
|
Using Kinetic Modelling to Infer Adaptations in Saccharomyces cerevisiae Carbohydrate Storage Metabolism to Dynamic Substrate Conditions. Metabolites 2023; 13:metabo13010088. [PMID: 36677014 PMCID: PMC9862193 DOI: 10.3390/metabo13010088] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 12/14/2022] [Accepted: 12/23/2022] [Indexed: 01/07/2023] Open
Abstract
Microbial metabolism is strongly dependent on the environmental conditions. While these can be well controlled under laboratory conditions, large-scale bioreactors are characterized by inhomogeneities and consequently dynamic conditions for the organisms. How Saccharomyces cerevisiae response to frequent perturbations in industrial bioreactors is still not understood mechanistically. To study the adjustments to prolonged dynamic conditions, we used published repeated substrate perturbation regime experimental data, extended it with proteomic measurements and used both for modelling approaches. Multiple types of data were combined; including quantitative metabolome, 13C enrichment and flux quantification data. Kinetic metabolic modelling was applied to study the relevant intracellular metabolic response dynamics. An existing model of yeast central carbon metabolism was extended, and different subsets of enzymatic kinetic constants were estimated. A novel parameter estimation pipeline based on combinatorial enzyme selection supplemented by regularization was developed to identify and predict the minimum enzyme and parameter adjustments from steady-state to dynamic substrate conditions. This approach predicted proteomic changes in hexose transport and phosphorylation reactions, which were additionally confirmed by proteome measurements. Nevertheless, the modelling also hints at a yet unknown kinetic or regulation phenomenon. Some intracellular fluxes could not be reproduced by mechanistic rate laws, including hexose transport and intracellular trehalase activity during substrate perturbation cycles.
Collapse
|
15
|
Lakrisenko P, Stapor P, Grein S, Paszkowski Ł, Pathirana D, Fröhlich F, Lines GT, Weindl D, Hasenauer J. Efficient computation of adjoint sensitivities at steady-state in ODE models of biochemical reaction networks. PLoS Comput Biol 2023; 19:e1010783. [PMID: 36595539 PMCID: PMC9838866 DOI: 10.1371/journal.pcbi.1010783] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2022] [Revised: 01/13/2023] [Accepted: 12/01/2022] [Indexed: 01/04/2023] Open
Abstract
Dynamical models in the form of systems of ordinary differential equations have become a standard tool in systems biology. Many parameters of such models are usually unknown and have to be inferred from experimental data. Gradient-based optimization has proven to be effective for parameter estimation. However, computing gradients becomes increasingly costly for larger models, which are required for capturing the complex interactions of multiple biochemical pathways. Adjoint sensitivity analysis has been pivotal for working with such large models, but methods tailored for steady-state data are currently not available. We propose a new adjoint method for computing gradients, which is applicable if the experimental data include steady-state measurements. The method is based on a reformulation of the backward integration problem to a system of linear algebraic equations. The evaluation of the proposed method using real-world problems shows a speedup of total simulation time by a factor of up to 4.4. Our results demonstrate that the proposed approach can achieve a substantial improvement in computation time, in particular for large-scale models, where computational efficiency is critical.
Collapse
Affiliation(s)
- Polina Lakrisenko
- Computational Health Center, Helmholtz Zentrum München Deutsches Forschungszentrum für Gesundheit und Umwelt (GmbH), Neuherberg, Germany
- Center for Mathematics, Technische Universität München, Garching, Germany
| | - Paul Stapor
- Computational Health Center, Helmholtz Zentrum München Deutsches Forschungszentrum für Gesundheit und Umwelt (GmbH), Neuherberg, Germany
- Center for Mathematics, Technische Universität München, Garching, Germany
| | - Stephan Grein
- University of Bonn, Life and Medical Sciences Institute, Bonn, Germany
| | | | - Dilan Pathirana
- University of Bonn, Life and Medical Sciences Institute, Bonn, Germany
| | - Fabian Fröhlich
- Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America
| | | | - Daniel Weindl
- Computational Health Center, Helmholtz Zentrum München Deutsches Forschungszentrum für Gesundheit und Umwelt (GmbH), Neuherberg, Germany
| | - Jan Hasenauer
- Computational Health Center, Helmholtz Zentrum München Deutsches Forschungszentrum für Gesundheit und Umwelt (GmbH), Neuherberg, Germany
- University of Bonn, Life and Medical Sciences Institute, Bonn, Germany
| |
Collapse
|
16
|
Estimation of Small-Scale Kinetic Parameters of Escherichia coli (E. coli) Model by Enhanced Segment Particle Swarm Optimization Algorithm ESe-PSO. Processes (Basel) 2023. [DOI: 10.3390/pr11010126] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
The ability to create “structured models” of biological simulations is becoming more and more commonplace. Although computer simulations can be used to estimate the model, they are restricted by the lack of experimentally available parameter values, which must be approximated. In this study, an Enhanced Segment Particle Swarm Optimization (ESe-PSO) algorithm that can estimate the values of small-scale kinetic parameters is described and applied to E. coli’s main metabolic network as a model system. The glycolysis, phosphotransferase system, pentose phosphate, the TCA cycle, gluconeogenesis, glyoxylate pathways, and acetate formation pathways of Escherichia coli are represented by the Differential Algebraic Equations (DAE) system for the metabolic network. However, this algorithm uses segments to organize particle movements and the dynamic inertia weight (ω) to increase the algorithm’s exploration and exploitation potential. As an alternative to the state-of-the-art algorithm, this adjustment improves estimation accuracy. The numerical findings indicate a good agreement between the observed and predicted data. In this regard, the result of the ESe-PSO algorithm achieved superior accuracy compared with the Segment Particle Swarm Optimization (Se-PSO), Particle Swarm Optimization (PSO), Genetic Algorithm (GA), and Differential Evolution (DE) algorithms. As a result of this innovative approach, it was concluded that small-scale and even entire cell kinetic model parameters can be developed.
Collapse
|
17
|
Fröhlich F. A Practical Guide for the Efficient Formulation and Calibration of Large, Energy- and Rule-Based Models of Cellular Signal Transduction. Methods Mol Biol 2023; 2634:59-86. [PMID: 37074574 DOI: 10.1007/978-1-0716-3008-2_3] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/20/2023]
Abstract
Aberrant signal transduction leads to complex diseases such as cancer. To rationally design treatment strategies with small molecule inhibitors, computational models have to be employed. Energy- and rule-based models allow the construction of mechanistic ordinary differential equation models based on structural insights. The detailed, energy-based description often generates large models, which are difficult to calibrate on experimental data. In this chapter, we provide a detailed, interactive protocol for the programmatic formulation and calibration of such large, energy- and rule-based models of cellular signal transduction based on an example model describing the action of RAF inhibitors on MAPK signaling. An interactive version of this chapter is available as Jupyter Notebook at github.com/FFroehlich/energy_modeling_chapter .
Collapse
Affiliation(s)
- Fabian Fröhlich
- Department of Systems Biology, Harvard Medical School, Boston, MA, USA.
| |
Collapse
|
18
|
Vittadello ST, Stumpf MPH. Open problems in mathematical biology. Math Biosci 2022; 354:108926. [PMID: 36377100 DOI: 10.1016/j.mbs.2022.108926] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2022] [Revised: 10/21/2022] [Accepted: 10/21/2022] [Indexed: 11/06/2022]
Abstract
Biology is data-rich, and it is equally rich in concepts and hypotheses. Part of trying to understand biological processes and systems is therefore to confront our ideas and hypotheses with data using statistical methods to determine the extent to which our hypotheses agree with reality. But doing so in a systematic way is becoming increasingly challenging as our hypotheses become more detailed, and our data becomes more complex. Mathematical methods are therefore gaining in importance across the life- and biomedical sciences. Mathematical models allow us to test our understanding, make testable predictions about future behaviour, and gain insights into how we can control the behaviour of biological systems. It has been argued that mathematical methods can be of great benefit to biologists to make sense of data. But mathematics and mathematicians are set to benefit equally from considering the often bewildering complexity inherent to living systems. Here we present a small selection of open problems and challenges in mathematical biology. We have chosen these open problems because they are of both biological and mathematical interest.
Collapse
Affiliation(s)
- Sean T Vittadello
- Melbourne Integrative Genomics, University of Melbourne, Australia; School of BioSciences, University of Melbourne, Australia
| | - Michael P H Stumpf
- Melbourne Integrative Genomics, University of Melbourne, Australia; School of BioSciences, University of Melbourne, Australia; School of Mathematics and Statistics, University of Melbourne, Australia.
| |
Collapse
|
19
|
Duong-Trung N, Born S, Kim JW, Schermeyer MT, Paulick K, Borisyak M, Cruz-Bournazou MN, Werner T, Scholz R, Schmidt-Thieme L, Neubauer P, Martinez E. When Bioprocess Engineering Meets Machine Learning: A Survey from the Perspective of Automated Bioprocess Development. Biochem Eng J 2022. [DOI: 10.1016/j.bej.2022.108764] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
|
20
|
Massonis G, Villaverde AF, Banga JR. Improving dynamic predictions with ensembles of observable models. Bioinformatics 2022; 39:6842325. [PMID: 36416122 PMCID: PMC9805594 DOI: 10.1093/bioinformatics/btac755] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Revised: 10/20/2022] [Accepted: 11/22/2022] [Indexed: 11/24/2022] Open
Abstract
MOTIVATION Dynamic mechanistic modelling in systems biology has been hampered by the complexity and variability associated with the underlying interactions, and by uncertain and sparse experimental measurements. Ensemble modelling, a concept initially developed in statistical mechanics, has been introduced in biological applications with the aim of mitigating those issues. Ensemble modelling uses a collection of different models compatible with the observed data to describe the phenomena of interest. However, since systems biology models often suffer from a lack of identifiability and observability, ensembles of models are particularly unreliable when predicting non-observable states. RESULTS We present a strategy to assess and improve the reliability of a class of model ensembles. In particular, we consider kinetic models described using ordinary differential equations with a fixed structure. Our approach builds an ensemble with a selection of the parameter vectors found when performing parameter estimation with a global optimization metaheuristic. This technique enforces diversity during the sampling of parameter space and it can quantify the uncertainty in the predictions of state trajectories. We couple this strategy with structural identifiability and observability analysis, and when these tests detect possible prediction issues we obtain model reparameterizations that surmount them. The end result is an ensemble of models with the ability to predict the internal dynamics of a biological process. We demonstrate our approach with models of glucose regulation, cell division, circadian oscillations and the JAK-STAT signalling pathway. AVAILABILITY AND IMPLEMENTATION The code that implements the methodology and reproduces the results is available at https://doi.org/10.5281/zenodo.6782638. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Gemma Massonis
- Computational Biology Lab, MBG-CSIC (Spanish National Research Council), Pontevedra, Galicia 36143, Spain
| | | | | |
Collapse
|
21
|
Kochen MA, Wiley HS, Feng S, Sauro HM. SBbadger: biochemical reaction networks with definable degree distributions. Bioinformatics 2022; 38:5064-5072. [PMID: 36111865 PMCID: PMC9665861 DOI: 10.1093/bioinformatics/btac630] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2022] [Revised: 07/24/2022] [Accepted: 09/15/2022] [Indexed: 12/24/2022] Open
Abstract
MOTIVATION An essential step in developing computational tools for the inference, optimization and simulation of biochemical reaction networks is gauging tool performance against earlier efforts using an appropriate set of benchmarks. General strategies for the assembly of benchmark models include collection from the literature, creation via subnetwork extraction and de novo generation. However, with respect to biochemical reaction networks, these approaches and their associated tools are either poorly suited to generate models that reflect the wide range of properties found in natural biochemical networks or to do so in numbers that enable rigorous statistical analysis. RESULTS In this work, we present SBbadger, a python-based software tool for the generation of synthetic biochemical reaction or metabolic networks with user-defined degree distributions, multiple available kinetic formalisms and a host of other definable properties. SBbadger thus enables the creation of benchmark model sets that reflect properties of biological systems and generate the kinetics and model structures typically targeted by computational analysis and inference software. Here, we detail the computational and algorithmic workflow of SBbadger, demonstrate its performance under various settings, provide sample outputs and compare it to currently available biochemical reaction network generation software. AVAILABILITY AND IMPLEMENTATION SBbadger is implemented in Python and is freely available at https://github.com/sys-bio/SBbadger and via PyPI at https://pypi.org/project/SBbadger/. Documentation can be found at https://SBbadger.readthedocs.io. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Michael A Kochen
- Department of Bioengineering, University of Washington, Seattle, WA 98105, USA
| | - H Steven Wiley
- Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, WA, 99354, USA
| | - Song Feng
- Biological Science Division, Pacific Northwest National Laboratory, Richland, WA 99354, USA
| | - Herbert M Sauro
- Department of Bioengineering, University of Washington, Seattle, WA 98105, USA
| |
Collapse
|
22
|
Kohjitani H, Koda S, Himeno Y, Makiyama T, Yamamoto Y, Yoshinaga D, Wuriyanghai Y, Kashiwa A, Toyoda F, Zhang Y, Amano A, Noma A, Kimura T. Gradient-based parameter optimization method to determine membrane ionic current composition in human induced pluripotent stem cell-derived cardiomyocytes. Sci Rep 2022; 12:19110. [PMID: 36351955 PMCID: PMC9646722 DOI: 10.1038/s41598-022-23398-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2022] [Accepted: 10/31/2022] [Indexed: 11/11/2022] Open
Abstract
Premature cardiac myocytes derived from human induced pluripotent stem cells (hiPSC-CMs) show heterogeneous action potentials (APs), probably due to different expression patterns of membrane ionic currents. We developed a method for determining expression patterns of functional channels in terms of whole-cell ionic conductance (Gx) using individual spontaneous AP configurations. It has been suggested that apparently identical AP configurations can be obtained using different sets of ionic currents in mathematical models of cardiac membrane excitation. If so, the inverse problem of Gx estimation might not be solved. We computationally tested the feasibility of the gradient-based optimization method. For a realistic examination, conventional 'cell-specific models' were prepared by superimposing the model output of AP on each experimental AP recorded by conventional manual adjustment of Gxs of the baseline model. Gxs of 4-6 major ionic currents of the 'cell-specific models' were randomized within a range of ± 5-15% and used as an initial parameter set for the gradient-based automatic Gxs recovery by decreasing the mean square error (MSE) between the target and model output. Plotting all data points of the MSE-Gx relationship during optimization revealed progressive convergence of the randomized population of Gxs to the original value of the cell-specific model with decreasing MSE. The absence of any other local minimum in the global search space was confirmed by mapping the MSE by randomizing Gxs over a range of 0.1-10 times the control. No additional local minimum MSE was obvious in the whole parameter space, in addition to the global minimum of MSE at the default model parameter.
Collapse
Affiliation(s)
- Hirohiko Kohjitani
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Shigeya Koda
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Yukiko Himeno
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Takeru Makiyama
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Yuta Yamamoto
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Daisuke Yoshinaga
- grid.258799.80000 0004 0372 2033Department Pediatrics, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Yimin Wuriyanghai
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Asami Kashiwa
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Futoshi Toyoda
- grid.410827.80000 0000 9747 6806Department of Physiology, Shiga University of Medical Science, Otsu, Japan
| | - Yixin Zhang
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Akira Amano
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Akinori Noma
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Takeshi Kimura
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| |
Collapse
|
23
|
Evaluation of the Influence of Rootstock Type on the Yield Parameters of Vines Using a Mathematical Model in Nontraditional Wine-Growing Conditions. APPLIED SCIENCES-BASEL 2022. [DOI: 10.3390/app12147293] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Great interest in viticulture in temperate climates results from the introduction of new interspecies hybrids of grapevines which are quite popular due to their high resistance to fungal diseases and lower temperature. However, the impact of rootstocks, common in vine cultivation, is little to not known, which makes setting up vineyards a challenge. This study aimed to evaluate the effect of the following six rootstock types: 101-14 Mgt, SORI, 161-49 C, 5 BB, SO4, 125 AA, and grapevines with their own roots on the yield quantity and berry quality (expressed by Brix extract) of Regent grapevines in temperate climates (southeastern Poland). A five-year experiment alongside a novel numerical model is applied to formulate precise and constructive findings about the rootstock impact in a temperate climate. Both the experimental and numerical part are supported by detailed statistical analysis. The five-year period of study indicates that the vines on rootstock 125 AA yielded the best, significantly. Shrubs grafted on rootstock 161-49 yielded the lowest, while the fruit extract content grafted on rootstock 101-14 was significantly lower among the evaluated ones. The parameters of own-rooted bushes and those grafted on SO4 rootstock did not differ significantly, except for the extract. The model leads to convergent conclusions with statistical analysis of raw experimental data. The 125 AA rootstock was the best for all nine tested case scenarios. On the other hand, 161-49 rootstock was the weakest, justified only in the most challenging conditions.
Collapse
|
24
|
Eriksson O, Bhalla US, Blackwell KT, Crook SM, Keller D, Kramer A, Linne ML, Saudargienė A, Wade RC, Hellgren Kotaleski J. Combining hypothesis- and data-driven neuroscience modeling in FAIR workflows. eLife 2022; 11:e69013. [PMID: 35792600 PMCID: PMC9259018 DOI: 10.7554/elife.69013] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2021] [Accepted: 05/13/2022] [Indexed: 12/22/2022] Open
Abstract
Modeling in neuroscience occurs at the intersection of different points of view and approaches. Typically, hypothesis-driven modeling brings a question into focus so that a model is constructed to investigate a specific hypothesis about how the system works or why certain phenomena are observed. Data-driven modeling, on the other hand, follows a more unbiased approach, with model construction informed by the computationally intensive use of data. At the same time, researchers employ models at different biological scales and at different levels of abstraction. Combining these models while validating them against experimental data increases understanding of the multiscale brain. However, a lack of interoperability, transparency, and reusability of both models and the workflows used to construct them creates barriers for the integration of models representing different biological scales and built using different modeling philosophies. We argue that the same imperatives that drive resources and policy for data - such as the FAIR (Findable, Accessible, Interoperable, Reusable) principles - also support the integration of different modeling approaches. The FAIR principles require that data be shared in formats that are Findable, Accessible, Interoperable, and Reusable. Applying these principles to models and modeling workflows, as well as the data used to constrain and validate them, would allow researchers to find, reuse, question, validate, and extend published models, regardless of whether they are implemented phenomenologically or mechanistically, as a few equations or as a multiscale, hierarchical system. To illustrate these ideas, we use a classical synaptic plasticity model, the Bienenstock-Cooper-Munro rule, as an example due to its long history, different levels of abstraction, and implementation at many scales.
Collapse
Affiliation(s)
- Olivia Eriksson
- Science for Life Laboratory, School of Electrical Engineering and Computer Science, KTH Royal Institute of TechnologyStockholmSweden
| | - Upinder Singh Bhalla
- National Center for Biological Sciences, Tata Institute of Fundamental ResearchBangaloreIndia
| | - Kim T Blackwell
- Department of Bioengineering, Volgenau School of Engineering, George Mason UniversityFairfaxUnited States
| | - Sharon M Crook
- School of Mathematical and Statistical Sciences, Arizona State UniversityTempeUnited States
| | - Daniel Keller
- Blue Brain Project, École Polytechnique Fédérale de LausanneLausanneSwitzerland
| | - Andrei Kramer
- Science for Life Laboratory, School of Electrical Engineering and Computer Science, KTH Royal Institute of TechnologyStockholmSweden
- Department of Neuroscience, Karolinska InstituteStockholmSweden
| | - Marja-Leena Linne
- Faculty of Medicine and Health Technology, Tampere UniversityTampereFinland
| | - Ausra Saudargienė
- Neuroscience Institute, Lithuanian University of Health SciencesKaunasLithuania
- Department of Informatics, Vytautas Magnus UniversityKaunasLithuania
| | - Rebecca C Wade
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS)HeidelbergGermany
- Center for Molecular Biology (ZMBH), ZMBH-DKFZ Alliance, University of HeidelbergHeidelbergGermany
- Interdisciplinary Center for Scientific Computing (IWR), Heidelberg UniversityHeidelbergGermany
| | - Jeanette Hellgren Kotaleski
- Science for Life Laboratory, School of Electrical Engineering and Computer Science, KTH Royal Institute of TechnologyStockholmSweden
- Department of Neuroscience, Karolinska InstituteStockholmSweden
| |
Collapse
|
25
|
Fröhlich F, Sorger PK. Fides: Reliable trust-region optimization for parameter estimation of ordinary differential equation models. PLoS Comput Biol 2022; 18:e1010322. [PMID: 35830470 PMCID: PMC9312381 DOI: 10.1371/journal.pcbi.1010322] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2022] [Revised: 07/25/2022] [Accepted: 06/21/2022] [Indexed: 11/18/2022] Open
Abstract
Ordinary differential equation (ODE) models are widely used to study biochemical reactions in cellular networks since they effectively describe the temporal evolution of these networks using mass action kinetics. The parameters of these models are rarely known a priori and must instead be estimated by calibration using experimental data. Optimization-based calibration of ODE models on is often challenging, even for low-dimensional problems. Multiple hypotheses have been advanced to explain why biochemical model calibration is challenging, including non-identifiability of model parameters, but there are few comprehensive studies that test these hypotheses, likely because tools for performing such studies are also lacking. Nonetheless, reliable model calibration is essential for uncertainty analysis, model comparison, and biological interpretation. We implemented an established trust-region method as a modular Python framework (fides) to enable systematic comparison of different approaches to ODE model calibration involving a variety of Hessian approximation schemes. We evaluated fides on a recently developed corpus of biologically realistic benchmark problems for which real experimental data are available. Unexpectedly, we observed high variability in optimizer performance among different implementations of the same mathematical instructions (algorithms). Analysis of possible sources of poor optimizer performance identified limitations in the widely used Gauss-Newton, BFGS and SR1 Hessian approximation schemes. We addressed these drawbacks with a novel hybrid Hessian approximation scheme that enhances optimizer performance and outperforms existing hybrid approaches. When applied to the corpus of test models, we found that fides was on average more reliable and efficient than existing methods using a variety of criteria. We expect fides to be broadly useful for ODE constrained optimization problems in biochemical models and to be a foundation for future methods development.
Collapse
Affiliation(s)
- Fabian Fröhlich
- Laboratory of Systems Pharmacology and Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Peter K. Sorger
- Laboratory of Systems Pharmacology and Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America
| |
Collapse
|
26
|
Liang E, Zhang S, Liu B, Qi B, Nie Y, Yuan Z. Mathematical Modeling for the Industrial 2-Mercaptobenzothiazole Batch Production Process. ACS OMEGA 2022; 7:6963-6977. [PMID: 35252688 PMCID: PMC8892915 DOI: 10.1021/acsomega.1c06646] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/23/2021] [Accepted: 02/08/2022] [Indexed: 06/14/2023]
Abstract
As an important chemical intermediate, 2-mercaptobenzothiazole (MBT) is widely used in various processes, especially in the rubber industry. However, there is no first-principles model that describes the synthetic process of MBT. This paper focuses on the formulation of a reliable mathematical model represented by a series of differential and algebraic equations for the industrial batch MBT production process. It is difficult to estimate all of the unknown parameters in the model because of the lack of sufficient industrial/experimental data. Thus, a reduced estimable parameter set is derived by performing estimability analysis on the original estimation problem. A multiple-starting-point strategy is then applied to numerically solve the non-convex parameter estimation problem with the weighted least-squares approach. Subsequently, a cross-validation technique is employed to evaluate the generalizability of the proposed model. Finally, it is confirmed that the proposed model produces encouraging prediction performance with regard to independent test data.
Collapse
Affiliation(s)
- Enzhi Liang
- State
Key Laboratory of Chemical Engineering, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Song Zhang
- China
Sunsine Chemical Holdings Ltd., Shanxian City, Shandong 274300, China
| | - Bin Liu
- China
Sunsine Chemical Holdings Ltd., Shanxian City, Shandong 274300, China
| | - Bujin Qi
- China
Sunsine Chemical Holdings Ltd., Shanxian City, Shandong 274300, China
| | - Yanpei Nie
- China
Sunsine Chemical Holdings Ltd., Shanxian City, Shandong 274300, China
| | - Zhihong Yuan
- State
Key Laboratory of Chemical Engineering, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| |
Collapse
|
27
|
Lai X, Keller C, Santos G, Schaft N, Dörrie J, Vera J. Multi-Level Computational Modeling of Anti-Cancer Dendritic Cell Vaccination Utilized to Select Molecular Targets for Therapy Optimization. Front Cell Dev Biol 2022; 9:746359. [PMID: 35186943 PMCID: PMC8847669 DOI: 10.3389/fcell.2021.746359] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2021] [Accepted: 12/23/2021] [Indexed: 01/18/2023] Open
Abstract
Dendritic cells (DCs) can be used for therapeutic vaccination against cancer. The success of this therapy depends on efficient tumor-antigen presentation to cytotoxic T lymphocytes (CTLs) and the induction of durable CTL responses by the DCs. Therefore, simulation of such a biological system by computational modeling is appealing because it can improve our understanding of the molecular mechanisms underlying CTL induction by DCs and help identify new strategies to improve therapeutic DC vaccination for cancer. Here, we developed a multi-level model accounting for the life cycle of DCs during anti-cancer immunotherapy. Specifically, the model is composed of three parts representing different stages of DC immunotherapy - the spreading and bio-distribution of intravenously injected DCs in human organs, the biochemical reactions regulating the DCs' maturation and activation, and DC-mediated activation of CTLs. We calibrated the model using quantitative experimental data that account for the activation of key molecular circuits within DCs, the bio-distribution of DCs in the body, and the interaction between DCs and T cells. We showed how such a data-driven model can be exploited in combination with sensitivity analysis and model simulations to identify targets for enhancing anti-cancer DC vaccination. Since other previous works show how modeling improves therapy schedules and DC dosage, we here focused on the molecular optimization of the therapy. In line with this, we simulated the effect in DC vaccination of the concerted modulation of combined intracellular regulatory processes and proposed several possibilities that can enhance DC-mediated immunogenicity. Taken together, we present a comprehensive time-resolved multi-level model for studying DC vaccination in melanoma. Although the model is not intended for personalized patient therapy, it could be used as a tool for identifying molecular targets for optimizing DC-based therapy for cancer, which ultimately should be tested in in vitro and in vivo experiments.
Collapse
Affiliation(s)
- Xin Lai
- Laboratory of Systems Tumor Immunology, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
- Deutsches Zentrum Immuntherapie and Comprehensive Cancer Center Erlangen-EMN, Erlangen, Germany
| | - Christine Keller
- Laboratory of Systems Tumor Immunology, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
| | - Guido Santos
- Laboratory of Systems Tumor Immunology, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
- Departament of Biochemistry, Microbiology, Cell Biology and Genetics, Faculty of Sciences, University of La Laguna, San Cristóbal de La Laguna, Spain
| | - Niels Schaft
- Deutsches Zentrum Immuntherapie and Comprehensive Cancer Center Erlangen-EMN, Erlangen, Germany
- RNA Group, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
| | - Jan Dörrie
- Deutsches Zentrum Immuntherapie and Comprehensive Cancer Center Erlangen-EMN, Erlangen, Germany
- RNA Group, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
| | - Julio Vera
- Laboratory of Systems Tumor Immunology, Department of Dermatology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) and Universitätsklinikum Erlangen, Erlangen, Germany
- Deutsches Zentrum Immuntherapie and Comprehensive Cancer Center Erlangen-EMN, Erlangen, Germany
| |
Collapse
|
28
|
Regularization and concave loss functions for estimation of chemical kinetic models. Appl Soft Comput 2022. [DOI: 10.1016/j.asoc.2021.108286] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
29
|
Villaverde AF, Pathirana D, Fröhlich F, Hasenauer J, Banga JR. A protocol for dynamic model calibration. Brief Bioinform 2022; 23:bbab387. [PMID: 34619769 PMCID: PMC8769694 DOI: 10.1093/bib/bbab387] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2021] [Revised: 08/06/2021] [Accepted: 08/29/2021] [Indexed: 12/23/2022] Open
Abstract
Ordinary differential equation models are nowadays widely used for the mechanistic description of biological processes and their temporal evolution. These models typically have many unknown and nonmeasurable parameters, which have to be determined by fitting the model to experimental data. In order to perform this task, known as parameter estimation or model calibration, the modeller faces challenges such as poor parameter identifiability, lack of sufficiently informative experimental data and the existence of local minima in the objective function landscape. These issues tend to worsen with larger model sizes, increasing the computational complexity and the number of unknown parameters. An incorrectly calibrated model is problematic because it may result in inaccurate predictions and misleading conclusions. For nonexpert users, there are a large number of potential pitfalls. Here, we provide a protocol that guides the user through all the steps involved in the calibration of dynamic models. We illustrate the methodology with two models and provide all the code required to reproduce the results and perform the same analysis on new models. Our protocol provides practitioners and researchers in biological modelling with a one-stop guide that is at the same time compact and sufficiently comprehensive to cover all aspects of the problem.
Collapse
Affiliation(s)
- Alejandro F Villaverde
- Universidade de Vigo, Department of Systems Engineering & Control, Vigo 36310, Galicia, Spain
| | - Dilan Pathirana
- Faculty of Mathematics and Natural Sciences, University of Bonn, Bonn 53115, Germany
| | - Fabian Fröhlich
- Institute of Computational Biology, Helmholtz Zentrum München, Neuherberg 85764, Germany
| | - Jan Hasenauer
- Center for Mathematics, Technische Universität München, Garching 85748, Germany
- Harvard Medical School, Cambridge, MA 02115, USA
| | - Julio R Banga
- Bioprocess Engineering Group, IIM-CSIC, Vigo 36208, Galicia, Spain
| |
Collapse
|
30
|
Lao-Martil D, Verhagen KJA, Schmitz JPJ, Teusink B, Wahl SA, van Riel NAW. Kinetic Modeling of Saccharomyces cerevisiae Central Carbon Metabolism: Achievements, Limitations, and Opportunities. Metabolites 2022; 12:74. [PMID: 35050196 PMCID: PMC8779790 DOI: 10.3390/metabo12010074] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2021] [Revised: 01/11/2022] [Accepted: 01/12/2022] [Indexed: 11/23/2022] Open
Abstract
Central carbon metabolism comprises the metabolic pathways in the cell that process nutrients into energy, building blocks and byproducts. To unravel the regulation of this network upon glucose perturbation, several metabolic models have been developed for the microorganism Saccharomyces cerevisiae. These dynamic representations have focused on glycolysis and answered multiple research questions, but no commonly applicable model has been presented. This review systematically evaluates the literature to describe the current advances, limitations, and opportunities. Different kinetic models have unraveled key kinetic glycolytic mechanisms. Nevertheless, some uncertainties regarding model topology and parameter values still limit the application to specific cases. Progressive improvements in experimental measurement technologies as well as advances in computational tools create new opportunities to further extend the model scale. Notably, models need to be made more complex to consider the multiple layers of glycolytic regulation and external physiological variables regulating the bioprocess, opening new possibilities for extrapolation and validation. Finally, the onset of new data representative of individual cells will cause these models to evolve from depicting an average cell in an industrial fermenter, to characterizing the heterogeneity of the population, opening new and unseen possibilities for industrial fermentation improvement.
Collapse
Affiliation(s)
- David Lao-Martil
- Department of Biomedical Engineering, Eindhoven University of Technology, Groene Loper 5, 5612 AE Eindhoven, The Netherlands;
| | - Koen J. A. Verhagen
- Lehrstuhl für Bioverfahrenstechnik, FAU Erlangen-Nürnberg, 91052 Erlangen, Germany; (K.J.A.V.); (S.A.W.)
| | - Joep P. J. Schmitz
- DSM Biotechnology Center, Alexander Fleminglaan 1, 2613 AX Delft, The Netherlands;
| | - Bas Teusink
- Systems Biology Lab, Amsterdam Institute of Molecular and Life Sciences, Vrije Universiteit Amsterdam, 1081 HZ Amsterdam, The Netherlands;
| | - S. Aljoscha Wahl
- Lehrstuhl für Bioverfahrenstechnik, FAU Erlangen-Nürnberg, 91052 Erlangen, Germany; (K.J.A.V.); (S.A.W.)
| | - Natal A. W. van Riel
- Department of Biomedical Engineering, Eindhoven University of Technology, Groene Loper 5, 5612 AE Eindhoven, The Netherlands;
- Amsterdam University Medical Center, University of Amsterdam, Meibergdreef 9, 1105 AZ Amsterdam, The Netherlands
| |
Collapse
|
31
|
Stapor P, Schmiester L, Wierling C, Merkt S, Pathirana D, Lange BMH, Weindl D, Hasenauer J. Mini-batch optimization enables training of ODE models on large-scale datasets. Nat Commun 2022; 13:34. [PMID: 35013141 PMCID: PMC8748893 DOI: 10.1038/s41467-021-27374-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2019] [Accepted: 11/11/2021] [Indexed: 11/09/2022] Open
Abstract
Quantitative dynamic models are widely used to study cellular signal processing. A critical step in modelling is the estimation of unknown model parameters from experimental data. As model sizes and datasets are steadily growing, established parameter optimization approaches for mechanistic models become computationally extremely challenging. Mini-batch optimization methods, as employed in deep learning, have better scaling properties. In this work, we adapt, apply, and benchmark mini-batch optimization for ordinary differential equation (ODE) models, thereby establishing a direct link between dynamic modelling and machine learning. On our main application example, a large-scale model of cancer signaling, we benchmark mini-batch optimization against established methods, achieving better optimization results and reducing computation by more than an order of magnitude. We expect that our work will serve as a first step towards mini-batch optimization tailored to ODE models and enable modelling of even larger and more complex systems than what is currently possible.
Collapse
Affiliation(s)
- Paul Stapor
- Helmholtz Zentrum München - German Research Center for Environmental Health, Institute of Computational Biology, 85764, Neuherberg, Germany
- Technische Universität München, Center for Mathematics, Chair of Mathematical Modeling of Biological Systems, 85748, Garching, Germany
| | - Leonard Schmiester
- Helmholtz Zentrum München - German Research Center for Environmental Health, Institute of Computational Biology, 85764, Neuherberg, Germany
- Technische Universität München, Center for Mathematics, Chair of Mathematical Modeling of Biological Systems, 85748, Garching, Germany
| | | | - Simon Merkt
- Universität Bonn, Faculty of Mathematics and Natural Sciences, 53115, Bonn, Germany
| | - Dilan Pathirana
- Universität Bonn, Faculty of Mathematics and Natural Sciences, 53115, Bonn, Germany
| | | | - Daniel Weindl
- Helmholtz Zentrum München - German Research Center for Environmental Health, Institute of Computational Biology, 85764, Neuherberg, Germany
| | - Jan Hasenauer
- Helmholtz Zentrum München - German Research Center for Environmental Health, Institute of Computational Biology, 85764, Neuherberg, Germany.
- Technische Universität München, Center for Mathematics, Chair of Mathematical Modeling of Biological Systems, 85748, Garching, Germany.
- Universität Bonn, Faculty of Mathematics and Natural Sciences, 53115, Bonn, Germany.
| |
Collapse
|
32
|
Reyes BC, Otero-Muras I, Petyuk VA. A numerical approach for detecting switch-like bistability in mass action chemical reaction networks with conservation laws. BMC Bioinformatics 2022; 23:1. [PMID: 34983366 PMCID: PMC8725470 DOI: 10.1186/s12859-021-04477-x] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2020] [Accepted: 11/11/2021] [Indexed: 12/24/2022] Open
Abstract
BACKGROUND Theoretical analysis of signaling pathways can provide a substantial amount of insight into their function. One particular area of research considers signaling pathways capable of assuming two or more stable states given the same amount of signaling ligand. This phenomenon of bistability can give rise to switch-like behavior, a mechanism that governs cellular decision making. Investigation of whether or not a signaling pathway can confer bistability and switch-like behavior, without knowledge of specific kinetic rate constant values, is a mathematically challenging problem. Recently a technique based on optimization has been introduced, which is capable of finding example parameter values that confer switch-like behavior for a given pathway. Although this approach has made it possible to analyze moderately sized pathways, it is limited to reaction networks that presume a uniterminal structure. It is this limited structure we address by developing a general technique that applies to any mass action reaction network with conservation laws. RESULTS In this paper we developed a generalized method for detecting switch-like bistable behavior in any mass action reaction network with conservation laws. The method involves (1) construction of a constrained optimization problem using the determinant of the Jacobian of the underlying rate equations, (2) minimization of the objective function to search for conditions resulting in a zero eigenvalue, (3) computation of a confidence level that describes if the global minimum has been found and (4) evaluation of optimization values, using either numerical continuation or directly simulating the ODE system, to verify that a bistability region exists. The generalized method has been tested on three motifs known to be capable of bistability. CONCLUSIONS We have developed a variation of an optimization-based method for the discovery of bistability, which is not limited to uniterminal chemical reaction networks. Successful completion of the method provides an S-shaped bifurcation diagram, which indicates that the network acts as a bistable switch for the given optimization parameters.
Collapse
Affiliation(s)
- Brandon C Reyes
- Advanced Computing, Math, and Data Division, Pacific Northwest National Laboratory, Richland, WA, 99352, USA
| | - Irene Otero-Muras
- BioProcess Engineering Group, IIM-CSIC (Spanish National Research Council), 36208, Vigo, Spain
| | - Vladislav A Petyuk
- Biological Sciences Division, Pacific Northwest National Laboratory, Richland, WA, 99352, USA.
| |
Collapse
|
33
|
Zhang Y, Wang H, Oliveira RHM, Zhao C, Popel AS. Systems biology of angiogenesis signaling: Computational models and omics. WIREs Mech Dis 2021; 14:e1550. [PMID: 34970866 PMCID: PMC9243197 DOI: 10.1002/wsbm.1550] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2021] [Revised: 12/03/2021] [Accepted: 12/06/2021] [Indexed: 01/10/2023]
Abstract
Angiogenesis is a highly regulated multiscale process that involves a plethora of cells, their cellular signal transduction, activation, proliferation, differentiation, as well as their intercellular communication. The coordinated execution and integration of such complex signaling programs is critical for physiological angiogenesis to take place in normal growth, development, exercise, and wound healing, while its dysregulation is critically linked to many major human diseases such as cancer, cardiovascular diseases, and ocular disorders; it is also crucial in regenerative medicine. Although huge efforts have been devoted to drug development for these diseases by investigation of angiogenesis‐targeted therapies, only a few therapeutics and targets have proved effective in humans due to the innate multiscale complexity and nonlinearity in the process of angiogenic signaling. As a promising approach that can help better address this challenge, systems biology modeling allows the integration of knowledge across studies and scales and provides a powerful means to mechanistically elucidate and connect the individual molecular and cellular signaling components that function in concert to regulate angiogenesis. In this review, we summarize and discuss how systems biology modeling studies, at the pathway‐, cell‐, tissue‐, and whole body‐levels, have advanced our understanding of signaling in angiogenesis and thereby delivered new translational insights for human diseases. This article is categorized under:Cardiovascular Diseases > Computational Models Cancer > Computational Models
Collapse
Affiliation(s)
- Yu Zhang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Rebeca Hannah M Oliveira
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Chen Zhao
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA.,School of Pharmacy, Nanjing Medical University, Nanjing, Jiangsu, China
| | - Aleksander S Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| |
Collapse
|
34
|
Schmiester L, Weindl D, Hasenauer J. Efficient gradient-based parameter estimation for dynamic models using qualitative data. BIOINFORMATICS (OXFORD, ENGLAND) 2021. [PMID: 34260697 DOI: 10.1101/2021.02.06.430039] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
MOTIVATION Unknown parameters of dynamical models are commonly estimated from experimental data. However, while various efficient optimization and uncertainty analysis methods have been proposed for quantitative data, methods for qualitative data are rare and suffer from bad scaling and convergence. RESULTS Here, we propose an efficient and reliable framework for estimating the parameters of ordinary differential equation models from qualitative data. In this framework, we derive a semi-analytical algorithm for gradient calculation of the optimal scaling method developed for qualitative data. This enables the use of efficient gradient-based optimization algorithms. We demonstrate that the use of gradient information improves performance of optimization and uncertainty quantification on several application examples. On average, we achieve a speedup of more than one order of magnitude compared to gradient-free optimization. In addition, in some examples, the gradient-based approach yields substantially improved objective function values and quality of the fits. Accordingly, the proposed framework substantially improves the parameterization of models from qualitative data. AVAILABILITY AND IMPLEMENTATION The proposed approach is implemented in the open-source Python Parameter EStimation TOolbox (pyPESTO). pyPESTO is available at https://github.com/ICB-DCM/pyPESTO. All application examples and code to reproduce this study are available at https://doi.org/10.5281/zenodo.4507613. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Leonard Schmiester
- Institute of Computational Biology, Helmholtz Zentrum München - German Research Center for Environmental Health, Neuherberg 85764, Germany
- Center for Mathematics, Technische Universität München, Garching 85748, Germany
| | - Daniel Weindl
- Institute of Computational Biology, Helmholtz Zentrum München - German Research Center for Environmental Health, Neuherberg 85764, Germany
| | - Jan Hasenauer
- Institute of Computational Biology, Helmholtz Zentrum München - German Research Center for Environmental Health, Neuherberg 85764, Germany
- Center for Mathematics, Technische Universität München, Garching 85748, Germany
- Faculty of Mathematics and Natural Sciences, University of Bonn, Bonn 53113, Germany
| |
Collapse
|
35
|
Liebermeister W, Noor E. Model Balancing: A Search for In-Vivo Kinetic Constants and Consistent Metabolic States. Metabolites 2021; 11:749. [PMID: 34822407 PMCID: PMC8621975 DOI: 10.3390/metabo11110749] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2021] [Revised: 10/08/2021] [Accepted: 10/09/2021] [Indexed: 11/16/2022] Open
Abstract
Enzyme kinetic constants in vivo are largely unknown, which limits the construction of large metabolic models. Given measured metabolic fluxes, metabolite concentrations, and enzyme concentrations, these constants may be inferred by model fitting, but the estimation problems are hard to solve if models are large. Here we show how consistent kinetic constants, metabolite concentrations, and enzyme concentrations can be determined from data if metabolic fluxes are known. The estimation method, called model balancing, can handle models with a wide range of rate laws and accounts for thermodynamic constraints between fluxes, kinetic constants, and metabolite concentrations. It can be used to estimate in-vivo kinetic constants, to complete and adjust available data, and to construct plausible metabolic states with predefined flux distributions. By omitting one term from the log posterior-a term for penalising low enzyme concentrations-we obtain a convex optimality problem with a unique local optimum. As a demonstrative case, we balance a model of E. coli central metabolism with artificial or experimental data and obtain a physically and biologically plausible parameterisation of reaction kinetics in E. coli central metabolism. The example shows what information about kinetic constants can be obtained from omics data and reveals practical limits to estimating in-vivo kinetic constants. While noise-free omics data allow for a reasonable reconstruction of in-vivo kcat and KM values, prediction from noisy omics data are worse. Hence, adjusting kinetic constants and omics data to obtain consistent metabolic models is the main application of model balancing.
Collapse
Affiliation(s)
| | - Elad Noor
- Department of Plant and Environmental Sciences, Weizmann Institute of Science, Rehovot 7610001, Israel;
| |
Collapse
|
36
|
Gaiewski MJ, Drewell RA, Dresch JM. Fitting thermodynamic-based models: Incorporating parameter sensitivity improves the performance of an evolutionary algorithm. Math Biosci 2021; 342:108716. [PMID: 34687735 DOI: 10.1016/j.mbs.2021.108716] [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: 02/03/2021] [Revised: 09/10/2021] [Accepted: 09/17/2021] [Indexed: 11/30/2022]
Abstract
A detailed comprehension of transcriptional regulation is critical to understanding the genetic control of development and disease across many different organisms. To more fully investigate the complex molecular interactions controlling the precise expression of genes, many groups have constructed mathematical models to complement their experimental approaches. A critical step in such studies is choosing the most appropriate parameter estimation algorithm to enable detailed analysis of the parameters that contribute to the models. In this study, we develop a novel set of evolutionary algorithms that use a pseudo-random Sobol Set to construct the initial population and incorporate parameter sensitivities into the adaptation of mutation rates, using local, global, and hybrid strategies. Comparison of the performance of these new algorithms to a number of current state-of-the-art global parameter estimation algorithms on a range of continuous test functions, as well as synthetic biological data representing models of gene regulatory systems, reveals improved performance of the new algorithms in terms of runtime, error and reproducibility. In addition, by analyzing the ability of these algorithms to fit datasets of varying quality, we provide the experimentalist with a guide to how the algorithms perform across a range of noisy data. These results demonstrate the improved performance of the new set of parameter estimation algorithms and facilitate meaningful integration of model parameters and predictions in our understanding of the molecular mechanisms of gene regulation.
Collapse
Affiliation(s)
- Michael J Gaiewski
- Department of Mathematics and Computer Science, Clark University, Worcester, MA, USA; Department of Mathematics, University of Connecticut, Storrs, CT, USA.
| | | | | |
Collapse
|
37
|
Benchmarking Optimisation Methods for Model Selection and Parameter Estimation of Nonlinear Systems. VIBRATION 2021. [DOI: 10.3390/vibration4030036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Characterisation and quantification of nonlinearities in the engineering structures include selecting and fitting a good mathematical model to a set of experimental vibration data with significant nonlinear features. These tasks involve solving an optimisation problem where it is difficult to choose a priori the best optimisation technique. This paper presents a systematic comparison of ten optimisation methods used to select the best nonlinear model and estimate its parameters through nonlinear system identification. The model selection framework fits the structure’s equation of motions using time-domain dynamic response data and takes into account couplings due to the presence of the nonlinearities. Three benchmark problems are used to evaluate the performance of two families of optimisation methods: (i) deterministic local searches and (ii) global optimisation metaheuristics. Furthermore, hybrid local–global optimisation methods are examined. All benchmark problems include a free play nonlinearity commonly found in mechanical structures. Multiple performance criteria are considered based on computational efficiency and robustness, that is, finding the best nonlinear model. Results show that hybrid methods, that is, the multi-start strategy with local gradient-based Levenberg–Marquardt method and the particle swarm with Levenberg–Marquardt method, lead to a successful selection of nonlinear models and an accurate estimation of their parameters within acceptable computational times.
Collapse
|
38
|
Introducing Parameter Clustering to the OED Procedure for Model Calibration of a Synthetic Inducible Promoter in S. cerevisiae. Processes (Basel) 2021. [DOI: 10.3390/pr9061053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
In recent years, synthetic gene circuits for adding new cell features have become one of the most powerful tools in biological and pharmaceutical research and development. However, because of the inherent non-linearity and noisy experimental data, the experiment-based model calibration of these synthetic parts is perceived as a laborious and time-consuming procedure. Although the optimal experimental design (OED) based on the Fisher information matrix (FIM) has been proved to be an effective means to improve the calibration efficiency, the required calculation increases dramatically with the model size (parameter number). To reduce the OED complexity without losing the calibration accuracy, this paper proposes two OED approaches with different parameter clustering methods and validates the accuracy of calibrated models with in-silico experiments. A model of an inducible synthetic promoter in S. cerevisiae is adopted for bench-marking. The comparison with the traditional off-line OED approach suggests that the OED approaches with both of the clustering methods significantly reduce the complexity of OED problems (for at least 49.0%), while slightly improving the calibration accuracy (11.8% and 19.6% lower estimation error in average for FIM-based and sensitivity-based approaches). This study implicates that for calibrating non-linear models of biological pathways, cluster-based OED could be a beneficial approach to improve the efficiency of optimal experimental design.
Collapse
|
39
|
|
40
|
van Riel NAW, Tiemann CA, Hilbers PAJ, Groen AK. Metabolic Modeling Combined With Machine Learning Integrates Longitudinal Data and Identifies the Origin of LXR-Induced Hepatic Steatosis. Front Bioeng Biotechnol 2021; 8:536957. [PMID: 33665185 PMCID: PMC7921164 DOI: 10.3389/fbioe.2020.536957] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2020] [Accepted: 12/16/2020] [Indexed: 11/23/2022] Open
Abstract
Temporal multi-omics data can provide information about the dynamics of disease development and therapeutic response. However, statistical analysis of high-dimensional time-series data is challenging. Here we develop a novel approach to model temporal metabolomic and transcriptomic data by combining machine learning with metabolic models. ADAPT (Analysis of Dynamic Adaptations in Parameter Trajectories) performs metabolic trajectory modeling by introducing time-dependent parameters in differential equation models of metabolic systems. ADAPT translates structural uncertainty in the model, such as missing information about regulation, into a parameter estimation problem that is solved by iterative learning. We have now extended ADAPT to include both metabolic and transcriptomic time-series data by introducing a regularization function in the learning algorithm. The ADAPT learning algorithm was (re)formulated as a multi-objective optimization problem in which the estimation of trajectories of metabolic parameters is constrained by the metabolite data and refined by gene expression data. ADAPT was applied to a model of hepatic lipid and plasma lipoprotein metabolism to predict metabolic adaptations that are induced upon pharmacological treatment of mice by a Liver X receptor (LXR) agonist. We investigated the excessive accumulation of triglycerides (TG) in the liver resulting in the development of hepatic steatosis. ADAPT predicted that hepatic TG accumulation after LXR activation originates for 80% from an increased influx of free fatty acids. The model also correctly estimated that TG was stored in the cytosol rather than transferred to nascent very-low density lipoproteins. Through model-based integration of temporal metabolic and gene expression data we discovered that increased free fatty acid influx instead of de novo lipogenesis is the main driver of LXR-induced hepatic steatosis. This study illustrates how ADAPT provides estimates for biomedically important parameters that cannot be measured directly, explaining (side-)effects of pharmacological treatment with LXR agonists.
Collapse
Affiliation(s)
- Natal A W van Riel
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, Netherlands.,Department of Vascular Medicine, Amsterdam UMC, Amsterdam, Netherlands.,Maastricht Centre for Systems Biology, Maastricht University, Maastricht, Netherlands
| | - Christian A Tiemann
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, Netherlands
| | - Peter A J Hilbers
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, Netherlands
| | - Albert K Groen
- Department of Vascular Medicine, Amsterdam UMC, Amsterdam, Netherlands.,Department of Laboratory Medicine, University Medical Center Groningen, Groningen, Netherlands
| |
Collapse
|
41
|
Städter P, Schälte Y, Schmiester L, Hasenauer J, Stapor PL. Benchmarking of numerical integration methods for ODE models of biological systems. Sci Rep 2021; 11:2696. [PMID: 33514831 PMCID: PMC7846608 DOI: 10.1038/s41598-021-82196-2] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2020] [Accepted: 01/08/2021] [Indexed: 11/09/2022] Open
Abstract
Ordinary differential equation (ODE) models are a key tool to understand complex mechanisms in systems biology. These models are studied using various approaches, including stability and bifurcation analysis, but most frequently by numerical simulations. The number of required simulations is often large, e.g., when unknown parameters need to be inferred. This renders efficient and reliable numerical integration methods essential. However, these methods depend on various hyperparameters, which strongly impact the ODE solution. Despite this, and although hundreds of published ODE models are freely available in public databases, a thorough study that quantifies the impact of hyperparameters on the ODE solver in terms of accuracy and computation time is still missing. In this manuscript, we investigate which choices of algorithms and hyperparameters are generally favorable when dealing with ODE models arising from biological processes. To ensure a representative evaluation, we considered 142 published models. Our study provides evidence that most ODEs in computational biology are stiff, and we give guidelines for the choice of algorithms and hyperparameters. We anticipate that our results will help researchers in systems biology to choose appropriate numerical methods when dealing with ODE models.
Collapse
Affiliation(s)
- Philipp Städter
- Institute of Computational Biology, Helmholtz Zentrum München - German Research Center for Environmental Health, 85764, Neuherberg, Germany
- Center for Mathematics, Technische Universität München, 85748, Garching, Germany
| | - Yannik Schälte
- Institute of Computational Biology, Helmholtz Zentrum München - German Research Center for Environmental Health, 85764, Neuherberg, Germany
- Center for Mathematics, Technische Universität München, 85748, Garching, Germany
| | - Leonard Schmiester
- Institute of Computational Biology, Helmholtz Zentrum München - German Research Center for Environmental Health, 85764, Neuherberg, Germany
- Center for Mathematics, Technische Universität München, 85748, Garching, Germany
| | - Jan Hasenauer
- Institute of Computational Biology, Helmholtz Zentrum München - German Research Center for Environmental Health, 85764, Neuherberg, Germany.
- Center for Mathematics, Technische Universität München, 85748, Garching, Germany.
- Faculty of Mathematics and Natural Sciences, University of Bonn, 53113, Bonn, Germany.
| | - Paul L Stapor
- Institute of Computational Biology, Helmholtz Zentrum München - German Research Center for Environmental Health, 85764, Neuherberg, Germany
- Center for Mathematics, Technische Universität München, 85748, Garching, Germany
| |
Collapse
|
42
|
Balsa-Canto E, Bandiera L, Menolascina F. Optimal Experimental Design for Systems and Synthetic Biology Using AMIGO2. Methods Mol Biol 2021; 2229:221-239. [PMID: 33405225 DOI: 10.1007/978-1-0716-1032-9_11] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Dynamic modeling in systems and synthetic biology is still quite a challenge-the complex nature of the interactions results in nonlinear models, which include unknown parameters (or functions). Ideally, time-series data support the estimation of model unknowns through data fitting. Goodness-of-fit measures would lead to the best model among a set of candidates. However, even when state-of-the-art measuring techniques allow for an unprecedented amount of data, not all data suit dynamic modeling.Model-based optimal experimental design (OED) is intended to improve model predictive capabilities. OED can be used to define the set of experiments that would (a) identify the best model or (b) improve the identifiability of unknown parameters. In this chapter, we present a detailed practical procedure to compute optimal experiments using the AMIGO2 toolbox.
Collapse
Affiliation(s)
- Eva Balsa-Canto
- (Bio)Process Engineering Group, IIM-CSIC (Spanish National Research Council), Vigo, Spain.
| | - Lucia Bandiera
- School of Engineering, Institute for Bioengineering, The University of Edinburgh, Edinburgh, UK.,SynthSys - Centre for Synthetic and Systems Biology, The University of Edinburgh, Edinburgh, UK
| | - Filippo Menolascina
- School of Engineering, Institute for Bioengineering, The University of Edinburgh, Edinburgh, UK.,SynthSys - Centre for Synthetic and Systems Biology, The University of Edinburgh, Edinburgh, UK
| |
Collapse
|
43
|
Yazdani A, Lu L, Raissi M, Karniadakis GE. Systems biology informed deep learning for inferring parameters and hidden dynamics. PLoS Comput Biol 2020; 16:e1007575. [PMID: 33206658 PMCID: PMC7710119 DOI: 10.1371/journal.pcbi.1007575] [Citation(s) in RCA: 48] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2019] [Revised: 12/02/2020] [Accepted: 10/11/2020] [Indexed: 01/23/2023] Open
Abstract
Mathematical models of biological reactions at the system-level lead to a set of ordinary differential equations with many unknown parameters that need to be inferred using relatively few experimental measurements. Having a reliable and robust algorithm for parameter inference and prediction of the hidden dynamics has been one of the core subjects in systems biology, and is the focus of this study. We have developed a new systems-biology-informed deep learning algorithm that incorporates the system of ordinary differential equations into the neural networks. Enforcing these equations effectively adds constraints to the optimization procedure that manifests itself as an imposed structure on the observational data. Using few scattered and noisy measurements, we are able to infer the dynamics of unobserved species, external forcing, and the unknown model parameters. We have successfully tested the algorithm for three different benchmark problems.
Collapse
Affiliation(s)
- Alireza Yazdani
- Division of Applied Mathematics, Brown University, Providence, Rhode Island, USA
| | - Lu Lu
- Department of Mathematics, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
| | - Maziar Raissi
- Department of Applied Mathematics, University of Colorado, Boulder, Colorado, USA
| | | |
Collapse
|
44
|
Schmiester L, Weindl D, Hasenauer J. Parameterization of mechanistic models from qualitative data using an efficient optimal scaling approach. J Math Biol 2020; 81:603-623. [PMID: 32696085 PMCID: PMC7427713 DOI: 10.1007/s00285-020-01522-w] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2019] [Revised: 05/28/2020] [Indexed: 12/21/2022]
Abstract
Quantitative dynamical models facilitate the understanding of biological processes and the prediction of their dynamics. These models usually comprise unknown parameters, which have to be inferred from experimental data. For quantitative experimental data, there are several methods and software tools available. However, for qualitative data the available approaches are limited and computationally demanding. Here, we consider the optimal scaling method which has been developed in statistics for categorical data and has been applied to dynamical systems. This approach turns qualitative variables into quantitative ones, accounting for constraints on their relation. We derive a reduced formulation for the optimization problem defining the optimal scaling. The reduced formulation possesses the same optimal points as the established formulation but requires less degrees of freedom. Parameter estimation for dynamical models of cellular pathways revealed that the reduced formulation improves the robustness and convergence of optimizers. This resulted in substantially reduced computation times. We implemented the proposed approach in the open-source Python Parameter EStimation TOolbox (pyPESTO) to facilitate reuse and extension. The proposed approach enables efficient parameterization of quantitative dynamical models using qualitative data.
Collapse
Affiliation(s)
- Leonard Schmiester
- Institute of Computational Biology, Helmholtz Zentrum München–German Research Center for Environmental Health, 85764 Neuherberg, Germany
- Center for Mathematics, Technische Universität München, 85748 Garching, Germany
| | - Daniel Weindl
- Institute of Computational Biology, Helmholtz Zentrum München–German Research Center for Environmental Health, 85764 Neuherberg, Germany
| | - Jan Hasenauer
- Institute of Computational Biology, Helmholtz Zentrum München–German Research Center for Environmental Health, 85764 Neuherberg, Germany
- Center for Mathematics, Technische Universität München, 85748 Garching, Germany
- Faculty of Mathematics and Natural Sciences, University of Bonn, 53113 Bonn, Germany
| |
Collapse
|
45
|
Hass H, Loos C, Raimúndez-Álvarez E, Timmer J, Hasenauer J, Kreutz C. Benchmark problems for dynamic modeling of intracellular processes. Bioinformatics 2020; 35:3073-3082. [PMID: 30624608 PMCID: PMC6735869 DOI: 10.1093/bioinformatics/btz020] [Citation(s) in RCA: 42] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2018] [Revised: 11/19/2018] [Accepted: 01/06/2019] [Indexed: 12/19/2022] Open
Abstract
Motivation Dynamic models are used in systems biology to study and understand cellular processes like gene regulation or signal transduction. Frequently, ordinary differential equation (ODE) models are used to model the time and dose dependency of the abundances of molecular compounds as well as interactions and translocations. A multitude of computational approaches, e.g. for parameter estimation or uncertainty analysis have been developed within recent years. However, many of these approaches lack proper testing in application settings because a comprehensive set of benchmark problems is yet missing. Results We present a collection of 20 benchmark problems in order to evaluate new and existing methodologies, where an ODE model with corresponding experimental data is referred to as problem. In addition to the equations of the dynamical system, the benchmark collection provides observation functions as well as assumptions about measurement noise distributions and parameters. The presented benchmark models comprise problems of different size, complexity and numerical demands. Important characteristics of the models and methodological requirements are summarized, estimated parameters are provided, and some example studies were performed for illustrating the capabilities of the presented benchmark collection. Availability and implementation The models are provided in several standardized formats, including an easy-to-use human readable form and machine-readable SBML files. The data is provided as Excel sheets. All files are available at https://github.com/Benchmarking-Initiative/Benchmark-Models, including step-by-step explanations and MATLAB code to process and simulate the models. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Helge Hass
- Center for Systems Biology (ZBSA), University of Freiburg, Freiburg 79104, Germany.,Institute of Physics, University of Freiburg, Freiburg 79104, Germany
| | - Carolin Loos
- Helmholtz Zentrum München - German Research Center for Environmental Health, Institute of Computational Biology, Neuherberg 85764, Germany.,Technische Universität München, Center for Mathematics, Chair of Mathematical Modeling of Biological Systems, Garching 85748, Germany
| | - Elba Raimúndez-Álvarez
- Helmholtz Zentrum München - German Research Center for Environmental Health, Institute of Computational Biology, Neuherberg 85764, Germany.,Technische Universität München, Center for Mathematics, Chair of Mathematical Modeling of Biological Systems, Garching 85748, Germany
| | - Jens Timmer
- Center for Systems Biology (ZBSA), University of Freiburg, Freiburg 79104, Germany.,Institute of Physics, University of Freiburg, Freiburg 79104, Germany.,Center for Data Analysis and Modelling (FDM), University of Freiburg, Freiburg 79104, Germany.,BIOSS Centre for Biological Signalling Studies, University of Freiburg, Freiburg 79104, Germany
| | - Jan Hasenauer
- Helmholtz Zentrum München - German Research Center for Environmental Health, Institute of Computational Biology, Neuherberg 85764, Germany.,Technische Universität München, Center for Mathematics, Chair of Mathematical Modeling of Biological Systems, Garching 85748, Germany
| | - Clemens Kreutz
- Center for Systems Biology (ZBSA), University of Freiburg, Freiburg 79104, Germany.,Institute of Physics, University of Freiburg, Freiburg 79104, Germany.,Center for Data Analysis and Modelling (FDM), University of Freiburg, Freiburg 79104, Germany
| |
Collapse
|
46
|
Rodman AD, Gerogiorgis DI. Parameter estimation and sensitivity analysis for dynamic modelling and simulation of beer fermentation. Comput Chem Eng 2020. [DOI: 10.1016/j.compchemeng.2019.106665] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
47
|
Evans MV, Eklund CR, Williams DN, Sey YM, Simmons JE. Global optimization of the Michaelis-Menten parameters using physiologically-based pharmacokinetic (PBPK) modeling and chloroform vapor uptake data in F344 rats. Inhal Toxicol 2020; 32:97-109. [PMID: 32241199 DOI: 10.1080/08958378.2020.1742818] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Abstract
Objective: To quantify metabolism, a physiologically based pharmacokinetic (PBPK) model for a volatile compound can be calibrated with the closed chamber (i.e. vapor uptake) inhalation data. Here, we introduce global optimization as a novel component of the predictive process and use it to illustrate a procedure for metabolic parameter estimation.Materials and methods: Male F344 rats were exposed in vapor uptake chambers to initial concentrations of 100, 500, 1000, and 3000 ppm chloroform. Chamber time-course data from these experiments, in combination with optimization using a chemical-specific PBPK model, were used to estimate Michaelis-Menten metabolic constants. Matlab® simulation software was used to integrate the mass balance equations and to perform the global optimizations using MEIGO (MEtaheuristics for systems biology and bIoinformatics Global Optimization - Version 64 bit, R2016A), a toolbox written for Matlab®. The cost function used the chamber time-course data and least squares to minimize the difference between data and simulation values.Results and discussion: The final values estimated for Vmax (maximum metabolic rate) and Km (affinity constant) were 1.2 mg/h and a range between 0.0005 and 0.6 mg/L, respectively. Also, cost function plots were used to analyze the dose-dependent capacity to estimate Vmax and Km within the experimental range used. Sensitivity analysis was used to assess identifiability for both parameters and show these kinetic data may not be sufficient to identify Km.Conclusion: In summary, this work should help toxicologists interested in optimization techniques understand the overall process employed when calibrating metabolic parameters in a PBPK model with inhalation data.
Collapse
Affiliation(s)
- Marina V Evans
- ORD, National Health and Environmental Effects Research Laboratory, ISTD, US Environmental Protection Agency, Research Triangle Park, NC, USA
| | - Christopher R Eklund
- ORD, National Health and Environmental Effects Research Laboratory, ISTD, US Environmental Protection Agency, Research Triangle Park, NC, USA
| | - David N Williams
- ORISE, Oak Ridge Institute for Science and Education, Oak Ridge, TN, USA
| | - Yusupha M Sey
- ORD, National Health and Environmental Effects Research Laboratory, ISTD, US Environmental Protection Agency, Research Triangle Park, NC, USA
| | - Jane Ellen Simmons
- ORD, National Health and Environmental Effects Research Laboratory, ISTD, US Environmental Protection Agency, Research Triangle Park, NC, USA
| |
Collapse
|
48
|
Raimúndez E, Keller S, Zwingenberger G, Ebert K, Hug S, Theis FJ, Maier D, Luber B, Hasenauer J. Model-based analysis of response and resistance factors of cetuximab treatment in gastric cancer cell lines. PLoS Comput Biol 2020; 16:e1007147. [PMID: 32119655 PMCID: PMC7067490 DOI: 10.1371/journal.pcbi.1007147] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2019] [Revised: 03/12/2020] [Accepted: 01/27/2020] [Indexed: 12/31/2022] Open
Abstract
Targeted cancer therapies are powerful alternatives to chemotherapies or can be used complementary to these. Yet, the response to targeted treatments depends on a variety of factors, including mutations and expression levels, and therefore their outcome is difficult to predict. Here, we develop a mechanistic model of gastric cancer to study response and resistance factors for cetuximab treatment. The model captures the EGFR, ERK and AKT signaling pathways in two gastric cancer cell lines with different mutation patterns. We train the model using a comprehensive selection of time and dose response measurements, and provide an assessment of parameter and prediction uncertainties. We demonstrate that the proposed model facilitates the identification of causal differences between the cell lines. Furthermore, our study shows that the model provides predictions for the responses to different perturbations, such as knockdown and knockout experiments. Among other results, the model predicted the effect of MET mutations on cetuximab sensitivity. These predictive capabilities render the model a basis for the assessment of gastric cancer signaling and possibly for the development and discovery of predictive biomarkers. Unraveling the causal differences between drug responders and non-responders is an important challenge. The information can help to understand molecular mechanisms and to guide the selection and design of targeted therapies. Here, we approach this problem for cetuximab treatment for gastric cancer using mechanistic mathematical modeling. The proposed model describes responder and non-responder gastric cancer cell lines and can predict the response in several validation experiments. Our analysis provides a differentiated view on mutations and explains, for instance, the relevance of MET mutations and the insignificance of PIK3CA mutation in the considered cell lines. The model might potentially provide the basis for understanding the recent failure of several clinical studies.
Collapse
Affiliation(s)
- Elba Raimúndez
- Helmholtz Zentrum München-German Research Center for Environmental Health, Institute of Computational Biology, Neuherberg, Germany
- Center for Mathematics, Technische Universität München, Garching, Germany
| | - Simone Keller
- Technical University of Munich, School of Medicine, Klinikum rechts der Isar, Institute of Pathology, Munich, Germany
| | - Gwen Zwingenberger
- Technical University of Munich, School of Medicine, Klinikum rechts der Isar, Institute of Pathology, Munich, Germany
| | - Karolin Ebert
- Technical University of Munich, School of Medicine, Klinikum rechts der Isar, Institute of Pathology, Munich, Germany
| | - Sabine Hug
- Helmholtz Zentrum München-German Research Center for Environmental Health, Institute of Computational Biology, Neuherberg, Germany
- Center for Mathematics, Technische Universität München, Garching, Germany
| | - Fabian J. Theis
- Helmholtz Zentrum München-German Research Center for Environmental Health, Institute of Computational Biology, Neuherberg, Germany
- Center for Mathematics, Technische Universität München, Garching, Germany
| | | | - Birgit Luber
- Technical University of Munich, School of Medicine, Klinikum rechts der Isar, Institute of Pathology, Munich, Germany
| | - Jan Hasenauer
- Helmholtz Zentrum München-German Research Center for Environmental Health, Institute of Computational Biology, Neuherberg, Germany
- Center for Mathematics, Technische Universität München, Garching, Germany
- Faculty of Mathematics and Natural Sciences, University of Bonn, Bonn, Germany
- * E-mail:
| |
Collapse
|
49
|
Kreutz C. Guidelines for benchmarking of optimization-based approaches for fitting mathematical models. Genome Biol 2019; 20:281. [PMID: 31842943 PMCID: PMC6915982 DOI: 10.1186/s13059-019-1887-9] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2019] [Accepted: 11/13/2019] [Indexed: 11/10/2022] Open
Abstract
Insufficient performance of optimization-based approaches for the fitting of mathematical models is still a major bottleneck in systems biology. In this article, the reasons and methodological challenges are summarized as well as their impact in benchmark studies. Important aspects for achieving an increased level of evidence for benchmark results are discussed. Based on general guidelines for benchmarking in computational biology, a collection of tailored guidelines is presented for performing informative and unbiased benchmarking of optimization-based fitting approaches. Comprehensive benchmark studies based on these recommendations are urgently required for the establishment of a robust and reliable methodology for the systems biology community.
Collapse
Affiliation(s)
- Clemens Kreutz
- Faculty of Medicine and Medical Center, Institute of Medical Biometry and Statistics, University of Freiburg, Stefan-Meier-Str. 26, Freiburg, 79104, Germany.
- CIBSS-Centre for Integrative Biological Signalling Studies, University of Freiburg, Freiburg, Germany.
| |
Collapse
|
50
|
Eisenkolb I, Jensch A, Eisenkolb K, Kramer A, Buchholz PCF, Pleiss J, Spiess A, Radde NE. Modeling of biocatalytic reactions: A workflow for model calibration, selection, and validation using Bayesian statistics. AIChE J 2019. [DOI: 10.1002/aic.16866] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Ina Eisenkolb
- Institute for Systems Theory and Automatic ControlUniversity of Stuttgart Stuttgart Germany
| | - Antje Jensch
- Institute for Systems Theory and Automatic ControlUniversity of Stuttgart Stuttgart Germany
| | - Kerstin Eisenkolb
- Institute for Systems Theory and Automatic ControlUniversity of Stuttgart Stuttgart Germany
| | - Andrei Kramer
- Science for Life LaboratoryKTH Royal Institute of Technology Stockholm Sweden
| | - Patrick C. F. Buchholz
- Institute for Biochemistry and Technical BiochemistryUniversity of Stuttgart Stuttgart Germany
| | - Jürgen Pleiss
- Institute for Biochemistry and Technical BiochemistryUniversity of Stuttgart Stuttgart Germany
| | - Antje Spiess
- Institute for Biochemical EngineeringTechnical University Braunschweig Braunschweig Germany
| | - Nicole E. Radde
- Institute for Systems Theory and Automatic ControlUniversity of Stuttgart Stuttgart Germany
| |
Collapse
|