Sequence heterogeneity accelerates protein search for targets on DNA

How to scan naked DNA using promiscuous recognition and no clamping: a model for pioneer transcription factors

Most DNA scanning proteins uniquely recognize their cognate sequence motif and slide on DNA assisted by some sort of clamping interface. The pioneer transcription factors that control cell fate in eukaryotes must forgo both elements to gain access to DNA in naked and chromatin forms; thus, whether or how these factors scan naked DNA is unknown. Here, we use single-molecule techniques to investigate naked DNA scanning by the Engrailed homeodomain (enHD) as paradigm of highly promiscuous recognition and open DNA binding interface. We find that enHD scans naked DNA quite effectively, and about 200000-fold faster than expected for a continuous promiscuous slide. To do so, enHD scans about 675 bp of DNA in 100 ms and then redeploys stochastically to another location 530 bp afar in just 10 ms. During the scanning phase enHD alternates between slow- and medium-paced modes every 3 and 40 ms, respectively. We also find that enHD binds nucleosomes and does so with enhanced affinity relative to naked DNA. Our results demonstrate that pioneer-like transcription factors can in principle do both, target nucleosomes and scan active DNA efficiently. The hybrid scanning mechanism used by enHD appears particularly well suited for the highly complex genomic signals of eukaryotic cells.

Why Are Nucleosome Breathing Dynamics Asymmetric?

In eukaryotic cells, DNA is bound to nucleosomes, but DNA segments occasionally unbind in the process known as nucleosome breathing. Although DNA can unwrap simultaneously from both ends of the nucleosome (symmetric breathing), experiments indicate that DNA prefers to dissociate from only one end (asymmetric breathing). However, the molecular origin of the asymmetry is not understood. We developed a new theoretical approach that gives microscopic explanations of asymmetric breathing. It is based on a stochastic description that leads to a comprehensive evaluation of dynamics by using effective free-energy landscapes. It is shown that asymmetric breathing follows the kinetically preferred pathways. In addition, it is also found that asymmetric breathing leads to a faster target search by transcription factors. Theoretical predictions, supported by computer simulations, agree with experiments. It is proposed that nature utilizes the symmetry of nucleosome breathing to achieve a better dynamic accessibility of chromatin for more efficient genetic regulation.

How Do Nucleosome Dynamics Regulate Protein Search on DNA?

Nearly three-fourths of all eukaryotic DNA is occupied by nucleosomes, protein-DNA complexes comprising octameric histone core proteins and ∼150 base pairs of DNA. In addition to acting as a DNA compaction vehicle, the dynamics of nucleosomes regulate the DNA site accessibility for the nonhistone proteins, thereby controlling regulatory processes involved in determining the cell identity and cell fate. Here, we propose an analytical framework to analyze the role of nucleosome dynamics on the target search process of transcription factors through a simple discrete-state stochastic description of the search process. By considering the experimentally determined kinetic rates associated with protein and nucleosome dynamics as the only inputs, we estimate the target search time of a protein via first-passage probability calculations separately during nucleosome breathing and sliding dynamics. Although both the nucleosome dynamics permit transient access to the DNA sites that are otherwise occluded by the histone proteins, our result suggests substantial differences between the protein search mechanism on a nucleosome performing breathing and sliding dynamics. Furthermore, we identify the molecular factors that influence the search efficiency and demonstrate how these factors together portray a highly dynamic landscape of gene regulation. Our analytical results are validated using extensive Monte Carlo simulations.

Understanding the Target Search by Multiple Transcription Factors on Nucleosomal DNA

The association of multiple Transcription Factors (TFs) in the cis-regulatory region is imperative for developmental changes in eukaryotes. The underlying process is exceedingly complex, and it is not at all clear what orchestrates the overall search process by multiple TFs. In this study, by developing a theoretical model based on a discrete-state stochastic approach, we investigated the target search mechanism of multiple TFs on nucleosomal DNA. Experimental kinetic rate constants of different TFs are taken as input to estimate the Mean-First-Passage time to recognize the binding motifs by two TFs on a dynamic nucleosome model. The theory systematically analyzes when the TFs search their binding motifs hierarchically and when simultaneously by proceeding via the formation of a protein-protein complex. Our results, validated by extensive Monte Carlo simulations, elucidate the molecular basis of the complex target search phenomenon of multiple TFs on nucleosomal DNA.

How Pioneer Transcription Factors Search for Target Sites on Nucleosomal DNA

All major biological processes start after protein molecules known as transcription factors detect specific regulatory sequences on DNA and initiate genetic expression by associating to them. But in eukaryotic cells, much of the DNA is covered by nucleosomes and other chromatin structures, preventing transcription factors from binding to their targets. At the same time, experimental studies show that there are several classes of proteins, called "pioneer transcription factors", that are able to reach the targets on nucleosomal DNA; however, the underlying microscopic mechanisms remain not well understood. We propose a new theoretical approach that might explain how pioneer transcription factors can find their targets. It is argued that pioneer transcription factors might weaken the interactions between the DNA and nucleosome by substituting them with similar interactions between transcription factors and DNA. Using this idea, we develop a discrete-state stochastic model that allows for exact calculations of target search dynamics on nucleosomal DNA using first-passage probabilities approach. It is found that the target search on nuclesomal DNA for pioneer transcription factors might be significantly accelerated while the search is slower on naked DNA in comparison with normal transcription factors. Our theoretical predictions are supported by Monte Carlo computer simulations, and they also agree with available experimental observations.

Dynamics of the Protein Search for Targets on DNA in Quorum Sensing Cells

Quorum sensing (QS) is a bacterial cell-cell communication process that regulates gene expression. The search and binding of the AHL-bound LuxR-type proteins to specific sites on DNA in quorum sensing cells in Gram-negative bacteria is a complex process and has been theoretically investigated based on a discrete-state stochastic approach. It is shown that several factors such as the rate of formation of the AHL-bound LuxR protein within the cells and its dissociation to freely diffusing autoinducer molecule (AHL), the diffusion of the latter in and out of the cells, positive feedback loops, and the cell population density play an important role in the protein target search and can control the gene regulation processes. Physical-chemical arguments to explain these observations are presented. Our calculations of the dynamic properties are also supplemented by Monte Carlo computer simulations. Our theoretical model provides physical insights into the complex mechanisms of protein target search in quorum sensing cells.

Influence of Nonspecific Interactions between Proteins and In Vivo Cytoplasmic Crowders in Facilitated Diffusion of Proteins: Theoretical Insights

The binding of proteins to their respective specific sites on the DNA through facilitated diffusion serves as the initial step of various important biological processes. While this search process has been thoroughly investigated via in vitro studies, the cellular environment is complex and may interfere with the protein's search dynamics. The cytosol is heavily crowded, which can potentially modify the search by nonspecifically interacting with the protein that has been mostly overlooked. In this work, we probe the target search dynamics in the presence of explicit crowding agents that have an affinity toward the protein. We theoretically investigate the role of such protein-crowder associations in the target search process using a discrete-state stochastic framework that allows for the analytical description of dynamic properties. It is found that stronger nonspecific associations between the crowder and proteins can accelerate the facilitated diffusion of proteins in comparison with a purely inert, rather weakly interacting cellular environment. This effect depends on how strong these associations are, the spatial positions of the target with respect to the crowders, and the size of the crowded region. Our theoretical results are also tested with Monte Carlo computer simulations. Our predictions are in qualitative agreement with existing experimental observations and computational studies.

Kinetic origin of nucleosome invasion by pioneer transcription factors

Recently, a cryo-electron microscopy study has captured different stages of nucleosome breathing dynamics that show partial unwrapping of DNA from histone core to permit transient access to the DNA sites by transcription factors. In practice, however, only a subset of transcription factors named pioneer factors can invade nucleosomes and bind to specific DNA sites to trigger essential DNA metabolic processes. We propose a discrete-state stochastic model that considers the interplay of nucleosome breathing and protein dynamics explicitly and estimate the mean time to search the target DNA sites. It is found that the molecular principle governing the search process on nucleosome is very different compared to that on naked DNA. The pioneer factors minimize their search times on nucleosomal DNA by compensating their nucleosome association rates by dissociation rates. A fine balance between the two presents a tradeoff between their nuclear mobility and error associated with the search process.

DNA Looping and DNA Conformational Fluctuations Can Accelerate Protein Target Search

Protein searching and binding to specific sites on DNA is a fundamentally important process that marks the beginning of all major cellular transformations. While the dynamics of protein-DNA interactions in in vitro settings is well investigated, the situation is much more complex for in vivo conditions because the DNA molecules in live cells are packed into chromosomal structures where they are undergoing strong dynamic and conformational fluctuations. In this work, we present a theoretical investigation on the role of DNA looping and DNA conformational fluctuations in the protein target search. It is based on a discrete-state stochastic analysis that allows for explicit calculations of dynamic properties, which is also supplemented by Monte Carlo computer simulations. It is found that for stronger nonspecific interactions between DNA and proteins the search occurs faster on the DNA looped conformation in comparison with the unlooped conformation, and the fastest search is observed when the loop is formed near the target site. It is also shown that DNA fluctuations between the looped and unlooped conformations influence the search dynamics, and this depends on the magnitude of conformational transition rates and on which conformation is more energetically stable. Physical-chemical arguments explaining these observations are presented. Our theoretical study suggests that the geometry and conformational changes in DNA are additional factors that might efficiently control the gene regulation processes.

Discrete-state stochastic kinetic models for target DNA search by proteins: Theory and experimental applications

To perform their functions, transcription factors and DNA-repair/modifying enzymes randomly search DNA in order to locate their specific targets on DNA. Discrete-state stochastic kinetic models have been developed to explain how the efficiency of the search process is influenced by the molecular properties of proteins and DNA as well as by other factors such as molecular crowding. These theoretical models not only offer explanations on the relation of microscopic processes to macroscopic behavior of proteins, but also facilitate the analysis and interpretation of experimental data. In this review article, we provide an overview on discrete-state stochastic kinetic models and explain how these models can be applied to experimental investigations using stopped-flow, single-molecule, nuclear magnetic resonance (NMR), and other biophysical and biochemical methods.

Effect of DNA Conformation on the Protein Search for Targets on DNA: A Theoretical Perspective

Protein-DNA interactions are important for all biological processes and involve the process of proteins searching for and recognizing specific binding sites on the DNA. Many aspects of the mechanism of the protein search for targets on DNA are not well understood. One important problem is the effect of DNA conformation on the protein search dynamics. Using a theoretical method based on a discrete-state stochastic approach, we obtained an analytical description of the dynamic properties. We investigated a system with two DNA conformers. It was found that the average search time on one DNA conformer via 1D or 3D motions depended on the dynamics of the search process on the other DNA conformer.

The effect of obstacles in multi-site protein target search with DNA looping

Many fundamental biological processes are regulated by protein-DNA complexes called synaptosomes, which possess multiple interaction sites. Despite the critical importance of synaptosomes, the mechanisms of their formation are not well understood. Because of the multisite nature of participating proteins, it is widely believed that their search for specific sites on DNA involves the formation and breaking of DNA loops and sliding in the looped configurations. In reality, DNA in live cells is densely covered by other biological molecules that might interfere with the formation of synaptosomes. In this work, we developed a theoretical approach to evaluate the role of obstacles in the target search of multisite proteins when the formation of DNA loops and the sliding in looped configurations are possible. Our theoretical method is based on analysis of a discrete-state stochastic model that uses a master equations approach and extensive computer simulations. It is found that the obstacle slows down the search dynamics in the system when DNA loops are long-lived, but the effect is minimal for short-lived DNA loops. In addition, the relative positions of the target and the obstacle strongly influence the target search kinetics. Furthermore, the presence of the obstacle might increase the noise in the system. These observations are discussed using physical-chemical arguments. Our theoretical approach clarifies the molecular mechanisms of formation of protein-DNA complexes with multiple interactions sites.

Eukaryotic transcription factors can track and control their target genes using DNA antennas

Eukaryotic transcription factors (TF) function by binding to short 6-10 bp DNA recognition sites located near their target genes, which are scattered through vast genomes. Such process surmounts enormous specificity, efficiency and celerity challenges using a molecular mechanism that remains poorly understood. Combining biophysical experiments, theory and bioinformatics, we dissect the interplay between the DNA-binding domain of Engrailed, a Drosophila TF, and the regulatory regions of its target genes. We find that Engrailed binding affinity is strongly amplified by the DNA regions flanking the recognition site, which contain long tracts of degenerate recognition-site repeats. Such DNA organization operates as an antenna that attracts TF molecules in a promiscuous exchange among myriads of intermediate affinity binding sites. The antenna ensures a local TF supply, enables gene tracking and fine control of the target site's basal occupancy. This mechanism illuminates puzzling gene expression data and suggests novel engineering strategies to control gene expression.

Facilitation of DNA loop formation by protein-DNA non-specific interactions

Complex DNA topological structures, including polymer loops, are frequently observed in biological processes when protein molecules simultaneously bind to several distant sites on DNA. However, the molecular mechanisms of formation of these systems remain not well understood. Existing theoretical studies focus only on specific interactions between protein and DNA molecules at target sequences. However, the electrostatic origin of primary protein-DNA interactions suggests that interactions of proteins with all DNA segments should be considered. Here we theoretically investigate the role of non-specific interactions between protein and DNA molecules on the dynamics of loop formation. Our approach is based on analyzing a discrete-state stochastic model via a method of first-passage probabilities supplemented by Monte Carlo computer simulations. It is found that depending on a protein sliding length during the non-specific binding event three different dynamic regimes of the DNA loop formation might be observed. In addition, the loop formation time might be optimized by varying the protein sliding length, the size of the DNA molecule, and the position of the specific target sequences on DNA. Our results demonstrate the importance of non-specific protein-DNA interactions in the dynamics of DNA loop formations.

Single-molecule DNA unzipping reveals asymmetric modulation of a transcription factor by its binding site sequence and context

Most functional transcription factor (TF) binding sites deviate from their ‘consensus’ recognition motif, although their sites and flanking sequences are often conserved across species. Here, we used single-molecule DNA unzipping with optical tweezers to study how Egr-1, a TF harboring three zinc fingers (ZF1, ZF2 and ZF3), is modulated by the sequence and context of its functional sites in the Lhb gene promoter. We find that both the core 9 bp bound to Egr-1 in each of the sites, and the base pairs flanking them, modulate the affinity and structure of the protein–DNA complex. The effect of the flanking sequences is asymmetric, with a stronger effect for the sequence flanking ZF3. Characterization of the dissociation time of Egr-1 revealed that a local, mechanical perturbation of the interactions of ZF3 destabilizes the complex more effectively than a perturbation of the ZF1 interactions. Our results reveal a novel role for ZF3 in the interaction of Egr-1 with other proteins and the DNA, providing insight on the regulation of Lhb and other genes by Egr-1. Moreover, our findings reveal the potential of small changes in DNA sequence to alter transcriptional regulation, and may shed light on the organization of regulatory elements at promoters.

Molecular search with conformational change: One-dimensional discrete-state stochastic model

Molecular search phenomena are observed in a variety of chemical and biological systems. During the search, the participating particles frequently move in complex inhomogeneous environments with random transitions between different dynamic modes. To understand the mechanisms of molecular search with alternating dynamics, we investigate the search dynamics with stochastic transitions between two conformations in a one-dimensional discrete-state stochastic model. It is explicitly analyzed using the first-passage time probability method to obtain a full dynamic description of the search process. A general dynamic phase diagram is developed. It is found that there are several dynamic regimes in the molecular search with conformational transitions, and they are determined by the relative values of the relevant length scales in the system. Theoretical predictions are fully supported by Monte Carlo computer simulations.

Mechanisms of Protein Search for Targets on DNA: Theoretical Insights

Protein-DNA interactions are critical for the successful functioning of all natural systems. The key role in these interactions is played by processes of protein search for specific sites on DNA. Although it has been studied for many years, only recently microscopic aspects of these processes became more clear. In this work, we present a review on current theoretical understanding of the molecular mechanisms of the protein target search. A comprehensive discrete-state stochastic method to explain the dynamics of the protein search phenomena is introduced and explained. Our theoretical approach utilizes a first-passage analysis and it takes into account the most relevant physical-chemical processes. It is able to describe many fascinating features of the protein search, including unusually high effective association rates, high selectivity and specificity, and the robustness in the presence of crowders and sequence heterogeneity.

Transcription pausing: biological significance of thermal fluctuations biased by repetitive genomic sequences

Transcription of DNA by RNA polymerase (RNAP) takes place in a cell environment dominated by thermal fluctuations. How are transcription reactions including initiation, elongation, and termination on genomic DNA so well-controlled during such fluctuations? A recent statistical mechanical approach using high-throughput sequencing data reveals that repetitive DNA sequence elements embedded into a genomic sequence provide the key mechanism to functionally bias the fluctuations of transcription elongation complexes. In particular, during elongation pausing, such repetitive sequence elements can increase the magnitude of one-dimensional diffusion of the RNAP enzyme on the DNA upstream of the pausing site, generating a large variation in the dwell times of RNAP pausing under the control of these genomic signals.

Mechanism of Genome Interrogation: How CRISPR RNA-Guided Cas9 Proteins Locate Specific Targets on DNA

The ability to precisely edit and modify a genome opens endless opportunities to investigate fundamental properties of living systems as well as to advance various medical techniques and bioengineering applications. This possibility is now close to reality due to a recent discovery of the adaptive bacterial immune system, which is based on clustered regularly interspaced short palindromic repeats (CRISPR)-associated proteins (Cas) that utilize RNA to find and cut the double-stranded DNA molecules at specific locations. Here we develop a quantitative theoretical approach to analyze the mechanism of target search on DNA by CRISPR RNA-guided Cas9 proteins, which is followed by a selective cleavage of nucleic acids. It is based on a discrete-state stochastic model that takes into account the most relevant physical-chemical processes in the system. Using a method of first-passage processes, a full dynamic description of the target search is presented. It is found that the location of specific sites on DNA by CRISPR Cas9 proteins is governed by binding first to protospacer adjacent motif sequences on DNA, which is followed by reversible transitions into DNA interrogation states. In addition, the search dynamics is strongly influenced by the off-target cutting. Our theoretical calculations allow us to explain the experimental observations and to give experimentally testable predictions. Thus, the presented theoretical model clarifies some molecular aspects of the genome interrogation by CRISPR RNA-guided Cas9 proteins.

Optimal Length of Conformational Transition Region in Protein Search for Targets on DNA

The starting point of many fundamental biological processes is associated with protein molecules finding and recognizing specific sites on DNA. However, despite a large number of experimental and theoretical studies on protein search for targets on DNA, many molecular aspects of underlying mechanisms are still not well understood. Experiments show that proteins bound to DNA can switch between slow recognition and fast search conformations. However, from a theoretical point of view, such conformational transitions should slow down the protein search for specific sites on DNA, in contrast to available experimental observations. In addition, experiments indicate that the nucleotide composition near the target site is more symmetrically homogeneous, leading to stronger effective interactions between proteins and DNA at these locations. However, as has been shown theoretically, this should also make the search less efficient, which is not observed. We propose a possible resolution of these problems by suggesting that conformational transitions occur only within a segment around the target where stronger interactions between proteins and DNA are observed. Two theoretical methods, based on continuum and discrete-state stochastic calculations, are developed, allowing us to obtain a comprehensive dynamic description for the protein search process in this system. The existence of an optimal length of the conformational transition zone with the shortest mean search time is predicted.

On the Mechanism of Homology Search by RecA Protein Filaments

Genetic stability is a key factor in maintaining, survival, and reproduction of biological cells. It relies on many processes, but one of the most important is a homologous recombination, in which the repair of breaks in double-stranded DNA molecules is taking place with a help of several specific proteins. In bacteria, this task is accomplished by RecA proteins that are active as nucleoprotein filaments formed on single-stranded segments of DNA. A critical step in the homologous recombination is a search for a corresponding homologous region on DNA, which is called a homology search. Recent single-molecule experiments clarified some aspects of this process, but its molecular mechanisms remain not well understood. We developed a quantitative theoretical approach to analyze the homology search. It is based on a discrete-state stochastic model that takes into account the most relevant physical-chemical processes in the system. Using a method of first-passage processes, a full dynamic description of the homology search is presented. It is found that the search dynamics depends on the degree of extension of DNA molecules and on the size of RecA nucleoprotein filaments, in agreement with experimental single-molecule measurements of DNA pairing by RecA proteins. Our theoretical calculations, supported by extensive Monte Carlo computer simulations, provide a molecular description of the mechanisms of the homology search.

The Role of DNA Looping in the Search for Specific Targets on DNA by Multisite Proteins

Many cellular processes involve simultaneous interactions between DNA and protein molecules at several locations. They are regulated and controlled by special protein-DNA complexes, which are known as synaptic complexes or synaptosomes. Because of the multisite nature of involved proteins, it was suggested that during the formation of synaptic complexes DNA loops might appear, but their role is unclear. We developed a theoretical model that allowed us to evaluate the effect of transient DNA loop formation. It is based on a discrete-state stochastic method that explicitly takes into account the free-energy contributions due to the appearance of DNA loops. The formation of the synaptic complexes is viewed as a search for a specific binding site on DNA by the protein molecule already bound to DNA at another location. It was found that the search might be optimized by varying the position of the target and the total length of DNA. Furthermore, the formation of transient DNA loops leads to faster dynamics if it is associated with favorable enthalpic contributions to nonspecific protein-DNA interactions. It is also shown that DNA looping might reduce stochastic noise in the system.

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.