Abstract

Tumor necrosis factor-α is a common cytokine that increases in inflammatory processes, slows the differentiation of bone formation, and induces osteodystrophy in the long-term inflammatory microenvironment. Our previous study confirmed that the Elongation protein 2 (ELP2) plays a significant role in osteogenesis and osteogenic differentiation, which is considered a drug discovery target in diseases related to bone formation and differentiation. In this study, we applied an in silico virtual screening method to select molecules that bind to the ELP2 protein from a chemical drug molecule library and obtained 95 candidates. Then, we included 11 candidates by observing the docking patterns and the noncovalent bonds. The binding affinity of the ELP2 protein with the candidate compounds was examined by SPR analysis, and 5 out of 11 compounds performed good binding affinity to the mouse ELP2 protein. After in vitro cell differentiation assay, candidates 2# and 5# were shown to reduce differentiation inhibition after tumor necrosis factor-α stimulation, allowing further optimization and development for potential clinical treatment of inflammation-mediated orthopedic diseases.

1. Introduction

Persistent bone tissue inflammation, such as bone fractures and rheumatoid arthritis (RA), destroys the bone formation and absorption balance, which reduces bone mass and results in significant impact on bone regeneration [1]. So far, no effective drugs can cope with this clinical problem. We have noticed that elevated tumor necrosis factor-a (TNF-α) is a major contributor to the bone pathophysiology through inhibition of osteoblasts function and stimulation of bone resorption activity of osteoclasts [2, 3]. Drug for lessening the potential damage caused by long-term TNF-α-associated inflammation exposure is urgent for clinical orthopedics.

Elongation protein 2 (ELP2) is a subordinate of the ELP123 complex and regulates the activity of elongation protein by integrating signals from various factors. It affects the inflammatory response through the JAK-STAT3 cascade and downstream pathways [4]. Our previous research identified the role of ELP2 in regulating TNF-α-induced osteoblast differentiation inhibition. A mechanistic study concluded that ELP2 blocked the osteogenic differentiation induced by BMP-2 through activation the STAT3 pathway and downregulating the utterance of BMPR2, thereby slowing the process of early bone tissue formation in the inflammatory microenvironment [5]. Thus, ELP2 was shown to be the potential drug target to develop novel therapeutic for inflammation-induced bone loss.

In silico screening methods, such as quantitative structure-activity relationships (QSAR), pharmacophore models, molecular dynamics (MD) simulations, and high-throughput molecular docking, have been used to find specific protein inhibitors in recent research [69]. In our study, we used the molecular docking method to screening novel structure-specific chemicals for ELP2 and obtained eleven chemical compounds, which showed effective binding affinity to ELP2 active pockets. Combined with label-free surface plasmon resonance (SPR) affinity analysis, five candidates showed obvious binding affinity to ELP2. After in vitro experiments, two of the candidate compounds (2 # and 5 #) were shown to block the retardation of osteoblast differentiation induced by TNF-α, leading a potent clinical treatment potential for inflammation-induced bone loss.

2. Materials and Method

2.1. Protein and Ligand Library Preparation

The target protein was obtained from the UniProtKB database (https://www.uniprot.org/) with an entry ID of Q91WG4 (Mus musculus), and the homocrystal structure was selected with the advanced sequence searching function aiming for high similarity. Then, the raw composition of the ELP2 was got from the database (https://www.rcsb.org/) with the identification number 5M2N. The homobuilding step was performed with Structure Prediction Wizard in the Schrodinger software (Maestro, version 2015, Schrödinger LLC, USA). The missing loops and mismatches were corrected using a knowledge-based algorithm provided with the wizard module. To ensure the correct starting structure, the initial structure of the protein 3D model was prepared using the Protein Preparation Wizard module. Hydrogen atoms were increase in structure consistently with physiological pH (7.0). Then, the structure was optimized by adjusting hydrogen bonding and removing atom collisions, and formal charges are added to the heteron and optimized at neutral pH. Finally, the redundant water molecules were removed, and the new structure was minimized using the OPLS-2005 force field optimization potential. The prepared protein structure was subjected to binding site calculation.

The binding sites were generated using the SiteMap tool of the Schrödinger comprehensive assessed according to the calculated attributes, including the size of the site, the degree of protein enclosure, the degree of solvent exposure, the degree of tightness of site interaction with the protein, and the hydrophilic and hydrophobic character of the site. The total site scores were then used as the main basis to select potential binding sites. Finally, we completed the receptor grid generation step ranked the sitemaps from high to low score.

The chemical ligand library, including 552,007 compounds, was retrieved from the InterBioScreen database (InterBioScreen ltd., Russia) in structure data file (.SDF) format (Mar 2018 version). The molecules were arranged to ligand preparation using the LigPrep module with default settings. The ligand preparation process involved saving the definite chiralities and generating five minimum low-energy stereoisomers per ligand with pH values ranging from 5.0 to 9.0.

2.2. Virtual Screening

During the docking stage, the ligands were docked to the grids by three-tire-docking [10], which started with HTVS, followed by SP and XP to improve the accuracy. Then, we screened and chose the ligand according to the ranking of the GLIDE score. Herein, the optimized formula of the GLIDE score was as follows: [11]. The ADME properties of the selected ligands were analyzed using the QikProp tool to conclude the natural properties of the drugs, such as hydrogen bond donor and acceptor even the number of rotatable bonds.

2.3. Surface Plasmon Resonance (SPR) Assay

To validate the binding affinity of the candidate compounds with ELP2 protein, we tested the binding affinity through SPR assay. The experiment was performed as previous reported [12]. Briefly, the candidate molecules are crossed onto the chip surface by hydrogen substitution reaction. During the experimental SPR, the chips were first primed with running buffer (1x PBS containing 5% DMSO). Mouse ELP2 protein used in the assay was diluted with running buffer at the concentrations from 200 nM to 3200 nM. To validate the reliability of our method, we selected biotin and rapamycin binding with FKBP12 as a system control. Additionally, the solvents DMSO for compounds and PBS () for proteins were measured individually as blank controls and background noise controls.

2.4. Cell Culture and Treatment

C2C12 cells were kept in DMEM complemented with the 10% fetal dairy serum and 1% penicillin/streptomycin and MC3T3-E1 cells in minimally needed culture medium (MEM, HyClone, GE Healthcare, USA) supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin, were grown individually, and then cultured at 37.0°C, 5% CO2. As previously reported [12], cells were cultured for differentiation induction in the presence of TNF-α with or without our candidate compounds.

2.5. Alkaline Phosphatase Measurement

Bone marrow transplantation was detected by ALP activity via the alkaline phosphatase assay kit followed by the instruction. A BioTek microplate spectrophotometer was used to assay the apparent density at the wavelength of 405 nm. ALP activity values are equally common in the density protein found in stem cells determined to use the BCA Protein Assay Kit.

2.6. Mineralization Assessment

As reported, in two differentiation cell models, the culture medium was replaced every 3 days till 35th day. Mineralization of the cell was found by staining using the Optimized Alizarin Red S Stain Kit by the kit instructions.

2.7. Real-Time Fluorescent Quantitative PCR

The total mRNA was extracted from cell culture using TRIzol reagent (Thermo Scientific USA). The concentration and quality of the RNA were checked using a NanoDrop™ 8000 spectrophotometer, and cDNA was then synthesized. Gene expression according to the recommendations of the reverse transcription kit was detected. Host gene β-actin was selected as the reference gene, and the gene expression level analysis was calculated by the 2-ΔΔCt quantification method [13].

2.8. Western Blot Assay

After treatment, using PBS to wash the cells, the cells were lysed with RIPA lysis buffer. The cell lysate was then separated by SDS-PAGE. All proteins were transferred to PVDF membranes (Bio-Rad, USA). The PVDF membrane was then blocked with 5% skim milk in Tris-buffered physiological saline -0.1% Tween-20 for 1 hour and incubated with the following primary antibody: anti-ALP antibody (rabbit), anti-COL-1 antibody (mouse), anti-OCN antibody (mouse), anti-BMP-2 antibody (mouse), anti-Runx2 antibody (mouse), Anti-Osterix (rabbit), and anti-β-actin (mouse), followed by suitable horseradish peroxidase- (HRP-) labeled secondary antibody, blocking buffer for goat anti-rabbit IgG. Images were captured using G&E imaging system.

2.9. Statistical Analysis

Multigroup comparison has been made with the two ANOVA methods. A single comparison was made using a one-sided Student’s test ( was counted as significance). Analysis was completed using SPSS 13.0 Statistics Package and Background 2017. The data is displayed as . The test was tripled and at least once.

3. Results

3.1. Modelling Preparation and Virtual Screening

The mouse ELP2 protein structure was homobuilt using PDB entry 5M2N. Missing loops were filled and mismatches were corrected through a knowledge-based algorithm. A preliminary model was then set up using the Protein Preparation Wizard. After further structure optimization, binding sites were generated and receptor grids were calculated by a grid generation module. The chemical ligand library consisting of 552,007 compounds was prepared with the LigPrep module, yielding five minimum low-energy stereoisomers per ligand with pH values ranging from 5.0 to 9.0. In the ligand-docking stage, the ELP2 model and the entire ligand library were first applied to perform a high-throughput virtual screening (HT-VS), and 11,216 compounds were passed. These screened compounds were further analyzed by the standard precision (SP) and extraprecision (XP) glide docking models. 397 hits belonging to 102 clusters were retrieved from the Glide SP docking with a docking score cut-off value of -6.0 kcal/mol. Subsequently, the Glide XP module was used to further reduce the hit set, and 95 hits belonging to 36 clusters were retrieved with a docking score cut-off value of -5.0 kcal/mol. Finally, these 95 ligand molecules were visually inspected based on docking poses and bonds interacting with protein residues. Finally, 11 hits were selected as candidate compounds. The structures of the 11 compound candidates are shown in Figure 1. The binding sites and docking surface structures are shown in Supplementary Figure 1A and 1B.

3.2. SPR Affinity Analysis

Under the technique of photo-cross-linker biosensor chip, the candidate compounds were immobilized onto the chemically modified chip surface in random orientations and without any attached label or linker. After the compounds were fixed onto the surface of the chip, ELP2 protein solved in PBS were passed over the chip surface. The original sensor gram information was collected in real time. The association rate constants (), dissociation rate constants (), and the equilibrium dissociation constant (KD) of samples are shown in Table 1. The binding curves are shown in Figure 2. According to the affinity measurement results, five among the eleven candidates possessed strong binding affinity to ELP2, and another two candidates (1# and 9#) showed much faster dissociation rates and were excluded due to presumable difficulty in maintaining adequate efficacy in vivo.

3.3. 2# and 5# Compounds Show Remarkable Mitigation against TNF-α-Induced ALP Activity Inhibition

To determine the effect of the candidate compounds on TNF-α-induced osteogenic differentiation inhibition, we examined the activity of osteogenic marker alkaline phosphatase in the cells. As mentioned in Materials and Method, cells were grown in an adequate osteoblast differentiation environment and tested the alkaline phosphatase activity. Staining observations and relative ALP activity are showed in Figures 3. The results indicated that compared with the osteogenic differentiation controls of the ID and SD differentiation models, TNF-α significantly reduced ALP activity and lightened the color after ALP staining, while the presence of candidates 2# and 5# brought a significant ease to TNF-α-induced inhibition of osteoblast differentiation in C2C12 (Figures 3(a) and 3(b)) and MC3T3-E1 cell (Figures 3(c) and 3(d)). The effects of candidates 7#, 8#, and 11# were not as strong as 2# and 5#, with only slightly easing of the inhibition resulting from TNF-α (data not shown). The half-effective concentration (EC50) of candidates 2# and 5# were measured through examining the ALP activity (Figure 4); the EC50 of 2# candidate was 14.71 μM in C2C12 cells of the ID model; at the same time, the EC50 in the MC3T3-E1 SD model was 29.68 μM. The EC50 values of 5# candidate were 12.15 μM and 18.24 μM in ID and SD model cells, respectively. Candidate 5# possessed better protective effect compared with candidate 2#, which was consistent with SPR binding affinity test result that 5# showed higher binding affinity to ELP2 protein (Figure 2).

3.4. 2# and 5# Compounds Reverse TNF-α-Induced Inhibition on Osteoblast Mineralization Activity

In order to test the function of the candidate compounds on the mineralization activity during osteoblasts differentiation, the self-differentiation (SD) model and inflammatory-differentiated (ID) model cells were stained with Alizarin Red S staining on day 35 after differentiation induction, and the mineralization activity of each group can be estimated by the number of stained spots. As Figure 5 has shown, extracellular mineralization was stained with Alizarin Red S staining; compared with the TNF-α control groups, candidates 2# and 5# increased the red spot on the cell surface, suggesting that cell mineralization can be reversed by the presence of candidates 2# and 5# in C2C12 and MC3T3-E1 cells (Figure 5), and we inferred that the TNF-α-induced osteoblast differentiation inhibition was reversed based on the appearance of mineralization activity changes in the cell models.

3.5. 2# and 5# Compounds Increased the Expression Level of Osteoblast Differentiation-Related Marker after TNF-α Stimulation

In order to investigate the changes in osteogenic differentiation associated molecules including COL-I, ALP, OCN, BMP-2, RUNX-2, and OSX genes, we performed Q-PCR and Western blot analysis on the differentiation cellular models on the 7th day after osteoblast cell differentiation are shown in Table 2. According to the difference of the gene expression levels, the differentiation control group was significantly upregulated compared with TNF-α treated cells on the day 7. The level of above gene expression was raised in the 2# and the 5# candidate-treated cell with cell morphology phenotype differentiation when compared to the cell stimulated with TNF-α (Figure 6). A further study of the expression of osteoblastic markers at the protein level was performed using Western blot. As shown in Figure 7, the expressions of the marker proteins were consistent with the mRNA expressions. Therefore, our above experiments indicated that the 2# and 5# ELP2 inhibitor candidates could effectively lighten the TNF-α-induced inhibition of primary osteoblast differentiation during bone inflammation.

4. Discussion

With different bone cells working together, it completes the role of absorbing the old bone and generating a new bone. On this way, the morphology and structure of the bone will be constantly destroyed and remodeled, and this depends on the number and proportion of cells, so the differentiation of osteoblasts can affect the morphology and structure of bone. Persistent inflammatory bone environment hinders osteoblastic differentiation, leading to bone dysplasia [1416]. TNF-α is a main contributor to persistent inflammation- and infection-induced bone regeneration inhibitions; the underline mechanism is still not clear [17]. A previous study [6] reported that ELP2 positively regulated the increased TNF-α-induced repression of osteoblast differentiation during persistent inflammation. In clinical studies, the repression of osteoblast differentiation will affect the regeneration and formation of bone. Therefore, ELP2 can inhibit bone dysplasia induced by TNF-α, which leads it playing a role in long-term chronic infections.

Molecular docking methods are an effective method to identify potential chemical structures for novel targets [18]. Compared with traditional high-throughput screening, the computational calculation method is faster and more cost-effective. The Schrödinger suit is also an efficient ligand screening program, which is widely used and has been proven in the field of ligand-seeking research [19, 20]. In this study, we performed virtual docking following the published method and using the IBS compound library and obtained 11 candidates that fit into the active site of the homobuilt mouse ELP2 protein. These candidate compounds were prepared for SPR affinity analysis with ELP2 protein, which is a technique that enables rapid detection of the interaction between two molecules together with their binding kinetic parameters in real time [21]. Up to now, no research has reported the epitope of ELP2 active center; therefore, we simulated some possible active centers using the simulation method with the protein structure building module. It is undeniable that such a method results in a large number of incorrect active centers and ineffective ligand candidates in the screening. We speculated that this was the reason only 2 of the 11 high-scored candidate compounds showed inhibition activity in the in vitro experiments.

During our experiment, we applied an in silico virtual screening to select compound molecules from the chemical drug library that bind to the ELP2 protein and obtained 95 candidates. We further filtered out 11 candidates by scoring the docking model and observing noncovalent bonds. Afterwards, the binding affinities of the ELP2 protein with candidate compounds were analyzed by SPR assay, and 5 of the 11 compounds possessed obvious binding affinity to ELP2 protein. Accordingly, these 5 potential compounds were used for the in vivo efficiency study to determine whether there was an associated potential drug effect. The candidate was tested in the differentiation model and verified, which are used for evaluation of osteoblast differentiation [22]. Further, Q-PCR and Western blot analyses were used to investigate the osteogenic differentiation markers which was affected by the two compounds; the results were consistent with previous, which indicated that 2 compounds could mitigate the effect of TNF-α-induced differentiation inhibition.

In summary, we obtained two competitive inhibitors of ELP2 protein; both of them have the potential to impede TNF-α-induced new bone-forming osteoblasts (OB) inhibition in two cells. The EC50 of candidate 2 # was 14.71 μM in the C2C12 cells of the ID model, while the EC50 in the MC3T3-E1 SD model was 29.68 μM. Candidate 5 # EC50 values were 12.15 μM and 18.24 μM in the ID and SD model cells, respectively. The results about the bone differentiation assay were in line with expectations. The results give hope that the 2 # and 5 # candidates will become new drugs in clinical orthopedics. Finally, further research will spare no effort to clarify the additional pharmacodynamics and metabolic kinetics of the 2# and 5# candidates, so as to develop them into clinical drugs that prevent the inhibition of bone regeneration caused by persistent inflammation.

Abbreviations

COL-I:Collagen alpha-1
EC50:Effective concentration 50
EC:Elongator complex
ELP123:Elongator subunits 1-3
ELP456:Elongator subunits 4-6
ELP2:Elongator protein 2
FBS:Fetal bovine serum
GLIDE:Grid-based ligand docking with energetics
HRP:Horseradish peroxidase
HTVS:High-throughput virtual screening
ID model:Inflammatory differentiation model
KD:Dissociation constant,
MD:Molecular dynamics
OCN:Osteocalcin
OSX:Osterix
Pol II:RNA polymerase II
QSAR:Quantitative structure-activity relationships
RA:Rheumatoid arthritis
RU:Response units
RUNX-2:Runt-related transcription factor 2
SD model:Self-differentiation model
SLE:Systemic lupus erythaematosus
SP:Standard precision
SPR:Surface plasmon resonance
StIP1:STAT3-interacting protein 1
XP:Extraprecision mode.

Data Availability

Once upon request, the corresponding author will offer the data included in this study.

Conflicts of Interest

The authors declared no conflicts of interests to this study.

Authors’ Contributions

Wen-Jiao Wu and Chang-Liang Xia contributed equally to this work. Wen-Jiao Wu and Chang-Liang Xia did the experiment and contributed equally to this research and should be regarded as co-first author. Chang-Peng Xu and Yong Qi contribute to the experimental plan design and funding application. The main work of Shuan-Ji Ou, Yang Yang, Yun-Fei Ma, and Yi-Long Hou is in silico docking and data analysis, including chart making and figure design. The work of Qing-Po Yang, Jun Zhang, and Jian-Wei Li was reporting of the work, including drafting, copywriting, and revising.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (No. 81972083), Science and Technology Planning Project of Guangzhou (Nos. 202102080052, 202102010057, and 201804010226), Science Foundation of Guangdong Second Provincial General Hospital (Nos. 3D-A2020004, 3D-A2020002, YQ2019-009, and 2019BSGZ012), Natural Science Foundation of Xinjiang Province (No. 2018D01C014), and Minority Science and Technology Talent Special Training Program Project of Xinjiang Province (No. 2019D03025).

Supplementary Materials

Supplementary Figure 1: the binding sites and docking surface structures of eleven candidates in docking procedure. (Supplementary Materials)