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

How nonspecifically DNA-binding proteins search for the target in crowded environments

We investigate how a tracer particle searches a target located in DNA modeled by a stiff chain in crowded environments using theoretical analysis and Langevin dynamics simulations. First, we show that the three-dimensional (3D) diffusion coefficient of the tracer only depends on the density of crowders ϕ, while its one-dimensional (1D) diffusion coefficient is affected by not only ϕ but also the nonspecific binding energy ε. With increasing ϕ and ε, no obvious change in the average 3D diffusion time is observed, while the average 1D sliding time apparently increases. We propose theoretically that the 1D sliding of the tracer along the chain could be well captured by the Kramers' law of escaping rather than the Arrhenius law, which is verified directly by the simulations. Finally, the average search time increases monotonously with an increase in ϕ while it has a minimum as a function of ε, which could be understood from the different behaviors of the average number of search rounds with the increasing ϕ or ε. These results provide a deeper understanding of the role of facilitated diffusion in target search of proteins on DNA in vivo.

Rad4 recognition-at-a-distance: Physical basis of conformation-specific anomalous diffusion of DNA repair proteins

Since Robert Brown's first observations of random walks by pollen particles suspended in solution, the concept of diffusion has been subject to countless theoretical and experimental studies in diverse fields from finance and social sciences, to physics and biology. Diffusive transport of macromolecules in cells is intimately linked to essential cellular functions including nutrient uptake, signal transduction, gene expression, as well as DNA replication and repair. Advancement in experimental techniques has allowed precise measurements of these diffusion processes. Mathematical and physical descriptions and computer simulations have been applied to model complicated biological systems in which anomalous diffusion, in addition to simple Brownian motion, was observed. The purpose of this review is to provide an overview of the major physical models of anomalous diffusion and corresponding experimental evidence on the target search problem faced by DNA-binding proteins, with an emphasis on DNA repair proteins and the role of anomalous diffusion in DNA target recognition.

Integrating Optical Tweezers, DNA Tightropes, and Single-Molecule Fluorescence Imaging: Pitfalls and Traps

Fluorescence imaging is one of the cornerstone techniques for understanding how single molecules search for their targets on DNA. By tagging individual proteins, it is possible to track their position with high accuracy. However, to understand how proteins search for targets, it is necessary to elongate the DNA to avoid protein localization ambiguities. Such structures known as "DNA tightropes" are tremendously powerful for imaging target location; however, they lack information about how force and load affect protein behavior. The use of optically trapped microstructures offers the means to apply and measure force effects. Here we describe a system that we recently developed to enable individual proteins to be directly manipulated on DNA tightropes. Proteins bound to DNA can be conjugated with Qdot fluorophores for visualization and also directly manipulated by an optically trapped, manufactured microstructure. Together this offers a new approach to understanding the physical environment of molecules, and the combination with DNA tightropes presents opportunities to study complex biological phenomena.

Tetrameric structure of the restriction DNA glycosylase R.PabI in complex with nonspecific double-stranded DNA

R.PabI is a type II restriction enzyme that recognizes the 5′-GTAC-3′ sequence and belongs to the HALFPIPE superfamily. Although most restriction enzymes cleave phosphodiester bonds at specific sites by hydrolysis, R.PabI flips the guanine and adenine bases of the recognition sequence out of the DNA helix and hydrolyzes the N-glycosidic bond of the flipped adenine in a similar manner to DNA glycosylases. In this study, we determined the structure of R.PabI in complex with double-stranded DNA without the R.PabI recognition sequence by X-ray crystallography. The 1.9 Å resolution structure of the complex showed that R.PabI forms a tetrameric structure to sandwich the double-stranded DNA and the tetrameric structure is stabilized by four salt bridges. DNA binding and DNA glycosylase assays of the R.PabI mutants showed that the residues that form the salt bridges (R70 and D71) are essential for R.PabI to find the recognition sequence from the sea of nonspecific sequences. R.PabI is predicted to utilize the tetrameric structure to bind nonspecific double-stranded DNA weakly and slide along it to find the recognition sequence.

Human Lin28 Forms a High-Affinity 1:1 Complex with the 106~363 Cluster miRNA miR-363

Lin28A is a post-transcriptional regulator of gene expression that interacts with and negatively regulates the biogenesis of let-7 family miRNAs. Recent data suggested that Lin28A also binds the putative tumor suppressor miR-363, a member of the 106~363 cluster of miRNAs. Affinity for this miRNA and the stoichiometry of the protein-RNA complex are unknown. Characterization of human Lin28's interaction with RNA has been complicated by difficulties in producing stable RNA-free protein. We have engineered a maltose binding protein fusion with Lin28, which binds let-7 miRNA with a Kd of 54.1 ± 4.2 nM, in agreement with previous data on a murine homologue. We show that human Lin28A binds miR-363 with a 1:1 stoichiometry and with a similar, if not higher, affinity (Kd = 16.6 ± 1.9 nM). Further analysis suggests that the interaction of the N-terminal cold shock domain of Lin28A with RNA is salt-dependent, supporting a model in which the cold shock domain allows the protein to sample RNA substrates through transient electrostatic interactions.

What have single-molecule studies taught us about gene expression?

The production of a single mRNA is the result of many sequential steps, from docking of transcription factors to polymerase initiation, elongation, splicing, and, finally, termination. Much of our knowledge about the fundamentals of RNA synthesis and processing come from ensemble in vitro biochemical measurements. Single-molecule approaches are very much in this same reductionist tradition but offer exquisite sensitivity in space and time along with the ability to observe heterogeneous behavior and actually manipulate macromolecules. These techniques can also be applied in vivo, allowing one to address questions in living cells that were previously restricted to reconstituted systems. In this review, we examine the unique insights that single-molecule techniques have yielded on the mechanisms of gene expression.

DNA binding proteins explore multiple local configurations during docking via rapid rebinding

Finding the target site and associating in a specific orientation are essential tasks for DNA-binding proteins. In order to make the target search process as efficient as possible, proteins should not only rapidly diffuse to the target site but also dynamically explore multiple local configurations before diffusing away. Protein flipping is an example of this second process that has been observed previously, but the underlying mechanism of flipping remains unclear. Here, we probed the mechanism of protein flipping at the single molecule level, using HIV-1 reverse transcriptase (RT) as a model system. In order to test the effects of long-range attractive forces on flipping efficiency, we varied the salt concentration and macromolecular crowding conditions. As expected, increased salt concentrations weaken the binding of RT to DNA while increased crowding strengthens the binding. Moreover, when we analyzed the flipping kinetics, i.e. the rate and probability of flipping, at each condition we found that flipping was more efficient when RT bound more strongly. Our data are consistent with a view that DNA bound proteins undergo multiple rapid re-binding events, or short hops, that allow the protein to explore other configurations without completely dissociating from the DNA.

Theory on the mechanism of site-specific DNA-protein interactions in the presence of traps

The speed of site-specific binding of transcription factor (TFs) proteins with genomic DNA seems to be strongly retarded by the randomly occurring sequence traps. Traps are those DNA sequences sharing significant similarity with the original specific binding sites (SBSs). It is an intriguing question how the naturally occurring TFs and their SBSs are designed to manage the retarding effects of such randomly occurring traps. We develop a simple random walk model on the site-specific binding of TFs with genomic DNA in the presence of sequence traps. Our dynamical model predicts that (a) the retarding effects of traps will be minimum when the traps are arranged around the SBS such that there is a negative correlation between the binding strength of TFs with traps and the distance of traps from the SBS and (b) the retarding effects of sequence traps can be appeased by the condensed conformational state of DNA. Our computational analysis results on the distribution of sequence traps around the putative binding sites of various TFs in mouse and human genome clearly agree well the theoretical predictions. We propose that the distribution of traps can be used as an additional metric to efficiently identify the SBSs of TFs on genomic DNA.

High-speed atomic force microscopy reveals strongly polarized movement of clostridial collagenase along collagen fibrils

Bacterial collagenases involved in donor infection are widely applied in many fields due to their high activity and specificity; however, little is known regarding the mechanisms by which bacterial collagenases degrade insoluble collagen in host tissues. Using high-speed atomic force microscopy, we simultaneously visualized the hierarchical structure of collagen fibrils and the movement of a representative bacterial collagenase, Clostridium histolyticum type I collagenase (ColG), to determine the relationship between collagen structure and collagenase movement. Notably, ColG moved ~14.5 nm toward the collagen N terminus in ~3.8 s in a manner dependent on a catalytic zinc ion. While ColG was engaged, collagen molecules were not only degraded but also occasionally rearranged to thicken neighboring collagen fibrils. Importantly, we found a similarity of relationship between the enzyme-substrate interface structure and enzyme migration in collagen-collagenase and DNA-nuclease systems, which share a helical substrate structure, suggesting a common strategy in enzyme evolution.

Protein search for multiple targets on DNA

Protein-DNA interactions are crucial for all biological processes. One of the most important fundamental aspects of these interactions is the process of protein searching and recognizing specific binding sites on DNA. A large number of experimental and theoretical investigations have been devoted to uncovering the molecular description of these phenomena, but many aspects of the mechanisms of protein search for the targets on DNA remain not well understood. One of the most intriguing problems is the role of multiple targets in protein search dynamics. Using a recently developed theoretical framework we analyze this question in detail. Our method is based on a discrete-state stochastic approach that takes into account most relevant physical-chemical processes and leads to fully analytical description of all dynamic properties. Specifically, systems with two and three targets have been explicitly investigated. It is found that multiple targets in most cases accelerate the search in comparison with a single target situation. However, the acceleration is not always proportional to the number of targets. Surprisingly, there are even situations when it takes longer to find one of the multiple targets in comparison with the single target. It depends on the spatial position of the targets, distances between them, average scanning lengths of protein molecules on DNA, and the total DNA lengths. Physical-chemical explanations of observed results are presented. Our predictions are compared with experimental observations as well as with results from a continuum theory for the protein search. Extensive Monte Carlo computer simulations fully support our theoretical calculations.

How motif environment influences transcription factor search dynamics: Finding a needle in a haystack

Transcription factors (TFs) have to find their binding sites, which are distributed throughout the genome. Facilitated diffusion is currently the most widely accepted model for this search process. Based on this model the TF alternates between one-dimensional sliding along the DNA, and three-dimensional bulk diffusion. In this view, the non-specific associations between the proteins and the DNA play a major role in the search dynamics. However, little is known about how the DNA properties around the motif contribute to the search. Accumulating evidence showing that TF binding sites are embedded within a unique environment, specific to each TF, leads to the hypothesis that the search process is facilitated by favorable DNA features that help to improve the search efficiency. Here, we review the field and present the hypothesis that TF-DNA recognition is dictated not only by the motif, but is also influenced by the environment in which the motif resides.

Twist-open mechanism of DNA damage recognition by the Rad4/XPC nucleotide excision repair complex

DNA damage repair starts with the recognition of damaged sites from predominantly normal DNA. In eukaryotes, diverse DNA lesions from environmental sources are recognized by the xeroderma pigmentosum C (XPC) nucleotide excision repair complex. Studies of Rad4 (radiation-sensitive 4; yeast XPC ortholog) showed that Rad4 "opens" up damaged DNA by inserting a β-hairpin into the duplex and flipping out two damage-containing nucleotide pairs. However, this DNA lesion "opening" is slow (˜5-10 ms) compared with typical submillisecond residence times per base pair site reported for various DNA-binding proteins during 1D diffusion on DNA. To address the mystery as to how Rad4 pauses to recognize lesions during diffusional search, we examine conformational dynamics along the lesion recognition trajectory using temperature-jump spectroscopy. Besides identifying the ˜10-ms step as the rate-limiting bottleneck towards opening specific DNA site, we uncover an earlier ˜100- to 500-μs step that we assign to nonspecific deformation (unwinding/"twisting") of DNA by Rad4. The β-hairpin is not required to unwind or to overcome the bottleneck but is essential for full nucleotide-flipping. We propose that Rad4 recognizes lesions in a step-wise "twist-open" mechanism, in which preliminary twisting represents Rad4 interconverting between search and interrogation modes. Through such conformational switches compatible with rapid diffusion on DNA, Rad4 may stall preferentially at a lesion site, offering time to open DNA. This study represents the first direct observation, to our knowledge, of dynamical DNA distortions during search/interrogation beyond base pair breathing. Submillisecond interrogation with preferential stalling at cognate sites may be common to various DNA-binding proteins.

Directly interrogating single quantum dot labelled UvrA2 molecules on DNA tightropes using an optically trapped nanoprobe

In this study we describe a new methodology to physically probe individual complexes formed between proteins and DNA. By combining nanoscale, high speed physical force measurement with sensitive fluorescence imaging we investigate the complex formed between the prokaryotic DNA repair protein UvrA₂ and DNA. This approach uses a triangular, optically-trapped “nanoprobe” with a nanometer scale tip protruding from one vertex. By scanning this tip along a single DNA strand suspended between surface-bound micron-scale beads, quantum-dot tagged UvrA₂ molecules bound to these ‘”DNA tightropes” can be mechanically interrogated. Encounters with UvrA₂ led to deflections of the whole nanoprobe structure, which were converted to resistive force. A force histogram from all 144 detected interactions generated a bimodal distribution centered on 2.6 and 8.1 pN, possibly reflecting the asymmetry of UvrA₂’s binding to DNA. These observations successfully demonstrate the use of a highly controllable purpose-designed and built synthetic nanoprobe combined with fluorescence imaging to study protein-DNA interactions at the single molecule level.

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.

Nonconsensus Protein Binding to Repetitive DNA Sequence Elements Significantly Affects Eukaryotic Genomes

Recent genome-wide experiments in different eukaryotic genomes provide an unprecedented view of transcription factor (TF) binding locations and of nucleosome occupancy. These experiments revealed that a large fraction of TF binding events occur in regions where only a small number of specific TF binding sites (TFBSs) have been detected. Furthermore, in vitro protein-DNA binding measurements performed for hundreds of TFs indicate that TFs are bound with wide range of affinities to different DNA sequences that lack known consensus motifs. These observations have thus challenged the classical picture of specific protein-DNA binding and strongly suggest the existence of additional recognition mechanisms that affect protein-DNA binding preferences. We have previously demonstrated that repetitive DNA sequence elements characterized by certain symmetries statistically affect protein-DNA binding preferences. We call this binding mechanism nonconsensus protein-DNA binding in order to emphasize the point that specific consensus TFBSs do not contribute to this effect. In this paper, using the simple statistical mechanics model developed previously, we calculate the nonconsensus protein-DNA binding free energy for the entire C. elegans and D. melanogaster genomes. Using the available chromatin immunoprecipitation followed by sequencing (ChIP-seq) results on TF-DNA binding preferences for ~100 TFs, we show that DNA sequences characterized by low predicted free energy of nonconsensus binding have statistically higher experimental TF occupancy and lower nucleosome occupancy than sequences characterized by high free energy of nonconsensus binding. This is in agreement with our previous analysis performed for the yeast genome. We suggest therefore that nonconsensus protein-DNA binding assists the formation of nucleosome-free regions, as TFs outcompete nucleosomes at genomic locations with enhanced nonconsensus binding. In addition, here we perform a new, large-scale analysis using in vitro TF-DNA preferences obtained from the universal protein binding microarrays (PBM) for ~90 eukaryotic TFs belonging to 22 different DNA-binding domain types. As a result of this new analysis, we conclude that nonconsensus protein-DNA binding is a widespread phenomenon that significantly affects protein-DNA binding preferences and need not require the presence of consensus (specific) TFBSs in order to achieve genome-wide TF-DNA binding specificity.

Interactions between proteins and DNA trigger many important biological processes. Therefore, to fully understand how the information encoded on the DNA transcribes into RNA, which in turn translates into proteins in the cell, we need to unravel the molecular design principles of protein-DNA interactions. It is known that many interactions occur when a protein is attracted to a specific short segment on the DNA called a specific protein-DNA binding motif. Strikingly, recent experiments revealed that many regulatory proteins reproducibly bind to different regions on the DNA lacking such specific motifs. This suggests that fundamental molecular mechanisms responsible for protein-DNA recognition specificity are not fully understood. Here, using high-throughput protein-DNA binding data obtained by two entirely different methods for ~100 TFs in each case, we show that DNA regions possessing certain repetitive sequence elements exert the statistical attractive potential on DNA-binding proteins, and as a result, such DNA regions are enriched in bound proteins. This is in agreement with our previous analysis performed for the yeast genome. We use the term nonconsensus protein-DNA binding in order to describe protein-DNA interactions that occur in the absence of specific protein-DNA binding motifs. Here we demonstrate that the identified nonconsensus effect is highly significant for a variety of organismal genomes and it affects protein-DNA binding preferences and nucleosome occupancy at the genome-wide level.

Molecular crowding effect on dynamics of DNA-binding proteins search for their targets

DNA-binding proteins locate and bind their target sequences positioned on DNA in crowded environments, but the molecular crowding effect on this search process is not clear. Using analytical techniques and Langevin dynamics simulations in two dimensions (2D), we find that the essential physics for facilitated diffusion in 2D search and 3D search is the same. We observe that the average search times have minima at the same optimal nonspecific binding energy for the cases with and without the crowding particle. Moreover, the molecular crowding increases the search time by increasing the average search rounds and the one-dimensional (1D) sliding time of a round, but almost not changing the average 2D diffusion time of a round. In addition, the fraction of 1D sliding time out of the total search time increases with increasing the concentration of crowders. For 2D diffusion, the molecular crowding decreases the jumping length and narrows its distribution due to the cage effect from crowders. These results shed light on the role of facilitated diffusion in DNA targeting kinetics in living cells.

Method for the analysis of contribution of sliding and hopping to a facilitated diffusion of DNA-binding protein: Application to in vivo data

DNA-binding protein searches for its target, a specific site on DNA, by means of diffusion. The search process consists of many recurrent steps of one-dimensional diffusion (sliding) along the DNA chain and three-dimensional diffusion (hopping) after dissociation of a protein from the DNA chain. Here we propose a computational method that allows extracting the contribution of sliding and hopping to the search process in vivo from the measurements of the kinetics of the target search by the lac repressor in Escherichia coli [P. Hammar et al., Science 336, 1595 (2012)]. The method combines lattice Monte Carlo simulations with the Brownian excursion theory and includes explicitly steric constraints for hopping due to the helical structure of DNA. The simulation results including all experimental data reveal that the in vivo target search is dominated by sliding. The short-range hopping to the same base pair interrupts one-dimensional sliding while long-range hopping does not contribute significantly to the kinetics of the search of the target in vivo.

Probing the target search of DNA-binding proteins in mammalian cells using TetR as model searcher

Many cellular functions rely on DNA-binding proteins finding and associating to specific sites in the genome. Yet the mechanisms underlying the target search remain poorly understood, especially in the case of the highly organized mammalian cell nucleus. Using as a model Tet repressors (TetRs) searching for a multi-array locus, we quantitatively analyse the search process in human cells with single-molecule tracking and single-cell protein–DNA association measurements. We find that TetRs explore the nucleus and reach their target by 3D diffusion interspersed with transient interactions with non-cognate sites, consistent with the facilitated diffusion model. Remarkably, nonspecific binding times are broadly distributed, underlining a lack of clear delimitation between specific and nonspecific interactions. However, the search kinetics is not determined by diffusive transport but by the low association rate to nonspecific sites. Altogether, our results provide a comprehensive view of the recruitment dynamics of proteins at specific loci in mammalian cells.

During transcription, replication and repair, DNA-binding proteins must find specific interaction sites hidden within a vast excess of genomic DNA. Here the authors use single-molecule tracking to quantitatively determine the contributions of the different processes that underlie target search in human cells.

Direct observation of TALE protein dynamics reveals a two-state search mechanism

Transcription activator-like effector (TALE) proteins are a class of programmable DNA-binding proteins for which the fundamental mechanisms governing the search process are not fully understood. Here we use single-molecule techniques to directly observe TALE search dynamics along DNA templates. We find that TALE proteins are capable of rapid diffusion along DNA using a combination of sliding and hopping behaviour, which suggests that the TALE search process is governed in part by facilitated diffusion. We also observe that TALE proteins exhibit two distinct modes of action during the search process—a search state and a recognition state—facilitated by different subdomains in monomeric TALE proteins. Using TALE truncation mutants, we further demonstrate that the N-terminal region of TALEs is required for the initial non-specific binding and subsequent rapid search along DNA, whereas the central repeat domain is required for transitioning into the site-specific recognition state.

TALEs are programmable DNA-binding proteins with practical use in genome engineering and synthetic biology. Here the authors use single-molecule fluorescence microscopy to establish that TALE proteins function using two distinct DNA-interaction modes during sequence-specific target search.

Mechanisms of diffusional search for specific targets by DNA-dependent proteins

To perform their functions, many DNA-dependent proteins have to quickly locate specific targets against the vast excess of nonspecific DNA. Although this problem was first formulated over 40 years ago, the mechanism of such search remains one of the unsolved fundamental problems in the field of protein-DNA interactions. Several complementary mechanisms have been suggested: sliding, based on one-dimensional random diffusion along the DNA contour; hopping, in which the protein "jumps" between the closely located DNA fragments; macroscopic association-dissociation of the protein-DNA complex; and intersegmental transfer. This review covers the modern state of the problem of target DNA search, theoretical descriptions, and methods of research at the macroscopic (molecule ensembles) and microscopic (individual molecules) levels. Almost all studied DNA-dependent proteins search for specific targets by combined three-dimensional diffusion and one-dimensional diffusion along the DNA contour.

Extracting physics of life at the molecular level: A review of single-molecule data analyses

Studying individual biomolecules at the single-molecule level has proved very insightful recently. Single-molecule experiments allow us to probe both the equilibrium and nonequilibrium properties as well as make quantitative connections with ensemble experiments and equilibrium thermodynamics. However, it is important to be careful about the analysis of single-molecule data because of the noise present and the lack of theoretical framework for processes far away from equilibrium. Biomolecular motion, whether it is free in solution, on a substrate, or under force, involves thermal fluctuations in varying degrees, which makes the motion noisy. In addition, the noise from the experimental setup makes it even more complex. The details of biologically relevant interactions, conformational dynamics, and activities are hidden in the noisy single-molecule data. As such, extracting biological insights from noisy data is still an active area of research. In this review, we will focus on analyzing both fluorescence-based and force-based single-molecule experiments and gaining biological insights at the single-molecule level. Inherently nonequilibrium nature of biological processes will be highlighted. Simulated trajectories of biomolecular diffusion will be used to compare and validate various analysis techniques.

Differential interaction kinetics of a bipolar structure-specific endonuclease with DNA flaps revealed by single-molecule imaging

As DNA repair enzymes are essential for preserving genome integrity, understanding their substrate interaction dynamics and the regulation of their catalytic mechanisms is crucial. Using single-molecule imaging, we investigated the association and dissociation kinetics of the bipolar endonuclease NucS from Pyrococcus abyssi (Pab) on 5′ and 3′-flap structures under various experimental conditions. We show that association of the PabNucS with ssDNA flaps is largely controlled by diffusion in the NucS-DNA energy landscape and does not require a free 5′ or 3′ extremity. On the other hand, NucS dissociation is independent of the flap length and thus independent of sliding on the single-stranded portion of the flapped DNA substrates. Our kinetic measurements have revealed previously unnoticed asymmetry in dissociation kinetics from these substrates that is markedly modulated by the replication clamp PCNA. We propose that the replication clamp PCNA enhances the cleavage specificity of NucS proteins by accelerating NucS loading at the ssDNA/dsDNA junctions and by minimizing the nuclease interaction time with its DNA substrate. Our data are also consistent with marked reorganization of ssDNA and nuclease domains occurring during NucS catalysis, and indicate that NucS binds its substrate directly at the ssDNA-dsDNA junction and then threads the ssDNA extremity into the catalytic site. The powerful techniques used here for probing the dynamics of DNA-enzyme binding at the single-molecule have provided new insight regarding substrate specificity of NucS nucleases.

Mechanistic stochastic model of histone modification pattern formation

BACKGROUND

The activity of a single gene is influenced by the composition of the chromatin in which it is embedded. Nucleosome turnover, conformational dynamics, and covalent histone modifications each induce changes in the structure of chromatin and its affinity for regulatory proteins. The dynamics of histone modifications and the persistence of modification patterns for long periods are still largely unknown.

RESULTS

In this study, we present a stochastic mathematical model that describes the molecular mechanisms of histone modification pattern formation along a single gene, with non-phenomenological, physical parameters. We find that diffusion and recruitment properties of histone modifying enzymes together with chromatin connectivity allow for a rich repertoire of stochastic histone modification dynamics and pattern formation. We demonstrate that histone modification patterns at a single gene can be established or removed within a few minutes through diffusion and weak recruitment mechanisms of histone modification spreading. Moreover, we show that strong synergism between diffusion and weak recruitment mechanisms leads to nearly irreversible transitions in histone modification patterns providing stable patterns. In the absence of chromatin connectivity spontaneous and dynamic histone modification boundaries can be formed that are highly unstable, and spontaneous fluctuations cause them to diffuse randomly. Chromatin connectivity destabilizes this synergistic system and introduces bistability, illustrating state switching between opposing modification states of the model gene. The observed bistable long-range and localized pattern formation are critical effectors of gene expression regulation.

CONCLUSION

This study illustrates how the cooperative interactions between regulatory proteins and the chromatin state generate complex stochastic dynamics of gene expression regulation.

Selective microbial genomic DNA isolation using restriction endonucleases

To improve the metagenomic analysis of complex microbiomes, we have repurposed restriction endonucleases as methyl specific DNA binding proteins. As an example, we use DpnI immobilized on magnetic beads. The ten minute extraction technique allows specific binding of genomes containing the DpnI Gm6ATC motif common in the genomic DNA of many bacteria including γ-proteobacteria. Using synthetic genome mixtures, we demonstrate 80% recovery of Escherichia coli genomic DNA even when only femtogram quantities are spiked into 10 µg of human DNA background. Binding is very specific with less than 0.5% of human DNA bound. Next Generation Sequencing of input and enriched synthetic mixtures results in over 100-fold enrichment of target genomes relative to human and plant DNA. We also show comparable enrichment when sequencing complex microbiomes such as those from creek water and human saliva. The technique can be broadened to other restriction enzymes allowing for the selective enrichment of trace and unculturable organisms from complex microbiomes and the stratification of organisms according to restriction enzyme enrichment.

Solution structure and intramolecular exchange of methyl-cytosine binding domain protein 4 (MBD4) on DNA suggests a mechanism to scan for mCpG/TpG mismatches

Unlike other members of the methyl-cytosine binding domain (MBD) family, MBD4 serves as a potent DNA glycosylase in DNA mismatch repair specifically targeting ^mCpG/TpG mismatches arising from spontaneous deamination of methyl-cytosine. The protein contains an N-terminal MBD (MBD4_MBD) and a C-terminal glycosylase domain (MBD4_GD) separated by a long linker. This arrangement suggests that the MBD4_MBD either directly augments enzymatic catalysis by the MBD4_GD or targets the protein to regions enriched for ^mCpG/TpG mismatches. Here we present structural and dynamic studies of MBD4_MBD bound to dsDNA. We show that MBD4_MBD binds with a modest preference for^mCpG as compared to mismatch, unmethylated and hydroxymethylated DNA. We find that while MBD4_MBD exhibits slow exchange between molecules of DNA (intermolecular exchange), the domain exhibits fast exchange between two sites in the same molecule of dsDNA (intramolecular exchange). Introducing a single-strand defect between binding sites does not greatly reduce the intramolecular exchange rate, consistent with a local hopping mechanism for moving along the DNA. These results support a model in which the MBD4_MBD4 targets the intact protein to ^mCpG islands and promotes scanning by rapidly exchanging between successive ^mCpG sites which facilitates repair of nearby ^mCpG/TpG mismatches by the glycosylase domain.

Microscopic mechanism of DNA damage searching by hOGG1

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

Type II restriction endonucleases--a historical perspective and more

This article continues the series of Surveys and Summaries on restriction endonucleases (REases) begun this year in Nucleic Acids Research. Here we discuss 'Type II' REases, the kind used for DNA analysis and cloning. We focus on their biochemistry: what they are, what they do, and how they do it. Type II REases are produced by prokaryotes to combat bacteriophages. With extreme accuracy, each recognizes a particular sequence in double-stranded DNA and cleaves at a fixed position within or nearby. The discoveries of these enzymes in the 1970s, and of the uses to which they could be put, have since impacted every corner of the life sciences. They became the enabling tools of molecular biology, genetics and biotechnology, and made analysis at the most fundamental levels routine. Hundreds of different REases have been discovered and are available commercially. Their genes have been cloned, sequenced and overexpressed. Most have been characterized to some extent, but few have been studied in depth. Here, we describe the original discoveries in this field, and the properties of the first Type II REases investigated. We discuss the mechanisms of sequence recognition and catalysis, and the varied oligomeric modes in which Type II REases act. We describe the surprising heterogeneity revealed by comparisons of their sequences and structures.

Positive and negative impacts of nonspecific sites during target location by a sequence-specific DNA-binding protein: origin of the optimal search at physiological ionic strength

The inducible transcription factor Egr-1, which recognizes a 9-bp target DNA sequence via three zinc-finger domains, rapidly activates particular genes upon cellular stimuli such as neuronal signals and vascular stresses. Here, using the stopped-flow fluorescence method, we measured the target search kinetics of the Egr-1 zinc-finger protein at various ionic strengths between 40 and 400 mM KCl and found the most efficient search at 150 mM KCl. We further investigated the kinetics of intersegment transfer, dissociation, and sliding of this protein on DNA at distinct concentrations of KCl. Our data suggest that Egr-1's kinetic properties are well suited for efficient scanning of chromosomal DNA in vivo. Based on a newly developed theory, we analyzed the origin of the optimal search efficiency at physiological ionic strength. Target association is accelerated by nonspecific binding to nearby sites and subsequent sliding to the target as well as by intersegment transfer. Although these effects are stronger at lower ionic strengths, such conditions also favor trapping of the protein at distant nonspecific sites, decelerating the target association. Our data demonstrate that Egr-1 achieves the optimal search at physiological ionic strength through a compromise between the positive and negative impacts of nonspecific interactions with DNA.

Ultrafast redistribution of E. coli SSB along long single-stranded DNA via intersegment transfer

Single-stranded DNA binding proteins (SSBs) selectively bind single-stranded DNA (ssDNA) and facilitate recruitment of additional proteins and enzymes to their sites of action on DNA. SSB can also locally diffuse on ssDNA, which allows it to quickly reposition itself while remaining bound to ssDNA. In this work, we used a hybrid instrument that combines single-molecule fluorescence and force spectroscopy to directly visualize the movement of Escherichia coli SSB on long polymeric ssDNA. Long ssDNA was synthesized without secondary structure that can hinder quantitative analysis of SSB movement. The apparent diffusion coefficient of E. coli SSB thus determined ranged from 70,000 to 170,000nt(2)/s, which is at least 600 times higher than that determined from SSB diffusion on short ssDNA oligomers, and is within the range of values reported for protein diffusion on double-stranded DNA. Our work suggests that SSB can also migrate via a long-range intersegment transfer on long ssDNA. The force dependence of SSB movement on ssDNA further supports this interpretation.

Obstacles may facilitate and direct DNA search by proteins

DNA recognition by DNA-binding proteins (DBPs), which is a pivotal event in most gene regulatory processes, is often preceded by an extensive search for the correct site. A facilitated diffusion process in which a DBP combines three-dimensional diffusion in solution with one-dimensional sliding along DNA has been suggested to explain how proteins can locate their target sites on DNA much faster than predicted by three-dimensional diffusion alone. Although experimental and theoretical studies have recently advanced understanding of the biophysical principles underlying the search mechanism, the process under in vivo cellular conditions is poorly understood. In this study, we used various computational approaches to explore how the presence of obstacle proteins on the DNA influences search efficiency. At a low obstacle occupancy (i.e., when few obstacles occupy sites on the DNA), sliding by the searching DBP may be confined, which may impair search efficiency. The obstacles, however, can be bypassed during hopping events, and the number of bypasses is larger for higher obstacle occupancies. Dynamism on the part of the obstacles may even further facilitate search kinetics. Our study shows that the nature and efficiency of the search process may be governed not only by the intrinsic properties of the DBP and the salt concentration of the medium, but also by the in vivo association of DNA with other macromolecular obstacles, their location, and occupancy.

The effects of transcription factor competition on gene regulation

Transcription factor (TF) molecules translocate by facilitated diffusion (a combination of 3D diffusion around and 1D random walk on the DNA). Despite the attention this mechanism received in the last 40 years, only a few studies investigated the influence of the cellular environment on the facilitated diffusion mechanism and, in particular, the influence of "other" DNA binding proteins competing with the TF molecules for DNA space. Molecular crowding on the DNA is likely to influence the association rate of TFs to their target site and the steady state occupancy of those sites, but it is still not clear how it influences the search in a genome-wide context, when the model includes biologically relevant parameters (such as: TF abundance, TF affinity for DNA and TF dynamics on the DNA). We performed stochastic simulations of TFs performing the facilitated diffusion mechanism, and considered various abundances of cognate and non-cognate TFs. We show that, for both obstacles that move on the DNA and obstacles that are fixed on the DNA, changes in search time are not statistically significant in case of biologically relevant crowding levels on the DNA. In the case of non-cognate proteins that slide on the DNA, molecular crowding on the DNA always leads to statistically significant lower levels of occupancy, which may confer a general mechanism to control gene activity levels globally. When the "other" molecules are immobile on the DNA, we found a completely different behavior, namely: the occupancy of the target site is always increased by higher molecular crowding on the DNA. Finally, we show that crowding on the DNA may increase transcriptional noise through increased variability of the occupancy time of the target sites.

Stopped-flow fluorescence kinetic study of protein sliding and intersegment transfer in the target DNA search process

Kinetic characterizations of protein translocation on DNA are nontrivial because the simultaneous presence of multiple different mechanisms makes it difficult to extract the information specific to a particular translocation mechanism. In this study, we have developed new approaches for the kinetic investigations of proteins' sliding and intersegment transfer (also known as "direct transfer") in the target DNA search process. Based on the analytical expression of the mean search time for the discrete-state stochastic model, we derived analytical forms of the apparent rate constant kapp for protein-target association in systems involving competitor DNA and the intersegment transfer mechanism. Our analytical forms of kapp facilitate the experimental determination of the kinetic rate constants for intersegment transfer and sliding in the target association process. Using stopped-flow fluorescence data for the target association kinetics along with the analytical forms of kapp, we have studied the translocation of the Egr-1 zinc-finger protein in the target DNA association process. Sliding was analyzed using the DNA-length-dependent kapp data. Using the dependence of kapp on the concentration of competitor DNA, we determined the second-order rate constant for intersegment transfer. Our results indicate that a major pathway in the target association process for the Egr-1 zinc-finger protein is the one involving intersegment transfer to a nonspecific site and the subsequent sliding to the target.

Asymmetric DNA-search dynamics by symmetric dimeric proteins

We focus on dimeric DNA-binding proteins from two well-studied families: orthodox type II restriction endonucleases (REs) and transcription factors (TFs). Interactions of the protein's recognition sites with the DNA and, particularly, the contribution of each of the monomers to one-dimensional (1D) sliding along nonspecific DNA were studied using computational tools. Coarse-grained molecular dynamics simulations of DNA scanning by various TFs and REs provide insights into how the symmetry of a homodimer can be broken while they nonspecifically interact with DNA. The characteristics of protein sliding along DNA, such as the average sliding length, partitioning between 1D and 3D search, and the one-dimensional diffusion coefficient D1, strongly depend on the salt concentration, which in turn affects the probability of the two monomers adopting a cooperative symmetric sliding mechanism. Indeed, we demonstrate that maximal DNA search efficiency is achieved when the protein adopts an asymmetric search mode in which one monomer slides while its partner hops. We find that proteins classified as TFs have a higher affinity for the DNA, longer sliding lengths, and an increased probability of symmetric sliding in comparison with REs. Moreover, TFs can perform their biological function over a much wider range of salt concentrations than REs. Our results demonstrate that the different biological functions of DNA-binding proteins are related to the different nonspecific DNA search mechanisms they adopt.

Extending single molecule fluorescence observation time by amplitude-modulated excitation

We present a hardware-based method that can improve single molecule fluorophore observation time by up to 1500% and super-localization by 47% for the experimental conditions used. The excitation was modulated using an acousto-optic modulator (AOM) synchronized to the data acquisition and inherent data conversion time of the detector. The observation time and precision in super-localization of four commonly used fluorophores were compared under modulated and traditional continuous excitation, including direct total internal reflectance excitation of Alexa 555 and Cy3, non-radiative Förster resonance energy transfer (FRET) excited Cy5, and direct epi-fluorescence wide field excitation of Rhodamine 6G. The proposed amplitude-modulated excitation does not perturb the chemical makeup of the system or sacrifice signal and is compatible with multiple types of fluorophores. Amplitude-modulated excitation has practical applications for any fluorescent study utilizing an instrumental setup with time-delayed detectors.

The helicase-like domains of type III restriction enzymes trigger long-range diffusion along DNA

Helicases are ubiquitous adenosine triphosphatases (ATPases) with widespread roles in genome metabolism. Here, we report a previously undescribed functionality for ATPases with helicase-like domains; namely, that ATP hydrolysis can trigger ATP-independent long-range protein diffusion on DNA in one dimension (1D). Specifically, using single-molecule fluorescence microscopy we show that the Type III restriction enzyme EcoP15I uses its ATPase to switch into a distinct structural state that diffuses on DNA over long distances and long times. The switching occurs only upon binding to the target site and requires hydrolysis of ~30 ATPs. We define the mechanism for these enzymes and show how ATPase activity is involved in DNA target site verification and 1D signaling, roles that are common in DNA metabolism: for example, in nucleotide excision and mismatch repair.

Proteins searching for their target on DNA by one-dimensional diffusion: overcoming the "speed-stability" paradox

The sequence dependence of DNA-protein interactions that allows proteins to find the correct reaction site also slows down the 1D diffusion of the protein along the DNA molecule, leading to the so-called "speed-stability paradox," wherein fast diffusion along the DNA molecule is seemingly incompatible with stable targeting of the reaction site. Here, we develop diffusion-reaction models that use discrete and continuous Gaussian random 1D diffusion landscapes with or without a high-energy cut-off, and two-state models with a transition to and from a "searching" mode in which the protein diffuses rapidly without recognizing the target. We show the conditions under which such considerations lead to a predicted speed-up of the targeting process, and under which the presence of a "searching" mode in a two-state model is nearly equivalent to the existence of a high-energy cut-off in a one-state model. We also determine the conditions under which the search is either diffusion-limited or reaction-limited, and develop quantitative expressions for the rate of successful targeting as a function of the site-specific reaction rate, the roughness of the DNA-protein interaction potential, and the presence of a "searching" mode. In general, we find that a rough landscape is compatible with a fast search if the highest energy barriers can be avoided by "hopping" or by the protein transitioning to a lower-energy "searching" mode. We validate these predictions with the results of Brownian dynamics, kinetic Metropolis, and kinetic Monte Carlo simulations of the diffusion and targeting process, and apply these concepts to the case of T7 RNA polymerase searching for its target site on T7 DNA.

Optical Methods to Study Protein-DNA Interactions in Vitro and in Living Cells at the Single-Molecule Level

The maintenance of intact genetic information, as well as the deployment of transcription for specific sets of genes, critically rely on a family of proteins interacting with DNA and recognizing specific sequences or features. The mechanisms by which these proteins search for target DNA are the subject of intense investigations employing a variety of methods in biology. A large interest in these processes stems from the faster-than-diffusion association rates, explained in current models by a combination of 3D and 1D diffusion. Here, we present a review of the single-molecule approaches at the forefront of the study of protein-DNA interaction dynamics and target search in vitro and in vivo. Flow stretch, optical and magnetic manipulation, single fluorophore detection and localization as well as combinations of different methods are described and the results obtained with these techniques are discussed in the framework of the current facilitated diffusion model.

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.

Roles for Helicases as ATP-Dependent Molecular Switches

On the basis of the familial name, a "helicase" might be expected to have an enzymatic activity that unwinds duplex polynucleotides to form single strands. A more encompassing taxonomy that captures alternative enzymatic roles has defined helicases as a sub-class of molecular motors that move directionally and processively along nucleic acids, the so-called "translocases". However, even this definition may be limiting in capturing the full scope of helicase mechanism and activity. Discussed here is another, alternative view of helicases-as machines which couple NTP-binding and hydrolysis to changes in protein conformation to resolve stable nucleoprotein assembly states. This "molecular switch" role differs from the classical view of helicases as molecular motors in that only a single catalytic NTPase cycle may be involved. This is illustrated using results obtained with the DEAD-box family of RNA helicases and with a model bacterial system, the ATP-dependent Type III restriction-modification enzymes. Further examples are discussed and illustrate the wide-ranging examples of molecular switches in genome metabolism.

RecA filament sliding on DNA facilitates homology search

During homologous recombination, RecA forms a helical filament on a single stranded (ss) DNA that searches for a homologous double stranded (ds) DNA and catalyzes the exchange of complementary base pairs to form a new heteroduplex. Using single molecule fluorescence imaging tools with high spatiotemporal resolution we characterized the encounter complex between the RecA filament and dsDNA. We present evidence in support of the 'sliding model' wherein a RecA filament diffuses along a dsDNA track. We further show that homology can be detected during sliding. Sliding occurs with a diffusion coefficient of approximately 8000 bp(2)/s allowing the filament to sample several hundred base pairs before dissociation. Modeling suggests that sliding can accelerate homology search by as much as 200 fold. Homology recognition can occur for as few as 6 nt of complementary basepairs with the recognition efficiency increasing for higher complementarity. Our data represents the first example of a DNA bound multi-protein complex which can slide along another DNA to facilitate target search.DOI:http://dx.doi.org/10.7554/eLife.00067.001.

AID and Apobec3G haphazard deamination and mutational diversity

Activation-induced deoxycytidine deaminase (AID) and Apobec 3G (Apo3G) cause mutational diversity by initiating mutations on regions of single-stranded (ss) DNA. Expressed in B cells, AID deaminates C → U in actively transcribed immunoglobulin (Ig) variable and switch regions to initiate the somatic hypermutation (SHM) and class switch recombination (CSR) that are essential for antibody diversity. Apo3G expressed in T cells catalyzes C deaminations on reverse transcribed cDNA causing HIV-1 retroviral inactivation. When operating properly, AID- and Apo3G-initiated mutations boost human fitness. Yet, both enzymes are potentially powerful somatic cell "mutators". Loss of regulated expression and proper genome targeting can cause human cancer. Here, we review well-established biological roles of AID and Apo3G. We provide a synopsis of AID partnering proteins during SHM and CSR, and describe how an Apo2 crystal structure provides "surrogate" insight for AID and Apo3G biochemical behavior. However, large gaps remain in our understanding of how dC deaminases search ssDNA to identify trinucleotide motifs to deaminate. We discuss two recent methods to analyze ssDNA scanning and deamination. Apo3G scanning and deamination is visualized in real-time using single-molecule FRET, and AID deamination efficiencies are determined with a random walk analysis. AID and Apo3G encounter many candidate deamination sites while scanning ssDNA. Generating mutational diversity is a principal aim of AID and an important ancillary property of Apo3G. Success seems likely to involve hit and miss deamination motif targeting, biased strongly toward miss.

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.