Successful cardiac resynchronization therapy reduces negative septal work in patient-specific models of dyssynchronous heart failure.

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.

Differential sensitivity to longitudinal and transverse stretch mediates transcriptional responses in mouse neonatal ventricular myocytes.

To identify how cardiomyocyte mechanosensitive signaling pathways are regulated by anisotropic stretch, micropatterned mouse neonatal cardiomyocytes were stretched primarily longitudinally or transversely to the myofiber axis. Four hours of static, longitudinal stretch induced differential expression of 557 genes, compared with 30 induced by transverse stretch, measured using RNA-seq. A logic-based ordinary differential equation model of the cardiac myocyte mechanosignaling network, extended to include the transcriptional regulation and expression of 784 genes, correctly predicted measured expression changes due to anisotropic stretch with 69% accuracy. The model also predicted published transcriptional responses to mechanical load in vitro or in vivo with 63-91% accuracy. The observed differences between transverse and longitudinal stretch responses were not explained by differential activation of specific pathways but rather by an approximately twofold greater sensitivity to longitudinal stretch than transverse stretch. In vitro experiments confirmed model predictions that stretch-induced gene expression is more sensitive to angiotensin II and endothelin-1, via RhoA and MAP kinases, than to the three membrane ion channels upstream of calcium signaling in the network. Quantitative cardiomyocyte gene expression differs substantially with the axis of maximum principal stretch relative to the myofilament axis, but this difference is due primarily to differences in stretch sensitivity rather than to selective activation of mechanosignaling pathways.NEW & NOTEWORTHY Anisotropic stretch applied to micropatterned neonatal mouse ventricular myocytes induced markedly greater acute transcriptional responses when the major axis of stretch was parallel to the myofilament axis than when it was transverse. Analysis with a novel quantitative network model of mechanoregulated cardiomyocyte gene expression suggests that this difference is explained by higher cell sensitivity to longitudinal loading than transverse loading than by the activation of differential signaling pathways.

Signaling network model of cardiomyocyte morphological changes in familial cardiomyopathy.

Familial cardiomyopathy is a precursor of heart failure and sudden cardiac death. Over the past several decades, researchers have discovered numerous gene mutations primarily in sarcomeric and cytoskeletal proteins causing two different disease phenotypes: hypertrophic (HCM) and dilated (DCM) cardiomyopathies. However, molecular mechanisms linking genotype to phenotype remain unclear. Here, we employ a systems approach by integrating experimental findings from preclinical studies (e.g., murine data) into a cohesive signaling network to scrutinize genotype to phenotype mechanisms. We developed an HCM/DCM signaling network model utilizing a logic-based differential equations approach and evaluated model performance in predicting experimental data from four contexts (HCM, DCM, pressure overload, and volume overload). The model has an overall prediction accuracy of 83.8%, with higher accuracy in the HCM context (90%) than DCM (75%). Global sensitivity analysis identifies key signaling reactions, with calcium-mediated myofilament force development and calcium-calmodulin kinase signaling ranking the highest. A structural revision analysis indicates potential missing interactions that primarily control calcium regulatory proteins, increasing model prediction accuracy. Combination pharmacotherapy analysis suggests that downregulation of signaling components such as calcium, titin and its associated proteins, growth factor receptors, ERK1/2, and PI3K-AKT could inhibit myocyte growth in HCM. In experiments with patient-specific iPSC-derived cardiomyocytes (MLP-W4R;MYH7-R723C iPSC-CMs), combined inhibition of ERK1/2 and PI3K-AKT rescued the HCM phenotype, as predicted by the model. In DCM, PI3K-AKT-NFAT downregulation combined with upregulation of Ras/ERK1/2 or titin or Gq protein could ameliorate cardiomyocyte morphology. The model results suggest that HCM mutations that increase active force through elevated calcium sensitivity could increase ERK activity and decrease eccentricity through parallel growth factors, Gq-mediated, and titin pathways. Moreover, the model simulated the influence of existing medications on cardiac growth in HCM and DCM contexts. This HCM/DCM signaling model demonstrates utility in investigating genotype to phenotype mechanisms in familial cardiomyopathy.

A deep learning approach for fully automated cardiac shape modeling in tetralogy of Fallot.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

Le Cœur en Sabot: shape associations with adverse events in repaired tetralogy of Fallot.

BACKGROUND: Maladaptive remodelling mechanisms occur in patients with repaired tetralogy of Fallot (rToF) resulting in a cycle of metabolic and structural changes. Biventricular shape analysis may indicate mechanisms associated with adverse events independent of pulmonary regurgitant volume index (PRVI). We aimed to determine novel remodelling patterns associated with adverse events in patients with rToF using shape and function analysis. METHODS: Biventricular shape and function were studied in 192 patients with rToF (median time from TOF repair to baseline evaluation 13.5 years). Linear discriminant analysis (LDA) and principal component analysis (PCA) were used to identify shape differences between patients with and without adverse events. Adverse events included death, arrhythmias, and cardiac arrest with median follow-up of 10 years. RESULTS: LDA and PCA showed that shape characteristics pertaining to adverse events included a more circular left ventricle (LV) (decreased eccentricity), dilated (increased sphericity) LV base, increased right ventricular (RV) apical sphericity, and decreased RV basal sphericity. Multivariate LDA showed that the optimal discriminative model included only RV apical ejection fraction and one PCA mode associated with a more circular and dilated LV base (AUC = 0.77). PRVI did not add value, and shape changes associated with increased PRVI were not predictive of adverse outcomes. CONCLUSION: Pathological remodelling patterns in patients with rToF are significantly associated with adverse events, independent of PRVI. Mechanisms related to incident events include LV basal dilation with a reduced RV apical ejection fraction.

Atlas-based methods for efficient characterization of patient-specific ventricular activation patterns.

AIMS: Ventricular 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 RESULTS: Atlases 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. CONCLUSION: Atlases 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.

Regional variations in ex-vivo diffusion tensor anisotropy are associated with cardiomyocyte remodeling in rats after left ventricular pressure overload.

BACKGROUND: Pressure 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. METHODS: We 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. RESULTS: We 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 < 0.001, R2 = 0.71). There was a transmural gradient in myocyte cross-sectional area in SHAM hearts that was diminished in the TAC group. We also found several regions of significantly altered DTI parameters in TAC LV compared to SHAM, especially in myocyte sheet angle dispersion and planar anisotropy. Among others, these parameters correlated significantly with directly measured myocyte aspect ratios. CONCLUSIONS: These results show that structural remodeling in pressure overload LV hypertrophy is regionally heterogeneous, especially transmurally, with a greater degree of remodeling in the sub-endocardium compared to the sub-epicardium. Additionally, several parameters derived from DTI correlated significantly with measurements of myocyte geometry from direct measurement in histological sections. We suggest that DTI may provide a non-invasive, comprehensive method to detect regional structural myocyte LV remodeling during disease.

Quantification of model and data uncertainty in a network analysis of cardiac myocyte mechanosignalling.

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'.

Predictive model identifies key network regulators of cardiomyocyte mechano-signaling.

Mechanical strain is a potent stimulus for growth and remodeling in cells. Although many pathways have been implicated in stretch-induced remodeling, the control structures by which signals from distinct mechano-sensors are integrated to modulate hypertrophy and gene expression in cardiomyocytes remain unclear. Here, we constructed and validated a predictive computational model of the cardiac mechano-signaling network in order to elucidate the mechanisms underlying signal integration. The model identifies calcium, actin, Ras, Raf1, PI3K, and JAK as key regulators of cardiac mechano-signaling and characterizes crosstalk logic imparting differential control of transcription by AT1R, integrins, and calcium channels. We find that while these regulators maintain mostly independent control over distinct groups of transcription factors, synergy between multiple pathways is necessary to activate all the transcription factors necessary for gene transcription and hypertrophy. We also identify a PKG-dependent mechanism by which valsartan/sacubitril, a combination drug recently approved for treating heart failure, inhibits stretch-induced hypertrophy, and predict further efficacious pairs of drug targets in the network through a network-wide combinatorial search.

Mechanotransduction in cardiac hypertrophy and failure.

Cardiac muscle cells have an intrinsic ability to sense and respond to mechanical load through a process known as mechanotransduction. In the heart, this process involves the conversion of mechanical stimuli into biochemical events that induce changes in myocardial structure and function. Mechanotransduction and its downstream effects function initially as adaptive responses that serve as compensatory mechanisms during adaptation to the initial load. However, under prolonged and abnormal loading conditions, the remodeling processes can become maladaptive, leading to altered physiological function and the development of pathological cardiac hypertrophy and heart failure. Although the mechanisms underlying mechanotransduction are far from being fully elucidated, human and mouse genetic studies have highlighted various cytoskeletal and sarcolemmal structures in cardiac myocytes as the likely candidates for load transducers, based on their link to signaling molecules and architectural components important in disease pathogenesis. In this review, we summarize recent developments that have uncovered specific protein complexes linked to mechanotransduction and mechanotransmission within the sarcomere, the intercalated disc, and at the sarcolemma. The protein structures acting as mechanotransducers are the first step in the process that drives physiological and pathological cardiac hypertrophy and remodeling, as well as the transition to heart failure, and may provide better insights into mechanisms driving mechanotransduction-based diseases.

Connexin defects underlie arrhythmogenic right ventricular cardiomyopathy in a novel mouse model.

Arrhythmogenic right ventricular cardiomyopathy (ARVC) termed a 'disease of the desmosome' is an inherited cardiomyopathy that recently underwent reclassification owing to the identification of left-dominant and biventricular disease forms. Homozygous loss-of-function mutations in the desmosomal component, desmoplakin, are found in patients exhibiting a biventricular form of ARVC; however, no models recapitulate the postnatal hallmarks of the disease as seen in these patients. To gain insights into the homozygous loss-of-function effects of desmoplakin in the heart, we generated cardiomyocyte-specific desmoplakin-deficient mice (DSP-cKO) using ventricular myosin light chain-2-Cre mice. Homozygous DSP-cKO mice are viable but display early ultrastructural defects in desmosomal integrity leading to a cardiomyopathy reminiscent of a biventricular form of ARVC, which includes cell death and fibro-fatty replacement within the ventricle leading to biventricular dysfunction, failure and premature death. DSP-cKO mice also exhibited ventricular arrhythmias that are exacerbated with exercise and catecholamine stimulation. Furthermore, DSP-cKO hearts exhibited right ventricular conduction defects associated with loss of connexin 40 expression and electrical wavefront propagation defects associated with loss of connexin 43 expression. Dose-dependent assessment of the effects of loss of desmoplakin in neonatal ventricular cardiomyocytes revealed primary loss of connexin 43 levels, phosphorylation and function independent of the molecular dissociation of the mechanical junction complex and fibro-fatty manifestation associated with ARVC, suggesting a role for desmoplakin as a primary stabilizer of connexin integrity. In summary, we provide evidence for a novel mouse model, which is reminiscent of the postnatal onset of ARVC while highlighting mechanisms underlying a biventricular form of human ARVC.

Non-invasive, model-based measures of ventricular electrical dyssynchrony for predicting CRT outcomes.

AIMS: Left ventricular activation delay due to left bundle branch block (LBBB) is an important determinant of the severity of dyssynchronous heart failure (DHF). We investigated whether patient-specific computational models constructed from non-invasive measurements can provide measures of baseline dyssynchrony and its reduction after CRT that may explain the degree of long-term reverse ventricular remodelling. METHODS AND RESULTS: LV end-systolic volume reduction (ΔESVLV) measured by 2D trans-thoracic echocardiography in eight patients following 6 months of CRT was significantly (P < 0.05) greater in responders (26 ± 20%, n = 4) than non-responders (11 ± 16%, n = 4). LV reverse remodelling did not correlate with baseline QRS duration or its change after biventricular pacing, but did correlate with baseline LV endocardial activation measured by electroanatomic mapping (R2 = 0.71, P < 0.01). Patient-specific models of LBBB ventricular activation with parameters obtained by matching model-computed vectorcardiograms (VCG) to those derived from standard patient ECGs yielded LV endocardial activation times that correlated well with those measured from endocardial maps (R2 = 0.90). Model-computed 3D LV activation times correlated strongly with the reduction in LVESV (R2 = 0.93, P < 0.001). Computed decreases due to simulated CRT in the time delay between LV septal and lateral activation correlated strongly with ΔESVLV (R2 = 0.92, P < 0.001). Models also suggested that optimizing VV delays may improve resynchronization by this measure of activation delay. CONCLUSIONS: Patient-specific computational models constructed from non-invasive measurements can compute estimates of LV dyssynchrony and their changes after CRT that may be as good as or better than electroanatomic mapping for predicting long-term reverse remodelling.

Atlas-Based Ventricular Shape Analysis for Understanding Congenital Heart Disease.

Congenital heart disease is associated with abnormal ventricular shape that can affect wall mechanics and may be predictive of long-term adverse outcomes. Atlas-based parametric shape analysis was used to analyze ventricular geometries of eight adolescent or adult single-ventricle CHD patients with tricuspid atresia and Fontans. These patients were compared with an "atlas" of non-congenital asymptomatic volunteers, resulting in a set of z-scores which quantify deviations from the control population distribution on a patient-by-patient basis. We examined the potential of these scores to: (1) quantify abnormalities of ventricular geometry in single ventricle physiologies relative to the normal population; (2) comprehensively quantify wall motion in CHD patients; and (3) identify possible relationships between ventricular shape and wall motion that may reflect underlying functional defects or remodeling in CHD patients. CHD ventricular geometries at end-diastole and end-systole were individually compared with statistical shape properties of an asymptomatic population from the Cardiac Atlas Project. Shape analysis-derived model properties, and myocardial wall motions between end-diastole and end-systole, were compared with physician observations of clinical functional parameters. Relationships between altered shape and altered function were evaluated via correlations between atlas-based shape and wall motion scores. Atlas-based shape analysis identified a diverse set of specific quantifiable abnormalities in ventricular geometry or myocardial wall motion in all subjects. Moreover, this initial cohort displayed significant relationships between specific shape abnormalities such as increased ventricular sphericity and functional defects in myocardial deformation, such as decreased long-axis wall motion. These findings suggest that atlas-based ventricular shape analysis may be a useful new tool in the management of patients with CHD who are at risk of impaired ventricular wall mechanics and chamber remodeling.

Modifying Ventricular Fibrillation by Targeted Rotor Substrate Ablation: Proof-of-Concept from Experimental Studies to Clinical VF.

INTRODUCTION: Recent work has suggested a role for organized sources in sustaining ventricular fibrillation (VF). We assessed whether ablation of rotor substrate could modulate VF inducibility in canines, and used this proof-of-concept as a foundation to suppress antiarrhythmic drug-refractory clinical VF in a patient with structural heart disease. METHODS AND RESULTS: In 9 dogs, we introduced 64-electrode basket catheters into one or both ventricles, used rapid pacing at a recorded induction threshold to initiate VF, and then defibrillated after 18±8 seconds. Endocardial rotor sites were identified from basket recordings using phase mapping, and ablation was performed at nonrotor (sham) locations (7 ± 2 minutes) and then at rotor sites (8 ± 2 minutes, P = 0.10 vs. sham); the induction threshold was remeasured after each. Sham ablation did not alter canine VF induction threshold (preablation 150 ± 16 milliseconds, postablation 144 ± 16 milliseconds, P = 0.54). However, rotor site ablation rendered VF noninducible in 6/9 animals (P = 0.041), and increased VF induction threshold in the remaining 3. Clinical proof-of-concept was performed in a patient with repetitive ICD shocks due to VF refractory to antiarrhythmic drugs. Following biventricular basket insertion, VF was induced and then defibrillated. Mapping identified 4 rotors localized at borderzone tissue, and rotor site ablation (6.3 ± 1.5 minutes/site) rendered VF noninducible. The VF burden fell from 7 ICD shocks in 8 months preablation to zero ICD therapies at 1 year, without antiarrhythmic medications. CONCLUSIONS: Targeted rotor substrate ablation suppressed VF in an experimental model and a patient with refractory VF. Further studies are warranted on the efficacy of VF source modulation.

Timing and magnitude of systolic stretch affect myofilament activation and mechanical work.

Dyssynchronous activation of the heart leads to abnormal regional systolic stretch. In vivo studies have suggested that the timing of systolic stretch can affect regional tension and external work development. In the present study, we measured the direct effects of systolic stretch timing on the magnitude of tension and external work development in isolated murine right ventricular papillary muscles. A servomotor was used to impose precisely timed stretches relative to electrical activation while a force transducer measured force output and strain was monitored using a charge-couple device camera and topical markers. Stretches taking place during peak intracellular Ca(2+) statistically increased peak tension up to 270%, whereas external work due to stretches in this interval reached values of 500 J/m. An experimental analysis showed that time-varying elastance overestimated peak tension by 100% for stretches occurring after peak isometric tension. The addition of the force-velocity relation explained some effects of stretches occurring before the peak of the Ca(2+) transient but had no effect in later stretches. An estimate of transient deactivation was measured by performing quick stretches to dissociate cross-bridges. The timing of transient deactivation explained the remaining differences between the model and experiment. These results suggest that stretch near the start of cardiac tension development substantially increases twitch tension and mechanical work production, whereas late stretches decrease external work. While the increased work can mostly be explained by the time-varying elastance of cardiac muscle, the decreased work in muscles stretched after the peak of the Ca(2+) transient is largely due to myofilament deactivation.

Biomechanics of cardiac electromechanical coupling and mechanoelectric feedback.

Cardiac mechanical contraction is triggered by electrical activation via an intracellular calcium-dependent process known as excitation-contraction coupling. Dysregulation of cardiac myocyte intracellular calcium handling is a common feature of heart failure. At the organ scale, electrical dyssynchrony leads to mechanical alterations and exacerbates pump dysfunction in heart failure. A reverse coupling between cardiac mechanics and electrophysiology is also well established. It is commonly referred as cardiac mechanoelectric feedback and thought to be an important contributor to the increased risk of arrhythmia during pathological conditions that alter regional cardiac wall mechanics, including heart failure. At the cellular scale, most investigations of myocyte mechanoelectric feedback have focused on the roles of stretch-activated ion channels, though mechanisms that are independent of ionic currents have also been described. Here we review excitation-contraction coupling and mechanoelectric feedback at the cellular and organ scales, and we identify the need for new multicellular tissue-scale model systems and experiments that can help us to obtain a better understanding of how interactions between electrophysiological and mechanical processes at the cell scale affect ventricular electromechanical interactions at the organ scale in the normal and diseased heart.

Successful Cardiac Resynchronization Therapy Reduces Negative Septal Work in Patient-Specific Models of Dyssynchronous Heart Failure.

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.

A universal biventricular coordinate system incorporating valve annuli: Validation in congenital heart disease.

Universal coordinate systems have been proposed to facilitate anatomic registration between three-dimensional images, data and models of the ventricles of the heart. However, current universal ventricular coordinate systems do not account for the outflow tracts and valve annuli where the anatomy is complex. Here we propose an extension to the 'Cobiveco' biventricular coordinate system that also accounts for the intervalvular bridges of the base and provides a tool for anatomically consistent registration between widely varying biventricular shapes. CobivecoX uses a novel algorithm to separate intervalvular bridges and assign new coordinates, including an inflow-outflow coordinate, to describe local positions in these regions uniquely and consistently. Anatomic consistency of registration was validated using curated three-dimensional biventricular shape models derived from cardiac MRI measurements in normal hearts and hearts from patients with congenital heart diseases. This new method allows the advantages of universal cardiac coordinates to be used for three-dimensional ventricular imaging data and models that include the left and right ventricular outflow tracts and valve annuli.

Novel role for vinculin in ventricular myocyte mechanics and dysfunction.

Vinculin (Vcl) plays a key structural role in ventricular myocytes that, when disrupted, can lead to contractile dysfunction and dilated cardiomyopathy. To investigate the role of Vcl in myocyte and myocardial function, cardiomyocyte-specific Vcl knockout mice (cVclKO) and littermate control wild-type mice were studied with transmission electron microscopy (TEM) and in vivo magnetic resonance imaging (MRI) tagging before the onset of global ventricular dysfunction. MRI revealed significantly decreased systolic strains transverse to the myofiber axis in vivo, but no changes along the muscle fibers or in fiber tension in papillary muscles from heterozygous global Vcl null mice. Myofilament lattice spacing from TEM was significantly greater in cVclKO versus wild-type hearts fixed in the unloaded state. AFM in Vcl heterozygous null mouse myocytes showed a significant decrease in membrane cortical stiffness. A multiscale computational model of ventricular mechanics incorporating cross-bridge geometry and lattice mechanics showed that increased transverse systolic stiffness due to increased lattice spacing may explain the systolic wall strains associated with Vcl deficiency, before the onset of ventricular dysfunction. Loss of cardiac myocyte Vcl may decrease systolic transverse strains in vivo by decreasing membrane cortical tension, which decreases transverse compression of the lattice thereby increasing interfilament spacing and stress transverse to the myofibers.

A novel mechanism involving four-and-a-half LIM domain protein-1 and extracellular signal-regulated kinase-2 regulates titin phosphorylation and mechanics.

Understanding mechanisms underlying titin regulation in cardiac muscle function is of critical importance given recent compelling evidence that highlight titin mutations as major determinants of human cardiomyopathy. We previously identified a cardiac biomechanical stress-regulated complex at the cardiac-specific N2B region of titin that includes four-and-a-half LIM domain protein-1 (Fhl1) and components of the mitogen-activated protein signaling cascade, which impacted muscle compliance in Fhl1 knock-out cardiac muscle. However, direct regulation of these molecular components in mediating titin N2B function remained unresolved. Here we identify Fhl1 as a novel negative regulator of titin N2B levels and phosphorylation-mediated mechanics. We specifically identify titin N2B as a novel substrate of extracellular signal regulated-kinase-2 (Erk2) and demonstrate that Fhl1 directly interferes with Erk2-mediated titin-N2B phosphorylation. We highlight the critical region in titin-N2B that interacts with Fhl1 and residues that are dependent on Erk2-mediated phosphorylation in situ. We also propose a potential mechanism for a known titin-N2B cardiomyopathy-causing mutation that involves this regulatory complex. These studies shed light on a novel mechanism regulating titin-N2B mechano-signaling as well as suggest that dysfunction of these pathways could be important in cardiac disease states affecting muscle compliance.