Thermodynamic analysis of the lactose repressor-operator DNA interaction

Bacterial chromatin proteins, transcription, and DNA topology: Inseparable partners in the control of gene expression

DNA in bacterial chromosomes is organized into higher-order structures by DNA-binding proteins called nucleoid-associated proteins (NAPs) or bacterial chromatin proteins (BCPs). BCPs often bind to or near DNA loci transcribed by RNA polymerase (RNAP) and can either increase or decrease gene expression. To understand the mechanisms by which BCPs alter transcription, one must consider both steric effects and the topological forces that arise when DNA deviates from its fully relaxed double-helical structure. Transcribing RNAP creates DNA negative (-) supercoils upstream and positive (+) supercoils downstream whenever RNAP and DNA are unable to rotate freely. This (-) and (+) supercoiling generates topological forces that resist forward translocation of DNA through RNAP unless the supercoiling is constrained by BCPs or relieved by topoisomerases. BCPs also may enhance topological stress and overall can either inhibit or aid transcription. Here, we review current understanding of how RNAP, BCPs, and DNA topology interplay to control gene expression.

Transcription of Bacterial Chromatin

Decades of research have probed the interplay between chromatin (genomic DNA associated with proteins and RNAs) and transcription by RNA polymerase (RNAP) in all domains of life. In bacteria, chromatin is compacted into a membrane-free region known as the nucleoid that changes shape and composition depending on the bacterial state. Transcription plays a key role in both shaping the nucleoid and organizing it into domains. At the same time, chromatin impacts transcription by at least five distinct mechanisms: (i) occlusion of RNAP binding; (ii) roadblocking RNAP progression; (iii) constraining DNA topology; (iv) RNA-mediated interactions; and (v) macromolecular demixing and heterogeneity, which may generate phase-separated condensates. These mechanisms are not mutually exclusive and, in combination, mediate gene regulation. Here, we review the current understanding of these mechanisms with a focus on gene silencing by H-NS, transcription coordination by HU, and potential phase separation by Dps. The myriad questions about transcription of bacterial chromatin are increasingly answerable due to methodological advances, enabling a needed paradigm shift in the field of bacterial transcription to focus on regulation of genes in their native state. We can anticipate answers that will define how bacterial chromatin helps coordinate and dynamically regulate gene expression in changing environments.

Price of disorder in the lac repressor hinge helix

The Lac system of genes has been pivotal in understanding gene regulation. When the lac repressor protein binds to the correct DNA sequence, the hinge region of the protein goes through a disorder to order transition. The structure of this region of the protein is well understood when it is in this bound conformation, but less so when it is not. Structural studies show that this region is flexible. Our simulations show this region is extremely flexible in solution; however, a high concentration of salt can help kinetically trap the hinge helix. Thermodynamically, disorder is more favorable without the DNA present.

The mechanism and high-free-energy transition state of lac repressor-lac operator interaction

Significant, otherwise-unavailable information about mechanisms and transition states (TS) of protein folding and binding is obtained from solute effects on rate constants. Here we characterize TS for lac repressor(R)–lac operator(O) binding by analyzing effects of RO-stabilizing and RO-destabilizing solutes on association (k_a) and dissociation (k_d) rate constants. RO-destabilizing solutes (urea, KCl) reduce k_a comparably (urea) or more than (KCl) they increase k_d, demonstrating that they destabilize TS relative to reactants and RO, and that TS exhibits most of the Coulombic interactions between R and O. Strikingly, three solutes which stabilize RO by favoring burial/dehydration of amide oxygens and anionic phosphate oxygens all reduce k_d without affecting k_a significantly. The lack of stabilization of TS by these solutes indicates that O phosphates remain hydrated in TS and that TS preferentially buries aromatic carbons and amide nitrogens while leaving amide oxygens exposed. In our proposed mechanism, DNA-binding-domains (DBD) of R insert in major grooves of O pre-TS, forming most Coulombic interactions of RO and burying aromatic carbons. Nucleation of hinge helices creates TS, burying sidechain amide nitrogens. Post-TS, hinge helices assemble and the DBD-hinge helix-O-DNA module docks on core repressor, partially dehydrating phosphate oxygens and tightening all interfaces to form RO.

Multiple LacI-mediated loops revealed by Bayesian statistics and tethered particle motion

The bacterial transcription factor LacI loops DNA by binding to two separate locations on the DNA simultaneously. Despite being one of the best-studied model systems for transcriptional regulation, the number and conformations of loop structures accessible to LacI remain unclear, though the importance of multiple coexisting loops has been implicated in interactions between LacI and other cellular regulators of gene expression. To probe this issue, we have developed a new analysis method for tethered particle motion, a versatile and commonly used in vitro single-molecule technique. Our method, vbTPM, performs variational Bayesian inference in hidden Markov models. It learns the number of distinct states (i.e. DNA–protein conformations) directly from tethered particle motion data with better resolution than existing methods, while easily correcting for common experimental artifacts. Studying short (roughly 100 bp) LacI-mediated loops, we provide evidence for three distinct loop structures, more than previously reported in single-molecule studies. Moreover, our results confirm that changes in LacI conformation and DNA-binding topology both contribute to the repertoire of LacI-mediated loops formed in vitro, and provide qualitatively new input for models of looping and transcriptional regulation. We expect vbTPM to be broadly useful for probing complex protein–nucleic acid interactions.

Precise quantification of transcription factors in a surface-based single-molecule assay

Biosensors have recognized a rapid development the last years in both industry and science. Recently, a single-molecule assay based on alternating laser excitation has been established for the quantitative detection of transcription factors. These proteins specifically recognize and bind DNA and play an important role in controlling gene expression. We implemented this assay format on a total internal reflection fluorescence microscope to detect transcription factors with immobilized single-molecule DNA biosensors. We quantify transcription factors via colocalization of the two halves of their binding site with immobilized single molecules of a two-color DNA biosensor. We could detect a model transcription factor, the bacterial lactose repressor, at different concentrations down to 150pM. We found that robust modeling of stoichiometry derived TIRF data is achieved with Student's t-distributions and nonlinear least-squares estimation with weights equal to the inverse of the expected number of bin entries. This significantly improved transcription factor concentration estimates with respect to distribution modeling with Gaussians without adding notable computational effort. The proposed model may enhance the precision of other single-molecule assays quantifying molecular distributions. Our measurements reliably confirm that the immobilized biosensor format is more sensitive than a previously published solution based approach.

Sequence dependence of transcription factor-mediated DNA looping

DNA is subject to large deformations in a wide range of biological processes. Two key examples illustrate how such deformations influence the readout of the genetic information: the sequestering of eukaryotic genes by nucleosomes and DNA looping in transcriptional regulation in both prokaryotes and eukaryotes. These kinds of regulatory problems are now becoming amenable to systematic quantitative dissection with a powerful dialogue between theory and experiment. Here, we use a single-molecule experiment in conjunction with a statistical mechanical model to test quantitative predictions for the behavior of DNA looping at short length scales and to determine how DNA sequence affects looping at these lengths. We calculate and measure how such looping depends upon four key biological parameters: the strength of the transcription factor binding sites, the concentration of the transcription factor, and the length and sequence of the DNA loop. Our studies lead to the surprising insight that sequences that are thought to be especially favorable for nucleosome formation because of high flexibility lead to no systematically detectable effect of sequence on looping, and begin to provide a picture of the distinctions between the short length scale mechanics of nucleosome formation and looping.

In vivo tests of thermodynamic models of transcription repressor function

One emphasis of the Gibbs Conference on Biothermodynamics is the value of thermodynamic measurements for understanding behaviors of biological systems. In this study, the correlation between thermodynamic measurements of in vitro DNA binding affinity with in vivo transcription repression was investigated for two transcription repressors. In the first system, which comprised an engineered LacI/GalR homolog, mutational changes altered the equilibrium constant for binding DNA. Changes correlated with altered repression, but estimates of in vivo repressor concentration suggest a ≥25-fold discrepancy with in vitro conditions. In the second system, changes in ligand binding to BirA altered dimerization and subsequent DNA occupancy. Again, these changes correlate with altered in vivo repression, but comparison with in vitro measurements reveals a ~10-fold discrepancy. Further analysis of each system suggests that the observed discrepancies between in vitro and in vivo results reflect the contributions of additional equilibria to the transcription repression process.

Single-molecule approaches to probe the structure, kinetics, and thermodynamics of nucleoprotein complexes that regulate transcription

Single-molecule experimentation has contributed significantly to our understanding of the mechanics of nucleoprotein complexes that regulate epigenetic switches. In this minireview, we will discuss the application of the tethered-particle motion technique, magnetic tweezers, and atomic force microscopy to (i) directly visualize and thermodynamically characterize DNA loops induced by the lac, gal, and lambda repressors and (ii) understand the mechanistic role of DNA-supercoiling and DNA-bending cofactors in both prokaryotic and eukaryotic systems.

The salt dependence of the interferon regulatory factor 1 DNA binding domain binding to DNA reveals ions are localized around protein and DNA

The equilibrium dissociation constant of the DNA binding domain of interferon regulatory factor 1 (IRF1 DBD) for its DNA binding site depends strongly on salt concentration and salt type. These dependencies are consistent with IRF1 DBD binding to DNA, resulting in the release of cations from the DNA and both release of anions from the protein and uptake of a cation by the protein. We demonstrated this by utilizing the fact that the release of fluoride from protein upon complex formation does not contribute to the salt concentration dependence of binding and by studying mutants in which charged residues in IRF1 DBD that form salt bridges with DNA phosphates are changed to alanine. The salt concentration dependencies of the dissociation constants of wild-type IRF1 DBD and the mutants R64A, D73A, K75A, and D73A/K75A were measured in buffer containing NaF, NaCl, or NaBr. The salt concentration and type dependencies of the mutants relative to wild-type IRF1 DBD provide evidence of charge neutralization by solution ions for R64 and by a salt bridge between D73 and K75 in buffer containing chloride or bromide salts. These data also allowed us to determine the number, type, and localization of condensed ions around both IRF1 DBD and its DNA binding site.

Use of urea and glycine betaine to quantify coupled folding and probe the burial of DNA phosphates in lac repressor-lac operator binding

Thermodynamic analysis of urea-biopolymer interactions and effects of urea on folding of proteins and alpha-helical peptides shows that urea interacts primarily with polar amide surface. Urea is therefore predicted to be a quantitative probe of coupled folding, remodeling, and other large-scale changes in the amount of water-accessible polar amide surface in protein processes. A parallel analysis indicates that glycine betaine [N,N,N-trimethylglycine (GB)] can be used to detect burial or exposure of anionic (carboxylate, phosphate) biopolymer surface. To test these predictions, we have investigated the effects of these solutes (0-3 m) on the formation of 1:1 complexes between lac repressor (LacI) and its symmetric operator site (SymL) at a constant KCl molality. Urea reduces the binding constant K(TO) [initial slope dlnK(TO)/dm(urea) = -1.7 +/- 0.2], and GB increases K(TO) [initial slope dlnK(TO)/dm(GB) = 2.1 +/- 0.2]. For both solutes, this derivative decreases with an increase in solute concentration. Analysis of these initial slopes predicts that (1.5 +/- 0.3) x 10(3) A2 of polar amide surface and (4.5 +/- 1.0) x 10(2) A2 of anionic surface are buried in the association process. Analysis of published structural data, together with modeling of unfolded regions of free LacI as extended chains, indicates that 1.5 x 10(3) A2 of polar amide surface and 6.3 x 10(2) A2 of anionic surface are buried in complexation. Quantitative agreement between structural and thermodynamic results is obtained for amide surface (urea); for anionic surface (GB), the experimental value is approximately 70% of the structural value. For LacI-SymL binding, two-thirds of the structurally predicted change in amide surface (1.0 x 10(3) A2) occurs outside the protein-DNA interface in protein-protein interfaces formed by folding of the hinge helices and interactions of the DNA binding domain (DBD) with the core of the repressor. Since urea interacts principally with amide surface, it is particularly well-suited to detect and quantify the extent of coupled folding and other large-scale remodeling events in the steps of protein-nucleic acid interactions and other protein associations.

Engineering a root-specific, repressor-operator gene complex

Strong, tissue-specific and genetically regulated expression systems are essential tools in plant biotechnology. An expression system tool called a 'repressor-operator gene complex' (ROC) has diverse applications in plant biotechnology fields including phytoremediation, disease resistance, plant nutrition, food safety, and hybrid seed production. To test this concept, we assembled a root-specific ROC using a strategy that could be used to construct almost any gene expression pattern. When a modified E. coli lac repressor with a nuclear localization signal was expressed from a rubisco small subunit expression vector, S1pt::lacIn, LacIn protein was localized to the nuclei of leaf and stem cells, but not to root cells. A LacIn repressible Arabidopsis actin expression vector A2pot was assembled containing upstream bacterial lacO operator sequences, and it was tested for organ and tissue specificity using beta-glucuronidase (GUS) and mercuric ion reductase (merA) gene reporters. Strong GUS enzyme expression was restricted to root tissues of A2pot::GUS/S1pt::lacIn ROC plants, while GUS activity was high in all vegetative tissues of plants lacking the repressor. Repression of shoot GUS expression exceeded 99.9% with no evidence of root repression, among a large percentage of doubly transformed plants. Similarly, MerA was strongly expressed in the roots, but not the shoots of A2pot::merA/S1pt::lacIn plants, while MerA levels remained high in both shoots and roots of plants lacking repressor. Plants with MerA expression restricted to roots were approximately as tolerant to ionic mercury as plants constitutively expressing MerA in roots and shoots. The superiority of this ROC over the previously described root-specific tobacco RB7 promoter is demonstrated.

The DNA trackwalkers: principles of lesion search and recognition by DNA glycosylases

DNA glycosylases, the pivotal enzymes in base excision repair, are faced with the difficult task of recognizing their substrates in a large excess of unmodified DNA. We present here a kinetic analysis of DNA glycosylase substrate specificity, based on the probability of error. This novel approach to this subject explains many features of DNA surveillance and catalysis of lesion excision by DNA glycosylases. This approach also is applicable to the general issue of substrate specificity. We discuss determinants of substrate specificity in damaged DNA and in the enzyme, as well as methods by which these determinants can be identified.

Effect of the detergent Tween-20 on the DNA affinity chromatography of Gal4, C/EBPalpha, and lac repressor with observations on column regeneration

C/EBPalpha, Gal4, and lac repressor, representing three different transcription factor homology families, were expressed as fusion proteins and used to characterize the effects of column aging, Mg2+, the nonionic detergent Tween-20, column loading, and bovine serum albumin on DNA-affinity chromatography. When lac-repressor-beta-galactosidase fusion protein is loaded onto a new DNA-Sepharose column, less elutes from a new column than one that has been used two or more times. Higher amounts of lac repressor, the Green Fluorescent Protein fusions with CAAT enhancer binding protein (C/EBPalpha) and Gal4, elute from the columns when 0.1% Tween-20 is added to the mobile phase. The amount of improvement found depends upon the transcription factor studied and the amount of the protein loaded on the column; lac repressor and Gal4 are eluted in higher amounts over a large range of protein loads while C/EBP shows the greatest effect at low protein loads. This detergent effect is seen when either Sepharose or silica is used for the stationary phase. Including bovine serum albumin in the mobile phase gives a similar though lesser improvement to that observed with Tween-20. Mg2+ or EDTA in the mobile phase gave similar chromatography for C/EBP; since EDTA protects columns from DNases, its inclusion in the mobile phase is preferred. After extended use, the DNA affinity columns no longer bind transcription factors and this is not due to losses of DNA from the columns. Two simple methods (sodium dodecylsulfate and KSCN) were developed to regenerate such worn out columns.

Methods for transcription factor separation

Recent advances in the separation of transcription factors (TFs) are reviewed in this article. An overview of the transcription factor families and their structure is discussed and a computer analysis of their sequences reveals that while they do not differ from other proteins in molecular mass or isoelectric pH, they do differ from other proteins in the abundance of certain amino acids. The chromatographic and electrophoretic methods which have been successfully used for purification and analysis are discussed and recent advances in stationary and mobile phase composition is discussed.

Detection and prevention of protein aggregation before, during, and after purification

The use of proteins for in vitro studies or as therapeutic agents is frequently hampered by protein aggregation during expression, purification, storage, or transfer into requisite assay buffers. A large number of potential protein stabilizers are available, but determining which are appropriate can take days or weeks. We developed a solubility assay to determine the best cosolvent for a given protein that requires very little protein and only a few hours to complete. This technique separates native protein from soluble and insoluble aggregates by filtration and detects both forms of protein by SDS-PAGE or Western blotting. Multiple buffers can be simultaneously screened to determine conditions that enhance protein solubility. The behavior of a single protein in mixtures and crude lysates can be analyzed with this technique, allowing testing prior to and throughout protein purification. Aggregated proteins can also be assayed for conditions that will stabilize native protein, which can then be used to improve subsequent purifications. This solubility assay was tested using both prokaryotic and eukaryotic proteins that range in size from 17 to 150 kDa and include monomeric and multimeric proteins. From the results presented, this technique can be applied to a variety of proteins.

Role of hydration in the binding of lac repressor to DNA

The osmotic stress technique was used to measure changes in macromolecular hydration that accompany binding of wild-type Escherichia coli lactose (lac) repressor to its regulatory site (operator O1) in the lac promoter and its transfer from site O1 to nonspecific DNA. Binding at O1 is accompanied by the net release of 260 +/- 32 water molecules. If all are released from macromolecular surfaces, this result is consistent with a net reduction of solvent-accessible surface area of 2370 +/- 550 A. This area is only slightly smaller than the macromolecular interface calculated for a crystalline repressor dimer-O1 complex but is significantly smaller than that for the corresponding complex with the symmetrical optimized O(sym) operator. The transfer of repressor from site O1 to nonspecific DNA is accompanied by the net uptake of 93 +/- 10 water molecules. Together these results imply that formation of a nonspecific complex is accompanied by the net release of 165 +/- 43 water molecules. The enhanced stabilities of repressor-DNA complexes with increasing osmolality may contribute to the ability of Escherichia coli cells to tolerate dehydration and/or high external salt concentrations.

Fine-tuning function: correlation of hinge domain interactions with functional distinctions between LacI and PurR

LacI and PurR are highly homologous proteins. Their functional units are homodimers, with an N-terminal DNA binding domain that comprises the helix-turn-helix (HTH), N-linker, and hinge regions from both monomers. Hinge structural changes are known to occur upon DNA dissociation but are difficult to monitor experimentally. The initial steps of hinge unfolding were therefore examined using molecular dynamics simulations, utilizing a truncated, chimeric protein comprising the LacI HTH/N-linker and PurR hinge. A terminal Gly-Cys-Gly was added to allow "dimerization" through disulfide bond formation. Simulations indicate that differences in LacI and PurR hinge primary sequence affect the quaternary structure of the hinge x hinge' interface. However, these alternate hinge orientations would be sterically restricted by the core domain. These results prompted detailed comparison of recently available DNA-bound structures for LacI and truncated LacI(1-62) with the PurR structure. Examination revealed that different N-linker and hinge contacts to the core domain of the partner monomer (which binds effector molecule) affect the juxtapositions of the HTH, N-linker, and hinge regions in the DNA binding domain. In addition, the two full-length repressors exhibit significant differences in the interactions between the core and the C-linker connection to the DNA binding domain. Both linkers and the hinge have been implicated in the allosteric response of these repressors. Intriguingly, one functional difference between these two proteins is that they exhibit opposite allosteric response to effector. Simulations and observed structural distinctions are correlated with mutational analysis and sequence information from the LacI/GalR family to formulate a mechanism for fine-tuning individual repressor function.

Temperature dependence of DNA affinity chromatography of transcription factors

Oligonucleotides bound by the CAAT enhancer binding protein (C/EBP), the lactose repressor, and Gal4 were chemically coupled to cyanogen bromide-activated Sepharose and the temperature dependence of transcription factor chromatography was characterized. Each transcription factor was applied to the appropriate column and eluted using a salt gradient at several temperatures. Each transcription factor showed a unique behavior. As temperature was increases, less salt was required to elute C/EBP, more salt was required to elute lac repressor, while Gal4 showed a biphasic dependency with the amount of salt first decreasing between 4 and 19 degrees C and then increasing above 19 degrees C. This temperature dependence is not due to protein or DNA unfolding but rather is a property of complex formation. By loading a column, washing it at a permissive temperature, and then rapidly changing the column temperature, highly selective elution can be obtained. The thermodynamics of this temperature effect are different for the binding of specific and nonspecific DNA sequences, making chromatography at different temperatures a potentially important way of purifying transcription factors.

Wrapping of flanking non-operator DNA in lac repressor-operator complexes: implications for DNA looping

In our studies of lac repressor tetramer (T)-lac operator (O) interactions, we observed that the presence of extended regions of non-operator DNA flanking a single lac operator sequence embedded in plasmid DNA produced large and unusual cooperative and anticooperative effects on binding constants (Kobs) and their salt concentration dependences for the formation of 1:1 (TO) and especially 1:2 (TO2) complexes. To explore the origin of this striking behavior we report and analyze binding data on 1:1 (TO) and 1:2 (TO2) complexes between repressor and a single O(sym) operator embedded in 40 bp, 101 bp, and 2514 bp DNA, over very wide ranges of [salt]. We find large interrelated effects of flanking DNA length and [salt] on binding constants (K(TO)obs, K(TO2)obs) and on their [salt]-derivatives, and quantify these effects in terms of the free energy contributions of two wrapping modes, designated local and global. Both local and global wrapping of flanking DNA occur to an increasing extent as [salt] decreases. Global wrapping of plasmid-length DNA is extraordinarily dependent on [salt]. We propose that global wrapping is driven at low salt concentration by the polyelectrolyte effect, and involves a very large number (>/similar 20) of coulombic interactions between DNA phosphates and positively charged groups on lac repressor. Coulombic interactions in the global wrap must involve both the core and the second DNA-binding domain of lac repressor, and result in a complex which is looped by DNA wrapping. The non-coulombic contribution to the free energy of global wrapping is highly unfavorable ( approximately +30-50 kcal mol(-1)), which presumably results from a significant extent of DNA distortion and/or entropic constraints. We propose a structural model for global wrapping, and consider its implications for looping of intervening non-operator DNA in forming a complex between a tetrameric repressor (LacI) and one multi-operator DNA molecule in vivo and in vitro. The existence of DNA wrapping in LacI-DNA interactions motivates the proposal that most if not all DNA binding proteins may have evolved the capability to wrap and thereby organize flanking regions of DNA.

Lactose repressor protein: functional properties and structure

The lactose repressor protein (LacI), the prototype for genetic regulatory proteins, controls expression of lactose metabolic genes by binding to its cognate operator sequences in E. coli DNA. Inducer binding elicits a conformational change that diminishes affinity for operator sequences with no effect on nonspecific binding. The release of operator is followed by synthesis of mRNA encoding the enzymes for lactose utilization. Genetic, chemical and physical studies provided detailed insight into the function of this protein prior to the recent completion of X-ray crystallographic structures. The structural information can now be correlated with the phenotypic data for numerous mutants. These structures also provide the opportunity for physical and chemical studies on mutants designed to examine various aspects of lac repressor structure and function. In addition to providing insight into protein structure-function correlations, LacI has been utilized in a wide variety of applications both in prokaryotic gene expression and in eukaryotic gene regulation and studies of mutagenesis.

Comparison of the DNA association kinetics of the Lac repressor tetramer, its dimeric mutant LacIadi, and the native dimeric Gal repressor

The rates of association of the tetrameric Lac repressor (LacI), dimeric LacIadi (a deletion mutant of LacI), and the native dimeric Gal repressor (GalR) to DNA restriction fragments containing a single specific site were investigated using a quench-flow DNase I "footprinting" technique. The dimeric proteins, LacIadi and GalR, and tetrameric LacI possess one and two DNA binding sites, respectively. The nanomolar protein concentrations used in these studies ensured that the state of oligomerization of each protein was predominantly either dimeric or tetrameric, respectively. The bimolecular association rate constants (ka) determined for the LacI tetramer exceed those of the dimeric proteins. The values of ka obtained for LacI, LacIadi, and GalR display different dependences on [KCl]. For LacIadi and GalR, they diminish as [KCl] increases from 25 mM to 200 mM, approaching rates predicted for three-dimensional diffusion. In contrast, the ka values determined for the tetrameric LacI remain constant up to 300 mM [KCl], the highest salt concentration that could be investigated by quench-flow footprinting. The enhanced rate of association of the tetramer relative to the dimeric proteins can be modeled by enhanced "sliding" (Berg, O. G., Winter, R. B., and von Hippel, P. H. (1981) Biochemistry 20, 6929-6948) of the LacI tetramer relative to the LacIadi dimer or a combination of enhanced sliding and the superimposition of "direct transfer" mediated by the bidentate DNA interactions of the tetramer.

Thermodynamics of sequence-specific protein-DNA interactions

The molecular recognition processes in sequence-specific protein-DNA interactions are complex. The only feature common to all sequence-specific protein-DNA structures is a large interaction interface, which displays a high degree of complementarity in terms of shape, polarity and electrostatics. Many molecular mechanisms act in concert to form the specific interface. These include conformational changes in DNA and protein, dehydration of surfaces, reorganization of ion atmospheres, and changes in dynamics. Here we review the current understanding of how different mechanisms contribute to the thermodynamics of the binding equilibrium and the stabilizing effect of the different types of noncovalent interactions found in protein-DNA complexes. The relation to the thermodynamics of small molecule-DNA binding and protein folding is also briefly discussed.

Modified nucleotides reveal the indirect role of the central base pairs in stabilizing the lac repressor-operator complex

Guanine residues in the lac operator were replaced by 2-aminopurine or purine analogues, pairing the modified nucleotides with C. The observed equilibrium dissociation constants for lac repressor binding to substituted operators were measured in 10 mM Tris, 150 mM KCl, 0.1 mM EDTA, 0.1 mM DTE, pH 7.6 at 25 degrees C. These measurements revealed five positions that destabilized the complex when substituted with either analogue. Two positions, which are related by a 2-fold symmetry, are in the major groove of the operator thought to directly interact with the protein. Three sites were in the central region of the operator. A purine analogue at a sixth site perturbed the local DNA structure and destabilized the complex. Alkylation interference experiments of the 2-aminopurine substituted operators demonstrated that, of the five affected, two substitutions displayed altered phosphate interference patterns at the phosphate adjacent to the substituted base. For these operators, complex formation was measured in different concentrations of KCl to assess the contribution of counterion release to the bimolecular process. The results indicated that both complexes were similar to wild-type, although minor changes were observed. The Kobs of the complex was then measured when 2-aminopurine or purine analogues were paired with uracil nucleotide, a base pair that serves to stabilize the DNA. The introduction of the new base pairs revealed two effects on the bimolecular interaction. For those operator sites that are thought to perturb the interaction directly, the affinity of the complex was weakened to levels observed for the singly-substituted operators. In contrast, the nucleotides of 2-aminopurine paired with uracil positioned in the central region of the operator served to enhance the stability of the complex. The purine-uracil base pair substitution on the other hand had a significant destabilizing effect on the interaction. We propose that the central base pairs modulate binding of the complex by altering the intrinsic properties of the DNA. Two specific attributes are required to achieve the lowest free energy of interaction. The DNA must have two interstrand hydrogen bonds to stabilize the duplex and it must have properties associated with directional bending or unwinding. This analysis does not rule out contributions by direct interactions between the protein and the central region of the operator but underscores how indirect effects play a major role in complex formation in this system.

Analysis of the thermodynamic linkage of DNA binding and ion binding for dimeric and tetrameric forms of the lac repressor

The salt concentration dependences of the observed association constants (Kobs) for the binding of wild-type lac repressor tetramer and the dimeric lacI-18 mutant repressor to lactose operator DNA were compared. For both proteins, the data are consistent with a class of linkage models in which ion binding by the protein is driven by differences in the ionic concentrations in bulk solution and in the volume near the DNA surface. The models that best agree with the data are those in which ion-binding reactions are cooperative. In spite of close agreement between these models and the data, the determination of ion stoichiometries and apparent ion-binding affinities requires additional mechanistic or structural information. The simplest ion-binding mechanism consistent with the data is compatible with a current structural model of the repressor-operator complex. At salt concentrations in excess of 50 mM, at which cation displacement from the DNA and anion displacement from the protein are expected to dominate, similar ion stoichiometries are found for the DNA binding of dimeric and tetrameric repressors. This supports the notion that the DNA contacts of these proteins are homologous. At lower salt concentrations, in which cation binding by the proteins is expected to be significant, the net ion stoichiometry of wild-type repressor binding is slightly greater than that of the lacI-18 mutant. This difference may reflect the availability of ion-binding sites in the distal subunits of tetramer that are not present in the dimer, or may be a consequence of the involvement of ion binding in the dimer/monomer equilibrium.

Coupling of local folding to site-specific binding of proteins to DNA

Thermodynamic studies have demonstrated the central importance of a large negative heat capacity change (delta C degree assoc) in site-specific protein-DNA recognition. Dissection of the large negative delta C degree assoc and the entropy change of protein-ligand and protein-DNA complexation provide a thermodynamic signature identifying processes in which local folding is coupled to binding. Estimates of the number of residues that fold on binding obtained from this analysis agree with structural data. Structural comparisons indicate that these local folding transitions create key parts of the protein-DNA interface. The energetic implications of this "induced fit" model for DNA site recognition are considered.

Interlocking of plasmid DNAs due to Lac repressor-operator interaction

The presence of a single lac repressor binding sequence on plasmid DNAs is shown to mediate the formation of interlocked dimers in E. coli. The presence of both homo- and hetero-interlocked dimers suggests that the lac repressor complex can bring together randomly two plasmid DNA molecules to facilitate gyrase-mediated interlocking. The exclusive formation of multiply intertwined dimers also suggest that the lac repressor complex may bind simultaneously to a pair of replicated daughter plasmid molecules prior to their segregation. The formation of interlocked plasmid DNAs can be indicative of interaction between two DNA bound proteins in vivo.

Thermodynamic stoichiometries of participation of water, cations and anions in specific and non-specific binding of lac repressor to DNA. Possible thermodynamic origins of the "glutamate effect" on protein-DNA interactions

The objective of this study is to quantify the contributions of cations, anions and water to stability and specificity of the interaction of lac repressor (lac R) protein with the strong-binding symmetric lac operator (Osym) DNA site. To this end, binding constants Kobs and their power dependences on univalent salt (MX) concentration (SKobs = d log Kobs/d log[MX]) have been determined for the interactions of lac R with Osym operator and with non-operator DNA using filter binding and DNA cellulose chromatography, respectively. For both specific and non-specific binding of lac R, Kobs at fixed salt concentration [KX] increases when chloride (Cl-) is replaced by the physiological anion glutamate (Glu-). At 0.25 M-KX, the increase in Kobs for Osym is observed to be approximately 40-fold, whereas for non-operator DNA the increase in Kobs is estimated by extrapolation to be approximately 300-fold. For non-operator DNA, SKobsRD is independent of salt concentration within experimental uncertainty, and is similar in KCl (SKobs,RDKCl = -9.8(+/- 1.0) between 0.13 M and 0.18 M-KCl) and KGlu (SKobs,RDKGlu = -9.3(+/- 0.7) between 0.23 M and 0.36 M-KGlu). For Osym DNA, SKobsRO varies significantly with the nature of the anion, and, at least in KGlu appears to decrease in magnitude with increasing [KGlu]. Average magnitudes of SKobsRO are less than SKobsRD, and, for specific binding decrease in the order [SKobsRO,KCl[>[SKobsRO,KAc[>[SKobsRO,KGlu[ . Neither KobsRO nor SKobsRO is affected by the choice of univalent cation M+ (Na+, K+, NH4+, or mixtures thereof, all as the chloride salt), and SKobsRO is independent of [MCl] in the range examined (0.125 to 0.3 M). This behavior of SKobsRO is consistent with that expected for a binding process with a large contribution from the polyelectrolyte effect. However, the lack of an effect of the nature of the cation on the magnitude of KobsRO at a fixed [MX] is somewhat unexpected, in view of the order of preference of cations for the immediate vicinity of DNA (NH4+ > K+ > Na+) observed by 23Na nuclear magnetic resonance. For both specific and non-specific binding, the large stoichiometry of cation release from the DNA polyelectrolyte is the dominant contribution to SKobs. To interpret these data, we propose that Glu- is an inert anion, whereas Ac- and Cl- compete with DNA phosphate groups in binding to lac repressor. A thermodynamic estimate of the minimum stoichiometry of water release from lac repressor and Osym operator (210(+/- 30) H2O) is determined from analysis of the apparently significant reduction in [SKobsRO,KGlu[ with increasing [KGlu] in the range 0.25 to 0.9 M. According to this analysis, SKobs values of specific and non-specific binding in KGlu differ primarily because of the release of water in specific binding. In KAc and KCl, we deduce that anion competition affects Kobs and SKobs to an extent which differs for different anions and for the different binding modes.

Thermodynamics of Cro protein-DNA interactions

Using a highly sensitive pulsed-flow microcalorimeter, we have measured the changes in enthalpy and determined the thermodynamic parameters delta H, delta S degree, delta G degree, and delta C(p) for Cro protein-DNA association reactions. The reactions studied include sequence-nonspecific DNA association and sequence-specific DNA associations involving single- and multiple-base alterations and/or single-amino acid alteration mutants. (i) The association of Cro protein with nonspecific DNA at 15 degrees C is characterized by delta H = +4.4 kcal.mol-1 (1 cal = 4.18J), delta S degrees = 49 cal.mol-1.K-1, delta G degrees = -9.7 kcal.mol-1, and delta Cp congruent to 0; the association with specific high-affinity operator OR3 DNA is characterized by delta H = +0.8 kcal.mol-1, delta S degree = 59 cal.mol-1.K-1, delta G degree = -16.1 kcal.mol-1, and delta Cp = -360 cal.mol-1.K-1, respectively. Both nonspecific and specific Cro-DNA associations are entropy-driven. (ii) Plots of delta H vs. delta Cp and delta S degree vs. delta Cp for the 20 association reactions studied fall into two correlation groups with linear slopes of +9.4 K and -20.5 K and of -0.03 and -0.14, respectively. These regression lines have common intercepts, at the delta H and delta S degree values of nonspecific association (where delta Cp congruent to 0). The results suggest that there are, at least, two distinct conformational subclasses in specific Cro-DNA complexes, stabilized by different combinations of enthalpic and entropic contributions. The delta G degree and delta Cp values form an approximately single linear correlation group as a consequence of compensatory contributions from delta H and delta S degree to delta G degree and to delta Cp. Cro protein-DNA associations share some similar thermodynamic properties with protein folding, but their overall energetics are quite different. Although the nonspecific complex is stabilized predominantly by electrostatic forces, it appears that H bonds, van der Waals contacts, hydrophobic effects, and charge interactions all contribute to the stability (delta G degree and delta Cp) of the specific complex. (iii) The variations in the values of the thermodynamic parameters are in general accord with our knowledge of the structure of the Cro-DNA complex.

DNA looping

DNA-looping mechanisms are part of networks that regulate all aspects of DNA metabolism, including transcription, replication, and recombination. DNA looping is involved in regulation of transcriptional initiation in prokaryotic operons, including ara, gal, lac, and deo, and in phage systems. Similarly, in eukaryotic organisms, the effects of enhancers appear to be mediated at least in part by loop formation, and examples of DNA looping by hormone receptor proteins and developmental regulatory proteins have been found. In addition, instances of looped structures have been found in replication and in recombination in both prokaryotes and eukaryotes. DNA loop formation may have different functions in different cellular contexts; in some cases, the loop itself is requisite for regulation, while in others the increase in the effective local concentration of protein may account for the effects observed. The ability of DNA to form loops is affected by the distance between binding sites; by the DNA sequence, which determines deformability and bendability; and by the presence of other proteins that exert an influence on the conformation of a particular sequence. Alteration of the stability of DNA loops and/or protein-DNA binding by extra- or intracellular signals provides responsivity to changing metabolic or environmental conditions. The fundamental property of site-specific protein binding to DNA can be combined with protein-protein and protein-ligand interaction to generate a broad range of physiological states.

31P NMR spectra of oligodeoxyribonucleotide duplex lac operator-repressor headpiece complexes: importance of phosphate ester backbone flexibility in protein-DNA recognition

The 31P NMR spectra of various 14-base-pair lac operators bound to both wild-type and mutant lac repressor headpiece proteins were analyzed to provide information on the backbone conformation in the complexes. The 31P NMR spectrum of a wild-type symmetrical operator, d(TGTGAGCGCTCACA)2, bound to the N-terminal 56-residue headpiece fragment of a Y7I mutant repressor was nearly identical to the spectrum of the same operator bound to the wild-type repressor headpiece. In contrast, the 31P NMR spectrum of the mutant operator, d(TATAGAGCGCTCATA)2, wild-type headpiece complex was significantly perturbed relative to the wild-type repressor-operator complex. The 31P chemical shifts of the phosphates of a second mutant operator, d(TGTGTGCGCACACA)2, showed small but specific changes upon complexation with either the wild-type or mutant headpiece. The 31P chemical shifts of the phosphates of a third mutant operator, d(TCTGAGCGCTCAGA)2, showed no perturbations upon addition of the wild-type headpiece. The 31P NMR results provide further evidence for predominant recognition of the 5'-strand of the 5'-TGTGA/3'-ACACT binding site in a 2:1 protein to headpiece complex. It is proposed that specific, strong-binding operator-protein complexes retain the inherent phosphate ester conformational flexibility of the operator itself, whereas the phosphate esters are conformationally restricted in the weak-binding operator-protein complexes. This retention of backbone torsional freedom in strong complexes is entropically favorable and provides a new (and speculative) mechanism for protein discrimination of different operator binding sites. It demonstrates the potential importance of phosphate geometry and flexibility on protein recognition and binding.

Analysis of equilibrium and kinetic measurements to determine thermodynamic origins of stability and specificity and mechanism of formation of site-specific complexes between proteins and helical DNA

The concentration and nature of the electrolyte are key factors determining (1) the equilibrium extent of binding of oligocations or proteins to DNA, (2) the distribution of bound protein between specific and nonspecific sites, and (3) the kinetics of association and dissociation of both specific and nonspecific complexes. Salt concentration may therefore be used to great advantage to probe the thermodynamic basis of stability and specificity of protein-DNA complexes, and the mechanisms of association and dissociation. Cation concentration serves as a thermodynamic probe of the contributions to stability and specificity from neutralization of DNA phosphate charges and/or reduction in phosphate charge density. Cation concentration also serves as a mechanistic probe of the kinetically significant steps in association and dissociation that involve cation uptake. In general, effects of electrolyte concentration on equilibrium constants (quantified by SKobs) and rate constants (quantified by Skobs) are primarily cation effects that result from the cation-exchange character of the interactions of proteins and oligocations with polyanionic DNA. The competitive effects of Mg2+ or polyamines on the equilibria and kinetics of protein-DNA interactions are interpretable in the context of the cation-exchange model. The nature of the anion often has a major effect on the magnitude of the equilibrium constant (Kobs) and rate constant (kobs) of protein-DNA interactions, but a minor effect on SKobs and Skobs, which are dominated by the cation stoichiometry. The order of effects of different anions generally follows the Hofmeister series and presumably reflects the relative extent of preferential accumulation or exclusion of these anions from the relevant surface regions of DNA-binding proteins. The question of which anion is most inert (i.e., neither accumulated nor excluded from the relevant regions of these proteins) remains unanswered. The characteristic effects of temperature on equilibrium constants and rate constants for protein-DNA interactions also serve as diagnostic probes of the thermodynamic origins of stability and specificity and of the mechanism of the interaction, since large changes in thermodynamic and activation heat capacities accompany processes with large changes in the amount of water-accessible nonpolar surface area.(ABSTRACT TRUNCATED AT 400 WORDS)

Lac repressor-operator interaction: DNA length dependence

The interaction of the E. coli lac operon repressor with its operator DNA has been directly examined as a function of the length of operator-containing DNA. The apparent bimolecular association rate constants were calculated as ka = (kd/KD), where the dissociation equilibrium constant, KD and the dissociation rate constant, kd, were measured by nitrocellulose filter adsorption assays. The values obtained for the overall association rate constants are compared with theoretical association rate curves for specific mechanisms. Association of the repressor with short operator containing DNA fragments (less than 70 base pairs) occurs at rates expected of three-dimensional diffusion. Our data also imply that at longer DNA lengths a combination of three-dimensional diffusion with one-dimensional sliding along with hopping and/or intersegment transfer must be involved to facilitate the repressor operator association.

lac repressor forms stable loops in vitro with supercoiled wild-type lac DNA containing all three natural lac operators

We have analyzed protein-DNA complexes formed between lac repressor and linear or differently supercoiled lac DNA (802 or 816 base-pairs in length), which carry all three natural lac operators (O1, O2 and O3) in their wild-type sequence context and spacing and compared them with constructs that contain specifically mutated "pseudo-operators" O2 or O3. We used gel retardation assays to identify the nature of the complexes according to their characteristic electrophoretic mobility and dissociation rate measurements to determine their stability. With linear DNA we found only indirect evidence for loop formation between O1 and O2. In covalently closed DNA minicircles the formation of a loop between O1 and O2 could be demonstrated by the observation that O1-O2 containing DNA with low negative supercoiling (sigma = -0.013 and less) is constricted by binding of lac repressor, resulting in an increased electrophoretic mobility. At elevated negative supercoiling (sigma = -0.025, -0.037, -0.05) O1-O2 containing DNA complexed with lac repressor migrates significantly slower than the corresponding O1-DNA, indicating loop formation. The dissociation of lac repressor-operator complexes is decreased with increasing negative supercoiling for all tested operator combinations of O1, O2 and O3. However, in the presence of at least two natural lac operators on the same DNA minicircle the enhancement of stability is particularly large. This indicates that a DNA loop is formed between these two lac operators, O1 and O2 as well as O1 and O3, since negative supercoiling is known specifically to promote the formation of looped structures. Additionally, we observe a dependence of dissociation rate on the spatial alignment of the operators as a result of changing helical periodicity in differently supercoiled DNA and consider this to be further evidence for loop formation between O1 and O2 as well as O1 and O3.

Interaction of the Escherichia coli Gal repressor protein with its DNA operators in vitro

The binding of Escherichia coli Gal repressor to linear DNA fragments containing two binding sites (OE and OI) within the gal operon was analyzed in vitro with quantitative footprint and mobility-shift techniques. In vivo analysis of the regulation of the gal operon [Haber, R., & Adhya, S. (1988) Proc. Natl. Acad. Sci. U.S.A. 85, 9683-9687] has suggested the role of a regulatory "looped complex" mediated by the association of Gal repressor dimers bound at OE and OI. The binding of Gal repressor to a single site can be described by a model in which monomer and dimer are in equilibrium and only the dimer binds to DNA. At pH 7.0, 25 mM KCl, and 20 degrees C, the binding and dimerization free energies are comparable, suggesting that the equilibrium governing the formation of dimers may be important to regulation. The two intrinsic binding constants, delta GI and delta GE, and a constant describing cooperativity, delta GIE, were determined by footprint titration analysis as a function of pH, [KCl], and temperature. Only at 4 and 0 degrees C was delta GIE negative, signifying cooperative binding. These results are thought to be due to a weak dimer to tetramer association interface. delta GE and delta GI had maximal values between pH 6 and pH 7. The dependence of these constants on [KCl] corresponded to the displacement of approximately 2 ion equiv. The temperature dependence could be described by a change in the heat capacity, delta Cp, of -2.3 kcal mol-1 deg-1. Mobility-shift titration experiments conducted at 20 and 0 degrees C yielded values for delta GIE that were consistent with those resolved from the footprint analysis. Unique values of delta GIE were determined by analysis of mobility-shift titrations of Gal repressor with wild-type operator subject to the constraint that delta GE = delta GI: a procedure that eliminates the need to simultaneously analyze wild-type titrations with titrations of OE- and OI- operators.

Role of the hydrophobic effect in stability of site-specific protein-DNA complexes

The site-specific binding interaction of lac repressor with a symmetric operator sequence and of EcoRI endonuclease with its specific recognition site both exhibit a characteristic dependence of equilibrium binding constant (Kobs) on temperature, in which Kobs attains a relative maximum in the physiologically relevant temperature range. This behavior, which appears to be quite general for site-specific protein-DNA interactions, is indicative of a large negative standard heat capacity change (delta C0P,obs) in the association process. By analogy with model compound transfer studies and protein folding data, we propose that this delta C0P,obs results primarily from the removal of non-polar surface from water in the association process. From delta C0P,obs we obtain semiquantitative information regarding the change in water-exposed non-polar surface area (delta Anp) and the corresponding hydrophobic driving force for association (delta G0hyd): delta G0hyd approximately equal to 8(+/- 1) x 10(1) delta C0P,obs approximately equal to -22(+/- 5) delta Anp. We propose that removal of non-polar surface from water (the hydrophobic effect) and release of cations (the polyelectrolyte effect) drive the thermodynamically unfavorable process (e.g. conformational distortions) necessary to achieve mutually complementary recognition surfaces (at a steric and functional-group level) in the specific complex.

Association kinetics of site-specific protein-DNA interactions: roles of nonspecific DNA sites and of the molecular location of the specific site

We have applied the formalism developed previously for the kinetics of domain-localized reactions [S. Mazur and M. T. Record, Jr. (1986) Biopolymers 25, 985-1008] to describe complex mechanisms of association of a protein with a specific site on a large DNA molecule also containing many nonspecific binding sites. These nonspecific sites participate in the mechanism of formation of the specific complex through competitive binding and the facilitating mechanisms of sliding and transfer. The effects of localizing the sites in a domain are represented by a simple algebraic expression, and the sequence of interactions within the domain are described by equations closely related to a conventional, homogeneous solution mechanism. We apply this formalism to examine the interplay between sliding and direct transfer in domain-localized interactions in general and in the lac repressor-lac operator interaction in particular. Experimental investigation of the effect of the molecular location of the specific site (e.g., end vs middle of the polymer chain) on the kinetics of association may allow the contributions of sliding and direct transfer to be resolved.

Chemical modification of arginine residues in the lactose repressor

The lactose repressor protein was chemically modified with 2,3-butanedione and phenylglyoxal. Arginine reaction was quantitated by either amino acid analysis or incorporation of 14C-labeled phenylglyoxal. Inducer binding activity was unaffected by the modification of arginine residues, while both operator and nonspecific DNA binding activities were diminished, although to differing degrees. The correlation of the decrease in DNA binding activities with the modification of approximately 1-2 equiv of arginine per monomer suggests increased reactivity of a functionally essential residue(s). For both reagents, operator DNA binding activity was protected by the presence of calf thymus DNA, and the extent of reaction with phenylglyoxal was simultaneously diminished. This protection presumably results from steric restriction of reagent access to an arginine(s) that is (are) essential for DNA binding interactions. These experiments suggest that there is (are) an essential reactive arginine(s) critical for repressor binding to DNA.

Recombination by resolvase is inhibited by lac repressor simultaneously binding operators between res sites

The Tn3 resolvase requires that the two recombination (res) sites be aligned as direct repeats on the same molecule for efficient recombination to occur. To test whether resolvase must contact the DNA between res sites as predicted by tracking models, we have determined the sensitivity of recombination to protein diffusion blockades. Recombination between two res sites is unaffected either by lac repressor or bacteriophage T7 RNA polymerase being bound between them. Yet recombination is inhibited by lac repressor if the res site is bounded by a lac operator on both sides. We demonstrate that lac repressor will bind to more than one DNA site under the conditions used to assay recombination. This result suggests that lac repressor can inhibit resolvase by forming a DNA loop that isolates a res site topologically. These results do not support a tracking model for resolvase but suggest that the structure and topology of the DNA substrate is important in the formation of a synapse between res sites.