NMR structure of the 5' splice site in the group IIB intron Sc.ai5γ--conformational requirements for exon-intron recognition

Multipolar electrostatics for hairpin and pseudoknots in RNA: Improving the accuracy of force field potential energy function

Molecular dynamics (MD) simulations that rely on force field methods has been widely used to explore the structure and function of RNAs. However, the current commonly used force fields are limited by the electrostatic description offered by atomic charge, dipole and at most quadrupole moments, failing to capture the anisotropic picture of electronic features. Actually, the distribution of electrons around atomic nuclei is not spherically symmetric but is geometry dependent. A multipolar electrostatic model based on high rank multipole moments is described in this work, which allows us to combine polarizability and anisotropy of electron density. RNA secondary structure was taken as a research system, and its substructures including stem, loops (hairpin loop, bulge loop, internal loop, and multi-branch loop), and pseudoknots (H-type and K-type) were investigated, respectively, as well as the hairpin. First, the atom-atom electrostatic properties derived from one chain of a duplex RNA 2MVY in our previous work (Ref. 58) were measured by the pilot RNA systems of hairpin, hairpin loop, stem, and H-type pseudoknot, respectively. The prediction results were not satisfactory. Consequently, to obtain a general set of electrostatic parameters for RNA force fields, the convergence behavior of the atom-atom electrostatic interactions in the pilot RNA systems was explored using high rank atomic multipole moments. The pilot RNA systems were cut into four types of different-sized molecular fragments, and the single nucleotide fragment and nucleotide-paired fragment proved to be the most reasonable systems for base-unpairing regions and base-pairing regions to investigate the convergence behavior of all types of atom-atom electrostatic interactions, respectively. Transferability of the electrostatic properties drawn from the pilot RNA systems to the corresponding test systems was also investigated. Furthermore, the convergence behavior of atomic electrostatic interactions in other substructures including bulge loop, internal loop, multi-branch loop, and K-type pseudoknot was expected to be modeled via the hairpin.

Metal ions and sugar puckering balance single-molecule kinetic heterogeneity in RNA and DNA tertiary contacts

The fidelity of group II intron self-splicing and retrohoming relies on long-range tertiary interactions between the intron and its flanking exons. By single-molecule FRET, we explore the binding kinetics of the most important, structurally conserved contact, the exon and intron binding site 1 (EBS1/IBS1). A comparison of RNA-RNA and RNA-DNA hybrid contacts identifies transient metal ion binding as a major source of kinetic heterogeneity which typically appears in the form of degenerate FRET states. Molecular dynamics simulations suggest a structural link between heterogeneity and the sugar conformation at the exon-intron binding interface. While Mg2+ ions lock the exon in place and give rise to long dwell times in the exon bound FRET state, sugar puckering alleviates this structural rigidity and likely promotes exon release. The interplay of sugar puckering and metal ion coordination may be an important mechanism to balance binding affinities of RNA and DNA interactions in general.

Contribution of Mobile Group II Introns to Sinorhizobium meliloti Genome Evolution

Mobile group II introns are ribozymes and retroelements that probably originate from bacteria. Sinorhizobium meliloti, the nitrogen-fixing endosymbiont of legumes of genus Medicago, harbors a large number of these retroelements. One of these elements, RmInt1, has been particularly successful at colonizing this multipartite genome. Many studies have improved our understanding of RmInt1 and phylogenetically related group II introns, their mobility mechanisms, spread and dynamics within S. meliloti and closely related species. Although RmInt1 conserves the ancient retroelement behavior, its evolutionary history suggests that this group II intron has played a role in the short- and long-term evolution of the S. meliloti genome. We will discuss its proposed role in genome evolution by controlling the spread and coexistence of potentially harmful mobile genetic elements, by ectopic transposition to different genetic loci as a source of early genomic variation and by generating sequence variation after a very slow degradation process, through intron remnants that may have continued to evolve, contributing to bacterial speciation.

An atomistic view on carbocyanine photophysics in the realm of RNA

Carbocyanine dyes have a long-standing tradition in fluorescence imaging and spectroscopy, due to their photostability and large spectral separation between individual dye species. Herein, we explore the versatility of cyanine dyes to probe the dynamics of nucleic acids and we report on the interrelation of fluorophores, RNA, and metal ions, namely K⁺ and Mg²⁺. Photophysical parameters including the fluorescence lifetime, quantum yield and dynamic anisotropy are monitored as a function of the nucleic acid composition, conformation, and metal ion abundance. Occasional excursions to a non-fluorescent cis-state hint at the remarkable sensitivity of carbocyanines to their local environment. Comparison of time-correlated single photon experiments with all-atom molecular dynamics simulations demonstrate that the propensity of photoisomerization is dictated by sterical constraints imposed on the fluorophore. Structural features in the vicinity of the dye play a crucial role in RNA recognition and have far-reaching implications on the mobility of the fluorescent probe. An atomic level description of the mutual interactions will ultimately benefit the quantitative interpretation of single-molecule FRET measurements on large RNA systems.

NMR Studies of the Structure and Function of the HIV-1 5'-Leader

The 5′-leader of the human immunodeficiency virus type 1 (HIV-1) genome plays several critical roles during viral replication, including differentially establishing mRNA versus genomic RNA (gRNA) fates. As observed for proteins, the function of the RNA is tightly regulated by its structure, and a common paradigm has been that genome function is temporally modulated by structural changes in the 5′-leader. Over the past 30 years, combinations of nucleotide reactivity mapping experiments with biochemistry, mutagenesis, and phylogenetic studies have provided clues regarding the secondary structures of stretches of residues within the leader that adopt functionally discrete domains. More recently, nuclear magnetic resonance (NMR) spectroscopy approaches have been developed that enable direct detection of intra- and inter-molecular interactions within the intact leader, providing detailed insights into the structural determinants and mechanisms that regulate HIV-1 genome packaging and function.

Improving NMR Structures of RNA

Here, we show that modern solution nuclear magnetic resonance (NMR) structures of RNA exhibit more steric clashes and conformational ambiguities than their crystallographic X-ray counterparts. To tackle these issues, we developed RNA-ff1, a new force field for structure calculation with Xplor-NIH. Using seven published NMR datasets, RNA-ff1 improves covalent geometry and MolProbity validation criteria for clashes and backbone conformation in most cases, relative to both the previous Xplor-NIH force field and the original structures associated with the experimental data. In addition, with smaller base-pair step rises in helical stems, RNA-ff1 structures enjoy more favorable base stacking. Finally, structural accuracy improves in the majority of cases, as supported by complete residual dipolar coupling cross-validation. Thus, the reported advances show great promise in bridging the quality gap that separates NMR and X-ray structures of RNA.

Structural features of a 3' splice site in influenza a

Influenza A is an RNA virus with a genome of eight negative sense segments. Segment 7 mRNA contains a 3′ splice site for alternative splicing to encode the essential M2 protein. On the basis of sequence alignment and chemical mapping experiments, the secondary structure surrounding the 3′ splice site has an internal loop, adenine bulge, and hairpin loop when it is in the hairpin conformation that exposes the 3′ splice site. We report structural features of a three-dimensional model of the hairpin derived from nuclear magnetic resonance spectra and simulated annealing with restrained molecular dynamics. Additional insight was provided by modeling based on ¹H chemical shifts. The internal loop containing the 3′ splice site has a dynamic guanosine and a stable imino (cis Watson–Crick/Watson–Crick) GA pair. The adenine bulge also appears to be dynamic with the A either stacked in the stem or forming a base triple with a Watson–Crick GC pair. The hairpin loop is a GAAA tetraloop closed by an AC pair.

Cation-induced kinetic heterogeneity of the intron-exon recognition in single group II introns

RNA is commonly believed to undergo a number of sequential folding steps before reaching its functional fold, i.e., the global minimum in the free energy landscape. However, there is accumulating evidence that several functional conformations are often in coexistence, corresponding to multiple (local) minima in the folding landscape. Here we use the 5'-exon-intron recognition duplex of a self-splicing ribozyme as a model system to study the influence of Mg(2+) and Ca(2+) on RNA tertiary structure formation. Bulk and single-molecule spectroscopy reveal that near-physiological M(2+) concentrations strongly promote interstrand association. Moreover, the presence of M(2+) leads to pronounced kinetic heterogeneity, suggesting the coexistence of multiple docked and undocked RNA conformations. Heterogeneity is found to decrease at saturating M(2+) concentrations. Using NMR, we locate specific Mg(2+) binding pockets and quantify their affinity toward Mg(2+). Mg(2+) pulse experiments show that M(2+) exchange occurs on the timescale of seconds. This unprecedented combination of NMR and single-molecule Förster resonance energy transfer demonstrates for the first time to our knowledge that a rugged free energy landscape coincides with incomplete occupation of specific M(2+) binding sites at near-physiological M(2+) concentrations. Unconventional kinetics in nucleic acid folding frequently encountered in single-molecule experiments are therefore likely to originate from a spectrum of conformations that differ in the occupation of M(2+) binding sites.

The role of Mg(II) in DNA cleavage site recognition in group II intron ribozymes: solution structure and metal ion binding sites of the RNA-DNA complex

Group II intron ribozymes catalyze the cleavage of (and their reinsertion into) DNA and RNA targets using a Mg2(+)-dependent reaction. The target is cleaved 3' to the last nucleotide of intron binding site 1 (IBS1), one of three regions that form base pairs with the intron's exon binding sites (EBS1 to -3).We solved the NMR solution structure of the d3' hairpin of the Sc.ai5γ intron containing EBS1 in its 11-nucleotide loop in complex with the dIBS1 DNA 7-mer and compare it with the analogous RNA-RNA contact. The EBS1-dIBS1 helix is slightly flexible and non-symmetric. NMR data reveal two major groove binding sites for divalent metal ions at the EBS1-dIBS1 helix, and surface plasmon resonance experiments show that low concentrations of Mg2(+) considerably enhance the affinity of dIBS1 for EBS1. Our results indicate that identification of both RNA and DNA IBS1 targets, presentation of the scissile bond, and stabilization of the structure by metal ions are governed by the overall structure of EBS1-dIBS1 and the surrounding loop nucleotides but are irrespective of different EBS1-(d)IBS1 geometries and interstrand affinities.

BOBA FRET: bootstrap-based analysis of single-molecule FRET data

Time-binned single-molecule Förster resonance energy transfer (smFRET) experiments with surface-tethered nucleic acids or proteins permit to follow folding and catalysis of single molecules in real-time. Due to the intrinsically low signal-to-noise ratio (SNR) in smFRET time traces, research over the past years has focused on the development of new methods to extract discrete states (conformations) from noisy data. However, limited observation time typically leads to pronounced cross-sample variability, i.e., single molecules display differences in the relative population of states and the corresponding conversion rates. Quantification of cross-sample variability is necessary to perform statistical testing in order to assess whether changes observed in response to an experimental parameter (metal ion concentration, the presence of a ligand, etc.) are significant. However, such hypothesis testing has been disregarded to date, precluding robust biological interpretation. Here, we address this problem by a bootstrap-based approach to estimate the experimental variability. Simulated time traces are presented to assess the robustness of the algorithm in conjunction with approaches commonly used in thermodynamic and kinetic analysis of time-binned smFRET data. Furthermore, a pair of functionally important sequences derived from the self-cleaving group II intron Sc.ai5γ (d3'EBS1*/IBS1*) is used as a model system. Through statistical hypothesis testing, divalent metal ions are shown to have a statistically significant effect on both thermodynamic and kinetic aspects of their interaction. The Matlab source code used for analysis (bootstrap-based analysis of smFRET data, BOBA FRET), as well as a graphical user interface, is available via http://www.aci.uzh.ch/rna/.