Cardiac length-dependent activation driven by force-dependent thick-filament dynamics.

The length-dependent activation (LDA) of maximum force and calcium sensitivity are established features of cardiac muscle contraction but the dominant underlying mechanisms remain to be fully clarified. Alongside the well-documented regulation of contraction via the thin filaments, experiments have identified an additional force-dependent thick-filament activation, whereby myosin heads parked in a so-called off state become available to generate force. This process produces a feedback effect that may potentially drive LDA. Using biomechanical modeling of a human left-ventricular myocyte, this study investigates the extent to which the off-state dynamics could, by itself, plausibly account for LDA, depending on the specific mathematical formulation of the feedback. We hypothesized four different models of the off-state regulatory feedback based on (A) total force, (B) active force, (C) sarcomere strain, and (D) passive force. We tested if these models could reproduce the isometric steady-state and dynamic LDA features predicted by an earlier published model of a human left-ventricle myocyte featuring purely phenomenological length dependences. The results suggest that only total-force feedback (A) is capable of reproducing the expected behaviors, but that passive tension could provide a length-dependent signal on which to initiate the feedback. Furthermore, by attributing LDA to off-state dynamics, our proposed model also qualitatively reproduces experimentally observed effects of the off-state-stabilizing drug mavacamten. Taken together, these results support off-state dynamics as a plausible primary mechanism underlying LDA.

Cell to whole organ global sensitivity analysis on a four-chamber heart electromechanics model using Gaussian processes emulators.

Cardiac pump function arises from a series of highly orchestrated events across multiple scales. Computational electromechanics can encode these events in physics-constrained models. However, the large number of parameters in these models has made the systematic study of the link between cellular, tissue, and organ scale parameters to whole heart physiology challenging. A patient-specific anatomical heart model, or digital twin, was created. Cellular ionic dynamics and contraction were simulated with the Courtemanche-Land and the ToR-ORd-Land models for the atria and the ventricles, respectively. Whole heart contraction was coupled with the circulatory system, simulated with CircAdapt, while accounting for the effect of the pericardium on cardiac motion. The four-chamber electromechanics framework resulted in 117 parameters of interest. The model was broken into five hierarchical sub-models: tissue electrophysiology, ToR-ORd-Land model, Courtemanche-Land model, passive mechanics and CircAdapt. For each sub-model, we trained Gaussian processes emulators (GPEs) that were then used to perform a global sensitivity analysis (GSA) to retain parameters explaining 90% of the total sensitivity for subsequent analysis. We identified 45 out of 117 parameters that were important for whole heart function. We performed a GSA over these 45 parameters and identified the systemic and pulmonary peripheral resistance as being critical parameters for a wide range of volumetric and hemodynamic cardiac indexes across all four chambers. We have shown that GPEs provide a robust method for mapping between cellular properties and clinical measurements. This could be applied to identify parameters that can be calibrated in patient-specific models or digital twins, and to link cellular function to clinical indexes.

Effect of scar and His-Purkinje and myocardium conduction on response to conduction system pacing.

INTRODUCTION: Conduction system pacing (CSP), in the form of His bundle pacing (HBP) or left bundle branch pacing (LBBP), is emerging as a valuable cardiac resynchronization therapy (CRT) delivery method. However, patient selection and therapy personalization for CSP delivery remain poorly characterized. We aim to compare pacing-induced electrical synchrony during CRT, HBP, LBBP, HBP with left ventricular (LV) epicardial lead (His-optimized CRT [HOT-CRT]), and LBBP with LV epicardial lead (LBBP-optimized CRT [LOT-CRT]) in patients with different conduction disease presentations using computational modeling. METHODS: We simulated ventricular activation on 24 four-chamber heart geometries, including His-Purkinje systems with proximal left bundle branch block (LBBB). We simulated septal scar, LV lateral wall scar, and mild and severe myocardium and LV His-Purkinje system conduction disease by decreasing the conduction velocity (CV) down to 70% and 35% of the healthy CV. Electrical synchrony was measured by the shortest interval to activate 90% of the ventricles (90% of biventricular activation time [BIVAT-90]). RESULTS: Severe LV His-Purkinje conduction disease favored CRT (BIVAT-90: HBP 101.5 ± 7.8 ms vs. CRT 93.0 ± 8.9 ms, p < .05), with additional electrical synchrony induced by HOT-CRT (87.6 ± 6.7 ms, p < .05) and LOT-CRT (73.9 ± 7.6 ms, p < .05). Patients with slow myocardium CV benefit more from CSP compared to CRT (BIVAT-90: CRT 134.5 ± 24.1 ms; HBP 97.1 ± 9.9 ms, p < .01; LBBP: 101.5 ± 10.7 ms, p < .01). Septal but not lateral wall scar made CSP ineffective, while CRT was able to resynchronize the ventricles in the presence of septal scar (BIVAT-90: baseline 119.1 ± 10.8 ms vs. CRT 85.1 ± 14.9 ms, p < .01). CONCLUSION: Severe LV His-Purkinje conduction disease attenuates the benefits of CSP, with additional improvements achieved with HOT-CRT and LOT-CRT. Septal but not lateral wall scars make CSP ineffective.

Cost-effectiveness analysis of leadless cardiac resynchronization therapy.

BACKGROUND: The Wireless Stimulation Endocardially for CRT (WiSE-CRT) system is a novel technology used to treat patients with dyssynchronous heart failure (HF) by providing leadless cardiac resynchronization therapy (CRT). Observational studies have demonstrated its safety and efficacy profile, however, the treatment cost-effectiveness has not previously been examined. METHODS: A cost-effectiveness evaluation of the WiSE-CRT System was performed using a cohort-based economic model adopting a "proportion in state" structure. In addition to the primary analysis, scenario analyses and sensitivity analyses were performed to test for uncertainty in input parameters. Outcomes were quantified in terms of quality-adjusted life year (QALY) differences. RESULTS: The primary analysis demonstrated that treatment with the WiSE-CRT system is likely to be cost-effective over a lifetime horizon at a QALY reimbursement threshold of £20 000, with a net monetary benefit (NMB) of £3781 per QALY. Cost-effectiveness declines at time horizons shorter than 10 years. Sensitivity analyses demonstrated that average system battery life had the largest impact on potential cost-effectiveness. CONCLUSION: Within the model limitations, these findings support the use of WiSE-CRT in indicated patients from an economic standpoint. However, improving battery technology should be prioritized to maximize cost-effectiveness in times when health services are under significant financial pressures.

Modelling the interaction between stem cells derived cardiomyocytes patches and host myocardium to aid non-arrhythmic engineered heart tissue design.

Application of epicardial patches constructed from human-induced pluripotent stem cell- derived cardiomyocytes (hiPSC-CMs) has been proposed as a long-term therapy to treat scarred hearts post myocardial infarction (MI). Understanding electrical interaction between engineered heart tissue patches (EHT) and host myocardium represents a key step toward a successful patch engraftment. EHT retain different electrical properties with respect to the host heart tissue due to the hiPSC-CMs immature phenotype, which may lead to increased arrhythmia risk. We developed a modelling framework to examine the influence of patch design on electrical activation at the engraftment site. We performed an in silico investigation of different patch design approaches to restore pre-MI activation properties and evaluated the associated arrhythmic risk. We developed an in silico cardiac electrophysiology model of a transmural cross section of host myocardium. The model featured an infarct region, an epicardial patch spanning the infarct region and a bath region. The patch is modelled as a layer of hiPSC-CM, combined with a layer of conductive polymer (CP). Tissue and patch geometrical dimensions and conductivities were incorporated through 10 modifiable model parameters. We validated our model against 4 independent experimental studies and showed that it can qualitatively reproduce their findings. We performed a global sensitivity analysis (GSA) to isolate the most important parameters, showing that the stimulus propagation is mainly governed by the scar depth, radius and conductivity when the scar is not transmural, and by the EHT patch conductivity when the scar is transmural. We assessed the relevance of small animal studies to humans by comparing simulations of rat, rabbit and human myocardium. We found that stimulus propagation paths and GSA sensitivity indices are consistent across species. We explored which EHT design variables have the potential to restore physiological propagation. Simulations predict that increasing EHT conductivity from 0.28 to 1-1.1 S/m recovered physiological activation in rat, rabbit and human. Finally, we assessed arrhythmia risk related to increasing EHT conductivity and tested increasing the EHT Na+ channel density as an alternative strategy to match healthy activation. Our results revealed a greater arrhythmia risk linked to increased EHT conductivity compared to increased Na+ channel density. We demonstrated that our modeling framework could capture the interaction between host and EHT patches observed in in vitro experiments. We showed that large (patch and tissue dimensions) and small (cardiac myocyte electrophysiology) scale differences between small animals and humans do not alter EHT patch effect on infarcted tissue. Our model revealed that only when the scar is transmural do EHT properties impact activation times and isolated the EHT conductivity as the main parameter influencing propagation. We predicted that restoring physiological activation by tuning EHT conductivity is possible but may promote arrhythmic behavior. Finally, our model suggests that acting on hiPSC-CMs low action potential upstroke velocity and lack of IK1 may restore pre-MI activation while not promoting arrhythmia.

Detection of focal source and arrhythmogenic substrate from body surface potentials to guide atrial fibrillation ablation.

Focal sources (FS) are believed to be important triggers and a perpetuation mechanism for paroxysmal atrial fibrillation (AF). Detecting FS and determining AF sustainability in atrial tissue can help guide ablation targeting. We hypothesized that sustained rotors during FS-driven episodes indicate an arrhythmogenic substrate for sustained AF, and that non-invasive electrical recordings, like electrocardiograms (ECGs) or body surface potential maps (BSPMs), could be used to detect FS and AF sustainability. Computer simulations were performed on five bi-atrial geometries. FS were induced by pacing at cycle lengths of 120-270 ms from 32 atrial sites and four pulmonary veins. Self-sustained reentrant activities were also initiated around the same 32 atrial sites with inexcitable cores of radii of 0, 0.5 and 1 cm. FS fired for two seconds and then AF inducibility was tested by whether activation was sustained for another second. ECGs and BSPMs were simulated. Equivalent atrial sources were extracted using second-order blind source separation, and their cycle length, periodicity and contribution, were used as features for random forest classifiers. Longer rotor duration during FS-driven episodes indicates higher AF inducibility (area under ROC curve = 0.83). Our method had accuracy of 90.6±1.0% and 90.6±0.6% in detecting FS presence, and 93.1±0.6% and 94.2±1.2% in identifying AF sustainability, and 80.0±6.6% and 61.0±5.2% in determining the atrium of the focal site, from BSPMs and ECGs of five atria. The detection of FS presence and AF sustainability were insensitive to vest placement (±9.6%). On pre-operative BSPMs of 52 paroxysmal AF patients, patients classified with initiator-type FS on a single atrium resulted in improved two-to-three-year AF-free likelihoods (p-value < 0.01, logrank tests). Detection of FS and arrhythmogenic substrate can be performed from ECGs and BSPMs, enabling non-invasive mapping towards mechanism-targeted AF treatment, and malignant ectopic beat detection with likely AF progression.

Correction: Linking statistical shape models and simulated function in the healthy adult human heart.

[This corrects the article DOI: 10.1371/journal.pcbi.1008851.].

Quantitative mapping of force-pCa curves to whole-heart contraction and relaxation.

The force-pCa (F-pCa) curve is used to characterize steady-state contractile properties of cardiac muscle cells in different physiological, pathological and pharmacological conditions. This provides a reduced preparation in which to isolate sarcomere mechanisms. However, it is unclear how changes in the F-pCa curve impact emergent whole-heart mechanics quantitatively. We study the link between sarcomere and whole-heart function using a multiscale mathematical model of rat biventricular mechanics that describes sarcomere, tissue, anatomy, preload and afterload properties quantitatively. We first map individual cell-level changes in sarcomere-regulating parameters to organ-level changes in the left ventricular function described by pressure-volume loop characteristics (e.g. end-diastolic and end-systolic volumes, ejection fraction and isovolumetric relaxation time). We next map changes in the sarcomere-regulating parameters to changes in the F-pCa curve. We demonstrate that a change in the F-pCa curve can be caused by multiple different changes in sarcomere properties. We demonstrate that changes in sarcomere properties cause non-linear and, importantly, non-monotonic changes in left ventricular function. As a result, a change in sarcomere properties yielding changes in the F-pCa curve that improve contractility does not guarantee an improvement in whole-heart function. Likewise, a desired change in whole-heart function (i.e. ejection fraction or relaxation time) is not caused by a unique shift in the F-pCa curve. Changes in the F-pCa curve alone cannot be used to predict the impact of a compound on whole-heart function. KEY POINTS: The force-pCa (F-pCa) curve is used to assess myofilament calcium sensitivity after pharmacological modulation and to infer pharmacological effects on whole-heart function. We demonstrate that there is a non-unique mapping from changes in F-pCa curves to changes in left ventricular (LV) function. The effect of changes in F-pCa on LV function depend on the state of the heart and could be different for different pathological conditions. Screening of compounds to impact whole-heart function by F-pCa should be combined with active tension and calcium transient measurements to predict better how changes in muscle function will impact whole-heart physiology.

Hypokalemia Promotes Arrhythmia by Distinct Mechanisms in Atrial and Ventricular Myocytes.

RATIONALE: Hypokalemia occurs in up to 20% of hospitalized patients and is associated with increased incidence of ventricular and atrial fibrillation. It is unclear whether these differing types of arrhythmia result from direct and perhaps distinct effects of hypokalemia on cardiomyocytes. OBJECTIVE: To investigate proarrhythmic mechanisms of hypokalemia in ventricular and atrial myocytes. METHODS AND RESULTS: Experiments were performed in isolated rat myocytes exposed to simulated hypokalemia conditions (reduction of extracellular [K+] from 5.0 to 2.7 mmol/L) and supported by mathematical modeling studies. Ventricular cells subjected to hypokalemia exhibited Ca2+ overload and increased generation of both spontaneous Ca2+ waves and delayed afterdepolarizations. However, similar Ca2+-dependent spontaneous activity during hypokalemia was only observed in a minority of atrial cells that were observed to contain t-tubules. This effect was attributed to close functional pairing of the Na+-K+ ATPase and Na+-Ca2+ exchanger proteins within these structures, as reduction in Na+ pump activity locally inhibited Ca2+ extrusion. Ventricular myocytes and tubulated atrial myocytes additionally exhibited early afterdepolarizations during hypokalemia, associated with Ca2+ overload. However, early afterdepolarizations also occurred in untubulated atrial cells, despite Ca2+ quiescence. These phase-3 early afterdepolarizations were rather linked to reactivation of nonequilibrium Na+ current, as they were rapidly blocked by tetrodotoxin. Na+ current-driven early afterdepolarizations in untubulated atrial cells were enabled by membrane hyperpolarization during hypokalemia and short action potential configurations. Brief action potentials were in turn maintained by ultra-rapid K+ current (IKur); a current which was found to be absent in tubulated atrial myocytes and ventricular myocytes. CONCLUSIONS: Distinct mechanisms underlie hypokalemia-induced arrhythmia in the ventricle and atrium but also vary between atrial myocytes depending on subcellular structure and electrophysiology.

In silico identification of potential calcium dynamics and sarcomere targets for recovering left ventricular function in rat heart failure with preserved ejection fraction.

Heart failure with preserved ejection fraction (HFpEF) is a complex disease associated with multiple co-morbidities, where impaired cardiac mechanics are often the end effect. At the cellular level, cardiac mechanics can be pharmacologically manipulated by altering calcium signalling and the sarcomere. However, the link between cellular level modulations and whole organ pump function is incompletely understood. Our goal is to develop and use a multi-scale computational cardiac mechanics model of the obese ZSF1 HFpEF rat to identify important biomechanical mechanisms that underpin impaired cardiac function and to predict how whole-heart mechanical function can be recovered through altering cellular calcium dynamics and/or cellular contraction. The rat heart was modelled using a 3D biventricular biomechanics model. Biomechanics were described by 16 parameters, corresponding to intracellular calcium transient, sarcomere dynamics, cardiac tissue and hemodynamics properties. The model simulated left ventricular (LV) pressure-volume loops that were described by 14 scalar features. We trained a Gaussian process emulator to map the 16 input parameters to each of the 14 outputs. A global sensitivity analysis was performed, and identified calcium dynamics and thin and thick filament kinetics as key determinants of the organ scale pump function. We employed Bayesian history matching to build a model of the ZSF1 rat heart. Next, we recovered the LV function, described by ejection fraction, peak pressure, maximum rate of pressure rise and isovolumetric relaxation time constant. We found that by manipulating calcium, thin and thick filament properties we can recover 34%, 28% and 24% of the LV function in the ZSF1 rat heart, respectively, and 39% if we manipulate all of them together. We demonstrated how a combination of biophysically based models and their derived emulators can be used to identify potential pharmacological targets. We predicted that cardiac function can be best recovered in ZSF1 rats by desensitising the myofilament and reducing the affinity to intracellular calcium concentration and overall prolonging the sarcomere staying in the active force generating state.

Linking statistical shape models and simulated function in the healthy adult human heart.

Cardiac anatomy plays a crucial role in determining cardiac function. However, there is a poor understanding of how specific and localised anatomical changes affect different cardiac functional outputs. In this work, we test the hypothesis that in a statistical shape model (SSM), the modes that are most relevant for describing anatomy are also most important for determining the output of cardiac electromechanics simulations. We made patient-specific four-chamber heart meshes (n = 20) from cardiac CT images in asymptomatic subjects and created a SSM from 19 cases. Nine modes captured 90% of the anatomical variation in the SSM. Functional simulation outputs correlated best with modes 2, 3 and 9 on average (R = 0.49 ± 0.17, 0.37 ± 0.23 and 0.34 ± 0.17 respectively). We performed a global sensitivity analysis to identify the different modes responsible for different simulated electrical and mechanical measures of cardiac function. Modes 2 and 9 were the most important for determining simulated left ventricular mechanics and pressure-derived phenotypes. Mode 2 explained 28.56 ± 16.48% and 25.5 ± 20.85, and mode 9 explained 12.1 ± 8.74% and 13.54 ± 16.91% of the variances of mechanics and pressure-derived phenotypes, respectively. Electrophysiological biomarkers were explained by the interaction of 3 ± 1 modes. In the healthy adult human heart, shape modes that explain large portions of anatomical variance do not explain equivalent levels of electromechanical functional variation. As a result, in cardiac models, representing patient anatomy using a limited number of modes of anatomical variation can cause a loss in accuracy of simulated electromechanical function.

A Quantitative Systems Pharmacology Perspective on the Importance of Parameter Identifiability.

There is an inherent tension in Quantitative Systems Pharmacology (QSP) between the need to incorporate mathematical descriptions of complex physiology and drug targets with the necessity of developing robust, predictive and well-constrained models. In addition to this, there is no "gold standard" for model development and assessment in QSP. Moreover, there can be confusion over terminology such as model and parameter identifiability; complex and simple models; virtual populations; and other concepts, which leads to potential miscommunication and misapplication of methodologies within modeling communities, both the QSP community and related disciplines. This perspective article highlights the pros and cons of using simple (often identifiable) vs. complex (more physiologically detailed but often non-identifiable) models, as well as aspects of parameter identifiability, sensitivity and inference methodologies for model development and analysis. The paper distills the central themes of the issue of identifiability and optimal model size and discusses open challenges.

Leadless left ventricular endocardial pacing for CRT upgrades in previously failed and high-risk patients in comparison with coronary sinus CRT upgrades.

AIMS: Cardiac resynchronization therapy (CRT) upgrades may be less likely to improve following intervention. Leadless left ventricular (LV) endocardial pacing has been used for patients with previously failed CRT or high-risk upgrades. We compared procedural and long-term outcomes in patients undergoing coronary sinus (CS) CRT upgrades with high-risk and previously failed CRT upgrades undergoing LV endocardial upgrades. METHOD AND RESULTS: Prospective consecutive CS upgrades between 2015 and 2019 were compared with those undergoing WiSE-CRT implantation. Cardiac resynchronization therapy response at 6 months was defined as improvement in clinical composite score (CCS) and a reduction in LV end-systolic volume (LVESV) ≥15%. A total of 225 patients were analysed; 121 CS and 104 endocardial upgrades. Patients receiving WiSE-CRT tended to have more comorbidities and were more likely to have previous cardiac surgery (30.9% vs. 16.5%; P = 0.012), hypertension (59.2% vs. 34.7%; P < 0.001), chronic obstructive airways disease (19.4% vs. 9.9%; P = 0.046), and chronic kidney disease (46.4% vs. 21.5%; P < 0.01) but similar LV ejection fraction (30.0 ± 8.3% vs. 29.5 ± 8.6%; P = 0.678). WiSE-CRT upgrades were successful in 97.1% with procedure-related mortality in 1.9%. Coronary sinus upgrades were successful in 97.5% of cases with a 2.5% rate of CS dissection and 5.6% lead malfunction/displacement. At 6 months, 91 WiSE-CRT upgrades and 107 CS upgrades had similar improvements in CCS (76.3% vs. 68.5%; P = 0.210) and reduction in LVESV ≥15% (54.2% vs. 56.3%; P = 0.835). CONCLUSION: Despite prior failed upgrades and high-risk patients with more comorbidities, WiSE-CRT upgrades had high rates of procedural success and similar improvements in CCS and LV remodelling with CS upgrades.

Tracking the motion of intracardiac structures aids the development of future leadless pacing systems.

BACKGROUND: Leadless pacemakers preclude the need for permanent leads to pace endocardium. However, it is yet to be determined whether a leadless pacemaker of a similar design to those manufactured for the right ventricle (RV) fits within the left ventricle (LV), without interfering with intracardiac structures. METHODS: Cardiac computed tomography scans were obtained from 30 patients indicated for cardiac resynchronisation therapy upgrade. The mitral valve annulus, chordae tendineae, papillary muscles and LV endocardial wall were marked in the end-diastolic frame. Intracardiac structures motions were tracked through the cardiac cycle. Two pacemaker designs similar to commercially manufactured leadless systems (Abbott's Nanostim LCP and Medtronic's Micra TPS) as well as theoretical designs with calculated optimal dimensions were evaluated. Pacemakers were virtually placed across the LV endocardial surface and collisions between them and intracardiac structures were detected throughout the cycle. RESULTS: Probability maps of LV intracardiac structures collisions on a 16-segment AHA model indicated possible placement for the Nanostim LCP, Micra TPS, and theoretical designs. Thresholding these maps at a 20% chance of collision revealed only about 36% of the endocardial surface remained collision-free with the deployment of Micra TPS design. The same threshold left no collision-free surface in the case of the Nanostim LCP. To reach at least half of the LV endocardium, the volume of Micra TPS, which is the smaller design, needed to be decreased by 41%. CONCLUSION: Due to the presence of intracardiac structures, placement of leadless pacemakers with dimensions similar to commercially manufactured RV systems would be limited to apical regions.

Gaussian process manifold interpolation for probabilistic atrial activation maps and uncertain conduction velocity.

In patients with atrial fibrillation, local activation time (LAT) maps are routinely used for characterizing patient pathophysiology. The gradient of LAT maps can be used to calculate conduction velocity (CV), which directly relates to material conductivity and may provide an important measure of atrial substrate properties. Including uncertainty in CV calculations would help with interpreting the reliability of these measurements. Here, we build upon a recent insight into reduced-rank Gaussian processes (GPs) to perform probabilistic interpolation of uncertain LAT directly on human atrial manifolds. Our Gaussian process manifold interpolation (GPMI) method accounts for the topology of the atrium, and allows for calculation of statistics for predicted CV. We demonstrate our method on two clinical cases, and perform validation against a simulated ground truth. CV uncertainty depends on data density, wave propagation direction and CV magnitude. GPMI is suitable for probabilistic interpolation of other uncertain quantities on non-Euclidean manifolds. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Balance of Active, Passive, and Anatomical Cardiac Properties in Doxorubicin-Induced Heart Failure.

Late-onset heart failure (HF) is a known side effect of doxorubicin chemotherapy. Typically, patients are diagnosed when already at an irreversible stage of HF, which allows few or no treatment options. Identifying the causes of compromised cardiac function in this patient group may improve early patient diagnosis and support treatment selection. To link doxorubicin-induced changes in cardiac cellular and tissue mechanical properties to overall cardiac function, we apply a multiscale biophysical biomechanics model of the heart to measure the plausibility of changes in model parameters representing the passive, active, or anatomical properties of the left ventricle for reproducing measured patient phenotypes. We create representative models of healthy controls (N = 10) and patients with HF induced by (N = 22) or unrelated to (N = 25) doxorubicin therapy. The model predicts that HF in the absence of doxorubicin is characterized by a 2- to 3-fold stiffness increase, decreased tension (0-20%), and ventricular dilation (of order 10-30%). HF due to doxorubicin was similar but showed stronger bias toward reduced active contraction (10-30%) and less dilation (0-20%). We find that changes in active, passive, and anatomical properties all play a role in doxorubicin-induced cardiotoxicity phenotypes. Differences in parameter changes between patient groups are consistent with doxorubicin cardiotoxicity having a greater dependence on reduced cellular contraction and less anatomical remodeling than HF not caused by doxorubicin.

Sex-Dependent QRS Guidelines for Cardiac Resynchronization Therapy Using Computer Model Predictions.

Cardiac resynchronization therapy (CRT) is an important treatment for heart failure. Low female enrollment in clinical trials means that current CRT guidelines may be biased toward males. However, females have higher response rates at lower QRS duration (QRSd) thresholds. Sex differences in the left ventricle (LV) size could provide an explanation for the improved female response at lower QRSd. We aimed to test if sex differences in CRT response at lower QRSd thresholds are explained by differences in LV size and hence predict sex-specific guidelines for CRT. We investigated the effect that LV size sex difference has on QRSd between male and females in 1093 healthy individuals and 50 CRT patients using electrophysiological computer models of the heart. Simulations on the healthy mean shape models show that LV size sex difference can account for 50-100% of the sex difference in baseline QRSd in healthy individuals. In the CRT patient cohort, model simulations predicted female-specific guidelines for CRT, which were 9-13 ms lower than current guidelines. Sex differences in the LV size are able to account for a significant proportion of the sex difference in QRSd and provide a mechanistic explanation for the sex difference in CRT response. Simulations accounting for the smaller LV size in female CRT patients predict 9-13 ms lower QRSd thresholds for female CRT guidelines.

A short history of the development of mathematical models of cardiac mechanics.

Cardiac mechanics plays a crucial role in atrial and ventricular function, in the regulation of growth and remodelling, in the progression of disease, and the response to treatment. The spatial scale of the critical mechanisms ranges from nm (molecules) to cm (hearts) with the fastest events occurring in milliseconds (molecular events) and the slowest requiring months (growth and remodelling). Due to its complexity and importance, cardiac mechanics has been studied extensively both experimentally and through mathematical models and simulation. Models of cardiac mechanics evolved from seminal studies in skeletal muscle, and developed into cardiac specific, species specific, human specific and finally patient specific calculations. These models provide a formal framework to link multiple experimental assays recorded over nearly 100 years into a single unified representation of cardiac function. This review first provides a summary of the proteins, physiology and anatomy involved in the generation of cardiac pump function. We then describe the evolution of models of cardiac mechanics starting with the early theoretical frameworks describing the link between sarcomeres and muscle contraction, transitioning through myosin-level models to calcium-driven systems, and ending with whole heart patient-specific models.

A comprehensive multi-index cardiac magnetic resonance-guided assessment of atrial fibrillation substrate prior to ablation: Prediction of long-term outcomes.

INTRODUCTION: Multiple cardiac magnetic resonance (CMR)-derived indices of atrial fibrillation (AF) substrate have been shown in isolation to predict long-term outcome following catheter ablation. Left atrial (LA) fibrosis, LA volume, LA ejection fraction (EF), left ventricular ejection fraction (LVEF), LA shape (sphericity) and pulmonary vein anatomy have all been shown to correlate with late AF recurrence. This study aimed to validate and assess the relative contribution of multiple indices in a long-term single-center study. METHODS AND RESULTS: Eighty-nine patients (53% paroxysmal AF, 73% male) underwent comprehensive CMR study before first-time AF ablation (median follow-up 726 days [IQR: 418-1010 days]). The 3D late gadolinium-enhanced acquisition (1.5T, 1.3 × 1.3 × 2 mm) was quantified for fibrosis; LA volume and sphericity were assessed on manual segmentation at atrial diastole; LAEF and LVEF were quantified on multislice cine imaging. AF recurred in 43 patients (48%) overall (31 at 1 year). In the recurrence group, LA fibrosis was higher (42% vs 29%; hazard ratio [HR]: 1.032; P = .002), left atrial ejection fraction (LAEF) lower (25% vs 34%; HR: 0.063; P = .016) and LVEF lower (57% vs 63%; HR: 0.011; P = .008). LA volume (135 vs 124 mL) and sphericity (0.819 vs 0.822) were similar. Multivariate Cox regression analysis was adjusted for age and sex (Model 1), additionally AF type (Model 2) and combined (Model 3). In Models 1 and 2, LA fibrosis, LAEF, and LVEF were independently associated with outcome, but only LA fibrosis was independent in Model 3 (HR: 1.021; P = .022). CONCLUSIONS: LAEF, LVEF, and LA fibrosis differed significantly in the AF recurrence cohort. However, on combined multivariate analysis only LA fibrosis remained independently associated with outcome.

Normoxic cells remotely regulate the acid-base balance of cells at the hypoxic core of connexin-coupled tumor growths.

ATP fuels the removal of metabolic end-products, including H+ ions that profoundly modulate biological activities. Energetic resources in hypoxic tumor regions are constrained by low-yielding glycolysis, and any means of reducing the cost of acid extrusion, without compromising pH homeostasis, would therefore be advantageous for cancer cells. Some cancers express connexin channels that allow solute exchange between cells, and we propose that, via this route, normoxic cells supply hypoxic neighbors with acid-neutralizing HCO3- ions. This hypothesis was tested by imaging cytoplasmic pH in spheroidal tissue growths of connexin43-positive pancreatic cancer Colo357 cells during light-controlled H+ uncaging at the hypoxic core. Cytoplasmic acid retention at the core was halved in the presence of CO2/HCO3-, but this process requires a restorative HCO3- flux. The effect of CO2/HCO3- was ablated by connexin43 inhibition or knockdown. In connexin-decoupled spheroids, 4,4'-diisothiocyano-2,2'-stilbenedisulfonic acid (DIDS), an inhibitor of HCO3- uptake, had no effect on cytoplasmic [H+] in the H+-uncaging region, indicating that DIDS-sensitive transport is not an adequate pH-regulatory strategy therein. With intact connexin-coupling, acid retention at the core increased upon DIDS treatment, indicating that HCO3- ions are taken up actively by peripheral cells and then transmitted passively to cells at the hypoxic core. Thus, the energetic burden of pH regulation is offloaded from hypoxic cells onto metabolically altruistic normoxic neighbors.-Dovmark, T. H., Hulikova, A., Niederer, S. A., Vaughan-Jones, R. D., Swietach, P. Normoxic cells remotely regulate the acid-base balance of cells at the hypoxic core of connexin-coupled tumor growths.