Search

The voltage-gated K + channel plays a key role in atrial excitability, conducting the ultra-rapid rectifier K + current (I Kur ) and contributing to the repolarization of the atrial action potential. In this study, we examine its regulation by hydrogen sulfide (H 2 S) in HL-1 cardiomyocytes and in HEK293 cells expressing human Kv1.5. Pacing induced remodeling resulted in shorting action potential duration, enhanced both Kv1.5 channel and H 2 S producing enzymes protein expression in HL-1 cardiomyocytes. H 2 S supplementation reduced these remodeling changes and restored action potential duration through inhibition of Kv1.5 channel. H 2 S also inhibited recombinant hKv1.5, lead to nitric oxide (NO) mediated S-nitrosylation and activated endothelial nitric oxide synthase (eNOS) by increased phosphorylation of Ser1177, prevention of NO formation precluded these effects. Regulation of I kur by H 2 S has important cardiovascular implications and represents a novel and potential therapeutic target., (© 2023. The Author(s).)

Chronic exposure to low levels of Carbon Monoxide is associated with an increased risk of cardiac arrhythmia. Microelectrode recordings from rat and guinea pig single isolated ventricular myocytes exposed to CO releasing molecule CORM-2 and excited at 0.2/s show repolarisation changes that develop over hundreds of seconds: action potential prolongation by delayed repolarisation, EADs, multiple EADs and oscillations around the plateau, leading to irreversible repolarisation failure. The measured direct effects of CO on currents in these cells, and ion channels expressed in mammalian systems showed an increase in prolonged late Na + , and a decrease in the maximal T- and L-type Ca ++ . peak and late Na + , ultra-rapid delayed, delayed rectifier, and the inward rectifier K + currents. Incorporation of these CO induced changes in maximal currents in ventricular cell models; (Gattoni et al., J. Physiol., 2016, 594, 4193-4224) (rat) and (Luo and Rudy, Circ. Res., 1994, 74, 1071-1096) (guinea-pig) and human endo-, mid-myo- and epi-cardial (O'Hara et al., PLoS Comput. Biol., 2011, 7, e1002061) models, by changes in maximal ionic conductance reproduces these repolarisation abnormalities. Simulations of cell populations with Gaussian distributions of maximal conductance parameters predict a CO induced increase in APD and its variability. Incorporation of these predicted CO induced conductance changes in human ventricular cell electrophysiology into ventricular tissue and wall models give changes in indices for the probability of the initiation of re-entrant arrhythmia., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Al-Owais, Steele, Holden and Benson.)

Hand held ECG recorders are transforming the way we detect and diagnose heart rhythm disorders. The Kardia 6 L was launched in 2019 to detect and diagnose heart rhythm disorders recording a six lead (limb lead) ECG. Recording and analysis of precordial leads are currently not supported by the Kardia 6 L. In this study we aim to assess if reliable chest lead data can be obtained using a simple modification to the recording system., Competing Interests: Declaration of Competing Interest No authors have anything to disclose specific to this manuscript., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Background: Using cardiovascular magnetic resonance imaging (CMR), it is possible to detect diffuse fibrosis of the left ventricle (LV) in patients with atrial fibrillation (AF), which may be independently associated with recurrence of AF after ablation. By conducting CMR, clinical, electrophysiology and biomarker assessment we planned to investigate LV myocardial fibrosis in patients undergoing AF ablation., Methods: LV fibrosis was assessed by T1 mapping in 31 patients undergoing percutaneous ablation for AF. Galectin-3, coronary sinus type I collagen C terminal telopeptide (ICTP), and type III procollagen N terminal peptide were measured with ELISA. Comparison was made between groups above and below the median for LV extracellular volume fraction (ECV), followed by regression analysis., Results: On linear regression analysis LV ECV had significant associations with invasive left atrial pressure (Beta 0.49, P = 0.008) and coronary sinus ICTP (Beta 0.75, P < 0.001), which remained significant on multivariable regression., Conclusion: LV fibrosis in patients with AF is associated with left atrial pressure and invasively measured levels of ICTP turnover biomarker.

Cardioversion and defibrillation by a single high energy shock applied by myocardial or body surface electrodes is painful, causes long term tissue damage, and is associated with worsening long term outcomes, but is almost always required for treatment of ventricular fibrillation . As a initial step towards developing methods that can terminate ventricular arrhythmias painlessly, we aim to determine if pacing stimuli at a rate of 5/s applied via an implantable cardiac defibrillator (ICD) can modify human ventricular fibrillation. In 8 patients undergoing defibrillation testing of a new/exchanged intracardiac defibrillator, five seconds of pacing at five stimuli per second was applied during the 10-20 seconds of induced ventricular fibrillation before the defibrillation shock was automatically applied, and the cardiac electrograms recorded and analyzed. The high frequency pacing did not entrain the ventricular fibrillation, but altered the dominant frequency in all 8 patients, and modulated the phase computed via the Hilbert Transform, in four of the patients. In this pilot study we demonstrate that high frequency pacing applied via ICD electrodes during VF can alter the dominant frequency and modulate the probability density of the phase of the electrogram of the ventricular fibrillation.

The human heart develops through complex mechanisms producing morphological and functional changes during gestation. We have recently demonstrated using diffusion tensor MRI that over the relatively short space of 40 days, between 100-140 days gestational age, the ventricular myocardium transforms from a disorganised tissue to the ordered structure characteristic of mature cardiac tissue. However, the genetic basis underpinning this maturation is unclear. Herein, we have used RNA-Seq to establish the developmentally-regulated transcriptome of gene expression in the developing human heart across three gestational ages in the first and second trimester. By comparing 9 weeks gestational age (WGA) with 12 WGA, we find 288 genes show significant differential expression. 305 genes were significantly altered comparing 12 and 16 WGA, and 806 genes differentially expressed between 9 and 16 WGA. Network analysis was used to identify genetic interactions, node properties and gene ontology categories. In summary, we present a comprehensive transcriptomic analysis of human heart development during early gestation, and identify differentially expressed genes during heart development between 9 and 16 weeks, overlapping the first and early second trimester.

Aberrant uterine myometrial activities in humans are major health issues. However, the cellular and tissue mechanism(s) that maintain the uterine myometrium at rest during gestation, and that initiate and maintain long-lasting uterine contractions during delivery are incompletely understood. In this study we construct a computational model for describing the electrical activity (simple and complex action potentials), intracellular calcium dynamics and mechanical contractions of isolated uterine myocytes from the pregnant rat. The model reproduces variant types of action potentials - from spikes with a smooth plateau, to spikes with an oscillatory plateau, to bursts of spikes - that are seen during late gestation under different physiological conditions. The effects of the hormones oestradiol (via reductions in calcium and potassium selective channel conductance), oxytocin (via an increase in intracellular calcium release) and the tocolytic nifedipine (via a block of L-type calcium channels currents) on action potentials and contractions are also reproduced, which quantitatively match to experimental data. All of these results validated the cell model development. In conclusion, the developed model provides a computational platform for further investigations of the ionic mechanism underlying the genesis and control of electrical and mechanical activities in the rat uterine myocytes.

Aims: The identification of arrhythmogenic right ventricular dysplasia (ARVD) from 12-channel standard electrocardiogram (ECG) is challenging. High density ECG data may identify lead locations and criteria with a higher sensitivity., Methods and Results: Eighty-channel ECG recording from patients diagnosed with ARVD and controls were quantified by magnitude and integral measures of QRS and T waves and by a measure (the average silhouette width) of differences in the shapes of the normalized ECG cycles. The channels with the best separability between ARVD patients and controls were near the right ventricular wall, at the third intercostal space. These channels showed pronounced differences in P waves compared to controls as well as the expected differences in QRS and T waves., Conclusion: Multichannel recordings, as in body surface mapping, add little to the reliability of diagnosing ARVD from ECGs. However, repositioning ECG electrodes to a high anterior position can improve the identification of ECG variations in ARVD. Additionally, increased P wave amplitude appears to be associated with ARVD.

Aims: To test the ability of four circulating biomarkers of fibrosis, and of low left atrial voltage, to predict recurrence of atrial fibrillation after catheter ablation., Background: Circulating biomarkers potentially may be used to improve patient selection for atrial fibrillation ablation. Low voltage areas in the left atrium predict arrhythmia recurrence when mapped in sinus rhythm. This study tested type III procollagen N terminal peptide (PIIINP), galectin-3 (gal-3), fibroblast growth factor 23 (FGF-23), and type I collagen C terminal telopeptide (ICTP), and whether low voltage areas in the left atrium predicted atrial fibrillation recurrence, irrespective of the rhythm during mapping., Methods: 92 atrial fibrillation ablation patients were studied. Biomarker levels in peripheral and intra-cardiac blood were measured with enzyme-linked immunosorbent assay. Low voltage (<0.5mV) was expressed as a proportion of the mapped left atrial surface area. Follow-up was one year. The primary endpoint was recurrence of arrhythmia. The secondary endpoint was a composite of recurrence despite two procedures, or after one procedure if no second procedure was undertaken., Results: The biomarkers were not predictive of either endpoint. After multivariate Cox regression analysis, high proportion of low voltage area in the left atrium was found to predict the primary endpoint in sinus rhythm mapping (hazard ratio 4.323, 95% confidence interval 1.337-13.982, p = 0.014) and atrial fibrillation mapping (hazard ratio 5.195, 95% confidence interval 1.032-26.141, p = 0.046). This effect was also apparent for the secondary endpoint., Conclusion: The studied biomarkers do not predict arrhythmia recurrence after catheter ablation. Left atrial voltage is an independent predictor of recurrence, whether the left atrium is mapped in atrial fibrillation or sinus rhythm.

Aims: Measurement of circulating biomarkers of fibrosis may have a role in selecting patients and treatment strategy for catheter ablation. Pro-collagen type III N-terminal pro-peptide (PIIINP), C-telopeptide of type I collagen (ICTP), fibroblast growth factor 23 (FGF-23), and galectin 3 (gal-3) have all been suggested as possible biomarkers for this indication, but studies assessing whether peripheral levels reflect intra-cardiac levels are scarce., Methods and Results: We studied 93 patients undergoing ablation for paroxysmal atrial fibrillation (AF) (n = 63) or non-paroxysmal AF (n = 30). Femoral venous, left and right atrial, and coronary sinus blood were analysed using ELISA to determine biomarker levels. Levels were compared with control patients (n = 36) and baseline characteristics, including left atrial voltage mapping data. C-telopeptide of type I collagen levels were higher in AF than in non-AF patients (P = 0.007). Peripheral ICTP levels were higher than all intra-cardiac levels (P < 0.001). Peripheral gal-3 levels were higher than left atrial levels (P = 0.001). Peripheral levels of FGF-23 and PIIINP were not significantly different from intra-cardiac levels. CS levels of ICTP were higher than right and left atrial levels (P < 0.001). gal-3 was higher in women vs. men (P ≤ 0.001) and with higher body mass index (P ≤ 0.001). ICTP levels increased with reducing ejection fraction (P ≤ 0.012)., Conclusions: Atrial fibrillation patients have higher levels of circulating ICTP than matched non-AF controls. In AF ablation patients, intra-cardiac sampling of FGF-23 or PIIINP gives no further information over peripheral sampling. For gal-3 and ICTP, intra-cardiac sampling may be necessary to assess their association with intra-cardiac processes. None of the biomarkers is related to fibrosis assessed by left atrial voltage., (Published on behalf of the European Society of Cardiology. All rights reserved. © The Author 2017. For Permissions, please email: journals.permissions@oup.com.)

Background: Prolongation of the QT interval of the electrocardiogram (ECG), underlain by prolongation of the action potential duration (APD) at the cellular level, is linked to increased vulnerability to cardiac arrhythmia. Pharmacological management of arrhythmia associated with QT prolongation is typically achieved through attempting to restore APD to control ranges, reversing the enhanced vulnerability to Ca 2+ -dependent afterdepolarisations (arrhythmia triggers) and increased transmural dispersion of repolarisation (arrhythmia substrate) associated with APD prolongation. However, such pharmacological modulation has been demonstrated to have limited effectiveness. Understanding the integrative functional impact of pharmacological modulation requires simultaneous investigation of both the trigger and substrate. Methods: We implemented a multi-scale (cell and tissue) in silico approach using a model of the human ventricular action potential, integrated with a model of stochastic 3D spatiotemporal Ca 2+ dynamics, and parameter modification to mimic prolonged QT conditions. We used these models to examine the efficacy of the hERG activator MC-II-157c in restoring APD to control ranges, examined its effects on arrhythmia triggers and substrates, and the interaction of these arrhythmia triggers and substrates. Results: QT prolongation conditions promoted the development of spontaneous release events underlying afterdepolarisations during rapid pacing. MC-II-157c applied to prolonged QT conditions shortened the APD, inhibited the development of afterdepolarisations and reduced the probability of afterdepolarisations manifesting as triggered activity in single cells. In tissue, QT prolongation resulted in an increased transmural dispersion of repolarisation, which manifested as an increased vulnerable window for uni-directional conduction block. In some cases, MC-II-157c further increased the vulnerable window through its effects on I Na . The combination of stochastic release event modulation and transmural dispersion of repolarisation modulation by MC-II-157c resulted in an integrative behavior wherein the arrhythmia trigger is reduced but the arrhythmia substrate is increased, leading to variable and non-linear overall vulnerability to arrhythmia. Conclusion: The relative balance of reduced trigger and increased substrate underlies a multi-dimensional role of MC-II-157c in modulation of cardiac arrhythmia vulnerability associated with prolonged QT interval.

The developmental timeline of the human heart remains elusive. The heart takes on its characteristic four chambered appearance by ~56 days gestational age (DGA). However, owing to the complexities (both technical and logistical) of exploring development in utero, we understand little of how the ventricular walls develop. To address this, we employed diffusion tensor magnetic resonance imaging to explore the architecture and tissue organization of the developing heart aged 95-143 DGA. We show that fractional anisotropy increases (from ~0.1 to ~0.5), diffusion coefficients decrease (from ~1 × 10 -3 mm 2 /sec to ~0.4 × 10 -3 mm 2 /sec), and fiber paths, extracted by tractography, increase linearly with gestation, indicative of the increasing organization of the ventricular myocytes. By 143 DGA, the developing heart has the classical helical organization observed in mature mammalian tissue. This was accompanied by an increase in connexin 43 and connexin 40 expression levels, suggesting their role in the development of the ventricular conduction system and that electrical propagation across the heart is facilitated in later gestation. Our findings highlight a key developmental window for the structural organization of the fetal heart.

Rhythm control of atrial fibrillation (AF) remains challenging, with modest long-term success rates. Atrial fibrosis has been associated with AF, but the clinical utility of assessment of this fibrosis has yet to be fully elucidated. In this paper we review the current state of understanding of the pathophysiology of atrial fibrosis in AF, and its impact upon the instigation and propagation of the arrhythmia. Fibrosis causes an increase in volume of dysfunctional extracellular matrix, and is associated with cellular alterations such as hypertrophy, apoptosis and membrane dysfunction within the atrial myocardium. In turn, these cause pathological alterations to atrial conduction, such as increased anisotropy, conduction block and re-entry, which can lead to AF. We review current methods of assessing atrial fibrosis and their impact upon the prediction of success of interventional rhythm control strategies such as ablation and cardioversion. We focus particularly on circulating biomarkers of fibrosis and scar formation; their role in the fibrotic process, and their value in the prediction of rhythm control success. We also review imaging and invasive electrocardiographic mapping techniques that may identify fibrosis, and again assess their potential predictive value. In this area there exist many unanswered questions, but further work will help to refine techniques to reliably identify and treat those patients who are most likely to benefit from rhythm control treatment strategies., (Copyright © 2016 Elsevier Ireland Ltd. All rights reserved.)

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., (Published on behalf of the European Society of Cardiology. All rights reserved. © The Author 2014. For permissions please email: journals.permissions@oup.com.)

Aims: We aim to engineer a computational model of propagation during normal sinus rhythm in the foetal human heart, by modifying models for adult cardiac tissue to match foetal electrocardiogram (fECG) characteristics. The model will be partially validated by fECG data, and applied to explore possible mechanisms of arrhythmogenesis in the foetal heart., Methods and Results: Foetal electrocardiograms have been recorded during pregnancy, with P- and T-waves, and the QRS complex, identified by averaging and signal processing. Intervals of the fECG are extracted and used to modify currently available human adult cardiomyocyte models. RR intervals inform models of the pacemaking cells by constraining their rate, the QT interval and its rate dependence constrain models of ventricular cells, and the width of the P-wave, the QR and PR intervals constrain propagation times, conduction velocities, and intercellular coupling. These cell models are coupled into a one-dimensional (1D) model of propagation during normal sinus rhythm in the human foetal heart. We constructed a modular, heterogeneous 1D model for propagation in the foetal heart, and predicted the effects of reduction in L-type Ca(++) current. These include bradycardia and atrioventricular conduction blocks. These may account quantitatively for congenital heart block produced by positive IgG antibodies., Conclusion: The fECG can be interpreted mechanistically and quantitatively by using a simple computational model for propagation. After further validation, by clinical recordings of the fECG and the electrophysiological experiments on foetal cardiac cells and tissues, the model may be used to predict the effects of maternally administered pharmaceuticals on the fECG.

The mechanisms leading to the initiation of normal, premature or dysfunctional human labour are poorly understood, as animal models are inappropriate, and experimental studies are limited. Computational modelling provides a means of linking non-invasive clinical data with the results of in vitro cell and tissue physiology. Nonlinear wave processes - propagation in an excitable medium - provides a quantitatively testable description of mechanisms of premature and full term labour, and a view of changes in uterine electrophysiology during gestation as a trajectory in excitation and intercellular coupling parameter space. Propagation phenomena can account for both premature and full term labour., (Copyright © 2013 Elsevier Ireland Ltd. All rights reserved.)

We construct the components for a family of computational models of the electrophysiology of the human foetal heart from 60 days gestational age (DGA) to full term. This requires both cell excitation models that reconstruct the myocyte action potentials, and datasets of cardiac geometry and architecture. Fast low-angle shot and diffusion tensor magnetic resonance imaging (DT-MRI) of foetal hearts provides cardiac geometry with voxel resolution of approximately 100 µm. DT-MRI measures the relative diffusion of protons and provides a measure of the average intravoxel myocyte orientation, and the orientation of any higher order orthotropic organization of the tissue. Such orthotropic organization in the adult mammalian heart has been identified with myocardial sheets and cleavage planes between them. During gestation, the architecture of the human ventricular wall changes from being irregular and isotropic at 100 DGA to an anisotropic and orthotropic architecture by 140 DGA, when it has the smooth, approximately 120° transmural change in myocyte orientation that is characteristic of the adult mammalian ventricle. The DT obtained from DT-MRI provides the conductivity tensor that determines the spread of potential within computational models of cardiac tissue electrophysiology. The foetal electrocardiogram (fECG) can be recorded from approximately 60 DGA, and RR, PR and QT intervals between the P, R, Q and T waves of the fECG can be extracted by averaging from approximately 90 DGA. The RR intervals provide a measure of the pacemaker rate, the QT intervals an index of ventricular action potential duration, and its rate-dependence, and so these intervals constrain and inform models of cell electrophysiology. The parameters of models of adult human sinostrial node and ventricular cells that are based on adult cell electrophysiology and tissue molecular mapping have been modified to construct preliminary models of foetal cell electrophysiology, which reproduce these intervals from fECG recordings. The PR and QR intervals provide an index of conduction times, and hence propagation velocities (approx. 1-10 cm s(-1), increasing during gestation) and so inform models of tissue electrophysiology. Although the developing foetal heart is small and the cells are weakly coupled, it can support potentially lethal re-entrant arrhythmia.

Micro-computed tomography (micro-CT) has been widely used to generate high-resolution 3-D tissue images from small animals nondestructively, especially for mineralized skeletal tissues. However, its application to the analysis of soft cardiovascular tissues has been limited by poor inter-tissue contrast. Recent ex vivo studies have shown that contrast between muscular and connective tissue in micro-CT images can be enhanced by staining with iodine. In the present study, we apply this novel technique for imaging of cardiovascular structures in canine hearts. We optimize the method to obtain high-resolution X-ray micro-CT images of the canine atria and its distinctive regions-including the Bachmann's bundle, atrioventricular node, pulmonary arteries and veins-with clear inter-tissue contrast. The imaging results are used to reconstruct and segment the detailed 3-D geometry of the atria. Structure tensor analysis shows that the arrangement of atrial fibers can also be characterized using the enhanced micro-CT images, as iodine preferentially accumulates within the muscular fibers rather than in connective tissues. This novel technique can be particularly useful in nondestructive imaging of 3-D cardiac architectures from large animals and humans, due to the combination of relatively high speed ( ~ 1 h/per scan of the large canine heart) and high voxel resolution (36 μm) provided. In summary, contrast micro-CT facilitates fast and nondestructive imaging and segmenting of detailed 3-D cardiovascular geometries, as well as measuring fiber orientation, which are crucial in constructing biophysically detailed computational cardiac models.

Computational models of human atrial cells, tissues and atria have been developed. Cell models, for atrial wall, crista terminalis, appendage, Bachmann's bundle and pectinate myocytes are characterised by action potentials, ionic currents and action potential duration (APD) restitution. The principal effect of the ion channel remodelling of persistent atrial fibrillation (AF), and a mutation producing familial AF, was APD shortening at all rates. Electrical alternans was abolished by the modelled action of Dronedarone. AF induced gap junctional remodelling slows propagation velocity at all rates. Re-entrant spiral waves in 2-D models are characterised by their frequency, wavelength, meander and stability. For homogenous models of normal tissue, spiral waves self-terminate, due to meander to inexcitable boundaries, and by dissipation of excitation. AF electrical remodelling in these homogenous models led to persistence of spiral waves, and AF fibrotic remodelling to their breakdown into fibrillatory activity. An anatomical model of the atria was partially validated by the activation times of normal sinus rhythm. The use of tissue geometry from clinical MRI, and tissue anisotropy from ex vivo diffusion tensor magnetic resonance imaging is outlined. In the homogenous model of normal atria, a single scroll breaks down onto spatio-temporal irregularity (electrical fibrillation) that is self-terminating; while in the AF remodelled atria the fibrillatory activity is persistent. The persistence of electrical AF can be dissected in the model in terms of ion channel and intercellular coupling processes, that can be modified pharmacologically; the effects of anatomy, that can be modified by ablation; and the permanent effects of fibrosis, that need to be prevented., (Copyright © 2011 Elsevier B.V. All rights reserved.)

Introduction: β-adrenergic stimulation increases the heart rate by accelerating the electrical activity of the pacemaker of the heart, the sinoatrial node (SAN). Ionic mechanisms underlying the actions of β-adrenergic stimulation are not yet fully understood. Isoprenaline (ISO), a β-adrenoceptor agonist, shifts voltage-dependent I(f) activation to more positive potentials resulting in an increase of I(f), which has been suggested to be the main mechanism underlying the effect of β-adrenergic stimulation. However, ISO has been found to increase the firing rate of rabbit SAN cells when I(f) is blocked. ISO also increases I(CaL), I(st), I(Kr), and I(Ks); and shifts the activation of I(Kr) to more negative potentials and increases the rate of its deactivation. ISO has also been reported to increase the intracellular Ca(2+) transient, which can contribute to chronotropy by modulating the "Ca(2+) clock." The aim of this study was to analyze the ionic mechanisms underlying the positive chronotropy of β-adrenergic stimulation using two distinct and well established computational models of the electrical activity of rabbit SAN cells., Methods and Results: We modified the Boyett et al. (2001) and Kurata et al. (2008) models of electrical activity for the central and peripheral rabbit SAN cells by incorporating equations for the known dose-dependent actions of ISO on various ionic channel currents (I(CaL), I(st), I(Kr), and I(Ks)), kinetics of I(Kr) and I(f), and the intracellular Ca(2+) transient. These equations were constructed from experimental data. To investigate the ionic basis of the effects of ISO, we simulated the chronotropic effect of a range of ISO concentrations when ISO exerted all its actions or just a subset of them., Conclusion: In both the Boyett et al. and Kurata et al. SAN models, the chronotropic effect of ISO was found to result from an integrated action of ISO on I(CaL), I(f), I(st), I(Kr), and I(Ks), among which an increase in the rate of deactivation of I(Kr) plays a prominent role, though the effect of ISO on I(f) and [Ca(2+)](i) also plays a role.

It has been shown by histology that cardiac myocytes are organized into laminae and this structure is important in function, both influencing the spread of electrical activation and enabling myocardial thickening in systole by laminar sliding. We have carried out high-spatial resolution three-dimensional MRI of the ventricular myolaminae of the entire volume of the isolated rat heart after contrast perfusion [dimeglumine gadopentate (Gd-DTPA)]. Four ex vivo rat hearts were perfused with Gd-DTPA and fixative and high-spatial resolution MRI was performed on a 9.4T MRI system. After MRI, cryosectioning followed by histology was performed. Images from MRI and histology were aligned, described, and quantitatively compared. In the three-dimensional MR images we directly show the presence of laminae and demonstrate that these are highly branching and are absent from much of the subepicardium. We visualized these MRI volumes to demonstrate laminar architecture and quantitatively demonstrated that the structural features observed are similar to those imaged in histology. We showed qualitatively and quantitatively that laminar architecture is similar in the four hearts. MRI can be used to image the laminar architecture of ex vivo hearts in three dimensions, and the images produced are qualitatively and quantitatively comparable with histology. We have demonstrated in the rat that: 1) laminar architecture is consistent between hearts; 2) myolaminae are absent from much of the subepicardium; and 3) although localized orthotropy is present throughout the myocardium, tracked myolaminae are branching structures and do not have a discrete identity.

We apply virtual tissue engineering to the full term human uterus with a view to reconstruction of the spatiotemporal patterns of electrical activity of the myometrium that control mechanical activity via intracellular calcium. The three-dimensional geometry of the gravid uterus has been reconstructed from segmented in vivo magnetic resonance imaging as well as ex vivo diffusion tensor magnetic resonance imaging to resolve fine scale tissue architecture. A late-pregnancy uterine smooth muscle cell model is constructed and bursting analysed using continuation algorithms. These cell models are incorporated into partial differential equation models for tissue synchronisation and propagation. The ultimate objective is to develop a quantitative and predictive understanding of the mechanisms that initiate and regulate labour., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Despite a vast amount of experimental and clinical data on the underlying ionic, cellular and tissue substrates, the mechanisms of common atrial arrhythmias (such as atrial fibrillation, AF) arising from the functional interactions at the whole atria level remain unclear. Computational modelling provides a quantitative framework for integrating such multi-scale data and understanding the arrhythmogenic behaviour that emerges from the collective spatio-temporal dynamics in all parts of the heart. In this study, we have developed a multi-scale hierarchy of biophysically detailed computational models for the human atria--the 3D virtual human atria. Primarily, diffusion tensor MRI reconstruction of the tissue geometry and fibre orientation in the human sinoatrial node (SAN) and surrounding atrial muscle was integrated into the 3D model of the whole atria dissected from the Visible Human dataset. The anatomical models were combined with the heterogeneous atrial action potential (AP) models, and used to simulate the AP conduction in the human atria under various conditions: SAN pacemaking and atrial activation in the normal rhythm, break-down of regular AP wave-fronts during rapid atrial pacing, and the genesis of multiple re-entrant wavelets characteristic of AF. Contributions of different properties of the tissue to mechanisms of the normal rhythm and arrhythmogenesis were investigated. Primarily, the simulations showed that tissue heterogeneity caused the break-down of the normal AP wave-fronts at rapid pacing rates, which initiated a pair of re-entrant spiral waves; and tissue anisotropy resulted in a further break-down of the spiral waves into multiple meandering wavelets characteristic of AF. The 3D virtual atria model itself was incorporated into the torso model to simulate the body surface ECG patterns in the normal and arrhythmic conditions. Therefore, a state-of-the-art computational platform has been developed, which can be used for studying multi-scale electrical phenomena during atrial conduction and AF arrhythmogenesis. Results of such simulations can be directly compared with electrophysiological and endocardial mapping data, as well as clinical ECG recordings. The virtual human atria can provide in-depth insights into 3D excitation propagation processes within atrial walls of a whole heart in vivo, which is beyond the current technical capabilities of experimental or clinical set-ups., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Atrial arrhythmias resulting from abnormally rapid focal activity in the atria may be reflected in an altered P-wave morphology (PWM) in the ECG. Although clinically important, detailed relationships between PWM and origins of atrial focal excitations have not been established. To study such relationships, we developed computational models of the human atria and torso. The model simulation results were used to evaluate an extant clinical algorithm for locating the origin of atrial focal points from the ECG. The simulations showed that the algorithm was practical and could predict the atrial focal locations with 85% accuracy. We proposed a further refinement of the algorithm to distinguish between focal locations within the large atrial bundles.

Mutations to hERG which result in changes to the rapid delayed rectifier current I(Kr) can cause long and short QT syndromes and are associated with an increased risk of cardiac arrhythmias. Experimental recordings of I(Kr) reveal the effects of mutations at the channel level, but how these changes translate to the cell and tissue levels remains unclear. We used computational models of human ventricular myocytes and tissues to predict and quantify the effects that de novo hERG mutations would have on cell and tissue electrophysiology. Mutations that decreased I(Kr) maximum conductance resulted in an increased cell and tissue action potential duration (APD) and a long QT interval on the electrocardiogram (ECG), whereas those that caused a positive shift in the inactivation curve resulted in a decreased APD and a short QT. Tissue vulnerability to re-entrant arrhythmias was correlated with transmural dispersion of repolarisation, and any change to this vulnerability could be inferred from the ECG QT interval or T wave peak-to-end time. Faster I(Kr) activation kinetics caused cell APD alternans to appear over a wider range of pacing rates and with a larger magnitude, and spatial heterogeneity in these cellular alternans resulted in discordant alternans at the tissue level. Thus, from channel kinetic data, we can predict the tissue-level electrophysiological effects of any hERG mutations and identify how the mutation would manifest clinically, as either a long or short QT syndrome with or without an increased risk of alternans and re-entrant arrhythmias.

Uterine contractions during labor are discretely regulated by rhythmic action potentials (AP) of varying duration and form that serve to determine calcium-dependent force production. We have employed a computational biology approach to develop a fuller understanding of the complexity of excitation-contraction (E-C) coupling of uterine smooth muscle cells (USMC). Our overall aim is to establish a mathematical platform of sufficient biophysical detail to quantitatively describe known uterine E-C coupling parameters and thereby inform future empirical investigations of physiological and pathophysiological mechanisms governing normal and dysfunctional labors. From published and unpublished data we construct mathematical models for fourteen ionic currents of USMCs: Ca2+ currents (L- and T-type), Na+ current, an hyperpolarization-activated current, three voltage-gated K+ currents, two Ca2+-activated K+ current, Ca2+-activated Cl current, non-specific cation current, Na+-Ca2+ exchanger, Na+-K+ pump and background current. The magnitudes and kinetics of each current system in a spindle shaped single cell with a specified surface area:volume ratio is described by differential equations, in terms of maximal conductances, electrochemical gradient, voltage-dependent activation/inactivation gating variables and temporal changes in intracellular Ca2+ computed from known Ca2+ fluxes. These quantifications are validated by the reconstruction of the individual experimental ionic currents obtained under voltage-clamp. Phasic contraction is modeled in relation to the time constant of changing [Ca2+]i. This integrated model is validated by its reconstruction of the different USMC AP configurations (spikes, plateau and bursts of spikes), the change from bursting to plateau type AP produced by estradiol and of simultaneous experimental recordings of spontaneous AP, [Ca2+]i and phasic force. In summary, our advanced mathematical model provides a powerful tool to investigate the physiological ionic mechanisms underlying the genesis of uterine electrical E-C coupling of labor and parturition. This will furnish the evolution of descriptive and predictive quantitative models of myometrial electrogenesis at the whole cell and tissue levels.

Reaction-diffusion computational models of cardiac electrophysiology require both dynamic excitation models that reconstruct the action potentials of myocytes as well as datasets of cardiac geometry and architecture that provide the electrical diffusion tensor D, which determines how excitation spreads through the tissue. We illustrate an experimental pipeline we have developed in our laboratories for constructing and validating such datasets. The tensor D changes with location in the myocardium, and is determined by tissue architecture. Diffusion tensor magnetic resonance imaging (DT-MRI) provides three eigenvectors e(i) and eigenvalues λ(i) at each voxel throughout the tissue that can be used to reconstruct this architecture. The primary eigenvector e(1) is a histologically validated measure of myocyte orientation (responsible for anisotropic propagation). The secondary and tertiary eigenvectors (e(2) and e(3)) specify the directions of any orthotropic structure if λ(2) is significantly greater than λ(3)-this orthotropy has been identified with sheets or cleavage planes. For simulations, the components of D are scaled in the fibre and cross-fibre directions for anisotropic simulations (or fibre, sheet and sheet normal directions for orthotropic tissues) so that simulated conduction velocities match values from optical imaging or plunge electrode experiments. The simulated pattern of propagation of action potentials in the models is partially validated by optical recordings of spatio-temporal activity on the surfaces of hearts. We also describe several techniques that enhance components of the pipeline, or that allow the pipeline to be applied to different areas of research: Q ball imaging provides evidence for multi-modal orientation distributions within a fraction of voxels, infarcts can be identified by changes in the anisotropic structure-irregularity in myocyte orientation and a decrease in fractional anisotropy, clinical imaging provides human ventricular geometry and can identify ischaemic and infarcted regions, and simulations in human geometries examine the roles of anisotropic and orthotropic architecture in the initiation of arrhythmias.

The electrical activity of cardiac and uterine tissues has been reconstructed by detailed computer models in the form of virtual tissues. Virtual tissues are biophysically and anatomically detailed, and represent quantitatively predictive models of the physiological and pathophysiological behaviours of tissue within an isolated organ. The cell excitation properties are quantitatively reproduced by equations that describe the kinetics of a few dozen proteins. These equations are derived from experimental measurements of membrane potentials, ionic currents, fluxes, and concentrations. Some of the measurements were taken from human cells and human ion channel proteins expressed in non-human cells, but they were mostly taken from cells of other animal species. Data on tissue geometry and architecture are obtained from the diffusion tensor magnetic resonance imaging of ex vivo or post mortem tissue, and are used to compute the spread of current in the tissue. Cardiac virtual tissues are well established and reproduce normal and pathological patterns of cardiac excitation within the atria or ventricles of the human heart. They have been applied to increase the understanding of normal cardiac electrophysiology, to evaluate the candidate mechanisms for re-entrant arrhythmias that lead to sudden cardiac death, and to predict the tissue level effects of mutant or pharmacologically-modified ion channels. The human full-term virtual uterus is still in development. This virtual tissue reproduces the in vitro behaviour of uterine tissue biopsies, and provides possible mechanisms for premature labour., (2010 FRAME.)

Intracellular Ca(2+) dynamics provides excitation-contraction coupling in cardiac myocytes. Under pathological conditions, spontaneous Ca(2+) release events can lead to intracellular Ca(2+) travelling waves, which can break, giving transitory or persistent intracellular re-entrant Ca(2+) scroll waves. Intracellular Ca(2+) waves can trigger cellular delayed after-depolarizations of membrane potential, which if they occur in a cluster of a few hundred neighbouring myocytes may lead to cardiac arrhythmia. Quantitative prediction of the initiation and propagation of intracellular Ca(2+) waves requires the dynamics of Ca(2+)-induced Ca(2+) release, and the intracellular spatial distribution of Ca(2+) release units (CRUs). The spatial distribution of ryanodine receptor clusters within a few sarcomeres was reconstructed directly from confocal imaging measurements. It was then embedded into a three-dimensional ventricular cell model, with a resting membrane potential and simple stochastic Ca(2+)-induced Ca(2+) release dynamics. Isotropic global Ca(2+) wave propagation can be produced within the anisotropic intracellular architecture, by isotropic local Ca(2+) diffusion, and the branching Z-disc structure providing inter Z-disc pathways for Ca(2+) propagation. The branching Z-disc provides a broader spatial distribution of ryanodine receptor clusters across Z-discs, which reduces the likelihood of wave initiation by spontaneous Ca(2+) releases. Intracellular Ca(2+) dynamics during catecholaminergic polymorphic ventricular tachycardia (CPVT) was simulated phenomenologically by increasing the Ca(2+) sensitivity factor of the CRU, which results in an increased rate of Ca(2+) release events. Flecainide has been shown to prevent arrhythmias in a murine model of CPVT and in patients. The modelled actions of flecainide on the time course of Ca(2+) release events prevented the initiation of Ca(2+) waves.

Over the last decade, magnetic resonance imaging (MRI) has become established as a useful referral diagnostic method in veterinary medicine that is widely used in small animal brain and spinal diseases, aural, nasal and orbital disorders, planning soft tissue surgery, oncology and small animal and equine orthopaedics. The use of MRI in these disciplines has grown due to its unparalleled capability to image soft tissue structures. This has been exploited in human cardiology where, despite the inherent difficulties in imaging a moving, contractile structure, cardiac MRI (CMRI) has become the optimal technique for the morphological assessment and quantification of ventricular function. Both CMRI hardware and software systems have developed rapidly in the last 10 years but although several preliminary veterinary CMRI studies have been reported, the technique's growth has been limited and is currently used primarily in clinical research. A review of published studies is presented with a description of CMRI technology and the potential of CMRI is discussed along with some of the reasons for its limited usage., (Copyright (c) 2008 Elsevier Ltd. All rights reserved.)

Fluctuations in the time interval between two consecutive R-waves of electrocardiogram during normal sinus rhythm may result from irregularities in the autonomic drive of the pacemaking sinoatrial node (SAN). We use a biophysically detailed mathematical model of the action potentials of rabbit SAN to quantify the effects of fluctuations in acetylcholine (ACh) on the pacemaker activity of the SAN and its variability. Fluctuations in ACh concentration model the effect of stochastic activity in the vagal parasympathetic fibers that innervate the SAN and produce varying rates of depolarization during the pacemaker potential, leading to fluctuations in cycle length (CL). Both the estimated maximal Lyapunov exponent and the noise limit of the resultant sequence of fluctuating CLs suggest chaotic dynamics. Apparently chaotic heart rate variability (HRV) seen in sinus rhythm can be produced by stochastic modulation of the SAN. The identification of HRV data as chaotic by use of time series measures such as a positive maximal Lyapunov exponent or positive noise limit requires both caution and a quantitative, predictive mechanistic model that is fully deterministic.

Quantitative understanding of the Ca(2+) handling in cardiac ventricular myocytes requires accurate knowledge of cardiac ultrastructure and protein distribution. We have therefore developed high-resolution imaging and analysis approaches to measure the three-dimensional distribution of immunolabelled proteins with confocal microscopy. Labelling of single rat cardiac myocytes with an antibody to the Z-line marker alpha-actinin revealed a complex architecture of sarcomere misalignment across single cells. Double immunolabelling was used to relate the Z-line structure to the distribution of ryanodine receptors (RyRs, the intracellular Ca(2+) release channels) and the transverse tubular system. Both RyR and transverse tubular system distributions exhibited frequent dislocations from the simple planar geometry generally assumed in existing mathematical models. To investigate potential effects of these irregularities on Ca(2+) dynamics, we determined the three-dimensional distribution of RyR clusters within an extended section of a single rat ventricular myocyte to construct a model of stochastic Ca(2+) dynamics with a measured Ca(2+) release unit (CRU) distribution. Calculations with this model were compared with a second model in which all CRUs were placed on flat planes. The model with a realistic CRU distribution supported Ca(2+) waves that spread axially along the cell at velocities of approximately 50 mum s(-1). By contrast, in the model with planar CRU distribution the axial wave spread was slowed roughly twofold and wave propagation often nearly faltered. These results demonstrate that spatial features of the CRU distribution on multiple length scales may significantly affect intracellular Ca(2+) dynamics and must be captured in detailed mechanistic models to achieve quantitative as well as qualitative insight.

Atrial fibrillation (AF) has been linked to increased inward rectifier potassium current, I(K1), either due to AF-induced electrical remodelling, or from functional changes due to the Kir2.1 V93I mutation. The aim of this simulation study was to identify at cell and tissue levels' mechanisms by which increased I(K1) facilitates and perpetuates AF. The Courtemanche et al. human atrial cell action potential (AP) model was modified to incorporate reported changes in I(K1) induced by the Kir2.1 V93I mutation in both heterozygous (Het) and homozygous (Hom) mutant forms. The modified models for wild type (WT), Het and Hom conditions were incorporated into homogeneous 1D, 2D and 3D tissue models. Restitution curves of AP duration (APD), effective refractory period (ERP) and conduction velocity (CV) were computed and both the temporal and the spatial vulnerability of atrial tissue to re-entry were measured. The lifespan and tip meandering pattern of re-entry were also characterised. For comparison, parallel simulations were performed by incorporating into the Courtmanche et al. model a linear increase in maximal I(K1) conductance. It was found that the gain-in-function of V93I 'mutant'I(K1) led to abbreviated atrial APs and flattened APD, ERP and CV restitution curves. It also hyperpolarised atrial resting membrane potential and slowed down intra-atrial conduction. V93I 'mutant'I(K1) reduced the tissue's temporal vulnerability but increased spatial vulnerability to initiate and sustain re-entry, resulting in an increased overall susceptibility of atrial tissue to arrhythmogenesis. In the 2D model, spiral waves self-terminated for WT (lifespan < 3.3 s) tissue, but persisted in Het and Hom tissues for the whole simulation period (lifespan > 10 s). The tip of the spiral wave meandered more in WT tissue than in Het and Hom tissues. Increased I(K1) due to augmented maximal conductance produced similar results to those of Het and Hom Kir2.1 V93I mutant conditions. In the 3D model the dynamic behaviour of scroll waves was stabilized by increased I(K1). In conclusion, increased I(K1) current, either by the Kir2.1 V93I mutation or by augmented maximal conductance, increases atrial susceptibility to arrhythmia by increasing the lifespan of re-entrant spiral waves and the stability of scroll waves in 3D tissue, thereby facilitating initiation and maintenance of re-entrant circuits.

Background: Dysfunction or pharmacologic inhibition of repolarizing cardiac ionic currents can lead to fatal arrhythmias. The hERG potassium channel underlies the repolarizing current I(Kr), and mutations therein can produce both long and short QT syndromes (LQT2 and SQT1). We previously reported on the diphenylurea compound NS1643, which acts on hERG channels in two distinct ways: by increasing overall conductance and by shifting the inactivation curve in the depolarized direction., Objective: The purpose of this study was to determine which of the two components contributes more to the antiarrhythmic effects of NS1643 under normokalemic and hypokalemic conditions., Methods: The study consisted of mathematical simulation of action potentials in a human ventricular ionic cell model in single cell and string of 100 cells., Results: Regardless of external potassium concentration or diastolic interval used, NS1643 decreases action potential duration and triangulation. For single cells, NS1643 increases the postrepolarization refractory time but shortens the absolute refractory period. In one dimensional simulations, NS1643 increases the vulnerable window for unidirectional block but suppresses the emergence of premature action potentials and unidirectional blocks around APD(90). During normokalemia, shifting the inactivation curve has greater impact than increasing conductance, whereas the opposite occurs during hypokalemia., Conclusion: Increased hERG conductance and the depolarizing shift of the inactivation curve both contribute to the antiarrhythmic actions of NS1643, with relative effects dependent on external K(+) concentration.

We have constructed computational models of canine ventricular cells and tissues, ultimately combining detailed tissue architecture and heterogeneous transmural electrophysiology. The heterogeneity is introduced by modifying the Hund-Rudy canine cell model in order to reproduce experimentally reported electrophysiological properties of endocardial, midmyocardial (M) and epicardial cells. These models are validated against experimental data for individual ionic current and action potential characteristics, and their rate dependencies. 1D and 3D heterogeneous virtual tissues are constructed, with detailed tissue architecture (anisotropy and orthotropy, due to fibre orientation and sheet structure) of the left ventricular wall wedge extracted from a diffusion tensor imaging data set. The models are used to study the effects of tissue heterogeneity and class III drugs on transmural propagation and tissue vulnerability to re-entry. We have determined relationships between the transmural dispersion of action potential duration (APD) and the vulnerable window in the 1D virtual ventricular wall, and demonstrated how changes in the transmural heterogeneity, and hence tissue vulnerability, can lead to generation of re-entry in the 3D ventricular wedge. Two class III drugs with opposite qualitative effects on transmural APD heterogeneity are considered: d-sotalol that increases transmural APD dispersion, and amiodarone that decreases it. Simulations with the 1D virtual ventricular wall show that under d-sotalol conditions the vulnerable window is substantially wider compared to amiodarone conditions, primarily in the epicardial region where unidirectional conduction block persists until the adjacent M cells are fully repolarised. Further simulations with the 3D ventricular wedge have shown that ectopic stimulation of the epicardial region results in generation of sustained re-entry under d-sotalol conditions, but not under amiodarone conditions or in control. Again, APD increase in M cells was identified as the major contributor to tissue vulnerability--re-entry was initiated primarily due to ectopic excitation propagating around the unidirectional conduction block in the M cell region. This suggests an electrophysiological mechanism for the anti- and proarrhythmic effects of the class III drugs: the relative safety of amiodarone in comparison to d-sotalol can be explained by relatively low transmural APD dispersion, and hence, a narrow vulnerable window and low probability of re-entry in the tissue.

Idiopathic short QT syndrome (SQTS) is a recently identified, genetically heterogeneous condition characterised by abbreviated QT intervals and an increased susceptibility to arrhythmia and sudden death. This simulation study identifies mechanisms by which cellular electrophysiological changes in the SQT2 (slow delayed rectifier, IKs, -linked) SQTS variant increases arrhythmia risk. The channel kinetics of the V307L mutation of the KCNQ1 subunit of the IKs channel were incorporated into human ventricular action potential (AP) models and into 1D and 2D transmural tissue simulations. Incorporating the V307L mutation into simulations reproduced defining features of the SQTS: abbreviation of the QT interval, and increases in T wave amplitude and Tpeak-Tend duration. In the single-cell model, the V307L mutation abbreviated ventricular cell AP duration at 90% repolarisation (APD90) and increased the maximal transmural voltage heterogeneity (deltaV) during APs; this resulted in augmented transmural heterogeneity of APD90 and of the effective refractory period (ERP). In the intact tissue model, the vulnerable window for unidirectional conduction block was also increased. In 2D tissue the V307L mutation facilitated and maintained reentrant excitation. Thus, in SQT2 increases in transmural heterogeneity of APD, deltaV, ERP and an increased vulnerable window for unidirectional conduction block generate an electrical substrate favourable to reentrant arrhythmia.

The architecture of the heart remains controversial despite extensive effort and recent advances in imaging techniques. Several opposing and non-mutually compatible models have been proposed to explain cardiac structure, and these models, although limited, have advanced the study and understanding of heart structure, function and development. We describe key areas of similarity and difference, highlight areas of contention and point to the important limitations of these models. Recent research in animal models on the nature, geometry and interaction of cardiac sheet structure allows unification of some seemingly conflicting features of the structural models. Intriguingly, evidence points to significant inter-individual structural variability (within constrained limits) in the canine, leading to the concept of a continuum (or distribution) of cardiac structures. This variability in heart structure partly explains the ongoing debate on myocardial architecture. These developments are used to construct an integrated description of cardiac structure unifying features of fibre, sheet and band architecture that provides a basis for (i) explaining cardiac electromechanics, (ii) computational simulations of cardiac physiology and (iii) designing interventions.

Cardiac ventricular cells and tissues are normally excitable, and are activated by propagating waves of excitation that are initiated in the specialized pacemaking region of the heart. However, isolated or repetitive activity can be initiated at abnormal (ectopic) sites in the ventricles. To trigger an endogenous ectopic beat, there must be a compact focus of cells with changed membrane excitation parameters and kinetics, which initiate activity by after-depolarizations triggered by propagating activity, or that have bifurcated into autorhythmicity. This ectopic focus needs to be large enough, and adequately coupled, to drive the surrounding tissue. We investigate the initiation of ectopic excitation in computational models of human ventricular cells triggered by after-depolarizations and by up/down regulation of specific membrane conductance systems, and the propagation and evolution of ectopic activity in homogeneous or heterogeneous and isotropic, anisotropic, or orthotropic tissues.

Voltage-dependent sodium (Na(+)) channels are heterogeneously distributed through the pacemaker of the heart, the sinoatrial node (SA node). The measured sodium channel current (i(Na)) density is higher in the periphery but low or zero in the center of the SA node. The functional roles of i(Na) in initiation and conduction of cardiac pacemaker activity remain uncertain. We evaluated the functional roles of i(Na) by computer modeling. A gradient model of the intact SA node and atrium of the rabbit heart was developed that incorporates both heterogeneities of the SA node electrophysiology and histological structure. Our computations show that a large i(Na) in the periphery helps the SA node to drive the atrial muscle. Removal i(Na) from the SA node slows down the pacemaking rate and increases the sinoatrial node-atrium conduction time. In some cases, reduction of the SA node i(Na) results in impairment of impulse initiation and conduction that leads to the SA node-atrium conduction exit block. Decrease in active SA node cell population has similar effects. Combined actions of reduced cell population and removal of i(Na) from the SA node have greater impacts on weakening the ability of the SA node to pace and drive the atrium.

Methods for the experimental and clinical investigation of cardiac arrhythmias are limited to inferring propagation within the myocardium, from surface measurements, or from electrodes at a few sites within the cardiac wall. Biophysically and anatomically detailed computational models of cardiac tissues offer a powerful way for studying the electrical propagation processes and arrhythmias within the virtual heart. We use virtual tissues to study and visualise the effects of patho- and physiological conditions, and pharmacological interventions on transmural propagation in the virtual ventricular walls. Class III drug actions are quantitatively explained by changes induced in the transmural dispersion of action potential duration. We illustrate the automated construction of a virtual anisotropic ventricle from Diffusion Tensor MRI for individual hearts, and use it to explore mechanisms leading to ventricular fibrillation. The virtual ventricular wall provides an effective tool for exploring, evaluating and visualising processes during the initiation and maintenance of ventricular arrhythmias.