Investigation of Estrogen Receptor (ESR1) for Breast Cancer from Traditional Chinese Medicine
Recently, an important topic of breast cancer had been published in 2013. In this report, estrogen receptor (ESR1) had defined the relation of hormone-cause breast cancer. The screening of traditional Chinese medicine (TCM) database has found the molecular compounds by simulating molecular docking and molecular dynamics to regulate ESR1. S-Allylmercaptocysteine and 5-hydroxy-L-tryptophan are selected according to the highest docking score than that of other TCM compounds and Raloxifene (control). The simulation from molecular dynamics is helpful in analyzing and detecting the protein-ligand interactions. After a comparing the control and the Apo form, then based on the docking poses, hydrophobic interactions, hydrogen bond and structure variations, this research postulates that S-allylmercaptocysteine may be more appropriate than other compounds for protein-ligand interaction.
There is an important topic of breast cancer that had been published in 2013. In this report, the estrogen receptor (ESR1) had defined the relation of breast cancer .
The breast cancer is a common cause of death in women and may originate either hormonally [2–5] from hereditary factors [6–9] and some other reasons . Although there are lots of methods to find breast cancer based on accurate prediction, the surgery treatment causes patient low emotion from losing breast. Thus the prevention and treatment early medicinal treatment is widely accepted.
The ESR1 is a ligand-activated transcription factor that mediates the biological effects of the steroid hormone estrogen [11, 12]. From the inhibition of ESR1, the hormone would be regulated then decreased the risk for hormone breast cancer.
The Raloxifene is an FDA approval drug for reduction in the risk of breast cancer (http://www.cancer.gov/cancertopics/druginfo/fda-raloxifene-hydrochloride). The reference  had reported that the Raloxifene could target ESR1. For the above reason, the Raloxifene was chosen as the control for this investigation.
The computer-aided drug design (CADD) is an in silico simulation technique used in the screening of compounds based on the structure and biological activity. The CADD has the advantages of both greater speed and lower cost than the traditional drug design. The two major application areas of CADD are structure-based drug design and ligand-based drug design [14–17]. CADD is used to investigate drug design application centered on structure-based drug design and molecular dynamics.
The personalized medicine and biomedicine  are a novel department knowledge which could analyze the mutation , and the cause for special disease . The traditional Chinese medicine (TCM) is defined a model case in this department. TCM has an important diagnosis culture in Asia, especially in China, Taiwan, Korea, and Japan. The TCM [email protected] (http://tcm.cmu.edu.tw/)  is the largest traditional Chinese medicine database in the world. Since it was established in 2011, there have been successful discoveries of novel lead compounds from the TCM [email protected] application for cancer treatment [22–25], memory disease prevention , pain relief , and antivirals [27–31]. Today, the screening of TCM compounds is possible from the application system of the website  and the cloud computing platform .
In this research, we screen a candidate compound against breast cancer from the TCM [email protected] The computational techniques of docking screening are used to select ligands prior to applying molecular simulation by molecular dynamics (MD) to investigate the variations in protein-ligand interactions that may contribute to the evaluation of the effects on ESR1 inhibition.
2. Materials and Methods
2.1. Data Set
A total of 61,000 TCM compounds were downloaded from the TCM database (http://tcm.cmu.edu.tw/). The human ESR1 (PDB ID: 1GWQ) crystal structure was obtained from RCSB Protein Data Bank [1, 13]. The Accelrys Discovery Studio 2.5 (DS 2.5) was used as the molecular simulations platform.
2.2. Disorder Protein Detection
The disorder region in protein plays an important role in drug design; thus we take the sequence to predict the disorder region by the Database of Protein Disorder (DisProt: http://www.disprot.org/) . The prediction decides the character of protein structure; then taking the comparison with the docking site could evaluate the efficacy of the drug during protein-ligand interaction.
2.3. Molecular Docking
The LigandFit module , a receptor-rigid docking algorithm program in Discovery Studio 2.5 (DS 2.5), was used for docking simulations of Raloxifene and TCM compounds to ESR1 in the CHARMm force field . The docking site of ESR1 was designed on the basis of the research [1, 13]. Through docking simulation, the top two compounds with the highest docking scores of the TCM compounds were selected to make the analysis of the hydrophobic interactions by Ligplus [37, 38].
2.4. Molecular Dynamics Simulation
The ligands of candidate complex must be reprepared before applying MD simulation by using SwissParam (http://swissparam.ch/)  based on the reference force field  of GROMACS 4.5.5 . The ESR1, with ligands, was placed in a simulation box in an appropriate buffer or other solution at a minimum distance of 1.2 Å from the complex. The solution for simulation was based on the TIP3P water model in which sodium and chloride ions were added to neutralize complex charges. Based on the Steepest Descent method for 5,000 steps to minimize the complex, the structure with the lowest energy was transferred to MD simulation. The electrostatic interactions were calculated on the basis of the particle-mesh Ewald (PME) method . The calculation with each time step was 2 fs and the numbers of steps were 5,000,000 times then the total simulation time of MD was 10,000 ps. The equilibration under the 100 ps constant temperature (PER ensemble) was based on the Berendsen weak thermal coupling method. The protocols in Gromacs used the MD data to calculate the MD trajectories, RMSD, energy variations, and eigenvector after MD.
3. Results and Discussion
3.1. The Detection of Disorder Protein
The disorder protein is defined as unstructured protein which makes the compounds dock to protein and stabilize the complex with difficultly while the docking site is a disorder region. The cited references [17, 43] indicate that the disorder region may have lower side effect than the widespread domain. Thus the disorder region can be defined as a challenge for drug design. In the prediction, the residue with the disposition greater than 0.5 is defined as disorder region (Figure 1). In this result, all the important amino acids of ESR1 are less than the threshold; thus the disorder protein has a weaker effect on docking and simulation.
3.2. Molecular Docking
Ranking the result of molecular docking by docking score, the two top TCM compounds, and the control, were selected (Table 1). These TCM compounds are S-allylmercaptocysteine and 5-hydroxy-L-tryptophan, extract from the TCM herbs Allium sativum and Mucuna pruriens, respectively. The top compound, S-allylmercaptocysteine, is defined as an hepatoprotective and anticancer compound [44–54] and the herb Allium sativum has antimicrobial properties [55–57], ameliorates tamoxifen-induced liver injury , and prevents cancer [45, 59–62]. The second ranked herb, Mucuna pruriens, has been identified as being able to reduce oxidation and prevent Parkinson’s disease [63, 64]. As mentioned above, the top-ranked compound could prevent or treat cancer. The second compound with antioxidation might ease symptoms of cancer. For the above reasons, we suggest that the selected compounds can have influence on ESR1.
The structure of control and the candidate compounds was selected after screening from the TCM database (Figure 2). The docking poses presents the ligands had interactions with different critical amino-acids in the protein. Figure 3 indicates the selected compounds could target and interact with amino acids around docking site.
The hydrophobic interaction can be analyzed by Ligplus (Figure 4). This result shows that the amino acids Glu353, Leu391, Arg394, and Phe404 can interact with all ligands through hydrophobic interactions or hydrogen bonds, indicating that these amino acids might be important in ligand-protein binding situations.
3.3. Molecular Dynamics Simulation
The RMSD and energy variation of a complex during MD simulation were recorded (Figure 5). The total energy is in the range between −744 and kJ/mol and tends to kJ/mol. From this figure, if both the amplitude of complex RMSD and energy are less; then the simulation may become balanced. S-allylmercaptocysteine has larger variation than other ligands but the complex RMSD of S-allylmercaptocysteine is the lowest. From this situation, we suggest S-allylmercaptocysteine can still interact, and this interaction makes the complex more stable.
The RMSF focus on each residue was analyzed, and on the variation of the whole protein, including with the ligand interaction (Figure 6). In this result, it can be seen that the regions of protein-ligand interaction are similar.
The reference-identified eigenvector was used to represent the protein variation . The first two eigenvectors were selected based on PCA (principal component analysis) calculation, and become the - and -axes. The comparison with apo (unbound protein) could find protein variation of first main character of protein (Figure 7). The upper subunit in these figures are the first eigenvector diffusion between apo and complex. The following is a matrix established from first two eigenvectors. After the comparison, we find that complex with S-allylmercaptocysteine is similar to apo then different from other compounds. This situation may not mean an absence variation, but the position of variation might be smaller.
The clustering is a result of the division of data into several groups based on RMSD variation; thus data in the same group will have the similar structure (position and composition) (Figure 8). In this result, S-allylmercaptocysteine has the least group and the largest one is in last. This situation means the complex with S-allylmercaptocysteine will tend to balance quickly. The largest group of apo and 5-hydroxy-L-tryptophan is not in last and it might present this simulation for apo and 5-hydroxy-L-tryptophan is not enough. Thus S-allylmercaptocysteine and control, Raloxifene, might have better effect on ESR1.
After the structure variation discussion is based on eigenvector and clustering, we should take focus on the structure variation during protein-ligand interaction (Figures 9 to 11). In Figure 9(a), there is high percentage (100%) of H bond occupancy in Glu353 which indicates that Glu353 may have a function in the inhibition of ESR1 from Raloxifene interaction. Figure 9(b) shows the variation between MD 0 ns and 10 ns which present the position and composition variation as control drug inhibit ESR1. As the variation of apo between 0 ns and 10 ns is smaller position variation, the variation of control might become a sample to detect the efficacy of compounds from the structure variation of simulation.
S-Allylmercaptocysteine has high H bond occupancy in both Glu353 and Arg394, with the variation in 2 and 5 being more variable than the control (Figure 10). This result might present that not only S-allylmercaptocysteine the efficacy as control, but also that the force might be stronger than control.
The 5-hydroxy-L-tryptophan complex interactions were recorded (Figure 11). Besides Glu353 and Arg394, Leu346 has been recorded in Figure 11(a). The structure variation of 5-hydroxy-L-tryptophan is more similar to S-allylmercaptocysteine than control. From this result, we suggest 5-hydroxy-L-tryptophan might also have the function for the inhibition of ESR1.
From these variations found, we suggest Glu353 might be important in inhibition and Arg394 might make the force stronger based on the H bond and structure variation.
The pathway definition is according to the calculation of caver 3.0 to determine the inter-path protein path during MD simulation . The pathway could help to determine the ligand moving and the pole provided from protein after interaction (Figure 12). In these results, we could find a path through ESR1. We suggest the ligands inhibit ESR1 while the interaction is in the protein.
Based on the above discussion, we can find that the top two TCM compounds S-allylmercaptocysteine and 5-hydroxy-L-tryptophan can have effect on ESR1 against breast cancer. Glu353 might have important role in inhibition based on high H bond occupancy in MD. Finally, according to the discussion from docking, interaction, and variation, we suggest that S-allylmercaptocysteine might be the best compound to inhibit ESR1 against breast cancer, even better than the control.
Conflict of Interests
The authors confirm that this paper’s content has no conflict of interests.
Tzu-Chieh Hung, Wen-Yuan Lee, and Kuen-Bao Chen contributed equally.
The research was supported by Grants from the National Science Council of Taiwan (NSC102-2325-B039-001, NSC102-2221-E-468-027-), Asia University (ASIA100-CMU-2, ASIA101-CMU-2, and 102-Asia-07), and China Medical University Hospital (DMR-103-058, DMR-103-001, and DMR-103-096). This study is also supported in part by Taiwan Department of Health Clinical Trial and Research Center of Excellence (DOH102-TD-B-111-004), Taiwan Department of Health Cancer Research Center of Excellence (MOHW103-TD-B-111-03), and CMU under the Aim for Top University Plan of the Ministry of Education, Taiwan. Finally, our gratitude goes to Dr. Tim Williams, Asia University.
W. Toy, Y. Shen, H. Won et al., “ESR1 ligand-binding domain mutations in hormone-resistant breast cancer,” Nature Genetics, vol. 45, no. 12, pp. 1439–1445, 2013.View at: Google Scholar
J. D. Yager and N. E. Davidson, “Estrogen carcinogenesis in breast cancer,” The New England Journal of Medicine, vol. 354, no. 3, pp. 228–282, 2006.View at: Publisher Site | Google Scholar
B. Park, K. Kim, M. Heo et al., “The changes of estrogen receptor-β variants expression in breast carcinogenesis: decrease of estrogen receptor-β2 expression is the key event in breast cancer development,” Journal of Surgical Oncology, vol. 93, no. 6, pp. 504–510, 2006.View at: Publisher Site | Google Scholar
D. Germain, “Estrogen carcinogenesis in breast cancer,” Endocrinology and Metabolism Clinics of North America, vol. 40, no. 3, pp. 473–484, 2011.View at: Publisher Site | Google Scholar
H. Al Kadri, S. Hassan, H. M. Al-Fozan, and A. Hajeer, “Hormone therapy for endometriosis and surgical menopause,” The Cochrane Database of Systematic Reviews, no. 1, Article ID CD005997, 2009.View at: Google Scholar
D. Mutch, L. Denny, and M. Quinn, “Hereditary gynecologic cancers,” International Journal of Gynaecology and Obstetrics, vol. 124, no. 3, pp. 189–192, 2014.View at: Google Scholar
A. Floquet, E. Stoeckle, S. Croce et al., “Hereditary ovarian carcinomas: clinico-biological features and treatment,” Bulletin du Cancer, vol. 101, no. 2, pp. 167–174, 2014.View at: Google Scholar
M. M. Koeneman, A. J. Kruse, S. J. Sep et al., “A family history questionnaire improves detection of women at risk for hereditary gynecologic cancer: a pilot study,” Familial Cancer, 2014.View at: Publisher Site | Google Scholar
E. de Geus, C. M. Aalfs, M. G. Verdam, H. C. de Haes, and E. M. Smets, “Informing relatives about their hereditary or familial cancer risk: study protocol for a randomized controlled trial,” Trials, vol. 15, no. 1, p. 86, 2014.View at: Google Scholar
V. Beral, D. Bull, R. Doll, R. Peto, and G. Reeves, “Breast cancer and breastfeeding: collaborative reanalysis of individual data from 47 epidemiological studies in 30 countries, including 50 302 women with breast cancer and 96 973 women without the disease,” The Lancet, vol. 360, no. 9328, pp. 187–195, 2002.View at: Publisher Site | Google Scholar
C. M. Klinge, “Estrogen receptor interaction with estrogen response elements,” Nucleic Acids Research, vol. 29, no. 14, pp. 2905–2919, 2001.View at: Google Scholar
D. Lopez, M. D. Sanchez, W. Shea-Eaton, and M. P. Mclean, “Estrogen activates the high-density lipoprotein receptor gene via binding to estrogen response elements and interaction with sterol regulatory element binding protein-1A,” Endocrinology, vol. 143, no. 6, pp. 2155–2168, 2002.View at: Publisher Site | Google Scholar
A. Wärnmark, E. Treuter, J. Gustafsson, R. E. Hubbard, A. M. Brzozowski, and A. C. W. Pike, “Interaction of transcriptional intermediary factor 2 nuclear receptor box peptides with the coactivator binding site of estrogen receptor α,” The Journal of Biological Chemistry, vol. 277, no. 24, pp. 21862–21868, 2002.View at: Publisher Site | Google Scholar
H. Huang, H. W. Yu, C. Chen et al., “Current developments of computer-aided drug design,” Journal of the Taiwan Institute of Chemical Engineers, vol. 41, no. 6, pp. 623–635, 2010.View at: Publisher Site | Google Scholar
W. I. Tou, S. Chang, C. Lee, and C. Y. Chen, “Drug design for neuropathic pain regulation from traditional Chinese medicine,” Scientific Reports, vol. 3, p. 844, 2013.View at: Google Scholar
C. Y. Chen, “A novel integrated framework and improved methodology of computer-aided drug design,” Current Topics in Medicinal Chemistry, vol. 13, no. 9, pp. 965–988, 2013.View at: Publisher Site | Google Scholar
C. Y. Chen and W. I. Tou, “How to design a drug for the disordered proteins?” Drug Discovery Today, vol. 18, pp. 910–915, 2013.View at: Google Scholar
W. Liao and F. Tsai, “Personalized medicine: a paradigm shift in healthcare,” BioMedicine, vol. 3, no. 2, pp. 66–72, 2013.View at: Publisher Site | Google Scholar
C.-C. Lee, C.-H. Tsai, L. Wan et al., “Increased incidence of Parkinsonism among Chinese with β-glucosidase mutation in central Taiwan,” BioMedicine, vol. 3, no. 2, pp. 92–94, 2013.View at: Google Scholar
I. C. Chou, W.-D. Lin, C.-H. Wang et al., “Möbius syndrome in a male with XX/XY mosaicism,” BioMedicine, vol. 3, no. 2, pp. 102–104, 2013.View at: Google Scholar
C. Y. Chen, “TCM [email protected]: the world's largest traditional Chinese medicine database for drug screening in silico,” PLoS ONE, vol. 6, no. 1, Article ID e15939, 2011.View at: Publisher Site | Google Scholar
H. Huang, K. Lee, H. W. Yu et al., “Structure-based and ligand-based drug design for HER 2 receptor,” Journal of Biomolecular Structure and Dynamics, vol. 28, no. 1, pp. 23–37, 2010.View at: Google Scholar
W. I. Tou and C. Y. Chen, “Traditional Chinese medicine as dual guardians against hypertension and cancer?” Journal of Biomolecular Structure and Dynamics, vol. 30, no. 3, pp. 299–317, 2012.View at: Publisher Site | Google Scholar
S. Yang, S. Chang, and C. Y. Chen, “Identifying HER2 inhibitors from natural products database,” PLoS ONE, vol. 6, no. 12, Article ID e28793, 2011.View at: Publisher Site | Google Scholar
C. Chen and C. Y. Chen, “Insights into designing the dual-targeted HER2/HSP90 inhibitors,” Journal of Molecular Graphics and Modelling, vol. 29, no. 1, pp. 21–31, 2010.View at: Publisher Site | Google Scholar
H.-J. Huang, C.-C. Lee, and C. Y.-C. Chen, “Pharmacological chaperone design for reducing risk factor of Parkinson's disease from traditional chinese medicine,” Evidence-Based Complementary and Alternative Medicine, vol. 2014, Article ID 830490, 12 pages, 2014.View at: Publisher Site | Google Scholar
S. Chang, H. Huang, and C. Y. Chen, “High performance screening, structural and molecular dynamics analysis to identify H1 inhibitors from TCM [email protected],” Molecular BioSystems, vol. 7, no. 12, pp. 3366–3374, 2011.View at: Publisher Site | Google Scholar
C. Lin, T. Chang, M. Sun et al., “Potent inhibitor design against H1N1 swine influenza: structure-based and molecular dynamics analysis for M2 inhibitors from traditional Chinese medicine database,” Journal of Biomolecular Structure and Dynamics, vol. 28, no. 4, pp. 471–482, 2011.View at: Google Scholar
S. Chang, H. Huang, and C. Y. Chen, “Two birds with one stone? Possible dual-targeting H1N1 inhibitors from traditional Chinese medicine,” PLoS Computational Biology, vol. 7, no. 12, Article ID e1002315, 2011.View at: Publisher Site | Google Scholar
C. Chen, Y. Chang, D. Bau et al., “Ligand-based dual target drug design for H1N1: swine flu—a preliminary first study,” Journal of Biomolecular Structure and Dynamics, vol. 27, no. 2, pp. 171–178, 2009.View at: Google Scholar
H. J. Huang, Y. R. Jian, and C. Y. Chen, “Traditional Chinese medicine application in HIV: an in silico study,” Journal of Biomolecular Structure and Dynamics, vol. 32, no. 1, pp. 1–12, 2014.View at: Google Scholar
T. Tsai, K. Chang, and C. Y. Chen, “IScreen: world's first cloud-computing web server for virtual screening and de novo drug design based on TCM [email protected],” Journal of Computer-Aided Molecular Design, vol. 25, no. 6, pp. 525–531, 2011.View at: Publisher Site | Google Scholar
K. Chang, T. Tsai, K. Chen et al., “iSMART: an integrated cloud computing web server for traditional Chinese medicine for online virtual screening, de novo evolution and drug design,” Journal of Biomolecular Structure and Dynamics, vol. 29, no. 1, pp. 243–250, 2011.View at: Google Scholar
B. Xue, R. L. Dunbrack, R. W. Williams, A. K. Dunker, and V. N. Uversky, “PONDR-FIT: a meta-predictor of intrinsically disordered amino acids,” Biochimica et Biophysica Acta—Proteins and Proteomics, vol. 1804, no. 4, pp. 996–1010, 2010.View at: Publisher Site | Google Scholar
C. M. Venkatachalam, X. Jiang, T. Oldfield, and M. Waldman, “LigandFit: a novel method for the shape-directed rapid docking of ligands to protein active sites,” Journal of Molecular Graphics and Modelling, vol. 21, no. 4, pp. 289–307, 2003.View at: Publisher Site | Google Scholar
B. R. Brooks, C. L. Brooks III, A. D. Mackerell Jr. et al., “CHARMM: the biomolecular simulation program,” Journal of Computational Chemistry, vol. 30, no. 10, pp. 1545–1614, 2009.View at: Publisher Site | Google Scholar
R. A. Laskowski and M. B. Swindells, “LigPlot+: multiple ligand-protein interaction diagrams for drug discovery,” Journal of Chemical Information and Modeling, vol. 51, no. 10, pp. 2778–2786, 2011.View at: Publisher Site | Google Scholar
A. C. Wallace, R. A. Laskowski, and J. M. Thornton, “LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions,” Protein Engineering, vol. 8, no. 2, pp. 127–134, 1995.View at: Google Scholar
V. Zoete, M. A. Cuendet, A. Grosdidier, and O. Michielin, “SwissParam: a fast force field generation tool for small organic molecules,” Journal of Computational Chemistry, vol. 32, no. 11, pp. 2359–2368, 2011.View at: Publisher Site | Google Scholar
U. D. Priyakumar and A. D. MacKerell, “Comparison of the CHARMM27, AMBER4. 1 and BMS nucleic acid force fields via free energy calculations of base flipping,” Abstracts of Papers of the American Chemical Society, vol. 230, pp. U1391–U1392, 2005.View at: Google Scholar
B. Hess, C. Kutzner, D. van der Spoel, and E. Lindahl, “GRGMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation,” Journal of Chemical Theory and Computation, vol. 4, no. 3, pp. 435–447, 2008.View at: Publisher Site | Google Scholar
T. A. Darden and L. G. Pedersen, “Molecular modeling: an experimental tool,” Environmental Health Perspectives, vol. 101, no. 5, pp. 410–412, 1993.View at: Google Scholar
W. I. Tou and C. Y. Chen, “May disordered protein cause serious drug side effect?” Drug Discovery Today, vol. 19, no. 4, pp. 367–372, 2014.View at: Google Scholar
J.-Y. Yan, F.-M. Tian, W.-N. Hu, J.-H. Zhang, H.-F. Cai, and N. Li, “Apoptosis of human gastric cancer cells line SGC 7901 induced by garlic-derived compound S-allylmercaptocysteine (SAMC),” European Review for Medical and Pharmacological Sciences, vol. 17, no. 6, pp. 745–751, 2013.View at: Google Scholar
W. Zhang, H. Xiao, and K. L. Parkin, “Apoptosis in MCF-7 breast cancer cells induced by S-alkenylmercaptocysteine (CySSR) species derived from Allium tissues in combination with sodium selenite,” Food and Chemical Toxicology C, vol. 68, pp. 1–10, 2014.View at: Google Scholar
D. Liang, Y. Qin, W. Zhao et al., “S-allylmercaptocysteine effectively inhibits the proliferation of colorectal cancer cells under in vitro and in vivo conditions,” Cancer Letters, vol. 310, no. 1, pp. 69–76, 2011.View at: Publisher Site | Google Scholar
H. Hu, X. P. Zhang, Y. L. Wang et al., “Identification of a novel function of Id-1 in mediating the anticancer responses of SAMC, a water-soluble garlic derivative, in human bladder cancer cells,” Molecular Medicine Reports, vol. 4, no. 1, pp. 9–16, 2011.View at: Publisher Site | Google Scholar
Y. Lee, “Induction of apoptosis by S-allylmercapto-L-cysteine, a biotransformed garlic derivative, on a human gastric cancer cell line,” International Journal of Molecular Medicine, vol. 21, no. 6, pp. 765–770, 2008.View at: Google Scholar
E. W. Howard, D. T. Lee, T. C. Yung, W. C. Chee, X. Wang, and C. W. Yong, “Evidence of a novel docetaxel sensitizer, garlic-derived S-allylmercaptocysteine, as a treatment option for hormone refractory prostate cancer,” International Journal of Cancer, vol. 122, no. 9, pp. 1941–1948, 2008.View at: Publisher Site | Google Scholar
J. O. Ban, D. Y. Yuk, K. S. Woo et al., “Inhibition of cell growth and induction of apoptosis via inactivation of NF-κB by a sulfurcompound isolated from garlic in human colon cancer cells,” Journal of Pharmacological Sciences, vol. 104, no. 4, pp. 374–383, 2007.View at: Publisher Site | Google Scholar
J. Xiao, Y. P. Ching, E. C. Liong, A. A. Nanji, M. L. Fung, and G. L. Tipoe, “Garlic-derived S-allylmercaptocysteine is a hepato-protective agent in non-alcoholic fatty liver disease in vivo animal model,” European Journal of Nutrition, vol. 52, no. 1, pp. 179–191, 2013.View at: Publisher Site | Google Scholar
J. Xiao, R. Guo, M. Fung et al., “Garlic-derived S-allylmercaptocysteine ameliorates nonalcoholic fatty liver disease in a rat model through inhibition of apoptosis and enhancing autophagy,” Evidence-Based Complementary and Alternative Medicine, vol. 2013, Article ID 642920, 11 pages, 2013.View at: Publisher Site | Google Scholar
J. Xiao, E. C. Liong, M. Ling, Y. Ching, M. Fung, and G. L. Tipoe, “S-allylmercaptocysteine reduces carbon tetrachloride-induced hepatic oxidative stress and necroinflammation via nuclear factor kappa B-dependent pathways in mice,” European Journal of Nutrition, vol. 51, no. 3, pp. 323–333, 2012.View at: Publisher Site | Google Scholar
Y. Lee, H. Kim, J. Lee, and K. Kim, “Anticancer activity of S-allylmercapto-L-cysteine on implanted tumor of human gastric cancer cell,” Biological and Pharmaceutical Bulletin, vol. 34, no. 5, pp. 677–681, 2011.View at: Publisher Site | Google Scholar
S. A. Seddiek, M. M. El-Shorbagy, H. F. Khater, and A. M. Ali, “The antitrichomonal efficacy of garlic and metronidazole against Trichomonas gallinae infecting domestic pigeons,” Parasitology Research, vol. 113, no. 4, pp. 1319–1329, 2014.View at: Google Scholar
Y. Liu, M. Song, T. M. Che et al., “Dietary plant extracts modulate gene expression profiles in ileal mucosa of weaned pigs after an Escherichia coli infection,” Journal of Animal Science, vol. 92, no. 5, pp. 2050–2062, 2014.View at: Google Scholar
R. Snowden, H. Harrington, K. Morrill et al., “A comparison of the anti-Staphylococcus aureus activity of extracts from commonly used medicinal plants,” Journal of Alternative and Complementary Medicine, vol. 20, no. 5, pp. 375–382, 2014.View at: Google Scholar
G. M. Suddek, “Allicin enhances chemotherapeutic response and ameliorates tamoxifen-induced liver injury in experimental animals,” Pharmaceutical Biology, 2014.View at: Publisher Site | Google Scholar
M. Jo, H. M. Yun, K. R. Park et al., “Anti-cancer effect of thiacremonone through down regulation of peroxiredoxin 6,” PloS ONE, vol. 9, no. 3, Article ID e91508, 2014.View at: Google Scholar
M. Dong, G. Yang, H. Liu et al., “Aged black garlic extract inhibits HT29 colon cancer cell growth via the PI3K/Akt signaling pathway,” Biomedical Reports, vol. 2, no. 2, pp. 250–254, 2014.View at: Google Scholar
G. Schafer and C. H. Kaschula, “The immunomodulation and anti-inflammatory effects of garlic organosulfur compounds in cancer chemoprevention,” Anti-Cancer Agents in Medicinal Chemistry, vol. 14, no. 2, pp. 233–240, 2014.View at: Google Scholar
E. R. Hahm and S. V. Singh, “Diallyl trisulfide inhibits estrogen receptor-alpha activity in human breast cancer cells,” Breast Cancer Research and Treatment, vol. 144, no. 1, pp. 47–57, 2014.View at: Google Scholar
S. K. Yadav, J. Prakash, S. Chouhan, and S. P. Singh, “Mucuna pruriens seed extract reduces oxidative stress in nigrostriatal tissue and improves neurobehavioral activity in paraquat-induced Parkinsonian mouse model,” Neurochemistry International, vol. 62, no. 8, pp. 1039–1047, 2013.View at: Publisher Site | Google Scholar
S. Ovallath and P. Deepa, “The history of parkinsonism: descriptions in ancient Indian medical literature,” Movement Disorders, vol. 28, no. 5, pp. 566–568, 2013.View at: Publisher Site | Google Scholar
J. H. Peters and B. L. de Groot, “Ubiquitin dynamics in complexes reveal molecular recognition mechanisms beyond induced fit and conformational selection,” PLoS Computational Biology, vol. 8, no. 10, Article ID e1002704, 2012.View at: Publisher Site | Google Scholar
E. Chovancova, A. Pavelka, P. Benes et al., “CAVER 3. 0: a tool for the analysis of transport pathways in dynamic protein structures,” PLoS Computational Biology, vol. 8, no. 10, Article ID e1002708, 2012.View at: Publisher Site | Google Scholar