DNA and RNA oligomer thermodynamics: the effect of mismatched bases on double-helix stability

How sequence alterations enhance the stability and delay expansion of DNA triplet repeat domains

DNA sequence alterations within DNA repeat domains inexplicably enhance the stability and delay the expansion of interrupted repeat domains. Here we propose mechanisms that rationalise such unanticipated outcomes. Specifically, we describe how interruption of a DNA repeat domain restricts the ensemble space available to dynamic, slip out, repeat bulge loops by introducing energetic barriers to loop migration. We explain how such barriers arise because some possible loop isomers result in energetically costly mismatches in the duplex portion of the repeat domain. We propose that the reduced ensemble space is the causative feature for the observed delay in repeat DNA expansion. We further posit that the observed loss of the interrupting repeat in some expanded DNAs reflects the transient occupation of loop isomer positions that result in a mismatch in the duplex stem due to 'leakiness' in the energy barrier. We propose that if the lifetime of such a low probability event allows for recognition by the mismatch repair system, then 'repair' of the repeat interruption can occur; thereby rationalising the absence of the interruption in the final expanded DNA 'product.' Our proposed mechanistic pathways provide reasoned explanations for what have been described as 'puzzling' observations, while also yielding insights into a biomedically important set of coupled genotypic phenomena that map the linkage between DNA origami thermodynamics and phenotypic disease states.

Nearest neighbor rules for RNA helix folding thermodynamics: improved end effects

Nearest neighbor parameters for estimating the folding stability of RNA secondary structures are in widespread use. For helices, current parameters penalize terminal AU base pairs relative to terminal GC base pairs. We curated an expanded database of helix stabilities determined by optical melting experiments. Analysis of the updated database shows that terminal penalties depend on the sequence identity of the adjacent penultimate base pair. New nearest neighbor parameters that include this additional sequence dependence accurately predict the measured values of 271 helices in an updated database with a correlation coefficient of 0.982. This refined understanding of helix ends facilitates fitting terms for base pair stacks with GU pairs. Prior parameter sets treated 5′GGUC3′ paired to 3′CUGG5′ separately from other 5′GU3′/3′UG5′ stacks. The improved understanding of helix end stability, however, makes the separate treatment unnecessary. Introduction of the additional terms was tested with three optical melting experiments. The average absolute difference between measured and predicted free energy changes at 37°C for these three duplexes containing terminal adjacent AU and GU pairs improved from 1.38 to 0.27 kcal/mol. This confirms the need for the additional sequence dependence in the model.

An asymmetric mesoscopic model for single bulges in RNA

Simple one-dimensional DNA or RNA mesoscopic models are of interest for their computational efficiency while retaining the key elements of the molecular interactions. However, they only deal with perfectly formed DNA or RNA double helices and consider the intra-strand interactions to be the same on both strands. This makes it difficult to describe highly asymmetric structures such as bulges and loops and, for instance, prevents the application of mesoscopic models to determine RNA secondary structures. Here we derived the conditions for the Peyrard-Bishop mesoscopic model to overcome these limitations and applied it to the calculation of single bulges, the smallest and simplest of these asymmetric structures. We found that these theoretical conditions can indeed be applied to any situation where stacking asymmetry needs to be considered. The full set of parameters for group I RNA bulges was determined from experimental melting temperatures using an optimization procedure, and we also calculated average opening profiles for several RNA sequences. We found that guanosine bulges show the strongest perturbation on their neighboring base pairs, considerably reducing the on-site interactions of their neighboring base pairs.

Synthesis of PNA oligoether conjugates

Several different approaches have been explored for conjugation of oligoethers to PNA with internally or N-terminal placed diaminopropionic acid residues. Single and double conjugation of 2-(2-(2-aminoethoxy)ethoxy)ethanol was obtained using carbonyldimidazole. Using a post PNA-assembly coupling procedure the building block 2-(2-(2-(benzoyloxy)ethoxy)ethoxy)acetic acid multiple attachment of 2-(2-(2-hydroxyethoxy)ethoxy)acetyl groups to both N-terminal and β-amino groups of inserted diaminopropionic acids residues was achieved. Use of a new oligoether functionalized amino acid allows inclusion of oligoether conjugates during on-line machine assisted synthesis which also allowed combination of methods for attachment of different oligoethers and co-conjugation of neocuproine as well as conjugation of an aminosugar.

Electrochemical measurements of DNA melting on surfaces

Thermal denaturation, or melting, measurements are a classic technique for analysis of thermodynamics of nucleic base driven associations in solution, as well as of interactions between nucleic acids and small molecule ligands such as drugs or carcinogens. Performed on surface-immobilized DNA films, this well-established technique can help understand how energetics of surface hybridization relate to those in solution, as well as provide high-throughput platforms for screening of small molecule ligands. Here we describe methods for measuring DNA melting transitions at solid/liquid interfaces with focus on the role of immobilization chemistry, including a common "immobilization-through-self-assembly" approach that is effective at moderate temperatures, and a thermo-stable approach based on polymer-supported DNA monolayers that can be used at elevated temperatures. We also discuss conditions necessary for reversible measurements, as signified by superimposition of the association (cooling) and dissociation (heating) transitions of immobilized DNA strands.

Electrochemically amplified molecular beacon biosensor for ultrasensitive DNA sequence-specific detection of Legionella sp

An electrochemically amplified molecular beacon (EAMB) biosensor is constructed using thiolated hairpin DNA-ferrocene probes on gold electrode. The switching from "on" to "off" states of individual probes in the presence of complementary DNA target influences the electrode potential, besides the current, owing to changes in surface density of the electroactive hairpin DNA-ferrocene probes. The EAMB biosensor demonstrates linear range over 8 orders of magnitude with ultrasensitive detection limit of 2.3 × 10(-14)M for the quantification of a 21-mer DNA sequence. Its applicability is tested against PCR amplicons derived from genomic DNA of live Legionella pneumophila. Excellent specificity down to one and three nucleotides mismatches in another strain of L. pneumophila and a different bacterium species, respectively, is demonstrated.

3D maps of RNA interhelical junctions

More than 50% of RNA secondary structure is estimated to be A-form helices, which are linked together by various junctions. Here we describe a protocol for computing three interhelical Euler angles describing the relative orientation of helices across RNA junctions. 5' and 3' helices, H1 and H2, respectively, are assigned based on the junction topology. A reference canonical helix is constructed using an appropriate molecular builder software consisting of two continuous idealized A-form helices (iH1 and iH2) with helix axis oriented along the molecular Z-direction running toward the positive direction from iH1 to iH2. The phosphate groups and the carbon and oxygen atoms of the sugars are used to superimpose helix H1 of a target interhelical junction onto the corresponding iH1 of the reference helix. A copy of iH2 is then superimposed onto the resulting H2 helix to generate iH2'. A rotation matrix R is computed, which rotates iH2' into iH2 and expresses the rotation parameters in terms of three Euler angles α(h), β(h) and γ(h). The angles are processed to resolve a twofold degeneracy and to select an overall rotation around the axis of the reference helix. The three interhelical Euler angles define clockwise rotations around the 5' (-γ(h)) and 3' (α(h)) helices and an interhelical bend angle (β(h)). The angles can be depicted graphically to provide a 'Ramachandran'-type view of RNA global structure that can be used to identify unusual conformations as well as to understand variations due to changes in sequence, junction topology and other parameters.

Stability of a surface-bound oligonucleotide duplex inferred from molecular dynamics: a study of single nucleotide defects using DNA microarrays

Microarray technology uses the sequence dependent hybridization (binding) affinity of surface-bound oligonucleotide strands for the quantification of complex nucleic acid mixtures. In spite of its huge potential in life science and medicine, microarray oligonucleotide hybridization remains far from being understood. Taking advantage of microarray combinatorial possibilities we show that, although surface bound, the hybridization affinities of single-base mismatched oligonucleotides can be derived from first principles using parameters from bulk.

Position dependent mismatch discrimination on DNA microarrays - experiments and model

Background

The propensity of oligonucleotide strands to form stable duplexes with complementary sequences is fundamental to a variety of biological and biotechnological processes as various as microRNA signalling, microarray hybridization and PCR. Yet our understanding of oligonucleotide hybridization, in particular in presence of surfaces, is rather limited. Here we use oligonucleotide microarrays made in-house by optically controlled DNA synthesis to produce probe sets comprising all possible single base mismatches and base bulges for each of 20 sequence motifs under study.

Results

We observe that mismatch discrimination is mostly determined by the defect position (relative to the duplex ends) as well as by the sequence context. We investigate the thermodynamics of the oligonucleotide duplexes on the basis of double-ended molecular zipper. Theoretical predictions of defect positional influence as well as long range sequence influence agree well with the experimental results.

Conclusion

Molecular zipping at thermodynamic equilibrium explains the binding affinity of mismatched DNA duplexes on microarrays well. The position dependent nearest neighbor model (PDNN) can be inferred from it. Quantitative understanding of microarray experiments from first principles is in reach.

Investigation of potential RNA bulge stabilizing elements

As a part of our interest in recognition and cleavage of RNA we carried out thermal melting studies with the aim of screening a number of simple oligonucleotide modifications for their potential as modifying elements for RNA bulge stabilizing oligonucleotides. A specific model system from our studies on oligonucleotide-based artificial nuclease (OBAN) systems was chosen and the bulge size was varied from one to five nucleotides. Introduction of single 2'-modified nucleoside moieties (2'-O-methyl, 2'-deoxy and 2'-deoxy-2'-amino) with different conformational preferences adjacent to the bulge revealed that a higher preference for the north conformers gave more stable bulges across the whole range of bulge sizes. Changing a bulge closing a G-U wobble base pair to a G-C pair resulted in the interesting observation that, although the fully complementary complex and small bulges were highly stabilized, there was little difference in the stability of the larger bulges. The wobble base pair even gave a slight stabilization of the 5 nt bulge system. Introduction of a uridine C-5 linker with a single ammonium group was clearly bulge stabilizing (DeltaT(m) + 4.6 to + 5.4 degrees C for the three most stabilized bulges), although with limited selectivity for different bulge sizes since the fully complementary duplex was also stabilized. Introduction of a naphthoyl group on a 2'-aminolinker mostly gave a destabilizing effect, while introduction of a 5-aminoneocuproine moiety on the same linker resulted in stabilization of all bulges, in particular those with two or four unpaired nucleotides (DeltaT(m) + 3.6 and + 2.9 degrees C respectively). The aromatic groups destabilize the fully complementary duplex, resulting in higher selectivity towards stabilization of bulges. A combination of the studied partial element should be suitable for future designs of modified oligonucleotides that, apart from standard base pairing, can also provide additional non-Watson-Crick recognition of RNA.

Stabilisation of RNA bulges by oligonucleotide complements containing an adenosine analogue

Incorporation of 2'-deoxy-2'-beta-(1-naphthylmethyl)tubercidin into an oligodeoxyribonucleotide mostly has little or a slightly negative effect on the T(m) values of complexes with DNA complements. With the same naphthylmethyl-substituted nucleoside at the 3'-end of a 2'-O-methyloligoribonucleotide, however, a stabilisation of 1-2 degrees C in the corresponding complexes with both DNA and RNA is observed. When the target sequence is an RNA fragment forming a two- or three-nucleotide bulge, complexes with (naphthylmethyl)tubercidin-modified oligodeoxyribonucleotides, as well as with the corresponding 2'-O-methyloligoribonucleotides, give stabilisations of 1-2 degrees C for the three-nucleotide bulge and of almost 4 degrees C for the two-nucleotide bulge. This stabilisation is specific to RNA, since the corresponding complexes with the DNA fragments do not display this effect. Thus, the (naphthylmethyl)tubercidin-containing oligonucleotides are the first reported oligonucleotide modifications that specifically stabilise bulged RNA.

Correlation between sequence-dependent glycosylase repair and the thermal stability of oligonucleotide duplexes containing 1, N6-ethenoadenine

Previous experiments on DNA sequence context reported that base modification, replication, and repair are affected by the nature of neighbor bases. We now report that repair by mammalian alkylpurine-DNA-N-glycosylases (APNG) of 15-mer oligonucleotides with a central 1,N6-ethenoadenine (epsilonA), flanked by 5' and 3' tandem bases, is also highly sequence dependent. Oligonucleotides with the central sequences -GGepsilonAGG- or -CCepsilonACC- are repaired 3-5-fold more efficiently than those containing -AAepsilonAAA- or -TTepsilonATT- when using human or mouse APNG. Melting curves of the same duplexes showed that oligomers with G.C/C. G neighbors were less denatured than those with A.T/T.A neighbors at 37 degreesC. This sequence-dependent difference in denaturation correlates with the relative thermodynamic stability of oligomers with G.C/C.G or A.T/T.A neighbors. The dependence of repair on thermal stability was confirmed by enzyme reactions performed over 0-45 degreesC. Under these conditions, repair of epsilonA flanked by G.C/C.G was dramatically increased at 37 degreesC with continuous increase up to 45 degreesC, in contrast to that with flanking A.T/T. A pairs, which was in agreement with the degree of denaturation of these duplexes. These results indicate that the thermodynamic stability conferred by base pairs flanking epsilonA plays an essential role in maintaining the integrity of the duplex structure which is necessary for repair.

Thermodynamic parameters for an expanded nearest-neighbor model for formation of RNA duplexes with Watson-Crick base pairs

Improved thermodynamic parameters for prediction of RNA duplex formation are derived from optical melting studies of 90 oligoribonucleotide duplexes containing only Watson-Crick base pairs. To test end or base composition effects, new sets of duplexes are included that have identical nearest neighbors, but different base compositions and therefore different ends. Duplexes with terminal GC pairs are more stable than duplexes with the same nearest neighbors but terminal AU pairs. Penalizing terminal AU base pairs by 0.45 kcal/mol relative to terminal GC base pairs significantly improves predictions of DeltaG degrees37 from a nearest-neighbor model. A physical model is suggested in which the differential treatment of AU and GC ends accounts for the dependence of the total number of Watson-Crick hydrogen bonds on the base composition of a duplex. On average, the new parameters predict DeltaG degrees37, DeltaH degrees, DeltaS degrees, and TM within 3.2%, 6.0%, 6.8%, and 1.3 degreesC, respectively. These predictions are within the limit of the model, based on experimental results for duplexes predicted to have identical thermodynamic parameters.

The ups and downs of nucleic acid duplex stability: structure-stability studies on chemically-modified DNA:RNA duplexes

In an effort to discover novel oligonucleotide modifications for antisense therapeutics, we have prepared oligodeoxyribonucleotides containing more than 200 different modifications and measured their affinities for complementary RNA. These include modifications to the heterocyclic bases, the deoxy-ribose sugar and the phosphodiester linkage. From these results, we have been able to determine structure-activity relationships that correlate hybridization affinity with changes in oligonucleotide structure. Data for oligonucleotides containing modified pyrimidine nucleotides are presented. In general, modifications that resulted in the most stable duplexes contained a heteroatom at the 2'-position of the sugar. Other sugar modifications usually led to diminished hybrid stability. Most backbone modifications that led to improved hybridization restricted backbone mobility and resulted in an A-type sugar pucker for the residue 5'to the modified internucleotide linkage. Among the heterocycles, C-5-substituted pyrimidines stood out as substantially increasing duplex stability.

Influence of defects on the electrophoretic, thermodynamic and dielectric properties of a 21 base pair DNA in solution

The thermodynamic and dielectric properties of a 21 base pair DNA have been evaluated and compared with those of samples with some defects. In particular, fragments in which the absence of a phosphate group (nick) or of two nucleotides (gap) causes chain interruptions were studied. Measurements of ultraviolet absorption as a function of temperature at different oligomer concentrations and at various ionic strengths were performed. Dielectric spectroscopy at radiofrequencies (1 MHz-1 GHz) was applied on aqueous solutions of the duplexes at 5 degrees C, where the solutes are thermally stable. Dielectric dispersions with 30-40 MHz characteristic frequencies were defined. The results of melting experiments indicate a thermal destabilization of the oligomers containing the defects. Electrophoretic data and the dielectric results show that the conformations of the nicked and control samples are very similar, while the oligomer with a gap is more compact with a different charge distribution at the ends.

Base sequence criteria and Cartesian coordinates for stable B/Z and B/Z/B junctions in relaxed DNA

It seems increasingly evident that if the Z form of DNA exists in the genome it must exist as short sections of alternating pyrimidine-purine sequences in the midst of very long sections of B-form DNA. We have determined the minimum length of a string of alternating CG base pairs that can go into the Z form in the middle of a long B form. Self-complimentary oligomers of the form T(M)(CG)(N)A(M) were synthesized. The conformation of the resulting duplex was determined in 6M aqueous NaCl solution by Raman scattering. We have found that 12 alternating CG base pairs is the minimum length required to form a stable Z form of DNA inside of a long B form section. Only the 4 center CG base pairs go into the Z form. These 4 CG base pairs in the Z form are flanked on each side by 4 CG base pairs in a non-Z (probably B) form as well as the ..TT.. ..AA.. sequences in the B form. We propose a model of the B/Z junction in which the double helix flips directly from the B form to the Z form so that there are no base pairs in the junction. In this model the B form is nucleated in the AT base pairs on each end and is propagated into the CG base pairs in the center. This model is supported by isotopic H/D exchange experiments that shows that the H/D exchange of the non-Z form CG base pairs is highly retarded and indicates that they remain in the B form. A Thermodynamic analysis of the concentration dependence of the melting point of the duplexes in both low and high salt, supports our model and rules out the possibility of hairpin formation. The enthalpy for the formation of a B/Z junction is determined to be about +16 kcal/junction. A comparison of these results with recent results on B/Z junctions in super-coiled DNA is given. Molecular modeling calculations permit us to obtain values for the coordinates and torsional angles of the oligomers showing both B/Z and B/Z/B junctions. The Cartesian coordinates for these oligomers as well as stereo figures of these models in color are available from the authors.

Flanking sequence effects within the pyrimidine triple-helix motif characterized by affinity cleaving

Nearest neighbor interactions affect the stabilities of triple-helical complexes. Within a pyrimidine triple-helical motif, the relative stabilities of natural base triplets T.AT, C + GC, and G.TA, as well as triplets, D3.TA and D3.CG, containing the nonnatural deoxyribonucleoside 1-(2-deoxy-beta-D-ribofuranosyl)-4-(3-benzamido)phenylimidazole (D3) were characterized by the affinity cleaving method in the context of different flanking triplets (T.AT, T.AT: T.AT, C + GC: C + GC, T.AT: G + GC, C + GC). The to be insensitive to substitutions in either the 3' or 5' directions, while the relative stabilities of triple helices containing C + GC triplets decreased as the number of adjacent C + GC triplets increased. Triple helices incorporating a G.TA interaction were most stable when this triplet was flanked by two T.AT triplets and were adversely affected when a C + GC triplet was placed in the adjacent 5' direction. Similarly, complexes containing D3.TA or D3.CG triplets were most stable when the triplet was flanked by two T.AT triplets but were destabilized when the adjacent 3' neighbor position was occupied with a C + GC triplet. This information regarding sequence composition effects in triple-helix formation establishes a set of guidelines for targeting sequences of double-helical DNA by the pyrimidine triple-helix motif.

Thermodynamic characterization of deoxyribooligonucleotide duplexes containing bulges

Ultraviolet absorption techniques were used to study the thermodynamics of duplex formation for a DNA decamer, d(GCGAAAAGCG).d(CGCTTTTCGC), and a series of related duplexes, each of which contains a bulged base centered in the A.T tract. Thermodynamic parameters were obtained from nonlinear least-squares fits of the melting curves and the concentration dependences of the melting temperatures. Duplexes containing a localized single-base bulge were found to be 3.5-4.6 kcal/mol less stable than the decamer at 37 degrees C. These results indicate that both the identity of the bulged base and the strand in which it is located may influence the amount by which the duplex is destabilized. Bulged bases located in the T-strand, d(CGCTTYTTCGC), in position Y, were observed to be slightly more destabilizing than those located in the A-strand, d(GCGAAXAAGCG), in position X. Bulged purines may be more destabilizing than bulged pyrimidines.

Circular dichroism studies of an oligodeoxyribonucleotide containing a hairpin loop made of a hexaethylene glycol chain: conformation and stability

An oligodeoxyribonucleotide, d(GCTCACAAT-X-ATTGTGAGC), where X represents a hexaethylene glycol chain, was studied using circular dichroism spectroscopy. Its conformation and conformational stability were compared to those of compounds where X was replaced by four thymines and to the duplex of same sequence without loop. The compound with the hexaethylene glycol chain can form a hairpin looped structure as well as a bulged duplex structure. In both cases the duplex region of the oligodeoxyribonucleotide exhibits the same conformation. In similar conditions the oligodeoxyribonucleotide with a four thymines loop forms exclusively a hairpin structure. Comparison between the thermodynamic parameters (delta H, delta S, delta G) associated with the formation of the structure of the three compounds are presented. In the case of the compound with the hexaethylene glycol chain it is shown that the large increase in its melting temperature (by about 35 degrees C in our experimental conditions) when compared to the non looped structure is mainly due to the fact that its melting process is intramolecular (monomolecular) whereas the other one is bimolecular.

The 6-thioguanine/5-methyl-2-pyrimidinone base pair

As part of a program to determine the physical possibility of expanding the number of types of base pairs in DNA, the pairing stabilities of the analog bases 6-thioguanine (GS) and 5-methyl-2-pyrimidinone (TH) in oligodeoxynucleotides were measured. Procedures were developed to synthesize oligodeoxynucleotides with the analog bases. The sequences of the synthesized oligomers were T-C-G-A-C-G-G-X-Y-C-C-G. An enzymatic procedure was developed to measure relative association constants of oligomer pairs with the self complementary reference oligomer, X = A and Y = T, K(T/A) = K. The results were K(C/G) = (5 +/- .5)K, K(TH/GS) = K/(1 +/- .5), K(T/G) = K/(9 +/- 3), K(TH/G) = K/(25 +/- 5), K(C/GS) less than K/30, K(TH/A) less than K/40, K(T/GS) less than K/40, K(C/A) less than K/40. The results with the standard bases are consistent with other methods of measurement. The stability of the base pair GS/TH is approximately the same as the standard base pair A/T.

Thermodynamic studies of base pairing involving 2,6-diaminopurine

The thermal stabilities of oligodeoxyribonucleotide duplexes containing 2,6-diaminopurine (D) matched with each of the four normal DNA bases were determined by optical melting techniques. Comparison of optical melting curves yielded relative stabilities for the D-containing standard base pairs in an otherwise identical base-pair sequence. The D:T pair was found to be more stable than the A:T pair in dC3DG3:dC3TG3, as stable as the A:T in dCT3DT3G:dCA3TA3G, and less stable than the A:T in dCA3DA3G:dCT7G. The order of stabilities for X:Y in the DNA duplex dCA3XA3G:dCT3YT3G is: (A:T) greater than (T:D) congruent to (D:T) greater than or equal to (T:A) greater than (C:D) congruent to (D:A) congruent to (D:G) greater than or equal to (D:C) congruent to (G:D) congruent to (D:D) greater than or equal to (A:D). Implications of these results for design of DNA oligonucleotide probes are discussed.

Ultraviolet-induced dimerization of non-adjacent pyrimidines. A potential mechanism for the targeted -1 frameshift mutation

The DNA photoproduct responsible for the ultraviolet (u.v.)-induced targeted -1 frameshift mutation is unknown. Based on mutagenesis studies by others, we surmised that this lesion might be found in high abundance in single-stranded DNA. u.v. irradiation of the single-stranded alternating copolymer poly[d(G-T)] yielded a photoproduct that was characterized in detail. It consists of a thymine-thymine cyclobutane dimer of predominantly cis-syn configuration occurring between non-adjacent thymidyl residues on the same strand. Its formation is strongly inhibited in double-stranded DNA. A similar u.v. photoproduct was obtained in higher yield from the polypyrimidine alternating copolymer poly[d(C-T)] under conditions in which it is single-stranded. It is proposed that replication across the unrepaired photoproduct: (formula; see text) is the cause of the targeted u.v.-induced -1 frameshift mutation.

Melting of a self-complementary DNA minicircle. Comparison of optical melting theory with exchange broadening of the nuclear magnetic resonance spectrum

Melting curves are calculated for the 16-base-pair duplex DNA sequence 5' GTATCCGTACGGATAC 3' linked on the ends by TTTT single-strand loops. The equilibrium statistical thermodynamic theory of DNA melting is modified to include effects of end-loops on the melting transition. An excellent fit of the experimental melting curve in 0.2 M-NaCl is obtained using two adjustable parameters, one for end-loop formation and the other for formation of the complete 40-base single-strand loop. The best-fit calculated melting curve permits evaluation of these parameters. The free energy to close a TTTT end-loop is 2.12 kcal/mol (1 cal = 4.184 J). A TTTT end-loop or hairpin loop is significantly more stable than an internal loop of comparable size sandwiched between two helical regions, even after allowing for the different stacking contributions. Reasons for this increased stability are presented. The loop free energy of the 40-base single-strand open minicircle is evaluated to be +1.27 kcal/mol, thus favoring the melting of two end-loops into the large open minicircle. The present results are compared with those of others for d(T-A) oligomers. The sequence TTTT forms a more stable end-loop, or hairpin, than TATA by about 2.0 kcal/mol. Theoretical rate constants for the proton-transfer step in the standard hydrogen-exchange model are calculated by extending the theory of diffusion-controlled reactions to take account of the electrostatic potential of the DNA. The predicted ratios of rate constants for different pairs of catalysts exchanging an A.T proton agree satisfactorily with the available experimental data for a 14-base-pair linear duplex, which confirms the diffusion-control of the proton-transfer step. Data presented here for the 16 base-pair duplex of the minicircle are consistent with catalysis-limited exchange in which the proton-transfer step is likewise diffusion-controlled. Under catalysis-limited conditions, the imino proton exchange rates are predicted from the catalytic rate constants, prevailing buffer catalyst concentrations, and the equilibrium constants to form the unstacked open state of optical melting theory. The observed exchange rates of the A.T base-pairs show no sign of the strong predicted end-melting trend, and exceed the predicted values by factors of 10 to 400. Moreover, the succession of "melting" in the nuclear magnetic resonance line-broadening deviates from that predicted by optical melting theory.(ABSTRACT TRUNCATED AT 400 WORDS)

Improved free-energy parameters for predictions of RNA duplex stability

Thermodynamic parameters for prediction of RNA duplex stability are reported. One parameter for duplex initiation and 10 parameters for helix propagation are derived from enthalpy and free-energy changes for helix formation by 45 RNA oligonucleotide duplexes. The oligomer sequences were chosen to maximize reliability of secondary structure predictions. Each of the 10 nearest-neighbor sequences is well-represented among the 45 oligonucleotides, and the sequences were chosen to minimize experimental errors in delta GO at 37 degrees C. These parameters predict melting temperatures of most oligonucleotide duplexes within 5 degrees C. This is about as good as can be expected from the nearest-neighbor model. Free-energy changes for helix propagation at dangling ends, terminal mismatches, and internal G X U mismatches, and free-energy changes for helix initiation at hairpin loops, internal loops, or internal bulges are also tabulated.

Predicting DNA duplex stability from the base sequence

We report the complete thermodynamic library of all 10 Watson-Crick DNA nearest-neighbor interactions. We obtained the relevant thermodynamic data from calorimetric studies on 19 DNA oligomers and 9 DNA polymers. We show how these thermodynamic data can be used to calculate the stability and predict the temperature-dependent behavior of any DNA duplex structure from knowledge of its base sequence. We illustrate our method of calculation by using the nearest-neighbor data to predict transition enthalpies and free energies for a series of DNA oligomers. These predicted values are in excellent agreement with the corresponding values determined experimentally. This agreement demonstrates that a DNA duplex structure thermodynamically can be considered to be the sum of its nearest-neighbor interactions. Armed with this knowledge and the nearest-neighbor thermodynamic data reported here, scientists now will be able to predict the stability (delta G degree) and the melting behavior (delta H degree) of any DNA duplex structure from inspection of its primary sequence. This capability should prove valuable in numerous applications, such as predicting the stability of a probe-gene complex; selecting optimal conditions for a hybridization experiment; deciding on the minimum length of a probe; predicting the influence of a specific transversion or transition on the stability of an affected DNA region; and predicting the relative stabilities of local domains within a DNA duplex.

The stiffness of dsRNA: hydrodynamic studies on fluorescence-labelled RNA segments of bovine rotavirus

The sedimentation coefficients of dsRNA segments of bovine rotavirus were determined in the analytical ultracentrifuge. The eleven segments were separated by preparative gel electrophoresis, and isolated by elution from gel pieces. The RNA was labelled by the intercalating fluorescent dye ethidium bromide at a ratio bound dye per base pair between 0.003 to 0.018. The analytical ultracentrifuge was equipped with a fluorescence recording optics. Sedimentation coefficients could be determined with amounts of RNA as little as 8 ng. All sedimentation coefficients were extrapolated to zero-concentration, zero-dye binding, and zero-impurities from the preparative gel electrophoresis. The hydrodynamic model of flexible cylinders was applied for the interpretation of the sedimentation coefficients. All dsRNA segments of rotavirus (663-3409 base pairs) and the dsRNA5 of cucumber mosaic virus (335 base pairs) fit the model of a "worm-like" or flexible cylinder with a persistence length of 1125 A and a hydrated diameter of 30 A. The results are compared with data from the literature on the persistence lengths of the B- and Z-forms of dsDNA and of viroids.

Temperature jump relaxation studies on the interactions between transfer RNAs with complementary anticodons. The effect of modified bases adjacent to the anticodon triplet

We have used the temperature-jump relaxation technique to determine the kinetic and thermodynamic parameters for the association between the following tRNAs pairs having complementary anticodons: tRNA(Ser) with tRNA(Gly), tRNA(Cys) with tRNA(Ala) and tRNA(Trp) with tRNA(Pro). The anticodon sequence of E. coli tRNA(Ser), GGA, is complementary to the U*CC anticodon of E. coli tRNA(Gly(2] (where U* is a still unknown modified uridine base) and A37 is not modified in none of these two tRNAs. E. coli tRNA(Ala) has a VGC anticodon (V is 5-oxyacetic acid uridine) while tRNA(Cys) has the complementary GCA anticodon with a modified adenine on the 3' side, namely 2-methylthio N6-isopentenyl adenine (mS2i6A37) in E. Coli tRNA(Cys) and N6-isopentenyl adenine (i6A37) in yeast tRNA(Cys). The brewer yeast tRNA(Trp) (anticodon CmCA) differs from the wild type E. coli tRNA(Trp) (anticodon CCA) in several positions of the nucleotide sequence. Nevertheless, in the anticodon loop, only two interesting differences are present: A37 is not modified while C34 at the first anticodon position is modified into a ribose 2'-O methyl derivative (Cm). The corresponding complementary tRNA is E.coli tRNA(Pro) with the VGG anticodon. Our results indicate a dominant effect of the nature and sequence of the anticodon bases and their nearest neighbor in the anticodon loop (particularly at position 37 on the 3' side); no detectable influence of modifications in the other tRNA stems has been detected. We found a strong stabilizing effect of the methylthio group on i6A37 as compared to isopentenyl modification of the same residue. We have not been able so far to assess the effect of isopentenyl modification alone in comparison to unmodified A37. The results obtained with the complex yeast tRNA(Trp)-E.coli tRNA(Pro) also suggest that a modification of C34 to Cm34 does not significantly increase the stability of tRNA(Trp) association with its complementary anticodon in tRNA(Pro). The observations are discussed in the light of inter- and intra-strand stacking interactions among the anticodon triplets and with the purine base adjacent to them, and of possible biological implications.

Consequences of methylation on the amino group of adenine. A proton two-dimensional NMR study of d(GGATATCC) and d(GGm6ATATCC)

A two-dimensional 500-MHz 1H-NMR study of two oligonucleotides, d(GGATATCC) and d(GGm6ATATCC), is presented in which we have investigated the effects of adenine methylation. The two-dimensional nuclear Overhauser spectra (NOESY) show that both oligonucleotides adopt a normal right-handed B-type helix and one-dimensional nuclear Overhauser enhancement (NOE) studies demonstrate that any difference in conformation must be small. However methylation drastically slows down the helix in equilibrium coil exchange which becomes slow on a proton NMR time scale. While d(GGATATCC) fits a two-site exchange model, d(GGm6ATATCC) does not and we invoke the presence of a third species which may be an intermediate in helix formation. NMR and ultraviolet spectroscopy show that methylation destabilizes the helix, measured by the melting temperature and enthalpy of dissociation.

Base-base mismatches. Thermodynamics of double helix formation for dCA3XA3G + dCT3YT3G (X, Y = A,C,G,T)

Thermodynamic parameters for double strand formation have been measured for the sixteen double helices of the sequence dCA3XA3G.dCT3YT3G, with each of the bases A, C, G and T at the positions labelled X and Y. The results are analyzed in terms of nearest-neighbors and are compared with thermodynamic parameters for RNA secondary structure. At room temperature the sequence (Formula: see text) is more stable than (Formula: see text) and is similar in stability to (Formula: see text) and (Formula: see text) are least stable. At higher temperatures the sequences containing a G.C base pair become more stable than those containing only A.T. All molecules containing mismatches are destabilized with respect to those with only Watson-Crick pairing, but there is a wide range of destabilization. At room temperature the most stable mismatches are those containing guanine (G.T, G.G, G.A); the least stable contain cytosine (C.A, C.C). At higher temperatures pyrimidine-pyrimidine mismatches become the least stable.

Deoxydodecanucleotide heteroduplex d(TTTTATAATAAA). d(TTTATTATAAAA) containing the promoter Pribnow sequence TATAAT. I. Double-helix stability by UV spectrophotometry and calorimetry

Thermally induced structural transition in the d(TTTTATAATAAA) d(TTTATTATAAAA) heteroduplex is characterized by UV-spectroscopy and differential scanning calorimetry. At low salt (less than 0.1 M) the occurrence of a cooperative transition in the lower temperature range, followed by a broad transition connected with small increase in absorbance is observed. At high salt (greater than or equal to 0.2 M) a single, monophasic transition appears. Linear dependence of the latter on log of salt concentration (dTm:dlogM = 14.2 degrees C) and of 1/Tm on log of oligomer concentration [derived therefrom delta H (v.H.) = 77.1 kcal/mole (duplex)] allows relating it to the melting of the heteroduplex helix. The non-cooperative transition, independent of oligomer concentration and similar to that of the single chain, was attributed to melting of short hairpin helices upon heteroduplex dissociation. Calorimetric enthalpy: 75.6 kcal/mole (duplex) proved significantly lower than predicted from known calorimetric data for poly[d(AT)] and poly d(A) X poly d(T).

Thermodynamic studies of RNA stability

Enthalpies and entropies of helix stabilization due to addition of 3' terminal unpaired nucleotides to a CCGG or GGCC core double helix are derived from UV melting studies. The results suggest stacking provides a significant fraction of the free energy of a terminal base pair. The effects of temperature, aggregation, and ionic strength on the determination of thermodynamic parameters are considered. Helix propagation parameters are revised and extended based on recent additions to the data set.

Escherichia coli polymerase I can use O2-methyldeoxythymidine or O4-methyldeoxythymidine in place of deoxythymidine in primed poly(dA-dT).poly(dA-dT) synthesis

O2-and O4-alkyldeoxythymidine are among the four O-alkyl base-modified derivatives produced by the reaction of N-nitroso alkylating agents with nucleic acids in vitro and in vivo. We find that both O2- and O4-methyl-dTTP can substitute for dTTP in alternating poly(dA-dT)-primed DNA synthesis. Up to 22% of the pyrimidines in the newly synthesized polymer were found by HPLC analysis to be O-methyldeoxythymidine. Little polymer synthesis was observed in the absence of dTTP. However, the O-methyl-dTTPs did not inhibit polymerization of dATP and dTTP. Polymers containing O2- or O4-methyldeoxythymidine were obtained in good yield, retaining the secondary structure of alternating poly(dA-dT). This was shown by the data for thermal transition under different conditions. In contrast, poly(dA-dT).poly(dA-dT) methylated or ethylated to less than 4% total modification by alkylnitrosoureas had a distinctly less stable structure. Neither O2- nor O4-methyldeoxythymidine can form more than one hydrogen bond with adenosine. The unchanged secondary structure of polymers containing these modified thymidines indicates that stacking interactions must play a major role in helix stabilization. O-Alkyldeoxythymidine may be formed by N-nitroso carcinogens that react intracellularly. We have shown that the triphosphates can be utilized by Escherichia coli DNA polymerase I as dTTP. The incorporated O4-methyl-dT causes misincorporation of G, both in transcription and synthesis. When O2-methyl-dT is present, less, but definite, misincorporation results.

Spectroscopic studies of (m5dC-dG)3: thermal stability of B- and Z-forms

The hexanucleoside pentaphosphate d(m5CpGpm5CpGpm5CpG) has been studied in solution by ultra-violet absorption, circular dichroism and 31P nuclear magnetic resonance under various experimental conditions. In 0.2 M NaClO4 at low temperature, an hexamer duplex is formed which has a B or B-like conformation. As the salt concentration is increased, a transition from a B-form to the Z-form occurs and is complete in 3 M NaClO4. In 3 M NaClO4, the behavior of the Z double helix is complex as a function of temperature. The variation of the circular dichroism at 295 nm is biphasic. A first transition occurs over a large range of temperature and corresponds to a conformational change due to a non-cooperative intramolecular process. Ultra-violet absorption and 31P nuclear magnetic resonance show that the new conformation arising from a distortion of the backbone is not similar to that observed in low salt conditions (B-form). At high hexanucleotide concentration, aggregates are formed. The second transition is cooperative and corresponds to the melting of a double stranded helix into single strands.