Isolation and characterization of LexA mutant repressors with enhanced DNA binding affinity

Single-molecule imaging of LexA degradation in Escherichia coli elucidates regulatory mechanisms and heterogeneity of the SOS response

The bacterial SOS response stands as a paradigm of gene networks controlled by a master transcriptional regulator. Self-cleavage of the SOS repressor, LexA, induces a wide range of cell functions that are critical for survival and adaptation when bacteria experience stress conditions¹, including DNA repair², mutagenesis^3,4, horizontal gene transfer^5–7, filamentous growth, and the induction of bacterial toxins^8–12, toxin-antitoxin systems¹³, virulence factors^6,14, and prophages^15–17. SOS induction is also implicated in biofilm formation and antibiotic persistence^11,18–20. Considering the fitness burden of these functions, it is surprising that the expression of LexA-regulated genes is highly variable across cells^10,21–23 and that cell subpopulations induce the SOS response spontaneously even in the absence of stress exposure^{9,11,12,16,24,25}. Whether this reflects a population survival strategy or a regulatory inaccuracy is unclear, as are the mechanisms underlying SOS heterogeneity. Here, we developed a single-molecule imaging approach based on a HaloTag fusion to directly monitor LexA inside live Escherichia coli cells, demonstrating the existence of 3 main states of LexA: DNA-bound stationary molecules, free LexA and degraded LexA species. These analyses elucidate the mechanisms by which DNA-binding and degradation of LexA regulate the SOS response in vivo. We show that self-cleavage of LexA occurs frequently throughout the population during unperturbed growth, rather than being restricted to a subpopulation of cells, which causes substantial cell-to-cell variation in LexA abundances. LexA variability underlies SOS gene expression heterogeneity and triggers spontaneous SOS pulses, which enhance bacterial survival in anticipation of stress.

The Parameter-Fitness Landscape of lexA Autoregulation in Escherichia coli

Feedback mechanisms are critical to control physiological responses. In gene regulation, one important example, termed negative autoregulation (NAR), occurs when a transcription factor (TF) inhibits its own production. NAR is common across the tree of life, enabling rapid homeostatic control of gene expression. NAR behavior can be described in accordance with its core biochemical parameters, but how constrained these parameters are by evolution is unclear. Here, we describe a model genetic network controlled by an NAR circuit within the bacterium Escherichia coli and elucidate these constraints by experimentally changing a key parameter and measuring its effect on circuit response and fitness. This analysis yielded a parameter-fitness landscape representing the genetic network, providing a window into what gene-environment conditions favor evolution of this regulatory strategy.

Feedback mechanisms are fundamental to the control of physiological responses. One important example in gene regulation, termed negative autoregulation (NAR), occurs when a transcription factor (TF) inhibits its own production through transcriptional repression. This enables more-rapid homeostatic control of gene expression. NAR circuits presumably evolve to limit the fitness costs of gratuitous gene expression. The key biochemical reactions of NAR can be parameterized using a mathematical model of promoter activity; however, this model of NAR has been studied mostly in the context of synthetic NAR circuits that are disconnected from the target genes of the TFs. Thus, it remains unclear how constrained NAR parameters are in a native circuit context, where the TF target genes can have fitness effects on the cell. To quantify these constraints, we created a panel of Escherichia coli strains with different lexA-NAR circuit parameters and analyzed the effect on SOS response function and bacterial fitness. Using a mathematical model for NAR, these experimental data were used to calculate NAR parameter values and derive a parameter-fitness landscape. Without feedback, survival of DNA damage was decreased due to high LexA concentrations and slower SOS “turn-on” kinetics. However, we show that, even in the absence of DNA damage, the lexA promoter is strong enough that, without feedback, high levels of lexA expression result in a fitness cost to the cell. Conversely, hyperfeedback can mimic lexA deletion, which is also costly. This work elucidates the lexA-NAR parameter values capable of balancing the cell’s requirement for rapid SOS response activation with limiting its toxicity.

IMPORTANCE Feedback mechanisms are critical to control physiological responses. In gene regulation, one important example, termed negative autoregulation (NAR), occurs when a transcription factor (TF) inhibits its own production. NAR is common across the tree of life, enabling rapid homeostatic control of gene expression. NAR behavior can be described in accordance with its core biochemical parameters, but how constrained these parameters are by evolution is unclear. Here, we describe a model genetic network controlled by an NAR circuit within the bacterium Escherichia coli and elucidate these constraints by experimentally changing a key parameter and measuring its effect on circuit response and fitness. This analysis yielded a parameter-fitness landscape representing the genetic network, providing a window into what gene-environment conditions favor evolution of this regulatory strategy.

Interconversion between bound and free conformations of LexA orchestrates the bacterial SOS response

The bacterial SOS response is essential for the maintenance of genomes, and also modulates antibiotic resistance and controls multidrug tolerance in subpopulations of cells known as persisters. In Escherichia coli, the SOS system is controlled by the interplay of the dimeric LexA transcriptional repressor with an inducer, the active RecA filament, which forms at sites of DNA damage and activates LexA for self-cleavage. Our aim was to understand how RecA filament formation at any chromosomal location can induce the SOS system, which could explain the mechanism for precise timing of induction of SOS genes. Here, we show that stimulated self-cleavage of the LexA repressor is prevented by binding to specific DNA operator targets. Distance measurements using pulse electron paramagnetic resonance spectroscopy reveal that in unbound LexA, the DNA-binding domains sample different conformations. One of these conformations is captured when LexA is bound to operator targets and this precludes interaction by RecA. Hence, the conformational flexibility of unbound LexA is the key element in establishing a co-ordinated SOS response. We show that, while LexA exhibits diverse dissociation rates from operators, it interacts extremely rapidly with DNA target sites. Modulation of LexA activity changes the occurrence of persister cells in bacterial populations.

Structure of the LexA-DNA complex and implications for SOS box measurement

The eubacterial SOS system is a paradigm of cellular DNA damage and repair, and its activation can contribute to antibiotic resistance. Under normal conditions, LexA represses the transcription of many DNA repair proteins by binding to SOS 'boxes' in their operators. Under genotoxic stress, accumulating complexes of RecA, ATP and single-stranded DNA (ssDNA) activate LexA for autocleavage. To address how LexA recognizes its binding sites, we determined three crystal structures of Escherichia coli LexA in complex with SOS boxes. Here we report the structure of these LexA-DNA complexes. The DNA-binding domains of the LexA dimer interact with the DNA in the classical fashion of a winged helix-turn-helix motif. However, the wings of these two DNA-binding domains bind to the same minor groove of the DNA. These wing-wing contacts may explain why the spacing between the two half-sites of E. coli SOS boxes is invariant.

Mifepristone-inducible LexPR system to drive and control gene expression in transgenic zebrafish

Effective transgenesis methods have been successfully employed in many organisms including zebrafish. However, accurate spatiotemporal control of transgene expression is still difficult to achieve. Here we describe a system for chemical-inducible gene expression and demonstrate its feasibility for generating transgenic driver lines in zebrafish. The key element of this system is a hybrid transcription factor engineered by fusion of the DNA-binding domain of the bacterial LexA repressor, a truncated ligand-binding domain of the human progesterone receptor, and the activation domain of the human NF-kappaB/p65 protein. This hybrid transcription factor (LexPR transactivator) binds to the synthetic steroid, mifepristone (RU-486), and functions in a ligand-dependent manner to induce expression of the gene(s) placed under the control of a synthetic operator-promoter sequence that harbors LexA binding sites. Transgene expression is strictly controlled and can be induced at any stage of the life cycle through administration of mifepristone in the water. To demonstrate the utility of this system, we generated stable transgenic lines which allow inducible tissue-specific expression of activated K-ras(V12). Combined with the Ac/Ds-mediated transgenesis, the LexPR expression system has many potential applications in the fields of genetics and biotechnology.

Thermodynamic and kinetic basis for the relaxed DNA sequence specificity of "promiscuous" mutant EcoRI endonucleases

Promiscuous mutant EcoRI endonucleases produce lethal to sublethal effects because they cleave Escherichia coli DNA despite the presence of the EcoRI methylase. Three promiscuous mutant forms, Ala138Thr, Glu192Lys and His114Tyr, have been characterized with respect to their binding affinities and first-order cleavage rate constants towards the three classes of DNA sites: specific, miscognate (EcoRI*) and non-specific. We have made the unanticipated and counterintuitive observations that the mutant restriction endonucleases that exhibit relaxed specificity in vivo nevertheless bind more tightly than the wild-type enzyme to the specific recognition sequence in vitro, and show even greater preference for binding to the cognate GAATTC site over miscognate sites. Binding preference for EcoRI* over non-specific DNA is also improved. The first-order cleavage rate constants of the mutant enzymes are normal for the cognate site GAATTC, but are greater than those of the wild-type enzyme at EcoRI* sites. Thus, the mutant enzymes use two mechanisms to partially bypass the multiple fail-safe mechanisms that protect against cleavage of genomic DNA in cells carrying the wild-type EcoRI restriction-modification system: (a) binding to EcoRI* sites is more probable than for wild-type enzyme because non-specific DNA is less effective as a competitive inhibitor; (b) the combination of increased affinity and elevated cleavage rate constants at EcoRI* sites makes double-strand cleavage of these sites a more probable outcome than it is for the wild-type enzyme. Semi-quantitative estimates of rates of EcoRI* site cleavage in vivo, predicted using the binding and cleavage constants measured in vitro, are in accord with the observed lethal phenotypes associated with the three mutations.

Improved model of a LexA repressor dimer bound to recA operator

A complete three dimensional model for the LexA repressor dimer bound to the recA operator site consistent with relevant biochemical and biophysical data for the repressor was proposed from our laboratory when no crystal structure of LexA was available. Subsequently, the crystal structures of four LexA mutants Delta(1-67) S119A, S119A, G85D and Delta(1-67) quadruple mutant in the absence of operator were reported. It is examined in this paper to what extent our previous model was correct and how, using the crystal structure of the operator-free LexA dimer we can predict an improved model of LexA dimer bound to recA operator. In our improved model, the C-domain dimerization observed repeatedly in the mutant operator-free crystals is retained but the relative orientation between the two domains within a LexA molecule changes. The crystal structure of wild type LexA with or without the recA operator cannot be solved as it autocleaves itself. We argue that the 'cleavable' cleavage site region found in the crystal structures is actually the more relevant form of the region in wild-type LexA since it agrees with the value of the pre-exponential Arrhenius factor for its autocleavage, absence of various types of trans-cleavages, difficulty in modifying the catalytic serine by diisopropyl flourophosphate and lack of cleavage at Arg 81 by trypsin; hence the concept of a 'conformational switch' inferred from the crystal structures is meaningless.

Model of a LexA repressor dimer bound to recA operator

A complete three dimensional model (RCSB000408; PDB code 1qaa) for the LexA repressor dimer bound to the recA operator site consistent with relevant biochemical and biophysical data for the repressor is proposed. A model of interaction of the N-terminal operator binding domain 1-72 with the operator was available. We have modelled residues 106-202 of LexA on the basis of the crystal structure of a homologous protein, UmuD'. Residues 70-105 have been modelled by us, residues 70-77 comprising the real hinge, followed by a beta-strand and an alpha-helix, both interacting with the rest of the C-domain. The preexponential Arrhenius factor for the LexA autocleavage is shown to be approximately 10(9) s(-1) at 298K whereas the exponential factor is approximately 2 x 10(-12), demanding that the autocleavage site is quite close to the catalytic site but reaction is slow due to an activation energy barrier. We propose that in the operator bound form, Ala 84- Gly 85 is about 7-10A from the catalytic groups, but the reaction does not occur as the geometry is not suitable for a nucleophilic attack from Ser 119 Ogamma, since Pro 87 is held in the cis conformation. When pH is elevated or under the action of activated RecA, cleavage may occur following a cis --> trans isomerization at Pro 87 and/or a rotation of the region beta9-beta10 about beta7-beta8 following the disruption of two hydrogen bonds. We show that the C-C interaction comprises the approach of two negatively charged surfaces neutralized by sodium ions, the C-domains of the monomers making a new beta barrel at the interface burying 710A2 of total surface area of each monomer.

Technical advance: An estrogen receptor-based transactivator XVE mediates highly inducible gene expression in transgenic plants

We have developed an estrogen receptor-based chemical-inducible system for use in transgenic plants. A chimeric transcription activator, XVE, was assembled by fusion of the DNA-binding domain of the bacterial repressor LexA (X), the acidic transactivating domain of VP16 (V) and the regulatory region of the human estrogen receptor (E; ER). The transactivating activity of the chimeric XVE factor, whose expression was controlled by the strong constitutive promoter G10-90, was strictly regulated by estrogens. In transgenic Arabidopsis and tobacco plants, estradiol-activated XVE can stimulate expression of a GFP reporter gene controlled by the target promoter, which consists of eight copies of the LexA operator fused upstream of the -46 35S minimal promoter. Upon induction by estradiol, GFP expression levels can be eightfold higher than that transcribed from a 35S promoter, whereas the uninduced controls have no detectable GFP transcripts, as monitored by Northern blot analysis. Neither toxic nor adverse physiological effects of the XVE system have been observed in transgenic Arabidopsis plants under all the conditions tested. The XVE system thus appears to be a reliable and efficient chemical-inducible system for regulating transgene expression in plants.

A LexA mutant repressor with a relaxed inter-domain linker

The LexA protein is part of a large family of prokaryotic transcriptional repressors that contain an amino-terminal DNA binding domain and a carboxy-terminal dimerization domain. These domains are separated by a linker or hinge region, which is generally considered to be rather flexible and unconstrained. So far, no structure of any of the full-length repressors is available. Here we show that a mutant LexA repressor harboring several point mutations in the hinge region gets sensitive to trypsin and Glu-C cleavage over a segment of at least 20 amino acids, whereas the LexA wild-type hinge region is resistant to these proteases. These data are not compatible with the hypothesis of an fully flexible and/or unstructured inter-domain linker and suggest that the LexA hinge region is, in fact, constrained by contacts with the carboxy-terminal domain and/or a fairly stable local structure of the linker region.

Preferential interactions of the Escherichia coli LexA repressor with anions and protons are coupled to binding the recA operator

The binding of Escherichia coli LexA repressor to the recA operator was examined as a function of the concentration of NaCl, KCl, NaF, and MgCl2 at pH 7.5, 21 degrees C. The effects of pH at 100 mM NaCl were also examined. Changes both in the qualitative appearance of the binding isotherms and in the magnitude of the apparent binding affinity with changes in solution conditions suggest that binding of anions and protons by LexA repressor is linked to oligomerization and/or operator binding. Binding of LexA repressor to the recA operator in the presence of NaCl ranging from 25 to 400 mM at picomolar DNA concentration showed a broad, apparently noncooperative, binding isotherm. Binding of LexA repressor in NaF at the same [DNA] yielded binding isotherms with a narrow transition, reflecting an apparently cooperative binding process. Also, the apparent binding affinity was weaker in NaF than in NaCl. Furthermore, the binding affinity and also the apparent binding mode, cooperative vs noncooperative, were pH dependent. The binding affinity of LexA repressor for operator was greatest near neutral pH. The apparent binding mode was noncooperative at pH 7-9 but was cooperative at pH 6 or 9.3. These observations suggest that the specific cation and anion composition and concentrations must be considered in understanding the details of regulation of the SOS system.

Binding of Escherichia coli LexA repressor to the RecA operator

Equilibrium binding of Escherichia coli LexA repressor to the recA operator was studied by the polyacrylamide gel mobility shift assay as a function of solution conditions. In the presence of NaCl at 20 degrees C, there was a significant salt dependence in binding to the recA operator, typical for protein-nucleic acid interactions with some electrostatic contribution to the binding free energy. In preliminary experiments in which the anion of the Na+ salt was changed from chloride to fluoride, little change was found with anion identity. This indicates that the salt effect on the binding interaction arises solely from the polyelectrolyte effect, not from anion binding or release by the protein upon complex formation. Increasing the temperature to 37 degrees C changed the binding affinity for complex formation at any given salt concentration and resulted in a change in the sensitivity of complex formation to NaCl concentration. Quantitative analysis of the data to obtain equilibrium binding constants is discussed.

Fos leucine zipper variants with increased association capacity

The Fos wild-type leucine zipper is unable to support homodimerization. This finding is generally explained by the negative net charge of the Fos zipper leading to the electrostatic repulsion of two monomers. Using a LexA-dependent in vivo assay in Escherichia coli, we show here that additional antideterminants for Fos zipper association are the residues in position a within the Fos zipper interface. If the wild-type Fos zipper is fused to the DNA binding domain of the LexA repressor (LexA-DBD), no excess repression is observed as compared with the LexA-DBD alone, in agreement with the incapacity of the wild-type Fos zipper to promote homodimerization. If hydrophobic amino acids (Ile, Leu, Val, Phe, Met) are inserted into the five a positions of a LexA-Fos zipper fusion protein, substantial transcriptional repression is recovered showing that Fos zipper homodimerization is not only limited by the repulsion of negatively charged residues but also by the nonhydrophobic nature of the a positions. The most efficient variants (harboring Ile or Leu in the five a positions) show an about 80-fold increase in transcriptional repression as compared with the wild-type Fos zipper fusion protein. In the case of multiple identical substitutions, the overall improvement is correlated with the hydrophobicity of the inserted side chains, i.e. Ile Leu > Val > Phe > Met. However at least for Val, Phe, and Met the impact of a given residue type on the association efficiency depends strongly on the heptad, i.e. on the local environment of the a residue. This is particularly striking for the second heptad of the Fos zipper, where Val is less well tolerated than Phe and Met. Most likely the a1 residue modulates the interhelical repulsion between two glutamic acid side chains in positions g1 and e2. Most of the hydrophobic Fos zipper variants are also improved in heteroassociation with a Jun leucine zipper, such that roughly half of the additional free energy of homodimerization is imported into the heterodimer. A few candidates (including the Fos wild-type zipper) deviate from this correlation, showing considerable excess heteroassociation.

A model for the LexA repressor DNA complex

A structural model for the interaction of the LexA repressor DNA binding domain (DBD) with operator DNA is derived by means of Monte Carlo docking. Protein-DNA complexes were generated by docking the LexA repressor DBD NMR solution structure onto both rigid and bent B-DNA structures while giving energy bonuses for contacts in agreement with experimental data. In the resulting complexes, helix III of the LexA repressor DBD is located in the major groove of the DNA and residues Asn-41, Glu-44, and Glu-45 form specific hydrogen bonds with bases of the CTGT DNA sequence. Ser-39, Ala-42, and Asn-41 are involved in a hydrophobic interaction with the methyl group of the first thymine base. Residues in the loop region connecting the two beta-sheet strands are involved in nonspecific contacts near the dyad axis of the operator. The contacts observed in the docked complexes cover the entire consensus CTGT half-site DNA operator, thus explaining the specificity of the LexA repressor for such sequences. In addition, a large number of nonspecific interactions between protein and DNA is observed. The agreement between the derived model for the LexA repressor DBD/DNA complex and experimental biochemical results is discussed.

Control of the LexA regulon by pH: evidence for a reversible inactivation of the LexA repressor during the growth cycle of Escherichia coli

The LexA repressor controls the expression of several genes, including lexA, recA, and sfiA, which are induced when exponentially growing bacteria are exposed to DNA-damaging agents. Induction of this so-called SOS response takes place while LexA is cleaved in a reaction that requires the RecA protein and damaged DNA. We have shown that large fluctuations in the cellular concentration of the LexA repressor and in the rate of transcription of the sfiA gene also occur spontaneously during bacterial growth in complex medium such as LB. The possibility that changes in external or internal pH may explain these fluctuations has been explored. A consistent pattern was established whereby conditions leading to either increased or decreased pH were associated with altered expression of the lexA and sfiA genes. These data can be explained by a model in which the LexA repressor exists in either of two forms in equilibrium: a form favoured at homeostatic internal pH, which has a low affinity for the operators of LexA-controlled genes; and a form accumulated in response to a transient decrease in internal pH, which has a high affinity for operators.

LEF-1 contains an activation domain that stimulates transcription only in a specific context of factor-binding sites

Lymphoid enhancer factor 1 (LEF-1) is a member of the high mobility group (HMG) family of proteins and participates in the regulation of the T cell receptor (TCR) alpha enhancer. We have previously shown that DNA binding by the HMG domain of LEF-1 induces a sharp bend in the DNA helix. Together with the dependence of LEF-1 on other factor-binding sites to regulate gene expression, DNA bending induced by the HMG domain suggested an 'architectural' role for LEF-1. In this study, we performed experiments to distinguish between a model in which the HMG domain is the only functional determinant of LEF-1 and a model in which additional domains of LEF-1 are involved in the regulation of gene expression. First, we show that the HMG domain alone is not sufficient to stimulate TCR alpha enhancer function. Second, we replaced the HMG domain of LEF-1 with the DNA-binding domain of the bacterial repressor LexA, which binds a specific nucleotide sequence without inducing a sharp bend in the DNA helix. The chimeric LEF-LexA protein increased the activity of a TCR alpha enhancer in which the LEF-1-binding site had been replaced with a LexA recognition sequence. Transcriptional stimulation by LEF-LexA, however, was less efficient than that observed with endogenous LEF-1. The LEF-LexA-mediated activation of gene expression was dependent upon an amino-terminal region of LEF-1 and a specific context of factor-binding sites in the TCR alpha enhancer. Neither multimerized LexA-binding sites, nor TCR alpha enhancers with altered spatial arrangements of factor-binding sites, were functional for regulation by LEF-LexA. Together, these data suggest that an aminoterminal region in LEF-1 contributes to the context-dependent regulation of the TCR alpha enhancer by LEF-1, presumably by interacting with other enhancer-bound proteins.