Diffusion of ring-shaped proteins along DNA: case study of sliding clamps

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.

The bacterial DNA sliding clamp, β-clamp: structure, interactions, dynamics and drug discovery

DNA replication is a tightly coordinated event carried out by a multiprotein replication complex. An essential factor in the bacterial replication complex is the ring-shaped DNA sliding clamp, β-clamp, ensuring processive DNA replication and DNA repair through tethering of polymerases and DNA repair proteins to DNA. β -clamp is a hub protein with multiple interaction partners all binding through a conserved clamp binding sequence motif. Due to its central role as a DNA scaffold protein, β-clamp is an interesting target for antimicrobial drugs, yet little effort has been put into understanding the functional interactions of β-clamp. In this review, we scrutinize the β-clamp structure and dynamics, examine how its interactions with a plethora of binding partners are regulated through short linear binding motifs and discuss how contexts play into selection. We describe the dynamic process of clamp loading onto DNA and cover the recent advances in drug development targeting β-clamp. Despite decades of research in β-clamps and recent landmark structural insight, much remains undisclosed fostering an increased focus on this very central protein.

Three-State Diffusion Model of DNA Glycosylase Translocation along Stretched DNA as Revealed by Free Energy Landscapes at the All-Atom Level

DNA glycosylases play key roles in the maintenance of genomic integrity. These enzymes effectively find rare damaged sites in DNA and participate in subsequent base excision repair. Single-molecule and ensemble experiments have revealed key aspects of this damage-site searching mechanism and the involvement of facilitated diffusion. In this study, we describe free energy landscapes of enzyme translocation along nonspecific DNA obtained using a fully atomistic molecular dynamics (MD) simulation of a well-known DNA glycosylase, human 8-oxoguanine DNA glycosylase 1 (hOGG1). Based on an analysis of simulated free energy profiles, we propose a three-state model for the damage-site searching mechanism. In the three states, named the L1, L2, and L3 states, the L1 state is a helical sliding mode in close contact with DNA, whereas the L2 state is a major- or minor-groove tracking mode in loose contact with DNA and the L3 state is a two-dimensional freely diffusing mode during which hOGG1 is somewhat removed from the DNA surface (∼24 Å away from the surface). This three-state model well describes key experimental findings obtained from single-molecule and ensemble experiments and provides a unified molecular picture of the DNA lesion-searching mechanism of hOGG1.

Structure and dynamics of double-stranded DNA rotaxanes

A DNA rotaxane, with its unique mechanically interlocked architecture consisting of a circular DNA molecule threaded onto a linear DNA axle, holds promise as a fundamental component for nanoscale functional devices. Nevertheless, its structural and dynamic behaviors, essential for advancing molecular machinery, remain largely unexplored. Using extensive all-atom molecular dynamics simulations, we investigated the behaviors of double-stranded DNA (dsDNA) rotaxanes, concentrating on the effects of shape distortion induced by torsional stress in small circular dsDNA containing 70-90 base pairs. We analyzed structural characteristics, including shape, intermolecular distances, and tilt angles, while also exploring dynamic properties such as translational diffusion and toroidal rotation. Our results indicate that shape distortion brings the circular and linear dsDNA components into closer proximity and causes a slight increase in translational diffusion yet a minor decrease in toroidal rotation. Nevertheless, there is no apparent evidence of coupling between translation and rotation. Overall, the insights from this study indicate that such shape distortion does not significantly alter their structure and dynamics. This finding provides flexibility for the design of DNA rotaxanes in nanoscale applications.

How Do Nucleosome Dynamics Regulate Protein Search on DNA?

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

The membrane-actin linker ezrin acts as a sliding anchor

Protein linkages to filamentous (F)-actin provide the cell membrane with mechanical stability and support intricate membrane architectures. However, the actin cytoskeleton is highly dynamic and undergoes rapid changes in shape during cell motility and other processes. The molecular mechanisms that generate a mechanically robust yet fluid connection between the membrane and actin cytoskeleton remain poorly understood. Here, we adapted a single-molecule optical trap assay to examine how the prototypical membrane-actin linker ezrin acts to anchor F-actin to the cell membrane. We find that ezrin forms a complex that slides along F-actin over micrometer distances while resisting detachment by forces oriented perpendicular to the filament axis. The ubiquity of ezrin and analogous proteins suggests that sliding anchors such as ezrin may constitute an important but overlooked element in the construction of the actin cytoskeleton.

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.

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.

Sequence Preference and Initiator Promiscuity for De Novo DNA Synthesis by Terminal Deoxynucleotidyl Transferase

The untemplated activity of terminal deoxynucleotidyl transferase (TdT) represents its most appealing feature. Its use is well established in applications aiming for extension of a DNA initiator strand, but a more recent focus points to its potential in enzymatic de novo synthesis of DNA. Whereas its low substrate specificity for nucleoside triphosphates has been studied extensively, here we interrogate how the activity of TdT is modulated by the nature of the initiating strands, in particular their length, chemistry, and nucleotide composition. Investigation of full permutational libraries of mono- to pentamers of d-DNA, l-DNA, and 2′O-methyl-RNA of differing directionality immobilized to glass surfaces, and generated via photolithographic in situ synthesis, shows that the efficiency of extension strongly depends on the nucleobase sequence. We also show TdT being catalytically active on a non-nucleosidic substrate, hexaethylene glycol. These results offer new perspectives on constraints and strategies for de novo synthesis of DNA using TdT regarding the requirements for initiation of enzymatic generation of DNA.

Evolutionary advantage of a dissociative search mechanism in DNA mismatch repair

Protein complexes involved in DNA mismatch repair diffuse along dsDNA as sliding clamps in order to locate a hemimethylated incision site. They have been observed to use a dissociative mechanism, in which two proteins, while continuously remaining attached to the DNA, sometimes associate into a single complex sliding on the DNA and sometimes dissociate into two independently sliding proteins. Here, we study the probability that these complexes locate a given target site via a semi-analytic, Monte Carlo calculation that tracks the association and dissociation of the sliding complexes. We compare such probabilities to those obtained using a nondissociative diffusive scan in the space of physically realistic diffusion constants, hemimethylated site distances, and total search times to determine the regions in which dissociative searching is more or less efficient than nondissociative searching. We conclude that the dissociative search mechanism is advantageous in the majority of the physically realistic parameter space, suggesting that the dissociative search mechanism confers an evolutionary advantage.

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.

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.

Improved Parameterization of Protein-DNA Interactions for Molecular Dynamics Simulations of PCNA Diffusion on DNA

As the field of molecular dynamics simulation utilizing the force fields is moving toward more complex systems, the accuracy of intermolecular interactions has become a central issue of the field. Here, we quantitatively evaluate the accuracy of the protein-DNA interactions in AMBER and CHARMM force fields by comparing experimental and simulated diffusion coefficients of proliferating cell nuclear antigen. We find that both force fields underestimate diffusion coefficients by at least an order of magnitude because the interactions between basic amino acids and DNA phosphate groups are too attractive. Then, we propose Lennard-Jones parameters optimized using the experimental osmotic pressure data of model chemicals, by using which one can reproduce the experimental diffusion coefficients. Newly optimized parameters will have a broad impact on general protein-DNA interactions.

Protein Diffusion on Charged Biopolymers: DNA versus Microtubule

Protein diffusion in lower-dimensional spaces is used for various cellular functions. For example, sliding on DNA is essential for proteins searching for their target sites, and protein diffusion on microtubules is important for proper cell division and neuronal development. On the one hand, these linear diffusion processes are mediated by long-range electrostatic interactions between positively charged proteins and negatively charged biopolymers and have similar characteristic diffusion coefficients. On the other hand, DNA and microtubules have different structural properties. Here, using computational approaches, we studied the mechanism of protein diffusion along DNA and microtubules by exploring the diffusion of both protein types on both biopolymers. We found that DNA-binding and microtubule-binding proteins can diffuse on each other's substrates; however, the adopted diffusion mechanism depends on the molecular properties of the diffusing proteins and the biopolymers. On the protein side, only DNA-binding proteins can perform rotation-coupled diffusion along DNA, with this being due to their higher net charge and its spatial organization at the DNA recognition helix. By contrast, the lower net charge on microtubule-binding proteins enables them to diffuse more quickly than DNA-binding proteins on both biopolymers. On the biopolymer side, microtubules possess intrinsically disordered, negatively charged C-terminal tails that interact with microtubule-binding proteins, thus supporting their diffusion. Thus, although both DNA-binding and microtubule-binding proteins can diffuse on the negatively charged biopolymers, the unique molecular features of the biopolymers and of their natural substrates are essential for function.

Rational design of an XNA ligase through docking of unbound nucleic acids to toroidal proteins

Xenobiotic nucleic acids (XNA) are nucleic acid analogues not present in nature that can be used for the storage of genetic information. In vivo XNA applications could be developed into novel biocontainment strategies, but are currently limited by the challenge of developing XNA processing enzymes such as polymerases, ligases and nucleases. Here, we present a structure-guided modelling-based strategy for the rational design of those enzymes essential for the development of XNA molecular biology. Docking of protein domains to unbound double-stranded nucleic acids is used to generate a first approximation of the extensive interaction of nucleic acid processing enzymes with their substrate. Molecular dynamics is used to optimise that prediction allowing, for the first time, the accurate prediction of how proteins that form toroidal complexes with nucleic acids interact with their substrate. Using the Chlorella virus DNA ligase as a proof of principle, we recapitulate the ligase's substrate specificity and successfully predict how to convert it into an XNA-templated XNA ligase.

p15PAF binding to PCNA modulates the DNA sliding surface

p15^PAF is an oncogenic intrinsically disordered protein that regulates DNA replication and lesion bypass by interacting with the human sliding clamp PCNA. In the absence of DNA, p15^PAF traverses the PCNA ring via an extended PIP-box that contacts the sliding surface. Here, we probed the atomic-scale structure of p15^PAF–PCNA–DNA ternary complexes. Crystallography and MD simulations show that, when p15^PAF occupies two subunits of the PCNA homotrimer, DNA within the ring channel binds the unoccupied subunit. The structure of PCNA-bound p15^PAF in the absence and presence of DNA is invariant, and solution NMR confirms that DNA does not displace p15^PAF from the ring wall. Thus, p15^PAF reduces the available sliding surfaces of PCNA, and may function as a belt that fastens the DNA to the clamp during synthesis by the replicative polymerase (pol δ). This constraint, however, may need to be released for efficient DNA lesion bypass by the translesion synthesis polymerase (pol η). Accordingly, our biochemical data show that p15^PAF impairs primer synthesis by pol η–PCNA holoenzyme against both damaged and normal DNA templates. In light of our findings, we discuss the possible mechanistic roles of p15^PAF in DNA replication and suppression of DNA lesion bypass.

The Role of Noncognate Sites in the 1D Search Mechanism of EcoRI

A one-dimensional (1D) search is an essential step in DNA target recognition. Theoretical studies have suggested that the sequence dependence of 1D diffusion can help resolve the competing demands of a fast search and high target affinity, a conflict known as the speed-selectivity paradox. The resolution requires that the diffusion energy landscape is correlated with the underlying specific binding energies. In this work, we report observations of a 1D search by quantum dot-labeled EcoRI. Our data supports the view that proteins search DNA via rotation-coupled sliding over a corrugated energy landscape. We observed that whereas EcoRI primarily slides along DNA at low salt concentrations, at higher concentrations, its diffusion is a combination of sliding and hopping. We also observed long-lived pauses at genomic star sites, which differ by a single nucleotide from the target sequence. To reconcile these observations with prior biochemical and structural data, we propose a model of search in which the protein slides over a sequence-independent energy landscape during fast search but rapidly interconverts with a "hemispecific" binding mode in which a half site is probed. This half site interaction stabilizes the transition to a fully specific mode of binding, which can then lead to target recognition.

Nucleosome Crowding in Chromatin Slows the Diffusion but Can Promote Target Search of Proteins

Dynamics of nuclear proteins in crowded chromatin has only been poorly understood. Here, we address the diffusion, target search, and structural dynamics of three proteins in a model chromatin using coarse-grained molecular simulations run on the K computer. We prepared two structures of chromatin made of 20 nucleosomes with different nucleosome densities and investigated dynamics of two transcription factors, HMGB1 and p53, and one signaling protein, ERK, embedded in the chromatin. We found fast and normal diffusion of the nuclear proteins in the low-density chromatins and slow and subdiffusional movements in the high-density chromatin. The diffusion of the largest transcription factor, p53, is slowed by high-density chromatin most markedly. The on rates and off rates for DNA binding are increased and decreased, respectively, in the high-density chromatin. To our surprise, the DNA sequence search was faster in chromatin with high nucleosome density, though the diffusion is slower. We also found that the three nuclear proteins preferred to bind on the linker DNA and the entry and exit regions of nucleosomal DNA. In addition to these regions, HMGB1 and p53 also bound to the dyad.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

Phosphorylation meets DNA mismatch repair

DNA mismatch repair (MMR) is a highly conserved process and ensures the removal of mispaired DNA bases and insertion-deletion loops right after replication. For this, a MutSα or MutSβ protein complex recognizes the DNA damage, MutLα nicks the erroneous strand, exonuclease 1 removes the wrong nucleotides, DNA polymerase δ refills the gap and DNA ligase I joins the fragments to seal the nicks and complete the repair process. The failure to accomplish these functions is associated with higher mutation rates and may lead to cancer, which highlights the importance of MMR by the maintenance of genomic stability. The post-replicative MMR implies that involved proteins are regulated at several levels, including posttranslational modifications (PTMs). Phosphorylation is one of the most common and major PTMs. Suitable with its regulatory force phosphorylation was shown to influence MMR factors thereby adjusting eukaryotic MMR activity. In this review, we summarized the current knowledge of the role of phosphorylation of MMR process involved proteins and their functional relevance.