Characterization of Escherichia coli dinJ-yafQ toxin-antitoxin system using insights from mutagenesis data

Mono- and multidomain defense toxins of the RelE/ParE superfamily

Toxin-antitoxin (TA) modules are widely distributed across prokaryotes, often existing in large numbers despite their associated fitness costs. Most type II TA modules are bicistronic operons encoding a monodomain toxin and its cognate protein antitoxin. The RelE/ParE superfamily encompasses toxins with a conserved Barnase-EndoU-ColicinE5/D-RelE (BECR) fold. Yet, their cellular targets differ remarkably: RelE toxins function as ribosome-dependent RNases, while ParE toxins act as DNA gyrase inhibitors. Using a comprehensive bioinformatics approach, this study analyzed 13 BECR-fold toxin families as classified in the Pfam database. Intriguingly, the ParE family was found to include a subcluster of mRNA-cleaving toxins, challenging its conventional role as solely DNA-targeting. This study identified a novel tripartite operon encoding a PtuA-like defense ATPase, a homolog of type IV restriction endonucleases, and a RelE homolog, suggesting a coordinated role in defense mechanisms. Multidomain BECR-fold toxins, including transmembrane variants, were also discovered, extending the functional repertoire of type II TA modules to membrane-associated systems. These findings clarify the evolutionary relationships and functional diversity within the RelE/ParE superfamily and discover novel, putative defense systems that can now be investigated experimentally.IMPORTANCEToxin-antitoxin modules play critical roles in prokaryotic survival and adaptation, contributing to genome stabilization and defense against phages and invading plasmids. The RelE/ParE superfamily exemplifies the structural and functional diversity of these systems, with members targeting distinct cellular processes, such as translation and DNA supercoiling. By elucidating the relationships among the 13 BECR-fold toxin families, this study enhances our understanding of microbial resistance mechanisms and reveals potential new opportunities for research into prokaryotic defense and regulation. These insights may have significant implications for medical and biotechnological applications, particularly in understanding bacterial responses to genetic invaders.

SProtFP: a machine learning-based method for functional classification of small ORFs in prokaryotes

Small proteins (≤100 amino acids) play important roles across all life forms, ranging from unicellular bacteria to higher organisms. In this study, we have developed SProtFP which is a machine learning-based method for functional annotation of prokaryotic small proteins into selected functional categories. SProtFP uses independent artificial neural networks (ANNs) trained using a combination of physicochemical descriptors for classifying small proteins into antitoxin type 2, bacteriocin, DNA-binding, metal-binding, ribosomal protein, RNA-binding, type 1 toxin and type 2 toxin proteins. We have also trained a model for identification of small open reading frame (smORF)-encoded antimicrobial peptides (AMPs). Comprehensive benchmarking of SProtFP revealed an average area under the receiver operator curve (ROC-AUC) of 0.92 during 10-fold cross-validation and an ROC-AUC of 0.94 and 0.93 on held-out balanced and imbalanced test sets. Utilizing our method to annotate bacterial isolates from the human gut microbiome, we could identify thousands of remote homologs of known small protein families and assign putative functions to uncharacterized proteins. This highlights the utility of SProtFP for large-scale functional annotation of microbiome datasets, especially in cases where sequence homology is low. SProtFP is freely available at http://www.nii.ac.in/sprotfp.html and can be combined with genome annotation tools such as ProsmORF-pred to uncover the functional repertoire of novel small proteins in bacteria.

Anaerobiosis, a neglected factor in phage-bacteria interactions

IMPORTANCE

Many parameters affect phage-bacteria interaction. Some of these parameters depend on the environment in which the bacteria are present. Anaerobiosis effect on phage infection in facultative anaerobic bacteria has not yet been studied. The absence of oxygen triggers metabolic changes in facultative bacteria and this affects phage infection and viral life cycle. Understanding how an anaerobic environment can alter the behavior of phages during infection is relevant for the phage therapy success.

Functional characterization and transcriptional repression by Lacticaseibacillus paracasei DinJ-YafQ

Abstract

DinJ-YafQ is a bacterial type II TA system formed by the toxin RNase YafQ and the antitoxin protein DinJ. The activity of YafQ and DinJ has been rigorously studied in Escherichia coli, but little has been reported about orthologous systems identified in different microorganisms. In this work, we report an in vitro and in vivo functional characterization of YafQ and DinJ identified in two different strains of Lacticaseibacillus paracasei and isolated as recombinant proteins. While DinJ is identical in both strains, the two YafQ orthologs differ only for the D72G substitution in the catalytic site. Both YafQ orthologs digest ribosomal RNA, albeit with different catalytic efficiencies, and their RNase activity is neutralized by DinJ. We further show that DinJ alone or in complex with YafQ can bind cooperatively to a 28-nt inverted repeat overlapping the −35 element of the TA operon promoter. Atomic force microscopy imaging of DinJ-YafQ in complex with DNA harboring the cognate site reveals the formation of different oligomeric states that prevent the binding of RNA polymerase to the promoter. A single amino acid substitution (R13A) within the RHH DNA-binding motif of DinJ is sufficient to abolish DinJ and DinJ-YafQ DNA binding in vitro. In vivo experiments confirm the negative regulation of the TA promoter by DinJ and DinJ-YafQ and unveil an unexpected high expression-related toxicity of the gfp reporter gene. A model for the binding of two YafQ-(DinJ)₂-YafQ tetramers to the promoter inverted repeat showing the absence of protein-protein steric clash is also presented.

Key points

• The RNase activity of L. paracasei YafQ toxin is neutralized by DinJ antitoxin.

• DinJ and DinJ-YafQ bind to an inverted repeat to repress their own promoter.

• The R13A mutation of DinJ abolishes DNA binding of both DinJ and DinJ-YafQ.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00253-022-12195-4.

Global Analysis of the Specificities and Targets of Endoribonucleases from Escherichia coli Toxin-Antitoxin Systems

Toxin-antitoxin systems are widely distributed genetic modules typically featuring toxins that can inhibit bacterial growth and antitoxins that can reverse inhibition. Although Escherichia coli encodes 11 toxins with known or putative endoribonuclease activity, the targets of most of these toxins remain poorly characterized. Using a new RNA sequencing (RNA-seq) pipeline that enables the mapping and quantification of RNA cleavage with single-nucleotide resolution, we characterized the targets and specificities of 9 endoribonuclease toxins from E. coli. We found that these toxins use low-information cleavage motifs to cut a significant proportion of mRNAs in E. coli, but not tRNAs or the rRNAs from mature ribosomes. However, all the toxins, including those that are ribosome dependent and cleave only translated RNA, inhibit ribosome biogenesis. This inhibition likely results from the cleavage of ribosomal protein transcripts, which disrupts the stoichiometry and biogenesis of new ribosomes and causes the accumulation of aberrant ribosome precursors. Collectively, our results provide a comprehensive, global analysis of endoribonuclease-based toxin-antitoxin systems in E. coli and support the conclusion that, despite their diversity, each disrupts translation and ribosome biogenesis.

Strategies to Investigate Membrane Damage, Nucleoid Condensation, and RNase Activity of Bacterial Toxin-Antitoxin Systems

A large number of bacterial toxin–antitoxin (TA) systems have been identified so far and different experimental approaches have been explored to investigate their activity and regulation both in vivo and in vitro. Nonetheless, a common feature of these methods is represented by the difficulty in cell transformation, culturing, and stability of the transformants, due to the expression of highly toxic proteins. Recently, in dealing with the type I Lpt/RNAII and the type II YafQ/DinJ TA systems, we encountered several of these problems that urged us to optimize methodological strategies to study the phenotype of recombinant Escherichia coli host cells. In particular, we have found conditions to tightly repress toxin expression by combining the pET expression system with the E. coli C41(DE3) pLysS strain. To monitor the RNase activity of the YafQ toxin, we developed a fluorescence approach based on Thioflavin-T which fluoresces brightly when complexed with bacterial RNA. Fluorescence microscopy was also applied to reveal loss of membrane integrity associated with the activity of the type I toxin Lpt, by using DAPI and ethidium bromide to selectively stain cells with impaired membrane permeability. We further found that atomic force microscopy can readily be employed to characterize toxin-induced membrane damages.

A minimal model for gene expression dynamics of bacterial type II toxin-antitoxin systems

Toxin-antitoxin (TA) modules are part of most bacteria's regulatory machinery for stress responses and general aspects of their physiology. Due to the interplay of a long-lived toxin with a short-lived antitoxin, TA modules have also become systems of interest for mathematical modelling. Here we resort to previous modelling efforts and extract from these a minimal model of type II TA system dynamics on a timescale of hours, which can be used to describe time courses derived from gene expression data of TA pairs. We show that this model provides a good quantitative description of TA dynamics for the 11 TA pairs under investigation here, while simpler models do not. Our study brings together aspects of Biophysics with its focus on mathematical modelling and Computational Systems Biology with its focus on the quantitative interpretation of 'omics' data. This mechanistic model serves as a generic transformation of time course information into kinetic parameters. The resulting parameter vector can, in turn, be mechanistically interpreted. We expect that TA pairs with similar mechanisms are characterized by similar vectors of kinetic parameters, allowing us to hypothesize on the mode of action for TA pairs still under discussion.

Comparative genome analysis of Salmonella enterica serovar Gallinarum biovars Pullorum and Gallinarum decodes strain specific genes

Salmonella enterica serovar Gallinarum biovar Pullorum (bvP) and biovar Gallinarum (bvG) are the etiological agents of pullorum disease (PD) and fowl typhoid (FT) respectively, which cause huge economic losses to poultry industry especially in developing countries including India. Vaccination and biosecurity measures are currently being employed to control and reduce the S. Gallinarum infections. High endemicity, poor implementation of hygiene and lack of effective vaccines pose challenges in prevention and control of disease in intensively maintained poultry flocks. Comparative genome analysis unravels similarities and dissimilarities thus facilitating identification of genomic features that aids in pathogenesis, niche adaptation and in tracing of evolutionary history. The present investigation was carried out to assess the genotypic differences amongst S.enterica serovar Gallinarum strains including Indian strain S. Gallinarum Sal40 VTCCBAA614. The comparative genome analysis revealed an open pan-genome consisting of 5091 coding sequence (CDS) with 3270 CDS belonging to core-genome, 1254 CDS to dispensable genome and strain specific genes i.e. singletons ranging from 3 to 102 amongst the analyzed strains. Moreover, the investigated strains exhibited diversity in genomic features such as virulence factors, genomic islands, prophage regions, toxin-antitoxin cassettes, and acquired antimicrobial resistance genes. Core genome identified in the study can give important leads in the direction of design of rapid and reliable diagnostics, and vaccine design for effective infection control as well as eradication. Additionally, the identified genetic differences among the S. enterica serovar Gallinarum strains could be used for bacterial typing, structure based inhibitor development by future experimental investigations on the data generated.

Comparative Genomics Suggests Mechanisms of Genetic Adaptation toward the Catabolism of the Phenylurea Herbicide Linuron in Variovorax

Biodegradation of the phenylurea herbicide linuron appears a specialization within a specific clade of the Variovorax genus. The linuron catabolic ability is likely acquired by horizontal gene transfer but the mechanisms involved are not known. The full-genome sequences of six linuron-degrading Variovorax strains isolated from geographically distant locations were analyzed to acquire insight into the mechanisms of genetic adaptation toward linuron metabolism. Whole-genome sequence analysis confirmed the phylogenetic position of the linuron degraders in a separate clade within Variovorax and indicated that they unlikely originate from a common ancestral linuron degrader. The linuron degraders differentiated from Variovorax strains that do not degrade linuron by the presence of multiple plasmids of 20–839 kb, including plasmids of unknown plasmid groups. The linuron catabolic gene clusters showed 1) high conservation and synteny and 2) strain-dependent distribution among the different plasmids. Most of them were bordered by IS1071 elements forming composite transposon structures, often in a multimeric array configuration, appointing IS1071 as a key element in the recruitment of linuron catabolic genes in Variovorax. Most of the strains carried at least one (catabolic) broad host range plasmid that might have been a second instrument for catabolic gene acquisition. We conclude that clade 1 Variovorax strains, despite their different geographical origin, made use of a limited genetic repertoire regarding both catabolic functions and vehicles to acquire linuron biodegradation.

12/111phiA Prophage Domestication Is Associated with Autoaggregation and Increased Ability to Produce Biofilm in Streptococcus agalactiae

CC17 Streptococcus agalactiae carrying group-A prophages is increasingly responsible for neonatal infections. To investigate the impact of the genetic features of a group-A prophage, we first conducted an in silico analysis of the genome of 12/111phiA, a group-A prophage carried by a strain responsible for a bloodstream infection in a parturient. This revealed a Restriction Modification system, suggesting a prophage maintenance strategy and five ORFs of interest for the host and encoding a type II toxin antitoxin system RelB/YafQ, an endonuclease, an S-adenosylmethionine synthetase MetK, and an StrP-like adhesin. Using the WT strain cured from 12/111phiA and constructing deleted mutants for the ORFs of interest, and their complemented mutants, we demonstrated an impact of prophage features on growth characteristics, cell morphology and biofilm formation. Our findings argue in favor of 12/111phiA domestication by the host and a role of prophage features in cell autoaggregation, glycocalyx and biofilm formation. We suggest that lysogeny may promote GBS adaptation to the acid environment of the vagina, consequently colonizing and infecting neonates.

An improved direct metamobilome approach increases the detection of larger-sized circular elements across kingdoms

Mobile genetic elements (MGEs) are instrumental in natural prokaryotic genome editing, permitting genome plasticity and allowing microbes to accumulate genetic diversity. MGEs serve as a vast communal gene pool and include DNA elements such as plasmids and bacteriophages (phages) among others. These mobile DNA elements represent a human health risk as they can introduce new traits, such as antibiotic resistance or virulence, to a bacterial strain. Sequencing libraries targeting environmental circular MGEs, referred to as metamobilomes, may broaden our current understanding of the mechanisms behind the mobility, prevalence and content of these elements. However, metamobilomics is affected by a severe bias towards small circular elements, introduced by multiple displacement amplification (MDA). MDA is typically used to overcome limiting DNA quantities after the removal of non-circular DNA during library preparations. By examining the relationship between sequencing coverage and the size of circular MGEs in paired metamobilome datasets with and without MDA, we show that larger circular elements are lost when using MDA. This study is the first to systematically demonstrate that MDA is detrimental to detecting larger-sized plasmids if small plasmids are present. It is also the first to show that MDA can be omitted when using enzyme-based DNA fragmentation and PCR in library preparation kits such as Nextera XT® from Illumina.

Rapid growth inhibitory activity of a YafQ-family endonuclease toxin of the Helicobacter pylori tfs4 integrative and conjugative element

Prokaryotic and archaeal chromosomes encode a diversity of toxin–antitoxin (TA) systems that contribute to a variety of stress-induced cellular processes in addition to stability and maintenance of mobile elements. Here, we find DinJ-YafQ family TA systems to be broadly distributed amongst diverse phyla, consistent with other ParE/RelE superfamily TAs, but more unusually occurring as a multiplicity of species-specific subtypes. In the gastric pathogen Helicobacter pylori we identify six distinct subtypes, of which three are predominantly associated with the mobilome, including the disease-associated integrative and conjugative element (ICE), tfs4. Whereas, the ICE-encoded proteins have characteristic features of DinJ-YafQ family Type II TA systems in general, the toxin component is distinguished by a broad metal-ion-dependent endonuclease activity with specificity for both RNA and DNA. We show that the remarkably rapid growth inhibitory activity of the ICE toxin is a correlate of a C-terminal lysine doublet which likely augments catalytic activity by increasing the positive electrostatic potential in the vicinity of the conserved active site. Our collective results reveal a structural feature of an ICE TA toxin that influences substrate catalysis and toxin function which may be relevant to specific TA-mediated responses in diverse genera of bacteria.

Evaluation of Putative Type II Toxin-Antitoxin Systems and Lon Protease Expression in Shigella flexneri Following Infection of Caco-2 Cells

: Shigella flexneri causes bacillary dysentery in developing countries. Due to recent reports regarding antimicrobial resistance in human S. flexneri, finding alternative therapeutics is of vital importance. Toxin-antitoxin (TA) systems have recently been introduced as antimicrobial targets owing to their involvement in bacterial survival in stress conditions and “persister” cell formation. In this study, the presence of four TA loci were studied in S. flexneri ATCC 12022. The presence of genes coding for the identified TA loci and Lon protease were confirmed by the PCR method using specific primers. Caco-2 cell lines were then infected with this standard strain, and 8 and 24 h post-infection, expression levels of genes coding for the studied TA loci, and Lon protease were evaluated using a real-time PCR method. Expression of mazF, GNAT (Gcn5-related N-acetyltransferase), yeeU, pfam13975, and Lon genes showed 5.4, 9.8, 2.3, 2.7, and 13.8-fold increase, respectively, 8 h after bacterial invasion of the Caco-2 cell line. In addition, the expression of the aforementioned genes showed 4.8, 10.8, 2.3, 3.7, and 16.8-fold increase after 24 h. The GNAT and lon genes showed significantly higher expression levels compared to the control (P value < 0.05). However, the increase in the expression level of yeeU was the same at 8 h and 24 h post-infection. In addition, mazF expression level showed a slight decrease at 24 h compared to 8h post-infection. Genes coding for GNAT and Lon protease showed a significantly higher expression after invading the Caco-2 cell line. Therefore, targeting GNAT or Lon protease can be taken into consideration for finding novel antimicrobial drug strategies. The exact functions and mechanisms of TA systems in S. flexneri isolates are suggested to be experimentally determined.

Complete genome sequence of Helicobacter pylori B128 7.13 and a single-step method for the generation of unmarked mutations

BACKGROUND

Helicobacter pylori represents an interesting model of bacterial pathogenesis given that most infections are asymptomatic, while a minority of infections cause severe gastric disease. H pylori strain B128 7.13 is used extensively to understand H pylori pathophysiology. Due to extensive restriction-modification systems, the fact that only some H pylori strains are naturally transformable, the inability of common plasmid and transposon vectors to replicate in this bacterium, as well as the limited number of antibiotic cassettes that are functional in H pylori, there are relatively few genetic tools for the mutagenesis of this bacterium.

MATERIALS AND METHODS

Here, we use PacBio and Illumina sequencing to reveal the complete genome sequence of H pylori B128 7.13. Furthermore, we describe a system to generate markerless and scarless mutations on the H pylori chromosome using the counter-selection marker, galactokinase from Escherichia coli.

RESULTS

We show that this mutagenesis strategy can be used to generate in-frame insertions, gene deletions, and multiple independent mutations in B128 7.13. Using the closed genome as a reference, we also report the absence of second site chromosomal mutations and/or rearrangements in our mutagenized strains. We compare the genome sequence of H pylori B128 7.13 with a closely related strain, H pylori B8, and reveal one notable region of difference, which is a 1430 bp insertion encoding a H pylori-specific DUF874 family protein of unknown function.

CONCLUSIONS

This article reports the closed genome of the important H pylori B128 7.13 strain and a mutagenesis method that can be adopted by researchers as an alternative strategy to generate isogenic mutants of H pylori in order to further our understanding of this bacterium.

Identification and first characterization of DinJ-YafQ toxin-antitoxin systems in Lactobacillus species of biotechnological interest

DinJ-YafQ is a type II TA system comprising the ribosome-dependent RNase YafQ toxin and the DinJ antitoxin protein. Although the module has been extensively characterized in Escherichia coli, little information is available for homologous systems in lactic acid bacteria. In this study, we employed bioinformatics tools to identify DinJ-YafQ systems in Lactobacillus casei, Lactobacillus paracasei and Lactobacillus rhamnosus species, commonly used in biotechnological processes. Among a total of nineteen systems found, two TA modules from Lactobacillus paracasei and two modules from Lactobacillus rhamnosus wild strains were isolated and their activity was verified by growth assays in Escherichia coli either in liquid and solid media. The RNase activity of the YafQ toxins was verified in vivo by probing mRNA dynamics and metabolism with single-cell Thioflavin T fluorescence. Our findings demonstrate that, albeit DinJ-YafQ TA systems are widely distributed in lactic acid bacteria, only few are fully functional, while others have lost toxicity even though they maintain high sequence identity with wild type YafQ and a likely functional antitoxin protein.

Modulation of Toxin-Antitoxin System Rnl AB Type II in Phage-Resistant Gammaproteobacteria Surviving Photodynamic Treatment

Type II toxin-antitoxin (TA) systems are the particular type of TA modules which take part in different kinds of cellular actions, such as biofilm formation, persistence, stress endurance, defense of the bacterial cell against multiple phage attacks, plasmid maintenance, and programmed cell death in favor of bacterial population. Although several bioinformatics and Pet lab studies have already been conducted to understand the functionality of already discovered TA systems, still, more work in this area is required. Rnl AB type II TA module, which is composed of RnlA toxin and RnlB antitoxin, is a newly discovered type II TA module which takes part in the defense mechanism against T4 bacteriophage attack in Escherichia coli K-12 strain MH1 that has not been widely studied in other bacteria. Because of the significant role of class Gammaproteobacteriacea in a diverse range of health problems, we chose here to focus on this class to survey the presence of the Rnl AB TA module. For better categorization and description of the distribution of this module in this class of bacteria, the corresponding phylogenetic trees are illustrated here. Neighbor-joining and the maximum parsimony methods were used in this study to take a look at the distribution of domains present in RnlA and RnlB proteins, among members of Gammaproteobacteria. Also, the possible roles of photodynamic therapy (PDT) in providing a substrate for better phage therapy are herein discussed.

Mechanisms for Differential Protein Production in Toxin-Antitoxin Systems

Toxin-antitoxin (TA) systems are key regulators of bacterial persistence, a multidrug-tolerant state found in bacterial species that is a major contributing factor to the growing human health crisis of antibiotic resistance. Type II TA systems consist of two proteins, a toxin and an antitoxin; the toxin is neutralized when they form a complex. The ratio of antitoxin to toxin is significantly greater than 1.0 in the susceptible population (non-persister state), but this ratio is expected to become smaller during persistence. Analysis of multiple datasets (RNA-seq, ribosome profiling) and results from translation initiation rate calculators reveal multiple mechanisms that ensure a high antitoxin-to-toxin ratio in the non-persister state. The regulation mechanisms include both translational and transcriptional regulation. We classified E. coli type II TA systems into four distinct classes based on the mechanism of differential protein production between toxin and antitoxin. We find that the most common regulation mechanism is translational regulation. This classification scheme further refines our understanding of one of the fundamental mechanisms underlying bacterial persistence, especially regarding maintenance of the antitoxin-to-toxin ratio.

The DinJ/RelE Toxin-Antitoxin System Suppresses Bacterial Proliferation and Virulence of Xylella fastidiosa in Grapevine

Xylella fastidiosa, the causal agent of Pierce's disease of grapes, is a slow-growing, xylem-limited, bacterial pathogen. Disease progression is characterized by systemic spread of the bacterium through xylem vessel networks, causing leaf-scorching symptoms, senescence, and vine decline. It appears to be advantageous to this pathogen to avoid excessive blockage of xylem vessels, because living bacterial cells are generally found in plant tissue with low bacterial cell density and minimal scorching symptoms. The DinJ/RelE toxin-antitoxin system is characterized here for a role in controlling bacterial proliferation and population size during plant colonization. The DinJ/RelE locus is transcribed from two separate promoters, allowing for coexpression of antitoxin DinJ with endoribonuclease toxin RelE, in addition to independent expression of RelE. The ratio of antitoxin/toxin expressed is dependent on bacterial growth conditions, with lower amounts of antitoxin present under conditions designed to mimic grapevine xylem sap. A knockout mutant of DinJ/RelE exhibits a hypervirulent phenotype, with higher bacterial populations and increased symptom development and plant decline. It is likely that DinJ/RelE acts to prevent excessive population growth, contributing to the ability of the pathogen to spread systemically without completely blocking the xylem vessels and increasing probability of acquisition by the insect vector.

Toxin YafQ Reduces Escherichia coli Growth at Low Temperatures

Toxin/antitoxin (TA) systems reduce metabolism under stress; for example, toxin YafQ of the YafQ/DinJ Escherichia coli TA system reduces growth by cleaving transcripts with in-frame 5'-AAA-G/A-3' sites, and antitoxin DinJ is a global regulator that represses its locus as well as controls levels of the stationary sigma factor RpoS. Here we investigated the influence on cell growth at various temperatures and found that deletion of the antitoxin gene, dinJ, resulted in both reduced metabolism and slower growth at 18°C but not at 37°C. The reduction in growth could be complemented by producing DinJ from a plasmid. Using a transposon screen to reverse the effect of the absence of DinJ, two mutations were found that inactivated the toxin YafQ; hence, the toxin caused the slower growth only at low temperatures rather than DinJ acting as a global regulator. Corroborating this result, a clean deletion of yafQ in the ΔdinJ ΔKmR strain restored both metabolism and growth at 18°C. In addition, production of YafQ was more toxic at 18°C compared to 37°C. Furthermore, by overproducing all the E. coli proteins, the global transcription repressor Mlc was found that counteracts YafQ toxicity only at 18°C. Therefore, YafQ is more effective at reducing metabolism at low temperatures, and Mlc is its putative target.

Toxin-Antitoxin Modules Are Pliable Switches Activated by Multiple Protease Pathways

Toxin-antitoxin (TA) modules are bacterial regulatory switches that facilitate conflicting outcomes for cells by promoting a pro-survival phenotypic adaptation and/or by directly mediating cell death, all through the toxin activity upon degradation of antitoxin. Intensive study has revealed specific details of TA module functions, but significant gaps remain about the molecular details of activation via antitoxin degradation used by different bacteria and in different environments. This review summarizes the current state of knowledge about the interaction of antitoxins with cellular proteases Lon and ClpP to mediate TA module activation. An understanding of these processes can answer long-standing questions regarding stochastic versus specific activation of TA modules and provide insight into the potential for manipulation of TA modules to alter bacterial growth.

Keeping the Wolves at Bay: Antitoxins of Prokaryotic Type II Toxin-Antitoxin Systems

In their initial stages of discovery, prokaryotic toxin-antitoxin (TA) systems were confined to bacterial plasmids where they function to mediate the maintenance and stability of usually low- to medium-copy number plasmids through the post-segregational killing of any plasmid-free daughter cells that developed. Their eventual discovery as nearly ubiquitous and repetitive elements in bacterial chromosomes led to a wealth of knowledge and scientific debate as to their diversity and functionality in the prokaryotic lifestyle. Currently categorized into six different types designated types I–VI, type II TA systems are the best characterized. These generally comprised of two genes encoding a proteic toxin and its corresponding proteic antitoxin, respectively. Under normal growth conditions, the stable toxin is prevented from exerting its lethal effect through tight binding with the less stable antitoxin partner, forming a non-lethal TA protein complex. Besides binding with its cognate toxin, the antitoxin also plays a role in regulating the expression of the type II TA operon by binding to the operator site, thereby repressing transcription from the TA promoter. In most cases, full repression is observed in the presence of the TA complex as binding of the toxin enhances the DNA binding capability of the antitoxin. TA systems have been implicated in a gamut of prokaryotic cellular functions such as being mediators of programmed cell death as well as persistence or dormancy, biofilm formation, as defensive weapons against bacteriophage infections and as virulence factors in pathogenic bacteria. It is thus apparent that these antitoxins, as DNA-binding proteins, play an essential role in modulating the prokaryotic lifestyle whilst at the same time preventing the lethal action of the toxins under normal growth conditions, i.e., keeping the proverbial wolves at bay. In this review, we will cover the diversity and characteristics of various type II TA antitoxins. We shall also look into some interesting deviations from the canonical type II TA systems such as tripartite TA systems where the regulatory role is played by a third party protein and not the antitoxin, and a unique TA system encoding a single protein with both toxin as well as antitoxin domains.

Molecular basis of ribosome recognition and mRNA hydrolysis by the E. coli YafQ toxin

Bacterial type II toxin-antitoxin modules are protein–protein complexes whose functions are finely tuned by rapidly changing environmental conditions. E. coli toxin YafQ is suppressed under steady state growth conditions by virtue of its interaction with its cognate antitoxin, DinJ. During stress, DinJ is proteolytically degraded and free YafQ halts translation by degrading ribosome-bound mRNA to slow growth until the stress has passed. Although structures of the ribosome with toxins RelE and YoeB have been solved, it is unclear what residues among ribosome-dependent toxins are essential for mediating both recognition of the ribosome and the mRNA substrate given their low sequence identities. Here we show that YafQ coordinates binding to the 70S ribosome via three surface-exposed patches of basic residues that we propose directly interact with 16S rRNA. We demonstrate that YafQ residues H50, H63, D67 and H87 participate in acid-base catalysis during mRNA hydrolysis and further show that H50 and H63 functionally complement as general bases to initiate the phosphodiester cleavage reaction. Moreover YafQ residue F91 likely plays an important role in mRNA positioning. In summary, our findings demonstrate the plasticity of ribosome-dependent toxin active site residues and further our understanding of which toxin residues are important for function.

Activation of Toxin-Antitoxin System Toxins Suppresses Lethality Caused by the Loss of σE in Escherichia coli

σ(E), an alternative σ factor that governs a major signaling pathway in envelope stress responses in Gram-negative bacteria, is essential for growth of Escherichia coli not only under stressful conditions, such as elevated temperature, but also under normal laboratory conditions. A mutational inactivation of the hicB gene has been reported to suppress the lethality caused by the loss of σ(E). hicB encodes the antitoxin of the HicA-HicB toxin-antitoxin (TA) system; overexpression of the HicA toxin, which exhibits mRNA interferase activity, causes cleavage of mRNAs and an arrest of cell growth, while simultaneous expression of HicB neutralizes the toxic effects of overproduced HicA. To date, however, how the loss of HicB rescues the cell lethality in the absence of σ(E) and, more specifically, whether HicA is involved in this process remain unknown. Here we showed that simultaneous disruption of hicA abolished suppression of the σ(E) essentiality in the absence of hicB, while ectopic expression of wild-type HicA, but not that of its mutant forms without mRNA interferase activity, restored the suppression. Furthermore, HicA and two other mRNA interferase toxins, HigB and YafQ, suppressed the σ(E) essentiality even in the presence of chromosomally encoded cognate antitoxins when these toxins were overexpressed individually. Interestingly, when the growth media were supplemented with low levels of antibiotics that are known to activate toxins, E. coli cells with no suppressor mutations grew independently of σ(E). Taken together, our results indicate that the activation of TA system toxins can suppress the σ(E) essentiality and affect the extracytoplasmic stress responses.

IMPORTANCE

σ(E) is an alternative σ factor involved in extracytoplasmic stress responses. Unlike other alternative σ factors, σ(E) is indispensable for the survival of E. coli even under unstressed conditions, although the exact reason for its essentiality remains unknown. Toxin-antitoxin (TA) systems are widely distributed in prokaryotes and are composed of two adjacent genes, encoding a toxin that exerts harmful effects on the toxin-producing bacterium itself and an antitoxin that neutralizes the cognate toxin. Curiously, it is known that inactivation of an antitoxin rescues the σ(E) essentiality, suggesting a connection between TA systems and σ(E) function. We demonstrate here that toxin activation is necessary for this rescue and suggest the possible involvement of TA systems in extracytoplasmic stress responses.

The bacterial translation stress response

Throughout their life, bacteria need to sense and respond to environmental stress. Thus, such stress responses can require dramatic cellular reprogramming, both at the transcriptional as well as the translational level. This review focuses on the protein factors that interact with the bacterial translational apparatus to respond to and cope with different types of environmental stress. For example, the stringent factor RelA interacts with the ribosome to generate ppGpp under nutrient deprivation, whereas a variety of factors have been identified that bind to the ribosome under unfavorable growth conditions to shut-down (RelE, pY, RMF, HPF and EttA) or re-program (MazF, EF4 and BipA) translation. Additional factors have been identified that rescue ribosomes stalled due to stress-induced mRNA truncation (tmRNA, ArfA, ArfB), translation of unfavorable protein sequences (EF-P), heat shock-induced subunit dissociation (Hsp15), or antibiotic inhibition (TetM, FusB). Understanding the mechanism of how the bacterial cell responds to stress will not only provide fundamental insight into translation regulation, but will also be an important step to identifying new targets for the development of novel antimicrobial agents.

Mechanisms of toxin inhibition and transcriptional repression by Escherichia coli DinJ-YafQ

Bacteria encounter environmental stresses that regulate a gene expression program required for adaptation and survival. Here, we report the 1.8-Å crystal structure of the Escherichia coli toxin-antitoxin complex YafQ-(DinJ)2-YafQ, a key component of the stress response. The antitoxin DinJ dimer adopts a ribbon-helix-helix motif required for transcriptional autorepression, and toxin YafQ contains a microbial RNase fold whose proposed active site is concealed by DinJ binding. Contrary to previous reports, our studies indicate that equivalent levels of transcriptional repression occur by direct interaction of either YafQ-(DinJ)2-YafQ or a DinJ dimer at a single inverted repeat of its recognition sequence that overlaps with the -10 promoter region. Surprisingly, multiple YafQ-(DinJ)2-YafQ complexes binding to the operator region do not appear to amplify the extent of repression. Our results suggest an alternative model for transcriptional autorepression that may be novel to DinJ-YafQ.

Structural and functional characterization of Escherichia coli toxin-antitoxin complex DinJ-YafQ

Toxin YafQ functions as a ribonuclease in the dinJ-yafQ toxin-antitoxin system of Escherichia coli. Antitoxin DinJ neutralizes YafQ-mediated toxicity by forming a stable protein complex. Here, crystal structures of the (DinJ)2-(YafQ)2 complex and the isolated YafQ toxin have been determined. The structure of the heterotetrameric complex (DinJ)2-(YafQ)2 revealed that the N-terminal region of DinJ folds into a ribbon-helix-helix motif and dimerizes for DNA recognition, and the C-terminal portion of each DinJ exclusively wraps around a YafQ molecule. Upon incorporation into the heterotetrameric complex, a conformational change of YafQ in close proximity to the catalytic site of the typical microbial ribonuclease fold was observed and validated. Mutagenesis experiments revealed that a DinJ mutant restored YafQ RNase activity in a tetramer complex in vitro but not in vivo. An electrophoretic mobility shift assay showed that one of the palindromic sequences present in the upstream intergenic region of DinJ served as a binding sequences for both the DinJ-YafQ complex and the antitoxin DinJ alone. Based on structure-guided and site-directed mutagenesis of DinJ-YafQ, we showed that two pairs of amino acids in DinJ were important for DNA binding; the R8A and K16A substitutions and the S31A and R35A substitutions in DinJ abolished the DNA binding ability of the DinJ-YafQ complex.

Crystal structure of toxin HP0892 from Helicobacter pylori with two Zn(II) at 1.8 Å resolution

Antibiotic resistance and microorganism virulence have been consistently exhibited by bacteria and archaea, which survive in conditions of environmental stress through toxin-antitoxin (TA) systems. The HP0892-HP0893 TA system is one of the two known TA systems belonging to Helicobacter pylori. The antitoxin, HP0893, binds and inhibits the HP0892 toxin and regulates the transcription of the TA operon. Here, we present the crystal structure of the zinc-bound HP0892 toxin at 1.8 Å resolution. Reorientation of residues at the mRNase active site was shown. The involved residues, namely E58A, H86A, and H58A/ H60A, were mutated and the binding affinity was monitored by ITC studies. Through the structural difference between the apo and the metal-bound state, and using a homology modeling tool, the involvement of the metal ion in mRNase active site could be identified. The most catalytically important residue, His86, reorients itself to exhibit RNase activity. His47, Glu58, and His60 are involved in metal binding where Glu58 acts as a general base and His47 and His60 may also act as a general acid in enzymatic activity. Glu58 and Asp64 are involved in substrate binding and specific sequence recognition. Arg83 is involved in phosphate binding and stabilization of the transition state, and Phe90 is involved in base packing and substrate orientation.

Post-transcriptional regulation of gene expression in bacterial pathogens by toxin-antitoxin systems

Toxin-antitoxin (TA) systems are small genetic elements ubiquitous in prokaryotic genomes that encode toxic proteins targeting various vital cellular functions. Typically, toxin activity is controlled by adjacently encoded protein or RNA antitoxins and unleashed as a consequence of genetic fluctuations or stressful conditions. Whereas some TA systems interfere with replication or cell wall synthesis, most of them influence transcriptional and post-transcriptional gene regulation. Antitoxin proteins often act as DNA binding transcriptional regulators and many TA toxins exhibit endoribonuclease activity to selectively degrade different RNA species and thus alter gene expression patterns. Some TA RNases cleave tRNA, tmRNAs or rRNAs, whereas most commonly mRNAs either in association with the ribosome or as free transcripts, are targeted. Examples are provided on how TA toxins differentially shape gene expression in bacterial pathogens by creating specialized ribosomes or by altering the transcriptome and how this may be tied in the control of pathogenicity factors.

Toxin-antitoxin systems as multilevel interaction systems

Toxin-antitoxin (TA) systems are small genetic modules usually composed of a toxin and an antitoxin counteracting the activity of the toxic protein. These systems are widely spread in bacterial and archaeal genomes. TA systems have been assigned many functions, ranging from persistence to DNA stabilization or protection against mobile genetic elements. They are classified in five types, depending on the nature and mode of action of the antitoxin. In type I and III, antitoxins are RNAs that either inhibit the synthesis of the toxin or sequester it. In type II, IV and V, antitoxins are proteins that either sequester, counterbalance toxin activity or inhibit toxin synthesis. In addition to these interactions between the antitoxin and toxin components (RNA-RNA, protein-protein, RNA-protein), TA systems interact with a variety of cellular factors, e.g., toxins target essential cellular components, antitoxins are degraded by RNAses or ATP-dependent proteases. Hence, TA systems have the capacity to interact with each other at different levels. In this review, we will discuss the different interactions in which TA systems are involved and their implications in TA system functions and evolution.

Light activation of Staphylococcus aureus toxin YoeBSa1 reveals guanosine-specific endoribonuclease activity

The Staphylococcus aureus chromosome harbors two homologues of the YefM-YoeB toxin-antitoxin (TA) system. The toxins YoeBSa1 and YoeBSa2 possess ribosome-dependent ribonuclease (RNase) activity in Escherichia coli. This activity is similar to that of the E. coli toxin YoeBEc, an enzyme that, in addition to ribosome-dependent RNase activity, possesses ribosome-independent RNase activity in vitro. To investigate whether YoeBSa1 is also a ribosome-independent RNase, we expressed YoeBSa1 using a novel strategy and characterized its in vitro RNase activity, sequence specificity, and kinetics. Y88 of YoeBSa1 was critical for in vitro activity and cell culture toxicity. This residue was mutated to o-nitrobenzyl tyrosine (ONBY) via unnatural amino acid mutagenesis. YoeBSa1-Y88ONBY could be expressed in the absence of the antitoxin YefMSa1 in E. coli. Photocaged YoeBSa1-Y88ONBY displayed UV light-dependent RNase activity toward free mRNA in vitro. The in vitro ribosome-independent RNase activity of YoeBSa1-Y88ONBY, YoeBSa1-Y88F, and YoeBSa1-Y88TAG was significantly reduced or abolished. In contrast to YoeBEc, which cleaves RNA at both adenosine and guanosine with a preference for adenosine, YoeBSa1 cleaved mRNA specifically at guanosine. Using this information, a fluorometric assay was developed and used to determine the kinetic parameters for ribosome-independent RNA cleavage by YoeBSa1.

2'-Phosphate cyclase activity of RtcA: a potential rationale for the operon organization of RtcA with an RNA repair ligase RtcB in Escherichia coli and other bacterial taxa

RNA terminal phosphate cyclase catalyzes the ATP-dependent conversion of a 3'-phosphate RNA end to a 2',3'-cyclic phosphate via covalent enzyme-(histidinyl-Nε)-AMP and RNA(3')pp(5')A intermediates. Here, we report that Escherichia coli RtcA (and its human homolog Rtc1) are capable of cyclizing a 2'-phosphate RNA end in high yield. The rate of 2'-phosphate cyclization by RtcA is five orders of magnitude slower than 3'-phosphate cyclization, notwithstanding that RtcA binds with similar affinity to RNA3'p and RNA2'p substrates. These findings expand the functional repertoire of RNA cyclase and suggest that phosphate geometry during adenylate transfer to RNA is a major factor in the kinetics of cyclization. RtcA is coregulated in an operon with an RNA ligase, RtcB, that splices RNA 5'-OH ends to either 3'-phosphate or 2',3'-cyclic phosphate ends. Our results suggest that RtcA might serve an end healing function in an RNA repair pathway, by converting RNA 2'-phosphates, which cannot be spliced by RtcB, to 2',3'-cyclic phosphates that can be sealed. The rtcBA operon is controlled by the σ(54) coactivator RtcR encoded by an adjacent gene. This operon arrangement is conserved in diverse bacterial taxa, many of which have also incorporated the RNA-binding protein Ro (which is implicated in RNA quality control under stress conditions) as a coregulated component of the operon.

Crystal structure of apo and copper bound HP0894 toxin from Helicobacter pylori 26695 and insight into mRNase activity

The toxin-antitoxin (TA) systems widely spread among bacteria and archaea are important for antibiotic resistance and microorganism virulence. The bacterial kingdom uses TA systems to adjust the global level of gene expression and translation through RNA degradation. In Helicobacter pylori, only two TA systems are known thus far. Our previous studies showed that HP0894-HP0895 acts as a TA system and that HP0894 exhibits intrinsic RNase activity. However, the precise molecular basis for interaction with substrate or antitoxin and the mechanism of mRNA cleavage remain unclear. Therefore, in an attempt to shed some light on the mechanism behind the TA system of HP0894-HP0895, here we present the crystal structures of apo- and metal-bound H. pylori 0894 at 1.28Å and 1.89Å, respectively. Through the combined approach of structural analysis and structural homology search, the amino acids involved in mRNase active site were monitored and the reorientations of different residues were discussed in detail. In the mRNase active site of HP0894 toxin, His84 acts as a catalytic residue and reorients itself to exhibit this type of activity, acting as a general acid in an acid-base catalysis reaction, while His47 and His60 stabilize the transition state. Lys52, Glu58, Asp64 and Arg80 have phosphate binding and specific sequence recognition. Glu58 also acts as a general base, and substrate reorientation is caused by Phe88. Based on experimental findings, a model for antitoxin binding could be suggested.

Identification and characterization of type II toxin-antitoxin systems in the opportunistic pathogen Acinetobacter baumannii

Acinetobacter baumannii is an opportunistic pathogen that causes nosocomial infections. Due to the ability to persist in the clinical environment and rapidly acquire antibiotic resistance, multidrug-resistant A. baumannii clones have spread in medical units in many countries in the last decade. The molecular basis of the emergence and spread of the successful multidrug-resistant A. baumannii clones is not understood. Bacterial toxin-antitoxin (TA) systems are abundant genetic loci harbored in low-copy-number plasmids and chromosomes and have been proposed to fulfill numerous functions, from plasmid stabilization to regulation of growth and death under stress conditions. In this study, we have performed a thorough bioinformatic search for type II TA systems in genomes of A. baumannii strains and estimated at least 15 possible TA gene pairs, 5 of which have been shown to be functional TA systems. Three of them were orthologs of bacterial and archaeal RelB/RelE, HicA/HicB, and HigB/HigA systems, and others were the unique SplT/SplA and CheT/CheA TA modules. The toxins of all five TA systems, when expressed in Escherichia coli, inhibited translation, causing RNA degradation. The HigB/HigA and SplT/SplA TA pairs of plasmid origin were highly prevalent in clinical multidrug-resistant A. baumannii isolates from Lithuanian hospitals belonging to the international clonal lineages known as European clone I (ECI) and ECII.

Relaxed cleavage specificity within the RelE toxin family

Bacterial type II toxin-antitoxin systems are widespread in bacteria. Among them, the RelE toxin family is one of the most abundant. The RelE(K-12) toxin of Escherichia coli K-12 represents the paradigm for this family and has been extensively studied, both in vivo and in vitro. RelE(K-12) is an endoribonuclease that cleaves mRNAs that are translated by the ribosome machinery as these transcripts enter the A site. Earlier in vivo reports showed that RelE(K-12) cleaves preferentially in the 5'-end coding region of the transcripts in a codon-independent manner. To investigate whether the molecular activity as well as the cleavage pattern are conserved within the members of this toxin family, RelE-like sequences were selected in Proteobacteria, Cyanobacteria, Actinobacteria, and Spirochaetes and tested in E. coli. Our results show that these RelE-like sequences are part of toxin-antitoxin gene pairs, and that they inhibit translation in E. coli by cleaving transcripts that are being translated. Primer extension analyses show that these toxins exhibit specific cleavage patterns in vivo, both in terms of frequency and location of cleavage sites. We did not observe codon-dependent cleavage but rather a trend to cleave upstream purines and between the second and third positions of codons, except for the actinobacterial toxin. Our results suggest that RelE-like toxins have evolved to rapidly and efficiently shut down translation in a large spectrum of bacterial species, which correlates with the observation that toxin-antitoxin systems are spreading by horizontal gene transfer.

Complete sequence of broad-host-range plasmid pNOR-2000 harbouring the metallo-β-lactamase gene blaVIM-2 from Pseudomonas aeruginosa

OBJECTIVES

Metallo-β-lactamases (MBLs) are increasingly reported not only in Enterobacteriaceae but also in Pseudomonas spp. These enzymes hydrolyse all β-lactams, including carbapenems, and are not inhibited by β-lactamase inhibitors. The aim of this study was to fully characterize a plasmid bearing the blaVIM-2 MBL gene identified in a Pseudomonas aeruginosa isolate.

METHODS

This plasmid was fully sequenced by high-density pyrosequencing and annotated using the GenDB version 2.0 annotation tool. The evaluation of the broad-host-range replication of the pNOR-2000 replication initiation gene was assessed using electro-transformation and conjugation assays and the distribution of this replicase gene was evaluated using an international collection of VIM-producing Pseudomonas spp.

RESULTS

Analysis of the 21 880 bp sequence of pNOR-2000 revealed a truncated and non-functional transfer operon, in addition to novel genes encoding a serine protease and toxin/antitoxin addiction systems. This broad-host-range plasmid shares high gene synteny with part of the mobile genomic island pKLC102 identified in P. aeruginosa strain C.

CONCLUSIONS

We report here the complete nucleotide sequence of plasmid pNOR-2000 from a P. aeruginosa clinical isolate harbouring the integron-located MBL gene blaVIM-2.

Identification of chromosomal HP0892-HP0893 toxin-antitoxin proteins in Helicobacter pylori and structural elucidation of their protein-protein interaction

Bacterial chromosomal toxin-antitoxin (TA) systems have been proposed not only to play an important role in the stress response, but also to be associated with antibiotic resistance. Here, we identified the chromosomal HP0892-HP0893 TA proteins in the gastric pathogen, Helicobacter pylori, and structurally characterized their protein-protein interaction. Previously, HP0892 protein was suggested to be a putative TA toxin based on its structural similarity to other RelE family TA toxins. In this study, we demonstrated that HP0892 binds to HP0893 strongly with a stoichiometry of 1:1, and HP0892-HP0893 interaction occurs mainly between the N-terminal secondary structure elements of HP0892 and the C-terminal region of HP0893. HP0892 cleaved mRNA in vitro, preferentially at the 5' end of A or G, and the RNase activity of HP0892 was inhibited by HP0893. In addition, heterologous expression of HP0892 in Escherichia coli cells led to cell growth arrest, and the cell toxicity of HP0892 was neutralized by co-expression with HP0893. From these results and a structural comparison with other TA toxins, it is concluded that HP0892 is a toxin with intrinsic RNase activity and HP0893 is an antitoxin against HP0892 from a TA system of H. pylori. It has been known that hp0893 gene and another TA antitoxin gene, hp0895, of H. pylori, are both genomic open reading frames that correspond to genes that are potentially expressed in response to interactions with the human gastric mucosa. Therefore, it is highly probable that TA systems of H. pylori are involved in virulence of H. pylori.