Spatial configuration of ordered polynucleotide chains. V. Conformational energy estimates of helical structure

Semiempirical potential energy functional used previously to account successfully for the mean-square unperturbed dimensions and nmr coupling constants of randomly coiling polynucleotides are used, after modifications, to account for base stacking and interstrand hydrogen bonding, and to evaluate the conformational energies of single- and double-stranded polynucleotide helices. Attention is focused upon the variety of A-genus helices with local backbone conformations resembling the known double-helical structures of RNA. Distinct structural differences between single- and double-stranded helices are predicted from the energy calculations. A second point of interest is the apparent failure of two conformationally identical left-handed polynucleotide chains to form a left-handed duplex. The third major observation of the study is the wide morphological variety of theoreticaly allowed right-handed polynucleotide duplexes. In addition to the familiar double helix stabilized by horizontal base stacking and hydrogen bonding, an unusual vertical double helix is predicted to form between complementary bases fixed in the unusual but not energetically forbidden high anti glycosyl conformation. Experimental results bearing upon the theoretical predictions are discussed.     

Base-Stacking Interactions in Double-Helical DNA Structures: Experiment Versus Theory

Abstract

Atom-atom potential energy calculations have been undertaken for deriving stacking energies in double-helical structures. A comparison between the energy patterns of A- and B-type double-helical fragments determined by single-crystal X-ray diffraction methods versus idealized uniform models based on fiber diffraction data shows that the van der Waals stacking energy is largely sensitive to local changes in the relative orientation of adjacent base pairs. The sequence-dependent conformational variability observed in the high-resolution structures appears to be a consequence of the equipartitioning of the stacking energy along the double helix. The large energy variations expected for a uniform structure are dampened considerably in the observed structures by means of local changes in conformational features such as helix rotation and roll angles between base pairs.     

Structural Studies of DNA Fragments: The G·T Wobble Base Pair in A,B and Z DNA; The G·A Base Pair in B-DNA

Abstract

The crystal structures of five double helical DNA fragments containing non-Watson-Crick complementary base pairs are reviewed. They comprise four fragments containing G·T base pairs: two deoxyoctamers d(GGGGCTCC) and d(GGGGTCCC) which crystallise as A type helices; a deoxydodecamer d(CGCGAATTTGCG) which crystallises in the B-DNA conformation; and the deoxyhexamer d(TGCGCG), which crystallises as a Z-DNA helix. In all four duplexes the G and T bases form wobble base pairs, with bases in the major tautomer forms and hydrogen bonds linking N1 of G with 02 of T and 06 of G with N3 of T. The X-ray analyses establish that the G·T wobble base pair can be accommodated in the A, B or Z double helix with minimal distortion of the global conformation. There are, however, changes in base stacking in the neighbourhood of the mismatched bases. The fifth structure, d(CGCGAATTAGCG), contains the purine purine mismatch G·A where G is in the anti and A in the syn conformation. The results represent the first direct structure determinations of base pair mismatches in DNA fragments and are discussed in relation to the fidelity of replication and mismatch recognition.     

Molecular-Mechanical studies of the mismatched base analogs of d(CGCGAATTCGCG)2:d(CGTGAATTCGCG)2, d(CGAGAATTCGCG)2, d(CGCGAATTCACG)2, d(CGCGAATTCTCG)2, and d(CGCAGAATTCGCG)·d(CGCGAATTCGCG)

Molecular-mechanical simulations have been carried out on “mismatched base” analogs of the DNA double-helical structure d(CGCGAATTCGCG)₂, in which the base pairs CG at the 3 and 10 positions have been replaced by CA, AG, TC, and TG base pairs, as well as an insertion analog in which an extra adenine has been incorporated into one strand of the above structure between bases 3 and 4. The results of these simulations (calculated relative stabilities, structures, and nmr ring-current shifts) have been compared with calorimetric and nmr data. The calculated relative stabilities of the double-helical parent dodecamer and the various “wobble” base pairs qualitatively correlate with the experimental melting temperatures. The base-pairing structure for the GT wobble pair is in agreement with that previously determined from nmr experiments. For the GA base pair, the structure with both bases anti has a slightly more favorable energy from base pairing and stacking than a structure with non-Watson-Crick H-bonding with adenine syn, in agreement with nmr experiments. The CA wobble base is calculated to favor an adenine 6NH₂ …? cytosine N3 H-bond over cytosine 4NH₂ …? adenine N1, again, in agreement with nmr experiments. There is no definitive experimental data on the TC base pair, but the existence of (somewhat long and weak) H-bonds involving cytosine 4NH₂ …? thymine 4CO and cytosine N3 …? thymine HN3 seems reasonable. We find a structure in which the extra adenine base of the insertion analogs sits “inside” the double helix.     

Base-stacking interactions in double-helical DNA structures: experiment versus theory

Atom-atom potential energy calculations have been undertaken for deriving stacking energies in double-helical structures. A comparison between the energy patterns of A- and B-type double-helical fragments determined by single-crystal X-ray diffraction methods versus idealized uniform models based on fibre diffraction data shows that the van der Waals stacking energy is largely sensitive to local changes in the relative orientation of adjacent base pairs. The sequence-dependent conformational variability observed in the high-resolution structures appears to be a consequence of the equipartitioning of the stacking energy along the double helix. The large energy variations expected for a uniform structure are dampened considerably in the observed structures by means of local changes in conformational features such as helix rotation and roll angles between base pairs.     

Calculations of the circular dichroism of double-helical nucleic acids. I. Effects involving pi leads to pi transitions

The circular dichroism of double-helical nucleic acids was calculated as a function of geometry using the theory of Tinoco and Johnson. This theory does not include contributions of near ultraviolet transitions that are not π → π* in nature. The calculated circular dichroism shows a strong dependence on the distance of base pairs from the helix axis and the tilt of the base pair. Smaller dependences are predicted for the propeller-like twist of a base pair and for variation of the angular increment per base pair. Moderately good agreement between calculated and many experimentally observed spectra could be generated.     

The crystal structure of an ‘All Locked’ nucleic acid duplex

‘Locked nucleic acids’ (LNAs) are known to introduce enhanced bio- and thermostability into natural nucleic acids rendering them powerful tools for diagnostic and therapeutic applications. We present the 1.9 Å X-ray structure of an ‘all LNA’ duplex containing exclusively modified β-d-2′-O-4′C-methylene ribofuranose nucleotides. The helix illustrates a new type of nucleic acid geometry that contributes to the understanding of the enhanced thermostability of LNA duplexes. A notable decrease of several local and overall helical parameters like twist, roll and propeller twist influence the structure of the LNA helix and result in a widening of the major groove, a decrease in helical winding and an enlarged helical pitch. A detailed structural comparison to the previously solved RNA crystal structure with the corresponding base pair sequence underlines the differences in conformation. The surrounding water network of the RNA and the LNA helix shows a similar hydration pattern.     

Proton and phosphorus NMR studies of d-CpG(pCpG)n duplexes in solution. Helix-coil transition and complex formation with actinomycin-D.

The Watson–Crick imino and amino exchangeable protons, the nonexchangeable base and sugar protons, and the backbone phosphates for d-CpG(pCpG)_n, n = 1 and 2, have been monitored by high-resolution nmr spectroscopy in aqueous solution over the temperature range 0°–90°C. The temperature dependence of the chemical shifts of the tetramer and hexamer resonances is consistent with the formation of stable duplexes at low temperature in solution. Comparison of the spectral characteristics of the tetranucleotide with those of the hexanucleotide with temperature permits the differentiation and assignment of the cytosine proton resonances on base pairs located at the end of the helix from those in an interior position. There is fraying at the terminal base pairs in the tetranucleotide and hexanucleotide duplexes. The Watson–Crick ring imino protons exchange at a faster rate than the Watson–Crick side-chain amino protons, with exchange occurring by transient opening of the double helix. The structure of the d-CpG(pCpG)_n double helices has been probed by proton relaxation time measurements, sugar proton coupling constants, and the proton chemical shift changes associated with the helix–coil transition. The experimental data support a structural model in solution, which incorporates an anti conformation about the glycosyl bonds, C(3) exo sugar ring pucker, and base overlap geometries similar to the B-DNA helix. Rotational correlation times of 1.7 and 0.9 × 10^?9 sec have been computed for the hexanucleotide and tetranucleotide duplexes in 0.1 M salt, D₂O, pH 6.25 at 27°C. The well-resolved ³¹P resonances for the internucleotide phosphates of the tetramer and hexamer sequences at superconducting fields shift upfield by 0.2–0.5 ppm on helix formation. These shifts reflect a conformational change about the ω,ω′ phosphodiester bonds from gauche-gauche in the duplex structure to a distribution of gauche-trans states in the coil structure. Significant differences are observed in the transition width and midpoint of the chemical shift versus temperature profiles plotted in differentiated form for the various base and sugar proton and internucleotide phosphorous resonances monitoring the d-CpG(pCpG)_n helix–coil transition. The twofold symmetry of the d-CpGpCpG duplex is removed on complex formation with the antibiotic actinomycin-D. Two phosphorous resonances are shifted downfield by ～2.6 ppm and ～1.6 ppm on formation of the 1:2 Act-D:d-CpGpCpG complex in solution. Model studies on binding of the antibiotic to dinucleotides of varying sequence indicate that intercalation of the actinomycin-D occurs at the GpC site in the d-CpGpCpG duplex and that the magnitude of the downfield shifts reflects strain at the O-P-O backbone angles and hydrogen bonding between the phenoxazone and the phosphate oxygens. Actinomycin-D is known to bind to nucleic acids that exhibit a B-DNA conformation; this suggests that the d-CpG(pCpG)_n duplexes exhibit a B-DNA conformation in solution.     

Helical coherence of DNA in crystals and solution

The twist, rise, slide, shift, tilt and roll between adjoining base pairs in DNA depend on the identity of the bases. The resulting dependence of the double helix conformation on the nucleotide sequence is important for DNA recognition by proteins, packaging and maintenance of genetic material, and other interactions involving DNA. This dependence, however, is obscured by poorly understood variations in the stacking geometry of the same adjoining base pairs within different sequence contexts. In this article, we approach the problem of sequence-dependent DNA conformation by statistical analysis of X-ray and NMR structures of DNA oligomers. We evaluate the corresponding helical coherence length—a cumulative parameter quantifying sequence-dependent deviations from the ideal double helix geometry. We find, e.g. that the solution structure of synthetic oligomers is characterized by 100–200 Å coherence length, which is similar to ~150 Å coherence length of natural, salmon-sperm DNA. Packing of oligomers in crystals dramatically alters their helical coherence. The coherence length increases to 800–1200 Å, consistent with its theoretically predicted role in interactions between DNA at close separations.     

The influence of complementary base pair interactions on secondary structure formation and conformational transitions in polynucleotides]

The calculations have been carried out of interaction energy between complementary base pairs of nucleic acids in the function of conformational parametres of double helix (Arnott's parameters) by the method of atom-atom potential functions. Interaction energy as a function of conformational parametres is valley-like and varies little along the bottom of the valley. The regions of interaction energy minima are compared with experimentally determined conformational parametres of nucleic acid double helices. On the basis of calculation results the pathways of conformational transitions between different forms of double-helical polynucleotides are discussed.     

Index to Authors,Volume 3

Abstract

Results of calculations using various empirical potentials suggest that base pair buckling, which commonly occurs in DNA crystal structures, is sufficient to eliminate the steric clash at CpG steps in B-DNA, originating from the base pair propeller twisting. The buckling is formed by an inclination of cytosines while deviations of guanines from a plane perpendicular to the double helix axis are unfavorable. The buckling is accompanied by an increased vertical separation of the base pair centers but the buckled arrangement of base pairs is at least as stable as when the vertical separation is normal and buckle zero. In addition, room is created by the increased vertical separation for the bases to propeller twist as is observed in DNA crystal structures. Further stabilization of base stacking is introduced into the buckled base pair arrangement by roll opening the base pairs into the double helix minor groove. The roll may lead to the double helix bending and liberation of guanines from the strictly perpendicular orientation to the double helix axis. The liberated guanines further contribute to the base pair buckling and stacking improvement. This work also suggests a characteristic very stable DNA structure promoted by nucleotide sequences in which runs of purines follow runs of pyrimidine bases.     

Interactions of molecules with nucleic acids. I. An algorithm to generate nucleic acid structures with an application to the B-DNA structure and a counterclockwise helix.

An algorithm is developed that enables the routine determination of backbone conformations of nucleic acids. All atomic positions including hydrogen are specified in accord with experimental bond lengths and angles but with theoretically determined conformational angles. For two Watson-Crick base pairs at a separation of 3.38 Å, and perpendicular to a common helical axis, minimum energy configurations are found for all 10 combinations at helical angles of α ～ 36°–38°, corresponding to the B-DNA structure with C(2′)-endo sugar puckers. Backbone configurations exist only within the range 35.5° ? α ? 42°, which suggests the origin of the 10-fold helix. Calculated stacking energies for the B-DNA structure increases for each of the clustered groups of base pairs: G·C with G·C, G·C with A·T, and A·T with A·T, and they are in approximate agreement with experimental observations. The counter-clockwise helix is examined, and physically meaningful structures are found only when the helical axes of successive base pairs are disjointed.     

Calculation of the circular dichroism of double-helical nucleic acids. 3. Dependence of the results on the choice of wave functions for purines and pyrimidines

The circular dichroism of double-helical nucleic acids was calculated using three different sets of wave functions for the purine and pyrimidine chromophores. Different wave functions give qualitatively the same types of spectra. Quantitatively, the differences are very substantial. However, the dependence of calculated circular dichroism spectra on nucleic acid geometry and base composition is very similar for the three different sets of wave functions. Each set of wave functions could be used to reproduce some, but not all experimentally observed spectra. The geometries required to reproduce experimentally observed spectra consistently require double-helical geometries similar to the B or C forms of DNA.     

Conformational and structural constraints in double-helical polynucleotides

Conformational analysis of antiparallel double-helical polynucleotides with Watson-Crick base pairing was reduced to a four-dimensional problem using original mathematical methods. In the four-dimensional conformational space the family of structures, characterized by the base-pair stacking with the most stable conformations in water solution as well as in the solid state, was localized. For the C′₂-endo sugar pucker, both right-handed and left-handed structures were found; right-handed structures only, however, seem to be allowed for the C′₃-endo pucker, the only possible one for ribonucleotides with base stacking.     

Conformational features of benzo‐homologated yDNA duplexes by molecular dynamics simulation

yDNA is a base‐modified nucleic acid duplex containing size‐expanded nucleobases. Base‐modified nucleic acids could expand the genetic alphabet and thereby enhance the functional potential of DNA. Unrestrained 100 ns MD simulations were performed in explicit solvent on the yDNA NMR sequence [5′(yA T yA yA T yA T T yA T)₂] and two modeled yDNA duplexes, [5′(yC yC G yC yC G G yC G G)₂] and [(yT5′ G yT A yC yG C yA yG T3′)?(yA5′ C T C yG C G yT A yC A3′)]. The force field parameters for the yDNA bases were derived in consistent with the well‐established AMBER force field. Our results show that DNA backbone can withstand the stretched size of the bases retaining the Watson‐Crick base pairing in the duplexes. The duplexes retained their double helical structure throughout the simulations accommodating the strain due to expanded bases in the backbone torsion angles, sugar pucker and helical parameters. The effect of the benzo‐expansion is clearly reflected in the extended C1′‐C1′ distances and enlarged groove widths. The size expanded base modification leads to reduction in base pair twist resulting in larger overlapping area between the stacked bases, enhancing inter and intra strand stacking interactions in yDNA in comparison with BDNA. This geometry could favour enhanced interactions with the groove binders and DNA binding proteins., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 55–64, 2016     

Adenine-guanine base pairing ribosomal RNA.

Analyses of secondary structures proposed for ribosomal RNA's show that, of the different kinds of base pairs directly adjoining the ends of postulated double-helical regions, only A-G with A at the 5' end significantly exceeds the number expected for a random base distribution. An A(syn)-G(trans) hydrogen-bonded basepair is proposed. This could fit at the end of an undistorted double helix, but would prevent further base stacking, thus favoring a break in the double helix to produce a non-linear tertiary structure.     

Mechanical properties of DNA-like polymers

The molecular structure of the DNA double helix has been known for 60 years, but we remain surprisingly ignorant of the balance of forces that determine its mechanical properties. The DNA double helix is among the stiffest of all biopolymers, but neither theory nor experiment has provided a coherent understanding of the relative roles of attractive base stacking forces and repulsive electrostatic forces creating this stiffness. To gain insight, we have created a family of double-helical DNA-like polymers where one of the four normal bases is replaced with various cationic, anionic or neutral analogs. We apply DNA ligase-catalyzed cyclization kinetics experiments to measure the bending and twisting flexibilities of these polymers under low salt conditions. Interestingly, we show that these modifications alter DNA bending stiffness by only 20%, but have much stronger (5-fold) effects on twist flexibility. We suggest that rather than modifying DNA stiffness through a mechanism easily interpretable as electrostatic, the more dominant effect of neutral and charged base modifications is their ability to drive transitions to helical conformations different from canonical B-form DNA.     

Conformation of 2-aminofluorene-modified DNA

Minimized potential-energy calculations were performed to determine the conformation of the 2-aminofluorene (AF) adduct to dCpdG at guanine C-8. The AF adduct has many low-energy conformers in both the anti and the syn domains of the guanine. This is in contrast with the acetylated adduct, (AAF), which greatly prefers the syn domain. Two types of low-energy guanine anti-conformations were obtained: (1) conformers that preserve guanine–cytidine stacking and (2) conformers with fluorene–cytidine stacking. Of special importance are conformers with ω′,ω,ψ = g^?, g^?, g⁺, characteristic of normal A- or B-helices, which are found in both groups. No conformers of this type were obtained for the acetylated AAF adduct. The guanine–cytidine stacked from with this conformation can be incorporated in the B-helix without any distortion, with the carcinogen situated at the helix exterior. The fluorene in this model can slide into the helix to yield a fluorene–cytidine stacked minimum-energy conformation. This requires no denaturation, although one base pair is unstacked and the helix axis is bent. Low-energy syn-conformations, similar to those obtained for the AAF adduct, were also computed. These were either guanine-cytidine stacked or fluorene–cytidine stacked. The syn froms are less likely to be observed in larger DNA polymers of the adduct, since they cause more distortion than the anti-conformations. However, they might well be observed in crystals of small subunits, and they should contribute significantly to the population in solution.     

Base pair buckling can eliminate the interstrand purine clash at the CpG steps in B-DNA caused by the base pair propeller twisting

Results of calculations using various empirical potentials suggest that base pair buckling, which commonly occurs in DNA crystal structures, is sufficient to eliminate the steric clash at CpG steps in B-DNA, originating from the base pair propeller twisting. The buckling is formed by an inclination of cytosines while deviations of guanines from a plane perpendicular to the double helix axis are unfavorable. The buckling is accompanied by an increased vertical separation of the base pair centers but the buckled arrangement of base pairs is at least as stable as when the vertical separation is normal and buckle zero. In addition, room is created by the increased vertical separation for the bases to propeller twist as is observed in DNA crystal structures. Further stabilization of base stacking is introduced into the buckled base pair arrangement by roll opening the base pairs into the double helix minor groove. The roll may lead to the double helix bending and liberation of guanines from the strictly perpendicular orientation to the double helix axis. The liberated guanines further contribute to the base pair buckling and stacking improvement. This work also suggests a characteristic very stable DNA structure promoted by nucleotide sequences in which runs of purines follow runs of pyrimidine bases.