Abstract

Objective. This study aimed to explore the therapeutic targets and related pathways of Kuntai capsules (KTCs) for premature ovarian insufficiency (POI) using network pharmacology and molecular docking. Methods. The active components and their targets of KTCs were retrieved from the Traditional Chinese Medicine Systems Pharmacology (TCMSP) website, and disease therapeutic targets of POI were downloaded from DisGeNET, GeneCards, and OMIM databases and combined with the disease differential genes of POI microarray dataset from the Gene Expression Omnibus (GEO) database. The intersecting genes of drug potential therapeutic targets and disease therapeutic targets were uploaded to the STRING database to form a protein-protein interaction network. Also, the possible pathway of KTCs in the treatment of POI was explored by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis; through the core therapeutic targets, the corresponding active ingredients of KTCs were found. Finally molecular docking was conducted to verify the accuracy of the drug action. Results. 120 potential therapeutic targets of KTCs for POI were found. The bioinformatics analysis revealed that KTCs may regulate the recruitment, growth, and development of follicles by controlling various pathways such as fluid shear stress, atherosclerosis, PI3K/AKT, and p53 signaling. They can inhibit granulosa cell atrophy and apoptosis, promote follicle maturation, regulate follicle sensitivity to follicle-stimulating hormone, and ultimately impact ovarian function. The core therapeutic targets were TP53, AKT1, and TNF, and the corresponding active ingredients were quercetin, kaempferol, and stigmasterol. The molecular docking results showed that all the root mean square deviations were less than 2. Conclusions. KTCs improve ovarian function probably by acting on regulating the recruitment of follicles, reducing the apoptosis of granulosa cells, promoting their growth and development, reducing oxidative stress damage, and improving the sensitivity of follicles to FSH.

1. Introduction

Premature ovarian insufficiency (POI) is ovarian failure in women prior to 40 years old. This generally manifests as menstrual abnormalities (amenorrhea, sparse, or frequent menstruation), high follicle-stimulating hormone levels (FSH >25 U/L), and decreased or volatile estrogen levels [1]. POI may be accompanied by a range of symptoms, including hot flushes, sweating, and irritability. POI is the depletion and quality decrease of ovarian follicles, leading to infertility and other complications, such as perimenopausal symptoms that can seriously affect a woman’s quality of life. Patients with reproductive needs often turn to assisted reproductive technology, but the results are often unsatisfactory. Kuntai capsule (KTC) is a patented Traditional Chinese Medicine (TCM) composed of Rehmanniae Radix Praeparata, Paeoniae Radix Alba, Coptidis Rhizoma, Asini Corii Colla, Scutellariae Radix, and Poria. In TCM terms, they have the effect of nourishing yin, clearing heat, transpiring nerves, and getting rid of annoyance. KTCs are recommended for various syndromes including both pre- and postmenopause. Previous clinical trials reported that these capsules can improve ovarian function and associated reproductive outcomes [24]. However, owing to the multiple components, targets, and pathways, it has not been able to fully determine the mechanism of KTCs in POI treatment through experimental research. Recently, the rise of network pharmacology approaches and molecular docking techniques has brought light to explore the mechanism of action of TCM [5, 6]. Based on the previously proposed TCM pathogenesis of POI, the potential therapeutic targets of KTCs in the treatment of POI were identified and validated using network pharmacology. In addition, GO and KEGG analyses were performed to predict the probable pathway, and the main components of KTCs were also molecularly docked with the core targets of POI, hoping to enrich the clinical application and translation of KTCs in POI. The workflow of the study is shown in Figure 1.

2. Material

The following databases and tools were utilized for the analysis: Traditional Chinese Medicine Systematic Pharmacology Database and Analysis Platform (TCMSP) (https://old.tcmsp-e.com/tcmsp.php), UniProt Database (https://www.uniprot.org), STRING Database (https://cn.string-db.org), GeneCards Database (https://www.genecards.org), OMIM Database (https://www.omim.org), DisGeNET Database (https://www.disgenet.org), GEO Database (https://www.ncbi.nlm.nih.gov/geo), Protein Data Bank (PDB, https://www.rcsb.org/), Venny 2.0 (https://bioinfogp.cnb.csic.es/tools/venny/index2.0.2.html), PubChem website (https://pubchem.ncbi.nlm.nih.gov/), as well as software tools like Cytoscape 3.7.2, PyMOL 2.5.2, AutoDockTool-1.5.6, Open Babel 2.4.1, Discovery Studio 2019 Client, Adobe Illustrator CS6, and R3.6.2.

3. Methods

3.1. Screening the Active Ingredients and Targets of KTCs

The chemical constituents of KTCs were retrieved in TCMSP and screened according to absorption, distribution, metabolism, and excretion (ADME) parameters based on the conditions of drug-like properties ≥0.18 and bioavailability ≥40%. Then, we searched for the targets of the active ingredients through the same database and converted the target names to the gene symbols through the UniProt database. Finally, the Cytoscape 3.7.2 software was used to construct the active ingredients of KTCs and their corresponding target genes’ interaction network. The period of time searched was May 2023.

3.2. Collecting the Therapeutic Targets of POI

The therapeutic targets were attained by searching DisGeNET, GeneCards, and OMIM with “POI,” “premature ovarian insufficiency,” “POF,” and “premature ovarian failure” as keywords. The period of time searched was May 2023. In GeneCards database, the higher the relevance score, the closer relationship between the targets and POI. As a rule of thumb, those with score value greater than the median were set as potential therapeutic targets if there are too many genes. Duplicate ones from the three disease databases were deleted after merging. Then, we converted the target names to the gene names by the UniProt database.

3.3. Screening Differentially Expressed Genes Related to POI

GEO is a public genomic data repository of high throughout gene expression data and microarrays. POI expression data were retrieved from GEO using the search terms “premature ovarian insufficiency,” “primary ovarian insufficiency,” “POI,” and “premature ovarian failure.” The period of time searched was from the establishment of the GEO database to May 2023. Then the R 3.6.2 “limma” package was further used for differentially expressed genes (DEGs) screening (|log 2 (fold change)| = 1 and adjust ), and volcanoes and heatmaps of DEGs were mapped.

3.4. Potential Therapeutic Targets of KTCs in the Treatment of POI

Therapeutic targets related to POI obtained from DisGeNET, GeneCards, and OMIM were combined with DEGs from POI dataset retrieved from the GEO database. These targets were then intersected with KTC therapeutic targets by Venny 2.0 to identify prospective KTC therapeutic targets for POI.

3.5. The Analysis of Protein-Protein Interaction (PPI) Network, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG)

Potential therapeutic targets of KTCs for POI were uploaded to the STRING database to obtain their interactions by setting species as “Homo sapiens,” the minimum interaction threshold value to 0.95, and hiding disconnected nodes in the network. The PPI network was illustrated by Cytoscape 3.7.2, and the core therapeutic targets were further screened according to the degree value. To investigate the probable molecular mechanisms of KTCs in the treatment of POI, R3.6.2 “clusterProfiler” and “org.Hs.eg.db” packages were used to perform GO and KEGG enrichment analysis ( value cutoff = 0.05).

3.6. Identifying the Core Ingredients of KTCs

The venny 2.0 website was used to intersect the targets of each ingredient in KTCs with the potential therapeutic targets of KTCs for POI. The components were then ranked based on the number of intersecting targets, from highest to lowest. The first 3 ones with the most overlapping targets were identified as the core ingredients. Targets corresponding to the main core ingredient were uploaded to the STRING database to obtain the PPI interaction network by setting species as “Homo sapiens,” the minimum interaction threshold value to 0.95, and hiding disconnected nodes in the network. The core therapeutic targets were further screened according to the degree value.

3.7. Molecular Docking

Molecular docking was performed using CB Dock server adopting blind approach [7, 8]. In our study, the three-dimensional (3D) structure of the main core active ingredients was obtained from the PubChem website. Their chemical structures were downloaded in ∗ mol2 format, processed using Mgltools1.5.6, hydrogenated, had their charges calculated, combined with non-polar hydrogen, and saved in ∗ pdbqt format. The core targets of KTCs in treating POI were obtained from the Protein Data Bank (PDB), processed with Mgltools1.5.6 by hydrogenation, charge calculation, and addition of non-polar hydrogen. Subsequently, the active sites were defined (Table 1). AutoDock Vina 1.1.2 was used to dock small molecules with proteins, and the conformation with the highest score was selected for analysis and mapping by Discovery Studio 2019. The binding energy between each active ingredient and the receptor protein was obtained, a lower binding energy between the ligand and receptor in docking indicates a stronger molecular binding and a higher likelihood of interaction [9]. Finally, the visualization of molecular docking was optimized with the “Adobe Illustrator CS6” software.

4. Results

4.1. The Active Ingredients and Targets of KTCs

According to the screening conditions, 41 active ingredients and 207 therapeutic targets of KTCs were obtained through the TCMSP database (Supplementary Tables 1 and 2). The regulatory network of the active ingredients and their corresponding targets were drawn by Cytoscape 3.7.2 (Figure 2). The middle diamonds represent genes, the circles around represent drugs, and different colours represent different drugs. The larger the pattern area and the darker the colour, the more genes the drug acts on. The degree value represented the number of edges connected to the node in the graph. In terms of degree value, MOL000098, MOL000422, and A2 were the top three pharmaceutical ingredients, corresponding to quercetin (CID: 5280343), kaempferol (CID: 5280863), and stigmasterol (CID: 5280794).

4.2. Therapeutic Targets for POI

We found 299, 1264, and 1172 therapeutic targets in the DisGeNET, GeneCards (relevance score = 7.86), and OMIM databases, respectively, by using keywords mentioned above (Supplementary Table 3).

4.3. DEGs Related to POI in the GEO Dataset

The qualified chip expression dataset GSE135697 was obtained from the GEO database. This study included 10 patients with biologically abnormal POI who were undergoing IVF/ICSI-ET at the Reproductive Medicine Center of Shandong University and 10 normal controls [10]. Biologically abnormal POI is defined by the basal serum FSH of ≥10 IU/L, age of <40 years, regular menstruation, and unilateral ovarian antral follicle counts of <5. Differential analysis was performed on mRNAs with |log2 (fold change)| ≥ 1 and adjust value <0.05. The fold change refers to the comparison of the mean expression level of serum mRNAs in POI patients with normal patients. We detected 60 DEGs from this dataset (Supplementary Table 4). As shown in Figure 3(a), volcano plot of DEGs in GSE135697 from POI patients was plotted using ggplot2 package in R 3.6.2; fold change ≥2 and adjust value <0.05 mean that the gene expression is upregulated, indicated by red dots; fold change ≤ −2 and adjust mean that the gene expression is downregulated, indicated by blue dots; −2 < fold change <2 and adjust mean that the gene expression is stable expression, indicated by grey dots. Heatmap for cluster analysis of DEGs between the two groups is shown in Figure 3(b); it was plotted using pheatmap function in R3.6.2. Each column represents a sample, and each row represents a differentially expressed gene (DEG). Red indicates that the DEG was upregulated compared with another group (). Green indicates that the DEG was downregulated compared with another group (). Black indicates that the DEG was intermediate expression.

4.4. Potential Therapeutic Targets of KTCs

The genes obtained from GSE135697 were combined with POI-related targets from the DisGeNET, GeneCards, and OMIM databases. Subsequently, all these genes were intersected with therapeutic targets of KTCs, resulting in 120 potential therapeutic targets for POI (Figures 4(a) and 4(b); Supplementary Table 5).

4.5. PPI Network and GO and KEGG Analyses of Potential Therapeutic Targets for POI

These 120 potential therapeutic targets were submitted to the STRING database. The PPI network is shown in Figure 5(a). Then, Cytoscape 3.7.2 was used to optimize the interaction (Figure 5(b)). The degree value represents the number of edges connected to the node in the graph. The core targets have the largest area and the darkest colour. According to the degree value, AKT1, TP53, HSP90AA1, TNF, and RELA were found to be the core therapeutic targets, with UniProt id and RCSB PDB protein structure database protein numbers of AKT1 (P31749, 3OS5), TP53 (P04637, 3TG5), HSP90AA1 (P07900, 3BM9), TNF (P01375, 5WJJ), and RELA (Q04206, 6HL6). The GO enrichment of the potential targets is shown in Figure 6(a), the vertical axis represents the GO term name, the horizontal axis represents the gene ratio, the size of the dot indicates the number of genes expressed in the GO term, and the colour of the dot corresponds to the adjust value range. The GO enrichment revealed that the targets were mostly enriched in growth factor, ubiquitin protein ligase and protease bindings as well as nuclear and steroid hormone receptor activities, etc. Further, KEGG enrichment analysis reported that the targets were mostly enriched in PI3K-Akt and p53 signalling pathways as well as endocrine resistance (Figure 6(b)). The vertical axis represents the KEGG term name, the horizontal axis represents the number of genes expressed in the KEGG term, and the colour of the histogram corresponds to the adjust value range. The PPI network and the result of GO and KEGG analyses are shown in Supplementary Table 6.

4.6. The Core Ingredients of KTCs

The 148 quercetin corresponding therapeutic targets from the TCMSP database have 102 intersection genes with the 120 potential targets of KTCs for POI treatment (Figure 7(a)), which shows that quercetin is the core ingredient of KTCs. Quercetin corresponding targets were uploaded to the STRING database, and the optimized PPI interaction network by Cytoscape 3.7.2 suggested that the core therapeutic targets of quercetin were AKT1, TP53, TNF, VEGFA, and IL 6, which are highly similar to KTCs (Figure 7(b)). Then, according to the same steps, we found that kaempferol and stigmasterol are the next two main ingredients.

4.7. The Results of Molecular Docking

We performed molecular docking according to the steps described in Methods section and obtained the corresponding data (Supplementary Table 7). We found that the root mean square deviations (RMSDs) between the conformation of the docked ligand and the conformation of the original crystal structure of all docking results were less than 2. The docking data between the main ingredients and the core therapeutic targets of KTCs in POI treatment are shown in Table 1. Also, images of docking results for quercetin-TNF, kaempferol-TNF, stigmasterol-AKT1, and stigmasterol-TP53 are shown in Figure 8.

5. Discussion

The mammalian ovary is the primary female reproductive organ involved in oocyte maturation and in the synthesis and secretion of sex steroid hormones, estrogens, and progesterone, which are crucial for female fertility. The cortex region of the human ovary possesses a finite pool of primordial follicles, whose number is set before birth [11]. A variety of factors, such as genetic, immune, iatrogenic, and environmental ones, lead to damage or premature depletion of primordial follicles, which is often accompanied by a decline in follicle quality, and ultimately lead to POI [12].

Patients with POI are young and may have fertility requirements. But so far, there is no recognized effective way for these patients undergoing IVF-ET to obtain ideal pregnancy rate by autologous eggs. In order to improve the outcomes of assisted reproduction, adding adjuvant drugs for pretreatment is considered to be a common method. These drugs, including DHEA, oral contraceptives, growth hormone, coenzyme Q10, and immunomodulator, are either expensive, ineffective, or have serious adverse effects [13, 14]. Several studies have shown that TCM has several advantages in improving ovarian function. KTC is a Chinese patent medicine for the treatment of perimenopausal syndrome. Recent studies have reported that KTCs can improve ovarian stimulation and pregnancy outcomes in this group, but the mechanism remains unestablished [1518].

By intersecting the therapeutic targets of KTCs and POI in the common databases, 120 potential therapeutic targets of KTCs for POI were finally obtained in our research. Then, they were uploaded to the STRING database to construct a PPI network, and finally 5 core targets were found, including TP53, AKT1, TNF, HSP90AA1, and RELA. TP53, a transcription factor, is involved in cell cycle regulation that acts to negatively regulate cell division by controlling a set of genes required for this process. Research studies have shown that cytotoxic drugs could induce cytotoxicity in the maturing oocyte by activating autophagy and apoptosis in a caspase-dependent manner and could induce oxidative stress by generating reactive oxygen species that elevated the mutated ovarian TP53 protein [19]. As an important regulator of ovarian function, AKT1 is one of the serine/threonine-protein kinases and regulates many processes including metabolism, proliferation, cell survival, growth, and angiogenesis through serine and/or threonine phosphorylation of a range of downstream substrates [20]. TNF is a cytokine secreted by macrophages with a wide range of biological activities. TNF-α is related to the uptake of glucose in tissues, which may lead to the decline of female fertility [21]. Heat shock protein 90AA1 (HSP90AA1) is a molecular chaperone that aids in protein folding. Functional HSP90 operates as dimer and has intrinsic ATPase activity. Interestingly, HSP90AA1 has been implicated in cytokine production and is involved in normal reproductive processes, including estrogen receptor α regulation and luteolysis [22, 23]. RELA (v-rel avian reticuloendotheliosis viral oncogene homolog A), formerly known as nuclear factor kappa light chain polypeptide gene enhancer, is a member of the NF-κB family and is involved in many biological processes such as inflammation, immunity, differentiation, cell growth, tumorigenesis, and apoptosis.

Our research findings indicate that KTCs influence the fluid shear stress, PI3K-Akt, and p53 signaling pathways based on GO and KEGG enrichment analysis. Previous studies have found that both oocytes and granulosa cells have PI3K/Akt signalling pathways. These are involved in regulating oocyte growth and primordial follicle development and play an important role in follicular recruitment [24, 25]. Follicle-stimulating hormones bind to specific receptors on the membrane of ovarian granulosa cells in order to activate upstream protein kinase A (PKA), which then activates the downstream PI3K/Akt pathway, thereby inducing aromatase expression, regulating the secretion of estrogen and progesterone, and promoting the maturation of ovarian granulosa cells [26, 27]. Abnormal regulation of the PI3K/Akt pathway can lead to premature recruitment of follicles, depleting the follicle pool [28] and resulting in premature ovarian failure. The activity of the p53 tumor suppressor protein can induce cell cycle arrest or apoptosis and is related to the apoptosis of granulosa cells. Several experiments have confirmed the close relationship between the p53 pathway and granulosa cell death [29]. Therefore, KTCs may regulate the recruitment, growth, and development of follicles by regulating PI3K/AKT and p53 signalling pathways, inhibit granulosa cell atresia and apoptosis, promote follicle maturation, regulate follicle sensitivity to follicle-stimulating hormone, and thus affect ovarian function.

By molecular docking, we found that quercetin-TNF, Kaempferol-TNF, stigmasterol-AKT1, and stigmasterol-TP53 have low docking energies which are less than −7 kcal/mol. The binding was better at binding energy ≤−5.0 kcal mol−1 and very good at binding energy ≤−7.5 kcal mol−1. Quercetin and kaempferol are flavonoids widely present in fruits, vegetables, and medicinal plants. They have attracted much attention due to their antioxidant, anti-inflammatory, and anticancer properties [30]. Quercetin works as an antioxidant by lowering free radical generation, preventing lipid peroxidation, and altering antioxidants, so as to protect the residual follicles of POI patients from oxidative stress damage and improve their ovarian function [31].

6. Conclusion

Using network pharmacology and molecular docking, the mechanism by which KTCs improve ovarian function in patients with POI was discovered. The core active ingredients of KTCs for POI were discovered to be quercetin, kaempferol, and stigmasterol, and its main therapeutic targets were AKT1, TP53, and TNF. Possible pathways related to endocrine resistance include the PI3K/AKT and p53 signaling pathways. They probably act by regulating the recruitment of follicles, reducing the apoptosis of granulosa cells, promoting their growth and development, reducing oxidative stress damage, and improving the sensitivity of follicles to FSH. However, there are some limitations. First, we only conducted bioinformatics mining, and our results should be verified in clinical samples. The second limitation is that only a small number of patients were included in the dataset GSE135697; we look forward to the release of more POI microarray data to refine our study.

Data Availability

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Ethical Approval

The data used in our study were obtained from public database; therefore, ethical approval was not required. No animals/humans were used in the studies that are the basis of this research.

Disclosure

The authors declare that all data were generated in-house and that no paper mill was used.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Ruihong Cai, Hailiang Wang, and Jinhua Zhang conceived and designed research. Jintuo Zhou edited images. Qiuping Lin analysed data. Ruihong Cai and Hailiang Wang wrote the manuscript. All authors have read and approved the manuscript.

Acknowledgments

The authors would like to thank https://www.enago.cn for providing manuscript insurance service.

Supplementary Materials

Supplementary Table 1: pharmacokinetic characteristics of active ingredients of KTCs. Supplementary Table 2: therapeutic targets corresponding to effective ingredients. Supplementary Table 3: therapeutic targets of POI disease database and GEO database. Supplementary Table 4: the DEGs related to POI in GSE135697. Supplementary Table 5: the intersection of KTCs’ therapeutic targets and POI-related databases. Supplementary Table 6: analysis of PPI network, GO, and KEGG of KTCs’ potential therapeutic targets. Supplementary Table 7: molecular docking results of TP53, AKT1, and TNF and their corresponding effective ingredients. (Supplementary Materials)