Volume 12 Supplement 10

## Proceedings of the Eighth Annual MCBIOS Conference. Computational Biology and Bioinformatics for a New Decade

# First-passage time analysis of a one-dimensional diffusion-reaction model: application to protein transport along DNA

- Michael L Mayo
^{1}Email author, - Edward J Perkins
^{1}and - Preetam Ghosh
^{2}

**12(Suppl 10)**:S18

**DOI: **10.1186/1471-2105-12-S10-S18

© Mayo et al; licensee BioMed Central Ltd. 2011

**Published: **18 October 2011

## Abstract

### Background

Proteins search along the DNA for targets, such as transcription initiation sequences, according to one-dimensional diffusion, which is interrupted by micro- and macro-hopping events and intersegmental transfers that occur under close packing conditions.

### Results

A one-dimensional diffusion-reaction model in the form of difference-differential equations is proposed to analyze the nonequilibrium protein sliding kinetics along a segment of bacterial DNA. A renormalization approach is used to derive an expression for the mean first-passage time to arrive at sites downstream of the origin from the occupation probabilities given by the individual transport equations. Monte Carlo simulations are employed to assess the validity of the proposed approach, and all results are interpreted within the context of bacterial transcription.

### Conclusions

Mean first-passage times decrease with increasing reaction rates, indicating that, on average, surviving proteins more rapidly locate downstream targets than their reaction-free counterparts, but at the price of increasing rarity. Two qualitatively different screening regimes are identified according to whether the search process operates under “small” or “large” values for the dissociation rate of the protein-DNA complex. Lower bounds are placed on the overall search time for varying reactive conditions. Good agreement with experimental estimates requires the reaction rate reside near the transition between both screening regimes, suggesting that biology balances a need for rapid searches against maximum exploration during each round of the sliding phase.

## Background

Information in the form of nucleotide sequences is processed into RNA during transcription, which are then processed into proteins responsible for regulating signaling pathways and other elements of cellular chemistry. Proteins, called transcription factors, are transported to the DNA by diffusion, wherein they attach to DNA “receptors,” forming semi-stable complexes. Initiation of transcription occurs after a search along DNA for specific terminating sequences; their presence attenuates the affinity of RNA polymerase to the promoter region of the targeted gene, increasing or decreasing transcriptional activity. Early experimental work on the *lac* repressor in *E. coli*[1, 2], and the seminal theoretical works of Winter, Berg, and von Hippel [3–5], laid the groundwork in the current understanding of protein-search kinetics. Targets (i.e., transcription initiation sequences) are located after conducting a one-dimensional search along the DNA, transported there by diffusion through the cytoplasm, wherein sliding is interrupted by micro- and macro-hopping events along single segments and intersegmental protein transfers across DNA strands in close proximity under close packing conditions. Such rounds of *jumping* and *sliding* across the DNA are now collectively referred to as “facilitated target location” [6]. In vivo, pure sliding is complicated by protein-phosphate backbone interactions [7], energetic and steric restrictions (e.g. DNA conformational changes [8, 9]), access limitations induced by diffusion gradients [10], macromolecular crowding [11], as well as proteins binding to multiple segments [12], to name just a few. Moreover, single molecule experiments conclude that protein association to the DNA dominates dissociation, wherein the resulting protein-DNA complex is relatively stable during the sliding process [13–15].

*t*

_{ s }, needed to find the transcription initiation sequence, which is the sum of the times spent sliding and jumping during each round:

*t*

_{ i }(Λ) is the time needed for the protein to traverse a length Λ (in units of base-pairs, bps) for each sliding event,

*t*

_{3D,i}is the time a protein spends free of the DNA beteween each round of sliding, and

*N*

_{ s }is the number of sliding rounds in the search. Equation 1 may be simplified by taking its average, , across all rounds of sliding and jumping [16], resulting in:

*t*

_{ i }(Λ)〉 =

*t*

_{Λ}labels the average time needed to explore a length Λ from the site of association, termed the exploration length [17]; 〈

*t*

_{3D,i}〉 =

*t*

_{3D}labels the average time a protein spends free of the DNA during the search. Furthermore, it has been argued that biology is optimized to minimize the overall search time [16], so that

*t*

_{Λ}=

*t*

_{3D}. Under this assumption, Eqn. 2 may be expressed as:

So, under suitable assumptions, the calculation of the total search time is reduced to an estimation of the length of the sequence actually searched, *L*, the exploration length, Λ, and the average time spent sliding during each round, *t*_{Λ}.

While others have investigated the transport properties of the protein-DNA complex in the presence of reversible binding kinetics [3], our analysis of the overall search times, based on Eq. 3, reveals new insight into the extraordinary speed of the target-location process. We frame our discussion in terms of a first-passage time analysis, from which we estimate *t*_{λ} for a one-dimensional diffusion-reaction model proposed below, and use it to show that the reaction rate serves to screen proteins from downstream exploration of the DNA-a fact that has implications for the timely response of mRNA transcription to the presence of transcription factors.

## Methods

Proteins are much smaller in length than the DNA they attach to. Indeed, ignoring the conformation of either protein or DNA, an *E. coli* genome is ~ 10^{5} times longer than a protein bound to it. A typical genome is ~ 10^{6} bases (~ 1*/* 3 mm, given that the distance between bases is Δ*l* = 1/3 nm), and the median “length” of an E. coli protein is approximately 278 amino acids [18], or ≈ 278*Å*; the ratio of protein to genome length is therefore approximately 8 × 10^{–5}. So, the DNA appears appears to be a long segment of nucleotide bases, as viewed from the perspective of DNA-bound proteins. It is natural, then, to model a long sequence as a one-dimensional lattice, wherein the lattice sites correspond to the base-pair positions of the underlying DNA.

*i*= 1, termed the origin, and terminates upon successful dissociation of the protein-DNA complex at any intermediate site

*i*= 2,3,…,

*N –*1, wherein

*N*is the number of bases in the modeled segment. The exact biological terminating mechanism of the target-location event is ignored: we merely assume that a particle trap is hosted at site

*i*=

*N*, which is implemented as a perfectly absorbing boundary condition (see below). The exact way in which the Hydrogen donor and acceptor bond-patterns destabilize the minimum interaction potential at the terminating site established during the protein-slide [19], is of no direct importance to the problem at hand.

### Energetic basis for protein diffusion along DNA

*i*(with respect to the cytoplasm),

*U*

_{ i }, is decomposed into two contributions: the specific, or sequence-dependent, component, , and the nonspecific, or sequence-independent, component,

*E*

_{ ns }. This decomposition gives:

wherein
is a binding sequence (*s*_{
i
} ∈ {*A*,*C*,*G*,*T*}) of length *l*.

wherein *μ* is the standard deviation of the distribution.

Although varying correlation lengths in the protein-DNA binding free energy may arise in several ways-such as directly, from the nonspecific energy [23], DNA curvature [8], and nonrandom sequences with high protein affinity (e.g. AT/GC-rich isochors [8]), or effectively, by mapping correlated micro-hops to a sliding process by “reducing” the free energy barriers to transport [9, 24]-on larger scales, Eqn. 5 provides evidence to support an assumption of an uncorrelated binding free energy landscape. This idea has formed the basis of other studies [16].

*μ*, which measures the value of an average fluctuation:

and estimated previously to be *μ* ~ *k*_{
B
}*T* under biologically relevant conditions (written in terms of the system temperature, *T*, and Boltzmann’s constant *k*_{
B
}) [8, 16].

*U*

_{i±1}–

*U*

_{ i }(Fig. 1). It can be shown that the transition probability to reach adjacent sites is related to the uninhibited transition rate,

*ν*, by an Arrhenius-type equation [16, 25]:

Equation 7 describes a DNA segment in thermal contact with cytoplasm of temperature *T*, 1/*β* = *k*_{
B
}*T* is the absolute temperature, and the factor of one half results from normalization. These transition probabilities are, in general, not symmetric; i.e., *ω*_{i+1,i} ≠ *ω*_{i,i+1}. Because the binding energies, *U*_{
i
}, are considered here to be Gaussian-distributed random variables, equation 7 implies that *ω*_{i+1,i} + *ω*_{i–1,i} ≤ *ν;* the difference being the probability per unit time for the protein not to jump to either site, termed the sojourn probability [26]. Furthermore, at any time during the walk, the protein-DNA complex may dissociate with rate *r*_{
i
} (Fig. 1), triggering extraneous protein macro- and micro-hopping events.

The biological problem of protein diffusion along the DNA under dissociative conditions is therefore mapped onto the generic problem of one-dimensional diffusion that proceeds according to randomly chosen transition probabilities, and subject to reactive conditions with the model lattice. In the next section we provide nonequilibrium transport equations that, collectively, serve as the mathematical foundation for the sliding process under study.

### Transport equations for the one-dimensional model

Given a protein that is currently associated to the lattice site *i*, it may move to adjacent sites according to the following “rules.” It may i) hop to site *i* + 1 with probability *α*_{i+1,i} = *ω*_{i+1,i}/*ν*; ii) hop to lattice site *i –* 1 with probability *α*_{
i
}_{–1,}_{
i
}; iii) react with, or dissociate from, the lattice site *i* with probability *σ*_{
i
} (1 – *α*_{i+1,i} – *α*_{i–1,i}), wherein *σ*_{
i
} = *r*_{
i
}/*ν* is the probability to react with the lattice at site *i* in time 1/*ν* . Dissociation from the DNA may occur, for example, by interacting with on-site obstacles during the slide [27].

*i*= 1 and terminating at site

*i*=

*n*, the occupation probability for any site

*i*along the segment,

*p*

_{i,1}(

*t*), at the elapsed time

*t*, can be found by formulating the diffusion rules i-iii) above as a set of difference-differential equations for

*n*lattice sites:

Note that reaction rates, transition probabilities, and times, are scaled in units of 1/*ν*, rendering the transport equations dimensionless.

*ν*~ 10

^{8}

*Hz–a*value which justifies the use of continuous-time differential equations for the protein diffusion problem. Additionally, the probability is sharp at the origin, so that initial data for equations 8-10 are given by:

wherein *δ*_{i 1}=0 (when *i* ≠ 1) or 1 (when *i* = 1) is the Kronecker delta.

Note that boundary conditions 8 and 10 must be specified to provide a unique set of occupation probabilities for the entire lattice; we have chosen to provide a reflecting boundary condition at the source of the diffusion. Although the diffusion properties are initially influenced by this boundary condition, boundary effects vanish for the late times and distances that are relevant at biological time (~ 1 second) and length scales (~ 100 bases), because of the uncorrelated nature of the diffusive trajectory.

### Renormalization of the first-passage time distribution

The mean first-passage time, *τ*_{
i
}, measures the average time for a diffusing object to arrive at a target site *i* for the first time, and is often used as an estimate for reaction rates of creation/descruction across a boundary, and transit times through a medium, among other phenomena of interest [28]. For a protein diffusing along the DNA, the mean first-passage time (MFPT) is the average *minimum* search time needed to locate a target site downstream from the site of association with the DNA, because transport is always terminated at the target sequence. As transcription factors at the promoter sequences of a gene either recruit or inhibit RNA polymerase binding, the rate at which they first arrive at the promoter sequence might also provide a rough estimate of a gene’s transcriptional activity.

Here, the term *p*_{i,i} (*t* – *t*′) denotes the probability of the protein transported to site *i*, to leave and return again by sliding to site *i* in time *t* – *t*′ ≥ 0; the protein may leave and return many times in this interval. Finally, the term *F*_{i,1} (*t*′) gives the probability per unit time that a protein will reach the site *i* for the first time, starting from the origin, at time *t*′ ≤ *t*.

*i*=

*n*, is equal to the rate in which probability enters into it:

because any walker that enters is restricted from leaving to contribute to occupation probabilities of any other sites, according to Eq. 10.

*n*–only the properties of those proteins surviving to this site are of interest. Equation 13 should therefore be renormalized to account for only those walkers surviving the entire transit to site

*n*:

*n*= 2, 3,…,

*N*, we have:

Note that the first-passage time for the origin is trivial, *τ*_{1} = 0.

### Monte Carlo simulations

Exact numerical simulations of the diffusion process, as defined above by the diffusion rules i-iii), were conducted by means of Monte Carlo simulation to study the diffusion-reaction process in terms of the mean first-passage time for a protein to reach a downstream target site. Individual trajectories were simulated, and the first-passage times were recorded for each site along the lattice. An arithmetic average of these first-passage times gives an approximation of the mean (equation 15).

*ν*). Figure 2 illustrates the pseudo-code for MCRW, logically outlined as follows:

- 1.
A 2D matrix,

*total_simulation_step*, records the number of jumps experienced by the protein for each unique trajectory. - 2.
A 2D matrix,

*fpt*, records the first-passage times for all the lattice sites visited for each trajectory. - 3.
A 1D array,

*visit*, records the particular site that is visited for each jump in any particular trajectory. This array is overwritten upon termination of each trajectory, and can be used to compute other statistical metrics, such as the root-mean-square displacement, among others.

This simulation terminates when a random walker survives to reach the target site, for a total of 1000 times, wherein a counter for complete trajectories is incremented and another unique trajectory begins. Similarly, if a reaction event returns TRUE, the current trajectory terminates and another is started again at the origin, wherein reflecting boundary conditions (equation 8) are given.

Note that while we have only discussed the measurement of the mean first-passage times using the algorithm MCRW, it is straightforward to compute the root-mean-square displacements for the different trajectories, as well as other statistical endpoints as they vary with the dissociation probability, random energy background, target site, and so on.

## Results

### Preferential selection of fast-moving proteins under reactive conditions

*μ*= 0.75

*k*

_{ B }

*T*), for a sequence of 100 bases in length, chosen because experiments on the BbvCI restriction enzyme suggest that Λ = 50 bases [30], but it is also long enough to minimize the effects of the reflecting boundary at the origin. It has been argued that the reaction rate depends almost entirely on nonspecific energy contributions to the total binding free energy [16],

*E*

_{ ns }, which is not dependent on the actual nucleotide sequence. We adopt this view by setting the reaction rates for each site equal to one another and dropping their indeces,

*r*=

*σν*. The influence of site-dependent reaction rates are considered below.

The mean first-passage time generally decreases with increasing reaction rate *σ*, demonstrating that only the “faster” trajectories-i.e., those clearing a maximum of bases in an equivalent amount of time-survive along the segment during each sliding phase. However, as *σ* increases, all trajectories become increasingly rare, in the sense that downstream sites are not visited as often as the ones associated with smaller *σ*. This indicates that transcription cannot entirely rely upon sliding as the main transport mechanism-an idea generally accepted since the early seminal works (see Background). If sliding were the only mechanism, a balance between minimizing search times while maximizing the length of individual slides might be required to most efficiently use the limited transcription factor proteins available in the cytoplasm (there may be as few as ~ 10-100 transcription factors present in the cytoplasm at any given time). However, since dissociation initiates the bulk diffusion involved in micro- and macro-hopping of the protein along the DNA, increasing distance covered and decreasing total search times, a larger reaction rate (but not too large) may actually support the search process by increasing the number of bulk excursions, allowing for the molecule to travel dominantly by 3D diffusion to distal DNA binding sites.

In the reactive cases, an increase in *σ* filters all but the superdiffusive trajectories (*x* ~ *t*^{
d
}, *d* > 1/2) from those of the reaction-free case, leaving those that scale directly with time: *x* ~ *t*. This speed-up is consistent with intuition, in that ever faster trajectories should increasingly acquire the properties of the maximal one. This clearly cannot hold in the limit *σ* → 1, as no molecules survive on the lattice to be transported anywhere downstream of the origin.

*σ*

_{ i }, were chosen according to the distribution shown within the figure inset. Here, the full width at half maximum of the distribution is chosen to be equal to the mean 〈

*σ*

_{ i }〉 = 10

^{–3}. Even under such a broad distribution (or equally, conditions of large fluctuations), the mean first-passage times are remarkably stable across the length of the shown segment, supporting the use of a site-independent reaction rate used throughout this article.

### Larger energetic fluctuations slow a purely sliding target search

*μ*, equation 6. Figure 5 illustrates the disorder-averaged mean first-passage times (y-axis) for several locations along the segment as a function of the parameters

*μ*and

*σ*. The first-passage time characteristics are roughly the same for all lattice sites: smaller energetic fluctuations allow for the fastest transits under most reaction conditions. As the distribution, equation 5, becomes sharper, the energetic barriers to transport vanish, and the probability for an adjacent transition tends to 1, i.e.

*α*

_{i+1,i}+

*α*

_{i–1,i}→ 1.

Interestingly, roughness values ranging from *μ* = 4*k*_{
B
}*T* to 16*k*_{
B
}*T* provide very similiar transport rates. As the distribution, equation 5, becomes wider, larger energetic barriers quench the transition frequency to adjacent sites, increasing residence times at the current site and slowing the sliding component of the overall search. Other studies suggest that a purely energetic terminating process would require *μ* ≥ 4*k*_{
B
}*T*[16], which, according to figure 5, generates slower searches under all reactive conditions. The “optimal” roughness reported in these studies, *μ* = 0.75*k*_{
B
}*T*, resides within a region of low variation in the mean first-passage times across a wide range of energetic conditions, wherein the roughness spans values ranging from *μ* = 1/2 to 1. Distrubtions adopting such values generate similar search times; such homogeneity under changing conditions provides a large degree of flexibility in the search process to many different types of proteins, and may provide tolerance against disruptions, such as temperature fluctuations.

### Survival times of proteins associated to DNA

*t*), quantifies this qualitative feature by giving the probability to find a protein anywhere along the segment as a function of time. Often of interest is the half-life of the protein-DNA comlex,

*T*

_{1/2}, which is the time required for the survival probability to be exactly 1

*/*2, and is found by inverting the following expression:

*μ*= 0.75

*k*

_{ B }

*T*), and the results are illustrated within Fig. 6 This figure suggests that for any downstream exploration to occur by sliding alone, the dissociation rate of the protein-DNA complex must be very small, because only the fastest molecules survive to significantly explore downstream sites for increasing reaction rates. Moreover, the power-law type scaling behavior,

*T*

_{1/2}~

*σ*

^{ γ }, transitions from

*γ*= 0 to

*γ*= –1 rather quickly at approximately

*σ*= 10

^{–4}, indicating that the influence of proteins trapped by the far boundary of the segment (

*i*=

*N*) washes out the contributions from any reactions with the medium for smaller

*σ*. However, dissociation events dominate the half-life behavior for larger values of

*σ*.

### Validation of the model and its biological implications

*L*and

*t*

_{λ}, given Λ = 50 bases. If we assume the span of the sequence that is actually searched is on the order of the length of the genome,

*L*= 10

^{6}bases, and we make the assumption that the search process terminates the first time a target sequence is found,

*t*

_{λ}=

*τ*

_{50}, then Eq. 3 reduces to:

^{ th }base (left axis) as a function of the dissociation/reaction rate for one particular realization of the binding free energy landscape (shown within the inset of Fig. 3). On the other axis (right axis), estimates for the minimum overall search time are provided, calculated according to equation 17, with the maximum hopping frequency given by

*ν*= 10

^{8}Hz.

*k*

_{ on }≈ 10

^{8}– 10

^{9}

*M*

^{–1}

*s*

^{–1}. Assuming a cytoplasm transcription factor concentration of ≈ 10

^{–9}

*M*, an order-of-magnitude estimate can be made for the time needed for a protein to bind one site [16]:

This value is in good agreement with the values predicted by the diffusion-reaction model presented here, as shown by the right axis of Fig. 7. Specifically, to agree with the lower bound of such estimates (i.e. ~ 1 second), the reaction rate should be bounded from above by *σ* ≤ 10^{–3}, placing it near the plateau region of figure 7.

There are two power-law type regimes here, *d* = 0 and *d* = –0.53, for scaling of type *τ*_{50} ~ *σ*^{
d
}. Here, the exponent of *d* ≈ -1*/* 2 means that a 100-fold increase in the reaction rate reduces the search time by 10-fold in the regime characterized by *σ* > 10^{–4}. Moreover, the transition the two regimes is abrupt, occuring across a decade or so in the reaction rate. The intersection of these power laws occurs at *σ* ≈ 10^{–4}, indicating that a significant number of molecules survive along the segment in biological conditions (refer to the half-life estimates of Fig. 6), so that dissociations of the protein-DNA complex during the search along nonspecific DNA are infrequent.

These two regimes can be interpreted as screening regimes, in that proteins are prevented from reaching downstream sites due to dissociation of the protein-DNA complex. The abrupt transition between these regimes can be attributed to a partial screening regime, wherein exploration is limited by the reaction rate, and representing a balance between maximizing the speed of the sliding phase while allowing for the greatest chance to arrive downstream. Biological advantage would appear to rest somewhere near this transition.

The plateau in the search times for *σ* < 10^{–4}, suggests that it may be experimentally difficult to measure reaction rates for these semi-stable comlexes, because the first-passage times vary very little across several orders of magnitude studied. In other words, the relationship between search times and reaction rate is not invertible in the region defined by *σ* < 10^{–4}. Nevertheless, the search time associated with the plateau, approximately *t*_{
min
} ≈ 2.4 seconds, sits comfortably within the bound of the experimental estimate provided above.

### Validity of the first-passage renormalization approach

## Conclusions

A one-dimensional diffusion-reaction model, in the form of a set of difference-differential equations has been proposed to analyze the nonequilibrium protein sliding kinetics along a segment of bacterial DNA. An iterative renormalization procedure was used to express the mean first-passage time to arrive at a downstream site in terms of the occupation probabilities provided by the transport equations. The validity of this approach was established by implementing a Monte Carlo type algorithm that directly compared the renormalization approach to first-passage times measured by the exact simulations.

We found that downstream targets are screened from access by upstream proteins according to the value of the dissociation rate of the protein-DNA complex; however, anomalous characteristics of the surviving associated proteins allow for faster search times as compared to their reaction-free counterparts. Moreover, access limitations induced by dissociation of the protein-DNA complex suggests that nature balances downstream exploration with the competing need for faster search times, possibly allowing for the most efficient use of resources, i.e., the few transcription factors typically present for each gene.

Indeed, experimental estimates of overall target-location times are in good agreement with a minimum bound presented here, indicating that the entire search process is flexible enough to handle a certain amount of variability associated with the increased rarity of DNA-bound proteins (depending on the actual reaction rate). It is clear from these analyses that sliding alone cannot univocally minimize the overall search times, and that other transport mechanisms, such as micro- and macro-hopping events and intersegment transfers, must assist in the search process to provide reliable delivery of proteins to target sites–a result that is consistent with current experimental and theoretical data. Further work is needed, however, to reveal the exact proportion in which these separate mechanisms contribute to the overall search process.

## Declarations

### Acknowledgements

Funding was provided by the US Army Environmental Quality and Installations 6.1 Basic Research Program. The Chief of Engineers has approved this material for publication.

This article has been published as part of *BMC Bioinformatics* Volume 12 Supplement 10, 2011: Proceedings of the Eighth Annual MCBIOS Conference. Computational Biology and Bioinformatics for a New Decade. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/12?issue=S10.

## Authors’ Affiliations

## References

- Riggs AD, Bourgeois S, Cohn M:
**The lac repressor-operator interaction. III. Kinetic studies.***J. Mol. Biol*1970,**53:**401–418. 10.1016/0022-2836(70)90074-4View ArticlePubMedGoogle Scholar - Riggs AD, Suzuki S, Bourgeois S:
**Lac repressor-operator interaction. I. Equilibrium studies.***J. Mol. Biol*1970,**48:**67–83. 10.1016/0022-2836(70)90219-6View ArticlePubMedGoogle Scholar - Berg OG, Winter RB, von Hippel PH:
**Diffusion-driven mechanisms of protein translocation on nucleic acids. 1. Models and theory.***Biochem*1981,**20:**6929–48. 10.1021/bi00527a028View ArticleGoogle Scholar - Winter RB, von Hippel PH:
**Diffusion-driven mechanisms of protein translocation on nucleic acids.***Biochem*1981,**20:**6948–6960. 10.1021/bi00527a029View ArticleGoogle Scholar - Winter RB, Berg OG, von Hippel PH:
**The Escherichia coli lac repressor-operator interaction: equilibrium measurements.***Biochem*1981,**20:**6961. 10.1021/bi00527a030View ArticleGoogle Scholar - von Hippel PH, Berg OG:
**Facilitated target location in biological systems.***J. Biol. Chem*1989,**264:**675–8.PubMedGoogle Scholar - Luscombe NM, Laskowski RA, Thornton JM:
**Amino acid-base interactions: a three-dimensional analysis of protein-DNA interactions at an atomic level.***Nucl. Acids Res*2001,**29:**2860–2874. 10.1093/nar/29.13.2860PubMed CentralView ArticlePubMedGoogle Scholar - Slutsky M, Kardar M, Mirny LA:
**Diffusion in correlated random potentials, with applications to DNA.***Phys. Rev. E*2004,**69:**061903.View ArticleGoogle Scholar - Wunderlich Z, Mirny LA:
**Spatial effects on the speed and reliability of protein-DNA search.***Nucl. Acids Res*2008,**36:**3570. 10.1093/nar/gkn173PubMed CentralView ArticlePubMedGoogle Scholar - Pfeifer P, Hagerty P:
**Screening transition in diffusion to and across fractal surfaces.**In*Fractals and Chaos in Chemical Engineering*. Edited by: Giona M, Biardi G. Singapore: World Scientific; 1997:151–164.Google Scholar - Li GW, Berg OG, Elf J:
**Effects of macromolecular crowding and DNA looping on gene regulation kinetics.***Nature Physics*2009,**5:**294–297. 10.1038/nphys1222View ArticleGoogle Scholar - Fried MG, Crothers DM:
**Kinetics and mechanism in the reaction of gene regulatory proteins with DNA.***J. Mol. Biol*1984,**172:**263–282. 10.1016/S0022-2836(84)80026-1View ArticlePubMedGoogle Scholar - Shimamoto N:
**One dimensional diffusion of proteins along DNA: its biological and chemical significance revealed by single-molecule measurements.***J. Biol. Chem*1999,**274:**15293–15296. 10.1074/jbc.274.22.15293View ArticlePubMedGoogle Scholar - Blainey PC, van Oijen AM, Banerjee A, Verdine GL, Xie XS:
**A base-excision DNA-repair protein finds intrahelical lesion bases by fast sliding in contact with DNA.***Proc. Natl. Acad. Sci. USA*2006,**103:**5752–5757. 10.1073/pnas.0509723103PubMed CentralView ArticlePubMedGoogle Scholar - Wang YM, Austin RH, Cox EC:
**Single-Molecule Measurements of Repressor Protein 1-D Diffusion on DNA.***Phys. Rev. Lett*2006,**97:**048302.View ArticlePubMedGoogle Scholar - Slutsky M, Mirny L:
**Kinetics of protein-DNA interaction: facilitated target location in sequence-dependent potential.***Biophys. J*2004,**87:**4021–35. 10.1529/biophysj.104.050765PubMed CentralView ArticlePubMedGoogle Scholar - Pfeifer P, Sapoval B:
**Optimization of diffusive transport to irregular surfaces with low sticking probability.***Mat. Res. Soc. Symp. Proc*1995,**366:**271–276.View ArticleGoogle Scholar - Brocchieri L, Karlin S:
**Protein length in eukaryotic and prokaryotic proteomes.***Nucl. Acids Res*2005,**33:**3390–3400. 10.1093/nar/gki615PubMed CentralView ArticlePubMedGoogle Scholar - Dahirel V, Paillusson F, Jardat M, Barbi M, Victor JM:
**Nonspecific DNA-protein interaction: why proteins can diffuse along DNA.***Phys. Rev. Lett*2009,**102:**228101.View ArticlePubMedGoogle Scholar - Berg OG, von Hippel PH:
**Selection of DNA binding sites by regulatory proteins. Statistical-mechanical theory and application to operators and promoters.***J. Mol. Biol*1987,**193:**723–750. 10.1016/0022-2836(87)90354-8View ArticlePubMedGoogle Scholar - Gerland U, Moroz JD, Hwa T:
**Physical constraints and functional characteristics of transcription factor-DNA interaction.***Proc. Natl. Acad. Sci. USA*2002,**99:**12015–12020. 10.1073/pnas.192693599PubMed CentralView ArticlePubMedGoogle Scholar - Stormo GD, Fields DS:
**Specificity, free energy and information content in protein-DNA interactions.***Trends Biochem. Sci*1998,**23:**109–113. 10.1016/S0968-0004(98)01187-6View ArticlePubMedGoogle Scholar - Das RK, Kolomeisky AB:
**Facilitated search of proteins on DNA: correlations are important.***Phys. Chem. Chem. Phys*2010,**12:**2999–3004.View ArticlePubMedGoogle Scholar - Loverdo C, Benichou O, Voituriez R, Biebricher A, Bonnet I, Desbiolles P:
**Quantifying hopping and jumping in facilitated diffusion of DNA-binding proteins.***Phys. Rev. Lett*2009,**102:**188101.View ArticlePubMedGoogle Scholar - Risken H:
*The Fokker-Planck Equation: Methods of Solution and Applications.*Berlin: Springer; 1989.View ArticleGoogle Scholar - Murthy KPN, Kehr KW:
**Mean first-passage time of random walks on a random lattice.***Phys. Rev. A*1989,**40:**2082–2087. 10.1103/PhysRevA.40.2082View ArticlePubMedGoogle Scholar - Zhou HX:
**A model for the mediation of processivity of DNA-targeting proteins by nonspecific binding: dependence on DNA length and presence of obstacles.***Biophys. J*2005,**88:**1608–1615. 10.1529/biophysj.104.052688PubMed CentralView ArticlePubMedGoogle Scholar - Redner S:
*A Guide to First-Passage Processes.*Cambridge: Cambridge University Press; 2001.View ArticleGoogle Scholar - Hughes B:
*Random walks and random environments.*Oxford: Oxford Science Publications; 1995.Google Scholar - Gower DM, Wilson GG, Halford SE:
**Measurement of the contributions of 1D and 3D pathways to the translocation of a protein along DNA.***Proc. Natl. Acad. Sci. USA*2005,**102:**15883–15888. 10.1073/pnas.0505378102View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.