Search

2'-Deoxy-ATP (dATP), a naturally occurring near analog of ATP, is a well-documented myosin activator that has been shown to increase contractile force, improve pump function, and enhance lusitropy in the heart. Calcium transients in cardiomyocytes with elevated levels of dATP show faster calcium decay compared with cardiomyocytes with basal levels of dATP, but the mechanisms behind this are unknown. Here, we design and utilize a multiscale computational modeling framework to test the hypothesis that dATP acts on the sarcoendoplasmic reticulum calcium-ATPase (SERCA) pump to accelerate calcium re-uptake into the sarcoplasmic reticulum during cardiac relaxation. Gaussian accelerated molecular dynamics simulations of human cardiac SERCA2A in the E1 apo, ATP-bound and dATP-bound states showed that dATP forms more stable contacts in the nucleotide binding pocket of SERCA and leads to increased closure of cytosolic domains. These structural changes ultimately lead to changes in calcium binding, which we assessed using Brownian dynamics simulations. We found that dATP increases calcium association rate constants to SERCA and that dATP binds to apo SERCA more rapidly than ATP. Using a compartmental ordinary differential equation model of human cardiomyocyte excitation-contraction coupling, we found that these increased association rate constants contributed to the accelerated rates of calcium transient decay observed experimentally. This study provides clear mechanistic evidence of enhancements in cardiac SERCA2A pump function due to interactions with dATP.

Cardiac shape modeling is a useful computational tool that has provided quantitative insights into the mechanisms underlying dysfunction in heart disease. The manual input and time required to make cardiac shape models, however, limits their clinical utility. Here we present an end-to-end pipeline that uses deep learning for automated view classification, slice selection, phase selection, anatomical landmark localization, and myocardial image segmentation for the automated generation of three-dimensional, biventricular shape models. With this approach, we aim to make cardiac shape modeling a more robust and broadly applicable tool that has processing times consistent with clinical workflows. Cardiovascular magnetic resonance (CMR) images from a cohort of 123 patients with repaired tetralogy of Fallot (rTOF) from two internal sites were used to train and validate each step in the automated pipeline. The complete automated pipeline was tested using CMR images from a cohort of 12 rTOF patients from an internal site and 18 rTOF patients from an external site. Manually and automatically generated shape models from the test set were compared using Euclidean projection distances, global ventricular measurements, and atlas-based shape mode scores. The mean absolute error (MAE) between manually and automatically generated shape models in the test set was similar to the voxel resolution of the original CMR images for end-diastolic models (MAE = 1.9 ± 0.5 mm) and end-systolic models (MAE = 2.1 ± 0.7 mm). Global ventricular measurements computed from automated models were in good agreement with those computed from manual models. The average mean absolute difference in shape mode Z-score between manually and automatically generated models was 0.5 standard deviations for the first 20 modes of a reference statistical shape atlas. Using deep learning, accurate three-dimensional, biventricular shape models can be reliably created. This fully automated end-to-end approach dramatically reduces the manual input required to create shape models, thereby enabling the rapid analysis of large-scale datasets and the potential to deploy statistical atlas-based analyses in point-of-care clinical settings. Training data and networks are available from cardiacatlas.org.

Current indications for pulmonary valve replacement (PVR) in repaired tetralogy of Fallot (rTOF) rely on cardiovascular magnetic resonance (CMR) image-based indices but are inconsistently applied, lead to mixed outcomes, and remain debated. This study aimed to test the hypothesis that specific markers of biventricular shape may discriminate differences between rTOF patients who did and did not require subsequent PVR better than standard imaging indices. In this cross-sectional retrospective study, biventricular shape models were customized to CMR images from 84 rTOF patients. A statistical atlas of end-diastolic shape was constructed using principal component analysis. Multivariate regression was used to quantify shape mode and imaging index associations with subsequent intervention status (PVR, n = 48 vs. No-PVR, n = 36), while accounting for confounders. Clustering analysis was used to test the ability of the most significant shape modes and imaging indices to discriminate PVR status as evaluated by a Matthews correlation coefficient (MCC). Geometric strain analysis was also conducted to assess shape mode associations with systolic function. PVR status correlated significantly with shape modes associated with right ventricular (RV) apical dilation and left ventricular (LV) dilation (p

Electrical impulse generation and its conduction within cells or cellular networks are the cornerstone of electrophysiology. However, the advancement of the field is limited by sensing accuracy and the scalability of current recording technologies. Here we describe a scalable platform that enables accurate recording of transmembrane potentials in electrogenic cells. The platform employs a three-dimensional high-performance field-effect transistor array for minimally invasive cellular interfacing that produces faithful recordings, as validated by the gold standard patch clamp. Leveraging the high spatial and temporal resolutions of the field-effect transistors, we measured the intracellular signal conduction velocity of a cardiomyocyte to be 0.182 m s-1, which is about five times the intercellular velocity. We also demonstrate intracellular recordings in cardiac muscle tissue constructs and reveal the signal conduction paths. This platform could provide new capabilities in probing the electrical behaviours of single cells and cellular networks, which carries broad implications for understanding cellular physiology, pathology and cell-cell interactions.

Dilated cardiomyopathy (DCM) is a life-threatening form of heart disease that is typically characterized by progressive thinning of the ventricular walls, chamber dilation, and systolic dysfunction. Multiple mutations in the gene encoding filamin C (FLNC), an actin-binding cytoskeletal protein in cardiomyocytes, have been found in patients with DCM. However, the mechanisms that lead to contractile impairment and DCM in patients with FLNC variants are poorly understood. To determine how FLNC regulates systolic force transmission and DCM remodeling, we used an inducible, cardiac-specific FLNC-knockout (icKO) model to produce a rapid onset of DCM in adult mice. Loss of FLNC reduced systolic force development in single cardiomyocytes and isolated papillary muscles but did not affect twitch kinetics or calcium transients. Electron and immunofluorescence microscopy showed significant defects in Z-disk alignment in icKO mice and altered myofilament lattice geometry. Moreover, a loss of FLNC induces a softening myocyte cortex and structural adaptations at the subcellular level that contribute to disrupted longitudinal force production during contraction. Spatially explicit computational models showed that these structural defects could be explained by a loss of inter-myofibril elastic coupling at the Z-disk. Our work identifies FLNC as a key regulator of the multiscale ultrastructure of cardiomyocytes and therefore plays an important role in maintaining systolic mechanotransmission pathways, the dysfunction of which may be key in driving progressive DCM.

Misregulated gene expression in human hearts can result in cardiovascular diseases that are leading causes of mortality worldwide. However, the limited information on the genomic location of candidate cis-regulatory elements (cCREs) such as enhancers and promoters in distinct cardiac cell types has restricted the understanding of these diseases. Here, we defined >287,000 cCREs in the four chambers of the human heart at single-cell resolution, which revealed cCREs and candidate transcription factors associated with cardiac cell types in a region-dependent manner and during heart failure. We further found cardiovascular disease-associated genetic variants enriched within these cCREs including 38 candidate causal atrial fibrillation variants localized to cardiomyocyte cCREs. Additional functional studies revealed that two of these variants affect a cCRE controlling KCNH2/HERG expression and action potential repolarization. Overall, this atlas of human cardiac cCREs provides the foundation for illuminating cell type-specific gene regulation in human hearts during health and disease.

AimsVentricular activation patterns can aid clinical decision-making directly by providing spatial information on cardiac electrical activation or indirectly through derived clinical indices. The aim of this work was to derive an atlas of the major modes of variation of ventricular activation from model-predicted 3D bi-ventricular activation time distributions and to relate these modes to corresponding vectorcardiograms (VCGs). We investigated how the resulting dimensionality reduction can improve and accelerate the estimation of activation patterns from surface electrogram measurements.Methods and resultsAtlases of activation time (AT) and VCGs were derived using principal component analysis on a dataset of simulated electrophysiology simulations computed on eight patient-specific bi-ventricular geometries. The atlases provided significant dimensionality reduction, and the modes of variation in the two atlases described similar features. Utility of the atlases was assessed by resolving clinical waveforms against them and the VCG atlas was able to accurately reconstruct the patient VCGs with fewer than 10 modes. A sensitivity analysis between the two atlases was performed by calculating a compact Jacobian. Finally, VCGs generated by varying AT atlas modes were compared with clinical VCGs to estimate patient-specific activation maps, and the resulting errors between the clinical and atlas-based VCGs were less than those from more computationally expensive method.ConclusionAtlases of activation and VCGs represent a new method of identifying and relating the features of these high-dimensional signals that capture the major sources of variation between patients and may aid in identifying novel clinical indices of arrhythmia risk or therapeutic outcome.

This document, published in the Journal of Applied Physics, discusses the importance of multiscale modeling and analysis in the field of biophysics. The authors emphasize the need for tools and analyses that integrate across physical scales and time scales in order to advance medical science. The document includes a collection of studies on various topics such as cardiac biomechanics, tissue mechanics, cerebral fluid flow, cell mechanics, protein modeling, and molecular fluorescence microscopy. The authors hope that readers will find these studies informative and inspiring. [Extracted from the article]

BackgroundPressure overload left ventricular (LV) hypertrophy is characterized by increased cardiomyocyte width and ventricle wall thickness, however the regional variation of this remodeling is unclear. Cardiovascular magnetic resonance (CMR) diffusion tensor imaging (DTI) may provide a non-invasive, comprehensive, and geometrically accurate method to detect regional differences in structural remodeling in hypertrophy. We hypothesized that DTI parameters, such as fractional and planar anisotropy, would reflect myocyte remodeling due to pressure overload in a regionally-dependent manner.MethodsWe investigated the regional distributions of myocyte remodeling in rats with or without transverse aortic constriction (TAC) via direct measurement of myocyte dimensions with confocal imaging of thick tissue sections, and correlated myocyte cross-sectional area and other geometric features with parameters of diffusivity from ex-vivo DTI in the same regions of the same hearts.ResultsWe observed regional differences in several parameters from DTI between TAC hearts and SHAM controls. Consistent with previous studies, helix angles from DTI correlated strongly with those measured directly from histological sections (p

BackgroundAdverse cardiac remodeling is an important precursor to anthracycline-related cardiac dysfunction, however conventional remodeling indices are limited. We sought to examine the utility of statistical atlas-derived measures of ventricular shape to improve the identification of adverse anthracycline-related remodeling in childhood cancer survivors.MethodsWe analyzed cardiac magnetic resonance imaging from a cross-sectional cohort of 20 childhood cancer survivors who were treated with low (

2'-deoxy-ATP (dATP) is a naturally occurring small molecule that has shown promise as a therapeutic because it significantly increases cardiac myocyte force development even at low dATP/ATP ratios. To investigate mechanisms by which dATP alters myosin crossbridge dynamics, we used Brownian dynamics simulations to calculate association rates between actin and ADP- or dADP-bound myosin. These rates were then directly incorporated in a mechanistic Monte Carlo Markov Chain model of cooperative sarcomere contraction. A unique combination of increased powerstroke and detachment rates was required to match experimental steady-state and kinetic data for dATP force production in rat cardiac myocytes when the myosin attachment rate in the model was constrained by the results of a Brownian dynamics simulation. Nearest-neighbor cooperativity was seen to contribute to, but not fully explain, the steep relationship between dATP/ATP ratio and steady-state force-development observed at lower dATP concentrations. Dynamic twitch simulations performed using measured calcium transients as inputs showed that the effects of dATP on the crossbridge alone were not sufficient to explain experimentally observed enhancement of relaxation kinetics by dATP treatment. Hence, dATP may also affect calcium handling even at low concentrations. By enabling the effects of dATP on sarcomere mechanics to be predicted, this multi-scale modeling framework may elucidate the molecular mechanisms by which dATP can have therapeutic effects on cardiac contractile dysfunction.

Dilated cardiomyopathy (DCM) is often associated with sarcomere protein mutations that confer reduced myofilament tension-generating capacity. We demonstrated that cardiac twitch tension-time integrals can be targeted and tuned to prevent DCM remodeling in hearts with contractile dysfunction. We employed a transgenic murine model of DCM caused by the D230N-tropomyosin (Tm) mutation and designed a sarcomere-based intervention specifically targeting the twitch tension-time integral of D230N-Tm hearts using multiscale computational models of intramolecular and intermolecular interactions in the thin filament and cell-level contractile simulations. Our models predicted that increasing the calcium sensitivity of thin filament activation using the cardiac troponin C (cTnC) variant L48Q can sufficiently augment twitch tension-time integrals of D230N-Tm hearts. Indeed, cardiac muscle isolated from double-transgenic hearts expressing D230N-Tm and L48Q cTnC had increased calcium sensitivity of tension development and increased twitch tension-time integrals compared with preparations from hearts with D230N-Tm alone. Longitudinal echocardiographic measurements revealed that DTG hearts retained normal cardiac morphology and function, whereas D230N-Tm hearts developed progressive DCM. We present a computational and experimental framework for targeting molecular mechanisms governing the twitch tension of cardiomyopathic hearts to counteract putative mechanical drivers of adverse remodeling and open possibilities for tension-based treatments of genetic cardiomyopathies.

RationaleZO-1 (Zonula occludens-1), a plasma membrane-associated scaffolding protein regulates signal transduction, transcription, and cellular communication. Global deletion of ZO-1 in the mouse is lethal by embryonic day 11.5. The function of ZO-1 in cardiac myocytes (CM) is largely unknown.ObjectiveTo determine the function of CM ZO-1 in the intact heart, given its binding to other CM proteins that have been shown instrumental in normal cardiac conduction and function.Methods and resultsWe generated ZO-1 CM-specific knockout (KO) mice using α-Myosin Heavy Chain-nuclear Cre (ZO-1cKO) and investigated physiological and electrophysiological function by echocardiography, surface ECG and conscious telemetry, intracardiac electrograms and pacing, and optical mapping studies. ZO-1cKO mice were viable, had normal Mendelian ratios, and had a normal lifespan. Ventricular morphometry and function were not significantly different between the ZO-1cKO versus control (CTL) mice, basally in young or aged mice, or even when hearts were subjected to hemodynamic loading. Atrial mass was increased in ZO-1cKO. Electrophysiological and optical mapping studies indicated high-grade atrioventricular (A-V) block in ZO-1cKO comparing to CTL hearts. While ZO-1-associated proteins such as vinculin, connexin 43, N-cadherin, and α-catenin showed no significant change with the loss of ZO-1, Connexin-45 and Coxsackie-adenovirus (CAR) proteins were reduced in atria of ZO-1cKO. Further, with loss of ZO-1, ZO-2 protein was increased significantly in ventricular CM in a presumed compensatory manner but was still not detected in the AV nodal myocytes. Importantly, the expression of the sodium channel protein NaV1.5 was altered in AV nodal cells of the ZO-1cKO versus CTL.ConclusionsZO-1 protein has a unique physiological role in cardiac nodal tissue. This is in alignment with its known interaction with CAR and Cx45, and a new function in regulating the expression of NaV1.5 in AV node. Uniquely, ZO-1 is dispensable for function of the working myocardium.

Cardiac myocytes transduce changes in mechanical loading into cellular responses via interacting cell signalling pathways. We previously reported a logic-based ordinary differential equation model of the myocyte mechanosignalling network that correctly predicts 78% of independent experimental results not used to formulate the original model. Here, we use Monte Carlo and polynomial chaos expansion simulations to examine the effects of uncertainty in parameter values, model logic and experimental validation data on the assessed accuracy of that model. The prediction accuracy of the model was robust to parameter changes over a wide range being least sensitive to uncertainty in time constants and most affected by uncertainty in reaction weights. Quantifying epistemic uncertainty in the reaction logic of the model showed that while replacing 'OR' with 'AND' reactions greatly reduced model accuracy, replacing 'AND' with 'OR' reactions was more likely to maintain or even improve accuracy. Finally, data uncertainty had a modest effect on assessment of model accuracy. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

We use Brownian-Langevin dynamics principles to derive a coarse-graining multiscale myofilament model that can describe the thin-filament activation process during contraction. The model links atomistic molecular simulations of protein-protein interactions in the thin-filament regulatory unit to sarcomere-level activation dynamics. We first calculate the molecular interaction energy between tropomyosin and actin surface using Brownian dynamics simulations. This energy profile is then generalized to account for the observed tropomyosin transitions between its regulatory stable states. The generalized energy landscape then served as a basis for developing a filament-scale model using Langevin dynamics. This integrated analysis, spanning molecular to thin-filament scales, is capable of tracking the events of the tropomyosin conformational changes as it moves over the actin surface. The tropomyosin coil with flexible overlap regions between adjacent tropomyosins is represented in the model as a system of coupled stochastic ordinary differential equations. The proposed multiscale approach provides a more detailed molecular connection between tropomyosin dynamics, the trompomyosin-actin interaction-energy landscape, and the generated force by the sarcomere.

The naturally occurring nucleotide 2-deoxy-adenosine 5'-triphosphate (dATP) can be used by cardiac muscle as an alternative energy substrate for myosin chemomechanical activity. We and others have previously shown that dATP increases contractile force in normal hearts and models of depressed systolic function, but the structural basis of these effects has remained unresolved. In this work, we combine multiple techniques to provide structural and functional information at the angstrom-nanometer and millisecond time scales, demonstrating the ability to make both structural measurements and quantitative kinetic estimates of weak actin-myosin interactions that underpin sarcomere dynamics. Exploiting dATP as a molecular probe, we assess how small changes in myosin structure translate to electrostatic-based changes in sarcomere function to augment contractility in cardiac muscle. Through Brownian dynamics simulation and computational structural analysis, we found that deoxy-hydrolysis products [2-deoxy-adenosine 5'-diphosphate (dADP) and inorganic phosphate (Pi)] bound to prepowerstroke myosin induce an allosteric restructuring of the actin-binding surface on myosin to increase the rate of cross-bridge formation. We then show experimentally that this predicted effect translates into increased electrostatic interactions between actin and cardiac myosin in vitro. Finally, using small-angle X-ray diffraction analysis of sarcomere structure, we demonstrate that the proposed increased electrostatic affinity of myosin for actin causes a disruption of the resting conformation of myosin motors, resulting in their repositioning toward the thin filament before activation. The dATP-mediated structural alterations in myosin reported here may provide insight into an improved criterion for the design or selection of small molecules to be developed as therapeutic agents to treat systolic dysfunction.

Nanoscale multipoint structure-function analysis is essential for deciphering the complexity of multiscale biological and physical systems. Atomic force microscopy (AFM) allows nanoscale structure-function imaging in various operating environments and can be integrated seamlessly with disparate probe-based sensing and manipulation technologies. Conventional AFMs only permit sequential single-point analysis; widespread adoption of array AFMs for simultaneous multipoint study is challenging owing to the intrinsic limitations of existing technological approaches. Here, we describe a prototype dispersive optics-based array AFM capable of simultaneously monitoring multiple probe-sample interactions. A single supercontinuum laser beam is utilized to spatially and spectrally map multiple cantilevers, to isolate and record beam deflection from individual cantilevers using distinct wavelength selection. This design provides a remarkably simplified yet effective solution to overcome the optical cross-talk while maintaining subnanometer sensitivity and compatibility with probe-based sensors. We demonstrate the versatility and robustness of our system on parallel multiparametric imaging at multiscale levels ranging from surface morphology to hydrophobicity and electric potential mapping in both air and liquid, mechanical wave propagation in polymeric films, and the dynamics of living cells. This multiparametric, multiscale approach provides opportunities for studying the emergent properties of atomic-scale mechanical and physicochemical interactions in a wide range of physical and biological networks.

Multi-scale computational modeling is a major branch of computational biology as evidenced by the US federal interagency Multi-Scale Modeling Consortium and major international projects. It invariably involves specific and detailed sequences of data analysis and simulation, often with multiple tools and datasets, and the community recognizes improved modularity, reuse, reproducibility, portability and scalability as critical unmet needs in this area. Scientific workflows are a well-recognized strategy for addressing these needs in scientific computing. While there are good examples if the use of scientific workflows in bioinformatics, medical informatics, biomedical imaging and data analysis, there are fewer examples in multi-scale computational modeling in general and cardiac electrophysiology in particular. Cardiac electrophysiology simulation is a mature area of multi-scale computational biology that serves as an excellent use case for developing and testing new scientific workflows. In this article, we develop, describe and test a computational workflow that serves as a proof of concept of a platform for the robust integration and implementation of a reusable and reproducible multi-scale cardiac cell and tissue model that is expandable, modular and portable. The workflow described leverages Python and Kepler-Python actor for plotting and pre/post-processing. During all stages of the workflow design, we rely on freely available open-source tools, to make our workflow freely usable by scientists.

In patients with dyssynchronous heart failure (DHF), cardiac conduction abnormalities cause the regional distribution of myocardial work to be non-homogeneous. Cardiac resynchronization therapy (CRT) using an implantable, programmed biventricular pacemaker/defibrillator, can improve the synchrony of contraction between the right and left ventricles in DHF, resulting in reduced morbidity and mortality and increased quality of life. Since regional work depends on wall stress, which cannot be measured in patients, we used computational methods to investigate regional work distributions and their changes after CRT. We used three-dimensional multi-scale patient-specific computational models parameterized by anatomic, functional, hemodynamic, and electrophysiological measurements in eight patients with heart failure and left bundle branch block (LBBB) who received CRT. To increase clinical translatability, we also explored whether streamlined computational methods provide accurate estimates of regional myocardial work. We found that CRT increased global myocardial work efficiency with significant improvements in non-responders. Reverse ventricular remodeling after CRT was greatest in patients with the highest heterogeneity of regional work at baseline, however the efficacy of CRT was not related to the decrease in overall work heterogeneity or to the reduction in late-activated regions of high myocardial work. Rather, decreases in early-activated regions of myocardium performing negative myocardial work following CRT best explained patient variations in reverse remodeling. These findings were also observed when regional myocardial work was estimated using ventricular pressure as a surrogate for myocardial stress and changes in endocardial surface area as a surrogate for strain. These new findings suggest that CRT promotes reverse ventricular remodeling in human dyssynchronous heart failure by increasing regional myocardial work in early-activated regions of the ventricles, where dyssynchrony is specifically associated with hypoperfusion, late systolic stretch, and altered metabolic activity and that measurement of these changes can be performed using streamlined approaches. Author summary: Cardiac resynchronization therapy (CRT) can improve the synchrony of contraction between the right and left ventricles in DHF and reduce morbidity and mortality and increase quality of life. However, regional work before CRT and changes with CRT are not well understood. We used three-dimensional multi-scale patient-specific computational models parameterized by anatomic, functional, hemodynamic, and electrophysiological measurements to assess changes in eight patients with heart failure and left bundle branch block (LBBB) who received CRT. To increase clinical translatability, we also explored the utility of streamlined methods to estimate regional myocardial work. We found that CRT increased global myocardial work efficiency with significant improvements in non-responders. Reverse ventricular remodeling after CRT was greatest in patients with the highest heterogeneity of regional work at baseline. Decreases in early-activated regions of myocardium performing negative myocardial work with CRT best explained patient variations in reverse remodeling. These findings were also observed when regional myocardial work was estimated using ventricular pressure as a surrogate for myocardial stress and changes in endocardial surface area as a surrogate for strain. These findings suggest that CRT promotes reverse ventricular remodeling in human dyssynchronous heart failure by increasing regional myocardial work in early-activated regions of the ventricles and that measurement of these changes can be performed using streamlined approaches. [ABSTRACT FROM AUTHOR]

Patient-specific models of the ventricular myocardium, combined with the computational power to run rapid simulations, are approaching the level where they could be used for personalized cardiovascular medicine. A major remaining challenge is determining model parameters from available patient data, especially for models of the Purkinje-myocardial junctions (PMJs): the sites of initial ventricular electrical activation. There are no non-invasive methods for localizing PMJs in patients, and the relationship between the standard clinical ECG and PMJ model parameters is underexplored. Thus, this study aimed to determine the sensitivity of the QRS complex of the ECG to the anatomical location and regional number of PMJs. The QRS complex was simulated using an image-based human torso and biventricular model, and cardiac electrophysiology was simulated using Cardioid. The PMJs were modeled as discrete current injection stimuli, and the location and number of stimuli were varied within initial activation regions based on published experiments. Results indicate that the QRS complex features were most sensitive to the presence or absence of four "seed" stimuli, and adjusting locations of nearby "regional" stimuli provided finer tuning. Decreasing number of regional stimuli by an order of magnitude resulted in virtually no change in the QRS complex. Thus, a minimal 12-stimuli configuration was identified that resulted in physiological excitation, defined by QRS complex feature metrics and ventricular excitation pattern. Overall, the sensitivity results suggest that parameterizing PMJ location, rather than number, be given significantly higher priority in future studies creating personalized ventricular models from patient-derived ECGs.

Ventricular fibrillation is a common life-threatening arrhythmia. The ECG of VF appears chaotic but may allow identification of sustaining mechanisms to guide therapy.We hypothesized that rotors and focal sources manifest distinct features on the ECG, and computational modeling may identify mechanisms of such features.VF induction was attempted in 31 patients referred for ventricular arrhythmia ablation. Simultaneous surface ECG and intracardiac electrograms were recorded using biventricular basket catheters. Endocardial phase maps were used to mechanistically classify each VF cycle as rotor or focally driven. ECGs were analyzed from patients demonstrating both mechanisms in the primary analysis and from all patients with induced VF in the secondary analysis. The ECG voltage variation during each mechanism was compared. Biventricular computer simulations of VF driven by focal sources or rotors were created and resulting ECGs of each VF mechanism were compared.Rotor-based VF exhibited greater voltage variation than focal source-based VF in both the primary analysis (n = 8, 110 ± 24% vs. 55 ± 32%, P = 0.02) and the secondary analysis (n = 18, 103 ± 30% vs. 67 ± 34%, P = 0.009). Computational VF simulations also revealed greater voltage variation in rotors compared to focal sources (110 ± 19% vs. 33 ± 16%, P = 0.001), and demonstrated that this variation was due to wavebreak, secondary rotor initiation, and rotor meander.Clinical and computational studies reveal that quantitative criteria of ECG voltage variation differ significantly between VF-sustaining rotors and focal sources, and provide insight into the mechanisms of such variation. Future studies should prospectively evaluate if these criteria can separate clinical VF mechanisms and guide therapy.

Cardiac fibrosis is a serious condition currently lacking effective treatments. It occurs as a result of cardiac fibroblast (CFB) activation and differentiation into myofibroblasts, characterized by proliferation, extracellular matrix (ECM) production and stiffening, and contraction due to the expression of smooth muscle α-actin. The mechanical properties of myocardium change regionally and over time after myocardial infarction (MI). Although mechanical cues are known to activate CFBs, it is unclear which specific mechanical stimuli regulate which specific phenotypic trait; thus we investigated these relationships using three in vitro models of CFB mechanical activation and found that 1) paracrine signaling from stretched cardiomyocytes induces CFB proliferation under mechanical conditions similar to those of the infarct border region; 2) direct stretch of CFBs mimicking the mechanical environment of the infarct region induces a synthetic phenotype with elevated ECM production; and 3) progressive matrix stiffening, modeling the mechanical effects of infarct scar maturation, causes smooth muscle α-actin fiber formation, up-regulation of collagen I, and down-regulation of collagen III. These findings suggest that myocyte stretch, fibroblast stretch, and matrix stiffening following MI may separately regulate different profibrotic traits of activated CFBs.

Cardiac fibrosis, the excessive accumulation of extracellular matrix (ECM), remains an unresolved problem in most forms of heart disease. In order to be successful in preventing, attenuating or reversing cardiac fibrosis, it is essential to understand the processes leading to ECM production and accumulation. Cardiac fibroblasts are the main producers of cardiac ECM, and harbor great phenotypic plasticity. They are activated by the disease-associated changes in mechanical properties of the heart, including stretch and increased tissue stiffness. Despite much remaining unknown, an interesting body of evidence exists on how mechanical forces are translated into transcriptional responses important for determination of fibroblast phenotype and production of ECM constituents. Such mechanotransduction can occur at multiple cellular locations including the plasma membrane, cytoskeleton and nucleus. Moreover, the ECM functions as a reservoir of pro-fibrotic signaling molecules that can be released upon mechanical stress. We here review the current status of knowledge of mechanotransduction signaling pathways in cardiac fibroblasts that culminate in pro-fibrotic gene expression.

2'-deoxy-ATP (dATP) improves cardiac function by increasing the rate of crossbridge cycling and Ca2+ transient decay. However, the mechanisms of these effects and how therapeutic responses to dATP are achieved when dATP is only a small fraction of the total ATP pool remain poorly understood. Here, we used a multiscale computational modeling approach to analyze the mechanisms by which dATP improves ventricular function. We integrated atomistic simulations of prepowerstroke myosin and actomyosin association, filament-scale Markov state modeling of sarcomere mechanics, cell-scale analysis of myocyte Ca2+ dynamics and contraction, organ-scale modeling of biventricular mechanoenergetics, and systems level modeling of circulatory dynamics. Molecular and Brownian dynamics simulations showed that dATP increases the actomyosin association rate by 1.9 fold. Markov state models predicted that dATP increases the pool of myosin heads available for crossbridge cycling, increasing steadystate force development at low dATP fractions by 1.3 fold due to mechanosensing and nearest-neighbor cooperativity. This was found to be the dominant mechanism by which small amounts of dATP can improve contractile function at myofilament to organ scales. Together with faster myocyte Ca2+ handling, this led to improved ventricular contractility, especially in a failing heart model in which dATP increased ejection fraction by 16% and the energy efficiency of cardiac contraction by 1%. This work represents a complete multiscale model analysis of a small molecule myosin modulator from single molecule to organ system biophysics and elucidates how the molecular mechanisms of dATP may improve cardiovascular function in heart failure with reduced ejection fraction. [ABSTRACT FROM AUTHOR]

Exercise‐induced muscle adaptations vary based on exercise modality and intensity. We constructed a signalling network model from 87 published studies of human or rodent skeletal muscle cell responses to endurance or resistance exercise in vivo or simulated exercise in vitro. The network comprises 259 signalling interactions between 120 nodes, representing eight membrane receptors and eight canonical signalling pathways regulating 14 transcriptional regulators, 28 target genes and 12 exercise‐induced phenotypes. Using this network, we formulated a logic‐based ordinary differential equation model predicting time‐dependent molecular and phenotypic alterations following acute endurance and resistance exercises. Compared with nine independent studies, the model accurately predicted 18/21 (85%) acute responses to resistance exercise and 12/16 (75%) acute responses to endurance exercise. Detailed sensitivity analysis of differential phenotypic responses to resistance and endurance training showed that, in the model, exercise regulates cell growth and protein synthesis primarily by signalling via mechanistic target of rapamycin, which is activated by Akt and inhibited in endurance exercise by AMP‐activated protein kinase. Endurance exercise preferentially activates inflammation via reactive oxygen species and nuclear factor κB signalling. Furthermore, the expected preferential activation of mitochondrial biogenesis by endurance exercise was counterbalanced in the model by protein kinase C in response to resistance training. This model provides a new tool for investigating cross‐talk between skeletal muscle signalling pathways activated by endurance and resistance exercise, and the mechanisms of interactions such as the interference effects of endurance training on resistance exercise outcomes. What is the central question of this study?How do the cell signalling pathways regulating skeletal myocyte responses to resistance and endurance exercise interact?What is the main finding and its importance?A new systems model of skeletal muscle signalling pathways activated by resistance and endurance training was developed with eight canonical signalling pathways regulating 14 transcriptional regulators, 28 target genes and 12 exercise‐induced phenotypes. The model accurately predicted 85% of independent measurements for resistance exercise and 75% for endurance training. Analysis revealed pathways regulating the preferential activation of protein synthesis and cell growth by resistance training and inflammation by endurance exercise. [ABSTRACT FROM AUTHOR]

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.

AimsCurrent mechanisms driving cardiac pacemaker function have focused on ion channel and gap junction channel function, which are essential for action potential generation and propagation between pacemaker cells. However, pacemaker cells also harbour desmosomes that structurally anchor pacemaker cells to each other in tissue, but their role in pacemaker function remains unknown.Methods and resultsTo determine the role of desmosomes in pacemaker function, we generated a novel mouse model harbouring cardiac conduction-specific ablation (csKO) of the central desmosomal protein, desmoplakin (DSP) using the Hcn4-Cre-ERT2 mouse line. Hcn4-Cre targets cells of the adult mouse sinoatrial node (SAN) and can ablate DSP expression in the adult DSP csKO SAN resulting in specific loss of desmosomal proteins and structures. Dysregulation of DSP via loss-of-function (adult DSP csKO mice) and mutation (clinical case of a patient harbouring a pathogenic DSP variant) in mice and man, respectively, revealed that desmosomal dysregulation is associated with a primary phenotype of increased sinus pauses/dysfunction in the absence of cardiomyopathy. Underlying defects in beat-to-beat regulation were also observed in DSP csKO mice in vivo and intact atria ex vivo. DSP csKO SAN exhibited migrating lead pacemaker sites associated with connexin 45 loss. In vitro studies exploiting ventricular cardiomyocytes that harbour DSP loss and concurrent early connexin loss phenocopied the loss of beat-to-beat regulation observed in DSP csKO mice and atria, extending the importance of DSP-associated mechanisms in driving beat-to-beat regulation of working cardiomyocytes.ConclusionWe provide evidence of a mechanism that implicates an essential role for desmosomes in cardiac pacemaker function, which has broad implications in better understanding mechanisms underlying beat-to-beat regulation as well as sinus node disease and dysfunction.

Subcellular compartmentation of the ubiquitous second messenger cAMP has been widely proposed as a mechanism to explain unique receptor-dependent functional responses. How exactly compartmentation is achieved, however, has remained a mystery for more than 40 years. In this study, we developed computational and mathematical models to represent a subcellular sarcomeric space in a cardiac myocyte with varying detail. We then used these models to predict the contributions of various mechanisms that establish subcellular cAMP microdomains. We used the models to test the hypothesis that phosphodiesterases act as functional barriers to diffusion, creating discrete cAMP signaling domains. We also used the models to predict the effect of a range of experimentally measured diffusion rates on cAMP compartmentation. Finally, we modeled the anatomical structures in a cardiac myocyte diad, to predict the effects of anatomical diffusion barriers on cAMP compartmentation. When we incorporated experimentally informed model parameters to reconstruct an in silico subcellular sarcomeric space with spatially distinct cAMP production sites linked to caveloar domains, the models predict that under realistic conditions phosphodiesterases alone were insufficient to generate significant cAMP gradients. This prediction persisted even when combined with slow cAMP diffusion. When we additionally considered the effects of anatomic barriers to diffusion that are expected in the cardiac myocyte dyadic space, cAMP compartmentation did occur, but only when diffusion was slow. Our model simulations suggest that additional mechanisms likely contribute to cAMP gradients occurring in submicroscopic domains. The difference between the physiological and pathological effects resulting from the production of cAMP may be a function of appropriate compartmentation of cAMP signaling. Therefore, understanding the contribution of factors that are responsible for coordinating the spatial and temporal distribution of cAMP at the subcellular level could be important for developing new strategies for the prevention or treatment of unfavorable responses associated with different disease states.

BackgroundWe have previously shown that non-equilibrium Na(+) current (INa) reactivation drives isoproterenol-induced phase-3 early afterdepolarizations (EADs) in mouse ventricular myocytes. In these cells, EAD initiation occurs secondary to potentiated sarcoplasmic reticulum Ca(2+) release and enhanced Na(+)/Ca(2+) exchange (NCX). This can be abolished by tetrodotoxin-blockade of INa, but not ranolazine, which selectively inhibits ventricular late INa.AimSince repolarization of human atrial myocytes is similar to mouse ventricular myocytes in that it is relatively rapid and potently modulated by Ca(2+), we investigated whether similar mechanisms can evoke EADs in human atrium. Indeed, phase-3 EADs have been shown to re-initiate atrial fibrillation (AF) during autonomic stimulation, which is a well-recognized initiator of AF.MethodsWe integrated a Markov model of INa gating in our human atrial myocyte model. To simulate experimental results, we rapidly paced this cell model at 10Hz in the presence of 0.1μM acetylcholine and 1μM isoproterenol, and assessed EAD occurrence upon return to sinus rhythm (1Hz).ResultsCellular Ca(2+) loading during fast pacing results in a transient period of hypercontractility after return to sinus rhythm. Here, fast repolarization and enhanced NCX facilitate INa reactivation via the canonical gating mode (i.e., not late INa burst mode), which drives EAD initiation. Simulating ranolazine administration reduces atrial peak INa and leads to faster repolarization, during which INa fails to reactivate and EADs are prevented.ConclusionsNon-equilibrium INa reactivation can critically contribute to arrhythmias, specifically in human atrial myocytes. Ranolazine might be beneficial in this context by blocking peak (not late) atrial INa.

Caveolin-3 (Cav-3) plays a critical role in organizing signaling molecules and ion channels involved in cardiac conduction and metabolism. Mutations in Cav-3 are implicated in cardiac conduction abnormalities and myopathies. Additionally, cardiac-specific overexpression of Cav-3 (Cav-3 OE) is protective against ischemic and hypertensive injury, suggesting a potential role for Cav-3 in basal cardiac electrophysiology and metabolism involved in stress adaptation. We hypothesized that overexpression of Cav-3 may alter baseline cardiac conduction and metabolism. We examined: (1) ECG telemetry recordings at baseline and during pharmacological interventions, (2) ion channels involved in cardiac conduction with immunoblotting and computational modeling, and (3) baseline metabolism in Cav-3 OE and transgene-negative littermate control mice. Cav-3 OE mice had decreased heart rates, prolonged PR intervals, and shortened QTc intervals with no difference in activity compared to control mice. Dobutamine or propranolol did not cause significant changes between experimental groups in maximal (dobutamine) or minimal (propranolol) heart rate. Cav-3 OE mice had an overall lower chronotropic response to atropine. The expression of Kv1.4 and Kv4.3 channels, Nav1.5 channels, and connexin 43 were increased in Cav-3 OE mice. A computational model integrating the immunoblotting results indicated shortened action potential duration in Cav-3 OE mice linking the change in channel expression to the observed electrophysiology phenotype. Metabolic profiling showed no gross differences in VO2, VCO2, respiratory exchange ratio, heat generation, and feeding or drinking. In conclusion, Cav-3 OE mice have changes in ECG intervals, heart rates, and cardiac ion channel expression. These findings give novel mechanistic insights into previously reported Cav-3 dependent cardioprotection.

The molecular pharmacology of the G protein-coupled receptors for sphingosine 1-phosphate (S1P) provides important insight into established and new therapeutic targets. A new, potent bitopic S1P3 antagonist, SPM-354, with in vivo activity, has been used, together with S1P3-knockin and S1P3-knockout mice to define the spatial and functional properties of S1P3 in regulating cardiac conduction. We show that S1P3 is a key direct regulator of cardiac rhythm both in vivo and in isolated perfused hearts. 2-Amino-2-[2-(4-octylphenyl)ethyl]propane-1,3-diol in vivo and S1P in isolated hearts induced a spectrum of cardiac effects, ranging from sinus bradycardia to complete heart block, as measured by a surface electrocardiogram in anesthetized mice and in volume-conducted Langendorff preparations. The agonist effects on complete heart block are absent in S1P3-knockout mice and are reversed in wild-type mice with SPM-354, as characterized and described here. Homologous knockin of S1P3-mCherry is fully functional pharmacologically and is strongly expressed by immunohistochemistry confocal microscopy in Hyperpolarization Activated Cyclic Nucleotide Gated Potassium Channel 4 (HCN4)-positive atrioventricular node and His-Purkinje fibers, with relative less expression in the HCN4-positive sinoatrial node. In Langendorff studies, at constant pressure, SPM-354 restored sinus rhythm in S1P-induced complete heart block and fully reversed S1P-mediated bradycardia. S1P3 distribution and function in the mouse ventricular cardiac conduction system suggest a direct mechanism for heart block risk that should be further studied in humans. A richer understanding of receptor and ligand usage in the pacemaker cells of the cardiac system is likely to be useful in understanding ventricular conduction in health, disease, and pharmacology.

Two hypertrophic cardiomyopathy-associated cardiac troponin I (cTnI) mutations, R146G and R21C, are located in different regions of cTnI, the inhibitory peptide and the cardiac-specific N terminus. We recently reported that these regions may interact when Ser-23/Ser-24 are phosphorylated, weakening the interaction of cTnI with cardiac TnC. Little is known about how these mutations influence the affinity of cardiac TnC for cTnI (KC-I) or contractile kinetics during β-adrenergic stimulation. Here, we tested how cTnI(R146G) or cTnI(R21C) influences contractile activation and relaxation and their response to protein kinase A (PKA). Both mutations significantly increased Ca(2+) binding affinity to cTn (KCa) and KC-I. PKA phosphorylation resulted in a similar reduction of KCa for all complexes, but KC-I was reduced only with cTnI(WT). cTnI(WT), cTnI(R146G), and cTnI(R21C) were complexed into cardiac troponin and exchanged into rat ventricular myofibrils, and contraction/relaxation kinetics were measured ± PKA phosphorylation. Maximal tension (Tmax) was maintained for cTnI(R146G)- and cTnI(R21C)-exchanged myofibrils, and Ca(2+) sensitivity of tension (pCa50) was increased. PKA phosphorylation decreased pCa50 for cTnI(WT)-exchanged myofibrils but not for either mutation. PKA phosphorylation accelerated the early slow phase relaxation for cTnI(WT) myofibrils, especially at Ca(2+) levels that the heart operates in vivo. Importantly, this effect was blunted for cTnI(R146G)- and cTnI(R21C)-exchanged myofibrils. Molecular dynamics simulations suggest both mutations inhibit formation of intra-subunit contacts between the N terminus and the inhibitory peptide of cTnI that is normally seen with WT-cTn upon PKA phosphorylation. Together, our results suggest that cTnI(R146G) and cTnI(R21C) blunt PKA modulation of activation and relaxation kinetics by prohibiting cardiac-specific N-terminal interaction with the cTnI inhibitory peptide.