Constitutive modeling of the non-linear visco-elasticity of the periodontal ligament

Functional non-uniformity of periodontal ligaments tunes mechanobiological stimuli across soft- and hard-tissue interfaces

The bone-periodontal ligament-tooth (BPT) complex is a unique mechanosensing soft-/hard-tissue interface, which governs the most rapid bony homeostasis in the body responding to external loadings. While the correlation between such loading and alveolar bone remodelling has been widely recognised, it has remained challenging to investigate the transmitted mechanobiological stimuli across such embedded soft-/hard-tissue interfaces of the BPT complex. Here, we propose a framework combining three distinct bioengineering techniques (i, ii, and iii below) to elucidate the innate functional non-uniformity of the PDL in tuning mechanical stimuli to the surrounding alveolar bone. The biphasic PDL mechanical properties measured via nanoindentation, namely the elastic moduli of fibres and ground substance at the sub-tissue level (i), were used as the input parameters in an image-based constitutive modelling framework for finite element simulation (ii). In tandem with U-net deep learning, the Gaussian mixture method enabled the comparison of 5195 possible pseudo-microstructures versus the innate non-uniformity of the PDL (iii). We found that the balance between hydrostatic pressure in PDL and the strain energy in the alveolar bone was maintained within a specific physiological range. The innate PDL microstructure ensures the transduction of favourable mechanobiological stimuli, thereby governing alveolar bone homeostasis. Our outcomes expand current knowledge of the PDL's mechanobiological roles and the proposed framework can be adopted to a broad range of similar soft-/hard- tissue interfaces, which may impact future tissue engineering, regenerative medicine, and evaluating therapeutic strategies. STATEMENT OF SIGNIFICANCE: A combination of cutting-edge technologies, including dynamic nanomechanical testing, high-resolution image-based modelling and machine learning facilitated computing, was used to elucidate the association between the microstructural non-uniformity and biomechanical competence of periodontal ligaments (PDLs). The innate PDL fibre network regulates mechanobiological stimuli, which govern alveolar bone remodelling, in different tissues across the bone-PDL-tooth (BPT) interfaces. These mechanobiological stimuli within the BPT are tuned within a physiological range by the non-uniform microstructure of PDLs, ensuring functional tissue homeostasis. The proposed framework in this study is also applicable for investigating the structure-function relationship in broader types of fibrous soft-/hard- tissue interfaces.

Construction of Human Periodontal Ligament Constitutive Model Based on Collagen Fiber Content

Periodontal ligament (PDL) is mainly composed of collagen fiber bundles, and the content of collagen fiber is an important factor affecting the mechanical properties of PDL. Based on this, the purpose of this study is to explore the effect of the PDL collagen fiber content on its viscoelastic mechanical behavior. Transverse and longitudinal samples of different regions of PDL were obtained from the human maxilla. The fiber content at different regions of human PDL was quantitatively measured using image processing software, and a new viscoelastic constitutive model was constructed based on the fiber content. The nano-indentation experiment was carried out with a loading rate of 0.5 mN·s-1, a peak load of 3 mN, and a holding time of 200 s, and the model parameters were obtained through the experiment data. The results showed that with the increase of fiber content, the deformation resistance of PDL also increased, and compared with the neck and middle region, the compressive strain in the apical region of PDL was the largest. The range of reduced elastic modulus of human PDL was calculated to be 0.39~5.08 MPa. The results of the experimental data and the viscoelastic constitutive model fit well, indicating that the model can well describe the viscoelastic behavior of human PDL.

Modelling and evaluating periodontal ligament mechanical behaviour and properties: A scoping review of current approaches and limitations

This scoping review is intended to synthesize the techniques proposed to model the tooth-periodontal ligament-bone complex (TPBC), while also evaluating the suggested periodontal ligament (PDL) material properties. It is concentrated on the recent advancements on the PDL and TPBC models, while identifying the advantages and limitations of the proposed approaches. Systematic searches were conducted up to December 2020 for articles that proposed PDL models to assess orthodontic tooth movement in Compendex, Web of Science, EMBASE, MEDLINE, PubMed, ScienceDirect, Google Scholar and Scopus databases. Although there have been many studies focused on the evaluation of PDL material properties through numerous modelling approaches, only a handful of approaches have been identified to investigate the interface properties of the PDL as a complete dynamical system (TPBC models). Past reviews on the analytical and experimental determination of the PDL properties already show a concerning range in reported output values-some nearly six orders of magnitude in difference-that strongly suggested the need for further investigation. Surprisingly, it has not yet been possible to determine a narrower range of values for the PDL material properties. Moreover, very few scientific approaches address the TPBC as an integrated complex system model. In consequence, current methods for capturing the PDL material behaviour in a clinical setting are limited and inconclusive. This synthesis encourages more systematic, pragmatic and phenomenological research in this area.

Tensile creep mechanical behavior of periodontal ligament: A hyper-viscoelastic constitutive model

OBJECTIVE

In orthodontic treatment, the biomechanical response of periodontal ligament (PDL) induces tooth movement. Coupling modeling of PDL can effectively reflect its biomechanical response. The nonlinear creep mechanical behavior of PDL was studied by uniaxial tensile creep test and a new hyper-viscoelastic constitutive model. Two coupling modeling methods with limitations were excluded.

METHODS

PDL specimens were prepared from the central incisors of pig mandible. The theoretical step function was replaced by static loading with a total loading time of 1 s. The creep loading with the constant stresses of 0.05, 0.1, and 0.15 MPa was selected and kept unchanged for 1000 s. The instantaneous hyperelastic mechanical behavior and time-dependent nonlinear viscoelastic mechanical behavior of PDL were characterized by coupled instantaneous third-order Ogden hyperelastic and time-dependent nonlinear creep models.

RESULTS

The results showed that the instantaneous elastic curve of PDL increases in the form of hyperelastic index. The creep strain and creep compliance curves increase rapidly before 200s, and then increase slowly in steady state. The creep strain increased with an increase in the constant stress; conversely, the creep compliance decreased with an increase in the constant stress. The results showed that the experimental data were highly consistent with the hyper-viscoelastic constitutive model (R²>0.97).

SIGNIFICANCE

We normalize the framework of hyper-viscoelastic coupling modeling (Instantaneous hyperelastic model + time-dependent nonlinear viscoelastic model). Which can be extended to other nonlinear viscoelastic biomaterials.

OpenMandible: An open-source framework for highly realistic numerical modelling of lower mandible physiology

OBJECTIVE

Computer modeling of lower mandible physiology remains challenging because prescribing realistic material characteristics and boundary conditions from medical scans requires advanced equipment and skill sets. The objective of this study is to provide a framework that could reduce simplifications made and inconsistency (in terms of geometry, materials, and boundary conditions) among further studies on the topic.

METHODS

The OpenMandible framework offers: 1) the first publicly available multiscale model of the mandible developed by combining cone beam computerized tomography (CBCT) and μCT imaging modalities, and 2) a C++ software tool for the generation of simulation-ready models (tet4 and hex8 elements). In addition to the application of conventional (Neumann and Dirichlet) boundary conditions, OpenMandible introduces a novel geodesic wave propagation - based approach for incorporating orthotropic micromechanical characteristics of cortical bone, and a unique algorithm for modeling muscles as uniformly directed vectors. The base intact model includes the mandible (spongy and compact bone), 14 teeth (comprising dentin, enamel, periodontal ligament, and pulp), simplified temporomandibular joints, and masticatory muscles (masseter, temporalis, medial, and lateral pterygoid).

RESULTS

The complete source code, executables, showcases, and sample data are freely available on the public repository: https://github.com/ArsoVukicevic/OpenMandible. It has been demonstrated that by slightly editing the baseline model, one can study different "virtual" treatments or diseases, including tooth restoration, placement of implants, mandible bone degradation, and others.

SIGNIFICANCE

OpenMandible eases the community to undertake a broad range of studies on the topic, while increasing their consistency and reproducibility. At the same time, the needs for dedicated equipment and skills for developing realistic simulation models are significantly reduced.

Mechanical effects of distributed fibre orientation in the periodontal ligament of an idealised geometry

In this study, we computationally assess the effects of the distributed fibre orientation in the periodontal ligament (PDL) on mechanical responses of the tooth-PDL complex. An idealised axial-symmetric geometry of a tooth-PDL complex was constructed. The fibre orientation in the PDL was modelled as a trigonometric function based on anatomical knowledge, and the PDL was modelled as a transversely isotropic hyperelastic material dependent on fibre orientations. Parametric studies of the fibre orientation on the mechanical responses of the tooth-PDL complex were conducted. Obtained results showed that the anatomically consistent fibre orientation functions as a supporting structure against not only vertical but also horizontal loads.

In silico study of cuspid' periodontal ligament damage under parafunctional and traumatic conditions of whole-mouth occlusions. A patient-specific evaluation

BACKGROUND AND OBJECTIVE

Although traumatic loading has been associated with periodontal ligament (PDL) damage and therefore with several oral disorders, the damage phenomena and the traumatic loads involved are still unclear. The complex composition and extremely thin size of the PDL make experimentation difficult, requiring computational studies that consider the macroscopic loading conditions, the microscopic composition and fine detailed geometry of the tissue. In this study, a new methodology to analyse the damage phenomena in the collagen network and the extracellular matrix of the PDL caused by parafunctional and traumatic occlusal forces was proposed.

METHODS

The entire human mandible and a portion thereof containing a full cuspid tooth were separately modelled using finite element analysis based on computed tomography and micro-computed tomography images, respectively. The first model was experimentally validated by occlusion analysis and subjected to the muscle loads produced during hard and soft chewing, traumatic cuspid occlusion, grinding, clenching, and simultaneous grinding and clenching. The occlusal forces computed by the first model were subsequently applied to the single tooth model to evaluate damage to the collagen network and the extracellular matrix of the PDL.

RESULTS

Early occlusal contact on the left cuspid tooth guided the mandible to the more occluded side (16.5% greater in the right side) and absorbed most of the lateral load. The intrusive occlusal loads on the posterior teeth were 0.77-13.3% greater than those on the cuspid. According to our findings, damage to the collagen network and the extracellular matrix of the PDL could occur in traumatic and grinding conditions, mainly due to fibre overstretching (>60%) and interstitial fluid overpressure (>4.7 kPa), respectively.

CONCLUSIONS

Our findings provide important biomechanical insights into the determination of damage mechanisms which are caused by mechanical loading and the key role of the porous-fibrous behaviour of the PDL in parafunctional and traumatic loading scenarios. Besides, the 3D loading conditions computed from occlusal contacts will help future studies in the design of new orthodontics appliances and encourage the application of computing methods in medical practice.

Physics of growing biological tissues: the complex cross-talk between cell activity, growth and resistance

For many organisms, shapes emerge from growth, which generates stresses, which in turn can feedback on growth. In this review, theoretical methods to analyse various aspects of morphogenesis are discussed with the aim to determine the most adapted method for tissue mechanics. We discuss the need to work at scales intermediate between cells and tissues and emphasize the use of finite elasticity for this. We detail the application of these ideas to four systems: active cells embedded in tissues, brain cortical convolutions, the cortex of Caenorhabditis elegans during elongation and finally the proliferation of epithelia on extracellular matrix. Numerical models well adapted to inhomogeneities are also presented. This article is part of the theme issue 'Rivlin's legacy in continuum mechanics and applied mathematics'.

A porous fibrous hyperelastic damage model for human periodontal ligament: Application of a microcomputerized tomography finite element model

The periodontal ligament (PDL) is a soft biological tissue that connects the tooth with the trabecular bone of the mandible. It plays a key role in load transmission and is primarily responsible for bone resorption and most common periodontal diseases. Although several numerical studies have analysed the biomechanical response of the PDL, most did not consider its porous fibrous structure, and only a few analysed damage to the PDL. This study presents an innovative numerical formulation of a porous fibrous hyperelastic damage material model for the PDL. The model considers two separate softening phenomena: fibre alignment during loading and fibre rupture. The parameters for the material model characterization were fitted using experimental data from the literature. Furthermore, the experimental tests used for characterization were computationally modelled to verify the material parameters. A finite element model of a portion of a human mandible, obtained by microcomputerized tomography, was developed, and the proposed constitutive model was implemented for the PDL. Our results confirm that damage to the PDL may occur mainly because of overpressure of the interstitial fluid, while large forces must be applied to damage the PDL fibrous network. Moreover, this study clarifies some aspects of the relationship between PDL damage and the bone remodelling process.

ORTHODONTIC PROCESS SAFETY EVALUATION BASED ON PERIODONTAL LIGAMENT CAPILLARY PRESSURE AND OGDEN MODEL

Taking the lower maxillary incisors as an example and the orthodontic forces along the near–far middle direction, the orthodontic forces along the crown–root direction and the orthodontic moment around the tongue–cheek direction as loading condition, the biomechanical simulation of the tooth is carried out by the method of finite element simulation in this paper. The CT images of the skull are segmented and denoised by Mimics. The solid models of teeth, periodontal ligament (PDL), alveolar bone and brackets are established by Gomagic and Solidworks. The material characteristics of the PDL are defined by the two-order Ogden hyperelastic model. Taking the PDL capillary pressure as a criterion for orthodontic safety, combined with the stress response of PDL, the safe orthodontic force range of mandibular central incisors is obtained by ANSYS finite element software.

Approach towards the porous fibrous structure of the periodontal ligament using micro-computerized tomography and finite element analysis

The periodontal ligament (PDL) is a porous and fibrous soft tissue situated around the tooth, which plays a key role in the transmission of loads from the tooth to the alveolar bone of the mandible. Although several studies have tried to characterize its mechanical properties, the behaviour of this tissue is not clear yet. In this study, a new simulation methodology based on a material model which considers the contribution of porous and fibrous structure with different material model formulations depending on the effort direction is proposed. The defined material model was characterized by a non-linear approximation of the porous fibrous matrix to experimental results obtained from samples of similar species and was validated by rigorous test simulations under tensile and compressive loads. The global PDL response was also validated using the parameters of the characterization in a finite element model of full human canine tooth obtained by micro-tomography. The results suggest that the porous contribution has high influence during compression because the bulk modulus of the material depends on the ability of interstitial fluid to drain. On the other hand, the collagen fibres running along the load direction are the main responsible of the ligament stiffness during tensile efforts. Thus, a material model with distinct responses depending of the load direction is proposed. Furthermore, the results suggest the importance of considering 3D finite element models based of the real morphology of human PDL for representing the irregular stress distribution caused by the coupling of complex material models and irregular morphologies.

RESEARCH OF MAXILLARY IMPACTED CANINE IN ORTHODONTIC TREATMENT BASED ON NANOINDENTATION EXPERIMENTS

This study selected the maxillary labial impacted canine as the research object to build the model of periodontal ligament (PDL) and simulate the process of orthodontic treatment. This paper obtained stress–strain curve by calculating and analyzing the data of nanoindentation experiments. The parameters were identified through curve fittings by ABAQUS. The fitting results show that the third-order Ogden model is in good agreement with the experimental data which demonstrate that the third-order Ogden model is able to reflect the material properties of the PDL. In this paper, orthodontic process of the maxillary labial impacted canine was simulated. The results show that inside and outside surfaces of PDL all have stress variation, the stress on the root apex and dental cervix of PDL is relatively large, the maximum appears at dental cervix and the minimum appears close to tooth impedance center.

MECHANICAL CHARACTERIZATION OF ANIMAL DERIVED GRAFTS FOR SURGICAL IMPLANTATION

Bioprostheses obtained from animal models are often adopted in abdominal surgery for repair and reconstruction. The functionality of these prosthetic implants is related also to their mechanical characteristics that are analyzed here. This work illustrates a constitutive model to describe the short-term mechanical response of Permacol[Formula: see text] bioprostheses. Experimental tests were developed on tissue samples to highlight mechanical non-linear characteristics and viscoelastic phenomena. Uni-axial tensile tests were developed to evaluate the strength and strain stiffening. Incremental uni-axial stress relaxation tests were carried out at nominal strain ranging from 10% to 20% and to monitor the stress relaxation process up to 400[Formula: see text]s. The constitutive model effectively describes the mechanical behavior found in experimental testing. The mechanical response appears to be independent on the loading direction, showing that the tissue can be considered as isotropic. The viscoelastic response of the tissue shows a strong decay of the stress in the first seconds of the relaxation process. The investigation performed is aimed at a general characterization of the biomechanical response and addresses the development of numerical models to evaluate the biomechanical performance of the graft with surrounding host tissues.

Expression analysis of α-smooth muscle actin and tenascin-C in the periodontal ligament under orthodontic loading or in vitro culture

α-smooth muscle actin (α-SMA) and tenascin-C are stress-induced phenotypic features of myofibroblasts. The expression levels of these two proteins closely correlate with the extracellular mechanical microenvironment. We investigated how the expression of α-SMA and tenascin-C was altered in the periodontal ligament (PDL) under orthodontic loading to indirectly reveal the intrinsic mechanical microenvironment in the PDL. In this study, we demonstrated the synergistic effects of transforming growth factor-β1 (TGF-β1) and mechanical tensile or compressive stress on myofibroblast differentiation from human periodontal ligament cells (hPDLCs). The hPDLCs under higher tensile or compressive stress significantly increased their levels of α-SMA and tenascin-C compared with those under lower tensile or compressive stress. A similar trend was observed in the tension and compression areas of the PDL under continuous light or heavy orthodontic load in rats. During the time-course analysis of expression, we observed that an increase in α-SMA levels was matched by an increase in tenascin-C levels in the PDL under orthodontic load in vivo. The time-dependent variation of α-SMA and tenascin-C expression in the PDL may indicate the time-dependent variation of intrinsic stress under constant extrinsic loading.

Viscoelasticity of periodontal ligament: an analytical model

BACKGROUND

Understanding of viscoelastic behaviour of a periodontal membrane under physiological conditions is important for many orthodontic problems. A new analytic model of a nearly incompressible viscoelastic periodontal ligament is suggested, employing symmetrical paraboloids to describe its internal and external surfaces.

METHODS

In the model, a tooth root is assumed to be a rigid body, with perfect bonding between its external surface and an internal surface of the ligament. An assumption of almost incompressible material is used to formulate kinematic relationships for a periodontal ligament; a viscoelastic constitutive equation with a fractional exponential kernel is suggested for its description.

RESULTS

Translational and rotational equations of motion are derived for ligament's points and special cases of translational displacements of the tooth root are analysed. Material parameters of the fractional viscoelastic function are assessed on the basis of experimental data for response of the periodontal ligament to tooth translation. A character of distribution of hydrostatic stresses in the ligament caused by vertical and horizontal translations of the tooth root is defined.

CONCLUSIONS

The proposed model allows generalization of the known analytical models of the viscoelastic periodontal ligament by introduction of instantaneous and relaxed elastic moduli, as well as the fractional parameter. The latter makes it possible to take into account different behaviours of the periodontal tissue under short- and long-term loads. The obtained results can be used to determine loads required for orthodontic tooth movements corresponding to optimal stresses, as well as to simulate bone remodelling on the basis of changes in stresses and strains in the periodontal ligament caused by such movements.

Experiment and hydro-mechanical coupling simulation study on the human periodontal ligament

In this paper, a new method involving an experiment in vivo and hydro-mechanical coupling simulations was proposed to investigate the biomechanical property of human periodontal ligament (PDL). Teeth were loaded and their displacements were measured in vivo. The finite element model of the experiment was built and hydro-mechanical coupling simulations were conducted to test some PDL's constitutive models. In the simulations, the linear elastic model, the hyperfoam model, and the Ogden model were assumed for the solid phase of the PDL coupled with a model of the fluid phase of the PDL. The displacements of the teeth derived from the simulations were compared with the experimental data to validate these constitutive models. The study shows that a proposed constitutive model of the PDL can be reliably tested by this method. Furthermore, the influence of species, areas, and the fluid volume ratio on PDL's mechanical property should be considered in the modeling and simulation of the mechanical property of the PDL.