Algorithm independent properties of RNA secondary structure predictions

Algorithms predicting RNA secondary structures based on different folding criteria – minimum free energies (mfe), kinetic folding (kin), maximum matching (mm) – and different parameter sets are studied systematically. Two base pairing alphabets were used: the binary GC and the natural four-letter AUGC alphabet. Computed structures and free energies depend strongly on both the algorithm and the parameter set. Statistical properties, such as mean number of base pairs, mean numbers of stacks, mean loop sizes, etc., are much less sensitive to the choice of parameter set and even of algorithm. Some features of RNA secondary structures, such as structure correlation functions, shape space covering and neutral networks, seem to depend only on the base pairing logic (GC or AUGC alphabet). Received: 16 May 1996 / Accepted: 10 July 1996     

Neutral adaptation of the genetic code to double-strand coding

Summary We lay new foundations to the hypothesis that the genetic code is adapted to evolutionary retention of information in the antisense strands of natural DNA/RNA sequences. In particular, we show that the genetic code exhibits, beyond the neutral replacement patterns of amino acid substitutions, optimal properties by favoring simultaneous evolution of proteins encoded in DNA/RNA sense-antisense strands. This is borne out in the sense-antisense transformations of the codons of every amino acid which target amino acids physicochemically similar to each other. Moreover, silent mutations in the sense strand generate conservative ones in its antisense counterpart and vice versa. Coevolution of proteins coded by complementary strands is shown to be a definite possibility, a result which does not depend on any physical interaction between the coevolving proteins. Likewise, the degree to which the present genetic code is dedicated to evolutionary sense-antisense tolerance is demonstrated by comparison with many randomized codes. Double-strand coding is quantified from an information-theoretical point of view.     

Automated identification of RNA 3D modules with discriminative power in RNA structural alignments

Recent progress in predicting RNA structure is moving towards filling the ‘gap’ in 2D RNA structure prediction where, for example, predicted internal loops often form non-canonical base pairs. This is increasingly recognized with the steady increase of known RNA 3D modules. There is a general interest in matching structural modules known from one molecule to other molecules for which the 3D structure is not known yet. We have created a pipeline, metaRNAmodules, which completely automates extracting putative modules from the FR3D database and mapping of such modules to Rfam alignments to obtain comparative evidence. Subsequently, the modules, initially represented by a graph, are turned into models for the RMDetect program, which allows to test their discriminative power using real and randomized Rfam alignments. An initial extraction of 22 495 3D modules in all PDB files results in 977 internal loop and 17 hairpin modules with clear discriminatory power. Many of these modules describe only minor variants of each other. Indeed, mapping of the modules onto Rfam families results in 35 unique locations in 11 different families. The metaRNAmodules pipeline source for the internal loop modules is available at http://rth.dk/resources/mrm.     

Endothelial Transient Receptor Potential Conical Channel (TRPC)-3 Activation Induces Vasogenic Edema Formation in the Rat Piriform Cortex Following Status Epilepticus

Transient receptor potential canonical channel (TRPC) is a nonselective cation channel permeable to Ca²⁺, which express in many cell types, including neurons. However the alterations in TRPC receptor expressions in response to status epilepticus (SE) have not been explored. Therefore, the present study was designated to elucidate the roles of TRPC3 in neuronal death and vasogenic edema within the rat piriform cortex (PC) following SE. In non-SE animals, TRPC3 immunoreactivity was abundantly detected in the PC. Following SE, TRPC3 immunoreactivity was increased in neurons. Furthermore, TRPC3 expression was detected in endothelial cells that did not contain it in non-SE animals. Loss of SMI-71 (a blood–brain barrier antigen) immunoreactivity was also observed in TRPC3 positive endothelial cells. In addition, FJB positive neurons and vasogenic edema were noticeably detected in the PC. To directly determine whether TRPC3 activation is correlated to SE-induced vasogenic edema formation and neuronal damages in the PC, the effect of Pyr-3 (a TRPC3 antagonist) on SE-induced insults were investigated. Pyr-3 infusion effectively attenuated vasogenic edema in the PC as compared to the vehicle. Therefore, our findings indicate that TRPC3 activation/overexpression induced by SE may involve BBB disruption and neuronal damages in the rat PC following SE. Therefore, the present study was TRPC3 may play an important role in SE-induced vasogenic edema formation through BBB disruptions in the rat PC.     

Structural and functional fingerprint of the mitochondrial ATP-binding cassette transporter Mdl1 from Saccharomyces cerevisiae

The ATP-binding cassette half-transporter Mdl1 from Saccharomyces cerevisiae has been proposed to be involved in the quality control of misassembled respiratory chain complexes by exporting degradation products generated by the m-AAA proteases from the matrix. Direct functional or structural data of the transport complex are, however, not known so far. After screening expression in various hosts, Mdl1 was overexpressed 100-fold to 1% of total mitochondrial membrane protein in S. cerevisiae. Based on detergent screens, Mdl1 was solubilized and purified to homogeneity. Mdl1 showed a high binding affinity for MgATP (Kd = 0.26 microm) and an ATPase activity with a Km of 0.86 mm (Hill coefficient of 0.98) and a turnover rate of 2.6 ATP/s. Mutagenesis of the conserved glutamate downstream of the Walker B motif (E599Q) or the conserved histidine of the H-loop (H631A) abolished ATP hydrolysis, whereas ATP binding was not affected. Mdl1 reconstituted into liposomes showed an ATPase activity similar to the solubilized complex. By single particle electron microscopy, a first three-dimensional structure of the mitochondrial ATP-binding cassette transporter was derived at 2.3-nm resolution, revealing a homodimeric complex in an open conformation.     

Alignment of RNA base pairing probability matrices

MOTIVATION: Many classes of functional RNA molecules are characterized by highly conserved secondary structures but little detectable sequence similarity. Reliable multiple alignments can therefore be constructed only when the shared structural features are taken into account. Since multiple alignments are used as input for many subsequent methods of data analysis, structure-based alignments are an indispensable necessity in RNA bioinformatics. RESULTS: We present here a method to compute pairwise and progressive multiple alignments from the direct comparison of base pairing probability matrices. Instead of attempting to solve the folding and the alignment problem simultaneously as in the classical Sankoff's algorithm, we use McCaskill's approach to compute base pairing probability matrices which effectively incorporate the information on the energetics of each sequences. A novel, simplified variant of Sankoff's algorithms can then be employed to extract the maximum-weight common secondary structure and an associated alignment. AVAILABILITY: The programs pmcomp and pmmulti described in this contribution are implemented in Perl and can be downloaded together with the example datasets from http://www.tbi.univie.ac.at/RNA/PMcomp/. A web server is available at http://rna.tbi.univie.ac.at/cgi-bin/pmcgi.pl     

Prediction of consensus RNA secondary structures including pseudoknots

Most functional RNA molecules have characteristic structures that are highly conserved in evolution. Many of them contain pseudoknots. Here, we present a method for computing the consensus structures including pseudoknots based on alignments of a few sequences. The algorithm combines thermodynamic and covariation information to assign scores to all possible base pairs, the base pairs are chosen with the help of the maximum weighted matching algorithm. We applied our algorithm to a number of different types of RNA known to contain pseudoknots. All pseudoknots were predicted correctly and more than 85 percent of the base pairs were identified.     

Vienna RNA secondary structure server

The Vienna RNA secondary structure server provides a web interface to the most frequently used functions of the Vienna RNA software package for the analysis of RNA secondary structures. It currently offers prediction of secondary structure from a single sequence, prediction of the consensus secondary structure for a set of aligned sequences and the design of sequences that will fold into a predefined structure. All three services can be accessed via the Vienna RNA web server at http://rna.tbi.univie.ac.at/.     

Horizontal Transfer of Archaeal Genes into the Deinococcaceae: Detection by Molecular and Computer-Based Approaches

Members of the Deinococcaceae (e.g., Thermus, Meiothermus, Deinococcus) contain A/V-ATPases typically found in Archaea or Eukaryotes which were probably acquired by horizontal gene transfer. Two methods were used to quantify the extent to which archaeal or eukaryotic genes have been acquired by this lineage. Screening of a Meiothermus ruber library with probes made against Thermoplasma acidophilum DNA yielded a number of clones which hybridized more strongly than background. One of these contained the prolyl tRNA synthetase (RS) gene. Phylogenetic analysis shows the M. ruber and D. radiodurans prolyl RS to be more closely related to archaeal and eukaryal forms of this gene than to the typical bacterial type. Using a bioinformatics approach, putative open reading frames (ORFs) from the prerelease version of the D. radiodurans genome were screened for genes more closely related to archaeal or eukaryotic genes. Putative ORFs were searched against representative genomes from each of the three domains using automated BLAST. ORFs showing the highest matches against archaeal and eukaryotic genes were collected and ranked. Among the top-ranked hits were the A/V-ATPase catalytic and noncatalytic subunits and the prolyl RS genes. Using phylogenetic methods, ORFs were analyzed and trees assessed for evidence of horizontal gene transfer. Of the 45 genes examined, 20 showed topologies in which D. radiodurans homologues clearly group with eukaryotic or archaeal homologues, and 17 additional trees were found to show probable evidence of horizontal gene transfer. Compared to the total number of ORFs in the genome, those that can be identified as having been acquired from Archaea or Eukaryotes are relatively few (approximately 1%), suggesting that interdomain transfer is rare.