The binary protein interactome of Treponema pallidum--the syphilis spirochete

Protein interaction networks shed light on the global organization of proteomes but can also place individual proteins into a functional context. If we know the function of bacterial proteins we will be able to understand how these species have adapted to diverse environments including many extreme habitats. Here we present the protein interaction network for the syphilis spirochete Treponema pallidum which encodes 1,039 proteins, 726 (or 70%) of which interact via 3,649 interactions as revealed by systematic yeast two-hybrid screens. A high-confidence subset of 991 interactions links 576 proteins. To derive further biological insights from our data, we constructed an integrated network of proteins involved in DNA metabolism. Combining our data with additional evidences, we provide improved annotations for at least 18 proteins (including TP0004, TP0050, and TP0183 which are suggested to be involved in DNA metabolism). We estimate that this "minimal" bacterium contains on the order of 3,000 protein interactions. Profiles of functional interconnections indicate that bacterial proteins interact more promiscuously than eukaryotic proteins, reflecting the non-compartmentalized structure of the bacterial cell. Using our high-confidence interactions, we also predict 417,329 homologous interactions ("interologs") for 372 completely sequenced genomes and provide evidence that at least one third of them can be experimentally confirmed.     

Identification of novel components of NAD-utilizing metabolic pathways and prediction of their biochemical functions

Nicotinamide adenine dinucleotide (NAD) is a ubiquitous cofactor participating in numerous redox reactions. It is also a substrate for regulatory modifications of proteins and nucleic acids via the addition of ADP-ribose moieties or removal of acyl groups by transfer to ADP-ribose. In this study, we use in-depth sequence, structure and genomic context analysis to uncover new enzymes and substrate-binding proteins in NAD-utilizing metabolic and macromolecular modification systems. We predict that Escherichia coli YbiA and related families of domains from diverse bacteria, eukaryotes, large DNA viruses and single strand RNA viruses are previously unrecognized components of NAD-utilizing pathways that probably operate on ADP-ribose derivatives. Using contextual analysis we show that some of these proteins potentially act in RNA repair, where NAD is used to remove 2'-3' cyclic phosphodiester linkages. Likewise, we predict that another family of YbiA-related enzymes is likely to comprise a novel NAD-dependent ADP-ribosylation system for proteins, in conjunction with a previously unrecognized ADP-ribosyltransferase. A similar ADP-ribosyltransferase is also coupled with MACRO or ADP-ribosylglycohydrolase domain proteins in other related systems, suggesting that all these novel systems are likely to comprise pairs of ADP-ribosylation and ribosylglycohydrolase enzymes analogous to the DraG-DraT system, and a novel group of bacterial polymorphic toxins. We present evidence that some of these coupled ADP-ribosyltransferases/ribosylglycohydrolases are likely to regulate certain restriction modification enzymes in bacteria. The ADP-ribosyltransferases found in these, the bacterial polymorphic toxin and host-directed toxin systems of bacteria such as Waddlia also throw light on the evolution of this fold and the origin of eukaryotic polyADP-ribosyltransferases and NEURL4-like ARTs, which might be involved in centrosomal assembly. We also infer a novel biosynthetic pathway that might be involved in the synthesis of a nicotinate-derived compound in conjunction with an asparagine synthetase and AMPylating peptide ligase. We use the data derived from this analysis to understand the origin and early evolutionary trajectories of key NAD-utilizing enzymes and present targets for future biochemical investigations.     

Sequence, structure and functional diversity of PD-(D/E)XK phosphodiesterase superfamily

Proteins belonging to PD-(D/E)XK phosphodiesterases constitute a functionally diverse superfamily with representatives involved in replication, restriction, DNA repair and tRNA–intron splicing. Their malfunction in humans triggers severe diseases, such as Fanconi anemia and Xeroderma pigmentosum. To date there have been several attempts to identify and classify new PD-(D/E)KK phosphodiesterases using remote homology detection methods. Such efforts are complicated, because the superfamily exhibits extreme sequence and structural divergence. Using advanced homology detection methods supported with superfamily-wide domain architecture and horizontal gene transfer analyses, we provide a comprehensive reclassification of proteins containing a PD-(D/E)XK domain. The PD-(D/E)XK phosphodiesterases span over 21 900 proteins, which can be classified into 121 groups of various families. Eleven of them, including DUF4420, DUF3883, DUF4263, COG5482, COG1395, Tsp45I, HaeII, Eco47II, ScaI, HpaII and Replic_Relax, are newly assigned to the PD-(D/E)XK superfamily. Some groups of PD-(D/E)XK proteins are present in all domains of life, whereas others occur within small numbers of organisms. We observed multiple horizontal gene transfers even between human pathogenic bacteria or from Prokaryota to Eukaryota. Uncommon domain arrangements greatly elaborate the PD-(D/E)XK world. These include domain architectures suggesting regulatory roles in Eukaryotes, like stress sensing and cell-cycle regulation. Our results may inspire further experimental studies aimed at identification of exact biological functions, specific substrates and molecular mechanisms of reactions performed by these highly diverse proteins.     

Comparative genomics reveals insights into cyanobacterial evolution and habitat adaptation

Cyanobacteria are photosynthetic prokaryotes that inhabit diverse aquatic and terrestrial environments. However, the evolutionary mechanisms involved in the cyanobacterial habitat adaptation remain poorly understood. Here, based on phylogenetic and comparative genomic analyses of 650 cyanobacterial genomes, we investigated the genetic basis of cyanobacterial habitat adaptation (marine, freshwater, and terrestrial). We show: (1) the expansion of gene families is a common strategy whereby terrestrial cyanobacteria cope with fluctuating environments, whereas the genomes of many marine strains have undergone contraction to adapt to nutrient-poor conditions. (2) Hundreds of genes are strongly associated with specific habitats. Genes that are differentially abundant in genomes of marine, freshwater, and terrestrial cyanobacteria were found to be involved in light sensing and absorption, chemotaxis, nutrient transporters, responses to osmotic stress, etc., indicating the importance of these genes in the survival and adaptation of organisms in specific habitats. (3) A substantial fraction of genes that facilitate the adaptation of Cyanobacteria to specific habitats are contributed by horizontal gene transfer, and such genetic exchanges are more frequent in terrestrial cyanobacteria. Collectively, our results further our understandings of the adaptations of Cyanobacteria to different environments, highlighting the importance of ecological constraints imposed by the environment in shaping the evolution of Cyanobacteria.Subject terms: Phylogenetics, Microbial ecology     

Comprehensive classification of nucleotidyltransferase fold proteins: identification of novel families and their representatives in human

This article presents a comprehensive review of large and highly diverse superfamily of nucleotidyltransferase fold proteins by providing a global picture about their evolutionary history, sequence-structure diversity and fulfilled functional roles. Using top-of-the-line homology detection method combined with transitive searches and fold recognition, we revised the realm of these superfamily in numerous databases of catalogued protein families and structures, and identified 10 new families of nucleotidyltransferase fold. These families include hundreds of previously uncharacterized and various poorly annotated proteins such as Fukutin/LICD, NFAT, FAM46, Mab-21 and NRAP. Some of these proteins seem to play novel important roles, not observed before for this superfamily, such as regulation of gene expression or choline incorporation into cell membrane. Importantly, within newly detected families we identified 25 novel superfamily members in human genome. Among these newly assigned members are proteins known to be involved in congenital muscular dystrophy, neurological diseases and retinal pigmentosa what sheds some new light on the molecular background of these genetic disorders. Twelve of new human nucleotidyltransferase fold proteins belong to Mab-21 family known to be involved in organogenesis and development. The determination of specific biological functions of these newly detected proteins remains a challenging task.     

The multifaceted roles of the HORMA domain in cellular signaling

The HORMA domain is a multifunctional protein–protein interaction module found in diverse eukaryotic signaling pathways including the spindle assembly checkpoint, numerous DNA recombination/repair pathways, and the initiation of autophagy. In all of these pathways, HORMA domain proteins occupy key signaling junctures and function through the controlled assembly and disassembly of signaling complexes using a stereotypical “safety belt” peptide interaction mechanism. A recent explosion of structural and functional work has shed new light on these proteins, illustrating how strikingly similar structural mechanisms give rise to radically different functional outcomes in each family of HORMA domain proteins.     

Multivariate selection shapes environment-dependent variation in the clonal morphology of a red seaweed

Within-individual strategies of variation (e.g., phenotypic plasticity) are particularly relevant to modular organisms, in which ramets of the same genetic individual may encounter diverse environments imposing diverse patterns of selection. Hence, measuring selection in heterogeneous environments is essential to understanding whether environment-dependent phenotypic change enhances the fitness of modular individuals. In sublittoral marine habitats, competition for light and space among modular taxa generates extreme patchiness in resource availability. Little is known, however, of the potential for plasticity within individuals to arise from spatially-variable selection in such systems. We tested whether plasticity enhances genet-level fitness in Asparagopsis armata, a clonal seaweed in which correlated traits mediate morphological responses to variation in light. Using the capacity for rapid, clonal growth to measure fitness, we identified aspects of ramet morphology targeted by selection in two contrasting light environments and compared patterns of selection across environments. We found that directional selection on single traits, coupled with linear and nonlinear selection on multi-trait interactions, shape ramet morphology within environments and favor different phenotypes in each. Evidence of environment-dependent, multivariate selection on correlated traits is novel for any marine modular organism and demonstrates that seaweeds, such as A. armata, may potentially adapt to environmental heterogeneity via plasticity in clonal morphology.     

Diverse heliorhodopsins detected via functional metagenomics in freshwater Actinobacteria,Chloroflexi and Archaea

The recently discovered rhodopsin family of heliorhodopsins (HeRs) is abundant in diverse microbial environments. So far, the functional and biological roles of HeRs remain unknown. To tackle this issue, we combined experimental and computational screens to gain some novel insights. Here, 10 readily expressed HeR genes were found using functional metagenomics on samples from two freshwater environments. These HeRs originated from diverse prokaryotic groups: Actinobacteria, Chloroflexi and Archaea. Heterologously expressed HeRs absorbed light in the green and yellow wavelengths (543–562 nm) and their photocycles exhibited diverse kinetic characteristics. To approach the physiological function of the HeRs, we used our environmental clones along with thousands of microbial genomes to analyze genes neighbouring HeRs. The strongest association was found with the DegV family involved in activation of fatty acids, which allowed us to hypothesize that HeRs might be involved in light-induced membrane lipid modifications.     

The Metagenomic Telescope

Next generation sequencing technologies led to the discovery of numerous new microbe species in diverse environmental samples. Some of the new species contain genes never encountered before. Some of these genes encode proteins with novel functions, and some of these genes encode proteins that perform some well-known function in a novel way. A tool, named the Metagenomic Telescope, is described here that applies artificial intelligence methods, and seems to be capable of identifying new protein functions even in the well-studied model organisms. As a proof-of-principle demonstration of the Metagenomic Telescope, we considered DNA repair enzymes in the present work. First we identified proteins in DNA repair in well–known organisms (i.e., proteins in base excision repair, nucleotide excision repair, mismatch repair and DNA break repair); next we applied multiple alignments and then built hidden Markov profiles for each protein separately, across well–researched organisms; next, using public depositories of metagenomes, originating from extreme environments, we identified DNA repair genes in the samples. While the phylogenetic classification of the metagenomic samples are not typically available, we hypothesized that some very special DNA repair strategies need to be applied in bacteria and Archaea living in those extreme circumstances. It is a difficult task to evaluate the results obtained from mostly unknown species; therefore we applied again the hidden Markov profiling: for the identified DNA repair genes in the extreme metagenomes, we prepared new hidden Markov profiles (for each genes separately, subsequent to a cluster analysis); and we searched for similarities to those profiles in model organisms. We have found well known DNA repair proteins, numerous proteins with unknown functions, and also proteins with known, but different functions in the model organisms.     

Hookworm SCP/TAPS protein structure--A key to understanding host-parasite interactions and developing new interventions

SCP/TAPS proteins are a diverse family of molecules in eukaryotes, including parasites. Despite their abundant occurrence in parasite secretomes, very little is known about their functions in parasitic nematodes, including blood-feeding hookworms. Current information indicates that SCP/TAPS proteins (called Ancylostoma-secreted proteins, ASPs) of the canine hookworm, Ancylostoma caninum, represent at least three distinct groups of proteins. This information, combined with comparative modelling, indicates that all known ASPs have an equatorial groove that binds extended structures, such as peptides or glycans. To elucidate structure-function relationships, we explored the three-dimensional crystal structure of an ASP (called Ac-ASP-7), which is highly up-regulated in expression in the transition of A. caninum larvae from a free-living to a parasitic stage. The topology of the N-terminal domain is consistent with pathogenesis-related proteins, and the C-terminal extension that resembles the fold of the Hinge domain. By anomalous diffraction, we identified a new metal binding site in the C-terminal extension of the protein. Ac-ASP-7 is in a monomer-dimer equilibrium, and crystal-packing analysis identified a dimeric structure which might resemble the homo-dimer in solution. The dimer interaction interface includes a novel binding site for divalent metal ions, and is proposed to serve as a binding site for proteins involved in the parasite-host interplay at the molecular level. Understanding this interplay and the integration of structural and functional data could lead to the design of new approaches for the control of parasitic diseases, with biotechnological outcomes.     

The SWIB and the MDM2 domains are homologous and share a common fold

Using a novel algorithm for protein sequence/structure analysis, we propose a probable homology between the SWIB domain (the conserved domain of the 60 kda subunit of the SWIB complex involved in chromatin remodelling) and the p53-binding domain of the MDM2 oncoprotein. The homology suggests that the SWIB domain would adopt a structure similar to that of the MDM2 domain and that these two families of proteins may share a similar functional mechanism.     

Metagenomics of microbial and viral life in terrestrial geothermal environments

Geothermally heated regions of Earth, such as terrestrial volcanic areas (fumaroles, hot springs, and geysers) and deep-sea hydrothermal vents, represent a variety of different environments populated by extremophilic archaeal and bacterial microorganisms. Since most of these microbes thriving in such harsh biotopes, they are often recalcitrant to cultivation; therefore, ecological, physiological and phylogenetic studies of these microbial populations have been hampered for a long time. More recently, culture-independent methodologies coupled with the fast development of next generation sequencing technologies as well as with the continuous advances in computational biology, have allowed the production of large amounts of metagenomic data. Specifically, these approaches have assessed the phylogenetic composition and functional potential of microbial consortia thriving within these habitats, shedding light on how extreme physico-chemical conditions and biological interactions have shaped such microbial communities. Metagenomics allowed to better understand that the exposure to an extreme range of selective pressures in such severe environments, accounts for genomic flexibility and metabolic versatility of microbial and viral communities, and makes extreme- and hyper-thermophiles suitable for bioprospecting purposes, representing an interesting source for novel thermostable proteins that can be potentially used in several industrial processes.     

Linking yeast Gcn5p catalytic function and gene regulation using a quantitative, graded dominant mutant approach

Establishing causative links between protein functional domains and global gene regulation is critical for advancements in genetics, biotechnology, disease treatment, and systems biology. This task is challenging for multifunctional proteins when relying on traditional approaches such as gene deletions since they remove all domains simultaneously. Here, we describe a novel approach to extract quantitative, causative links by modulating the expression of a dominant mutant allele to create a function-specific competitive inhibition. Using the yeast histone acetyltransferase Gcn5p as a case study, we demonstrate the utility of this approach and (1) find evidence that Gcn5p is more involved in cell-wide gene repression, instead of the accepted gene activation associated with HATs, (2) identify previously unknown gene targets and interactions for Gcn5p-based acetylation, (3) quantify the strength of some Gcn5p-DNA associations, (4) demonstrate that this approach can be used to correctly identify canonical chromatin modifications, (5) establish the role of acetyltransferase activity on synthetic lethal interactions, and (6) identify new functional classes of genes regulated by Gcn5p acetyltransferase activity--all six of these major conclusions were unattainable by using standard gene knockout studies alone. We recommend that a graded dominant mutant approach be utilized in conjunction with a traditional knockout to study multifunctional proteins and generate higher-resolution data that more accurately probes protein domain function and influence.     

Identification of novel restriction endonuclease-like fold families among hypothetical proteins

Restriction endonucleases and other nucleic acid cleaving enzymes form a large and extremely diverse superfamily that display little sequence similarity despite retaining a common core fold responsible for cleavage. The lack of significant sequence similarity between protein families makes homology inference a challenging task and hinders new family identification with traditional sequence-based approaches. Using the consensus fold recognition method Meta-BASIC that combines sequence profiles with predicted protein secondary structure, we identify nine new restriction endonuclease-like fold families among previously uncharacterized proteins and predict these proteins to cleave nucleic acid substrates. Application of transitive searches combined with gene neighborhood analysis allow us to confidently link these unknown families to a number of known restriction endonuclease-like structures and thus assign folds to the uncharacterized proteins. Finally, our method identifies a novel restriction endonuclease-like domain in the C-terminus of RecC that is not detected with structure-based searches of the existing PDB database.