Yilin Liao
Yiling Ding1- 1Department of Obstetrics and Gynecology, The Second Xiangya Hospital of Central South University, Changsha, Hunan, China
- 2Department of Orthopedics, The Second Xiangya Hospital of Central South University, Changsha, China
Background: Pregnancy-induced Hypertension (PIH) is a disease that causes serious maternal and fetal morbidity and mortality. Alisma Orientale (AO) has a long history of use as traditional Chinese medicine therapy for PIH. This study explores its potential mechanism and biosafety based on network pharmacology, network toxicology, molecular docking and molecular dynamics simulation.
Methods: Compounds of AO were screened in TCMSP, TCM-ID, TCM@Taiwan, BATMAN, TOXNET and CTD database; PharmMapper and SwissTargetPrediction, GeneCards, DisGeNET and OMIM databases were used to predict the targets of AO anti-PIH. The protein-protein interaction analysis and the KEGG/GO enrichment analysis were applied by STRING and Metascape databases, respectively. Then, we constructed the “herb-compound-target-pathway-disease” map in Cytoscape software to show the core regulatory network. Finally, molecular docking and molecular dynamics simulation were applied to analyze binding affinity and reliability. The same procedure was conducted for network toxicology to illustrate the mechanisms of AO hepatotoxicity and nephrotoxicity.
Results: 29 compounds with 78 potential targets associated with the therapeutic effect of AO on PIH, 10 compounds with 117 and 111 targets associated with AO induced hepatotoxicity and nephrotoxicity were obtained, respectively. The PPI network analysis showed that core therapeutic targets were IGF, MAPK1, AKT1 and EGFR, while PPARG and TNF were toxicity-related targets. Besides, GO/KEGG enrichment analysis showed that AO might modulate the PI3K-AKT and MAPK pathways in treating PIH and mainly interfere with the lipid and atherosclerosis pathways to induce liver and kidney injury. The “herb-compound-target-pathway-disease” network showed that triterpenoids were the main therapeutic compounds, such as Alisol B 23-Acetate and Alisol C, while emodin was the main toxic compounds. The results of molecular docking and molecular dynamics simulation also showed good binding affinity between core compounds and targets.
Conclusion: This research illustrated the mechanism underlying the therapeutic effects of AO against PIH and AO induced hepato-nephrotoxicity. However, further experimental verification is warranted for optimal use of AO during clinical practice.
Introduction
Pregnancy-induced Hypertension (PIH) encompasses a series of serious pregnancy complications in the second trimester (after 20 weeks of gestation), including pre-eclampsia, eclampsia, chronic hypertension with pre-eclampsia and chronic hypertension with pregnancy (Chen et al., 2019). Little is currently known about the exact etiology, although it has been established that genetic, immunological and oxidative stress are risk factors (Magalhães et al., 2020). The pathological features include endothelial cell dysfunction, multiple cytokine stimulation and activation of the coagulation system with vasospasm and increased vascular reactivity, which lead to inadequate spiral arterioles and decreased invasion of trophoblasts into the maternal decidua of the uterus, leading to a series of clinical symptoms, such as generalized edema, vomiting, blurred vision, proteinuria and even disseminated intravascular coagulation (DIC) (Chaiworapongsa et al., 2014). Furthermore, it causes fetal growth retardation, premature birth and perinatal stillborn fetus (Al Khalaf et al., 2022). Current evidence suggests that PIH is a significant threat to more than 70,000 gravidas and half a million fetuses worldwide yearly, with a high maternal mortality rate ranging from 10 to 16% (Bruno et al., 2022).
The complexity of the interactions also limits the efficacy of treatment. Currently available management consists of symptomatic treatment, including spasmolytics (magnesium sulfate), anticoagulants (aspirin), antihypertensive drugs, and termination of pregnancy as a last resort (Zhang et al., 2022a). Nevertheless, the clinical curative efficacy is highly heterogeneous and has been reported to cause maternal and fetal adverse effects (Ni et al., 2022).
Traditional Chinese medicine (TCM) therapy is an ancient practice used for more than 2000 years in China (Tian et al., 2021). The application of TCM therapy in modern society provides a novel approach for complementary and alternative medicine (CAM) treatment and is recognized by more and more people worldwide (Chu et al., 2022). Zhong-jing Zhang originally described PIH in “Synopsis in the Golden Chamber” as “nausea in pregnancy” “eclampsia” and “edema during pregnancy” (Xue et al., 2016). According to the syndrome differentiation theory, PIH is caused by Yin deficiency and Yang excess, deficiency of both spleen and kidney and blood stasis (Yeh et al., 2009). Over the years, many therapeutic methods and prescriptions of TCM have been used to treat PIH, and countless patients have benefited from it (Perry et al., 2018).
Alisma Orientale (AO), also called Ze Xie in Chinese, is a high-efficacy and low-toxicity herb medicine primarily used in Southeast Asian countries. Its efficacy is mainly mediated by removing dampness and promoting water metabolism (Wang et al., 2022), accounting for its efficacy in treating oliguria, edema and hypertension. “Women’s prescription” also recorded that AO could relieve systemic edema of gravida. In recent years, pharmacological research has demonstrated that AO compounds exhibit diuretic, anti-atherosclerotic and immunomodulatory activities, which all suitable for the treatment of PIH (Loh et al., 2017). Chen et al. reported that the natural compound of AO, alisol B 23-acetate, could inhibit Ang II-induced RAS/Wnt/β-catenin axis and attenuate podocyte injury (Chen et al., 2017). Moreover, it could inhibit collagen I, vimentin and α-smooth muscle actin at the mRNA and protein levels in rats (Chen et al., 2020). Notwithstanding that hundreds of active compounds have been isolated from AO, the specific compounds and correlative mechanisms have been largely understudied. Therefore, a novel method is urgently needed to reveal the complex treatment network and identify compounds underlying the therapeutic effect of AO against PIH. Besides, toxicology studies on AO have revealed that chronic administration may induce mild nephrotoxicity and hepatotoxicity (Yuen et al., 2006), emphasizing the need to verify the toxicity and safety profile of AO compounds and their molecular mechanisms.
The concept of network pharmacology was first officially reported by the British Pharmacologist Hopkins in 2007 (Hopkins, 2008), while network toxicology was originally proposed by Academician Liu of Chinese in 2011 (Li et al., 2019a). Both methods were based on the theories of “multi-target” and “multi-pathway” between the drug and disease, which are consistent with the characteristics of TCM therapy (Li et al., 2022).
In this study, we applied network pharmacology to screen the main therapeutic compounds of AO and predicted the potential mechanism. In addition, we analyzed the potential hepato-nephrotoxicity mechanism of AO compounds by network toxicology during the treatment of PIH. Finally, molecular docking and molecular dynamics simulation analyses were used to estimate the binding stability. Overall, this study improves our current understanding of the potential pharmacodynamic and toxicity mechanism of AO in PIH patients and explores how to enhance efficacy and reduce toxicity. The detailed research content of this study is presented in the flow chart (Figure 1).
FIGURE 1. The flow chart of this study to explore the potential molecular mechanism of AO in the treatment of PIH and toxic compounds induced hepato-nephrotoxicity.
Methods
Obtaining therapeutic and toxicity-related compounds
The Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP http://lsp.nwu.edu.cn/tcmsp.php) database (Ru et al., 2014), Traditional Chinese Medicine Information Database (Kang et al., 2013) (TCM-ID http://bidd.group), BATMAN-TCM (Liu et al., 2016) (a bioinformatics analysis Tool for molecular mechanism of traditional Chinese medicine http://bionet.ncpsb.org.cn/) and TCM@Taiwa (Chen, 2011) (http://tcm.cmu.edu.tw/zh-tw/) database were used to acquire detailed information on identified AO compounds. Oral bioavailability (OB) refers to the dosage of drugs that enter the blood circulation through oral absorption and act on local tissues and organs to produce corresponding pharmacological effects. Drug-likeness (DL) is defined as the similarity of the chemical structure of a drug compared with known drugs. Both characteristics are important to assess the potential medicinal value of compounds (Dai et al., 2022). We selected the biologically active compounds of AO in the TCMSP database using the screening criteria: OB ≥ 30 and DL ≥ 0.18 (Gao et al., 2022). Similarly, compounds from other database sources were screened in SwissADME (Daina et al., 2017) (http://www.swissadme.ch) database according to their pharmacokinetic parameters. The toxicity compound of AO was further searched in Toxicology data network (TOXNET, http://toxnet.nlm.nih.gov/index.html) database (Fowler and Schnall, 2014) and Comparative Toxicogenomics Database (CTD, https://toxnet.nlm.nih.gov/newtoxnet/ctd.htm) to analyze the biosafety profile of AO in the treatment of PIH (Stanic et al., 2021).
Potential targets related to AO compounds, PIH, liver and kidney injury
The structural formulas of therapeutic and toxic compounds obtained in the previous step were input into the PharmMapper (Liu et al., 2010) (http://www.lilab-ecust.cn/pharmmapper/) and SwissTargetPrediction database (Daina et al., 2019) (http://www.swisstargetprediction.ch/) to predict the possible targets based on the spatial conformation. The uniport ID of the target was converted to the standardized gene name by Uniport (Holzhüter and Geertsma, 2022) (https://www.uniprot.org/) database. MeSH database (https://www.ncbi.nlm.nih.gov/mesh) was utilized to verify the standard name of the disease as “Pregnancy-induced Hypertension” “drug-induced liver injury” (DILI) and “drug-induced kidney injury” (DIKI) (Zouhal et al., 2021). Then, GeneCards database (Barshir et al., 2021) (https://www.genecards.org/), DisGeNET database (Piñero et al., 2020) (https://www.disgenet.org) and Online Mendelian Inheritance in Man (Li et al., 2012) (OMIM, https://omim.org) database were retrieved to obtain potential targets using the keywords “Pregnancy-induced Hypertension” “drug-induced liver injury” and “drug-induced kidney injury”. After merging and removing duplicate targets, the intersected targets in the Venn plot were defined as potential targets of OA compounds during PIH treatment and toxic compounds associated with hepato-nephrotoxicity.
Constructing PPI network to screen core target
A protein-protein interaction (PPI) network was generated to identify interacting proteins. STRING (Szklarczyk et al., 2021) (http://string-db.org, Version 11.5) database was used to construct a PPI network of the anti-PIH effect of AO, with species limited to “Homo sapiens” and interaction score >0.4. The Cytoscape plug-ins Cytohubba and MCODE (molecular complex detection) were used to analyze the topological parameters and select the core targets of AO anti-HIP with high accuracy and exhibit the complex relationship between disease and herb (Ye et al., 2022). The same steps were used to construct the PPI network of network toxicology.
GO and KEGG pathway enrichment analysis
Functional enrichment analysis was conducted to better understand the functions of the screened genes. Gene Ontology (GO) enables the analysis of gene function based on the cellular compound (CC), molecular function (MF) and biological process (BP). Kyoto Encyclopedia of Genes and Genomes (KEGG) enables an understanding of the biological pathways associated with genes. Both were used to illustrate the core pathway and mechanism of AO anti-HIP. Targets were inputted in the Metascape database (Zhou et al., 2019) (http://www.metascape.org/), and the cut-off p-value, min overlap and enrichment value were set to 0.01, 3 and 1.5, respectively. And q-value < 0.05 (Benjamini–Hochberg method) was used to remove the false positive enrichment results (Zou et al., 2016). The R package available on the bioinformatics website (http://www.bioinformatics.com.cn/) was used to visualize the enriched results as bar and bubble plots. Then, the KEGG (Kanehisa and Sato, 2020) (http://www.kegg.jp.org/) database was used to map and color the detailed target messages in the most significantly enriched pathway. Finally, the complex regulatory network of AO in the treatment of PIH was shown by the “herb-compound-target-pathway-disease” network, incorporating the elements of AO, therapeutic compound, potential target, and top 10 enriched pathway, and visualized by Cytoscape software (v.3.9.1, https://cytoscape.org/). The same steps were used to analyze the core pathways of toxic AO compounds associated with hepato-nephrotoxicity.
Molecular dock verified the binding affinity
The Swiss dock platform (Grosdidier et al., 2011) (http://www.swissdock.ch/) is an online molecular docking (MD) tool that enables the calculation of the binding affinity of each binding site between small molecule ligand and receptor protein. The X-ray diffraction of the protein crystal structure of core targets was downloaded from the Protein Data Bank (PDB) database (www.rcsb.org) (Nakamura et al., 2022). The results were ranked by the binding affinity score, and the lowest binding affinity score value corresponded to the best binding site. The binding affinity < -7 kcal/mol indicates a strong combination possibility, which was visualized in Discovery Studio 2019 software (https://www.3ds.com/) (Sultana et al., 2022).
Molecular dynamics simulation verified the binding stability
Molecular Dynamics simulations (MDS) were further conducted to investigate the stability of small molecule ligands in proteins by the “Standard Dynamics Cascade” unit of Discovery Studio2019 software for ligand-protein complexes with the lowest binding affinity after molecular docking (Hu et al., 2022). The ligand-protein complex was placed in a solvent chamber, filled with water molecules and stabilized the electrically neutral system with Cl− and Na+. After balancing the system by the NPT ensemble (fixed the pressure, temperature and number of particles), the simulation time value was set to 500 ps and heating, equilibrium and production phases were conducted. Finally, the result was derived by trajectories analyzed with Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF) and hydrogen bond properties.
Result
Target prediction of AO compounds, PIH, liver and kidney injury
There were 77 compounds of AO obtained from TCMSP, TCM-ID, TCM@Taiwan and BATMAN database (Figure 3A) (Details were listed in Supplementary Table S1). Based on the screening criteria OB ≥ 30 and DL ≥ 0.18, 10 eligible active compounds were obtained from the TCMSP database, including sitosterol (MOL000359), Alisol B (MOL000830), Alisol B monoacetate (MOL000831), alisol b 23 acetate (MOL000832), 16β-methoxyalisol B monoacetate (MOL000849), alisol B (MOL000853), alisol C (MOL000854), alisol C monoacetate (MOL000856), 1-Monolinolein (MOL002464) and Alisol B acetate (MOL000862). There were 19 bioactive compounds from other database sources, including 13Β,17Β-Epoxyalisol A, 24-Deacetylalisol O, 25-Anhydroalisol F, Alismol, Alisol B 23-Acetate, Alisol E 23-Acetate, Alisol E 24-Acetate, Alizexol A, Neoalisol, Oriediterpenol, Oriediterpenoside and other eight compounds without formally name. The 10 toxic compounds retrieved from TOXNET and CTD databases were choline (MOL000394), emodin (MOL000472), 5-hydroxymethylfurfural (HMF, MOL000748), 1 h-indole-3-carboxylic acid (MOL000823), sucrose (MOL000842), nicotinamide (NCA, MOL000857), stearic acid (MOL000860), Healip (MOL000861), 2-Furaldehyde and 2-Furancarboxylic acid (details shown in Table1, chemistry structure shown in Figure 2).
TABLE 1. The chemical characteristics of 39 AO bioactive compounds obtained in TCMSP, TCM-ID, TCM@Taiwan, and BATMAN-TCM database.
FIGURE 2. The chemical structure of 39 AO bioactive compounds obtained from TCMSP, TCM-ID, TCM@Taiwan and BATMAN-TCM database (A), The chemical structure of top 10 therapeutic compounds (B) and (C), The chemical structure of remaining 19 therapeutic compounds (D), The chemical structure of 10 toxic compounds.
After screening and removing duplicates, we identified 138 potential targets associated with 29 therapeutic compounds and 160 targets associated with 10 AO toxic compounds. Besides, 1136 putative targets of PIH, 2188 targets of drug-induced liver injury and 2094 targets of drug-induced kidney injury were retrieved from the GeneCards, DisGeNET and OMIM databases. 78 targets overlapped between therapeutic compounds of AO and PIH, while 117 and 111 targets overlapped between AO toxic compounds and liver and kidney injury were obtained by Venn (detailed targets provided in Table 2).
TABLE 2. The detail gene name of potential targets of AO in the treatment of PIH and AO toxic compounds induced hepato-nephrotoxicity.
PPI network analysis
The 78 targets of AO in the treatment of PIH were incorporated into the STRING database to construct the PPI network. 78 nodes and 771 edges of protein-protein interaction were visually displayed in Cytoscape software, topological analysis including degree centrality (DC) and combined score, and the maximal clique centrality (MCC) algorithm was further used to screen the core targets. The top 10 targets were AKT1, BCL2L1, CASP3, EGFR, ESR1, HSP90AA1, IGF1, MAPK1, MAPK8 and SRC (Figure 3). Besides, after removing the disconnect node, the PPI network of toxic compounds-induced liver injury included 116 nodes and 1176 edges, the top 10 targets were CASP3, EGF, EGFR, ESR1, HSP90AA1, IL1B, MYC, PPARG, TNF and TP53 (Figures 5A,B). And the PPI network of toxic compounds-induced kidney injury included 109 nodes and 1130 edges, the top 10 targets were CASP3, EGF, EGFR, FOS, HSP90AA1, MMP9, MYC, PTGS2, TNF and TP53 (Supplementary Figure S1A,B).
FIGURE 3. The PPI network of AO active compounds in the treatment of PIH (A), Venn diagram of AO identified compounds from TCMSP, TCM-ID, TCM@Taiwan and BATMAN-TCM database (B), Venn diagram of potential targets to the AO treat PIH (C), The degree value of compound potential targets (D), The 78 targets PPI network gain from STRING database (E), Plug-in of Cytoscape to screen core targets.
GO and KEGG pathway enrichment analysis
A total of 1022 biological process, 50 cellular compound and 93 molecular function items were obtained from the Metascape database, and 151 KEGG pathways were significantly enriched (p < 0.01, q < 0.05), demonstrating the characteristic of multi-pathway of AO in the treatment of PIH. The top 10 biological processes included the response to hormone, regulation of kinase activity, enzyme-linked receptor protein signaling pathway, cellular response to hormone stimulus, regulation of MAPK cascade cellular response to lipid, cellular response to nitrogen compound, positive regulation of protein phosphorylation, transmembrane receptor protein tyrosine kinase and positive regulation of cell migration (Figures 4A–C).
FIGURE 4. GO/KEGG enrichment analysis of AO treat PIH (A), The bar graph of top 10 GO items, including biological process (BP), cellular compound (CC) and molecular function (MF) (B), The bubble diagram of the top 10 KEGG pathway (C), GO chord chart (D), The core pathway, PI3K-AKT cascade mapped by KEGG Mapper database, whereby therapeutic targets were colored in red, targets of AO but not in the treatment of PIH were in blue, the other targets of PIH were yellow (E), The H-C-T network map showed the importance of each AO compounds, whereby 29 compounds ranked by degree value (F), The H-C-T-P-D network map showed the complex regulation of AO treat PIH. The deep green diamond represents herb (AO), the light green hexagon represents the therapeutic compound of AO, the red and orange circles represent targets based on the degree value, the blue V icon represents the pathway involved, and the yellow rectangle is the disease (PIH).
The top 10 significantly enriched KEGG pathways included Pathways in cancer, MAPK signaling pathway, PI3K-Akt signaling pathway, Proteoglycans in cancer, Prostate cancer, Lipid and atherosclerosis, Ras signaling pathway, Endocrine resistance, Chemical carcinogenesis -receptor activation and reactive oxygen species. Based on the enrichment ratio of pathways and their correlation with diseases, the PI3K-Akt signaling pathway was the most important. We further colored the detailed targets of the PI3K-Akt pathway in the KEGG mapper database, whereby therapeutic targets were colored in red, targets of AO but not in the treatment of PIH were in blue, the other targets of PIH were yellow (Figure 4D). Besides, a total of 1262 biological process, 55 cellular compound and 117 molecular function and 168 KEGG pathway items were significantly enriched (p < 0.01, q < 0.05), and the lipid and atherosclerosis pathway was significantly associated with AO induced liver injury (Figures 5C–E). And a total of 1260BP, 52 CC, 115 MF items, and 161 KEGG pathway were enriched (p < 0.01, q < 0.05), the lipid and atherosclerosis pathway was also significantly associated with AO induced kidney injury (Supplementary Figures S1C,E).
FIGURE 5. The PPI network and GO/KEGG enrichment analysis of AO-induced liver injury (A), Venn diagram of potential targets to the AO-induced hepatotoxicity (B), The 116 targets PPI network and plug-in of Cytoscape to screen core targets (C), The bar graph of top 10 BP, CC and MF items (D), The bubble diagram of the top 10 KEGG pathway (E), The lipid and atherosclerosis pathway mapped by KEGG Mapper database (F), The H-C-T-P-D network showed the regulation mechanism of AO induced liver injury.
Herb-compound-target-pathway-disease network analysis
After identifying the core targets and the top 10 pathways from previous analysis, a “herb-compound-target-pathway-disease” (H-C-T-P-D) network was constructed by Cytoscape 3.9.1 to reveal the complex molecular mechanism and multi-compound, multi-target and multi-pathway characteristic of AO in the treatment of PIH (Figures 4E,F). The major active compounds from the network were Alisol B 23-Acetate, Alisol E 24-Acetate and alisol C. The deep green diamond represents herb (AO), the light green hexagon represents the therapeutic compound of AO, the red and orange circles represent targets based on the degree value, the blue V icon represents the pathway involved, and the yellow rectangle is the disease (PIH). Similarly, the “H-C-T-P-D” network was used to illustrate the mechanism of AO toxic compounds induced liver injury (Figure 5F) and kidney injury (Supplementary Figure S1F), whereby the green diamond, the light green hexagon, the red circle, the blue Vicon and the yellow rectangle represent the herb, the toxic compounds of AO, targets, pathways, and disease (AO induced hepato-nephrotoxicity). The main toxic compound was emodin (MOL000472).
Molecular docking
Supplementary Table S2 exhibits the binding affinity of the top 10 therapeutic compounds and emodin with their core targets, respectively. The results showed that the active compounds had good binding activities with all targets. A binding affinity < -5.0 kcal/mol indicated that the small molecule ligand has a good binding activity with the receptor protein, and binding affinity < -7.0 kcal/mol indicated stronger binding activity. Among them, we found a binding affinity < -7.0 kcal/mol between all compounds and IGF1 as well as emodin with TNF and PPARG, indicating that IGF1 may be an important target mediating the therapeutic effect against PIH, and PPARG and TNF may be the main targets of AO-induced liver injury. Figure 6A was the binding affinity heatmap of ligand-protein complexes, and Figures 6B–D shows the details of molecular docking between Alisol B monoacetate and IGF1 (lowest binding affinity); the intermolecular forces include alkyl bonding, π-alkyl bonding, conventional hydrogen bonding and carbon-hydrogen bonding. Figure 7 shows the results of emodin binding with PPARG and TNF.
FIGURE 6. The molecular docking results of AO therapeutic compounds with core targets (A), The binding affinity heatmap of ligand-receptor complex (B), 3D structures of the Alisol B monoacetate-IGF1 complex (C), spatial structure of Alisol B monoacetate-IGF1 show binding details (D), 2D structure of Alisol B monoacetate-IGF1 complex show intermolecular force types.
FIGURE 7. The molecular docking results of AO toxic compounds with corresponding core targets (A), emodin docking with TNF (B), emodin docking with PPARG (From left to right: 3D, spatial and 2D structure).
Molecular dynamics simulation
Molecular dynamics simulation represents an important technology for analyzing the ligand-protein complex’s conformational change and stability after docking. To explore the binding stability of Alisol B monoacetate with IGF1, the RMSD, RMSF curves, energies alternate tendency and hydrogen bonding heat map were calculated. Figure 8 shows that after 50 ps, the trajectories of all molecules and energy levels tend to stabilized. The RMSD, RMSF curve and hydrogen bond heat map exhibited good stability.
FIGURE 8. The molecular dynamic simulations results of Alisol B monoacetate binding with IGF1 (A), RMSD and (B), RMSF curves of MDS (C), Hydrogen bond heat map of Alisol B monoacetate binding with IGF1.
Discussion
Pregnancy-induced hypertension is a common disease during pregnancy, with a 10–12% incidence rate (Lan et al., 2022). Elevated arterial pressure and proteinuria are the most prominent clinical manifestations, and disease progression can cause reversible pathological damage to multiple systems and organs. Given the importance of the safety profile of drugs indicated for pregnant women and the unknown pathogenesis of PIH, significant emphasis has been placed on better understanding the underlying therapeutic and toxic mechanisms of these drugs (Reddy et al., 2022).
As mentioned above, oxidative stress, immune dysregulation, endothelial cell malfunction and other perplex factors have been associated with the pathogenesis of PIH. Unlike Western medicine which adopts a targeted approach, TCM fosters a more holistic approach based on the multi-compound, multi-target and multi-pathway of TCM drugs and has huge prospects for use as complementary and alternative medicine.
Alisma Orientale, also known as “Ze Xie” in Chinese, is the most important herbal compound in Ze Xie decoction indicated for treating hypertension, according to ancient records. Modern studies have isolated more than 100 active compounds, mainly terpenoids and small amounts of flavonoids, alkaloids, phytosterols and fatty acids. It is worth mentioning that protostane triterpenoids are the characteristic compounds, including alisol A-I and their derivatives (Shu et al., 2016). It has been established that Alisol B 23-acetate is an important compound with excellent bioactivities used to characterize AO in the “Pharmacopoeia of the People’s Republic of China” (Li et al., 2019b). Our findings showed that compounds with therapeutic effects included Alisol B 23-acetate, Alisol E 24-acetate and other derivatives, consistent with the “H-C-T-P-D” network.
Previous studies have demonstrated that AO plants, especially triterpene compounds, have diuretic activity. Zhang et al. reported that Alisol B, alisol B 23-acetate and other derivatives could interfere with the Na+, K+, and Cl− co-transport carrier on the luminal membrane of the Henle and sodium–chloride co-transporter in the renal distal convoluting tubule to exert diuretic action (Zhang et al., 2017). It could also significantly affect K+ excretion by competitively binding the receptor site in the collecting tube to influence sodium, potassium exchange and acid absorption. Besides, the immunomodulatory activity of the methanol extract of triterpenoids from the AO plant (alisol A, B and their acetate derivatives) significantly inhibited paw swelling in rats showing type I-IV hypersensitivity (Shao et al., 2018). In addition, alisol B and alisol B monoacetate could inhibit the complement system through the antigen–antibody-mediated process, exhibiting a good therapeutic effect against immune-related diseases (Lee et al., 2003). Further studies illustrated that alisol B 23-acetate could prevent lipid peroxidation and regulate inducible nitric oxide synthase (iNOS) expression, with significant reticuloendothelial system-potentiating activities (Matsuda et al., 1999). In conclusion, the active compounds of AO possess diuretic, immunomodulatory and vascular endothelial modulatory activities, exhibiting multi-compound, multi-target and cooperation characteristics as well as multi-compound amplification effects in the treatment of PIH.
Few studies have hitherto assessed the toxicity and side effects of AO and its compounds to determine its biosafety profile (Tian et al., 2014). Accordingly, network toxicology was applied to explore the potential mechanism in this research. Our findings showed that emodin belonging to flavonoids was the main compound that induced hepato-nephrotoxicity, and other toxic compounds with smaller proportions, including fatty acids and carbohydrates, such as stearic acid and sucrose. PPI network and GO/KEGG pathway enrichment analysis showed that AO toxic compounds mainly targeted PPARG and TNF and regulated the lipid and atherosclerosis signaling pathway. PPARG, also known as peroxisome proliferator-activated receptor gamma, is a nuclear receptor binding peroxisome proliferator that regulates adipocyte differentiation and controls the peroxisomal beta-oxidation pathway of fatty acids (Małodobra-Mazur et al., 2022). Tumor necrosis factor (TNF) is mainly secreted by macrophages and can induce insulin resistance through inhibition of Insulin receptor substrate 1 (IRS1), tyrosine phosphorylation and G kinase-anchoring protein 1 (GKAP42) degradation in adipocytes (Ando et al., 2015). Taken together, we provide preliminary evidence that the potential mechanism in AO-induced liver injury may mediate lipid metabolism via network toxicology.
PPI network analysis showed that IGF, AKT1, EGFR and MAPK1 were the main targets of AO in the treatment of PIH; the molecular docking and molecular dynamics simulation results also exhibited good binding affinity and stability. Besides, KEGG pathway enrichment analysis showed that MAPK and PI3K-AKT signal pathways were significantly enriched. Insulin-like growth factor I (IGF1) can bind to the alpha subunit of IGF1 receptor and modulate the tyrosine kinase activity to initiate the down-stream signaling events activation of PI3K-AKT and Ras-MAPK pathways (Mohindra et al., 2022). It has been reported that IGF1 is differentially expressed between PIH and normal pregnancy newborns, and IGF1 could stimulate extra-villous trophoblast (EVT) cell migration and invasion, given that insufficient invasion into the uterine endometrium may be the crucial reason for PIH (Seedorf et al., 2020). Epidermal growth factor receptor (EGFR) and RAC-alpha serine/threonine-protein kinase (AKT1) are important hub genes that regulate a series of biological functions, including angiogenesis (Zhang et al., 2022b). Mitogen-activated protein kinase 1 (MAPK1) has been associated with the transduction of endothelial inflammatory response and inflammatory factor expression action (Galganska et al., 2021). Wang Z et al. indicated that overexpressing miR-106 and downregulating the MAPK signaling cascade could attenuate oxidative stress injury and inflammatory response in the liver of mice with PIH (Wang et al., 2019). Li X et al. found significant Cx43 protein overexpression on human umbilical vein endothelial cells (HUVECs) in PIH patients, which may lead to monocyte-endothelial adhesion increase by PI3K-AKT signaling pathway activation and upregulate the downstream genes of VCAM-1 and ICAM-1 (Li et al., 2019c).
In conclusion, we found that the main compounds of AO in treating PIH are triterpenoids, especially Alisol B, Alisol C and their derivatives, which target IGF, AKT1, EGFR and MAPK1 and mediate the MAPK and PI3K-AKT signal pathways to produce diuretic and immunomodulatory activities and regulate the function of endothelial cells, anti-vascular inflammation and oxidative stress. In addition, our study preliminarily indicated that the toxic compounds in AO were mainly emodin and small amounts of fatty acids, which may produce hepato-nephrotoxicity by interfering with lipid metabolism. Furthermore, our study suggests that an optimal therapeutic effect against PIH may be achieved by increasing the triterpenoid content and reducing toxic compounds through appropriate methods. However, some shortcomings in the present study affect the robustness of our findings to some extent. First, our study was based on the identified and main bioactive compounds of AO, which may have some selection bias. Second, and crucially, further in vivo and in vitro basic research is needed to confirm and extend our conclusions. And this was exactly what we need to do next, to verify the mechanisms of AO core compounds with therapeutic efficacy in PIH and hepato-nephrotoxicity in future research.
Conclusion
This study was based on the bioinformatics analysis of network pharmacology, network toxicology, molecular docking and molecular dynamics simulation. The finding shows that the triterpenoids were the main therapeutic compounds and emodin was the main toxic compound of Alisma Orientale. And also illustrated the potential mechanism underlying the therapeutic effects of AO against PIH and AO induced hepato-nephrotoxicity. In short, our research could guide directions for subsequent basic experimental verification and provide points of view for enhance efficacy and reduce toxicity in clinical.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
YL designed the study and wrote the manuscript, MY derived the datasets and performed statistical analyses; YD revised the manuscript; LY and CX draw all figures and prepared the supplementary materials.
Acknowledgments
We sincerely appreciate the databases of TCMSP, TCM-ID, TCM@Taiwan, BATMAN-TCM, TOXNET, CTD, PharmMapper, Swiss Institute of Bioinformatics (SIB, including SwissADME, SwissTargetPrediction and Swiss Dock), Uniport, GeneCards, DisGeNET, OMIM, STRING, Metascape, bioinformatics and Cytoscape, Discovery Studio 2019 software for provided the data and making the charts.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.1027112/full#supplementary-material
References
Al Khalaf, S., Khashan, A. S., Chappell, L. C., O'Reilly É, J., and McCarthy, F. P. (2022). Role of antihypertensive treatment and blood pressure control in the occurrence of adverse pregnancy outcomes: A population-based study of linked electronic health records. Hypertension 79, 1548–1558. doi:10.1161/HYPERTENSIONAHA.122.18920
Ando, Y., Shinozawa, Y., Iijima, Y., Yu, B. C., Sone, M., Ooi, Y., et al. (2015). Tumor necrosis factor (TNF)-α-induced repression of GKAP42 protein levels through cGMP-dependent kinase (cGK)-Iα causes insulin resistance in 3T3-L1 adipocytes. J. Biol. Chem. 290, 5881–5892. doi:10.1074/jbc.M114.624759
Barshir, R., Fishilevich, S., Iny-Stein, T., Zelig, O., Mazor, Y., Guan-Golan, Y., et al. (2021). GeneCaRNA: A comprehensive gene-centric database of human non-coding RNAs in the GeneCards suite. J. Mol. Biol. 433, 166913. doi:10.1016/j.jmb.2021.166913
Bruno, A. M., Allshouse, A. A., Metz, T. D., and Theilen, L. H. (2022). Trends in hypertensive disorders of pregnancy in the United States from 1989 to 2020. Obstet. Gynecol. 140, 83–86. doi:10.1097/AOG.0000000000004824
Chaiworapongsa, T., Chaemsaithong, P., Yeo, L., and Romero, R. (2014). Pre-eclampsia part 1: Current understanding of its pathophysiology. Nat. Rev. Nephrol. 10, 466–480. doi:10.1038/nrneph.2014.102
Chen, C. Y. (2011). TCM Database@Taiwan: The world's largest traditional Chinese medicine database for drug screening in silico, PLoS One 6, e15939. doi:10.1371/journal.pone.0015939
Chen, H., Wang, M. C., Chen, Y. Y., Chen, L., Wang, Y. N., Vaziri, N. D., et al. (2020). Alisol B 23-acetate attenuates CKD progression by regulating the renin-angiotensin system and gut-kidney axis. Ther. Adv. Chronic Dis. 11, 2040622320920025. doi:10.1177/2040622320920025
Chen, L., Chen, D. Q., Wang, M., Liu, D., Chen, H., Dou, F., et al. (2017). Role of RAS/Wnt/β-catenin axis activation in the pathogenesis of podocyte injury and tubulo-interstitial nephropathy. Chem. Biol. Interact. 273, 56–72. doi:10.1016/j.cbi.2017.05.025
Chen, S., Li, N., Mei, Z., Ye, R., Li, Z., Liu, J., et al. (2019). Micronutrient supplementation during pregnancy and the risk of pregnancy-induced hypertension: A randomized clinical trial. Clin. Nutr. 38, 146–151. doi:10.1016/j.clnu.2018.01.029
Chu, H., Moon, S., Park, J., Bak, S., Ko, Y., and Youn, B. Y. (2022). The use of artificial intelligence in complementary and alternative medicine: A systematic scoping review. Front. Pharmacol. 13, 826044. doi:10.3389/fphar.2022.826044
Dai, W., Chen, C., Dong, G., Li, G., Peng, W., Liu, X., et al. (2022). Alleviation of fufang fanshiliu decoction on type II diabetes mellitus by reducing insulin resistance: A comprehensive network prediction and experimental validation. J. Ethnopharmacol. 294, 115338. doi:10.1016/j.jep.2022.115338
Daina, A., Michielin, O., and Zoete, V. (2017). SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 7, 42717. doi:10.1038/srep42717
Daina, A., Michielin, O., and Zoete, V. (2019). SwissTargetPrediction: Updated data and new features for efficient prediction of protein targets of small molecules. Nucleic Acids Res. 47, W357-W364–w364. doi:10.1093/nar/gkz382
Fowler, S., and Schnall, J. G. (2014). Toxnet: Information on toxicology and environmental health. Am. J. Nurs. 114, 61–63. doi:10.1097/01.NAJ.0000443783.75162.79
Galganska, H., Jarmuszkiewicz, W., and Galganski, L. (2021). Carbon dioxide inhibits COVID-19-type proinflammatory responses through extracellular signal-regulated kinases 1 and 2, novel carbon dioxide sensors. Cell. Mol. Life Sci. 78, 8229–8242. doi:10.1007/s00018-021-04005-3
Gao, Y., Ji, W., Lu, M., Wang, Z., Jia, X., Wang, D., et al. (2022). Systemic pharmacological verification of Guizhi Fuling decoction in treating endometriosis-associated pain. J. Ethnopharmacol. 297, 115540. doi:10.1016/j.jep.2022.115540
Grosdidier, A., Zoete, V., and Michielin, O. (2011). SwissDock, a protein-small molecule docking web service based on EADock DSS. Nucleic Acids Res. 39, W270–W277. doi:10.1093/nar/gkr366
Holzhüter, K., and Geertsma, E. R. (2022). Uniport, not proton-symport, in a non-mammalian SLC23 transporter. J. Mol. Biol. 434, 167393. doi:10.1016/j.jmb.2021.167393
Hopkins, A. L. (2008). Network pharmacology: The next paradigm in drug discovery. Nat. Chem. Biol. 4, 682–690. doi:10.1038/nchembio.118
Hu, Y., Wu, Y., Jiang, C., Wang, Z., Shen, C., Zhu, Z., et al. (2022). Investigative on the molecular mechanism of licorice flavonoids anti-melanoma by network pharmacology, 3D/2D-QSAR, molecular docking, and molecular dynamics simulation. Front. Chem. 10, 843970. doi:10.3389/fchem.2022.843970
Kanehisa, M., and Sato, Y. (2020). KEGG Mapper for inferring cellular functions from protein sequences. Protein Sci. 29, 28–35. doi:10.1002/pro.3711
Kang, H., Tang, K., Liu, Q., Sun, Y., Huang, Q., Zhu, R., et al. (2013). HIM-herbal ingredients in-vivo metabolism database. J. Cheminform. 5, 28. doi:10.1186/1758-2946-5-28
Lan, X., Guo, L., Zhu, S., Cao, Y., Niu, Y., Han, S., et al. (2022). First-trimester serum cytokine profile in pregnancies conceived after assisted reproductive technology (ART) with subsequent pregnancy-induced hypertension. Front. Immunol. 13, 930582. doi:10.3389/fimmu.2022.930582
Lee, S. M., Kim, J. H., Zhang, Y., An, R. B., Min, B. S., Joung, H., et al. (2003). Anti-complementary activity of protostane-type triterpenes from Alismatis rhizoma. Arch. Pharm. Res. 26, 463–465. doi:10.1007/BF02976863
Li, S., Wang, L., Du, Z., Jin, S., Song, C., Jia, S., et al. (2019). Identification of the lipid-lowering component of triterpenes from Alismatis rhizoma based on the MRM-based characteristic chemical profiles and support vector machine model. Anal. Bioanal. Chem. 411, 3257–3268. doi:10.1007/s00216-019-01818-x
Li, X., Chen, H., Yang, H., Liu, J., Li, Y., Dang, Y., et al. (2022). Study on the potential mechanism of tonifying kidney and removing dampness formula in the treatment of postmenopausal dyslipidemia based on network pharmacology, molecular docking and experimental evidence. Front. Endocrinol. (Lausanne). 13, 918469. doi:10.3389/fendo.2022.918469
Li, X., Zhang, Q., Zhang, R., Cheng, N., Guo, N., Liu, Y., et al. (2019). Down-regulation of Cx43 expression on PIH-HUVEC cells attenuates monocyte-endothelial adhesion. Thromb. Res. 179, 104–113. doi:10.1016/j.thromres.2019.05.009
Li, Y., Zhang, Y., Wang, Y., Li, Y., Yang, F., Zhang, P., et al. (2019). A strategy for the discovery and validation of toxicity quality marker of Chinese medicine based on network toxicology. Phytomedicine 54, 365–370. doi:10.1016/j.phymed.2018.01.018
Li, Z., Ying, B., Liu, X., Zhang, X., and Yu, H. (2012). An examination of the OMIM database for associating mutation to a consensus reference sequence. Protein Cell. 3, 198–203. doi:10.1007/s13238-012-2037-2
Liu, X., Ouyang, S., Yu, B., Liu, Y., Huang, K., Gong, J., et al. (2010). PharmMapper server: A web server for potential drug target identification using pharmacophore mapping approach. Nucleic Acids Res. 38, W609–W614. doi:10.1093/nar/gkq300
Liu, Z., Guo, F., Wang, Y., Li, C., Zhang, X., Li, H., et al. (2016). BATMAN-TCM: A bioinformatics analysis tool for molecular mechANism of traditional Chinese medicine. Sci. Rep. 6, 21146. doi:10.1038/srep21146
Loh, Y. C., Tan, C. S., Ch'ng, Y. S., Ahmad, M., Asmawi, M. Z., and Yam, M. F. (2017). Vasodilatory effects of combined traditional Chinese medicinal herbs in optimized ratio. J. Med. Food 20, 265–278. doi:10.1089/jmf.2016.3836
Magalhães, E., Méio, M., Peixoto-Filho, F. M., Gonzalez, S., da Costa, A. C. C., and Moreira, M. E. L. (2020). Pregnancy-induced hypertension, preterm birth, and cord blood adipokine levels. Eur. J. Pediatr. 179, 1239–1246. doi:10.1007/s00431-020-03586-8
Małodobra-Mazur, M., Cierzniak, A., Ryba, M., Sozański, T., Piórecki, N., and Kucharska, A. Z. (2022). Cornus mas L. Increases glucose uptake and the expression of PPARG in insulin-resistant adipocytes. Nutrients 14, 2307. doi:10.3390/nu14112307
Matsuda, H., Kageura, T., Toguchida, I., Murakami, T., Kishi, A., and Yoshikawa, M. (1999). Effects of sesquiterpenes and triterpenes from the rhizome of Alisma orientale on nitric oxide production in lipopolysaccharide-activated macrophages: Absolute stereostructures of alismaketones-B 23-acetate and -C 23-acetate. Bioorg. Med. Chem. Lett. 9, 3081–3086. doi:10.1016/s0960-894x(99)00536-3
Mohindra, V., Chowdhury, L. M., Chauhan, N., Maurya, R. K., and Jena, J. K. (2022). Transcriptome analysis revealed hub genes for muscle growth in Indian major carp, Catla catla (Hamilton, 1822). Genomics 114, 110393. doi:10.1016/j.ygeno.2022.110393
Nakamura, R., Ogawa, S., Takahashi, Y., and Fujishiro, T. (2022). Cycloserine enantiomers inhibit PLP-dependent cysteine desulfurase SufS via distinct mechanisms. Febs J. 289, 5947–5970. doi:10.1111/febs.16455
Ni, Z., Mei, S., You, S., Lin, Y., Cheng, W., Zhou, L., et al. (2022). Adverse effects of polycystic ovarian syndrome on pregnancy outcomes in women with frozen-thawed embryo transfer: Propensity score-matched study. Front. Endocrinol. 13, 878853. doi:10.3389/fendo.2022.878853
Perry, H., Khalil, A., and Thilaganathan, B. (2018). Preeclampsia and the cardiovascular system: An update. Trends cardiovasc. Med. 28, 505–513. doi:10.1016/j.tcm.2018.04.009
Piñero, J., Ramírez-Anguita, J. M., Saüch-Pitarch, J., Ronzano, F., Centeno, E., Sanz, F., et al. (2020). The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 48, D845-D855–d855. doi:10.1093/nar/gkz1021
Reddy, M., Palmer, K., Rolnik, D. L., Wallace, E. M., Mol, B. W., and Da Silva Costa, F. (2022). Role of placental, fetal and maternal cardiovascular markers in predicting adverse outcome in women with suspected or confirmed pre-eclampsia. Ultrasound Obstet. Gynecol. 59, 596–605. doi:10.1002/uog.24851
Ru, J., Li, P., Wang, J., Zhou, W., Li, B., Huang, C., et al. (2014). Tcmsp: A database of systems pharmacology for drug discovery from herbal medicines. J. Cheminform. 6, 13. doi:10.1186/1758-2946-6-13
Seedorf, G., Kim, C., Wallace, B., Mandell, E. W., Nowlin, T., Shepherd, D., et al. (2020). rhIGF-1/BP3 preserves lung growth and prevents pulmonary hypertension in experimental bronchopulmonary dysplasia. Am. J. Respir. Crit. Care Med. 201, 1120–1134. doi:10.1164/rccm.201910-1975OC
Shao, C., Fu, B., Ji, N., Pan, S., Zhao, X., Zhang, Z., et al. (2018). Alisol B 23-acetate inhibits IgE/Ag-mediated mast cell activation and allergic reaction. Int. J. Mol. Sci. 19, 4092. doi:10.3390/ijms19124092
Shu, Z., Pu, J., Chen, L., Zhang, Y., Rahman, K., Qin, L., et al. (2016). Alisma orientale: Ethnopharmacology, phytochemistry and pharmacology of an important traditional Chinese medicine. Am. J. Chin. Med. 44, 227–251. doi:10.1142/S0192415X16500142
Stanic, B., Samardzija Nenadov, D., Fa, S., Pogrmic-Majkic, K., and Andric, N. (2021). Integration of data from the cell-based ERK1/2 ELISA and the Comparative Toxicogenomics Database deciphers the potential mode of action of bisphenol A and benzo[a]pyrene in lung neoplasm. Chemosphere 285, 131527. doi:10.1016/j.chemosphere.2021.131527
Sultana, N., Chung, H. J., Emon, N. U., Alam, S., Taki, M. T. I., Rudra, S., et al. (2022). Biological functions of Dillenia pentagyna roxb. Against pain, inflammation, fever, diarrhea, and thrombosis: Evidenced from in vitro, in vivo, and molecular docking study. Front. Nutr. 9, 911274. doi:10.3389/fnut.2022.911274
Szklarczyk, D., Gable, A. L., Nastou, K. C., Lyon, D., Kirsch, R., Pyysalo, S., et al. (2021). The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 49, D605–d612. doi:10.1093/nar/gkaa1074
Tian, N. N., Zheng, Y. B., Li, Z. P., Zhang, F. W., and Zhang, J. F. (2021). Histone methylatic modification mediates the tumor-suppressive activity of curcumol in hepatocellular carcinoma via an Hotair/EZH2 regulatory axis. J. Ethnopharmacol. 280, 114413. doi:10.1016/j.jep.2021.114413
Tian, T., Chen, H., and Zhao, Y. Y. (2014). Traditional uses, phytochemistry, pharmacology, toxicology and quality control of Alisma orientale (sam.) juzep: A review. J. Ethnopharmacol. 158 Pt A, 373–387. doi:10.1016/j.jep.2014.10.061
Wang, B., Peng, Q., Wang, R., Yu, S., Li, Q., and Huang, C. (2022). Efficient Microcystis removal and sulfonamide-resistance gene propagation mitigation by constructed wetlands and functional genes analysis. Chemosphere 292, 133481. doi:10.1016/j.chemosphere.2021.133481
Wang, Z., Bao, X., Song, L., Tian, Y., and Sun, P. (2019). Role of miR-106-mediated mitogen-activated protein kinase signaling pathway in oxidative stress injury and inflammatory infiltration in the liver of the mouse with gestational hypertension. J. Cell. Biochem. 122, 958–968. doi:10.1002/jcb.29552
Xue, B., Yu, Y., Zhang, Z., Guo, F., Beltz, T. G., Thunhorst, R. L., et al. (2016). Leptin mediates high-fat diet sensitization of angiotensin II-elicited hypertension by upregulating the brain renin-angiotensin system and inflammation. Hypertension 67, 970–976. doi:10.1161/HYPERTENSIONAHA.115.06736
Ye, X. W., Wang, H. L., Cheng, S. Q., Xia, L. J., Xu, X. F., and Li, X. R. (2022). Network pharmacology-based strategy to investigate the pharmacologic mechanisms of coptidis rhizoma for the treatment of alzheimer's disease. Front. Aging Neurosci. 14, 890046. doi:10.3389/fnagi.2022.890046
Yeh, H. Y., Chen, Y. C., Chen, F. P., Chou, L. F., Chen, T. J., and Hwang, S. J. (2009). Use of traditional Chinese medicine among pregnant women in Taiwan. Int. J. Gynaecol. Obstet. 107, 147–150. doi:10.1016/j.ijgo.2009.07.024
Yuen, M. F., Tam, S., Fung, J., Wong, D. K., Wong, B. C., and Lai, C. L. (2006). Traditional Chinese medicine causing hepatotoxicity in patients with chronic Hepatitis B infection: A 1-year prospective study. Aliment. Pharmacol. Ther. 24, 1179–1186. doi:10.1111/j.1365-2036.2006.03111.x
Zhang, M., Li, J., Fu, X., Zhang, Y., Zhang, T., Wu, B., et al. (2022). Endometrial thickness is an independent risk factor of hypertensive disorders of pregnancy: A retrospective study of 13, 458 patients in frozen-thawed embryo transfers. Reprod. Biol. Endocrinol. 20, 93. doi:10.1186/s12958-022-00965-8
Zhang, X., Li, X. Y., Lin, N., Zhao, W. L., Huang, X. Q., Chen, Y., et al. (2017). Diuretic activity of compatible triterpene components of alismatis rhizoma. Molecules 22, 1459. doi:10.3390/molecules22091459
Zhang, X., Liu, L., Liu, D., Li, Y., He, J., and Shen, L. (2022). 17β-Estradiol promotes angiogenesis of bone marrow mesenchymal stem cells by upregulating the PI3K-Akt signaling pathway. Comput. Struct. Biotechnol. J. 20, 3864–3873. doi:10.1016/j.csbj.2022.07.028
Zhou, Y., Zhou, B., Pache, L., Chang, M., Khodabakhshi, A. H., Tanaseichuk, O., et al. (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 10, 1523. doi:10.1038/s41467-019-09234-6
Zou, X., Holmes, E., Nicholson, J. K., and Loo, R. L. (2016). Automatic spectroscopic data categorization by clustering analysis (asclan): A data-driven approach for distinguishing discriminatory metabolites for phenotypic subclasses. Anal. Chem. 88, 5670–5679. doi:10.1021/acs.analchem.5b04020
Keywords: pregnancy induced hypertension, Alisma orientale, network pharmacology, network toxicology, molecular docking
Citation: Liao Y, Ding Y, Yu L, Xiang C and Yang M (2022) Exploring the mechanism of Alisma orientale for the treatment of pregnancy induced hypertension and potential hepato-nephrotoxicity by using network pharmacology, network toxicology, molecular docking and molecular dynamics simulation. Front. Pharmacol. 13:1027112. doi: 10.3389/fphar.2022.1027112
Received: 09 September 2022; Accepted: 18 October 2022;
Published: 15 November 2022.
Edited by:
Swaran J. S. Flora, National Institute of Pharmaceutical Education and Research, IndiaReviewed by:
YuFeng Zhang, Nanjing University of Chinese Medicine, ChinaJaved Ahmad, Najran University, Saudi Arabia
Huachong Xu, Jinan University, China
Copyright © 2022 Liao, Ding, Yu, Xiang and Yang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mengyuan Yang, yangmengyuan_xy@csu.edu.cn