A Score of the Ability of a Three-Dimensional Protein Model to Retrieve Its Own Sequence as a Quantitative Measure of Its Quality and Appropriateness

BACKGROUND: Despite the remarkable progress of bioinformatics, how the primary structure of a protein leads to a three-dimensional fold, and in turn determines its function remains an elusive question. Alignments of sequences with known function can be used to identify proteins with the same or similar function with high success. However, identification of function-related and structure-related amino acid positions is only possible after a detailed study of every protein. Folding pattern diversity seems to be much narrower than sequence diversity, and the amino acid sequences of natural proteins have evolved under a selective pressure comprising structural and functional requirements acting in parallel. PRINCIPAL FINDINGS: The approach described in this work begins by generating a large number of amino acid sequences using ROSETTA [Dantas G et al. (2003) J Mol Biol 332:449-460], a program with notable robustness in the assignment of amino acids to a known three-dimensional structure. The resulting sequence-sets showed no conservation of amino acids at active sites, or protein-protein interfaces. Hidden Markov models built from the resulting sequence sets were used to search sequence databases. Surprisingly, the models retrieved from the database sequences belonged to proteins with the same or a very similar function. Given an appropriate cutoff, the rate of false positives was zero. According to our results, this protocol, here referred to as Rd.HMM, detects fine structural details on the folding patterns, that seem to be tightly linked to the fitness of a structural framework for a specific biological function. CONCLUSION: Because the sequence of the native protein used to create the Rd.HMM model was always amongst the top hits, the procedure is a reliable tool to score, very accurately, the quality and appropriateness of computer-modeled 3D-structures, without the need for spectroscopy data. However, Rd.HMM is very sensitive to the conformational features of the models' backbone

MediPlEx - a tool to combine in silico & experimental gene expression profiles of the model legume Medicago truncatula

Henckel K, Küster H, Stutz L, Goesmann A. MediPlEx - a tool to combine in silico and experimental gene expression profiles of the model legume Medicago truncatula. BMC Research Notes. 2010;3(1): 262.BACKGROUND:Expressed Sequence Tags (ESTs) are in general used to gain a first insight into gene activities from a species of interest. Subsequently, and typically based on a combination of EST and genome sequences, microarray-based expression analyses are performed for a variety of conditions. In some cases, a multitude of EST and microarray experiments are conducted for one species, covering different tissues, cell states, and cell types. Under these circumstances, the challenge arises to combine results derived from the different expression profiling strategies, with the goal to uncover novel information on the basis of the integrated datasets.FINDINGS:Using our new application, MediPlEx (MEDIcago truncatula multiPLe EXpression analysis), expression data from EST experiments, oligonucleotide microarrays and Affymetrix GeneChips can be combined and analyzed, leading to a novel approach to integrated transcriptome analysis. We have validated our tool via the identification of a set of well-characterized AM-specific and AM-induced marker genes, identified by MediPlEx on the basis of in silico and experimental gene expression profiles from roots colonized with AM fungi.CONCLUSIONS:MediPlEx offers an integrated analysis pipeline for different sets of expression data generated for the model legume Medicago truncatula. As expected, in silico and experimental gene expression data that cover the same biological condition correlate well. The collection of differentially expressed genes identified via MediPlEx provides a starting point for functional studies in plant mutants. MediPlEx can freely be used at http://www.cebitec.uni-bielefeld.de/mediplex

Improved accuracy of multiple ncRNA alignment by incorporating structural information into a MAFFT-based framework

<p>Abstract</p> <p>Background</p> <p>Structural alignment of RNAs is becoming important, since the discovery of functional non-coding RNAs (ncRNAs). Recent studies, mainly based on various approximations of the Sankoff algorithm, have resulted in considerable improvement in the accuracy of pairwise structural alignment. In contrast, for the cases with more than two sequences, the practical merit of structural alignment remains unclear as compared to traditional sequence-based methods, although the importance of multiple structural alignment is widely recognized.</p> <p>Results</p> <p>We took a different approach from a straightforward extension of the Sankoff algorithm to the multiple alignments from the viewpoints of accuracy and time complexity. As a new option of the MAFFT alignment program, we developed a multiple RNA alignment framework, X-INS-i, which builds a multiple alignment with an iterative method incorporating structural information through two components: (1) pairwise structural alignments by an external pairwise alignment method such as SCARNA or LaRA and (2) a new objective function, Four-way Consistency, derived from the base-pairing probability of every sub-aligned group at every multiple alignment stage.</p> <p>Conclusion</p> <p>The BRAliBASE benchmark showed that X-INS-i outperforms other methods currently available in the sum-of-pairs score (SPS) criterion. As a basis for predicting common secondary structure, the accuracy of the present method is comparable to or rather higher than those of the current leading methods such as RNA Sampler. The X-INS-i framework can be used for building a multiple RNA alignment from any combination of algorithms for pairwise RNA alignment and base-pairing probability. The source code is available at the webpage found in the Availability and requirements section.</p

A Mathematical Framework for Protein Structure Comparison

Comparison of protein structures is important for revealing the evolutionary relationship among proteins, predicting protein functions and predicting protein structures. Many methods have been developed in the past to align two or multiple protein structures. Despite the importance of this problem, rigorous mathematical or statistical frameworks have seldom been pursued for general protein structure comparison. One notable issue in this field is that with many different distances used to measure the similarity between protein structures, none of them are proper distances when protein structures of different sequences are compared. Statistical approaches based on those non-proper distances or similarity scores as random variables are thus not mathematically rigorous. In this work, we develop a mathematical framework for protein structure comparison by treating protein structures as three-dimensional curves. Using an elastic Riemannian metric on spaces of curves, geodesic distance, a proper distance on spaces of curves, can be computed for any two protein structures. In this framework, protein structures can be treated as random variables on the shape manifold, and means and covariance can be computed for populations of protein structures. Furthermore, these moments can be used to build Gaussian-type probability distributions of protein structures for use in hypothesis testing. The covariance of a population of protein structures can reveal the population-specific variations and be helpful in improving structure classification. With curves representing protein structures, the matching is performed using elastic shape analysis of curves, which can effectively model conformational changes and insertions/deletions. We show that our method performs comparably with commonly used methods in protein structure classification on a large manually annotated data set

Heterotopic Ossifications in a Mouse Model of Albright Hereditary Osteodystrophy

Albright hereditary osteodystrophy (AHO) is characterized by short stature, brachydactyly, and often heterotopic ossifications that are typically subcutaneous. Subcutaneous ossifications (SCO) cause considerable morbidity in AHO with no effective treatment. AHO is caused by heterozygous inactivating mutations in those GNAS exons encoding the α-subunit of the stimulatory G protein (Gαs). When inherited maternally, these mutations are associated with obesity, cognitive impairment, and resistance to certain hormones that mediate their actions through G protein-coupled receptors, a condition termed pseudohypoparathyroidism type 1a (PHP1a). When inherited paternally, GNAS mutations cause only AHO but not hormonal resistance, termed pseudopseudohypoparathyroidism (PPHP). Mice with targeted disruption of exon 1 of Gnas (GnasE1−/+) replicate human PHP1a or PPHP phenotypically and hormonally. However, SCO have not yet been reported in GnasE1+/− mice, at least not those that had been analyzed by us up to 3 months of age. Here we now show that GnasE1−/+ animals develop SCO over time. The ossified lesions increase in number and size and are uniformly detected in adult mice by one year of age. They are located in both the dermis, often in perifollicular areas, and the subcutis. These lesions are particularly prominent in skin prone to injury or pressure. The SCO comprise mature bone with evidence of mineral deposition and bone marrow elements. Superficial localization was confirmed by radiographic and computerized tomographic imaging. In situ hybridization of SCO lesions were positive for both osteonectin and osteopontin. Notably, the ossifications were much more extensive in males than females. Because GnasE1−/+ mice develop SCO features that are similar to those observed in AHO patients, these animals provide a model system suitable for investigating pathogenic mechanisms involved in SCO formation and for developing novel therapeutics for heterotopic bone formation. Moreover, these mice provide a model with which to investigate the regulatory mechanisms of bone formation

The Essentials of Protein Import in the Degenerate Mitochondrion of Entamoeba histolytica

Several essential biochemical processes are situated in mitochondria. The metabolic transformation of mitochondria in distinct lineages of eukaryotes created proteomes ranging from thousands of proteins to what appear to be a much simpler scenario. In the case of Entamoeba histolytica, tiny mitochondria known as mitosomes have undergone extreme reduction. Only recently a single complete metabolic pathway of sulfate activation has been identified in these organelles. The E. histolytica mitosomes do not produce ATP needed for the sulfate activation pathway and for three molecular chaperones, Cpn60, Cpn10 and mtHsp70. The already characterized ADP/ATP carrier would thus be essential to provide cytosolic ATP for these processes, but how the equilibrium of inorganic phosphate could be maintained was unknown. Finally, how the mitosomal proteins are translocated to the mitosomes had remained unclear. We used a hidden Markov model (HMM) based search of the E. histolytica genome sequence to discover candidate (i) mitosomal phosphate carrier complementing the activity of the ADP/ATP carrier and (ii) membrane-located components of the protein import machinery that includes the outer membrane translocation channel Tom40 and membrane assembly protein Sam50. Using in vitro and in vivo systems we show that E. histolytica contains a minimalist set up of the core import components in order to accommodate a handful of mitosomal proteins. The anaerobic and parasitic lifestyle of E. histolytica has produced one of the simplest known mitochondrial compartments of all eukaryotes. Comparisons with mitochondria of another amoeba, Dictystelium discoideum, emphasize just how dramatic the reduction of the protein import apparatus was after the loss of archetypal mitochondrial functions in the mitosomes of E. histolytica

Extensive microbial and functional diversity within the chicken cecal microbiome

Chickens are major source of food and protein worldwide. Feed conversion and the health of chickens relies on the largely unexplored complex microbial community that inhabits the chicken gut, including the ceca. We have carried out deep microbial community profiling of the microbiota in twenty cecal samples via 16S rRNA gene sequences and an in-depth metagenomics analysis of a single cecal microbiota. We recovered 699 phylotypes, over half of which appear to represent previously unknown species. We obtained 648,251 environmental gene tags (EGTs), the majority of which represent new species. These were binned into over two-dozen draft genomes, which included Campylobacter jejuni and Helicobacter pullorum. We found numerous polysaccharide- and oligosaccharide-degrading enzymes encoding within the metagenome, some of which appeared to be part of polysaccharide utilization systems with genetic evidence for the co-ordination of polysaccharide degradation with sugar transport and utilization. The cecal metagenome encodes several fermentation pathways leading to the production of short-chain fatty acids, including some with novel features. We found a dozen uptake hydrogenases encoded in the metagenome and speculate that these provide major hydrogen sinks within this microbial community and might explain the high abundance of several genera within this microbiome, including Campylobacter, Helicobacter and Megamonas

Reconstruction of the Core and Extended Regulons of Global Transcription Factors

The processes underlying the evolution of regulatory networks are unclear. To address this question, we used a comparative genomics approach that takes advantage of the large number of sequenced bacterial genomes to predict conserved and variable members of transcriptional regulatory networks across phylogenetically related organisms. Specifically, we developed a computational method to predict the conserved regulons of transcription factors across α-proteobacteria. We focused on the CRP/FNR super-family of transcription factors because it contains several well-characterized members, such as FNR, FixK, and DNR. While FNR, FixK, and DNR are each proposed to regulate different aspects of anaerobic metabolism, they are predicted to recognize very similar DNA target sequences, and they occur in various combinations among individual α-proteobacterial species. In this study, the composition of the respective FNR, FixK, or DNR conserved regulons across 87 α-proteobacterial species was predicted by comparing the phylogenetic profiles of the regulators with the profiles of putative target genes. The utility of our predictions was evaluated by experimentally characterizing the FnrL regulon (a FNR-type regulator) in the α-proteobacterium Rhodobacter sphaeroides. Our results show that this approach correctly predicted many regulon members, provided new insights into the biological functions of the respective regulons for these regulators, and suggested models for the evolution of the corresponding transcriptional networks. Our findings also predict that, at least for the FNR-type regulators, there is a core set of target genes conserved across many species. In addition, the members of the so-called extended regulons for the FNR-type regulators vary even among closely related species, possibly reflecting species-specific adaptation to environmental and other factors. The comparative genomics approach we developed is readily applicable to other regulatory networks

Comparative Genomics of Gardnerella vaginalis Strains Reveals Substantial Differences in Metabolic and Virulence Potential

Gardnerella vaginalis is described as a common vaginal bacterial species whose presence correlates strongly with bacterial vaginosis (BV). Here we report the genome sequencing and comparative analyses of three strains of G. vaginalis. Strains 317 (ATCC 14019) and 594 (ATCC 14018) were isolated from the vaginal tracts of women with symptomatic BV, while Strain 409-05 was isolated from a healthy, asymptomatic individual with a Nugent score of 9.Substantial genomic rearrangement and heterogeneity were observed that appeared to have resulted from both mobile elements and substantial lateral gene transfer. These genomic differences translated to differences in metabolic potential. All strains are equipped with significant virulence potential, including genes encoding the previously described vaginolysin, pili for cytoadhesion, EPS biosynthetic genes for biofilm formation, and antimicrobial resistance systems, We also observed systems promoting multi-drug and lantibiotic extrusion. All G. vaginalis strains possess a large number of genes that may enhance their ability to compete with and exclude other vaginal colonists. These include up to six toxin-antitoxin systems and up to nine additional antitoxins lacking cognate toxins, several of which are clustered within each genome. All strains encode bacteriocidal toxins, including two lysozyme-like toxins produced uniquely by strain 409-05. Interestingly, the BV isolates encode numerous proteins not found in strain 409-05 that likely increase their pathogenic potential. These include enzymes enabling mucin degradation, a trait previously described to strongly correlate with BV, although commonly attributed to non-G. vaginalis species.Collectively, our results indicate that all three strains are able to thrive in vaginal environments, and therein the BV isolates are capable of occupying a niche that is unique from 409-05. Each strain has significant virulence potential, although genomic and metabolic differences, such as the ability to degrade mucin, indicate that the detection of G. vaginalis in the vaginal tract provides only partial information on the physiological potential of the organism