Conformational Effects of DNA Stretching

DNA is an extensible molecule, and an extended conformation of DNA is involved in some biological processes. We have examined the effect of elongation stress on the conformational properties of DNA base pairs by conformational analysis. The calculations show that stretching does significantly affect the conformational properties and flexibilities of base pairs. In particular, we have found that the propeller twist in base pairs reverses its sign upon stretching. The energy profile analysis indicates that electrostatic interactions make a major contribution to the stabilization of the positive-propeller-twist configuration in stretched DNA. This stretching also results in a monotonic decrease in the helical twist angle, tending to unwind the double helix. Fluctuations in most variables initially increase upon stretching, because of unstacking of base pairs, but then the fluctuations decrease as DNA is stretched further, owing to the formation of specific interactions between base pairs induced by the positive propeller twist. Thus, the stretching of DNA has particularly significant effects upon DNA flexibility. These changes in both the conformation and flexibility of base pairs probably have a role in functional interactions with proteins.     

Conformational effects of DNA stretching

DNA is an extensible molecule, and an extended conformation of DNA is involved in some biological processes. We have examined the effect of elongation stress on the conformational properties of DNA base pairs by conformational analysis. The calculations show that stretching does significantly affect the conformational properties and flexibilities of base pairs. In particular, we have found that the propeller twist in base pairs reverses its sign upon stretching. The energy profile analysis indicates that electrostatic interactions make a major contribution to the stabilization of the positive-propeller-twist configuration in stretched DNA. This stretching also results in a monotonic decrease in the helical twist angle, tending to unwind the double helix. Fluctuations in most variables initially increase upon stretching, because of unstacking of base pairs, but then the fluctuations decrease as DNA is stretched further, owing to the formation of specific interactions between base pairs induced by the positive propeller twist. Thus, the stretching of DNA has particularly significant effects upon DNA flexibility. These changes in both the conformation and flexibility of base pairs probably have a role in functional interactions with proteins.     

Raman spectroscopy of supercoiled and nicked ColE1 plasmid

Native supercoiled and nicked ColE1 DNA were examined using laser Raman spectroscopy. ColE1 contains 6646 base pairs (bp) and, when supercoiled, approximately 47 negative supercoils. An analytical buoyant density gradient centrifugation technique developed by Burke and Bauer was scaled to preparative quantities, and used to isolate the supercoiled plasmid fraction from its nicked counterpart. This procedure allowed enriched fractions of the supercoiled plasmid to be extracted without the use of the optical contaminant ethidium bromide. The intensities of several Raman bands were altered between the spectra of the two topological forms. Notably absent were any changes in bands arising from cytosine and guanine vibrations. The observed changes are interpreted in terms of the polymorphic structures which have been observed in many DNA structural studies. The results of this study suggest that accommodation of supercoiling takes place chiefly in A-T base pairs and backbone moieties, without substantial modification of G-C base-pair structure. Premelting effects may account for the observed changes, including a slight shift to lower frequency of a band known to be responsive to base-pair disruption. Heteronomous ribose sugar pucker is evident in both supercoiled and nicked plasmid species. No gross conformational transitions were detected for native supercoiled DNA, and consequently, subtle rearrangements appear sufficient to absorb the supercoiling deformations.     

Influence of sequence on the conformation of the B-DNA helix.

We have tried to ascertain whether the variability found in the conformational features of the 10 base steps in B-DNA is mainly due to the flanking sequences or to interactions with the environment. From an analysis of the twist parameter of the base-pair steps available from crystals of oligonucleotides and protein/oligonucleotide complexes, we conclude that in most cases the flanking sequences show little influence: the conformation of a DNA region results from the combination of the independent intrinsic features of each base step (average conformation and intrinsic variability), modulated by their interactions with the environment. Only in some cases (YR steps, in particular CG and CA/TG) does it appear that flanking sequences have an influence on the conformation of the central base step. The values obtained allow an approximation to the parameters expected for repetitive DNA sequences. In particular, it is found that poly[d(AG/CT)] should have a strongly alternating conformation, in agreement with recently reported oligonucleotide structures.     

Atomic structures of excited state A–T Hoogsteen base pairs in duplex DNA by combining NMR relaxation dispersion,mutagenesis, and chemical shift calculations

NMR relaxation dispersion studies indicate that in canonical duplex DNA, Watson–Crick base pairs (bps) exist in dynamic equilibrium with short-lived low abundance excited state Hoogsteen bps. N1-methylated adenine (m¹A) and guanine (m¹G) are naturally occurring forms of damage that stabilize Hoogsteen bps in duplex DNA. NMR dynamic ensembles of DNA duplexes with m¹A–T Hoogsteen bps reveal significant changes in sugar pucker and backbone angles in and around the Hoogsteen bp, as well as kinking of the duplex towards the major groove. Whether these structural changes also occur upon forming excited state Hoogsteen bps in unmodified duplexes remains to be established because prior relaxation dispersion probes provided limited information regarding the sugar-backbone conformation. Here, we demonstrate measurements of C3′ and C4′ spin relaxation in the rotating frame (R1ρ) in uniformly ¹³C/¹⁵N labeled DNA as sensitive probes of the sugar-backbone conformation in DNA excited states. The chemical shifts, combined with structure-based predictions using an automated fragmentation quantum mechanics/molecular mechanics method, show that the dynamic ensemble of DNA duplexes containing m¹A–T Hoogsteen bps accurately model the excited state Hoogsteen conformation in two different sequence contexts. Formation of excited state A–T Hoogsteen bps is accompanied by changes in sugar-backbone conformation that allow the flipped syn adenine to form hydrogen-bonds with its partner thymine and this in turn results in overall kinking of the DNA toward the major groove. Results support the assignment of Hoogsteen bps as the excited state observed in canonical duplex DNA, provide an atomic view of DNA dynamics linked to formation of Hoogsteen bps, and lay the groundwork for a potentially general strategy for solving structures of nucleic acid excited states.     

The crystal structure of d(CCCCGGGG): a new A-form variant with an extended backbone conformation

The crystal structure of d(CCCCGGGG) has been determined at a resolution of 2.25 A. The oligomers crystallize as A-DNA duplexes occupying crystallographic two-fold axes. The backbone conformation is, in general, similar to that observed in previously reported crystal structures of A-DNA fragments, except for the central linkage, where it adopts an extended structure resulting from all trans conformation at the P-O5'-C5'-C4' bonds. This type of conformation facilitates interstrand stacking between the guanines at the C-G site. The local helix twist at this step is very small (25 degrees) compared to an overall average of 33.5 degrees. The unique structure of the C-G base-pair step, namely the extended backbone and the distinct stacking geometry, may be an important feature in the recognition mechanism between double-stranded DNA molecules and restriction endonucleases such as Msp I, which cuts the sequence CCGG very specifically with a rate unaffected by neighboring base pairs.     

Constructing optimal backbone segments for joining fixed DNA base pairs.

A method is presented to link a sequence of space-fixed base pairs by the sugar-phosphate segments of single nucleotides and to evaluate the effects in the backbone caused by this positioning of the bases. The entire computational unit comprises several nucleotides that are energy-minimized, subject to constraints imposed by the sugar-phosphate backbone segments being anchored to space-fixed base pairs. The minimization schemes are based on two stages, a conjugate gradient method followed by a Newton-Raphson algorithm. Because our purpose is to examine the response, or relaxation, of an artificially stressed backbone, it is essential to be able to obtain, as closely as possible, a lowest minimum energy conformation of the backbone segment in conformational space. For this purpose, an algorithm is developed that leads to the generation of an assembly of many local energy minima. From these sets of local minima, one conformation corresponding to the one with the lowest minimum is then selected and designated to represent the backbone segment at its minimum. The effective electrostatic potential of mean force is expressed in terms of adjustable parameters that incorporate solvent screening action in the Coulombic interactions between charged backbone atoms; these parameters are adjusted to obtain the best fit of the nearest-neighbor phosphorous atoms in an x-ray structure.     

Sequence dependence of DNA conformational flexibility

By using conformational free energy calculations, we have studied the sequence dependence of flexibility and its anisotropy along various conformational variables of DNA base pairs. The results show the AT base step to be very flexible along the twist coordinate. On the other hand, homonucleotide steps, GG(CC) and AA(TT), are among the most rigid sequences. For the roll motion that would correspond to a bend, the TA step is most flexible, while the GG(CC) step is least flexible. The flexibility of roll is quite anisotropic; the ratio of fluctuations toward the major and minor grooves is the largest for the GC step and the smallest for the AA(TT) and CG steps. Propeller twisting of base pairs is quite flexible, especially of A.T base pairs; propeller twist can reach 19 degrees by thermal fluctuation. We discuss the effect of electrostatic parameters, comparison with available experimental results, and biological relevance of these results.     

An ab initio molecular orbital study on the sequence-dependency of DNA conformation: an evaluation of intra- and inter-strand stacking interaction energy

The various nearest neighbor stacking interaction energies of stacked base pairs in the DNA double helix are calculated for both A- and B-type conformations using an ab initio molecular orbital method. It is demonstrated that the sequence-dependent conformational preference for A- or B-type results from the stacking interaction. In particular, the base sequence showing the highest preference for an A-type conformation is revealed as GC/GC, and the one with the next highest preference, AT/AT; for a B-type conformation, the respective sequences are CG/CG and CA/TG. The overall conformation of a DNA fragment is not determined by these particular sequences only but is influenced by all base pair steps. An intrinsically favorable conformation is predicted from the constituent stacking interaction.     

Refinement of the solution structure of the B DNA hexamer 5'd(C-G-T-A-C-G)2 on the basis of inter-proton distance data

A restrained least-squares refinement of the solution structure of the self-complementary B DNA hexamer 5'd(C-G-T-A-C-G)2 is presented. The structure is refined on the basis of 190 inter-proton distances determined by pre-steady-state nuclear Overhauser enhancement measurements. Two refinements were carried out starting from two initial B DNA structures differing by an overall root-mean-square (r.m.s.) difference of 0.32 A. In both cases, the final r.m.s. difference between the experimental and calculated inter-proton distances was 0.12 A compared to 0.61 A and 0.58 A for the two initial structures. The difference between the two refined structures is small, with an overall r.m.s. difference of 0.16 A, and represents the error in the refined co-ordinates. The refined structures have a B-type conformation with local structural variations in backbone and glycosidic bond torsion angles, and base-pair propellor twist, base roll, base tilt and local helical twist angles.     

Sequence-dependent conformation of an A-DNA double helix. The crystal structure of the octamer d(G-G-T-A-T-A-C-C)

The crystal structures of the synthetic self-complementary octamer d(G-G-T-A-T-A-C-C) and its 5-bromouracil-containing analogue have been refined to R values of 20% and 14% at resolutions of 1.8 and 2.25 A, respectively. The molecules adopt and A-DNA type double-helical conformation, which is minimally affected by crystal forces. A detailed analysis of the structure shows a considerable influence of the nucleotide sequence on the base-pair stacking patterns. In particular, the electrostatic stacking interactions between adjacent guanine and thymine bases produce symmetric bending of the double helix and a major-groove widening. The sugar-phosphate backbone appears to be only slightly affected by the base sequence. The local variations in the base-pair orientation are brought about by correlated adjustments in the backbone torsion angles and the glycosidic orientation. Sequence-dependent conformational variations of the type observed here may contribute to the specificity of certain protein-DNA interactions.     

Conformational characteristics and correlations in crystal structures of nucleic acid oligonucleotides: evidence for sub-states

Sugar phosphate backbone conformations are a structural element inextricably involved in a complete understanding of specific recognition nucleic acid ligand interactions, from early stage discrimination of the correct target to complexation per se, including any structural adaptation on binding. The collective results of high resolution DNA, RNA and protein/DNA crystal structures provide an opportunity for an improved and enhanced statistical analysis of standard and unusual sugar-phosphate backbone conformations together with corresponding dinucleotide sequence effects as a basis for further exploration of conformational effects on binding. In this study, we have analyzed the conformations of all relevant crystal structures in the nucleic acids data base, determined the frequency distribution of all possible epsilon, zeta, alpha, beta and gamma backbone angle arrangements within four nucleic acid categories (A-RNA and A-DNA, free and bound B-DNA) and explored the relationships between backbone angles, sugar puckers and selected helical parameters. The trends in the correlations are found to be similar regardless of the nucleic acid category. It is interesting that specific structural effects exhibited by the different unusual backbone sub-states are in some cases contravariant. Certain alpha/gamma changes are accompanied by C3' endo (north) sugars, small twist angles and positive values of base pair roll, and favor a displacement of nucleotide bases towards the minor groove compared to that of canonical B form structures. Unusual epsilon/zeta combinations occur with C2' (south) sugars, high twist angles, negative values of base pair roll, and base displacements towards the major groove. Furthermore, any unusual backbone correlates with a reduced dispersion of equilibrium structural parameters of the whole double helix, as evidenced by the reduced standard deviations of almost all conformational parameters. Finally, a strong sequence effect is displayed in the free oligomers, but reduced somewhat in the ligand bound forms. The most variable steps are GpA and CpA, and, to a lesser extent, their partners TpC and TpG. The results provide a basis for considering if the variable and non-variable steps within a biological active sequence precisely determine morphological structural features as the curvature direction, the groove depth, and the accessibility of base pair for non covalent associations.     

Solution structure of RNA duplexes containing alternating CG base pairs: NMR study of r(CGCGCG)2 and 2'-O-Me(CGCGCG)2 under low salt conditions.

Structures of r(CGCGCG)2 and 2'-O-Me(CGCGCG)2 have been determined by NMR spectroscopy under low salt conditions. All protons and phosphorus nuclei resonances have been assigned. Signals of H5'/5" have been assigned stereospecifically. All 3JH,H and 3JP,H coupling constants have been measured. The structures were determined and refined using an iterative relaxation matrix procedure (IRMA) and the restrained MD simulation. Both duplexes form half-turn, right-handed helices with several conformational features which deviate significantly from a canonical A-RNA structure. Duplexes are characterised as having C3'-endo sugar pucker, very low base-pair rise and high helical twist and inclination angles. Helices are overwound with <10 bp per turn. There is limited inter-strand guanine stacking for CG steps. Within CG steps of both duplexes, the planes of the inter-strand cytosines are not parallel while guanines are almost parallel. For the GC steps this pattern is reversed. The 2'-O-methyl groups are spatially close to the 5'-hydrogens of neighbouring residues from the 3'-side and are directed towards the minor groove of 2'-O-Me(CGCGCG)2 forming a hydrophobic layer. Solution structures of both duplexes are similar; the effect of 2'-O-methylation on the parent RNA structure is small. This suggests that intrinsic properties imposed by alternating CG base pairs govern the overall conformation of both duplexes.     

Sequence dependence of the B-A conformational transition of DNA

We have studied, by conformational analysis, the sequence dependence of DNA conformational transition between B- and A-forms. We have considered intramolecular interactions between base pairs, without backbone, to examine their role in the conformational transition between B- and A-forms, and found that base pairs themselves usually have intrinsic conformational preferences for the B- or A-form. Calculation of all ten possible base steps shows that the base combinations, CC (or GG), GC, AT, and TA, have tendencies to assume the A-conformation. Results show that it is particularly easy to slide along the long axis of the base pair for these steps, with AT and CC showing especially flat energies. These calculations show that a preference for the B- or A-conformation depends on the electrostatic energy parameters, in particular, on dielectric and shielding constants; the A-conformation is preferred for low dielectric constant or low shielding. Both the A- and B-conformations are mainly stabilized by electrostatic interactions between favorably juxtaposed atomic charges on base pairs; however, the B-conformation generally has more favorable van der Waals interactions than the A-form. These sequence-dependent conformational preference and environmental effects agree roughly with experimental observations, suggesting that the origin of the conformational polymorphism is attributable to the intrinsic conformational preference of base pairs.     

Sequence-dependent motions of DNA: a normal mode analysis at the base-pair level

Computer simulation of the dynamic structure of DNA can be carried out at various levels of resolution. Detailed high resolution information about the motions of DNA is typically collected for the atoms in a few turns of double helix. At low resolution, by contrast, the sequence-dependence features of DNA are usually neglected and molecules with thousands of base pairs are treated as ideal elastic rods. The present normal mode analysis of DNA in terms of six base-pair "step" parameters per chain residue addresses the dynamic structure of the double helix at intermediate resolution, i.e., the mesoscopic level of a few hundred base pairs. Sequence-dependent effects are incorporated into the calculations by taking advantage of "knowledge-based" harmonic energy functions deduced from the mean values and dispersion of the base-pair "step" parameters in high-resolution DNA crystal structures. Spatial arrangements sampled along the dominant low frequency modes have end-to-end distances comparable to those of exact polymer models which incorporate all possible chain configurations. The normal mode analysis accounts for the overall bending, i.e., persistence length, of the double helix and shows how known discrepancies in the measured twisting constants of long DNA molecules could originate in the deformability of neighboring base-pair steps. The calculations also reveal how the natural coupling of local conformational variables affects the global motions of DNA. Successful correspondence of the computed stretching modulus with experimental data requires that the DNA base pairs be inclined with respect to the direction of stretching, with chain extension effected by low energy transverse motions that preserve the strong van der Waals' attractions of neighboring base-pair planes. The calculations further show how one can "engineer" the macroscopic properties of DNA in terms of dimer deformability so that polymers which are intrinsically straight in the equilibrium state exhibit the mesoscopic bending anisotropy essential to DNA curvature and loop formation.     

Crystal structure at 1.5-A resolution of d(CGCICICG), an octanucleotide containing inosine, and its comparison with d(CGCG) and d(CGCGCG) structures.

The octadeoxyribonucleotide d(CGCICICG) has been crystallized in space group P(6)5(22) with unit cell dimensions of a = b = 31.0 A and c = 43.7 A, and X-ray diffraction data have been collected to 1.5-A resolution. Precession photographs and the self-Patterson function indicate that 12 base pairs of Z-conformation DNA stack along the c-axis, and the double helices pack in a hexagonal array similar to that seen in other crystals of Z-DNA. The structure has been solved by both Patterson deconvolution and molecular replacement methods and refined in space group P(6)5 to an R factor of 0.225 using 2503 unique reflections greater than 3.0 sigma (F). Comparison of the molecules within the hexagonal lattice with highly refined crystal structures of other Z-DNA reveals only minor conformational differences, most notably in the pucker of the deoxyribose of the purine residues. The DNA has multiple occupancy of C:I and C:G base pairs, and C:I base pairs adopt a conformation similar to that of C:G base pairs.     

The solution structure of a B-DNA undecamer comprising a portion of the specific target site for the cAMP receptor protein in the gal operon. Refinement on the basis of interproton distance data.

A restrained least squares refinement of the solution structure of the double-stranded DNA undecamer 5'd(AAGTGT-GACAT).5'd(ATGTCACACTT) comprising a portion of the specific target site of the cAMP receptor protein in the gal operon is presented. The structure is refined on the basis of both distance and planarity restraints, 2331 in all. The distance restraints comprise 150 interproton distances determined from pre-steady state nuclear Overhauser enhancement measurements and 2159 other interatomic distances derived from idealized geometry (i.e., distances between covalently bonded atoms, between atoms defining fixed bond angles, and between atoms defining hydrogen bonding in AT and GC base pairs). Two refinements were carried out and in both cases the final RMS difference between the experimental and calculated interproton distances was 0.2 A. The difference between the two refined structures is small (overall RMS difference of 0.23 A) and represents the error in the refined coordinates. Although the refined structures have an overall B-type conformation there are large variations in many of the local conformational parameters including backbone and glycosidic bond torsion angles, helical twist and propellor twist, base roll and base tilt angles.