Na/K pump regulation of cardiac repolarization: insights from a systems biology approach

Repurposing Metformin for the Treatment of Atrial Fibrillation: Current Insights

Metformin is an orally effective anti-hyperglycemic drug that despite being introduced over 60 years ago is still utilized by an estimated 120 to 150 million people worldwide for the treatment of type 2 diabetes (T2D). Metformin is used off-label for the treatment of polycystic ovary syndrome (PCOS) and for pre-diabetes and weight loss. Metformin is a safe, inexpensive drug with side effects mostly limited to gastrointestinal issues. Prospective clinical data from the United Kingdom Prospective Diabetes Study (UKPDS), completed in 1998, demonstrated that metformin not only has excellent therapeutic efficacy as an anti-diabetes drug but also that good glycemic control reduced the risk of micro- and macro-vascular complications, especially in obese patients and thereby reduced the risk of diabetes-associated cardiovascular disease (CVD). Based on a long history of clinical use and an excellent safety record metformin has been investigated to be repurposed for numerous other diseases including as an anti-aging agent, Alzheimer's disease and other dementias, cancer, COVID-19 and also atrial fibrillation (AF). AF is the most frequently diagnosed cardiac arrythmia and its prevalence is increasing globally as the population ages. The argument for repurposing metformin for AF is based on a combination of retrospective clinical data and in vivo and in vitro pre-clinical laboratory studies. In this review, we critically evaluate the evidence that metformin has cardioprotective actions and assess whether the clinical and pre-clinical evidence support the use of metformin to reduce the risk and treat AF.

Bioelectric stimulation controls tissue shape and size

Epithelial tissues sheath organs and electro-mechanically regulate ion and water transport to regulate development, homeostasis, and hydrostatic organ pressure. Here, we demonstrate how external electrical stimulation allows us to control these processes in living tissues. Specifically, we electrically stimulate hollow, 3D kidneyoids and gut organoids and find that physiological-strength electrical stimulation of ∼ 5 - 10 V/cm powerfully inflates hollow tissues; a process we call electro-inflation. Electro-inflation is mediated by increased ion flux through ion channels/transporters and triggers subsequent osmotic water flow into the lumen, generating hydrostatic pressure that competes against cytoskeletal tension. Our computational studies suggest that electro-inflation is strongly driven by field-induced ion crowding on the outer surface of the tissue. Electrically stimulated tissues also break symmetry in 3D resulting from electrotaxis and affecting tissue shape. The ability of electrical cues to regulate tissue size and shape emphasizes the role and importance of the electrical micro-environment for living tissues.

Quantitative and qualitative comparison of Rubidium-82 and Oxygen-15 water cardiac PET

BACKGROUND

Differences in tracer characteristics may influence the interpretation of positron emission tomography myocardial perfusion imaging (MPI). We compare the reading of MPIs with a low-extraction retention tracer (82Rb) and a high-extraction non-retention tracer (15O-water) in a selected cohort of patients with known coronary artery disease (CAD).

METHODS

Thirty-nine patients with known CAD referred to 82Rb MPI due to angina underwent rest and stress imaging with both tracers and experienced MPI readers provided blinded consensus reads of all studies. In addition, a comparison of regional and global quantitative measures of perfusion was performed.

RESULTS

The results showed 74 % agreement in the reading of 82Rb and 15O-water MPI for regional reversible ischemia and global disease, and 82 % agreement for regional irreversible ischemia. The 15O-water MPI identified more cases of global disease (n = 12 (15O-water) vs n = 4 (82Rb), p = 0.03), whereas differences in reversible ischemia (n = 22 vs n = 16, p = 0.11) and, irreversible ischemia (n = 8 vs n = 11, p = 0.45) were not significant. The correlation between myocardial blood flow measured using the two tracers was similar to previous studies (R2 = 0.78) with wide limits of agreement (-0.93 to 0.84 ml/g/min).

CONCLUSIONS

Agreement between consensus readings of 82Rb and 15O-water MPI was good in patients with known CAD. In this limited size study, no significant differences in the identification of reversible and irreversible ischemia found, whereas 15O-water MPI had a higher positive rate for suspected global disease.

Computational Modeling of Cardiac Electrophysiology

Mathematical modeling and simulation are well-established and powerful tools to integrate experimental data of individual components of cardiac electrophysiology, excitation-contraction coupling, and regulatory signaling pathways, to gain quantitative and mechanistic insight into pathophysiological processes and guide therapeutic strategies. Here, we briefly describe the processes governing cardiac myocyte electrophysiology and Ca2+ handling and their regulation, as well as action potential propagation in tissue. We discuss the models and methods used to describe these phenomena, including procedures for model parameterization and validation, in addition to protocols for model interrogation and analysis and techniques that account for phenotypic variability and parameter uncertainty. Our objective is to provide a summary of basic concepts and approaches as a resource for scientists training in this discipline and for all researchers aiming to gain an understanding of cardiac modeling studies.

How Merkel cells transduce mechanical stimuli: A biophysical model of Merkel cells

Merkel cells combine with Aβ afferents, producing slowly adapting type 1(SA1) responses to mechanical stimuli. However, how Merkel cells transduce mechanical stimuli into neural signals to Aβ afferents is still unclear. Here we develop a biophysical model of Merkel cells for mechanical transduction by incorporating main ingredients such as Ca2+ and K+ voltage-gated channels, Piezo2 channels, internal Ca2+ stores, neurotransmitters release, and cell deformation. We first validate our model with several experiments. Then we reveal that Ca2+ and K+ channels on the plasma membrane shape the depolarization of membrane potentials, further regulating the Ca2+ transients in the cells. We also show that Ca2+ channels on the plasma membrane mainly inspire the Ca2+ transients, while internal Ca2+ stores mainly maintain the Ca2+ transients. Moreover, we show that though Piezo2 channels are rapidly adapting mechanical-sensitive channels, they are sufficient to inspire sustained Ca2+ transients in Merkel cells, which further induce the release of neurotransmitters for tens of seconds. Thus our work provides a model that captures the membrane potentials and Ca2+ transients features of Merkel cells and partly explains how Merkel cells transduce the mechanical stimuli by Piezo2 channels.

What determines the optimal pharmacological treatment of atrial fibrillation? Insights from in silico trials in 800 virtual atria

The best pharmacological treatment for each atrial fibrillation (AF) patient is unclear. We aim to exploit AF simulations in 800 virtual atria to identify key patient characteristics that guide the optimal selection of anti-arrhythmic drugs. The virtual cohort considered variability in electrophysiology and low voltage areas (LVA) and was developed and validated against experimental and clinical data from ionic currents to ECG. AF sustained in 494 (62%) atria, with large inward rectifier K+ current (IK1 ) and Na+ /K+ pump (INaK ) densities (IK1 0.11 ± 0.03 vs. 0.07 ± 0.03 S mF-1 ; INaK 0.68 ± 0.15 vs. 0.38 ± 26 S mF-1 ; sustained vs. un-sustained AF). In severely remodelled left atrium, with LVA extensions of more than 40% in the posterior wall, higher IK1 (median density 0.12 ± 0.02 S mF-1 ) was required for AF maintenance, and rotors localized in healthy right atrium. For lower LVA extensions, rotors could also anchor to LVA, in atria presenting short refractoriness (median L-type Ca2+ current, ICaL , density 0.08 ± 0.03 S mF-1 ). This atrial refractoriness, modulated by ICaL and fast Na+ current (INa ), determined pharmacological treatment success for both small and large LVA. Vernakalant was effective in atria presenting long refractoriness (median ICaL density 0.13 ± 0.05 S mF-1 ). For short refractoriness, atria with high INa (median density 8.92 ± 2.59 S mF-1 ) responded more favourably to amiodarone than flecainide, and the opposite was found in atria with low INa (median density 5.33 ± 1.41 S mF-1 ). In silico drug trials in 800 human atria identify inward currents as critical for optimal stratification of AF patient to pharmacological treatment and, together with the left atrial LVA extension, for accurately phenotyping AF dynamics. KEY POINTS: Atrial fibrillation (AF) maintenance is facilitated by small L-type Ca2+ current (ICaL ) and large inward rectifier K+ current (IK1 ) and Na+ /K+ pump. In severely remodelled left atrium, with low voltage areas (LVA) covering more than 40% of the posterior wall, sustained AF requires higher IK1 and rotors localize in healthy right atrium. For lower LVA extensions, rotors can also anchor to LVA, if the atria present short refractoriness (low ICaL ) Vernakalant is effective in atria presenting long refractoriness (high ICaL ). For short refractoriness, atria with fast Na+ current (INa ) up-regulation respond more favourably to amiodarone than flecainide, and the opposite is found in atria with low INa . The inward currents (ICaL and INa ) are critical for optimal stratification of AF patient to pharmacological treatment and, together with the left atrial LVA extension, for accurately phenotyping AF dynamics.

The role of β-adrenergic stimulation in QT interval adaptation to heart rate during stress test

The adaptation lag of the QT interval after heart rate (HR) has been proposed as an arrhythmic risk marker. Most studies have quantified the QT adaptation lag in response to abrupt, step-like changes in HR induced by atrial pacing, in response to tilt test or during ambulatory recordings. Recent studies have introduced novel methods to quantify the QT adaptation lag to gradual, ramp-like HR changes in stress tests by evaluating the differences between the measured QT series and an estimated, memoryless QT series obtained from the instantaneous HR. These studies have observed the QT adaptation lag to progressively reduce when approaching the stress peak, with the underlying mechanisms being still unclear. This study analyzes the contribution of β-adrenergic stimulation to QT interval rate adaptation in response to gradual, ramp-like HR changes. We first quantify the QT adaptation lag in Coronary Artery Disease (CAD) patients undergoing stress test. To uncover the involved mechanisms, we use biophysically detailed computational models coupling descriptions of human ventricular electrophysiology and β-adrenergic signaling, from which we simulate ventricular action potentials and ECG signals. We characterize the adaptation of the simulated QT interval in response to the HR time series measured from each of the analyzed CAD patients. We show that, when the simulated ventricular tissue is subjected to a time-varying β-adrenergic stimulation pattern, with higher stimulation levels close to the stress peak, the simulated QT interval presents adaptation lags during exercise that are more similar to those measured from the patients than when subjected to constant β-adrenergic stimulation. During stress test recovery, constant and time-varying β-adrenergic stimulation patterns render similar adaptation lags, which are generally shorter than during exercise, in agreement with results from the patients. In conclusion, our findings support the role of time-varying β-adrenergic stimulation in contributing to QT interval adaptation to gradually increasing HR changes as those seen during the exercise phase of a stress test.

Polarized NHE1 and SWELL1 regulate migration direction, efficiency and metastasis

Cell migration regulates diverse (patho)physiological processes, including cancer metastasis. According to the Osmotic Engine Model, polarization of NHE1 at the leading edge of confined cells facilitates water uptake, cell protrusion and motility. The physiological relevance of the Osmotic Engine Model and the identity of molecules mediating cell rear shrinkage remain elusive. Here, we demonstrate that NHE1 and SWELL1 preferentially polarize at the cell leading and trailing edges, respectively, mediate cell volume regulation, cell dissemination from spheroids and confined migration. SWELL1 polarization confers migration direction and efficiency, as predicted mathematically and determined experimentally via optogenetic spatiotemporal regulation. Optogenetic RhoA activation at the cell front triggers SWELL1 re-distribution and migration direction reversal in SWELL1-expressing, but not SWELL1-knockdown, cells. Efficient cell reversal also requires Cdc42, which controls NHE1 repolarization. Dual NHE1/SWELL1 knockdown inhibits breast cancer cell extravasation and metastasis in vivo, thereby illustrating the physiological significance of the Osmotic Engine Model.

Effect of Cell Cycle on Cell Surface Expression of Voltage-Gated Sodium Channels and Na+,K+-ATPase

Voltage-gated sodium channels (VGSCs) are the target for many therapies. Variation in membrane potential occurs throughout the cell cycle, yet little attention has been devoted to the role of VGSCs and Na⁺,K⁺-ATPases. We hypothesized that in addition to doubling DNA and cell membrane in anticipation of cell division, there should be a doubling of VGSCs and Na⁺,K⁺-ATPase compared to non-dividing cells. We tested this hypothesis in eight immortalized cell lines by correlating immunocytofluorescent labeling of VGSCs or Na⁺,K⁺-ATPase with propidium iodide or DAPI fluorescence using flow cytometry and imaging. Cell surface expression of VGSCs during phases S through M was double that seen during phases G0–G1. By contrast, Na⁺,K⁺-ATPase expression increased only 1.5-fold. The increases were independent of baseline expression of channels or pumps. The variation in VGSC and Na⁺,K⁺-ATPase expression has implications for both our understanding of sodium’s role in controlling the cell cycle and variability of treatments targeted at these components of the Na⁺ handling system.

Syringic Acid Ameliorates Cardiac, Hepatic, Renal and Neuronal Damage Induced by Chronic Hyperglycaemia in Wistar Rats: A Behavioural, Biochemical and Histological Analysis

This study investigated the effects of syringic acid (SA) on renal, cardiac, hepatic, and neuronal diabetic complications in streptozotocin-induced neonatal (nSTZ) diabetic rats. STZ (110 mg/kg i.p) was injected into Wistar rat neonates as a split dose (second and third postnatal day). Diabetes mellitus was diagnosed in adults by measuring fasting blood glucose levels, urine volume, and food and water intake. The treatment of SA (25 mg/kg, 50 mg/kg p.o) was given from the 8th to 18th postnatal week. To assess the development of diabetic complications and the effect of therapy, biochemical indicators in serum and behavioural parameters were recorded at specific intervals during the study period. SA (25 mg/kg, 50 mg/kg p.o) treatment reduced hyperglycaemia, polydipsia, polyphagia, polyuria, relative organ weight, cardiac hypertrophic indices, inflammatory markers, cell injury markers, glycated haemoglobin, histopathological score, and oxidative stress, and increased Na/K ATPase activity. These findings suggest that SA might significantly alleviate diabetic complications and/or renal, neuronal, cardiac, and hepatic damage in nSTZ diabetic rats.

Allopurinol and valproic acid improve cardiac triglyceride and Na+-K+-ATPase activity independent of circulating aldosterone in female rats with glucose intolerance

Context: Studies have shown that cardiac triglyceride accumulation and impaired Na⁺-K⁺-ATPase activity are linked to diabetes- related cardiovascular disease, particularly in women.Objectives: We hypothesised that allopurinol (ALL) and valproic acid (VPA) treatment would improve cardiac triglyceride and Na⁺-K⁺-ATPase activity independent of circulating aldosterone in Combined Oral Contraceptive (COC)-induced dysglycemiaMaterials and methods: Rats received COC (1.0 μg ethinylestradiol and 5.0 μg levonorgestrel; po) with or without ALL (1 mg; po) and VPA (20 mg; po) for 6 weeks.Results: COC-treatment led to impaired glucose tolerance, accumulated abdominal fat, dyslipidemia, elevated plasma MDA, PAI-1 and aldosterone levels and also reduced plasma nitric oxide bioavailability and cardiac Na⁺-K⁺-ATPase activity. However, either ALL or VPA treatment ameliorated these alterations comparably independent of elevated aldosterone levelDiscussion and conclusion: Our results suggest that either ALL or VPA would improve cardiac TG and Na⁺-K⁺-ATPase activity comparably in COC-treated rats, regardless of circulating aldosterone level.

In-silico drug trials for precision medicine in atrial fibrillation: From ionic mechanisms to electrocardiogram-based predictions in structurally-healthy human atria

Atrial fibrillation (AF) inducibility, sustainability and response to pharmacological treatment of individual patients are expected to be determined by their ionic current properties, especially in structurally-healthy atria. Mechanisms underlying AF and optimal cardioversion are however still unclear. In this study, in-silico drug trials were conducted using a population of human structurally-healthy atria models to 1) identify key ionic current properties determining AF inducibility, maintenance and pharmacological cardioversion, and 2) compare the prognostic value for predicting individual AF cardioversion of ionic current properties and electrocardiogram (ECG) metrics. In the population of structurally-healthy atria, 477 AF episodes were induced in ionic current profiles with both steep action potential duration (APD) restitution (eliciting APD alternans), and high excitability (enabling propagation at fast rates that transformed alternans into discordant). High excitability also favored 211 sustained AF episodes, so its decrease, through prolonged refractoriness, explained pharmacological cardioversion. In-silico trials over 200 AF episodes, 100 ionic profiles and 10 antiarrhythmic compounds were consistent with previous clinical trials, and identified optimal treatments for individual electrophysiological properties of the atria. Algorithms trained on 211 simulated AF episodes exhibited >70% accuracy in predictions of cardioversion for individual treatments using either ionic current profiles or ECG metrics. In structurally-healthy atria, AF inducibility and sustainability are enabled by discordant alternans, under high excitability and steep restitution conditions. Successful pharmacological cardioversion is predicted with 70% accuracy from either ionic or ECG properties, and it is optimal for treatments maximizing refractoriness (thus reducing excitability) for the given ionic current profile of the atria.

Interpretable machine learning of action potential duration restitution kinetics in single-cell models of atrial cardiomyocytes

Action potential duration (APD) restitution curve and its maximal slope (Smax) reflect single cell-level dynamic instability for inducing chaotic heart rhythms. However, conventional parameter sensitivity analysis often fails to describe nonlinear relationships between ion channel parameters and electrophysiological phenotypes, such as Smax. We explored the parameter-phenotype mapping in a population of 5000 single-cell atrial cell models through interpretable machine learning (ML) approaches. Parameter sensitivity analyses could explain the linear relationships between parameters and electrophysiological phenotypes, including APD₉₀, resting membrane potential, Vmax, refractory period, and APD/calcium alternans threshold, but not for Smax. However, neural network models had better prediction performance for Smax. To interpret the ML model, we evaluated the parameter importance at the global and local levels by computing the permutation feature importance and the local interpretable model-agnostic explanations (LIME) values, respectively. Increases in I_CaL, I_NCX, and I_Kr, and decreases in I_K1, I_b,Cl, I_Kur, I_SERCA, and I_to are correlated with higher Smax values. The LIME algorithm determined that I_NaK plays a significant role in determining Smax as well as I_to and I_Kur. The atrial cardiomyocyte population was hierarchically clustered into three distinct groups based on the LIME values and the single-cell simulation confirmed that perturbations in I_NaK resulted in different behaviors of APD restitution curves in three clusters. Our combined top-down interpretable ML and bottom-up mechanistic simulation approaches uncovered the role of I_NaK in heterogeneous behaviors of Smax in the atrial cardiomyocyte population.

Advancement of Nanofibrous Mats and Common Useful Drug Delivery Applications

Electrospinning enables simple and cost-effective production of polymer nanofibers from different polymer materials. Drug delivery systems are capable of achieving maximum drug treatment benefits by significantly reducing adverse complications. Electrospun nanofibers have recently attracted considerable attention owing to their distinctive properties, including flexibility and biocompatibility. The implementation of functional constituents within nanostructure fibers blends is an effective technique for the administration of a variety of drugs in animal research, broadening the nanofiber capability and reliability. The nanofibrous mesh and its various application purposes are discussed in terms of a summary of recent research, emphasizing the ease of streaming and a large number of combinations of this approach, which could lead to a breakthrough in targeted therapy.

Evaluation of the influence of diet supplementation with conjugated linoleic acid isomers on elemental composition in the cardio-oncological nutritional programming rat' model

BACKGROUND

The 'developmental origin of health and disease' (DOHaD) hypothesis assumes that due to the action of some stimuli during fetal life the long-term physiological changes occurs and may affects the risk of various diseases. The main aim of this study was to assess impact of supplementation of maternal' and early postnatal diet with conjugated linoleic acids (CLA) isomers on selected elements content in hearts of their female offspring with chemically induced breast cancer.

METHODS

Elemental composition was determined by quadrupole mass spectrometer with inductively coupled plasma ionization (ICP-MS).

RESULTS

The effect of maternal' diet on the elements content was more pronounced than the progeny diet modifications. Significant correlations among determined elements, especially macroelements, and lipidomic parameters, in the experimental factors dependent manner were observed. It can be concluded that supplementation of maternal and progeny diets with CLA isomers may significantly influence the contents of both macro- and microelements in cardiac tissue of newborns.

CONCLUSION

Our results also indicate, that dynamic and intricate balance among various elements in body may be affected by the lipid dietary supplements also in the pathological state. Utility of cardio-oncological approach in developmental programming study was confirmed.

Effect of ivabradine in heart failure: a meta-analysis of heart failure patients with reduced versus preserved ejection fraction

In clinical trials of heart failure reduced ejection fraction (HFrEF), ivabradine seemed to be an effective heart rate lowering agent associated with lower risk of cardiovascular death. In contrast, ivabradine failed to improve cardiovascular outcomes in heart failure preserved ejection fraction (HFpEF) despite the significant effect on heart rate. This meta-analysis is the first to compare the effects of ivabradine on heart rate and mortality parameters in HFpEF versus HFrEF. We screened three databases: PubMed, Embase, and Cochrane Library. The outcomes of these studies were mortality, reduction in heart rate, and left ventricular function improvement. We compared the efficacy of ivabradine treatment in HFpEF versus HFrEF. Heart rate analysis of pooled data showed decrease in both HFrEF (-17.646 beats/min) and HFpEF (-11.434 beats/min), and a tendency to have stronger bradycardic effect in HFrEF (p = 0.094) in randomized clinical trials. Left ventricular ejection fraction analysis revealed significant improvement in HFrEF (5.936, 95% CI: [4.199-7.672], p < 0.001) when compared with placebo (p < 0.001). We found that ivabradine significantly improves left ventricular performance in HFrEF, at the same time it exerts a tendency to have improved bradycardic effect in HFrEF. These disparate effects of ivabradine and the higher prevalence of non-cardiac comorbidities in HFpEF may explain the observed beneficial effects in HFrEF and the unchanged outcomes in HFpEF patients after ivabradine treatment.

Hydrogen, Bicarbonate, and Their Associated Exchangers in Cell Volume Regulation

Cells lacking a stiff cell wall, e.g., mammalian cells, must actively regulate their volume to maintain proper cell function. On the time scale that protein production is negligible, water flow in and out of the cell determines the cell volume variation. Water flux follows hydraulic and osmotic gradients; the latter is generated by various ion channels, transporters, and pumps in the cell membrane. Compared to the widely studied roles of sodium, potassium, and chloride in cell volume regulation, the effects of proton and bicarbonate are less understood. In this work, we use mathematical models to analyze how proton and bicarbonate, combined with sodium, potassium, chloride, and buffer species, regulate cell volume upon inhibition of ion channels, transporters, and pumps. The model includes several common, widely expressed ion transporters and focuses on obtaining generic outcomes. Results show that the intracellular osmolarity remains almost constant before and after cell volume change. The steady-state cell volume does not depend on water permeability. In addition, to ensure the stability of cell volume and ion concentrations, cells need to develop redundant mechanisms to maintain homeostasis, i.e., multiple ion channels or transporters are involved in the flux of the same ion species. These results provide insights for molecular mechanisms of cell volume regulation with additional implications for water-driven cell migration.

Electrospinning for drug delivery applications: A review

Drug delivery devices are promising tools in the pharmaceutical field, as they are able to maximize the therapeutic effects of the delivered drug while minimizing the undesired side effects. In the past years, electrospun nanofibers attracted rising attention due to their unique features, like biocompatibility and broad flexibility. Incorporation of active principles in nanofibrous meshes proved to be an efficient method for in situ delivery of a wide range of drugs, expanding the possibility and applicability of those devices. In this review, the principle of electrospinning and different fields of applications are treated to give an overview of the recent literature, underlining the easy tuning and endless combination of this technique, that in the future could be the new frontier of personalized medicine.

The influence of different liquid environments on the atomic force microscopy detection of living bEnd.3 cells

Atomic force microscopy (AFM) is one of the most important tools in the field of biomedical science, and it can be used to perform the high-resolution three-dimensional imaging of samples in liquid environments to obtain their physical properties (such as surface potentials and mechanical properties). The influence of the liquid environment on the image quality of the sample and the detection results cannot be ignored. In this work, quantitative imaging (QI) mode AFM imaging and mechanical detection were performed on mouse brain microvascular endothelial (bEnd.3) cells in different liquid environments. The gray-level variance product (SMD2) function was used to evaluate the imaging quality of the cells in liquids with different physical properties, and the variations in cell mechanical properties were quantitatively analyzed. An AFM detection liquid containing less ions and organics compared with the traditional culture medium, which is beneficial for improving the imaging quality, is introduced, and it shows similar mechanical detection results within 3 h. This can greatly reduce the detection costs and could have positive significance in the field of AFM living-cell detection.

Immediate and Delayed Response of Simulated Human Atrial Myocytes to Clinically-Relevant Hypokalemia

Although plasma electrolyte levels are quickly and precisely regulated in the mammalian cardiovascular system, even small transient changes in K⁺, Na⁺, Ca²⁺, and/or Mg²⁺ can significantly alter physiological responses in the heart, blood vessels, and intrinsic (intracardiac) autonomic nervous system. We have used mathematical models of the human atrial action potential (AP) to explore the electrophysiological mechanisms that underlie changes in resting potential (V_r) and the AP following decreases in plasma K⁺, [K⁺]_o, that were selected to mimic clinical hypokalemia. Such changes may be associated with arrhythmias and are commonly encountered in patients (i) in therapy for hypertension and heart failure; (ii) undergoing renal dialysis; (iii) with any disease with acid-base imbalance; or (iv) post-operatively. Our study emphasizes clinically-relevant hypokalemic conditions, corresponding to [K⁺]_o reductions of approximately 1.5 mM from the normal value of 4 to 4.5 mM. We show how the resulting electrophysiological responses in human atrial myocytes progress within two distinct time frames: (i) Immediately after [K⁺]_o is reduced, the K⁺-sensing mechanism of the background inward rectifier current (I_K1) responds. Specifically, its highly non-linear current-voltage relationship changes significantly as judged by the voltage dependence of its region of outward current. This rapidly alters, and sometimes even depolarizes, V_r and can also markedly prolong the final repolarization phase of the AP, thus modulating excitability and refractoriness. (ii) A second much slower electrophysiological response (developing 5-10 minutes after [K⁺]_o is reduced) results from alterations in the intracellular electrolyte balance. A progressive shift in intracellular [Na⁺]_i causes a change in the outward electrogenic current generated by the Na⁺/K⁺ pump, thereby modifying V_r and AP repolarization and changing the human atrial electrophysiological substrate. In this study, these two effects were investigated quantitatively, using seven published models of the human atrial AP. This highlighted the important role of I_K1 rectification when analyzing both the mechanisms by which [K⁺]_o regulates V_r and how the AP waveform may contribute to "trigger" mechanisms within the proarrhythmic substrate. Our simulations complement and extend previous studies aimed at understanding key factors by which decreases in [K⁺]_o can produce effects that are known to promote atrial arrhythmias in human hearts.

Simple models including energy and spike constraints reproduce complex activity patterns and metabolic disruptions

In this work, we introduce new phenomenological neuronal models (eLIF and mAdExp) that account for energy supply and demand in the cell as well as the inactivation of spike generation how these interact with subthreshold and spiking dynamics. Including these constraints, the new models reproduce a broad range of biologically-relevant behaviors that are identified to be crucial in many neurological disorders, but were not captured by commonly used phenomenological models. Because of their low dimensionality eLIF and mAdExp open the possibility of future large-scale simulations for more realistic studies of brain circuits involved in neuronal disorders. The new models enable both more accurate modeling and the possibility to study energy-associated disorders over the whole time-course of disease progression instead of only comparing the initially healthy status with the final diseased state. These models, therefore, provide new theoretical and computational methods to assess the opportunities of early diagnostics and the potential of energy-centered approaches to improve therapies.

Neurons, even “at rest”, require a constant supply of energy to function. They cannot sustain high-frequency activity over long periods because of regulatory mechanisms, such as adaptation or sodium channels inactivation, and metabolic limitations. These limitations are especially severe in many neuronal disorders, where energy can become insufficient and make the neuronal response change drastically, leading to increased burstiness, network oscillations, or seizures. Capturing such behaviors and impact of energy constraints on them is an essential prerequisite to study disorders such as Parkinson’s disease and epilepsy. However, energy and spiking constraints are not present in any of the standard neuronal models used in computational neuroscience. Here we introduce models that provide a simple and scalable way to account for these features, enabling large-scale theoretical and computational studies of neurological disorders and activity patterns that could not be captured by previously used models. These models provide a way to study energy-associated disorders over the whole time-course of disease progression, and they enable a better assessment of energy-centered approaches to improve therapies.

Cardiac transmembrane ion channels and action potentials: cellular physiology and arrhythmogenic behavior

Cardiac arrhythmias are among the leading causes of mortality. They often arise from alterations in the electrophysiological properties of cardiac cells and their underlying ionic mechanisms. It is therefore critical to further unravel the pathophysiology of the ionic basis of human cardiac electrophysiology in health and disease. In the first part of this review, current knowledge on the differences in ion channel expression and properties of the ionic processes that determine the morphology and properties of cardiac action potentials and calcium dynamics from cardiomyocytes in different regions of the heart are described. Then the cellular mechanisms promoting arrhythmias in congenital or acquired conditions of ion channel function (electrical remodeling) are discussed. The focus is on human-relevant findings obtained with clinical, experimental, and computational studies, given that interspecies differences make the extrapolation from animal experiments to human clinical settings difficult. Deepening the understanding of the diverse pathophysiology of human cellular electrophysiology will help in developing novel and effective antiarrhythmic strategies for specific subpopulations and disease conditions.

Elucidation of Cellular Mechanisms That Regulate the Sustained Contraction and Relaxation of the Mammalian Iris

Purpose

In mammals, pupil constriction and dilation form the pupillary light reflex (PLR), which is mediated by both brain-regulated (parasympathetic) and local iris-driven reflexes. To better understand the cellular mechanisms that regulate pupil physiological dynamics via central and local photoreception, we have examined the regulation of the PLR via parasympathetic and local activation, respectively.

Methods

In this study, the PLR was examined in mouse enucleated eyes ex vivo in real-time under different ionic conditions in response to acetylcholine and/or blue light (480 nm). The use of pupillometry recordings captured the relaxation, contraction, and pupil escape (redilation) processes for 10 minutes up to 1 hour.

Results

Among others, our results show that ryanodine receptor channels are the main driver for iridal stimulation-contraction coupling, in which extracellular influx of Ca2+ is required for amplification of pupil constriction. Both local and parasympathetic iridal activations are necessary, but not sufficient for sustained pupil constriction. Moreover, the degree of membrane potential repolarization in the dark is correlated with the latency and velocity of iridal constriction. Furthermore, pupil escape is driven by membrane potential hyperpolarization where voltage-gated potassium channels play a crucial role.

Conclusions

Together, this study presents new mechanisms regulating synchronized pupil dilation and contraction, sustained pupil constriction, iridal stimulation-contraction coupling, and pupil escape.

Calibration of ionic and cellular cardiac electrophysiology models

Cardiac electrophysiology models are among the most mature and well-studied mathematical models of biological systems. This maturity is bringing new challenges as models are being used increasingly to make quantitative rather than qualitative predictions. As such, calibrating the parameters within ion current and action potential (AP) models to experimental data sets is a crucial step in constructing a predictive model. This review highlights some of the fundamental concepts in cardiac model calibration and is intended to be readily understood by computational and mathematical modelers working in other fields of biology. We discuss the classic and latest approaches to calibration in the electrophysiology field, at both the ion channel and cellular AP scales. We end with a discussion of the many challenges that work to date has raised and the need for reproducible descriptions of the calibration process to enable models to be recalibrated to new data sets and built upon for new studies. This article is categorized under: Analytical and Computational Methods > Computational Methods Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Cellular Models.

Mathematical model of the ventricular action potential and effects of isoproterenol-induced cardiac hypertrophy in rats

Mathematical action potential (AP) modeling is a well-established but still-developing area of research to better understand physiological and pathological processes. In particular, changes in AP mechanisms in the isoproterenol (ISO) -induced hypertrophic heart model are incompletely understood. Here we present a mathematical model of the rat AP based on recordings from rat ventricular myocytes. In our model, for the first time, all channel kinetics are defined with a single type of function that is simple and easy to apply. The model AP and channels dynamics are consistent with the APs recorded from rats for both Control (absence of ISO) and ISO-treated cases. Our mathematical model helps us to understand the reason for the prolongation in AP duration after ISO application while ISO treatment helps us to validate our mathematical model. We reveal that the smaller density and the slower gating kinetics of the transient K⁺ current help explain the prolonged AP duration after ISO treatment and the increasing amplitude of the rapid and the slow inward rectifier currents also contribute to this prolongation alongside the flux in Ca²⁺ currents. ISO induced an increase in the density of the Na⁺ current that can explain the faster upstroke. We believe that AP dynamics from rat ventricular myocytes can be reproduced very well with this mathematical model and that it provides a powerful tool for improved insights into the underlying dynamics of clinically important AP properties such as ISO application.

Expansion of the Transporter-Opsin-G protein-coupled receptor superfamily with five new protein families

Here we provide bioinformatic evidence that the Organo-Arsenical Exporter (ArsP), Endoplasmic Reticulum Retention Receptor (KDELR), Mitochondrial Pyruvate Carrier (MPC), L-Alanine Exporter (AlaE), and the Lipid-linked Sugar Translocase (LST) protein families are members of the Transporter-Opsin-G Protein-coupled Receptor (TOG) Superfamily. These families share domains homologous to well-established TOG superfamily members, and their topologies of transmembranal segments (TMSs) are compatible with the basic 4-TMS repeat unit characteristic of this Superfamily. These repeat units tend to occur twice in proteins as a result of intragenic duplication events, often with subsequent gain/loss of TMSs in many superfamily members. Transporters within the ArsP family allow microbial pathogens to expel toxic arsenic compounds from the cell. Members of the KDELR family are involved in the selective retrieval of proteins that reside in the endoplasmic reticulum. Proteins of the MPC family are involved in the transport of pyruvate into mitochondria, providing the organelle with a major oxidative fuel. Members of family AlaE excrete L-alanine from the cell. Members of the LST family are involved in the translocation of lipid-linked glucose across the membrane. These five families substantially expand the range of substrates of transport carriers in the superfamily, although KDEL receptors have no known transport function. Clustering of protein sequences reveals the relationships among families, and the resulting tree correlates well with the degrees of sequence similarity documented between families. The analyses and programs developed to detect distant relatedness, provide insights into the structural, functional, and evolutionary relationships that exist between families of the TOG superfamily, and should be of value to many other investigators.

A Heart for Diversity: Simulating Variability in Cardiac Arrhythmia Research

In cardiac electrophysiology, there exist many sources of inter- and intra-personal variability. These include variability in conditions and environment, and genotypic and molecular diversity, including differences in expression and behavior of ion channels and transporters, which lead to phenotypic diversity (e.g., variable integrated responses at the cell, tissue, and organ levels). These variabilities play an important role in progression of heart disease and arrhythmia syndromes and outcomes of therapeutic interventions. Yet, the traditional in silico framework for investigating cardiac arrhythmias is built upon a parameter/property-averaging approach that typically overlooks the physiological diversity. Inspired by work done in genetics and neuroscience, new modeling frameworks of cardiac electrophysiology have been recently developed that take advantage of modern computational capabilities and approaches, and account for the variance in the biological data they are intended to illuminate. In this review, we outline the recent advances in statistical and computational techniques that take into account physiological variability, and move beyond the traditional cardiac model-building scheme that involves averaging over samples from many individuals in the construction of a highly tuned composite model. We discuss how these advanced methods have harnessed the power of big (simulated) data to study the mechanisms of cardiac arrhythmias, with a special emphasis on atrial fibrillation, and improve the assessment of proarrhythmic risk and drug response. The challenges of using in silico approaches with variability are also addressed and future directions are proposed.

Reproducible model development in the cardiac electrophysiology Web Lab

The modelling of the electrophysiology of cardiac cells is one of the most mature areas of systems biology. This extended concentration of research effort brings with it new challenges, foremost among which is that of choosing which of these models is most suitable for addressing a particular scientific question. In a previous paper, we presented our initial work in developing an online resource for the characterisation and comparison of electrophysiological cell models in a wide range of experimental scenarios. In that work, we described how we had developed a novel protocol language that allowed us to separate the details of the mathematical model (the majority of cardiac cell models take the form of ordinary differential equations) from the experimental protocol being simulated. We developed a fully-open online repository (which we termed the Cardiac Electrophysiology Web Lab) which allows users to store and compare the results of applying the same experimental protocol to competing models. In the current paper we describe the most recent and planned extensions of this work, focused on supporting the process of model building from experimental data. We outline the necessary work to develop a machine-readable language to describe the process of inferring parameters from wet lab datasets, and illustrate our approach through a detailed example of fitting a model of the hERG channel using experimental data. We conclude by discussing the future challenges in making further progress in this domain towards our goal of facilitating a fully reproducible approach to the development of cardiac cell models.

Electromechanics and Volume Dynamics in Nonexcitable Tissue Cells

Cell volume regulation is fundamentally important in phenomena such as cell growth, proliferation, tissue homeostasis, and embryogenesis. How the cell size is set, maintained, and changed over a cell's lifetime is not well understood. In this work we focus on how the volume of nonexcitable tissue cells is coupled to the cell membrane electrical potential and the concentrations of membrane-permeable ions in the cell environment. Specifically, we demonstrate that a sudden cell depolarization using the whole-cell patch clamp results in a 50% increase in cell volume, whereas hyperpolarization results in a slight volume decrease. We find that cell volume can be partially controlled by changing the chloride or the sodium/potassium concentrations in the extracellular environment while maintaining a constant external osmotic pressure. Depletion of external chloride leads to a volume decrease in suspended HN31 cells. Introducing cells to a high-potassium solution causes volume increase up to 50%. Cell volume is also influenced by cortical tension: actin depolymerization leads to cell volume increase. We present an electrophysiology model of water dynamics driven by changes in membrane potential and the concentrations of permeable ions in the cells surrounding. The model quantitatively predicts that the cell volume is directly proportional to the intracellular protein content.

Mechanism of doxorubicin cardiotoxicity evaluated by integrating multiple molecular effects into a biophysical model

BACKGROUND AND PURPOSE

Doxorubicin (DOX) is an effective cancer therapeutic agent but causes therapy-limiting cardiotoxicity. The effects of DOX and its metabolite doxorubicinol (DOXL) on individual channels have been well characterized in isolation. However, it is unknown how the action and interaction of affected channels combine to generate the phenotypic cardiotoxic outcome. We sought to develop an in silico model that links drug effects on channels to action potential duration (APD) and intracellular Ca²⁺ concentration in order to address this gap in knowledge.

EXPERIMENTAL APPROACH

We first propose two methods to obtain, from published values, consensus drug effects on the currents of individual channels, transporters and pumps. Separately, we obtained equivalent values for APD and Ca²⁺ concentration (the readouts used as surrogates for cardiotoxicity). Once derived, the consensus effects on the currents were incorporated into established biophysical models of the cardiac myocyte and were refined adjusting the sarcoplasmic reticulum Ca²⁺ leak current (I_Leak ) until the consensus effects on APD and Ca²⁺ dynamics were replicated. Using factorial analysis, we then quantified the relative contribution of each channel to DOX and DOXL cardiotoxicity.

KEY RESULTS

The factorial analysis identified the rapid delayed rectifying K⁺ current, the L-type Ca²⁺ current and the sarcoplasmic reticulum I_Leak as the targets primarily responsible for the cardiotoxic effects on APD and Ca²⁺ dynamics.

CONCLUSIONS AND IMPLICATIONS

This study provides insight into the mechanisms of DOX-induced cardiotoxicity and a framework for the development of future diagnostic and therapeutic strategies.

Slow Adaptation of Ventricular Repolarization as a Cause of Arrhythmia?

Introduction: This article is part of the Focus Theme of Methods of Information in Medicine on “Biosignal Interpretation: Advanced Methods for Studying Cardiovascular and Respiratory Systems”.

Background: Adaptation of the QT-interval to changes in heart rate reflects on the body-surface electrocardiogram the adaptation of action potential duration (APD) at the cellular level. The initial fast phase of APD adaptation has been shown to modulate the arrhythmia substrate. Whether the slow phase is potentially proarrhythmic remains unclear.

Objectives: To analyze in-vivo human data and use computer simulations to examine effects of the slow APD adaptation phase on dispersion of repolarization and reentry in the human ventricle.

Methods: Electrograms were acquired from 10 left and 10 right ventricle (LV/RV) endocardial sites in 15 patients with normal ventricles during RV pacing. Activation-recovery intervals, as a surrogate for APD, were measured during a sustained increase in heart rate. Observed dynamics were studied using computer simulations of human tissue electrophysiology.

Results: Spatial heterogeneity of rate adaptation was observed in all patients. Inhomogeneity in slow APD adaptation time constants (ΔTs) was greater in LV than RV (ΔTs LV = 31.8 ± 13.2, ΔTs RV = 19.0 ± 12.8 s, P < 0.01). Simulations showed that altering local slow time constants of adaptation was sufficient to convert partial wavefront block to block with successful reentry.

Conclusions: Using electrophysiological data acquired in-vivo in human and computer simulations, we identify heterogeneity in the slow phase of APD adaptation as an important component of arrhythmogenesis.

Sensitivity analysis revealing the effect of modulating ionic mechanisms on calcium dynamics in simulated human heart failure

Abnormal intracellular Ca²⁺ handling is the major contributor to the depressed cardiac contractility observed in heart failure. The electrophysiological remodeling associated with this pathology alters both the action potential and the Ca²⁺ dynamics, leading to a defective excitation-contraction coupling that ends in mechanical dysfunction. The importance of maintaining a correct intracellular Ca²⁺ concentration requires a better understanding of its regulation by ionic mechanisms. To study the electrical activity and ionic homeostasis of failing myocytes, a modified version of the O’Hara et al. human action potential model was used, including electrophysiological remodeling. The impact of the main ionic transport mechanisms was analyzed using single-parameter sensitivity analyses, the first of which explored the modulation of electrophysiological characteristics related to Ca²⁺ exerted by the remodeled parameters. The second sensitivity analysis compared the potential consequences of modulating individual channel conductivities, as one of the main effects of potential drugs, on Ca²⁺ dynamic properties under both normal conditions and in heart failure. The first analysis revealed the important contribution of the sarcoplasmic reticulum Ca²⁺-ATPase (SERCA) dysfunction to the altered Ca²⁺ homeostasis, with the Na⁺/Ca²⁺ exchanger (NCX) and other Ca²⁺ cycling proteins also playing a significant role. Our results highlight the importance of improving the SR uptake function to increase Ca²⁺ content and restore Ca²⁺ homeostasis and contractility. The second sensitivity analysis highlights the different response of the failing myocyte versus the healthy myocyte to potential pharmacological actions on single channels. The result of modifying the conductances of the remodeled proteins such as SERCA and NCX in heart failure has less impact on Ca²⁺ modulation. These differences should be taken into account when designing drug therapies.

Action Potential Shortening and Impairment of Cardiac Function by Ablation of Slc26a6

BACKGROUND

Intracellular pH (pH_i) is critical to cardiac excitation and contraction; uncompensated changes in pH_i impair cardiac function and trigger arrhythmia. Several ion transporters participate in cardiac pH_i regulation. Our previous studies identified several isoforms of a solute carrier Slc26a6 to be highly expressed in cardiomyocytes. We show that Slc26a6 mediates electrogenic Cl^-/HCO₃^- exchange activities in cardiomyocytes, suggesting the potential role of Slc26a6 in regulation of not only pH_i, but also cardiac excitability.

METHODS AND RESULTS

To test the mechanistic role of Slc26a6 in the heart, we took advantage of Slc26a6 knockout (Slc26a6^-/- ) mice using both in vivo and in vitro analyses. Consistent with our prediction of its electrogenic activities, ablation of Slc26a6 results in action potential shortening. There are reduced Ca²⁺ transient and sarcoplasmic reticulum Ca²⁺ load, together with decreased sarcomere shortening in Slc26a6^-/- cardiomyocytes. These abnormalities translate into reduced fractional shortening and cardiac contractility at the in vivo level. Additionally, pH_i is elevated in Slc26a6^-/- cardiomyocytes with slower recovery kinetics from intracellular alkalization, consistent with the Cl^-/HCO₃^- exchange activities of Slc26a6. Moreover, Slc26a6^-/- mice show evidence of sinus bradycardia and fragmented QRS complex, supporting the critical role of Slc26a6 in cardiac conduction system.

CONCLUSIONS

Our study provides mechanistic insights into Slc26a6, a unique cardiac electrogenic Cl^-/HCO₃^- transporter in ventricular myocytes, linking the critical roles of Slc26a6 in regulation of pH_i, excitability, and contractility. pH_i is a critical regulator of other membrane and contractile proteins. Future studies are needed to investigate possible changes in these proteins in Slc26a6^-/- mice.

Slow [Na+]i dynamics impacts arrhythmogenesis and spiral wave reentry in cardiac myocyte ionic model

Accumulation of intracellular Na⁺ is gaining recognition as an important regulator of cardiac myocyte electrophysiology. The intracellular Na⁺ concentration can be an important determinant of the cardiac action potential duration, can modulate the tissue-level conduction of excitation waves, and can alter vulnerability to arrhythmias. Mathematical models of cardiac electrophysiology often incorporate a dynamic intracellular Na⁺ concentration, which changes much more slowly than the remaining variables. We investigated the dependence of several arrhythmogenesis-related factors on [Na⁺]_i in a mathematical model of the human atrial action potential. In cell simulations, we found that [Na⁺]_i accumulation stabilizes the action potential duration to variations in several conductances and that the slow dynamics of [Na⁺]_i impacts bifurcations to pro-arrhythmic afterdepolarizations, causing intermittency between different rhythms. In long-lasting tissue simulations of spiral wave reentry, [Na⁺]_i becomes spatially heterogeneous with a decreased area around the spiral wave rotation center. This heterogeneous region forms a functional anchor, resulting in diminished meandering of the spiral wave. Our findings suggest that slow, physiological, rate-dependent variations in [Na⁺]_i may play complex roles in cellular and tissue-level cardiac dynamics.

Atrial Fibrillation Dynamics and Ionic Block Effects in Six Heterogeneous Human 3D Virtual Atria with Distinct Repolarization Dynamics

Atrial fibrillation (AF) usually manifests as reentrant circuits propagating through the whole atria creating chaotic activation patterns. Little is yet known about how differences in electrophysiological and ionic properties between patients modulate reentrant patterns in AF. The goal of this study is to quantify how variability in action potential duration (APD) at different stages of repolarization determines AF dynamics and their modulation by ionic block using a set of virtual whole-atria human models. Six human whole-atria models are constructed based on the same anatomical structure and fiber orientation, but with different electrophysiological phenotypes. Membrane kinetics for each whole-atria model are selected with distinct APD characteristics at 20, 50, and 90% repolarization, from an experimentally calibrated population of human atrial action potential models, including AF remodeling and acetylcholine parasympathetic effects. Our simulations show that in all whole-atria models, reentrant circuits tend to organize around the pulmonary veins and the right atrial appendage, thus leading to higher dominant frequency (DF) and more organized activation in the left atrium than in the right atrium. Differences in APD in all phases of repolarization (not only APD₉₀) yielded quantitative differences in fibrillation patterns with long APD associated with slower and more regular dynamics. Long APD₅₀ and APD₂₀ were associated with increased interatrial conduction block and interatrial differences in DF and organization index, creating reentry instability and self-termination in some cases. Specific inhibitions of I_K1, I_NaK, or I_Na reduce DF and organization of the arrhythmia by enlarging wave meandering, reducing the number of secondary wavelets, and promoting interatrial block in all six virtual patients, especially for the phenotypes with short APD at 20, 50, and/or 90% repolarization. This suggests that therapies aiming at prolonging the early phase of repolarization might constitute effective antiarrhythmic strategies for the pharmacological management of AF. In summary, simulations report significant differences in atrial fibrillatory dynamics resulting from differences in APD at all phases of repolarization.

The Electrogenic Na+/K+ Pump Is a Key Determinant of Repolarization Abnormality Susceptibility in Human Ventricular Cardiomyocytes: A Population-Based Simulation Study

Background: Cellular repolarization abnormalities occur unpredictably due to disease and drug effects, and can occur even in cardiomyocytes that exhibit normal action potentials (AP) under control conditions. Variability in ion channel densities may explain differences in this susceptibility to repolarization abnormalities. Here, we quantify the importance of key ionic mechanisms determining repolarization abnormalities following ionic block in human cardiomyocytes yielding normal APs under control conditions.

Methods and Results: Sixty two AP recordings from non-diseased human heart preparations were used to construct a population of human ventricular models with normal APs and a wide range of ion channel densities. Multichannel ionic block was applied to investigate susceptibility to repolarization abnormalities. I_Kr block was necessary for the development of repolarization abnormalities. Models that developed repolarization abnormalities over the widest range of blocks possessed low Na⁺/K⁺ pump conductance below 50% of baseline, and I_CaL conductance above 70% of baseline. Furthermore, I_NaK made the second largest contribution to repolarizing current in control simulations and the largest contribution under 75% I_Kr block. Reversing intracellular Na⁺ overload caused by reduced I_NaK was not sufficient to prevent abnormalities in models with low Na⁺/K⁺ pump conductance, while returning Na⁺/K⁺ pump conductance to normal substantially reduced abnormality occurrence, indicating I_NaK is an important repolarization current.

Conclusions: I_NaK is an important determinant of repolarization abnormality susceptibility in human ventricular cardiomyocytes, through its contribution to repolarization current rather than homeostasis. While we found I_Kr block to be necessary for repolarization abnormalities to occur, I_NaK decrease, as in disease, may amplify the pro-arrhythmic risk of drug-induced I_Kr block in humans.

Computer modelling reveals new conformers of the ATP binding loop of Na+/K+-ATPase involved in the transphosphorylation process of the sodium pump

Hydrolysis of ATP by Na⁺/K⁺-ATPase, a P-Type ATPase, catalyzing active Na⁺ and K⁺ transport through cellular membranes leads transiently to a phosphorylation of its catalytical α-subunit. Surprisingly, three-dimensional molecular structure analysis of P-type ATPases reveals that binding of ATP to the N-domain connected by a hinge to the P-domain is much too far away from the Asp³⁶⁹ to allow the transfer of ATP’s terminal phosphate to its aspartyl-phosphorylation site. In order to get information for how the transfer of the γ-phosphate group of ATP to the Asp³⁶⁹ is achieved, analogous molecular modeling of the M₄–M₅ loop of ATPase was performed using the crystal data of Na⁺/K⁺-ATPase of different species. Analogous molecular modeling of the cytoplasmic loop between Thr³³⁸ and Ile⁷⁶⁰ of the α₂-subunit of Na⁺/K⁺-ATPase and the analysis of distances between the ATP binding site and phosphorylation site revealed the existence of two ATP binding sites in the open conformation; the first one close to Phe⁴⁷⁵ in the N-domain, the other one close to Asp³⁶⁹ in the P-domain. However, binding of Mg²⁺•ATP to any of these sites in the “open conformation” may not lead to phosphorylation of Asp³⁶⁹. Additional conformations of the cytoplasmic loop were found wobbling between “open conformation” <==> “semi-open conformation <==> “closed conformation” in the absence of 2Mg²⁺•ATP. The cytoplasmic loop’s conformational change to the “semi-open conformation”—characterized by a hydrogen bond between Arg⁵⁴³ and Asp⁶¹¹—triggers by binding of 2Mg²⁺•ATP to a single ATP site and conversion to the “closed conformation” the phosphorylation of Asp³⁶⁹ in the P-domain, and hence the start of Na⁺/K⁺-activated ATP hydrolysis.

Cell mechanics: a dialogue

Under the microscope, eukaryotic animal cells can adopt a variety of different shapes and sizes. These cells also move and deform, and the physical mechanisms driving these movements and shape changes are important in fundamental cell biology, tissue mechanics, as well as disease biology. This article reviews some of the basic mechanical concepts in cells, emphasizing continuum mechanics description of cytoskeletal networks and hydrodynamic flows across the cell membrane. We discuss how cells can generate movement and shape changes by controlling mass fluxes at the cell boundary. These mass fluxes can come from polymerization/depolymerization of actin cytoskeleton, as well as osmotic and hydraulic pressure-driven flow of water across the cell membrane. By combining hydraulic pressure control with force balance conditions at the cell surface, we discuss a quantitative mechanism of cell shape and volume control. The broad consequences of this model on cell mechanosensation and tissue mechanics are outlined.

Potassium currents in the heart: functional roles in repolarization, arrhythmia and therapeutics

This is the second of the two White Papers from the fourth UC Davis Cardiovascular Symposium Systems Approach to Understanding Cardiac Excitation-Contraction Coupling and Arrhythmias (3-4 March 2016), a biennial event that brings together leading experts in different fields of cardiovascular research. The theme of the 2016 symposium was 'K⁺ channels and regulation', and the objectives of the conference were severalfold: (1) to identify current knowledge gaps; (2) to understand what may go wrong in the diseased heart and why; (3) to identify possible novel therapeutic targets; and (4) to further the development of systems biology approaches to decipher the molecular mechanisms and treatment of cardiac arrhythmias. The sessions of the Symposium focusing on the functional roles of the cardiac K⁺ channel in health and disease, as well as K⁺ channels as therapeutic targets, were contributed by Ye Chen-Izu, Gideon Koren, James Weiss, David Paterson, David Christini, Dobromir Dobrev, Jordi Heijman, Thomas O'Hara, Crystal Ripplinger, Zhilin Qu, Jamie Vandenberg, Colleen Clancy, Isabelle Deschenes, Leighton Izu, Tamas Banyasz, Andras Varro, Heike Wulff, Eleonora Grandi, Michael Sanguinetti, Donald Bers, Jeanne Nerbonne and Nipavan Chiamvimonvat as speakers and panel discussants. This article summarizes state-of-the-art knowledge and controversies on the functional roles of cardiac K⁺ channels in normal and diseased heart. We endeavour to integrate current knowledge at multiple scales, from the single cell to the whole organ levels, and from both experimental and computational studies.

Balance between sodium and calcium currents underlying chronic atrial fibrillation termination: An in silico intersubject variability study

BACKGROUND

Atrial remodeling as a result of long-standing persistent atrial fibrillation (AF) induces substrate modifications that lead to different perpetuation mechanisms than in paroxysmal AF and a reduction in the efficacy of antiarrhythmic treatments.

OBJECTIVE

The purpose of this study was to identify the ionic current modifications that could destabilize reentries during chronic AF and serve to personalize antiarrhythmic strategies.

METHODS

A population of 173 mathematical models of remodeled human atrial tissue with realistic intersubject variability was developed based on action potential recordings of 149 patients diagnosed with AF. The relationship of each ionic current with AF maintenance and the dynamics of functional reentries (rotor meandering, dominant frequency) were evaluated by means of 3-dimensional simulations.

RESULTS

Self-sustained reentries were maintained in 126 (73%) of the simulations. AF perpetuation was associated with higher expressions of I_Na and I_CaL (P <.01), with no significant differences in the remaining currents. I_CaL blockade promoted AF extinction in 30% of these 126 models. The mechanism of AF termination was related with collisions between rotors because of an increase in rotor meandering (1.71 ± 2.01cm²) and presented an increased efficacy in models with a depressed I_Na (P <.01).

CONCLUSION

Mathematical simulations based on a population of models representing intersubject variability allow the identification of ionic mechanisms underlying rotor dynamics and the definition of new personalized pharmacologic strategies. Our results suggest that the underlying mechanism of the diverging success of I_CaL block as an antiarrhythmic strategy is dependent on the basal availability of sodium and calcium ion channel conductivities.

Mechanisms of pro-arrhythmic abnormalities in ventricular repolarisation and anti-arrhythmic therapies in human hypertrophic cardiomyopathy

INTRODUCTION

Hypertrophic cardiomyopathy (HCM) is a cause of sudden arrhythmic death, but the understanding of its pro-arrhythmic mechanisms and an effective pharmacological treatment are lacking. HCM electrophysiological remodelling includes both increased inward and reduced outward currents, but their role in promoting repolarisation abnormalities remains unknown. The goal of this study is to identify key ionic mechanisms driving repolarisation abnormalities in human HCM, and to evaluate anti-arrhythmic effects of single and multichannel inward current blocks.

METHODS

Experimental ionic current, action potential (AP) and Ca(2+)-transient (CaT) recordings were used to construct populations of human non-diseased and HCM AP models (n=9118), accounting for inter-subject variability. Simulations were conducted for several degrees of selective and combined inward current block.

RESULTS

Simulated HCM cardiomyocytes exhibited prolonged AP and CaT, diastolic Ca(2+) overload and decreased CaT amplitude, in agreement with experiments. Repolarisation abnormalities in HCM models were consistently driven by L-type Ca(2+) current (ICaL) re-activation, and ICaL block was the most effective intervention to normalise repolarisation and diastolic Ca(2+), but compromised CaT amplitude. Late Na(+) current (INaL) block partially abolished repolarisation abnormalities, with small impact on CaT. Na(+)/Ca(2+) exchanger (INCX) block effectively restored repolarisation and CaT amplitude, but increased Ca(2+) overload. Multichannel block increased efficacy in normalising repolarisation, AP biomarkers and CaT amplitude compared to selective block.

CONCLUSIONS

Experimentally-calibrated populations of human AP models identify ICaL re-activation as the key mechanism for repolarisation abnormalities in HCM, and combined INCX, INaL and ICaL block as effective anti-arrhythmic therapies also able to partially reverse the HCM electrophysiological phenotype.

Multi-scale Modeling of the Cardiovascular System: Disease Development, Progression, and Clinical Intervention

Cardiovascular diseases (CVDs) are the leading cause of death in the western world. With the current development of clinical diagnostics to more accurately measure the extent and specifics of CVDs, a laudable goal is a better understanding of the structure-function relation in the cardiovascular system. Much of this fundamental understanding comes from the development and study of models that integrate biology, medicine, imaging, and biomechanics. Information from these models provides guidance for developing diagnostics, and implementation of these diagnostics to the clinical setting, in turn, provides data for refining the models. In this review, we introduce multi-scale and multi-physical models for understanding disease development, progression, and designing clinical interventions. We begin with multi-scale models of cardiac electrophysiology and mechanics for diagnosis, clinical decision support, personalized and precision medicine in cardiology with examples in arrhythmia and heart failure. We then introduce computational models of vasculature mechanics and associated mechanical forces for understanding vascular disease progression, designing clinical interventions, and elucidating mechanisms that underlie diverse vascular conditions. We conclude with a discussion of barriers that must be overcome to provide enhanced insights, predictions, and decisions in pre-clinical and clinical applications.

Flow-Driven Cell Migration under External Electric Fields

Electric fields influence many aspects of cell physiology, including various forms of cell migration. Many cells are sensitive to electric fields, and they can migrate toward a cathode or an anode, depending on the cell type. In this Letter, we examine an actomyosin-independent mode of cell migration under electrical fields. Our theory considers a one-dimensional cell with water and ionic fluxes at the cell boundary. Water fluxes through the membrane are governed by the osmotic pressure difference across the cell membrane. Fluxes of cations and anions across the cell membrane are determined by the properties of the ion channels as well as the external electric field. Results show that without actin polymerization and myosin contraction, electric fields can also drive cell migration, even when the cell is not polarized. The direction of migration with respect to the electric field direction is influenced by the properties of ion channels, and are cell-type dependent.

Cell-specific cardiac electrophysiology models

The traditional cardiac model-building paradigm involves constructing a composite model using data collected from many cells. Equations are derived for each relevant cellular component (e.g., ion channel, exchanger) independently. After the equations for all components are combined to form the composite model, a subset of parameters is tuned, often arbitrarily and by hand, until the model output matches a target objective, such as an action potential. Unfortunately, such models often fail to accurately simulate behavior that is dynamically dissimilar (e.g., arrhythmia) to the simple target objective to which the model was fit. In this study, we develop a new approach in which data are collected via a series of complex electrophysiology protocols from single cardiac myocytes and then used to tune model parameters via a parallel fitting method known as a genetic algorithm (GA). The dynamical complexity of the electrophysiological data, which can only be fit by an automated method such as a GA, leads to more accurately parameterized models that can simulate rich cardiac dynamics. The feasibility of the method is first validated computationally, after which it is used to develop models of isolated guinea pig ventricular myocytes that simulate the electrophysiological dynamics significantly better than does a standard guinea pig model. In addition to improving model fidelity generally, this approach can be used to generate a cell-specific model. By so doing, the approach may be useful in applications ranging from studying the implications of cell-to-cell variability to the prediction of intersubject differences in response to pharmacological treatment.

Mathematical models of cardiac cell electrophysiology are widely used as predictive and illuminatory tools, but have been developed for decades using a suboptimal process. The models are typically constructed by manual adjustment of parameters to fit simple data and therefore often underperform when used to predict complex behavior such as arrhythmias. We present a novel method of model parameterization using automated optimization and dynamically rich fitting data and then demonstrate that this approach is better at finding the “real” model of a cell. Application of the method to cardiac myocytes leads to cell-specific models, which may enable well-controlled studies of both cellular- and subject-level population heterogeneity in disease propensity and response to therapies.

Abnormal calcium homeostasis in heart failure with preserved ejection fraction is related to both reduced contractile function and incomplete relaxation: an electromechanically detailed biophysical modeling study

Heart failure with preserved ejection fraction (HFpEF) accounts for about 50% of heart failure cases. It has features of incomplete relaxation and increased stiffness of the left ventricle. Studies from clinical electrophysiology and animal experiments have found that HFpEF is associated with impaired calcium homeostasis, ion channel remodeling and concentric left ventricle hypertrophy (LVH). However, it is still unclear how the abnormal calcium homeostasis, ion channel and structural remodeling affect the electro-mechanical dynamics of the ventricles. In this study we have developed multiscale models of the human left ventricle from single cells to the 3D organ, which take into consideration HFpEF-induced changes in calcium handling, ion channel remodeling and concentric LVH. Our simulation results suggest that at the cellular level, HFpEF reduces the systolic calcium level resulting in a reduced systolic contractile force, but elevates the diastolic calcium level resulting in an abnormal residual diastolic force. In our simulations, these abnormal electro-mechanical features of the ventricular cells became more pronounced with the increase of the heart rate. However, at the 3D organ level, the ejection fraction of the left ventricle was maintained due to the concentric LVH. The simulation results of this study mirror clinically observed features of HFpEF and provide new insights toward the understanding of the cellular bases of impaired cardiac electromechanical functions in heart failure.

Investigation of therapeutic potential and molecular mechanism of vitamin P and digoxin in I/R-induced myocardial infarction in rat

Ischemic-reperfusion (I/R) is a major event in the pathogenesis of ischemic heart disease that leads to higher rate of mortality. The study has been designed to investigate the therapeutic potential and molecular mechanism of vitamin P and digoxin in I/R-induced myocardial infarction in isolated rat heart preparation by using Langendorff apparatus. The animals were treated with vitamin P (50 and 100 mg/kg; p.o.) and digoxin (500 μg/kg) for 5 consecutive days. Digoxin served as a positive control in the present study. On the sixth day, the heart was harvested and induced to 30 min of global ischemia followed by 120 min of reperfusion using Langendorff apparatus. The coronary effluent was collected at different time intervals (i.e. basal, 1, 15, 30, 45, 60 and 120 min.) for the assessment of myocardial contractility function. In addition, creatine kinase-M and B subunits (CK-MB), lactate dehydrogenase (LDH1) and Na(+)-K(+)-ATPase activity along with oxidative tissue biomarkers (i.e. thio-barbituric acid reactive substances (TBARS) and reduced glutathione (GSH)) changes were estimated. The I/R of myocardium produced decrease in coronary flow rate; increase in CK-MB, LDH1 and Na(+)-K(+)-ATPase activity along with increase in TBARS and decrease in GSH levels as compared to normal group. The treatment with vitamin P (100 mg/kg) and digoxin (500 μg/kg) have produced a significant (p < 0.05) ameliorative effect against I/R induced above functional, metabolic and tissue biomarkers changes. Vitamin P has an ameliorative potential against I/R induced myocardial functional changes. It may be due to its free radical scavenging and anti-infarct property via inhibition of Na(+)-K(+)-ATPase activity. Therefore, it can be used as a potential therapeutic medicine for the management of cardiovascular disorders.

Inter-subject variability in human atrial action potential in sinus rhythm versus chronic atrial fibrillation

AIMS

Human atrial electrophysiology exhibits high inter-subject variability in both sinus rhythm (SR) and chronic atrial fibrillation (cAF) patients. Variability is however rarely investigated in experimental and theoretical electrophysiological studies, thus hampering the understanding of its underlying causes but also its implications in explaining differences in the response to disease and treatment. In our study, we aim at investigating the ability of populations of human atrial cell models to capture the inter-subject variability in action potential (AP) recorded in 363 patients both under SR and cAF conditions.

METHODS AND RESULTS

Human AP recordings in atrial trabeculae (n = 469) from SR and cAF patients were used to calibrate populations of computational SR and cAF atrial AP models. Three populations of over 2000 sampled models were generated, based on three different human atrial AP models. Experimental calibration selected populations of AP models yielding AP with morphology and duration in range with experimental recordings. Populations using the three original models can mimic variability in experimental AP in both SR and cAF, with median conductance values in SR for most ionic currents deviating less than 30% from their original peak values. All cAF populations show similar variations in G(K1), G(Kur) and G(to), consistent with AF-related remodeling as reported in experiments. In all SR and cAF model populations, inter-subject variability in I(K1) and I(NaK) underlies variability in APD90, variability in I(Kur), I(CaL) and I(NaK) modulates variability in APD50 and combined variability in Ito and I(Kur) determines variability in APD20. The large variability in human atrial AP triangulation is mostly determined by I(K1) and either I(NaK) or I(NaCa) depending on the model.

CONCLUSION

Experimentally-calibrated human atrial AP models populations mimic AP variability in SR and cAF patient recordings, and identify potential ionic determinants of inter-subject variability in human atrial AP duration and morphology in SR versus cAF.

Comprehensive analyses of ventricular myocyte models identify targets exhibiting favorable rate dependence

Reverse rate dependence is a problematic property of antiarrhythmic drugs that prolong the cardiac action potential (AP). The prolongation caused by reverse rate dependent agents is greater at slow heart rates, resulting in both reduced arrhythmia suppression at fast rates and increased arrhythmia risk at slow rates. The opposite property, forward rate dependence, would theoretically overcome these parallel problems, yet forward rate dependent (FRD) antiarrhythmics remain elusive. Moreover, there is evidence that reverse rate dependence is an intrinsic property of perturbations to the AP. We have addressed the possibility of forward rate dependence by performing a comprehensive analysis of 13 ventricular myocyte models. By simulating populations of myocytes with varying properties and analyzing population results statistically, we simultaneously predicted the rate-dependent effects of changes in multiple model parameters. An average of 40 parameters were tested in each model, and effects on AP duration were assessed at slow (0.2 Hz) and fast (2 Hz) rates. The analysis identified a variety of FRD ionic current perturbations and generated specific predictions regarding their mechanisms. For instance, an increase in L-type calcium current is FRD when this is accompanied by indirect, rate-dependent changes in slow delayed rectifier potassium current. A comparison of predictions across models identified inward rectifier potassium current and the sodium-potassium pump as the two targets most likely to produce FRD AP prolongation. Finally, a statistical analysis of results from the 13 models demonstrated that models displaying minimal rate-dependent changes in AP shape have little capacity for FRD perturbations, whereas models with large shape changes have considerable FRD potential. This can explain differences between species and between ventricular cell types. Overall, this study provides new insights, both specific and general, into the determinants of AP duration rate dependence, and illustrates a strategy for the design of potentially beneficial antiarrhythmic drugs.

Several drugs intended to treat cardiac arrhythmias have failed because of unfavorable rate-dependent properties. That is, the drugs fail to alter electrical activity at fast heart rates, where this would be beneficial, but they do affect electrical activity at slow rates, where this is unwanted. In targeted studies, several agents have been shown to exhibit these unfavorable properties, suggesting that these rate-dependent responses may be intrinsic to ventricular muscle. To determine whether drugs with desirable rate-dependent properties could be rationally designed, we performed comprehensive and systematic analyses of several heart cell models. These analyses calculated the rate-dependent properties of changes in any model parameter, thereby generating simultaneously a large number of model predictions. The analyses showed that targets with favorable rate-dependent properties could indeed be identified, and further simulations uncovered the mechanisms underlying these behaviors. Moreover, a quantitative comparison of results obtained in different models provided new insight in why a given drug applied to different species, or to different tissue types, might produce different rate-dependent behaviors. Overall this study shows how a comprehensive and systematic approach to heart cell models can both identify novel targets and produce more general insight into rate-dependent alterations to cardiac electrical activity.