Joyce Chuang

Patient-specific models of cardiac biomechanics

Patient-specific models of cardiac function have the potential to improve diagnosis and management of heart disease by integrating medical images with heterogeneous clinical measurements subject to constraints imposed by physical first principles and prior experimental knowledge. We describe new methods for creating three-dimensional patient-specific models of ventricular biomechanics in the failing heart. Three-dimensional bi-ventricular geometry is segmented from cardiac CT images at end-diastole from patients with heart failure. Human myofiber and sheet architecture is modeled using eigenvectors computed from diffusion tensor MR images from an isolated, fixed human organ-donor heart and transformed to the patient-specific geometric model using large deformation diffeomorphic mapping. Semi-automated methods were developed for optimizing the passive material properties while simultaneously computing the unloaded reference geometry of the ventricles for stress analysis. Material properties of active cardiac muscle contraction were optimized to match ventricular pressures measured by cardiac catheterization, and parameters of a lumped-parameter closed-loop model of the circulation were estimated with a circulatory adaptation algorithm making use of information derived from echocardiography. These components were then integrated to create a multi-scale model of the patient-specific heart. These methods were tested in five heart failure patients from the San Diego Veterans Affairs Medical Center who gave informed consent. The simulation results showed good agreement with measured echocardiographic and global functional parameters such as ejection fraction and peak cavity pressures.

Mouse and computational models link Mlc2v dephosphorylation to altered myosin kinetics in early cardiac disease.

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.

Patient-Specific Modeling of Dyssynchronous Heart Failure: A Case Study

The development and clinical use of patient-specific models of the heart is now a feasible goal. Models have the potential to aid in diagnosis and support decision-making in clinical cardiology. Several groups are now working on developing multi-scale models of the heart for understanding therapeutic mechanisms and better predicting clinical outcomes of interventions such as cardiac resynchronization therapy. Here we describe the methodology for generating a patient-specific model of the failing heart with a myocardial infarct and left ventricular bundle branch block. We discuss some of the remaining challenges in developing reliable patient-specific models of cardiac electromechanical activity, and identify some of the main areas for focusing future research efforts. Key challenges include: efficiently generating accurate patient-specific geometric meshes and mapping regional myofiber architecture to them; modeling electrical activation patterns based on cellular alterations in human heart failure, and estimating regional tissue conductivities based on clinically available electrocardiographic recordings; estimating unloaded ventricular reference geometry and material properties for biomechanical simulations; and parameterizing systemic models of circulatory dynamics from available hemodynamic measurements.

Increased infarct wall thickness by a bio-inert material is insufficient to prevent negative left ventricular remodeling after myocardial infarction.

Background Several injectable materials have been shown to preserve or improve cardiac function as well as prevent or slow left ventricular (LV) remodeling post-myocardial infarction (MI). However, it is unclear as to whether it is the structural support or the bioactivity of these polymers that lead to beneficial effects. Herein, we examine how passive structural enhancement of the LV wall by an increase in wall thickness affects cardiac function post-MI using a bio-inert, non-degradable synthetic polymer in an effort to better understand the mechanisms by which injectable materials affect LV remodeling. Methods and Results Poly(ethylene glycol) (PEG) gels of storage modulus G′ = 0.5±0.1 kPa were injected and polymerized in situ one week after total occlusion of the left coronary artery in female Sprague Dawley rats. The animals were imaged using magnetic resonance imaging (MRI) at 7±1 day(s) post-MI as a baseline and again post-injection 49±4 days after MI. Infarct wall thickness was statistically increased in PEG gel injected vs. control animals (p<0.01). However, animals in the polymer and control groups showed decreases in cardiac function in terms of end diastolic volume, end systolic volume and ejection fraction compared to baseline (p<0.01). The cellular response to injection was also similar in both groups. Conclusion The results of this study demonstrate that passive structural reinforcement alone was insufficient to prevent post-MI remodeling, suggesting that bioactivity and/or cell infiltration due to degradation of injectable materials are likely playing a key role in the preservation of cardiac function, thus providing a deeper understanding of the influencing properties of biomaterials necessary to prevent post-MI negative remodeling.

Myosin Light Chain Phosphorylation is Critical for Adaptation to Cardiac Stress

Background— Cardiac hypertrophy is a common response to circulatory or neurohumoral stressors as a mechanism to augment contractility. When the heart is under sustained stress, the hypertrophic response can evolve into decompensated heart failure, although the mechanism(s) underlying this transition remain largely unknown. Because phosphorylation of cardiac myosin light chain 2 (MLC2v), bound to myosin at the head-rod junction, facilitates actin-myosin interactions and enhances contractility, we hypothesized that phosphorylation of MLC2v plays a role in the adaptation of the heart to stress. We previously identified an enzyme that predominantly phosphorylates MLC2v in cardiomyocytes, cardiac myosin light-chain kinase (cMLCK), yet the role(s) played by cMLCK in regulating cardiac function in health and disease remain to be determined. Methods and Results— We found that pressure overload induced by transaortic constriction in wild-type mice reduced phosphorylated MLC2v levels by ≈40% and cMLCK levels by ≈85%. To examine how a reduction in cMLCK and the corresponding reduction in phosphorylated MLC2v affect function, we generated Mylk3 gene-targeted mice and transgenic mice overexpressing cMLCK specifically in cardiomyocytes. Pressure overload led to severe heart failure in cMLCK knockout mice but not in mice with cMLCK overexpression in which cMLCK protein synthesis exceeded degradation. The reduction in cMLCK protein during pressure overload was attenuated by inhibition of ubiquitin-proteasome protein degradation systems. Conclusions— Our results suggest the novel idea that accelerated cMLCK protein turnover by the ubiquitin-proteasome system underlies the transition from compensated hypertrophy to decompensated heart failure as a result of reduced phosphorylation of MLC2v.

Determination of three‐dimensional ventricular strain distributions in gene‐targeted mice using tagged MRI

A model‐based method for calculating three‐dimensional (3D) cardiac wall strain distributions in the mouse has been developed and tested in a genetically engineered mouse model of dilated cardiomyopathy. Data from MR tagging and harmonic phase (HARP) tracking were used to measure material point displacements, and 3D Lagrangian strains were calculated throughout the entire left ventricle (LV) with a deformable parametric model. A mouse model where cardiomyocytes are specifically made deficient in vinculin (VclKO) were compared to wild‐type (WT) littermates. 3D strain analysis revealed differences in LV wall mechanics between WT and VclKO mice at 8 weeks of age when systolic function had just begun to decline. Most notably, end‐systolic radial strain and torsional shear were reduced in VclKO hearts which contributed to regional mechanical dysfunction. This study demonstrates the feasibility of using MRI tagging methods to detect alterations in 3D myocardial strain distributions in genetically engineered mouse models of cardiovascular disease. Magn Reson Med, 2010.

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.

Accuracy of Seattle Heart Failure Model and HeartMate II Risk Score in Non–Inotrope-Dependent Advanced Heart Failure Patients: Insights From the ROADMAP Study (Risk Assessment and Comparative Effectiveness of Left Ventricular Assist Device and Medical Management in Ambulatory Heart Failure Patients)

Background— Timing of left ventricular assist device (LVAD) implantation in advanced heart failure patients not on inotropes is unclear. Relevant prediction models exist (SHFM [Seattle Heart Failure Model] and HMRS [HeartMate II Risk Score]), but use in this group is not established. Methods and Results— ROADMAP (Risk Assessment and Comparative Effectiveness of Left Ventricular Assist Device and Medical Management in Ambulatory Heart Failure Patients) is a prospective, multicenter, nonrandomized study of 200 advanced heart failure patients not on inotropes who met indications for LVAD implantation, comparing the effectiveness of HeartMate II support versus optimal medical management. We compared SHFM-predicted versus observed survival (overall survival and LVAD-free survival) in the optimal medical management arm (n=103) and HMRS-predicted versus observed survival in all LVAD patients (n=111) using Cox modeling, receiver–operator characteristic (ROC) curves, and calibration plots. In the optimal medical management cohort, the SHFM was a significant predictor of survival (hazard ratio=2.98; P <0.001; ROC area under the curve=0.71; P <0.001) but not LVAD-free survival (hazard ratio=1.41; P =0.097; ROC area under the curve=0.56; P =0.314). SHFM showed adequate calibration for survival but overestimated LVAD-free survival. In the LVAD cohort, the HMRS had marginal discrimination at 3 (Cox P =0.23; ROC area under the curve=0.71; P =0.026) and 12 months (Cox P =0.036; ROC area under the curve=0.62; P =0.122), but calibration was poor, underestimating survival across time and risk subgroups. Conclusions— In non–inotrope-dependent advanced heart failure patients receiving optimal medical management, the SHFM was predictive of overall survival but underestimated the risk of clinical worsening and LVAD implantation. Among LVAD patients, the HMRS had marginal discrimination and underestimated survival post–LVAD implantation. Clinical Trial Registration— URL: . Unique identifier: [NCT01452802][1]. [1]: /lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT01452802&atom=%2Fcirchf%2F10%2F5%2Fe003745.atomBackground— Timing of left ventricular assist device (LVAD) implantation in advanced heart failure patients not on inotropes is unclear. Relevant prediction models exist (SHFM [Seattle Heart Failure Model] and HMRS [HeartMate II Risk Score]), but use in this group is not established. Methods and Results— ROADMAP (Risk Assessment and Comparative Effectiveness of Left Ventricular Assist Device and Medical Management in Ambulatory Heart Failure Patients) is a prospective, multicenter, nonrandomized study of 200 advanced heart failure patients not on inotropes who met indications for LVAD implantation, comparing the effectiveness of HeartMate II support versus optimal medical management. We compared SHFM-predicted versus observed survival (overall survival and LVAD-free survival) in the optimal medical management arm (n=103) and HMRS-predicted versus observed survival in all LVAD patients (n=111) using Cox modeling, receiver–operator characteristic (ROC) curves, and calibration plots. In the optimal medical management cohort, the SHFM was a significant predictor of survival (hazard ratio=2.98; P<0.001; ROC area under the curve=0.71; P<0.001) but not LVAD-free survival (hazard ratio=1.41; P=0.097; ROC area under the curve=0.56; P=0.314). SHFM showed adequate calibration for survival but overestimated LVAD-free survival. In the LVAD cohort, the HMRS had marginal discrimination at 3 (Cox P=0.23; ROC area under the curve=0.71; P=0.026) and 12 months (Cox P=0.036; ROC area under the curve=0.62; P=0.122), but calibration was poor, underestimating survival across time and risk subgroups. Conclusions— In non–inotrope-dependent advanced heart failure patients receiving optimal medical management, the SHFM was predictive of overall survival but underestimated the risk of clinical worsening and LVAD implantation. Among LVAD patients, the HMRS had marginal discrimination and underestimated survival post–LVAD implantation. Clinical Trial Registration— URL: http://www.clinicaltrials.gov. Unique identifier: NCT01452802.

Patient-Reported Health-Related Quality of Life Is a Predictor of Outcomes in Ambulatory Heart Failure Patients Treated with Left Ventricular Assist Device Compared with Medical Management: Results from the ROADMAP Study (Risk Assessment and Comparative Effectiveness of Left Ventricular Assist Device and Medical Management)

Background— The prospective observational ROADMAP study (Risk Assessment and Comparative Effectiveness of Left Ventricular Assist Device and Medical Management) demonstrated that ambulatory advanced heart failure patients selected for left ventricular assist device (LVAD) were more likely to be alive at 1 year on original therapy with ≥75-m improvement in 6-minute walk distance compared with patients assigned to optimal medical management. Whether baseline health-related quality of life (hrQoL) resulted in a heterogeneity of this treatment benefit is unknown. Methods and Results— Patient-reported hrQoL was assessed with EuroQol questionnaire and visual analogue scale (VAS). We aimed to identify predictors of event-free survival and survival with acceptable hrQoL (VAS≥60). LVAD patients had significant improvement in 3 of 5 EuroQol dimensions (P<0.05), but no significant changes were observed with optimal medical management. Among patients with baseline VAS<55, survival on original treatment was lower for optimal medical management patients compared with those assigned to LVAD (58±7% versus 82±5%; P=0.004). No such difference was seen if baseline VAS was ≥55 (70±7% versus 75±9%; P=0.79). Survival on original therapy with acceptable quality of life was also more likely with LVAD versus optimal medical management if baseline VAS was <55, whereas outcomes in patients with higher baseline VAS scores were similar regardless of treatment assignment (P=0.046 for treatment arm and baseline VAS interaction). Conclusions— LVAD therapy resulted in improvement of patient health status in heart failure patients with low self-reported hrQoL, but not in patients with acceptable quality of life at the time of LVAD implantation. Patient-reported hrQoL should be integrated into decision making concerning the use and timing of LVAD therapy in heart failure patients who are symptom limited but remain ambulatory. Clinical Trial Registration— URL: http://www.ClinicalTrials.gov. Unique identifier: NCT01452802.

Multicenter Evaluation of Octreotide as Secondary Prophylaxis in Patients With Left Ventricular Assist Devices and Gastrointestinal Bleeding

Background Gastrointestinal (GI) bleeding is one of the most common complications after continuous-flow left ventricular assist device implantation. More than one third of patients with incident bleed go on to develop recurrent GI bleeding. Octreotide, a somatostatin analog, is proposed to reduce the risk of recurrent GI bleeding in this population. Methods and Results This multicenter, retrospective analysis evaluated 51 continuous-flow left ventricular assist device patients who received secondary prophylaxis with octreotide after their index GI bleed from 2009 to 2015. All patients had a hospitalization for GI bleed and received octreotide after discharge. Patient demographics, medical and medication history, and clinical characteristics of patients who rebled after receiving octreotide were compared with non–rebleeders. These data were also compared with matched historical control patients previously enrolled in the HMII (HeartMate II) clinical trials, none of whom received octreotide, to provide a context for the bleeding rates. Twelve patients (24%) who received secondary octreotide prophylaxis developed another GI bleed, whereas 39 (76%) did not. There were similar intergroup demographics; however, significantly more bleeders had a previous GI bleeding history before left ventricular assist device placement (33% versus 5%; P=0.02) and greater frequency of angiodysplasia confirmed during endoscopy (58% versus 23%; P=0.03). Fewer patients in this study experienced a recurrent GI bleed compared with a matched historical control group that did not receive octreotide (24% versus 43%; P=0.04). Conclusions Patients with continuous-flow left ventricular assist device receiving secondary prophylaxis with octreotide had a significantly lower GI bleed recurrence compared with historical controls not treated with octreotide. Additional prospective studies are needed to confirm these data.