The conformation of single-strand polynucleotides in solution: sedimentation studies of apurinic acid

Physical chemistry of nucleic acids and their complexes

Studies of the physical properties of nucleic acids began almost immediately following the discovery of the DNA structure. Early investigations focused on the stability and specificity of multi-strand polynucleotide complexes, then gradually on their interaction with other molecules, particularly proteins. As molecular and structural biology expanded to provide detailed information about biochemical mechanisms, physical studies eventually acquired the additional constraint that they should be relevant to functioning biological systems. We describe work in our laboratory that began with investigations of relatively simple questions about the role of electrostatic interactions in the stabilization of multi-strand nucleic acid structures, and evolved to studies of chromatin structure in vitro and within the nucleus.

Moving beyond Watson-Crick models of coarse grained DNA dynamics

DNA produces a wide range of structures in addition to the canonical B-form of double-stranded DNA. Some of these structures are stabilized by Hoogsteen bonds. We developed an experimentally parameterized, coarse-grained model that incorporates such bonds. The model reproduces many of the microscopic features of double-stranded DNA and captures the experimental melting curves for a number of short DNA hairpins, even when the open state forms complicated secondary structures. We demonstrate the utility of the model by simulating the folding of a thrombin aptamer, which contains G-quartets, and strand invasion during triplex formation. Our results highlight the importance of including Hoogsteen bonding in coarse-grained models of DNA.

Moments of the End-to-End Vector of a Chain Molecule, Its Persistence and Distribution

The persistence vector a is defined as the configurational average of the chain vector r connecting the ends of the molecule and expressed in a reference frame fixed with respect to the first two skeletal bonds. Moments of second and higher orders in the components of r may readily be calculated for real chains in the rotational isomeric state approximation, and from them the corresponding moments of the vector rho = r - a measured from the terminus of a. Development of the density distribution of rho about a is proposed as an alternative to the customary treatment of the density distribution of r about r = 0 on the assumption that this latter distribution should be (approximately) symmetric. Past difficulties in the analysis of cyclization equilibria involving rings of moderate size, such as occur in single strands of polynucleotide chains, conceivably may be overcome by adoption of this alternative.

Detailed scaling analysis of low-force polyelectrolyte elasticity

Single-molecule force-extension data are typically compared to ideal models of polymer behavior that ignore the effects of self-avoidance. Here, we demonstrate a link between single-molecule data and the scaling pictures of a real polymer. We measure a low-force elasticity regime where the extension L of chemically denatured single-stranded DNA grows as a power law with force f : L approximately f;{gamma} , with gamma approximately 0.60-0.69 . This compares favorably with the "tensile-blob" model of a self-avoiding polymer, which predicts gamma=2/3 . We show that the transition out of the low-force regime is highly salt dependent, and use the tensile-blob model to relate this effect to the salt dependence of the polymer's Kuhn length and excluded-volume parameter. We find that, contrary to the well-known Odijk-Skolnick-Fixman theory, the Kuhn length of single-stranded DNA is linearly proportional to the Debye length of the solution. Finally, we show that the low-force elasticity becomes linear (gamma=1) at approximately 3 M salt, and interpret this as a Theta point of the polymer. At this point, the force-extension data is best described by the wormlike chain model, from which we estimate the bare (nonelectrostatic) persistence length of the polymer to be approximately 0.6 nm .

Mechanism of thermal renaturation and hybridization of nucleic acids: Kramers' process and universality in Watson-Crick base pairing

Renaturation and hybridization reactions lead to the pairing of complementary single-stranded nucleic acids. We present here a theoretical investigation of the mechanism of these reactions in vitro under thermal conditions (dilute solutions of single-stranded chains, in the presence of molar concentrations of monovalent salts and at elevated temperatures). The mechanism follows a Kramers' process, whereby the complementary chains overcome a potential barrier through Brownian motion. The barrier originates from a single rate-limiting nucleation event in which the first complementary base pairs are formed. The reaction then proceeds through a fast growth of the double helix. For the DNA of bacteriophages T7, T4, and phiX174, as well as for Escherichia coli DNA, the bimolecular rate k2 of the reaction increases as a power law of the average degree of polymerization <N> of the reacting single-strands: k2 is proportional to <N> alpha. This relationship holds for 100 < or = <N> < or = 50,000 with an experimentally determined exponent alpha = 0.51 +/- 0.01. The length dependence results from a thermodynamic excluded-volume effect. The reacting single-stranded chains are predicted to be in universal good solvent conditions, and the scaling law is determined by the relevant equilibrium monomer contact probability. The value theoretically predicted for the exponent is alpha = 1 - nutheta2, where nu is Flory's swelling exponent (nu approximately equal 0.588), and theta2 is a critical exponent introduced by des Cloizeaux (theta2 approximately equal 0.82), yielding alpha = 0.52 +/- 0.01, in agreement with the experimental results.

The flexibility of locally melted DNA

Protein-bound duplex DNA is often bent or kinked. Yet, quantification of intrinsic DNA bending that might lead to such protein interactions remains enigmatic. DNA cyclization experiments have indicated that DNA may form sharp bends more easily than predicted by the established worm-like chain (WLC) model. One proposed explanation suggests that local melting of a few base pairs introduces flexible hinges. We have expanded this model to incorporate sequence and temperature dependence of the local melting, and tested it for three sequences at temperatures from 23°C to 42°C. We find that small melted bubbles are significantly more flexible than double-stranded DNA and can alter DNA flexibility at physiological temperatures. However, these bubbles are not flexible enough to explain the recently observed very sharp bends in DNA.

Melting of persistent double-stranded polymers

Motivated by recent DNA-pulling experiments, we revisit the Poland-Scheraga model of melting a double-stranded polymer. We include distinct bending rigidities for both the double-stranded segments and the single-stranded segments forming a bubble. There is also bending stiffness at the branch points between the two segment types. The transfer matrix technique for single persistent chains is generalized to describe the branching bubbles. Properties of spherical harmonics are then exploited in truncating and numerically solving the resulting transfer matrix. This allows efficient computation of phase diagrams and force-extension curves (isotherms). While the main focus is on exposition of the transfer matrix technique, we provide general arguments for a reentrant melting transition in stiff double strands. Our theoretical approach can also be extended to study polymers with bubbles of any number of strands, with potential applications to molecules such as collagen.

Simulation of polymer translocation through protein channels

A modeling algorithm is presented to compute simultaneously polymer conformations and ionic current, as single polymer molecules undergo translocation through protein channels. The method is based on a combination of Langevin dynamics for coarse-grained models of polymers and the Poisson-Nernst-Planck formalism for ionic current. For the illustrative example of ssDNA passing through the alpha-hemolysin pore, vivid details of conformational fluctuations of the polymer inside the vestibule and beta-barrel compartments of the protein pore, and their consequent effects on the translocation time and extent of blocked ionic current are presented. In addition to yielding insights into several experimentally reported puzzles, our simulations offer experimental strategies to sequence polymers more efficiently.

Assessment of adenyl residue reactivity within model nucleic acids by surface enhanced Raman spectroscopy

We rank the reactivity of the adenyl residues (A) of model DNA and RNA molecules with electropositive subnano size [Ag]n+ sites as a function of nucleic acid primary sequences and secondary structures and in the presence of biological amounts of Cl- and Na+ or Mg2+ ions. In these conditions A is markedly more reactive than any other nucleic acid bases. A reactivity is higher in ribo (r) than in deoxyribo (d) species [pA>pdA and (pA)n>>(pdA)n]. Base pairing decreases A reactivity in corresponding duplexes but much less in r than in d. In linear single and paired dCAG or dGAC loci, base stacking inhibits A reactivity even if A is bulged or mispaired (A.A). dA tracts are highly reactive only when dilution prevents self-association and duplex structures. In d hairpins the solvent-exposed A residues are reactive in CAG and GAC triloops and even more in ATC loops. Among the eight rG1N2R3A4 loops, those bearing a single A (A4) are the least reactive. The solvent-exposed A2 is reactive, but synergistic structural transitions make the initially stacked A residues of any rGNAA loop much more reactive. Mg2+ cross-bridging single strands via phosphates may screen A reactivity. In contrast d duplexes cross-bridging enables "A flipping" much more in rA.U pairs than in dA.T. Mg2+ promotes A reactivity in unpaired strands. For hairpins Mg2+ binding stabilizes the stems, but according to A position in the loops, A reactivity may be abolished, reduced, or enhanced. It is emphasized that not only accessibility but also local flexibility, concerted docking, and cation and anion binding control A reactivity.

Mechanical studies of single ribosome/mRNA complexes

Methodology was developed for specifically anchoring Escherichia coli 70S ribosomes onto a chemically modified, cysteine-reactive glass surface. Immobilized ribosomes maintain the capability of binding a polyuridylic acid (poly(U)) template, enabling investigation of mechanical properties of individual ribosome-poly(U) complexes using laser tweezers. Streptavidin-coated polystyrene microspheres bound specifically to the biotinylated 3' end of long (up to 10,000 bases) poly(U) strands. A novel optical method was built to control the position of the laser trap along the microscope optical axis at 2 nm resolution, facilitating measurement of the force-extension relationship for poly(U). Some immobilized ribosome-poly(U) complexes supported 100 pN of force applied at the 3' end of the mRNA. Binding of N-acetylated Phe-tRNA(Phe), an analog of the initiator fMet-tRNA(Met), enhanced the population of complexes that could withstand high forces. The persistence length of poly(U) RNA homopolymer, modeled as a worm-like chain, was found to be 0.79 +/- 0.05 nm and the backbone elasticity was 900 +/- 140 pN, similar to values for single-stranded DNA.

Gapped DNA and cyclization of short DNA fragments

We use the cyclization of small DNA molecules, approximately 200 bp in length, to study conformational properties of DNA fragments with single-stranded gaps. The approach is extremely sensitive to DNA conformational properties and, being complemented by computations, allows a very accurate determination of the fragment's conformational parameters. Sequence-specific nicking endonucleases are used to create the 4-nt-long gap. We determined the bending rigidity of the single-stranded region in the gapped DNA. We found that the gap of 4 nt in length makes all torsional orientations of DNA ends equally probable. Our results also show that the gap has isotropic bending rigidity. This makes it very attractive to use gapped DNA in the cyclization experiments to determine DNA conformational properties, since the gap eliminates oscillations of the cyclization efficiency with the DNA length. As a result, the number of measurements is greatly reduced in the approach, and the analysis of the data is greatly simplified. We have verified our approach on DNA fragments containing well-characterized intrinsic bends caused by A-tracts. The obtained experimental results and theoretical analysis demonstrate that gapped-DNA cyclization is an exceedingly sensitive and accurate approach for the determination of DNA bending.

Probing single-stranded DNA conformational flexibility using fluorescence spectroscopy

Single-stranded DNA (ssDNA) is an essential intermediate in various DNA metabolic processes and interacts with a large number of proteins. Due to its flexibility, the conformations of ssDNA in solution can only be described using statistical approaches, such as flexibly jointed or worm-like chain models. However, there is limited data available to assess such models quantitatively, especially for describing the flexibility of short ssDNA and RNA. To address this issue, we performed FRET studies of a series of oligodeoxythymidylates, (dT)N, over a wide range of salt concentrations and chain lengths (10 < or = N < or = 70 nucleotides), which provide systematic constraints for testing theoretical models. Unlike in mechanical studies where available ssDNA conformations are averaged out during the time it takes to perform measurements, fluorescence lifetimes may act here as an internal clock that influences fluorescence signals depending on how fast the ssDNA conformations fluctuate. A reasonably good agreement could be obtained between our data and the worm-like chain model provided that limited relaxations of the ssDNA conformations occur within the fluorescence lifetime of the donor. The persistence length thus estimated ranges from 1.5 nm in 2 M NaCl to 3 nm in 25 mM NaCl.

Stereoselectivity in DNA-templated organic synthesis and its origins

DNA-templated synthesis is a surprisingly general strategy for controlling chemical reactivity that enables synthetic products to be manipulated in ways previously available only to biological macromolecules. The chiral nature of the DNA template raises the possibility that DNA-templated synthesis can proceed stereoselectively. Here, we show that DNA-templated substitution reactions can exhibit stereoselectivity without the assistance of chiral groups other than those present in DNA. By characterizing changes in stereoselectivity as a result of systematic changes in the structure of the template-reagent complexes, we begin to reveal the origins of the observed stereoselectivity. We propose that the conformations of nucleotides adjacent to the reactants are largely responsible for stereoselectivity. Indeed, template and reagent sequences that can adopt either a left-handed Z-form DNA helix or a normal right-handed B-form DNA helix generate opposite stereoselectivities in the Z-form and B-form even though they share the same covalent structure and the same absolute stereochemistry. Our findings establish ways in which the chirality of an information carrier can be transmitted to the stereochemistry of encoded products through templated synthesis.

A physical origin for functional domain structure in nucleic acids as evidenced by cross-linking entropy: I

A global strategy for estimating the entropy of long sequences of RNA is proposed to help improve the predictive capacity of RNA secondary structure dynamic programming algorithm (DPA) free energy (FE) minimization methods. These DPA strategies only consider the effects that occur in the immediate (nearest neighbor) vicinity of a given base pair (bp) in a secondary structure plot. They are therefore defined as nearest-neighbor secondary structure (NNSS) strategies. The new approach utilizes the statistical properties of the Gaussian polymer chain model to introduce both local and global contributions to the entropy of a given secondary structure. These entropic contributions are primarily a function of the persistence length of the RNA. Limits on the domain size are strongly suggested by this model and these limits are a function of both the length and the percentage of bp enclosed within a given domain. The model generalizes the penalties found in the NNSS algorithms. The approach considers the importance of flexibility in the folding and stability of RNA by considering the role of the persistence length in a biopolymer structure. The theory also suggests that molecular machinery may also take advantage of this global entropic effect to bring about catalytic effects. The applications can also be extended to protein structure calculations with some additional considerations.

Synthesis of abasic locked nucleic acid and two seco-LNA derivatives and evaluation of their hybridization properties compared with their more flexible DNA counterparts

To investigate the structural basis of the unique hybridization properties of LNA (locked nucleic acid) three novel LNA derivatives with modified carbohydrate parts were synthesized and evaluated with respect to duplex stabilities. The abasic LNA monomer (X(L), Figure 1) with the rigid carbohydrate moiety of LNA but no nucleobase attached showed no enhanced duplex stabilities compared to its more flexible abasic DNA counterpart (X, Figure 1). These results suggest that the exceptional hybridization properties of LNA primarily originate from improved intrastrand nucleobase stacking and not backbone preorganization. Two monocyclic seco-LNA derivatives, obtained by cleavage of the C1'-O4' bond of an LNA monomer or complete removal of the O4'-furanose oxygen atom (Z(L) and dZ(L), respectively, Figure 1), were compared to their acyclic DNA counterpart (Z, Figure 1). Even though they are more constrained than Z, the seco-LNA derivatives Z(L) and dZ(L) destabilize duplex formation even more than the flexible seco-DNA monomer Z.

Driven polymer translocation through a narrow pore

Motivated by experiments in which a polynucleotide is driven through a proteinaceous pore by an electric field, we study the diffusive motion of a polymer threaded through a narrow channel with which it may have strong interactions. We show that there is a range of polymer lengths in which the system is approximately translationally invariant, and we develop a coarse-grained description of this regime. From this description, general features of the distribution of times for the polymer to pass through the pore may be deduced. We also introduce a more microscopic model. This model provides a physically reasonable scenario in which, as in experiments, the polymer's speed depends sensitively on its chemical composition, and even on its orientation in the channel. Finally, we point out that the experimental distribution of times for the polymer to pass through the pore is much broader than expected from simple estimates, and speculate on why this might be.

Flexibility of single-stranded DNA: use of gapped duplex helices to determine the persistence lengths of poly(dT) and poly(dA)

The forces responsible for the formation and stabilization of secondary and higher-order nucleic acid structure can be more fully understood once the sequence-dependent properties (e.g. intrinsic rigidity, effective rise) of the component single-stranded species are well-defined. Knowledge of the conformations of the single-stranded polymers is also important for the development of better polyelectrolyte models for various structural or strand-dissociation reactions. However, there is at present little quantitative information regarding the sequence dependence of either rise or rigidity in single-stranded DNA or RNA polymers. To address this issue, we describe a form of transient electric birefringence (TEB) measurement in which the rotational decay times (taugap) of DNA molecules possessing central, single-stranded regions (gaps) are compared with the corresponding times (taudplx) for duplex control molecules of the same length (in nucleotides per strand) as the continuous strand in the "gapped duplex". For magnesium ion concentrations above 1-2 mM, the tau ratios ( identical withtaugap/taudplx) for the gapped duplexes reach plateau values, above which no further change in tau ratio is observed; values for the persistence length (P) and internucleotide spacing (h) of the gap sequences are obtained from the experimental tau ratios. For dTn, the permissible ranges of P and h are 20-30 A and 5-7 A per nucleotide (nt), respectively, with optimal values of 31 A (P) and 5. 2 A/nt (h). For dAn, the persistence length for the low temperature (4 degreesC), stacked form is 78 (+/-8) A for a helix rise of 3.2 A/nt. One significant advantage of the current method over previous approaches is the use of short (80-100 nt) molecules, thus facilitating the production of various gap sequences through synthetic means.

Overstretching B-DNA: the elastic response of individual double-stranded and single-stranded DNA molecules

Single molecules of double-stranded DNA (dsDNA) were stretched with force-measuring laser tweezers. Under a longitudinal stress of approximately 65 piconewtons (pN), dsDNA molecules in aqueous buffer undergo a highly cooperative transition into a stable form with 5.8 angstroms rise per base pair, that is, 70% longer than B form dsDNA. When the stress was relaxed below 65 pN, the molecules rapidly and reversibly contracted to their normal contour lengths. This transition was affected by changes in the ionic strength of the medium and the water activity or by cross-linking of the two strands of dsDNA. Individual molecules of single-stranded DNA were also stretched giving a persistence length of 7.5 angstroms and a stretch modulus of 800 pN. The overstretched form may play a significant role in the energetics of DNA recombination.

A genetic algorithm based molecular modeling technique for RNA stem-loop structures

A new modeling technique for arriving at the three dimensional (3-D) structure of an RNA stem-loop has been developed based on a conformational search by a genetic algorithm and the following refinement by energy minimization. The genetic algorithm simultaneously optimizes a population of conformations in the predefined conformational space and generates 3-D models of RNA. The fitness function to be optimized by the algorithm has been defined to reflect the satisfaction of known conformational constraints. In addition to a term for distance constraints, the fitness function contains a term to constrain each local conformation near to a prepared template conformation. The technique has been applied to the two loops of tRNA, the anticodon loop and the T-loop, and has found good models with small root mean square deviations from the crystal structure. Slightly different models have also been found for the anticodon loop. The analysis of a collection of alternative models obtained has revealed statistical features of local variations at each base position.

Electrophoretic evidence that single-stranded regions of one or more nucleotides dramatically increase the flexibility of DNA

The influence of single-stranded nicks and gaps on the flexibility of DNA has been investigated by subjecting to gel electrophoresis sets of molecules containing single-stranded regions of defined position and length. The DNA molecules were produced by ligating together synthetic oligomers that contained either nicks or single-stranded gaps of 1-4 nucleotides; the oligomer repeat lengths were 20, 21, 22, 23, or 26 bp, in order to produce nicks or gaps that were either in- or out-of-phase with the helix repeat of DNA. Nick-containing DNA molecules displayed nearly normal electrophoretic behavior, with maximum reductions in gel mobility (41 degrees C; 12% polyacrylamide gels) of approximately 10% for 230-bp molecules containing 10 nicks. In contrast, molecules containing gaps of 2-4 nucleotides demonstrated dramatic reductions in mobility, approaching one-half of the values of their full-duplex counterparts; molecules containing 1-nucleotide gaps displayed intermediate behavior. The observed (relative) mobilities of molecules containing gaps of more than 1 nucleotide were remarkably insensitive to temperature and to the presence of magnesium ions in the electrophoresis buffer. The central conclusion of the current study is that single-stranded gaps represent points of swivel-like character, whereas nicks retain much of the rigid character of double-helical DNA.

Molecular modelling of the 3-D structure of RNA tetraloops with different nucleotide sequences

One surprisingly common element of RNA secondary structure consists of a hairpin capped by a four-base loop (or the tetraloop). Recently the 3-D structures of two RNA-tetraloops have been determined by NMR-studies. Both structures have a similar architecture: the first and the last bases of the loop form a hydrogen bonded pair which is stacked on the stem base pair. We have analysed the ability of tetraloops, with the other combinations of the first and the fourth bases, to adopt such a 'diloop' conformation using computer modelling. The analysis has shown that the 'diloop' conformation has many covalent and steric constraints which give a possibility for reliable structural predictions. As a result, a set of the tetraloop 3-D structures in which hydrogen bonded pairing of the first and the last bases does not cause covalent and steric hindrances has been selected. In most cases several predicted 3-D structures corresponded to one tetraloop sequence. Taking into consideration the folding pathway of RNA hairpins we have resolved this ambiguity and predicted the most probable 3-D structure for every possible nucleotide sequence of the tetraloop. On the basis of these results a conclusion has been drawn on the possible reasons of the tetraloop phylogenetic preference.

Reactivity of antibodies to guanosine modified by the carcinogen N-acetoxy-N-2-acetylaminofluorene

N-(guanosin-8-yl) acetylaminofluorene (Guo-AAF) was prepared by the reaction of N-acetoxy-N-2-acetylaminofluorene (AAAF) and guanosine. Antibodies to Guo-AAF were elicited in rabbits by immunization with bovine serum albumin-Guo-AAF conjugate. The antibodies were purified by affinity chromatography on a Sepharose-Guo-AAF column. The reactivity of these antibodies towards several ligands was studied by radioimmunoassay. The antibodies have the same affinity for double stranded DNA-AAF and single stranded DNA-AAF. Thus the geometry of the regions of DNA substituted by AAF residues is the same in native and denatured DNA. The affinity of the antibodies is smaller for DNA-AAF than for Guo-AAF. This can be due in part to the stacking of AAF residues with the adjacent bases as shown by the study of the interactions between the antibodies and AAF-oligonucleotides. The circular dichroism spectra of AAF-oligonucleotides bound to the antibodies are reported.

A model for the single stranded random coil form of polydeoxyadenylic acid from minimum energy conformations of the dimeric subunit

The minimum energy conformations of dApdA have been examined for their suitability as buildings blocks of the single stranded coil form of polynucleotides. Calculations of the characteristic ratio C difference = less than ro greater than 2/n liter2 were made for a polymer generated from all the low energy conformers, as well as for selected combinations. A polymer composed of a conformer with omega', omega = t*,g+,(skewed) psi = t, C-(2)-endo type pucker, in combination with the 'B' form, has a C difference equal to that observed in coils of apurinic acid (6) when the fraction of 'B' form conformers is approximately 25% and approximately 91%. The t*,g+ conformer is the second lowest energy form in the C-(2)-endo puckering domain, following the 'B' form.

RNA structure

This review is concerned primarily with the physical structure and changes in the structure of RNA molecules. It will be evident that we have not attempted comprehensive coverage of what amounts to a vast literature. We have tried to stay away from particular areas that have been recently reviewed elsewhere. Citations to and information from them are included, however, so that access to the literature is available. Much of what we treat in depth deals with the crystal structures and solution behaviour of model RNA compounds, including synthetic polymers and molecular fragments such as dinucleoside phosphates. Sequence data on natural RNA are cited, but not in detail. Similarly, apart from tRNA, natural RNAs the structural determinations of which are presently not so far advanced, are not dwelt upon. We have tried to present in detail the available structural data with scaled drawings that permit facile comparisons of molecular geometries.