A molecular docking study of phytochemical estrogen mimics from dietary herbal supplements

In Silico Pharmacol. 2015; 3: 4. Published online 2015 Mar 22. doi: 10.1186/s40203-015-0008-z PMCID: PMC4397262 Chelsea N Powers and William N Setzercorresponding author Department of Chemistry, University of Alabama in Huntsville, Huntsville, AL 35899 USA Chelsea N Powers, Email: ude.hau@7000pnc. Contributor Information. corresponding authorCorresponding author. This article has been cited by other articles in PMC. Go to: Abstract Purpose The purpose of this study is to use a molecular docking approach to identify potential estrogen mimics or anti-estrogens in phytochemicals found in popular dietary herbal supplements. Methods In this study, 568 phytochemicals found in 17 of the most popular herbal supplements sold in the United States were built and docked with two isoforms of the estrogen receptor, ERα and ERβ (a total of 27 different protein crystal structures). Results The docking results revealed six strongly docking compounds in Echinacea, three from milk thistle (Silybum marianum), three from Gingko biloba, one from Sambucus nigra, none from maca (Lepidium meyenii), five from chaste tree (Vitex agnus-castus), two from fenugreek (Trigonella foenum-graecum), and two from Rhodiola rosea. Notably, of the most popular herbal supplements for women, there were numerous compounds that docked strongly with the estrogen receptor: Licorice (Glycyrrhiza glabra) had a total of 26 compounds strongly docking to the estrogen receptor, 15 with wild yam (Dioscorea villosa), 11 from black cohosh (Actaea racemosa), eight from muira puama (Ptychopetalum olacoides or P. uncinatum), eight from red clover (Trifolium pratense), three from damiana (Turnera aphrodisiaca or T. diffusa), and three from dong quai (Angelica sinensis). Of possible concern were the compounds from men’s herbal supplements that exhibited strong docking to the estrogen receptor: Gingko biloba had three compounds, gotu kola (Centella asiatica) had two, muira puama (Ptychopetalum olacoides or P. uncinatum) had eight, and Tribulus terrestris had six compounds. Conclusions This molecular docking study has revealed that almost all popular herbal supplements contain phytochemical components that may bind to the human estrogen receptor and exhibit selective estrogen receptor modulation. As such, these herbal supplements may cause unwanted side effects related to estrogenic activity. Keywords: Molecular docking, Estrogen receptor, Herbal supplements Go to: Background The use of alternative medicines in the United States, particularly herbal supplements, has dramatically increased since the beginning of the 21st century (Figure 1). Filling American minds with promises of enhanced beauty, sharper senses, and optimum organ functions, herbal supplements claim to increase, or improve almost all issues a person could have with their body. Without a doubt it is appealing to have problems solved by simply swallowing a pill or drinking a tea, not much effort required, however it has been widely ignored the consecutive consequences these supplements can provide (Cupp 1999). Figure 1 Figure 1 Relationship between herbal supplement purchases in the United States and the year. Two major factors play a part in the ongoing, unnoticed herbal supplement crisis: Regulations for herbal supplements and uneducated consumers. Beginning with the first, the United States does not classify herbal supplements as drugs, and therefore supplements are not required to undergo the extensive testing that pharmaceutical drugs do before put on the market. Courtesy of the “Dietary Supplement Health and Education Act of 1994”, herbal supplements are not evaluated by the Food and Drug Administration (Calixto 2000) making it easy for supplement companies to rapidly introduce new supplements to consumers, with or without the knowledge of possible harmful side effects. Unspecified drugs, contaminations, toxins, and/or heavy metals (Au et al. 2000) can be included in an herbal supplement, and since companies are not required to subject their products to quality analysis, this spectrum of harmful compounds could be digested by a consumer and induce adverse effects. As for the second, biologically uneducated consumers do not understand or simply do not consider the concept that plants are not always beneficial. They believe anything that is natural must be good for their health and safe to consume (Stonemetz 2008), which is far from the truth. Plants contain hundreds of phytochemicals, some of which are indeed toxic to the human body. One class of phytochemicals of major concern, which is the focus of this study, phytoestrogens, can interfere and react with the human estrogen receptors, which regulate neural, skeletal, cardiovascular, and reproductive tissues. This interference, however, is not always adverse. For example, some phytoestrogens can promote carcinogenic growth, while others can inhibit the growth. The purpose of this study was to identify potential estrogen mimics or anti-estrogens in phytochemicals found in popular dietary herbal supplements. The data gathered can only suggest the possibility of a phytochemical to be an anti-estrogen or a mimic, not confirm its estrogenic properties. It is our hope that the discoveries made during this study can help to identify the estrogenic activity of the phytochemicals examined. This information can then lead to the health benefits or hazards associated with the phytochemicals, which in turn could greatly affect the increasingly popular herbal supplement movement. Go to: Methods Literature survey A literature survey on herbal supplements was carried out to identify the most popular general [Echinacea, milk thistle (Silybum marianum), Ginkgo biloba, Sambucus nigra, maca (Lepidium meyenii), chaste tree (Vitex agnus-castus), fenugreek (Trigonella foenum-graecum), and Rhodiola rosea], women’s [damiana leaf (Turnera aphrodisiaca, T. diffusa), muira puama (Ptychopetalum olacoides, P. uncinatum), black cohosh (Actaea racemosa = Cimifuga racemosa), licorice root (Glycyrrhiza glabra), wild yam (Dioscorea villosa), dong quai (Angelica sinensis) and red clover (Trifolium pretense)], and men’s [Gingko biloba, gotu kola (Centella asiatica), muira puama (Ptychopetalum olacoides, P. uncinatum), and Tribulus terrestris] herbal supplements advertised and used in the United States. A survey of the literature, including the Dictionary of Natural Products (2014) and Duke’s Phytochemical Database (1998), was carried out to determine the phytochemical constituents of each herb. Molecular modeling of phytochemicals Each phytochemical ligand structure (see Figures 2, ​,3,3, ​,4,4, ​,5,5, ​,6,6, ​,7,7, ​,8,8, ​,9,9, ​,10,10, ​,11,11, ​,12,12, ​,13,13, ​,14,14, ​,15,15, ​,16,16, ​,17,17, ​,18,18, ​,19,19, ​,20,20, ​,21,21, ​,22,22, ​,23,23, ​,24,24, ​,25,25, ​,26,26, ​,27,27, ​,28,28, ​,29,29, ​,30,30, ​,31,31, ​,32,32, ​,33,33, ​,34,34, ​,35,35, ​,36,36, ​,37,37, ​,38,38, ​,39,39, ​,40,40, ​,41,41, ​,42,42, and ​and43)43) was built using Spartan ’14 for Windows (2013). For each ligand, a conformational search and geometry optimization was carried out using the MMFF force field (Halgren 1996). Figure 2 Figure 2 Alkaloid ligands examined in this work. Figure 3 Figure 3 Chalcone ligands examined in this work. Figure 4 Figure 4 Coumarin ligands examined in this work. Figure 5 Figure 5 Additional coumarin ligands examined in this work. Figure 6 Figure 6 Diterpenoid ligands examined in this work. Figure 7 Figure 7 Additional diterpenoid ligands examined in this work. Figure 8 Figure 8 Additional diterpenoid ligands examined in this work. Figure 9 Figure 9 Additional diterpenoid ligands examined in this work. Figure 10 Figure 10 Flavonoid ligands examined in this work. Figure 11 Figure 11 Additional flavonoid ligands examined in this work. Figure 12 Figure 12 Additional flavonoid ligands examined in this work. Figure 13 Figure 13 Additional flavonoid ligands examined in this work. Figure 14 Figure 14 Additional flavonoid ligands examined in this work. Figure 15 Figure 15 Additional flavonoid ligands examined in this work. Figure 16 Figure 16 Additional flavonoid ligands examined in this work. Figure 17 Figure 17 Additional flavonoid ligands examined in this work. Figure 18 Figure 18 Isoflavonoid ligands examined in this work. Figure 19 Figure 19 Additional isoflavonoid ligands examined in this work. Figure 20 Figure 20 Additional isoflavonoid ligands examined in this work. Figure 21 Figure 21 Lignan ligands examined in this work. Figure 22 Figure 22 Phenanthrenoid ligands examined in this work. Figure 23 Figure 23 Miscellaneous phenolic compounds examined in this work. Figure 24 Figure 24 Additional miscellaneous phenolic compounds examined in this work. Figure 25 Figure 25 Additional miscellaneous phenolic compounds examined in this work. Figure 26 Figure 26 Sesquiterpenoid ligands examined in this work. Figure 27 Figure 27 Steroid ligands examined in this work. Figure 28 Figure 28 Additional steroid ligands examined in this work. Figure 29 Figure 29 Additional steroid ligands examined in this work. Figure 30 Figure 30 Additional steroid ligands examined in this work. Figure 31 Figure 31 Additional steroid ligands examined in this work. Figure 32 Figure 32 Additional steroid ligands examined in this work. Figure 33 Figure 33 Additional steroid ligands examined in this work. Figure 34 Figure 34 Additional steroid ligands examined in this work. Figure 35 Figure 35 Stilbenoid ligands examined in this work. Figure 36 Figure 36 Additional stilbenoid ligands examined in this work. Figure 37 Figure 37 Triterpenoid ligands examined in this work. Figure 38 Figure 38 Additional triterpenoid ligands examined in this work. Figure 39 Figure 39 Additional triterpenoid ligands examined in this work. Figure 40 Figure 40 Additional triterpenoid ligands examined in this work. Figure 41 Figure 41 Additional triterpenoid ligands examined in this work. Figure 42 Figure 42 Miscellaneous phytochemical ligands examined in this work. Figure 43 Figure 43 Additional miscellaneous phytochemical ligands examined in this work. Molecular docking Protein-ligand docking studies were carried out based on the crystal structures of human estrogen receptor α [ERα: PDB 1X7E (Manas et al. 2004a), PDB 1X7R (Manas et al. 2004b), and PDB 3ERD (Shiau et al. 1998)] and human estrogen receptor β [ERβ: PDB 1U3Q, 1U3R, 1U3S (Malamas et al. 2004), 1U9E, 1X7B, 1X76, 1X78 (Manas et al. 2004a), and 1X7J (Manas et al. 2004b)]. Prior to docking all solvent molecules and the co-crystallized ligands were removed from the structures. Molecular docking calculations for all compounds with each of the proteins were undertaken using Molegro Virtual Docker v. 6.0 (2013). Potential binding sites in the protein structures were identified using the grid-based cavity prediction algorithm of the Molegro Virtual Docker (2013) program. The location of the volume used by the docking search algorithm was positioned at the center of the cavity and a sphere (15 Å radius) large enough to encompass the entire cavity of the binding site of each protein structure was selected in order to allow each ligand to search. If a co-crystallized inhibitor or substrate was present in the structure, then that site was chosen as the binding site. If no co-crystallized ligand was present, then suitably sized (>50 Å3) cavities were used as potential binding sites. The docking searches were constrained to those cavities. Standard protonation states of the proteins based on neutral pH were used in the docking studies. Each protein was used as a rigid model structure; no relaxation of the protein was performed. Assignments of charges on each protein were based on standard templates as part of the Molegro Virtual Docker (2013) program (Thomsen and Christensen 2006); no other charges were necessary to be set. Flexible ligand models were used in the docking and subsequent optimization scheme. As a test of docking accuracy and for docking energy comparison, co-crystallized ligands were re-docked into the protein structures (see Table 1). Additionally, as positive controls, the known estrogenic compounds 17β-estradiol and α-zearalenone were docked with each protein structure in order to compare docking energies with the herbal phytochemicals. Different orientations of the ligands were searched and ranked based on their energy scores. The RMSD threshold for multiple cluster poses was set at <1.00 Å. The docking algorithm was set at maximum iterations of 1500 with a simplex evolution population size of 50 and a minimum of 30 runs for each ligand. Each binding site of oligomeric structures was searched with each ligand. The lowest-energy (strongest-docking) poses for each ligand in each protein target are summarized in Tables 2, ​,3,3, ​,4,4, ​,5,5, ​,6,6, ​,7,7, ​,8,8, ​,9,9, ​,10,10, ​,11,11, ​,12,12, ​,13,13, ​,1414 and ​and1515. Table 1 Table 1 MolDock docking energies of co-crystallized ligands and root-mean-squared deviations between the co-crystallized ligand and the re-docked poses of the co-crystallized ligand with human estrogen receptors α and β Table 2 Table 2 MolDock molecular docking energies (kJ/mol) for alkaloids with human estrogen receptors α and β Table 3 Table 3 MolDock molecular docking energies (kJ/mol) for chalcones with human estrogen receptors α and β Table 4 Table 4 MolDock molecular docking energies (kJ/mol) for coumarins with human estrogen receptors α and β Table 5 Table 5 MolDock molecular docking energies (kJ/mol) for diterpenoids with human estrogen receptors α and β Table 6 Table 6 MolDock molecular docking energies (kJ/mol) for flavonoids with human estrogen receptors α and β Table 7 Table 7 MolDock molecular docking energies (kJ/mol) for isoflavonoids with human estrogen receptors α and β Table 8 Table 8 MolDock molecular docking energies (kJ/mol) for lignans with human estrogen receptors α and β Table 9 Table 9 MolDock molecular docking energies (kJ/mol) for phenanthrenoids with human estrogen receptors α and β Table 10 Table 10 MolDock molecular docking energies (kJ/mol) for miscellaneous phenolic ligands with human estrogen receptors α and β Table 11 Table 11 MolDock molecular docking energies (kJ/mol) for sesquiterpenoids with human estrogen receptors α and β Table 12 Table 12 MolDock molecular docking energies (kJ/mol) for steroids with human estrogen receptors α and β Table 13 Table 13 MolDock molecular docking energies (kJ/mol) for stilbenoids with human estrogen receptors α and β Table 14 Table 14 MolDock molecular docking energies (kJ/mol) for triterpenoids with human estrogen receptors α and β Table 15 Table 15 MolDock molecular docking energies (kJ/mol) for miscellaneous phytochemicals with human estrogen receptors α and β Go to: Results Alkaloids The alkaloid ligands examined in this study are shown in Figure 2. The molecular docking results for the alkaloids are summarized in Table 2. Of the alkaloids examined in this study, cis- and trans-clovamide, with docking energies of −119.8 and −113.6 kJ/mol, respectively, and N-trans-feruloyltyramine (Edock = −103.1 kJ/mol) were found to dock well with ERα. Their docking energies were more exothermic than those of estradiol, −92.0 kJ/mol and the corresponding co-crystallized ligand genistein, −93.4 kJ/mol, and the clovamides were more exothermic than zearalenone (Edock = −104.1 kJ/mol). The co-crystallized ligand, genistein, and the clovamide and feruloylyramine ligands have similar positions in the binding site (Figure 44). Phe 404, Leu 525, Leu 346, Leu 387, and Leu 391 form a hydrophobic pocket around the docked alkaloids. Phe 404 exhibited edge-to-face π–π interactions between the phenyl substituent of Phe with the caffeic or ferulic substituents of the alkaloids and with the hydroxyphenyl substituent of genistein. Notable hydrogen bonds in the lowest-energy docked pose of cis-clovamide were the 3-OH and 4-OH of the cis-caffeic moiety with the carboxylate residue of Glu 353 and the 3-OH group with the guanidine residue of Arg 394 (Figure 45). The docked trans-clovamide had hydrogen bonds between the 4-OH of the caffeate with the guanidine of Arg 394 and the carbonyl group of Leu 387 and the 3-OH group with the carboxylate of Glu 353. Hydrogen bonds were formed between the 4-OH group on the ferulyl substituent of N-trans-feruloyltyramine and carbonyl group of Leu 387, the guanidine group of Arg 394, and the carboxylate of Glu 353. Figure 44 Figure 44 Lowest-energy docked poses of alkaloids [ N-trans -feruloyltyramine (aqua), cis- clovamide (red), and trans- clovamide (blue) along with the co-crystallized ligand, genistein (green)] with ERα (PDB 1X7R). A: Docked poses showing the entire ribbon ... Figure 45 Figure 45 Lowest-energy docked pose of cis- clovamide with ERα (PDB 1X7R) showing the principle amino acid contacts in the binding site. Hydrogen bonds are indicated by blue dashed lines. Similarly, cis-clovamide, trans-clovamide, and N-trans-feruloyltyramine were the alkaloids that docked well with ERβ. Their docking energies (−124.9, −122.0, and −113.8 kJ/mol, respectively) were more exothermic than those of estradiol, −100.0 kJ/mol, zearalenone, −104.9 kJ/mol, and the corresponding co-crystallized ligand 2-(3-fluoro-4-hydroxyphenyl)-7-vinyl-1,3-benzoxazol-5-ol, −107.9 kJ/mol. The alkaloid and the co-crystallized ligand occupied similar positions in the binding site, a hydrophobic pocket formed by Leu 298, Phe 356, Leu 339, and His 475. Phe 356 exhibited edge-to-face π–π interactions with the caffeic or ferulic substituents of the docked alkaloid ligands as well as with the hydroxyphenyl substituent of the co-crystallized ligand. There were two notable hydrogen bonds formed between the 4-OH group on the ferulyl substituent of N-trans-feruloyltyramine and the guanidine group of Arg 346, and the carboxylate of Glu 305 (Figure 46). These same two residues formed hydrogen bonds with the 4-hydroxyphenyl group of the co-crystallized ligand. The caffeoyl group of cis-clovamide formed hydrogen bonds with Glu 305 and Leu 298. trans-Clovamide, however, formed hydrogen bonds with Leu 339, Arg 346, and Glu 305. Figure 46 Figure 46 Lowest-energy docked pose of N-trans- feruloyltyramine with ERβ (PDB 1X7B) showing the principle amino acid contacts in the binding site. Hydrogen bonds are indicated by blue dashed lines. Chalcones The structures of the chalcones are shown in Figure 3, while the docking energies are summarized in Table 3. Xanthohumol was the strongest docking chalcone with ERα. Its docking energy, −116.8 kJ/mol, is more exothermic than those of estradiol, −92.0 kJ/mol, zearalenone, −104.1 kJ/mol, and the corresponding co-crystallized ligand genistein, −93.4 kJ/mol. The 4-OH group of xanthohumol forms three hydrogen-bonds with the protein (the carboxylate of Glu 353, the guanidine of Arg 394, and the carbonyl oxygen of Phe 404). The 4′-OH group of xanthohumol forms hydrogen-bonds with the imidazole N-H of His 524 and the carbonyl oxygen of Gly 521. Kanzonol Y (Edock = −111.2 kJ/mol) and licochalcone B (Edock = −107.8 kJ/mol) were the only other chalcone ligands to dock well with ERα. Of the chalcone ligands examined, kanzonol Y (Edock = −122.4 kJ/mol), xanthohumol (Edock = −116.8 kJ/mol), and licoagrochalcone A (Edock = −115.5 kJ/mol), docked best with ERβ. Their docking energies were decidedly more exothermic than those of estradiol, zearalenone, and the corresponding co-crystallized ligand 2-(5-hydroxy-naphthalen-1-yl)-1,3-benzooxazol-6-ol (Edock = −109.2 kJ/mol). Apparently, the hydrophobic prenyl groups allow for stronger docking. Thus, kanzonol Y docked to ERβ much better than the non-prenylated α,2′,4,4′-tetrahydroxydihydrochalcone (Edock = −105.0 kJ/mol). In the lowest-energy docked pose of kanzonol Y, the 3-prenyl group is sandwiched between the hydrophobic residues of Phe 356 and Leu 339, while the 5′-prenyl group is sandwiched between Leu 476 and Thr 299 (Figure 47). It has been shown that prenylation of flavonoids and related compounds does alter the estrogenic activity and often results in antiestrogenic activity (Kretzschmar et al. 2010; Simons et al. 2012). Figure 47 Figure 47 Lowest-energy docked pose of the prenylated chalcone kanzonol Y with ERβ (PDB 1X7J) showing the hydrophobic amino acid contacts with the isoprenyl groups of the ligand. Coumarins The MolDock docking energies of the coumarins are summarized in Table 4, and the structures of the coumarins are shown in Figures 4 and ​and5.5. Glabrene and pratenol B were the strongest docking coumarins with ERα (Edock = −104.8 kJ/mol), more exothermic than either estradiol (−92.0 kJ/mol) or genistein (−93.4 kJ/mol), and comparable to zearalenone (−104.1 kJ/mol). The 7-hydroxychromene moieties of both glabrene and pratenol B are held in a hydrophobic pocket formed by Leu 387, Phe 404, Met 388, and Leu 391. Furthermore, there are edge-to-face π–π interactions between Phe 404 and the chromene benzene rings of the ligands, as well as hydrogen bonds between the 7-OH group of the chromene and the guanidine group of Arg 394 and the carboxylate of Glu 356, analogous to the co-crystallized ligand genistein. Additionally, one of the carboxylates of pratenol B forms a hydrogen bond with imidazole substituent NH group of His 475. In addition to glabrene (Edock = −114.9 kJ/mol) and pratenol B (Edock = −112.0 kJ/mol), mirificoumestan docked strongly with ERβ with a docking energy of −113.0 kJ/mol. These compounds docked more exothermically than estradiol, zearalenone, or the co-crystallized ligand [5-hydroxy-2-(4-hydroxyphenyl)-1-benzofuran-7-yl]acetonitrile (Edock = −107.7 kJ/mol). Glabrene docked into the hydrophobic pocket formed by Leu 298, Leu 339, and Phe 356. Phe 356 exhibited edge-to-face π–π interactions with the hydroxychromene substituent of glabrene. There was a hydrogen bond between the imidazole substituent NH group of His 475 and the 5′-OH group of glabrene, and three hydrogen bonding interactions were seen between the 7-OH group of glabrene and the carbonyl group of Leu 339, the guanidine moiety of Arg 346, and the carboxylate of Glu 305. Diterpenoids The structures for the diterpenoid ligands examined in this work are shown in Figures 6, ​,7,7, ​,8,8, and ​and9,9, and the docking energies are listed in Table 5. The strongest docking diterpenoids with ERα were diosbulbin F, diosbulbin K, and diosbulbin L (Edock = −111.2, −112.1, and −110.8 kJ/mol, respectively). Each of these ligands docked in the hydrophobic binding pocket with the furan group forming a hydrogen bond to the guanidine of Arg 394. The methyl ester of diosbulbin F and the carboxylate of diosbulbin L also formed hydrogen bonds to the imizadole of His 524. These three diterpenoids also docked strongly to ERβ. In addition, diosbulbin F, diosbulbin H, ptycho-6α,7α-diol, and ptycholide IV all had docking energies more exothermic (Edock = −114.8, −114.1, −122.9, and −114.7 kJ/mol) than the co-crystallized ligand, 2-(5-hydroxynaphthalen-1-yl)-1,3-benzooxazol-6-ol (Edock = −111.3 kJ/mol). Ptycho-6α,7α-diol docked with ERβ with hydrogen bonds between the lactone carbonyl of the ligand and Arg 346 and the 6-OH group of the ligand with His 475. Flavonoids The structures of the flavonoid ligands examined in this work are shown in Figures 10, ​,11,11, ​,12,12, ​,13,13, ​,14,14, ​,15,15, ​,16,16, and ​and17,17, and the docking energies are listed in Table 6. Of the flavonoids, luteolin-8-propenoic acid docked the strongest to ERα, with a docking energy of −113.127 kJ/mol, more exothermic than those of estradiol, zearalenone, or genistein. A common docking orientation for phenolic ligands in ERα is the hydrophobic pocket of Leu 387, Phe 404, Met 388, and Leu 391, along with edge-to-face π–π interactions with Phe 404, and hydrogen bonds between the phenolic –OH group and the guanidine group of Arg 394 and the carboxylate of Glu 356. The 7-OH group of the ligand made an additional hydrogen bond with the carbonyl oxygen of Gly 521. No other flavonoid ligands showed notably strong docking with ERα. Luteolin-8-propenoic acid was also the strongest docking flavonoid with ERβ (Edock = −123.1 kJ/mol), is far more exothermic than estradiol, zearalenone, and the co-crystallized ligand, 2-(3-fluoro-4-hydroxyphenyl)-7-vinyl-1,3-benzoxazol-5-ol (Edock = −106.2 kJ/mol). As observed in other phenolic ligands with ERβ, luteolin-8-propenoic acid occupied the hydrophobic pocket formed by Leu 298, Leu 339, and Phe 356; edge-to-face π–π interactions of the phenolic ligand with Phe 356 and hydrogen boding of the phenolic –OH group with the carbonyl group of Leu 339, the guanidine group of Arg 346, and the carboxylate of Glu 305. The 7-OH group of the ligand made additional hydrogen bonds with the carbonyl oxygen of Gly 472 and His 475. Casticin (Edock = −106.4 kJ/mol), gonzalitosin (Edock = −106.3 kJ/mol), gossypetin (Edock = −105.7 kJ/mol), laricitrin (Edock = −107.3 kJ/mol), myricetin (Edock = −106.2 kJ/mol), quercetin (Edock = −106.0 kJ/mol), santin (Edock = −108.0 kJ/mol), and tricetin (Edock = −106.2 kJ/mol), all had more exothermic docking energies than estradiol or zearalenone but were less exothermic than the co-crystallized ligand, 2-(5-hydroxynaphthalen-1-yl)-1,3-benzooxazol-6-ol (Edock = −111.3 kJ/mol). Isoflavonoids The docking energies of the isoflavonoids are summarized in Table 7 and the structures are shown in Figure 18, ​,19,19, and ​and20.20. Genistein is the quintessential estrogenic isoflavonoid, but it is a weaker docking ligand than estradiol or zearalenone for either ERα or ERβ. The strongest docking isoflavonoid with ERα was the aglycone of licoagroside A (Edock = −100.5 kJ/mol), but this was weaker than zearalenone. On the other hand, several isoflavonoid ligands docked to ERβ more strongly than zearalenone: licoagroside A aglycone (Edock = −107.7 kJ/mol), 1-methoxyphaseollin (Edock = −110.5 kJ/mol), pratensein (Edock = −106.4 kJ/mol), and 3′,5,7-trihydroxy-5′-methoxyisoflavone (Edock = −107.9 kJ/mol), but none of these docked more strongly than the synthetic co-crystallized ligand, 2-(5-hydroxynaphthalen-1-yl)-1,3-benzooxazol-6-ol (Edock = −111.3 kJ/mol). Lignans The structures and the docking energies of the lignans are shown in Figure 21 and Table 8, respectively. Nortrachelogenin and 7′-hydroxymatairesinol were the strongest docking lignans to ERα (Edock = −112.0 and −112.3 kJ/mol, respectively). Nortrachelogenin was also the strongest docking lignan to ERβ (Edock = −125.4 kJ/mol). Sesamin showed notable selectivity for ERβ over ERα (Edock = −121.8 and −99.1 kJ/mol, respectively). Nortrachelogenin occupied the same orientation and hydrogen-bonding pattern in both ERα and ERβ, with one of the phenolic –OH groups hydrogen bonded to argenine in the binding pocket (Arg 394 in ERα; Arg 346 in ERβ) and the other phenolic –OH group hydrogen bonded to the histidine (His 524 in ERα; His 475 in ERβ). Phenanthrenoids The docking energies and structures of phenanthrenoids are shown in Table 9 and Figure 22, respectively. None of the phenanthrenoids examined in this work showed docking energies lower than estradiol or zearalenone for ERα or ERβ. Miscellaneous phenolics The docking energies for miscellaneous herbal phenolic compounds are listed in Table 10, and the structures are shown in Figures 23, ​,24,24, and ​and25.25. The strongest docking ligands of the miscellaneous phenolic compounds for ERα was cimicifugic acid F (Edock = −126.2 kJ/mol), and this ligand also docked strongly with ERβ (Edock = −125.2 kJ/mol). Several other phenolic ligands docked with very exothermic energies to ERβ: the aglycone of agnucastoside C (Edock = −130.0 kJ/mol), caffeoyl-p-coumaroyl tartaric acid (Edock = −129.8 kJ/mol), cimiracemate B (Edock = −127.3 kJ/mol), cimiracemate D (Edock = −128.5 kJ/mol), and fukinolic acid (Edock = −127.3 kJ/mol). Analogous to other phenolic compounds (see above), the cinnamate moiety of the lowest-energy pose of cimicifugic acid F in ERα shows edge-to-face π–π interactions with Phe 404, and hydrogen bonding between the –OCH3 group and the guanidine group of Arg 394. Additionally, the 3-carboxylate group of the ligand is hydrogen-bonded to His 524, and the phenolic group fits into a hydrophobic pocket formed by Met 388, Met 421, and Ile 424. The lowest-energy docked pose of the aglycone of agnucastoside C with ERβ, as with other phenolic compounds (see above), has the p-coumarate phenolic –OH group hydrogen bonded to the carbonyl group of Leu 339 and the guanidine group of Arg 346, and edge-to-face π–π interactions of the phenolic ligand with Phe 356. The cyclic hemiacetal group is hydrogen bonded to His 475. Sesquiterpenoids Of the sesquiterpenoids, only cinnamoylechinadiol gave a notable docking energy (−120.8 kJ/mol) with ERβ (Figure 26, Table 11). Steroids The steroidal ligands examined in this study are shown in Figures 27, ​,28,28, ​,29,29, ​,30,30, ​,31,31, ​,32,32, ​,33,33, and ​and34,34, and their docking energies are listed in Table 12. Several pregnane steroids exhibited docking energies less than the co-crystallized ligand 2-(5-hydroxynaphthalen-1-yl)-1,3-benzooxazol-6-ol (Edock = −111.3 kJ/mol) for ERβ: 2,3-Dihydroxypregn-16-en-20-one (Edock = −116.6 kJ/mol), 3,16-dihydroxypregn-5-en-20-one (Edock = −116.6 kJ/mol), 3,21-dihydroxypregna-5,16-dien-20-one (Edock = −121.4 kJ/mol), and pregnadienolone (Edock = −115.4 kJ/mol). The lowest-energy docked poses of the pregnane ligands show them all to adopt the same orientation (Figure 48) with key hydrogen-bonding interactions of the 3-OH group of the steroids with the guanidine of Arg 346 and the amide carbonyl of Leu 339, and the 20-ketone group of the ligand with the imidazole N-H of His 475. Two ligands, 3,16-dihydroxypregn-5-en-20-one and 3,21-dihydroxypregna-5,16-dien-20-one, showed selectivity for ERβ over ERα (24.8 and 18.4 kJ/mol, respectively). Figure 48 Figure 48 Lowest-energy docked poses of pregnane steroids [2,3-dihydroxypregn-16-en-20-one (magenta), 3,16-dihydroxypregn-5-en-20-one (dark green), 3,21-dihydroxypregna-5,16-dien-20-one (white), and pregnadienolone (bright green)] with ERβ (PDB 1U3R). ... Stilbenoids Structures and docking energies for the stilbenoid ligands are shown in Figures 35 and ​and36,36, and Table 13, respectively. Several stilbenoid ligands showed notably strong docking energies; lower than estradiol, zearalenone, or the respective co-crystallized ligands: 3-acetoxy-4′,5-dihydroxy-3′-prenyldihydrostilbene (Edock to ERα = −119.0 kJ/mol), licoagrodione (Edock to ERβ = −116.9 kJ/mol), 3,3′,4,5′-tetrahydroxy-4′,5-diprenylbibenzyl (Edock to ERβ = −117.1 kJ/mol), 3,3′,4,5′-tetrahydroxy-5-prenylbibenzyl (Edock to ERβ = −115.0 kJ/mol), and uralstilbene (Edock to ERβ = −122.1 kJ/mol). Prenylation of stilbenoids seems to improve docking energies by about 20 kJ/mol. Triterpenoids Docking energies are presented in Table 14 and structures of triterpenoid ligands are illustrated in Figures 37, ​,38,38, ​,39,39, ​,40,40, and ​and41.41. Unlike the pregnane steroids, there were no triterpenoid ligands that showed good docking with either ERα or ERβ. Miscellaneous phytochemicals Several miscellaneous phytochemicals found in herbal supplements were included in this study (Table 15, Figures 42 and ​and43).43). Of these ligands, orobanchyl acetate gave excellent docking energies for both ERα and ERβ (Edock = −111.3 and −122.8 kJ/mol, respectively). Go to: Discussion Angelica sinensis Dong quai (Angelica sinensis root) has been used in Chinese traditional medicine for thousands of years for various female health conditions (e.g., dysmenorrhea, pelvic pain, symptoms of menopause) (Chye 2006; Al-Bareeq et al. 2010; Fang et al. 2012). In spite of its history, dong quai provided no clinical relief of menopausal symptoms (Hirata et al. 1997). In fact, dong quai has been shown to stimulate the growth of MCF-7 (ER+ human mammary carcinoma) cells (Lau et al. 2005), but does not bind either ERα or ERβ (Liu et al. 2001). The plant contains several miscellaneous phytochemicals, only two of which have notable docking energies, 10-angeloylbutylphthalide (−107.1 kJ/mol with ERβ) and angeliferulate (−110.7 and −121.5 kJ/mol with ERα and ERβ, respectively). Centella asiatica Centella asiatica (gotu kola) has been used in Ayurvedic traditional medicine for cognitive enhancement (Rao et al. 2005), to alleviate symptoms of anxiety and promote relaxation (Wijeweera et al. 2006), as well as for headache, body ache, asthma, ulcers, and wound healing (Kumar and Gupta 2002). Animal studies have revealed cognitive enhancement (Kumar and Gupta 2002; Rao et al. 2005), neuroprotective (Subathra et al. 2005), and anxiolytic (Wijeweera et al. 2006) effects. To our knowledge, there have been no reports on the estrogenic activity of C. asiatica. The plant contains the flavonoids castillicetin, castilliferol, kaempferol, and quercetin; the triterpenoids 2,3,20,23-tetrahydroxy-28-ursanoic acid, 2,3,23-trihydroxy-20-ursen-28-oic acid, 2,3-dihydroxy-5-(hydroxymethyl)-24-norolean-12-en-28-oic acid, 3,6,23-trihydroxy-12-ursen-28-oic acid, 6β-hydroxymaslinic acid, asiatic acid, asiaticoside G, betulafolienetriol, centellasapogenol A, centelloside A, corosolic acid, isothankunic acid, madasiatic acid, madecassic acid, quasipanaxadiol, terminolic acid, uncargenin C and zemoside A; the steroid campesterol, and the miscellaneous compounds asiaticin, homosilphiperfoloic acid, and irbic acid. Both quercetin and asiaticin had notable exothermic docking energies with ERβ (−106.0 and −109.0 kJ/mol, respectively). Quercetin has shown preferential binding to ERβ (Kuiper et al. 1998). Cimicifuga racemosa (syn. Actaea racemosa) Although black cohosh extracts have demonstrated clinical efficacy against some symptoms of menopause (Lieberman 1998; McKenna et al. 2001; Liske et al. 2002; Pockaj et al. 2004; Wuttke et al. 2006), several studies have demonstrated little or no estrogenic activity (Liu et al. 2001; Kronenberg 2003; Lupu et al. 2003; Mahady 2003). The efficacy of C. racemosa extracts on post-menopausal symptoms has been attributed to partial agonism of the serotonin receptor (Burdette et al. 2003) and the μ-opiate receptor (Rhyu et al. 2006). C. racemosa extracts have revealed several triterpenoids (Shao et al. 2000), phenylpropanoids (Chen et al. 2002) and caffeic acid derivatives (Li et al. 2003). Very few of the C. racemosa triterpenoids showed negative docking energies and are, therefore, unlikely estrogen receptor binding agents. Several C. racemosa phenolic compounds did show remarkable docking affinities for both ERα and ERβ: cimicifugic acid A, cimicifugic acid B, cimicifugic acid G, cimiciphenol, cimiciphenone, cimiracemate A, cimiracemate B, cimiracemate C, cimiracemate D, and fukinolic acid. It is likely that any estrogenic activity of C. racemosa extract (Seidlová-Wuttke et al. 2003a) is due to the presence of phenolic components rather than triterpenoids. Dioscorea villosa The rhizomes of wild yam, Dioscorea villosa, have been used to treat symptoms associated with menopause and premenstrual syndrome (PMS) as well as to relieve labor pains and sooth dysmenorrhea (Dutta 2015). The genus contains numerous steroidal glycosides (Sautour et al. 2006; Sautour et al. 2007; Ali et al. 2013). In this work, we have carried out in-silico screening of phytochemicals from the genus Dioscorea (Dictionary of Natural Products 2014). Of these, the diosbulbins D, F, H, J, K, and L (diterpenoids from D. bulbifera (Komori 1997; Liu et al. 2010) gave remarkable docking energies with both ERα and ERβ while the pregnane steroids 3,16-dihydroxypregn-5-en-20-one, 3,21-dihydroxypregna-5,16-dien-20-one, ergost-5-ene-3,26-diol, and pregnadienolone, showed selective docking with ERβ. Interestingly, neither furostane nor the spirostane steroids, common in Dioscorea spp. docked well with the estrogen receptors. It is worth noting that clinical studies have shown D. villosa to have little effect on menopausal symptoms (Komesaroff et al. 2001). Echinacea spp. Echinacea (E. angustifolia, E. pallida, and E. purpurea) is one of the most popular herbal supplements sold in the United States and has been used as a treatment for the common cold, coughs, bronchitis, upper respiratory infections, and inflammatory conditions (Percival 2000). Recent studies have demonstrated Echinacea to exhibit immune-system-stimulating activity (Block and Mead 2003). Phytochemicals that have been isolated from Echinacea spp. include chicoric acid and monomethyl and dimethyl ethers, trichocarpinine, cinnamoylechinadiol, cinnamoylechinaxanthol, cinnamoylepoxyechinadiol, cinnamoyldihydroxynardol, caftaric acid, caffeoyl-p-coumaroyltartaric acid, burkinabin A, burkinabin B, kaempferol, luteolin, and quercetin. Although Echinacea has not shown estrogenic activity (Zava et al. 1998), six phytochemicals were identified in this docking study that showed strong docking to the estrogen receptor: the flavonoid quercetin; the phenolic compounds caffeoyl-p-coumaroyltartaric acid, caftaric acid, and chicoric acid; and the sesquiterpenoids cinnamoylechinadiol and cinnamoylepoxyechinadiol. Gingko biloba G. biloba is commonly used as a supplement to improve cognitive abilities (Kennedy et al. 2000), and for women specifically, it has been used to treat some of the side effects accompanying menopause (Oh and Chung 2004). The extracts of G. biloba have previously been shown to exhibit feeble estrogenic effects, and act as selective estrogen receptor modulators (SERMs) with the α and β estrogen receptors (Oh and Chung 2006). Phytochemical analyses have revealed G. biloba extracts to contain the flavonoids (2R,3S,4S)-3,3′,4,4′,5,5′,7-heptahydroxyflavan, 8-(5-carboxy-2-methoxyphenyl)-5,7-dihydroxy-4′-methoxyflavone, acacetin, amentoflavone, apigenin, bilobetin, 5′-methoxybilobetin, epigallocatechin, ginkgetin, isoginkgetin, kaempferol, luteolin, quercetin, sciadopitysin, and tricetin, as well as the lignan sesamin, the sesquiterpenoids bilobanol, and the steroid globosterol. Of these, quercetin, tricetin, and sesamin gave large negative docking energies. Sesamin, in particular, docked strongly with ERβ. Glycyrrhiza glabra The phytochemistry of Glycyrrhiza glabra has been well studied and numerous compounds have been isolated and identified, including aurones (licoagroaurone and licoagrone), chalcones (1,2-dihydroparatocarpin A, 2,4,4′-trihydroxychalcone, 4-hydroxychalcone, cordifolin, isoliquiriteginin, licoagrochalcones A-D, licochalcones A and B, α,2′,4,4′-tetrahydroxydihydrochalcone, and kanzonol Y), coumarins (2′-O-methylglabridin, 3,4-didehydroglabridin, 3′-hydroxy-4′-methoxyglabridin, 4′-O-methylgrabridin, 4′-O-methylkanzonol W, bergapten, gancaonin F, glabrene, glabrocoumarin, glycycoumarin, glyinflanin H, hispaglabridin A, hispaglabridin B, isoglycycoumarin, isoglycyrol, kanzonol U, kanzonol V, kanzonol W, and licocoumarin A), flavonoids (3-hydroxyglabrol, 6-prenyleriodictyol, 6-prenylpinocembrin, folerogenin, glabranin, glabrol, isolicoflavonol, isoschaftoside, isoviolanthin, kaempferol, kumatakenin, licoagrodin, licoflavanone, naringenin, norwogonin, pinocembrin, quercetin, shinflavanone, and vitexin), isoflavonoids (1-methoxyphaseollin, 2′,4′,5,7-tetrahydroxy-3′,8-diprenylisoflavanone, 7-acetoxy-2-methylisoflavone, 7-hydroxy-2-methylisoflavone, 7-methoxy-2-methylisoflavone, 8-prenylphaseollinisoflavan, genistein, glabraisoflavanone A, glabraisoflavanone B, glabroisoflavanone A, glabridin, glabroisoflavanone B, glabrone, glyasperin B, glyasperin K, glyzaglabrin, glyzarin, isoderrone, isoglabrone, isomucronulatol, kanzonol R, kanzonol T, kanzonol X, licoagroside A, licoricidin, lupiwighteone, phaseollinisoflavan, prunetin, shinpterocarpin, tetrapterol G, and wighteone), the lignan licoagrocarpin, stilbenoids (3,3′,4,5′-tetrahydroxy-4′,5-diprenylbibenzyl, 3,3′,4,5′-tetrahydroxy-5-prenylbibenzyl, 3,3′,5′-trihydroxy-4-methoxy-5-prenylbibenzyl, 3,3′,5′-trihydroxy-4-methoxybibenzyl, 3,4′,5-trihydroxy-3′,4-diprenylbibenzyl, 3,4′,5-trihydroxy-3′-prenyldihydrostilbene, 3-acetoxy,4′,5-dihydroxy-3′-prenyldihydrostilbene, licoagrodione, and uralstilbene), and triterpenoids (11-deoxoglycyrrhetic acid, 18α-glycyrrhetic acid, 18α-hydroxyglycyrrhetic acid, 21-hydroxyisoglabrolide, 24-hydroxyglycyrrhetic acid, 24-hydroxyliquiritic acid, 28-hydroxyglycyrrhetic acid, 3,24-dihydroxy-11,13(18)-oleanadien-30-oic acid methyl ester, 3,24-dihydroxy-9(11),12-oleanadien-30-oic acid, 9(11)-dehydroglycyrrhetic acid, β-amyrin, betulinic acid, desoxoglabrolide, glabric acid, glabrolide, glycyrrhetic acid, glycyrrhetol, isoglabrolide, lanosta-5,24-dien-3-ol, liquiridiolic acid, liquiritic acid, and liquoric acid). Licorice (Glycyrrhiza glabra) root has been used for thousands of years by different cultures and for a variety of reasons (Fenwick et al. 1990). Although licorice root has been suggested as a treatment for symptoms of menopause (Ojeda 2003), G. glabra root extracts have been shown to be inactive in terms of ERα or ERβ binding (Liu et al. 2001). Nevertheless, however, fractionation of G. glabra extracts has revealed several ER-modulating components (Khalaf et al. 2010; Simons et al. 2011). Glabrene binds to human ER and shows estrogenic activity (Tamir et al. 2001; Simons et al. 2011). Our in-silico docking study shows glabrene to be a strongly docking ligand to both ERα and ERβ (−104.8 and −114.9 kJ/mol, respectively). Glabridin, on the other hand, displayed ERα-selective antagonism (Simons et al. 2011), in contrast to the docking results that showed glabridin to have ERβ docking selectivity (−15.8 and −92.9 kJ/mol, respectively). The chalcones isoliquiritigenin (Tamir et al. 2001; Maggiolini et al. 2002) and licochalcone A (Rafi et al. 2000), the flavonoid quercetin (Kuiper et al. 1998), and the isoflavonoid genistein (Ososki and Kennelly 2003) also bind to human ERα and show estrogenic activity. These ligands all show negative docking energies with ERα (range from −93.2 to −99.9 kJ/mol) and ERβ (−98.9 to −107.8 kJ/mol). Lepidium meyenii Maca (Lepidium meyenii) is native to the central Andes of Peru (3500–4500 m asl) (Wang et al. 2007). The root has been used by native Amerindians to improve fertility, as an aphrodisiac for both men and women. Maca was found to increase sperm counts and gonadal mass in a rat model (Chung et al. 2005), to improve copulatory performance of male mice and rats (Zheng et al. 2000; Cicero et al. 2001), and to increase litter size (Ruiz-Luna et al. 2005) and pregnancy rates in female mice (Kuo et al. 2003). In adult human males, maca treatment led to increased semen volume and sperm count (Gonzales et al. 2001) and increased sexual desire (Gonzales et al. 2002; Stone et al. 2009). In addition, maca reduced sexual dysfunction in postmenopausal women (Brooks et al. 2008) and inhibited estrogen-deficient osteoporosis in ovariectomized rats (Zhang et al. 2006), but maca extracts have not shown estrogenic activity (Brooks et al. 2008). None of the L. meyenii phytochemicals investigated in this in-silico study showed remarkable docking energies, consistent with the non-estrogenic activity previously reported. Ptychopetalum olacoides, P. uncinatum Muira puama (bark and root extracts of P. olacoides or P. uncinatum) has been used in Amazonian Brazil during highly stressful periods, to treat CNS-related ailments, neuromuscular problems, “nervous weakness”, sexual debility, frigidity, impotence, and rheumatism (Schultes and Raffauf 1990; Siqueira et al. 1998; Duke et al. 2009). Consistent with these traditional uses, P. olacoides ethanol root extract has shown memory retrieval improvement in young and aging mice (da Silva et al. 2004), in-vitro acetylcholine esterase inhibitory activity (Siqueira et al. 2003), and prevention of stress-induced hypothalamic-pituitary-adrenal hyperactivity (Piato et al. 2008). In addition, Muira puama formulations have demonstrated efficacy in treating male erectile dysfunction and low libido (Waynberg 1994) and low sex drive in women (Waynberg and Brewer 2000). A number of clerodane diterpenoids have been isolated from P. olacoides bark (Tang et al. 2008; Tang et al. 2009; Tang et al. 2011). Several of these have given excellent docking energies with the estrogen receptor, but two in particular, ptycho-6α,7α-diol and ptycholide IV had remarkable docking to ERβ (−122.9 and −114.7 kJ/mol, respectively). To our knowledge, the estrogenic effects of Muira puama have not been investigated. Rhodiola rosea R. rosea is reputed to strengthen the nervous system, fight depression, enhance memory, and improve energy levels (Brown et al. 2002), which has been attributed to adaptogenic properties of the herb (Spasov et al. 2000; Darbinyan et al. 2000). The flavonoids gossypetin, herbacetin, and rhodiolin, and the lignan (+)-lariciresinol, have been identified in R. rosea. Lariciresinol showed strong docking to both ERα and ERβ. There are conflicting reports on the potential estrogenic effects of R. rosea, however (Eagon et al. 2004; Kim et al. 2005). Sambucus nigra The bark, leaves, flowers, fruit, and root extracts of black elderberry (Sambucus nigra) have been used traditionally to treat respiratory ailments such as bronchitis, cough, influenza, and upper respiratory infections (Zakay-Rones et al. 1995; Krawitz et al. 2011). Compounds identified in S. nigra extracts include α-amyrin, β-amyrin, α-amyrone, betulin, campesterol, cycloartenol, lupeol, oleanolic acid, quercetin, ursolic acid, and dihydrodehydrodiconiferyl alcohol (9-acetate). Of these, only the flavonoid quercetin showed good docking to the estrogen receptor. To our knowledge, there are no reports on the estrogenic activity of S. nigra extracts. Silybum marianum Milk thistle, Silybum marianum, extracts (silymarin) have been used for centuries to treat liver diseases (Flora et al. 1998). Silymarin contains the flavonoids apigenin, chrysoeriol, cisilandrin, eriodictyol, isocisilandrin, isosilandrin A, isosilandrin B, isosilybin A, isosilybin B, isosilybin C, isosilybin D, isosilychristin, kaempferol, naringenin, neosilyhermin A, neosilyhermin B, quercetin, silandrin A, silandrin B, silyamandin, silybin A, silybin B 2,3-dehydrosilybin, silychristin, 2,3-dehydrosilychristin, silychristin B, silydianin, silyhermin, silymonin, taxifolin; the lignan (Z)-dehydrodiconiferyl alcohol; the steroids 24-methylenelanost-8-ene-3,25,28-triol and marianine; and the triterpenoids silymin A and silymin B. Silybin B (also called silibinin) has been shown to selectively bind to the ERβ receptor rather than ERα (Seidlová-Wuttke et al. 2003b) and previous docking studies have shown selective docking to ERβ over ERα (El-Shitany et al. 2010). In contrast, our current docking study revealed that neither silybin A nor silybin B gave negative docking energies with either ERα or ERα. Quercetin and taxifolin, on the other hand, gave docking energies comparable to zearalenone with ERβ (−106.0, −104.6, and −104.5 kJ/mol, respectively), and these flavonoids have also shown estrogenic activity (Plíšková et al. 2005). Silymarin modulation of ERβ may be responsible for the estrogenic effects of the extract (Seidlová-Wuttke et al. 2003b; Plíšková et al. 2005; El-Shitany et al. 2010). Tribulus terrestris Tribulus terrestris has been used to contribute to physical and sexual strength (De Combarieu et al. 2003; Neychev and Mitev 2005). The plant is rich in steroidal glycosides (Wu et al. 1996; Yan et al. 1996; De Combarieu et al. 2003; Dinchev et al. 2008) and T. terrestris extracts have shown androgenic effects in animal models (Gauthaman et al. 2002; Gauthaman and Ganesan 2008), but had no influence on androgen production (Neychev and Mitev 2005) or gains in strength or muscle mass (Rogerson et al. 2007) in young men. To our knowledge, there have been no reports on the estrogenic effects of T. terrestris. In addition to steroids, several alkaloids (N-trans-feruloyltyramine, perlolyrine, terresoxazine, terrestriamide, tribulusamide A and tribulusamide B) and flavonoids (isorhamnetin, kaempferol, and quercetin) have been found in T. terrestris. In-silico molecular docking has shown that N-trans-feruloyltyramine docks strongly to both ERα and ERβ, perlolyrine docks strongly to ERβ, terrestribisamide docks strongly to ERα, quercetin docks strongly to ERβ, and the steroid 2,3-dihydroxypregn-16-en-20-one docks strongly to both ERα and ERβ. Trifolium pratense The alkaloids cis- and trans-clovamide, several coumarins, flavonoids, and isoflavonoids have been identified in red clover (Trifolium pratense). Although there have been no ethnobotanical reports to support it, the presence of isoflavones has led to suggest that red clover may serve as a phytochemical alternative to post-menopausal hormone replacement therapy (Coon et al. 2007). Red clover extract has been shown to exhibit weakly estrogenic effects in a rat model (Burdette et al. 2002) and does show in-vitro ERα and ERβ binding ability (Dornstauder et al. 2001; Beck et al. 2003; Overk et al. 2005). The clinical effectiveness of red clover has not, however, been demonstrated (Fugh-Berman and Kronenberg 2001; Booth et al. 2006). Biochanin A, daidzein, formononetin, and genistein have been identified as T. pratense isoflavonoids with ERα and ERβ binding activity (Beck et al. 2003; Pfitscher et al. 2008) and these compounds did show a binding preference for ERβ. Molecular docking of these ligands also showed preference for ERβ. They were not, however, the best docking ligands of T. pratense phytochemicals. The alkaloids cis- and trans-clovamide, and the coumarin pratenol B showed good docking to both ERα and ERβ. Of the T. pratense isoflavonoids, calycosin and pseudobaptigenin had more exothermic docking energies to ERβ than did biochanin A, daidzein, formononetin, or genistein. Trigonella foenum-graecum Fenugreek (Trigonella foenum-graecum) is used as an antidiabetic (Abdel-Barry et al. 1997) and for lowering blood lipid and cholesterol levels (Prasanna 2000). Phytochemicals identified in T. foenum-graecum include the alkaloid gentianine; the coumarins 7-acetoxy-4-methylcoumarin, trigocoumarin, and trigoforin; the flavonoids 2″-O-p-coumaroylvitexin, 2″-O-p-coumaroylorientin, 6,8-digalactosylapigenin, 8-galactopyranosyl-6-quinovopyranosylapigenin, 8-β-D-galactopyranosyl-6-β-D-xylopyranosylapigenin, isoorientin, isovitexin, kaempferol, luteolin, neocorymboside, orientin, quercetin, trigraecum, vicenin 1, vicenin 2, and vicenin 3; the isoflavonoid 3′,5,7-trihydroxy-5′-methoxyisoflavone; p-coumaric acid; and the steroids 25R-spirosta-3,5-diene (= Δ3,5-deoxyneotigogenin), diosgenin, furostane-2,3,22,26-tetrol (= trigoneoside aglycone), gitogenin, neotigogenin, tigogenin, and yamogenin. Quercetin and 3′,5,7-trihydroxy-5′-methoxyisoflavone were the strongest docking ligands for ERβ. Fenugreek seed extract has shown estrogenic activity (Sreeja et al. 2010). Turnera aphrodisiaca (syn. T. diffusa) Damiana leaf has been used in traditional medicine in Neotropical cultures as a stimulant and aphrodisiac (Morton 1981), and the herb is marketed as a sexual enhancer for men and women. Extracts of T. aphrodisiaca have shown aphrodisiac activity in mouse (Helmrick and Reiser 2000; Kumar et al. 2009) and rat (Estrada-Reyes et al. 2009; Estrada-Reyes et al. 2013) models, as well as anxiolytic activity in mice (Kumar and Sharma 2005). Phytochemical investigations have revealed damiana to be rich in flavonoids, including acacetin, apigenin, gonzalitosin, laricitrin, luteolin-8-propenoic acid, orientin, orientin-3″-ketone, pinocembrin, and syringetin (Piacente et al. 2002; Zhao et al. 2007). The extract has shown anti-aromatase (due primarily to acacetin and pinocembrin) and estrogenic activity (due primarily to apigenin, Z-echinacin, and pinocembrin) (Zhao et al. 2008). In contrast to these experimental results, molecular docking revealed the strongest docking Turnera compound to be luteolin-8-propenoic acid (−113.1, −123.1 kJ/mol for ERα and ERβ, respectively) but this compound was inactive in the aromatase and estrogen assays. In contrast, pinocembrin, which was active in both experimental assays, had ERα and ERβ docking energies of −81.4 and −87.9 kJ/mol, respectively. Vitex agnus-castus Vitex agnus-castus, “chaste tree”, has been used as a tonic for female reproductive disorders, including menstrual disorders (amenorrhea, dysmenorrhea), premenstrual syndrome (PMS, corpus luteum insufficiency, hyperprolactinemia, infertility, menopause, and disrupted lactation (Daniele et al. 2005; van Die et al. 2013). V. agnus-castus extracts contain several flavonoids (casticin, isoorientin, 6″-caffeoylisoorientin, 6″-caffeoylisoorientin-4″-methyl ether, isovitexin, luteolin, orientin, santin, 5-O-demethyltangeretin, and vitexin), diterpenoids (8,14-labdadiene-6,7,13-triol-6,7-diacetate, viteagnuside A, viteagnusin A, viteagnusin B, viteagnusin D, viteagnusin E, viteagnusin F, viteagnusin G, viteagnusin H, viteagnusin I, viteagnusin J, and vitexlactam A), the isoflavonoid vitexcarpan, and the phenolic compounds agnucastoside C and agnuside. Based upon docking energies with ERβ, casticin, santin, vitexcarpan, and the aglycones of agnucastoside C and agnuside, could be expected to exhibit ER modulation. Indeed, V. agnus-castus extracts have been shown to bind to ERα (Liu et al. 2001) and ERβ (Jarry et al. 2003). It has been suggested that V. agnus-castus is a source of 3-ketosteroids with progesterone-like activity (Bruneton 1999; Dweck 2006). Go to: Conclusions This molecular docking study has revealed that almost all popular herbal supplements contain phytochemical components that may bind to the human estrogen receptor and exhibit selective estrogen receptor modulation. As such, these herbal supplements may cause unwanted side effects related to estrogenic activity. For example, estrogenic agents may be effective and potent growth stimulators of estrogen-receptor positive tumors and pose a hazard to patients with breast cancer who have ER-positive tumors and who are being treated with antiestrogens. The strongest docking (most exothermic docking energies) phytochemical ligands were phenolic compounds and the weakest docking ligands were triterpenoids. A common binding motif for phenolic ligands in ERα is the hydrophobic pocket of Leu 387, Phe 404, Met 388, and Leu 391, along with edge-to-face π–π interactions with Phe 404, and hydrogen bonds between the phenolic –OH group and the guanidine group of Arg 394 and the carboxylate of Glu 356. Similarly, interactions of phenolic ligands with ERβ include binding in a hydrophobic pocket formed by Leu 298, Leu 339, and Phe 356; edge-to-face π–π interactions of the phenolic ligand with Phe 356 and hydrogen boding of the phenolic –OH group with the carbonyl group of Leu 339, the guanidine group of Arg 346, and the carboxylate of Glu 305. Common hydrogen-bonding residues in the binding sites of ERα are the guanidine group of Arg 394, the imidazole group of His 524. Hydrogen-bonding residues of ERβ are Arg 346, His 475. There are several limitations to these docking results: Some of the herbal phytochemicals examined may not be bioavailable due to limited solubility, membrane permeability; This docking study has only examined docking of the natural ligands (or their aglycones) and does not take into account in-vivo hydrolysis or other metabolic derivatization; The docking studies do not account for synergism in the estrogen receptor binding; The molecular docking method itself suffers from inherent limitations (e.g., the protein is modeled as a rigid structure without flexibility, solvation of the binding site and the ligand is excluded, and free-energy estimation of protein-ligand complexes is largely ignored) (Yuriev et al. 2011; Yuriev and Ramsland 2013). Docking energies do not provide information about whether strongly binding ligands may function as agonists or antagonists of the estrogen receptor. Go to: Acknowledgments CNP is grateful to the University of Alabama in Huntsville Department of Chemistry for an undergraduate summer research fellowship. Go to: Abbreviations ERα Estrogen receptor α ERβ Estrogen receptor β ER+ Estrogen receptor positive MMFF Merck molecular force field PDB Protein data bank PMS Pre-menstrual syndrom RMSD Root mean square deviation SERM Selective estrogen receptor modulator Go to: Footnotes Competing interests The authors declare that they have no competing interests. Authors’ contributions CNP and WNS conceived the study, carried out molecular mechanics and molecular docking, and drafted and revised the manuscript. Both authors read and approved the final manuscript. Go to: Contributor Information Chelsea N Powers, Email: ude.hau@7000pnc. William N Setzer, Email: ude.hau.yrtsimehc@reztesw. Go to: References Abdel-Barry JA, Abdel-Hassan IA, Al-Hakiem MHH. Hypoglycaemic and antihyperglycaemic effects of Trigonella foenum-graecum leaf in normal and alloxan induced diabetic rats. J Ethnopharmacol. 1997;58:149–155. doi: 10.1016/S0378-8741(97)00101-3. [PubMed] [Cross Ref] Al-Bareeq RJ, Ray AA, Nott L, Pautler SE, Razvi H. Dong quai (Angelica sinensis) in the treatment of hot flashes for men on androgen deprivation therapy: results of a randomized double-blind placebo controlled trial. Can Urol Assoc J. 2010;4:49–53. [PMC free article] [PubMed] Ali Z, Smillie TJ, Khan IA. Cholestane steroid glycosides from the rhizomes of Dioscorea villosa. Carbohydr Res. 2013;370:86–91. doi: 10.1016/j.carres.2012.12.022. [PMC free article] [PubMed] [Cross Ref] Au AM, Ko R, Boo FO, Hsu R, Perez G, Yang Z. Screening methods for drugs and heavy metals in Chinese patent medicines. Bull Environ Contam Toxicol. 2000;65:112–119. doi: 10.1007/s001280000102. [PubMed] [Cross Ref] Beck V, Unterrieder E, Krenn L, Kubelka W, Jungbauer A. Comparison of hormonal activity (estrogen, androgen, and progestin) of standardized plant extracts for large scale use in hormone replacement therapy. Mol Biol. 2003;84:259–268. [PubMed] Block KI, Mead MN. Immune system effects of echinacea, ginseng, and astragalus: a review. Integr Cancer Ther. 2003;2:247–267. doi: 10.1177/1534735403256419. [PubMed] [Cross Ref] Booth NL, Piersen CE, Banuvar S, Geller SE, Shulman LP, Farnsworth NR. Clinical studies of red clover (Trifolium pratense) dietary supplements in menopause: a literature review. Menopause. 2006;13:251–264. doi: 10.1097/01.gme.0000198297.40269.f7. [PubMed] [Cross Ref] Brooks NA, Wilcox G, Walker KZ, Ashton JF, Cox MB, Stojanovska L. Beneficial effects of Lepidium meyenii (Maca) on psychological symptoms and measures of sexual dysfunction in postmenopausal women are not related to estrogen or androgen content. Menopause. 2008;15:1157–1162. doi: 10.1097/gme.0b013e3181732953. [PubMed] [Cross Ref] Brown RP, Gerbarg PL, Ramazanov Z. Rhodiola rosea A phytomedicinal overview. Herbal Gram. 2002;56:40–52. Bruneton J. Pharmacognosy Phytochemistry Medicinal Plants. 2. Paris: Lavoisier; 1999. Burdette JE, Liu J, Lantvit D, Lim E, Booth N, Bhat KPL, Hedayat S, Van Breemen RB, Constantinou AI, Pezzuto JM, Farnsworth NR, Bolton JL. Trifolium pratense (red clover) exhibits estrogenic effects in vivo in ovariectomized Sprague–Dawley rats. J Nutr. 2002;132:27–30. [PubMed] Burdette JE, Liu J, Chen SN, Fabricant DS, Piersen CE, Barker EL, Pezzuto JM, Mesecar A, van Breemen RB, Farnsworth NR, Bolton JL. Black cohosh acts as a mixed competitive ligand and partial agonist of the serotonin receptor. J Agric Food Chem. 2003;51:5661–5670. doi: 10.1021/jf034264r. [PubMed] [Cross Ref] Calixto JB. Efficacy, safety, quality control, marketing and regulatory guidelines for herbal medicines (phytotherapeutic agents) Braz J Med Biol Res. 2000;33:179–189. doi: 10.1590/S0100-879X2000000200004. [PubMed] [Cross Ref] Chen SN, Fabricant DS, Lu ZZ, Zhang H, Fong HHS, Farnsworth NR. Cimiracemates A-D, phenylpropanoids esters from the rhizomes of Cimicifuga racemosa. Phytochemistry. 2002;61:409–413. doi: 10.1016/S0031-9422(02)00209-1. [PubMed] [Cross Ref] Chung F, Rubio J, Gonzales C, Gasco M, Gonzales GF. Dose–response effects of Lepidium meyenii (Maca) aqueous extract on testicular function and weight of different organs in adult rats. J Ethnopharmacol. 2005;98:143–147. doi: 10.1016/j.jep.2005.01.028. [PubMed] [Cross Ref] Chye PLH. Traditional Asian folklore medicines in sexual health. Indian J Urol. 2006;22:241–245. doi: 10.4103/0970-1591.27632. [Cross Ref] Cicero AFG, Bandieri E, Arletti R. Lepidium meyenii Walp. improves sexual behavior in male rats independently from its action on spontaneous locomotor activity. J Ethnopharmacol. 2001;75:225–229. doi: 10.1016/S0378-8741(01)00195-7. [PubMed] [Cross Ref] Coon JT, Pittler MH, Ernst E. Trifolium pratense isoflavones in the treatment of menopausal hot flushes: a systematic review and meta-analysis. Phytomedicine. 2007;14:153–159. doi: 10.1016/j.phymed.2006.12.009. [PubMed] [Cross Ref] Cupp MJ. Herbal remedies: adverse effects and drug interactions. Am Fam Physician. 1999;59:1239–1244. [PubMed] da Silva AL, Piato ALS, Bardini S, Netto CA, Nunes DS, Elisabetsky E. Memory retrieval improvement by Ptychopetalum olacoides in young and aging mice. J Ethnopharmacol. 2004;95:199–203. doi: 10.1016/j.jep.2004.07.019. [PubMed] [Cross Ref] Daniele C, Thompson Coon J, Pittler MH, Ernst E. Vitex agnus castus. Drug Saf. 2005;28:319–322. doi: 10.2165/00002018-200528040-00004. [PubMed] [Cross Ref] Darbinyan V, Kteyan A, Panossian A, Gabrielian E, Wikman G, Wagner H. Rhodiola rosea in stress induced fatigue – A double blind cross-over study of a standardized extract SHR-5 with a repeated low-dose regimen on the mental performance of healthy physicians during night duty. Phytomedicine. 2000;7:365–371. doi: 10.1016/S0944-7113(00)80055-0. [PubMed] [Cross Ref] De Combarieu E, Fuzzati N, Lovati M, Mercalli E. Furostanol saponins from Tribulus terrestris. Fitoterapia. 2003;74:583–591. doi: 10.1016/S0367-326X(03)00152-7. [PubMed] [Cross Ref] Dictionary of Natural Products on DVD v.23:1 (2014) CRC Press, Boca Raton, Florida, USA Dinchev D, Janda B, Evstatieva L, Oleszek W, Aslani MR, Kostova I. Distribution of steroidal saponins in Tribulus terrestris from different geographical regions. Phytochemistry. 2008;69:176–186. doi: 10.1016/j.phytochem.2007.07.003. [PubMed] [Cross Ref] Dornstauder E, Jisa E, Unterrieder I, Krenn L, Kubelka W, Jungbauer A. Estrogenic activity of two standardized red clover extracts (Menoflavon®) intended for large scale use in hormone replacement therapy. J Steroid Biochem Mol Biol. 2001;78:67–75. doi: 10.1016/S0960-0760(01)00075-9. [PubMed] [Cross Ref] Duke J, Bogenschutz MJ (1998) Dr. Duke’s phytochemical and ethnobotanical databases. http://www.ars-grin.gov/duke/ Duke JA, Bogenschutz-Godwin MJ, Ottesen AR. Duke’s Handbook of Medicinal Plants of Latin America. Boca Raton, Florida: CRC Press; 2009. pp. 587–588. Dutta B. Food and medicinal values of certain species of Dioscorea with special reference to Assam. J Pharmacog Phytochem. 2015;3:15–18. Dweck AC. Isoflavones, phytohormones and phytosterols. J Appl Costetol. 2006;24:17–33. Eagon PK, Elm MS, Gerbarg PL, Brown RP, Check JJ, Diorio GJ, Houghton F. Evaluation of the medicinal botanical Rhodiola rosea for estrogenicity. Proc Am Assoc Cancer Res. 2004;45:663. El-Shitany NA, Hegazy S, El-desoky K. Evidences for antiosteoporotic and selective estrogen receptor modulator activity of silymarin compared with ethinyl estradiol in ovariectomized rats. Phytomedicine. 2010;17:116–125. doi: 10.1016/j.phymed.2009.05.012. [PubMed] [Cross Ref] Estrada-Reyes R, Ortiz-López P, Gutiérrez-Ortíz J, Martínez-Mota L. Turnera diffusa Wild (Turneraceae) recovers sexual behavior in sexually exhausted males. J Ethnopharmacol. 2009;123:423–429. doi: 10.1016/j.jep.2009.03.032. [PubMed] [Cross Ref] Estrada-Reyes R, Carro-Juárez M, Martínez-Mota L. Pro-sexual effects of Turnera diffusa Wild (Turneraceae) in male rats involves the nitric oxide pathway. J Ethnopharmacol. 2013;146:164–172. doi: 10.1016/j.jep.2012.12.025. [PubMed] [Cross Ref] Fang L, Xiao XF, Liu CX, He X. Recent advance in studies on Angelica sinensis. Chinese Herbal Med. 2012;4:12–25. Fenwick GR, Lutomski J, Nieman C. Liquorice, Glycyrrhiza glabra L. – Composition, uses and analysis. Food Chem. 1990;38:119–143. doi: 10.1016/0308-8146(90)90159-2. [Cross Ref] Flora K, Hahn M, Rosen H, Benner K. Milk thistle (Silybum marianum) for the therapy of liver disease. Am J Gastroenterol. 1998;93:139–143. doi: 10.1111/j.1572-0241.1998.00139.x. [PubMed] [Cross Ref] Fugh-Berman A, Kronenberg F. Red clover (Trifolium pratense) for menopausal women: current state of knowledge. Menopause. 2001;8:333–337. doi: 10.1097/00042192-200109000-00007. [PubMed] [Cross Ref] Gauthaman K, Ganesan AP. The hormonal effects of Tribulus terrestris and its role in the management of male erectile dysfunction – an evaluation using primates, rabbit and rat. Phytomedicine. 2008;15:44–54. doi: 10.1016/j.phymed.2007.11.011. [PubMed] [Cross Ref] Gauthaman K, Adaikan PG, Prasad RNV. Aphrodisiac properties of Tribulus terrestris extract (Protodioscin) in normal and castrated rats. Life Sci. 2002;71:1385–1396. doi: 10.1016/S0024-3205(02)01858-1. [PubMed] [Cross Ref] Gonzales GF, Córdova A, Gonzales C, Chung A, Vega K, Villena A. Lepidium meyenii (maca) improved semen parameters in adult men. Asian J Androl. 2001;3:301–304. [PubMed] Gonzales GF, Córdova A, Vega K, Chung A, Villena A, Góñez C, Castillo S. Effect of Lepidium meyenii (MACA) on sexual desire and its absent relationship with serum testosterone levels in adult healthy men. Andrologia. 2002;34:367–372. doi: 10.1046/j.1439-0272.2002.00519.x. [PubMed] [Cross Ref] Halgren TA. Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF 94. J Comput Chem. 1996;17:490–519. doi: 10.1002/(SICI)1096-987X(199604)17:5/6<490::AID-JCC1>3.0.CO;2-P. [Cross Ref] Helmrick L, Reiser C. Aphrodisiac properties of Turnera diffusa. J Undergrad Res. 2000;3:109–113. Hirata JD, Swiersz LM, Zell B, Small R, Ettinger B. Does dong quai have estrogenic effects in postmenopausal women? A double-blind, placebo-controlled trial. Fertil Steril. 1997;68:981–986. doi: 10.1016/S0015-0282(97)00397-X. [PubMed] [Cross Ref] Jarry H, Spengler B, Porzel A, Schmidt J, Wuttke W, Christoffel V. Evidence for estrogen receptor β-selective activity of Vitex agnus-castus and isolated flavones. Planta Med. 2003;69:945–947. doi: 10.1055/s-2003-45105. [PubMed] [Cross Ref] Kennedy DO, Scholey AB, Wesnes KA. The dose-dependent cognitive effects of acute administration of Ginkgo biloba to healthy young volunteers. Psychopharmacology (Berl) 2000;151:416–423. doi: 10.1007/s002130000501. [PubMed] [Cross Ref] Khalaf I, Vlase L, Lazăr D, Corciovă A, Ivãnescu B, Lazăr MI. HPLC-MS study of phytoestrogens from Glycyrrhiza glabra. Farmacia. 2010;58:89–94. Kim JB, Price RK, Stein RC, O’Hara MJ. Risks of the complementary medicine Rhodiola rosea in breast cancer. J Nutr (Suppl) 2005;135:3048–3049. Komesaroff PA, Black CVS, Cable V, Sudhir K. Effects of wild yam extract on menopausal symptoms, lipids and sex hormones in healthy menopausal women. Climacteric. 2001;4:144–150. doi: 10.1080/cmt.4.2.144.150. [PubMed] [Cross Ref] Komori T. Glycosides from Dioscorea bulbifera. Toxicon. 1997;35:1531–1536. doi: 10.1016/S0041-0101(97)00032-9. [PubMed] [Cross Ref] Krawitz C, Abu Mraheil M, Stein M, Imirzalioglu C, Domann E, Pleschka S, Hain T. Inhibitory activity of a standardized elderberry liquid extract against clinically-relevant human respiratory bacterial pathogens and influenza A and B viruses. BMC Complement Altern Med. 2011;11:16. doi: 10.1186/1472-6882-11-16. [PMC free article] [PubMed] [Cross Ref] Kretzschmar G, Zierau O, Wober J, Tischer S, Metz P, Vollmer G. Prenylation has a compound specific effect on the estrogenicity of naringenin and genistein. J Steroid Biochem Mol Biol. 2010;118:1–6. doi: 10.1016/j.jsbmb.2009.08.005. [PubMed] [Cross Ref] Kronenberg F. Black cohosh, a menopausal remedy, does not have estrogenic activity and does not promote breast cancer cell growth. Int J Oncol. 2003;23:1407–1412. [PubMed] Kuiper GGJM, Lemmen JG, Carlsson B, Corton JC, Safe SH, van der Saag PT, van der Burg B, Gustafsson JÅ. Interaction of estrogenic chemicals and phytoestrogens with estrogen receptor β Endocrinology. 1998;139:4252–4263. [PubMed] Kumar MHV, Gupta YK. Effect of different extracts of Centella asiatica on cognition and markers of oxidative stress in rats. J Ethnopharmacol. 2002;79:253–260. doi: 10.1016/S0378-8741(01)00394-4. [PubMed] [Cross Ref] Kumar S, Sharma A. Anti-anxiety activity studies on homeopathic formulations of Turnera aphrodisiaca Ward. Evid Based Complement Alternat Med. 2005;2:117–119. doi: 10.1093/ecam/neh069. [PMC free article] [PubMed] [Cross Ref] Kumar S, Madaan R, Sharma A. Evaluation of aphrodisiac activity of Turnera aphrodisiaca. Int J Pharmacog Phytochem Res. 2009;1:1–4. Kuo TF, Chang MH, Liau MY. Effects of Lepidium meyenii Walp. (maca) on fecundity and puppy growth in mice. Taiwan Vet J. 2003;29:1–8. Lau CB, Ho TC, Chan TW, Kim SC. Use of dong quai (Angelica sinensis) to treat peri- or postmenopausal symptoms in women with breast cancer: is it appropriate? Menopause. 2005;12:734–740. doi: 10.1097/01.gme.0000184419.65943.01. [PubMed] [Cross Ref] Li W, Sun Y, Liang W, Fitzloff JF, van Breeman RB. Identification of caffeic acid derivatives in Actaea racemosa (Cimicifuga racemosa, black cohosh) by liquid chromatography / tandem mass spectrometry. Rapid Commun Mass Spectrom. 2003;17:978–982. doi: 10.1002/rcm.1008. [PubMed] [Cross Ref] Lieberman S. A review of the effectiveness of Cimicifuga racemosa (black cohosh) for the symptoms of menopause. J Womens Health. 1998;7:525–529. doi: 10.1089/jwh.1998.7.525. [PubMed] [Cross Ref] Liske E, Hänggi W, Henneicke-von Zepelin HH, Boblitz N, Wüstenberg P, Rahlfs VW. Physiological investigation of a unique extract of black cohosh (Cimicifugae racemosae rhozoma): A 6-month clinical study demonstrates no systemic estrogeni effect. J Womens Health Gend Based Med. 2002;11:163–174. doi: 10.1089/152460902753645308. [PubMed] [Cross Ref] Liu J, Burdette JE, Xu H, Gu C, van Breemen RB, Bhat KPL, Booth N, Constantinou AI, Pezzuto JM, Fong HHS, Farnsworth NR, Bolton JL. Evaluation of estrogenic activity of plant extracts for the potential treatment of menopausal symptoms. J Agric Food Chem. 2001;49:2472–2479. doi: 10.1021/jf0014157. [PubMed] [Cross Ref] Liu H, Chou GX, Guo YL, Ji LL, Wang JM, Wang ZT. Norclerodane diterpenoids from rhizomes of Dioscorea bulbifera. Phytochemistry. 2010;71:1174–1180. doi: 10.1016/j.phytochem.2010.04.002. [PubMed] [Cross Ref] Lupu R, Mehmi I, Atlas E, Tsai MS, Pisha E, Oketch-Rabah HA, Nuntanakorn P, Kennelly EJ, Kronenberg F. Black cohosh, a menopausal remedy, does not have estrogenic activity and does not promote breast cancer cell growth. Int J Oncol. 2003;23:1407–1412. [PubMed] Maggiolini M, Statti G, Vivacqua A, Gabriele S, Rago V, Loizzo M, Menichini F, Amdò S. Estrogenic and antiproliferative activities of isoliquiritigenin in MCF7 breast cancer cells. J Steroid Biochem Mol Biol. 2002;82:315–322. doi: 10.1016/S0960-0760(02)00230-3. [PubMed] [Cross Ref] Mahady GB. Is black cohosh estrogenic? Nutr Rev. 2003;61:183–186. doi: 10.1301/nr.2003.may.183-186. [PubMed] [Cross Ref] Malamas MS, Manas ES, McDevitt RE, Gunawan I, Xu ZB, Collini MD, Miller CP, Dinh T, Henderson RA, Keith JC, Jr, Harris HA. Design and synthesis of aryl diphenolic azoles as potent and selective estrogen receptor-beta ligands. J Med Chem. 2004;47:5021–5040. doi: 10.1021/jm049719y. [PubMed] [Cross Ref] Manas ES, Unwalla RJ, Xu ZB, Malamas MS, Miller CP, Harris HA, Hsiao C, Akopian T, Hum WT, Malakian K, Wolfrom S, Bapat A, Bhat RA, Stahl ML, Somers WS, Alvarez JC. Structure-based design of estrogen receptor-beta selective ligands. J Am Chem Soc. 2004;126:15106–15119. doi: 10.1021/ja047633o. [PubMed] [Cross Ref] Manas ES, Xu ZB, Unwalla RJ, Somers WS. Understanding the selectivity of genistein for human estrogen receptor-beta using X-ray crystallography and computational methods. Structure. 2004;12:2197–2207. doi: 10.1016/j.str.2004.09.015. [PubMed] [Cross Ref] McKenna DJ, Jones K, Humphrey S, Hughes K. Black cohosh: efficacy, safety, and use in clinical and preclinical applications. Altern Ther Health Med. 2001;7:93–100. [PubMed] Molegro Virtual Docker version 6.0.0 (2013) Molegro ApS, Aarhus, Denmark Morton JF. Atlas of Medicinal Plants of Middle America. Vol I. Springfield, Illinois: Charles C. Thomas Publisher; 1981. p. 584. Neychev VK, Mitev VI. The aphrodisiac herb Tribulus terrestris does not influence the androgen production in young men. J Ethnopharmacol. 2005;101:319–323. doi: 10.1016/j.jep.2005.05.017. [PubMed] [Cross Ref] Oh SM, Chung KH. Estrogenic effects of Ginkgo biloba extracts. Life Sci. 2004;74:1325–1335. doi: 10.1016/j.lfs.2003.06.045. [PubMed] [Cross Ref] Oh SM, Chung KH. Antiestrogenic activities of Ginkgo biloba extracts. J Steroid Biochem Mol Biol. 2006;100:167–176. doi: 10.1016/j.jsbmb.2006.04.007. [PubMed] [Cross Ref] Ojeda L. Menopause without Medicine. 5. Alameda, California: Hunter House; 2003. Ososki AL, Kennelly EJ. Phytoestrogens: a review of the present state of research. Phytother Res. 2003;17:845–869. doi: 10.1002/ptr.1364. [PubMed] [Cross Ref] Overk CR, Yao P, Chadwick LR, Nikolic D, Sun Y, Cuendet MA, Deng Y, Hedayat AS, Pauli GF, Farnsworth NR, van Breemen RB, Bolton JL. Comparison of the in vitro estrogenic activities of compounds from hops (Humulus lupulus) and red clover (Trifolium pratense) J Agric Food Chem. 2005;53:6246–6253. doi: 10.1021/jf050448p. [PMC free article] [PubMed] [Cross Ref] Percival SS. Use of Echinacea in medicine. Biochem Pharmacol. 2000;60:155–158. doi: 10.1016/S0006-2952(99)00413-X. [PubMed] [Cross Ref] Pfitscher A, Reiter E, Jungbauer A. Receptor binding and transactivation activities of red clover isoflavones and their metabolites. J Steroid Biochem Mol Biol. 2008;112:87–94. doi: 10.1016/j.jsbmb.2008.08.007. [PubMed] [Cross Ref] Piacente S, Camargo EES, Zampelli A, Gracioso JS, Souza Brito AR, Pizza C, Vilegas W. Flavonoids and arbutin from Turnera diffusa. Z Naturforsch. 2002;57c:983–985. [PubMed] Piato AL, Detanico BC, Jesus JF, Rodrigues Lhullier FL, Nunes DS, Elisabetsky E. Effects of Marapuama in the chronic mild stress model: further indication of antidepressant properties. J Ethnopharmacol. 2008;118:300–304. doi: 10.1016/j.jep.2008.04.018. [PubMed] [Cross Ref] Plíšková M, Vondráček J, Křen V, Gažák R, Sedmera P, Walterová D, Psotová J, Šimánek V, Machala M. Effects of silymarin flavonolignans and synthetic silybin derivatives on estrogen and aryl hydrocarbon receptor activation. Toxicology. 2005;215:80–89. doi: 10.1016/j.tox.2005.06.020. [PubMed] [Cross Ref] Pockaj BA, Loprinzi CL, Sloan JA, Novotny PJ, Barton DL, Hagenmaier A, Zhang H, Lambert GH, Reeser KA, Wisbey JA. Pilot evaluation of black cohosh for the treatment of hot flashes in women. Cancer Invest. 2004;22:515–521. doi: 10.1081/CNV-200026394. [PubMed] [Cross Ref] Prasanna M. Hypolipidemic effect of fenugreek: a clinical study. Indian J Pharmacol. 2000;32:34–36. Rafi MM, Rosen RT, Vassil A, Ho CT, Zhang H, Ghai G, Lambert G, DiPaola RS. Modulation of bcl-2 and cytotoxicity of licochalcone-A, a novel estrogenic flavonoid. Anticancer Res. 2000;20:2653–2658. [PubMed] Rao SB, Chetana M, Devi PU. Centella asiatica treatment during postnatal period enhances learning and memory in mice. Physiol Behav. 2005;86:449–457. doi: 10.1016/j.physbeh.2005.07.019. [PubMed] [Cross Ref] Rhyu MR, Lu J, Webster DE, Fabricant DS, Farnsworth NR, Wang ZJ. Black cohosh (Actaea racemosa, Cimicifuga racemosa) behaves as a mixed competitive ligand and partial agonist at the human μ opiate receptor. J Agric Food Chem. 2006;54:9852–9857. doi: 10.1021/jf062808u. [PMC free article] [PubMed] [Cross Ref] Rogerson S, Riches CJ, Jennings C, Weatherby RP, Meir RA, Marshall-Gradisnik SM. The effect of five weeks of Tribulus terrestris supplementation on muscle strength and body composition during preseason training in elite rugby league players. J Strength Cond Res. 2007;21:348–353. [PubMed] 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:16. doi: 10.1186/1477-7827-3-16. [PMC free article] [PubMed] [Cross Ref] Sautour M, Miyamoto T, Lacaille-Dibois MA. Steroidal saponins and flavan-3-ol glycosides from Dioscorea villosa. Biochem System Ecol. 2006;34:60–63. doi: 10.1016/j.bse.2005.07.007. [Cross Ref] Sautour M, Mitaine-Offer AC, Lacaille-Dubois MA. The Dioscorea genus: a review of bioactive steroid saponins. J Nat Med. 2007;61:91–101. doi: 10.1007/s11418-006-0126-3. [Cross Ref] Schultes RE, Raffauf RF. The Healing Forest. Portland, Oregon: Dioscorides Press; 1990. p. 343. Seidlová-Wuttke D, Hesse O, Jarry H, Christoffel V, Spengler B, Becker T, Wuttke W. Evidence for selective estrogen receptor modulator activity in a black cohosh (Cimicifuga racemosa) extract: comparison with estradiol-17β Eur J Endocrinol. 2003;149:351–362. doi: 10.1530/eje.0.1490351. [PubMed] [Cross Ref] Seidlová-Wuttke D, Becker T, Christoffel V, Jarry H, Wuttke W. Silymarin is a selective estrogen receptor β (ERβ) agonist and has estrogenic effects in the metaphysis of the femur but no or antiestrogenic effects in the uterus of ovariectomized (ovx) rats. J Steroid Biochem Mol Biol. 2003;86:179–188. doi: 10.1016/S0960-0760(03)00270-X. [PubMed] [Cross Ref] Shao Y, Harris A, Wang M, Zhang H, Cordell GA, Bowman M, Lemmo E. Triterpene glycosides from Cimicifuga racemosa. J Nat Prod. 2000;63:905–910. doi: 10.1021/np000047y. [PubMed] [Cross Ref] Shiau AK, Barstad D, Loria PM, Cheng L, Kushner PJ, Agard DA, Greene GL. The structural basis of estrogen receptor/coactivator recognition and the antagonism of this interaction by tamoxifen. Cell. 1998;95:927–937. doi: 10.1016/S0092-8674(00)81717-1. [PubMed] [Cross Ref] Simons R, Vincken JP, Mol LAM, The SAM, Bovee TFH, Luijendijk TJC, Verbruggen MA, Gruppen H. Agonistic and antagonistic estrogens in licorice root (Glycyrrhiza glabra) Anal Bioanal Chem. 2011;401:305–313. doi: 10.1007/s00216-011-5061-9. [PMC free article] [PubMed] [Cross Ref] Simons R, Gruppen H, Bovee TFH, Verbruggen MA, Vincken JP. Prenylated isoflavonoids from plants as selective estrogen receptor modulators (phytoSERMs) Food Funct. 2012;3:810–827. doi: 10.1039/c2fo10290k. [PubMed] [Cross Ref] Siqueira IR, Lara DR, Silva D, Gaieski FS, Nunes DS, Elisabetsky E. Psychopharmacological properties of Ptychopetalum olacoides Bentham (Olacaceae) Pharmaceut Biol. 1998;36:327–334. doi: 10.1076/phbi.36.5.327.4657. [Cross Ref] Siqueira IR, Fochesatto C, da Silva AL, Nunes DS, Battastini AM, Netto CA, Elisabetsky E. Ptychopetalum olacoides, a traditional Amazonian “nerve tonic”, possesses anticholinesterase activity. Pharmacol Biochem Behav. 2003;75:645–650. doi: 10.1016/S0091-3057(03)00113-8. [PubMed] [Cross Ref] Spartan ’14 v.1.1.2 (2013) Wavefunction, Inc., Irvine, California, USA Spasov AA, Wikman GK, Mandrikov VB, Mironova IA, Neumoin VV. A double-blind, placebo-controlled pilot study of the stimulating and adaptogenic effect of Rhodiola rosea SHR-5 extract on the fatigue of students caused by stress during an examination period with a repeated low-dose regimen. Phytomedicine. 2000;7:85–89. doi: 10.1016/S0944-7113(00)80078-1. [PubMed] [Cross Ref] Sreeja S, Anju VS, Sreeja S. In-vitro estrogenic activities of fenugreek Trigonella foenum graecum seeds. Indian J Med Res. 2010;131:814–819. [PubMed] Stone M, Ibarra A, Roller M, Zangara A, Stevenson E. A pilot investigation into the effect of maca supplementation on physical activity and sexual desire in sportsmen. J Ethnopharmacol. 2009;126:574–576. doi: 10.1016/j.jep.2009.09.012. [PubMed] [Cross Ref] Stonemetz D. A review of the clinical efficacy of evening primrose. Holist Nurs Pract. 2008;22:171–174. doi: 10.1097/01.HNP.0000318026.45527.07. [PubMed] [Cross Ref] Subathra M, Shila S, Devi MA, Panneerselvam C. Emerging role of Centella asiatica in improving age-related neurological antioxidant status. Exptl Gerontol. 2005;40:707–715. doi: 10.1016/j.exger.2005.06.001. [PubMed] [Cross Ref] Tamir S, Eizenberg M, Somjen D, Izrael S, Vaya J. Estrogen-like activity of glabrene and other constituents isolated from licorice root. J Steroid Biochem Mol Biol. 2001;78:291–298. doi: 10.1016/S0960-0760(01)00093-0. [PubMed] [Cross Ref] Tang W, Hioki H, Harada K, Kubo M, Fukuyama Y. Clerodane diterpenoids with NGF-potentiating activity from Ptychopetalum olacoides. J Nat Prod. 2008;71:1760–1763. doi: 10.1021/np8004002. [PubMed] [Cross Ref] Tang W, Kubo M, Harada H, Hioki H, Fukuyama Y. Novel NGF-potentiating diterpenoids from a Brazilian medicinal plant, Ptychopetalum olacoides. Bioorg Med Chem Lett. 2009;19:882–886. doi: 10.1016/j.bmcl.2008.11.100. [PubMed] [Cross Ref] Tang W, Harada K, Kubo M, Hioki H, Fukuyama Y. Eight new clerodane diterpenoids from the bark of Ptychopetalum olacoides. Nat Prod Commun. 2011;6:327–332. [PubMed] Thomsen R, Christensen MH. MolDock: a new technique for high-accuracy molecular docking. J Med Chem. 2006;49:3315–3321. doi: 10.1021/jm051197e. [PubMed] [Cross Ref] van Die MD, Burger HG, Teede HJ, Bone KM. Vitex agnus-castus extracts for female reproductive disorders: a systematic review of clinical trials. Planta Med. 2013;79:562–575. [PubMed] Wang Y, Wang Y, McNeil B, Harvey LM. Maca: an Andean crop with multi-pharmacological functions. Food Res Int. 2007;40:783–792. doi: 10.1016/j.foodres.2007.02.005. [Cross Ref] Waynberg J. Yohimbine vs Muira puama in the treatment of sexual dysfunction. Am J Nat Med. 1994;1:8–9. Waynberg J, Brewer S. Effects of Herbal vX on libido and sexual activity in premenopausal and postmenopausal women. Adv Ther. 2000;17:255–262. doi: 10.1007/BF02853164. [PubMed] [Cross Ref] Wijeweera P, Arnason JT, Koszycki D, Merali Z. Evaluation of anxiolytic properties of Gotukola – (Centella asiatica) extracts and asiaticoside in rat behavioral models. Phytomedicine. 2006;13:668–676. doi: 10.1016/j.phymed.2006.01.011. [PubMed] [Cross Ref] Wu G, Jiang S, Jiang F, Zhu D, Wu H, Jiang S. Steroidal glycosides from Tribulus terrestris. Phytochemistry. 1996;42:1677–1681. doi: 10.1016/0031-9422(96)00182-3. [PubMed] [Cross Ref] Wuttke W, Gorkow C, Seidlová-Wuttke D. Effects of black cohosh (Cimicifuga racemosa) on bone turnover, vaginal mucosa, and various blood parameters in postmenopausal women: a double-blind, placebo-controlled, and conjugated estrogens-controlled study. Menopause. 2006;13:185–196. doi: 10.1097/01.gme.0000174470.44822.57. [PubMed] [Cross Ref] Yan W, Ohtani K, Kasai R, Yamasaki K. Steroidal saponins from fruits of Tribulus terrestris. Phytochemistry. 1996;42:1417–1422. doi: 10.1016/0031-9422(96)00131-8. [PubMed] [Cross Ref] Yuriev E, Ramsland PA. Latest developments in molecular docking: 2010–2011 in review. J Mol Recog. 2013;26:215–239. doi: 10.1002/jmr.2266. [PubMed] [Cross Ref] Yuriev E, Agostino M, Ramsland PA. Challenges and advances in computational docking: 2009 in review. J Mol Recog. 2011;24:149–164. doi: 10.1002/jmr.1077. [PubMed] [Cross Ref] Zakay-Rones Z, Varsano N, Zlotnik M, Manor O, Regev L, Schlesinger M, Mumcuoglu M. Inhibition of several strains of influenza virus in vitro and reduction of symptoms by an elderberry extract (Sambucus nigra L.) during an outbreak of influenza B Panama. J Altern Complement Med. 1995;1:361–369. doi: 10.1089/acm.1995.1.361. [PubMed] [Cross Ref] Zava DT, Dollbaum CM, Blen M. Estrogen and progestin bioactivity of foods, herbs, and spices. Proc Soc Exp Biol Med. 1998;217:369–378. doi: 10.3181/00379727-217-44247. [PubMed] [Cross Ref] Zhang Y, Yu L, Ao M, Jin W. Effect of ethanol extract of Lepidium meyenii Walp. on osteoporosis in ovariectomized rat. J Ethnopharmacol. 2006;105:274–279. doi: 10.1016/j.jep.2005.12.013. [PubMed] [Cross Ref] Zhao J, Pawar RS, Ali Z, Khan IA. Phytochemical investigation of Turnera diffusa. J Nat Prod. 2007;70:289–292. doi: 10.1021/np060253r. [PubMed] [Cross Ref] Zhao J, Dasmahapatra AK, Khan SI, Khan IA. Anti-aromatase activity of the constituents from damiana (Turnera diffusa) J Ethnopharmacol. 2008;120:387–393. doi: 10.1016/j.jep.2008.09.016. [PubMed] [Cross Ref] Zheng BL, He K, Kim CH, Rogers L, 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. doi: 10.1016/S0090-4295(99)00549-X. [PubMed] [Cross Ref] Articles from In Silico Pharmacology are provided here courtesy of Springer