Extended molecular dynamics of a c-kit promoter quadruplex

The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 μs long, ∼50 μs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal μs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.     

Molecular dynamics studies of the 3D structure and planar ligand binding of a quadruplex dimer

G-rich sequences can fold into a four-stranded structure called a G-quadruplex, and sequences with short loops are able to aggregate to form stable quadruplex multimers. Few studies have characterized the properties of this variety of quadruplex multimers. Using molecular modeling and molecular dynamics simulations, the present study investigated a dimeric G-quadruplex structure formed from a simple sequence of d(GGGTGGGTGGGTGGGT) (G1), and its interactions with a planar ligand of a perylene derivative (Tel03). A series of analytical methods, including free energy calculations and principal components analysis (PCA), was used. The results show that a dimer structure with stacked parallel monomer structures is maintained well during the entire simulation. Tel03 can bind to the dimer efficiently through end stacking, and the binding mode of the ligand stacked with the 3′-terminal thymine base is most favorable. PCA showed that the dominant motions in the free dimer occur on the loop regions, and the presence of the ligand reduces the flexibility of the loops. Our investigation will assist in understanding the geometric structure of stacked G-quadruplex multimers and may be helpful as a platform for rational drug design.     

Conformational dynamics of the human propeller telomeric DNA quadruplex on a microsecond time scale

The human telomeric DNA sequence with four repeats can fold into a parallel-stranded propeller-type topology. NMR structures solved under molecular crowding experiments correlate with the crystal structures found with crystal-packing interactions that are effectively equivalent to molecular crowding. This topology has been used for rationalization of ligand design and occurs experimentally in a number of complexes with a diversity of ligands, at least in the crystalline state. Although G-quartet stems have been well characterized, the interactions of the TTA loop with the G-quartets are much less defined. To better understand the conformational variability and structural dynamics of the propeller-type topology, we performed molecular dynamics simulations in explicit solvent up to 1.5 μs. The analysis provides a detailed atomistic account of the dynamic nature of the TTA loops highlighting their interactions with the G-quartets including formation of an A:A base pair, triad, pentad and hexad. The results present a threshold in quadruplex simulations, with regards to understanding the flexible nature of the sugar-phosphate backbone in formation of unusual architecture within the topology. Furthermore, this study stresses the importance of simulation time in sampling conformational space for this topology.     

Water spines and networks in G-quadruplex structures

Quadruplex DNAs can fold into a variety of distinct topologies, depending in part on loop types and orientations of individual strands, as shown by high-resolution crystal and NMR structures. Crystal structures also show associated water molecules. We report here on an analysis of the hydration arrangements around selected folded quadruplex DNAs, which has revealed several prominent features that re-occur in related structures. Many of the primary-sphere water molecules are found in the grooves and loop regions of these structures. At least one groove in anti-parallel and hybrid quadruplex structures is long and narrow and contains an extensive spine of linked primary-sphere water molecules. This spine is analogous to but fundamentally distinct from the well-characterized spine observed in the minor groove of A/T-rich duplex DNA, in that every water molecule in the continuous quadruplex spines makes a direct hydrogen bond contact with groove atoms, principally phosphate oxygen atoms lining groove walls and guanine base nitrogen atoms on the groove floor. By contrast, parallel quadruplexes do not have extended grooves, but primary-sphere water molecules still cluster in them and are especially associated with the loops, helping to stabilize loop conformations.     

Distinctive features in the structure and dynamics of the DNA repeat sequence GGCGGG

G-rich DNA has been known to form a variety of folded and multistranded structures, with even single base modifications causing important structural changes. But, very little is known about the dynamic characteristics of the structures, which may play crucial roles in facilitating the structural transitions. In this background, we report here NMR investigations on the structure and dynamics of a DNA repeat sequence GGCGGG in aqueous solution containing Na+ ions at neutral pH. The chosen sequence d-TGGCGGGT forms a parallel quadruplex with a C-tetrad in the middle, formed by symmetrical pairing of four Cs in a plane via NH2-O2 H-bonds. 13C relaxation measurements at natural abundance for C' sugar carbons provided valuable insight into the sequence specific dynamism of G and C-tetrads in the quadruplex. The C4 tetrad seems to introduce high conformational dynamism at milli- to micro-second time scale in the quadruplex. Concomitantly, there is a decrease in the pico-second time scale dynamics. Interestingly, these effects are seen more prominently at the G-tetrads on the 3' end of C-tetrad than on its 5' end. These observations would have important implications for the roles the tetrads may play in many biological functions.     

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)₃] 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.     

Relative stability of G‐quadruplex structures: Interactions between the human Bcl2 promoter region and derivatives of carbazole and diphenylamine

The bcl2 promoter region forms a G‐quadruplex structure, which is a crucial target for anticancer drug development. In this study, we provide theoretical predictions of the stability of different G‐quadruplex folds of the 23‐mer bcl2 promoter region and G‐quadruplex ligand. We take into account the whole G‐quadruplex structure, including bound‐cations and solvent effects, in order to compute the ligand binding free energy using molecular dynamics simulation. Two series of the carbazole and diphenylamine derivatives are used to screen for the most potent drug in terms of stabilization. The energy analysis identifies the predominant energy components affecting the stability of the various different G‐quadruplex folds. The energy associated with the stability of the G‐quadruplex‐K⁺ structures obtained displays good correlation with experimental T_m measurements. We found that loop orientation has an intrinsic influence on G‐quadruplex stability and that the basket structure is the most stable. Furthermore, parallel loops are the most effective drug binding site. Our studies also demonstrate that rigidity and planarity are the key structural elements of a drug that stabilizes the G‐quadruplex structure. BMVC‐4 is the most potential G‐quadruplex ligand. This approach demonstrates significant promise and should benefit drug design. © 2014 Wiley Periodicals, Inc. Biopolymers 101: 1038–1050, 2014.     

EFFECTS OF A 8-OXOADENOSINE INCORPORATION ON QUADRUPLEX STRUCTURES: THERMAL STABILITIES AND STRUCTURAL STUDIES

The effects of incorporation of 8-oxoadenosine in two different truncations of human telomeric sequence forming quadruplex structures are reported. In order to characterise their structures, a combination of NMR and UV spectroscopy and computational techniques were used. Both oligonucleotides have been found to form fourfold symmetric quadruplex structures. As a tautomeric equilibrium between keto and enol forms of 8-oxoadenosine may establish in solution and intrinsic stabilities effects, such as internal H-bonds, for example, may determine the predominance of some particular tautomer, molecular modelling studies were performed on quadruplex structures containing both the tautomeric forms. Both molecules resulted to be thermally less stable than the natural.     

Effects of a 8-oxoadenosine incorporation on quadruplex structures: thermal stabilities and structural studies

The effects of incorporation of 8-oxoadenosine in two different truncations of human telomeric sequence forming quadruplex structures are reported. In order to characterise their structures, a combination of NMR and UV spectroscopy and computational techniques were used. Both oligonucleotides have been found to form fourfold symmetric quadruplex structures. As a tautomeric equilibrium between keto and enol forms of 8-oxoadenosine may establish in solution and intrinsic stabilities effects, such as internal H-bonds, for example, may determine the predominance of some particular tautomer, molecular modelling studies were performed on quadruplex structures containing both the tautomeric forms. Both molecules resulted to be thermally less stable than the natural.     

Molecular dynamics and principal components analysis of human telomeric quadruplex multimers

Guanine-rich DNA repeat sequences located at the terminal ends of chromosomal DNA can fold in a sequence-dependent manner into G-quadruplex structures, notably the terminal 150-200 nucleotides at the 3′ end, which occur as a single-stranded DNA overhang. The crystal structures of quadruplexes with two and four human telomeric repeats show an all-parallel-stranded topology that is readily capable of forming extended stacks of such quadruplex structures, with external TTA loops positioned to potentially interact with other macromolecules. This study reports on possible arrangements for these quadruplex dimers and tetramers, which can be formed from 8 or 16 telomeric DNA repeats, and on a methodology for modeling their interactions with small molecules. A series of computational methods including molecular dynamics, free energy calculations, and principal components analysis have been used to characterize the properties of these higher-order G-quadruplex dimers and tetramers with parallel-stranded topology. The results confirm the stability of the central G-tetrads, the individual quadruplexes, and the resulting multimers. Principal components analysis has been carried out to highlight the dominant motions in these G-quadruplex dimer and multimer structures. The TTA loop is the most flexible part of the model and the overall multimer quadruplex becoming more stable with the addition of further G-tetrads. The addition of a ligand to the model confirms the hypothesis that flat planar chromophores stabilize G-quadruplex structures by making them less flexible.     

Folding Topology of a Bimolecular DNA Quadruplex Containing a Stable Mini-hairpin Motif within the Diagonal Loop

We describe the NMR structural characterisation of a bimolecular anti-parallel DNA quadruplex d(G₃ACGTAGTG₃)₂ containing an autonomously stable mini-hairpin motif inserted within the diagonal loop. A folding topology is identified that is different from that observed for the analogous d(G₃T₄G₃)₂ dimer with the two structures differing in the relative orientation of the diagonal loops. This appears to reflect specific base stacking interactions at the quadruplex-duplex interface that are not present in the structure with the T₄-loop sequence. A truncated version of the bimolecular quadruplex d(G₂ACGTAGTG₂)₂, with only two core G-tetrads, is less stable and forms a heterogeneous mixture of three 2-fold symmetric quadruplexes with different loop arrangements. We demonstrate that the nature of the loop sequence, its ability to form autonomously stable structure, the relative stabilities of the hairpin loop and core quadruplex, and the ability to form favourable stacking interactions between these two motifs are important factors in controlling DNA G-quadruplex topology.     

NMR structure refinement and dynamics of the K+-[d(G3T4G3)]2 quadruplex via particle mesh Ewald molecular dynamics simulations.

The solution structure and dynamical properties of the potassium-stabilized, hairpin dimer quadruplex formed by the oligonucleotide d(G3T4G3) have been elucidated by a combination of high-resolution NMR and molecular dynamics simulations. Refinement calculations were carried out both in vacuo, without internally coordinated K+ cations, and in explicit water, with internally coordinated K+ cations. In the latter case, the electrostatic interactions were calculated using the particle mesh Ewald (PME) method. The NMR restraints indicate that the K+ quadruplex has a folding arrangement similar to that formed by the same oligonucleotide in the presence of sodium, but with significant local differences. Unlike the Na+ quadruplex, the thymine loops found in K+ exhibit considerable flexibility, and appear to interconvert between two preferred conformations. Furthermore, the NMR evidence points toward K+-stabilized guanine quartets of slightly larger diameter relative to the Na+-stabilized structure. The characteristics of the quartet stem are greatly affected by the modeling technique employed: caged cations alter the size and symmetry of the quartets, and explicit water molecules form hydration spines within the grooves. These results provide insight into those factors that determine the overall stability of hairpin dimer quadruplexes and the effects of different cations in modulating the relative stability of the dimeric hairpin and linear, four-stranded, quadruplex forms.     

Determinants of DNA quadruplex structural type: sequence and potassium binding

There are DNA sequences which adopt the same quadruplex structural type in the presence of sodium as in the presence of sodium and potassium. There are also sequences that appear to have a requirement for the presence of potassium for the adoption of a particular quadruplex structural type. Information about the basis for these potassium effects has been obtained by examining the structures of a set of DNAs with differing numbers of loop residues and different lengths of runs of dG residues in the presence of sodium alone and in the presence of potassium and sodium. On the basis of the results, obtained primarily via solution-state NMR, it appears that very small loops favor parallel stranded quartet structures which do not require the presence of potassium. DNAs with loops of two to four residues and runs of two dG residues can form quadruplex structures of the "edge" or "chair" type in the presence of potassium but not in the presence of sodium alone. When all of the loops contain four residues, a "crossover" or "basket" type structure can be formed in the presence of sodium as well as in the presence of sodium and potassium. Structures with runs of three or four dG residues and with loops from two to four residues can form basket or crossover type structures in the absence of potassium. The presence of a purine in a loop can block both potassium binding and formation of chair type structures. Modeling of the interactions of cations with these quadruplex structures indicates that the potassium ions required for chair type structures interact with a terminal quartet and residues in the adjacent loop.     

Crystal structure of a c-kit promoter quadruplex reveals the structural role of metal ions and water molecules in maintaining loop conformation

We report here the 1.62 Å crystal structure of an intramolecular quadruplex DNA formed from a sequence in the promoter region of the c-kit gene. This is the first reported crystal structure of a promoter quadruplex and the first observation of localized magnesium ions in a quadruplex structure. The structure reveals that potassium and magnesium ions have an unexpected yet significant structural role in stabilizing particular quadruplex loops and grooves that is distinct from but in addition to the role of potassium ions in the ion channel at the centre of all quadruplex structures. The analysis also shows how ions cluster together with structured water molecules to stabilize the quadruplex arrangement. This particular quadruplex has been previously studied by NMR methods, and the present X-ray structure is in accord with the earlier topology assignment. However, as well as the observations of potassium and magnesium ions, the crystal structure has revealed a highly significant difference in the dimensions of the large cleft in the structure, which is a plausible target for small molecules. This difference can be understood by the stabilizing role of structured water networks.     

A comparison of the properties and the solution structure for RNA and DNA quadruplexes which contain two GGAGG sequences joined with a tetranucleotide linker

We have determined solution structure of r(GGAGGUUUUGGAGG) (R14) by NMR; the RNA 14-mer forms an intra-strand parallel quadruplex with a G-tetrad and a hexad, in which a G-tetrad core is augmented by association of two A residues. The quadruplex further forms a dimer through stacking interaction between the hexads. In order to obtain insight into the difference between RNA and DNA quadruplexes, we synthesized the corresponding DNA 14-mer, d(GGAGGTTTTGGAGG) (D14), and examined its properties and structure by CD, gel electrophoresis, and NMR. K+ ions increased the thermal stability of both R14 and D14 structures. The binding affinity of K+ ions to R14 was much higher than that to D14. The CD and gel electrophoretic studies suggest that D14 forms a quadruplex entirely different from that of R14 in the presence of K+ ions; two molecules of D14 form a quadruplex with both antiparallel and parallel strand alignments and with diagonal loops at both ends of the stacked G-tetrads. The NMR study also gave results that are consistent with such structure: alternate glycosidic conformation, 5'G(syn)-G(anti)3', and characteristic chemical shift data observed for many quadruplexes containing diagonal TTTT loops.     

Molecular dynamics simulation for probing the flexibility of the 35 nucleotide SL1 sequence kissing complex from HIV-1Lai genomic RNA

The SL1 stem-loop located in the encapsidation domain is responsible for initiating the dimerisation of HIV-1 genomic RNA by means of a loop-loop interaction known as Kissing Complex (KC). The SL1 secondary structure has been predicted as a 35 nucleotides [K. G. Murti, M. Bondurant, and A. Tereba. J Virol 37, 411-419 (1981)] stem-loop composed of a 4 base pairs (bp) terminal duplex, a 4 nt asymmetrical internal loop, a 7 bp internal duplex, and a 9 nt apical loop. Several high resolution structures of the monomer and of KC of a 23 nt sequence containing only the internal duplex and the apical loop of SL1 are available in the literature. No experimental high resolution structure of the complete native SL1 sequence has been reported so far, either for the monomer or for KC. The asymmetrical internal loop has been described from NMR studies of different monomeric hairpin sequences, leading to divergent results, which suggests its high flexibility. In this work, we built a SL1(35) KC model which was submitted to a 31 ns molecular dynamics simulation (MD). Our results allows to describe the internal dynamics of SL1(35) KC and the differences of behavior of the different parts of the dimer. Thus, we could show the stability of the interactions between the two apical loops and of the terminal duplexes, the destabilization of the internal duplexes and the high flexibility of the asymmetrical internal loops.     

Crystal structure of the potassium form of an Oxytricha nova G-quadruplex

The crystal structures of the potassium-containing quadruplex formed from the Oxytricha nova sequence d(GGGGTTTTGGGG) are reported, in two space groups, the orthorhombic P2(1)2(1)2(1) and the trigonal P3(2)21, which diffract to 2.0 A and 1.49 A, respectively. The orthorhombic form contains two independent quadruplexes in the asymmetric unit, and the trigonal form contains one. All three of these quadruplexes adopt an identical fold, with two strands forming an antiparallel diagonal arrangement. This is identical with that observed previously in NMR studies of the native sodium and potassium forms, and a crystallographic analysis of it complexed with an O. nova protein. The present analysis demonstrates that the native structure is the same in solution and in the crystalline state and, moreover, that the nature of the counter-ion does not affect the overall fold of this quadruplex. The analysis corrects an earlier crystallographic study of this quadruplex. The conformation of the tetra-thymine loop is described in detail, which involves the third thymine base folding back to interact with the first thymine base. The water networks in the grooves and loops are described and, in particular, the ability of water molecules to form a continuous spine of hydration in the narrow groove is detailed. Each quadruplex has five potassium ions organised in a linear channel, with square antiprismatic coordination to each ion from oxygen atoms.     

Predictive modelling of topology and loop variations in dimeric DNA quadruplex structures

We have used a combination of simulated annealing (SA), molecular dynamics (MD) and locally enhanced sampling (LES) methods in order to predict the favourable topologies and loop conformations of dimeric DNA quadruplexes with T2 or T3 loops. This follows on from our previous MD simulation studies on the influence of loop lengths on the topology of intramolecular quadruplex structures [P. Hazel et al. (2004) J. Am. Chem. Soc., 126, 16405–16415], which provided results consistent with biophysical data. The recent crystal structures of d(G4T3G4)2 and d(G4BrUT2G4) (P. Hazel et al. (2006) J. Am. Chem. Soc., in press) and the NMR-determined topology of d(TG4T2G4T)2 [A.T. Phan et al. (2004) J. Mol. Biol., 338, 93–102] have been used in the present study for comparison with simulation results. These together with MM-PBSA free-energy calculations indicate that lateral T3 loops are favoured over diagonal loops, in accordance with the experimental structures; however, distinct loop conformations have been predicted to be favoured compared to those found experimentally. Several lateral and diagonal loop conformations have been found to be similar in energy. The simulations suggest an explanation for the distinct patterns of observed dimer topology for sequences with T3 and T2 loops, which depend on the loop lengths, rather than only on G-quartet stability.     

The effect on quadruplex stability of North-nucleoside derivatives in the loops of the thrombin-binding aptamer

Modified thrombin-binding aptamers (TBAs) carrying uridine (U), 2'-deoxy-2'-fluorouridine (FU) and North-methanocarbathymidine (NT) residues in the loop regions were synthesized and analyzed by UV thermal denaturation experiments and CD spectroscopy. The replacement of thymidines in the TGT loop by U and FU results in an increased stability of the antiparallel quadruplex structure described for the TBA while the presence of NT residues in the same positions destabilizes the antiparallel structure. The substitution of the thymidines in the TT loops for U, FU and NT induce a destabilization of the antiparallel quadruplex, indicating the crucial role of these positions. NMR studies on TBAs modified with uridines at the TGT loop also confirm the presence of the antiparallel quadruplex structure. Nevertheless, replacement of two Ts in the TT loops by uridine gives a more complex scenario in which the antiparallel quadruplex structure is present along with other partially unfolded species or aggregates.