Skip to main content

A study on the anti-osteoporosis mechanism of isopsoralen based on network pharmacology and molecular experiments



Osteoporosis (OP) is a disease caused by multiple factors. Studies have pointed out that isopsoralen (IPRN) is one of the most effective drugs for the treatment of OP. Based on network pharmacological and molecular experimental analysis, the molecular mechanism of IPRN in osteoporosis is clarified.


IPRN target genes and OP-related genes were predicted from the databases. Intersections were obtained and visualized. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed on target genes, which was confirmed by experiments internal and external experiments. Molecular docking was used to verify the binding between IPRN and target proteins. Molecular dynamics (MD) simulates the binding affinity of protein targets and active compounds.


87 IPRN target genes and 242 disease-related targets were predicted. The protein–protein interaction (PPI) network identified 18 IPRN target proteins for the treatment of OP. GO analysis indicated that target genes were involved in biological processes. KEGG analysis showed that pathways such as PI3K/AKT/mTOR were associated with OP. Cell experiments (qPCR and WB) found that the expressions of PI3K, AKT, and mTOR in MC3T3-E1 cells at 10 μM, 20 μM, and 50 μM IPRN concentrations, especially at 20 μM IPRN treatment, were higher than those in the control group at 48 h. Animal experiments also showed that compared with the control group, 40 mg/kg/time IPRN could promote the expression of the PI3K gene in chondrocytes of SD rats.


This study predicted the target genes of IPRN in the treatment of OP and preliminarily verified that IPRN plays an anti-OP role through the PI3K/AKT/mTOR pathway, which provides a new drug for the treatment of OP.


Osteoporosis (OP) is a chronic disease with impaired bone strength caused by low bone density and bone microstructure damage, which mainly occurs in adults, especially postmenopausal women [1, 2]. Since bone density and bone quality are specific manifestations of bone strength, OP will damage bone mineral density, leading to increased fragility and further fracture [3]. The causes of OP are mainly related to low estrogen content, insufficient calcium intake and vitamin D supplementation, and excessive use of glucocorticoids [4, 5]. The literature indicates that corticosteroid-induced osteoporosis (CIO) is the most common type of secondary osteoporosis, which can lead to fracture and increase morbidity and mortality [6]. Statistical research found that through meta-analysis, the prevalence of OP in countries around the world was 18.3% (95% CI: 16.2–20.7), and the prevalence of OP in Africa was the highest at 39.5% (95% CI: 22.3–59.7) [7]. In addition, it pointed out that in a cross-sectional study of 20,416 people, the prevalence of OP among Chinese adults aged 40 and above in the mainland was 5.0% in males and 20.6% in females [8]. Obviously, with such a serious aging population in China, the incidence of OP is increasing year by year. Therefore, researching the mechanism and treatment of OP has become a popular direction. At present, its specific pathogenesis and treatment are still under in-depth research.

Currently, On the one hand, the popular biochemical markers of bone turnover (BTM), such as serum bone alkaline phosphatase (bALP) and type I procollagen N propeptide (PINP), provide specific and dynamic indicators of bone turnover mechanism in the balance between bone formation and bone absorption, and can be used as predictive factors of OP occurrence and monitoring indicators of OP patients [9]. On the other hand, hormone replacement therapy (HRT) is used to treat OP in clinical practice. Raloxifene, for example, is a selective estrogen receptor modulator (SERM) for the treatment and prevention of OP. It has been shown to reduce the risk of vertebral fractures by 36% in large clinical trials. Therefore, it greatly improves estrogen's protection of bones from fractures [10]. In addition, it also includes other anti-absorption prescription drugs for OP, such as bisphosphonate (BP) (alendronate, risedronate, and so on) and RANK ligand inhibitors (such as denosumab) [11, 12]. However, clinical findings have suggested that although commonly used anti-OP drugs are effective, severe side effects and loss of potency may limit the long-term use of single drugs [13]. And the long-term use of hormone therapy in patients will increase the risk of cardiovascular disease, coronary heart disease, breast cancer, and other diseases to varying degrees [14].

Therefore, we try to find a suitable and effective traditional Chinese medicine (TCM), to reduce the harm of side effects to the body as much as possible. Previous studies have shown that Epimedium [15], Angelica Sinensis root [16], and Du-Zhong (Eucommia ulmoides) cortex extracts (DZCE) [17] in TCM can prevent or delay the onset of OP. In addition, another TCM, Psoralea corylifolia Linn (PCL) also plays an important role in the treatment of OP. The study has found that PCL is a kind of coumarin-type TCM compound to prevent OP, and its main active components are psoralen (PRN) and isopsoralen (IPRN) (Fig. 1) [18]. IPRN, a component of psoralen seeds isolated from previous studies, is the same as the angelica in the root of Angelica sinensis. PCL has been widely used in various TCM formulations for the treatment of cardiovascular disease, osteoporosis, cancer, and other diseases [18]. Previous literature has concluded that PRN and IPRN have biological activities such as anti-tumor [19, 20], increased levels of estrogen [21], osteoblast differentiation [22], and so on. At present, there are many studies on PRN, but few studies on the mechanism of IPRN in preventing and treating OP. Therefore, IPRN is also a good choice of TCM when studying OP caused by hormone deficiency.

Fig. 1
figure 1

Pharmacological and molecular properties of isopspralen by SwissADME

At this stage, the development of network pharmacology is particularly rapid and its application is extensive. This is a new discipline based on the theory of systems biology, which analyzes the network of biological systems and selects specific signal nodes to construct a molecular network of multi-target drugs and multi-way regulation [23]. Improving the therapeutic effect of drugs and reducing toxic and side effects, thereby increasing the success rate of new drugs in clinical trials and saving the cost of drug research and development are emphasized in network pharmacology. It can be seen that the key to the modernization of TCM is to find suitable models and methods for the study of the TCM system and effective therapeutic drugs.

Therefore, based on network pharmacology, this study explored the role of IPRN in OP through different molecular targets and different pathways and was verified by cell experiments and animal experiments, aiming to provide a new direction for the treatment of OP.

Materials and methods

Query and visualization IPRN target genes

"Isopsoralen" was entered and searched in the PubChem ( Compound database, and the retrieved target genes were exported and saved in Excel. Then the STRING ( database and DAVID ( database were used to unify gene names and IDs [24, 25]. The obtained target genes were imported into Cytoscape 3.7.0 software to visualize the network diagram of drugs and target genes [26].

Acquisition of OP target genes

"Osteoporosis" was sequentially entered into PubMed (, GeneCards (, and OMIM (Online Mendelian Inheritance in Man) ( databases for retrieval, and the obtained genes were integrated. The duplicated genes were removed, and thus, the target-related genes of OP were finally obtained. According to the method of 2.1, Cytoscape 3.7.0 software was still used to visualize the network diagram of the disease.

Screening of IPRN-OP intersected genes and PPI network construction and visualization

Drug target genes and OP target genes were input by Excel to make a Venn diagram. The crossover genes obtained in this section served as potential targets of IPRN to treat OP in later studies.

The common targets of IPRN-OP were input into the STRING database for retrieval. The protein type was set as "Homo sapiens", and the minimum interaction threshold was set as 0.4. Then, the network relationship data of the target interaction was obtained, and the PPI network was drawn [24]. Finally, Cytoscape software was used for visualization.

GO annotation analysis and KEGG enrichment analysis

The Metascape platform was used to perform Gene Ontology (GO) annotation analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the targets. And the biological processes, cellular components, molecular functions, and the key signaling pathways of the targets were obtained. The analytical results were then visualized by GraphPad Prism 8.0, with P < 0.05 and false discovery rate (FDR) < 0.05 as the screening criteria.

Molecular docking verification of key targets and drug components

The predicted target is evaluated by molecular docking. First, the crystal structures of candidate components and targets are obtained from the TCMSP database and RCSB protein database respectively ( Download. Then, AutoDockTools 1.5.6 software is used to read candidate molecules and target protein files, hydrogenate, calculate charges, add atomic types (candidate molecules as ligands, and proteins as receptors), and saves them in pdbqt format. Then, import AutoDockTools 1.5.6 to build docking grid boxes for each target. The docking is completed through Autodock Vina 1.1.2 software. The lower the vina score is, the more stable the binding between ligand and receptor is, which is used to preliminarily evaluate the binding activity of compounds to targets. Finally, Pymol software was used to output the protein and docking results as complexes, and discovery Studio was used to analyze the force.

Molecular dynamics simulation

To analyze the binding affinities of protein targets and active compounds obtained by molecular docking, molecular dynamics (MD) simulations were carried out by GROMACS 2020.3 software. The MD simulation workflow included four steps: energy minimization, heating, equilibrium, and production dynamics simulation. The simulation box size was optimized with the distance between each atom of the protein and the box greater than 1.0 nm. To make the simulation system electrically neutral, the water molecules were replaced with Cl and Na+ ions. Following the steepest descent method, energy optimization of 5.0 × 104 steps was performed to minimize the energy consumption of the entire system, and finally to reduce the unreasonable contact or atom overlap in the entire system. After energy minimization, first-phase equilibration was performed with the NVT ensemble at 300 K for 100 ps to stabilize the temperature of the system. Second-phase equilibration was simulated with the NPT ensemble at 1 bar and 100 ps. The primary objective of the simulation is to optimize the interaction between the target protein and the solvent and ions so that the simulation system is fully pre-equilibrated. All MD simulations were performed for 100 ns under an isothermal and isostatic ensemble with a temperature of 300 K and a pressure of 1 atmosphere. Lennard–Jones function was used to calculate the Van der Waals force, and the nonbond truncation distance was set to 1.4 nm. The long-range electrostatic interaction was calculated by the particle Mesh-Ewald method with the Fourier spacing of 0.16 nm.

Molecular biological verification

Drug-regulated cell experiments

IPRN was purchased from MCE Co., Ltd. (MedChemExpress, Monmouth Junction, NJ, USA). The product number of IPRN was HY-N0763, and its purity is 99.86%. The mouse osteoblast cell line (MC3T3-E1) was obtained from Procell Life Science&Technology Co., Ltd. (Wuhan, China).

MC3T3-E1 cells were cultured in DMEM medium containing 10% fetal bovine serum (FBS) and then were passaged and cryopreserved. The prepared MC3T3-E1 cell suspensions were divided into the control group and the experimental group. In the experimental group, three gradients of IPRN, 10 μM (low concentration), 20 μM (medium concentration), and 50 μM (high concentration), were sequentially added to the complete medium according to the conditions of the previous experiments. At 48 h after administration, cell RNA and protein samples were collected, and GAPDH was used as an internal control to detect the expression of PI3K, AKT, mTOR genes and their phosphorylated proteins by quantitative reverse transcription polymerase chain reaction (qRT-PCR) and Western blot (WB). The primer sequence is shown in Table 1.

Table 1 Primer sequences of GAPDH, PI3K, AKT, and mTOR genes

Drug-regulated animal experiments

A total of 9 two-month-old Sprague–Dawley (SD) rats which were SPF grade rats, half male, and half female and weighed 200 ± 20 g, were purchased from Hubei Provincial Laboratory Animal Research Center, license number SCXK (E) 2020–0018. They were fed at the Experimental Animal Center of Xizang Minzu University. Animal experiments complied with ethics and principles.

SD rats were randomly divided into the control group and the experimental group (3 rats in each group) according to the random number table method. The control group was given the same amount of normal saline every day for 7 consecutive days. The experimental group was given IPRN (40 mg/kg/time) continuously for 7 days. At 24 h after the last gavage, the knee joint cartilage tissues of rats in the two groups were collected, and RNA was extracted. Rat β-actin was used as an internal reference, and the relative expression levels of the PI3K gene were calculated according to 2−∆∆Ct.

Statistical analysis

The control group and the experimental group were repeated three times each group, and the statistical analysis was carried out in turn, and finally, the effect of IPRN on MC3T3-E1 cells and SD rats was obtained. The mean ± standard deviation (SD) of all experimental data was taken. Statistical analyses were performed by GraphPad Prism 8.0 statistical software. The data were compared between the two groups by Student’s t test, and P < 0.05 was considered to be statistically significant.


Target search and network construction of IPRN

As shown in Fig. 2, we found 87 IPRN-related target proteins and constructed a network diagram of the drug and target proteins by Cytoscape software. It can be seen from the figure that the purple inverted triangle on the left represents IPRN, and the rectangles, ellipses, and parallelograms radiating from the IPRN node represent the target proteins.

Fig. 2
figure 2

Network map of isopsoralen-related target proteins

Visualization of OP-target network

The number of OP-related genes retrieved from the PubMed, GeneCards, and OMIM databases was 317, 4066, and 34, respectively. The results from the three databases were exported and then merged in Excel to remove duplicates and take the intersected genes, and finally, 242 genes associated with OP were obtained. As shown in Fig. 3, the red triangle represents OP, and the rounded rectangles and ellipses radiating from the OP node represent OP-related genes.

Fig. 3
figure 3

Network diagram OP targets genes

Integration of the key IPRN targets in the treatment of OP

A total of 19 common target genes were obtained from the Venn diagram (Fig. 4). These target proteins were input into the STRING online software for analysis and retrieval (Fig. 5). The target protein FDPS didn't interact with the other 18 cross-targets, so it was deleted. The results were imported into Cytoscape software for visualization (Fig. 6), and the degree values of protein interaction network nodes were analyzed in ascending order (Fig. 7). The higher the degree value reaches, the longer the histogram is, indicating that the gene is at the core of the network. As can be seen in Fig. 6, AKT serum/triamino kinase 1 (AKT1) is the core of the network. Moreover, phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) and mechanistic target of rapamycin (mTOR) are closely related to AKT1. Therefore, AKT1, PIK3CA, and mTOR may be the key targets of IPRN in the treatment of OP, and other key overlapping target proteins are interleukin 6 (IL6), Runt-related transcription factor 2 (RUNX2), bone gamma-carboxyglutamate protein (BGLAP), SRY-box 9 (SOX9), collagen type I alpha 1 chain (COL1A1), estrogen receptor alpha (ESR1), tumor necrosis factor (TNF) and albumin (ALB). These targets may show important biological effects of IPRN in the treatment of OP. Details of the above target proteins are presented in Table 2.

Fig. 4
figure 4

Intersection target Venn diagram

Fig. 5
figure 5

PPI network diagram of cross-targets

Fig. 6
figure 6

Network construction diagram of overlapping target proteins of drugs and diseases

Fig. 7
figure 7

Degree values of the targets

Table 2 Basic information on overlapping target proteins of drugs and diseases

GO functional annotation scores and KEGG enrichment analysis of common targets

As shown in Fig. 8, the functional distribution of 18 target proteins was explored by GO functional analysis. The GO analytical results of the predicted key IPRN targets for OP were shown in this figure. The results of this study showed that these targets can interact with cytokine receptors, β-catenin, DNA-binding transcription factors, cell proliferation, etc. These findings further indicated that IPRN may cause OP by participating in these biological processes.

Fig. 8
figure 8

GO enrichment analysis of isopsoralen

We obtained the results of the KEGG enrichment analysis in Figs. 9 and 10. From these results, we concluded that the anti-OP IPRN is mainly involved in the PI3K-Akt signaling pathway, AGE-RAGE signaling pathway in diabetic complications, proteoglycans in cancer, TNF signaling pathway, and other biological pathways.

Fig. 9
figure 9

Histogram of KEGG functional enrichment analysis

Fig. 10
figure 10

Key signaling pathways of overlapping drug-disease targets

Molecular docking

According to the final prediction results of network pharmacology and literature review, first, from the RSCB PDB database ( download the AKT1 protein structure (PDB ID: 3OS5), PIK3CA protein structure (PDB ID: 4JPS) and mTOR protein structure (PDB ID: 4O75). Then we use AutoDock Vina 1.1.2 software to dock these target proteins (AKT1, PIK3CA and mTOR) with IPRN. The larger the absolute value of the Docking Score is, the stronger the bond between the compound and the target becomes, and the more stable the molecular conformation becomes. Therefore, we concluded that PIK3CA has the strongest binding force to IPRN, and its molecular structure is the most stable one (Fig. 11). The second is AKT1 (Fig. 12A), and the weakest one is mTOR (Fig. 12B). The binding energy of IPRN and the molecular docking of related targets are shown in Table 3.

Fig. 11
figure 11

Schematic diagram of the molecular docking between isopsoralen and its target (PI3KCA)

Fig. 12
figure 12

Schematic diagram of the molecular docking between isopsoralen and its target (A: AKT1; B: mTOR)

Table 3 Molecular docking of Isopsoralen with related targets

Molecular dynamics simulation

MD simulation provides insight into the stability of protein–ligand complexes. In this study, based on the docking results, PI3K(4JPS)-IPRN, AKT(3OS5)-IPRN, and mTOR(4O75)- IPRN were subjected to molecular dynamics simulation analyses for 100 ns to evaluate the motion, trajectory, structural features, binding potential, and conformational changes of molecules. The root means square deviation (RMSD) is a good measure of the conformational stability of protein and ligands and is a measure of the extent of deviation in the position of atoms from the starting position. A lower deviation indicates better conformational stability. The changes in RMSD values of the complexes were analyzed in this study. As shown in Fig. 13A, the RMSD of the AKT-IPRN composite system reached equilibrium after 10 ns. In the same way, the RMSD of the mTOR- IPRN composite system reached equilibrium after 7 ns (Fig. 13B), and the fluctuation value was small during the whole simulation process; Although the RMSD of PI3K- IPRN composite system fluctuated, it tended to be stable after 20 ns (Fig. 13C). This shows that the conformation of the protein does not change significantly after the combination of a small molecular ligand and protein, and the combination is relatively stable.

Fig. 13
figure 13

Root-mean-square deviation (RMSD) plots of AKT(3OS5)-IPRN (A), mTOR(4O75)-IPRN (B), PI3K(4JPS)-IPRN (C) during of molecular dynamics simulation

Results of molecular experimental validation

Cell experiments

In the early stage of this experiment, by exploring the administration time, 24, 48, and 72 h were selected, and it was found that the difference in gene expression was more obvious when the drug was administered at 48 h. Therefore, 48 h was chosen as the experimental condition. Figure 14 shows the gene expression of the MC3T3-E1 cell line at 48 h and different drug concentrations. As can be seen from the figure, compared with the control group, 10, 20, and 50 μM concentrations in the 48 h group can significantly increase the expression level of PI3K, but with the increase of the dose, the effect of promoting the expression of PI3K does not increase significantly, and even appears downward trend (P < 0.001). In addition, compared with the control group, IPRN at concentrations of 10, 20, and 50 μM in the 48 h group could significantly increase the expression levels of AKT and mTOR, and it was positively correlated with the dose (P ˂ 0.05). More importantly, WB results showed that the expression of PI3K/p-PI3K, AKT/p-AKT, and mTOR/p-mTOR was significantly increased when 20 μM IPRN was added compared with the control group (Fig. 15A–C). In a word, the above results suggest that IPRN may affect the growth of MC3T3-E1 by participating in the PI3K-AKT-mTOR signaling pathway.

Fig. 14
figure 14

The expression of PI3K (A), AKT (B), mTOR (C) genes. *P < 0.05, **P < 0.01, ***P < 0.001

Fig. 15
figure 15

The expression of PI3K/p-PI3K (A), AKT/p-AKT (B) and mTOR/p-mTOR (C). *P < 0.05, **P < 0.01, ***P < 0.001

Animal experiments

Knee articular cartilage tissues were collected from rats in the experimental group and the control group, and qRT-PCR was used to detect related genes, and then the Ct values and expression levels of β-actin and PI3K in each group were obtained (Fig. 16). After the calculation of 2−ΔΔCt, they conformed to the normal distribution. The results showed that 40 mg/kg/time IPRN could increase the expression level of the PI3K gene in chondrocytes of SD rats after gavage for 7 days, which was significantly higher in the experimental group than that in the control group (P < 0.01). This result suggested that IPRN may exert an anti-OP effect by activating the PI3K gene.

Fig. 16
figure 16

The expression of PI3K and β-actin. *P < 0.05, **P < 0.01, ***P < 0.001


OP is a degenerative disease characterized by decreased bone mass, degeneration of bone microarchitecture, and increased risk of bone fragility, which results in fractures and seriously endangers the health of middle-aged and elderly people [27]. At this stage, the side effects of Western medicine in the treatment of OP have caused considerable controversy and therefore, the researchers pay more attention to TCM which emphasizes the regulation of systemic function. As a medicinal ingredient, IPRN has potential advantages in the multi-targeted treatment of OP. Therefore, based on network pharmacology, this study systematically explored the pharmacological effects and potential targets of IPRN and discussed the potential molecular mechanism of IPRN in the treatment and prevention of OP through molecular experiments. The results of this study demonstrated that the anti-OP effect of IPRN is achieved through multiple targets and multiple pathways.

As we have learned, the bone tissue of the human body will be continuously updated with the progression of osteoblast-mediated bone formation and osteoclast-mediated bone resorption. In old age, bone resorption is significantly greater than formation, which causes bone loss. Previous studies have shown that IPRN can promote osteoblast differentiation and mineralization, increase calcium nodule levels and alkaline phosphatase (ALP) activity, and upregulate osteoblast markers including ALP, RUNX2, and COL1A1 [21]. Furthermore, IPRN can also alleviate oxidative stress-induced osteoblast injury through the Wnt/β-catenin signaling pathway [28]. The latest literature pointed out that the eight best differential genes, including AKT1, have a high diagnostic effect on the risk of OP [29]. In addition, tumor necrosis factor-α (TNF-α) can promote osteoclast differentiation to enhance bone resorption and inhibit osteoblast differentiation to impair bone formation [30]. In this study, a network diagram of IPRN-OP interacting proteins was constructed, with a total of 18 overlapping genes, including RUNX2, COL1A1, AKT1, and TNF. On the one hand, the results reflected that the best interaction point of IPRN was found by network pharmacology. On the other hand, it has been confirmed that IPRN can promote the expression of PI3K, AKT, and mTOR genes in MC3T3-E1 cells after 48 h. And in vivo experiments have also confirmed that 40 mg/kg/time IPRN can promote the PI3K gene expression in SD rat chondrocytes. These results expanded the idea of studying the specific pathogenesis of OP.

From another perspective, the anti-OP mechanism of prescriptions of TCM involves multiple regulatory pathways, such as PI3K/AKT/m TOR [31], Wnt/β-catenin [32], BMP/Smad [33], MAPK pathway [34], and so on. Based on the KEGG pathway enrichment analysis, we found that PIK3CA, AKT1, and mTOR were all enriched in the PI3K/AKT/mTOR signaling pathway. According to the literature, this pathway plays a crucial role in the pathological process of OP, and the activation of the PI3K/AKT/mTOR pathway can promote the differentiation of bone component cells including osteoblasts, chondrocytes, myoblasts, and adipocytes [35,36,37]; Naringenin exerts a protective effect on glucocorticoid-induced OP by regulating the PI3K/AKT/mTOR signaling pathway [38]. IGF-1 stimulates osteoblast precursors to form osteoblasts through the PI3K/AKT/mTOR pathway and promotes osteogenic differentiation [39]. In this study, the key genes of this pathway were verified in the follow-up experiments, and it was further concluded that IPRN may exert its anti-OP effect through the PI3K/AKT/mTOR pathway. In addition, the literature pointed out that the functional enrichment analysis of the targets of Liuwei Dihuang Pills in diabetic nephropathy-related OP found that the AGE-RAGE signaling pathway in diabetic complications was related to the main targets [40]. Similarly, this study enriched the same pathway by IPRN. The results suggested that IPRN has a positive effect in the treatment of diabetic nephropathy-related OP and can be used as a therapeutic drug.

Although this study has obtained the above results, it still has its shortcomings. First of all, the research on animal models of OP with TCM preparations is mainly to find new drugs or to further prove efficacy. Animal experiments in this study are relatively simple, just to verify whether the effect of IPRN on OP is obvious. Human OP is divided into primary and secondary, and animal models may not be completely consistent with the disease, so it is still necessary to improve and refine this aspect in subsequent research. Second, the verification of relevant protein expression is lacking in animal experiments. In the future, we are going to add WB, immunohistochemistry (IHC), and other methods to detect the expression of proteins. Despite the above shortcomings, our study still found that IPRN can make an anti-OP difference, which further complements the molecular mechanism of TCM in the treatment of OP.


Based on network pharmacology, this study predicted the target genes of IPRN in the treatment of OP and preliminarily verified by in vitro and in vivo experiments that IPRN exerts an anti-OP effect by regulating the PI3K/AKT/mTOR pathway.


  1. Möckel L, Bartneck M, Möckel C. Risk of falls in postmenopausal women treated with romosozumab: Preliminary indices from a meta-analysis of randomized, controlled trials. Osteoporosis Sarcopenia. 2020;6:20–6.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Migliorini F, Colarossi G, Baroncini A, Eschweiler J, Tingart M, Maffulli N. Pharmacological management of postmenopausal osteoporosis: a level I evidence based—expert opinion. Expert Rev Clin Pharmacol. 2021;14:105–19.

    Article  CAS  PubMed  Google Scholar 

  3. Migliorini F, Giorgino R, Hildebrand F, Spiezia F, Peretti GM, Alessandri-Bonetti M, Eschweiler J, Maffulli N. Fragility fractures: risk factors and management in the elderly. Medicina (Kaunas). 2021;57:1119.

    Article  PubMed  Google Scholar 

  4. Gennari C, Agnusdei D. Calcitonin, estrogens and the bone. J Steroid Biochem Mol Biol. 1990;37:451–5.

    Article  CAS  PubMed  Google Scholar 

  5. Moschonis G, Katsaroli I, Lyritis GP, Manios Y. The effects of a 30-month dietary intervention on bone mineral density: the Postmenopausal Health Study. Br J Nutr. 2010;104:100–7.

    Article  CAS  PubMed  Google Scholar 

  6. Migliorini F, Colarossi G, Eschweiler J, Oliva F, Driessen A, Maffulli N. Antiresorptive treatments for corticosteroid-induced osteoporosis: a Bayesian network meta-analysis. Br Med Bull. 2022;143:46–56.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Oumer KS, Liu Y, Yu Q, Wu F, Yang S. Awareness of osteoporosis among 368 residents in China: a cross-sectional study. BMC Musculoskelet Disord. 2020;21:197.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Wang L, Yu W, Yin X, Cui L, Tang S, Jiang N, Cui L, Zhao N, Lin Q, Chen L, Lin H, Jin X, Dong Z, Ren Z, Hou Z, Zhang Y, Zhong J, Cai S, Liu Y, Meng R, Deng Y, Ding X, Ma J, Xie Z, Shen L, Wu W, Zhang M, Ying Q, Zeng Y, Dong J, Cummings SR, Li Z, Xia W. Prevalence of osteoporosis and fracture in China: the China osteoporosis prevalence study. JAMA Netw Open. 2021;4:e2121106.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Migliorini F, Maffulli N, Colarossi G, Eschweiler J, Tingart M, Betsch M. Effect of drugs on bone mineral density in postmenopausal osteoporosis: a Bayesian network meta-analysis. J Orthop Surg Res. 2021;16:533.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Osteoporosis prevention, diagnosis, and therapy. Jama. 2001;(285):785–95.

  11. Migliorini F, Maffulli N, Spiezia F, Peretti GM, Tingart M, Giorgino R. Potential of biomarkers during pharmacological therapy setting for postmenopausal osteoporosis: a systematic review. J Orthop Surg Res. 2021;16:351.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Migliorini F, Maffulli N, Spiezia F, Tingart M, Maria PG, Riccardo G. Biomarkers as therapy monitoring for postmenopausal osteoporosis: a systematic review. J Orthop Surg Res. 2021;16:318.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Li H, Xiao Z, Quarles LD, Li W. Osteoporosis: mechanism, molecular target and current status on drug development. Curr Med Chem. 2021;28:1489–507.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Rossouw JE, Anderson GL, Prentice RL, LaCroix AZ, Kooperberg C, Stefanick ML, Jackson RD, Beresford SA, Howard BV, Johnson KC, Kotchen JM, Ockene J. Risks and benefits of estrogen plus progestin in healthy postmenopausal women: principal results from the Women’s Health Initiative randomized controlled trial. JAMA. 2002;288:321–33.

    Article  CAS  PubMed  Google Scholar 

  15. Indran IR, Liang RL, Min TE, Yong EL. Preclinical studies and clinical evaluation of compounds from the genus Epimedium for osteoporosis and bone health. Pharmacol Ther. 2016;162:188–205.

    Article  CAS  PubMed  Google Scholar 

  16. Lim DW, Kim YT. Anti-osteoporotic effects of Angelica sinensis (Oliv.) Diels extract on ovariectomized rats and its oral toxicity in rats. Nutrients. 2014;6:4362–72.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Pan Y, Niu Y, Li C, Zhai Y, Zhang R, Guo X, Mei Q. Du-zhong (Eucommia ulmoides) prevents disuse-induced osteoporosis in hind limb suspension rats. Am J Chin Med. 2014;42:143–55.

    Article  PubMed  Google Scholar 

  18. Zhang X, Zhao W, Wang Y, Lu J, Chen X. The chemical constituents and bioactivities of Psoralea corylifolia Linn.: a review. Am J Chin Med. 2016;44:35–60.

    Article  PubMed  Google Scholar 

  19. Wang Y, Hong C, Zhou C, Xu D, Qu HB. Screening antitumor compounds psoralen and isopsoralen from Psoralea corylifolia L. seeds. Evid Based Complement Alternat Med. 2011;2011:363052.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Lu H, Zhang L, Liu D, Tang P, Song F. Isolation and purification of psoralen and isopsoralen and their efficacy and safety in the treatment of osteosarcoma in nude rats. Afr Health Sci. 2014;14:641–7.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Ge L, Cui Y, Cheng K, Han J. Isopsoralen enhanced osteogenesis by targeting AhR/ERα. Molecules. 2018;23:2600.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Tang DZ, Yang F, Yang Z, Huang J, Shi Q, Chen D, Wang YJ. Psoralen stimulates osteoblast differentiation through activation of BMP signaling. Biochem Biophys Res Commun. 2011;405:256–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Hopkins AL. Network pharmacology: the next paradigm in drug discovery. Nat Chem Biol. 2008;4:682–90.

    Article  CAS  PubMed  Google Scholar 

  24. Szklarczyk D, Gable AL, Nastou KC, Lyon D, Kirsch R, Pyysalo S, Doncheva NT, Legeay M, Fang T, Bork P, Jensen LJ, von Mering C. The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021;49:D605–12.

    Article  CAS  PubMed  Google Scholar 

  25. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Article  CAS  PubMed  Google Scholar 

  26. Su G, Morris JH, Demchak B, Bader GD. Biological network exploration with Cytoscape 3. Curr Protoc Bioinformat. 2014;47:8–13.

    Article  Google Scholar 

  27. Ye G, Huang Y, Yin L, Wang J, Huang X, Bin X. Association between LEPR polymorphism and susceptibility of osteoporosis in Chinese Mulao people. Artif Cells Nanomed Biotechnol. 2022;50:10–6.

    Article  CAS  PubMed  Google Scholar 

  28. Li YP, Wu B, Liang J, Li F. Isopsoralen ameliorates H2O2-induced damage in osteoblasts via activating the Wnt/β-catenin pathway. Exp Ther Med. 2019;18:1899–906.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Tang H, Han Q, Yin Y. Screening of Important markers in peripheral blood mononuclear cells to predict female osteoporosis risk using LASSO regression algorithm and SVM method. Evol Bioinform Online. 2022;18:11769343221075014.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Liu Z, Li C, Huang P, Hu F, Jiang M, Xu X, Li B, Deng L, Ye T, Guo L. CircHmbox1 targeting miRNA-1247-5p is involved in the regulation of bone metabolism by TNF-α in postmenopausal osteoporosis. Front Cell Dev Biol. 2020;8:594785.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Cai P, Lu Y, Yin Z, Wang X, Zhou X, Li Z. Baicalein ameliorates osteoporosis via AKT/FOXO1 signaling. Aging. 2021;13:17370–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Wang X, Tian Y, Liang X, Yin C, Huai Y, Zhao Y, Huang Q, Chu X, Wang W, Qian A. Bergamottin promotes osteoblast differentiation and bone formation via activating the Wnt/β-catenin signaling pathway. Food Funct. 2022;13:2913–24.

    Article  CAS  PubMed  Google Scholar 

  33. Lin GL, Hankenson KD. Integration of BMP, Wnt, and notch signaling pathways in osteoblast differentiation. J Cell Biochem. 2011;112:3491–501.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Wang Y, Cui K, Zhao H, Li D, Wang W, Zhu Y. Bushen Ningxin Decoction pharmacological serum promotes the proliferation and suppresses the apoptosis of murine osteoblasts through MAPK pathway. J Ethnopharmacol. 2009;122:221–6.

    Article  PubMed  Google Scholar 

  35. Zhou H, Jiao G, Dong M, Chi H, Wang H, Wu W, Liu H, Ren S, Kong M, Li C, Zhang L, Chen Y. Orthosilicic acid accelerates bone formation in human osteoblast-like cells through the PI3K-Akt-mTOR pathway. Biol Trace Elem Res. 2019;190:327–35.

    Article  CAS  PubMed  Google Scholar 

  36. Ghosh-Choudhury N, Abboud SL, Nishimura R, Celeste A, Mahimainathan L, Choudhury GG. Requirement of BMP-2-induced phosphatidylinositol 3-kinase and Akt serine/threonine kinase in osteoblast differentiation and Smad-dependent BMP-2 gene transcription. J Biol Chem. 2002;277:33361–8.

    Article  CAS  PubMed  Google Scholar 

  37. Kaliman P, Viñals F, Testar X, Palacín M, Zorzano A. Phosphatidylinositol 3-kinase inhibitors block differentiation of skeletal muscle cells. J Biol Chem. 1996;271:19146–51.

    Article  CAS  PubMed  Google Scholar 

  38. Ge X, Zhou G. Protective effects of naringin on glucocorticoid-induced osteoporosis through regulating the PI3K/Akt/mTOR signaling pathway. Am J Transl Res. 2021;13:6330–41.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. Bakker AD, Jaspers RT. IL-6 and IGF-1 signaling within and between muscle and bone: How important is the mTOR pathway for bone metabolism? Curr Osteoporos Rep. 2015;13:131–9.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Liu MM, Lv NN, Geng R, Hua Z, Ma Y, Huang GC, Cheng J, Xu HY. Uncovering the key miRNAs and targets of the Liuwei Dihuang pill in diabetic nephropathy-related osteoporosis based on weighted gene co-expression network and network pharmacology analysis. Endocr Metab Immune Disord Drug Targets. 2021;22:274–89.

    Article  Google Scholar 

Download references


We thank all participants and volunteers in this study. We also thank Inner Mongolia People's Hospital for helping with sample collection.


This research was supported by the Medical and Health Science and Technology Plan Project of the Inner Mongolia Autonomous Region Health Commission (202201004), a school-level joint project of Inner Mongolia Medical University (YKD2021LH039), and the Inner Mongolia Autonomous Region Science and Technology Plan Project (2022YFSH0006).

Author information

Authors and Affiliations



Jian Wang and Tianyu Chen drafted the work or revised it critically for important content; Xiang Li and Yu Zhang performed the experiments; Shuang Fu and Ruikun Huo analyzed the data and prepared figures and tables; Yan Duan conceived and designed the experiments. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Yan Duan.

Ethics declarations

Ethics approval

All animal experiments are conducted in accordance with the ethical standards of Inner Mongolia People’s Hospital.

Competing interests

The authors declare that they have no known competing financial interests or personal relationships that could influence the work of this study.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, J., Chen, T., Li, X. et al. A study on the anti-osteoporosis mechanism of isopsoralen based on network pharmacology and molecular experiments. J Orthop Surg Res 18, 304 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Isopsoralen
  • Osteoporosis
  • Network pharmacology
  • PI3K/AKT/mTOR pathway
  • Molecular docking