Skip to main content

Advertisement

In silico profiling for secondary metabolites from Lepidium meyenii (maca) by the pharmacophore and ligand-shape-based joint approach

Article metrics

Abstract

Background

Lepidium meyenii Walpers (maca) is an herb known as a traditional nutritional supplement and widely used in Peru, North America, and Europe to enhance human fertility and treat osteoporosis. The secondary metabolites of maca, namely, maca alkaloids, macaenes, and macamides, are bioactive compounds, but their targets are undefined.

Methods

The pharmacophore-based PharmaDB targets database screening joint the ligand shape similarity-based WEGA validation approach is proposed to predict the targets of these unique constituents and was performed using Discovery Studio 4.5 and PharmaDB. A compounds–targets–diseases network was established using Cytoscape 3.2. These suitable targets and their genes were calculated and analyzed using ingenuity pathway analysis and GeneMANIA.

Results

Certain targets were identified in osteoporosis (8 targets), prostate cancer (9 targets), and kidney diseases (11 targets). This was the first study to identify the targets of these bioactive compounds in maca for cardiovascular diseases (29 targets). The compound with the most targets (46) was an amide alkaloid (MA-24).

Conclusion

In silico target fishing identified maca’s traditional effects on treatment and prevention of osteoporosis, prostate cancer, and kidney diseases, and its potential function of treating cardiovascular diseases, as the most important of this herb’s possible activities.

Background

Lepidium meyenii Walpers (maca) belongs to the brassica (mustard) family and the Lepidium genus, which grows robustly only at altitudes over 4000 m [1]. Maca has three major phenotypes, yellow, red and black, based on its hypocotyl and stem coloration [2]. The underground part of the maca is consumed as a food and as a folk medicine to enhance fertility and sexual behaviors and has multiple bioactivities [3]. Currently, maca is used in nutrition and health care products sold from Peru to North America and Europe [4, 5]. Maca contains abundant valuable nutritional ingredients [6], such as maca alkaloids, macaenes, glucosinolates, sterols, and polyphenols, and other secondary metabolites. The maca alkaloids, especially macamides and macaenes, are the main functional constituents of maca [7, 8]. To date, 31 maca alkaloids and four macaenes have been isolated from L. meyenii. Their structures are shown below (Table 1). All of the macamides, which are found only in maca, are N-benzylamides. Wu [9] synthesized 11 of the reported macamides as well as a series of structurally related amides that resemble macamides (Table 2). These synthesized compounds were collected in this study and used in our experiments.

Table 1 Structures of alkaloids and macaenes isolated from L. meyenii
Table 2 Structures of the synthetic amides resembling macamides

The biological and pharmacological effects of maca have been investigated in experimental animals such as rats [10], mice [11], fish [12], and bulls [13]. Maca can enhance sexual behavior and increase sperm count [14], improve prostate function related to testosterone enanthate (TE) levels [15], and improve the quality of embryos [16, 17]. Maca also has beneficial effects on learning and memory in scopolamine-induced memory impairment mice [18]. Moreover, maca aqueous extract scavenges free radicals and protects cells against oxidative stress [5].

Several in vitro/in vivo animal experimental studies have shown that secondary metabolites extracted from maca have bioactivities that help treat osteoporosis and enhance prostate function and sexual function [1921]. However, most in vitro and in vivo experiments have not specified the molecular target of these secondary metabolites and the mechanisms of the functions of the compounds obtained from maca are unclear. The pharmacophore model is reliable for parallel screening to predict and mimic the binding situation of compounds and targets [22, 23]. This study aimed to investigate the network involved in the mechanisms of action of secondary metabolites of maca. We used the pharmacophore-based method in combination with a novel ligand shape similarity strategy and used the weighted Gaussian algorithm (WEGA).

Methods

Collection of chemical constituents

The natural constituents of maca were collected from the literature [9, 2430] using the search terms “lepidium or maca” combined with “constitutes, compounds, chemical or metabolites.” Traditionally used maca contains a dominant pattern of secondary metabolites, particularly alkaloids and macaenes [31]. The secondary metabolite constituents from maca were evaluated to precisely predict the active compounds. A total of 47 alkaloids extracted from maca and synthetic amides were categorized into classes.

Conformer generation

All chemical structures were prepared in SD format, converted from a 2D cdx file format to 3D models, using Open Babel GUI [32] version 2.3.2 (OpenBableGUI; Chris Morley, USA). Molecular energy was minimized using the Energy Minimization module of Discovery Studio version 4.5 (DS 4.5; Accelrys Inc., San Diego, CA, USA) under the chemistry at Harvard Macromolecular Mechanics (CHARMM) force field. This survey led to the construction of the 3D multi-conformational maca compounds molecular structure database (i.e., maca-DB), which was generated by a Monte Carlo-based conformational analysis (FAST mode). These compounds are rigid; the number of conformers for each compound is much less than 255. The maca-DB contains a total of 47 constituents and 9976 conformations.

Pharmacophore model collection

In silico profiling of the maca-DB was performed using the generated 3D chemical feature-based pharmacophore models. The pharmacophore models were used to represent the binding mode of particular compounds to specific drug targets [33]. Each pharmacophore model contained several convictive chemical features that determine the chemical functionalities of a certain ligand: H-bond donors or acceptors, hydrophobic groups, aromatic nuclei, and positive or negative ionizable moieties [34]. Unlike common docking methods, pharmacophore-based virtual screening outlines the specific compounds and their multiple pharmacologic targets and determines novel actions of these compounds.

PharmaDB

PharmaDB is the only pharmacophore database implemented in DS 4.5. A total of 68,000 pharmacophores were derived from 8000 protein–ligand complexes in the sc-PDB dataset. sc-PDB is designed to identify binding sites suitable for the docking of a drug-like ligand, and 9276 three-dimensional structures of binding sites were identified using the Protein Data Bank (PDB) [35].

Parameters

In this study, PharmaDB was used for profiling. All pharmacophore models with the shape of the binding pocket were selected for virtual screening using the default settings of the Ligand Profiler module of DS 4.5. In principle, each alkaloid that mapped to a chemical feature of the respective pharmacophore model was counted as one hit. The screening was conducted using the default settings and with a minimal inter-feature distance of 0.5 Å.

Binding mode refinement

All the poses of the ligands mapped to the pharmacophore were preserved. A series of target–ligand pairs were selected for further examination. The selection was based on compatibility with previously reported pharmacological activities and the traditional use of maca. Further refinement was carried in Molecular Operating Environment (MOE; Chemical Computing Group Inc., Canada) to identify the protein–ligand binding modes. Energy minimization was performed by conjugated gradient minimization with the Merck Molecular Force Field 94× (MMFF94×) until a root-mean-square deviation of 0.1 kcal mol−1 Ǻ−1 was reached.

WEGA validation

The WEGA is an accurate shape-based virtual screening method [36]. In this research, we validated the reliability of the binding model by calculating the binding efficiencies of the compounds and the original ligands of the hit targets using the shape similarity calculations function of WEGA.

The sc-PDB also provides separate MOL2 files for the ligand, its binding site, and the corresponding protein chain(s). Ions and cofactors at the vicinity of the ligand are included in the protein. This helps to evaluate the influence of ligand binding on binding site diversity for docking. MOL2 files of the hit-target protein ligands were selected to create the target-ligands database (tl-DB). WEGA validation was performed by comparing the contents of the maca-DB to those of the tl-DB.

Network construction

The Table 3 showing interactions between all mapped compounds and hit targets shows the ligand profiling results. For each target, the protein name, gene name, and pathway information were collected from the PDB, Kyoto Encyclopedia of Genes and Genomes (KEGG) [37], and Cell Signaling Technology (CST) [38]. The target–target interactions were mapped using GeneMANIA [39]; all targets were analyzed using Ingenuity Pathway Analysis (IPA®; Qiagen, Redwood City, CA, USA). All diseases related to the targets were retrieved from the Therapeutic Target Database (TTD) [40] and DrugBank database [41]. The overall compound–target–pathway networks were generated using Cytoscape, version 3.2 (Cytoscape Consortium, USA). In the graphical networks, nodes represent the compounds, targets, and related diseases. The edges linking the compound-target and target-diseases represent their relationships and are marked with two types of lines. The related diseases are marked with different colors at the nodes. The targeted diseases pathway was mapped using the KEGGscape plugin of Cytoscape, version 3.2.

Table 3 Four categories of disease targets of the selected compounds from profiling

Results and discussion

Evaluation of constituents

In modern drug discovery, large compound libraries are compared, and the diversity of these libraries must be analyzed [42]. The constituents collected and synthesized from maca could be divided into eight compound classes (Fig. 1). The 40 compounds examined in this research were fished by targets (Fig. 2). The compounds with higher degree values were distributed across different categories, such as amide alkaloids (MA-24; 25), macaenes (MA-32; 33), and synthetic amides (MA-43; 44). Compounds that participate in more interactions than other components have a higher bioactivity value.

Fig. 1
figure1

The diversity of compounds analyzed by the scaffold-based classification approach (SCA). CID means compound class ID, categories by complexity; the SCA also outputs the following structural descriptor values: (1) Cyclicity side chain value; (2) AE average electronegativity; (3) HD number of H-bond donors; (3) HA number of H-bond acceptors; (4) AB number of aromatic bonds; (5) ATMS number of non-H atoms; (6) BNDS number of non-H-involved bonds; (7) SSSRS number of the smallest set of smallest rings; (8) AZ average atomic numbers; (9) RB number of rotating bonds; (10) MW molecular weight

Fig. 2
figure2

Overview of number of hit compounds of the parallel pharmacophore screening against druggable targets in the field of diseases

Network analysis

In total, 950 models were selected for in silico screening of the maca-DB. These models belong to 125 protein targets; 87 of those targets were validated by the TTD database, were involved in 60 pathways, and were targets of 41 maca constituents. As shown in Fig. 2, we chose three major disease areas (prostate cancer, osteoporosis, and kidney diseases) to validate the traditional medical action and the fished maca compounds.

Interpreting the mechanisms of action

An array of well-defined in-house structure and ligand-based pharmacophore models was selected from PharmaDB. For the profiling results, all biological functions of hit targets were annotated from TTD and DrugBank. The identified targets had variable pharmacologic usages, such as the treatment of osteoporosis (8 targets), prostate disease (9 targets), and kidney disease (11 targets). Some targets were related to cardiovascular diseases (CVD), such as hypertension, myocardial infarction, ischemic heart disease, and dyslipidemia (29 targets). Figure 3 and Table 3 provides an overview of the selected targets in the four categories of disease mentioned above. A total of 125 targets were mapped, and the IPA analysis indicated that 107 of them have been used to make drugs. The druggable list is presented in Additional file 1.

Fig. 3
figure3

The major pharmacologic network of maca. Hexagon, compounds; circle, targets (red: osteoporosis; light blue: kidney diseases; deep blue: prostate cancer)

WEGA validation

The WEGA is suitable for large-scale parallel screening of a series of bioactive compounds; regardless of the conformations of the compounds, their targets can be experimentally determined. The results determined by shape showed that all binding models obtained with DS 4.5 had a high ligand-receptor structure binding value: all scores were above 0.5 (Additional file 2). Important molecular superimposition images are shown in Fig. 4.

Fig. 4
figure4

Compounds from maca align with natural ligands from PDB structure (light blue) by WEGA (a matrix metalloproteinases; b androgen receptor; c carbonic anhydrase II; d estrogen receptor)

Selected targets related to prostate cancer

All the gene interactions of these targets were analyzed using GeneMANIA (Fig. 5; Additional file 3). Epidemiological studies have found that consumption of maca could reduce the risk of prostate cancer, which might be associated with aromatic glucosinolate content [7, 8]. Animal experiments in mice [43] and rats [7, 44, 45] showed that maca reduced TE level in a dose-dependent manner and induced prostatic hyperplasia. Red maca aqueous extracts can also reduce ventral prostate size in normal and TE-treated rats [7].

Fig. 5
figure5

Network images for the interaction of targets related to four disease categories fished by profiling (a prostate cancer; b kidney diseases; c osteoporosis; d cardiovascular diseases)

In this study, nine targets annotated by the TTD database were related to prostate cancer. One of the most common targets used in the treatment of prostate diseases is the androgen receptor (AR). In some cell types, testosterone is converted by 5-alpha-reductase into dihydrotestosterone, which is an even more potent agonist for AR activation than testosterone [46]. AR is a sequence-specific DNA-binding protein involved in cellular proliferation in prostate cancer and in the development of secondary sexual characteristics through activation by dihydrotestosterone [47]. AR is also involved in the regulation of the adhesion of prostate cancer cells to the extracellular matrix and to the invasion of prostate cancer cells through its influence on the expression of specific integrin subunits [48]. There is increasing evidence that the genus Lepidium could reduce the risk of prostate cancer development [49, 50]. Research also suggests that other cruciferous plants from the genus Lepidium could be used as important alternative treatments for prostate diseases [51]. Growth of the prostate is a hormone-mediated phenomenon regulated by both androgens and estrogens [52]. A recent report indicates that maca’s effect on ventral prostate size may partly be a result of the action of glucosinolate metabolites on AR [28]. In this study, AR was fished by MA-4, MA-24, and MA-44 (Fig. 5a).

Stromelysin-1, also known as matrix metalloproteinase-3 (MMP-3), is an enzyme that activates other targeted matrix metalloproteinases (MMPs): MMP-9, MMP-12, and MMP-13 [53]. The expression of MMPs is primarily regulated at the level of transcription, where the promoter of the gene responds to various stimuli, including growth factors, cytokines, tumor promoters, and oncogene products [54]. MMPs are associated with various physiological and pathological processes, such as morphogenesis, angiogenesis, tissue repair, cirrhosis, arthritis, and metastasis [55]. MMPs also play a significant role in the development and metastasis of prostate cancer [56]. In particular, MMP-9 has been shown experimentally to be involved in prostate cancer [57]. These four targets share similar protein domains and were fished by MA-18, MA-25, MA-44, and MA-45 (Fig. 5b).

Proto-oncogene tyrosine-protein kinase Src (c-Src) is a non-receptor tyrosine kinase protein that may be involved in cancer progression by promoting other signals [58]. C-Src is highly expressed in malignant prostate cells [59]. When primary prostate cells were treated with a c-Src inhibitor in vitro, their proliferation, migration, and invasive potential were reduced [60].

As c-Src is a solid druggable target, several c-Src tyrosine kinase inhibitors have been utilized therapeutically [61]. Dasatinib has been approved for the treatment of chronic myeloid leukemia [62], which is an Src family inhibitor. Based on the binding results of those compounds that fished c-Src, we can predict that MA-24 and MA-25 may be pharmacologically similar to Dasatinib.

Selected targets related to kidney diseases

Renal tubular acidosis (RTA) is a metabolic acidosis caused by impaired excretion of acid by the kidney. Carbonic anhydrase II (CA2) is one of the 14 forms of human α carbonic anhydrases and the one with the highest catalytic activity. The physiological functions of CA2 include pH regulation, CO2 and H2CO3 transport, and maintenance of H2O and electrolyte balance. CA2 deficiency syndrome can lead to osteoporosis, RTA, and cerebral calcification.

In inherited CA2 deficiency, isolated proximal RTA presents with osteoporosis (owing to impaired osteoclast function), cerebral calcification, and variable levels of mental retardation. Although this form of inherited RTA is clinically more proximal, it can also present with a mixed proximal and distal phenotype, which reflects the presence of CA2 in cells all along the renal tubule. CA2 was fished by MA-18, MA-24, MA-38, MA-44, and MA-45 (Fig. 5c). Kidney status directly affects the reproductive function, especially sexual behavior. Oral administration of a purified lipidic extract from maca could enhance sexual behavior by increasing the number of complete intromissions in normal mice and decreasing the latent period of erection in erectile dysfunction male rats [63].

Selected targets related to osteoporosis

Osteoporosis is a skeletal fragility disorder and is common in elderly people. Its prevalence is increasing as more individuals are developing low bone mineral density [64]. The edible part of maca, the hypocotyl, has been widely used to treat osteoporosis [65]. Ethanol extract of maca has anti-osteoporotic activity and indicated maca alkaloids, steroids, glucosinolates, isothiocyanates, and macamides are probably responsible for its biological functions [19].

Estrogen receptor α (ER-α) binds to estrogens and regulates bone homeostasis and prevents postmenopausal bone loss [66, 67]. Estrogen deficiency is a major determinant of bone loss in postmenopausal women [68, 69]. In one ovariectomized rats experiment, ER-α was the predominant ER form expressed in mesenchymal stem cells [70]. Co-expression of ER-α with other genes indicates its activator function in the osteogenic differentiation of mesenchymal stem cells, which causes osteoporosis [71]. Estradiol, estrone, and raloxifene bind to the alpha receptor. However, because the ER’s helix 12 domain plays a crucial role in determining its interactions with co-activators and co-repressors, different ER combinations may respond differently to various ligands, which may translate into tissue selective agonistic and antagonistic effects [72, 73]. For example, tamoxifen is an antagonist in the breast and is used as a breast cancer treatment, but it is an ER agonist in bone and therefore prevents osteoporosis [69]. Recent studies have suggested that maca contains phyto-estrogens, which may have estrogenic activity [74, 75]. We found that three compounds were connected to ER: MA-19, MA-24, and MA-25 (Fig. 5d).

Enzymatic cleavage by MMPs is involved in the destruction of articular cartilage, and the high expression of MMP-9 and MMP-13 could be detected in pathologic synovium and cartilage samples [76]. Several natural substances containing maca extract tested in vitro are effective agents, as evidenced by the strong regulation of MMP-9 and MMP-13 [77]. In osteoclast migration, MMPs control the cell–matrix interactions required in the model of osteoclast recruitment in primitive long bones [78]. As classical anti-osteoporotic agents, bisphosphonates are involved in the inhibition of the functions of several MMPs (MMP-3, -9, -12, and -13), which were mapped in this virtual screening.

Prediction of the function of maca compounds in the treatment of cardiovascular diseases

A total of 29 targets related to CVD were mapped. Maca could be used in the treatment of CVD characterized by atherogenic lipoprotein profile, and showed relevant angiotensin I-converting enzyme inhibitory activities, indicating potential anti-hypertension activity; however, the mechanisms of these activities are still to be clarified. This result indicated that maca might have significant potential for the treatment of CVD.

Kinase activity mediated by mitogen-activated protein kinase 14 (MAPK14), also called p38α, has been identified in many tissues [79]. p38α is mainly activated through MAPK kinase kinase cascades and exerts its biological function via downstream substrate phosphorylation [80]. Pharmacological and genetic inhibition of p38α has revealed its biological significance regarding physiological functions and its potential for targeting p38α in human diseases, especially CVD [8183]. MAPK14 activity regulates myocyte cytokinesis and promotes cell-cycle exit during maturation in the newborn mouse heart [84]. MAPK14 has also been associated with cell-cycle arrest in mammalian cardiomyocytes [85], and its inhibition might be a strategy to promote cardiac regeneration in response to injury [86]. Furthermore, MAPK14 promoted myocyte apoptosis and cardiomyocyte hypertrophy, and targeted IRS-1-mediated Akt signaling and promoted myocyte death under chronic insulin stimulation in vitro [87, 88].

Peroxisome proliferator-activated receptors (PPARs) are a group of nuclear receptor proteins [89], including PPARα, PPARδ, and PPARγ, whose ligand and DNA-binding domains share 60–80 % homology [90, 91]. PPARs are widely expressed in the vasculature, myocardium, and the immune cells, such as monocytes and macrophages [92]. Additionally, PPAR-retinoid X receptor heterodimers repress CLOCK/BMAL1 gene expression [93]. Hence, PPARs could regulate the expression of a series of genes involved in metabolism that impact cardiovascular physiology [94]. Different PPAR isoforms are observed in various cardiovascular pathologies, such as atherosclerosis, hypertension, and cardiac hypertrophy [95]. Both PPARα and PPARγ are expressed in endothelial cells, vascular smooth muscle cells, and monocytes/macrophages [96, 97]. In atherosclerosis, activation of these two proteins reduces leukocyte recruitment and cell adhesion [98]. Both regulate cytokine-induced genes (such as VCAM-1 and tissue factor), and PPARα and PPARγ inhibit the expressions of tumor necrosis factor-α and MCP-1, respectively [99]. PPARδ activation decreases the expressions of MCP-1, ICAM-1, and inflammatory cytokines and attenuates atherosclerosis development [100].

The potential use of PPAR agonists and dual PPAR agonists, including PPARα/γ, PPARα/δ, and PPARδ/γ dual agonists, in the treatment of CVD has recently received attention [101]. Compounds that are capable of targeting more than one PPAR isotype and are effective at treating CVD have emerged as an interesting and efficient treatment approach. Both MAPK14 and PPARs are related to a series of maca compounds (Table 4).

Table 4 Compounds fished by MAPK14 and PPARs

Conclusion

In silico target fishing identified maca’s traditional effects on treatment and prevention of osteoporosis, prostate cancer, and kidney diseases, and its potential function of treating cardiovascular diseases, as the most important of this herb’s possible activities.

Abbreviations

AR:

androgen receptor

BMD:

bone mineral density

CA2:

carbonic anhydrase II

CHARMM:

chemistry at Harvard macromolecular mechanics

c-Src:

tyrosine-protein kinase Src

CST:

cell signaling technology

CVD:

cardiovascular disease

DS:

discovery studio

ER-α:

estrogen receptor α

IPA:

ingenuity pathway analysis

KEGG:

Kyoto encyclopedia of genes and genomes

LPE:

latent period of erection

MAPK14:

mitogen-activated protein kinase 14

MMFF94x:

Merck molecular force field 94×

MMPs:

matrix metalloproteinases

MOE:

molecular operating environment

MSCs:

mesenchymal stem cells

PDB:

protein data bank

PPARs:

peroxisome proliferator-activated receptors

RMSD:

root mean square deviation

RTA:

renal tubular acidosis

TE:

testosterone enanthate

TTD:

therapeutic target database

WEGA:

weighted Gaussian algorithm

References

  1. 1.

    Hermann M, Heller J. Andean roots and tubers: ahipa, arracacha, maca and yacon. Rome: Inernational Plant Genetic Resources Insititute; 1997. p. 1–199.

  2. 2.

    Esparza E, Hadzich A, Kofer W, Mithöfer A, Cosio EG. Bioactive maca (Lepidium meyenii) alkamides are a result of traditional Andean postharvest drying practices. Phytochemistry. 2015;116:138–48.

  3. 3.

    León J. The, “Maca” (Lepidium meyenii), A little known food plant of Peru. Econ Bot. 1964;18:122–7.

  4. 4.

    Gonzales GF, Singh VK, Govil JN, Ahmad K, Sharma RK. Biological effects of Lepidium meyenii, maca, a plant from the highlands of Peru. Nat Prod. 2007;1:209–34.

  5. 5.

    Sandoval M, Okuhama NN, Angeles FM, Melchor VV, Condezo LA, Lao J, Miller MJS. Antioxidant activity of the cruciferous vegetable Maca (Lepidium meyenii). Food Chem. 2002;79:207–13.

  6. 6.

    Piacente S, Carbone V, Plaza A, Zampelli A, Pizza C. Investigation of the tuber constituents of maca (Lepidium meyenii Walp.). J Agric Food Chem. 2002;50:5621–5.

  7. 7.

    Gonzales GF, Miranda S, Nieto J, Fernández G, Yucra S, Rubio J, Yi P, Gasco M. Red maca (Lepidium meyenii) reduced prostate size in rats. Reprod Biol Endocrinol. 2005;3:1–16.

  8. 8.

    Večeřa R, Orolin J, Škottová N, Kazdová L, Oliyarnik O, Ulrichová J, Šimánek V. The influence of maca (Lepidium meyenii) on antioxidant status, lipid and glucose metabolism in rat. Plant Food Hum Nutr. 2007;62:59–63.

  9. 9.

    Wu H, Kelley CJ, Pino-Figueroa A, Vu HD, Maher TJ. Macamides and their synthetic analogs: evaluation of in vitro FAAH inhibition. Bioorg Med Chem. 2013;21:5188–97.

  10. 10.

    Cicero AFG, Bandieri E, Arletti R. Lepidium meyenii Walp. improves sexual behaviour in male rats independently from its action on spontaneous locomotor activity. J Ethnopharmacol. 2001;75:225–9.

  11. 11.

    Ruiz-Luna AC, Salazar S, Aspajo NJ, Rubio J, Gasco M, Gonzales GF. Lepidium meyenii (Maca) increases litter size in normal adult female mice. Reprod Biol Endocrinol. 2005;3:1–16.

  12. 12.

    Lee KJ, Dabrowski K, Rinchard J, Gomez C, Guz L, Vilchez C. Supplementation of maca (Lepidium meyenii) tuber meal in diets improves growth rate and survival of rainbow trout Oncorhynchus mykiss (Walbaum) alevins and juveniles. Aquac Res. 2004;35:215–23.

  13. 13.

    Clément C, Kneubühler J, Urwyler A, Witschi U, Kreuzer M. Effect of maca supplementation on bovine sperm quantity and quality followed over two spermatogenic cycles. Theriogenology. 2010;74:173–83.

  14. 14.

    Gonzales C, Rubio J, Gasco M, Nieto J, Yucra S, Gonzales GF. Effect of short-term and long-term treatments with three ecotypes of Lepidium meyenii (MACA) on spermatogenesis in rats. J Ethnopharmacol. 2006;103:448–54.

  15. 15.

    Gasco M, Villegas L, Yucra S, Rubio J, Gonzales GF. Dose-response effect of Red Maca (Lepidium meyenii) on benign prostatic hyperplasia induced by testosterone enanthate. Phytomedicine. 2007;14:460–4.

  16. 16.

    Lembè DM, Gasco M, Gonzales GF. Fertility and estrogenic activity of Turraeanthus africanus in combination with Lepidium meyenii (Black maca) in female mice. Eur J Integr Med. 2012;4:e345–51.

  17. 17.

    Rubio J, Qiong W, Liu XM, Jiang Z, Dang HX, Chen SL, Gonzales GF. Aqueous extract of black maca (Lepidium meyenii) on memory impairment induced by ovariectomy in mice. Evid Based Complement Alternat Med. 2011;2011:253958.

  18. 18.

    Rubio J, Dang H, Gong M, Liu XM, Chen SL, Gonzales GF. Aqueous and hydroalcoholic extracts of Black Maca (Lepidium meyenii) improve scopolamine-induced memory impairment in mice. Food Chem Toxicol. 2007;45:1882–90.

  19. 19.

    Rosales-Hartshorn M. Maca: botanical medicine from the Andes. Adv Food Tech Nutr Sci Open J. 2015;1:e1–6.

  20. 20.

    Lentz A, Gravitt K, Carson CC, Marson L. Acute and chronic dosing of Lepidium meyenii (Maca) on male rat sexual behavior. J Sex Med. 2007;4:332–40.

  21. 21.

    Shin BC, Lee MS, Yang EJ, Lim HS, Ernst E. Maca (L. meyenii) for improving sexual function: a systematic review. BMC Complement Altern Med. 2010;10:44–50.

  22. 22.

    Dror O, Shulman-Peleg A, Nussinov R, Wolfson HJ. Predicting molecular interactions in silico: I. A guide to pharmacophore identification and its applications to drug design. Curr Med Chem. 2004;11:71–90.

  23. 23.

    Taha MO, Bustanji Y, Al-Bakri AG, Yousef AM, Zalloum WA, Al-Masri IM, Atallah N. Discovery of new potent human protein tyrosine phosphatase inhibitors via pharmacophore and QSAR analysis followed by in silico screening. J Mol Graph Model. 2007;25:870–84.

  24. 24.

    Muhammad I, Zhao JP, Dunbar DC, Khan IA. Constituents of Lepidium meyenii ‘maca’. Phytochemistry. 2002;59:105–10.

  25. 25.

    Cui BL, Zheng BL, He K, Zheng QY. Imidazole alkaloids from Lepidium meyenii. J Nat Prod. 2003;66:1101–3.

  26. 26.

    Zhao JP, Muhammad I, Dunbar DC, Mustafa J, Khan IA. New alkamides from maca (Lepidium meyenii). J Agric Food Chem. 2005;53:690–3.

  27. 27.

    McCollom MM, Villinski JR, McPhail KL, Craker LE, Gafner S. Analysis of macamides in samples of Maca (Lepidium meyenii) by HPLC-UV-MS/MS. Phytochem Anal. 2005;16:463–9.

  28. 28.

    Tellez MR, Khan IA, Kobaisy M, Schrader KK, Dayan FE, Osbrink W. Composition of the essential oil of Lepidium meyenii (Walp.). Phytochemistry. 2002;61:149–55.

  29. 29.

    Gonzales GF, Gonzales-Castañeda C. The methyltetrahydro-β-carbolines in Maca (Lepidium meyenii). Evid Based Complement Alternat Med. 2009;6:315–6.

  30. 30.

    Ganzera M, Zhao J, Muhammad I, Khan IA. Chemical profiling and standardization of Lepidium meyenii (Maca) by reversed phase high performance liquid chromatography. Chem Pharm Bull. 2002;50:988–91.

  31. 31.

    Gonzales GF. Ethnobiology and ethnopharmacology of Lepidium meyenii (Maca), a plant from the Peruvian highlands. Evid Based Complement Alternat Med. 2012;2012:193496.

  32. 32.

    O’Boyle NM, Banck M, James CA, Morley C, Vandermeersch T, Hutchison GR. Open babel: an open chemical toolbox. J Cheminform. 2011;3:33–46.

  33. 33.

    Guner OF. Pharmacophore perception, development, and use in drug design, vol. 2. La Jolla: International Univeristy Line; 2000.

  34. 34.

    Giannakakou P, Gussio R, Nogales E, Downing KH, Zaharevitz D, Bollbuck B, Poy G, Sackett D, Nicolaou K, Fojo T. A common pharmacophore for epothilone and taxanes: molecular basis for drug resistance conferred by tubulin mutations in human cancer cells. Proc Natl Acad Sci USA. 2000;97:2904–9.

  35. 35.

    http://www.rcsb.org/pdb/home/home.do.

  36. 36.

    Yan X, Li J, Liu Z, Zheng M, Ge H, Xu J. Enhancing molecular shape comparison by weighted Gaussian functions. J Chem Inf Model. 2013;53:1967–78.

  37. 37.

    http://www.kegg.jp.

  38. 38.

    http://www.cellsignal.com/contents/science/cst-pathways/science-pathways.

  39. 39.

    Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT, Maitland A, Mostafavi S, Montojo J, Shao Q, Wright G, Bader GD, Morris Q. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38:W214–20.

  40. 40.

    Zhu F, Shi Z, Qin C, Tao L, Liu X, Xu F, Zhang L, Song Y, Liu XH, Zhang JX, Han BC, Zhang P, Chen YZ. Therapeutic target database update 2012: a resource for facilitating target-oriented drug discovery. Nucleic Acids Res. 2012;40:D1128–36.

  41. 41.

    Knox C, Law V, Jewison T, Liu P, Ly S, Frolkis A, Pon A, Banco K, Mak C, Neveu V, Djoumbou Y, Eisner R, Guo AC, Wishart DS. DrugBank 3.0: a comprehensive resource for ‘omics’ research on drugs. Nucleic Acids Res. 2011;39:D1035–41.

  42. 42.

    Xu J. A new approach to finding natural chemical structure classes. J Med Chem. 2002;45:5311–20.

  43. 43.

    Gonzales GF, Gasco M, Malheiros-Pereira A, Gonzales-Castañeda C. Antagonistic effect of Lepidium meyenii (red maca) on prostatic hyperplasia in adult mice. Andrologia. 2008;40:179–85.

  44. 44.

    Gonzales GF, Vasquez V, Rodriguez D, Maldonado C, Mormontoy J, Portella J, Pajuelo M, Villegas L, Gasco M, Mathur P. Effect of two different extracts of red maca in male rats with testosterone-induced prostatic hyperplasia. Asian J Androl. 2007;9:245–51.

  45. 45.

    Gasco M, Villegas L, Yucra S, Rubio J, Gonzales G. Dose-response effect of Red Maca (Lepidium meyenii) on benign prostatic hyperplasia induced by testosterone enanthate. Phytomedicine. 2007;14:460–4.

  46. 46.

    Davis SR, Worsley R. Androgen treatment of postmenopausal women. J Steroid Biochem. 2014;142:107–14.

  47. 47.

    Fu M, Wang CG, Reutens AT, Wang J, Angeletti RH, Siconolfi-Baez L, Ogryzko V, Avantaggiati ML, Pestell RG. p300 and p300/cAMP-response element-binding protein-associated factor acetylate the androgen receptor at sites governing hormone-dependent transactivation. J Biol Chem. 2000;275:20853–60.

  48. 48.

    Nagakawa O, Akashi T, Hayakawa Y, Junicho A, Koizumi K, Fujiuchi Y, Furuya Y, Matsuda T, Fuse H, Saiki I. Differential expression of integrin subunits in DU-145/AR prostate cancer cells. Oncol Rep. 2004;12:837–41.

  49. 49.

    Azimi H, Khakshur AA, Aghdasi I, Fallah-Tafti M, Abdollahi M. A review of animal and human studies for management of benign prostatic hyperplasia with natural products: perspective of new pharmacological agents. Inflamm Allergy Drug Targets. 2012;11:207–21.

  50. 50.

    Shrivastava A, Gupta VB. Various treatment options for benign prostatic hyperplasia: a current update. J Midlife Health. 2012;3:10–9.

  51. 51.

    Marti S, Fernandez CC, Perez-Fernandez R. Effect of an integral suspension of Lepidium latifolium on prostate hyperplasia in rats. Fitoterapia. 2004;75:187–91.

  52. 52.

    Lund TD, Munson DJ, Adlercreutz H, Handa RJ, Lephart ED. Androgen receptor expression in the rat prostate is down-regulated by dietary phytoestrogens. Reprod Biol Endocrinol. 2004;2:1–5.

  53. 53.

    Egeblad M, Werb Z. New functions for the matrix metalloproteinases in cancer progression. Nat Rev Cancer. 2002;2:161–74.

  54. 54.

    Matrisian LM. Metalloproteinases and their inhibitors in matrix remodeling. Trends Genet. 1990;6:121–5.

  55. 55.

    Eguchi T, Kubota S, Kawata K, Mukudai Y, Uehara J, Ohgawara T, Ibaragi S, Sasaki A, Kuboki T, Takigawa M. Novel transcription factor-like function of human matrix metalloproteinase 3 regulating the CTGF/CCN2 gene. Mol Cell Biol. 2008;28:2391–413.

  56. 56.

    Lokeshwar BL. MMP inhibition in prostate cancer. Ann NY Acad Sci. 1999;878:271–89.

  57. 57.

    Kong DJ, Li YW, Wang ZW, Banerjee S, Sarkar FH. Inhibition of angiogenesis and invasion by 3,3′-diindolylmethane is mediated by the nuclear factor-κB downstream target genes MMP-9 and uPA that regulated bioavailability of vascular endothelial growth factor in prostate cancer. Cancer Res. 2007;67:3310–9.

  58. 58.

    Devary Y, Gottlieb RA, Smeal T, Karin M. The mammalian ultraviolet response is triggered by activation of Src tyrosine kinases. Cell. 1992;71:1081–91.

  59. 59.

    Wheeler DL, Iida M, Dunn EF. The role of Src in solid tumors. Oncologist. 2009;14:667–78.

  60. 60.

    Park SI, Zhang J, Phillips KA, Araujo JC, Najjar AM, Volgin AY, Gelovani JG, Kim SJ, Wang Z, Gallick GE. Targeting Src family kinases inhibits growth and lymph node metastases of prostate cancer in an orthotopic nude mouse model. Cancer Res. 2008;68:3323–33.

  61. 61.

    Musumeci F, Schenone S, Brullo C, Botta M. An update on dual Src/Abl inhibitors. Future Med Chem. 2012;4:799–822.

  62. 62.

    Breccia M, Salaroli A, Molica M, Alimena G. Systematic review of dasatinib in chronic myeloid leukemia. Onco Targets Ther. 2013;6:257–65.

  63. 63.

    Zheng BL, He K, Kim CH, Rogers LL, Shao Y, Huang ZY, Lu Y, Yan SJ, Qien LC, Zheng QY. Effect of a lipidic extract from lepidium meyenii on sexual behavior in mice and rats. Urology. 2000;55:598–602.

  64. 64.

    Hodgson SF, Watts NB, Bilezikian JP, Clarke BL, Gray TK, Harris DW, Johnston CC, Kleerekoper M, Lindsay R, Luckey MM, McClung MR, Nankin HR, Petak SM, Recker RR. American association of clinical endocrinologists medical guidelines for clinical practice for the prevention and treatment of postmenopausal osteoporosis: 2001 edition, with selected updates for 2003. Endocr Pract. 2003;9:544–64.

  65. 65.

    Zhang YZ, Yu LJ, Ao MZ, Jin WW. Effect of ethanol extract of Lepidium meyenii Walp. on osteoporosis in ovariectomized rat. J Ethnopharmacol. 2006;105:274–9.

  66. 66.

    Dahlman-Wright K, Cavailles V, Fuqua SA, Jordan VC, Katzenellenbogen JA, Korach KS, Maggi A, Muramatsu M, Parker MG, Gustafsson JA. International union of pharmacology LXIV. Estrogen receptors. Pharmacol Rev. 2006;58:773–81.

  67. 67.

    Eriksen EF, Colvard DS, Berg NJ, Graham ML, Mann KG, Spelsberg TC, Riggs BL. Evidence of estrogen receptors in normal human osteoblast-like cells. Science. 1988;241:84–6.

  68. 68.

    Sims NA, Dupont S, Krust A, Clement-Lacroix P, Minet D, Resche-Rigon M, Gaillard-Kelly M, Baron R. Deletion of estrogen receptors reveals a regulatory role for estrogen receptors-β in bone remodeling in females but not in males. Bone. 2002;30:18–25.

  69. 69.

    Deroo BJ, Korach KS. Estrogen receptors and human disease. J Clin Invest. 2006;116:561–70.

  70. 70.

    Clark CR, MacLusky NJ, Parsons B, Naftolin F. Effects of estrogen deprivation on brain estrogen and progestin receptor levels and the activation of female sexual behavior. Horm Behav. 1981;15:289–98.

  71. 71.

    Zhou S, Zilberman Y, Wassermann K, Bain SD, Sadovsky Y, Gazit D. Estrogen modulates estrogen receptor α and β expression, osteogenic activity, and apoptosis in mesenchymal stem cells (MSCs) of osteoporotic mice. J Cell Biochem Suppl. 2001;81:144–55.

  72. 72.

    Bakas P, Liapis A, Vlahopoulos S, Giner M, Logotheti S, Creatsas G, Meligova AK, Alexis MN, Zoumpourlis V. Estrogen receptor α and β in uterine fibroids: a basis for altered estrogen responsiveness. Fertil Steril. 2008;90:1878–85.

  73. 73.

    Shang Y, Brown M. Molecular determinants for the tissue specificity of SERMs. Science. 2002;295:2465–8.

  74. 74.

    Rochira V, Balestrieri A, Madeo B, Baraldi E, Faustini-Fustini M, Granata AR, Carani C. Congenital estrogen deficiency: in search of the estrogen role in human male reproduction. Mol Cell Endocrinol. 2001;178:107–15.

  75. 75.

    Cutillo DM. Method for metabolic correction. US Patent 20,150,313,949. 2015.

  76. 76.

    Oh H, Yang S, Park M, Chun JS. Matrix metalloproteinase (MMP)-12 regulates MMP-9 expression in interleukin-1β-treated articular chondrocytes. J Cell Biochem. 2008;105:1443–50.

  77. 77.

    Akhtar N, Miller MJ, Haqqi TM. Effect of a Herbal-Leucine mix on the IL-1β-induced cartilage degradation and inflammatory gene expression in human chondrocytes. BMC Complement Altern Med. 2011;11:66–76.

  78. 78.

    Delaissé JM, Andersen TL, Engsig MT, Henriksen K, Troen T, Blavier L. Matrix metalloproteinases (MMP) and cathepsin K contribute differently to osteoclastic activities. Microsc Res Tech. 2003;61:504–13.

  79. 79.

    Lee JC, Laydon JT, McDonnell PC, Gallagher TF, Kumar S, Green D, McNulty D, Blumenthal MJ, Heys JR, Landvatter SW. A protein kinase involved in the regulation of inflammatory cytokine biosynthesis. Nature. 1994;372:739–46.

  80. 80.

    Han J, Lee JD, Bibbs L, Ulevitch RJ. A MAP kinase targeted by endotoxin and hyperosmolarity in mammalian cells. Science. 1994;265:808–11.

  81. 81.

    Alexa A. Characterisation of mapk14/mapk11 mediated cardioprotection using RNAi. Circulation. 2013;128:A12734.

  82. 82.

    Waterworth DM, Li L, Scott R, Warren L, Gillson C, Aponte J, Sarov-Blat L, Sprecher D, Dupuis J, Reiner A, Psaty BM, Tracy RP, Lin H, McPherson R, Chissoe S, Wareham N, Ehm MG. A low-frequency variant in MAPK14 provides mechanistic evidence of a link with myeloperoxidase: a prognostic cardiovascular risk marker. J Am Heart Assoc. 2014;3:e001074.

  83. 83.

    Muslin AJ. MAPK signalling in cardiovascular health and disease: molecular mechanisms and therapeutic targets. Clin Sci. 2008;115:203–18.

  84. 84.

    Liu Z, Fan HY, Wang Y, Richards JS. Targeted disruption of Mapk14 (p38 MAPK α) 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:1794–804.

  85. 85.

    Wagner EF, Nebreda ÁR. Signal integration by JNK and p38 MAPK pathways in cancer development. Nat Rev Cancer. 2009;9:537–49.

  86. 86.

    Jopling C, Boue S, Izpisua Belmonte JC. Dedifferentiation, transdifferentiation and reprogramming: three routes to regeneration. Nat Rev Mol Cell Biol. 2011;12:79–89.

  87. 87.

    Knowlton DL, Tang K, Henstock PV, Subramanian RR. miRNA alterations modify kinase activation in the IGF-1 pathway and correlate with colorectal cancer stage and progression in patients. J Cancer. 2011;2:490–502.

  88. 88.

    Földes G, Mioulane M, Wright JS, Liu AQ, Novak P, Merkely B, Gorelik J, Schneider MD, Ali NN, Harding SE. Modulation of human embryonic stem cell-derived cardiomyocyte growth: a testbed for studying human cardiac hypertrophy? J Mol Cell Cardiol. 2011;50:367–76.

  89. 89.

    Michalik L, Auwerx J, Berger JP, Chatterjee VK, Glass CK, Gonzalez FJ, Grimaldi PA, Kadowaki T, Lazar MA, O’Rahilly S, Palmer CN, Plutzy J, Reddy JK, Spiegelman BM, Staels B, Wahli W. International Union of Pharmacology LXI. Peroxisome proliferator-activated receptors. Pharmacol Rev. 2006;58:726–41.

  90. 90.

    Dreyer C, Krey G, Keller H, Givel F, Helftenbein G, Wahli W. Control of the peroxisomal β-oxidation pathway by a novel family of nuclear hormone receptors. Cell. 1992;68:879–87.

  91. 91.

    Ziouzenkova O, Perrey S, Marx N, Bacqueville D, Plutzky J. Peroxisome proliferator-activated receptors. Curr Atheroscler Rep. 2002;4:59–64.

  92. 92.

    Abdelrahman M, Sivarajah A, Thiemermann C. Beneficial effects of PPAR-γ ligands in ischemia-reperfusion injury, inflammation and shock. Cardiovasc Res. 2005;65:772–81.

  93. 93.

    Nakamura K, Inoue I, Takahashi S, Komoda T, Katayama S. Cryptochrome and period proteins are regulated by the CLOCK/BMAL1 Gene: crosstalk between the PPARs/RXRα-regulated and CLOCK/BMAL1-regulated system. PPAR Res. 2008;2008:348610.

  94. 94.

    Chandra V, Huang P, Hamuro Y, Raghuram S, Wang Y, Burris TP, Rastinejad F. Structure of the intact PPAR-γ-RXR-α nuclear receptor complex on DNA. Nature. 2008;456:350–6.

  95. 95.

    Asakawa M, Takano H, Nagai T, Uozumi H, Hasegawa H, Kubota N, Saito T, Masuda Y, Kadowaki T, Komuro I. Peroxisome proliferator-activated receptor γ plays a critical role in inhibition of cardiac hypertrophy in vitro and in vivo. Circulation. 2002;105:1240–6.

  96. 96.

    Marx N, Sukhova GK, Collins T, Libby P, Plutzky J. PPARα activators inhibit cytokine-induced vascular cell adhesion molecule-1 expression in human endothelial cells. Circulation. 1999;99:3125–31.

  97. 97.

    Marx N, Bourcier T, Sukhova GK, Libby P, Plutzky J. PPARγ activation in human endothelial cells increases plasminogen activator inhibitor type-1 expression PPARγ as a potential mediator in vascular disease. Arterioscler Thromb Vasc Biol. 1999;19:546–51.

  98. 98.

    Kersten S, Desvergne B, Wahli W. Roles of PPARs in health and disease. Nature. 2000;405:421–4.

  99. 99.

    Tham DM, Martin-McNulty B, Wang YX, Wilson DW, Vergona R, Sullivan ME, Dole W, Rutledge JC. Angiotensin II is associated with activation of NF-κB-mediated genes and downregulation of PPARs. Physiol Genomics. 2002;11:21–30.

  100. 100.

    Rival Y, Benéteau N, Taillandier T, Pezet M, Dupont-Passelaigue E, Patoiseau JF, Junquéro D, Colpaert FC, Delhon A. PPARα and PPARδ activators inhibit cytokine-induced nuclear translocation of NF-κB and expression of VCAM-1 in EAhy926 endothelial cells. Eur J Pharmacol. 2002;435:143–51.

  101. 101.

    Henke BR. Peroxisome proliferator-activated receptor α/γ dual agonists for the treatment of type 2 diabetes. J Med Chem. 2004;47:4118–27.

Download references

Authors' contributions

HBL conceived and designed the study. FY, XLT and XY performed the experiments. FY analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We are very grateful to Dr. Jun Xu (Research Center for Drug Discovery, School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou) for his valuable suggestions. We would like to thank Mr. Ke-hao Wu (College of Life Science, Hunan Normal University, Changsha) for assistance with data processing.

Competing interests

The authors declare that they have no competing interests.

Author information

Correspondence to Hai-bo Liu.

Additional files

13020_2016_112_MOESM1_ESM.xlsx

Additional file 1. IPA druggable datasheet.

13020_2016_112_MOESM2_ESM.txt

Additional file 2. WEGA validation scores.

13020_2016_112_MOESM3_ESM.pdf

Additional file 3. Total GeneMANIA report of maca.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Prostate Cancer
  • Osteoporosis
  • Androgen Receptor
  • Dasatinib
  • Virtual Screening