Multi-scale computational modelling in biology and physiology

Abstract Recent advances in biotechnology and the availability of ever more powerful computers have led to the formulation of increasingly complex models at all levels of biology. One of the main aims of systems biology is to couple these together to produce integrated models across multiple spatial scales and physical processes. In this review, we formulate a definition of multi-scale in terms of levels of biological organisation and describe the types of model that are found at each level. Key issues that arise in trying to formulate and solve multi-scale and multi-physics models are considered and examples of how these issues have been addressed are given for two of the more mature fields in computational biology: the molecular dynamics of ion channels and cardiac modelling. As even more complex models are developed over the coming few years, it will be necessary to develop new methods to model them (in particular in coupling across the interface between stochastic and deterministic processes) and new techniques will be required to compute their solutions efficiently on massively parallel computers. We outline how we envisage these developments occurring.

[1]  F. Yin,et al.  A multiaxial constitutive law for mammalian left ventricular myocardium in steady-state barium contracture or tetanus. , 1998, Journal of biomechanical engineering.

[2]  G. Voth,et al.  Hybrid Ab-Initio/Empirical Molecular Dynamics: Combining the ONIOM Scheme with the Atom-Centered Density Matrix Propagation (ADMP) Approach , 2004 .

[3]  D DiFrancesco,et al.  A model of cardiac electrical activity incorporating ionic pumps and concentration changes. , 1985, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[4]  P. J. Hunter,et al.  Generation of an Anatomically Based Geometric Coronary Model , 2004, Annals of Biomedical Engineering.

[5]  Elizabeth M Cherry,et al.  Suppression of alternans and conduction blocks despite steep APD restitution: electrotonic, memory, and conduction velocity restitution effects. , 2004, American journal of physiology. Heart and circulatory physiology.

[6]  G. W. Beeler,et al.  Reconstruction of the action potential of ventricular myocardial fibres , 1977, The Journal of physiology.

[7]  N Trayanova,et al.  Termination of Spiral Waves with Biphasic Shocks: , 2000, Journal of cardiovascular electrophysiology.

[8]  José Jalife,et al.  Dynamics of rotating vortices in the Beeler-Reuter model of cardiac tissue , 1995 .

[9]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[10]  C. Fraser,et al.  Transmission Dynamics of the Etiological Agent of SARS in Hong Kong: Impact of Public Health Interventions , 2003, Science.

[11]  Steve McKeever,et al.  On the application of partial evaluation to the optimisation of cardiac electrophysiological simulations , 2006, PEPM '06.

[12]  D. Noble Cardiac Action and Pacemaker Potentials based on the Hodgkin-Huxley Equations , 1960, Nature.

[13]  D Peter Tieleman,et al.  COMPUTER SIMULATIONS OF TRANSPORT THROUGH MEMBRANES: PASSIVE DIFFUSION, PORES, CHANNELS AND TRANSPORTERS , 2006, Clinical and experimental pharmacology & physiology.

[14]  Ivet Bahar,et al.  Common mechanism of pore opening shared by five different potassium channels. , 2006, Biophysical journal.

[15]  R. Winslow,et al.  Electrophysiological modeling of cardiac ventricular function: from cell to organ. , 2000, Annual review of biomedical engineering.

[16]  Andrew J. Pullan,et al.  An Anatomically Based Model of Transient Coronary Blood Flow in the Heart , 2002, SIAM J. Appl. Math..

[17]  Catherine M Lloyd,et al.  CellML: its future, present and past. , 2004, Progress in biophysics and molecular biology.

[18]  D. Noble,et al.  Improved guinea-pig ventricular cell model incorporating a diadic space, IKr and IKs, and length- and tension-dependent processes. , 1998, The Canadian journal of cardiology.

[19]  C. Luo,et al.  A model of the ventricular cardiac action potential. Depolarization, repolarization, and their interaction. , 1991, Circulation research.

[20]  Robert Plonsey,et al.  A planar slab bidomain model for cardiac tissue , 2006, Annals of Biomedical Engineering.

[21]  Y. Rudy,et al.  Ionic Current Basis of Electrocardiographic Waveforms: A Model Study , 2002, Circulation research.

[22]  G. Neufeld,et al.  Numerical and experimental study of steady-state CO2 and inert gas washout. , 1988, Journal of applied physiology.

[23]  Jonathan P. Whiteley,et al.  An Efficient Numerical Technique for the Solution of the Monodomain and Bidomain Equations , 2006, IEEE Transactions on Biomedical Engineering.

[24]  C. Hahn,et al.  The effects of ventilation pattern on carbon dioxide transfer in three computer models of the airways , 2002, Respiratory Physiology & Neurobiology.

[25]  Denis Noble,et al.  Models of cardiac ventricular action potentials: iterative interaction between experiment and simulation , 2001, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[26]  K B Campbell,et al.  Nonlinear myofilament regulatory processes affect frequency-dependent muscle fiber stiffness. , 2001, Biophysical journal.

[27]  Edmund J. Crampin,et al.  Multiscale computational modelling of the heart , 2004, Acta Numerica.

[28]  Andrew D. McCulloch,et al.  Effect of Laminar Orthotropic Myofiber Architecture on Regional Stress and Strain in the Canine Left Ventricle , 2000 .

[29]  R L Winslow,et al.  Generation and propagation of ectopic beats induced by spatially localized Na–K pump inhibition in atrial network models , 1993, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[30]  G Plank,et al.  Computational tools for modeling electrical activity in cardiac tissue. , 2003, Journal of electrocardiology.

[31]  A Garfinkel,et al.  Scroll wave dynamics in a three-dimensional cardiac tissue model: roles of restitution, thickness, and fiber rotation. , 2000, Biophysical journal.

[32]  Valentina Tozzini,et al.  Coarse-grained models for proteins. , 2005, Current opinion in structural biology.

[33]  Gunnar Seemann,et al.  Quantitative Reconstruction of Cardiac Electromechanics in Human Myocardium: , 2003, Journal of cardiovascular electrophysiology.

[34]  R. H. Mitchell,et al.  Cellular automaton model of ventricular fibrillation , 1992, IEEE Transactions on Biomedical Engineering.

[35]  I. Shrivastava,et al.  Simulations of ion permeation through a potassium channel: molecular dynamics of KcsA in a phospholipid bilayer. , 2000, Biophysical journal.

[36]  Shin-Ho Chung,et al.  Tests of continuum theories as models of ion channels. II. Poisson-Nernst-Planck theory versus brownian dynamics. , 2000, Biophysical journal.

[37]  Klaus Schulten,et al.  Molecular dynamics simulations of proteins in lipid bilayers. , 2005, Current opinion in structural biology.

[38]  A. Atilgan,et al.  Direct evaluation of thermal fluctuations in proteins using a single-parameter harmonic potential. , 1997, Folding & design.

[39]  P. Hunter,et al.  Integration from proteins to organs: the Physiome Project , 2003, Nature Reviews Molecular Cell Biology.

[40]  Alfio Quarteroni,et al.  Complex Systems in Biomedicine , 2006 .

[41]  E. Werner Meeting Report: The Future and Limits of Systems Biology , 2005, Science's STKE.

[42]  Thomas F. Miller,et al.  Quantum free energies of the conformers of glycine on an ab initio potential energy surfaceElectronic supplementary information (ESI) available: Calculated harmonic frequencies of the glycine conformers. See http://www.rsc.org/suppdata/cp/b3/b314644h/ , 2004 .

[43]  B. Roux,et al.  Molecular dynamics of the KcsA K(+) channel in a bilayer membrane. , 2000, Biophysical journal.

[44]  J. Ruppersberg Ion Channels in Excitable Membranes , 1996 .

[45]  J. Clark,et al.  Mathematical model of an adult human atrial cell: the role of K+ currents in repolarization. , 1998, Circulation research.

[46]  Richard H Clayton,et al.  Computational framework for simulating the mechanisms and ECG of re-entrant ventricular fibrillation , 2002, Physiological measurement.

[47]  H. Murakoshi,et al.  ION CONDUCTION AND SELECTIVITY IN K + CHANNELS , 2006 .

[48]  M R Franz,et al.  Effect of Sustained Load on Dispersion of Ventricular Repolarization and Conduction Time in the Isolated Intact Rabbit Heart , 1996, Journal of cardiovascular electrophysiology.

[49]  L. Atlas,et al.  Multiscale Modeling of Cardiac Cellular Energetics , 2005, Annals of the New York Academy of Sciences.

[50]  Susan Dumps,et al.  A model study. , 1988, Nursing standard (Royal College of Nursing (Great Britain) : 1987).

[51]  Kazuo Nakazawa,et al.  Spiral Wave Control by a Localized Stimulus: , 2004, Journal of cardiovascular electrophysiology.

[52]  Roy C. P. Kerckhoffs,et al.  Computational Methods for Cardiac Electromechanics , 2006, Proceedings of the IEEE.

[53]  Natalia Trayanova,et al.  Cardiac vulnerability to electric shocks during phase 1A of acute global ischemia. , 2004, Heart rhythm.

[54]  Toshiaki Hisada,et al.  Analysis of Fluid-Structure Interaction Problems with Structural Buckling and Large Domain Changes by ALE Finite Element Method. , 2001 .

[55]  Ling Xia,et al.  Analysis of cardiac ventricular wall motion based on a three-dimensional electromechanical biventricular model , 2005, Physics in medicine and biology.

[56]  A V Panfilov,et al.  Generation of Reentry in Anisotropic Myocardium , 1993, Journal of cardiovascular electrophysiology.

[57]  P. Hunter,et al.  Mathematical model of geometry and fibrous structure of the heart. , 1991, The American journal of physiology.

[58]  A. McCulloch,et al.  Computational model of three-dimensional cardiac electromechanics , 2002 .

[59]  Andrew D. McCulloch,et al.  Computational Methods for Cardiac Electrophysiology , 2004 .

[60]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

[61]  H. Kitano Systems Biology: A Brief Overview , 2002, Science.

[62]  Philip C Biggin,et al.  Ion channel gating: insights via molecular simulations , 2003, FEBS letters.

[63]  A. McCulloch,et al.  Relationship Between Regional Shortening and Asynchronous Electrical Activation in a Three‐Dimensional Model of Ventricular Electromechanics , 2003, Journal of cardiovascular electrophysiology.

[64]  James P. Keener,et al.  Mathematical physiology , 1998 .

[65]  Natalia A. Trayanova,et al.  Computational techniques for solving the bidomain equations in three dimensions , 2002, IEEE Transactions on Biomedical Engineering.

[66]  Denis Noble,et al.  Cellular Open Resource (COR): a Public CellML Based Environment for Modeling Biological Function , 2003, Int. J. Bifurc. Chaos.

[67]  J L Puglisi,et al.  LabHEART: an interactive computer model of rabbit ventricular myocyte ion channels and Ca transport. , 2001, American journal of physiology. Cell physiology.

[68]  D. Noble Modeling the Heart--from Genes to Cells to the Whole Organ , 2002, Science.

[69]  H M Byrne,et al.  The influence of growth-induced stress from the surrounding medium on the development of multicell spheroids , 2001, Journal of mathematical biology.

[70]  D. Noble,et al.  Reconstruction of the electrical activity of cardiac Purkinje fibres. , 1975, The Journal of physiology.

[71]  Toshiaki Hisada,et al.  Multiphysics simulation of left ventricular filling dynamics using fluid-structure interaction finite element method. , 2004, Biophysical journal.

[72]  N. El-Sherif,et al.  A logical state model of reentrant ventricular activation , 1990, IEEE Transactions on Biomedical Engineering.

[73]  Jeppe Olsen,et al.  A coupled cluster and full configuration interaction study of CN and CN , 2004 .

[74]  A. Hodgkin,et al.  A quantitative description of membrane current and its application to conduction and excitation in nerve , 1990 .

[75]  Tamar Schlick,et al.  Molecular Modeling and Simulation: An Interdisciplinary Guide , 2010 .

[76]  E. McVeigh,et al.  Electromechanics of paced left ventricle simulated by straightforward mathematical model: comparison with experiments. , 2005, American journal of physiology. Heart and circulatory physiology.

[77]  A V Panfilov,et al.  Phase singularities and filaments: simplifying complexity in computational models of ventricular fibrillation. , 2006, Progress in biophysics and molecular biology.

[78]  J. Humphrey Continuum biomechanics of soft biological tissues , 2003 .

[79]  J. Bonfield,et al.  Finishing the euchromatic sequence of the human genome , 2004, Nature.

[80]  Arthur C. Guyton,et al.  Handbook of Physiology—The Cardiovascular System , 1985 .

[81]  P. Hunter,et al.  Stretch-induced changes in heart rate and rhythm: clinical observations, experiments and mathematical models. , 1999, Progress in biophysics and molecular biology.

[82]  D. Levitt,et al.  Modeling of Ion Channels , 1999, The Journal of general physiology.

[83]  G S Kassab,et al.  Morphometry of pig coronary arterial trees. , 1993, The American journal of physiology.

[84]  Tianhai Tian,et al.  A multi-scaled approach for simulating chemical reaction systems. , 2004, Progress in biophysics and molecular biology.

[85]  P. Hunter,et al.  New developments in a strongly coupled cardiac electromechanical model. , 2005, Europace : European pacing, arrhythmias, and cardiac electrophysiology : journal of the working groups on cardiac pacing, arrhythmias, and cardiac cellular electrophysiology of the European Society of Cardiology.

[86]  Yoram Rudy,et al.  Computational biology in the study of cardiac ion channels and cell electrophysiology , 2006, Quarterly Reviews of Biophysics.

[87]  Peter Kohl,et al.  Myocardial ischemia lowers precordial thump efficacy: an inquiry into mechanisms using three-dimensional simulations. , 2006, Heart rhythm.

[88]  Leo K. Cheng,et al.  Cardiac electrical activity--from heart to body surface and back again. , 2003, Journal of electrocardiology.

[89]  G S Kassab,et al.  Topology and dimensions of pig coronary capillary network. , 1994, The American journal of physiology.

[90]  Qiang Shi,et al.  Mixed atomistic and coarse-grained molecular dynamics: simulation of a membrane-bound ion channel. , 2006, The journal of physical chemistry. B.

[91]  Colleen E. Clancy,et al.  Na+ Channel Mutation That Causes Both Brugada and Long-QT Syndrome Phenotypes: A Simulation Study of Mechanism , 2002, Circulation.

[92]  A. McCulloch,et al.  Three-dimensional analysis of regional cardiac function: a model of rabbit ventricular anatomy. , 1998, Progress in biophysics and molecular biology.

[93]  A V Holden,et al.  Computer simulation of re-entry sources in myocardium in two and three dimensions. , 1993, Journal of theoretical biology.

[94]  Denis Noble,et al.  Computational modelling of biological systems: tools and visions , 2000, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[95]  Satoshi Matsuoka,et al.  Role of individual ionic current systems in ventricular cells hypothesized by a model study. , 2003, The Japanese journal of physiology.

[96]  Serdar Kuyucak,et al.  Models of permeation in ion channels , 2001 .

[97]  Peter Kohl,et al.  Induction of ventricular arrhythmias following mechanical impact: A simulation study in 3D , 2004, Journal of Molecular Histology.

[98]  F. Fenton,et al.  Vortex dynamics in three-dimensional continuous myocardium with fiber rotation: Filament instability and fibrillation. , 1998, Chaos.

[99]  Andrew D McCulloch,et al.  An ionic model of stretch-activated and stretch-modulated currents in rabbit ventricular myocytes. , 2005, Europace : European pacing, arrhythmias, and cardiac electrophysiology : journal of the working groups on cardiac pacing, arrhythmias, and cardiac cellular electrophysiology of the European Society of Cardiology.

[100]  Peter J Hunter,et al.  Sarcomere length changes in a 3D mathematical model of the pig ventricles. , 2003, Progress in biophysics and molecular biology.

[101]  Antonio Fasano,et al.  Mathematical modelling of tumour growth and treatment , 2006 .

[102]  Richard H. Clayton,et al.  Filament behavior in a computational model of ventricular fibrillation in the canine heart , 2004, IEEE Transactions on Biomedical Engineering.

[103]  D. P. Nickerson,et al.  A model of cardiac cellular electromechanics , 2001, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[104]  P. Hunter,et al.  Modelling the mechanical properties of cardiac muscle. , 1998, Progress in biophysics and molecular biology.

[105]  Eric J. Sorin,et al.  How well can simulation predict protein folding kinetics and thermodynamics? , 2005, Annual review of biophysics and biomolecular structure.

[106]  E. Jakobsson,et al.  Hierarchical approach to predicting permeation in ion channels. , 2001, Biophysical journal.

[107]  David Harel,et al.  A Turing-like test for biological modeling , 2005, Nature Biotechnology.

[108]  N. G. Sepulveda,et al.  Current injection into a two-dimensional anisotropic bidomain. , 1989, Biophysical journal.

[109]  F. Fenton,et al.  Modeling wave propagation in realistic heart geometries using the phase-field method. , 2005, Chaos.

[110]  M. Lai,et al.  An Immersed Boundary Method with Formal Second-Order Accuracy and Reduced Numerical Viscosity , 2000 .

[111]  Ron Elber,et al.  Long-timescale simulation methods. , 2005, Current opinion in structural biology.

[112]  S. Yoshizawa,et al.  An Active Pulse Transmission Line Simulating Nerve Axon , 1962, Proceedings of the IRE.

[113]  Natalia A Trayanova,et al.  Differences Between Left and Right Ventricular Chamber Geometry Affect Cardiac Vulnerability to Electric Shocks , 2005, Circulation research.

[114]  E. Sobie,et al.  Nonuniform responses of transmembrane potential during electric field stimulation of single cardiac cells. , 1999, American journal of physiology. Heart and circulatory physiology.

[115]  D. Noble,et al.  Rectifying Properties of Heart Muscle , 1960, Nature.

[116]  P. Hunter,et al.  Computational physiology and the physiome project , 2004, Experimental physiology.

[117]  Rodrigo Weber dos Santos,et al.  Algebraic Multigrid Preconditioner for the Cardiac Bidomain Model , 2007, IEEE Transactions on Biomedical Engineering.

[118]  S. Cowin,et al.  Biomechanics: Mechanical Properties of Living Tissues, 2nd ed. , 1994 .

[119]  Benoît Roux,et al.  Ion conduction and selectivity in K(+) channels. , 2005, Annual review of biophysics and biomolecular structure.

[120]  C. Luo,et al.  A dynamic model of the cardiac ventricular action potential. II. Afterdepolarizations, triggered activity, and potentiation. , 1994, Circulation research.

[121]  S. Hanash,et al.  Challenges in deriving high-confidence protein identifications from data gathered by a HUPO plasma proteome collaborative study , 2006, Nature Biotechnology.

[122]  N. Trayanova,et al.  Cardiac Defibrillation and the Role of Mechanoelectric Feedback in Postshock Arrhythmogenesis , 2006, Annals of the New York Academy of Sciences.

[123]  Y Rudy,et al.  Action potential and contractility changes in [Na(+)](i) overloaded cardiac myocytes: a simulation study. , 2000, Biophysical journal.

[124]  P. J. Hunter,et al.  Generation of an Anatomically Based Three-Dimensional Model of the Conducting Airways , 2000, Annals of Biomedical Engineering.

[125]  Robert F Gilmour,et al.  Ionic mechanism of electrical alternans. , 2002, American journal of physiology. Heart and circulatory physiology.

[126]  W. Krassowska,et al.  Electroporation and Shock-Induced Transmembrane Potential in a Cardiac Fiber During Defibrillation Strength Shocks , 1998, Annals of Biomedical Engineering.

[127]  W. Rheinboldt,et al.  A COMPUTER MODEL OF ATRIAL FIBRILLATION. , 1964, American heart journal.

[128]  P. Hunter,et al.  Laminar structure of the heart: ventricular myocyte arrangement and connective tissue architecture in the dog. , 1995, The American journal of physiology.

[129]  Shawn T. Brown,et al.  Advances in methods and algorithms in a modern quantum chemistry program package. , 2006, Physical chemistry chemical physics : PCCP.

[130]  Serdar Kuyucak,et al.  Recent advances in ion channel research. , 2002, Biochimica et biophysica acta.

[131]  B. Chait,et al.  The structure of the potassium channel: molecular basis of K+ conduction and selectivity. , 1998, Science.

[132]  P. Schleyer Encyclopedia of computational chemistry , 1998 .

[133]  Luigi Preziosi,et al.  Cancer Modelling and Simulation , 2003 .

[134]  Jay D. Humphrey,et al.  Review Paper: Continuum biomechanics of soft biological tissues , 2003, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[135]  R. FitzHugh Impulses and Physiological States in Theoretical Models of Nerve Membrane. , 1961, Biophysical journal.

[136]  J. Whiteley Some factors affecting pulmonary oxygen transport. , 2006, Mathematical biosciences.

[137]  James P. Keener,et al.  An Immersed Interface Method for Solving Anisotropic Elliptic Boundary Value Problems in Three Dimensions , 2003, SIAM J. Sci. Comput..

[138]  C. Luo,et al.  A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes. , 1994, Circulation research.

[139]  A. Hodgkin,et al.  A quantitative description of membrane current and its application to conduction and excitation in nerve , 1952, The Journal of physiology.

[140]  C. Peskin,et al.  A three-dimensional computer model of the human heart for studying cardiac fluid dynamics , 2000, SIGGRAPH 2000.

[141]  Shin-Ho Chung,et al.  Permeation of ions across the potassium channel: Brownian dynamics studies. , 1999, Biophysical journal.

[142]  Natalia Trayanova,et al.  Defibrillation of the heart: insights into mechanisms from modelling studies , 2006, Experimental physiology.

[143]  Gregory A Voth,et al.  A multiscale coarse-graining method for biomolecular systems. , 2005, The journal of physical chemistry. B.

[144]  Andrew J. Pullan,et al.  A Finite Element Method for an Eikonal Equation Model of Myocardial Excitation Wavefront Propagation , 2002, SIAM J. Appl. Math..

[145]  Peter J. Hunter,et al.  Modeling RBC and Neutrophil Distribution Through an Anatomically Based Pulmonary Capillary Network , 2004, Annals of Biomedical Engineering.

[146]  Alexander D. MacKerell,et al.  CHARMM: The Energy Function and Its Parameterization , 2002 .

[147]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[148]  Denis Noble,et al.  The IUPS human physiome project , 2002, Pflügers Archiv.

[149]  L. Clerc Directional differences of impulse spread in trabecular muscle from mammalian heart. , 1976, The Journal of physiology.

[150]  J. Bonfield,et al.  Finishing the euchromatic sequence of the human genome , 2004, Nature.

[151]  Roy C. P. Kerckhoffs,et al.  Homogeneity of Cardiac Contraction Despite Physiological Asynchrony of Depolarization: A Model Study , 2003, Annals of Biomedical Engineering.

[152]  J A Negroni,et al.  A cardiac muscle model relating sarcomere dynamics to calcium kinetics. , 1996, Journal of molecular and cellular cardiology.

[153]  P. Maini,et al.  Multiscale modeling in biology , 2007 .

[154]  Tirion,et al.  Large Amplitude Elastic Motions in Proteins from a Single-Parameter, Atomic Analysis. , 1996, Physical review letters.

[155]  R. Plonsey Bioelectric sources arising in excitable fibers (Alza lecture) , 2006, Annals of Biomedical Engineering.

[156]  G. R. Smith,et al.  Simulation approaches to ion channel structure–function relationships , 2001, Quarterly Reviews of Biophysics.

[157]  S. Chung,et al.  Tests of continuum theories as models of ion channels. I. Poisson-Boltzmann theory versus Brownian dynamics. , 2000, Biophysical journal.