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/.

Under stressful conditions, bacterial RelA-SpoT Homolog (RSH) enzymes synthesize the alarmone (p)ppGpp, a nucleotide second messenger. (p)ppGpp rewires bacterial transcription and metabolism to cope with stress, and, at high concentrations, inhibits the process of protein synthesis and bacterial growth to save and redirect resources until conditions improve. Single-domain small alarmone synthetases (SASs) are RSH family members that contain the (p)ppGpp synthesis (SYNTH) domain, but lack the hydrolysis (HD) domain and regulatory C-terminal domains of the long RSHs such as Rel, RelA, and SpoT. We asked whether analysis of the genomic context of SASs can indicate possible functional roles. Indeed, multiple SAS subfamilies are encoded in widespread conserved bicistronic operon architectures that are reminiscent of those typically seen in toxin-antitoxin (TA) operons. We have validated five of these SASs as being toxic (toxSASs), with neutralization by the protein products of six neighboring antitoxin genes. The toxicity of Cellulomonas marina toxSAS FaRel is mediated by the accumulation of alarmones ppGpp and ppApp, and an associated depletion of cellular guanosine triphosphate and adenosine triphosphate pools, and is counteracted by its HD domain-containing antitoxin. Thus, the ToxSAS-antiToxSAS system with its multiple different antitoxins exemplifies how ancient nucleotide-based signaling mechanisms can be repurposed as TA modules during evolution, potentially multiple times independently.

PoxtA and OptrA are ATP binding cassette (ABC) proteins of the F subtype (ABCF). They confer resistance to oxazolidinone and phenicol antibiotics, such as linezolid and chloramphenicol, which stall translating ribosomes when certain amino acids are present at a defined position in the nascent polypeptide chain. These proteins are often encoded on mobile genetic elements, facilitating their rapid spread amongst Gram-positive bacteria, and are thought to confer resistance by binding to the ribosome and dislodging the bound antibiotic. However, the mechanistic basis of this resistance remains unclear. Here we refine the PoxtA spectrum of action, demonstrate alleviation of linezolid-induced context-dependent translational stalling, and present cryo-electron microscopy structures of PoxtA in complex with the Enterococcus faecalis 70S ribosome. PoxtA perturbs the CCA-end of the P-site tRNA, causing it to shift by ∼4 Å out of the ribosome, corresponding to a register shift of approximately one amino acid for an attached nascent polypeptide chain. We postulate that the perturbation of the P-site tRNA by PoxtA thereby alters the conformation of the attached nascent chain to disrupt the drug binding site.

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.

Many Gram-positive pathogenic bacteria employ ribosomal protection proteins (RPPs) to confer resistance to clinically important antibiotics. In Bacillus subtilis, the RPP VmlR confers resistance to lincomycin (Lnc) and the streptogramin A (SA) antibiotic virginiamycin M (VgM). VmlR is an ATP-binding cassette (ABC) protein of the F type, which, like other antibiotic resistance (ARE) ABCF proteins, is thought to bind to antibiotic-stalled ribosomes and promote dissociation of the drug from its binding site. To investigate the molecular mechanism by which VmlR confers antibiotic resistance, we have determined a cryo-electron microscopy (cryo-EM) structure of an ATPase-deficient B. subtilis VmlR-EQ2 mutant in complex with a B. subtilis ErmDL-stalled ribosomal complex (SRC). The structure reveals that VmlR binds within the E site of the ribosome, with the antibiotic resistance domain (ARD) reaching into the peptidyltransferase center (PTC) of the ribosome and a C-terminal extension (CTE) making contact with the small subunit (SSU). To access the PTC, VmlR induces a conformational change in the P-site tRNA, shifting the acceptor arm out of the PTC and relocating the CCA end of the P-site tRNA toward the A site. Together with microbiological analyses, our study indicates that VmlR allosterically dissociates the drug from its ribosomal binding site and exhibits specificity to dislodge VgM, Lnc, and the pleuromutilin tiamulin (Tia), but not chloramphenicol (Cam), linezolid (Lnz), nor the macrolide erythromycin (Ery).

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.