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

Analyzing paramagnetic NMR data on target DNA search by proteins using a discrete-state kinetic model for translocation

Before reaching their targets, sequence-specific DNA-binding proteins nonspecifically bind to DNA through electrostatic interactions and stochastically change their locations on DNA. Investigations into the dynamics of DNA-scanning by proteins are nontrivial due to the simultaneous presence of multiple translocation mechanisms and many sites for the protein to nonspecifically bind to DNA. Nuclear magnetic resonance (NMR) spectroscopy can provide information about the target DNA search processes at an atomic level. Paramagnetic relaxation enhancement (PRE) is particularly useful to study how the proteins scan DNA in the search process. Previously, relatively simple two-state or three-state exchange models were used to explain PRE data reflecting the target search process. In this work, using more realistic discrete-state stochastic kinetics models embedded into an NMR master equation, we analyzed the PRE data for the HoxD9 homeodomain interacting with DNA. The kinetic models that incorporate sliding, dissociation, association, and intersegment transfer can reproduce the PRE profiles observed at some different ionic strengths. The analysis confirms the previous interpretation of the PRE data and shows that the protein's probability distribution among nonspecific sites is nonuniform during the target DNA search process.

Spatio-temporal dynamics of the DNA glycosylase OGG1 in finding and processing 8-oxoguanine

OGG1 is the DNA glycosylase responsible for the removal of the oxidative lesion 8-oxoguanine (8-oxoG) from DNA. The recognition of this lesion by OGG1 is a complex process that involves scanning the DNA for the presence of 8-oxoG, followed by recognition and lesion removal. Structural data have shown that OGG1 evolves through different stages of conformation onto the DNA, corresponding to elementary steps of the 8-oxoG recognition and extrusion from the double helix. Single-molecule studies of OGG1 on naked DNA have shown that OGG1 slides in persistent contact with the DNA, displaying different binding states probably corresponding to the different conformation stages. However, in cells, the DNA is not naked and OGG1 has to navigate into a complex and highly crowded environment within the nucleus. To ensure rapid detection of 8-oxoG, OGG1 alternates between 3D diffusion and sliding along the DNA. This process is regulated by the local chromatin state but also by protein co-factors that could facilitate the detection of oxidized lesions. We will review here the different methods that have been used over the last years to better understand how OGG1 detects and process 8-oxoG lesions.

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

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

Diverse role of decoys on emergence and precision of oscillations in a biomolecular clock

Biomolecular clocks are key drivers of oscillatory dynamics in diverse biological processes including cell-cycle regulation, circadian rhythms, and pattern formation during development. A minimal clock implementation is based on the classical Goodwin oscillator, in which a repressor protein inhibits its own synthesis via time-delayed negative feedback. Clock motifs, however, do not exist in isolation; its components are open to interacting with the complex environment inside cells. For example, there are ubiquitous high-affinity binding sites along the genome, known as decoys, where transcription factors such as repressor proteins can potentially interact. This binding affects the availability of transcription factors and has often been ignored in theoretical studies. How does such genomic decoy binding impact the clock's robustness and precision? To address this question, we systematically analyze deterministic and stochastic models of the Goodwin oscillator in the presence of reversible binding of the repressor to a finite number of decoy sites. Our analysis reveals that the relative stability of decoy-bound repressors compared to the free repressor plays distinct roles on the emergence and precision of oscillations. Interestingly, active degradation of the bound repressor can induce sustained oscillations that are otherwise absent without decoys. In contrast, decoy abundances can kill oscillation dynamics if the bound repressor is protected from degradation. Taking into account low copy-number fluctuations in clock components, we show that the degradation of the bound repressors enhances precision by attenuating noise in both the amplitude and period of oscillations. Overall, these results highlight the versatile role of otherwise hidden decoys in shaping the stochastic dynamics of biological clocks and emphasize the importance of synthetic decoys in designing robust clocks.

Single-Molecule Microscopy Meets Molecular Dynamics Simulations for Characterizing the Molecular Action of Proteins on DNA and in Liquid Condensates

DNA-binding proteins trigger various cellular functions and determine cellular fate. Before performing functions such as transcription, DNA repair, and DNA recombination, DNA-binding proteins need to search for and bind to their target sites in genomic DNA. Under evolutionary pressure, DNA-binding proteins have gained accurate and rapid target search and binding strategies that combine three-dimensional search in solution, one-dimensional sliding along DNA, hopping and jumping on DNA, and intersegmental transfer between two DNA molecules. These mechanisms can be achieved by the unique structural and dynamic properties of these proteins. Single-molecule fluorescence microscopy and molecular dynamics simulations have characterized the molecular actions of DNA-binding proteins in detail. Furthermore, these methodologies have begun to characterize liquid condensates induced by liquid-liquid phase separation, e.g., molecular principles of uptake and dynamics in droplets. This review discusses the molecular action of DNA-binding proteins on DNA and in liquid condensate based on the latest studies that mainly focused on the model protein p53.

Testing mechanisms of DNA sliding by architectural DNA-binding proteins: dynamics of single wild-type and mutant protein molecules in vitro and in vivo

Architectural DNA-binding proteins (ADBPs) are abundant constituents of eukaryotic or bacterial chromosomes that bind DNA promiscuously and function in diverse DNA reactions. They generate large conformational changes in DNA upon binding yet can slide along DNA when searching for functional binding sites. Here we investigate the mechanism by which ADBPs diffuse on DNA by single-molecule analyses of mutant proteins rationally chosen to distinguish between rotation-coupled diffusion and DNA surface sliding after transient unbinding from the groove(s). The properties of yeast Nhp6A mutant proteins, combined with molecular dynamics simulations, suggest Nhp6A switches between two binding modes: a static state, in which the HMGB domain is bound within the minor groove with the DNA highly bent, and a mobile state, where the protein is traveling along the DNA surface by means of its flexible N-terminal basic arm. The behaviors of Fis mutants, a bacterial nucleoid-associated helix-turn-helix dimer, are best explained by mobile proteins unbinding from the major groove and diffusing along the DNA surface. Nhp6A, Fis, and bacterial HU are all near exclusively associated with the chromosome, as packaged within the bacterial nucleoid, and can be modeled by three diffusion modes where HU exhibits the fastest and Fis the slowest diffusion.

A Bit Stickier, a Bit Slower, a Lot Stiffer: Specific vs. Nonspecific Binding of Gal4 to DNA

Transcription factors regulate gene activity by binding specific regions of genomic DNA thanks to a subtle interplay of specific and nonspecific interactions that is challenging to quantify. Here, we exploit Reflective Phantom Interface (RPI), a label-free biosensor based on optical reflectivity, to investigate the binding of the N-terminal domain of Gal4, a well-known gene regulator, to double-stranded DNA fragments containing or not its consensus sequence. The analysis of RPI-binding curves provides interaction strength and kinetics and their dependence on temperature and ionic strength. We found that the binding of Gal4 to its cognate site is stronger, as expected, but also markedly slower. We performed a combined analysis of specific and nonspecific binding—equilibrium and kinetics—by means of a simple model based on nested potential wells and found that the free energy gap between specific and nonspecific binding is of the order of one kcal/mol only. We investigated the origin of such a small value by performing all-atom molecular dynamics simulations of Gal4–DNA interactions. We found a strong enthalpy–entropy compensation, by which the binding of Gal4 to its cognate sequence entails a DNA bending and a striking conformational freezing, which could be instrumental in the biological function of Gal4.

Nucleic acid binding by SAMHD1 contributes to the antiretroviral activity and is enhanced by the GpsN modification

SAMHD1 impedes infection of myeloid cells and resting T lymphocytes by retroviruses, and the enzymatic activity of the protein-dephosphorylation of deoxynucleotide triphosphates (dNTPs)-implicates enzymatic dNTP depletion in innate antiviral immunity. Here we show that the allosteric binding sites of the enzyme are plastic and can accommodate oligonucleotides in place of the allosteric activators, GTP and dNTP. SAMHD1 displays a preference for oligonucleotides containing phosphorothioate bonds in the Rp configuration located 3' to G nucleotides (GpsN), the modification pattern that occurs in a mechanism of antiviral defense in prokaryotes. In the presence of GTP and dNTPs, binding of GpsN-containing oligonucleotides promotes formation of a distinct tetramer with mixed occupancy of the allosteric sites. Mutations that impair formation of the mixed-occupancy complex abolish the antiretroviral activity of SAMHD1, but not its ability to deplete dNTPs. The findings link nucleic acid binding to the antiretroviral activity of SAMHD1, shed light on the immunomodulatory effects of synthetic phosphorothioated oligonucleotides and raise questions about the role of nucleic acid phosphorothioation in human innate immunity.

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.

How Retroactivity Affects the Behavior of Incoherent Feedforward Loops

An incoherent feedforward loop (IFFL) is a network motif known for its ability to accelerate responses and generate pulses. It remains an open question to understand the behavior of IFFLs in contexts with high levels of retroactivity, where an upstream transcription factor binds to numerous downstream binding sites. Here we study the behavior of IFFLs by simulating and comparing ODE models with different levels of retroactivity. We find that increasing retroactivity in an IFFL can increase, decrease, or keep the network's response time and pulse amplitude constant. This suggests that increasing retroactivity, traditionally considered an impediment to designing robust synthetic systems, could be exploited to improve the performance of IFFLs. In contrast, we find that increasing retroactivity in a negative autoregulated circuit can only slow the response. The ability of an IFFL to flexibly handle retroactivity may have contributed to its significant abundance in both bacterial and eukaryotic regulatory networks.

Racemic phosphorothioate as a tool for NMR investigations of protein-DNA complexes

A major driving force for protein-nucleic acid association is electrostatic interactions via ion pairs of the positively charged basic side chains and negatively charged phosphates. For a better understanding of how proteins scan DNA and recognize particular signatures, it is important to gain atomic-level insight into the behavior of basic side chains at the protein-DNA interfaces. NMR spectroscopy is a powerful tool for investigating the structural, dynamic, and kinetic aspects of protein-DNA interactions. However, resonance assignment of basic side-chain cationic moieties at the molecular interfaces remains to be a major challenge. Here, we propose a fast, robust, and inexpensive approach that greatly facilitates resonance assignment of interfacial moieties and also allows for kinetic measurements of protein translocation between two DNA duplexes. This approach utilizes site-specific incorporation of racemic phosphorothioate at the position of a phosphate that interacts with a protein side chain. This modification retains the electric charge of phosphate and therefore is mild, but causes significant chemical shift perturbations for the proximal protein side chains, which facilitates resonance assignment. Due to the racemic nature of the modification, two different chemical shifts are observed for the species with different diastereomers R_P and S_P of the incorporated phosphorothioate group. Kinetic information on the exchange of the protein molecule between R_P and S_P DNA duplexes can be obtained by ¹⁵N_z exchange spectroscopy. We demonstrate the applications of this approach to the Antennapedia homeodomain-DNA complex and the CREB1 basic leucine-zipper (bZIP)-DNA complex.

Enhancement of gene expression noise from transcription factor binding to genomic decoy sites

The genome contains several high-affinity non-functional binding sites for transcription factors (TFs) creating a hidden and unexplored layer of gene regulation. We investigate the role of such “decoy sites” in controlling noise (random fluctuations) in the level of a TF that is synthesized in stochastic bursts. Prior studies have assumed that decoy-bound TFs are protected from degradation, and in this case decoys function to buffer noise. Relaxing this assumption to consider arbitrary degradation rates for both bound/unbound TF states, we find rich noise behaviors. For low-affinity decoys, noise in the level of unbound TF always monotonically decreases to the Poisson limit with increasing decoy numbers. In contrast, for high-affinity decoys, noise levels first increase with increasing decoy numbers, before decreasing back to the Poisson limit. Interestingly, while protection of bound TFs from degradation slows the time-scale of fluctuations in the unbound TF levels, the decay of bound TFs leads to faster fluctuations and smaller noise propagation to downstream target proteins. In summary, our analysis reveals stochastic dynamics emerging from nonspecific binding of TFs and highlights the dual role of decoys as attenuators or amplifiers of gene expression noise depending on their binding affinity and stability of the bound TF.

Transcription Factor Binding in Embryonic Stem Cells Is Constrained by DNA Sequence Repeat Symmetry

Transcription factor (TF) recognition is dictated by the underlying DNA motif sequence specific for each TF. Here, we reveal that DNA sequence repeat symmetry plays a central role in defining TF-DNA-binding preferences. In particular, we find that different TFs bind similar symmetry patterns in the context of different developmental layers. Most TFs possess dominant preferences for similar DNA repeat symmetry types. However, in some cases, preferences of specific TFs are changed during differentiation, suggesting the importance of information encoded outside of known motif regions. Histone modifications also exhibit strong preferences for similar DNA repeat symmetry patterns unique to each type of modification. Next, using an in vivo reporter assay, we show that gene expression in embryonic stem cells can be positively modulated by the presence of genomic and computationally designed DNA oligonucleotides containing identified nonconsensus-repetitive sequence elements. This supports the hypothesis that certain nonconsensus-repetitive patterns possess a functional ability to regulate gene expression. We also performed a solution NMR experiment to probe the stability of double-stranded DNA via imino proton resonances for several double-stranded DNA sequences characterized by different repetitive patterns. We suggest that such local stability might play a key role in determining TF-DNA binding preferences. Overall, our findings show that despite the enormous sequence complexity of the TF-DNA binding landscape in differentiating embryonic stem cells, this landscape can be quantitatively characterized in simple terms using the notion of DNA sequence repeat symmetry.

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.

How p53 Molecules Solve the Target DNA Search Problem: A Review

Interactions between DNA and DNA-binding proteins play an important role in many essential cellular processes. A key function of the DNA-binding protein p53 is to search for and bind to target sites incorporated in genomic DNA, which triggers transcriptional regulation. How do p53 molecules achieve "rapid" and "accurate" target search in living cells? The search dynamics of p53 were expected to include 3D diffusion in solution, 1D diffusion along DNA, and intersegmental transfer between two different DNA strands. Single-molecule fluorescence microscopy enabled the tracking of p53 molecules on DNA and the characterization of these dynamics quantitatively. Recent intensive single-molecule studies of p53 succeeded in revealing each of these search dynamics. Here, we review these studies and discuss the target search mechanisms of p53.

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.

Structural Basis of Enhanced Facilitated Diffusion of DNA-Binding Protein in Crowded Cellular Milieu

Although the fast association between DNA-binding proteins (DBPs) and DNA is explained by a facilitated diffusion mechanism, in which DBPs adopt a weighted combination of three-dimensional diffusion and one-dimensional (1D) sliding and hopping modes of transportation, the role of cellular environment that contains many nonspecifically interacting proteins and other biomolecules is mostly overlooked. By performing large-scale computational simulations with an appropriately tuned model of protein and DNA in the presence of nonspecifically interacting bulk and DNA-bound crowders (genomic crowders), we demonstrate the structural basis of the enhanced facilitated diffusion of DBPs inside a crowded cellular milieu through, to our knowledge, novel 1D scanning mechanisms. In this one-dimensional scanning mode, the protein can float along the DNA under the influence of nonspecific interactions of bulk crowder molecules. The search mode is distinctly different compared to usual 1D sliding and hopping dynamics in which protein diffusion is regulated by the DNA electrostatics. In contrast, the presence of genomic crowders expedites the target search process by transporting the protein over DNA segments through the formation of a transient protein-crowder bridged complex. By analyzing the ruggedness of the associated potential energy landscape, we underpin the molecular origin of the kinetic advantages of these search modes and show that they successfully explain the experimentally observed acceleration of facilitated diffusion of DBPs by molecular crowding agents and crowder-concentration-dependent enzymatic activity of transcription factors. Our findings provide crucial insights into gene regulation kinetics inside the crowded cellular milieu.

Rational design of a DNA sequence-specific modular protein tag by tuning the alkylation kinetics

Sequence-selective chemical modification of DNA by synthetic ligands has been a long-standing challenge in the field of chemistry. Even when the ligand consists of a sequence-specific DNA binding domain and reactive group, sequence-selective reactions by these ligands are often accompanied by off-target reactions. A basic principle to design DNA modifiers that react at specific sites exclusively governed by DNA sequence recognition remains to be established. We have previously reported selective DNA modification by a self-ligating protein tag conjugated with a DNA-binding domain, termed as a modular adaptor, and orthogonal application of modular adaptors by relying on the chemoselectivity of the protein tag. The sequence-specific crosslinking reaction by the modular adaptor is thought to proceed in two steps: the first step involves the formation of a DNA-protein complex, while in the second step, a proximity-driven intermolecular crosslinking occurs. According to this scheme, the specific crosslinking reaction of a modular adaptor would be driven by the DNA recognition process only when the dissociation rate of the DNA complex is much higher than the rate constant for the alkylation reaction. In this study, as a proof of principle, a set of combinations for modular adaptors and their substrates were utilized to evaluate the reactions. Three types of modular adaptors consisting of a single type of self-ligating tag and three types of DNA binding proteins fulfill the kinetic requirements for the reaction of the self-ligating tag with a substrate and the dissociation of the DNA-protein complex. These modular adaptors actually undergo sequence-specific crosslinking reactions exclusively driven by the recognition of a specific DNA sequence. The design principle of sequence-specific modular adaptors based on the kinetic aspects of complex formation and chemical modification is applicable for developing recognition-driven selective modifiers for proteins and other biological macromolecules.

Intrinsically disordered domain of tumor suppressor p53 facilitates target search by ultrafast transfer between different DNA strands

Intersegmental transfer (IST) is an important strategy in the target search used by sequence-specific DNA-binding proteins (DBPs), enabling DBPs to search for targets between multiple DNA strands without dissociation. We examined the IST of the tumor suppressor p53 using ensemble stopped-flow and single-molecule fluorescence measurements. The ensemble measurements demonstrated that p53 exhibits very fast IST, whose rate constant was ∼108 M-1 s-1. To determine the domains of p53 responsible for IST, two mutants with deletions of one of its two DNA binding domains were generated. The mutant lacking the disordered C-terminal (CT) domain (the CoreTet mutant) abolished IST, whereas the mutant lacking the structured core domain (the TetCT mutant) maintained IST, clearly demonstrating the importance of the CT domain. Single-molecule fluorescence measurements further demonstrated the transfer of p53 between two tethered DNA strands. The pseudo-wild-type p53 and the TetCT mutant showed significant transfer efficiencies, whereas the transfer efficiency for the CoreTet mutant was zero. These results suggest that ultrafast IST might be promoted by four copies of the CT domain, by binding to two DNA strands simultaneously. Such ultrafast IST might be important to avoid nearby-bound DBPs during the target search process of p53 in nucleus.

Role for first zinc finger of WT1 in DNA sequence specificity: Denys-Drash syndrome-associated WT1 mutant in ZF1 enhances affinity for a subset of WT1 binding sites

Wilms tumor protein (WT1) is a Cys2-His2 zinc-finger transcription factor vital for embryonic development of the genitourinary system. The protein contains a C-terminal DNA binding domain with four tandem zinc-fingers (ZF1-4). An alternative splicing of Wt1 can add three additional amino acids-lysine (K), threonine (T) and serine (S)-between ZF3 and ZF4. In the -KTS isoform, ZF2-4 determine the sequence-specificity of DNA binding, whereas the function of ZF1 remains elusive. Three X-ray structures are described here for wild-type -KTS isoform ZF1-4 in complex with its cognate DNA sequence. We observed four unique ZF1 conformations. First, like ZF2-4, ZF1 can be positioned continuously in the DNA major groove forming a 'near-cognate' complex. Second, while ZF2-4 make base-specific interactions with one DNA molecule, ZF1 can interact with a second DNA molecule (or, presumably, two regions of the same DNA molecule). Third, ZF1 can intercalate at the joint of two tail-to-head DNA molecules. If such intercalation occurs on a continuous DNA molecule, it would kink the DNA at the ZF1 binding site. Fourth, two ZF1 units can dimerize. Furthermore, we examined a Denys-Drash syndrome-associated ZF1 mutation (methionine at position 342 is replaced by arginine). This mutation enhances WT1 affinity for a guanine base. X-ray crystallography of the mutant in complex with its preferred sequence revealed the interactions responsible for this affinity change. These results provide insight into the mechanisms of action of WT1, and clarify the fact that ZF1 plays a role in determining sequence specificity of this critical transcription factor.

Influence of Oligonucleotide Grafting Density on Surface-Mediated DNA Transport and Hybridization

Adsorption of soluble DNA to surfaces decorated with complementary DNA plays an important role in many bionanotechnology applications, and previous studies have reported complex dependencies of the surface density of immobilized DNA on hybridization. While these effects have been speculatively ascribed to steric or electrostatic effects, the influence of surface-mediated molecular transport (i.e., intermittent "hopping diffusion") has not been fully appreciated. Here, single-molecule tracking and Förster resonance energy transfer (FRET) were employed to characterize the mobility and the hybridization efficiency of adsorbed ssDNA oligonucleotides ("target") at solid-liquid interfaces exhibiting surface-immobilized ssDNA ("probe") over a wide range of surface grafting densities. Two distinct regimes were observed, with qualitatively different transport and hybridization behaviors. At dilute grafting density, only 1-3% of target molecules were observed to associate with probes (i.e., to hybridize). Adsorbing target molecules often searched unsuccessfully and "flew", via desorption-mediated diffusion, to secondary locations before hybridizing. In contrast, at high probe grafting density, approximately 20% of target DNA hybridized to immobilized probes, and almost always in the vicinity of initial adsorption. Moreover, following a dehybridization event, target molecules rehybridized at high probe density, but rehybridization was infrequent in the dilute density regime. Interestingly, the intermittent interfacial transport of mobile target molecules was suppressed by the presence of immobilized probe DNA, presumably due to an increased probability of readsorption following each "hop". Together, these findings suggested that many salient effects of grafting density on surface-mediated DNA hybridization can be directly related to the mechanisms of surface-mediated intermittent diffusion.

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.

DNA scanning by base excision repair enzymes and implications for pathway coordination

Site-specific DNA binding proteins must search the genome to locate their target sites, and many DNA modifying enzymes have the ability to scan along DNA in search of their substrates. This process is termed processive searching, and it serves to decrease the search time by effectively increasing the DNA binding footprint of a protein. The repertoire of proteins capable of processive searching is expanding, highlighting the need to understand the governing principles behind this fundamental process. Many of the enzymes in the base excision DNA repair pathway are capable of processive searching. Here, we briefly summarize methodology for determining if a protein can scan DNA and highlight the discovery that the base excision repair DNA polymerase β performs a processive search. Elucidation of physical models for DNA searching has also provided a plausible mechanism for pathway coordination during repair. The ability of BER enzymes to transiently sample adjacent DNA sites while bound to their product confers accessibility to downstream enzymes and does not require protein-protein interactions for coordination.

Facilitated Diffusion Mechanisms in DNA Base Excision Repair and Transcriptional Activation

Preservation of the coding potential of the genome and highly regulated gene expression over the life span of a human are two fundamental requirements of life. These processes require the action of repair enzymes or transcription factors that efficiently recognize specific sites of DNA damage or transcriptional regulation within a restricted time frame of the cell cycle or metabolism. A failure of these systems to act results in accumulated mutations, metabolic dysfunction, and disease. Despite the multifactorial complexity of cellular DNA repair and transcriptional regulation, both processes share a fundamental physical requirement that the proteins must rapidly diffuse to their specific DNA-binding sites that are embedded within the context of a vastly greater number of nonspecific DNA-binding sites. Superimposed on the needle-in-the-haystack problem is the complex nature of the cellular environment, which contains such high concentrations of macromolecules that the time frame for diffusion is expected to be severely extended as compared to dilute solution. Here we critically review the mechanisms for how these proteins solve the needle-in-the-haystack problem and how the effects of cellular macromolecular crowding can enhance facilitated diffusion processes. We restrict the review to human proteins that use stochastic, thermally driven site-recognition mechanisms, and we specifically exclude systems involving energy cofactors or circular DNA clamps. Our scope includes ensemble and single-molecule studies of the past decade or so, with an emphasis on connecting experimental observations to biological function.

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.

NMR-based investigations into target DNA search processes of proteins

To perform their function, transcription factors and DNA-repair/modifying enzymes must first locate their targets in the vast presence of nonspecific, but structurally similar sites on genomic DNA. Before reaching their targets, these proteins stochastically scan DNA and dynamically move from one site to another on DNA. Solution NMR spectroscopy provides unique atomic-level insights into the dynamic DNA-scanning processes, which are difficult to gain by any other experimental means. In this review, we provide an introductory overview on the NMR methods for the structural, dynamic, and kinetic investigations of target DNA search by proteins. We also discuss advantages and disadvantages of these NMR methods over other methods such as single-molecule techniques and biochemical approaches.

Acceleration of bursty multiprotein target search kinetics on DNA by colocalisation

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

Measurement of nanoscale DNA translocation by uracil DNA glycosylase in human cells

DNA ‘sliding’ by human repair enzymes is considered to be important for DNA damage detection. Here, we transfected uracil-containing DNA duplexes into human cells and measured the probability that nuclear human uracil DNA glycosylase (hUNG2) excised two uracil lesions spaced 10–80 bp apart in a single encounter without escaping the micro-volume containing the target sites. The two-site transfer probabilities were 100% and 54% at a 10 and 40 bp spacing, but dropped to only 10% at 80 bp. Enzyme trapping experiments suggested that site transfers over 40 bp followed a DNA ‘hopping’ pathway in human cells, indicating that authentic sliding does not occur even over this short distance. The transfer probabilities were much greater than observed in aqueous buffers, but similar to in vitro measurements in the presence of polymer crowding agents. The findings reveal a new role for the crowded nuclear environment in facilitating DNA damage detection.

Enzyme cycling contributes to efficient induction of genome mutagenesis by the cytidine deaminase APOBEC3B

The single-stranded DNA cytidine deaminases APOBEC3B, APOBEC3H haplotype I, and APOBEC3A can contribute to cancer through deamination of cytosine to form promutagenic uracil in genomic DNA. The enzymes must access single-stranded DNA during the dynamic processes of DNA replication or transcription, but the enzymatic mechanisms enabling this activity are not known. To study this, we developed a method to purify full length APOBEC3B and characterized it in comparison to APOBEC3A and APOBEC3H on substrates relevant to cancer mutagenesis. We found that the ability of an APOBEC3 to cycle between DNA substrates determined whether it was able to efficiently deaminate single-stranded DNA produced by replication and single-stranded DNA bound by replication protein A (RPA). APOBEC3 deaminase activity during transcription had a size limitation that inhibited APOBEC3B tetramers, but not APOBEC3A monomers or APOBEC3H dimers. Altogether, the data support a model in which the availability of single-stranded DNA is necessary, but alone not sufficient for APOBEC3-induced mutagenesis in cells because there is also a dependence on the inherent biochemical properties of the enzymes. The biochemical properties identified in this study can be used to measure the mutagenic potential of other APOBEC enzymes in the genome.

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.

Potential role of DNA methylation as a facilitator of target search processes for transcription factors through interplay with methyl-CpG-binding proteins

Eukaryotic genomes contain numerous non-functional high-affinity sequences for transcription factors. These sequences potentially serve as natural decoys that sequester transcription factors. We have previously shown that the presence of sequences similar to the target sequence could substantially impede association of the transcription factor Egr-1 with its targets. In this study, using a stopped-flow fluorescence method, we examined the kinetic impact of DNA methylation of decoys on the search process of the Egr-1 zinc-finger protein. We analyzed its association with an unmethylated target site on fluorescence-labeled DNA in the presence of competitor DNA duplexes, including Egr-1 decoys. DNA methylation of decoys alone did not affect target search kinetics. In the presence of the MeCP2 methyl-CpG-binding domain (MBD), however, DNA methylation of decoys substantially (∼10-30-fold) accelerated the target search process of the Egr-1 zinc-finger protein. This acceleration did not occur when the target was also methylated. These results suggest that when decoys are methylated, MBD proteins can block them and thereby allow Egr-1 to avoid sequestration in non-functional locations. This effect may occur in vivo for DNA methylation outside CpG islands (CGIs) and could facilitate localization of some transcription factors within regulatory CGIs, where DNA methylation is rare.

Discrete-State Kinetics Model for NMR-Based Analysis of Protein Translocation on DNA at Equilibrium

In the target DNA search process, sequence-specific DNA-binding proteins first nonspecifically bind to DNA and stochastically move from one site to another before reaching their targets. To rigorously assess how the translocation process influences NMR signals from proteins interacting with nonspecific DNA, we incorporated a discrete-state kinetic model for protein translocation on DNA into the McConnell equation. Using this equation, we simulated line shapes of NMR signals from proteins undergoing translocations on DNA through sliding, dissociation/reassociation, and intersegment transfer. Through this analysis, we validated an existing NMR approach for kinetic investigations of protein translocation on DNA, which utilizes NMR line shapes of two nonspecific DNA-protein complexes and their mixture. We found that, despite its use of simplistic two-state approximation neglecting the presence of many microscopic states, the previously proposed NMR approach provides accurate kinetic information on the intermolecular translocations of proteins between two DNA molecules. Interestingly, our results suggest that the same NMR approach can also provide qualitative information about the one-dimensional diffusion coefficient for proteins sliding on DNA.

Processive searching ability varies among members of the gap-filling DNA polymerase X family

DNA repair proteins must locate rare damaged sites within the genome. DNA polymerase β (Pol β), a member of the DNA polymerase X family that is involved in base excision repair, uses a processive hopping search mechanism to locate substrates. This effectively enhances its search footprint on DNA, increasing the probability of locating damaged sites. Processive searching has been reported or proposed for many DNA-binding proteins, raising the question of how widespread or specific to certain enzymes the ability to perform this function is. To provide insight into this question, we compared the ability of three homologous DNA Pol X family members to perform a processive search for 1-nucleotide gaps in DNA using a previously developed biochemical assay. We found that at near-predicted physiological ionic strengths, the intramolecular searching ability of Pol β is at least 4-fold higher than that of Pol μ and ∼2-fold higher than that of Pol λ. Pol β also was able to perform intersegmental transfer with the intersegmental searching ability of Pol β being at least 6- and ∼2-fold higher than that of Pols μ and λ, respectively. Mutational analysis suggested that differences in the N-terminal domains of these polymerases are responsible for the varying degrees of searching competence. Of note, the differences in processive searching ability observed among the DNA Pol X family members correlated with their proposed biological functions in base excision repair and nonhomologous end joining.

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.

Thermodynamic Additivity for Impacts of Base-Pair Substitutions on Association of the Egr-1 Zinc-Finger Protein with DNA

The transcription factor Egr-1 specifically binds as a monomer to its 9 bp target DNA sequence, GCGTGGGCG, via three zinc fingers and plays important roles in the brain and cardiovascular systems. Using fluorescence-based competitive binding assays, we systematically analyzed the impacts of all possible single-nucleotide substitutions in the target DNA sequence and determined the change in binding free energy for each. Then, we measured the changes in binding free energy for sequences with multiple substitutions and compared them with the sum of the changes in binding free energy for each constituent single substitution. For the DNA variants with two or three nucleotide substitutions in the target sequence, we found excellent agreement between the measured and predicted changes in binding free energy. Interestingly, however, we found that this thermodynamic additivity broke down with a larger number of substitutions. For DNA sequences with four or more substitutions, the measured changes in binding free energy were significantly larger than predicted. On the basis of these results, we analyzed the occurrences of high-affinity sequences in the genome and found that the genome contains millions of such sequences that might functionally sequester Egr-1.

Sequence heterogeneity accelerates protein search for targets on DNA

The process of protein search for specific binding sites on DNA is fundamentally important since it marks the beginning of all major biological processes. We present a theoretical investigation that probes the role of DNA sequence symmetry, heterogeneity, and chemical composition in the protein search dynamics. Using a discrete-state stochastic approach with a first-passage events analysis, which takes into account the most relevant physical-chemical processes, a full analytical description of the search dynamics is obtained. It is found that, contrary to existing views, the protein search is generally faster on DNA with more heterogeneous sequences. In addition, the search dynamics might be affected by the chemical composition near the target site. The physical origins of these phenomena are discussed. Our results suggest that biological processes might be effectively regulated by modifying chemical composition, symmetry, and heterogeneity of a genome.

Regulation of transcription factors via natural decoys in genomic DNA

Eukaryotic genomic DNA contains numerous high-affinity sites for transcription factors. Only a small fraction of these sites directly regulates target genes. Other high-affinity sites can serve as naturally present decoys that sequester transcription factors. Such natural decoys in genomic DNA may provide novel regulatory mechanisms for transcription factors.

Single-molecule fluorescence studies on DNA looping

Structure and dynamics of DNA impact how the genetic code is processed and maintained. In addition to its biological importance, DNA has been utilized as building blocks of various nanomachines and nanostructures. Thus, understanding the physical properties of DNA is of fundamental importance to basic sciences and engineering applications. DNA can undergo various physical changes. Among them, DNA looping is unique in that it can bring two distal sites together, and thus can be used to mediate interactions over long distances. In this paper, we introduce a FRET-based experimental tool to study DNA looping at the single molecule level. We explain the connection between experimental measurables and a theoretical concept known as the J factor with the intent of raising awareness of subtle theoretical details that should be considered when drawing conclusions. We also explore DNA looping-assisted protein diffusion mechanism called intersegmental transfer using protein induced fluorescence enhancement (PIFE). We present some preliminary results and future outlooks.

Crowding on DNA in Protein Search for Targets

Proteins searching and recognizing specific sites on DNA is required for initiating all major biological processes. While the details of the protein search for targets on DNA in purified in vitro systems are reasonably well understood, the situation in real cells is much less clear. The presence of other types of molecules on DNA should prevent reaching the targets, but experiments show that, surprisingly, the molecular crowding on DNA influences the search dynamics much less than expected. We develop a theoretical method that allowed us to clarify the mechanisms of the protein search on DNA in the presence of crowding. It is found that the dimensionality of the search trajectories specifies whether the crowding will affect the target finding. For 3D search pathways it is minimal, while the strongest effect is for 1D search pathways when the crowding particle can block the search. In addition, for 1D search we determined that the critical parameter is a mobility of crowding agents: highly mobile molecules do not affect the search dynamics, while the slow particles can significantly slow down the process. Physical-chemical explanations of the observed phenomena are presented. Our theoretical predictions thus explain the experimental observations, and they are also supported by extensive numerical simulations.

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.

Changes in conformational dynamics of basic side chains upon protein-DNA association

Basic side chains play major roles in recognition of nucleic acids by proteins. However, dynamic properties of these positively charged side chains are not well understood. In this work, we studied changes in conformational dynamics of basic side chains upon protein–DNA association for the zinc-finger protein Egr-1. By nuclear magnetic resonance (NMR) spectroscopy, we characterized the dynamics of all side-chain cationic groups in the free protein and in the complex with target DNA. Our NMR order parameters indicate that the arginine guanidino groups interacting with DNA bases are strongly immobilized, forming rigid interfaces. Despite the strong short-range electrostatic interactions, the majority of the basic side chains interacting with the DNA phosphates exhibited high mobility, forming dynamic interfaces. In particular, the lysine side-chain amino groups exhibited only small changes in the order parameters upon DNA-binding. We found a similar trend in the molecular dynamics (MD) simulations for the free Egr-1 and the Egr-1–DNA complex. Using the MD trajectories, we also analyzed side-chain conformational entropy. The interfacial arginine side chains exhibited substantial entropic loss upon binding to DNA, whereas the interfacial lysine side chains showed relatively small changes in conformational entropy. These data illustrate different dynamic characteristics of the interfacial arginine and lysine side chains.

Residence Times of Molecular Complexes in Solution from NMR Data of Intermolecular Hydrogen-Bond Scalar Coupling

The residence times of molecular complexes in solution are important for understanding biomolecular functions and drug actions. We show that NMR data of intermolecular hydrogen-bond scalar couplings can yield information on the residence times of molecular complexes in solution. The molecular exchange of binding partners via the breakage and reformation of a complex causes self-decoupling of intermolecular hydrogen-bond scalar couplings, and this self-decoupling effect depends on the residence time of the complex. For protein-DNA complexes, we investigated the salt concentration dependence of intermolecular hydrogen-bond scalar couplings between the protein side-chain (15)N and DNA phosphate (31)P nuclei, from which the residence times were analyzed. The results were consistent with those obtained by (15)Nz-exchange spectroscopy. This self-decoupling-based kinetic analysis is unique in that it does not require any different signatures for the states involved in the exchange, whereas such conditions are crucial for kinetic analyses by typical NMR and other methods.

Influence of quasi-specific sites on kinetics of target DNA search by a sequence-specific DNA-binding protein

Functions of transcription factors require formation of specific complexes at particular sites in cis-regulatory elements of genes. However, chromosomal DNA contains numerous sites that are similar to the target sequences recognized by transcription factors. The influence of such “quasi-specific” sites on functions of the transcription factors is not well understood at present by experimental means. In this work, using fluorescence methods, we have investigated the influence of quasi-specific DNA sites on the efficiency of target location by the zinc finger DNA-binding domain of the inducible transcription factor Egr-1, which recognizes a 9 bp sequence. By stopped-flow assays, we measured the kinetics of Egr-1’s association with a target site on 143 bp DNA in the presence of various competitor DNAs, including nonspecific and quasi-specific sites. The presence of quasi-specific sites on competitor DNA significantly decelerated the target association by the Egr-1 protein. The impact of the quasi-specific sites depended strongly on their affinity, their concentration, and the degree of their binding to the protein. To quantitatively describe the kinetic impact of the quasi-specific sites, we derived an analytical form of the apparent kinetic rate constant for the target association and used it for fitting to the experimental data. Our kinetic data with calf thymus DNA as a competitor suggested that there are millions of high-affinity quasi-specific sites for Egr-1 among the 3 billion bp of genomic DNA. This study quantitatively demonstrates that naturally abundant quasi-specific sites on DNA can considerably impede the target search processes of sequence-specific DNA-binding proteins.

Physicochemical Properties of Ion Pairs of Biological Macromolecules

Ion pairs (also known as salt bridges) of electrostatically interacting cationic and anionic moieties are important for proteins and nucleic acids to perform their function. Although numerous three-dimensional structures show ion pairs at functionally important sites of biological macromolecules and their complexes, the physicochemical properties of the ion pairs are not well understood. Crystal structures typically show a single state for each ion pair. However, recent studies have revealed the dynamic nature of the ion pairs of the biological macromolecules. Biomolecular ion pairs undergo dynamic transitions between distinct states in which the charged moieties are either in direct contact or separated by water. This dynamic behavior is reasonable in light of the fundamental concepts that were established for small ions over the last century. In this review, we introduce the physicochemical concepts relevant to the ion pairs and provide an overview of the recent advancement in biophysical research on the ion pairs of biological macromolecules.

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

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