The role of rigidity in DNA looping-unlooping by AraC

A Career's Work, the l-Arabinose Operon: How It Functions and How We Learned It

Very few labs have had the good fortune to have been able to focus for more than 50 years on a relatively narrow research topic and to be in a field in which both basic knowledge and the research technology and methods have progressed as rapidly as they have in molecular biology. My research group, first at Brandeis University and then at Johns Hopkins University, has had this opportunity. In this review, therefore, I will describe largely the work from my laboratory that has spanned this period and which was carried out by 40 plus graduate students, several postdoctoral associates, my technician, and me. In addition to presenting the scientific findings or results, I will place many of the topics in scientific context and, because we needed to develop a good many of the experimental methods behind our findings, I will also describe some of these methods and their importance. Also included will be occasional comments on how the research community or my research group functioned. Because a wide variety of approaches were used throughout our work, no ideal organization of this review is apparent. Therefore, I have chosen to use a hybrid structure in which there are six sections. Within each of the sections, experiments and findings will be described roughly in chronological order. Frequent cross references between parts and sections will be made because some findings and experimental approaches could logically have been described in more than one place.

Single-molecule evidence for a chemical ratchet in binding between the cam repressor and its operator

The affinity for regulator-operator binding on DNA sometimes depends on the length of the DNA harboring the operator, which is known as the antenna effect. One-dimensional diffusion along DNA has been suggested to be the cause, but this may contradict the binding affinity independent of the reaction pathways, which is derived from the detailed balance of the reaction at equilibrium. Recently, the chemical ratchet was proposed to solve this contradiction by suggesting a stationary state containing microscopic non-equilibrium. In a single-molecule observation, P. putida CamR molecules associate with their operator via one-dimensional diffusion along the DNA, while they mostly dissociated from the operator without the diffusion. Consistently, the observed overall association rate was dependent on the DNA length, while the overall dissociation rate was not, leading to an antenna effect. E. coli RNA polymerase did not show this behavior, and thus it is a specific property of a protein. The bipartite interaction domains containing the helix-turn-helix motif are speculated to be one of the possible causes. The biological significance of the chemical ratchet and a model for its microscopic mechanism are also discussed.

Co-evolution of drug resistance and broadened substrate recognition in HIV protease variants isolated from an Escherichia coli genetic selection system

A genetic selection system for activity of HIV protease is described that is based on a synthetic substrate constructed as a modified AraC regulatory protein that when cleaved stimulate L-arabinose metabolism in an Escherichia coli araC strain. Growth stimulation on selective plates was shown to depend on active HIV protease and the scissile bond in the substrate. In addition, the growth of cells correlated well with the established cleavage efficiency of the sites in the viral polyprotein, Gag, when these sites were individually introduced into the synthetic substate of the selection system. Plasmids encoding protease variants selected based on stimulation of cell growth in the presence of saquinavir or cleavage of a site not cleaved by wild-type protease, were indistinguishable with respect to both phenotypes. Also, both groups of selected plasmids encoded side chain substitutions known from clinical isolates or displayed different side chain substitutions but at identical positions. One highly frequent side chain substitution, E34V, not regarded as a major drug resistance substitution was found in variants obtained under both selective conditions and is suggested to improve protease processing of the synthetic substrate. This substitution is away from the substrate-binding cavity and together with other substitutions in the selected reading frames supports the previous suggestion of a substrate-binding site extended from the active site binding pocket itself.

Prediction of whole-cell transcriptional response with machine learning

MOTIVATION

Applications in synthetic and systems biology can benefit from measuring whole-cell response to biochemical perturbations. Execution of experiments to cover all possible combinations of perturbations is infeasible. In this paper, we present the host response model (HRM), a machine learning approach that maps response of single perturbations to transcriptional response of the combination of perturbations.

RESULTS

The HRM combines high-throughput sequencing with machine learning to infer links between experimental context, prior knowledge of cell regulatory networks, and RNASeq data to predict a gene's dysregulation. We find that the HRM can predict the directionality of dysregulation to a combination of inducers with an accuracy of >90% using data from single inducers. We further find that the use of prior, known cell regulatory networks doubles the predictive performance of the HRM (an R2 from 0.3 to 0.65). The model was validated in two organisms, Escherichia coli and Bacillus subtilis, using new experiments conducted after training. Finally, while the HRM is trained with gene expression data, the direct prediction of differential expression makes it possible to also conduct enrichment analyses using its predictions. We show that the HRM can accurately classify >95% of the pathway regulations. The HRM reduces the number of RNASeq experiments needed as responses can be tested in silico prior to the experiment.

AVAILABILITY AND IMPLEMENTATION

The HRM software and tutorial are available at https://github.com/sd2e/CDM and the configurable differential expression analysis tools and tutorials are available at https://github.com/SD2E/omics_tools.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

One-dimensional diffusion of TrpR along DNA enhances its affinity for the operator by chemical ratchet mechanism

Several DNA-binding proteins show the affinities for their specific DNA sites that positively depend on the length of DNA harboring the sites, i. e. antenna effect. DNA looping can cause the effect for proteins with two or more DNA binding sites, i. e. the looping mechanism. One-dimensional diffusion also has been suggested to cause the effect for proteins with single DNA sites, the diffusion mechanism, which could violate detailed balance. We addressed which mechanism is possible for E. coli TrpR showing 10⁴-fold antenna effect with a single DNA binding site. When a trpO-harboring DNA fragment was connected to a nonspecific DNA with biotin-avidin connection, the otherwise sevenfold antenna effect disappeared. This result denies the looping mechanism with an unknown second DNA binding site. The 3.5-fold repression by TrpR in vivo disappeared when a tight LexA binding site was introduced at various sites near the trpO, suggesting that the binding of LexA blocks one-dimensional diffusion causing the antenna effect. These results are consistent with the chemical ratchet recently proposed for TrpR-trpO binding to solve the deviation from detailed balance, and evidence that the antenna effect due to one-dimensional diffusion exists in cells.

RNA Polymerase and Transcription Mechanisms: The Forefront of Physicochemical Studies of Chemical Reactions

The study of transcription and its regulation is an interdisciplinary field that is closely connected with genetics, structural biology, and reaction theory. Among these, although less attention has been paid to reaction theory, it is becoming increasingly useful for research on transcription. Rate equations are commonly used to describe reactions involved in transcription, but they tend to be used unaware of the timescales of relevant physical processes. In this review, we discuss the limitation of rate equation for describing three-dimensional diffusion and one-dimensional diffusion along DNA. We then introduce the chemical ratchet mechanism recently proposed for explaining the antenna effect, an enhancement of the binding affinity to a specific site on longer DNA, which deviates from a thermodynamic rule. We show that chemical ratchet cannot be described with a single set of rate equations but alternative sets of rate equations that temporally switch no faster than the binding reaction.

Arabinose Alters Both Local and Distal H-D Exchange Rates in the Escherichia coli AraC Transcriptional Regulator

In the absence of arabinose, the dimeric Escherichia coli regulatory protein of the l-arabinose operon, AraC, represses expression by looping the DNA between distant half-sites. Binding of arabinose to the dimerization domains forces AraC to preferentially bind two adjacent DNA half-sites, which stimulates RNA polymerase transcription of the araBAD catabolism genes. Prior genetic and biochemical studies hypothesized that arabinose allosterically induces a helix-coil transition of a linker between the dimerization and DNA binding domains that switches the AraC conformation to an inducing state [Brown, M. J., and Schleif, R. F. (2019) Biochemistry, preceding paper in this issue (DOI: 10.1021/acs.biochem.9b00234)]. To test this hypothesis, hydrogen-deuterium exchange mass spectrometry was utilized to identify structural regions involved in the conformational activation of AraC by arabinose. Comparison of the hydrogen-deuterium exchange kinetics of individual dimeric dimerization domains and the full-length dimeric AraC protein in the presence and absence of arabinose reveals a prominent arabinose-induced destabilization of the amide hydrogen-bonded structure of linker residues (I₁₆₇ and N₁₆₈). This destabilization is demonstrated to result from an increased probability to form a helix capping motif at the C-terminal end of the dimerizing α-helix of the dimerization domain that preceeds the interdomain linker. These conformational changes could allow for quaternary repositioning of the DNA binding domains required for induction of the araBAD promoter through rotation of peptide backbone dihedral angles of just a couple of residues. Subtle changes in exchange rates are also visible around the arabinose binding pocket and in the DNA binding domain.

The mechanics of DNA loops bridged by proteins unveiled by single-molecule experiments

Protein-induced DNA bridging and looping is a common mechanism for various and essential processes in bacterial chromosomes. This mechanism is preserved despite the very different bacterial conditions and their expected influence on the thermodynamic and kinetic characteristics of the bridge formation and stability. Over the last two decades, single-molecule techniques carried out on in vitro DNA systems have yielded valuable results which, in combination with theoretical works, have clarified the effects of different parameters of nucleoprotein complexes on the protein-induced DNA bridging and looping process. In this review, I will outline the features that can be measured for such processes with various single-molecule techniques in use in the field. I will then describe both the experimental results and the theoretical models that illuminate the contribution of the DNA molecule itself as well as that of the bridging proteins in the DNA looping mechanism at play in the nucleoid of E. coli.

Modulation of DNA binding by gene-specific transcription factors

The transcription of many genes, particularly in prokaryotes, is controlled by transcription factors whose activity can be modulated by controlling their DNA binding affinity. Understanding the molecular mechanisms by which DNA binding affinity is regulated is important, but because forming definitive conclusions usually requires detailed structural information in combination with data from extensive biophysical, biochemical, and sometimes genetic experiments, little is truly understood about this topic. This review describes the biological requirements placed upon DNA binding transcription factors and their consequent properties, particularly the ways that DNA binding affinity can be modulated and methods for its study. What is known and not known about the mechanisms modulating the DNA binding affinity of a number of prokaryotic transcription factors, including CAP and lac repressor, is provided.

Understanding the basis of a class of paradoxical mutations in AraC through simulations

Most mutations at position 15 in the N-terminal arm of the regulatory protein AraC leave the protein incapable of responding to arabinose and inducing the proteins required for arabinose catabolism. Mutations at other positions of the arm do not have this behavior. Simple energetic analysis of the interactions between the arm and bound arabinose do not explain the uninducibility of AraC with mutations at position 15. Extensive molecular dynamics (MD) simulations, carried out largely on the Open Science Grid, were done of the wild-type protein with and without bound arabinose and of all possible mutations at position 15, many of which were constructed and measured for this work. Good correlation was found for deviation of arm position during the simulations and inducibility as measured in vivo of the same mutant proteins. Analysis of the MD trajectories revealed that preservation of the shape of the arm is critical to inducibility. To maintain the correct shape of the arm, the strengths of three interactions observed to be strong in simulations of the wild-type AraC protein need to be preserved. These interactions are between arabinose and residue 15, arabinose and residues 8-9, and residue 13 and residue 15. The latter interaction is notable because residues L9, Y13, F15, W95, and Y97 form a hydrophobic cluster which needs to be preserved for retention of the correct shape.

Mutational analysis of the N-terminal domain of UreR, the positive transcriptional regulator of urease gene expression

The Escherichia coli plasmid-encoded urease, a virulence factor in human and animal infections of the urinary and gastroduodenal tracts, is induced when the substrate urea is present in the growth medium. Urea-dependent urease expression is mediated at the transcriptional level by the AraC-like activator UreR. Previous work has shown that a peptide representing the N-terminal 194 amino-acid residues of UreR binds urea at a single site, full-length UreR forms an oligomer, and the oligomerization motif is thought to reside in the N-terminal portion of the molecule. The C-terminal domain of UreR contains two helix-turn-helix motifs presumed to be necessary for DNA binding. In this study, we exploited mutational analyses at the N-terminal domain of UreR to determine if this domain dimerizes similar to other AraC family members. UreR mutants were analyzed for the ability to activate transcription of lacZ from an ureDp-lacZ transcriptional fusion. A construct encoding the N-terminal 194 amino acids of UreR, eluted as an oligomer by gel filtration and had a dominant negative phenotype over the wild-type ureR allele. We hypothesize that this dominant negative phenotype results from the formation of inactive heterodimers between wild-type and truncated UreR. Dominant negative analysis and cross-linking assays demonstrated that E. coli UreR is active as a dimer and dimerization occurs within the first 180 residues.

Active role of the interdomain linker of AraC

Mutations in the interdomain linker of the gene for the AraC regulatory protein of Escherichia coli that severely interfere with the protein's ability either to repress or to activate transcription have been found. These mutations have relatively small effects on the dimerization domain's ability to bind arabinose or to dimerize the protein or on the DNA-binding domain's affinity for a single DNA half-site. The linker mutations, however, dramatically change the affinity of AraC for binding to two direct-repeat DNA half-sites. Less dramatically, the induction-deficient linker variants also display altered DNA sequence selectivity. These results show that changing the sequence of the interdomain linker can profoundly affect the dimerization domain-DNA-binding domain interactions in AraC. The smaller effects on the functions of the individual domains could be the direct result of the linker alterations but more likely are the indirect result of the altered dimerization domain-DNA-binding domain interactions. In summary, the linker does not simply function as a passive and flexible connector between the domains of AraC but, instead, is more directly involved in the protein's dimerization domain-DNA-binding domain interactions.

Understanding apparent DNA flexibility enhancement by HU and HMGB architectural proteins

Understanding and predicting the mechanical properties of protein/DNA complexes are challenging problems in biophysics. Certain architectural proteins bind DNA without sequence specificity and strongly distort the double helix. These proteins rapidly bind and unbind, seemingly enhancing the flexibility of DNA as measured by cyclization kinetics. The ability of architectural proteins to overcome DNA stiffness has important biological consequences, but the detailed mechanism of apparent DNA flexibility enhancement by these proteins has not been clear. Here, we apply a novel Monte Carlo approach that incorporates the precise effects of protein on DNA structure to interpret new experimental data for the bacterial histone-like HU protein and two eukaryotic high-mobility group class B (HMGB) proteins binding to ∼200-bp DNA molecules. These data (experimental measurement of protein-induced increase in DNA cyclization) are compared with simulated cyclization propensities to deduce the global structure and binding characteristics of the closed protein/DNA assemblies. The simulations account for all observed (chain length and concentration dependent) effects of protein on DNA behavior, including how the experimental cyclization maxima, observed at DNA lengths that are not an integral helical repeat, reflect the deformation of DNA by the architectural proteins and how random DNA binding by different proteins enhances DNA cyclization to different levels. This combination of experiment and simulation provides a powerful new approach to resolve a long-standing problem in the biophysics of protein/DNA interactions.

Gene repression by minimal lac loops in vivo

The inflexibility of double-stranded DNA with respect to bending and twisting is well established in vitro. Understanding apparent DNA physical properties in vivo is a greater challenge. Here, we exploit repression looping with components of the Escherichia coli lac operon to monitor DNA flexibility in living cells. We create a minimal system for testing the shortest possible DNA repression loops that contain an E. coli promoter, and compare the results to prior experiments. Our data reveal that loop-independent repression occurs for certain tight operator/promoter spacings. When only loop-dependent repression is considered, fits to a thermodynamic model show that DNA twisting limits looping in vivo, although the apparent DNA twist flexibility is 2- to 4-fold higher than in vitro. In contrast, length-dependent resistance to DNA bending is not observed in these experiments, even for the shortest loops constraining <0.4 persistence lengths of DNA. As observed previously for other looping configurations, loss of the nucleoid protein heat unstable (HU) markedly disables DNA looping in vivo. Length-independent DNA bending energy may reflect the activities of architectural proteins and the structure of the DNA topological domain. We suggest that the shortest loops are formed in apical loops rather than along the DNA plectonemic superhelix.

AraC protein, regulation of the l-arabinose operon in Escherichia coli, and the light switch mechanism of AraC action

This review covers the physiological aspects of regulation of the arabinose operon in Escherichia coli and the physical and regulatory properties of the operon's controlling gene, araC. It also describes the light switch mechanism as an explanation for many of the protein's properties. Although many thousands of homologs of AraC exist and regulate many diverse operons in response to many different inducers or physiological states, homologs that regulate arabinose-catabolizing genes in response to arabinose were identified. The sequence similarities among them are discussed in light of the known structure of the dimerization and DNA-binding domains of AraC.

Computational predictions of the mutant behavior of AraC

An algorithm implemented in Rosetta correctly predicts the folding capabilities of the 17-residue N-terminal arm of the AraC gene regulatory protein when arabinose is bound to the protein and the dramatically different structure of this arm when arabinose is absent. The transcriptional activity of 43 mutant AraC proteins with alterations in the arm sequences was measured in vivo and compared with their predicted folding properties. Seventeen of the mutants possessed regulatory properties that could be directly compared with folding predictions. Sixteen of the 17 mutants were correctly predicted. The algorithm predicts that the N-terminal arm sequences of AraC homologs fold to the Escherichia coli AraC arm structure. In contrast, it predicts that random sequences of the same length and many partially randomized E. coli arm sequences do not fold to the E. coli arm structure. The high level of success shows that relatively "simple" computational methods can in some cases predict the behavior of mutant proteins with good reliability.

Solution structure of the DNA binding domain of AraC protein

We report the solution structure of the DNA binding domain of the Escherichia coli regulatory protein AraC determined in the absence of DNA. The 20 lowest energy structures, determined on the basis of 1507 unambiguous nuclear Overhauser restraints and 180 angle restraints, are well resolved with a pair wise backbone root mean square deviation of 0.7 A. The protein, free of DNA, is well folded in solution and contains seven helices arranged in two semi-independent sub domains, each containing one helix-turn-helix DNA binding motif, joined by a 19 residue central helix. This solution structure is discussed in the context of extensive biochemical and physiological data on AraC and with respect to the DNA-bound structures of the MarA and Rob homologs.

A DNA-assisted binding assay for weak protein-protein interactions

We describe a new method used for quantitating weak interactions between proteins in which the weak interaction is "assisted" by a known DNA-DNA interaction. Oligonucleotides, which are conjugated to proteins of interest, contain short complementary DNA sequences that provide additional binding energy for protein-protein interactions. A stretch of unpaired bases links the protein to the hybridizing DNA sequence to allow formation of both protein-protein and DNA-DNA interactions with minimal structural interference. We validated the DNA-assisted binding method using heterodimerizing coiled-coil proteins. The method was then used to measure the predicted weak interaction between two domains of the Escherichia coli L-arabinose operon regulatory protein AraC. The interaction between domains has the expected magnitude (K(d)=0.37 mM) in the absence of arabinose. Upon addition of arabinose, we detected a weaker and unexpected interaction, which may necessitate modification of the proposed mechanism of AraC. The DNA-assisted binding method may also prove useful in the study of other weak protein-protein interactions.

Opposite allosteric mechanisms in TetR and CAP

Regulation of the DNA binding affinity of an oligomeric protein can be considered to consist of an intrinsic component, in which the affinity of an individual DNA-binding domain is modulated in response to effector binding, and an extrinsic component, in which the relative position of the protein's two DNA-binding domains are altered so that they can or cannot contact both half-site operators simultaneously. We demonstrated directly that the TetR repressor utilizes an extrinsic mechanism and CAP, the catabolite activator protein, utilizes an intrinsic mechanism.

Constitutive mutations in the Escherichia coli AraC protein

The Escherichia coli AraC protein represses and induces the araBAD operon in response to the absence or presence of l-arabinose. Constitutive mutations in the AraC gene no longer require the presence of l-arabinose to convert AraC from its repressing to its inducing state. Such mutations were isolated directly by virtue of their constitutivity or by their resistance to the nonmetabolizable arabinose analog, d-fucose. The majority of the constitutive mutations lie within the same residues of the N-terminal regulatory arm of AraC. Two, however, were found in the core of the dimerization domain. As predicted by the light switch mechanism of AraC, constitutive mutations increase the susceptibility of the N-terminal arms to digestion by trypsin or chymotrypsin, suggesting that these mutations weaken or disrupt the arm structure required for repression by AraC. Fluorescence, circular dichroism, and cysteine reactivity measurements show that the constitutive mutations in the core of the dimerization domain lead to a weakening of the support for the arms and reduce the stability of the minus-arabinose arm structure. These mutations also weaken the interaction between the two-helix bundle and the beta-barrel subdomains of the dimerization domain and reduce the structural stability of the beta-barrels.

Functional modes of the regulatory arm of AraC

One of the two crystal structures of the arm-dimerization domain determined in the absence of arbinose fails to show the arm, whereas the other structure does show it. The two structures lead to different pictures for the regulatory behavior of the arms. Trypsin digestion, fluorescence anisotropy, and NMR experiments presented here were designed to resolve the issue and show that in arm-dimerization domain, the arms are structured, although differently, in the presence and absence of arabinose. The arms have also been shown to interact with the DNA binding domains of AraC by their requirement for the immobilization of the DNA binding domains that is necessary for DNA looping and repression. The binding of arabinose has been shown to release the DNA binding domains and looping ceases. The picture resulting from the new experiments and the crystal structures of the arm-dimerization domain is that in the absence of arabinose, the arm adopts one structure on the dimerization domain and that the DNA binding domain then binds to this complex. Upon binding arabinose, the arm restructures and as a result, no longer serves as a gasket between the DNA binding domain and dimerization domain. The DNA binding domain is then released, subject only to the constraints imposed by the flexible linker connecting it to dimerization domain, and the protein relocates on the DNA and activates transcription.

Eukaryotic HMGB proteins as replacements for HU in E. coli repression loop formation

DNA looping is important for gene repression and activation in Escherichia coli and is necessary for some kinds of gene regulation and recombination in eukaryotes. We are interested in sequence-nonspecific architectural DNA-binding proteins that alter the apparent flexibility of DNA by producing transient bends or kinks in DNA. The bacterial heat unstable (HU) and eukaryotic high-mobility group B (HMGB) proteins fall into this category. We have exploited a sensitive genetic assay of DNA looping in living E. coli cells to explore the extent to which HMGB proteins and derivatives can complement a DNA looping defect in E. coli lacking HU protein. Here, we show that derivatives of the yeast HMGB protein Nhp6A rescue DNA looping in E. coli lacking HU, in some cases facilitating looping to a greater extent than is observed in E. coli expressing normal levels of HU protein. Nhp6A-induced changes in the DNA length-dependence of repression efficiency suggest that Nhp6A alters DNA twist in vivo. In contrast, human HMGB2-box A derivatives did not rescue looping.

Differences in the mechanism of the allosteric l-rhamnose responses of the AraC/XylS family transcription activators RhaS and RhaR

Proteins in the largest subset of AraC/XylS family transcription activators, including RhaS and RhaR, have C-terminal domains (CTDs) that mediate DNA-binding and transcription activation, and N-terminal domains (NTDs) that mediate dimerization and effector binding. The mechanism of the allosteric effector response in this family has been identified only for AraC. Here, we investigated the mechanism by which RhaS and RhaR respond to their effector, l-rhamnose. Unlike AraC, N-terminal truncations suggested that RhaS and RhaR do not use an N-terminal arm to inhibit activity in the absence of effector. We used random mutagenesis to isolate RhaS and RhaR variants with enhanced activation in the absence of l-rhamnose. NTD substitutions largely clustered around the predicted l-rhamnose-binding pockets, suggesting that they mimic the structural outcome of effector binding to the wild-type proteins. RhaS-CTD substitutions clustered in the first HTH motif, and suggested that l-rhamnose induces improved DNA binding. In contrast, RhaR-CTD substitutions clustered at a single residue in the second HTH motif, at a position consistent with improved RNAP contacts. We propose separate allosteric mechanisms for the two proteins: Without l-rhamnose, RhaS does not effectively bind DNA while RhaR does not effectively contact RNAP. Upon l-rhamnose binding, both proteins undergo structural changes that enable transcription activation.

Geometry of mediating protein affects the probability of loop formation in DNA

Recent single molecule experiments have determined the probability of loop formation in DNA as a function of the DNA contour length for different types of looping proteins. The optimal contour length for loop formation as well as the probability density functions have been found to be strongly dependent on the type of looping protein used. We show, using Monte Carlo simulations and analytical calculations, that these observations can be replicated using the wormlike-chain model for double-stranded DNA if we account for the nonzero size of the looping protein. The simulations have been performed in two dimensions so that bending is the only mode of deformation available to the DNA while the geometry of the looping protein enters through a single variable which is representative of its size. We observe two important effects that seem to directly depend on the size of the enzyme: 1), the overall propensity of loop formation at any given value of the DNA contour length increases with the size of the enzyme; and 2), the contour length corresponding to the first peak as well as the first well in the probability density functions increases with the size of the enzyme. Additionally, the eigenmodes of the fluctuating shape of the looped DNA calculated from simulations and theory are in excellent agreement, and reveal that most of the fluctuations in the DNA occur in regions of low curvature.

DNA tape measurements of AraC

A new method for measuring distances between points in the AraC–DNA complex was developed and applied. It utilizes variable lengths of single-stranded DNA that connect double-stranded regions containing the two half-site binding sequences of AraC. These distances plus the protein interdomain linker distances are compatible with two classes of structure for the dimeric AraC gene regulatory protein. In one class, the N-terminal regulatory arm of one dimerization domain is capable of interacting with the DNA-binding domain on the same polypeptide chain for a cis interaction. In the other class, the possible arm-DNA-binding domain interaction is trans, where it adds to the dimerization interface.

Specific interactions by the N-terminal arm inhibit self-association of the AraC dimerization domain

Deletion of the regulatory N-terminal arms of the AraC protein from its dimerization domain fragments increases the susceptibility of the dimerization domain to form a series of higher order polymers by indefinite self-association. We investigated how the normal presence of the arm inhibits this self-association. One possibility is that arms can act as an entropic bristles to interfere with the approach of other macromolecules, thereby decreasing collision frequencies. We examined the repulsive effect of flexible arms by measuring the rate of trypsin cleavage of a specially constructed ubiquitin-arm protein. Adding an arm to ubiquitin or increasing its length produced only a modest repulsive effect. This suggests that arms such as the N-terminal arm of AraC do not reduce self-association by entropic exclusion. We consequently tested the hypothesis that the arm on AraC reduces self-association by binding to the core of the dimerization domain even in the absence of arabinose. The behaviors of dimerization domain mutants containing deletions or alterations in the N-terminal arms substantiate this hypothesis. Apparently, interactions between the N-terminal arm and the dimerization domain core position the arm to interfere with the protein-protein contacts necessary for self-association.

Structure and properties of a truely apo form of AraC dimerization domain

The arabinose-binding pockets of wild type AraC dimerization domains crystallized in the absence of arabinose are occupied with the side chains of Y31 from neighboring domains. This interaction leads to aggregation at high solution concentrations and prevents determination of the structure of truely apo AraC. In this work we found that the aggregation does not significantly occur at physiological concentrations of AraC. We also found that the Y31V mutation eliminates the self-association, but does not affect regulation properties of the protein. At the same time, the mutation allows crystallization of the dimerization domain of the protein with only solvent in the arabinose-binding pocket. Using a distance difference method suitable for detecting and displaying even minor structural variation among large groups of similar structures, we find that there is no significant structural change in the core of monomers of the AraC dimerization domain resulting from arabinose, fucose, or tyrosine occupancy of the ligand-binding pocket. A slight change is observed in the relative orientation of monomers in the dimeric form of the domain upon the binding of arabinose but its significance cannot yet be assessed.

Mutational analysis of the Escherichia coli melR gene suggests a two-state concerted model to explain transcriptional activation and repression in the melibiose operon

Transcription of the Escherichia coli melAB operon is regulated by the MelR protein, an AraC family member whose activity is modulated by the binding of melibiose. In the absence of melibiose, MelR is unable to activate the melAB promoter but autoregulates its own expression by repressing the melR promoter. Melibiose triggers MelR-dependent activation of the melAB promoter and relieves MelR-dependent repression of the melR promoter. Twenty-nine single amino acid substitutions in MelR that result in partial melibiose-independent activation of the melAB promoter have been identified. Combinations of different substitutions result in almost complete melibiose-independent activation of the melAB promoter. MelR carrying each of the single substitutions is less able to repress the melR promoter, while MelR carrying some combinations of substitutions is completely unable to repress the melR promoter. These results argue that different conformational states of MelR are responsible for activation of the melAB promoter and repression of the melR promoter. Supporting evidence for this is provided by the isolation of substitutions in MelR that block melibiose-dependent activation of the melAB promoter while not changing melibiose-independent repression of the melR promoter. Additional experiments with a bacterial two-hybrid system suggest that interactions between MelR subunits differ according to the two conformational states.

Single-molecule spectroscopic determination of lac repressor-DNA loop conformation

The Escherichia coli lactose repressor protein (LacI) provides a classic model for understanding protein-induced DNA looping. LacI has a C-terminal four-helix bundle tetramerization domain that may act as a flexible hinge. In previous work, several DNA constructs, each containing two lac operators bracketing a sequence-induced bend, were designed to stabilize different possible looping geometries. The resulting hyperstable LacI-DNA loops exist as both a compact "closed" form with a V-shaped repressor and also a more "open" form with an extended hinge. The "9C14" construct was of particular interest because footprinting, electrophoretic mobility shift, and ring closure experiments suggested that it forms both geometries. Previous fluorescence resonance energy transfer (FRET) measurements gave an efficiency of energy transfer (ET) of 70%, confirming the existence of a closed form. These measurements could not determine whether open form or intermediate geometries are populated or the timescale of interconversion. We have now applied single-molecule FRET to Cy3, Cy5 double-labeled LacI-DNA loops diffusing freely in solution. By using multiple excitation wavelengths and by carefully examining the behavior of the zero-ET peak during titration with LacI, we show that the LacI-9C14 loop exists exclusively in a single closed form exhibiting essentially 100% ET.

Bacterial repression loops require enhanced DNA flexibility

The Escherichia coli lac operon provides a classic paradigm for understanding regulation of gene transcription. It is now appreciated that lac promoter repression involves cooperative binding of the bidentate lac repressor tetramer to pairs of lac operators via DNA looping. We have adapted components of this system to create an artificial assay of DNA flexibility in E.coli. This approach allows for systematic study of endogenous and exogenous proteins as architectural factors that enhance apparent DNA flexibility in vivo. We show that inducer binding does not completely remove repression loops but it does alter their geometries. Deletion of the E.coli HU protein drastically destabilizes small repression loops, an effect that can be partially overcome by expression of a heterologous mammalian HMG protein. These results emphasize that the inherent torsional inflexibility of DNA restrains looping and must be modulated in vivo.

DNA twisting flexibility and the formation of sharply looped protein-DNA complexes

Gene-regulatory complexes often require that pairs of DNA-bound proteins interact by looping-out short (often approximately 100-bp) stretches of DNA. The loops can vary in detailed length and sequence and, thus, in total helical twist, which radically alters their geometry. How this variability is accommodated structurally is not known. Here we show that the inherent twistability of 89- to 105-bp DNA circles exceeds theoretical expectation by up to 400-fold. These results can be explained only by greatly enhanced DNA flexibility, not by permanent bends. They invalidate the use of classic theories of flexibility for understanding sharp DNA looping but support predictions of two recent theories. Our findings imply an active role for DNA flexibility in loop formation and suggest that variability in the detailed helical twist of regulatory loops is accommodated naturally by the inherent twistability of the DNA.

LeuO protein delimits the transcriptionally active and repressive domains on the bacterial chromosome

LeuO protein relieves bacterial gene silencer AT8-mediated transcriptional repression as part of a promoter relay mechanism found in the ilvIH-leuO-leuABCD gene cluster. The gene silencing activity has recently been characterized as a nucleoprotein filament initiated at the gene silencer. In this gene locus, the nucleoprotein filament cis-spreads toward the target leuO promoter and results in the repression of the leuO gene. Although the cis-spreading nature of the transcriptionally repressive nucleoprotein filament has been revealed, the mechanism underlying LeuO-mediated gene silencing relief remains unknown. We have demonstrated here that LeuO functions analogously to the eukaryotic boundary element that delimits the transcriptionally active and repressive domains on the chromosome by blocking the cis-spreading pathway of the transcriptionally repressive heterochromatin. Given that one LeuO-binding site is positioned between the gene silencer and the target promoter, the simultaneous presence of a second LeuO-binding site synergistically enhances the blockade, resulting in a cooperative increase in LeuO-mediated gene silencing relief. A known DNA loop-forming protein, the lac repressor (LacI), was used to confirm that cooperative protein binding via DNA looping is responsible for the blocking synergy. Indeed, a distal LeuO site located downstream cooperates with the LeuO sites located upstream of the leuO gene, resulting in synergistic relief for the repressed leuO gene via looping out the intervening DNA between LeuO sites in the ilvIH-leuO-leuABCD gene cluster.

Geometric and dynamic requirements for DNA looping, wrapping and unwrapping in the activation of E.coli glnAp2 transcription by NtrC

Transcriptional activation by the E.coli NtrC protein can occur via DNA looping between a DNA-bound activator and the target sigma(54) RNA polymerase. NtrC forms an octamer on DNA that is capable of binding two DNA molecules. Its ATPase activity is required for open complex formation. Geometric requirements for activation were assessed using a library of DNA bending sequences created by random ligation of A-tract oligonucleotides, as well as several designed sequences. Thirty random or designed sequences with a variety of DNA lengths and bending geometries were cloned in plasmids, and the library was used to replace the spacer between the NtrC binding sites and the core glnAp2 promoter. The activity of each promoter construct under nitrogen limitation was determined in vivo, in a lambda phage lacZ reporter system integrated as a single-copy lysogen to avoid titrating NtrC or polymerase. A wide variety of bending geometries was found to support a similar level of transcriptional activation ( approximately 3-4-fold). Computer modeling of the DNA trajectories suggests that the most inactive promoters have short spacer DNA and the NtrC sites on the opposite side of the helix as the wild-type sites; otherwise, the loop can form effectively. Flexibility and multivalency of the NtrC-Esigma(54) interaction apparently provides substantial independence from DNA stiffness constraints, and in general activation requires less efficient looping than repression. However, none of the random templates were as active as wild-type promoter. Subsidiary activator binding sites in the wild-type were found to be required for full activity, but, surprisingly, these sites could not be functionally replaced by strong binding sites. This suggests that one or more protomers in the NtrC octamer must form and then release contacts with DNA in order to complete the ATPase cycle and act as an AAA(+) activator of the Esigma(54). This dynamic DNA wrapping around the NtrC octamer is proposed to be necessary for efficient activation, and the wrapping may also reduce adventitious activation of other promoters.

Biochemical and physiological properties of the DNA binding domain of AraC protein

Intact AraC protein is poorly soluble and difficult to purify, whereas its dimerization domain is the opposite. Unexpectedly, the DNA binding domain of AraC proved also to be soluble in cells when overproduced and is easily purified to homogeneity. The DNA binding affinity of the DNA binding domain for its binding site could not be measured by electrophoretic mobility shift because of its rapid association and dissociation rates, but its affinity could be measured with a fluorescence assay and was found to have a dissociation constant of 1 x 10(-8)M in 100 mM KCl. The binding of monomers of the DNA binding domain to adjacent half-sites occurs without substantial positive or negative cooperativity. A simple analysis relates the DNA binding affinities of monomers of DNA binding domain and normal dimeric AraC protein.

Site-directed mutagenesis of conserved inverted repeat sequences in the xylanase C promoter region from Streptomyces sp. EC3

Streptomyces sp. EC3, a strain which was originally isolated from cattle manure compost, was shown to possess a strong xylanolytic activity. One of the genes responsible for this activity, xlnC, encodes a secreted xylanase. In the native strain, as in the heterologous host S. lividans, expression of xlnC was detectable in the presence of xylan but not in the presence of glucose. Induction by xylan was shown to take place at the transcriptional level. The transcriptional start site of xlnC was mapped and likely -35 (5'-TTGACA-3') and -10 (5'-GAGAAC-3') motifs were identified. In order to localise putative conserved regulatory sequences, the promoter regions of xylanase-encoding genes from various Streptomyces species were aligned. This alignment revealed the existence of three sets of quite well conserved palindromic AT rich sequences called boxes 1, 2 and 3. Box 3 (5'-CGAAA N TTTCG-3') is the farthest away from the promoter region (150-200 bp). A shorter version of this palindrome (5'-GAAA NN TTTC-3') or (5'-CGAAA-3') constitutes box 1, which is located just upstream of the putative -35 promoter sequence. Box 2, located 5-7 bp upstream of box 1, comprises a shorter palindrome than box 3, with inverted polarity [5'-(G/C)TTTC (N) GAAA(G/C)-3']. The putative regulatory role of the conserved inverted repeats in boxes 2 and 3 in the promoter region of the xlnC gene from Streptomyces sp. EC3, was assessed. These boxes were modified by site-directed mutagenesis, and the mutant promoter regions, as well as the wild-type promoter region, were separately fused to a beta-lactamase reporter gene. Analysis of the expression patterns of these fusions in cultures grown in the presence of glucose, xylan or both carbon sources demonstrated that these motifs were cis -acting negative regulatory elements, each playing a specific role in the regulation of xlnC expression. Box 3 was shown to be critical for the establishment of repression of xlnC expression by glucose, whereas box 2 was shown to play an important role in the induction of xlnC expression by xylan.

Phage N4 RNA polymerase II recruitment to DNA by a single-stranded DNA-binding protein

Transcription of bacteriophage N4 middle genes is carried out by a phage-coded, heterodimeric RNA polymerase (N4 RNAPII), which belongs to the family of T7-like RNA polymerases. In contrast to phage T7-RNAP, N4 RNAPII displays no activity on double-stranded templates and low activity on single-stranded templates. In vivo, at least one additional N4-coded protein (p17) is required for N4 middle transcription. We show that N4 ORF2 encodes p17 (gp2). Characterization of purified gp2revealed that it is a single-stranded DNA-binding protein that activates N4 RNAPII transcription on single-stranded DNA templates through specific interaction with N4 RNAPII. On the basis of the properties of the proteins involved in N4 RNAPII transcription and of middle promoters, we propose a model for N4 RNAPII promoter recognition, in which gp2plays two roles, stabilization of a single-stranded region at the promoter and recruitment of N4 RNAPII through gp2-N4 RNAPII interactions. Furthermore, we discuss our results in the context of transcription initiation by mitochondrial RNA polymerases.

AraC protein: a love-hate relationship

In the bacterium Escherichia coli, the AraC protein positively and negatively regulates expression of the proteins required for the uptake and catabolism of the sugar L-arabinose. This essay describes how work from my laboratory on this system spanning more than thirty years has aided our understanding of positive regulation, revealed DNA looping (a mechanism that explains many action-at-a-distance phenomena) and, more recently, has uncovered the mechanism by which arabinose shifts AraC from a state where it prefers to bind to two well-separated DNA half-sites and form a DNA loop to a state where it binds to two adjacent half-sites and activates transcription. This work required learning how to assay, purify, and work with a protein possessing highly uncooperative biochemical properties. Present work is focussed on understanding arabinose-responsive mechanism in atomic detail and is also directed towards understanding protein structure and function well enough to be able to engineer the allosteric mechanism seen in AraC onto other proteins.

Inducible bacteriocin production in Lactobacillus is regulated by differential expression of the pln operons and by two antagonizing response regulators, the activity of which is enhanced upon phosphorylation

Expression of the five (pln) operons involved in the bacteriocin production of Lactobacillus plantarum C11 is regulated by a so-called pheromone-based signal-transducing network, in which the peptide pheromone (PlnA) induces bacteriocin production through the action of a histidine protein kinase (PlnB) and two antagonizing response regulators (PlnC as an activator and PlnD as a negative regulator). All pln-regulated promoters contain a conserved pair of direct repeats that serve as binding sites for PlnC and PlnD. In the present work, we show that the five PlnA-responsive operons are differentially expressed with regard to both timing and strength, and that the pheromone triggers a strong autoactivating loop of the regulatory unit (plnABCD) during an early stage of induction that gradually leads to enhanced activation of the other operons. The transport operon (plnGHSTUV), which is involved in the secretion of the pheromone and bacteriocins, is also expressed relatively early upon induction, but is quickly turned off soon after peak expression. Further investigation of the various promoters revealed that, although subtle differences within the promoter regions could account for the observed differential regulation, the presence of a downstream promoter-proximal sequence in one promoter was found to cause delayed peak activity. How phosphorylation regulates the activity of the pln response regulators was also accessed by direct mutagenesis at their phosphorylation sites. It was found that the two response regulators exert activity at two different levels: a low level when they are not phosphorylated and an elevated level when they are phosphorylated. The present data demonstrate that bacteriocin production in L. plantarum C11 is a highly regulated process, in which different regulatory mechanisms are applied to fine tune the timing and strength of expression of the five pln operons.

Subdivision of the helix-turn-helix GntR family of bacterial regulators in the FadR, HutC, MocR, and YtrA subfamilies

Haydon and Guest (Haydon, D. J, and Guest, J. R. (1991) FEMS Microbiol. Lett. 63, 291-295) first described the helix-turn-helix GntR family of bacterial regulators. They presented them as transcription factors sharing a similar N-terminal DNA-binding (d-b) domain, but they observed near-maximal divergence in the C-terminal effector-binding and oligomerization (E-b/O) domain. To elucidate this C-terminal heterogeneity, structural, phylogenetic, and functional analyses were performed on a family that now comprises about 270 members. Our comparative study first focused on the C-terminal E-b/O domains and next on DNA-binding domains and palindromic operator sequences, has classified the GntR members into four subfamilies that we called FadR, HutC, MocR, and YtrA. Among these subfamilies a degree of similarity of about 55% was observed throughout the entire sequence. Structure/function associations were highlighted although they were not absolutely stringent. The consensus sequences deduced for the DNA-binding domain were slightly different for each subfamily, suggesting that fusion between the D-b and E-b/O domains have occurred separately, with each subfamily having its own D-b domain ancestor. Moreover, the compilation of the known or predicted palindromic cis-acting elements has highlighted different operator sequences according to our subfamily subdivision. The observed C-terminal E-b/O domain heterogeneity was therefore reflected on the DNA-binding domain and on the cis-acting elements, suggesting the existence of a tight link between the three regions involved in the regulating process.

Biophysical evidence of arm-domain interactions in AraC

We report development of a method for the direct measurement of the interaction between the N-terminal arm and the remainder of the dimerization domain in the Escherichia coli AraC protein, the regulator of the l-arabinose operon. The interaction was measured using surface plasmon resonance to monitor the association between the immobilized peptide arm and the dimerization domain, truncated of its arm, in solution. As expected from genetic and physiological data, the interaction is strongly stimulated by l-arabinose and is insensitive to sugars like d-glucose or d-galactose. Alterations in the sequence of the arm which physiological experiments predict either to strengthen or weaken the arm produce the expected responses.

Mapping arm-DNA-binding domain interactions in AraC

AraC protein, the regulator of the l-arabinose operon in Escherichia coli has been postulated to function by a light switch mechanism. According to this mechanism, it should be possible to find mutations in the DNA-binding domain of AraC that result in weaker arm-DNA-binding domain interactions and which make the protein constitutive, that is, it no longer requires arabinose to activate transcription. We isolated such mutations by randomizing three contiguous leucine residues in the DNA-binding domain, and then by systematically scanning surface residues of the DNA-binding domain with alanine and glutamic acid. As a result, a total of 20 constitutive mutations were found at ten different positions. They form a contiguous trail on the DNA-distal face of the DNA-binding domain, and likely define the region where the N-terminal arm that extends from the N-terminal dimerization domain contacts the C-terminal DNA-binding domain.

The AraC transcriptional activators

The AraC family of bacterial transcriptional activators regulate diverse genetic systems. Recent X-ray diffraction studies show that the monomeric MarA and Rob activators bind to their asymmetric degenerate DNA sites via two different helix-turn-helix elements. Activation by MarA, SoxS or Rob requires a particular orientation of the asymmetric binding sequence (and hence the activator), depending on its distance from the -10 RNAP signal. Genetic studies are beginning to clarify how the activators interact with RNAP. Growing evidence suggests that for the sugar metabolism activators, multiple binding sites upstream of the promoter anchor the activator in a repressing or nonactivating configuration. By interaction with the sugar and/or CRP, the activator is allosterically altered so it can bind a new set of sites that enable it to activate the promoter. Surprisingly, the virulence activator, Rns, must bind to both upstream and downstream sites in order to activate the rns promoter.