PubMed HealthSearch

Biomedical subjects

J SantaLucia

Publications and source records attributed to J SantaLucia.

At least 19 recordsLinked to original sources

Nearest-neighbor thermodynamics and NMR of DNA sequences with internal A.A, C.C, G.G, and T.T mismatches.

Thermodynamic measurements are reported for 51 DNA duplexes with A.A, C.C, G.G, and T.T single mismatches in all possible Watson-Crick contexts. These measurements were used to test the applicability of the nearest-neighbor model and to calculate the 16 unique nearest-neighbor parameters for the 4 single like with like base mismatches next to a Watson-Crick pair. The observed trend in stabilities of mismatches at 37 degrees C is G.G > T.T approximately A.A > C.C. The observed stability trend for the closing Watson-Crick pair on the 5' side of the mismatch is G.C >/= C.G >/= A.T >/= T.A. The mismatch contribution to duplex stability ranges from -2.22 kcal/mol for GGC.GGC to +2.66 kcal/mol for ACT.ACT. The mismatch nearest-neighbor parameters predict the measured thermodynamics with average deviations of DeltaG degrees 37 = 3.3%, DeltaH degrees = 7. 4%, DeltaS degrees = 8.1%, and TM = 1.1 degrees C. The imino proton region of 1-D NMR spectra shows that G.G and T.T mismatches form hydrogen-bonded structures that vary depending on the Watson-Crick context. The data reported here combined with our previous work provide for the first time a complete set of thermodynamic parameters for molecular recognition of DNA by DNA with or without single internal mismatches. The results are useful for primer design and understanding the mechanism of triplet repeat diseases.

Base Pair Mismatch

NMR solution structure of a DNA dodecamer containing single G.T mismatches.

The three-dimensional solution structure of the self-complementary DNA dodecamer (CGT_GACGT_TACG above GCAT_TGCAG_TGC] which contains the thermodynamically destabilizing [TG_A above AT_T] motif was determined using two-dimensional NMR spectroscopy and simulated annealing protocols. Relaxation matrix analysis methods were used to yield accurate NOE derived distance restraints. Scalar coupling constants for the sugar protons were determined by quantitative simulations of DQF-COSY cross-peaks and used to determine sugar pucker populations. Twenty refined structures starting from random geometries converged to an average pairwise root mean square deviation of 0.49 A. Back calculated NOEs give Rc and Rx factors of 0.38 and 0.088, respectively. The final structure shows that each of the single G@T mismatches form a wobble pair with two hydrogen bonds where the guanine projects into the minor groove and the thymine projects into the major groove. The incorporation of the destabilizing [TG_A above AT_T] motif has little effect on the backbone torsion angles and helical parameters compared to standard B-form duplexes, which may explain why G.T mismatches are among the most commonly observed in DNA. The structure shows that perturbations caused by a G.T mismatch extend only to its neighboring Watson-Crick base pair, thus providing a structural basis for the applicability of the nearest-neighbor model to the thermodynamics of internal G.T mismatches.

Base Pair Mismatch

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.

Base Pairing

Nearest-neighbor thermodynamics of internal A.C mismatches in DNA: sequence dependence and pH effects.

Thermodynamics of 27 oligonucleotides with internal A.C mismatches at two different pHs were determined from UV absorbance versus temperature melting profiles. The data were combined with four literature values and used to derive nearest-neighbor parameters for all 16 trimer sequences with internal A.C mismatches at pH 7.0 and 5. 0. The results indicate that the contribution of single A.C mismatches to duplex stability is strongly dependent on the solution pH and the nearest-neighbor context. On average, the protonation of an internal A.C mismatch stabilizes the duplex by 1.39 kcal/mol for DeltaG degrees37 and 7.0 degreesC for the TM. The nearest-neighbor parameters predict DeltaG degrees37, DeltaH degrees, DeltaS degrees, and TM of oligonucleotides presented in this study with average deviations of 6.3%, 11.0%, 12.2%, and 1.8 degreesC, respectively, at pH 7.0 and 4.7%, 5.9%, 6.1%, and 1.3 degreesC, respectively, at pH 5. 0. At pH 7.0, the contribution of single A.C mismatches to helix stability ranges from 2.25 kcal/mol for TCA/AAT to 1.22 kcal/mol for GCG/CAC. At pH 5.0, however, the contribution of A+.C mismatches ranges from 1.09 kcal/mol for TCT/AAA to -0.43 kcal/mol for GCC/CAG. Implications of the results for replication fidelity and mismatch repair are discussed.

Adenine

Effect of mismatched complementary strands and 5'-change in sequence context on the thermodynamics and structure of benzo[a]pyrene-modified oligonucleotides.

Benzo[a]pyrene (B[a]P) is a well-studied environmental carcinogen that when activated can react with DNA to form four major adducts: (+)-trans-, (-)-trans-, (+)-cis-, and (-)-cis-anti-B[a]P-dG. In this study, two oligonucleotides (5'-dCCATT-GB[a]P-CTACC-3' and 5'-dCCATC-GB[a]P-CTACC-3') were prepared, each containing the four isomeric adducts, and these were hybridized to either complementary sequences or to sequences containing an A, G, or T opposite the adducted guanine. Thermal melting curves, CD, and UV spectra of each duplex were measured and compared with the unmodified counterpart. The raw and relative thermodynamic measurements were then compared which indicated that differences occur that are both adduct and sequence dependent. These differences were next compared with the in vitro DNA polymerase incorporation data and were found to be strikingly correlated. Most significantly, for all four B[a]P isomers a mismatch of an A across from the adduct resulted in the least amount of relative destabilization, while the Watson-Crick complement C showed the most; in vitro studies showed that A is the preferred base incorporated across from each isomer, while C was incorporated least often. This observed correlation suggests that one factor contributing to misincorporation at an adduct site is the thermodynamic stability of the incorporated base. Structurally, the effect of sequence context and mismatched complementary strands were also compared, suggesting that all adducts tend to intercalate within the helix when they are complemented with a mismatched complementary strand. In addition, the level of this intercalation seems to be both sequence and stereoisomer dependent.

Base Sequence

Thermodynamics of internal C.T mismatches in DNA.

Thermodynamics of 23 oligonucleotides with internal single C.T mismatches were obtained by measuring UV absorbance as a function of temperature. Results from these 23 duplexes were combined with three measurements from the literature to derive nearest-neighbor thermodynamic parameters for seven linearly independent trimer sequences with internal C.T mismatches. The data show that the nearest-neighbor model is adequate for predicting thermodynamics of oligonucleotides with internal C.T with average deviations for Delta G degrees37, Delta H degrees, Delta S degrees and T m of 6.4%, 9.9%, 10.6%, and 1.9 degreesC respectively. C.T mismatches destabilize the duplex in all sequence contexts. The thermodynamic contribution of C. T mismatches to duplex stability varies weakly depending on the orientation of the mismatch and its context and ranges from +1.02 kcal/mol for GCG/CTC and CCG/GTC to +1.95 kcal/mol for TCC/ATG.

Cytidine

Nearest neighbor thermodynamic parameters for internal G.A mismatches in DNA.

Thermodynamics of 22 oligonucleotides with internal single G.A mismatches dissolved in 1 M NaCl were determined from absorbance versus temperature melting curves. These data, combined with five literature sequences, were used to derive nearest-neighbor thermodynamic parameters for seven linearly independent trimer sequences with internal G.A mismatches and Watson-Crick flanking base pairs. The G.A mismatch parameters predict DeltaG degrees 37, DeltaH degrees, DeltaS degrees, and TM with average deviations of 4.4%, 7.4%, 8.0%, and 1.5 degrees C, respectively. The nearest-neighbor parameters show that G.A mismatch stability is strongly context dependent, and DeltaG degrees 37 ranges from +1.16 kcal/mol for TGA/AAT to -0.78 kcal/mol for GGC/CAG. In addition, one-dimensional 1H NMR spectra show that the G.A pairing geometry is pH and context dependent.

Base Composition

A unified view of polymer, dumbbell, and oligonucleotide DNA nearest-neighbor thermodynamics.

A unified view of polymer, dumbbell, and oligonucleotide nearest-neighbor (NN) thermodynamics is presented. DNA NN DeltaG degrees 37 parameters from seven laboratories are presented in the same format so that careful comparisons can be made. The seven studies used data from natural polymers, synthetic polymers, oligonucleotide dumbbells, and oligonucleotide duplexes to derive NN parameters; used different methods of data analysis; used different salt concentrations; and presented the NN thermodynamics in different formats. As a result of these differences, there has been much confusion regarding the NN thermodynamics of DNA polymers and oligomers. Herein I show that six of the studies are actually in remarkable agreement with one another and explanations are provided in cases where discrepancies remain. Further, a single set of parameters, derived from 108 oligonucleotide duplexes, adequately describes polymer and oligomer thermodynamics. Empirical salt dependencies are also derived for oligonucleotides and polymers.

DNA

The crystal structure of an RNA oligomer incorporating tandem adenosine-inosine mismatches.

The X-ray crystallographic structure of the RNA duplex [r(CGCAIGCG)]2 has been refined to 2.5 A. It shows a symmetric internal loop of two non-Watson-Crick base pairs which form in the middle of the duplex. The tandem A-I/I-A pairs are related by a crystallographic two-fold axis. Both A(anti)-I(anti) mismatches are in a head-to-head conformation forming hydrogen bonds using the Watson-Crick positions. The octamer duplexes stack above one another in the cell forming a pseudo-infinite helix throughout the crystal. A hydrated calcium ion bridges between the 3'-terminal of one molecule and the backbone of another. The tandem A-I mismatches are incorporated with only minor distortion to the backbone. This is in contrast to the large helical perturbations often produced by sheared G-A pairs in RNA oligonucleotides.

Adenosine

Thermodynamics and NMR of internal G.T mismatches in DNA.

Thermodynamics of 39 oligonucleotides with internal G.T mismatches dissolved in 1 M NaCl were determined from UV absorbance versus temperature profiles. These data were combined with literature values of six sequences to derive parameters for 10 linearly independent trimer and tetramer sequences with G.T mismatches and Watson-Crick base pairs. The G.T mismatch parameters predict DeltaG degrees 37, DeltaH degrees , DeltaS degrees , and TM with average deviations of 5.1%, 7.5%, 8.0%, and 1.4 degrees C, respectively. These predictions are within the limits of what can be expected for a nearest-neighbor model. The data show that the contribution of a single G.T mismatch to helix stability is context dependent and ranges from +1.05 kcal/mol for AGA/TTT to -1.05 kcal/mol for CGC/GTG. Several tests of the applicability of the nearest-neighbor model to G.T mismatches are described. Analysis of imino proton chemical shifts show that structural perturbations from the G.T mismatches are highly localized. One-dimensional NOE difference spectra demonstrate that G.T mismatches form stable hydrogen-bonded wobble pairs in diverse contexts. Refined nearest-neighbor parameters for Watson-Crick base pairs are also presented.

Base Composition

In vivo determination of RNA structure-function relationships: analysis of the 790 loop in ribosomal RNA.

The 790 loop is a conserved hairpin located between positions 786 and 796 of Escherichia coli 16 S rRNA that is required for ribosome function. Using a novel genetic approach, all positions in the loop were simultaneously mutated and functional mutant sequences were selected in vivo. This "instant evolution" experiment revealed that approximately 190 of the 262,144 possible mutant sequences were functional. Analysis of functional mutant sequences allowed discrimination between nucleotides directly involved in protein synthesis and those involved primarily in loop structure. Among the functional mutant sequences, positions 789 and 791 were invariant and extensive covariation was observed among the nucleotides at the base of the loop at positions 787, 788, 794 and 795. NMR and thermodynamic analyses of model 790 hairpins in vitro revealed weak pairing interactions between positions 787 and 795 and between positions 788 and 794 consistent with the in vivo mutational analysis. Functional analysis of site-directed mutants containing all possible nucleotide combinations at positions 787 and 795 in vivo showed that stable base-pairs at these positions prevent subunit association.

Chloramphenicol

Solution structure of (rGGCAGGCC)2 by two-dimensional NMR and the iterative relaxation matrix approach.

The three-dimensional solution structure of the RNA self-complementary duplex [formula in text] was derived from two-dimensional NMR and the iterative relaxation matrix approach. Each GA mismatch forms two hydrogen bonds: A-NH6 to G-O6 and A-N1 to G-NH1 (imino). This GA structure differs from the sheared tandem GA structure in [formula in text] which also has two hydrogen bonds: A-N7 to G-NH2 and A-NH6 to G-N3 [SantaLucia, J., Jr & Turner, D. H. (1993) Biochemistry 32, 12612-12623], although the only difference between the two sequences is the order of the two GA mismatches. Inspection of three-dimensional structures indicates that substituting [formula in text] for [formula in text] makes GA mismatches unable to have stable sheared conformations. This may explain why the [formula in text] motif is rarely observed in nature, whereas [formula in text] common.

Crystallography, X-Ray

Measuring the thermodynamics of RNA secondary structure formation.

The thermodynamics of RNA secondary structure formation in small model systems provides a database for predicting RNA structure from sequence. Methods for making these measurements are reviewed with emphasis on optical methods and treatment of experimental errors. Analysis of experimental results in terms of simple nearest-neighbor models is presented. Some measured sequence dependences of non-Watson-Crick motifs are discussed.

Base Composition

Improved nearest-neighbor parameters for predicting DNA duplex stability.

Thermodynamic data were determined from UV absorbance vs temperature profiles of 23 oligonucleotides. These data were combined with data from the literature for 21 sequences to derive improved parameters for the 10 Watson-Crick nearest neighbors. The observed trend in nearest-neighbor stabilities at 37 degrees C is GC > CG > GG > GA approximately GT approximately CA > CT > AA > AT > TA (where only the top strand is shown for each nearest neighbor). This trend suggests that both sequence and base composition are important determinants of DNA duplex stability. On average, the improved parameters predict deltaG degrees(37), deltaH degrees, deltaS degrees, and T(m) within 4%, 7%, 8%, and 2 degrees C, respectively. The parameters are optimized for the prediction of oligonucleotides dissolved in 1 M NaC1.

Base Sequence

Synthesis and NMR of RNA with selective isotopic enrichment in the bases.

Efficient syntheses of pyrimidine and purine nucleosides and nucleotides with selective 13C enrichment in the base moieties are described. Uridine and cytidine are labeled at position C6 and adenosine and guanosine are labeled at position C8. The selectively labeled nucleosides were converted to nucleoside triphosphates and used with in vitro transcription to synthesize labeled RNA. Isotope-edited 12C and 13C sub-spectra of a omega 1-1/2-X-filtered NOESY experiment are demonstrated to be useful for making resonance assignments and for deriving structural information in large (> 20 nt) RNA molecules. The labeled RNAs also allow heteronuclear J-couplings and relaxation parameters to be measured without complications from 13C-13C J-couplings.

Adenine

Structure of (rGGCGAGCC)2 in solution from NMR and restrained molecular dynamics.

The duplex (rGGCGAGCC)2 contains tandem G x A mismatches--a common motif in the secondary structures of biological RNAs. The three-dimensional structure of (rGGCGAGCC)2 was derived using molecular dynamics and energy minimization with NMR-derived restraints for 78 interproton distances (per strand), 18 hydrogen bonds for the six Watson-Crick G x C pairs, and 26 dihedral angles (per strand). The G x A mismatch structures are similar to those observed in a DNA duplex [Li, Y., Zon, G., & Wilson, W. D. (1991) Proc. Natl. Acad. Sci. U.S.A. 80, 26-30] and an RNA hairpin [Heus, H. A., & Pardi, A. (1991) Science 253, 191-193], with hydrogen bonds from guanine 2-amino and N3 to adenine N7 and 6-amino, respectively. The other G 2-amino and A 6-amino protons are within hydrogen-bonding distance of a phosphate oxygen and 2'-oxygen, respectively. Strong interstrand A-A and G-G stacking is observed between the G x A mismatches. This contrasts with the poor stacking observed between the G x A mismatches and closing G x C base pairs. The stems are basically A-form with all bases in the anti conformation and all nonterminal sugars in the C3'-endo conformation. The structure rationalizes previous thermodynamic, circular dichroism, and imino proton NMR results and suggests tandem G x A mismatches in RNA may provide a contact site for tertiary interactions.

Base Sequence

Context dependence of hydrogen bond free energy revealed by substitutions in an RNA hairpin.

Prediction and modeling of RNA structure requires knowledge of the free energy contributions of various interactions. Many unusual hydrogen bonds were recently proposed in the structure of a GCAA hairpin determined from nuclear magnetic resonance. The contributions of these hydrogen bonds to the folding stability of the hairpin formed by rGGCGCAAGCC have now been investigated through the use of functional group substitutions. These and previous results suggest a strong context dependence for the free energy of hydrogen bond formation. The results also suggest that the phylogenetic preference for GNRA (where N = A, C, G, or U and R = A or G) tetraloops may have a functional rather than thermodynamic basis.

Base Sequence