Systematic Understanding of Anti-Aging Effect of Coenzyme Q10 on Oocyte Through a Network Pharmacology Approach

Abstract

Background:

Maternal oocyte aging is strongly contributing to age-related decline in female fertility. Coenzyme Q10 (CoQ10) exerts positive effects in improving aging-related deterioration of oocyte quality, but the exact mechanism is unclear.

Objective:

To reveal the system-level mechanism of CoQ10’s anti-aging effect on oocytes based on network pharmacology.

Methods:

This study adopted a systems network pharmacology approach, including target identification, data integration, network and module construction, bioinformatics analysis, molecular docking, and molecular dynamics simulation.

Result:

A total of 27 potential therapeutic targets were screened out. Seven hub targets (PPARA, CAT, MAPK14, SQSTM1, HMOX1, GRB2, and GSR) were identified. Functional and pathway enrichment analysis indicated that these 27 putative targets exerted therapeutic effects on oocyte aging by regulating signaling pathways (e.g., PPAR, TNF, apoptosis, necroptosisn, prolactin, and MAPK signaling pathway), and are involved oxidation-reduction process, mitochondrion, enzyme binding, reactive oxygen species metabolic process, ATP binding, among others. In addition, five densely linked functional modules revealed the potential mechanisms of CoQ10 in improving aging-related deterioration of oocyte quality are closely related to antioxidant, mitochondrial function enhancement, autophagy, anti-apoptosis, and immune and endocrine system regulation. The molecular docking study reveals that seven hub targets have a good binding affinity towards CoQ10, and molecular dynamics simulation confirms the stability of the interaction between the hub targets and the CoQ10 ligand.

Conclusion:

This network pharmacology study revealed the multiple mechanisms involved in the anti-aging effect of CoQ10 on oocytes. The molecular docking and molecular dynamics simulation provide evidence that CoQ10 may act on these hub targets to fight against oocytes aging.

Introduction

Aging leads to a gradual and continuous decline of physiological processes. Oogenesis is one of the first processes to fail during human aging. Human oocyte development begins in the early fetal life, reaches its peak at 20 weeks of gestation, and remains dormant before complete meiosis I (1). Over time, the quality and quantity of these oocytes gradually declined, concurrently with the reduced reproductive ability and increased occurrence of aneuploidy, miscarriage, and birth defects (2–4). In contrast, the aging-related decline in fecundity was negated in patients using oocytes obtained from healthy young donors (5). These data suggest that the age-related decline in fertility is strongly attributed to a decline in the oocyte quality.

Oocytes accumulate a large pool of maternal RNA and protein during meiotic maturation to support early embryogenesis before embryonic genome activation (6, 7). Therefore, oocyte quality is a critical prerequisite for fertilization and subsequent embryo development (8). In the clinic, the paradigm that aged oocytes represent low‐quality oocytes is widely confirmed. Data from a reproductive center in China shows that among 20,687 women receiving IVF, the complete full cycle live birth rate of women under 31 years old was 63.81%, while the rate of women over 40 years old decreased to 4.71% (3). However, due to education, social, economic, and other factors, the average age of first-time mothers is rising worldwide, which poses a great challenge to reproductive medicine (9, 10).

Various interventions have been proposed to improve oocyte quality, thereby improving pregnancy outcomes in oocyte aging-related infertility. As the third most consumed dietary supplement, coenzyme Q10 (CoQ10) has attracted global interest due to its crucial role in mitochondrial bioenergy, anti-oxidation, anti-aging, and immune system regulation. In the meantime, the role of CoQ10 in human reproduction has received widespread attention, especially in improving oocyte quality and counteracting oocyte aging (11–14). Turi et al. confirmed for the first time that CoQ10 is present in human follicular fluid, and its level was significantly higher in the mature oocytes versus dysmorphic oocytes; correspondingly, the CoQ10 level in grade I-II embryos was significantly higher than that of III-IV embryos (15). A further clinical study conducted by Giannubilo et al. indicated that oral supplements with CoQ10 increased the content of CoQ10 in follicular fluid, and CoQ10 in follicular fluid was positively correlated with oocytes quality (12). In addition, pretreatment with CoQ10 significantly increased the number of retrieved oocytes, fertilization rate, and the number of high-quality embryos in women with poor ovarian response (POR) during the IVF-ICSI cycle (13). Rodent studies have also confirmed that CoQ10 counteracts physiological ovarian aging and improves ovarian reserve through multiple mechanisms (16–19).

Currently, although anti-oocyte aging activities exerted by CoQ10 have been reported, details of biomarkers and the biological pathways through which CoQ10 counteracts oocyte aging are limited. In previous studies, a network pharmacology-based approach was successfully used to reveal all candidate targets, functions, and potential therapeutic mechanisms of drugs (20–22). In the present study, system network pharmacology combined with molecular docking was performed to reveal the potential targets and therapeutic mechanisms of CoQ10 against oocyte aging. The workflow is shown in Figure 1.

Figure 1

Materials and Methods

Selection of Oocyte Aging-Associated Genes from scRNA-Seq Results

Oocyte aging-associated genes were searched in the GEO database (https://www.ncbi.nlm.nih.gov/geo/), and the series (GSE96477 and GSE155179) were selected. The GSE95477 and GSE155179 include single-cell RNA sequencing (scRNA-seq) data of MII human oocytes from younger (≤30 years) and older (≥40 years) patients. DESeq2 software (version 1.32.0) was used to determine differentially-expressed genes (DEGs) in all datasets. Genes with a P value <0.05 and |log 2 (fold change) | >1 were considered to be DEGs, namely, oocyte aging-related genes.

In addition, we searched Pubmed databases (https://pubmed.ncbi.nlm.nih.gov) using human oocyte aging and scRNA-seq as keywords. We include Zhang’s research which isolated mRNA from human oocytes obtained from old and young patients for scRNA-seq and subsequently identified differentially expressed genes (23). In summary, three scRNA-seq data sets were included for differential gene analysis. All differentially expressed genes were integrated into oocyte aging genes.

Identification of Putative CoQ10 Targets

Simplified molecular-input line-entry specification (SMILES) information and 3D structure of CoQ10 were obtained from PubChem database (https://pubchem.ncbi.nlm.nih.gov/). PharmMapper server (https://lilabecust.cn/pharmmapper/), Drug bank (https://go.drugbank.com/),TargetNet (http://targetnet.scbdd.com/), and Comparative Toxicogenomics Database (CTD, https://ctdbase.org/) were used to identify possible targets for CoQ10.

Construction of Protein-Protein Interaction Network

The oocyte aging-rated genes and the CoQ10 action targets were intersected to obtain the genes of CoQ10 against oocyte aging. The Venn diagram was drawn by using the bioinformatics platform Jvenn (http://jvenn.toulouse.inra.fr/app/index.html) and the protein-protein interaction (PPI) network of the common targets was obtained by using the STRING database (https://string-db.org/). The species was limited to Homo sapiens, and we selected the combined score >0.4. The Cytoscape 3.8.2 was used to construct the PPI network.

GO and KEGG Pathway Enrichment

The OmicShare tools (https://www.omicshare.com/), a free data analysis platform, were used to analyze the GO and KEGG enrichment of the potential therapeutic targets. The Go analysis includes biological processes (BP), cellular components (CC), molecular function (MF).

Construction of Network

Two networks were constructed: (1) the PPI network was constructed using the genes of CoQ10 against oocyte aging; (2) the CoQ10-targets-pathways network, this network was divided into different functional modules using the Community Cluster (Glay) algorithm of clusterMaker2.

Molecular Docking of the Active Components-Targets

The 3D structures of protein receptors were obtained from the RCSB Protein Data Bank (https://www.rcsb.org/) (24). PyMol software (v2.4.0) was applied to remove the heteroatom and water molecules from the proteins. The 2D structure of CoQ10, used as a ligand, was downloaded from the PubChem website (https://pubchem.ncbi.nlm.nih.gov/) (25). AutoDockTools (v1.5.7) software was used to convert the above structure into a pdbqt file. The prediction of binding sites were obtained from Deepsite website (https://www.playmolecule.org/deepsite/). Molecular docking was completed by Autodock vina software (v1.1.2) (26). The complexes were then observed and plotted using PyMol software.

Molecular Dynamic Simulation

The conformations with the highest scores obtained from the molecular docking (protein-ligand complexes) were further subjected to molecular dynamic (MD) simulations using Gromacs 2020.6. (27). The CHARMM36 force field was used to generate the parameters of the protein (28). The Acpype Server (https://www.bio2byte.be/acpype/) (29) was used for ligand topology. The protein-ligand complex system was built in a dodecahedron box and filled with three-point water model. Sodium ions or chloride ions neutralize the charge of the system, making it electrically neutral. The system was energy-minimized for 50,000 steps at the temperature of 300k and equilibrated by 100 ps of NVT and NPT simulations subsequently. Then, a production run of 30 ns was performed with a step length of 2.0 fs, and the coordinate track file was saved every 10 ps. The root-mean-square deviation (RMSD) and root-mean-square fluctuation (RMSF) throughout the trajectory were analyzed using the gmx rms and rmsf, respectively. The g_mmpbsa program was used to evaluated the binding free energy of the complexes.

Results

Oocyte Aging-Related Genes

The GSE95477 data set derived from the GEO database was analyzed and 61 DEGs were obtained according to the cut-off criteria. 1455 DEGs of GSE155179 date are obtained via the same method.

In addition, we downloaded the supplementary data of Zhang’s research (23). Then, we re-screened DEGs according to the cut-off criteria and obtained 1439 DEGs. Finally, all the obtained DEGs were merged, 2807 oocyte aging-related genes were identified (Supplementary 1).

CoQ10−Putative Targets

A total of 242 CoQ10 targets were obtained after removing duplications from the PharmMapper, DrugBank, TargetNet and CTD databases (Supplementary 2).

PPI Network of the Potential Therapeutic Targets

Based on the results above, 31 common target genes of CoQ10-putative targets and oocyte aging-related genes were identified after removing EDN1which was up-regulated in GSE155179 but down-regulated in Zhang’s research (Figure 2A).

Figure 2

The PPI network of common target genes consisted of 27 nodes and 51 edges, and 4 free nodes (targets) were removed.

The 27 nodes of the PPI network were identified as potential therapeutic targets (Figure 2B). The CytoNCA, an app plugin in Cytoscape, was used to analyze the topological parameters of the network (Supplementary 3).

To find the hub targets in this complex biological network, the topological parameters including degree centrality (DC), betweenness centrality (BC), closeness centrality (CC), eigenvector centrality (EC), local average connectivity-based method (LAC), network centrality (NC) were analyzed. As a result, a total of seven genes (PPARA, CAT, MAPK14, SQSTM1, HMOX1, GRB2, and GSR), were obtained according to the median values of DC, BC, CC, EC, LAC, NC (Supplementary 4). The seven genes were used for the subsequent molecular docking and molecular dynamic simulation.

GO Enrichment Analysis

To further study the molecular mechanism of CoQ10 against oocyte aging, GO enrichment of the potential therapeutic targets was performed based on BP, CC, and MF (Figure 3A). In terms of BP, the anti-aging effect of CoQ10 on oocyte mainly involves the regulation of reactive oxygen species, participation of immune process, and response to hormone etc. The MF mainly involves nuclear transcription, cytokine activity, and ATP binding etc. The CC is mainly manifested in cell membrane, cytoplasm, and mitochondrion etc. (Figure 3B).

Figure 3

KEGG Enrichment Analysis

The KEGG pathway enrichment was identified based on the potential therapeutic targets. 53 pathways with P value < 0.05 were obtained. After data screening, 10 significant pathways were selected (Figure 4A and Table 1).

Figure 4

Then, the CoQ10-targets-pathways network was constructed after sorting out the genes enriched in each pathway (Figure 4B). To further investigate the mechanism of CoQ10 against oocyte aging, five functional modules were constructed using the Glay algorithm of clusterMaker2 in the targets-pathways network (Figure 4C). Modules 1 consisted of six pathways, including the TNF signaling pathway (has04668), apoptosis (has04210), necroptosis (has04217), IL-17 signaling pathway (has04657), Toll-like receptor signaling pathway (has04620), and NOD-like receptor signaling pathway (has04621). Modules 2 comprised 4 pathway, including the prolactin signaling pathway (has04917), T cell receptor signaling pathway (has04660), FaxO signaling pathway (has04068), and GnRH signaling pathway (has04912). Module 3 consisted of 2 pathways, including the Th1 and Th2 cell differential (has04658) and Jak-STAT signaling pathway (has04630). Module 4 comprised PPAR signaling pathway (has03320). Module 5 comprised 2 pathways, including the MAPK signaling pathway (hsa04010) and Ras signaling pathway (hsa04014).

Molecular Docking

The binding energy can be calculated to predict their affinity. Generally, the lower the binding energy of the ligand and the receptor, the more stable the binding conformation. As presented in Table 2, all binding affinities are lower than -5 kcal/mol, indicating that CoQ10 has a strong binding activity to the seven hub genes, and the binding conformation is stable. Figure 5 shows the binding mode of CoQ10 with the hub targets.

Figure 5

MD Simulation

To further validate these molecular docking simulation results, MD simulations were performed. The RMSD and RMSF of the backbone were used to scrutinize the stability of the structure model of the CoQ10 and its targets complex. The trajectory files in MD simulation were made as movies (Supplementary 5).

As shown in Figure 6, the RMSD value fluctuates within 0.2 nm after 20 ns, and the system gradually reaches equilibrium. In addition, the average fluctuations of the protein residues were evaluated as the RMSF during the MD simulation (Figure 7).

Figure 6

Figure 7

According to the RMSD values, we extracted the 20-ns to 30-ns trajectory for analysis. Van Der Waals interaction, electronic energy, polar solvation energy, solvent-accessible, and surface area (SASA) energy between CoQ10 and targets were calculated using g-MMPB/SA software. The results were shown in Table 3. It is observed that the van der Waals interaction, electronic energy, and SASA energy promote the binding of CoQ10 to the target protein, followed by electronic energy and SASA energy, while the polar solvation energy have the opposite effect.

To further explore the conformation changes of CoQ10 and its targets complex, we constructed the free energy landscapes of the complex and defined two reaction coordinates of the free landscape map: one is RMSD, which shows the stability of the protein structure; the other is the radius of gyrate (Rg), which reflects whether the complex is folded stably. The energy landscape was calculated along these two reaction coordinates using the gmx sham tools (Figure 8).

Figure 8

As shown in Figure 8, the blue area indicates the minimum value of free energy, while the cyan and green areas indicate the metastable conformational state. For all complexes, a higher blue color is observed, which represents that these complexes were stable.

The MD simulation study further evaluated the stability of the interaction between the hub targets and the CoQ10 ligand predicted by the docking experiment. In addition, free energy calculation provides an assessment of the binding affinity of hub targets to CoQ10. In our study, MD simulations showed that the structures of all complexes are very stable. Among them, the PPARA-CoQ10 complex has the lowest binding free energy, indicating that its binding affinity is the highest and may play a major role in CoQ10’s anti-aging effect on oocytes.

Discussion

The reproductive potential in women gradually declines but significantly beginning approximately in their early 30s. However, many young women choose to postpone marriage and/or childbirth driven by social trends. Aging-related deterioration of oocyte quality is one of the great challenges in reproductive medicine. As mentioned previously, increasing evidence has shown that CoQ10 effectively improve oocyte quality and counteract oocyte aging, but the mechanisms underlying these benefits are still not fully understood. In recent years, the emergence of scRNA-seq technology has dramatically expanded our knowledge of normal and disease-related physiological processes (30, 31). Therefore, in the present study, for the first time, we integrated scRNA-seq data from younger (≤30 years), and older (≥40 years) female oocytes to obtain oocyte aging-related genes. Next, a comprehensive network pharmacology method was utilized to elucidate the mechanism of action of CoQ10 against oocyte aging. According to the network pharmacology results, PPARA, CAT,MAPK14,SQSTM1,HMOX1,GRB2, and GSR play vital roles in protecting against oocyte aging via CoQ10. In addition, the molecular docking of the hub genes and CoQ10 exhibited high affinities. Moreover, MD simulations showed that the structures of all complexes (seven hub genes and CoQ10 complexes) are very stable, implying that the seven hub genes may be highly correlated in against oocyte aging with CoQ10.

The Hub Targets of CoQ10’s in Countering Oocyte Aging

PPARA (also known as PPARα) and two additional PPAR isotypes, PPARβ and PPARγ, are a family of nuclear hormone receptors belonging to the steroid receptor superfamily. PPARα is the first PPAR identified by Issemann and Green in 1990 (32). Since then, it was learned that all the three PPAR isotypes are expressed in the hypothalamic-pituitary-gonadal axis and play a critical role in regulating the proliferation of the different ovarian cells, gametogenesis, steroidogenesis, ovulation, and corpus luteum regression (33, 34). For example, a study using PPARα (−/−) knockout mice has shown that PPARα-dependent pathways play a key role in mediating DEHP’s regulation of ovarian ESR1 expression (35). In addition, the pan-agonist of PPARα, bezafibrate, prevents porcine oocytes (in in vitro-matured) aging through its antioxidant properties and mitochondrial protection (25). Tiefenbach, et al. found that CoQ10 and its derivatives serve as direct partial agonists for both PPARγ and PPARα (36). Consistently, our molecular docking results showed that PPARα has a good binding ability with CoQ10. Further MD simulations showed that the PPARA-CoQ10 complex has the lowest binding free energy, indicating that its binding affinity is the highest, and may play a major role in CoQ10’s anti-aging effect on oocytes.

CAT, which encodes catalase, is a critical endogenous antioxidant enzyme defence against oxidative stress (37, 38). Oxidative stress is caused by the imbalance between the production and destruction of reactive oxygen species (ROS), directly damages the cellular proteins, lipids, and DNA, and thus accelerates oocyte apoptosis and deteriorates oocyte quality (39–42). CoQ10 is one of the representative antioxidants. Several studies have shown that CoQ10 treatment significantly increased the CAT activity, improved preantral follicles viability, and protected ovarian reserve by inhibiting oxidative stress (43–45).

MAPK14 (also known as p38MAPKα) is a member of the MAP kinase family. It is known that activated MAPK can regulate the expression of gonadotropin synthesis genes LHβ and FSHβ, and induce follicle growth and ovulation (67, 68).However, numerous studies have shown that specific MAPK14 inhibitors could block FSH-induced MAPK 3/1 phosphorylation in granulosa cells (GCs). Consistently, the MAPK14 inhibitor SB203580, blocks cumulus cell-oocyte (COC) expansion and inhibits the induction of genes (FSH, AREG, and PGE2) critical for this process in a dose-dependent manner (46). In summary, MAPK14 plays an essential role in COC expansion, oocyte maturation, and steroidogenesis.

Autophagy is a ubiquitous dysfunctional cell components degradation process that plays a critical role in cellular homeostasis (47). SQSTM1, also known as p62, is a classical receptor of autophagy, and its role in regulating oxidative stress and promoting oocyte longevity has attracted widespread attention (48–50). Recent evidence indicates that endogenous SQSTM1/p62 is present as gel-like structures, which serves as a platform for autophagosome formation and anti-oxidative stress response (49). In the experiment of spermidine-induced cytoprotective autophagy of female germline stem cells (FGSCs), Yuan et al. found that spermidine, via up-regulating the expression of SQSTM1/p62, ameliorate the cellular senescence of FGSCs caused by OS (51).

GCs have long been known to provide the physical support and microenvironment required for oocyte ​development. Cumulative evidence has shown that HMOX1, commonly HO-1, has important antioxidant and antiapoptotic in human and animal ovarian GCs (52, 53). GRB2 is a critical molecule in intracellular signal transduction that participates in organogenesis, development, and induction of meiosis in oocytes (54, 55). GSR is a key member of the glutathione antioxidant defense system. Studies have shown that the down-regulation of antioxidant genes, including GSR, is a unique aging feature of early-stage oocytes (56, 57).

In summary, CoQ10 may act on these hub targets to fight against oocytes aging as verified by the molecular docking and MD Simulation.

Important Pathways and Functional Modules of CoQ10’s Putative Targets

The molecular mechanisms of oocyte aging are greatly complicated, involving various biological processes and pathways. The 27 therapeutic targets of CoQ10 against oocyte aging screened in this study mainly participate in the oxidation-reduction process, mitochondrion, enzyme binding, reactive oxygen species metabolic process, ATP binding, electron transfer activity, and membrane raft et al. Notably, of all biological processes, oxidative stress and mitochondrial-related biological functions seem to be the most relevant.

Mitochondria are the most abundant organelles in mammalian oocytes. They act as energy factories in the cells and provide energy for oocyte meiosis maturation, fertilization, and embryonic development through oxidative phosphorylation (OXPHOS). In addition, the mitochondrial genome only comes from the mother, mitochondrial DNA (mtDNA) copy number can be used as a potential biomarker to predict oocyte quality and embryo viability (58, 59). Therefore, reproductive specialists proposed that mitochondrial microinjection is the most promising therapy to restore the quality of oocytes and overcome age-related infertility (60). As shown in Figure 9A, OXPHOS involves the action of the mitochondrial electron transport chain (ETC), which consists of five protein complexes located on the inner mitochondrial membrane. As an essential component of the ETC and ATP synthesis, CoQ10 acts as an electron and proton transporter in the ETC, transferring electrons derived from complexes I and II to complex III. Complex III pumps protons across the membrane and transfers to Cyt c. Then, complex IV transfers the electrons received from Cyt c to oxygen, consuming 2 protons on the matrix side and transferring 2 protons to the membrane space simultaneously, generating an electrochemical proton gradient. At last, ATP synthase uses this proton gradient to form ATP from ADP and phosphate, providing energy for all oocyte events.

Figure 9

According to the free radical theory, oxidative stress caused by elevated intracellular levels of ROS is the most significant contributor to oocyte aging. However, mitochondria are the principal source of cellular ROS, especially through electron leakage from Complexes I and III, and the mitochondria themselves are the primary targets of ROS detrimental effects. In addition, mtDNA is particularly vulnerable to ROS-mediated oxidative damage due to the lack of histone protection. ROS is mainly produced in the matrix side of the inner mitochondrial membrane and overlaps with mtDNA positions. This condition promotes the formation of mtDNA-protein cross-linking mediated by ROS, which is one of the most serious forms of DNA damage. In turn, mtDNA damage directly impairs electron transfer chain and ROS production in the mitochondria, forming a “ROS vicious circle” (Figure 9B). In the clinic, age-related increases in ROS have been found in the follicular fluid of women, while the levels of antioxidants (CAT, SOD, and GPX) in ovaries are reduced with female age (61, 62). The imbalance between ROS and antioxidants leads to a decline in oocyte quality, an essential factor affecting oocyte quality. Fortunately, CoQ10, as a lipid-soluble antioxidant, can scavenge ROS and inhibit DNA oxidation.

Interestingly, there is a tendency to CoQ10 in a variety of mammalian tissues to be reduced with age (especially after the age of 30 in humans) (63). This age-related decline in the CoQ10 level seems to coincide with the time when female fertility began to decline, suggesting a contribution of the reduced CoQ10 level to ovarian aging. To date, lots of studies have revealed that CoQ10 supplements may increase ATP content in oocytes, improve oocyte quality, and significantly delays reproductive aging (16).

Furthermore, the KEGG pathway analysis indicated that the PPAR signaling pathway was the most enriched. Other related pathways include the TNF signaling pathway, Apoptosis, Necroptosis, Prolactin signaling pathway, MAPK signaling pathway, etc. Previous studies have shown that the PPAR signaling pathway is involved in activating the most abundant fertility reserve in female mammals-dormant primordial follicles (64). In addition, between the delicate and intricate game of equilibrium of antithetic peroxidation and antioxidation forces, the PPAR signaling pathway stands out as a central player devoted to against oxidative stress (65). Furthermore, the PPAR signaling pathway also plays a key role in mitochondrial function. In the diabetes-induced atrial, the agonist of PPAR-γ significantly increases mitochondrial biogenesis-related transcription factors (PGC1α, NRF1–2, TFAM) and mtDNA copy number and reduces mitochondrial ROS production (66).

To further understand the CoQ10 mechanisms in improving oocyte quality, five densely linked functional modules were identified based on the target pathway network, as shown in Figure 6C. These modules may be particularly relevant for CoQ10 against oocyte aging, as proteins in a given module may collectively participate in specific biological functions. The first module consists of pathways in the immune system, cell growth and death, and related signaling pathways. The second module consists of pathways in the endocrine and immune system. The third module includes pathways in signal transduction and the immune system, and the fourth mainly involves the PPAR signal pathway.

The last module includes pathways in the signal transduction related to GCs proliferation and oocyte meiosis. Among them, three-fifths of the modules involve the immune system. Specifically, immune system-associated cells and molecules exist in the hypothalamus-pituitary–ovarian axis, regulating the development of fetal and adult primordial follicles. Moreover, Bukovsky et al. demonstrated that immune cells and their secreting factors are an important part of the ovarian germline stem cells (OSC) niche. They participate in the asymmetric division of OSCs and promote the division, proliferation and migration of germ cells and the growth of follicles (67, 68). These modules reflected CoQ10’s effects on endocrine and immune regulation, anti-apoptosis, and GCs proliferation and oocyte meiosis.

Conclusions

In the present work, we found that the potential mechanisms of CoQ10 in improving aging-related deterioration of oocyte quality are closely related to antioxidant, mitochondrial function enhancement, autophagy, endocrine and immune system regulation, and anti-apoptosis (Figure 10). These biological behaviors are regulated by the PPAR, TNF, apoptosis, necroptosis, prolactin, and MAPK signaling pathway, among others. The molecular docking studies and the dynamic simulations verified the potential of CoQ10 as a promising anti-aging agent on oocytes.

Figure 10

Funding

This work was financially supported through grants from the National Natural Science Foundation of China (82004003), the Project of Zhejiang Province Scientific Research Foundation (2020ZA078), Science and Technology Projects of Zhejiang Province (2019C03086), the Project of Zhejiang Province Scientific Research Foundation (2022ZQ061), the Medical and Health Science and Technology Plan of Zhejiang Province (2021KY920), and Zhejiang Zhangqin Famous Traditional Chinese Medicine Expert Inheritance Studio Project (GZS2017014).

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.

References

• 1

  AmanvermezRTosunM. An Update on Ovarian Aging and Ovarian Reserve Tests. Int J Fertil Steril (2016) 9(4):411–5. doi: 10.22074/ijfs.2015.4591

  • CrossRef
  • Google Scholar

• 2

  WangJChenLZhouCWangLXieHXiaoYet al. Prospective Chromosome Analysis of 3429 Amniocentesis Samples in China Using Copy Number Variation Sequencing. Am J Obstet Gynecol (2018) 219(3):287.e281–287.e218. doi: 10.1016/j.ajog.2018.05.030

  • CrossRef
  • Google Scholar

• 3

  ZhuQChenQWangLLuXLyuQWangYet al. Live Birth Rates in the First Complete IVF Cycle Among 20 687 Women Using a Freeze-All Strategy. Hum Reprod (2018) 33(5):924–9. doi: 10.1093/humrep/dey044

  • CrossRef
  • Google Scholar

• 4

  ChenYWangXLiLLuSZhangZ. New Cut-Off Values for Screening of Trisomy 21, 18 and Open Neural Tube Defects (ONTD) During the Second Trimester in Pregnant Women With Advanced Maternal Age. BMC Pregnancy Childbirth (2020) 20(1):776. doi: 10.1186/s12884-020-03464-z

  • CrossRef
  • Google Scholar

• 5

  SunderamSKissinDMCrawfordSBFolgerSGJamiesonDJBarfieldWD. Assisted Reproductive Technology Surveillance–United State. MMWR Surveill Summ (2014) 63(10):1–28. doi: 10.15585/mmwr.ss6606a1

  • CrossRef
  • Google Scholar

• 6

  SmithHLStevensAMinogueBSneddonSShawLWoodLet al. Systems Based Analysis of Human Embryos and Gene Networks Involved in Cell Lineage Allocation. BMC Genomics (2019) 20(1):171. doi: 10.1186/s12864-019-5558-8

  • CrossRef
  • Google Scholar

• 7

  XuRLiCLiuXGaoS. Insights Into Epigenetic Patterns in Mammalian Early Embryos. Protein Cell (2021) 12(1):7–28. doi: 10.1007/s13238-020-00757-z

  • CrossRef
  • Google Scholar

• 8

  YanezLZHanJBehrBBPeraRARCamarilloDB. Human Oocyte Developmental Potential is Predicted by Mechanical Properties Within Hours After Fertilization. Nat Commun (2016) 7(1):10809. doi: 10.1038/ncomms10809

  • CrossRef
  • Google Scholar

• 9

  BeaujouanE. Latest-Late Fertility? Decline and Resurgence of Late Parenthood Across the Low-Fertility Countries. Population Dev Rev (2020) 46(2):219–47. doi: 10.1111/padr.12334

  • CrossRef
  • Google Scholar

• 10

  Ben MessaoudKBouyerJde la RochebrochardE. Infertility Treatment in Franc-2017: A Challenge of Growing Treatment Needs at Older Ages. Am J Public Health (2020) 110(9):1418–20. doi: 10.2105/ajph.2020.305781

  • CrossRef
  • Google Scholar

• 11

  BentovYHannamTJurisicovaAEsfandiariNCasperRF. Coenzyme Q10 Supplementation and Oocyte Aneuploidy in Women Undergoing IVF-ICSI Treatment. Clin Med Insights Reprod Health (2014) 8:31–6. doi: 10.4137/cmrh.S14681

  • CrossRef
  • Google Scholar

• 12

  GiannubiloSROrlandoPSilvestriSCirilliIMarcheggianiFCiavattiniAet al. CoQ10 Supplementation in Patients Undergoing IVF-ET: The Relationship With Follicular Fluid Content and Oocyte Maturity. Antioxid (Basel) (2018) 7(10):141. doi: 10.3390/antiox7100141

  • CrossRef
  • Google Scholar

• 13

  XuYNisenblatVLuCLiRQiaoJZhenXet al. Pretreatment With Coenzyme Q10 Improves Ovarian Response and Embryo Quality in Low-Prognosis Young Women With Decreased Ovarian Reserve: A Randomized Controlled Trial. Reprod Biol Endocrinol (2018) 16(1):29. doi: 10.1186/s12958-018-0343-0

  • CrossRef
  • Google Scholar

• 14

  ZhangYZhangCShuJGuoJChangH-MLeungPCKet al. Adjuvant Treatment Strategies in Ovarian Stimulation for Poor Responders Undergoing IVF: A Systematic Review and Network Meta-Analysis. Hum Reprod Update (2020) 26(2):247–63. doi: 10.1093/humupd/dmz046

  • CrossRef
  • Google Scholar

• 15

  TuriAGiannubiloSRBrugèFPrincipiFBattistoniSSantoniFet al. Coenzyme Q10 Content in Follicular Fluid and its Relationship With Oocyte Fertilization and Embryo Grading. Arch Gynecol Obstet (2012) 285(4):1173–6. doi: 10.1007/s00404-011-2169-2

  • CrossRef
  • Google Scholar

• 16

  Ben-MeirABursteinEBorrego-AlvarezAChongJWongEYavorskaTet al. Coenzyme Q10 Restores Oocyte Mitochondrial Function and Fertility During Reproductive Aging. Aging Cell (2015) 14(5):887–95. doi: 10.1111/acel.12368

  • CrossRef
  • Google Scholar

• 17

  BootsCEBoudouresAZhangWDruryAMoleyKH. Obesity-Induced Oocyte Mitochondrial Defects are Partially Prevented and Rescued by Supplementation With Co-Enzyme Q10 in a Mouse Model. Hum Reprod (2016) 31(9):2090–7. doi: 10.1093/humrep/dew181

  • CrossRef
  • Google Scholar

• 18

  ZhangMShiYangXZhangYMiaoYChenYCuiZet al. Coenzyme Q10 Ameliorates the Quality of Postovulatory Aged Oocytes by Suppressing DNA Damage and Apoptosis. Free Radic Biol Med (2019) 143:84–94. doi: 10.1016/j.freeradbiomed.2019.08.002

  • CrossRef
  • Google Scholar

• 19

  XingXZhangJZhangJWangYWangJKangJet al. Coenzyme Q10 Supplement Rescues Postovulatory Oocyte Aging by Regulating SIRT4 Expression. Curr Mol Pharmacol (2021) 5:190–203. doi: 10.2174/1874467214666210420112819

  • CrossRef
  • Google Scholar

• 20

  LiRGuoCLiYLiangXYangLHuangW. Therapeutic Target and Molecular Mechanism of Vitamin C-Treated Pneumonia: A Systematic Study of Network Pharmacology. Food Funct (2020) 11(5):4765–72. doi: 10.1039/d0fo00421a

  • CrossRef
  • Google Scholar

• 21

  XieBGengQXuJLuHLuoHHuYet al. The Multi-Targets Mechanism of Hydroxychloroquine in the Treatment of Systemic Lupus Erythematosus Based on Network Pharmacology. Lupus (2020) 29(13):1704–11. doi: 10.1177/0961203320952541

  • CrossRef
  • Google Scholar

• 22

  YangLXuHChenYMiaoCZhaoYXingYet al. Melatonin: Multi-Target Mechanism Against Diminished Ovarian Reserve Based on Network Pharmacology. Front Endocrinol (Lausanne) (2021) 12:630504. doi: 10.3389/fendo.2021.630504

  • CrossRef
  • Google Scholar

• 23

  ZhangJJLiuXChenLZhangSZhangXHaoCet al. Advanced Maternal Age Alters Expression of Maternal Effect Genes That are Essential for Human Oocyte Quality. Aging (Albany NY) (2020) 12(4):3950–61. doi: 10.18632/aging.102864

  • CrossRef
  • Google Scholar

• 24

  BurleySKBhikadiyaCBiCBittrichSChenLCrichlowGVet al. RCSB Protein Data Bank: Powerful New Tools for Exploring 3D Structures of Biological Macromolecules for Basic and Applied Research and Education in Fundamental Biology, Biomedicine, Biotechnology, Bioengineering and Energy Sciences. Nucleic Acids Res (2021) 49(D1):D437–51. doi: 10.1093/nar/gkaa1038

  • CrossRef
  • Google Scholar

• 25

  KimJYZhouDCuiXS. Bezafibrate Prevents Aging in In Vitro-Matured Porcine Oocytes. J Anim Sci Technol (2021) 63(4):766–77. doi: 10.5187/jast.2021.e64

  • CrossRef
  • Google Scholar

• 26

  TrottOOlsonAJ. AutoDock Vina: Improving the Speed and Accuracy of Docking With a New Scoring Function, Efficient Optimization, and Multithreading. J Comput Chem (2010) 3(2):455–61. doi: 10.1002/jcc.21334

  • CrossRef
  • Google Scholar

• 27

  PállSZhmurovABauerPAbrahamMLundborgMGrayAet al. Heterogeneous Parallelization and Acceleration of Molecular Dynamics Simulations in GROMACS. J Chem Phys (2020) 153(13):134110. doi: 10.1063/5.0018516

  • CrossRef
  • Google Scholar

• 28

  LeeJChengXSwailsJMYeomMSEastmanPKLemkulJAet al. CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field. J Chem Theory Comput (2016) 12(1):405–13. doi: 10.1021/acs.jctc.5b00935

  • CrossRef
  • Google Scholar

• 29

  Sousa da SilvaAWVrankenWF. ACPYPE - AnteChamber PYthon Parser InterfacE. BMC Res Notes (2012) 5:367. doi: 10.1186/1756-0500-5-367

  • CrossRef
  • Google Scholar

• 30

  TangXHuangYLeiJLuoHZhuX. The Single-Cell Sequencing: New Developments and Medical Applications. Cell Biosci (2019) 9(1):53. doi: 10.1186/s13578-019-0314-y

  • CrossRef
  • Google Scholar

• 31

  LaHYooHLeeEJThangNXChoiHJOhJet al. Insights From the Applications of Single-Cell Transcriptomic Analysis in Germ Cell Development and Reproductive Medicine. Int J Mol Sci (2021) 22(2):863. doi: 10.3390/ijms22020823

  • CrossRef
  • Google Scholar

• 32

  IssemannIGreenS. Activation of a Member of the Steroid Hormone Receptor Superfamily by Peroxisome Proliferators. Nature (1990) 347(6294):645–50. doi: 10.1038/347645a0

  • CrossRef
  • Google Scholar

• 33

  KomarCM. Peroxisome Proliferator-Activated Receptors (PPARs) and Ovarian Function–Implications for Regulating Steroidogenesis, Differentiation, and Tissue Remodeling. Reprod Biol Endocrinol (2005) 3:41. doi: 10.1186/1477-7827-3-41

  • CrossRef
  • Google Scholar

• 34

  VittiMDi EmidioGDi CarloMCartaGAntonosanteAArtiniPGet al. Peroxisome Proliferator-Activated Receptors in Female Reproduction and Fertility. PPAR Res (2016) 2016:4612306. doi: 10.1155/2016/4612306

  • CrossRef
  • Google Scholar

• 35

  KawanoMQinXYYoshidaMFukudaTNansaiHHayashiYet al. Peroxisome Proliferator-Activated Receptor α Mediates Di-(2-Ethylhexyl) Phthalate Transgenerational Repression of Ovarian Esr1 Expression in Female Mice. Toxicol Lett (2014) 228(3):235–40. doi: 10.1016/j.toxlet.2014.04.019

  • CrossRef
  • Google Scholar

• 36

  TiefenbachJMagomedovaLLiuJReunovAATsaiREappenNSet al. Idebenone and Coenzyme Q(10) are Novel Pparα/γ Ligands, With Potential for Treatment of Fatty Liver Diseases. Dis Model Mech (2018) 11(9):1–10. doi: 10.1242/dmm.034801

  • CrossRef
  • Google Scholar

• 37

  IghodaroOMAkinloyeOA. First Line Defence Antioxidants-Superoxide Dismutase (SOD), Catalase (CAT) and Glutathione Peroxidase (GPX): Their Fundamental Role in the Entire Antioxidant Defence Grid. Alexandria J Med (2018) 54(4):287–93. doi: 10.1016/j.ajme.2017.09.001

  • CrossRef
  • Google Scholar

• 38

  NandiAYanLJJanaCKDasN. Role of Catalase in Oxidative Stress- and Age-Associated Degenerative Diseases. Oxid Med Cell Longev (2019) 2019:9613090. doi: 10.1155/2019/9613090

  • CrossRef
  • Google Scholar

• 39

  WangSHeGChenMZuoTXuWLiuX. The Role of Antioxidant Enzymes in the Ovaries. Oxid Med Cell Longev (2017) 2017:4371714. doi: 10.1155/2017/4371714

  • CrossRef
  • Google Scholar

• 40

  SohelMMHAkyuzBKoncaYArslanKSariozkanSCinarMU. Oxidative Stress Modulates the Expression of Apoptosis-Associated microRNAs in Bovine Granulosa Cells. Vitro Cell Tissue Res (2019) 376(2):295–308. doi: 10.1007/s00441-019-02990-3

  • CrossRef
  • Google Scholar

• 41

  KowalskaMPiekutTPrendeckiMSodelAKozubskiWDorszewskaJ. Mitochondrial and Nuclear DNA Oxidative Damage in Physiological and Pathological Aging. DNA Cell Biol (2020) 39(8):1410–20. doi: 10.1089/dna.2019.5347

  • CrossRef
  • Google Scholar

• 42

  YangLChenYLiuYXingYMiaoCZhaoYet al. The Role of Oxidative Stress and Natural Antioxidants in Ovarian Aging. Front Pharmacol (2021) 11:617843(2364). doi: 10.3389/fphar.2020.617843

  • CrossRef
  • Google Scholar

• 43

  HosseinzadehEZavarehSLashkarboloukiT. Antioxidant Properties of Coenzyme Q10-Pretreated Mouse Pre-Antral Follicles Derived From Vitrified Ovaries. J Obstet Gynaecol Res (2017) 43(1):140–8. doi: 10.1111/jog.13173

  • CrossRef
  • Google Scholar

• 44

  SangsefidiZSYaghoubiFHajiahmadiSHosseinzadehM. The Effect of Coenzyme Q10 Supplementation on Oxidative Stress: A Systematic Review and Meta-Analysis of Randomized Controlled Clinical Trials. Food Sci Nutr (2020) 8(4):1766–76. doi: 10.1002/fsn3.1492

  • CrossRef
  • Google Scholar

• 45

  KiremitliTKiremitliSAkselimBYilmazBMammadovRTorIHet al. Protective Effect of Coenzyme Q10 on Oxidative Ovarian and Uterine Damage Induced by Methotrexate in Rats. Hum Exp Toxicol (2021) 40(9):1537–44. doi: 10.1177/09603271211002891

  • CrossRef
  • Google Scholar

• 46

  LiuZFanHYWangYRichardsJS. Targeted Disruption of Mapk14 (P38mapkalpha) in Granulosa Cells and Cumulus Cells Causes Cell-Specific Changes in Gene Expression Profiles That Rescue COC Expansion and Maintain Fertility. Mol Endocrinol (2010) 24(9):1794–804. doi: 10.1210/me.2010-0086

  • CrossRef
  • Google Scholar

• 47

  GalluzziLBaehreckeEHBallabioABoyaPBravo-San PedroJMCecconiFet al. Molecular Definitions of Autophagy and Related Processes. EMBO J (2017) 36(13):1811–36. doi: 10.15252/embj.201796697

  • CrossRef
  • Google Scholar

• 48

  LinF-HZhangW-LLiHTianX-DZhangJLiXet al. Role of Autophagy in Modulating Post-Maturation Aging of Mouse Oocytes. Cell Death Dis (2018) 9(3):308. doi: 10.1038/s41419-018-0368-5

  • CrossRef
  • Google Scholar

• 49

  KageyamaSGudmundssonSRSouY-SIchimuraYTamuraNKazunoSet al. P62/SQSTM1-Droplet Serves as a Platform for Autophagosome Formation and Anti-Oxidative Stress Response. Nat Commun (2021) 12(1):16. doi: 10.1038/s41467-020-20185-1

  • CrossRef
  • Google Scholar

• 50

  PetersAECabanSJMcLaughlinEARomanSDBromfieldEGNixonBet al. The Impact of Aging on Macroautophagy in the Pre-Ovulatory Mouse Oocyte. Front Cell Dev Biol (2021) 9:691826. doi: 10.3389/fcell.2021.691826

  • CrossRef
  • Google Scholar

• 51

  YuanXTianGGPeiXHuXWuJ. Spermidine Induces Cytoprotective Autophagy of Female Germline Stem Cells In Vitro and Ameliorates Aging Caused by Oxidative Stress Through Upregulated Sequestosome-1/P62 Expression. Cell Biosci (2021) 11(1):107. doi: 10.1186/s13578-021-00614-4

  • CrossRef
  • Google Scholar

• 52

  WangYYangCElsheikhNAHLiCYangFWangGet al. HO-1 Reduces Heat Stress-Induced Apoptosis in Bovine Granulosa Cells by Suppressing Oxidative Stress. Aging (Albany NY) (2019) 11(15):5535–47. doi: 10.18632/aging.102136

  • CrossRef
  • Google Scholar

• 53

  ZouLChengGXuCLiuHWangYLiNet al. Copper Nanoparticles Induce Oxidative Stress via the Heme Oxygenase 1 Signaling Pathway In Vitro Studies. Int J Nanomed (2021) 16:1565–73. doi: 10.2147/ijn.S292319

  • CrossRef
  • Google Scholar

• 54

  CailliauKBrowaeys-PolyEBroutin-L'HermiteINiochePGarbayCDucruixAet al. Grb2 Promotes Reinitiation of Meiosis in Xenopus Oocytes. Cell Signal (2001) 13(1):51–5. doi: 10.1016/s0898-6568(00)00138-8

  • CrossRef
  • Google Scholar

• 55

  WangCQinHZhaoCYangLYuTZhangYet al. Whole-Genome Re-Sequencing and Transcriptome Reveal Oogenesis-Related Genes in Autotetraploid Carassius Auratus. Mar Biotechnol (NY) (2021) 23(2):233–41. doi: 10.1007/s10126-021-10018-7

  • CrossRef
  • Google Scholar

• 56

  WangSZhengYLiJYuYZhangWSongMet al. Single-Cell Transcriptomic Atlas of Primate Ovarian Aging. Cell (2020) 180(3):585–600.e519. doi: 10.1016/j.cell.2020.01.009

  • CrossRef
  • Google Scholar

• 57

  ZhangSZhuDMeiXLiZLiJXieMet al. Advances in Biomaterials and Regenerative Medicine for Primary Ovarian Insufficiency Therapy. Bioact Mater (2021) 6(7):1957–72. doi: 10.1016/j.bioactmat.2020.12.008

  • CrossRef
  • Google Scholar

• 58

  WellsD. Mitochondrial DNA Quantity as a Biomarker for Blastocyst Implantation Potential. Fertil Steril (2017) 108(5):742–7. doi: 10.1016/j.fertnstert.2017.10.007

  • CrossRef
  • Google Scholar

• 59

  CecchinoGNGarcia-VelascoJA. Mitochondrial DNA Copy Number as a Predictor of Embryo Viability. Fertil Steril (2019) 111(2):205–11. doi: 10.1016/j.fertnstert.2018.11.021

  • CrossRef
  • Google Scholar

• 60

  LabartaEde Los SantosMJEscribáMJPellicerAHerraizS. Mitochondria as a Tool for Oocyte Rejuvenation. Fertil Steril (2019) 111(2):219–26. doi: 10.1016/j.fertnstert.2018.10.036

  • CrossRef
  • Google Scholar

• 61

  OyawoyeOAbdel GadirAGarnerAConstantinoviciNPerrettCHardimanP. Antioxidants and Reactive Oxygen Species in Follicular Fluid of Women Undergoing IVF: Relationship to Outcome. Hum Reprod (2003) 18(11):2270–4. doi: 10.1093/humrep/deg450

  • CrossRef
  • Google Scholar

• 62

  TatoneCCarboneMCFaloneSAimolaPGiardinelliACasertaDet al. Age-Dependent Changes in the Expression of Superoxide Dismutases and Catalase are Associated With Ultrastructural Modifications in Human Granulosa Cells. Mol Hum Reprod (2006) 12(11):655–60. doi: 10.1093/molehr/gal080

  • CrossRef
  • Google Scholar

• 63

  MilesMVHornPSTangPHMorrisonJAMilesLDeGrauwTet al. Age-Related Changes in Plasma Coenzyme Q10 Concentrations and Redox State in Apparently Healthy Children and Adults. Clin Chim Acta (2004) 347(1-2):139–44. doi: 10.1016/j.cccn.2004.04.003

  • CrossRef
  • Google Scholar

• 64

  YoonSYKimRJangHShinDHLeeJISeolDet al. Peroxisome Proliferator-Activated Receptor Gamma Modulator Promotes Neonatal Mouse Primordial Follicle Activation. In Vitro Int J Mol Sci (2020) 21(9):3120. doi: 10.3390/ijms21093120

  • CrossRef
  • Google Scholar

• 65

  LeeC. Collaborative Power of Nrf2 and Pparγ Activators Against Metabolic and Drug-Induced Oxidative Injury. Oxid Med Cell Longev (2017) 2017:1378175. doi: 10.1155/2017/1378175

  • CrossRef
  • Google Scholar

• 66

  ZhangZZhangXMengLGongMLiJShiWet al. Pioglitazone Inhibits Diabetes-Induced Atrial Mitochondrial Oxidative Stress and Improves Mitochondrial Biogenesis, Dynamics, and Function Through the PPAR-γ/PGC-1α Signaling Pathway. Front Pharmacol (2021) 12:658362. doi: 10.3389/fphar.2021.658362

  • CrossRef
  • Google Scholar

• 67

  BukovskyA. Cell Commitment by Asymmetric Division and Immune System Involvement. Prog Mol Subcell Biol (2007) 45:179–204. doi: 10.1007/978-3-540-69161-7_8

  • CrossRef
  • Google Scholar

• 68

  BukovskyACaudleMR. Immunoregulation of Follicular Renewal, Selection, POF, and Menopause In Vivo, vs. Neo-Oogenesis In Vitro, POF and Ovarian Infertility Treatment, and a Clinical Trial. Reprod Biol Endocrinol (2012) 10:97. doi: 10.1186/1477-7827-10-97

  • CrossRef
  • Google Scholar