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'.
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.
1 Biomedical Sciences Graduate Program and Departments of 2 Pharmacology and 3 Medicine, University of California, San Diego, La Jolla, California; and 4 Department of Molecular Biology, The Scripps ...Research Institute, La Jolla, California
Submitted 5 December 2006
; accepted in final form 6 February 2007
Sphingosine 1-phosphate (S1P) is released at sites of tissue injury and effects cellular responses through activation of G protein-coupled receptors. The role of S1P in regulating cardiomyocyte survival following in vivo myocardial ischemia-reperfusion (I/R) injury was examined by using mice in which specific S1P receptor subtypes were deleted. Mice lacking either S1P 2 or S1P 3 receptors and subjected to 1-h coronary occlusion followed by 2 h of reperfusion developed infarcts equivalent to those of wild-type (WT) mice. However, in S1P 2,3 receptor double-knockout mice, infarct size following I/R was increased by >50%. I/R leads to activation of ERK, JNK, and p38 MAP kinases; however, these responses were not diminished in S1P 2,3 receptor knockout compared with WT mice. In contrast, activation of Akt in response to I/R was markedly attenuated in S1P 2,3 receptor knockout mouse hearts. Neither S1P 2 nor S1P 3 receptor deletion alone impaired I/R-induced Akt activation, which suggests redundant signaling through these receptors and is consistent with the finding that deletion of either receptor alone did not increase I/R injury. The involvement of cardiomyocytes in S1P 2 and S1P 3 receptor mediated activation of Akt was tested by using cells from WT and S1P receptor knockout hearts. Akt was activated by S1P, and this was modestly diminished in cardiomyocytes from S1P 2 or S1P 3 receptor knockout mice and completely abolished in the S1P 2,3 receptor double-knockout myocytes. Our data demonstrate that activation of S1P 2 and S1P 3 receptors plays a significant role in protecting cardiomyocytes from I/R damage in vivo and implicate the release of S1P and receptor-mediated Akt activation in this process.
cardioprotection; mitogen-activated kinase; G protein-coupled receptors; infarct
Address for reprint requests and other correspondence: J. H. Brown, Dept. of Pharmacology, Univ. of California, San Diego, 9500 Gilman Dr., La Jolla, CA, 92093-0636 (e-mail: jhbrown{at}ucsd.edu )
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.
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.
Actin-myosin interactions provide the driving force underlying each heartbeat. The current view is that actin-bound regulatory proteins play a dominant role in the activation of calcium-dependent ...cardiac muscle contraction. In contrast, the relevance and nature of regulation by myosin regulatory proteins (for example, myosin light chain-2 MLC2) in cardiac muscle remain poorly understood. By integrating gene-targeted mouse and computational models, we have identified an indispensable role for ventricular Mlc2 (Mlc2v) phosphorylation in regulating cardiac muscle contraction. Cardiac myosin cycling kinetics, which directly control actin-myosin interactions, were directly affected, but surprisingly, Mlc2v phosphorylation also fed back to cooperatively influence calcium-dependent activation of the thin filament. Loss of these mechanisms produced early defects in the rate of cardiac muscle twitch relaxation and ventricular torsion. Strikingly, these defects preceded the left ventricular dysfunction of heart disease and failure in a mouse model with nonphosphorylatable Mlc2v. Thus, there is a direct and early role for Mlc2 phosphorylation in regulating actin-myosin interactions in striated muscle contraction, and dephosphorylation of Mlc2 or loss of these mechanisms can play a critical role in heart failure.
The maturation of cardiac myocytes during the immediate prenatal period coincides with changes in the mechanical properties of the extracellular matrix. We investigated the effects of extracellular ...stiffness on cardiomyocyte maturation in neonatal rat ventricular myocytes grown on collagen‐coated gels. Cells on 10‐kPa substrates developed aligned sarcomeres, while cells on stiffer substrates had unaligned sarcomeres and stress fibers. Cells generated greater mechanical force on gels with stiffness similar to that of the native myocardium than on stiffer or softer substrates. To investigate the differentiation of myocyte progenitors, we used clonal expansion of engineered human embryonic stem cells. Puromycin‐selected cardiomyocytes exhibited a gene expression profile similar to that of adult human cardiomyocytes and generated force and action potentials consistent with normal fetal cardiomyocytes. These results suggest that extracellular stiffness significantly affects maturation and differentiation of immature ventricular myocytes.
Mechanics-based cardiac growth models can now predict changes in mass, chamber size, and wall thickness in response to perturbations such as pressure overload (PO), volume overload, and myocardial ...infarction with a single set of growth parameters. As these models move toward clinical applications, many of the most interesting applications involve predictions of whether or how a patient’s heart will reverse its growth after an intervention. In the case of PO, significant regression in wall thickness is observed both experimentally and clinically following relief of overload, for example following replacement of a stenotic aortic valve. Therefore, the objective of this work was to evaluate the ability of a published cardiac growth model that captures forward growth in multiple situations to predict growth reversal following relief of PO. Using a finite element model of a beating canine heart coupled to a circuit model of the circulation, we quantitatively matched hemodynamic data from a canine study of aortic banding followed by unbanding. Surprisingly, although the growth model correctly predicted the time course of PO-induced hypertrophy, it predicted only limited growth reversal given the measured unbanding hemodynamics, contradicting experimental and clinical observations. We were able to resolve this discrepancy only by incorporating an evolving homeostatic setpoint for the governing growth equations. Furthermore, our analysis suggests that many strain- and stress-based growth laws using the traditional volumetric growth framework will have similar difficulties capturing regression following the relief of PO unless growth setpoints are allowed to evolve.
Although a truly complete understanding of whole heart activation, contraction, and deformation is well beyond our current reach, a significant amount of effort has been devoted to discovering and ...understanding the mechanisms by which myocardial structure determines cardiac function to better treat patients with cardiac disease. Several experimental studies have shown that transmural fiber strain is relatively uniform in both diastole and systole, in contrast to predictions from traditional mechanical theory. Similarly, mathematical models have largely predicted uniform fiber stress across the wall. The development of this uniform pattern of fiber stress and strain during filling and ejection is due to heterogeneous transmural distributions of several myocardial structures. This review summarizes these transmural gradients, their contributions to fiber mechanics, and the potential functional effects of their remodeling during pressure overload hypertrophy.
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.
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.
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.
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.
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.
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.
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, R
= 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.
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.