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

Post-infarct evolution of ventricular and myocardial function

Adverse ventricular remodeling following acute myocardial infarction (MI) may induce ventricular dilation, fibrosis, and loss of global contractile function, possibly resulting in heart failure (HF). Understanding the relation between the time-dependent changes in material properties of the myocardium and the contractile function of the heart may further our understanding of the development of HF post-MI and guide the development of novel therapies. A finite element model of cardiac mechanics was used to model MI in a thick-walled truncated ellipsoidal geometry. Infarct core and border zone comprised 9.6 and 8.1% of the LV wall volume, respectively. Acute MI was modeled by inhibiting active stress generation. Chronic MI was modeled by the additional effect of infarct material stiffening, wall thinning and fiber reorientation. In acute MI, stroke work decreased by 25%. In the infarct core, fiber stress was reduced but fiber strain was increased, depending on the degree of infarct stiffening. Fiber work density was equal to zero. Healthy tissue adjacent to the infarct showed decreased work density depending on the degree of infarct stiffness and the orientation of the myofibers with respect to the infarct region. Thinning of the wall partially restored this loss in work density while the effects of fiber reorientation were minimal. We found that the relative loss in pump function in the infarcted heart exceeds the relative loss in healthy myocardial tissue due to impaired mechanical function in healthy tissue adjacent to the infarct. Infarct stiffening, wall thinning and fiber reorientation did not affect pump function but did affect the distribution of work density in tissue adjacent to the infarct.

A Lumped Two-Compartment Model for Simulation of Ventricular Pump and Tissue Mechanics in Ischemic Heart Disease

Introduction: Computational modeling of cardiac mechanics and hemodynamics in ischemic heart disease (IHD) is important for a better understanding of the complex relations between ischemia-induced heterogeneity of myocardial tissue properties, regional tissue mechanics, and hemodynamic pump function. We validated and applied a lumped two-compartment modeling approach for IHD integrated into the CircAdapt model of the human heart and circulation.

Methods: Ischemic contractile dysfunction was simulated by subdividing a left ventricular (LV) wall segment into a hypothetical contractile and noncontractile compartment, and dysfunction severity was determined by the noncontractile volume fraction (NCVF). Myocardial stiffness was determined by the zero-passive stress length (Ls0,pas) and nonlinearity (kECM) of the passive stress-sarcomere length relation of the noncontractile compartment. Simulated end-systolic pressure volume relations (ESPVRs) for 20% acute ischemia were qualitatively compared between a two- and one-compartment simulation, and parameters of the two-compartment model were tuned to previously published canine data of regional myocardial deformation during acute and prolonged ischemia and reperfusion. In six patients with myocardial infarction (MI), the NCVF was automatically estimated using the echocardiographic LV strain and volume measurements obtained acutely and 6 months after MI. Estimated segmental NCVF values at the baseline and 6-month follow-up were compared with percentage late gadolinium enhancement (LGE) at 6-month follow-up.

Results: Simulation of 20% of NCVF shifted the ESPVR rightward while moderately reducing the slope, while a one-compartment simulation caused a leftward shift with severe reduction in the slope. Through tuning of the NCVF, Ls0,pas, and kECM, it was found that manipulation of the NCVF alone reproduced the deformation during acute ischemia and reperfusion, while additional manipulations of Ls0,pas and kECM were required to reproduce deformation during prolonged ischemia and reperfusion. Out of all segments with LGE>25% at the follow-up, the majority (68%) had higher estimated NCVF at the baseline than at the follow-up. Furthermore, the baseline NCVF correlated better with percentage LGE than NCVF did at the follow-up.

Conclusion: We successfully used a two-compartment model for simulation of the ventricular pump and tissue mechanics in IHD. Patient-specific optimizations using regional myocardial deformation estimated the NCVF in a small cohort of MI patients in the acute and chronic phase after MI, while estimated NCVF values closely approximated the extent of the myocardial scar at the follow-up. In future studies, this approach can facilitate deformation imaging–based estimation of myocardial tissue properties in patients with cardiovascular diseases.

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 study of border zone formation in ischemic heart using electro-chemical coupled computational model

Myocardial infarction (MI) is the most common cause of a heart failure, which occurs due to myocardial ischemia leading to left ventricular (LV) remodeling. LV remodeling particularly occurs at the ischemic area and the region surrounds it, known as the border zone. The role of the border zone in initiating LV remodeling process urges the investigation on the correlation between early border zone changes and remodeling outcome. Thus, this study aims to simulate a preliminary conceptual work of the border zone formation and evolution during onset of MI and its effect towards early LV remodeling processes by incorporating the oxygen concentration effect on the electrophysiology of an idealized three-dimensional LV through electro-chemical coupled mathematical model. The simulation result shows that the region of border zone, represented by the distribution of electrical conductivities, keeps expanding over time. Based on this result, the border zone is also proposed to consist of three sub-regions, namely mildly, moderately, and seriously impaired conductivity regions, which each region categorized depending on its electrical conductivities. This division could be used as a biomarker for classification of reversible and irreversible myocardial injury and will help to identify the different risks for the survival of patient. Larger ischemic size and complete occlusion of the coronary artery can be associated with an increased risk of developing irreversible injury, in particular if the reperfusion treatment is delayed. Increased irreversible injury area can be related with cardiovascular events and will further deteriorate the LV function over time.

The study of myocardial ischemia-reperfusion treatment through computational modelling

Reperfusion of the blood flow to ischemic myocardium is the standard treatment for patients suffering myocardial infarction. However, the reperfusion itself can also induce myocardial injury, in which the actual mechanism and its risk factors remain unclear. This work aims to study the mechanism of ischemia-reperfusion treatment using a three-dimensional (3D) oxygen diffusion model. An electrical model is then coupled to an oxygen model to identify the possible region of myocardial damage. Our findings show that the value of oxygen exceeds its optimum (>1.0) at the ischemic area during early reperfusion period. This complication was exacerbated in a longer ischemic period. While a longer reperfusion time causes a continuous excessive oxygen supply to the ischemic area throughout the reperfusion time. This work also suggests the use of less than 0.8 of initial oxygen concentration in the reperfusion treatment to prevent undesired upsurge at the early reperfusion period and further myocardial injury. We also found the region at risk for myocardial injury is confined in the ischemic vicinity revealed by its electrical conductivity impairment. Although there is a risk that reperfusion leads to myocardial injury for excessive oxygen accumulation, the reperfusion treatment is helpful in reducing the infarct size.

The role of end-diastolic myocardial fibre stretch on infarct extension

Myocardial infarct extension, a process involving the enlargement of infarct and border zone, leads to progressive degeneration of left ventricular (LV) function and eventually gives rise to heart failure. Despite carrying a high risk, the causation of infarct extension is still a subject of much speculation. In this study, patient-specific LV models were developed to investigate the correlation between infarct extension and impaired regional mechanics. Subsequently, sensitivity analysis was performed to examine the causal factors responsible for the impaired regional mechanics observed in regions surrounding the infarct and border zone. From our simulations, fibre strain, fibre stress and fibre stress-strain loop (FSSL) were the key biomechanical variables affected in these regions. Among these variables, only FSSL was correlated with infarct extension, as reflected in its work density dissipation (WDD) index value, with high WDD indices recorded at regions with infarct extension. Impaired FSSL is caused by inadequate contraction force generation during the isovolumic contraction and ejection phases. Our further analysis revealed that the inadequacy in contraction force generation is not necessarily due to impaired myocardial intrinsic contractility, but at least in part, due to inadequate muscle fibre stretch at end-diastole, which depresses the ability of myocardium to generate adequate contraction force in the subsequent systole (according to the Frank-Starling law). Moreover, an excessively stiff infarct may cause its neighbouring myocardium to be understretched at end-diastole, subsequently depressing the systolic contractile force of the neighbouring myocardium, which was found to be correlated with infarct extension.

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.

Impact of myocardial infarction on intraventricular vortex and flow energetics assessed using computational simulations

Flow energetics have been proposed as early indicators of progressive left ventricular (LV) functional impairment in patients with myocardial infarction (MI), but its correlation with individual MI parameters has not been fully explored. Using electro-fluid-structure interaction LV models, this study investigated the correlation between four MI parameters: infarct size, infarct multiplicity, regional enhancement of contractility at the viable myocardium area (RECVM), and LV mechanical dyssynchrony (LVMD) with intraventricular vortex and flow energetics. In LV with small infarcts, our results showed that infarct appearance amplified the energy dissipation index (DI), where substantial viscous energy loss was observed in areas with high flow velocity and near the infarct-vortex interface. The LV with small multiple infarcts and RECVM showed remarkable DI increment during systole and diastole. In correlation analysis, the systolic kinetic energy fluctuation index (E') was positively related to ejection fraction (EF) (R²  = 0.982) but negatively correlated with diastolic E' (R²  = 0.970). Diastolic E' was inversely correlated with vortex kinetic energy (R²  = 0.960) and vortex depth (R²  = 0.876). We showed an excessive systolic DI could differentiate infarcted LV with normal EF from healthy LV. Strong flow acceleration, LVMD, and vortex-infarct interactions were predominant factors that induced excessive DI in infarcted LVs. Instead of causing undesired flow turbulence, high systolic E' suggested the existence of energetic flow acceleration, while high diastolic E' implied an inefficient diastolic filling. Thus, systolic E' is not a suitable early indicator for progressive LV dysfunction in MI patients, while diastolic E' may be a useful index to indicate diastolic impairment in these patients.

Post-infarct morphine treatment mitigates left ventricular remodeling and dysfunction in a rat model of ischemia-reperfusion

Following myocardial infarction, the heart undergoes a series of dramatic compensations which may later form a maladaptive picture characterized by ventricular dilation and pump failure. Among several opioid agents, morphine has been shown to confer protection against reperfusion injury and infarct size. Here, we sought to study the cardioprotective effect of post-infarct morphine treatment against left ventricular adverse remodeling. We induced myocardial infarction in male Sprague - Dawley rats by ligating left anterior descending artery and then, treated these animals with three different doses of morphine -0.3, 3 and 10 mg/kg (i.p.). The echocardiographic evaluation depicted improved cardiac performance and lesser chamber dilation in the animals that had received 3 mg/kg of morphine. Next, we studied the effect of 3 mg/kg morphine administration on left ventricular hemodynamics, infarct size, tissue architecture, changes in lung and heart weight, circulating TNF-α level and post-MI mRNA expression of collagen-1, collagen-3, TGF-β, TNF-α, MMP-2 and MMP-9. Five-day morphine administration markedly improved LV function, and also reduced infarct size, myocyte hypertrophy, fibrosis, index of infarct expansion, heart weight and serum TNF-α level. Moreover, morphine alleviated MI-induced increase in wet and dry lung weight. Morphine also altered the mRNA expression of fibrosis-related genes, TNF-α, MMP-2 and MMP-9. In conclusion, post-infarct morphine treatment can mitigate adverse remodeling and cardiac dysfunction after MI. Beside analgesic effect, we may be able to harvest benefits from the antifibrotic and anti-remodeling action of morphine in patients with the acute coronary syndrome.

Gaussian Process Regressions for Inverse Problems and Parameter Searches in Models of Ventricular Mechanics

Patient specific models of ventricular mechanics require the optimization of their many parameters under the uncertainties associated with imaging of cardiac function. We present a strategy to reduce the complexity of parametric searches for 3-D FE models of left ventricular contraction. The study employs automatic image segmentation and analysis of an image database to gain geometric features for several classes of patients. Statistical distributions of geometric parameters are then used to design parametric studies investigating the effects of: (1) passive material properties during ventricular filling, and (2) infarct geometry on ventricular contraction in patients after a heart attack. Gaussian Process regression is used in both cases to build statistical models trained on the results of biophysical FEM simulations. The first statistical model estimates unloaded configurations based on either the intraventricular pressure or the end-diastolic fiber strain. The technique provides an alternative to the standard fixed-point iteration algorithm, which is more computationally expensive when used to unload more than 10 ventricles. The second statistical model captures the effects of varying infarct geometries on cardiac output. For training, we designed high resolution models of non-transmural infarcts including refinements of the border zone around the lesion. This study is a first effort in developing a platform combining HPC models and machine learning to investigate cardiac function in heart failure patients with the goal of assisting clinical diagnostics.