|
|
||||||||
J Thorac Cardiovasc Surg 2008;135:1094-1102
© 2008 The American Association for Thoracic Surgery
Surgery for Acquired Cardiovascular Disease |
a Joint Graduate Group in Bioengineering, University of California, San Francisco and Berkeley, Calif
b Department of Surgery, University of California, San Francisco, Calif
c Department of Anesthesia, University of California, San Francisco, Calif
d Department of Radiology, University of California, San Francisco, Calif
e Department of Veterans Affairs Medical Center, San Francisco, Calif
f Department of Biomedical Engineering, Duke University, Durham, NC
Received for publication May 10, 2007; revisions received October 27, 2007; accepted for publication November 15, 2007. * Address for reprints: Julius M. Guccione, PhD, Division of Surgical Services (112D), VA Medical Center, 4150 Clement Street, San Francisco, CA 94121. (Email: GuccioneJ{at}surgery.ucsf.edu).
| Abstract |
|---|
|
|
|---|
Methods: Nine sheep underwent myocardial infarction and aneurysm repair with a linear repair (13.4 ± 2.3 weeks postmyocardial infarction). Satisfactory magnetic resonance imaging examinations were obtained in 6 sheep (6.6 ± 0.5 weeks postrepair). Finite element models were constructed from in vivo magnetic resonance imaging-based cardiac geometry and postmortem measurement of myofiber helix angles using diffusion tensor magnetic resonance imaging. Material properties were iteratively determined by comparing the finite element model output with systolic tagged magnetic resonance imaging strain measurements.
Results: At the mid-wall, fiber stress in the border zone decreased by 39% (sham = 32.5 ± 2.5 kPa, repair = 19.7 ± 3.6 kPa, P = .001) to the level of remote regions after repair. In the septum, however, border zone fiber stress remained high (sham = 31.3 ± 5.4 kPa, repair = 23.8 ± 5.8 kPa, P = .29). Cross-fiber stress at the mid-wall decreased by 41% (sham = 13.0 ± 1.5 kPa, repair = 7.7 ± 2.1 kPa, P = .01), but cross-fiber stress in the un-excluded septal infarct was 75% higher in the border zone than remote regions (remote = 5.9 ± 1.9 kPa, border zone = 10.3 ± 3.6 kPa, P < .01). However, end-diastolic fiber and cross-fiber stress were not reduced in the remote myocardium after plication.
Conclusion: With the exception of the retained septal infarct, end-systolic stress is reduced in all areas of the left ventricle after infarct plication. Consequently, we expect the primary positive effect of infarct plication to be in the infarct border zone. However, the amount of stress reduction necessary to halt or reverse nonischemic infarct extension in the infarct border zone and eccentric hypertrophy in the remote myocardium is unknown.
| Introduction |
|---|
|
|
|---|
A previous study from our laboratory analyzed the end-diastolic and end-systolic pressure–volume relationships in a cohort of sheep treated with LV aneurysm plication.5
The results showed acute decreases in LV volume and diastolic compliance, an increase in end-systolic elastance, and an acute decrease in the Starling relationship. After 6 weeks, however, these parameters returned to preplication values. This redilation was hypothesized to be the result of high wall stress in the repair border zone and remote myocardium. More recently, we analyzed end-systolic 3-dimensional strain in another cohort of sheep treated with LV aneurysm plication and found that repair of LV aneurysm significantly increases systolic circumferential shortening at the border zone in sheep.6
Because there is no available method to directly measure LV wall stress,7
it must be estimated computationally or from simplified analytic expressions. At present, there are various equations available to compute approximate LV stress;8
however, the most versatile estimates use the finite element (FE) method.9-13
By using FE analysis, a recent study from our laboratory calculated mid-wall stresses in 4 sheep LVs with anteroapical aneurysms.14
Material properties in FE models were iteratively determined by comparing calculated with experimentally measured strains. FE models also included the most recent estimates of aneurysm material properties15
and myofiber helix angles.16
Stress measurements in the infarct border zone showed a 24% increase in fiber stress and a 115% increase in cross-fiber stress relative to noninfarct tissue. Moreover, it was found that a transverse contracting stress equal to 40% of that along the muscle fiber improved agreement with experimental measurements by 24%.
By using the same methodology, the present study used FE analysis to compute mid-wall stress in sheep hearts that have undergone aneurysm plication. Material properties were determined that allowed the FE models to reproduce end-systolic 3-dimensional strain measurements made with tagged magnetic resonance imaging (MRI). FE models included myocyte fiber angle measurements made for the first time after this procedure.16
The present study was undertaken to test the hypothesis that end-systolic stress is reduced in all areas of the LV after infarct plication.
| Materials and Methods |
|---|
|
|
|---|
Myocardial Infarction
The infarct procedure has been described in detail.17
In brief, 8 adult sheep underwent anteroapical MI. Male Dorsett sheep (40–50 kg) were anesthetized (ketamine, 33 mg/kg (intramuscularly), maintained with isoflurane (2%–4% inspired), and mechanically ventilated (tidal volume 20 mL/kg; model 309-0612-800, Ohio Medical Products, Madison, Wis). A left thoracotomy was performed using a sterile technique, and the left anterior descending and second left anterior descending diagonal coronary arteries were ligated at a point 40% of the distance from the apex to the base as previously described.17
If present, branches of the posterior descending artery that also perfused this region were ligated 20% of the distance from the apex to the base. The thoracotomy was closed in layers, and the sheep recovered from anesthesia.
Aneurysm Plication
After development of a ventricular aneurysm, aneurysm plication was performed without cardiopulmonary bypass (13.4 ± 2.3 weeks post-MI). A partial lower sternotomy was performed under general anesthesia, and pericardial adhesions were divided. The transition between infarcted and noninfarcted myocardium was palpated, and the LV aneurysm was plicated between 2 strips of Dacron felt. Polypropylene sutures (2-0 Prolene, MH needle, Ethicon Inc, Sommerville, NJ) were passed through the felt, through and through the aneurysm at its border, and through the opposite strip of felt in a horizontal mattress fashion. The sternotomy was closed, and the sheep recovered from anesthesia.
Experimental Measurements
Six weeks after the plication procedure (20.0 ± 2.9 weeks post-MI, 6.6 ± 0.5 weeks postplication), systolic strains were measured with tagged MRI as described in greater detail previously.6
In summary, the sheep was intubated and anesthetized (1% inspired isoflurane [Forane]). Under fluoroscopic guidance, a nonferromagnetic transducer-tipped pressure catheter (model SPC-320, Millar Instruments, Houston, Tex) was introduced into the LV via sterile neck incisions. The animal was transported to the MRI scanner and positioned on its left side in the magnet (Siemens Symphony MRI scanner, 1.5 Tesla with Quantum gradients; Siemens Medical Systems, Iselin, NJ) with its chest centered in a Helmholtz coil. A series of tagged magnetic resonance images were acquired in orthogonal short- and long-axis planes synchronized to the R wave of the electrocardiogram signal at 40-ms intervals for the approximate completion of the cardiac cycle. LV pressure was recorded intermittently between slice acquisitions.
Image acquisition was synchronized to the R wave of the electrocardiogram signal. During the actual image data acquisition, the ventilator (Hallowell 2000, Hallowell EMC, Pittsfield, Mass) was stopped for approximately 10 seconds at maximum inspiration (to minimize respiratory motion and the associated motion artefacts in our magnetic resonance images). During this time period, a series of images were acquired at 50-ms intervals until the approximate completion of the entire cardiac cycle. A segmented cine gradient echo sequence was used, the cardiac cycle was divided into 50-ms windows, and 7 phase encoding lines were acquired in each window. Image data were collected during 16 heart cycles providing 99 phase encoding steps. The number of cardiac phases collected depended on the specific animal's R-R interval, with 10 cardiac phases being typical. The gradient echo sequence had an echo time of 4 ms, an excitation angle (
) of 20 degrees, and an acquisition matrix of 256 x 102. The field of view was set to 280 x 230 mm2 and 280 x 230 mm2 for the short- and long-axis images, respectively. Total MRI study time was approximately 45 minutes. Raw data were interpolated to 256 x 256 pixel images, and these were transferred to a Silicon Graphics workstation (Silicon Graphics Inc, Mountain View, Calif) and converted from Siemens format to SGI format using custom software. Systolic myocardial strains were calculated from tag-line deformation18
in cardiac coordinates (ie, circumferential, longitudinal, and radial) at the mid-wall in 12 equally spaced "sectors" around the circumference in each short-axis slice.
Fiber Angles
Helix angle measurements made using diffusion tensor magnetic resonance imaging (DTMRI) in the same hearts16
were incorporated into the FE model by aligning geometry from the 2 data sets and sampling the nearest points from the DTMRI data set, as described previously.14
For the data set alignment, we aligned the valve plane, LV long-axis, and right ventricular insertion points. All alignments were rigid body motions. The long axes were determined from a least-squares fit to the centroid of the epicardial short-axis contours. Fiber angles in the septal aneurysm were set to 0 degrees. The aneurysmal fiber angle is assumed to be zero because the anisotropy of the aneurysm is referred to "cardiac" coordinates. Fiber angles of zero and 90 degrees correspond to the circumferential and longitudinal directions, respectively. During biaxial mechanical testing of sheep LV aneurysms, the specimens are stretched in the circumferential and longitudinal directions. The results from these tests indicate that sheep LV aneurysms are anisotropic (ie, significantly stiffer in the longitudinal direction).
Finite Element Analysis
The FE analysis followed the same methodology as described previously14
and summarized here. FE models replicated in vivo geometry at end-isovolumic relaxation and contained 18 elements longitudinally, 12 elements circumferentially, and 1 element transmurally with linear interpolation in prolate spheroidal coordinates. Remote, border zone, and aneurysmal (in the septum) regions were determined from ventricular wall thickness. The border zone region was defined at the steep transition in wall thickness between the remote and aneurysm regions. As in previous studies,13,14
azimuthal displacements at all nodes at the base and apex and circumferential displacement of the basal epicardial nodes were constrained. In prolate spheroidal coordinates, the azimuthal direction is aligned with the longitudinal or LV apex-to-base direction only at the equator. At the apex, the azimuthal direction is aligned with the circumferential direction. Thus, these azimuthal constraints do not prevent apex-to-base shortening.
Material Law
The material laws have been described.13,14
In short, the passive myocardium was modeled by a strain energy function, W, that is anisotropic relative to the local fiber direction:
|
| (1) |
|
| (2) |
Active contraction was simulated by adding stress in the muscle fiber direction defined by a time-varying elastance model.19
For end systole, this reduced to the following:
|
| (3) |
|
| (4) |
|
| (5) |
Determination of Material Parameters
In the septal aneurysm, material parameters determined from biaxial stretching experiments were used.15
In the border zone, Tmax from equation 3 was reduced by 50% in accordance with results from Guccione and colleagues13
and Walker and colleagues.14
Because of a lack of experimental data of aneurysm tissue under shear loading, bfs was kept constant across all regions. In accordance with biaxial stretching experiments of Lin and Yin22
and results from our previous study,14
an in-plane, cross-fiber stress was added equal to 40% of that along the muscle fiber direction.
By following an iterative approach described previously14
and illustrated in
Figure 1, transversely isotropic material parameters were systematically varied and Tmax scaled to reproduce measured ventricular volumes and myocardial strains. The FE model was first inflated to the measured offset end-diastolic pressure (
Table 1) with initial material parameters from previously published studies;23-25
material parameter C from equation 1 was scaled in nonaneurysmal regions to obtain the correct end-diastolic volume. Once the end-diastolic volume was within a few percent, the model was inflated to the offset end-systolic pressure and systolic contraction simulated in the remote region and Tmax was scaled accordingly to obtain the end-systolic volume within a few percent.14
|
|
Statistical Analysis
Stress values were taken from the central integration point of each element, weighted by the element volume, and averaged over regions. Statistical analysis was performed using paired t tests within groups and unpaired t tests between groups. Stress calculations addressed 4 fundamental questions (1 and 2: Is ventricular stress (fiber and cross-fiber) reduced after plication? 3 and 4: Is ventricular stress (fiber and cross-fiber) elevated in the border zone relative to remote regions after plication?); consequently, a Bonferroni adjustment for 4 comparisons was applied to all end-systolic stress statistics. All other statistics presented were for exploratory purposes and were not adjusted. All data are presented as mean ± standard deviation.
| Results |
|---|
|
|
|---|
|
|
|
| Discussion |
|---|
|
|
|---|
Although there have been several FE studies of MI,9,13
to our knowledge none have been performed on LV aneurysm plication. Savage and coworkers27
computed wall stress in the ovine LV aneurysm model before and after plication using a closed form equation developed by Janz.8
Comparison of circumferential stress from Savage and colleagues' study with mid-wall fiber stress (which is predominantly circumferential at the mid-wall28
) from the present study indicates serious discrepancies, as much as 122% difference in the anterior aneurysm preplication and 60% difference in the anterior border zone postplication (data taken from Table 2,27
assuming the units posted are dynes/cm2 and not dynes/sec/cm2). In addition, large standard deviations exist in their stress estimates, suggesting sensitivity in the method. Taniguchi and colleagues2
used the same closed-form equation for their study and also present large variations in their stress estimates. Because of the simplifying assumptions made to develop the closed form equation and the regional validation to tagged MRI presented in the current study, we expect that our stress calculations are more accurate.
The assumption in DTMRI applied to myocardium is that the direction of the muscle fibers is collinear with the direction of the eigenvector associated with the largest eigen diffusivity.29
We believe our results were made more accurate by incorporating 3-dimensional variations in myofiber orientation in our FE models. However, a uniform transmural myofiber variation gave comparable results.
The relatively wide variance in some of the material parameters that are expressed in Table E1 is typical of soft tissues.30
Incorporating this variance allows us to equally fit the measured end-systolic strain with our FE models. Currently, we are developing and validating a user-friendly software tool that performs a formal optimization of myocardial material properties.
Stress in the Myocardium
The reason that infarct plication fails to stabilize LV volume must be that stress reduction is not sufficiently great. Stress is important because it determines myocardial energy expenditure31
and coronary blood flow,32
causes hypertrophy,10,33
and may lead to apoptosis34
and nonischemic infarct extension (NIIE), a process in which normally perfused segments adjacent to the infarct increase in size over time in response to high systolic stress.35,36
Even if the effect of a surgical ventricular remodeling operation on function is neutral, the effect on stress may make the operation worthwhile. Specifically, a reduction in systolic border zone stress may stop or reduce NIIE in the border zone, and a reduction in diastolic stress in the remote myocardium may stop or reduce eccentric hypertrophy after MI. Our FE-based calculations show that end-systolic fiber and cross-fiber stress in the border zone and residual aneurysm were reduced after infarct plication. However, preliminary end-diastolic stress calculations suggest that end-diastolic fiber and cross-fiber stresses were not reduced in the remote myocardium after plication. Consequently, we expect the primary positive effect of infarct plication to be in the infarct border zone.
However, the amount of stress reduction necessary to halt or reverse NIIE in the infarct BZ and eccentric hypertrophy in the remote myocardium is unknown. In a recent study, Moustakidis and coworkers37
reported that peak systolic stress (circumferential) in normal sheep was approximately 12 kPa (Figure 4, A). It should be noted that plication, which reduces end-systolic fiber stress to 19.2 ± 4.6 kPa, does not succeed in "normalizing" stress.
Moustakidis and colleagues37
did not calculate diastolic stress, so there is no way to say whether plication "normalizes" diastolic fiber stress. Moreover a strict comparison cannot be performed because Moustakidis and colleagues did not report the pressures in their models or validate their models against tagged MRI measurements.37
Further evaluation of end-systolic and end-diastolic stress levels in normal sheep using the methods in this study may be warranted.
Other Analytic Approaches
In the past, analytic estimates of stress have used equations based on force balance and LaPlace's law.8,10
However, these equations cannot account for varying myofiber orientation or regional material properties (eg, a stiffer aneurysm). In recent literature, ventricular stress calculations have been presented from FE models with detailed geometry and loading conditions but a linear elastic, infinitesimal deformation analysis.37,38
The advantage to this FE approach is that the computation time is much shorter because it does not require incremental loading or iterative solutions. However, this does not reflect the nonlinear behavior of myocardium or account for the large deformations of the beating heart. More than 20 years ago, Janz and Grimm39
showed circumferential strains can be off by as much as 90% with linear elasticity. Furthermore, solutions were off by as much as 100% when not accounting for large deformations.40
The FE method in the present study incorporates a large-deformation, nonlinear elastic analysis that is anisotropic relative to the local myofiber orientation (measured in detail with DTMRI). Moreover, because of the rigorous comparison with systolic strain measurements made from tagged MRI, we expect our calculations to be the best estimates of end-systolic stress in the LV after a linear aneurysm repair to date. Diastolic strain will be measured to further validate the diastolic model in the future but is not currently available.
Persistent systolic stress in the retained septal infarct may have been responsible for ongoing NIIE and the redilation seen in our previous study of aneurysm plication.5
Of note, the Dor operation41
excludes the septal infarct completely, and therefore ongoing NIIE is not possible. Nevertheless, there may still be stress concentrations around the inserted patch. FE analysis of this and similar operations is currently being pursued in our laboratory.
We chose to intervene at 13 weeks (and 16 weeks in our Dor study) because the LV has remodeled significantly at that point. From a practical standpoint, the hearts are large enough to perform cardiac operations with (Dor) and without (plication) cardiopulmonary bypass. We thought that the 2- and 6-week postplication time points were reasonable early and intermediate times after the intervention. However, long-term observation after plication and the Dor procedure in this animal model may be beneficial.
Limitations
Because the methods used are the same as in our previous study, the same limitations apply. In summary, long computation times prohibit a formal optimization of material constants25
and greater mesh refinement. Software and hardware upgrades, currently being pursued in our laboratory, will address both of these issues. Second, limited accuracy of measurements of radial strain (Err) have been acknowledged by multiple researchers as a limitation of tagged MRI26
but would be valuable for further validation of the FE analysis. It is important to note here that we calculate Err (and the other 5 strain components) in our FE model, but we do not try to fit model material properties to the experimental Err data (only to the other 5 components). This does not result in significant inaccuracy because Err can be determined from the other 5 strain components (assuming incompressibility). Third, a stress-free state in the beating heart is difficult to obtain, especially when complicated with a dyskinetic infarct. Kresh and Wechsler42
have argued that the dilated LV is devoid of residual stress. Our previous study using a cylindric model of the equatorial region of the normal canine LV23
suggests that the passive material properties would be different but not overly so if the unloaded model configuration is not stress free. Our previous study of partial left ventriculectomy43
suggests that fiber stress is not significantly affected by residual stress. Fourth, for simplicity we did not include secondary effects, such as those associated with scarring around the surgical repair (there is a lack of model restraint in the area of the aneurysm plication), Teflon felt (we did not include different material properties for this), and the right ventricle (we did not load the opposite side of the septum with right ventricular pressures or include the right ventricular attachments to the septum). Fifth, the transmural change in myofiber helix angle is rapid, and small spatial misalignments can cause inaccurate extraction of myofiber data. This is a limitation of applying ex vivo data to in vivo data. However, it is at least as good if not better than histologic sectioning in which the heart is also removed and the measured angles mapped back to in vivo locations. Last, our FE model includes only mechanical effects, whereas aneurysm plication probably has significant effects on LV electrophysiology, metabolic activity, and so forth.
| Conclusions |
|---|
|
|
|---|
| Figure E1 |
|---|
|
|
|---|
|
| Table E1 |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||
SD, Standard deviation. C = 0.290, bf = 33.17, and bt = 41.56 in the aneurysm; bfs was constant across all regions. Tmax was reduced by one half in the border zone.
| Acknowledgments |
|---|
| Footnotes |
|---|
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. F. Wenk, P. Eslami, Z. Zhang, C. Xu, E. Kuhl, J. H. Gorman III, J. D. Robb, M. B. Ratcliffe, R. C. Gorman, and J. M. Guccione A Novel Method for Quantifying the In-Vivo Mechanical Effect of Material Injected Into a Myocardial Infarction Ann. Thorac. Surg., September 1, 2011; 92(3): 935 - 941. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Zhong, Y. Su, L. Gobeawan, S. Sola, R.- S. Tan, J. L. Navia, D. N. Ghista, T. Chua, J. Guccione, and G. S. Kassab Impact of surgical ventricular restoration on ventricular shape, wall stress, and function in heart failure patients Am J Physiol Heart Circ Physiol, May 1, 2011; 300(5): H1653 - H1660. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. F. Wenk, Z. Zhang, G. Cheng, D. Malhotra, G. Acevedo-Bolton, M. Burger, T. Suzuki, D. A. Saloner, A. W. Wallace, J. M. Guccione, et al. First Finite Element Model of the Left Ventricle With Mitral Valve: Insights Into Ischemic Mitral Regurgitation Ann. Thorac. Surg., May 1, 2010; 89(5): 1546 - 1553. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Yu, K. L. Christman, E. Chin, R. E. Sievers, M. Saeed, and R. J. Lee Restoration of left ventricular geometry and improvement of left ventricular function in a rodent model of chronic ischemic cardiomyopathy. J. Thorac. Cardiovasc. Surg., January 1, 2009; 137(1): 180 - 187. [Abstract] [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 |