Unusual DNA conformations

Impact of geometry optimization on base-base stacking interaction energies in the canonical A- and B-forms of DNA

Base stacking is known to make an important contribution to the stability of DNA and RNA, and accordingly, significant efforts are ongoing to calculate stacking energies using ab initio quantum mechanical methods. To date, impressive improvements have been made in the model chemistries used to perform stacking energy calculations, including extensions that include robust treatments of electron correlation with extended basis sets, as required to treat interactions where dispersion makes a significant contribution. However, those efforts typically use rigid monomer geometries when calculating the interaction energies. To overcome this, in the present work, we describe a novel internal coordinate definition that allows the relative, intermolecular orientation of stacked base monomers to be constrained during geometry optimizations while allowing full optimization of the intramolecular degrees of freedom. Use of the novel reference frame to calculate the impact of full geometry optimization versus constraining the bases to be planar on base monomer stacking energies, combined with density-fitted, spin-component scaling MP2 treatment of electron correlation, shows that full optimization makes the average stacking energy more favorable by -3.4 and -1.5 kcal/mol for the canonical A and B conformations of the 16 5' to 3' base stacked monomers. Thus, treatment of geometry optimization impacts the stacking energies to an extent similar to or greater than the impact of current state of the art increases in the rigor of the model chemistry itself used to treat base stacking. Results also indicate that stacking favors the B-form of DNA, though the average difference versus the A-form decreases from -2.6 to -0.6 kcal/mol when the intramolecular geometry is allowed to fully relax. However, stacking involving cytosine is shown to favor the A-form of DNA, with that contribution generally larger in the fully optimized bases. The present results show the importance of allowing geometry optimization, as well as properly treating the appropriate model chemistry, in studies of nucleic acid base stacking.

A quantitative measure of chirality inside nucleic acid databank

We show the capability of a chirality index (Pietropaolo et al., Proteins 2008;70:667-677) to investigate nucleic acid structures because of its high sensitivity to helical conformations. By analyzing selected structures of DNA and RNA, we have found that sequences rich in cytosine and guanine have a tendency to left-handed chirality, in contrast to regions rich in adenine or thymine which show strong negative, right-handed, chirality values. We also analyze RNA structures, where specific loops and hairpin motifs are characterized by a well-defined chirality value. We find that in nucleosome the chirality is exalted, whereas in ribosome it is reduced. Our results illustrate the sensitivity of this descriptor for nucleic acid conformations.

Effect of gamma-hydroxypropano deoxyguanosine, the major acrolein-derived adduct, on monomolecular quadruplex structure of telomeric repeat d(TTAGGG)(4)

The three oligodeoxyribonucleotides (ODNs) a-c, having the telomeric repeat d(TTAGGG)(4) sequence and incorporating gamma-hydroxypropano deoxyguanosine at different positions, were synthesized. Gel electrophoresis and CD analyses indicated that the ODNs assume monomolecular quadruplex structures in Na+ and in K+ buffers. The T(m) values, obtained by CD melting experiments, showed that the presence of the acrolein-dG adduct into the ODN b decreases the thermal stability of the monomolecular quadruplex structure in Na+ solution, whereas for a and c no significant effect could be detected in the same experimental conditions. On the contrary, all ODNs a-d show the same behaviour in K+ buffer. These findings are briefly discussed.

How sequence defines structure: a crystallographic map of DNA structure and conformation

The fundamental question of how sequence defines conformation is explicitly answered if the structures of all possible sequences of a macromolecule are determined. We present here a crystallographic screen of all permutations of the inverted repeat DNA sequence d(CCnnnN6N7N8GG), where N6, N7, and N8 are any of the four naturally occurring nucleotides. At this point, 63 of the 64 possible permutations have been crystallized from a defined set of solutions. When combined with previous work, we have assembled a data set of 37 single-crystal structures from 29 of the sequences in this motif, representing three structural classes of DNA (B-DNA, A-DNA, and four-stranded Holliday junctions). This data set includes a unique set of amphimorphic sequence, those that crystallize in two different conformations and serve to bridge the three structural phases. We have thus constructed a map of DNA structures that can be walked through in single nucleotide steps. Finally, the resulting data set allows us to dissect in detail the stabilization of and conformational variations within structural classes and identify significant conformational deviations within a particular structural class that result from sequence rather than crystal or crystallization effects.

Empirical force fields for biological macromolecules: overview and issues

Empirical force field-based studies of biological macromolecules are becoming a common tool for investigating their structure-activity relationships at an atomic level of detail. Such studies facilitate interpretation of experimental data and allow for information not readily accessible to experimental methods to be obtained. A large part of the success of empirical force field-based methods is the quality of the force fields combined with the algorithmic advances that allow for more accurate reproduction of experimental observables. Presented is an overview of the issues associated with the development and application of empirical force fields to biomolecular systems. This is followed by a summary of the force fields commonly applied to the different classes of biomolecules; proteins, nucleic acids, lipids, and carbohydrates. In addition, issues associated with computational studies on "heterogeneous" biomolecular systems and the transferability of force fields to a wide range of organic molecules of pharmacological interest are discussed.

Selective recognition of G-qQuadruplex telomeric DNA by a bis(quinacridine) macrocycle

The interaction of G-quadruplex DNA with the macrocyclic compound BOQ1, which possesses two dibenzophenanthroline (quinacridine) subunits, has been investigated by a variety of methods. The oligonucleotide 5'-A(GGGT(2)A)(3)G(3), which mimics the human telomeric repeat sequence and forms an intramolecular quadruplex, was used as one model system. Equilibrium binding constants measured by biosensor surface plasmon resonance (SPR) methods indicate a high affinity of the macrocycle for the quadruplex conformation (K > 1 x 10(7) M(-)(1)) with two equivalent binding sites. The affinity of BOQ1 for DNA duplexes is at least 1 order of magnitude lower. In addition, the macrocycle is more selective than the monomeric control compound (MOQ2), which is not able to discriminate between the two DNA structures (K(duplex) approximately K(quadruplex) approximately 10(6) M(-)(1)). Strong binding of BOQ1 to G4 DNA sequences was confirmed by fluorometric titrations with a tetraplex-forming oligonucleotide. Competition dialysis experiments with a panel of different DNA structures, from single strands to quadruplexes, clearly established the quadruplex binding specificity of BOQ1. Fluorescence resonance energy transfer (FRET) T(m) experiments with a doubly labeled oligonucleotide also revealed a strong stabilization of the G4 conformation in the presence of BOQ1 (DeltaT(m) = +28 degrees C). This DeltaT(m) value is one of the highest values measured for a G-quadruplex ligand and is significantly higher than observed for the monomer control compounds (DeltaT(m) = +10-12 degrees C). Gel mobility shift assays indicated that the macrocycle efficiently induces the formation of G-tetraplexes. Strong inhibition of telomerase was observed in the submicromolar range (IC(50) = 0.13 microM). These results indicate that macrocycles represent an exciting new development opportunity for targeting DNA quadruplexes.

Multivariate curve resolution: a powerful tool for the analysis of conformational transitions in nucleic acids

A successful application is reported of the multivariate curve resolution alternating least-squares method (MCR-ALS) for the analysis of nucleic acid melting and salt-induced transitions. Under conditions where several structures co-exist in a conformational equilibrium, MCR-ALS analysis of the UV and circular dichroism (CD) spectra at different temperatures, ionic strength and oligonucleotide concentration allows for the resolution of concentration profiles and pure spectra of the different species. The methodology is illustrated by the case of the cyclic oligonucleotide d<pTGCTCGCT>. The melting transition of this molecule at different oligonucleotide concentrations was studied at 0, 2 and 10 mM MgCl2 by UV and CD spectroscopy. In addition, salt titration experiments were carried out at 21.0 and 54.0 degrees C. The MCR-ALS analysis indicates that three different conformations of this molecule co-exist in solution. In agreement with previous NMR studies, these conformations were assigned to a monomeric dumbbell-like structure, a dimeric four-stranded conformation and a disordered (random coil) structure. The MCR-ALS methodology allows for a detailed analysis of how this equilibrium is affected by temperature, salt and oligonucleotide concentration.

ADAPT: a molecular mechanics approach for studying the structural properties of long DNA sequences

We describe an original approach to determining sequence-structure relationships for DNA. This approach, termed ADAPT, combines all-atom molecular mechanics with a multicopy algorithm to build nucleotides that contain all four standard bases in variable proportions. These nucleotides enable us to search very rapidly for base sequences that energetically favor chosen types of DNA deformation or chosen DNA-protein or DNA-ligand interactions. Sequences satisfying the chosen criteria can be found by energy minimization, combinatorial sequence searching, or genome scanning, in a manner similar to the threading approaches developed for protein structure prediction. In the latter case, we are able to analyze roughly 2000 base pairs per second. Applications of the method to DNA allomorphic transitions, DNA deformation, and specific DNA interactions are presented.

Intrinsic conformational energetics associated with the glycosyl torsion in DNA: a quantum mechanical study

The glycosyl torsion (chi) in nucleic acids has long been recognized to be a major determinant of their conformational properties. chi torsional energetics were systematically mapped in deoxyribonucleosides using high-level quantum mechanical methods, for north and south sugar puckers and with gamma in the g(+) and trans conformations. In all cases, the syn conformation is found higher in energy than the anti. When gamma is changed from g(+) to trans, the anti orientation of the base is strongly destabilized, and the energy difference and barrier between anti and syn are significantly decreased. The barrier between anti and syn in deoxyribonucleosides is found to be less than 10 kcal/mol and tends to be lower with purines than with pyrimidines. With gamma = g(+)/chi = anti, a south sugar yields a significantly broader energy well than a north sugar with no energy barrier between chi values typical of A or B DNA. Contrary to the prevailing view, the syn orientation is not more stable with south puckers than with north puckers. The syn conformation is significantly more energetically accessible with guanine than with adenine in 5-nucleotides but not in nucleosides. Analysis of nucleic acid crystal structures shows that gamma = trans/chi = anti is a minor but not negligible conformation. Overall, chi appears to be a very malleable structural parameter with the experimental chi distributions reflecting, to a large extent, the associated intrinsic torsional energetics.

Metal-modified nucleobase sextet: joining four linear metal fragments (trans-a2PtII) and six model nucleobases to an exceedingly stable entity

Crosslinking of three different model nucleobases (9-ethyladenine, 9-EtA; 9-ethylguanine, 9-EtGH; 1-methyluracil, 1-MeU) by two linear trans-aPtII (a = NH3 or CH3NH2) entities leads to a flat metal-modified base triplet, trans,trans-[(NH3)2Pt(1-MeU-N3)(mu-9-EtA-N7,N1)Pt(CH3NH2)2(9-EtGH-N7)]3+ (4b). Upon hemideprotonation of the 9-ethylguanine base at the N1 position. 4b spontaneously dimerizes to the metalated nucleobase sextet 5, [(4b)(triple bond)(4b-H)]5+. In this dimeric structure a neutral and an anionic guanine ligand, which are complementary to each other, are joined through three H bonds and additionally by two H bonds between guanine and uracil nucleobases. Four additional interbase H bonds maintain the approximate coplanarity of all six bases. The two base triplets form an exceedingly stable entity (KD = 500 +/- 150 M(-1) in DMSO), which is unprecedented in nucleobase chemistry. The precursor of 4b and several related complexes are described and their structures and solution properties are reported.

Twisting and stretching single DNA molecules

The elastic properties of DNA are essential for its biological function. They control its bending and twisting as well as the induction of structural modifications in the molecule. These can affect its interaction with the cell machinery. The response of a single DNA molecule to a mechanical stress can be precisely determined in single-molecule experiments which give access to an accurate measurement of the elastic parameters of DNA.

Stress-induced structural transitions in DNA and proteins

The ability to manipulate, stretch and twist biomolecules opens the way to an understanding of their structural transitions. We review some of the recently discovered stress-induced structural transitions in DNA as well as the application of single molecule manipulation techniques to DNA unzipping and to the study of protein folding/unfolding transitions.

Fluorescence energy transfer analysis of DNA structures containing several bulges and their interaction with CAP

DNA molecules with three bulges separated by double-stranded helical sections of B-DNA were constructed to be used as substrates for DNA-protein binding assays. Fluorescence resonance energy transfer (FRET) between dye molecules attached to the 5'-ends of the DNA molecules is used to monitor the protein binding. The A5 bulge, which consists of five unpaired adenine nucleotides, alters the direction of the helical axis by approximately 80 to 90 at every bulge site. Computer molecular modeling facilitated a pre-selection of suitable helix lengths that bring the labeled ends of the three-bulge DNA molecules (60 to 70 base-pairs long) into close proximity. The FRET experiments verified that the labeled ends of the helices of these long molecules were indeed close. A series of FRET experiments was carried out with two A5 and two A7 bulge molecules. The relative positions of the bulges were varied along the central helical DNA sequence (between the bulges) in order to determine the relative angular juxtapositions of the outlying helical arms flanking the central helical region. The global structural features of the DNA molecules are manifested in the FRET data. The FRET experiments, especially those of the two-bulge series, could be interpreted remarkably well with molecular models based on the NMR structure of the A5 bulge. These models assume that the DNA molecules do not undergo large torsional conformational fluctuations at the bulge sites. The magnitude of the FRET efficiency attests to a relatively rigid structure for many of the long 5'-end-labeled molecules. The changes in the FRET efficiency of three-bulge structures containing the specific binding sequence of the catabolite activator protein (CAP) demonstrated significant deformation of the DNA upon binding of CAP. No direct interaction of CAP with the dyes was observed.

Modeling DNA deformations

Recent developments have been made in modeling double-helical DNA at four levels of three-dimensional structure: the all-atom level, whereby an oligonucleotide duplex is surrounded by a shroud of solvent molecules; the base-pair level, with explicit backbone atoms; the mesoscopic level, that is, a few hundred base pairs, with the local duplex conformation described by knowledge-based harmonic energy functions; and the scale of several thousand nucleotides, with the duplex described as an ideal elastic rod. Predictions of the sequence-dependent bending and twisting of the double helix, as well as solvent- and force-induced B-->A and over-stretching conformational transitions, are compared with experimental data. These subtle conformational changes are critical to the functioning of the double helix, including its packaging in the close confines of the cell, the mutual fit of DNA and protein in nucleoprotein complexes, and the effective recognition of base pairs in recombination and transcription.

Nucleic acids: theory and computer simulation, Y2K

Molecular dynamics simulations on DNA and RNA that include solvent are now being performed under realistic environmental conditions of water activity and salt. Improvements to force-fields and treatments of long-range interactions have significantly increased the reliability of simulations. New studies of sequence effects, axis bending, solvation and conformational transitions have appeared.

Fold-back tetraplex DNA species in DNase I-resistant DNA isolated from HeLa cells

A DNase I-resistant DNA species has been isolated and purified from HeLa cells by gel electrophoresis. Our studies indicate that the DNase I-resistant DNA species was about 40-60 bp fragment sizes responding to double-strand DNA marker and has higher guanine content. The image of AFM showed that this species has been assumed to be tetraplex structure according to its apparent width and height. Its CD, UV spectrum also exhibited characteristics similar to some tetraplex structure, which was different from the standard duplex DNA. 32P-labeled probes (TTAGGG)4 and 5'-TGGGGAGGGTGGGGAGGGTGGGGAAGG-3' could be hybridized to purified DNase I-resistant species. All results suggest that the DNase I-resistant DNA species have at least two components, which adopt an intrastrand fold-back DNA tetraplex. Their sequences were similar to human telomere and human c-myc locus (NHE), respectively.

Association of double-stranded DNA fragments into multistranded DNA structures

We have previously observed that double-stranded DNA fragments containing a tract of the tandemly repeated sequence poly(CA). poly(TG) can associate in vitro to form stable complexes of low electrophoretic mobility, which are recognized with high specificity by proteins HMG1 and HMG2. The formation of such complexes has since been observed to depend on interactions of DNA with polypropylene surfaces, with the suggestion that the formation of low mobility complexes might be the result of strand dissociation followed by misaligned reassociation of the repetitive sequences. The data presented here show that at high ionic strength the interactions of DNA with polypropylene are sufficiently strong for DNA to remain bound to the polypropylene surface, which suggests that DNA might also be involved in interactions with hydrophobic molecules in vivo. Under such conditions, low-mobility complexes are found only in the material adsorbed to the polypropylene surface, and all DNA fragments are able to form low-mobility structures, whether or not they contain repetitive sequences. Preventing the separation of strands by ligating hairpin loop oligonucleotides at both ends of the fragments does not prevent the formation of low-mobility complexes. Our results suggest two different pathways for the formation of complexes. In the first, dissociation is followed by misaligned reassociation of repetitive sequences, yielding duplexes with single-stranded end regions that associate to form multimeric complexes. In the second, repetitive as well as nonrepetitive DNA molecules bound to polypropylene adopt a conformation with locally unwound regions, which allows interactions between neighboring duplexes adsorbed on the surface, resulting in the formation of low-mobility complexes.

Collective variable modelling of nucleic acids

Collective variable models continue to contribute to our knowledge of nucleic acids. The past year has seen considerable progress both in modelling sequence-dependent effects on nucleic acid conformation and in understanding how proteins or external stresses influence nucleic acid structure. Algorithmic developments have also allowed collective models to be applied to studies of thermal fluctuations and dynamics. For larger systems, models with varying degrees of resolution are being refined and applied to nucleic acids containing hundreds or thousands of nucleotides.

Spectroscopic comparison of different DNA structures formed by oligonucleotides

Six different nucleic acid structures including duplex, triplex and quadruplex are formed by oligonucleotides. Their structural properties are studied in detail by four spectroscopic techniques, i.e. CD, UV, NMR and fluorescence. Results are: CD Spectra: The common characteristics is a negative band at 240 nm, and the spectra are different from each other in the range 260-300 nm. Many factors such as chain direction, sugar puckering, orientation of the glycosyl bond, base stacking and sequence can effect their conformation and then show diversity and complexity in the spectra. UV Spectra: The UV spectra of all forms are quite similar, all of them exhibit a sharp positive peak around 210 nm and a broad positive band in the region of 240-280 nm. Although the bands are different in absorbance, the spectra are not characteristic enough to distinguish these forms. In addition, their thermal denaturation is also observed by UV spectrum, different melting curves and points are shown and some thermodynamic information is provided. NMR Spectra: Since the G residues in the six samples all participate in hydrogen bond, the imino proton can not exchange with the solvent freely so as to allow an observable resonance to arise. The resonance number and chemical shift will vary with the change in base-pairing number and mode as well as the whole geometry of its molecule. Fluorescence Spectra: The interaction mechanisms between EB and these structures are different. B type duplex and triplex adopt an intercalative mode in which the efficiency of energy transfer is relatively high and the fluorescence of EB can not be quenched easily. While for the parallel duplex, outside binding is predominant in which energy transfer can hardly happen and most of its fluorescence can be quenched. As for the quadruplex, groove binding is possible, so the efficiency of energy transfer is higher than that in outside binding, but lower than that in intercalative binding, and fluorescence is quenched partly.

Progressive DNA bending is made possible by gradual changes in the torsion angle of the glycosyl bond

Structural comparisons have led to the suggestion that the conformational rearrangement that would be required to change A-DNA into the TA-DNA form of DNA observed in the complex with the TATA box binding protein (TBP) could be completed by modifying only the value of the glycosyl bond chi by approximately 45 degrees. The lack of a high number of crystal structures of this type makes it difficult to conclude whether a smooth transition from A-DNA to TA-DNA can occur without disrupting at any point either the Watson-Crick base pairing or the A-DNA conformation of the backbone. To explore the possibility of such a smooth transition, constrained molecular dynamics simulations were carried out for the double-stranded dodecamer d(GGTATATAAAAC), in which a transition from A-DNA to TA-DNA was induced by modifying only the chi angle values. The results demonstrate the feasibility of a continuous path in the A-DNA to TA-DNA transition. Varying extents of DNA curvature are also attainable, by maintaining the A-DNA backbone structure and Watson-Crick hydrogen bonding while changing the chi angle value smoothly from that in A-DNA to one corresponding to B-DNA.