ALBERT

Notes: There are 10 unique dinucleotides of double-stranded DNA, but only 8 independent nearest-neighbor energies that occur in circular DNA, as shown by D. M. Gray and I. Tinoco [(1970) Biopolymers 9, 223-244]. We extend that analysis to include end effects, and show that the number of unique dinuoleotide pairs (including ends) is 14. but there are only 12 independent energies. We discuss how these 12 energies (or spectra or any other pairwise additive property) can be measured and displayed, and how they should and should not be compared between experimenters. As an example, we analyzed the recently reported melting curves [M. J. Doktycz et al. (1992) Biopolymers, 32, 849-864] of 16 DNA dumbbells in two different Na+ environments. This analysis reveals a new means for evaluating end effects and the emergence of longer than nearest-neighbor interactions at low salt concentration. © 1992 John Wiley & Sons, Inc.

Notes: CD spectra and melting curves were collected for a 28 base-pair DNA fragment in the form of a DNA dumbbell (linked on both ends by T4 single-strand loops) and the same DNA sequence in the linear form (without end loops). The central 16 base pairs (bp) of the 28-bp duplex region is the poly(pu) sequence: 5′-AGGAAGGAGGAAAGAG-3′. Mixtures of the dumbbell and linear DNAs with the 16-base single-strand sequence 5′-TCCTTCCTCCTTTCTC-3′ were also prepared and studied. At 22°C, CD measurements of the mixtures in 950 mM NaCl, 10 mM sodium acetate, 1 mM EDTA, pH 5.5, at a duplex concentration of 1.8 μM, provided evidence for triplex formation. Spectroscopic features of the triplexes formed with either a dumbbell or linear substrate were quite similar. Melting curves of the duplex molecules alone and in mixtures with the third strand were collected as a function of duplex concentration from 0.16 to 2.15 μM. Melting curves of the dumbbell alone and mixtures with the third strand were entirely independent of DNA concentration. In contrast, melting curves of the linear duplex alone or mixed with the third strand were concentration dependent. At identical duplex concentrations, the dumbbell alone melts ∼20°C higher than the linear duplex. The curve of the linear duplex displayed a significant pretransition probably due to end fraying.On melting curves of mixtures of the dumbbell or linear duplex with the third strand, a low temperature transition with much lower relative hyperchromicity change (∼ 5%) was observed. This transition was attributed to the melting of a new molecular species, e.g., the triplex formed between the duplex and single-strand DNA molecules. In the case of the dumbbell/single-strand mixture, these melting transitions of the triplex and the dumbbell were entirely resolvable. In contrast, the melting transitions of the linear duplex and the triplex overlapped, thereby preventing their clear distinction. To analyze the data, a three-state equilibrium model is presented. The analysis utilizes differences in relative absorbance vs temperature curves of dumbbells (or linear molecules) alone and in mixtures with the third strand. From the model analysis a straightforward derivation of fT(T), the fraction of triplex as a function of temperature, was obtained. Analysis of fT vs temperature curves, in effect melting curves of the triplexes, provided evaluation of thermodynamic parameters of the melting transition. For the triplex formed with the dumbbell substrate, the total transition enthalpy is ΔHT = 118.4 ± 12.8 kcal/mol (7.4 ± 0.8 kcal/mol per triplet unit) and the total transition entropy is ΔST = 344 ± 36.8 cal/K · mol (eu) (21.5 ± 2.3 eu per triple unit). The transition curves of the triplex formed with the linear duplex substrate displayed two distinct regions. A broad pretransition region from fT = 0 to 0.55 and a higher, sharper transition above fT = 0.55. The transition parameters derived from the lower temperature region of the curve are ΔH′T = 44.8 ± 9.6 kcal/mol and ΔS′T = 112 ± 33.6 eu (or ΔH′ = 2.8 ± 0.6 kcal/mol and ΔS′ = 7.0 ± 2.1 eu per triplet). These values are probably too small to correspond to actual melting of the triplex but instead likely reveal effects of end fraying of the duplex substrate on triplex stability. Transition parameters of the upper transition are ΔH′T = 128.0 ± 2.3 kcal/mol and ΔS′T = 379.2 ± 6.4 eu (ΔH′ = 8.0 ± 0.2 kcal/mol and ΔS′ = 23.7 ± 0.4 eu per triplet) in good agreement (within experimental error) with the transition parameters of the triplex formed with the dumbbell substrate. Supposing this upper transition reflects actual dissociation of the third strand from the linear duplex substrate this triplex is comparable in thermodynamic stability to the triplex formed with a dumbbell substrate. Even so, the biphasic melting character of the linear triplex obscures the whole analysis, casting doubt on its absolute reliability. Apparently triplexes formed with a dumbbell substrate offer technical advantages over triplexes formed from linear or hairpin duplex substrates for studies of DNA triplex stability. © 1993 John Wiley & Sons, Inc.

Notes: The preparation and characterization of DNA dumbbells that contain the 16 base-pair duplex sequences 5′G-C-A-T-A-G-A-T-G-A-G-A-A-T-G-C3′ (set 1) and 5′G-C-A-T-C-A-T-C-G-A-T-G-A-T-G-C3 ′(set 2) are reported. The dumbbells of set 1 have the duplex stem nucleated on both ends by Tn (n = 2, 3, 4, 6, 8, 10, and 14) loops. The dumbbells of set 2 have Tn (n = 2, 4, 8, 10) end loops. For the molecules of set 1, effects of end loop size on the electrophoretic mobility, CD and UV absorbance spectra, and cleavage by restriction enzymes, were investigated. Effects of loop size on the CD spectra and restriction enzyme cleavage of the molecules of set 2 were also examined. Optical melting curves of the molecules of set 1 were collected as a function of sodium ion concentration from 30 to 120 mM. These investigations revealed that as loop size decreases, the electrophoretic mobilities, rates of enzyme cleavage, and optical melting temperatures increase. For end loops with at least three T's the observed increases are inversely proportional to loop size. The behavior of the dumbbell with T2 end loops departs from this linear dependence and is anomalous in every experimental context. For molecules with end loops comprised of at least four T's CD spectra were virtually indistinguishable. However, these spectra differed considerably from the CD spectrum of the T2-looped molecule. The CD spectrum of the dumbbell with T3 end loops displayed features common to the dumbbells with larger loops and T2 end loops. Thermodynamic evidence that the terminal G · C base pairs (bps) nucleating the T2 end loops were intact was obtained from a comparison of the melting temperature of this molecule with that of a DNA dumbbell containing the 14 central bps of the set 1 duplex sequence linked instead by end loops comprised of the four base sequence, C-T-T-C. The tm of this latter molecule was determined to be 9°C less than that of the former dumbbell assumed to contain a 16-bp stem and T2 end loops. One- and two-dimensional proton nmr measurements performed on the dumbbell of set 2 with T2 end loops and the same sequence with unpaired T's on the ends also verified formation of the T2 end loops, indicated the duplex residues of the dumbbell were in the B form, and showed that the first 3 bps adjoining the loops are distorted. Results are compared with several recently published relevant studies.

Notes: The preparation and melting of a 16 base-pair duplex DNA linked on both ends by C12H24 (dodecyl) chains is described. Absorbance vs temperature curves (optical melting curves) were measured for the dodecyl-linked molecule and the same duplex molecule linked on the ends instead by T4 loops. Optical melting curves of both molecules were measured in 25, 55, and 85 mM Na+ and revealed, regardless of [Na +], the duplex linked by dodecyl loops is more stable by at least 6°C than the same duplex linked by T4 loops. Experimental curves in each salt environment were analyzed in terms of the two-state and multistate theoretical models. In the two-state, or van't Hoff analysis, the melting transition is assumed to occur in an all-or-none manner. Thus, the only possible states accessible to the molecule throughout the melting transition are the completely intact duplex and the completely melted duplex or minicircle. In the multistate analysis no assumptions regarding the melting transition are required and the statistical occurrence of every possible partially melted state of the duplex is explicitly considered. Results of the analysis revealed the melting transitions of both the dodecyl-linked molecule and the dumbbell with T4 end loops are essentially two state in 25 and 55 mM Na+. In contrast, significant deviations from two-state behavior were observed in 85 m MNa+. From our previously published melting data of DNA dumbbells with Tn end loops where n = 2, 3, 4, 6, 8, 10, 14 [T. M. Paner, M. Amaratunga, and A. S. Benight, (1992) Biopolymers, Vol. 32, pp. 881-892] and the dumbbell with T4 end loops of this study, a plot of d(Tm)/d ln [Na+] was constructed. Extrapolation of this data to n = 1 intersects with the value of d (Tm)/d ln [Na+] obtained for the alkyl-linked dumbbell, suggesting the salt-dependent stability of the alkyl-linked molecule behaves as though the duplex of this molecule were linked by end loops comprised of a single T residue. © 1993 John Wiley & Sons, Inc.

Notes: Many important applications of DNA sequence-dependent hybridization reactions have recently emerged. This has sparked a renewed interest in analytical calculations of sequence-dependent melting stability of duplex DNA. In particular, for many applications it is often desirable to accurately predict the transition temperature, or tm, of short duplex DNA oligomers (∼ 20 base pairs or less) from their sequence and concentration. The thermodynamic analytical method underlying these predictive calculations is based on the nearest-neighbor model. At least 11 sets of nearest-neighbor sequence-dependent thermodynamic parameters for DNA have been published. These sets are compared. Use of the nearest-neighbor sets in predicting tm from the DNA sequence is demonstrated, and the ability of the nearest-neighbor parameters to provide accurate predictions of experimental tm's of short duplex DNA oligomers is assessed. © 1998 John Wiley & Sons, Inc. Biopoly 44: 217-239, 1997

Notes: Theoretical melting curves were calculated for four DNA restriction fragments, 157-257 base pairs (bp), and a series of hypothetical block DNAs with sequences d(C2xAxC2x). d(C2xTxG2x), 5 ≤ x ≤ 40. These DNAs provided a mixture of A·T/G·C sequence distributions with which to investigate the effects of parameters and base-pair changes on the melting of short DNAs. The sensitivity of DNA melting curves to changes in internal loop melting parameters σ and κ was examined. As Expected, theoretical melting curves of short DNAs with a quasirandom base-pair sequence vary little with changes in internal loop parameters. End melting dominates the transition behaviour of these moleucles. This was also observed for the block DNAs up to x = 22. Beyond this length, melting curves are highly sensitive to the internal loop parameters. Sensitivity is also predicted for a 157-bp fragment with a block distribution of A·T and G·C pairs. These results indicate that accurate evaluation of internal loop parameters is possible with short DNAs (100-200 bp) containing a G·C/A·T/G·C block distribution with at least 22 bp in each block. Duplex-to-single-strands dissociation parameters were reevaluated form experimental melting curve data of eight DNA fragments using a least squares fit approach. This analysis confirmed parameter values previously found with a simplified dissociation model. A Priori predictions are made on the effects of base-pair changes on the melting curves of three characterized DNA restriction fragments. Single base-pair changes are predicted to induce small but measurable changes in the melting curves. The characteristics of the altered melting curves depend on the location of the base-pair change.

Notes: Dynamic and static light scattering, CD, and optical melting experiments have been conducted on M13mp19 viral circular single-strand DNA as a function of NaCl concentration. Over the 10,000-fold range in concentration from 100 μM to 1.0M NaCl, the melting curves and CD spectra indicate an increase in base stacking and stability of stacked regions with increased salt concentration. Analysis of dynamic light scattering measurements of the single-strand DNA solutions as a function of K2 from 1.56 to 20 × 1010 cm-2 indicates the collected autocorrelation functions are biexponential, thus revealing the presence of two decaying dynamic components. These components are taken to correspond to (1) global translational motions of the molecular center of mass and (2) motions of the internal molecular subunits. From the evaluated relaxation rates of these components, diffusion coefficients D0 and Dplat are determined. The center of mass translational diffusion coefficient D0, varies in a nonmonotonic manner, by 10%, from 3.75 × 10-8 to 3.39 × 10-8 cm2/s over the NaCl concentration range from 100 μM to 1.0M. Likewise, the radius of gyration RG, obtained from static light scattering experiments, varies by 15% from 699 to 830 Å over the same NaCl range. Dplat, the diffusion coefficient of the internal subunits, displays a different dependence on the NaCl concentration and decreases, by nearly 22% in a titratable fashion, from 12.46 × 10-8 to 10.26 × 10-8 cm2/s, when the salt is increased from 100 μM to 1.0M. A semiquantitative interpretation of these results is provided by analysis of the light scattering data in terms of the circular Rouse-Zimm chain. Rouse-Zimm model parameters are estimated from the experimental results, assuming the circular chains are composed of a fixed number of Gaussian segments, N + 1 = 15. The rms displacement of the internal segments, b, is estimated to be the smallest (442 Å) in 100 mM NaCl. Increases of b to 467 Å in 100 μM and 524 Å in 1.0M NaCl are observed. Meanwhile, the hypothetical friction factor of the internal subunits, f, progressively increases as the NaCl concentration is raised. It is inferred from the evaluated Rouse-Zimm model parameters that both the static flexibility of the circular chain and diffusive displacements of the internal subunits decrease with increases in NaCl concentration from 100 mM to 1.0M. These decreases directly contract the salt-dependent behavior of double-stranded DNA, where greater flexibility is observed when the Na+ concentration is increased. The melting and CD measurements indicate the decrease in flexibility and internal motions is due to the formation of nucleotide stacking in the higher NaCl environments. In 100 μM NaCl, where stacking is highly unfavored, a significant electrostatic contribution to the persistence length likely acts to stiffen the molecule. It appears the observable changes in the internal dynamics of M13mp19 single-strand DNA are associated with increases in base stacking that occur from 100 μM to 1.0M NaCl, which apparently induce relatively small perturbations in the overall global tertiary conformation of the DNA.

Notes: Optical melting transitions of the short DNA hairpins formed from the self-complementary DNA oligomers d[GGATACX4GTATCC] where X = A, T, G, or C measured in 100 mM NaCl are presented. A significant dependence of the melting transitions on loop sequence is observed and transition temperatures, tm, of the hairpins vary from 58.3°C for the T4 loop hairpin to 55.3°C for the A4 loop. A nearest-neighbor sequence-dependent theoretical algorithm for calculating melting curves of DNA hairpins is presented and employed to analyze the experimental melting transitions. Experimental melting curves were fit by adjustment of a single theoretical parameter, Fend(n), the weighting function for a hairpin loop comprised of n single-strand bases. Empirically determined values of Fend(n) provide an evaluation of the free-energy of hairpin loop formation and stability. Effects of heterogeneous nearest-neighbor sequence interactions in the duplex stem on hairpin loop for mation were investigated by evaluating Fend(n) in individual fitting procedures using two of the published sets of nearest-neighbor stacking interactions in DNA evaluated in 100 mM NaCl and given by Wartell and Benight, 1985. In all cases, evaluated values of Fend(n) were obtained that provided exact theoretical predictions of the experimental transitions.Results of the evaluations indicate: (1) Evaluated free-energies of hairpin loop formation are only slightly dependent on loop sequences examined. At the transition temperature, Tm, the free-energy of forming a loop of four bases is approximately equal for T4, G4, or C4 loops and varies from 3.9 to 4.8 kcal/mole depending on the set of nearest-neighbor interactions employed in the evaluations. This result suggests, in light of the observed differences in stability between the T4, G4, and C4 loop hairpins, that sequence-dependent interactions between base residues of the loop are most likely not the source of the enhanced stability of a T4 loop. In contrast, the evaluated free-energy of forming an A4 loop is approximately 400 cal/mole higher for each nearest-neighbor set indicating unfavorable interactions between A bases in a loop-affect loop formation and overall hairpin stability, (2) The absolute value for the free-energy of loop formation at the Tm of each hairpin varies by about 1 kcal/mole depending on the set of nearest-neighbor interactions employed and the relative hierarchy of stability for each loop is conserved for different nearest neighbor sets, (3) The melting process of each hairpin deviates from strict two-state behavior in the order according to loop sequence of T 〉 A 〉 G 〉 C, (4) Results of our analysis are compared with the early work of Scheffler et al., 1970 on the hairpins formed from the copolymer sequences d(T - A)q where q = 9-21. Comparisons with the more recent works a DNA dumbbell (Benight et al., 1988) and the very similar DNA hairpins studied by Senior et al., 1988 are presented.

Notes: Seventeen DNA dumbbells were constructed that have duplex sequences ranging in length from 14 to 18 base pairs linked on the ends by T4 single-strand loops. Fifteen of the molecules have the core duplexes with the sequences 5′G-T-A-T-C-C-(W-X-Y-Z)-G-G-A-T-A-C3′, where (W-X-Y-Z) represents a unique combination of A · T, T · A, G · C, and C · G base pairs. The remaining two molecules have the central sequences (W-X-Y-Z) = A-C and A-C-A-C-A-C. These duplex sequences were designed such that the central sequences include different combinations of the 10 possible nearest-neighbor (n-n) stacks in DNA. In this sense the set of molecules is complete and serves as a model system for evaluating sequence-dependent local stability of DNA. Optical melting curves of the samples were collected in 25, 55, 85, and 115 mM [Na +], and showed, regardless of solvent ionic strength, that the transition temperatures of the dumbbells vary by as much as 14° for different molecules of the set.Results of melting experiments analyzed in terms of a n-n sequence-dependent model allowed evaluation of nine independent linear combinations of the n-n stacking interactions in DNA as a function of solvent ionic strength. Although there are in principle 10 possible different n-n interactions in DNA, these 10 are not linearly independent and therefore can not be uniquely determined. For molecules with ends, there are 9 linearly independent combinations, as opposed to circular or semiinfinite repeating copolymers where only 8 linear combinations of the 10 possible n-n interactions are linearly independent. The n-n interactions are presented as combinations of the deviations from average stacking for the 5′-3′ base-pair doublets, δGi, and reveal several interesting features: (1) Titratable changes in the values of δGi, with changing salt environment are observed. In all salts the most stable unique combination is δG4 = (δGGpC + δGCpG,)/2, and the least stable is the GpG/CpC stack, δG2 = δGGpG/CpC. (2) The χ2 values of the fits of the evaluated δGi's to experimental data increased with decreasing [Na +], suggesting that significant interactions beyond nearest neighbors become more pronounced, particularly at 25 mM Na +. (3) In 85 and 115 mM Na +, where the n-n approximation seems to be most valid, the absolute value of δGi for any n-n stack or average of two n-n stacks is not more than ∼ 220 cal/mole, indicating that deviations from average stacking due to n-n interactions represent about 15% of the total stability of a base pair. The overall thermodynamic stability of DNA is predominantly determined by the sequence content (%G · C). Even though the contribution of n-n interactions to overall stability are intrinsically small, reliable predictions of DNA transition temperatures de novo from sequence can be significantly compromised by cumulative errors in the δGi's. (4) Comparisons of our set of n-n linear combinations evaluated in 115 mM Na+ with various published sets evaluated from melting experiments of long restriction fragments, synthetic polymers, and short oligomers, and those obtained from a reanalysis of published melting data of synthetic polymers, are presented. The analysis reveals a major consensus agreement between n-n free energies evaluated from melting data of restriction fragments and long synthetic repeating copolymers. In contrast, only a minor consensus agreement is obtained between our n-n set and these values or those obtained from melting analysis of a combination of short oligomers and long polymers or those theoretically calculated. Results of these comparisons suggest the values of n-n interactions evaluated from DNA melting curves depend on the length of the melted duplex regions of the DNA molecules that comprise the sample set.

Notes: Dynamic and static intensity light scattering techniques were employed to identify conditions allowing preparation of homogeneous solutions of distinct oligomeric states of RecA protein. These hydrodynamically distinguishable oligomer populations of RecA protein were obtained in homogeneous pure quantities sufficient for physical studies. Results indicate two fairly narrow distributions of RecA oligomers comprised on average of 42 ± 3 and 18 ± 1 RecA monomers. These structures, denoted RecA42 and RecA18, respectively, could be obtained reproducibly in milligram quantities and were stable for at least one week. This enabled reliable characterizations of their hydrodynamic properties by dynamic and total intensity light scattering. These measurements revealed RecA42 had an average translational diffusion coefficient, D20(L) = 8 ± 2 × 10-8 cm2/s, molecular weight, Mr = 1.6 ± 0.1 × 106, and radius of gyration, RG = 465 ± 29 Å. The smaller aggregate, RecA18, had D20(S) = 20.5 ± 2.5 × 10-8 cm2/s. Mr = 7.0 ± 0.4 × 105, and RG = 300 ± 20 Å. Heating RecA18 at 37°C overnight resulted in conversion to a species with hydrodynamic properties indistinguishable from RecA42, called RecA18/42. Conversion of RecA42 to RecA18 occurred almost instantaneously by 50% dilution at 38°C or very slowly with incubation at 4°C for at least 39 days.Self-association reactions of the three starting oligomeric states (RecA18, RecA42, and RecA18/42) induced by MgCl2 were monitored at several temperatures by dynamic light scattering. Results of these experiments provided evaluations of kinetic activation parameters of the self-association reactions. The activation parameters found for each starting oligomeric state of the protein were significantly different, revealing the variable influence of MgCl2 on the activation barriers to RecA self-association. Highly aggregated equilibrium solutions that ultimately form in solutions of each starting oligomeric species, incubated in MgCl2 at 38°C for four days, were characterized by total intensity light scattering. Interpretations of these data in terms of characteristic behavior of random polymers suggests the surface morphologies of these highly associated equilibrium states formed from RecA42 and RecA18/42 are similar but contrast with that of RecA18. Calculated values of the translational diffusion coefficient D0 were obtained for oligomeric structures modeled as helical arrays of connected monomer spheres. Best agreement with experimentally determined diffusion coefficients required that constituent monomer spheres of RecA42 have radii 33-40% larger than the monomer spheres of RecA18. Results suggest the hydrodynamically distinct oligomeric forms of RecA may reside in conformational states with different surface exposure of hydrophobic residues, which results in substantial differences in local solvation/hydration. © 1996 John Wiley & Sons, Inc.