Abstract
Telomeric ends of chromosomes, which comprise noncoding repeat sequences of guanine-rich DNA, which are the fundamental in protecting the cell from recombination and degradation. Telomeric DNA sequences can form four stranded quadruplex structures, which are involved in the structure of telomere ends. The formation and stabilization of telomeric quadruplexes has been shown to inhibit the activity of telomerase, thus establishing telomeric DNA quadrulex as an attractive target for cancer therapeutic intervention. Molecular dynamic simulation offers the prospects of detailed description of the dynamical structure with ion and water at molecular level. In this work we have taken a oligomeric part of human telomeric DNA, d(TAGGGT) to form different monomeric quadruplex structures d(TAGGGT)4. Here we report the relative stabilities of these structures under K+ ion conditions and binding interaction between the strands, as determined by molecular dynamic simulations followed by energy calculation. We have taken locked nucleic acid (LNA) in this study. The free energy molecular mechanics Poission Boltzman surface area calculations are performed for the determination of most stable complex structure between all modified structures. We calculated binding free energy for the combination of different strands as the ligand and receptor for all structures. The energetic study shows that, a mixed hybrid type quadruplex conformation in which two parallel strands are bind with other two antiparallel strands, are more stable than other conformations. The possible mechanism for the inhibition of the cancerous growth has been discussed. Such studies may be helpful for the rational drug designing.
Similar content being viewed by others
References
Feldser DM, Hackett JA, Greider CW (2003) Telomere dysfunction and the initiation of genome instability. Nat Rev Cancer 3:623–627
Cech TR (2000) Life at the end of the chromosome: telomeres and telomerase. Angew Chem Int Ed 39:34–43
Wright WE, Tesmer VM, Huffman KE, Levene SD, Shay JW (1997) Normal human chromosomes have long G-rich telomeric overhangs at one end. Genes Dev 11:2801–2809
Organisian L, Bryan TM (2007) Physiological relevance of telomeric G-quadruplex formation: a potential drug target. Bioessays 29:155–165
Kim NW, Piatyszek MA, Prowse KR, Harley CB, West MD, Ho PL, Coviello GM, Wright WE, Weinrich SL et al (1994) Specific association of human telomerase activity with immortal cells and cancer. Science 266:2011–2015
Hanahan D, Weinberg RA (2000) The hallmarks of cancer. Cell 100:57–70
Mergny JL, Helene C (1998) G-quadruplex DNA: a target for drug design. Nat Med 4:1366–1367
Sun DY, Hurley LH (2001) Methods in enzymology, drug-nucleic acid interactions. Academic Press, Inc 340:573–592
Hurley LH (2002) DNA and its associated processes as targets for cancer therapy. Nat Rev Cancer 2:188–200
Neidle S, Parkinson G (2002) Telomere maintenance as a target for anticancer drug discovery. Nat Rev Drug Discov 1:383–393
Bodnar AG, Ouellette M, Frolkis M, Holt SE, Chiu CP, Morin GB, Harley CB, Shay JW, Lichtsteiner S et al (1998) Extension of life-span by introduction of telomerase into normal human cells. Science 279:349–352
Harley CB (1991) Telomere loss: mitotic clock or genetic time bomb? Mutat Res 256:271–282
Sun H, Karow JK, Hickson ID, Maizels N (1998) The Bloom’s syndrome helicase unwinds G4 DNA. J Biol Chem 273:27587–27592
Sun H, Bennett RJ, Maizels N (1999) The saccharomyces cerevisiae Sgs1 helicase esciently unwinds G-G paired DNAs. Nucleic Acids Res 27:1978–1984
Hackett JA, Feldser DM, Greider CW (2001) Telomere dysfunction increases mutation rate and genomic instability. Cell 106:275–286
Salazar M, Thompson BD, Kerwin SM, Hurley LH (1996) Thermally induced DNA:RNA hybrid to G-quadruplex transitions: possible implications for telomere synthesis by telomerase. Biochemistry 35:16110–16115
Hurley LH (2001) Secondary DNA structures as molecular targets for cancer therapeutics. Biochem Soc Trans 29:692–696
Hurley LH, Wheelhouse RT, Sun D, Kerwin SM, Salazar M, Fedoro OY, Han FX, Han HY, Izbicka E et al (2000) G-quadruplexes as targets for drug design. Pharmacol Ther 85:141–158
Neidle S, Read MA (2000) G-quadruplexes as therapeutic targets. Biopolymers 56:195–208
Mergny J-L, Helene C (1998) G-quadruplex DNA: a target for drug design. Nat Genet 4:1366–1367
Bearss DJ, Hurley LH, Von Hoff DD (2000) Telomere maintenance mechanisms as a target for drug development. Oncogene 19:6632–6641
Gowan SM et al (2002) A G-Quadruplex-interactive potent small-molecule inhibitor of telomerase exhibiting in vitro and in vivo antitumor activity. Mol Pharmacol 61:1154–1162
Read MA et al (2001) Structure-based design of selective and potent G quadruplex-mediated telomerase inhibitors. Proc Natl Acad Sci USA 98(9):4844–4849
Parkinson GN, Lee MPH, Neidle S (2002) Crystal structure of parallel quadruplexes from human telomeric DNA. Nature 417:876–880
Ambrus A, Chen D, Dai JX, Bialis T, Jones RA, Yang DZ (2006) Human telomeric sequence forms a hybrid-type intra- molecular G-quadruplex structure with mixed parallel/antiparallel strands in potassium solution. Nucleic Acids Res 34:2723–2735
Xu Y, Noguchi Y, Sugiyama H (2006) The new models of the human telomere d[AGGG(TTAGGG)(3)] in K+ solution. Bioorg Med Chem 14:5584–5591
Luu KN, Phan AT, Kuryavyi V, Lacroix L, Patel DJ (2006) Structure of the human telomere in K+ solution: an intramolecular (3 + 1) G-quadruplex scafold. J Am Chem Soc 128:9963–9970
Wang Y, Patel DJ (1993) Solution structure of the human telomeric repeat d[AG3(T2AG3)3] G-tetraplex. Structure 1:263–282
He Y, Neumann RD, Panyutin IG (2004) Intramolecular quadruplex conformation of human telomeric DNA assessed with 125I-radioprobing. Nucleic Acid Res 32:5359–5367
Redon S, Bombard S, Elizondo-Riojas MA, Chottard JC (2003) Platinum cross-linking of adenines and guanines on the quadruplex structures of the AG3(T2AG3)3 and (T2AG3)4 human telomere sequences in Na+ and K+ solutions. Nucleic Acid Res 31:1605–1613
Qi J, Shafer RH (2005) Covalent ligation studies on the human telomere quadruplex. Nucleic Acid Res 33:3185–3192
Li J, Correia JJ, Wang L, Trent JO, Chaires JB (2005) Not so crystal clear: the structure of the human telomere G-quadruplex in solution differs from that present in a crystal. Nucleic Acid Res 33:4649–4659
Crooke ST (2004) Progress in antisence technology. Annu Rev Med 55:61–95
Bondensgaard K, Petersen M, Singh SK, Rajwanshi VK, Kumar R, Wengel J, Jacobsen JP (2000) Structural studies of LNA: RNA duplexes by NMR: conformations and implications for RNase H activity. Chemistry 6:2687–2695
Nielsen KE, Singh SK, Wengel J, Jacobsen JP (2000) Solution structure of an LNA hybridized to DNA: NMR study of the d(CT(L)GCT(L)T(L)CT(L)GC): d(GCAGAAGCAG) duplex containing four locked nucleotides. Bioconjug Chem 11:228–238
Petersen M, Bondensgaard K, Wengel J, Jacobsen JP (2002) Locked nucleic acid (LNA) recognition of RNA: NMR solution structures of LNA: RNA hybrids. J Am Chem Soc 124:5974–5982
Randazzo A, Esposito V, Ohlenschlager O, Ramachandran R, Mayol L (2004) NMR solution structure of a parallel LNA quadruplex. Nucleic Acids Res 32:3083–3092
Sacca B, Lacroix L, Mergny J-L (2005) The effect of chemical modifications on the thermal stability of different G-quadruplex-forming oligonucleotides. Nucleic Acids Res 33:1182–1192
Ivanova A, Rosch N (2007) The structure of LNA: DNA hybrids from molecular dynamics simulations: the effect of locked nucleotides. J Phys Chem A 111(38):9307–9319
Pande V, Nilsson L (2008) Insights into structure, dynamics and hydration of locked nucleic acid (LNA) strand-based duplexes from molecular dynamics simulation. Nucleic Acids Res 36(5):1508–1516
Cornell WD, Cieplak P, Baylyl CI, Gould IR, Merz JKM et al (1995) A second generation force field for the simulation of proteins, nucleic acids and organic molecules. J Am Chem Soc 117:5179–5197
Price DJ, Brooks CL 3rd (2004) A modified TIP3P water potential for simulation with Ewald summation. J Chem Phys 121:10096–10103
Duan Y, Wu C, Chowdhury S, Lee MC, Xiong G, Zhang W, Yang R, Cieplak P, Luo R, Lee T (2003) A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations. J Comput Chem 24:1999–2012
Lee MC, Duan Y (2004) Distinguish protein decoys by using a scoring function based on a new Amber force field, short molecular dynamics simulations, and the generalized born solvent model. Proteins 55:620–634
Yang L, Tan C, Hsieh M-J, Wang J, Duan Y, Cieplak P, Caldwell J, Kollman PA, Luo R (2006) New-generation Amber united-atom force field. J Phys Chem B 110:13166–13176
Case DA, Cheathem TE, Daren T, Gohlke H, Luo R, Merz KM, Onufriev A et al (2005) The AMBER biomolecular simulation programs. J Comput Chem 26:1668–1688
Ryckaert JP, Ciccotti G, Berendsen HJC (1977) Numerical integration of the Cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. J Comput Phys 23:327–341
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE (2004) UCSF Chimera—a visualization system for exploratory research and analysis. J Comput Chem 25:1605–1612
Gohlke H, Case DA (2004) Converging free energy estimates: MM-PB(GB)SA studies on the protein-protein complex Ras-Raf. J Comput Chem 25:238–250
Dubey KD, Ojha RP (2011) Binding free energy calculations with hybrid QM/MM methods for Abl kinase inhibitors. J Biol Phys 37:69–78
Hud NV, Schultze P, Sklenar V, Feigon J (1999) Binding sites and dynamics of ammonium ions in a telomere repeat DNA quadruplex. J Mol Biol 285:233–243
Spackova N, Berger I, Sponer J (1999) Nanosecond molecular dynamics simulation of parallel and antiparallel guanine quadruplex DNA molcules. J Am Chem Soc 121:5519–5534
Chowdhury S, Bansal M (2001) G-quadruplex structure can be stable with only some coordination sites being occupied by cations: a six-nanosecond molecular dynamics study. J Phys Chem B 105:7572–7578
Cavallari N, Calzolarii A, Garbesi A, Di Felice R (2006) Stability and migration of metal ions in G4-wires by molecular dynamics simulations. J Phys Chem B 110:26337–26348
Agrawal S, Ojha RP, Maiti S (2008) Energetics of the human Tel-22 quadruplex—telomestatin interaction: a molecular dynamics study. J Phys Chem B 112(22):6828–6836
Gu J, Leszczynski J, Bansal M (1999) A new insight into the structure and stability of Hoogsteen hydrogen-bonded G-tetrad: an ab initio SCF study. Chem Phys Lett 311:209
Arnott S, Chandrasekaran R, Marttila CM (1974) Structures for polyinosinic acid and polyguanylic acid. Biochem J 141:537–543
Phillips K, Dauter Z, Murchie AIH, Lilley DMJ, Luisi B (1997) The crystal structure of a parallel-stranded guanine Tetraplex at 0.95 Å resolution. J Mol Biol 273:171–182
Acknowledgments
The authors are thankful to DST, New Delhi for computational facility in the form of FIST scheme. We thankfully acknowledge the partial computational work at BRAF of C-DAC, Pune, India. KDD acknowledges CSIR for senior research fellowship.
Author information
Authors and Affiliations
Corresponding author
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
About this article
Cite this article
Chaubey, A.K., Dubey, K.D. & Ojha, R.P. Stability and free energy calculation of LNA modified quadruplex: a molecular dynamics study. J Comput Aided Mol Des 26, 289–299 (2012). https://doi.org/10.1007/s10822-012-9548-z
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10822-012-9548-z