Not so crystal clear: the structure of the human telomere G-quadruplex in solution differs from that present in a crystal

The structure of human telomere DNA is of intense interest because of its role in the biology of both cancer and aging. The sequence [5′-AGGG(TTAGGG)(3)] has been used as a model for telomere DNA in both NMR and X-ray crystallographic studies, the results of which show dramatically different structures. In Na(+) solution, NMR revealed an antiparallel G-quadruplex structure that featured both diagonal and lateral TTA loops. Crystallographic studies in the presence of K(+) revealed a flattened, propeller-shaped structure featuring a parallel-stranded G-quadruplex with symmetrical external TTA loops. We report the results of biophysical experiments in solution and computational studies that are inconsistent with the reported crystal structure, indicating that a different structure exists in K(+) solutions. Sedimentation coefficients were determined experimentally in both Na(+) and K(+) solutions and were compared with values calculated using bead models for the reported NMR and crystal structures. Although the solution NMR structure accurately predicted the observed S-value in Na(+) solution, the crystal structure predicted an S-value that differed dramatically from that experimentally observed in K(+) solution. The environments of loop adenines were probed by quantitative fluorescence studies using strategic and systematic single-substitutions of 2-aminopurine for adenine bases. Both fluorescence intensity and quenching experiments in K(+) yielded results at odds with quantitative predictions from the reported crystal structure. Circular dichroism and fluorescence quenching studies in the presence of the crowding agent polyethylene glycol showed dramatic changes in the quadruplex structure in K(+) solutions, but not in Na(+) solutions, suggesting that the crystal environment may have selected for a particular conformational form. Molecular dynamics simulations were performed to yield model structures for the K(+) quadruplex form that are consistent with our biophysical results and with previously reported chemical modification studies. These models suggest that the biologically relevant structure of the human telomere quadruplex in K(+) solution is not the one determined in the published crystalline state

Electrostatic free energy calculations using the generalized solvent boundary potential method

Copyright 2002 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics. The following article appeared in The Journal of Chemical Physics and may be found at http://dx.doi.org/10.1063/1.1507108.Free energyperturbation (FEP) calculations using all-atom molecular dynamics simulations with a large number of explicit solvent molecules are a powerful approach to study ligand–macromolecule association processes at the atomic level. One strategy to carry out FEP calculations efficiently and reduce computational time is to consider the explicit dynamics of only a small number of atoms in a localized region around the ligand. Such an approximation is motivated by the observation that the factors governing binding specificity are dominated by interactions in the vicinity of the ligand. However, a straightforward truncation of the system may yield inaccurate results as the influence exerted by the remote regions of the macromolecule and the surrounding solvent through long-range electrostatic effects may be significant. To obtain meaningful results, it is important to incorporate the influence of the remote regions of the ligand–macromolecule complex implicitly using some effective potential. The generalized solvent boundary potential (GSBP) that was developed recently [W. Im, S. Bernèche, and B. Roux, J. Chem. Phys. 114, 2924 (2001)] is an efficient computational method to represent the long-range electrostaticinteractions arising from remote (outer) regions in simulations of a localized (inner) region with a small number of explicit atoms. In the present work, FEP calculations combined with GSBP are used to illustrate the importance of these long-range electrostatic factors in estimation of the charging free energy of an aspartate ligand bound to the aspartyl-tRNA synthetase. Calculations with explicit spherical simulation inner regions of different radii are used to test the accuracy of the GSBP method and also illustrate the importance of explicit protein and solvent dynamics in the free energy estimation. The influence of the represented outer region is tested using separate simulations in which the reaction field and/or the protein static field are excluded. Both components are shown to be essential to obtain quantitatively meaningful results. The ability of implicitly treating the influence of protein fluctuations in the outer region using a protein dielectric constant is examined. It is shown that accurate charging free energy calculations can be performed for this system with a spherical region of 15 to 20 Å radius, which roughly corresponds to 1500–3500 moving atoms. The results indicate that GSBP in combination with FEP calculations is a precise and efficient approach to include long-range electrostatic effects in the study of ligand binding to large macromolecules

Ranking ligand affinity for the DNA minor groove by experiment and simulation

The structural and thermodynamic basis for the strength and selectivity of the interactions of minor-groove binders (MGBs) with DNA is not fully understood. In 2003 we reported the first example of a thiazole containing MGB that bound in a phase shifted pattern that spanned 6 base-pairs rather than the usual 4 (for tricyclic distamycin-like compounds). Since then, using DNA footprinting, nuclear magnetic resonance spectroscopy, isothermal titration calorimetry and molecular dynamics, we have established that the flanking bases around the central 4 being read by the ligand have subtle effects on recognition. We have investigated the effect of these flanking sequences on binding and the reasons for the differences and established a computational method to rank ligand affinity against varying DNA sequences

Microcalorimetry and spectroscopic studies on the binding of dye janus green blue to deoxyribonucleic acid

The interaction of the phenazinium dye janus green blue (JGB) with deoxyribonucleic acid was investigated using isothermal titration calorimetry and thermal melting experiments. The calorimetric data were supplemented by spectroscopic studies. Calorimetry results suggested the binding affinity of the dye to DNA to be of the order of 105 M-1. The binding was predominantly entropy driven with a small negative favorable enthalpy contribution to the standard molar Gibbs energy change.The binding became weaker as the temperature and salt concentration was raised. The temperature dependence of the standard molar enthalpy changes yielded negative values of standard molar heat capacity change for the complexation revealing substantial hydrophobic contribution in the DNA binding. An enthalpy–entropy compensation behavior was also observed in the system. The salt dependence of the binding yielded the release of 0.69 number of cations on binding of each dye molecule. The non-polyelectrolytic contribution was found to be the predominant force in the binding interaction. Thermal melting studies revealed that the DNA helix was stabilized against denaturation by the dye. The binding was also characterized by absorbance, resonance light scattering and circular dichroism spectral measurements. The binding constants from the spectral results were close to those obtained from the calorimetric data. The energetic aspects of the interaction of the dye JGB to double stranded DNA are supported by strong binding revealed from the spectral data

Characterization of a K+-induced conformational switch in a human telomeric DNA oligonucleotide using 2-aminopurine fluorescence

Human telomeric DNA consists of tandem repeats of the DNA sequence d(GGGTTA). Oligodeoxynucletotide telomere models such as d[A(GGGTTA)(3)GGG] (Tel22) fold in a cation-dependent manner into quadruplex structures consisting of stacked G-quartets linked by d(TTA) loops. NMR has shown that in Na(+) solutions Tel22 forms a ‘basket’ topology of four antiparallel strands; in contrast, Tel22 in K(+) solutions consists of a mixture of unknown topologies. Our previous studies on the mechanism of folding of Tel22 and similar telomere analogs utilized changes in UV absorption between 270 and 325 nm that report primarily on G-quartet formation and stacking showed that quadruplex formation occurs within milliseconds upon mixing with an appropriate cation. In the current study, we assessed the dynamics and equilibria of folding of specific loops by using Tel22 derivatives in which the dA residues were serially substituted with the fluorescent reporter base, 2-aminopurine (2-AP). Tel22 folding induced by Na(+) or K(+) assessed by changes in 2-AP fluorescence consists of at least three kinetic steps with time constants spanning a range of ms to several hundred seconds. Na(+)-dependent equilibrium titrations of Tel22 folding could be approximated as a cooperative two-state process. In contrast, K(+)-dependent folding curves were biphasic, revealing that different conformational ensembles are present in 1 mM and 30 mM K(+). This conclusion was confirmed by (1)H NMR. Molecular dynamics simulations revealed a K(+) binding pocket in Tel22 located near dA1 that is specific for the so-called hybrid-1 conformation in which strand 1 is in a parallel arrangement. The possible presence of this topologically specific binding site suggests that K(+) may play an allosteric role in regulating telomere conformation and function by modulating quadruplex tertiary structure

Cellular uptake, cytotoxicity and DNA-binding studies of the novel imidazoacridinone antineoplastic agent C1311

C1311 is a novel therapeutic agent with potent activity against experimental colorectal cancer that has been selected for entry into clinical trial. The compound has previously been shown to have DNA-binding properties and to inhibit the catalytic activity of topoisomerase II. In this study, cellular uptake and mechanisms by which C1311 interacts with DNA and exerts cytotoxic effects in intact colon carcinoma cells were investigated. The HT29 colon cancer cell line was chosen to follow cellular distribution of C1311 over a time course of 24 h at drug concentrations that just inhibited cell proliferation by 50% or 100%. Nuclear uptake of C1311 and co-localization with lysosomal or mitochondrial dyes was examined by fluorescence microscopy and effects on these cellular compartments were determined by measurement of acid phosphatase levels, rhodamine 123 release or DNA-binding behaviour. The strength and mode of DNA binding was established by thermal melting stabilization, direct titration and viscometric studies of host duplex length. The onset of apoptosis was followed using a TUNEL assay and DNA-fragmentation to determine a causal relationship of cell death. Growth inhibition of HT29 cells by C1311 was concomitant with rapid drug accumulation in nuclei and in this context we showed that the compound binds to duplex DNA by intercalation, with likely A/T sequence-preferential binding. Drug uptake was also seen in lysosomes, leading to lysosomal rupture and a marked increase of acid phosphatase activity 8 h after exposure to C1311 concentrations that effect total growth inhibition. Moreover, at these concentrations lysosomal swelling and breakdown preceded apoptosis, which was not evident up to 24 h after exposure to drug. Thus, the lysosomotropic effect of C1311 appears to be a novel feature of this anticancer agent. As it is unlikely that C1311-induced DNA damage alone would be sufficient for cytotoxic activity, lysosomal rupture may be a critical component for therapeutic efficacy. © 1999 Cancer Research Campaig

Unprocessed Viral DNA Could Be the Primary Target of the HIV-1 Integrase Inhibitor Raltegravir

Integration of HIV DNA into host chromosome requires a 3′-processing (3′-P) and a strand transfer (ST) reactions catalyzed by virus integrase (IN). Raltegravir (RAL), commonly used in AIDS therapy, belongs to the family of IN ST inhibitors (INSTIs) acting on IN-viral DNA complexes (intasomes). However, studies show that RAL fails to bind IN alone, but nothing has been reported on the behaviour of RAL toward free viral DNA. Here, we assessed whether free viral DNA could be a primary target for RAL, assuming that the DNA molecule is a receptor for a huge number of pharmacological agents. Optical spectroscopy, molecular dynamics and free energy calculations, showed that RAL is a tight binder of both processed and unprocessed LTR (long terminal repeat) ends. Complex formation involved mainly van der Waals forces and was enthalpy driven. Dissociation constants (Kds) revealed that RAL affinity for unbound LTRs was stronger than for bound LTRs. Moreover, Kd value for binding of RAL to LTRs and IC50 value (half concentration for inhibition) were in same range, suggesting that RAL binding to DNA and ST inhibition are correlated events. Accommodation of RAL into terminal base-pairs of unprocessed LTR is facilitated by an extensive end fraying that lowers the RAL binding energy barrier. The RAL binding entails a weak damping of fraying and correlatively of 3′-P inhibition. Noteworthy, present calculated RAL structures bound to free viral DNA resemble those found in RAL-intasome crystals, especially concerning the contacts between the fluorobenzyl group and the conserved 5′C4pA33′ step. We propose that RAL inhibits IN, in binding first unprocessed DNA. Similarly to anticancer drug poisons acting on topoisomerases, its interaction with DNA does not alter the cut, but blocks the subsequent joining reaction. We also speculate that INSTIs having viral DNA rather IN as main target could induce less resistance

Interaction of Pyrrolobenzodiazepine (PBD) Ligands with Parallel Intermolecular G-Quadruplex Complex Using Spectroscopy and ESI-MS

Studies on ligand interaction with quadruplex DNA, and their role in stabilizing the complex at concentration prevailing under physiological condition, has attained high interest. Electrospray ionization mass spectrometry (ESI-MS) and spectroscopic studies in solution were used to evaluate the interaction of PBD and TMPyP4 ligands, stoichiometry and selectivity to G-quadruplex DNA. Two synthetic ligands from PBD family, namely pyrene-linked pyrrolo[2,1-c][1,4]benzodiazepine hybrid (PBD1), mixed imine-amide pyrrolobenzodiazepine dimer (PBD2) and 5,10,15,20-tetrakis(N-methyl-4-pyridyl)porphyrin (TMPyP4) were studied. G-rich single-stranded oligonucleotide d(5′GGGGTTGGGG3′) designated as d(T2G8), from the telomeric region of Tetrahymena Glaucoma, was considered for the interaction with ligands. ESI-MS and spectroscopic methods viz., circular dichroism (CD), UV-Visible, and fluorescence were employed to investigate the G-quadruplex structures formed by d(T2G8) sequence and its interaction with PBD and TMPyP4 ligands. From ESI-MS spectra, it is evident that the majority of quadruplexes exist as d(T2G8)2 and d(T2G8)4 forms possessing two to ten cations in the centre, thereby stabilizing the complex. CD band of PBD1 and PBD2 showed hypo and hyperchromicity, on interaction with quadruplex DNA, indicating unfolding and stabilization of quadruplex DNA complex, respectively. UV-Visible and fluorescence experiments suggest that PBD1 bind externally where as PBD2 intercalate moderately and bind externally to G-quadruplex DNA. Further, melting experiments using SYBR Green indicate that PBD1 unfolds and PBD2 stabilizes the G-quadruplex complex. ITC experiments using d(T2G8) quadruplex with PBD ligands reveal that PBD1 and PBD2 prefer external/loop binding and external/intercalative binding to quadruplex DNA, respectively. From experimental results it is clear that the interaction of PBD2 and TMPyP4 impart higher stability to the quadruplex complex