Population of computational rabbit-specific ventricular action potential models for investigating sources of variability in cellular repolarisation

Variability is observed at all levels of cardiac electrophysiology. Yet, the underlying causes and importance of this variability are generally unknown, and difficult to investigate with current experimental techniques. The aim of the present study was to generate populations of computational ventricular action potential models that reproduce experimentally observed intercellular variability of repolarisation (represented by action potential duration) and to identify its potential causes. A systematic exploration of the effects of simultaneously varying the magnitude of six transmembrane current conductances (transient outward, rapid and slow delayed rectifier K(+), inward rectifying K(+), L-type Ca(2+), and Na(+)/K(+) pump currents) in two rabbit-specific ventricular action potential models (Shannon et al. and Mahajan et al.) at multiple cycle lengths (400, 600, 1,000 ms) was performed. This was accomplished with distributed computing software specialised for multi-dimensional parameter sweeps and grid execution. An initial population of 15,625 parameter sets was generated for both models at each cycle length. Action potential durations of these populations were compared to experimentally derived ranges for rabbit ventricular myocytes. 1,352 parameter sets for the Shannon model and 779 parameter sets for the Mahajan model yielded action potential duration within the experimental range, demonstrating that a wide array of ionic conductance values can be used to simulate a physiological rabbit ventricular action potential. Furthermore, by using clutter-based dimension reordering, a technique that allows visualisation of multi-dimensional spaces in two dimensions, the interaction of current conductances and their relative importance to the ventricular action potential at different cycle lengths were revealed. Overall, this work represents an important step towards a better understanding of the role that variability in current conductances may play in experimentally observed intercellular variability of rabbit ventricular action potential repolarisation.

A compliant, banded outflow cannula for decreased afterload sensitivity of rotary right ventricular assist devices

Biventricular support with dual rotary ventricular assist devices (VADs) has been implemented clinically with restriction of the right VAD (RVAD) outflow cannula to artificially increase afterload and, therefore, operate within recommended design speed ranges. However, the low preload and high afterload sensitivity of these devices increase the susceptibility of suction events. Active control systems are prone to sensor drift or inaccurate inferred (sensor-less) data, therefore an alternative solution may be of benefit. This study presents the in vitro evaluation of a compliant outflow cannula designed to passively decrease the afterload sensitivity of rotary RVADs and minimize left-sided suction events. A one-way fluid-structure interaction model was initially used to produce a design with suitable flow dynamics and radial deformation. The resultant geometry was cast with different initial cross-sectional restrictions and concentrations of a softening diluent before evaluation in a mock circulation loop. Pulmonary vascular resistance (PVR) was increased from 50 dyne s/cm5 until left-sided suction events occurred with each compliant cannula and a rigid, 4.5 mm diameter outflow cannula for comparison. Early suction events (PVR ∼ 300 dyne s/cm5) were observed with the rigid outflow cannula. Addition of the compliant section with an initial 3 mm diameter restriction and 10% diluent expanded the outflow restriction as PVR increased, thus increasing RVAD flow rate and preventing left-sided suction events at PVR levels beyond 1000 dyne s/cm5. Therefore, the compliant, restricted outflow cannula provided a passive control system to assist in the prevention of suction events with rotary biventricular support while maintaining pump speeds within normal ranges of operation.

Mechanism design of ventricular assist device

This research presents an innovative design approach for the development of high efficiency Ventricular assist device that can be used for long-term support a heart failure patient. Computational fluid dynamics (CFD) techniques were applied to the development and intensive analysis to improve the performance and reliability of the pump. From the CFD analysis, a prototype pump was created and evaluated on the mock circulation loop that simulate the human circulatory system environment to evaluate its performance in support varying heart conditions.

Ventricular repolarization during cardiovascular autonomic function testing

The autonomic nervous system is an important modulator of ventricular repolarization and arrhythmia vulnerability. This study explored the effects of cardiovascular autonomic function tests on repolarization and its heterogeneity, with a special reference to congenital arrhythmogenic disorders typically associated with stress-induced fatal ventricular arrhythmias. The first part explored the effects of standardized autonomic tests on QT intervals in a 12-lead electrocardiogram and in multichannel magnetocardiography in 10 healthy adults. The second part studied the effects of deep breathing, Valsalva manouvre, mental stress, sustained handgrip and mild exercise on QT intervals in asymptomatic patients with LQT1 subtype of the hereditary long QT syndrome (n=9) and in patients with arrhythmogenic right ventricular dysplasia (ARVD, n=9). Even strong sympathetic activation had no effects on spatial QT interval dispersion in healthy subjects, but deep respiratory efforts and Valsalva influenced it in ways that were opposite in electrocardiographic and magnetocardiographic recordings. LQT1 patients showed blunted QT interval and sinus nodal responses to sympathetic challenge, as well as an exaggerated QT prolongation during the recovery phases. LQT1 patients showed a QT interval recovery overshoot in 2.4 ± 1.7 tests compared with 0.8 ± 0.7 in healthy controls (P = 0.02). Valsalva strain prolonged the T wave peak to T wave end interval only in the LQT1 patients, considered to reflect the arrhythmogenic substrate in this syndrome. ARVD patients showed signs of abnormal repolarization in the right ventricle, modulated by abrupt sympathetic activation. An electrocardiographic marker reflecting interventricular dispersion of repolarization was introduced. It showed that LQT1 patients exhibit a repolarization gradient from the left ventricle towards the right ventricle, significantly larger than in controls. In contrast, ARVD patients showed a repolarization gradient from the right ventricle towards the left. Valsalva strain amplified the repolarization gradient in LQT1 patients whereas it transiently reversed it in patients with ARVD. In conclusion, intrathoracic volume and pressure changes influence regional electrocardiographic and magnetocardiographic QT interval measurements differently. Especially recovery phases of standard cardiovascular autonomic functions tests and Valsalva manoeuvre reveal the abnormal repolarization in asymptomatic LQT1 patients. Both LQT1 and ARVD patients have abnormal interventricular repolarization gradients, modulated by abrupt sympathetic activation. Autonomic testing and in particular the Valsalva manoeuvre are potentially useful in unmasking abnormal repolarization in these syndromes.

Ion channel gene variants predisposing to severe human ventricular arrhythmias

Congenital long QT syndrome (LQTS) with an estimated prevalence of 1:2000-1:10 000 manifests with prolonged QT interval on electrocardiogram and risk for ventricular arrhythmias and sudden death. Several ion channel genes and hundreds of mutations in these genes have been identified to underlie the disorder. In Finland, four LQTS founder mutations of potassium channel genes account for up to 40-70% of genetic spectrum of LQTS. Acquired LQTS has similar clinical manifestations, but often arises from usage of QT-prolonging medication or electrolyte disturbances. A prolonged QT interval is associated with increased morbidity and mortality not only in clinical LQTS but also in patients with ischemic heart disease and in the general population. The principal aim of this study was to estimate the actual prevalence of LQTS founder mutations in Finland and to calculate their effect on QT interval in the Finnish background population. Using a large population-based sample of over 6000 Finnish individuals from the Health 2000 Survey, we identified LQTS founder mutations KCNQ1 G589D (n=8), KCNQ1 IVS7-2A>G (n=1), KCNH2 L552S (n=2), and KCNH2 R176W (n=16) in 27 study participants. This resulted in a weighted prevalence estimate of 0.4% for LQTS in Finland. Using a linear regression model, the founder mutations resulted in a 22- to 50-ms prolongation of the age-, sex-, and heart rate-adjusted QT interval. Collectively, these data suggest that one of 250 individuals in Finland may be genetically predisposed to ventricular arrhythmias arising from the four LQTS founder mutations. A KCNE1 D85N minor allele with a frequency of 1.4% was associated with a 10-ms prolongation in adjusted QT interval and could thus identify individuals at increased risk of ventricular arrhythmias at the population level. In addition, the previously reported associations of KCNH2 K897T, KCNH2 rs3807375, and NOS1AP rs2880058 with QT interval duration were confirmed in the present study. In a separate study, LQTS founder mutations were identified in a subgroup of acquired LQTS, providing further evidence that congenital LQTS gene mutations may underlie acquired LQTS. Catecholaminergic polymorphic ventricular tachycardia (CPVT) is characterized by exercise-induced ventricular arrhythmias in a structurally normal heart and results from defects in the cardiac Ca2+ signaling proteins, mainly ryanodine receptor type 2 (RyR2). In a patient population of typical CPVT, RyR2 mutations were identifiable in 25% (4/16) of patients, implying that noncoding variants or other genes are involved in CPVT pathogenesis. A 1.1 kb RyR2 exon 3 deletion was identified in two patients independently, suggesting that this region may provide a new target for RyR2-related molecular genetic studies. Two novel RyR2 mutations showing a gain-of-function defect in vitro were identified in three victims of sudden cardiac death. Extended pedigree analyses revealed some surviving mutation carriers with mild structural abnormalities of the heart and resting ventricular arrhythmias suggesting that not all RyR2 mutations lead to a typical CPVT phenotype, underscoring the relevance of tailored risk stratification of a RyR2 mutation carrier.

Electrocardiographic parameters of ventricular repolarization : modifiers and the prognostic value

Electric activity of the heart consists of repeated cardiomyocyte depolarizations and repolarizations. Abnormalities in repolarization predispose to ventricular arrhythmias. In body surface electrocardiogram, ventricular repolarization generates the T wave. Several electrocardiographic measures have been developed both for clinical and research purposes to detect repolarization abnormalities. The study aim was to investigate modifiers of ventricular repolarization with the focus on the relationship of the left ventricular mass, antihypertensive drugs, and common gene variants, to electrocardiographic repolarization parameters. The prognostic value of repolarization parameters was also assessed. The study subjects originated from a population of more than 200 middle-aged hypertensive men attending the GENRES hypertension study, and from an epidemiological survey, the Health 2000 Study, including more than 6000 participants. Ventricular repolarization was analysed from digital standard 12-lead resting electrocardiograms with two QT-interval based repolarization parameters (QT interval, T-wave peak to T-wave end interval) and with a set of four T-wave morphology parameters. The results showed that in hypertensive men, a linear change in repolarization parameters is present even in the normal range of left ventricular mass, and that even mild left ventricular hypertrophy is associated with potentially adverse electrocardiographic repolarization changes. In addition, treatments with losartan, bisoprolol, amlodipine, and hydrochlorothiazide have divergent short-term effects on repolarization parameters in hypertensive men. Analyses of the general population sample showed that single nucleotide polymorphisms in KCNH2, KCNE1, and NOS1AP genes are associated with changes in QT-interval based repolarization parameters but not consistently with T-wave morphology parameters. T-wave morphology parameters, but not QT interval or T-wave peak to T-wave end interval, provided independent prognostic information on mortality. The prognostic value was specifically related to cardiovascular mortality. The results indicate that, in hypertension, altered ventricular repolarization is already present in mild left ventricular mass increase, and that commonly used antihypertensive drugs may relatively rapidly and treatment-specifically modify electrocardiographic repolarization parameters. Common variants in cardiac ion channel genes and NOS1AP gene may also modify repolarization-related arrhythmia vulnerability. In the general population, T-wave morphology parameters may be useful in the risk assessment of cardiovascular mortality.

Diabetic cardiomyopathy and post-infarct ventricular remodelling : Effects of levosimendan in a rodent model of type II diabetes

Type 2 diabetes is a risk factor for the development of cardiovascular disease. Recently, the term diabetic cardiomyopathy has been proposed to describe the changes in the heart that occur in response to chronic hyperglycemia and insulin resistance. Ventricular remodelling in diabetic cardiomyopathy includes left ventricular hypertrophy, increased interstitial fibrosis, apoptosis and diastolic dysfunction. Mechanisms behind these changes are increased oxidative stress and renin-angiotensin system activation. The diabetic Goto-Kakizaki rat is a non-obese model of type 2 diabetes that exhibits defective insulin signalling. Recently two interconnected stress response pathways have been discovered that link insulin signalling, longevity, apoptosis and cardiomyocyte hypertrophy. The insulin-receptor PI3K/Ak pathway inhibits proapoptotic FOXO3a in response to insulin signalling and the nuclear Sirt1 deacetylase inhibits proapoptotic p53 and modulates FOXO3a in favour of survival and growth. --- Levosimendan is a calcium sensitizing agent used for the management of acute decompensated heart failure. Levosimendan acts as a positive inotrope by sensitizing cardiac troponin C to calcium and exerts vasodilation by opening mitochondrial and sarcolemmal ATP-sensitive potassium channels. Levosimendan has been described to have beneficial effects in ventricular remodelling after myocardial infarction. The aims of the study were to characterize whether diabetic cardiomyopathy associates with cardiac dysfunction, cardiomyocyte apoptosis, hypertrophy and fibrosis in spontaneously diabetic Goto-Kakizaki (GK) rats, which were used to model type 2 diabetes. Protein expression and activation of the Akt FOXO3a and Sirt1 p53 pathways were examined in the development of ventricular remodelling in GK rats with and without myocardial infarction (MI). The third and fourth studies examined the effects of levosimendan on ventricular remodelling and gene expression in post-MI GK rats. The results demonstrated that diabetic GK rats develop both modest hypertension and features similar to diabetic cardiomyopathy including cardiac dysfunction, LV hypertrophy and fibrosis and increased apoptotic signalling. MI induced a sustained increase in cardiomyocyte apoptosis in GK rats together with aggravated LV hypertrophy and fibrosis. The GK rat myocardium exhibited decreased Akt- FOXO3a phosphorylation and increased nuclear translocation of FOXO3a and overproduction of the Sirt1 protein. Treatment with levosimendan decreased cardiomyocyte apoptosis, senescence and LV hypertrophy and altered the gene expression profile in GK rat myocardium. The findings indicate that impaired cardioprotection via Akt FOXO3a and p38 MAPK is associated with increased apoptosis, whereas Sirt1 functions in counteracting apoptosis and the development of LV hypertrophy in the GK rat myocardium. Overall, levosimendan treatment protects against post-MI ventricular remodelling and alters the gene expression profile in the GK rat myocardium.

Modification of ventricular repolarization in silent long QT syndrome mutation carriers

Congenital long QT syndrome (LQTS) is a familial disorder characterized by ventricular repolarization that makes carriers vulnerable to malignant ventricular tachycardia and sudden cardiac death. The three main subtypes (LQT1, LQT2 and LQT3) constitute 95% of cases. The disorder is characterized by a prolonged QT interval in electrocardiograms (ECG), but a considerable portion are silent carriers presenting normal (QTc < 440 ms) or borderline (QTc < 470 ms) QT interval. Genetic testing is available only for 60-70% of patients. A number of pharmaceutical compounds also affect ventricular repolarization, causing a clinically similar disorder called acquired long QT syndrome. LQTS carriers - who already have impaired ventricular repolarization - are especially vulnerable. In this thesis, asymptomatic genotyped LQTS mutation carriers with non-diagnostic resting ECG were studied. The body surface potential mapping (BSPM) system was utilized for ECG recording, and signals were analyzed with an automated analysis program. QT interval length, and the end part of the T wave, the Tpe interval, was studied during exercise stress testing and an epinephrine bolus test. In the latter, T wave morphology was also analyzed. The effect of cetirizine was studied in LQTS carriers and also with supra- therapeutic dose in healthy volunteers. At rest, LQTS mutation carriers had a slightly longer heart rate adjusted QTc interval than healthy subjects (427 ± 31 ms and 379 ± 26 ms; p<0.001), but significant overlapping existed. LQT2 mutation carriers had a conspicuously long Tpe-interval (113 ± 24 ms; compared to 79 ± 11 ms in LQT1, 81 ± 17 ms in LQT3 and 78 ± 10 ms in controls; p<0.001). In exercise stress tests, LQT1 mutation carriers exhibit a long QT interval at high heart rates and during recovery, whereas LQT2 mutation carriers have a long Tpe interval at the beginning of exercise and at the end of recovery at low heart rates. LQT3 mutation carriers exhibit prominent shortening of both QT and Tpe intervals during exercise. A small epinephrine bolus revealed disturbed repolarization, especially in LQT2 mutation carriers, who developed prolonged Tpe intervals. A higher epinephrine bolus caused abnormal T waves with a different T wave profile in LQTS mutation carriers compared to healthy controls. These effects were seen in LQT3 as well, a group that may easily escape other provocative tests. In the cetirizine test, the QT and Tpe intervals were not prolonged in LQTS mutation carriers or in healthy controls. Subtype-specific findings in exercise test and epinephrine bolus test help to diagnose silent LQTS mutation carriers and to guide subtype-specific treatments. The Tpe interval, which signifies the repolarization process, seems to be a sensitive marker of disturbed repolarization along with the QT interval, which signifies the end of repolarization. This method may be used in studying compounds that are suspected to affect repolarization. Cetirizine did not adversely alter ventricular repolarization and would not be pro-arrhythmic in common LQT1 and LQT2 subtypes when used at its recommended doses.

The Mathematical modeling of inhomogeneites in ventricular tissue

Cardiac arrhythmias such as ventricular tachycardia (VT) or ventricular fibrillation (VF) are the leading cause of death in the industrialised world. There is a growing consensus that these arrhythmias arise because of the formation of spiral waves of electrical activation in cardiac tissue; unbroken spiral waves are associated with VT and broken ones with VF. Several experimental studies have been carried out to determine the effects of inhomogeneities in cardiac tissue on such arrhythmias. We give a brief overview of such experiments, and then an introduction to partial-differential-equation models for ventricular tissue. We show how different types of inhomogeneities can be included in such models, and then discuss various numerical studies, including our own, of the effects of these inhomogeneities on spiral-wave dynamics. The most remarkable qualitative conclusion of our studies is that the spiral-wave dynamics in such systems depends very sensitively on the positions of these inhomogeneities.

Spiral-Wave Dynamics in a Mathematical Model of Human Ventricular Tissue with Myocytes and Fibroblasts

Cardiac fibroblasts, when coupled functionally with myocytes, can modulate the electrophysiological properties of cardiac tissue. We present systematic numerical studies of such modulation of electrophysiological properties in mathematical models for (a) single myocyte-fibroblast (MF) units and (b) two-dimensional (2D) arrays of such units; our models build on earlier ones and allow for zero-, one-, and two-sided MF couplings. Our studies of MF units elucidate the dependence of the action-potential (AP) morphology on parameters such as E-f, the fibroblast resting-membrane potential, the fibroblast conductance G(f), and the MF gap-junctional coupling G(gap). Furthermore, we find that our MF composite can show autorhythmic and oscillatory behaviors in addition to an excitable response. Our 2D studies use (a) both homogeneous and inhomogeneous distributions of fibroblasts, (b) various ranges for parameters such as G(gap), G(f), and E-f, and (c) intercellular couplings that can be zero-sided, one-sided, and two-sided connections of fibroblasts with myocytes. We show, in particular, that the plane-wave conduction velocity CV decreases as a function of G(gap), for zero-sided and one-sided couplings; however, for two-sided coupling, CV decreases initially and then increases as a function of G(gap), and, eventually, we observe that conduction failure occurs for low values of G(gap). In our homogeneous studies, we find that the rotation speed and stability of a spiral wave can be controlled either by controlling G(gap) or E-f. Our studies with fibroblast inhomogeneities show that a spiral wave can get anchored to a local fibroblast inhomogeneity. We also study the efficacy of a low-amplitude control scheme, which has been suggested for the control of spiral-wave turbulence in mathematical models for cardiac tissue, in our MF model both with and without heterogeneities.

A Study of Early Afterdepolarizations in a Model for Human Ventricular Tissue

Sudden cardiac death is often caused by cardiac arrhythmias. Recently, special attention has been given to a certain arrhythmogenic condition, the long-QT syndrome, which occurs as a result of genetic mutations or drug toxicity. The underlying mechanisms of arrhythmias, caused by the long-QT syndrome, are not fully understood. However, arrhythmias are often connected to special excitations of cardiac cells, called early afterdepolarizations (EADs), which are depolarizations during the repolarizing phase of the action potential. So far, EADs have been studied mainly in isolated cardiac cells. However, the question on how EADs at the single-cell level can result in fibrillation at the tissue level, especially in human cell models, has not been widely studied yet. In this paper, we study wave patterns that result from single-cell EAD dynamics in a mathematical model for human ventricular cardiac tissue. We induce EADs by modeling experimental conditions which have been shown to evoke EADs at a single-cell level: by an increase of L-type Ca currents and a decrease of the delayed rectifier potassium currents. We show that, at the tissue level and depending on these parameters, three types of abnormal wave patterns emerge. We classify them into two types of spiral fibrillation and one type of oscillatory dynamics. Moreover, we find that the emergent wave patterns can be driven by calcium or sodium currents and we find phase waves in the oscillatory excitation regime. From our simulations we predict that arrhythmias caused by EADs can occur during normal wave propagation and do not require tissue heterogeneities. Experimental verification of our results is possible for experiments at the cell-culture level, where EADs can be induced by an increase of the L-type calcium conductance and by the application of I-Kr blockers, and the properties of the emergent patterns can be studied by optical mapping of the voltage and calcium.