G-quadruplexes can maintain their structure in the gas phase

Several very extended (0.5-1 micros) molecular dynamics (MD) simulations of parallel and antiparallel G-quadruplex DNA strongly suggest that in the presence of suitable cations the quadruplex not only remains stable in the gas phase, but also displays a structure that closely resembles that found in extended (25-ns long) trajectories in aqueous solution. In the absence of the crucial cations, the trajectories become unstable and in general the quadruplex structure is lost. To our knowledge, this is the first physiologically relevant structure of DNA for which very large MD simulations suggest that the structure in water and in the gas phase are indistinguishable.     

G-quadruplex DNA assemblies: loop length, cation identity, and multimer formation

G-rich DNA sequences are able to fold into structures called G-quadruplexes. To obtain general trends in the influence of loop length on the structure and stability of G-quadruplex structures, we studied oligodeoxynucleotides with random bases in the loops. Sequences studied are dGGGW(i)GGGW(j)GGGW(k)GGG, with W = thymine or adenine with equal probability, and i, j, and k comprised between 1 and 4. All were studied by circular dichroism, native gel electrophoresis, UV-monitored thermal denaturation, and electrospray mass spectrometry, in the presence of 150 mM potassium, sodium, or ammonium cations. Parallel conformations are favored by sequences with short loops, but we also found that sequences with short loops form very stable multimeric quadruplexes, even at low strand concentration. Mass spectrometry reveals the formation of dimers and trimers. When the loop length increases, preferred quadruplex conformations tend to be more intramolecular and antiparallel. The nature of the cation also has an influence on the adopted structures, with K(+) inducing more parallel multimers than NH4(+) and Na(+). Structural possibilities are discussed for the new quadruplex higher-order assemblies.     

Spectroscopic detection of DNA quadruplexes by vibrational circular dichroism

The four-stranded G-quadruplex motif is a conformation frequently adopted by guanine-rich nucleic acids that plays an important role in biology, medicine, and nanotechnology. Although vibrational spectroscopy has been widely used to investigate nucleic acid structure, association of particular spectral features with the quadruplex structure has to date been ambiguous. In this work, experimental IR absorption and vibrational circular dichroism (VCD) spectra of the model quadruplex systems d(G)(8) and deoxyguanosine-5'-monophosphate (5'-dGMP) were analyzed using molecular dynamics (MD) and quantum-chemical modeling. The experimental spectra were unambiguously assigned to the quadruplex DNA arrangement, and several IR and VCD bands related to this structural motif were determined. Involvement of MD in the modeling was essential for realistic simulation of the spectra. The VCD signal was found to be more sensitive to dynamical structural variations than the IR signal. The combination of the spectroscopic techniques with multiscale simulations provides extended information about nucleic acid conformations and their dynamics.     

Stability and cations coordination of DNA and RNA 14-mer G-quadruplexes: a multiscale computational approach

Molecular dynamics simulations have been used to study the differences between two DNA and RNA 14-mer quadruplexes of analogous sequences. Their structures present a completely different fold: DNA forms a bimolecular quadruplex containing antiparallel strands and diagonal loops; RNA forms an intrastrand parallel quadruplex containing a G-tetrad and an hexad, which dimerizes by hexad stacking. We used a multiscale computational approach combining classical Molecular dynamics simulations and density functional theory calculations to elucidate the difference in stability of the 2-folds and their ability in coordinating cations. The presence of 2'-OH groups in the RNA promotes the formation of a large number of intramolecular hydrogen bonds that account for the difference in fold and stability of the two 14-mers. We observe that the adenines in the RNA quadruplex play a key role in conserving the geometry of the hexad. We predict the cation coordination mode of the two quadruplexes, not yet observed experimentally, and we offer a rationale for the corresponding binding energies involved.     

Topology variation and loop structural homology in crystal and simulated structures of a bimolecular DNA quadruplex

The topology of DNA quadruplexes depends on the nature and number of the nucleotides linking G-quartet motifs. To assess the effects of a three-nucleotide TTT linker, the crystal structure of the DNA sequence d(G(4)T(3)G(4)) has been determined at 1.5 A resolution, together with that of the brominated analogue d(G(4)(Br)UTTG(4)) at 2.4 A resolution. Both sequences form bimolecular intermolecular G-quadruplexes with lateral loops. d(G(4)(Br)UTTG(4)) crystallized in the monoclinic space group P2(1) with three quadruplex molecules in the asymmetric unit, two associating together as a head-to-head stacked dimer, and the third as a single head-to-tail dimer. The head-to-head dimers have two lateral loops on the same G-quadruplex face and form an eight-G-quartet stack, with a linear array of seven K(+) ions between the quartets. d(G(4)T(3)G(4)) crystallized in the orthorhombic space group C222 and has a structure very similar to the head-to-tail dimer in the P2(1) unit cell. The sequence studied here is able to form several different folds; however, all four quadruplexes in the two structures have lateral loops, in contrast to the diagonal loops reported for the analogous quadruplex with T(4) loops. A total of seven independent T(3) loops were observed in the two structures. These can be classified into two discrete conformational classes, suggesting that these represent preferred loop conformations that are independent of crystal-packing forces.     

Insight into G-DNA structural polymorphism and folding from sequence and loop connectivity through free energy analysis

The lengths of G-tracts and their connecting loop sequences determine G-quadruplex folding and stability. Complete understanding of the sequence-structure relationships remains elusive. Here, single-loop G-quadruplexes were investigated using explicit solvent molecular dynamics (MD) simulations to characterize the effect of loop length, loop sequence, and G-tract length on the folding topologies and stability of G-quadruplexes. Eight loop types, including different variants of lateral, diagonal, and propeller loops, and six different loop sequences [d0 (i.e., no intervening residues in the loop), dT, dT(2), dT(3), dTTA, and dT(4)] were considered through MD simulation and free energy analysis. In most cases the free energetic estimates agree well with the experimental observations. The work also provides new insight into G-quadruplex folding and stability. This includes reporting the observed instability of the left propeller loop, which extends the rules for G-quadruplex folding. We also suggest a plausible explanation why human telomere sequences predominantly form hybrid-I and hybrid-II type structures in K(+) solution. Overall, our calculation results indicate that short loops generally are less stable than longer loops, and we hypothesize that the extreme stability of sequences with very short loops could possibly derive from the formation of parallel multimers. The results suggest that free energy differences, estimated from MD and free energy analysis with current force fields and simulation protocols, are able to complement experiment and to help dissect and explain loop sequence, loop length, and G-tract length and orientation influences on G-quadruplex structure.     

Molecular dynamics and thermodynamics of protein-RNA interactions: mutation of a conserved aromatic residue modifies stacking interactions and structural adaptation in the U1A-stem loop 2 RNA complex.

Molecular dynamics (MD) simulations and free energy component analysis have been performed to evaluate the molecular origins of the 5.5 kcal/mol destabilization of the complex formed between the N-terminal RNP domain of U1A and stem loop 2 of U1 snRNA upon mutation of a conserved aromatic residue, Phe56, to Ala. MD simulations, including counterions and water, have been carried out on the wild type and Phe56Ala peptide-stem loop 2 RNA complexes, the free wild type and Phe56Ala peptides, and the free stem loop 2 RNA. The MD structure of the Phe56Ala-stem loop 2 complex is similar to that of the wild type complex except the stacking interaction between Phe56 and A6 of stem loop 2 is absent and loop 3 of the peptide is more dynamic. However, the MD simulations predict large changes in the structure and dynamics of helix C and increased dynamic range of loop 3 for the free Phe56Ala peptide compared to the wild type peptide. Since helix C and loop 3 are highly variable regions of RNP domains, this indicates that a significant contribution to the reduced affinity of the Phe56Ala peptide for RNA results from cooperation between highly conserved and highly variable regions of the RNP domain of U1A. Surprisingly, these structural effects, which are manifested as cooperative free energy changes, occur in the free peptide, rather than in the complex, and are revealed only by study of both the initial and final states of the complexation process. Free energy component analysis correctly accounts for the destabilization of the Phe56Ala-stem loop 2 complex, and indicates that approximately 80% of the destabilization is due to the loss of the stacking interaction and approximately 20% is due to differences in U1A adaptation.     

Comparison between self-guided Langevin dynamics and molecular dynamics simulations for structure refinement of protein loop conformations

This article presents a comparative analysis of two replica‐exchange simulation methods for the structure refinement of protein loop conformations, starting from low‐resolution predictions. The methods are self‐guided Langevin dynamics (SGLD) and molecular dynamics (MD) with a Nosé–Hoover thermostat. We investigated a small dataset of 8‐ and 12‐residue loops, with the shorter loops placed initially from a coarse‐grained lattice model and the longer loops from an enumeration assembly method (the Loopy program). The CHARMM22 + CMAP force field with a generalized Born implicit solvent model (molecular‐surface parameterized GBSW2) was used to explore conformational space. We also assessed two empirical scoring methods to detect nativelike conformations from decoys: the all‐atom distance‐scaled ideal‐gas reference state (DFIRE‐AA) statistical potential and the Rosetta energy function. Among the eight‐residue loop targets, SGLD out performed MD in all cases, with a median of 0.48 Å reduction in global root‐mean‐square deviation (RMSD) of the loop backbone coordinates from the native structure. Among the more challenging 12‐residue loop targets, SGLD improved the prediction accuracy over MD by a median of 1.31 Å, representing a substantial improvement. The overall median RMSD for SGLD simulations of 12‐residue loops was 0.91 Å, yielding refinement of a median 2.70 Å from initial loop placement. Results from DFIRE‐AA and the Rosetta model applied to rescoring conformations failed to improve the overall detection calculated from the CHARMM force field. We illustrate the advantage of SGLD over the MD simulation model by presenting potential‐energy landscapes for several loop predictions. Our results demonstrate that SGLD significantly outperforms traditional MD in the generation and populating of nativelike loop conformations and that the CHARMM force field performs comparably to other empirical force fields in identifying these conformations from the resulting ensembles. Published 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

NMR Structure of a Triangulenium‐Based Long‐Lived Fluorescence Probe Bound to a G‐Quadruplex

An NMR structural study of the interaction between a small‐molecule optical probe (DAOTA‐M2) and a G‐quadruplex from the promoter region of the c‐myc oncogene revealed that they interact at 1:2 binding stoichiometry. NMR‐restrained structural calculations show that binding of DAOTA‐M2 occurs mainly through π–π stacking between the polyaromatic core of the ligand and guanine residues of the outer G‐quartets. Interestingly, the binding affinities of DAOTA‐M2 differ by a factor of two for the outer G‐quartets of the unimolecular parallel G‐quadruplex under study. Unrestrained MD calculations indicate that DAOTA‐M2 displays significant dynamic behavior when stacked on a G‐quartet plane. These studies provide molecular guidelines for the design of triangulenium derivatives that can be used as optical probes for G‐quadruplexes.     

Structural basis of telomeric RNA quadruplex--acridine ligand recognition

Human telomeric DNA is now known to be transcribed into noncoding RNA sequences, termed TERRA. These sequences, which are believed to play roles in the regulation of telomere function, can form higher-order quadruplex structures and may themselves be the target of therapeutic intervention. The crystal structure of a TERRA quadruplex-acridine small-molecule complex at a resolution of 2.60 ?, is reported here and contrasts remarkably with the structure of the analogous DNA quadruplex complex. The bimolecular RNA complex has a parallel-stranded topology with propeller-like UUA loops. These loops are held in particular conformations by multiple hydrogen bonds involving the O2' hydroxyl groups of the ribonucleotide sugars and play an active role in binding the acridine molecules to the RNA quadruplex. By contrast, the analogous DNA quadruplex complex has simpler 1:1 acridine binding, with no loop involvement. There are significant loop conformational changes in the RNA quadruplex compared to the native TERRA quadruplex (Collie, G. W.; Haider, S. M.; Neidle, S.; Parkinson, G. N. Nucleic Acids Res. 2010, 38, 5569 - 5580), which have implications for the future design of small molecules targeting TERRA quadruplexes, and RNA quadruplexes more generally.     

Quadruplex-based molecular beacons as tunable DNA probes

Molecular beacons (MBs) are fluorescent nucleic acid probes with a hairpin-shaped structure in which the 5' and 3' ends are self-complementary. Due to a change in their emissive properties upon recognition with complementary sequences, MBs allow the diagnosis of single-stranded DNA or RNA with high mismatch discrimination, in vitro and in vivo. Whereas the stems of MB hairpins usually rely on the formation of a Watson-Crick duplex, we demonstrate in this report that the preceding structure can be replaced by a G-quadruplex motif (G4). Intramolecular quadruplexes may still be formed with a central loop composed of 12 to 21 bases, therefore extending the sequence repertoire of quadruplex formation. G4-MB can efficiently be used for oligonucleotide discrimination: in the presence of a complementary sequence, the central loop hybridizes and forms a duplex that causes opening of the quadruplex stem. The corresponding G4-MB unfolding can be detected by a change in its fluorescence emission. We discuss the thermodynamic and kinetic opportunities that are provided by using G4-MB instead of traditional MB. In particular, the intrinsic feature of the quadruplex motif facilitates the design of functional molecular beacons by independently varying the concentration of monovalent or divalent cations in the medium.     

Effect of loop orientation on quadruplex-TMPyP4 interaction

G-quadruplexes are believed to be potential targets for therapeutic intervention and this has resulted in designing of various quadruplex interacting ligands. Moreover, reports about existence of quadruplex forming sequences across the genome have propelled greater interest in understanding their interaction with small molecules. An intramolecular quadruplex sequence can adopt different conformations, owing to different orientation of loops in the structure. The differences in the loop orientation can affect their molecular recognition. Herein, we have studied the interaction of 5,10,15,20-tetrakis(1-methyl-4-pyridyl)-21H, 23H-porphine (TMPyP4), a well-known G quadruplex binding ligand with three DNA quadruplexes differing in loop orientations. Results obtained from UV, ITC, and SPR studies have coherently revealed that the TMPyP4 molecule shows preferential binding to parallel G-quadruplex ( c-myc and c-kit) over its antiparallel counterpart (human telomeric). The binding affinity for parallel quadruplex was (10(7)) 1 order of magnitude higher than that for antiparallel DNA quadruplex (10 ). The study shows two binding modes, stronger binding (10(7)) of TMPyP4 involving end stacking and a weaker external binding (10 ), while TMPyP4 shows only one binding mode with duplex with a binding affinity of the order of 10(6). Overall, the study emphasizes that differences in the loop orientation give rise to different conformations of quadruplex, which in turn govern its binding to small molecules, and thereby play a pivotal role in molecular recognition.     

The stability of intramolecular DNA quadruplexes with extended loops forming inter- and intra-loop duplexes

We have examined the formation of intramolecular quadruplex DNA structures in which the loops have been extended so as to generate short DNA duplexes. Fluorescence melting and DNase I cleavage studies show that duplexes can be formed within each loop, but that duplexes between the loops are not stable.     

Modeling and molecular dynamics simulation of the human gonadotropin-releasing hormone receptor in a lipid bilayer

In the present study, a model for the human gonadotropin-releasing hormone receptor embedded in an explicit lipid bilayer was developed. The final conformation was obtained by extensive molecular dynamics simulations of a homology model based on the bovine rhodopsin crystal structure. The analysis of the receptor structure allowed us to detect a number of specific contacts between different amino acid residues, as well as water- and lipid-mediated interactions. These interactions were stable in six additional independent 35 ns long simulations at 310 and 323 K, which used the refined model as the starting structure. All loops, particularly the extracellular loop 2 and the intracellular loop 3, exhibited high fluctuations, whereas the transmembrane helices were more static. Although other models of this receptor have been previously developed, none of them have been subjected to extensive molecular dynamics simulations, and no other three-dimensional structure is publicly available. Our results suggest that the presence of ions as well as explicit solvent and lipid molecules are critical for the structure of membrane protein models, and that molecular dynamics simulations are certainly useful for their refinement.     

A tale of two ions: the conformational landscapes of bis(trifluoromethanesulfonyl)amide and N,N-dialkylpyrrolidinium

The conformational landscapes of two commonly used ionic liquid ions, the anion bis(trifluoromethanesulfonyl)amide (Ntf2) and the cations N-propyl- and N-butyl-N-methylpyrrolidinium, were investigated using data obtained from Raman spectroscopy, molecular dynamics, and ab initio techniques. In the case of Ntf2, the plotting of three-dimensional potential energy surfaces (PES) and the corresponding molecular dynamics (MD) simulations confirmed the existence of two stable isomers (each existing as a pair of enantiomers) and evidenced the nature of the anion as a flexible, albeit hindered, molecule capable of interconversion between conformers in the liquid state, a result confirmed by the Raman data. In the case of the N,N-dialkylpyrrolidinium cations, the PES show a much more limited conformational behavior of the pyrrolidinium ring (pseudorotation). Nevertheless, such pseudorotation produces two stable isomers with the propyl and butyl side chains in completely different positions (axial-envelope and equatorial-envelope conformations). This result was also confirmed by Raman spectra analyses and MD simulations in the liquid phase. The implications of the conformational behavior of the two types of ions are discussed in terms of the solvation properties of the corresponding ionic liquids.     

Molecular Simulation of Interlayer Structure and Dynamics in 12.4 Å Cs-Smectite Hydrates

A detailed understanding of hydrated Cs-smectites is necessary to predict the permeability of clay liners to radiocesium cations at nuclear waste containment facilities. Monte Carlo (MC) and molecular dynamics (MD) modeling techniques were applied to three representative Cs-smectites to interpret a variety of experimental data on interlayer structure and dynamics. Spectroscopic and surface chemistry methods that attempt to differentiate interlayer water from water residing in micropores have provided data suggesting that, in stable 12.4 ? Cs-smectite hydrates, the interlamellar water content is less than one-half monolayer. Convergence profiles in MC simulations predicted stable hydrates at interlayer water contents of 1/3 or possibly 2/3 water monolayer. Radial distribution functions and coordination number data illustrated the ability of Cs(+) to organize water molecules into partial hydration shells and displayed the distortions of water structure induced by the clay surface. Molecular dynamics simulations of the MC-stable Cs-smectites revealed interlayer Cs(+) to be strongly bound as innersphere surface complexes, in agreement with published bulk diffusion coefficients. The strongly adsorbed Cs(+) can be associated with one of the species identified in (133)Cs NMR spectroscopic studies of hydrated Cs-smectites. These cations typically exhibited jump diffusion, whereas continuous diffusion of H(2)O occurred. Copyright 2001 Academic Press.