The influence of fibre orientation, extracted from different segments of the human left ventricle, on the activation and repolarization sequence: a simulation study

Computational analysis of arrhythmogenesis in KCNH2 T618I mutation-associated short QT syndrome and the pharmacological effects of quinidine and sotalol

Short QT syndrome (SQTS) is a rare but dangerous genetic disease. In this research, we conducted a comprehensive in silico investigation into the arrhythmogenesis in KCNH2 T618I-associated SQTS using a multi-scale human ventricle model. A Markov chain model of IKr was developed firstly to reproduce the experimental observations. It was then incorporated into cell, tissue, and organ models to explore how the mutation provided substrates for ventricular arrhythmias. Using this T618I Markov model, we explicitly revealed the subcellular level functional alterations by T618I mutation, particularly the changes of ion channel states that are difficult to demonstrate in wet experiments. The following tissue and organ models also successfully reproduced the changed dynamics of reentrant spiral waves and impaired rate adaptions in hearts of T618I mutation. In terms of pharmacotherapy, we replicated the different effects of a drug under various conditions using identical mathematical descriptions for drugs. This study not only simulated the actions of an effective drug (quinidine) at various physiological levels, but also elucidated why the IKr inhibitor sotalol failed in SQT1 patients through profoundly analyzing its mutation-dependent actions.

Spatial distribution of physiologic 12-lead QRS complex

The normal physiologic range of QRS complex duration spans between 80 and 125 ms with known differences between females and males which cannot be explained by the anatomical variations of heart sizes. To investigate the reasons for the sex differences as well as for the wide range of normal values, a technology is proposed based on the singular value decomposition and on the separation of different orthogonal components of the QRS complex. This allows classification of the proportions of different components representing the 3-dimensional representation of the electrocardiographic signal as well as classification of components that go beyond the 3-dimensional representation and that correspond to the degree of intricate convolutions of the depolarisation sequence. The technology was applied to 382,019 individual 10-s ECG samples recorded in 639 healthy subjects (311 females and 328 males) aged 33.8 ± 9.4 years. The analyses showed that QRS duration was mainly influenced by the proportions of the first two orthogonal components of the QRS complex. The first component demonstrated statistically significantly larger proportion of the total QRS power (expressed by the absolute area of the complex in all independent ECG leads) in females than in males (64.2 ± 11.6% vs 59.7 ± 11.9%, p < 0.00001—measured at resting heart rate of 60 beats per minute) while the second component demonstrated larger proportion of the QRS power in males compared to females (33.1 ± 11.9% vs 29.6 ± 11.4%, p < 0.001). The analysis also showed that the components attributable to localised depolarisation sequence abnormalities were significantly larger in males compared to females (2.85 ± 1.08% vs 2.42 ± 0.87%, p < 0.00001). In addition to the demonstration of the technology, the study concludes that the detailed convolution of the depolarisation waveform is individual, and that smoother and less intricate depolarisation propagation is the mechanism likely responsible for shorter QRS duration in females.

Wide-area low-energy surface stimulation of large mammalian ventricular tissue

The epicardial and endocardial surfaces of the heart are attractive targets to administer antiarrhythmic electrotherapies. Electrically stimulating wide areas of the surfaces of small mammalian ventricles is straightforward given the relatively small scale of their myocardial dimensions compared to the tissue space constant and electrical field. However, it has yet to be proven for larger mammalian hearts with tissue properties and ventricular dimensions closer to humans. Our goal was to address the feasibility and impact of wide-area electrical stimulation on the ventricular surfaces of large mammalian hearts at different stimulus strengths. This was accomplished by placing long line electrodes on the ventricular surfaces of pig hearts that span wide areas, and activating them individually. Stimulus efficacy was assessed and compared between surfaces, and tissue viability was evaluated. Activation time was dependent on stimulation strength and location, achieving uniform linear stimulation at 9x threshold strength. Endocardial stimulation activated more tissue transmurally than epicardial stimulation, which could be considered a potential target for future cardiac electrotherapies. Overall, our results indicate that electrically stimulating wide areas of the ventricular surfaces of large mammals is achievable with line electrodes, minimal tissue damage, and energies under the human pain threshold (100 mJ).

Tissue Anisotropy Modeling Using Soft Composite Materials

Soft tissues in general exhibit anisotropic mechanical behavior, which varies in three dimensions based on the location of the tissue in the body. In the past, there have been few attempts to numerically model tissue anisotropy using composite-based formulations (involving fibers embedded within a matrix material). However, so far, tissue anisotropy has not been modeled experimentally. In the current work, novel elastomer-based soft composite materials were developed in the form of experimental test coupons, to model the macroscopic anisotropy in tissue mechanical properties. A soft elastomer matrix was fabricated, and fibers made of a stiffer elastomer material were embedded within the matrix material to generate the test coupons. The coupons were tested on a mechanical testing machine, and the resulting stress-versus-stretch responses were studied. The fiber volume fraction (FVF), fiber spacing, and orientations were varied to estimate the changes in the mechanical responses. The mechanical behavior of the soft composites was characterized using hyperelastic material models such as Mooney-Rivlin's, Humphrey's, and Veronda-Westmann's model and also compared with the anisotropic mechanical behavior of the human skin, pelvic tissues, and brain tissues. This work lays the foundation for the experimental modelling of tissue anisotropy, which combined with microscopic studies on tissues can lead to refinements in the simulation of localized fiber distribution and orientations, and enable the development of biofidelic anisotropic tissue phantom materials for various tissue engineering and testing applications.

EFFECTS OF ACUTE GLOBAL ISCHEMIA ON RE-ENTRANT ARRHYTHMOGENESIS: A SIMULATION STUDY

Sudden cardiac death is mainly caused by arrhythmogenesis. For a functional abnormal heart, such as an ischemic heart, the probability of arrhythmia occurring is greatly increased. During myocardial ischemia, re-entry is prone to degenerate into ventricular fibrillation (VF). Therefore it has important meaning to investigate the intricate mechanisms underlying VF under an ischemic condition in order to better facilitate therapeutic interventions. In this paper, to analyze the functional influence of acute global ischemia on cardiac electrical activity and subsequently on re-entrant arrhythmogenesis, we take into account three main pathophysiological consequences of ischemia: hyperkalaemia, acidosis, and anoxia, and develop a 3D human ventricular ischemic model that combines a detailed biophysical description of the excitation kinetics of human ventricular cells with an integrated geometry of human ventricular tissue which incorporates fiber direction anisotropy and the stimulation activation sequence. The results show that under acute global ischemia, the tissue excitability and the slope of ventricular cellular action potential duration restitution (APDR) are greatly decreased. As a result, the complexity of VF activation patterns is reduced. For the three components of ischemia, hyperkalaemia is the dominant contributor to the stability of re-entry under acute global ischemia. Increasing [K+]o acts to prolong the cell refractory period, reduce the tissue excitability and slow the conduction velocity. Our results also show that VF can be eliminated by decreasing cellular excitability, primarily by elevating the concentration value of extracellular K+.

Comparison of the relation between left ventricular anatomy and QRS duration in patients with cardiomyopathy with versus without left bundle branch block

QRS duration (QRSd) is used to diagnose left bundle branch block (LBBB) and is important to determine cardiac resynchronization therapy eligibility. The same QRSd thresholds established decades ago are used for all patients. However, significant interpatient variability of normal QRSd exists, and individualized QRSd thresholds might improve diagnosis and intervention strategies. Previous work reported left ventricular (LV) mass and papillary muscle location predicted QRSd in healthy subjects, but the relation in diseased ventricles is unknown. The aim of the present study was to determine the association between LV anatomy and QRSd in patients with cardiomyopathy. Patients referred for primary prevention implantable defibrillators (n = 166) received cardiac magnetic resonance imaging, and those with normal conduction (without bundle branch or fascicular block) and LBBB were studied. The LV mass, length, internal diameter, LV end-diastolic volume, septal and lateral wall thicknesses, and papillary muscle location were measured. In patients with normal conduction, LV length (r = 0.35, p <0.001), mass (r = 0.32, p <0.001), diameter (r = 0.20, p = 0.03), and septal wall thickness (r = 0.20, p = 0.03) had positive correlations with QRSd. In patients with LBBB, LV length (r = 0.32, p = 0.03), mass (r = 0.39, p = 0.01), diameter (r = 0.34, p = 0.02), and LV end-diastolic volume (r = 0.32, p = 0.04) had positive correlations with QRSd. Contrary to previous studies in healthy subjects, papillary muscle angle (location) was not associated with QRSd in cardiomyopathy patients with normal conduction or LBBB. In conclusion, increasing LV anatomical measurements were associated with increasing QRSd in patients with cardiomyopathy. Future work should investigate the use of LV anatomical measurements in developing individualized QRSd thresholds for diagnosing conduction abnormalities such as LBBB and identifying candidates for cardiac resynchronization therapy.

Effects of pacing site and stimulation history on alternans dynamics and the development of complex spatiotemporal patterns in cardiac tissue

Alternans of action potential duration has been associated with T wave alternans and the development of arrhythmias because it produces large gradients of repolarization. However, little is known about alternans dynamics in large mammalian hearts. Using optical mapping to record electrical activations simultaneously from the epicardium and endocardium of 9 canine right ventricles, we demonstrate novel arrhythmogenic complex spatiotemporal dynamics. (i) Alternans predominantly develops first on the endocardium. (ii) The postulated simple progression from normal rhythm to concordant to discordant alternans is not always observed; concordant alternans can develop from discordant alternans as the pacing period is decreased. (iii) In contrast to smaller tissue preparations, multiple stationary nodal lines may exist and need not be perpendicular to the pacing site or to each other. (iv) Alternans has fully three-dimensional dynamics and the epicardium and endocardium can show significantly different dynamics: multiple nodal surfaces can be transmural or intramural and can form concave/convex surfaces resulting in islands of discordant alternans. (v) The complex spatiotemporal patterns observed during alternans are very sensitive to both the site of stimulation and the stimulation history. Alternans in canine ventricles not only exhibit larger amplitudes and persist for longer cycle length regimes compared to those found in smaller mammalian hearts, but also show novel dynamics not previously described that enhance dispersion and show high sensitivity to initial conditions. This indicates some underlying predisposition to chaos and can help to guide the design of new drugs and devices controlling and preventing arrhythmic events.

Concurrent optimization of timing delays and electrode positioning in biventricular pacing based on a computer heart model assuming 17 left ventricular segments

BACKGROUND

The efficacy of cardiac resynchronization therapy through biventricular pacing (BVP) has been demonstrated by numerous studies in patients suffering from congestive heart failure. In order to achieve a guideline for optimal treatment with BVP devices, an automated non-invasive strategy based on a computer model of the heart is presented.

MATERIALS AND METHODS

The presented research investigates an off-line optimization algorithm regarding electrode positioning and timing delays. The efficacy of the algorithm is demonstrated in four patients suffering from left bundle branch block (LBBB) and myocardial infarction (MI). The computer model of the heart was used to simulate the LBBB in addition to several MI allocations according to the different left ventricular subdivisions introduced by the American Heart Association. Furthermore, simulations with reduced interventricular conduction velocity were performed in order to model interventricular excitation conduction delay. More than 800,000 simulations were carried out by adjusting a variety of 121 pairs of atrioventricular and interventricular delays and 36 different electrode positioning set-ups. Additionally, three different conduction velocities were examined. The optimization measures included the minimum root mean square error (E(RMS)) between physiological, pathological and therapeutic excitation, and also the difference of QRS-complex duration. Both of these measures were computed automatically.

RESULTS

Depending on the patient's pathology and conduction velocity, a reduction of E(RMS) between physiological and therapeutic excitation could be reached. For each patient and pathology, an optimal pacing electrode pair was determined. The results demonstrated the importance of an individual adjustment of BVP parameters to the patient's anatomy and pathology.

CONCLUSION

This work proposes a novel non-invasive optimization algorithm to find the best electrode positioning sites and timing delays for BVP in patients with LBBB and MI. This algorithm can be used to plan an optimal therapy for an individual patient.

Simulation of the QRS complex using papillary muscle positions as the site of early activation in human subjects

BACKGROUND

Simulation of the electrical activation of the heart and its comparison with real in vivo activation is a promising method in testing potential determinants of excitation. Simulation of the electrical activity of the human heart is now emerging as a step forward for understanding and predicting electrophysiologic patterns in humans. Initial points of excitation and the manner in which the activation spreads from these points are important variables determining QRS complex characteristics. It has been suggested that in humans, the initial excitation of the left ventricle is a primary determinant of QRS complex characteristics, and that excitation begins at the papillary muscles and septum, where the fascicles of the left bundle branch insert. The aim of this study is to test the hypothesis that QRS duration and direction of QRS axis in the frontal plane have excellent agreement between real QRS and simulated QRS using papillary muscle position to indicate the border of the origin of early ventricular activation.

METHODS

Fourteen healthy adult volunteers were included in the study. Magnetic resonance imaging data were obtained to assess the papillary muscle positions. Twelve-lead electrocardiographic (ECG) recordings were used to obtain real ECG data for assessment of QRS duration and QRS axis in each subject. Simulation software developed by ECG-TECH Corp (Huntington, NY) was used to simulate the ECG of each subject to determine simulated QRS duration and QRS frontal plane axis. QRS duration and QRS axis data were compared between simulated and real ECG and agreement between these variables was calculated.

RESULTS

Seventy-nine percent of subjects had a difference of the QRS duration between real and simulated ECG of less than 10 milliseconds. The calculated strength of agreement between simulated and real QRS duration was 71% and considered as "good" (kappa statistics). In 70% of subjects, the difference in the QRS axis was less than 10 degrees . The calculated strength of agreement between simulated and real QRS axis was 80% and considered as "excellent" (kappa statistics).

CONCLUSIONS

The results of this study suggest that the sites of the initiation of electrical activity in the left ventricle, as assessed by the positions of papillary muscles, may be considered as primary determinants of the QRS duration and QRS axis in humans. This knowledge may help in predicting normal QRS characteristic on a patient-specific basis. In this study, simulation of the QRS complex was based on papillary muscles from human hearts.

Innovation focus: the patient with arrhythmia

Great strides have been made over the last two decades in the management of patients with rhythm disorders. Despite this, however, the remaining critical problems of stroke related to atrial fibrillation or as a result of radiofrequency ablation require innovative solutions to fully realize the potential of these recent advances. Similarly, implanted cardiac devices have revolutionized the care of patients with bradyrhythmias and tachyarrhythmias. Dyssynchronus ventricular pacing associated with present devices; however, results in heart failure, tricuspid regurgitation, and inappropriate device therapy once again create a demand for creative solutions. While not technically an arrhythmia, epilepsy management today is riddled with many of the problems that plagued cardiac arrhythmia management previously, and thus an appreciation of the similarities in requirement for investigative solutions may yield groundbreaking solutions. In this paper, we describe some novel methods to reduce complications associated with rhythm disorders and their treatment and apply the lessons learned from cardiovascular arrhythmia management to the brain. These include: a method to reduce coagulum formation and thus subsequent thromboembolism with indwelling catheters specifically during radiofrequency ablation procedures; a technique to ligate the left atrial appendage through percutaneous subxiphoid pericardial access; development and testing of a novel intramyocardial pace-sense lead, particularly used in a unique anatomic location (the atrioventricular septum) to allow pacing the ventricles in a relatively synchronous manner without crossing the tricuspid valve or entering the coronary sinus; finally, novel modifications of the cardiovascular mapping and ablation techniques used for the management of the central nervous system disorders primarily via the venous drainage of the brain. Innovative and potential solutions to treat the patient with arrhythmia are presented.