Resilience of biochemical activity in protein domains in the face of structural divergence

The Lipocone Superfamily: A Unifying Theme In Metabolism Of Lipids, Peptidoglycan And Exopolysaccharides, Inter-Organismal Conflicts And Immunity

Wnt proteins are critical signaling molecules in developmental processes across animals. Despite intense study, their evolutionary roots have remained enigmatic. Using sensitive sequence analysis and structure modeling, we establish that the Wnts are part of a vast assemblage of domains, the Lipocone superfamily, defined here for the first time. It includes previously studied enzymatic domains like the phosphatidylserine synthases (PTDSS1/2) and the TelC toxin domain from Streptococcus intermedius , the enigmatic VanZ proteins, the animal Serum Amyloid A (SAA) and a further host of uncharacterized proteins in a total of 30 families. Though the metazoan Wnts are catalytically inactive, we present evidence for a conserved active site across this superfamily, versions of which are consistently predicted to operate on head groups of either phospholipids or polyisoprenoid lipids, catalyzing transesterification and phosphate-containing head group severance reactions. We argue that this superfamily originated as membrane proteins, with one branch (including Wnt and SAA) evolving into soluble versions. By comprehensively analyzing contextual information networks derived from comparative genomics, we establish that they act in varied functional contexts, including regulation of membrane lipid composition, extracellular polysaccharide biosynthesis, and biogenesis of bacterial outer-membrane components, like lipopolysaccharides. On multiple occasions, members of this superfamily, including the bacterial progenitors of Wnt and SAA, have been recruited as effectors in biological conflicts spanning inter-organismal interactions and anti-viral immunity in both prokaryotes and eukaryotes. These findings establish a unifying theme in lipid biochemistry, explain the origins of Wnt signaling and provide new leads regarding immunity across the tree of life.

New biochemistry in the Rhodanese-phosphatase superfamily: emerging roles in diverse metabolic processes, nucleic acid modifications, and biological conflicts

The protein-tyrosine/dual-specificity phosphatases and rhodanese domains constitute a sprawling superfamily of Rossmannoid domains that use a conserved active site with a cysteine to catalyze a range of phosphate-transfer, thiotransfer, selenotransfer and redox activities. While these enzymes have been extensively studied in the context of protein/lipid head group dephosphorylation and various thiotransfer reactions, their overall diversity and catalytic potential remain poorly understood. Using comparative genomics and sequence/structure analysis, we comprehensively investigate and develop a natural classification for this superfamily. As a result, we identified several novel clades, both those which retain the catalytic cysteine and those where a distinct active site has emerged in the same location (e.g. diphthine synthase-like methylases and RNA 2' OH ribosyl phosphate transferases). We also present evidence that the superfamily has a wider range of catalytic capabilities than previously known, including a set of parallel activities operating on various sugar/sugar alcohol groups in the context of NAD+-derivatives and RNA termini, and potential phosphate transfer activities involving sugars and nucleotides. We show that such activities are particularly expanded in the RapZ-C-DUF488-DUF4326 clade, defined here for the first time. Some enzymes from this clade are predicted to catalyze novel DNA-end processing activities as part of nucleic-acid-modifying systems that are likely to function in biological conflicts between viruses and their hosts.

Ribonuclease T2 represents a distinct circularly permutated version of the BECR RNases

Detection of homologous relationships among proteins and understanding their mechanisms of diversification are major topics in the fields of protein science, bioinformatics, and phylogenetics. Recent developments in sequence/profile-based and structural similarity-based methods have greatly facilitated the unification and classification of many protein families into superfamilies or folds, yet many proteins remain unclassified in current protein databases. As one of the three earliest identified RNases in biology, ribonuclease T2, also known as RNase I in Escherichia coli, RNase Rh in fungi, or S-RNase in plant, is thought to be an ancient RNase family due to its widespread distribution and distinct structure. In this study, we present evidence that RNase T2 represents a circularly permutated version of the BECR (Barnase-EndoU-Colicin E5/D-RelE) fold RNases. This subtle relationship cannot be detected by traditional methods such as sequence/profile-based comparisons, structure-similarity searches, and circular permutation detections. However, we were able to identify the structural similarity using rational reconstruction of a theoretical RNase T2 ancestor via a reverse circular permutation process, followed by structural modeling using AlphaFold2, and structural comparisons. This relationship is further supported by the fact that RNase T2 and other typical BECR RNases, namely Colicin D, RNase A, and BrnT, share similar catalytic site configurations, all involving an analogous set of conserved residues on the α0 helix and the β4 strand of the BECR fold. This study revealed a hidden root of RNase T2 in bacterial toxin systems and demonstrated that reconstruction and modeling of ancestral topology is an effective strategy to identify remote relationship between proteins.

Functional and Structural Diversity of Bacterial Contact-Dependent Growth Inhibition Effectors

Bacteria live in complex communities and environments, competing for space and nutrients. Within their niche habitats, bacteria have developed various inter-bacterial mechanisms to compete and communicate. One such mechanism is contact-dependent growth inhibition (CDI). CDI is found in many Gram-negative bacteria, including several pathogens. These CDI+ bacteria encode a CdiB/CdiA two-partner secretion system that delivers inhibitory toxins into neighboring cells upon contact. Toxin translocation results in the growth inhibition of closely related strains and provides a competitive advantage to the CDI+ bacteria. CdiB, an outer-membrane protein, secretes CdiA onto the surface of the CDI+ bacteria. When CdiA interacts with specific target-cell receptors, CdiA delivers its C-terminal toxin region (CdiA-CT) into the target-cell. CdiA-CT toxin proteins display a diverse range of toxic functions, such as DNase, RNase, or pore-forming toxin activity. CDI+ bacteria also encode an immunity protein, CdiI, that specifically binds and neutralizes its cognate CdiA-CT, protecting the CDI+ bacteria from auto-inhibition. In Gram-negative bacteria, toxin/immunity (CdiA-CT/CdiI) pairs have highly variable sequences and functions, with over 130 predicted divergent toxin/immunity complex families. In this review, we will discuss biochemical and structural advances made in the characterization of CDI. This review will focus on the diverse array of CDI toxin/immunity complex structures together with their distinct toxin functions. Additionally, we will discuss the most recent studies on target-cell recognition and toxin entry, along with the discovery of a new member of the CDI loci. Finally, we will offer insights into how these diverse toxin/immunity complexes could be harnessed to fight human diseases.

BioS2Net: Holistic Structural and Sequential Analysis of Biomolecules Using a Deep Neural Network

BACKGROUND

For decades, the rate of solving new biomolecular structures has been exceeding that at which their manual classification and feature characterisation can be carried out efficiently. Therefore, a new comprehensive and holistic tool for their examination is needed.

METHODS

Here we propose the Biological Sequence and Structure Network (BioS2Net), which is a novel deep neural network architecture that extracts both sequential and structural information of biomolecules. Our architecture consists of four main parts: (i) a sequence convolutional extractor, (ii) a 3D structure extractor, (iii) a 3D structure-aware sequence temporal network, as well as (iv) a fusion and classification network.

RESULTS

We have evaluated our approach using two protein fold classification datasets. BioS2Net achieved a 95.4% mean class accuracy on the eDD dataset and a 76% mean class accuracy on the F184 dataset. The accuracy of BioS2Net obtained on the eDD dataset was comparable to results achieved by previously published methods, confirming that the algorithm described in this article is a top-class solution for protein fold recognition.

CONCLUSIONS

BioS2Net is a novel tool for the holistic examination of biomolecules of known structure and sequence. It is a reliable tool for protein analysis and their unified representation as feature vectors.

Photoglobin, a distinct family of non-heme binding globins, defines a potential photosensor in prokaryotic signal transduction systems

Globins constitute an ancient superfamily of proteins, exhibiting enormous structural and functional diversity, as demonstrated by many heme-binding families and two non-heme binding families that were discovered in bacterial stressosome component RsbR and in light-harvesting phycobiliproteins (phycocyanin) in cyanobacteria and red algae. By comprehensively exploring the globin repertoire using sensitive computational analyses of sequences, structures, and genomes, we present the identification of the third family of non-heme binding globins—the photoglobin. By conducting profile-based comparisons, clustering analyses, and structural modeling, we demonstrate that photoglobin is related to, but distinct from, the phycocyanin family. Photoglobin preserves a potential ligand-binding pocket, whose residue configuration closely resembles that of phycocyanin, indicating that photoglobin potentially binds to a comparable linear tetrapyrrole. By exploring the contextual information provided by the photoglobin’s domain architectures and gene-neighborhoods, we found that photoglobin is frequently associated with the B12-binding light sensor domain and many domains typical of prokaryotic signal transduction systems. Structural modeling using AlphaFold2 demonstrated that photoglobin and B12-binding domains form a structurally conserved hub among different domain architecture contexts. Based on these strong associations, we predict that the coupled photoglobin and B12-binding domains act as a light-sensing regulatory bundle, with each domain sensing different wavelengths of light resulting in switch-like regulation of downstream signaling effectors. Thus, based on the above lines of evidence, we present a distinct non-heme binding globin family and propose that it may define a new type of light sensor, by means of a linear tetrapyrrole, in complex prokaryotic signal transduction systems.

Genomic discovery and structural dissection of a novel type of polymorphic toxin system in gram-positive bacteria

Bacteria have developed several molecular conflict systems to facilitate kin recognition and non-kin competition to gain advantages in the acquisition of growth niches and of limited resources. One such example is a large class of so-called polymorphic toxin systems (PTSs), which comprise a variety of the toxin proteins secreted via T2SS, T5SS, T6SS, T7SS and many others. These systems are highly divergent in terms of sequence/structure, domain architecture, toxin-immunity association, and organization of the toxin loci, which makes it difficult to identify and characterize novel systems using traditional experimental and bioinformatic strategies. In recent years, we have been developing and utilizing unique genome-mining strategies and pipelines, based on the organizational principles of both domain architectures and genomic loci of PTSs, for an effective and comprehensive discovery of novel PTSs, dissection of their components, and prediction of their structures and functions. In this study, we present our systematic discovery of a new type of PTS (S8-PTS) in several gram-positive bacteria. We show that the S8-PTS contains three components: a peptidase of the S8 family (subtilases), a polymorphic toxin, and an immunity protein. We delineated the typical organization of these polymorphic toxins, in which a N-terminal signal peptide is followed by a potential receptor binding domain, BetaH, and one of 16 toxin domains. We classified each toxin domain by the distinct superfamily to which it belongs, identifying nine BECR ribonucleases, one Restriction Endonuclease, one HNH nuclease, two novel toxin domains homologous to the VOC enzymes, one toxin domain with the Frataxin-like fold, and several other unique toxin families such as Ntox33 and HicA. Accordingly, we identified 20 immunity families and classified them into different classes of folds. Further, we show that the S8-PTS-associated peptidases are analogous to many other processing peptidases found in T5SS, T7SS, T9SS, and many proprotein-processing peptidases, indicating that they function to release the toxin domains during secretion. The S8-PTSs are mostly found in animal and plant-associated bacteria, including many pathogens. We propose S8-PTSs will facilitate the competition of these bacteria with other microbes or contribute to the pathogen-host interactions.

Comparative Phylogenomic Analysis Reveals Evolutionary Genomic Changes and Novel Toxin Families in Endophytic Liberibacter Pathogens

Liberibacter pathogens are the causative agents of several severe crop diseases worldwide, including citrus Huanglongbing and potato zebra chip. These bacteria are endophytic and nonculturable, which makes experimental approaches challenging and highlights the need for bioinformatic analysis in advancing our understanding about Liberibacter pathogenesis. Here, we performed an in-depth comparative phylogenomic analysis of the Liberibacter pathogens and their free-living, nonpathogenic, ancestral species, aiming to identify major genomic changes and determinants associated with their evolutionary transitions in living habitats and pathogenicity. Using gene neighborhood analysis and phylogenetic classification, we systematically uncovered, annotated, and classified all prophage loci into four types, including one previously unrecognized group. We showed that these prophages originated through independent gene transfers at different evolutionary stages of Liberibacter and only the SC-type prophage was associated with the emergence of the pathogens. Using ortholog clustering, we vigorously identified two additional sets of genomic genes, which were either lost or gained in the ancestor of the pathogens. Consistent with the habitat change, the lost genes were enriched for biosynthesis of cellular building blocks. Importantly, among the gained genes, we uncovered several previously unrecognized toxins, including new toxins homologous to the EspG/VirA effectors, a YdjM phospholipase toxin, and a secreted endonuclease/exonuclease/phosphatase (EEP) protein. Our results substantially extend the knowledge of the evolutionary events and potential determinants leading to the emergence of endophytic, pathogenic Liberibacter species, which will facilitate the design of functional experiments and the development of new methods for detection and blockage of these pathogens. IMPORTANCELiberibacter pathogens are associated with several severe crop diseases, including citrus Huanglongbing, the most destructive disease to the citrus industry. Currently, no effective cure or treatments are available, and no resistant citrus variety has been found. The fact that these obligate endophytic pathogens are not culturable has made it extremely challenging to experimentally uncover the genes/proteins important to Liberibacter pathogenesis. Further, earlier bioinformatics studies failed to identify key genomic determinants, such as toxins and effector proteins, that underlie the pathogenicity of the bacteria. In this study, an in-depth comparative genomic analysis of Liberibacter pathogens along with their ancestral nonpathogenic species identified the prophage loci and several novel toxins that are evolutionarily associated with the emergence of the pathogens. These results shed new light on the disease mechanism of Liberibacter pathogens and will facilitate the development of new detection and blockage methods targeting the toxins.

Resolving the Rules of Robustness and Resilience in Biology Across Scales

Why do some biological systems and communities persist while others fail? Robustness, a system's stability, and resilience, the ability to return to a stable state, are key concepts that span multiple disciplines within and outside the biological sciences. Discovering and applying common rules that govern the robustness and resilience of biological systems is a critical step toward creating solutions for species survival in the face of climate change, as well as the for the ever-increasing need for food, health, and energy for human populations. We propose that network theory provides a framework for universal scalable mathematical models to describe robustness and resilience and the relationship between them, and hypothesize that resilience at lower organization levels contribute to robust systems. Insightful models of biological systems can be generated by quantifying the mechanisms of redundancy, diversity, and connectivity of networks, from biochemical processes to ecosystems. These models provide pathways towards understanding how evolvability can both contribute to and result from robustness and resilience under dynamic conditions. We now have an abundance of data from model and non-model systems and the technological and computational advances for studying complex systems. Several conceptual and policy advances will allow the research community to elucidate the rules of robustness and resilience. Conceptually, a common language and data structure that can be applied across levels of biological organization needs to be developed. Policy advances such as cross-disciplinary funding mechanisms, access to affordable computational capacity, and the integration of network theory and computer science within the standard biological science curriculum will provide the needed research environments. This new understanding of biological systems will allow us to derive ever more useful forecasts of biological behaviors and revolutionize the engineering of biological systems that can survive changing environments or disease, navigate the deepest oceans, or sustain life throughout the solar system.

Unification and extensive diversification of M/Orf3-related ion channel proteins in coronaviruses and other nidoviruses

The coronavirus, Severe Acute Respiratory Syndrome (SARS)-CoV-2, responsible for the ongoing coronavirus disease 2019 (COVID-19) pandemic, has emphasized the need for a better understanding of the evolution of virus-host interactions. ORF3a in both SARS-CoV-1 and SARS-CoV-2 are ion channels (viroporins) implicated in virion assembly and membrane budding. Using sensitive profile-based homology detection methods, we unify the SARS-CoV ORF3a family with several families of viral proteins, including ORF5 from MERS-CoVs, proteins from beta-CoVs (ORF3c), alpha-CoVs (ORF3b), most importantly, the Matrix (M) proteins from CoVs, and more distant homologs from other nidoviruses. We present computational evidence that these viral families might utilize specific conserved polar residues to constitute an aqueous pore within the membrane-spanning region. We reconstruct an evolutionary history of these families and objectively establish the common origin of the M proteins of CoVs and Toroviruses. We also show that the divergent ORF3 clade (ORF3a/ORF3b/ORF3c/ORF5 families) represents a duplication stemming from the M protein in alpha- and beta-CoVs. By phyletic profiling of major structural components of primary nidoviruses, we present a hypothesis for their role in virion assembly of CoVs, ToroVs, and Arteriviruses. The unification of diverse M/ORF3 ion channel families in a wide range of nidoviruses, especially the typical M protein in CoVs, reveal a conserved, previously under-appreciated role of ion channels in virion assembly and membrane budding. We show that M and ORF3 are under different evolutionary pressures; in contrast to the slow evolution of M as core structural component, the ORF3 clade is under selection for diversification, which suggests it might act at the interface with host molecules and/or immune attack.

An RNA Repair Operon Regulated by Damaged tRNAs

Many bacteria contain an RNA repair operon, encoding the RtcB RNA ligase and the RtcA RNA cyclase, that is regulated by the RtcR transcriptional activator. Although RtcR contains a divergent version of the CARF (CRISPR-associated Rossman fold) oligonucleotide-binding regulatory domain, both the specific signal that regulates operon expression and the substrates of the encoded enzymes are unknown. We report that tRNA fragments activate operon expression. Using a genetic screen in Salmonella enterica serovar Typhimurium, we find that the operon is expressed in the presence of mutations that cause tRNA fragments to accumulate. RtcA, which converts RNA phosphate ends to 2′, 3′-cyclic phosphate, is also required. Operon expression and tRNA fragment accumulation also occur upon DNA damage. The CARF domain binds 5′ tRNA fragments ending in cyclic phosphate, and RtcR oligomerizes upon binding these ligands, a prerequisite for operon activation. Our studies reveal a signaling pathway involving broken tRNAs and implicate the operon in tRNA repair.

Hughes et al. demonstrate that a bacterial RNA repair operon, containing the RtcB RNA ligase and the RtcA RNA cyclase, is regulated by binding of 5′ tRNA halves ending in 2′, 3′-cyclic phosphate to the RtcR transcriptional activator. These studies show how tRNA fragments can regulate bacterial gene expression.

Structural bioinformatics predicts that the Retinitis Pigmentosa-28 protein of unknown function FAM161A is a homologue of the microtubule nucleation factor Tpx2

Background: FAM161A is a microtubule-associated protein conserved widely across eukaryotes, which is mutated in the inherited blinding disease Retinitis Pigmentosa-28. FAM161A is also a centrosomal protein, being a core component of a complex that forms an internal skeleton of centrioles. Despite these observations about the importance of FAM161A, current techniques used to examine its sequence reveal no homologies to other proteins.

Methods: Sequence profiles derived from multiple sequence alignments of FAM161A homologues were constructed by PSI-BLAST and HHblits, and then used by the profile-profile search tool HHsearch, implemented online as HHpred, to identify homologues. These in turn were used to create profiles for reverse searches and pair-wise searches. Multiple sequence alignments were also used to identify amino acid usage in functional elements.

Results: FAM161A has a single homologue: the targeting protein for Xenopus kinesin-like protein-2 (Tpx2), which is a strong hit across more than 200 residues. Tpx2 is also a microtubule-associated protein, and it has been shown previously by a cryo-EM molecular structure to nucleate microtubules through two small elements: an extended loop and a short helix. The homology between FAM161A and Tpx2 includes these elements, as FAM161A has three copies of the loop, and one helix that has many, but not all, properties of the one in Tpx2.

Conclusions: FAM161A and ­its homologues are predicted to be a previously unknown variant of Tpx2, and hence bind microtubules in the same way. This prediction allows precise, testable molecular models to be made of FAM161A-microtubule complexes.

A malaria parasite subtilisin propeptide-like protein is a potent inhibitor of the egress protease SUB1

Subtilisin-like serine peptidases (subtilases) play important roles in the life cycle of many organisms, including the protozoan parasites that are the causative agent of malaria, Plasmodium spp. As with other peptidases, subtilase proteolytic activity has to be tightly regulated in order to prevent potentially deleterious uncontrolled protein degradation. Maturation of most subtilases requires the presence of an N-terminal propeptide that facilitates folding of the catalytic domain. Following its proteolytic cleavage, the propeptide acts as a transient, tightly bound inhibitor until its eventual complete removal to generate active protease. Here we report the identification of a stand-alone malaria parasite propeptide-like protein, called SUB1-ProM, encoded by a conserved gene that lies in a highly syntenic locus adjacent to three of the four subtilisin-like genes in the Plasmodium genome. Template-based modelling and ab initio structure prediction showed that the SUB1-ProM core structure is most similar to the X-ray crystal structure of the propeptide of SUB1, an essential parasite subtilase that is discharged into the parasitophorous vacuole (PV) to trigger parasite release (egress) from infected host cells. Recombinant Plasmodium falciparum SUB1-ProM was found to be a fast-binding, potent inhibitor of P. falciparum SUB1, but not of the only other essential blood-stage parasite subtilase, SUB2, or of other proteases examined. Mass-spectrometry and immunofluorescence showed that SUB1-ProM is expressed in the PV of blood stage P. falciparum, where it may act as an endogenous inhibitor to regulate SUB1 activity in the parasite.

Small protein folds at the root of an ancient metabolic network

Life on Earth is driven by electron transfer reactions catalyzed by a suite of enzymes that comprise the superfamily of oxidoreductases (Enzyme Classification EC1). Most modern oxidoreductases are complex in their structure and chemistry and must have evolved from a small set of ancient folds. Ancient oxidoreductases from the Archean Eon between ca. 3.5 and 2.5 billion years ago have been long extinct, making it challenging to retrace evolution by sequence-based phylogeny or ancestral sequence reconstruction. However, three-dimensional topologies of proteins change more slowly than sequences. Using comparative structure and sequence profile-profile alignments, we quantify the similarity between proximal cofactor-binding folds and show that they are derived from a common ancestor. We discovered that two recurring folds were central to the origin of metabolism: ferredoxin and Rossmann-like folds. In turn, these two folds likely shared a common ancestor that, through duplication, recruitment, and diversification, evolved to facilitate electron transfer and catalysis at a very early stage in the origin of metabolism.

Highly regulated, diversifying NTP-dependent biological conflict systems with implications for the emergence of multicellularity

Social cellular aggregation or multicellular organization pose increased risk of transmission of infections through the system upon infection of a single cell. The generality of the evolutionary responses to this outside of Metazoa remains unclear. We report the discovery of several thematically unified, remarkable biological conflict systems preponderantly present in multicellular prokaryotes. These combine thresholding mechanisms utilizing NTPase chaperones (the MoxR-vWA couple), GTPases and proteolytic cascades with hypervariable effectors, which vary either by using a reverse transcriptase-dependent diversity-generating system or through a system of acquisition of diverse protein modules, typically in inactive form, from various cellular subsystems. Conciliant lines of evidence indicate their deployment against invasive entities, like viruses, to limit their spread in multicellular/social contexts via physical containment, dominant-negative interactions or apoptosis. These findings argue for both a similar operational 'grammar' and shared protein domains in the sensing and limiting of infections during the multiple emergences of multicellularity.

Convergent Evolution of the Barnase/EndoU/Colicin/RelE (BECR) Fold in Antibacterial tRNase Toxins

Contact-dependent growth inhibition (CDI) is a form of interbacterial competition mediated by CdiB-CdiA two-partner secretion systems. CdiA effector proteins carry polymorphic C-terminal toxin domains (CdiA-CT), which are neutralized by specific CdiI immunity proteins to prevent self-inhibition. Here, we present the crystal structures of CdiA-CT⋅CdiI complexes from Klebsiella pneumoniae 342 and Escherichia coli 3006. The toxins adopt related folds that resemble the ribonuclease domain of colicin D, and both are isoacceptor-specific tRNases that cleave the acceptor stem of deacylated tRNA_GAU^Ile. Although the toxins are similar in structure and substrate specificity, CdiA-CT^Kp342 activity requires translation factors EF-Tu and EF-Ts, whereas CdiA-CT^EC3006 is intrinsically active. Furthermore, the corresponding immunity proteins are unrelated in sequence and structure. CdiI^Kp342 forms a dimeric β sandwich, whereas CdiI^EC3006 is an α-solenoid monomer. Given that toxin-immunity genes co-evolve as linked pairs, these observations suggest that the similarities in toxin structure and activity reflect functional convergence.

Different Ways of Doing the Same: Variations in the Two Last Steps of the Purine Biosynthetic Pathway in Prokaryotes

The last two steps of the purine biosynthetic pathway may be catalyzed by different enzymes in prokaryotes. The genes that encode these enzymes include homologs of purH, purP, purO and those encoding the AICARFT and IMPCH domains of PurH, here named purV and purJ, respectively. In Bacteria, these reactions are mainly catalyzed by the domains AICARFT and IMPCH of PurH. In Archaea, these reactions may be carried out by PurH and also by PurP and PurO, both considered signatures of this domain and analogous to the AICARFT and IMPCH domains of PurH, respectively. These genes were searched for in 1,403 completely sequenced prokaryotic genomes publicly available. Our analyses revealed taxonomic patterns for the distribution of these genes and anticorrelations in their occurrence. The analyses of bacterial genomes revealed the existence of genes coding for PurV, PurJ, and PurO, which may no longer be considered signatures of the domain Archaea. Although highly divergent, the PurOs of Archaea and Bacteria show a high level of conservation in the amino acids of the active sites of the protein, allowing us to infer that these enzymes are analogs. Based on the results, we propose that the gene purO was present in the common ancestor of all living beings, whereas the gene encoding PurP emerged after the divergence of Archaea and Bacteria and their isoforms originated in duplication events in the common ancestor of phyla Crenarchaeota and Euryarchaeota. The results reported here expand our understanding of the diversity and evolution of the last two steps of the purine biosynthetic pathway in prokaryotes.

Oxidative opening of the aromatic ring: Tracing the natural history of a large superfamily of dioxygenase domains and their relatives

A diverse collection of enzymes comprising the protocatechuate dioxygenases (PCADs) has been characterized in several extradiol aromatic compound degradation pathways. Structural studies have shown a relationship between PCADs and the more broadly-distributed, functionally enigmatic Memo domain linked to several human diseases. To better understand the evolution of this PCAD-Memo protein superfamily, we explored their structural and functional determinants to establish a unified evolutionary framework, identifying 15 clearly-delineable families, including a previously-underappreciated diversity in five Memo clade families. We place the superfamily's origin within the greater radiation of the nucleoside phosphorylase/hydrolase-peptide/amidohydrolase fold prior to the last universal common ancestor of all extant organisms. In addition to identifying active-site residues across the superfamily, we describe three distinct, structurally-variable regions emanating from the core scaffold often housing conserved residues specific to individual families. These were predicted to contribute to the active-site pocket, potentially in substrate specificity and allosteric regulation. We also identified several previously-undescribed conserved genome contexts, providing insight into potentially novel substrates in PCAD clade families. We extend known conserved contextual associations for the Memo clade beyond previously-described associations with the AMMECR1 domain and a radical S-adenosylmethionine family domain. These observations point to two distinct yet potentially overlapping contexts wherein the elusive molecular function of the Memo domain could be finally resolved, thereby linking it to nucleotide base and aliphatic isoprenoid modification. In total, this report throws light on the functions of large swaths of the experimentally-uncharacterized PCAD-Memo families.

Antimicrobial Peptides, Polymorphic Toxins, and Self-Nonself Recognition Systems in Archaea: an Untapped Armory for Intermicrobial Conflicts

Diverse and highly variable systems involved in biological conflicts and self-versus-nonself discrimination are ubiquitous in bacteria but much less studied in archaea. We performed comprehensive comparative genomic analyses of the archaeal systems that share components with analogous bacterial systems and propose an approach to identify new systems that could be involved in these functions. We predict polymorphic toxin systems in 141 archaeal genomes and identify new, archaea-specific toxin and immunity protein families. These systems are widely represented in archaea and are predicted to play major roles in interactions between species and in intermicrobial conflicts. This work is expected to stimulate experimental research to advance the understanding of poorly characterized major aspects of archaeal biology.

Numerous, diverse, highly variable defense and offense genetic systems are encoded in most bacterial genomes and are involved in various forms of conflict among competing microbes or their eukaryotic hosts. Here we focus on the offense and self-versus-nonself discrimination systems encoded by archaeal genomes that so far have remained largely uncharacterized and unannotated. Specifically, we analyze archaeal genomic loci encoding polymorphic and related toxin systems and ribosomally synthesized antimicrobial peptides. Using sensitive methods for sequence comparison and the “guilt by association” approach, we identified such systems in 141 archaeal genomes. These toxins can be classified into four major groups based on the structure of the components involved in the toxin delivery. The toxin domains are often shared between and within each system. We revisit halocin families and substantially expand the halocin C8 family, which was identified in diverse archaeal genomes and also certain bacteria. Finally, we employ features of protein sequences and genomic locus organization characteristic of archaeocins and polymorphic toxins to identify candidates for analogous but not necessarily homologous systems among uncharacterized protein families. This work confidently predicts that more than 1,600 archaeal proteins, currently annotated as “hypothetical” in public databases, are components of conflict and self-versus-nonself discrimination systems.

Unexpected Evolution of Lesion-Recognition Modules in Eukaryotic NER and Kinetoplast DNA Dynamics Proteins from Bacterial Mobile Elements

The provenance of several components of major uniquely eukaryotic molecular machines are increasingly being traced back to prokaryotic biological conflict systems. Here, we demonstrate that the N-terminal single-stranded DNA-binding domain from the anti-restriction protein ArdC, deployed by bacterial mobile elements against their host, was independently acquired twice by eukaryotes, giving rise to the DNA-binding domains of XPC/Rad4 and the Tc-38-like proteins in the stem kinetoplastid. In both instances, the ArdC-N domain tandemly duplicated forming an extensive DNA-binding interface. In XPC/Rad4, the ArdC-N domains (BHDs) also fused to the inactive transglutaminase domain of a peptide-N-glycanase ultimately derived from an archaeal conflict system. Alongside, we delineate several parallel acquisitions from conjugative elements/bacteriophages that gave rise to key components of the kinetoplast DNA (kDNA) replication apparatus. These findings resolve two outstanding questions in eukaryote biology: (1) the origin of the unique DNA lesion-recognition component of NER and (2) origin of the unusual, plasmid-like features of kDNA.

Evolutionary convergence and divergence in archaeal chromosomal proteins and Chromo-like domains from bacteria and eukaryotes

SH3-fold-β-barrel domains of the chromo-like superfamily recognize epigenetic marks in eukaryotic proteins. Their provenance has been placed either in archaea, based on apparent structural similarity to chromatin-compacting Sul7d and Cren7 proteins, or in bacteria based on the presence of sequence homologs. Using sequence and structural evidence we establish that the archaeal Cren7/Sul7 proteins emerged from a zinc ribbon (ZnR) ancestor. Further, we show that the ancestral eukaryotic chromo-like domains evolved from bacterial versions, likely acquired from early endosymbioses, which already possessed an aromatic cage for recognition of modified amino-groups. These bacterial versions are part of a radiation of secreted SH3-fold domains, which spawned both chromo-like domains and classical SH3 domains in the context of peptide-recognition in the peptidoglycan or the extracellular matrix. This establishes that Cren7/Sul7 converged to a “SH3”-like state from a ZnR precursor via the loss of metal-chelation and acquisition of stronger hydrophobic interactions; it is unlikely to have participated in the evolution of the chromo-like domains. We show that archaea possess several Cren7/Sul7-related proteins with intact Zn-chelating ligands, which we predict to play previously unstudied roles in chromosome segregation during cell-division comparable to the PRC barrel and CdvA domain proteins.

The ColM Family, Polymorphic Toxins Breaching the Bacterial Cell Wall

Bacteria host an arsenal of antagonism-mediating molecules to combat for ecologic space. Bacteriocins represent a pivotal group of secreted antibacterial peptides and proteins assisting in this fight, mainly eliminating relatives. Colicin M, a model for peptidoglycan-interfering bacteriocins in Gram-negative bacteria, appears to be part of a set of polymorphic toxins equipped with such a catalytic domain (ColM) targeting lipid II. Diversifying recombination has enabled parasitism of different receptors and has also given rise to hybrid bacteriocins in which ColM is associated with another toxin module. Remarkably, ColM toxins have recruited a diverse array of immunity partners, comprising cytoplasmic membrane-associated proteins with different topologies. Together, these findings suggest that different immunity mechanisms have evolved for ColM, in contrast to bacteriocins with nuclease activities.

Polyvalent Proteins, a Pervasive Theme in the Intergenomic Biological Conflicts of Bacteriophages and Conjugative Elements

Intense biological conflicts between prokaryotic genomes and their genomic parasites have resulted in an arms race in terms of the molecular “weaponry” deployed on both sides. Using a recursive computational approach, we uncovered a remarkable class of multidomain proteins with 2 to 15 domains in the same polypeptide deployed by viruses and plasmids in such conflicts. Domain architectures and genomic contexts indicate that they are part of a widespread conflict strategy involving proteins injected into the host cell along with parasite DNA during the earliest phase of infection. Their unique feature is the combination of domains with highly disparate biochemical activities in the same polypeptide; accordingly, we term them polyvalent proteins. Of the 131 domains in polyvalent proteins, a large fraction are enzymatic domains predicted to modify proteins, target nucleic acids, alter nucleotide signaling/metabolism, and attack peptidoglycan or cytoskeletal components. They further contain nucleic acid-binding domains, virion structural domains, and 40 novel uncharacterized domains. Analysis of their architectural network reveals both pervasive common themes and specialized strategies for conjugative elements and plasmids or (pro)phages. The themes include likely processing of multidomain polypeptides by zincin-like metallopeptidases and mechanisms to counter restriction or CRISPR/Cas systems and jump-start transcription or replication. DNA-binding domains acquired by eukaryotes from such systems have been reused in XPC/RAD4-dependent DNA repair and mitochondrial genome replication in kinetoplastids. Characterization of the novel domains discovered here, such as RNases and peptidases, are likely to aid in the development of new reagents and elucidation of the spread of antibiotic resistance.

IMPORTANCE This is the first report of the widespread presence of large proteins, termed polyvalent proteins, predicted to be transmitted by genomic parasites such as conjugative elements, plasmids, and phages during the initial phase of infection along with their DNA. They are typified by the presence of multiple domains with disparate activities combined in the same protein. While some of these domains are predicted to assist the invasive element in replication, transcription, or protection of their DNA, several are likely to target various host defense systems or modify the host to favor the parasite's life cycle. Notably, DNA-binding domains from these systems have been transferred to eukaryotes, where they have been incorporated into DNA repair and mitochondrial genome replication systems.

RNA damage in biological conflicts and the diversity of responding RNA repair systems

RNA is targeted in biological conflicts by enzymatic toxins or effectors. A vast diversity of systems which repair or ‘heal’ this damage has only recently become apparent. Here, we summarize the known effectors, their modes of action, and RNA targets before surveying the diverse systems which counter this damage from a comparative genomics viewpoint. RNA-repair systems show a modular organization with extensive shuffling and displacement of the constituent domains; however, a general ‘syntax’ is strongly maintained whereby systems typically contain: a RNA ligase (either ATP-grasp or RtcB superfamilies), nucleotidyltransferases, enzymes modifying RNA-termini for ligation (phosphatases and kinases) or protection (methylases), and scaffold or cofactor proteins. We highlight poorly-understood or previously-uncharacterized repair systems and components, e.g. potential scaffolding cofactors (Rot/TROVE and SPFH/Band-7 modules) with their respective cognate non-coding RNAs (YRNAs and a novel tRNA-like molecule) and a novel nucleotidyltransferase associating with diverse ligases. These systems have been extensively disseminated by lateral transfer between distant prokaryotic and microbial eukaryotic lineages consistent with intense inter-organismal conflict. Components have also often been ‘institutionalized’ for non-conflict roles, e.g. in RNA-splicing and in RNAi systems (e.g. in kinetoplastids) which combine a distinct family of RNA-acting prim-pol domains with DICER-like proteins.

Antibiotic Binding Drives Catalytic Activation of Aminoglycoside Kinase APH(2″)-Ia

APH(2″)-Ia is a widely disseminated resistance factor frequently found in clinical isolates of Staphylococcus aureus and pathogenic enterococci, where it is constitutively expressed. APH(2″)-Ia confers high-level resistance to gentamicin and related aminoglycosides through phosphorylation of the antibiotic using guanosine triphosphate (GTP) as phosphate donor. We have determined crystal structures of the APH(2″)-Ia in complex with GTP analogs, guanosine diphosphate, and aminoglycosides. These structures collectively demonstrate that aminoglycoside binding to the GTP-bound kinase drives conformational changes that bring distant regions of the protein into contact. These changes in turn drive a switch of the triphosphate cofactor from an inactive, stabilized conformation to a catalytically competent active conformation. This switch has not been previously reported for antibiotic kinases or for the structurally related eukaryotic protein kinases. This catalytic triphosphate switch presents a means by which the enzyme can curtail wasteful hydrolysis of GTP in the absence of aminoglycosides, providing an evolutionary advantage to this enzyme.

Transposons to toxins: the provenance, architecture and diversification of a widespread class of eukaryotic effectors

Enzymatic effectors targeting nucleic acids, proteins and other cellular components are the mainstay of conflicts across life forms. Using comparative genomics we identify a large class of eukaryotic proteins, which include effectors from oomycetes, fungi and other parasites. The majority of these proteins have a characteristic domain architecture with one of several N-terminal 'Header' domains, which are predicted to play a role in trafficking of these effectors, including a novel version of the Ubiquitin fold. The Headers are followed by one or more diverse C-terminal domains, such as restriction endonuclease (REase), protein kinase, HNH endonuclease, LK-nuclease (a RNase) and multiple distinct peptidase domains, which are predicted to carry their toxicity determinants. The most common types of these proteins appear to have originated from prokaryotic transposases (e.g. TN7 and Mu) and combine a CDC6/ORC1-STAND clade NTPase domain with a C-terminal REase domain. Other than the so-called Crinkler effectors of oomycetes and fungi, these effectors are encoded by other eukaryotic parasites such as trypanosomatids (the RHS proteins) and the rhizarian Plasmodiophora, and symbionts like Capsaspora Remarkably, we also find these proteins in free-living eukaryotes, including several viridiplantae, fungi, amoebozoans and animals. These versions might either still be transposons or function in other poorly understood eukaryote-specific inter-organismal and inter-genomic conflicts. These include the Medea1 selfish element of Tribolium that spreads via post-zygotic killing. We present a unified mechanism for the recombination-dependent diversification and action of this widespread class of molecular weaponry deployed across diverse conflicts ranging from parasitic to free-living forms.

Comparative genomic analyses reveal a vast, novel network of nucleotide-centric systems in biological conflicts, immunity and signaling

Cyclic di- and linear oligo-nucleotide signals activate defenses against invasive nucleic acids in animal immunity; however, their evolutionary antecedents are poorly understood. Using comparative genomics, sequence and structure analysis, we uncovered a vast network of systems defined by conserved prokaryotic gene-neighborhoods, which encode enzymes generating such nucleotides or alternatively processing them to yield potential signaling molecules. The nucleotide-generating enzymes include several clades of the DNA-polymerase β-like superfamily (including Vibrio cholerae DncV), a minimal version of the CRISPR polymerase and DisA-like cyclic-di-AMP synthetases. Nucleotide-binding/processing domains include TIR domains and members of a superfamily prototyped by Smf/DprA proteins and base (cytokinin)-releasing LOG enzymes. They are combined in conserved gene-neighborhoods with genes for a plethora of protein superfamilies, which we predict to function as nucleotide-sensors and effectors targeting nucleic acids, proteins or membranes (pore-forming agents). These systems are sometimes combined with other biological conflict-systems such as restriction-modification and CRISPR/Cas. Interestingly, several are coupled in mutually exclusive neighborhoods with either a prokaryotic ubiquitin-system or a HORMA domain-PCH2-like AAA+ ATPase dyad. The latter are potential precursors of equivalent proteins in eukaryotic chromosome dynamics. Further, components from these nucleotide-centric systems have been utilized in several other systems including a novel diversity-generating system with a reverse transcriptase. We also found the Smf/DprA/LOG domain from these systems to be recruited as a predicted nucleotide-binding domain in eukaryotic TRPM channels. These findings point to evolutionary and mechanistic links, which bring together CRISPR/Cas, animal interferon-induced immunity, and several other systems that combine nucleic-acid-sensing and nucleotide-dependent signaling.

Two independently folding units of Plasmodium profilin suggest evolution via gene fusion

Gene fusion is a common mechanism of protein evolution that has mainly been discussed in the context of multidomain or symmetric proteins. Less is known about fusion of ancestral genes to produce small single-domain proteins. Here, we show with a domain-swapped mutant Plasmodium profilin that this small, globular, apparently single-domain protein consists of two foldons. The separation of binding sites for different protein ligands in the two halves suggests evolution via an ancient gene fusion event, analogous to the formation of multidomain proteins. Finally, the two fragments can be assembled together after expression as two separate gene products. The possibility to engineer both domain-swapped dimers and half-profilins that can be assembled back to a full profilin provides perspectives for engineering of novel protein folds, e.g., with different scaffolding functions.

An interbacterial NAD(P)(+) glycohydrolase toxin requires elongation factor Tu for delivery to target cells

Type VI secretion (T6S) influences the composition of microbial communities by catalyzing the delivery of toxins between adjacent bacterial cells. Here, we demonstrate that a T6S integral membrane toxin from Pseudomonas aeruginosa, Tse6, acts on target cells by degrading the universally essential dinucleotides NAD(+) and NADP(+). Structural analyses of Tse6 show that it resembles mono-ADP-ribosyltransferase proteins, such as diphtheria toxin, with the exception of a unique loop that both excludes proteinaceous ADP-ribose acceptors and contributes to hydrolysis. We find that entry of Tse6 into target cells requires its binding to an essential housekeeping protein, translation elongation factor Tu (EF-Tu). These proteins participate in a larger assembly that additionally directs toxin export and provides chaperone activity. Visualization of this complex by electron microscopy defines the architecture of a toxin-loaded T6S apparatus and provides mechanistic insight into intercellular membrane protein delivery between bacteria.

Contact-Dependent Growth Inhibition (CDI) and CdiB/CdiA Two-Partner Secretion Proteins

Bacteria have developed several strategies to communicate and compete with one another in complex environments. One important mechanism of inter-bacterial competition is contact-dependent growth inhibition (CDI), in which Gram-negative bacteria use CdiB/CdiA two-partner secretion proteins to suppress the growth of neighboring target cells. CdiB is an Omp85 outer-membrane protein that exports and assembles CdiA exoproteins onto the inhibitor cell surface. CdiA binds to receptors on susceptible bacteria and subsequently delivers its C-terminal toxin domain (CdiA-CT) into the target cell. CDI systems also encode CdiI immunity proteins, which specifically bind to the CdiA-CT and neutralize its toxin activity, thereby protecting CDI(+) cells from auto-inhibition. Remarkably, CdiA-CT sequences are highly variable between bacteria, as are the corresponding CdiI immunity proteins. Variations in CDI toxin/immunity proteins suggest that these systems function in bacterial self/non-self recognition and thereby play an important role in microbial communities. In this review, we discuss recent advances in the biochemistry, structural biology and physiology of CDI.

The eukaryotic translation initiation regulator CDC123 defines a divergent clade of ATP-grasp enzymes with a predicted role in novel protein modifications

Abstract

Deciphering the origin of uniquely eukaryotic features of sub-cellular systems, such as the translation apparatus, is critical in reconstructing eukaryogenesis. One such feature is the highly conserved, but poorly understood, eukaryotic protein CDC123, which regulates the abundance of the eukaryotic translation initiation eIF2 complex and binds one of its components eIF2γ. We show that the eukaryotic protein CDC123 defines a novel clade of ATP-grasp enzymes distinguished from all other members of the superfamily by a RAGNYA domain with two conserved lysines (henceforth the R2K clade). Combining the available biochemical and genetic data on CDC123 with the inferred enzymatic function, we propose that the eukaryotic CDC123 proteins are likely to function as ATP-dependent protein-peptide ligases which modify proteins by ribosome-independent addition of an oligopeptide tag. We also show that the CDC123 family emerged first in bacteria where it appears to have diversified along with the two other families of the R2K clade. The bacterial CDC123 family members are of two distinct types, one found as part of type VI secretion systems which deliver polymorphic toxins and the other functioning as potential effectors delivered to amoeboid eukaryotic hosts. Representatives of the latter type have also been independently transferred to phylogenetically unrelated amoeboid eukaryotes and their nucleo-cytoplasmic large DNA viruses. Similarly, the two other prokaryotic R2K clade families are also proposed to participate in biological conflicts between bacteriophages and their hosts. These findings add further evidence to the recently proposed hypothesis that the horizontal transfer of enzymatic effectors from the bacterial endosymbionts of the stem eukaryotes played a fundamental role in the emergence of the characteristically eukaryotic regulatory systems and sub-cellular structures.

Reviewers

This article was reviewed by Michael Galperin and Sandor Pongor.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0053-x) contains supplementary material, which is available to authorized users.

New players in the toxin field: polymorphic toxin systems in bacteria

Bacteria have evolved numerous strategies to increase their competitiveness and fight against each other. Indeed, a large arsenal of antibacterial weapons is available in order to inhibit the proliferation of competitor cells. Polymorphic toxin systems (PTS), recently identified by bioinformatics in all major bacterial lineages, correspond to such a system primarily involved in conflict between related bacterial strains. They are typically composed of a secreted multidomain toxin, a protective immunity protein, and multiple cassettes encoding alternative toxic domains. The C-terminal domains of polymorphic toxins carry the toxic activity, whereas the N-terminal domains are related to the trafficking mode. In silico analysis of PTS identified over 150 distinct toxin domains, including putative nuclease, deaminase, or peptidase domains. Immunity genes found immediately downstream of the toxin genes encode small proteins that protect bacteria against their own toxins or against toxins secreted by neighboring cells. PTS encompass well-known colicins and pyocins, contact-dependent growth inhibition systems which include CdiA and Rhs toxins and some effectors of type VI secretion systems. We have recently characterized the MafB toxins, a new family of PTS deployed by pathogenic Neisseria spp. Many other putative PTS have been identified by in silico predictions but have yet to be characterized experimentally. However, the high number of these systems suggests that PTS have a fundamental role in bacterial biology that is likely to extend beyond interbacterial competition.

One cannot rule them all: Are bacterial toxins-antitoxins druggable?

Type II (proteic) toxin–antitoxin (TA) operons are widely spread in bacteria and archaea. They are organized as operons in which, usually, the antitoxin gene precedes the cognate toxin gene. The antitoxin generally acts as a transcriptional self-repressor, whereas the toxin acts as a co-repressor, both proteins constituting a harmless complex. When bacteria encounter a stressful environment, TAs are triggered. The antitoxin protein is unstable and will be degraded by host proteases, releasing the free toxin to halt essential processes. The result is a cessation of cell growth or even death. Because of their ubiquity and the essential processes targeted, TAs have been proposed as good candidates for development of novel antimicrobials. We discuss here the possible druggability of TAs as antivirals and antibacterials, with focus on the potentials and the challenges that their use may find in the ‘real’ world. We present strategies to develop TAs as antibacterials in view of novel technologies, such as the use of very small molecules (fragments) as inhibitors of protein–protein interactions. Appropriate fragments could disrupt the T:A interfaces leading to the release of the targeted TA pair. Possible ways of delivery and formulation of Tas are also discussed.

We consider various approaches to develop the toxins of the type II family as possible candidates to drug discovery; druggability of toxins-antitoxins could be possible as antivirals. As antibacterials, they might be considered as druggable but delivery and formulation may not be simple so far.

The enzymatic nature of an anonymous protein sequence cannot reliably be inferred from superfamily level structural information alone

As the largest fraction of any proteome does not carry out enzymatic functions, and in order to leverage 3D structural data for the annotation of increasingly higher volumes of sequence data, we wanted to assess the strength of the link between coarse grained structural data (i.e., homologous superfamily level) and the enzymatic versus non-enzymatic nature of protein sequences. To probe this relationship, we took advantage of 41 phylogenetically diverse (encompassing 11 distinct phyla) genomes recently sequenced within the GEBA initiative, for which we integrated structural information, as defined by CATH, with enzyme level information, as defined by Enzyme Commission (EC) numbers. This analysis revealed that only a very small fraction (about 1%) of domain sequences occurring in the analyzed genomes was found to be associated with homologous superfamilies strongly indicative of enzymatic function. Resorting to less stringent criteria to define enzyme versus non-enzyme biased structural classes or excluding highly prevalent folds from the analysis had only modest effect on this proportion. Thus, the low genomic coverage by structurally anchored protein domains strongly associated to catalytic activities indicates that, on its own, the power of coarse grained structural information to infer the general property of being an enzyme is rather limited.