Systems approach to understanding electromechanical activity in the human heart: a national heart, lung, and blood institute workshop summary

The physics of heart rhythm disorders

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.

Diffuse fibrosis and repolarization disorders explain ventricular arrhythmias in Brugada syndrome: a computational study

In this work, we reported a computational study to quantitatively determine the individual contributions of three candidate arrhythmic factors associated with Brugada Syndrome. In particular, we focused our analysis on the role of structural abnormalities, dispersion of repolarization, and size of the diseased region. We developed a human phenomenological model capable of replicating the action potential characteristics both in Brugada Syndrome and in healthy conditions. Inspired by physiological observations, we employed the phenomenological model in a 2D geometry resembling the pathological RVOT coupled with healthy epicardial tissue. We assessed the insurgence of sustained reentry as a function of electrophysiological and structural abnormalities. Our computational study indicates that both structural and repolarization abnormalities are essential to induce sustained reentry. Furthermore, our results suggest that neither dispersion of repolarization nor structural abnormalities are sufficient on their own to induce sustained reentry. It should be noted how our study seems to explain an arrhythmic mechanism that unifies the classic repolarization and depolarization hypotheses of the pathophysiology of the Brugada Syndrome. Finally, we believe that this work may offer a new perspective on the computational and clinical investigation of Brugada Syndrome and its arrhythmic behaviour.

Characteristics of Bisphenol Cardiotoxicity: Impaired Excitability, Contractility, and Relaxation

Bisphenol a (BPA) is a high production volume chemical that is frequently used to manufacture epoxy resins and polycarbonate plastics. BPA-containing products are now pervasive, and as a result, biomonitoring studies report widespread exposure in > 90% of adults, adolescents, and children. Both epidemiological and experimental studies have reported associations between BPA exposure and adverse cardiovascular health outcomes. With increasing concerns regarding BPA exposure, a few structurally similar bisphenol chemicals have been introduced as replacements, including bisphenol s (BPS) and bisphenol f (BPF). In accordance with the recently established "Key characteristics of cardiovascular toxicants", we reviewed the literature to highlight the immediate effects of bisphenol chemicals on (1) cardiac excitability, and (2) contractility and relaxation. BPA inhibits key cardiac ion channels, impairs cardiac excitability, and acts as a more potent inhibitor as compared to BPF and BPS. Through the inhibition of calcium current, some studies report that bisphenol chemicals can act as negative inotropic agents. Yet, others suggest that low dose exposures may increase contractility and precipitate triggered arrhythmias via the phosphorylation of key calcium handling proteins. Accordingly, we propose additional considerations for future work to comprehensively address the cardiac safety profile of BPA, as compared to replacement chemicals.

Quantitative cross-species translators of cardiac myocyte electrophysiology: Model training, experimental validation, and applications

Animal experimentation is key in the evaluation of cardiac efficacy and safety of novel therapeutic compounds. However, interspecies differences in the mechanisms regulating excitation-contraction coupling can limit the translation of experimental findings from animal models to human physiology and undermine the assessment of drugs’ efficacy and safety. Here, we built a suite of translators for quantitatively mapping electrophysiological responses in ventricular myocytes across species. We trained these statistical operators using a broad dataset obtained by simulating populations of our biophysically detailed computational models of action potential and Ca2+ transient in mouse, rabbit, and human. We then tested our translators against experimental data describing the response to stimuli, such as ion channel block, change in beating rate, and β-adrenergic challenge. We demonstrate that this approach is well suited to predicting the effects of perturbations across different species or experimental conditions and suggest its integration into mechanistic studies and drug development pipelines.

A computationally efficient dynamic model of human epicardial tissue

We present a new phenomenological model of human ventricular epicardial cells and we test its reentry dynamics. The model is derived from the Rogers-McCulloch formulation of the FitzHugh-Nagumo equations and represents the total ionic current divided into three contributions corresponding to the excitatory, recovery and transient outward currents. Our model reproduces the main characteristics of human epicardial tissue, including action potential amplitude and morphology, upstroke velocity, and action potential duration and conduction velocity restitution curves. The reentry dynamics is stable, and the dominant period is about 270 ms, which is comparable to clinical values. The proposed model is the first phenomenological model able to accurately resemble human experimental data by using only 3 state variables and 17 parameters. Indeed, it is more computationally efficient than existing models (i.e., almost two times faster than the minimal ventricular model). Beyond the computational efficiency, the low number of parameters facilitates the process of fitting the model to the experimental data.

Transgenic Rabbit Models in Proarrhythmia Research

Drug-induced proarrhythmia constitutes a potentially lethal side effect of various drugs. Most often, this proarrhythmia is mechanistically linked to the drug's potential to interact with repolarizing cardiac ion channels causing a prolongation of the QT interval in the ECG. Despite sophisticated screening approaches during drug development, reliable prediction of proarrhythmia remains very challenging. Although drug-induced long-QT-related proarrhythmia is often favored by conditions or diseases that impair the individual's repolarization reserve, most cellular, tissue, and whole animal model systems used for drug safety screening are based on normal, healthy models. In recent years, several transgenic rabbit models for different types of long QT syndromes (LQTS) with differences in the extent of impairment in repolarization reserve have been generated. These might be useful for screening/prediction of a drug's potential for long-QT-related proarrhythmia, particularly as different repolarizing cardiac ion channels are impaired in the different models. In this review, we summarize the electrophysiological characteristics of the available transgenic LQTS rabbit models, and the pharmacological proof-of-principle studies that have been performed with these models—highlighting the advantages and disadvantages of LQTS models for proarrhythmia research. In the end, we give an outlook on potential future directions and novel models.

Closing the 'phenotype gap' in precision medicine: improving what we measure to understand complex disease mechanisms

The central concept underlying precision medicine is a mechanistic understanding of each disease and its response to therapy sufficient to direct a specific intervention. To execute on this vision requires parsing incompletely defined disease syndromes into discrete mechanistic subsets and developing interventions to precisely address each of these etiologically distinct entities. This will require substantial adjustment of traditional paradigms which have tended to aggregate high-level phenotypes with very different etiologies. In the current environment, where diagnoses are not mechanistic, drug development has become so expensive that it is now impractical to imagine the cost-effective creation of new interventions for many prevalent chronic conditions. The vision of precision medicine also argues for a much more seamless integration of research and development with clinical care, where shared taxonomies will enable every clinical interaction to inform our collective understanding of disease mechanisms and drug responses. Ideally, this would be executed in ways that drive real-time and real-world discovery, innovation, translation, and implementation. Only in oncology, where at least some of the biology is accessible through surgical excision of the diseased tissue or liquid biopsy, has "co-clinical" modeling proven feasible. In most common germline disorders, while genetics often reveal the causal mutations, there still remain substantial barriers to efficient disease modeling. Aggregation of similar disorders under single diagnostic labels has directly contributed to the paucity of etiologic and mechanistic understanding by directly reducing the resolution of any subsequent studies. Existing clinical phenotypes are typically anatomic, physiologic, or histologic, and result in a substantial mismatch in information content between the phenomes in humans or in animal 'models' and the variation in the genome. This lack of one-to-one mapping of discrete mechanisms between disease and animal models causes a failure of translation and is one form of 'phenotype gap.' In this review, we will focus on the origins of the phenotyping deficit and approaches that may be considered to bridge the gap, creating shared taxonomies between human diseases and relevant models, using cardiovascular examples.

Successful knock-in of Hypertrophic Cardiomyopathy-mutation R723G into the MYH7 gene mimics HCM pathology in pigs

Familial Hypertrophic Cardiomyopathy (HCM) is the most common inherited cardiac disease. About 30% of the patients are heterozygous for mutations in the MYH7 gene encoding the ß-myosin heavy chain (MyHC). Hallmarks of HCM are cardiomyocyte disarray and hypertrophy of the left ventricle, the symptoms range from slight arrhythmias to sudden cardiac death or heart failure. To gain insight into the underlying mechanisms of the diseases' etiology we aimed to generate genome edited pigs with an HCM-mutation. We used TALEN-mediated genome editing and successfully introduced the HCM-point mutation R723G into the MYH7 gene of porcine fibroblasts and subsequently cloned pigs that were heterozygous for the HCM-mutation R723G. No off-target effects were determined in the R723G-pigs. Surprisingly, the animals died within 24 h post partem, probably due to heart failure as indicated by a shift in the a/ß-MyHC ratio in the left ventricle. Most interestingly, the neonatal pigs displayed features of HCM, including mild myocyte disarray, malformed nuclei, and MYH7-overexpression. The finding of HCM-specific pathology in neonatal R723G-piglets suggests a very early onset of the disease and highlights the importance of novel large animal models for studying causative mechanisms and long-term progression of human cardiac diseases.

Larger rate dependence of late sodium current in cardiac Purkinje cells: A potential link to arrhythmogenesis

BACKGROUND

Purkinje cells (PCs) have a steeper rate dependence of repolarization and are more susceptible to arrhythmic activity than do ventricular myocytes (VMs). Late sodium current (I_NaL) is rate dependent and contributes to rate dependence of repolarization.

OBJECTIVE

This study sought to test our hypothesis that PCs have a larger rate dependence of I_NaL, contributing to their steeper rate dependence of repolarization and higher susceptibility to arrhythmic activity, than do VMs.

METHODS

I_NaL was recorded in isolated rabbit PCs and VMs with the whole-cell patch clamp technique. Action potential was examined using the microelectrode technique.

RESULTS

Compared with VMs, PCs exhibited a significantly larger rate dependence of I_NaL with a larger I_NaL to basic cycle length (BCL) slope. Moreover, PCs had a larger rate dependence of I_NaL decay and slower recovery kinetics. Interestingly, the larger rate dependence of I_NaL matched to a steeper rate dependence of action potential duration (APD) in PCs. The I_NaL blocker tetrodotoxin significantly blunted, while the I_NaL enhancer anemone toxin (ATX-II) significantly increased, the rate dependence of I_NaL and APD in PCs and VMs. In the presence of ATX-II, the rate dependence of I_NaL in PCs was markedly larger than that in VMs, causing a much steeper rate dependence of APD in PCs. Accordingly, PCs exhibited greater rate-dependent electrical instability and were more prone to ATX-II-induced early afterdepolarizations, which were completely inhibited by the I_NaL inhibitor ranolazine.

CONCLUSION

PCs have a significantly larger rate dependence of I_NaL than do VMs because of distinctive I_NaL decay and recovery kinetics, which contributes to their larger rate adaptation, and simultaneously predisposes them to a higher risk of arrhythmogenesis.

Myocardial Infarct Segmentation From Magnetic Resonance Images for Personalized Modeling of Cardiac Electrophysiology

Accurate representation of myocardial infarct geometry is crucial to patient-specific computational modeling of the heart in ischemic cardiomyopathy. We have developed a methodology for segmentation of left ventricular (LV) infarct from clinically acquired, two-dimensional (2D), late-gadolinium enhanced cardiac magnetic resonance (LGE-CMR) images, for personalized modeling of ventricular electrophysiology. The infarct segmentation was expressed as a continuous min-cut optimization problem, which was solved using its dual formulation, the continuous max-flow (CMF). The optimization objective comprised of a smoothness term, and a data term that quantified the similarity between image intensity histograms of segmented regions and those of a set of training images. A manual segmentation of the LV myocardium was used to initialize and constrain the developed method. The three-dimensional geometry of infarct was reconstructed from its segmentation using an implicit, shape-based interpolation method. The proposed methodology was extensively evaluated using metrics based on geometry, and outcomes of individualized electrophysiological simulations of cardiac dys(function). Several existing LV infarct segmentation approaches were implemented, and compared with the proposed method. Our results demonstrated that the CMF method was more accurate than the existing approaches in reproducing expert manual LV infarct segmentations, and in electrophysiological simulations. The infarct segmentation method we have developed and comprehensively evaluated in this study constitutes an important step in advancing clinical applications of personalized simulations of cardiac electrophysiology.

Image-based reconstruction of three-dimensional myocardial infarct geometry for patient-specific modeling of cardiac electrophysiology

PURPOSE

Accurate three-dimensional (3D) reconstruction of myocardial infarct geometry is crucial to patient-specific modeling of the heart aimed at providing therapeutic guidance in ischemic cardiomyopathy. However, myocardial infarct imaging is clinically performed using two-dimensional (2D) late-gadolinium enhanced cardiac magnetic resonance (LGE-CMR) techniques, and a method to build accurate 3D infarct reconstructions from the 2D LGE-CMR images has been lacking. The purpose of this study was to address this need.

METHODS

The authors developed a novel methodology to reconstruct 3D infarct geometry from segmented low-resolution (Lo-res) clinical LGE-CMR images. Their methodology employed the so-called logarithm of odds (LogOdds) function to implicitly represent the shape of the infarct in segmented image slices as LogOdds maps. These 2D maps were then interpolated into a 3D image, and the result transformed via the inverse of LogOdds to a binary image representing the 3D infarct geometry. To assess the efficacy of this method, the authors utilized 39 high-resolution (Hi-res) LGE-CMR images, including 36 in vivo acquisitions of human subjects with prior myocardial infarction and 3 ex vivo scans of canine hearts following coronary ligation to induce infarction. The infarct was manually segmented by trained experts in each slice of the Hi-res images, and the segmented data were downsampled to typical clinical resolution. The proposed method was then used to reconstruct 3D infarct geometry from the downsampled images, and the resulting reconstructions were compared with the manually segmented data. The method was extensively evaluated using metrics based on geometry as well as results of electrophysiological simulations of cardiac sinus rhythm and ventricular tachycardia in individual hearts. Several alternative reconstruction techniques were also implemented and compared with the proposed method.

RESULTS

The accuracy of the LogOdds method in reconstructing 3D infarct geometry, as measured by the Dice similarity coefficient, was 82.10% ± 6.58%, a significantly higher value than those of the alternative reconstruction methods. Among outcomes of electrophysiological simulations with infarct reconstructions generated by various methods, the simulation results corresponding to the LogOdds method showed the smallest deviation from those corresponding to the manual reconstructions, as measured by metrics based on both activation maps and pseudo-ECGs.

CONCLUSIONS

The authors have developed a novel method for reconstructing 3D infarct geometry from segmented slices of Lo-res clinical 2D LGE-CMR images. This method outperformed alternative approaches in reproducing expert manual 3D reconstructions and in electrophysiological simulations.

How computer simulations of the human heart can improve anti-arrhythmia therapy

Over the last decade, the state-of-the-art in cardiac computational modelling has progressed rapidly. The electrophysiological function of the heart can now be simulated with a high degree of detail and accuracy, opening the doors for simulation-guided approaches to anti-arrhythmic drug development and patient-specific therapeutic interventions. In this review, we outline the basic methodology for cardiac modelling, which has been developed and validated over decades of research. In addition, we present several recent examples of how computational models of the human heart have been used to address current clinical problems in cardiac electrophysiology. We will explore the use of simulations to improve anti-arrhythmic pacing and defibrillation interventions; to predict optimal sites for clinical ablation procedures; and to aid in the understanding and selection of arrhythmia risk markers. Together, these studies illustrate how the tremendous advances in cardiac modelling are poised to revolutionize medical treatment and prevention of arrhythmia.

TRPM4 non-selective cation channels influence action potentials in rabbit Purkinje fibres

KEY POINTS

The transient receptor potential melastatin 4 (TRPM4) inhibitor 9-phenanthrol reduces action potential duration in rabbit Purkinje fibres but not in ventricle. TRPM4-like single channel activity is observed in isolated rabbit Purkinje cells but not in ventricular cells. The TRPM4-like current develops during the notch and early repolarization phases of the action potential in Purkinje cells.

ABSTRACT

Transient receptor potential melastatin 4 (TRPM4) Ca(2+)-activated non-selective cation channel activity has been recorded in cardiomyocytes and sinus node cells from mammals. In addition, TRPM4 gene mutations are associated with human diseases of cardiac conduction, suggesting that TRPM4 plays a role in this aspect of cardiac function. Here we evaluate the TRPM4 contribution to cardiac electrophysiology of Purkinje fibres. Ventricular strips with Purkinje fibres were isolated from rabbit hearts. Intracellular microelectrodes recorded Purkinje fibre activity and the TRPM4 inhibitor 9-phenanthrol was applied to unmask potential TRPM4 contributions to the action potential. 9-Phenanthrol reduced action potential duration measured at the point of 50 and 90% repolarization with an EC50 of 32.8 and 36.1×10(-6) mol l(-1), respectively, but did not modulate ventricular action potentials. Inside-out patch-clamp recordings were used to monitor TRPM4 activity in isolated Purkinje cells. TRPM4-like single channel activity (conductance = 23.8 pS; equal permeability for Na(+) and K(+); sensitivity to voltage, Ca(2+) and 9-phenanthrol) was observed in 43% of patches from Purkinje cells but not from ventricular cells (0/16). Action potential clamp experiments performed in the whole-cell configuration revealed a transient inward 9-phenanthrol-sensitive current (peak density = -0.65 ± 0.15 pA pF(-1); n = 5) during the plateau phases of the Purkinje fibre action potential. These results show that TRPM4 influences action potential characteristics in rabbit Purkinje fibres and thus could modulate cardiac conduction and be involved in triggering arrhythmias.

Acute effects of ethanol on action potential and intracellular Ca(2+) transient in cardiac ventricular cells: a simulation study

Alcohol consumption may result in electrocardiographic changes and arrhythmias, at least partly due to effects of ethanol on cardiac ionic currents. Contractility and intracellular Ca(2+) dynamics seem to be altered as well. In this study, we integrated the available (mostly animal) experimental data into previously published models of the rat and human ventricular myocytes to assess the share of ionic current components in ethanol-induced changes in AP configuration and cytosolic Ca(2+) transient in ventricular cardiomyocytes. The rat model reproduced well the experimentally observed changes in AP duration (APD) under ethanol (slight prolongation at 0.8 mM and shortening at ≥8 mM). These changes were almost exclusively caused by the ethanol-induced alterations of I K1. The cytosolic Ca(2+) transient decreased gradually with the increasing ethanol concentration as a result of the ethanol-induced inhibition of I Ca. In the human model, ethanol produced a dose-dependent APD lengthening, dominated by ethanol effect on I Kr, the key repolarising current in human ventricles. This effect might contribute to the clinically observed proarrhythmic effects of ethanol in predisposed individuals.

The adverse cardiac effects of Di(2-ethylhexyl)phthalate and Bisphenol A

The ubiquitous nature of plastics has raised concerns pertaining to continuous exposure to plastic polymers and human health risks. Of particular concern is the use of endocrine-disrupting chemicals in plastic production, including di(2-ethylhexyl)phthalate (DEHP) and bisphenol A (BPA). Widespread and continuous exposure to DEHP and BPA occurs through dietary intake, inhalation, dermal and intravenous exposure via consumer products and medical devices. This article reviews the literature examining the relationship between DEHP and BPA exposure and cardiac toxicity. In vitro and in vivo experimental reports are outlined, as well as epidemiological studies which examine the association between these chemicals and cardiovascular outcomes. Gaps in our current knowledge are also discussed, along with future investigative endeavors that may help resolve whether DEHP and/or BPA exposure has a negative impact on cardiovascular physiology.

Ion Channels in the Heart

Optimal cardiac function depends on proper timing of excitation and contraction in various regions of the heart, as well as on appropriate heart rate. This is accomplished via specialized electrical properties of various components of the system, including the sinoatrial node, atria, atrioventricular node, His-Purkinje system, and ventricles. Here we review the major regionally determined electrical properties of these cardiac regions and present the available data regarding the molecular and ionic bases of regional cardiac function and dysfunction. Understanding these differences is of fundamental importance for the investigation of arrhythmia mechanisms and pharmacotherapy.

Computational modeling of cardiac optogenetics: Methodology overview & review of findings from simulations

Cardiac optogenetics is emerging as an exciting new potential avenue to enable spatiotemporally precise control of excitable cells and tissue in the heart with low-energy optical stimuli. This approach involves the expression of exogenous light-sensitive proteins (opsins) in target heart tissue via viral gene or cell delivery. Preliminary experiments in optogenetically-modified cells, tissue, and organisms have made great strides towards demonstrating the feasibility of basic applications, including the use of light stimuli to pace or disrupt reentrant activity. However, it remains unknown whether techniques based on this intriguing technology could be scaled up and used in humans for novel clinical applications, such as pain-free optical defibrillation or dynamic modulation of action potential shape. A key step towards answering such questions is to explore potential optogenetics-based therapies using sophisticated computer simulation tools capable of realistically representing opsin delivery and light stimulation in biophysically detailed, patient-specific models of the human heart. This review provides (1) a detailed overview of the methodological developments necessary to represent optogenetics-based solutions in existing virtual heart platforms and (2) a survey of findings that have been derived from such simulations and a critical assessment of their significance with respect to the progress of the field.

Exposure to phthalates affects calcium handling and intercellular connectivity of human stem cell-derived cardiomyocytes

Background

The pervasive nature of plastics has raised concerns about the impact of continuous exposure to plastic additives on human health. Of particular concern is the use of phthalates in the production of flexible polyvinyl chloride (PVC) products. Di-2-ethylhexyl-phthalate (DEHP) is a commonly used phthalate ester plasticizer that imparts flexibility and elasticity to PVC products. Recent epidemiological studies have reported correlations between urinary phthalate concentrations and cardiovascular disease, including an increased risk of high blood pressure and coronary risk. Yet, there is little direct evidence linking phthalate exposure to adverse effects in human cells, including cardiomyocytes.

Methods and Results

The effect of DEHP on calcium handling was examined using monolayers of gCAMP3 human embryonic stem cell-derived cardiomyocytes, which contain an endogenous calcium sensor. Cardiomyocytes were exposed to DEHP (5 – 50 μg/mL), and calcium transients were recorded using a Zeiss confocal imaging system. DEHP exposure (24 – 72 hr) had a negative chronotropic and inotropic effect on cardiomyocytes, increased the minimum threshold voltage required for external pacing, and modified connexin-43 expression. Application of Wy-14,643 (100 μM), an agonist for the peroxisome proliferator-activated receptor alpha, did not replicate DEHP’s effects on calcium transient morphology or spontaneous beating rate.

Conclusions

Phthalates can affect the normal physiology of human cardiomyocytes, including DEHP elicited perturbations in cardiac calcium handling and intercellular connectivity. Our findings call for additional studies to clarify the extent by which phthalate exposure can alter cardiac function, particularly in vulnerable patient populations who are at risk for high phthalate exposure.

Alternating membrane potential/calcium interplay underlies repetitive focal activity in a genetic model of calcium-dependent atrial arrhythmias

KEY POINTS

Atrial fibrillation is often initiated and perpetuated by abnormal electrical pulses repetitively originating from regions outside the heart's natural pacemaker. In this study we examined the causal role of abnormal calcium releases from the sarcoplasmic reticulum in producing repetitive electrical discharges in atrial cells and tissues. Calsequestrin2 is a protein that stabilizes the closed state of calcium release channels, i.e. the ryanodine receptors. In the atria from mice predisposed to abnormal calcium releases secondary to the absence of calsequestrin2, we observed abnormal repetitive electrical discharges that may lead to atrial fibrillation. Here, we report a novel pathological rhythm generator. Specifically, abnormal calcium release leads to electrical activation, which in turn results in another abnormal calcium release. This process repeats itself and thus sustains the repetitive electrical discharges. These results suggest that improving the stability of ryanodine receptors might be useful to treat atrial fibrillation.

ABSTRACT

Aberrant diastolic calcium (Ca) release due to leaky ryanodine receptors (RyR2s) has been recently associated with atrial fibrillation (AF) and catecholaminergic polymorphic ventricular tachycardia (CPVT). However, it remains unclear how diastolic Ca release contributes to the rising of rapid repetitive focal activity, which is considered as a common AF triggering mechanism. To address this question, we conducted simultaneous voltage/Ca optical mapping in atrial tissue and one-/two-dimensional confocal imaging in atrial tissue and myocytes from wild-type (WT, n = 15) and CPVT mice lacking calsequestrin 2 (Casq2(-/-), n = 45), which promotes diastolic Ca release. During β-adrenergic stimulation (100 nM isoproterenol), only Casq2(-/-) atrial myocytes showed pacing-induced self-sustained repetitive activity (31 ± 21 s vs. none in WT). Importantly, in atrial tissue, this repetitive activity could translate to Ca-dependent focal arrhythmia. Ectopic action potential (AP) firing during repetitive activity occurred only when diastolic Ca release achieved a sufficient level of synchronization. The AP, in turn, synchronized subsequent diastolic Ca release by temporally aligning multiple sources of Ca waves both within individual myocytes and throughout the atrial tissue. This alternating interplay between AP and diastolic Ca release perpetuates the self-sustaining repetitive activity. In fact, pharmacological disruption of synchronized diastolic Ca release (by ryanodine) prevented aberrant APs; and vice versa, the inhibition of AP (by TTX or 0 Na, 0 Ca solution) de-synchronized diastolic Ca release. Taken together, these results suggest that a cyclical interaction between synchronized diastolic Ca release and AP forms a pathological rhythm generator that is involved in Ca-dependent atrial arrhythmias in CPVT.

Exploring susceptibility to atrial and ventricular arrhythmias resulting from remodeling of the passive electrical properties in the heart: a simulation approach

Under diseased conditions, remodeling of the cardiac tissue properties (“passive properties”) takes place; these are aspects of electrophysiological behavior that are not associated with active ion transport across cell membranes. Remodeling of the passive electrophysiological properties most often results from structural remodeling, such as gap junction down-regulation and lateralization, fibrotic growth infiltrating the myocardium, or the development of an infarct scar. Such structural remodeling renders atrial or ventricular tissue as a major substrate for arrhythmias. The current review focuses on these aspects of cardiac arrhythmogenesis. Due to the inherent complexity of cardiac arrhythmias, computer simulations have provided means to elucidate interactions pertinent to this spatial scale. Here we review the current state-of-the-art in modeling atrial and ventricular arrhythmogenesis as arising from the disease-induced changes in the passive tissue properties, as well as the contributions these modeling studies have made to our understanding of the mechanisms of arrhythmias in the heart. Because of the rapid advance of structural imaging methodologies in cardiac electrophysiology, we chose to present studies that have used such imaging methodologies to construct geometrically realistic models of cardiac tissue, or the organ itself, where the regional remodeling properties of the myocardium can be represented in a realistic way. We emphasize how the acquired knowledge can be used to pave the way for clinical applications of cardiac organ modeling under the conditions of structural remodeling.

Interrogation of living myocardium in multiple static deformation states with diffusion tensor and diffusion spectrum imaging

Diffusion tensor magnetic resonance imaging (MRI) reveals valuable insights into tissue histo-anatomy and microstructure, and has steadily gained traction in the cardiac community. Its wider use in small animal cardiac imaging in vivo has been constrained by its extreme sensitivity to motion, exaggerated by the high heart rates usually seen in rodents. Imaging of the isolated heart eliminates respiratory motion and, if conducted on arrested hearts, cardiac pulsation. This serves as an important intermediate step for basic and translational studies. However, investigating the micro-structural basis of cardiac deformation in the same heart requires observations in different deformation states.

Here, we illustrate the imaging of isolated rat hearts in three mechanical states mimicking diastole (cardioplegic arrest), left-ventricular (LV) volume overload (cardioplegic arrest plus LV balloon inflation), and peak systole (lithium-induced contracture). An optimised MRI-compatible Langendorff perfusion setup with the radio-frequency (RF) coil integrated into the wet chamber was developed for use in a 9.4T horizontal bore scanner. Signal-to-noise ratio improved significantly, by 75% compared to a previous design with external RF coil, and stability tests showed no significant changes in mean T₁, T₂ or LV wall thickness over a 170 min period. In contracture, we observed a significant reduction in mean fractional anisotropy from 0.32 ± 0.02 to 0.28 ± 0.02, as well as a significant rightward shift in helix angles with a decrease in the proportion of left-handed fibres, as referring to the locally prevailing cell orientation in the heart, from 24.9% to 23.3%, and an increase in the proportion of right-handed fibres from 25.5% to 28.4%. LV overload, in contrast, gave rise to a decrease in the proportion of left-handed fibres from 24.9% to 21.4% and an increase in the proportion of right-handed fibres from 25.5% to 26.0%. The modified perfusion and coil setup offers better performance and control over cardiac contraction states.

We subsequently performed high-resolution diffusion spectrum imaging (DSI) and 3D whole heart fibre tracking in fixed ex vivo rat hearts in slack state and contracture. As a model-free method, DSI augmented the measurements of water diffusion by also informing on multiple intra-voxel diffusion orientations and non-Gaussian diffusion. This enabled us to identify the transition from right- to left-handed fibres from the subendocardium to the subepicardium, as well as voxels in apical regions that were traversed by multiple fibres. We observed that both the mean generalised fractional anisotropy and mean kurtosis were lower in hearts in contracture compared to the slack state, by 23% and 9.3%, respectively. While its heavy acquisition burden currently limits the application of DSI in vivo, ongoing work in acceleration techniques may enable its use in live animals and patients. This would provide access to the as yet unexplored dimension of non-Gaussian diffusion that could serve as a highly sensitive marker of cardiac micro-structural integrity.

A compartmentalized mathematical model of the β1-adrenergic signaling system in mouse ventricular myocytes

The β₁-adrenergic signaling system plays an important role in the functioning of cardiac cells. Experimental data shows that the activation of this system produces inotropy, lusitropy, and chronotropy in the heart, such as increased magnitude and relaxation rates of [Ca²⁺]_i transients and contraction force, and increased heart rhythm. However, excessive stimulation of β₁-adrenergic receptors leads to heart dysfunction and heart failure. In this paper, a comprehensive, experimentally based mathematical model of the β₁-adrenergic signaling system for mouse ventricular myocytes is developed, which includes major subcellular functional compartments (caveolae, extracaveolae, and cytosol). The model describes biochemical reactions that occur during stimulation of β₁-adrenoceptors, changes in ionic currents, and modifications of Ca²⁺ handling system. Simulations describe the dynamics of major signaling molecules, such as cyclic AMP and protein kinase A, in different subcellular compartments; the effects of inhibition of phosphodiesterases on cAMP production; kinetics and magnitudes of phosphorylation of ion channels, transporters, and Ca²⁺ handling proteins; modifications of action potential shape and duration; magnitudes and relaxation rates of [Ca²⁺]_i transients; changes in intracellular and transmembrane Ca²⁺ fluxes; and [Na⁺]_i fluxes and dynamics. The model elucidates complex interactions of ionic currents upon activation of β₁-adrenoceptors at different stimulation frequencies, which ultimately lead to a relatively modest increase in action potential duration and significant increase in [Ca²⁺]_i transients. In particular, the model includes two subpopulations of the L-type Ca²⁺ channels, in caveolae and extracaveolae compartments, and their effects on the action potential and [Ca²⁺]_i transients are investigated. The presented model can be used by researchers for the interpretation of experimental data and for the developments of mathematical models for other species or for pathological conditions.

The holistic effects of acupuncture treatment

Traditional Chinese Medicine (TCM), as a complex medical science which reflects philosophical principles and embodies large dialectical thought, is used to place the human body into a large system for observation. Acupuncture as a vital part of TCM, has been practiced to treat various diseases and symptoms. However, acupuncture is also facing severe challenges resulted from insufficient modern scientific research. Nowadays, the holistic effects of acupuncture can be researched by some modern approaches, such as the systems biology and fMRI technique. It is believed that having a better understand will greatly promote acupuncture research and be beneficial to scientization and modernization of acupuncture.

Induced pluripotent stem cells as cardiac arrhythmic in vitro models and the impact for drug discovery

INTRODUCTION

The development of new antiarrhythmic agents is challenging and is hampered by high attrition rate of novel drug candidates. One of the reasons for this is limited predictability of existing preclinical models for drug assessment. Cardiomyocytes (CMs) derived from disease-specific induced pluripotent stem cells (iPSC) represent a novel in vitro cellular model of cardiac arrhythmias with an unprecedented potential for generating new mechanistic insight into disease pathophysiology and improving the process of drug development.

AREAS COVERED

This review outlines recent studies demonstrating the suitability and limitations of iPSC-derived CMs (iPS-CMs) for in vitro modeling inherited arrhythmias and drug testing. The authors focus on channelopathies and outline the properties of iPS-CMs, highlighting their utility and limitations for investigating the mechanism of cardiac arrhythmias and drug discovery.

EXPERT OPINION

The iPS-CMs represent a valuable addition to the already existing armamentarium of cardiac arrhythmic models. However, the superiority of iPS-CMs over other arrhythmia models has not yet been rigorously established and the limitations of the model must be overcome before its full potential for antiarrhythmic drug discovery can be realized. Nevertheless, iPS cell-based platforms hold a great potential for increasing our knowledge about cellular arrhythmia mechanisms and improving the drug discovery process.

Heart failure gene therapy: the path to clinical practice

Gene therapy, aimed at the correction of key pathologies being out of reach for conventional drugs, bears the potential to alter the treatment of cardiovascular diseases radically and thereby of heart failure. Heart failure gene therapy refers to a therapeutic system of targeted drug delivery to the heart that uses formulations of DNA and RNA, whose products determine the therapeutic classification through their biological actions. Among resident cardiac cells, cardiomyocytes have been the therapeutic target of numerous attempts to regenerate systolic and diastolic performance, to reverse remodeling and restore electric stability and metabolism. Although the concept to intervene directly within the genetic and molecular foundation of cardiac cells is simple and elegant, the path to clinical reality has been arduous because of the challenge on delivery technologies and vectors, expression regulation, and complex mechanisms of action of therapeutic gene products. Nonetheless, since the first demonstration of in vivo gene transfer into myocardium, there have been a series of advancements that have driven the evolution of heart failure gene therapy from an experimental tool to the threshold of becoming a viable clinical option. The objective of this review is to discuss the current state of the art in the field and point out inevitable innovations on which the future evolution of heart failure gene therapy into an effective and safe clinical treatment relies.

Clinical Applications of Omics Technologies on ZHENG Differentiation Research in Traditional Chinese Medicine

Traditional Chinese medicine (TCM) ZHENG is the basic concept of TCM theory. The effectiveness of TCM treatment depends on the accuracy of ZHENG differentiation. ZHENG differentiation, using the "four diagnostic methods," has the drawbacks of subjectivity and variability. Following development of omics technologies, which study the functional activities of human body from a system-wide perspective, it has been more and more applied in study of objectivity differentiating TCM ZHENG and understanding its biological mechanisms. This paper reviewed the literatures of clinical TCM ZHENG differentiation researches, underlying omics technologies, and indicated the increased trends of related articles with four kinds of omics technologies, including genomics, transcriptomics, proteomics and metabolomics, and the correlations between ZHENG differentiation and findings in omics studies. Moreover, the paper summarized the typical omics application in common studied diseases and TCM ZHENGs and discussed the main problems and countermeasure of ZHENG differentiation researches. The work here may provide a reference for further research of TCM ZHENG differentiation using omics technologies.

Mining disease fingerprints from within genetic pathways

Mining biological networks can be an effective means to uncover system level knowledge out of micro level associations, such as encapsulated in genetic pathways. Analysis of human disease genetic pathways can lead to the identification of major mechanisms that may underlie disorders at an abstract functional level. The focus of this study was to develop an approach for structural pattern analysis and classification of genetic pathways of diseases. A probabilistic model was developed to capture characteristic components ('fingerprints') of functionally annotated pathways. A probability estimation procedure of this model searched for fingerprints in each disease pathway while improving probability estimates of model parameters. The approach was evaluated on data from the Kyoto Encyclopedia of Genes and Genomes (consisting of 56 pathways across seven disease categories). Based on the achieved average classification accuracy of up to ~77%, the findings suggest that these fingerprints may be used for classification and discovery of genetic pathways.

Systems biology and cardiac arrhythmias

During the past few years, the development of effective, empirical technologies for treatment of cardiac arrhythmias has exceeded the pace at which detailed knowledge of the underlying biology has accumulated. As a result, although some clinical arrhythmias can be cured with techniques such as catheter ablation, drug treatment and prediction of the risk of sudden death remain fairly primitive. The identification of key candidate genes for monogenic arrhythmia syndromes shows that to bring basic biology to the clinic is a powerful approach. Increasingly sophisticated experimental models and methods of measurement, including stem cell-based models of human cardiac arrhythmias, are being deployed to study how perturbations in several biologic pathways can result in an arrhythmia-prone heart. The biology of arrhythmia is largely quantifiable, which allows for systematic analysis that could transform treatment strategies that are often still empirical into management based on molecular evidence.

Systems biology of meridians, acupoints, and chinese herbs in disease

Meridians, acupoints, and Chinese herbs are important components of traditional Chinese medicine (TCM). They have been used for disease treatment and prevention and as alternative and complementary therapies. Systems biology integrates omics data, such as transcriptional, proteomic, and metabolomics data, in order to obtain a more global and complete picture of biological activity. To further understand the existence and functions of the three components above, we reviewed relevant research in the systems biology literature and found many recent studies that indicate the value of acupuncture and Chinese herbs. Acupuncture is useful in pain moderation and relieves various symptoms arising from acute spinal cord injury and acute ischemic stroke. Moreover, Chinese herbal extracts have been linked to wound repair, the alleviation of postmenopausal osteoporosis severity, and anti-tumor effects, among others. Different acupoints, variations in treatment duration, and herbal extracts can be used to alleviate various symptoms and conditions and to regulate biological pathways by altering gene and protein expression. Our paper demonstrates how systems biology has helped to establish a platform for investigating the efficacy of TCM in treating different diseases and improving treatment strategies.

Markov models of use-dependence and reverse use-dependence during the mouse cardiac action potential

The fast component of the cardiac transient outward current, I_Ktof, is blocked by a number of drugs. The major molecular bases of I_Ktof are Kv4.2/Kv4.3 voltage-gated potassium channels. Drugs with similar potencies but different blocking mechanisms have differing effects on action potential duration (APD). We used in silico analysis to determine the effect of I_Ktof-blocking drugs with different blocking mechanisms on mouse ventricular myocytes. We used our existing mouse model of the action potential, and developed 4 new Markov formulations for I_Ktof, I_Ktos, I_Kur, I_Ks. We compared effects of theoretical I_Ktof-specific channel blockers: (1) a closed state, and (2) an open channel blocker. At concentrations lower or close to IC₅₀, the drug which bound to the open state always had a much greater effect on APD than the drug which bound to the closed state. At concentrations much higher than IC₅₀, both mechanisms had similar effects at very low pacing rates. However, an open state binding drug had a greater effect on APD at faster pacing rates, particularly around 10 Hz. In summary, our data indicate that drug effects on APD are strongly dependent not only on IC₅₀, but also on the drug binding state.

At the heart of computational modelling

The link between experimental data and biophysically based mathematical models is key to computational simulation meeting its potential to provide physiological insight. However, despite the importance of this link, scrutiny and analysis of the processes by which models are parameterised from data are currently lacking. While this situation is common to many areas of physiological modelling, to provide a concrete context, we use examples drawn from detailed models of cardiac electro-mechanics. Using this biophysically detailed cohort of models we highlight the specific issues of model parameterization and propose this process can be separated into three stages: observation, fitting and validation. Finally, future research challenges and directions in this area are discussed.

Simulation of the undiseased human cardiac ventricular action potential: model formulation and experimental validation

Cellular electrophysiology experiments, important for understanding cardiac arrhythmia mechanisms, are usually performed with channels expressed in non myocytes, or with non-human myocytes. Differences between cell types and species affect results. Thus, an accurate model for the undiseased human ventricular action potential (AP) which reproduces a broad range of physiological behaviors is needed. Such a model requires extensive experimental data, but essential elements have been unavailable. Here, we develop a human ventricular AP model using new undiseased human ventricular data: Ca(2+) versus voltage dependent inactivation of L-type Ca(2+) current (I(CaL)); kinetics for the transient outward, rapid delayed rectifier (I(Kr)), Na(+)/Ca(2+) exchange (I(NaCa)), and inward rectifier currents; AP recordings at all physiological cycle lengths; and rate dependence and restitution of AP duration (APD) with and without a variety of specific channel blockers. Simulated APs reproduced the experimental AP morphology, APD rate dependence, and restitution. Using undiseased human mRNA and protein data, models for different transmural cell types were developed. Experiments for rate dependence of Ca(2+) (including peak and decay) and intracellular sodium ([Na(+)](i)) in undiseased human myocytes were quantitatively reproduced by the model. Early afterdepolarizations were induced by I(Kr) block during slow pacing, and AP and Ca(2+) alternans appeared at rates >200 bpm, as observed in the nonfailing human ventricle. Ca(2+)/calmodulin-dependent protein kinase II (CaMK) modulated rate dependence of Ca(2+) cycling. I(NaCa) linked Ca(2+) alternation to AP alternans. CaMK suppression or SERCA upregulation eliminated alternans. Steady state APD rate dependence was caused primarily by changes in [Na(+)](i), via its modulation of the electrogenic Na(+)/K(+) ATPase current. At fast pacing rates, late Na(+) current and I(CaL) were also contributors. APD shortening during restitution was primarily dependent on reduced late Na(+) and I(CaL) currents due to inactivation at short diastolic intervals, with additional contribution from elevated I(Kr) due to incomplete deactivation.

Human pluripotent stem cells in pharmacological and toxicological screening: new perspectives for personalized medicine

Human stem cells provide an important novel tool for generating in vitro pharmacological and toxicological test systems. In the development of new targeted therapies, as well as in critical safety issues, including hepato-, neuro- and cardio-toxicity, animal-based tests are mostly unsatisfactory, whereas the use of in vitro model systems is limited by the unavailability of relevant human tissues. Human embryonic stem cell lines may fill this gap and offer an advantage over primary cultures as well as tissue-derived (adult) stem cells. Human embryonic stem cells represent an unlimited source for the production of differentiated somatic progenies and allow various stable genetic manipulations. As a new opening in personalized medicine test systems, the generation of induced pluripotent stem cell lines and their derivatives can provide patient- and disease-specific cellular assays for drug development and safety assessments. This article reviews promising human stem cell applications in pharmacological and toxicological screenings, focusing on the implications for personalized medicine.

Simulating human cardiac electrophysiology on clinical time-scales

In this study, the feasibility of conducting in silico experiments in near-realtime with anatomically realistic, biophysically detailed models of human cardiac electrophysiology is demonstrated using a current national high-performance computing facility. The required performance is achieved by integrating and optimizing load balancing and parallel I/O, which lead to strongly scalable simulations up to 16,384 compute cores. This degree of parallelization enables computer simulations of human cardiac electrophysiology at 240 times slower than real time and activation times can be simulated in approximately 1 min. This unprecedented speed suffices requirements for introducing in silico experimentation into a clinical workflow.

Are computational models of any use to psychiatry?

Mathematically rigorous descriptions of key hypotheses and theories are becoming more common in neuroscience and are beginning to be applied to psychiatry. In this article two fictional characters, Dr. Strong and Mr. Micawber, debate the use of such computational models (CMs) in psychiatry. We present four fundamental challenges to the use of CMs in psychiatry: (a) the applicability of mathematical approaches to core concepts in psychiatry such as subjective experiences, conflict and suffering; (b) whether psychiatry is mature enough to allow informative modelling; (c) whether theoretical techniques are powerful enough to approach psychiatric problems; and (d) the issue of communicating clinical concepts to theoreticians and vice versa. We argue that CMs have yet to influence psychiatric practice, but that they help psychiatric research in two fundamental ways: (a) to build better theories integrating psychiatry with neuroscience; and (b) to enforce explicit, global and efficient testing of hypotheses through more powerful analytical methods. CMs allow the complexity of a hypothesis to be rigorously weighed against the complexity of the data. The paper concludes with a discussion of the path ahead. It points to stumbling blocks, like the poor communication between theoretical and medical communities. But it also identifies areas in which the contributions of CMs will likely be pivotal, like an understanding of social influences in psychiatry, and of the co-morbidity structure of psychiatric diseases.

Sodium channel molecular conformations and antiarrhythmic drug affinity

Class I cardiac antiarrhythmic drugs, for example, lidocaine, mexiletine, flecainide, quinidine, and procainamide, continue to play an important role in the therapy for cardiac arrhythmias because of the presence of use-dependent block. Lidocaine, as well as related drugs such as mepivacaine, bupivacaine, and cocaine, also belong to the class of medications referred to as local anesthetics. In this review, we will consider lidocaine as the prototypical antiarrhythmic drug because it continues to be widely used both as an antiarrhythmic drug (first used as an antiarrhythmic drug in 1950) as well as a local anesthetic agent. Both of these clinical uses depend upon block of sodium current (I(Na)), but it is the presence of use-dependent I(Na) block, that is, an increasing amount of block at faster heart rates, which enables a local anesthetic agent to be a useful antiarrhythmic drug. Although many early studies investigated the action of antiarrhythmic drugs on Na currents, the availability of site-directed mutant Na channels has enabled for major advances in understanding their mechanisms of action based upon molecular conformations of the Na channel.

S100A1 gene therapy for heart failure: a novel strategy on the verge of clinical trials

Representing the common endpoint of various cardiovascular disorders, heart failure (HF) shows a dramatically growing prevalence. As currently available therapeutic strategies are not capable of terminating the progress of the disease, HF is still associated with a poor clinical prognosis. Among the underlying molecular mechanisms, the loss of cardiomyocyte Ca(2+) cycling integrity plays a key role in the pathophysiological development and progression of the disease. The cardiomyocyte EF-hand Ca(2+) sensor protein S100A1 emerged as a regulator both of sarcoplasmic reticulum (SR), sarcomere and mitochondrial function implicating a significant role in cardiac physiology and dysfunction. In this review, we aim to recapitulate the translation of S100A1-based investigation from first clinical observations over basic research experiments back to a near-clinical setting on the verge of clinical trials today. We also address needs for further developments towards "second-generation" gene therapy and discuss the therapeutic potential of S100A1 gene therapy for HF as a promising novel strategy for future cardiologists. This article is part of a Special Section entitled "Special Section: Cardiovascular Gene Therapy".

Short-term disruption in regional left ventricular electrical conduction patterns increases interstitial matrix metalloproteinase activity

Increased matrix metalloproteinase (MMP) abundance occurs with adverse left ventricular (LV) remodeling in a number of cardiac disease states, including those induced by long-standing arrhythmias. However, whether regionally contained aberrant electrical activation of the LV, with consequent dyskinesia, alters interstitial MMP activation remained unknown. Electrical activation of the LV of pigs (n = 10, 30-35 kg) was achieved by pacing (150 beats/min) at left atrial and LV sites such that normal atrioventricular activation (60 min) was followed by regional early LV activation for 60 min within 1.5 cm of the paced site and restoration of normal atrioventricular pacing for 120 min. Regional shortening (piezoelectric crystals) and interstitial MMP activity (microdialysis with MMP fluorogenic substrate) at the LV pacing site and a remote LV site were monitored at 30-min intervals. During aberrant electrical stimulation, interstitial MMP activity at the paced site was increased (122 +/- 4%) compared with the remote region (100%, P < 0.05). Restoration of atrioventricular pacing after the 60-min period of aberrant electrical activation normalized segmental shortening (8.5 +/- 0.4%), but MMP activity remained elevated (121 +/- 6%, P < 0.05). This study demonstrates that despite the restoration of mechanical function, disturbances in electrical conduction, in and of itself, can cause acute increases in regional in vivo MMP activation and, therefore, contribute to myocardial remodeling.

Why does troponin I have so many phosphorylation sites? Fact and fancy

We discuss a current controversy regarding the relative role of phosphorylation sites on cardiac troponin I (cTnI) (Fig. 1) in physiological and patho-physiological cardiac function. Studies with mouse models and in vitro studies indicate that multi-site phosphorylations are involved in both control of maximum tension and sarcomeric responsiveness to Ca(2+). Thus one hypothesis is that cardiac function reflects a balance of cTnI phosphorylations and a tilt in this balance may be maladaptive in acquired and genetic disorders of the heart. Studies on human heart samples taken mainly at end-stage heart failure, and in depth proteomic analysis of human and rat heart samples demonstrate that Ser23/Ser24 are the major and perhaps the only sites likely to be relevant to control cardiac function. Thus functional significance of Ser23/Ser24 phosphorylation is taken as fact, whereas the function of some other sites is treated as fancy. Maybe the extremes will meet: in any case we both agree that further work needs to be carried out with relatively large mammals and with determination of the time course of changes in phosphorylation to identify transient modifications that may be relevant at a beat-to-beat basis. Moreover, we agree that the changes and effects of cTnI phosphorylation need to be fully integrated into the effects of other phosphorylations in the cardiac myocyte.

Current progress in patient-specific modeling

We present a survey of recent advancements in the emerging field of patient-specific modeling (PSM). Researchers in this field are currently simulating a wide variety of tissue and organ dynamics to address challenges in various clinical domains. The majority of this research employs three-dimensional, image-based modeling techniques. Recent PSM publications mostly represent feasibility or preliminary validation studies on modeling technologies, and these systems will require further clinical validation and usability testing before they can become a standard of care. We anticipate that with further testing and research, PSM-derived technologies will eventually become valuable, versatile clinical tools.

Cardiac repolarization: insights from mathematical modeling and electrocardiographic imaging (ECGI)

Cardiac repolarization is a complex rate dependent process. At the cellular level, it depends on a delicate dynamic balance of ion channel currents. At the heart level, it is spatially heterogeneous, leading to spatial gradients of potential and excitability. This article provides insights into the molecular mechanisms of the delayed rectifiers I(Kr) (rapid) and I(Ks) (slow) that underlie effective function of these channels as repolarizing currents during the cardiac action potential (AP). We also provide noninvasive images of cardiac repolarization in humans. Methodologically, computational biology is used to simulate ion channel function and to incorporate it into a model of the cardiac cell. ECG imaging (ECGI) is applied to normal subjects and Wolff-Parkinson-White (WPW) patients to obtain epicardial maps of repolarization. The simulations demonstrate that I(Kr) and I(Ks) generate their peak current late during the AP, where they effectively participate in repolarization. I(Kr) maximizes the current by fast inactivation and gradual recovery during the AP. I(Ks) does so by generating an available reserve of channels in closed states from which the channels can open rapidly. ECGI shows that in the human heart, normal repolarization epicardial potential maps are static with 40 ms dispersion between RV and LV. In WPW, ECGI located the accessory pathway(s) and showed a large base-to-apex repolarization gradient that resolved to normal one month post-ablation, demonstrating presence of "cardiac memory". We conclude that computational biology can provide a mechanistic link across scales, from the molecular functioning of ion channels to the cellular AP. ECGI can noninvasively image human cardiac repolarization and its alteration by disease and interventions. This property makes it a potential tool for noninvasive risk stratification and evaluation of treatment by drugs and devices.

Towards predictive modelling of the electrophysiology of the heart

The simulation of cardiac electrical function is an example of a successful integrative multiscale modelling approach that is directly relevant to human disease. Today we stand at the threshold of a new era, in which anatomically detailed, tomographically reconstructed models are being developed that integrate from the ion channel to the electromechanical interactions in the intact heart. Such models hold high promise for interpretation of clinical and physiological measurements, for improving the basic understanding of the mechanisms of dysfunction in disease, such as arrhythmias, myocardial ischaemia and heart failure, and for the development and performance optimization of medical devices. The goal of this article is to present an overview of current state-of-art advances towards predictive computational modelling of the heart as developed recently by the authors of this article. We first outline the methodology for constructing electrophysiological models of the heart. We then provide three examples that demonstrate the use of these models, focusing specifically on the mechanisms for arrhythmogenesis and defibrillation in the heart. These include: (1) uncovering the role of ventricular structure in defibrillation; (2) examining the contribution of Purkinje fibres to the failure of the shock; and (3) using magnetic resonance imaging reconstructed heart models to investigate the re-entrant circuits formed in the presence of an infarct scar.