*Corresponding Author:
Y. Dong
Department of Traditional Chinese Medicine, College of Pharmacy and Biotechnology, Tianjin Medical College, Tianjin 300222, China
E-mail:
zxsd411@126.com
Date of Received 10 July 2021
Date of Revision 19 August 2022
Date of Acceptance 03 February 2023
Indian J Pharm Sci 2023;85(1):233-240  

This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 3.0 License, which allows others to remix, tweak, and build upon the work non-commercially, as long as the author is credited and the new creations are licensed under the identical terms

Abstract

Qing’E pills are a classical formula in traditional Chinese medicine that has the effect of tonifying the kidneys and strengthening the bones, first appearing in the tonic formulae in prescriptions of peaceful benevolent dispensary (10th century CE). To preliminarily elucidate the potential mechanisms of Qing’E pills in the treatment of postmenopausal osteoporosis based on network pharmacology. The principal components and corresponding protein targets of Qing’E pills were searched on traditional Chinese medicine systems pharmacology, bioinformatics analysis tool for molecular mechanism-traditional Chinese medicine, SymMap and encyclopedia of traditional Chinese medicine and the compound-target network was constructed by Cytoscape 3.9.1. The targets of postmenopausal osteoporosis were searched on online mendelian inheritance in man, GeneCards, histidine triad nucleotide binding protein and DisGeNET databases. The intersection of compound target and disease target was obtained and the coincidence target was imported into the search tool for retrieval of interacting genes database to construct a protein-protein interaction network. We further performed gene ontology and Kyoto encyclopedia of genes and genomes enrichment analyses on the targets. The 33 active compounds in Qing’E pills corresponded to 433 targets. 1647 targets were associated with postmenopausal osteoporosis disease. The active compounds with high degree values were quercetin, stearic acid and marine. The proteins with higher values in the protein-protein interaction network were mitogenactivated protein kinase 1, protein kinase B, phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha, phosphoinositide-3-kinase regulatory 1, amyloid beta precursor protein, Src proto-oncogene, nonreceptor tyrosine kinase, mitogen-activated protein kinase 8, heat shock protein 90 alpha family class A member 1, epidermal growth factor receptor, and Janus kinase 2. Gene ontology and Kyoto encyclopedia of genes and genomes analysis showed that Qing’E pills treatment of postmenopausal osteoporosis involves signaling pathways that mainly include the tumor necrosis factor signaling pathway, osteoclast differentiation, Wnt signaling pathway and estrogen signaling pathway. We revealed the active ingredients and potential molecular mechanism of Qing’E pills in the treatment of postmenopausal osteoporosis and provided a reference for subsequent basic research

Keywords

Qing’E pills, network pharmacology, traditional Chinese medicine, postmenopausal osteoporosis

Postmenopausal Osteoporosis (PMO) is a metabolic disorder caused by estrogen deficiency and characterized by reduced Bone Mineral Density (BMD) and deterioration of bone microarchitecture[1,2]. Previous studies have shown a high incidence of PMO[3]. Patients with PMO have an increased susceptibility to hypo traumatic fractures, which are often accompanied by severe painful manifestations after fracture, which severely reduces the quality of life of patients. According to epidemiological studies, PMO affects half the women in the 6th and 7th decade of life. More than 57 % of women at risk in 27 European Union countries (EU27) do not receive bone-specific treatment[4,5]. As the population ages, aging comorbidities such as PMO pose a huge economic and health burden on society[6], and have become a growing public health problem[7]. The International Osteoporosis Foundation (IOF) and the European Federation of Pharmaceutical Industries Association (EFPIA) report that the cost of PMO treatment will increase by 25 % by 2025. The search for rational and effective treatments has become an important issue for clinicians to address.

In the late 1980s, estrogen replacement, perhaps calcitonin, along with calcium and vitamin D supplementation were used as the main treatment strategy for PMO[7]. In 2017, selective estrogen receptor modulators[8] (raloxifene), bisphosphonates[9] (alendronate, risedronate, ibandronate and zoledronic acid), human anti-Receptor Activator of Nuclear Factor Kappa-Β Ligand (RANKL) antibody[10] (ANKL; denosumab) and parathyroid hormone analog[11] (teriparatide) are included in the range of therapeutic agents. Despite significant advances in the development of drugs for the prevention and treatment of fractures, there are still major challenges in implementing treatment such as adverse effects, for example, osteonecrosis of the jaw and atypical femoral fractures have been reported following teriparatide treatment[3].

Qing’E Pills (QEP) is a mixture of Traditional Chinese Medicine (TCM) consisting of Eucommiae Cortex (Eucommia ulmoides (E. ulmoides), Psoraleae Fructus (Psoralea corylifolia), Juglandis Semen (Juglans regia (J. regia)) and Garlic Rhizoma[12-14]. However, the mechanism by which QEP improves the balance of bone metabolism in PMO patients is not known. As an emerging discipline, network pharmacology can accelerate the clinical translation of precise and effective therapeutic interventions[15]. The targets of the main active components of QEP and the targets of sepsis were collected by using the network pharmacology technology. The QEP component target network was built, as well as the Protein-Protein Interaction (PPI) network of chemical and disease interaction targets. The signaling pathways of biological function and efficacy were then investigated using Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis.

Materials and Methods

Data sources:

Procurement and screening of active compounds in QEP: The compounds and the corresponding targets of Eucommiae Cortex (E. ulmoides), Psoralea Fructus (Psoralea corylifolia), Juglandis Semen (J. regia) in QEP were obtained from the following databases, namely Traditional Chinese Medicine Systems Pharmacology (TCMSP) (tcmsp)[16], Bioinformatics Analysis Tool for Molecular Mechanism (BATMAN)-TCM (batman-tcm) [17], SymMap (http://www.symmap.org/)[18] and Encyclopedia of Traditional Chinese Medicine (ETCM) (http://www.tcmip.cn/ETCM/index.php/ Home/Index/)[19]. Oral Bioavailability (OB) refers to the amount of medicine that enters the circulatory system. Drug similarity (DL) is the similarity between a chemical and a known drug[16]. In the TCMSP system, OB ≥30 % and DL ≥0.18 were set as the filters for retrieving pharmacokinetic data.

Target prediction for the active compounds of QEP: The 2-dimensional structure files of the active ingredients were obtained by accessing the organic small molecule bioactivity database (https://pubchem.ncbi.nlm.nih.gov/, https:/pubchem.ncbi.nlm.nih.gov/)[20]. The structures of the active ingredients were cross-validated for ligand similarity and analytically aligned via the SwissTargetPrediction platform (http://www.swisstargetprediction.ch/)[21] to complement the TCMSP system with compound components not included and their corresponding targets of action. The active ingredient targets of the QEP were screened using the criterion of "probability value >0". The top 10 related targets (related target proteins) of each compound were selected and uploaded to the PharmMapper database (http://lilab.ecust.edu.cn/)[22]. The "Homo sapiens" was selected for the identification of potential target proteins (hereafter referred to as "targets").

Target prediction: Searching through GeneCards database (https://www.genecards.org/)[23], Online Mendelian Inheritance in Man (OMIM) database (https://omim.org/)[24], DisGeNET database (https://www.disgenet.org/)[25] for PMO-related target retrieval. The potential targets for PMO treatment with QEP were mapped through the ImageGP platform (http://www.ehbio.com/ImageGP/).

PPI network construction:

To further study the interaction targets, the protein identifier or protein sequence provided by Search Tool for Retrieval of Interacting Genes (STRING) (https://string-db.org/)[26] was used for analysis. Select “Homo sapiens” in the organism column and the minimum required interaction score is “highest confidence (0.9)”. Then input the correlation information between the targets into Cytoscape 3.7.2 for visualization and construct a PPI network.

Functional enrichment analysis of QEP'S therapeutic targets for PMO:

The Meta scape system (http://metascape.org/, version 3.0) was used to study the biological function of potential targets in PMO[27]. We collected terms with p<0.01, a minimum count of 3 and an enrichment factor >1.5. The Cytoscape 3.8.0 software package (Boston, MA, United States of America (USA)) was used to explain the complicated relationships among compounds, genes, pathways and diseases.

Results and Discussion

After removing duplications of the initial ingredients, a total of 33 initial chemical components of QEP were obtained through the TCMSP system and the SwissTarget system retrieval and manual retrieval. E. ulmoides Oliv. contains 25 active compounds, Psoralea corylifolia L contains 5 active compounds, J. regia L. and Allium sativum L include 2, 1 active compounds respectively. The 433 targets of 33 active ingredients were selected as common targets with PMO (1647 targets) (fig. 1) to construct a drug-active ingredient-target network for QEP for PMO (fig. 2). The active compounds with high QEP median values are shown in Table 1.

Type Name Degree Stress Radiality
Eucommia ulmoides Oliv. Quercetin 127 1394918 0.99
Psoralea corylifolia L Stearic acid 115 638198 0.99
Allium sativum L Marine 101 889700 0.99
Psoralea corylifolia L Xanthotoxin 81 1101216 0.99
Juglans regia L. Ellagic acid 58 781208 0.989
Eucommia ulmoides Oliv. Kaempferol 54 536792 0.99
Eucommia ulmoides Oliv. Beta-sitosterol 36 251990 0.99
Eucommia ulmoides Oliv. Cinchonan-9-al 30 204364 0.98
Eucommia ulmoides Oliv. Tabernemontanine 28 133726 0.99
Eucommia ulmoides Oliv. Beta-carotene 18 67074 0.98

Table 1: Active Compounds with High Moderate QEP Values

IJPS-Venn

Fig. 1: Venn diagram of QEP and PMO targets

IJPS-compounds

Fig. 2: Drug-active compounds-target diagram of QEP for PMO

We performed PPI network construction for the QEP-PMO common target based on Cytoscape software results (fig. 3). If the network contains between 3 and 500 proteins, the molecular complex detection algorithm was applied to identify densely connected network components. The nodes in fig. 3 represent target proteins and the edges represent two proteins connected, the larger the node, the larger the position in the network. The top 10 nodes are Mitogen-Activated Protein Kinase (MAPK) 1, Protein Kinase B (AKT1), Phosphatidylinositol-4, 5-Bisphosphate 3-Kinase Catalytic Subunit Alpha (PIK3CA), Phosphoinositide-3-Kinase Regulatory 1 (PIK3R1), Amyloid Beta Precursor Protein (APP), SRC Proto-Oncogene, Non-Receptor Tyrosine Kinase (SRC), MAPK8, Heat Shock Protein 90 Alpha Family Class A Member 1 (HSP90AA1), Epidermal Growth Factor Receptor (EGFR) and Janus Kinase 2 (JAK2), which are presumed to be important in the QEP treatment of PMO.

IJPS-relationship

Fig. 3: PPI network relationship diagram

During the GO enrichment analysis of the targets, 150 biological processes, 28 cellular components and 38 molecular functions were identified, indicating that the biological functions of QEP interfering with PMO were mainly in the areas of estrogen regulation, regulation of bone mineralization and regulation of factor binding. We screened 15 biological functions, cellular components and molecular functional pathways that had a significant impact in GO analysis for visualization (fig. 4), which related to positive regulation of transcription from nuclear receptor activity, transcription factor binding, protein domain specific binding, heme binding, protein homodimerization activity, response to hormone, response to inorganic substance, response to xenobiotic stimulus, cellular response to nitrogen compound, response to extracellular stimulus, membrane raft, transcription regulator complex, side of membrane and extracellular matrix.

IJPS-drug

Fig. 4: GO analysis of common drug-disease targets

In order to explore the signal pathway mechanism of QEP in the treatment of PMO, we performed KEGG enrichment analysis (fig. 5). We found that the targets were mainly related to 36 key pathways, including the calcium signaling pathway, the estrogen signaling pathway and the cyclic Guanosine Monophosphate-activates Protein Kinase G (cGMP-PKG) signaling pathway. At the same time the pathways in cancer, lipid and atherosclerosis, chemical carcinogenesis-receptor activation, prostate cancer, Phosphoinositide-3-Kinase (PI3K)-Akt signaling pathway, transcriptional misregulation in cancer, chemical carcinogenesis-reactive oxygen species, MAPK signaling pathway and tuberculosis were of great importance in the QEP efficacy research.

IJPS-analysis

Fig. 5: KEGG analysis of common drug-disease targets

As a recommended remedy in the PMO Chinese Medicine Guideline, QEP can promote bone calcium absorption and thus regulate the balance between calcium and phosphorus. Pharmacological studies have shown that QEP can inhibit osteoclast activity, maintain stable bone metabolism and regulate estrogen metabolism, thereby preventing and treating PMO. However, there is a lack of research on the mechanism of QEP efficacy. We investigate the potential mechanism of QEP in the treatment of PMO through a predictive study, with a view to providing directional options for the clinical management of PMO.

Thirty-three chemical components were collected from the TCMSP database, which can act in combination with other active components at multiple targets in the network. Six of the ten active ingredients with high degree values were from Eucommia tulipifera (E. tulipifera) and the associated components of E. tulipifera, such as quercetin and kaempferol, were connected to drug-disease common targets with high density, suggesting that E. tulipifera may play a good role in the treatment of PMO. The lignans and phenylpropanoids in Eucommia are immune-enhancing, and the sterol in bone marrow is used as a raw material for the production of vitamin D3, so it is assumed that QEP may prevent PMO exacerbation to some extent[28-34]. The topological analysis of the main compounds and the major intersecting targets revealed the highest binding of the target CYP3A4 and hirsutine. CYP3A4 (cytochrome P450 3A4) is a hemoglobin thioredoxin that completes a variety of oxidative reactions that affect bone loss[35]. Evellyn et al. showed that stimulation of the estrogen signaling pathway could regulate the induction of osteogenesis. Evellyn et al. also found that stimulation of the estrogen signaling pathway could regulate the body's secretion of estrogen. It suggested that QEP may improve bone loss due to estrogen reduction by modulating the calcium signaling pathway, estrogen signaling pathway and cGMP-PKG signaling pathway, which are the key effector pathways[36].

Among the high value targets, MAPK1 and MAKE8 are members of the MAPK signaling pathway, which intervenes in osteoblast apoptosis via the mitochondrial pathway[37-42]. AKT1 was a member of the AKT gene family[43]. AKT1 was found to have limited skeletal development in rats treated with AKT1 knockout, suggesting that AKTI may play a role in promoting osteogenesis[44]. AKT1, PIK3CA and PIK3R1 are members of the PI3K-Akt signaling pathway, through which AKTI plays an anti-inflammatory role in the metabolic process of the body[45-49]. AK2, a member of the JAK/signal transducer and activator of transcription signaling pathway, is involved in the development of bone destruction through inflammatory responses[50]. The key targets, including AKT1, PIK3CA and PIK3R1 are presumed by GO enrichment analysis to function mainly at the plasma membrane, membrane rafts and nucleus and to perform molecular functions such as kinase activity, protein binding and protein kinase binding through positive regulation of nitric oxide biosynthesis, signal transduction, protein phosphorylation and other biological processes. Although the cancer signaling pathway was ranked first in the enrichment analysis of this study, there are few studies on this pathway in relation to PMO and the mechanisms mediating it need to be investigated in subsequent studies.

There are still some limitations in this study. We only used modern bioinformatics methods to predict the role of QEP in PMO through network pharmacology. However, the current network information technology needs to be further improved and the accuracy and timeliness of the database data needs to be scientifically validated. Also, unidentified and undocumented compounds or targets may not be included in our analysis. Secondly, the bioactive components in this study are not fully representative of QEP. We believe that the pharmacodynamics and molecular biology experiments should be considered to further validate our results.

We conducted a preliminary exploration of the main chemical components and potential mechanisms of QEP for the treatment of PMO based on network pharmacology techniques. The study suggested that dulcimer may be the main drug for the treatment of PMO. Quercetin, stearic acid and marine may be the main active component of QEP. MAPK1, AKT1, PIK3CA, PIK3R1, APP, SRC, MAPK8, HSP90AA1, EGFR and JAK2 may be the potential therapeutic targets of QEP for PMO. Several signaling pathways such as tumor necrosis factor signaling pathway, osteoclast differentiation, Wnt signaling pathway and estrogen signaling pathway may be involved in the mechanism of QEP efficacy.

Conflict of interests:

The authors declared no conflict of interests.

References