Classification of 3D structural character of RNA by hydrogen bond and base stacking.

We are developing a computational system to classify RNA structures by its structural character. Here, an improved grouping algorithm was introduced to the system and the base-stacking pattern (BSP) is used as a criterion for the classification in addition to hydrogen-bond pattern (HBP). 279 conformers of 15 mer RNA hairpin were classified into 89 and 36 groups by HBP and BSP, respectively, suggesting that HBP represents conformational character better than BSP.     

Force field based conformational analysis of RNA structural motifs: GNRA tetraloops and their pyrimidine relatives

The protocol of conformational analysis applied here to ribonucleotide oligomers combines conformational search in the space of torsion angles and energy minimization using the AMBER4.1 force field with a continuum treatment of electrostatic solute-solvent interactions. RNA fragments with 5′-GGGCGNNAGCCU-3′ sequences commonly fold into hairpins with four-membered loops. The combinatorial search for acceptable conformations using the MC-SYM program was restricted to loop nucleotides and yielded roughly 1500 structures being compatible with a double-stranded stem. After energy minimization by the JUMNA program (without applying any experimental constraints), these structures converged into an ensemble of 74 different conformers including 26 structures which contained the sheared G-A base pair observed in experimental studies of GNRA tetraloops. Energetic analysis shows that inclusion of solvent electrostatic effects is critically important for the selection of conformers that agree with experimentally determined structures. The continuum model accounts for solvent polarization by means of the electrostatic reaction field. In the case of GNRA loop sequences, the contributions of the reaction field shift relative stabilities towards conformations showing most of the structural features derived from NMR studies. The agreement of computed conformations with the experimental structures of GAAA, GCAA, and GAGA tetraloops suggests that the continuum treatment of the solvent represents a definitive improvement over methods using simple damping models in electrostatic energy calculations. Application of the procedure described here to the evaluation of the relative stabilities of conformers resulting from searching the conformational space of RNA structural motifs provides some progress in (non-homology based) RNA 3D-structure prediction. Received: 20 January 1999 / Revised version: 4 June 1999 / Accepted: 10 June 1999     

Escherichia coli 5S RNA A and B conformers. Characterisation by enzymatic and chemical methods

The structures of the two stable conformers of Escherichia coli 5 S RNA, the and B form, were compared. Information about the structures were obtained using the methods of limited enzymatic hydrolysis and chemical modification of accessible nucleotides. Base-specific modifications were performed for adenosines and cytidines using diethylpyrocarbonate and dimethylsulfate in combination with a strand-scission reaction at the modified site. Base-specific (RNase T1) as well as conformation-specific (nuclease S1, cobra venom nuclease) enzymes were employed for the limited enzymatic hydrolysis. Clear differences in the accessibility of the two 5 S RNA conformers to the enzymes and the chemical reagents were established and the regions with altered reactivities were localized in the 5 S RNA structure. The results are consistent with the disruption of the secondary structural interactions in helix II and partly in helices III and IV during the transition from the A to the B form. (The numbering of the helices is according to the generally accepted Fox and Woese model.) In addition some regions presumably involved in the tertiary structure are distorted. There is evidence, however, for the new formation of structural regions between two distant sites in the 5 S RNA B form. The results enable us to refine the existing 5 S RNA A-form model and provide insight into the structural dynamics that lead to the formation of the 5 S RNA B form.     

In-gel probing of individual RNA conformers within a mixed population reveals a dimerization structural switch in the HIV-1 leader

Definitive secondary structural mapping of RNAs in vitro can be complicated by the presence of more than one structural conformer or multimerization of some of the molecules. Until now, probing a single structure of conformationally flexible RNA molecules has typically relied on introducing stabilizing mutations or adjusting buffer conditions or RNA concentration. Here, we present an in-gel SHAPE (selective 2′OH acylation analysed by primer extension) approach, where a mixed structural population of RNA molecules is separated by non-denaturing gel electrophoresis and the conformers are individually probed within the gel matrix. Validation of the technique using a well-characterized RNA stem-loop structure, the HIV-1 trans-activation response element, showed that authentic structure was maintained and that the method was accurate and highly reproducible. To further demonstrate the utility of in-gel SHAPE, we separated and examined monomeric and dimeric species of the HIV-1 packaging signal RNA. Extensive differences in acylation sensitivity were seen between monomer and dimer. The results support a recently proposed structural switch model of RNA genomic dimerization and packaging, and demonstrate the discriminatory power of in-gel SHAPE.     

Two distinct conformations of rat liver ribosomal 5S RNA.

Three different conformers of rat liver 5S ribosomal RNA were investigated by partial nuclease cleavage technique using S1 nuclease and cobra venom endoribonuclease (CVE) as conformational probes. Urea-treated and renatured 5S RNA co-migrate on non-denaturing gels, but exhibit distinct differences in their nuclease cleavage patterns. The most prominent differences in S1 nuclease and CVE accessibility of these conformers are located in region 30-50 and around nucleotides 70 and 90. The third form of 5S RNA with higher electrophoretic mobility was generated by EDTA treatment. The cleavage patterns of this 5S RNA conformer are similar to that characteristic for the renatured 5S RNA. The results demonstrate the difference in secondary structure and possibly different tertiary base-pairing interactions of 5S RNA conformers.     

Lipid-dependent Generation of Dual Topology for a Membrane Protein

The mechanism by which membrane proteins exhibit structural and functional duality in the same membrane or different membranes is unknown. We posit that such duality is determined by both the protein sequence and the membrane lipid composition wherein a spatial or temporal change in the latter can result in a post-assembly change in protein structure and function. To investigate whether co-existence of multiple topological conformers is dependent on the membrane lipid composition, we determined the topological organization of lactose permease in an Escherichia coli model cell system in which phosphatidylethanolamine membrane content can be systematically varied. At intermediate levels of phosphatidylethanolamine a mixture of native and topologically mis-oriented conformers co-existed. There was no threshold level of phosphatidylethanolamine determining a sharp transition from one conformer to the other. Co-existing conformers were not in rapid equilibrium at a static lipid composition indicating that duality of topology is established during an early folding step. Depletion of intermediate levels of phosphatidylethanolamine after final protein assembly resulted in complete mis-orientation of the native conformer. Combined with previous results, such topological dynamics are reversible in both directions. We propose a thermodynamically based model for how lipid-protein interactions can result in a mixed topological organization and how changes in lipid composition can result in changes in the ratio of topologically distinct conformers of proteins. These observations demonstrate a potential lipid-dependent biological switch for generating dynamic structural and functional heterogeneity for a protein within the same membrane or between different membranes in more complex eukaryotic cells.     

A complete conformational map for RNA.

A simple stereochemical framework for understanding RNA structure has remained elusive to date. We present a comprehensive conformational map for two nucleoside-5',3'-diphosphates and for a truncated dinucleotide derived from a grid search of all potential conformers using hard sphere steric exclusion criteria to define allowed conformers. The eight-dimensional conformational space is presented as a series of two-dimensional projections. These projections reveal several well-defined allowed and disallowed regions which correlate well with data obtained from X-ray crystallography of both large and small RNA molecules. Furthermore, the two-dimensional projections show that consecutive and ribose ring-proximal torsion angles are interdependent, while more distant torsion angles are not. Remarkably, using steric criteria alone, it is possible to generate a predictive conformational map for RNA.     

Semi-automated update and cleanup of structural RNA alignment databases.

We have developed a series of programs which assist in maintenance of structural RNA databases. A main program BLASTs the RNA database against GenBank and automatically extends and realigns the sequences to include the entire range of the RNA query sequences. After manual update of the database, other programs can examine base pair consistency and phylogenetic support. The output can be applied iteratively to refine the structural alignment of the RNA database. Using these tools, the number of potential misannotations per sequence was reduced from 20 to 3 in the Signal Recognition Particle RNA database. AVAILABILITY: A quick-server and programs are available at http://www.bioinf.au.dk/rnadbtool/     

Hyperreactivity of adenines and conformational flexibility of a translational repression site.

We have used a diethylpyrocarbonate (DEPC) modification [(1976) Prog. Nucl. Acids Res. 16, 189-262] to probe the accessibility of adenines essential for coat protein binding in the MS2 translational operator [(1983) Biochemistry 22, 2601-2610, 2610-2615, 4723-4730; (1987) Biochemistry 26, 1563-1568]. The essential adenines are apparently hyperreactive with this reagent relative to other sites within the same molecule. Variation of ionic strength, pH and divalent cation concentrations reveal the existence of two distinct conformers of the RNA operator as judged by DEPC reactivity. We propose that the hyperreactivity observed is due to the participation of neighbouring bases in the DEPC modification reaction and can be used as a novel structural probe.     

NMR evidence for the existence of two native conformations of 5S RNA

NMR spectra of the non-exchangeable protons in 5S RNA from E. coli show the existence of two distinct conformers of the molecule which meet the operational definition of "A form" or native 5S RNA. Both are easily distinguished spectroscopically from denatured, "B form" 5S RNA. The conditions which interconvert the two A form conformers strongly suggest that the transition between them gives rise to the low temperature optical melting transition first reported in 5S RNA by Kao and Crothers (1).     

Non-canonical base pairs and higher order structures in nucleic acids: crystal structure database analysis

Non-canonical base pairs, mostly present in the RNA, often play a prominent role towards maintaining their structural diversity. Higher order structures like base triples are also important in defining and stabilizing the tertiary folded structure of RNA. We have developed a new program BPFIND to analyze different types of canonical and non-canonical base pairs and base triples involving at least two direct hydrogen bonds formed between polar atoms of the bases or sugar O2' only. We considered 104 possible types of base pairs, out of which examples of 87 base pair types are found to occur in the available RNA crystal structures. Analysis indicates that approximately 32.7% base pairs in the functional RNA structures are non-canonical, which include different types of GA and GU Wobble base pairs apart from a wide range of base pair possibilities. We further noticed that more than 10.4% of these base pairs are involved in triplet formation, most of which play important role in maintaining long-range tertiary contacts in the three-dimensional folded structure of RNA. Apart from detection, the program also gives a quantitative estimate of the conformational deformation of detected base pairs in comparison to an ideal planar base pair. This helps us to gain insight into the extent of their structural variations and thus assists in understanding their specific role towards structural and functional diversity.     

RNABase: an annotated database of RNA structures

RNABase is a unified database of all three-dimensional structures containing RNA deposited in either the Protein Data Bank (PDB) or Nucleic Acid Data Base (NDB). For each structure, RNABase contains a brief summary as well as annotation of conformational parameters, identification of possible model errors, Ramachandran-style conformational maps and classification of ribonucleotides into conformers. These same analyses can also be performed on structures submitted by users. To facilitate access, structures are automatically placed into a variety of functional and structural categories, including: ribozymes, pseudoknots, etc. RNABase can be freely accessed on the web at http://www.rnabase.org. We are committed to maintaining this database indefinitely.     

Conformational analysis of octa- and tetrabromo tetraphenylporphyrins and their Ni(II) and Tb(III) complexes

Molecular mechanics (MM) calculations were used to analyze the puckering of metalloporphyrins as a function of metal ion size and the position of substituents on the porphyrin periphery, on a three series of octa- and tetrabromo tetraphenylporphyrins: without metal, and with Ni(II), and Tb(III) as representative small and large metal ions, respectively. Molecular energy optimization calculations were carried out using the Consistent Force Field (CFF) program, with the parameters developed previously and new parameters for bromine atom. Normal-coordinate structural decomposition (NSD) analysis was performed on the equilibrium structures obtained by MM calculations. The conformers are also stereochemically characterized, compared with available X-ray structures and with the conformers obtained in our previous MM study using chloro instead of bromo beta-pyrrole substituents.     

Structural versatility of bovine ribonuclease A. Distinct conformers of trimeric and tetrameric aggregates of the enzyme.

Lyophilization of bovine ribonuclease A (RNase A; Sigma, type XII-A) from 40% acetic acid solutions leads to the formation of approximately 14 aggregated species that can be separated by ion-exchange chromatography. Several aggregates were identified, including two variously deamidated dimeric subspecies, two distinct trimeric and two distinct tetrameric RNase A conformers, besides the two forms of dimer characterized previously [Gotte, G. & Libonati, M. (1998) Two different forms of aggregated dimers of ribonuclease A. Biochim. Biophys. Acta 1386, 106-112]. We also have possible evidence for the existence of two forms of pentameric RNase A. The two forms of trimers and tetramers are characterized by: (a) slightly different gel filtration patterns; (b) different retention times in ion-exchange chromatography; and (c) different mobilities in cathodic gel electrophoresis under nondenaturing conditions. Therefore, they appear to have distinct structural organizations responsible for a different availability of their positively charged amino acid residues. All RNase A oligomers, in particular the two distinct trimeric and tetrameric conformers, degrade poly(A).poly(U), viral double-stranded RNA and polyadenylate with a catalytic efficiency that is in general higher for the more basic species. On the contrary, the activity of the RNase A oligomers, from dimer to pentamer, on yeast RNA and poly(C) (Kunitz assay) is lower than that of monomeric RNase A.     

In vitro assembly of Ebola virus nucleocapsidlike complex expressed in E. coli

Ebola virus (EBOV) harbors an RNA genome encapsidated by nucleoprotein (NP) along with other viral proteins to form a nucleocapsid complex. Previous Cryoeletron tomography and biochemical studies have shown the helical structure of EBOV nucleocapsid at nanometer resolution and the first 450 amino-acid of NP (NPΔ451–739) alone is capable of forming a helical nucleocapsid-like complex (NLC). However, the structural basis for NP-NP interaction and the dynamic procedure of the nucleocapsid assembly is yet poorly understood. In this work, we, by using an E. coli expression system, captured a series of images of NPΔ451–739 conformers at different stages of NLC assembly by negative-stain electron microscopy, which allowed us to picture the dynamic procedure of EBOV nucleocapsid assembly. Along with further biochemical studies, we showed the assembly of NLC is salt-sensitive, and also established an indispensible role of RNA in this process. We propose the diverse modes of NLC elongation might be the key determinants shaping the plasticity of EBOV virions. Our findings provide a new model for characterizing the self-oligomerization of viral nucleoproteins and studying the dynamic assembly process of viral nucleocapsid in vitro.