Model for RNA Binding and the Catalytic Site of the RNase Kid of the Bacterial parD Toxin–Antitoxin System

PemK's Arg24 is a crucial residue for PemIK toxin-antitoxin system to induce the persistence of Weissella cibaria against ciprofloxacin stress

The toxin-antitoxin (TA) system plays a key role in bacteria escaping antibiotic stress with persistence, however, the mechanisms by which persistence is controlled remain poorly understood. Weissella cibaria, a novel probiotic, can enters a persistent state upon encountering ciprofloxacin stress. Conversely, it resumes from the persistence when ciprofloxacin stress is relieved or removed. Here, it was found that PemIK TA system played a role in transitioning between these two states. And the PemIK was consisted of PemK, an endonuclease toxic to mRNA, and antitoxin PemI which neutralized its toxicity. The PemK specifically cleaved the U↓AUU in mRNA encoding enzymes involved in glycolysis, TCA cycle and respiratory chain pathways. This cleavage event subsequently disrupted the crucial cellular processes such as hydrogen transfer, electron transfer, NADH and FADH2 synthesis, ultimately leading to a decrease in ATP levels and an increase in membrane depolarization and persister frequency. Notably, Arg24 was a critical active residue for PemK, its mutation significantly reduced the mRNA cleavage activity and the adverse effects on metabolism. These insights provided a clue to comprehensively understand the mechanism by which PemIK induced the persistence of W. cibaria to escape ciprofloxacin stress, thereby highlighting another novel aspect PemIK respond for antibiotic stress.

OUP accepted manuscript

Staphylococcus aureus is a notorious and globally distributed pathogenic bacterium. New strategies to develop novel antibiotics based on intrinsic bacterial toxin–antitoxin (TA) systems have been recently reported. Because TA systems are present only in bacteria and not in humans, these distinctive systems are attractive targets for developing antibiotics with new modes of action. S. aureus PemIK is a type II TA system, comprising the toxin protein PemK and the labile antitoxin protein PemI. Here, we determined the crystal structures of both PemK and the PemIK complex, in which PemK is neutralized by PemI. Our biochemical approaches, including fluorescence quenching and polarization assays, identified Glu20, Arg25, Thr48, Thr49, and Arg84 of PemK as being important for RNase function. Our study indicates that the active site and RNA-binding residues of PemK are covered by PemI, leading to unique conformational changes in PemK accompanied by repositioning of the loop between β1 and β2. These changes can interfere with RNA binding by PemK. Overall, PemK adopts particular open and closed forms for precise neutralization by PemI. This structural and functional information on PemIK will contribute to the discovery and development of novel antibiotics in the form of peptides or small molecules inhibiting direct binding between PemI and PemK.

Applications of phage-derived RNA-based technologies in synthetic biology

As the most abundant biological entities with incredible diversity, bacteriophages (also known as phages) have been recognized as an important source of molecular machines for the development of genetic-engineering tools. At the same time, phages are crucial for establishing and improving basic theories of molecular biology. Studies on phages provide rich sources of essential elements for synthetic circuit design as well as powerful support for the improvement of directed evolution platforms. Therefore, phages play a vital role in the development of new technologies and central scientific concepts. After the RNA world hypothesis was proposed and developed, novel biological functions of RNA continue to be discovered. RNA and its related elements are widely used in many fields such as metabolic engineering and medical diagnosis, and their versatility led to a major role of RNA in synthetic biology. Further development of RNA-based technologies will advance synthetic biological tools as well as provide verification of the RNA world hypothesis. Most synthetic biology efforts are based on reconstructing existing biological systems, understanding fundamental biological processes, and developing new technologies. RNA-based technologies derived from phages will offer abundant sources for synthetic biological components. Moreover, phages as well as RNA have high impact on biological evolution, which is pivotal for understanding the origin of life, building artificial life-forms, and precisely reprogramming biological systems. This review discusses phage-derived RNA-based technologies terms of phage components, the phage lifecycle, and interactions between phages and bacteria. The significance of RNA-based technology derived from phages for synthetic biology and for understanding the earliest stages of biological evolution will be highlighted.

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.

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.

tRNA is a new target for cleavage by a MazF toxin

Toxin-antitoxin (TA) systems play key roles in bacterial persistence, biofilm formation and stress responses. The MazF toxin from the Escherichia coli mazEF TA system is a sequence- and single-strand-specific endoribonuclease, and many studies have led to the proposal that MazF family members exclusively target mRNA. However, recent data indicate some MazF toxins can cleave specific sites within rRNA in concert with mRNA. In this report, we identified the repertoire of RNAs cleaved by Mycobacterium tuberculosis toxin MazF-mt9 using an RNA-seq-based approach. This analysis revealed that two tRNAs were the principal targets of MazF-mt9, and each was cleaved at a single site in either the tRNA(Pro14) D-loop or within the tRNA(Lys43) anticodon. This highly selective target discrimination occurs through recognition of not only sequence but also structural determinants. Thus, MazF-mt9 represents the only MazF family member known to target tRNA and to require RNA structure for recognition and cleavage. Interestingly, the tRNase activity of MazF-mt9 mirrors basic features of eukaryotic tRNases that also generate stable tRNA-derived fragments that can inhibit translation in response to stress. Our data also suggest a role for tRNA distinct from its canonical adapter function in translation, as cleavage of tRNAs by MazF-mt9 downregulates bacterial growth.

Coupling between the basic replicon and the Kis-Kid maintenance system of plasmid R1: modulation by Kis antitoxin levels and involvement in control of plasmid replication

kis-kid, the auxiliary maintenance system of plasmid R1 and copB, the auxiliary copy number control gene of this plasmid, contribute to increase plasmid replication efficiency in cells with lower than average copy number. It is thought that Kis antitoxin levels decrease in these cells and that this acts as the switch that activates the Kid toxin; activated Kid toxin reduces copB-mRNA levels and this increases RepA levels that increases plasmid copy number. In support of this model we now report that: (i) the Kis antitoxin levels do decrease in cells containing a mini-R1 plasmid carrying a repA mutation that reduces plasmid copy number; (ii) kid-dependent replication rescue is abolished in cells in which the Kis antitoxin levels or the CopB levels are increased. Unexpectedly we found that this coordination significantly increases both the copy number of the repA mutant and of the wt mini-R1 plasmid. This indicates that the coordination between plasmid replication functions and kis-kid system contributes significantly to control plasmid R1 replication.

An RNA-seq method for defining endoribonuclease cleavage specificity identifies dual rRNA substrates for toxin MazF-mt3

Toxin-antitoxin (TA) systems are widespread in prokaryotes. Among these, the mazEF TA system encodes an endoribonucleolytic toxin, MazF, that inhibits growth by sequence-specific cleavage of single-stranded RNA. Defining the physiological targets of a MazF toxin first requires the identification of its cleavage specificity, yet the current toolkit is antiquated and limited. We describe a rapid genome-scale approach, MORE (Mapping by Overexpression of an RNase in Escherichia coli) RNA-seq, for defining the cleavage specificity of endoribonucleolytic toxins. Application of MORE RNA-seq to MazF-mt3 from Mycobacterium tuberculosis reveals two critical ribosomal targets — the essential, evolutionarily conserved helix/loop 70 of 23S rRNA and the anti-Shine-Dalgarno (aSD) sequence of 16S rRNA. Our findings support an emerging model where both rRNA and mRNA are principal targets of MazF toxins and suggest that, as in E. coli, removal of the aSD sequence by a MazF toxin modifies ribosomes to selectively translate leaderless mRNAs in M. tuberculosis.

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.

Structural basis of mRNA recognition and cleavage by toxin MazF and its regulation by antitoxin MazE in Bacillus subtilis

MazF is an mRNA interferase, which, upon activation during stress conditions, cleaves mRNAs in a sequence-specific manner, resulting in cellular growth arrest. During normal growth conditions, the MazF toxin is inactivated through binding to its cognate antitoxin, MazE. How MazF specifically recognizes its mRNA target and carries out cleavage and how the formation of the MazE-MazF complex inactivates MazF remain unclear. We present crystal structures of MazF in complex with mRNA substrate and antitoxin MazE in Bacillus subtilis. The structure of MazF in complex with uncleavable UUdUACAUAA RNA substrate defines the molecular basis underlying the sequence-specific recognition of UACAU and the role of residues involved in the cleavage through site-specific mutational studies. The structure of the heterohexameric (MazF)2-(MazE)2-(MazF)2 complex in Bacillus subtilis, supplemented by mutational data, demonstrates that the positioning of the C-terminal helical segment of MazE within the RNA-binding channel of the MazF dimer prevents mRNA binding and cleavage by MazF.

ACA-specific RNA sequence recognition is acquired via the loop 2 region of MazF mRNA interferase

MazF is an mRNA interferase that cleaves mRNAs at a specific RNA sequence. MazF from E. coli (MazF-ec) cleaves RNA at A^CA. To date, a large number of MazF homologs that cleave RNA at specific three- to seven-base sequences have been identified from bacteria to archaea. MazF-ec forms a dimer, in which the interface between the two subunits is known to be the RNA substrate-binding site. Here, we investigated the role of the two loops in MazF-ec, which are closely associated with the interface of the MazF-ec dimer. We examined whether exchanging the loop regions of MazF-ec with those from other MazF homologs, such as MazF from Myxococcus xanthus (MazF-mx) and MazF from Mycobacterium tuberculosis (MazF-mt3), affects RNA cleavage specificity. We found that exchanging loop 2 of MazF-ec with loop 2 regions from either MazF-mx or MazF-mt3 created a new cleavage sequence at (A/U)(A/U)AA^C in addition to the original cleavage site, A^CA, whereas exchanging loop 1 did not alter cleavage specificity. Intriguingly, exchange of loop 2 with 8 or 12 consecutive Gly residues also resulted in a new RNA cleavage site at (A/U)(A/U)AA^C. The present study suggests a method for expanding the RNA cleavage repertoire of mRNA interferases, which is crucial for potential use in the regulation of specific gene expression and for biotechnological applications.

Multiple fuzzy interactions in the moonlighting function of thymosin-β4

Thymosine β4 (Tß4) is a 43 amino acid long intrinsically disordered protein (IDP), which was initially identified as an actin-binding and sequestering molecule. Later it was described to have multiple other functions, such as regulation of endothelial cell differentiation, blood vessel formation, wound repair, cardiac cell migration, and survival.1 The various functions of Tβ4 are mediated by interactions with distinct and structurally unrelated partners, such as PINCH, ILK, and stabilin-2, besides the originally identified G-actin. Although the cellular readout of these interactions and the formation of these complexes have been thoroughly described, no attempt was made to study these interactions in detail, and to elucidate the thermodynamic, kinetic, and structural underpinning of this range of moonlighting functions. Because Tβ4 is mostly disordered, and its 4 described partners are structurally unrelated (the CTD of stabilin-2 is actually fully disordered), it occurred to us that this system might be ideal to characterize the structural adaptability and ensuing moonlighting functions of IDPs. Unexpectedly, we found that Tβ4 engages in multiple weak, transient, and fuzzy interactions, i.e., it is capable of mediating distinct yet specific interactions without adapting stable folded structures.

Selectivity and self-assembly in the control of a bacterial toxin by an antitoxic noncoding RNA pseudoknot

Bacterial small RNAs perform numerous regulatory roles, including acting as antitoxic components in toxin-antitoxin systems. In type III toxin-antitoxin systems, small processed RNAs directly antagonize their toxin protein partners, and in the systems characterized the toxin and antitoxin components together form a trimeric assembly. In the present study, we sought to define how the RNA antitoxin, ToxI, inhibits its potentially lethal protein partner, ToxN. We show through cross-inhibition experiments with the ToxIN systems from Pectobacterium atrosepticum (ToxIN(Pa)) and Bacillus thuringiensis (ToxIN(Bt)) that ToxI RNAs are highly selective enzyme inhibitors. Both systems have an "addictive" plasmid maintenance phenotype. We demonstrate that ToxI(Pa) can inhibit ToxN(Pa) in vitro both in its processed form and as a repetitive precursor RNA, and this inhibition is linked to the self-assembly of the trimeric complex. Inhibition and self-assembly are both mediated entirely by the ToxI(Pa) RNA, with no requirement for cellular factors or exogenous energy. Finally, we explain the origins of ToxI antitoxin selectivity through our crystal structure of the ToxIN(Bt) complex. Our results show how a processed RNA pseudoknot can inhibit a deleterious protein with exquisite molecular specificity and how these self-contained and addictive RNA-protein pairs can confer different adaptive benefits in their bacterial hosts.

A common origin for the bacterial toxin-antitoxin systems parD and ccd, suggested by analyses of toxin/target and toxin/antitoxin interactions

Bacterial toxin-antitoxin (TA) systems encode two proteins, a potent inhibitor of cell proliferation (toxin) and its specific antidote (antitoxin). Structural data has revealed striking similarities between the two model TA toxins CcdB, a DNA gyrase inhibitor encoded by the ccd system of plasmid F, and Kid, a site-specific endoribonuclease encoded by the parD system of plasmid R1. While a common structural fold seemed at odds with the two clearly different modes of action of these toxins, the possibility of functional crosstalk between the parD and ccd systems, which would further point to their common evolutionary origin, has not been documented. Here, we show that the cleavage of RNA and the inhibition of protein synthesis by the Kid toxin, two activities that are specifically counteracted by its cognate Kis antitoxin, are altered, but not inhibited, by the CcdA antitoxin. In addition, Kis was able to inhibit the stimulation of DNA gyrase-mediated cleavage of DNA by CcdB, albeit less efficiently than CcdA. We further show that physical interactions between the toxins and antitoxins of the different systems do occur and define the stoichiometry of the complexes formed. We found that CcdB did not degrade RNA nor did Kid have any reproducible effect on the tested DNA gyrase activities, suggesting that these toxins evolved to reach different, rather than common, cellular targets.

Toxin-antitoxin systems in bacteria and archaea

Almost all bacteria and many archaea contain genes whose expression inhibits cell growth and may lead to cell death when overproduced, reminiscent of apoptotic genes in higher systems. The cellular targets of these toxins are quite diverse and include DNA replication, mRNA stability, protein synthesis, cell-wall biosynthesis, and ATP synthesis. These toxins are co-expressed and neutralized with their cognate antitoxins from a TA (toxin-antitoxin) operon in normally growing cells. Antitoxins are more labile than toxins and are readily degraded under stress conditions, allowing the toxins to exert their toxic effect. Presence of at least 33 TA systems in Escherichia coli and more than 60 TA systems in Mycobacterium tuberculosis suggests that the TA systems are involved not only in normal bacterial physiology but also in pathogenicity of bacteria. The elucidation of their cellular function and regulation is thus crucial for our understanding of bacterial physiology under various stress conditions.

Kis antitoxin couples plasmid R1 replication and parD (kis,kid) maintenance modules

The coupling between the replication and parD (kis, kid) maintenance modules of R1 has been revisited here by the isolation of a significant collection of conditional replication mutants in the pKN1562 mini-R1 plasmid, and in its derivative, pJLV01, specifically affected in the RNase activity of the Kid toxin. This new analysis aims to identify key factors in this coupling. For this purpose we have quantified and characterized the restriction introduced by parD to isolate conditional replication mutants of this plasmid, a signature of the modular coupling. This restriction depends on the RNase activity of the Kid toxin and it is relieved by either over-expression of the Kis antitoxin or by preventing its degradation by Lon and ClpAP proteases. Based on these data and on the correlation between copy numbers and parD transcriptional levels obtained in the different mutants, it is proposed that a reduction of Kis antitoxin levels in response to inefficient plasmid replication is the key factor for coupling plasmid replication and parD modules.

Nucleotide docking: prediction of reactant state complexes for ribonuclease enzymes

Ribonuclease enzymes (RNases) play key roles in the maturation and metabolism of all RNA molecules. Computational simulations of the processes involved can help to elucidate the underlying enzymatic mechanism and is often employed in a synergistic approach together with biochemical experiments. Theoretical calculations require atomistic details regarding the starting geometries of the molecules involved, which, in the absence of crystallographic data, can only be achieved from computational docking studies. Fortunately, docking algorithms have improved tremendously in recent years, so that reliable structures of enzyme-ligand complexes can now be successfully obtained from computation. However, most docking programs are not particularly optimized for nucleotide docking. In order to assist our studies on the cleavage of RNA by the two most important ribonuclease enzymes, RNase A and RNase H, we evaluated four docking tools-MOE2009, Glide 5.5, QXP-Flo+0802, and Autodock 4.0-for their ability to simulate complexes between these enzymes and RNA oligomers. To validate our results, we analyzed the docking results with respect to the known key interactions between the protein and the nucleotide. In addition, we compared the predicted complexes with X-ray structures of the mutated enzyme as well as with structures obtained from previous calculations. In this manner, we were able to prepare the desired reaction state complex so that it could be used as the starting structure for further DFT/B3LYP QM/MM reaction mechanism studies.

Crystal structures of Phd-Doc, HigA, and YeeU establish multiple evolutionary links between microbial growth-regulating toxin-antitoxin systems

Bacterial toxin-antitoxin (TA) systems serve a variety of physiological functions including regulation of cell growth and maintenance of foreign genetic elements. Sequence analyses suggest that TA families are linked by complex evolutionary relationships reflecting likely swapping of functional domains between different TA families. Our crystal structures of Phd-Doc from bacteriophage P1, the HigA antitoxin from Escherichia coli CFT073, and YeeU of the YeeUWV systems from E. coli K12 and Shigella flexneri confirm this inference and reveal additional, unanticipated structural relationships. The growth-regulating Doc toxin exhibits structural similarity to secreted virulence factors that are toxic for eukaryotic target cells. The Phd antitoxin possesses the same fold as both the YefM and NE2111 antitoxins that inhibit structurally unrelated toxins. YeeU, which has an antitoxin-like activity that represses toxin expression, is structurally similar to the ribosome-interacting toxins YoeB and RelE. These observations suggest extensive functional exchanges have occurred between TA systems during bacterial evolution.

The relBE2Spn toxin-antitoxin system of Streptococcus pneumoniae: role in antibiotic tolerance and functional conservation in clinical isolates

Type II (proteic) chromosomal toxin-antitoxin systems (TAS) are widespread in Bacteria and Archaea but their precise function is known only for a limited number of them. Out of the many TAS described, the relBE family is one of the most abundant, being present in the three first sequenced strains of Streptococcus pneumoniae (D39, TIGR4 and R6). To address the function of the pneumococcal relBE2Spn TAS in the bacterial physiology, we have compared the response of the R6-relBE2Spn wild type strain with that of an isogenic derivative, Delta relB2Spn under different stress conditions such as carbon and amino acid starvation and antibiotic exposure. Differences on viability between the wild type and mutant strains were found only when treatment directly impaired protein synthesis. As a criterion for the permanence of this locus in a variety of clinical strains, we checked whether the relBE2Spn locus was conserved in around 100 pneumococcal strains, including clinical isolates and strains with known genomes. All strains, although having various types of polymorphisms at the vicinity of the TA region, contained a functional relBE2Spn locus and the type of its structure correlated with the multilocus sequence type. Functionality of this TAS was maintained even in cases where severe rearrangements around the relBE2Spn region were found. We conclude that even though the relBE2Spn TAS is not essential for pneumococcus, it may provide additional advantages to the bacteria for colonization and/or infection.

Structural and thermodynamic characterization of Vibrio fischeri CcdB

CcdB(Vfi) from Vibrio fischeri is a member of the CcdB family of toxins that poison covalent gyrase-DNA complexes. In solution CcdB(Vfi) is a dimer that unfolds to the corresponding monomeric components in a two-state fashion. In the unfolded state, the monomer retains a partial secondary structure. This observation correlates well with the crystal and NMR structures of the protein, which show a dimer with a hydrophobic core crossing the dimer interface. In contrast to its F plasmid homologue, CcdB(Vfi) possesses a rigid dimer interface, and the apparent relative rotations of the two subunits are due to structural plasticity of the monomer. CcdB(Vfi) shows a number of non-conservative substitutions compared with the F plasmid protein in both the CcdA and the gyrase binding sites. Although variation in the CcdA interaction site likely determines toxin-antitoxin specificity, substitutions in the gyrase-interacting region may have more profound functional implications.

PemK toxin of Bacillus anthracis is a ribonuclease: an insight into its active site, structure, and function

Bacillus anthracis genome harbors a toxin-antitoxin (TA) module encoding pemI (antitoxin) and pemK (toxin). This study describes the rPemK as a potent ribonuclease with a preference for pyrimidines (C/U), which is consistent with our previous study that demonstrated it as a translational attenuator. The in silico structural modeling of the PemK in conjunction with the site-directed mutagenesis confirmed the role of His-59 and Glu-78 as an acid-base couple in mediating the ribonuclease activity. The rPemK is shown to form a complex with the rPemI, which is in line with its function as a TA module. This rPemI-rPemK complex becomes catalytically inactive when both the proteins interact in a molar stoichiometry of 1. The rPemI displays vulnerability to proteolysis but attains conformational stability only upon rPemK interaction. The pemI-pemK transcript is shown to be up-regulated upon stress induction with a concomitant increase in the amount of PemK and a decline in the PemI levels, establishing the role of these modules in stress. The artificial perturbation of TA interaction could unleash the toxin, executing bacterial cell death. Toward this end, synthetic peptides are designed to disrupt the TA interaction. The peptides are shown to be effective in abrogating TA interaction in micromolar range in vitro. This approach can be harnessed as a potential antibacterial strategy against anthrax in the future.

A mutagenic analysis of the RNase mechanism of the bacterial Kid toxin by mass spectrometry

Kid, the toxin of the parD (kis, kid) maintenance system of plasmid R1, is an endoribonuclease that preferentially cleaves RNA at the 5' of A in the core sequence 5'-UA(A/C)-3'. A model of the Kid toxin interacting with the uncleavable mimetic 5'-AdUACA-3' is available. To evaluate this model, a significant collection of mutants in some of the key residues proposed to be involved in RNA binding (T46, A55, T69 and R85) or RNA cleavage (R73, D75 and H17) were analysed by mass spectrometry in RNA binding and cleavage assays. A pair of substrates, 5'-AUACA-3', and its uncleavable mimetic 5'-AdUACA-3', used to establish the model and structure of the Kid-RNA complex, were used in both the RNA cleavage and binding assays. A second RNA substrate, 5'-UUACU-3' efficiently cleaved by Kid both in vivo and in vitro, was also used in the cleavage assays. Compared with the wild-type protein, mutations in the residues of the catalytic site abolished RNA cleavage without substantially altering RNA binding. Mutations in residues proposed to be involved in RNA binding show reduced binding efficiency and a corresponding decrease in RNA cleavage efficiency. The cleavage profiles of the different mutants were similar with the two substrates used, but RNA cleavage required much lower protein concentrations when the 5'-UUACU-3' substrate was used. Protein synthesis and growth assays are consistent with there being a correlation between the RNase activity of Kid and its inhibitory potential. These results give important support to the available models of Kid RNase and the Kid-RNA complex.

Significant bias against the ACA triplet in the tmRNA sequence of Escherichia coli K-12

The toxin MazF in Escherichia coli cleaves single-stranded RNAs specifically at ACA sequences. MazF overexpression virtually eliminates all cellular mRNAs to completely block protein synthesis. However, protein synthesis can continue on an mRNA that is devoid of ACA triplets. The finding that ribosomal RNAs remain intact in the face of complete translation arrest suggested a purpose for such preservation. We therefore examined the sequences of all transcribed RNAs to determine if there was any statistically significant bias against ACA. While ACA motifs are absent from tmRNA, 4.5S RNA, and seven of the eight 5S rRNAs, statistical analysis revealed that only for tmRNA was the absence nonrandom. The introduction of single-strand ACAs makes tmRNA highly susceptible to MazF cleavage. Furthermore, analysis of tmRNA sequences from 442 bacteria showed that the discrimination against ACA in tmRNAs was seen mostly in enterobacteria. We propose that the unusual bias against ACA in tmRNA may have coevolved with the acquisition of MazF.

Insights into the DNA cleavage mechanism of human LINE-1 retrotransposon endonuclease

The human LINE-1 endonuclease (L1-EN) contributes in defining the genomic integration sites of the abundant human L1 and Alu retrotransposons. LINEs have been considered as possible vehicles for gene delivery and understanding the mechanism of L1-EN could help engineering them as genetic tools. We tested the in vitro activity of point mutants in three L1-EN residues--Asp145, Arg155, Ile204--that are key for DNA cleavage, and determined their crystal structures. The L1-EN structure remains overall unaffected by the mutations, which change the enzyme activity but leave DNA cleavage sequence specificity mostly unaffected. To better understand the mechanism of L1-EN, we performed molecular dynamics simulations using as model the structures of wild type EN-L1, of two betaB6-betaB5 loop exchange mutants we have described previously to be important for DNA recognition, of the R155A mutant from this study, and of the homologous TRAS1 endonuclease: all confirm a rigid scaffold. The simulations crucially indicate that the betaB6-betaB5 loop shows an anticorrelated motion with the surface loops betaA6-betaA5 and betaB3-alphaB1. The latter loop harbors N118, a residue that alters DNA cleavage specificity in homologous endonucleases, and implies that the plasticity and correlated motion of these loops has a functional importance in DNA recognition and binding. To further explore how these loops are possibly involved in DNA binding, we docked computationally two DNA substrates to our structure, one involving a flipped-out nucleotide downstream the scissile phosphodiester; and one not. The models for both scenarios are feasible and agree with the hypotheses derived from the dynamic simulations. The reduced cleavage activity we have observed for the I204Y mutant above however, favors the flipped out nucleotide model.

Structural analysis of MED-1 reveals unexpected diversity in the mechanism of DNA recognition by GATA-type zinc finger domains

MED-1 is a member of a group of divergent GATA-type zinc finger proteins recently identified in several species of Caenorhabditis. The med genes are transcriptional regulators that are involved in the specification of the mesoderm and endoderm precursor cells in nematodes. Unlike other GATA-type zinc fingers that recognize the consensus sequence (A/C/T)GATA(A/G), the MED-1 zinc finger (MED1zf) binds the larger and atypical site GTATACT(T/C)(3). We have examined the basis for this unusual DNA specificity using a range of biochemical and biophysical approaches. Most strikingly, we show that although the core of the MED1zf structure is similar to that of GATA-1, the basic tail C-terminal to the zinc finger unexpectedly adopts an alpha-helical structure upon binding DNA. This additional helix appears to contact the major groove of the DNA, making contacts that explain the extended DNA consensus sequence observed for MED1zf. Our data expand the versatility of DNA recognition by GATA-type zinc fingers and perhaps shed new light on the DNA-binding properties of mammalian GATA factors.

Influence of operator site geometry on transcriptional control by the YefM-YoeB toxin-antitoxin complex

YefM-YoeB is among the most prevalent and well-characterized toxin-antitoxin complexes. YoeB toxin is an endoribonuclease whose activity is inhibited by YefM antitoxin. The regions 5' of yefM-yoeB in diverse bacteria possess conserved sequence motifs that mediate transcriptional autorepression. The yefM-yoeB operator site arrangement is exemplified in Escherichia coli: a pair of palindromes with core hexamer motifs and a center-to-center distance of 12 bp overlap the yefM-yoeB promoter. YefM is an autorepressor that initially recognizes a long palindrome containing the core hexamer, followed by binding to a short repeat. YoeB corepressor greatly enhances the YefM-operator interaction. Scanning mutagenesis demonstrated that the short repeat is crucial for correct interaction of YefM-YoeB with the operator site in vivo and in vitro. Moreover, altering the relative positions of the two palindromes on the DNA helix abrogated YefM-YoeB cooperative interactions with the repeats: complex binding to the long repeat was maintained but was perturbed to the short repeat. Although YefM lacks a canonical DNA binding motif, dual conserved arginine residues embedded in a basic patch of the protein are crucial for operator recognition. Deciphering the molecular basis of toxin-antitoxin transcriptional control will provide key insights into toxin-antitoxin activation and function.

Novel Escherichia coli RF1 mutants with decreased translation termination activity and increased sensitivity to the cytotoxic effect of the bacterial toxins Kid and RelE

Novel mutations in prfA, the gene for the polypeptide release factor RF1 of Escherichia coli, were isolated using a positive genetic screen based on the parD (kis, kid) toxin–antitoxin system. This original approach allowed the direct selection of mutants with altered translational termination efficiency at UAG codons. The isolated prfA mutants displayed a ∼10-fold decrease in UAG termination efficiency with no significant changes in RF1 stability in vivo. All three mutations, G121S, G301S and R303H, were situated close to the nonsense codon recognition site in RF1:ribosome complexes. The prfA mutants displayed increased sensitivity to the RelE toxin encoded by the relBE system of E. coli, thus providing in vivo support for the functional interaction between RF1 and RelE. The prfA mutants also showed increased sensitivity to the Kid toxin. Since this toxin can cleave RNA in a ribosome-independent manner, this result was not anticipated and provided first evidence for the involvement of RF1 in the pathway of Kid toxicity. The sensitivity of the prfA mutants to RelE and Kid was restored to normal levels upon overproduction of the wild-type RF1 protein. We discuss these results and their utility for the design of novel antibacterial strategies in the light of the recently reported structure of ribosome-bound RF1.

Towards deciphering the principles underlying an mRNA recognition code

Messenger RNAs interact with a number of different molecules that determine the fate of each transcript and contribute to the overall pattern of gene expression. These interactions are governed by specific mRNA signals, which in principle could represent a special mRNA recognition 'code'. Both, small molecules and proteins demonstrate a diversity of mRNA binding modes often dependent on the structural context of the regions surrounding specific target sequences. In this review, we have highlighted recent structural studies that illustrate the diversity of recognition principles used by mRNA binders for timely and specific targeting and processing of the message.

Structural and biophysical analysis of the DNA binding properties of myelin transcription factor 1

Zinc binding domains, or zinc fingers (ZnFs), form one of the most numerous and most diverse superclasses of protein structural motifs in eukaryotes. Although our understanding of the functions of several classes of these domains is relatively well developed, we know much less about the molecular mechanisms of action of many others. Myelin transcription factor 1 (MyT1) type ZnFs are found in organisms as diverse as nematodes and mammals and are found in a range of sequence contexts. MyT1, one of the early transcription factors expressed in the developing central nervous system, contains seven MyT1 ZnFs that are very highly conserved both within the protein and between species. We have used a range of biophysical techniques, including NMR spectroscopy and data-driven macromolecular docking, to investigate the structural basis for the interaction between MyT1 ZnFs and DNA. Our data indicate that MyT1 ZnFs recognize the major groove of DNA in a way that appears to differ from other known zinc binding domains.

Interactions between the toxin Kid of the bacterial parD system and the antitoxins Kis and MazE

The proteins Kid and Kis are the toxin and antitoxin, respectively, encoded by the parD operon of Escherichia coli plasmid R1. Kis prevents the inhibition of E. coli cell growth caused by the RNA cleavage activity of Kid. Overproduction of MazE, the chromosome-encoded homologue of Kis, has been demonstrated to neutralize Kid toxicity to a certain extent in the absence of native Kis. Here, we show that a high structural similarity exists between these antitoxins, using NMR spectroscopy. We report about the interactions between Kid and Kis that are responsible for neutralization of Kid toxicity and enhance autoregulation of parD transcription. Native macromolecular mass spectrometry data demonstrate that Kid and Kis form multiple complexes. At Kis:Kid ratios equal to or exceeding 1:1, as found in vivo in a plasmid-containing cell, various complexes are present, ranging from Kid(2)-Kis(2) tetramer up to Kis(2)-Kid(2)-Kis(2)-Kid(2)-Kis(2) decamer. When Kid is in excess of Kis, corresponding to an in vivo situation immediately after loss of the plasmid, the Kid(2)-Kis(2)-Kid(2) heterohexamer is the most abundant species. NMR chemical shift and intensity perturbations in the (1)H (15)N HSQC spectra of Kid and Kis, observed when titrating the partner protein, show that the interaction sites of Kid and Kis resemble those within the previously reported MazF(2)-MazE(2)-MazF(2) complex. Furthermore, we demonstrate that Kid(2)-MazE(2) tetramers can be formed via weak interactions involving a limited part of the Kis-binding residues of Kid. The functional roles of the identified Kid-Kis and Kid-MazE interaction sites and complexes in toxin neutralization and repression of transcription are discussed.

An archaeal protein with homology to the eukaryotic translation initiation factor 5A shows ribonucleolytic activity

To identify proteins that are involved in RNA degradation and processing in Halobacterium sp. NRC-1, we purified proteins with RNA-degrading activity by classical biochemical techniques. One of these proteins showed strong homology to the eukaryotic initiation factor 5A (eIF-5A) and was accordingly named archaeal initiation factor 5A (aIF-5A). Eukaryotic IF-5A is known to be involved in mRNA turnover and to bind RNA. Hypusination of eIF-5A is required for sequence-specific binding of RNA. This unique post-translational modification is restricted to Eukarya and Archaea. The exact function of eIF-5A in RNA turnover remained obscure. Here we show for the first time that aIF-5A from Halobacterium sp. NRC-1 exhibits RNA cleavage activity, preferentially cleaving adjacent to A nucleotides. Detectable RNA binding could be shown for aIF-5A purified from Halobacterium sp. NRC-1 but not from Escherichia coli, while both proteins possess RNA cleavage activity, indicating that hypusination of aIF-5A is required for RNA binding but not for its RNA cleavage activity. Furthermore, we show that the hypusinated form of eIF-5A also shows RNase activity while the unmodified protein does not. Charged amino acids in the N-terminal domain of aIF-5A as well as in the C-terminal domain, which is highly similar to the cold shock protein A (CspA), an RNA chaperone of E. coli, are important for RNA cleavage activity. Moreover our results reveal that activity of aIF-5A depends on its oligomeric state.

Interactions of Kid-Kis toxin-antitoxin complexes with the parD operator-promoter region of plasmid R1 are piloted by the Kis antitoxin and tuned by the stoichiometry of Kid-Kis oligomers

The parD operon of Escherichia coli plasmid R1 encodes a toxin–antitoxin system, which is involved in plasmid stabilization. The toxin Kid inhibits cell growth by RNA degradation and its action is neutralized by the formation of a tight complex with the antitoxin Kis. A fascinating but poorly understood aspect of the kid–kis system is its autoregulation at the transcriptional level. Using macromolecular (tandem) mass spectrometry and DNA binding assays, we here demonstrate that Kis pilots the interaction of the Kid–Kis complex in the parD regulatory region and that two discrete Kis-binding regions are present on parD. The data clearly show that only when the Kis concentration equals or exceeds the Kid concentration a strong cooperative effect exists between strong DNA binding and Kid₂–Kis₂–Kid₂–Kis₂ complex formation. We propose a model in which transcriptional repression of the parD operon is tuned by the relative molar ratio of the antitoxin and toxin proteins in solution. When the concentration of the toxin exceeds that of the antitoxin tight Kid₂–Kis₂–Kid₂ complexes are formed, which only neutralize the lethal activity of Kid. Upon increasing the Kis concentration, (Kid₂–Kis₂)_n complexes repress the kid–kis operon.

Fungal ribotoxins: molecular dissection of a family of natural killers

RNase T1 is the best known representative of a large family of ribonucleolytic proteins secreted by fungi, mostly Aspergillus and Penicillium species. Ribotoxins stand out among them by their cytotoxic character. They exert their toxic action by first entering the cells and then cleaving a single phosphodiester bond located within a universally conserved sequence of the large rRNA gene, known as the sarcin-ricin loop. This cleavage leads to inhibition of protein biosynthesis, followed by cellular death by apoptosis. Although no protein receptor has been found for ribotoxins, they preferentially kill cells showing altered membrane permeability, such as those that are infected with virus or transformed. Many steps of the cytotoxic process have been elucidated at the molecular level by means of a variety of methodological approaches and the construction and purification of different mutant versions of these ribotoxins. Ribotoxins have been used for the construction of immunotoxins, because of their cytotoxicity. Besides this activity, Aspf1, a ribotoxin produced by Aspergillus fumigatus, has been shown to be one of the major allergens involved in allergic aspergillosis-related pathologies. Protein engineering and peptide synthesis have been used in order to understand the basis of these pathogenic mechanisms as well as to produce hypoallergenic proteins with potential diagnostic and immunotherapeutic applications.

Toxin-antitoxin systems are ubiquitous and plasmid-encoded in vancomycin-resistant enterococci

Vancomycin-resistant enterococci (VRE) are common hospital pathogens that are resistant to most major classes of antibiotics. The incidence of VRE is increasing rapidly, to the point where over one-quarter of enterococcal infections in intensive care units are now resistant to vancomycin. The exact mechanism by which VRE maintains its plasmid-encoded resistance genes is ill-defined, and novel targets for the treatment of VRE are lacking. In an effort to identify novel protein targets for the treatment of VRE infections, we probed the plasmids obtained from 75 VRE isolates for the presence of toxin-antitoxin (TA) gene systems. Remarkably, genes for one particular TA pair, the mazEF system (originally identified on the Escherichia coli chromosome), were present on plasmids from 75/75 (100%) of the isolates. Furthermore, mazEF was on the same plasmid as vanA in the vast majority of cases (>90%). Plasmid stability tests and RT-PCR raise the possibility that this plasmid-encoded mazEF is indeed functional in enterococci. Given this ubiquity of mazEF in VRE and the deleterious activity of the MazF toxin, disruption of mazEF with pharmacological agents is an attractive strategy for tailored antimicrobial therapy.

Information-driven protein-DNA docking using HADDOCK: it is a matter of flexibility

Intrinsic flexibility of DNA has hampered the development of efficient protein−DNA docking methods. In this study we extend HADDOCK (High Ambiguity Driven DOCKing) [C. Dominguez, R. Boelens and A. M. J. J. Bonvin (2003) J. Am. Chem. Soc.125, 1731–1737] to explicitly deal with DNA flexibility. HADDOCK uses non-structural experimental data to drive the docking during a rigid-body energy minimization, and semi-flexible and water refinement stages. The latter allow for flexibility of all DNA nucleotides and the residues of the protein at the predicted interface. We evaluated our approach on the monomeric repressor−DNA complexes formed by bacteriophage 434 Cro, the Escherichia coli Lac headpiece and bacteriophage P22 Arc. Starting from unbound proteins and canonical B-DNA we correctly predict the correct spatial disposition of the complexes and the specific conformation of the DNA in the published complexes. This information is subsequently used to generate a library of pre-bent and twisted DNA structures that served as input for a second docking round. The resulting top ranking solutions exhibit high similarity to the published complexes in terms of root mean square deviations, intermolecular contacts and DNA conformation. Our two-stage docking method is thus able to successfully predict protein−DNA complexes from unbound constituents using non-structural experimental data to drive the docking.

Shutdown decay of mRNA

Although plasmid-borne and chromosomal toxin-antitoxin (TA) operons have been known for some time, the recent identification of mRNA as the target of at least two different classes of toxins has led to a dramatic renewal of interest in these systems as mediators of stress responses. Members of the MazF/PemK family, the so-called mRNA interferases, are ribonucleases that inhibit translation by destroying cellular mRNAs under stress conditions, while the founder member of the RelE family promotes cleavage of mRNAs through the ribosome. Detailed structures of these enzymes, often in complex with their inhibitors, have provided vital clues to their mechanisms of action. The primary role and regulation of these systems has been the subject of some controversy. One model suggests they play a beneficial role by wiping the slate clean and preventing wasteful energy consumption by the translational apparatus during adaptation to stress conditions, while another favours the idea that their main function is programmed cell death. The two models might not be mutually exclusive if a side-effect of prolonged exposure to toxic RNase activity without de novo synthesis of the inhibitor were a state of dormancy for which we do not yet understand the key to recovery. In this review, I discuss the recent developments in the rapidly expanding field of what I refer to as bacterial shutdown decay.