Andrew C. Doxey

Ordered surface carbons distinguish antifreeze proteins and their ice-binding regions

Antifreeze proteins (AFPs) are found in cold-adapted organisms and have the unusual ability to bind to and inhibit the growth of ice crystals. However, the underlying molecular basis of their ice-binding activity is unclear because of the difficulty of studying the AFP-ice interaction directly and the lack of a common motif, domain or fold among different AFPs. We have formulated a generic ice-binding model and incorporated it into a physicochemical pattern-recognition algorithm. It successfully recognizes ice-binding surfaces for a diverse range of AFPs, and clearly discriminates AFPs from other structures in the Protein Data Bank. The algorithm was used to identify a novel AFP from winter rye, and the antifreeze activity of this protein was subsequently confirmed. The presence of a common and distinct physicochemical pattern provides a structural basis for unifying AFPs from fish, insects and plants.

Modular Evolution and the Origins of Symmetry: Reconstruction of a Three-Fold Symmetric Globular Protein

The high frequency of internal structural symmetry in common protein folds is presumed to reflect their evolutionary origins from the repetition and fusion of ancient peptide modules, but little is known about the primary sequence and physical determinants of this process. Unexpectedly, a sequence and structural analysis of symmetric subdomain modules within an abundant and ancient globular fold, the β-trefoil, reveals that modular evolution is not simply a relic of the ancient past, but is an ongoing and recurring mechanism for regenerating symmetry, having occurred independently in numerous existing β-trefoil proteins. We performed a computational reconstruction of a β-trefoil subdomain module and repeated it to form a newly three-fold symmetric globular protein, ThreeFoil. In addition to its near perfect structural identity between symmetric modules, ThreeFoil is highly soluble, performs multivalent carbohydrate binding, and has remarkably high thermal stability. These findings have far-reaching implications for understanding the evolution and design of proteins via subdomain modules.

Cold-Active Winter Rye Glucanases with Ice-Binding Capacity

Extracellular pathogenesis-related proteins, including glucanases, are expressed at cold temperatures in winter rye (Secale cereale) and display antifreeze activity. We have characterized recombinant cold-induced glucanases from winter rye to further examine their roles and contributions to cold tolerance. Both basic β-1,3-glucanases and an acidic β-1,3;1,4-glucanase were expressed in Escherichia coli, purified, and assayed for their hydrolytic and antifreeze activities in vitro. All were found to be cold active and to retain partial hydrolytic activity at subzero temperatures (e.g. 14%–35% at −4°C). The two types of glucanases had antifreeze activity as measured by their ability to modify the growth of ice crystals. Structural models for the winter rye β-1,3-glucanases were developed on which putative ice-binding surfaces (IBSs) were identified. Residues on the putative IBSs were charge conserved for each of the expressed glucanases, with the exception of one β-1,3-glucanase recovered from nonacclimated winter rye in which a charged amino acid was present on the putative IBS. This protein also had a reduced antifreeze activity relative to the other expressed glucanases. These results support the hypothesis that winter rye glucanases have evolved to inhibit the formation of large, potentially fatal ice crystals, in addition to having enzymatic activity with a potential role in resisting infection by psychrophilic pathogens. Glucanases of winter rye provide an interesting example of protein evolution and adaptation aimed to combat cold and freezing conditions.

Aquatic metagenomes implicate Thaumarchaeota in global cobalamin production.

Cobalamin (vitamin B12) is a complex metabolite and essential cofactor required by many branches of life, including most eukaryotic phytoplankton. Algae and other cobalamin auxotrophs rely on environmental cobalamin supplied from a relatively small set of cobalamin-producing prokaryotic taxa. Although several Bacteria have been implicated in cobalamin biosynthesis and associated with algal symbiosis, the involvement of Archaea in cobalamin production is poorly understood, especially with respect to the Thaumarchaeota. Based on the detection of cobalamin synthesis genes in available thaumarchaeotal genomes, we hypothesized that Thaumarchaeota, which are ubiquitous and abundant in aquatic environments, have an important role in cobalamin biosynthesis within global aquatic ecosystems. To test this hypothesis, we examined cobalamin synthesis genes across sequenced thaumarchaeotal genomes and 430 metagenomes from a diverse range of marine, freshwater and hypersaline environments. Our analysis demonstrates that all available thaumarchaeotal genomes possess cobalamin synthesis genes, predominantly from the anaerobic pathway, suggesting widespread genetic capacity for cobalamin synthesis. Furthermore, although bacterial cobalamin genes dominated most surface marine metagenomes, thaumarchaeotal cobalamin genes dominated metagenomes from polar marine environments, increased with depth in marine water columns, and displayed seasonality, with increased winter abundance observed in time-series datasets (e.g., L4 surface water in the English Channel). Our results also suggest niche partitioning between thaumarchaeotal and cyanobacterial ribosomal and cobalamin synthesis genes across all metagenomic datasets analyzed. These results provide strong evidence for specific biogeographical distributions of thaumarchaeotal cobalamin genes, expanding our understanding of the global biogeochemical roles played by Thaumarchaeota in aquatic environments.

Botulinum neurotoxin homologs in non-Clostridium species.

Clostridial neurotoxins (CNTs) are the deadliest toxins known and the causative agents of botulism and tetanus. Despite their structural and functional complexity, no CNT homologs are currently known outside Clostridium. Here, we report the first homologs of Clostridium CNTs within the genome of the rice fermentation organism Weissella oryzae SG25. One gene in W. oryzae S25 encodes a protein with a four‐domain architecture and HExxH protease motif common to botulinum neurotoxins (BoNTs). An adjacent gene with partial similarity to CNTs is also present, and both genes seem to have been laterally transferred into the W. oryzae genome from an unknown source. Identification of mobile, CNT‐related genes outside of Clostridium has implications for our understanding of the evolution of this important toxin family.

Prediction of molecular mimicry candidates in human pathogenic bacteria

Molecular mimicry of host proteins is a common strategy adopted by bacterial pathogens to interfere with and exploit host processes. Despite the availability of pathogen genomes, few studies have attempted to predict virulence-associated mimicry relationships directly from genomic sequences. Here, we analyzed the proteomes of 62 pathogenic and 66 non-pathogenic bacterial species, and screened for the top pathogen-specific or pathogen-enriched sequence similarities to human proteins. The screen identified approximately 100 potential mimicry relationships including well-characterized examples among the top-scoring hits (e.g., RalF, internalin, yopH, and others), with about 1/3 of predicted relationships supported by existing literature. Examination of homology to virulence factors, statistically enriched functions, and comparison with literature indicated that the detected mimics target key host structures (e.g., extracellular matrix, ECM) and pathways (e.g., cell adhesion, lipid metabolism, and immune signaling). The top-scoring and most widespread mimicry pattern detected among pathogens consisted of elevated sequence similarities to ECM proteins including collagens and leucine-rich repeat proteins. Unexpectedly, analysis of the pathogen counterparts of these proteins revealed that they have evolved independently in different species of bacterial pathogens from separate repeat amplifications. Thus, our analysis provides evidence for two classes of mimics: complex proteins such as enzymes that have been acquired by eukaryote-to-pathogen horizontal transfer, and simpler repeat proteins that have independently evolved to mimic the host ECM. Ultimately, computational detection of pathogen-specific and pathogen-enriched similarities to host proteins provides insights into potentially novel mimicry-mediated virulence mechanisms of pathogenic bacteria.

Structure-aided prediction of mammalian transcription factor complexes in conserved non-coding elements

Mapping the DNA-binding preferences of transcription factor (TF) complexes is critical for deciphering the functions of cis-regulatory elements. Here, we developed a computational method that compares co-occurring motif spacings in conserved versus unconserved regions of the human genome to detect evolutionarily constrained binding sites of rigid TF complexes. Structural data were used to estimate TF complex physical plausibility, explore overlapping motif arrangements seldom tackled by non-structure-aware methods, and generate and analyse three-dimensional models of the predicted complexes bound to DNA. Using this approach, we predicted 422 physically realistic TF complex motifs at 18% false discovery rate, the majority of which (326, 77%) contain some sequence overlap between binding sites. The set of mostly novel complexes is enriched in known composite motifs, predictive of binding site configurations in TF–TF–DNA crystal structures, and supported by ChIP-seq datasets. Structural modelling revealed three cooperativity mechanisms: direct protein–protein interactions, potentially indirect interactions and ‘through-DNA’ interactions. Indeed, 38% of the predicted complexes were found to contain four or more bases in which TF pairs appear to synergize through overlapping binding to the same DNA base pairs in opposite grooves or strands. Our TF complex and associated binding site predictions are available as a web resource at http://bejerano.stanford.edu/complex.

Horizontal Gene Transfer Contributed to the Evolution of Extracellular Surface Structures: The Freshwater Polyp Hydra Is Covered by a Complex Fibrous Cuticle Containing Glycosaminoglycans and Proteins of the PPOD and SWT (Sweet Tooth) Families

The single-cell layered ectoderm of the fresh water polyp Hydra fulfills the function of an epidermis by protecting the animals from the surrounding medium. Its outer surface is covered by a fibrous structure termed the cuticle layer, with similarity to the extracellular surface coats of mammalian epithelia. In this paper we have identified molecular components of the cuticle. We show that its outermost layer contains glycoproteins and glycosaminoglycans and we have identified chondroitin and chondroitin-6-sulfate chains. In a search for proteins that could be involved in organising this structure we found PPOD proteins and several members of a protein family containing only SWT (sweet tooth) domains. Structural analyses indicate that PPODs consist of two tandem β-trefoil domains with similarity to carbohydrate-binding sites found in lectins. Experimental evidence confirmed that PPODs can bind sulfated glycans and are secreted into the cuticle layer from granules localized under the apical surface of the ectodermal epithelial cells. PPODs are taxon-specific proteins which appear to have entered the Hydra genome by horizontal gene transfer from bacteria. Their acquisition at the time Hydra evolved from a marine ancestor may have been critical for the transition to the freshwater environment.

Insights into the evolutionary origins of clostridial neurotoxins from analysis of the Clostridium botulinum strain A neurotoxin gene cluster

BackgroundClostridial neurotoxins (CNTs) are the most deadly toxins known and causal agents of botulism and tetanus neuroparalytic diseases. Despite considerable progress in understanding CNT structure and function, the evolutionary origins of CNTs remain a mystery as they are unique to Clostridium and possess a sequence and structural architecture distinct from other protein families. Uncovering the origins of CNTs would be a significant contribution to our understanding of how pathogens evolve and generate novel toxin families.ResultsThe C. botulinum strain A genome was examined for potential homologues of CNTs. A key link was identified between the neurotoxin and the flagellin gene (CBO0798) located immediately upstream of the BoNT/A neurotoxin gene cluster. This flagellin sequence displayed the strongest sequence similarity to the neurotoxin and NTNH homologue out of all proteins encoded within C. botulinum strain A. The CBO0798 gene contains a unique hypervariable region, which in closely related flagellins encodes a collagenase-like domain. Remarkably, these collagenase-containing flagellins were found to possess the characteristic HEXXH zinc-protease motif responsible for the neurotoxins endopeptidase activity. Additional links to collagenase-related sequences and functions were detected by further analysis of CNTs and surrounding genes, including sequence similarities to collagen-adhesion domains and collagenases. Furthermore, the neurotoxins HCRn domain was found to exhibit both structural and sequence similarity to eukaryotic collagen jelly-roll domains.ConclusionMultiple lines of evidence suggest that the neurotoxin and adjacent genes evolved from an ancestral collagenase-like gene cluster, linking CNTs to another major family of clostridial proteolytic toxins. Duplication, reshuffling and assembly of neighboring genes within the BoNT/A neurotoxin gene cluster may have lead to the neurotoxins unique architecture. This work provides new insights into the evolution of C. botulinum neurotoxins and the evolutionary mechanisms underlying the origins of virulent genes.

MetAnnotate: function-specific taxonomic profiling and comparison of metagenomes.

BackgroundMetagenomes provide access to the taxonomic composition and functional capabilities of microbial communities. Although metagenomic analysis methods exist for estimating overall community composition or metabolic potential, identifying specific taxa that encode specific functions or pathways of interest can be more challenging. Here we present MetAnnotate, which addresses the common question: “which organisms perform my function of interest within my metagenome(s) of interest?” MetAnnotate uses profile hidden Markov models to analyze shotgun metagenomes for genes and pathways of interest, classifies retrieved sequences either through a phylogenetic placement or best hit approach, and enables comparison of these profiles between metagenomes.ResultsBased on a simulated metagenome dataset, the tool achieves high taxonomic classification accuracy for a broad range of genes, including both markers of community abundance and specific biological pathways. Lastly, we demonstrate MetAnnotate by analyzing for cobalamin (vitamin B12) synthesis genes across hundreds of aquatic metagenomes in a fraction of the time required by the commonly used Basic Local Alignment Search Tool top hit approach.ConclusionsMetAnnotate is multi-threaded and installable as a local web application or command-line tool on Linux systems. Metannotate is a useful framework for general and/or function-specific taxonomic profiling and comparison of metagenomes.