|
|
||||||||
a Department of Surgery, University of California, San Francisco, California
b Department of Bioengineering, University of California, San Francisco, California
c Department of Anesthesia, University of California, San Francisco, California
d Department of Radiology, University of California, San Francisco, California
e Veterans Affairs Medical Center, San Francisco, California
f Livermore Software Technology Corporation, Livermore, California
Accepted for publication February 12, 2010.
* Address correspondence to Dr Ratcliffe, Surgical Service (112), San Francisco Veterans Affairs Medical Center, 4150 Clement St, San Francisco, CA 94121 (Email: mark.ratcliffe{at}va.gov).
| Abstract |
|---|
|
|
|---|
Methods: A finite element model of the left ventricle, mitral apparatus, and chordae tendineae was created from magnetic resonance images from a sheep that developed moderate mitral regurgitation after posterobasal myocardial infarction. Each region of the model was characterized by a specific constitutive law that captured the material response when subjected to physiologic pressure loading.
Results: The model simulation produced a gap between the posterior and anterior leaflets, just above the infarcted posterior papillary muscle, which is indicative of mitral regurgitation. When the stiffness of the infarct region was reduced, this caused the wall to distend and the gap area between the leaflets to increase by 33%. Additionally, the stress in the leaflets increased around the chordal connection points near the gap.
Conclusions: The methodology outlined in this work will allow a finite element model of both the left ventricle and mitral valve to be generated using noninvasive techniques.
| Introduction |
|---|
|
|
|---|
Ischemic MR is usually managed with mitral annuloplasty alone [3, 4]. However, the need for surgical therapy and the type of surgical therapy remains unclear. For instance, the type of annuloplasty (flexible versus rigid, complete versus incomplete) is controversial [3, 4] and the role of surgical remodeling with devices including the Acorn CorCap (Acorn Cardiovascular, St Paul, MN) [5] and Coapsys (Myocor, Maple Grove, MN) [6] is unclear.
Finite element (FE) modeling of the heart and heart valves is becoming more and more common. For instance, there have been numerous FE modeling studies of the aortic valve, either as a bioprosthesis or as native tissue. In contrast, the mitral valve is structurally more complex and includes the contracting left ventricular (LV) wall as a structural element. As a consequence, most FE models of the mitral valve assume that the papillary muscles are fixed in space. Examples include the work of Kunzelman and colleagues [7, 8], who generated the first FE model of an isolated mitral valve. Finite element models of the isolated mitral valve have now been used to determine the effect of artificial chordae on valve deformation [9], the effects of leaflet and annular curvature on leaflet stress [10], and the effects of innovative annuloplasty ring shape [11, 12]. More recently, modeling of the mitral leaflets has become more realistic [13] and fluid structure interaction has been added [14].
On the other hand, FE models of the LV that incorporate myocardial contractility have been described and used to determine the effect of myocardial infarction on structure and function [15–17]. Additionally, the efficacy of various surgical procedures has been simulated with the FE method [18, 19]. Most recently, a method for determining myocardial material properties noninvasively was developed using FE modeling and optimization [20]. To date, however, FE models of the LV have not included the mitral valve.
In the current study, we expand our previous FE models of the LV [20] to incorporate the leaflets and chordae of the mitral valve. The model is based on magnetic resonance imaging (MRI) data from a sheep that has developed moderate ischemic MR after posterobasal myocardial infarction [21, 22]. We demonstrate the utility and power of the method by using the FE model to test the hypothesis that a reduction in infarct stiffness will increase dyskenesis of the posterior LV wall, increase the displacement of the posterior papillary muscle, and thereby increase ischemic mitral regurgitation.
| Material and Methods |
|---|
|
|
|---|
Myocardial Infarction
A single adult sheep underwent posterobasal myocardial infarction as previously described [21, 22]. In brief, a castrated male Dorsett sheep was anesthetized (ketamine, 1 gm, intramuscularly; isofluorane, 2 to 4% inspired) and mechanically ventilated (Model 309-0612-800; Ohio Medical Products, Madison, WI). The exhaled partial pressure of carbon dioxide was continuously monitored and tidal volume and respiratory rate adjusted to maintain a partial pressure of carbon dioxide less than 40 mm Hg. During a left thoracotomy, ligatures were placed around the second and third obtuse marginal coronary arteries immediately above the atrioventricular groove. After the thoracotomy was closed, the sheep recovered from anesthesia. Eight weeks after myocardial infarction, the sheep underwent transdiaphragmatic echo and magnetic resonance imaging (MRI) as previously described [23].
Echocardiography
A midline sub-xyphoid incision was made and subdiaphragmatic two-dimensional long axis images were obtained using an Acuson Sequoia C256 with a 3V2c (2 to 3 MHz) probe. The severity of MR was assessed quantitatively as the area of the regurgitant jet as a percentage of left atrial area in the apical four-chamber view. The following grading was used: grade 1, less than 20%; grade 2, 20 to 40%, grade 3, 40 to 60%; and grade 4, greater than 60% of the left atrial area.
Magnetic Resonance Imaging
The LV pressure was measured with a nonferromagnetic transducer-tipped pressure catheter (model SPC-350; Millar Instruments, Houston, TX) as described previously by Guccione and colleagues [24]. A series of orthogonal short-axis and long-axis MRIs were acquired (Fig 1) as described previously in Reference 24.
|
A customized version of the MRI postprocessing software, FindTags (Laboratory of Cardiac Energetics, National Institutes of Health, Bethesda, MD) was used to contour the endocardial and epicardial LV surfaces (Fig 1), and also to segment the systolic tags for each image slice [25]. The mitral valve annulus and leaflets were identified from the magnetic resonance images using the custom software Find Contours (Cardiac Biomechanics Laboratory, San Francisco, CA). Within each long-axis slice, 3 points were selected along the posterior leaflet and 4 points were selected along the anterior leaflet. The annulus was constructed from a total of 12 points around the perimeter of the mitral valve. The chordae tendineae could not be identified on the MRIs and thus could not be contoured. However, the computational reconstruction of the chordae was approximated from anatomic images of a mitral valve from an excised heart, which included the attached chordae and papillary muscles [26].
Finite Element Model
A FE model was created. Infarct, remote, and border zone regions of the LV were determined based on the LV wall thickness using Find Contours. Specifically, the border zone region was defined as the steep transition in wall thickness between the remote and infarcted regions [27]. Surface meshes were then created from the LV contours to replicate the in vivo geometry (Rapidform; INUS Technology, Inc, Sunnyvale, CA), as shown in Figure 2(A). The spaces between the endocardial and epicardial surfaces were filled with 8-node trilinear brick elements, with a single integration point for computational efficiency, to generate a volumetric mesh that was 4 elements thick transmurally (Truegrid; XYZ Scientific Applications, Inc, Livermore, CA) (Figs 3 A; B). Each region, infarct, remote, and border zone was assigned different material properties.
|
|
The chordae tendineae were composed of 2-node beam elements. The strut (basal) and edge (marginal) chordae were assigned different material properties and were connected directly between the papillary muscles and the mitral leaflets, without branching. More specifically, 24 edge chordae were attached to the nodes on the free edge of each leaflet to replicate the splay at the connection points seen in vivo, while 8 strut chordae were distributed over two small regions in the midsection of each leaflet. The connection points of the chordae on the papillary muscles were chosen to distribute the load, as seen in vivo, rather than concentrating them at a single point (Figs 3B and 4B).
|
Loading Conditions
The basal nodes of the LV were extended above the plane of the mitral leaflets and material was added in between in order to reduce the interaction of the boundary constraints and the leaflets. The extended basal nodes were fully constrained in all directions. The endocardial wall was loaded to the measured in vivo end-diastolic and end-systolic LV pressures. The leaflets were loaded in order to represent the pressure gradient between the atrium and ventricle, which is given by an inward pressure from the atrium during diastole and an outward pressure from the ventricle during systole.
Constitutive Models
Nearly incompressible, transversely isotropic, hyperelastic constitutive laws for passive [29] and active myocardium [30] were modeled in a user-defined material subroutine in the explicit FE solver, LS-DYNA (Livermore Software Technology Corporation, Livermore, CA). The key equations can be found in prior work by Walker and colleagues [31]. The mitral valve leaflets were also modeled with a nearly incompressible, transversely isotropic, hyperelastic material response. The soft tissue material model (*MAT_091), which is built into LS-DYNA [32] and is based on the constitutive law developed by Quapp and Weiss [33], was used to model the leaflet response. An overview of several widely used material models for the mitral valve is given by Weinberg and Mofrad [34]. The chordae tendineae were modeled with a cable element formulation (*MAT_071), which allows a load to be carried only in tension and does not sustain bending [32]. Rather than using a strain energy function, the uniaxial stress-strain curve was used directly in the model as a look-up table.
The material properties of the myocardium were averaged from previous studies [17, 20] and are shown Table 1. The infarct was assigned a value roughly 10 times greater than the remote, but with no contractility. The contractility in the border zone was reduced by 50%, relative to the remote region [17, 20]. The material properties for the mitral leaflets were estimated by comparing the stress-strain curves from experimental data [35] and model simulations performed with the soft tissue model using LS-DYNA. The values were determined manually, rather than using a formal optimization, and are shown in Table 2. Unfortunately, it was not possible to exactly reproduce the curves. The stress-strain relation for the chordae tendineae was based on data from Kunzelman and Cochran [36].
|
|
| Results |
|---|
|
|
|---|
The severity of MR was quantified from the echocardiograph images (Fig 5). The grade of MR was found to be 2+ in the animal used in this study. It was not possible to ascertain the presence of a gap in the images, due to low resolution. However, the presence of a regurgitant jet is indicative of an open passage between the LV and atrium.
|
|
The Effect of Papillary Muscle Displacement on Simulated Results
The infarct material properties were varied to investigate the effects on deformation in the mitral valve leaflets. When the infarct stiffness, C, was reduced from 3.6 kPa to 0.3 kPa, the wall deformation increased significantly (Fig 7). The area of the gap between the leaflets was estimated by placing a cut plane through the gap and measuring the area of the bounded region. The area of the gap was found to increase by roughly 33%, from 3.74 mm2 to 4.96 mm2, when the stiffness of the infarct was reduced, as seen in Figure 8. The stress in the leaflets also increased due to the increased wall deformation. The average effective stresses at several locations in the leaflets are given in Table 3. It is worth noting that the stress in the center of the posterior leaflet did not increase significantly during this simulation.
|
|
|
| Comment |
|---|
|
|
|---|
The Sheep Model of Ischemic MR
Ligation of the second and third obtuse marginal coronary arteries in the sheep causes MR to gradually develop over 6 weeks until it reaches a moderate (2 to 3+) level [21, 22]. Ischemic MR in the sheep is the mechanical result of leaflet teathering (Carpentier type 3B mitral leaflet motion [40]) and annular dilation. As a consequence of LV remodeling after the posterobasal myocardial infarction, the posterior papillary muscle moves laterally [41, 42]. This causes both anterior and posterior mitral leaflets to be tethered and ischemic MR to occur. The mitral annulus also dilates. For instance, Gorman and colleagues [43] found that the area of the annulus increased 9.2% at end-systole. While there are differences between human and sheep mitral anatomy [44], the mechanism of ischemic MR in the sheep model is thought to be similar to the mechanism of ischemic MR in humans.
Our model of the LV with mitral valve was able to demonstrate the presence of ischemic MR. For instance, a gap was seen between the anterior and posterior leaflet in the A3P3 region at end-systole. In addition, the model was faithful to the other anatomic characteristics of ischemic MR. Specifically, the shape of the mitral annulus in the model was relatively flat. This has been observed by Ryan and colleagues [45] in sheep with MR. The deformed shape of the leaflets showed a slight billowing effect in the center of the leaflets, while the strut chordae attachment points on the leaflets were restrained [42].
The Effect of Infarct Stiffness
We used the finite element model of LV with mitral valve to test the hypothesis that a reduction in infarct stiffness will increase ischemic mitral regurgitation. Specifically, we were able to demonstrate that a reduction in infarct stiffness caused further displacement of the posterior papillary muscle at end-systole. As a consequence, there was an increase in the gap between the anterior and posterior leaflets in the A3P3 region.
The ability of lateral papillary muscle displacement to decrease mitral leaflet coaptation and increase mitral regurgitation has been previously shown using an isolated mitral valve preparation [46] and a finite element model of the mitral valve [11]. In both of those studies, papillary muscle location was manually prescribed, whereas in our model the location of the papillary muscle was solved for by simulating the motion of the myocardial wall.
Leaflet Stress
The stress distribution varied significantly over the leaflets in the present study, due to the nonsymmetric geometry and tethering. Leaflet stress in isolated mitral valve models has been reported [13], but most of these studies apply a ventricular pressure of 120 mm Hg to models that are symmetric, whereas in the present study we apply a ventricular pressure of 91.46 mm Hg. However, it was found that the magnitude of stress in the center of the leaflets was roughly the same order of magnitude as that shown in recent studies. For example, the effective stress in the center of the anterior and posterior leaflets was reported by Prot and colleagues [13] to be roughly 220 kPa and 130 kPa, respectively. In another study, Salgo and colleagues [10] reported much higher stresses in a mitral valve model that were on the order of 4 to 400 MPa. However, it should be noted that the study was conducted with a linear orthotropic material law.
Limitations
The geometry of the mitral valve does not include the clefts in the posterior leaflet. This could affect the deformation of the leaflets near the endpoints of the free edges. Second, the leaflets were not modeled to include contact, thus there is some penetration of the anterior leaflet into the posterior. The inclusion of contact would most likely change the shape of the noncoapting region. However, the result of the model simulation would still produce a gap, indicative of MR. Third, the material model for the leaflets does not capture the rapid stiffening effect from the fibers at the critical stretch, but rather has a more gradual stiffening effect. This most likely contributes to the slight billowing in the leaflet deformation. Finally, the material properties used for the myocardium have not been validated, but are based on average values.
Conclusions and Future Directions
The methodology presented here provides a systematic framework for generating an accurate model of the LV with mitral valve. The model presented in this study was able to capture the presence of MR, which was seen as a gap between the leaflets, and to show the effect of displacement of the papillary muscle on ischemic MR.
We plan to improve the model of the LV with mitral valve in a number of ways. For instance, future models will be made more anatomically correct with the use of coregistered three-dimensional echocardiography and MRIs. Second, a contact algorithm to prevent penetration between the leaflets and implementation of a more widely used material model for leaflet stiffness will be integrated into the model. Third, the model will be expanded to incorporate the effect of blood flow and flow structure interaction. This will generate a more accurate profile of the loading applied on the endocardium and mitral valve. Last, we will employ our material optimization techniques to calculate actual myocardial and mitral annular material properties.
The success of this study opens the way for further finite element-based studies of therapy for ischemic MR. First and most obvious would be a more detailed study of infarct stiffness which would determine the infarct stiffness that optimally reduces MR and mitral leaflet stress. Models of passive cardiac constraint are being developed in our laboratory and a logical next step would be to simulate the effect of the Acorn CorCap on ischemic MR. Beyond that there are numerous possibilities, including the determination of optimal annuloplasty ring material properties and shape and the effect of Coapsys-like intraventricular struts on LV and mitral valve function.
Ultimately, the goal is to use the model to determine the optimal surgical repair of ischemic MR. We believe that as this model improves it will be a powerful tool for planning future clinical and animal trials, thus reducing the need for expensive and time consuming studies.
| Acknowledgments |
|---|
|
|
|---|
| References |
|---|
|
|
|---|
Related Article
This article has been cited by other articles:
![]() |
R. Carrick, L. Ge, L. C. Lee, Z. Zhang, R. Mishra, L. Axel, J. M. Guccione, E. A. Grossi, and M. B. Ratcliffe Patient-Specific Finite Element-Based Analysis of Ventricular Myofiber Stress After Coapsys: Importance of Residual Stress Ann. Thorac. Surg., June 1, 2012; 93(6): 1964 - 1971. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. M. Wong, J. F. Wenk, Z. Zhang, G. Cheng, G. Acevedo-Bolton, M. Burger, D. A. Saloner, A. W. Wallace, J. M. Guccione, M. B. Ratcliffe, et al. The Effect of Mitral Annuloplasty Shape in Ischemic Mitral Regurgitation: A Finite Element Simulation Ann. Thorac. Surg., March 1, 2012; 93(3): 776 - 782. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Braun Invited Commentary Ann. Thorac. Surg., May 1, 2010; 89(5): 1554 - 1554. [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 |