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

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

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

Skipping events impose repeated binding attempts: profound kinetic implications of protein-DNA conformational changes

The kinetics of protein-DNA recognition, along with its thermodynamic properties, including affinity and specificity, play a central role in shaping biological function. Protein-DNA recognition kinetics are characterized by two key elements: the time taken to locate the target site amid various nonspecific alternatives; and the kinetics involved in the recognition process, which may necessitate overcoming an energetic barrier. In this study, we developed a coarse-grained (CG) model to investigate interactions between a transcription factor called the sex-determining region Y (SRY) protein and DNA, in order to probe how DNA conformational changes affect SRY-DNA recognition and binding kinetics. We find that, not only does a requirement for such a conformational DNA transition correspond to a higher energetic barrier for binding and therefore slower kinetics, it may further impede the recognition kinetics by increasing unsuccessful binding events (skipping events) where the protein partially binds its DNA target site but fails to form the specific protein-DNA complex. Such skipping events impose the need for additional cycles protein search of nonspecific DNA sites, thus significantly extending the overall recognition time. Our results highlight a trade-off between the speed with which the protein scans nonspecific DNA and the rate at which the protein recognizes its specific target site. Finally, we examine molecular approaches potentially adopted by natural systems to enhance protein-DNA recognition despite its intrinsically slow kinetics.

Nonspecific vs. specific DNA binding free energetics of a transcription factor domain protein

Transcription factor (TF) proteins regulate gene expression by binding to specific sites on the genome. In the facilitated diffusion model, an optimized search process is achieved by the TF alternating between 3D diffusion in the bulk and 1D diffusion along DNA. While undergoing 1D diffusion, the protein can switch from a search mode for fast diffusion along nonspecific DNA to a recognition mode for stable binding to specific DNA. It was recently noticed that, for a small TF domain protein, reorientations on DNA happen between the nonspecific and specific DNA binding. We here conducted all-atom molecular dynamics simulations with steering forces to reveal the protein-DNA binding free energetics, confirming that the search and recognition modes are distinguished primarily by protein orientations on the DNA. As the binding free energy difference between the specific and nonspecific DNA system slightly deviates from that being estimated directly from dissociation constants on 15-bp DNA constructs, we hypothesize that the discrepancy can come from DNA sequences flanking the 6-bp central binding sites that impact on the dissociation kinetics measurements. The hypothesis is supported by a simplified spherical protein-DNA model along with stochastic simulations and kinetic modeling.

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.

Can molecular dynamics be used to simulate biomolecular recognition?

There are many problems in biochemistry that are difficult to study experimentally. Simulation methods are appealing due to direct availability of atomic coordinates as a function of time. However, direct molecular simulations are challenged by the size of systems and the time scales needed to describe relevant motions. In theory, enhanced sampling algorithms can help to overcome some of the limitations of molecular simulations. Here, we discuss a problem in biochemistry that offers a significant challenge for enhanced sampling methods and that could, therefore, serve as a benchmark for comparing approaches that use machine learning to find suitable collective variables. In particular, we study the transitions LacI undergoes upon moving between being non-specifically and specifically bound to DNA. Many degrees of freedom change during this transition and that the transition does not occur reversibly in simulations if only a subset of these degrees of freedom are biased. We also explain why this problem is so important to biologists and the transformative impact that a simulation of it would have on the understanding of DNA regulation.

Importin alpha family NAAT/IBB domain: Functions of a pleiotropic long chameleon sequence

Nuclear transport is essential for eukaryotic cell survival and regulates the movement of functional molecules in and out of the nucleus via the nuclear pore. Transport is facilitated by protein-protein interactions between cargo and transport receptors, which contribute to the expression and regulation of downstream genetic information. This chapter focuses on the molecular basis of the multifunctional nature of the importin α family, the representative transport receptors that bring proteins into the nucleus. Importin α performs multiple functions during the nuclear transport cycle through interactions with multiple molecules by a single domain called the IBB domain. This domain is a long chameleon sequence, which can change its conformation and binding mode depending on the interaction partners. By considering the evolutionarily conserved biochemical/physicochemical propensities of the amino acids constituting the functional complex interfaces, together with their structural properties, the mechanisms of switching between multiple complexes formed via IBB and the regulation of downstream functions are examined in detail. The mechanism of regulation by IBB indicates that the time has come for a paradigm shift in the way we view the molecular mechanisms by which proteins regulate downstream functions through their interactions with other molecules.

Conformational Change of Transcription Factors from Search to Specific Binding: A lac Repressor Case Study

In a process known as facilitated diffusion, DNA-binding proteins find their target sites by combining three-dimensional diffusion and one-dimensional scanning of the DNA. Following the trade-off between speed and stability, agile exploration of DNA requires loose binding, whereas, at the DNA target site, the searching protein needs to establish tight interactions with the DNA. To enable both efficient search and stable binding, DNA-binding proteins and DNA often switch conformations upon recognition. Here, we study the one-dimensional diffusion and DNA binding of the dimeric lac repressor (LacI), which was reported to adopt two different conformations when binding different conformations of DNA. Using coarse-grained molecular dynamic simulations, we studied the diffusion and the sequence-specific binding of these conformations of LacI, as well as their truncated or monomeric variants, with two DNA conformations: straight and bent. The simulations were compared to experimental observables. This study supports that linear diffusion along DNA combines tight rotation-coupled groove tracking and rotation-decoupled hopping, where the protein briefly dissociates and reassociates just a few base pairs away. Tight groove tracking is crucial for target-site recognition, while hopping speeds up the overall search process. We investigated the diffusion of different LacI conformations on DNA and show how the flexibility of LacI's hinge regions ensures agility on DNA as well as faithful groove tracking. If the hinge regions instead form α-helices at the protein-DNA interface, tight groove tracking is not possible. On the contrary, the helical hinge region is essential for tight binding to bent, specific DNA, for the formation of the specific complex. Based on our study of different encounter complexes, we argue that the conformational change in LacI and DNA bending are somewhat coupled. Our findings underline the importance of two distinct protein conformations for facilitated diffusion and specific binding, respectively.

Function and regulation of transcription factors during mitosis-to-G1 transition

During cell division, drastic cellular changes characteristic of mitosis result in the inactivation of the transcriptional machinery, and global downregulation of transcription. Sequence-specific transcription factors (TFs) have thus been considered mere bystanders, devoid of any regulatory function during mitosis. This view changed significantly in recent years, upon the conclusion that many TFs associate with condensed chromosomes during cell division, even occupying a fraction of their genomic target sites in mitotic chromatin. This finding was at the origin of the concept of mitotic bookmarking by TFs, proposed as a mechanism to propagate gene regulatory information across cell divisions, by facilitating the reactivation of specific bookmarked genes. While the underlying mechanisms and biological significance of this model remain elusive, recent developments in this fast-moving field have cast new light into TF activity during mitosis, beyond a bookmarking role. Here, we start by reviewing the most recent findings on the complex nature of TF-chromatin interactions during mitosis, and on mechanisms that may regulate them. Next, and in light of recent reports describing how transcription is reinitiated in temporally distinct waves during mitosis-to-G1 transition, we explore how TFs may contribute to defining this hierarchical gene expression process. Finally, we discuss how TF activity during mitotic exit may impact the acquisition of cell identity upon cell division, and propose a model that integrates dynamic changes in TF-chromatin interactions during this cell-cycle period, with the execution of cell-fate decisions.

Investigating the Conformational Dynamics of a Y-Family DNA Polymerase during Its Folding and Binding to DNA and a Nucleotide

During DNA polymerization, the Y-family DNA polymerases are capable of bypassing various DNA damage, which can stall the replication fork progression. It has been well acknowledged that the structures of the Y-family DNA polymerases have been naturally evolved to undertake this vital task. However, the mechanisms of how these proteins utilize their unique structural and conformational dynamical features to perform the translesion DNA synthesis are less understood. Here, we developed structure-based models to study the precatalytic DNA polymerization process, including DNA and nucleotide binding to DPO4, a paradigmatic Y-family polymerase from Sulfolobus solfataricus. We studied the interplay between the folding and the conformational dynamics of DPO4 and found that DPO4 undergoes first unraveling (unfolding) and then folding for accomplishing the functional "open-to-closed" conformational transition. DNA binding dynamically modulates the conformational equilibrium in DPO4 during the stepwise binding through different types of interactions, leading to different conformational distributions of DPO4 at different DNA binding stages. We observed that nucleotide binding induces modulation of a few contacts surrounding the active site of the DPO4-DNA complex associated with a high free energy barrier. Our simulation results resonate with the experimental evidence that the conformational change at the active site led by nucleotide is the rate-limiting step of nucleotide incorporation. In combination with localized frustration analyses, we underlined the importance of DPO4 conformational dynamics and fluctuations in facilitating DNA and nucleotide binding. Our findings offer mechanistic insights into the processes of DPO4 conformational dynamics associated with the substrate binding and contribute to the understanding of the "structure-dynamics-function" relationship in the Y-family DNA polymerases.

Overview of the Biomolecular Association and Dynamics session at the 20th IUPAB congress, 45th Brazilian congress of SBBF, and the 50th annual meeting of SBBq

In this session, experts in molecular biophysics described the dynamics of biopolymers across a wide range of length and time scales. This discussion highlighted numerous techniques that span from highly detailed simulations, to coarse-grained theoretical models, as well as high-resolution structural analysis. The topics were equally diverse, where there was discussion of biological processes at small (individual atoms), intermediate (assemblies) and very large scales (phase separation).

Dynamics of Proteins and Macromolecular Machines in Escherichia coli

Proteins are major contributors to the composition and the functions in the cell. They often assemble into larger structures, macromolecular machines, to carry out intricate essential functions. Although huge progress in understanding how macromolecular machines function has been made by reconstituting them in vitro, the role of the intracellular environment is still emerging. The development of fluorescence microscopy techniques in the last 2 decades has allowed us to obtain an increased understanding of proteins and macromolecular machines in cells. Here, we describe how proteins move by diffusion, how they search for their targets, and how they are affected by the intracellular environment. We also describe how proteins assemble into macromolecular machines and provide examples of how frequent subunit turnover is used for them to function and to respond to changes in the intracellular conditions. This review emphasizes the constant movement of molecules in cells, the stochastic nature of reactions, and the dynamic nature of macromolecular machines.

Direct visualization of the effect of DNA structure and ionic conditions on HU-DNA interactions

Architectural DNA–binding proteins are involved in many important DNA transactions by virtue of their ability to change DNA conformation. Histone-like protein from E. coli strain U93, HU, is one of the most studied bacterial architectural DNA–binding proteins. Nevertheless, there is still a limited understanding of how the interactions between HU and DNA are affected by ionic conditions and the structure of DNA. Here, using optical tweezers in combination with fluorescent confocal imaging, we investigated how ionic conditions affect the interaction between HU and DNA. We directly visualized the binding and the diffusion of fluorescently labelled HU dimers on DNA. HU binds with high affinity and exhibits low mobility on the DNA in the absence of Mg²⁺; it moves 30-times faster and stays shorter on the DNA with 8 mM Mg²⁺ in solution. Additionally, we investigated the effect of DNA tension on HU–DNA complexes. On the one hand, our studies show that binding of HU enhances DNA helix stability. On the other hand, we note that the binding affinity of HU for DNA in the presence of Mg²⁺ increases at tensions above 50 pN, which we attribute to force-induced structural changes in the DNA. The observation that HU diffuses faster along DNA in presence of Mg²⁺ compared to without Mg²⁺ suggests that the free energy barrier for rotational diffusion along DNA is reduced, which can be interpreted in terms of reduced electrostatic interaction between HU and DNA, possibly coinciding with reduced DNA bending.

Importin α2 association with chromatin: Direct DNA binding via a novel DNA-binding domain

The nuclear transport of proteins is important for facilitating appropriate nuclear functions. The importin α family proteins play key roles in nuclear transport as transport receptors for copious nuclear proteins. Additionally, these proteins possess other functions, including chromatin association and gene regulation. However, these nontransport functions of importin α are not yet fully understood, especially their molecular-level mechanisms and consequences for functioning with chromatin. Here, we report the novel molecular characteristics of importin α binding to diverse DNA sequences in chromatin. We newly identified and characterized a DNA-binding domain-the Nucleic Acid Associating Trolley pole domain (NAAT domain)-in the N-terminal region of importin α within the conventional importin β binding (IBB) domain that is necessary for nuclear transport of cargo proteins. Furthermore, we found that the DNA binding of importin α synergistically coupled the recruitment of its cargo protein to DNA. This is the first study to delineate the interaction between importin α and chromatin DNA via the NAAT domain, indicating the bifunctionality of the importin α N-terminal region for nuclear transport and chromatin association.

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.

Regulator-dependent temporal dynamics of a restriction-modification system's gene expression upon entering new host cells: single-cell and population studies

Restriction-modification (R-M) systems represent a first line of defense against invasive DNAs, such as bacteriophage DNAs, and are widespread among bacteria and archaea. By acquiring a Type II R-M system via horizontal gene transfer, the new hosts generally become more resistant to phage infection, through the action of a restriction endonuclease (REase), which cleaves DNA at or near specific sequences. A modification methyltransferase (MTase) serves to protect the host genome against its cognate REase activity. The production of R-M system components upon entering a new host cell must be finely tuned to confer protective methylation before the REase acts, to avoid host genome damage. Some type II R-M systems rely on a third component, the controller (C) protein, which is a transcription factor that regulates the production of REase and/or MTase. Previous studies have suggested C protein effects on the dynamics of expression of an R-M system during its establishment in a new host cell. Here, we directly examine these effects. By fluorescently labelling REase and MTase, we demonstrate that lack of a C protein reduces the delay of REase production, to the point of being simultaneous with, or even preceding, production of the MTase. Single molecule tracking suggests that a REase and a MTase employ different strategies for their target search within host cells, with the MTase spending much more time diffusing in proximity to the nucleoid than does the REase. This difference may partially ameliorate the toxic effects of premature REase expression.

What Are the Molecular Requirements for Protein Sliding along DNA?

DNA-binding proteins rely on linear diffusion along the longitudinal DNA axis, supported by their nonspecific electrostatic affinity for DNA, to search for their target recognition sites. One may therefore expect that the ability to engage in linear diffusion along DNA is universal to all DNA-binding proteins, with the detailed biophysical characteristics of that diffusion differing between proteins depending on their structures and functions. One key question is whether the linear diffusion mechanism is defined by translation coupled with rotation, a mechanism that is often termed sliding. We conduct coarse-grained and atomistic molecular dynamics simulations to investigate the minimal requirements for protein sliding along DNA. We show that coupling, while widespread, is not universal. DNA-binding proteins that slide along DNA transition to uncoupled translation–rotation (i.e., hopping) at higher salt concentrations. Furthermore, and consistently with experimental reports, we find that the sliding mechanism is the less dominant mechanism for some DNA-binding proteins, even at low salt concentrations. In particular, the toroidal PCNA protein is shown to follow the hopping rather than the sliding mechanism.

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.

Protein-DNA target search relies on quantum walk

Protein-DNA interactions play a fundamental role in all life systems. A critical issue of such interactions is given by the strategy of protein search for specific targets on DNA. The mechanisms by which the protein are able to find relatively small cognate sequences, typically 15-20 base pairs (bps) for repressors, and 4-6 bps for restriction enzymes among the millions of bp of non-specific chromosomal DNA have hardly engaged researchers for decades. Recent experimental studies have generated new insights on the basic processes of protein-DNA interactions evidencing the underlying complex dynamic phenomena involved, which combine three-dimensional and one-dimensional motion along the DNA chain. It has been demonstrated that protein molecules have an extraordinary ability to find the target very quickly on the DNA chain, in some cases, with two orders of magnitude faster than the diffusion limit. This unique property of protein-DNA search mechanism is known as facilitated diffusion. Several theoretical mechanisms have been suggested to describe the origin of facilitated diffusion. However, none of such models currently has the ability to fully describe the protein search strategy. In this paper, we suggest that the ability of proteins to identify consensus sequences on DNA is based on the entanglement of π-π electrons between DNA nucleotides and protein amino acids. The π-π entanglement is based on Quantum Walk (QW), through Coin-position entanglement (CPE). First, the protein identifies a dimer belonging to the consensus sequence, and localize a π on such dimer, hence, the other π electron scans the DNA chain until the sequence is identified. Focusing on the example of recognition of consensus sequences of EcoRV or EcoRI, we will describe the quantum features of QW on protein-DNA complexes during the search strategy, such as walker quadratic spreading on a coherent superposition of different vertices and environment-supported long-time survival probability of the walker. We will employ both discrete- or continuous-time versions of QW. Biased and unbiased classical Random Walk (CRW) have been used for a long time to describe the Protein-DNA search strategy. QW, the quantum version of CRW, has been widely studied for its applications in quantum information applications. In our biological application, the walker (the protein) resides at a vertex in a graph (the DNA structural topology). Differently to CRW, where the walker moves randomly, the quantum walker can hop along the edges in the graph to reach other vertices entering coherently a superposition across different vertices spreading quadratically faster than CRW analogous evidencing the typical speed up features of the QW. When applied to a protein-DNA target search problem, QW gives the possibility to achieve the experimental diffusional motion of proteins over diffusion classical limits experienced along DNA chains exploiting quantum features such as CPE and long-time survival probability supported by the environment. In turn, we come to the conclusion that, under quantum picture, the protein search strategy does not distinguish between one-dimensional (1D) and three-dimensional (3D) cases.

Prerecognition Diffusion Mechanism of Human DNA Mismatch Repair Proteins along DNA: Msh2-Msh3 versus Msh2-Msh6

DNA mismatch repair (MMR) is an important postreplication process that eliminates mispaired or unpaired nucleotides to ensure genomic replication fidelity. In humans, Msh2-Msh6 and Msh2-Msh3 are the two mismatch repair initiation factors that recognize DNA lesions. While X-ray crystal structures exist for these proteins in complex with DNA lesions, little is known about their structures during the initial search along nonspecific double-stranded DNA, because they are short-lived and difficult to determine experimentally. In this study, various computational approaches were used to sidestep these difficulties. All-atom and coarse-grained simulations based on the crystal structures of Msh2-Msh3 and Msh2-Msh6 showed no translation along the DNA, suggesting that the initial search conformation differs from the lesion-bound crystal structure. We modeled probable search-mode structures of MSH proteins and showed, using coarse-grained molecular dynamics simulations, that they can perform rotation-coupled diffusion on DNA, which is a suitable and efficient search mechanism for their function and one predicted earlier by fluorescence resonance energy transfer and fluorescence microscopy studies. This search mechanism is implemented by electrostatic interactions among the mismatch-binding domain (MBD), the clamp domains, and the DNA backbone. During simulations, their diffusion rate did not change significantly with an increasing salt concentration, which is consistent with observations from experimental studies. When the gap between their DNA-binding clamps was increased, Msh2-Msh3 diffused mostly via the clamp domains while Msh2-Msh6 still diffused using the MBD, reproducing the experimentally measured lower diffusion coefficient of Msh2-Msh6. Interestingly, Msh2-Msh3 was capable of dissociating from the DNA, whereas Msh2-Msh6 always diffused on the DNA duplex. This is consistent with the experimental observation that Msh2-Msh3, unlike Msh2-Msh6, can overcome obstacles such as nucleosomes. Our models provide a molecular picture of the different mismatch search mechanisms undertaken by Msh2-Msh6 and Msh2-Msh3, despite the similarity of their structures.

A lattice model on the rate of in vivo site-specific DNA-protein interactions

We develop a lattice model of site-specific DNA-protein interactions under in vivo conditions where DNA is modelled as a self-avoiding random walk that is embedded in a cubic lattice box resembling the living cell. The protein molecule searches for its cognate site on DNA via a combination of three dimensional (3D) and one dimensional (1D) random walks. Hopping and intersegmental transfers occur depending on the conformational state of DNA. Results show that the search acceleration ratio (= search time in pure 3D route/search time in 3D and 1D routes) asymptotically increases towards a limiting value as the dilution factor of DNA (= volume of the cell/the volume of DNA) tends towards infinity. When the dilution ratio is low, then hopping and intersegmental transfers significantly enhance the search efficiency over pure sliding. At high dilution ratio, hopping does not enhance the search efficiency much since under such situation DNA will be in a relaxed conformation that favors only sliding. In the absence of hopping and intersegmental transfers, there exists an optimum sliding time at which the search acceleration ratio attains a maximum in line with the current theoretical results. However, existence of such optimum sliding length disappears in the presence of hopping. When the DNA is confined in a small volume inside the cell resembling a natural cell system, then there exists an optimum dilution and compression ratios (= total cell volume/volume in which DNA is confined) at which the search acceleration factor attains a maximum especially in the presence of hopping and intersegmental transfers. These optimum values are consistent with the values observed in the Escherichia coli cell system. In the absence of confinement of DNA, position of the specific binding site on the genomic DNA significantly influences the search acceleration. However, such position dependent changes in the search acceleration ratio will be nullified in the presence of hopping and intersegmental transfers especially when the DNA is confined in a small volume that is embedded in an outer cell.

Lost in the Crowd: How Does Human 8-Oxoguanine DNA Glycosylase 1 (OGG1) Find 8-Oxoguanine in the Genome?

The most frequent DNA lesion resulting from an oxidative stress is 7,8-dihydro-8-oxoguanine (8-oxoG). 8-oxoG is a premutagenic base modification due to its capacity to pair with adenine. Thus, the repair of 8-oxoG is critical for the preservation of the genetic information. Nowadays, 8-oxoG is also considered as an oxidative stress-sensor with a putative role in transcription regulation. In mammalian cells, the modified base is excised by the 8-oxoguanine DNA glycosylase (OGG1), initiating the base excision repair (BER) pathway. OGG1 confronts the massive challenge that is finding rare occurrences of 8-oxoG among a million-fold excess of normal guanines. Here, we review the current knowledge on the search and discrimination mechanisms employed by OGG1 to find its substrate in the genome. While there is considerable data from in vitro experiments, much less is known on how OGG1 is recruited to chromatin and scans the genome within the cellular nucleus. Based on what is known of the strategies used by proteins searching for rare genomic targets, we discuss the possible scenarios allowing the efficient detection of 8-oxoG by OGG1.

Value of Temporal Information When Analyzing Reaction Coordinates

Reaction coordinates chart pathways from reactants to products of chemical reactions. Determination of reaction coordinates from ensembles of molecular trajectories has thus been the focus of many studies. A widely used and insightful choice of a reaction coordinate is the committor function, defined as the probability that a trajectory will reach the product before the reactant. Here, we consider alternatives to the committor function that add useful mechanistic information, the mean first passage time, and the exit time to the product. We further derive a simple relationship between the functions of the committor, the mean first passage time, and the exit time. We illustrate the diversity of mechanisms predicted by alternative reaction coordinates with several toy problems and with a simple model of protein searching for a specific DNA motif.

C1188D mutation abolishes specific recognition between MLL1-CXXC domain and CpG site by inducing conformational switch of flexible N-terminal

Mixed lineage leukemia protein (MLL1 protein) recognizes the CpG site via its CXXC domain and is frequently associated with leukemia. The specific recognition is abolished by C1188D mutation, which also prevents MLL-related leukemia. In this paper, multiple molecular dynamic (MD) simulations were performed to investigate the mechanism of recognition and influences of C1188D mutation. Started from fully dissociated DNA and MLL1-CXXC domain, remarkably, the center of mass (COM) of MLL1-CXXC domain quickly concentrates on the vicinity of the CpG site in all 53 short MD simulations. Extended simulations of the wild type showed that the native complex formed in 500 ns among 4 of 53 simulations. In contrast, the C1188D mutant COM distributed broadly around the DNA and the native complex was not observed in any of the extended simulations. Simulations on the apo MLL1-CXXC domain further suggest that the wild type protein remained predominantly in an open form that closely resembles its structure in the native complex whereas C1188D mutant formed predominantly compact structures in which the N- terminal bends to D1188. This conformational switch hinders the formation of encounter complex, thus abolishes the recognition. Our study also provides clues to the study mechanism of recognition, by the CXXC domain from proteins like DNA methyltransferase and ten-eleven translocation enzymes.

Transcription Factors and DNA Play Hide and Seek

Transcription factors (TFs) bind to specific DNA motifs to regulate the expression of target genes. To reach their binding sites, TFs diffuse in 3D and perform local motions such as 1D sliding, hopping, or intersegmental transfer. TF-DNA interactions depend on multiple parameters, such as the chromatin environment, TF partitioning into distinct subcellular regions, and cooperativity with other DNA-binding proteins. In this review, how current understanding of the search process has initially been shaped by prokaryotic studies is discussed, as well as what is known about the parameters regulating TF search efficiency in the context of the complex eukaryotic chromatin landscape.