Spatial configuration of ordered polynucleotide chains: a novel double helix

Biophysical Reviews' "Meet the Editors Series"-a profile of Wilma K. Olson

As one of the five Executive Editors of Biophysical Reviews I have been asked to provide this short biographical sketch for the readers of the journal. I have been a member of the Editorial Board since the inception of the journal in 2008 and an Executive Editor since 2014. I hold a B.S. degree in Chemistry from the University of Delaware and a Ph.D. in Physical Chemistry from Stanford University. Except for a year as a Damon Runyon Postdoctoral Fellow at Columbia University, I have spent my entire professional career at Rutgers, the State University of New Jersey, where I am currently the Mary I. Bunting Professor of Chemistry and Chemical Biology. I served for many years as Founding Director of the Rutgers University Center for Molecular Biophysics and Biophysical Chemistry, and have trained undergraduate, graduate, and postdoctoral students from a variety of academic disciplines.

Atomistic simulations reveal bubbles, kinks and wrinkles in supercoiled DNA

Although DNA is frequently bent and supercoiled in the cell, much of the available information on DNA structure at the atomistic level is restricted to short linear sequences. We report atomistic molecular dynamics (MD) simulations of a series of DNA minicircles containing between 65 and 110 bp which we compare with a recent biochemical study of structural distortions in these tight DNA loops. We have observed a wealth of non-canonical DNA structures such as kinks, denaturation bubbles and wrinkled conformations that form in response to bending and torsional stress. The simulations show that bending alone is sufficient to induce the formation of kinks in circles containing only 65 bp, but we did not observe any defects in simulations of larger torsionally relaxed circles containing 110 bp over the same MD timescales. We also observed that under-winding in minicircles ranging in size from 65 to 110 bp leads to the formation of single stranded bubbles and wrinkles. These calculations are used to assess the ability of atomistic MD simulations to determine the structure of bent and supercoiled DNA.

Backbone-base inclination as a fundamental determinant of nucleic acid self- and cross-pairing

The crystal structure of the duplex formed by oligo(2′,3′-dideoxy-β-d-glucopyranosyl)nucleotides (homo-DNA) revealed strongly inclined backbone and base-pair axes [Egli,M., Pallan,P.S., Pattanayek,R., Wilds,C.J., Lubini,P., Minasov,G., Dobler,M., Leumann,C.J. and Eschenmoser,A. (2006) Crystal structure of homo-DNA and nature's choice of pentose over hexose in the genetic system. J. Am. Chem. Soc., 128, 10847–10856]. This inclination is easily perceived because homo-DNA exhibits only a modest helical twist. Conversely, the tight coiling of strands conceals that the backbone-base inclinations for A- (DNA and RNA) and B-form (DNA) duplexes differ considerably. We have defined a parameter η_B that corresponds to the local inclination between sugar-phosphate backbone and base plane in nucleic acid strands. Here, we show its biological significance as a predictive measure for the relative strand polarities (antiparallel, aps, or parallel, ps) in duplexes of DNA, RNA and artificial nucleic acid pairing systems. The potential of formation of ps duplexes between complementary 16-mers with eight A and U(T) residues each was investigated with DNA, RNA, 2′-O-methylated RNA, homo-DNA and p-RNA, the ribopyranosyl isomer of RNA. The thermodynamic stabilities of the corresponding aps duplexes were also measured. As shown previously, DNA is capable of forming both ps and aps duplexes. However, all other tested systems are unable to form stable ps duplexes with reverse Watson–Crick (rWC) base pairs. This observation illustrates the handicap encountered by nucleic acid systems with inclinations η_B that differ significantly from 0° to form a ps rWC paired duplex. Accordingly, RNA with a backbone-base inclination of −30°, pairs strictly in an aps fashion. On the other hand, the more or less perpendicular orientation of backbone and bases in DNA allows it to adopt a ps rWC paired duplex. In addition to providing a rationalization of relative strand polarity with nucleic acids, the backbone-base inclination parameter is also a determinant of cross-pairing. Thus, systems with strongly deviating η_B angles will not pair with each other. Nucleic acid pairing systems with significant backbone-base inclinations can also be expected to display different stabilities depending on which terminus carries unpaired nucleotides. The negative inclination of RNA is consistent with the higher stability of duplexes with 3′- compared to those with 5′-dangling ends.

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.

DNA stretching and compression: large-scale simulations of double helical structures

Computer-simulated elongation and compression of A - and B -DNA structures beyond the range of thermal fluctuations provide new insights into high energy "activated" forms of DNA implicated in biochemical processes, such as recombination and transcription. All-atom potential energy studies of regular poly(dG).poly(dC) and poly(dA).poly(dT) double helices, stretched from compressed states of 2.0 A per base-pair step to highly extended forms of 7.0 A per residue, uncover four different hyperfamilies of right-handed structures that differ in mutual base-pair orientation and sugar-phosphate backbone conformation. The optimized structures embrace all currently known right-handed forms of double-helical DNA identified in single crystals as well as non-canonical forms, such as the original "Watson-Crick" duplex with trans conformations about the P-O5' and C5'-C4' backbone bonds. The lowest energy minima correspond to canonical A and B -form duplexes. The calculations further reveal a number of unusual helical conformations that are energetically disfavored under equilibrium conditions but become favored when DNA is highly stretched or compressed. The variation of potential energy versus stretching provides a detailed picture of dramatic conformational changes that accompany the transitions between various families of double-helical forms. In particular, the interchanges between extended canonical and non-canonical states are reminiscent of the cooperative transitions identified by direct stretching experiments. The large-scale, concerted changes in base-pair inclination, brought about by changes in backbone and glycosyl torsion angles, could easily give rise to the observed sharp increase in force required to stretch single DNA molecules more than 1.6-1.65 times their canonical extension. Our extended duplexes also help to tie together a number of previously known structural features of the RecA-DNA complex and offer a self-consistent stereochemical model for the single-stranded/duplex DNA recognition brought in register by recombination proteins. The compression of model duplexes, by contrast, yields non-canonical structures resembling the deformed steps in crystal complexes of DNA with the TATA-box binding protein (TBP). The crystalline TBP-bound DNA steps follow the calculated compression-elongation pattern of an unusual "vertical" duplex with base planes highly inclined with respect to the helical axis, exposed into the minor groove, and accordingly accessible for recognition.Significantly, the double helix can be stretched by a factor of two and compressed roughly in half before its computed internal energy rises sharply. The energy profiles show that DNA extension-compression is related not only to the variation of base-pair Rise but also to concerted changes of Twist, Roll, and Slide. We suggest that the high energy "activated" forms calculated here are critical for DNA processing, e.g. nucleo-protein recognition, DNA/RNA synthesis, and strand exchange.

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.

Local DNA stretching mimics the distortion caused by the TATA box-binding protein

X-ray structures of the TATA box-binding protein complexed with its DNA target show that the nucleic acid is severely bent away from the protein and also strongly unwound. We have used molecular mechanics and energy mapping to understand how such an unusual conformation can be induced. The results show that simple deformation pathways involving local stretching or unwinding of DNA reproduce many features of the experimental structure. Notably, kinked junctions with the flanking B-DNA regions occur without the need for any specific local

Metropolis Monte Carlo calculations of DNA structure using internal coordinates and NMR distance restraints: an alternative method for generating a high-resolution solution structure

A new method, a restrained Monte Carlo (rMC) calculation, is demonstrated for generating high-resolution structures of DNA oligonucleotides in solution from interproton distance restraints and bounds derived from complete relaxation matrix analysis of two-dimensional nuclear Overhauser effect (NOE) spectral peak intensities. As in the case of restrained molecular dynamics (rMD) refinement of structures, the experimental distance restraints and bounds are incorporated as a pseudo-energy term (or penalty function) into the mathematical expression for the molecular energy. However, the use of generalized helical parameters, rather than Cartesian coordinates, to define DNA conformation increases efficiency by decreasing by an order of magnitude the number of parameters needed to describe a conformation and by simplifying the potential energy profile. The Metropolis Monte Carlo method is employed to simulate an annealing process. The rMC method was applied to experimental 2D NOE data from the octamer duplex d(GTATAATG).d(CATTATAC). Using starting structures from different locations in conformational space (e.g. A-DNA and B-DNA), the rMC calculations readily converged, with a root-mean-square deviation (RMSD) of < 0.3 A between structures generated using different protocols and starting structures. Theoretical 2D NOE peak intensities were calculated for the rMC-generated structures using the complete relaxation matrix program CORMA, enabling a comparison with experimental intensities via residual indices. Simulation of the vicinal proton coupling constants was carried out for the structures generated, enabling a comparison with the experimental deoxyribose ring coupling constants, which were not utilized in the structure determination in the case of the rMC simulations. Agreement with experimental 2D NOE and scalar coupling data was good in all cases. The rMC structures are quite similar to that refined by a traditional restrained MD approach (RMSD < 0.5 A) despite the different force fields used and despite the fact that MD refinement was conducted with additional restraints imposed on the endocyclic torsion angles of deoxyriboses. The computational time required for the rMC and rMD calculations is about the same. A comparison of structural parameters is made and some limitations of both methods are discussed with regard to the average nature of the experimental restraints used in the refinement.

Paranemic structures of DNA and their role in DNA unwinding

A DNA structure is defined as paranemic if the participating strands can be separated without mutual rotation of the opposite strands. The experimental methods employed to detect paranemic, unwound, DNA regions is described, including probing by single-strand specific nucleases (SNN), conformation-specific chemical probes, topoisomer analysis, NMR, and other physical methods. The available evidence for the following paranemic structures is surveyed: single-stranded DNA, slippage structures, cruciforms, alternating B-Z regions, triplexes (H-DNA), paranemic duplexes and RNA, protein-stabilized paranemic DNA. The problem of DNA unwinding during gene copying processes is analyzed; the possibility that extended paranemic DNA regions are transiently formed during replication, transcription, and recombination is considered, and the evidence supporting the participation of paranemic DNA forms in genes committed to or undergoing copying processes is summarized.

Structural and conformational studies on deoxyguanosyl-3',5'-deoxyadenosine monophosphate and its ethyl phosphotriester analogs--left-handed dimers

The mode of base-base stacking, the handedness and the sugar(dGpA)phosphate backbone conformation of deoxyguanosyl 3'-5' deoxyadenosine and its diastereomeric ethyl phosphotriester analogs were studied by 1H NMR, UV and CD spectroscopy. The results indicate the three dimers are left-handed, while the sugar phosphate backbone is comprised predominantly of C2-endo,gg(C4-C5) and g'g (C5-O) conformers. The two bases are extensively stacked and interact about 90 degrees along the dyad axes. The extent of base overlap in dGpA is slightly greater than in either ethyl phosphotriester analog. The absolute configurations of the two ethyl phosphotriester diastereoisomers of dGpA can be assigned by one-dimensional and two-dimensional 1H NMR nuclear Overhauser enhancement experiments.

Conformational studies of trinucleoside bisphosphates. 2. Potential energy calculations

Semi-empirical potential energy calculations were carried out for the conformations of ApApA, a trinucleoside bisphosphate (trimer), as proposed in the preceding paper in this journal. Energy minimals were obtained for these conformers, i. e. conformers I-I, I-II, III-I' and I'-I' (the first Roman numeral represents the conformation of the ApAp- moiety and the second one the -pApA moiety of ApApA), of the nearest-neighbor interactions, and conformers B1, B2 and B3 of bulged structures. The results show that the detailed conformations of the component dinuclotide (dimer) moieties in I-I, I-II, III-I' and I'-I' are close (+/- 10 degrees) to the corresponding stable conformations of dinucleotides in solution proposed by Lee and coworkers. The bulged conformations B1, B2 and B3 have open forms for each dinucleotide moiety. Their preferred local conformations are anti, gamma +, beta +, epsilon + (except for Ap- in B3, epsilon-), and 3' endo (except for Ap- in B3, 2' endo) for the glycosidic bond, C4'-C5', C5'-O5', C3'-O3' and the ribose moiety, respectively. The conformation of the phosphodiester linkage in -pApA moiety is zeta t2, alpha t3 for all the three bulged conformers, while that in ApAp- is zeta- 1, alpha t2 for B1, zeta+ 1 alpha- 2 for B2, and zeta+ 1, alpha+ 2 for B3. The calculated composite chemical shifts of the base protons from ring-current effects agree with the corresponding observed trimerization shifts. Conformers I-I and I'-I' are helix-like structures and correspond to the right-handed and the 'vertical-stacked' helices, respectively. The structure of I-II is similar to that observed in the ApApA single crystals. The bulged conformations may serve as models for frameshift mutations.

Conformational studies of 13 trinucleoside bisphosphates by 360-MHz 1H-NMR spectroscopy. 1. Ribose protons

The ribose protons of 13 trinucleoside bisphosphates (trimers) were studied, using 360-MHz proton nuclear magnetic resonance spectroscopy. Complete assignments and analyses of the NMR signals of these protons were carried out by the methods of homonuclear decoupling and computer line-shape simulations. It was shown that the trinucleotides preferred the anti, 3' endo, gamma +, beta t and epsilon t/epsilon- conformations for the glycosidic torsions, the ribose rings, the C4'-C5' bonds, the C5'-O5' bonds, and the C3'-O3' bonds, respectively. It was also found that the trimers, especially those which had noticeable population of 'bulged' structures, did not necessarily have a higher population of these preferred local conformations than their component dimers. The overall conformations of the trinucleotides are classified into two categories. The conformations in the first category involve the nearest-neighbor interactions. Each dinucleotide moiety can assume one of the four stable conformations (I, I', II and III) or the open forms of dinucleoside monophosphates. However, due to steric hindrance, there are only four cases in which both dinucleotide moieties can assume one of the four stable conformations at the same time. These four combinations of conformations are I-I, I'-I', I-II and III-I', where the first Roman numeral represents the conformation of the NpN'p-moiety and the second one, that of the -pN'pN'' moiety of the trimers. Among them, I-I and I'-I' are helical structures, capable of forming a double helix. The second category contains conformations with bulged structures which have the two dinucleotide moieties in open forms (i.e. no nearest-neighbor interactions) and the bases of the two terminal residues stacking on each other while the middle residue is bulged out. These bulged conformations may serve as structural models for frame-shift mutations.

Determination of handedness of DNA double helices from NOE difference spectra: the structure of poly(dG-dC).poly(dG-dC) in low salt

Combined use of shielding constant computations, measurements of chemical shifts and NOE studies reveal that poly(dG-dC).(poly)dG-dC) in low salt solutions exist as a right-handed B-DNA double helix described by Gupta, Dhingra, Sarma, Sarma, Rajagopalan and Sasisekharan, J. Biomole. Str. Dyn. 1. 395, 1983. We present a simple and direct method to determine the handedness of DNA double helices from NOE difference spectra. This method takes advantage of the NOE between base protons and the H2'H2" sugar protons; and in the difference NOE spectra in the H2'H2" region the signatures of the right and left-handed helices become imprinted.

Alternate g-g- conformation for dinucleoside phosphates in solution

An alternative g-g- conformation (conformer I') for dinucleosides in solution has been deduced, based on potential energy calculations and nuclear magnetic resonance spectroscopy. This conformation is characterized by larger glycosidic torsional angles (chi = 94-111 degrees) than those of conformer I (chi = 8-35 degrees), although the other torsional angles are similar. There are thus four stable conformers (I, I', II and III) for dinucleosides in equilibrium with the open forms. The structure of conformer I' supports that of the 'vertical' double helix constructed by Olson (W.K. Olson, Proc. Natl. Acad. Sci. U.S.A. 74, (1977) 1775). Our data may suggest the possibility of interconversion between the vertical double helix and the regular double helix of A-form DNA, RNA or A'-form RNA.

Interferon induction: a conformational hypothesis

The ability of polynucleotides or polynucleotide duplexes such as poly(I).poly(C) to induce interferon production is proposed to depend on the existence of certain stable glycosidic orientations. It appears that a slight increase in instability of 1--3 kcal/mole (1 cal = 4.184 J) in the conformational regions near 20 degrees, 80 degrees, and 160 degrees leads to a loss of potency with respect to interferon induction. Thus, it is proposed that, for a polynucleotide to exist in the overall conformation necessary for interferon induction, stability of glycosidic orientations near 20 degrees, 80 degrees, and 160 degrees may be necessary to confer flexibility and activity on polynucleotide structures. This proposed conformational triad of stable conformational regions essential to interferon induction is based on the results of conformational energy calculations of the glycoside rotational profiles of adenosine, 7-deazaadenosine, inosine, and 7-deazainosine, as well as the conformational properties of other purine nucleoside analogs, and on inferences derived from calculations about the conformational effect in polynucleotides of removing the 2'-OH group.