title: Fitting local repolarization parameters in cardiac reaction-diffusion models in the presence of electrotonic coupling.
abstract: Repolarization gradients contribute to arrhythmogenicity. In reaction-diffusion models of cardiac tissue, heterogeneities in action potential duration (APD) can be created by locally modifying an intrinsic membrane kinetics parameter. Electrotonic coupling, however, acts as a confounding factor that modulates APD dispersion.
  
  
  
  We developed an algorithm based on a quasi-Newton method that iteratively adjusts the spatial distribution of a membrane parameter to reproduce a pre-defined target APD map in a coupled tissue. The method assumes that the relation between the adjustable parameter and APD is bijective in an isolated cell. Each iteration of the algorithm involved simulating the cardiac reaction-diffusion system with the updated parameter profile for one beat and extracting the APD map. The algorithm was extended to simultaneous estimation of two parameter profiles based on two APD maps at different repolarization thresholds.
  
  
  
  The method was validated in 1D, 2D and 3D atrial tissues using synthetic target APD maps with controllable total variation and maximum APD gradient. The adjustable parameter was local acetylcholine concentration. The iterations converged provided that APD gradients were not too steep. Convergence was found to be faster 2-5 iterations) when the maximal gradient was less steep, when APD range was smaller and when tissue conductivity was reduced.
  
  
  
  This algorithm provides a tool to automatically generate arrhythmogenic substrates with controllable repolarization gradients and possibly incorporate experimental APD maps into computer models.
authors.affiliation: Université de Montréal, Département de Pharmacologie et Physiologie , Institut de Génie Biomédical, Montréal, Canada; Hôpital du Sacré-Coeur de Montréal, Centre de Recherche, Montréal, Canada.
authors.forename: Akshay Kota
authors.initials: AK
authors.lastname: Aswath Kumar
authors.affiliation: Université de Montréal, Département de Pharmacologie et Physiologie , Institut de Génie Biomédical, Montréal, Canada; Hôpital du Sacré-Coeur de Montréal, Centre de Recherche, Montréal, Canada.
authors.forename: Angelina
authors.initials: A
authors.lastname: Drahi
authors.affiliation: Université de Montréal, Département de Pharmacologie et Physiologie , Institut de Génie Biomédical, Montréal, Canada; Hôpital du Sacré-Coeur de Montréal, Centre de Recherche, Montréal, Canada. Electronic address: vincent.jacquemet@umontreal.ca.
authors.forename: Vincent
authors.initials: V
authors.lastname: Jacquemet
mesh_terms: D000200:Action Potentials / Q000502:physiology*; D000465:Algorithms; D000818:Animals; D018780:Body Surface Potential Mapping / Q000379:methods*; D003198:Computer Simulation; D006329:Heart Conduction System / Q000502:physiology*; D006801:Humans; D008564:Membrane Potentials / Q000502:physiology*; D008955:Models, Cardiovascular; D032383:Myocytes, Cardiac / Q000502:physiology*; D009716:Numerical Analysis, Computer-Assisted
keywords: Action potential duration; Cardiac electrophysiology; Cell coupling; Computer modeling; Parameter estimation