1
|
Ni H, Morotti S, Grandi E. A Heart for Diversity: Simulating Variability in Cardiac Arrhythmia Research. Front Physiol 2018; 9:958. [PMID: 30079031 PMCID: PMC6062641 DOI: 10.3389/fphys.2018.00958] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2018] [Accepted: 06/29/2018] [Indexed: 12/31/2022] Open
Abstract
In cardiac electrophysiology, there exist many sources of inter- and intra-personal variability. These include variability in conditions and environment, and genotypic and molecular diversity, including differences in expression and behavior of ion channels and transporters, which lead to phenotypic diversity (e.g., variable integrated responses at the cell, tissue, and organ levels). These variabilities play an important role in progression of heart disease and arrhythmia syndromes and outcomes of therapeutic interventions. Yet, the traditional in silico framework for investigating cardiac arrhythmias is built upon a parameter/property-averaging approach that typically overlooks the physiological diversity. Inspired by work done in genetics and neuroscience, new modeling frameworks of cardiac electrophysiology have been recently developed that take advantage of modern computational capabilities and approaches, and account for the variance in the biological data they are intended to illuminate. In this review, we outline the recent advances in statistical and computational techniques that take into account physiological variability, and move beyond the traditional cardiac model-building scheme that involves averaging over samples from many individuals in the construction of a highly tuned composite model. We discuss how these advanced methods have harnessed the power of big (simulated) data to study the mechanisms of cardiac arrhythmias, with a special emphasis on atrial fibrillation, and improve the assessment of proarrhythmic risk and drug response. The challenges of using in silico approaches with variability are also addressed and future directions are proposed.
Collapse
Affiliation(s)
| | | | - Eleonora Grandi
- Department of Pharmacology, University of California, Davis, Davis, CA, United States
| |
Collapse
|
2
|
Images as drivers of progress in cardiac computational modelling. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2014; 115:198-212. [PMID: 25117497 PMCID: PMC4210662 DOI: 10.1016/j.pbiomolbio.2014.08.005] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/24/2014] [Accepted: 08/02/2014] [Indexed: 11/28/2022]
Abstract
Computational models have become a fundamental tool in cardiac research. Models are evolving to cover multiple scales and physical mechanisms. They are moving towards mechanistic descriptions of personalised structure and function, including effects of natural variability. These developments are underpinned to a large extent by advances in imaging technologies. This article reviews how novel imaging technologies, or the innovative use and extension of established ones, integrate with computational models and drive novel insights into cardiac biophysics. In terms of structural characterization, we discuss how imaging is allowing a wide range of scales to be considered, from cellular levels to whole organs. We analyse how the evolution from structural to functional imaging is opening new avenues for computational models, and in this respect we review methods for measurement of electrical activity, mechanics and flow. Finally, we consider ways in which combined imaging and modelling research is likely to continue advancing cardiac research, and identify some of the main challenges that remain to be solved.
Collapse
|
3
|
Roberts BN, Yang PC, Behrens SB, Moreno JD, Clancy CE. Computational approaches to understand cardiac electrophysiology and arrhythmias. Am J Physiol Heart Circ Physiol 2012; 303:H766-83. [PMID: 22886409 DOI: 10.1152/ajpheart.01081.2011] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Cardiac rhythms arise from electrical activity generated by precisely timed opening and closing of ion channels in individual cardiac myocytes. These impulses spread throughout the cardiac muscle to manifest as electrical waves in the whole heart. Regularity of electrical waves is critically important since they signal the heart muscle to contract, driving the primary function of the heart to act as a pump and deliver blood to the brain and vital organs. When electrical activity goes awry during a cardiac arrhythmia, the pump does not function, the brain does not receive oxygenated blood, and death ensues. For more than 50 years, mathematically based models of cardiac electrical activity have been used to improve understanding of basic mechanisms of normal and abnormal cardiac electrical function. Computer-based modeling approaches to understand cardiac activity are uniquely helpful because they allow for distillation of complex emergent behaviors into the key contributing components underlying them. Here we review the latest advances and novel concepts in the field as they relate to understanding the complex interplay between electrical, mechanical, structural, and genetic mechanisms during arrhythmia development at the level of ion channels, cells, and tissues. We also discuss the latest computational approaches to guiding arrhythmia therapy.
Collapse
Affiliation(s)
- Byron N Roberts
- Tri-Institutional MD-PhD Program, Physiology, Biophysics and Systems Biology Graduate Program, Weill Cornell Medical College/The Rockefeller University/Sloan-Kettering Cancer Institute, Weill Medical College of Cornell University, New York, New York, USA
| | | | | | | | | |
Collapse
|
4
|
Reumann M, Fitch BG, Rayshubskiy A, Pitman MC, Rice JJ. Orthogonal recursive bisection as data decomposition strategy for massively parallel cardiac simulations. BIOMED ENG-BIOMED TE 2011; 56:129-45. [PMID: 21657987 DOI: 10.1515/bmt.2011.100] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
We present the orthogonal recursive bisection algorithm that hierarchically segments the anatomical model structure into subvolumes that are distributed to cores. The anatomy is derived from the Visible Human Project, with electrophysiology based on the FitzHugh-Nagumo (FHN) and ten Tusscher (TT04) models with monodomain diffusion. Benchmark simulations with up to 16,384 and 32,768 cores on IBM Blue Gene/P and L supercomputers for both FHN and TT04 results show good load balancing with almost perfect speedup factors that are close to linear with the number of cores. Hence, strong scaling is demonstrated. With 32,768 cores, a 1000 ms simulation of full heart beat requires about 6.5 min of wall clock time for a simulation of the FHN model. For the largest machine partitions, the simulations execute at a rate of 0.548 s (BG/P) and 0.394 s (BG/L) of wall clock time per 1 ms of simulation time. To our knowledge, these simulations show strong scaling to substantially higher numbers of cores than reported previously for organ-level simulation of the heart, thus significantly reducing run times. The ability to reduce runtimes could play a critical role in enabling wider use of cardiac models in research and clinical applications.
Collapse
Affiliation(s)
- Matthias Reumann
- IBM T. J. Watson Research Center, Yorktown Heights, NY 10598, USA.
| | | | | | | | | |
Collapse
|
5
|
Reumann M, Fitch BG, Rayshubskiy A, Keller DUJ, Seemann G, Dossel O, Pitman MC, Rice JJ. Orthogonal recursive bisection data decomposition for high performance computing in cardiac model simulations: dependence on anatomical geometry. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2009; 2009:2799-2802. [PMID: 19964263 DOI: 10.1109/iembs.2009.5333803] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
Orthogonal recursive bisection (ORB) algorithm can be used as data decomposition strategy to distribute a large data set of a cardiac model to a distributed memory supercomputer. It has been shown previously that good scaling results can be achieved using the ORB algorithm for data decomposition. However, the ORB algorithm depends on the distribution of computational load of each element in the data set. In this work we investigated the dependence of data decomposition and load balancing on different rotations of the anatomical data set to achieve optimization in load balancing. The anatomical data set was given by both ventricles of the Visible Female data set in a 0.2 mm resolution. Fiber orientation was included. The data set was rotated by 90 degrees around x, y and z axis, respectively. By either translating or by simply taking the magnitude of the resulting negative coordinates we were able to create 14 data set of the same anatomy with different orientation and position in the overall volume. Computation load ratios for non - tissue vs. tissue elements used in the data decomposition were 1:1, 1:2, 1:5, 1:10, 1:25, 1:38.85, 1:50 and 1:100 to investigate the effect of different load ratios on the data decomposition. The ten Tusscher et al. (2004) electrophysiological cell model was used in monodomain simulations of 1 ms simulation time to compare performance using the different data sets and orientations. The simulations were carried out for load ratio 1:10, 1:25 and 1:38.85 on a 512 processor partition of the IBM Blue Gene/L supercomputer. Th results show that the data decomposition does depend on the orientation and position of the anatomy in the global volume. The difference in total run time between the data sets is 10 s for a simulation time of 1 ms. This yields a difference of about 28 h for a simulation of 10 s simulation time. However, given larger processor partitions, the difference in run time decreases and becomes less significant. Depending on the processor partition size, future work will have to consider the orientation of the anatomy in the global volume for longer simulation runs.
Collapse
Affiliation(s)
- Matthias Reumann
- Computational Biology Center, IBM TJ Watson Research Center, Yorktown Heights, 1101 Kitchawan Road, Route 134, NY 10598, USA.
| | | | | | | | | | | | | | | |
Collapse
|
6
|
Reumann M, Fitch BG, Rayshubskiy A, Keller DUJ, Seemann G, Dossel O, Pitman MC, Rice JJ. Strong scaling and speedup to 16,384 processors in cardiac electro-mechanical simulations. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2009; 2009:2795-8. [PMID: 19964262 DOI: 10.1109/iembs.2009.5333802] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
High performance computing is required to make feasible simulations of whole organ models of the heart with biophysically detailed cellular models in a clinical setting. Increasing model detail by simulating electrophysiology and mechanical models increases computation demands. We present scaling results of an electro - mechanical cardiac model of two ventricles and compare them to our previously published results using an electrophysiological model only. The anatomical data-set was given by both ventricles of the Visible Female data-set in a 0.2 mm resolution. Fiber orientation was included. Data decomposition for the distribution onto the distributed memory system was carried out by orthogonal recursive bisection. Load weight ratios for non-tissue vs. tissue elements used in the data decomposition were 1:1, 1:2, 1:5, 1:10, 1:25, 1:38.85, 1:50 and 1:100. The ten Tusscher et al. (2004) electrophysiological cell model was used and the Rice et al. (1999) model for the computation of the calcium transient dependent force. Scaling results for 512, 1024, 2048, 4096, 8192 and 16,384 processors were obtained for 1 ms simulation time. The simulations were carried out on an IBM Blue Gene/L supercomputer. The results show linear scaling from 512 to 16,384 processors with speedup factors between 1.82 and 2.14 between partitions. The most optimal load ratio was 1:25 for on all partitions. However, a shift towards load ratios with higher weight for the tissue elements can be recognized as can be expected when adding computational complexity to the model while keeping the same communication setup. This work demonstrates that it is potentially possible to run simulations of 0.5 s using the presented electro-mechanical cardiac model within 1.5 hours.
Collapse
Affiliation(s)
- Matthias Reumann
- Computational Biology Center, IBM TJ Watson Research Center, Yorktown Heights, 1101 Kitchawan Road, Route 134, NY 10598, USA.
| | | | | | | | | | | | | | | |
Collapse
|
7
|
Reumann M, Gurev V, Rice JJ. Computational modeling of cardiac disease: potential for personalized medicine. Per Med 2009; 6:45-66. [DOI: 10.2217/17410541.6.1.45] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Cardiovascular diseases are leading causes of death, reduce life quality and consume almost half a trillion dollars in healthcare expenses in the USA alone. Cardiac modeling and simulation technologies hold promise as important tools to improve cardiac care and are already in use to elucidate the fundamental mechanisms of cardiac physiology and pathophysiology. However, the emphasis has been on simulating average or exemplar cases. This report describes two classes of cardiac modeling efforts. First, electrophysiological models of channelopathies simulate the altered electrical activity that is thought to promote arrhythmias. Second, electromechanical models attempt to capture both the electrophysiological and mechanical aspects of heart function. One goal of the community is to develop models with sufficient patient customization to assist in personalized treatment planning. Some model aspects can be customized with existing data collection techniques to more closely represent individual patients while other model aspects will likely remain based on generic data. Despite important challenges, cardiac models hold promise to be important enablers of personalized medicine.
Collapse
Affiliation(s)
- Matthias Reumann
- Functional Genomics and Systems Biology, IBM T.J. Watson Research Center, PO Box 218, Yorktown Heights, NY 10598, USA
| | - Viatcheslav Gurev
- Department of Biomedical Engineering and Institute for Computational Medicine, The Johns Hopkins University, MD, USA
| | - John Jeremy Rice
- Functional Genomics and Systems Biology, IBM T.J. Watson Research Center, PO Box 218, Yorktown Heights, NY 10598, USA
- Department of Biomedical Engineering and Institute for Computational Medicine, The Johns Hopkins University, MD, USA
| |
Collapse
|
8
|
Fenner JW, Brook B, Clapworthy G, Coveney PV, Feipel V, Gregersen H, Hose DR, Kohl P, Lawford P, McCormack KM, Pinney D, Thomas SR, Van Sint Jan S, Waters S, Viceconti M. The EuroPhysiome, STEP and a roadmap for the virtual physiological human. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2008; 366:2979-99. [PMID: 18559316 DOI: 10.1098/rsta.2008.0089] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
Biomedical science and its allied disciplines are entering a new era in which computational methods and technologies are poised to play a prevalent role in supporting collaborative investigation of the human body. Within Europe, this has its focus in the virtual physiological human (VPH), which is an evolving entity that has emerged from the EuroPhysiome initiative and the strategy for the EuroPhysiome (STEP) consortium. The VPH is intended to be a solution to common infrastructure needs for physiome projects across the globe, providing a unifying architecture that facilitates integration and prediction, ultimately creating a framework capable of describing Homo sapiens in silico. The routine reliance of the biomedical industry, biomedical research and clinical practice on information technology (IT) highlights the importance of a tailor-made and robust IT infrastructure, but numerous challenges need to be addressed if the VPH is to become a mature technological reality. Appropriate investment will reap considerable rewards, since it is anticipated that the VPH will influence all sectors of society, with implications predominantly for improved healthcare, improved competitiveness in industry and greater understanding of (patho)physiological processes. This paper considers issues pertinent to the development of the VPH, highlighted by the work of the STEP consortium.
Collapse
Affiliation(s)
- J W Fenner
- Department of Medical Physics and Clinical Engineering, University of Sheffield, I Floor, Royal Hallamshire Hospital, Glossop Road, Sheffield S10 2JF, UK.
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
9
|
Clayton RH, Panfilov AV. A guide to modelling cardiac electrical activity in anatomically detailed ventricles. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2007; 96:19-43. [PMID: 17825362 DOI: 10.1016/j.pbiomolbio.2007.07.004] [Citation(s) in RCA: 88] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
One of the most recent trends in cardiac electrophysiology is the development of integrative anatomically accurate models of the heart, which include description of cardiac activity from sub-cellular and cellular level to the level of the whole organ. In order to construct this type of model, a researcher needs to collect a wide range of information from books and journal articles on various aspects of biology, physiology, electrophysiology, numerical mathematics and computer programming. The aim of this methodological article is to survey recent developments in integrative modelling of electrical activity in the ventricles of the heart, and to provide a practical guide to the resources and tools that are available for work in this exciting and challenging area.
Collapse
Affiliation(s)
- R H Clayton
- Department of Computer Science, University of Sheffield, Regent Court, 211 Portobello Street, Sheffield, S1 4DP, UK.
| | | |
Collapse
|
10
|
Burton RAB, Plank G, Schneider JE, Grau V, Ahammer H, Keeling SL, Lee J, Smith NP, Gavaghan D, Trayanova N, Kohl P. Three-dimensional models of individual cardiac histoanatomy: tools and challenges. Ann N Y Acad Sci 2007; 1080:301-19. [PMID: 17132791 PMCID: PMC3313659 DOI: 10.1196/annals.1380.023] [Citation(s) in RCA: 81] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
There is a need for, and utility in, the acquisition of data sets of cardiac histoanatomy, with the vision of reconstructing individual hearts on the basis of noninvasive imaging, such as MRI, enriched by reference to detailed atlases of serial histology obtained from representative samples. These data sets would be useful not only as a repository of knowledge regarding the specifics of cardiac histoanatomy, but could form the basis for generation of individualized high-resolution cardiac structure-function models. The current article presents a step in this general direction: it illustrates how whole-heart noninvasive imaging can be combined with whole-heart histology in an approach to achieve automated construction of histoanatomically detailed models of cardiac 3D structure and function at hitherto unprecedented resolution and accuracy (based on 26.4 x 26.4 x 24.4 microm MRI voxel size, and enriched by histological detail). It provides an overview of the tools used in this quest and outlines challenges posed by the approach in the light of applications that may benefit from the availability of such data and tools.
Collapse
Affiliation(s)
- Rebecca A B Burton
- Department of Physiology, Anatomy and Genetics, University of Oxford, Oxford OX1 3PT, UK
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
11
|
Welsh E, Jirotka M, Gavaghan D. Post-genomic science: cross-disciplinary and large-scale collaborative research and its organizational and technological challenges for the scientific research process. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2006; 364:1533-49. [PMID: 16766359 DOI: 10.1098/rsta.2006.1785] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
We examine recent developments in cross-disciplinary science and contend that a 'Big Science' approach is increasingly evident in the life sciences-facilitated by a breakdown of the traditional barriers between academic disciplines and the application of technologies across these disciplines. The first fruits of 'Big Biology' are beginning to be seen in, for example, genomics, (bio)-nanotechnology and systems biology. We suggest that this has profound implications for the research process and presents challenges both in technological design, in the provision of infrastructure and training, in the organization of research groups, and in providing suitable research funding mechanisms and reward systems. These challenges need to be addressed if the promise of this approach is to be fully realized. In this paper, we will draw on the work of social scientists to understand how these developments in science and technology relate to organizational culture, organizational change and the context of scientific work. We seek to learn from previous technological developments that seemed to offer similar potential for organizational and social change.
Collapse
Affiliation(s)
- Elaine Welsh
- Department of Sociology, University of Warwick, Coventry CV4 7AL, UK.
| | | | | |
Collapse
|
12
|
Seemann G, Höper C, Sachse FB, Dössel O, Holden AV, Zhang H. Heterogeneous three-dimensional anatomical and electrophysiological model of human atria. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2006; 364:1465-81. [PMID: 16766355 DOI: 10.1098/rsta.2006.1781] [Citation(s) in RCA: 134] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
Investigating the mechanisms underlying the genesis and conduction of electrical excitation in the atria at physiological and pathological states is of great importance. To provide knowledge concerning the mechanisms of excitation, we constructed a biophysical detailed and anatomically accurate computer model of human atria that incorporates both structural and electrophysiological heterogeneities. The three-dimensional geometry was extracted from the visible female dataset. The sinoatrial node (SAN) and atrium, including crista terminalis (CT), pectinate muscles (PM), appendages (APG) and Bachmann's bundle (BB) were segmented in this work. Fibre orientation in CT, PM and BB was set to local longitudinal direction. Descriptions for all used cell types were based on modifications of the Courtemanche et al. model of a human atrial cell. Maximum conductances of Ito, IKr and ICa,L were modified for PM, CT, APG and atrioventricular ring to reproduce measured action potentials (AP). Pacemaker activity in the human SAN was reproduced by removing IK1, but including If, ICa,T, and gradients of channel conductances as described in previous studies for heterogeneous rabbit SAN. Anisotropic conduction was computed with a monodomain model using the finite element method. The transversal to longitudinal ratio of conductivity for PM, CT and BB was 1:9. Atrial working myocardium (AWM) was set to be isotropic. Simulation of atrial electrophysiology showed initiation of APs in the SAN centre. The excitation spread afterwards to the periphery near to the region of the CT and preferentially towards the atrioventricular region. The excitation extends over the right atrium along PM. Both CT and PM activated the right AWM. Earliest activation of the left atrium was through BB and excitation spread over to the APG. The conduction velocities were 0.6ms-1 for AWM, 1.2ms-1 for CT, 1.6ms-1 for PM and 1.1ms-1 for BB at a rate of 63bpm. The simulations revealed that bundles form dominant pathways for atrial conduction. The preferential conduction towards CT and along PM is comparable with clinical mapping. Repolarization is more homogeneous than excitation due to the heterogeneous distribution of electrophysiological properties and hence the action potential duration.
Collapse
Affiliation(s)
- Gunnar Seemann
- Institute of Biomedical Engineering, University Karlsruhe (TH), Kaiserstrasse 12, 76128 Karlsruhe, Germany.
| | | | | | | | | | | |
Collapse
|
13
|
Shim EB, Leem CH, Abe Y, Noma A. A new multi-scale simulation model of the circulation: from cells to system. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2006; 364:1483-500. [PMID: 16766356 DOI: 10.1098/rsta.2006.1782] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
We developed a comprehensive cell model that simulates the sequential cellular events from membrane excitation to contraction in the human ventricle. By combining this ventricular cell model with a lumped circulation model, we examined how blood pressure dynamics in the ventricle and aorta are related to the cellular processes. To convert cell contraction into ventricular pressure using Laplace's law, we introduced a simple geometric model of a ventricle: one shaped like a thin-walled hemisphere. The force of contraction of a single cell induces tension in the hemispheric ventricular wall, which generates the ventricular and aortic pressures in the lumped circulation model. The time courses of the hemodynamic properties, as well as the volume-pressure trajectory of the left ventricle, were well reproduced. Our multi-scale cardiovascular model, which covers from cardiac cells to the circulatory system, simulates the typical characteristics of heart mechanics, such as the pressure-volume relationship, stroke volume and the effect of the increased maximum free calcium concentration on cardiovascular hemodynamics. To test the cell-circulation coupling characteristics of the model, we simulated the effects of a decrease in L-type calcium channel conductance (cell level) on left ventricular pressure (system level). The variation due to different pacing frequencies for myocyte excitation was also investigated to assess the effects of heart rate on cardiac cells and the circulatory system.
Collapse
Affiliation(s)
- Eun Bo Shim
- Department of Mechanical and Biomedical Engineering, Kangwon National University, Hyoja-dong, Chuncheon, Kangwon 200-701, South Korea
| | | | | | | |
Collapse
|
14
|
Gavaghan D, Garny A, Maini PK, Kohl P. Mathematical models in physiology. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2006; 364:1099-106. [PMID: 16608698 DOI: 10.1098/rsta.2006.1757] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
Computational modelling of biological processes and systems has witnessed a remarkable development in recent years. The search-term (modelling OR modeling) yields over 58000 entries in PubMed, with more than 34000 since the year 2000: thus, almost two-thirds of papers appeared in the last 5-6 years, compared to only about one-third in the preceding 5-6 decades. The development is fuelled both by the continuously improving tools and techniques available for bio-mathematical modelling and by the increasing demand in quantitative assessment of element inter-relations in complex biological systems. This has given rise to a worldwide public domain effort to build a computational framework that provides a comprehensive theoretical representation of integrated biological function-the Physiome. The current and next issues of this journal are devoted to a small sub-set of this initiative and address biocomputation and modelling in physiology, illustrating the breadth and depth of experimental data-based model development in biological research from sub-cellular events to whole organ simulations.
Collapse
Affiliation(s)
- David Gavaghan
- Oxford University Computing Laboratory, Wolfson Building, Parks Road, Oxford OX1 3QD, UK.
| | | | | | | |
Collapse
|