Target search of N sliding proteins on a DNA

Role of Shape Deformation of DNA-Binding Sites in Regulating the Efficiency and Specificity in Their Recognition by DNA-Binding Proteins

Accurate transcription of genetic information is crucial, involving precise recognition of the binding motifs by DNA-binding proteins. While some proteins rely on short-range hydrophobic and hydrogen bonding interactions at binding sites, others employ a DNA shape readout mechanism for specific recognition. In this mechanism, variations in DNA shape at the binding motif resulted from either inherent flexibility or binding of proteins at adjacent sites are sensed and capitalized by the searching proteins to locate them specifically. Through extensive computer simulations, we investigate both scenarios to uncover the underlying mechanism and origin of specificity in the DNA shape readout mechanism. Our findings reveal that deformation in shape at the binding motif creates an entropy funnel, allowing information about altered shapes to manifest as fluctuations in minor groove widths. This signal enhances the efficiency of nonspecific search of nearby proteins by directing their movement toward the binding site, primarily driven by a gain in entropy. We propose this as a generic mechanism for DNA shape readout, where specificity arises from the alignment between the molecular frustration of the searching protein and the ruggedness of the entropic funnel governed by molecular features of the protein and arrangement of the DNA bases at the binding site, respectively.

Dynamic interplay between target search and recognition for a Type I CRISPR-Cas system

CRISPR-Cas effector complexes enable the defense against foreign nucleic acids and have recently been exploited as molecular tools for precise genome editing at a target locus. To bind and cleave their target, the CRISPR-Cas effectors have to interrogate the entire genome for the presence of a matching sequence. Here we dissect the target search and recognition process of the Type I CRISPR-Cas complex Cascade by simultaneously monitoring DNA binding and R-loop formation by the complex. We directly quantify the effect of DNA supercoiling on the target recognition probability and demonstrate that Cascade uses facilitated diffusion for its target search. We show that target search and target recognition are tightly linked and that DNA supercoiling and limited 1D diffusion need to be considered when understanding target recognition and target search by CRISPR-Cas enzymes and engineering more efficient and precise variants.

Computational investigations on target-site searching and recognition mechanisms by thymine DNA glycosylase during DNA repair process

DNA glycosylase, as one member of DNA repair machineries, plays an essential role in correcting mismatched/damaged DNA nucleotides by cleaving the N-glycosidic bond between the sugar and target nucleobase through the base excision repair (BER) pathways. Efficient corrections of these DNA lesions are critical for maintaining genome integrity and preventing premature aging and cancers. The target-site searching/recognition mechanisms and the subsequent conformational dynamics of DNA glycosylase, however, remain challenging to be characterized using experimental techniques. In this review, we summarize our recent studies of sequential structural changes of thymine DNA glycosylase (TDG) during the DNA repair process, achieved mostly by molecular dynamics (MD) simulations. Computational simulations allow us to reveal atomic-level structural dynamics of TDG as it approaches the target-site, and pinpoint the key structural elements responsible for regulating the translocation of TDG along DNA. Subsequently, upon locating the lesions, TDG adopts a base-flipping mechanism to extrude the mispaired nucleobase into the enzyme active-site. The constructed kinetic network model elucidates six metastable states during the base-extrusion process and suggests an active role of TDG in flipping the intrahelical nucleobase. Finally, the molecular mechanism of product release dynamics after catalysis is also summarized. Taken together, we highlight to what extent the computational simulations advance our knowledge and understanding of the molecular mechanism underlying the conformational dynamics of TDG, as well as the limitations of current theoretical work.

Time averaging and emerging nonergodicity upon resetting of fractional Brownian motion and heterogeneous diffusion processes

How different are the results of constant-rate resetting of anomalous-diffusion processes in terms of their ensemble-averaged versus time-averaged mean-squared displacements (MSDs versus TAMSDs) and how does stochastic resetting impact nonergodicity? We examine, both analytically and by simulations, the implications of resetting on the MSD- and TAMSD-based spreading dynamics of particles executing fractional Brownian motion (FBM) with a long-time memory, heterogeneous diffusion processes (HDPs) with a power-law space-dependent diffusivity D(x)=D_{0}|x|^{γ} and their "combined" process of HDP-FBM. We find, inter alia, that the resetting dynamics of originally ergodic FBM for superdiffusive Hurst exponents develops disparities in scaling and magnitudes of the MSDs and mean TAMSDs indicating weak ergodicity breaking. For subdiffusive HDPs we also quantify the nonequivalence of the MSD and TAMSD and observe a new trimodal form of the probability density function. For reset FBM, HDPs and HDP-FBM we compute analytically and verify by simulations the short-time MSD and TAMSD asymptotes and long-time plateaus reminiscent of those for processes under confinement. We show that certain characteristics of these reset processes are functionally similar despite a different stochastic nature of their nonreset variants. Importantly, we discover nonmonotonicity of the ergodicity-breaking parameter EB as a function of the resetting rate r. For all reset processes studied we unveil a pronounced resetting-induced nonergodicity with a maximum of EB at intermediate r and EB∼(1/r)-decay at large r. Alongside the emerging MSD-versus-TAMSD disparity, this r-dependence of EB can be an experimentally testable prediction. We conclude by discussing some implications to experimental systems featuring resetting dynamics.

Atomic resolution of short-range sliding dynamics of thymine DNA glycosylase along DNA minor-groove for lesion recognition

Thymine DNA glycosylase (TDG), as a repair enzyme, plays essential roles in maintaining the genome integrity by correcting several mismatched/damaged nucleobases. TDG acquires an efficient strategy to search for the lesions among a vast number of cognate base pairs. Currently, atomic-level details of how TDG translocates along DNA as it approaches the lesion site and the molecular mechanisms of the interplay between TDG and DNA are still elusive. Here, by constructing the Markov state model based on hundreds of molecular dynamics simulations with an integrated simulation time of ∼25 μs, we reveal the rotation-coupled sliding dynamics of TDG along a 9 bp DNA segment containing one G·T mispair. We find that TDG translocates along DNA at a relatively faster rate when distant from the lesion site, but slows down as it approaches the target, accompanied by deeply penetrating into the minor-groove, opening up the mismatched base pair and significantly sculpturing the DNA shape. Moreover, the electrostatic interactions between TDG and DNA are found to be critical for mediating the TDG translocation. Notably, several uncharacterized TDG residues are identified to take part in regulating the conformational switches of TDG occurred in the site-transfer process, which warrants further experimental validations.

Exact solution to the first-passage problem for a particle with a dichotomous diffusion coefficient

We consider the problem of first-passage time for reaching a boundary of a particle which diffuses in one dimension and is confined to the region x∈(0,L), with a diffusion coefficient that switches randomly between two states, having diffusivities that are different. Exact analytical expressions are found for the survival probability of the particle as a function of time. The survival probability has a multiexponential decay, and to characterize it, we use the average rate constant k, as well as the instantaneous rate r(t). Our approach can easily be extended to the case where the diffusion coefficient takes n different values. The model should be of interest to biological processes, in which a reactant searches for a target in a heterogeneous environment, making the diffusion coefficient a random function of time. The best example for this is a protein searching for a target site on the DNA.

Motion of a Brownian particle in the presence of reactive boundaries

We study the one-dimensional motion of a Brownian particle inside a confinement described by two reactive boundaries which can partially reflect or absorb the particle. Understanding the effects of such boundaries is important in physics, chemistry, and biology. We compute the probability density of the particle displacement exactly, from which we derive expressions for the survival probability and the mean absorption time as a function of the reactive coefficients. Furthermore, using the Feynman-Kac formalism, we investigate the local time profile, which is the fluctuating time spent by the particle at a given location, both till a fixed observation time and till the absorption time. Our analytical results are compared to numerical simulations, showing perfect agreement.

Quantifying the two-state facilitated diffusion model of protein-DNA interactions

The current report extends the facilitated diffusion model to account for conflict between the search and recognition binding modes adopted by DNA-binding proteins (DBPs) as they search DNA and subsequently recognize and bind to their specific binding site. The speed of the search dynamics is governed by the energetic ruggedness of the protein-DNA landscape, whereas the rate for the recognition process is mostly dictated by the free energy barrier for the transition between the DBP's search and recognition binding modes. We show that these two modes are negatively coupled, such that fast 1D sliding and rapid target site recognition probabilities are unlikely to coexist. Thus, a tradeoff occurs between optimizing the timescales for finding and binding the target site. We find that these two kinetic properties can be balanced to produce a fast timescale for the total target search and recognition process by optimizing frustration. Quantification of the facilitated diffusion model by including a frustration term enables it to explain several experimental observations concerning search and recognition speeds. The extended model captures experimental estimate of the energetic ruggedness of the protein-DNA landscape and predicts how various molecular properties of protein-DNA binding affect recognition kinetics. Particularly, point mutations may change the frustration and so affect protein association with DNA, thus providing a means to modulate protein-DNA affinity by manipulating the protein's association or dissociation reactions.

Acceleration of bursty multiprotein target search kinetics on DNA by colocalisation

Proteins are capable of locating specific targets on DNA by employing a facilitated diffusion process with intermittent 1D and 3D search steps. We here uncover the implications of colocalisation of protein production and DNA binding sites via computer simulations.

Bulk-mediated surface transport in the presence of bias

Surface transport, when the particle is allowed to leave the surface, travel in the bulk for some time, and then return to the surface, is referred to as bulk-mediated surface transport. Recently, we proposed a formalism that significantly simplifies analysis of bulk-mediated surface diffusion [A. M. Berezhkovskii, L. Dagdug, and S. M. Bezrukov, J. Chem. Phys. 143, 084103 (2015)]. Here this formalism is extended to bulk-mediated surface transport in the presence of bias, i.e., when the particle has arbitrary drift velocities on the surface and in the bulk. A key advantage of our approach is that the transport problem reduces to that of a two-state problem of the particle transitions between the surface and the bulk. The latter can be solved with relative ease. The formalism is used to find the Laplace transforms of the first two moments of the particle displacement over the surface in time t at arbitrary values of the particle drift velocities and diffusivities on the surface and in the bulk. This allows us to analyze in detail the time dependence of the effective drift velocity of the particle on the surface, which can be highly nontrivial.

Facilitated diffusion in the presence of obstacles on the DNA

Biological functions of DNA depend on the sequence-specific binding of DNA-binding proteins to their corresponding binding sites. Binding of these proteins to their binding sites occurs through a facilitated diffusion process that combines three-dimensional diffusion in the cytoplasm with one-dimensional diffusion (sliding) along the DNA. In this work, we use a lattice model of facilitated diffusion to study how the dynamics of binding of a protein to a specific site (e.g., binding of an RNA polymerase to a promoter or of a transcription factor to its operator site) is affected by the presence of other proteins bound to the DNA, which act as 'obstacles' in the sliding process. Different types of these obstacles with different dynamics are implemented. While all types impair facilitated diffusion, the extent of the hindrance depends on the type of obstacle. As a consequence of hindrance by obstacles, more excursions into the cytoplasm are required for optimal target binding compared to the case without obstacles.

Dynamics of the Protein Search for Targets on DNA in the Presence of Traps

Protein search for specific binding sites on DNA is a fundamental biological phenomenon associated with the beginning of most major biological processes. It is frequently found that proteins find and recognize their specific targets quickly and efficiently despite the complex nature of protein-DNA interactions in living cells. Although significant experimental and theoretical efforts were made in recent years, the mechanisms of these processes remain not well-clarified. We present a theoretical study of the protein target search dynamics in the presence of semispecific binding sites which are viewed as traps. Our theoretical approach employs a discrete-state stochastic method that accounts for the most important physical and chemical processes in the system. It also leads to a full analytical description for all dynamic properties of the protein search. It is found that the presence of traps can significantly modify the protein search dynamics. This effect depends on the spatial positions of the targets and traps, on distances between them, on the average sliding length of the protein along the DNA, and on the total length of DNA. Theoretical predictions are discussed using simple physical-chemical arguments, and they are also validated with extensive Monte Carlo computer simulations.

A new approach to the problem of bulk-mediated surface diffusion

This paper is devoted to bulk-mediated surface diffusion of a particle which can diffuse both on a flat surface and in the bulk layer above the surface. It is assumed that the particle is on the surface initially (at t = 0) and at time t, while in between it may escape from the surface and come back any number of times. We propose a new approach to the problem, which reduces its solution to that of a two-state problem of the particle transitions between the surface and the bulk layer, focusing on the cumulative residence times spent by the particle in the two states. These times are random variables, the sum of which is equal to the total observation time t. The advantage of the proposed approach is that it allows for a simple exact analytical solution for the double Laplace transform of the conditional probability density of the cumulative residence time spent on the surface by the particle observed for time t. This solution is used to find the Laplace transform of the particle mean square displacement and to analyze the peculiarities of its time behavior over the entire range of time. We also establish a relation between the double Laplace transform of the conditional probability density and the Fourier-Laplace transform of the particle propagator over the surface. The proposed approach treats the cases of both finite and infinite bulk layer thicknesses (where bulk-mediated surface diffusion is normal and anomalous at asymptotically long times, respectively) on equal footing.

Predicting the influence of long-range molecular interactions on macroscopic-scale diffusion by homogenization of the Smoluchowski equation

The macroscopic diffusion constant for a charged diffuser is in part dependent on (1) the volume excluded by solute "obstacles" and (2) long-range interactions between those obstacles and the diffuser. Increasing excluded volume reduces transport of the diffuser, while long-range interactions can either increase or decrease diffusivity, depending on the nature of the potential. We previously demonstrated [P. M. Kekenes-Huskey et al., Biophys. J. 105, 2130 (2013)] using homogenization theory that the configuration of molecular-scale obstacles can both hinder diffusion and induce diffusional anisotropy for small ions. As the density of molecular obstacles increases, van der Waals (vdW) and electrostatic interactions between obstacle and a diffuser become significant and can strongly influence the latter's diffusivity, which was neglected in our original model. Here, we extend this methodology to include a fixed (time-independent) potential of mean force, through homogenization of the Smoluchowski equation. We consider the diffusion of ions in crowded, hydrophilic environments at physiological ionic strengths and find that electrostatic and vdW interactions can enhance or depress effective diffusion rates for attractive or repulsive forces, respectively. Additionally, we show that the observed diffusion rate may be reduced independent of non-specific electrostatic and vdW interactions by treating obstacles that exhibit specific binding interactions as "buffers" that absorb free diffusers. Finally, we demonstrate that effective diffusion rates are sensitive to distribution of surface charge on a globular protein, Troponin C, suggesting that the use of molecular structures with atomistic-scale resolution can account for electrostatic influences on substrate transport. This approach offers new insight into the influence of molecular-scale, long-range interactions on transport of charged species, particularly for diffusion-influenced signaling events occurring in crowded cellular environments.

First-passage times, mobile traps, and Hopf bifurcations

For a random walk on a confined one-dimensional domain, we consider mean first-passage times (MFPT) in the presence of a mobile trap. The question we address is whether a mobile trap can improve capture times over a stationary trap. We consider two scenarios: a randomly moving trap and an oscillating trap. In both cases, we find that a stationary trap actually performs better (in terms of reducing expected capture time) than a very slowly moving trap; however, a trap moving sufficiently fast performs better than a stationary trap. We explicitly compute the thresholds that separate the two regimes. In addition, we find a surprising relation between the oscillating trap problem and a moving-sink problem that describes reduced dynamics of a single spike in a certain regime of the Gray-Scott model. Namely, the above-mentioned threshold corresponds precisely to a Hopf bifurcation that induces oscillatory motion in the location of the spike. We use this correspondence to prove the uniqueness of the Hopf bifurcation.

An integrated model of transcription factor diffusion shows the importance of intersegmental transfer and quaternary protein structure for target site finding

We present a computational model of transcription factor motion that explains both the observed rapid target finding of transcription factors, and how this motion influences protein and genome structure. Using the Smoldyn software, we modelled transcription factor motion arising from a combination of unrestricted 3D diffusion in the nucleoplasm, sliding along the DNA filament, and transferring directly between filament sections by intersegmental transfer. This presents a fine-grain picture of the way in which transcription factors find their targets two orders of magnitude faster than 3D diffusion alone allows. Eukaryotic genomes contain sections of nucleosome free regions (NFRs) around the promoters; our model shows that the presence and size of these NFRs can be explained as their acting as antennas on which transcription factors slide to reach their targets. Additionally, our model shows that intersegmental transfer may have shaped the quaternary structure of transcription factors: sequence specific DNA binding proteins are unusually enriched in dimers and tetramers, perhaps because these allow intersegmental transfer, which accelerates target site finding. Finally, our model shows that a ‘hopping’ motion can emerge from 3D diffusion on small scales. This explains the apparently long sliding lengths that have been observed for some DNA binding proteins observed in vitro. Together, these results suggest that transcription factor diffusion dynamics help drive the evolution of protein and genome structure.

Single molecule technologies

Single molecule technologies provide an alternative set of approaches to conventional techniques and promise to deliver fundamentally new information about biological processes at the level of molecular movement, dynamics and function. As several of these mature and become more accessible for routine use in molecular biology laboratories, the potential impact on drug discovery research and development should be significant.:

Breathing dynamics based parameter sensitivity analysis of hetero-polymeric DNA

We study the parameter sensitivity of hetero-polymeric DNA within the purview of DNA breathing dynamics. The degree of correlation between the mean bubble size and the model parameters is estimated for this purpose for three different DNA sequences. The analysis leads us to a better understanding of the sequence dependent nature of the breathing dynamics of hetero-polymeric DNA. Out of the 14 model parameters for DNA stability in the statistical Poland-Scheraga approach, the hydrogen bond interaction ε(hb)(AT) for an AT base pair and the ring factor ξ turn out to be the most sensitive parameters. In addition, the stacking interaction ε(st)(TA-TA) for an TA-TA nearest neighbor pair of base-pairs is found to be the most sensitive one among all stacking interactions. Moreover, we also establish that the nature of stacking interaction has a deciding effect on the DNA breathing dynamics, not the number of times a particular stacking interaction appears in a sequence. We show that the sensitivity analysis can be used as an effective measure to guide a stochastic optimization technique to find the kinetic rate constants related to the dynamics as opposed to the case where the rate constants are measured using the conventional unbiased way of optimization.

Physical constraints determine the logic of bacterial promoter architectures

Site-specific transcription factors (TFs) bind to their target sites on the DNA, where they regulate the rate at which genes are transcribed. Bacterial TFs undergo facilitated diffusion (a combination of 3D diffusion around and 1D random walk on the DNA) when searching for their target sites. Using computer simulations of this search process, we show that the organization of the binding sites, in conjunction with TF copy number and binding site affinity, plays an important role in determining not only the steady state of promoter occupancy, but also the order at which TFs bind. These effects can be captured by facilitated diffusion-based models, but not by standard thermodynamics. We show that the spacing of binding sites encodes complex logic, which can be derived from combinations of three basic building blocks: switches, barriers and clusters, whose response alone and in higher orders of organization we characterize in detail. Effective promoter organizations are commonly found in the E. coli genome and are highly conserved between strains. This will allow studies of gene regulation at a previously unprecedented level of detail, where our framework can create testable hypothesis of promoter logic.

Intracellular facilitated diffusion: searchers, crowders, and blockers

In bacteria, regulatory proteins search for a specific DNA-binding target via "facilitated diffusion": a series of rounds of three-dimensional diffusion in the cytoplasm, and one-dimensional (1D) linear diffusion along the DNA contour. Using large scale Brownian dynamics simulations we find that each of these steps is affected differently by crowding proteins, which can either be bound to the DNA acting as a road block to the 1D diffusion, or freely diffusing in the cytoplasm. Macromolecular crowding can strongly affect mechanistic features such as the balance between three-dimensional and 1D diffusion, but leads to surprising robustness of the total search time.

Effect of DNA conformation on facilitated diffusion

Within a living cell, site-specific DNA-binding proteins need to search the whole genome to find a target of ~10-20 bp. That they find the target, and do so quickly, is vital for the correct functioning of the DNA, and of the cell as a whole. The current understanding is that this search is performed via facilitated diffusion, i.e. by combining three-dimensional bulk diffusion within the cytoplasm or nucleoplasm, with one-dimensional diffusion along the DNA backbone, to which the protein binds non-specifically. After reviewing the standard theory of facilitated diffusion, we discuss in the present article the still rather rare direct computer simulations of this process, focusing on the three-dimensional part of the search, and the effect of DNA looping and the general DNA conformation on its efficiency. We close by highlighting some open questions in this field.

Lévy noise induced switch in the gene transcriptional regulatory system

The study of random fluctuations in a gene transcriptional regulatory system is extended to the case of non-Gaussian Lévy noise, which can describe unpredictable jump changes of the random environment. The stationary probability densities are given to explore the key roles of Lévy noise in a gene transcriptional regulatory system. The results demonstrate that the parameters of Lévy noise, including noise intensity, stability index, and skewness parameter, can induce switches between distinct gene-expression states. A further concern is the switching time (from the high concentration state to the low concentration one or from the low concentration state to the high concentration one), which is a random variable and often referred to as the mean first passage time. The effects of Lévy noise on expression and degradation time are studied by computing the mean first passage time in two directions and a number of different peculiarities of non-Gaussian Lévy noise compared with Gaussian noise are observed.

Speed-selectivity paradox in the protein search for targets on DNA: is it real or not?

Protein search for targets on DNA starts all major biological processes. Although significant experimental and theoretical efforts have been devoted to investigation of these phenomena, mechanisms of protein-DNA interactions during the search remain not fully understood. One of the most surprising observations is known as a speed-selectivity paradox. It suggests that experimentally observed fast findings of targets require smooth protein-DNA binding potentials, while the stability of the specific protein-DNA complex imposes a large energy gap which should significantly slow down the protein molecule. We developed a discrete-state stochastic approach that allowed us to investigate explicitly target search phenomena and to analyze the speed-selectivity paradox. A general dynamic phase diagram for different search regimes is constructed. The effect of the target position on search dynamics is investigated. Using experimentally observed parameters, it is found that slow protein diffusion on DNA does not lead to an increase in the search times. Thus, our theory resolves the speed-selectivity paradox by arguing that it does not exist. It is just an artifact of using approximate continuum theoretical models for analyzing protein search in the region of the parameter space beyond the range of validity of these models. In addition, the presented method, for the first time, provides an explanation for fast target search at the level of single protein molecules. Our theoretical predictions agree with all available experimental observations, and extensive Monte Carlo computer simulations are performed to support analytical calculations.

In vivo facilitated diffusion model

Under dilute in vitro conditions transcription factors rapidly locate their target sequence on DNA by using the facilitated diffusion mechanism. However, whether this strategy of alternating between three-dimensional bulk diffusion and one-dimensional sliding along the DNA contour is still beneficial in the crowded interior of cells is highly disputed. Here we use a simple model for the bacterial genome inside the cell and present a semi-analytical model for the in vivo target search of transcription factors within the facilitated diffusion framework. Without having to resort to extensive simulations we determine the mean search time of a lac repressor in a living E. coli cell by including parameters deduced from experimental measurements. The results agree very well with experimental findings, and thus the facilitated diffusion picture emerges as a quantitative approach to gene regulation in living bacteria cells. Furthermore we see that the search time is not very sensitive to the parameters characterizing the DNA configuration and that the cell seems to operate very close to optimal conditions for target localization. Local searches as implied by the colocalization mechanism are only found to mildly accelerate the mean search time within our model.

Computational models for large-scale simulations of facilitated diffusion

The binding of site-specific transcription factors to their genomic target sites is a key step in gene regulation. While the genome is huge, transcription factors belong to the least abundant protein classes in the cell. It is therefore fascinating how short the time frame is that they require to home in on their target sites. The underlying search mechanism is called facilitated diffusion and assumes a combination of three-dimensional diffusion in the space around the DNA combined with one-dimensional random walk on it. In this review, we present the current understanding of the facilitated diffusion mechanism and identify questions that lack a clear or detailed answer. One way to investigate these questions is through stochastic simulation and, in this manuscript, we support the idea that such simulations are able to address them. Finally, we review which biological parameters need to be included in such computational models in order to obtain a detailed representation of the actual process.

Bullied no more: when and how DNA shoves proteins around

The predominant protein-centric perspective in protein-DNA-binding studies assumes that the protein drives the interaction. Research focuses on protein structural motifs, electrostatic surfaces and contact potentials, while DNA is often ignored as a passive polymer to be manipulated. Recent studies of DNA topology, the supercoiling, knotting, and linking of the helices, have shown that DNA has the capability to be an active participant in its transactions. DNA topology-induced structural and geometric changes can drive, or at least strongly influence, the interactions between protein and DNA. Deformations of the B-form structure arise from both the considerable elastic energy arising from supercoiling and from the electrostatic energy. Here, we discuss how these energies are harnessed for topology-driven, sequence-specific deformations that can allow DNA to direct its own metabolism.

Generalized facilitated diffusion model for DNA-binding proteins with search and recognition states

Transcription factors (TFs) such as the lac repressor find their target sequence on DNA at remarkably high rates. In the established Berg-von Hippel model for this search process, the TF alternates between three-dimensional diffusion in the bulk solution and one-dimensional sliding along the DNA chain. To overcome the so-called speed-stability paradox, in similar models the TF was considered as being present in two conformations (search state and recognition state) between which it switches stochastically. Combining both the facilitated diffusion model and alternating states, we obtain a generalized model. We explicitly treat bulk excursions for rodlike chains arranged in parallel and consider a simplified model for coiled DNA. Compared to previously considered facilitated diffusion models, corresponding to limiting cases of our generalized model, we surprisingly find a reduced target search rate. Moreover, at optimal conditions there is no longer an equipartition between the time spent by the protein on and off the DNA chain.

Effective surface motion on a reactive cylinder of particles that perform intermittent bulk diffusion

In many biological and small scale technological applications particles may transiently bind to a cylindrical surface. In between two binding events the particles diffuse in the bulk, thus producing an effective translation on the cylindrical surface. We here derive the effective motion on the surface allowing for additional diffusion on the cylindrical surface itself. We find explicit solutions for the number of adsorbed particles at one given instant, the effective surface displacement, as well as the surface propagator. In particular sub- and superdiffusive regimes are found, as well as an effective stalling of diffusion visible as a plateau in the mean squared displacement. We also investigate the corresponding first passage problem.

Theory on the dynamic memory in the transcription-factor-mediated transcription activation

We develop a theory to explain the origin of the static and dynamical memory effects in transcription-factor-mediated transcription activation. Our results suggest that the following inequality conditions should be satisfied to observe such memory effects: (a) τ(L)≫max(τ(R),τ(E)), (b) τ(LT)≫τ(T), and (c) τ(I)≥(τ(EL)+τ(TR)) where τ(L) is the average time required for the looping-mediated spatial interactions of enhancer-transcription-factor complex with the corresponding promoter--RNA-polymerase or eukaryotic RNA polymerase type II (PolII in eukaryotes) complex that is located L base pairs away from the cis-acting element, (τ(R),τ(E)) are respectively the search times required for the site-specific binding of the RNA polymerase and the transcription factor with the respective promoter and the cis-regulatory module, τ(LT) is the time associated with the relaxation of the looped-out segment of DNA that connects the cis-acting site and promoter, τ(T) is the time required to generate a complete transcript, τ(I) is the transcription initiation time, τ(EL) is the elongation time, and τ(TR) is the termination time. We have theoretically derived the expressions for the various searching, looping, and loop-relaxation time components. Using the experimentally determined values of various time components we further show that the dynamical memory effects cannot be experimentally observed whenever the segment of DNA that connects the cis-regulatory element with the promoter is not loaded with bulky histone bodies. Our analysis suggests that the presence of histone-mediated compaction of the connecting segment of DNA can result in higher values of looping and loop-relaxation times, which is the origin of the static memory in the transcription activation that is mediated by the memory gene loops in eukaryotes.

Physics of protein-DNA interactions: mechanisms of facilitated target search

One of the most critical aspects of protein-DNA interactions is the ability of protein molecules to quickly find and recognize specific target sequences on DNA. Experimental measurements indicate that the corresponding association rates to few specific sites among large number of non-specific sites are typically large. For some proteins they might be even larger than maximal allowed three-dimensional diffusion rates. Although significant progress in understanding protein search and recognition of targets on DNA has been achieved, detailed mechanisms of these processes are still strongly debated. Here we present a critical review of current theoretical approaches and some experimental observations in this area. Specifically, the role of lowering dimensionality, non-specific interactions, diffusion along the DNA molecules, protein and target sites concentrations, and electrostatic effects are critically analyzed. Possible future directions and outstanding problems are also presented and discussed.

Theory of site-specific DNA-protein interactions in the presence of conformational fluctuations of DNA binding domains

We develop a theory that explains how the thermally driven conformational fluctuations in the DNA binding domains (DBDs) of the DNA binding proteins (DBPs) are effectively coupled to the one-dimensional searching dynamics of DBPs for their cognate sites on DNA. We show that the rate gammaopt, associated with the flipping of conformational states of DBDs beyond which the maximum search efficiency of DBPs is achieved, varies with the one-dimensional sliding length L as gammaopt proportional, L(-2) and with the number of roadblock protein molecules present on the same DNA m as gammaopt proportional, m2. The required free energy barrier ERTO associated with this flipping transition seems to be varying with L as ERTO proportional, variant ln L2. When the barrier height associated with the conformational flipping of DBDs is comparable with that of the thermal free energy, then the possible value of L under in vivo conditions seems to be L<or=70 bps.

Dynamic strategies for target-site search by DNA-binding proteins

Gene regulatory proteins find their target sites on DNA remarkably quickly; the experimental binding rate for lac repressor is orders-of-magnitude higher than predicted by free diffusion alone. It has been proposed that nonspecific binding aids the search by allowing proteins to slide and hop along DNA. We develop a reaction-diffusion theory of protein translocation that accounts for transport both on and off the strand and incorporates the physical conformation of DNA. For linear DNA modeled as a wormlike chain, the distribution of hops available to a protein exhibits long, power-law tails that make the long-time displacement along the strand superdiffusive. Our analysis predicts effective superdiffusion coefficients for given nonspecific binding and unbinding rate parameters. Translocation rate exhibits a maximum at intermediate values of the binding rate constant, while search efficiency is optimized at larger binding rate constant values. Thus, our theory predicts a region of values of the nonspecific binding and unbinding rate parameters that balance the protein translocation rate and the efficiency of the search. Published data for several proteins falls within this predicted region of parameter values.

Biophysical characterization of DNA binding from single molecule force measurements

Single molecule force spectroscopy is a powerful method that uses the mechanical properties of DNA to explore DNA interactions. Here we describe how DNA stretching experiments quantitatively characterize the DNA binding of small molecules and proteins. Small molecules exhibit diverse DNA binding modes, including binding into the major and minor grooves and intercalation between base pairs of double-stranded DNA (dsDNA). Histones bind and package dsDNA, while other nuclear proteins such as high mobility group proteins bind to the backbone and bend dsDNA. Single-stranded DNA (ssDNA) binding proteins slide along dsDNA to locate and stabilize ssDNA during replication. Other proteins exhibit binding to both dsDNA and ssDNA. Nucleic acid chaperone proteins can switch rapidly between dsDNA and ssDNA binding modes, while DNA polymerases bind both forms of DNA with high affinity at distinct binding sites at the replication fork. Single molecule force measurements quantitatively characterize these DNA binding mechanisms, elucidating small molecule interactions and protein function.

A model of anomalous enzyme-catalyzed gel degradation kinetics

We show that a model of target location involving n noninteracting particles moving subdiffusively along a line segment (a generalization of a model introduced by Sokolov et al. [Biophys. J. 2005, 89, 895.]) provides a basis for understanding recent experiments by Pelta et al. [Phys. Rev. Lett. 2007, 98, 228302.] on the kinetics of diffusion-limited gel degradation. These experiments find that the time t(c) taken by the enzyme thermolysin to completely hydrolyze a gel varies inversely as roughly the 3/2 power of the initial enzyme concentration [E]. In general, however, this time would be expected to vary either as [E](-1) or as [E](-2), depending on whether the brownian diffusion of the enzyme to the site of cleavage took place along the network chains (1-d diffusion) or through the pore spaces (3-d diffusion). In our model, the unusual dependence of t(c) on [E] is explained in terms of a reaction-diffusion equation that is formulated in terms of fractional rather than ordinary time derivatives.

Microscopic mechanism for experimentally observed anomalous elasticity of DNA in two dimensions

By exploring a recent model in which DNA bending elasticity, described by the wormlike chain model, is coupled to basepair denaturation, we demonstrate that small denaturation bubbles lead to anomalies in the flexibility of DNA at the nanometric scale, when confined in two dimensions (2D), as reported in atomic-force microscopy experiments. Our model yields very good fits to experimental data and quantitative predictions that can be tested experimentally. Although such anomalies exist when DNA fluctuates freely in three dimensions (3D), they are too weak to be detected. Interactions between bases in the helical double-stranded DNA are modified by electrostatic adsorption on a 2D substrate, which facilitates local denaturation. This work reconciles the apparent discrepancy between observed 2D and 3D DNA elastic properties and points out that conclusions about the 3D properties of DNA (and its companion proteins and enzymes) do not directly follow from 2D experiments by atomic-force microscopy.

Interaction of bacteriophage T4 and T7 single-stranded DNA-binding proteins with DNA

Bacteriophages T4 and T7 are well-studied model replication systems, which have allowed researchers to determine the roles of many proteins central to DNA replication, recombination and repair. Here we summarize and discuss the results from two recently developed single-molecule methods to determine the salt-dependent DNA-binding kinetics and thermodynamics of the single-stranded DNA (ssDNA)-binding proteins (SSBs) from these systems. We use these methods to characterize both the equilibrium double-stranded DNA (dsDNA) and ssDNA binding of the SSBs T4 gene 32 protein (gp32) and T7 gene 2.5 protein (gp2.5). Despite the overall two-orders-of-magnitude weaker binding of gp2.5 to both forms of DNA, we find that both proteins exhibit four-orders-of-magnitude preferential binding to ssDNA relative to dsDNA. This strong preferential ssDNA binding as well as the weak dsDNA binding is essential for the ability of both proteins to search dsDNA in one dimension to find available ssDNA-binding sites at the replication fork.

Packaging effects on site-specific DNA-protein interactions

We show that the rate of site-specific association of a protein molecule of interest with the DNA chain can be approximately 10(2) times higher than that of the three-dimensional diffusion-controlled collision rate limit approximately 10(8) mol(-1) s(-1) only when the protein molecule of interest searches for its specific site on the DNA chain in a reduced dimensional space with a dimensionality dr of dr<1. Upon considering the concurrent dynamics of the linear DNA chain that is embedded in a d-dimensional space along with the one-dimensional diffusion dynamics of the nonspecifically bound protein molecule on the DNA chain, we derive the generalized scaling law epsilon approximately 2(3(2-d)+3), where epsilon is the number of times by which the rate of site-specific association of the protein molecule with the DNA chain can be enhanced over the three-dimensional diffusion-controlled collision rate limit and d is the dimensionality of the reduced search space. Using the analogy between the self-intersection loop length in the theory of random walks and the ring-closure events in the theory of site specific interactions of a protein molecule with the DNA chain, we further show that the extent of packaging and volume compression of the genomic DNA inside the living cell is designed in such a way that the efficiency of the protein molecule in the process of searching for its specific site on the genomic DNA is a maximum. Our simulation results suggest that the volume compression factor theta which is the ratio between the total volume of the living cell and the volume occupied by the DNA chain along with all the other bound protein molecules should be such that theta>or=100 for an efficient site specific interaction of a protein molecule of interest with the linear DNA chain that is embedded in a three-dimensional space. Our theoretical and simulation results agree well with the E. coli cellular system.

Using the bias from flow to elucidate single DNA repair protein sliding and interactions with DNA

We perform single-molecule spatial tracking measurements of a DNA repair protein, the C-terminal domain of Ada (C-Ada) from Escherichia coli, moving on DNA extended by flow. The trajectories of single proteins labeled with a fluorophore are constructed. We analyze single-protein dwell times on DNA for different flow rates and conclude that sliding (with essentially no hopping) is the mechanism of C-Ada motion along stretched DNA. We also analyze the trajectory results with a drift-diffusion Langevin equation approach to elucidate the influence of flow on the protein motion; systematic variation of the flow enables one to estimate the microscopic friction. We integrate the step-size probability distribution to obtain a version of the fluctuation theorem that articulates the relation between the entropy production and consumption under the adjustable drag (i.e., bias) from the flow. This expression allows validation of the Langevin equation description of the motion. Comparison of the rate of sliding with recent computer simulations of DNA repair suggests that C-Ada could conduct its repair function while moving at near the one-dimensional diffusion limit.

Facilitated diffusion with DNA coiling

When DNA-binding proteins search for their specific binding site on a DNA molecule they alternate between linear 1-dimensional diffusion along the DNA molecule, mediated by nonspecific binding, and 3-dimensional volume excursion events between successive dissociation from and rebinding to DNA. If the DNA molecule is kept in a straight configuration, for instance, by optical tweezers, these 3-dimensional excursions may be divided into long volume excursions and short hops along the DNA. These short hops correspond to immediate rebindings after dissociation such that a rebinding event to the DNA occurs at a site that is close to the site of the preceding dissociation. When the DNA molecule is allowed to coil up, immediate rebinding may also lead to so-called intersegmental jumps, i.e., immediate rebindings to a DNA segment that is far away from the unbinding site when measured in the chemical distance along the DNA, but close by in the embedding 3-dimensional space. This effect is made possible by DNA looping. The significance of intersegmental jumps was recently demonstrated in a single DNA optical tweezers setup. Here we present a theoretical approach in which we explicitly take the effect of DNA coiling into account. By including the spatial correlations of the short hops we demonstrate how the facilitated diffusion model can be extended to account for intersegmental jumping at varying DNA densities. It is also shown that our approach provides a quantitative interpretation of the experimentally measured enhancement of the target location by DNA-binding proteins.

Diffusion, dimensionality, and noise in transcriptional regulation

The precision of biochemical signaling is limited by randomness in the diffusive arrival of molecules at their targets. For proteins binding to specific sites on DNA and regulating transcription, the ability of the proteins to diffuse in one dimension by sliding along the length of the DNA, in addition to their diffusion in bulk solution, would seem to generate a larger target for DNA binding, consequently reducing the noise in the occupancy of the regulatory site. Here we show that this effect is largely canceled by the enhanced temporal correlations in one-dimensional diffusion. With realistic parameters, sliding along DNA has surprisingly little effect on the physical limits to the precision of transcriptional regulation.

Bubble merging in breathing DNA as a vicious walker problem in opposite potentials

We investigate the coalescence of two DNA bubbles initially located at weak domains and separated by a more stable barrier region in a designed construct of double-stranded DNA. In a continuum Fokker-Planck approach, the characteristic time for bubble coalescence and the corresponding distribution are derived, as well as the distribution of coalescence positions along the barrier. Below the melting temperature, we find a Kramers-type barrier crossing behavior, while at high temperatures, the bubble corners perform drift diffusion toward coalescence. In the calculations, we map the bubble dynamics on the problem of two vicious walkers in opposite potentials. We also present a discrete master equation approach to the bubble coalescence problem. Numerical evaluation and stochastic simulation of the master equation show excellent agreement with the results from the continuum approach. Given that the coalesced state is thermodynamically stabilized against a state where only one or a few of the base pairs of the barrier region are re-established, it appears likely that this type of setup could be useful for the quantitative investigation of thermodynamic DNA stability data as well as the rate constants involved in the unzipping and zipping dynamics of DNA in single molecule fluorescence experiments.

Directional dependent dynamics of protein molecules on DNA

We demonstrate that a protein molecule of interest undergoing the one-dimensional Brownian dynamics along DNA can exhibit a directional dependent net transport either toward or away from its target site depending on the distribution of the initial positions of the other classes of protein molecules present on the same DNA. Directionality arises as a consequence of the confinement of the search space and dynamic reflections by other protein molecules present on the same DNA chain. Energy cost for such directionality comes from the free energy spent on setting the initial positions of the other protein molecules. In the mechanism of action of cis-acting elements on the initiation of transcription, such free-energy inputs are derived from the site-specific binding affinities of the inflowing transcriptional factors toward their cis-acting elements. If the initial distribution of other protein molecules is a random one, then the protein molecule of interest exhibits a net transport away from its target site. This directionality originates from unequal natures of enhancing and retarding effects of the randomly distributed other classes of protein molecules. The protein molecule of interest overcomes the retarding effects of the other classes of protein molecules in a dynamical manner by increasing the number of dissociation-association events when it is far away from its target site and then by switching back to the sliding dynamics due to increase in the enhancing effects as it moves closer to its target site.

Single DNA denaturation and bubble dynamics

While the Watson-Crick double-strand is the thermodynamically stable state of DNA in a wide range of temperature and salt conditions, even at physiological conditions local denaturation bubbles may open up spontaneously due to thermal activation. By raising the ambient temperature, titration, or by external forces in single molecule setups bubbles proliferate until full denaturation of the DNA occurs. Based on the Poland-Scheraga model we investigate both the equilibrium transition of DNA denaturation and the dynamics of the denaturation bubbles with respect to recent single DNA chain experiments for situations below, at, and above the denaturation transition. We also propose a new single molecule setup based on DNA constructs with two bubble zones to measure the bubble coalescence and extract the physical parameters relevant to DNA breathing. Finally we consider the interplay between denaturation bubbles and selectively single-stranded DNA binding proteins.

Optimizing intermittent reaction paths

Various examples of biochemical reactions in cells, such as DNA/protein interactions, reveal that in extremely diluted regimes reaction paths are not always simple brownian trajectories. They can rather be qualified as intermittent, since they combine slow diffusion phases on one hand and a second mode of faster transport on the other hand, which can be either a faster diffusion mode, as in the case of DNA-binding proteins, or a ballistic mode powered by molecular motors in the case of intracellular transport. In this article, we introduce simple theoretical models which permit to calculate explicitly the reaction rates for reactions limited by intermittent transport. This approach shows quantitatively that intermittent reaction pathways are actually very efficient, since they permit to significantly increase the reaction rates, which could explain why they are observed so often. Moreover, we give theoretical arguments which suggest that intermittent transport could also be useful for in vitro chemistry. Indeed, we show that intermittent transport naturally pops up in the context of reaction at interfaces, where reactants combine surface diffusion phases and bulk excursions, and could permit to enhance reactivity. In this case, adjusting chemically the affinity of reactants with the interface makes possible to optimize the reaction rate.

How DNA coiling enhances target localization by proteins

Many genetic processes depend on proteins interacting with specific sequences on DNA. Despite the large excess of nonspecific DNA in the cell, proteins can locate their targets rapidly. After initial nonspecific binding, they are believed to find the target site by 1D diffusion ("sliding") interspersed by 3D dissociation/reassociation, a process usually referred to as facilitated diffusion. The 3D events combine short intrasegmental "hops" along the DNA contour, intersegmental "jumps" between nearby DNA segments, and longer volume "excursions." The impact of DNA conformation on the search pathway is, however, still unknown. Here, we show direct evidence that DNA coiling influences the specific association rate of EcoRV restriction enzymes. Using optical tweezers together with a fast buffer exchange system, we obtained association times of EcoRV on single DNA molecules as a function of DNA extension, separating intersegmental jumping from other search pathways. Depending on salt concentration, targeting rates almost double when the DNA conformation is changed from fully extended to a coiled configuration. Quantitative analysis by an extended facilitated diffusion model reveals that only a fraction of enzymes are ready to bind to DNA. Generalizing our results to the crowded environment of the cell we predict a major impact of intersegmental jumps on target localization speed on DNA.

Kinetics and thermodynamics of salt-dependent T7 gene 2.5 protein binding to single- and double-stranded DNA

Bacteriophage T7 gene 2.5 protein (gp2.5) is a single-stranded DNA (ssDNA)-binding protein that has essential roles in DNA replication, recombination and repair. However, it differs from other ssDNA-binding proteins by its weaker binding to ssDNA and lack of cooperative ssDNA binding. By studying the rate-dependent DNA melting force in the presence of gp2.5 and its deletion mutant lacking 26 C-terminal residues, we probe the kinetics and thermodynamics of gp2.5 binding to ssDNA and double-stranded DNA (dsDNA). These force measurements allow us to determine the binding rate of both proteins to ssDNA, as well as their equilibrium association constants to dsDNA. The salt dependence of dsDNA binding parallels that of ssDNA binding. We attribute the four orders of magnitude salt-independent differences between ssDNA and dsDNA binding to nonelectrostatic interactions involved only in ssDNA binding, in contrast to T4 gene 32 protein, which achieves preferential ssDNA binding primarily through cooperative interactions. The results support a model in which dimerization interactions must be broken for DNA binding, and gp2.5 monomers search dsDNA by 1D diffusion to bind ssDNA. We also quantitatively compare the salt-dependent ssDNA- and dsDNA-binding properties of the T4 and T7 ssDNA-binding proteins for the first time.