Home Join Contact
 

Research Article
Open Access
Geraniol and Limonene Interaction with 3-hydroxy-3- methylglutaryl-CoA (HMG-CoA) Reductase for their Role as Cancer Chemo-preventive Agents
Madhumita Pattanayak1,2, P K Seth3, Suchi Smita4, Shailendra K Gupta1,5*
1Indian Institute of Toxicology Research, Lucknow (CSIR) India
2Centre for Cellular and Molecular Biology, Hyderabad (CSIR) India
3Biotech Park, Lucknow India
4Department of Biological Sciences, University of Rostock, Rostock Germnay
5System Biology & Bioinformatics Group, University of Rostock, Rostock Germnay
*Corresponding author:
Dr. Shailendra K Gupta,
System Biology & Bioinformatics Group,
University of Rostock, 18051, Rostock, Germany,
Tel       : +49 381 4987578,
Fax     : +49 381 4987572,
E-mail: shailendra.gupta@uni-rostock.de
Received September 26, 2009; Accepted November 23, 2009; Published November 24, 2009
Citation: Pattanayak M, Seth PK, Smita S, Gupta SK (2009) Geraniol and Limonene Interaction with 3-hydroxy-3-methylglutaryl-CoA (HMGCoA) Reductase for their Role as Cancer Chemo-preventive Agents. J Proteomics Bioinform 2: 466-474. doi:10.4172/jpb.1000107
 
Copyright: © 2009 Pattanayak M, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
 
Abstract
Recent studies have shown that monoterpenes exhibit antitumor activities and suggest that these compounds are a new class of cancer chemo-preventive agents. Limonene, a main constituent of orange and citrus peel oils has been reported to exert antitumor activity against mammary gland, lung, liver, stomach and skin cancers in rodents whereas, geraniol, a principal constituent of Geranium and Ocimum inhibits the growth of human colon cancer cells. Prenylation of proteins is essential for progression of cells into the S phase and involves post-translational covalent attachment of a lipophilic farnesyl or geranylgeranyl isoprenoid group to numerous proteins. Suppression of prenylation of proteins leads to inhibition of DNA synthesis. Further, epidemiologic evidences suggest that suppression of hydrophilic 3-hydroxy-3-methylglutaryl-CoA (HMG-CoA) reductase, a key enzyme of mevalonate biosynthesis, leads to reduction of the mevalonate pool and
thus limits protein isoprenylation.

Geraniol and limonene inhibit the activity of HMG-CoA reductase subsequently reducing the possibility of cancer growth. In the present work, we analyzed binding affinity of limonene and geraniol with HMG-CoA and explored mechanism of interaction using in silico approaches. The binding positions were verified according to their energy, PMF (Potential of Mean Force) value, PLP (Piecewise Linear Potential) value and Ligand Internal energy. It was found that limonene had greater binding affinity with the receptor suggesting better antitumor agent in comparison to geraniol.

Keywords
Limonene; Geraniol; HMG CoA reductase; Cancer chemo-preventive agents; docking

Introduction
Essential oils are highly concentrated volatile aromatic essences of plants. They are mainstay of aromatherapy but are also used in flavoring, perfumes and even as solvents. Terpenes, aldehydes, esters, ketones, alcohol, phenol and oxides are major components of essential oils. Monoterpenes function physiologically as chemo-attractants or chemo-repellents, and they are largely responsible for the distinctive fragrance of many plants (McGarvey et al., 1995). Significant scientific evidences
are there to suggest that nutritive and non-nutritive plant-based dietary factors can inhibit the process of carcinogenesis effectively (Singletery, 2000). Monoterpenes are non nutritive dietary
components found in the essential oils of plants having antitumor activity, exhibiting not only the ability to prevent the formation or progression of cancer, but also regress existing malignant tumors (Crowell, 1999). The human exposure to monoterpenes through the diet or environment is widespread.

Major monoterpenes includes limonene, pinenene, menthol, geraniol, camphene, sabinene, cadinine. Monoterpenes consist of two isoprene units with the molecular formula C10H16. Monoterpenes may be linear (acyclic) or contain rings. These 10 carbon isoprenoids are derived from the mevalonate pathway in plants but are not produced by mammals, fungi or other species (Loza-Tavera, 1999). Citrus fruit, orange and peppermint are the main sources of d-limonene i.e. p-mentha-1,8-diene (Kodama et al., 1977). d-limonene (Figure 1) is a prevalent flavoring agent and because of its pleasant citrus fragrance, it is commonly added to cosmetics, soa psand other cleaning products. It is a cyclic monoterpene and formed by the cyclization of geranylpyrophosphate in a reaction catalyzed by limonene synthase (Alonso et al., 1992; Kjonaas et al., 1983). Limonene has well-established chemo-preventive activity against many cancer types. Limonene has been shown to inhibit the development of spontaneous neoplasms in mice at the dose of 1200 mg/kg orally (National Toxicology Program, 1990). Dietary limonene also reduces the incidence of spontaneous lymphomas in p53-/- mice (Salim et al., 2003). When administered either in pure form or as orange peel oil (95% d-limonene), limonene inhibits the development of chemically induced rodent mammary (Asamoto et al., 2002), skin (Elegbede et al., 1986), liver (Lu et al., 2004), lung and stomach (Raphael and Kuttan, 2003) cancers. In rat mammary carcinogenesis models, the chemo-preventive effects of limonene are evident during the initiation phase of 7-12- dimethylbenz[a]anthracene (DMBA)-induced cancer (Elson et al., 1988) and during the promotion phase of both DMBA- and nitrosomethylurea (NMU)-induced cancers (Chander et al., 1994). Dietary limonene also inhibits the development of ras oncogene–induced mammary carcinomas in rats (Gould et al., 1994). Development of azoxymethane-induced aberrantcrypt f oci in the colon of rats was significantly reduced when they were given 0.5% limonene in the drinking water (Kawamori et al., 1996).

Figure1: Chemical structure of d-limonene (p-mentha-1,8-diene ).

M.F: C10H16, CAS No: 5989-27-5.

Main sources of geraniol i.e. trans-3,7-Dimethyl-2,6-octadien- 1-ol (Figure 2) are bergamot, carrot, coriander, lavender, lemon, lime, nutmeg, orange, rose, blueberry, basil and blackberry. It is mainly used in perfumery and flavouring industries. Geraniol synthase is involved in the terpene biosynthetic pathway converting geranyl diphosphate to geraniol (Iijima et al., 2004). Geraniol, an acyclic monoterpene, has antitumor activity against mu rine leukemia, hepatoma and melanoma cells in vivo when administered before and after tumor cell transplantation. It has antiproliferative effects on hepatoma and melanoma cell growth (Polo and de Bravo, 2006). Geraniol (400 μM) caused a 70% inhibition of cell growth in human colon cancer cell lines. Geraniol has shown anti-tumoral efficacy on TC-118 human tumors transplanted in Swiss nu/nu mice. Geraniol (150 μM) has been identified to reduce thymidylate synthase and thymidine kinase expression in cancer cells. In nude mice, the combined administration of 5-fluorouracil (20 mg/kg) and geraniol (150 mg/kg) caused a 53% reduction of the tumor volume, whereas a 26% reduction was obtained with geraniol alone (Carnesecchi et al., 2004).

Figure2: Chemical structure of geraniol (trans-3,7-Dimethyl- 2,6-octadien-1-ol).

M.F: C10H18O, CAS No: 106-24-1.

HMG-CoA reductase (Figure 3) is a polytopic, transmembrane protein that catalyzes a key step in the mevalonate pathway (conversion of HMG-CoA to mevalonate). Mevalonate is necessary for cell growth (Swanson and Hohl, 2006) and is involved in the synthesis of sterols, isoprenoids and other lipids. HMG-CoA reductase is the rate-limiting step in cholesterol synthesis and represents the sole major drug target for contemporary cholesterol- lowering drugs (Genser et al., 2008). HMG-CoA reductase is also an important developmental enzyme. Limonene and geraniol suppress HMG-CoA reductase synthesis in mammalian cells by decreasing the translational efficiency of HMG-CoA reductase transcripts (Peffley and Gayen, 2003) and thus reduce mevalonate production. Terpenoids reduce cancer formation by the simple reduction of synthesis of chlolesterol and ubiquinone and other cholesterol derivatives that are necessary for the cell proliferation. It is speculated that mevalonate is probably involved in the post-translational modification of proteins involved in cell turnover. The reduction of the mevalonate pool limits protein isoprenylation, which involves the post-translational covalent attachment of a lipophilic farnesyl or geranylgeranyl isoprenoid group to numerous proteins (Clarke, 1992).

Figure3: Rotatable bonds shown in green circle (A) Geraniol: with 5 rotatable bonds (B) Limonene: with 1 rotatable bond.

Materials and Methods

Preparation of the Receptor and the Ligands
The structure file of HMG-CoA reductase complexed with atorvastatin (an inhibitor) was downloaded from Protein Data Bank (PDB id: 1HWK). Structure was resolved using x-ray crystallography experiment at 2.22 Å resolution with R-value 0.212 from Homo sapiens (Istvan and Deisenhofer, 2001). To study the interaction of HMG-CoA reductase with geraniol and limonene, water molecules and non-protein residues were deleted from the complex. CHARMm forcefield was applied to study the molecular dynamics. CHARMm uses a flexible and comprehensive empirical energy function that is a summation of many individual energy terms. The energy function is based on separable internal coordinate terms and pairwise nonbond interaction terms (Brooks et al., 1983). The total energy is expressed by equation 1.

  E = Eb + Eθ + Eφ + Eω + EvdW + Eel + Ehb + ECr + ECφ
(1)

where, E is the total energy; Eb(bond potential), Eθ(bond angle potential), Eφ(Dihedral angle potential), Eω(improper torsions) are internal energy terms, EvdW(Van der Waals interactions), Eel(Electrostatic potential), Ehb(hydrogen bond energy) are nonbonded internal/external interactions energy terms, ECr(constraints) and ECφ(user defined energy function) are special energy terms.

Identification of Binding Cavity on Receptor Surface
After energy minimization, the binding pockets of the receptor were determined by using “eraser” algorithm using Accelrys Discovery Studio. This algorithm is first used to remove all grid points outside the receptor. The boundary between the “inside” and “outside” region is determined by the “site opening” parameter. For the remaining grid points (i.e., those “inside” the site), a flood-filling algorithm is employed to find contiguous regions consisting of unoccupied, connected grid points. Each such region is identified as a possible site. A user-specified size cutoff used to remove sites smaller than the specified volume for further consideration (Venkatachalam et al., 2003).

Interaction Protocol and Scoring Functions for Docking
The interaction of the receptor and the ligand was performed using “LigandFit” protocol on Accelrys Discovery Studio. In the first phase of LigandFit docking procedure, binding sites were indentified on the receptor surface. Site partitioning approach was followed to sample different parts of the larger binding site for docking. In the second phase, docking between receptor and ligand was performed in the specified site.

Docking ligands to the specified sites has different approaches like conformational search to generate candidate ligand conformations for docking, ligand/site shape matching to select ligand
conformations that are similar to the shape of site or site partitions. Candidate ligand poses in the binding site are evaluated and prioritized according to the DockScore function on the basis of forcefield approximation (equation 2), Piecewise Linear Potential function (PLP) (equation 3), LigScore1, LigScore2, Potential of Mean Force (PMF) and Jain scores.

  DockScore(forcefield) = - (ligand/receptor interaction energy + ligand internal energy)
(2)

  DockScore(PLP) = - (PLP potential)
(3)

As shown in Eq. 2, this version of DockScore contain two energy terms, these are internal energy of the ligand and the interaction energy of the ligand with the receptor. The interaction energy is taken as the sum of the van der Waals energy and electrostatic energy. To reduce the time needed for the computation of the interaction energy, a grid-based estimation of the ligand/receptor interaction energy is employed. PLP is a fast, simple, docking function that has been shown to correlate well with protein-ligand binding affinities. PLP scores are measured in arbitrary units, with negative PLP scores reported in order to make them suitable for subsequent use in consensus score calculations.
Higher PLP scores indicate stronger receptor-ligand binding (larger pKi values). LigScore1 is a scoring function for predicting receptor-ligand binding affinities. vdW, C+pol and TotPol^2 descriptors are used to calculate LigScore1 (equation 4, 5), which is computed in units of pKi (-log Ki). When scoring
ligands, the individual contributions of these descriptors may also be provided along with the overall LigScore1 value. Two slightly different equations are used in the calculation of LigScore1 depending on the forcefield (Dreiding or CFF) employed for the calculation of the vdW descriptor and the corresponding charge model (Gasteiger or CFF) used to assign atoms as polar or nonpolar.

LigScore2 is another fast and simple scoring function for predicting receptor-ligand binding affinities. vdW, C+pol, and BuryPol^2 descriptors are used to calculate LigScore2 (equation 6,7), which is computed in units of pKi (-log Ki). When scoring ligands, the individual contributions of these descriptors may also be provided along with the overall LigScore2 value. Two slightly different equations are used in the calculation of LigScore2 depending on the forcefield (Dreiding or CFF) employed for the calculation of the vdW descriptor.

  LigScore1_CFF = 0.4896 - 0.04551*vdW + 0.1439*C+pol - 0.001010*TotPol^2
(4)

  LigScore1_Dreiding = -0.3498 - 0.04673*vdW + 0.1653*C+pol -0.001132*TotPol^2
(5)

  LigScore2_CFF = 1.900 - 0.0730*vdW + 0.06246*C+pol - 0.00007324*BuryPol^2
(6)

  LigScore2_Dreiding = 1.539 - 0.07622*vdW + 0.6501*C+pol - 0.00007821*BuryPol^2
(7)

where the coefficients were obtained through regression analysis of the binding affinities of a series of protein-ligand complexes (Krammer et al., 2005).

The PMF scoring function (Muegge et al., 2005) is based on statistical analysis of the 3D structures of protein-ligand complexes. They were found to correlate well with protein-ligand binding free energies while being fast and simple to calculate. The scores are calculated by summing pairwise interaction terms over all interatomic pairs of the receptor-ligand complex.

The Jain score is a sum of five interaction terms (Jain, 1996). These are Lipophilic interactions, Polar attractive interactions, Polar repulsive interactions, Solvation of the protein and ligand and an entropy term for the ligand. Only proximate protein-ligand atoms are considered for the pairwise interaction terms. The lipophilic and polar interaction terms are each represented by a weighted sum of a Gaussian and a sigmoidal function. This functional form is short-ranged with a pronounced maximum that occurs at close surface contacts. It also incurs a significant penalty for short contacts between protein and ligand atoms.

Parameters for Docking Study
For docking study, the Energy Grid Force Field parameter was set to Dreiding, for computing ligand-protein interaction energy. The Energy Grid parameters control the grid bases docking used in the initial evaluation of the poses. In the Dreiding force field the Gasteiger charging method is employed to calculate the partial charges of ligands and proteins. The Energy Grid Extension from site was set to 5.0 Å. The Conformation search Number of Monte Carlo Trial was set to “0” to perform a rigid
docking. Maximum poses for ligand in the receptor cavity was set to 10. Ligand poses in the receptor cavity were evaluated using LigScore1, LigScore2, PLP1, PLP2, PMF, Jain, Dock Score empirical scoring functions.

Result and Discussion
Molecular properties of genaniol and limonene were analysed, to identify, if they are satisfying Lipinski rule of 5. According to Lipinski rule of 5, for any druggable compound, molecular weight should be less than 500; number of H-donors less than 5; number of H-acceptor less than 10; and octanol-water partition coefficient (ALogP) value should be less than 5. Calculated molecular properties values of geraniol and limonene are shown in Table 1. Rotatable bonds of genaniol and limonene are shown in Figure 3. Geraniol contains total 5 rotatable bonds, while limonene has only 1 rotatable bond. Ligand conformations were generated using search small molecule confirmation tools available in Accelrys discovery studio. Systematic search method was used with energy threshold 20 kcal/mol to generate total 56 conformation poses of geraniol (Table 2). Energy plot
of all 56 confirmation poses of geraniol is shown in Figure 4. In limonene, only 1 rotatable bond was present, thus only 3 conformation poses were generated using systematic search conformation generation method at energy threshold 20 kcal/mol (Table 3). Energy plot of all 3 confirmation poses of limonene is shown in Figure 5. HMG-CoA reductase chain A was analysed for all the possible binding sites (Figure 6). These active sites were selected from the receptor according to their volume of the binding cavity. Docking was performed by selecting one site at a time. Binding site 2 with 326 interacting points and 40.750 Å ^3 volume showed interaction with both geraniol and limonene (Figure 7). Best poses for each geraniol and limonene with HMG-CoA reductase were analysed for different energy parameters.

Table 1: Molecular properties of Geraniol and Limonene.

Table 2: Confirmation poses of Geraniol generated using systematic search with energy threshold 20 kcal/mol.

Table 3: Confirmation poses of Limonene generated using systematic search with energy threshold 20 kcal/mol.

Figure4: Energy plot of all 56 confirmation poses of Geraniol.

Figure5: Energy plot of all 3confirmation poses of Limonene.

Figure6: 3-hydroxy-3-methylglutaryl-CoA (HMG-CoA) reductase (chain A). Molecular surface is colored based on calculated interpolated charges. Protein back bone is displayed as solid ribbon and colored by secondary structure type. Binding site 2 is shown with green dots.

Figure7: Detail view of binding site in 3-hydroxy-3- methylglutaryl-CoA (HMG-CoA) reductase (chain A). All amino acids, surrounding the binding sites are labeled. Surface is colored based on calculated interpolated charges.

HMG-CoA reductase -geraniol Interaction
Docking was performed with all 56 conformation poses of geraniol and top 8 poses were analysed in the binding cavity of HMG-CoA reductase (Table 4). The best pose of geraniol (with dock score = 9.448) interacting with threonine 809, aspartic acid 767 and glycine 765 of HMG-CoA reductase. The hydrogen atom at position 29 of geraniol interacts with hydrogen atom at position 22 of threonine present at position 809 of HMG-CoA reductase. Same hydrogen atom at position 29 of geraniol interacts with oxygen of the C=O of glycine present at position 765 of the receptor molecule. Hydrogen at position 20 of geraniol interacts with hydrogen beta 1 of aspartic acid at position 767.
Total ligand internal energy for the best post is calculated to 7.642. (Figure 8).

Table 4: Conformation poses of geraniol and limonene with different scoring functions. Poses are arranged with the descending dock score value.

Figure8: Geraniol interaction with HMG-CoA reductase. All interacting amino acids with the ligand are labeled.

HMG-CoA reductase -limonene Interaction
For limonene and HMG-CoA reductase interaction, we generated three conformations using small molecule conformation generation method. All three poses of limonene were analyzed in the binding cavity of HMG-CoA reductase. The best pose of limonene has dock score as high as 10.593, much more than the best dock score in case of geraniol (Table 4). Best pose of limonene has three atoms interacting withHMG-CoA reductase. These are hydrogen at position 20, 17 and carbon at position 3. Hydrogen at position 20 and carbon at position 3 form bond with hydrogen atom of glycine amino terminus at position 808. The second interaction is in between hydrogen at position 17 of limonene and H-β1 of aspartic acid present at position 767 of HMG-CoA reductase (Figure 9).

Figure9: Limonene interaction with HMG-CoA reductase. All interacting amino acids with the ligand are labeled.

Conclusion
Molecular docking studies provide lead to determine the potential of ligand interaction in the binding cavities of receptor molecules. Considering the high dock score and low ligand internal energy, it can be concluded that limonene has greater binding affinity with HMG-CoA reductase and thus having better antitumor activity in comparison to geraniol. Aspartic acid at position 767 of HMG-CoA reductase is interacting with both geraniol and limonene. This amino acid acts as a major anchor
point for the ligands to interact with the receptor molecule for their anti-tumor activities.

Acknowledgements
SKG acknowledge Director, IITR, Lucknow for his generous support and CSIR NWP-17 project for providing necessary infrastructure for the work.

References
  1. Alonso WR, Rajaonarivony JI, Gershenzon J, Croteau R (1992) Purification of 4S-limonene synthase, a monoterpene cyclase from the glandular trichomes of peppermint (Mentha x piperita) and spearmint (Mentha spicata). J Biol Chem 267: 7582-7. »  CrossRef  »  PubMed  »  Google Scholar

  2. Asamoto M, Ota T, Toriyama-Baba H, Hokaiwado N, Naito A, et al. (2002) Mammary carcinomas induced in human c- Ha-ras proto-oncogene transgenic rats are estrogen-independent, but responsive to d-limonene treatment. Jpn J Cancer Res 93: 32-5. »  CrossRef  »  PubMed  »  Google Scholar

  3. Brooks BR, Bruccoleri RE, Olafson BD, States DJ, Swaminathan S, et al. (1983) CHARMM: A program for macromolecular energy minimization and dynamics calculations. J Comput Chem 4: 187-217. »  CrossRef  »  Google Scholar

  4. Carnesecchi S, Bras-Gonçalves R, Bradaia A, Zeisel M, Gossé F, et al. (2004) Geraniol, a component of plant essential oils, modulates DNA synthesis and potentiates 5-fluorouracil efficacy on human colon tumor xenografts. Cancer Lett 215: 53-9. »  CrossRef  »  PubMed  »  Google Scholar

  5. Chander SK, Lansdown AG, Luqmani YA, Gomm JJ, Coope RC, et al. (1994) Effectiveness of combined limonene and 4-hydroxyandrostenedione in the treatment of NMU-induced rat mammary tumours. Br J Cancer 69: 879-82. »  PubMed  »  Google Scholar

  6. Clarke S (1992) Protein isoprenylation and methylation at carboxyl-terminal cysteine residues. Annu Rev Biochem 61: 355-86. »  CrossRef  »  PubMed  »  Google Scholar

  7. Crowell PL (1999) Prevention and therapy of cancer by dietary monoterpenes. J Nutr 129: 775S-778S. »  CrossRef  »  PubMed  »  Google Scholar

  8. Elegbede JA, Maltzman TH, Verma AK, Tanner MA, Elson CE, et al. (1986) Mouse skin tumor promoting activity of orange peel oil and d-limonene: a re-evaluation. Carcinogenesis 7: 2047-9. »  CrossRef  »  PubMed  »  Google Scholar

  9. Elson CE, Maltzman TH, Boston JL, Tanner MA, Gould MN (1988) Anti-carcinogenic activity of d-limonene during the initiation and promotion/progression stages of DMBA-induced rat mammary carcinogenesis. Carcinogenesis 9: 331- 2. »  CrossRef  »  PubMed  »  Google Scholar

  10. Genser B, Grammer TB, Stojakovic T, Siekmeier R, März W (2008) Effect of HMG CoA reductase inhibitors on lowdensity lipoprotein cholesterol and C-reactive protein: systematic review and meta-analysis. Int J Clin Pharmacol Ther 46: 497-510.  »  PubMed  »  Google Scholar

  11. Gould MN, Moore CJ, Zhang R, Wang B, Kennan WS, et al. (1994) Limonene chemoprevention of mammary carcinoma induction following direct in situ transfer of v-Ha-ras. Cancer Res 54: 3540-3543. »  CrossRef  »  PubMed  »  Google Scholar

  12. Iijima Y, Davidovich-Rikanati R, Fridman E, Gang DR, Bar E, et al. (2004) The Biochemical and molecular basis for the divergent patterns in the biosynthesis of terpenes and phenylpropenes in the peltate glands of three cultivars of Basil. Plant Physiol 136: 3724-3736. »  CrossRef  »  PubMed  »  Google Scholar

  13. Istvan ES, Deisenhofer J (2001) Structural mechanism for statin inhibition of HMG-CoA reductase. Science 292: 1160- 4. »  CrossRef  »  PubMed  »  Google Scholar

  14. Jain AN (1996) Scoring noncovalent protein-ligand interactions: a continuous differentiable function tuned to compute binding affinities. J Comput Aided Mol Des 10: 427-40. »  CrossRef  »  PubMed  »  Google Scholar

  15. Kawamori T, Tanaka T, Hirose Y, Ohnishi M, Mori H (1996) Inhibitory effects of d-limonene on the development of colonic aberrant crypt foci induced by azoxymethane in F344 rats. Carcinogenesis. 17: 369-372. »  CrossRef  »  PubMed  »  Google Scholar

  16. Kjonaas R, Croteau R (1983) Demonstration that limonene is the first cyclic intermediate in the biosynthesis of oxygenated p-menthane monoterpenes in Mentha Piperita and other Mentha species. Arch Biochem Biophys 220: 79-89. »  CrossRef  »  PubMed  »  Google Scholar

  17. Kodama R, Yano T, Furukawa K, Noda K, Ide H (1976) Studies on the metabolism of d-limonene (p-mentha-1,8-diene). IV. Isolation and characterization of new metabolites and species differences in metabolism. Xenobiotica 6: 377- 89. »  CrossRef  »  PubMed  »  Google Scholar

  18. Krammer A, Kirchhoff PD, Jiang X, Venkatachalam CM, Waldman M (2005) LigScore: a novel scoring function for predicting binding affinities. J Mol Graph Model 23: 395- 407. »  CrossRef  »  PubMed  »  Google Scholar

  19. Loza-Tavera H (1999) Monoterpenes in essential oils. Biosynthesis and properties. Adv Exp Med Biol 464: 49-62. »  PubMed  »  Google Scholar

  20. Lu XG, Zhan LB, Feng BA, Qu MY, Yu LH, et al. (2004) Inhibition of growth and metastasis of human gastric cancer implanted in nude mice by d-limonene. World J Gastroenterol 10: 2140-4. »  CrossRef  »  PubMed  »  Google Scholar

  21. McGarvey DJ, Croteau R (1995) Terpenoid metabolism. Plant Cell 7: 1015-1026. »  CrossRef  »  PubMed  »  Google Scholar

  22. Muegge I, Martin YC (1999) A general and fast scoring function for protein-ligand interactions: a simplified potential approach. J Med Chem 42: 791-804. »  CrossRef  »  PubMed  »  Google Scholar

  23. National Toxicology Program (1990) NTP Toxicology and Carcinogenesis Studies of d-Limonene (CAS No. 5989-27- 5) in F344/N Rats and B6C3F1 Mice (Gavage Studies). Natl Toxicol Program Tech Rep Ser 347: 1-165.  »  PubMed  

  24. Peffley DM, Gayen AK (2003) Plant-derived monoterpenes suppress hamster kidney cell 3-hydroxy-3-methylglutaryl coenzyme a reductase synthesis at the post-transcriptional level. J Nutr 133: 38-44. »  CrossRef  »  PubMed  »  Google Scholar

  25. Polo MP, de Bravo MG (2006) Effect of geraniol on fattyacid and mevalonate metabolism in the human hepatoma cell line Hep G2. Biochem Cell Biol 84: 102-11. »  CrossRef  »  PubMed  »  Google Scholar

  26. Raphael TJ, Kuttan G (2003) Effect of naturally occurring monoterpenes carvone, limonene and perillic acid in the inhibition of experimental lung metastasis induced by B16F- 10 melanoma cells. J Exp Clin Cancer Res 22: 419-24.  »  PubMed  »  Google Scholar

  27. Salim EI, Wanibuchi H, Morimura K, Wei M, Mitsuhashi M, et al. (2003) Carcinogenicity of dimethylarsinic acid in p53 heterozygous knockout and wild-type C57BL/6J mice. Carcinogenesis 24: 335-42. »  CrossRef  »  PubMed  »  Google Scholar

  28. Singletery K (2000) Diet, natural products and cancer chemoprevention. J Nutr 130: 465S-466S. »  CrossRef  »  PubMed  »  Google Scholar

  29. Swanson KM, Hohl RJ (2006) Anti-cancer therapy: targeting the mevalonate pathway. Curr Cancer Drug Targets 6: 15-37. »  CrossRef  »  PubMed  »  Google Scholar

  30. Venkatachalam CM, Jiang X, Oldfield T, Waldman M (2003) LigandFit: a novel method for the shape-directed rapid docking of ligands to protein active sites. J Mol Graph Model 21: 289-307. »  CrossRef  »  PubMed  »  Google Scholar
 
This Article
DOWNLOAD
» XML (76 KB)
» PDF (2, 450 KB)
» Citation
CONTRIBUTE

SHARE

EXPLORE
Related Article at