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 new approach to the study of sequence-dependent properties of polynucleotides

A sequence-dependent properly of any long polynucleotide can be related to the properties of a limited number of other polynucleotides. The properties of the constituent sequence subunits need not be individually known. Only the properties of a set of linearly independent sequence combinations in polymers need be known. The number of independent sequences depends on the number of neighboring bases (or base pairs) contributing to the property and on the number of different bases (or base pairs) allowed. General formulae are derived. If only nearest neighbors contribute to a property, there are 13 linearly independent single-strand sequences with four different bases and there are 8 linearly independent double-strand sequences with two different base pairs. The study of independent sequences in polymers should be especially useful with double-stranded polynucleotides. It should be possible by this polymer approach to estimate nearest-neighbor frequencies from the circular dichroism of double strands and also to investigate the relation between solution conformation and double-strand sequence.     

Dimer of 2,7-diamino-1,8-naphthyridine for the detection of mismatches formed by pyrimidine nucleotide bases

Discrimination of base mismatches from normal Watson-Crick base pairs in duplex DNA constitutes a key approach to the detection of single nucleotide polymorphisms (SNPs). We have developed a sensor for a surface plasmon resonance (SPR) assay system to detect G-G, A-A, and C-C mismatch duplexes by employing a surface upon which mismatch-binding ligands (MBLs) are immobilized. We synthesized a new MBL consisting of 2,7-diamino-1,8-naphthyridine (damND) and immobilized it onto a CM5 sensor chip to carry out the SPR assay of DNA duplexes containing a single-base mismatch. The SPR sensor with damND revealed strong responses to all C-C mismatches, and sequence-dependent C-T and T-T mismatches. Compared to ND- and naphthyridine-azaquinolone hybrid (NA)-immobilized sensor surfaces, with affinity to mismatches composed of purine nucleotide bases, the damND-immobilized surface was useful for the detection of the mismatches composed of pyrimidine nucleotide bases.     

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)     

Mispair specificity of methyl-directed DNA mismatch correction in vitro

To evaluate the substrate specificity of methyl-directed mismatch repair in Escherichia coli extracts, we have constructed a set of DNA heteroduplexes, each of which contains one of the eight possible single base pair mismatches and a single hemimethylated d(GATC) site. Although all eight mismatches were located at the same position within heteroduplex molecules and were embedded within the same sequence environment, they were not corrected with equal efficiencies in vitro. G-T was corrected most efficiently, with A-C, C-T, A-A, T-T, and G-G being repaired at rates 40-80% of that of the G-T mispair. Correction of each of these six mispairs occurred in a methyl-directed manner in a reaction requiring mutH, mutL, and mutS gene products. C-C and A-G mismatches showed different behavior. C-C was an extremely poor substrate for correction while repair of A-G was anomalous. Although A-G was corrected to A-T by the mutHLS-dependent, methyl-directed pathway, repair of A-G to C-G occurred largely by a pathway that is independent of the methylation state of the heteroduplex and which does not require mutH, mutL, or mutS gene products. Similar results were obtained with a second A-G mismatch in a different sequence environment suggesting that a novel pathway may exist for processing A-G mispairs to C-G base pairs. As judged by DNase I footprint analysis, MutS protein is capable of recognizing each of the eight possible base-base mismatches. Use of this method to estimate the apparent affinity of MutS protein for each of the mispairs revealed a rough correlation between MutS affinity and efficiency of correction by the methyl-directed pathway. However, the A-C mismatch was an exception in this respect indicating that interactions other than mismatch recognition may contribute to the efficiency of repair.     

Polymer-supported RNA synthesis and its application to test the nearest-neighbor model for duplex stability

A solid-phase method using a phosphoramidite approach is described for synthesis of oligoribonucleotides. The method was used to synthesize pairs of oligomers with identical nearest neighbors but different sequences. Comparison of thermodynamic parameters for these pairs provides a test of the nearest-neighbor hypothesis for prediction of helix stability. In general, pairs of sequences with identical nearest neighbors have enthalpy and entropy changes for helix formation that differ by 8% on average, delta Go37 that differ by 6% on average, and melting temperatures within 0-5 degrees C of each other. These limits are typical of the accuracy that should be expected from nearest-neighbor predictions of RNA helix stability. UCAUGA and UGAUCA have the same nearest neighbors but melting temperatures that differ by 7 degrees C. This suggests some sequences will not be approximated well by the nearest-neighbor model.     

A comparison of the binding to polynucleotides of complementary and noncomplementary oligonucleotides.

The binding of oligonucleotides to synthetic polynucleotides has been studied as a control for investigations of the binding of oligonucleotides to natural RNA molecules. Only combinations that involved A-U, G-C, and G-G pairs were found to be significantly stable under the experimental conditions used here. The stability of the oligomerpolymer pairing increased with the length of the region paired and with its G + C content. Further, some different sequence isomers of the same G + C content exhibited quite different binding constants. This variability is consistent with certain sequence differences in the double-strand stacking interactions stabilizing the oligomer-polymer association. Oligomer binding was also shown to depend upon the identity of the polymer residues neighboring the binding site, indicating the effect upon oligomer binding of small changes in the single-strand conformation of the binding site. These observations validate the criteria that allow one to decide if an observed association constant of an oligomer to an RNA molecule reflects a complete complementarity between the two or not. This improves the basis for using oligonucleotide binding constants to RNA of known sequences to map secondary structure.     

The binding affinity of Ff gene 5 protein depends on the nearest-neighbor composition of the ssDNA substrate

The Ff gene 5 protein (g5p) is considered to be a nonspecific single-stranded DNA binding protein, because it binds cooperatively to and saturates the Ff bacteriophage single-stranded DNA genome and other single-stranded polynucleotides. However, the binding affinity Komega (the intrinsic binding constant times a cooperativity factor) differs by over an order of magnitude for binding to single-stranded polynucleotides such as poly[d(A)] and poly[d(C)]. A polynucleotide that is more stacked, like poly[d(A)], binds more weakly than one that is less stacked, like poly[d(C)]. To test the hypothesis that DNA base stacking, a nearest-neighbor property, is involved in the binding affinity of the Ff g5p for different DNA sequences, Komega values were determined as a function of NaCl concentration for binding to six synthetic sequences 48 nucleotides in length: dA48, dC48, d(AAC)16, d(ACC)16, d(AACC)12, and d(AAACC)9A3. The binding affinities of the protein for these sequences were indeed found to be related to the nearest-neighbor compositions of the sequences, rather than to simple base compositions. That is, the g5p binding site, which is spanned by four nucleotides, discriminates among these sequences on the basis of the relative numbers of nearest neighbors (AA, CC, and AC plus CA) in the sequence. The results support the hypothesis that the extent of base stacking/unstacking of the free, nonbound ssDNA plays an important role in the binding affinity of the Ff gene 5 protein.     

Analysis of local helix geometry in three B-DNA decamers and eight dodecamers

Local variations in B-DNA helix structure are compared among three decamers and eight dodecamers, which contain examples of all ten base-pair step types. All pairwise combinations of helix parameters are compared by linear regression analysis, in a search for internal relationships as well as correlations with base sequence. The primary conclusions are: (1) Three-center hydrogen bonds between base-pairs occur frequently in the major groove at C-C, C-A, A-A and A-C steps, but are less convincing at C-C and C-T steps in the minor groove. The requirements for large base-pair propeller are (1) that the base-pair should be A.T rather than G.C, and (2) that it be involved in a major groove three-center hydrogen bond with the following base-pair. Either condition alone is insufficient. Hence, a large propeller is expected at the leading base-pair of A-A and A-C steps, but not at A-T, T-A, C-A or C-C steps. (2) A systematic and quantitative linkage exists between helix variables twist, rise, cup and roll, of such strength that the rise between base-pairs can hardly be described as an independent variable at all. Two typical patterns of behavior are observed at steps from one base-pair to the next: high twist profile (HTP), characterized by high twist, low rise, positive cup and negative roll, and low twist profile (LTP), marked by low twist, high rise; negative cup and positive roll. Examples of HTP are steps G-C, G-A and Y-C-A-R, where Y is pyrimidine and R is purine. Examples of LTP steps are C-G, G-G, A-G and C-A steps other than Y-C-A-R. (3) The minor groove is especially narrow across the two base-pairs of the following steps: A-T, T-A, A-A and G-A. (4) In general, base step geometry cannot be correlated solely with the bases that define the step in question; the two flanking steps also must be taken into account. Hence, local helix structure must be studied in the context, not of two base-pairs: A-B, but of four: x-A-B-y. Calladine's rules, although too simple in detail, were correct in defining the length of sequence over which a given perturbation is expressed. Whereas ten different two-base steps are possible, allowing for the identity of complementary sequences, there are 136 different four-base steps. Only 33 of these 136 four-base steps are represented in the decamer and dodecamer structures solved to date, and hence it is premature to try to set up detailed structural algorithms. (5) The sugar-phosphate backbone chains of B-DNA place strong limits on sequence-induced structural variation, damping down most variables within four or five base-pairs, and preventing purine-purine anti-anti mismatches from causing bulges in the double helix. Hence, although short-range sequence-induced deformations (or deformability) are observed, long-range deformations propagated down the helix are not to be expected.     

Peptide-induced parallel DNA duplexes for oligopyrimidines. Stereospecificity in complexation for oligo(L-lysine) and oligo(L-ornithine)

It is shown that the cationic oligopeptides octadeca(L-lysine) (Lys18) and octadeca(L-ornithine) (Orn18) can induce a parallel duplex for the natural DNA oligomer dT10 with thymine-thymine base pairs. Complexation of the ammonium groups in the peptide side chains with the DNA phosphates leads to diminished electrostatic phosphate-phosphate repulsions, which allows this T-T base pair formation. From combined NOESY 1H NMR and molecular mechanics studies, it follows that the parallel duplex is right-handed, with the peptide located in the groove of the duplex. For the natural DNA oligomers dC10, d(C6T6), and d(T6C2T2), only Lys18 is able to induce the formation of parallel duplexes with C-C and T-T base pairs. It is shown that, for Orn18, a complexation must occur with one of the nonbonded oxygen atoms in the phosphate groups (OR) in such a way that unfavorable steric interactions are present with the C-C base pairs, which have a larger propellor twist angle than T-T base pairs. An analogy is presented between peptide complexation with the phosphates and the neutralization of the phosphate groups by methylation, which is known to lead to parallel duplexes with T-T base pairs (for both the Sp and Rp configurations) and C-C base pairs (only for the Sp configuration).     

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.     

The nature of preferred hairpin structures in 16S-like rRNA variable regions.

Variable length hairpins in 16S-like rRNA show a predominance for tetra-loops, its degree correlates with the protein content of the ribosome. The number of base-pairs adjacent to the loop (the tip size) and the nearest neighbor composition contribute to the stability of hairpin structures. The average tip size in length variable hairpins correlates with the thermophilicity of the organism, i.e. in temperate environments less stable stem structures are tolerated or even necessary. The most abundant loop families UUCG, GCAA, and CUUG occur most frequently at loop sizes 3, 2, and 7, respectively. Short tips of size less than or equal to 4 generally prefer nearest-neighbor combinations that result in CCC-GGG. Loop-specific tipmost nearest neighbors are revealed at longer tips: CUC(UUCG)GAG, GUA(GCAA)UAC with a maximum at tip sizes 5-6, and GWG(CUUG)CWC. Conserved hairpins, however, prefer variants of the UUCG and GCAA motifs with additional purines. Minor loop families and single motifs such as UUUA, UUUU, CUUGU, UUCGG, and UUU are investigated for preferable tip sizes and nearest-neighbor composition. Specific features are revealed for prominent hexa-loops.     

Effects of ssDNA sequences on non-sequence-specific protein binding

The circular dichroism (CD) spectra of single-stranded DNAs (ssDNAs) are significantly perturbed by the binding of single-stranded DNA binding proteins such as the Ff bacteriophage gene 5 protein (g5p) and the A domain of the 70 kDa subunit of human replication protein A (RPA70-A). These two proteins have similar OB-fold secondary structures, although their CD spectra at wavelengths below 250 nm differ greatly. The spectrum of g5p is dominated by a tyrosyl L(a) band at 229 nm, while that of RPA70-A is dominated by its beta secondary structure. Despite differences in their inherent spectral properties, these two proteins similarly perturb the spectra of bound nucleic acid oligomers. CD spectra of free, non-protein-bound ssDNAs are dependent on interactions of the nearest-neighboring nucleotides in the sequence. The CD spectra (per mol of nucleotide) of simple repetitive sequences 48 nucleotides in length and containing simple combinations of A and C are related by nearest-neighbor equations. For example, 3 x Deltaepsilon[d(AAC)(16)] = 3 x Deltaepsilon[d(ACC)(16)] + Deltaepsilon[d(A)(48)] - Deltaepsilon[d(C)(48)]. Moreover, nearest-neighbor equations relate the spectra of ssDNAs when they are bound by g5p, indicating that each type of perturbed nearest neighbor has a similar average structure within the binding site of the protein.     

Evaluation of nearest-neighbor methods for detection of chimeric small-subunit rRNA sequences.

Detection of chimeric artifacts formed when PCR is used to retrieve naturally occurring small-subunit (SSU) rRNA sequences may rely on demonstrating that different sequence domains have different phylogenetic affiliations. We evaluated the CHECK_CHIMERA method of the Ribosomal Database Project and another method which we developed, both based on determining nearest neighbors of different sequence domains, for their ability to discern artificially generated SSU rRNA chimeras from authentic Ribosomal Database Project sequences. The reliability of both methods decreases when the parental sequences which contribute to chimera formation are more than 82 to 84% similar. Detection is also complicated by the occurrence of authentic SSU rRNA sequences that behave like chimeras. We developed a naive statistical test based on CHECK_CHIMERA output and used it to evaluate previously reported SSU rRNA chimeras. Application of this test also suggests that chimeras might be formed by retrieving SSU rRNAs as cDNA. The amount of uncertainty associated with nearest-neighbor analyses indicates that such tests alone are insufficient and that better methods are needed.     

Sequence specificity of DNA adenine methylase in the protozoan Tetrahymena thermophila.

The sequence specificity of the Tetrahymena DNA-adenine methylase was determined by nearest-neighbor analyses of in vivo and in vitro methylated DNA. In vivo all four common bases were found to the 5' side of N6-methyladenine, but only thymidine was 3'. Homologous DNA already methylated in vivo and heterologous Micrococcus luteus DNA were methylated in vitro by a partially purified DNA-adenine methylase activity isolated from Tetrahymena macronuclei. The in vitro-methylated sequence differed from the in vivo sequence in that both thymidine and cytosine were 3' nearest neighbors of N6-methyladenine.     

Thermodynamics of single strand DNA base stacking

The thermodynamics of the stacking to unstacking transitions of 24 single-stranded DNA sequences (ssDNA), 10-12 bases in length, in sodium phosphate buffer were determined from 10 to 95 degrees C, using differential scanning calorimetry (DSC). An additional 22 ssDNA sequences did not exhibit an S<-->U transition in this temperature range. The transition properties of the ssDNA sequences with 相似文献   

The analysis of stress-induced duplex destabilization in long genomic DNA sequences.

We present a method for calculating predicted locations and extents of stress-induced DNA duplex destabilization (SIDD) as functions of base sequence and stress level in long DNA molecules. The base pair denaturation energies are assigned individually, so the influences of near neighbors, methylated bases, adducts, or lesions can be included. Sample calculations indicate that copolymeric energetics give results that are close to those derived when full near-neighbor energetics are used; small but potentially informative differences occur only in the calculated SIDD properties of moderately destabilized regions. The method presented here for analyzing long sequences calculates the destabilization properties within windows of fixed length N, with successive windows displaced by an offset distance d(o). The final values of the relevant destabilization parameters for each base pair are calculated as weighted averages of the values computed for each window in which that base pair appears. This approach implicitly assumes that the strength of the direct coupling between remote base pairs that is induced by the imposed stress attenuates with their separation distance. This strategy enables calculations of the destabilization properties of DNA sequences of any length, up to and including complete chromosomes. We illustrate its utility by calculating the destabilization properties of the entire E. coli genomic DNA sequence. A preliminary analysis of the results shows that promoters are associated with SIDD regions in a highly statistically significant manner, suggesting that SIDD attributes may prove useful in the computational prediction of promoter locations in prokaryotes.     

DNA interstrand crosslink formation by mechlorethamine at a cytosine-cytosine mismatch pair: kinetics and sequence dependence

Expansion of the triplet repeat DNA sequence d[CGG]n.d[CCG]n is a characteristic of Fragile X syndrome, a human neurodegenerative disease. Stable intrastrand conformations formed by both d[CGG]n and d[CCG]n, and involving G-G and C-C mismatch pairs, respectively, are believed to be of importance in the development of the disease. We have shown previously that C-C mismatch pairs can be crosslinked covalently by mechlorethamine, a nitrogen mustard alkylating agent, and hence this reaction may be of value as a probe for conformers of d[CCG]n. To characterize the mechlorethamine C-C crosslink reaction further, here we report the kinetics and sequence dependence of formation of the crosslink species, using a series of model duplexes. The rate of reaction depends on the base sequence proximal to the C-C mismatch pair. Hence, in 19mer duplexes containing a central d[M4M3M2M1Cn1n2n3n4].d[N4N3N2N1Cm1m2m3m4] sequence, where M-m and N-n are complementary base pairs, the amount of crosslink increased with increasing G-C content of the eight base pairs neighboring the C-C mismatch and with the proximity of the G-C pairs to the C-C mismatch. Molecular dynamics simulations of the solvated duplexes provided an explanation of these data. Hence, for a C-C pair flanked by G-C base pairs the mismatched cytosine bases remain stacked within the duplex, but for a C-C pair flanked by A-T base pairs, the simulations suggested local opening of the duplex around the C-C pair, making it a less effective target for mechlorethamine.     

Lipid chains and cholesterol in model membranes: a Monte Carlo Study

The Monte Carlo method has been employed to study the equilibrium properties of a planar array of hydrocarbon chains interacting with a cholesterol molecule. The chains are arranged to model one monolayer of a lipid bilayer and within this monolayer are allowed to move laterally and change conformations by gauche rotations. In the simulation cell there are 90 lipid chains and a single cholesterol molecule. Periodic boundary conditions are imposed upon the cell. The primary results of the calculations are order parameter profiles for the C-C bonds. These are calculated for (i) all chains, (ii) the 6 chains which are nearest neighbors to the cholesterol, and (iii) the 12 chains which are next-nearest neighbors to the cholesterol. Calculations are carried out for C-14, C-16, and C-18 chains. The results show that cholesterol strongly affects the upper portions of the chains, leaving them less able to change conformations. For C-16 and C-18 chains, the chain termini of the cholesterol neighbors are more disordered than the bulk chain termini. The magnitude of the effect depends strongly on the chain length. The results suggest that the changes in the lipid phase transition caused by cholesterol are a consequence of each cholesterol hindering the rotameric freedom of five to seven lipid chains.     

Statistical analysis of pair-wise compatibility of spatially nearest neighbor and adjacent residues in alpha-helix and beta-strands: application to a minimal model for secondary structure prediction

Secondary structural elements like alpha-helix and beta-strands possess distinctly different structural features and thus the relative positioning of the nearest neighbor residues, and also the sequence-wise adjacent residues is important in determining the structural preference. In the present work we have statistically examined the pair-wise compatibility pattern of physically nearest neighbors and separately the adjacent residue pairs along the sequence in between the nearest neighbor partners in alpha-helices and beta-strands. It has been demonstrated that the patterns and hence, the physical basis of the compatibility of adjacent residue pairs and the spatially nearest neighbors are significantly different in most cases. The influence of tertiary contacts on the pair-wise compatibility is shown to be significant for beta-strands while it is small for alpha-helices. Based on the compatibility of physically nearest neighbors and the sequence-wise adjacent residue pairs, a minimal model has been constructed to predict the alpha-helices, beta-strands and coils of a protein from its sequence. Application of this method to 100 sequences shows that it has a predictive capability comparable to that of other more sophisticated statistical methods.