Stochastic analysis of the SOS response in Escherichia coli

A self-activating electron transfer antibacterial strategy: Co3O4/TiO2 P-N heterojunctions combined with photothermal therapy

Implant-associated infections are significant impediments to successful surgical outcomes, often resulting from persistent bacterial contamination. It has been hypothesized that bacteria can transfer electrons to semiconductors with comparable potential to the biological redox potential (BRP). Building on this concept, we developed an antibiotic-free bactericidal system, Co3O4/TiO2-Ti, capable of achieving real-time and sustainable bactericidal effects. Our study demonstrated that Co3O4/TiO2-Ti, possessing an appropriately set valence band, initiated charge transfer, reactive oxygen species (ROS) production, and membrane damage in adherent Staphylococcus aureus (S. aureus). Notably, in vivo experiments illustrated the remarkable antibacterial activity of Co3O4/TiO2-Ti, while promoting soft-tissue reconstruction and demonstrating excellent cytocompatibility. Transcriptomic analysis further revealed a down-regulation of aerobic respiration-associated genes and an up-regulation of ROS-associated genes in S. aureus in the presence of Co3O4/TiO2-Ti compared to Ti. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis and gene set enrichment analysis (GSEA) identified alterations in respiratory metabolism, oxidative phosphorylation, and the synthesis of amino acid in S. aureus cultured on Co3O4/TiO2-Ti. Furthermore, when combined with near-infrared (NIR) irradiation and photothermal therapy (PTT), Co3O4/TiO2-Ti eliminated 95.71% of floating and adherent S. aureus in vitro. The findings suggest that this antibiotic-free strategy holds substantial promise in enhancing implant sterilization capabilities, thereby contributing to the prevention and treatment of bacterial infections through bandgap engineering of implants and NIR irradiation.

An efficient metal-organic framework-based drug delivery platform for synergistic antibacterial activity and osteogenesis

Bone implants for clinical application should be endowed with antibacterial activity, biocompatibility, and even osteogenesis-promoting properties. In this work, metal-organic framework (MOF) based drug delivery platform was used to modify titanium implants for improved clinical applicability. Methyl Vanillate@Zeolitic Imidazolate Framework-8 (MV@ZIF-8) was immobilized on the polydopamine (PDA) modified titanium. The sustainable release of the Zn²⁺ and MV causes substantial oxidative damage to Escherichia coli (E. coli) and Staphylococcus aureus (S. aureus). The increased reactive oxygen species (ROS) significantly up-regulates the expression of oxidative stress and DNA damage response genes. Meanwhile, the structural disruption of lipid membranes caused by the ROS, the damage caused by Zinc active sites and the damage accelerated by the MV are both involved in inhibiting bacterial proliferation. The up-regulated expression of the osteogenic-related genes and proteins indicated that the MV@ZIF-8 could effectively promote the osteogenic differentiation of the human bone mesenchymal stem cells (hBMSCs). RNA sequencing and Western blotting analysis revealed that the MV@ZIF-8 coating activates the canonical Wnt/β-catenin signaling pathway through the regulation of tumor necrosis factor (TNF) pathway, thereby promoting the osteogenic differentiation of the hBMSCs. This work demonstrates a promising application of the MOF-based drug delivery platform in bone tissue engineering.

Modeling endonuclease colicin-like bacteriocin operons as 'genetic arms' in plasmid-genome conflicts

Plasmids are acellular propagating entities that depend on bacteria, as molecular parasites, for propagation. A 'tussle' between bacteria and plasmid ensues; bacteria for riddance of the plasmid and plasmid for persistence within a live host. Plasmid-maintenance systems such as endonuclease Colicin-Like Bacteriocins (CLBs) ensure plasmid propagation within the population; (i) the plasmid-cured cells are killed by the CLBs; (ii) damaged cells lyse and release the CLBs that eliminate the competitors, and (iii) the released plasmids invade new bacteria. Surprisingly, endonuclease CLB operons occur on bacterial genomes whose significance is unknown. Here, we study genetics, eco-evolutionary drive, and physiological relevance of genomic endonuclease CLB operons. We investigated plasmidic and genomic endonuclease CLB operons using sequence analyses from an eco-evolutionary perspective. We found 1266 genomic and plasmidic endonuclease CLB operons across 30 bacterial genera. Although 51% of the genomes harbor endonuclease CLB operons, the majority of the genomic endonuclease CLB operons lacked a functional lysis gene, suggesting the negative selection of lethal genes. The immunity gene of the endonuclease CLB operon protects the plasmid-cured host, eliminating the metabolic burden. We show mutual exclusivity of endonuclease CLB operons on genomes and plasmids. We propose an anti-addiction hypothesis for genomic endonuclease CLB operons. Using a stochastic hybrid agent-based model, we show that the endonuclease CLB operons on genomes confer an advantage to the host genome in terms of immunity to the toxin and elimination of plasmid burden. The conflict between bacterial genome and plasmids allows the emergence of 'genetic arms' such as CLB operons that regulate the ecological interplay of bacterial genomes and plasmids.

The vulnerable versatility of Salmonella antibiotic persisters during infection

Tolerance and persistence are superficially similar phenomena by which bacteria survive bactericidal antibiotics. It is assumed that the same physiology underlies survival of individual tolerant and persistent bacteria. However, by comparing tolerance and persistence during Salmonella Typhimurium infection, we reveal that these two phenomena are underpinned by different bacterial physiologies. Multidrug-tolerant mutant Salmonella enter a near-dormant state protected from immune-mediated genotoxic damages. However, the numerous tolerant cells, optimized for survival, lack the capabilities necessary to initiate infection relapse following antibiotic withdrawal. In contrast, persisters retain an active state. This leaves them vulnerable to accumulation of macrophage-induced dsDNA breaks but concurrently confers the versatility to initiate infection relapse if protected by RecA-mediated DNA repair. Accordingly, recurrent, invasive, non-typhoidal Salmonella clinical isolates display hallmarks of persistence rather than tolerance during antibiotic treatment. Our study highlights the complex trade-off that antibiotic-recalcitrant Salmonella balance to act as a reservoir for infection relapse.

Covalent grafting of hyperbranched poly-L-lysine on Ti-based implants achieves dual functions of antibacteria and promoted osteointegration in vivo

The dual functional implants of antibacteria and osteointegration are highly demanded in orthopedic and dentistry, especially for patients who suffer from diabetes or osteoporosis simultaneously. However, there is lack of the facile and robust method to produce clinically applicable implants with this dual function although coatings possessing single function have been extensively developed. Herein, hyperbranched poly-L-lysine (HBPL) polymers were covalently immobilized onto the alkali-heat treated titanium (Ti) substrates and implants by using 3-glycidyloxypropyltrimethoxysilane (GPTMS) as the coupling agent, which displayed excellent antibacterial activity against S. aureus and E. coli with an efficiency as high as 89.4% and 92.2% in vitro, respectively. The HBPL coating also significantly promoted the adhesion, spreading, proliferation and osteogenic differentiation of MC3T3-E1 cells in vitro. Furthermore, the results of a S. aureus infection rat model in vivo ulteriorly verified that the HBPL-modified screws had good antibacterial and anti-inflammatory abilities at an early stage of implantation and better osteointegration compared with the control Ti screws.

WD40 domain of RqkA regulates its kinase activity and role in extraordinary radioresistance of D. radiodurans

RqkA, a DNA damage responsive serine/threonine kinase, is characterized for its role in DNA repair and cell division in D. radiodurans. It has a unique combination of a kinase domain at N-terminus and a WD40 type domain at C-terminus joined through a linker. WD40 domain is comprised of eight β-propeller repeats held together via 'tryptophan-docking motifs' and forming a typical 'velcro' closure structure. RqkA mutants lacking the WD40 region (hereafter referred to as WD mutant) could not complement RqkA loss in γ radiation resistance in D. radiodurans and lacked γ radiation-mediated activation of kinase activity in vivo. WD mutants failed to phosphorylate its cognate substrate (e.g. DrRecA) in surrogate E. coli cells. Unlike wild-type enzyme, the kinase activity of its WD40 mutants was not stimulated by pyrroloquinoline quinine (PQQ) indicating the role of the WD motifs in PQQ interaction and stimulation of its kinase activity. Together, results highlighted the importance of the WD40 domain in the regulation of RqkA kinase signaling functions in vivo, and thus, the role of WD40 domain in the regulation of any STPK is first time demonstrated in bacteria.Communicated by Ramaswamy H. Sarma.

The PafBC-mediated response sensitizes a bistable DNA damage response in Mycobacteria

Due to the genotoxically challenging environments in which they live in, Mycobacteria have a complex DNA damage repair system that is governed by two major DNA damage responses, namely, the LexA/RecA-dependent response and the newly characterized PafBC-mediated response (Müller et al., 2018). The LexA/RecA-dependent response is a well-known bistable response found in different types of bacteria, and the Mycobacteria-specific PafBC-mediated response interacts with and modifies the LexA/RecA-dependent response (Müller et al., 2018). The interaction between the LexA/RecA-dependent response and the PafBC-mediated response has not been characterized mathematically. Our analysis shows that the addition of the PafBC-mediated response sensitizes the overall DNA damage response, effectively lowering the DNA damage rate threshold for activation.

A Roadblock-and-Kill Mechanism of Action Model for the DNA-Targeting Antibiotic Ciprofloxacin

Fluoroquinolones, antibiotics that cause DNA damage by inhibiting DNA topoisomerases, are clinically important, but their mechanism of action is not yet fully understood. In particular, the dynamical response of bacterial cells to fluoroquinolone exposure has hardly been investigated, although the SOS response, triggered by DNA damage, is often thought to play a key role. Here, we investigated the growth inhibition of the bacterium Escherichia coli by the fluoroquinolone ciprofloxacin at low concentrations. We measured the long-term and short-term dynamical response of the growth rate and DNA production rate to ciprofloxacin at both the population and single-cell levels. We show that, despite the molecular complexity of DNA metabolism, a simple roadblock-and-kill model focusing on replication fork blockage and DNA damage by ciprofloxacin-poisoned DNA topoisomerase II (gyrase) quantitatively reproduces long-term growth rates in the presence of ciprofloxacin. The model also predicts dynamical changes in the DNA production rate in wild-type E. coli and in a recombination-deficient mutant following a step-up of ciprofloxacin. Our work highlights that bacterial cells show a delayed growth rate response following fluoroquinolone exposure. Most importantly, our model explains why the response is delayed: it takes many doubling times to fragment the DNA sufficiently to inhibit gene expression. We also show that the dynamical response is controlled by the timescale of DNA replication and gyrase binding/unbinding to the DNA rather than by the SOS response, challenging the accepted view. Our work highlights the importance of including detailed biophysical processes in biochemical-systems models to quantitatively predict the bacterial response to antibiotics.

Dynamics of ColicinE2 production and release determine the competitive success of a toxin-producing bacterial population

The release of toxins is one mechanism used by bacterial species to establish dominance over competitors, but how the dynamics of toxin expression determine the competitive success of a toxin-producing population is largely unknown. Here, we investigate how the expression dynamics of ColicinE2 - a toxic bacteriocin - affect competition between toxin-producing and toxin-sensitive strains of Escherichia coli. We demonstrate that, in addition to genetic modifications in the toxin expression system, alterations of the growth medium can be used to modulate the timing of toxin production and the amount of toxin released. Thus cells that release the toxin at later times can accumulate more colicin. In experiments, we found that delaying toxin release does not significantly alter competition outcome. However, our theoretical analysis allowed us to assess the relative contributions of release time and toxin level to the competitive success of the producer strain, that might counteract each other in experiments. The results reveal that the importance of delaying toxin release lies in increasing the toxin amount. This is a more effective strategy for the toxin-producing strain than prompt discharge of the colicin. In summary, our study shows how the toxin release dynamics influence the competitive success of the toxin-producing bacterial population.

Gene expression noise in a complex artificial toxin expression system

Gene expression is an intrinsically stochastic process. Fluctuations in transcription and translation lead to cell-to-cell variations in mRNA and protein levels affecting cellular function and cell fate. Here, using fluorescence time-lapse microscopy, we quantify noise dynamics in an artificial operon in Escherichia coli, which is based on the native operon of ColicinE2, a toxin. In the natural system, toxin expression is controlled by a complex regulatory network; upon induction of the bacterial SOS response, ColicinE2 is produced (cea gene) and released (cel gene) by cell lysis. Using this ColicinE2-based operon, we demonstrate that upon induction of the SOS response noise of cells expressing the operon is significantly lower for the (mainly) transcriptionally regulated gene cea compared to the additionally post-transcriptionally regulated gene cel. Likewise, we find that mutations affecting the transcriptional regulation by the repressor LexA do not significantly alter the population noise, whereas specific mutations to post-transcriptionally regulating units, strongly influence noise levels of both genes. Furthermore, our data indicate that global factors, such as the plasmid copy number of the operon encoding plasmid, affect gene expression noise of the entire operon. Taken together, our results provide insights on how noise in a native toxin-producing operon is controlled and underline the importance of post-transcriptional regulation for noise control in this system.

A Quantitative Model Explains Single-Cell Dynamics of the Adaptive Response in Escherichia coli

DNA damage caused by alkylating chemicals induces an adaptive response in Escherichia coli that increases the tolerance of cells to further damage. Signaling of the response occurs through irreversible methylation of the Ada protein, which acts as a DNA repair protein and damage sensor. Methylated Ada induces its own gene expression through a positive feedback loop. However, random fluctuations in the abundance of Ada jeopardize the reliability of the induction signal. I developed a quantitative model to test how gene expression noise and feedback amplification affect the fidelity of the adaptive response. A remarkably simple model accurately reproduced experimental observations from single-cell measurements of gene expression dynamics in a microfluidic device. Stochastic simulations showed that delays in the adaptive response are a direct consequence of the very low number of Ada molecules present to signal DNA damage. For cells that have zero copies of Ada, response activation becomes a memoryless process that is dictated by an exponential waiting time distribution between basal Ada expression events. Experiments also confirmed the model prediction that the strength of the adaptive response drops with an increasing growth rate of cells.

Locality of interactions in three-strain bacterial competition in E. coli

The population dynamics that determine the composition and stability of ecosystems ultimately emerge from interactions between individual organisms. One well-studied system is the three-strain E. coli interaction of a heterogeneously toxin-producing C strain that interacts with a toxin-sensitive S and a toxin-resistant R strain. Here, we employ a multi-scale fluorescence microscopy approach, that has been proven useful in identifying previously unknown or underestimated stochastic effects in C-S competition. This approach allows us to investigate the microscopic interaction of the R strain and to quantify the role of stochastic effects in the spatially structured C-R-S interaction. We show that the early colony patterning at 12 h and at small length scales (near single cell level) is characterized by a number of microscopic variables (the number of C and R cell clusters and the area occupied by S) and is subject to random processes in positioning and toxin production. Then, in a second competition phase, mainly deterministic processes such as bacterial growth and global toxin action determine the following population dynamics. Consequently, together with environmental factors, the microscopic variables were predictive of the competition outcome. However, interactions of neighboring R and C clusters could amplify local variations. If R clusters originated near a C cell cluster, R could profit from the toxin produced by C without bearing the cost of production-a mechanism called cheating. By combining information from the micro- and macro-scale dynamics, we can estimate the distance at which the cheating interaction significantly changes to be in the order of 250 µm. In summary, after an initial phase influenced by stochastic patterning, largely deterministic growth dynamics follow, which are additionally affected by local interactions of neighboring clusters. As such, the results underline the importance of stochasticity and local effects in the context of ecological interactions.

Non-equilibrium repressor binding kinetics link DNA damage dose to transcriptional timing within the SOS gene network

Biochemical pathways are often genetically encoded as simple transcription regulation networks, where one transcription factor regulates the expression of multiple genes in a pathway. The relative timing of each promoter’s activation and shut-off within the network can impact physiology. In the DNA damage repair pathway (known as the SOS response) of Escherichia coli, approximately 40 genes are regulated by the LexA repressor. After a DNA damaging event, LexA degradation triggers SOS gene transcription, which is temporally separated into subsets of ‘early’, ‘middle’, and ‘late’ genes. Although this feature plays an important role in regulating the SOS response, both the range of this separation and its underlying mechanism are not experimentally defined. Here we show that, at low doses of DNA damage, the timing of promoter activities is not separated. Instead, timing differences only emerge at higher levels of DNA damage and increase as a function of DNA damage dose. To understand mechanism, we derived a series of synthetic SOS gene promoters which vary in LexA-operator binding kinetics, but are otherwise identical, and then studied their activity over a large dose-range of DNA damage. In distinction to established models based on rapid equilibrium assumptions, the data best fit a kinetic model of repressor occupancy at promoters, where the drop in cellular LexA levels associated with higher doses of DNA damage leads to non-equilibrium binding kinetics of LexA at operators. Operators with slow LexA binding kinetics achieve their minimal occupancy state at later times than operators with fast binding kinetics, resulting in a time separation of peak promoter activity between genes. These data provide insight into this remarkable feature of the SOS pathway by demonstrating how a single transcription factor can be employed to control the relative timing of each gene’s transcription as a function of stimulus dose.

As the precise timing of gene expression is critical for cells to respond and adapt to new environments, it is important to understand the underlying mechanisms which control this timing. In this report, we studied the timing of transcription for genes in the bacterial DNA damage repair pathway (known as the SOS response), a regulatory system where each gene is controlled by the same transcriptional repressor, LexA. By specifically isolating the role of the LexA binding interaction at SOS gene promoters, we found a relationship between the amount of DNA damage incurred by the cell, LexA binding kinetics at a promoter, and the timing of promoter activation. Our data fit a kinetic model that reveals how a disequilibrium between the LexA-operator binding reaction and cellular LexA concentrations causes timing differences between genes to emerge only at higher doses of DNA damage. Taken together, we show that non-equilibrium DNA binding kinetics is the mechanism by which a single transcription factor can modulate timing differences across an entire network of genes as a function of stimulus dose.

Effects of stochasticity and division of labor in toxin production on two-strain bacterial competition in Escherichia coli

In phenotypically heterogeneous microbial populations, the decision to adopt one or another phenotype is often stochastically regulated. However, how this stochasticity affects interactions between competing microbes in mixed communities is difficult to assess. One example of such an interaction system is the competition of an Escherichia coli strain C, which performs division of labor between reproducers and self-sacrificing toxin producers, with a toxin-sensitive strain S. The decision between reproduction or toxin production within a single C cell is inherently stochastic. Here, combining experimental and theoretical approaches, we demonstrate that this stochasticity in the initial phase of colony formation is the crucial determinant for the competition outcome. In the initial phase (t < 12h), stochasticity influences the formation of viable C clusters at the colony edge. In the subsequent phase, the effective fitness differences (set primarily by the degree of division of labor in the C strain population) dictate the deterministic population dynamics and consequently competition outcome. In particular, we observe that competitive success of the C strain is only found if (i) a C edge cluster has formed at the end of the initial competition phase and (ii) the beneficial and detrimental effects of toxin production are balanced, which is the case at intermediate toxin producer fractions. Our findings highlight the importance of stochastic processes during the initial phase of colony formation, which might be highly relevant for other microbial community interactions in which the random choice between phenotypes can have long-lasting consequences for community fate.

Competition is the dominant interaction type between species of bacteria. Bacterial toxin-mediated competition is often accompanied by a division of labor between toxin-producing cells and reproducers within a species. In populations with large cell numbers, the stochastic dynamics that determine if an individual cell switches into the toxin-producing state are often unnoticed. Consequently, we know little about how stochastic effects influence bacterial competition. Here, combining experimental and theoretical efforts, we study the competition of a toxin-producing strain with a toxin-sensitive strain. By correlating the initial conditions—at near single-cell level—to the macroscopic competition outcome, we investigate both the importance of the division of labor as well as the influence of the stochastic toxin production dynamics on competition outcome. Our results highlight the impact of the initial phase of competition as a major determinant for the success of the toxin-producing population.

Hierarchical Post-transcriptional Regulation of Colicin E2 Expression in Escherichia coli

Post-transcriptional regulation of gene expression plays a crucial role in many bacterial pathways. In particular, the translation of mRNA can be regulated by trans-acting, small, non-coding RNAs (sRNAs) or mRNA-binding proteins, each of which has been successfully treated theoretically using two-component models. An important system that includes a combination of these modes of post-transcriptional regulation is the Colicin E2 system. DNA damage, by triggering the SOS response, leads to the heterogeneous expression of the Colicin E2 operon including the cea gene encoding the toxin colicin E2, and the cel gene that codes for the induction of cell lysis and release of colicin. Although previous studies have uncovered the system's basic regulatory interactions, its dynamical behavior is still unknown. Here, we develop a simple, yet comprehensive, mathematical model of the colicin E2 regulatory network, and study its dynamics. Its post-transcriptional regulation can be reduced to three hierarchically ordered components: the mRNA including the cel gene, the mRNA-binding protein CsrA, and an effective sRNA that regulates CsrA. We demonstrate that the stationary state of this system exhibits a pronounced threshold in the abundance of free mRNA. As post-transcriptional regulation is known to be noisy, we performed a detailed stochastic analysis, and found fluctuations to be largest at production rates close to the threshold. The magnitude of fluctuations can be tuned by the rate of production of the sRNA. To study the dynamics in response to an SOS signal, we incorporated the LexA-RecA SOS response network into our model. We found that CsrA regulation filtered out short-lived activation peaks and caused a delay in lysis gene expression for prolonged SOS signals, which is also seen in experiments. Moreover, we showed that a stochastic SOS signal creates a broad lysis time distribution. Our model thus theoretically describes Colicin E2 expression dynamics in detail and reveals the importance of the specific regulatory components for the timing of toxin release.

Intracellular d-Serine Accumulation Promotes Genetic Diversity via Modulated Induction of RecA in Enterohemorrhagic Escherichia coli

We recently discovered that exposure of enterohemorrhagic Escherichia coli (EHEC) to d-serine resulted in accumulation of this unusual amino acid, induction of the SOS regulon, and downregulation of the type III secretion system that is essential for efficient colonization of the host. Here, we have investigated the physiological relevance of this elevated SOS response, which is of particular interest given the presence of Stx toxin-carrying lysogenic prophages on the EHEC chromosome that are activated during the SOS response. We found that RecA elevation in response to d-serine, while being significant, was heterogeneous and not capable of activating stx expression or stx phage transduction to a nonlysogenic recipient. This "SOS-like response" was, however, capable of increasing the mutation frequency associated with low-level RecA activity, thus promoting genetic diversity. Furthermore, this response was entirely dependent on RecA and enhanced in the presence of a DNA-damaging agent, indicating a functional SOS response, but did not result in observable cleavage of the LexA repressor alone, indicating a controlled mechanism of induction. This work demonstrates that environmental factors not usually associated with DNA damage are capable of promoting an SOS-like response. We propose that this modulated induction of RecA allows EHEC to adapt to environmental insults such as d-serine while avoiding unwanted phage-induced lysis.

IMPORTANCE

The SOS response is a global stress network that is triggered by the presence of DNA damage due to breakage or stalled replication forks. Activation of the SOS response can trigger the replication of lytic bacteriophages and promote genetic diversification through error-prone polymerases. We have demonstrated that the host-associated metabolite d-serine contributes to Escherichia coli niche specification and accumulates inside cells that cannot catabolize it. This results in a modulated activation of the SOS antirepressor RecA that is insufficient to trigger lytic bacteriophage but capable of increasing the SOS-associated mutation frequency. These findings describe how relevant signals not normally associated with DNA damage can hijack the SOS response, promoting diversity as E. coli strains adapt while avoiding unwanted phage lysis.

Studying the organization of DNA repair by single-cell and single-molecule imaging

•

Single-cell experiments to study stochastic events and heterogeneity in DNA repair.

•

Quantifying DNA repair protein concentration, diffusion, and localization in cells.

•

Direct observation of DNA repair using photoactivated single-molecule tracking.

DNA repair safeguards the genome against a diversity of DNA damaging agents. Although the mechanisms of many repair proteins have been examined separately in vitro, far less is known about the coordinated function of the whole repair machinery in vivo. Furthermore, single-cell studies indicate that DNA damage responses generate substantial variation in repair activities across cells. This review focuses on fluorescence imaging methods that offer a quantitative description of DNA repair in single cells by measuring protein concentrations, diffusion characteristics, localizations, interactions, and enzymatic rates. Emerging single-molecule and super-resolution microscopy methods now permit direct visualization of individual proteins and DNA repair events in vivo. We expect much can be learned about the organization of DNA repair by linking cell heterogeneity to mechanistic observations at the molecular level.

DNA damage responses in prokaryotes: regulating gene expression, modulating growth patterns, and manipulating replication forks

Recent advances in the area of bacterial DNA damage responses are reviewed here. The SOS pathway is still the major paradigm of bacterial DNA damage response, and recent studies have clarified the mechanisms of SOS induction and key physiological roles of SOS including a very major role in genetic exchange and variation. When considering diverse bacteria, it is clear that SOS is not a uniform pathway with one purpose, but rather a platform that has evolved for differing functions in different bacteria. Relating in part to the SOS response, the field has uncovered multiple apparent cell-cycle checkpoints that assist cell survival after DNA damage and remarkable pathways that induce programmed cell death in bacteria. Bacterial DNA damage responses are also much broader than SOS, and several important examples of LexA-independent regulation will be reviewed. Finally, some recent advances that relate to the replication and repair of damaged DNA will be summarized.

Structure-function study of deinococcal serine/threonine protein kinase implicates its kinase activity and DNA repair protein phosphorylation roles in radioresistance of Deinococcus radiodurans

The DR2518 (RqkA) a eukaryotic type serine/threonine protein kinase in Deinococcus radiodurans was characterized for its role in bacterial response to oxidative stress and DNA damage. The K42A, S162A, T169A and S171A mutation in RqkA differentially affected its kinase activity and functional complementation for γ radiation resistance in Δdr2518 mutant. For example, K42A mutant was completely inactive and showed no complementation while S171A, T169A and T169A/S171A mutants were less active and complemented proportionally to different levels as compared to wild type. Amongst, different DNA binding proteins that purified RqkA could phosphorylate, PprA a DNA repair protein, phosphorylation had improved its affinity to DNA by 4 fold and could enhance its supportive role in intermolecular ligation by T4 DNA ligase. RqkA phosphorylates PprA at threonine 72 (T72), serine 112 (S112) and threonine 144 (T144) in vitro with the majority of it goes to T72 site. Unlike wild type PprA and single mutants of T72, S112 and T144 residues, the T72AS112A double and T72AS112AT144A triple mutant derivatives of PprA did not phosphorylate in vivo and also failed to complement PprA loss in D. radiodurans. Deletion of rqkA in pprA::cat background enhanced radiosensitivity of pprA mutant, which became nearly similar to ΔrqkA resistance to γ radiation. These results suggested that K42 of RqkA is essential for catalytic functions and the kinase activity of RqkA as well as phosphorylation of PprA have roles in γ radiation resistance of D. radiodurans.

Evolutionary tradeoffs between economy and effectiveness in biological homeostasis systems

Biological regulatory systems face a fundamental tradeoff: they must be effective but at the same time also economical. For example, regulatory systems that are designed to repair damage must be effective in reducing damage, but economical in not making too many repair proteins because making excessive proteins carries a fitness cost to the cell, called protein burden. In order to see how biological systems compromise between the two tasks of effectiveness and economy, we applied an approach from economics and engineering called Pareto optimality. This approach allows calculating the best-compromise systems that optimally combine the two tasks. We used a simple and general model for regulation, known as integral feedback, and showed that best-compromise systems have particular combinations of biochemical parameters that control the response rate and basal level. We find that the optimal systems fall on a curve in parameter space. Due to this feature, even if one is able to measure only a small fraction of the system's parameters, one can infer the rest. We applied this approach to estimate parameters in three biological systems: response to heat shock and response to DNA damage in bacteria, and calcium homeostasis in mammals.

Many systems in the cell work to keep homeostasis, or balance. For example, damage repair systems make special repair proteins to resolve damage. These systems typically have many biochemical parameters such as biochemical rate constants, and it is not clear how much of the huge parameter space is filled by actual biological systems. We examined how natural selection acts on these systems when there are two important tasks: effectiveness – rapidly repairing damage, and economy – avoiding excessive production of repair proteins. We find that this multi-task optimization situation leads to natural selection of circuits that lie on a curve in parameter space. Thus, most of parameter space is empty. Estimating only a few parameters of the circuit is enough to predict the rest. This approach allowed us to estimate parameters for bacterial heat shock and DNA repair systems, and for a mammalian hormone system responsible for calcium homeostasis.

Current knowledge on regulatory RNAs and their machineries in Staphylococcus aureus

Staphylococcus aureus is one of the major human pathogens, which causes numerous community-associated and hospital-acquired infections. The regulation of the expression of numerous virulence factors is coordinated by complex interplays between two component systems, transcriptional regulatory proteins, and regulatory RNAs. Recent studies have identified numerous novel RNAs comprising cis-acting regulatory RNAs, antisense RNAs, small non coding RNAs and small mRNAs encoding peptides. We present here several examples of RNAs regulating S. aureus pathogenicity and describe various aspects of antisense regulation.

A comparative proteomic analysis of the early response to compatible symbiotic bacteria in the roots of a supernodulating soybean variety

To reveal the processes involved in the early stages of symbiosis between soybean plants and root nodule bacteria, we conducted a proteomic analysis of the response to bacterial inoculation in the roots of supernodulating (En-b0-1) and non-nodulating (En1282) varieties, and their parental normal-nodulating variety (Enrei). A total of 56 proteins were identified from 48 differentially expressed protein spots in normal-nodulating variety after bacterial inoculation. Among 56 proteins, metabolism- and energy production-related proteins were upregulated in supernodulating and downregulated in non-nodulating varieties compared to normal-nodulating variety. The supernodulating and non-nodulating varieties responded oppositely to bacterial inoculation with respect to the expression of 11 proteins. Seven proteins of these proteins was downregulated in supernodulating varieties compared to non-nodulating variety, but expression of proteasome subunit alpha type 6, gamma glutamyl hydrolase, glucan endo-1,3-beta glucosidase, and nodulin 35 was upregulated. The expression of seven proteins mirrored the degree of nodule formation. At the transcript level, expression of stem 31kDa glycoprotein, leucine aminopeptidase, phosphoglucomutase, and peroxidase was downregulated in the supernodulating variety compared to the non-nodulating variety, and their expression in the normal-nodulating variety was intermediate. These results suggest that suppression of the autoregulatory mechanism in the supernodulating variety might be due to negative regulation of defense and signal transduction-related processes.

Multi-scale stochastic simulation of diffusion-coupled agents and its application to cell culture simulation

Many biological systems consist of multiple cells that interact by secretion and binding of diffusing molecules, thus coordinating responses across cells. Techniques for simulating systems coupling extracellular and intracellular processes are very limited. Here we present an efficient method to stochastically simulate diffusion processes, which at the same time allows synchronization between internal and external cellular conditions through a modification of Gillespie's chemical reaction algorithm. Individual cells are simulated as independent agents, and each cell accurately reacts to changes in its local environment affected by diffusing molecules. Such a simulation provides time-scale separation between the intra-cellular and extra-cellular processes. We use our methodology to study how human monocyte-derived dendritic cells alert neighboring cells about viral infection using diffusing interferon molecules. A subpopulation of the infected cells reacts early to the infection and secretes interferon into the extra-cellular medium, which helps activate other cells. Findings predicted by our simulation and confirmed by experimental results suggest that the early activation is largely independent of the fraction of infected cells and is thus both sensitive and robust. The concordance with the experimental results supports the value of our method for overcoming the challenges of accurately simulating multiscale biological signaling systems.

Optimality in DNA repair

DNA within cells is subject to damage from various sources. Organisms have evolved a number of mechanisms to repair DNA damage. The activity of repair enzymes carries its own risk, however, because the repair of two nearby lesions may lead to the breakup of DNA and result in cell death. We propose a mathematical theory of the damage and repair process in the important scenario where lesions are caused in bursts. We use this model to show that there is an optimum level of repair enzymes within cells which optimises the cell's response to damage. This optimal level is explained as the best trade-off between fast repair and a low probability of causing double-stranded breaks. We derive our results analytically and test them using stochastic simulations, and compare our predictions with current biological knowledge.

Characterization of the role of the RadS/RadR two-component system in the radiation resistance of Deinococcus radiodurans

Deinococcus radiodurans shows extraordinary tolerance to DNA damage, and exhibits differential gene expression and protein recycling. A putative response regulator, the DRB0091 (RadR) ORF, was identified from a pool of DNA-binding proteins induced in response to gamma radiation in this bacterium. radR is located upstream of drB0090, which encodes a putative sensor histidine kinase (RadS) on the megaplasmid. Deletion of these genes both individually and together resulted in hypersensitivity to DNA-damaging agents and a delayed or altered double-strand break repair. A ΔradRradS double mutant and a ΔradR single mutant showed nearly identical responses to gamma radiation and UVC. Wild-type RadR and RadS complemented the corresponding mutant strains, but also exhibited significant cross-complementation, albeit at lower doses of gamma radiation. The radS transcript was not detected in the ΔradR mutant, suggesting the existence of a radRS operon. Recombinant RadS was autophosphorylated and could catalyse the transfer of γ phosphate from ATP to RadR in vitro. These results indicated the functional interaction of RadS and RadR, and suggested a role for the RadS/RadR two-component system in the radiation resistance of this bacterium.

Diversity and distribution of transcription factors: their partner domains play an important role in regulatory plasticity in bacteria

The ability of bacteria to deal with diverse environmental changes depends on their repertoire of genes and their ability to regulate their expression. In this process, DNA-binding transcription factors (TFs) have a fundamental role because they affect gene expression positively and/or negatively depending on operator context and ligand-binding status. Here, we show an exhaustive analysis of winged helix-turn-helix domains (wHTHs), a class of DNA-binding TFs. These proteins were identified in high proportions and widely distributed in bacteria, representing around half of the total TFs identified so far. In addition, we evaluated the repertoire of wHTHs in terms of their partner domains (PaDos), identifying a similar trend, as with TFs, i.e. they are abundant and widely distributed in bacteria. Based on the PaDos, we defined three main groups of families: (i) monolithic, those families with little PaDo diversity, such as LysR; (ii) promiscuous, those families with a high PaDo diversity; and (iii) monodomain, with families of small sizes, such as MarR. These findings suggest that PaDos have a very important role in the diversification of regulatory responses in bacteria, probably contributing to their regulatory complexity. Thus, the TFs discriminate over longer regions on the DNA through their diverse DNA-binding domains. On the other hand, the PaDos would allow a great flexibility for transcriptional regulation due to their ability to sense diverse stimuli through a variety of ligand-binding compounds.

Role of cell-to-cell variability in activating a positive feedback antiviral response in human dendritic cells

In the first few hours following Newcastle disease viral infection of human monocyte-derived dendritic cells, the induction of IFNB1 is extremely low and the secreted type I interferon response is below the limits of ELISA assay. However, many interferon-induced genes are activated at this time, for example DDX58 (RIGI), which in response to viral RNA induces IFNB1. We investigated whether the early induction of IFNBI in only a small percentage of infected cells leads to low level IFN secretion that then induces IFN-responsive genes in all cells. We developed an agent-based mathematical model to explore the IFNBI and DDX58 temporal dynamics. Simulations showed that a small number of early responder cells provide a mechanism for efficient and controlled activation of the DDX58-IFNBI positive feedback loop. The model predicted distributions of single cell responses that were confirmed by single cell mRNA measurements. The results suggest that large cell-to-cell variation plays an important role in the early innate immune response, and that the variability is essential for the efficient activation of the IFNB1 based feedback loop.

Noise-induced modulation of the relaxation kinetics around a non-equilibrium steady state of non-linear chemical reaction networks

Stochastic effects from correlated noise non-trivially modulate the kinetics of non-linear chemical reaction networks. This is especially important in systems where reactions are confined to small volumes and reactants are delivered in bursts. We characterise how the two noise sources confinement and burst modulate the relaxation kinetics of a non-linear reaction network around a non-equilibrium steady state. We find that the lifetimes of species change with burst input and confinement. Confinement increases the lifetimes of all species that are involved in any non-linear reaction as a reactant. Burst monotonically increases or decreases lifetimes. Competition between burst-induced and confinement-induced modulation may hence lead to a non-monotonic modulation. We quantify lifetime as the integral of the time autocorrelation function (ACF) of concentration fluctuations around a non-equilibrium steady state of the reaction network. Furthermore, we look at the first and second derivatives of the ACF, each of which is affected in opposite ways by burst and confinement. This allows discriminating between these two noise sources. We analytically derive the ACF from the linear Fokker–Planck approximation of the chemical master equation in order to establish a baseline for the burst-induced modulation at low confinement. Effects of higher confinement are then studied using a partial-propensity stochastic simulation algorithm. The results presented here may help understand the mechanisms that deviate stochastic kinetics from its deterministic counterpart. In addition, they may be instrumental when using fluorescence-lifetime imaging microscopy (FLIM) or fluorescence-correlation spectroscopy (FCS) to measure confinement and burst in systems with known reaction rates, or, alternatively, to correct for the effects of confinement and burst when experimentally measuring reaction rates.

Growth parameter components of adaptive specificity during experimental evolution of the UVR-inducible mutator Pseudomonas cichorii 302959

BACKGROUND

Mutagenic DNA repair (MDR) transiently increases mutation rate through the activation of low-fidelity repair polymerases in response to specific, DNA-damaging environmental stress conditions such as ultraviolet radiation (UVR) exposure. These repair polymerases also confer UVR tolerance, intimately linking mutability and survival in bacteria that colone habitats subject to regular UVR exposure.

METHODOLOGY/PRINCIPAL FINDINGS

Here, we investigate adaptive specificity in experimental lineages of the highly UVR-mutable epiphytic plant pathogen Pseudomonas cichorii 302959. Relative fitness measurements of isolates and population samples from replicate lineages indicated that adaptive improvements emerged early in all lineages of our evolution experiment and specific increases in relative fitness correlated with distinct improvements in doubling and lag times. Adaptive improvements gained under UVR and non-UVR conditions were acquired preferentially, and differentially contributed to relative fitness under varied growth conditions.

CONCLUSIONS

These results support our earlier observations that MDR activation may contribute to gains in relative fitness without impeding normal patterns of adaptive specificity in P. cichorii 302959.

Characterization of a DNA damage-inducible membrane protein kinase from Deinococcus radiodurans and its role in bacterial radioresistance and DNA strand break repair

Deinococcus radiodurans mutant lacking pyrroloquinoline-quinone (PQQ) synthesis shows sensitivity to γ-rays and impairment of DNA double strand break repair. The genome of this bacterium encodes five putative proteins having multiple PQQ binding motifs. The deletion mutants of corresponding genes were generated, and their response to DNA damage was monitored. Only the Δdr2518 mutant exhibited higher sensitivity to DNA damage. Survival of these cells decreased by 3-log cycle both at 6 kGy γ-rays and 1200 Jm(-2) UV (254 nm) radiation, and 2.5-log cycle upon 14 days desiccation at 5% humidity. The Δdr2518 mutant showed complete inhibition of DSB repair until 24 h PIR and disappearance of a few phosphoproteins. The Δdr2518pqqE:cat double mutant showed γ-ray sensitivity similar to Δdr2518 indicating functional interaction of these genes in D. radiodurans. DR2518 contains a eukaryotic type Ser/Thr kinase domain and structural topology suggesting stress responsive transmembrane protein. Its autokinase activity in solution was stimulated by nearly threefold with PQQ and twofold with linear DNA, but not with circular plasmid DNA. More than 15-fold increase in dr2518 transcription and several-fold enhanced in vivo phosphorylation of DR2518 were observed in response to γ irradiation. These results suggest that DR2518 as a DNA damage-responsive protein kinase plays an important role in radiation resistance and DNA strand break repair in D. radiodurans.

Plato's cave algorithm: inferring functional signaling networks from early gene expression shadows

Improving the ability to reverse engineer biochemical networks is a major goal of systems biology. Lesions in signaling networks lead to alterations in gene expression, which in principle should allow network reconstruction. However, the information about the activity levels of signaling proteins conveyed in overall gene expression is limited by the complexity of gene expression dynamics and of regulatory network topology. Two observations provide the basis for overcoming this limitation: a. genes induced without de-novo protein synthesis (early genes) show a linear accumulation of product in the first hour after the change in the cell's state; b. The signaling components in the network largely function in the linear range of their stimulus-response curves. Therefore, unlike most genes or most time points, expression profiles of early genes at an early time point provide direct biochemical assays that represent the activity levels of upstream signaling components. Such expression data provide the basis for an efficient algorithm (Plato's Cave algorithm; PLACA) to reverse engineer functional signaling networks. Unlike conventional reverse engineering algorithms that use steady state values, PLACA uses stimulated early gene expression measurements associated with systematic perturbations of signaling components, without measuring the signaling components themselves. Besides the reverse engineered network, PLACA also identifies the genes detecting the functional interaction, thereby facilitating validation of the predicted functional network. Using simulated datasets, the algorithm is shown to be robust to experimental noise. Using experimental data obtained from gonadotropes, PLACA reverse engineered the interaction network of six perturbed signaling components. The network recapitulated many known interactions and identified novel functional interactions that were validated by further experiment. PLACA uses the results of experiments that are feasible for any signaling network to predict the functional topology of the network and to identify novel relationships.