Stretch Harmonizes Sarcomere Strain Across the Cardiomyocyte.

BACKGROUND: Increasing cardiomyocyte contraction during myocardial stretch serves as the basis for the Frank-Starling mechanism in the heart. However, it remains unclear how this phenomenon occurs regionally within cardiomyocytes, at the level of individual sarcomeres. We investigated sarcomere contractile synchrony and how intersarcomere dynamics contribute to increasing contractility during cell lengthening. METHODS: Sarcomere strain and Ca2+ were simultaneously recorded in isolated left ventricular cardiomyocytes during 1 Hz field stimulation at 37 °C, at resting length and following stepwise stretch. RESULTS: We observed that in unstretched rat cardiomyocytes, differential sarcomere deformation occurred during each beat. Specifically, while most sarcomeres shortened during the stimulus, ≈10% to 20% of sarcomeres were stretched or remained stationary. This nonuniform strain was not traced to regional Ca2+ disparities but rather shorter resting lengths and lower force production in systolically stretched sarcomeres. Lengthening of the cell recruited additional shortening sarcomeres, which increased contractile efficiency as less negative, wasted work was performed by stretched sarcomeres. Given the known role of titin in setting sarcomere dimensions, we next hypothesized that modulating titin expression would alter intersarcomere dynamics. Indeed, in cardiomyocytes from mice with titin haploinsufficiency, we observed greater variability in resting sarcomere length, lower recruitment of shortening sarcomeres, and impaired work performance during cell lengthening. CONCLUSIONS: Graded sarcomere recruitment directs cardiomyocyte work performance, and harmonization of sarcomere strain increases contractility during cell stretch. By setting sarcomere dimensions, titin controls sarcomere recruitment, and its lowered expression in haploinsufficiency mutations impairs cardiomyocyte contractility.

Frequent blood flow restricted training not to failure and to failure induces similar gains in myonuclei and muscle mass.

The purpose of the present study was to compare the effects of short-term high-frequency failure vs non-failure blood flow-restricted resistance exercise (BFRRE) on changes in satellite cells (SCs), myonuclei, muscle size, and strength. Seventeen untrained men performed four sets of BFRRE to failure (Failure) with one leg and not to failure (Non-failure; 30-15-15-15 repetitions) with the other leg using knee-extensions at 20% of one repetition maximum (1RM). Fourteen sessions were distributed over two 5-day blocks, separated by a 10-day rest period. Muscle samples obtained before, at mid-training, and 10-day post-intervention (Post10) were analyzed for muscle fiber area (MFA), myonuclei, and SC. Muscle size and echo intensity of m.rectus femoris (RF) and m.vastus lateralis (VL) were measured by ultrasonography, and knee extension strength with 1RM and maximal isometric contraction (MVC) up until Post24. Both protocols increased myonuclear numbers in type-1 (12%-17%) and type-2 fibers (20%-23%), and SC in type-1 (92%-134%) and type-2 fibers (23%-48%) at Post10 (p < 0.05). RF and VL size increased by 5%-10% in both legs at Post10 to Post24, whereas the MFA of type-1 fibers in Failure was decreased at Post10 (-10 ± 16%; p = 0.02). Echo intensity increased by ~20% in both legs during Block1 (p < 0.001) and was ~8 to 11% below baseline at Post24 (p = 0.001-0.002). MVC and 1RM decreased by 5%-10% after Block1, but increased in both legs by 6%-11% at Post24 (p < 0.05). In conclusion, both short-term high-frequency failure and non-failure BFRRE induced increases in SCs, in myonuclei content, muscle size, and strength, concomitant with decreased echo intensity. Intriguingly, the responses were delayed and peaked 10-24 days after the training intervention. Our findings may shed light on the mechanisms involved in resistance exercise-induced overreaching and supercompensation.

Arrhythmogenic Current Generation by Myofilament-Triggered Ca2+ Release and Sarcomere Heterogeneity.

Heterogeneous mechanical dyskinesis has been implicated in many arrhythmogenic phenotypes. Strain-dependent perturbations to cardiomyocyte electrophysiology may contribute to this arrhythmogenesis through processes referred to as mechanoelectric feedback. Although the role of stretch-activated ion currents has been investigated using computational models, experimental studies suggest that mechanical strain may also promote arrhythmia by facilitating calcium wave propagation. To investigate whether strain-dependent changes in calcium affinity to the myofilament may promote arrhythmogenic intracellular calcium waves, we modified a mathematical model of rabbit excitation-contraction coupling coupled to a model of myofilament activation and force development. In a one-dimensional compartmental analysis, we bidirectionally coupled 50 sarcomere models in series to model calcium diffusion and stress transfer between adjacent sarcomeres. These considerations enabled the model to capture 1) the effects of mechanical feedback on calcium homeostasis at the sarcomeric level and 2) the combined effects of mechanical and calcium heterogeneities at the cellular level. The results suggest that in conditions of calcium overload, the vulnerable window of stretch-release to trigger suprathreshold delayed afterdepolarizations can be affected by heterogeneity in sarcomere length. Furthermore, stretch and sarcomere heterogeneity may modulate the susceptibility threshold for delayed afterdepolarizations and the aftercontraction wave propagation velocity.

Computational quantification of patient-specific changes in ventricular dynamics associated with pulmonary hypertension.

Pulmonary arterial hypertension (PAH) causes an increase in the mechanical loading imposed on the right ventricle (RV) that results in progressive changes to its mechanics and function. Here, we quantify the mechanical changes associated with PAH by assimilating clinical data consisting of reconstructed three-dimensional geometry, pressure, and volume waveforms, as well as regional strains measured in patients with PAH (n = 12) and controls (n = 6) within a computational modeling framework of the ventricles. Modeling parameters reflecting regional passive stiffness and load-independent contractility as indexed by the tissue active tension were optimized so that simulation results matched the measurements. The optimized parameters were compared with clinical metrics to find usable indicators associated with the underlying mechanical changes. Peak contractility of the RV free wall (RVFW) γRVFW,max was found to be strongly correlated and had an inverse relationship with the RV and left ventricle (LV) end-diastolic volume ratio (i.e., RVEDV/LVEDV) (RVEDV/LVEDV)+ 0.44, R2 = 0.77). Correlation with RV ejection fraction (R2 = 0.50) and end-diastolic volume index (R2 = 0.40) were comparatively weaker. Patients with with RVEDV/LVEDV > 1.5 had 25% lower γRVFW,max (P < 0.05) than that of the control. On average, RVFW passive stiffness progressively increased with the degree of remodeling as indexed by RVEDV/LVEDV. These results suggest a mechanical basis of using RVEDV/LVEDV as a clinical index for delineating disease severity and estimating RVFW contractility in patients with PAH.NEW & NOTEWORTHY This article presents patient-specific data assimilation of a patient cohort and physical description of clinical observations.

A computational framework for testing arrhythmia marker sensitivities to model parameters in functionally calibrated populations of atrial cells.

Models of cardiac cell electrophysiology are complex non-linear systems which can be used to gain insight into mechanisms of cardiac dynamics in both healthy and pathological conditions. However, the complexity of cardiac models can make mechanistic insight difficult. Moreover, these are typically fitted to averaged experimental data which do not incorporate the variability in observations. Recently, building populations of models to incorporate inter- and intra-subject variability in simulations has been combined with sensitivity analysis (SA) to uncover novel ionic mechanisms and potentially clarify arrhythmogenic behaviors. We used the Koivumäki human atrial cell model to create two populations, representing normal Sinus Rhythm (nSR) and chronic Atrial Fibrillation (cAF), by varying 22 key model parameters. In each population, 14 biomarkers related to the action potential and dynamic restitution were extracted. Populations were calibrated based on distributions of biomarkers to obtain reasonable physiological behavior, and subjected to SA to quantify correlations between model parameters and pro-arrhythmia markers. The two populations showed distinct behaviors under steady state and dynamic pacing. The nSR population revealed greater variability, and more unstable dynamic restitution, as compared to the cAF population, suggesting that simulated cAF remodeling rendered cells more stable to parameter variation and rate adaptation. SA revealed that the biomarkers depended mainly on five ionic currents, with noted differences in sensitivities to these between nSR and cAF. Also, parameters could be selected to produce a model variant with no alternans and unaltered action potential morphology, highlighting that unstable dynamical behavior may be driven by specific cell parameter settings. These results ultimately suggest that arrhythmia maintenance in cAF may not be due to instability in cell membrane excitability, but rather due to tissue-level effects which promote initiation and maintenance of reentrant arrhythmia.

Methodology for Morphometric Analysis of Modern Human Contralateral Premolars.

OBJECTIVE: This study aimed at developing a standard methodology for morphometric analysis and comparison of contralateral human premolar pulp space using microcomputed tomography (micro-CT) and semiautomated software. The primary objective was to establish a method to compare the complex and minute morphological internal volumes of contralateral premolar pulp spaces and determine their degree of similarity. The secondary aim was to introduce new methodology for selecting contralateral premolars for the study of biomaterials and techniques. METHODS: Forty-one intact human premolar pairs (n = 82) extracted from 28 patients were scanned with micro-CT. Quantitative comparative evaluation was performed through geometric morphometric deviation analysis of the pulp spaces after mirroring, automatic alignment, and coregistration with semiautomated software. Geometric parameters compared included volume, surface, and surface over volume. Shape deviation analysis of transformed mean distances and root mean square errors was conducted. RESULTS: The geometric parameters of the contralateral premolar pulp spaces had significantly higher similarity coefficients than random pairs (P < 0.001). The shape deviation analysis and transformed mean distances had significantly lower values for contralaterals compared with random pairs (P < 0.001). The contralateral geometries had a statistically significant narrower distribution in deviation when compared with random pairs (P < 0.001). CONCLUSIONS: We present a methodology that sets a new standard for internal validation of the teeth to be used in ex vivo testing of endodontic materials and techniques. It also shows that the resolution of the CT scan is crucial and that studies using cone beam CT cannot be used for such studies.

Advancing the understanding and treatment of post-traumatic stress disorder with computational modelling.

The existing theories of post-traumatic stress disorder (PTSD) have inspired large volumes of research and have contributed substantially to our current knowledge base. However, most of the theories are of a qualitative and verbal nature, and may be difficult to evaluate and compare with each other. In this paper, we propose that one way forward is to use computational modelling to formulate more precise theories of PTSD that can be evaluated by (1) assessing whether the model can explain fundamental phenomena related to PTSD, and (2) comparing simulated outcomes with real data. Computational modelling can force us to describe processes more precisely and achieve stronger theories that are viable for testing. Establishing the theoretical groundwork before undertaking empirical studies can help us to avoid doing research with low probability of valid results, and counteract the replicability crisis in psychology. In conclusion, computational modelling is a promising avenue for advancing the understanding and treatment of PTSD.

Computational modelling can help us to specify the psychological processes involved in PTSD, which may increase our understanding of how best to help people to recover after traumatic events.With computational models of PTSD, we can simulate the consequences of the theoretical principles and make sure to design research studies that are theoretically well grounded.To validate the computational models, high-quality empirical data are still needed.

A computational study of right ventricular mechanics in a rat model of pulmonary arterial hypertension.

Pulmonary arterial hypertension (PAH) presents a significant challenge to right ventricular (RV) function due to progressive pressure overload, necessitating adaptive remodeling in the form of increased wall thickness, enhanced myocardial contractility and stiffness to maintain cardiac performance. However, the impact of these remodeling mechanisms on RV mechanics in not clearly understood. In addition, there is a lack of quantitative understanding of how each mechanism individually influences RV mechanics. Utilizing experimental data from a rat model of PAH at three distinct time points, we developed biventricular finite element models to investigate how RV stress and strain evolved with PAH progression. The finite element models were fitted to hemodynamic and morphological data to represent different disease stages and used to analyze the impact of RV remodeling as well as the altered RV pressure. Furthermore, we performed a number of theoretical simulation studies with different combinations of morphological and physiological remodeling, to assess and quantify their individual impact on overall RV load and function. Our findings revealed a substantial 4-fold increase in RV stiffness and a transient 2-fold rise in contractility, which returned to baseline by week 12. These changes in RV material properties in addition to the 2-fold increase in wall thickness significantly mitigated the increase in wall stress and strain caused by the progressive increase in RV afterload. Despite the PAH-induced cases showing increased wall stress and strain at end-diastole and end-systole compared to the control, our simulations suggest that without the observed remodeling mechanisms, the increase in stress and strain would have been much more pronounced. Our model analysis also indicated that while changes in the RV's material properties-particularly increased RV stiffness - have a notable effect on its mechanics, the primary compensatory factor limiting the stress and strain increase in the early stages of PAH was the significant increase in wall thickness. These findings underscore the importance of RV remodeling in managing the mechanical burden on the right ventricle due to pressure overload.

Impact of left atrial wall motion assumptions in fluid simulations on proposed predictors of thrombus formation.

Atrial fibrillation (AF) poses a significant risk of stroke due to thrombus formation, which primarily occurs in the left atrial appendage (LAA). Medical image-based computational fluid dynamics (CFD) simulations can provide valuable insight into patient-specific hemodynamics and could potentially enhance personalized assessment of thrombus risk. However, the importance of accurately representing the left atrial (LA) wall dynamics has not been fully resolved. In this study, we compared four modeling scenarios; rigid walls, a generic wall motion based on a reference motion, a semi-generic wall motion based on patient-specific motion, and patient-specific wall motion based on medical images. We considered a LA geometry acquired from 4D computed tomography during AF, systematically performed convergence tests to assess the numerical accuracy of our solution strategy, and quantified the differences between the four approaches. The results revealed that wall motion had no discernible impact on LA cavity hemodynamics, nor on the markers that indicate thrombus formation. However, the flow patterns within the LAA deviated significantly in the rigid model, indicating that the assumption of rigid walls may lead to errors in the estimated risk factors. In contrast, the generic, semi-generic, and patient-specific cases were qualitatively similar. The results highlight the crucial role of wall motion on hemodynamics and predictors of thrombus formation, and also demonstrate the potential of using a generic motion model as a surrogate for the more complex patient-specific motion. While the present study considered a single case, the employed CFD framework is entirely open-source and designed for adaptability, allowing for integration of additional models and generic motions.

A verified and validated moving domain computational fluid dynamics solver with applications to cardiovascular flows.

Computational fluid dynamics (CFD) in combination with patient-specific medical images has been used to correlate flow phenotypes with disease initiation, progression and outcome, in search of a prospective clinical tool. A large number of CFD software packages are available, but are typically based on rigid domains and low-order finite volume methods, and are often implemented in massive low-level C++ libraries. Furthermore, only a handful of solvers have been appropriately verified and validated for their intended use. Our goal was to develop, verify and validate an open-source CFD solver for moving domains, with applications to cardiovascular flows. The solver is an extension of the CFD solver Oasis, which is based on the finite element method and implemented using the FEniCS open source framework. The new solver, named OasisMove, extends Oasis by expressing the Navier-Stokes equations in the arbitrary Lagrangian-Eulerian formulation, which is suitable for handling moving domains. For code verification we used the method of manufactured solutions for a moving 2D vortex problem, and for validation we compared our results against existing high-resolution simulations and laboratory experiments for two moving domain problems of varying complexity. Verification results showed that the L 2 error followed the theoretical convergence rates. The temporal accuracy was second-order, while the spatial accuracy was second- and third-order using ℙ 1 / ℙ 1 and ℙ 2 / ℙ 1 finite elements, respectively. Validation results showed good agreement with existing benchmark results, by reproducing lift and drag coefficients with less than 1% error, and demonstrating the solver's ability to capture vortex patterns in transitional and turbulent-like flow regimes. In conclusion, we have shown that OasisMove is an open-source, accurate and reliable solver for cardiovascular flows in moving domains.

A cell-based framework for modeling cardiac mechanics.

Cardiomyocytes are the functional building blocks of the heart-yet most models developed to simulate cardiac mechanics do not represent the individual cells and their surrounding matrix. Instead, they work on a homogenized tissue level, assuming that cellular and subcellular structures and processes scale uniformly. Here we present a mathematical and numerical framework for exploring tissue-level cardiac mechanics on a microscale given an explicit three-dimensional geometrical representation of cells embedded in a matrix. We defined a mathematical model over such a geometry and parametrized our model using publicly available data from tissue stretching and shearing experiments. We then used the model to explore mechanical differences between the extracellular and the intracellular space. Through sensitivity analysis, we found the stiffness in the extracellular matrix to be most important for the intracellular stress values under contraction. Strain and stress values were observed to follow a normal-tangential pattern concentrated along the membrane, with substantial spatial variations both under contraction and stretching. We also examined how it scales to larger size simulations, considering multicellular domains. Our work extends existing continuum models, providing a new geometrical-based framework for exploring complex cell-cell and cell-matrix interactions.

Regional Left Ventricular Fiber Stress Analysis for Cardiac Resynchronization Therapy Response.

Cardiac resynchronization therapy (CRT) is an effective treatment for a subgroup of heart failure (HF) patients, but more than 30% of those selected do not improve after CRT implantation. Imperfect pre-procedural criteria for patient selection and optimization are the main causes of the high non-response rate. In this study, we evaluated a novel measure for assessing CRT response. We used a computational modeling framework to calculate the regional stress of the left ventricular wall of seven CRT patients and seven healthy controls. The standard deviation of regional wall stress at the time of mitral valve closure (SD_MVC) was used to quantify dyssynchrony and compared between patients and controls and among the patients. The results show that SD_MVC is significantly lower in controls than patients and correlates with long-term response in patients, based on end-diastolic volume reduction. In contrast to our initial hypothesis, patients with lower SD_MVC respond better to therapy. The patient with the highest SD_MVC was the only non-responder in the patient cohort. The distribution of fiber stress at the beginning of the isovolumetric phase seems to correlate with the degree of response and the use of this measurement could potentially improve selection criteria for CRT implantation. Further studies with a larger cohort of patients are needed to validate these results.

Electromechanical feedback with reduced cellular connectivity alters electrical activity in an infarct injured left ventricle: a finite element model study.

Myocardial infarction (MI) significantly alters the structure and function of the heart. As abnormal strain may drive heart failure and the generation of arrhythmias, we used computational methods to simulate a left ventricle with an MI over the course of a heartbeat to investigate strains and their potential implications to electrophysiology. We created a fully coupled finite element model of myocardial electromechanics consisting of a cellular physiological model, a bidomain electrical diffusion solver, and a nonlinear mechanics solver. A geometric mesh built from magnetic resonance imaging (MRI) measurements of an ovine left ventricle suffering from a surgically induced anteroapical infarct was used in the model, cycled through the cardiac loop of inflation, isovolumic contraction, ejection, and isovolumic relaxation. Stretch-activated currents were added as a mechanism of mechanoelectric feedback. Elevated fiber and cross fiber strains were observed in the area immediately adjacent to the aneurysm throughout the cardiac cycle, with a more dramatic increase in cross fiber strain than fiber strain. Stretch-activated channels decreased action potential (AP) dispersion in the remote myocardium while increasing it in the border zone. Decreases in electrical connectivity dramatically increased the changes in AP dispersion. The role of cross fiber strain in MI-injured hearts should be investigated more closely, since results indicate that these are more highly elevated than fiber strain in the border of the infarct. Decreases in connectivity may play an important role in the development of altered electrophysiology in the high-stretch regions of the heart.

Computational models of ventricular mechanics and adaptation in response to right-ventricular pressure overload.

Pulmonary arterial hypertension (PAH) is associated with substantial remodeling of the right ventricle (RV), which may at first be compensatory but at a later stage becomes detrimental to RV function and patient survival. Unlike the left ventricle (LV), the RV remains understudied, and with its thin-walled crescent shape, it is often modeled simply as an appendage of the LV. Furthermore, PAH diagnosis is challenging because it often leaves the LV and systemic circulation largely unaffected. Several treatment strategies such as atrial septostomy, right ventricular assist devices (RVADs) or RV resynchronization therapy have been shown to improve RV function and the quality of life in patients with PAH. However, evidence of their long-term efficacy is limited and lung transplantation is still the most effective and curative treatment option. As such, the clinical need for improved diagnosis and treatment of PAH drives a strong need for increased understanding of drivers and mechanisms of RV growth and remodeling (G&R), and more generally for targeted research into RV mechanics pathology. Computational models stand out as a valuable supplement to experimental research, offering detailed analysis of the drivers and consequences of G&R, as well as a virtual test bench for exploring and refining hypotheses of growth mechanisms. In this review we summarize the current efforts towards understanding RV G&R processes using computational approaches such as reduced-order models, three dimensional (3D) finite element (FE) models, and G&R models. In addition to an overview of the relevant literature of RV computational models, we discuss how the models have contributed to increased scientific understanding and to potential clinical treatment of PAH patients.

Computational cardiac physiology for new modelers: Origins, foundations, and future.

Mathematical models of the cardiovascular system have come a long way since they were first introduced in the early 19th century. Driven by a rapid development of experimental techniques, numerical methods, and computer hardware, detailed models that describe physical scales from the molecular level up to organs and organ systems have been derived and used for physiological research. Mathematical and computational models can be seen as condensed and quantitative formulations of extensive physiological knowledge and are used for formulating and testing hypotheses, interpreting and directing experimental research, and have contributed substantially to our understanding of cardiovascular physiology. However, in spite of the strengths of mathematics to precisely describe complex relationships and the obvious need for the mathematical and computational models to be informed by experimental data, there still exist considerable barriers between experimental and computational physiological research. In this review, we present a historical overview of the development of mathematical and computational models in cardiovascular physiology, including the current state of the art. We further argue why a tighter integration is needed between experimental and computational scientists in physiology, and point out important obstacles and challenges that must be overcome in order to fully realize the synergy of experimental and computational physiological research.

A Computational Study of the Effects of Tachycardia-Induced Remodeling on Calcium Wave Propagation in Rabbit Atrial Myocytes.

In atrial cardiomyocytes without a well-developed T-tubule system, calcium diffuses from the periphery toward the center creating a centripetal wave pattern. During atrial fibrillation, rapid activation of atrial myocytes induces complex remodeling in diffusion properties that result in failure of calcium to propagate in a fully regenerative manner toward the center; a phenomenon termed "calcium silencing." This has been observed in rabbit atrial myocytes after exposure to prolonged rapid pacing. Although experimental studies have pointed to possible mechanisms underlying calcium silencing, their individual effects and relative importance remain largely unknown. In this study we used computational modeling of the rabbit atrial cardiomyocyte to query the individual and combined effects of the proposed mechanisms leading to calcium silencing and abnormal calcium wave propagation. We employed a population of models obtained from a newly developed model of the rabbit atrial myocyte with spatial representation of intracellular calcium handling. We selected parameters in the model that represent experimentally observed cellular remodeling which have been implicated in calcium silencing, and scaled their values in the population to match experimental observations. In particular, we changed the maximum conductances of ICaL, INCX, and INaK, RyR open probability, RyR density, Serca2a density, and calcium buffering strength. We incorporated remodeling in a population of 16 models by independently varying parameters that reproduce experimentally observed cellular remodeling, and quantified the resulting alterations in calcium dynamics and wave propagation patterns. The results show a strong effect of ICaL in driving calcium silencing, with INCX, INaK, and RyR density also resulting in calcium silencing in some models. Calcium alternans was observed in some models where INCX and Serca2a density had been changed. Simultaneously incorporating changes in all remodeled parameters resulted in calcium silencing in all models, indicating the predominant role of decreasing ICaL in the population phenotype.

Nationwide rollout reveals efficacy of epidemic control through digital contact tracing.

Fuelled by epidemiological studies of SARS-CoV-2, contact tracing by mobile phones has been put to use in many countries. Over a year into the pandemic, we lack conclusive evidence on its effectiveness. To address this gap, we used a unique real world contact data set, collected during the rollout of the first Norwegian contact tracing app in the Spring of 2020. Our dataset involves millions of contacts between 12.5% of the adult population, which enabled us to measure the real-world app performance. The technological tracing efficacy was measured at 80%, and we estimated that at least 11.0% of the discovered close contacts could not have been identified by manual contact tracing. Our results also indicated that digital contact tracing can flag individuals with excessive contacts, which can help contain superspreading related outbreaks. The overall effectiveness of digital tracing depends strongly on app uptake, but significant impact can be achieved for moderate uptake numbers. Used as a supplement to manual tracing and other measures, digital tracing can be instrumental in controlling the pandemic. Our findings can thus help informing public health policies in the coming months.

Computational Modeling Studies of the Roles of Left Ventricular Geometry, Afterload, and Muscle Contractility on Myocardial Strains in Heart Failure with Preserved Ejection Fraction.

Global longitudinal strain and circumferential strain are found to be reduced in HFpEF, which some have interpreted that the global left ventricular (LV) contractility is impaired. This finding is, however, contradicted by a preserved ejection fraction (EF) and confounded by changes in LV geometry and afterload resistance that may also affect the global strains. To reconcile these issues, we used a validated computational framework consisting of a finite element LV model to isolate the effects of HFpEF features in affecting systolic function metrics. Simulations were performed to quantify the effects on myocardial strains due to changes in LV geometry, active tension developed by the tissue, and afterload. We found that only a reduction in myocardial contractility and an increase in afterload can simultaneously reproduce the blood pressures, EF and strains measured in HFpEF patients. This finding suggests that it is likely that the myocardial contractility is reduced in HFpEF patients. Graphical abstract.

A Novel Computational Model of the Rabbit Atrial Cardiomyocyte With Spatial Calcium Dynamics.

Models of cardiac electrophysiology are widely used to supplement experimental results and to provide insight into mechanisms of cardiac function and pathology. The rabbit has been a particularly important animal model for studying mechanisms of atrial pathophysiology and atrial fibrillation, which has motivated the development of models for the rabbit atrial cardiomyocyte electrophysiology. Previously developed models include detailed representations of membrane currents and intracellular ionic concentrations, but these so-called "common-pool" models lack a spatially distributed description of the calcium handling system, which reflects the detailed ultrastructure likely found in cells in vivo. Because of the less well-developed T-tubular system in atrial compared to ventricular cardiomyocytes, spatial gradients in intracellular calcium concentrations may play a more significant role in atrial cardiomyocyte pathophysiology, rendering common-pool models less suitable for investigating underlying electrophysiological mechanisms. In this study, we developed a novel computational model of the rabbit atrial cardiomyocyte incorporating detailed compartmentalization of intracellular calcium dynamics, in addition to a description of membrane currents and intracellular processes. The spatial representation of calcium was based on dividing the intracellular space into eighteen different compartments in the transversal direction, each with separate systems for internal calcium storage and release, and tracking ionic fluxes between compartments in addition to the dynamics driven by membrane currents and calcium release. The model was parameterized employing a population-of-models approach using experimental data from different sources. The parameterization of this novel model resulted in a reduced population of models with inherent variability in calcium dynamics and electrophysiological properties, all of which fall within the range of observed experimental values. As such, the population of models may represent natural variability in cardiomyocyte electrophysiology or inherent uncertainty in the underlying experimental data. The ionic model population was also able to reproduce the U-shaped waveform observed in line-scans of triggered calcium waves in atrial cardiomyocytes, characteristic of the absence of T-tubules, resulting in a centripetal calcium wave due to subcellular calcium diffusion. This novel spatial model of the rabbit atrial cardiomyocyte can be used to integrate experimental findings, offering the potential to enhance our understanding of the pathophysiological role of calcium-handling abnormalities under diseased conditions, such as atrial fibrillation.

Multisite pacing and myocardial scars: a computational study.

Cardiac resynchronization therapy (CRT) is a frequently effective treatment modality for dyssynchronous heart failure, however, 30% of patients do not respond, usually due to suboptimal activation of the left ventricle (LV). Multisite pacing (MSP) may increase the response rate, but its effect in the presence of myocardial scars is not fully understood. We use a computational model to study the outcome of MSP in an LV with scars in two different locations and of two different sizes. The LV was stimulated from anterior, posterior and lateral locations individually and in pairs, while a septal stimulation site represented right ventricular (RV) pacing. Intraventricular pressures were measured, and outcomes evaluated in terms of maximum LV pressure gradient (dP/dtmax)- change compared to isolated RV pacing. The best result obtained using various LV pacing locations included a combination of sites remote from scars and the septum. The highest dP/dtmax increase was achieved, regardless of scar size, using MSP with one pacing site located on the LV free wall opposite to the scar and one site opposite to the septum. These in silico modelling results suggest that making placement of pacing electrodes dependent on location of scarring, may alter acute haemodynamics and that such modelling may contribute to future CRT optimization.