New evidence that the hydrophobic effect and dispersion are not major driving forces for nucleotide base stacking

"End-to-end" stacking of small dsRNA

PELDOR (pulsed electron-electron double resonance) is an established method to study intramolecular distances and can give evidence for conformational changes and flexibilities. However, it can also be used to study intermolecular interactions as for example oligerimization. Here, we used PELDOR to study the "end-to-end" stacking of small double-stranded (ds) RNAs. For this study, the dsRNA molecules were only singly labeled with the spin label TPA to avoid multispin effects and to measure only the intermolecular stacking interactions. It can be shown that small dsRNAs tend to assemble to rod-like structures due to π-π interactions between the base pairs at the end of the strands. On the one hand, these interactions can influence or complicate measurements aimed at the determining of the structure and dynamics of the dsRNA molecule itself. On the other hand, it can be interesting to study such intermolecular stacking interactions in more detail, as for example their dependence on ion concentration. We quantitatively determined the stacking probability as a function of the monovalent NaCl salt and the dsRNA concentration. From these data, the dissociation constant K_d was deduced and found to depend on the ratio between the NaCl salt and dsRNA concentrations. Additionally, the distances and distance distributions obtained predict a model for the stacking geometry of dsRNAs. Introducing a nucleotide overhangs at one end of the dsRNA molecule restricts the stacking to the other end, leading only to dimer formations. Introducing such an overhang at both ends of the dsRNA molecule fully suppresses stacking, as we demonstrate by PELDOR experiments quantitatively.

Geometric Patterns for Neighboring Bases Near the Stacked State in Nucleic Acid Strands

Structural variation in base stacking has been analyzed frequently in isolated double helical contexts for nucleic acids, but not as often in nonhelical geometries or in complex biomolecular environments. In this study, conformations of two neighboring bases near their stacked state in any environment are comprehensively characterized for single-strand dinucleotide (SSD) nucleic acid crystal structure conformations. An ensemble clustering method is used to identify a reduced set of representative stacking geometries based on pairwise distances between select atoms in consecutive bases, with multiple separable conformational clusters obtained for categories divided by nucleic acid type (DNA/RNA), SSD sequence, stacking face orientation, and the presence or absence of a protein environment. For both DNA and RNA, SSD conformations are observed that are either close to the A-form, or close to the B-form, or intermediate between the two forms, or further away from either form, illustrating the local structural heterogeneity near the stacked state. Among this large variety of distinct conformations, several common stacking patterns are observed between DNA and RNA, and between nucleic acids in isolation or in complex with proteins, suggesting that these might be stable stacking orientations. Noncanonical face/face orientations of the two bases are also observed for neighboring bases in the same strand, but their frequency is much lower, with multiple SSD sequences across categories showing no occurrences of such unusual stacked conformations. The resulting reduced set of stacking geometries is directly useful for stacking-energy comparisons between empirical force fields, prediction of plausible localized variations in single-strand structures near their canonical states, and identification of analogous stacking patterns in newly solved nucleic acid containing structures.

Unraveling Base Stacking Driving Forces in DNA

Base stacking is a key determinant of nucleic acid structures, but the precise origin of the thermodynamic driving force behind the stacking of nucleobases remains open. The rather mild stacking free energy measured experimentally, roughly a kcal/mol depending on the identity of the bases, is physiologically significant because while base stacking confers stability to the genome in its double helix form, the duplex also has to be unwound in order to be replicated or transcribed. A stacking free energy that is either too high or too low will over- or understabilize the genome, impacting the storage of genetic information and also its retrieval. While the molecular origin of stacking driving force has been attributed to many different sources including dispersion, electrostatics, and solvent hydrogen bonding, here we show via a systematic decomposition of the stacking free energy using large-scale computer simulations that the dominant driving force stabilizing base stacking is nonhydrophobic solvent entropy. Counteracting this is the conformational entropic penalty on the sugar-phosphate backbone against stacking, while solvent hydrogen-bonding, charge-charge interactions, and dispersive forces produce only secondary perturbations. Solvent entropic forces and DNA backbone conformational strains therefore work against each other, leading to a very mild composite stacking free energy in agreement with experiments.

Probing phenylalanine environments in oligomeric structures with pentafluorophenylalanine and cyclohexylalanine

Stabilization of protein structures and protein-protein interactions are critical in the engineering of industrially useful enzymes and in the design of pharmaceutically valuable ligands. Hydrophobic interactions involving phenylalanine residues play crucial roles in protein stability and protein-protein/peptide interactions. To establish an effective method to explore the hydrophobic environments of phenylalanine residues, we present a strategy that uses pentafluorophenylalanine (F5Phe) and cyclohexylalanine (Cha). In this study, substitution of F5Phe or Cha for three Phe residues at positions 328, 338, and 341 in the tetramerization domain of the tumor suppressor protein p53 was performed. These residues are located at the interfaces of p53-p53 interactions and are important in the stabilization of the tetrameric structure. The stability of the p53 tetrameric structure did not change significantly when F5Phe-containing peptides at positions Phe328 or Phe338 were used. In contrast, the substitution of Cha for Phe341 in the hydrophobic core enhanced the stability of the tetrameric structure with a T(m) value of 100 degrees C. Phe328 and Phe338 interact with each other through pi-interactions, whereas Phe341 is buried in the surrounding alkyl side-chains of the hydrophobic core of the p53 tetramerization domain. Furthermore, high pressure-assisted denaturation analysis indicated improvement in the occupancy of the hydrophobic core. Considerable stabilization of the p53 tetramer was achieved by filling the identified cavity in the hydrophobic core of the p53 tetramer. The results indicate the status of the Phe residues, indicating that the "pair substitution" of Cha and F5Phe is highly suitable for probing the environments of Phe residues.

Predicting helical coaxial stacking in RNA multibranch loops

The hypothesis that RNA coaxial stacking can be predicted by free energy minimization using nearest-neighbor parameters is tested. The results show 58.2% positive predictive value (PPV) and 65.7% sensitivity for accuracy of the lowest free energy configuration compared with crystal structures. The probability of each stacking configuration can be predicted using a partition function calculation. Based on the dependence of accuracy on the calculated probability of the stacks, a probability threshold of 0.7 was chosen for predicting coaxial stacks. When scoring these likely stacks, the PPV was 66.7% at a sensitivity of 51.9%. It is observed that the coaxial stacks of helices that are not separated by unpaired nucleotides can be predicted with a significantly higher accuracy (74.0% PPV, 66.1% sensitivity) than the coaxial stacks mediated by noncanonical base pairs (55.9% PPV, 36.5% sensitivity). It is also shown that the prediction accuracy does not show any obvious trend with multibranch loop complexity as measured by three different parameters.

Structural specificity conferred by a group I RNA peripheral element

Like proteins, structured RNAs must specify a native conformation that is more stable than all other possible conformations. Local structure is much more stable for RNA than for protein, so it is likely that the principal challenge for RNA is to stabilize the native structure relative to misfolded and partially folded intermediates rather than unfolded structures. Many structured RNAs contain peripheral structural elements, which surround the core elements. Although it is clear that peripheral elements stabilize structure within RNAs that contain them, it has not yet been explored whether they specifically stabilize the native states relative to alternative folds. A two-piece version of the group I intron RNA from Tetrahymena is used here to show that the peripheral element P5abc binds to the native conformation of the rest of the RNA 50,000 times more tightly than it binds to a long-lived misfolded conformation. Thus, P5abc stabilizes the native conformation by approximately 6 kcal/mol relative to this misfolded conformation. Further, activity measurements show that for the RNA lacking P5abc, the native conformation is only marginally preferred over the misfolded conformation (<0.5 kcal/mol), indicating that the peripheral structure of this RNA is required to achieve a significant thermodynamic preference for the native state. Such "structural specificity" may be a general function of RNA peripheral domains.

Analysis of the pi-pi stacking interactions between the aminoglycoside antibiotic kinase APH(3')-IIIa and its nucleotide ligands

A key contact in the active site of an aminoglycoside phosphotransferase enzyme (APH(3')-IIIa) is a pi-pi stacking interaction between Tyr42 and the adenine ring of bound nucleotides. We investigated the prevalence of similar Tyr-adenine contacts and found that many different protein systems employ Tyr residues in the recognition of the adenine ring. The geometry of these stacking interactions suggests that electrostatics play a role in the attraction between these aromatic systems. Kinetic and calorimetric experiments on wild-type and mutant forms of APH(3')-IIIa yielded further experimental evidence of the importance of electrostatics in the adenine binding region and suggested that the stacking interaction contributes approximately 2 kcal/mol of binding energy. This type of information concerning the forces that govern nucleotide binding in APH(3')-IIIa will facilitate inhibitor design strategies that target the nucleotide binding site of APH-type enzymes.

Contribution of aromatic interactions to alpha-helix stability

The influence of natural and unnatural i, i + 4 aromatic side chain-side chain interactions on alpha-helix stability was determined in Ala-Lys host peptides by circular dichroism (CD). All interactions investigated provided some stability to the helix; however, phenylalanine-phenylalanine (F-F) and phenylalanine-pentafluorophenylalanine (F-f5F) interactions resulted in the greatest enhancement in helicity, doubling the helical content over i, i + 5 control peptides at internal positions. Quantification of these interactions using AGADIR multistate helix-coil algorithm revealed that the F-F and F-f5F interaction energies are equivalent at internal positions in the sequence (deltaGF-F = deltaGF-f5F = -0.27 kcal/mol), despite the differences in their expected geometries. As the strength of a face-to-face stacked phenyl-pentafluorophenyl interaction should surpass an edge-to-face or offset-stacked phenyl-phenyl interaction, we believe this result reflects the inability of the side chains in F-f5F to attain a fully stacked geometry within the context of an alpha-helix. Positioning the interactions at the C-terminus led to much stronger interactions (deltaGF-F = -0.8 kcal/mol; deltaGF-f5F = -0.55 kcal/mol) likely because of favorable chi(1) rotameric preferences for aromatic residues at C-capping regions of alpha-helices, suggesting that aromatic side chain-side chain interactions are an effective alpha-helix C-capping method.

Unexpected substituent effects in offset pi-pi stacked interactions in water

We have measured the rotational barriers of meta- and para-substituted N-benzyl-2-(2-fluorophenyl)pyridinium bromides in aqueous solution by dynamic NMR as a model system for offset-stacking interactions in proteins. Because the benzyl ring can stack with the 2-fluorophenyl ring in the offset conformation in the ground state, but not the transition state, the rotational barrier reflects the magnitude of the stacking interaction. Only a small (0.1 kcal/mol) change in rotational barrier was found for para substituents relative to hydrogen. A much larger energy difference was found for electronegative meta substituents (up to 0.66 kcal/mol for CF3). Evidence suggests that this is due at least in part to an electrostatic interaction between electron-poor hydrogens on one ring with the electronegative substituents on the other ring.

Parallel multiplex thermodynamic analysis of coaxial base stacking in DNA duplexes by oligodeoxyribonucleotide microchips

Parallel thermodynamic analysis of the coaxial stacking effect of two bases localized in one strand of DNA duplexes has been performed. Oligonucleotides were immobilized in an array of three-dimensional polyacrylamide gel pads of microchips (MAGIChips'). The stacking effect was studied for all combinations of two bases and assessed by measuring the increase in melting temperature and in the free energy of duplexes formed by 5mers stacked to microchip-immobilized 10mers. For any given interface, the effect was studied for perfectly paired bases, as well as terminal mismatches, single base overlaps, single and double gaps, and modified terminal bases. Thermodynamic parameters of contiguous stacking determined by using microchips closely correlated with data obtained in solution. The extension of immobilized oligonucleotides with 5,6-dihydroxyuridine, a urea derivative of deoxyribose, or by phosphate, decreased the stacking effect moderately, while extension with FITC or Texas Red virtually eliminated stacking. The extension of the immobilized oligonucleotides with either acridine or 5-nitroindole increased stacking to mispaired bases and in some GC-rich interfaces. The measurements of stacking parameters were performed in different melting buffers. Although melting temperatures of AT- and GC-rich oligonucleotides in 5 M tetramethylammonium chloride were equalized, the energy of stacking interaction was significantly diminished.

Stability and structure of RNA duplexes containing isoguanosine and isocytidine

Isoguanosine (iG) and isocytidine (iC) differ from guanosine (G) and cytidine (C), respectively, in that the amino and carbonyl groups are transposed. The thermodynamic properties of a set of iG, iC containing RNA duplexes have been measured by UV optical melting. It is found that iG-iC replacements usually stabilize duplexes, and the stabilization per iG-iC pair is sequence-dependent. The sequence dependence can be fit to a nearest-neighbor model in which the stabilities of iG--iC pairs depend on the adjacent iG--iC or G--C pairs. For 5'-CG-3'/3'-GC-5' and 5'-GG-3'/3'-CC-5' nearest neighbors, the free energy differences upon iG-iC replacement are smaller than 0.2 kcal/mol at 37 degrees C, regardless of the number of replacements. For 5'-GC-3'/3'-CG-5', however, each iG--iC replacement adds 0.6 kcal/mol stabilizing free energy at 37 degrees C. Stacking propensities of iG and iC as unpaired nucleotides at the end of a duplex are similar to those of G and C. An NMR structure is reported for r(CiGCGiCG)(2) and found to belong to the A-form family. The structure has substantial deviations from standard A-form but is similar to published NMR and/or crystal structures for r(CGCGCG)(2) and 2'-O-methyl (CGCGCG)(2). These results provide benchmarks for theoretical calculations aimed at understanding the fundamental physical basis for the thermodynamic stabilities of nucleic acid duplexes.

The physical basis of nucleic acid base stacking in water

It has been argued that the stacking of adenyl groups in water must be driven primarily by electrostatic interactions, based upon NMR data showing stacking for two adenyl groups joined by a 3-atom linker but not for two naphthyl groups joined by the same linker. In contrast, theoretical work has suggested that adenine stacking is driven primarily by nonelectrostatic forces, and that electrostatic interactions actually produce a net repulsion between adenines stacking in water. The present study provides evidence that the experimental data for the 3-atom-linked bis-adenyl and bis-naphthyl compounds are consistent with the theory indicating that nonelectrostatic interactions drive adenine stacking. First, a theoretical conformational analysis is found to reproduce the observed ranking of the stacking tendencies of the compounds studied experimentally. A geometric analysis identifies two possible reasons, other than stronger electrostatic interactions, why the 3-atom-linked bis-adenyl compounds should stack more than the bis-naphthyl compounds. First, stacked naphthyl groups tend to lie further apart than stacked adenyl groups, based upon both quantum calculations and crystal structures. This may prevent the bis-naphthyl compound from stacking as extensively as the bis-adenyl compound. Second, geometric analysis shows that more stacked conformations are sterically accessible to the bis-adenyl compound than to the bis-naphthyl compound because the linker is attached to the sides of the adenyl groups, but to the ends of the naphthyl groups. Finally, ab initio quantum mechanics calculations and energy decompositions for relevant conformations of adenine and naphthalene dimers support the view that stacking in these compounds is driven primarily by nonelectrostatic interactions. The present analysis illustrates the importance of considering all aspects of a molecular system when interpreting experimental data, and the value of computer models as an adjunct to chemical intuition.

Thermodynamic parameters for DNA sequences with dangling ends

The thermodynamic contributions to duplex formation of all 32 possible single-nucleotide dangling ends on a Watson-Crick pair are reported. In most instances, dangling ends are stabilizing with free energy contributions ranging from +0.48 (GT(A)) to-0.96 kcal/mol (). In comparison, Watson-Crick nearest-neighbor increments range from -0. 58 (TA/AT) to -2.24 (GC/CG) kcal/mol. Hence, in some cases, a dangling end contributes as much to duplex stability as a Watson-Crick A-T base pair. The implications of these results for DNA probe design are discussed. Analysis of the sequence dependence of dangling-end stabilities show that the nature of the closing base pair largely determines the stabilization. For a given closing base pair, however, adenine dangling ends are always more or equally as stable as the other dangling nucleotides. Moreover, 5' dangling ends are more or equally as stabilizing as their 3' counterparts. Comparison of DNA with RNA dangling-end motifs shows that DNA motifs with 5' dangling ends contribute to stability equally or more than their RNA counterparts. Conversely, RNA 3' dangling ends contribute to stability equally or more than their DNA counterparts. This data set has been incorporated into a DNA secondary structure prediction algorithm (DNA MFOLD) (http://mfold2.wustl.edu/mfold/dna/for m1.cgi) as well as a DNA hybridization prediction algorithm (HYTHERtrade mark) (http://jsl1.chem.wayne.edu/Hyther/hythermenu .html).

Factors Contributing to Aromatic Stacking in Water: Evaluation in the Context of DNA

We report the use of thermodynamic measurements in a self-complementary DNA duplex (5'-dXCGCGCG)(2), where X is an unpaired natural or nonnatural deoxynucleoside, to study the forces that stabilize aqueous aromatic stacking in the context of DNA. Thermal denaturation experiments show that the core duplex (lacking X) is formed with a free energy (37 °C) of -8.1 kcal·mol(-1) in a pH 7.0 buffer containing 1 M Na(+). We studied the effects of adding single dangling nucleosides (X) where the aromatic "base" is adenine, guanine, thymine, cytosine, pyrrole, benzene, 4-methylindole, 5-nitroindole, trimethylbenzene, difluorotoluene, naphthalene, phenanthrene, and pyrene. Adding these dangling residues is found to stabilize the duplex by an additional -0.8 to -3.4 kcal·mol(-1). At 5 μM DNA concentration, T(m) values range from 41.7 °C (core sequence) to 64.1 °C (with dangling pyrene residues). For the four natural bases, the order of stacking ability is A > G ≥ T = C. The nonpolar analogues stack more strongly in general than the more polar natural bases. The stacking geometry was confirmed in two cases (X = adenine and pyrene) by 2-D NOESY experiments. Also studied is the effect of ethanol cosolvent on the stacking of natural bases and pyrene. Stacking abilities were compared to calculated values for hydrophobicity, dipole moment, polarizability, and surface area. In general, hydrophobic effects are found to be larger than other effects stabilizing stacking (electrostatic effects, dispersion forces); however, the natural DNA bases are found to be less dependent on hydrophobic effects than are the more nonpolar compounds. The results also point out strategies for the design nucleoside analogues that stack considerably more strongly than the natural bases; such compounds may be useful in stabilizing designed DNA structures and complexes.