Assignments of 31P NMR resonances in oligodeoxyribonucleotides: origin of sequence-specific variations in the deoxyribose phosphate backbone conformation and the 31P chemical shifts of double-helical nucleic acids

It is now possible to unambiguously assign all 31P resonances in the 31P NMR spectra of oligonucleotides by either two-dimensional NMR techniques or site-specific 17O labeling of the phosphoryl groups. Assignment of 31P signals in tetradecamer duplexes, (dTGTGAGCGCTCACA)2, (dTAT-GAGCGCTCATA)2, (dTCTGAGCGCTCAGA)2, and (dTGTGTGCGCACACA)2, and the dodecamer duplex d(CGTGAATTCGCG)2 containing one base-pair mismatch, combined with additional assignments in the literature, has allowed an analysis of the origin of the sequence-specific variation in 31P chemical shifts of DNA. The 31P chemical shifts of duplex B-DNA phosphates correlate reasonably well with some aspects of the Dickerson/Calladine sum function for variation in the helical twist of the oligonucleotides. Correlations between experimentally measured P-O and C-O torsional angles and results from molecular mechanics energy minimization calculations show that these results are consistent with the hypothesis that sequence-specific variations in 31P chemical shifts are attributable to sequence-specific changes in the deoxyribose phosphate backbone. The major structural variation responsible for these 31P shift perturbations appears to be P-O and C-O backbone torsional angles which respond to changes in the local helical structure. Furthermore, 31P chemical shifts and JH3'-P coupling constants both indicate that these backbone torsional angle variations are more permissive at the ends of the double helix than in the middle. Thus 31P NMR spectroscopy and molecular mechanics energy minimization calculations appear to be able to support sequence-specific structural variations along the backbone of the DNA in solution.     

Triple helical DNA in a duplex context and base pair opening

It is fundamental to explore in atomic detail the behavior of DNA triple helices as a means to understand the role they might play in vivo and to better engineer their use in genetic technologies, such as antigene therapy. To this aim we have performed atomistic simulations of a purine-rich antiparallel triple helix stretch of 10 base triplets flanked by canonical Watson–Crick double helices. At the same time we have explored the thermodynamic behavior of a flipping Watson–Crick base pair in the context of the triple and double helix. The third strand can be accommodated in a B-like duplex conformation. Upon binding, the double helix changes shape, and becomes more rigid. The triple-helical region increases its major groove width mainly by oversliding in the negative direction. The resulting conformations are somewhere between the A and B conformations with base pairs remaining almost perpendicular to the helical axis. The neighboring duplex regions maintain a B DNA conformation. Base pair opening in the duplex regions is more probable than in the triplex and binding of the Hoogsteen strand does not influence base pair breathing in the neighboring duplex region.     

The response of DNA length and twist to changes in ionic strength

We examine twist‐stretch coupling of unconstrained DNA using polyelectrolyte theory as applied to a line‐charge model along with published data on the ionic‐strength dependence of the twist angle. We conclude that twist‐stretch coupling is negative: environmental changes that stretch free DNA, unconstrained by externally applied pulling or twisting forces, are accompanied by unwinding of the double helix. We also analyze a helical model and conclude that the observed unwinding of the DNA helix when ionic strength is decreased is driven by radial swelling of the helix. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 223–226, 2015.     

Sequence-dependent twist-stretch coupling in DNA

Recent single-molecule micromanipulation experiments on DNA subject to small distortion revealed positive coupling between DNA stretching and twisting—for instance, DNA elongates when overtwisted. Here we propose a method to calculate the twist-stretch coupling constant specific to a DNA fragment of a given sequence. The method employs a sequence-dependent dinucleotide force field and is based on constrained minimization of the fragment's deformation energy. Using a force field inferred from atomistic molecular dynamics simulations, we obtain the twist-stretch coupling for random sequence to be 0.30 nm/turn, close to experimental values. An exhaustive calculation for all oligomers of nine basepairs yields values between 0.14 and 0.45 nm/turn, positively correlated with the contents of pyrimidine-purine steps in the sequence. Our method is simple to use and allows one to explore the hypothesis that some sequences may be optimized for twist-stretch coupling.     

Molecular structure of deoxycytidyl-3''-methylphosphonate (RP) 5''-deoxyguanidine, d[Cp(CH3)G]. A neutral dinucleotide with Watson-Crick base pairing and a right handed helical twist.

The crystal structure of d[Cp(CH3)G] has been determined as part of a project to study the mechanism of the B----Z transition in DNA. The asymmetric unit contains two dinucleotides and the equivalent of 7.5 water molecules, partially disordered over 12 definable positions. The two symmetry-independent dinucleotides form a duplex with Watson-Crick base-pairing and a right-handed helical sense. Comparison with previously determined structures of the B and A conformation showed that this duplex is closer to B than to A but significantly different from B. It corresponds to a stretched out helix with a 4 A rise per base pair and a helical twist of 32 degrees. This structure may serve as a model for the bending of DNA in certain situations. The configuration at the methyl phosphonate is RP, and a mechanism, based on this assignment, is presented for the B----Z transition in DNA.     

The mechanics of winding and unwinding helices in recombination: torsional stress associated with strand transfer promoted by RecA protein

Homologous recombination usually involves the production of heteroduplex DNA, DNA containing strands contributed from two different duplexes. RecA protein of E. coli can promote the formation of heteroduplex DNA in vitro by the exchange of DNA strands between two helical structures, duplex DNA and a helical recA nucleoprotein filament containing a single strand of DNA. Complete unwinding of the parental duplex and the rewinding of one strand with a new complement requires rotation of the helical structures about one another, or about their respective longitudinal axes. The observations described here demonstrate an association of torsional stress with strand exchange, and suggest that exchange is accomplished principally by concomitant rotation of duplex DNA and the recA nucleoprotein filament, each about its longitudinal axis.     

Modeling DNA Thermodynamics under Torsional Stress

Negatively twisted DNA is essential to many biological functions. Due to torsional stress, duplex DNA can have local, sequence-dependent structural defects. In this work, a thermodynamic model of DNA was built to qualitatively predict the local sequence-dependent mechanical instabilities under torsional stress. The results were compared to both simulation of a coarse-grained model and experiment results. By using the Kirkwood superposition approximation, we built an analytical model to represent the free energy difference ΔW of a hydrogen-bonded basepair between the B-form helical state and the basepair opened (or locally melted) state, within a given sequence under torsional stress. We showed that ΔW can be well approximated by two-body interactions with its nearest-sequence-neighbor basepairs plus a free energy correction due to long-range correlations. This model is capable of rapidly predicting the position and thermodynamics of local defects in a given sequence. The result qualitatively matches with an in vitro experiment for a long DNA sequence (>4000 basepairs). The 12 parameters used in this model can be further quantitatively refined when more experimental data are available.     

Bulged adenosine influence on the RNA duplex conformation in solution

The RNA single bulge motif is an unpaired residue within a strand of several complementary base pairs. To gain insight into structural changes induced by the presence of the adenosine bulge on RNA duplex, the solution structures of RNA duplex containing a single adenine bulge (5'-GCAGAAGAGCG-3'/5'-CGCUCUCUGC-3') and a reference duplex with all Watson-Crick base pairs (5'-GCAGAGAGCG-3'/5'-CGCUCUCUGC-3') have been determined by NMR spectroscopy. The reference duplex structure is a regular right-handed helix with all of the attributes of an A-type helix. In the bulged duplex, single adenine bulge stacks into the helix, and the bulge region forms a well-defined structure. Both structures were analyzed by the use of calculated helical parameters. Distortions induced by the accommodation of unpaired residue into the helical structure propagate over the entire structure and are manifested as the reduced base pairs inclination and x-displacement. Intrahelical position of bulged adenine A5 is stabilized by efficient stacking with 5'-neighboring residues G4.     

Blind Predictions of DNA and RNA Tweezers Experiments with Force and Torque

Single-molecule tweezers measurements of double-stranded nucleic acids (dsDNA and dsRNA) provide unprecedented opportunities to dissect how these fundamental molecules respond to forces and torques analogous to those applied by topoisomerases, viral capsids, and other biological partners. However, tweezers data are still most commonly interpreted post facto in the framework of simple analytical models. Testing falsifiable predictions of state-of-the-art nucleic acid models would be more illuminating but has not been performed. Here we describe a blind challenge in which numerical predictions of nucleic acid mechanical properties were compared to experimental data obtained recently for dsRNA under applied force and torque. The predictions were enabled by the HelixMC package, first presented in this paper. HelixMC advances crystallography-derived base-pair level models (BPLMs) to simulate kilobase-length dsDNAs and dsRNAs under external forces and torques, including their global linking numbers. These calculations recovered the experimental bending persistence length of dsRNA within the error of the simulations and accurately predicted that dsRNA''s “spring-like” conformation would give a two-fold decrease of stretch modulus relative to dsDNA. Further blind predictions of helix torsional properties, however, exposed inaccuracies in current BPLM theory, including three-fold discrepancies in torsional persistence length at the high force limit and the incorrect sign of dsRNA link-extension (twist-stretch) coupling. Beyond these experiments, HelixMC predicted that ‘nucleosome-excluding’ poly(A)/poly(T) is at least two-fold stiffer than random-sequence dsDNA in bending, stretching, and torsional behaviors; Z-DNA to be at least three-fold stiffer than random-sequence dsDNA, with a near-zero link-extension coupling; and non-negligible effects from base pair step correlations. We propose that experimentally testing these predictions should be powerful next steps for understanding the flexibility of dsDNA and dsRNA in sequence contexts and under mechanical stresses relevant to their biology.     

Modeling DNA Thermodynamics under Torsional Stress

Negatively twisted DNA is essential to many biological functions. Due to torsional stress, duplex DNA can have local, sequence-dependent structural defects. In this work, a thermodynamic model of DNA was built to qualitatively predict the local sequence-dependent mechanical instabilities under torsional stress. The results were compared to both simulation of a coarse-grained model and experiment results. By using the Kirkwood superposition approximation, we built an analytical model to represent the free energy difference ΔW of a hydrogen-bonded basepair between the B-form helical state and the basepair opened (or locally melted) state, within a given sequence under torsional stress. We showed that ΔW can be well approximated by two-body interactions with its nearest-sequence-neighbor basepairs plus a free energy correction due to long-range correlations. This model is capable of rapidly predicting the position and thermodynamics of local defects in a given sequence. The result qualitatively matches with an in vitro experiment for a long DNA sequence (>4000 basepairs). The 12 parameters used in this model can be further quantitatively refined when more experimental data are available.     

Conformational studies of (2''-5'') polynucleotides: theoretical computations of energy, base morphology, helical structure, and duplex formation.

A detailed theoretical analysis has been carried out to probe the conformational characteristics of (2'-5') polynucleotide chains. Semi-empirical energy calculations are used to estimate the preferred torsional combinations of the monomeric repeating unit. The resulting morphology of adjacent bases and the tendency to form regular single-stranded structures are determined by standard computational procedures. The torsional preferences are in agreement with available nmr measurements on model compounds. The tendencies to adopt base stacked and intercalative geometries are markedly depressed compared to those in (3'-5') chains. Very limited families of regular monomerically repeating single-stranded (2'-5') helices are found. Base stacking, however, can be enhanced (but helix formation is at the same time depressed) in mixed puckered chains. Constrained (2'-5') duplex structures have been constructed from a search of all intervening glycosyl and sugar conformations that form geometrically feasible phosphodiester linkages. Both A- and B-type base stacking are found to generate non-standard backbone torsions and mixed glycosyl/sugar combinations. The 2'- and 5'-residues are locked in totally different arrangements and are thereby prevented from generating long helical structures.     

A Renaissance Man: In Memoriam of Jon Widom (1955–2011)

Abstract

Assignment of the ¹H and ³¹P resonances of a decamer DNA duplex, d(CGCTTAAGCG)₂ was determined by two-dimensional COSY, NOESY and ¹H- ³¹P Pure Absorption phase Constant time (PAC) heteronuclear correlation spectroscopy. The solution structure of the decamer was calculated by an iterative hybrid relaxation matrix method combined with NOESY-distance restrained molecular dynamics. The distances from the 2D NOESY spectra were calculated from the relaxation rate matrix which were evaluated from a hybrid NOESY volume matrix comprising elements from the experiment and those calculated from an initial structure. The hybrid matrix-derived distances were then used in a restrained molecular dynamics procedure to obtain a new structure that better approximates the NOESY spectra. The resulting partially refined structure was then used to calculate an improved theoretical NOESY volume matrix which is once again merged with the experimental matrix until refinement is complete. J_H3′-P coupling constants for each of the phosphates of the decamer were obtained from ¹H-³¹P J-resolved selective proton flip 2D spectra. By using a modified Karplus relationship the C4′-C3′-03′-P torsional angles (?) were obtained. Comparison of the ³¹P chemical shifts and J_H3′-P coupling constants of this sequence has allowed a greater insight into the various factors responsible for ³¹P chemical shift variations in oligonucleotides. It also provides an important probe of the sequence-dependent structural variation of the deoxyribose phosphate backbone of DNA in solution. These correlations are consistent with the hypothesis that changes in local helical structure perturb the deoxyribose phosphate backbone. The variation of the ³¹P chemical shift, and the degree of this variation from one base step to the next is proposed as a potential probe of local helical conformation within the DNA double helix. The pattern of calculated ? and ζ torsional angles from the restrained molecular dynamics refinement agrees quite well with the measured J_H3′-P coupling constants. Thus, the local helical parameters determine the length of the phosphodiester backbone which in turn constrains the phosphate in various allowed conformations.     

Crystal structure of a 14 bp RNA duplex with non-symmetrical tandem GxU wobble base pairs

Adjacent GxU wobble base pairs are frequently found in rRNA. Atomic structures of small RNA motifs help to provide a better understanding of the effects of various tandem mismatches on duplex structure and stability, thereby providing better rules for RNA structure prediction and validation. The crystal structure of an RNA duplex containing the sequence r(GGUAUUGC-GGUACC)2 has been solved at 2.1 A resolution using experimental phases. Novel refinement strategies were needed for building the correct solvent model. At present, this is the only short RNA duplex structure containing 5'-U-U-3'/3'-G-G-5' non-symmetric tandem GxU wobble base pairs. In the 14mer duplex, the six central base pairs are all displaced away from the helix axis, yielding significant changes in local backbone conformation, helix parameters and charge distribution that may provide specific recognition sites for biologically relevant ligand binding. The greatest deviations from A-form helix occur where the guanine of a wobble base pair stacks over a purine from the opposite strand. In this vicinity, the intra-strand phosphate distances increase significantly, and the major groove width increases up to 3 A. Structural comparisons with other short duplexes containing symmetrical tandem GxU or GxT wobble base pairs show that nearest-neighbor sequence dependencies govern helical twist and the occurrence of cross-strand purine stacks.     

Binding of ethidium bromide to a DNA triple helix. Evidence for intercalation

The interaction of ethidium, a DNA intercalator, with the poly(dA).poly(dT) duplex and the poly (dA).2poly(dT) triplex has been investigated by a variety of spectrophotometric and hydrodynamic techniques. The fluorescence of ethidium is increased when either the duplex or triplex form is present. Binding constants, determined from absorbance measurements, indicate that binding to the triple helical form is substantially stronger than to the duplex, with a larger binding site size (2.8 base triplets compared to 2.4 base pairs). Furthermore, while binding to poly(dA).poly(dT) shows strong positive cooperativity, binding to the triplex is noncooperative. Thermal denaturation experiments demonstrate that ethidium stabilizes the triple helix. Binding to either form induces a weak circular dichroism band in the visible wavelength region, while in the region around 310 nm, there is a band that is strongly dependent on the degree of saturation of the duplex, and which is positive for the duplex but negative for the triplex. Both fluorescence energy transfer and quenching studies provide evidence of intercalation of ethidium in both duplex and triplex complexes. Binding of ethidium leads to an initial decrease in viscosity for both the duplex and triplex structures, followed by an increase, which is greater for the duplex. Taken together, these results strongly suggest that ethidium binds to the poly (dA).2poly(dT) triple helix via an intercalative mechanism.     

Torsional tension in the DNA double helix measured with trimethylpsoralen in living E. coli cells: analogous measurements in insect and human cells

The rate of covalent photobinding of trimethylpsoralen to DNA is greater when the DNA is wound with negative superhelical tension than when it is relaxed. In vitro the rate of photobinding is directly proportional to the negative superhelical density of the DNA. Thus measurement of the rate of photobinding provides an assay for probing in vivo unrestrained tension in the winding of the DNA double helix. This approach has been applied to measure torsional tension in DNA as it is packaged in living E. coli. Drosophila and HeLa cells. A method is described for measuring the rate of photobinding to intracellular DNA and rRNA, and for using the latter measurement as an internal control of the rate of me3-psoralen photobinding in vivo. This permits more accurate and reproducible measurement of changes in the DNA-psoralen photobinding reaction. The me3-psoralen probe interacts with intracellular bacterial DNA as expected for a purified DNA duplex wound with superhelical density sigma = -0.05 +/- 0.01. This superhelical tension is relaxed in cells when multiple single-strand breaks are introduced into the chromosomal DNA by gamma-irradiation. Similar relaxation occurs when cells are treated with the DNA gyrase inhibitor coumermycin. The results suggest that the DNA double helix is wound with torsional tension in vivo and that DNA supercoils which are equilibrated with this tension are not completely restrained in nucleosome-like structures. Torsional tension in the DNA of eucaryotic cells is not detectable in analogous measurements of the packaged DNA of HeLa and Drosophila cells. The simplest interpretation of this finding is that, within the limits of detection, all superhelical turns in the DNA are restrained in nucleosomes or nucleosome-like structures in these eucaryotic cells.     

The hydrogen-bonding structure in parallel-stranded duplex DNA is reverse Watson-Crick

Raman spectra of the parallel-stranded duplex formed from the deoxyoligonucleotides 5'-d-[(A)10TAATTTTAAATATTT]-3' (D1) and 5'-d[(T)10ATTAAAATTTATAAA]-3' (D2) in H2O and D2O have been acquired. The spectra of the parallel-stranded DNA are then compared to the spectra of the antiparallel double helix formed from the deoxyoligonucleotides D1 and 5'-d(AAATATTTAAAATTA-(T)10]-3' (D3). The Raman spectra of the antiparallel-stranded (aps) duplex are reminiscent of the spectra of poly[d(A)].poly[d(T)] and a B-form structure similar to that adopted by the homopolymer duplex is assigned to the antiparallel double helix. The spectra of the parallel-stranded (ps) and antiparallel-stranded duplexes differ significantly due to changes in helical organization, i.e., base pairing, base stacking, and backbone conformation. Large changes observed in the carbonyl stretching region (1600-1700 cm-1) implicate the involvement of the C(2) carbonyl of thymine in base pairing. The interaction of adenine with the C(2) carbonyl of thymine is consistent wtih formation of reverse Watson-Crick base pairing in parallel-stranded DNA. Phosphate-furanose vibrations similar to those observed for B-form DNA of heterogenous sequence and high A,T content are observed at 843 and 1092 cm-1 in the spectra of the parallel-stranded duplex. The 843-cm-1 band is due to the presence of a sizable population of furanose rings in the C2'-endo conformation. Significant changes observed in the regions from 1150 to 1250 cm-1 and from 1340 to 1400 cm-1 in the spectra of the parallel-stranded duplex are attributed to variations in backbone torsional and glycosidic angles and base stacking.     

Energetics of the hairpin to mismatched duplex transition of d(GCCGCAGC) on NaCl solution

We report Potential of Mean Force studies to describe the relative thermodynamic stabilities of d(GCCGCAGC) in a mismatched duplex and a hairpin monomer conformation in NaCl solution. The PMF calculations are combined with previous molecular mechanics and normal mode analysis in order to estimate the role of different components of the free energy in determining the relative stability of the duplex and hairpin structures. The high entropy associated with the loop region and the lack of minor groove phosphate-phosphate interactions in the hairpin compete against the gain in enthalpic contribution to the free energy due to base pairing in the mismatched duplex. The combined free energy calculations show that the hairpin is the most stable conformation at low salt and that a hairpin to duplex transition takes place at approximately 0.47 M NaCl. In addition, we studied the hairpin to partially stacked single helical conformation equilibrium at low salt. We found a small variation in transition temperature in salt concentration, delta Tm/delta log10(cs) approximately 2-3 degrees K/decade, in contrast to the duplex to hairpin or duplex to partially stacked single helix transition where the transition temperature exhibited marked dependence on salt concentration. This is in qualitative agreement with experimental data. Based on the Potential of Mean Force free energy calculation, the order of relative stability of the three-conformations studied varies with salt concentration. We observed the following orders of stability: stacked single helix greater than hairpin greater than duplex for cs less than 0.77 M NaCl; single helix greater than duplex greater than hairpin for 0.77 less than Cs less than 2.1 M; and duplex greater than hairpin greater than single strand for cs greater than 2.1 M. From the calculated PMF free energy curves in the NaCl concentration range, 0.012 less than cs less than 5.0 M, we can assign upper and lower bounds for the non-ionic differences in free energy between the duplex, hairpin, and stacked single helical states (at standard conditions: cs = 1.0 M, T = 25 degrees C, and 1 M oligomer concentration). We found that for delta G duplex single helix = G duplex - 2 x G single helix less than -7.38 Kcal/mol, the single helix is the least stable state. For the duplex-to-hairpin free energy difference in the range, -1.87 less than delta G duplex-hairpin less than 0.03 Kcal/mol, there will always be a salt-induced hairpin-to-duplex transition for 0.01 less than cs less than 1.6 M NaCl. If delta G duplex-hairpin less than -1.87, the duplex is always more stable than the hairpin; and for delta G duplex-hairpin greater than Kcal/mol, the hairpin state is always more stable than the duplex, for all salt concentrations.     

Interaction of the Escherichia coli HU protein with DNA. Evidence for formation of nucleosome-like structures with altered DNA helical pitch

Nuclease digestion studies of DNA bound to the histone-like protein HU show that cuts in each strand of the DNA double helix are made with a periodicity of 8.5 base-pairs. By contrast, similar digestions of DNA in eukaryotic nucleosomes show a repeat of 10.4 base-pairs. This and other results (including circular dichroism studies) are consistent with the proposal that the pitch of the DNA double helix in the HU complex is reduced from a repeat length of 10.5 to 8.5 base-pairs per helical turn. Simultaneously, the DNA in the HU-DNA complex containing two dimers of HU per 60 base-pairs has its linking number decreased by 1.0 turn per 290 base-pairs. From these changes it is calculated that HU imposes a DNA writhe of 1.0 per three to four monomers of HU. The results suggest a model in which DNA is coiled in left-handed toroidal supercoils on the HU complex, having a stoichiometry resembling that of the half-nucleosome of eukaryotic chromatin. An important distinction is that HU complexes can restrain the same number of DNA superhelical turns as eukaryotic nucleosomes, yet the DNA retains more negative torsional tension, just as is observed in prokaryotic chromosomes in vivo. Another distinction is that HU-DNA complexes are less stable, having a dissociation half-life of 0.6 min in 50 mM-NaCl. This last property may explain prior difficulties in detecting prokaryotic nucleosome-like structures.     

On the use of phasing experiments to measure helical repeat and bulge loop-associated twist in RNA.

In a phasing experiment, two bends are introduced into a long duplex RNA or DNA and the number of base pairs between them varied. When electrophoresed in a gel, the set of molecules may show a periodic variation in mobility that contains information about the twist associated with the bends and the intervening helix. We show how a set of three phasing experiments can be used to extract this information, and apply it to an RNA helix bend at the bulge sequence A2. The bulge introduces a negative (left-handed) twist of approximately 30 degrees; at low temperatures, it is mostly confined to the 5' side of the bulge. The apparent helical repeat of random sequence RNA measured in these experiments was 10.2 +/- 0.1 base pairs, an unexpectedly low value. It is likely that moderate curvative of the RNA helix axis (30-40 degrees over 80 bp) has affected the measurement.