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.     

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.     

Mechanism of promoter repression by Lac repressor–DNA loops

The Escherichia coli lactose (lac) operon encodes the first genetic switch to be discovered, and lac remains a paradigm for studying negative and positive control of gene expression. Negative control is believed to involve competition of RNA polymerase and Lac repressor for overlapping binding sites. Contributions to the local Lac repressor concentration come from free repressor and repressor delivered to the operator from remote auxiliary operators by DNA looping. Long-standing questions persist concerning the actual role of DNA looping in the mechanism of promoter repression. Here, we use experiments in living bacteria to resolve four of these questions. We show that the distance dependence of repression enhancement is comparable for upstream and downstream auxiliary operators, confirming the hypothesis that repressor concentration increase is the principal mechanism of repression loops. We find that as few as four turns of DNA can be constrained in a stable loop by Lac repressor. We show that RNA polymerase is not trapped at repressed promoters. Finally, we show that constraining a promoter in a tight DNA loop is sufficient for repression even when promoter and operator do not overlap.     

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.     

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.     

Fluorescence Resonance Energy Transfer over ∼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, ~130 basepairs apart. For 9C14, efficient LacI-induced energy transfer (~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, ~52% transfer was observed. Time-resolved FRET suggests that this complex exists in both closed- and open form populations. Less efficient transfer, ~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.     

Direct demonstration and quantification of long-range DNA looping by the λ bacteriophage repressor

Recently, it was proposed that DNA looping by the λ repressor (CI protein) strengthens repression of lytic genes during lysogeny and simultaneously ensures efficient switching to lysis. To investigate this hypothesis, tethered particle motion experiments were performed and dynamic CI-mediated looping of single DNA molecules containing the λ repressor binding sites separated by 2317 bp (the wild-type distance) was quantitatively analyzed. DNA containing all three intact operators or with mutated o3 operators were compared. Modeling the thermodynamic data established the free energy of CI octamer-mediated loop formation as 1.7 kcal/mol, which decreased to –0.7 kcal/mol when supplemented by a tetramer (octamer+tetramer-mediated loop). These results support the idea that loops secured by an octamer of CI bound at oL1, oL2, oR1 and oR2 operators must be augmented by a tetramer of CI bound at the oL3 and oR3 to be spontaneous and stable. Thus the o3 sites are critical for loops secured by the CI protein that attenuate cI expression.     

The Effect of Nonspecific Binding of Lambda Repressor on DNA Looping Dynamics

The λ repressor (CI) protein-induced DNA loop maintains stable lysogeny, yet allows efficient switching to lysis. Herein, the kinetics of loop formation and breakdown has been characterized at various concentrations of protein using tethered particle microscopy and a novel, to our knowledge, method of analysis. Our results show that a broad distribution of rate constants and complex kinetics underlie loop formation and breakdown. In addition, comparison of the kinetics of looping in wild-type DNA and DNA with mutated o3 operators showed that these sites may trigger nucleation of nonspecific binding at the closure of the loop. The average activation energy calculated from the rate constant distribution is consistent with a model in which nonspecific binding of CI between the operators shortens their effective separation, thereby lowering the energy barrier for loop formation and broadening the rate constant distribution for looping. Similarly, nonspecific binding affects the kinetics of loop breakdown by increasing the number of loop-securing protein interactions, and broadens the rate constant distribution for this reaction. Therefore, simultaneous increase of the rate constant for loop formation and reduction of that for loop breakdown stabilizes lysogeny. Given these simultaneous changes, the frequency of transitions between the looped and the unlooped state remains nearly constant. Although the loop becomes more stable thermodynamically with increasing CI concentration, it still opens periodically, conferring sensitivity to environmental changes, which may require switching to lytic conditions.     

Computational analysis of looping of a large family of highly bent DNA by LacI

Sequence-dependent intrinsic curvature of DNA influences looping by regulatory proteins such as LacI and NtrC. Curvature can enhance stability and control shape, as observed in LacI loops formed with three designed sequences with operators bracketing an A-tract bend. We explore geometric, topological, and energetic effects of curvature with an analysis of a family of highly bent sequences, using the elastic rod model from previous work. A unifying straight-helical-straight representation uses two phasing parameters to describe sequences composed of two straight segments that flank a common helically supercoiled segment. We exercise the rod model over this two-dimensional space of phasing parameters to evaluate looping behaviors. This design space is found to comprise two subspaces that prefer parallel versus anti-parallel binding topologies. The energetic cost of looping varies from 4 to 12 kT. Molecules can be designed to yield distinct binding topologies as well as hyperstable or hypostable loops and potentially loops that can switch conformations. Loop switching could be a mechanism for control of gene expression. Model predictions for linking numbers and sizes of LacI-DNA loops can be tested using multiple experimental approaches, which coupled with theory could address whether proteins or DNA provide the observed flexibility of protein-DNA loops.     

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.     

Bead size effects on protein-mediated DNA looping in tethered-particle motion experiments

Tethered particle motion (TPM) has become an important tool for single-molecule studies of biomolecules; however, concerns remain that the method may alter the dynamics of the biophysical process under study. We investigate the effect of the attached microsphere on an illustrative biological example: the formation and breakdown of protein-mediated DNA loops in the lac repressor system. By comparing data from a conventional TPM experiment with 800 nm polystyrene beads and dark-field TPM using 50 nm Au nanoparticles, we found that the lifetimes of the looped and unlooped states are only weakly modified, less than two-fold, by the presence of the large bead. This is consistent with our expectation of weak excluded-volume effects and hydrodynamic surface interactions from the cover glass and microsphere.     

DNA-binding properties of a lac repressor mutant incapable of forming tetramers

The interaction of proteins bound to sites widely separated on the genome is a recurrent motif in both prokaryotic and eukaryotic regulatory systems. Lac repressor mediates the formation of "DNA loops" by the simultaneous interaction of a single protein tetramer with two DNA-binding sites. The DNA-binding properties of a Lac repressor mutant (LacIadi) deficient in the association of protein dimers to tetramers was investigated. The results of quantitative footprint and gel mobility-shift titrations suggest that the wild-type Lac repressor (LacI+) binds cooperatively to two operator sites separated by 11 helical turns on a linear DNA restriction fragment by the formation of a "looped complex." LacIadi binds to this two-site operator non-cooperatively and without formation of a looped complex. These results demonstrate that the dimer-tetramer association of LacI+ is directly responsible for its cooperative binding and its ability to mediate formation of a looped complex. The Iadi mutation disrupts the monomer-dimer as well as eliminating the dimer-tetramer association equilibria while the DNA binding affinity of LacIadi to a single site is unchanged relative to the wild-type protein. These results suggest that DNA binding and dimer-tetramer association are functionally unlinked. The similarity of the DNA-binding properties of LacIadi and Gal repressor, a protein believed to function by mediating the formation of a looped complex, are discussed.