ORIGINAL RESEARCH article

Front. Physiol., 29 August 2016

Sec. Craniofacial Biology and Dental Research

Volume 7 - 2016 | https://doi.org/10.3389/fphys.2016.00364

Theoretical Considerations and a Mathematical Model for the Analysis of the Biomechanical Response of Human Keratinized Oral Mucosa

  • 1. Department of Prosthodontics, Faculty of Health Sciences, School of Dentistry, Aristotle University Thessaloniki, Greece

  • 2. Laboratory of Thin Films-Nanosystems and Nanometrology, School of Physics, Aristotle University Thessaloniki, Greece

  • 3. Department of Engineering, Cambridge Graphene Centre, Cambridge University Cambridge, UK

  • 4. Department of Electrical and Computer Engineering, School of Engineering, Aristotle University Thessaloniki, Greece

  • 5. Division of Graduate and Postgraduate Prosthodontics, Tufts University School of Dental Medicine Boston, MA, USA

Abstract

Removable complete and partial dentures are supported by the residual alveolar ridges consisting of mucosa, submucosa, periosteum, and bone. An understanding of the biomechanical behavior of the oral mucosa is essential in order to improve the denture-bearing foundations for complete and partially edentulous patients. The purpose of this paper was to examine the biomechanical behavior of the soft tissues supporting a removable denture and develop a model for that reason. Keratinized oral mucosa blocks with their underlying bone were harvested from the maxillary palatal area adjacent to the edentulous ridges of a cadaver. The compressive response of the oral mucosa was tested by using atomic force microscopy. The specimens were first scanned in order their topography to be obtained. The mechanical properties of the specimens were tested using a single crystal silicon pyramidal tip, which traversed toward the keratinized oral mucosa specimens. Loading-unloading cycles were registered and four mathematical models were tested using MATLAB to note which one approximates the force-displacement curve as close as possible: a. spherical, b. conical, c. third order polynomial, d. Murphy (fourth order polynomial, non-linear Hertzian based). The third order polynomial model showed the best accuracy in representing the force-displacement data of the tested specimens. A model was developed in order to analyze the biomechanical behavior of the human oral keratinized mucosa and obtain information about its mechanical properties.

Introduction

Removable complete and partial dentures are supported by the residual alveolar ridges consisting of mucosa, submucosa, periosteum, and bone. This supporting area has been calculated to be 23 cm2 for the maxilla, and about 12.25 cm2 for the edentulous mandible (Hobkirk and Zarb, 2013). It should be mentioned however that these figures vary, depending on the size of the maxilla or mandible and the amount of resorption after teeth extraction. Nevertheless, these numbers are substantially lower than the supporting mechanism of the teeth—i.e., the periodontal ligament—which is about 45 cm2, for each dental arch (Hobkirk and Zarb, 2013).

Besides the smaller area of the denture bearing surface in comparison with that of the teeth, there are some other distinct differences too. These include the involved sensory mechanisms and the anatomical features of each structure. The periodontal ligament is connective tissue with a thickness of 0.15–0.35 mm, consisting of collagen, oxytalan and eulanin fibers, glycosaminoglycans and blood vessels ranging between 4 and 47% of the total tissue volume (Jonas and Riede, 1980; Blaushild et al., 1992; Johnson and Pylypas, 1992; Embery et al., 1995; Sloan and Carter, 1995; Michalakis et al., 2012). The periodontal ligament is organized into six different groups of fibers which are not unidirectionally distributed: 1. transeptal, 2. alveolar crest, 3. horizontal, 4. oblique, 5. apical, and 6. interradicular, which are present only between the roots of multirooted teeth (Carranza, 1990; Lindhe et al., 2003). The existence of many blood vessels into the periodontal ligament and the hemodynamic pressure that these vessels exert affects its biomechanical response (Kristiansen and Heyreaas, 1989; Sims, 1995; Ioi et al., 2002a,b). The biomechanical response of the periodontal ligament to occlusal loads is not clear (Caputo and Standlee, 1987). Three focal hypotheses have been made in the past, in order to describe the way in which the periodontal ligament supports the tooth: (i) the tensional mechanism model, supporting the idea that the fibers have a wavy configuration and consequently load transmission from the tooth to the neighboring alveolar bone occurs through a gradual unfolding of these fibers (Mühlemann, 1951; Picton, 1965, 1969), (ii) the viscoelastic model, considering that tooth movement within the socket is controlled more by the vascular elements than by the fibers (Bien, 1966; Fung, 1973; Natali et al., 2004), (iii) the collagenous thixotropic model, supporting the notion that tooth support is possible because of the periodontal ligament's thixotropic gel properties (Kardos and Simpson, 1979, 1980).

The oral mucosa covering the hard palate and the attached gingiva is termed masticatory mucosa and consists of the stratified squamous epithelium at the surface and the lamina propria which lies deeper. The stratified squamous epithelium consists of four layers, which—from most superficial to deepest—are: (a) stratum corneum, (b) stratum granulosum, (c) stratum spinosum, (d) stratum basale (Nanci, 2013). Lamina propria is connective tissue, which is composed of cells, mainly fibroblasts, and an extracellular matrix, consisting of a ground substance and fibers. Fibroblasts are responsible for the secretion of collagen and other elements of the extracellular matrix. The ground substance of the matrix is composed of glycoproteins, glysosaminoglycans and proteoglycans, while the fibers are mainly collagen, providing tensile strength and flexibility to the tissue and elastic fibers, contributing resiliency. Below the oral mucosa of the attached gingiva and the hard palate lies the mucoperiosteum with dense collagenous connective tissue attaching directly to the periosteum. The mucoperiosteum contains fat and salivary glands (Slavkin and Bavetta, 1972; Newcomb, 1981; Clausen et al., 1983; Meyer et al., 1984; Dahllöf et al., 1986; Bourke et al., 2000). An understanding of the biomechanical behavior of the oral mucosa is essential in order to improve the denture-bearing foundations for complete and partially edentulous patients, by better managing traumatized tissues and giving instructions to patients regarding the time which is required for tissues to recover, after applying occlusal loads during daytime. Furthermore, finite element analysis models studying the deformation of oral mucosa under occlusal loading require use of an equation, which unfortunately is not supplied by the bibliography. Nevertheless, finite element analyses of the oral mucosa have been performed in the past and several material models have been adopted in order mucosal behavior to be interpreted. These include linear elastic, biphasic, multi-phasic elastic, and hyperelastic models (Chen et al., 2015). Additionally, knowledge of oral mucosa biomechanics can be helpful in fabricating dental materials with similar or complementary behavior to that of oral tissues (Saitoh et al., 2010; Hong et al., 2012).

Although, numerous articles have been published on the biomechanics of the periodontal ligament (Kurashima, 1965; Komatsu and Viidik, 1966; Daly et al., 1974; Wills et al., 1976; Atkinson and Ralph, 1977; Wills and Picton, 1978; Dorow et al., 2003; Natali et al., 2003; Bergomi et al., 2010, 2011), the research on the biomechanics of keratinized oral mucosa is scarce (Keilig et al., 2009; Goktas et al., 2011).

The purpose of this study was to examine the biomechanical behavior of the soft tissues supporting a removable complete denture and develop a model for that reason.

Materials and methods

Tissue preparation

This in vitro study was performed in accordance with the guidelines of the Declaration of Helsinki and the research protocol was approved by the Ethics Committee of the Aristotle University (256/06-07-2011), prior to initiation.

Eight 12 × 8 × 8 mm keratinized oral mucosa blocks with their underlying bone were provided by the Laboratory of Anatomy of the Medical Faculty of the Health Sciences School of the Aristotle University of Thessaloniki (Goktas et al., 2011; Herris et al., 2013). The specimens were harvested from maxillary edentulous areas, by using a low speed 0.2 mm thickness diamond disc (Thin Flex X929-7 TP; Abrasive Technology Inc, Lewis Center, OH, USA) under continuous saline irrigation (Figure 1). To prevent dehydration, the specimens were then stored in a 10% neutral buffered formalin solution (water 91.9–92%, formaldehyde 4%, methyl alcohol 1–2%, sodium phosphate dibasic 0.65%, sodium phosphate monobasic, monohydrate 0.4%) until the mechanical analysis testing, which took place 1 h after.

Figure 1

The compressive response of the oral mucosa was the only biomechanical characteristic tested, by using an Atomic Force Microscope (Solver P47H; NT-MDT Co., Moscow, Russia). A standard square based pyramidal single crystal silicon (Si) tip (NSG 10; NT-MDT Co., Moscow, Russia) with a < 10 nm typical curvature tip radius and a lateral surface of 500 μm was used. A nominal spring constant of 12.9 0.06 N/m was used after calculating it by using the Sader method (Sader et al., 1999). The length of the cantilever was 95 ± 5 μm, the width was 30 ± 3 μm, while the thickness was 2 ± 0.5 μm.

Before initiation of the experimental procedures a calibration of the cantilever took place. First the deflection was converted to force using the Hooke's law and then the response of the cantilever was subtracted from the measurement. For this purpose a force-displacement curve in a hard surface was acquired, followed by a force-displacement curve on the sample. The difference of these measurements supplied the indentation depth.

The specimens were then fixed on glass slabs with Histoacryl topical tissue adhesive (B. Braun Corp., Melsungen, Germany), which was used according to manufacturer's instructions. The fixing luting agent was left to set for a period of 1 min before initiation of the measurements. For each specimen a time-period of 3 min was required for the preparation and testing procedures.

Specimens topography

The specimens were first scanned in order their topography to be obtained. This can be achieved in two ways: (1) the contact mode and, (2) the tapping (semicontact) mode (Ethier and Simmons, 2008).

In the contact mode the probe exerts a constant force to the specimens, which has as a result the application of large lateral forces and therefore a possible deformation of the specimen (Ethier and Simmons, 2008). Thus, the tapping mode, in which the cantilever is either magnetically or acoustically driven, was selected. (Figure 2) During the scanning period the oscillating tip of the cantilever was moving in very close proximity to the surface of the specimen, touching it regularly.

Figure 2

The oscillations of the cantilever are due to repulsive and attracting forces, and they have been discussed previously in the literature (Goodman and Garcia, 1991; García and San Paulo, 1999). The first ones are short range forces with an exponential decaying and can be considered as Pauli's exclusion principle interaction, electron-electron Coulomb interaction and hard sphere repulsion. The second ones are long range forces, including van der Waals interactions, electrostatic and chemical forces. For the interaction between the tip of the cantilever and the specimen's surface, the van der Waals forces and van der Waals potential obey to the following (Hamaker, 1937):

Where, A is the Hamaker constant related to the material, R is the sphere radius and d is the separation between the tip and the half-space surface (Argento and French, 1996).

The scanning rate was set at ~5,6 μm/s step 11 nm Hz, and the specimens' surface roughness was determined by the two following equations:

Mean absolute deviation surface roughness (Ra)

Where Zmean represents the mean height, as this was calculated over the entire area of the biologic specimen, discretized in the grid of Z(i,j), i = 1,…,Nx and j = 1,…,Ny

The Root Mean Square (RMS) surface roughness

which represents the average deviation between the height and the mean surface.

Additionally, Ry (peak to valley) measurements were recorded. Two- and three-dimensional images of the specimens' topography were obtained (Figure 3).

Figure 3

Indentation test

The mechanical properties of the specimens were then tested. The Si pyramidal tip traversed vertically toward the keratinized oral mucosa specimens and the deflection of the cantilever was measured. As the rigid tip moved toward the soft biologic specimen the latter deflected around the probe. This problem of contact mechanics is based on the Hertz theory and the deflection of the cantilever arm is given by the following (Haga et al., 2002; Ethier and Simmons, 2008):

Where z is the vertical deflection, zo denotes the height of the probe where the force F applied to the biologic specimen becomes non-zero, kc is the stiffness of the cantilever, E is the elastic modulus of the biologic material in N/m−2, v is the Poisson's ratio of the material (due to high water content Poisson's ratio for most biological specimens is considered to be 0.5), and α is the face angle for the silicon-nitride cantilever (Dimitriadis et al., 2002; Ethier and Simmons, 2008; Figure 4).

Figure 4

The force applied by a four sided pyramidal tip is given by:

A vertical oscillating frequency of ~331 kHz was used and data were recorded at multiple sites of the tested specimens (Weisenhorn et al., 1992, 1993a,b; Stolz et al., 2004; Rahmat and Hubert, 2010; Figure 5). Each loading-unloading cycle lasted 0.3 s.

Figure 5

Force calibration and displacement data

Identification of the first contact between the tip of the cantilever and the material tested presents a challenge in the indentation tests of soft biologic tissues, when using atomic force microscopy (Stolz et al., 2004; Herris et al., 2013). The abrupt change in the force-development curve, due to the generation of repulsive forces between the two interacting bodies, was taken as the initial contact point (z0), and the corresponding force (F0) at that point was zero.

The indentation depth was given by the following (Rahmat and Hubert, 2010; Herris et al., 2013):

where, z denotes the displacement of the cantilever base, d0 represents the deflection of the cantilever at initial contact and d is the deflection of the cantilever. It should be pointed out that the cantilever's deflection was given by

where, F represents the force and k denotes the stiffness of the cantilever.

Room temperature (21 ± 1°C) and relative humidity (50 ± 10%) were recorded throughout the experiments.

Results

The indentation tests were performed only on the coronal surface of the specimens. Indentation tests on sagittal planes were not performed due to their low clinical application.

Data from the indentation tests of the tested specimens was collected and force-displacement diagrams have been obtained (Figure 5). It should be mentioned that the negative values in the curve correspond to the initial phase of the testing procedure, when the tip lands on the surface of specimen.

Five loading-unloading cycles were registered. The simplest fitting approach was selected. That consisted of a visual inspection of the force-displacement curve and identification of the initial contact point (zo, do), as already presented earlier (Lin et al., 2007). The area of the curve representing the noncontact region was ignored and for the contact region of the curve four mathematical models were tested using MATLAB (Mathworks Inc; Natick, MA, USA) to note which one approximates the force-displacement curve, as close as possible: a. spherical, b. conical, c. third order polynomial, d. Murphy (fourth order polynomial, non-linear Hertzian based), which is a fourth order polynomial model (Murphy et al., 2013; Table 1, Figure 6).

Table 1

Fit nameNo.1No.2No.3No.4No.5
exp = 3/20.1264030.0099390.0211020.0034260.004355
exp = 20.0474430.0853100.1174850.0608890.066963
Murphy0.0003800.0023660.0046340.0026580.01629
Polynomial 30.0003960.0000710.0001420.0001890.000023

Comparison of the numerical fit results including the goodness-of-fit statistic of the normalized mean square error, where the normalization is by the sample variance.

Figure 6

The goodness-of-fit measure presented in Table 1 is the normalized mean square error (NMSE), and was defined as follows:

where xi is the i-th observation of variable X, is the corresponding estimated value by the fitted model, is the sample mean of the n observations of X.

The spherical and the conical models did not approximate the force-displacement curves The model presented by Murphy et al. (2013), with the form

performed well for only one specimen.

The third order polynomial model of the form

approximated the curve very closely for all tested specimens.

Therefore, a proposition of a model which fits the experimental data better than the aforementioned models and the Hertz model (Herz, 1881), which is fully elastic, is attempted:

Murphy et al. (2013) have speculated that the Young's modulus of biological materials varies with displacement and is given by the following second order polynomial:

The authors hypothesized that the triad set k1, k2, and Eb governs the modulus of elasticity of the material they tested. Specifically k1 and k2 represent the non-linear region, while Eb represents the contact stiffness, which is in the elastic region of the force-displacement curve.

The coefficients Γ1, Γ2, and Γ3 of their fourth order polynomial (11) are:

Accordingly, in the present study E(δ) is given by the following equation:

where k1, k2, and Eb describe the oral mucosa's and mucoperiosteum's properties as a triad set, and

Eb represents the mucoperiosteum's contact stiffness and δ is the indentation depth.

A similar approach has been adopted by Murphy et al. (2013). Other researchers have come independently to the same conclusion (Herris et al., 2013).

The constant term P4 of (12) is omitted, as when F = 0, δ = 0.

Then the general polynomial model (12) becomes:

Therefore, by substituting E(δ) from (17) and (18) to (6), the latter becomes:

where F is the applied force, δ is the indentation depth, α is the face angle for the silicon-nitride cantilever, v is the Poisson's ratio of the mucosa, and k1, k2, and k3 are the moduli of elasticity for the stratified squamous epithelium, lamina propria and mucoperiosteum respectively.

Discussion

In the present ex vivo study AFM was used to test the mechanical properties of human keratinized oral mucosa. The specimens were harvested from the edentulous areas which support the maxillary complete denture (Hobkirk and Zarb, 2013). The force-development curves were analyzed and a third order polynomial model different than the classic elastic Hertz model (Herz, 1881) was developed in order to describe its biomechanical response.

In the past 25 years AFM has been used to study the mechanical properties of materials. This concept has also been used in the early 90's for “force-displacement” data collection from biological materials, including tendons, ligaments, muscle tissues and menisci (Burnham and Colton, 1989; Radmacher et al., 1992; Tao et al., 1992; Weiss et al., 2002; Sweigart et al., 2004; Yin and Elliot, 2004; Van Loocke et al., 2006, 2008; Villegas et al., 2007; Cheng and Gan, 2008; Ciarletta et al., 2008). A nano-indenter could have been used in this study, as well. However, nano-indenters have a resolution of approximately 100 nN, while the forces applied by the AFM can range from pico-Newtons to several hundreds of micro-Newtons by changing the stiffness of the cantilever. Thus, the sensitivity and the versatility of the AFM makes it an ideal tool for mechanical properties testing of biological materials (Stolz et al., 2009; Notbohm et al., 2012). Additionally, use of the tapping mode of the AFM prevents the distortion of the biologic specimen, as it is presented later in the discussion. In most of the published research, standard manufacturers' cantilevers and pyramidal tips have been employed, while modifications with attachment of silica microspheres have also been reported (Mahaffy et al., 2000; Dimitriadis et al., 2002). The solution for the microsphere tips is supplied by the original Hertz model which deals with the shallow contact between two spherical bodies (Herz, 1881). The original Hertzian theory has been used by many researchers who studied contact deformation and many modifications have been made in order to account for large deformations (Gao and Gao, 2000), viscoelasticity (Gillies et al., 2002; Yang et al., 2004; Attard, 2007; Chen et al., 2013), anisotropy (Batra and Jiang, 2008), multi-layered structures (Ai et al., 2002), and adhesive interactions (Cao et al., 2005; Yang, 2006; Ebenstein, 2011; Chen et al., 2013; Kohn and Ebenstein, 2013). In the present study a sharp pyramidal tip was employed. It should be mentioned that the Bilodeau solution applies for this case (Bilodeau, 1992).

The model developed in the present study is probably valid only when the forces are applied to the coronal part of the oral keratinized mucosa, as it has been proven that tissues are both inhomogeneous and anisotropic. Biological specimens' anisotropic properties have been demonstrated with nonlinear laser scanning microscopy, which has been used for elastin and collagen distribution—within the specimens—imaging (Herris et al., 2013). Additionally their mechanical properties are site and direction-specific dependent (Stolz et al., 2004). It has been shown in the past that the indentation elastic moduli of biological specimens differ depending on the depth, increasing from the superficial to deep layers (Herris et al., 2013). Furthermore, it has been demonstrated that both the Young's modulus and Poisson's ratio change during different development stages of the cell (Zhang et al., 2009). Three different moduli of elasticity were incorporated to the model developed in this study, with k1 being the modulus of elasticity of stratified squamous epithelium, k2 of lamina propria and k3 of mucoperiosteum. This is in accordance with the present experimental results and previous research findings, which have demonstrated that the elastic moduli increase from the superficial to the deep layer (Herris et al., 2013; Murphy et al., 2013). Furthermore, the cells and the extracellular matrix present different elastic modulus and may contribute to the observed inhomogeneity (Trickey et al., 2006; Han et al., 2011). Third-degree polynomial models have also proven to be valid for mechanical forces acting on cell biomembranes, as well as in other anatomical structures as the descending aorta (Stefanadis et al., 2000; Zhang et al., 2009). The third degree polynomial model without a constant term which was tested in the present study not only approximated better the force-displacement curve of the keratinized oral mucosa, but it also provided a simpler model than the fourth order polynomial presented by Murphy et al. (2013). It should be mentioned however that different biological tissues have been examined in these two studies.

As already mentioned, before testing the specimens were kept in a buffered 10% formalin solution until testing, which occurred within an hour. The 10% formalin solution is actually a 4% formaldehyde, which is by definition 1.3 molar. A totally unbuffered formaldehyde solution exerts an osmotic pressure of about 1300 mO. Isotonic salt solutions present osmolarities of 250–350 mO. Therefore, it is expected that formalin diffuses into tissues faster. Formaldehyde has a molecular weight of 30 and it is expected to penetrate the tissues fast. Nevertheless, fixation actually takes a relatively long time. It has been estimated that a time period of 6–16 h, depending on the specimen, is needed. Furthermore, since all specimens received the same treatment simultaneously it can be assumed that formalin penetration was uniform in all specimens. Macroscopically, no swelling was noticed within the 1-h period between specimen harvesting and testing (Bono et al., 2001; Thavarajah et al., 2012).

The human oral keratinized tissue tested in the present study was bonded with cyanoacrylate cement to a glass slab, since there is scientific evidence that force-displacement results depend strongly on whether or not the specimens are attached to the substrate (Yang, 1998). A glass substrate was used as its modulus of elasticity is 50 GPa, which is much higher than that of biologic specimens, while the modulus of elasticity of the pyramidal tip employed was 150 GPa (Grafström et al., 1993; Weisenhorn et al., 1993a). Therefore, the deformation of both the tip and the substrate can be considered as negligible.

The loading-unloading curves of the specimens tested did not coincide, suggesting a viscoelastic behavior (Lakes, 1999). Furthermore, the hysteresis of the loading-unloading curves, indicates dissipation of energy. During loading, the area beneath represents the energy stored, while during unloading, this area represents the energy recovered. Quantification of the hysteresis can be performed by introduction of the plasticity index η. In experiments involving solid objects this parameter characterizes the elastic/plastic behavior of the material, when external forces are applied. The value of plasticity index can range between 0, indicating a fully elastic material, and 1, when the material displays a fully plastic behavior (Briscoe et al., 1998; Klymenko et al., 2009).

The tapping mode was selected for the present study, as it presents more advantages than the contact and the non-contact imaging methods. The tapping mode eliminates the lateral, frictional forces transmitted from the tip of the cantilever to the specimen's surface. In this mode the probe oscillates with sufficient amplitude to prevent it being trapped by adhesive meniscus forces from the contaminant layer (e.g., water), and it makes intermittent contact with the surface of the specimen (Manning et al., 2003; Rogers et al., 2004). As a result, tissue damage is minimal, if any, and recovery is guaranteed. The oscillation frequency usually ranges between 50,000 and 500,000 cycles per second. An oscillation frequency of 311,000 cycles was used in the present study. With this frequency the tip-specimen adhesion forces are overcome and only vertical forces are applied. The high frequency oscillations of the tip may be regarded as a disadvantage of the method, as they are not comparable to the human chewing frequency. It should be mentioned however that, neither the forces applied by the tip to the specimen are comparable to mastication forces. Atomic force microscopy, like all laboratory techniques, has some limitations.

It should be pointed out however that, in studies like the present one, a fundamental assumption is made: the mechanical response of biological materials relies on contribution of different structures which act in sequence. In this way interpretation of the mechanical properties of materials which present a hierarchy can be made (Bonilla et al., 2015).

Development of a model describing the behavior of oral mucosa under mechanical forces is instrumental for the knowledge of its mechanical properties, such as the Young's modulus, and an understanding of how masticatory function is affected by mechanical interactions. Moreover, this modeling assists in fabricating biomaterials (e.g., tissue conditioners) which will act in a similar way with, or in a complimentary way to the oral keratinized mucosa.

It should be mentioned that ex-vivo studies performed with AFM present certain drawbacks, including the identification of the most appropriate area for force application, the high frequency oscillation of the cantilever, the possible minor dehydration of the tissue, the molecular roughness of the pyramidal tip and the uncertainty of the first contact between the tip and the biological specimen. These have also been acknowledged by other authors and possibly affect measurement's accuracy (Vinckier and Semenza, 1998).

Further studies are needed to determine the contribution of each layer, as well as, that of the cells and the extracellular matrix in the biomechanical behavior of the oral mucosa.

Conclusion

Within the limitations of the present study, the following conclusions can be made regarding the force-displacement data of the human oral keratinized mucosa:

  • The third order polynomial model examined in the present study showed a very good accuracy.

  • The Murphy model (fourth order) performed well for only one specimen.

  • The spherical and the conical models did not approximate the force-displacement curves.

  • A mathematical model for the analysis of the biomechanical response of human keratinized oral mucosa was developed.

Author note

The results of this paper were presented at the IADR General Session, March 20–23, 2013, Seattle, WA, USA, and at the European Society of Biomechanics, August 25–28, 2013, Patras, Greece.

Statements

Author contributions

All authors (AT, PK, MS, SK, DK, SL, ON, AP, and KM) have made substantial contributions to the conception and/or design of the work; or the acquisition, analysis, or interpretation of data for the work; and greatly assisted in drafting the work or revising it critically for important intellectual content; and they approved the submitted version; and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    AiZ. Y.YueZ. Q.ThamL. G.YangM. (2002). Extended Sneddon and Muki solutions for multilayered elastic materials. Int. J. Eng. Sci.40, 14531483. 10.1016/S0020-7225(02)00022-8

  • 2

    ArgentoC.FrenchR. H. (1996). Parametric tip model and force-distance relation for Hamaker constant determination from atomic force microscopy. J. Appl. Phys.80, 60816090. 10.1063/1.363680

  • 3

    AtkinsonH. F.RalphW. J. (1977). In vitro strength of the human periodontal ligament. J. Dent. Res.56, 4852. 10.1177/00220345770560011001

  • 4

    AttardP. (2007). Measurement and interpretation of elastic and viscoelastic properties with the atomic force microscope. J. Phys. Condens. Matt.19, 47320161. 10.1088/0953-8984/19/47/473201

  • 5

    BatraR. C.JiangW. (2008). Analytical solution of the contact problem of a rigid indenterand an anisotropic linear elastic layer. Int. J. Solids Struct.45, 58145830. 10.1016/j.ijsolstr.2008.06.016

  • 6

    BergomiM.CugnomiJ.BotsisJ.BelserU. C.WiskottH. W. (2010). The role of the fluid phase in the viscous response of bovine periodontal ligament. J. Biomech.43, 11461152. 10.1016/j.jbiomech.2009.12.020

  • 7

    BergomiM.CugnomiJ.GalliM.BotsisJ.BelserU. C.WiskotH. W. (2011). Hydro-mechanical coupling in the periodontal ligament: a porohyperelastic finite element model. J. Biomech.44, 343810.1016/j.jbiomech.2010.08.019

  • 8

    BienS. M. (1966). Hydrodynamic damping of tooth movement. J. Dent. Res.45, 907914. 10.1177/00220345660450036701

  • 9

    BilodeauG. (1992). Regular pyramid punch problem. J Appl Mech59, 519523. 10.1115/1.2893754

  • 10

    BlaushildN.MihaeliY.SteigmanS. (1992). Histomorphometric study of the periodontal vasculature of the rat incisor. J. Dent. Res.71, 19081912. 10.1177/00220345920710121001

  • 11

    BonillaM. R.StokesJ. R.GidleyM. J.YakubovG. E. (2015). Interpreting atomic force microscopy nanoindentation of hierarchical biological materials using multi-regime analysis. Soft Matt.11, 12811292. 10.1039/C4SM02440K

  • 12

    BonoC.RenardR.SabatinoC.LevineR.TornettaP. (2001). The effects of varied concentrations of formalin on the tensile strength of cortical bone: should embalmed bone ever be used for biomechanical testing, in Orthopedics Research Society, 47th Annual Meeting February (San Francisco, CA).

  • 13

    BourkeK. A.HaasseH.LiH.DaleyT.BartoldP. M. (2000). Distribution and synthesis of elastin in porcine gingiva and alveolar mucosa. J. Periodont. Res.35, 361368. 10.1034/j.1600-0765.2000.035006361.x

  • 14

    BriscoeB. J.FioriL.PelilloE. (1998). Nano-indentation of polymeric surfaces. J. Phys. D Appl. Phys.31, 23952405. 10.1088/0022-3727/31/19/006

  • 15

    BurnhamN. A.ColtonR. J. (1989). Measuring the nanomechanical properties and surface forces of materials using the AFM. J. Vac. Sci. Technol.7, 29062913. 10.1116/1.576168

  • 16

    CaoY. F.YangD. H.SoboyejoyW. (2005). Nanoindentation method for determining the initial contact and adhesion characteristics of soft polydimethylsiloxane. J. Mater. Res.20, 20042011. 10.1557/JMR.2005.0256

  • 17

    CaputoA. A.StandleeJ. P. (1987). Biomechanics in Clinical Dentistry. Chicago, IL: Quintessence.

  • 18

    CarranzaF. A. (1990). Glickman's Clinical Periodontology.Philadelphia, PA: W. B. Saunders.

  • 19

    ChenJ.AhmadR.LiW.SwainM.LiQ. (2015). Biomechanics of oral mucosa. J. R. Interface12:20150325. 10.1098/rsif.2015.0325

  • 20

    ChenZ.DiebelsS.PeterN. J.SchneiderA. S. (2013). Identification of finite viscoelasticity and adhesion effects in nanoindentation of a soft polymer by inverse method. Comput. Mater. Sci.72, 127139. 10.1016/j.commatsci.2013.01.040

  • 21

    ChengT.GanR. Z. (2008). Mechanical properties of anterior malleolar ligament from experimental measurement and material modelling analysis. Biomech. Model. Mechanobiol.7, 387394. 10.1007/s10237-007-0094-x

  • 22

    CiarlettaP.DarioP.MiceraS. (2008). Pseudo-hyperelastic model of tendon hysteresis from adaptive recruitment of collagen type I fibrils. Biomaterials29, 764770. 10.1016/j.biomaterials.2007.10.020

  • 23

    ClausenH.VedtofteP.MoeD.DabelsteenE. (1983). Keratin pattern in human and buccal and hard palate mucosa. Scand. J. Dent. Res.91, 411413. 10.1111/j.1600-0722.1983.tb00838.x

  • 24

    DahllöfG.MondèerT.ReinholtF. P.WilkströmB.HjerpeA. (1986). Proteoglycans and glycosaminoglycans in phenytoin –induced gingival overgrowth. J. Periodont. Res.21, 132110.1111/j.1600-0765.1986.tb01432.x

  • 25

    DalyC. H.NichollsJ. I.KyddW. L.NansenP. D. (1974). The response of the human periodontal ligament to torsional loading- I. Experimental methods. J. Biomech.7, 517522. 10.1016/0021-9290(74)90086-4

  • 26

    DimitriadisE. K.HorkayF.MarescaJ.KacharB.ChadwickR. S. (2002). Determination of elastic moduli of thin layers of soft material using the atomic force microscope. Biophys. J.82, 27982810. 10.1016/S0006-3495(02)75620-8

  • 27

    DorowC.KristinN.SanderF. G. (2003). Determination of the mechanical properties of the periodontal ligament in a uniaxial tensional experiment. J. Orofac. Orthop.64, 100107. 10.1007/s00056-003-0225-7

  • 28

    EbensteinD. M. (2011). Nano-JKR force curve method overcomes challenges of surface detection and adhesion for nanoindentation of a compliant polymer in air and water. J. Mater. Res.26, 10261035. 10.1557/jmr.2011.42

  • 29

    EmberyG.WaddingtonR.HallR. (1995). The ground substance of the periodontal ligament, in The Periodontal Ligament in Health and Disease, eds BerkovitzB. K. B.MoxhamB. J.NewmanH. N. (London: Mosby-Wolfe), 83106.

  • 30

    EthierC. R.SimmonsC. A. (2008). Introductory Biomechanics. From Cells to Organisms.Cambridge: University Press.

  • 31

    FungY. C. (1973). Biorheology of soft tissues. Biorheology10, 139155.

  • 32

    GaoY. C.GaoT. J. (2000). Large deformation contact of a rubber notch with a rigid wedge. Int. J. Solids Struct.37, 43194334. 10.1016/S0020-7683(99)00191-2

  • 33

    GarcíaR.San PauloA. (1999). Attractive and repulsive tip-sample interaction regimes in tapping-mode atomic force microscopy. Phys. Rev. B.60, 49614967. 10.1103/PhysRevB.60.4961

  • 34

    GilliesG.PrestidgeC. A.AttardP. (2002). An AFM study of the deformation and nanorheology of cross-linked PDMS droplets. Langmuir18, 16741679. 10.1021/la011461g

  • 35

    GoktasS.DmytrykJ. J.McFetridgeP. S. (2011). Biomechanical behavior of oral soft tissues. J. Periodontol.82, 11781186. 10.1902/jop.2011.100573

  • 36

    GoodmanF. O.GarciaN. (1991). Roles of the attractive and repulsive forces in atomic-force microscopy. Phys. Rev. B.43, 47284731. 10.1103/PhysRevB.43.4728

  • 37

    GrafströmS.NeitzertM.HagenT.AckermannJ.NeumanR.PtobstO.et al. (1993). The role of topography and friction for the image contrast in lateral force microscopy. Nanotechnology4, 143151. 10.1088/0957-4484/4/3/003

  • 38

    HagaH.SasakiS.KawabataK.ItoE.UshikiT.SambongiT. (2002). Elasticity mapping of living fibroblasts by AFM and immunofluresence observation of the cytoskeleton. Ultramicroscopy82, 253258. 10.1016/S0304-3991(99)00157-6

  • 39

    HamakerH. C. (1937). The London-van der Waals attraction between spherical particles. Physica4, 10581072. 10.1016/S0031-8914(37)80203-7

  • 40

    HanL.GrodzinskyA. J.OrtizC. (2011). Nanomechanics of the cartilage extracellular matrix. Annu. Rev. Mater. Res.41, 133168. 10.1146/annurev-matsci-062910-100431

  • 41

    HerrisH. K.MiriA. K.TripathyU.BarthelatF.MongeauL. (2013). Indentation of poroviscoleastic vocal fold tissue using an atomic force microscope. J. Mech. Behav. Biomed. Mater.28, 383392. 10.1016/j.jmbbm.2013.05.026

  • 42

    HerzH. (1881). Über die berührung fester elastischer körper. J. die reine und Angewandte Mathematik92, 156171.

  • 43

    HobkirkJ. A.ZarbG. (2013). The edentulous state, in Prosthodontic Treatment for Edentulous Patients. Complete Dentures and Implant-Supported Prostheses, 13th Edn., eds ZarbG.HobkirkJ. A.EckertS. E.JacobR. F. (St. Louis, MO: Elsevier Mosby), 127.

  • 44

    HongG.MaedaT.MurataH.SasakiK. (2012). The dynamic viscoelasticity and plasticizer leachability of tissue conditioners. Gerodontology29, 284291. 10.1111/j.1741-2358.2012.00639.x

  • 45

    IoiH.MorishitaT.NakataS.NakashimaA.NandaR. S. (2002a). Evaluation of physiological tooth movements within clinically normal periodontal tissues by means of periodontal pulsation measurements. J. Periodont. Res.37, 110117. 10.1034/j.1600-0765.2001.00665.x

  • 46

    IoiH.NakataS.NakashimaA.CountsA. L.NandaR. S. (2002b). Changes in tooth position in humans in relation to arterial blood pressure. Arch. Oral. Biol.47, 219226. 10.1016/S0003-9969(01)00110-8

  • 47

    JohnsonR. B.PylypasS. P. (1992). A re-evaluation of the distribution of the elastic meshwork within the periodontal ligament of the mouse. J. Periodont. Res.27, 239249. 10.1111/j.1600-0765.1992.tb01674.x

  • 48

    JonasI. E.RiedeU. N. (1980). Reaction of oxytalan fibers in human periodontium to mechanical stress. A combined histochemical and morphometric analysis. J. Histochem. Cytochem.28, 211216. 10.1177/28.3.7354216

  • 49

    KardosT. B.SimpsonL. O. (1979). A theoretical consideration of the periodontal membrane as a collagenous thixotropic system and its relationship to tooth eruption. J. Periodont. Res.14, 444451. 10.1111/j.1600-0765.1979.tb00243.x

  • 50

    KardosT. B.SimpsonL. O. (1980). A new periodontal membrane biology based on thixotropic concepts. Am. J. Orthod.77, 508515. 10.1016/0002-9416(80)90130-X

  • 51

    KeiligL.StarkH.BayerS.UtzK. H.StrazzaM.GrünnerM.et al. (2009). Numerical investigation of the mechanical loading of supporting soft tissue for partial dentures. Int. J. Prosthodont.22, 201203.

  • 52

    KlymenkoO.Wiltowska-ZuberJ.LekkaM.KwiatekW. M. (2009). Energy dissipation in the AFM elasticity measurements. Acta Phys. Pol. A.115, 548551. 10.12693/APhysPolA.115.548

  • 53

    KohnJ. C.EbensteinD. M. (2013). Eliminating adhesion errors in nanoindentation of compliant polymers and hydrogels. J. Mech. Behav. Biomed. Mater.20, 316326. 10.1016/j.jmbbm.2013.02.002

  • 54

    KomatsuK.ViidikA. (1966). Changes in the fibre arrangement of the rat incisor periodontal ligament in relation to various loading levels in vitro. Arch. Oral Biol.41, 147159. 10.1016/0003-9969(95)00114-X

  • 55

    KristiansenA. B.HeyreaasK. J. (1989). Micropuncture measurements of interstitial fluid pressure in the rat periodontal ligament. Proc. Finn. Dent. Soc.85, 295300.

  • 56

    KurashimaK. (1965). Viscoelastic properties of periodontal tissue. Bull. Tokyo Med. Dent. Univ.36, 148153.

  • 57

    LakesR. S. (1999). Viscoelastic Solids. Boca Raton, FL: CRC Press.

  • 58

    LinD. C.DimitriadisE. K.HorkayF. (2007). Robust strategies for automated AFM force curve analysis-I. Non-adhesive indentation of soft inhomogenous materials. J. Biomech. Eng.129, 430440. 10.1115/1.2720924

  • 59

    LindheJ.KarringT.LangN. P. (2003). Clinical Periodontology and Implant Dentistry. Copenhagen: Blackwell Munksgaard.

  • 60

    MahaffyR. E.ShihC. K.MacKintoshF. C.KäsJ. (2000). Scanning probe-based frequency-dependent microrheology of polymer gels and biological cells. Phys. Rev. Lett.85, 880883. 10.1103/PhysRevLett.85.880

  • 61

    ManningL.RogersB.JonesM.AdamsJ. D.FusteJ. L.MinneS. C. (2003). Self-oscillating tapping mode atomic force microscopy. Rev. Sci. Instrum.74, 42204222. 10.1063/1.1602935

  • 62

    MeyerJ.SquierC. A.GersonS. J. (1984). The Structure and Function of Oral Mucosa. Oxford: Pergamon Press.

  • 63

    MichalakisK. X.CalvaniP.HirayamaH. (2012). Biomechanical considerations on tooth-implant supported fixed partial dentures. J. Dent. Biomech.3:175873601246202510.1177/1758736012462025

  • 64

    MühlemannH. R. (1951). Periodontometry, a method for measuring tooth mobility. Oral. Surg. Oral. Med. Oral. Path4, 12201233. 10.1016/0030-4220(51)90080-1

  • 65

    MurphyM. F.LilleyF.LalorM. J.CrosbyS. R.MaddenG.JohnstonG.et al (2013). Evaluation of a nonlinear Hertzian-based model reveals prostate cancer cells respond differently to force than normal prostate cells. Microsc. Res. Tech.76, 3641. 10.1002/jemt.22132

  • 66

    NanciA. (2013). Ten Cate's Oral Histology. St. Louis, MO: Elsevier Mosby.

  • 67

    NataliA. N.PavanP. G.CarnielE. L.DorowC. (2003). A transversally isotropic elasto-damage constitutive model for the periodontal ligament. Comput. Methods Biomech. Biomed. Engin.6, 329336. 10.1080/10255840310001639840

  • 68

    NataliA. N.PavanP. G.CarnielE. L.DorowC. (2004). Visco-elastic response of the periodontal ligament: an experimental-numerical analysis. Connect. Tissue Res.45, 222230. 10.1080/03008200490885742

  • 69

    NewcombG. M. (1981). An ultrastructural study of epithelial specialization at the porcine mucogingival junction. J. Periodont. Res.16, 5165. 10.1111/j.1600-0765.1981.tb00949.x

  • 70

    NotbohmJ.PoonB.RavichandranG. (2012). Analysis of nanoindentation of soft materials with an atomic force microscope. J. Mater. Res.27, 229237. 10.1557/jmr.2011.252

  • 71

    PictonD. C. (1965). On the part played by the socket on tooth support. Arch. Oral Biol.6, 945955. 10.1016/0003-9969(65)90088-9

  • 72

    PictonD. C. (1969). The effect of external forces on the periodontium, in Biology of the Periodontium, eds MelcherA. H.BowenW. H. (New York, NY: Academic Press), 363419.

  • 73

    RadmacherM.TillmanR. W.FritzM.GaubH. E. (1992). From molecules to cells: imaging soft samples with the atomic force microsope. Science257, 19001905. 10.1126/science.1411505

  • 74

    RahmatB.HubertP. (2010). Interaction stress measurement using atomic force microscopy: a stepwise discretization method. J. Phys. Chem. C.114, 1502915035. 10.1021/jp104993f

  • 75

    RogersB.ManningL.SulchekT.AdamsJ. D. (2004). Improving tapping mode atomic force microscopy with piezoelectric cantilevers. Utramicroscopy100, 267276. 10.1016/j.ultramic.2004.01.016

  • 76

    SaderJ. E.ChonJ. W. M.MulvaneyP. (1999). Calibration of rectangular atomic force microscope cantilevers. Rev. Sci. Instrum.70, 39673969. 10.1063/1.1150021

  • 77

    SaitohS.SasakiK.NezuT.TairaM. (2010). Viscoelastic behavior of commercially available tissue conditioners under compression. Dent. Mater. J.29, 461468. 10.4012/dmj.2009-130

  • 78

    SimsM. (1995). The morphology of the vasculature of the periodontal ligament, in The Periodontal Ligament in Health and Disease, eds BerkovitzB. K. B.MoxhamB. J.NewmanH. N. (London: Mosby-Wolfe), 10720.

  • 79

    SlavkinH. C.BavettaL. A. (1972). Developmental Aspects of Oral Biology.New York, NY: Academic Press.

  • 80

    SloanP.CarterD. H. (1995). Structural organization of the fibres of the periodontal ligament, in The Periodontal Ligament in Health and Disease, eds BerkovitzB. K. B.MoxhamB. J.NewmanH. N. (London: Mosby-Wolfe).

  • 81

    StefanadisC.DernellisJ.TsiamisE.DiamantopoulosL.MichaelidesA.ToutouzasP. (2000). Assessment of aortic line of elasticity using polynomial regression analysis. Circulation101, 18191825. 10.1161/01.CIR.101.15.1819

  • 82

    StolzM.GottardiR.RaitteriR.MiotS.MartinI.ImerR.et al. (2009). Early detection of aging cartilageand osteoarthritis in mice and patient samples using atomic force microscopy. Nat. Nanotechnol.4, 186192. 10.1038/nnano.2008.410

  • 83

    StolzM.RaiteriR.DanielsA. U.VanLandinghamM. R.BaschongW.AebiU. (2004). Dynamic elastic modulus of porcine articular cartilage determined at two different levels of tissue organization by indentation-type atomic force microscopy. Biophys. J.86, 32693283. 10.1016/S0006-3495(04)74375-1

  • 84

    SweigartM. A.ZhuC. F.BurtD. M.DeHollP. D.AgrawalC. M.ClantonT. O.et al. (2004). Intraspecies and interspecies comparison of the compressive properties of the medial meniscus. Ann. Biomed. Eng.32, 15691579. 10.1114/B:ABME.0000049040.70767.5c

  • 85

    TaoN. J.LindsayS. M.LeesS. (1992). Measuring the microelastic properties of biological materials. Biophys. J.63, 11651169. 10.1016/S0006-3495(92)81692-2

  • 86

    ThavarajahR.MudimbaimannarV. K.ElizabethJ.RaoU. K.RanganathamK. (2012). Chemical and physical basics of routine formaldehyde fixation. J. Oral Maxillofac. Pathol.16, 400405. 10.4103/0973-029X.102496

  • 87

    TrickeyW. R.BaaijensF. P.LaursenT. A.AlexopoulosL. G.GuilakF. (2006). Determination of the Poisson's ratio of the cell: recovery properties of chondrocytes after release from complete micropipette aspiration. J. Biomech.39, 7887. 10.1016/j.jbiomech.2004.11.006

  • 88

    Van LoockeM.LyonsC. G.SimmsC. K. (2006). A validated model of passive muscle in compression. J. Biomech.39, 29993009. 10.1016/j.jbiomech.2005.10.016

  • 89

    Van LoockeM.LyonsC. G.SimmsC. K. (2008). Viscoelastic properties of skeletal muscle in compression: stress-relaxation behavior and constitutive modelling. J. Biomech.41, 15551566. 10.1016/j.jbiomech.2008.02.007

  • 90

    VillegasD. F.MaesJ. A.MageeS. D.DonahueT. L. (2007). Failure properties and strain distribution analysis of meniscal attachments. J. Biomech.40, 26552662. 10.1016/j.jbiomech.2007.01.015

  • 91

    VinckierA.SemenzaG. (1998). Measuring elasticity of biological materials by atomic force microscopy. FEBS Lett.430, 1216. 10.1016/S0014-5793(98)00592-4

  • 92

    WeisenhornA. L.KasasS.SollettiJ. M.KhorsandiM.GotzosV.RomerD. U.et al. (1993a). Deformation observed on soft surfaces with an AFM. Proc. SPIE18552634. 10.1117/12.146382

  • 93

    WeisenhornA. L.KhorsandiM.KasasS.GotzosV.ButtH. J. (1993b). Deformation and height anomaly of soft surfaces studied with an AFM. Nanotechnology4, 106113. 10.1088/0957-4484/4/2/006

  • 94

    WeisenhornA. L.MaivaldP.ButtH. J.HansmaP. K. (1992). Measuring adhesion, attraction, and repulsion between surfaces in liquids with an atomic-force microscope. Phys. Rev. B. Condensed. Matt.45, 1122611232. 10.1103/PhysRevB.45.11226

  • 95

    WeissJ. A.GardinerJ. C.Bonifasi-ListaC. (2002). Ligament material behavior is nonlinear, viscoelastic and rate-independent under shear loading. J. Biomech.35, 943950. 10.1016/S0021-9290(02)00041-6

  • 96

    WillsD. J.PictonD. C. (1978). Changes in the mobility and resting position of incisor teeth in macaque monkeys. Arch. Oral. Biol.23, 225229. 10.1016/0003-9969(78)90221-2

  • 97

    WillsD. J.PictonD. C.DaviesW. I. R. (1976). A study of the fluid systems of the periodontium in macaque monkeys. Arch. Oral. Biol.21, 175185. 10.1016/0003-9969(76)90127-8

  • 98

    YangF. (1998). Indentation of an incompressible elastic film. Mech. Mater.30, 275286. 10.1016/S0167-6636(98)00035-0

  • 99

    YangF. (2006). Effect of adhesion energy on the contact stiffness in nanoindentation. J. Mater. Res.21, 26832688. 10.1557/jmr.2006.0331

  • 100

    YangS.ZhangY. W.ZengK. Y. (2004). Analysis of nanoidentation creep for polymeric materials. J. Appl. Phys.95, 36553666. 10.1063/1.1651341

  • 101

    YinL.ElliotD. M. (2004). A biphasic and transversely isotropic mechanical model for tendon: application to mouse tail fascicles in uniaxial tension. J. Biomech.37, 907916. 10.1016/j.jbiomech.2003.10.007

  • 102

    ZhangY. L.HanM. L.VidyalakshmiJ.SheeC. Y.AngW. T. (2009). Automatic control of mechanical forces acting on cell biomembranes using a vision-guided microrobotic system in computer microscopy. J. Microsc.236, 7078. 10.1111/j.1365-2818.2009.03209.x

Summary

Keywords

human oral keratinized mucosa, atomic force microscopy, oral tissue mechanics, contact mechanics, polynomial model, mathematical model, curve fitting

Citation

Tsaira A, Karagiannidis P, Sidira M, Kassavetis S, Kugiumtzis D, Logothetidis S, Naka O, Pissiotis A and Michalakis K (2016) Theoretical Considerations and a Mathematical Model for the Analysis of the Biomechanical Response of Human Keratinized Oral Mucosa. Front. Physiol. 7:364. doi: 10.3389/fphys.2016.00364

Received

28 May 2016

Accepted

09 August 2016

Published

29 August 2016

Volume

7 - 2016

Edited by

Thimios Mitsiadis, University of Zurich, Switzerland

Reviewed by

Jean-Christophe Farges, Claude Bernard University Lyon 1, France; Christian Hellmich, Vienna University of Technology, Austria; Nenad Filipovic, University of Kragujevac, Serbia

Updates

Copyright

*Correspondence: Konstantinos Michalakis ; ;

This article was submitted to Craniofacial Biology, a section of the journal Frontiers in Physiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics