Sliding of Proteins Non-specifically Bound to DNA: Brownian Dynamics Studies with Coarse-Grained Protein and DNA Models

DNA binding proteins efficiently search for their cognitive sites on long genomic DNA by combining 3D diffusion and 1D diffusion (sliding) along the DNA. Recent experimental results and theoretical analyses revealed that the proteins show a rotation-coupled sliding along DNA helical pitch. Here, we performed Brownian dynamics simulations using newly developed coarse-grained protein and DNA models for evaluating how hydrodynamic interactions between the protein and DNA molecules, binding affinity of the protein to DNA, and DNA fluctuations affect the one dimensional diffusion of the protein on the DNA. Our results indicate that intermolecular hydrodynamic interactions reduce 1D diffusivity by 30%. On the other hand, structural fluctuations of DNA give rise to steric collisions between the CG-proteins and DNA, resulting in faster 1D sliding of the protein. Proteins with low binding affinities consistent with experimental estimates of non-specific DNA binding show hopping along the CG-DNA. This hopping significantly increases sliding speed. These simulation studies provide additional insights into the mechanism of how DNA binding proteins find their target sites on the genome.     

Association kinetics with coupled diffusional flows: Special application to the lac repressor-operator system

The time development of the association of the lac repressor to the operator is considered in a model where the repressor is allowed to bind unspecifically to DNA and move along the DNA chain in a one-dimensional diffusion. The coupling to the three-dimensional diffusion outside the chain is introduced by letting the repressor associate and dissociate from the chain until it is finally bound to the operator. All distance correlations along the chain are included. The mean time of association is calculated and through a comparison with experimental data the molecular parameters are determined. The one-dimensional diffusion constant is found to lie of the order of 10^?9cm² s^?1. The model is sufficiently general to tie applicable to other similar systems.     

Kinetics and mechanism in the reaction of gene regulatory proteins with DNA

We have measured the kinetic properties of the Escherichia coli cAMP receptor protein (CAP) and lac repressor interacting with lac promoter restriction fragments. Under our reaction conditions (10 mM-Tris X HCl (pH 8.0 at 21 degrees C), 1 mM-EDTA, 10 microM-cAMP, 50 micrograms bovine serum albumin/ml, 5% glycerol), the association of CAP is at least a two-step process, with an initial, unstable complex formed with rate constant kappa a = 5(+/- 2.5) X 10(7) M-1 s-1. Subsequent formation of a stable complex occurs with an apparent bimolecular rate constant kappa a = 6.7 X 10(6) M-1 s-1. At low total DNA concentration, the dissociation rate constant for the specific CAP-DNA complex is 1.2 X 10(-4) s-1. The ratio of formation and dissociation rate constants yields an estimate of the equilibrium constant, Keq = 5 X 10(10) M-1, in good agreement with static results. We observed that the dissociation rate constant of both CAP-DNA and repressor-DNA complexes is increased by adding non-specific "catalytic" DNA to the reaction mixture. CAP dissociation by the concentration-dependent pathway is second-order in added non-specific DNA, consistent with either the simultaneous or the sequential participation of two DNA molecules in the reaction mechanism. The results imply a role for distal DNA in assembly-disassembly of specific CAP-DNA complexes, and are consistent with a model in which the subunits in the CAP dimer separate in the assembly-disassembly process. The dissociation of lac repressor-operator complexes was found to be DNA concentration-dependent as well, although in contrast to CAP, the reaction is first-order in catalytic DNA. Added excess operator-rich DNA gave more rapid dissociation than equivalent concentrations of non-specific DNA, indicating that the sequence content of the competing DNA influences the rate of repressor dissociation. The simplest interpretation of these observations is that lac repressor can be transferred directly from one DNA molecule to another. A comparison of the translocation rates calculated for direct transfer with those predicted by the one-dimensional sliding model indicates that direct transfer may play a role in the binding site search of lac repressor.     

Microscopic mechanism of DNA damage searching by hOGG1

The DNA backbone is often considered a track that allows long-range sliding of DNA repair enzymes in their search for rare damage sites in DNA. A proposed exemplar of DNA sliding is human 8-oxoguanine (^oG) DNA glycosylase 1 (hOGG1), which repairs mutagenic ^oG lesions in DNA. Here we use our high-resolution molecular clock method to show that macroscopic 1D DNA sliding of hOGG1 occurs by microscopic 2D and 3D steps that masquerade as sliding in resolution-limited single-molecule images. Strand sliding was limited to distances shorter than seven phosphate linkages because attaching a covalent chemical road block to a single DNA phosphate located between two closely spaced damage sites had little effect on transfers. The microscopic parameters describing the DNA search of hOGG1 were derived from numerical simulations constrained by the experimental data. These findings support a general mechanism where DNA glycosylases use highly dynamic multidimensional diffusion paths to scan DNA.     

Local and translational dynamics in DNA-lipid assemblies monitored by solid-state and diffusion NMR

The influence of electrostatic interactions on the dynamic properties of complexes containing DNA and mixtures of cationic- (DDA) and zwitterionic (DLPC) lipids are studied by means of NMR. The systems are arranged in lamellar membrane stacks intercalated by DNA molecules. This is confirmed by ³¹P-NMR, where a superposition of an axially symmetric powder pattern arising from the phospholipid membrane and an asymmetric tensor due to DNA can be fitted to the experimentally observed lineshape. The local mobility and order is assessed using two solid-state NMR techniques applicable to samples with natural isotopic abundance: WIdeline SEparation (WISE) and Separated Local Field (SLF) spectroscopy. Both experiments yield highly resolved ¹³C spectra in the direct dimension. The indirect dimension contains information about molecular dynamics through the ¹H dipolar linewidth (WISE) or the ¹H-¹³C dipolar coupling constant (SLF). The experiments suggest that DNA is static while it induces an increased disorder in the hydrocarbon chains as compared to the parent lipid case. DDA chain order is more affected than DLPC due to the attractive electrostatic interaction between DNA and the cationic lipid. Translational dynamics of the lipids and the water was measured with the Pulsed Field Gradient STimulated Echo (PFG STE) technique. The influence of lamellar domain size and the angular dependence of the diffusion coefficients and nuclear relaxation times on the results of the PFG STE experiments are discussed. The local water diffusion coefficient is reduced by a factor four from the value of bulk water, and increases as the DLPC content is increased. We observe two lipid components with an order of magnitude difference in diffusion coefficients in the DNA:DDA:DLPC precipitate and these are assigned to DLPC (fast) and DDA (slow). Cationic lipid (DDA) diffusion is decreasing a factor of 2 when DLPC is added to the pure DNA:DDA system, indicating DNA-induced lipid segregation within the bilayer and the transition from locally 2D to 1D diffusion of the DDA. The results show that DNA-lipid electrostatic interactions reduce the long-range lipid mobility but locally enhance the hydrocarbon chain dynamics by perturbing the preferred lipid packing.     

High Free-Energy Barrier of 1D Diffusion Along DNA by Architectural DNA-Binding Proteins

Architectural DNA-binding proteins function to regulate diverse DNA reactions and have the defining property of significantly changing DNA conformation. Although the 1D movement along DNA by other types of DNA-binding proteins has been visualized, the mobility of architectural DNA-binding proteins on DNA remains unknown. Here, we applied single-molecule fluorescence imaging on arrays of extended DNA molecules to probe the binding dynamics of three structurally distinct architectural DNA-binding proteins: Nhp6A, HU, and Fis. Each of these proteins was observed to move along DNA, and the salt concentration independence of the 1D diffusion implies sliding with continuous contact to DNA. Nhp6A and HU exhibit a single sliding mode, whereas Fis exhibits two sliding modes. Based on comparison of the diffusion coefficients and sizes of many DNA binding proteins, the architectural proteins are categorized into a new group distinguished by an unusually high free-energy barrier for 1D diffusion. The higher free-energy barrier for 1D diffusion by architectural proteins can be attributed to the large DNA conformational changes that accompany binding and impede rotation-coupled movement along the DNA grooves.     

Improved NMR spectra of a protein–DNA complex through rational mutagenesis and the application of a sensitivity optimized isotope-filtered NOESY experiment

The NMR spectra of the complex between the DNA-binding domain of the Dead ringer protein (DRI-DBD, Gly262-Gly398) and its DNA binding site (DRI-DBD:DNA, 26 kDa) have been optimized by biochemical and spectroscopic means. First, we demonstrate the utility of a modified 2D [F1,F2] ¹³C-filtered NOESY experiment that employs a ¹J_HC versus chemical shift optimized adiabatic ¹³C inversion pulse [Zwahlen, C. et al. (1997) J. Am. Chem. Soc., 119, 6711–6721]. The new sequence is shown to be more sensitive than previously published pulse schemes (up to 40% in favorable cases) and its utility is demonstrated using two protein–DNA complexes. Second, we demonstrate that the targeted replacement of an interfacial aromatic residue in the DRI-DBD:DNA complex substantially reduces line broadening within its NMR spectra. The spectral changes are dramatic, salvaging a protein–DNA complex that was originally ill suited for structural analysis by NMR. This biochemical approach is not a general method, but may prove useful in the spectral optimization of other protein complexes that suffer from interfacial line broadening caused by dynamic changes in proximal aromatic rings.     

Diffusion of the Restriction Nuclease EcoRI along DNA

Many specific sequence DNA binding proteins locate their target sequence by first binding to DNA nonspecifically, then by linearly diffusing or hopping along DNA until either the protein dissociates from the DNA or it finds the recognition sequence. We have devised a method for measuring one-dimensional diffusion along DNA based on the ratio of the dissociation rate of protein from DNA fragments containing one specific binding site to the dissociation rate from DNA fragments containing two specific binding sites. Our extensive measurements of dissociation rates and specific-nonspecific relative binding constants of the restriction nuclease EcoRI enable us to determine the diffusion rate of nonspecifically bound protein along the DNA. By varying the distance between the two binding sites, we confirm a linear diffusion mechanism. The sliding rate is relatively insensitive to salt concentration and osmotic pressure, indicating that the protein moves smoothly along the DNA probably following the helical phosphate-sugar backbone of DNA. We calculate a diffusion coefficient for EcoRI of 3 × 10⁴ bp² s^− 1 EcoRI is able to diffuse ∼ 150 bp, on average, along the DNA in 1 s. This diffusion rate is about 2000-fold slower than the diffusion of free protein in solution. A factor of 40-50 can be accounted for by rotational friction resulting from following the helical path of the DNA backbone. Two possibilities could account for the remaining activation energy: salt bridges between the DNA and the protein are transiently broken, or the water structure at the protein-DNA interface is disrupted as the two surfaces move past each other.     

Lac repressor-operator interaction: DNA length dependence

The interaction of the E. coli lac operon repressor with its operator DNA has been directly examined as a function of the length of operator-containing DNA. The apparent bimolecular association rate constants were calculated as ka = (kd/KD), where the dissociation equilibrium constant, KD and the dissociation rate constant, kd, were measured by nitrocellulose filter adsorption assays. The values obtained for the overall association rate constants are compared with theoretical association rate curves for specific mechanisms. Association of the repressor with short operator containing DNA fragments (less than 70 base pairs) occurs at rates expected of three-dimensional diffusion. Our data also imply that at longer DNA lengths a combination of three-dimensional diffusion with one-dimensional sliding along with hopping and/or intersegment transfer must be involved to facilitate the repressor operator association.     

Sliding and jumping of single EcoRV restriction enzymes on non-cognate DNA

The restriction endonuclease EcoRV can rapidly locate a short recognition site within long non-cognate DNA using 'facilitated diffusion'. This process has long been attributed to a sliding mechanism, in which the enzyme first binds to the DNA via nonspecific interaction and then moves along the DNA by 1D diffusion. Recent studies, however, provided evidence that 3D translocations (hopping/jumping) also help EcoRV to locate its target site. Here we report the first direct observation of sliding and jumping of individual EcoRV molecules along nonspecific DNA. Using fluorescence microscopy, we could distinguish between a slow 1D diffusion of the enzyme and a fast translocation mechanism that was demonstrated to stem from 3D jumps. Salt effects on both sliding and jumping were investigated, and we developed numerical simulations to account for both the jump frequency and the jump length distribution. We deduced from our study the 1D diffusion coefficient of EcoRV, and we estimated the number of jumps occurring during an interaction event with nonspecific DNA. Our results substantiate that sliding alternates with hopping/jumping during the facilitated diffusion of EcoRV and, furthermore, set up a framework for the investigation of target site location by other DNA-binding proteins.     

How do site-specific DNA-binding proteins find their targets?

Essentially all the biological functions of DNA depend on site-specific DNA-binding proteins finding their targets, and therefore ‘searching’ through megabases of non-target DNA. In this article, we review current understanding of how this sequence searching is done. We review how simple diffusion through solution may be unable to account for the rapid rates of association observed in experiments on some model systems, primarily the Lac repressor. We then present a simplified version of the ‘facilitated diffusion’ model of Berg, Winter and von Hippel, showing how non-specific DNA–protein interactions may account for accelerated targeting, by permitting the protein to sample many binding sites per DNA encounter. We discuss the 1-dimensional ‘sliding’ motion of protein along non-specific DNA, often proposed to be the mechanism of this multiple site sampling, and we discuss the role of short-range diffusive ‘hopping’ motions. We then derive the optimal range of sliding for a few physical situations, including simple models of chromosomes in vivo, showing that a sliding range of ~100 bp before dissociation optimizes targeting in vivo. Going beyond first-order binding kinetics, we discuss how processivity, the interaction of a protein with two or more targets on the same DNA, can reveal the extent of sliding and we review recent experiments studying processivity using the restriction enzyme EcoRV. Finally, we discuss how single molecule techniques might be used to study the dynamics of DNA site-specific targeting of proteins.     

Unique 31P Spectral Response to the Formation of a Specific Restriction Enzyme–DNA Complex

Protein-induced distortion is a dramatic but not universally observed feature of sequence-specific DNA interactions. This is illustrated by the crystal structures of restriction enzyme–DNA complexes: While some of these structures exhibit DNA distortion, others do not. Among the latter is PvuII endonuclease, a small enzyme that is also amenable to NMR spectroscopic studies. Here ³¹P NMR spectroscopy is applied to demonstrate the unique spectral response of DNA to sequence-specific protein interactions. The ³¹P NMR spectrum of a noncognate DNA exhibits only spectral broadening upon the addition of enzyme. However, when enzyme is added to target DNA, a number of ³¹P resonances shift dramatically. The magnitudes of the chemical shifts (2–3 ppm) are among the largest observed. Site-specific substitution with phosphoramidates and phosphorothioates are used analyze these effects. While such spectral features have been interpreted as indicative of DNA backbone distortions, FRET analysis indicates that this does not occur in PvuII-cognate DNA complexes in solution. The distinct ³¹P spectral signature observed for cognate DNA mirrors that observed for the enzyme, underscoring the unique features of cognate complex formation.     

The rotation-coupled sliding of EcoRV

It has been proposed that certain type II restriction enzymes (REs), such as EcoRV, track the helical pitch of DNA as they diffuse along DNA, a so-called rotation-coupled sliding. As of yet, there is no direct experimental observation of this phenomenon, but mounting indirect evidence gained from single-molecule imaging of RE-DNA complexes support the hypothesis. We address this issue by conjugating fluorescent labels of varying size (organic dyes, proteins and quantum dots) to EcoRV, and by fusing it to the engineered Rop protein scRM6. Single-molecule imaging of these modified EcoRVs sliding along DNA provides us with their linear diffusion constant (D(1)), revealing a significant size dependency. To account for the dependence of D(1) on the size of the EcoRV label, we have developed four theoretical models describing different types of motion along DNA and find that our experimental results are best described by rotation-coupled sliding of the protein. The similarity of EcoRV to other type II REs and DNA binding proteins suggests that this type of motion could be widely preserved in other biological contexts.     

Search by proteins for their DNA target site: 1. The effect of DNA conformation on protein sliding

The recognition of DNA-binding proteins (DBPs) to their specific site often precedes by a search technique in which proteins slide, hop along the DNA contour or perform inter-segment transfer and 3D diffusion to dissociate and re-associate to distant DNA sites. In this study, we demonstrated that the strength and nature of the non-specific electrostatic interactions, which govern the search dynamics of DBPs, are strongly correlated with the conformation of the DNA. We tuned two structural parameters, namely curvature and the extent of helical twisting in circular DNA. These two factors are mutually independent of each other and can modulate the electrostatic potential through changing the geometry of the circular DNA conformation. The search dynamics for DBPs on circular DNA is therefore markedly different compared with linear B-DNA. Our results suggest that, for a given DBP, the rotation-coupled sliding dynamics is precluded in highly curved DNA (as well as for over-twisted DNA) because of the large electrostatic energy barrier between the inside and outside of the DNA molecule. Under such circumstances, proteins prefer to hop in order to explore interior DNA sites. The change in the balance between sliding and hopping propensities as a function of DNA curvature or twisting may result in different search efficiency and speed.     

The Role of Noncognate Sites in the 1D Search Mechanism of EcoRI

A one-dimensional (1D) search is an essential step in DNA target recognition. Theoretical studies have suggested that the sequence dependence of 1D diffusion can help resolve the competing demands of a fast search and high target affinity, a conflict known as the speed-selectivity paradox. The resolution requires that the diffusion energy landscape is correlated with the underlying specific binding energies. In this work, we report observations of a 1D search by quantum dot-labeled EcoRI. Our data supports the view that proteins search DNA via rotation-coupled sliding over a corrugated energy landscape. We observed that whereas EcoRI primarily slides along DNA at low salt concentrations, at higher concentrations, its diffusion is a combination of sliding and hopping. We also observed long-lived pauses at genomic star sites, which differ by a single nucleotide from the target sequence. To reconcile these observations with prior biochemical and structural data, we propose a model of search in which the protein slides over a sequence-independent energy landscape during fast search but rapidly interconverts with a “hemispecific” binding mode in which a half site is probed. This half site interaction stabilizes the transition to a fully specific mode of binding, which can then lead to target recognition.     

Motion of a DNA sliding clamp observed by single molecule fluorescence spectroscopy

DNA sliding clamps attach to polymerases and slide along DNA to allow rapid, processive replication of DNA. These clamps contain many positively charged residues that could curtail the sliding due to attractive interactions with the negatively charged DNA. By single-molecule spectroscopy we have observed a fluorescently labeled sliding clamp (polymerase III beta subunit or beta clamp) loaded onto freely diffusing, single-stranded M13 circular DNA annealed with fluorescently labeled DNA oligomers of up to 90 bases. We find that the diffusion constant for the beta clamp diffusing along DNA is on the order of 10(-14) m(2)/s, at least 3 orders of magnitude less than that for diffusion through water alone. We also find evidence that the beta clamp remains at the 3' end in the presence of Escherichia coli single-stranded-binding protein. These results may imply that the clamp not only acts to hold the polymerase on the DNA but also prevents excessive drifting along the DNA.