Direct observation of the temperature-induced melting process of the Salmonella fourU RNA thermometer at base-pair resolution

In prokaryotes, RNA thermometers regulate a number of heat shock and virulence genes. These temperature sensitive RNA elements are usually located in the 5′-untranslated regions of the regulated genes. They repress translation initiation by base pairing to the Shine–Dalgarno sequence at low temperatures. We investigated the thermodynamic stability of the temperature labile hairpin 2 of the Salmonella fourU RNA thermometer over a broad temperature range and determined free energy, enthalpy and entropy values for the base-pair opening of individual nucleobases by measuring the temperature dependence of the imino proton exchange rates via NMR spectroscopy. Exchange rates were analyzed for the wild-type (wt) RNA and the A8C mutant. The wt RNA was found to be stabilized by the extraordinarily stable G14–C25 base pair. The mismatch base pair in the wt RNA thermometer (A8–G31) is responsible for the smaller cooperativity of the unfolding transition in the wt RNA. Enthalpy and entropy values for the base-pair opening events exhibit linear correlation for both RNAs. The slopes of these correlations coincide with the melting points of the RNAs determined by CD spectroscopy. RNA unfolding occurs at a temperature where all nucleobases have equal thermodynamic stabilities. Our results are in agreement with a consecutive zipper-type unfolding mechanism in which the stacking interaction is responsible for the observed cooperativity. Furthermore, remote effects of the A8C mutation affecting the stability of nucleobase G14 could be identified. According to our analysis we deduce that this effect is most probably transduced via the hydration shell of the RNA.     

Thermodynamics of RNA melting,one base pair at a time

The melting of base pairs is a ubiquitous feature of RNA structural transitions, which are widely used to sense and respond to cellular stimuli. A recent study employing solution nuclear magnetic resonance (NMR) imino proton exchange spectroscopy provides a rare base-pair-specific view of duplex melting in the Salmonella FourU RNA thermosensor, which regulates gene expression in response to changes in temperature at the translational level by undergoing a melting transition. The authors observe “microscopic” enthalpy–entropy compensation—often seen “macroscopically” across a series of related molecular species—across base pairs within the same RNA. This yields variations in base-pair stabilities that are an order of magnitude smaller than corresponding variations in enthalpy and entropy. A surprising yet convincing link is established between the slopes of enthalpy–entropy correlations and RNA melting points determined by circular dichroism (CD), which argues that unfolding occurs when base-pair stabilities are equalized. A single AG-to-CG mutation, which enhances the macroscopic hairpin thermostability and folding cooperativity and renders the RNA thermometer inactive in vivo, spreads its effect microscopically throughout all base pairs in the RNA, including ones far removed from the site of mutation. The authors suggest that an extended network of hydration underlies this long-range communication. This study suggests that the deconstruction of macroscopic RNA unfolding in terms of microscopic unfolding events will require careful consideration of water interactions.     

Use of ultra stable UNCG tetraloop hairpins to fold RNA structures: thermodynamic and spectroscopic applications.

RNA molecules of > 20 nucleotides have been the focus of numerous recent NMR structural studies. Several investigators have used the UNCG family of hairpins to ensure proper folding. We show that th UUCG hairpin has a minimum requirement of a two base-pair stem. Hairpins with a CG loop closing base pair and an initial 5'CG or 5'GC base pair have a melting temperature approximately 55 degrees C in 10 mM sodium phosphate. The high stability of even such small hairpins suggests that the hairpin can serve as a nucleation site for folding. For high resolution NMR work, the UNCG loop family (UACG in particular) provides excellent spectroscopic markers in one-dimensional exchangeable spectra, in two-dimensional COSY spectra and in NOESY spectra that clearly define it as forming a hairpin. This allows straightforward initiation of chemical shift assignments.     

Bacterial RNA thermometers: molecular zippers and switches

Bacteria use complex strategies to coordinate temperature-dependent gene expression. Many genes encoding heat shock proteins and virulence factors are regulated by temperature-sensing RNA sequences, known as RNA thermometers (RNATs), in their mRNAs. For these genes, the 5' untranslated region of the mRNA folds into a structure that blocks ribosome access at low temperatures. Increasing the temperature gradually shifts the equilibrium between the closed and open conformations towards the open structure in a zipper-like manner, thereby increasing the efficiency of translation initiation. Here, we review the known molecular principles of RNAT action and the hierarchical RNAT cascade in Escherichia coli. We also discuss RNA-based thermosensors located upstream of cold shock and other genes, translation of which preferentially occurs at low temperatures and which thus operate through a different, more switch-like mechanism. Finally, we consider the potential biotechnological applications of natural and synthetic RNATs.     

Integration Preferences of Wildtype AAV-2 for Consensus Rep-Binding Sites at Numerous Loci in the Human Genome

Adeno-associated virus type 2 (AAV) is known to establish latency by preferential integration in human chromosome 19q13.42. The AAV non-structural protein Rep appears to target a site called AAVS1 by simultaneously binding to Rep-binding sites (RBS) present on the AAV genome and within AAVS1. In the absence of Rep, as is the case with AAV vectors, chromosomal integration is rare and random. For a genome-wide survey of wildtype AAV integration a linker-selection-mediated (LSM)-PCR strategy was designed to retrieve AAV-chromosomal junctions. DNA sequence determination revealed wildtype AAV integration sites scattered over the entire human genome. The bioinformatic analysis of these integration sites compared to those of rep-deficient AAV vectors revealed a highly significant overrepresentation of integration events near to consensus RBS. Integration hotspots included AAVS1 with 10% of total events. Novel hotspots near consensus RBS were identified on chromosome 5p13.3 denoted AAVS2 and on chromsome 3p24.3 denoted AAVS3. AAVS2 displayed seven independent junctions clustered within only 14 bp of a consensus RBS which proved to bind Rep in vitro similar to the RBS in AAVS3. Expression of Rep in the presence of rep-deficient AAV vectors shifted targeting preferences from random integration back to the neighbourhood of consensus RBS at hotspots and numerous additional sites in the human genome. In summary, targeted AAV integration is not as specific for AAVS1 as previously assumed. Rather, Rep targets AAV to integrate into open chromatin regions in the reach of various, consensus RBS homologues in the human genome.     

NMR structure of stem-loop SL2 of the HIV-1 psi RNA packaging signal reveals a novel A-U-A base-triple platform

The genome of the human immunodeficiency virus type-1 (HIV-1) contains a stretch of approximately 120 nucleotides known as the psi-site that is essential for RNA packaging during virus assembly. These nucleotides have been proposed to form four stem-loops (SL1-SL4) that have both independent and overlapping functions. Stem-loop SL2 is important for efficient recognition and packaging of the full-length, unspliced viral genome, and also contains the major splice-donor site (SD) for mRNA splicing. We have determined the structure of the 19-residue SL2 oligoribonucleotide by heteronuclear NMR methods. The structure is generally consistent with the most recent of two earlier secondary structure predictions, with residues G1-G2-C3-G4 and C6-U7 forming standard Watson Crick base-pairs with self-complementary residues C16-G17-C18-C19 and A12-G13, respectively. However, residue A15, which is located near the center of the stem, does not form a predicted bulge, and residues A5 and U14 do not form an expected Watson-Crick base-pair. Instead, these residues form a novel A5-U14-A15 base-triple that appears to be stabilized by hydrogen bonds from A15-H61 and -H62 to A5-N1 and U14-O2, respectively; from A5-H61 to U14-O2, and from C16-H42 to U14-O2'. A kink in the backbone allows the aromatic rings of the sequential U14-A15 residues to be approximately co-planar, adopting a stable "platform motif" that is structurally similar to the A-A (adenosine) platforms observed in the P4-P6 ribozyme domain of the Tetrahymena group I intron. Platform motifs generally function in RNA by mediating long-range interactions, and it is therefore possible that the A-U-A base-triple platform mediates long-range interactions that either stabilize the psi-RNA or facilitate splicing and/or packaging. Residue G8 of the G8-G9-U10-G11 tetraloop is stacked above the U7-A12 base-pair, and the remaining tetraloop residues are disordered and available for potential interactions with either other RNA or protein components.     

A nonlinear dynamic model of DNA with a sequence-dependent stacking term

No simple model exists that accurately describes the melting behavior and breathing dynamics of double-stranded DNA as a function of nucleotide sequence. This is especially true for homogenous and periodic DNA sequences, which exhibit large deviations in melting temperature from predictions made by additive thermodynamic contributions. Currently, no method exists for analysis of the DNA breathing dynamics of repeats and of highly G/C- or A/T-rich regions, even though such sequences are widespread in vertebrate genomes. Here, we extend the nonlinear Peyrard–Bishop–Dauxois (PBD) model of DNA to include a sequence-dependent stacking term, resulting in a model that can accurately describe the melting behavior of homogenous and periodic sequences. We collect melting data for several DNA oligos, and apply Monte Carlo simulations to establish force constants for the 10 dinucleotide steps (CG, CA, GC, AT, AG, AA, AC, TA, GG, TC). The experiments and numerical simulations confirm that the GG/CC dinucleotide stacking is remarkably unstable, compared with the stacking in GC/CG and CG/GC dinucleotide steps. The extended PBD model will facilitate thermodynamic and dynamic simulations of important genomic regions such as CpG islands and disease-related repeats.     

Thermal unfolding of small proteins with SH3 domain folding pattern

The thermal unfolding of three SH3 domains of the Tec family of tyrosine kinases was studied by differential scanning calorimetry and CD spectroscopy. The unfolding transition of the three protein domains in the acidic pH region can be described as a reversible two-state process. For all three SH3 domains maximum stability was observed in the pH region 4.5 < pH < 7.0 where these domains unfold at temperatures of 353K (Btk), 342K (Itk), and 344K (Tec). At these temperatures an enthalpy change of 196 kJ/mol, 178 kJ/mol, and 169 kJ/mol was measured for Btk-, Itk-, and Tec-SH3 domains, respectively. The determined changes in heat capacity between the native and the denatured state are in an usual range expected for small proteins. Our analysis revealed that all SH3 domains studied are only weakly stabilized and have free energies of unfolding which do not exceed 12–16 kJ/mol but show quite high melting temperatures. Comparing unfolding free energies measured for eukaryotic SH3 domains with those of the topologically identical Sso7d protein from the hyperthermophile Sulfolobus solfataricus, the increased melting temperature of the thermostable protein is due to a broadening as well as a significant lifting of its stability curve. However, at their physiological temperatures, 310K for mesophilic SH3 domains and 350K for Sso7d, eukaryotic SH3 domains and Sso7d show very similar stabilities. Proteins 31:309–319, 1998. © 1998 Wiley-Liss, Inc.     

Translational Entropy and DNA Duplex Stability

Investigation of folding/unfolding DNA duplexes of various size and composition by superprecise calorimetry has revised several long-held beliefs concerning the forces responsible for the formation of the double helix. It was established that: 1) the enthalpy and the entropy of duplex unfolding are temperature dependent, increasing with temperature rise and having the same heat capacity increment for CG and AT pairs; 2) the enthalpy of AT melting is greater than that of the CG pair, so the stabilizing effect of the CG pair in comparison with AT results not from its larger enthalpic contribution (as expected from its extra hydrogen bond), but from the larger entropic contribution of the AT pair that results from its ability to fix ordered water in the minor groove and release it upon duplex unfolding; 3) the translation entropy, resulting from the appearance of a new kinetic unit on duplex dissociation, determines the dependence of duplex stability on its length and its concentration (it is an order-of-magnitude smaller than predicted from the statistical mechanics of gases and is fully expressed by the stoichiometric correction term); 4) changes in duplex stability on reshuffling the sequence (the “nearest-neighbor effect”) result from the immobilized water molecules fixed by AT pairs in the minor groove; and 5) the evaluated thermodynamic components permit a quantitative expression of DNA duplex stability.     

Rabbit liver tRNAVal1: II. Unusual secondary structure of T ψ C stem and loop due to a U54:A60 base pair

In contrast to all other known tRNAs, mammalian tRNA^Val₁ contains two adenosines A₅₉ and A₆₀, opposite to U₅₄ and ψ₅₅ in the UψCG sequence of the TψC loop, which could form unusual A:U (or A:ψ) pairs in addition to the five “normal” G:C pairs. In order to measure the number of G:C and A:U (A:ψ) pairs in the TψC stem, we prepared the 30 nucleotide long 3′-terminal fragment of this tRNA by “m⁷G-cleavage”. From differentiated melting curves and temperature jump experiments it was concluded that the TψC stem in this fragment is in fact extended by an additional A₆₀:U₅₄ pair. A dimer of this fragment with 14 base pairs was characterized by gel electrophoresis and by the same physical methods. An additional A:U pair in the tRNA^Val₁ fragment does not necessarily mean that this is also true for intact tRNA. However, we showed that U₅₄ is far less available for enzymatic methylation in mammalian tRNA^Val₁ compared to tRNA from T⁻E. coli. This clear difference in U₅₄ reactivity, together with the identification of an extra A₆₀:U₅₄ pair in the UψCG containing fragment suggests the presence of a 6 base pair TψC stem and a 5 nucleotide TψC loop in this tRNA.     

Thermodynamic differences among homologous thermophilic and mesophilic proteins.

Here, we analyze the thermodynamic parameters and their correlations in families containing homologous thermophilic and mesophilic proteins which show reversible two-state folding <--> unfolding transitions between the native and the denatured states. For the proteins in these families, the melting temperatures correlate with the maximal protein stability change (between the native and the denatured states) as well as with the enthalpic and entropic changes at the melting temperature. In contrast, the heat capacity change is uncorrelated with the melting temperature. These and additional results illustrate that higher melting temperatures are largely obtained via an upshift and broadening of the protein stability curves. Both thermophilic and mesophilic proteins are maximally stable around room temperature. However, the maximal stabilities of thermophilic proteins are considerably greater than those of their mesophilic homologues. At the living temperatures of their respective source organisms, homologous thermophilic and mesophilic proteins have similar stabilities. The protein stability at the living temperature of the source organism does not correlate with the living temperature of the protein. We tie thermodynamic observations to microscopics via the hydrophobic effect and a two-state model of the water structure. We conclude that, to achieve higher stability and greater resistance to high and low temperatures, specific interactions, particularly electrostatic, should be engineered into the protein. The effect of these specific interactions is largely reflected in an increased enthalpy change at the melting temperature.     

Why water-soluble, compact, globular proteins have similar specific enthalpies of unfolding at 110 degrees C.

The changes in free energy, enthalpy, and entropy of unfolding have been measured for many water-soluble, compact, globular proteins by a number of workers. In principle, a wide range in stability could be achieved by proteins, as measured by the free energy of unfolding; in practice, evolution only allows a narrow range in this quantity. Proteins are only marginally stable at room temperature for many possible reasons, including ensuring that folding is reversible and polypeptide chains are not trapped in incorrectly folded structures. Many of these proteins have approximately the same values of enthalpy of unfolding around 110 degrees C. We show here that this arises because the change in entropy of unfolding at room temperature and the change in heat capacity on unfolding, which governs the temperature variation of the enthalpy and entropy, both vary with the magnitude of the hydrophobic effect in the protein. As all these proteins have evolved to achieve similar stabilities at room temperature, the enthalpy of unfolding will also vary with the size of the hydrophobic effect in the protein. A consequence of this is that curves of the specific unfolding enthalpy against temperature for different proteins intersect around 110 degrees C. A similar conclusion, on the basis of similar melting points rather than similar free energies of unfolding, has been reached independently by Baldwin and Muller (R. L. Baldwin, personal communication).     

Cytarabine-Induced Destabilization and Bending of a Model Okazaki Fragment

Abstract

The effects of cytarabine on the structural and thermodynamic properties of an Okazaki fragment were investigated using UV hyperchromicity and 2D ¹H NMR. Cytarabine significantly decreased the stability of this model Okazaki fragment, decreasing the melting temperature from 46.8 °C to 42.4 °C at 1.33 × 10^?5M. Cytarabine also markedly increased the bend angle of the Okazaki fragment duplex from 20° to 42°. Changes to the structures and stabilities of Okazaki fragments may cause the biological effects of cytarabine.     

Equilibrium unfolding pathway of an H-type RNA pseudoknot which promotes programmed -1 ribosomal frameshifting.

The equilibrium unfolding pathway of a 41-nucleotide frameshifting RNA pseudoknot from the gag-pro junction of mouse intracisternal A-type particles (mIAP), an endogenous retrovirus, has been determined through analysis of dual optical wavelength, equilibrium thermal melting profiles and differential scanning calorimetry. The mIAP pseudoknot is an H-type pseudoknot proposed to have structural features in common with the gag-pro frameshifting pseudoknots from simian retrovirus-1 (SRV-1) and mouse mammary tumor virus (MMTV). In particular, the mIAP pseudoknot is proposed to contain an unpaired adenosine base at the junction of the two helical stems (A15), as well as one in the middle of stem 2 (A35). A mutational analysis of stem 1 hairpins and compensatory base-pair substitutions incorporated into helical stem 2 was used to assign optical melting transitions to molecular unfolding events. The optical melting profile of the wild-type RNA is most simply described by four sequential two-state unfolding transitions. Stem 2 melts first in two closely coupled low-enthalpy transitions at low tmin which the stem 3' to A35, unfolds first, followed by unfolding of the remainder of the helical stem. The third unfolding transition is associated with some type of stacking interactions in the stem 1 hairpin loop not present in the pseudoknot. The fourth transition is assigned to unfolding of stem 1. In all RNAs investigated, DeltaHvH approximately DeltaHcal, suggesting that DeltaCpfor unfolding is small. A35 has the thermodynamic properties expected for an extrahelical, unpaired nucleotide. Deletion of A15 destabilizes the stem 2 unfolding transition in the context of both the wild-type and DeltaA35 mutant RNAs only slightly, by DeltaDeltaG degrees approximately 1 kcal mol-1(at 37 degrees C). The DeltaA15 RNA is considerably more susceptible to thermal denaturation in the presence of moderate urea concentrations than is the wild-type RNA, further evidence of a detectable global destabilization of the molecule. Interestingly, substitution of the nine loop 2 nucleotides with uridine residues induces a more pronounced destabilization of the molecule (DeltaDeltaG degrees approximately 2.0 kcal mol-1), a long-range, non-nearest neighbor effect. These findings provide the thermodynamic basis with which to further refine the relationship between efficient ribosomal frameshifting and pseudoknot structure and stability.