Dynamic Strategies for Target-Site Search by DNA-Binding Proteins

Gene regulatory proteins find their target sites on DNA remarkably quickly; the experimental binding rate for lac repressor is orders-of-magnitude higher than predicted by free diffusion alone. It has been proposed that nonspecific binding aids the search by allowing proteins to slide and hop along DNA. We develop a reaction-diffusion theory of protein translocation that accounts for transport both on and off the strand and incorporates the physical conformation of DNA. For linear DNA modeled as a wormlike chain, the distribution of hops available to a protein exhibits long, power-law tails that make the long-time displacement along the strand superdiffusive. Our analysis predicts effective superdiffusion coefficients for given nonspecific binding and unbinding rate parameters. Translocation rate exhibits a maximum at intermediate values of the binding rate constant, while search efficiency is optimized at larger binding rate constant values. Thus, our theory predicts a region of values of the nonspecific binding and unbinding rate parameters that balance the protein translocation rate and the efficiency of the search. Published data for several proteins falls within this predicted region of parameter values.     

Specific and non-specific interactions of ParB with DNA: implications for chromosome segregation

The segregation of many bacterial chromosomes is dependent on the interactions of ParB proteins with centromere-like DNA sequences called parS that are located close to the origin of replication. In this work, we have investigated the binding of Bacillus subtilis ParB to DNA in vitro using a variety of biochemical and biophysical techniques. We observe tight and specific binding of a ParB homodimer to the parS sequence. Binding of ParB to non-specific DNA is more complex and displays apparent positive co-operativity that is associated with the formation of larger, poorly defined, nucleoprotein complexes. Experiments with magnetic tweezers demonstrate that non-specific binding leads to DNA condensation that is reversible by protein unbinding or force. The condensed DNA structure is not well ordered and we infer that it is formed by many looping interactions between neighbouring DNA segments. Consistent with this view, ParB is also able to stabilize writhe in single supercoiled DNA molecules and to bridge segments from two different DNA molecules in trans. The experiments provide no evidence for the promotion of non-specific DNA binding and/or condensation events by the presence of parS sequences. The implications of these observations for chromosome segregation are discussed.     

Conformational selection in protein binding and function

Protein binding and function often involves conformational changes. Advanced nuclear magnetic resonance (NMR) experiments indicate that these conformational changes can occur in the absence of ligand molecules (or with bound ligands), and that the ligands may “select” protein conformations for binding (or unbinding). In this review, we argue that this conformational selection requires transition times for ligand binding and unbinding that are small compared to the dwell times of proteins in different conformations, which is plausible for small ligand molecules. Such a separation of timescales leads to a decoupling and temporal ordering of binding/unbinding events and conformational changes. We propose that conformational‐selection and induced‐change processes (such as induced fit) are two sides of the same coin, because the temporal ordering is reversed in binding and unbinding direction. Conformational‐selection processes can be characterized by a conformational excitation that occurs prior to a binding or unbinding event, while induced‐change processes exhibit a characteristic conformational relaxation that occurs after a binding or unbinding event. We discuss how the ordering of events can be determined from relaxation rates and effective on‐ and off‐rates determined in mixing experiments, and from the conformational exchange rates measured in advanced NMR or single‐molecule fluorescence resonance energy transfer experiments. For larger ligand molecules such as peptides, conformational changes and binding events can be intricately coupled and exhibit aspects of conformational‐selection and induced‐change processes in both binding and unbinding direction.     

Escape of a Small Molecule from Inside T4 Lysozyme by Multiple Pathways

The T4 lysozyme L99A mutant is often used as a model system to study small-molecule binding to proteins, but pathways for ligand entry and exit from the buried binding site and the associated protein conformational changes have not been fully resolved. Here, molecular dynamics simulations were employed to model benzene exit from its binding cavity using the weighted ensemble (WE) approach to enhance sampling of low-probability unbinding trajectories. Independent WE simulations revealed four pathways for benzene exit, which correspond to transient tunnels spontaneously formed in previous simulations of apo T4 lysozyme. Thus, benzene unbinding occurs through multiple pathways partially created by intrinsic protein structural fluctuations. Motions of several α-helices and side chains were involved in ligand escape from metastable microstates. WE simulations also provided preliminary estimates of rate constants for each exit pathway. These results complement previous works and provide a semiquantitative characterization of pathway heterogeneity for binding of small molecules to proteins.     

Concentration-dependent exchange accelerates turnover of proteins bound to double-stranded DNA

The multistep kinetics through which DNA-binding proteins bind their targets are heavily studied, but relatively little attention has been paid to proteins leaving the double helix. Using single-DNA stretching and fluorescence detection, we find that sequence-neutral DNA-binding proteins Fis, HU and NHP6A readily exchange with themselves and with each other. In experiments focused on the Escherichia coli nucleoid-associated protein Fis, only a small fraction of protein bound to DNA spontaneously dissociates into protein-free solution. However, if Fis is present in solution, we find that a concentration-dependent exchange reaction occurs which turns over the bound protein, with a rate of k_exch = 6 × 10⁴ M⁻¹s⁻¹. The bacterial DNA-binding protein HU and the yeast HMGB protein NHP6A display the same phenomenon of protein in solution accelerating dissociation of previously bound labeled proteins as exchange occurs. Thus, solvated proteins can play a key role in facilitating removal and renewal of proteins bound to the double helix, an effect that likely plays a major role in promoting the turnover of proteins bound to DNA in vivo and, therefore, in controlling the dynamics of gene regulation.     

Entamoeba histolytica Dmc1 Catalyzes Homologous DNA Pairing and Strand Exchange That Is Stimulated by Calcium and Hop2-Mnd1

Meiosis depends on homologous recombination (HR) in most sexually reproducing organisms. Efficient meiotic HR requires the activity of the meiosis-specific recombinase, Dmc1. Previous work shows Dmc1 is expressed in Entamoeba histolytica, a eukaryotic parasite responsible for amoebiasis throughout the world, suggesting this organism undergoes meiosis. Here, we demonstrate Dmc1 protein is expressed in E. histolytica. We show that purified ehDmc1 forms presynaptic filaments and catalyzes ATP-dependent homologous DNA pairing and DNA strand exchange over at least several thousand base pairs. The DNA pairing and strand exchange activities are enhanced by the presence of calcium and the meiosis-specific recombination accessory factor, Hop2-Mnd1. In combination, calcium and Hop2-Mnd1 dramatically increase the rate of DNA strand exchange activity of ehDmc1. The biochemical system described herein provides a basis on which to better understand the role of ehDmc1 and other HR proteins in E. histolytica.     

Operational Principles for the Dynamics of the In Vitro ParA-ParB System

In many bacteria the ParA-ParB protein system is responsible for actively segregating DNA during replication. ParB proteins move by interacting with DNA bound ParA-ATP, stimulating their unbinding by catalyzing hydrolysis, that leads to rectified motion due to the creation of a wake of depleted ParA. Recent in vitro experiments have shown that a ParB covered magnetic bead can move with constant speed over a DNA covered substrate that is bound by ParA. It has been suggested that the formation of a gradient in ParA leads to diffusion-ratchet like motion of the ParB bead but how it forms and generates a force is still a matter of exploration. Here we develop a deterministic model for the in vitro ParA-ParB system and show that a ParA gradient can spontaneously form due to any amount of initial spatial noise in bound ParA. The speed of the bead is independent of this noise but depends on the ratio of the range of ParA-ParB force on the bead to that of removal of surface bound ParA by ParB. We find that at a particular ratio the speed attains a maximal value. We also consider ParA rebinding (including cooperativity) and ParA surface diffusion independently as mechanisms for ParA recovery on the surface. Depending on whether the DNA covered surface is undersaturated or saturated with ParA, we find that the bead can accelerate persistently or potentially stall. Our model highlights key requirements of the ParA-ParB driving force that are necessary for directed motion in the in vitro system that may provide insight into the in vivo dynamics of the ParA-ParB system.     

RAD54 controls access to the invading 3′-OH end after RAD51-mediated DNA strand invasion in homologous recombination in Saccharomyces cerevisiae

Rad51 is a key protein in homologous recombination performing homology search and DNA strand invasion. After DNA strand exchange Rad51 protein is stuck on the double-stranded heteroduplex DNA product of DNA strand invasion. This is a problem, because DNA polymerase requires access to the invading 3′-OH end to initiate DNA synthesis. Here we show that, the Saccharomyces cerevisiae dsDNA motor protein Rad54 solves this problem by dissociating yeast Rad51 protein bound to the heteroduplex DNA after DNA strand invasion. The reaction required species-specific interaction between both proteins and the ATPase activity of Rad54 protein. This mechanism rationalizes the in vivo requirement of Rad54 protein for the turnover of Rad51 foci and explains the observed dependence of the transition from homologous pairing to DNA synthesis on Rad54 protein in vegetative and meiotic yeast cells.     

The differential extension in dsDNA bound to Rad51 filaments may play important roles in homology recognition and strand exchange

RecA and Rad51 proteins play an important role in DNA repair and homologous recombination. For RecA, X-ray structure information and single molecule force experiments have indicated that the differential extension between the complementary strand and its Watson–Crick pairing partners promotes the rapid unbinding of non-homologous dsDNA and drives strand exchange forward for homologous dsDNA. In this work we find that both effects are also present in Rad51 protein. In particular, pulling on the opposite termini (3′ and 5′) of one of the two DNA strands in a dsDNA molecule allows dsDNA to extend along non-homologous Rad51-ssDNA filaments and remain stably bound in the extended state, but pulling on the 3′5′ ends of the complementary strand reduces the strand-exchange rate for homologous filaments. Thus, the results suggest that differential extension is also present in dsDNA bound to Rad51. The differential extension promotes rapid recognition by driving the swift unbinding of dsDNA from non-homologous Rad51-ssDNA filaments, while at the same time, reducing base pair tension due to the transfer of the Watson–Crick pairing of the complementary strand bases from the highly extended outgoing strand to the slightly less extended incoming strand, which drives strand exchange forward.     

Structural and Functional Characterisation of a Conserved Archaeal RadA Paralog with Antirecombinase Activity

DNA recombinases (RecA in bacteria, Rad51 in eukarya and RadA in archaea) catalyse strand exchange between homologous DNA molecules, the central reaction of homologous recombination, and are among the most conserved DNA repair proteins known. RecA is the sole protein responsible for this reaction in bacteria, whereas there are several Rad51 paralogs that cooperate to catalyse strand exchange in eukaryotes. All archaea have at least one (and as many as four) RadA paralog, but their function remains unclear. Herein, we show that the three RadA paralogs encoded by the Sulfolobus solfataricus genome are expressed under normal growth conditions and are not UV inducible. We demonstrate that one of these proteins, Sso2452, which is representative of the large archaeal RadC subfamily of archaeal RadA paralogs, functions as an ATPase that binds tightly to single-stranded DNA. However, Sso2452 is not an active recombinase in vitro and inhibits D-loop formation by RadA. We present the high-resolution crystal structure of Sso2452, which reveals key structural differences from the canonical RecA family recombinases that may explain its functional properties. The possible roles of the archaeal RadA paralogs in vivo are discussed.     

Increase in incidence of chromosome instability and non-conservative recombination between repeats in Saccharomyces cerevisiae hpr1Δ strains

Null hprl Δ strains show a large increase (up to 2000-fold) over wild type in the frequency of occurrence of deletions between direct repeats on three different chromosomes. However, we show that hprl Δ mutations have little or no effect on reciprocal exchange, gene conversion or unequal sister chromatid exchange, as determined using intrachromosomal, interchromosomal and plasmid-chromosome assay systems. A novel intrachromosomal recombination system has allowed us to determine that over 95% of deletions in hpr1 Δ strains do not occur by reciprocal exchange. On the other hand, hpr1 Δ strains show chromosome loss frequencies of up to 100 times the wild-type level. Our results suggest that yeast cells have a very efficient non-conservative recombination mechanism, dependent on RADI and RAD52, that causes deletions between direct DNA repeats, and this mechanism is strongly stimulated in hpr1 Δ strains. The results indicate that the Hpr1 protein is required for stability of DNA repeats and chromosomes. We propose that in the absence of the Hprl protein the cell destabilizes the genome by allowing the initiation of events that lead to deletions of sequences between repeats, and to chromosome instability. We discuss the roles that proteins such as Hprl have in maintaining direct repeats and in preventing non-conservative recombination and consider the importance of these functions for chromosome stability.     

Staphylococcus aureus protein SAUGI acts as a uracil-DNA glycosylase inhibitor

DNA mimic proteins are unique factors that control the DNA binding activity of target proteins by directly occupying their DNA binding sites. The extremely divergent amino acid sequences of the DNA mimics make these proteins hard to predict, and although they are likely to be ubiquitous, to date, only a few have been reported and functionally analyzed. Here we used a bioinformatic approach to look for potential DNA mimic proteins among previously reported protein structures. From ∼14 candidates, we selected the Staphylococcus conserved hypothetical protein SSP0047, and used proteomic and structural approaches to show that it is a novel DNA mimic protein. In Staphylococcus aureus, we found that this protein acts as a uracil-DNA glycosylase inhibitor, and therefore named it S. aureus uracil-DNA glycosylase inhibitor (SAUGI). We also determined and analyzed the complex structure of SAUGI and S. aureus uracil-DNA glycosylase (SAUDG). Subsequent BIAcore studies further showed that SAUGI has a high binding affinity to both S. aureus and human UDG. The two uracil-DNA glycosylase inhibitors (UGI and p56) previously known to science were both found in Bacillus phages, and this is the first report of a bacterial DNA mimic that may regulate SAUDG’s functional roles in DNA repair and host defense.     

Inferring Diffusion Dynamics from FCS in Heterogeneous Nuclear Environments

Fluorescence correlation spectroscopy (FCS) is a noninvasive technique that probes the diffusion dynamics of proteins down to single-molecule sensitivity in living cells. Critical mechanistic insight is often drawn from FCS experiments by fitting the resulting time-intensity correlation function, G(t), to known diffusion models. When simple models fail, the complex diffusion dynamics of proteins within heterogeneous cellular environments can be fit to anomalous diffusion models with adjustable anomalous exponents. Here, we take a different approach. We use the maximum entropy method to show—first using synthetic data—that a model for proteins diffusing while stochastically binding/unbinding to various affinity sites in living cells gives rise to a G(t) that could otherwise be equally well fit using anomalous diffusion models. We explain the mechanistic insight derived from our method. In particular, using real FCS data, we describe how the effects of cell crowding and binding to affinity sites manifest themselves in the behavior of G(t). Our focus is on the diffusive behavior of an engineered protein in 1) the heterochromatin region of the cell’s nucleus as well as 2) in the cell’s cytoplasm and 3) in solution. The protein consists of the basic region-leucine zipper (BZip) domain of the CCAAT/enhancer-binding protein (C/EBP) fused to fluorescent proteins.     

Ligand-Induced Changes of the Apparent Transition-State Position in?Mechanical Protein Unfolding

Force-spectroscopic measurements of ligand-receptor systems and the unfolding/folding of nucleic acids or proteins reveal information on the underlying energy landscape along the pulling coordinate. The slope Δx^‡ of the force-dependent unfolding/unbinding rates is interpreted as the distance from the folded/bound state to the transition state for unfolding/unbinding and, hence, often related to the mechanical compliance of the sample molecule. Here we show that in ligand-binding proteins, the experimentally inferred Δx^‡ can depend on the ligand concentration, unrelated to changes in mechanical compliance. We describe the effect in single-molecule, force-spectroscopy experiments of the calcium-binding protein calmodulin and explain it in a simple model where mechanical unfolding and ligand binding occur on orthogonal reaction coordinates. This model predicts changes in the experimentally inferred Δx^‡, depending on ligand concentration and the associated shift of the dominant barrier between the two reaction coordinates. We demonstrate quantitative agreement between experiments and simulations using a realistic six-state kinetic scheme using literature values for calcium-binding kinetics and affinities. Our results have important consequences for the interpretation of force-spectroscopic data of ligand-binding proteins.     

Rad51 gain-of-function mutants that exhibit high affinity DNA binding cause DNA damage sensitivity in the absence of Srs2

We previously identified several rad51 gain-of-function alleles that partially suppress the requirement for RAD55 and RAD57 in DNA repair. To gain further insight into the mechanism of action of these alleles, we compared the activities of Rad51-V328A, Rad51-P339S and Rad51-I345T with wild-type Rad51, for DNA binding, filament stability, strand exchange and interaction with the antirecombinase helicase, Srs2. These alleles were chosen because they show the highest activity in suppression of ionizing radiation sensitivity of the rad57 mutant, and Val 328 and Ile 345 are conserved in the human Rad51 protein. All three mutant proteins exhibited higher affinity for single-stranded DNA (ssDNA) and showed more robust strand exchange activity with oligonucleotide substrates than wild-type Rad51, with the Rad51-I345T and Rad51-V328A proteins displaying higher activity than Rad51-P339S. However, the Srs2 antirecombinase was able to disrupt Rad51–ssDNA complexes formed with all the mutant proteins. In vivo, the rad51-I345T mutant strain exhibited high resistance to methyl methane sulfonate that was dependent on functional SRS2. These results suggest the Srs2 translocase is able to disrupt Rad51–ssDNA complexes at stalled replication forks, but in the absence of Srs2 the enhanced DNA binding of the Rad51-I345T protein is detrimental to cell survival.     

Human Rad51 filaments on double- and single-stranded DNA: correlating regular and irregular forms with recombination function

Recombinase proteins assembled into helical filaments on DNA are believed to be the catalytic core of homologous recombination. The assembly, disassembly and dynamic rearrangements of this structure must drive the DNA strand exchange reactions of homologous recombination. The sensitivity of eukaryotic recombinase activity to reaction conditions in vitro suggests that the status of bound nucleotide cofactors is important for function and possibly for filament structure. We analyzed nucleoprotein filaments formed by the human recombinase Rad51 in a variety of conditions on double-stranded and single-stranded DNA by scanning force microscopy. Regular filaments with extended double-stranded DNA correlated with active in vitro recombination, possibly due to stabilizing the DNA products of these assays. Though filaments formed readily on single-stranded DNA, they were very rarely regular structures. The irregular structure of filaments on single-stranded DNA suggests that Rad51 monomers are dynamic in filaments and that regular filaments are transient. Indeed, single molecule force spectroscopy of Rad51 filament assembly and disassembly in magnetic tweezers revealed protein association and disassociation from many points along the DNA, with kinetics different from those of RecA. The dynamic rearrangements of proteins and DNA within Rad51 nucleoprotein filaments could be key events driving strand exchange in homologous recombination.     

Regulation of Deinococcus radiodurans RecA Protein Function via Modulation of Active and Inactive Nucleoprotein Filament States

The RecA protein of Deinococcus radiodurans (DrRecA) has a central role in genome reconstitution after exposure to extreme levels of ionizing radiation. When bound to DNA, filaments of DrRecA protein exhibit active and inactive states that are readily interconverted in response to several sets of stimuli and conditions. At 30 °C, the optimal growth temperature, and at physiological pH 7.5, DrRecA protein binds to double-stranded DNA (dsDNA) and forms extended helical filaments in the presence of ATP. However, the ATP is not hydrolyzed. ATP hydrolysis of the DrRecA-dsDNA filament is activated by addition of single-stranded DNA, with or without the single-stranded DNA-binding protein. The ATPase function of DrRecA nucleoprotein filaments thus exists in an inactive default state under some conditions. ATPase activity is thus not a reliable indicator of DNA binding for all bacterial RecA proteins. Activation is effected by situations in which the DNA substrates needed to initiate recombinational DNA repair are present. The inactive state can also be activated by decreasing the pH (protonation of multiple ionizable groups is required) or by addition of volume exclusion agents. Single-stranded DNA-binding protein plays a much more central role in DNA pairing and strand exchange catalyzed by DrRecA than is the case for the cognate proteins in Escherichia coli. The data suggest a mechanism to enhance the efficiency of recombinational DNA repair in the context of severe genomic degradation in D. radiodurans.     

Extending Bell's model: how force transducer stiffness alters measured unbinding forces and kinetics of molecular complexes

Forced unbinding of complementary macromolecules such as ligand-receptor complexes can reveal energetic and kinetic details governing physiological processes ranging from cellular adhesion to drug metabolism. Although molecular-level experiments have enabled sampling of individual ligand-receptor complex dissociation events, disparities in measured unbinding force F_R among these methods lead to marked variation in inferred binding energetics and kinetics at equilibrium. These discrepancies are documented for even the ubiquitous ligand-receptor pair, biotin-streptavidin. We investigated these disparities and examined atomic-level unbinding trajectories via steered molecular dynamics simulations, as well as via molecular force spectroscopy experiments on biotin-streptavidin. In addition to the well-known loading rate dependence of F_R predicted by Bell's model, we find that experimentally accessible parameters such as the effective stiffness of the force transducer k can significantly perturb the energy landscape and the apparent unbinding force of the complex for sufficiently stiff force transducers. Additionally, at least 20% variation in unbinding force can be attributed to minute differences in initial atomic positions among energetically and structurally comparable complexes. For force transducers typical of molecular force spectroscopy experiments and atomistic simulations, this energy barrier perturbation results in extrapolated energetic and kinetic parameters of the complex that depend strongly on k. We present a model that explicitly includes the effect of k on apparent unbinding force of the ligand-receptor complex, and demonstrate that this correction enables prediction of unbinding distances and dissociation rates that are decoupled from the stiffness of actual or simulated molecular linkers.