Target search of N sliding proteins on a DNA

At low to moderate ambient salt concentrations, DNA-binding proteins bind relatively tightly to DNA, and only very rarely detach. Intersegmental transfer due to DNA-looping can be excluded by applying an external pulling force to the DNA molecule. Under such conditions, we explore the targeting dynamics of N proteins sliding diffusively along DNA in search of their specific target sequence. At lower densities of binding proteins, we find a reduction of the characteristic search time proportional to N(-2), with corrections at higher concentrations. Rates for detachment and attachment of binding proteins are incorporated in the model. Our findings are in agreement with recent single molecule studies in the presence of bacteriophage T4 gene 32 protein for which the unbinding rate is much lower than the specific binding rate.     

Searching DNA via a “Monkey Bar” Mechanism: The Significance of Disordered Tails

The search through nonspecific DNA for a specific site by proteins is known to be facilitated by sliding, hopping, and intersegment transfer between separate DNA strands, yet the driving forces of these protein dynamics from the molecular perspective are unclear. In this study, molecular features of the DNA search mechanism were explored for three homologous proteins (the HoxD9, Antp, and NK-2 homeodomains) using a simple computational model in which protein-DNA interactions are represented solely by electrostatic forces. In particular, we studied the impact that disordered N-terminal tails (N-tails), which are more common in DNA-binding proteins than in other proteins, have on the efficiency of DNA search. While the three homeodomain proteins were found to use similar binding interfaces in specific and nonspecific interactions with DNAs, their different electrostatic potentials affect the nature of their sliding dynamics. The different lengths and net charges of the N-tails of the homeodomains affect their motion along the DNA. The presence of an N-tail increases sliding propensity but slows linear diffusion along the DNA. When the search is performed in the presence of two parallel DNA molecules, a direct transfer, which is facilitated by the protein tail, from one nonspecific DNA to another occurs. The tailed proteins jump between two DNA molecules through an intermediate in which the recognition helix of the protein is adsorbed to one DNA fragment and the N-tail is adsorbed to the second, suggesting a “monkey bar” mechanism. Our study illustrates how the molecular architecture of proteins controls the efficiency of DNA scanning.     

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.     

Accelerated Search Kinetics Mediated by Redox Reactions of DNA Repair Enzymes

A charge transport (CT) mechanism has been proposed in several articles to explain the localization of base excision repair (BER) enzymes to lesions on DNA. The CT mechanism relies on redox reactions of iron-sulfur cofactors that modify the enzyme's binding affinity. These redox reactions are mediated by the DNA strand and involve the exchange of electrons between BER enzymes along DNA. We propose a mathematical model that incorporates enzyme binding/unbinding, electron transport, and enzyme diffusion along DNA. Analysis of our model within a range of parameter values suggests that the redox reactions can increase desorption of BER enzymes not already bound to lesions, allowing the enzymes to be recycled—thus accelerating the overall search process. This acceleration mechanism is most effective when enzyme copy numbers and enzyme diffusivity along the DNA are small. Under such conditions, we find that CT BER enzymes find their targets more quickly than simple passive enzymes that simply attach to the DNA without desorbing.     

Human Rad54 protein stimulates DNA strand exchange activity of hRad51 protein in the presence of Ca2+

Rad51 and Rad54 proteins play a key role in homologous recombination in eukaryotes. Recently, we reported that Ca2+ is required in vitro for human Rad51 protein to form an active nucleoprotein filament that is important for the search of homologous DNA and for DNA strand exchange, two critical steps of homologous recombination. Here we find that Ca2+ is also required for hRad54 protein to effectively stimulate DNA strand exchange activity of hRad51 protein. This finding identifies Ca2+ as a universal cofactor of DNA strand exchange promoted by mammalian homologous recombination proteins in vitro. We further investigated the hRad54-dependent stimulation of DNA strand exchange. The mechanism of stimulation appeared to include specific interaction of hRad54 protein with the hRad51 nucleoprotein filament. Our results show that hRad54 protein significantly stimulates homology-independent coaggregation of dsDNA with the filament, which represents an essential step of the search for homologous DNA. The results obtained indicate that hRad54 protein serves as a dsDNA gateway for the hRad51-ssDNA filament, promoting binding and an ATP hydrolysis-dependent translocation of dsDNA during the search for homologous sequences.     

Protein Sliding along DNA: Dynamics and Structural Characterization

Efficient search of DNA by proteins is fundamental to the control of cellular regulatory processes. It is currently believed that protein sliding, hopping, and transfer between adjacent DNA segments, during which the protein nonspecifically interacts with DNA, are central to the speed of their specific recognition. In this study, we focused on the structural and dynamic features of proteins when they scan the DNA. Using a simple computational model that represents protein-DNA interactions by electrostatic forces, we identified that the protein makes use of identical binding interfaces for both nonspecific and specific DNA interactions. Accordingly, in its one-dimensional diffusion along the DNA, the protein is bound at the major groove and performs a helical motion, which is stochastic and driven by thermal diffusion. A microscopic structural insight into sliding from our model, which is governed by electrostatic forces, corroborates previous experimental studies suggesting that the active site of some regulatory proteins continually faces the interior of the DNA groove while sliding along sugar-phosphate rails. The diffusion coefficient of spiral motion along the major groove of the DNA is not affected by salt concentration, but the efficiency of the search can be significantly enhanced by increasing salt concentration due to a larger number of hopping events. We found that the most efficient search comprises ∼ 20% sliding along the DNA and ∼ 80% hopping and three-dimensional diffusion. The presented model that captures various experimental features of facilitated diffusion has the potency to address other questions regarding the nature of DNA search, such as the sliding characteristics of oligomeric and multidomain DNA-binding proteins that are ubiquitous in the cell.     

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.     

Analytic binding isotherms describing competitive interactions of a protein ligand with specific and nonspecific sites on the same DNA oligomer

Many studies of specific protein-nucleic acid binding use short oligonucleotides or restriction fragments, in part to minimize the potential for nonspecific binding of the protein. However, when the specificity ratio is low, multiple nonspecifically bound proteins may occupy the region of DNA corresponding to one specific site; this situation was encountered in our recent calorimetric study of binding of integration host factor (IHF) protein to its specific 34-bp H' DNA site. Here, beginning from the analytical McGhee and von Hippel infinite-lattice nonspecific binding isotherm, we derive a novel analytic isotherm for nonspecific binding of a ligand to a finite lattice. This isotherm is an excellent approximation to the exact factorial-based Epstein finite lattice isotherm even for short lattices and therefore is of great practical significance for analysis of experimental data and for analytic theory. Using this isotherm, we develop an analytic treatment of the competition between specific and nonspecific binding of a large ligand to the same finite lattice (i.e., DNA oligomer) containing one specific and multiple overlapping nonspecific binding sites. Analysis of calorimetric data for IHF-H' DNA binding using this treatment yields enthalpies and binding constants for both specific and nonspecific binding and the nonspecific site size. This novel analysis demonstrates the potential contribution of nonspecific binding to the observed thermodynamics of specific binding, even with very short DNA oligomers, and the need for reverse (constant protein) titrations or titrations with nonspecific DNA to resolve specific and nonspecific contributions. The competition treatment is useful in analyzing low-specificity systems, including those where specificity is weakened by mutations or the absence of specificity factors.     

A model for the mediation of processivity of DNA-targeting proteins by nonspecific binding: dependence on DNA length and presence of obstacles

A physical and mathematical model is presented to explain processivity of proteins on DNA. In this model, a DNA-targeting protein such as a restriction enzyme can diffuse to the DNA surface and nonspecifically bind to it. Once on the DNA surface it will either move along the DNA or equilibrate with the surrounding region. Owing to the nonspecific binding, the search for a specific site on the DNA occurs in a reduced dimensionality, and the protein appears processive when moving from one specific site to another. The simplest version of this nonspecific-binding-facilitated diffusion model is solved and the results quantitatively explain experimentally observed dependence of the processivity ratio on the intervening DNA length between two specific sites.     

The core and carboxyl-terminal domains of the integrase protein of human immunodeficiency virus type 1 each contribute to nonspecific DNA binding.

The integrase protein of human immunodeficiency virus type 1 removes two nucleotides from the 3' ends of reverse-transcribed human immunodeficiency virus type 1 DNA (3' processing) and covalently inserts the processed ends into a target DNA (DNA strand transfer). Mutant integrase proteins that lack the amino-and/or carboxyl-terminal domains are incapable of catalyzing 3' processing and DNA strand transfer but are competent for an apparent reversal of the DNA strand transfer reaction (disintegration) in vitro. Here, we investigate the binding of integrase to DNA by UV cross-linking. Cross-linked complexes form with a variety of DNA substrates independent of the presence of divalent metal ion. Analysis with amino- and carboxyl-terminal deletion mutant proteins shows that residues 213 to 266 of the 288-residue protein are required for efficient cross-linking in the absence of divalent metal ion. Carboxyl-terminal deletion mutants that lack this region efficiently cross-link only to the branched disintegration DNA substrate, and this reaction is dependent on the presence of metal ion. Both the core and C-terminal domains of integrase therefore contribute to nonspecific DNA binding.     

Stopped-Flow Fluorescence Kinetic Study of Protein Sliding and Intersegment Transfer in the Target DNA Search Process

Kinetic characterizations of protein translocation on DNA are nontrivial because the simultaneous presence of multiple different mechanisms makes it difficult to extract the information specific to a particular translocation mechanism. In this study, we have developed new approaches for the kinetic investigations of proteins' sliding and intersegment transfer (also known as “direct transfer”) in the target DNA search process. Based on the analytical expression of the mean search time for the discrete-state stochastic model, we derived analytical forms of the apparent rate constant k_app for protein-target association in systems involving competitor DNA and the intersegment transfer mechanism. Our analytical forms of k_app facilitate the experimental determination of the kinetic rate constants for intersegment transfer and sliding in the target association process. Using stopped-flow fluorescence data for the target association kinetics along with the analytical forms of k_app, we have studied the translocation of the Egr-1 zinc-finger protein in the target DNA association process. Sliding was analyzed using the DNA-length-dependent k_app data. Using the dependence of k_app on the concentration of competitor DNA, we determined the second-order rate constant for intersegment transfer. Our results indicate that a major pathway in the target association process for the Egr-1 zinc-finger protein is the one involving intersegment transfer to a nonspecific site and the subsequent sliding to the target.     

Two steps forward,one step back: Determining XPD helicase mechanism by single-molecule fluorescence and high-resolution optical tweezers

XPD-like helicases constitute a prominent DNA helicase family critical for many aspects of genome maintenance. These enzymes share a unique structural feature, an auxiliary domain stabilized by an iron-sulphur (FeS) cluster, and a 5′–3′ polarity of DNA translocation and duplex unwinding. Biochemical analyses alongside two single-molecule approaches, total internal reflection fluorescence microscopy and high-resolution optical tweezers, have shown how the unique structural features of XPD helicase and its specific patterns of substrate interactions tune the helicase for its specific cellular function and shape its molecular mechanism. The FeS domain forms a duplex separation wedge and contributes to an extended DNA binding site. Interactions within this site position the helicase in an orientation to unwind the duplex, control the helicase rate, and verify the integrity of the translocating strand. Consistent with its cellular role, processivity of XPD is limited and is defined by an idiosyncratic stepping kinetics. DNA duplex separation occurs in single base pair steps punctuated by frequent backward steps and conformational rearrangements of the protein–DNA complex. As such, the helicase in isolation mainly stabilizes spontaneous base pair opening and exhibits a limited ability to unwind stable DNA duplexes. The presence of a cognate ssDNA binding protein converts XPD into a vigorous helicase by destabilizing the upstream dsDNA as well as by trapping the unwound strands. Remarkably, the two proteins can co-exist on the same DNA strand without competing for binding. The current model of the XPD unwinding mechanism will be discussed along with possible modifications to this mechanism by the helicase interacting partners and unique features of such bio-medically important XPD-like helicases as FANCJ (BACH1), RTEL1 and CHLR1 (DDX11).     

Differential ATP requirements distinguish the DNA translocation and DNA unwinding activities of the Escherichia coli PRI A protein

The Escherichia coli primosome is a mobile multiprotein DNA replication-priming apparatus that assembles at a specific site (termed a primosome assembly site (PAS] on single-stranded DNA-binding protein-coated single-stranded DNA. The PRI A protein (factor Y, protein n') is a PAS sequence-specific (d)ATPase as well as a DNA helicase and is believed to direct the assembly of the primosome at a PAS. In this report, the PRI A DNA helicase reaction is dissected in vitro, by use of a strand displacement assay, into three steps with distinct ATP requirements. First, the PRI A protein gains entry to the DNA via an ATP-independent, PAS sequence-specific binding event. Second, the PRI A protein translocates along the single-stranded DNA in the 3'----5' direction at a maximal rate of 90 nucleotides/s. DNA translocation requires ATP hydrolysis. The ATP concentration required to support half of the maximal translocation rate is 100 microM, which is identical to the Km for ATP of the PRI A protein DNA-dependent ATPase activity. Finally, the PRI A protein unwinds duplex DNA. The ATP concentration required for duplex DNA unwinding depends upon the length of the duplex region to be unwound. Displacement of a 24-nucleotide long oligomer required no more ATP than that required for the translocation of PRI A protein along single-stranded DNA, whereas displacement of a 390-nucleotide long DNA fragment required a 10-fold higher concentration of ATP than that required for oligomer displacement.     

Mechanism of translocation and kinetics of DNA unwinding by the helicase RecG

RecG is a DNA helicase involved in the repair of damage at a replication fork and catalyzes the reversal of the fork to a point beyond the damage in the template strand. It unwinds duplex DNA in reactions that are coupled to ATP hydrolysis. The kinetic mechanism of duplex DNA unwinding by RecG was analyzed using a quantitative fluorescence assay based on the process of contact quenching between Cy3 and Dabcyl groups attached to synthetic three-way DNA junctions. The data show that the protein moves at a rate of 26 bp s(-1) along the duplex DNA during the unwinding process. RecG ATPase activity during translocation indicates a constant rate of 7.6 s(-1), measured using a fluorescent phosphate sensor, MDCC-PBP. These two rates imply a movement of approximately 3 bp per ATP hydrolyzed. We demonstrate in several trapping experiments that RecG remains attached to DNA after translocation to the end of the arm of the synthetic DNA junction. ATPase activity continues after translocation is complete. Dissociation of RecG from the product DNA occurs only very slowly, suggesting strong interactions between them. The data support the idea that interactions of the duplex template arm with the protein are the major sites of binding and production of translocation.     

Multiple-binding-site mechanism explains concentration-dependent unbinding rates of DNA-binding proteins

Recent work has demonstrated concentration-dependent unbinding rates of proteins from DNA, using fluorescence visualization of the bacterial nucleoid protein Fis [Graham et al. (2011) (Concentration-dependent exchange accelerates turnover of proteins bound to double-stranded DNA. Nucleic Acids Res., 39:2249)]. The physical origin of this concentration-dependence is unexplained. We use a combination of coarse-grained simulation and theory to demonstrate that this behavior can be explained by taking into account the dimeric nature of the protein, which permits partial dissociation and exchange with other proteins in solution. Concentration-dependent unbinding is generated by this simple model, quantitatively explaining experimental data. This effect is likely to play a major role in determining binding lifetimes of proteins in vivo where there are very high concentrations of solvated molecules.     

Effects of various single-stranded-DNA-binding proteins on reactions promoted by RecA protein

To relate the roles of Escherichia coli SSB in recombination in vivo and in vitro, we have studied the mutant proteins SSB-1 and SSB-113, the variant SSBc produced by chymotryptic cleavage, the partially homologous variant F SSB (encoded by the E. coli sex factor), and the protein encoded by gene 32 of bacteriophage T4. All of these, with the exception of SSB-1, augmented both the initial rate of homologous pairing and strand exchange promoted by RecA protein. From these and related observations, we conclude that SSB stimulates the initial formation of joint molecules by nonspecifically promoting the binding of RecA protein to single-stranded DNA; that SSB plays no role in synapsis of the RecA nucleoprotein filament with duplex DNA; that stimulation of strand exchange by SSB is similarly nonspecific; and that all members of the class of proteins represented by SSB, F SSB, and gene 32 protein may play equivalent roles in making single-stranded DNA more accessible to RecA protein.     

Forces and energetics of hapten-antibody dissociation: a biased molecular dynamics simulation study

The unbinding of fluorescein from the single-chain Fv fragment of the 4D5Flu antibody is investigated by biased molecular dynamics with an implicit solvation model. To obtain statistically meaningful results, a large number of unbinding trajectories are calculated; they involve a total simulation time of more than 200 ns. Simulations are carried out with a time-dependent perturbation and in the presence of a constant force. The two techniques, which provide complementary information, induce unbinding by favoring an increase in the distance between the ligand and the antibody. This distance is an appropriate progress variable for the dissociation reaction and permits direct comparison of the unbinding forces in the simulations with data from atomic force microscopy (AFM). The time-dependent perturbation generates unfolding pathways that are close to equilibrium and can be used to reconstruct the mean force; i.e. the derivative of the potential of mean force, along the reaction coordinate. This is supported by an analysis of the overall unbinding profile and the magnitude of the mean force, which are similar to those of the unbinding force (i.e. the external force due to the time-dependent perturbation) averaged over several unbinding events.The multiple simulations show that unbinding proceeds along a rather well-defined pathway for a broad range of effective pulling speeds. Initially, there is a distortion of the protein localized in the C-terminal region followed by the fluorescein exit from the binding site. This occurs in steps that involve breaking of specific electrostatic and van der Waals interactions. It appears that the simulations do not explore the same barriers as those measured in the AFM experiments because of the much higher unfolding speed in the former. The dependence of the force on the logarithm of the loading rate is linear and the slope is higher than in the AFM, in agreement with experiment in other systems, where different slopes were observed for different regimes. Based on the unbinding events, mutations in the 4D5Flu antigen binding site are predicted to result in significant changes in the unbinding force.     

Uncoupling DNA translocation and helicase activity in PcrA: direct evidence for an active mechanism

DNA footprinting and nuclease protection studies of PcrA helicase complexed with a 3'-tailed DNA duplex reveal a contact region that covers a significant region of the substrate both in the presence and absence of a non-hydrolysable analogue of ATP, ADPNP. However, details of the interactions of the enzyme with the duplex region are altered upon binding of nucleotide. By combining this information with that obtained from crystal structures of PcrA complexed with a similar DNA substrate, we have designed mutant proteins that are defective in helicase activity but that leave the ATPase and single-stranded DNA translocation activities intact. These mutants are all located in domains 1B and 2B, which interact with the duplex portion of the DNA substrate. Taken together with the crystal structures, these data support an 'active' mechanism for PcrA that involves two distinct ATP-dependent processes: destabilization of the duplex DNA ahead of the enzyme that is coupled to DNA translocation along the single strand product.     

Effect of mutations in the C-terminal domain of Mu B on DNA binding and interactions with Mu A transposase

Bacteriophage Mu transposition requires two phage-encoded proteins, the transposase, Mu A, and an accessory protein, Mu B. Mu B is an ATP-dependent DNA-binding protein that is required for target capture and target immunity and is an allosteric activator of transpososome function. The recent NMR structure of the C-terminal domain of Mu B (Mu B223-312) revealed that there is a patch of positively charged residues on the solvent-exposed surface. This patch may be responsible for the nonspecific DNA binding activity displayed by the purified Mu B223-312 peptide. We show that mutations of three lysine residues within this patch completely abolish nonspecific DNA binding of the C-terminal peptide (Mu B223- 312). To determine how this DNA binding activity affects transposition we mutated these lysine residues in the full-length protein. The full-length protein carrying all three mutations was deficient in both strand transfer and allosteric activation of transpososome function but retained ATPase activity. Peptide binding studies also revealed that this patch of basic residues within the C-terminal domain of Mu B is within a region of the protein that interacts directly with Mu A. Thus, we conclude that this protein segment contributes to both DNA binding and protein-protein contacts with the Mu transposase.