Critical Volume of Human Myocardium Necessary to Maintain Ventricular Fibrillation

Human sodium current voltage-dependence at physiological temperature measured by coupling a patch-clamp experiment to a mathematical model

Voltage-gated Na+ channels are crucial to action potential propagation in excitable tissues. Because of the high amplitude and rapid activation of the Na+ current, voltage-clamp measurements are very challenging and are usually performed at room temperature. In this study, we measured Na+ current voltage-dependence in stem cell-derived cardiomyocytes at physiological temperature. While the apparent activation and inactivation curves, measured as the dependence of current amplitude on voltage, fall within the range reported in previous studies, we identified a systematic error in our measurements. This error is caused by the deviation of the membrane potential from the command potential of the amplifier. We demonstrate that it is possible to account for this artifact using computer simulation of the patch-clamp experiment. We obtained surprising results through patch-clamp model optimization: a half-activation of -11.5 mV and a half-inactivation of -87 mV. Although the half-activation deviates from previous research, we demonstrate that this estimate reproduces the conduction velocity dependence on extracellular potassium concentration. KEY POINTS: Voltage-gated Na+ currents play a crucial role in excitable tissues including neurons, cardiac and skeletal muscle. Measurement of Na+ current is challenging because of its high amplitude and rapid kinetics, especially at physiological temperature. We have used the patch-clamp technique to measure human Na+ current voltage-dependence in human induced pluripotent stem cell-derived cardiomyocytes. The patch-clamp data were processed by optimization of the model accounting for voltage-clamp experiment artifacts, revealing a large difference between apparent parameters of Na+ current and the results of the optimization. We conclude that actual Na+ current activation is extremely depolarized in comparison to previous studies. The new Na+ current model provides a better understanding of action potential propagation; we demonstrate that it explains propagation in hyperkalaemic conditions.

Age-related structural and functional changes of the intracardiac nervous system

BACKGROUND

Although aging is known to be associated with an increased incidence of both atrial and ventricular arrhythmias, there is limited knowledge about how Schwann cells (SC) and the intracardiac nervous system (iCNS) remodel with age. Here we investigate the differences in cardiac SC, parasympathetic nerve fibers, and muscarinic acetylcholine receptor M2 (M2R) expression in young and old mice. Additionally, we examine age-related changes in cardiac responses to sympathomimetic and parasympathomimetic drugs.

METHODS AND RESULTS

Lower SC density, lower SC proliferation and fewer parasympathetic nerve fibers were observed in cardiac and, as a control sciatic nerves from old (20-24 months) compared to young mice (2-3 months). In old mice, chondroitin sulfate proteoglycan 4 (CSPG4) was increased in sciatic but not cardiac nerves. Expression of M2R was lower in ventricular myocardium and ventricular conduction system from old mice compared to young mice, while no significant difference was seen in M2R expression in sino-atrial or atrio-ventricular node pacemaker tissue. Heart rate was slower and PQ intervals were longer in Langendorff-perfused hearts from old mice. Ventricular tachycardia and fibrillation were more frequently observed in response to carbachol administration in hearts from old mice versus those from young mice.

CONCLUSIONS

On the background of reduced presence of SC and parasympathetic nerve fibers, and of lower M2R expression in ventricular cardiomyocytes and conduction system of aged hearts, the propensity of ventricular arrhythmogenesis upon parasympathomimetic drug application is increased. Whether this is caused by an increase in heterogeneity of iCNS structure and function remains to be elucidated.

Numerical study of the drift of scroll waves by optical feedback in cardiac tissue

Nonlinear waves were found in various types of physical, chemical, and biological excitable media, e.g., in heart muscle. They can form three-dimensional (3D) vortices, called scroll waves, that are of particular significance in the heart, as they underlie lethal cardiac arrhythmias. Thus controlling the behavior of scroll waves is interesting and important. Recently, the optical feedback control procedure for two-dimensional vortices, called spiral waves, was developed. It can induce directed linear drift of spiral waves in optogenetically modified cardiac tissue. However, the extension of this methodology to 3D scroll waves is nontrivial, as optogenetic signals only penetrate close to the surface of cardiac tissue. Here we present a study of this extension in a two-variable reaction-diffusion model and in a detailed model of cardiac tissue. We show that the success of the control procedure is determined by the tension of the scroll wave filament. In tissue with positive filament tension the control procedure works in all cases. However, in the case of negative filament tension for a sufficiently large medium, instabilities occur and make drift and control of scroll waves impossible. Because in normal cardiac tissue the filament tension is assumed to be positive, we conclude that the proposed optical feedback scheme can be a robust method in inducing the linear drift of scroll waves that can control their positions in cardiac tissue.

Hardware-Mappable Cellular Neural Networks for Distributed Wavefront Detection in Next-Generation Cardiac Implants

Artificial intelligence algorithms are being adopted to analyze medical data, promising faster interpretation to support doctors' diagnostics. The next frontier is to bring these powerful algorithms to implantable medical devices. Herein, a closed-loop solution is proposed, where a cellular neural network is used to detect abnormal wavefronts and wavebrakes in cardiac signals recorded in human tissue is trained to achieve >96% accuracy, >92% precision, >99% specificity, and >93% sensitivity, when floating point precision weights are assumed. Unfortunately, the current hardware technologies for floating point precision are too bulky or energy intensive for compact standalone applications in medical implants. Emerging device technologies, such as memristors, can provide the compact and energy-efficient hardware fabric to support these efforts and can be reliably embedded with existing sensor and actuator platforms in implantable devices. A distributed design that considers the hardware limitations in terms of overhead and limited bit precision is also discussed. The proposed distributed solution can be easily adapted to other medical technologies that require compact and efficient computing, like wearable devices and lab-on-chip platforms.

Electrophysiology and Arrhythmogenesis in the Human Right Ventricular Outflow Tract

BACKGROUND

Right ventricular outflow tract (RVOT) is a common source of ventricular tachycardia, which often requires ablation. However, the mechanisms underlying the RVOT's unique arrhythmia susceptibility remain poorly understood due to lack of detailed electrophysiological and molecular studies of the human RVOT.

METHODS

We conducted optical mapping studies in 16 nondiseased donor human RVOT preparations subjected to pharmacologically induced adrenergic and cholinergic stimulation to evaluate susceptibility to arrhythmias and characterize arrhythmia dynamics.

RESULTS

We found that under control conditions, RVOT has shorter action potential duration at 80% repolarization relative to the right ventricular apical region. Treatment with isoproterenol (100 nM) shortened action potential duration at 80% repolarization and increased incidence of premature ventricular contractions (P=0.003), whereas acetylcholine (100 μM) stimulation alone had no effect on action potential duration at 80% repolarization or premature ventricular contractions. However, acetylcholine treatment after isoproterenol stimulation reduced the incidence of premature ventricular contractions (P=0.034) and partially reversed action potential duration at 80% repolarization shortening (P=0.029). Immunolabeling of RVOT (n=4) confirmed the presence of cholinergic marker VAChT (vesicular acetylcholine transporter) in the region. Rapid pacing revealed RVOT susceptibility to both concordant and discordant alternans. Investigation into transmural arrhythmia dynamics showed that arrhythmia wave fronts and phase singularities (rotors) were relatively more organized in the endocardium than in the epicardium (P=0.006). Moreover, there was a weak but positive spatiotemporal autocorrelation between epicardial and endocardial arrhythmic wave fronts and rotors. Transcriptome analysis (n=10 hearts) suggests a trend that MAPK (mitogen-activated protein kinase) signaling, calcium signaling, and cGMP-PKG (protein kinase G) signaling are among the pathways that may be enriched in the male RVOT, whereas pathways of neurodegeneration may be enriched in the female RVOT.

CONCLUSIONS

Human RVOT electrophysiology is characterized by shorter action potential duration relative to the right ventricular apical region. Cholinergic right ventricular stimulation attenuates the arrhythmogenic effects of adrenergic stimulation, including increase in frequency of premature ventricular contractions and shortening of wavelength. Right ventricular arrhythmia is characterized by positive spatial-temporal autocorrelation between epicardial-endocardial arrhythmic wave fronts and rotors that are relatively more organized in the endocardium.

Cardiac Conduction Velocity, Remodeling and Arrhythmogenesis

Cardiac electrophysiological disorders, in particular arrhythmias, are a key cause of morbidity and mortality throughout the world. There are two basic requirements for arrhythmogenesis: an underlying substrate and a trigger. Altered conduction velocity (CV) provides a key substrate for arrhythmogenesis, with slowed CV increasing the probability of re-entrant arrhythmias by reducing the length scale over which re-entry can occur. In this review, we examine methods to measure cardiac CV in vivo and ex vivo, discuss underlying determinants of CV, and address how pathological variations alter CV, potentially increasing arrhythmogenic risk. Finally, we will highlight future directions both for methodologies to measure CV and for possible treatments to restore normal CV.

Theoretical Models and Computational Analysis of Action Potential Dispersion for Cardiac Arrhythmia Risk Stratification

Reentrant cardiac arrhythmias such as atrial fibrillation (AF) and ventricular fibrillation (VF) are common cardiac arrhythmias that account for substantial morbidity and mortality throughout the world. However, the mechanisms and optimal ablation treatment strategies for such arrhythmias are still unclear. Using 2D optical mapping of a mouse model with AF and VF, we have identified regional heterogeneity of the action potential duration (APD) in the atria and ventricles of the heart as key drivers for the initiation and persistence of reentry. The purpose of this paper is to discuss theoretical patterns of dispersion, demonstrate patterns of dispersion seen in our mouse model and discuss the computational analysis of APD dispersion patterns. These analyses and discussions may lead to better understanding of dispersion patterns in patients with these arrhythmias, as well as help comprehend whether and how reducing dispersion can lead to arrhythmia risk stratification and treatment strategies for arrhythmias.

Prolonged ursodeoxycholic acid administration reduces acute ischaemia-induced arrhythmias in adult rat hearts

Acute myocardial ischaemia and reperfusion (I-R) are major causes of ventricular arrhythmias in patients with a history of coronary artery disease. Ursodeoxycholic acid (UDCA) has previously been shown to be antiarrhythmic in fetal hearts. This study was performed to investigate if UDCA protects against ischaemia-induced and reperfusion-induced arrhythmias in the adult myocardium, and compares the effect of acute (perfusion only) versus prolonged (2 weeks pre-treatment plus perfusion) UDCA administration. Langendorff-perfused adult Sprague-Dawley rat hearts were subjected to acute regional ischaemia by ligation of the left anterior descending artery (10 min), followed by reperfusion (2 min), and arrhythmia incidence quantified. Prolonged UDCA administration reduced the incidence of acute ischaemia-induced arrhythmias (p = 0.028), with a reduction in number of ventricular ectopic beats during the ischaemic phase compared with acute treatment (10 ± 3 vs 58 ± 15, p = 0.036). No antiarrhythmic effect was observed in the acute UDCA administration group. Neither acute nor prolonged UDCA treatment altered the incidence of reperfusion arrhythmias. The antiarrhythmic effect of UDCA may be partially mediated by an increase in cardiac wavelength, due to the attenuation of conduction velocity slowing (p = 0.03), and the preservation of Connexin43 phosphorylation during acute ischaemia (p = 0.0027). The potential antiarrhythmic effects of prolonged UDCA administration merit further investigation.

Catheter-integrated soft multilayer electronic arrays for multiplexed sensing and actuation during cardiac surgery

The rigidity and relatively primitive modes of operation of catheters equipped with sensing or actuation elements impede their conformal contact with soft-tissue surfaces, limit the scope of their uses, lengthen surgical times and increase the need for advanced surgical skills. Here, we report materials, device designs and fabrication approaches for integrating advanced electronic functionality with catheters for minimally invasive forms of cardiac surgery. By using multiphysics modelling, plastic heart models and Langendorff animal and human hearts, we show that soft electronic arrays in multilayer configurations on endocardial balloon catheters can establish conformal contact with curved tissue surfaces, support high-density spatiotemporal mapping of temperature, pressure and electrophysiological parameters and allow for programmable electrical stimulation, radiofrequency ablation and irreversible electroporation. Integrating multimodal and multiplexing capabilities into minimally invasive surgical instruments may improve surgical performance and patient outcomes.

Genetic algorithm-based personalized models of human cardiac action potential

We present a novel modification of genetic algorithm (GA) which determines personalized parameters of cardiomyocyte electrophysiology model based on set of experimental human action potential (AP) recorded at different heart rates. In order to find the steady state solution, the optimized algorithm performs simultaneous search in the parametric and slow variables spaces. We demonstrate that several GA modifications are required for effective convergence. Firstly, we used Cauchy mutation along a random direction in the parametric space. Secondly, relatively large number of elite organisms (6-10% of the population passed on to new generation) was required for effective convergence. Test runs with synthetic AP as input data indicate that algorithm error is low for high amplitude ionic currents (1.6±1.6% for IKr, 3.2±3.5% for IK1, 3.9±3.5% for INa, 8.2±6.3% for ICaL). Experimental signal-to-noise ratio above 28 dB was required for high quality GA performance. GA was validated against optical mapping recordings of human ventricular AP and mRNA expression profile of donor hearts. In particular, GA output parameters were rescaled proportionally to mRNA levels ratio between patients. We have demonstrated that mRNA-based models predict the AP waveform dependence on heart rate with high precision. The latter also provides a novel technique of model personalization that makes it possible to map gene expression profile to cardiac function.

Granger Causality-Based Analysis for Classification of Fibrillation Mechanisms and Localization of Rotational Drivers

BACKGROUND

The mechanisms sustaining myocardial fibrillation remain disputed, partly due to a lack of mapping tools that can accurately identify the mechanism with low spatial resolution clinical recordings. Granger causality (GC) analysis, an econometric tool for quantifying causal relationships between complex time-series, was developed as a novel fibrillation mapping tool and adapted to low spatial resolution sequentially acquired data.

METHODS

Ventricular fibrillation (VF) optical mapping was performed in Langendorff-perfused Sprague-Dawley rat hearts (n=18), where novel algorithms were developed using GC-based analysis to (1) quantify causal dependence of neighboring signals and plot GC vectors, (2) quantify global organization with the causality pairing index, a measure of neighboring causal signal pairs, and (3) localize rotational drivers (RDs) by quantifying the circular interdependence of neighboring signals with the circular interdependence value. GC-based mapping tools were optimized for low spatial resolution from downsampled optical mapping data, validated against high-resolution phase analysis and further tested in previous VF optical mapping recordings of coronary perfused donor heart left ventricular wedge preparations (n=12), and adapted for sequentially acquired intracardiac electrograms during human persistent atrial fibrillation mapping (n=16).

RESULTS

Global VF organization quantified by causality pairing index showed a negative correlation at progressively lower resolutions (50% resolution: P=0.006, R2=0.38, 12.5% resolution, P=0.004, R2=0.41) with a phase analysis derived measure of disorganization, locations occupied by phase singularities. In organized VF with high causality pairing index values, GC vector mapping characterized dominant propagating patterns and localized stable RDs, with the circular interdependence value showing a significant difference in driver versus nondriver regions (0.91±0.05 versus 0.35±0.06, P=0.0002). These findings were further confirmed in human VF. In persistent atrial fibrillation, a positive correlation was found between the causality pairing index and presence of stable RDs (P=0.0005,R2=0.56). Fifty percent of patients had RDs, with a low incidence of 0.9±0.3 RDs per patient.

CONCLUSIONS

GC-based fibrillation analysis can measure global fibrillation organization, characterize dominant propagating patterns, and map RDs using low spatial resolution sequentially acquired data.

Optocardiography and Electrophysiology Studies of Ex Vivo Langendorff-perfused Hearts

Small animal models are most commonly used in cardiovascular research due to the availability of genetically modified species and lower cost compared to larger animals. Yet, larger mammals are better suited for translational research questions related to normal cardiac physiology, pathophysiology, and preclinical testing of therapeutic agents. To overcome the technical barriers associated with employing a larger animal model in cardiac research, we describe an approach to measure physiological parameters in an isolated, Langendorff-perfused piglet heart. This approach combines two powerful experimental tools to evaluate the state of the heart: electrophysiology (EP) study and simultaneous optical mapping of transmembrane voltage and intracellular calcium using parameter sensitive dyes (RH237, Rhod2-AM). The described methodologies are well suited for translational studies investigating the cardiac conduction system, alterations in action potential morphology, calcium handling, excitation-contraction coupling and the incidence of cardiac alternans or arrhythmias.

Optical Mapping

Optical mapping of electrical activity in the heart is based on voltage-sensitive and lipophilic fluorescence dyes. Optical signals recorded from cardiac cells correlate well with their transmembrane potentials. High spatiotemporal resolution, wide field mapping, and high sensitivity to transmembrane potential enable detailed characterization of action potential initiation and propagation. Optical mapping is used to study complex patterns of excitation propagation, including propagation across the sinoatrial and atrioventricular nodes and during atrial and ventricular arrhythmias.Optical mapping is used to study the role of reentrant activity in atrial and ventricular fibrillation.

Spontaneous termination of chaotic spiral wave dynamics in human cardiac ion channel models

Chaotic spiral or scroll wave dynamics can be found in diverse systems. In cardiac dynamics, spiral or scroll waves of electrical excitation determine the dynamics during life-threatening arrhythmias like ventricular fibrillation. In numerical studies it was found that chaotic episodes of spiral and scroll waves can be transient, thus they terminate spontaneously. We show in this study that this behavior can also be observed using models which describe the ion channel dynamics of human cardiomyocytes (Bueno-Orovio-Cherry-Fenton model and the Ten Tusscher-Noble-Noble-Panfilov model). For both models we find that the average lifetime of the chaotic transients grows exponentially with the system size. With this behavior, we classify the systems into the group of type-II supertransients. We observe a significant difference of the breakup behavior between the models, which results in a distinct dynamics during the final phase just before the termination. The observation of a (temporally) stable single-spiral state affects the prevailing description of the dynamics of type-II supertransients as being “quasi-stationary” and also the feasibility of predicting the spontaneous termination of the spiral wave dynamics. In the long term, the relation between the breakup behavior of spiral waves and properties of chaotic transients like predictability or average transient lifetime may contribute to an improved understanding and classification of cardiac arrhythmias.