Filament tension and phase-locked drift of meandering scroll waves

Rotating scroll waves are self-organising patterns which are found in many oscillating or excitable systems. Here we show that quasi-periodic (meandering) scroll waves, which include the rotors that organise cardiac arrhythmias, exhibit filament tension when averaged over the meander cycle. With strong filament curvature or medium thickness gradients, however, scroll wave dynamics are governed by phase-locked drift instead of filament tension. Our results are validated in computational models of cycloidal meander and a cardiac tissue model with linear core.Comment: accepted for publication in Physical Review Letters (December 2017

Noise Effects on the Complex Patterns of Abnormal Heartbeats

Patients at high risk for sudden death often exhibit complex heart rhythms in which abnormal heartbeats are interspersed with normal heartbeats. We analyze such a complex rhythm in a single patient over a 12-hour period and show that the rhythm can be described by a theoretical model consisting of two interacting oscillators with stochastic elements. By varying the magnitude of the noise, we show that for an intermediate level of noise, the model gives best agreement with key statistical features of the dynamics.Comment: 4 pages, 4 figures, RevTe

Conditions for propagation and block of excitation in an asymptotic model of atrial tissue

Detailed ionic models of cardiac cells are difficult for numerical simulations because they consist of a large number of equations and contain small parameters. The presence of small parameters, however, may be used for asymptotic reduction of the models. Earlier results have shown that the asymptotics of cardiac equations are non-standard. Here we apply such a novel asymptotic method to an ionic model of human atrial tissue in order to obtain a reduced but accurate model for the description of excitation fronts. Numerical simulations of spiral waves in atrial tissue show that wave fronts of propagating action potentials break-up and self-terminate. Our model, in particular, yields a simple analytical criterion of propagation block, which is similar in purpose but completely different in nature to the `Maxwell rule' in the FitzHugh-Nagumo type models. Our new criterion agrees with direct numerical simulations of break-up of re-entrant waves.Comment: Revised manuscript submitted to Biophysical Journal (30 pages incl. 10 figures

Remodelling of human atrial K+ currents but not ion channel expression by chronic β-blockade

Chronic β-adrenoceptor antagonist (β-blocker) treatment in patients is associated with a potentially anti-arrhythmic prolongation of the atrial action potential duration (APD), which may involve remodelling of repolarising K+ currents. The aim of this study was to investigate the effects of chronic β-blockade on transient outward, sustained and inward rectifier K+ currents (ITO, IKSUS and IK1) in human atrial myocytes and on the expression of underlying ion channel subunits. Ion currents were recorded from human right atrial isolated myocytes using the whole-cell-patch clamp technique. Tissue mRNA and protein levels were measured using real time RT-PCR and Western blotting. Chronic β-blockade was associated with a 41% reduction in ITO density: 9.3 ± 0.8 (30 myocytes, 15 patients) vs 15.7 ± 1.1 pA/pF (32, 14), p < 0.05; without affecting its voltage-, time- or rate dependence. IK1 was reduced by 34% at −120 mV (p < 0.05). Neither IKSUS, nor its increase by acute β-stimulation with isoprenaline, was affected by chronic β-blockade. Mathematical modelling suggested that the combination of ITO- and IK1-decrease could result in a 28% increase in APD90. Chronic β-blockade did not alter mRNA or protein expression of the ITO pore-forming subunit, Kv4.3, or mRNA expression of the accessory subunits KChIP2, KChAP, Kvβ1, Kvβ2 or frequenin. There was no reduction in mRNA expression of Kir2.1 or TWIK to account for the reduction in IK1. A reduction in atrial ITO and IK1 associated with chronic β-blocker treatment in patients may contribute to the associated action potential prolongation, and this cannot be explained by a reduction in expression of associated ion channel subunits

Competition and quality indicators in the health care sector: empirical evidence from the Dutch hospital sector

textabstractThere is much debate about the effect of competition in healthcare and especially the effect of competition on the quality of healthcare, although empirical evidence on this subject is mixed. The Netherlands provides an interesting case in this debate. The Dutch system could be characterized as a system involving managed competition and mandatory healthcare insurance. Information about the quality of care provided by hospitals has been publicly available since 2008. In this paper, we evaluate the relationship between quality scores for three diagnosis groups and the market power indicators of hospitals. We estimate the impact of competition on quality in an environment of liberalized pricing. For this research, we used unique price and production data relating to three diagnosis groups (cataract, adenoid and tonsils, bladder tumor) produced by Dutch hospitals in the period 2008–2011. We also used the quality indicators relating to these diagnosis groups. We reveal a negative relationship between market share and quality score for two of the three diagnosis groups studied, meaning that hospitals in competitive markets have better quality scores than those in concentrated markets. We therefore conclude that more competition is associated with higher quality scores

Effects of Electrical and Structural Remodeling on Atrial Fibrillation Maintenance: A Simulation Study

Atrial fibrillation, a common cardiac arrhythmia, often progresses unfavourably: in patients with long-term atrial fibrillation, fibrillatory episodes are typically of increased duration and frequency of occurrence relative to healthy controls. This is due to electrical, structural, and contractile remodeling processes. We investigated mechanisms of how electrical and structural remodeling contribute to perpetuation of simulated atrial fibrillation, using a mathematical model of the human atrial action potential incorporated into an anatomically realistic three-dimensional structural model of the human atria. Electrical and structural remodeling both shortened the atrial wavelength - electrical remodeling primarily through a decrease in action potential duration, while structural remodeling primarily slowed conduction. The decrease in wavelength correlates with an increase in the average duration of atrial fibrillation/flutter episodes. The dependence of reentry duration on wavelength was the same for electrical vs. structural remodeling. However, the dynamics during atrial reentry varied between electrical, structural, and combined electrical and structural remodeling in several ways, including: (i) with structural remodeling there were more occurrences of fragmented wavefronts and hence more filaments than during electrical remodeling; (ii) dominant waves anchored around different anatomical obstacles in electrical vs. structural remodeling; (iii) dominant waves were often not anchored in combined electrical and structural remodeling. We conclude that, in simulated atrial fibrillation, the wavelength dependence of reentry duration is similar for electrical and structural remodeling, despite major differences in overall dynamics, including maximal number of filaments, wave fragmentation, restitution properties, and whether dominant waves are anchored to anatomical obstacles or spiralling freely

Impact of Sarcoplasmic Reticulum Calcium Release on Calcium Dynamics and Action Potential Morphology in Human Atrial Myocytes: A Computational Study

Electrophysiological studies of the human heart face the fundamental challenge that experimental data can be acquired only from patients with underlying heart disease. Regarding human atria, there exist sizable gaps in the understanding of the functional role of cellular Ca2+ dynamics, which differ crucially from that of ventricular cells, in the modulation of excitation-contraction coupling. Accordingly, the objective of this study was to develop a mathematical model of the human atrial myocyte that, in addition to the sarcolemmal (SL) ion currents, accounts for the heterogeneity of intracellular Ca2+ dynamics emerging from a structurally detailed sarcoplasmic reticulum (SR). Based on the simulation results, our model convincingly reproduces the principal characteristics of Ca2+ dynamics: 1) the biphasic increment during the upstroke of the Ca2+ transient resulting from the delay between the peripheral and central SR Ca2+ release, and 2) the relative contribution of SL Ca2+ current and SR Ca2+ release to the Ca2+ transient. In line with experimental findings, the model also replicates the strong impact of intracellular Ca2+ dynamics on the shape of the action potential. The simulation results suggest that the peripheral SR Ca2+ release sites define the interface between Ca2+ and AP, whereas the central release sites are important for the fire-diffuse-fire propagation of Ca2+ diffusion. Furthermore, our analysis predicts that the modulation of the action potential duration due to increasing heart rate is largely mediated by changes in the intracellular Na+ concentration. Finally, the results indicate that the SR Ca2+ release is a strong modulator of AP duration and, consequently, myocyte refractoriness/excitability. We conclude that the developed model is robust and reproduces many fundamental aspects of the tight coupling between SL ion currents and intracellular Ca2+ signaling. Thus, the model provides a useful framework for future studies of excitation-contraction coupling in human atrial myocytes

Switches, Excitable Responses and Oscillations in the Ring1B/Bmi1 Ubiquitination System

In an active, self-ubiquitinated state, the Ring1B ligase monoubiquitinates histone H2A playing a critical role in Polycomb-mediated gene silencing. Following ubiquitination by external ligases, Ring1B is targeted for proteosomal degradation. Using biochemical data and computational modeling, we show that the Ring1B ligase can exhibit abrupt switches, overshoot transitions and self-perpetuating oscillations between its distinct ubiquitination and activity states. These different Ring1B states display canonical or multiply branched, atypical polyubiquitin chains and involve association with the Polycomb-group protein Bmi1. Bistable switches and oscillations may lead to all-or-none histone H2A monoubiquitination rates and result in discrete periods of gene (in)activity. Switches, overshoots and oscillations in Ring1B catalytic activity and proteosomal degradation are controlled by the abundances of Bmi1 and Ring1B, and the activities and abundances of external ligases and deubiquitinases, such as E6-AP and USP7