1
|
Tönsing C, Steiert B, Timmer J, Kreutz C. Likelihood-ratio test statistic for the finite-sample case in nonlinear ordinary differential equation models. PLoS Comput Biol 2023; 19:e1011417. [PMID: 37738254 PMCID: PMC10550180 DOI: 10.1371/journal.pcbi.1011417] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Revised: 10/04/2023] [Accepted: 08/08/2023] [Indexed: 09/24/2023] Open
Abstract
Likelihood ratios are frequently utilized as basis for statistical tests, for model selection criteria and for assessing parameter and prediction uncertainties, e.g. using the profile likelihood. However, translating these likelihood ratios into p-values or confidence intervals requires the exact form of the test statistic's distribution. The lack of knowledge about this distribution for nonlinear ordinary differential equation (ODE) models requires an approximation which assumes the so-called asymptotic setting, i.e. a sufficiently large amount of data. Since the amount of data from quantitative molecular biology is typically limited in applications, this finite-sample case regularly occurs for mechanistic models of dynamical systems, e.g. biochemical reaction networks or infectious disease models. Thus, it is unclear whether the standard approach of using statistical thresholds derived for the asymptotic large-sample setting in realistic applications results in valid conclusions. In this study, empirical likelihood ratios for parameters from 19 published nonlinear ODE benchmark models are investigated using a resampling approach for the original data designs. Their distributions are compared to the asymptotic approximation and statistical thresholds are checked for conservativeness. It turns out, that corrections of the likelihood ratios in such finite-sample applications are required in order to avoid anti-conservative results.
Collapse
Affiliation(s)
- Christian Tönsing
- Institute of Physics, University of Freiburg, Germany
- CIBSS Centre for Integrative Biological Signalling Studies, University of Freiburg, Germany
- FDM Freiburg Center for Data Analysis and Modeling, University of Freiburg, Germany
| | | | - Jens Timmer
- Institute of Physics, University of Freiburg, Germany
- CIBSS Centre for Integrative Biological Signalling Studies, University of Freiburg, Germany
- FDM Freiburg Center for Data Analysis and Modeling, University of Freiburg, Germany
| | - Clemens Kreutz
- CIBSS Centre for Integrative Biological Signalling Studies, University of Freiburg, Germany
- FDM Freiburg Center for Data Analysis and Modeling, University of Freiburg, Germany
- Faculty of Medicine and Medical Center, Institute of Medical Biometry and Statistics, University of Freiburg, Germany
| |
Collapse
|
2
|
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
|
3
|
Varshney RK, Bohra A, Roorkiwal M, Barmukh R, Cowling WA, Chitikineni A, Lam HM, Hickey LT, Croser JS, Bayer PE, Edwards D, Crossa J, Weckwerth W, Millar H, Kumar A, Bevan MW, Siddique KHM. Fast-forward breeding for a food-secure world. Trends Genet 2021; 37:1124-1136. [PMID: 34531040 DOI: 10.1016/j.tig.2021.08.002] [Citation(s) in RCA: 61] [Impact Index Per Article: 20.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2021] [Revised: 08/03/2021] [Accepted: 08/04/2021] [Indexed: 10/20/2022]
Abstract
Crop production systems need to expand their outputs sustainably to feed a burgeoning human population. Advances in genome sequencing technologies combined with efficient trait mapping procedures accelerate the availability of beneficial alleles for breeding and research. Enhanced interoperability between different omics and phenotyping platforms, leveraged by evolving machine learning tools, will help provide mechanistic explanations for complex plant traits. Targeted and rapid assembly of beneficial alleles using optimized breeding strategies and precise genome editing techniques could deliver ideal crops for the future. Realizing desired productivity gains in the field is imperative for securing an adequate future food supply for 10 billion people.
Collapse
Affiliation(s)
- Rajeev K Varshney
- Centre of Excellence in Genomics and Systems Biology, International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad 502324, India; State Agricultural Biotechnology Centre, Centre for Crop and Food Innovation, Food Futures Institute, Murdoch University, Murdoch WA 6150, Western Australia, Australia; The UWA Institute of Agriculture, The University of Western Australia, Perth, WA 6009, Australia.
| | - Abhishek Bohra
- ICAR-Indian Institute of Pulses Research (IIPR), Kanpur, India
| | - Manish Roorkiwal
- Centre of Excellence in Genomics and Systems Biology, International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad 502324, India; The UWA Institute of Agriculture, The University of Western Australia, Perth, WA 6009, Australia
| | - Rutwik Barmukh
- Centre of Excellence in Genomics and Systems Biology, International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad 502324, India
| | - Wallace A Cowling
- The UWA Institute of Agriculture, The University of Western Australia, Perth, WA 6009, Australia
| | - Annapurna Chitikineni
- Centre of Excellence in Genomics and Systems Biology, International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad 502324, India
| | - Hon-Ming Lam
- School of Life Sciences and Center for Soybean Research of the State Key Laboratory of Agrobiotechnology, The Chinese University of Hong Kong, Shatin, Hong Kong, China
| | - Lee T Hickey
- Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, QLD, Australia
| | - Janine S Croser
- The UWA Institute of Agriculture, The University of Western Australia, Perth, WA 6009, Australia
| | - Philipp E Bayer
- The UWA Institute of Agriculture, The University of Western Australia, Perth, WA 6009, Australia; School of Biological Sciences, The University of Western Australia, Crawley, WA, Australia
| | - David Edwards
- The UWA Institute of Agriculture, The University of Western Australia, Perth, WA 6009, Australia; School of Biological Sciences, The University of Western Australia, Crawley, WA, Australia
| | - José Crossa
- International Maize and Wheat Improvement Center (CIMMYT), Texcoco, Mexico
| | - Wolfram Weckwerth
- Department of Ecogenomics and Systems Biology, Vienna Metabolomics Center (VIME), University of Vienna, Vienna, Austria
| | - Harvey Millar
- The UWA Institute of Agriculture, The University of Western Australia, Perth, WA 6009, Australia; ARC Centre of Excellence in Plant Energy Biology, The University of Western Australia, Crawley, WA, Australia
| | - Arvind Kumar
- Deputy Director General's Office, International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad 502324, India
| | | | - Kadambot H M Siddique
- The UWA Institute of Agriculture, The University of Western Australia, Perth, WA 6009, Australia
| |
Collapse
|
4
|
Kordes C, Bock HH, Reichert D, May P, Häussinger D. Hepatic stellate cells: current state and open questions. Biol Chem 2021; 402:1021-1032. [PMID: 34008380 DOI: 10.1515/hsz-2021-0180] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2021] [Accepted: 05/03/2021] [Indexed: 01/14/2023]
Abstract
This review article summarizes 20 years of our research on hepatic stellate cells within the framework of two collaborative research centers CRC575 and CRC974 at the Heinrich Heine University. Over this period, stellate cells were identified for the first time as mesenchymal stem cells of the liver, and important functions of these cells in the context of liver regeneration were discovered. Furthermore, it was determined that the space of Disse - bounded by the sinusoidal endothelium and hepatocytes - functions as a stem cell niche for stellate cells. Essential elements of this niche that control the maintenance of hepatic stellate cells have been identified alongside their impairment with age. This article aims to highlight previous studies on stellate cells and critically examine and identify open questions and future research directions.
Collapse
Affiliation(s)
- Claus Kordes
- Clinic of Gastroenterology, Hepatology, and Infectious Diseases, Heinrich Heine University, Moorenstraße 5, D-40225 Düsseldorf, Germany
| | - Hans H Bock
- Clinic of Gastroenterology, Hepatology, and Infectious Diseases, Heinrich Heine University, Moorenstraße 5, D-40225 Düsseldorf, Germany
| | - Doreen Reichert
- Clinic of Gastroenterology, Hepatology, and Infectious Diseases, Heinrich Heine University, Moorenstraße 5, D-40225 Düsseldorf, Germany
| | - Petra May
- Clinic of Gastroenterology, Hepatology, and Infectious Diseases, Heinrich Heine University, Moorenstraße 5, D-40225 Düsseldorf, Germany
| | - Dieter Häussinger
- Clinic of Gastroenterology, Hepatology, and Infectious Diseases, Heinrich Heine University, Moorenstraße 5, D-40225 Düsseldorf, Germany
| |
Collapse
|
5
|
Adlung L, Stapor P, Tönsing C, Schmiester L, Schwarzmüller LE, Postawa L, Wang D, Timmer J, Klingmüller U, Hasenauer J, Schilling M. Cell-to-cell variability in JAK2/STAT5 pathway components and cytoplasmic volumes defines survival threshold in erythroid progenitor cells. Cell Rep 2021; 36:109507. [PMID: 34380040 DOI: 10.1016/j.celrep.2021.109507] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2020] [Revised: 04/01/2021] [Accepted: 07/19/2021] [Indexed: 12/25/2022] Open
Abstract
Survival or apoptosis is a binary decision in individual cells. However, at the cell-population level, a graded increase in survival of colony-forming unit-erythroid (CFU-E) cells is observed upon stimulation with erythropoietin (Epo). To identify components of Janus kinase 2/signal transducer and activator of transcription 5 (JAK2/STAT5) signal transduction that contribute to the graded population response, we extended a cell-population-level model calibrated with experimental data to study the behavior in single cells. The single-cell model shows that the high cell-to-cell variability in nuclear phosphorylated STAT5 is caused by variability in the amount of Epo receptor (EpoR):JAK2 complexes and of SHP1, as well as the extent of nuclear import because of the large variance in the cytoplasmic volume of CFU-E cells. 24-118 pSTAT5 molecules in the nucleus for 120 min are sufficient to ensure cell survival. Thus, variability in membrane-associated processes is sufficient to convert a switch-like behavior at the single-cell level to a graded population-level response.
Collapse
Affiliation(s)
- Lorenz Adlung
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany; Department of Medicine, University Medical Center Hamburg-Eppendorf, 20246 Hamburg, Germany; Hamburg Center for Translational Immunology (HCTI), University Medical Center Hamburg-Eppendorf, 20246 Hamburg, Germany
| | - 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
| | - Christian Tönsing
- Institute of Physics and Freiburg Center for Data Analysis and Modelling (FDM), University of Freiburg, 79104 Freiburg, Germany; CIBSS-Centre for Integrative Biological Signalling Studies, University of Freiburg, 79104 Freiburg, 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
| | - Luisa E Schwarzmüller
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Lena Postawa
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Dantong Wang
- 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
| | - Jens Timmer
- Institute of Physics and Freiburg Center for Data Analysis and Modelling (FDM), University of Freiburg, 79104 Freiburg, Germany; CIBSS-Centre for Integrative Biological Signalling Studies, University of Freiburg, 79104 Freiburg, Germany.
| | - Ursula Klingmüller
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany; Translational Lung Research Center (TLRC), Member of the German Center for Lung Research (DZL), 69120 Heidelberg, 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; Faculty of Mathematics and Natural Sciences, University of Bonn, 53113 Bonn, Germany.
| | - Marcel Schilling
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany.
| |
Collapse
|
6
|
Lee S, Clinedinst L. Mathematical Biology: Expand, Expose, and Educate! Bull Math Biol 2020; 82:120. [PMID: 32910273 DOI: 10.1007/s11538-020-00796-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2020] [Accepted: 08/21/2020] [Indexed: 10/23/2022]
Abstract
Mathematical biology has made significant contributions and advancements in the biological sciences. Recruitment efforts focus on encouraging students, especially those who are underrepresented and underserved, to pursue the field of mathematical biology, regardless of their undergraduate institution type, and raise awareness about the countless professional and academic possibilities provided by this specialized training. This article examines the need to expand, expose, and educate others about mathematical biology. To support field expansion, we give several recommendations of ways to integrate mathematics applied curricula to attract broader student interest. With this exposure-whether it is led by an individual, a department, a university, or researchers in mathematical biology-each can help to promote a base knowledge and appreciation of the field. In order to encourage the next generation of researchers to consider mathematical biology, we highlight current interdisciplinary programs, share popular mathematical tools, and present some thoughts on ways to support a thriving and inclusive mathematical biology community for years to come.
Collapse
Affiliation(s)
- Shernita Lee
- Graduate School, Virginia Tech, Blacksburg, VA, USA.
| | | |
Collapse
|
7
|
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
|