Thermodynamic parameters based on a nearest-neighbor model for DNA sequences with a single-bulge loop

All 64 possible thermodynamic parameters for a single-bulge loop in the middle of a sequence were derived from optical melting studies. The relative stability of a single bulge depended on both the type of bulged base and its flanking base pairs. The contribution of the single bulge to helix stability ranged from 3.69 kcal/mol for a TAT bulge to -1.05 kcal/mol for an ACC bulge. Thermodynamics for 10 sequences with a GTG bulge were determined to test the applicability of the nearest-neighbor model to a single-bulge loop. Thermodynamic parameters for the GTG bulge and Watson-Crick base pairs predict, DeltaH degrees, DeltaS degrees, and T(M)(50 microM) values with average deviations of 3.0%, 4.3%, 4.7%, and 0.9 degrees C, respectively. The prediction accuracy was within the limits of what can be expected for a nearest-neighbor model. This certified that the thermodynamics for single-bulge loops can be estimated adequately using a nearest-neighbor model.     

Stabilization factors affecting duplex formation of peptide nucleic acid with DNA.

Peptide nucleic acid (PNA) is an oligonucleotide analogue in which the sugar-phosphate backbone is replaced by an N-(2-aminoethyl)glycine unit to which the nucleobases are attached. We investigated the thermodynamic behavior of PNA/DNA hybrid duplexes with identical nearest neighbors but with different sequences and chain lengths (5, 6, 7, 8, 10, 12, and 16 mers) to reveal whether the nearest-neighbor model is valid for the PNA/DNA duplex stability. CD spectra of 6, 7, and 8 mer PNA/DNA duplexes showed similar signal, while 10, 12, and 16 mer duplexes did not. The average difference in Delta G degrees (37) for short PNA/DNA duplexes with identical nearest-neighbor pairs was only 3.5%, whereas that of longer duplexes (10, 12, and 16 mers) was 16.4%. Therefore, the nearest-neighbor model seems to be useful at least for the short PNA/DNA duplexes. Thermodynamics of PNA/DNA duplexes containing 1--3 bulge residues were also studied. While the stability of the 12 mer DNA/DNA duplex decreased as the number of bulge bases increases, the number of bulge bases in PNA/DNA unchanged the duplex stability. Thus, the influence of bulge insertion in the PNA/DNA duplexes is different from that of a DNA/DNA duplex. This might be due to the different base geometry in a helix which may potentially make hydrogen bonds in a base pair and stacking interaction unfavorable compared with DNA/DNA duplexes.     

Kinetics and thermodynamics of double-helix formation of self-complementary deoxyribooctanucleotides d(TCTATAGA) and d(TAGATCTA).

Double-helix formation of self-complementary deoxyribooctanucleotides, d(TCTATAGA) and d(TAGATCTA), with identical nearest neighbor base pairs has been studied by means of UV melting and temperature-jump techniques. The self-complementary duplexes of both octanucleotides with identical nearest neighbors had similar stabilities: The stabilization energies of the octanucleotides at 25 degrees C were 5.8 kcal mol-1 for d(TCTATAGA) and 6.7 kcal mol-1 for d(TAGATCTA). On the kinetic curve the melting reactions finished within 20 ms for d(TCTATAGA) and 40 ms for d(TAGATCTA) at 20 degrees C. For both octanucleotides the rate constants of dissociation increased and the rate constants of association decreased with increasing temperature.     

The thermal stability of DNA fragments with tandem mismatches at a d(CXYG).d(CY'X'G) site.

Temperature-Gradient Gel Electrophoresis (TGGE) was employed to determine the thermal stabilities of 28 DNA fragments, 373 bp long, with two adjacent mismatched base pairs, and eight DNAs with Watson-Crick base pairs at the same positions. Heteroduplex DNAs containing two adjacent mismatches were formed by melting and reannealing pairs of homologous 373 bp DNA fragments differing by two adjacent base pairs. Product DNAs were separated based on their thermal stability by parallel and perpendicular TGGE. The polyacrylamide gel contained 3.36 M urea and 19.2 % formamide to lower the DNA melting temperatures. The order of stability was determined in the sequence context d(CXYG).d(CY'X'G) where X.X' and Y.Y" represent the mismatched or Watson-Crick base pairs. The identity of the mismatched bases and their stacking interactions influence DNA stability. Mobility transition melting temperatures (T u) of the DNAs with adjacent mismatches were 1.0-3.6 degrees C (+/-0.2 degree C) lower than the homoduplex DNA with the d(CCAG).d(CTGG) sequence. Two adjacent G.A pairs, d(CGAG).d(CGAG), created a more stable DNA than DNAs with Watson-Crick A.T pairs at the same sites. The d(GA).d(GA) sequence is estimated to be 0.4 (+/-30%) kcal/mol more stable in free energy than d(AA).d(TT) base pairs. This result confirms the unusual stability of the d(GA).d(GA) sequence previously observed in DNA oligomers. All other DNAs with adjacent mismatched base pairs were less stable than Watson-Crick homoduplex DNAs. Their relative stabilities followed an order expected from previous results on single mismatches. Two homoduplex DNAs with identical nearest neighbor sequences but different next-nearest neighbor sequences had a small but reproducible difference in T u value. This result indicates that sequence dependent next neighbor stacking interactions influence DNA stability.     

Improved free energies for G.C base-pairs

Thermodynamic parameters of helix formation are reported for seven oligoribonucleotides containing only G.C pairs. These data are used with the nearest-neighbor model to calculate enthalpies and free energies of base-pair formation for G.C pairs. For helix initiation, the free energy change at 37 degrees C, delta G(0)37, is +3.9 kcal/mol; for helix propagation, the delta G(0)37 values are -2.3, -3.2 and -3.3 kcal/mol for C-G, G-G and G-C neighbors, respectively.     

Studies of DNA dumbbells. I. Melting curves of 17 DNA dumbbells with different duplex stem sequences linked by T4 endloops: evaluation of the nearest-neighbor stacking interactions in DNA.

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 sequence (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 degrees 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, delta Gi, and reveal several interesting features: (1) Titratable changes in the values of delta Gi with changing salt environment are observed. In all salts the most stable unique combination is delta G4 = (delta GGpC+delta GCpG)/2, and the least stable is the GpG/CpC stack, delta G2 = delta GGpG/CpC. (2) The chi 2 values of the fits of the evaluated delta Gi's to experimental data increased with decreasing [Na+], suggesting that significant interactions beyond nearest neighbors become more pronounced, particularly at 25 nM Na+.(ABSTRACT TRUNCATED AT 400 WORDS)     

Solution structure and thermal stability of ribosomal protein L30e from hyperthermophilic archaeon Thermococcus celer

To understand the structural basis of thermostability, we have determined the solution structure of a thermophilic ribosomal protein L30e from Thermococcus celer by NMR spectroscopy. The conformational stability of T. celer L30e was measured by guanidine and thermal-induced denaturation, and compared with that obtained for yeast L30e, a mesophilic homolog. The melting temperature of T. celer L30e was 94 degrees C, whereas the yeast protein denatured irreversibly at temperatures >45 degrees C. The two homologous proteins also differ greatly in their stability at 25 degrees C: the free energy of unfolding was 45 kJ/mole for T. celer L30e and 14 kJ/mole for the yeast homolog. The solution structure of T. celer L30e was compared with that of the yeast homolog. Although the two homologous proteins do not differ significantly in their number of hydrogen bonds and the amount of solvent accessible surface area buried with folding, the thermophilic T. celer L30e was found to have more long-range ion pairs, more proline residues in loops, and better helix capping residues in helix-1 and helix-4. A K9A variant of T. celer L30e was created by site-directed mutagenesis to examine the role of electrostatic interactions on protein stability. Although the melting temperatures of the K9A variant is approximately 8 degrees C lower than that of the wild-type L30e, their difference in T(m) is narrowed to approximately 4.2 degrees C at 0.5 M NaCl. This salt-dependency of melting temperatures strongly suggests that electrostatic interactions contribute to the thermostability of T. celer L30e.     

Non-nearest-neighbor dependence of the stability for RNA bulge loops based on the complete set of group I single-nucleotide bulge loops

Fifty-nine RNA duplexes containing single-nucleotide bulge loops were optically melted in 1 M NaCl, and the thermodynamic parameters DeltaH degrees, DeltaS degrees, DeltaG 37 degrees, and TM for each sequence were determined. Sequences from this study were combined with sequences from previous studies [Longfellow, C. E., et al. (1990) Biochemistry 29, 278-285; Znosko, B. M., et al. (2002) Biochemistry 41, 10406-10417], thus examining all possible group I single-nucleotide bulge loop and nearest-neighbor sequence combinations. The free energy increments at 37 degrees C for the introduction of a group I single-nucleotide bulge loop range between 1.3 and 5.2 kcal/mol. The combined data were used to develop a model for predicting the free energy of a RNA duplex containing a single-nucleotide bulge. For bulge loops with adjacent Watson-Crick base pairs, neither the identity of the bulge nor the nearest-neighbor base pairs had an effect on the influence of the bulge loop on duplex stability. The proposed model for prediction of the stability of a duplex containing a bulged nucleotide was primarily affected by non-nearest-neighbor interactions. The destabilization of the duplex by the bulge was related to the stability of the stems adjacent to the bulge. Specifically, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. The stability of a duplex containing a bulged nucleotide adjacent to a wobble base pair also was primarily affected by non-nearest-neighbor interactions. Again, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. However, when one or both of the bulge nearest neighbors was a wobble base pair, the free energy increment for insertion of a bulge loop is dependent upon the position and orientation of the wobble base pair relative the bulged nucleotide. Bulge sequences of the type ((5'UBX)(3'GY)), ((5'GBG)(3'UU)) and ((5'UBU)(3'GG)) are less destabilizing by 0.6 kcal/mol, and bulge sequences of the type ((5'GBX)(3'UY)) and ((5'XBU)(3'YG)) are more destabilizing by 0.4 kcal/mol than bulge loops adjacent to Watson-Crick base pairs.     

An estimate of the nearest neighbor base-pair content of 5S RNA using CD and absorption spectroscopy.

We analyzed the CD and uv absorption spectra of 5S RNA from Escherichia coli using the method developed in the preceding paper. The analysis of spectra of 5S RNA at 20 degrees C in 0.1M NaClO4, 2.5 mM Na+ (phosphate), pH 7.0, and 0.5 mM MgSO4 gave 7 +/- 3.6 A.U base pairs, 25 +/- 3.6 G.C base pairs, and 7.5 +/- 3.6 G.U base pairs. Estimates of nearest neighbor base pairs were more consistent with the Pieler-Erdmann and the Gewirth-Moore secondary structure models than with the Fox-Woese or the Burns-Luoma-Marshall models. We also examined the structure of 5S RNA as a function of temperature. The melting profile exhibited two transitions--one at about 35 degrees C and one above 50 degrees C. Our spectral data showed that helices I and II were stable during the first transition, and agreed with other data that helix III was the most likely helix to have melted. The results from this in-depth study of 5S RNA indicate that our method of analysis should be useful for studying the secondary structures of other small, unmodified RNAs.     

N(4)C-ethyl-N(4)C cross-linked DNA: synthesis and characterization of duplexes with interstrand cross-links of different orientations.

The preparation and physical properties of short DNA duplexes that contain a N(4)C-ethyl-N(4)C interstrand cross-link are described. Duplexes that contain an interstrand cross-link between mismatched C-C residues and duplexes in which the C residues of a -CG- or -GC- step are linked to give "staggered" interstrand cross-links were prepared using a novel N(4)C-ethyl-N(4)C phosphoramidite reagent. Duplexes with the C-C mismatch cross-link have UV thermal transition temperatures that are 25 degrees C higher than the melting temperatures of control duplexes in which the cross-link is replaced with a G-C base pair. It appears that this cross-link stabilizes adjacent base pairs and does not perturb the structure of the helix, a conclusion that is supported by the CD spectrum of this duplex and by molecular models. An even higher level of stabilization, 49 degrees C, is seen with the duplex that contains a -CG- staggered cross-link. Molecular models suggest that this cross-link may induce propeller twisting in the cross-linked base pairs, and the CD spectrum of this duplex exhibits an unusual negative band at 298 nm, although the remainder of the spectrum is similar to that of B-form DNA. Mismatched C-C or -CG- staggered cross-linked duplexes that have complementary overhanging ends can undergo self-ligation catalyzed by T4 DNA ligase. Analysis of the ligated oligomers by nondenaturing polyacrylamide gel electrophoresis shows that the resulting oligomers migrate in a manner similar to that of a mixture of non-cross-linked control oligomers and suggests that these cross-links do not result in significant bending of the helix. However, the orientation of the staggered cross-link can have a significant effect on the structure and stability of the cross-linked duplex. Thus, the thermal stability of the duplex that contains a -GC- staggered cross-link is 10 degrees C lower than the melting temperature of the control, non-cross-linked duplex. Unlike the -CG- staggered cross-link, in which the cross-linked base pairs can still maintain hydrogen bond contacts, molecular models suggest that formation of the -GC- staggered cross-link disrupts hydrogen bonding and may also perturb adjacent base pairs leading to an overall reduction in helix stability. Duplexes with specifically positioned and oriented cross-links can be used as substrates to study DNA repair mechanisms.     

A set of nearest neighbor parameters for predicting the enthalpy change of RNA secondary structure formation

A complete set of nearest neighbor parameters to predict the enthalpy change of RNA secondary structure formation was derived. These parameters can be used with available free energy nearest neighbor parameters to extend the secondary structure prediction of RNA sequences to temperatures other than 37°C. The parameters were tested by predicting the secondary structures of sequences with known secondary structure that are from organisms with known optimal growth temperatures. Compared with the previous set of enthalpy nearest neighbor parameters, the sensitivity of base pair prediction improved from 65.2 to 68.9% at optimal growth temperatures ranging from 10 to 60°C. Base pair probabilities were predicted with a partition function and the positive predictive value of structure prediction is 90.4% when considering the base pairs in the lowest free energy structure with pairing probability of 0.99 or above. Moreover, a strong correlation is found between the predicted melting temperatures of RNA sequences and the optimal growth temperatures of the host organism. This indicates that organisms that live at higher temperatures have evolved RNA sequences with higher melting temperatures.     

Analysis of melting transitions of the DNA hairpins formed from the oligomer sequences d[GGATAC(X)4GTATCC] (X = A, T, G, C)

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 degrees C for the T4 loop hairpin to 55.3 degrees 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 formation 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.(ABSTRACT TRUNCATED AT 400 WORDS)     

Purine-purine mismatches in RNA helices: evidence for protonated G.A pairs and next-nearest neighbor effects.

Thermodynamic parameters are presented for 12 different RNA duplexes containing A.A, A.G, G.A and G.G mismatches flanked by C-G base pairs. UV melting studies were conducted under three different buffer conditions in order to evaluate the effects of salt concentration and pH on the stability of each mismatch-containing duplex. The main findings are: (i) the mismatches have a wide range of effects on duplex stability, decreasing delta G degrees 37 of denaturation by approximately 0-7 kcal/mol; (ii) the nearest-neighbor assumption commonly used to calculate helix stability breaks down for G.A mismatches; and (iii) G.A mismatches separated by 2 bp form a protonated structure.     

NMR evidence of the stabilisation by the carcinogen N-2-acetylaminofluorene of a frameshift mutagenesis intermediate.

Two heteroduplexes d(C1A2C3T4C5G6C7A8C9A10C11)-d (G12T13G14T15G16G17A18G19T20G21) containing a bulged guanine either unmodified or modified with the carcinogen N-2-acetylaminofluorene (AAF) have been studied by nuclear magnetic resonance (NMR) as models of slipped mutagenic intermediates (SMI). Conformational equilibria are observed in both the unmodified and the AAF-modified heteroduplexes. The major conformation of the unmodified duplex is one where the extra guanine is stacked in the helix and the major conformation of the AAF-modified heteroduplex is one where the AAF is external to the helix. Unusual sugar proton chemical shifts of C5- and G6-AAF indicate that the AAF ring is pointing out in the 5' direction. A strong increase in the modified heteroduplex melting temperature (+15 degrees C) is observed. Moreover, in contrast to the unmodified heteroduplex, which shows extensive melting in the vicinity of the bulged guanine, the base pairs around the bulge in the AAF-modified heteroduplex remain paired at temperatures up to 30 degrees C. This exceptional stability of the site around the bulged modified guanine is suggested to be responsible for the high rate of -1 mutation induced by AAF at repetitive sequences.     

Derivation of nearest-neighbor properties from data on nucleic acid oligomers. II. Thermodynamic parameters of DNA.RNA hybrids and DNA duplexes

Using nearest-neighbor models consisting of independent short sequence combinations of nearest neighbors (ISS models), values of thermodynamic parameters for sets of independent sequences are derived from published oligomer data for DNA.RNA hybrids [N. Sugimoto, S. Nakano, M. Katoh, A. Matsumura, H. Nakamuta, T. Ohmichi, M. Yoneyama, and M. Sasaki (1995) Biochemistry, Vol. 34, pp. 11211-11216] and dsDNA duplexes [J. SantaLucia, Jr., H. T. Allawi, and P. A. Seneviratne (1996) Biochemistry, Vol. 35, pp. 3555-3562]. The results are compared with those from models that assign values of thermodynamic parameters to individual nearest neighbors (INN models). Differences in the use of ISS and INN models are also illustrated in an appendix, which shows examples of analyses for values of a fictitious nearest-neighbor property. INN models that include an initiation parameter contain an implicit assumption that combinations of end neighbors have the same value of a property. It is found that combinations of end neighbors (e.g., base pairs neighboring solvent) in oligomers can have significant and different apparent values of thermodynamic properties, so that the assumption inherent in INN models is not always correct. Even though ISS models do not allow the assignment of values to individual nearest neighbors, except for the like neighbors [such as d(AA)/r(UU), etc., for hybrids and d(AA)/d(TT) and d(GG)/d(CC) for DNA duplexes], they do provide physically meaningful values for the like neighbors, for sequence combinations, and for specified combinations of end neighbors.     

A parallel stranded linear DNA duplex incorporating dG.dC base pairs

DNA oligonucleotides with appropriately designed complementary sequences can form a duplex in which the two strands are paired in a parallel orientation and not in the conventional antiparallel double helix of B-DNA. All parallel stranded (ps) molecules reported to date have consisted exclusively of dA.dT base pairs. We have substituted four dA.dT base pairs of a 25-nt parallel stranded linear duplex (ps-D1.D2) with dG.dC base pairs. The two strands still adopt a duplex structure with the characteristic spectroscopic properties of the ps conformation but with a reduced thermodynamic stability. Thus, the melting temperature of the ps duplex with four dG.dC base pairs (ps-D5.D6) is 10-16 degrees C lower and the van't Hoff enthalpy difference delta HvH for the helix-coil transition is reduced by 20% (in NaCl) and 10% (in MgCl2) compared to that of ps-D1.D2. Based on energy minimizations of a ps-[d(T5GA5).d(A5CT5)] duplex using force field calculations we propose a model for the conformation of a trans dG.dC base pair in a ps helix.     

Base-pair neutral homozygotes can be discriminated by calibrated high-resolution melting of small amplicons

Genotyping by high-resolution melting analysis of small amplicons is homogeneous and simple. However, this approach can be limited by physical and chemical components of the system that contribute to intersample melting variation. It is challenging for this method to distinguish homozygous G::C from C::G or A::T from T::A base-pair neutral variants, which comprise approximately 16% of all human single nucleotide polymorphisms (SNPs). We used internal oligonucleotide calibrators and custom analysis software to improve small amplicon (42-86 bp) genotyping on the LightScanner. Three G/C (PAH c.1155C>G, CHK2 c.1-3850G>C and candidate gene BX647987 c.261+22,290C>G) and three T/A (CPS1 c.3405-29A>T, OTC c.299-8T>A and MSH2 c.1511-9A>T) human single nucleotide variants were analyzed. Calibration improved homozygote genotyping accuracy from 91.7 to 99.7% across 1105 amplicons from 141 samples for five of the six targets. The average T(m) standard deviations of these targets decreased from 0.067 degrees C before calibration to 0.022 degrees C after calibration. We were unable to generate a small amplicon that could discriminate the BX647987 c.261+22,290C>G (rs1869458) SNP, despite reducing standard deviations from 0.086 degrees C to 0.032 degrees C. Two of the sites contained symmetric nearest neighbors adjacent to the SNPs. Unexpectedly, we were able to distinguish these homozygotes by T(m) even though current nearest neighbor models predict that the two homozygous alleles would be identical.     

Stability of 3' double nucleotide overhangs that model the 3' ends of siRNA

Thermodynamic parameters are reported for duplex formation in 1 M NaCl for 16 RNA sequences, each containing a core tetramer duplex, GGCC, and a 3' overhang consisting of two bases. The results indicate additional double-helical stability is conferred by the double 3' terminal overhang relative to the single 3' terminal overhang. A nearest-neighbor analysis of the data indicates that the free energy contribution at 37 degrees C of the second base in the double 3' terminal overhang varies from 0 to 0.7 kcal/mol. The second base in the 3' double overhang can contribute nearly the same stability to a duplex as a base pair or a 3' dangling overhang. Stability contribution of a dangling base, two nucleotides removed from the 3' end of a duplex, is dependent upon both the identity of the base as well as that of the dangling base that it neighbors. A second dangling base only increases the stability of the duplex when it is neighboring a 3' purine dangling nucleotide. Furthermore, a second dangling pyrimidine provides a greater contribution to duplex stability than a purine. A nearest-neighbor model was developed to predict the influence of 3' double overhang on the stability of duplex formation. The model improves the prediction of free energy and melting temperature when tested against six sequences with different core duplexes.     

Comparison of thermodynamic stabilities between PNA (peptide nucleic acid)/DNA hybrid duplexes and DNA/DNA duplexes.

We have examined quantitatively stabilities of PNA/DNA hybrid duplexes with identical nearest-neighbor base pairs and compared stabilities between PNA/DNA and DNA/DNA. The average difference of stabilization energy of the short PNA/DNA was 0.9 kcal mol(-1), which suggests that the stability of the hybrids with identical nearest-neighbor base pairs can be predicted with the nearest-neighbor model as well as those of nucleic acid duplexes.     

Derivation of nearest-neighbor properties from data on nucleic acid oligomers. I. Simple sets of independent sequences and the influence of absent nearest neighbors

The constraints on combinations of nearest neighbors in nucleic acid sequences and the numbers of independent sequences needed to describe nearest-neighbor properties of oligomers and polymers are derived and summarized. It has been pointed out in previous work [D. M. Gray and I. Tinoco, Jr. (1970) Biopolymers, Vol. 9, pp. 223-244; R. F. Goldstein and A. S. Benight (1992) Biopolymers, Vol. 32, pp. 1679-1693] that these constraints restrict the information available from measurements of properties of sequence combinations. The emphasis in this paper is on the properties of oligomer sequences that vary in length, where each nucleotide or base pair at the end of the sequence makes a significant contribution to the measured property by interacting with its boundary of fixed sequence or solvent. In such cases it is not be possible to determine values of properties of individual nearest neighbors, except for the like neighbors [e.g., d(A-A), d(G-G), d(T-T), and d(C-C) nucleotide neighbors in single-stranded DNA or d(A-A)/d(T-T) and d(G-G)/d(C-C) base pair neighbors in double-stranded DNA], solely from measurements of properties of different sequences. Even values for properties of the like neighbors cannot be determined from such oligomeric sequences if the sequences are all of the same length. Nearest-neighbor properties of oligomer sequences that vary in length can be summarized in terms of the values for independent sets of sequences that are nearest neighbors and monomers all with boundaries of the fixed sequence or solvent. Straightforward combinations of the values for the independent sequences will give the values of the property for any dependent sequence, without explicit knowledge of the individual nearest-neighbor values. These considerations have important consequences for the derivation of widely used thermodynamic parameters, as discussed in the following paper.