Urea destabilizes RNA by forming stacking interactions and multiple hydrogen bonds with nucleic acid bases

Single-molecule fluorescence resonance energy transfer studies of the human telomerase RNA pseudoknot: temperature-/urea-dependent folding kinetics and thermodynamics

The ribonucleoprotein telomerase is an RNA-dependent DNA polymerase that catalyzes the repetitive addition of a short, species-specific, DNA sequence to the ends of linear eukaryotic chromosomes. The single RNA component of telomerase contains both the template sequence for DNA synthesis and a functionally critical pseudoknot motif, which can also exist as a less stable hairpin. Here we use a minimal version of the human telomerase RNA pseudoknot to study this hairpin-pseudoknot structural equilibrium using temperature-controlled single-molecule fluorescence resonance energy transfer (smFRET) experiments. The urea dependence of these experiments aids in determination of the folding kinetics and thermodynamics. The wild-type pseudoknot behavior is compared and contrasted to a mutant pseudoknot sequence implicated in a genetic disorder-dyskeratosis congenita. These findings clearly identify that this 2nt noncomplementary mutation destabilizes the folding of the wild-type pseudoknot by substantially reducing the folding rate constant (≈ 400-fold) while only nominally increasing the unfolding rate constant (≈ 5-fold). Furthermore, the urea dependence of the equilibrium and rate constants is used to develop a free energy landscape for this unimolecular equilibrium and propose details about the structure of the transition state. Finally, the urea-dependent folding experiments provide valuable physical insights into the mechanism for destabilization of RNA pseudoknots by such chemical denaturants.

Molecular crowding favors reactivity of a human ribozyme under physiological ionic conditions

In an effort to relate RNA folding to function under cellular-like conditions, we monitored the self-cleavage reaction of the human hepatitis delta virus-like CPEB3 ribozyme in the background of physiological ionic concentrations and various crowding and cosolute agents. We found that at physiological free Mg(2+) concentrations (∼0.1-0.5 mM), both crowders and cosolutes stimulate the rate of self-cleavage, up to ∼6-fold, but that in 10 mM Mg(2+) (conditions widely used for in vitro ribozyme studies) these same additives have virtually no effect on the self-cleavage rate. We further observe a dependence of the self-cleavage rate on crowder size, wherein the level of rate stimulation is diminished for crowders larger than the size of the unfolded RNA. Monitoring effects of crowding and cosolute agents on rates in biological amounts of urea revealed additive-promoted increases at both low and high Mg(2+) concentrations, with a maximal stimulation of more than 10-fold and a rescue of the rate to its urea-free values. Small-angle X-ray scattering experiments reveal a structural basis for this stimulation in that higher-molecular weight crowding agents favor a more compact form of the ribozyme in 0.5 mM Mg(2+) that is essentially equivalent to the form under standard ribozyme conditions of 10 mM Mg(2+) without a crowder. This finding suggests that at least a portion of the rate enhancement arises from favoring the native RNA tertiary structure. We conclude that cellular-like crowding supports ribozyme reactivity by favoring a compact form of the ribozyme, but only under physiological ionic and cosolute conditions.

Self-assembly of nonionic surfactant Tween 20@2β-CD inclusion complexes in dilute solution

It has long been considered that the addition of cyclodextrins (CDs) disfavors the self-assembly of surfactants in dilute solutions since the hydrophobic effect is destroyed upon the formation of the hydrophiphilic CD/surfactant inclusion complex. However, in this work, we found that β-CD/nonionic surfactant inclusion complexes are able to self-assemble into vesicles in dilute solutions, namely in solutions with concentration lower than the CMC of surfactants. When using Tween 20 as a model surfactant, HNMR and MS measurements indicate that the building block for the vesicles is the channel type Tween 20@2β-CD inclusion complex. Structure and IR analysis suggests that the self-assembly of hydrophilic Tween 20@2β-CD is driven by H-bonds between both the headgroup of Tween 20 and the hydroxyl groups of β-CD. The self-assembly of the inclusion complex between the β-CD and the nonionic surfactant in dilute solution is found to be a general phenomenon. Undoubtedly, surfactant@2β-CD inclusion complex can be a novel building block for nonamphiphilic self-assembly, which provides a new physical insight for the influence of cyclodextrins on the self-assembly of surfactants.

Urea-induced denaturation of preQ1-riboswitch

Urea, a polar molecule with a large dipole moment, not only destabilizes folded RNA structures but can also enhance the folding rates of large ribozymes. Unlike the mechanism of urea-induced unfolding of proteins, which is well understood, the action of urea on RNA has barely been explored. We performed extensive all-atom molecular dynamics simulations to determine the molecular underpinnings of urea-induced RNA denaturation. Urea displays its denaturing power in both secondary and tertiary motifs of the riboswitch structure. Our simulations reveal that the denaturation of RNA structures is mainly driven by the hydrogen-bonding and stacking interactions of urea with the bases. Through detailed studies of the simulation trajectories, we found that geminate pairs between urea and bases due to hydrogen bonds and stacks persist only ~0.1-1 ns, which suggests that the urea-base interaction is highly dynamic. Most importantly, the early stage of base-pair disruption is triggered by penetration of water molecules into the hydrophobic domain between the RNA bases. The infiltration of water into the narrow space between base pairs is critical in increasing the accessibility of urea to transiently disrupted bases, thus allowing urea to displace inter-base hydrogen bonds. This mechanism--water-induced disruption of base pairs resulting in the formation of a "wet" destabilized RNA followed by solvation by urea--is the exact opposite of the two-stage denaturation of proteins by urea. In the latter case, initial urea penetration creates a dry globule, which is subsequently solvated by water, leading to global protein unfolding. Our work shows that the ability to interact with both water and polar or nonpolar components of nucleotides makes urea a powerful chemical denaturant for nucleic acids.

Translocating kilobase RNA through the Staphylococcal α-hemolysin nanopore

The electrophoretic translocation of polynucleotides through nanopores may permit direct single-molecule nucleic acid sequencing. Here we describe the translocation of ssRNA heteropolymers (91-6083 bases) through the α-hemolysin nanopore. Translocation of these long ssRNAs is characterized by surprisingly long, almost complete ionic current blockades with durations averaging milliseconds per base (at +180 mV). The event durations decrease exponentially with increased transmembrane potential but are largely unaffected by the presence of urea. When the ssRNA is coupled at the 3' end to streptavidin, which cannot translocate through the pore, permanent blockades are observed, supporting our conclusion that the transient blockades arise from ssRNA translocation.

Quantifying functional group interactions that determine urea effects on nucleic acid helix formation

Urea destabilizes helical and folded conformations of nucleic acids and proteins, as well as protein-nucleic acid complexes. To understand these effects, extend previous characterizations of interactions of urea with protein functional groups, and thereby develop urea as a probe of conformational changes in protein and nucleic acid processes, we obtain chemical potential derivatives (μ23 = dμ2/dm3) quantifying interactions of urea (component 3) with nucleic acid bases, base analogues, nucleosides, and nucleotide monophosphates (component 2) using osmometry and hexanol-water distribution assays. Dissection of these μ23 values yields interaction potentials quantifying interactions of urea with unit surface areas of nucleic acid functional groups (heterocyclic aromatic ring, ring methyl, carbonyl and phosphate O, amino N, sugar (C and O); urea interacts favorably with all these groups, relative to interactions with water. Interactions of urea with heterocyclic aromatic rings and attached methyl groups (as on thymine) are particularly favorable, as previously observed for urea-homocyclic aromatic ring interactions. Urea m-values determined for double helix formation by DNA dodecamers near 25 °C are in the range of 0.72-0.85 kcal mol(-1)m(-1) and exhibit little systematic dependence on nucleobase composition (17-42% GC). Interpretation of these results using the urea interaction potentials indicates that extensive (60-90%) stacking of nucleobases in the separated strands in the transition region is required to explain the m-value. Results for RNA and DNA dodecamers obtained at higher temperatures, and literature data, are consistent with this conclusion. This demonstrates the utility of urea as a quantitative probe of changes in surface area (ΔASA) in nucleic acid processes.

Single-molecule chemical denaturation of riboswitches

To date, single-molecule RNA science has been developed almost exclusively around the effect of metal ions as folding promoters and stabilizers of the RNA structure. Here, we introduce a novel strategy that combines single-molecule Förster resonance energy transfer (FRET) and chemical denaturation to observe and manipulate RNA dynamics. We demonstrate that the competing interplay between metal ions and denaturant agents provides a platform to extract information that otherwise will remain hidden with current methods. Using the adenine-sensing riboswitch aptamer as a model, we provide strong evidence for a rate-limiting folding step of the aptamer domain being modulated through ligand binding, a feature that is important for regulation of the controlled gene. In the absence of ligand, the rate-determining step is dominated by the formation of long-range key tertiary contacts between peripheral stem-loop elements. In contrast, when the adenine ligand interacts with partially folded messenger RNAs, the aptamer requires specifically bound Mg²⁺ ions, as those observed in the crystal structure, to progress further towards the native form. Moreover, despite that the ligand-free and ligand-bound states are indistinguishable by FRET, their different stability against urea-induced denaturation allowed us to discriminate them, even when they coexist within a single FRET trajectory; a feature not accessible by existing methods.

Introductory lecture: interpreting and predicting Hofmeister salt ion and solute effects on biopolymer and model processes using the solute partitioning model

Understanding how Hofmeister salt ions and other solutes interact with proteins, nucleic acids, other biopolymers and water and thereby affect protein and nucleic acid processes as well as model processes (e.g. solubility of model compounds) in aqueous solution is a longstanding goal of biophysical research. Empirical Hofmeister salt and solute "m-values" (derivatives of the observed standard free energy change for a model or biopolymer process with respect to solute or salt concentration m3) are equal to differences in chemical potential derivatives: m-value = delta(dmu2/dm3) = delta mu23, which quantify the preferential interactions of the solute or salt with the surface of the biopolymer or model system (component 2) exposed or buried in the process. Using the solute partitioning model (SPM), we dissect mu23 values for interactions of a solute or Hofmeister salt with a set of model compounds displaying the key functional groups of biopolymers to obtain interaction potentials (called alpha-values) that quantify the interaction of the solute or salt per unit area of each functional group or type of surface. Interpreted using the SPM, these alpha-values provide quantitative information about both the hydration of functional groups and the competitive interaction of water and the solute or salt with functional groups. The analysis corroborates and quantifies previous proposals that the Hofmeister anion and cation series for biopolymer processes are determined by ion-specific, mostly unfavorable interactions with hydrocarbon surfaces; the balance between these unfavorable nonpolar interactions and often-favorable interactions of ions with polar functional groups determine the series null points. The placement of urea and glycine betaine (GB) at opposite ends of the corresponding series of nonelectrolytes results from the favorable interactions of urea, and unfavorable interactions of GB, with many (but not all) biopolymer functional groups. Interaction potentials and local-bulk partition coefficients quantifying the distribution of solutes (e.g. urea, glycine betaine) and Hofmeister salt ions in the vicinity of each functional group make good chemical sense when interpreted in terms of competitive noncovalent interactions. These interaction potentials allow solute and Hofmeister (noncoulombic) salt effects on protein and nucleic acid processes to be interpreted or predicted, and allow the use of solutes and salts as probes of

Denaturation of RNA secondary and tertiary structure by urea: simple unfolded state models and free energy parameters account for measured m-values

To investigate the mechanism by which urea destabilizes RNA structure, urea-induced unfolding of four different RNA secondary and tertiary structures was quantified in terms of an m-value, the rate at which the free energy of unfolding changes with urea molality. From literature data and our osmometric study of a backbone analogue, we derived average interaction potentials (per square angstrom of solvent accessible surface) between urea and three kinds of RNA surfaces: phosphate, ribose, and base. Estimates of the increases in solvent accessible surface areas upon RNA denaturation were based on a simple model of unfolded RNA as a combination of helical and single-strand segments. These estimates, combined with the three interaction potentials and a term to account for interactions of urea with released ions, yield calculated m-values that are in good agreement with experimental values (200 mm monovalent salt). Agreement was obtained only if single-stranded RNAs were modeled in a highly stacked, A-form conformation. The primary driving force for urea-induced denaturation is the strong interaction of urea with the large surface areas of bases that become exposed upon denaturation of either RNA secondary or tertiary structure, though interactions of urea with backbone and released ions may account for up to a third of the m-value. Urea m-values for all four RNAs are salt-dependent, which we attribute to an increased extension (or decreased charge density) of unfolded RNAs with an increased urea concentration. The sensitivity of the urea m-value to base surface exposure makes it a potentially useful probe of the conformations of RNA unfolded states.

Balancing target flexibility and target denaturation in computational fragment-based inhibitor discovery

Accounting for target flexibility and selecting "hot spots" most likely to be able to bind an inhibitor continue to be challenges in the field of structure-based drug design, especially in the case of protein-protein interactions. Computational fragment-based approaches using molecular dynamics (MD) simulations are a promising emerging technology having the potential to address both of these challenges. However, the optimal MD conditions permitting sufficient target flexibility while also avoiding fragment-induced target denaturation remain ambiguous. Using one such technology (Site Identification by Ligand Competitive Saturation, SILCS), conditions were identified to either prevent denaturation or identify and exclude trajectories in which subtle but important denaturation was occurring. The target system used was the well-characterized protein cytokine IL-2, which is involved in a protein-protein interface and, in its unliganded crystallographic form, lacks surface pockets that can serve as small-molecule binding sites. Nonetheless, small-molecule inhibitors have previously been discovered that bind to two "cryptic" binding sites that emerge only in the presence of ligand binding, highlighting the important role of IL-2 flexibility. Using the above conditions, SILCS with hydrophobic fragments was able to identify both sites based on favorable fragment binding while avoiding IL-2 denaturation. An important additional finding was that acetonitrile, a water-miscible fragment, fails to identify either site yet can induce target denaturation, highlighting the importance of fragment choice.

Exploring the strength, mode, dynamics, and kinetics of binding interaction of a cationic biological photosensitizer with DNA: implication on dissociation of the drug-DNA complex via detergent sequestration

The present study aims at exploring a detailed characterization of the binding interaction of a promising cancer cell photosensitizer, harmane (HM), with DNA extracted from herring sperm. The polarity-sensitive prototropic transformation of HM, a naturally occurring, fluorescent, drug-binding alkaloid, β-carboline, is remarkably modified upon interaction with DNA and is manifested through significant modulations on the absorption and emission profiles of HM. From the series of studies undertaken in the present program, for example, absorption; steady-state emission; the effect of chaotrope (urea); iodide ion-induced steady-state fluorescence quenching; circular dichroism (CD); and helix melting from absorption spectroscopy; the mode of binding of HM into the DNA helix has been substantiated to be principally intercalative. Concomitantly, a discernible dependence of the photophysics of the DNA-bound drug on the medium ionic strength indicates that electrostatic attraction should not be ignored in the interaction. Efforts have also been delivered to delineate the dynamical aspects of the interaction, such as modulation in time-resolved fluorescence decay and rotational relaxation dynamics of the drug within the DNA environment. In view of the prospective biological applications of HM, the issue of facile dissociation of intercalated HM from the DNA helix also comprises a crucial prerequisite for the functioning as an effective therapeutic agent. In this context, our results imply that the concept of detergent-sequestered dissociation of the drug from the drug-DNA complex can be a prospective strategy through an appropriate choice of the detergent molecule. The utility of the present work resides in exploring the potential applicability of the fluorescence property of HM for studying its interactions with a relevant biological target, for example, DNA. In addition, the methods and techniques used in the present work can also be exploited to study the interaction of HM with other biological, biomimicking assemblies and drug delivery vehicles, and so forth.

Effects of a protecting osmolyte on the ion atmosphere surrounding DNA duplexes

Osmolytes are small, chemically diverse, organic solutes that function as an essential component of cellular stress response. Protecting osmolytes enhance protein stability via preferential exclusion, and nonprotecting osmolytes, such as urea, destabilize protein structures. Although much is known about osmolyte effects on proteins, less is understood about osmolyte effects on nucleic acids and their counterion atmospheres. Nonprotecting osmolytes destabilize nucleic acid structures, but effects of protecting osmolytes depend on numerous factors including the type of nucleic acid and the complexity of the functional fold. To begin quantifying protecting osmolyte effects on nucleic acid interactions, we used small-angle X-ray scattering (SAXS) techniques to monitor DNA duplexes in the presence of sucrose. This protecting osmolyte is a commonly used contrast matching agent in SAXS studies of protein-nucleic acid complexes; thus, it is important to characterize interaction changes induced by sucrose. Measurements of interactions between duplexes showed no dependence on the presence of up to 30% sucrose, except under high Mg(2+) conditions where stacking interactions were disfavored. The number of excess ions associated with DNA duplexes, reported by anomalous small-angle X-ray scattering (ASAXS) experiments, was sucrose independent. Although protecting osmolytes can destabilize secondary structures, our results suggest that ion atmospheres of individual duplexes remain unperturbed by sucrose.

Equilibrium study of protein denaturation by urea

Though urea is commonly used to denature proteins, the molecular mechanism of its denaturing ability is still a subject of considerable debate. Previous molecular dynamics simulation studies have sought to elucidate the mechanism of urea denaturation by focusing on the pathway of denaturation rather than examining the effect of urea on the folding/unfolding equilibrium, which is commonly measured in experiment. Here we report the reversible folding/unfolding equilibrium of Trp-cage miniprotein in the presence of urea, over a broad range of urea concentrations, using all-atom Replica exchange MD simulations. The simulations capture the experimentally observed linear dependence of unfolding free energy on urea concentration. We find that the denaturation is driven by favorable direct interaction of urea with the protein through both electrostatic and van der Waals forces and quantify their contribution. Though the magnitude of direct electrostatic interaction of urea is larger than van der Waals, the difference between unfolded and folded ensembles is dominated by the van der Waals interaction. We also find that hydrogen bonding of urea to the peptide backbone does not play a dominant role in denaturation. The unfolded ensemble sampled depends on urea concentration, with greater urea concentration favoring conformations with greater solvent exposure. The m-value is predicted to increase with temperature and more strongly so with pressure.