Abstract
Objective: To evaluate the risk of root resorption, individual finite element models (FEMs) of extracted human maxillary first premolars were created, and the distribution of the hydrostatic pressure in the periodontal ligament (PDL) of these models was simulated.
Materials and Methods: A continuous lingual torque of 3 Nmm and 6 Nmm respectively was applied in vivo to the aforementioned teeth. After extraction, FEMs of these double-rooted teeth were created based on high-resolution microcomputed tomographics (micro CT, voxel size: 35 microns). This high volumetric resolution made the recognition of very small resorption lacunae possible. Scanning electron micrographs of the root surfaces were created as well. This enabled the investigation of advantages and disadvantages of the different imaging techniques from the viewpoint of the examination of root resorption. Using the FEMs, the same loading conditions as applied in vivo were simulated.
Results: The results of clinical examination and simulations were compared using the identical roots of the teeth. The regions that showed increased hydrostatic pressure (>0.0047 MPa) correlated well with the locations of root resorption for each tooth. Increased torque resulted in increased high-pressure areas and increased magnitudes of hydrostatic pressure, correlating with the experiments.
Conclusion: If hydrostatic pressure exceeds typical human capillary blood pressure in the PDL, the risk of root resorption increases.
INTRODUCTION
Minimization of the risk of root resorption is important during orthodontic treatment. The biologic mechanisms of this process have not been fully studied but are known to be strongly dependent on individual factors of the patient. A multitude of factors influencing root resorption have been discussed. Genetic factors,1 the mineral content of cementum,2 and orthodontic forces and moments, as well as their duration during application, are believed to influence susceptibility to root resorption. Recent studies show that the extent of root resorption depends on the applied force system.3–6 Schwarz7,8 stated that root resorption can occur if capillary blood pressure is exceeded. In the literature, the range of capillary blood pressure is reported to be between 15 and 35 mm Hg (equivalent to 0.0020 to 0.0047 MPa).9,10
Casa et al4 applied torque to human premolars in vivo, and after extracting them, they investigated the root resorption that occurred. They observed resorption lacunae on the linguoapical portion of the root and on the buccocervical third of the root.
The goal of the present study was to simulate these experiments based on the same teeth. The distribution of the hydrostatic pressure in the periodontal ligament (PDL) of these double-rooted premolars was compared in a finite element model (FEM) with the results of the in vivo experiments of the same teeth for validation of the capillary blood pressure theory of Schwarz. It was assumed that a stress greater than blood pressure could make the capillary blood vessels in the PDL collapse.11
Stress in the PDL has been analyzed in the literature.9,10,12–14 In single-rooted teeth stress maxima were found near the apex. To our knowledge, there are no existing studies using the identical teeth for simulation of orthodontic torque and comparing the results with root resorption on these teeth.
To date, root resorption and the mineralized structure of cementum have been examined only by scanning electron microscopy (SEM).3–5,15–17 This is a rather difficult means of examining root resorption, as sputtering is necessary and the desired image view must be chosen before making an image. In this study, root resorption lacunae and their morphology were visualized by SEM and microcomputed tomography (micro CT) to show the advantages and disadvantages of these different visualization techniques.
MATERIALS AND METHODS
The extracted teeth used for this study came from an in vivo study performed by Casa et al.4 They investigated root resorption in human maxillary first premolars after application of continuous torque moments (3 Nmm and 6 Nmm), using superelastic nickel-titanium–stainless steel springs for the application of the torque moments.
The continuity of the moments over a treatment time of 4 weeks was induced by temperature conditioning the wire at the Department of Orthodontics, Paulista University, São Paulo, Brazil.
After 4 weeks, the premolars were extracted. Clinically, the patients required the extraction of these teeth as part of their orthodontic treatment. All patients were informed about the procedures and their written consent was obtained. The study was authorized by the Medical Ethics Committee of the University of São Paulo.
The root surfaces were examined by SEM.4 SEM is a rather difficult means of examining root resorption. Therefore, micro CT datasets (Institute of Orthopaedic Research and Biomechanics, Ulm, Germany) were created for the aforementioned teeth (3 Nmm torque and 6 Nmm torque) and compared with the SEM results. The resolution of the dataset was 35 microns. The density information in the micro CT dataset was used to create isosurface representations of the surface of the tooth roots. Figure 1 shows such an isosurface view of the tooth treated by a torque of 6 Nmm in the experiment.
Model generation and preprocessing for finite element analysis was performed by AMIRA (Mercury Computer Systems, Berlin, Germany), SURFACER (Ann Arbor, Mich) and ANSYS (Canonsburg, Penn), according to the method published by Clement et al.18
The PDL was constructed as a layer with a thickness of 0.3 mm around the tooth root. A sensitivity analysis with a single-rooted premolar showed that this was a justifiable approximation. The number of elements forming the tooth roots is shown in Table 1. This number differs between both tooth roots because individual meshes were created for the 3-Nmm tooth root and for the 6-Nmm tooth root.
Final preprocessing, finite element calculations, and postprocessing were performed in ABAQUS (Providence, RI). A linearly elastic constitutive equation was chosen for the PDL and the tooth (Table 2). Because resorption generally begins after a few weeks of force application, the time factor did not have to be considered. Thus, the goal is not to investigate the instantaneous reaction of the tooth to loading, but the stress distribution in the position of equilibrium. Figure 2 illustrates the loading conditions of the numeric simulation of a torque also successfully performed by Dorow and Sander.9
Two forces working in opposite directions were applied to two nodes on the lingual and labial sides of the root near the alveolar crest. They resulted in a total lingual moment of My = 3 Nmm or My = 6 Nmm, respectively, and were used to simulate lingual torque. This approach was justified because only the relative position of the forces is important for simulation of a torque.
RESULTS
Figure 3 shows an SEM image of the tooth under a torque of 6 Nmm. This high-resolution representation of the apical root part gives detailed information about the microscopic structure of the resorption lacunae. Very large resorption lacunae are visible in the lingual surfaces of the apical parts of the root.
Figure 1 shows an isosurface view of the apical part of the same tooth based on micro CT. The morphology of the roots with resorption is clearly visible, but the resolution is too low for representation of the microscopic structure and the mineralized structure of cementum.
In this study, the hydrostatic pressure (p), which is equal to- σH = (σ1 + σ2 + σ3)/3 with hydrostatic stress σH and the principal stresses σ1, σ2 and σ3, was calculated. In Figures 4 (3-Nmm tooth) and 5 (6-Nmm tooth), the pressure on the surface of the PDL is shown. This pressure also represents the pressure within the PDL. A test showed that pressure does not change significantly within the PDL in directions tangential to the surface of the PDL.
After 4 weeks of torque application (3 Nmm), Casa et al4 observed resorption only at the lingual part of the two apical root parts. They also found some resorption lacunae at the buccocervical third of the root.
This corresponds to the pressure zones (greater than 0.0047 MPa, maximum literature value for capillary blood pressure9) resulting from the simulation of the same torque application to the same root (light grey areas in Figure 4).
The finite element calculations with the same root show regions of pressure exceeding 0.0047 MPa at the buccocervical part of the root and overall lingual surfaces of the two apical root parts. There were no other resorption zones found in the experiments, and no other pressure zones exceeding 0.0047 MPa resulted from the simulations.
For a lingual torque of 6 Nmm, more and deeper resorption lacunae were observed linguoapically and buccocervically.
Figure 3 shows an SEM of the apical part of the root after 4 weeks of application of 6 Nmm of torque. Figure 5 shows the results of the finite element calculations with the same root and the same torque. Pressures above 0.0160 MPa are represented in grey. Pressures exceeding 0.0047 MPa can be observed on the lingual sides of the two apical root parts. Buccocervically, the simulations also showed an area with a pressure that exceeded 0.0047 MPa.
These high-pressure zones correspond to the zones where root resorption occurred in the experiment. The area of the buccocervical high-pressure zone is greater than the areas of the linguoapical high-pressure zones. The pressure maximum in the buccocervical high-pressure zone is also greater than the pressure maxima of the linguoapical high-pressure zones.
DISCUSSION
SEM is a rather difficult method of examining tooth roots for resorption because the desired image view must be chosen before an image is created. By contrast, during micro CT imaging, a dataset of the whole tooth is recorded and sputtering is not necessary. During computer-aided image processing, the investigator can choose any image view desired. Furthermore, it is possible to investigate the resorption lacunae volumetrically, as described by Geiger et al.19
Although SEM gives much more detailed information about the resorption lacunae and the mineralized structure (Figure 3), the micro CT dataset gives the morphology of the root with resorption lacunae, which makes identification of the major resorption lacunae easier (Figure 1). If only morphologic information about the resorption lacunae is required, micro CT can be an interesting alternative.
In regions where the capillary blood pressure was not exceeded, no root resorption occurred. This indicated that there seems to be no risk of root resorption in regions where the capillary blood pressure is not exceeded. In regions of the tooth root where the capillary blood pressure exceeded 0.0047 MPa (buccocervically and linguoapically), root resorption was observed in this study.
After application of a torque of 6 Nmm, the root resorption was greater than for the teeth that received 3 Nmm of torque in the in vivo experiments. This agrees with the simulations, in which more extended high-pressure zones were observed for the 6-Nmm teeth than for the 3-Nmm teeth. The hydrostatic pressure in the PDL obviously depends on the magnitude of the force being applied. This contradicts the findings of Owman-Moll et al20,21 and Maltha et al22 but is in agreement with the findings of Faltin et al,3,6 Casa et al,4 and Chan and Darendeliler,5 who observed that the extent of root resorption depended on the magnitude of the applied force system.
Simulations performed by Dorow and Sander9 of a single-rooted mandibular first premolar with a lingual torque show that increased hydrostatic pressure can be observed at the linguoapical and the buccocervical parts of the PDL. This observation was confirmed in this study of double-rooted teeth. There were no other high-pressure zones apart from those shown in the simulations.
For lingual torque, the highest pressures and the largest high-pressure areas were observed at the buccal part of the PDL near the alveolar crest. One would expect root resorption that is at least as severe as in the apical part, but this was not observed in vivo. This indicates that a difference in dentoclast activity may reduce the sensitivity of the root to resorption in the cervical third.
The different blood sources for vascularization of the PDL may also play a part, depending on the position in the PDL considered. The cervical third is also connected with the supraperiosteal arteries, while the apical third of the PDL is near the arteria dentalis.23 The distribution of capillary blood pressure along the PDL is not known but might also contribute to the differences.
Apical cementum is softer than cervical cementum and there are fewer Sharpey's fibers, which may play a part in resorption susceptibility.24,25 Henry and Weinmann26 discovered that the apical third of the root is more susceptible to root resorption than the cervical third. Because the cervical acellular cementum is more mineralized than the apical cellular cementum, it is possible that the mineral content influences susceptibility to resorption.
An increase in hydrostatic pressure and increases in high-pressure areas (higher than 0.0047 MPa) corresponded with more severe root resorption. Preceding calculations simulating intrusion showed the same results.
Further studies are necessary to determine whether the increased magnitude of hydrostatic pressure or the increased area of the high-pressure zones is responsible for the more severe root resorption during increased torque. The dysfunction of blood supply in larger PDL regions may be a cause of more severe root resorption. An improvement of the calculation method is required for a better understanding of root resorption.
In furcations of the roots, the simulations as well as the in vivo experiments showed no sign of local pressure maxima or resorption lacunae, respectively. This is in contradiction to the intrusion of teeth, during which an increased risk for root resorption near the furcation of the roots3 is present.
Areas with stress in the simulations (buccoapically and linguocervically) did not show root resorption in vivo. This is in disagreement with Han et al,16 who observed root resorption during the extrusion of teeth. An explanation could be the different application systems and application periods used in each study.
The individual resorption lacunae and their typical hemispheric shape17 could not be found by these simulations, but pure geometric causes can be excluded as an explanation. Perhaps micromechanical properties resulting from distribution of Sharpey's fibers in root cementum play a part. Maxima of hydrostatic pressure were observed in the apical thirds of the roots, in agreement with the findings of Dorow and Sander9 and Rudolph et al27 (studies with single-rooted teeth).
CONCLUSIONS
It is possible to examine root resorption by micro CT imaging if only morphologic information is required.
If the hydrostatic pressure rises above 0.0047 MPa, a decisive condition for root resorption is fulfilled, but other factors such as the location of the pressure on the root may also influence the severity of resorption.
Hydrostatic pressure is one of the indicators of the regions in which root resorption can occur.
Individual factors have an additional effect on root resorption.
Acknowledgments
The authors thank the Deutsche Forschungsgemeinschaft for financial support of project 272/6-1 and 272/7-1.