Cardiac electrical dynamics: maximizing dynamical heterogeneity

Short-Long Heart Rate Variation Increases Dispersion of Action Potential Duration in Long QT Type 2 Transgenic Rabbit Model

The initiation of polymorphic ventricular tachycardia in long QT syndrome type 2 (LQT2) has been associated with a characteristic ECG pattern of short-long RR intervals. We hypothesize that this characteristic pattern increases APD dispersion in LQT2, thereby promoting arrhythmia. We investigated APD dispersion and its dependence on two previous cycle lengths (CLs) in transgenic rabbit models of LQT2, LQT1, and their littermate controls (LMC) using random stimulation protocols. The results show that the short-long RR pattern was associated with a larger APD dispersion in LQT2 but not in LQT1 rabbits. The multivariate analyses of APD as a function of two previous CLs (APD_n = C + α₁CL_n−1 + α₂CL_n−2) showed that α₁ (APD restitution slope) is largest and heterogeneous in LQT2 but uniform in LQT1, enhancing APD dispersion under long CL_n−1 in LQT2. The α₂ (short-term memory) was negative in LQT2 while positive in LQT1, and the spatial pattern of α₁ was inversely correlated to α₂ in LQT2, which explains why a short-long combination causes a larger APD dispersion in LQT2 but not in LQT1 rabbits. In conclusion, short-long RR pattern increased APD dispersion only in LQT2 rabbits through heterogeneous APD restitution and the short-term memory, underscoring the genotype-specific triggering of arrhythmias in LQT syndrome.

Laser-Etched Designs for Molding Hydrogel-Based Engineered Tissues

Rapid prototyping and fabrication of elastomeric molds for sterile culture of engineered tissues allow for the development of tissue geometries that can be tailored to different in vitro applications and customized as implantable scaffolds for regenerative medicine. Commercially available molds offer minimal capabilities for adaptation to unique conditions or applications versus those for which they are specifically designed. Here we describe a replica molding method for the design and fabrication of poly(dimethylsiloxane) (PDMS) molds from laser-etched acrylic negative masters with ∼0.2 mm resolution. Examples of the variety of mold shapes, sizes, and patterns obtained from laser-etched designs are provided. We use the patterned PDMS molds for producing and culturing engineered cardiac tissues with cardiomyocytes derived from human-induced pluripotent stem cells. We demonstrate that tight control over tissue morphology and anisotropy results in modulation of cell alignment and tissue-level conduction properties, including the appearance and elimination of reentrant arrhythmias, or circular electrical activation patterns. Techniques for handling engineered cardiac tissues during implantation in vivo in a rat model of myocardial infarction have been developed and are presented herein to facilitate development and adoption of surgical techniques for use with hydrogel-based engineered tissues. In summary, the method presented herein for engineered tissue mold generation is straightforward and low cost, enabling rapid design iteration and adaptation to a variety of applications in tissue engineering. Furthermore, the burden of equipment and expertise is low, allowing the technique to be accessible to all.

Effect of mental stress on dynamic electrophysiological properties of the endocardium and epicardium in humans

BACKGROUND

Striking temporal associations exist between ventricular arrhythmia and acute mental stress, for example, during natural disasters, or defibrillator shocks associated with stressful events. We hypothesized that electrophysiological changes in response to mental stress may be exaggerated at short coupling intervals and hence relevant to arrhythmia initiation.

OBJECTIVE

The aim of this study was to determine the dynamic response in human electrophysiology during mental stress.

METHODS

Patients with normal hearts and supraventricular tachycardia underwent electrophysiological studies avoiding sedation. Conditions of relaxation and stress were induced with standardized psychometric protocols (mental arithmetic and anger recall) during decremental S1S2 right ventricular (RV) pacing. Unipolar electrograms were acquired simultaneously from the RV endocardium, left ventricular (LV) endocardium (LV endo), and epicardium (LV epi), and activation-recovery intervals (ARIs) computed.

RESULTS

Twelve patients ( 9 women; median age 34 years) were studied. During stress, effective refractory period (ERP) reduced from 228 ± 23 to 221 ± 21 ms (P < .001). ARIs reduced during mental stress (P < .001), with greater reductions in LV endocardium than in the epicardium or RV endocardium (LV endo -8 ms; LV epi -5 ms; RV endo -4 ms; P < .001). Mental stress depressed the entire electrical restitution curve, with minimal effect on slope. A substantial reduction in minimal ARIs on the restitution curve in LV endo occurred, commensurate with the reduction in ERP (LV endo ARI 195 ± 31 ms at rest to 182 ± 32 ms during mental stress; P < .001). Dispersion of repolarization increased sharply at coupling intervals approaching ERP during stress but not at rest.

CONCLUSION

Mental stress induces significant electrophysiological changes. The increase in dispersion of repolarization at short coupling intervals may be relevant to observed phenomena of arousal-associated arrhythmia.

Absolute beat-to-beat variability and instability parameters of ECG intervals: biomarkers for predicting ischaemia-induced ventricular fibrillation

BACKGROUND AND PURPOSE

Predicting lethal arrhythmia liability from beat-to-beat variability and instability (BVI) of the ECG intervals is a useful technique in drug assessment. Most investigators use only arrhythmia-free ECGs for this. Recently, it was shown that drug-induced torsades de pointes (TdP) liability can be predicted more accurately from BVI measured irrespective of rhythm, even during arrhythmias (absolute BVI). The present study tested the broader applicability of this assessment by examining whether absolute BVI parameters predict another potential lethal arrhythmia, ischaemia-induced ventricular fibrillation (VF).

EXPERIMENTAL APPROACH

Langendorff-perfused rat hearts were subjected to regional ischaemia for 15 min. Absolute BVI parameters were derived from ECG intervals measured in 40 consecutive ventricular complexes (irrespective of rhythm) immediately preceding VF onset and compared with time-matched values in hearts not expressing VF.

KEY RESULTS

Increased frequency of non-sinus beats and 'R on T' arrhythmic beats, shortened mean RR and electrical diastolic intervals, and increased BVI of cycle length and repolarization predicted VF occurrence. Absolute BVI parameters that quantify variability of repolarization (e.g. 'short-term variability' of QT interval) had the best predictive power with high sensitivity and specificity. In contrast, VF was not predicted by any BVI parameter derived from the last arrhythmia-free interlude before VF.

CONCLUSIONS AND IMPLICATIONS

The novel absolute BVI parameters that predicted TdP in rabbit also predict ischaemia-induced VF in rat, indicating a diagnostic and mechanistic congruence. Repolarization inhomogeneity represents a pivotal biomarker of ischaemia-induced VF. The newly validated biomarkers could serve as surrogates for VF in pre-clinical drug investigations.

Effects of human atrial ionic remodelling by β-blocker therapy on mechanisms of atrial fibrillation: a computer simulation

AIMS

Atrial anti-arrhythmic effects of β-adrenoceptor antagonists (β-blockers) may involve both a suppression of pro-arrhythmic effects of catecholamines, and an adaptational electrophysiological response to chronic β-blocker use; so-called 'pharmacological remodelling'. In human atrium, such remodelling decreases the transient outward (Ito) and inward rectifier (IK1) K(+) currents, and increases the cellular action potential duration (APD) and effective refractory period (ERP). However, the consequences of these changes on mechanisms of genesis and maintenance of atrial fibrillation (AF) are unknown. Using mathematical modelling, we tested the hypothesis that the long-term adaptational decrease in human atrial Ito and IK1 caused by chronic β-blocker therapy, i.e. independent of acute electrophysiological effects of β-blockers, in an otherwise un-remodelled atrium, could suppress AF.

METHODS AND RESULTS

Contemporarily, biophysically detailed human atrial cell and tissue models were used to investigate effects of the β-blocker-based pharmacological remodelling. Chronic β-blockade remodelling prolonged atrial cell APD and ERP. The incidence of small amplitude APD alternans in the CRN model was reduced. At the 1D tissue level, β-blocker remodelling decreased the maximum pacing rate at which APs could be conducted. At the three-dimensional organ level, β-blocker remodelling reduced the life span of re-entry scroll waves.

CONCLUSION

This study improves our understanding of the electrophysiological mechanisms of AF suppression by chronic β-blocker therapy. Atrial fibrillation suppression may involve a reduced propensity for maintenance of re-entrant excitation waves, as a consequence of increased APD and ERP.

Bridging the gap between computation and clinical biology: validation of cable theory in humans

Introduction: Computerized simulations of cardiac activity have significantly contributed to our understanding of cardiac electrophysiology, but techniques of simulations based on patient-acquired data remain in their infancy. We sought to integrate data acquired from human electrophysiological studies into patient-specific models, and validated this approach by testing whether electrophysiological responses to sequential premature stimuli could be predicted in a quantitatively accurate manner.

Methods: Eleven patients with structurally normal hearts underwent electrophysiological studies. Semi-automated analysis was used to reconstruct activation and repolarization dynamics for each electrode. This S₂ extrastimuli data was used to inform individualized models of cardiac conduction, including a novel derivation of conduction velocity restitution. Activation dynamics of multiple premature extrastimuli were then predicted from this model and compared against measured patient data as well as data derived from the ten-Tusscher cell-ionic model.

Results: Activation dynamics following a premature S₃ were significantly different from those after an S₂. Patient specific models demonstrated accurate prediction of the S₃ activation wave, (Pearson's R² = 0.90, median error 4%). Examination of the modeled conduction dynamics allowed inferences into the spatial dispersion of activation delay. Further validation was performed against data from the ten-Tusscher cell-ionic model, with our model accurately recapitulating predictions of repolarization times (R² = 0.99).

Conclusions: Simulations based on clinically acquired data can be used to successfully predict complex activation patterns following sequential extrastimuli. Such modeling techniques may be useful as a method of incorporation of clinical data into predictive models.

A novel surface electrocardiogram-based marker of ventricular arrhythmia risk in patients with ischemic cardiomyopathy

Background

Better sudden cardiac death risk markers are needed in ischemic cardiomyopathy (ICM). Increased heterogeneity of electrical restitution is an important mechanism underlying the risk of ventricular arrhythmia (VA). Our aim was to develop and test a novel quantitative surface electrocardiogram–based measure of VA risk in patients with ICM: the Regional Restitution Instability Index (R2I2).

Methods and Results

R2I2, the mean of the standard deviation of residuals from the mean gradient for each ECG lead at a range of diastolic intervals, was measured retrospectively from high-resolution 12-lead ECGs recorded during an electrophysiology study. Patient groups were as follows: Study group, 26 patients with ICM being assessed for implantable defibrillator; Control group, 29 patients with supraventricular tachycardia undergoing electrophysiology study; and Replication group, 40 further patients with ICM. R2I2 was significantly higher in the Study patients than in Controls (mean ± standard error of the mean: 1.09±0.06 versus 0.63±0.04, P<0.001). Over a median follow-up period of 23 months, 6 of 26 Study group patients had VA or death. R2I2 predicted VA or death independently of demographic factors, electrophysiology study result, left ventricular ejection fraction, or QRS duration (Cox model, P=0.029). R2I2 correlated with peri-infarct zone as assessed by cardiac magnetic resonance imaging (r=0.51, P=0.024). The findings were replicated in the Replication group: R2I2 was significantly higher in 11 of 40 Replication patients experiencing VA (1.18±0.10 versus 0.92±0.05, P=0.019). In combined analysis of ICM cohorts, R2I2 ≥1.03 identified subjects with significantly higher risk of VA or death (43%) compared with R2I2 <1.03 (11%) (P=0.004).

Conclusions

In this pilot study, we have developed a novel VA risk marker, R2I2, and have shown that it correlated with a structural measure of arrhythmic risk and predicted risk of VA or death in patients with ICM. R2I2 may improve risk stratification and merits further evaluation. (J Am Heart Assoc. 2012;1:e001552 doi: 10.1161/JAHA.112.001552.)

Computational approaches to understand cardiac electrophysiology and arrhythmias

Cardiac rhythms arise from electrical activity generated by precisely timed opening and closing of ion channels in individual cardiac myocytes. These impulses spread throughout the cardiac muscle to manifest as electrical waves in the whole heart. Regularity of electrical waves is critically important since they signal the heart muscle to contract, driving the primary function of the heart to act as a pump and deliver blood to the brain and vital organs. When electrical activity goes awry during a cardiac arrhythmia, the pump does not function, the brain does not receive oxygenated blood, and death ensues. For more than 50 years, mathematically based models of cardiac electrical activity have been used to improve understanding of basic mechanisms of normal and abnormal cardiac electrical function. Computer-based modeling approaches to understand cardiac activity are uniquely helpful because they allow for distillation of complex emergent behaviors into the key contributing components underlying them. Here we review the latest advances and novel concepts in the field as they relate to understanding the complex interplay between electrical, mechanical, structural, and genetic mechanisms during arrhythmia development at the level of ion channels, cells, and tissues. We also discuss the latest computational approaches to guiding arrhythmia therapy.

Unstable QT interval dynamics precedes ventricular tachycardia onset in patients with acute myocardial infarction: a novel approach to detect instability in QT interval dynamics from clinical ECG

BACKGROUND

Instability in ventricular repolarization in the presence of premature activations (PA) plays an important role in arrhythmogenesis. However, such instability cannot be detected clinically. This study developed a methodology for detecting QT interval (QTI) dynamics instability from the ECG and explored the contribution of PA and QTI instability to ventricular tachycardia (VT) onset.

METHODS AND RESULTS

To examine the contribution of PAs and QTI instability to VT onset, ECGs of 24 patients with acute myocardial infarction, 12 of whom had sustained VT (VT) and 12 nonsustained VT (NSVT), were used. From each patient ECG, 2 10-minute-long ECG recordings were extracted, 1 right before VT onset (onset epoch) and 1 at least 1 hour before it (control epoch). To ascertain how PA affects QTI dynamics stability, pseudo-ECGs were calculated from an MRI-based human ventricular model. Clinical and pseudo-ECGs were subdivided into 1-minute recordings (minECGs). QTI dynamics stability of each minECG was assessed with a novel approach. Frequency of PAs (f(PA)) and the number of minECGs with unstable QTI dynamics (N(us)) were determined for each patient. In the VT group, f(PA) and N(us) of the onset epoch were larger than in control. Positive regression relationships between f(PA) and N(us) were identified in both groups. The simulations showed that both f(PA) and the PA degree of prematurity contribute to QTI dynamics instability.

CONCLUSIONS

Increased PA frequency and QTI dynamics instability precede VT onset in patients with acute myocardial infarction, as determined by novel methodology for detecting instability in QTI dynamics from clinical ECGs.

A novel methodology for assessing the bounded-input bounded-output instability in QT interval dynamics: application to clinical ECG with ventricular tachycardia

The goal of this paper is to present a new methodology for assessing the bounded-input bounded-output (BIBO) stability in QT interval (QTI) dynamics from clinical ECG. The ECG recordings were collected from 15 patients who experienced ventricular tachycardia (VT). Ten-minute-long ECG recordings extracted immediately before the onset of a chosen VT, one per patient, were assembled into a VT group, while the control group comprised 10-min-long ECGs extracted 1 h before VT onset and at least 1 h after any prior arrhythmic event. Each 10-min recording was subdivided into 1-min ECG recordings (minECGs). The QTI dynamics of each minECG was defined as a function of several prior QTIs and RR intervals; the BIBO stability of this function was then assessed in the z -domain. The number of minECGs with unstable QTI dynamics (N (us)) and the frequency of premature activations (PA),  f (PA) , were counted for each ECG recording and were compared between the VT and control groups. The results show that the present methodology successfully captured the instability in QTI dynamics leading to VT onset in the studied population. Significantly larger N (us) was found in the VT group compared against the control and a positive correlation between N (us) and f (PA) was identified in both groups.

Kinematics of spiral waves under feedback-related spatial gradients

The kinematics of spiral waves with artificially constructed spatial excitability is numerically investigated in the Oregonator model. On an assumption that the rotation center of spiral's tip drifts at angle δ to the direction of a local gradient, a kinematic formula of motion of spiral's tip is derived. To test the formula, we have presented two forms of feedback-related spatial fields with radial gradients (RGs) and concentric circular gradients (CGs) both centering on a reference point. It is found that both rigidly rotating and meandering spiral waves are attracted to the reference point of an inward RG and a clockwise CG perturbation but moved away from it under an outward RG and a counterclockwise CG. Simulations of the drift-velocity formulas provide a quantitative testing of the numerical results.

Boundary-induced reentry in homogeneous excitable tissue

Heterogeneity of cardiac electrical properties can lead to heart rhythm disorders. Numerical studies have shown that stimuli chosen to maximize dynamic heterogeneity terminate wave propagation. However, experimental investigations suggest that similar sequences induce fragmentation of the wave fronts, rather than complete wave block. In this paper we show that an insulating boundary in an otherwise homogeneous medium can disrupt dynamically induced wave block by breaking a symmetry in the spatial pattern of action potential duration, leading to unidirectional block and reentrant activation.