1
|
Dorešić D, Grein S, Hasenauer J. Efficient parameter estimation for ODE models of cellular processes using semi-quantitative data. Bioinformatics 2024; 40:i558-i566. [PMID: 38940161 PMCID: PMC11211815 DOI: 10.1093/bioinformatics/btae210] [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] [Indexed: 06/29/2024] Open
Abstract
MOTIVATION Quantitative dynamical models facilitate the understanding of biological processes and the prediction of their dynamics. The parameters of these models are commonly estimated from experimental data. Yet, experimental data generated from different techniques do not provide direct information about the state of the system but a nonlinear (monotonic) transformation of it. For such semi-quantitative data, when this transformation is unknown, it is not apparent how the model simulations and the experimental data can be compared. RESULTS We propose a versatile spline-based approach for the integration of a broad spectrum of semi-quantitative data into parameter estimation. We derive analytical formulas for the gradients of the hierarchical objective function and show that this substantially increases the estimation efficiency. Subsequently, we demonstrate that the method allows for the reliable discovery of unknown measurement transformations. Furthermore, we show that this approach can significantly improve the parameter inference based on semi-quantitative data in comparison to available methods. AVAILABILITY AND IMPLEMENTATION Modelers can easily apply our method by using our implementation in the open-source Python Parameter EStimation TOolbox (pyPESTO) available at https://github.com/ICB-DCM/pyPESTO.
Collapse
Affiliation(s)
- Domagoj Dorešić
- Life and Medical Sciences (LIMES) Institute, University of Bonn, 53113 Bonn, Germany
- Institute of Computational Biology, Helmholtz Zentrum München – German Research Center for Environmental Health, 85764 Neuherberg, Germany
| | - Stephan Grein
- Life and Medical Sciences (LIMES) Institute, University of Bonn, 53113 Bonn, Germany
| | - Jan Hasenauer
- Life and Medical Sciences (LIMES) Institute, University of Bonn, 53113 Bonn, Germany
- 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
|
2
|
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
|
3
|
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
|
4
|
Faida P, Attiogbe MKI, Majeed U, Zhao J, Qu L, Fan D. Lung cancer treatment potential and limits associated with the STAT family of transcription factors. Cell Signal 2023:110797. [PMID: 37423343 DOI: 10.1016/j.cellsig.2023.110797] [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: 04/06/2023] [Revised: 06/19/2023] [Accepted: 07/04/2023] [Indexed: 07/11/2023]
Abstract
Lung cancer is one of the mortal cancers and the leading cause of cancer-related mortality, with a cancer survival rate of fewer than 5% in developing nations. This low survival rate can be linked to things like late-stage detection, quick postoperative recurrences in patients receiving therapy, and chemoresistance developing against various lung cancer treatments. Signal transducer and activator of transcription (STAT) family of transcription factors are involved in lung cancer cell proliferation, metastasis, immunological control, and treatment resistance. By interacting with specific DNA sequences, STAT proteins trigger the production of particular genes, which in turn result in adaptive and incredibly specific biological responses. In the human genome, seven STAT proteins have been discovered (STAT1 to STAT6, including STAT5a and STAT5b). Many external signaling proteins can activate unphosphorylated STATs (uSTATs), which are found inactively in the cytoplasm. When STAT proteins are activated, they can increase the transcription of several target genes, which leads to unchecked cellular proliferation, anti-apoptotic reactions, and angiogenesis. The effects of STAT transcription factors on lung cancer are variable; some are either pro- or anti-tumorigenic, while others maintain dual, context-dependent activities. Here, we give a succinct summary of the various functions that each member of the STAT family plays in lung cancer and go into more detail about the advantages and disadvantages of pharmacologically targeting STAT proteins and their upstream activators in the context of lung cancer treatment.
Collapse
Affiliation(s)
- Paison Faida
- Shaanxi Key Laboratory of Degradable Biomedical Materials and Shaanxi R&D Center of Biomaterials and Fermentation Engineering, School of Chemical Engineering, Northwest University, Taibai North Road 229, Xi'an, Shaanxi 710069, China; Biotech. & Biomed. Research Institute, Northwest University, Taibai North Road 229, Xi'an, Shaanxi 710069, China
| | - Mawusse K I Attiogbe
- Department of Pharmacology, School of Basic Medical Sciences, Xi'an Jiaotong University, Xi'an, Shaanxi 710061, China
| | - Usman Majeed
- College of Food Science and Technology, Northwest University, Xi'an, Shaanxi 710069, China
| | - Jing Zhao
- Shaanxi Key Laboratory of Degradable Biomedical Materials and Shaanxi R&D Center of Biomaterials and Fermentation Engineering, School of Chemical Engineering, Northwest University, Taibai North Road 229, Xi'an, Shaanxi 710069, China; Biotech. & Biomed. Research Institute, Northwest University, Taibai North Road 229, Xi'an, Shaanxi 710069, China
| | - Linlin Qu
- Shaanxi Key Laboratory of Degradable Biomedical Materials and Shaanxi R&D Center of Biomaterials and Fermentation Engineering, School of Chemical Engineering, Northwest University, Taibai North Road 229, Xi'an, Shaanxi 710069, China; Biotech. & Biomed. Research Institute, Northwest University, Taibai North Road 229, Xi'an, Shaanxi 710069, China
| | - Daidi Fan
- Shaanxi Key Laboratory of Degradable Biomedical Materials and Shaanxi R&D Center of Biomaterials and Fermentation Engineering, School of Chemical Engineering, Northwest University, Taibai North Road 229, Xi'an, Shaanxi 710069, China; Biotech. & Biomed. Research Institute, Northwest University, Taibai North Road 229, Xi'an, Shaanxi 710069, China.
| |
Collapse
|
5
|
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
|
6
|
Chakraborty S, Andrieux G, Kastl P, Adlung L, Altamura S, Boehm ME, Schwarzmüller LE, Abdullah Y, Wagner MC, Helm B, Gröne HJ, Lehmann WD, Boerries M, Busch H, Muckenthaler MU, Schilling M, Klingmüller U. Erythropoietin-driven dynamic proteome adaptations during erythropoiesis prevent iron overload in the developing embryo. Cell Rep 2022; 40:111360. [PMID: 36130519 DOI: 10.1016/j.celrep.2022.111360] [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: 03/14/2022] [Revised: 06/22/2022] [Accepted: 08/23/2022] [Indexed: 11/26/2022] Open
Abstract
Erythropoietin (Epo) ensures survival and proliferation of colony-forming unit erythroid (CFU-E) progenitor cells and their differentiation to hemoglobin-containing mature erythrocytes. A lack of Epo-induced responses causes embryonic lethality, but mechanisms regulating the dynamic communication of cellular alterations to the organismal level remain unresolved. By time-resolved transcriptomics and proteomics, we show that Epo induces in CFU-E cells a gradual transition from proliferation signature proteins to proteins indicative for differentiation, including heme-synthesis enzymes. In the absence of the Epo receptor (EpoR) in embryos, we observe a lack of hemoglobin in CFU-E cells and massive iron overload of the fetal liver pointing to a miscommunication between liver and placenta. A reduction of iron-sulfur cluster-containing proteins involved in oxidative phosphorylation in these embryos leads to a metabolic shift toward glycolysis. This link connecting erythropoiesis with the regulation of iron homeostasis and metabolic reprogramming suggests that balancing these interactions is crucial for protection from iron intoxication and for survival.
Collapse
Affiliation(s)
- Sajib Chakraborty
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany; Systems Cell-Signalling Laboratory, Department of Biochemistry and Molecular Biology, University of Dhaka, Dhaka 1000, Bangladesh
| | - Geoffroy Andrieux
- Institute of Medical Bioinformatics and Systems Medicine, Medical Center-University of Freiburg, Faculty of Medicine, University of Freiburg, 79110 Freiburg, Germany; German Cancer Consortium (DKTK), Freiburg, Germany and German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Philipp Kastl
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Lorenz Adlung
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany; Department of Medicine & Hamburg Center for Translational Immunology, University Medical Center Hamburg-Eppendorf, 20246 Hamburg, Germany
| | - Sandro Altamura
- Center for Translational Biomedical Iron Research (CeTBI), Department of Pediatric Hematology, Oncology and Immunology, Heidelberg University, 69120 Heidelberg, Germany
| | - Martin E Boehm
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Luisa E Schwarzmüller
- Division Molecular Genome Analysis, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Yomn Abdullah
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Marie-Christine Wagner
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Barbara Helm
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Hermann-Josef Gröne
- Division Cellular and Molecular Pathology, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Wolf D Lehmann
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Melanie Boerries
- Institute of Medical Bioinformatics and Systems Medicine, Medical Center-University of Freiburg, Faculty of Medicine, University of Freiburg, 79110 Freiburg, Germany; German Cancer Consortium (DKTK), Freiburg, Germany and German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany; Comprehensive Cancer Center Freiburg (CCCF), Medical Center-University of Freiburg, University of Freiburg, 79106 Freiburg im Breisgau, Germany.
| | - Hauke Busch
- Institute of Medical Bioinformatics and Systems Medicine, Medical Center-University of Freiburg, Faculty of Medicine, University of Freiburg, 79110 Freiburg, Germany; Institute of Experimental Dermatology, University of Lübeck, 23562 Lübeck, Germany.
| | - Martina U Muckenthaler
- Center for Translational Biomedical Iron Research (CeTBI), Department of Pediatric Hematology, Oncology and Immunology, Heidelberg University, 69120 Heidelberg, Germany; Translational Lung Research Center Heidelberg (TLRC), German Center for Lung Research (DZL), 69120 Heidelberg, Germany; German Center for Cardiovascular Research, Partner Site Heidelberg/Mannheim, 69120 Heidelberg, Germany.
| | - Marcel Schilling
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany.
| | - Ursula Klingmüller
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany; Translational Lung Research Center Heidelberg (TLRC), German Center for Lung Research (DZL), 69120 Heidelberg, Germany.
| |
Collapse
|
7
|
Fleury M, Vazquez-Mateo C, Hernandez-Escalante J, Dooms H. Partial STAT5 signaling is sufficient for CD4 + T cell priming but not memory formation. Cytokine 2022; 150:155770. [PMID: 34839177 PMCID: PMC8761165 DOI: 10.1016/j.cyto.2021.155770] [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: 08/23/2021] [Accepted: 11/06/2021] [Indexed: 02/03/2023]
Abstract
Signal transducer and activator of transcription 5 (STAT5) plays an important role in regulating gene expression in response to cytokines of the common (γc) chain family. In this capacity, STAT5 promotes CD8+ effector and memory T cell survival and regulatory T cell development. However, its function in conventional CD4+ T cells is less clear. In this study, the requirement of intact STAT5 signaling for CD4+ effector and memory T cell generation and maintenance was investigated by using DO11.10 TCR transgenic T cells that are genetically deficient in STAT5A or B, as well as by transducing DO11 T cells with a dominant-negative STAT5 to temporally block STAT5 function. We found that the presence of STAT5A or B alone was sufficient for primary CD4+ effector T cell generation, but not for establishing a long-lived memory cell population. Similarly, blocking STAT5 signaling during priming did not prevent initial T cell activation, but inhibited the generation of memory cells. Surprisingly, blocking STAT5 post-priming did not impact the long-term survival of CD4+ memory T cells in vivo. Mechanistically, intact STAT5B, but not STAT5A, was required for IL-7Rα re-expression in activated T cells, which is an important cytokine receptor for CD4+ memory generation. These data show that fully functional STAT5 is essential to deliver an early, non-redundant signal for memory programming during the primary CD4+ T cell response, while partial STAT5 signaling is sufficient for effector differentiation. Our results have implications for the precise use of STAT5 inhibitors to timely inhibit memory T cell responses.
Collapse
Affiliation(s)
- Michelle Fleury
- Arthritis and Autoimmune Diseases Research Center, Rheumatology Section, Department of Medicine, Boston University School of Medicine, Boston MA 02118, United States; Department of Microbiology, Boston University School of Medicine, Boston MA 02118, United States
| | - Cristina Vazquez-Mateo
- Arthritis and Autoimmune Diseases Research Center, Rheumatology Section, Department of Medicine, Boston University School of Medicine, Boston MA 02118, United States
| | - Jaileene Hernandez-Escalante
- Arthritis and Autoimmune Diseases Research Center, Rheumatology Section, Department of Medicine, Boston University School of Medicine, Boston MA 02118, United States; Department of Microbiology, Boston University School of Medicine, Boston MA 02118, United States
| | - Hans Dooms
- Arthritis and Autoimmune Diseases Research Center, Rheumatology Section, Department of Medicine, Boston University School of Medicine, Boston MA 02118, United States; Department of Microbiology, Boston University School of Medicine, Boston MA 02118, United States; Department of Pathology, University of California San Francisco, San Francisco CA 94143, United States.
| |
Collapse
|
8
|
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
|
9
|
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
|
10
|
Schmiester L, Weindl D, Hasenauer J. Efficient gradient-based parameter estimation for dynamic models using qualitative data. Bioinformatics 2021; 37:4493-4500. [PMID: 34260697 PMCID: PMC8652033 DOI: 10.1093/bioinformatics/btab512] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2021] [Revised: 07/02/2021] [Accepted: 07/08/2021] [Indexed: 11/22/2022] Open
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
|
11
|
Borisov I, Metelkin E. Confidence intervals by constrained optimization-An algorithm and software package for practical identifiability analysis in systems biology. PLoS Comput Biol 2020; 16:e1008495. [PMID: 33347435 PMCID: PMC7785248 DOI: 10.1371/journal.pcbi.1008495] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Revised: 01/05/2021] [Accepted: 11/06/2020] [Indexed: 11/30/2022] Open
Abstract
Practical identifiability of Systems Biology models has received a lot of attention in recent scientific research. It addresses the crucial question for models’ predictability: how accurately can the models’ parameters be recovered from available experimental data. The methods based on profile likelihood are among the most reliable methods of practical identification. However, these methods are often computationally demanding or lead to inaccurate estimations of parameters’ confidence intervals. Development of methods, which can accurately produce parameters’ confidence intervals in reasonable computational time, is of utmost importance for Systems Biology and QSP modeling. We propose an algorithm Confidence Intervals by Constraint Optimization (CICO) based on profile likelihood, designed to speed-up confidence intervals estimation and reduce computational cost. The numerical implementation of the algorithm includes settings to control the accuracy of confidence intervals estimates. The algorithm was tested on a number of Systems Biology models, including Taxol treatment model and STAT5 Dimerization model, discussed in the current article. The CICO algorithm is implemented in a software package freely available in Julia (https://github.com/insysbio/LikelihoodProfiler.jl) and Python (https://github.com/insysbio/LikelihoodProfiler.py). Differential equations-based models are widely used in Systems Biology and Quantitative Systems Pharmacology and play a significant role in the discovery of new disease-directed drugs. Complexity of models is a trade off from their employment to crucial fields of biology and medicine. These areas of application require large non-linear models with many unknown parameters. How accurately can the parameters of a model be recovered from experimental data? What is the identifiable subset of parameters? Can the model be reduced or reparameterized to become identifiable? All those questions of identifiability analysis are essential for model’s predictability and reliability. That explains why the topic of identifiability of Systems Biology models has received a lot of attention in recent scientific research. However, existing numerical methods of identifiability analysis are computationally demanding or often lead to inaccurate estimations. Development of methods, which can accurately produce parameters’ confidence intervals in reasonable computational time, is of utmost importance for Systems Biology and QSP modeling. We propose an algorithm and a software package to test identifiability of Systems Biology models, designed to speed-up confidence intervals estimation and reduce computational cost. The software package was tested on a number of Systems Biology models, including Taxol treatment model and STAT5 Dimerization model, discussed in the current article.
Collapse
|
12
|
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
|
13
|
Schälte Y, Hasenauer J. Efficient exact inference for dynamical systems with noisy measurements using sequential approximate Bayesian computation. Bioinformatics 2020; 36:i551-i559. [PMID: 32657404 PMCID: PMC7355286 DOI: 10.1093/bioinformatics/btaa397] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
MOTIVATION Approximate Bayesian computation (ABC) is an increasingly popular method for likelihood-free parameter inference in systems biology and other fields of research, as it allows analyzing complex stochastic models. However, the introduced approximation error is often not clear. It has been shown that ABC actually gives exact inference under the implicit assumption of a measurement noise model. Noise being common in biological systems, it is intriguing to exploit this insight. But this is difficult in practice, as ABC is in general highly computationally demanding. Thus, the question we want to answer here is how to efficiently account for measurement noise in ABC. RESULTS We illustrate exemplarily how ABC yields erroneous parameter estimates when neglecting measurement noise. Then, we discuss practical ways of correctly including the measurement noise in the analysis. We present an efficient adaptive sequential importance sampling-based algorithm applicable to various model types and noise models. We test and compare it on several models, including ordinary and stochastic differential equations, Markov jump processes and stochastically interacting agents, and noise models including normal, Laplace and Poisson noise. We conclude that the proposed algorithm could improve the accuracy of parameter estimates for a broad spectrum of applications. AVAILABILITY AND IMPLEMENTATION The developed algorithms are made publicly available as part of the open-source python toolbox pyABC (https://github.com/icb-dcm/pyabc). SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Yannik Schälte
- Helmholtz Zentrum München, Institute of Computational Biology, Neuherberg 85764, Germany
- Department of Mathematics, Chair of Mathematical Modeling of Biological Systems, Technical University Munich, Garching 85748, Germany
| | - Jan Hasenauer
- Helmholtz Zentrum München, Institute of Computational Biology, Neuherberg 85764, Germany
- Department of Mathematics, Chair of Mathematical Modeling of Biological Systems, Technical University Munich, Garching 85748, Germany
- Research Unit Biomathematics, University of Bonn, Bonn 53113, Germany
| |
Collapse
|
14
|
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
|
15
|
Maurer B, Kollmann S, Pickem J, Hoelbl-Kovacic A, Sexl V. STAT5A and STAT5B-Twins with Different Personalities in Hematopoiesis and Leukemia. Cancers (Basel) 2019; 11:E1726. [PMID: 31690038 PMCID: PMC6895831 DOI: 10.3390/cancers11111726] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2019] [Revised: 10/25/2019] [Accepted: 11/01/2019] [Indexed: 12/14/2022] Open
Abstract
The transcription factors STAT5A and STAT5B have essential roles in survival and proliferation of hematopoietic cells-which have been considered largely redundant. Mutations of upstream kinases, copy number gains, or activating mutations in STAT5A, or more frequently in STAT5B, cause altered hematopoiesis and cancer. Interfering with their activity by pharmacological intervention is an up-and-coming therapeutic avenue. Precision medicine requests detailed knowledge of STAT5A's and STAT5B's individual functions. Recent evidence highlights the privileged role for STAT5B over STAT5A in normal and malignant hematopoiesis. Here, we provide an overview on their individual functions within the hematopoietic system.
Collapse
Affiliation(s)
- Barbara Maurer
- Institute of Pharmacology and Toxicology, University of Veterinary Medicine, 1210 Vienna, Austria.
| | - Sebastian Kollmann
- Institute of Pharmacology and Toxicology, University of Veterinary Medicine, 1210 Vienna, Austria
| | - Judith Pickem
- Institute of Pharmacology and Toxicology, University of Veterinary Medicine, 1210 Vienna, Austria
| | - Andrea Hoelbl-Kovacic
- Institute of Pharmacology and Toxicology, University of Veterinary Medicine, 1210 Vienna, Austria
| | - Veronika Sexl
- Institute of Pharmacology and Toxicology, University of Veterinary Medicine, 1210 Vienna, Austria
| |
Collapse
|
16
|
Mitra ED, Suderman R, Colvin J, Ionkov A, Hu A, Sauro HM, Posner RG, Hlavacek WS. PyBioNetFit and the Biological Property Specification Language. iScience 2019; 19:1012-1036. [PMID: 31522114 PMCID: PMC6744527 DOI: 10.1016/j.isci.2019.08.045] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2019] [Revised: 06/21/2019] [Accepted: 08/22/2019] [Indexed: 02/07/2023] Open
Abstract
In systems biology modeling, important steps include model parameterization, uncertainty quantification, and evaluation of agreement with experimental observations. To help modelers perform these steps, we developed the software PyBioNetFit, which in addition supports checking models against known system properties and solving design problems. PyBioNetFit introduces Biological Property Specification Language (BPSL) for the formal declaration of system properties. BPSL allows qualitative data to be used alone or in combination with quantitative data. PyBioNetFit performs parameterization with parallelized metaheuristic optimization algorithms that work directly with existing model definition standards: BioNetGen Language (BNGL) and Systems Biology Markup Language (SBML). We demonstrate PyBioNetFit's capabilities by solving various example problems, including the challenging problem of parameterizing a 153-parameter model of cell cycle control in yeast based on both quantitative and qualitative data. We demonstrate the model checking and design applications of PyBioNetFit and BPSL by analyzing a model of targeted drug interventions in autophagy signaling.
Collapse
Affiliation(s)
- Eshan D Mitra
- Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, USA
| | - Ryan Suderman
- Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, USA
| | - Joshua Colvin
- Department of Biological Sciences, Northern Arizona University, Flagstaff, AZ, USA
| | - Alexander Ionkov
- Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, USA
| | - Andrew Hu
- Department of Bioengineering, University of Washington, Seattle, WA, USA
| | - Herbert M Sauro
- Department of Bioengineering, University of Washington, Seattle, WA, USA
| | - Richard G Posner
- Department of Biological Sciences, Northern Arizona University, Flagstaff, AZ, USA
| | - William S Hlavacek
- Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, USA.
| |
Collapse
|
17
|
Kreutz C. An easy and efficient approach for testing identifiability. Bioinformatics 2019; 34:1913-1921. [PMID: 29365095 DOI: 10.1093/bioinformatics/bty035] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2017] [Accepted: 01/19/2018] [Indexed: 11/12/2022] Open
Abstract
Motivation The feasibility of uniquely estimating parameters of dynamical systems from observations is a widely discussed aspect of mathematical modelling. Several approaches have been published for analyzing this so-called identifiability of model parameters. However, they are typically computationally demanding, difficult to perform and/or not applicable in many application settings. Results Here, an approach is presented which enables quickly testing of parameter identifiability. Numerical optimization with a penalty in radial direction enforcing displacement of the parameters is used to check whether estimated parameters are unique, or whether the parameters can be altered without loss of agreement with the data indicating non-identifiability. This Identifiability-Test by Radial Penalization (ITRP) can be employed for every model where optimization-based parameter estimation like least-squares or maximum likelihood is feasible and is therefore applicable for all typical systems biology models. The approach is illustrated and tested using 11 ordinary differential equation (ODE) models. Availability and implementation The presented approach can be implemented without great efforts in any modelling framework. It is available within the free Matlab-based modelling toolbox Data2Dynamics. Source code is available at https://github.com/Data2Dynamics. Contact ckreutz@fdm.uni-freiburg.de. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Clemens Kreutz
- Center for Systems Biology (ZBSA), Habsburger Str. 49, University of Freiburg, 79104 Freiburg, Germany.,Center for Data Analysis and Modelling (FDM), Eckerstr. 1, University of Freiburg, 79104 Freiburg, Germany
| |
Collapse
|
18
|
Lucarelli P, Schilling M, Kreutz C, Vlasov A, Boehm ME, Iwamoto N, Steiert B, Lattermann S, Wäsch M, Stepath M, Matter MS, Heikenwälder M, Hoffmann K, Deharde D, Damm G, Seehofer D, Muciek M, Gretz N, Lehmann WD, Timmer J, Klingmüller U. Resolving the Combinatorial Complexity of Smad Protein Complex Formation and Its Link to Gene Expression. Cell Syst 2018; 6:75-89.e11. [DOI: 10.1016/j.cels.2017.11.010] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2016] [Revised: 06/23/2017] [Accepted: 11/14/2017] [Indexed: 12/11/2022]
|
19
|
Adlung L, Kar S, Wagner MC, She B, Chakraborty S, Bao J, Lattermann S, Boerries M, Busch H, Wuchter P, Ho AD, Timmer J, Schilling M, Höfer T, Klingmüller U. Protein abundance of AKT and ERK pathway components governs cell type-specific regulation of proliferation. Mol Syst Biol 2017; 13:904. [PMID: 28123004 PMCID: PMC5293153 DOI: 10.15252/msb.20167258] [Citation(s) in RCA: 57] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Signaling through the AKT and ERK pathways controls cell proliferation. However, the integrated regulation of this multistep process, involving signal processing, cell growth and cell cycle progression, is poorly understood. Here, we study different hematopoietic cell types, in which AKT and ERK signaling is triggered by erythropoietin (Epo). Although these cell types share the molecular network topology for pro‐proliferative Epo signaling, they exhibit distinct proliferative responses. Iterating quantitative experiments and mathematical modeling, we identify two molecular sources for cell type‐specific proliferation. First, cell type‐specific protein abundance patterns cause differential signal flow along the AKT and ERK pathways. Second, downstream regulators of both pathways have differential effects on proliferation, suggesting that protein synthesis is rate‐limiting for faster cycling cells while slower cell cycles are controlled at the G1‐S progression. The integrated mathematical model of Epo‐driven proliferation explains cell type‐specific effects of targeted AKT and ERK inhibitors and faithfully predicts, based on the protein abundance, anti‐proliferative effects of inhibitors in primary human erythroid progenitor cells. Our findings suggest that the effectiveness of targeted cancer therapy might become predictable from protein abundance.
Collapse
Affiliation(s)
- Lorenz Adlung
- Division of Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Sandip Kar
- Division of Theoretical Systems Biology, German Cancer Research Center (DKFZ), Heidelberg, Germany.,BioQuant Center, University of Heidelberg, Heidelberg, Germany.,Department of Chemistry, Indian Institute of Technology, Mumbai, India
| | - Marie-Christine Wagner
- Division of Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Bin She
- Division of Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Sajib Chakraborty
- Division of Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Jie Bao
- Systems Biology of the Cellular Microenvironment Group, IMMZ, ALU, Freiburg, Germany
| | - Susen Lattermann
- Division of Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Melanie Boerries
- Systems Biology of the Cellular Microenvironment Group, IMMZ, ALU, Freiburg, Germany.,German Cancer Consortium (DKTK), Freiburg, Germany.,German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Hauke Busch
- Systems Biology of the Cellular Microenvironment Group, IMMZ, ALU, Freiburg, Germany.,German Cancer Consortium (DKTK), Freiburg, Germany.,German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Patrick Wuchter
- Department of Medicine V, University of Heidelberg, Heidelberg, Germany.,Institute for Transfusion Medicine and Immunology, University of Heidelberg, Mannheim, Germany
| | - Anthony D Ho
- Department of Medicine V, University of Heidelberg, Heidelberg, Germany
| | - Jens Timmer
- Center for Biological Signaling Studies (BIOSS), Institute of Physics, University of Freiburg, Freiburg, Germany
| | - Marcel Schilling
- Division of Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Thomas Höfer
- Division of Theoretical Systems Biology, German Cancer Research Center (DKFZ), Heidelberg, Germany .,BioQuant Center, University of Heidelberg, Heidelberg, Germany
| | - Ursula Klingmüller
- Division of Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), Heidelberg, Germany .,Translational Lung Research Center (TLRC), Member of the German Center for Lung Research (DZL), Heidelberg, Germany
| |
Collapse
|
20
|
Lehmann WD. A timeline of stable isotopes and mass spectrometry in the life sciences. MASS SPECTROMETRY REVIEWS 2017; 36:58-85. [PMID: 26919394 DOI: 10.1002/mas.21497] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/07/2015] [Accepted: 01/21/2016] [Indexed: 06/05/2023]
Abstract
This review retraces the role of stable isotopes and mass spectrometry in the life sciences. The timeline is divided into four segments covering the years 1920-1950, 1950-1980, 1980-2000, and 2000 until today. For each period methodic progress and typical applications are discussed. Application of stable isotopes is driven by improvements of mass spectrometry, chromatography, and related fields in sensitivity, mass accuracy, structural specificity, complex sample handling ability, data output, and data evaluation. We currently experience the vision of omics-type analyses, that is, the comprehensive identification and quantification of a complete compound class within one or a few analytical runs. This development is driven by stable isotopes without competition by radioisotopes. In metabolic studies as classic field of isotopic tracer experiments, stable isotopes and radioisotopes were competing solutions, with stable isotopes as the long-term junior partner. Since the 1990s the number of metabolic studies with radioisotopes decreases, whereas stable isotope studies retain their slow but stable upward tendency. Unique fields of stable isotopes are metabolic tests in newborns, metabolic experiments in healthy controls, newborn screening for inborn errors, quantification of drugs and drug metabolites in doping control, natural isotope fractionation in geology, ecology, food authentication, or doping control, and more recently the field of quantitative omics-type analyses. There, cells or whole organisms are systematically labeled with stable isotopes to study proteomic differences or specific responses to stimuli or genetic manipulation. The duo of stable isotopes and mass spectrometry will probably continue to grow in the life sciences, since it delivers reference-quality quantitative data with molecular specificity, often combined with informative isotope effects. © 2016 Wiley Periodicals, Inc. Mass Spec Rev 36:58-85, 2017.
Collapse
Affiliation(s)
- Wolf D Lehmann
- German Cancer Research Center (DKFZ), D-69120 Heidelberg, Germany
| |
Collapse
|
21
|
Merkle R, Steiert B, Salopiata F, Depner S, Raue A, Iwamoto N, Schelker M, Hass H, Wäsch M, Böhm ME, Mücke O, Lipka DB, Plass C, Lehmann WD, Kreutz C, Timmer J, Schilling M, Klingmüller U. Identification of Cell Type-Specific Differences in Erythropoietin Receptor Signaling in Primary Erythroid and Lung Cancer Cells. PLoS Comput Biol 2016; 12:e1005049. [PMID: 27494133 PMCID: PMC4975441 DOI: 10.1371/journal.pcbi.1005049] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2016] [Accepted: 07/05/2016] [Indexed: 01/23/2023] Open
Abstract
Lung cancer, with its most prevalent form non-small-cell lung carcinoma (NSCLC), is one of the leading causes of cancer-related deaths worldwide, and is commonly treated with chemotherapeutic drugs such as cisplatin. Lung cancer patients frequently suffer from chemotherapy-induced anemia, which can be treated with erythropoietin (EPO). However, studies have indicated that EPO not only promotes erythropoiesis in hematopoietic cells, but may also enhance survival of NSCLC cells. Here, we verified that the NSCLC cell line H838 expresses functional erythropoietin receptors (EPOR) and that treatment with EPO reduces cisplatin-induced apoptosis. To pinpoint differences in EPO-induced survival signaling in erythroid progenitor cells (CFU-E, colony forming unit-erythroid) and H838 cells, we combined mathematical modeling with a method for feature selection, the L1 regularization. Utilizing an example model and simulated data, we demonstrated that this approach enables the accurate identification and quantification of cell type-specific parameters. We applied our strategy to quantitative time-resolved data of EPO-induced JAK/STAT signaling generated by quantitative immunoblotting, mass spectrometry and quantitative real-time PCR (qRT-PCR) in CFU-E and H838 cells as well as H838 cells overexpressing human EPOR (H838-HA-hEPOR). The established parsimonious mathematical model was able to simultaneously describe the data sets of CFU-E, H838 and H838-HA-hEPOR cells. Seven cell type-specific parameters were identified that included for example parameters for nuclear translocation of STAT5 and target gene induction. Cell type-specific differences in target gene induction were experimentally validated by qRT-PCR experiments. The systematic identification of pathway differences and sensitivities of EPOR signaling in CFU-E and H838 cells revealed potential targets for intervention to selectively inhibit EPO-induced signaling in the tumor cells but leave the responses in erythroid progenitor cells unaffected. Thus, the proposed modeling strategy can be employed as a general procedure to identify cell type-specific parameters and to recommend treatment strategies for the selective targeting of specific cell types. A major challenge in the development of therapeutic interventions is the selective inhibition of a signal transduction pathway in one cell type such as a cancer cell leaving the other cell type such as a healthy cell as unaffected as possible. Here, we propose a new approach that combines mathematical modeling based on quantitative experimental data with statistical methods. We demonstrate based on simulated data that our approach can determine which parameters are the same and which parameters differ in two exemplary cell types. We compare a lung cancer cell line to the precursor cells of red blood cells. We show that the same signal transduction network induced by erythropoietin (EPO), a hormone that is frequently employed to treat anemia in cancer patients, regulates survival of both cell types. Based on our experimental data in combination with our computational approach, we identify seven cell type-specific differences in this signaling pathway. Our strategy allows predicting therapeutic targets that could be inhibited to interfere with survival of lung cancer cells while leaving production of red blood cells unaffected.
Collapse
Affiliation(s)
- Ruth Merkle
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), INF 280, Heidelberg, Germany
- Translational Lung Research Center (TLRC), German Center for Lung Research (DZL), Heidelberg, Germany
| | - Bernhard Steiert
- Institute of Physics, University of Freiburg, Germany & BIOSS Centre for Biological Signalling Studies, University of Freiburg, Germany
| | - Florian Salopiata
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), INF 280, Heidelberg, Germany
- Translational Lung Research Center (TLRC), German Center for Lung Research (DZL), Heidelberg, Germany
| | - Sofia Depner
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), INF 280, Heidelberg, Germany
- Translational Lung Research Center (TLRC), German Center for Lung Research (DZL), Heidelberg, Germany
| | - Andreas Raue
- Institute of Physics, University of Freiburg, Germany & BIOSS Centre for Biological Signalling Studies, University of Freiburg, Germany
| | - Nao Iwamoto
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), INF 280, Heidelberg, Germany
| | - Max Schelker
- Institute of Physics, University of Freiburg, Germany & BIOSS Centre for Biological Signalling Studies, University of Freiburg, Germany
| | - Helge Hass
- Institute of Physics, University of Freiburg, Germany & BIOSS Centre for Biological Signalling Studies, University of Freiburg, Germany
| | - Marvin Wäsch
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), INF 280, Heidelberg, Germany
- Translational Lung Research Center (TLRC), German Center for Lung Research (DZL), Heidelberg, Germany
| | - Martin E. Böhm
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), INF 280, Heidelberg, Germany
| | - Oliver Mücke
- Division Epigenomics and Cancer Risk Factors, German Cancer Research Center (DKFZ), INF 280, Heidelberg, Germany
| | - Daniel B. Lipka
- Regulation of Cellular Differentiation Group, Division Epigenomics and Cancer Risk Factors, German Cancer Research Center (DKFZ), INF 280, Heidelberg, Germany
| | - Christoph Plass
- Division Epigenomics and Cancer Risk Factors, German Cancer Research Center (DKFZ), INF 280, Heidelberg, Germany
| | - Wolf D. Lehmann
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), INF 280, Heidelberg, Germany
| | - Clemens Kreutz
- Institute of Physics, University of Freiburg, Germany & BIOSS Centre for Biological Signalling Studies, University of Freiburg, Germany
| | - Jens Timmer
- Institute of Physics, University of Freiburg, Germany & BIOSS Centre for Biological Signalling Studies, University of Freiburg, Germany
- * E-mail: (JT); (MS); (UK)
| | - Marcel Schilling
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), INF 280, Heidelberg, Germany
- * E-mail: (JT); (MS); (UK)
| | - Ursula Klingmüller
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), INF 280, Heidelberg, Germany
- Translational Lung Research Center (TLRC), German Center for Lung Research (DZL), Heidelberg, Germany
- * E-mail: (JT); (MS); (UK)
| |
Collapse
|
22
|
Verma R, Green JM, Schatz PJ, Wojchowski DM. A dimeric peptide with erythropoiesis-stimulating activity uniquely affects erythropoietin receptor ligation and cell surface expression. Exp Hematol 2016; 44:765-769.e1. [PMID: 27174804 DOI: 10.1016/j.exphem.2016.04.015] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2016] [Revised: 04/14/2016] [Accepted: 04/20/2016] [Indexed: 10/21/2022]
Abstract
Erythropoiesis-stimulating agents (ESAs) that exert long-acting antianemia effects have been developed recently, but their mechanisms are poorly understood. Analyses reveal unique erythropoietin receptor (EPOR)-binding properties for one such ESA, the synthetic EPOR agonist peginesatide. Compared with recombinant human EPO and darbepoietin, peginesatide exhibited a slow on rate, but sustained EPOR residency and resistant displacement. In EPO-dependent human erythroid progenitor UT7epo cells, culture in peginesatide unexpectedly upmodulated endogenous cell surface EPOR levels with parallel increases in full-length EPOR-68K levels. These unique properties are suggested to contribute to the durable activity of this (and perhaps additional) dimeric peptide hematopoietic growth factor receptor agonist.
Collapse
Affiliation(s)
- Rakesh Verma
- Molecular Medicine Division, Maine Medical Center Research Institute, Scarborough, ME
| | | | | | - Don M Wojchowski
- Molecular Medicine Division, Maine Medical Center Research Institute, Scarborough, ME; Department of Medicine, Tufts University Medical Center, Boston, MA.
| |
Collapse
|
23
|
Kreutz C. New Concepts for Evaluating the Performance of Computational Methods * *The author acknowledge financial support by the by the German Ministry of Education and Research (BMBF) via e:Bio Grant No. 031L0080. ACTA ACUST UNITED AC 2016. [DOI: 10.1016/j.ifacol.2016.12.104] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|