|
|
||||||||
a Department of Mechanical Engineering, University of Kentucky, Lexington, Kentucky
b Department of Surgery, University of Kentucky, Lexington, Kentucky
c Department of Surgery, University of California, San Francisco, California
d Department of Bioengineering, University of California, San Francisco, California
e Department of Radiology and Biomedical Imaging, University of California, San Francisco, California
f Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
g Gorman Cardiovascular Research Group, University of Pennsylvania, Philadelphia, Pennsylvania
Accepted for publication December 19, 2011.
* Address correspondence to Dr Guccione, UCSF/VA Medical Center (112D), 4150 Clement St, San Francisco, CA 94121 (Email: guccionej{at}surgery.ucsf.edu).
| Abstract |
|---|
|
|
|---|
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 patient's 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.
| Introduction |
|---|
|
|
|---|
An important feature of this investigation is the combined use of three-dimensional complementary spatial modulation of magnetization for measuring regional myocardial deformation and a realistic mathematical (FE) model of the infarcted LV for computing regional myocardial force development. The technique has been previously validated by direct ex vivo active force measurements in skinned fiber preparations [5]. Direct measurement of such forces in an intact LV has been unreliable in animals and is not possible in patients.
| Material and Methods |
|---|
|
|
|---|
Experimental Measurements
A 62-year-old man experienced an MI in 1985. In 2010, he underwent coronary artery bypass grafting. Before coronary artery bypass grafting, persantine thallium stress testing demonstrated reversible defects that were all amenable to surgical revascularization. The patient underwent three-dimensional complementary spatial modulation of magnetization (Fig 1
A) 6 weeks after coronary artery bypass grafting, at which time there was no evidence of ongoing ischemia. In addition, the patient underwent magnetic resonance delayed hyperenhancement to precisely delineate the infarct region. Left ventricular pressure was continuously measured by means of cardiac catheterization. The endocardial and epicardial LV surfaces were contoured from the magnetic resonance images, and the systolic tags were segmented to compute the systolic myocardial strain at the midwall of the LV [5, 6].
|
Boundary and Loading Conditions
The boundary conditions of the LV were assigned to fully constrain the displacement at the epicardial-basal edge, while allowing the remaining nodes at the base to move in the circumferential-radial plane. The inner endocardial wall of the LV was loaded with the clinically measured pressures to simulate the end-diastolic (ED) and end-systolic (ES) states. It should be noted that the pressure boundary conditions were applied by quickly ramping the load to the measured value and then holding the load constant as the simulation reached a steady-state solution at ED and ES.
Material Properties of the Heart
Nearly incompressible, transversely isotropic, hyperelastic constitutive laws for passive [8] and active myocardium [9] were modeled with a user-defined material subroutine in the explicit FE solver, LS-DYNA (Livermore Software Technology Corporation, Livermore, CA). Passive material properties were represented by the strain energy function:
|
| (1) |
where E 11 is fiber strain, E 22 is cross-fiber strain, E 33 is radial strain, E 23 is shear strain in the transverse plane, and E 12 and E 13 are shear strain in the fiber-cross fiber and fiber-radial planes, respectively. Values for the material constants b f , b t , and b fs were chosen as 24.63, 9.63, and 8.92, respectively, which were scaled (ie, a value of 6 was added to each exponential material constant) from previous studies of canine myocardium [8]. To validate the models, the material constant C was adjusted until the LV ED volumes matched the experimentally measured values. As was the case in our previous FE models of the infarcted ovine LV [2, 5, 6], we assume that C has the same value in the BZ and remote myocardium but is 10 times stiffer in the infarct.
Active contraction was modeled by defining total stress as the sum of the passive stress derived from the strain energy function and an active fiber directional component, T
0, which is a function of time, t, peak intracellular calcium concentration, Ca
0, sarcomere length, l, and maximal isometric tension achieved at the longest sarcomere length, T
max [5]:
|
| (2) |
where S is the second Piola-Kirchoff stress tensor, p is the hydrostatic pressure introduced as the Lagrange multiplier needed to ensure incompressibility, J is the Jacobian of the deformation gradient tensor, C is the right Cauchy-Green deformation tensor, Dev is the deviatoric projection operator:
|
| (3) |
and
is the deviatoric contribution of the strain energy function,
(Eq 1). The assumption of near incompressibility of the myocardium requires the decoupling of the strain energy function into dilational and deviatoric components:
|
| (4) |
where U is the dilatational (volumetric) contribution [10].
The active fiber-directed stress component is defined by a time-varying elastance model, which at end-systole, is reduced to [11]:
|
| (5) |
with m and b as constants, and the length-dependent calcium sensitivity, ECa
50, is given by:
|
| (6) |
where B is a constant, (Ca 0)max is the maximal peak intracellular calcium concentration, l 0 is the sarcomere length at which no active tension develops, and l R is the stress-free sarcomere length. The material constants for active contraction were taken to be [12] Ca 0, 4.35 μmol/L; (Ca 0)max, 4.35 μmol/L; B, 4.75 μm–1; l 0, 1.58 μm; m, 1.0489 s/μm; b, –1.429 seconds, and l R was set at 1.85 μm, the sarcomere length in the unloaded reference configuration. On the basis of biaxial stretching experiments [13] and FE analyses [14, 15], cross-fiber, in-plane stress equivalent to 40% of that along the myocardial fiber direction was added. To further validate the models, the parameter T max was formally optimized.
Material Parameter Optimization
The objective function for the optimization was taken to be the mean squared error (MSE) [1]. The passive material parameters were determined such that the FE model-predicted LV ED volume matched the patient-specific in vivo measured value. The systolic material parameters in the remote (T
max_R) and BZ (T
max_B) regions were estimated by minimizing the errors between FE model-predicted and in vivo measured systolic strains and LV ES volume. The goal of the optimization is to minimize the MSE:
|
| (7) |
where n is the in vivo strain point, N is the total number of in vivo strain points,
is the computed FE strain at each strain point, and V
ED and V
ES are the computed FE LV ED and ES volumes, respectively. The overbar represents the experimental in vivo measurements.
| Results |
|---|
|
|
|---|
The optimization displayed good convergence, with a 90% confidence interval of 18%, which only took 10 iterations to reach (Fig 2). However, the optimization procedure was allowed to run 5 extra iterations to ensure convergence. The final MSE value of 8.53 was calculated using 816 strain components and two volumes. The MSE value indicates generally good agreement between the FE model-predicted systolic strains and the patient-specific in vivo measured strains, and is similar in magnitude to our previous studies, which were validated with infarcted ovine hearts [5, 6].
|
|
|
| Comment |
|---|
|
|
|---|
Myocardial tissue tagging using complementary spatial modulation of magnetization allows detailed assessment of myocardial motion. To capture the complex three-dimensional cardiac motion pattern, multiple two-dimensional tagged slices are usually acquired in different orientations. These approaches are prone to slice misregistration and associated with long acquisition times. In this work, we applied the accelerated three-dimensional tagging acquisition method of Rutz and coworkers [17], which enabled assessment of three-dimensional motion information with whole heart coverage in three short breathholds. They found hypokinetic regions in patients with an MI corresponded well with regions exhibiting hyperenhancement after contrast injection. However, Rutz and colleagues did not attempt to quantify the forces or stresses responsible for that hypokinesis.
Using an MRI-based FE stress analysis of a clinically relevant large animal preparation, Guccione and associates [2] suggested that the mechanism underlying mechanical dysfunction in the BZ of LV aneurysm is primarily the result of myocardial contractile dysfunction rather than increased wall stress in this region. Then, Walker and coworkers [15, 18] used a fixed ratio of 2:1 in remote versus BZ region contractility to compute regional myocardial material properties and stress in 6 animals [18] after linear repair of LV aneurysm and in 5 animals that underwent a sham operation [15]. More recently, Sun and colleagues [5] developed a computationally efficient formal optimization that allows regional myocardial contractility to be quantified without enforcing the fixed ratio mentioned above. Moreover, in vivo estimates of regional myocardial contractility were validated using ex vivo direct force measurements in skinned fiber preparations. Then, that formal optimization was applied in 6 animals 2 weeks before and 2 and 6 weeks after patch repair of LV aneurysm [16], as well as in an animal with posterobasal MI [6]. Most important, in every single animal and time included in those previous studies (32 instances), myocardial contractility in the BZ of the MI was significantly less than that in regions remote from the MI.
In this paper we present two distinct values for the contractility parameters for BZ (T max_B) and remote myocardium (T max_R). In a recent study of an ovine model of LV aneurysm [19], we found that our corresponding FE model was better able to reproduce the experimental LV pressure versus myocardial strain data when we allowed T max_B to have values that vary linearly from 0 at the boundary between the aneurysm and BZ to T max_R at the boundary between the BZ and remote myocardium. We expect a similar gradient in T max_B to exist in the BZ of other types of MIs, including the MI of the patient in the present study.
At present, the greatest challenge in quantifying regional myocardial contractility in patients is the very sophisticated MRI hardware required to quantify regional myofiber orientation in vivo. The diffusion imaging approach proposed by Gamper and associates [20] appears to be the state of the art, but it requires an MRI scanner with uniquely large gradients (of at least 80 mT/m). Because we did not have access to such a scanner, we instead repeated our quantification of regional myocardial contractility in a patient with an MI using a range of myofiber angle distributions. Case 2 of the fiber angle sensitivity study corresponds to the in vivo diffusion MRI data of Wu and coworkers [21, 22] concerning patients with an MI in which the percentage of left-handed helical fibers (negative subepicardial helix angles) increased from the remote zone to the BZ and infarct zone, and the percentage of right-handed helical fibers (positive subendocardial helix angles) decreased from the remote zone to the BZ and infarct zone. The findings of Wu and colleagues [21, 22] are at variance with the ex vivo observations of Chen and associates [23]. Recent ex vivo diffusion MRI data from infarcted rat hearts indicate that BZ helix angles are preserved [24, 25].
This investigation demonstrates the substantial depression of myocardial contractility in the BZ of a human MI relative to that in regions remote from the MI. The current study evaluates one patient at one time very remote from the index MI; however, the results of this study support the central role of the BZ myocardium in the post-MI remodeling process in patients that has been demonstrated in large animal infarct models. To fully establish this important conclusion, serial evaluations beginning at the time of the index MI will need to be performed in a large cohort of patients. The current study supports the importance and demonstrates the feasibility of such large and long-term studies.
| Acknowledgments |
|---|
|
|
|---|
| References |
|---|
|
|
|---|
Related Article
This article has been cited by other articles:
![]() |
J. Braun Invited Commentary Ann. Thorac. Surg., April 1, 2012; 93(4): 1193 - 1194. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ANN THORAC SURG | ASIAN CARDIOVASC THORAC ANN | EUR J CARDIOTHORAC SURG |
| J THORAC CARDIOVASC SURG | ICVTS | ALL CTSNet JOURNALS |