Stem Cells International
Volume 2018 (2018), Article ID 3208285, 17 pages
https://doi.org/10.1155/2018/3208285

Review Article

In Vitro Weight-Loaded Cell Models for Understanding Mechanodependent Molecular Pathways Involved in Orthodontic Tooth Movement: A Systematic Review

Mila Janjic,¹ Denitsa Docheva,² Olivera Trickovic Janjic,³ Andrea Wichelhaus,¹ and Uwe Baumert¹

¹Department of Orthodontics and Dentofacial Orthopedics, University Hospital, LMU Munich, 80336 Munich, Germany
²Experimental Trauma Surgery, Department of Trauma Surgery, University Regensburg Medical Centre, 93053 Regensburg, Germany
³Department of Preventive and Pediatric Dentistry, Faculty of Medicine, University of Niš, 18000 Niš, Serbia

Correspondence should be addressed to Uwe Baumert; uwe.baumert@med.uni-muenchen.de

Received 26 January 2018; Accepted 9 May 2018; Published 31 July 2018

Academic Editor: Andrea Ballini

Copyright © 2018 Mila Janjic et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Cells from the mesenchymal lineage in the dental area, including but not limited to PDL fibroblasts, osteoblasts, and dental stem cells, are exposed to mechanical stress in physiological (e.g., chewing) and nonphysiological/therapeutic (e.g., orthodontic tooth movement) situations. Close and complex interaction of these different cell types results in the physiological and nonphysiological adaptation of these tissues to mechanical stress. Currently, different in vitro loading models are used to investigate the effect of different types of mechanical loading on the stress adaptation of these cell types. We performed a systematic review according to the PRISMA guidelines to identify all studies in the field of dentistry with focus on mechanobiology using in vitro loading models applying uniaxial static compressive force. Only studies reporting on cells from the mesenchymal lineage were considered for inclusion. The results are summarized regarding gene expression in relation to force duration and magnitude, and the most significant signaling pathways they take part in are identified using protein-protein interaction networks.

1. Introduction

The aim of orthodontics is to move an abnormally positioned tooth through the application of a continuous force on its surface. This force stimulates bone remodelling in the surrounding tissue, namely, the periodontal ligament (PDL) and the alveolar bone, resulting in the bone removal in the direction of the tooth movement and bone apposition in the opposite direction (Figure 1). Thus, the underlying mechanism of orthodontic tooth movement (OTM) is the stimulation of bone remodelling by the application of an orthodontic force [1].

(a)

(b)

Figure 1: Bone remodelling during orthodontic tooth movement. (a) Initial displacement of the tooth due to stretching of the fibres within the PDL on the tension side and compression on the opposite with the application of the orthodontic force. (b) Bone apposition on the tension side and resorption on the compression side as the result of the long-term force application.

Histologically, the effects of orthodontic force on the tooth and its surrounding tissues are now well understood and the underlying stages in OTM are identified [2]. Human periodontal ligament cells (hPDLCs) and human osteoblasts (hOBs) are recognized as the cell types originating from the mesenchymal lineage, which play the most dominant role during OTM. Unlike hOBs, which represent well a characterized cell type, hPDLCs represent a mixed population of mostly fibroblast-like cells [3]. Among them, mesenchymal stem cells are of special importance as the source of progenitors responsible for the regeneration and remodulation of not only PDL itself but also alveolar bone [4].

In order to better understand morphological changes during OTM, it is important to elucidate molecular and cellular signaling mechanisms between and within these cell types. The complex in vivo structure of the tissues involved makes it impossible to investigate force sensing and cellular communication of individual cells. Therefore, in vitro models using cells isolated from the PDL or from alveolar bone were established and different types of forces mimicking those found during OTM were applied [5]. These in vitro models are used to answer open questions including but not limited to how cells sense force, how they convert mechanical stress into molecular signals, and how these molecular signals influence the specific response of these cells to that specific force.

On the basis of the most commonly used approaches to apply mechanical stress on cells, present in vitro loading models can be classified into those using substrate deformation-based approaches, hydrostatic pressure approach, centrifugation approach, fluid flow approach, vibration approach, and weight approach [6]. Also, there has been increasing interest in moving from conventional monolayer, two-dimensional (2D) in vitro loading models to three-dimensional (3D) in vitro loading models.

Weight-based in vitro loading models have been successfully used over several years to investigate the effect of static, compressive, unidirectional force on the cells. In models using 2D cell cultures, cells are precultured in cell culture dishes (e.g., 6-well plates). After reaching the desired confluency, the cells are subjected to weight-based compression. In most cases, a glass slide is laid on top of the cell monolayer. Then, a weight is applied by positioning a glass cylinder filled with lead granules on top of this slide. The glass slide is used to secure even distribution of the force [7]. Increasing or reducing the number of granules in the glass cylinder adjusts the level of compressive force (Figure 2(a)). The same type of force is applied by slight modifications of this model: some authors used a stack of glass slides of different heights (e.g., [8]) or glass discs of different thicknesses (e.g., [9]) replacing the glass cylinder filled with lead granules. This in vitro loading model can also be used to apply static compressive force on 3D cell cultures. In this case, the same principle is used, except that the cells are embedded in a 3D matrix that is then compressed in the described manner (Figure 2(b)). Yang et al. [6] coined the term “weight approach”-based (WAB) for this in vitro model. To refer to this specific setup, we will also use WAB throughout this publication.

(a)

(b)

Figure 2: Schematic illustration of the static 2D (a) and 3D (b) in vitro loading model based on the weight approach applied in the literature (details are found in the text).

The primary aim of this review was to identify all articles related to the field of orthodontics using either a 2D or 3D WAB in vitro loading model and provide an overview of the details of their use: the most commonly used loading durations, force magnitudes, and scaffolds and their findings regarding gene expression and substance secretion in relation to force application. The secondary objective was to discover most commonly examined genes and to identify important pathways in OTM that most of the identified genes from these studies are involved in, focusing especially on hPDLCs.

2. Materials and Methods

To conduct this review, the “Preferred Reporting Items for Systematic Review and Meta-Analysis Protocols” (PRISMA-P) 2015 statement was consulted [10].

2.1. Defining the Eligibility Criteria

Inclusion criteria were as follows:

(i)Studies in the field of dentistry that examined the effect of mechanical stress on tooth surrounding tissues(ii)Application of the 2D or 3D WAB in vitro loading model…(iii)…on hPDLCs, hOBs, or all bone-like cell types/lines of human or animal origin(iv)Only studies written in English language, identified on the PubMed database until 01.12.2017, were taken into consideration

2.2. Literature Search and Study Selection Process

Separate search strategies were created for studies using either the 2D or the 3D in vitro setup for mechanical cell loading (Supplement 1). Searches were performed in the PubMed database following these predefined search strategies.

After identification of relevant studies in the PubMed database, the downloaded records from each search were imported into the bibliographic software EndNote X8 (Clarivate Analytics, Philadelphia, Pennsylvania, USA). All records were examined by two reviewers independently (MJ and UB), according to predefined inclusion and exclusion criteria (see above): first by title, then by abstract. If the abstract was not available, the full text of the report was obtained. Records that were obviously irrelevant were excluded, and the full texts of all remaining records were acquired. After the full-text assessment, the final list of included articles was generated. Any disagreements during this process were dissolved through discussion with another review author (DD) until reaching a consensus. The articles that did not meet all inclusion criteria after full-text assessment were excluded from further examination. Additional relevant studies were further identified through forward and backward reference chaining and hand-search of specific journals. Study quality assessment of the included studies was not performed, since the goal of this article was to provide an overview of all findings in the field only.

2.3. Data Extraction

The following information was extracted from each study obtained in full length: author, journal, year of publication, and used cell type. Force magnitude and duration, examined genes or substances, gene expression, or substance secretion details were recorded only if their response was directly connected to mechanical force stimulus. Gene symbols were used in the tables whenever possible. In case the identity or variant of a gene was doubtful or not clear primer sequences were examined using Primer-BLAST (URL: https://www.ncbi.nlm.nih.gov/tools/primer-blast/) [11]. If Western blot, ELISA, or inhibition experiments were reported, we tried to verify the antibodies and/or inhibitor specificity to determine the exact protein species (variant). Additionally, the method used for evaluation of the gene/substance expression was recorded. Data regarding the used scaffolds were collected for studies applying 3D WAB in vitro setups.

The following tables were prepared to summarize the findings: (1) studies applying the 2D WAB in vitro loading model on human primary cells from the orofacial region (i.e., hPDLCs, hOBs, and human oral bone marrow cells), (2) studies applying the 2D WAB in vitro loading model on human and nonhuman cells and cell lines not included in the first table, and (3) studies applying the 3D WAB in vitro loading model on human and nonhuman cells and cell lines.

2.4. STRING Analysis

The examined genes and metabolites using the 2D approach were summarized in two separate lists: one for hPDLFs and one for hOBs and other human bone-derived cell lines. Protein-protein interaction (PPI) networks were generated for both lists separately using the STRING database (10.5, URL: https://string-db.org/) [12]. From within STRING, the KEGG database [13] was queried to identify the main pathways involved. Only pathways with a false discovery rate below 1.00E−05 were considered.

3. Results

3.1. Study Selection Process

Figure 3 summarises the results of both 2D and 3D searches using a flow chart according to PRISMA. Separate searches were conducted for the studies applying either the 2D or 3D (Supplement 1) WAB in vitro loading models.

Figure 3: PRISMA flow diagram of the review process.

The search formula applied to identify 2D WAB in vitro loading studies is shown in Supplement 1. Altogether, 2284 abstracts were identified in the PubMed database (Figure 3).

Additionally, 7 articles were identified through forward and backward reference chaining and hand-search of specific journals. After reading the titles and abstracts of all identified studies, we excluded 2184. The remaining 107 articles were then checked by full-text reading. Fifty-six of them meet our inclusion criteria and were included for further analysis. The remaining did not meet the inclusion criteria. Reasons for their exclusion are listed in Supplement 1.

The search formula applied to identify 3D WAB in vitro loading studies is shown in Supplement 1. We identified a total of 1038 articles in PubMed (Figure 3). Additional 4 articles were discovered through forward and backward reference chaining and hand-search of specific journals. After initial screening, we excluded 992 articles and proceeded with full-text reading of the 50 articles. Finally, 17 of them meet our inclusion criteria. The remaining articles were excluded from further analysis. Reasons for their exclusion are summarized in Supplement 1.

All studies fulfilling the inclusion criteria were organised into three different supplementary tables: Supplement 2 summarises 2D WAB in vitro loading studies using human primary cells from the orofacial region. In Supplement 3, the two-dimensional WAB in vitro loading studies using human nonorofacial-derived cells and animal cells and cell lines are found. Supplement 4 summarises the 3D WAB in vitro loading studies.

3.2. Force Durations and Force Magnitudes Used in the Studies

3.2.1. 2D WAB In Vitro Loading Model

In these studies, compression forces ranging from 0.25 g/cm² to 5 g/cm² were applied on cells in 2D culture. The most commonly used compressive force was 2 g/cm², irrespectively which cell type was used in the study. In most of the studies, the force was applied for 24 h (Supplements 2 and 3).

3.2.2. 3D WAB In Vitro Loading Model

Force duration and magnitude depended on the scaffold used (Supplement 4). In most of the studies, scaffolds made from collagen gel and the polylactic-co-glycolic acid (PLGA) were applied. One of the studies [14] used a hydrophilically modified poly-L-lactide (PLLA) matrix. Collagen gel scaffolds were used with force magnitudes varying between 0.5 g/cm² and 9.5 g/cm²; the most commonly used force was 6 g/cm². Force was applied for 0.5 to 72 h. Most commonly used force application periods were 12 and 24 h. Force levels between 5 and 35 g/cm² were applied to cells embedded in PLGA scaffolds. The most commonly applied force was 25 g/cm². The duration of force application was from 3 to 72 h. The study using the hydrophilically modulated PLLA matrix [14] applied force magnitudes from 5 to 35 g/cm². The duration of force application varied between one day and 14 days.

3.3. Cell Types Used in the Studies

3.3.1. 2D WAB In Vitro Loading Model

Forty of these studies used human primary cells isolated from the tooth surrounding tissues (Supplement 2): hPDLCs, hOBs, and human orofacial bone marrow-derived cells (hOBMC). The remaining studies used other cells and cell lines from human and animal sources: MG63, RAW264.7, ST-2, Saos-2, OCCM-30, MC3T3-E1, C2C12, U2OS, rat-derived PDLCs, or bone marrow-derived osteoblasts and the cementoblast cell line HCEM-SV40 (Supplement 3).

3.3.2. 3D WAB In Vitro Loading Model

hPDLCs and human gingival fibroblasts were used in 13 studies (Supplement 4). The remaining two studies used cell types and lines from the nonoral region or nonhuman origin (Supplement 4): the murine cell line MC3T3-E1 and murine osteoblasts.

Taken together, the most commonly used cells were hPDLCs. They were used in total 51 studies (2D: 38; 3D: 13) (Supplements 2 and 4). According to the isolation method applied, we distinguished between the following sources: “explant method” [15, 16] (2D: 18; 3D: 4), “enzyme digestion method” [4] (2D: 9; 3D: 6), commercial sources (2D: 3; 3D: 1), or “no detailed information of isolation available” (2D: 8; 3D: 2).

3.4. Genes and Substances Examined in the Studies

A complete overview of genes and metabolites examined in 2D and 3D WAB studies and details of their expression can be found in Supplements 2 and 3 (2D) and Supplement 4 (3D).

In this review, special attention was paid to hPDLCs as the most examined cell type among studies and their prominent role in OTM. The most examined genes and metabolites in relation to hPDLCs were TNF superfamily member 11 (TNFSF11), TNF receptor superfamily member 11B (TNFRSF11B), prostaglandin-endoperoxide synthase 2 (PTGS2), and prostaglandin E₂ (PGE₂). In Table 1, details regarding their expression/secretion, including the information at which time points or force magnitudes the highest/lowest value was reached, is summarized.

Table 1: Top four examined genes or substances in studies applying 2D or 3D in vitro WAB loading model on hPDLCs. For each gene or substance, cell culture type and literature reference are given. Additionally, examined force durations and magnitudes are summarized. Force effect on gene expression/substance secretion was evaluated (increase/decrease/no change). In all cases, the most prominent changes (increase/decrease) or “no change” are noted. For each change of expression/secretion, the corresponding maximum (increase/decrease) of force duration and magnitudes are additionally provided.

Gene symbol or metabolite
Cell culture
Reference
Examined force applied
Gene expression
Substance secretion

Duration (h)
Magnitude (g/cm²)
Increase/decrease/no change
Change in relation to force duration (h)
Change in relation to force magnitude (g/cm²)
Increase/decrease/no change
Change in relation to force duration (h)
Change in relation to force magnitude (g/cm²)

PGE₂
2D
Benjakul et al. in press [17]
48
1.5
na

Increase (qPCR: GAPDH)
48
1.5

Jin et al. 2015 [18]
0; 0.5; 3; 6; 12
2.0
na

Increase (ELISA)
12
2.0

Kang et al. 2010 [19]
0.5; 2; 6; 24; 48
2.0
na

Increase (ELISA)
48
2.0

Kanzaki et al. 2002 [20]
0.5; 1.5; 6; 24; 48 (+ELISA: 60)
0.5; 1.0; 2.0; 3.0; 4.0 (ELISA: 2.0)
na

Increase (ELISA)
60
2.0

Kirschneck et al. 2015 [21]
24
2.0
na

Not explicitly stated (ELISA)

Liu et al. 2006 [22]
48
2.0
na

Increase (ELISA)
48
2.0

Mayahara et al. 2007 [23]
3; 6; 12; 24; 48
2
na

Increase (ELISA)
48
2

Premaraj et al. 2013 [24]
0.5; 1; 3; 6
5.0
na

Increase (ELISA)
1
5.0

Proff et al. 2014 [9]
24
2
na

Increase (ELISA)
24
2

Römer et al. 2013 [25]
24
2
na

Increase (ELISA)
24
2

3D (Coll. gel)
de Araujo et al. 2007 [26]
3; 12; 24; 48; 72
6.0

Increase (EIA)
72
6.0

3D (PLGA)
Li et al. 2016 [27]
6; 24; 72
5.0; 15.0; 25.0
na

Increase (ELISA)
24
15.0…25.0

Yi et al. 2016 [28]
24
25.0

Increase (ELISA)
24
25.0

PTGS2
2D
Jin et al. 2015 [18]
0; 0.5; 3; 6; 12
2.0
Increase (qPCR: GAPDH)
12
2.0

Kang et al. 2010 [19]
0.5; 2; 6; 24; 48
2.0
Increase (qPCR: GAPDH)
48

Kanzaki et al. 2002 [20]
0.5; 1.5; 6; 24; 48
0.5; 1.0; 2.0; 3.0; 4.0
Increase (sqPCR: ACTNB)
6
2.0

Kirschneck et al. 2015 [21]
24
2.0
Increase (qPCR: POL2RA)
24
2.0

Liu et al. 2006 [22]
48
2.0
Increase (sqPCR: ACTNB)
48
2.0

Mayahara et al. 2007 [23]
3; 6; 12; 24; 48
2
Increase (qPCR: GAPDH)
48
2

Mayahara et al. 2010 [29]
3; 6; 12; 24; 48
2.0
Increase (qPCR: GAPDH)
48
2

Premaraj et al. 2013 [24]
6
0.2; 2.2; 5.0
nd

Increase (WB)
6
5.0

Proff et al. 2014 [9]
24
2
Increase (qPCR: POL2RA)
24
2
Increase (WB)
24
2

Römer et al. 2013 [25]
24
2
Increase (qPCR: POL2RA)
24
2

Wongkhantee et al. 2007 [30]
24
0; 1.25; 2.5
Increase (sqPCR: GAPDH)
24
2.5

3D (Coll. gel)
de Araujo et al. 2007 [26]
1; 3; 6; 12; 24; 48; 72
3.6; 6.0; 7.1; 9.5
Increase (sqPCR: GAPDH)
6
7.1

3D (PLGA)
Li et al. 2016 [31]
6; 24; 72
25.0
Increase (qPCR: GAPDH)
6
25.0

Li et al. 2013 [32]
6; 24; 72
25.0
Increase (qPCR: GAPDH)
6
25.0

Li et al. 2016 [27]
6; 24; 72
5.0; 15.0; 25.0
Increase (qPCR: GAPDH)
6
25.0

Li et al. 2011 [33]
6
5; 15; 25; 35
Increase (qPCR: GAPDH)
6
35.0

Yi et al. 2016 [28]
24
25.0
Increase (qPCR: GAPDH)
24
25.0
Increase (WB)
24
25.0

TNFRSF11B
2D
Benjakul et al. in press [17]
48
1.5
No change (qPCR: GAPDH)

No change

Jin et al. 2015 [18]
0; 0.5; 3; 6; 12
2.0
No change (qPCR: GAPDH)

Kanzaki et al. 2002 [20]
0.5; 1.5; 6; 24; 48
0.5; 1.0; 2.0; 3.0; 4.0
No change (sqPCR: ACTNB)

Kim et al. 2013 [8]
0.5; 2; 6; 24; 48
2.0
Transitory downregulated. (qPCR: GAPDH)
6
2.0
Transitory downregulation (ELISA)
6
2.0

Kirschneck et al. 2015 [21]
24
2.0
No change (qPCR: POL2RA)

Lee et al. 2015 [34]
0; 2; 4; 8; 24; 48
2.5
No change (qPCR: ACTNB)

Liu et al. 2017 [35]
6; 12; 24
0.5; 1.0; 1.5
nd

Decrease (WB)
n. g.
1.5

Luckprom et al. 2011 [36]
2; 4
2.5
No change (sqPCR: GAPDH)

Mitsuhashi et al. 2011 [37]
1; 3; 6; 9; 12; 24
4.0
No change (qPCR: ACTNB)

Nakajima et al. 2008 [38]
0; 1; 3; 6; 9; 12; 24
0.5; 1.0; 2.0; 3.0; 4.0
nd

Increase (ELISA)
24
0.5

Nishijima et al. 2006 [39]
48
0; 0.5; 1.0; 2.0; 3.0
nd

Decrease (ELISA)
48
2.0

Römer et al. 2013 [25]
24
2
No change (qPCR: RNA-polymerase-2-polypeptide A)

Yamada et al. 2013 [40]
12
4.0
Decrease (qPCR: GAPDH)
12
4.0
Decrease (ELISA)
12
4.0

Yamaguchi et al. 2006 [41]
0; 3; 6; 9; 12; 24; 48
0.5; 1.0; 2.0; 3.0
n. d.

Decrease (ELISA)
12…48
2.0

3D (Coll. gel)
Kaku et al. 2016 [42]
12; 24
0.5; 1.0; 2.0
Increase (qPCR: GAPDH)
12
1.0

3D (PLLA modif.)
Liao et al. 2016 [14]
1 d; 3 d; 7 d; 14 d
5.0; 15.0; 25.0; 35.0
No change (qPCR: GAPDH)

3D (PLGA)
Jianru et al. 2015 [43]
3; 6; 12 (WB: 12)
25.0
Decrease followed by increase (qPCR: GAPDH)
3 (decrease)
12 (increase)
25.0
Increase (WB)
12
25.0

Li et al. 2016 [31]
6; 24; 72
25.0
Decrease followed by Increase (qPCR: GAPDH)
6 (decrease)
72 (increase)
25.0

Li et al. 2016 [27]
6; 24; 72
5.0; 15.0; 25.0
Decrease followed by increase (qPCR: GAPDH)
6 (decrease)
72 (increase)
15.0 (decrease)
25.0 (increase)
Decrease followed by Increase (qPCR: GAPDH)
6 (decrease)
72 (increase)
25.0 (decrease)
25.0 (increase)

Li et al. 2011 [33]
6; 24; 72
25
Decrease followed by increase (qPCR: GAPDH)
6 (decrease)
72 (increase)
25.0

Yi et al. 2016 [28]
24
25.0
Decrease (qPCR: GAPDH)
24
25.0
No change (WB)

TNFSF11
2D
Benjakul et al. in press [17]
48
1.5
Increase (qPCR: GAPDH)
48
1.5
Increase (qPCR: GAPDH)
48
1.5

Jin et al. 2015 [18]
0; 0.5; 3; 6; 12
2.0
Increase (qPCR: GAPDH)
12
2.0

Kang et al. 2013 [44]
2; 48
2.0
Increase (qPCR: GAPDH)
48
2.0

Kanzaki et al. 2002 [20]
0.5; 1.5; 6; 24; 48
0.5; 1.0; 2.0; 3.0; 4.0
Increase (sqPCR: ACTNB)
48
2.0
Increase (WB): 40-kDa+ 55-kDa
48
2.0

Kikuta et al. 2015 [45]
1; 3; 6; 9; 12; 24 (+ELISA: 48)
4.0
Increase (qPCR: GAPDH)
12
4.0
Increase (ELISA)
24
4.0

Kim et al. 2013 [8]
0.5; 2; 6; 24; 48
2.0 ++
Increase (qPCR: GAPDH)
24
2.0
Increase (ELISA)
48
2.0

Kirschneck et al. 2015 [21]
24
2.0
Increase (qPCR: POL2RA)
24
2.0

Lee et al. 2015 [34]
0; 2; 4; 8; 24; 48
2.5
Increase (qPCR: ACTNB)
24
2.5

Liu et al. 2017 [35]
6, 12, 24
0.5; 1.0; 1.5
nd

Increase (WB: GAPDH)
ng
1.5

Liu et al. 2006 [22]
48
2.0
Increase (sqPCR: ACTNB)
48
2.0

Luckprom et al. 2011 [36]
2; 4
2.5
Increase (sqPCR: GAPDH)
2
2.5
Increase (WB)
4
2.5

Mitsuhashi et al. 2011 [37]
1; 3; 6; 9; 12; 24
4.0
Temporary increase (qPCR: ACTNB)
6…9
4.0

Nakajima et al. 2008 [38]
0; 1; 3; 6; 9; 12; 24
0.5; 1.0; 2.0; 3.0; 4.0
nd

Increase (ELISA)
24
4.0

Nishijima et al. 2006 [39]
48
0; 0.5; 1.0; 2.0; 3.0
nd

Increase (ELISA)
12…48
2.0

Römer et al. 2013 [25]
24
2
Increase (qPCR: RNA-polymerase-2-polypeptide A)
24
2

Wongkhantee et al. 2007 [30]
24
0; 1.25; 2.5
Increase (sqPCR: GAPDH)
24
2.5
Increase (WB; ACTNB)
24
2.5

Yamada et al. 2013 [40]
12
4.0
Increase (qPCR: GAPDH)
12
4.0
Increase (ELISA)
12
4.0

Yamaguchi et al. 2006 [41]
0; 3; 6; 9; 12; 24; 48
0.5; 1.0; 2.0; 3.0
nd

Increase (ELISA): sRANKL
Increase (WB)
12…48
12
2.0
2.0

3D (Coll. gel)
Kang et al. 2013 [44]
2; 48
2.0
Increase (qPCR: GAPDH)
2
2.0

3D (PLLA modif.)
Liao et al. 2016 [14]
1 d; 3 d; 7 d; 14 d
5.0; 15.0; 25.0; 35.0
Increase (qPCR: GAPDH)
Day 14
35.0

3D (PLGA)
Jianru et al. 2015 [43]
3; 6; 12 (WB: 12)
25.0
Increase (qPCR: GAPDH)
6
25.0
Increase (WB)
12
25.0

Li et al. 2016 [31]
6; 24; 72
25.0
Increase (qPCR: GAPDH)
6
25.0

Li et al. 2016 [27]
6; 24; 72
5.0; 15.0; 25.0
Increase (qPCR: GAPDH)
6
25.0
Decrease (ELISA)
72
25.0

Li et al. 2011 [33]
6; 24; 72
5; 15; 25; 35
Increase (qPCR: GAPDH)
Increase followed by no change (qPCR: GAPDH)
6
6 (increase)
72 (no change)
25…35.0
25
25

Yi et al. 2016 [28]
24
25.0
Increase (qPCR: GAPDH)
24
25.0
Increase (WB)
24
25.0

2D: two-dimensional cell culture; 3D (Coll. gel): three-dimensional cell culture, collagen gel; 3D (PLGA): three-dimensional cell culture using PLGA scaffolds; 3D (PLLA modif.): three-dimensional cell culture, hydrophilically modified PLLA scaffolds; qPCR: quantitative polymerase chain reaction (e.g., real-time PCR); sqPCR: semiquantitative polymerase chain reaction, followed by reference gene used; nr: not reported; na: not applicable; ELISA: enzyme-linked immune absorbent assay; WB: Western blot; IF: immunofluorescence; FLM: fluorescence microscopy; EIA: enzyme immunoassay.

3.5. STRING Analysis and KEGG Pathways

3.5.1. Construction of Protein-Protein Interaction (PPI) Network

In order to elucidate the molecular mechanisms of OTM and the role of the hPDLCs and bone cells in this process, we used STRING to construct PPI networks. Two separate gene lists were compiled from those studies using hPDLCs (“hPDLC list”; data from Supplement 3) and from those using hOBs or human bone-cells and cell lines (“hOB list”; data from Supplements 2 and 3). The hPDLC list contained 48 different genes (Figure 4(a)) and the hOB list 51 different genes (Figure 4(b)).

(a)

(b)

Figure 4: Protein-protein interaction networks for the (a) “hPDLC list” and the (b) “hOB list”. The gene lists are shown in the lower left part of each subfigure. Those genes with the highest number of interactions (“top 10”) are given in tables in the lower right part of each subfigure.

Two separate PPI networks were obtained, based on the interactions with a high level of confidence (>0.700) (Figure 4). Nodes in the networks represent the proteins produced by a single protein-coding gene locus; edges represent protein-protein interaction. Based on the colour of the edge, eight different interactions based on “gene neighbourhood,” “gene fusion,” “cooccurrence,” “coexpression,” “experiments,” “databases,” and “text mining” can be differentiated [12]. The top 10 nodes with the highest degree of connections from each of the two gene lists are also shown in Figure 4. PPI enrichment values for each constructed network were calculated in STRING. These show that both PPI networks had significantly more interactions than expected and that the nodes are not random (PP enrichment value < 1.0E–16).

3.5.2. Identification of KEGG Pathways

According to our STRING analysis, KEGG pathways relevant for OTM for each set of genes are listed in Table 2.

Table 2: KEGG pathways relevant for OTM with false discovery rates below 1.00E − 05 derived from STRING analysis using the set of examined genes from human periodontal ligament cells (“hPDLC list”; top panel) and human bone and bone-related cells and cell lines (“hOB list”; bottom panel). “X”, gene involved in that specific pathway.

(a)

KEGG ID
4060
4668
4510
4620
4370
4062
4380
4010
4064

KEGG name
Cytokine-cytokine receptor interaction
TNF signaling pathway
Focal adhesion
Toll-like receptor signaling pathway
VEGF signaling pathway
Chemokine signaling pathway
Osteoclast differentiation
MAPK signaling pathway
NF-kappa B signaling pathway

False discovery rate
2.62E–15
2.06E–12
3.90E–11
2.04E–09
9.47E–08
1.33E–07
2.29E–07
1.42E–06
1.86E–05

ADRB2

AKT1

X
X
X
X
X
X
X

ALPL

BGLAP

CBS

CCL2
X
X

X

CCL3
X

X

X

CCL5
X
X

X

X

CCND1

X

CCR5
X

X

CDH11

COL1A1

X

COL3A1

X

COL5A1

X

CSF1
X
X

X

CTNNB1

X

CTSB

CTSL

CXCL8 (= IL8)
X

X

X

X

FGF2

X

GJA1

GSK3b

X

X

HMGB1

HSP90AA1

HSPA4

HSPB1

X

X

IGF1

X

IL17A
X

IL1B
X
X

X

X
X
X

IL6
X
X

X

JAG1

X

LGALS3BP

MMP13

MMP3

X

PIEZO1

PLA2G4A

X

X

POSTN

PTGS1

PTGS2

X

X

X

PTK2

X

X
X

RUNX2

SPP1

X
X

TGFB1
X

X

TGFB3
X

X
X

TNF
X
X

X

X
X
X

TNFRSF11B
X

X

TNFSF11
X

X

X

VEGFA
X

X

X

(b)

KEGG ID
4350
4060
4064
4390
4668
4210
4380
4620
4066

KEGG name
TGF-beta signaling pathway
Cytokine-cytokine receptor interaction
NF-kappa B signaling pathway
Hippo signaling pathway
TNF signaling pathway
Apoptosis
Osteoclast differentiation
Toll-like receptor signaling pathway
HIF-1 signaling pathway

False discovery rate
8.33E–23
2.37E–21
8.32E–11
5.07E–09
1.01E–08
6.26E–08
1.02E–05
6.79E–05
7.16E–05

ACVR1
X
X

ACVR2A
X
X

ACVR2B
X
X

ALPL

BAX

X

BCL2

X

X

X

BGLAP

BMP2
X
X

X

BMP4
X

X

BMP6
X

X

BMP7
X
X

X

BMPR1A
X
X

X

BMPR1B
X
X

X

BMPR2
X
X

X

Casp3

X
X

CHRD
X

CXCR1

X

FST
X

GREM1

IBSP

IL11

X

IL11RA

IL1b

X
X

X
X
X
X

IL1r1

X
X

X
X

IL6

X

X
X
X

IL6R

X

X

IL8

X
X

X
X

MKI67

MMP1

MMP13

MMP14

MMP2

MMP3

NOG
X

PLAT

PLAU

X

PTGS2

X

RUNX2

SERPINE1

X

X

SMAD1
X

X

SP7

SPP1

X
X

TIMP1

X

TIMP2

TIMP3

TIMP4

TNF
X
X
X

X
X
X
X

TNFRSF11B

X

TNFRSF1A

X
X

X
X

TNFSF11

X
X

ZNF354C

4. Discussion

In vivo bone remodelling during OTM represents a complex biological process, triggered by mechanical stimuli. OTM involves numerous events, spatially and temporary orchestrated and coordinated by different cell types, signaling factors, and networks [1]. Systematic breakdown and analysis of individual components of this complex process is the key for understanding its molecular background and a possible way to accelerate and improve it. Therefore, a variety of in vitro mechanical loading models have been established [5, 6]. The in vitro loading model based on the weight approach has been considered as the most appropriate loading model for the stimulation of the orthodontic force on the compressive site [6].

4.1. Characteristics of 2D and 3D WAB In Vitro Loading Models

4.1.1. Conventional 2D WAB

In vitro loading model, initially described by Kanai et al. [7], has been used for more than two decades for studying the compression-induced osteoclastogenesis and is still considered as the gold standard. It represents a simple and effective method for application of static compressive, unidirectional force to a cell monolayer.

The advantages of WAB in vitro loading model are the following:

(i)It reduces the need for animal studies, which are costly and time consuming.(ii)It enables the analysis of specific cell types independently or in cocultures with other cells of interest.(iii)Human primary cells can be used for better approximation to clinical situation.

From our point of view, the main disadvantage is its missing impact of the natural surrounding environment. There has been an increasing interest in the development of the 3D cell culture WAB in vitro loading model during the last years, in order to approximate the in vitro situation to the in vivo situation.

4.1.2. 3D WAB In Vitro Loading Model

During the last years, more studies have been using cells incorporated into biological scaffolds instead of monolayer cultures. This is due to the demand of mimicking an extracellular matrix, which is beneficial for cell behaviour, instead of growing cells on artificial plastic cell culture surface [46]. According to our data, three types of scaffolds have been used so far in combination with the 3D WAB in vitro loading model. The first identified studies used collagen I scaffolds [26, 47, 48]. Although the collagen gels are still widely used for this purpose, there is the increasing interest in the development of scaffolds composed of synthetic polymers. In 2011, Li et al. [33] introduced the PLGA scaffolds that had a higher stiffness in comparison to collagen gels and an elastic modulus very close to that of human PDL. The only disadvantage was that cells growing in PLGA displayed a disordered grow pattern that differs from the one in natural PDL [33]. Liao et al. [14] went one step further and introduced a hydrophilically modified PLLA matrix. This matrix displayed several advantages: higher nutrient and oxygen permeability and a better cell attachment, making it more suitable for long-term force application [14].

4.2. Force Magnitude Used in the Studies

According to Schwarz [49], optimal orthodontic force (OOF) in clinical orthodontics should be equal to capillary blood vessel pressure (≈25 g/cm²) [49]. On a tissue level, OOF should enable the desired clinical outcome without causing the unwanted side effects, for example, root resorption. On the cellular level, it should evoke best biologic cellular response without inhibiting the cell proliferation significantly [27]. Optimal orthodontic force in vitro varies between different models. Estimation of OOF for each in vitro model is of crucial importance for their successful application in OTM simulation [20, 33].

In 2D cell culture WAB in vitro loading models, applied forces varied between 0.2 and 5.0 g/cm². Our data suggest that 2.0 g/cm² was the most commonly used force magnitude in the studies so far. According to Kanzaki et al. [20], this force magnitude proved to induce the best cellular response. Few studies reported a decrease in cell viability in a force-dependent manner, especially with the application of 4 g/cm² force [20, 37, 50, 51].

In studies applying the 3D WAB in vitro loading models, the force magnitude used was chosen depending on the stiffness of the scaffold. Studies using collagen gel scaffolds most commonly applied 6 g/cm² force onto their in vitro models. According to Araujo et al. [47], this force was corresponding to the therapeutic orthodontic force, giving the best cellular response. For PLGA scaffolds, the force magnitude showing the best performance was 25 g/cm² (range: 5–35 g/cm²). The same range of forces were applied in the study of Liao et al. [14] using a hydrophilically modified PLLA scaffold matrix. This range also corresponds to the one used in clinical settings, which indicates that these scaffolds are closest to the mechanical properties of in vivo PDL [14, 33]. This qualifies them also as a suitable model for investigation of light and heavy forces, which are considered as a cause of orthodontic treatment failure.

4.3. Duration of the Force Application

The length of the force application in the studies rarely exceeded 72 h. In most of the cases, force was applied up to 24 and 48 h. Considering the fact that the first 10 days are of crucial importance for OTM ([52], p. 303), the duration of force application in most of the conducted studies is insufficient to fully understand the molecular background of OTM. Additionally, we would like to point out that only a few studies observed cell viability during the experiment. Most of them confirmed a reduction of cell viability, not only due to the force level but also depending on time [19, 50, 51]. We assume that one of the limitations, especially in the 2D WAB in vitro models, is compromised nutrient and oxygen supply in the pressure area. To overcome especially the time limitation of previous models, Liao et al. [14] introduced the hydrophilically modified PLLA matrix as a new scaffold for 3D cultures. They have shown that this scaffold can be used for up to 14 days without affecting cell viability, claiming that it provides good perfusion of the nutrients and oxygen over longer periods of time [14]. Establishing an in vitro model suitable for long-term force application (up to or more than 10 days) is beneficial for progress in this research field.

4.4. Role of PDL and hPDLCs in OTM

Due to lack of PDL, ankylosed teeth and implants cannot undergo OTM, which depict best PDL’s key role in transmitting the mechanical stimulus and initiating the process of bone remodelling [1, 53]. Beside its mechanotransduction properties, it also contributes to tissue homoeostasis and repair, mostly due to the presence of mesenchymal stem cells which are an important part in the normal hPDLC population [4]. This portion of hPDLCs is known to be present in a higher extent in hPDLCs isolated with the “enzyme digestion method” [54], commonly used among the studies in this review, especially in the 3D group.

4.5. Most Examined Genes in the Studies That Used hPDLCs

To explain the contribution of hPDLCs in OTM on the molecular level, we summarised all data regarding the most commonly examined genes and substances in this cell type (Table 1). These were TNFSF11, PTGS2, and PGE₂, known as osteoclastogenesis inducers, and TNFRSF11B, known as an osteoclastogenesis inhibitor.

TNFSF11 (also known as “RANKL”) [55] plays a crucial role in bone resorption on the compression side during OTM, inducing the osteoclast formation. TNFSF11 showed an increased gene expression in all studies that used the 2D WAB in vitro loading model (Table 1). In most of the studies using this model, TNFSF11 gene expression, as well as protein secretion, was positively correlated with both force duration and magnitude reaching the maximum expression level after 12–24 hours of force application. Studies using the 3D WAB in vitro loading model also reported an increase in the TNFSF11 secretion, most of them after 6 hours of force application (Table 1). In cells grown in PLGA scaffolds, a positive correlation between force magnitude and gene expression but a negative correlation between force duration and gene expression was noticed.

TNFRSF11B, also referred to as osteoprotegerin (OPG), is TNFSF11’s antagonist that inhibits osteoclastogenesis [55]. Most of the studies applying the 2D WAB in vitro loading model reported no observed change in gene expression (), with exception of two studies that reported downregulation [40] or transitory downregulation [8] (Table 1). Considering protein secretion, results were contradictory. Most studies, however, reported a decrease in protein secretion or did not report any change. Results from studies using 3D WAB in vitro loading were also contrary, depending on the scaffold used. In a study using collagen gel scaffolds, an increase in TNFRSF11B gene expression was observed [26]. In all studies applying PLGA scaffolds, a decrease in TNFRSF11B secretion was positively correlated with force magnitude and negatively correlated with force duration [27, 28, 31, 33, 43]. With one exception [28], a comparison of TNFSF11 and TNFRSF11B gene expression in the aforementioned studies showed that a rapid down/regulation of TNFRSF11B appears parallel to a rapid upregulation of TNFSF11 in 3D WAB in vitro loading. Since both genes represent antagonists in bone turnover regulation, this was explained as a good representation of the cyclic changes in the bone metabolism on the compression side during OTM [31, 33]. It was also suggested that downregulation of TNFSF11 in later stages might have something to do with other inducers for prolonged osteoclastogenesis promotion [33].

Gene expression of PTGS2 was increased upon force application in both 2D and 3D studies. In most of the 2D WAB studies, PTGS2 showed a positive correlation between the duration of the experiment and gene expression (Table 1). In those studies, using the 3D WAB in vitro loading model, PTGS2 seemed to be negatively correlated with force duration and positively correlated with force magnitude. On the other hand, PTGS2 protein quantity was shown to be in positive correlation with both duration and force magnitude using Western blotting (Table 1). Since PTGS2 is involved in prostaglandin E₂ metabolism, an upregulation of PTGS2 gene expression (maximum at 24 to 48 h after force application) is correlated with an upregulation of PGE₂ secretion (maximum at 48 h after force application) in all studies (Table 1).

Taken together, there seems to be some inconsistency between studies using the 2D and the 3D WAB in vitro loading model. The results within the 2D WAB group of studies are quite similar and comparable. However, a noticeable higher heterogeneity among those studies using the 3D WAB in vitro loading model is recognizable. This heterogeneity can be related to the type of scaffolds used.

4.6. STRING PPI Analysis

We performed STRING PPI analysis for two selected sets of genes (“hPDLC list” and “hOB list”). PPI enrichment values obtained from both PPI networks (Figure 4) had significantly more interactions than expected. This implicates that the genes examined in the studies were not chosen randomly. From our point of view, this is not surprising, since most of the studies were selecting “the genes of interest” for their analysis, all previously known or suspected to be involved in bone metabolism. Just a few of the studies performed microarray analysis in order to identify all genes responding to force application [26, 32, 44, 48].

In addition, KEGG pathways relevant for OTM, identified for each set of genes in STRING analysis (Table 2), can be useful source for discovering new genes that might influence OTM.

5. Conclusions

In summary, the WAB in vitro loading model represents a simple and very efficient way to investigate molecular events during OTM. The purpose of this review was to provide an overview of all used forms of the WAB in vitro loading model (2D and 3D in combination with different scaffolds), present all current findings, and point out at certain questions for their further improvement.

3D WAB in vitro loading models have shown to be promising for use in future research by bringing a more real environment in in vitro setups. However, unlike well-established 2D models that provide comparable results, 3D models show inconsistency in results. Obviously, there is a need for further improvement in order to establish standardised in vitro models that will provide comparable results. Also, there is a need to elucidate molecular events during longer periods of force application. Therefore, the future goal is to establish both 2D and 3D loading models that will allow us to conduct long-term investigations. The study of Liao et al. [14] is a good example for this, and there should be more research in that direction.

Abbreviations

2D: | Two-dimensional
3D: | Three-dimensional
ATP: | Adenosine triphosphate
cAMP: | Cyclic adenosine monophosphate
ECM: | Extracellular matrix
ELISA: | Enzyme-linked immunosorbent assay
H₂S: | Hydrogen sulfide
hOBMCs: | Human oral bone marrow cells
hOBs: | Human osteoblasts
hPDLCs: | Human periodontal ligament cells
KEGG: | Kyoto encyclopedia of genes and genomes
NO: | Nitric oxide
OOF: | Optimal orthodontic force
OPG: | Osteoprotegerin
OTM: | Orthodontic tooth movement
PDL: | Periodontal ligament
PGE₂: | Prostaglandin E₂
PLGA: | Polylactic-co-glycolic acid
PLLA: | Poly-L-lactide acid
PPI: | Protein-protein interaction
PTGS2: | Prostaglandin-endoperoxide synthase 2
RANKL: | Receptor activator of nuclear factor kappa-Β ligand
ROS: | Reactive oxygen species
STRING: | Search tool for the retrieval of interacting genes/proteins
TNF: | Tumor necrosis factor
TNFRSF11B: | TNF receptor superfamily member 11b
TNFSF11: | TNF superfamily member 11
WAB: | Weight approach based.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this manuscript.

Acknowledgments

Mila Janjic received a study grant from BAYHOST (Bayerisches Hochschulzentrum für Mittel-, Ost- und Südosteuropa, Regensburg, Germany) and from the Fund for Young Talents of the Republic of Serbia (Government of the Republic of Serbia, Ministry of Youth and Sports, Belgrade, Serbia).

Supplementary Materials

Supplementary 1. Search strategy designed for the studies applying the in vitro loading model based on a weight approach on cells in 2D or 3D cell culture and lists the excluded studies after full-text reading with reasons.

Supplementary 2. Studies applying the 2D weight approach on human primary cells from the orofacial region, that is, human periodontal ligament cells (hPDLC), human oral bone marrow cells (hOBMC), and human alveolar bone osteoblasts (hOB). For each gene or metabolite force magnitude and force duration, the change in gene expression or substance secretion (increase, decrease, and no change) and the techniques applied are given.

Supplementary 3. Studies applying the 2D weight approach on human and nonhuman cells and cell lines not included in Supplement 2 (i.e., human primary cells from the orofacial region). For each gene or metabolite force magnitude and force duration, the change in gene expression or substance secretion (increase, decrease, and no change) and the techniques applied are given.

Supplementary 4. Studies applying the 3D weight approach on human and nonhuman cells and cell lines. For each gene or metabolite force magnitude and force duration, the change in gene expression or substance secretion (increase, decrease, and no change) and the techniques applied are given.

References

A. Wichelhaus, Orthodontic Therapy - Fundamental Treatment Concepts, Georg Thieme, New York, 2017.

Z. Davidovitch and V. Krishnan, “Biological basis of orthodontic tooth movement: an historical perspective,” in Biological Mechanisms of Tooth Movement, V. Krishnan and Z. Davidovitch Eds, pp. 3–14, Wiley, Chichester, West Sussex, UK, 2015.

Z. Davidovitch, “Tooth movement,” Critical Reviews in Oral Biology and Medicine, vol. 2, no. 4, pp. 411–450, 1991.

B. M. Seo, M. Miura, S. Gronthos et al., “Investigation of multipotent postnatal stem cells from human periodontal ligament,” The Lancet, vol. 364, no. 9429, pp. 149–155, 2004.

U. Baumert, I. Golan, B. Becker et al., “Pressure simulation of orthodontic force in osteoblasts: a pilot study,” Orthodontics and Craniofacial Research, vol. 7, no. 1, pp. 3–9, 2004.

L. Yang, Y. Yang, S. Wang, Y. Li, and Z. Zhao, “In vitro mechanical loading models for periodontal ligament cells: from two-dimensional to three-dimensional models,” Archives of Oral Biology, vol. 60, no. 3, pp. 416–424, 2015.

K. Kanai, H. Nohara, and K. Hanada, “Initial effects of continuously applied compressive stress to human periodontal ligament fibroblasts,” The Journal of Japan Orthodontic Society, vol. 51, pp. 153–163, 1992.

S. J. Kim, K. H. Park, Y. G. Park, S. W. Lee, and Y. G. Kang, “Compressive stress induced the up-regulation of M-CSF, RANKL, TNF-α expression and the down-regulation of OPG expression in PDL cells via the integrin-FAK pathway,” Archives of Oral Biology, vol. 58, no. 6, pp. 707–716, 2013.

P. Proff, C. Reicheneder, A. Faltermeier, D. Kubein-Meesenburg, and P. Römer, “Effects of mechanical and bacterial stressors on cytokine and growth-factor expression in periodontal ligament cells,” Journal of Orofacial Orthopedics / Fortschritte der Kieferorthopädie, vol. 75, no. 3, pp. 191–202, 2014.

D. Moher, PRISMA-P Group, L. Shamseer et al., “Preferred reporting items for systematic review and meta-analysis protocols (PRISMA-P) 2015 statement,” Systematic Reviews, vol. 4, no. 1, 2015.

J. Ye, G. Coulouris, I. Zaretskaya, I. Cutcutache, S. Rozen, and T. L. Madden, “Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction,” BMC Bioinformatics, vol. 13, no. 1, p. 134, 2012.

D. Szklarczyk, J. H. Morris, H. Cook et al., “The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible,” Nucleic Acids Research, vol. 45, no. D1, pp. D362–D368, 2016.

J. Du, Z. Yuan, Z. Ma, J. Song, X. Xie, and Y. Chen, “KEGG-PATH: Kyoto encyclopedia of genes and genomes-based pathway analysis using a path analysis model,” Molecular BioSystems, vol. 10, no. 9, pp. 2441–2447, 2014.

W. Liao, M. Okada, K. Inami, Y. Hashimoto, and N. Matsumoto, “Cell survival and gene expression under compressive stress in a three-dimensional in vitro human periodontal ligament-like tissue model,” Cytotechnology, vol. 68, no. 2, pp. 249–260, 2016.

B. Ragnarsson, G. Carr, and J. C. Daniel, “Basic biological sciences isolation and growth of human periodontal ligament cells in vitro,” Journal of Dental Research, vol. 64, no. 8, pp. 1026–1030, 1985.

M. J. Somerman, S. Y. Archer, G. R. Imm, and R. A. Foster, “A comparative study of human periodontal ligament cells and gingival fibroblasts in vitro,” Journal of Dental Research, vol. 67, no. 1, pp. 66–70, 1988.

S. Benjakul, S. Jitpukdeebodintra, and C. Leethanakul, “Effects of low magnitude high frequency mechanical vibration combined with compressive force on human periodontal ligament cells in vitro,” European Journal of Orthodontics, , In press.

Y. Jin, J. Li, Y. Wang et al., “Functional role of mechanosensitive ion channel Piezo1 in human periodontal ligament cells,” The Angle Orthodontist, vol. 85, no. 1, pp. 87–94, 2015.

Y. G. Kang, J. H. Nam, K. H. Kim, and K. S. Lee, “FAK pathway regulates PGE2 production in compressed periodontal ligament cells,” Journal of Dental Research, vol. 89, no. 12, pp. 1444–1449, 2010.

H. Kanzaki, M. Chiba, Y. Shimizu, and H. Mitani, “Periodontal ligament cells under mechanical stress induce osteoclastogenesis by receptor activator of nuclear factor κB ligand up-regulation via prostaglandin E2 synthesis,” Journal of Bone and Mineral Research, vol. 17, no. 2, pp. 210–220, 2002.

C. Kirschneck, P. Proff, M. Maurer, C. Reicheneder, and P. Römer, “Orthodontic forces add to nicotine-induced loss of periodontal bone : An in vivo and in vitro study,” Journal of Orofacial Orthopedics, vol. 76, no. 3, pp. 195–212, 2015.

L. Liu, K. Igarashi, H. Kanzaki, M. Chiba, H. Shinoda, and H. Mitani, “Clodronate inhibits PGE2 production in compressed periodontal ligament cells,” Journal of Dental Research, vol. 85, no. 8, pp. 757–760, 2006.

K. Mayahara, Y. Kobayashi, K. Takimoto, N. Suzuki, N. Mitsui, and N. Shimizu, “Aging stimulates cyclooxygenase-2 expression and prostaglandin E₂ production in human periodontal ligament cells after the application of compressive force,” Journal of Periodontal Research, vol. 42, no. 1, pp. 8–14, 2007.

S. Premaraj, I. Souza, and T. Premaraj, “Focal adhesion kinase mediates β-catenin signaling in periodontal ligament cells,” Biochemical and Biophysical Research Communications, vol. 439, no. 4, pp. 487–492, 2013.

P. Römer, J. Köstler, V. Koretsi, and P. Proff, “Endotoxins potentiate COX-2 and RANKL expression in compressed PDL cells,” Clinical Oral Investigations, vol. 17, no. 9, pp. 2041–2048, 2013.

R. M. S. de Araujo, Y. Oba, and K. Moriyama, “Identification of genes related to mechanical stress in human periodontal ligament cells using microarray analysis,” Journal of Periodontal Research, vol. 42, no. 1, pp. 15–22, 2007.

M. Li, J. Yi, Y. Yang, W. Zheng, Y. Li, and Z. Zhao, “Investigation of optimal orthodontic force at the cellular level through three-dimensionally cultured periodontal ligament cells,” The European Journal of Orthodontics, vol. 38, no. 4, pp. 366–372, 2016.

J. Yi, B. Yan, M. Li et al., “Caffeine may enhance orthodontic tooth movement through increasing osteoclastogenesis induced by periodontal ligament cells under compression,” Archives of Oral Biology, vol. 64, pp. 51–60, 2016.

K. Mayahara, A. Yamaguchi, M. Sakaguchi, Y. Igarashi, and N. Shimizu, “Effect of Ga-Al-As laser irradiation on COX-2 and cPLA₂-α expression in compressed human periodontal ligament cells,” Lasers in Surgery and Medicine, vol. 42, no. 6, pp. 489–493, 2010.

S. Wongkhantee, T. Yongchaitrakul, and P. Pavasant, “Mechanical stress induces osteopontin expression in human periodontal ligament cells through rho kinase,” Journal of Periodontology, vol. 78, no. 6, pp. 1113–1119, 2007.

M. Le Li, J. Yi, Y. Yang et al., “Compression and hypoxia play independent roles while having combinative effects in the osteoclastogenesis induced by periodontal ligament cells,” The Angle Orthodontist, vol. 86, no. 1, pp. 66–73, 2016.

Y. Li, M. Li, L. Tan et al., “Analysis of time-course gene expression profiles of a periodontal ligament tissue model under compression,” Archives of Oral Biology, vol. 58, no. 5, pp. 511–522, 2013.

Y. Li, W. Zheng, J. S. Liu et al., “Expression of osteoclastogenesis inducers in a tissue model of periodontal ligament under compression,” Journal of Dental Research, vol. 90, no. 1, pp. 115–120, 2010.

S. Y. Lee, H. I. Yoo, and S. H. Kim, “CCR5-CCL axis in PDL during orthodontic biophysical force application,” Journal of Dental Research, vol. 94, no. 12, pp. 1715–1723, 2015.

F. Liu, F. Wen, D. He et al., “Force-induced H₂S by PDLSCs modifies osteoclastic activity during tooth movement,” Journal of Dental Research, vol. 96, no. 6, pp. 694–702, 2017.

P. Luckprom, K. Kanjanamekanant, and P. Pavasant, “Role of connexin43 hemichannels in mechanical stress-induced ATP release in human periodontal ligament cells,” Journal of Periodontal Research, vol. 46, no. 5, pp. 607–615, 2011.

M. Mitsuhashi, M. Yamaguchi, T. Kojima, R. Nakajima, and K. Kasai, “Effects of HSP70 on the compression force-induced TNF-α and RANKL expression in human periodontal ligament cells,” Inflammation Research, vol. 60, no. 2, pp. 187–194, 2011.

R. Nakajima, M. Yamaguchi, T. Kojima, M. Takano, and K. Kasai, “Effects of compression force on fibroblast growth factor-2 and receptor activator of nuclear factor kappa B ligand production by periodontal ligament cells in vitro,” Journal of Periodontal Research, vol. 43, no. 2, pp. 168–173, 2008.

Y. Nishijima, M. Yamaguchi, T. Kojima, N. Aihara, R. Nakajima, and K. Kasai, “Levels of RANKL and OPG in gingival crevicular fluid during orthodontic tooth movement and effect of compression force on releases from periodontal ligament cells in vitro,” Orthodontics and Craniofacial Research, vol. 9, no. 2, pp. 63–70, 2006.

K. Yamada, M. Yamaguchi, M. Asano, S. Fujita, R. Kobayashi, and K. Kasai, “Th17-cells in atopic dermatitis stimulate orthodontic root resorption,” Oral Diseases, vol. 19, no. 7, pp. 683–693, 2013.

M. Yamaguchi, N. Aihara, T. Kojima, and K. Kasai, “RANKL increase in compressed periodontal ligament cells from root resorption,” Journal of Dental Research, vol. 85, no. 8, pp. 751–756, 2006.

M. Kaku, J. M. Rosales Rocabado, M. Kitami et al., “Mechanical loading stimulates expression of collagen cross-linking associated enzymes in periodontal ligament,” Journal of Cellular Physiology, vol. 231, no. 4, pp. 926–933, 2016.

Y. I. Jianru, L. I. MeiLe, Y. YANG, W. ZHENG, L. I. Yu, and Z. ZHAO, “Static compression regulates OPG expression in periodontal ligament cells via the CAMK II pathway,” Journal of Applied Oral Science, vol. 23, no. 6, pp. 549–554, 2015.

K. L. Kang, S. W. Lee, Y. S. Ahn, S. H. Kim, and Y. G. Kang, “Bioinformatic analysis of responsive genes in two-dimension and three-dimension cultured human periodontal ligament cells subjected to compressive stress,” Journal of Periodontal Research, vol. 48, no. 1, pp. 87–97, 2013.

J. Kikuta, M. Yamaguchi, M. Shimizu, T. Yoshino, and K. Kasai, “Notch signaling induces root resorption via RANKL and IL-6 from hPDL cells,” Journal of Dental Research, vol. 94, no. 1, pp. 140–147, 2014.

M. Weinreb and C. E. Nemcovsky, “In vitro models for evaluation of periodontal wound healing/regeneration,” Periodontology, vol. 68, no. 1, pp. 41–54, 2015.

R. M. Santos de Araujo, Y. Oba, and K. Moriyama, “Role of regulator of G-protein signaling 2 (RGS2) in periodontal ligament cells under mechanical stress,” Cell Biochemistry and Function, vol. 25, no. 6, pp. 753–758, 2007.

Y. H. Lee, D. S. Nahm, Y. K. Jung et al., “Differential gene expression of periodontal ligament cells after loading of static compressive force,” Journal of Periodontology, vol. 78, no. 3, pp. 446–452, 2007.

A. M. Schwarz, “Tissue changes incidental to orthodontic tooth movement,” International Journal of Orthodontia, Oral Surgery and Radiography, vol. 18, no. 4, pp. 331–352, 1932.

K. Kanjanamekanant, P. Luckprom, and P. Pavasant, “Mechanical stress-induced interleukin-1beta expression through adenosine triphosphate/P2X7 receptor activation in human periodontal ligament cells,” Journal of Periodontal Research, vol. 48, no. 2, pp. 169–176, 2013.

P. Tripuwabhrut, K. Mustafa, P. Brudvik, and M. Mustafa, “Initial responses of osteoblasts derived from human alveolar bone to various compressive forces,” European Journal of Oral Sciences, vol. 120, no. 4, pp. 311–318, 2012.

W. R. Proffit and H. W. Fields Jr., Contemporary Orthodontics, Mosby, St. Louis, 2000.

G. J. King, S. D. Keeling, and T. J. Wronski, “Histomorphometric study of alveolar bone turnover in orthodontic tooth movement,” Bone, vol. 12, no. 6, pp. 401–409, 1991.

K. Tanaka, K. Iwasaki, K. E. Feghali, M. Komaki, I. Ishikawa, and Y. Izumi, “Comparison of characteristics of periodontal ligament cells obtained from outgrowth and enzyme-digested culture methods,” Archives of Oral Biology, vol. 56, no. 4, pp. 380–388, 2011.

M. Yamaguchi, “RANK/RANKL/OPG during orthodontic tooth movement,” Orthodontics & Craniofacial Research, vol. 12, no. 2, pp. 113–119, 2009.