Detecting protein-induced folding of the U4 snRNA kink-turn by single-molecule multiparameter FRET measurements

Connecting genotype and phenotype in minor spliceosome diseases

Minor spliceosome is responsible for recognizing and excising a specific subset of divergent introns during the pre-mRNA splicing process. Mutations in the unique snRNA and protein components of the minor spliceosome are increasingly being associated with a variety of germline and somatic human disorders, collectively termed as minor spliceosomopathies. Understanding the mechanistic basis of these diseases has been challenging due to limited functional information on many minor spliceosome components. However, recently published cryo-electron microscopy (cryo-EM) structures of various minor spliceosome assembly intermediates have marked a significant advancement in elucidating the roles of these components during splicing. These structural breakthroughs have not only enhanced our comprehension of the minor spliceosome's functionality but also shed light on how disease-associated mutations disrupt its functions. Consequently, research focus is now shifting toward investigating how these splicing defects translate into broader pathological processes within gene expression pathways. Here we outline the current structural and functional knowledge of the minor spliceosome, explore the mechanistic consequences of its mutations, and discuss emerging challenges in connecting molecular dysfunctions to clinical phenotypes.

RNA kink-turns are highly anisotropic with respect to lateral displacement of the flanking stems

Kink-turns are highly bent internal loop motifs commonly found in the ribosome and other RNA complexes. They frequently act as binding sites for proteins and mediate tertiary interactions in larger RNA structures. Kink-turns have been a topic of intense research, but their elastic properties in the folded state are still poorly understood. Here we use extensive all-atom molecular dynamics simulations to parameterize a model of kink-turn in which the two flanking helical stems are represented by effective rigid bodies. Time series of the full set of six interhelical coordinates enable us to extract minimum energy shapes and harmonic stiffness constants for kink-turns from different RNA functional classes. The analysis suggests that kink-turns exhibit isotropic bending stiffness but are highly anisotropic with respect to lateral displacement of the stems. The most flexible lateral displacement mode is perpendicular to the plane of the static bend. These results may help understand the structural adaptation and mechanical signal transmission by kink-turns in complex natural and artificial RNA structures.

The kink-turn in the structural biology of RNA

The kink-turn (k-turn) is a widespread structural motif found in functional RNA species. It typically comprises a three-nucleotide bulge followed by tandem trans sugar edge-Hoogsteen G:A base pairs. It introduces a sharp kink into the axis of duplex RNA, juxtaposing the minor grooves. Cross-strand H-bonds form at the interface, accepted by the conserved adenine nucleobases of the G:A basepairs. Alternative acceptors for one of these divides the k-turns into two conformational classes N3 and N1. The base pair that follows the G:A pairs (3b:3n) determines which conformation is adopted by a given k-turn. k-turns often mediate tertiary contacts in folded RNA species and frequently bind proteins. Common k-turn binding proteins include members of the L7Ae family, such as the human 15·5k protein. A recognition helix within these proteins binds in the widened major groove on the outside of the k-turn, that makes specific H-bonds with the conserved guanine nucleobases of the G:A pairs. L7Ae binds with extremely high affinity, and single-molecule data are consistent with folding by conformational selection. The standard, simple k-turn can be elaborated in a variety of ways, that include the complex k-turns and the k-junctions. In free solution in the absence of added metal ions or protein k-turns do not adopt the tightly-kinked conformation. They undergo folding by the binding of proteins, by the formation of tertiary contacts, and some (but not all) will fold on the addition of metal ions. Whether or not folding occurs in the presence of metal ions depends on local sequence, including the 3b:3n position, and the -1b:-1n position (5' to the bulge). In most cases -1b:-1n = C:G, so that the 3b:3n position is critical since it determines both folding properties and conformation. In general, the selection of these sequence matches a given k-turn to its biological requirements. The k-turn structure is now very well understood, to the point at which they can be used as a building block for the formation of RNA nano-objects, including triangles and squares.

Methodologies for studying the spliceosome's RNA dynamics with single-molecule FRET

The spliceosome is an extraordinarily dynamic molecular machine in which significant changes in composition as well as protein and RNA conformation are required for carrying out pre-mRNA splicing. Single-molecule fluorescence resonance energy transfer (smFRET) can be used to elucidate these dynamics both in well-characterized model systems and in entire spliceosomes. These types of single-molecule data provide novel information about spliceosome components and can be used to identify sub-populations of molecules with unique behaviors. When smFRET is combined with single-molecule fluorescence colocalization, conformational dynamics can be further linked to the presence or absence of a given spliceosome component. Here, we provide a description of experimental considerations, approaches, and workflows for smFRET with an emphasis on applications for the splicing machinery.

RNAComposer and RNA 3D structure prediction for nanotechnology

RNAs adopt specific, stable tertiary architectures to perform their activities. Knowledge of RNA tertiary structure is fundamental to understand RNA functions beginning with transcription and ending with turnover. Contrary to advanced RNA secondary structure prediction algorithms, which allow good accuracy when experimental data are integrated into the prediction, tertiary structure prediction of large RNAs still remains a significant challenge. However, the field of RNA tertiary structure prediction is rapidly developing and new computational methods based on different strategies are emerging. RNAComposer is a user-friendly and freely available server for 3D structure prediction of RNA up to 500 nucleotide residues. RNAComposer employs fully automated fragment assembly based on RNA secondary structure specified by the user. Importantly, this method allows incorporation of distance restraints derived from the experimental data to strengthen the 3D predictions. The potential and limitations of RNAComposer are discussed and an application to RNA design for nanotechnology is presented.

The solution structural ensembles of RNA kink-turn motifs and their protein complexes

With the growing number of crystal structures of RNA and RNA-protein complexes, a critical next step is understanding the dynamic solution behavior of these entities in terms of conformational ensembles and energy landscapes. To this end, we have used X-ray scattering interferometry (XSI) to probe the ubiquitous RNA kink-turn motif and its complexes with the canonical kink-turn binding protein L7Ae. XSI revealed that the folded kink-turn is best described as a restricted conformational ensemble. The ions present in solution alter the nature of this ensemble, and protein binding can perturb the kink-turn ensemble without collapsing it to a unique state. This study demonstrates how XSI can reveal structural and ensemble properties of RNAs and RNA-protein complexes and uncovers the behavior of an important RNA-protein motif. This type of information will be necessary to understand, predict and engineer the behavior and function of RNAs and their protein complexes.

The Kink Turn, a Key Architectural Element in RNA Structure

Kink turns (k-turns) are widespread structural elements that introduce an axial bend into duplex RNA with an included angle of 50°. These mediate key tertiary interactions and bind specific proteins including members of the L7Ae family. The standard k-turn comprises a three-nucleotide bulge followed by G·A and A·G pairs. The RNA kinks by an association of the two minor grooves, stabilized by the formation of a number of key cross-strand hydrogen bonds mostly involving the adenine bases of the G·A and A·G pairs. The k-turns may be divided into two conformational classes, depending on the receptor for one of these hydrogen bonds. k-turns become folded by one of three different processes. Some, but not all, k-turns become folded in the presence of metal ions. Whether or not a given k-turn is folded under these conditions is determined by its sequence. We present a set of rules for the prediction of folding properties and the structure adopted on local sequence.

Automated RNA 3D Structure Prediction with RNAComposer

RNAs adopt specific structures to perform their activities and these are critical to virtually all RNA-mediated processes. Because of difficulties in experimentally assessing structures of large RNAs using NMR, X-ray crystallography, or cryo-microscopy, there is currently great demand for new high-resolution 3D structure prediction methods. Recently we reported on RNAComposer, a knowledge-based method for the fully automated RNA 3D structure prediction from a user-defined secondary structure. RNAComposer method is especially suited for structural biology users. Since our initial report in 2012, both servers, freely available at http://rnacomposer.ibch.poznan.pl and http://rnacomposer.cs.put.poznan.pl have been often visited. Therefore this chapter provides guidance for using RNAComposer and discusses points that should be considered when predicting 3D RNA structure. An application example presents current scope and limitations of RNAComposer.

Resolution of Two Sub-Populations of Conformers and Their Individual Dynamics by Time Resolved Ensemble Level FRET Measurements

Most active biopolymers are dynamic structures; thus, ensembles of such molecules should be characterized by distributions of intra- or intermolecular distances and their fast fluctuations. A method of choice to determine intramolecular distances is based on Förster resonance energy transfer (FRET) measurements. Major advances in such measurements were achieved by single molecule FRET measurements. Here, we show that by global analysis of the decay of the emission of both the donor and the acceptor it is also possible to resolve two sub-populations in a mixture of two ensembles of biopolymers by time resolved FRET (trFRET) measurements at the ensemble level. We show that two individual intramolecular distance distributions can be determined and characterized in terms of their individual means, full width at half maximum (FWHM), and two corresponding diffusion coefficients which reflect the rates of fast ns fluctuations within each sub-population. An important advantage of the ensemble level trFRET measurements is the ability to use low molecular weight small-sized probes and to determine nanosecond fluctuations of the distance between the probes. The limits of the possible resolution were first tested by simulation and then by preparation of mixtures of two model peptides. The first labeled polypeptide was a relatively rigid Pro₇ and the second polypeptide was a flexible molecule consisting of (Gly-Ser)₇ repeats. The end to end distance distributions and the diffusion coefficients of each peptide were determined. Global analysis of trFRET measurements of a series of mixtures of polypeptides recovered two end-to-end distance distributions and associated intramolecular diffusion coefficients, which were very close to those determined from each of the pure samples. This study is a proof of concept study demonstrating the power of ensemble level trFRET based methods in resolution of subpopulations in ensembles of flexible macromolecules.

Assembly and dynamics of the U4/U6 di-snRNP by single-molecule FRET

In large ribonucleoprotein machines, such as ribosomes and spliceosomes, RNA functions as an assembly scaffold as well as a critical catalytic component. Protein binding to the RNA scaffold can induce structural changes, which in turn modulate subsequent binding of other components. The spliceosomal U4/U6 di-snRNP contains extensively base paired U4 and U6 snRNAs, Snu13, Prp31, Prp3 and Prp4, seven Sm and seven LSm proteins. We have studied successive binding of all protein components to the snRNA duplex during di-snRNP assembly by electrophoretic mobility shift assay and accompanying conformational changes in the U4/U6 RNA 3-way junction by single-molecule FRET. Stems I and II of the duplex were found to co-axially stack in free RNA and function as a rigid scaffold during the entire assembly, but the U4 snRNA 5' stem-loop adopts alternative orientations each stabilized by Prp31 and Prp3/4 binding accounting for altered Prp3/4 binding affinities in presence of Prp31.

AMBER-DYES: Characterization of Charge Fluctuations and Force Field Parameterization of Fluorescent Dyes for Molecular Dynamics Simulations

Recent advances in single molecule fluorescence experiments and theory allow a direct comparison and improved interpretation of experiment and simulation. To this end, force fields for a larger number of dyes are required which are compatible with and can be integrated into existing biomolecular force fields. Here, we developed, characterized, and implemented AMBER-DYES, a modular fluorescent label force field, for a set of 22 fluorescent dyes and their linkers from the Alexa, Atto, and Cy families, which are in common use for single molecule spectroscopy experiments. The force field is compatible with the AMBER protein force fields and the GROMACS molecular dynamics simulation program. The high electronic polarizability of the delocalized π-electron orbitals, as found in many fluorescent dyes, poses a particular challenge to point charge based force fields such as AMBER. To quantify the charge fluctuations due to the electronic polarizability, we simulated the 22 dyes in explicit solvent and sampled the charge fluctuations using QM/MM simulations at the B3LYP/6-31G*//TIP3P level of theory. The analysis of the simulations enabled us to derive ensemble fitted RESP charges from the solvated charge distributions of multiple trajectories. We observed broad, single peaked charge distributions for the conjugated ring atoms with well-defined mean values. The charge fitting procedure was validated against published charges of the dyelike amino acid tryptophan, which showed good agreement with existing tryptophan parameters from the AMBER, CHARMM, and OPLS force field families. A principal component analysis of the charge fluctuations revealed that a small number of collective coordinates suffices to describe most of the in-plane dye polarizability. The AMBER-DYES force field allows the rapid preparation of all atom molecular dynamics simulations of fluorescent systems for state of the art multi microsecond trajectories.

Biochemical defects in minor spliceosome function in the developmental disorder MOPD I

Biallelic mutations of the human RNU4ATAC gene, which codes for the minor spliceosomal U4atac snRNA, cause the developmental disorder, MOPD I/TALS. To date, nine separate mutations in RNU4ATAC have been identified in MOPD I patients. Evidence suggests that all of these mutations lead to abrogation of U4atac snRNA function and impaired minor intron splicing. However, the molecular basis of these effects is unknown. Here, we use a variety of in vitro and in vivo assays to address this question. We find that only one mutation, 124G>A, leads to significantly reduced expression of U4atac snRNA, whereas four mutations, 30G>A, 50G>A, 50G>C and 51G>A, show impaired binding of essential protein components of the U4atac/U6atac di-snRNP in vitro and in vivo. Analysis of MOPD I patient fibroblasts and iPS cells homozygous for the most common mutation, 51G>A, shows reduced levels of the U4atac/U6atac.U5 tri-snRNP complex as determined by glycerol gradient sedimentation and immunoprecipitation. In this report, we establish a mechanistic basis for MOPD I disease and show that the inefficient splicing of genes containing U12-dependent introns in patient cells is due to defects in minor tri-snRNP formation, and the MOPD I-associated RNU4ATAC mutations can affect multiple facets of minor snRNA function.

The K-turn motif in riboswitches and other RNA species

The kink turn is a widespread structure motif that introduces a tight bend into the axis of duplex RNA. This generally functions to mediate tertiary interactions, and to serve as a specific protein binding site. K-turns or closely related structures are found in at least seven different riboswitch structures, where they function as key architectural elements that help generate the ligand binding pocket. This article is part of a Special Issue entitled: Riboswitches.

Kink turn sRNA folding upon L7Ae binding using molecular dynamics simulations

The kink-turn sRNA motif in archaea, whose combination with protein L7Ae initializes the assembly of small ribonucleoprotein particles (sRNPs), plays a key role in ribosome maturation and the translation process. Although many studies have been reported on this motif, the mechanism of sRNA folding coupled with protein binding is still poorly understood. Here, room and high temperature molecular dynamics (MD) simulations were performed on the complex of 25-nt kink-turn sRNA and L7Ae. The average RMSD values between the bound and corresponding apo structures and Kolmogorov-Smirnov P test analysis indicate that sRNA may follow an induced fit mechanism upon binding with L7Ae, both locally and globally. These conclusions are further supported by high-temperature unfolding kinetic analysis. Principal component analysis (PCA) found both closing and opening motions of the kink-turn sRNA. This might play a key role in the sRNP assembly and methylation catalysis. These combined computational methods can be used to study the specific recognition of other sRNAs and proteins.

The plasticity of a structural motif in RNA: structural polymorphism of a kink turn as a function of its environment

The k-turn is a widespread structural motif that introduces a tight kink into the helical axis of double-stranded RNA. The adenine bases of consecutive G•A pairs are directed toward the minor groove of the opposing helix, hydrogen bonding in a typical A-minor interaction. We show here that the available structures of k-turns divide into two classes, depending on whether N3 or N1 of the adenine at the 2b position accepts a hydrogen bond from the O2' at the -1n position. There is a coordinated structural change involving a number of hydrogen bonds between the two classes. We show here that Kt-7 can adopt either the N3 or N1 structures depending on environment. While it has the N1 structure in the ribosome, on engineering it into the SAM-I riboswitch, it changes to the N3 structure, resulting in a significant alteration in the trajectory of the helical arms.

Single-molecule observation of the induction of k-turn RNA structure on binding L7Ae protein

The k-turn is a commonly occurring structural motif that introduces a tight kink into duplex RNA. In free solution, it can exist in an extended form, or by folding into the kinked structure. Binding of proteins including the L7Ae family can induce the formation of the kinked geometry, raising the question of whether this occurs by passive selection of the kinked structure, or a more active process in which the protein manipulates the RNA structure. We have devised a single-molecule experiment whereby immobilized L7Ae protein binds Cy3-Cy5-labeled RNA from free solution. We find that all bound RNA is in the kinked geometry, with no evidence for transitions to an extended form at the millisecond timescale of the camera. Furthermore, real-time binding experiments provide no evidence for a more extended intermediate even at the earliest times, at a time resolution of 16 ms. The data support a passive conformational selection model by which the protein selects a fraction of RNA that is already in the kinked conformation, thereby drawing the equilibrium into this form.

Probing nucleic acid interactions and pre-mRNA splicing by Förster Resonance Energy Transfer (FRET) microscopy

Förster resonance energy transfer (FRET) microscopy is a powerful technique routinely used to monitor interactions between biomolecules. Here, we focus on the techniques that are used for investigating the structure and interactions of nucleic acids (NAs). We present a brief overview of the most commonly used FRET microscopy techniques, their advantages and drawbacks. We list experimental approaches recently used for either in vitro or in vivo studies. Next, we summarize how FRET contributed to the understanding of pre-mRNA splicing and spliceosome assembly.

The structure and folding of kink turns in RNA

The kink turn (k-turn) is a widespread structural motif that introduces a tight kink into the axis of double-stranded RNA, with an included angle ∼60°. A standard k-turn comprises a three-nucleotide bulge followed on the 3' side by a G•A pair, an A•G pair, and usually further non-Watson-Crick pairs. The kinked conformation may be stabilized by three processes. These are the addition of metal ions, the binding of proteins such as the L7Ae family, and by the formation of tertiary interactions. The structure is characterized by specific A-minor interactions with the adenine nucleobases of the G•A pairs, and some very well-conserved hydrogen bonds involving 2'-hydroxyl groups. We can identify two classes of k-turns, that differ in the manner of the hydrogen bonding at the adenine of the bulge-distal G•A pair.

Structure and folding of a rare, natural kink turn in RNA with an A*A pair at the 2b*2n position

The kink turn (k-turn) is a frequently occurring motif, comprising a bulge followed by G•A and A•G pairs that introduces a sharp axial bend in duplex RNA. Natural k-turn sequences exhibit significant departures from the consensus, including the A•G pairs that form critical interactions stabilizing the core of the structure. Kt-23 found in the small ribosomal subunit differs from the consensus in many organisms, particularly in the second A•G pair distal to the bulge (2b•2n). Analysis of many Kt-23 sequences shows that the frequency of occurrence at the 2n position (i.e., on the nonbulged strand, normally G in standard k-turns) is U>C>G>A. Less than 1% of sequences have A at the 2n position, but one such example occurs in Thelohania solenopsae Kt-23. This sequence folds only weakly in the presence of Mg²⁺ ions but is induced to fold normally by the binding of L7Ae protein. Introduction of this sequence into the SAM-I riboswitch resulted in normal binding of SAM ligand, indicating that tertiary RNA contacts have resulted in k-turn folding. X-ray crystallography shows that the T. solenopsae Kt-23 adopts a standard k-turn geometry, making the key, conserved hydrogen bonds in the core and orienting the 1n (of the bulge-proximal A•G pair) and 2b adenine nucleobases in position facing the opposing minor groove. The 2b and 2n adenine nucleobases are not directly hydrogen bonded, but each makes hydrogen bonds to their opposing strands.

RNA tertiary interactions in a riboswitch stabilize the structure of a kink turn

The kink turn is a widespread RNA motif that introduces an acute kink into the axis of duplex RNA, typically comprising a bulge followed by a G⋅A and A⋅G pairs. The kinked conformation is stabilized by metal ions, or the binding of proteins including L7Ae. We now demonstrate a third mechanism for the stabilization of k-turn structure, involving tertiary interactions within a larger RNA structure. The SAM-I riboswitch contains an essential standard k-turn sequence that kinks a helix so that its terminal loop can make a long-range interaction. We find that some sequence variations in the k-turn within the riboswitch do not prevent SAM binding, despite preventing the folding of the k-turn in isolation. Furthermore, two crystal structures show that the sequence-variant k-turns are conventionally folded within the riboswitch. This study shows that the folded structure of the k-turn can be stabilized by tertiary interactions within a larger RNA structure.

Understanding RNA Flexibility Using Explicit Solvent Simulations: The Ribosomal and Group I Intron Reverse Kink-Turn Motifs

Reverse kink-turn is a recurrent elbow-like RNA building block occurring in the ribosome and in the group I intron. Its sequence signature almost matches that of the conventional kink-turn. However, the reverse and conventional kink-turns have opposite directions of bending. The reverse kink-turn lacks basically any tertiary interaction between its stems. We report unrestrained, explicit solvent molecular dynamics simulations of ribosomal and intron reverse kink-turns (54 simulations with 7.4 μs of data in total) with different variants (ff94, ff99, ff99bsc0, ff99χOL, and ff99bsc0χOL) of the Cornell et al. force field. We test several ion conditions and two water models. The simulations characterize the directional intrinsic flexibility of reverse kink-turns pertinent to their folded functional geometries. The reverse kink-turns are the most flexible RNA motifs studied so far by explicit solvent simulations which are capable at the present simulation time scale to spontaneously and reversibly sample a wide range of geometries from tightly kinked ones through flexible intermediates up to extended, unkinked structures. A possible biochemical role of the flexibility is discussed. Among the tested force fields, the latest χOL variant is essential to obtaining stable trajectories while all force field versions lacking the χ correction are prone to a swift degradation toward senseless ladder-like structures of stems, characterized by high-anti glycosidic torsions. The type of explicit water model affects the simulations considerably more than concentration and the type of ions.

Biological mechanisms, one molecule at a time

The last 15 years have witnessed the development of tools that allow the observation and manipulation of single molecules. The rapidly expanding application of these technologies for investigating biological systems of ever-increasing complexity is revolutionizing our ability to probe the mechanisms of biological reactions. Here, we compare the mechanistic information available from single-molecule experiments with the information typically obtained from ensemble studies and show how these two experimental approaches interface with each other. We next present a basic overview of the toolkit for observing and manipulating biology one molecule at a time. We close by presenting a case study demonstrating the impact that single-molecule approaches have had on our understanding of one of life's most fundamental biochemical reactions: the translation of a messenger RNA into its encoded protein by the ribosome.

Structural heterogeneity and quantitative FRET efficiency distributions of polyprolines through a hybrid atomistic simulation and Monte Carlo approach

Förster Resonance Energy Transfer (FRET) experiments probe molecular distances via distance dependent energy transfer from an excited donor dye to an acceptor dye. Single molecule experiments not only probe average distances, but also distance distributions or even fluctuations, and thus provide a powerful tool to study biomolecular structure and dynamics. However, the measured energy transfer efficiency depends not only on the distance between the dyes, but also on their mutual orientation, which is typically inaccessible to experiments. Thus, assumptions on the orientation distributions and averages are usually made, limiting the accuracy of the distance distributions extracted from FRET experiments. Here, we demonstrate that by combining single molecule FRET experiments with the mutual dye orientation statistics obtained from Molecular Dynamics (MD) simulations, improved estimates of distances and distributions are obtained. From the simulated time-dependent mutual orientations, FRET efficiencies are calculated and the full statistics of individual photon absorption, energy transfer, and photon emission events is obtained from subsequent Monte Carlo (MC) simulations of the FRET kinetics. All recorded emission events are collected to bursts from which efficiency distributions are calculated in close resemblance to the actual FRET experiment, taking shot noise fully into account. Using polyproline chains with attached Alexa 488 and Alexa 594 dyes as a test system, we demonstrate the feasibility of this approach by direct comparison to experimental data. We identified cis-isomers and different static local environments as sources of the experimentally observed heterogeneity. Reconstructions of distance distributions from experimental data at different levels of theory demonstrate how the respective underlying assumptions and approximations affect the obtained accuracy. Our results show that dye fluctuations obtained from MD simulations, combined with MC single photon kinetics, provide a versatile tool to improve the accuracy of distance distributions that can be extracted from measured single molecule FRET efficiencies.

Structural dynamics of the box C/D RNA kink-turn and its complex with proteins: the role of the A-minor 0 interaction, long-residency water bridges, and structural ion-binding sites revealed by molecular simulations

Kink-turns (K-turns) are recurrent elbow-like RNA motifs that participate in protein-assisted RNA folding and contribute to RNA dynamics. We carried out a set of molecular dynamics (MD) simulations using parm99 and parmbsc0 force fields to investigate structural dynamics of the box C/D RNA and its complexes with two proteins: native archaeal L7ae protein and human 15.5 kDa protein, originally bound to very similar structure of U4 snRNA. The box C/D RNA forms K-turn with A-minor 0 tertiary interaction between its canonical (C) and noncanonical (NC) stems. The local K-turn architecture is thus different from the previously studied ribosomal K-turns 38 and 42 having A-minor I interaction. The simulations reveal visible structural dynamics of this tertiary interaction involving altogether six substates which substantially contribute to the elbow-like flexibility of the K-turn. The interaction can even temporarily shift to the A-minor I type pattern; however, this is associated with distortion of the G/A base pair in the NC-stem of the K-turn. The simulations show reduction of the K-turn flexibility upon protein binding. The protein interacts with the apex of the K-turn and with the NC-stem. The protein-RNA interface includes long-residency hydration sites. We have also found long-residency hydration sites and major ion-binding sites associated with the K-turn itself. The overall topology of the K-turn remains stable in all simulations. However, in simulations of free K-turn, we observed instability of the key C16(O2')-A7(N1) H-bond, which is a signature interaction of K-turns and which was visibly more stable in simulations of K-turns possessing A-minor I interaction. It may reflect either some imbalance of the force field or it may be a correct indication of early stages of unfolding since this K-turn requires protein binding for its stabilization. Interestingly, the 16(O2')-7(N1) H- bond is usually not fully lost since it is replaced by a water bridge with a tightly bound water, which is adenine-specific similarly as the original interaction. The 16(O2')-7(N1) H-bond is stabilized by protein binding. The stabilizing effect is more visible with the human 15.5 kDa protein, which is attributed to valine to arginine substitution in the binding site. The behavior of the A-minor interaction is force-field-dependent because the parmbsc0 force field attenuates the A-minor fluctuations compared to parm99 simulations. Behavior of other regions of the box C/D RNA is not sensitive to the force field choice. Simulation with net-neutralizing Na(+) and 0.2 M excess salt conditions appear in all aspects equivalent. The simulations show loss of a hairpin tetraloop, which is not part of the K-turn. This was attributed to force field limitations.

A structural database for k-turn motifs in RNA

The kink-turn (k-turn) is a common structural motif in RNA that introduces a tight kink into the helical axis. k-turns play an important architectural role in RNA structures and serve as binding sites for a number of proteins. We have created a database of known and postulated k-turn sequences and three-dimensional (3D) structures, available via the internet. This site provides (1) a database of sequence and structure, as a resource for the RNA community, and (2) a tool to enable the manipulation and comparison of 3D structures where known.

Ion-induced folding of a kink turn that departs from the conventional sequence

Kink turns (k-turns) are important structural motifs that create a sharp axial bend in RNA. Most conform to a consensus in which a three-nucleotide bulge is followed by consecutive G*A and A*G base pairs, and when these G*A pairs are modified in vitro this generally leads to a failure to adopt the k-turn conformation. Kt-23 in the 30S ribosomal subunit of Thermus thermophilus is a rare exception in which the bulge-distal A*G pair is replaced by a non-Watson-Crick A*U pair. In the context of the ribosome, Kt-23 adopts a completely conventional k-turn geometry. We show here that this sequence is induced to fold into a k-turn structure in an isolated RNA duplex by Mg(2+) or Na(+) ions. Therefore, the Kt-23 is intrinsically stable despite lacking the key A*G pair; its formation requires neither tertiary interactions nor protein binding. Moreover, the Kt-23 k-turn is stabilized by the same critical hydrogen-bonding interactions within the core of the structure that are found in more conventional sequences such as the near-consensus Kt-7. T. thermophilus Kt-23 has two further non-Watson-Crick base pairs within the non-canonical helix, three and four nucleotides from the bulge, and we find that the nature of these pairs influences the ability of the RNA to adopt k-turn conformation, although the base pair adjacent to the A*U pair is more important than the other.

Single-molecule Förster resonance energy transfer studies of RNA structure, dynamics and function

Single-molecule fluorescence microscopy experiments on RNA molecules brought to light the highly complex dynamics of key biological processes, including RNA folding, catalysis of ribozymes, ligand sensing of riboswitches and aptamers, and protein synthesis in the ribosome. By using highly advanced biophysical spectroscopy techniques in combination with sophisticated biochemical synthesis approaches, molecular dynamics of individual RNA molecules can be observed in real time and under physiological conditions in unprecedented detail that cannot be achieved with bulk experiments. Here, we review recent advances in RNA folding and functional studies of RNA and RNA-protein complexes addressed by using single-molecule Förster (fluorescence) resonance energy transfer (smFRET) technique.

Single-molecule FRET measures bends and kinks in DNA

We present advances in the use of single-molecule FRET measurements with flexibly linked dyes to derive full 3D structures of DNA constructs based on absolute distances. The resolution obtained by this single-molecule approach harbours the potential to study in detail also protein- or damage-induced DNA bending. If one is to generate a geometric structural model, distances between fixed positions are needed. These are usually not experimentally accessible because of unknown fluorophore-linker mobility effects that lead to a distribution of FRET efficiencies and distances. To solve this problem, we performed studies on DNA double-helices by systematically varying donor acceptor distances from 2 to 10 nm. Analysis of dye-dye quenching and fluorescence anisotropy measurements reveal slow positional and fast orientational fluorophore dynamics, that results in an isotropic average of the FRET efficiency. We use a nonlinear conversion function based on MD simulations that allows us to include this effect in the calculation of absolute FRET distances. To obtain unique structures, we performed a quantitative statistical analysis for the conformational search in full space based on triangulation, which uses the known helical nucleic acid features. Our higher accuracy allowed the detection of sequence-dependent DNA bending by 16 degrees . For DNA with bulged adenosines, we also quantified the kink angles introduced by the insertion of 1, 3 and 5 bases to be 32 degrees +/- 6 degrees , 56 degrees +/- 4 degrees and 73 +/- 2 degrees , respectively. Moreover, the rotation angles and shifts of the helices were calculated to describe the relative orientation of the two arms in detail.

The importance of G.A hydrogen bonding in the metal ion- and protein-induced folding of a kink turn RNA

The kink turn (K-turn) is a common motif in RNA structure, found in many RNA species important in translation, RNA modification and splicing, and the control of gene expression. In general the K-turn comprises a three nucleotide bulge followed by trans sugar-Hoogsteen G.A pairs. The RNA adopts a tightly kinked conformation, and is a common target for binding proteins, exemplified by the L7Ae family. We have measured the rates of association and dissociation for the binding of L7Ae to the Kt-7 kink turn, from which we calculate an affinity of K(D)=10 pM. This high affinity is consistent with the role of this binding as the first stage in the assembly of key functional nucleoproteins such as box C/D snoRNP. Kink-turn RNA undergoes a two-state transition between the kinked conformation, and a more extended structure, and folding into the kinked form is induced by divalent metal ions, or by binding of proteins of the L7Ae class. The K-turn provides an excellent, simple model for RNA folding, which can be dissected at the atomic level. We have analyzed the contributions of the hydrogen bonds that form the G.A pairs to the ion- and protein-induced folding of the K-turn. We find that all four hydrogen bonds are important to the stability of the kinked form of the RNA, and we can now define all the important hydrogen bonding interactions that stabilize the K-turn. The high affinity of L7Ae binding is coupled to the induced folding of the K-turn, allowing some sub-optimal variants to adopt the kinked geometry. However, in all such cases the affinity is lowered, and the results underline the importance of both G.A pairs to the stability of the K-turn.

Fluorescence from diffusing single molecules illuminates biomolecular structure and dynamics

This review provides an account of single-molecule fluorescence methodologies for freely diffusing molecules applied to a diverse array of biological problems pertaining to biomolecular folding and assembly. We describe the principles of confocal fluorescence microscopy to detect and analyze fluorescence bursts from diffusing single molecules. These methods including single-molecule fluorescence resonance energy transfer, coincidence, correlations and polarization offer a powerful means to uncover hidden information about conformational sub-populations and interconversion dynamics of biomolecular systems in a wide range of timescales. We offer several key examples to illustrate how these methodologies have been extremely useful in teasing out structural and dynamical aspects of many important biomolecular systems.

Staufen1 regulates diverse classes of mammalian transcripts

It is currently unknown how extensively the double-stranded RNA-binding protein Staufen (Stau)1 is utilized by mammalian cells to regulate gene expression. To date, Stau1 binding to the 3'-untranslated region (3'-UTR) of ADP ribosylation factor (ARF)1 mRNA has been shown to target ARF1 mRNA for Stau1-mediated mRNA decay (SMD). ARF1 SMD depends on translation and recruitment of the nonsense-mediated mRNA decay factor Upf1 to the ARF1 3'-UTR by Stau1. Here, we demonstrate that Stau1 binds to a complex structure within the ARF1 3'-UTR. We also use microarrays to show that 1.1 and 1.0% of the 11 569 HeLa-cell transcripts that were analyzed are upregulated and downregulated, respectively, at least two-fold upon Stau1 depletion in three independently performed experiments. We localize the Stau1 binding site to the 3'-UTR of four mRNAs that we define as natural SMD targets. Additionally, we provide evidence that the efficiency of SMD increases during the differentiation of C2C12 myoblasts to myotubes. We propose that Stau1 influences the expression of a wide variety of physiologic transcripts and metabolic pathways.

An improved definition of the RNA-binding specificity of SECIS-binding protein 2, an essential component of the selenocysteine incorporation machinery

By binding to SECIS elements located in the 3′-UTR of selenoprotein mRNAs, the protein SBP2 plays a key role in the assembly of the selenocysteine incorporation machinery. SBP2 contains an L7Ae/L30 RNA-binding domain similar to that of protein 15.5K/Snu13p, which binds K-turn motifs with a 3-nt bulge loop closed by a tandem of G.A and A.G pairs. Here, by SELEX experiments, we demonstrate the capacity of SBP2 to bind such K-turn motifs with a protruding U residue. However, we show that conversion of the bulge loop into an internal loop reinforces SBP2 affinity and to a greater extent RNP stability. Opposite variations were found for Snu13p. Accordingly, footprinting assays revealed strong contacts of SBP2 with helices I and II and the 5′-strand of the internal loop, as opposed to the loose interaction of Snu13p. Our data also identifies new determinants for SBP2 binding which are located in helix II. Among the L7Ae/L30 family members, these determinants are unique to SBP2. Finally, in accordance with functional data on SECIS elements, the identity of residues at positions 2 and 3 in the loop influences SBP2 affinity. Altogether, the data provide a very precise definition of the SBP2 RNA specificity.

The role of specific 2'-hydroxyl groups in the stabilization of the folded conformation of kink-turn RNA

The role of 2'-hydroxyl groups in stabilizing the tightly kinked geometry of the kink-turn (K-turn) has been investigated. Individual 2'-OH groups have been removed by chemical synthesis, and the kinking of the RNA has been studied by gel electrophoresis and fluorescence resonance energy transfer. The results have been analyzed by reference to a database of 11 different crystallographic structures of K-turns. The potential hydrogen bonds fall into several classes. The most important are those in the core of the turn and ribose-phosphate interactions around the bulge. Of these the single most important hydrogen bond is one donated from the 2'-OH of the 5' nucleotide of the bulge to the N1 of the adenine of the kink-proximal A*G pair. This is present in all known K-turn structures, and removal of the 2'-OH completely prevents metal ion-induced folding. Hydrogen bonds formed in the minor grooves of the helical stems are less important, and removal of the participating 2'-OH groups leads to reduced impairment of folding. These interactions are generally more polymorphic, and hydrogen bonds probably form where possible, as permitted by the global structure.

Analysis of sequence and structural features that identify the B/C motif of U3 small nucleolar RNA as the recognition site for the Snu13p-Rrp9p protein pair

The eukaryal Snu13p/15.5K protein binds K-turn motifs in U4 snRNA and snoRNAs. Two Snu13p/15.5K molecules bind the nucleolar U3 snoRNA required for the early steps of preribosomal processing. Binding of one molecule on the C'/D motif allows association of proteins Nop1p, Nop56p, and Nop58p, whereas binding of the second molecule on the B/C motif allows Rrp9p recruitment. To understand how the Snu13p-Rrp9p pair recognizes the B/C motif, we first improved the identification of RNA determinants required for Snu13p binding by experiments using the systematic evolution of ligands by exponential enrichment. This demonstrated the importance of a U.U pair stacked on the sheared pairs and revealed a direct link between Snu13p affinity and the stability of helices I and II. Sequence and structure requirements for efficient association of Rrp9p on the B/C motif were studied in yeast cells by expression of variant U3 snoRNAs and immunoselection assays. A G-C pair in stem II, a G residue at position 1 in the bulge, and a short stem I were found to be required. The data identify the in vivo function of most of the conserved residues of the U3 snoRNA B/C motif. They bring important information to understand how different K-turn motifs can recruit different sets of proteins after Snu13p association.

mRNPs take shape by CLIPPING and PAIRING

The interaction of RNA-binding proteins (RBPs) with RNA is a crucial aspect of normal cellular metabolism. Yet, the diverse number of RBPs and RNA motifs to which they bind, the wide range of interaction strengths and the fact that RBPs associate in dynamic complexes have made it challenging to determine whether a particular RNA-binding protein binds a particular RNA. Recent work by three different laboratories has led to the development of new tools to query such interactions in the more physiological environs of cultured cells. The use of these methods has led to insights into (1) the networks of RNAs regulated by a particular protein, (2) the identification of new protein partners within messenger ribonucleoprotein particles and (3) the flux of RNA-binding proteins on an mRNA throughout its lifecycle. Here, I examine these new methods and discuss their relative strengths and current limitations.

Assembly of the archaeal box C/D sRNP can occur via alternative pathways and requires temperature-facilitated sRNA remodeling

Archaeal dual-guide box C/D small nucleolar RNA-like RNAs (sRNAs) bind three core proteins in sequential order at both terminal box C/D and internal C'/D' motifs to assemble two ribonuclear protein (RNP) complexes active in guiding nucleotide methylation. Experiments have investigated the process of box C/D sRNP assembly and the resultant changes in sRNA structure or "remodeling" as a consequence of sRNP core protein binding. Hierarchical assembly of the Methanocaldococcus jannaschii sR8 box C/D sRNP is a temperature-dependent process with binding of L7 and Nop56/58 core proteins to the sRNA requiring elevated temperature to facilitate necessary RNA structural dynamics. Circular dichroism (CD) spectroscopy and RNA thermal denaturation revealed an increased order and stability of sRNA folded structure as a result of L7 binding. Subsequent binding of the Nop56/58 and fibrillarin core proteins to the L7-sRNA complex further remodeled sRNA structure. Assessment of sR8 guide region accessibility using complementary RNA oligonucleotide probes revealed significant changes in guide region structure during sRNP assembly. A second dual-guide box C/D sRNA from M. jannaschii, sR6, also exhibited RNA remodeling during temperature-dependent sRNP assembly, although core protein binding was affected by sR6's distinct folded structure. Interestingly, the sR6 sRNP followed an alternative assembly pathway, with both guide regions being continuously exposed during sRNP assembly. Further experiments using sR8 mutants possessing alternative guide regions demonstrated that sRNA folded structure induced by specific guide sequences impacted the sRNP assembly pathway. Nevertheless, assembled sRNPs were active for sRNA-guided methylation independent of the pathway followed. Thus, RNA remodeling appears to be a common and requisite feature of archaeal dual-guide box C/D sRNP assembly and formation of the mature sRNP can follow different assembly pathways in generating catalytically active complexes.

Protein-protein and protein-RNA contacts both contribute to the 15.5K-mediated assembly of the U4/U6 snRNP and the box C/D snoRNPs

The k-turn-binding protein 15.5K is unique in that it is essential for the hierarchical assembly of three RNP complexes distinct in both composition and function, namely, the U4/U6 snRNP, the box C/D snoRNP, and the RNP complex assembled on the U3 box B/C motif. 15.5K interacts with the cognate RNAs via an induced fit mechanism, which results in the folding of the surrounding RNA to create a binding site(s) for the RNP-specific proteins. However, it is possible that 15.5K also mediates RNP formation via protein-protein interactions with the complex-specific proteins. To investigate this possibility, we created a series of 15.5K mutations in which the surface properties of the protein had been changed. We assessed their ability to support the formation of the three distinct RNP complexes and found that the formation of each RNP requires a distinct set of regions on the surface of 15.5K. This implies that protein-protein contacts are essential for RNP formation in each complex. Further supporting this idea, direct protein-protein interaction could be observed between hU3-55K and 15.5K. In conclusion, our data suggest that the formation of each RNP involves the direct recognition of specific elements in both 15.5K protein and the specific RNA.