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.     

Exclusion of repetitive DNA elements from gnathostome Hox clusters

Despite their homology and analogous function, the Hox gene clusters of vertebrates and invertebrates are subject to different constraints on their structural organization. This is demonstrated by a drastically different distribution of repetitive DNA elements in the Hox cluster regions. While gnathostomes have a strong tendency to exclude repetitive DNA elements from the inside of their Hox clusters, no such trend can be detected in the Hox gene clusters of protostomes. Repeats "invade" the gnathostome Hox clusters from the 5' and 3' ends while the core of the clusters remains virtually free of repetitive DNA. This invasion appears to be correlated with relaxed constraints associated with gene loss after cluster duplications.     

Insights into the genome of the enteric bacterium Escherichia blattae: cobalamin (B12) biosynthesis, B12-dependent reactions, and inactivation of the gene region encoding B12-dependent glycerol dehydratase by a new mu-like prophage

The enteric bacterium Escherichia blattae has been analyzed for the presence of cobalamin (B12) biosynthesis and B12-dependent pathways. Biochemical studies revealed that E. blattae synthesizes B12 de novo aerobically and anaerobically. Genes exhibiting high similarity to all genes of Salmonella enterica serovar Typhimurium, which are involved in the oxygen-independent route of B12 biosynthesis, were present in the genome of E. blattae DSM 4481. The dha regulon encodes the key enzymes for the anaerobic conversion of glycerol to 1,3-propanediol, including coenzyme B12-dependent glycerol dehydratase. E. blattae DSM 4481 lacked glycerol dehydratase activity and showed no anaerobic growth with glycerol, but the genome of E. blattae DSM 4481 contained a dha regulon. The E. blattaedha regulon is unusual, since it harbors genes for two types of dihydroxyacetone kinases. The major difference to dha regulons of other enteric bacteria is the inactivation of the dehydratase-encoding gene region by insertion of a 33,339-bp prophage (MuEb). Sequence analysis revealed that MuEb belongs to the Mu family of bacteriophages. The E. blattae strains ATCC 33429 and ATCC 33430 did not contain MuEb. Accordingly, both strains harbored an intact dehydratase-encoding gene region and fermented glycerol. The properties of the glycerol dehydratases and the correlating genes (dhaBCE) of both strains were similar to other B12-dependent glycerol and diol dehydratases, but both dehydratases exhibited the highest affinity for glycerol of all B12-dependent dehydratases characterized so far. In addition to the non-functional genes encoding B12-dependent glycerol dehydratase, the genome of E. blattae DSM 4481 contained the genes for only one other B12-dependent enzyme, the methylcobalamin-dependent methionine synthase.     

Inner kinetochore of the pathogenic yeast Candida glabrata

The human pathogenic yeast Candida glabrata is the second most common Candida pathogen after Candida albicans, causing both bloodstream and mucosal infections. The centromere (CEN) DNA of C. glabrata (CgCEN), although structurally very similar to that of Saccharomyces cerevisiae, is not functional in S. cerevisiae. To further examine the structure of the C. glabrata inner kinetochore, we isolated several C. glabrata homologs of S. cerevisiae inner kinetochore protein genes, namely, genes for components of the CBF3 complex (Ndc10p, Cep3p, and Ctf13p) and genes for the proteins Mif2p and Cse4p. The amino acid sequence identities of these proteins were 32 to 49% relative to S. cerevisiae. CgNDC10, CgCEP3, and CgCTF13 are required for growth in C. glabrata and are specifically found at CgCEN, as demonstrated by chromatin immunoprecipitation experiments. Cross-complementation experiments revealed that the isolated genes, with the exception of CgCSE4, are species specific and cannot functionally substitute for the corresponding genes in S. cerevisiae deletion strains. Likewise, the S. cerevisiae CBF3 genes NDC10, CEP3, and CTF13 cannot functionally replace their homologs in C. glabrata CBF3 deletion strains. Two-hybrid analysis revealed several interactions between these proteins, all of which were previously reported for the inner kinetochore proteins of S. cerevisiae. Our findings indicate that although many of the inner kinetochore components have evolved considerably between the two closely related species, the organization of the C. glabrata inner kinetochore is similar to that in S. cerevisiae.     

Mal de Meleda (MDM) caused by mutations in the gene for SLURP-1 in patients from Germany,Turkey, Palestine,and the United Arab Emirates

Mal de Meleda (MDM) or keratosis palmoplantaris transgrediens of Siemens is an autosomal recessive skin disorder characterized by diffuse palmoplantar keratoderma (PPK) and transgressive keratosis with an onset in early infancy. There is no associated involvement of other organs; however, a spectrum of clinical presentations with optional and variable features has been described. Mutations in the ARS (component B)-81/s gene ( LY6LS) on chromosome 8q24-qter, which encodes SLURP-1, have recently been identified in patients with MDM. Here, we have analyzed four MDM families for mutations in SLURP-1. In a large Palestinian pedigree with multiple consanguinity, patients are homozygous for a new mutation that substitutes an arginine for a conserved glycine residue at position 86. A different mutation in Turkish patients results in the same amino acid exchange. Some remarkable similarities are seen in the clinical picture of patients from both families. Patients of an Emirati Bedouin family have a homozygous alteration of the translation initiation codon. In a German family with no known consanguinity, we have shown pseudodominant inheritance. Three affected children and their affected mother are homozygous for the missense mutation W15R. Our findings indicate that the MDM type of transgressive PPK is caused by SLURP-1 mutations in patients from various origins and demonstrate allelic heterogeneity for mutations in SLURP-1.     

Serum-free generation of antigen presenting cells from acute myeloid leukaemic blasts for active specific immunisation

PURPOSE: Immunotherapy holds promise as a new strategy for the eradication of residual cells in acute myeloid leukaemia (AML). Leukaemic antigen presenting cells (APCs) combining optimal antigen presentation and tumour antigenicity could be used as potent T cell activators. For clinical purposes it is desirable to culture APCs under serum-free conditions. Therefore, we compared morphological, immunophenotypical and functional outcome of the serum-free culture of AML-APCs to their serum-enriched culture. METHODS: AML blasts (n=19) were cultured in the presence of either a cytokine mix or calcium ionophore (CI) for 14 and 2 days, respectively, in FCS-containing medium (FCS), StemSpan serum-free medium (SP) and CellGro serum-free medium (CG). After culture relative yields were calculated and immunophenotypic analysis of APC markers was performed. The mixed leukocyte reaction (MLR) was used to determine T cell stimulating capacity. RESULTS: Serum-free culture of AML-APCs resulted in comparable morphology, relative yields and immunophenotype to serum-enriched culture. By comparing both serum-free media we observed a trend towards a more mature phenotype of CI-cultured AML-APCs in SP. MLR showed that serum-free cultured cells have equal T cell stimulatory capacity in comparison with serum-enriched culture. CONCLUSION: These data show that the serum-free culture of AML-APCs is feasible and that these APCs are comparable to serum-enriched cultured AML-APCs with regard to morphological, immunophenotypical and functional characteristics. These AML-APCs are suitable for the development of active specific immunisation protocols which meet the criteria for good clinical practise (GCP).     

Protein-protein docking predictions for the CAPRI experiment

We predicted structures for all seven targets in the CAPRI experiment using a new method in development at the time of the challenge. The technique includes a low-resolution rigid body Monte Carlo search followed by high-resolution refinement with side-chain conformational changes and rigid body minimization. Decoys (approximately 10(6) per target) were discriminated using a scoring function including van der Waals and solvation interactions, hydrogen bonding, residue-residue pair statistics, and rotamer probabilities. Decoys were ranked, clustered, manually inspected, and selected. The top ranked model for target 6 predicted the experimental structure to 1.5 A RMSD and included 48 of 65 correct residue-residue contacts. Target 7 was predicted at 5.3 A RMSD with 22 of 37 correct residue-residue contacts using a homology model from a known complex structure. Using a preliminary version of the protocol in round 1, target 1 was predicted within 8.8 A although few contacts were correct. For targets 2 and 3, the interface locations and a small fraction of the contacts were correctly identified.     

Crystal structure of the avilamycin resistance-conferring methyltransferase AviRa from Streptomyces viridochromogenes

The emergence of antibiotic-resistant bacterial strains is a widespread problem in contemporary medical practice and drug design. It is therefore important to elucidate the underlying mechanism in each case. The methyltransferase AviRa from Streptomyces viridochromogenes mediates resistance to the antibiotic avilamycin, which is closely related to evernimicin, an oligosaccharide antibiotic that has been used in medical studies. The structure of AviRa was determined by X-ray diffraction at 1.5A resolution. Phases were obtained from one selenomethionine residue introduced by site-directed mutagenesis. The chain-fold is similar to that of most methyltransferases, although AviRa contains two additional helices as a specific feature. A putative-binding site for the cofactor S-adenosyl-L-methionine was derived from homologous structures. It agrees with the conserved pattern of interacting amino acid residues. AviRa methylates a specific guanine base within the peptidyltransferase loop of the 23S ribosomal RNA. Guided by the target, the enzyme was docked to the cognate ribosomal surface, where it fit well into a deep cleft without contacting any ribosomal protein. The two additional alpha-helices of AviRa filled a depression in the surface. Since the transferred methyl group of the cofactor is in a pocket beneath the enzyme surface, the targeted guanine base has to flip out for methylation.