Deciphering the Therapeutic Mechanisms of Wuzi Ershen Decoction in Treating Oligoasthenozoospermia through the Network Pharmacology Approach

Background Infertility affects approximately 15% of couples around the world, and male factors are accounted for 40–50%. Oligoasthenozoospermia is the most common reason for male infertility. Unfortunately, effective drug therapy is still lacking except for assisted reproductive technology (ART). Previous researchers found that Wuzi Ershen decoction (WZESD) can increase sperm count, enhance sperm vitality, and improve semen quality. However, the pharmacological mechanisms remain unclear. Methods In this study, we screened compounds and predicted the targets of WZESD based on the TCMSP and BATMAN-TCM database combined with literature searching in the PubMed database. We obtained proteins related to oligoasthenozoospermia through GeneCards and submitted them to STRING to obtain the protein-protein interaction (PPI) network. Potential targets of WZESD were mapped to the network, and the hub targets were screened by topology. We used online platform Metascape and Enrichr for GO and KEGG enrichment analyses. AutoDock Vina was utilized for further verification of the binding mode between compounds and targets. Results Totally, 276 bioactive compounds were obtained and targeted 681 proteins. 446 oligoasthenozoospermia disease-specific proteins were acquired, and further bioinformatics analysis found that they were mainly involved in the formation of gametes, meiosis, and sperm differentiation. Protein interaction network analysis revealed that target proteins of WZESD were associated with oligoasthenozoospermia disease-specific proteins. The 79 targets of disease-specific proteins, which were anchored by WZESD, mainly participate in the cellular response to the organic cyclic compound, regulation of the apoptotic process, nitricoxide biosynthetic and metabolic process, oxidative stress, and protein phosphorylation regulation, which are the causes for oligoasthenozoospermia. Molecular docking simulation further validated that bioactive compounds originated from WZESD with targeted proteins showed high binding efficiency. Conclusions This study uncovers the therapeutic mechanisms of WZESD for oligoasthenozoospermia treatment from the perspective of network pharmacology and may provide a valuable reference for further experimental research studies and clinical applications.


Introduction
Infertility, defined as the inability to achieve a clinical pregnancy within 12 months of regular unprotected intercourse, is estimated to affect 15% of all couples globally [1]. Male factor is believed to account for 50% of infertile couples. It has been exponentially increasing in recent years due to a comprehensive evaluation of reproductive male function and improvements in diagnostic tools [2]. More than 40% of infertile men are diagnosed with oligoasthenozoospermia, and their poor sperm quality is considered as one of the major causes of infertility [3]. Genetic abnormity (Klinefelter syndrome, microdeletion on the Y chromosome), endocrine disorders (primary or secondary hypogonadism), testicular dysfunction (heat, drug, or radiation therapy), and infections (mycoplasma and ureaplasma) are In this paper, an integrated approach including active ingredient screening, target prediction, network construction, and analysis with molecular docking was used to reveal potential drug targets related to oligoasthenozoospermia, active compounds from WZESD, and their pharmacological mechanisms of action for oligoasthenozoospermia therapy. e present work may provide a valuable reference for further experiment research studies and clinical applications of WZESD for oligoasthenozoospermia treatment.

Database Construction.
Candidate bioactive ingredients of 12 herbs in WZESD were screened from the Traditional Chinese Medicine Systems Pharmacology database (TCMSP, http://lsp.nwu.edu.cn/tcmsp.php) [25], the BAT-MAN-TCM [26,27] (http://bionet.ncpsb.org/batman-tcm/) database, and literature searching in the PubMed database. TCMSP database is a unique systems pharmacology platform designed for herbal medicines. e newly developed TCMSP provides up-to-date, quantitative, and system information about TCM ingredients, ADME-related properties, targets, and diseases [28]. e newest version of TCMSP comprises 510 effective herbal entries registered in the Chinese Pharmacopoeia with more than 31,000 ingredients, which spread over 18 different drug classes [29]. BATMAN-TCM database is the first online bioinformatics analysis tool specially designed for the research of molecular mechanisms of TCM based on TCM ingredients' target prediction and subsequent network pharmacology analyses [26]. e structures of these ingredients were saved as a mol2 format for further analysis. Chem3D Pro 14.0 was employed to optimize these molecules and minimize the energy.
Active ingredient screening including absorption, distribution, metabolism, and excretion (ADME) evaluations of drugs is crucial in drug development and discovery [30]. Because biological experiments are time-consuming and of high cost, the identification of ADME properties by computational methods has now become an inevitable choice in pharmaceutical research. Here, 2 ADME-related models, including the evaluation of oral bioavailability (OB) and drug-likeness (DL), were calculated to screen the potential bioactive ingredients of WZESD. OB represents the percentage of an orally administered dose of unchanged drug that reaches the systemic circulation and is one of the most important pharmacokinetic parameters [31]. Poor OB is usually the main reason for the failure of drug discovery, especially for TCM which in most cases are oral administration. DL index is a qualitative concept used in drug design for estimating the "druggability" of a substance calculated using the Tanimoto coefficient [32,33]. e evaluation of DL can facilitate screening for excellent compounds and increase the hit rate for the candidate drugs. For screening in the TCMSP database, OB ≥ 30% (a suggested criterion by the TCMSP database) and DL ≥ 0.18 (mean DL value for all DrugBank compounds) were regarded as the threshold for screening possible candidate drugs presently. At the same time, we also screened the BATMAN-TCM database and pieces of literature in the PubMed database with the Latin name of each herb for the candidate compounds.

Target Prediction.
Drug-target mappings were obtained from three sources. Experimental validated drug-target pairs were retrieved from the HIT database [34]. For those ingredients without validated targets, an in-house developed model SysDT was used which efficiently integrated the chemical, genomic, and pharmacological information for drug targeting and discovery with the random forest (RF) and support vector machine (SVM) algorithm [35]. is method displays the incredible performance of prediction for drug-target interactions, with a concordance of 82.83%, a sensitivity of 81.33%, and a specificity of 93.62%, respectively [28]. e information for the targets based on the above 2 methods was collected from the TCMSP database. A similarity-based method on the basis of the BATMAN-TCM database was employed to predict potential targets of TCM ingredients, the core idea of which was to rank potential drug-target interactions based on their similarity to the known drug-target interactions [26]. UniProt [36] (http:// www.uniprot.org/) was applied to obtain the official name of the predicted targets.

Oligoasthenozoospermia-Specific Protein Collection.
Information on oligoasthenozoospermia-associated target genes was collected from GeneCards: the human gene database [37] (http://www.genecards.org/, ver. 4.8.2), and only "Homo sapiens" proteins linked to oligoasthenozoospermia were selected. GeneCards is a searchable, integrative database that automatically integrates gene-centric data from ∼125 web sources, including genomic, transcriptomic, proteomic, clinical, and functional information. Furthermore, these target genes were submitted to STRING [38,39] (https://string-db.org/) to generate the protein interacting network. e STRING database aims to collect, score, and integrate all publicly available sources of protein-protein interaction information and to complement these with computational predictions and then achieve a comprehensive and objective global network, including direct (physical) as well as indirect (functional) interactions [40]. [41] was used in this study to evaluate the potential binding mode between bioactive compounds and putative targets; Discovery Studio 4.5 Client, BIOVIA [42], was applied to analyze the docked structures. e crystal structure of the target proteins of WZESD for oligoasthenozoospermia therapy was downloaded from the RCSB Protein Data Bank [43] (http://www. rcsb.org). e 3D chemical structures of bioactive ingredients were downloaded from the PubChem Compound database [44] or TCMSP database and submitted to minimize the energy by using the molecular mechanics-2 (MM2) force field in Chem3D Pro. e protein-ligand docking active site was defined by the location of the original ligand. Dimensions of the grids were set at 30 × 30 × 30 Ǻ in the x-, y-, and z-directions, with a spacing of 0.375 Ǻ between the grid points and the center placed at the active site of the original ligand crystallographic structures. All other docking and consequent scoring parameters used were kept at their default settings. e compound was regarded to be an efficient drug if the binding affinity was higher than that of the original ligand.

Network Construction and Analysis.
To characterize the "multicomponent", "multitarget" therapeutic mechanisms of herbal medicine for the treatment of oligoasthenozoospermia from a network target perspective, network construction was performed as follows: (1) Candidate compounds and candidate targets of WZESD were used to construct a candidate compound-candidate target (cC-cT) network (2) e PPI data obtained above were used to establish the oligoasthenozoospermia-specific protein interaction network (3) Potential compounds and putative targets from WZESD for oligoasthenozoospermia therapy were used to build a potential compound-potential target (pC-pT) network (4) Compounds and targets through molecular docking validation were used to construct a compound-target (C-T) network All networks were generated and analyzed by an open source of bioinformatics package for biological network analysis and visualization, Cytoscape 3.5.1 [45]. Two topological parameters, degree and betweenness centrality, were calculated for the obtained networks which imply the significance of a node. e top 10 compounds and targets based on the 2 parameters were considered as key components and core proteins.
2.6. Bioinformatics Analysis. Gene Ontology (GO) analysis is a major bioinformatics tool for annotating genes and gene products, and it aims to unify the representation of gene and gene product attributes across all species [46]. e Kyoto Encyclopedia of Genes and Genomes (KEGG) [47] is a collection of databases containing advanced functional information for the systematic analysis of gene functions, biological pathways, diseases, drugs, and chemical substances. For GO and KEGG pathway analyses, the Metascape online tool [48] (http://metascape.org) and "Enrichr" platform [49] (http://amp.pharm.mssm.edu/Enrichr/) were used to identify the predominant biological processes (BP) and KEGG pathways regulated by WZESD.

Results
ousands of years' clinical practices in herbal medicine have proven the in vivo efficacy and safety of herbal medicines [50], despite their mechanisms of action being generally unknown [51]. Fortunately, advances in systems biology and medicine have allowed the application of the systems pharmacology approach in the study of the herbal medicine. Unlike conventional medicine in which drugs are studied and used isolated, herbal medicine typically integrates several medicinal herbs which contain multiple chemical compounds that show more efficient therapeutic Evidence-Based Complementary and Alternative Medicine effects than an isolated single constituent [52]. In our work, a novel systems pharmacology approach integrated with polypharmacology and network biology was applied to uncover the therapeutic mechanisms of WZESD from a systematic level.

Bioactive Compound Screening and Target Prediction for WZESD.
A total of 276 candidate compounds were identified in WZESD, including 66 in SC, 10 in PS, 69 in LF, 15 in CS, 6 in RF, 77 in RS, 11 in FR, 3 in DRR, 5 in RRP, 14 in PR, 6 in OJ, and 23 in EH, respectively (Table S1). e 276 candidate compounds yield 681 candidate targets (Table S2), and the connections between them reach up to 5055. e candidate compounds and candidate targets were submitted to generate the cC-cT network (Figure 1), and details are depicted in the following part. Figure 1, the network consists of 998 nodes and 5055 edges, including 12 herbs (squares), 276 compounds (triangles), and 681 targets (squares). e numbers of candidate targets in SC, PS, LF, RF, CS, RS, FR, DRR, RRP, PR, OJ, and EH were 204, 162, 388, 215, 178, 238, 178, 20, 20, 134, 41, and 216, respectively (Table S3). Although the number of targets in each herb was different, they overlapped dramatically in the 12 herbs indicating that different ingredients in WZESD shared common or similar targets with synergistic effects. e duplicated target proteins were removed, and 681 unique terms were retained. Much of the compounds possessed more than one target, and a target was anchored by more than one compound, revealing the "multicompound", "multitarget" therapeutic mechanism of WZESD. Two centrality indicators, degree and betweenness, were calculated in this network. Different centralities reflect different importance of nodes in a network from different angles. e top 10 compounds and targets through network analysis are shown in Tables 1 and 2. Quercetin, L-asparagine, kaempferol, and stigmasterol were predicted as the major active compounds of WZESD. e proteins including PTGS2, NCOA2, PGR, PTGS1, and HSP90 were predicted as essential pharmacological proteins for the therapeutic effects of WZESD.

Analyses on the Oligoasthenozoospermia-Based Specific
Protein Interaction Network. 446 target genes related to oligoasthenozoospermia were searched in the GeneCards database (Table S4) and further submitted to the STRING database (34 proteins were not recognized) to construct the protein-protein interaction (PPI) network (Figure 2(a)). is network contains 412 nodes and 4297 edges. e size and color of the node are proportional to the value of betweenness and degree (Table S5), respectively. e core proteins through network analysis were ALB, TP53, INS, AKT1, and MAPK1 (Table 3) which indicated the essential role of these proteins in the pathophysiology process of oligoasthenozoospermia. Seventy-nine proteins in this network were anchored by WZESD and further discussed in the following part ( Figure 2(b)).

Co-Bioinformatics Analysis for Targets of WZESD and
Oligoasthenozoospermia-Specific Proteins. To explore the pharmacology mechanism of WZESD for oligoasthenozoospermia treatment, the relationship between WZESD targets and oligoasthenozoospermia-specific proteins was analyzed ( Figure 3). 79 proteins were overlapped in the 2 groups of protein lists (Figure 3(a)). All the selected targets of WZESD and oligoasthenospermia-specific proteins were collected for enrichment network analysis, and most of the targets of WZESD could fall into GO items with the same statistical significance as oligoasthenozoospermia-specific proteins (Figures 3(b)-3(d)), indicating the strong function association between the two comparison cohorts. Gamete generation (GO: 0007276), the developmental process involved in reproduction (GO: 0003006), regulation of hormone levels (GO: 0010817), etc., were the top 10 enriched GO clusters.

Potential Compound-Potential Target (pC-pT) Network Construction and GO and KEGG Analyses.
To further investigate the therapeutic mechanism of WZESD for oligoasthenozoospermia treatment, the 79 overlapped target proteins were clarified in detail, and a pC-pT network was constructed ( Figure 4). is network consists of 330 nodes (12 herbs, 213 potential compounds, and 79 potential targets) and 1390 edges. Compounds were manually divided into 2 groups of lists: compounds originated from WZYZ pill (LightSlateBlue) and compounds from other herbs of WZESD (red). e Wuzi Yanzong (WZYZ) pill is one of the most commonly prescribed Chinese herbal formulas for male infertility treatment. It was first recorded in the book called "Xuan Jie Lu" in AD 733 and spread around Northeast Asia over the next thousand years. were added to enhance the efficacy. So, we compared the compounds and targets between the two groups. Seven compounds, including quercetin, campesterol, kaempferol, beta-sitosterol, sitosterol, gamma-sitosterol, and stigmasterol, were overlapped between the 2 groups (DeepSkyBlue triangles). 108 specific compounds (LightSlateBlue triangles) and 7 overlapped compounds (DeepSkyBlue triangles) from WZYZ pill anchored 70 potential targets (including 12 specific). However, 98 components (7 overlapped included) in the additional group targeted 67 potential proteins (including 9 specific proteins: EDN1, PCNA, DNMT1, SLC6A14, VDAC1, VDAC3, VDAC2, BGLAP, and NFKB1). Quercetin, kaempferol, luteolin, beta-sitosterol, and        Figure 3: Co-bioinformatics analysis for targets of WZESD and oligoasthenozoospermia-specific proteins. (a) Circos plot for the 2 groups of proteins. Purple lines link the same proteins that are shared by multiple lists. Blue lines link different proteins where they fall into the same ontology term. (b) Enrichment network analysis for the union of the 2 groups of proteins. Each term is represented by a circle node, where its size is proportional to the number of input proteins that fall into that term, and its color represents its cluster identity. Terms with a similarity score >0.3 are linked by an edge (the thickness of the edge represents the similarity score). (c) e same enrichment network has its nodes displayed as pies. Each pie sector is proportional to the number of hits originated from the 2 lists. (d) e same enrichment network has its nodes coloured by the p value, as shown in the legend. e darker the color, the more statistically significant the node. 8 Evidence-Based Complementary and Alternative Medicine stigmasterol were the crucial potential active compounds, while PTGS1/2, PGR, NCOA2, ESR1, and AR were the core targets of WZESD for oligoasthenozoospermia therapy through network analysis (Tables 4 and 5). For example, quercetin is a flavonoid that has been reported to possess strong antioxidant properties and is an effective free radical scavenger, and it showed positive effects on multiple functional parameters of spermatozoa, including viability and motility [53]. Kaempferol is another kind of flavonoid, and both quercetin and kaempferol can reduce the DNA damaging in the human sperm caused by reactive oxygen species (ROS) [54]. Luteolin can protect mice from severe acute pancreatitis through anti-inflammatory and antioxidant effects [55]. For target analysis, prostaglandin G/H synthase 1/2 (PTGS1/2), namely, cyclooxygenase-1/2 (COX-1/2), is the key enzyme in the conversion of polyunsaturated fatty acids and arachidonic acid to prostaglandin (PG). COX-2 has been demonstrated to be upregulated in infertility men [56], and treatment with the COX-2 inhibitor can improve sperm motility and morphology as well as increase the pregnancy rate in infertile males [57]. Estrogen receptor 1 (ESR1) is one of the main nuclear hormone receptors for estrogen and plays an important role in the regulation of spermatogenesis [58]. : Potential compound-potential target (pC-pT) network of WZESD for treating oligoasthenozoospermia. e squares, triangles, and circles represent the herbs, potential compounds, and potential targets, respectively. LightSlateBlue: potential compounds and targets originated from the WZYZ pill. Red: potential compounds and targets derived from other herbs for WZESD. DeepSkyBlue: common compounds and targets in 2 comparison cohorts. e size of the node is proportional to the value of the degree. Efferent ductules and epididymal functions are dependent on estrogen signaling through ESR1, whose loss impaired ion transport and water reabsorption, resulting in an abnormal sperm [59]. Androgen receptor (AR) is essential for normal male reproductive development and function. AR signaling is required for the maintenance of spermatogonial numbers, blood-testis barrier integrity, completion of meiosis, adhesion of spermatids, and spermiation [60].
GO and KEGG analyses were conducted via the online "Enrichr" platform to further explain the pharmacological mechanisms [61] (Figure 5). e 79 potential targets mainly participate in nitric oxide biosynthetic and metabolic process, apoptosis, and regulation of reactive oxygen biosynthetic process, which were strongly associated with oligoasthenozoospermia ( Figure 5(a)). Reactive oxygen species (ROS) are an integral component of sperm developmental physiology, capacitation, and function. Elevated ROS levels, from processes such as infection or inflammation, can be associated with aberrations of sperm development, function, and fertilizing capacity [62]. e oxidative damage targets all cell components, reducing sperm motility and mitochondrial activity [63]. erapy with ROS scavengers can significantly improve the in vitro function of human spermatozoa and may be beneficial in infertility patients [64,65]. Apoptosis plays an important role in regulating spermatogenesis of various mammalian species, including humans [66]. is process of regulated cell death serves several important functions in the testis, a few of which include maintaining appropriate germ cells to Sertoli cell ratio, removing defective germ cells, and maintenance of overall quality control in sperm production [67]. However, high rates of apoptosis have been reported in testicular biopsies from infertile men [68]. Inflammatory disease has been established to affect male reproductive function and fertility. Male accessory gland infections account for almost 15% of all cases of male infertility seen in infertility clinics [69]. e top 10 enriched KEGG pathways were ranked in the descending order using a combined score as follows: pathways in cancer, prostate cancer, HIF-1 signaling pathway, etc. (Figure 5(b)). We found that pathways in cancer was the most significant enrichment term for the 79 target proteins; others such as hepatitis B, toxoplasmosis, and apoptosis share several same target genes with cancerrelated pathways. is implies the crosstalk among different pathways in the complex molecular biological network. Many pathways such as the Wnt signaling pathway, PI3K-AKT signaling pathway, and MAPK signaling pathway are the main part of the cancer signaling pathway, which results in sustained angiogenesis, apoptosis evading, and cell proliferation. To a certain degree, there are several common pathways between the two different pathophysiological processes of spermatogenesis and tumorigenesis. e regulation of these signaling pathways by WZESD may be the therapeutic mechanism for oligoasthenozoospermia treatment.

Molecular Docking Validation.
To evaluate the reliability of the above drug-target interactions and to further predict the accurate binding modes, a molecular docking method with AutoDock Vina [41] was employed. Molecular docking is a computational method that attempts to predict the noncovalent binding of macromolecules or, more frequently, of a macromolecule (receptor) and a small molecule (ligand) efficiently. e crystal structures of the potential targets were downloaded from the RCSB Protein Data Bank (http://www.rcsb.org). e original ligand was split and redocked to the proteins, the same as the potential compounds. A compound was considered to be an efficient drug only when the binding affinity was higher than that of the original ligand. 39 potential targets and 180 potential compounds were employed for molecular docking (Table S6). Finally, 61 potential compounds showed better binding affinity with 21 potential targets and formed 120 pairs of compound-target interactions ( Figure 6). AR,   Table S7. ree proteins, including PTGS2, XDH, and NOS2 with corresponding compounds that have the highest binding affinity, were used to display the exact binding mode (Figure 7). 1,2,5,6-Tetrahydrotanshinone was located in the binding pocket of PTGS2 and formed one conventional hydrogen bond by interacting with the key amino acid SER-530. Additionally, van der Waals with TYR-348, PHE-381, and ARG-120, pi-sigma with VAL-349, LEU-352, and ALA-527, and amide-pi stacked with GLY-526 were found in the active site which helped in the stabilization of the compound at the binding site (Figure 7(a), the global view; Figure 7(b), 3D partial view; Figure 7(c), 2D view for the interaction). In the active site of XDH, SER-876, ARG-880, SER-1008, and THR1010 formed 4 solid conventional hydrogen bonds with luteolin. Besides, other noncovalent bonds such as van der Waals, carbon-hydrogen bond, and amide-pi stacking were also formed in the binding site and helped in the stability of the complex (Figure 7(d), the global view; Figure 7(e), 3D partial view; Figure 7(f ), 2D view for the interaction). However, there were no hydrogen bonds found in the interaction between NOS2 and hypolaetin, but van der Waals with ILE-195, LEU-203, SER-236, etc., and amide-pi stacked with TRP-188 and PHE-363 were recognized in the active site and mediated the pharmacologic action (Figure 7(g), the global view; Figure 7(h), 3D partial view; Figure 7(i), 2D view for the interaction). Based on these findings, we consider that noncovalent bonds such as hydrogen bond, van der Waals, carbon-hydrogen bond, and amide-pi stacking play key roles in the protein-ligand recognition and stability, which may help determine the pharmacology activities.

Discussion
Infertility pertains to approximately 15% of sexually active couples, and male factor is present in approximately 40% of infertility cases [70]. Many factors such as age, infection, environmental pollutant, and endocrine disorder are the main cause for male infertility. e most common cause is oligoasthenozoospermia. A variety of medications have been developed to improve the sperm quality and treat male infertility such as follicle-stimulating hormone (FSH), antiestrogens, L-carnitine, and antioxidants, but showed limited curative effect, and made it still a controversial issue [71].
us, men with infertility seek for one or more alternative therapies such as herbal medicines or acupuncture   [72]. Wuzi Yanzong (WZYZ) pill, a kidney-reinforcing Chinese herbal formula, has been widely used to treat the syndrome of kidney deficiency, including impotence, spermatorrhea, premature ejaculation, and male sterility. Many research studies have reported that WZYZ pill can treat oligoasthenozoospermia by increasing sperm concentration, improving sperm motility, enhancing the activity of the acrosomal enzyme, and decreasing the sperm DNA fragmentation index [12,73]. However, the efficacy of WZYZ pill in treating male infertility remains a matter of debate [12].

Evidence-Based Complementary and Alternative Medicine
WZESD, a TCM formula development from the WZYZ pill according to TCM theory, showed a more efficient therapeutic effect for male infertility [13], but the underlying therapeutic mechanism is still unclear. Furthermore, due to the complex nature of TCM, it remains difficult to unveil such holistic medicine by the current reductionism research strategies [74]. In the current study, a network pharmacology-based method, including active ingredients' screening, multiple drug-target prediction and validation, and network construction and analysis with molecular docking, was used to clarify the pharmacological mechanisms of WZESD for treating oligoasthenozoospermia. 115 compounds (including 7 overlapped) of WZYZ pill anchored 70 potential targets (including 12 specific proteins). However, 98 components (including 7 overlapped) in the additional Van   group targeted 67 potential proteins (including 9 specific proteins). e 79 potential targets mainly participate in nitric oxide biosynthetic and metabolic process, apoptosis, and regulation of reactive oxygen biosynthetic process, which were strongly associated with oligoasthenozoospermia. Previous limited studies reported that WZYZ Pill can improve spermatogenesis, possibly through modulating the expression of the secretory proteins of Sertoli cells [75], intervene in Tip60-mediated apoptosis [73], and alleviate testicular damage via its antioxidation [76], which were consistent with our study results. Clinical trial studies have also confirmed that WZYZ Pill has obvious effects in improving the quality of male semen [77,78] and elevating the semen volume and sperm density in infertility patients with low semen counts [79], but the active compounds, direct targets, and molecule mechanisms were not mentioned in their research studies. Our study can suggest more comprehensive mechanisms of therapeutic effects of WZESD in terms of compounds, target proteins, and pathways than manual reviews of the literature.
WZYZ pill is an ancient formula and has been widely used for the treatment of male infertility. However, changes with the environment, climate, and human through hundreds of years made it less efficient. "Syndrome differentiation" and "human-environment interrelation" are 2 of the most important concepts in TCM theory. So, we need to modify the herb composition to deal with the changes. WZESD, based on WZYZ pill, added 7 commonly used herb medicines, which played the role of tonifying the liver and kidney, nourishing Yin and promoting fluid, accelerating blood circulation, removing blood stasis and dredging collaterals, and oozing wet and drenching, better accommodated the changes, and showed enhanced therapeutic effects in clinical research. In the molecule level, more active compounds (91) and proteins (9) apart from WZYZ pill were found in WZESD. We also found that quercetin, kaempferol, luteolin, beta-sitosterol, and stigmasterol were the crucial potential active compounds, while PTGS1/2, PGR, NCOA2, ESR1, and AR were the core targets of WZESD for oligoasthenozoospermia therapy (Tables 4 and 5). In another network pharmacology study, 39 bioactive compounds and 40 targets of the WZYZ pill associated with spermatogenesis disorder were obtained, and protein-protein interaction identified TP53, TNF, AKT1, Bcl-XL, Bcl-2, and IκBA as hub targets [79]. Compared with two studies, WZESD could provide more compounds and targets; it was more likely to play a role in the treatment of oligoasthenospermia through more pathways. e "multicompound", "multitarget", and "multipathway" therapeutic mechanisms were further validated and may explain the efficacy of WZYZ pill as well as why WZESD showed better effects than WZYZ pill at the molecule level.
e dominant paradigm in drug discovery is the concept of designing maximally selective ligands to act on individual drug targets. e increase in the rate of drugs failing in latestage clinical development over the past decade has been concurrent with the dominance of the assumption that the goal of drug discovery is to design exquisitely selective ligands that act on a single disease target [16]. Many common diseases such as cancer and rheumatoid arthritis as well as male infertility are complex biological systems caused by multiple molecular abnormalities. Drugs that modulate a single target might not always yield the desired outcome even if they completely interdict the functions of their direct targets [80]. With the rapid progress of bioinformatics, systems biology, and polypharmacology, network-based drug discovery is considered a promising approach toward more cost-effective drug development [17]. Coinciding with the holistic and systemic characteristics of TCM, network pharmacology is expected to bridge the gap between TCM and modern medicine. Our study successfully deciphers the efficacy of WZESD for the treatment of oligoasthenozoospermia at a molecule level and from a holistic perspective via a network pharmacology approach, which may provide a valuable reference for further experiment research studies and clinical usage.
Two main limitations of this study should be pointed out: (1) experimental investigations of WZESD for oligoasthenozoospermia treatment were not conducted; (2) the ingredients of herbs screened by our method may also be partially missing out compounds that still have effects.

Conclusions
e therapeutic effects of herbal formulas in disease management have been demonstrated by clinical practice over thousands of years, and the efficacy of WZESD for treating oligoasthenozoospermia has been proved by clinical research studies. Our work uncovers the therapeutic mechanisms of WZESD for oligoasthenozoospermia treatment from the perspective of network pharmacology and may provide a valuable reference for further experiment research studies and clinical applications.

Data Availability
All the data generated and analyzed during this study are included within this published article and its supplementary information files. e datasets supporting the conclusion in this study are available in a public database from TCMSP, BATMAN-TCM, PubChem Compound, GeneCards, Uni-Prot, STRING, Metascape, and Enrichr.

Conflicts of Interest
e authors declare that they have no conflicts of interest.

Authors' Contributions
Mingrui Hu, Yuanyuan Zhong, Jiekun Luo, and Yang Wang designed the study, collected and analyzed the data, and drafted the manuscript. Wei Xiao, Tao Tang, Shunshun Wang, Hanjin Cui, and Teng Li reviewed and corrected the manuscript. All authors read and approved the final manuscript.
14 Evidence-Based Complementary and Alternative Medicine Table S1: candidate compounds of each herb in WZESD. Table S2: detailed information for the 681 candidate targets of WZESD. Table S3: candidate compounds and their candidate targets for WZESD. Table S4: 446 target genes related to oligoasthenozoospermia obtained from Gene-Cards. Table S5: topological parameter for the generated oligoasthenozoospermia-specific protein interaction network. Table S6: the proteins and compounds performed for molecular docking.