Uncovering new families and folds in the natural protein universe.

We are now entering a new era in protein sequence and structure annotation, with hundreds of millions of predicted protein structures made available through the AlphaFold database1. These models cover nearly all proteins that are known, including those challenging to annotate for function or putative biological role using standard homology-based approaches. In this study, we examine the extent to which the AlphaFold database has structurally illuminated this 'dark matter' of the natural protein universe at high predicted accuracy. We further describe the protein diversity that these models cover as an annotated interactive sequence similarity network, accessible at https://uniprot3d.org/atlas/AFDB90v4 . By searching for novelties from sequence, structure and semantic perspectives, we uncovered the ß-flower fold, added several protein families to Pfam database2 and experimentally demonstrated that one of these belongs to a new superfamily of translation-targeting toxin-antitoxin systems, TumE-TumA. This work underscores the value of large-scale efforts in identifying, annotating and prioritizing new protein families. By leveraging the recent deep learning revolution in protein bioinformatics, we can now shed light into uncharted areas of the protein universe at an unprecedented scale, paving the way to innovations in life sciences and biotechnology.

The structural basis of hyperpromiscuity in a core combinatorial network of type II toxin-antitoxin and related phage defense systems.

Toxin-antitoxin (TA) systems are a large group of small genetic modules found in prokaryotes and their mobile genetic elements. Type II TAs are encoded as bicistronic (two-gene) operons that encode two proteins: a toxin and a neutralizing antitoxin. Using our tool NetFlax (standing for Network-FlaGs for toxins and antitoxins), we have performed a large-scale bioinformatic analysis of proteinaceous TAs, revealing interconnected clusters constituting a core network of TA-like gene pairs. To understand the structural basis of toxin neutralization by antitoxins, we have predicted the structures of 3,419 complexes with AlphaFold2. Together with mutagenesis and functional assays, our structural predictions provide insights into the neutralizing mechanism of the hyperpromiscuous Panacea antitoxin domain. In antitoxins composed of standalone Panacea, the domain mediates direct toxin neutralization, while in multidomain antitoxins the neutralization is mediated by other domains, such as PAD1, Phd-C, and ZFD. We hypothesize that Panacea acts as a sensor that regulates TA activation. We have experimentally validated 16 NetFlax TA systems and used domain annotations and metabolic labeling assays to predict their potential mechanisms of toxicity (such as membrane disruption, and inhibition of cell division or protein synthesis) as well as biological functions (such as antiphage defense). We have validated the antiphage activity of a RosmerTA system encoded by Gordonia phage Kita, and used fluorescence microscopy to confirm its predicted membrane-depolarizing activity. The interactive version of the NetFlax TA network that includes structural predictions can be accessed at http://netflax.webflags.se/.

A hyperpromiscuous antitoxin protein domain for the neutralization of diverse toxin domains.

Toxin-antitoxin (TA) gene pairs are ubiquitous in microbial chromosomal genomes and plasmids as well as temperate bacteriophages. They act as regulatory switches, with the toxin limiting the growth of bacteria and archaea by compromising diverse essential cellular targets and the antitoxin counteracting the toxic effect. To uncover previously uncharted TA diversity across microbes and bacteriophages, we analyzed the conservation of genomic neighborhoods using our computational tool FlaGs (for flanking genes), which allows high-throughput detection of TA-like operons. Focusing on the widespread but poorly experimentally characterized antitoxin domain DUF4065, our in silico analyses indicated that DUF4065-containing proteins serve as broadly distributed antitoxin components in putative TA-like operons with dozens of different toxic domains with multiple different folds. Given the versatility of DUF4065, we have named the domain Panacea (and proteins containing the domain, PanA) after the Greek goddess of universal remedy. We have experimentally validated nine PanA-neutralized TA pairs. While the majority of validated PanA-neutralized toxins act as translation inhibitors or membrane disruptors, a putative nucleotide cyclase toxin from a Burkholderia prophage compromises transcription and translation as well as inducing RelA-dependent accumulation of the nucleotide alarmone (p)ppGpp. We find that Panacea-containing antitoxins form a complex with their diverse cognate toxins, characteristic of the direct neutralization mechanisms employed by Type II TA systems. Finally, through directed evolution, we have selected PanA variants that can neutralize noncognate TA toxins, thus experimentally demonstrating the evolutionary plasticity of this hyperpromiscuous antitoxin domain.

Toxins MazF and MqsR cleave Escherichia coli rRNA precursors at multiple sites.

The endoribonuclease toxins of the E. coli toxin-antitoxin systems arrest bacterial growth and protein synthesis by targeting cellular mRNAs. As an exception, E. coli MazF was reported to cleave also 16S rRNA at a single site and separate an anti-Shine-Dalgarno sequence-containing RNA fragment from the ribosome. We noticed extensive rRNA fragmentation in response to induction of the toxins MazF and MqsR, which suggested that these toxins can cleave rRNA at multiple sites. We adapted differential RNA-sequencing to map the toxin-cleaved 5'- and 3'-ends. Our results show that the MazF and MqsR cleavage sites are located within structured rRNA regions and, therefore, are not accessible in assembled ribosomes. Most of the rRNA fragments are located in the aberrant ribosomal subunits that accumulate in response to toxin induction and contain unprocessed rRNA precursors. We did not detect MazF- or MqsR-cleaved rRNA in stationary phase bacteria and in assembled ribosomes. Thus, we conclude that MazF and MqsR cleave rRNA precursors before the ribosomes are assembled and potentially facilitate the decay of surplus rRNA transcripts during stress.

Mechanism of phage sensing and restriction by toxin-antitoxin-chaperone systems.

Toxin-antitoxins (TAs) are prokaryotic two-gene systems composed of a toxin neutralized by an antitoxin. Toxin-antitoxin-chaperone (TAC) systems additionally include a SecB-like chaperone that stabilizes the antitoxin by recognizing its chaperone addiction (ChAD) element. TACs mediate antiphage defense, but the mechanisms of viral sensing and restriction are unexplored. We identify two Escherichia coli antiphage TAC systems containing host inhibition of growth (HigBA) and CmdTA TA modules, HigBAC and CmdTAC. HigBAC is triggered through recognition of the gpV major tail protein of phage λ. Chaperone HigC recognizes gpV and ChAD via analogous aromatic molecular patterns, with gpV outcompeting ChAD to trigger toxicity. For CmdTAC, the CmdT ADP-ribosyltransferase toxin modifies mRNA to halt protein synthesis and limit phage propagation. Finally, we establish the modularity of TACs by creating a hybrid broad-spectrum antiphage system combining the CmdTA TA warhead with a HigC chaperone phage sensor. Collectively, these findings reveal the potential of TAC systems in broad-spectrum antiphage defense.

Transcriptional cross-activation between toxin-antitoxin systems of Escherichia coli.

BACKGROUND: Bacterial toxin-antitoxin (TA) systems are formed by potent regulatory or suicide factors (toxins) and their short-lived inhibitors (antitoxins). Antitoxins are DNA-binding proteins and auto-repress transcription of TA operons. Transcription of multiple TA operons is activated in temporarily non-growing persister cells that can resist killing by antibiotics. Consequently, the antitoxin levels of persisters must have been dropped and toxins are released of inhibition. RESULTS: Here, we describe transcriptional cross-activation between different TA systems of Escherichia coli. We find that the chromosomal relBEF operon is activated in response to production of the toxins MazF, MqsR, HicA, and HipA. Expression of the RelE toxin in turn induces transcription of several TA operons. We show that induction of mazEF during amino acid starvation depends on relBE and does not occur in a relBEF deletion mutant. Induction of TA operons has been previously shown to depend on Lon protease which is activated by polyphospate accumulation. We show that transcriptional cross-activation occurs also in strains deficient for Lon, ClpP, and HslV proteases and polyphosphate kinase. Furthermore, we find that toxins cleave the TA mRNA in vivo, which is followed by degradation of the antitoxin-encoding fragments and selective accumulation of the toxin-encoding regions. We show that these accumulating fragments can be translated to produce more toxin. CONCLUSION: Transcriptional activation followed by cleavage of the mRNA and disproportionate production of the toxin constitutes a possible positive feedback loop, which can fire other TA systems and cause bistable growth heterogeneity. Cross-interacting TA systems have a potential to form a complex network of mutually activating regulators in bacteria.

Fragmentation of Escherichia coli mRNA by MazF and MqsR.

MazEF and MqsRA are toxin-antitoxin systems, where the toxins MazF and MqsR sequence-specifically cleave single-stranded RNA, thereby shutting down protein synthesis and cell growth. However, it has been proposed that MazF functions in a highly specific pathway, where it truncates the 5' ends of a set of E. coli transcripts (the MazF regulon), which are then translated under stress conditions by specialized ribosomes. We mapped the cleavage sites of MazF and MqsR throughout the E. coli transcriptome. Our results show that both toxins cleave mRNA independently of the recognition site position and MazF freely cleaves transcripts of the proposed MazF regulon within coding sequences. Proteome analysis indicated that MazF expression leads to overall inhibition of protein synthesis and the putative MazF regulon proteins are not selectively synthesized in response to the toxin. Our results support a simpler role for endoribonuclease TA systems as indifferent destroyers of unstructured RNA.