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     

A probability matrix for the identification of vibrios

A probability matrix for computer-assisted identification of vibrios has been constructed, based on the API 20E system. Data were gathered from 173 strains representing 31 taxa of vibrios and related organisms, from a variety of sources. The matrix was tested internally by four statistical programs. Program OVERMAT tested the separation and program MOSTTYP the discretion and homogeneity of the taxa. Most of the taxa were satisfactory but a few were less so; reasons for this are discussed. Program CHARSEP and program DIACHAR tested the separation and diagnostic values, respectively, of the characters used. The overall test error was 4.5%. The matrix was assessed externally by its performance in the identification of vibrio-like strains isolated from freshwater. Of 243 wild strains, 79.4% were identified with ten taxa, with a Willcox score of greater than or equal to 0.99.     

A probability matrix for the identification of vibrios

A probability matrix for computer-assisted identification of vibrios has been constructed, based on the API 20E system. Data were gathered from 173 strains representing 31 taxa of vibrios and related organisms, from a variety of sources. The matrix was tested internally by four statistical programs. Program OVERMAT tested the separation and program MOSTTYP the discretion and homogeneity of the taxa. Most of the taxa were satisfactory but a few were less so; reasons for this are discussed. Program CHARSEP and program DIACHAR tested the separation and diagnostic values, respectively, of the characters used. The overall test error was 4–5%. The matrix was assessed externally by its performance in the identification of vibrio-like strains isolated from freshwater. Of 243 wild strains, 79–4% were identified with ten taxa, with a Willcox score of ζ 0–99.     

Using substitution matrices to estimate probability distributions for biological sequences.

Accurately estimating probabilities from observations is important for probabilistic-based approaches to problems in computational biology. In this paper we present a biologically-motivated method for estimating probability distributions over discrete alphabets from observations using a mixture model of common ancestors. The method is an extension of substitution matrix-based probability estimation methods. In contrast to previous such methods, our method has a simple Bayesian interpretation and has the advantage over Dirichlet mixtures that it is both effective and simple to compute for large alphabets. The method is applied to estimate amino acid probabilities based on observed counts in an alignment and is shown to perform comparably to previous methods. The method is also applied to estimate probability distributions over protein families and improves protein classification accuracy.     

Software for the development and evaluation of probabilistic identification matrices

A number of programs are described for the development and evaluationof probabilistic identification matrices for use with computer-assistedidentification. The program BEST reads an initial matrix ofper cent probabilities for all the binary characters examinedduring a cluster analysis and determines the most useful setof tests for distinguishing taxa in the matrix. Program RESORTcreates from the initial matrix an identification matrix inwhich the order of the tests may be different or the numberof tests reduced. A printed version of a matrix can be producedby MATPRJNT, which creates tables giving the per cent probabilityof a positive test result and the test results presented as‘–’, ‘v’ and ‘+’ dependingon a user-specified threshold. Program IDSC evaluates an identificationmatrix by calculating the best identification score for eachtaxon in the matrix using the expected result for each test.The programs were written in FORTRAN 77 and can be run on anymicrocomputer using the PC/MS-DOS operating system. Received on April 27, 1990; accepted on November 5, 1990     

New polyacrylamide matrices for drift-free isoelectric focusing

The major cause for pH gradient decay and cathodic drift during isoelectric focusing in polyacrylamide gels has been found to be electroendo-osmotic flow generated by fixed charges in the gel matrix. These charges have the following causes: (a) trace impurities of acrylic acid in the co-monomers; (b) covalent incorporation of catalysts (persulphate and riboflavin 5'-phosphate) as terminal groups in polyacrylamide chains; (c) hydrolysis of amide groups to acrylic acid in the gel layer underneath the cathodic filter paper strip. The result of these fixed negative charges in the matrix is a movement of counter-ions with hydration water towards the cathode (i.e. electroendo-osmosis) with concomitant drift of pH gradient and focused protein zones in the same direction. It is impossible to cure the cathodic drift by increasing the pH of the anolyte, or decreasing the pH of the catholyte, or both, as previously suggested. One way to reduce the cathodic drift efficiently is to incorporate covalently in the matrix tertiary or quaternary groups (3-dimethylaminopropylmethacrylamide) in stoichiometric amounts as compared with the negative charges.This ‘balanced’ polyacrylamide displays zero drift for at least 5000V·h, which is considered to be an ample time for equilibrium separation of protein species in isoelectric focusing.     

New Metabolites of Streptomyces alboniger

A mutant unable to grow under acidic conditions was isolated from the acidophilic heterotrophic bacterium Acidiphilium facilis 24R. The growth of the mutant could be fully restored by the addition of spermidine or lysine at the concentration of 100 μm. The HPLC analysis of polyamine composition showed that spermidine and putrescine were major polyamine components in the parental strain. In the mutant strain, putrescine was replaced by cadaverine. It was found that some polyamines in the cells were conjugated with the other cell components. The growth of the bacterium in the medium below pH 4.5 was inhibited in the presence of α-methylornithine or methylglyoxal-bis(guanylhydrazone), which are inhibitors of rate-limiting enzymes involved in the biosynthesis of polyamines. The growth of the bacterium that had been inhibited in the presence of inhibitors could be fully restored by the addition of putrescine or spermidine. On the basis of these results, it was concluded that polyamines have a significant role in the growth of Acidiphilium facilis 24R under acidic conditions.     

Rapid identification of Streptomyces isolates by MALDI-TOF MS

The recent emergence of multidrug-resistant bacteria over the last decade has led to a renewal in the discovery of new antimicrobial drugs. Streptomyces members are practically unlimited sources of new antibiotics. However, the identification of Streptomyces species is difficult and time-consuming. Therefore, there is a need for alternative methods for their rapid identification. In this study, an efficient protocol of identification using Matrix-Assisted Laser Desorption Ionization Time-Of-Flight Mass Spectrometry (MALDI-TOF MS) was developed and applied for the rapid identification of Streptomyces isolates from the El Kala lakes in northeastern Algeria. A collection of 48 Streptomyces isolates were used for this study. The optimized procedure allowed us to obtain specific and reproducible protein spectra for each Streptomyces isolate tested. The spectra generated were used to build a preliminary local database based on their initial 16S rRNA identification. The blind test used for the identification of 20 Streptomyces strains already available in our created database and 20 unknown Streptomyces isolates showed that all (100%) of the Streptomyces strains listed in the database were rapidly (<30 min) identified with high scores of up to 2.8. Here, for the first time we showed that MALDI-TOF MS could be used as a cost-effective tool for the rapid identification of Streptomyces isolates.     

Genome analyses of Streptomyces peucetius ATCC 27952 for the identification and comparison of cytochrome P450 complement with other Streptomyces

We have determined the genome sequence of 8.7 Mb chromosome of Streptomyces peucetius ATCC 27952, which produces clinically important anthracycline chemotherapeutic agents of the polyketide class of antibiotics, daunorubicin and doxorubicin. The cytochrome P450 (CYP) superfamily is represented by 19 sequences in the S. peucetius. Among those, 15 code for functional genes, whereas the remaining four are pseudo genes. CYPs from S. peucetius are phylogenetically close to those of Streptomyces amermitilis. Four CYPs are associated with modular PKS of avermectin and two with doxorubicin biosynthetic gene cluster. CYP252A1 is the new family found in S. peucetius, which shares 38% identity to CYP51 from Streptomyces coelicolor A3 (2). Nine CYPs from S. peucetius are found in the cluster containing various regulatory genes including rar operon, conserved in S. coelicolor A3 (2) and Streptomyces griseus. Although two ferredoxins and four ferredoxin reductases have been identified so far, only one ferredoxin reductase was found in the cluster of CYP147F1 in S. peucetius. To date, 174 CYPs have been described from 45 Streptomyces species in all searchable databases. However, only 18 CYPs are clustered with ferredoxin. The comparative study of cytochrome P450s, ferredoxins, and ferredoxin reductases should be useful for the future development and manipulation of antibiotic biosynthetic pathways.     

A probability matrix for the identification of species of Vibrio and related genera

A matrix for the probabilistic identification of species of Vibrio and related genera has been constructed using the data from 1091 strains collected throughout the world and classified. Thirty-eight phenons are included in the matrix, 31 of these represent previously identified species or biovars and seven represent phenons which could not be identified and may represent new species. The identification matrix incorporates 81 characters although a subset of 30 tests can be used to distinguish the 38 phenons from each other. The additional 51 tests were included to assist the identification of some strains for which the initial 30 tests were inadequate. No significant cluster overlap was found at the 5% level and the identification score for the Hypothetical Median Organism of each cluster exceeded 0.9999 in all cases.     

A probability matrix for the identification of species of Vibrio and related genera

A matrix for the probabilistic identification of species of Vibrio and related genera has been constructed using the data from 1091 strains collected throughout the world and classified. Thirty-eight phenons are included in the matrix, 31 of these represent previously identified species or biovars and seven represent phenons which could not be identified and may represent new species. The identification matrix incorporates 81 characters although a subset of 30 tests can be used to distinguish the 38 phenons from each other. The additional 51 tests were included to assist the identification of some strains for which the initial 30 tests were inadequate. No significant cluster overlap was found at the 5% level and the identification score for the Hypothetical Median Organism of each cluster exceeded 0.9999 in all cases.     

A probability matrix for the identification of campylobacters, helicobacters and allied taxa

A probabilistic identification matrix for campylobacteria, comprising 67 phenotypic characters and 37 taxa, is described. The accuracy and integrity of the matrix was evaluated using established computer-assisted methods. Certain taxa (for example, Campylobacter concisus and Camp. gracilis ) demonstrated significant phenotypic diversity; previous data corroborated these findings. Differentiation between a few pairs of taxa proved difficult, although discriminatory characteristics were noted in each of these cases. The results indicate that most campylobacteria can be identified accurately and objectively with phenotypic tests when probabilistic methods of data assessment are employed.     

New matrices for the purification of pectinases by affinity chromatography

Polygalacturonic acid was used as a ligand in the affinity technique for pectinases purification from the filtrate of Aspergillus niger 71 culture. For this purpose four matrices were examined, namely, alkylamine controlled porous glass (CPG), alkylamine silica gel as well as keratin or polyamide coated silica gel. Good results of pectinase purification was obtained on silanized CPG or keratin coated silica gel supports.     

Modelling and parameter identification for batch fermentations with Streptomyces tendae under phosphate limitation

Summary A simple structured model for the dynamics of phosphate-limited batch fermentations with Streptomyces tendae is presented. The model describes the influence of intracellular phosphate storage upon the growth behav our of the culture. The development of the model takes into account the possible internal regulatory processes of phosphate metabolism. These complex biochemical pathways are summarized with regard to rate-limiting steps to obtain relatively simple model equations. The model parameters are fitted to the experimental data with an identification programme based on the sequential quadratic programming algorithm. Modifications in this algorithm yield a good performance for this application. With respect to the sensitivity of the model parameters, a feedback on the modelling is given. After several loops of modelling and identification, a model was achieved that fits to a set of batch fermentations. Furthermore the simulations show that RNA measurements of some recent fermentations can be interpreted by the simulated internal state variable and that there is evidence for RNA as an intracellular phosphate reserve. Offprint requests to: K.-P. Kuhn