RNAPack: an integrated NMR approach to RNA structure determination

Over the last decade, a vast number of useful nuclear magnetic resonance (NMR) experiments have been developed and successfully employed to determine the structure and dynamics of RNA oligonucleotides. Despite this progress, high-resolution RNA structure determination by NMR spectroscopy still remains a lengthy process and requires programming and extensive calibrations to perform NMR experiments successfully. To accelerate RNA structure determination by NMR spectroscopy, we have designed and programmed a package of RNA NMR experiments, called RNAPack. The user-friendly package contains a set of semiautomated single, double, and triple resonance NMR experiments, which are fully optimized for high-resolution RNA solution structure determination on Varian NMR spectrometers. RNAPack provides an autocalibration feature that allows rapid calibration of all NMR experiments in a single step and thereby speeds up the NMR data collection and eliminates user errors. In our laboratory, we have successfully employed this technology to solve RNA solution structures of domains of the internal ribosome entry site of the genomic hepatitis C viral RNA in less than 3 months. RNAPack therefore makes NMR spectroscopy an attractive and rapid structural tool and allows integration of atomic resolution structural information into biochemical studies of large RNA systems.     

Mechanisms of allosteric gene regulation by NMR quantification of microsecond-millisecond protein dynamics

The trp RNA-binding attenuation protein (TRAP) is a paradigmatic allosteric protein that regulates the tryptophan biosynthetic genes associated with the trp operon in bacilli. The ring-shaped 11-mer TRAP is activated for recognition of a specific trp-mRNA target by binding up to 11 tryptophan molecules. To characterize the mechanisms of tryptophan-induced TRAP activation, we have performed methyl relaxation dispersion (MRD) nuclear magnetic resonance (NMR) experiments that probe the time-dependent structure of TRAP in the microsecond-to-millisecond "chemical exchange" time window. We find significant side chain flexibility localized to the RNA and tryptophan binding sites of the apo protein and that these dynamics are dramatically reduced upon ligand binding. Analysis of the MRD NMR data provides insights into the structural nature of transiently populated conformations sampled in solution by apo TRAP. The MRD data are inconsistent with global two-state exchange, indicating that conformational sampling in apo TRAP is asynchronous. These findings imply a temporally heterogeneous population of structures that are incompatible with RNA binding and substantiate the study of TRAP as a paradigm for probing and understanding essential dynamics in allosteric, regulatory proteins.     

Structure, dynamics and function of the outer membrane protein A (OmpA) and influenza hemagglutinin fusion domain in detergent micelles by solution NMR

Recent progress from our laboratories to determine structures of small membrane proteins (up to 20 kDa) in detergent micelles by solution nuclear magnetic resonance (NMR) is reviewed. NMR opens a new window to also study, for the first time, the dynamics of membrane proteins. We report on recent attempts to correlate dynamic measurements on OmpA with the ion channel function of this protein. We also summarize how NMR and spin-label electron paramagnetic resonance spectroscopy and selective mutagenesis can be combined to provide a structural basis towards understanding the mechanism of influenza hemagglutinin-mediated membrane fusion.     

Assignment of the L30-mRNA complex using selective isotopic labeling and RNA mutants.

The helix-loop-helix structure formed in the pre-mRNA and the mRNA of L30, a ribosomal protein from the yeast Saccharomyces cerevisiae, serves as an auto-regulatory binding site for the protein to suppress the L30 synthesis upon overproduction. Using a 33-nucleotide model RNA, the structures of the L30 binding site RNA in the presence and absence of the protein were investigated using nuclear magnetic resonance (NMR) spectroscopy. Homonuclear and(13)C/(15)N-based resonance assignments and spectral comparisons indicated that the purine-rich internal loop is dynamic in the free RNA but becomes ordered in the presence of L30 protein. Although the resonances in the loop region are sharper and more disperse in the bound RNA, their assignment was extremely challenging, due to spectral complexity and broadened resonances caused by local dynamics. Two strategies, namely selective(13)C/(15)N-labeling and NMR analyses of five complexes with RNA mutants, were used to overcome these difficulties. Only using these approaches could assignment of the internal loop resonances and identification of the unusual NOEs and nucleotide conformations within the internal loop be made. In the case of structural determination of the L30-mRNA complex, it was critical to be able to take advantage of the available biochemical information in order to complete the structure determination.     

Recent developments in structural proteomics for protein structure determination

The major challenges in structural proteomics include identifying all the proteins on the genome-wide scale, determining their structure-function relationships, and outlining the precise three-dimensional structures of the proteins. Protein structures are typically determined by experimental approaches such as X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy. However, the knowledge of three-dimensional space by these techniques is still limited. Thus, computational methods such as comparative and de novo approaches and molecular dynamic simulations are intensively used as alternative tools to predict the three-dimensional structures and dynamic behavior of proteins. This review summarizes recent developments in structural proteomics for protein structure determination; including instrumental methods such as X-ray crystallography and NMR spectroscopy, and computational methods such as comparative and de novo structure prediction and molecular dynamics simulations.     

Conjoined Use of EM and NMR in RNA Structure Refinement

More than 40% of the RNA structures have been determined using nuclear magnetic resonance (NMR) technique. NMR mainly provides local structural information of protons and works most effectively on relatively small biomacromolecules. Hence structural characterization of large RNAs can be difficult for NMR alone. Electron microscopy (EM) provides global shape information of macromolecules at nanometer resolution, which should be complementary to NMR for RNA structure determination. Here we developed a new energy term in Xplor-NIH against the density map obtained by EM. We conjointly used NMR and map restraints for the structure refinement of three RNA systems — U2/U6 small-nuclear RNA, genome-packing motif (Ψ^CD)₂ from Moloney murine leukemia virus, and ribosome-binding element from turnip crinkle virus. In all three systems, we showed that the incorporation of a map restraint, either experimental or generated from known PDB structure, greatly improves structural precision and accuracy. Importantly, our method does not rely on an initial model assembled from RNA duplexes, and allows full torsional freedom for each nucleotide in the torsion angle simulated annealing refinement. As increasing number of macromolecules can be characterized by both NMR and EM, the marriage between the two techniques would enable better characterization of RNA three-dimensional structures.     

The model structure of the hammerhead ribozyme formed by RNAs of reciprocal chirality

RNA-based tools are frequently used to modulate gene expression in living cells. However, the stability and effectiveness of such RNA-based tools is limited by cellular nuclease activity. One way to increase RNA’s resistance to nucleases is to replace its D-ribose backbone with L-ribose isomers. This modification changes chirality of an entire RNA molecule to L-form giving it more chance of survival when introduced into cells. Recently, we have described the activity of left-handed hammerhead ribozyme (L-Rz, L-HH) that can specifically hydrolyse RNA with the opposite chirality at a predetermined location. To understand the structural background of the RNA specific cleavage in a heterochiral complex, we used circular dichroism (CD) and nuclear magnetic resonance (NMR) spectroscopy as well as performed molecular modelling and dynamics simulations of homo- and heterochiral RNA complexes. The active ribozyme-target heterochiral complex showed a mixed chirality as well as low field imino proton NMR signals. We modelled the 3D structures of the oligoribonucleotides with their ribozyme counterparts of reciprocal chirality. L- or D-ribozyme formed a stable, homochiral helix 2, and two short double heterochiral helixes 1 and 3 of D- or L-RNA strand thorough irregular Watson–Crick base pairs. The formation of the heterochiral complexes is supported by the result of simulation molecular dynamics. These new observations suggest that L-catalytic nucleic acids can be used as tools in translational biology and diagnostics.     

Assessing and refining molecular dynamics simulations of proteins with nuclear magnetic resonance data

The sophistication of the force fields, algorithms and hardware used for molecular dynamics (MD) simulations of proteins is continuously increasing. No matter how advanced the methodology, however, it is essential to evaluate the appropriateness of the structures sampled in a simulation by comparison with quantitative experimental data. Solution nuclear magnetic resonance (NMR) data are particularly useful for checking the quality of protein simulations, as they provide both structural and dynamic information on a variety of temporal and spatial scales. Here, various features and implications of using NMR data to validate and bias MD simulations are outlined, including an overview of the different types of NMR data that report directly on structural properties and of relevant simulation techniques. The focus throughout is on how to properly account for conformational averaging, particularly within the context of the assumptions inherent in the relationships that link NMR data to structural properties.     

Unraveling the structural complexity in a single-stranded RNA tail: implications for efficient ligand binding in the prequeuosine riboswitch

Single-stranded RNAs (ssRNAs) are ubiquitous RNA elements that serve diverse functional roles. Much of our understanding of ssRNA conformational behavior is limited to structures in which ssRNA directly engages in tertiary interactions or is recognized by proteins. Little is known about the structural and dynamic behavior of free ssRNAs at atomic resolution. Here, we report the collaborative application of nuclear magnetic resonance (NMR) and replica exchange molecular dynamics (REMD) simulations to characterize the 12 nt ssRNA tail derived from the prequeuosine riboswitch. NMR carbon spin relaxation data and residual dipolar coupling measurements reveal a flexible yet stacked core adopting an A-form-like conformation, with the level of order decreasing toward the terminal ends. An A-to-C mutation within the polyadenine tract alters the observed dynamics consistent with the introduction of a dynamic kink. Pre-ordering of the tail may increase the efficacy of ligand binding above that achieved by a random-coil ssRNA. The REMD simulations recapitulate important trends in the NMR data, but suggest more internal motions than inferred from the NMR analysis. Our study unmasks a previously unappreciated level of complexity in ssRNA, which we believe will also serve as an excellent model system for testing and developing computational force fields.     

核磁共振波谱应用于结构生物学的研究进展

综述了核磁共振波谱在结构生物学研究中的进展。在溶液中测定生物大分子的结构,分子大小的限制正被减少,尽管新结构的测定仍然需要付出比较大的努力。核磁共振是一个有效的手段,可用于研究在许多细胞过程中存在的弱的或者瞬态的蛋白质-蛋白质相互作用。结构的柔性在蛋白质分子功能中起了中心作用。由于最近方法学的发展,使NMR可以表征蛋白质的动力学,从而可以对分子机制有新的认识。核磁共振波谱可以在原子分辨率下表征无序的蛋白质系统,可以研究折叠路径。跨膜蛋白在细胞中起了关键作用,这使它们成为药物的靶标。应用液体和固体核磁共振技术已经成功测定了跨膜蛋白质的结构。     

Local folding coupled to RNA binding in the yeast ribosomal protein L30.

The ribosomal protein L30 from yeast Saccharomyces cerevisiae auto-regulates its own synthesis by binding to a structural element in both its pre-mRNA and its mRNA. The three-dimensional structures of L30 in the free (f L30) and the pre-mRNA bound (b L30) forms have been solved by nuclear magnetic resonance spectroscopy. Both protein structures contain four alternating alpha-helices and four beta-strands segments and adopt an overall topology that is an alphabetaalpha three-layer sandwich, representing a unique fold. Three loops on one end of the alphabetaalpha sandwich have been mapped as the RNA binding site on the basis of structural comparison, chemical shift perturbation and the inter-molecular nuclear Overhauser effects to the RNA. The structural and dynamic comparison of f L30 and b L30 reveals that local dynamics may play an important role in the RNA binding. The fourth helix in b L30 is longer than in f L30, and is stabilized by RNA binding. The exposed hydrophobic surface that is buried upon RNA binding may provide the energy necessary to drive secondary structure formation, and may account for the increased stability of b L30.     

Three-dimensional structure of bradykinin in SDS micelles. Study using nuclear magnetic resonance, distance geometry, and restrained molecular mechanics and dynamics

The conformational properties of bradykinin in five molar excess sodium dodecyl sulfate (SDS) micelles have been examined by two-dimensional nuclear magnetic resonance (NMR) techniques at 500 MHz. Detailed structural information for bradykinin in SDS was obtained from quantitative 2-D nuclear Overhauser enhancement (n.O.e.) analyses, distance geometry, and restrained molecular mechanics and dynamics calculations. The conformation of bradykinin in SDS micelles, as determined by these methods, is characterized by a beta-turn-like structure at residues 6-9. A detailed comparison of the structures derived from distance geometry and restrained molecular mechanics and dynamics is also presented.     

Structural model of the HIV-1 Tat(46-58)-TAR complex

The trans-activator protein (Tat) of human immunodeficiency virus type 1 (HIV-1) binds to an uridine-rich bulge of an RNA target (TAR; trans-activation responsive element) predominantly via its basic sequence domain. The structure of the Tat(46-58)-TAR complex has been determined by a novel modeling approach relying on structural information about one crucial arginine residue and crosslink data. The strategy described here solely uses this experimental data without additional "modeling" assumptions about the structure of the complex in order to avoid human bias. Model building was performed in a fashion similar to structure calculations from nuclear magnetic resonance (NMR)-spectroscopic data using restrained molecular dynamics. The resulting set of structures of Tat(46-58) in its complex with TAR reveals that all models have converged to a common fold, showing a backbone root mean square deviation (RMSD) of 1.36A. Analysis of the calculated structures suggests that HIV-I Tat forms a hairpin loop in its complex with TAR that shares striking similarity to the hairpin formed by the structure of the bovine immunodeficiency virus Tat protein after TAR binding as determined by NMR studies. The outlined approach is not limited to the Tat-TAR complex modeling, but is also applicable to all molecular complexes with sufficient biochemical and biophysical data available.     

High-resolution NMR structure of an RNA model system: the 14-mer cUUCGg tetraloop hairpin RNA

We present a high-resolution nuclear magnetic resonance (NMR) solution structure of a 14-mer RNA hairpin capped by cUUCGg tetraloop. This short and very stable RNA presents an important model system for the study of RNA structure and dynamics using NMR spectroscopy, molecular dynamics (MD) simulations and RNA force-field development. The extraordinary high precision of the structure (root mean square deviation of 0.3 Å) could be achieved by measuring and incorporating all currently accessible NMR parameters, including distances derived from nuclear Overhauser effect (NOE) intensities, torsion-angle dependent homonuclear and heteronuclear scalar coupling constants, projection-angle-dependent cross-correlated relaxation rates and residual dipolar couplings. The structure calculations were performed with the program CNS using the ARIA setup and protocols. The structure quality was further improved by a final refinement in explicit water using OPLS force field parameters for non-bonded interactions and charges. In addition, the 2′-hydroxyl groups have been assigned and their conformation has been analyzed based on NOE contacts. The structure currently defines a benchmark for the precision and accuracy amenable to RNA structure determination by NMR spectroscopy. Here, we discuss the impact of various NMR restraints on structure quality and discuss in detail the dynamics of this system as previously determined.     

Protein dynamics from solution NMR

Solution nuclear magnetic resonance (NMR) spectroscopy is unique in its ability to elucidate the details of atomic-level structural and dynamical properties of biological macromolecules under native-like conditions. Recent advances in NMR techniques and protein sample preparation now allow comprehensive investigation of protein dynamics over timescales ranging 14 orders of magnitude at nearly every atomic site. Thus, solution NMR is poised to reveal aspects of the physico-chemical properties that govern the ensemble distribution of protein conformers and the dynamics of their interconversion. We review these advances as well as their recent application to the study of proteins.     

Protein dynamics from solution NMR: theory and applications

Solution nuclear magnetic resonance (NMR) spectroscopy is unique in its ability to elucidate the details of atomic-level structural and dynamical properties of biological macromolecules under native-like conditions. Recent advances in NMR techniques and protein sample preparation now allow comprehensive investigation of protein dynamics over timescales ranging 14 orders of magnitude at nearly every atomic site. Thus, solution NMR is poised to reveal aspects of the physico-chemical properties that govern the ensemble distribution of protein conformers and the dynamics of their interconversion. We review these advances as well as their recent application to the study of proteins.     

NMR of membrane proteins in micelles and bilayers: the FXYD family proteins

Determining the atomic resolution structures of membrane proteins is of particular interest in contemporary structural biology. Helical membrane proteins constitute one-third of the expressed proteins encoded in a genome, many drugs have membrane-bound proteins as their receptors, and mutations in membrane proteins result in human diseases. Although integral membrane proteins provide daunting technical challenges for all methods of protein structure determination, nuclear magnetic resonance (NMR) spectroscopy can be an extremely versatile and powerful method for determining their structures and characterizing their dynamics, in lipid environments that closely mimic the cell membranes. Once milligram amounts of isotopically labeled protein are expressed and purified, micelle samples can be prepared for solution NMR analysis, and lipid bilayer samples can be prepared for solid-state NMR analysis. The two approaches are complementary and can provide detailed structural and dynamic information. This paper describes the steps for membrane protein structure determination using solution and solid-state NMR. The methods for protein expression and purification, sample preparation and NMR experiments are described and illustrated with examples from the FXYD proteins, a family of regulatory subunits of the Na,K-ATPase.     

Integrating experimental data with molecular simulations to investigate RNA structural dynamics

Conformational dynamics is crucial for ribonucleic acid (RNA) function. Techniques such as nuclear magnetic resonance, cryo-electron microscopy, small- and wide-angle X-ray scattering, chemical probing, single-molecule Förster resonance energy transfer, or even thermal or mechanical denaturation experiments probe RNA dynamics at different time and space resolutions. Their combination with accurate atomistic molecular dynamics (MD) simulations paves the way for quantitative and detailed studies of RNA dynamics. First, experiments provide a quantitative validation tool for MD simulations. Second, available data can be used to refine simulated structural ensembles to match experiments. Finally, comparison with experiments allows for improving MD force fields that are transferable to new systems for which data is not available. Here we review the recent literature and provide our perspective on this field.