Influence of cardiac tissue anisotropy on re-entrant activation in computational models of ventricular fibrillation

The aim of this study was to establish the role played by anisotropic diffusion in (i) the number of filaments and epicardial phase singularities that sustain ventricular fibrillation in the heart, (ii) the lifetimes of filaments and phase singularities, and (iii) the creation and annihilation dynamics of filaments and phase singularities. A simplified monodomain model of cardiac tissue was used, with membrane excitation described by a simplified 3-variable model. The model was configured so that a single re-entrant wave was unstable, and fragmented into multiple re-entrant waves. Re-entry was then initiated in tissue slabs with varying anisotropy ratio. The main findings of this computational study are: (i) anisotropy ratio influenced the number of filaments Sustaining simulated ventricular fibrillation, with more filaments present in simulations with smaller values of transverse diffusion coefficient, (ii) each re-entrant filament was associated with around 0.9 phase singularities on the surface of the slab geometry, (iii) phase singularities were longer lived than filaments, and (iv) the creation and annihilation of filaments and phase singularities were linear functions of the number of filaments and phase singularities, and these relationships were independent of the anisotropy ratio. This study underscores the important role played by tissue anisotropy in cardiac ventricular fibrillation

Mechanisms of exercise-induced improvements in the contractile apparatus of the mammalian myocardium

One of the main outcomes of aerobic endurance exercise training is the improved maximal oxygen uptake, and this is pivotal to the improved work capacity that follows the exercise training. Improved maximal oxygen uptake in turn is at least partly achieved because exercise training increases the ability of the myocardium to produce a greater cardiac output. In healthy subjects, this has been demonstrated repeatedly over many decades. It has recently emerged that this scenario may also be true under conditions of an initial myocardial dysfunction. For instance, myocardial improvements may still be observed after exercise training in post-myocardial infarction heart failure. In both health and disease, it is the changes that occur in the individual cardiomyocytes with respect to their ability to contract that by and large drive the exercise training-induced adaptation to the heart. Here, we review the evidence and the mechanisms by which exercise training induces beneficial changes in the mammalian myocardium, as obtained by means of experimental and clinical studies, and argue that these changes ultimately alter the function of the whole heart and contribute to the changes in whole-body function

Cardiac cell modelling: Observations from the heart of the cardiac physiome project

In this manuscript we review the state of cardiac cell modelling in the context of international initiatives such as the IUPS Physiome and Virtual Physiological Human Projects, which aim to integrate computational models across scales and physics. In particular we focus on the relationship between experimental data and model parameterisation across a range of model types and cellular physiological systems. Finally, in the context of parameter identification and model reuse within the Cardiac Physiome, we suggest some future priority areas for this field

Multiple mechanisms of spiral wave breakup in a model of cardiac electrical activity

It has become widely accepted that the most dangerous cardiac arrhythmias are due to re- entrant waves, i.e., electrical wave(s) that re-circulate repeatedly throughout the tissue at a higher frequency than the waves produced by the heart's natural pacemaker (sinoatrial node). However, the complicated structure of cardiac tissue, as well as the complex ionic currents in the cell, has made it extremely difficult to pinpoint the detailed mechanisms of these life-threatening reentrant arrhythmias. A simplified ionic model of the cardiac action potential (AP), which can be fitted to a wide variety of experimentally and numerically obtained mesoscopic characteristics of cardiac tissue such as AP shape and restitution of AP duration and conduction velocity, is used to explain many different mechanisms of spiral wave breakup which in principle can occur in cardiac tissue. Some, but not all, of these mechanisms have been observed before using other models; therefore, the purpose of this paper is to demonstrate them using just one framework model and to explain the different parameter regimes or physiological properties necessary for each mechanism (such as high or low excitability, corresponding to normal or ischemic tissue, spiral tip trajectory types, and tissue structures such as rotational anisotropy and periodic boundary conditions). Each mechanism is compared with data from other ionic models or experiments to illustrate that they are not model-specific phenomena. The fact that many different breakup mechanisms exist has important implications for antiarrhythmic drug design and for comparisons of fibrillation experiments using different species, electromechanical uncoupling drugs, and initiation protocols.Comment: 128 pages, 42 figures (29 color, 13 b&w

A multiscale model for collagen alignment in wound healing

It is thought that collagen alignment plays a significant part in scar tissue formation during dermal wound healing. We present a multiscale model for collagen deposition and alignment during this process. We consider fibroblasts as discrete units moving within an extracellular matrix of collagen and fibrin modelled as continua. Our model includes flux induced alignment of collagen by fibroblasts, and contact guidance of fibroblasts by collagen fibres. We can use the model to predict the effects of certain manipulations, such as varying fibroblast speed, or placing an aligned piece of tissue in the wound. We also simulate experiments which alter the TGF-β concentrations in a healing dermal wound and use the model to offer an explanation of the observed influence of this growth factor on scarring

Distinct properties of the atrial voltage-gated sodium channel

Sodium handling properties of left atria (LA) were compared to those of left ventricle (LV), with the aim to understand the atria’s susceptibility to arrhythmia and improve therapy. Mouse LA sodium channels displayed distinctive activation, inactivation and recovery kinetics compared to LV sodium channels. Distinctive voltage dependence of LA sodium channel inactivation was instrumental in reducing INa in LA compared to LV, when initiated from physiological holding potentials. Flecainide sodium channel inhibition was greater in LA than LV, likely also due to differences in kinetic properties of the sodium channels between chambers. Additionally, the greater inhibitory effect of flecainide at more positive membrane potentials could result in even greater LA sodium channel inhibition in vivo. Activation and inactivation distinctions observed between LA and LV sodium channels were conserved between chambers in the Plako+/- mouse. However, there was no difference in physiological INa density, sodium channel recovery or flecainide inhibition between Plako+/- LA and LV chambers. The novel Langendorff-free isolation method produced high yields of viable mouse cardiomyocytes comparable in morphology, signalling, calcium handling and sodium channel electrophysiology to cardiomyocytes isolated using the traditional Langendorff method. This maintained that injection isolation is a valuable method for obtaining cardiomyocytes for cardiac research

Cancer modelling: Getting to the heart of the problem

Paradoxically, improvements in healthcare that have enhanced the life expectancy of humans in the Western world have, indirectly, increased the prevalence of certain types of cancer such as prostate and breast. It remains unclear whether this phenomenon should be attributed to the ageing process itself or the cumulative effect of prolonged exposure to harmful environmental stimuli such as ultraviolet light, radiation and carcinogens (Franks and Teich, 1988). Equally, there is also compelling evidence that certain genetic abnormalities can predispose individuals to specific cancers (Ilyas et al., 1999). The variety of factors that have been implicated in the development of solid tumours stems, to a large extent, from the fact that ‘cancer’ is a generic term, often used to characterize a series of disorders that share common features. At this generic level of description, cancer may be viewed as a cellular disease in which controls that usually regulate growth and maintain homeostasis are disrupted. Cancer is typically initiated by genetic mutations that lead to enhanced mitosis of a cell lineage and the formation of an avascular tumour. Since it receives nutrients by diffusion from the surrounding tissue, the size of an avascular tumour is limited to several millimeters in diameter. Further growth relies on the tumour acquiring the ability to stimulate the ingrowth of a new, circulating blood supply from the host vasculature via a process termed angiogenesis (Folkman, 1974). Once vascularised, the tumour has access to a vast nutrient source and rapid growth ensues. Further, tumour fragments that break away from the primary tumour, on entering the vasculature, may be transported to other organs in which they may establish secondary tumours or metastases that further compromise the host. Invasion is another key feature of solid tumours whereby contact with the tissue stimulates the production of enzymes that digest the tissue, liberating space into which the tumour cells migrate. Thus, cancer is a complex, multiscale process. The spatial scales of interest range from the subcellular level, to the cellular and macroscopic (or tissue) levels while the timescales may vary from seconds (or less) for signal transduction pathways to months for tumour doubling times The variety of phenomena involved, the range of spatial and temporal scales over which they act and the complex way in which they are inter-related mean that the development of realistic theoretical models of solid tumour growth is extremely challenging. While there is now a large literature focused on modelling solid tumour growth (for a review, see, for example, Preziosi, 2003), existing models typically focus on a single spatial scale and, as a result, are unable to address the fundamental problem of how phenomena at different scales are coupled or to combine, in a systematic manner, data from the various scales. In this article, a theoretical framework will be presented that is capable of integrating a hierarchy of processes occurring at different scales into a detailed model of solid tumour growth (Alarcon et al., 2004). The model is formulated as a hybrid cellular automaton and contains interlinked elements that describe processes at each spatial scale: progress through the cell cycle and the production of proteins that stimulate angiogenesis are accounted for at the subcellular level; cell-cell interactions are treated at the cellular level; and, at the tissue scale, attention focuses on the vascular network whose structure adapts in response to blood flow and angiogenic factors produced at the subcellular level. Further coupling between the different spatial scales arises from the transport of blood-borne oxygen into the tissue and its uptake at the cellular level. Model simulations will be presented to illustrate the effect that spatial heterogeneity induced by blood flow through the vascular network has on the tumour’s growth dynamics and explain how the model may be used to compare the efficacy of different anti-cancer treatment protocols