Lik Chuan Lee

Distribution of normal human left ventricular myofiber stress at end diastole and end systole: a target for in silico design of heart failure treatments

Ventricular wall stress is believed to be responsible for many physical mechanisms taking place in the human heart, including ventricular remodeling, which is frequently associated with heart failure. Therefore, normalization of ventricular wall stress is the cornerstone of many existing and new treatments for heart failure. In this paper, we sought to construct reference maps of normal ventricular wall stress in humans that could be used as a target for in silico optimization studies of existing and potential new treatments for heart failure. To do so, we constructed personalized computational models of the left ventricles of five normal human subjects using magnetic resonance images and the finite-element method. These models were calibrated using left ventricular volume data extracted from magnetic resonance imaging (MRI) and validated through comparison with strain measurements from tagged MRI (950 ± 170 strain comparisons/subject). The calibrated passive material parameter values were C0 = 0.115 ± 0.008 kPa and B0 = 14.4 ± 3.18; the active material parameter value was Tmax = 143 ± 11.1 kPa. These values could serve as a reference for future construction of normal human left ventricular computational models. The differences between the predicted and the measured circumferential and longitudinal strains in each subject were 3.4 ± 6.3 and 0.5 ± 5.9%, respectively. The predicted end-diastolic and end-systolic myofiber stress fields for the five subjects were 2.21 ± 0.58 and 16.54 ± 4.73 kPa, respectively. Thus these stresses could serve as targets for in silico design of heart failure treatments.

Algisyl-LVR™ with coronary artery bypass grafting reduces left ventricular wall stress and improves function in the failing human heart.

BACKGROUND Left ventricular (LV) wall stress reduction is a cornerstone in treating heart failure. Large animal models and computer simulations indicate that adding non-contractile material to the damaged LV wall can potentially reduce myofiber stress. We sought to quantify the effects of a novel implantable hydrogel (Algisyl-LVR™) treatment in combination with coronary artery bypass grafting (i.e. Algisyl-LVR™+CABG) on both LV function and wall stress in heart failure patients. METHODS AND RESULTS Magnetic resonance images obtained before treatment (n=3), and at 3 months (n=3) and 6 months (n=2) afterwards were used to reconstruct the LV geometry. Cardiac function was quantified using end-diastolic volume (EDV), end-systolic volume (ESV), regional wall thickness, sphericity index and regional myofiber stress computed using validated mathematical modeling. The LV became more ellipsoidal after treatment, and both EDV and ESV decreased substantially 3 months after treatment in all patients; EDV decreased from 264 ± 91 ml to 146 ± 86 ml and ESV decreased from 184 ± 85 ml to 86 ± 76 ml. Ejection fraction increased from 32 ± 8% to 47 ± 18% during that period. Volumetric-averaged wall thickness increased in all patients, from 1.06 ± 0.21 cm (baseline) to 1.3 ± 0.26 cm (3 months). These changes were accompanied by about a 35% decrease in myofiber stress at end-of-diastole and at end-of-systole. Post-treatment myofiber stress became more uniform in the LV. CONCLUSIONS These results support the novel concept that Algisyl-LVR™+CABG treatment leads to decreased myofiber stress, restored LV geometry and improved function.

First Evidence of Depressed Contractility in the Border Zone of a Human Myocardial Infarction

BACKGROUND The temporal progression in extent and severity of regional myofiber contractile dysfunction in normally perfused border zone (BZ) myocardium adjacent to a myocardial infarction (MI) has been shown to be an important pathophysiologic feature of the adverse remodeling process in large animal models. We sought, for the first time, to document the presence of impaired contractility of the myofibers in the human BZ myocardium. METHODS A 62-year-old man who experienced an MI in 1985 and had recently had complete revascularization was studied. Myofiber systolic contractile stress developed in the normally perfused BZ adjacent to the MI (T(max_B)) and that developed in regions remote from the MI (T(max_R)) were quantified using cardiac catheterization, magnetic resonance imaging, and mathematical modeling. RESULTS The resulting finite element model of the patients beating left ventricle was able to simulate the reduced systolic strains measured using magnetic resonance imaging at matching left ventricular pressures and volumes. The T(max_B) (73.1 kPa) was found to be greatly reduced relative to T(max_R) (109.5 kPa). These results were found to be independent of assumptions relating to BZ myofiber orientation. CONCLUSIONS The results of this study document the presence of impaired contractility of the myofibers in the BZ myocardium and support its role in the post-MI remodeling process in patients. To fully establish this important conclusion serial evaluations beginning at the time of the index MI will need to be performed in a cohort of patients. The current study supports the importance and demonstrates the feasibility of larger and longer-term studies.

Heterogeneous growth-induced prestrain in the heart

Even when entirely unloaded, biological structures are not stress-free, as shown by Y.C. Fung׳s seminal opening angle experiment on arteries and the left ventricle. As a result of this prestrain, subject-specific geometries extracted from medical imaging do not represent an unloaded reference configuration necessary for mechanical analysis, even if the structure is externally unloaded. Here we propose a new computational method to create physiological residual stress fields in subject-specific left ventricular geometries using the continuum theory of fictitious configurations combined with a fixed-point iteration. We also reproduced the opening angle experiment on four swine models, to characterize the range of normal opening angle values. The proposed method generates residual stress fields which can reliably reproduce the range of opening angles between 8.7±1.8 and 16.6±13.7 as measured experimentally. We demonstrate that including the effects of prestrain reduces the left ventricular stiffness by up to 40%, thus facilitating the ventricular filling, which has a significant impact on cardiac function. This method can improve the fidelity of subject-specific models to improve our understanding of cardiac diseases and to optimize treatment options.

Patient-specific finite element based analysis of ventricular myofiber stress after Coapsys: Importance of residual stress.

BACKGROUND We sought to determine regional myofiber stress after Coapsys device (Myocor, Inc, Maple Grove, MN) implantation using a finite element model of the left ventricle (LV). Chronic ischemic mitral regurgitation is caused by LV remodeling after posterolateral myocardial infarction. The Coapsys device consists of a single trans-LV chord placed below the mitral valve such that when tensioned it alters LV shape and decreases chronic ischemic mitral regurgitation. METHODS Finite element models of the LV were based on magnetic resonance images obtained before (preoperatively) and after (postoperatively) coronary artery bypass grafting with Coapsys implantation in a single patient. To determine the effect of Coapsys and LV before stress, virtual Coapsys was performed on the preoperative model. Diastolic and systolic material variables in the preoperative, postoperative, and virtual Coapsys models were adjusted so that model LV volume agreed with magnetic resonance imaging data. Chronic ischemic mitral regurgitation was abolished in the postoperative models. In each case, myofiber stress and pump function were calculated. RESULTS Both postoperative and virtual Coapsys models shifted end-systolic and end-diastolic pressure-volume relationships to the left. As a consequence and because chronic ischemic mitral regurgitation was reduced after Coapsys, pump function was unchanged. Coapsys decreased myofiber stress at end-diastole and end-systole in both the remote and infarct regions of the myocardium. However, knowledge of Coapsys and LV prestress was necessary for accurate calculation of LV myofiber stress, especially in the remote zone. CONCLUSIONS Coapsys decreases myofiber stress at end-diastole and end-systole. The improvement in myofiber stress may contribute to the long-term effect of Coapsys on LV remodeling.

Modeling Pathologies of Diastolic and Systolic Heart Failure

AbstractChronic heart failure is a medical condition that involves structural and functional changes of the heart and a progressive reduction in cardiac output. Heart failure is classified into two categories: diastolic heart failure, a thickening of the ventricular wall associated with impaired filling; and systolic heart failure, a dilation of the ventricles associated with reduced pump function. In theory, the pathophysiology of heart failure is well understood. In practice, however, heart failure is highly sensitive to cardiac microstructure, geometry, and loading. This makes it virtually impossible to predict the time line of heart failure for a diseased individual. Here we show that computational modeling allows us to integrate knowledge from different scales to create an individualized model for cardiac growth and remodeling during chronic heart failure. Our model naturally connects molecular events of parallel and serial sarcomere deposition with cellular phenomena of myofibrillogenesis and sarcomerogenesis to whole organ function. Our simulations predict chronic alterations in wall thickness, chamber size, and cardiac geometry, which agree favorably with the clinical observations in patients with diastolic and systolic heart failure. In contrast to existing single- or bi-ventricular models, our new four-chamber model can also predict characteristic secondary effects including papillary muscle dislocation, annular dilation, regurgitant flow, and outflow obstruction. Our prototype study suggests that computational modeling provides a patient-specific window into the progression of heart failure with a view towards personalized treatment planning.

Human Cardiac Function Simulator for the Optimal Design of a Novel Annuloplasty Ring with a Sub-valvular Element for Correction of Ischemic Mitral Regurgitation

AbstractIschemic mitral regurgitation is associated with substantial risk of death. We sought to: (1) detail significant recent improvements to the Dassault Systèmes human cardiac function simulator (HCFS); (2) use the HCFS to simulate normal cardiac function as well as pathologic function in the setting of posterior left ventricular (LV) papillary muscle infarction; and (3) debut our novel device for correction of ischemic mitral regurgitation. We synthesized two recent studies of human myocardial mechanics. The first study presented the robust and integrative finite element HCFS. Its primary limitation was its poor diastolic performance with an LV ejection fraction below 20% caused by overly stiff ex vivo porcine tissue parameters. The second study derived improved diastolic myocardial material parameters using in vivo MRI data from five normal human subjects. We combined these models to simulate ischemic mitral regurgitation by computationally infarcting an LV region including the posterior papillary muscle. Contact between our novel device and the mitral valve apparatus was simulated using Dassault Systèmes SIMULIA software. Incorporating improved cardiac geometry and diastolic myocardial material properties in the HCFS resulted in a realistic LV ejection fraction of 55%. Simulating infarction of posterior papillary muscle caused regurgitant mitral valve mechanics. Implementation of our novel device corrected valve dysfunction. Improvements in the current study to the HCFS permit increasingly accurate study of myocardial mechanics. The first application of this simulator to abnormal human cardiac function suggests that our novel annuloplasty ring with a sub-valvular element will correct ischemic mitral regurgitation.

Utility of high-resolution electroanatomic mapping of the left ventricle using a multispline basket catheter in a swine model of chronic myocardial infarction.

BACKGROUND Standard electroanatomic mapping systems use a single catheter to perform left ventricular substrate mapping. A new mapping system uses a 64-electrode mini-basket catheter to perform rapid automated acquisition of chamber geometry, voltage, and activation. OBJECTIVE The aim of this study was to compare the accuracy of electroanatomic mapping using the basket catheter with that of mapping using a standard linear catheter in a swine model of chronic myocardial infarction. METHODS Ten swine underwent left anterior descending coronary artery occlusion to create an anteroseptal myocardial infarction. Animals underwent delayed-enhancement magnetic resonance imaging (MRI) and then detailed left ventricular voltage mapping with both the basket and the linear catheter. Map characteristics and scar area were compared between the basket catheter, linear catheter, and MRI. Induced ventricular tachycardia (VT) was mapped with the basket catheter. RESULTS More points were acquired with the basket catheter than with the standard catheter (8762 ± 3164 vs 1712 ± 551; P < .001). The fifth percentile of normal bipolar voltage distribution with the basket catheter was 1.54 mV. Using a bipolar voltage cutoff of 1.5 mV, the total infarct areas measured using the basket catheter, linear catheter, and MRI were similar (17.8 cm(2) vs 20.9 cm(2) vs 17.5 cm(2); P = .69); however, the correlation between MRI and catheter scar area measurement was best for the basket catheter (basket vs linear: r = .76 vs r = .71). In 3 animals, sustained poorly tolerated VT was initiated and the circuit mapped successfully with the basket catheter in <5 minutes. CONCLUSION Rapid substrate and activation mapping using a 64-electrode mini-basket catheter allows detailed voltage and activation mapping in postinfarction cardiomyopathy. This system may be useful for substrate and VT mapping in human postinfarction cardiomyopathy.

Applications of Computational Modeling in Cardiac Surgery

Although computational modeling is common in many areas of science and engineering, only recently have advances in experimental techniques and medical imaging allowed this tool to be applied in cardiac surgery. Despite its infancy in cardiac surgery, computational modeling has been useful in calculating the effects of clinical devices and surgical procedures. In this review, we present several examples that demonstrate the capabilities of computational cardiac modeling in cardiac surgery. Specifically, we demonstrate its ability to simulate surgery, predict myofiber stress and pump function, and quantify changes to regional myocardial material properties. In addition, issues that would need to be resolved in order for computational modeling to play a greater role in cardiac surgery are discussed. doi: 10.1111/jocs.12332 (J Card Surg 2014;29:293–302)

A computational model that predicts reverse growth in response to mechanical unloading

Ventricular growth is widely considered to be an important feature in the adverse progression of heart diseases, whereas reverse ventricular growth (or reverse remodeling) is often considered to be a favorable response to clinical intervention. In recent years, a number of theoretical models have been proposed to model the process of ventricular growth while little has been done to model its reverse. Based on the framework of volumetric strain-driven finite growth with a homeostatic equilibrium range for the elastic myofiber stretch, we propose here a reversible growth model capable of describing both ventricular growth and its reversal. We used this model to construct a semi-analytical solution based on an idealized cylindrical tube model, as well as numerical solutions based on a truncated ellipsoidal model and a human left ventricular model that was reconstructed from magnetic resonance images. We show that our model is able to predict key features in the end-diastolic pressure–volume relationship that were observed experimentally and clinically during ventricular growth and reverse growth. We also show that the residual stress fields generated as a result of differential growth in the cylindrical tube model are similar to those in other nonidentical models utilizing the same geometry.