Predicting left ventricular contractile function via Gaussian process emulation in aortic-banded rats.

Cardiac contraction is the result of integrated cellular, tissue and organ function. Biophysical in silico cardiac models offer a systematic approach for studying these multi-scale interactions. The computational cost of such models is high, due to their multi-parametric and nonlinear nature. This has so far made it difficult to perform model fitting and prevented global sensitivity analysis (GSA) studies. We propose a machine learning approach based on Gaussian process emulation of model simulations using probabilistic surrogate models, which enables model parameter inference via a Bayesian history matching (HM) technique and GSA on whole-organ mechanics. This framework is applied to model healthy and aortic-banded hypertensive rats, a commonly used animal model of heart failure disease. The obtained probabilistic surrogate models accurately predicted the left ventricular pump function (R2 = 0.92 for ejection fraction). The HM technique allowed us to fit both the control and diseased virtual bi-ventricular rat heart models to magnetic resonance imaging and literature data, with model outputs from the constrained parameter space falling within 2 SD of the respective experimental values. The GSA identified Troponin C and cross-bridge kinetics as key parameters in determining both systolic and diastolic ventricular function. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Cardiomyocyte substructure reverts to an immature phenotype during heart failure.

KEY POINTS: As reactivation of the fetal gene program has been implicated in pathological remodelling during heart failure (HF), we examined whether cardiomyocyte subcellular structure and function revert to an immature phenotype during this disease. Surface and internal membrane structures appeared gradually during development, and returned to a juvenile state during HF. Similarly, dyadic junctions between the cell membrane and sarcoplasmic reticulum were progressively 'packed' with L-type Ca2+ channels and ryanodine receptors during development, and 'unpacked' during HF. Despite similarities in subcellular structure, dyads were observed to be functional from early developmental stages, but exhibited an impaired ability to release Ca2+ in failing cardiomyocytes. Thus, while immature and failing cardiomyocytes share similarities in subcellular structure, these do not fully account for the marked impairment of Ca2+ homeostasis observed in HF. ABSTRACT: Reactivation of the fetal gene programme has been implicated as a driver of pathological cardiac remodelling. Here we examined whether pathological remodelling of cardiomyocyte substructure and function during heart failure (HF) reflects a reversion to an immature phenotype. Using scanning electron microscopy, we observed that Z-grooves and t-tubule openings at the cell surface appeared gradually during cardiac development, and disappeared during HF. Confocal and super-resolution imaging within the cell interior revealed similar structural parallels; disorganization of t-tubules in failing cells was strikingly reminiscent of the late stages of postnatal development, with fewer transverse elements and a high proportion of longitudinal tubules. Ryanodine receptors (RyRs) were observed to be laid down in advance of developing t-tubules and similarly 'orphaned' in HF, although RyR distribution along Z-lines was relatively sparse. Indeed, nanoscale imaging revealed coordinated packing of L-type Ca2+ channels and RyRs into dyadic junctions during development, and orderly unpacking during HF. These findings support a 'last in, first out' paradigm, as the latest stages of dyadic structural development are reversed during disease. Paired imaging of t-tubules and Ca2+ showed that the disorganized arrangement of dyads in immature and failing cells promoted desynchronized and slowed Ca2+ release in these two states. However, while developing cells exhibited efficient triggering of Ca2+ release at newly formed dyads, dyadic function was impaired in failing cells despite similar organization of Ca2+ handling proteins. Thus, pathologically deficient Ca2+ homeostasis during HF is only partly linked to the re-emergence of immature subcellular structure, and additionally reflects lost dyadic functionality.

Evidence for heterogeneous subsarcolemmal Na+ levels in rat ventricular myocytes.

The intracellular Na+ concentration ([Na+]) regulates cardiac contractility. Previous studies have suggested that subsarcolemmal [Na+] is higher than cytosolic [Na+] in cardiac myocytes, but this concept remains controversial. Here, we used electrophysiological experiments and mathematical modeling to test whether there are subsarcolemmal pools with different [Na+] and dynamics compared with the bulk cytosol in rat ventricular myocytes. A Na+ dependency curve for Na+-K+-ATPase (NKA) current was recorded with symmetrical Na+ solutions, i.e., the same [Na+] in the superfusate and internal solution. This curve was used to estimate [Na+] sensed by NKA in other experiments. Three experimental observations suggested that [Na+] is higher near NKA than in the bulk cytosol: 1) when extracellular [Na+] was high, [Na+] sensed by NKA was ~6 mM higher than the internal solution in quiescent cells; 2) long trains of Na+ channel activation almost doubled this gradient; compared with an even intracellular distribution of Na+, the increase of [Na+] sensed by NKA was 10 times higher than expected, suggesting a local Na+ domain; and 3) accumulation of Na+ near NKA after trains of Na+ channel activation dissipated very slowly. Finally, mathematical models assuming heterogeneity of [Na+] between NKA and the Na+ channel better reproduced experimental data than the homogeneous model. In conclusion, our data suggest that NKA-sensed [Na+] is higher than [Na+] in the bulk cytosol and that there are differential Na+ pools in the subsarcolemmal space, which could be important for cardiac contractility and arrhythmogenesis. NEW & NOTEWORTHY Our data suggest that the Na+-K+-ATPase-sensed Na+ concentration is higher than the Na+ concentration in the bulk cytosol and that there are differential Na+ pools in the subsarcolemmal space, which could be important for cardiac contractility and arrhythmogenesis. Listen to this article's corresponding podcast at https://ajpheart.podbean.com/e/heterogeneous-sodium-in-ventricular-myocytes/ .

Hypokalaemia induces Ca²âº overload and Ca²âº waves in ventricular myocytes by reducing Na⁺,K⁺-ATPase α2 activity.

KEY POINTS: Hypokalaemia is a risk factor for development of ventricular arrhythmias. In rat ventricular myocytes, low extracellular K(+) (corresponding to clinical moderate hypokalaemia) increased Ca(2+) wave probability, Ca(2+) transient amplitude, sarcoplasmic reticulum (SR) Ca(2+) load and induced SR Ca(2+) leak. Low extracellular K(+) reduced Na(+),K(+)-ATPase (NKA) activity and hyperpolarized the resting membrane potential in ventricular myocytes. Both experimental data and modelling indicate that reduced NKA activity and subsequent Na(+) accumulation sensed by the Na(+), Ca(2+) exchanger (NCX) lead to increased Ca(2+) transient amplitude despite concomitant hyperpolarization of the resting membrane potential. Low extracellular K(+) induced Ca(2+) overload by lowering NKA α2 activity. Triggered ventricular arrhythmias in patients with hypokalaemia may therefore be attributed to reduced NCX forward mode activity linked to an effect on the NKA α2 isoform. ABSTRACT: Hypokalaemia is a risk factor for development of ventricular arrhythmias. The aim of this study was to determine the cellular mechanisms leading to triggering of arrhythmias in ventricular myocytes exposed to low Ko. Low Ko, corresponding to moderate hypokalaemia, increased Ca(2+) transient amplitude, sarcoplasmic reticulum (SR) Ca(2+) load, SR Ca(2+) leak and Ca(2+) wave probability in field stimulated rat ventricular myocytes. The mechanisms leading to Ca(2+) overload were examined. Low Ko reduced Na(+),K(+)-ATPase (NKA) currents, increased cytosolic Na(+) concentration and increased the Na(+) level sensed by the Na(+), Ca(2+) exchanger (NCX). Low Ko also hyperpolarized the resting membrane potential (RMP) without significant alterations in action potential duration. Experiments in voltage clamped and field stimulated ventricular myocytes, along with mathematical modelling, suggested that low Ko increases the Ca(2+) transient amplitude by reducing NKA activity despite hyperpolarization of the RMP. Selective inhibition of the NKA α2 isoform by low dose ouabain abolished the ability of low Ko to reduce NKA currents, to increase Na(+) levels sensed by NCX and to increase the Ca(2+) transient amplitude. We conclude that low Ko, within the range of moderate hypokalaemia, increases Ca(2+) levels in ventricular myocytes by reducing the pumping rate of the NKA α2 isoform with subsequent Na(+) accumulation sensed by the NCX. These data highlight reduced NKA α2 -mediated control of NCX activity as a possible mechanism underlying triggered ventricular arrhythmias in patients with hypokalaemia.

Impaired left ventricular mechanical and energetic function in mice after cardiomyocyte-specific excision of Serca2.

Sarco(endo)plasmic reticulum Ca2+ -ATPase (SERCA)2 transports Ca2+ from the cytosol into the sarcoplasmic reticulum of cardiomyocytes and is essential for maintaining myocardial Ca2+ handling and thus the mechanical function of the heart. SERCA2 is a major ATP consumer in excitation-contraction coupling but is regarded to contribute to energetically efficient Ca2+ handling in the cardiomyocyte. Previous studies using cardiomyocyte-specific SERCA2 knockout (KO) mice have demonstrated that decreased SERCA2 activity reduces the Ca2+ transient amplitude and induces compensatory Ca2+ transport mechanisms that may lead to more inefficient Ca2+ transport. In this study, we examined the relationship between left ventricular (LV) function and myocardial O2 consumption (MVo2) in ex vivo hearts from SERCA2 KO mice to directly measure how SERCA2 elimination influences mechanical and energetic features of the heart. Ex vivo hearts from SERCA2 KO hearts developed mechanical dysfunction at 4 wk and demonstrated virtually no working capacity at 7 wk. In accordance with the reported reduction in Ca2+ transient amplitude in cardiomyocytes from SERCA2 KO mice, work-independent MVo2 was decreased due to a reduced energy cost of excitation-contraction coupling. As these hearts also showed a marked impairment in the efficiency of chemomechanical energy transduction (contractile efficiency, i.e, work-dependent MVo2), hearts from SERCA2 KO mice were found to be mechanically inefficient. This ex vivo evaluation of mechanical and energetic function in hearts from SERCA2 KO mice brings together findings from previous experimental and mathematical modeling-based studies and demonstrates that reduced SERCA2 activity not only leads to mechanical dysfunction but also to energetic dysfunction.

Effects of serum-free storage on morphology, phenotype, and viability of ex vivo cultured human conjunctival epithelium.

The use of amniotic membrane (AM) represents one of the major developments in ocular surface reconstruction. However, in a study on patients with primary pterygium, transplantation of AM with ex vivo expanded human conjunctival epithelial cells (HCjE) promoted earlier epithelialization than AM alone. We previously showed that cultured human limbal epithelial cells maintain their morphology, phenotype, and viability for one week when stored at 23°C. The current study investigates the feasibility of storing HCjE in HEPES-MEM and Optisol-GS at 23°C for 4 and 7 days, respectively. The five experimental groups were analyzed by light microscopy, immunohistochemistry, transmission electron microscopy, and a viability assay. The ultrastructural integrity of cultured HCjE was well preserved following 4 days of storage, however, 7 days of storage resulted in some loss of cell-cell contacts and epithelial detachment from the amniotic membrane. The number of microvilli in cultured HCjE not subjected to storage was 2.03±0.38 microvilli/µm. In comparison, after 4 and 7 days of HEPES-MEM storage this number was 1.69±0.54 microvilli/µm; P=0.98 and 0.89±1.0 microvilli/µm; P=0.28, respectively. After Optisol-GS storage for 4 and 7 days, the mean number of microvilli was 1.07±1.0 microvilli/µm; P=0.47 and 0.07±0.07 microvilli/µm; P=0.03, respectively. The number of cell layers in cultured HCjE not subjected to storage was 4.4±0.3 cell layers, as opposed to 4.0±0.9 cell layers; P=0.89 after 4 days of HEPES-MEM storage and 2.8±0.6 cell layers; P=0.01 after 7 days of storage in HEPES-MEM. The number of cell layers after 4 and 7 days of storage in Optisol-GS was 3.7±0.2 cell layers; P=0.46 and 3.4±0.4 cell layers; P=0.18, respectively. The expression of markers for undifferentiated cells (ΔNp63α, ABCG2 and p63), proliferating cells (Ki67 and PCNA), goblet cells (Ck7 and MUC5AC), stratified squamous epithelial cells (Ck4), and apoptotic cells (caspase-3) in cultured HCjE appeared to be unchanged after 4 and 7 days of HEPES-MEM and Optisol-GS storage. The percentage of viable cells in cultured HCjE not subjected to storage (91.4%±3.2%) was sustained after 4 and 7 days of storage in HEPES-MEM (94.1%±4.5%; P=0.99 and 85.1%±13.7%; P=0.87, respectively) as well as after 4 and 7 days of storage in Optisol-GS (87.7%±15.2%; P=0.97 and 79.8%±15.7%; P=0.48, respectively). We conclude that cultured HCjE may be stored for at least 4 days in serum-free conditions at 23°C while maintaining the phenotype and viability. HEPES-MEM appears to be comparable to Optisol-GS for serum-free storage with preservation of the ultrastructure for at least 4 days.

Calcium dynamics in the ventricular myocytes of SERCA2 knockout mice: A modeling study.

We describe a simulation study of Ca²(+) dynamics in mice with cardiomyocyte-specific conditional excision of the sarco(endo)plasmic reticulum calcium ATPase (SERCA) gene, using an experimental data-driven biophysically-based modeling framework. Previously, we reported a moderately impaired heart function measured in mice at 4 weeks after SERCA2 gene deletion (knockout (KO)), along with a >95% reduction in the level of SERCA2 protein. We also reported enhanced Ca²(+) flux through the L-type Ca²(+) channels and the Na(+)/Ca²(+) exchanger in ventricular myocytes isolated from these mice, compared to the control Serca2(flox/flox) mice (flox-flox (FF)). In the current study, a mathematical model-based analysis was applied to enable further quantitative investigation into changes in the Ca²(+) handling mechanisms in these KO cardiomyocytes. Model parameterization based on a wide range of experimental measurements showed a 67% reduction in SERCA activity and an over threefold increase in the activity of the Na(+)/Ca²(+) exchanger. The FF and KO models were then validated against experimentally measured [Ca²(+)](i) transients and experimentally estimated sarco(endo)plasmic reticulum (SR) function. Simulation results were in quantitative agreement with experimental measurements, confirming that sustained [Ca²(+)](i) transients could be maintained in the KO cardiomyocytes despite severely impaired SERCA function. In silico analysis shows that diastolic [Ca²(+)](i) rises sharply with progressive reductions in SERCA activity at physiologically relevant pacing frequencies. Furthermore, an analysis of the roles of the compensatory mechanisms revealed that the major combined effect of the compensatory mechanisms is to lower diastolic [Ca²(+)](i). Finally, by using a comprehensive sensitivity analysis of the role of all cellular calcium handling mechanisms, we show that the combination of upregulation of the Na(+)/Ca²(+) exchanger and increased L-type Ca²(+) current is the most effective means to maintain diastolic and systolic calcium levels after loss of SERCA function.

Losartan improves recovery of contraction and inhibits transient inward current in a cellular model of cardiac ischemia and reperfusion.

Losartan, a selective angiotensin II (AII) type I receptor antagonist, may protect against myocardial stunning and arrhythmia in ischemia and reperfusion. To examine the cellular basis for these protective actions, we studied effects of losartan and AII on contractile and electrical activity of ventricular myocytes exposed to simulated ischemia and reperfusion. Ionic currents were measured with voltage-clamp techniques and contractions were measured with a video edge detector. After 10 min of superfusion with Tyrode's solution at 37 degrees C, cells were exposed to simulated ischemia (hypoxia, acidosis, hyperkalemia, hypercapnia, lactate accumulation, and substrate deprivation) for 30 min followed by 25 min of reperfusion with normal Tyrode's solution. During ischemia, drug-treated cells were exposed to either 0.1 microM AII, 10 microM losartan, or both simultaneously. In reperfusion, contractions were depressed to 42% of preischemic levels in untreated cells. Losartan treatment significantly improved contractile recovery to 84% (P <. 05) of preischemic levels. AII-treated cells showed contractile recovery similar to untreated cells (40%), whereas cells treated with losartan plus AII recovered to 101% of preischemic levels. Cells exposed to losartan or losartan plus AII also exhibited reduced incidence of transient inward current (I(TI)) (20%, P <.05; 36%) relative to untreated cells (60%). However, I(TI) incidence was not altered by treatment with AII alone (57%). Treatment with exogenous agonist did not potentiate contractile depression or I(TI) incidence, and losartan exerted protective effects in the presence and absence of AII. Thus, losartan may have effects that are independent of AII receptor blockade.