Cooperative and anticooperative effects in binding of the first and second plasmid Osym operators to a LacI tetramer: evidence for contributions of non-operator DNA binding by wrapping and looping

Biomolecular Assemblies: Moving from Observation to Predictive Design

Biomolecular assembly is a key driving force in nearly all life processes, providing structure, information storage, and communication within cells and at the whole organism level. These assembly processes rely on precise interactions between functional groups on nucleic acids, proteins, carbohydrates, and small molecules, and can be fine-tuned to span a range of time, length, and complexity scales. Recognizing the power of these motifs, researchers have sought to emulate and engineer biomolecular assemblies in the laboratory, with goals ranging from modulating cellular function to the creation of new polymeric materials. In most cases, engineering efforts are inspired or informed by understanding the structure and properties of naturally occurring assemblies, which has in turn fueled the development of predictive models that enable computational design of novel assemblies. This Review will focus on selected examples of protein assemblies, highlighting the story arc from initial discovery of an assembly, through initial engineering attempts, toward the ultimate goal of predictive design. The aim of this Review is to highlight areas where significant progress has been made, as well as to outline remaining challenges, as solving these challenges will be the key that unlocks the full power of biomolecules for advances in technology and medicine.

Proteins mediating DNA loops effectively block transcription

Loops are ubiquitous topological elements formed when proteins simultaneously bind to two noncontiguous DNA sites. While a loop‐mediating protein may regulate initiation at a promoter, the presence of the protein at the other site may be an obstacle for RNA polymerases (RNAP) transcribing a different gene. To test whether a DNA loop alters the extent to which a protein blocks transcription, the lac repressor (LacI) was used. The outcome of in vitro transcription along templates containing two LacI operators separated by 400 bp in the presence of LacI concentrations that produced both looped and unlooped molecules was visualized with scanning force microscopy (SFM). An analysis of transcription elongation complexes, moving for 60 s at an average of 10 nt/s on unlooped DNA templates, revealed that they more often surpassed LacI bound to the lower affinity O2 operator than to the highest affinity Os operator. However, this difference was abrogated in looped DNA molecules where LacI became a strong roadblock independently of the affinity of the operator. Recordings of transcription elongation complexes, using magnetic tweezers, confirmed that they halted for several minutes upon encountering a LacI bound to a single operator. The average pause lifetime is compatible with RNAP waiting for LacI dissociation, however, the LacI open conformation visualized in the SFM images also suggests that LacI could straddle RNAP to let it pass. Independently of the mechanism by which RNAP bypasses the LacI roadblock, the data indicate that an obstacle with looped topology more effectively interferes with transcription.

Aspects of protein-DNA interactions: a review of quantitative thermodynamic theory for modelling synthetic circuits utilising LacI and CI repressors, IPTG and the reporter gene lacZ

Protein-DNA interactions are central to the control of gene expression across all forms of life. The development of approaches to rigorously model such interactions has often been hindered both by a lack of quantitative binding data and by the difficulty in accounting for parameters relevant to the intracellular situation, such as DNA looping and thermodynamic non-ideality. Here, we review these considerations by developing a thermodynamically based mathematical model that attempts to simulate the functioning of an Escherichia coli expression system incorporating two of the best characterised prokaryotic DNA binding proteins, Lac repressor and lambda CI repressor. The key aim was to reproduce experimentally observed reporter gene activities arising from the expression of either wild-type CI repressor or one of three positive-control CI mutants. The model considers the role of several potentially important, but sometimes neglected, biochemical features, including DNA looping, macromolecular crowding and non-specific binding, and allowed us to obtain association constants for the binding of CI and its variants to a specific operator sequence.

Hypothetical biomolecular probe based on a genetic switch with tunable symmetry and stability

Background

Genetic switches are ubiquitous in nature, frequently associated with the control of cellular functions and developmental programs. In the realm of synthetic biology, it is of great interest to engineer genetic circuits that can change their mode of operation from monostable to bistable, or even to multistable, based on the experimental fine-tuning of readily accessible parameters. In order to successfully design robust, bistable synthetic circuits to be used as biomolecular probes, or understand modes of operation of such naturally occurring circuits, we must identify parameters that are key in determining their characteristics.

Results

Here, we analyze the bistability properties of a general, asymmetric genetic toggle switch based on a chemical-reaction kinetic description. By making appropriate approximations, we are able to reduce the system to two coupled differential equations. Their deterministic stability analysis and stochastic numerical simulations are in excellent agreement. Drawing upon this general framework, we develop a model of an experimentally realized asymmetric bistable genetic switch based on the LacI and TetR repressors. By varying the concentrations of two synthetic inducers, doxycycline and isopropyl β-D-1-thiogalactopyranoside, we predict that it will be possible to repeatedly fine-tune the mode of operation of this genetic switch from monostable to bistable, as well as the switching rates over many orders of magnitude, in an experimental setting. Furthermore, we find that the shape and size of the bistability region is closely connected with plasmid copy number.

Conclusions

Based on our numerical calculations of the LacI-TetR asymmetric bistable switch phase diagram, we propose a generic work-flow for developing and applying biomolecular probes: Their initial state of operation should be specified by controlling inducer concentrations, and dilution due to cellular division would turn the probes into memory devices in which information could be preserved over multiple generations. Additionally, insights from our analysis of the LacI-TetR system suggest that this particular system is readily available to be employed in this kind of probe.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0279-y) contains supplementary material, which is available to authorized users.

Simplified mechanistic models of gene regulation for analysis and design

Simplified mechanistic models of gene regulation are fundamental to systems biology and essential for synthetic biology. However, conventional simplified models typically have outputs that are not directly measurable and are based on assumptions that do not often hold under experimental conditions. To resolve these issues, we propose a ‘model reduction’ methodology and simplified kinetic models of total mRNA and total protein concentration, which link measurements, models and biochemical mechanisms. The proposed approach is based on assumptions that hold generally and include typical cases in systems and synthetic biology where conventional models do not hold. We use novel assumptions regarding the ‘speed of reactions’, which are required for the methodology to be consistent with experimental data. We also apply the methodology to propose simplified models of gene regulation in the presence of multiple protein binding sites, providing both biological insights and an illustration of the generality of the methodology. Lastly, we show that modelling total protein concentration allows us to address key questions on gene regulation, such as efficiency, burden, competition and modularity.

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.

FRET studies of a landscape of Lac repressor-mediated DNA loops

DNA looping mediated by the Lac repressor is an archetypal test case for modeling protein and DNA flexibility. Understanding looping is fundamental to quantitative descriptions of gene expression. Systematic analysis of LacI•DNA looping was carried out using a landscape of DNA constructs with lac operators bracketing an A-tract bend, produced by varying helical phasings between operators and the bend. Fluorophores positioned on either side of both operators allowed direct Förster resonance energy transfer (FRET) detection of parallel (P1) and antiparallel (A1, A2) DNA looping topologies anchored by V-shaped LacI. Combining fluorophore position variant landscapes allows calculation of the P1, A1 and A2 populations from FRET efficiencies and also reveals extended low-FRET loops proposed to form via LacI opening. The addition of isopropyl-β-d-thio-galactoside (IPTG) destabilizes but does not eliminate the loops, and IPTG does not redistribute loops among high-FRET topologies. In some cases, subsequent addition of excess LacI does not reduce FRET further, suggesting that IPTG stabilizes extended or other low-FRET loops. The data align well with rod mechanics models for the energetics of DNA looping topologies. At the peaks of the predicted energy landscape for V-shaped loops, the proposed extended loops are more stable and are observed instead, showing that future models must consider protein flexibility.

Anti-cooperative and cooperative protein-protein interactions between TetR isoforms on synthetic enhancers

Protein-protein interactions play an important role in determining the regulatory output of cis regulatory regions. In this work, we revisit the regulatory output functions recorded for the synthetic enhancers that contain binding sites for TetR. We use our thermodynamic model as an analysis tool to infer that two different types of interactions may take place between the TetR molecules. First, a strong mutually exclusive anti-cooperative interaction precludes the synthetic enhancer from being occupied by more than one AT (the aTc bound TetR isoform) protein, and a second weak cooperative interaction exists between the aTc-free TetR isoform (T). Consequently, this work highlights the power of the synthetic enhancer approach as a tool for studying protein-protein interactions via an experimentally verifiable prediction for the general mode of binding of the TetR repressor.

Crystal structure of a 1.6-hexanediol bound tetrameric form of Escherichia coli lac-repressor refined to 2.1 A resolution

We report the structure of a novel tetrameric form of the lactose repressor (LacI) protein from Escherichia coli refined to 2.1 A resolution. The tetramer is bound to 1.6-hexanediol present in the crystallization solution and the final R(free) for the structure is 0.201. The structure confirms previously reported structures on the monomer level. However, the tetramer is much more densely packed. This adds a new level of complexity to the interpretation of mutational effects and challenges details in the current model for LacI function. Several amino acids, previously associated with changes in function but unexplained at the structural level, appear in a new structural context in this tetramer which provides new implications for their function.

First-principles calculation of DNA looping in tethered particle experiments

We calculate the probability of DNA loop formation mediated by regulatory proteins such as Lac repressor (LacI), using a mathematical model of DNA elasticity. Our model is adapted to calculating quantities directly observable in tethered particle motion (TPM) experiments, and it accounts for all the entropic forces present in such experiments. Our model has no free parameters; it characterizes DNA elasticity using information obtained in other kinds of experiments. It assumes a harmonic elastic energy function (or wormlike chain type elasticity), but our Monte Carlo calculation scheme is flexible enough to accommodate arbitrary elastic energy functions. We show how to compute both the 'looping J factor' (or equivalently, the looping free energy) for various DNA construct geometries and LacI concentrations, as well as the detailed probability density function of bead excursions. We also show how to extract the same quantities from recent experimental data on TPM, and then compare to our model's predictions. In particular, we present a new method to correct observed data for finite camera shutter time and other experimental effects. Although the currently available experimental data give large uncertainties, our first-principles predictions for the looping free energy change are confirmed to within about 1 k(B)T, for loops of length around 300 basepairs. More significantly, our model successfully reproduces the detailed distributions of bead excursion, including their surprising three-peak structure, without any fit parameters and without invoking any alternative conformation of the LacI tetramer. Indeed, the model qualitatively reproduces the observed dependence of these distributions on tether length (e.g., phasing) and on LacI concentration (titration). However, for short DNA loops (around 95 basepairs) the experiments show more looping than is predicted by the harmonic-elasticity model, echoing other recent experimental results. Because the experiments we study are done in vitro, this anomalously high looping cannot be rationalized as resulting from the presence of DNA-bending proteins or other cellular machinery. We also show that it is unlikely to be the result of a hypothetical 'open' conformation of the LacI tetramer.

Concentration and length dependence of DNA looping in transcriptional regulation

In many cases, transcriptional regulation involves the binding of transcription factors at sites on the DNA that are not immediately adjacent to the promoter of interest. This action at a distance is often mediated by the formation of DNA loops: Binding at two or more sites on the DNA results in the formation of a loop, which can bring the transcription factor into the immediate neighborhood of the relevant promoter. These processes are important in settings ranging from the historic bacterial examples (bacterial metabolism and the lytic-lysogeny decision in bacteriophage), to the modern concept of gene regulation to regulatory processes central to pattern formation during development of multicellular organisms. Though there have been a variety of insights into the combinatorial aspects of transcriptional control, the mechanism of DNA looping as an agent of combinatorial control in both prokaryotes and eukaryotes remains unclear. We use single-molecule techniques to dissect DNA looping in the lac operon. In particular, we measure the propensity for DNA looping by the Lac repressor as a function of the concentration of repressor protein and as a function of the distance between repressor binding sites. As with earlier single-molecule studies, we find (at least) two distinct looped states and demonstrate that the presence of these two states depends both upon the concentration of repressor protein and the distance between the two repressor binding sites. We find that loops form even at interoperator spacings considerably shorter than the DNA persistence length, without the intervention of any other proteins to prebend the DNA. The concentration measurements also permit us to use a simple statistical mechanical model of DNA loop formation to determine the free energy of DNA looping, or equivalently, the for looping.

Network dynamics

Probably one of the most characteristic features of a living system is its continual propensity to change as it juggles the demands of survival with the need to replicate. Internally these adjustments are manifest as changes in metabolite, protein, and gene activities. Such changes have become increasingly obvious to experimentalists, with the advent of high-throughput technologies. In this chapter we highlight some of the quantitative approaches used to rationalize the study of cellular dynamics. The chapter focuses attention on the analysis of quantitative models based on differential equations using biochemical control theory. Basic pathway motifs are discussed, including straight chain, branched, and cyclic systems. In addition, some of the properties conferred by positive and negative feedback loops are discussed, particularly in relation to bistability and oscillatory dynamics.

Intrinsic curvature of DNA influences LacR-mediated looping

Protein-mediated DNA looping is a common mechanism for regulating gene expression. Loops occur when a protein binds to two operators on the same DNA molecule. The probability of looping is controlled, in part, by the basepair sequence of inter-operator DNA, which influences its structural properties. One structural property is the intrinsic or stress-free curvature. In this article, we explore the influence of sequence-dependent intrinsic curvature by exercising a computational rod model for the inter-operator DNA as applied to looping of the LacR-DNA complex. Starting with known sequences for the inter-operator DNA, we first compute the intrinsic curvature of the helical axis as input to the rod model. The crystal structure of the LacR (with bound operators) then defines the requisite boundary conditions needed for the dynamic rod model that predicts the energetics and topology of the intervening DNA loop. A major contribution of this model is its ability to predict a broad range of published experimental data for highly bent (designed) sequences. The model successfully predicts the loop topologies known from fluorescence resonance energy transfer measurements, the linking number distribution known from cyclization assays with the LacR-DNA complex, the relative loop stability known from competition assays, and the relative loop size known from gel mobility assays. In addition, the computations reveal that highly curved sequences tend to lower the energetic cost of loop formation, widen the energy distribution among stable and meta-stable looped states, and substantially alter loop topology. The inclusion of sequence-dependent intrinsic curvature also leads to nonuniform twist and necessitates consideration of eight distinct binding topologies from the known crystal structure of the LacR-DNA complex.

Analysis of kinetics in noisy systems: application to single molecule tethered particle motion

In the tethered particle motion method the length of a DNA molecule is monitored by measuring the range of diffusion of a microsphere tethered to the surface of a microscope coverslip through the DNA molecule itself. Looping of DNA (induced by binding of a specific protein) can be detected with this method and the kinetics of the looping/unlooping processes can be measured at the single molecule level. The microsphere's position variance represents the experimental variable reporting on the polymer length. Therefore, data windowing is required to obtain position variance from raw position data. Due to the characteristic diffusion time of the microsphere, the low-pass filtering required to attain a good signal/noise ratio (S/N) in the discrimination of looped versus unlooped state impacts significantly the measurement's time resolution. Here we present a method for measuring lifetimes based on half-amplitude thresholding and then correcting the kinetic measurements, taking into account low S/N (leading to false events) and limited time resolution (leading to missed events). This method allows an accurate and unbiased estimation of the kinetic parameters under investigation, independently of the choice of the window used for variance calculation, with potential applications to other single molecule measurements with low S/N.

Synthetic tetracycline-inducible regulatory networks: computer-aided design of dynamic phenotypes

BACKGROUND

Tightly regulated gene networks, precisely controlling the expression of protein molecules, have received considerable interest by the biomedical community due to their promising applications. Among the most well studied inducible transcription systems are the tetracycline regulatory expression systems based on the tetracycline resistance operon of Escherichia coli, Tet-Off (tTA) and Tet-On (rtTA). Despite their initial success and improved designs, limitations still persist, such as low inducer sensitivity. Instead of looking at these networks statically, and simply changing or mutating the promoter and operator regions with trial and error, a systematic investigation of the dynamic behavior of the network can result in rational design of regulatory gene expression systems. Sophisticated algorithms can accurately capture the dynamical behavior of gene networks. With computer aided design, we aim to improve the synthesis of regulatory networks and propose new designs that enable tighter control of expression.

RESULTS

In this paper we engineer novel networks by recombining existing genes or part of genes. We synthesize four novel regulatory networks based on the Tet-Off and Tet-On systems. We model all the known individual biomolecular interactions involved in transcription, translation, regulation and induction. With multiple time-scale stochastic-discrete and stochastic-continuous models we accurately capture the transient and steady state dynamics of these networks. Important biomolecular interactions are identified and the strength of the interactions engineered to satisfy design criteria. A set of clear design rules is developed and appropriate mutants of regulatory proteins and operator sites are proposed.

CONCLUSION

The complexity of biomolecular interactions is accurately captured through computer simulations. Computer simulations allow us to look into the molecular level, portray the dynamic behavior of gene regulatory networks and rationally engineer novel ones with useful applications. We are able to propose, test and accept or reject design principles for each network. Guided by simulations, we develop a set of design principles for novel tetracycline-inducible networks.

Optimization of a stochastically simulated gene network model via simulated annealing

By rearranging naturally occurring genetic components, gene networks can be created that display novel functions. When designing these networks, the kinetic parameters describing DNA/protein binding are of great importance, as these parameters strongly influence the behavior of the resulting gene network. This article presents an optimization method based on simulated annealing to locate combinations of kinetic parameters that produce a desired behavior in a genetic network. Since gene expression is an inherently stochastic process, the simulation component of simulated annealing optimization is conducted using an accurate multiscale simulation algorithm to calculate an ensemble of network trajectories at each iteration of the simulated annealing algorithm. Using the three-gene repressilator of Elowitz and Leibler as an example, we show that gene network optimizations can be conducted using a mechanistically realistic model integrated stochastically. The repressilator is optimized to give oscillations of an arbitrary specified period. These optimized designs may then provide a starting-point for the selection of genetic components needed to realize an in vivo system.

Lac repressor hinge flexibility and DNA looping: single molecule kinetics by tethered particle motion

The tethered particle motion (TPM) allows the direct detection of activity of a variety of biomolecules at the single molecule level. First pioneered for RNA polymerase, it has recently been applied also to other enzymes. In this work we employ TPM for a systematic investigation of the kinetics of DNA looping by wild-type Lac repressor (wt-LacI) and by hinge mutants Q60G and Q60 + 1. We implement a novel method for TPM data analysis to reliably measure the kinetics of loop formation and disruption and to quantify the effects of the protein hinge flexibility and of DNA loop strain on such kinetics. We demonstrate that the flexibility of the protein hinge has a profound effect on the lifetime of the looped state. Our measurements also show that the DNA bending energy plays a minor role on loop disruption kinetics, while a strong effect is seen on the kinetics of loop formation. These observations substantiate the growing number of theoretical studies aimed at characterizing the effects of DNA flexibility, tension and torsion on the kinetics of protein binding and dissociation, strengthening the idea that these mechanical factors in vivo may play an important role in the modulation of gene expression regulation.

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.

Long-range interactions between transcription factors

We discuss a method for analysing the number of GFP-LacI fusion transcription factors bound to a construct of 256 contiguous LacI binding sites using photon bleaching statics. We show by using a combination of imaging of the construct in nanochannels, photon statistics and addition of IGFP that the binding coefficient of the LacI decreases with increasing occupation of the construct, with a binding coefficient of 10(-6) M when only 15 of the 256 possible sites are occupied. We model this effect by assuming that the GFP-LacI dimer introduces elastic strain into the helix by generalized deformations, and that this strain propagates over distances at least as large as the persistence length.

Single-molecule studies of repressor-DNA interactions show long-range interactions

We have performed single-molecule studies of GFP-LacI repressor proteins bound to bacteriophage lambda DNA containing a 256 tandem lac operator insertion confined in nanochannels. An integrated photon molecular counting method was developed to determine the number of proteins bound to DNA. By using this method, we determined the saturated mean occupancy of the 256 tandem lac operators to be 13, which constitutes only 2.5% of the available sites. This low occupancy level suggests that the repressors influence each other even when they are widely separated, at distances on the order of 200 nm, or several DNA persistence lengths.

Coarse-grained model of entropic allostery

Many signaling functions in molecular biology require proteins to bind to substrates such as DNA in response to environmental signals such as the simultaneous binding to a small molecule. Examples are repressor proteins which may transmit information via a conformational change in response to the ligand binding. An alternative entropic mechanism of "allostery" suggests that the inducer ligand changes the intramolecular vibrational entropy, not just the mean static structure. We present a quantitative, coarse-grained model of entropic allostery, which suggests design rules for internal cohesive potentials in proteins employing this effect. It also addresses the issue of how the signal information to bind or unbind is transmitted through the protein. The model may be applicable to a wide range of repressors and also to signaling in trans-membrane proteins.

Thermal and urea-induced unfolding of the marginally stable lac repressor DNA-binding domain: a model system for analysis of solute effects on protein processes

Thermodynamic and structural evidence indicates that the DNA binding domains of lac repressor (lacI) exhibit significant conformational adaptability in operator binding, and that the marginally stable helix-turn-helix (HTH) recognition element is greatly stabilized by operator binding. Here we use circular dichroism at 222 nm to quantify the thermodynamics of the urea- and thermally induced unfolding of the marginally stable lacI HTH. Van't Hoff analysis of the two-state unfolding data, highly accurate because of the large transition breadth and experimental access to the temperature of maximum stability (T(S); 6-10 degrees C), yields standard-state thermodynamic functions (deltaG(o)(obs), deltaH(o)(obs), deltaS(o)(obs), deltaC(o)(P,obs)) over the temperature range 4-40 degrees C and urea concentration range 0 </= C(3) </= 6 M. For unfolding the HTH, deltaG(o)(obs) decreases linearly with increasing C(3) at all temperatures examined, which directly confirms the validity of the linear extrapolation method (LEM) to obtain the intrinsic stability of this protein. At 25 degrees C (pH 7.3 and 50 mM K(+)), both linear extrapolation and extrapolation via the local-bulk domain model (LBDM) to C(3) = 0 yield deltaG(o)(obs) = 1.23 +/- 0.05 kcal mol(-)(1), in agreement with direct measurement (1.24 +/- 0.30 kcal mol(-)(1)). Like deltaG(o)(obs), both deltaH(o)(obs) and deltaS(o)(obs) decrease linearly with increasing C(3); the derivatives with respect to C(3) of deltaG(o)(obs), deltaH(o)(obs) and TdeltaS(o)(obs) (in cal mol(-)(1) M(-)(1)) are -449 +/- 11, -661 +/- 90, and -203 +/- 91 at 25 degrees C, indicating that the effect of urea on deltaG(o)(obs) is primarily enthalpic. The deltaC(o)(P,obs) of unfolding (0.63 +/- 0.05 kcal mol(-)(1) K(-)(1)) is not detectibly dependent on C(3) or temperature. The urea m-value of the lacI HTH (-d deltaG(o)(obs),/dC(3) = 449 +/- 11 cal mol(-)(1) M(-)(1) at 25 degrees C) is independent of C(3) up to at least 6 M. Use of the LBDM to fit the C(3)-dependence of deltaG(o)(obs) yields the local-bulk partition coefficient for accumulation of urea at the protein surface exposed upon denaturation: K(P) = 1.103 +/- 0.002 at 25 degrees C. This partition coefficient is the same within uncertainty as those previously determined by LBDM analysis of osmometric data for solutions of urea and native (folded) bovine serum albumin, as well as LBDM analysis of the proportionality of m-values to changes in water accessible surface area upon protein unfolding. From the correspondence between values of K(P), we conclude that the average local urea concentration at both folded and unfolded protein surface exceeds the bulk by approximately 10% at 25 degrees C. The observed decrease in m-value for the lacI HTH with increasing temperature, together with the observed reductions in both deltaH(o)(obs) and deltaS(o)(obs) of unfolding with increasing urea concentration, demonstrate that K(P) for urea decreases with increasing temperature and that transfer of urea from the bulk solution to the local domain at the protein surface exposed on denaturation is enthalpically driven and entropically unfavorable.

Fluorescence resonance energy transfer over approximately 130 basepairs in hyperstable lac repressor-DNA loops

Lac repressor (LacI) binds two operator DNA sites, looping the intervening DNA. DNA molecules containing two lac operators bracketing a sequence-directed bend were previously shown to form hyperstable LacI-looped complexes. Biochemical studies suggested that orienting the operators outward relative to the bend direction (in construct 9C14) stabilizes a positively supercoiled closed form, with a V-shaped LacI, but that the most stable loop construct (11C12) is a more open form. Here, fluorescence resonance energy transfer (FRET) is measured on DNA loops, between fluorescein and TAMRA attached near the two operators, approximately 130 basepairs apart. For 9C14, efficient LacI-induced energy transfer ( approximately 74% based on donor quenching) confirms that the designed DNA shape can force the looped complex into a closed form. From enhanced acceptor emission, correcting for observed donor-dependent quenching of acceptor fluorescence, approximately 52% transfer was observed. Time-resolved FRET suggests that this complex exists in both closed- and open form populations. Less efficient transfer, approximately 10%, was detected for DNA-LacI sandwiches and 11C12-LacI, consistent with an open form loop. This demonstration of long-range FRET in large DNA loops confirms that appropriate DNA design can control loop geometry. LacI flexibility may allow it to maintain looping with other proteins bound or under different intracellular conditions.

Protein surface salt bridges and paths for DNA wrapping

The organization of large regions of DNA on the surface of proteins is critical to many DNA 'transactions', including replication, transcription, recombination and repair, as well as the packaging of chromosomal DNA. Recent thermodynamic and structural studies of DNA binding by integration host factor indicate that the disruption of protein surface salt bridges (dehydrated ion pairs) dominates the observed thermodynamics of integration host factor binding and, more generally, allows the wrapping of DNA on protein surfaces. The proposed thermodynamic signature of wrapping with coupled salt bridge disruption includes large negative salt-concentration-dependent enthalpy, entropy and heat capacity changes and smaller than expected magnitudes of the observed binding constant and its power dependence on salt concentration. Examination of the free structures of proteins recently shown to wrap DNA leads us to hypothesize that a pattern of surface salt bridges interspersed with cationic sidechains provides a structural signature for wrapping and that the number and organization of salt bridges and cationic groups dictate the thermodynamics and topology of DNA wrapping, which in turn are critical to function.

Ion concentration and temperature dependence of DNA binding: comparison of PurR and LacI repressor proteins

Purine repressor (PurR) binding to specific DNA is enhanced by complexing with purines, whereas lactose repressor (LacI) binding is diminished by interaction with inducer sugars despite 30% identity in their protein sequences and highly homologous tertiary structures. Nonetheless, in switching from low- to high-affinity DNA binding, these proteins undergo a similar structural change in which the hinge region connecting the DNA and effector binding domains folds into an alpha-helix and contacts the DNA minor groove. The differences in response to effector for these proteins should be manifest in the polyelectrolyte effect which arises from cations displaced from DNA by interaction with positively charged side chains on a protein and is quantitated by measurement of DNA binding affinity as a function of ion concentration. Consistent with structural data for these proteins, high-affinity operator DNA binding by the PurR-purine complex involved approximately 15 ion pairs, a value significantly greater than that for the corresponding state of LacI (approximately 6 ion pairs). For both proteins, however, conversion to the low-affinity state results in a decrease of approximately 2-fold in the number of cations released per dimeric DNA binding site. Heat capacity changes (DeltaC(p)) that accompany DNA binding, derived from buried apolar surface area, coupled folding, and restriction of motional freedom of polar groups in the interface, also reflect the differences between these homologous repressor proteins. DNA binding of the PurR-guanine complex is accompanied by a DeltaC(p) (-2.8 kcal mol(-1) K(-1)) more negative than that observed previously for LacI (-0.9 to -1.5 kcal mol(-1) K(-1)), suggesting that more extensive protein folding and/or enhanced structural rigidity may occur upon DNA binding for PurR compared to DNA binding for LacI. The differences between these proteins illustrate plasticity of function despite high-level sequence and structural homology and undermine efforts to predict protein behavior on the basis of such similarities.

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.

Designed hyperstable Lac repressor.DNA loop topologies suggest alternative loop geometries

Lac repressor (LacI) forms DNA loops which are critical for efficient operator binding and transcriptional repression. Designed DNA loops formed on three constructs with lac operators bracketing phased A-tract bends were characterized by mobility shift, footprinting, and DNA cyclization and topology. Operator dyad axes point either in or out relative to the sequence-induced curvature. Possible conformations suggested from X-ray structures of LacI and LacI.DNA include "wrapping away" (WA), "simple loop" (SL), and "wrapping toward" (WT) models. The WA loop should be preferentially stabilized by the outward operators, the SL and WT loops by the inward operators. Competition experiments demonstrated increased loop stability for all the bent constructs, with the SL/WT construct supporting hyperstable loops (t1/2 of days). This offers a general approach to stabilizing multi-protein DNA complexes on short DNA. DNA cyclization of loops gave minicircle products with altered topologies. WA constructs afforded relaxed and positive topoisomers, and cyclization kinetics indicated slow interconversion of precursors to the two topoisomers. The SL/WT construct gave a relaxed topoisomer, with a small amount of negative supercoil. These results suggest that while it is possible to force the WA loop to form (as in a model proposed from the LacI.DNA structure), the most stable loop geometry is different, probably a U-shape around an extended LacI tetramer. The topological results show how a protein-induced positive supercoil can be reconciled with LacI's preference for binding negatively supercoiled DNA. We suggest that looping proteins can affect the assembly of subsequent proteins by controlling loop topology.

Dimerisation mutants of Lac repressor. I. A monomeric mutant, L251A, that binds Lac operator DNA as a dimer

Dimer formation between monomers of the Escherichia coli Lac repressor is substantially specificed by the interactions between three alpha-helices in each monomer which form a hydrophobic interface. As a first step in analysing the specificity of this interaction, we examined the mutant L251A. LacR bearing this mutation in a background lacking the C-terminal heptad repeats is completely incapable of forming dimers in solution, with a dimer-monomer equilibrium dissociation constant, or Kd, higher than 10(-5)M. This correlates with a 200-fold decrease in its ability to repress the lac operon in vivo compared to dimeric LacR. Surprisingly, the mutant is still capable of forming dimers upon binding to short operator DNA in vitro. Analysis of the kinetic parameters of binding of the mutant to operator DNA reveals a 2000 to 3000-fold increase in the equilibrium dissociation constant (Kd) of the mutant-DNA complex in comparison to dimeric LacR-operator complexes, with the change almost entirely due to a greater than 1000-fold decrease in association rate. The dissociation rate varies only by a factor of about two, in comparison to dimeric LacR. This change reflects a kinetic pathway in which dimer formation, in solution or on DNA, is the rate-limiting step. These findings have implications for the specificity and stability of the protein-protein interface in question.

Thermodynamic analysis of unfolding and dissociation in lactose repressor protein

Lactose repressor protein, regulator of lac enzyme expression in Escherichia coli, maintains its structure and function at extremely low protein concentrations (<10(-)12 M). To examine the unfolding and dissociation of this tetrameric protein, structural transitions in the presence of varying concentrations of urea were monitored by fluorescence and circular dichroism spectroscopy, analytical ultracentrifugation, and functional activities. The spectroscopic data demonstrated a single cooperative transition with no evidence of folded dimeric or monomeric species of this protein. These spectroscopic transitions were reversible provided a long incubation step was employed in the refolding reaction at approximately 3 M urea. The refolded repressor protein possessed the same functional and structural properties as wild-type repressor protein. The absence of concentration dependence expected for tetramer dissociation to unfolded monomer (M4 <--> 4U) in the spectral transitions indicates that the disruption of the monomer-monomer interface and monomer unfolding are a concerted reaction (M4 <--> U4) that may occur prior to the dissociation of the dimer-dimer interface. Thus, we propose that the unfolded monomers remain associated at the C-terminus by the 4-helical coiled-coil structure that forms the dimer-dimer interface and that this intermediate is the end point detected in the spectral transitions. Efforts to confirm the existence of this species by ultracentrifugation were inhibited by the aggregation of this intermediate. Based upon these observations, the wild-type fluorescence and CD data were fit to a model, M4 <--> U4, which resulted in an overall DeltaG degrees for unfolding of 40 kcal/mol. Using a mutant protein, K84L, in which the monomer-monomer interface is stabilized, sedimentation equilibrium results demonstrated that the dimer-dimer interface of lac repressor could persist at higher levels of urea than the monomer-monomer interface. The tetramer-dimer transition monitored using this mutant repressor yields a DeltaG degrees of 20.4 kcal/mol. Using this free energy value for the dissociation process of U4 <--> 4U, an overall free energy change of approximately 60 kcal/mol was calculated for dissociation of all interfaces and unfolding of the tetrameric lac repressor, reflecting the exceptional stability of this protein.

The function of auxiliary operators

Gene regulation by control of transcription has been analysed in great detail both in prokaryotes and in eukaryotes. The frequency of transcription may be decreased by repressors or increased by activators. A repressor may work by decreasing the concentration of RNA polymerase at a promoter capable of forming an open complex. An activator may work by increasing the concentration of RNA polymerase at a promoter capable of forming an open complex. For this purpose, a strategy is used over and over again. It is called increase in local concentration. How Escherichia coli uses this strategy efficiently is discussed.

Operator search by mutant Lac repressors

The Escherichia coli Lac and Gal repressors are two members of a large family of bacterial repressor proteins that share significant sequence and structural homology. Efficient repression by all family members requires specific binding to a site or sites close to the transcriptional start of the genes regulated. Both LacR and GalR have to bind to at least two sites for efficient repression, yet they differ in one important respect: LacR is a homotetramer whereas GalR is a homodimer. In an attempt to understand this difference, we studied the operator binding activity of a LacR variant that has the DNA-binding specificity of GalR (LacR-V17A18). A tetrameric version of this protein shows a 30-fold decrease in association rate to operator located on a long (lambda) DNA molecule, in comparison to wild-type LacR, while a dimeric version of this protein shows an unaltered association rate in comparison to dimeric LacR. This reduction in association rate correlates with a broadened DNA-binding specificity for base-pairs 4 and 5 of the operator: examination of an additional LacR variant with an even broader DNA-binding specificity indicates that a tetrameric version also shows a 30-fold decrease in association rate in comparison to wild-type LacR, while a dimeric version again shows an unaltered association rate in comparison to dimeric LacR. This difference in association rate in vitro correlates with whether a tetrameric or dimeric variant of LacR of a given DNA-binding specificity will repress lacZ under control of a single operator more efficiently in vivo. We therefore propose that the formation of stable homotetramers becomes a distinct disadvantage unless a high degree of DNA-binding specificity is also present, and demonstrate that this in indeed the case for GalR-mediated repression of the gal operon. This functional constraint seems to have influenced the evolution of the LacI-GalR family of repressors, most of which have a relatively broad specificity of DNA-binding and most of which form only stable homodimers.

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 the interactions of lac repressor with variants of the symmetric lac operator: effects of converting a consensus site to a non-specific site

What are the thermodynamic consequences of the stepwise conversion of a highly specific (consensus) protein-DNA interface to one that is nonspecific? How do the magnitudes of key favorable contributions to complex stability (burial of hydrophobic surfaces and reduction of DNA phosphate charge density) change as the DNA sequence of the specific site is detuned? To address these questions we investigated the binding of lac repressor (LacI) to a series of 40 bp fragments carrying symmetric (consensus) and variant operator sequences over a range of temperatures and salt concentrations. Variant DNA sites contained symmetrical single and double base-pair substitutions at positions 4 and/or 5 [sequence: see text] in each 10 bp half site of the symmetric lac operator (Osym). Non-specific interactions were examined using a 40 bp non-operator DNA fragment. Disruption of the consensus interface by a single symmetrical substitution reduces the observed equilibrium association constant (K(obs)) for Osym by three to four orders of magnitude; double symmetrical substitutions approach the six orders in magnitude difference between specific and non-specific binding to a 40 bp fragment. At these adjacent positions in the consensus site, the free energy effects of multiple substitutions are non-additive: the first reduces /deltaG(obs)o/ by 3 to 5 kcal mol(-1), approximately halfway to the non-specific level, whereas the second is less deleterious, reducing /deltaG(obs)o/ by less than 3 kcal mol(-1). Variant-specific dependences of K(obs) on temperature and salt concentration characterize these LacI-operator interactions. In general, binding constants and standard free energies of binding both exhibit characteristic extrema near 290 K. As a consequence, both the enthalpic and entropic contributions to stability of Osym and variant complexes change from positive (i.e. entropy driven) at lower temperatures to negative (i.e. enthalpy driven) at higher temperatures, indicating that the heat capacity change upon binding, deltaC(obs)o, is large and negative. In general, /deltaC(obs)o/ decreases as the specificity and stability of the variant complex decreases. Stabilities of complexes of LacI with Osym and all variant operators are strongly [salt]-dependent. Binding constants for the variant complexes exhibit a power-dependence on [salt] that is larger in magnitude (i.e. more negative) than for Osym, but no obvious trend relates changes in contributions from the polyelectrolyte effect and the observed reductions in stability (delta deltaG(obs)o). These variant-specific thermodynamic signatures provide novel insights into the consequences of converting a consensus interface to a less specific one; such insights are not obtained from comparisons at the level of delta deltaG(obs)o. We propose that this variant-specific behavior arises from a strong effect of operator sequence on the extent of induced conformational changes in the protein (and possibly also in the DNA site) which accompany binding.

Lac repressor-operator complex

For many years the lac operon of Escherichia coli has been the paradigm for gene regulation. Recently, the structures of the lac repressor core bound to isopropyl-beta-D-1-thiogalactoside (IPTG), the intact apo lac repressor, the intact lac repressor complexes with IPTG and a 21-base-pair symmetric operator, and the refined headpiece of the repressor have been determined. These structures have provided a framework for understanding a wealth of biochemical and genetic information. An analysis of these structures, as well as a description of their function and a comparison to homologous proteins, is now possible.

Escherichia coli lac repressor-lac operator interaction and the influence of allosteric effectors

The wild type E. coli lac operator is embedded in a 35 base-pair DNA sequence containing extensive 2-fold symmetry, suggesting a symmetric repressor operator complex. However, deviations from strict 2-fold symmetry occur at the central base-pair and at three additional base-pairs. Using an operator fragment binding analysis we have determined: (a) a relative contribution each pair provides to the lac repressor-lac operator DNA complex, (b) the operator DNA length necessary for maximum binding to lac repressor; and (c) the contribution of the several non-symmetric base in the wild-type operator to the binding affinity. Since lac repressor-lac operator DNA interaction is reduced upon binding of the gratuitous inducer, isopropyl-beta-D-galactoside (IPTG), the same DNA fragment binding analysis was performed with the low affinity form of lac repressor. In the presence of inducer, the affinity for the left half site of the wild-type lac operator is reduced without significant reduction on the right half of the operator. Conversely, the anti-inducer orthonitrophenylfucoside (ONPF) which stabilizes the lac repressor-lac operator complex increases the binding affinity, particularly to the right half of the operator.