Rapid Genetic Algorithm Optimization of a Mouse Computational Model: Benefits for Anthropomorphization of Neonatal Mouse Cardiomyocytes

While the mouse presents an invaluable experimental model organism in biology, its usefulness in cardiac arrhythmia research is limited in some aspects due to major electrophysiological differences between murine and human action potentials (APs). As previously described, these species-specific traits can be partly overcome by application of a cell-type transforming clamp (CTC) to anthropomorphize the murine cardiac AP. CTC is a hybrid experimental-computational dynamic clamp technique, in which a computationally calculated time-dependent current is inserted into a cell in real-time, to compensate for the differences between sarcolemmal currents of that cell (e.g., murine) and the desired species (e.g., human). For effective CTC performance, mismatch between the measured cell and a mathematical model used to mimic the measured AP must be minimal. We have developed a genetic algorithm (GA) approach that rapidly tunes a mathematical model to reproduce the AP of the murine cardiac myocyte under study. Compared to a prior implementation that used a template-based model selection approach, we show that GA optimization to a cell-specific model results in a much better recapitulation of the desired AP morphology with CTC. This improvement was more pronounced when anthropomorphizing neonatal mouse cardiomyocytes to human-like APs than to guinea pig APs. CTC may be useful for a wide range of applications, from screening effects of pharmaceutical compounds on ion channel activity, to exploring variations in the mouse or human genome. Rapid GA optimization of a cell-specific mathematical model improves CTC performance and may therefore expand the applicability and usage of the CTC technique.

[1]  Eve Marder,et al.  Structure and visualization of high-dimensional conductance spaces. , 2006, Journal of neurophysiology.

[2]  Alan Garfinkel,et al.  Arrhythmogenic consequences of myofibroblast-myocyte coupling. , 2012, Cardiovascular research.

[3]  Erik De Schutter,et al.  Complex Parameter Landscape for a Complex Neuron Model , 2006, PLoS Comput. Biol..

[4]  Alan Garfinkel,et al.  Shaping a new Ca2+ conductance to suppress early afterdepolarizations in cardiac myocytes , 2011, The Journal of physiology.

[5]  Ronald Wilders,et al.  UvA-DARE ( Digital Academic Repository ) HERG channel ( dys ) function revealed by dynamic action potential clamp technique , 2004 .

[6]  Thomas O'Hara,et al.  Quantitative comparison of cardiac ventricular myocyte electrophysiology and response to drugs in human and nonhuman species. , 2012, American journal of physiology. Heart and circulatory physiology.

[7]  A. V. van Ginneken,et al.  Long‐QT syndrome‐related sodium channel mutations probed by the dynamic action potential clamp technique , 2006, The Journal of physiology.

[8]  R. Winslow,et al.  A computational model of the human left-ventricular epicardial myocyte. , 2004, Biophysical journal.

[9]  D. Christini,et al.  Instability in action potential morphology underlies phase 2 reentry: a mathematical modeling study. , 2009, Heart rhythm.

[10]  K H W J Ten Tusscher,et al.  Cell model for efficient simulation of wave propagation in human ventricular tissue under normal and pathological conditions , 2006, Physics in medicine and biology.

[11]  T. Colatsky,et al.  Inhibition of cardiac delayed rectifier K+ current by overexpression of the long-QT syndrome HERG G628S mutation in transgenic mice. , 1998, Circulation research.

[12]  John A. White,et al.  Effects of imperfect dynamic clamp: Computational and experimental results , 2008, Journal of Neuroscience Methods.

[13]  N. Rosenthal,et al.  Analysis of cardiac myocyte biology in transgenic mice: a protocol for preparation of neonatal mouse cardiac myocyte cultures. , 2010, Methods in molecular biology.

[14]  J. Mertz,et al.  Control of Local Intracellular Calcium Concentration with Dynamic-Clamp Controlled 2-Photon Uncaging , 2011, PloS one.

[15]  David J. Christini,et al.  Real-time Experiment Interface for biological control applications , 2010, 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology.

[16]  S. Verheule,et al.  Cardiac electrophysiology in mice: a matter of size , 2012, Front. Physio..

[17]  S. Grandy,et al.  Postnatal development has a marked effect on ventricular repolarization in mice. , 2007, American journal of physiology. Heart and circulatory physiology.

[18]  John A. White,et al.  GenNet: A Platform for Hybrid Network Experiments , 2011, Front. Neuroinform..

[19]  D. Christini,et al.  Anthropomorphizing the mouse cardiac action potential via a novel dynamic clamp method. , 2009, Biophysical journal.

[20]  C. Henriquez,et al.  An integrative model of mouse cardiac electrophysiology from cell to torso. , 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.

[21]  Gerald Friedland,et al.  Revisiting a basic function on current CPUs : A fast logarithm implementation with adjustable accuracy , 2007 .

[22]  E. Marbán,et al.  Calcium cycling and contractile activation in intact mouse cardiac muscle , 1998, The Journal of physiology.

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

[24]  Nicol N. Schraudolph,et al.  A Fast, Compact Approximation of the Exponential Function , 1999, Neural Computation.

[25]  Francis A. Ortega,et al.  Cardiac myocyte model parameter sensitivity analysis and model transformation using a genetic algorithm , 2011, GECCO.

[26]  Henry Markram,et al.  A Novel Multiple Objective Optimization Framework for Constraining Conductance-Based Neuron Models by Experimental Data , 2007, Front. Neurosci..

[27]  James M. Bower,et al.  A Comparative Survey of Automated Parameter-Search Methods for Compartmental Neural Models , 1999, Journal of Computational Neuroscience.

[28]  Victor Yodaiken,et al.  A Real-Time Linux , 2000 .

[29]  Dieter Jaeger,et al.  The use of automated parameter searches to improve ion channel kinetics for neural modeling , 2011, Journal of Computational Neuroscience.

[30]  E. Neher Correction for liquid junction potentials in patch clamp experiments. , 1992, Methods in enzymology.

[31]  Robert J Butera,et al.  MRCI: a flexible real-time dynamic clamp system for electrophysiology experiments , 2004, Journal of Neuroscience Methods.

[32]  D. Bers Calcium cycling and signaling in cardiac myocytes. , 2008, Annual review of physiology.

[33]  W. Giles,et al.  A mathematical model of action potential heterogeneity in adult rat left ventricular myocytes. , 2001, Biophysical journal.

[34]  Christopher R. Johnson,et al.  Three-dimensional Propagation in Mathematic Models: Integrative Model of the Mouse Heart , 2004 .

[35]  Gavin C. Cawley,et al.  On a Fast, Compact Approximation of the Exponential Function , 2000, Neural Computation.

[36]  Ronald Wilders,et al.  Dynamic clamp: a powerful tool in cardiac electrophysiology , 2006, The Journal of physiology.

[37]  J. Nerbonne Studying cardiac arrhythmias in the mouse--a reasonable model for probing mechanisms? , 2004, Trends in cardiovascular medicine.

[38]  Deanna M. Church,et al.  Genome Reference Consortium , 2013 .

[39]  Edward J. Vigmond,et al.  Atrial cell action potential parameter fitting using genetic algorithms , 2005, Medical and Biological Engineering and Computing.

[40]  Eve Marder,et al.  The dynamic clamp comes of age , 2004, Trends in Neurosciences.

[41]  David J. Christini,et al.  Real-Time Linux Dynamic Clamp: A Fast and Flexible Way to Construct Virtual Ion Channels in Living Cells , 2001, Annals of Biomedical Engineering.

[42]  N. Trayanova,et al.  A Computational Model to Predict the Effects of Class I Anti-Arrhythmic Drugs on Ventricular Rhythms , 2011, Science Translational Medicine.