Towards structural genomics of RNA: rapid NMR resonance assignment and simultaneous RNA tertiary structure determination using residual dipolar couplings

We report a new residual dipolar couplings (RDCs) based NMR procedure for rapidly determining RNA tertiary structure demonstrated on a uniformly (15)N/(13)C-labeled 27 nt variant of the trans-activation response element (TAR) RNA from HIV-I. In this procedure, the time-consuming nuclear Overhauser enhancement (NOE)-based sequential assignment step is replaced by a fully automated RDC-based assignment strategy. This approach involves examination of all allowed sequence-specific resonance assignment permutations for best-fit agreement between measured RDCs and coordinates for sub-structures in a target RNA. Using idealized A-form geometries to model Watson-Crick helices and coordinates from a previous X-ray structure to model a hairpin loop in TAR, the best-fit RDC assignment solutions are determined very rapidly (相似文献   

Impact of static and dynamic A-form heterogeneity on the determination of RNA global structural dynamics using NMR residual dipolar couplings

We examined how static and dynamic deviations from the idealized A-form helix propagate into errors in the principal order tensor parameters determined using residual dipolar couplings (rdcs). A 20-ns molecular dynamics (MD) simulation of the HIV-1 transactivation response element (TAR) RNA together with a survey of spin relaxation studies of RNA dynamics reveals that pico-to-nanosecond local motions in non-terminal Watson-Crick base-pairs will uniformly attenuate base and sugar one bond rdcs by approximately 7%. Gaussian distributions were generated for base and sugar torsion angles through statistical comparison of 40 RNA X-ray structures solved to <3.0 A resolution. For a typical number (>or=11) of one bond C-H base and sugar rdcs, these structural deviations together with rdc uncertainty (1.5 Hz) lead to average errors in the magnitude and orientation of the principal axis of order that are <9% and <4 degrees, respectively. The errors decrease to <5% and <4 degrees for >or=17 rdcs. A protocol that allows for estimation of error in A-form order tensors due to both angular deviations and rdc uncertainty (Aform-RDC) is validated using theoretical simulations and used to analyze rdcs measured previously in TAR in the free state and bound to four distinct ligands. Results confirm earlier findings that the two TAR helices undergo large changes in both their mean relative orientation and dynamics upon binding to different targets.     

3D maps of RNA interhelical junctions

More than 50% of RNA secondary structure is estimated to be A-form helices, which are linked together by various junctions. Here we describe a protocol for computing three interhelical Euler angles describing the relative orientation of helices across RNA junctions. 5' and 3' helices, H1 and H2, respectively, are assigned based on the junction topology. A reference canonical helix is constructed using an appropriate molecular builder software consisting of two continuous idealized A-form helices (iH1 and iH2) with helix axis oriented along the molecular Z-direction running toward the positive direction from iH1 to iH2. The phosphate groups and the carbon and oxygen atoms of the sugars are used to superimpose helix H1 of a target interhelical junction onto the corresponding iH1 of the reference helix. A copy of iH2 is then superimposed onto the resulting H2 helix to generate iH2'. A rotation matrix R is computed, which rotates iH2' into iH2 and expresses the rotation parameters in terms of three Euler angles α(h), β(h) and γ(h). The angles are processed to resolve a twofold degeneracy and to select an overall rotation around the axis of the reference helix. The three interhelical Euler angles define clockwise rotations around the 5' (-γ(h)) and 3' (α(h)) helices and an interhelical bend angle (β(h)). The angles can be depicted graphically to provide a 'Ramachandran'-type view of RNA global structure that can be used to identify unusual conformations as well as to understand variations due to changes in sequence, junction topology and other parameters.     

Effect of single-residue bulges on RNA double-helical structures: crystallographic database analysis and molecular dynamics simulation studies

Asymmetric bulge loop motifs are widely dispersed in all types of functional RNAs. They are frequently occurring structural motifs in folded RNA structures and appear commonly in pre-microRNA and ribosomes, where they are involved in specific RNA–RNA and RNA–protein interactions. It is therefore necessary to understand such motifs from a structural point of view. We analyzed all available RNA structures and identified quite a few fragments of double helices that contain bulges. We found that these discontinuities often introduce kinks into the double helices, which also affects the stacking overlap between the base pairs across the irregularity. In order to understand the influence of these bulges on stability and flexibility, we carried out molecular dynamics simulations of three different single-residue bulge-containing RNA helices using the CHARMM36 force field. The structural variability at the junctions of RNA bulges is expected to differ from that in continuous double-helical stretches. The structural features of the junction region were observed to vary noticeably depending on the orientation of the bulge residue. When the base of the bulge residue is looped out, the RNA stretch behaves like a standard long A-form RNA double helix, whereas the entire RNA behaves differently when the base of the bulge residue is intercalated between base pairs inside the RNA stem. Such single-base intercalation was found to introduce a permanent kink into the composite double helix, which could be a recognition element for Dicer during the maturation of miRNA.     

Coplanar and coaxial orientations of RNA bases and helices

Electrostatic interactions, base-pairing, and especially base-stacking dominate RNA three-dimensional structures. In an A-form RNA helix, base-stacking results in nearly perfect parallel orientations of all bases in the helix. Interestingly, when an RNA structure containing multiple helices is visualized at the atomic level, it is often possible to find an orientation such that only the edges of most bases are visible. This suggests that a general aspect of higher level RNA structure is a coplanar arrangement of base-normal vectors. We have analyzed all solved RNA crystal structures to determine the degree to which RNA base-normal vectors are globally coplanar. Using a statistical test based on the Watson-Girdle distribution, we determined that 330 out of 331 known RNA structures show statistically significant (p < 0.05; false discovery rate [FDR] = 0.05) coplanar normal vector orientations. Not surprisingly, 94% of the helices in RNA show bipolar arrangements of their base-normal vectors (p < 0.05). This allows us to compute a mean axis for each helix and compare their orientations within an RNA structure. This analysis revealed that 62% (208/331) of the RNA structures exhibit statistically significant coaxial packing of helices (p < 0.05, FDR = 0.08). Further analysis reveals that the bases in hairpin loops and junctions are also generally planar. This work demonstrates coplanar base orientation and coaxial helix packing as an emergent behavior of RNA structure and may be useful as a structural modeling constraint.     

Comparison of alignment tensors generated for native tRNAVal using magnetic fields and liquid crystalline media

Residual dipolar couplings (RDCs) complement standard NOE distance and J-coupling torsion angle data to improve the local and global structure of biomolecules in solution. One powerful application of RDCs is for domain orientation studies, which are especially valuable for structural studies of nucleic acids, where the local structure of a double helix is readily modeled and the orientations of the helical domains can then be determined from RDC data. However, RDCs obtained from only one alignment media generally result in degenerate solutions for the orientation of multiple domains. In protein systems, different alignment media are typically used to eliminate this orientational degeneracy, where the combination of RDCs from two (or more) independent alignment tensors can be used to overcome this degeneracy. It is demonstrated here for native E. coli tRNA^Val that many of the commonly used liquid crystalline alignment media result in very similar alignment tensors, which do not eliminate the 4-fold degeneracy for orienting the two helical domains in tRNA. The intrinsic magnetic susceptibility anisotropy (MSA) of the nucleobases in tRNA^Val was also used to obtain RDCs for magnetic alignment at 800 and 900 MHz. While these RDCs yield a different alignment tensor, the specific orientation of this tensor combined with the high rhombicity for the tensors in the liquid crystalline media only eliminates two of the four degenerate orientations for tRNA^Val. Simulations are used to show that, in optimal cases, the combination of RDCs obtained from liquid crystalline medium and MSA-induced alignment can be used to obtain a unique orientation for the two helical domains in tRNA^Val. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Measurement of 1H�C15N and 1H�C13C residual dipolar couplings in nucleic acids from TROSY intensities

Analogous to the recently introduced ARTSY method for measurement of one-bond (1)H-(15)N residual dipolar couplings (RDCs) in large perdeuterated proteins, we introduce methods for measurement of base (13)C-(1)H and (15)N-(1)H RDCs in protonated nucleic acids. Measurements are based on quantitative analysis of intensities in (1)H-(15)N and (13)C-(1)H TROSY-HSQC spectra, and are illustrated for a 71-nucleotide adenine riboswitch. Results compare favorably with those of conventional frequency-based measurements in terms of completeness and convenience of use. The ARTSY method derives the size of the coupling from the ratio of intensities observed in two TROSY-HSQC spectra recorded with different dephasing delays, thereby minimizing potential resonance overlap problems. Precision of the RDC measurements is limited by the signal-to-noise ratio, S/N, achievable in the 2D TROSY-HSQC reference spectrum, and is approximately given by 30/(S/N) Hz for (15)N-(1)H and 65/(S/N) Hz for (13)C-(1)H. The signal-to-noise ratio of both (1)H-(15)N and (1)H-(13)C spectra greatly benefits when water magnetization during the experiments is not perturbed, such that rapid magnetization transfer from bulk water to the nucleic acid, mediated by rapid amino and hydroxyl hydrogen exchange coupled with (1)H-(1)H NOE transfer, allows for fast repetition of the experiment. RDCs in the mutated helix 1 of the riboswitch are compatible with nucleotide-specifically modeled, idealized A-form geometry and a static orientation relative to the helix 2/3 pair, which differs by ca 6° relative to the X-ray structure of the native riboswitch.     

Spectroscopic studies of chimeric DNA-RNA and RNA 29-base intramolecular triple helices.

Fourier transform infrared (FTIR), UV absorption and exchangeable proton NMR spectroscopies have been used to study the formation and stability of two intramolecular pH-dependent triple helices composed by a chimeric 29mer DNA-RNA (DNA double strand and RNA third strand) or by the analogous 29mer RNA. In both cases decrease of pH induces formation of a triple helical structure containing either rU*dA.dT and rC+*dG.dC or rU*rA.rU and rC+*rG.rC triplets. FTIR spectroscopy shows that exclusively N-type sugars are present in the triple helix formed by the 29mer RNA while both N- and S-type sugars are detected in the case of the chimeric 29mer DNA-RNA triple helix. Triple helix formation with the third strand RNA and the duplex as DNA appears to be associated with the conversion of the duplex part from a B-form secondary structure to one which contains partly A-form sugars. Thermal denaturation experiments followed by UV spectroscopy show that a major stabilization occurs upon formation of the triple helices. Monophasic melting curves indicate a simultaneous disruption of the Hoogsteen and Watson-Crick hydrogen bonds in the intramolecular triplexes when the temperature is increased. This is in agreement with imino proton NMR spectra recorded as a function of temperature. Comparison with experiments concerning intermolecular triplexes of identical base and sugar composition shows the important role played by the two tetrameric loops in the stabilization of the intramolecular triple helices studied.     

Characterisation of the conformational properties of urea-unfolded Im7: implications for the early stages of protein folding

The colicin immunity protein Im7 folds from its unfolded state in 6 M urea to its native four-helix structure through an on-pathway intermediate that lacks one of the helices of the native structure (helix III). In order to further characterize the folding mechanism of Im7, we have studied the conformational properties of the protein unfolded in 6 M urea in detail using heteronuclear NMR. Triple-resonance experiments with 13C/15N-labelled Im7 in 6 M urea provided almost complete resonance assignments for the backbone nuclei, and measurement of backbone 15N relaxation parameters allowed dynamic ordering of the unfolded polypeptide chain to be investigated. Reduced spectral density mapping and fitting backbone R2 relaxation rates to a polymer dynamics model identified four clusters of interacting residues, each predicted by the average area buried upon folding for each residue. Chemical shift analyses and measurement of NOEs detected with a long mixing-time 1H-1H-15N NOESY-HSQC spectrum confirmed the formation of four clusters. Each cluster of interacting side-chains in urea-unfolded Im7 occurs in a region of the protein that forms a helix in the protein, with the largest clusters being associated with the three long helices that are formed in the on-pathway folding intermediate, whilst the smallest cluster forms a helix only in the native state. NMR studies of a Phe15Ala Im7 variant and a protein in which residues 51-56 are replaced by three glycine residues (H3G3 Im7*), indicated that the clusters do not interact with each other, possibly because they are solvated by urea, as indicated by analysis of NOEs between the protein and the solvent. Based on these data, we suggest that dilution of the chaotrope to initiate refolding will result in collapse of the clusters, leading to the formation of persistent helical structure and the generation of the three-helix folding intermediate.     

Both helix topology and counterion distribution contribute to the more effective charge screening in dsRNA compared with dsDNA

The recent discovery of the RNA interference mechanism emphasizes the biological importance of short, isolated, double-stranded (ds) RNA helices and calls for a complete understanding of the biophysical properties of dsRNA. However, most previous studies of the electrostatics of nucleic acid duplexes have focused on DNA. Here, we present a comparative investigation of electrostatic effects in RNA and DNA. Using resonant (anomalous) and non-resonant small-angle X-ray scattering, we characterized the charge screening efficiency and counterion distribution around short (25 bp) dsDNA and RNA molecules of comparable sequence. Consistent with theoretical predictions, we find counterion mediated screening to be more efficient for dsRNA than dsDNA. Furthermore, the topology of the RNA A-form helix alters the spatial distribution of counterions relative to B-form DNA. The experimental results reported here agree well with ion-size-corrected non-linear Poisson–Boltzmann calculations. We propose that differences in electrostatic properties aid in selective recognition of different types of short nucleic acid helices by target binding partners.     

Solution studies of the dimerization initiation site of HIV-1 genomic RNA.

Dimerization of HIV-1 genomic RNA is an essential step of the viral cycle, initiated at a conserved stem-loop structure which forms a 'kissing complex' involving loop-loop interactions (dimerization initiation site, DIS). A 19mer RNA oligonucleotide (DIS-19) has been synthesized which forms a stable symmetrical dimer in solution at millimolar concentrations. Dimerization does not depend on addition of Mg2+. RNA ligation experiments unambiguously indicate that the formed dimer is a stable kissing complex under the NMR experimental conditions.1H NMR resonance assignments were obtained for DIS-19 at 290 K and pH 6.5. Analysis of the pattern of NOE connectivities reveals that the helix formed by loop-loop base pairing is stacked onto the two terminal stems. Non-canonical base pairs between two essential and conserved adenines are found at the junctions between the two intramolecular and the single intramolecular helices.     

Helix conformations in 7TM membrane proteins determined using oriented-sample solid-state NMR with multiple residue-specific 15N labeling

Oriented solid-state NMR in combination with multiple-residue-specific ¹⁵N labeling and extensive numerical spectral analysis is proposed to determine helix conformations of large membrane proteins in native membranes. The method is demonstrated on uniaxially oriented samples of ¹⁵N-methionine, -valine, and -glycine-labeled bacteriorhopsin in native purple membranes. Experimental two-dimensional ¹H-¹⁵N dipole-dipole coupling versus ¹⁵N chemical shift spectra for all samples are analyzed numerically to establish combined constraints on the orientation of the seven transmembrane helices relative to the membrane bilayer normal. Since the method does not depend on specific resonance assignments and proves robust toward nonidealities in the sample alignment, it may be generally feasible for the study of conformational arrangement and function-induced conformation changes of large integral membrane proteins.     

Major groove width variations in RNA structures determined by NMR and impact of <Superscript>13</Superscript>C residual chemical shift anisotropy and <Superscript>1</Superscript>H–<Superscript>13</Superscript>C residual dipolar coupling on refinement

Ribonucleic acid structure determination by NMR spectroscopy relies primarily on local structural restraints provided by ¹H– ¹H NOEs and J-couplings. When employed loosely, these restraints are broadly compatible with A- and B-like helical geometries and give rise to calculated structures that are highly sensitive to the force fields employed during refinement. A survey of recently reported NMR structures reveals significant variations in helical parameters, particularly the major groove width. Although helical parameters observed in high-resolution X-ray crystal structures of isolated A-form RNA helices are sensitive to crystal packing effects, variations among the published X-ray structures are significantly smaller than those observed in NMR structures. Here we show that restraints derived from aromatic ¹H– ¹³C residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs) can overcome NMR restraint and force field deficiencies and afford structures with helical properties similar to those observed in high-resolution X-ray structures.     

Top-down approach in protein RDC data analysis: <Emphasis Type="Italic">de novo</Emphasis> estimation of the alignment tensor

In solution NMR spectroscopy the residual dipolar coupling (RDC) is invaluable in improving both the precision and accuracy of NMR structures during their structural refinement. The RDC also provides a potential to determine protein structure de novo. These procedures are only effective when an accurate estimate of the alignment tensor has already been made. Here we present a top–down approach, starting from the secondary structure elements and finishing at the residue level, for RDC data analysis in order to obtain a better estimate of the alignment tensor. Using only the RDCs from N–H bonds of residues in α-helices and CA–CO bonds in β-strands, we are able to determine the offset and the approximate amplitude of the RDC modulation-curve for each secondary structure element, which are subsequently used as targets for global minimization. The alignment order parameters and the orientation of the major principal axis of individual helix or strand, with respect to the alignment frame, can be determined in each of the eight quadrants of a sphere. The following minimization against RDC of all residues within the helix or strand segment can be carried out with fixed alignment order parameters to improve the accuracy of the orientation. For a helical protein Bax, the three components A _xx , A _yy and A _zz , of the alignment order can be determined with this method in average to within 2.3% deviation from the values calculated with the available atomic coordinates. Similarly for β-sheet protein Ubiquitin they agree in average to within 8.5%. The larger discrepancy in β-strand parameters comes from both the diversity of the β-sheet structure and the lower precision of CA–CO RDCs. This top-down approach is a robust method for alignment tensor estimation and also holds a promise for providing a protein topological fold using limited sets of RDCs.     

Validation of helical tilt angles in the solution NMR structure of the Z domain of Staphylococcal protein A by combined analysis of residual dipolar coupling and NOE data

Staphylococcal protein A (SpA) is a virulence factor from Staphylococcus aureus that is able to bind to immunoglobulins. The 3D structures of its immunoglobulin (Ig) binding domains have been extensively studied by NMR and X-ray crystallography, and are often used as model structures in developing de novo or ab initio strategies for predicting protein structure. These small three-helix-bundle structures, reported in free proteins or Ig-bound complexes, have been determined previously using medium- to high-resolution data. Although the location and relative orientation of the three helices in most of these published 3D domain structures are consistent, there are significant differences among the reported structures regarding the tilt angle of the first helix (helix 1). We have applied residual dipolar coupling data, together with nuclear Overhauser enhancement and scalar coupling data, in refining the NMR solution structure of an engineered IgG-binding domain (Z domain) of SpA. Our results demonstrate that the three helices are almost perfectly antiparallel in orientation, with the first helix tilting slightly away from the other two helices. We propose that this high-accuracy structure of the Z domain of SpA is a more suitable target for theoretical predictions of the free domain structure than previously published lower-accuracy structures of protein A domains.     

A polynomial-time nuclear vector replacement algorithm for automated NMR resonance assignments.

High-throughput NMR structural biology can play an important role in structural genomics. We report an automated procedure for high-throughput NMR resonance assignment for a protein of known structure, or of a homologous structure. These assignments are a prerequisite for probing protein-protein interactions, protein-ligand binding, and dynamics by NMR. Assignments are also the starting point for structure determination and refinement. A new algorithm, called Nuclear Vector Replacement (NVR) is introduced to compute assignments that optimally correlate experimentally measured NH residual dipolar couplings (RDCs) to a given a priori whole-protein 3D structural model. The algorithm requires only uniform( 15)N-labeling of the protein and processes unassigned H(N)-(15)N HSQC spectra, H(N)-(15)N RDCs, and sparse H(N)-H(N) NOE's (d(NN)s), all of which can be acquired in a fraction of the time needed to record the traditional suite of experiments used to perform resonance assignments. NVR runs in minutes and efficiently assigns the (H(N),(15)N) backbone resonances as well as the d(NN)s of the 3D (15)N-NOESY spectrum, in O(n(3)) time. The algorithm is demonstrated on NMR data from a 76-residue protein, human ubiquitin, matched to four structures, including one mutant (homolog), determined either by x-ray crystallography or by different NMR experiments (without RDCs). NVR achieves an assignment accuracy of 92-100%. We further demonstrate the feasibility of our algorithm for different and larger proteins, using NMR data for hen lysozyme (129 residues, 97-100% accuracy) and streptococcal protein G (56 residues, 100% accuracy), matched to a variety of 3D structural models. Finally, we extend NVR to a second application, 3D structural homology detection, and demonstrate that NVR is able to identify structural homologies between proteins with remote amino acid sequences using a database of structural models.     

Structure and dynamics of a membrane protein in micelles from three solution NMR experiments

Three solution NMR experiments on a uniformly ¹⁵N labeled membrane protein in micelles provide sufficient information to describe the structure, topology, and dynamics of its helices, as well as additional information that characterizes the principal features of residues in terminal and inter-helical loop regions. The backbone amide resonances are assigned with an HMQC-NOESY experiment and the backbone dynamics are characterized by a ¹H-¹⁵N heteronuclear NOE experiment, which clearly distinguishes between the structured helical residues and the more mobile residues in the terminal and interhelical loop regions of the protein. The structure and topology of the helices are described by Dipolar waves and PISA wheels derived from experimental measurements of residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs). The results show that the membrane-bound form of Pf1 coat protein has a 20-residue trans-membrane hydrophobic helix with an orientation that differs by about 90° from that of an 8-residue amphipathic helix. This combination of three-experiments that yields Dipolar waves and PISA wheels has the potential to contribute to high-throughput structural characterizations of membrane proteins.     

Protein side-chain resonance assignment and NOE assignment using RDC-defined backbones without TOCSY data

One bottleneck in NMR structure determination lies in the laborious and time-consuming process of side-chain resonance and NOE assignments. Compared to the well-studied backbone resonance assignment problem, automated side-chain resonance and NOE assignments are relatively less explored. Most NOE assignment algorithms require nearly complete side-chain resonance assignments from a series of through-bond experiments such as HCCH-TOCSY or HCCCONH. Unfortunately, these TOCSY experiments perform poorly on large proteins. To overcome this deficiency, we present a novel algorithm, called Nasca (NOE Assignment and Side-Chain Assignment), to automate both side-chain resonance and NOE assignments and to perform high-resolution protein structure determination in the absence of any explicit through-bond experiment to facilitate side-chain resonance assignment, such as HCCH-TOCSY. After casting the assignment problem into a Markov Random Field (MRF), Nasca extends and applies combinatorial protein design algorithms to compute optimal assignments that best interpret the NMR data. The MRF captures the contact map information of the protein derived from NOESY spectra, exploits the backbone structural information determined by RDCs, and considers all possible side-chain rotamers. The complexity of the combinatorial search is reduced by using a dead-end elimination (DEE) algorithm, which prunes side-chain resonance assignments that are provably not part of the optimal solution. Then an A* search algorithm is employed to find a set of optimal side-chain resonance assignments that best fit the NMR data. These side-chain resonance assignments are then used to resolve the NOE assignment ambiguity and compute high-resolution protein structures. Tests on five proteins show that Nasca assigns resonances for more than 90% of side-chain protons, and achieves about 80% correct assignments. The final structures computed using the NOE distance restraints assigned by Nasca have backbone RMSD 0.8–1.5 Å from the reference structures determined by traditional NMR approaches.     

Crystal structure of an RNA tertiary domain essential to HCV IRES-mediated translation initiation

The hepatitis C virus (HCV) internal ribosome entry site (IRES) RNA drives internal initiation of viral protein synthesis during host cell infection. In the tertiary structure of the IRES RNA, two helical junctions create recognition sites for direct binding of the 40S ribosomal subunit and eukaryotic initiation factor 3 (eIF3). The 2.8 A resolution structure of the IIIabc four-way junction, which is critical for binding eIF3, reveals how junction nucleotides interact with an adjacent helix to position regions directly involved in eIF3 recognition. Two of the emergent helices stack to form a nearly continuous A-form duplex, while stacking of the other two helices is interrupted by the insertion of junction residues into the helix minor groove. This distorted stack probably serves as an important recognition surface for the translational machinery.