DNA modeling reveals an extended lac repressor conformation in classic in vitro binding assays

Protein-mediated DNA looping, such as that induced by the lactose repressor (LacI) of Escherichia coli, is a well-known gene regulation mechanism. Although researchers have given considerable attention to DNA looping by LacI, many unanswered questions about this mechanism, including the role of protein flexibility, remain. Recent single-molecule observations suggest that the two DNA-binding domains of LacI are capable of splaying open about the tetramerization domain into an extended conformation. We hypothesized that if recent experiments were able to reveal the extended conformation, it is possible that such structures occurred in previous studies as well. In this study, we tested our hypothesis by reevaluating two classic in vitro binding assays using a computational rod model of DNA. The experiments and computations evaluate the looping of both linear DNA and supercoiled DNA minicircles over a broad range of DNA interoperator lengths. The computed energetic minima align well with the experimentally observed interoperator length for optimal loop stability. Of equal importance, the model reveals that the most stable loops for linear DNA occur when LacI adopts the extended conformation. In contrast, for DNA minicircles, optimal stability may arise from either the closed or the extended protein conformation depending on the degree of supercoiling and the interoperator length.     

Looping charged elastic rods: applications to protein-induced DNA loop formation

We analyze looping of thin charged elastic filaments under applied torques and end forces, using the solution of linear elasticity theory equations. In application to DNA, we account for its polyelectrolyte character and charge renormalization, calculating electrostatic energies stored in the loops. We argue that the standard theory of electrostatic persistence is only valid when the loop’s radius of curvature and close-contact distance are much larger than the Debye screening length. We predict that larger twist rates are required to trigger looping of charged rods as compared with neutral ones. We then analyze loop shapes formed on charged filaments of finite length, mimicking DNA looping by proteins with two DNA-binding domains. We find optimal loop shapes at different salt amounts, minimizing the sum of DNA elastic, DNA electrostatic, and protein elastic energies. We implement a simple model where intercharge repulsions do not affect the loop shape directly but can choose the energy-optimized shape from the allowed loop types. At low salt concentrations more open loops are favored due to enhanced repulsion of DNA charges, consistent with the results of computer simulations on formation of DNA loops by lac repressor. Then, we model the precise geometry of DNA binding by the lac tetramer and explore loop shapes, varying the confined DNA length and protein opening angle. The characteristics of complexes obtained, such as the total loop energy, stretching forces required to maintain its shape, and the reduction of electrostatic energy with increment of salt, are in good agreement with the outcomes of more elaborate numerical calculations for lac-repressor-induced DNA looping.     

Statistical-mechanical theory of DNA looping

The lack of a rigorous analytical theory for DNA looping has caused many DNA-loop-mediated phenomena to be interpreted using theories describing the related process of DNA cyclization. However, distinctions in the mechanics of DNA looping versus cyclization can have profound quantitative effects on the thermodynamics of loop closure. We have extended a statistical mechanical theory recently developed for DNA cyclization to model DNA looping, taking into account protein flexibility. Notwithstanding the underlying theoretical similarity, we find that the topological constraint of loop closure leads to the coexistence of multiple classes of loops mediated by the same protein structure. These loop topologies are characterized by dramatic differences in twist and writhe; because of the strong coupling of twist and writhe within a loop, DNA looping can exhibit a complex overall helical dependence in terms of amplitude, phase, and deviations from uniform helical periodicity. Moreover, the DNA-length dependence of optimal looping efficiency depends on protein elasticity, protein geometry, and the presence of intrinsic DNA bends. We derive a rigorous theory of loop formation that connects global mechanical and geometric properties of both DNA and protein and demonstrates the importance of protein flexibility in loop-mediated protein-DNA interactions.     

Multimerization of the Drosophila zeste protein is required for efficient DNA binding.

The Drosophila zeste protein forms multimeric species in vitro through its C-terminal domain. Multimerization is required for efficient binding to DNA containing multiple recognition sequences and increasing the number of binding sites stimulates binding in a cooperative manner. Mutants that can only form dimers still bind to a dimeric site, but with lower affinity. Mutations or progressive deletions from the C-terminal show that when even dimer formation is prevented, DNA-binding activity is lost. Surprisingly, binding activity is regained with larger deletions that leave only the DNA-binding domain. Additional protein sequences apparently inhibit DNA binding unless they permit multimerization. The DNA-binding domain peptides bind strongly even to isolated recognition sequences and they bind as monomers. The ability of various zeste peptides to stimulate white gene expression in vivo shows that multimeric forms are the functional species of the zeste product in vivo. The DNA-binding domain peptide binds well to DNA in vitro, but it cannot stimulate white gene expression in vivo. This failure may reflect the need for an activation domain or it may be caused by indiscriminate binding of this peptide to non-functional isolated sites. Multimerization increases binding specificity, selecting only sites with multiple recognition sequences.     

Quantitative comparison of DNA looping in vitro and in vivo: chromatin increases effective DNA flexibility at short distances

The probability that two sites on a linear DNA molecule will contact each other by looping depends on DNA flexibility. Although the flexibility of naked DNA in vitro is well characterized, looping in chromatin is poorly understood. By extending existing theory, we present a single equation that describes DNA looping over all distances. We also show that DNA looping in vitro can be measured accurately by FLP recombination between sites from 74 bp to 15 kb apart. In agreement with previous work, a persistence length of 50 nm was determined. FLP recombination of the same substrates in mammalian cells showed that chromatin increases the flexibility of DNA at short distances, giving an apparent persistence length of 27 nm.     

Real-time observation of DNA looping dynamics of Type IIE restriction enzymes NaeI and NarI

Many restriction enzymes require binding of two copies of a recognition sequence for DNA cleavage, thereby introducing a loop in the DNA. We investigated looping dynamics of Type IIE restriction enzymes NaeI and NarI by tracking the Brownian motion of single tethered DNA molecules. DNA containing two endonuclease recognition sites spaced a few 100 bp apart connect small polystyrene beads to a glass surface. The position of a bead is tracked through video microscopy. Protein-mediated looping and unlooping is then observed as a sudden specific change in Brownian motion of the bead. With this method we are able to directly follow DNA looping kinetics of single protein–DNA complexes to obtain loop stability and loop formation times. We show that, in the absence of divalent cations, NaeI induces DNA loops of specific size. In contrast, under these conditions NarI mainly creates non-specific loops, resulting in effective DNA compaction for higher enzyme concentrations. Addition of Ca²⁺ increases the NaeI-DNA loop lifetime by two orders of magnitude and stimulates specific binding by NarI. Finally, for both enzymes we observe exponentially distributed loop formation times, indicating that looping is dominated by (re)binding the second recognition site.     

Identification of the DNA-binding domain of the FLP recombinase

We have subjected the FLP protein of the 2-micron plasmid to partial proteolysis by proteinase K and have found that FLP can be digested into two major proteinase K-resistant peptides of 21 and 13 kDa, respectively. The 21-kDa peptide contains a site-specific DNA-binding domain that binds to the FLP recognition target (FRT) site with an affinity similar to that observed for the native FLP protein. This peptide can induce DNA bending upon binding to a DNA fragment containing the FRT site, but the angle of the bend (approximately 24 degrees) is smaller in magnitude than that induced by the native FLP protein (60 degrees). The additional DNA bending induced by the interaction between two native FLP molecules bound to the FRT site is not observed with the 21-kDa DNA-binding peptide. Amino-terminal sequencing has been used to map this peptide to an internal region of FLP that begins at residue Leu-148. It is likely that the DNA-binding peptide includes the catalytic site of the FLP protein.     

Dynamics of the Mrf-2 DNA-binding domain free and in complex with DNA

Mrf-2 is a member of a new class of DNA-binding proteins known as the AT-rich interaction domain family or ARID. Chemical shift indices and characteristic NOE values indicate that the three-dimensional structure of the Mrf-2 ARID in complex with DNA is nearly identical to that of the free protein. The backbone dynamics of the Mrf-2 domain free and in complex with DNA have been characterized by (15)N NMR relaxation measurements and model-free analysis. Chemical shift perturbations and dynamic studies suggest that two flexible interhelical loops, the flexible C-terminal tail, and one alpha-helix are involved in DNA recognition, indicating the importance of protein dynamics in DNA binding. Some well-structured regions, in particular the putative DNA-contacting helix, in Mrf-2 show a decrease in the order parameters (S(2)) upon complex formation. The less well-structured loops and the unstructured C-terminus show reduced flexibility upon DNA binding. In addition, the model-free analysis indicates motions on the picosecond to nanosecond and micro- to millisecond time scales at the DNA-binding surface of the bound Mrf-2 ARID, suggesting a model where interactions between the protein and DNA are highly dynamic.     

Interplay of Protein and DNA Structure Revealed in Simulations of the lac Operon

The E. coli Lac repressor is the classic textbook example of a protein that attaches to widely spaced sites along a genome and forces the intervening DNA into a loop. The short loops implicated in the regulation of the lac operon suggest the involvement of factors other than DNA and repressor in gene control. The molecular simulations presented here examine two likely structural contributions to the in-vivo looping of bacterial DNA: the distortions of the double helix introduced upon association of the highly abundant, nonspecific nucleoid protein HU and the large-scale deformations of the repressor detected in low-resolution experiments. The computations take account of the three-dimensional arrangements of nucleotides and amino acids found in crystal structures of DNA with the two proteins, the natural rest state and deformational properties of protein-free DNA, and the constraints on looping imposed by the conformation of the repressor and the orientation of bound DNA. The predicted looping propensities capture the complex, chain-length-dependent variation in repression efficacy extracted from gene expression studies and in vitro experiments and reveal unexpected chain-length-dependent variations in the uptake of HU, the deformation of repressor, and the folding of DNA. Both the opening of repressor and the presence of HU, at levels approximating those found in vivo, enhance the probability of loop formation. HU affects the global organization of the repressor and the opening of repressor influences the levels of HU binding to DNA. The length of the loop determines whether the DNA adopts antiparallel or parallel orientations on the repressor, whether the repressor is opened or closed, and how many HU molecules bind to the loop. The collective behavior of proteins and DNA is greater than the sum of the parts and hints of ways in which multiple proteins may coordinate the packaging and processing of genetic information.     

Monovalent cations regulate DNA sequence recognition by 434 repressor

The bacteriophage 434 repressor distinguishes between its six naturally occurring binding sites using indirect readout. In indirect readout, sequence-dependent differences in the structure and flexibility of non-contacted bases in a protein's DNA-binding site modulate the affinity of DNA for protein. The conformation and flexibility of a DNA sequence can be influenced by the interaction of the DNA bases or backbone with solution components. We examined the effect of changing the cation-type present in solution on the stability and structure of 434 repressor complexes with wild-type and mutant OR1 and OR3, binding sites that differ in their contacted and non-contacted base sequences. We find that the affinity of repressor for OR1, but not for OR3, depends remarkably on the type and concentration of monovalent cation. Moreover, the formation of a stable, specific repressor-OR1 complex requires the presence of monovalent cations; however, repressor-OR3 complex formation has no such requirement. Changing monovalent cation type alters the ability of repressor to protect OR1, but not OR3, from *OH radical cleavage. Altering the relative length of the poly(dA) x poly(dT) tract in the non-contacted regions of the OR1 and OR3 can reverse the cation sensitivity of repressor's affinities for these two sites. Taken together these findings show that cation-dependent alterations in DNA structure underlies indirect readout of DNA sequence by 434 repressor and perhaps other proteins.     

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

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