1
|
Kumar N, Mim MS, Dowling A, Zartman JJ. Reverse engineering morphogenesis through Bayesian optimization of physics-based models. NPJ Syst Biol Appl 2024; 10:49. [PMID: 38714708 PMCID: PMC11076624 DOI: 10.1038/s41540-024-00375-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2023] [Accepted: 04/17/2024] [Indexed: 05/10/2024] Open
Abstract
Morphogenetic programs coordinate cell signaling and mechanical interactions to shape organs. In systems and synthetic biology, a key challenge is determining optimal cellular interactions for predicting organ shape, size, and function. Physics-based models defining the subcellular force distribution facilitate this, but it is challenging to calibrate parameters in these models from data. To solve this inverse problem, we created a Bayesian optimization framework to determine the optimal cellular force distribution such that the predicted organ shapes match the experimentally observed organ shapes. This integrative framework employs Gaussian Process Regression, a non-parametric kernel-based probabilistic machine learning modeling paradigm, to learn the mapping functions relating to the morphogenetic programs that maintain the final organ shape. We calibrated and tested the method on Drosophila wing imaginal discs to study mechanisms that regulate epithelial processes ranging from development to cancer. The parameter estimation framework successfully infers the underlying changes in core parameters needed to match simulation data with imaging data of wing discs perturbed with collagenase. The computational pipeline identifies distinct parameter sets mimicking wild-type shapes. It enables a global sensitivity analysis to support the regulation of actomyosin contractility and basal ECM stiffness to generate and maintain the curved shape of the wing imaginal disc. The optimization framework, combined with experimental imaging, identified that Piezo, a mechanosensitive ion channel, impacts fold formation by regulating the apical-basal balance of actomyosin contractility and elasticity of ECM. This workflow is extensible toward reverse-engineering morphogenesis across organ systems and for real-time control of complex multicellular systems.
Collapse
Affiliation(s)
- Nilay Kumar
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, IN, 46556, USA
| | - Mayesha Sahir Mim
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, IN, 46556, USA
- Bioengineering Graduate Program, University of Notre Dame, Notre Dame, IN, 46556, USA
| | - Alexander Dowling
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, IN, 46556, USA
| | - Jeremiah J Zartman
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, IN, 46556, USA.
- Bioengineering Graduate Program, University of Notre Dame, Notre Dame, IN, 46556, USA.
| |
Collapse
|
2
|
Cain JY, Evarts JI, Yu JS, Bagheri N. Incorporating temporal information during feature engineering bolsters emulation of spatio-temporal emergence. BIOINFORMATICS (OXFORD, ENGLAND) 2024; 40:btae131. [PMID: 38444088 PMCID: PMC10957516 DOI: 10.1093/bioinformatics/btae131] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/26/2023] [Revised: 02/08/2024] [Accepted: 03/01/2024] [Indexed: 03/07/2024]
Abstract
MOTIVATION Emergent biological dynamics derive from the evolution of lower-level spatial and temporal processes. A long-standing challenge for scientists and engineers is identifying simple low-level rules that give rise to complex higher-level dynamics. High-resolution biological data acquisition enables this identification and has evolved at a rapid pace for both experimental and computational approaches. Simultaneously harnessing the resolution and managing the expense of emerging technologies-e.g. live cell imaging, scRNAseq, agent-based models-requires a deeper understanding of how spatial and temporal axes impact biological systems. Effective emulation is a promising solution to manage the expense of increasingly complex high-resolution computational models. In this research, we focus on the emulation of a tumor microenvironment agent-based model to examine the relationship between spatial and temporal environment features, and emergent tumor properties. RESULTS Despite significant feature engineering, we find limited predictive capacity of tumor properties from initial system representations. However, incorporating temporal information derived from intermediate simulation states dramatically improves the predictive performance of machine learning models. We train a deep-learning emulator on intermediate simulation states and observe promising enhancements over emulators trained solely on initial conditions. Our results underscore the importance of incorporating temporal information in the evaluation of spatio-temporal emergent behavior. Nevertheless, the emulators exhibit inconsistent performance, suggesting that the underlying model characterizes unique cell populations dynamics that are not easily replaced. AVAILABILITY AND IMPLEMENTATION All source codes for the agent-based model, emulation, and analyses are publicly available at the corresponding DOIs: 10.5281/zenodo.10622155, 10.5281/zenodo.10611675, 10.5281/zenodo.10621244, respectively.
Collapse
Affiliation(s)
- Jason Y Cain
- Department of Chemical Engineering, University of Washington, Seattle, WA 98195, United States
| | - Jacob I Evarts
- Department of Biology, University of Washington, Seattle, WA 98195, United States
| | - Jessica S Yu
- Department of Biology, University of Washington, Seattle, WA 98195, United States
| | - Neda Bagheri
- Department of Chemical Engineering, University of Washington, Seattle, WA 98195, United States
- Department of Biology, University of Washington, Seattle, WA 98195, United States
| |
Collapse
|
3
|
Kumar N, Dowling A, Zartman J. Reverse engineering morphogenesis through Bayesian optimization of physics-based models. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.21.553928. [PMID: 37662294 PMCID: PMC10473585 DOI: 10.1101/2023.08.21.553928] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/05/2023]
Abstract
Morphogenetic programs direct the cell signaling and nonlinear mechanical interactions between multiple cell types and tissue layers to define organ shape and size. A key challenge for systems and synthetic biology is determining optimal combinations of intra- and inter-cellular interactions to predict an organ's shape, size, and function. Physics-based mechanistic models that define the subcellular force distribution facilitate this, but it is extremely challenging to calibrate parameters in these models from data. To solve this inverse problem, we created a Bayesian optimization framework to determine the optimal cellular force distribution such that the predicted organ shapes match the desired organ shapes observed within the experimental imaging data. This integrative framework employs Gaussian Process Regression (GPR), a non-parametric kernel-based probabilistic machine learning modeling paradigm, to learn the mapping functions relating to the morphogenetic programs that generate and maintain the final organ shape. We calibrated and tested the method on cross-sections of Drosophila wing imaginal discs, a highly informative model organ system, to study mechanisms that regulate epithelial processes that range from development to cancer. As a specific test case, the parameter estimation framework successfully infers the underlying changes in core parameters needed to match simulation data with time series imaging data of wing discs perturbed with collagenase. Unexpectedly, the framework also identifies multiple distinct parameter sets that generate shapes similar to wild-type organ shapes. This platform enables an efficient, global sensitivity analysis to support the necessity of both actomyosin contractility and basal ECM stiffness to generate and maintain the curved shape of the wing imaginal disc. The optimization framework, combined with fixed tissue imaging, identified that Piezo, a mechanosensitive ion channel, impacts fold formation by regulating the apical-basal balance of actomyosin contractility and elasticity of ECM. This framework is extensible toward reverse-engineering the morphogenesis of any organ system and can be utilized in real-time control of complex multicellular systems.
Collapse
|
4
|
Kondo Y, Achouri NL, Falou HA, Atar L, Aumann T, Baba H, Boretzky K, Caesar C, Calvet D, Chae H, Chiga N, Corsi A, Delaunay F, Delbart A, Deshayes Q, Dombrádi Z, Douma CA, Ekström A, Elekes Z, Forssén C, Gašparić I, Gheller JM, Gibelin J, Gillibert A, Hagen G, Harakeh MN, Hirayama A, Hoffman CR, Holl M, Horvat A, Horváth Á, Hwang JW, Isobe T, Jiang WG, Kahlbow J, Kalantar-Nayestanaki N, Kawase S, Kim S, Kisamori K, Kobayashi T, Körper D, Koyama S, Kuti I, Lapoux V, Lindberg S, Marqués FM, Masuoka S, Mayer J, Miki K, Murakami T, Najafi M, Nakamura T, Nakano K, Nakatsuka N, Nilsson T, Obertelli A, Ogata K, de Oliveira Santos F, Orr NA, Otsu H, Otsuka T, Ozaki T, Panin V, Papenbrock T, Paschalis S, Revel A, Rossi D, Saito AT, Saito TY, Sasano M, Sato H, Satou Y, Scheit H, Schindler F, Schrock P, Shikata M, Shimizu N, Shimizu Y, Simon H, Sohler D, Sorlin O, Stuhl L, Sun ZH, Takeuchi S, Tanaka M, Thoennessen M, Törnqvist H, Togano Y, Tomai T, Tscheuschner J, Tsubota J, Tsunoda N, Uesaka T, Utsuno Y, Vernon I, Wang H, Yang Z, Yasuda M, Yoneda K, Yoshida S. First observation of 28O. Nature 2023; 620:965-970. [PMID: 37648757 PMCID: PMC10630140 DOI: 10.1038/s41586-023-06352-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2022] [Accepted: 06/21/2023] [Indexed: 09/01/2023]
Abstract
Subjecting a physical system to extreme conditions is one of the means often used to obtain a better understanding and deeper insight into its organization and structure. In the case of the atomic nucleus, one such approach is to investigate isotopes that have very different neutron-to-proton (N/Z) ratios than in stable nuclei. Light, neutron-rich isotopes exhibit the most asymmetric N/Z ratios and those lying beyond the limits of binding, which undergo spontaneous neutron emission and exist only as very short-lived resonances (about 10-21 s), provide the most stringent tests of modern nuclear-structure theories. Here we report on the first observation of 28O and 27O through their decay into 24O and four and three neutrons, respectively. The 28O nucleus is of particular interest as, with the Z = 8 and N = 20 magic numbers1,2, it is expected in the standard shell-model picture of nuclear structure to be one of a relatively small number of so-called 'doubly magic' nuclei. Both 27O and 28O were found to exist as narrow, low-lying resonances and their decay energies are compared here to the results of sophisticated theoretical modelling, including a large-scale shell-model calculation and a newly developed statistical approach. In both cases, the underlying nuclear interactions were derived from effective field theories of quantum chromodynamics. Finally, it is shown that the cross-section for the production of 28O from a 29F beam is consistent with it not exhibiting a closed N = 20 shell structure.
Collapse
Affiliation(s)
- Y Kondo
- Department of Physics, Tokyo Institute of Technology, Tokyo, Japan.
- RIKEN Nishina Center, Saitama, Japan.
| | - N L Achouri
- LPC Caen UMR6534, Université de Caen Normandie, ENSICAEN, CNRS/IN2P3, Caen, France
| | - H Al Falou
- Lebanese University, Beirut, Lebanon
- Lebanese-French University of Technology and Applied Sciences, Deddeh, Lebanon
| | - L Atar
- Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, Germany
| | - T Aumann
- Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, Germany
- GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany
- Helmholtz Research Academy Hesse for FAIR, Darmstadt, Germany
| | - H Baba
- RIKEN Nishina Center, Saitama, Japan
| | - K Boretzky
- GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany
| | - C Caesar
- Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, Germany
- GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany
| | - D Calvet
- Irfu, CEA, Université Paris-Saclay, Gif-sur-Yvette, France
| | - H Chae
- Institute for Basic Science, Daejeon, Republic of Korea
| | - N Chiga
- RIKEN Nishina Center, Saitama, Japan
| | - A Corsi
- Irfu, CEA, Université Paris-Saclay, Gif-sur-Yvette, France
| | - F Delaunay
- LPC Caen UMR6534, Université de Caen Normandie, ENSICAEN, CNRS/IN2P3, Caen, France
| | - A Delbart
- Irfu, CEA, Université Paris-Saclay, Gif-sur-Yvette, France
| | - Q Deshayes
- LPC Caen UMR6534, Université de Caen Normandie, ENSICAEN, CNRS/IN2P3, Caen, France
| | | | - C A Douma
- ESRIG, University of Groningen, Groningen, The Netherlands
| | - A Ekström
- Institutionen för Fysik, Chalmers Tekniska Högskola, Göteborg, Sweden
| | | | - C Forssén
- Institutionen för Fysik, Chalmers Tekniska Högskola, Göteborg, Sweden
| | - I Gašparić
- RIKEN Nishina Center, Saitama, Japan
- Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, Germany
- Ruđer Bošković Institute, Zagreb, Croatia
| | - J-M Gheller
- Irfu, CEA, Université Paris-Saclay, Gif-sur-Yvette, France
| | - J Gibelin
- LPC Caen UMR6534, Université de Caen Normandie, ENSICAEN, CNRS/IN2P3, Caen, France
| | - A Gillibert
- Irfu, CEA, Université Paris-Saclay, Gif-sur-Yvette, France
| | - G Hagen
- Physics Division, Oak Ridge National Laboratory, Oak Ridge, TN, USA
- Department of Physics and Astronomy, University of Tennessee, Knoxville, TN, USA
| | - M N Harakeh
- GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany
- ESRIG, University of Groningen, Groningen, The Netherlands
| | - A Hirayama
- Department of Physics, Tokyo Institute of Technology, Tokyo, Japan
| | - C R Hoffman
- Physics Division, Argonne National Laboratory, Argonne, IL, USA
| | - M Holl
- Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, Germany
- GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany
| | - A Horvat
- GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany
| | - Á Horváth
- Eötvös Loránd University, Budapest, Hungary
| | - J W Hwang
- Center for Exotic Nuclear Studies, Institute for Basic Science, Daejeon, Republic of Korea
- Department of Physics and Astronomy, Seoul National University, Seoul, Republic of Korea
| | - T Isobe
- RIKEN Nishina Center, Saitama, Japan
| | - W G Jiang
- Institutionen för Fysik, Chalmers Tekniska Högskola, Göteborg, Sweden
| | - J Kahlbow
- RIKEN Nishina Center, Saitama, Japan
- Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, Germany
| | | | - S Kawase
- Department of Advanced Energy Engineering Science, Kyushu University, Fukuoka, Japan
| | - S Kim
- Center for Exotic Nuclear Studies, Institute for Basic Science, Daejeon, Republic of Korea
- Department of Physics and Astronomy, Seoul National University, Seoul, Republic of Korea
| | | | - T Kobayashi
- Department of Physics, Tohoku University, Miyagi, Japan
| | - D Körper
- GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany
| | - S Koyama
- Department of Physics, The University of Tokyo, Tokyo, Japan
| | - I Kuti
- Atomki, Debrecen, Hungary
| | - V Lapoux
- Irfu, CEA, Université Paris-Saclay, Gif-sur-Yvette, France
| | - S Lindberg
- Institutionen för Fysik, Chalmers Tekniska Högskola, Göteborg, Sweden
| | - F M Marqués
- LPC Caen UMR6534, Université de Caen Normandie, ENSICAEN, CNRS/IN2P3, Caen, France
| | - S Masuoka
- Center for Nuclear Study, The University of Tokyo, Saitama, Japan
| | - J Mayer
- Institut für Kernphysik, Universität zu Köln, Köln, Germany
| | - K Miki
- Department of Physics, Tohoku University, Miyagi, Japan
| | - T Murakami
- Department of Physics, Kyoto University, Kyoto, Japan
| | - M Najafi
- ESRIG, University of Groningen, Groningen, The Netherlands
| | - T Nakamura
- Department of Physics, Tokyo Institute of Technology, Tokyo, Japan
- RIKEN Nishina Center, Saitama, Japan
| | - K Nakano
- Department of Advanced Energy Engineering Science, Kyushu University, Fukuoka, Japan
| | - N Nakatsuka
- Department of Physics, Kyoto University, Kyoto, Japan
| | - T Nilsson
- Institutionen för Fysik, Chalmers Tekniska Högskola, Göteborg, Sweden
| | - A Obertelli
- Irfu, CEA, Université Paris-Saclay, Gif-sur-Yvette, France
| | - K Ogata
- Department of Physics, Kyushu University, Fukuoka, Japan
- Research Center for Nuclear Physics, Osaka University, Osaka, Japan
- Department of Physics, Osaka City University, Osaka, Japan
| | - F de Oliveira Santos
- Grand Accélérateur National d'Ions Lourds (GANIL), CEA/DRF-CNRS/IN2P3, Caen, France
| | - N A Orr
- LPC Caen UMR6534, Université de Caen Normandie, ENSICAEN, CNRS/IN2P3, Caen, France
| | - H Otsu
- RIKEN Nishina Center, Saitama, Japan
| | - T Otsuka
- RIKEN Nishina Center, Saitama, Japan
- Department of Physics, The University of Tokyo, Tokyo, Japan
| | - T Ozaki
- Department of Physics, Tokyo Institute of Technology, Tokyo, Japan
| | - V Panin
- RIKEN Nishina Center, Saitama, Japan
| | - T Papenbrock
- Physics Division, Oak Ridge National Laboratory, Oak Ridge, TN, USA
- Department of Physics and Astronomy, University of Tennessee, Knoxville, TN, USA
| | - S Paschalis
- Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, Germany
| | - A Revel
- LPC Caen UMR6534, Université de Caen Normandie, ENSICAEN, CNRS/IN2P3, Caen, France
- Grand Accélérateur National d'Ions Lourds (GANIL), CEA/DRF-CNRS/IN2P3, Caen, France
| | - D Rossi
- Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, Germany
| | - A T Saito
- Department of Physics, Tokyo Institute of Technology, Tokyo, Japan
| | - T Y Saito
- Department of Physics, The University of Tokyo, Tokyo, Japan
| | - M Sasano
- RIKEN Nishina Center, Saitama, Japan
| | - H Sato
- RIKEN Nishina Center, Saitama, Japan
| | - Y Satou
- Department of Physics and Astronomy, Seoul National University, Seoul, Republic of Korea
| | - H Scheit
- Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, Germany
| | - F Schindler
- Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, Germany
| | - P Schrock
- Center for Nuclear Study, The University of Tokyo, Saitama, Japan
| | - M Shikata
- Department of Physics, Tokyo Institute of Technology, Tokyo, Japan
| | - N Shimizu
- Center for Computational Sciences, University of Tsukuba, Ibaraki, Japan
| | - Y Shimizu
- RIKEN Nishina Center, Saitama, Japan
| | - H Simon
- GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany
| | | | - O Sorlin
- Grand Accélérateur National d'Ions Lourds (GANIL), CEA/DRF-CNRS/IN2P3, Caen, France
| | - L Stuhl
- RIKEN Nishina Center, Saitama, Japan
- Center for Exotic Nuclear Studies, Institute for Basic Science, Daejeon, Republic of Korea
| | - Z H Sun
- Physics Division, Oak Ridge National Laboratory, Oak Ridge, TN, USA
- Department of Physics and Astronomy, University of Tennessee, Knoxville, TN, USA
| | - S Takeuchi
- Department of Physics, Tokyo Institute of Technology, Tokyo, Japan
| | - M Tanaka
- Department of Physics, Osaka University, Osaka, Japan
| | - M Thoennessen
- Facility for Rare Isotope Beams, Michigan State University, East Lansing, MI, USA
| | - H Törnqvist
- Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, Germany
- GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany
| | - Y Togano
- Department of Physics, Tokyo Institute of Technology, Tokyo, Japan
- Department of Physics, Rikkyo University, Tokyo, Japan
| | - T Tomai
- Department of Physics, Tokyo Institute of Technology, Tokyo, Japan
| | - J Tscheuschner
- Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, Germany
| | - J Tsubota
- Department of Physics, Tokyo Institute of Technology, Tokyo, Japan
| | - N Tsunoda
- Center for Nuclear Study, The University of Tokyo, Saitama, Japan
| | - T Uesaka
- RIKEN Nishina Center, Saitama, Japan
| | - Y Utsuno
- Advanced Science Research Center, Japan Atomic Energy Agency, Ibaraki, Japan
| | - I Vernon
- Department of Mathematical Sciences, Durham University, Durham, UK
| | - H Wang
- RIKEN Nishina Center, Saitama, Japan
| | - Z Yang
- RIKEN Nishina Center, Saitama, Japan
| | - M Yasuda
- Department of Physics, Tokyo Institute of Technology, Tokyo, Japan
| | - K Yoneda
- RIKEN Nishina Center, Saitama, Japan
| | - S Yoshida
- Liberal and General Education Center, Institute for Promotion of Higher Academic Education, Utsunomiya University, Tochigi, Japan
| |
Collapse
|
5
|
Dokuchaev A, Kursanov A, Balakina-Vikulova NA, Katsnelson LB, Solovyova O. The importance of mechanical conditions in the testing of excitation abnormalities in a population of electro-mechanical models of human ventricular cardiomyocytes. Front Physiol 2023; 14:1187956. [PMID: 37362439 PMCID: PMC10285544 DOI: 10.3389/fphys.2023.1187956] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2023] [Accepted: 05/25/2023] [Indexed: 06/28/2023] Open
Abstract
Background: Populations of in silico electrophysiological models of human cardiomyocytes represent natural variability in cell activity and are thoroughly calibrated and validated using experimental data from the human heart. The models have been shown to predict the effects of drugs and their pro-arrhythmic risks. However, excitation and contraction are known to be tightly coupled in the myocardium, with mechanical loads and stretching affecting both mechanics and excitation through mechanisms of mechano-calcium-electrical feedback. However, these couplings are not currently a focus of populations of cell models. Aim: We investigated the role of cardiomyocyte mechanical activity under different mechanical conditions in the generation, calibration, and validation of a population of electro-mechanical models of human cardiomyocytes. Methods: To generate a population, we assumed 11 input parameters of ionic currents and calcium dynamics in our recently developed TP + M model as varying within a wide range. A History matching algorithm was used to generate a non-implausible parameter space by calibrating the action potential and calcium transient biomarkers against experimental data and rejecting models with excitation abnormalities. The population was further calibrated using experimental data on human myocardial force characteristics and mechanical tests involving variations in preload and afterload. Models that passed the mechanical tests were validated with additional experimental data, including the effects of drugs with high or low pro-arrhythmic risk. Results: More than 10% of the models calibrated on electrophysiological data failed mechanical tests and were rejected from the population due to excitation abnormalities at reduced preload or afterload for cell contraction. The final population of accepted models yielded action potential, calcium transient, and force/shortening outputs consistent with experimental data. In agreement with experimental and clinical data, the models demonstrated a high frequency of excitation abnormalities in simulations of Dofetilide action on the ionic currents, in contrast to Verapamil. However, Verapamil showed a high frequency of failed contractions at high concentrations. Conclusion: Our results highlight the importance of considering mechanoelectric coupling in silico cardiomyocyte models. Mechanical tests allow a more thorough assessment of the effects of interventions on cardiac function, including drug testing.
Collapse
Affiliation(s)
- Arsenii Dokuchaev
- Laboratory of Mathematical Physiology, Institute of Immunology and Physiology, Ural Branch of Russian Academy of Sciences, Ekaterinburg, Russia
| | - Alexander Kursanov
- Laboratory of Mathematical Physiology, Institute of Immunology and Physiology, Ural Branch of Russian Academy of Sciences, Ekaterinburg, Russia
- Laboratory of Mathematical Modeling in Physiology and Medicine Based on Supercomputers, Ural Federal University, Ekaterinburg, Russia
| | - Nathalie A. Balakina-Vikulova
- Laboratory of Mathematical Physiology, Institute of Immunology and Physiology, Ural Branch of Russian Academy of Sciences, Ekaterinburg, Russia
- Laboratory of Mathematical Modeling in Physiology and Medicine Based on Supercomputers, Ural Federal University, Ekaterinburg, Russia
| | - Leonid B. Katsnelson
- Laboratory of Mathematical Physiology, Institute of Immunology and Physiology, Ural Branch of Russian Academy of Sciences, Ekaterinburg, Russia
- Laboratory of Mathematical Modeling in Physiology and Medicine Based on Supercomputers, Ural Federal University, Ekaterinburg, Russia
| | - Olga Solovyova
- Laboratory of Mathematical Physiology, Institute of Immunology and Physiology, Ural Branch of Russian Academy of Sciences, Ekaterinburg, Russia
- Laboratory of Mathematical Modeling in Physiology and Medicine Based on Supercomputers, Ural Federal University, Ekaterinburg, Russia
| |
Collapse
|
6
|
Scarponi D, Iskauskas A, Clark RA, Vernon I, McKinley TJ, Goldstein M, Mukandavire C, Deol A, Weerasuriya C, Bakker R, White RG, McCreesh N. Demonstrating multi-country calibration of a tuberculosis model using new history matching and emulation package - hmer. Epidemics 2023; 43:100678. [PMID: 36913805 DOI: 10.1016/j.epidem.2023.100678] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2022] [Revised: 02/23/2023] [Accepted: 03/06/2023] [Indexed: 03/09/2023] Open
Abstract
Infectious disease models are widely used by epidemiologists to improve the understanding of transmission dynamics and disease natural history, and to predict the possible effects of interventions. As the complexity of such models increases, however, it becomes increasingly challenging to robustly calibrate them to empirical data. History matching with emulation is a calibration method that has been successfully applied to such models, but has not been widely used in epidemiology partly due to the lack of available software. To address this issue, we developed a new, user-friendly R package hmer to simply and efficiently perform history matching with emulation. In this paper, we demonstrate the first use of hmer for calibrating a complex deterministic model for the country-level implementation of tuberculosis vaccines to 115 low- and middle-income countries. The model was fit to 9-13 target measures, by varying 19-22 input parameters. Overall, 105 countries were successfully calibrated. Among the remaining countries, hmer visualisation tools, combined with derivative emulation methods, provided strong evidence that the models were misspecified and could not be calibrated to the target ranges. This work shows that hmer can be used to simply and rapidly calibrate a complex model to data from over 100 countries, making it a useful addition to the epidemiologist's calibration tool-kit.
Collapse
Affiliation(s)
- Danny Scarponi
- Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, London School of Hygiene and Tropical Medicine, UK.
| | | | - Rebecca A Clark
- Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, London School of Hygiene and Tropical Medicine, UK
| | - Ian Vernon
- Department of Mathematical Sciences, Durham University, UK
| | | | | | - Christinah Mukandavire
- Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, London School of Hygiene and Tropical Medicine, UK
| | - Arminder Deol
- Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, London School of Hygiene and Tropical Medicine, UK
| | - Chathika Weerasuriya
- Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, London School of Hygiene and Tropical Medicine, UK
| | - Roel Bakker
- Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, London School of Hygiene and Tropical Medicine, UK
| | - Richard G White
- Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, London School of Hygiene and Tropical Medicine, UK
| | - Nicky McCreesh
- Centre for Mathematical Modelling of Infectious Diseases, Department of Infectious Disease Epidemiology, London School of Hygiene and Tropical Medicine, UK
| |
Collapse
|
7
|
Nieto Ramos A, Fenton FH, Cherry EM. Bayesian inference for fitting cardiac models to experiments: estimating parameter distributions using Hamiltonian Monte Carlo and approximate Bayesian computation. Med Biol Eng Comput 2023; 61:75-95. [PMID: 36322242 DOI: 10.1007/s11517-022-02685-y] [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/15/2022] [Accepted: 10/02/2022] [Indexed: 01/07/2023]
Abstract
Customization of cardiac action potential models has become increasingly important with the recognition of patient-specific models and virtual patient cohorts as valuable predictive tools. Nevertheless, developing customized models by fitting parameters to data poses technical and methodological challenges: despite noise and variability associated with real-world datasets, traditional optimization methods produce a single "best-fit" set of parameter values. Bayesian estimation methods seek distributions of parameter values given the data by obtaining samples from the target distribution, but in practice widely known Bayesian algorithms like Markov chain Monte Carlo tend to be computationally inefficient and scale poorly with the dimensionality of parameter space. In this paper, we consider two computationally efficient Bayesian approaches: the Hamiltonian Monte Carlo (HMC) algorithm and the approximate Bayesian computation sequential Monte Carlo (ABC-SMC) algorithm. We find that both methods successfully identify distributions of model parameters for two cardiac action potential models using model-derived synthetic data and an experimental dataset from a zebrafish heart. Although both methods appear to converge to the same distribution family and are computationally efficient, HMC generally finds narrower marginal distributions, while ABC-SMC is less sensitive to the algorithmic settings including the prior distribution.
Collapse
Affiliation(s)
- Alejandro Nieto Ramos
- School of Mathematical Sciences, Rochester Institute of Technology, 1 Lomb Memorial Drive, 14623, Rochester, NY, USA.,Epilepsy Center, Neurological Institute, Cleveland Clinic, 9500 Euclid Avenue, Cleveland, OH, 44195, USA
| | - Flavio H Fenton
- School of Physics, Georgia Institute of Technology, 837 State Street NW, 30332, Atlanta, GA, USA
| | - Elizabeth M Cherry
- School of Computational Science and Engineering, Georgia Institute of Technology, 756 West Peachtree Street, 30308, Atlanta, GA, USA.
| |
Collapse
|
8
|
Vernon I, Owen J, Aylett-Bullock J, Cuesta-Lazaro C, Frawley J, Quera-Bofarull A, Sedgewick A, Shi D, Truong H, Turner M, Walker J, Caulfield T, Fong K, Krauss F. Bayesian emulation and history matching of JUNE. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20220039. [PMID: 35965471 PMCID: PMC9376712 DOI: 10.1098/rsta.2022.0039] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/18/2022] [Accepted: 06/07/2022] [Indexed: 05/21/2023]
Abstract
We analyze JUNE: a detailed model of COVID-19 transmission with high spatial and demographic resolution, developed as part of the RAMP initiative. JUNE requires substantial computational resources to evaluate, making model calibration and general uncertainty analysis extremely challenging. We describe and employ the uncertainty quantification approaches of Bayes linear emulation and history matching to mimic JUNE and to perform a global parameter search, hence identifying regions of parameter space that produce acceptable matches to observed data, and demonstrating the capability of such methods. This article is part of the theme issue 'Technical challenges of modelling real-life epidemics and examples of overcoming these'.
Collapse
Affiliation(s)
- I. Vernon
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Department of Mathematical Sciences, Durham University, Durham DH13LE, UK
| | - J. Owen
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Department of Mathematical Sciences, Durham University, Durham DH13LE, UK
| | - J. Aylett-Bullock
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Institute for Particle Physics Phenomenology, Durham University, Durham DH13LE, UK
| | - C. Cuesta-Lazaro
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH13LE, UK
| | - J. Frawley
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Advanced Research Computing, Durham University, Durham DH13LE, UK
| | - A. Quera-Bofarull
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH13LE, UK
| | - A. Sedgewick
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Centre for Extragalactic Astronomy, Durham University, Durham DH13LE, UK
| | - D. Shi
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH13LE, UK
| | - H. Truong
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Institute for Particle Physics Phenomenology, Durham University, Durham DH13LE, UK
| | - M. Turner
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Advanced Research Computing, Durham University, Durham DH13LE, UK
| | - J. Walker
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Institute for Particle Physics Phenomenology, Durham University, Durham DH13LE, UK
| | - T. Caulfield
- Department of Computer Science, Durham University, Durham DH13LE, UK
| | - K. Fong
- Department of Science, Technology, Engineering and Public Policy, University College London, London WC1E6BT, UK
- Department of Anaesthesia, University College London Hospital, London NW12BU, UK
| | - F. Krauss
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Institute for Particle Physics Phenomenology, Durham University, Durham DH13LE, UK
| |
Collapse
|
9
|
Vernon I, Owen J, Aylett-Bullock J, Cuesta-Lazaro C, Frawley J, Quera-Bofarull A, Sedgewick A, Shi D, Truong H, Turner M, Walker J, Caulfield T, Fong K, Krauss F. Bayesian emulation and history matching of JUNE. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20210039. [PMID: 35965471 DOI: 10.1098/rsta.2021.0039] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/01/2021] [Accepted: 12/06/2021] [Indexed: 05/21/2023]
Abstract
We analyze JUNE: a detailed model of COVID-19 transmission with high spatial and demographic resolution, developed as part of the RAMP initiative. JUNE requires substantial computational resources to evaluate, making model calibration and general uncertainty analysis extremely challenging. We describe and employ the uncertainty quantification approaches of Bayes linear emulation and history matching to mimic JUNE and to perform a global parameter search, hence identifying regions of parameter space that produce acceptable matches to observed data, and demonstrating the capability of such methods. This article is part of the theme issue 'Technical challenges of modelling real-life epidemics and examples of overcoming these'.
Collapse
Affiliation(s)
- I Vernon
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Department of Mathematical Sciences, Durham University, Durham DH13LE, UK
| | - J Owen
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Department of Mathematical Sciences, Durham University, Durham DH13LE, UK
| | - J Aylett-Bullock
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Institute for Particle Physics Phenomenology, Durham University, Durham DH13LE, UK
| | - C Cuesta-Lazaro
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH13LE, UK
| | - J Frawley
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Advanced Research Computing, Durham University, Durham DH13LE, UK
| | - A Quera-Bofarull
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH13LE, UK
| | - A Sedgewick
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Centre for Extragalactic Astronomy, Durham University, Durham DH13LE, UK
| | - D Shi
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH13LE, UK
| | - H Truong
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Institute for Particle Physics Phenomenology, Durham University, Durham DH13LE, UK
| | - M Turner
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Advanced Research Computing, Durham University, Durham DH13LE, UK
| | - J Walker
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Institute for Particle Physics Phenomenology, Durham University, Durham DH13LE, UK
| | - T Caulfield
- Department of Computer Science, Durham University, Durham DH13LE, UK
| | - K Fong
- Department of Science, Technology, Engineering and Public Policy, University College London, London WC1E6BT, UK
- Department of Anaesthesia, University College London Hospital, London NW12BU, UK
| | - F Krauss
- Institute for Data Science, Durham University, Durham DH13LE, UK
- Institute for Particle Physics Phenomenology, Durham University, Durham DH13LE, UK
| |
Collapse
|
10
|
Du YH, Wang MY, Yang LH, Tong LL, Guo DS, Ji XJ. Optimization and Scale-Up of Fermentation Processes Driven by Models. Bioengineering (Basel) 2022; 9:bioengineering9090473. [PMID: 36135019 PMCID: PMC9495923 DOI: 10.3390/bioengineering9090473] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2022] [Revised: 09/05/2022] [Accepted: 09/09/2022] [Indexed: 11/16/2022] Open
Abstract
In the era of sustainable development, the use of cell factories to produce various compounds by fermentation has attracted extensive attention; however, industrial fermentation requires not only efficient production strains, but also suitable extracellular conditions and medium components, as well as scaling-up. In this regard, the use of biological models has received much attention, and this review will provide guidance for the rapid selection of biological models. This paper first introduces two mechanistic modeling methods, kinetic modeling and constraint-based modeling (CBM), and generalizes their applications in practice. Next, we review data-driven modeling based on machine learning (ML), and highlight the application scope of different learning algorithms. The combined use of ML and CBM for constructing hybrid models is further discussed. At the end, we also discuss the recent strategies for predicting bioreactor scale-up and culture behavior through a combination of biological models and computational fluid dynamics (CFD) models.
Collapse
Affiliation(s)
- Yuan-Hang Du
- School of Food Science and Pharmaceutical Engineering, Nanjing Normal University, Nanjing 210023, China
| | - Min-Yu Wang
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Biotechnology and Pharmaceutical Engineering, Nanjing Tech University, Nanjing 211816, China
| | - Lin-Hui Yang
- School of Food Science and Pharmaceutical Engineering, Nanjing Normal University, Nanjing 210023, China
| | - Ling-Ling Tong
- School of Food Science and Pharmaceutical Engineering, Nanjing Normal University, Nanjing 210023, China
| | - Dong-Sheng Guo
- School of Food Science and Pharmaceutical Engineering, Nanjing Normal University, Nanjing 210023, China
- Correspondence: (D.-S.G.); (X.-J.J.)
| | - Xiao-Jun Ji
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Biotechnology and Pharmaceutical Engineering, Nanjing Tech University, Nanjing 211816, China
- Correspondence: (D.-S.G.); (X.-J.J.)
| |
Collapse
|
11
|
Hu B, Jiang W, Miyagi T, Sun Z, Ekström A, Forssén C, Hagen G, Holt JD, Papenbrock T, Stroberg SR, Vernon I. Ab initio predictions link the neutron skin of 208Pb to nuclear forces. NATURE PHYSICS 2022; 18:1196-1200. [PMID: 36217363 PMCID: PMC9537109 DOI: 10.1038/s41567-022-01715-8] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/23/2021] [Accepted: 07/11/2022] [Indexed: 05/14/2023]
Abstract
Heavy atomic nuclei have an excess of neutrons over protons, which leads to the formation of a neutron skin whose thickness is sensitive to details of the nuclear force. This links atomic nuclei to properties of neutron stars, thereby relating objects that differ in size by orders of magnitude. The nucleus 208Pb is of particular interest because it exhibits a simple structure and is experimentally accessible. However, computing such a heavy nucleus has been out of reach for ab initio theory. By combining advances in quantum many-body methods, statistical tools and emulator technology, we make quantitative predictions for the properties of 208Pb starting from nuclear forces that are consistent with symmetries of low-energy quantum chromodynamics. We explore 109 different nuclear force parameterizations via history matching, confront them with data in select light nuclei and arrive at an importance-weighted ensemble of interactions. We accurately reproduce bulk properties of 208Pb and determine the neutron skin thickness, which is smaller and more precise than a recent extraction from parity-violating electron scattering but in agreement with other experimental probes. This work demonstrates how realistic two- and three-nucleon forces act in a heavy nucleus and allows us to make quantitative predictions across the nuclear landscape.
Collapse
Affiliation(s)
- Baishan Hu
- TRIUMF, Vancouver, British Columbia Canada
| | - Weiguang Jiang
- Department of Physics, Chalmers University of Technology, Gothenburg, Sweden
| | - Takayuki Miyagi
- TRIUMF, Vancouver, British Columbia Canada
- Department of Physics, Technische Universität Darmstadt, Darmstadt, Germany
- ExtreMe Matter Institute EMMI, GSI Helmholtzzentrum für Schwerionenforschung GmbH, Darmstadt, Germany
| | - Zhonghao Sun
- Department of Physics and Astronomy, University of Tennessee, Knoxville, TN USA
- Physics Division, Oak Ridge National Laboratory, Oak Ridge, TN USA
| | - Andreas Ekström
- Department of Physics, Chalmers University of Technology, Gothenburg, Sweden
| | - Christian Forssén
- Department of Physics, Chalmers University of Technology, Gothenburg, Sweden
| | - Gaute Hagen
- TRIUMF, Vancouver, British Columbia Canada
- Department of Physics and Astronomy, University of Tennessee, Knoxville, TN USA
- Physics Division, Oak Ridge National Laboratory, Oak Ridge, TN USA
| | - Jason D. Holt
- TRIUMF, Vancouver, British Columbia Canada
- Department of Physics, McGill University, Montreal, Quebec Canada
| | - Thomas Papenbrock
- Department of Physics and Astronomy, University of Tennessee, Knoxville, TN USA
- Physics Division, Oak Ridge National Laboratory, Oak Ridge, TN USA
| | - S. Ragnar Stroberg
- Department of Physics, University of Washington, Seattle, WA USA
- Physics Division, Argonne National Laboratory, Lemont, IL USA
| | - Ian Vernon
- Department of Mathematical Sciences, Durham University, Durham, UK
| |
Collapse
|
12
|
Dunne M, Mohammadi H, Challenor P, Borgo R, Porphyre T, Vernon I, Firat EE, Turkay C, Torsney-Weir T, Goldstein M, Reeve R, Fang H, Swallow B. Complex model calibration through emulation, a worked example for a stochastic epidemic model. Epidemics 2022; 39:100574. [PMID: 35617882 PMCID: PMC9109972 DOI: 10.1016/j.epidem.2022.100574] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2021] [Revised: 04/22/2022] [Accepted: 04/29/2022] [Indexed: 12/03/2022] Open
Abstract
Uncertainty quantification is a formal paradigm of statistical estimation that aims to account for all uncertainties inherent in the modelling process of real-world complex systems. The methods are directly applicable to stochastic models in epidemiology, however they have thus far not been widely used in this context. In this paper, we provide a tutorial on uncertainty quantification of stochastic epidemic models, aiming to facilitate the use of the uncertainty quantification paradigm for practitioners with other complex stochastic simulators of applied systems. We provide a formal workflow including the important decisions and considerations that need to be taken, and illustrate the methods over a simple stochastic epidemic model of UK SARS-CoV-2 transmission and patient outcome. We also present new approaches to visualisation of outputs from sensitivity analyses and uncertainty quantification more generally in high input and/or output dimensions.
Collapse
Affiliation(s)
- Michael Dunne
- College of Engineering, Mathematics and Physical Sciences, University of Exeter, Exeter, UK
| | - Hossein Mohammadi
- College of Engineering, Mathematics and Physical Sciences, University of Exeter, Exeter, UK
| | - Peter Challenor
- College of Engineering, Mathematics and Physical Sciences, University of Exeter, Exeter, UK
| | - Rita Borgo
- Department of Informatics, King's College London, London, UK
| | - Thibaud Porphyre
- Laboratoire de Biométrie et Biologie Evolutive, VetAgro Sup, Marcy l'Etoile, France
| | - Ian Vernon
- Department of Mathematical Sciences, Durham University, Durham, UK
| | - Elif E Firat
- Department of Computer Science, University of Nottingham, Nottingham, UK
| | - Cagatay Turkay
- Centre for Interdisciplinary Methodologies, University of Warwick, Coventry, UK
| | - Thomas Torsney-Weir
- VRVis Zentrum für Virtual Reality und Visualisierung Forschungs-GmbH, Vienna, Austria
| | | | - Richard Reeve
- Boyd Orr Centre for Population and Ecosystem Health, Institute of Biodiversity, Animal Health and Comparative Medicine, College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow, UK
| | - Hui Fang
- Department of Computer Science, Loughborough University, Loughborough, UK
| | - Ben Swallow
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK.
| |
Collapse
|
13
|
Schöpping M, Gaspar P, Neves AR, Franzén CJ, Zeidan AA. Identifying the essential nutritional requirements of the probiotic bacteria Bifidobacterium animalis and Bifidobacterium longum through genome-scale modeling. NPJ Syst Biol Appl 2021; 7:47. [PMID: 34887435 PMCID: PMC8660834 DOI: 10.1038/s41540-021-00207-4] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2021] [Accepted: 11/03/2021] [Indexed: 12/14/2022] Open
Abstract
Although bifidobacteria are widely used as probiotics, their metabolism and physiology remain to be explored in depth. In this work, strain-specific genome-scale metabolic models were developed for two industrially and clinically relevant bifidobacteria, Bifidobacterium animalis subsp. lactis BB-12® and B. longum subsp. longum BB-46, and subjected to iterative cycles of manual curation and experimental validation. A constraint-based modeling framework was used to probe the metabolic landscape of the strains and identify their essential nutritional requirements. Both strains showed an absolute requirement for pantethine as a precursor for coenzyme A biosynthesis. Menaquinone-4 was found to be essential only for BB-46 growth, whereas nicotinic acid was only required by BB-12®. The model-generated insights were used to formulate a chemically defined medium that supports the growth of both strains to the same extent as a complex culture medium. Carbohydrate utilization profiles predicted by the models were experimentally validated. Furthermore, model predictions were quantitatively validated in the newly formulated medium in lab-scale batch fermentations. The models and the formulated medium represent valuable tools to further explore the metabolism and physiology of the two species, investigate the mechanisms underlying their health-promoting effects and guide the optimization of their industrial production processes.
Collapse
Affiliation(s)
- Marie Schöpping
- Systems Biology, Discovery, Chr. Hansen A/S, 2970, Hørsholm, Denmark
- Division of Industrial Biotechnology, Department of Biology and Biological Engineering, Chalmers University of Technology, 41296, Gothenburg, Sweden
| | - Paula Gaspar
- Systems Biology, Discovery, Chr. Hansen A/S, 2970, Hørsholm, Denmark
| | - Ana Rute Neves
- Systems Biology, Discovery, Chr. Hansen A/S, 2970, Hørsholm, Denmark
- Arla Foods Ingredients Group P/S, 6920, Videbæk, Denmark
| | - Carl Johan Franzén
- Division of Industrial Biotechnology, Department of Biology and Biological Engineering, Chalmers University of Technology, 41296, Gothenburg, Sweden
| | - Ahmad A Zeidan
- Systems Biology, Discovery, Chr. Hansen A/S, 2970, Hørsholm, Denmark.
| |
Collapse
|
14
|
Aylett-Bullock J, Cuesta-Lazaro C, Quera-Bofarull A, Icaza-Lizaola M, Sedgewick A, Truong H, Curran A, Elliott E, Caulfield T, Fong K, Vernon I, Williams J, Bower R, Krauss F. June: open-source individual-based epidemiology simulation. ROYAL SOCIETY OPEN SCIENCE 2021. [PMID: 34295529 DOI: 10.5281/zenodo.4925939] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
We introduce June, an open-source framework for the detailed simulation of epidemics on the basis of social interactions in a virtual population constructed from geographically granular census data, reflecting age, sex, ethnicity and socio-economic indicators. Interactions between individuals are modelled in groups of various sizes and properties, such as households, schools and workplaces, and other social activities using social mixing matrices. June provides a suite of flexible parametrizations that describe infectious diseases, how they are transmitted and affect contaminated individuals. In this paper, we apply June to the specific case of modelling the spread of COVID-19 in England. We discuss the quality of initial model outputs which reproduce reported hospital admission and mortality statistics at national and regional levels as well as by age strata.
Collapse
Affiliation(s)
- Joseph Aylett-Bullock
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Particle Physics Phenomenology, Durham University, Durham DH1 3LE, UK
| | - Carolina Cuesta-Lazaro
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH1 3LE, UK
| | - Arnau Quera-Bofarull
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH1 3LE, UK
| | - Miguel Icaza-Lizaola
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH1 3LE, UK
| | - Aidan Sedgewick
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Centre for Extragalactic Astronomy, Durham University, Durham DH1 3LE, UK
| | - Henry Truong
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Particle Physics Phenomenology, Durham University, Durham DH1 3LE, UK
| | - Aoife Curran
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH1 3LE, UK
| | - Edward Elliott
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH1 3LE, UK
| | - Tristan Caulfield
- Department of Computer Science, University College London, London WC1E 6BT, UK
| | - Kevin Fong
- Department of Science, Technology, Engineering and Public Policy, University College London, London WC1E 6BT, UK
- Department of Anaesthesia, University College London Hospital, London NW1 2BU, UK
| | - Ian Vernon
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Department of Mathematical Sciences, Durham University, Durham DH1 3LE, UK
| | - Julian Williams
- Institute of Hazard, Risk and Resilience, Durham University, Durham DH1 3LE, UK
| | - Richard Bower
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH1 3LE, UK
| | - Frank Krauss
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Particle Physics Phenomenology, Durham University, Durham DH1 3LE, UK
| |
Collapse
|
15
|
Aylett-Bullock J, Cuesta-Lazaro C, Quera-Bofarull A, Icaza-Lizaola M, Sedgewick A, Truong H, Curran A, Elliott E, Caulfield T, Fong K, Vernon I, Williams J, Bower R, Krauss F. June: open-source individual-based epidemiology simulation. ROYAL SOCIETY OPEN SCIENCE 2021; 8:210506. [PMID: 34295529 PMCID: PMC8261230 DOI: 10.1098/rsos.210506] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/25/2021] [Accepted: 06/22/2021] [Indexed: 05/09/2023]
Abstract
We introduce June, an open-source framework for the detailed simulation of epidemics on the basis of social interactions in a virtual population constructed from geographically granular census data, reflecting age, sex, ethnicity and socio-economic indicators. Interactions between individuals are modelled in groups of various sizes and properties, such as households, schools and workplaces, and other social activities using social mixing matrices. June provides a suite of flexible parametrizations that describe infectious diseases, how they are transmitted and affect contaminated individuals. In this paper, we apply June to the specific case of modelling the spread of COVID-19 in England. We discuss the quality of initial model outputs which reproduce reported hospital admission and mortality statistics at national and regional levels as well as by age strata.
Collapse
Affiliation(s)
- Joseph Aylett-Bullock
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Particle Physics Phenomenology, Durham University, Durham DH1 3LE, UK
| | - Carolina Cuesta-Lazaro
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH1 3LE, UK
| | - Arnau Quera-Bofarull
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH1 3LE, UK
| | - Miguel Icaza-Lizaola
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH1 3LE, UK
| | - Aidan Sedgewick
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Centre for Extragalactic Astronomy, Durham University, Durham DH1 3LE, UK
| | - Henry Truong
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Particle Physics Phenomenology, Durham University, Durham DH1 3LE, UK
| | - Aoife Curran
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH1 3LE, UK
| | - Edward Elliott
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH1 3LE, UK
| | - Tristan Caulfield
- Department of Computer Science, University College London, London WC1E 6BT, UK
| | - Kevin Fong
- Department of Science, Technology, Engineering and Public Policy, University College London, London WC1E 6BT, UK
- Department of Anaesthesia, University College London Hospital, London NW1 2BU, UK
| | - Ian Vernon
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Department of Mathematical Sciences, Durham University, Durham DH1 3LE, UK
| | - Julian Williams
- Institute of Hazard, Risk and Resilience, Durham University, Durham DH1 3LE, UK
| | - Richard Bower
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Computational Cosmology, Durham University, Durham DH1 3LE, UK
| | - Frank Krauss
- Institute for Data Science, Durham University, Durham DH1 3LE, UK
- Institute for Particle Physics Phenomenology, Durham University, Durham DH1 3LE, UK
| |
Collapse
|
16
|
Volodina V, Challenor P. The importance of uncertainty quantification in model reproducibility. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2021; 379:20200071. [PMID: 33775141 PMCID: PMC8059558 DOI: 10.1098/rsta.2020.0071] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Accepted: 09/21/2020] [Indexed: 06/12/2023]
Abstract
Many computer models possess high-dimensional input spaces and substantial computational time to produce a single model evaluation. Although such models are often 'deterministic', these models suffer from a wide range of uncertainties. We argue that uncertainty quantification is crucial for computer model validation and reproducibility. We present a statistical framework, termed history matching, for performing global parameter search by comparing model output to the observed data. We employ Gaussian process (GP) emulators to produce fast predictions about model behaviour at the arbitrary input parameter settings allowing output uncertainty distributions to be calculated. History matching identifies sets of input parameters that give rise to acceptable matches between observed data and model output given our representation of uncertainties. Modellers could proceed by simulating computer models' outputs of interest at these identified parameter settings and producing a range of predictions. The variability in model results is crucial for inter-model comparison as well as model development. We illustrate the performance of emulation and history matching on a simple one-dimensional toy model and in application to a climate model. This article is part of the theme issue 'Reliability and reproducibility in computational science: implementing verification, validation and uncertainty quantification in silico'.
Collapse
Affiliation(s)
- Victoria Volodina
- The Alan Turing Institute, British Library, 96 Euston Road, London, NW1 2DB, UK
| | - Peter Challenor
- The Alan Turing Institute, British Library, 96 Euston Road, London, NW1 2DB, UK
- College of Engineering, Mathematics and Physical Sciences, University of Exeter, Exeter EX4 4QE, UK
| |
Collapse
|
17
|
Sangha GS, Goergen CJ, Prior SJ, Ranadive SM, Clyne AM. Preclinical techniques to investigate exercise training in vascular pathophysiology. Am J Physiol Heart Circ Physiol 2021; 320:H1566-H1600. [PMID: 33385323 PMCID: PMC8260379 DOI: 10.1152/ajpheart.00719.2020] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Atherosclerosis is a dynamic process starting with endothelial dysfunction and inflammation and eventually leading to life-threatening arterial plaques. Exercise generally improves endothelial function in a dose-dependent manner by altering hemodynamics, specifically by increased arterial pressure, pulsatility, and shear stress. However, athletes who regularly participate in high-intensity training can develop arterial plaques, suggesting alternative mechanisms through which excessive exercise promotes vascular disease. Understanding the mechanisms that drive atherosclerosis in sedentary versus exercise states may lead to novel rehabilitative methods aimed at improving exercise compliance and physical activity. Preclinical tools, including in vitro cell assays, in vivo animal models, and in silico computational methods, broaden our capabilities to study the mechanisms through which exercise impacts atherogenesis, from molecular maladaptation to vascular remodeling. Here, we describe how preclinical research tools have and can be used to study exercise effects on atherosclerosis. We then propose how advanced bioengineering techniques can be used to address gaps in our current understanding of vascular pathophysiology, including integrating in vitro, in vivo, and in silico studies across multiple tissue systems and size scales. Improving our understanding of the antiatherogenic exercise effects will enable engaging, targeted, and individualized exercise recommendations to promote cardiovascular health rather than treating cardiovascular disease that results from a sedentary lifestyle.
Collapse
Affiliation(s)
- Gurneet S Sangha
- Fischell Department of Bioengineering, University of Maryland, College Park, Maryland
| | - Craig J Goergen
- Weldon School of Biomedical Engineering, Purdue University, West Lafayette, Indiana.,Purdue University Center for Cancer Research, Purdue University, West Lafayette, Indiana
| | - Steven J Prior
- Department of Kinesiology, University of Maryland School of Public Health, College Park, Maryland.,Baltimore Veterans Affairs Geriatric Research, Education, and Clinical Center, Baltimore, Maryland
| | - Sushant M Ranadive
- Department of Kinesiology, University of Maryland School of Public Health, College Park, Maryland
| | - Alisa M Clyne
- Fischell Department of Bioengineering, University of Maryland, College Park, Maryland
| |
Collapse
|
18
|
Rowe JH, Jones AM. Focus on biosensors: Looking through the lens of quantitative biology. QUANTITATIVE PLANT BIOLOGY 2021; 2:e12. [PMID: 37077214 PMCID: PMC10095858 DOI: 10.1017/qpb.2021.10] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/01/2021] [Revised: 07/27/2021] [Accepted: 07/27/2021] [Indexed: 05/02/2023]
Abstract
In recent years, plant biologists interested in quantifying molecules and molecular events in vivo have started to complement reporter systems with genetically encoded fluorescent biosensors (GEFBs) that directly sense an analyte. Such biosensors can allow measurements at the level of individual cells and over time. This information is proving valuable to mathematical modellers interested in representing biological phenomena in silico, because improved measurements can guide improved model construction and model parametrisation. Advances in synthetic biology have accelerated the pace of biosensor development, and the simultaneous expression of spectrally compatible biosensors now allows quantification of multiple nodes in signalling networks. For biosensors that directly respond to stimuli, targeting to specific cellular compartments allows the observation of differential accumulation of analytes in distinct organelles, bringing insights to reactive oxygen species/calcium signalling and photosynthesis research. In conjunction with improved image analysis methods, advances in biosensor imaging can help close the loop between experimentation and mathematical modelling.
Collapse
Affiliation(s)
- James H. Rowe
- Sainsbury Laboratory, Cambridge University, Cambridge, United Kingdom
| | - Alexander M. Jones
- Sainsbury Laboratory, Cambridge University, Cambridge, United Kingdom
- Author for correspondence: Alexander M. Jones, E-mail:
| |
Collapse
|
19
|
Tong H, Madison I, Long TA, Williams CM. Computational solutions for modeling and controlling plant response to abiotic stresses: a review with focus on iron deficiency. CURRENT OPINION IN PLANT BIOLOGY 2020; 57:8-15. [PMID: 32619968 DOI: 10.1016/j.pbi.2020.05.006] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/30/2020] [Revised: 05/15/2020] [Accepted: 05/23/2020] [Indexed: 06/11/2023]
Abstract
Computational solutions enable plant scientists to model protein-mediated stress responses and characterize novel gene functions that coordinate responses to a variety of abiotic stress conditions. Recently, density functional theory was used to study proteins active sites and elucidate enzyme conversion mechanisms involved in iron deficiency responsive signaling pathways. Computational approaches for protein homology modeling and the kinetic modeling of signaling pathways have also resolved the identity and function in proteins involved in iron deficiency signaling pathways. Significant changes in gene relationships under other stress conditions, such as heat or drought stress, have been recently identified using differential network analysis, suggesting that stress tolerance is achieved through asynchronous control. Moreover, the increasing development and use of statistical modeling and systematic modeling of transcriptomic data have provided significant insight into the gene regulatory mechanisms associated with abiotic stress responses. These types of in silico approaches have facilitated the plant science community's future goals of developing multi-scale models of responses to iron deficiency stress and other abiotic stress conditions.
Collapse
Affiliation(s)
- Haonan Tong
- Electrical and Computer Engineering, North Carolina State University, Raleigh, USA
| | - Imani Madison
- Plant and Microbial Biology, North Carolina State University, Raleigh, USA
| | - Terri A Long
- Plant and Microbial Biology, North Carolina State University, Raleigh, USA.
| | - Cranos M Williams
- Electrical and Computer Engineering, North Carolina State University, Raleigh, USA.
| |
Collapse
|
20
|
Abstract
This review aims to reflect upon the major developments in PARP14 research from late 2017 to early 2020. In doing so, this report will focus on the continual elucidation of PARP14's function including an emerging role in viral replication. This is in addition to other functional developments in cancer and inflammation, along with reflecting upon the leads in inhibitor design, including the increased attention toward the macrodomain. This report will also include a brief recap on contemporary poly(ADP-ribose) polymerase inhibitors and reflect upon the development surrounding the other poly(ADP-ribose) polymerases to overall give a succinct update to assist the development of selective PARP14 inhibitors.
Collapse
Affiliation(s)
- Amanda L Tauber
- Faculty of Health Sciences & Medicine, Bond University, Gold Coast 4229, Queensland, Australia
| | - Stephan M Levonis
- Faculty of Health Sciences & Medicine, Bond University, Gold Coast 4229, Queensland, Australia
| | - Stephanie S Schweiker
- Faculty of Health Sciences & Medicine, Bond University, Gold Coast 4229, Queensland, Australia
| |
Collapse
|
21
|
Marucci L, Barberis M, Karr J, Ray O, Race PR, de Souza Andrade M, Grierson C, Hoffmann SA, Landon S, Rech E, Rees-Garbutt J, Seabrook R, Shaw W, Woods C. Computer-Aided Whole-Cell Design: Taking a Holistic Approach by Integrating Synthetic With Systems Biology. Front Bioeng Biotechnol 2020; 8:942. [PMID: 32850764 PMCID: PMC7426639 DOI: 10.3389/fbioe.2020.00942] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2020] [Accepted: 07/21/2020] [Indexed: 01/03/2023] Open
Abstract
Computer-aided design (CAD) for synthetic biology promises to accelerate the rational and robust engineering of biological systems. It requires both detailed and quantitative mathematical and experimental models of the processes to (re)design biology, and software and tools for genetic engineering and DNA assembly. Ultimately, the increased precision in the design phase will have a dramatic impact on the production of designer cells and organisms with bespoke functions and increased modularity. CAD strategies require quantitative models of cells that can capture multiscale processes and link genotypes to phenotypes. Here, we present a perspective on how whole-cell, multiscale models could transform design-build-test-learn cycles in synthetic biology. We show how these models could significantly aid in the design and learn phases while reducing experimental testing by presenting case studies spanning from genome minimization to cell-free systems. We also discuss several challenges for the realization of our vision. The possibility to describe and build whole-cells in silico offers an opportunity to develop increasingly automatized, precise and accessible CAD tools and strategies.
Collapse
Affiliation(s)
- Lucia Marucci
- Department of Engineering Mathematics, University of Bristol, Bristol, United Kingdom.,School of Cellular and Molecular Medicine, University of Bristol, Bristol, United Kingdom.,Bristol Centre for Synthetic Biology (BrisSynBio), University of Bristol, Bristol, United Kingdom
| | - Matteo Barberis
- Systems Biology, School of Biosciences and Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford, United Kingdom.,Centre for Mathematical and Computational Biology, CMCB, University of Surrey, Guildford, United Kingdom.,Synthetic Systems Biology and Nuclear Organization, Swammerdam Institute for Life Sciences, University of Amsterdam, Amsterdam, Netherlands
| | - Jonathan Karr
- Icahn Institute for Data Science and Genomic Technology, Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States
| | - Oliver Ray
- Department of Computer Science, University of Bristol, Bristol, United Kingdom
| | - Paul R Race
- Bristol Centre for Synthetic Biology (BrisSynBio), University of Bristol, Bristol, United Kingdom.,School of Biochemistry, University of Bristol, Bristol, United Kingdom
| | - Miguel de Souza Andrade
- Brazilian Agricultural Research Corporation/National Institute of Science and Technology - Synthetic Biology, Brasília, Brazil.,Department of Cell Biology, Institute of Biological Sciences, University of Brasília, Brasília, Brazil
| | - Claire Grierson
- Bristol Centre for Synthetic Biology (BrisSynBio), University of Bristol, Bristol, United Kingdom.,School of Biological Sciences, University of Bristol, Bristol, United Kingdom
| | - Stefan Andreas Hoffmann
- Manchester Institute of Biotechnology, The University of Manchester, Manchester, United Kingdom
| | - Sophie Landon
- Department of Engineering Mathematics, University of Bristol, Bristol, United Kingdom.,Bristol Centre for Synthetic Biology (BrisSynBio), University of Bristol, Bristol, United Kingdom
| | - Elibio Rech
- Brazilian Agricultural Research Corporation/National Institute of Science and Technology - Synthetic Biology, Brasília, Brazil
| | - Joshua Rees-Garbutt
- Bristol Centre for Synthetic Biology (BrisSynBio), University of Bristol, Bristol, United Kingdom.,School of Biological Sciences, University of Bristol, Bristol, United Kingdom
| | - Richard Seabrook
- Elizabeth Blackwell Institute for Health Research (EBI), University of Bristol, Bristol, United Kingdom
| | - William Shaw
- Department of Bioengineering, Imperial College London, London, United Kingdom
| | - Christopher Woods
- Bristol Centre for Synthetic Biology (BrisSynBio), University of Bristol, Bristol, United Kingdom.,School of Chemistry, University of Bristol, Bristol, United Kingdom
| |
Collapse
|
22
|
Jackson SE, Vernon I, Liu J, Lindsey K. Understanding hormonal crosstalk in Arabidopsis root development via emulation and history matching. Stat Appl Genet Mol Biol 2020; 19:sagmb-2018-0053. [PMID: 32649296 DOI: 10.1515/sagmb-2018-0053] [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: 09/07/2018] [Accepted: 05/12/2020] [Indexed: 11/15/2022]
Abstract
A major challenge in plant developmental biology is to understand how plant growth is coordinated by interacting hormones and genes. To meet this challenge, it is important to not only use experimental data, but also formulate a mathematical model. For the mathematical model to best describe the true biological system, it is necessary to understand the parameter space of the model, along with the links between the model, the parameter space and experimental observations. We develop sequential history matching methodology, using Bayesian emulation, to gain substantial insight into biological model parameter spaces. This is achieved by finding sets of acceptable parameters in accordance with successive sets of physical observations. These methods are then applied to a complex hormonal crosstalk model for Arabidopsis root growth. In this application, we demonstrate how an initial set of 22 observed trends reduce the volume of the set of acceptable inputs to a proportion of 6.1 × 10-7 of the original space. Additional sets of biologically relevant experimental data, each of size 5, reduce the size of this space by a further three and two orders of magnitude respectively. Hence, we provide insight into the constraints placed upon the model structure by, and the biological consequences of, measuring subsets of observations.
Collapse
Affiliation(s)
- Samuel E Jackson
- Southampton Statistical Sciences Research Institute, University of Southampton, Southampton, UK
| | - Ian Vernon
- Department of Mathematical Sciences, Durham University, Durham, UK
| | - Junli Liu
- School of Biological and Biomedical Sciences, Durham University, Durham, UK
| | - Keith Lindsey
- School of Biological and Biomedical Sciences, Durham University, Durham, UK
| |
Collapse
|
23
|
Reanalysis and integration of public microarray datasets reveals novel host genes modulated in leprosy. Mol Genet Genomics 2020; 295:1355-1368. [PMID: 32661593 DOI: 10.1007/s00438-020-01705-6] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2020] [Accepted: 07/01/2020] [Indexed: 01/24/2023]
Abstract
Due to multiple hypothesis testing with often limited sample size, microarrays and other-omics technologies can sometimes produce irreproducible findings. Complementary to better experimental design, reanalysis and integration of gene expression datasets may help overcome reproducibility issues by identifying consistent differentially expressed genes from independent studies. In this work, after a systematic search, nine microarray datasets evaluating host gene expression in leprosy were reanalyzed and the information was integrated to strengthen evidence of differential expression for several genes. Our results are relevant in prioritizing genes and pathways for further investigation, whether in functional studies or in biomarker discovery. Reanalysis of individual datasets revealed several differentially expressed genes (DEGs) in accordance with original reports. Then, five integration methods (P value and effect size based) were tested. In the end, random-effects model and ratio association were selected as the main methods to pinpoint DEGs. Overall, classic pathways were found corroborating previous findings and validating this approach. Also, we identified some novel DEG involved especially with skin development processes (AQP3, AKR1C3, CYP27B1, LTB, VDR) and keratinocyte biology (CSTA, DSG1, KRT14, KRT5, PKP1, IVL), both still poorly understood in leprosy context. In addition, here we provide aggregated evidence towards some gene candidates that should be prioritized in further leprosy research, as they are likely important in immunopathogenesis. Altogether, these data are useful in better understanding host responses to the disease and, at the same time, provide a list of potential host biomarkers that could be useful in complementing leprosy diagnosis based on transcriptional levels.
Collapse
|
24
|
Longobardi S, Lewalle A, Coveney S, Sjaastad I, Espe EKS, Louch WE, Musante CJ, Sher A, Niederer SA. Predicting left ventricular contractile function via Gaussian process emulation in aortic-banded rats. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2020; 378:20190334. [PMID: 32448071 PMCID: PMC7287330 DOI: 10.1098/rsta.2019.0334] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
Cardiac contraction is the result of integrated cellular, tissue and organ function. Biophysical in silico cardiac models offer a systematic approach for studying these multi-scale interactions. The computational cost of such models is high, due to their multi-parametric and nonlinear nature. This has so far made it difficult to perform model fitting and prevented global sensitivity analysis (GSA) studies. We propose a machine learning approach based on Gaussian process emulation of model simulations using probabilistic surrogate models, which enables model parameter inference via a Bayesian history matching (HM) technique and GSA on whole-organ mechanics. This framework is applied to model healthy and aortic-banded hypertensive rats, a commonly used animal model of heart failure disease. The obtained probabilistic surrogate models accurately predicted the left ventricular pump function (R2 = 0.92 for ejection fraction). The HM technique allowed us to fit both the control and diseased virtual bi-ventricular rat heart models to magnetic resonance imaging and literature data, with model outputs from the constrained parameter space falling within 2 SD of the respective experimental values. The GSA identified Troponin C and cross-bridge kinetics as key parameters in determining both systolic and diastolic ventricular function. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.
Collapse
Affiliation(s)
- S Longobardi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - A Lewalle
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - S Coveney
- Insigneo Institute for in-silico Medicine and Department of Computer Science, University of Sheffield, Sheffield, UK
| | - I Sjaastad
- Institute for Experimental Medical Research and KG Jebsen Center for Cardiac Research, University of Oslo, Oslo, Norway
| | - E K S Espe
- Institute for Experimental Medical Research and KG Jebsen Center for Cardiac Research, University of Oslo, Oslo, Norway
| | - W E Louch
- Institute for Experimental Medical Research and KG Jebsen Center for Cardiac Research, University of Oslo, Oslo, Norway
| | - C J Musante
- Pfizer Worldwide Research, Development and Medical, Cambridge, MA, USA
| | - A Sher
- Pfizer Worldwide Research, Development and Medical, Cambridge, MA, USA
| | - S A Niederer
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| |
Collapse
|
25
|
Sanchez-Taltavull D, Perkins TJ, Dommann N, Melin N, Keogh A, Candinas D, Stroka D, Beldi G. Bayesian correlation is a robust gene similarity measure for single-cell RNA-seq data. NAR Genom Bioinform 2020; 2:lqaa002. [PMID: 33575552 PMCID: PMC7671344 DOI: 10.1093/nargab/lqaa002] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2019] [Revised: 11/30/2019] [Accepted: 01/09/2020] [Indexed: 02/07/2023] Open
Abstract
Assessing similarity is highly important for bioinformatics algorithms to determine correlations between biological information. A common problem is that similarity can appear by chance, particularly for low expressed entities. This is especially relevant in single-cell RNA-seq (scRNA-seq) data because read counts are much lower compared to bulk RNA-seq. Recently, a Bayesian correlation scheme that assigns low similarity to genes that have low confidence expression estimates has been proposed to assess similarity for bulk RNA-seq. Our goal is to extend the properties of the Bayesian correlation in scRNA-seq data by considering three ways to compute similarity. First, we compute the similarity of pairs of genes over all cells. Second, we identify specific cell populations and compute the correlation in those populations. Third, we compute the similarity of pairs of genes over all clusters, by considering the total mRNA expression. We demonstrate that Bayesian correlations are more reproducible than Pearson correlations. Compared to Pearson correlations, Bayesian correlations have a smaller dependence on the number of input cells. We show that the Bayesian correlation algorithm assigns high similarity values to genes with a biological relevance in a specific population. We conclude that Bayesian correlation is a robust similarity measure in scRNA-seq data.
Collapse
Affiliation(s)
- Daniel Sanchez-Taltavull
- Visceral Surgery and Medicine, Inselspital, Bern University Hospital, Department for BioMedical Research, University of Bern, Murtenstrasse 35, 3008 Bern, Switzerland
| | - Theodore J Perkins
- Regenerative Medicine Program, Ottawa Hospital Research Institute, Ottawa, Ontario, ON K1H8L6, Canada.,Department of Biochemistry, Microbiology and Immunology, University of Ottawa, Ottawa, Ontario, ON K1H8L6, Canada
| | - Noelle Dommann
- Visceral Surgery and Medicine, Inselspital, Bern University Hospital, Department for BioMedical Research, University of Bern, Murtenstrasse 35, 3008 Bern, Switzerland
| | - Nicolas Melin
- Visceral Surgery and Medicine, Inselspital, Bern University Hospital, Department for BioMedical Research, University of Bern, Murtenstrasse 35, 3008 Bern, Switzerland
| | - Adrian Keogh
- Visceral Surgery and Medicine, Inselspital, Bern University Hospital, Department for BioMedical Research, University of Bern, Murtenstrasse 35, 3008 Bern, Switzerland
| | - Daniel Candinas
- Visceral Surgery and Medicine, Inselspital, Bern University Hospital, Department for BioMedical Research, University of Bern, Murtenstrasse 35, 3008 Bern, Switzerland
| | - Deborah Stroka
- Visceral Surgery and Medicine, Inselspital, Bern University Hospital, Department for BioMedical Research, University of Bern, Murtenstrasse 35, 3008 Bern, Switzerland
| | - Guido Beldi
- Visceral Surgery and Medicine, Inselspital, Bern University Hospital, Department for BioMedical Research, University of Bern, Murtenstrasse 35, 3008 Bern, Switzerland
| |
Collapse
|
26
|
Alden K, Cosgrove J, Coles M, Timmis J. Using Emulation to Engineer and Understand Simulations of Biological Systems. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2020; 17:302-315. [PMID: 29994223 DOI: 10.1109/tcbb.2018.2843339] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Modeling and simulation techniques have demonstrated success in studying biological systems. As the drive to better capture biological complexity leads to more sophisticated simulators, it becomes challenging to perform statistical analyses that help translate predictions into increased understanding. These analyses may require repeated executions and extensive sampling of high-dimensional parameter spaces: analyses that may become intractable due to time and resource limitations. Significant reduction in these requirements can be obtained using surrogate models, or emulators, that can rapidly and accurately predict the output of an existing simulator. We apply emulation to evaluate and enrich understanding of a previously published agent-based simulator of lymphoid tissue organogenesis, showing an ensemble of machine learning techniques can reproduce results obtained using a suite of statistical analyses within seconds. This performance improvement permits incorporation of previously intractable analyses, including multi-objective optimization to obtain parameter sets that yield a desired response, and Approximate Bayesian Computation to assess parametric uncertainty. To facilitate exploitation of emulation in simulation-focused studies, we extend our open source statistical package, spartan, to provide a suite of tools for emulator development, validation, and application. Overcoming resource limitations permits enriched evaluation and refinement, easing translation of simulator insights into increased biological understanding.
Collapse
|
27
|
Leng J, Shoura M, McLeish TCB, Real AN, Hardey M, McCafferty J, Ranson NA, Harris SA. Securing the future of research computing in the biosciences. PLoS Comput Biol 2019; 15:e1006958. [PMID: 31095554 PMCID: PMC6521984 DOI: 10.1371/journal.pcbi.1006958] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
Improvements in technology often drive scientific discovery. Therefore, research requires sustained investment in the latest equipment and training for the researchers who are going to use it. Prioritising and administering infrastructure investment is challenging because future needs are difficult to predict. In the past, highly computationally demanding research was associated primarily with particle physics and astronomy experiments. However, as biology becomes more quantitative and bioscientists generate more and more data, their computational requirements may ultimately exceed those of physical scientists. Computation has always been central to bioinformatics, but now imaging experiments have rapidly growing data processing and storage requirements. There is also an urgent need for new modelling and simulation tools to provide insight and understanding of these biophysical experiments. Bioscience communities must work together to provide the software and skills training needed in their areas. Research-active institutions need to recognise that computation is now vital in many more areas of discovery and create an environment where it can be embraced. The public must also become aware of both the power and limitations of computing, particularly with respect to their health and personal data.
Collapse
Affiliation(s)
- Joanna Leng
- School of Computing, University of Leeds, Leeds, United Kingdom
| | - Massa Shoura
- School of Pathology, Stanford University, Palo Alto, California, United States of America
| | | | - Alan N. Real
- Advanced Research Computing, University of Durham, Durham, United Kingdom
| | - Mariann Hardey
- Advanced Research Computing, University of Durham, Durham, United Kingdom
- School of Business, University of Durham, Durham, United Kingdom
| | | | - Neil A. Ranson
- Astbury Centre for Structural and Molecular Biology, University of Leeds, Leeds, United Kingdom
- School of Molecular and Cellular Biology, University of Leeds, Leeds, United Kingdom
| | - Sarah A. Harris
- Astbury Centre for Structural and Molecular Biology, University of Leeds, Leeds, United Kingdom
- School of Physics and Astronomy, University of Leeds, Leeds, United Kingdom
- * E-mail:
| |
Collapse
|
28
|
Abstract
Alterations in the human gut microbiota play an important role in disease pathogenesis. Although next-generation sequencing has provided observational evidence linking shifts in gut microbiota composition to alterations in the human host, underlying mechanisms remain elusive. Metabolites generated within complex microbial communities and at the crossroads with host cells may be able to explain the impact of the gut microbiome on human homeostasis. Emerging technologies including novel culturing protocols, microfluidic systems, engineered organoids, and single-cell imaging approaches are providing new perspectives from which the gut microbiome can be studied paving the way to new diagnostic markers and personalized therapeutic interventions.
Collapse
Affiliation(s)
- Paola Brun
- Department of Molecular Medicine, University of Padova, Padova, Italy
| |
Collapse
|
29
|
Rutten JP, Ten Tusscher K. In Silico Roots: Room for Growth. TRENDS IN PLANT SCIENCE 2019; 24:250-262. [PMID: 30665820 DOI: 10.1016/j.tplants.2018.11.005] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/31/2018] [Revised: 11/14/2018] [Accepted: 11/19/2018] [Indexed: 06/09/2023]
Abstract
Computational models are invaluable tools for understanding the hormonal and genetic control of root development. Thus far, models have focused on the crucial roles that auxin transport and metabolism play in determining the auxin signaling gradient that controls the root meristem. Other hormones such as cytokinins, gibberellins, and ethylene have predominantly been considered as modulators of auxin dynamics, but their underlying patterning mechanisms are currently unresolved. In addition, the effects of cell- and tissue-level growth dynamics, which induce dilution and displacement of signaling molecules, have remained unexplored. Elucidating these additional mechanisms will be essential to unravel how root growth is patterned in a robust and self-organized manner. Models incorporating growth will thus be crucial in unraveling the underlying logic of root developmental decision making.
Collapse
Affiliation(s)
- Jacob Pieter Rutten
- Computational Developmental Biology Group, Faculty of Science, Utrecht University, Utrecht, The Netherlands
| | - Kirsten Ten Tusscher
- Computational Developmental Biology Group, Faculty of Science, Utrecht University, Utrecht, The Netherlands.
| |
Collapse
|
30
|
Coveney S, Clayton RH. Fitting two human atrial cell models to experimental data using Bayesian history matching. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2018; 139:43-58. [PMID: 30145156 DOI: 10.1016/j.pbiomolbio.2018.08.001] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/19/2018] [Revised: 08/02/2018] [Accepted: 08/06/2018] [Indexed: 12/18/2022]
Abstract
Cardiac cell models are potentially valuable tools for applications such as quantitative safety pharmacology, but have many parameters. Action potentials in real cardiac cells also vary from beat to beat, and from one cell to another. Calibrating cardiac cell models to experimental observations is difficult, because the parameter space is large and high-dimensional. In this study we have demonstrated the use of history matching to calibrate the maximum conductance of ion channels and exchangers in two detailed models of the human atrial action potential against measurements of action potential biomarkers. History matching is an approach developed in other modelling communities, based on constructing fast-running Gaussian process emulators of the model. Emulators were constructed from a small number of model runs (around 102), and then run many times (>106) at low computational cost, each time with a different set of model parameters. Emulator outputs were compared with experimental biomarkers using an implausibility measure, which took into account experimental variance as well as emulator variance. By repeating this process, the region of non-implausible parameter space was iteratively reduced. Both cardiac cell models were successfully calibrated to experimental datasets, resulting in sets of parameters that could be sampled to produce variable action potentials. However, model parameters did not occupy a small range of values. Instead, the history matching process exposed inputs that can co-vary across a wide range and still be consistent with a particular biomarker. We also found correlations between some biomarkers, indicating a need for better descriptors of action potential shape.
Collapse
Affiliation(s)
- Sam Coveney
- Insigneo Institute for in-silico Medicine and Department of Computer Science, University of Sheffield, UK.
| | - Richard H Clayton
- Insigneo Institute for in-silico Medicine and Department of Computer Science, University of Sheffield, UK.
| |
Collapse
|
31
|
Babtie AC, Stumpf MPH. How to deal with parameters for whole-cell modelling. J R Soc Interface 2017; 14:20170237. [PMID: 28768879 PMCID: PMC5582120 DOI: 10.1098/rsif.2017.0237] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2017] [Accepted: 06/22/2017] [Indexed: 11/12/2022] Open
Abstract
Dynamical systems describing whole cells are on the verge of becoming a reality. But as models of reality, they are only useful if we have realistic parameters for the molecular reaction rates and cell physiological processes. There is currently no suitable framework to reliably estimate hundreds, let alone thousands, of reaction rate parameters. Here, we map out the relative weaknesses and promises of different approaches aimed at redressing this issue. While suitable procedures for estimation or inference of the whole (vast) set of parameters will, in all likelihood, remain elusive, some hope can be drawn from the fact that much of the cellular behaviour may be explained in terms of smaller sets of parameters. Identifying such parameter sets and assessing their behaviour is now becoming possible even for very large systems of equations, and we expect such methods to become central tools in the development and analysis of whole-cell models.
Collapse
Affiliation(s)
- Ann C Babtie
- Department of Life Sciences, Imperial College London, London, UK
| | | |
Collapse
|