Next Article in Journal
Title: Is Celiac Trunk Revascularization Necessary After High-Flow Pancreaticoduodenal Arterial Arcades Aneurysm Retrograde Embolization?
Previous Article in Journal
The Effect of Bio-Electromagnetic Energy Regulation Therapy on Erectile Dysfunction in Patients with Multiple Sclerosis: A Triple-Blind Randomized Clinical Trial
Previous Article in Special Issue
Effectiveness of Activated Sodium Hypochlorite Irrigation by Shock Wave-Enhanced Emission Photoacoustic Streaming, Sonic and Ultrasonic Devices in Removing Enterococcus faecalis Biofilm From Root Canal System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Periodontal Breakdown, Orthodontic Movements and Pulpal Ischemia Correlations—A Comparison Between Five Study Methods

by
Radu-Andrei Moga
1,*,
Cristian Doru Olteanu
2,* and
Ada Gabriela Delean
1
1
Department of Cariology, Endodontics and Oral Pathology, School of Dental Medicine, University of Medicine and Pharmacy Iuliu Hatieganu, Str. Motilor 33, 400001 Cluj-Napoca, Romania
2
Department of Orthodontics, School of Dental Medicine, University of Medicine and Pharmacy Iuliu Hatieganu, Str. Avram Iancu 31, 400083 Cluj-Napoca, Romania
*
Authors to whom correspondence should be addressed.
J. Clin. Med. 2024, 13(23), 7062; https://doi.org/10.3390/jcm13237062
Submission received: 7 November 2024 / Revised: 20 November 2024 / Accepted: 21 November 2024 / Published: 22 November 2024
(This article belongs to the Special Issue Surgical and Non-surgical Endodontics in 2024 and Beyond)

Abstract

:
Background/Objectives: This study assessed the biomechanical behavior of dental pulp and the neuro-vascular bundle/NVB as well as the ischemic risks during orthodontic movements in a gradual horizontal periodontal breakdown, using five methods and aiming to identify the most accurate one. Methods: Seventy-two models of second lower premolar (from nine patients) were subjected to 3 N of intrusion, extrusion, rotation, tipping, and translation. Five numerical methods, Tresca, Von Mises/VM, Maximum and Minimum Principal, and hydrostatic pressure were used in a total of 1800 numerical simulations. The results were color-coded projections of the stress areas that were then correlated with maximum physiological hydrostatic pressure/MHP and known clinical biomechanical behavior. Results: During periodontal breakdown, all five methods displayed, for all movements, quantitative stresses lower than MHP, suggesting that 3 N are not inducing any local tissular ischemic risks for the healthy intact tissues. All five methods displayed rotation as the most stressful movement during periodontal breakdown, while translation was the least. The NVB was more exposed to ischemic risks than dental pulp during the periodontal breakdown due to constant tissular deformations. Only VM and Tresca methods showed translation as more prone to expose dental pulp (both coronal and radicular) to ischemic risks (than the other movements) during the periodontal breakdown simulation. However, all five methods showed intrusion and extrusion as more prone to expose the NVB to higher ischemic risks than the other movements during the periodontal breakdown simulation. Conclusions: During periodontal breakdown, Tresca and Von Mises were more accurate, with Tresca being the most accurate of all.

1. Introduction

There is a continuous, increasing interest in the individualized study of each dental tissular component, manifested through new numerical studies [1,2,3,4]. Dental pulp and the neuro-vascular bundle/NVB are one of the least studied among the dental tissues despite their great importance during orthodontic treatment [5,6,7,8,9,10,11,12,13]. Under orthodontic loadings, various amounts of local circulatory disturbances appear, triggering the orthodontic movements [5,6,7,8,9,10,11,12,13,14]. The physiological tissular local maximal hydrostatic pressure/MHP is 16–22 KPa (about 80% of the systolic pressure) [15,16,17,18,19,20,21,22]. Any pressure exceeding this, induced by the orthodontic loads, will lead to a variable degree of changes in the pressure amount, with a partial collapse of circulatory vessels and localized ischemia (triggering the orthodontic movement) [23,24,25,26,27,28,29,30,31,32,33]. However, if the amount of force is too high or applied for longer periods of time it will locally induce prolonged ischemia followed by necrosis and degenerative-resorptive mechanisms [23,24,25,26,27,28,29,30,31,32,33]. Despite being clinically accepted that light orthodontic forces of up to 1 N are safe for the intact periodontium [15,16,17,18,34], there are multiple studies regarding the optimal amount of force for intact periodontium, with various contradicting results, reporting variable 1–3 N and various movements as more tissular stressful (e.g., rotation vs. intrusion), creating more confusion [35,36,37,38,39,40] and no consensus. However, there is no available data regarding the amount of force safely applied in various amounts of reduced periodontium, despite being clinically quite frequent.
It must be emphasized that healthy intact tissues hold the ability to adapt to these local changes and naturally avoid any necrosis or resorptive processes [23,24,25,26,27,28,29,30,31,32,33]. However, if the living tissues were previously traumatized/injured (e.g., occlusal trauma for the NVB [9,10,11,12,13] and direct/indirect pulp capping for dental pulp [23,24,25,26,27,28,29,30,31,32,33]), some morphological structural changes occur with diminishing the tissular adaptability potential [9,10,11,12,13,33,41,42,43,44,45,46,47,48]. Thus, an orthodontic force that in healthy intact tissues will pose no ischemic risks, in a diminished tissular adaptability situation could induce ischemia, followed by necrosis and further tissular loss. Unfortunately, these previous injuries and functional tissular changes are not clinically directly visible (but can be predicted); thus, in the oral cavity only the ischemic results will be seen. Thus, an accurate method to anticipate these ischemic risks is of extreme importance for planning orthodontic treatment and predicting its outcome.
Moreover, the biomechanical tissular individual behavior is variable at various levels of periodontal breakdown [9,10,11,36,37,39,40,42,43,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64] with great influence over the tissular stress distribution and directly influencing the localized circulatory disturbances in pulp and the NVB. Most of the orthodontic patients have various amounts of bone loss (usually of a few mm), with direct influence over the amount and local orthodontic stress distribution. There are no studies to be found regarding the biomechanical behavior of dental pulp and the NVB during the periodontal breakdown. Moreover, there are no comparative studies about the accuracy of study methods. We previously approached some of these above-mentioned issues but only for the light orthodontic forces, but no data are available for larger forces [15,16,17,18,19,20,21,22].
Clinical studies are not able to individually study each dental component but can provide a general picture of the biomechanical behavior of the entire tissular structure (i.e., tooth with pulp, the NVB, periodontal ligament and bone) [15,16,17,18,19,20,21,22]. The only method enabling the individual analysis of each structure is the numerical method based on radiological CBCT images (cone beam computed tomography) [1,2,3,4]. Nevertheless, due to the small dimensions of these dental tissues (pulp and the NVB), the radiological image must have a small voxel size of less than 0.01 mm [15,16,17,18,19,20,21,22].
The numerical methods used in dental studies were originally designed for the engineering field, where they are still extremely accurate and reliable (and widely used) [64]. However, since they were imported into the dental field, their accuracy was dramatically reduced (even for the new studies as well [1,2,3,4]), providing reports that often contradicted the clinical data, varying from one report to another and creating debatable results and a mistrust of the method [15,16,17,18,19,20,21,22]. The most studied were the PDL and dental implant–bone interfaces, while the newer studies are extended to the entire oral cavity [1,2,3,4]. These accuracy issues (never addressed and resolved) were directly linked (as reported by our previous reports [15,16,17,18,19,20,21,22]) to the misunderstanding of the numerical studies methodological protocol. The numerical study employs a certain methodology based on the analyzed type of material (brittle, ductile and liquid/gas), an extremely accurate model, as well as proper boundary assumptions [15,16,17,18,19,20,21,22,64]. None of these were approached in the current numerical dental studies research flow [15,16,17,18,19,20,21,22,64]. Most of the studies [35,36,37,38,39,40,49,50,51,52,53,54,56,57,58,59,60,61,62,63,64,65,66] employed (still continuing the same pattern [1,2,3,4]) a numerical method without discussing if it is suitable for the investigated tissues [15,16,17,18,19,20,21,22]. Moreover, correlations with the physiological MHP values are rarely seen [15,16,17,18,19,20,21,22]. The still employed numerical methods [35,36,37,38,39,40,49,50,51,52,53,54,56,57,58,59,60,61,62,63,64,65,66] are Maximum (tensile) and Minimum (compressive) Principal, Von Mises/VM (overall), Tresca (shear, only recently introduced) and hydrostatic pressure (liquids/gas), despite the general acceptance that dental tissues are of ductile resemblance with a certain brittle flow mode [15,16,17,18,19,20,21,22]. Thus, the reported results contained multiple inconsistencies [35,36,37,38,39,40] that contradicted the clinical reality [5,6,7,8]. In previous research, Tresca along with Von Mises was reported for light orthodontic forces to be the most accurate [15,16,17,18,19,20,21,22]. No similar data are available for large forces and various levels of bone loss [15,16,17,18,19,20,21,22]. Additionally, some concerns were also raised due to reports regarding the poor quality of multiple in vivo studies [5,7,42] used for indirect validation.
Since numerical methods are the only possibility to analyze dental tissue [15,16,17,18,19,20,21,22,64] (e.g., pulp and NVB structures), an accurate working method is needed, but only based on a scientific comparison between multiple methods.
Our aim was to assess the biomechanical behavior of dental pulp and the neuro-vascular bundle as well as the induced ischemic risks during five orthodontic movements and 3 N in a gradual horizontal periodontal breakdown of 1–8 mm, using five numerical methods. Additionally, another objective was to analyze the differences between the reports of these five numerical methods and to identify the most accurate of all in the study of pulp and the NVB.

2. Materials and Methods

This study is part of a larger research project conducted in a stepwise manner (with clinical protocol 158/02.04.2018) and studying the biomechanical behavior of dental tissues in intact and reduced periodontium [15,16,17,18,19,20,21,22]. This study focused on assessing the most suitable method of study to obtain the most accurate results.
Our study involved nine patients (4 males/5 females, mean age 29.81 ± 1.45) and totaled 1800 numerical simulations. The simulations were performed on 72 3D models of a healthy intact second lower premolar.
The inclusion criteria were as follows: no malposition and missing teeth in the investigated region, intact healthy second lower premolar, limited to 1–2 mm bone loss around the second lower premolar, orthodontic treatment indication and non-inflamed periodontium associated with proper oral hygiene. The exclusion criteria were as follows: more than 2 mm bone loss, particular root geometry, various crown and root abnormalities, radiologically visible bone defects and poor oral hygiene after inclusion.
The region of interest held the lower first molar and two premolars and was radiologically investigated by a CBCT (cone beam computed tomography—ProMax 3DS, Planmeca, Helsinki, Finland; voxel size of 0.075 mm). For manual tissular reconstructions, the AMIRA 5.4.0 (Visage Imaging Inc., Andover, MA, USA) was employed. The manual reconstruction process identified on DICOM images all tissular components (enamel, dentine, dental pulp, the neuro-vascular bundle, the periodontal ligament/PDL and trabecular and cortical bone). All above-mentioned tissues were assembled into 3D models (Figure 1). The PDL had a physiological variable thickness of 0.15–0.225 mm. The cementum component could not be properly separated from dentine; thus, based on the physical similarity between their properties, cementum was reconstructed as dentine (Table 1). The molar and the first premolar were not reconstructed, while their alveolar sockets were filled with trabecular and cortical bone. The missing bone and PDL were reconstructed for obtaining intact periodontium models. A base of a stainless-steel bracket was reconstructed on the vestibular side of the enamel.
For obtaining a gradual horizontal periodontal breakdown, each of the nine intact periodontium models were subjected to a gradual loss of 1 mm, obtaining thus various bone loss models (1–8 mm loss).
The mesh models had 5.06–6.05 million C3D4 tetrahedral elements, 0.97–1.07 million nodes and global element size of 0.08–0.116 mm (Figure 1). There were no mesh errors, but only a reduced number of element warnings (e.g., for one of the models: pulp–NVB mesh with 4 element warnings, meaning 0.0158% for a total of 25,252 elements—Figure 1K, and tooth mesh with 39 element warnings, meaning 0.00589% for a total of 661,137 elements—Figure 1L). A mesh element warning means a small non-essential areas discontinuity, but with no impact over the biomechanical behavior of the model.
The numerical analysis employed the ABAQUS 6.13-1 (Dassault Systèmes Simulia Corp., Maastricht, The Netherlands) and used the following five most used numerical methods/failure criteria in dental studies: Von Mises (maximum overall), Tresca (maximum shear), Maximum Principal (maximum tensile), Minimum Principal (maximum compressive) and hydrostatic pressure (liquids/gas). As boundary assumptions, the base of the models had zero displacements, and isotropy, linearity, and homogeneity/non-homogeneity (as in most of the current numerical studies). Five orthodontic movements were simulated: extrusion, intrusion, rotation, tipping and translation and an orthodontic force of 3 N. The results were displayed as color-coded projections of various intensities (red-orange for high, yellow-green for moderate and blue for low) showing the tissular deformations induced by the orthodontic movements. Correlations with the maximum physiological hydrostatic pressure of 16–22 KPa as well as with the known clinical biomechanical behavior were performed to assess the ischemic risks.

3. Results

All five methods analyzing the pulp and NVB biomechanical behavior (Figure 2, Figure 3 and Figure 4 and Table 2) during the five orthodontic movements and gradual horizontal periodontal breakdown showed that the maximum amount of stress was present in the NVB, while the pulpal stress was extremely small. Thus, the NVB stress was ranging between 1.17 and 8.87 KPa, while the pulpal stress was of 0.16–0.48 KPa, for 8 mm bone loss (Table 2), lower than the 16–22 KPa of the local maximum physiological hydrostatic pressure.
It seems that 3 N of force induces no ischemic risks during orthodontic movements and periodontal breakdown (up to 8 mm tissular loss) for healthy intact teeth. However, since the NVB stress was 7.3–18.47 times higher than the pupal stress (for 8 mm loss), it is more exposed to ischemic risks than the pulp. The highest amount of stress was present (in each of the five methods) at 8 mm of loss, showing that the periodontal breakdown induces a direct increase in tissular stress (especially visible in the NVB).
Rotation movements (in each of the five methods) displayed the highest amount of stress (Table 2), seeming to be the most stressful movement, closely followed by tipping and intrusion/extrusion, while the translation was the least of them.
The hydrostatic pressure method/failure criteria displayed the highest amounts of stress, closely followed by the Maximum (tensile) and Minimum (compressive) Principal, while the Von Mises (overall) and Tresca (shear) were the least of them.
In all movements and bone loss levels, the coronal and radicular pulp displayed similar qualitative stress (with two exceptions for rotation and translation for VM and Tresca methods), suggesting that the orthodontic force induced limited tissular stress.
The VM and Tresca methods are the only methods to display localized areas of coronal pulpal stress (more pronounced for translation) (Figure 2, Figure 3 and Figure 4). For the 1 mm bone loss level, during translation, the coronal pulpal stress displayed the highest spread on vestibular, mesial and distal sides of coronal pulp (Figure 2D,E) but with no radicular pulp stress. At a 4 mm loss, the stress areas extend to the radicular pulp (Figure 3D,E) in the cervical and middle thirds of the root. An 8 mm periodontal loss displayed limited coronal pulpal stress (Figure 4D,E), with an increase in the middle third of the root (correlated with the bone loss level). The rotation movement displayed a localized reduced coronal pulpal stress at the 1 mm periodontal breakdown, completely disappearing at 4 mm loss. Based on the above, despite being the least stressful among the five movements (safe for healthy intact tissues), the translation seems to directly influence the coronal pulp, inducing potential ischemic risks at this level (especially on previously injured tissues—direct/indirect pulp capping).
Quantitatively, all five methods displayed a similar pattern of increased stress (Table 2), the only difference being their amounts. Thus, hydrostatic pressure and Maximum and Minimum Principal methods displayed quantitative results 7–18 times higher than Von Mises and Tresca, but with no possible scientific evaluation of their accuracy (since all were lower than MHP).
However, qualitatively, there are significant differences between their biomechanical behavior stress displays, which can provide an orientation to their accuracy (Figure 2D,E, Figure 3D,E and Figure 4D,E). Thus, Tresca and Von Mises methods/criteria are the only ones to display the same constant stress distribution pattern during the entire periodontal breakdown and for all five movements (constant extremely small intensity blue shades amounts of stress for pulp and various amounts of higher/medium intensity stress for the NVB). This stress display pattern is correct when compared with the known clinical tissular biomechanical behavior, and it is due to the correct design (based on the analyzed type of material) of these two methods for ductile-like materials (dental tissues being materials of ductile resemblance). When comparing the two criteria, another aspect must be addressed, respectively, the homogeneity/non-homogeneity nature of the analyzed material. The Von Mises method was designed for homogenous materials and the Tresca method for non-homogenous materials, while the living dental tissues are non-homogenous materials. Thus, despite qualitatively showing no differences, quantitatively the Tresca amounts of stress are a little bit higher when compared with Von Mises (1.12–1.17 times, falling into the literature average interval of 1.15—Table 2). Based on the above, Tresca seems quantitatively more accurate than Von Mises.
The hydrostatic pressure method/criteria displayed various inconsistencies (Figure 2A, Figure 3A and Figure 4A) in the intensity of stress display with significant variables from one movement to another and from one bone loss level to another. This biomechanical behavioral stress display is due to the failure criteria specially designed for liquids and gas where there are no shear stresses and an internal micro-architecture fundamentally different from that of the studied dental tissues (despite the high percentage of contained water/liquid, which is the main reason for it to be employed in the first place).
The Maximum (tensile) and Minimum (compressive) stress methods/failure criteria displayed the same qualitative inconsistencies (Figure 2B,C, Figure 3B,C and Figure 4B,C) as hydrostatic pressure, with various intensities of color-coded stress displays, different from one movement to another and from one periodontal breakdown level to the other. This behavioral display is because both criteria were specially designed for brittle-like materials while the dental tissues are known to have a ductile resemblance.
These three methods displayed both qualitative biomechanical behavioral pattern inconsistencies (variable mixes of negative–positive stresses, with no scientific explanation, since no compressive-tensile stresses should be present in these enclosed/protected areas) contradicting the physical known structural clinical behavior, as well as higher amounts of stress (when compared to the first two methods), seeming to be less accurate for the study of pulp and the NVB. The similar color-coded stress display in both pulp and the NVB suggests comparable tissular ischemic risks despite the anatomical topography. Biomechanically, the pulp is protected by dentine (pulp chamber and radicular root canals), while the NVB is held in the apical third of the periodontal ligament (directly exposed to stresses). Thus, it seems that the qualitative display contradicts the quantitative display (with the NVB higher than pulpal stress—Table 2) as well as clinical knowledge.
Nevertheless, despite the above-mentioned differences, all five methods similarly displayed the NVB tissular deformations induced by the orthodontic movements during the periodontal breakdown, seeming to provide a correct picture of the clinical biomechanical behavior for larger tissular deformations/effects. Despite all quantitative results being lower than MHP (i.e., 3 N are safely applied in healthy intact tissues), the displayed NVB tissular deformations could have a potential ischemic risk on previously injured tissues (e.g., occlusal trauma).

4. Discussion

Despite the importance of the NVB and dental pulp during orthodontic treatment, there is a limited number of numerical studies (some of them belonging to our team [15,16,17,18,19,20,21,22]) with a focus on these tissues, due to segmentation and 3D reconstruction difficulties. Most of these studies are for intact periodontium, and only one focuses on periodontal breakdown [18]. No other study performing a correlation between the five methods addressing the tissular biomechanical behavior issues was found in the current research flow (except ours); thus, this study is the first. However, since the NVB is held in the apical third of the periodontal ligament (absorption-dissipation function [15,16,17,18,19,20,21,22]), correlations with numerical PDL studies as well as with MHP (16–22 KPa) can be made.
Our numerical study assessed the biomechanical behavior of dental pulp and NVB during 1–8 mm periodontal breakdown totaling 1800 numerical simulations. Tresca, Von Mises, Maximum and Minimum Principal and hydrostatic pressure methods were all used by the current dental numerical studies [1,2,3,4,35,36,37,38,39,40,49,50,51,52,53,54,55,56,65,66]. The Tresca method was introduced in the current research flow by our previous research [15,16,17,18,19,20,21,22], in order to develop and provide an accurate method for the study of dental tissues [64].
Each of the five above-mentioned methods have been developed in the engineering field for the study of different types of materials (with acknowledged accuracy), an issue of extreme importance when studying dental tissues [64]. These methods had been specifically designed for a certain type of studied material [64]. The Maximum (tensile) and Minimum (compressive) Principal had been designed for brittle-like materials, that when subjected to loads suffer from non-recoverable form deformations, leading to small deformations and rapidly followed by cracking and destruction (e.g., glass and concrete) [15,16,17,18,19,20,21,22]. The Tresca (shear) and Von Mises (overall) were designed for ductile-like materials, that when subjected to loads suffer recoverable form deformations, totally recovering the form when the load stops (e.g., rubber and steel) [15,16,17,18,19,20,21,22]. Hydrostatic pressure was specially designed for liquids and gas (no shear stress), materials with no internal structure [15,16,17,18,19,20,21,22].
In numerical studies, the employed method is selected based on the analyzed material type (ductile, brittle or liquid/gas) to provide accurate results [64]. Dental tissues are ductile-like materials with a certain degree of brittle flow mode [15,16,17,18,19,20,21,22,64]. None of the previous numerical studies (except ours) had this approach, thus, reporting results with serious accuracy issues. Another issue is the natural non-homogeneity of dental tissues; thus, the employed method must fulfill this request (e.g., only Tresca was specifically designed for non-homogenous materials) [15,16,17,18,19,20,21,22,64]. On the other side, dental studies [1,2,3,4,35,36,37,38,39,40,49,50,51,52,53,54,55,56,65,66] used forces up to 6 N (extremely small when compared with the engineering field), and the tissular deformations and displacements were extremely small during the orthodontic movements; thus, differences when employing one of the five methods are less visible than when displaying the biomechanical results side by side (as in this study). The dental pulp is held in the dental pulp chamber, an incompressible hard tissue cavity made of dentine, with a protection role, where there are virtually no deformations and displacements. On the other hand, the NVB is held in the apical third of the PDL, an extremely deformable and tensile-compressive tissue, with nutritional, metabolic and absorption-dissipation roles. Thus, the NVB is extremely prone to deformations and displacements, within the PDL deformations limits, and subjected to ischemic risks. It must be emphasized that this study is the first to address the benefits of this approach in the current research flow.
When the biomechanical behavioral results of the five methods are displayed side-by-side and compared, some general patterns (in agreement with known clinical data) are similar: the NVB stress is always higher than pulpal stress (quantitatively), while the NVB tissular deformations induced by the movements in various bone loss levels are clearly visible, in agreement with known clinical tissular biomechanics. However, when further investigating the side-by-side qualitative display, the accuracy issues became more and more visible (all related to the misused method) [15,16,17,18,19,20,21,22]. Based only on the quantitative results (all lower than MHP, see Table 2), it is difficult to distinguish a scientifically based difference (although two of the criteria showed very low values compared to the other three—in line with other studies [15,16,17,18,19,20,21,22]).
Qualitatively, Tresca and VM displayed a constant, similar pattern of color-coded blue small intensity stress for dental pulp for all five movements during the periodontal breakdown simulation, which is in agreement with the known tissular clinical biomechanical behavior and other reports [5,7,9,10,11,12,33,34,42,43,44,45,46,47,48,64], as well as our previous research [15,16,17,18,19,20,21,22]. The quantitative values were also the smallest among the five methods, 5.2–145.5 times under MHP, posing no ischemic risks for the 3 N of applied force, which is also in agreement with previous reports and the above-mentioned clinical data. A previous study [18], with similar methodology using 0.5 N of force, reported similar correlations as this study, with VM and Tresca being more accurate, with Tresca better suited for dental tissues. However, translation, despite being the least stressful movement among the five, displayed coronal pulp and radicular stress in both methods (Figure 2D,E, Figure 3D,E and Figure 4D,E) during the entire horizontal periodontal breakdown simulation. Despite being quantitatively lower than MHP, with no ischemic risks for healthy tissues, if pulpal tissues suffered of an earlier injury (i.e., direct/indirect pulp capping [9,10,11,12,33,42,43,44,45,46,47,48]), ischemic risks are expected (in agreement with clinical data [13,41]).
The tissular NVB deformations reported by all five methods, are in line with known clinical biomechanical behavior and other reports [5,7,9,10,11,12,33,34,42,43,44,45,46,47,48,64]. All five methods reported quantitative amounts of stress lower than MHP, thus posing no ischemic risks for healthy intact tissues. However, there are reports of ischemia and necrosis in previously injured tissues (i.e., occlusal trauma [9,10,11,23,24,25,26,27,28,29,30,31,32,33]); thus, our study is in agreement with those studies [23,24,25,26,27,28,29,30,31,32,33]. Generally, all five methods displayed NVB deformations and lower (than MHP) quantitative amounts of stress (despite other visible differences), due to the extremely small tissular biomechanical deformations and displacements imposed by the anatomical micro-architecture. However, if the local anatomical topographic conditions would have allowed larger deformations and tissue displacements, these qualitative–quantitative differences would have been even larger, confirming the above-mentioned.
Rotation was reported to be the most stressful by all five methods during the periodontal breakdown, in agreement with our previous research [15,16,17,18,19,20,21,22], other research [23,24,25,26,27,28,29,30,31,32,33] and reports by Wu et al. [35,36,37]. Intrusion and extrusion were also seen as stressful movements (Table 2), in agreement with reports by Minch et al. [38] and Hofman et al. [39,40] (but for the intact periodontium).
The other three methods (Maximum and Minimum Principal and hydrostatic pressure), despite displaying quantitative data lower than MHP and confirming the safety of 3 N appliance during periodontal breakdown, displayed qualitative stress patterns with various inconsistencies (negative/positive compression-tension in dental pulp where this was not clinically present and bizarre variations in color-coded stress display from one movement to another and between bone loss levels) with no biomechanical scientific explanations. Since the dental pulp is held in an incompressible cavity made entirely of solid material, no deformations or displacements clinically occur, and thus the color-coded projections should display the same color of low/no stress intensity during the entire periodontal breakdown. However, the three methods displayed various stress intensities contradicting the mechanical physical principles. Thus, the only possible explanation is that the lower accuracy is due to design issues (brittle/liquid type of material).
The analysis of other numerical reports [1,2,3,4,35,36,37,38,39,40,49,50,51,52,53,54,55,56,65,66] confirmed the above-mentioned differences. Multiple intact periodontium PDL numerical studies have assessed a few orthodontic movements, one/two [39,40,56] anatomically simplified 3D models (i.e., upper first premolar [36,37,39,40,56]; canine [35,36,37]; first molar [65] and incisor [37,50,51,52,53,54]), intact periodontium [35,36,37,39,40,56] or variable levels of reduced periodontium [52,53,65] and employed failure criterions (hydrostatic pressure [35,36,37,39,40,56,66], Maximum and/or Minimum Principal [49,50,51,52,53] and VM [49,54,56]) without any reference to the type of material studied, homogeneity/non-homogeneity or MHP [1,49,50,51,52,53,54], employing non-linearity [49,50,51]/linearity [52,53,54].
Some reported results exceeded physiological MHP [49,50,51,52,53,66] when using 1–6 N loads [35,36,37,39,40,66], and with hydrostatic pressure [39,40,56] being the only method for the study of PDL (despite not performing any comparison/correlation with other methods), this contradicts clinical knowledge. Wu et al. reported as optimal a rotation of 2.1–2.9 N [35,36,37], contradicting Profit et al. [34] and clinical knowledge (optimal should be light to avoid any ischemic and necrotic risks).
These major issues were due to tissular models since PDL studies [1,35,36,37,38,39,40,49,50,51,52,53,54,56,65,66] did not reconstruct the NVB structure. Thus, despite reporting various amounts of stress, the accuracy is debatable (e.g., the PDL apical third amounts of stress exceeding MHP for light orthodontic forces [35,36,37,39,40] are clinically incorrect). More biomechanically bizarre is the qualitative stress display of extremely high PDL apical third stress associated with insignificant cervical third stress [35,36,37], which is biomechanically fundamentally incorrect [9,10,11,34,42,43,50,51]. From the mechanical physics point of view, correlated with tissular anatomical micro-architecture, the first major stress should be displayed in the cervical third, where the PDL surface is important, and the absorption-dissipation ability is high. Additionally, in the apical third from the same mechanical point of view, the stress should be at a minimum since most of the stresses are absorbed and dissipated in the middle and cervical thirds. Moreover, the vascularization is much better in the apical third than in cervical and middle thirds, while the natural tissular adaptation ability is higher. Other research (with no NVB tissular reconstruction) reported the opposite, with high cervical stress and almost no apical third stress [1]. The above-mentioned reports [1,35,36,37,38,39,40,49,50,51,52,53,54,55,56,65,66] suffer from the same methodological issues: the employed method lacked a material-based method selection and no correlations/comparations with other methods.
The numerical method limits are related to the employed proper failure criteria and the anatomical accuracy of the 3D model, as well as multiple simulations on the same model for indirect validation. Due to the importance of pulp and NVB functionality in orthodontic treatment [5,6,7,8,9,10,11,12,13,14], and the need for an individualized biomechanical behavioral study, numerical studies are the only option. Numerical studies are not able to completely mimic the clinical conditions (usually there are no pure movements and only associations). Thus, clinically, the quantitative amounts of stress could be a little bit lower than in this study, but with no impact over the report’s accuracy.
The proper method/failure criteria selection is related to the brittleness–ductileness of the analyzed tissues [64]. There are many new numerical studies employing brittle-like material methods for the dental tissues of ductile resemblance (i.e., the same as the above-mentioned flows) reporting results with accuracy issues [1,2,3,4]. In the engineering field, there are significant mechanical behavioral differences between brittle and ductile materials. In dentistry, due to limited deformations and displacements, these differences are seen only if the behavioral results are put side by side. Nevertheless, the differences are significant from the mechanical point of view when describing the tissular biomechanics. The anatomical accuracy of the models is directly related to the number of elements and nodes and the small global element size (i.e., a variable PDL thickness of 0.15–0.225 mm vs. a constant thickness of 0.25 mm [1]; 6 million tetrahedral elements and 1 million nodes, global element size of 0.08–0.116 mm, 40–12,731 times more elements, 4.4–1463 times more nodes than previous studies [9,10,11,39,40,42,43,49,50,51,52,53,54,55,56]). The CBCT image accuracy (small voxel size) also influences the anatomical accuracy (a voxel size of 0.075 mm herein vs. 0.3 mm [1]). A larger sample size also enhances the result accuracy (nine in this study vs. one in current studies [1,50,51,54,57,58,59,60,61,62,63,64]). Multiple simulations also increase the reports’ accuracy (1800 in this study vs. a few in most of the other studies [35,36,37,38,39,40,49,50,51,52,53,54,56,57,58,59,60,61,62,63,64,65,66]). Since there are so many numerical studies with accuracy issues (due to a methodological lack of understanding) [5,7,42,64], there is a need for numerical studies to address the above-mentioned issues. The limitations of the numerical studies in terms of clinical practice are mostly related to the fact that they cannot accurately reproduce the clinical situations, meaning that the quantitative results herein provided are lower than clinical ones. Qualitatively, since the associations of movements are often used (instead of pure movements as in this study), this means that the stress display areas are in clinical reality merged. Nevertheless, the results from this study do not lose their accuracy, since each time the results are clinically lower, while this study provided the maximum amount of stress and tissular deformations induced by the movements.

5. Conclusions

  • During the periodontal breakdown, all five methods displayed, for all five movements, quantitative amounts of stress lower than MHP, suggesting that 3 N are not inducing any local tissular ischemic risks for the healthy intact NVB and pulp.
  • During the periodontal breakdown, among the five methods, Tresca and Von Mises were more accurate than the other three, with Tresca being the most accurate of all.
  • All five methods displayed rotation as the most stressful movement for 1–8 mm of bone loss, while translation was the least stressful.
  • NVB was more exposed to ischemic risks than dental pulp during the periodontal breakdown due to constant tissular deformations during the 1–8 mm bone loss simulation, visible for all five movements.
  • Only VM and Tresca showed the translation movement as more prone to exposing dental pulp (both coronal and radicular) to more ischemic risks than the other five movements during the periodontal breakdown simulation.
  • All five methods showed intrusion and extrusion as more prone to exposing the NVB to higher ischemic risks than the other movements during the periodontal breakdown simulation.

6. Practical Implications

This research is the first to perform a comparison of five methods during horizontal periodontal breakdown and to highlight the methodological issues found in the current research flow with an impact on the current reports’ accuracy. It is also the first to study the biomechanical behavior of pulp and the NVB under 3 N of orthodontic force, during five movements in various levels of periodontal loss up to 8 mm loss. It is of clinically significant importance to know that 3 N of orthodontic force has almost no ischemic risks over the pulp and the NVB in healthy intact tissues, but this should be carefully used if previous injury/trauma are present in the conditions of bone loss presence. It also provides a complete set of biomechanical data extremely useful for the orthodontic planning phase. This study provides a clear comparative analysis, which is extremely okuseful when studying dental tissues. It also provides knowledge to better understand the nature of numerical studies as well as the requirements for accurate results.

Author Contributions

Conceptualization: R.-A.M.; Methodology: R.-A.M.; Software: R.-A.M.; Validation: R.-A.M. and C.D.O.; Formal analysis: R.-A.M.; Investigation: R.-A.M.; Resources: R.-A.M.; Data Curation: R.-A.M.; Writing-original draft preparation: R.-A.M.; Writing-review and editing: R.-A.M., A.G.D. and C.D.O.; Visualization- Supervision- Project administration: R.-A.M., A.G.D. and C.D.O.; Funding acquisition: R.-A.M., A.G.D. and C.D.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. The authors were the funders of this research project.

Institutional Review Board Statement

The research protocol has been approved by the Ethical Committee of the University of Medicine (158/2.04.2018).

Informed Consent Statement

Informed written consent was obtained from all subjects involved in the study.

Data Availability Statement

The original contributions presented in this study are included in the article. Further inquiries can be directed to the corresponding author(s).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

  1. Xing, J.; Zhang, G.; Sun, M.; Pan, H.; Zhang, C.; Liu, Y.; Li, K.; He, Z.; Zhang, K.; Wang, J.; et al. Clinical insights into tooth extraction via torsion method: A biomechanical analysis of the tooth-periodontal ligament complex. Front. Bioeng. Biotechnol. 2024, 12, 1479751. [Google Scholar] [CrossRef] [PubMed]
  2. Zhu, Y.; Zhu, J.; Yin, D.; Liu, Y. Improved stomatognathic model for highly realistic finite element analysis of temporomandibular joint biomechanics. J. Mech. Behav. Biomed. Mater. 2024, 160, 106780. [Google Scholar] [CrossRef] [PubMed]
  3. Zhu, Y.; Zheng, F.; Gong, Y.; Zhu, J.; Yin, D.; Liu, Y. Effect of occlusal contact on TMJ loading during occlusion: An in silico study. Comput. Biol. Med. 2024, 178, 108725. [Google Scholar] [CrossRef] [PubMed]
  4. Zheng, F.; Gong, Y.; Zhu, Y.; Yin, D.; Liu, Y. Three-dimensional theoretical model for effectively describing the effect of craniomaxillofacial structural factors on loading situation in the temporomandibular joint. J. Mech. Behav. Biomed. Mater. 2024, 151, 106371. [Google Scholar] [CrossRef] [PubMed]
  5. Yamaguchi, M.; Fukasawa, S. Is Inflammation a Friend or Foe for Orthodontic Treatment?: Inflammation in Orthodontically Induced Inflammatory Root Resorption and Accelerating Tooth Movement. Int. J. Mol. Sci. 2021, 22, 2388. [Google Scholar] [CrossRef]
  6. Weissheimer, T.; Silva, E.; Pinto, K.P.; Só, G.B.; Rosa, R.A.; Só, M.V.R. Do orthodontic tooth movements induce pulp necrosis? A systematic review. Int. Endod. J. 2021, 54, 1246–1262. [Google Scholar] [CrossRef]
  7. Vitali, F.C.; Cardoso, I.V.; Mello, F.W.; Flores-Mir, C.; Andrada, A.C.; Dutra-Horstmann, K.L.; Duque, T.M. Effect of orthodontic force on dental pulp histomorphology and tissue factor expression. Angle Orthod. 2021, 91, 830–842. [Google Scholar] [CrossRef]
  8. Vermiglio, G.; Centofanti, A.; Matarese, G.; Militi, A.; Matarese, M.; Arco, A.; Nicita, F.; Cutroneo, G. Human Dental Pulp Tissue during Orthodontic Tooth Movement: An Immunofluorescence Study. J. Funct. Morphol. Kinesiol. 2020, 5, 65. [Google Scholar] [CrossRef]
  9. Bauss, O.; Rohling, J.; Meyer, K.; Kiliaridis, S. Pulp vitality in teeth suffering trauma during orthodontic therapy. Angle Orthod. 2009, 79, 166–171. [Google Scholar] [CrossRef]
  10. Bauss, O.; Rohling, J.; Rahman, A.; Kiliaridis, S. The effect of pulp obliteration on pulpal vitality of orthodontically intruded traumatized teeth. J. Endod. 2008, 34, 417–420. [Google Scholar] [CrossRef]
  11. Bauss, O.; Schäfer, W.; Sadat-Khonsari, R.; Knösel, M. Influence of orthodontic extrusion on pulpal vitality of traumatized maxillary incisors. J. Endod. 2010, 36, 203–207. [Google Scholar] [CrossRef] [PubMed]
  12. Patro, S.; Meto, A.; Mohanty, A.; Chopra, V.; Miglani, S.; Das, A.; Luke, A.M.; Hadi, D.A.; Meto, A.; Fiorillo, L.; et al. Diagnostic Accuracy of Pulp Vitality Tests and Pulp Sensibility Tests for Assessing Pulpal Health in Permanent Teeth: A Systematic Review and Meta-Analysis. Int J Env. Res Public Health 2022, 19, 9599. [Google Scholar] [CrossRef] [PubMed]
  13. Cărămizaru, M.; Pleşea, I.E.; Dragomir, L.P.; Popescu, M.R.; Uscatu, C.D.; Şerbănescu, M.S.; Alexandru, D.O.; Comănescu, T.M. Quantitative assessment of morphological changes of dental pulp components of teeth affected by occlusal trauma. Rom. J. Morphol. Embryol. = Rev. Roum. Morphol. Embryol. 2018, 59, 729–740. [Google Scholar]
  14. França, C.M.; Riggers, R.; Muschler, J.L.; Widbiller, M.; Lococo, P.M.; Diogenes, A.; Bertassoni, L.E. 3D-Imaging of Whole Neuronal and Vascular Networks of the Human Dental Pulp via CLARITY and Light Sheet Microscopy. Sci. Rep. 2019, 9, 10860. [Google Scholar] [CrossRef] [PubMed]
  15. Moga, R.A.; Olteanu, C.D.; Botez, M.D.; Buru, S.M. Assessment of the Orthodontic External Resorption in Periodontal Breakdown—A Finite Elements Analysis (Part I). Healthcare 2023, 11, 1447. [Google Scholar] [CrossRef]
  16. Moga, R.A.; Delean, A.G.; Buru, S.M.; Botez, M.D.; Olteanu, C.D. Orthodontic Internal Resorption Assessment in Periodontal Breakdown—A Finite Elements Analysis (Part II). Healthcare 2023, 11, 2622. [Google Scholar] [CrossRef]
  17. Moga, R.A.; Olteanu, C.D.; Botez, M.; Buru, S.M. Assessment of the Maximum Amount of Orthodontic Force for Dental Pulp and Apical Neuro-Vascular Bundle in Intact and Reduced Periodontium on Bicuspids (Part II). Int. J. Env. Res. Public Health 2023, 20, 1179. [Google Scholar] [CrossRef]
  18. Moga, R.A.; Buru, S.M.; Olteanu, C.D. Assessment of the Best FEA Failure Criteria (Part II): Investigation of the Biomechanical Behavior of Dental Pulp and Apical-Neuro-Vascular Bundle in Intact and Reduced Periodontium. Int. J. Env. Res. Public Health 2022, 19, 15635. [Google Scholar] [CrossRef]
  19. Moga, R.A.; Buru, S.M.; Chiorean, C.G. Overall stress in periodontal ligament under orthodontic movement during a periodontal breakdown. Am. J. Orthod. Dentofac. Orthop. Off. Publ. Am. Assoc. Orthod. Its Const. Soc. Am. Board Orthod. 2022, 161, e127–e135. [Google Scholar] [CrossRef]
  20. Moga, R.A.; Olteanu, C.D.; Buru, S.M.; Botez, M.D.; Delean, A.G. Cortical and Trabecular Bone Stress Assessment during Periodontal Breakdown—A Comparative Finite Element Analysis of Multiple Failure Criteria. Medicina 2023, 59, 1462. [Google Scholar] [CrossRef]
  21. Moga, R.A.; Olteanu, C.D.; Buru, S.M.; Botez, M.D.; Delean, A.G. Finite Elements Analysis of Biomechanical Behavior of the Bracket in a Gradual Horizontal Periodontal Breakdown—A Comparative Analysis of Multiple Failure Criteria. Appl. Sci. 2023, 13, 9480. [Google Scholar] [CrossRef]
  22. Moga, R.-A.; Olteanu, C.D.; Delean, A.G. Investigating the Ability of the Tooth and Surrounding Support Tissues to Absorb and Dissipate Orthodontic Loads during Periodontal Breakdown—Finite Elements Analysis. Appl. Sci. 2024, 14, 1041. [Google Scholar] [CrossRef]
  23. Cardenas-Duque, L.M.; Yoshida, M.; Goto, G. Pulpal response to different pulp capping methods after pulp exposure by air abrasion. J. Clin. Pediatr. Dent. 2002, 26, 269–273. [Google Scholar] [CrossRef]
  24. Murray, P.E.; Hafez, A.A.; Windsor, L.J.; Smith, A.J.; Cox, C.F. Comparison of pulp responses following restoration of exposed and non-exposed cavities. J. Dent. 2002, 30, 213–222. [Google Scholar] [CrossRef] [PubMed]
  25. Mjör, I.A. Pulp-dentin biology in restorative dentistry. Part 7: The exposed pulp. Quintessence Int. 2002, 33, 113–135. [Google Scholar]
  26. Murray, P.E.; Hafez, A.A.; Smith, A.J.; Cox, C.F. Identification of hierarchical factors to guide clinical decision making for successful long-term pulp capping. Quintessence Int. 2003, 34, 61–70. [Google Scholar]
  27. Kitasako, Y.; Murray, P.E.; Tagami, J.; Smith, A.J. Histomorphometric analysis of dentinal bridge formation and pulpal inflammation. Quintessence Int. 2002, 33, 600–608. [Google Scholar]
  28. Medina, V.O., 3rd; Shinkai, K.; Shirono, M.; Tanaka, N.; Katoh, Y. Histopathologic study on pulp response to single-bottle and self-etching adhesive systems. Oper. Dent. 2002, 27, 330–342. [Google Scholar]
  29. Suzuki, M.; Katsumi, A.; Watanabe, R.; Shirono, M.; Katoh, Y. Effects of an experimentally developed adhesive resin system and CO2 laser irradiation on direct pulp capping. Oper. Dent. 2005, 30, 702–718. [Google Scholar]
  30. Kitasako, Y.; Ikeda, M.; Tagami, J. Pulpal responses to bacterial contamination following dentin bridging beneath hard-setting calcium hydroxide and self-etching adhesive resin system. Dent. Traumatol. Off. Publ. Int. Assoc. Dent. Traumatol. 2008, 24, 201–206. [Google Scholar] [CrossRef]
  31. Schuurs, A.H.; Gruythuysen, R.J.; Wesselink, P.R. Pulp capping with adhesive resin-based composite vs. calcium hydroxide: A review. Endod. Dent. Traumatol. 2000, 16, 240–250. [Google Scholar] [CrossRef] [PubMed]
  32. Cox, C.F.; Hafez, A.A. Biocomposition and reaction of pulp tissues to restorative treatments. Dent. Clin. N. Am. 2001, 45, 31–48. [Google Scholar] [CrossRef]
  33. Farughi, A.; Rouhani, A.; Shahmohammadi, R.; Jafarzadeh, H. Clinical comparison of sensitivity and specificity between sensibility and vitality tests in determining the pulp vitality of mandibular premolars. Aust. Endod. J. J. Aust. Soc. Endodontol. Inc. 2021, 47, 474–479. [Google Scholar] [CrossRef]
  34. Proffit, W.R.; Fields, H.W., Jr.; Sarver, D.M. Contemporary Orthodontics, 5th ed.; Elsevier: St. Louis, MO, USA, 2012. [Google Scholar]
  35. Wu, J.L.; Liu, Y.F.; Peng, W.; Dong, H.Y.; Zhang, J.X. A biomechanical case study on the optimal orthodontic force on the maxillary canine tooth based on finite element analysis. J. Zhejiang Univ. Sci. B 2018, 7, 535–546. [Google Scholar] [CrossRef]
  36. Wu, J.; Liu, Y.; Wang, D.; Zhang, J.; Dong, X.; Jiang, X.; Xu, X. Investigation of effective intrusion and extrusion force for maxillary canine using finite element analysis. Comput. Methods Biomech. Biomed. Eng. 2019, 22, 1294–1302. [Google Scholar] [CrossRef] [PubMed]
  37. Wu, J.; Liu, Y.; Li, B.; Wang, D.; Dong, X.; Sun, Q.; Chen, G. Numerical simulation of optimal range of rotational moment for the mandibular lateral incisor, canine and first premolar based on biomechanical responses of periodontal ligaments: A case study. Clin. Oral Investig. 2021, 25, 1569–1577. [Google Scholar] [CrossRef] [PubMed]
  38. Minch, L.E.; Sarul, M.; Nowak, R.; Kawala, B.; Antoszewska-Smith, J. Orthodontic intrusion of periodontally-compromised maxillary incisors: 3-dimensional finite element method analysis. Adv. Clin. Exp. Med. Off. Organ Wroc. Med. Univ. 2017, 26, 829–833. [Google Scholar] [CrossRef]
  39. Hohmann, A.; Wolfram, U.; Geiger, M.; Boryor, A.; Kober, C.; Sander, C.; Sander, F.G. Correspondences of hydrostatic pressure in periodontal ligament with regions of root resorption: A clinical and a finite element study of the same human teeth. Comput. Methods Programs Biomed. 2009, 93, 155–161. [Google Scholar] [CrossRef]
  40. Hohmann, A.; Wolfram, U.; Geiger, M.; Boryor, A.; Sander, C.; Faltin, R.; Faltin, K.; Sander, F.G. Periodontal ligament hydrostatic pressure with areas of root resorption after application of a continuous torque moment. Angle Orthod. 2007, 77, 653–659. [Google Scholar] [CrossRef]
  41. Rusu Olaru, A.; Popescu, M.R.; Pleşea, I.E.; Şerbănescu, M.S.; Pleşea, R.M.; Cojocaru, M.O.; Coculescu, E.C. Abrasion and dental pulp morphological changes in occlusal dysfunction. Rom. J. Morphol. Embryol. = Rev. Roum. Morphol. Embryol. 2024, 65, 279–290. [Google Scholar] [CrossRef]
  42. Javed, F.; Al-Kheraif, A.A.; Romanos, E.B.; Romanos, G.E. Influence of orthodontic forces on human dental pulp: A systematic review. Arch. Oral Biol. 2015, 60, 347–356. [Google Scholar] [CrossRef] [PubMed]
  43. Bauss, O.; Röhling, J.; Sadat-Khonsari, R.; Kiliaridis, S. Influence of orthodontic intrusion on pulpal vitality of previously traumatized maxillary permanent incisors. Am. J. Orthod. Dentofac. Orthop. Off. Publ. Am. Assoc. Orthod. Its Const. Soc. Am. Board Orthod. 2008, 134, 12–17. [Google Scholar] [CrossRef] [PubMed]
  44. Strobl, H.; Haas, M.; Norer, B.; Gerhard, S.; Emshoff, R. Evaluation of pulpal blood flow after tooth splinting of luxated permanent maxillary incisors. Dent. Traumatol. Off. Publ. Int. Assoc. Dent. Traumatol. 2004, 20, 36–41. [Google Scholar] [CrossRef] [PubMed]
  45. Emshoff, R.; Emshoff, I.; Moschen, I.; Strobl, H. Diagnostic characteristics of pulpal blood flow levels associated with adverse outcomes of luxated permanent maxillary incisors. Dent. Traumatol. Off. Publ. Int. Assoc. Dent. Traumatol. 2004, 20, 270–275. [Google Scholar] [CrossRef] [PubMed]
  46. Chen, E.; Abbott, P.V. Dental pulp testing: A review. Int. J. Dent. 2009, 2009, 365785. [Google Scholar] [CrossRef]
  47. Balevi, B. Cold pulp testing is the simplest and most accurate of all dental pulp sensibility tests. Evid.-Based Dent. 2019, 20, 22–23. [Google Scholar] [CrossRef]
  48. Mainkar, A.; Kim, S.G. Diagnostic Accuracy of 5 Dental Pulp Tests: A Systematic Review and Meta-analysis. J. Endod. 2018, 44, 694–702. [Google Scholar] [CrossRef]
  49. Toms, S.R.; Eberhardt, A.W. A nonlinear finite element analysis of the periodontal ligament under orthodontic tooth loading. Am. J. Orthod. Dentofac. Orthop. Off. Publ. Am. Assoc. Orthod. Its Const. Soc. Am. Board Orthod. 2003, 123, 657–665. [Google Scholar] [CrossRef]
  50. Hemanth, M.; Deoli, S.; Raghuveer, H.P.; Rani, M.S.; Hegde, C.; Vedavathi, B. Stress Induced in the Periodontal Ligament under Orthodontic Loading (Part I): A Finite Element Method Study Using Linear Analysis. J. Int. Oral Health JIOH 2015, 7, 129–133. [Google Scholar]
  51. Hemanth, M.; Deoli, S.; Raghuveer, H.P.; Rani, M.S.; Hegde, C.; Vedavathi, B. Stress Induced in Periodontal Ligament under Orthodontic Loading (Part II): A Comparison of Linear Versus Non-Linear Fem Study. J. Int. Oral Health JIOH 2015, 7, 114–118. [Google Scholar]
  52. Geramy, A. Initial stress produced in the periodontal membrane by orthodontic loads in the presence of varying loss of alveolar bone: A three-dimensional finite element analysis. Eur. J. Orthod. 2002, 24, 21–33. [Google Scholar] [CrossRef]
  53. Geramy, A.; Faghihi, S. Secondary trauma from occlusion: Three-dimensional analysis using the finite element method. Quintessence Int. 2004, 35, 835–843. [Google Scholar] [PubMed]
  54. Shaw, A.M.; Sameshima, G.T.; Vu, H.V. Mechanical stress generated by orthodontic forces on apical root cementum: A finite element model. Orthod. Craniofacial Res. 2004, 7, 98–107. [Google Scholar] [CrossRef]
  55. Gupta, M.; Madhok, K.; Kulshrestha, R.; Chain, S.; Kaur, H.; Yadav, A. Determination of stress distribution on periodontal ligament and alveolar bone by various tooth movements—A 3D FEM study. J. Oral Biol. Craniofacial Res. 2020, 10, 758–763. [Google Scholar] [CrossRef]
  56. Roscoe, M.G.; Cattaneo, P.M.; Dalstra, M.; Ugarte, O.M.; Meira, J.B.C. Orthodontically induced root resorption: A critical analysis of finite element studies’ input and output. Am. J. Orthod. Dentofac. Orthop. Off. Publ. Am. Assoc. Orthod. Its Const. Soc. Am. Board Orthod. 2021, 159, 779–789. [Google Scholar] [CrossRef]
  57. Prados-Privado, M.; Martínez-Martínez, C.; Gehrke, S.A.; Prados-Frutos, J.C. Influence of Bone Definition and Finite Element Parameters in Bone and Dental Implants Stress: A Literature Review. Biology 2020, 9, 224. [Google Scholar] [CrossRef]
  58. Yamanishi, Y.; Yamaguchi, S.; Imazato, S.; Nakano, T.; Yatani, H. Effects of the implant design on peri-implant bone stress and abutment micromovement: Three-dimensional finite element analysis of original computer-aided design models. J. Periodontol. 2014, 85, e333–e338. [Google Scholar] [CrossRef] [PubMed]
  59. Pérez-Pevida, E.; Brizuela-Velasco, A.; Chávarri-Prado, D.; Jiménez-Garrudo, A.; Sánchez-Lasheras, F.; Solaberrieta-Méndez, E.; Diéguez-Pereira, M.; Fernández-González, F.J.; Dehesa-Ibarra, B.; Monticelli, F. Biomechanical Consequences of the Elastic Properties of Dental Implant Alloys on the Supporting Bone: Finite Element Analysis. BioMed Res. Int. 2016, 2016, 1850401. [Google Scholar] [CrossRef] [PubMed]
  60. Tawara, D.; Nagura, K. Predicting changes in mechanical properties of trabecular bone by adaptive remodeling. Comput. Methods Biomech. Biomed. Eng. 2017, 20, 415–425. [Google Scholar] [CrossRef]
  61. Aunmeungtong, W.; Khongkhunthian, P.; Rungsiyakull, P. Stress and strain distribution in three different mini dental implant designs using in implant retained overdenture: A finite element analysis study. ORAL Implantol. 2016, 9, 202–212. [Google Scholar]
  62. Merdji, A.; Bachir Bouiadjra, B.; Achour, T.; Serier, B.; Ould Chikh, B.; Feng, Z.O. Stress analysis in dental prosthesis. Comput. Mater. Sci. 2010, 49, 126–133. [Google Scholar] [CrossRef]
  63. Field, C.; Ichim, I.; Swain, M.V.; Chan, E.; Darendeliler, M.A.; Li, W.; Li, Q. Mechanical responses to orthodontic loading: A 3-dimensional finite element multi-tooth model. Am. J. Orthod. Dentofac. Orthop. Off. Publ. Am. Assoc. Orthod. Its Const. Soc. Am. Board Orthod. 2009, 135, 174–181. [Google Scholar] [CrossRef] [PubMed]
  64. Perez-Gonzalez, A.; Iserte-Vilar, J.L.; Gonzalez-Lluch, C. Interpreting finite element results for brittle materials in endodontic restorations. Biomed. Eng. Online 2011, 10, 44. [Google Scholar] [CrossRef] [PubMed]
  65. Hussein Mahmood Ghuloom, K.; Mascarenhas, R.; Parveen, S.; Husain, A. Finite element analysis of orthodontically induced stress in the periodontal ligament of the maxillary first molar with simulated bone loss. J. Comput. Methods Sci. Eng. 2017, 17, 243–252. [Google Scholar] [CrossRef]
  66. Zhong, J.; Chen, J.; Weinkamer, R.; Darendeliler, M.A.; Swain, M.V.; Sue, A.; Zheng, K.; Li, Q. In vivo effects of different orthodontic loading on root resorption and correlation with mechanobiological stimulus in periodontal ligament. J. R. Soc. Interface 2019, 16, 20190108. [Google Scholar] [CrossRef]
Figure 1. Mesh model (one of the nine 3D models): (A) 2nd lower right premolar model with 1 mm loss periodontium, (B) 4 mm loss model, (C) 8 mm loss model, (D) 1 mm loss models with applied extrusion and restraining movements base vectors, (E) second lower premolar with the base of the bracket and the NVB and (F) dental pulp and NVB. Applied vectors: (G) intrusion, (H) rotation, (I) tipping, (J) translation, (K) dental pulp and NVB mesh with elements warnings and (L) second lower premolar with mesh elements warnings.
Figure 1. Mesh model (one of the nine 3D models): (A) 2nd lower right premolar model with 1 mm loss periodontium, (B) 4 mm loss model, (C) 8 mm loss model, (D) 1 mm loss models with applied extrusion and restraining movements base vectors, (E) second lower premolar with the base of the bracket and the NVB and (F) dental pulp and NVB. Applied vectors: (G) intrusion, (H) rotation, (I) tipping, (J) translation, (K) dental pulp and NVB mesh with elements warnings and (L) second lower premolar with mesh elements warnings.
Jcm 13 07062 g001
Figure 2. Comparative stress distribution displayed by the five methods for 3 N applied force (of one of the nine 3D models), in 1 mm periodontal breakdown and five movements (extrusion, intrusion, rotation, tipping and translation): (A) hydrostatic pressure, (B) Minimum Principal, (C) Maximum Principal, (D) Tresca and (E) Von Mises.
Figure 2. Comparative stress distribution displayed by the five methods for 3 N applied force (of one of the nine 3D models), in 1 mm periodontal breakdown and five movements (extrusion, intrusion, rotation, tipping and translation): (A) hydrostatic pressure, (B) Minimum Principal, (C) Maximum Principal, (D) Tresca and (E) Von Mises.
Jcm 13 07062 g002
Figure 3. Comparative stress distribution displayed by the five methods for 3 N applied force (of one of the nine 3D models), in 4 mm periodontal breakdown and five movements (extrusion, intrusion, rotation, tipping and translation): (A) hydrostatic pressure, (B) Minimum Principal, (C) Maximum Principal, (D) Tresca, (E) Von Mises.
Figure 3. Comparative stress distribution displayed by the five methods for 3 N applied force (of one of the nine 3D models), in 4 mm periodontal breakdown and five movements (extrusion, intrusion, rotation, tipping and translation): (A) hydrostatic pressure, (B) Minimum Principal, (C) Maximum Principal, (D) Tresca, (E) Von Mises.
Jcm 13 07062 g003
Figure 4. Comparative stress distribution displayed by the five methods for 3 N applied force (of one of the nine 3D models), in 8 mm periodontal breakdown and five movements (extrusion, intrusion, rotation, tipping and translation): (A) hydrostatic pressure, (B) Minimum Principal, (C) Maximum Principal, (D) Tresca and (E) Von Mises.
Figure 4. Comparative stress distribution displayed by the five methods for 3 N applied force (of one of the nine 3D models), in 8 mm periodontal breakdown and five movements (extrusion, intrusion, rotation, tipping and translation): (A) hydrostatic pressure, (B) Minimum Principal, (C) Maximum Principal, (D) Tresca and (E) Von Mises.
Jcm 13 07062 g004
Table 1. Physical properties of materials.
Table 1. Physical properties of materials.
MaterialsYoung’s Modulus, E (GPa)Poisson Ratio, ʋRefs.
Enamel800.33[15,16,17,18,19,20,21,22]
Dentin/Cementum 18.60.31[15,16,17,18,19,20,21,22]
Pulp and NVB0.00210.45[15,16,17,18,19,20,21,22]
PDL0.06670.49[15,16,17,18,19,20,21,22]
Cortical bone14.50.323[15,16,17,18,19,20,21,22]
Trabecular bone1.370.3[15,16,17,18,19,20,21,22]
Stainless streel bracket (Cr-Co)2180.33[15,16,17,18,19,20,21,22]
Table 2. Maximum stress average values (KPa) produced by 3 N in the NVB and coronal pulp during periodontal breakdown and five orthodontic movements.
Table 2. Maximum stress average values (KPa) produced by 3 N in the NVB and coronal pulp during periodontal breakdown and five orthodontic movements.
Resorption (mm) 12345678
IntrusionTrescaNVB0.840.991.131.281.501.731.952.17
c0.080.090.100.110.130.150.170.18
VMNVB0.710.840.981.121.311.551.641.88
c0.070.080.090.100.110.130.150.16
PressureNVB2.963.273.593.914.515.105.696.28
c0.270.310.350.400.580.770.961.14
S1NVB−2.56−2.90−3.24−3.57−3.99−4.40−4.80−5.22
c0.270.290.310.330.390.440.490.55
S3NVB−2.99−3.22−3.45−3.68−4.17−4.67−5.16−5.66
c0.200.200.210.210.270.320.380.44
ExtrusionTrescaNVB0.840.991.131.281.501.731.952.17
c0.080.090.100.110.130.150.170.18
VMNVB0.710.840.981.121.311.551.641.88
c0.070.080.090.100.110.130.150.16
PressureNVB−2.96−3.27−3.59−3.91−4.51−5.10−5.69−6.28
c−0.27−0.31−0.35−0.40−0.58−0.77−0.96−1.14
S1NVB2.983.443.904.354.695.035.365.70
c−0.55−0.61−0.670.720.740.770.810.83
S3NVB3.053.293.523.754.294.825.365.89
c−0.26−0.28−0.30−0.32−0.36−0.41−0.44−0.48
TranslationTrescaNVB0.640.720.800.880.981.081.191.29
c0.110.120.140.150.160.170.170.18
VMNVB0.550.620.690.760.860.961.071.17
c0.130.140.140.150.150.150.150.16
PressureNVB2.432.542.652.763.233.714.194.66
c−0.35−0.36−0.370.380.450.510.580.64
S1NVB−2.11−2.24−2.39−2.51−2.94−3.36−3.79−4.22
c0.280.320.350.380.420.470.510.55
S3NVB−2.57−2.71−2.85−2.99−3.52−4.05−4.58−5.10
c−0.10−0.11−0.110.110.130.151.760.20
RotationTrescaNVB0.901.211.512.072.322.572.823.07
3 N c0.160.170.180.190.210.240.260.28
VMNVB0.971.191.411.631.741.851.952.06
c0.140.150.160.170.190.220.240.27
PressureNVB3.984.535.085.636.337.037.738.42
c1.691.761.831.902.132.362.592.82
S1NVB−3.98−4.67−5.33−6.00−6.72−7.44−8.16−8.87
c0.881.201.511.832.082.322.572.81
S3NVB−4.28−4.95−5.61−6.28−7.05−7.83−8.60−9.37
c0.470.560.65−0.74−0.82−0.91−0.99−1.08
TippingTrescaNVB0.911.091.261.441.641.842.042.24
c0.130.140.140.140.160.170.180.19
VMNVB0.790.941.101.251.431.621.801.98
c0.110.110.110.110.130.140.160.17
PressureNVB3.093.554.014.474.945.425.896.36
c1.171.211.251.291.331.371.401.44
S1NVB−2.80−3.05−3.29−4.18−4.57−4.95−5.34−5.72
c−0.69−0.88−1.071.261.541.822.092.37
S3NVB−3.19−3.56−3.92−4.28−4.87−5.46−6.05−6.64
c−0.30−0.35−0.40−0.440.450.470.470.48
NVB—neuro-vascular bundle, c—coronal pulp.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Moga, R.-A.; Olteanu, C.D.; Delean, A.G. Periodontal Breakdown, Orthodontic Movements and Pulpal Ischemia Correlations—A Comparison Between Five Study Methods. J. Clin. Med. 2024, 13, 7062. https://doi.org/10.3390/jcm13237062

AMA Style

Moga R-A, Olteanu CD, Delean AG. Periodontal Breakdown, Orthodontic Movements and Pulpal Ischemia Correlations—A Comparison Between Five Study Methods. Journal of Clinical Medicine. 2024; 13(23):7062. https://doi.org/10.3390/jcm13237062

Chicago/Turabian Style

Moga, Radu-Andrei, Cristian Doru Olteanu, and Ada Gabriela Delean. 2024. "Periodontal Breakdown, Orthodontic Movements and Pulpal Ischemia Correlations—A Comparison Between Five Study Methods" Journal of Clinical Medicine 13, no. 23: 7062. https://doi.org/10.3390/jcm13237062

APA Style

Moga, R.-A., Olteanu, C. D., & Delean, A. G. (2024). Periodontal Breakdown, Orthodontic Movements and Pulpal Ischemia Correlations—A Comparison Between Five Study Methods. Journal of Clinical Medicine, 13(23), 7062. https://doi.org/10.3390/jcm13237062

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop