Conditional cooperativity of toxin - antitoxin regulation can mediate bistability between growth and dormancy

Quantifying Dynamic Phenotypic Heterogeneity in Resistant Escherichia coli under Translation-Inhibiting Antibiotics

Understanding the phenotypic heterogeneity of antibiotic-resistant bacteria following treatment and the transitions between different phenotypes is crucial for developing effective infection control strategies. The study expands upon previous work by explicating chloramphenicol-induced phenotypic heterogeneities in growth rate, gene expression, and morphology of resistant Escherichia coli using time-lapse microscopy. Correlating the bacterial growth rate and cspC expression, four interchangeable phenotypic subpopulations across varying antibiotic concentrations are identified, surpassing the previously described growth rate bistability. Notably, bacterial cells exhibiting either fast or slow growth rates can concurrently harbor subpopulations characterized by high and low gene expression levels, respectively. To elucidate the mechanisms behind this enhanced heterogeneity, a concise gene expression network model is proposed and the biological significance of the four phenotypes is further explored. Additionally, by employing Hidden Markov Model fitting and integrating the non-equilibrium landscape and flux theory, the real-time data encompassing diverse bacterial traits are analyzed. This approach reveals dynamic changes and switching kinetics in different cell fates, facilitating the quantification of observable behaviors and the non-equilibrium dynamics and thermodynamics at play. The results highlight the multi-dimensional heterogeneous behaviors of antibiotic-resistant bacteria under antibiotic stress, providing new insights into the compromised antibiotic efficacy, microbial response, and associated evolution processes.

Integrative biology of persister cell formation: molecular circuitry, phenotypic diversification and fitness effects

Microbial populations often contain persister cells, which reduce the extinction risk upon sudden stresses. Persister cell formation is deeply intertwined with physiology. Due to this complexity, it cannot be satisfactorily understood by focusing only on mechanistic, physiological or evolutionary aspects. In this review, we take an integrative biology perspective to identify common principles of persister cell formation, which might be applicable across evolutionary-distinct microbes. Persister cells probably evolved to cope with a fundamental trade-off between cellular stress and growth tasks, as any biosynthetic resource investment in growth-supporting proteins is at the expense of stress tasks and vice versa. Natural selection probably favours persister cell subpopulation formation over a single-phenotype strategy, where each cell is prepared for growth and stress to a suboptimal extent, since persister cells can withstand harsher environments and their coexistence with growing cells leads to a higher fitness. The formation of coexisting phenotypes requires bistable molecular circuitry. Bistability probably emerges from growth-modulated, positive feedback loops in the cell's growth versus stress control network, involving interactions between sigma factors, guanosine pentaphosphate and toxin-antitoxin (TA) systems. We conclude that persister cell formation is most likely a response to a sudden reduction in growth rate, which can be achieved by antibiotic addition, nutrient starvation, sudden stresses, nutrient transitions or activation of a TA system.

Bacterial toxin-antitoxin modules: classification, functions, and association with persistence

•

Ubiquitously present bacterial Toxin-Antitoxin (TA) modules consist of stable toxin associated with labile antitoxin.

•

Classification of TAs modules based on inhibition of toxin through antitoxin in 8 different classes.

•

Variety of specific toxin targets and the abundance of TA modules in various deadly pathogens.

•

Specific role of TAs modules in conservation of the resistant genes, emergence of persistence & biofilm formation.

•

Proposed antibacterial strategies involving TA modules for elimination of multi-drug resistance.

Toxin-antitoxin (TA) modules are ubiquitous gene loci among bacteria and are comprised of a toxin part and its cognate antitoxin part. Under normal physiological conditions, antitoxin counteracts the toxicity of the toxin whereas, during stress conditions, TA modules play a crucial role in bacterial physiology through involvement in the post-segregational killing, abortive infection, biofilms, and persister cell formation. Most of the toxins are proteinaceous that affect translation or DNA replication, although some other intracellular molecular targets have also been described. While antitoxins may be a protein or RNA, that generally neutralizes its cognate toxin by direct interaction or with the help of other signaling elements and thus helps in the TA module regulation. In this review, we have discussed the current state of the multifaceted TA (type I–VIII) modules by highlighting their classification and specific targets. We have also discussed the presence of TA modules in the various pathogens and their role in antibiotic persistence development as well as biofilm formation, by influencing the different cellular processes. In the end, assembling knowledge about ubiquitous TA systems from pathogenic bacteria facilitated us to propose multiple novel antibacterial strategies involving artificial activation of TA modules.

Multiple Copies of flhDC in Paraburkholderia unamae Regulate Flagellar Gene Expression, Motility, and Biofilm Formation

FlhDC is a heterohexameric complex that acts as a master regulator of flagellar biosynthesis genes in numerous bacteria. Previous studies have identified a single flhDC operon encoding this complex. However, we found that two flhDC loci are present throughout Paraburkholderia, and two additional flhC copies are also present in Paraburkholderia unamae. Systematic deletion analysis in P. unamae of the different flhDC copies showed that one of the operons, flhDC1, plays the predominant role, with deletion of its genes resulting in a severe inhibition of motility and biofilm formation. Expression analysis using promoter-lacZ fusions and real-time quantitative PCR support the primary role of flhDC1 in flagellar gene regulation, with flhDC2 a secondary contributor. Phylogenetic analysis shows the presence of the flhDC1 and flhDC2 operons throughout Paraburkholderia. In contrast, Burkholderia and other bacteria only carry the copy syntenous with flhDC2. The variations in impact each copy of flhDC has on downstream processes indicate that regulation of FlhDC in P. unamae, and likely other Paraburkholderia species, is regulated at least in part by the presence of multiple copies of these genes. IMPORTANCE Motility is important in the colonization of plant roots by beneficial and pathogenic bacteria, with flagella playing essential roles in host cell adhesion, entrance, and biofilm formation. Flagellar biosynthesis is energetically expensive. Its complex regulation by the FlhDC master regulator is well studied in peritrichous flagella expressing enterics. We report the unique presence throughout Paraburkholderia of multiple copies of flhDC. In P. unamae, the flhDC1 copy showed higher expression and a greater effect on swim motility, flagellar development, and regulation of downstream genes, than the flhDC2 copy that is syntenous to flhDC in Escherichia coli and pathogenic Burkholderia spp. The flhDC genes have evolved differently in these plant-growth-promoting bacteria, giving an additional layer of complexity in gene regulation by FlhDC.

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.

Impacts of type II toxin-antitoxin systems on cell physiology and environmental behavior in acetic acid bacteria

Acetic acid bacteria (AAB) are a group of Gram-negative and strictly aerobic microorganisms widely used in vinegar industry, especially the species belonging to the genera Acetobacter and Komagataeibacter. The environments inhabited by AAB during the vinegar fermentation, in particular those natural traditional bioprocesses, are complex and dynamically changed, usually accompanied by diverse microorganisms, bacteriophages, and the increasing acetic acid concentration. For this reason, how AAB survive to such harsh niches has always been an interesting research field. Previous omic analyses (e.g., genomics, proteomics, and transcriptomics) have provided abundant clues for the metabolic pathways and bioprocesses indispensable for the acid stress adaptation of AAB. Nevertheless, it is far from fully understanding what factors regulate these modular mechanisms overtly and covertly upon shifting environments. Bacterial toxin-antitoxin systems (TAS), usually consisting of a pair of genes encoding a stable toxin and an unstable antitoxin that is capable of counteracting the toxin, have been uncovered to have a variety of biological functions. Recent studies focusing on the role of TAS in Acetobacter pasteurianus suggest that TAS contribute substantially to the acid stress resistance. In this mini review, we discuss the biological functions of type II TAS in the context of AAB with regard to the acid stress resistance, persister formation and resuscitation, genome stability, and phage immunity. KEY POINTS: • Type II TAS act as regulators in the acid stress resistance of AAB. • Type II TAS are implicated in the formation of acid-tolerant persister cells in AAB. • Type II TAS are potential factors responsible for phage immunity and genome stability.

A systematic search for switch-like behavior in type II toxin-antitoxin systems

Bistable switch-like behavior is a ubiquitous feature of gene regulatory networks with decision-making capabilities. Type II toxin-antitoxin (TA) systems are hypothesized to facilitate a bistable switch in toxin concentration that influences the dormancy transition in persister cells. However, a series of recent retractions has raised fundamental questions concerning the exact mechanism of toxin propagation in persister cells and the relationship between type II TA systems and cellular dormancy. Through a careful modeling search, we identify how sp: bistablilty can emerge in type II TA systems by systematically modifying a basic model for the RelBE system with other common biological mechanisms. Our systematic search uncovers a new combination of mechanisms influencing bistability in type II TA systems and explores how toxin bistability emerges through synergistic interactions between paired type II TA systems. Our analysis also illustrates how Descartes' rule of signs and the resultant can be used as a powerful delineator of bistability in mathematical systems regardless of application.

Figure 1 Theory Meets Figure 2 Experiments in the Study of Gene Expression

It is tempting to believe that we now own the genome. The ability to read and rewrite it at will has ushered in a stunning period in the history of science. Nonetheless, there is an Achilles' heel exposed by all of the genomic data that has accrued: We still do not know how to interpret them. Many genes are subject to sophisticated programs of transcriptional regulation, mediated by DNA sequences that harbor binding sites for transcription factors, which can up- or down-regulate gene expression depending upon environmental conditions. This gives rise to an input-output function describing how the level of expression depends upon the parameters of the regulated gene-for instance, on the number and type of binding sites in its regulatory sequence. In recent years, the ability to make precision measurements of expression, coupled with the ability to make increasingly sophisticated theoretical predictions, has enabled an explicit dialogue between theory and experiment that holds the promise of covering this genomic Achilles' heel. The goal is to reach a predictive understanding of transcriptional regulation that makes it possible to calculate gene expression levels from DNA regulatory sequence. This review focuses on the canonical simple repression motif to ask how well the models that have been used to characterize it actually work. We consider a hierarchy of increasingly sophisticated experiments in which the minimal parameter set learned at one level is applied to make quantitative predictions at the next. We show that these careful quantitative dissections provide a template for a predictive understanding of the many more complex regulatory arrangements found across all domains of life.

Non-monotonic auto-regulation in single gene circuits

We theoretically study the effects of non-monotonic response curves in genetic auto-regulation by exploring the possible dynamical behaviors for such systems. Our motivation is twofold: we aim at conceiving the simplest genetic circuits for synthetic biology and at understanding the natural auto-regulation of the LrpB protein of the Sulfolobus solfataricus archaeon which exhibits non-monotonicity. We analyzed three toy models, based on mass-action kinetics, with increasing complexity and sought for oscillations and (fast) bistable switching. We performed large parameter scans and sensitivity analyses, and quantified the quality of the oscillators and switches by computing relative volumes in parameter space reproducing the sought dynamical behavior. All single gene systems need finely tuned parameters in order to oscillate, but bistable switches are more robust against parameter changes. We expected non-monotonic switches to be faster than monotonic ones, however solutions combining both auto-activation and repression in the physiological range to obtain fast switches are scarce. Our analysis shows that the Ss-LrpB system can not provide a bistable switch and that robust oscillations are unlikely. Gillespie simulations suggest that the function of the natural Ss-LrpB system is sensing via a spiking behavior, which is in line with the fact that this protein has a metabolic regulatory function and binds to a ligand.

Birth and Resuscitation of (p)ppGpp Induced Antibiotic Tolerant Persister Cells

Transient antibiotic treatment typically eradicates most sensitive bacteria except a few survivors called persisters. The second messenger (p)ppGpp plays a key role in persister formation in Escherichia coli populations but the underlying mechanisms have remained elusive. In this study we induced (p)ppGpp synthesis by modulating tRNA charging and then directly observed the stochastic appearance, antibiotic tolerance, and resuscitation of persister cells using live microscopy. Different physiological parameters of persister cells as well as their regularly growing ancestors and sisters were continuously monitored using fluorescent reporters. Our results confirmed previous findings that high (p)ppGpp levels are critical for persister formation, but the phenomenon remained strikingly stochastic without any correlation between (p)ppGpp levels and antibiotic tolerance on the single-cell level. We could not confirm previous notions that persisters exhibit markedly low concentrations of intracellular ATP or were linked to post-transcriptional effects of (p)ppGpp through the activation of small genetic elements known as toxin-antitoxin (TA) modules. Instead, we suggest that persister cell formation under regular conditions is driven by the transcriptional response to increased (p)ppGpp levels.

A dual role in regulation and toxicity for the disordered N-terminus of the toxin GraT

Bacterial toxin-antitoxin (TA) modules are tightly regulated to maintain growth in favorable conditions or growth arrest during stress. A typical regulatory strategy involves the antitoxin binding and repressing its own promoter while the toxin often acts as a co-repressor. Here we show that Pseudomonas putida graTA-encoded antitoxin GraA and toxin GraT differ from other TA proteins in the sense that not the antitoxin but the toxin possesses a flexible region. GraA auto-represses the graTA promoter: two GraA dimers bind cooperatively at opposite sides of the operator sequence. Contrary to other TA modules, GraT is a de-repressor of the graTA promoter as its N-terminal disordered segment prevents the binding of the GraT₂A₂ complex to the operator. Removal of this region restores operator binding and abrogates Gr aT toxicity. GraTA represents a TA module where a flexible region in the toxin rather than in the antitoxin controls operon expression and toxin activity.

The Pseudomonas putida toxin GraT and antitoxin GraA form a type II toxin-antoxin module. Here the authors present the crystal structures of the GraA dimer, GraTA and GraA-DNA complexes and show that GraT contains a functionally important N-terminal intrinsic disordered region that prevents the binding of the GraTA complex to the operator.

Excitable dynamics through toxin-induced mRNA cleavage in bacteria

Toxin-antitoxin (TA) systems in bacteria and archaea are small genetic elements consisting of the genes coding for an intracellular toxin and an antitoxin that can neutralize this toxin. In various cases, the toxins cleave the mRNA. In this theoretical work we use deterministic and stochastic modeling to explain how toxin-induced cleavage of mRNA in TA systems can lead to excitability, allowing large transient spikes in toxin levels to be triggered. By using a simplified network where secondary complex formation and transcriptional regulation are not included, we show that a two-dimensional, deterministic model captures the origin of such toxin excitations. Moreover, it allows to increase our understanding by examining the dynamics in the phase plane. By systematically comparing the deterministic results with Gillespie simulations we demonstrate that even though the real TA system is intrinsically stochastic, toxin excitations can be accurately described deterministically. A bifurcation analysis of the system shows that the excitable behavior is due to a nearby Hopf bifurcation in the parameter space, where the system becomes oscillatory. The influence of stress is modeled by varying the degradation rate of the antitoxin and the translation rate of the toxin. We find that stress increases the frequency of toxin excitations. The inclusion of secondary complex formation and transcriptional regulation does not fundamentally change the mechanism of toxin excitations. Finally, we show that including growth rate suppression and translational inhibition can lead to longer excitations, and even cause excitations in cases when the system would otherwise be non-excitable. To conclude, the deterministic model used in this work provides a simple and intuitive explanation of toxin excitations in TA systems.

Autoregulation of bacterial gene expression: lessons from the MazEF toxin-antitoxin system

Autoregulation is the direct modulation of gene expression by the product of the corresponding gene. Autoregulation of bacterial gene expression has been mostly studied at the transcriptional level, when a protein acts as the cognate transcriptional repressor. A recent study investigating dynamics of the bacterial toxin-antitoxin MazEF system has shown how autoregulation at both the transcriptional and post-transcriptional levels affects the heterogeneity of Escherichia coli populations. Toxin-antitoxin systems hold a crucial but still elusive part in bacterial response to stress. This perspective highlights how these modules can also serve as a great model system for investigating basic concepts in gene regulation. However, as the genomic background and environmental conditions substantially influence toxin activation, it is important to study (auto)regulation of toxin-antitoxin systems in well-defined setups as well as in conditions that resemble the environmental niche.

Modeling slow-processing of toxin messenger RNAs in type-I toxin-antitoxin systems: post-segregational killing and noise filtering

In type-I toxin-antitoxin (TA) systems, the action of growth-inhibiting toxin proteins is counteracted by the antitoxin small RNAs (sRNAs) that prevent the translation of toxin messenger RNAs (mRNAs). When a TA module is encoded on a plasmid, the short lifetime of antitoxin sRNA compared to toxin mRNAs mediates post-segregational killing (PSK) that contribute the plasmid maintenance, while some of the chromosomal encoded TA loci have been reported to contribute to persister formation in response to a specific upstream signal. Some of the well studied type-I TA systems such as hok/sok are known to have a rather complex regulatory mechanism. Transcribed full-length toxin mRNAs fold such that the ribosome binding site is not accessible and hence cannot be translated. The mRNAs are slowly processed by RNases, and the truncated mRNAs can be either translated or bound by antitoxin sRNA to be quickly degraded. We analyze the role of this extra processing by a mathematical model. We first consider the PSK scenario, and demonstrate that the extra processing compatibly ensures the high toxin expression upon complete plasmid loss, without inducing toxin expression upon acquisition of a plasmid or decrease of plasmid number to a non-zero number. We further show that the extra processing help filtering the transcription noise, avoiding random activation of toxins in transcriptionally regulated TA systems as seen in chromosomal ones. The present model highlights impacts of the slow processing reaction, offering insights on why the slow processing reactions are commonly identified in multiple type-I TA systems.

Toxins targeting transfer RNAs: Translation inhibition by bacterial toxin-antitoxin systems

Prokaryotic toxin-antitoxin (TA) systems are composed of a protein toxin and its cognate antitoxin. These systems are abundant in bacteria and archaea and play an important role in growth regulation. During favorable growth conditions, the antitoxin neutralizes the toxin's activity. However, during conditions of stress or starvation, the antitoxin is inactivated, freeing the toxin to inhibit growth and resulting in dormancy. One mechanism of growth inhibition used by several TA systems results from targeting transfer RNAs (tRNAs), either through preventing aminoacylation, acetylating the primary amino group, or endonucleolytic cleavage. All of these mechanisms inhibit translation and result in growth arrest. Many of these toxins only act on a specific tRNA or a specific subset of tRNAs; however, more work is necessary to understand the specificity determinants of these toxins. For the toxins whose specificity has been characterized, both sequence and structural components of the tRNA appear important for recognition by the toxin. Questions also remain regarding the mechanisms used by dormant bacteria to resume growth after toxin induction. Rescue of stalled ribosomes by transfer-messenger RNAs, removal of acetylated amino groups from tRNAs, or ligation of cleaved RNA fragments have all been implicated as mechanisms for reversing toxin-induced dormancy. However, the mechanisms of resuming growth after induction of the majority of tRNA targeting toxins are not yet understood. This article is categorized under: Translation > Translation Regulation RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition.

Reassessing the Role of Type II Toxin-Antitoxin Systems in Formation of Escherichia coli Type II Persister Cells

Persistence is a reversible and low-frequency phenomenon allowing a subpopulation of a clonal bacterial population to survive antibiotic treatments. Upon removal of the antibiotic, persister cells resume growth and give rise to viable progeny. Type II toxin-antitoxin (TA) systems were assumed to play a key role in the formation of persister cells in Escherichia coli based on the observation that successive deletions of TA systems decreased persistence frequency. In addition, the model proposed that stochastic fluctuations of (p)ppGpp levels are the basis for triggering activation of TA systems. Cells in which TA systems are activated are thought to enter a dormancy state and therefore survive the antibiotic treatment. Using independently constructed strains and newly designed fluorescent reporters, we reassessed the roles of TA modules in persistence both at the population and single-cell levels. Our data confirm that the deletion of 10 TA systems does not affect persistence to ofloxacin or ampicillin. Moreover, microfluidic experiments performed with a strain reporting the induction of the yefM-yoeB TA system allowed the observation of a small number of type II persister cells that resume growth after removal of ampicillin. However, we were unable to establish a correlation between high fluorescence and persistence, since the fluorescence of persister cells was comparable to that of the bulk of the population and none of the cells showing high fluorescence were able to resume growth upon removal of the antibiotic. Altogether, these data show that there is no direct link between induction of TA systems and persistence to antibiotics.IMPORTANCE Within a growing bacterial population, a small subpopulation of cells is able to survive antibiotic treatment by entering a transient state of dormancy referred to as persistence. Persistence is thought to be the cause of relapsing bacterial infections and is a major public health concern. Type II toxin-antitoxin systems are small modules composed of a toxic protein and an antitoxin protein counteracting the toxin activity. These systems were thought to be pivotal players in persistence until recent developments in the field. Our results demonstrate that previous influential reports had technical flaws and that there is no direct link between induction of TA systems and persistence to antibiotics.

Autoregulation of mazEF expression underlies growth heterogeneity in bacterial populations

The MazF toxin sequence-specifically cleaves single-stranded RNA upon various stressful conditions, and it is activated as a part of the mazEF toxin-antitoxin module in Escherichia coli. Although autoregulation of mazEF expression through the MazE antitoxin-dependent transcriptional repression has been biochemically characterized, less is known about post-transcriptional autoregulation, as well as how both of these autoregulatory features affect growth of single cells during conditions that promote MazF production. Here, we demonstrate post-transcriptional autoregulation of mazF expression dynamics by MazF cleaving its own transcript. Single-cell analyses of bacterial populations during ectopic MazF production indicated that two-level autoregulation of mazEF expression influences cell-to-cell growth rate heterogeneity. The increase in growth rate heterogeneity is governed by the MazE antitoxin, and tuned by the MazF-dependent mazF mRNA cleavage. Also, both autoregulatory features grant rapid exit from the stress caused by mazF overexpression. Time-lapse microscopy revealed that MazF-mediated cleavage of mazF mRNA leads to increased temporal variability in length of individual cells during ectopic mazF overexpression, as explained by a stochastic model indicating that mazEF mRNA cleavage underlies temporal fluctuations in MazF levels during stress.

Division rate, cell size and proteome allocation: impact on gene expression noise and implications for the dynamics of genetic circuits

The cell division rate, size and gene expression programmes change in response to external conditions. These global changes impact on average concentrations of biomolecule and their variability or noise. Gene expression is inherently stochastic, and noise levels of individual proteins depend on synthesis and degradation rates as well as on cell-cycle dynamics. We have modelled stochastic gene expression inside growing and dividing cells to study the effect of division rates on noise in mRNA and protein expression. We use assumptions and parameters relevant to Escherichia coli, for which abundant quantitative data are available. We find that coupling of transcription, but not translation rates to the rate of cell division can result in protein concentration and noise homeostasis across conditions. Interestingly, we find that the increased cell size at fast division rates, observed in E. coli and other unicellular organisms, buffers noise levels even for proteins with decreased expression at faster growth. We then investigate the functional importance of these regulations using gene regulatory networks that exhibit bi-stability and oscillations. We find that network topology affects robustness to changes in division rate in complex and unexpected ways. In particular, a simple model of persistence, based on global physiological feedback, predicts increased proportion of persister cells at slow division rates. Altogether, our study reveals how cell size regulation in response to cell division rate could help controlling gene expression noise. It also highlights that understanding circuits' robustness across growth conditions is key for the effective design of synthetic biological systems.

Division rate, cell size and proteome allocation: impact on gene expression noise and implications for the dynamics of genetic circuits

The cell division rate, size and gene expression programmes change in response to external conditions. These global changes impact on average concentrations of biomolecule and their variability or noise. Gene expression is inherently stochastic, and noise levels of individual proteins depend on synthesis and degradation rates as well as on cell-cycle dynamics. We have modelled stochastic gene expression inside growing and dividing cells to study the effect of division rates on noise in mRNA and protein expression. We use assumptions and parameters relevant to Escherichia coli, for which abundant quantitative data are available. We find that coupling of transcription, but not translation rates to the rate of cell division can result in protein concentration and noise homeostasis across conditions. Interestingly, we find that the increased cell size at fast division rates, observed in E. coli and other unicellular organisms, buffers noise levels even for proteins with decreased expression at faster growth. We then investigate the functional importance of these regulations using gene regulatory networks that exhibit bi-stability and oscillations. We find that network topology affects robustness to changes in division rate in complex and unexpected ways. In particular, a simple model of persistence, based on global physiological feedback, predicts increased proportion of persister cells at slow division rates. Altogether, our study reveals how cell size regulation in response to cell division rate could help controlling gene expression noise. It also highlights that understanding circuits' robustness across growth conditions is key for the effective design of synthetic biological systems.

Toxins, Targets, and Triggers: An Overview of Toxin-Antitoxin Biology

Bacterial toxin-antitoxin (TA) modules are abundant genetic elements that encode a toxin protein capable of inhibiting cell growth and an antitoxin that counteracts the toxin. The majority of toxins are enzymes that interfere with translation or DNA replication, but a wide variety of molecular activities and cellular targets have been described. Antitoxins are proteins or RNAs that often control their cognate toxins through direct interactions and, in conjunction with other signaling elements, through transcriptional and translational regulation of TA module expression. Three major biological functions of TA modules have been discovered, post-segregational killing ("plasmid addiction"), abortive infection (bacteriophage immunity through altruistic suicide), and persister formation (antibiotic tolerance through dormancy). In this review, we summarize the current state of the field and highlight how multiple levels of regulation shape the conditions of toxin activation to achieve the different biological functions of TA modules.

Synchronized switching of multiple toxin-antitoxin modules by (p)ppGpp fluctuation

Toxin–antitoxin (TA) loci are widespread in bacteria including important pathogenic species. Recent studies suggest that TA systems play a key role in persister formation. However, the persistence phenotype shows only weak dependence on the number of TA systems, i.e. they are functionally redundant. We use a mathematical model to investigate the interaction of multiple TA systems in the switching between growth and persistence. We explore two scenarios: (i) TA systems are bistable and each TA system experiences its own noise and (ii) the noise in the level of common stress signal (e.g. (p)ppGpp) coordinates all TA systems simultaneously. We find that in the first scenario the exit from the persister state strongly depends on the number of TA systems. However in the second case, we could reproduce the weak dependence. The duration of the high (p)ppGpp state was found to be the key parameter for persistence. The (p)ppGpp-driven synchronized transition of all TA systems results in the redundancy.

MazF activation promotes translational heterogeneity of the grcA mRNA in Escherichia coli populations

Bacteria adapt to adverse environmental conditions by altering gene expression patterns. Recently, a novel stress adaptation mechanism has been described that allows Escherichia coli to alter gene expression at the post-transcriptional level. The key player in this regulatory pathway is the endoribonuclease MazF, the toxin component of the toxin-antitoxin module mazEF that is triggered by various stressful conditions. In general, MazF degrades the majority of transcripts by cleaving at ACA sites, which results in the retardation of bacterial growth. Furthermore, MazF can process a small subset of mRNAs and render them leaderless by removing their ribosome binding site. MazF concomitantly modifies ribosomes, making them selective for the translation of leaderless mRNAs. In this study, we employed fluorescent reporter-systems to investigate mazEF expression during stressful conditions, and to infer consequences of the mRNA processing mediated by MazF on gene expression at the single-cell level. Our results suggest that mazEF transcription is maintained at low levels in single cells encountering adverse conditions, such as antibiotic stress or amino acid starvation. Moreover, using the grcA mRNA as a model for MazF-mediated mRNA processing, we found that MazF activation promotes heterogeneity in the grcA reporter expression, resulting in a subpopulation of cells with increased levels of GrcA reporter protein.

Toxin-Antitoxin Systems: Implications for Plant Disease

Toxin-antitoxin (TA) systems are gene modules that are ubiquitous in free-living prokaryotes. Diverse in structure, cellular function, and fitness roles, TA systems are defined by the presence of a toxin gene that suppresses bacterial growth and a toxin-neutralizing antitoxin gene, usually encoded in a single operon. Originally viewed as DNA maintenance modules, TA systems are now thought to function in many roles, including bacterial stress tolerance, virulence, phage defense, and biofilm formation. However, very few studies have investigated the significance of TA systems in the context of plant-microbe interactions. This review discusses the potential impact and application of TA systems in plant-associated bacteria, guided by insights gained from animal-pathogenic model systems.

HicA toxin of Escherichia coli derepresses hicAB transcription to selectively produce HicB antitoxin

Antitoxins encoded by type II toxin - antitoxin (TA) modules neutralize cognate toxins by direct protein - protein contact and in addition, regulate TA operon transcription by binding to operators in the promoter regions. On top of the simple negative feed-back regulation, canonical type II TA operons are regulated by a mechanism called 'Conditional Cooperativity'(CC). In CC, the cellular toxin:antitoxin (T:A) ratio controls the transcription-rate such that low T:A ratios favour repression and high T:A ratios favour de-repression of TA operon transcription. Here a new molecular mechanism that secures selective synthesis of antitoxin in the presence of excess toxin was unravelled. The hicAB locus of E. coli K-12 encodes HicA mRNase and HicB antitoxin. It was shown that hicAB is transcribed by two promoters, an upstream one that is activated by CRP-cAMP and competence factor Sxy and a downstream one that is autorepressed solely by HicB. Excess HicA destabilizes the HicB•operator complex in vitro and consistently, activates hicAB transcription in vivo. Remarkably, the hicAB transcript synthesized from the HicB-controlled promoter produces HicB but not HicA. Thus, the HicA-mediated derepression of hicAB transcription provides a mechanism that conditionally and selectively stimulates synthesis of HicB antitoxin under conditions of excess HicA toxin.

Bacterial persistence from a system-level perspective

In recent years, our understanding about bacterial persistence has significantly advanced: we comprehend the persister phenotype better, more triggers for persistence entry have been found, and more insights in the involvement and role of toxin-antitoxin systems and other molecular mechanisms have been unravelled. In this review, we attempt to put these findings into an integrated, system-level perspective. From this point of view, persistence can be seen as a response to a strong perturbation of metabolic homeostasis, either triggered environmentally, or by means of intracellular stochasticity. Metabolic-flux-regulated resource allocation ensures stress protection, and several feedback mechanisms stabilize the cells in this protected state. We hope that this novel view can advance our understanding about persistence.

Molecular mechanisms and clinical implications of bacterial persistence

Any bacterial population harbors a small number of phenotypic variants that survive exposure to high concentrations of antibiotic. Importantly, these so-called 'persister cells' compromise successful antibiotic therapy of bacterial infections and are thought to contribute to the development of antibiotic resistance. Intriguingly, drug-tolerant persisters have also been identified as a factor underlying failure of chemotherapy in tumor cell populations. Recent studies have begun to unravel the complex molecular mechanisms underlying persister formation and revolve around stress responses and toxin-antitoxin modules. Additionally, in vitro evolution experiments are revealing insights into the evolutionary and adaptive aspects of this phenotype. Furthermore, ever-improving experimental techniques are stimulating efforts to investigate persisters in their natural, infection-associated, in vivo environment. This review summarizes recent insights into the molecular mechanisms of persister formation, explains how persisters complicate antibiotic treatment of infections, and outlines emerging strategies to combat these tolerant cells.

Proteolytic crosstalk in multi-protease networks

Processive proteases, such as ClpXP in E. coli, are conserved enzyme assemblies that can recognize and rapidly degrade proteins. These proteases are used for a number of purposes, including degrading mistranslated proteins and controlling cellular stress response. However, proteolytic machinery within the cell is limited in capacity and can lead to a bottleneck in protein degradation, whereby many proteins compete ('queue') for proteolytic resources. Previous work has demonstrated that such queueing can lead to pronounced statistical relationships between different protein counts when proteins compete for a single common protease. However, real cells contain many different proteases, e.g. ClpXP, ClpAP, and Lon in E. coli, and it is not clear how competition between proteins for multiple classes of protease would influence the dynamics of cellular networks. In the present work, we theoretically demonstrate that a multi-protease proteolytic bottleneck can substantially couple the dynamics for both simple and complex (oscillatory) networks, even between substrates with substantially different affinities for protease. For these networks, queueing often leads to strong positive correlations between protein counts, and these correlations are strongest near the queueing theoretic point of balance. Furthermore, we find that the qualitative behavior of these networks depends on the relative size of the absolute affinity of substrate to protease compared to the cross affinity of substrate to protease, leading in certain regimes to priority queue statistics.

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.

Computational Methods to Model Persistence

Bacterial persister cells are dormant cells, tolerant to multiple antibiotics, that are involved in several chronic infections. Toxin-antitoxin modules play a significant role in the generation of such persister cells. Toxin-antitoxin modules are small genetic elements, omnipresent in the genomes of bacteria, which code for an intracellular toxin and its neutralizing antitoxin. In the past decade, mathematical modeling has become an important tool to study the regulation of toxin-antitoxin modules and their relation to the emergence of persister cells. Here, we provide an overview of several numerical methods to simulate toxin-antitoxin modules. We cover both deterministic modeling using ordinary differential equations and stochastic modeling using stochastic differential equations and the Gillespie method. Several characteristics of toxin-antitoxin modules such as protein production and degradation, negative autoregulation through DNA binding, toxin-antitoxin complex formation and conditional cooperativity are gradually integrated in these models. Finally, by including growth rate modulation, we link toxin-antitoxin module expression to the generation of persister cells.

Function does not follow form in gene regulatory circuits

Gene regulatory circuits are to the cell what arithmetic logic units are to the chip: fundamental components of information processing that map an input onto an output. Gene regulatory circuits come in many different forms, distinct structural configurations that determine who regulates whom. Studies that have focused on the gene expression patterns (functions) of circuits with a given structure (form) have examined just a few structures or gene expression patterns. Here, we use a computational model to exhaustively characterize the gene expression patterns of nearly 17 million three-gene circuits in order to systematically explore the relationship between circuit form and function. Three main conclusions emerge. First, function does not follow form. A circuit of any one structure can have between twelve and nearly thirty thousand distinct gene expression patterns. Second, and conversely, form does not follow function. Most gene expression patterns can be realized by more than one circuit structure. And third, multifunctionality severely constrains circuit form. The number of circuit structures able to drive multiple gene expression patterns decreases rapidly with the number of these patterns. These results indicate that it is generally not possible to infer circuit function from circuit form, or vice versa.

Connecting growth with gene expression: of noise and numbers

Growth is a dynamic process whereby cells accumulate mass. Growth rates of single cells are connected to RNA and protein synthesis rates, and therefore with biomolecule numbers. Noise in gene expression depends on these numbers, and is thus linked with cellular growth. Whether these global attributes of the cell participate in gene regulation is still largely unexplored. New experimental and modelling studies suggest that systemic variations in biomolecule numbers can coordinate cellular processes, including growth itself, through global regulatory feedback that acts in addition to genetic regulatory networks. Here, we review these findings and speculate on possible implications of this less appreciated layer of gene regulation for cellular physiology and adaptation to changing environments.

Cut to the chase--Regulating translation through RNA cleavage

Activation of toxin-antitoxin (TA) systems provides an important mechanism for bacteria to adapt to challenging and ever changing environmental conditions. Known TA systems are classified into five families based on the mechanisms of antitoxin inhibition and toxin activity. For type II TA systems, the toxin is inactivated in exponentially growing cells by tightly binding its antitoxin partner protein, which also serves to regulate cellular levels of the complex through transcriptional auto-repression. During cellular stress, however, the antitoxin is degraded thus freeing the toxin, which is then able to regulate central cellular processes, primarily protein translation to adjust cell growth to the new conditions. In this review, we focus on the type II TA pairs that regulate protein translation through cleavage of ribosomal, transfer, or messenger RNA.

Expression, purification, and functional analysis of novel RelE operon from X. nematophila

Bacterial toxin-antitoxin (TA) complexes induce programmed cell death and also function to relieve cell from stress by various response mechanisms. Escherichia coli RelB-RelE TA complex consists of a RelE toxin functionally counteracted by RelB antitoxin. In the present study, a novel homolog of RelE toxin designated as Xn-relE toxin from Xenorhabdus nematophila possessing its own antitoxin designated as Xn-relEAT has been identified. Expression and purification of recombinant proteins under native conditions with GST and Ni-NTA chromatography prove the existence of novel TA module. The expression of recombinant Xn-relE under tightly regulated ara promoter in E. coli Top 10 cells confirms its toxic nature in endogenous toxicity assay. The neutralization activity in endogenous toxicity assay by Xn-relEAT antitoxin confirms its antidote nature when studying the whole TA operon under ara regulated promoter. This study promotes newly discovered TA module to be regarded as important as other proteins of type II toxin-antitoxin system.

The bacterial antitoxin HipB establishes a ternary complex with operator DNA and phosphorylated toxin HipA to regulate bacterial persistence

Nearly all bacteria exhibit a type of phenotypic growth described as persistence that is thought to underlie antibiotic tolerance and recalcitrant chronic infections. The chromosomally encoded high-persistence (Hip) toxin–antitoxin proteins HipA_SO and HipB_SO from Shewanella oneidensis, a proteobacterium with unusual respiratory capacities, constitute a type II toxin–antitoxin protein module. Here we show that phosphorylated HipA_SO can engage in an unexpected ternary complex with HipB_SO and double-stranded operator DNA that is distinct from the prototypical counterpart complex from Escherichia coli. The structure of HipB_SO in complex with operator DNA reveals a flexible C-terminus that is sequestered by HipA_SO in the ternary complex, indicative of its role in binding HipA_SO to abolish its function in persistence. The structure of HipA_SO in complex with a non-hydrolyzable ATP analogue shows that HipA_SO autophosphorylation is coupled to an unusual conformational change of its phosphorylation loop. However, HipA_SO is unable to phosphorylate the translation factor Elongation factor Tu, contrary to previous reports, but in agreement with more recent findings. Our studies suggest that the phosphorylation state of HipA is an important factor in persistence and that the structural and mechanistic diversity of HipAB modules as regulatory factors in bacterial persistence is broader than previously thought.

Latent phenotypes pervade gene regulatory circuits

BACKGROUND

Latent phenotypes are non-adaptive byproducts of adaptive phenotypes. They exist in biological systems as different as promiscuous enzymes and genome-scale metabolic reaction networks, and can give rise to evolutionary adaptations and innovations. We know little about their prevalence in the gene expression phenotypes of regulatory circuits, important sources of evolutionary innovations.

RESULTS

Here, we study a space of more than sixteen million three-gene model regulatory circuits, where each circuit is represented by a genotype, and has one or more functions embodied in one or more gene expression phenotypes. We find that the majority of circuits with single functions have latent expression phenotypes. Moreover, the set of circuits with a given spectrum of functions has a repertoire of latent phenotypes that is much larger than that of any one circuit. Most of this latent repertoire can be easily accessed through a series of small genetic changes that preserve a circuit's main functions. Both circuits and gene expression phenotypes that are robust to genetic change are associated with a greater number of latent phenotypes.

CONCLUSIONS

Our observations suggest that latent phenotypes are pervasive in regulatory circuits, and may thus be an important source of evolutionary adaptations and innovations involving gene regulation.

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.