NMR structure of an intracellular loop peptide derived from prostaglandin EP3alpha receptor

We found that a peptide (EP3a: TIKALVSRCRAKAAV) corresponding to the N-terminal site of the intracellular third loop of human prostaglandin EP3α receptor could activate G protein α-subunit directly. The activity was almost same as Mastoparan-X, a G protein activating peptide from wasp venom. The three-dimensional molecular structure of the peptide in SDS-d₂₅ micelles was determined by 2D ¹H NMR spectroscopy. The structure of EP3a consists of a positive charge cluster on the C-terminal helical site. The cluster was also found in several corresponding receptor peptides. Therefore, the positive charge cluster on the helical structure might play a crucial role in activation of G protein.     

The UAA/GAN internal loop motif: a new RNA structural element that forms a cross-strand AAA stack and long-range tertiary interactions

Analysis of aligned RNA sequences and high-resolution crystal structures has revealed a new RNA structural element, termed the UAA/GAN motif. Found in internal loops of the 23 S rRNA, as well as in RNase P RNA and group I and II introns, this six-nucleotide motif adopts a distinctive local structure that includes two base-pairs with non-canonical conformations and three conserved adenine bases, which form a cross-strand AAA stack in the minor groove. Most importantly, the motif invariably forms long-range tertiary contacts, as the AAA stack typically forms A-minor interactions and the flipped-out N nucleotide forms additional contacts that are specific to the structural context of each loop. The widespread presence of this motif and its propensity to form long-range contacts suggest that it plays a critical role in defining the architectures of structured RNAs.     

Structure of the pyrimidine-rich internal loop in the poliovirus 3'-UTR: the importance of maintaining pseudo-2-fold symmetry in RNA helices containing two adjacent non-canonical base-pairs

Formation of non-canonical base-pairs in RNA often plays a very important functional role. In addition they frequently serve as factors in stabilizing the secondary structure elements that provide the frame of large compact RNA structures. Here we describe the structure of an internal loop containing a 5'CU3'/5'UU3' non-canonical tandem base-pair motif, which is conserved within the 3'-UTR of poliovirus-like enteroviruses. Structural details reveal striking regularities of the local helix geometry, resulting from alternating geometrical adjustments, which are important for understanding and predicting stabilities and configurations of tandem non-canonical base-pairs. The C-U and U-U base-pairs severely contract the minor groove of the sugar-phosphate backbone, which might be important for protein recognition or binding to other RNA elements.     

RNA helical packing in solution: NMR structure of a 30 kDa GAAA tetraloop-receptor complex

Tertiary interactions are critical for proper RNA folding and ribozyme catalysis. RNA tertiary structure is often condensed through long-range helical packing interactions mediated by loop-receptor motifs. RNA structures displaying helical packing by loop-receptor interactions have been solved by X-ray crystallography, but not by NMR. Here, we report the NMR structure of a 30 kDa GAAA tetraloop-receptor RNA complex. In order to stabilize the complex, we used a modular design in which the RNA was engineered to form a homodimer, with each subunit containing a GAAA tetraloop phased one helical turn apart from its cognate 11-nucleotide receptor domain. The structure determination utilized specific isotopic labeling patterns (2H, 13C and 15N) and refinement against residual dipolar couplings. We observe a unique and highly unusual chemical shift pattern for an adenosine platform interaction that reveals a spectroscopic fingerprint for this motif. The structure of the GAAA tetraloop-receptor interaction is well defined solely from experimental NMR data, shows minor deviations from previously solved crystal structures, and verifies the previously inferred hydrogen bonding patterns within this motif. This work demonstrates the feasibility of using engineered homodimers as modular systems for the determination of RNA tertiary interactions by NMR.     

Solution structure of an alternate conformation of helix27 from Escherichia coli16S rRNA

Helix (H)27 of 16S ribosomal (r)RNA from Escherichia coli was dubbed the "switch helix" when mutagenesis suggested that two alternative base pair registers may have distinct functional roles in the bacterial ribosome. Although more recent genetic analyses suggest that H27 conformational switching is not required for translation, previous solution studies demonstrated that the isolated E. coli H27 can dynamically convert between the 885 and 888 conformations. Here, we have solved the nuclear magnetic resonance solution structure of a locked 888 conformation. NOE and residual dipolar coupling restraints reveal an architecture that markedly differs from that of the 885 conformation found in crystal structures of the bacterial ribosome. In place of the loop E motif that characterizes the 885 conformer and that the 888 conformer cannot adopt, we find evidence for an asymmetrical A-rich internal loop stabilized by stacking interactions among the unpaired A's. Comparison of the isolated H27 888 solution structure with the 885 crystal structure within the context of the ribosome suggests a difference in overall length of H27 that presents one plausible reason for the absence of H27 conformational switching within the sterically confining ribosome.     

Application of dipolar coupling data to the refinement of the solution structure of the sarcin-ricin loop RNA

Residual dipolar couplings can provide the long-range information that most NMR solution structures lack. The use of such data in protein structure determinations is now fairly routine, but even though these data should be much more useful for nucleic acids, their application to nucleic acid structure determination is still in its infancy. Here we present a method for producing accurate, dipolar-refined structures of nucleic acids that is more efficient than those used previously, and apply it to E73, a 29 nucleotide RNA that includes the sarcin-ricin loop from rat 28S rRNA. The results enable us to address the differences between the crystal structure of E73 and the solution structure proposed for it previously.     

Functional stabilization of an RNA recognition motif by a noncanonical N-terminal expansion

RNA recognition motifs (RRMs) constitute versatile macromolecular interaction platforms. They are found in many components of spliceosomes, in which they mediate RNA and protein interactions by diverse molecular strategies. The human U11/U12-65K protein of the minor spliceosome employs a C-terminal RRM to bind hairpin III of the U12 small nuclear RNA (snRNA). This interaction comprises one side of a molecular bridge between the U11 and U12 small nuclear ribonucleoprotein particles (snRNPs) and is reminiscent of the binding of the N-terminal RRMs in the major spliceosomal U1A and U2B″ proteins to hairpins in their cognate snRNAs. Here we show by mutagenesis and electrophoretic mobility shift assays that the β-sheet surface and a neighboring loop of 65K C-terminal RRM are involved in RNA binding, as previously seen in canonical RRMs like the N-terminal RRMs of the U1A and U2B″ proteins. However, unlike U1A and U2B″, some 30 residues N-terminal of the 65K C-terminal RRM core are additionally required for stable U12 snRNA binding. The crystal structure of the expanded 65K C-terminal RRM revealed that the N-terminal tail adopts an α-helical conformation and wraps around the protein toward the face opposite the RNA-binding platform. Point mutations in this part of the protein had only minor effects on RNA affinity. Removal of the N-terminal extension significantly decreased the thermal stability of the 65K C-terminal RRM. These results demonstrate that the 65K C-terminal RRM is augmented by an N-terminal element that confers stability to the domain, and thereby facilitates stable RNA binding.     

miRPlant: an integrated tool for identification of plant miRNA from RNA sequencing data

Background

Small RNA sequencing is commonly used to identify novel miRNAs and to determine their expression levels in plants. There are several miRNA identification tools for animals such as miRDeep, miRDeep2 and miRDeep*. miRDeep-P was developed to identify plant miRNA using miRDeep’s probabilistic model of miRNA biogenesis, but it depends on several third party tools and lacks a user-friendly interface. The objective of our miRPlant program is to predict novel plant miRNA, while providing a user-friendly interface with improved accuracy of prediction.

Result

We have developed a user-friendly plant miRNA prediction tool called miRPlant. We show using 16 plant miRNA datasets from four different plant species that miRPlant has at least a 10% improvement in accuracy compared to miRDeep-P, which is the most popular plant miRNA prediction tool. Furthermore, miRPlant uses a Graphical User Interface for data input and output, and identified miRNA are shown with all RNAseq reads in a hairpin diagram.

Conclusions

We have developed miRPlant which extends miRDeep* to various plant species by adopting suitable strategies to identify hairpin excision regions and hairpin structure filtering for plants. miRPlant does not require any third party tools such as mapping or RNA secondary structure prediction tools. miRPlant is also the first plant miRNA prediction tool that dynamically plots miRNA hairpin structure with small reads for identified novel miRNAs. This feature will enable biologists to visualize novel pre-miRNA structure and the location of small RNA reads relative to the hairpin. Moreover, miRPlant can be easily used by biologists with limited bioinformatics skills.miRPlant and its manual are freely available at http://www.australianprostatecentre.org/research/software/mirplant or http://sourceforge.net/projects/mirplant/.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-275) contains supplementary material, which is available to authorized users.     

Sheared Aanti.Aanti base pairs in a destabilizing 2 x 2 internal loop: the NMR structure of 5'(rGGCAAGCCU)2

The 5'(rGGCAAGCCU)(2) duplex contains tandem A.A pairs. The three-dimensional structure of the 5'(rGGCAAGCCU)(2) duplex was modeled by molecular dynamics and energy minimization with NMR-derived distance and dihedral angle restraints. Although the 5'(rCAAG)(2) loop is thermodynamically destabilizing by 1.1 kcal/mol, the tandem A.A pairs adopt a predominant conformation: a sheared anti-anti (A.A trans Hoogsteen/Sugar-edge) alignment similar to that observed in the crystal structure of the P4-P6 domain of the Tetrahymena thermophila intron [Cate, J. H., Gooding, A. R., Podell, E., Zhou, K., Golden, B. L., Kundrot, C. E., Cech, T. R., and Doudna, J. A. (1996) Science 273, 1678-1685]. The NMR-derived structure of the 5'(rGGCAAGCCU)(2) duplex exhibits cross-strand hydrogen bonds from N3 of A4 to an amino hydrogen of A5 and from the 2' oxygen of the A4 sugar to the other amino hydrogen of A5. An intrastrand hydrogen bond is formed from the 2' OH hydrogen of A4 to O5' of A5. The cross-strand A5 bases are stacked. The Watson-Crick G-C regions are essentially A-form. The sheared anti-anti (A.A trans Hoogsteen/Sugar-edge) alignment provides potential contact sites for tertiary interactions and, therefore, is a possible target site for therapeutics. Thus, thermodynamically destabilizing internal loops can be preorganized for tertiary interactions or ligand binding.     

A systematic comparison of three structure determination methods from NMR data: Dependence upon quality and quantity of data

Summary We have systematically examined how the quality of NMR protein structures depends on (1) the number of NOE distance constraints. (2) their assumed precision, (3) the method of structure calculation and (4) the size of the protein. The test sets of distance constraints have been derived from the crystal structures of crambin (5 kDa) and staphylococcal nuclease (17 kDa). Three methods of structure calculation have been compared: Distance Geometry (DGEOM), Restrained Molecular Dynamics (XPLOR) and the Double Iterated Kalman Filter (DIKF). All three methods can reproduce the general features of the starting structure under all conditions tested. In many instances the apparent precision of the calculated structure (as measured by the RMS dispersion from the average) is greater than its accuracy (as measured by the RMS deviation of the average structure from the starting crystal structure). The global RMS deviations from the reference structures decrease exponentially as the number of constraints is increased, and after using about 30% of all potential constraints, the crrors asymptotically approach a limiting value. Increasing the assumed precision of the constraints has the same qualitative effect as increasing the number of constraints. For comparable numbers of constraints/residue, the precision of the calculated structure is less for the larger than for the smaller protein, regardless of the method of calculation. The accuracy of the average structure calculated by Restrained Molecular Dynamics is greater than that of structures obtained by purely geometric methods (DGEOM and DIKF).     

Towards a structural classification of phosphate binding sites in protein-nucleotide complexes: an automated all-against-all structural comparison using geometric matching

A method is described for the rapid comparison of protein binding sites using geometric matching to detect similar three-dimensional structure. The geometric matching detects common atomic features through identification of the maximum common sub-graph or clique. These features are not necessarily evident from sequence or from global structural similarity giving additional insight into molecular recognition not evident from current sequence or structural classification schemes. Here we use the method to produce an all-against-all comparison of phosphate binding sites in a number of different nucleotide phosphate-binding proteins. The similarity search is combined with clustering of similar sites to allow a preliminary structural classification. Clustering by site similarity produces a classification of binding sites for the 476 representative local environments producing ten main clusters representing half of the representative environments. The similarities make sense in terms of both structural and functional classification schemes. The ten main clusters represent a very limited number of unique structural binding motifs for phosphate. These are the structural P-loop, di-nucleotide binding motif [FAD/NAD(P)-binding and Rossman-like fold] and FAD-binding motif. Similar classification schemes for nucleotide binding proteins have also been arrived at independently by others using different methods.     

Crystal structure of hypothetical protein TTHB192 from Thermus thermophilus HB8 reveals a new protein family with an RNA recognition motif-like domain

We have determined the crystal structure of hypothetical protein TTHB192 from Thermus thermophilus HB8 at 1.9 A resolution. This protein is a member of the Escherichia coli ygcH sequence family, which contains approximately 15 sequence homologs of bacterial origin. These homologs have a high isoelectric point. The crystal structure reveals that TTHB192 consists of two independently folded domains, and that each domain exhibits a ferredoxin-like fold with a four-stranded antiparallel beta-sheet packed on one side by alpha-helices. These two tandem domains face each other to generate a beta-sheet platform. TTHB192 displays overall structural similarity to Sex-lethal protein and poly(A)-binding protein fragments. These proteins have RNA binding activity which is supported by a beta-sheet platform formed by two tandem repeats of an RNA recognition motif domain with signature sequence motifs on the beta-sheet surface. Although TTHB192 does not have the same signature sequence motif as the RNA recognition motif domain, the presence of an evolutionarily conserved basic patch on the beta-sheet platform could be functionally relevant for nucleic acid-binding. This report shows that TTHB192 and its sequence homologs adopt an RNA recognition motif-like domain and provides the first testable functional hypothesis for this protein family.     

Full 1H NMR assignment of a 24-nucleotide RNA hairpin: Application of the 1H 3D-NOE/2QC experiment

The subject RNA models the binding site for the coat protein of the R17 virus, as well as the ribosome recognition sequence for the R17 replicase gene. With an RNA of this size, overlaps among the sugar protons complicate assignments of the ¹H NMR spectrum. The cross peaks that overlap significantly in 2D-NOE spectra can frequently be resolved by introducing a third, in our approach the double-quantum, frequency axis. In particular the planes in a 3D-NOE/2QC spectrum perpendicular to the 2Q axis are extremely useful, showing a highly informative repeating NOE-2Q pattern. In this experiment substantial J-coupling confers special advantages. This always occurs for geminal pairs (H5/H5 for RNA plus H2/H2 for DNA), as well as for H5/H6, for H3/H4 in sugars with substantial populations of the N-pucker, for H1/H2 for S-puckered sugars, and usually for H2/H3. For the 24-mer RNA hairpin the additional information from the 3D-NOE/2QC spectrum allowed assignment of all of the non-exchangeable protons, eliminating the need for stable-isotope labeling.     

NMR structure determination of alpha-conotoxin BuIA, a novel neuronal nicotinic acetylcholine receptor antagonist with an unusual 4/4 disulfide scaffold

We have determined a high-resolution three-dimensional structure of alpha-conotoxin BuIA, a 13-residue peptide toxin isolated from Conus bullatus. Despite its unusual 4/4 disulfide bond layout alpha-conotoxin BuIA exhibits strong antagonistic activity at alpha6/alpha3beta2beta3, alpha3beta2, and alpha3beta4 nAChR subtypes like some alpha4/7 conotoxins. alpha-Conotoxin BuIA lacks the C-terminal beta-turn present within the second disulfide loop of alpha4/7 conotoxins, having only a "pseudo omega-shaped" molecular topology. Nevertheless, it contains a functionally critical two-turn helix motif, a feature ubiquitously found in alpha4/7 conotoxins. Such an aspect seems mainly responsible for similarities in the receptor recognition profile of alpha-conotoxin BuIA to alpha4/7 conotoxins. Structural comparison of alpha-conotoxin BuIA with alpha4/7 conotoxins and alpha4/3 conotoxin ImI suggests that presence of the second helical turn portion of the two-turn helix motif in alpha4/7 and alpha4/4 conotoxins may be important for binding to the alpha3 and/or alpha6 subunit of nAChR.     

The NMR solution structure and function of RPA3313: a putative ribosomal transport protein from Rhodopseudomonas palustris

Protein function elucidation often relies heavily on amino acid sequence analysis and other bioinformatics approaches. The reliance is extended to structure homology modeling for ligand docking and protein–protein interaction mapping. However, sequence analysis of RPA3313 exposes a large, unannotated class of hypothetical proteins mostly from the Rhizobiales order. In the absence of sequence and structure information, further functional elucidation of this class of proteins has been significantly hindered. A high quality NMR structure of RPA3313 reveals that the protein forms a novel split ββαβ fold with a conserved ligand binding pocket between the first β‐strand and the N‐terminus of the α‐helix. Conserved residue analysis and protein–protein interaction prediction analyses reveal multiple protein binding sites and conserved functional residues. Results of a mass spectrometry proteomic analysis strongly point toward interaction with the ribosome and its subunits. The combined structural and proteomic analyses suggest that RPA3313 by itself or in a larger complex may assist in the transportation of substrates to or from the ribosome for further processing. Proteins 2016; 85:93–102. © 2016 Wiley Periodicals, Inc.