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

Designed architectural proteins that tune DNA looping in bacteria

Architectural proteins alter the shape of DNA. Some distort the double helix by introducing sharp kinks. This can serve to relieve strain in tightly-bent DNA structures. Here, we design and test artificial architectural proteins based on a sequence-specific Transcription Activator-like Effector (TALE) protein, either alone or fused to a eukaryotic high mobility group B (HMGB) DNA-bending domain. We hypothesized that TALE protein binding would stiffen DNA to bending and twisting, acting as an architectural protein that antagonizes the formation of small DNA loops. In contrast, fusion to an HMGB domain was hypothesized to generate a targeted DNA-bending architectural protein that facilitates DNA looping. We provide evidence from Escherichia coli Lac repressor gene regulatory loops supporting these hypotheses in living bacteria. Both data fitting to a thermodynamic DNA looping model and sophisticated molecular modeling support the interpretation of these results. We find that TALE protein binding inhibits looping by stiffening DNA to bending and twisting, while the Nhp6A domain enhances looping by bending DNA without introducing twisting flexibility. Our work illustrates artificial approaches to sculpt DNA geometry with functional consequences. Similar approaches may be applicable to tune the stability of small DNA loops in eukaryotes.

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.

Exploring LacI-DNA dynamics by multiscale simulations using the SIRAH force field

The lac repressor protein (LacI) together with its target regulatory sequence are a common model for studying DNA looping and its implications on transcriptional control in bacteria. Owing to the molecular size of this system, standard all-atom (AA) simulations are prohibitive for achieving relevant biological time scales. As an alternative, multiscale models, which combine AA descriptions at particular regions with coarse-grained (CG) representations of the remaining components, were used to address this computational challenge while preserving the relevant details of the system. In this work, we implement a new multiscale approach based on the SIRAH force field to gain deeper insights into the dynamics of the LacI-DNA system. Our methodology allows for a dual resolution treatment of the solute and solvent, explicitly representing the protein, DNA, and solvent environment without compromising the AA region. Starting from the P1 loop configuration in an undertwisted conformation, we were able to observe the transition to the more stable overtwisted state. Additionally, a detailed characterization of the conformational space sampled by the DNA loop was done. In agreement with experimental and theoretical evidence, we observed the transient formation of kinks at the loop, which were stabilized by the presence of counterions at the minor groove. We also show that the loop's intrinsic flexibility can account for reported FRET measurements and bent conformations required to bind the CAP transcription factor.

Weak operator binding enhances simulated Lac repressor-mediated DNA looping

The 50th anniversary of Biopolymers coincides closely with the like celebration of the discovery of the Escherichia coli (lac) lactose operon, a classic genetic system long used to illustrate the influence of biomolecular structure on function. The looping of DNA induced by the binding of the Lac repressor protein to sequentially distant operator sites on DNA continues to serve as a paradigm for understanding long-range genomic communication. Advances in analyses of DNA structures and in incorporation of proteins in computer simulations of DNA looping allow us to address long-standing questions about the role of protein-mediated DNA loop formation in transcriptional control. Here we report insights gained from studies of the sequence-dependent contributions of the natural lac operators to Lac repressor-mediated DNA looping. Novel superposition of the ensembles of protein-bound operator structures derived from NMR measurements reveals variations in DNA folding missed in conventional structural alignments. The changes in folding affect the predicted ease with which the repressor induces loop formation and the ways that DNA closes between the protein headpieces. The peeling of the auxiliary operators away from the repressor enhances the formation of loops with the 92-bp wildtype spacing and hints of a structural reason behind their weak binding.

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.

Probing DNA topology using tethered particle motion

Transcription factors mediate the formation of nucleoprotein complexes that are critical for efficient regulation of epigenetic switches. In these complexes, DNA is frequently bent or looped by the protein; other times, strong interactions lead the DNA to fully wrap the regulatory protein(s). The equilibrium between the bending, looping, full and partial wrapping of DNA governs the level of transcriptional regulation and is tuned by biophysical parameters. Characterization of the structure, kinetics, and thermodynamics of formation of such nucleoprotein complexes is fundamental to the understanding of the molecular mechanisms that underlie the operation of the genetic switches controlled by them. Here, we describe in detail how to perform tethered particle motion experiments aimed at understanding how protein-DNA interactions influence the formation and breakdown of these regulatory complexes.

FLEXIBILITY AND CONTROL OF PROTEIN–DNA LOOPS

Protein–DNA loops are essential for efficient transcriptional repression and activation. The geometry and stability of the archetypal Lac repressor tetramer (LacI)–DNA loop were investigated using designed hyperstable loops containing lac operators bracketing a sequence-directed bend. Electrophoretic mobility shift assays, DNA cyclization, and bulk and single-molecule fluorescence resonance energy transfer (FRET) demonstrate that the DNA sequence controls whether the LacI–DNA loop forms a compact loop with positive writhe or an open loop with little writhe. Monte Carlo methods for simulation of DNA ring closure were extended to DNA loops, including treatment of variable protein hinge angles. The observed distribution of topoisomer products upon cyclization provides a strong constraint on possible models. The experiments and modeling imply that LacI–DNA can adopt a wide range of geometries but has a strong intrinsic preference for an open form. The flexibility of LacI helps explain in vivo observations that DNA looping is less sensitive to DNA length and shape than that expected from the physical properties of DNA. While DNA cyclization suggests two pools of precursor loops for the 9C14 construct, single-molecule FRET demonstrates a single population. This discrepancy suggests that the LacI–DNA structure is strongly influenced by flanking DNA.

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.

Wrapped-around models for the lac operon complex

The protein-DNA complex, involved in the lac operon of enteric bacteria, is paradigmatic in understanding the extent of DNA bending and plasticity due to interactions with protein assemblies acting as DNA regulators. For the lac operon, two classes of structures have been proposed: 1), with the protein tetramer lying away from the DNA loop (wrapped-away model); and 2), with the protein tetramer lying inside the DNA loop (wrapped-around model). A recently developed electrostatic analytical model shows that the size and net charge of the Lac protein tetramer allow the bending of DNA, which is consistent with another wrapped-around model from the literature. Coarse-grained models, designed based on this observation, are extensively investigated and show three kinds of wrapped-around arrangements of DNA and a lower propensity for wrapped-away configurations. Molecular dynamics simulations of an all-atom model, built on the basis of the most tightly collapsed coarse-grained model, show that most of the DNA double-helical architecture is maintained in the region between O3 and O1 DNA operators, that the DNA distortion is concentrated in the chain beyond the O1 operator, and that the protein tetramer can adapt the N-terminal domains to the DNA tension.

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.

Interconvertible lac repressor-DNA loops revealed by single-molecule experiments

At many promoters, transcription is regulated by simultaneous binding of a protein to multiple sites on DNA, but the structures and dynamics of such transcription factor-mediated DNA loops are poorly understood. We directly examined in vitro loop formation mediated by Escherichia coli lactose repressor using single-molecule structural and kinetics methods. Small (∼150 bp) loops form quickly and stably, even with out-of-phase operator spacings. Unexpectedly, repeated spontaneous transitions between two distinct loop structures were observed in individual protein–DNA complexes. The results imply a dynamic equilibrium between a novel loop structure with the repressor in its crystallographic “V” conformation and a second structure with a more extended linear repressor conformation that substantially lessens the DNA bending strain. The ability to switch between different loop structures may help to explain how robust transcription regulation is maintained even though the mechanical work required to form a loop may change substantially with metabolic conditions.

Some proteins that regulate DNA transcription do so by binding simultaneously to two separated sites on the DNA molecule, forming a DNA loop. Although such loops are common, many of their features are poorly characterized. Of particular interest is the question of how some proteins accommodate the formation of loops of different sizes, particularly when the loops are small and thus require strong bending (and, in some cases, twisting) of the DNA to form. We observed the shape and behavior of individual DNA molecules bent into tight loops by Lac repressor, a transcription-regulating protein from the bacterium Escherichia coli. Loops were formed in DNA molecules with repressor-binding sites on opposite faces of the DNA double helix almost as readily as in those with sites on the same side, suggesting that the repressor is highly flexible. The DNA can switch back and forth between a tighter and a looser loop structure “on the fly” during the lifetime of a single loop, further evidence that Lac repressor is capable of adopting different shapes that may serve to minimize DNA bending or twisting in loops. The ability of the repressor to readily switch between different loop shapes may allow it to maintain effective control of transcription across situations in which the difficulty of bending or twisting DNA changes substantially.

A large-scale conformational change in a transcription factor protein allows DNA loops to dynamically switch between alternative conformations that may contribute to robust transcription regulation.

Formation of a wrapped DNA-protein interface: experimental characterization and analysis of the large contributions of ions and water to the thermodynamics of binding IHF to H' DNA

To characterize driving forces and driven processes in formation of a large-interface, wrapped protein-DNA complex analogous to the nucleosome, we have investigated the thermodynamics of binding the 34-base pair (bp) H' DNA sequence to the Escherichia coli DNA-remodeling protein integration host factor (IHF). Isothermal titration calorimetry and fluorescence resonance energy transfer are applied to determine effects of salt concentration [KCl, KF, K glutamate (KGlu)] and of the excluded solute glycine betaine (GB) on the binding thermodynamics at 20 degrees C. Both the binding constant K(obs) and enthalpy Delta H degrees (obs) depend strongly on [salt] and anion identity. Formation of the wrapped complex is enthalpy driven, especially at low [salt] (e.g., Delta H(o)(obs)=-20.2 kcal x mol(-1) in 0.04 M KCl). Delta H degrees (obs) increases linearly with [salt] with a slope (d Delta H degrees (obs)/d[salt]), which is much larger in KCl (38+/-3 kcal x mol(-1) M(-1)) than in KF or KGlu (11+/-2 kcal x mol(-1) M(-1)). At 0.33 M [salt], K(obs) is approximately 30-fold larger in KGlu or KF than in KCl, and the [salt] derivative SK(obs)=dlnK(obs)/dln[salt] is almost twice as large in magnitude in KCl (-8.8+/-0.7) as in KF or KGlu (-4.7+/-0.6). A novel analysis of the large effects of anion identity on K(obs), SK(obs) and on Delta H degrees (obs) dissects coulombic, Hofmeister, and osmotic contributions to these quantities. This analysis attributes anion-specific differences in K(obs), SK(obs), and Delta H degrees (obs) to (i) displacement of a large number of water molecules of hydration [estimated to be 1.0(+/-0.2)x10(3)] from the 5340 A(2) of IHF and H' DNA surface buried in complex formation, and (ii) significant local exclusion of F(-) and Glu(-) from this hydration water, relative to the situation with Cl(-), which we propose is randomly distributed. To quantify net water release from anionic surface (22% of the surface buried in complexation, mostly from DNA phosphates), we determined the stabilizing effect of GB on K(obs): dlnK(obs)/d[GB]=2.7+/-0.4 at constant KCl activity, indicating the net release of ca. 150 H(2)O molecules from anionic surface.

Simultaneous binding of minor groove binder and intercalator to dodecamer DNA: importance of relative orientation of donor and acceptor in FRET

In the present study, steady-state, picosecond time-resolved fluorescence and polarization gated anisotropy have been used to establish simultaneous binding of an intercalator (ethidium bromide, EtBr) and a minor groove binder (Hoeschst 33258, H258) to a dodecamer DNA of specific sequence. The Förster resonance energy transfer (FRET) studies between the dyes H258 (donor) and EtBr (acceptor) in the dodecamer, where the ligands have a particular relative orientation of the transition dipoles, in contrast to the cases in sodium dodecyl sulfate (SDS) micelles and larger genomic DNA, where the orientations are random, reveal the effect of the binding geometry of the ligands in the constrained environment. Our study establishes that reconsideration of the value of the orientation factor (kappa2) is crucial for correct estimation of the donor-acceptor distance when the ligands are simultaneously bound to a specific region of biological macromolecules.

Effect of supercoiling on formation of protein-mediated DNA loops

DNA loop formation is one of several mechanisms used by organisms to regulate genes. The free energy of forming a loop is an important factor in determining whether the associated gene is switched on or off. In this paper we use an elastic rod model of DNA to determine the free energy of forming short (50-100 basepair), protein mediated DNA loops. Superhelical stress in the DNA of living cells is a critical factor determining the energetics of loop formation, and we explicitly account for it in our calculations. The repressor protein itself is regarded as a rigid coupler; its geometry enters the problem through the boundary conditions it applies on the DNA. We show that a theory with these ingredients is sufficient to explain certain features observed in modulation of in vivo gene activity as a function of the distance between operator sites for the lac repressor. We also use our theory to make quantitative predictions for the dependence of looping on superhelical stress, which may be testable both in vivo and in single-molecule experiments such as the tethered particle assay and the magnetic bead assay.

Compaction of DNA on nanoscale three-dimensional templates

There exist several important in vivo examples, where a DNA chain is compacted on interacting with nanoscale objects such as proteins, thereby forming complexes with a well defined molecular architecture. One of the well known manifestations of such a natural organization of a semi-flexible DNA chain on nanoscale objects is hierarchical DNA molecule assembly into a chromosome, which is mediated by cationic histone proteins at the first stages of compaction. The biological importance of this and other natural nanostructural organizations of the DNA molecule has inspired many theoretical and numerical studies to gain physical insight into this problem. On the other hand, the experimental model systems containing DNA and nanoobjects, which are important to extend our knowledge beyond natural systems, were almost unavailable until the last decade. Accelerating progress in nanoscale chemistry and materials science has brought about various nanoscale three-dimensional structures such as dendrimers, nanoparticles, and nanotubes, and thus has provided a basis for the next important step in creating novel DNA-containing nanostructures, modelling of natural DNA compaction, and verification of accumulated theoretical predictions on the interaction between DNA and nanoscale templates. This review is written to highlight this early stage of nano-inspired progress and it is focused on physico-chemical and biophysical experimental investigations as well as theoretical and numerical studies dedicated to the compaction of DNA on nanoscale three-dimensional templates.

Modeling the Lac repressor-operator assembly: the influence of DNA looping on Lac repressor conformation

Repression of transcription of the Escherichia coli Lac operon by the Lac repressor (LacR) is accompanied by the simultaneous binding of LacR to two operators and the formation of a DNA loop. A recently developed theory of sequence-dependent DNA elasticity enables one to relate the fine structure of the LacR-DNA complex to a wide range of heretofore-unconnected experimental observations. Here, that theory is used to calculate the configuration and free energy of the DNA loop as a function of its length and base-pair sequence, its linking number, and the end conditions imposed by the LacR tetramer. The tetramer can assume two types of conformations. Whereas a rigid V-shaped structure is observed in the crystal, EM images show extended forms in which two dimer subunits are flexibly joined. Upon comparing our computed loop configurations with published experimental observations of permanganate sensitivities, DNase I cutting patterns, and loop stabilities, we conclude that linear DNA segments of short-to-medium chain length (50-180 bp) give rise to loops with the extended form of LacR and that loops formed within negatively supercoiled plasmids induce the V-shaped structure.

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.

Energetic differences between the specific binding of a 40bp DNA duplex and the lac promoter to lac repressor protein

The energetics of LRP binding to a 104 bp lac promoter determined from ITC measurements were compared to the energetics of binding to a shorter 40 bp DNA duplex with the 21 bp promoter binding site sequence. The promoter binding affinity of 2.47 +/- 0.0 1x 10(7) M(-1) was higher than the DNA binding affinity of 1.81 +/- 0.67 x 10(7) M(-1) while the binding enthalpy of -804 +/- 41 kJ mol(-1) was lower than that of the DNA binding enthalpy of -145 +/- 16 kJ mol(-1) at 298.15 K. Both the promoter and DNA binding reactions were exothermic in phosphate buffer but endothermic in Tris buffer that showed the transfer of four protons to LRP in the former reaction but only two in the latter. A more complicated dependence of these parameters on temperature was observed for promoter binding. These energetic differences are attributable to additional LRP-promoter interactions from wrapping of the promoter around the LRP.

DNA trajectory in the Gal repressosome

The Gal repressosome is a higher-order nucleoprotein complex that represses transcription of the gal operon in Escherichia coli. During the repressosome assembly, a DNA loop is formed by the interaction of two GalR dimers, bound to two spatially separated operators, OE and OI, flanking the gal promoters. Structure-based genetic analysis indicated that GalR homodimers interact directly and form a V-shaped stacked tetramer in repressosome, further stabilized by HU binding to an architecturally critical position on the DNA. In this scheme of GalR tetramerization, the alignment of the operators in the DNA loop could be in either parallel (PL) or antiparallel (AL) mode. As each mode can have two alternative geometries differing in the mutual stacking of the OE- and OI-bound GalR dimers, it is possible to have four different DNA trajectories in the repressosome. Feasibilities of these trajectories were tested by in vitro transcription repression assays, first by isolating GalR mutants with altered operator specificity and then by constructing four different potential loops with mutant GalR heterodimers bound to specifically designed hybrid operators in such a way as to give rise to only one of the four putative trajectories. Results show that OE and OI adopt a mutual antiparallel orientation in an under-twisted DNA loop, consistent with the energetically optimal structural model. In this structure the center of the HU-binding site is located at the apex of the DNA loop. The approach reported here can be used to distinguish between otherwise indistinguishable DNA trajectories in complex nucleoprotein machines.

Modelling DNA loops using continuum and statistical mechanics

The classical Kirchhoff elastic-rod model applied to DNA is extended to account for sequence-dependent intrinsic twist and curvature, anisotropic bending rigidity, electrostatic force interactions, and overdamped Brownian motion in a solvent. The zero-temperature equilibrium rod model is then applied to study the structural basis of the function of the lac repressor protein in the lac operon of Escherichia coli. The structure of a DNA loop induced by the clamping of two distant DNA operator sites by lac repressor is investigated and the optimal geometries for the loop of length 76 bp are predicted. Further, the mimicked binding of catabolite gene activator protein (CAP) inside the loop provides solutions that might explain the experimentally observed synergy in DNA binding between the two proteins. Finally, a combined Monte Carlo and Brownian dynamics solver for a worm-like chain model is described and a preliminary analysis of DNA loop-formation kinetics is presented.

DNA spools under tension

DNA spools, structures in which DNA is wrapped and helically coiled onto itself or onto a protein core, are ubiquitous in nature. We develop a general theory describing the nonequilibrium behavior of DNA spools under linear tension. Two puzzling and seemingly unrelated recent experimental findings, the sudden quantized unwrapping of nucleosomes and that of DNA toroidal condensates under tension, are theoretically explained and shown to be of the same origin. The study provides new insights into nucleosome and chromatin fiber stability and dynamics.

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.

Genetic analysis of GalR tetramerization in DNA looping during repressosome assembly

The Gal repressosome is a nucleoprotein complex consisting of 2 GalR dimers, 1 HU, and 1 DNA loop, which represses the transcription of the gal operon. We have adopted a structure-based genetic approach to complement ongoing physical studies of the complex. Homology-based and subsequent alanine-scanning mutageneses suggest that five residues in the DNA-distal subdomain of GalR dimer are important for repressosome formation. A further analysis of these and intragenic suppressors of looping-defective GalR mutants as well as gain-of-function mutants that permit repressosome assembly in the absence of HU show that GalR dimers contact each other in the repressosome in a partially stacked configuration.

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.

Charge neutralization and DNA bending by the Escherichia coli catabolite activator protein

We are interested in the role of asymmetric phosphate neutralization in DNA bending induced by proteins. We describe an experimental estimate of the actual electrostatic contribution of asymmetric phosphate neutralization to the bending of DNA by the Escherichia coli catabolite activator protein (CAP), a prototypical DNA-bending protein. Following assignment of putative electrostatic interactions between CAP and DNA phosphates based on X-ray crystal structures, appropriate phosphates in the CAP half-site DNA were chemically neutralized by methylphosphonate substitution. DNA shape was then evaluated using a semi-synthetic DNA electrophoretic phasing assay. Our results confirm that the unmodified CAP DNA half-site sequence is intrinsically curved by 26 degrees in the direction enhanced in the complex with protein. In the absence of protein, neutralization of five appropriate phosphates increases DNA curvature to 32 degrees (approximately 23% increase), in the predicted direction. Shifting the placement of the neutralized phosphates changes the DNA shape, suggesting that sequence-directed DNA curvature can be modified by the asymmetry of phosphate neutralization. We suggest that asymmetric phosphate neutralization contributes favorably to DNA bending by CAP, but cannot account for the full DNA deformation.

Radiolysis of lac repressor by gamma-rays and heavy ions: a two-hit model for protein inactivation

Upon gamma-ray or argon ion irradiation of the lac repressor protein, its peptide chain is cleaved and the protein loses its lac operator-binding activity, as shown respectively by polyacrylamide gel electrophoresis and retardation gel electrophoresis. We developed phenomenological models that satisfactorily account for the experimental results: the peptide chain cleavage model considers that the average number of chain breaks per protomer is proportional to the irradiation dose and that the distribution of the number of breaks per protomer obeys Poisson's law. The repressor inactivation model takes into account the quaternary structure (a dimer of dimer) and the organization of the repressor in domains (two DNA binding sites, one per dimer). A protomer is inactivated by at least two different radiation-induced damages. A dimer is inactivated when at least one of the two protomers is inactivated. A tetramer is inactivated when both dimers are inactivated. From the combination of both models, we can deduce that chain cleavage cannot account for the protein inactivation, which should mainly result from oxidation of amino acid side chains. Indeed, particularly oxidizable and accessible amino acids (Tyr, His) are involved in the DNA binding process.

The type IIs restriction endonuclease BspMI is a tetramer that acts concertedly at two copies of an asymmetric DNA sequence

Type IIs endonucleases recognize asymmetric DNA sequences and cleave both strands at fixed positions downstream of the sequence. Many type IIs enzymes, including BspMI, cleave substrates with two sites more rapidly than those with one site. They usually act sequentially on DNA with two sites, but BspMI converted such a substrate directly to the final products cut at both sites. The BspMI endonuclease was found to be a tetramer, in contrast to the monomeric structures for many type IIs enzymes. No change in subunit association occurred during the BspMI reaction. Plasmids with two BspMI sites were cleaved in cis, in reactions spanning sites in the same DNA, even when the sites were separated by just 38 bp. Plasmids with one BspMI site were cleaved in trans, with the enzyme bridging sites in separate DNA molecules: these slow reactions could be accelerated by adding a second DNA with the recognition sequence. Thus, whereas many type IIs enzymes dimerize before cleaving DNA, a process facilitated by two recognition sites in cis, the BspMI tetramer binds two copies of its recognition sequence before cleaving the DNA in both strands at both sites.

Specific and non-specific interactions of integration host factor with DNA: thermodynamic evidence for disruption of multiple IHF surface salt-bridges coupled to DNA binding

Site-specific DNA binding of architectural protein integration host factor (IHF) is involved in formation of functional multiprotein-DNA assemblies in Escherichia coli, while non-specific binding of IHF and other histone-like proteins serves to structure the nucleoid. Here, we report an isothermal titration calorimetry study of the thermodynamics of binding IHF to a 34 bp fragment composed entirely of the specific H' site from lambda-phage DNA. At low to moderate [K(+)] (60-100 mM), strong competition is observed between specific and non-specific binding as a result of a low specificity ratio (approximately 10(2)) and a very small non-specific site size. In this [K(+)] range, both specific and non-specific binding are enthalpy-driven, with large negative enthalpy, entropy and heat capacity changes and binding constants that are insensitive to [K(+)]. Above 100 mM K(+), only specific binding is observed, and both the binding constant and the magnitudes of enthalpy, entropy and heat capacity changes all decrease strongly with increasing [K(+)]. When interpreted in the context of the structure of the specific complex, the thermodynamics provide compelling evidence for a previously unrecognized design principle by which proteins that form extensive binding interfaces with nucleic acids control binding constants, binding site sizes and effects of temperature and ion concentrations on stability and specificity. We propose that up to 22 of the 23 IHF cationic side-chains that are located within 6 A of DNA phosphate oxygen atoms in the complex, are masked in the absence of DNA by pairing with anionic carboxylate groups in intramolecular salt-bridges (dehydrated ion-pairs). These salt-bridges increase in stability with increasing temperature and decreasing [K(+)]. To explain the unusual thermodynamics of IHF-DNA interactions, we propose that both specific and non-specific binding at low [K(+)] require disruption of salt-bridges (as many as 18 for specific binding) whereupon many of the unmasked charged groups hydrate and the cationic groups interact with DNA. From structural or thermodynamic parallels with IHF, we propose that large-scale coupling of disruption of protein salt-bridges to DNA binding is significant for other large-interface DNA wrapping proteins including the nucleosome, lac repressor core tetramer, RNA polymerase core protein, HU and SSB.

Strong DNA binding by covalently linked dimeric Lac headpiece: evidence for the crucial role of the hinge helices

The combined structural and biochemical studies on Lac repressor bound to operator DNA have demonstrated the central role of the hinge helices in operator bending and the induction mechanism. We have constructed a covalently linked dimeric Lac-headpiece that binds DNA with four orders of magnitude higher affinity as compared with the monomeric form. This enabled a detailed biochemical and structural study of Lac binding to its cognate wild-type and selected DNA operators. The results indicate a profound contribution of hinge helices to the stability of the protein-DNA complex and highlight their central role in operator recognition. Furthermore, protein-DNA interactions in the minor groove appear to modulate hinge helix stability, thus accounting for affinity differences and protein-induced DNA bending among the various operator sites. Interestingly, the in vitro DNA-binding affinity of the reported dimeric Lac construct can de readily modulated by simple adjustment of redox conditions, thus rendering it a potential artificial gene regulator.

Alternative geometries of DNA looping: an analysis using the SfiI endonuclease

Many processes are governed by proteins that bind to separate sites in DNA and loop out the intervening DNA, but the geometries of the loops have seldom been determined. The SfiI endonuclease cleaves DNA after interacting with two recognition sites, and is a favourable system for the analysis of DNA looping. A gel-shift assay was used here to examine the binding of SfiI to a series of linear DNA molecules containing two SfiI sites separated by 109-170 base-pairs. The complexes in which SfiI trapped a loop by binding to two sites in the same DNA were separated from the complexes containing SfiI bound to separate DNA molecules. Step-wise changes in the inter-site spacing generated two forms of the looped complex with different electrophoretic mobilities. The yields of each looped complex and the complexes from intermolecular synapses all varied cyclically with the inter-site spacing, with similar periodicities ( approximately 10.5 base-pairs) but with different phases. One looped complex predominated whenever the DNA between the sites needed to be underwound in order to produce the correct helical orientation of the binding sites. The other looped complex predominated whenever the intervening DNA needed to be overwound. We conclude that the former has trapped a right-handed loop with a negative node and the latter a left-handed loop with a positive node.