Dynamics of the IRE RNA hairpin loop probed by 2-aminopurine fluorescence and stochastic dynamics simulations

Computational Modeling of RNA Aptamers: Structure Prediction of the Apo State

RNA aptamers are single-stranded oligonucleotides that bind to specific molecular targets with high affinity and specificity. To design aptamers for new applications, it is critical to understand the ligand binding mechanism in terms of the structure and dynamics of the ligand-bound and apo states. The problem is that most of the NMR or X-ray crystal structures available for RNA aptamers are for ligand-bound states. Available apo state structures, mostly characterized by crystallization under nonphysiological conditions or probed by low resolution techniques, might fail to represent the diverse structural variations of the apo state in solution. Here, we develop an approach to obtain a representative ensemble of apo structures that are based on in silico RNA 3D structure prediction and in vitro experiments that characterize base stacking. Using the neomycin-B aptamer as a case study, an ensemble of structures for the aptamer in the apo (unbound) state are validated and then used to investigate the ligand-binding mechanism for the aptamer in complex with neomycin-B.

Hybridization kinetics of out-of-equilibrium mixtures of short RNA oligonucleotides

Hybridization and strand displacement kinetics determine the evolution of the base paired configurations of mixtures of oligonucleotides over time. Although much attention has been focused on the thermodynamics of DNA and RNA base pairing in the scientific literature, much less work has been done on the time dependence of interactions involving multiple strands, especially in RNA. Here we provide a study of oligoribonucleotide interaction kinetics and show that it is possible to calculate the association, dissociation and strand displacement rates displayed by short oligonucleotides (5nt–12nt) that exhibit no expected secondary structure as simple functions of oligonucleotide length, CG content, ΔG of hybridization and ΔG of toehold binding. We then show that the resultant calculated kinetic parameters are consistent with the experimentally observed time dependent changes in concentrations of the different species present in mixtures of multiple competing RNA strands. We show that by changing the mixture composition, it is possible to create and tune kinetic traps that extend by orders of magnitude the typical sub-second hybridization timescale of two complementary oligonucleotides. We suggest that the slow equilibration of complex oligonucleotide mixtures may have facilitated the nonenzymatic replication of RNA during the origin of life.

The Iron Maiden. Cytosolic Aconitase/IRP1 Conformational Transition in the Regulation of Ferritin Translation and Iron Hemostasis

Maintaining iron homeostasis is fundamental for almost all living beings, and its deregulation correlates with severe and debilitating pathologies. The process is made more complicated by the omnipresence of iron and by its role as a fundamental component of a number of crucial metallo proteins. The response to modifications in the amount of the free-iron pool is performed via the inhibition of ferritin translation by sequestering consensus messenger RNA (mRNA) sequences. In turn, this is regulated by the iron-sensitive conformational equilibrium between cytosolic aconitase and IRP1, mediated by the presence of an iron-sulfur cluster. In this contribution, we analyze by full-atom molecular dynamics simulation, the factors leading to both the interaction with mRNA and the conformational transition. Furthermore, the role of the iron-sulfur cluster in driving the conformational transition is assessed by obtaining the related free energy profile via enhanced sampling molecular dynamics simulations.

Single-molecule fluorescence detection of a tricyclic nucleoside analogue

Fluorescent nucleobase surrogates capable of Watson–Crick hydrogen bonding are essential probes of nucleic acid structure and dynamics, but their limited brightness and short absorption and emission wavelengths have rendered them unsuitable for single-molecule detection. Aiming to improve on these properties, we designed a new tricyclic pyrimidine nucleoside analogue with a push–pull conjugated system and synthesized it in seven sequential steps. The resulting C-linked 8-(diethylamino)benzo[b][1,8]naphthyridin-2(1H)-one nucleoside, which we name ABN, exhibits ε₄₄₂ = 20 000 M⁻¹ cm⁻¹ and Φ_em,540 = 0.39 in water, increasing to Φ_em = 0.50–0.53 when base paired with adenine in duplex DNA oligonucleotides. Single-molecule fluorescence measurements of ABN using both one-photon and two-photon excitation demonstrate its excellent photostability and indicate that the nucleoside is present to > 95% in a bright state with count rates of at least 15 kHz per molecule. This new fluorescent nucleobase analogue, which, in duplex DNA, is the brightest and most red-shifted known, is the first to offer robust and accessible single-molecule fluorescence detection capabilities.

Fluorescent nucleoside analogue ABN is readily detected at the single-molecule level and retains a quantum yield >50% in duplex DNA oligonucleotides.

Divalent ions tune the kinetics of a bacterial GTPase center rRNA folding transition from secondary to tertiary structure

Folding of an RNA from secondary to tertiary structure often depends on divalent ions for efficient electrostatic charge screening (nonspecific association) or binding (specific association). To measure how different divalent cations modify folding kinetics of the 60 nucleotide Ecoli rRNA GTPase center, we combined stopped-flow fluorescence in the presence of Mg²⁺, Ca²⁺, or Sr²⁺ together with time-resolved small angle X-ray scattering (SAXS) in the presence of Mg²⁺ to observe the folding process. Immediately upon addition of each divalent ion, the RNA undergoes a transition from an extended state with secondary structure to a more compact structure. Subsequently, specific divalent ions modulate populations of intermediates in conformational ensembles along the folding pathway with transition times longer than 10 msec. Rate constants for the five folding transitions act on timescales from submillisecond to tens of seconds. The sensitivity of RNA tertiary structure to divalent cation identity affects all but the fastest events in RNA folding, and allowed us to identify those states that prefer Mg²⁺ The GTPase center RNA appears to have optimized its folding trajectory to specifically utilize this most abundant intracellular divalent ion.

Quantitative sampling of conformational heterogeneity of a DNA hairpin using molecular dynamics simulations and ultrafast fluorescence spectroscopy

Molecular dynamics (MD) simulations and time resolved fluorescence (TRF) spectroscopy were combined to quantitatively describe the conformational landscape of the DNA primary binding sequence (PBS) of the HIV-1 genome, a short hairpin targeted by retroviral nucleocapsid proteins implicated in the viral reverse transcription. Three 2-aminopurine (2AP) labeled PBS constructs were studied. For each variant, the complete distribution of fluorescence lifetimes covering 5 orders of magnitude in timescale was measured and the populations of conformers experimentally observed to undergo static quenching were quantified. A binary quantification permitted the comparison of populations from experimental lifetime amplitudes to populations of aromatically stacked 2AP conformers obtained from simulation. Both populations agreed well, supporting the general assumption that quenching of 2AP fluorescence results from pi-stacking interactions with neighboring nucleobases and demonstrating the success of the proposed methodology for the combined analysis of TRF and MD data. Cluster analysis of the latter further identified predominant conformations that were consistent with the fluorescence decay times and amplitudes, providing a structure-based rationalization for the wide range of fluorescence lifetimes. Finally, the simulations provided evidence of local structural perturbations induced by 2AP. The approach presented is a general tool to investigate fine structural heterogeneity in nucleic acid and nucleoprotein assemblies.

Mode of DNA binding with γ-butyrolactone receptor protein CprB from Streptomyces coelicolor revealed by site-specific fluorescence dynamics

BACKGROUND

The γ-butyrolactone (GBL) binding transcription factors in Streptomyces species are known for their involvement in quorum sensing where they control the expression of various genes initiating secondary metabolic pathways. The structurally characterized member of this family CprB from Streptomyces coelicolor had earlier been demonstrated to bind a multitude of sequences containing a specific binding signature. Though structural breakthrough has been obtained for its complex with a consensus DNA sequence there is, however a dearth of information regarding the overall and site specific dynamics of protein-DNA interaction.

METHODS

To delineate the effect of CprB on the bound DNA, changes in motional dynamics of the fluorescent probe 2-aminopurine were monitored at three conserved base positions (5th, 12th and 23rd) for two DNA sequences: the consensus and the biologically relevant cognate element, on complex formation.

RESULTS

The changes in lifetime and generalized order parameter revealed a similarity in the binding pattern of the protein to both sequences with greater dynamic restriction at the end positions, 5th and 23rd, as compared to the middle 12th position. Also differences within this pattern demonstrated the influence of even small changes in sequence on protein interactions.

CONCLUSIONS

Here the study of motional dynamics was instrumental in establishing a structural footprint for the cognate DNA sequence and explaining the dynamics for the consensus DNA from structural correspondence.

GENERAL SIGNIFICANCE

Motional dynamics can be a powerful tool to efficiently study the mode of DNA binding to proteins that interact differentially with a plethora of DNA sequences, even in the absence of structural breakthrough.

Refinement of Generalized Born Implicit Solvation Parameters for Nucleic Acids and Their Complexes with Proteins

The Generalized Born (GB) implicit solvent model has undergone significant improvements in accuracy for modeling of proteins and small molecules. However, GB still remains a less widely explored option for nucleic acid simulations, in part because fast GB models are often unable to maintain stable nucleic acid structures or they introduce structural bias in proteins, leading to difficulty in application of GB models in simulations of protein-nucleic acid complexes. Recently, GB-neck2 was developed to improve the behavior of protein simulations. In an effort to create a more accurate model for nucleic acids, a similar procedure to the development of GB-neck2 is described here for nucleic acids. The resulting parameter set significantly reduces absolute and relative energy error relative to Poisson-Boltzmann for both nucleic acids and nucleic acid-protein complexes, when compared to its predecessor GB-neck model. This improvement in solvation energy calculation translates to increased structural stability for simulations of DNA and RNA duplexes, quadruplexes, and protein-nucleic acid complexes. The GB-neck2 model also enables successful folding of small DNA and RNA hairpins to near native structures as determined from comparison with experiment. The functional form and all required parameters are provided here and also implemented in the AMBER software.

2-Aminopurine Fluorescence as a Probe of Local RNA Structure and Dynamics and Global Folding

The biology of an RNA is encoded in its structure and dynamics, whether that be binding to a protein, binding to another RNA, enzymatic catalysis, or becoming a substrate. In solution, most RNA molecules are sampling conformations, and their structures are best described as conformational ensembles. For larger RNAs, experiments that can describe the conformations of their domains can be particularly daunting, especially when the RNA is novel and not well characterized. Here, we explain how we have used site-specific 2-aminopurine as a fluorescent probe of the secondary and tertiary structures of a 60 nucleotide RNA, and what new findings we have about its Mg(2+)-dependent conformational changes. We focus on this RNA from prokaryotic ribosome as a proof of concept as well as a research project. Its tertiary structure is known from a cocrystal, and its secondary structure is modeled from phylogenetic conservation, but there are virtually no data describing the motions of its nucleotides in solution, or its folding kinetics. It is a perfect system to illustrate the unique information that comes from a comprehensive fluorescence study of this intricate RNA.

Dynamics of riboswitches: Molecular simulations

Riboswitch RNAs play key roles in bacterial metabolism and represent a promising new class of antibiotic targets for treatment of infectious disease. While many studies of riboswitches have been performed, the exact mechanism of riboswitch operation is still not fully understood at the atomistic level of detail. Molecular dynamics simulations are useful for interpreting existing experimental data and producing predictions for new experiments. Here, a wide range of computational studies on riboswitches is reviewed. By elucidating the key principles of riboswitch operation, computation may aid in the effort to design more specific antibiotics with affinities greater than those of the native ligand. Such a detailed understanding may be required to improve efficacy and reduce side effects. These studies are laying the groundwork for understanding the action mechanism of new compounds that inhibit riboswitch activity. Future directions such as magnesium effects, large-scale conformational changes, expression platforms and co-transcriptional folding are also discussed. This article is part of a Special Issue entitled: Riboswitches.

Structural and thermodynamic signatures that define pseudotriloop RNA hairpins

Pseudotriloop (PTL) structures in RNAs have been recognized as essential elements in RNA folding and recognition of proteins. PTL structures are derived from hexaloops by formation of a cross-loop base pair leaving a triloop and 3' bulged out residue. Despite their common presence and functional importance, insufficient structural and thermodynamic data are available that can be used to predict formation of PTLs from sequence alone. Using NMR spectroscopy and UV-melting data we established factors that contribute to the formation and stability of PTL structures derived from hepatitis B virus and human foamy virus. The NMR data show that, besides the cross-loop base pair, also a 3' pyrimidine bulge and a G-C loop-closing base pair are primary determinants of PTL formation. By changing the G-C closing base pair into C-G, the PTL switches into a hexaloop. Comparison of these rules with regular triloop hairpins and PTLs from other sources is discussed as well as the conservation of a PTL in human foamy virus and other spumaretroviruses.

Sequence-dependent structural dynamics of primate adenosine-to-inosine editing substrates

Humans have the highest level of adenosine-to-inosine (A-to-I) editing amongst primates, yet the reasons for this difference remain unclear. Sequence analysis of the Alu Sg elements (A-to-I RNA substrates) corresponding to the Nup50 gene in human, chimp, and rhesus reveals subtle sequence variations surrounding the edit sites. We have developed three constructs that represent human (HuAp5), chimp (ChAp5), and rhesus (RhAp5) Nup50 Alu Sg A-to-I editing substrates. Here, 2-aminopurine (2-Ap) was substituted for edited adenosine (A5) so as to monitor the fluorescence intensity with respect to temperature. UV and steady-state fluorescence (SSF) T(M) plots indicate that local and global unfolding are coincident, with the human construct displaying a T(M) of approximately 70°C, compared to 60°C for chimp and 54°C for rhesus. However, time-resolved fluorescence (TRF) resolves three different fluorescence lifetimes that we assign to folded, intermediate(s), and unfolded states. The TRF data fit well to a two-intermediate model, whereby both intermediates (M, J) are in equilibrium with each other, and the folded/unfolded states. Our model suggests that, at 37°C, human state J and the folded state will be the most heavily populated in comparison to the other primate constructs. In order for adenosine deaminase acting on RNA (ADAR) to efficiently dock, a stable duplex must be present that corresponds to the human construct, globally. Next, the enzyme must "flip out" the base of interest to facilitate the A-to-I conversion; a nucleotide in an intermediate-like position would enhance this conformational change. Our experiments demonstrate that subtle variations in RNA sequence might contribute to the high A-to-I editing levels found in humans.

GTP binding leads to narrowing of the conformer population while preserving the structure of the RNA aptamer: a site-specific time-resolved fluorescence dynamics study

In this study, we employed a combination of steady-state and time-resolved fluorescence spectroscopy and studied the site-specific dynamics in a GTP aptamer using 2-aminopurine as a fluorescent probe. We compared the dynamics of the GTP-bound aptamer with that of the free aptamer as well as when it is denatured. GTP binding leads to an overall compaction of structure in the aptamer. The general pattern of fluorescence lifetimes and correlation times scanned across several locations in the aptamer does not seem to change following GTP binding. However, a remarkable narrowing of the lifetime distribution of the aptamer ensues following its compaction by GTP binding. Interestingly, such a "conformational narrowing" is evident from the lifetime readouts of the nucleotide belonging to the stem as well as the "bulge" part of the aptamer, independent of whether it is directly interacting with GTP. Taken together, these results underscore the importance of an overall intrinsic structure associated with the free aptamer that is further modulated following GTP binding. This work provides strong support for the "conformational selection" hypothesis of ligand binding.

Spectroscopic probes of RNA structure and dynamics

A single technique is insufficient to characterize the properties of an RNA molecule so this Guide provides advice and suggestions for use of several spectroscopic methods applied to RNA. It begins with a discussion of design features to synthesize a suitable molecule for study, assuming that the reader is familiar with in vitro enzymatic and chemical synthesis of RNA. With the RNA in hand, the application of UV melting studies to characterize its folding/unfolding transitions is presented, followed by instructions for fluorescence experiments to augment the UV data. Since RNAs are flexible molecules, it is important to characterize their mobility, and so time-resolved fluorescence data are analyzed. A brief exposition of the power of NMR spectroscopy to identify ion-binding sites is provided. RNA examples are described for each method to give a perspective on what can and cannot be learned.

Applying 6-methylisoxanthopterin-enhanced fluorescence to examine protein-DNA interactions in the picomolar range

Incorporation of fluorescent nucleoside analogues into duplex DNA usually leads to a reduction in quantum yield, which significantly limits their potential use and application. We have identified two pentamer DNA sequences containing 6-methylisoxanthopterin (6-MI) (ATFAA and AAFTA, where F is 6-MI) that exhibit significant enhancement of fluorescence upon formation of duplex DNA with quantum yields close to that of monomeric 6-MI. The enhanced fluorescence dramatically increases the utility and sensitivity of the probe and is used to study protein-DNA interactions of nanomolar specificity in this work. The increased sensitivity of 6-MI allows anisotropy binding measurements to be performed at DNA concentrations of 1 nM and fluorescence intensity measurements at 50 pM DNA. The ATFAA sequence was incorporated into DNA constructs to measure the binding affinity of four different protein-DNA interactions that exhibit sequence-specific and non-sequence-specific recognition. In all cases, the K(d) values obtained were consistent with previously reported values measured by other methods. Time-resolved and steady-state fluorescence measurements demonstrate that 6-MI fluorescence is very sensitive to local distortion and reports on different degrees of protein-induced perturbations with single-base resolution, where the largest changes occur at the site of protein binding.

Fe2+ binds iron responsive element-RNA, selectively changing protein-binding affinities and regulating mRNA repression and activation

Iron increases synthesis rates of proteins encoded in iron-responsive element (IRE)-mRNAs; metabolic iron ("free," "labile") is Fe(2+). The noncoding IRE-RNA structure, approximately 30 nt, folds into a stem loop to control synthesis of proteins in iron trafficking, cell cycling, and nervous system function. IRE-RNA riboregulators bind specifically to iron-regulatory proteins (IRP) proteins, inhibiting ribosome binding. Deletion of the IRE-RNA from an mRNA decreases both IRP binding and IRP-independent protein synthesis, indicating effects of other "factors." Current models of IRE-mRNA regulation, emphasizing iron-dependent degradation/modification of IRP, lack answers about how iron increases IRE-RNA/IRP protein dissociation or how IRE-RNA, after IRP dissociation, influences protein synthesis rates. However, we observed Fe(2+) (anaerobic) or Mn(2+) selectively increase the IRE-RNA/IRP K(D). Here we show: (i) Fe(2+) binds to the IRE-RNA, altering its conformation (by 2-aminopurine fluorescence and ethidium bromide displacement); (ii) metal ions increase translation of IRE-mRNA in vitro; (iii) eukaryotic initiation factor (eIF)4F binds specifically with high affinity to IRE-RNA; (iv) Fe(2+) increased eIF4F/IRE-RNA binding, which outcompetes IRP binding; (v) exogenous eIF4F rescued metal-dependent IRE-RNA translation in eIF4F-depeleted extracts. The regulation by metabolic iron binding to IRE-RNA to decrease inhibitor protein (IRP) binding and increase activator protein (eIF4F) binding identifies IRE-RNA as a riboregulator.

Comparative structural effects of HIV-1 Gag and nucleocapsid proteins in binding to and unwinding of the viral RNA packaging signal

The major RNA binding region of the HIV-1 Gag polyprotein is the nucleocapsid (NC) domain, which is responsible for the specific capture of the genomic RNA genome during viral assembly. The Gag polyprotein has other RNA chaperone functions, which are mirrored by the isolated NC protein after physiological cleavage from Gag. Gag, however, is suggested to have superior nucleic acid chaperone activity. Here we investigate the interaction of Gag and NC with the core RNA structure of the HIV-1 packaging signal (Ψ), using 2-aminopurine substitution to create a series of modified RNAs based on the Ψ helix loop structure. The effects of 2-aminopurine substitution on the physical and structural properties of the viral Ψ were characterized. The fluorescence properties of the 2-aminopurine substitutions showed features consistent with the native GNAR tetraloop. Dissociation constants (K(d)) of the two viral proteins, measured by fluorescence polarization (FP), were similar, and both NC and Gag affected the 2-aminopurine fluorescence of bases close to the loop binding region in a similar fashion. However, the influence of Gag on the fluorescence of the 2-aminopurine nucleotides at the base of the helix implied a much more potent helix destabilizing action on the RNA stem loop (SL) versus that seen with NC. This was further supported when the viral Ψ SL was tagged with a 5' fluorophore and 3' quencher. In the absence of any viral protein, minimal fluorescence was detected; addition of NC yielded a slight increase in fluorescence, while addition of the Gag protein yielded a large change in fluorescence, further suggesting that, compared to NC, the Gag protein has a greater propensity to affect RNA structure and that Ψ helix unwinding may be an intrinsic step in RNA encapsidation.

Characterization of a K+-induced conformational switch in a human telomeric DNA oligonucleotide using 2-aminopurine fluorescence

Human telomeric DNA consists of tandem repeats of the DNA sequence d(GGGTTA). Oligodeoxynucleotide telomere models such as d[A(GGGTTA)(3)GGG] (Tel22) fold in a cation-dependent manner into quadruplex structures consisting of stacked G-quartets linked by d(TTA) loops. NMR has shown that in Na(+) solutions Tel22 forms a "basket" topology of four antiparallel strands; in contrast, Tel22 in K(+) solutions consists of a mixture of unknown topologies. Our previous studies on the mechanism of folding of Tel22 and similar telomere analogues utilized changes in UV absorption between 270 and 325 nm that report primarily on G-quartet formation and stacking showed that quadruplex formation occurs within milliseconds upon mixing with an appropriate cation. In this study, we assess the dynamics and equilibria of folding of specific loops by using Tel22 derivatives in which the dA residues were serially substituted with the fluorescent reporter base, 2-aminopurine (2-AP). Tel22 folding induced by Na(+) or K(+) assessed by changes in 2-AP fluorescence consists of at least three kinetic steps with time constants spanning a range from milliseconds to several hundred seconds. Na(+)-dependent equilibrium titrations of Tel22 folding could be approximated as a cooperative two-state process. In contrast, K(+)-dependent folding curves were biphasic, revealing that different conformational ensembles are present in 1 and 30 mM K(+). This conclusion was confirmed by (1)H NMR. Molecular dynamics simulations revealed a K(+) binding pocket in Tel22 located near dA1 that is specific for the so-called hybrid-1 conformation in which strand 1 is in a parallel arrangement. The possible presence of this topologically specific binding site suggests that K(+) may play an allosteric role in regulating telomere conformation and function by modulating quadruplex tertiary structure.

Fast folding of an RNA tetraloop on a rugged energy landscape detected by a stacking-sensitive probe

We investigate the microsecond-timescale kinetics of the RNA hairpin ga*cUUCGguc. The fluorescent nucleotide 2-aminopurine (a*) reports mainly on base stacking. Ten kinetic traces and the temperature denaturation curve are globally fitted to four-state models of the free-energy surface. In the best-fitting sequential model, the hairpin unfolds over successively larger barriers in at least three stages: stem fraying and increased base-stacking fluctuations; concerted loss of hydrogen bonding and partial unstacking; and additional unstacking of single strands at the highest temperatures. Parallel and trap models also provide adequate fits: such pathways probably also play a role in the complete free-energy surface of the hairpin. To interpret the model states structurally, 200 ns of molecular dynamics, including six temperature-jump simulations, were run. Although the sampling is by no means comprehensive, five different states were identified using hydrogen bonding and base stacking as reaction coordinates. The four to five states required to explain the experiments or simulations set a lower limit on the complexity of this small RNA hairpin's energy landscape.

A DNA nanoswitch incorporating the fluorescent base analogue 2-aminopurine detects single nucleotide mismatches in unlabelled targets

DNA nanoswitches can be designed to detect unlabelled nucleic acid targets and have been shown to discriminate between targets which differ in the identity of only one base. This paper demonstrates that the fluorescent base analogue 2-aminopurine (AP) can be used to discriminate between nanoswitches with and without targets and to discriminate between matched and mismatched targets. In particular, we have used both steady-state and time-resolved fluorescence spectroscopy to determine differences in AP environment at the branchpoint of nanoswitches assembled using complementary targets and targets which incorporate single base mismatches.

Dissecting structural transitions in the HIV-1 dimerization initiation site RNA using 2-aminopurine fluorescence

A highly conserved 35 nucleotide RNA stem-loop, the dimerization initiation site (DIS), in the 5' untranslated region (UTR) of the human immunodeficiency virus type I (HIV-1) genome has been identified as the sequence primarily responsible for initiation of viral genome dimerization. The DIS initiates viral genome dimerization through a loop-loop 'kissing' interaction and is converted from an intermediate 'kissing' to a more thermodynamically stable extended duplex dimer in a conformational rearrangement that is chaperoned by the HIV-1 nucleocapsid protein (NCp7). Here we describe fluorescence methods designed to probe local RNA dynamics and structural transitions associated with the DIS dimer formation and its NCp7 chaperoned structural conversion. These methods take advantage of the exquisite sensitivity of the quantum yield of the fluorescent nucleotide base analog, 2-aminopurine (2-AP), to its immediate structural and dynamic environment. The 2-AP fluorescence methods described allow a detailed kinetic and thermodynamic examination of this type of RNA-RNA interaction, as well as an analysis of the molecular mechanism of NCp7 chaperone activity.

2-aminopurine as a probe of RNA conformational transitions

2-aminopurine (2AP) is a fluorescent nucleobase that provides the means to probe structure and dynamics of RNA molecules. Because 2AP can base pair with Uridine, it can replace normal A:U pairs without substantial deformation of duplexes. It is best used as a probe of ostensibly single-stranded regions: its fluorescence intensity reports on base stacking and its fluorescence decay lifetimes report on its conformational dynamics. Three examples of its use are described here, illustrating how 2AP fluorescence has been used to probe RNA folding and hairpin loop dynamics.

RNA folding: conformational statistics, folding kinetics, and ion electrostatics

RNA folding is a remarkably complex problem that involves ion-mediated electrostatic interaction, conformational entropy, base pairing and stacking, and noncanonical interactions. During the past decade, results from a variety of experimental and theoretical studies pointed to (a) the potential ion correlation effect in Mg2+-RNA interactions, (b) the rugged energy landscapes and multistate RNA folding kinetics even for small RNA systems such as hairpins and pseudoknots, (c) the intraloop interactions and sequence-dependent loop free energy, and (d) the strong nonadditivity of chain entropy in RNA pseudoknot and other tertiary folds. Several related issues, which have not been thoroughly resolved, require combined approaches with thermodynamic and kinetic experiments, statistical mechanical modeling, and all-atom computer simulations.

Time-resolved fluorescent properties of 8-vinyl-deoxyadenosine and 2-amino-deoxyribosylpurine exhibit different sensitivity to their opposite base in duplexes

8-Vinyl-deoxyadenosine (8VA) has been recently introduced as a fluorescent analogue of adenosine that is less perturbing and less quenched than the well-established 2-amino-deoxyribosylpurine (2AP) probe when inserted in oligonucleotides. To further validate 8VA as a fluorescent substitute of A, we compared the ability of 8VA and 2AP in sequences of the type d(CGT TTT XNX TTT TGC) (with N=8VA or 2AP and X=T and C) to discriminate the nature of the opposite base (Y) in duplexes. For both probes, systematic variations in the amplitudes of the short- and long-lived lifetimes of the fluorescence intensity decays as well as in the amplitude of the fast rotational correlation time of the fluorescence anisotropy decays were observed as a function of the nature of Y. From these parameters, we inferred a stability order 8VA-T > 8VA-G > 8VA-A > 8VA-C, similar to the stability order with the native A base, but different from the stability order with 2AP. Using a combination of molecular mechanics and ab initio calculations, we found that the time-resolved parameters of 8VA, but not the 2AP ones, correlate well with the geometry and the strength of the A-Y base-pairing interaction. This may be rationalized by the smaller structural and electronic perturbations induced by the vinyl group in position 8 as compared to the amino group at position 2. As a consequence, substitution of A by 8VA in a base pair was found to only minimally modify the structure and interaction energy of the base pair. Thus, 8VA can be used as a fluorescent substitute of the natural A, to straightforwardly discriminate the nature of the opposite base. This may find interesting applications notably in the elucidation of the mechanisms and dynamics of the DNA mismatch repair system.

Site-specific variations in RNA folding thermodynamics visualized by 2-aminopurine fluorescence

The fluorescent base analogue 2-aminopurine (2-AP) is commonly used to study specific conformational and protein binding events involving nucleic acids. Here, combinations of steady-state and time-resolved fluorescence spectroscopy of 2-AP were employed to monitor conformational transitions within a model hairpin RNA from diverse structural perspectives. RNA substrates adopting stable, unambiguous secondary structures were labeled with 2-AP at an unpaired base, within the loop, or inside the base-paired stem. Steady-state fluorescence was monitored as the RNA hairpins made the transitions between folded and unfolded conformations using thermal denaturation, urea titration, and cation-mediated folding. Unstructured control RNA substrates permitted the effects of higher-order RNA structures on 2-AP fluorescence to be distinguished from stimulus-dependent changes in intrinsic 2-AP photophysics and/or interactions with adjacent residues. Thermodynamic parameters describing local conformational changes were thus resolved from multiple perspectives within the model RNA hairpin. These data provided energetic bases for construction of folding mechanisms, which varied among different folding-unfolding stimuli. Time-resolved fluorescence studies further revealed that 2-AP exhibits characteristic signatures of component fluorescence lifetimes and respective fractional contributions in different RNA structural contexts. Together, these studies demonstrate localized conformational events contributing to RNA folding and unfolding that could not be observed by approaches monitoring only global structural transitions.

Structural heterogeneity in DNA: temperature dependence of 2-aminopurine fluorescence in dinucleotides

The fluorescent base analogue 2-aminopurine is a sensitive probe for local dynamics of DNA. Its fluorescence is quenched by interaction with the neighboring bases, but the underlying mechanisms are still under investigation. We studied 2-aminopurine fluorescence in dinucleotides with each of the natural bases. Consistently, two of the four fluorescence-decay components depend strongly on temperature. Our results indicate that these components are due to the excited-state dynamics of a single conformational state. We propose a variation of the gating model in which transient unstacking occurs in the excited state.

Solution structure of the apical stem-loop of the human hepatitis B virus encapsidation signal

Hepatitis B virus (HBV) replication is initiated by HBV RT binding to the highly conserved encapsidation signal, epsilon, at the 5′ end of the RNA pregenome. Epsilon contains an apical stem–loop, whose residues are either totally conserved or show rare non-disruptive mutations. Here we present the structure of the apical stem–loop based on NOE, RDC and ¹H chemical shift NMR data. The ¹H chemical shifts proved to be crucial to define the loop conformation. The loop sequence 5′-CUGUGC-3′ folds into a UGU triloop with a CG closing base pair and a bulged out C and hence forms a pseudo-triloop, a proposed protein recognition motif. In the UGU loop conformations most consistent with experimental data, the guanine nucleobase is located on the minor groove face and the two uracil bases on the major groove face. The underlying helix is disrupted by a conserved non-paired U bulge. This U bulge adopts multiple conformations, with the nucleobase being located either in the major groove or partially intercalated in the helix from the minor groove side, and bends the helical stem. The pseudo-triloop motif, together with the U bulge, may represent important anchor points for the initial recognition of epsilon by the viral RT.

Molecular control of vertebrate iron homeostasis by iron regulatory proteins

Both deficiencies and excesses of iron represent major public health problems throughout the world. Understanding the cellular and organismal processes controlling iron homeostasis is critical for identifying iron-related diseases and in advancing the clinical treatments for such disorders of iron metabolism. Iron regulatory proteins (IRPs) 1 and 2 are key regulators of vertebrate iron metabolism. These RNA binding proteins post-transcriptionally control the stability or translation of mRNAs encoding proteins involved in iron homeostasis thereby controlling the uptake, utilization, storage or export of iron. Recent evidence provides insight into how IRPs selectively control the translation or stability of target mRNAs, how IRP RNA binding activity is controlled by iron-dependent and iron-independent effectors, and the pathological consequences of dysregulation of the IRP system.

Structural Fluctuations and Excitation Transfer between Adenine and 2-Aminopurine in Single-Stranded Deoxytrinucleotides

Steady-state fluorescence measurements on the deoxytrinucleotides (5')dTp2APpA(3') and (5')dAp2APpA(3') show a temperature-dependence and a viscosity-dependence for energy transfer that qualitatively differ from those seen in our previous study of charge transfer (CT) in these systems. Time-resolved anisotropy studies and molecular dynamics simulations are presented that provide a detailed characterization of the structural dynamics of these systems and how these fluctuations modulate the electronic interaction between 2AP and its neighbors. To gain quantitative insight into the interplay of conformational fluctuations and stacking-induced energy transfer, we present results from a new hybrid quantum-classical simulation method for computing the A --> 2AP energy transfer rate that makes use of the full three-dimensional nature of the donor and acceptor transition densities. Analysis of the results shows that the standard transition dipole-transition dipole approximation for the Coulombic coupling substantially overestimates the transfer rate and that the nearest neighbor energy transfer from adenine to 2AP occurs on a much faster time scale than that for CT. This suggests that, unlike the CT dynamics where conformational "gating" plays a critical role, the large amplitude fluctuations that modulate the process are largely "frozen" out on the energy transfer time scale.

Exploring the complex folding kinetics of RNA hairpins: I. General folding kinetics analysis

Depending on the nucleotide sequence, the temperature, and other conditions, RNA hairpin-folding kinetics can be very complex. The complexity with a wide range of cooperative and noncooperative kinetic behaviors arises from the interplay between the formation of the loops, the disruption of the misfolded states, and the formation of the rate-limiting base stacks. With a rate constant model and a kinetic-cluster theory, we explore the broad landscape for RNA hairpin-folding kinetics. The model is validated through direct tests against several experimental measurements. The general kinetic folding mechanisms and the predicted great variety of folding kinetics are directly applicable and quantitatively testable in experiments. The results from this study suggest that 1), previous experimental findings based on the individual hairpins revealed only a small fraction of much broader and more complex RNA hairpin-folding landscapes; 2), even for structures as simple as hairpins, universal folding timescales and pathways do not exist; and 3), to treat the loop size as the sole factor to determine the hairpin-folding rate is an oversimplification.

Exploring the complex folding kinetics of RNA hairpins: II. Effect of sequence, length, and misfolded states

The complexity of RNA hairpin folding arises from the interplay between the loop formation, the disruption of the slow-breaking misfolded states, and the formation of the slow-forming native base stacks. We investigate the general physical mechanism for the dependence of the RNA hairpin folding kinetics on the sequence and the length of the hairpin loop and the helix stem. For example, 1), the folding would slow down when a stable GC basepair moves to the middle of the stem; 2), hairpin with GC basepair near the loop would fold/unfold faster than the one with GC near the tail of the stem; 3), within a certain range of the stem length, a longer stem can cause faster folding; and 4), certain misfolded states can assist folding through the formation of scaffold structures to lower the entropic barrier for the folding. All our findings are directly applicable and quantitatively testable in experiments. In addition, our results can be useful for molecular design to achieve desirable fast/slow-folding hairpins, hairpins with/without specific misfolded intermediates, and hairpins that fold along designed pathways.

Site-directed spin labeling studies reveal solution conformational changes in a GAAA tetraloop receptor upon Mg(2+)-dependent docking of a GAAA tetraloop

The Mg(2+)-dependent GAAA tetraloop interaction with its 11 nucleotide receptor is one of the most frequently occurring long-range tertiary interactions in RNAs. To explore conformational changes in the receptor during tetraloop docking, nitroxide spin labels were attached at each of four uridine bases, one at a time, within an RNA molecule containing the receptor sequence. In the presence of Mg2+ and the tetraloop, the electron paramagnetic resonance (EPR) spectrum of one of the labeled bases reflected a large increase in mobility, indicating unstacking of the base upon tetraloop docking. This provides direct evidence that base unstacking is an intrinsic feature of the solution tetraloop-receptor complex formed in the presence of Mg2+. Additional evidence suggests that in solution the bound receptor conformation is similar to that observed in the crystal structure of a group I intron ribozyme domain. In Mg2+ alone, a receptor conformation with an unstacked base was not detectable, suggesting that this conformation is of higher standard state free energy than that of the free receptor. This leads to the conclusion that the extensive RNA-RNA interactions observed in the crystal structure of the tetraloop-receptor complex provide larger interaction energy than the measured apparent affinity between the tetraloop and the free receptor. This is compatible with a high specificity of the tetraloop-receptor interaction.

Iron responsive element RNA flexibility described by NMR and isotropic reorientational eigenmode dynamics

The first example of the application of reorientational eigenmode dynamics (RED) to RNA is shown here for the small and floppy Iron Responsive Element (IRE) RNA hairpin. Order parameters calculated for bases and riboses from a 12 ns molecular dynamics trajectory are compared to experimentally determined order parameters from 13C-1H NMR relaxation experiments, and shown to be in qualitative agreement. Given the small size of the IRE hairpin and its very flexible loop, isotropic RED (iRED) was also used to analyze the trajectory in order to describe its dynamic motions. iRED analysis shows that the global and internal dynamics of the IRE are not rigorously separable, which will result in inaccurate experimental order parameters. In addition, the iRED analysis described the many correlated motions that comprise the dynamics of the IRE RNA. The combined use of NMR relaxation, RED, and iRED provide a uniquely detailed description of IRE RNA dynamics.