Thermodynamic analysis of conserved loop-stem interactions in P1-P2 frameshifting RNA pseudoknots from plant Luteoviridae

The RNA genomes of plant luteovirids beet western yellows virus (BWYV), potato leaf roll virus (PLRV), and pea enation mosaic virus (PEMV RNA1; PEMV-1) contain a short mRNA pseudoknotted motif overlapping the P1 and P2 open reading frames required for programmed -1 mRNA ribosomal frameshifting. The relationship between structure, stability, and function is poorly understood in these RNA systems. A m(5)-C(8)-substituted BWYV RNA is employed to establish that the BWYV P1-P2 pseudoknot is protonated at cytidine 8 in loop L1 (delta(N(3)H)+ = 12.98 ppm), which stabilizes a C(+.)(G-C) major groove base triple by Delta(DeltaG(37))(protonation) = 3.1 (+/-0.4) kcal mol(-1). The stabilities of both the PLRV and PEMV-1 P1-P2 pseudoknots are also strongly pH-dependent, with Delta(DeltaG(37))(protonation) = 2.1 (+/-0.2) kcal mol(-1) for the PEMV-1 pseudoknot despite a distinct structural context. As previously found for the BWYV pseudoknot [Nixon and Giedroc (2000) J. Mol. Biol. 296, 659], both the PLRV and PEMV-1 RNAs are stabilized by DeltaH > or = 30 kcal mol(-)(1) in excess of secondary structure predictions, attributed to loop L2-stem S1 minor groove triplex interactions. BWYV RNAs containing single 2'-deoxy or A --> G substitutions that disrupt L2-S1 hydrogen bonding are strongly destabilized with Delta(DeltaG(37))(folding) (pH = 7.0) ranging from approximately 1.8 (+/-0.3) to > or =4.0 kcal mol(-1), relative to the wild-type BWYV RNA. These findings suggest that each member of this family of pseudoknots adopts a tightly folded structure that maximizes the cooperativity and complementarity of L1-S2 and L2-S1 loop-stem interactions required in part to offset the low intrinsic stability of the short three base pair pseudoknot stem S2.     

The GAAA tetraloop-receptor interaction contributes differentially to folding thermodynamics and kinetics for the P4-P6 RNA domain

Tetraloops with the generic sequence GNRA are commonly found in RNA secondary structure, and interactions of such tetraloops with "receptors" elsewhere in RNA play important roles in RNA structure and folding. However, the contributions of tetraloop-receptor interactions specifically to the kinetics of RNA tertiary folding, rather than the thermodynamics of maintaining tertiary structure once folded, have not been reported. Here we investigate the role of the key GAAA tetraloop-receptor motif in folding of the P4-P6 domain of the Tetrahymena group I intron RNA. Insertions of one or more nucleotides into the tetraloop significantly disrupt the thermodynamics of tertiary folding; single-nucleotide insertions shift the folding free energy by 2-4 kcal/mol (DeltaDeltaG(o)'). The folding kinetics of several modified P4-P6 domains were determined by stopped-flow fluorescence spectroscopy, using an internally incorporated pyrene residue as the chromophore. In contrast to the thermodynamic results, the kinetics of Mg(2+)-induced folding were barely affected by the tetraloop modifications, with a DeltaDeltaG(++) of 0.2-0.4 kcal/mol and a Phi value (ratio of the kinetic and thermodynamic contributions) of <0.1. These data indicate an early transition state for folding of P4-P6 with respect to forming the tetraloop-receptor contact, consistent with previous results for modifications elsewhere in P4-P6. We conclude that the GAAA tetraloop-receptor motif contributes little to the stabilization of the transition state for Mg(2+)-induced P4-P6 folding. Rather, the tetraloop-receptor motif acts to clamp the RNA once folding has occurred. This is the first report to correlate the kinetic and thermodynamic contributions of an important RNA tertiary motif, the GNRA tetraloop-receptor. The results are related to possible models for the Mg(2+)-induced folding of the P4-P6 RNA, including a model invoking rapid nonspecific electrostatic collapse.     

Contribution of hydrogen bonding to protein stability estimated from isotope effects

An unresolved issue in structural biology concerns the relative contribution of H bonds to protein stability. We use the small molecules 4-acetamidobenzoic acid and N-acetylanthranilic acid as model compounds to relate the energetic contribution from hydrogen bonds (H bonds) to the deuterium/hydrogen amide isotope effect. N-Acetylanthranilic acid models carbonyl-amide H bonds formed during protein folding; 4-acetamidobenzoic acid models the unfolded state in which the amide H bonds to water. NMR is used to measure shifts in the pK(a) of the ionizable carboxyl group when the amides of the compounds are either protonated or deuterated. From the pK(a) shift, we obtain a quantitative scale factor: SF = partial partial differential(DeltaG(HB))/partial partial differential(RT ln Phi), where DeltaG(HB) is the change in free energy of an H bond upon isotope substitution and Phi is the fractionation factor. Isotope effect data also are reported for a small globular protein, lambda repressor, using the "C(m) experiment". The protein's isotope effect, which reports on the shape of the energy well, is converted to H-bonding free energy by applying the scale factor. We estimate that amide-related H bonds (amide-carbonyl and amide-water) contribute favorably to protein stability by approximately 30-50 kcal/mol in lambda repressor, GCN4 coiled coil, and cytochrome c but unfavorably by approximately 6 kcal/mol in ubiquitin. The results indicate that H-bond strength varies from one protein to another and presumably at different sites within the same protein.     

Sequence and structural conservation in RNA ribose zippers

The "ribose zipper", an important element of RNA tertiary structure, is characterized by consecutive hydrogen-bonding interactions between ribose 2'-hydroxyls from different regions of an RNA chain or between RNA chains. These tertiary contacts have previously been observed to also involve base-backbone and base-base interactions (A-minor type). We searched for ribose zipper tertiary interactions in the crystal structures of the large ribosomal subunit RNAs of Haloarcula marismortui and Deinococcus radiodurans, and the small ribosomal subunit RNA of Thermus thermophilus and identified a total of 97 ribose zippers. Of these, 20 were found in T. thermophilus 16 S rRNA, 44 in H. marismortui 23 S rRNA (plus 2 bridging 5 S and 23 S rRNAs) and 30 in D. radiodurans 23 S rRNA (plus 1 bridging 5 S and 23 S rRNAs). These were analyzed in terms of sequence conservation, structural conservation and stability, location in secondary structure, and phylogenetic conservation. Eleven types of ribose zippers were defined based on ribose-base interactions. Of these 11, seven were observed in the ribosomal RNAs. The most common of these is the canonical ribose zipper, originally observed in the P4-P6 group I intron fragment. All ribose zippers were formed by antiparallel chain interactions and only a single example extended beyond two residues, forming an overlapping ribose zipper of three consecutive residues near the small subunit A-site. Almost all ribose zippers link stem (Watson-Crick duplex) or stem-like (base-paired), with loop (external, internal, or junction) chain segments. About two-thirds of the observed ribose zippers interact with ribosomal proteins. Most of these ribosomal proteins bridge the ribose zipper chain segments with basic amino acid residues hydrogen bonding to the RNA backbone. Proteins involved in crucial ribosome function and in early stages of ribosomal assembly also stabilize ribose zipper interactions. All ribose zippers show strong sequence conservation both within these three ribosomal RNA structures and in a large database of aligned prokaryotic sequences. The physical basis of the sequence conservation is stacked base triples formed between consecutive base-pairs on the stem or stem-like segment with bases (often adenines) from the loop-side segment. These triples have previously been characterized as Type I and Type II A-minor motifs and are stabilized by base-base and base-ribose hydrogen bonds. The sequence and structure conservation of ribose zippers can be directly used in tertiary structure prediction and may have applications in molecular modeling and design.     

Energetics of a strongly pH dependent RNA tertiary structure in a frameshifting pseudoknot

Retroviruses employ -1 translational frameshifting to regulate the relative concentrations of structural and non-structural proteins critical to the viral life cycle. The 1.6 A crystal structure of the -1 frameshifting pseudoknot from beet western yellows virus reveals, in addition to Watson-Crick base-pairing, many loop-stem RNA tertiary structural interactions and a bound Na(+). Investigation of the thermodynamics of unfolding of the beet western yellows virus pseudoknot reveals strongly pH-dependent loop-stem tertiary structural interactions which stabilize the molecule, contributing a net of DeltaH approximately -30 kcal mol(-1) and DeltaG degrees (37) of -3.3 kcal mol(-1) to a total DeltaH and DeltaG degrees (37) of -121 and -16 kcal mol(-1), respectively, at pH 6.0, 0.5 M K(+) by DSC. Characterization of mutant RNAs supports the presence of a C8(+).G12-C26 loop 1-stem 2 base-triple (pK(a)=6.8), protonation of which contributes nearly -3.5 kcal mol(-1) in net stability in the presence of a wild-type loop 2. Substitution of the nucleotides in loop 2 with uridine bases, which would eliminate the minor groove triplex, destroys pseudoknot formation. An examination of the dependence of the monovalent ion and type on melting profiles suggests that tertiary structure unfolding occurs in a manner quantitatively consistent with previous studies on the stabilizing effects of K(+), NH(4)(+) and Na(+) on other simple duplex and pseudoknotted RNAs.     

Multiple folding pathways for the P4-P6 RNA domain

We recently described site-specific pyrene labeling of RNA to monitor Mg(2+)-dependent equilibrium formation of tertiary structure. Here we extend these studies to follow the folding kinetics of the 160-nucleotide P4-P6 domain of the Tetrahymena group I intron RNA, using stopped-flow fluorescence with approximately 1 ms time resolution. Pyrene-labeled P4-P6 was prepared using a new phosphoramidite that allows high-yield automated synthesis of oligoribonucleotides with pyrene incorporated at a specific 2'-amino-2'-deoxyuridine residue. P4-P6 forms its higher-order tertiary structure rapidly, with k(obs) = 15-31 s(-1) (t(1/2) approximately 20-50 ms) at 35 degrees C and [Mg(2+)] approximately 10 mM in Tris-borate (TB) buffer. The folding rate increases strongly with temperature from 4 to 45 degrees C, demonstrating a large activation enthalpy DeltaH(double dagger) approximately 26 kcal/mol; the activation entropy DeltaS(double dagger) is large and positive. In low ionic strength 10 mM sodium cacodylate buffer at 35 degrees C, a slow (t(1/2) approximately 1 s) folding component is also observed. The folding kinetics are both ionic strength- and temperature-dependent; the slow phase vanishes upon increasing [Na(+)] in the cacodylate buffer, and the kinetics switch completely from fast at 30 degrees C to slow at 40 degrees C. Using synchrotron hydroxyl radical footprinting, we confirm that fluorescence monitors the same kinetic events as hydroxyl radical cleavage, and we show that the previously reported slow P4-P6 folding kinetics apply only to low ionic strength conditions. One model to explain the fast and slow folding kinetics postulates that some tertiary interactions are present even without Mg(2+) in the initial state. The fast kinetic phase reflects folding that is facilitated by these interactions, whereas the slow kinetics are observed when these interactions are disrupted at lower ionic strength and higher temperature.     

Energetics of hydrogen bond networks in RNA: hydrogen bonds surrounding G+1 and U42 are the major determinants for the tertiary structure stability of the hairpin ribozyme

The hairpin ribozyme, a small catalytic RNA consisting of two helix-loop-helix motifs, serves as a paradigm for RNA folding. In the active conformer, the ribozyme is docked into a compact structure via loop-loop interactions. The crystal structure of the docked hairpin ribozyme shows an intricate network of hydrogen bonding interactions at the docking interface, mediated by the base, sugar, and phosphate groups of U42 and G+1 [Rupert, P. B., and Ferre-D'Amare, A. R. (2001) Nature 410, 780-786]. To elucidate the determinants for tertiary structure stability in the hairpin ribozyme, we evaluated the energetic contributions of hydrogen bonds surrounding U42 and G+1 by time-resolved fluorescence resonance energy transfer using modified ribozymes that lack one or more of the individual interactions. Elimination of a single tertiary hydrogen bond consistently resulted in a net destabilization of approximately 2 kJ/mol. The results of double- and triple-mutant cycles suggest that individual hydrogen bonds surrounding G+1 or U42 act cooperatively and form extended hydrogen bond networks that stabilize the docked ribozyme. These results demonstrate that RNAs, similar to proteins, can exploit coupled hydrogen bond networks to stabilize the docking of distant structural domains.     

Structural insight into the kinetics and DeltaCp of interactions between TEM-1 beta-lactamase and beta-lactamase inhibitory protein (BLIP)

In a previous study, we examined thermodynamic parameters for 20 alanine mutants in beta-lactamase inhibitory protein (BLIP) for binding to TEM-1 beta-lactamase. Here we have determined the structures of two thermodynamically distinctive complexes of BLIP mutants with TEM-1 beta-lactamase. The complex BLIP Y51A-TEM-1 is a tight binding complex with the most negative binding heat capacity change (DeltaG = approximately -13 kcal mol(-1) and DeltaCp = approximately -0.8 kcal mol(-1) K(-1)) among all of the mutants, whereas BLIP W150A-TEM-1 is a weak complex with one of the least negative binding heat capacity changes (DeltaG = approximately -8.5 kcal mol(-1) and DeltaCp = approximately -0.27 kcal mol(-1) K(-1)). We previously determined that BLIP Tyr51 is a canonical and Trp150 an anti-canonical TEM-1-contact residue, where canonical refers to the alanine substitution resulting in a matched change in the hydrophobicity of binding free energy. Structure determination indicates a rearrangement of the interactions between Asp49 of the W150A BLIP mutant and the catalytic pocket of TEM-1. The Asp49 of W150A moves more than 4 angstroms to form two new hydrogen bonds while losing four original hydrogen bonds. This explains the anti-canonical nature of the Trp150 to alanine substitution, and also reveals a strong long distance coupling between Trp150 and Asp49 of BLIP, because these two residues are more than 25 angstroms apart. Kinetic measurements indicate that the mutations influence the dissociation rate but not the association rate. Further analysis of the structures indicates that an increased number of interface-trapped water molecules correlate with poor interface packing in a mutant. It appears that the increase of interface-trapped water molecules is inversely correlated with negative binding heat capacity changes.     

Metal-activated histidine carbon donor hydrogen bonds contribute to metalloprotein folding and function

Carbon donor hydrogen bonds are typically weak interactions that contribute less than 2 kcal/mol, and provide only modest stabilization in proteins. One exception is the class of hydrogen bonds donated by heterocyclic side chain carbons. Histidine is capable of particularly strong interactions through the Cε¹ and Cδ² carbons when the imidazole is protonated or bound to metal. Given the frequent occurrence of metal-bound histidines in metalloproteins, we characterized the energies of these interactions through DFT calculations on model compounds. Imidazole-water hydrogen bonding could vary from −11.0 to −17.0 kcal/mol, depending on the metal identity and oxidation state. A geometric search of metalloprotein structures in the PDB identified a number of candidate His C-H···O hydrogen bonds which may be important for folding or function. DFT calculations on model complexes of superoxide reductase show a carbon donor hydrogen bond positioning a water molecule above the active site.     

Contribution of a buried hydrogen bond to lambda repressor folding kinetics.

A hydrogen bond between the buried residues Asp 14 and Ser 77 in monomeric lambda repressor has been removed by mutation of these residues to alanine. Double mutant cycles show that the interaction stabilizes the native state of the protein by 1.5 kcal/mol. Removal of the interaction affects mainly the unfolding rates and not the folding rates, suggesting that this hydrogen bond is not substantially formed in the rate-limiting steps in the folding pathways of the protein. Mutations in two versions of lambda6-85, wild type and the faster folding G46A/G48A (WT), show similar effects. Diffusion-collision correctly predicts the behavior of WT but not of wild type. Our analysis suggests that folding of helix 3 is a crucial slow step along the various folding pathways and generally occurs before the formation of the 14-77 hydrogen bond. Experiments removing tertiary interactions, combined with experiments altering helical stability and diffusion-collision calculations, provide a strategy to unravel the folding mechanisms of small helical proteins.     

An early transition state for folding of the P4-P6 RNA domain

Tertiary folding of the 160-nt P4-P6 domain of the Tetrahymena group I intron RNA involves burying of substantial surface area, providing a model for the folding of other large RNA domains involved in catalysis. Stopped-flow fluorescence was used to monitor the Mg2+-induced tertiary folding of pyrene-labeled P4-P6. At 35 degrees C with [Mg2+] approximately 10 mM, P4-P6 folds on the tens of milliseconds timescale with k(obs) = 15-31 s(-1). From these values, an activation free energy deltaG(double dagger) of approximately 8-16 kcal/mol is calculated, where the large range for deltaG(double dagger) arises from uncertainty in the pre-exponential factor relating k(obs) and delta G(double dagger). The folding rates of six mutant P4-P6 RNAs were measured and found to be similar to that of the wild-type RNA, in spite of significant thermodynamic destabilization or stabilization. The ratios of the kinetic and thermodynamic free energy changes phi = delta deltaG(double dagger)/delta deltaG(o') are approximately 0, implying a folding transition state in which most of the native-state tertiary contacts are not yet formed (an early folding transition state). The k(obs) depends on the Mg2+ concentration, and the initial slope of k(obs) versus [Mg2+] suggests that only approximately 1 Mg2+ ion is bound in the rate-limiting folding step. This is consistent with an early folding transition state, because folded P4-P6 binds many Mg2+ ions. The observation of a substantial deltaG(double dagger) despite an early folding transition state suggests that a simple two-state folding diagram for Mg2+-induced P4-P6 folding is incomplete. Our kinetic data are some of the first to provide quantitative values for an activation barrier and location of a transition state for tertiary folding of an RNA domain.     

Linkage of monovalent and divalent ion binding in the folding of the P4-P6 domain of the Tetrahymena ribozyme

We have explored the linkage of monovalent and divalent ion binding in the folding of the P4-P6 domain of Tetrahymena thermophila ribozyme by examining the Mg2+-induced folding and the urea-induced denaturation of the folded state as a function of Na+ under equilibrium folding conditions using hydroxyl radical footprinting. These studies allowed a thermodynamic examination of eight discrete protection sites within P4-P6 that are involved in several tertiary structure contacts. Monovalent ions compete with Mg2+ ions in mediating P4-P6 folding. The urea denaturation isotherms demonstrated DeltaDeltaG values of >2 kcal x mol(-1) in experiments conducted in 10 versus 200 mM NaCl at a constant 10 mM MgCl2. However, the individual-site isotherms reported by footprinting revealed that larger than average changes in DeltaG values were localized to specific sites within the Mg2+-rich A-bulge. The competitive effects of monovalent ions were less when K+ rather than Na+ was the monovalent cation present. This result indicates the importance of the specific K+ binding sites that are associated with AA-platform structures to P4-P6 folding and stability. These site-specific footprinting data provide quantitative and site-specific measurements of the ion-linked stability for P4-P6 that are interpreted with respect to crystallographic data.     

Contribution of intra- and intermolecular hydrogen bonds to the conformational stability of human lysozyme(,).

In globular proteins, there are intermolecular hydrogen bonds between protein and water molecules, and between water molecules, which are bound with the proteins, in addition to intramolecular hydrogen bonds. To estimate the contribution of these hydrogen bonds to the conformational stability of a protein, the thermodynamic parameters for denaturation and the crystal structures of five Thr to Val and five Thr to Ala mutant human lysozymes were determined. The denaturation Gibbs energy (DeltaG) of Thr to Val and Thr to Ala mutant proteins was changed from 4.0 to -5.6 kJ/mol and from 1.6 to -6.3 kJ/mol, respectively, compared with that of the wild-type protein. The contribution of hydrogen bonds to the stability (DeltaDeltaG(HB)) of the Thr and other mutant human lysozymes previously reported was extracted from the observed stability changes (DeltaDeltaG) with correction for changes in hydrophobicity and side chain conformational entropy between the wild-type and mutant structures. The estimation of the DeltaDeltaG(HB) values of all mutant proteins after removal of hydrogen bonds, including protein-water hydrogen bonds, indicates a favorable contribution of the intra- and intermolecular hydrogen bonds to the protein stability. The net contribution of an intramolecular hydrogen bond (DeltaG(HB[pp])), an intermolecular one between protein and ordered water molecules (DeltaG(HB[pw])), and an intermolecular one between ordered water molecules (DeltaG(HB[ww])) could be estimated to be 8. 5, 5.2, and 5.0 kJ/mol, respectively, for a 3 A long hydrogen bond. This result shows the different contributions to protein stability of intra- and intermolecular hydrogen bonds. The entropic cost due to the introduction of a water molecule (DeltaG(H)()2(O)) could be also estimated to be about 8 kJ/mol.     

Tertiary structure stability of the hairpin ribozyme in its natural and minimal forms: different energetic contributions from a ribose zipper motif

The hairpin catalytic motif in tobacco ringspot virus satellite RNA consists of two helix-loop-helix elements on two adjacent arms of a four-way helical junction. The bases essential for catalytic activity are located in the loops that are brought into proximity by a conformational change as a prerequisite for catalysis. The two loops interact via a ribose zipper motif involving the 2'-hydroxyls of A10, G11, A24, and C25 [Rupert, P. B., and Ferre d'Amare, A. R. (2001) Nature 401, 780-786]. To quantify the energetic importance of the ribose zipper hydrogen bonds, we have incorporated deoxy modifications at these four positions and determined the resulting destabilization of the docked conformer by means of time-resolved fluorescence resonance energy transfer. In a minimal form of the ribozyme, in which the loops are placed on the arms of a two-way helical junction, all modifications lead to a significant loss in tertiary structure stability and altered Mg2+ binding. Surprisingly, no significant destabilization was seen with the natural four-way junction ribozyme, suggesting that hydrogen bonding interactions involving the 2'-hydroxyls do not contribute to the stability of the docked conformer. These results suggest that the energetic contributions of ribose zipper hydrogen bonds are highly context dependent and differ significantly for the minimal and natural forms of the ribozyme.     

Contribution of cation-pi interactions to protein stability

Calculations predict that cation- interactions make an important contribution to protein stability. While there have been some attempts to experimentally measure strengths of cation-pi interactions using peptide model systems, much less experimental data are available for globular proteins. We have attempted to determine the magnitude of cation-pi interactions of Lys with aromatic amino acids in four different proteins (LIVBP, MBP, RBP, and Trx). In each case, Lys was replaced with Gln and Met. In a separate series of experiments, the aromatic amino acid in each cation-pi pair was replaced by Leu. Stabilities of wild-type (WT) and mutant proteins were characterized by both thermal and chemical denaturation. Gln and aromatic --> Leu mutants were consistently less stable than corresponding Met mutants, reflecting the nonisosteric nature of these substitutions. The strength of the cation-pi interaction was assessed by the value of the change in the free energy of unfolding [DeltaDeltaG(degrees) = DeltaG(degrees)(Met) - DeltaG(degrees)(WT)]. This ranged from +1.1 to -1.9 kcal/mol (average value -0.4 kcal/mol) at 298 K and +0.7 to -2.6 kcal/mol (average value -1.1 kcal/mol) at the Tm of each WT. It therefore appears that the strength of cation-pi interactions increases with temperature. In addition, the experimentally measured values are appreciably smaller in magnitude than calculated values with an average difference /DeltaG(degrees)expt - DeltaG(degrees)calc/av of 2.9 kcal/mol. At room temperature, the data indicate that cation-pi interactions are at best weakly stabilizing and in some cases are clearly destabilizing. However, at elevated temperatures, close to typical Tm's, cation-pi interactions are generally stabilizing.     

Sheared Aanti.Aanti base pairs in a destabilizing 2 x 2 internal loop: the NMR structure of 5'(rGGCAAGCCU)2

The 5'(rGGCAAGCCU)(2) duplex contains tandem A.A pairs. The three-dimensional structure of the 5'(rGGCAAGCCU)(2) duplex was modeled by molecular dynamics and energy minimization with NMR-derived distance and dihedral angle restraints. Although the 5'(rCAAG)(2) loop is thermodynamically destabilizing by 1.1 kcal/mol, the tandem A.A pairs adopt a predominant conformation: a sheared anti-anti (A.A trans Hoogsteen/Sugar-edge) alignment similar to that observed in the crystal structure of the P4-P6 domain of the Tetrahymena thermophila intron [Cate, J. H., Gooding, A. R., Podell, E., Zhou, K., Golden, B. L., Kundrot, C. E., Cech, T. R., and Doudna, J. A. (1996) Science 273, 1678-1685]. The NMR-derived structure of the 5'(rGGCAAGCCU)(2) duplex exhibits cross-strand hydrogen bonds from N3 of A4 to an amino hydrogen of A5 and from the 2' oxygen of the A4 sugar to the other amino hydrogen of A5. An intrastrand hydrogen bond is formed from the 2' OH hydrogen of A4 to O5' of A5. The cross-strand A5 bases are stacked. The Watson-Crick G-C regions are essentially A-form. The sheared anti-anti (A.A trans Hoogsteen/Sugar-edge) alignment provides potential contact sites for tertiary interactions and, therefore, is a possible target site for therapeutics. Thus, thermodynamically destabilizing internal loops can be preorganized for tertiary interactions or ligand binding.     

Nine hydrophobic side chains are key determinants of the thermodynamic stability and oligomerization status of tumour suppressor p53 tetramerization domain.

The contribution of almost each amino acid side chain to the thermodynamic stability of the tetramerization domain (residues 326-353) of human p53 has been quantitated using 25 mutants with single-residue truncations to alanine (or glycine). Truncation of either Leu344 or Leu348 buried at the tetramer interface, but not of any other residue, led to the formation of dimers of moderate stability (8-9 kcal/mol of dimer) instead of tetramers. One-third of the substitutions were moderately destabilizing (<3.9 kcal/mol of tetramer). Truncations of Arg333, Asn345 or Glu349 involved in intermonomer hydrogen bonds, Ala347 at the tetramer interface or Thr329 were more destabilizing (4.1-5.7 kcal/mol). Strongly destabilizing (8.8- 11.7 kcal/mol) substitutions included those of Met340 at the tetramer interface and Phe328, Arg337 and Phe338 involved peripherally in the hydrophobic core. Truncation of any of the three residues involved centrally in the hydrophobic core of each primary dimer either prevented folding (Ile332) or allowed folding only at high protein concentration or low temperature (Leu330 and Phe341). Nine hydrophobic residues per monomer constitute critical determinants for the stability and oligomerization status of this p53 domain.     

pH-dependent interactions and the stability and folding kinetics of the N-terminal domain of L9. Electrostatic interactions are only weakly formed in the transition state for folding

The role of electrostatic interactions in the stability and the folding of the N-terminal domain of the ribosomal protein L9 (NTL9) was investigated by determining the effects of varying the pH conditions. Urea denaturations and thermal unfolding experiments were used to measure the free energy of folding, DeltaG degrees, at 18 different pH values, ranging from pH 1.1 to pH 10.5. Folding rates were measured at 19 pH values between pH 2.1 and pH 9.5, and unfolding rates were determined at 15 pH values in this range using stopped-flow fluorescence experiments. The protein is maximally stable between pH 5.5 and 7.5 with a value of DeltaG degrees =4.45 kcal mol(-1). The folding rate reaches a maximum at pH 5.5, however the change in folding rates with pH is relatively modest. Over the pH range of 2.1 to 5.5 there is a small increase in folding rates, ln (k(f)) changes from 5.1 to 6.8. However, the change in stability is more dramatic, with a difference of 2.6 kcal mol(-1) between pH 2.0 and pH 5.4. The change in stability is largely due to the smaller barrier for unfolding at low pH values. The natural log of the unfolding rates varies by approximately four units between pH 2.1 and pH 5.5. The stability of the protein decreases above pH 7.5 and again the change is largely due to changes in the unfolding rate. ln (k(f)) varies by less than one unit between pH 5.5 and pH 9.5 while DeltaG degrees decreases by 2.4 kcal mol(-1) over the range of pH 5. 4 to pH 10.0, which corresponds to a change in ln K(eq) of 4.0. These studies show that pH-dependent interactions contribute significantly to the overall stability of the protein but have only a small effect upon the folding kinetics, indicating that electrostatic interactions are weakly formed in the transition state for folding.     

Sequence dependence for the energetics of dangling ends and terminal base pairs in ribonucleic acid

Stability increments of terminal unpaired nucleotides (dangling ends) and terminal base pairs on the core helixes AUGCAU and UGCGCA are reported. Enthalpy, entropy, and free energy changes of helix formation were measured spectrophotometrically for 18 oligoribonucleotides containing the core sequences. The results indicate 3' dangling purines add more stability than 3' dangling pyrimidines. In most cases, the additional stability from a 3' dangling end on an AU base pair is less than that on a GC base pair [Freier, S.M., Burger, B.J., Alkema, D., Neilson, T., & Turner, D.H. (1985) Biochemistry 22, 6198-6206]. The sequence dependence provides a test for the importance of dangling ends for various RNA interactions. Correlations are suggested with codon context effects and with the three-dimensional structure of yeast phenylalanine transfer RNA. In the latter case, all terminal unpaired nucleotides having stability increments more favorable than -1 kcal/mol are stacked on the adjacent base pair. All terminal unpaired nucleotides having stability increments less favorable than -0.3 kcal/mol are not stacked on the adjacent base pair. In several cases, this lack of stacking is associated with a turn in the sugar-phosphate backbone. This suggests stability increments measured on oligoribonucleotides may be useful for predicting tertiary structure in large RNA molecules. Comparison of the stability increments for terminal dangling ends and base pairs, and of terminal GC and AU base pairs, indicates the free energy increment associated with forming a hydrogen bond can be about -1 kcal/mol of hydrogen bond.