Transcription factors (TFs) efficiently locate their target DNA sequences by combining three-dimensional diffusion and one-dimensional sliding on nonspecific DNA. To balance rapid sliding with strong specific binding, TFs were proposed to switch between search and recognition conformations. For Escherichia coli lac repressor (LacI), the folding of the hinge helices has been implicated in the conformational switch. Here, we tested how mutations in the hinge region impact the search speed and binding stability. Based on molecular dynamics simulations, we selected two LacI mutants favoring either search or recognition conformation. We measured the binding kinetics of the mutants both in vitro on DNA microarrays with 2479 different Lac operators and in vivo via single-molecule experiments. We identified a mutation that enhances the specificity but reduces binding strength globally, and another mutation that makes the operator binding stronger but also reduces the specificity. However, the altered specificity impacts the search time less than expected. Instead, the major effect was impaired dissociation in response to Isopropyl β-D-1-thiogalactopyranoside (IPTG) induction for the strongly binding mutant. Together with earlier reports of affinity-inducibility trade-offs in LacI, our data support the model in which the trade-off is between binding stability and inducibility rather than between speed and binding stability.

The molecular details governing transcription factor (TF) binding and the formation of accessible chromatin are not yet quantitatively understood—including how sequence context modulates affinity, how TFs search DNA, the kinetics of TF occupancy, and how motif grammars coordinate binding. To resolve these questions for a human TF, erythroid Krüppel-like factor (eKLF/KLF1), we quantitatively compare, in high throughput, in vitro TF binding rates and affinities with in vivo single-molecule TF and nucleosome occupancies and in vivo-derived deep learning models. We find that 40-fold flanking sequence effects on affinity are consistent with distal flanks tuning TF search parameters and captured by a linear energy model. Motif recognition probability, rather than time in the bound state, drives affinity changes, and in vitro and in nuclei measurements exhibit consistent, minutes-long TF residence times. Finally, in vitro biophysical parameters predict in vivo sequence preferences and single-molecule chromatin states for unseen motif grammars.

DNA folding thermodynamics are central to many biological processes and biotechnological applications involving base-pairing. Current methods for predicting stability from DNA sequence use nearest-neighbor models that struggle to accurately capture the diverse sequence dependence of secondary structural motifs beyond Watson-Crick base pairs, likely due to insufficient experimental data. In this work, we introduce a massively parallel method, Array Melt, that uses fluorescence-based quenching signals to measure the equilibrium stability of millions of DNA hairpins simultaneously on a repurposed Illumina sequencing flow cell. By leveraging this dataset of 27,732 sequences with two-state melting behaviors, we derive a NUPACK-compatible model (dna24), a rich parameter model that exhibits higher accuracy, and a graph neural network (GNN) model that identifies relevant interactions within DNA beyond nearest neighbors. All models show improved accuracy in predicting DNA folding thermodynamics, enabling more effective in silico design of qPCR primers, oligo hybridization probes, and DNA origami.

Transcription factors (TFs) efficiently locate their target DNA sequences by combining three-dimensional diffusion and one-dimensional sliding on nonspecific DNA. To balance rapid sliding with strong specific binding, TFs were proposed to switch between search and recognition conformations. For E. coli lac repressor (LacI), the folding of the hinge helices has been implicated in the conformational switch. Here, we tested how mutations in the hinge region impact the search speed and binding stability. Based on molecular dynamics simulations, we selected two LacI mutants favoring either search or recognition conformation. We measured the binding kinetics of the mutants both in vitro on DNA microarrays with 2,479 different Lac operators and in vivo via single-molecule experiments. We conclude that a hinge region mutation causing less helix propensity enhances the specificity but reduces binding strength globally, while a hinge region mutation causing higher helix propensity has opposite effects. However, altered specificity impacts the search time less than expected. Instead, the major effect was impaired dissociation in response to IPTG induction for the strongly binding mutant. Together with earlier reports of affinity-inducibility trade-offs in LacI, our data support the model in which the hinge region governs a trade-off between binding stability and inducibility rather than between speed and binding stability.

The binding of multiple transcription factors (TFs) to genomic enhancers activates gene expression in mammalian cells. However, the molecular details that link enhancer sequence to TF binding, promoter state, and gene expression levels remain opaque. We applied single-molecule footprinting (SMF) to measure the simultaneous occupancy of TFs, nucleosomes, and components of the transcription machinery on engineered enhancer/promoter constructs with variable numbers of TF binding sites for both a synthetic and an endogenous TF. We find that activation domains enhance a TF’s capacity to compete with nucleosomes for binding to DNA in a BAF-dependent manner, TF binding on nucleosome-free DNA is consistent with independent binding between TFs, and average TF occupancy linearly contributes to promoter activation rates. We also decompose TF strength into separable binding and activation terms, which can be tuned and perturbed independently. Finally, we develop thermodynamic and kinetic models that quantitatively predict both the binding microstates observed at the enhancer and subsequent time-dependent gene expression. This work provides a template for quantitative dissection of distinct contributors to gene activation, including the activity of chromatin remodelers, TF activation domains, chromatin acetylation, TF concentration, TF binding affinity, and TF binding site configuration.

The binding of multiple transcription factors (TFs) to genomic enhancers drives gene expression in mammalian cells1. However, the molecular details that link enhancer sequence to TF binding, promoter state and transcription levels remain unclear. Here we applied single-molecule footprinting2,3 to measure the simultaneous occupancy of TFs, nucleosomes and other regulatory proteins on engineered enhancer-promoter constructs with variable numbers of TF binding sites for both a synthetic TF and an endogenous TF involved in the type I interferon response. Although TF binding events on nucleosome-free DNA are independent, activation domains recruit cofactors that destabilize nucleosomes, driving observed TF binding cooperativity. Average TF occupancy linearly determines promoter activity, and we decompose TF strength into separable binding and activation terms. Finally, we develop thermodynamic and kinetic models that quantitatively predict both the enhancer binding microstates and gene expression dynamics. This work provides a template for the quantitative dissection of distinct contributors to gene expression, including TF activation domains, concentration, binding affinity, binding site configuration and recruitment of chromatin regulators.

Short tandem repeats (STRs) are enriched in eukaryotic cis-regulatory elements and alter gene expression, yet how they regulate transcription remains unknown. We found that STRs modulate transcription factor (TF)-DNA affinities and apparent on-rates by about 70-fold by directly binding TF DNA-binding domains, with energetic impacts exceeding many consensus motif mutations. STRs maximize the number of weakly preferred microstates near target sites, thereby increasing TF density, with impacts well predicted by statistical mechanics. Confirming that STRs also affect TF binding in cells, neural networks trained only on in vivo occupancies predicted effects identical to those observed in vitro. Approximately 90% of TFs preferentially bound STRs that need not resemble known motifs, providing a cis-regulatory mechanism to target TFs to genomic sites.

RNAs are central to fundamental biological processes in all known organisms. The set of possible intramolecular interactions of RNA nucleotides defines the range of alternative structural conformations of a specific RNA that can coexist, and these structures enable functional catalytic properties of RNAs and/or their productive intermolecular interactions with other RNAs or proteins. However, the immense combinatorial space of potential RNA sequences has precluded predictive mapping between RNA sequence and molecular structure and function. Recent advances in high-throughput approaches in vitro have enabled quantitative thermodynamic and kinetic measurements of RNA-RNA and RNA-protein interactions, across hundreds of thousands of sequence variations. In this Review, we explore these techniques, how they can be used to understand RNA function and how they might form the foundations of an accurate model to predict the structure and function of an RNA directly from its nucleotide sequence. The experimental techniques and modelling frameworks discussed here are also highly relevant for the sampling of sequence-structure-function space of DNAs and proteins.

Sequence-specific binding of proteins to DNA is essential for accessing genetic information. Here, we derive a simple equation for target-site recognition, which uncovers a previously unrecognized coupling between the macroscopic association and dissociation rates of the searching protein. Importantly, this relationship makes it possible to recover the relevant microscopic rates from experimentally determined macroscopic ones. We directly test the equation by observing the binding and unbinding of individual lac repressor (LacI) molecules during target search. We find that LacI dissociates from different target sequences with essentially identical microscopic dissociation rates. Instead, sequence specificity is determined by the efficiency with which the protein recognizes different targets, effectively reducing its risk of being retained on a non-target sequence. Our theoretical framework also accounts for the coupling between off-target binding and unbinding of the catalytically inactive Cas9 (dCas9), showing that the binding pathway can be obtained from macroscopic data. One Sentence Summary Association and dissociation rates are anti-correlated for reactions that include a nonspecific probing step.

Sequence-specific binding of proteins to DNA is essential for accessing genetic information. We derive a model that predicts an anticorrelation between the macroscopic association and dissociation rates of DNA binding proteins. We tested the model for thousands of different lac operator sequences with a protein binding microarray and by observing kinetics for individual lac repressor molecules in single-molecule experiments. We found that sequence specificity is mainly governed by the efficiency with which the protein recognizes different targets. The variation in probability of recognizing different targets is at least 1.7 times as large as the variation in microscopic dissociation rates. Modulating the rate of binding instead of the rate of dissociation effectively reduces the risk of the protein being retained on nontarget sequences while searching.

Many proteins that bind specific DNA sequences search the genome by combining three-dimensional diffusion with one-dimensional sliding on nonspecific DNA1–5. Here we combine resonance energy transfer and fluorescence correlation measurements to characterize how individual lac repressor (LacI) molecules explore the DNA surface during the one-dimensional phase of target search. To track the rotation of sliding LacI molecules on the microsecond timescale, we use real-time single-molecule confocal laser tracking combined with fluorescence correlation spectroscopy (SMCT–FCS). The fluctuations in fluorescence signal are accurately described by rotation-coupled sliding, in which LacI traverses about 40 base pairs (bp) per revolution. This distance substantially exceeds the 10.5-bp helical pitch of DNA; this suggests that the sliding protein frequently hops out of the DNA groove, which would result in the frequent bypassing of target sequences. We directly observe such bypassing using single-molecule fluorescence resonance energy transfer (smFRET). A combined analysis of the smFRET and SMCT–FCS data shows that LacI hops one or two grooves (10–20 bp) every 200–700 μs. Our data suggest a trade-off between speed and accuracy during sliding: the weak nature of nonspecific protein–DNA interactions underlies operator bypassing, but also speeds up sliding. We anticipate that SMCT–FCS, which monitors rotational diffusion on the microsecond timescale while tracking individual molecules with millisecond resolution, will be applicable to the real-time investigation of many other biological interactions and will effectively extend the accessible time regime for observing these interactions by two orders of magnitude. Single-molecule fluorescence resonance energy transfer and real-time confocal laser tracking with fluorescence correlation spectroscopy together characterize how individual lac repressor molecules bypass operator sites while exploring the DNA surface at microsecond timescales.

Many proteins that bind specific DNA sequences search the genome by combining three dimensional (3D) diffusion in the cytoplasm with one dimensional (1D) sliding on non-specific regions of the DNA. It is however not known how sliding proteins are oriented with respect to DNA in order to recognize specific sequences. Here we measure the polarization of fluorescence emission from single fluorescently labeled lac repressor (LacI) molecules sliding on stretched DNA. Real-time feedback-coupled confocal single-particle tracking allows us to measure fluorescence correlation of the sliding molecules. We find that the fluctuations in the fluorescence signal on the μs timescale are accurately described by rotation-coupled sliding on DNA. On average, LacI moves ~50 base pairs per revolution, which is significantly longer than the 10.5 bp helical periodicity of DNA. Our data support a facilitated diffusion model where the transcription factor (TF) scans the DNA grooves for hydrogen bonding opportunities in a pre-aligned orientation with occasional slippage out of the groove.

The lactose operon repressor protein LacI has long served as a paradigm of the bacterial transcription factors. However, the mechanisms whereby LacI rapidly locates its cognate binding site on the bacterial chromosome are still elusive. Single-molecule fluorescence imaging approaches are well suited for the study of these mechanisms but rely on a functionally compatible fluorescence labeling of LacI. Particularly attractive for protein fluorescence labeling are synthetic fluorophores due to their small size and favorable photophysical characteristics. Synthetic fluorophores are often conjugated to natively occurring cysteine residues using maleimide chemistry. For a site-specific and functionally compatible labeling with maleimide fluorophores, the target protein often needs to be redesigned to remove unwanted native cysteines and to introduce cysteines at locations better suited for fluorophore attachment. Biochemical screens can then be employed to probe for the functional activity of the redesigned protein both before and after dye labeling. Here, we report a mutagenesis-based redesign of LacI to enable a functionally compatible labeling with maleimide fluorophores. To provide an easily accessible labeling site in LacI, we introduced a single cysteine residue at position 28 in the DNA-binding headpiece of LacI and replaced two native cysteines with alanines where derivatization with bulky substituents is known to compromise the protein’s activity. We find that the redesigned LacI retains a robust activity in vitro and in vivo, provided that the third native cysteine at position 281 is retained in LacI. In a total internal reflection microscopy assay, we observed individual Cy3-labeled LacI molecules bound to immobilized DNA harboring the cognate O1 operator sequence, indicating that the dye-labeled LacI is functionally active. We have thus been able to generate a functional fluorescently labeled LacI that can be used to unravel mechanistic details of LacI target search at the single molecule level.

The polarization of the light emitted from a fluorophore is determined by its dipole moment and hence by the direction of the molecule. We are using this phenomena to measure the orientation of labeled transcription factor molecules. We are studying the lac repressor, labeled with a fluorophore that is predicted to be rigid and oriented along the DNA when the protein is associated to the DNA substrate. In our experiments we measure polarization of the light emitted from single transcription factor molecules when they are bound specifically and sliding on flow stretched DNA. The polarization data acquired is shifted towards the direction of DNA with relatively small changes in the signal over several seconds when using an integration time of ∼100 ms. This shows that the transcription factor has a preferred and maintained orientation while interacting with the DNA substrate on this timescale. For transcription factors sliding on DNA, very little polarization preference is observed in the direction that would be averaged out by rotation of the fluorophore around the DNA. Taken together these observations strongly support the model of rotational sliding for the lac repressor. By using a simplified model of dipole excitation and emission we have also been able to quantify the fluorophore orientation and wobbliness for individual transcription factor molecules sliding on DNA. This analysis shows that around 40 % of the labeled molecules have an average dipole orientation close to parallel to DNA (<20°) with a standard deviation in the orientation of 40°-60°. The remaining 60 % have a more undefined orientation with larger standard deviation (>60°).

DNA binding proteins are very efficient in recognizing/searching target sites over long stretches of non-specific DNA sequences. Current theories suggests that the search efficiency comes partly from non-specific interactions with DNA. During these DNA interactions the search dimensionality is lowered from a three-dimensional diffusion process to a one-dimensional. When this occurs the protein is for a short time sliding along the DNA contour. Indirect evidence suggest that the protein is rotating along the DNA spiral but the underlying dynamics of these sliding events are still not well understood. To experimentally investigate sliding dynamics, over time scales comparable to expected protein rotation times, we have developed a single molecule tracking polarization microscope. The platform is based on the scanning confocal principle and has capability of real time positioning feedback, polarization measurement and single photon time tagging. Photon acquisition and positioning controls are all implemented on an FPGA allowing for fast feedback without losing single photon timing information. With a reaction time close to a millisecond, for the positioning feedback, and a 200 Mhz clock for time tagging we study sliding of the lac repressor. Our lac repressors are each labeled with a fluorophore that is predicted to be in a rigid conformation. This allows us to use polarization as an indicator of the protein orientation. Using flow stretched DNA we are able to capture DNA-protein sliding events. Preliminary results reveals a polarization shift that is consistent with the expected dipole orientation and flow direction and thus support the rotational hypothesis.

Human ALC1 is an oncogene-encoded chromatin-remodeling enzyme required for DNA repair that possesses a poly(ADP-ribose) (PAR)-binding macro domain. Its engagement with PARylated PARP1 activates ALC1 at sites of DNA damage, but the underlying-mechanism remains unclear. Here, we establish a dual role for the macro domain in autoinhibition of ALC1 ATPase activity and coupling to nucleosome mobilization. In the absence of DNA damage, an inactive conformation of the ATPase is maintained by juxtaposition of the macro domain against predominantly the C-terminal ATPase lobe through conserved electrostatic interactions. Mutations within this interface displace the macro domain, constitutively activate the ALC1 ATPase independent of PARylated PARP1, and alter the dynamics of ALC1 recruitment at DNA damage sites. Upon DNA damage, binding of PARylated PARP1 by the macro domain induces a conformational change that relieves autoinhibitory interactions with the ATPase motor, which selectively activates ALC1 remodeling upon recruitment to sites of DNA damage.

The helical structure of DNA imposes constraints on the rate of diffusionlimited protein binding. Here we solve the reaction–diffusion equations for DNA-like geometries and extend with simulations when necessary. We find that the helical structure can make binding to the DNA more than twice as fast compared to a case where DNA would be reactive only along one side. We also find that this rate advantage remains when the contributions from steric constraints and rotational diffusion of the DNA-binding protein are included. Furthermore, we find that the association rate is insensitive to changes in the steric constraints on the DNA in the helix geometry, while it is much more dependent on the steric constraints on the DNA-binding protein. We conclude that the helical structure of DNA facilitates the nonspecific binding of transcription factors and structural DNA-binding proteins in general. S Online supplementary data available from stacks.iop.org/JPA/49/364002/ mmedia

The helical structure of DNA imposes constraints on the rate of diffusion-limited protein binding. Here we solve the reaction-diffusion equations for DNA-like geometries and extend with simulations when necessary. We find that the helical structure can make binding to the DNA more than twice as fast compared to a case where DNA would be reactive only along one side. We also find that this rate advantage remains when the contributions from steric constraints and rotational diffusion of the DNA-binding protein are included. Furthermore, we find that the association rate is insensitive to changes in the steric constraints on the DNA in the helix geometry, while it is much more dependent on the steric constraints on the DNA-binding protein. We conclude that the helical structure of DNA facilitates the nonspecific binding of transcription factors and structural DNA-binding proteins in general.