MetaMicrobesOnline: phylogenomic analysis of microbial communities.

The metaMicrobesOnline database (freely available at http://meta.MicrobesOnline.org) offers phylogenetic analysis of genes from microbial genomes and metagenomes. Gene trees are constructed for canonical gene families such as COG and Pfam. Such gene trees allow for rapid homologue analysis and subfamily comparison of genes from multiple metagenomes and comparisons with genes from microbial isolates. Additionally, the genome browser permits genome context comparisons, which may be used to determine the closest sequenced genome or suggest functionally associated genes. Lastly, the domain browser permits rapid comparison of protein domain organization within genes of interest from metagenomes and complete microbial genomes.

Tracing determinants of dual substrate specificity in glycoside hydrolase family 5.

Enzymes are traditionally viewed as having exquisite substrate specificity; however, recent evidence supports the notion that many enzymes have evolved activities against a range of substrates. The diversity of activities across glycoside hydrolase family 5 (GH5) suggests that this family of enzymes may contain numerous members with activities on multiple substrates. In this study, we combined structure- and sequence-based phylogenetic analysis with biochemical characterization to survey the prevalence of dual specificity for glucan- and mannan-based substrates in the GH5 family. Examination of amino acid profile differences between the subfamilies led to the identification and subsequent experimental confirmation of an active site motif indicative of dual specificity. The motif enabled us to successfully discover several new dually specific members of GH5, and this pattern is present in over 70 other enzymes, strongly suggesting that dual endoglucanase-mannanase activity is widespread in this family. In addition, reinstatement of the conserved motif in a wild type member of GH5 enhanced its catalytic efficiency on glucan and mannan substrates by 175 and 1,600%, respectively. Phylogenetic examination of other GH families further indicates that the prevalence of enzyme multispecificity in GHs may be greater than has been experimentally characterized. Single domain multispecific GHs may be exploited for developing improved enzyme cocktails or facile engineering of microbial hosts for consolidated bioprocessing of lignocellulose.

GLAMM: Genome-Linked Application for Metabolic Maps.

The Genome-Linked Application for Metabolic Maps (GLAMM) is a unified web interface for visualizing metabolic networks, reconstructing metabolic networks from annotated genome data, visualizing experimental data in the context of metabolic networks and investigating the construction of novel, transgenic pathways. This simple, user-friendly interface is tightly integrated with the comparative genomics tools of MicrobesOnline [Dehal et al. (2010) Nucleic Acids Research, 38, D396-D400]. GLAMM is available for free to the scientific community at glamm.lbl.gov.

Metagenome-assembled genome extraction and analysis from microbiomes using KBase.

Uncultivated Bacteria and Archaea account for the vast majority of species on Earth, but obtaining their genomes directly from the environment, using shotgun sequencing, has only become possible recently. To realize the hope of capturing Earth's microbial genetic complement and to facilitate the investigation of the functional roles of specific lineages in a given ecosystem, technologies that accelerate the recovery of high-quality genomes are necessary. We present a series of analysis steps and data products for the extraction of high-quality metagenome-assembled genomes (MAGs) from microbiomes using the U.S. Department of Energy Systems Biology Knowledgebase (KBase) platform ( http://www.kbase.us/ ). Overall, these steps take about a day to obtain extracted genomes when starting from smaller environmental shotgun read libraries, or up to about a week from larger libraries. In KBase, the process is end-to-end, allowing a user to go from the initial sequencing reads all the way through to MAGs, which can then be analyzed with other KBase capabilities such as phylogenetic placement, functional assignment, metabolic modeling, pangenome functional profiling, RNA-Seq and others. While portions of such capabilities are available individually from other resources, the combination of the intuitive usability, data interoperability and integration of tools in a freely available computational resource makes KBase a powerful platform for obtaining MAGs from microbiomes. While this workflow offers tools for each of the key steps in the genome extraction process, it also provides a scaffold that can be easily extended with additional MAG recovery and analysis tools, via the KBase software development kit (SDK).

Manipulation of the carbon storage regulator system for metabolite remodeling and biofuel production in Escherichia coli.

BACKGROUND: Microbial engineering strategies that elicit global metabolic perturbations have the capacity to increase organism robustness for targeted metabolite production. In particular, perturbations to regulators of cellular systems that impact glycolysis and amino acid production while simultaneously decreasing fermentation by-products such as acetate and CO(2) make ideal targets. Intriguingly, perturbation of the Carbon Storage Regulator (Csr) system has been previously implicated in large changes in central carbon metabolism in E. coli. Therefore, we hypothesized that perturbation of the Csr system through the CsrA-CsrB ribonucleoprotein complex might increase production of biofuels and their intermediates from heterologous pathways. RESULTS: We engaged the CsrA-CsrB ribonucleoprotein complex of E. coli via overexpression of CsrB. CsrB is a 350-nucleotide non-coding RNA that antagonizes CsrA, an RNA-binding protein that regulates translation of specific mRNA targets. By using shotgun proteomics and targeted metabolomics we established that elevation of CsrB levels leads to alterations in metabolite and protein levels in glycolysis, the TCA cycle and amino acid levels. Consequently, we show that such changes can be suitably applied to improve the production of desired compounds through the native fatty acid and heterologous n-butanol and isoprenoid pathways by up to two-fold. We also observed concomitant decreases in undesirable fermentation by-products such as acetate and CO(2). CONCLUSIONS: We have demonstrated that simple engineering of the RNA-based Csr global regulatory system constitutes a novel approach to obtaining pathway-independent improvements within engineered hosts. Additionally, since Csr is conserved across most prokaryotic species, this approach may also be amenable to a wide variety of production hosts.

MicrobesOnline: an integrated portal for comparative and functional genomics.

Since 2003, MicrobesOnline (http://www.microbesonline.org) has been providing a community resource for comparative and functional genome analysis. The portal includes over 1000 complete genomes of bacteria, archaea and fungi and thousands of expression microarrays from diverse organisms ranging from model organisms such as Escherichia coli and Saccharomyces cerevisiae to environmental microbes such as Desulfovibrio vulgaris and Shewanella oneidensis. To assist in annotating genes and in reconstructing their evolutionary history, MicrobesOnline includes a comparative genome browser based on phylogenetic trees for every gene family as well as a species tree. To identify co-regulated genes, MicrobesOnline can search for genes based on their expression profile, and provides tools for identifying regulatory motifs and seeing if they are conserved. MicrobesOnline also includes fast phylogenetic profile searches, comparative views of metabolic pathways, operon predictions, a workbench for sequence analysis and integration with RegTransBase and other microbial genome resources. The next update of MicrobesOnline will contain significant new functionality, including comparative analysis of metagenomic sequence data. Programmatic access to the database, along with source code and documentation, is available at http://microbesonline.org/programmers.html.

Optimization of a heterologous mevalonate pathway through the use of variant HMG-CoA reductases.

Expression of foreign pathways often results in suboptimal performance due to unintended factors such as introduction of toxic metabolites, cofactor imbalances or poor expression of pathway components. In this study we report a 120% improvement in the production of the isoprenoid-derived sesquiterpene, amorphadiene, produced by an engineered strain of Escherichia coli developed to express the native seven-gene mevalonate pathway from Saccharomyces cerevisiae (Martin et al. 2003). This substantial improvement was made by varying only a single component of the pathway (HMG-CoA reductase) and subsequent host optimization to improve cofactor availability. We characterized and tested five variant HMG-CoA reductases obtained from publicly available genome databases with differing kinetic properties and cofactor requirements. The results of our in vitro and in vivo analyses of these enzymes implicate substrate inhibition of mevalonate kinase as an important factor in optimization of the engineered mevalonate pathway. Consequently, the NADH-dependent HMG-CoA reductase from Delftia acidovorans, which appeared to have the optimal kinetic parameters to balance HMG-CoA levels below the cellular toxicity threshold of E. coli and those of mevalonate below inhibitory concentrations for mevalonate kinase, was identified as the best producer for amorphadiene (54% improvement over the native pathway enzyme, resulting in 2.5mM or 520 mg/L of amorphadiene after 48 h). We further enhanced performance of the strain bearing the D. acidovorans HMG-CoA reductase by increasing the intracellular levels of its preferred cofactor (NADH) using a NAD(+)-dependent formate dehydrogenase from Candida boidinii, along with formate supplementation. This resulted in an overall improvement of the system by 120% resulting in 3.5mM or 700 mg/L amorphadiene after 48 h of fermentation. This comprehensive study incorporated analysis of several key parameters for metabolic design such as in vitro and in vivo kinetic performance of variant enzymes, intracellular levels of protein expression, in-pathway substrate inhibition and cofactor management to enable the observed improvements. These metrics may be applied to a broad range of heterologous pathways for improving the production of biologically derived compounds.

A genomic catalog of Earth's microbiomes.

The reconstruction of bacterial and archaeal genomes from shotgun metagenomes has enabled insights into the ecology and evolution of environmental and host-associated microbiomes. Here we applied this approach to >10,000 metagenomes collected from diverse habitats covering all of Earth's continents and oceans, including metagenomes from human and animal hosts, engineered environments, and natural and agricultural soils, to capture extant microbial, metabolic and functional potential. This comprehensive catalog includes 52,515 metagenome-assembled genomes representing 12,556 novel candidate species-level operational taxonomic units spanning 135 phyla. The catalog expands the known phylogenetic diversity of bacteria and archaea by 44% and is broadly available for streamlined comparative analyses, interactive exploration, metabolic modeling and bulk download. We demonstrate the utility of this collection for understanding secondary-metabolite biosynthetic potential and for resolving thousands of new host linkages to uncultivated viruses. This resource underscores the value of genome-centric approaches for revealing genomic properties of uncultivated microorganisms that affect ecosystem processes.

Superfamily assignments for the yeast proteome through integration of structure prediction with the gene ontology.

Saccharomyces cerevisiae is one of the best-studied model organisms, yet the three-dimensional structure and molecular function of many yeast proteins remain unknown. Yeast proteins were parsed into 14,934 domains, and those lacking sequence similarity to proteins of known structure were folded using the Rosetta de novo structure prediction method on the World Community Grid. This structural data was integrated with process, component, and function annotations from the Saccharomyces Genome Database to assign yeast protein domains to SCOP superfamilies using a simple Bayesian approach. We have predicted the structure of 3,338 putative domains and assigned SCOP superfamily annotations to 581 of them. We have also assigned structural annotations to 7,094 predicted domains based on fold recognition and homology modeling methods. The domain predictions and structural information are available in an online database at http://rd.plos.org/10.1371_journal.pbio.0050076_01.

Contribution of mobile genetic elements to Desulfovibrio vulgaris genome plasticity.

The genome of Desulfovibrio vulgaris strain DePue, a sulfate-reducing Deltaproteobacterium isolated from heavy metal-impacted lake sediment, was completely sequenced and compared with the type strain D. vulgaris Hildenborough. The two genomes share a high degree of relatedness and synteny, but harbour distinct prophage and signatures of past phage encounters. In addition to a highly variable phage contribution, the genome of strain DePue contains a cluster of open-reading frames not found in strain Hildenborough coding for the production and export of a capsule exopolysaccharide, possibly of relevance to heavy metal resistance. Comparative whole-genome microarray analysis on four additional D. vulgaris strains established greater interstrain variation within regions associated with phage insertion and exopolysaccharide biosynthesis.

Homology modeling using parametric alignment ensemble generation with consensus and energy-based model selection.

The accuracy of a homology model based on the structure of a distant relative or other topologically equivalent protein is primarily limited by the quality of the alignment. Here we describe a systematic approach for sequence-to-structure alignment, called 'K*Sync', in which alignments are generated by dynamic programming using a scoring function that combines information on many protein features, including a novel measure of how obligate a sequence region is to the protein fold. By systematically varying the weights on the different features that contribute to the alignment score, we generate very large ensembles of diverse alignments, each optimal under a particular constellation of weights. We investigate a variety of approaches to select the best models from the ensemble, including consensus of the alignments, a hydrophobic burial measure, low- and high-resolution energy functions, and combinations of these evaluation methods. The effect on model quality and selection resulting from loop modeling and backbone optimization is also studied. The performance of the method on a benchmark set is reported and shows the approach to be effective at both generating and selecting accurate alignments. The method serves as the foundation of the homology modeling module in the Robetta server.

Assessment of predictions submitted for the CASP7 domain prediction category.

This paper details the assessment process and evaluation results for the Critical Assessment of Protein Structure Prediction (CASP7) domain prediction category. Domain predictions were assessed using the Normalized Domain Overlap score introduced in CASP6 and the accuracy of prediction of domain break points. The results of the analysis clearly demonstrate that the best methods are able to make consistently reliable predictions when the target has a structural template, although they are less good when the domain break occurs in a region not covered by a template. The conditions of the experiment meant that it was impossible to draw any conclusions about domain prediction for free modeling targets and it was also difficult to draw many distinctions between the best groups. Two thirds of the targets submitted were single domains and hence regarded as easy to predict. Even those targets defined as having multiple domains always had at least one domain with a similar template structure.

Structure prediction for CASP7 targets using extensive all-atom refinement with Rosetta@home.

We describe predictions made using the Rosetta structure prediction methodology for both template-based modeling and free modeling categories in the Seventh Critical Assessment of Techniques for Protein Structure Prediction. For the first time, aggressive sampling and all-atom refinement could be carried out for the majority of targets, an advance enabled by the Rosetta@home distributed computing network. Template-based modeling predictions using an iterative refinement algorithm improved over the best existing templates for the majority of proteins with less than 200 residues. Free modeling methods gave near-atomic accuracy predictions for several targets under 100 residues from all secondary structure classes. These results indicate that refinement with an all-atom energy function, although computationally expensive, is a powerful method for obtaining accurate structure predictions.

Genomics for environmental microbiology.

The utilization of natural microbial diversity in biotechnology is hindered by our inability to culture the vast majority of microorganisms and the observation that laboratory engineered bacteria rarely function in the wild. It is now clear that an understanding of the community structure, function and evolution of bacteria in their natural environments is required to meet the promise of microbial biotechnology. To meet these new challenges, microbiologists are applying the tools of genomics and related high-throughput technologies to both cultured microbes and environmental samples. This work will lead to new views on ecosystems and biological function together with the biotechnology enabled by this science.

Protein structure prediction and analysis using the Robetta server.

The Robetta server (http://robetta.bakerlab.org) provides automated tools for protein structure prediction and analysis. For structure prediction, sequences submitted to the server are parsed into putative domains and structural models are generated using either comparative modeling or de novo structure prediction methods. If a confident match to a protein of known structure is found using BLAST, PSI-BLAST, FFAS03 or 3D-Jury, it is used as a template for comparative modeling. If no match is found, structure predictions are made using the de novo Rosetta fragment insertion method. Experimental nuclear magnetic resonance (NMR) constraints data can also be submitted with a query sequence for RosettaNMR de novo structure determination. Other current capabilities include the prediction of the effects of mutations on protein-protein interactions using computational interface alanine scanning. The Rosetta protein design and protein-protein docking methodologies will soon be available through the server as well.

Automated prediction of domain boundaries in CASP6 targets using Ginzu and RosettaDOM.

Domain boundary prediction is an important step in both experimental and computational protein structure characterization. We have developed two fully automated domain parsing methods: the first, Ginzu, which we have described previously, utilizes information from homologous sequences and structures, while the second, RosettaDOM, which has not been described previously, uses only information in the query sequence. Ginzu iteratively assigns domains by homology to structures and sequence families using successively less confident methods. RosettaDOM uses the Rosetta de novo structure prediction method to build three-dimensional models, and then applies Taylor's structure based domain assignment method to parse the models into domains. Domain boundaries observed repeatedly in the models are predicted to be domain boundaries for the protein. Interestingly, RosettaDOM produced quite good domain predictions for proteins of a size typically considered to be beyond the reach of de novo structure prediction methods. For remote fold recognition targets and new folds, both Ginzu and RosettaDOM produced promising results, and in some cases where one method failed to detect the correct domain boundary, it was correctly identified by the other method. We describe here the successes and failures using both methods, and address the possibility of incorporating both protocols into an improved hybrid method.

Prediction of CASP6 structures using automated Robetta protocols.

The Robetta server and revised automatic protocols were used to predict structures for CASP6 targets. Robetta is a publicly available protein structure prediction server (http://robetta.bakerlab.org/ that uses the Rosetta de novo and homology modeling structure prediction methods. We incorporated some of the lessons learned in the CASP5 experiment into the server prior to participating in CASP6. We additionally tested new ideas that were amenable to full-automation with an eye toward improving the server. We find that the Robetta server shows the greatest promise for the more challenging targets. The most significant finding from CASP5, that automated protocols can be roughly comparable in ability with the better human-intervention predictors, is repeated here in CASP6.