Structure Determination of Protein⧸RNA Complexes by NMR

An RNA-centric historical narrative around the Protein Data Bank

Some of the amazing contributions brought to the scientific community by the Protein Data Bank (PDB) are described. The focus is on nucleic acid structures with a bias toward RNA. The evolution and key roles in science of the PDB and other structural databases for nucleic acids illustrate how small initial ideas can become huge and indispensable resources with the unflinching willingness of scientists to cooperate globally. The progress in the understanding of the molecular interactions driving RNA architectures followed the rapid increase in RNA structures in the PDB. That increase was consecutive to improvements in chemical synthesis and purification of RNA molecules, as well as in biophysical methods for structure determination and computer technology. The RNA modeling efforts from the early beginnings are also described together with their links to the state of structural knowledge and technological development. Structures of RNA and of its assemblies are physical objects, which, together with genomic data, allow us to integrate present-day biological functions and the historical evolution in all living species on earth.

A solid-state NMR tool box for the investigation of ATP-fueled protein engines

Motor proteins are involved in a variety of cellular processes. Their main purpose is to convert the chemical energy released during adenosine triphosphate (ATP) hydrolysis into mechanical work. In this review, solid-state Nuclear Magnetic Resonance (NMR) approaches are discussed allowing studies of structures, conformational events and dynamic features of motor proteins during a variety of enzymatic reactions. Solid-state NMR benefits from straightforward sample preparation based on sedimentation of the proteins directly into the Magic-Angle Spinning (MAS) rotor. Protein resonance assignment is the crucial and often time-limiting step in interpreting the wealth of information encoded in the NMR spectra. Herein, potentials, challenges and limitations in resonance assignment for large motor proteins are presented, focussing on both biochemical and spectroscopic approaches. This work highlights NMR tools available to study the action of the motor domain and its coupling to functional processes, as well as to identify protein-nucleotide interactions during events such as DNA replication. Arrested protein states of reaction coordinates such as ATP hydrolysis can be trapped for NMR studies by using stable, non-hydrolysable ATP analogues that mimic the physiological relevant states as accurately as possible. Recent advances in solid-state NMR techniques ranging from Dynamic Nuclear Polarization (DNP), ³¹P-based heteronuclear correlation experiments, ¹H-detected spectra at fast MAS frequencies >100 kHz to paramagnetic NMR are summarized and their applications to the bacterial DnaB helicase from Helicobacter pylori are discussed.

Interaction proteomics by using in-cell NMR spectroscopy

A synopsis of in-cell NMR spectroscopic approaches to study interaction proteomics in prokaryotic and eukaryotic cells is presented. We describe the use of in-cell NMR spectroscopy to resolve high resolution protein structures, discuss methodologies for determining and analyzing high and low affinity protein-target structural interactions, including intrinsically disordered proteins, and detail important functional interactions that result from these interactions. SIGNIFICANCE: The ultimate goal of structural and biochemical research is to understand how macromolecular interactions give rise to and regulate biological activity in living cells. The challenge is formidable due to the complexity that arises not only from the number of proteins (genes) expressed by the organism, but also from the combinatorial interactions between them. Despite ongoing efforts to decipher the complex nature of protein interactions, new methods for structurally characterizing protein complexes are needed to fully understand molecular networks. With the onset of in-cell NMR spectroscopy, molecular structures and interactions can be studied under physiological conditions shedding light on the structural underpinning of biological activity.

Partner-specific prediction of RNA-binding residues in proteins: A critical assessment

RNA-protein interactions play essential roles in regulating gene expression. While some RNA-protein interactions are "specific", that is, the RNA-binding proteins preferentially bind to particular RNA sequence or structural motifs, others are "non-RNA specific." Deciphering the protein-RNA recognition code is essential for comprehending the functional implications of these interactions and for developing new therapies for many diseases. Because of the high cost of experimental determination of protein-RNA interfaces, there is a need for computational methods to identify RNA-binding residues in proteins. While most of the existing computational methods for predicting RNA-binding residues in RNA-binding proteins are oblivious to the characteristics of the partner RNA, there is growing interest in methods for partner-specific prediction of RNA binding sites in proteins. In this work, we assess the performance of two recently published partner-specific protein-RNA interface prediction tools, PS-PRIP, and PRIdictor, along with our own new tools. Specifically, we introduce a novel metric, RNA-specificity metric (RSM), for quantifying the RNA-specificity of the RNA binding residues predicted by such tools. Our results show that the RNA-binding residues predicted by previously published methods are oblivious to the characteristics of the putative RNA binding partner. Moreover, when evaluated using partner-agnostic metrics, RNA partner-specific methods are outperformed by the state-of-the-art partner-agnostic methods. We conjecture that either (a) the protein-RNA complexes in PDB are not representative of the protein-RNA interactions in nature, or (b) the current methods for partner-specific prediction of RNA-binding residues in proteins fail to account for the differences in RNA partner-specific versus partner-agnostic protein-RNA interactions, or both.

Ribosome Mediated Quinary Interactions Modulate In-Cell Protein Activities

Ribosomes are present inside bacterial cells at micromolar concentrations and occupy up to 20% of the cell volume. Under these conditions, even weak quinary interactions between ribosomes and cytosolic proteins can affect protein activity. By using in-cell and in vitro NMR spectroscopy, and biophysical techniques, we show that the enzymes, adenylate kinase and dihydrofolate reductase, and the respective coenzymes, ATP and NADPH, bind to ribosomes with micromolar affinity, and that this interaction suppresses the enzymatic activities of both enzymes. Conversely, thymidylate synthase, which works together with dihydrofolate reductase in the thymidylate synthetic pathway, is activated by ribosomes. We also show that ribosomes impede diffusion of green fluorescent protein in vitro and contribute to the decrease in diffusion in vivo. These results strongly suggest that ribosome-mediated quinary interactions contribute to the differences between in vitro and in vivo protein activities and that ribosomes play a previously under-appreciated nontranslational role in regulating cellular biochemistry.

RNA structure refinement using NMR solvent accessibility data

NMR spectroscopy is a powerful technique to study ribonucleic acids (RNAs) which are key players in a plethora of cellular processes. Although the NMR toolbox for structural studies of RNAs expanded during the last decades, they often remain challenging. Here, we show that solvent paramagnetic relaxation enhancements (sPRE) induced by the soluble, paramagnetic compound Gd(DTPA-BMA) provide a quantitative measure for RNA solvent accessibility and encode distance-to-surface information that correlates well with RNA structure and improves accuracy and convergence of RNA structure determination. Moreover, we show that sPRE data can be easily obtained for RNAs with any isotope labeling scheme and is advantageous regarding sample preparation, stability and recovery. sPRE data show a large dynamic range and reflect the global fold of the RNA suggesting that they are well suited to identify interaction surfaces, to score structural models and as restraints in RNA structure determination.

Sequence-Based Prediction of RNA-Binding Residues in Proteins

Identifying individual residues in the interfaces of protein-RNA complexes is important for understanding the molecular determinants of protein-RNA recognition and has many potential applications. Recent technical advances have led to several high-throughput experimental methods for identifying partners in protein-RNA complexes, but determining RNA-binding residues in proteins is still expensive and time-consuming. This chapter focuses on available computational methods for identifying which amino acids in an RNA-binding protein participate directly in contacting RNA. Step-by-step protocols for using three different web-based servers to predict RNA-binding residues are described. In addition, currently available web servers and software tools for predicting RNA-binding sites, as well as databases that contain valuable information about known protein-RNA complexes, RNA-binding motifs in proteins, and protein-binding recognition sites in RNA are provided. We emphasize sequence-based methods that can reliably identify interfacial residues without the requirement for structural information regarding either the RNA-binding protein or its RNA partner.

NMR Methods for Characterization of RNA Secondary Structure

Knowledge of RNA secondary structure is often sufficient to identify relationships between the structure of RNA and processing pathways, and the design of therapeutics. Nuclear magnetic resonance (NMR) can identify types of nucleotide base pairs and the sequence, thus limiting possible secondary structures. Because NMR experiments, like chemical mapping, are performed in solution, not in single crystals, experiments can be initiated as soon as the biomolecule is expressed and purified. This chapter summarizes NMR methods that permit rapid identification of RNA secondary structure, information that can be used as supplements to chemical mapping, and/or as preliminary steps required for 3D structure determination. The primary aim is to provide guidelines to enable a researcher with minimal knowledge of NMR to quickly extract secondary structure information from basic datasets. Instrumental and sample considerations that can maximize data quality are discussed along with some details for optimal data acquisition and processing parameters. Approaches for identifying base pair types in both unlabeled and isotopically labeled RNA are covered. Common problems, such as missing signals and overlaps, and approaches to address them are considered. Programs under development for merging NMR data with structure prediction algorithms are briefly discussed.

A fully enzymatic method for site-directed spin labeling of long RNA

Site-directed spin labeling is emerging as an essential tool to investigate the structural and dynamical features of RNA. We propose here an enzymatic method, which allows the insertion of a paramagnetic center at a specific position in an RNA molecule. The technique is based on a segmental approach using a ligation protocol with T4 RNA ligase 2. One transcribed acceptor RNA is ligated to a donor RNA in which a thio-modified nucleotide is introduced at its 5′-end by in vitro transcription with T7 RNA polymerase. The paramagnetic thiol-specific reagent is subsequently attached to the RNA ligation product. This novel strategy is demonstrated by introducing a paramagnetic probe into the 55 nucleotides long RNA corresponding to K-turn and Specifier Loop domains from the Bacillus subtilis tyrS T-Box leader RNA. The efficiency of the coupling reaction and the quality of the resulting spin-labeled RNA were assessed by Mass Spectrometry, Electron Paramagnetic Resonance (EPR) and Nuclear Magnetic Resonance (NMR). This method enables various combinations of isotopic segmental labeling and spin labeling schemes, a strategy that will be of particular interest to investigate the structural and dynamical properties of large RNA complexes by NMR and EPR spectroscopies.

RNABindRPlus: a predictor that combines machine learning and sequence homology-based methods to improve the reliability of predicted RNA-binding residues in proteins

Protein-RNA interactions are central to essential cellular processes such as protein synthesis and regulation of gene expression and play roles in human infectious and genetic diseases. Reliable identification of protein-RNA interfaces is critical for understanding the structural bases and functional implications of such interactions and for developing effective approaches to rational drug design. Sequence-based computational methods offer a viable, cost-effective way to identify putative RNA-binding residues in RNA-binding proteins. Here we report two novel approaches: (i) HomPRIP, a sequence homology-based method for predicting RNA-binding sites in proteins; (ii) RNABindRPlus, a new method that combines predictions from HomPRIP with those from an optimized Support Vector Machine (SVM) classifier trained on a benchmark dataset of 198 RNA-binding proteins. Although highly reliable, HomPRIP cannot make predictions for the unaligned parts of query proteins and its coverage is limited by the availability of close sequence homologs of the query protein with experimentally determined RNA-binding sites. RNABindRPlus overcomes these limitations. We compared the performance of HomPRIP and RNABindRPlus with that of several state-of-the-art predictors on two test sets, RB44 and RB111. On a subset of proteins for which homologs with experimentally determined interfaces could be reliably identified, HomPRIP outperformed all other methods achieving an MCC of 0.63 on RB44 and 0.83 on RB111. RNABindRPlus was able to predict RNA-binding residues of all proteins in both test sets, achieving an MCC of 0.55 and 0.37, respectively, and outperforming all other methods, including those that make use of structure-derived features of proteins. More importantly, RNABindRPlus outperforms all other methods for any choice of tradeoff between precision and recall. An important advantage of both HomPRIP and RNABindRPlus is that they rely on readily available sequence and sequence-derived features of RNA-binding proteins. A webserver implementation of both methods is freely available at http://einstein.cs.iastate.edu/RNABindRPlus/.

Present and future of NMR for RNA-protein complexes: a perspective of integrated structural biology

Nucleic acids are gaining enormous importance as key molecules in almost all biological processes. Most nucleic acids do not act in isolation but are generally associated with proteins to form high-molecular-weight nucleoprotein complexes. In this perspective article I focus on the structural studies of supra-molecular ribonucleoprotein (RNP) assemblies in solution by a combination of state-of-the-art TROSY-based NMR experiments and other structural biology techniques. I discuss ways how to combine sparse NMR data with low-resolution structural information from small-angle scattering, fluorescence and electron paramagnetic resonance spectroscopy to obtain the structure of large RNP particles by an integrated structural biology approach. In the last section I give a perspective for the study of RNP complexes by solid-state NMR.

Intrinsic dynamics of an extended hydrophobic core in the S. cerevisiae RNase III dsRBD contributes to recognition of specific RNA binding sites

The Saccharomyces cerevisiae RNase III enzyme Rnt1p preferentially binds to double-stranded RNA hairpin substrates with a conserved (A/u)GNN tetraloop fold, via shape-specific interactions by its double-stranded RNA-binding domain (dsRBD) helix α1 to the tetraloop minor groove. To investigate whether conformational flexibility in the dsRBD regulates the binding specificity, we determined the backbone dynamics of the Rnt1p dsRBD in the free and AGAA hairpin-bound states using NMR spin-relaxation experiments. The intrinsic microsecond-to-millisecond timescale dynamics of the dsRBD suggests that helix α1 undergoes conformational sampling in the free state, with large dynamics at some residues in the α1-β1 loop (α1-β1 hinge). To correlate free dsRBD dynamics with structural changes upon binding, we determined the solution structure of the free dsRBD used in the previously determined RNA-bound structures. The Rnt1p dsRBD has an extended hydrophobic core comprising helix α1, the α1-β1 loop, and helix α3. Analysis of the backbone dynamics and structures of the free and bound dsRBD reveals that slow-timescale dynamics in the α1-β1 hinge are associated with concerted structural changes in the extended hydrophobic core that govern binding of helix α1 to AGAA tetraloops. The dynamic behavior of the dsRBD bound to a longer AGAA hairpin reveals that dynamics within the hydrophobic core differentiate between specific and nonspecific sites. Mutations of residues in the α1-β1 hinge result in changes to the dsRBD stability and RNA-binding affinity and cause defects in small nucleolar RNA processing invivo. These results reveal that dynamics in the extended hydrophobic core are important for binding site selection by the Rnt1p dsRBD.

Protein-RNA footprinting: an evolving tool

As more RNA molecules with important cellular functions are discovered, there is a strong need to characterize their structures, functions, and interactions. Chemical and enzymatic footprinting methods are used to map RNA secondary and tertiary structure, to monitor ligand interactions and conformational changes, and in the study of protein-RNA interactions. These methods provide data at single-nucleotide resolution that nicely complements the structural information available from X-ray diffraction, nuclear magnetic resonance spectroscopy (NMR), or cryo-electron microscopy. Footprinting methods also complement the dynamic information derived from single-molecule Förster resonance energy transfer. RNA footprinting tools have been used for decades, but we have recently seen spectacular advances, for instance, the use in combination with massive parallel sequencing techniques. Large libraries of RNA molecules (small or large in size) can now be probed in high-throughput manner when RNA footprinting methods are combined with fluorescent probe technologies and automation. In this article, after a brief historical overview, we summarize recent advances in RNA-protein footprinting methodologies that now integrate tools for massive parallel analysis.

NMR studies of protein-RNA interactions

This chapter describes the preparation of NMR quantities of RNA purified to single-nucleotide resolution for protein-RNA interaction studies. The protocol is easily modified to make nucleotide-specific isotopically labeled RNAs or uniformly labeled RNA fragments for ligation to generate segmentally labeled RNAs.

Structure of a yeast RNase III dsRBD complex with a noncanonical RNA substrate provides new insights into binding specificity of dsRBDs

dsRBDs often bind dsRNAs with some specificity, yet the basis for this is poorly understood. Rnt1p, the major RNase III in Saccharomyces cerevisiae, cleaves RNA substrates containing hairpins capped by A/uGNN tetraloops, using its dsRBD to recognize a conserved tetraloop fold. However, the identification of a Rnt1p substrate with an AAGU tetraloop raised the question of whether Rnt1p binds to this noncanonical substrate differently than to A/uGNN tetraloops. The solution structure of Rnt1p dsRBD bound to an AAGU-capped hairpin reveals that the tetraloop undergoes a structural rearrangement upon binding to Rnt1p dsRBD to adopt a backbone conformation that is essentially the same as the AGAA tetraloop, and indicates that a conserved recognition mode is used for all Rnt1p substrates. Comparison of free and RNA-bound Rnt1p dsRBD reveals that tetraloop-specific binding requires a conformational change in helix α1. Our findings provide a unified model of binding site selection by this dsRBD.

Preparation of selective and segmentally labeled single-stranded DNA for NMR by self-primed PCR and asymmetrical endonuclease double digestion

We demonstrate a new, efficient and easy-to-use method for enzymatic synthesis of (stereo-)specific and segmental (13)C/(15)N/(2)H isotope-labeled single-stranded DNA in amounts sufficient for NMR, based on the highly efficient self-primed PCR. To achieve this, new approaches are introduced and combined. (i) Asymmetric endonuclease double digestion of tandem-repeated PCR product. (ii) T4 DNA ligase mediated ligation of two ssDNA segments. (iii) In vitro dNTP synthesis, consisting of in vitro rNTP synthesis followed by enzymatic stereo-selective reduction of the C2' of the rNTP, and a one-pot add-up synthesis of dTTP from dUTP. The method is demonstrated on two ssDNAs: (i) a 36-nt three-way junction, selectively (13)C(9)/(15)N(3)/(2)H((1',2'',3',4',5',5''))-dC labeled and (ii) a 39-nt triple-repeat three-way junction, selectively (13)C(9)/(15)N(3)/(2)H((1',2'',3',4',5',5''))-dC and (13)C(9)/(15)N(2)/(2)H((1',2'',3',4',5',5''))-dT labeled in segment C20-C39. Their NMR spectra show the spectral simplification, while the stereo-selective (2)H-labeling in the deoxyribose of the dC-residues, straightforwardly provided assignment of their C1'-H2' and C2'-H2' resonances. The labeling protocols can be extended to larger ssDNA molecules and to more than two segments.

The 3'-flap pocket of human flap endonuclease 1 is critical for substrate binding and catalysis

Flap endonuclease 1 (FEN1) proteins, which are present in all kingdoms of life, catalyze the sequence-independent hydrolysis of the bifurcated nucleic acid intermediates formed during DNA replication and repair. How FEN1s have evolved to preferentially cleave flap structures is of great interest especially in light of studies wherein mice carrying a catalytically deficient FEN1 were predisposed to cancer. Structural studies of FEN1s from phage to human have shown that, although they share similar folds, the FEN1s of higher organisms contain a 3'-extrahelical nucleotide (3'-flap) binding pocket. When presented with 5'-flap substrates having a 3'-flap, archaeal and eukaryotic FEN1s display enhanced reaction rates and cleavage site specificity. To investigate the role of this interaction, a kinetic study of human FEN1 (hFEN1) employing well defined DNA substrates was conducted. The presence of a 3'-flap on substrates reduced Km and increased multiple- and single turnover rates of endonucleolytic hydrolysis at near physiological salt concentrations. Exonucleolytic and fork-gap-endonucleolytic reactions were also stimulated by the presence of a 3'-flap, and the absence of a 3'-flap from a 5'-flap substrate was more detrimental to hFEN1 activity than removal of the 5'-flap or introduction of a hairpin into the 5'-flap structure. hFEN1 reactions were predominantly rate-limited by product release regardless of the presence or absence of a 3'-flap. Furthermore, the identity of the stable enzyme product species was deduced from inhibition studies to be the 5'-phosphorylated product. Together the results indicate that the presence of a 3'-flap is the critical feature for efficient hFEN1 substrate recognition and catalysis.

Multiple segmental and selective isotope labeling of large RNA for NMR structural studies

Multiple segmental and selective isotope labeling of RNA with three segments has been demonstrated by introducing an RNA segment, selectively labeled with ¹³C₉/¹⁵N₂/²H_{(1′, 3′, 4′, 5′, 5′′)}-labeled uridine residues, into the central position of the 20 kDa ε-RNA of Duck Hepatitis B Virus. The RNA molecules were produced via two efficient protocols: a two-step protocol, which uses T4 DNA ligase and T4 RNA ligase 1, and a one-pot protocol, which uses T4 RNA ligase 1 alone. With T4 RNA ligase 1 all not-to-be-ligated termini are usually protected to prevent formation of side products. We show that such labor-intensive protection of termini is not required, provided segmentation sites can be chosen such that the segments fold into the target structure or target-like structures and thus are not trapped into stable alternate structures. These sites can be reliably predicted via DINAMelt. The simplified NMR spectrum provided evidence for the presence of a U28 H₃-imino resonance, previously obscured in the fully labeled sample, and thus of the non-canonical base pair U28:C37. The demonstrated multiple segmental labeling protocols are generally applicable to large RNA molecules and can be extended to more than three segments.

DMS footprinting of structured RNAs and RNA-protein complexes

We describe a protocol in which dimethyl sulfate (DMS) modification of the base-pairing faces of unpaired adenosine and cytidine nucleotides is used for structural analysis of RNAs and RNA-protein complexes (RNPs). The protocol is optimized for RNAs of small to moderate size (< or = 500 nt). The RNA or RNP is first exposed to DMS under conditions that promote formation of the folded structure or complex, as well as 'control' conditions that do not allow folding or complex formation. The positions and extents of modification are then determined by primer extension, polyacrylamide gel electrophoresis and quantitative analysis. From changes in the extent of modification upon folding or protein binding (appearance of a 'footprint'), it is possible to detect local changes in the secondary and tertiary structure of RNA, as well as the formation of RNA-protein contacts. This protocol takes 1.5-3 d to complete, depending on the type of analysis used.

RNA structure determination by NMR

This chapter reviews the methodologies for RNA structure determination by liquid-state nuclear magnetic resonance (NMR). The routine production of milligram quantities of isotopically labeled RNA remains critical to the success of NMR-based structure studies. The standard method for the preparation of isotopically labeled RNA for structural studies in solution is in vitro transcription from DNA oligonucleotide templates using T7 RNA polymerase and unlabeled or isotopically labeled nucleotide triphosphates (NTPs). The purification of the desired RNA can be performed by either denaturing polyacrylamide gel electrophoresis (PAGE) or anion-exchange chromatography. Our basic strategy for studying RNA in solution by NMR is outlined. The topics covered include RNA resonance assignment, restraint collection, and the structure calculation process. Selected examples of NMR spectra are given for a correctly folded 30 nucleotide-containing RNA.

Short, synthetic and selectively 13C-labeled RNA sequences for the NMR structure determination of protein-RNA complexes

We report an optimized synthesis of all canonical 2'-O-TOM protected ribonucleoside phosphoramidites and solid supports containing [13C5]-labeled ribose moieties, their sequence-specific introduction into very short RNA sequences and their use for the structure determination of two protein-RNA complexes. These specifically labeled sequences facilitate RNA resonance assignments and are essential to assign a high number of sugar-sugar and intermolecular NOEs, which ultimately improve the precision and accuracy of the resulting structures. This labeling strategy is particularly useful for the study of protein-RNA complexes with single-stranded RNA in solution, which is rapidly an increasingly relevant research area in biology.

Specific recognition of HIV TAR RNA by the dsRNA binding domains (dsRBD1-dsRBD2) of PKR

PKR (double-stranded RNA-dependent protein kinase) is an important component of host defense to virus infection. Binding of dsRNA to two dsRBDs (double-stranded RNA binding domains) of PKR modulates its own kinase activation. How structural features of natural target RNAs, such as bulges and loops, have an effect on the binding to two dsRBDs of PKR still remains unclear. By using ITC and NMR, we show here that both the bulge and loop of TAR RNA are necessary for the high affinity binding to dsRBD1-dsRBD2 of PKR with 1:1 stoichiometry. The binding site for the dsRBD1-dsRBD2 spans from upper bulge to lower stem of the TAR RNA, based on chemical shift mapping. The backbone resonances in the 40 kDa TAR.dsRBD1-dsRBD2 were assigned. NMR chemical shift perturbation data suggest that the beta1-beta2 loop of the dsRBD1 interacts with the TAR RNA, whereas that of the dsRBD2 is less involved in the TAR RNA recognition. In addition, the residues of the interdomain linker between the dsRBD1 and the dsRBD2 also show large chemical perturbations indicating that the linker is involved in the recognition of TAR RNA. The results presented here provide the biophysical and spectroscopic basis for high-resolution structural studies, and show how local RNA structural features modulate recognition by dsRBDs.