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.

Rationally designed coiled-coil DNA looping peptides control DNA topology

Artificial DNA looping peptides were engineered to study the roles of protein and DNA flexibility in controlling the geometry and stability of protein-mediated DNA loops. These LZD (leucine zipper dual-binding) peptides were derived by fusing a second, C-terminal, DNA-binding region onto the GCN4 bZip peptide. Two variants with different coiled-coil lengths were designed to control the relative orientations of DNA bound at each end. Electrophoretic mobility shift assays verified formation of a sandwich complex containing two DNAs and one peptide. Ring closure experiments demonstrated that looping requires a DNA-binding site separation of 310 bp, much longer than the length needed for natural loops. Systematic variation of binding site separation over a series of 10 constructs that cyclize to form 862-bp minicircles yielded positive and negative topoisomers because of two possible writhed geometries. Periodic variation in topoisomer abundance could be modeled using canonical DNA persistence length and torsional modulus values. The results confirm that the LZD peptides are stiffer than natural DNA looping proteins, and they suggest that formation of short DNA loops requires protein flexibility, not unusual DNA bendability. Small, stable, tunable looping peptides may be useful as synthetic transcriptional regulators or components of protein–DNA nanostructures.

LacI-DNA-IPTG loops: equilibria among conformations by single-molecule FRET

The E. coli Lac repressor (LacI) tetramer binds simultaneously to a promoter-proximal DNA binding site (operator) and an auxiliary operator, resulting in a DNA loop, which increases repression efficiency. Induction of the lac operon by allolactose reduces the affinity of LacI for DNA, but induction does not completely prevent looping in vivo. Our previous work on the conformations of LacI loops used a hyperstable model DNA construct, 9C14, that contains a sequence directed bend flanked by operators. Single-molecule fluorescence resonance energy transfer (SM-FRET) on a dual fluorophore-labeled LacI-9C14 loop showed that it adopts a single, stable, high-FRET V-shaped LacI conformation. Ligand-induced changes in loop geometry can affect loop stability, and the current work assesses loop population distributions for LacI-9C14 complexes containing the synthetic inducer IPTG. SM-FRET confirms that the high-FRET LacI-9C14 loop is only partially destabilized by saturating IPTG. LacI titration experiments and FRET fluctuation analysis suggest that the addition of IPTG induces loop conformational dynamics and re-equilibration between loop population distributions that include a mixture of looped states that do not exhibit high-efficiency FRET. The results show that repression by looping even at saturating IPTG should be considered in models for regulation of the operon. We propose that persistent DNA loops near the operator function biologically to accelerate rerepression upon exhaustion of inducer.

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.

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.

Gene repression by minimal lac loops in vivo

The inflexibility of double-stranded DNA with respect to bending and twisting is well established in vitro. Understanding apparent DNA physical properties in vivo is a greater challenge. Here, we exploit repression looping with components of the Escherichia coli lac operon to monitor DNA flexibility in living cells. We create a minimal system for testing the shortest possible DNA repression loops that contain an E. coli promoter, and compare the results to prior experiments. Our data reveal that loop-independent repression occurs for certain tight operator/promoter spacings. When only loop-dependent repression is considered, fits to a thermodynamic model show that DNA twisting limits looping in vivo, although the apparent DNA twist flexibility is 2- to 4-fold higher than in vitro. In contrast, length-dependent resistance to DNA bending is not observed in these experiments, even for the shortest loops constraining <0.4 persistence lengths of DNA. As observed previously for other looping configurations, loss of the nucleoid protein heat unstable (HU) markedly disables DNA looping in vivo. Length-independent DNA bending energy may reflect the activities of architectural proteins and the structure of the DNA topological domain. We suggest that the shortest loops are formed in apical loops rather than along the DNA plectonemic superhelix.

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.