1
|
Quinn KN, Abbott MC, Transtrum MK, Machta BB, Sethna JP. Information geometry for multiparameter models: new perspectives on the origin of simplicity. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2022; 86:10.1088/1361-6633/aca6f8. [PMID: 36576176 PMCID: PMC10018491 DOI: 10.1088/1361-6633/aca6f8] [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: 11/03/2021] [Accepted: 11/29/2022] [Indexed: 05/20/2023]
Abstract
Complex models in physics, biology, economics, and engineering are oftensloppy, meaning that the model parameters are not well determined by the model predictions for collective behavior. Many parameter combinations can vary over decades without significant changes in the predictions. This review uses information geometry to explore sloppiness and its deep relation to emergent theories. We introduce themodel manifoldof predictions, whose coordinates are the model parameters. Itshyperribbonstructure explains why only a few parameter combinations matter for the behavior. We review recent rigorous results that connect the hierarchy of hyperribbon widths to approximation theory, and to the smoothness of model predictions under changes of the control variables. We discuss recent geodesic methods to find simpler models on nearby boundaries of the model manifold-emergent theories with fewer parameters that explain the behavior equally well. We discuss a Bayesian prior which optimizes the mutual information between model parameters and experimental data, naturally favoring points on the emergent boundary theories and thus simpler models. We introduce a 'projected maximum likelihood' prior that efficiently approximates this optimal prior, and contrast both to the poor behavior of the traditional Jeffreys prior. We discuss the way the renormalization group coarse-graining in statistical mechanics introduces a flow of the model manifold, and connect stiff and sloppy directions along the model manifold with relevant and irrelevant eigendirections of the renormalization group. Finally, we discuss recently developed 'intensive' embedding methods, allowing one to visualize the predictions of arbitrary probabilistic models as low-dimensional projections of an isometric embedding, and illustrate our method by generating the model manifold of the Ising model.
Collapse
Affiliation(s)
- Katherine N Quinn
- Center for the Physics of Biological Function, Princeton University, Princeton, NJ, United States of America
| | - Michael C Abbott
- Department of Physics, Yale University, New Haven, CT, United States of America
| | - Mark K Transtrum
- Department of Physics and Astronomy, Brigham Young University, Provo, UT, United States of America
| | - Benjamin B Machta
- Department of Physics and Systems Biology Institute, Yale University, New Haven, CT, United States of America
| | - James P Sethna
- Department of Physics, Cornell University, Ithaca, NY, United States of America
| |
Collapse
|
2
|
Abstract
The global burden caused by cardiovascular disease is substantial, with heart disease representing the most common cause of death around the world. There remains a need to develop better mechanistic models of cardiac function in order to combat this health concern. Heart rhythm disorders, or arrhythmias, are one particular type of disease which has been amenable to quantitative investigation. Here we review the application of quantitative methodologies to explore dynamical questions pertaining to arrhythmias. We begin by describing single-cell models of cardiac myocytes, from which two and three dimensional models can be constructed. Special focus is placed on results relating to pattern formation across these spatially-distributed systems, especially the formation of spiral waves of activation. Next, we discuss mechanisms which can lead to the initiation of arrhythmias, focusing on the dynamical state of spatially discordant alternans, and outline proposed mechanisms perpetuating arrhythmias such as fibrillation. We then review experimental and clinical results related to the spatio-temporal mapping of heart rhythm disorders. Finally, we describe treatment options for heart rhythm disorders and demonstrate how statistical physics tools can provide insights into the dynamics of heart rhythm disorders.
Collapse
Affiliation(s)
- Wouter-Jan Rappel
- Department of Physics, University of California San Diego, La Jolla, CA 92037
| |
Collapse
|
3
|
Whittaker DG, Wang J, Shuttleworth JG, Venkateshappa R, Kemp JM, Claydon TW, Mirams GR. Ion channel model reduction using manifold boundaries. JOURNAL OF THE ROYAL SOCIETY, INTERFACE 2022; 19:20220193. [PMID: 35946166 PMCID: PMC9363999 DOI: 10.1098/rsif.2022.0193] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Mathematical models of voltage-gated ion channels are used in basic research, industrial and clinical settings. These models range in complexity, but typically contain numerous variables representing the proportion of channels in a given state, and parameters describing the voltage-dependent rates of transition between states. An open problem is selecting the appropriate degree of complexity and structure for an ion channel model given data availability. Here, we simplify a model of the cardiac human Ether-à-go-go related gene (hERG) potassium ion channel, which carries cardiac IKr, using the manifold boundary approximation method (MBAM). The MBAM approximates high-dimensional model-output manifolds by reduced models describing their boundaries, resulting in models with fewer parameters (and often variables). We produced a series of models of reducing complexity starting from an established five-state hERG model with 15 parameters. Models with up to three fewer states and eight fewer parameters were shown to retain much of the predictive capability of the full model and were validated using experimental hERG1a data collected in HEK293 cells at 37°C. The method provides a way to simplify complex models of ion channels that improves parameter identifiability and will aid in future model development.
Collapse
Affiliation(s)
- Dominic G Whittaker
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| | - Jiahui Wang
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| | - Joseph G Shuttleworth
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| | | | - Jacob M Kemp
- Department of Biomedical Physiology and Kinesiology, Simon Fraser University, Burnaby, Canada
| | - Thomas W Claydon
- Department of Biomedical Physiology and Kinesiology, Simon Fraser University, Burnaby, Canada
| | - Gary R Mirams
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| |
Collapse
|
4
|
Karanam A, Rappel WJ. Boolean modelling in plant biology. QUANTITATIVE PLANT BIOLOGY 2022; 3:e29. [PMID: 37077966 PMCID: PMC10095905 DOI: 10.1017/qpb.2022.26] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/17/2022] [Revised: 10/24/2022] [Accepted: 11/16/2022] [Indexed: 05/03/2023]
Abstract
Signalling and genetic networks underlie most biological processes and are often complex, containing many highly connected components. Modelling these networks can provide insight into mechanisms but is challenging given that rate parameters are often not well defined. Boolean modelling, in which components can only take on a binary value with connections encoded by logic equations, is able to circumvent some of these challenges, and has emerged as a viable tool to probe these complex networks. In this review, we will give an overview of Boolean modelling, with a specific emphasis on its use in plant biology. We review how Boolean modelling can be used to describe biological networks and then discuss examples of its applications in plant genetics and plant signalling.
Collapse
Affiliation(s)
- Aravind Karanam
- Department of Physics, University of California, San Diego, La Jolla, California92093, USA
| | - Wouter-Jan Rappel
- Department of Physics, University of California, San Diego, La Jolla, California92093, USA
- Author for correspondence: W.-J. Rappel, E-mail:
| |
Collapse
|
5
|
Simulating Notch-Dome Morphology of Action Potential of Ventricular Cell: How the Speeds of Positive and Negative Feedbacks on Transmembrane Voltage Can Influence the Health of a Cell? BIOMED RESEARCH INTERNATIONAL 2020; 2020:5169241. [PMID: 32953882 PMCID: PMC7487097 DOI: 10.1155/2020/5169241] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/27/2019] [Revised: 03/23/2020] [Accepted: 07/18/2020] [Indexed: 11/24/2022]
Abstract
Ventricular action potential is well-known because of its plateau phase with a spike-notch-dome morphology. As such, the morphology of action potential is necessary for ensuring a correct heart functioning. Any distraction from normal notch-dome morphology may trigger a circus movement reentry in the form of lethal ventricular fibrillation. When the epicardial action potential dome propagates from a site where it is maintained to regions where it has been lost, it gives rise to the proposed mechanism for the Brugada syndrome. Despite the impact of notch-dome dynamics on the heart function, no independent and explicit research has been performed on the simulation of notch-dome dynamics and morphology. In this paper, using a novel mathematical approach, a three-state variable model is proposed; we show that our proposed model not only can simulate morphology of action potential of ventricular cells but also can propose a biological reasonable tool for controlling of the morphology of action potential spike-notch-dome. We show that the processes of activation and inactivation of ionic gating variables (as positive or negative feedbacks on the voltage of cell membrane) and the ratio of their speeds (time constants) can be treated as a reasonable biological tool for simulating ventricular cell notch-dome. This finding may led to a new insight to the quantification of the health of a ventricular cell and may also propose a new drug therapy strategy for cardiac diseases.
Collapse
|
6
|
Whittaker DG, Clerx M, Lei CL, Christini DJ, Mirams GR. Calibration of ionic and cellular cardiac electrophysiology models. WILEY INTERDISCIPLINARY REVIEWS. SYSTEMS BIOLOGY AND MEDICINE 2020; 12:e1482. [PMID: 32084308 PMCID: PMC8614115 DOI: 10.1002/wsbm.1482] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Revised: 01/17/2020] [Accepted: 01/18/2020] [Indexed: 12/30/2022]
Abstract
Cardiac electrophysiology models are among the most mature and well-studied mathematical models of biological systems. This maturity is bringing new challenges as models are being used increasingly to make quantitative rather than qualitative predictions. As such, calibrating the parameters within ion current and action potential (AP) models to experimental data sets is a crucial step in constructing a predictive model. This review highlights some of the fundamental concepts in cardiac model calibration and is intended to be readily understood by computational and mathematical modelers working in other fields of biology. We discuss the classic and latest approaches to calibration in the electrophysiology field, at both the ion channel and cellular AP scales. We end with a discussion of the many challenges that work to date has raised and the need for reproducible descriptions of the calibration process to enable models to be recalibrated to new data sets and built upon for new studies. This article is categorized under: Analytical and Computational Methods > Computational Methods Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Cellular Models.
Collapse
Affiliation(s)
- Dominic G. Whittaker
- Centre for Mathematical Medicine & Biology, School of Mathematical SciencesUniversity of NottinghamNottinghamUK
| | - Michael Clerx
- Computational Biology & Health Informatics, Department of Computer ScienceUniversity of OxfordOxfordUK
| | - Chon Lok Lei
- Computational Biology & Health Informatics, Department of Computer ScienceUniversity of OxfordOxfordUK
| | | | - Gary R. Mirams
- Centre for Mathematical Medicine & Biology, School of Mathematical SciencesUniversity of NottinghamNottinghamUK
| |
Collapse
|
7
|
Houston C, Marchand B, Engelbert L, Cantwell CD. Reducing complexity and unidentifiability when modelling human atrial cells. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2020. [PMID: 32448063 DOI: 10.5061/dryad.p2ngf1vmc] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
Mathematical models of a cellular action potential (AP) in cardiac modelling have become increasingly complex, particularly in gating kinetics, which control the opening and closing of individual ion channel currents. As cardiac models advance towards use in personalized medicine to inform clinical decision-making, it is critical to understand the uncertainty hidden in parameter estimates from their calibration to experimental data. This study applies approximate Bayesian computation to re-calibrate the gating kinetics of four ion channels in two existing human atrial cell models to their original datasets, providing a measure of uncertainty and indication of potential issues with selecting a single unique value given the available experimental data. Two approaches are investigated to reduce the uncertainty present: re-calibrating the models to a more complete dataset and using a less complex formulation with fewer parameters to constrain. The re-calibrated models are inserted back into the full cell model to study the overall effect on the AP. The use of more complete datasets does not eliminate uncertainty present in parameter estimates. The less complex model, particularly for the fast sodium current, gave a better fit to experimental data alongside lower parameter uncertainty and improved computational speed. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.
Collapse
Affiliation(s)
- C Houston
- ElectroCardioMaths Programme, Centre for Cardiac Engineering, Imperial College, London, UK
- Department of Aeronautics, Imperial College, London, UK
| | - B Marchand
- Department of Aeronautics, Imperial College, London, UK
| | - L Engelbert
- Department of Aeronautics, Imperial College, London, UK
| | - C D Cantwell
- ElectroCardioMaths Programme, Centre for Cardiac Engineering, Imperial College, London, UK
- Department of Aeronautics, Imperial College, London, UK
| |
Collapse
|
8
|
Houston C, Marchand B, Engelbert L, Cantwell CD. Reducing complexity and unidentifiability when modelling human atrial cells. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2020; 378:20190339. [PMID: 32448063 PMCID: PMC7287336 DOI: 10.1098/rsta.2019.0339] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/07/2023]
Abstract
Mathematical models of a cellular action potential (AP) in cardiac modelling have become increasingly complex, particularly in gating kinetics, which control the opening and closing of individual ion channel currents. As cardiac models advance towards use in personalized medicine to inform clinical decision-making, it is critical to understand the uncertainty hidden in parameter estimates from their calibration to experimental data. This study applies approximate Bayesian computation to re-calibrate the gating kinetics of four ion channels in two existing human atrial cell models to their original datasets, providing a measure of uncertainty and indication of potential issues with selecting a single unique value given the available experimental data. Two approaches are investigated to reduce the uncertainty present: re-calibrating the models to a more complete dataset and using a less complex formulation with fewer parameters to constrain. The re-calibrated models are inserted back into the full cell model to study the overall effect on the AP. The use of more complete datasets does not eliminate uncertainty present in parameter estimates. The less complex model, particularly for the fast sodium current, gave a better fit to experimental data alongside lower parameter uncertainty and improved computational speed. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.
Collapse
Affiliation(s)
- C. Houston
- ElectroCardioMaths Programme, Centre for Cardiac Engineering, Imperial College, London, UK
- Department of Aeronautics, Imperial College, London, UK
- e-mail:
| | - B. Marchand
- Department of Aeronautics, Imperial College, London, UK
| | - L. Engelbert
- Department of Aeronautics, Imperial College, London, UK
| | - C. D. Cantwell
- ElectroCardioMaths Programme, Centre for Cardiac Engineering, Imperial College, London, UK
- Department of Aeronautics, Imperial College, London, UK
| |
Collapse
|
9
|
Zhang J, De-Oliveira-Ceciliato P, Takahashi Y, Schulze S, Dubeaux G, Hauser F, Azoulay-Shemer T, Tõldsepp K, Kollist H, Rappel WJ, Schroeder JI. Insights into the Molecular Mechanisms of CO 2-Mediated Regulation of Stomatal Movements. Curr Biol 2019; 28:R1356-R1363. [PMID: 30513335 DOI: 10.1016/j.cub.2018.10.015] [Citation(s) in RCA: 45] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Plants must continually balance the influx of CO2 for photosynthesis against the loss of water vapor through stomatal pores in their leaves. This balance can be achieved by controlling the aperture of the stomatal pores in response to several environmental stimuli. Elevation in atmospheric [CO2] induces stomatal closure and further impacts leaf temperatures, plant growth and water-use efficiency, and global crop productivity. Here, we review recent advances in understanding CO2-perception mechanisms and CO2-mediated signal transduction in the regulation of stomatal movements, and we explore how these mechanisms are integrated with other signaling pathways in guard cells.
Collapse
Affiliation(s)
- Jingbo Zhang
- Division of Biological Sciences, Cell and Developmental Biology Section, University of California San Diego, La Jolla, CA 92093, USA
| | - Paulo De-Oliveira-Ceciliato
- Division of Biological Sciences, Cell and Developmental Biology Section, University of California San Diego, La Jolla, CA 92093, USA
| | - Yohei Takahashi
- Division of Biological Sciences, Cell and Developmental Biology Section, University of California San Diego, La Jolla, CA 92093, USA
| | - Sebastian Schulze
- Division of Biological Sciences, Cell and Developmental Biology Section, University of California San Diego, La Jolla, CA 92093, USA
| | - Guillaume Dubeaux
- Division of Biological Sciences, Cell and Developmental Biology Section, University of California San Diego, La Jolla, CA 92093, USA
| | - Felix Hauser
- Division of Biological Sciences, Cell and Developmental Biology Section, University of California San Diego, La Jolla, CA 92093, USA
| | - Tamar Azoulay-Shemer
- Division of Biological Sciences, Cell and Developmental Biology Section, University of California San Diego, La Jolla, CA 92093, USA
| | - Kadri Tõldsepp
- Institute of Technology, University of Tartu, Tartu 50411, Estonia
| | - Hannes Kollist
- Institute of Technology, University of Tartu, Tartu 50411, Estonia
| | - Wouter-Jan Rappel
- Physics Department, University of California San Diego, La Jolla, CA 92093, USA
| | - Julian I Schroeder
- Division of Biological Sciences, Cell and Developmental Biology Section, University of California San Diego, La Jolla, CA 92093, USA.
| |
Collapse
|
10
|
Lombardo DM, Rappel WJ. Chaotic tip trajectories of a single spiral wave in the presence of heterogeneities. Phys Rev E 2019; 99:062409. [PMID: 31330597 DOI: 10.1103/physreve.99.062409] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2017] [Indexed: 11/07/2022]
Abstract
Spiral waves have been observed in a variety of physical, chemical, and biological systems. They play a major role in cardiac arrhythmias, including fibrillation, where the observed irregular activation patterns are generally thought to arise from the continuous breakup of multiple unstable spiral waves. Using spatially extended simulations of different electrophysiological models of cardiac tissue, we show that a single spiral wave in the presence of heterogeneities can display chaotic tip trajectories, consistent with fibrillation. We also show that the simulated spiral tip dynamics, including chaotic trajectories, can be captured by a simple particle model which only describes the dynamics of the spiral tip. This shows that spiral wave breakup, or interactions with other waves, are not necessary to initiate chaos in spiral waves.
Collapse
Affiliation(s)
- Daniel M Lombardo
- Department of Physics, University of California San Diego, San Diego, California 92093, USA
| | - Wouter-Jan Rappel
- Department of Physics, University of California San Diego, San Diego, California 92093, USA
| |
Collapse
|
11
|
Jæger KH, Wall S, Tveito A. Detecting undetectables: Can conductances of action potential models be changed without appreciable change in the transmembrane potential? CHAOS (WOODBURY, N.Y.) 2019; 29:073102. [PMID: 31370420 DOI: 10.1063/1.5087629] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/03/2019] [Accepted: 06/12/2019] [Indexed: 05/23/2023]
Abstract
Mathematical models describing the dynamics of the cardiac action potential are of great value for understanding how changes to the system can disrupt the normal electrical activity of cells and tissue in the heart. However, to represent specific data, these models must be parameterized, and adjustment of the maximum conductances of the individual contributing ionic currents is a commonly used method. Here, we present a method for investigating the uniqueness of such resulting parameterizations. Our key question is: Can the maximum conductances of a model be changed without giving any appreciable changes in the action potential? If so, the model parameters are not unique and this poses a major problem in using the models to identify changes in parameters from data, for instance, to evaluate potential drug effects. We propose a method for evaluating this uniqueness, founded on the singular value decomposition of a matrix consisting of the individual ionic currents. Small singular values of this matrix signify lack of parameter uniqueness and we show that the conclusion from linear analysis of the matrix carries over to provide insight into the uniqueness of the parameters in the nonlinear case. Using numerical experiments, we quantify the identifiability of the maximum conductances of well-known models of the cardiac action potential. Furthermore, we show how the identifiability depends on the time step used in the observation of the currents, how the application of drugs may change identifiability, and, finally, how the stimulation protocol can be used to improve the identifiability of a model.
Collapse
Affiliation(s)
| | - Samuel Wall
- Simula Research Laboratory, 1325 Lysaker, Norway
| | - Aslak Tveito
- Simula Research Laboratory, 1325 Lysaker, Norway
| |
Collapse
|