First evidence of depressed contractility in the border zone of a human myocardial infarction

Bi-ventricular elastic material parameters estimation using 3D CMR myocardial strains in rheumatic heart disease patients

Understanding the elastic material behavior of myocardium during the diastolic phase is critical for evaluating cardiac function and improving treatments for diastolic abnormalities. This study introduces a novel multi-objective optimization framework that incorporates both strain and volume measurements to enhance the accuracy of myocardial property assessments in Rheumatic Heart Disease (RHD) patients and healthy controls. By employing global volume and strain measurements instead of segmented strains from the sixteen AHA regions, we achieve a robust alignment with the Klotz curve across all groups, indicating an accurate simulation of end-diastolic pressure-volume relationships (EDPVRs). Our approach uniquely integrates combinations of longitudinal, circumferential, and radial strains, resulting in an unprecedented reduction in errors between clinical and simulated strain values, with less than one percent difference for targeted parameters. The results demonstrate that the alignment between computational predictions and clinical measurements depends significantly on the choice of optimization target. The study reveals significant differences in tissue mechanics between RHD patients and healthy controls, with notable variations in ventricular stiffness and fiber orientations across optimization targets, confirmed through rigorous statistical analyses. The observed variations in fiber angles, particularly the smaller angles for longitudinal strains and steeper angles for circumferential strains, underscore the intricate relationship between myocardial fiber architecture and cardiac deformation, offering deeper insights into ventricular biomechanics. By presenting qualitative and quantitative differences in stress and strain distributions, this research advances the understanding of myocardial mechanics, highlighting the clinical relevance of fiber orientation and material properties in modeling cardiac mechanics and distinguishing diseased from healthy myocardial behavior.

Effect of excitation sequence of myocardial contraction on the mechanical response of the left ventricle

In the past two decades there has been rapid development in the field of computational cardiac models. These have included either (i) mechanical models that assumed simultaneous myocardial activation, or (ii) electromechanical models that assumed time-varying myocardial activation. The influence of these modelling assumptions of myocardial activation on clinically relevant metrics, like myocardial strain, commonly used for validation of cardiac models has yet to be systematically examined, leading to uncertainty over their influence on the predictions of these models. This study examined the effects of simultaneous (mechanical), uniform endocardial, 3-patch endocardial (simulating the fascicles of the His-Purkinje system) and 1-patch endocardial (simulating the atrioventricular node) excitation sequences on the mechanical response of a synthetic human left ventricular model. The influence of the duration of the activation and time-to-peak contraction was also investigated. The electromechanical and mechanical models produced different strain distributions in early systole. However, these differences decayed as systole progressed. Using the same activation duration (74 ms) the average peak-systolic circumferential strain difference between the models was 0.65±0.37 %. A slightly prolonged activation duration (134 ms) resulted in no substantial difference increase (0.76±0.47 %). Differences up to 3.5 % were observed for prolonged activation durations (200 ms). Endocardial excitation produced non-physiological cumulative activation time distributions compared to the other models. Septal 1-patch excitation resulted in early systolic strain response that resembled pathological left bundle branch block. Decreased time-to-peak contraction exaggerated the effects of electrophysiology. The study found that excitation sequence minimally affects strain distributions at peak systole for physiological and even slightly pathological activation durations. However, electromechanical models with (patho)physiologically informed activation sequences are important for the accurate prediction of early systolic and pathological late systolic responses.

Simulated Effects of Acute Left Ventricular Myocardial Infarction on Mitral Regurgitation in an Ovine Model

Ischemic mitral regurgitation (IMR) occurs from incomplete coaptation of the mitral valve (MV) after myocardial infarction (MI), typically worsened by continued remodeling of the left ventricular (LV). The importance of LV remodeling is clear as IMR is induced by the post-MI dual mechanisms of mitral annular dilation and leaflet tethering from papillary muscle (PM) distension via the MV chordae tendineae (MVCT). However, the detailed etiology of IMR remains poorly understood, in large part due to the complex interactions of the MV and the post-MI LV remodeling processes. Given the patient-specific anatomical complexities of the IMR disease processes, simulation-based approaches represent an ideal approach to improve our understanding of this deadly disease. However, development of patient-specific models of left ventricle-mitral valve (LV-MV) interactions in IMR are complicated by the substantial variability and complexity of the MR etiology itself, making it difficult to extract underlying mechanisms from clinical data alone. To address these shortcomings, we developed a detailed ovine LV-MV finite element (FE) model based on extant comprehensive ovine experimental data. First, an extant ovine LV FE model (Sci. Rep. 2021 Jun 29;11(1):13466) was extended to incorporate the MV using a high fidelity ovine in vivo derived MV leaflet geometry. As it is not currently possible to image the MVCT in vivo, a functionally equivalent MVCT network was developed to create the final LV-MV model. Interestingly, in pilot studies, the MV leaflet strains did not agree well with known in vivo MV leaflet strain fields. We then incorporated previously reported MV leaflet prestrains (J. Biomech. Eng. 2023 Nov 1;145(11):111002) in the simulations. The resulting LV-MV model produced excellent agreement with the known in vivo ovine MV leaflet strains and deformed shapes in the normal state. We then simulated the effects of regional acute infarctions of varying sizes and anatomical locations by shutting down the local myocardial contractility. The remaining healthy (noninfarcted) myocardium mechanical behaviors were maintained, but allowed to adjust their active contractile patterns to maintain the prescribed pressure-volume loop behaviors in the acute post-MI state. For all cases studied, the LV-MV simulation demonstrated excellent agreement with known LV and MV in vivo strains and MV regurgitation orifice areas. Infarct location was shown to play a critical role in resultant MV leaflet strain fields. Specifically, extensional deformations of the posterior leaflets occurred in the posterobasal and laterobasal infarcts, while compressive deformations of the anterior leaflet were observed in the anterobasal infarct. Moreover, the simulated posterobasal infarct induced the largest MV regurgitation orifice area, consistent with experimental observations. The present study is the first detailed LV-MV simulation that reveals the important role of MV leaflet prestrain and functionally equivalent MVCT for accurate predictions of LV-MV interactions. Importantly, the current study further underscored simulation-based methods in understanding MV function as an integral part of the LV.

Heart-on-a-chip systems: disease modeling and drug screening applications

Cardiovascular disease (CVD) is the leading cause of death worldwide, casting a substantial economic footprint and burdening the global healthcare system. Historically, pre-clinical CVD modeling and therapeutic screening have been performed using animal models. Unfortunately, animal models oftentimes fail to adequately mimic human physiology, leading to a poor translation of therapeutics from pre-clinical trials to consumers. Even those that make it to market can be removed due to unforeseen side effects. As such, there exists a clinical, technological, and economical need for systems that faithfully capture human (patho)physiology for modeling CVD, assessing cardiotoxicity, and evaluating drug efficacy. Heart-on-a-chip (HoC) systems are a part of the broader organ-on-a-chip paradigm that leverages microfluidics, tissue engineering, microfabrication, electronics, and gene editing to create human-relevant models for studying disease, drug-induced side effects, and therapeutic efficacy. These compact systems can be capable of real-time measurements and on-demand characterization of tissue behavior and could revolutionize the drug development process. In this review, we highlight the key components that comprise a HoC system followed by a review of contemporary reports of their use in disease modeling, drug toxicity and efficacy assessment, and as part of multi-organ-on-a-chip platforms. We also discuss future perspectives and challenges facing the field, including a discussion on the role that standardization is expected to play in accelerating the widespread adoption of these platforms.

The dependency of fetal left ventricular biomechanics function on myocardium helix angle configuration

The helix angle configuration of the myocardium is understood to contribute to the heart function, as finite element (FE) modeling of postnatal hearts showed that altered configurations affected cardiac function and biomechanics. However, similar investigations have not been done on the fetal heart. To address this, we performed image-based FE simulations of fetal left ventricles (LV) over a range of helix angle configurations, assuming a linear variation of helix angles from epicardium to endocardium. Results showed that helix angles have substantial influence on peak myofiber stress, cardiac stroke work, myocardial deformational burden, and spatial variability of myocardial strain. A good match between LV myocardial strains from FE simulations to those measured from 4D fetal echo images could only be obtained if the transmural variation of helix angle was generally between 110 and 130°, suggesting that this was the physiological range. Experimentally discovered helix angle configurations from the literature were found to produce high peak myofiber stress, high cardiac stroke work, and a low myocardial deformational burden, but did not coincide with configurations that would optimize these characteristics. This may suggest that the fetal development of myocyte orientations depends concurrently on several factors rather than a single factor. We further found that the shape, rather than the size of the LV, determined the manner at which helix angles influenced these characteristics, as this influence changed significantly when the LV shape was varied, but not when a heart was scaled from fetal to adult size while retaining the same shape. This may suggest that biomechanical optimality would be affected during diseases that altered the geometric shape of the LV.

A myocardial infarct border-zone-on-a-chip demonstrates distinct regulation of cardiac tissue function by an oxygen gradient

After a myocardial infarction, the boundary between the injured, hypoxic tissue and the adjacent viable, normoxic tissue, known as the border zone, is characterized by an oxygen gradient. Yet, the impact of an oxygen gradient on cardiac tissue function is poorly understood, largely due to limitations of existing experimental models. Here, we engineered a microphysiological system to controllably expose engineered cardiac tissue to an oxygen gradient that mimics the border zone and measured the effects of the gradient on electromechanical function and the transcriptome. The gradient delayed calcium release, reuptake, and propagation; decreased diastolic and peak systolic stress; and increased expression of inflammatory cascades that are hallmarks of myocardial infarction. These changes were distinct from those observed in tissues exposed to uniform normoxia or hypoxia, demonstrating distinct regulation of cardiac tissue phenotypes by an oxygen gradient. Our border-zone-on-a-chip model advances functional and mechanistic insight into oxygen-dependent cardiac tissue pathophysiology.

Left Ventricle Biomechanics of Low-Flow, Low-Gradient Aortic Stenosis: A Patient-Specific Computational Model

This study aimed to create an imaging-derived patient-specific computational model of low-flow, low-gradient (LFLG) aortic stenosis (AS) to obtain biomechanics data about the left ventricle. LFLG AS is now a commonly recognized sub-type of aortic stenosis. There remains much controversy over its management, and investigation into ventricular biomechanics may elucidate pathophysiology and better identify patients for valve replacement. ECG-gated cardiac computed tomography images from a patient with LFLG AS were obtained to provide patient-specific geometry for the computational model. Surfaces of the left atrium, left ventricle (LV), and outflow track were segmented. A previously validated multi-scale, multi-physics computational human heart model was adapted to the patient-specific geometry, yielding a model consisting of 91,000 solid elements. This model was coupled to a virtual circulatory system and calibrated to clinically measured parameters from echocardiography and cardiac catheterization data. The simulation replicated key physiologic parameters within 10% of their clinically measured values. Global LV systolic myocardial stress was 7.1 ± 1.8 kPa. Mean stress of the basal, middle, and apical segments were 7.7 ± 1.8 kPa, 9.1 ± 3.8 kPa, and 6.4 ± 0.4 kPa, respectively. This is the first patient-specific computational model of LFLG AS based on clinical imaging. Low myocardial stress correlated with low ejection fraction and eccentric LV remodeling. Further studies are needed to understand how alterations in LV biomechanics correlates with clinical outcomes of AS.

Effects of dispersed fibres in myocardial mechanics, Part II: active response

This work accompanies the first part of our study "effects of dispersed fibres in myocardial mechanics: Part I passive response" with a focus on myocardial active contraction. Existing studies have suggested that myofibre architecture plays an important role in myocardial active contraction. Following the first part of our study, we firstly study how the general fibre architecture affects ventricular pump function by varying the mean myofibre rotation angles, and then the impact of fibre dispersion along the myofibre direction on myocardial contraction in a left ventricle model. Dispersed active stress is described by a generalised structure tensor method for its computational efficiency. Our results show that both the myofibre rotation angle and its dispersion can significantly affect cardiac pump function by redistributing active tension circumferentially and longitudinally. For example, larger myofibre rotation angle and higher active tension along the sheet-normal direction can lead to much reduced end-systolic volume and higher longitudinal shortening, and thus a larger ejection fraction. In summary, these two studies together have demonstrated that it is necessary and essential to include realistic fibre structures (both fibre rotation angle and fibre dispersion) in personalised cardiac modelling for accurate myocardial dynamics prediction.

The role of regional myocardial topography post-myocardial infarction on infarct extension

Infarct extension involves necrosis of healthy myocardium in the border zone (BZ), progressively enlarging the infarct zone (IZ) and recruiting the remote zone (RZ) into the BZ, eventually leading to heart failure. The mechanisms underlying infarct extension remain unclear, but myocyte stretching has been suggested as the most likely cause. Using human patient-specific left-ventricular (LV) numerical simulations established from cardiac magnetic resonance imaging (MRI) of myocardial infarction (MI) patients, the correlation between infarct extension and regional mechanics abnormality was investigated by analysing the fibre stress-strain loops (FSSLs). FSSL abnormality was characterised using the directional regional external work (DREW) index, which measures FSSL area and loop direction. Sensitivity studies were also performed to investigate the effect of infarct stiffness on regional myocardial mechanics and potential for infarct extension. We found that infarct extension was correlated to severely abnormal FSSL in the form of counter-clockwise loop at the RZ close to the infarct, as indicated by negative DREW values. In regions demonstrating negative DREW values, we observed substantial fibre stretching in the isovolumic relaxation (IVR) phase accompanied by a reduced rate of systolic shortening. Such stretching in IVR phase in part of the RZ was due to its inability to withstand the high LV pressure that was still present and possibly caused by regional myocardial stiffness inhomogeneity. Further analysis revealed that the occurrence of severely abnormal FSSL due to IVR fibre stretching near the RZ-BZ boundary was due to a large amount of surrounding infarcted tissue, or an excessively stiff IZ.

The impact of myocardial compressibility on organ-level simulations of the normal and infarcted heart

Myocardial infarction (MI) rapidly impairs cardiac contractile function and instigates maladaptive remodeling leading to heart failure. Patient-specific models are a maturing technology for developing and determining therapeutic modalities for MI that require accurate descriptions of myocardial mechanics. While substantial tissue volume reductions of 15-20% during systole have been reported, myocardium is commonly modeled as incompressible. We developed a myocardial model to simulate experimentally-observed systolic volume reductions in an ovine model of MI. Sheep-specific simulations of the cardiac cycle were performed using both incompressible and compressible tissue material models, and with synchronous or measurement-guided contraction. The compressible tissue model with measurement-guided contraction gave best agreement with experimentally measured reductions in tissue volume at peak systole, ventricular kinematics, and wall thickness changes. The incompressible model predicted myofiber peak contractile stresses approximately double the compressible model (182.8 kPa, 107.4 kPa respectively). Compensatory changes in remaining normal myocardium with MI present required less increase of contractile stress in the compressible model than the incompressible model (32.1%, 53.5%, respectively). The compressible model therefore provided more accurate representation of ventricular kinematics and potentially more realistic computed active contraction levels in the simulated infarcted heart. Our findings suggest that myocardial compressibility should be incorporated into future cardiac models for improved accuracy.

Novel Methodology for Measuring Regional Myocardial Efficiency

Our approach differs from the usual global measure of cardiac efficiency by using PET/MRI to measure efficiency of small pieces of cardiac tissue whose limiting size is equal to the spatial resolution of the PET scanner. We initiated a dynamic cardiac PET study immediately prior to the injection of 15.1 mCi of 11C-acetate acquiring data for 25 minutes while simultaneously acquiring MRI cine data. 1) A 3D finite element (FE) biomechanical model of the imaged heart was constructed by utilizing nonrigid deformable image registration to alter the Dassault Systèmes FE Living Heart Model (LHM) to fit the geometry in the cardiac MRI cine data. The patient specific FE cardiac model with estimates of stress, strain, and work was transformed into PET/MRI format. 2) A 1-tissue compartment model was used to calculate wash-in (K1) and the linear portion of the decay in the PET 11C-acetate time activity curve (TAC) was used to calculate the wash-out k2(mono) rate constant. K1 was used to calculate blood flow and k2(mono) was used to calculate myocardial volume oxygen consumption ( MVO2 ). 3) Estimates of stress and strain were used to calculate Myocardial Equivalent Minute Work ( MEMW ) and Cardiac Efficiency = MEMW/MVO2 was then calculated for 17 tissue segments of the left ventricle. The global MBF was 0.96 ± 0.15 ml/min/gm and MVO2 ranged from 8 to 17 ml/100gm/min. Six central slices of the MRI cine data provided a range of MEMW of 0.1 to 0.4 joules/gm/min and a range of Cardiac Efficiency of 6 to 18%.

Biomechanics of Human Fetal Hearts with Critical Aortic Stenosis

Critical aortic stenosis (AS) of the fetal heart causes a drastic change in the cardiac biomechanical environment. Consequently, a substantial proportion of such cases will lead to a single-ventricular birth outcome. However, the biomechanics of the disease is not well understood. To address this, we performed Finite Element (FE) modelling of the healthy fetal left ventricle (LV) based on patient-specific 4D ultrasound imaging, and simulated various disease features observed in clinical fetal AS to understand their biomechanical impact. These features included aortic stenosis, mitral regurgitation (MR) and LV hypertrophy, reduced contractility, and increased myocardial stiffness. AS was found to elevate LV pressures and myocardial stresses, and depending on severity, can drastically decrease stroke volume and myocardial strains. These effects are moderated by MR. AS alone did not lead to MR velocities above 3 m/s unless LV hypertrophy was included, suggesting that hypertrophy may be involved in clinical cases with high MR velocities. LV hypertrophy substantially elevated LV pressure, valve flow velocities and stroke volume, while reducing LV contractility resulted in diminished LV pressure, stroke volume and wall strains. Typical extent of hypertrophy during fetal AS in the clinic, however, led to excessive LV pressure and valve velocity in the FE model, suggesting that reduced contractility is typically associated with hypertrophy. Increased LV passive stiffness, which might represent fibroelastosis, was found to have minimal impact on LV pressures, stroke volume, and wall strain. This suggested that fibroelastosis could be a by-product of the disease progression and does not significantly impede cardiac function. Our study demonstrates that FE modelling is a valuable tool for elucidating the biomechanics of congenital heart disease and can calculate parameters which are difficult to measure, such as intraventricular pressure and myocardial stresses.

In vivo estimation of normal left ventricular stiffness and contractility based on routine cine MR acquisition

Post-myocardial infarction remodeling process is known to alter the mechanical properties of the heart. Biomechanical parameters, such as tissue stiffness and contractility, would be useful for clinicians to better assess the severity of the diseased heart. However, these parameters are difficult to obtain in the current clinical practice. In this paper, we estimated subject-specific in vivo myocardial stiffness and contractility from 21 healthy volunteers, based on left ventricle models constructed from data acquired from routine cardiac MR acquisition only. The subject-specific biomechanical parameters were quantified using an inverse finite-element modelling approach. The personalized models were evaluated against relevant clinical metrics extracted from the MR data, such as circumferential strain, wall thickness and fractional thickening. We obtained the ranges of healthy biomechanical indices of 1.60 ± 0.22 kPa for left ventricular stiffness and 95.13 ± 14.56 kPa for left ventricular contractility. These reference normal values can be used for future model-based investigation on the stiffness and contractility of ischemic myocardium.

Upregulated lncRNA Pvt1 may be important for cardiac remodeling at the infarct border zone

Myocardial infarction (MI) is a leading cause of mortality due to progression to ventricular arrhythmias (VAs) or heart failure (HF). Cardiac remodeling at the infarct border zone (IBZ) is the primary contributor for VAs or HF. Therefore, genes involved in IBZ remodeling may be potential targets for the treatment of MI, but the mechanism remains unclear. The present study aimed to explain the molecular mechanisms of IBZ remodeling based on the roles of long non-coding RNAs (lncRNAs). After downloading miRNA (GSE76592) and mRNA/lncRNA (GSE52313) datasets from the Gene Expression Omnibus database, 23 differentially expressed miRNAs (DEMs), 2,563 genes (DEGs) and 168 lncRNAs (DELs) were identified between IBZ samples of MI mice and sham controls. A total of 483 DEGs were predicted to be regulated by 23 DEMs, among which Itgam, Met and TNF belonged to hub genes after five topological parameters were calculated for genes in the protein-protein interaction network. These hub genes-associated DEMs (mmu-miR-181a, mmu-miR-762) can also interact with six DELs (Gm15832, Gas5, Gm6634, Pvt1, Gm14636 and A330023F24Rik) to constitute the competing endogenous RNA (ceRNA) axes. Furthermore, a co-expression network was constructed based on the co-expression pairs between 44 DELs and 297 DEGs, in which Pvt1 and Bst1 were overlapped with the ceRNA network. Thus, Bst1-associated ceRNA (Pvt1-mmu-miR-181a-Bst1) and co-expression (Pvt-Bst1) axes were also pivotal for MI. Accordingly, Pvt1 may be a crucial lncRNA for modification of cardiac remodeling in the IBZ after MI and may function by acting as a ceRNA for miR-181a to regulate TNF/Met/Itgam/Bst1 or by co-expressing with Bst1.

Intra-myocardial alginate hydrogel injection acts as a left ventricular mid-wall constraint in swine

Despite positive initial outcomes emerging from preclinical and early clinical investigation of alginate hydrogel injection therapy as a treatment for heart failure, the lack of knowledge about the mechanism of action remains a major shortcoming that limits the efficacy of treatment design. To identify the mechanism of action, we examined previously unobtainable measurements of cardiac function from in vivo, ex vivo, and in silico states of clinically relevant heart failure (HF) in large animals. High-resolution ex vivo magnetic resonance imaging and histological data were used along with state-of-the-art subject-specific computational model simulations. Ex vivo data were incorporated in detailed geometric computational models for swine hearts in health (n = 5), ischemic HF (n = 5), and ischemic HF treated with alginate hydrogel injection therapy (n = 5). Hydrogel injection therapy mitigated elongation of sarcomere lengths (1.68 ± 0.10μm [treated] vs. 1.78 ± 0.15μm [untreated], p<0.001). Systolic contractility in treated animals improved substantially (ejection fraction = 43.9 ± 2.8% [treated] vs. 34.7 ± 2.7% [untreated], p<0.01). The in silico models realistically simulated in vivo function with >99% accuracy and predicted small myofiber strain in the vicinity of the solidified hydrogel that was sustained for up to 13 mm away from the implant. These findings suggest that the solidified alginate hydrogel material acts as an LV mid-wall constraint that significantly reduces adverse LV remodeling compared to untreated HF controls without causing negative secondary outcomes to cardiac function. STATEMENT OF SIGNIFICANCE: Heart failure is considered a growing epidemic and hence an important health problem in the US and worldwide. Its high prevalence (5.8 million and 23 million, respectively) is expected to increase by 25% in the US alone by 2030. Heart failure is associated with high morbidity and mortality, has a 5-year mortality rate of 50%, and contributes considerably to the overall cost of health care ($53.1 billion in the US by 2030). Despite positive initial outcomes emerging from preclinical and early clinical investigation of alginate hydrogel injection therapy as a treatment for heart failure, the lack of knowledge concerning the mechanism of action remains a major shortcoming that limits the efficacy of treatment design. To understand the mechanism of action, we combined high-resolution ex vivo magnetic resonance imaging and histological data in swine with state-of-the-art subject-specific computational model simulations. The in silico models realistically simulated in vivo function with >99% accuracy and predicted small myofiber strain in the vicinity of the solidified hydrogel that was sustained for up to 13 mm away from the implant. These findings suggest that the solidified alginate hydrogel material acts as a left ventricular mid-wall constraint that significantly reduces adverse LV remodeling compared to untreated heart failure controls without causing negative secondary outcomes to cardiac function. Moreover, if the hydrogel can be delivered percutaneously rather than via the currently used open-chest procedure, this therapy may become routine for heart failure treatment. A minimally invasive procedure would be in the best interest of this patient population; i.e., one that cannot tolerate general anesthesia and surgery, and it would be significantly more cost-effective than surgery.

Effect of myofibre architecture on ventricular pump function by using a neonatal porcine heart model: from DT-MRI to rule-based methods

Myofibre architecture is one of the essential components when constructing personalized cardiac models. In this study, we develop a neonatal porcine bi-ventricle model with three different myofibre architectures for the left ventricle (LV). The most realistic one is derived from ex vivo diffusion tensor magnetic resonance imaging, and other two simplifications are based on rule-based methods (RBM): one is regionally dependent by dividing the LV into 17 segments, each with different myofibre angles, and the other is more simplified by assigning a set of myofibre angles across the whole ventricle. Results from different myofibre architectures are compared in terms of cardiac pump function. We show that the model with the most realistic myofibre architecture can produce larger cardiac output, higher ejection fraction and larger apical twist compared with those of the rule-based models under the same pre/after-loads. Our results also reveal that when the cross-fibre contraction is included, the active stress seems to play a dual role: its sheet-normal component enhances the ventricular contraction while its sheet component does the opposite. We further show that by including non-symmetric fibre dispersion using a general structural tensor, even the most simplified rule-based myofibre model can achieve similar pump function as the most realistic one, and cross-fibre contraction components can be determined from this non-symmetric dispersion approach. Thus, our study highlights the importance of including myofibre dispersion in cardiac modelling if RBM are used, especially in personalized models.

Method for Calibration of Left Ventricle Material Properties using 3D Echocardiography Endocardial Strains

We sought to calibrate mechanical properties of left ventricle (LV) based on 3D speckle tracking echocardiographic imaging data recorded from 16 segments defined by American Heart Association (AHA). The in vivo data were used to create finite element (FE) LV and biventricular (BV) models. The orientation of the fibers in the LV model was rule-based, but diffusion tensor magnetic resonance imaging (DT-MRI) data were used for the fiber directions in the BV model. A nonlinear fiber-reinforced constitutive equation was used to describe the passive behavior of the myocardium, whereas the active tension was described by a model based on tissue contraction (Tmax). Isight was used for optimization, which used Abaqus as the forward solver (Simulia, Providence, USA). The calibration of passive properties based on the end diastolic pressure volume relation (ED PVR) curve resulted in relatively good agreement (mean error = -0.04 ml). The difference between experimental and computational strains decreased after segmental strain metrics, rather than global metrics, were used for calibration: for the LV model, the mean difference reduced from 0.129 to 0.046 (circumferential) and from 0.076 to 0.059 (longitudinal); for the BV model, the mean difference nearly did not change in the circumferential direction (0.061) but reduced in the longitudinal direction from 0.076 to 0.055. The calibration of mechanical properties for myocardium can be improved using segmental strain metrics. The importance of realistic fiber orientation and geometry for modeling of the LV was shown.

Computational Modelling of the Effect of Infarct Stiffness on Regional Myocardial Mechanics

Ventricular remodeling after myocardial infarction increases the rate of mortality and is highly associated with the extent of infarct transmurality. It is hypothesized that infarct stiffness alters regional mechanics and affects the likelihood of human ventricular remodeling. However, this is yet to be studied in detail. In this paper, we present simulations from an actively-contracting left ventricular model to investigate the effects of transmural infarct stiffness on myofiber regional mechanics. Results show that higher infarct stiffness reduces systolic stress at the infarct and border zones, minimizing infarct bulging but increasing the diastolic stress at the endocardial border zone. Determining a proper amount of infarct stiffness is required to achieve a balanced regional mechanics across the cardiac cycle that may be useful in therapy, such as myocardial hydrogel injection to adjust its stiffness and reduce stress to prevent ventricular remodeling.

Use of Multifactorial Treatments to Address the Challenge of Translating Experimental Myocardial Infarct Reduction Strategies

Myocardial tissue damage that occurs during an ischemic event leads to a spiraling deterioration of cardiac muscle structural and functional integrity. Reperfusion is the only known efficacious strategy and is the most commonly used treatment to reduce injury and prevent remodeling. However, timing is critical, and the procedure is not always feasible for a variety of reasons. The complex molecular basis for cardioprotection has been studied for decades but formulation of a viable therapeutic that can significantly attenuate myocardial injury remains elusive. In this review, we address barriers to the development of a fruitful approach that will substantially improve the prognosis of those suffering from this widespread and largely unmitigated disease. Furthermore, we proffer that ephrinA1, a candidate molecule that satisfies many of the important criteria discussed, possesses robust potential to overcome these hurdles and thus offers protection that surpasses the limitations currently observed.

Organ-level validation of a cross-bridge cycling descriptor in a left ventricular finite element model: effects of ventricular loading on myocardial strains

Although detailed cell‐based descriptors of cross‐bridge cycling have been applied in finite element (FE) heart models to describe ventricular mechanics, these multiscale models have never been tested rigorously to determine if these descriptors, when scaled up to the organ‐level, are able to reproduce well‐established organ‐level physiological behaviors. To address this void, we here validate a left ventricular (LV) FE model that is driven by a cell‐based cross‐bridge cycling descriptor against key organ‐level heart physiology. The LV FE model was coupled to a closed‐loop lumped parameter circulatory model to simulate different ventricular loading conditions (preload and afterload) and contractilities. We show that our model is able to reproduce a linear end‐systolic pressure volume relationship, a curvilinear end‐diastolic pressure volume relationship and a linear relationship between myocardial oxygen consumption and pressure–volume area. We also show that the validated model can predict realistic LV strain‐time profiles in the longitudinal, circumferential, and radial directions. The predicted strain‐time profiles display key features that are consistent with those measured in humans, such as having similar peak strains, time‐to‐peak‐strain, and a rapid change in strain during atrial contraction at late‐diastole. Our model shows that the myocardial strains are sensitive to not only LV contractility, but also to the LV loading conditions, especially to a change in afterload. This result suggests that caution must be exercised when associating changes in myocardial strain with changes in LV contractility. The methodically validated multiscale model will be used in future studies to understand human heart diseases.

Construction and Validation of Subject-Specific Biventricular Finite-Element Models of Healthy and Failing Swine Hearts From High-Resolution DT-MRI

Predictive computational modeling has revolutionized classical engineering disciplines and is in the process of transforming cardiovascular research. This is particularly relevant for investigating emergent therapies for heart failure, which remains a leading cause of death globally. The creation of subject-specific biventricular computational cardiac models has been a long-term endeavor within the biomedical engineering community. Using high resolution (0.3 × 0.3 × 0.8 mm) ex vivo data, we constructed a precise fully subject-specific biventricular finite-element model of healthy and failing swine hearts. Each model includes fully subject-specific geometries, myofiber architecture and, in the case of the failing heart, fibrotic tissue distribution. Passive and active material properties are prescribed using hyperelastic strain energy functions that define a nearly incompressible, orthotropic material capable of contractile function. These materials were calibrated using a sophisticated multistep approach to match orthotropic tri-axial shear data as well as subject-specific hemodynamic ventricular targets for pressure and volume to ensure realistic cardiac function. Each mechanically beating heart is coupled with a lumped-parameter representation of the circulatory system, allowing for a closed-loop definition of cardiovascular flow. The circulatory model incorporates unidirectional fluid exchanges driven by pressure gradients of the model, which in turn are driven by the mechanically beating heart. This creates a computationally meaningful representation of the dynamic beating of the heart coupled with the circulatory system. Each model was calibrated using subject-specific experimental data and compared with independent in vivo strain data obtained from echocardiography. Our methods produced highly detailed representations of swine hearts that function mechanically in a remarkably similar manner to the in vivo subject-specific strains on a global and regional comparison. The degree of subject-specificity included in the models represents a milestone for modeling efforts that captures realism of the whole heart. This study establishes a foundation for future computational studies that can apply these validated methods to advance cardiac mechanics research.

Evaluation of a Novel Finite Element Model of Active Contraction in the Heart

Finite element (FE) modeling is becoming a widely used approach for the investigation of global heart function. In the present study, a novel model of cellular-level systolic contraction, which includes both length- and velocity-dependence, was implemented into a 3D non-linear FE code. To validate this new FE implementation, an optimization procedure was used to determine the contractile parameters, associated with sarcomeric function, by comparing FE-predicted pressure and strain to experimental measures collected with magnetic resonance imaging and catheterization in the ventricles of five healthy rats. The pressure-volume relationship generated by the FE models matched well with the experimental data. Additionally, the regional distribution of end-systolic strains and circumferential-longitudinal shear angle exhibited good agreement with experimental results overall, with the main deviation occurring in the septal region. Moreover, the FE model predicted a heterogeneous distribution of sarcomere re-lengthening after ventricular ejection, which is consistent with previous in vivo studies. In conclusion, the new FE active contraction model was able to predict the global performance and regional mechanical behaviors of the LV during the entire cardiac cycle. By including more accurate cellular-level mechanisms, this model could provide a better representation of the LV and enhance cardiac research related to both systolic and diastolic dysfunction.

Changes and classification in myocardial contractile function in the left ventricle following acute myocardial infarction

In this research, we hypothesized that novel biomechanical parameters are discriminative in patients following acute ST-segment elevation myocardial infarction (STEMI). To identify these biomechanical biomarkers and bring computational biomechanics ‘closer to the clinic’, we applied state-of-the-art multiphysics cardiac modelling combined with advanced machine learning and multivariate statistical inference to a clinical database of myocardial infarction. We obtained data from 11 STEMI patients (ClinicalTrials.gov NCT01717573) and 27 healthy volunteers, and developed personalized mathematical models for the left ventricle (LV) using an immersed boundary method. Subject-specific constitutive parameters were achieved by matching to clinical measurements. We have shown, for the first time, that compared with healthy controls, patients with STEMI exhibited increased LV wall active tension when normalized by systolic blood pressure, which suggests an increased demand on the contractile reserve of remote functional myocardium. The statistical analysis reveals that the required patient-specific contractility, normalized active tension and the systolic myofilament kinematics have the strongest explanatory power for identifying the myocardial function changes post-MI. We further observed a strong correlation between two biomarkers and the changes in LV ejection fraction at six months from baseline (the required contractility (r = − 0.79, p < 0.01) and the systolic myofilament kinematics (r = 0.70, p = 0.02)). The clinical and prognostic significance of these biomechanical parameters merits further scrutinization.

Computational Investigation of Transmural Differences in Left Ventricular Contractility

Myocardial contractility of the left ventricle (LV) plays an essential role in maintaining normal pump function. A recent ex vivo experimental study showed that cardiomyocyte force generation varies across the three myocardial layers of the LV wall. However, the in vivo distribution of myocardial contractile force is still unclear. The current study was designed to investigate the in vivo transmural distribution of myocardial contractility using a noninvasive computational approach. For this purpose, four cases with different transmural distributions of maximum isometric tension (Tmax) and/or reference sarcomere length (lR) were tested with animal-specific finite element (FE) models, in combination with magnetic resonance imaging (MRI), pressure catheterization, and numerical optimization. Results of the current study showed that the best fit with in vivo MRI-derived deformation was obtained when Tmax assumed different values in the subendocardium, midmyocardium, and subepicardium with transmurally varying lR. These results are consistent with recent ex vivo experimental studies, which showed that the midmyocardium produces more contractile force than the other transmural layers. The systolic strain calculated from the best-fit FE model was in good agreement with MRI data. Therefore, the proposed noninvasive approach has the capability to predict the transmural distribution of myocardial contractility. Moreover, FE models with a nonuniform distribution of myocardial contractility could provide a better representation of LV function and be used to investigate the effects of transmural changes due to heart disease.

Advances in computational modelling for personalised medicine after myocardial infarction

Myocardial infarction (MI) is a leading cause of premature morbidity and mortality worldwide. Determining which patients will experience heart failure and sudden cardiac death after an acute MI is notoriously difficult for clinicians. The extent of heart damage after an acute MI is informed by cardiac imaging, typically using echocardiography or sometimes, cardiac magnetic resonance (CMR). These scans provide complex data sets that are only partially exploited by clinicians in daily practice, implying potential for improved risk assessment. Computational modelling of left ventricular (LV) function can bridge the gap towards personalised medicine using cardiac imaging in patients with post-MI. Several novel biomechanical parameters have theoretical prognostic value and may be useful to reflect the biomechanical effects of novel preventive therapy for adverse remodelling post-MI. These parameters include myocardial contractility (regional and global), stiffness and stress. Further, the parameters can be delineated spatially to correspond with infarct pathology and the remote zone. While these parameters hold promise, there are challenges for translating MI modelling into clinical practice, including model uncertainty, validation and verification, as well as time-efficient processing. More research is needed to (1) simplify imaging with CMR in patients with post-MI, while preserving diagnostic accuracy and patient tolerance (2) to assess and validate novel biomechanical parameters against established prognostic biomarkers, such as LV ejection fraction and infarct size. Accessible software packages with minimal user interaction are also needed. Translating benefits to patients will be achieved through a multidisciplinary approach including clinicians, mathematicians, statisticians and industry partners.

Estimating prognosis in patients with acute myocardial infarction using personalized computational heart models

Biomechanical computational models have potential prognostic utility in patients after an acute ST-segment-elevation myocardial infarction (STEMI). In a proof-of-concept study, we defined two groups (1) an acute STEMI group (n = 6, 83% male, age 54 ± 12 years) complicated by left ventricular (LV) systolic dysfunction; (2) an age- and sex- matched hyper-control group (n = 6, 83% male, age 46 ± 14 years), no prior history of cardiovascular disease and normal systolic blood pressure (SBP < 130 mmHg). Cardiac MRI was performed in the patients (2 days & 6 months post-STEMI) and the volunteers, and biomechanical heart models were synthesized for each subject. The candidate parameters included normalized active tension (AT norm) and active tension at the resting sarcomere length (T req, reflecting required contractility). Myocardial contractility was inversely determined from personalized heart models by matching CMR-imaged LV dynamics. Compared with controls, patients with recent STEMI exhibited increased LV wall active tension when normalized by SBP. We observed a linear relationship between T req 2 days post-MI and global longitudinal strain 6 months later (r = 0.86; p = 0.03). T req may be associated with changes in LV function in the longer term in STEMI patients complicated by LV dysfunction. Further studies seem warranted.

Injectable Bioengineered Hydrogel Therapy in the Treatment of Ischemic Cardiomyopathy

OPINION STATEMENT

Over the past two decades, the field of cardiovascular medicine has seen the rapid development of multiple different modalities for the treatment of ischemic myocardial disease. Most research efforts have focused on strategies aimed at coronary revascularization, with significant technological advances made in percutaneous coronary interventions as well as coronary artery bypass graft surgery. However, recent research efforts have shifted towards ways to address the downstream effects of myocardial infarction on both cellular and molecular levels. To this end, the broad application of injectable hydrogel therapy after myocardial infarction has stimulated tremendous interest. In this article, we will review what hydrogels are, how they can be bioengineered in unique ways to optimize therapeutic potential, and how they can be used as part of a treatment strategy after myocardial infarction.

The role of infarct transmural extent in infarct extension: A computational study

Infarct extension, a process involving progressive extension of the infarct zone (IZ) into the normally perfused border zone (BZ), leads to continuous degradation of the myocardial function and adverse remodelling. Despite carrying a high risk of mortality, detailed understanding of the mechanisms leading to BZ hypoxia and infarct extension remains unexplored. In the present study, we developed a 3D truncated ellipsoidal left ventricular model incorporating realistic electromechanical properties and fibre orientation to examine the mechanical interaction among the remote, infarct and BZs in the presence of varying infarct transmural extent (TME). Localized highly abnormal systolic fibre stress was observed at the BZ, owing to the simultaneous presence of moderately increased stiffness and fibre strain at this region, caused by the mechanical tethering effect imposed by the overstretched IZ. Our simulations also demonstrated the greatest tethering effect and stress in BZ regions with fibre direction tangential to the BZ-remote zone boundary. This can be explained by the lower stiffness in the cross-fibre direction, which gave rise to a greater stretching of the IZ in this direction. The average fibre strain of the IZ, as well as the maximum stress in the sub-endocardial layer, increased steeply from 10% to 50% infarct TME, and slower thereafter. Based on our stress-strain loop analysis, we found impairment in the myocardial energy efficiency and elevated energy expenditure with increasing infarct TME, which we believe to place the BZ at further risk of hypoxia. Copyright © 2016 John Wiley & Sons, Ltd.

Effects of using the unloaded configuration in predicting the in vivo diastolic properties of the heart

Computational models are increasingly being used to investigate the mechanical properties of cardiac tissue. While much insight has been gained from these studies, one important limitation associated with computational modeling arises when using in vivo images of the heart to generate the reference state of the model. An unloaded reference configuration is needed to accurately represent the deformation of the heart. However, it is rare for a beating heart to actually reach a zero-pressure state during the cardiac cycle. To overcome this, a computational technique was adapted to determine the unloaded configuration of an in vivo porcine left ventricle (LV). In the current study, in vivo measurements were acquired using magnetic resonance images (MRI) and synchronous pressure catheterization in the LV (N = 5). The overall goal was to quantify the effects of using early-diastolic filling as the reference configuration (common assumption used in modeling) versus using the unloaded reference configuration for predicting the in vivo properties of LV myocardium. This was accomplished by using optimization to minimize the difference between MRI measured and finite element predicted strains and cavity volumes. The results show that when using the unloaded reference configuration, the computational method predicts material properties for LV myocardium that are softer and less anisotropic than when using the early-diastolic filling reference configuration. This indicates that the choice of reference configuration could have a significant impact on capturing the realistic mechanical response of the heart.

Fibrous scaffolds for building hearts and heart parts

Extracellular matrix (ECM) structure and biochemistry provide cell-instructive cues that promote and regulate tissue growth, function, and repair. From a structural perspective, the ECM is a scaffold that guides the self-assembly of cells into distinct functional tissues. The ECM promotes the interaction between individual cells and between different cell types, and increases the strength and resilience of the tissue in mechanically dynamic environments. From a biochemical perspective, factors regulating cell-ECM adhesion have been described and diverse aspects of cell-ECM interactions in health and disease continue to be clarified. Natural ECMs therefore provide excellent design rules for tissue engineering scaffolds. The design of regenerative three-dimensional (3D) engineered scaffolds is informed by the target ECM structure, chemistry, and mechanics, to encourage cell infiltration and tissue genesis. This can be achieved using nanofibrous scaffolds composed of polymers that simultaneously recapitulate 3D ECM architecture, high-fidelity nanoscale topography, and bio-activity. Their high porosity, structural anisotropy, and bio-activity present unique advantages for engineering 3D anisotropic tissues. Here, we use the heart as a case study and examine the potential of ECM-inspired nanofibrous scaffolds for cardiac tissue engineering. We asked: Do we know enough to build a heart? To answer this question, we tabulated structural and functional properties of myocardial and valvular tissues for use as design criteria, reviewed nanofiber manufacturing platforms and assessed their capabilities to produce scaffolds that meet our design criteria. Our knowledge of the anatomy and physiology of the heart, as well as our ability to create synthetic ECM scaffolds have advanced to the point that valve replacement with nanofibrous scaffolds may be achieved in the short term, while myocardial repair requires further study in vitro and in vivo.

MRI evaluation of injectable hyaluronic acid-based hydrogel therapy to limit ventricular remodeling after myocardial infarction

Injectable biomaterials are an attractive therapy to attenuate left ventricular (LV) remodeling after myocardial infarction (MI). Although studies have shown that injectable hydrogels improve cardiac structure and function in vivo, temporal changes in infarct material properties after treatment have not been assessed. Emerging imaging and modeling techniques now allow for serial, non-invasive estimation of infarct material properties. Specifically, cine magnetic resonance imaging (MRI) assesses global LV structure and function, late-gadolinium enhancement (LGE) MRI enables visualization of infarcted tissue to quantify infarct expansion, and spatial modulation of magnetization (SPAMM) tagging provides passive wall motion assessment as a measure of tissue strain, which can all be used to evaluate infarct properties when combined with finite element (FE) models. In this work, we investigated the temporal effects of degradable hyaluronic acid (HA) hydrogels on global LV remodeling, infarct thinning and expansion, and infarct stiffness in a porcine infarct model for 12 weeks post-MI using MRI and FE modeling. Hydrogel treatment led to decreased LV volumes, improved ejection fraction, and increased wall thickness when compared to controls. FE model simulations demonstrated that hydrogel therapy increased infarct stiffness for 12 weeks post-MI. Thus, evaluation of myocardial tissue properties through MRI and FE modeling provides insight into the influence of injectable hydrogel therapies on myocardial structure and function post-MI.

Electromechanics modeling of the effects of myocardial infarction on left ventricular function

Ventricular remodeling may occur following myocardial infarction of the left ventricle (LV) and such remodeling has been shown to be correlated with increased patient morbidity and mortality. It is thus important to estimate the likelihood of remodeling from the state of the infarcted LV. In this paper, we present simulations from an actively-contracting truncated ellipsoid LV model, incorporating realistic fiber orientation and electromechanical properties, to investigate the effects of infarct size and transmural extent (TME) on myofiber regional mechanics. Results showed that transmural infarcts greatly elevated both the myofiber stress and strain at the border zone during end systole, making the LV more susceptible to structural remodeling. It was found that TME rather than infarct size was more predictive of LV remodeling.

A Novel Method for Quantifying Smooth Regional Variations in Myocardial Contractility Within an Infarcted Human Left Ventricle Based on Delay-Enhanced Magnetic Resonance Imaging

Heart failure is increasing at an alarming rate, making it a worldwide epidemic. As the population ages and life expectancy increases, this trend is not likely to change. Myocardial infarction (MI)-induced adverse left ventricular (LV) remodeling is responsible for nearly 70% of heart failure cases. The adverse remodeling process involves an extension of the border zone (BZ) adjacent to an MI, which is normally perfused but shows myofiber contractile dysfunction. To improve patient-specific modeling of cardiac mechanics, we sought to create a finite element model of the human LV with BZ and MI morphologies integrated directly from delayed-enhancement magnetic resonance (DE-MR) images. Instead of separating the LV into discrete regions (e.g., the MI, BZ, and remote regions) with each having a homogeneous myocardial material property, we assumed a functional relation between the DE-MR image pixel intensity and myocardial stiffness and contractility--we considered a linear variation of material properties as a function of DE-MR image pixel intensity, which is known to improve the accuracy of the model's response. The finite element model was then calibrated using measurements obtained from the same patient--namely, 3D strain measurements-using complementary spatial modulation of magnetization magnetic resonance (CSPAMM-MR) images. This led to an average circumferential strain error of 8.9% across all American Heart Association (AHA) segments. We demonstrate the utility of our method for quantifying smooth regional variations in myocardial contractility using cardiac DE-MR and CSPAMM-MR images acquired from a 78-yr-old woman who experienced an MI approximately 1 yr prior. We found a remote myocardial diastolic stiffness of C(0) = 0.102 kPa, and a remote myocardial contractility of T(max) = 146.9 kPa, which are both in the range of previously published normal human values. Moreover, we found a normalized pixel intensity range of 30% for the BZ, which is consistent with the literature. Based on these regional myocardial material properties, we used our finite element model to compute patient-specific diastolic and systolic LV myofiber stress distributions, which cannot be measured directly. One of the main driving forces for adverse LV remodeling is assumed to be an abnormally high level of ventricular wall stress, and many existing and new treatments for heart failure fundamentally attempt to normalize LV wall stress. Thus, our noninvasive method for estimating smooth regional variations in myocardial contractility should be valuable for optimizing new surgical or medical strategies to limit the chronic evolution from infarction to heart failure.

Numerical evaluation of myofiber orientation and transmural contractile strength on left ventricular function

The left ventricle (LV) of the heart is composed of a complex organization of cardiac muscle fibers, which contract to generate force and pump blood into the body. It has been shown that both the orientation and contractile strength of these myofibers vary across the ventricular wall. The hypothesis of the current study is that the transmural distributions of myofiber orientation and contractile strength interdependently impact LV pump function. In order to quantify these interactions a finite element (FE) model of the LV was generated, which incorporated transmural variations. The influences of myofiber orientation and contractile strength on the Starling relationship and the end-systolic (ES) apex twist of the LV were assessed. The results suggest that reductions in contractile strength within a specific transmural layer amplified the effects of altered myofiber orientation in the same layer, causing greater changes in stroke volume (SV). Furthermore, when the epicardial myofibers contracted the strongest, the twist of the LV apex was greatest, regardless of myofiber orientation. These results demonstrate the important role of transmural distribution of myocardial contractile strength and its interplay with myofiber orientation. The coupling between these two physiologic parameters could play a critical role in the progression of heart failure.

Quasi-static image-based immersed boundary-finite element model of left ventricle under diastolic loading

Finite stress and strain analyses of the heart provide insight into the biomechanics of myocardial function and dysfunction. Herein, we describe progress toward dynamic patient-specific models of the left ventricle using an immersed boundary (IB) method with a finite element (FE) structural mechanics model. We use a structure-based hyperelastic strain-energy function to describe the passive mechanics of the ventricular myocardium, a realistic anatomical geometry reconstructed from clinical magnetic resonance images of a healthy human heart, and a rule-based fiber architecture. Numerical predictions of this IB/FE model are compared with results obtained by a commercial FE solver. We demonstrate that the IB/FE model yields results that are in good agreement with those of the conventional FE model under diastolic loading conditions, and the predictions of the LV model using either numerical method are shown to be consistent with previous computational and experimental data. These results are among the first to analyze the stress and strain predictions of IB models of ventricular mechanics, and they serve both to verify the IB/FE simulation framework and to validate the IB/FE model. Moreover, this work represents an important step toward using such models for fully dynamic fluid-structure interaction simulations of the heart. © 2014 The Authors. International Journal for Numerical Methods in Engineering published by John Wiley & Sons, Ltd.

Dynamic finite-strain modelling of the human left ventricle in health and disease using an immersed boundary-finite element method

Detailed models of the biomechanics of the heart are important both for developing improved interventions for patients with heart disease and also for patient risk stratification and treatment planning. For instance, stress distributions in the heart affect cardiac remodelling, but such distributions are not presently accessible in patients. Biomechanical models of the heart offer detailed three-dimensional deformation, stress and strain fields that can supplement conventional clinical data. In this work, we introduce dynamic computational models of the human left ventricle (LV) that are derived from clinical imaging data obtained from a healthy subject and from a patient with a myocardial infarction (MI). Both models incorporate a detailed invariant-based orthotropic description of the passive elasticity of the ventricular myocardium along with a detailed biophysical model of active tension generation in the ventricular muscle. These constitutive models are employed within a dynamic simulation framework that accounts for the inertia of the ventricular muscle and the blood that is based on an immersed boundary (IB) method with a finite element description of the structural mechanics. The geometry of the models is based on data obtained non-invasively by cardiac magnetic resonance (CMR). CMR imaging data are also used to estimate the parameters of the passive and active constitutive models, which are determined so that the simulated end-diastolic and end-systolic volumes agree with the corresponding volumes determined from the CMR imaging studies. Using these models, we simulate LV dynamics from enddiastole to end-systole. The results of our simulations are shown to be in good agreement with subject-specific CMR-derived strain measurements and also with earlier clinical studies on human LV strain distributions.

Bioinjection treatment: effects of post-injection residual stress on left ventricular wall stress

Injection of biomaterials into diseased myocardium has been associated with decreased myofiber stress, restored left ventricular (LV) geometry and improved LV function. However, its exact mechanism(s) of action remained unclear. In this work, we present the first patient-specific computational model of biomaterial injection that accounts for the possibility of residual strain and stress introduced by this treatment. We show that the presence of residual stress can create more heterogeneous regional myofiber stress and strain fields. Our simulation results show that the treatment generates low stress and stretch areas between injection sites, and high stress and stretch areas between the injections and both the endocardium and epicardium. Globally, these local changes are translated into an increase in average myofiber stress and its standard deviation (from 6.9 ± 4.6 to 11.2 ± 48.8 kPa and 30 ± 15 to 35.1 ± 50.9 kPa at end-diastole and end-systole, respectively). We also show that the myofiber stress field is sensitive to the void-to-size ratio. For a constant void size, the myofiber stress field became less heterogeneous with decreasing injection volume. These results suggest that the residual stress and strain possibly generated by biomaterial injection treatment can have large effects on the regional myocardial stress and strain fields, which may be important in the remodeling process.

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.

Patient-specific finite element modeling of the Cardiokinetix Parachute(®) device: effects on left ventricular wall stress and function

The Parachute(®) (Cardiokinetix, Inc., Menlo Park, California) is a catheter-based device intended to reverse left ventricular (LV) remodeling after antero-apical myocardial infarction. When deployed, the device partitions the LV into upper and lower chambers. To simulate its mechanical effects, we created a finite element LV model based on computed tomography (CT) images from a patient before and 6 months after Parachute(®) implantation. Acute mechanical effects were determined by in silico device implantation (VIRTUAL-Parachute). Chronic effects of the device were determined by adjusting the diastolic and systolic material parameters to better match the 6-month post-implantation CT data and LV pressure data at end-diastole (ED) (POST-OP). Regional myofiber stress and pump function were calculated in each case. The principal finding is that VIRTUAL-Parachute was associated with a 61.2 % reduction in the lower chamber myofiber stress at ED. The POST-OP model was associated with a decrease in LV diastolic stiffness and a larger reduction in myofiber stress at the upper (27.1%) and lower chamber (78.4%) at ED. Myofiber stress at end-systole and stroke volume was little changed in the POST-OP case. These results suggest that the primary mechanism of Parachute(®) is a reduction in ED myofiber stress, which may reverse eccentric post-infarct LV hypertrophy.

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.

Reduction in left ventricular wall stress and improvement in function in failing hearts using Algisyl-LVR

Injection of Algisyl-LVR, a treatment under clinical development, is intended to treat patients with dilated cardiomyopathy. This treatment was recently used for the first time in patients who had symptomatic heart failure. In all patients, cardiac function of the left ventricle (LV) improved significantly, as manifested by consistent reduction of the LV volume and wall stress. Here we describe this novel treatment procedure and the methods used to quantify its effects on LV wall stress and function. Algisyl-LVR is a biopolymer gel consisting of Na(+)-Alginate and Ca(2+)-Alginate. The treatment procedure was carried out by mixing these two components and then combining them into one syringe for intramyocardial injections. This mixture was injected at 10 to 19 locations mid-way between the base and apex of the LV free wall in patients. Magnetic resonance imaging (MRI), together with mathematical modeling, was used to quantify the effects of this treatment in patients before treatment and at various time points during recovery. The epicardial and endocardial surfaces were first digitized from the MR images to reconstruct the LV geometry at end-systole and at end-diastole. Left ventricular cavity volumes were then measured from these reconstructed surfaces. Mathematical models of the LV were created from these MRI-reconstructed surfaces to calculate regional myofiber stress. Each LV model was constructed so that 1) it deforms according to a previously validated stress-strain relationship of the myocardium, and 2) the predicted LV cavity volume from these models matches the corresponding MRI-measured volume at end-diastole and end-systole. Diastolic filling was simulated by loading the LV endocardial surface with a prescribed end-diastolic pressure. Systolic contraction was simulated by concurrently loading the endocardial surface with a prescribed end-systolic pressure and adding active contraction in the myofiber direction. Regional myofiber stress at end-diastole and end-systole was computed from the deformed LV based on the stress-strain relationship.