Monte Carlo applications to thermal and chemical denaturation experiments of nucleic acids and proteins

The Role of XPB/Ssl2 dsDNA Translocase Processivity in Transcription Start-site Scanning

The general transcription factor TFIIH contains three ATP-dependent catalytic activities. TFIIH functions in nucleotide excision repair primarily as a DNA helicase and in Pol II transcription initiation as a dsDNA translocase and protein kinase. During initiation, the XPB/Ssl2 subunit of TFIIH couples ATP hydrolysis to dsDNA translocation facilitating promoter opening and the kinase module phosphorylates Pol II to facilitate the transition to elongation. These functions are conserved between metazoans and yeast; however, yeast TFIIH also drives transcription start-site scanning in which Pol II scans downstream DNA to locate productive start-sites. The ten-subunit holo-TFIIH from S. cerevisiae has a processive dsDNA translocase activity required for scanning and a structural role in scanning has been ascribed to the three-subunit TFIIH kinase module. Here, we assess the dsDNA translocase activity of ten-subunit holo- and core-TFIIH complexes (i.e. seven subunits, lacking the kinase module) from both S. cerevisiae and H. sapiens. We find that neither holo nor core human TFIIH exhibit processive translocation, consistent with the lack of start-site scanning in humans. Furthermore, in contrast to holo-TFIIH, the S. cerevisiae core-TFIIH also lacks processive translocation and its dsDNA-stimulated ATPase activity was reduced ~5-fold to a level comparable to the human complexes, potentially explaining the reported upstream shift in start-site observed in vitro in the absence of the S. cerevisiae kinase module. These results suggest that neither human nor S. cerevisiae core-TFIIH can translocate efficiently, and that the S. cerevisiae kinase module functions as a processivity factor to allow for robust transcription start-site scanning.

Escherichia coli DnaK Allosterically Modulates ClpB between High- and Low-Peptide Affinity States

ClpB and DnaKJE provide protection to Escherichia coli cells during extreme environmental stress. Together, this co-chaperone system can resolve protein aggregates, restoring misfolded proteins to their native form and function in solubilizing damaged proteins for removal by the cell's proteolytic systems. DnaK is the component of the KJE system that directly interacts with ClpB. There are many hypotheses for how DnaK affects ClpB-catalyzed disaggregation, each with some experimental support. Here, we build on our recent work characterizing the molecular mechanism of ClpB-catalyzed polypeptide translocation by developing a stopped-flow FRET assay that allows us to detect ClpB's movement on model polypeptide substrates in the absence or presence of DnaK. We find that DnaK induces ClpB to dissociate from the polypeptide substrate. We propose that DnaK acts as a peptide release factor, binding ClpB and causing the ClpB conformation to change to a low-peptide affinity state. Such a role for DnaK would allow ClpB to rebind to another portion of an aggregate and continue nonprocessive translocation to disrupt the aggregate.

Molecular motor translocation kinetics: Application of Monte Carlo computer simulations to determine microscopic kinetic parameters

Methods for studying the translocation of motor proteins along a filament (e.g., nucleic acid and polypeptide) typically monitor the total production of ADP, the arrival/departure of the motor protein at/from a particular location (often one end of the filament), or the dissociation of the motor protein from the filament. The associated kinetic time courses are often analyzed using a simple sequential uniform n-step mechanism to estimate the macroscopic kinetic parameters (e.g., translocation rate and processivity) and the microscopic kinetic parameters (e.g., kinetic step-size and the rate constant for the rate-limiting step). These sequential uniform n-step mechanisms assume repetition of uniform and irreversible rate-limiting steps of forward motion along the filament. In order to determine how the presence of non-uniform motion (e.g., backward motion, random pauses, or jumping) affects the estimates of parameters obtained from such analyses, we evaluated computer simulated translocation time courses containing non-uniform motion using a simple sequential uniform n-step model. By comparing the kinetic parameters estimated from the analysis of the data generated by these simulations with the input parameters of the simulations, we were able to determine which of the kinetic parameters were likely to be over/under estimated due to non-uniform motion of the motor protein.

Low processivity for DNA translocation by the ISWI molecular motor

The motor protein ISWI (Imitation SWItch) is the conserved catalytic ATPase domain of the ISWI family of chromatin remodelers. Members of the ISWI family are involved in regulating the structure of cellular chromatin during times of transcription, translation, and repair. Current models for the nucleosome repositioning activity of ISWI and other chromatin remodelers require the translocation of the remodeling protein along double-stranded DNA through an ATP-dependent mechanism. Here we report results from spectrofluorometric stopped-flow experiments which demonstrate that ISWI displays very low processivity for free DNA translocation. By combining these results with those from experiments monitoring the DNA stimulated ATPase activity of ISWI we further demonstrate that the DNA translocation by ISWI is tightly coupled to ATP hydrolysis. The calculated coupling efficiency of 0.067±0.018 ATP/ISWI/bp is seemingly quite low in comparison to similar DNA translocases and we present potential models to account for this. Nevertheless, the tight coupling of ATP hydrolysis to DNA translocation suggests that DNA translocation is not energetically rate limiting for nucleosome repositioning by ISWI.

Translocation of Saccharomyces cerevisiae Pif1 helicase monomers on single-stranded DNA

In Saccharomyces cerevisiae Pif1 participates in a wide variety of DNA metabolic pathways both in the nucleus and in mitochondria. The ability of Pif1 to hydrolyse ATP and catalyse unwinding of duplex nucleic acid is proposed to be at the core of its functions. We recently showed that upon binding to DNA Pif1 dimerizes and we proposed that a dimer of Pif1 might be the species poised to catalysed DNA unwinding. In this work we show that monomers of Pif1 are able to translocate on single-stranded DNA with 5′ to 3′ directionality. We provide evidence that the translocation activity of Pif1 could be used in activities other than unwinding, possibly to displace proteins from ssDNA. Moreover, we show that monomers of Pif1 retain some unwinding activity although a dimer is clearly a better helicase, suggesting that regulation of the oligomeric state of Pif1 could play a role in its functioning as a helicase or a translocase. Finally, although we show that Pif1 can translocate on ssDNA, the translocation profiles suggest the presence on ssDNA of two populations of Pif1, both able to translocate with 5′ to 3′ directionality.

Kinetic mechanism of DNA translocation by the RSC molecular motor

ATP-dependent nucleosome repositioning by chromatin remodeling enzymes requires the translocation of these enzymes along the nucleosomal DNA. Using a fluorescence stopped-flow assay we monitored DNA translocation by a minimal RSC motor and through global analysis of these time courses we have determined that this motor has a macroscopic translocation rate of 2.9 bp/s with a step size of 1.24 bp. From the complementary quantitative analysis of the associated time courses of ATP consumption during DNA translocation we have determined that this motor has an efficiency of 3.0 ATP/bp, which is slightly less that the efficiency observed for several genetically related DNA helicases and which likely results from random pausing by the motor during translocation. Nevertheless, this motor is able to exert enough force during translocation to displace streptavidin from biotinylated DNA. Taken together these results are the necessary first step for quantifying both the role of DNA translocation in nucleosome repositioning by RSC and the efficiency at which RSC couples ATP binding and hydrolysis to nucleosome repositioning.

Asymmetric regulation of bipolar single-stranded DNA translocation by the two motors within Escherichia coli RecBCD helicase

Repair of double-stranded DNA breaks in Escherichia coli is initiated by the RecBCD helicase that possesses two superfamily-1 motors, RecB (3' to 5' translocase) and RecD (5' to 3' translocase), that operate on the complementary DNA strands to unwind duplex DNA. However, it is not known whether the RecB and RecD motors act independently or are functionally coupled. Here we show by directly monitoring ATP-driven single-stranded DNA translocation of RecBCD that the 5' to 3' rate is always faster than the 3' to 5' rate on DNA without a crossover hotspot instigator site and that the translocation rates are coupled asymmetrically. That is, RecB regulates both 3' to 5' and 5' to 3' translocation, whereas RecD only regulates 5' to 3' translocation. We show that the recently identified RecBC secondary translocase activity functions within RecBCD and that this contributes to the coupling. This coupling has implications for how RecBCD activity is regulated after it recognizes a crossover hotspot instigator sequence during DNA unwinding.

Single-stranded DNA translocation of E. coli UvrD monomer is tightly coupled to ATP hydrolysis

Escherichia coli UvrD is an SF1A (superfamily 1 type A) helicase/translocase that functions in several DNA repair pathways. A UvrD monomer is a rapid and processive single-stranded DNA (ssDNA) translocase but is unable to unwind DNA processively in vitro. Based on data at saturating ATP (500 μM), we proposed a nonuniform stepping mechanism in which a UvrD monomer translocates with biased (3' to 5') directionality while hydrolyzing 1 ATP per DNA base translocated, but with a kinetic step size of 4-5 nt/step, suggesting that a pause occurs every 4-5 nt translocated. To further test this mechanism, we examined UvrD translocation over a range of lower ATP concentrations (10-500 μM ATP), using transient kinetic approaches. We find a constant ATP coupling stoichiometry of ∼1 ATP/DNA base translocated even at the lowest ATP concentration examined (10 μM), indicating that ATP hydrolysis is tightly coupled to forward translocation of a UvrD monomer along ssDNA with little slippage or futile ATP hydrolysis during translocation. The translocation kinetic step size remains constant at 4-5 nt/step down to 50 μM ATP but increases to ∼7 nt/step at 10 μM ATP. These results suggest that UvrD pauses more frequently during translocation at low ATP but with little futile ATP hydrolysis.

Dynamic light scattering to study allosteric regulation

The Escherichia coli ClpA protein, like many AAA+ motor proteins, is allosterically regulated by nucleotide binding. We have combined analytical ultracentrifugation approaches with dynamic light scattering (DLS) to examine the self-association properties and the allosteric linkage of assembly to nucleotide binding. Here we present a protocol for the rapid and precise determination of the diffusion coefficient using DLS measurements in a model-independent fashion. When combined with analytical ultracentrifugation experiments, such an approach can yield a more complete understanding of the hydrodynamic and thermodynamic properties of the system.

Allosteric interactions of DNA and nucleotides with S. cerevisiae RSC

RSC (remodel the structure of chromatin) is an essential chromatin remodeler of Saccharomyces cerevisiae that has been shown to have DNA translocase properties. We studied the DNA binding properties of a "trimeric minimal RSC" (RSCt) of the RSC chromatin remodeling complex and the effect of nucleotides on this interaction using fluorescence anisotropy. RSCt binds to 20 bp fluorescein-labeled double-stranded DNA with a K(d) of ∼100 nM. The affinity of RSCt for DNA is reduced in the presence of AMP-PNP and ADP in a concentration-dependent manner with the addition of AMP-PNP having more pronounced effect. These differences in the magnitude at which the binding of ADP and AMP-PNP affects the affinity of DNA binding by RSCt suggest that the physical movement of the enzyme along DNA begins between the binding of ATP and its subsequent hydrolysis. Furthermore, the fact that the highest affinity for DNA binding by RSCt occurs in the absence of bound nucleotide offers a mechanistic explanation for the apparent low processivity of DNA translocation by the enzyme.

Ensemble methods for monitoring enzyme translocation along single stranded nucleic acids

We review transient kinetic methods developed to study the mechanism of translocation of nucleic acid motor proteins. One useful stopped-flow fluorescence method monitors arrival of the translocase at the end of a fluorescently labeled nucleic acid. When conducted under single-round conditions the time courses can be analyzed quantitatively using n-step sequential models to determine the kinetic parameters for translocation (rate, kinetic step size and processivity). The assay and analysis discussed here can be used to study enzyme translocation along a linear lattice such as ssDNA or ssRNA. We outline the methods for experimental design and two approaches, along with their limitations, that can be used to analyze the time courses. Analysis of the full time courses using n-step sequential models always yields an accurate estimate of the translocation rate. An alternative semi-quantitative "time to peak" analysis yields accurate estimates of translocation rates only if the enzyme initiates translocation from a unique site on the nucleic acid. However, if initiation occurs at random sites along the nucleic acid, then the "time to peak" analysis can yield inaccurate estimates of even the rates of translocation depending on the values of other kinetic parameters, especially the rate of dissociation of the translocase. Thus, in those cases analysis of the full time course is needed to obtain accurate estimates of translocation rates.

The macroscopic rate of nucleic acid translocation by hepatitis C virus helicase NS3h is dependent on both sugar and base moieties

The nonstructural protein 3 helicase (NS3h) of hepatitis C virus is a 3'-to-5' superfamily 2 RNA and DNA helicase that is essential for the replication of hepatitis C virus. We have examined the kinetic mechanism of the translocation of NS3h along single-stranded nucleic acid with bases uridylate (rU), deoxyuridylate (dU), and deoxythymidylate (dT), and have found that the macroscopic rate of translocation is dependent on both the base moiety and the sugar moiety of the nucleic acid, with approximate macroscopic translocation rates of 3 nt s(-1) (oligo(dT)), 35 nt s(-1) (oligo(dU)), and 42 nt s(-1) (oligo(rU)), respectively. We found a strong correlation between the macroscopic translocation rates and the binding affinity of the translocating NS3h protein for the respective substrates such that weaker affinity corresponded to faster translocation. The values of K(0.5) for NS3h translocation at a saturating ATP concentration are as follows: 3.3+/-0.4 microM nucleotide (poly(dT)), 27+/-2 microM nucleotide (poly(dU)), and 36+/-2 microM nucleotide (poly(rU)). Furthermore, results of the isothermal titration of NS3h with these oligonucleotides suggest that differences in TDeltaS(0) are the principal source of differences in the affinity of NS3h binding to these substrates. Interestingly, despite the differences in macroscopic translocation rates and binding affinities, the ATP coupling stoichiometries for NS3h translocation were identical for all three substrates (approximately 0.5 ATP molecule consumed per nucleotide translocated). This similar periodicity of ATP consumption implies a similar mechanism for NS3h translocation along RNA and DNA substrates.

Investigation of translocation, DNA unwinding, and protein displacement by NS3h, the helicase domain from the hepatitis C virus helicase

Helicases are motor proteins that are involved in DNA and RNA metabolism, replication, recombination, transcription, and repair. The motors are powered by ATP binding and hydrolysis. Hepatitis C virus encodes a helicase called nonstructural protein (NS3). NS3 possesses protease and helicase activities on its N-terminal and C-terminal domains, respectively. The helicase domain of NS3 is termed NS3h. In vitro, NS3h catalyzes RNA and DNA unwinding in a 3'-5' direction. The directionality of unwinding is thought to arise in part from the enzyme's ability to translocate along DNA, but translocation has not been shown explicitly. We examined the DNA translocase activity of NS3h by using single-stranded oligonucleotide substrates containing a fluorescent probe on the 5' end. NS3h can bind to the ssDNA and in the presence of ATP move toward the 5' end. When the enzyme encounters the fluorescent probe, a fluorescence change is observed that allows translocation to be characterized. Under conditions that favor binding of one NS3h per DNA substrate (100 nM NS3h and 200 nM oligonucleotide), we find that NS3h translocates on ssDNA at a rate of 46 +/- 5 nucleotides/s, and that it can move for 230 +/- 60 nucleotides before dissociating from the DNA. The translocase activity of some helicases is responsible for displacing proteins that are bound to DNA. We studied protein displacement by using a ssDNA oligonucleotide covalently linked to biotin on the 5' end. Upon addition of streptavidin, a "protein block" was placed in the pathway of the helicase. Interestingly, NS3h was unable to displace streptavidin from the end of the oligonucleotide, despite its ability to translocate along the DNA. The DNA unwinding activity of NS3h was examined using a 22 bp duplex DNA substrate under conditions that were identical to those used to study translocation. NS3h exhibited little or no DNA unwinding under single-cycle conditions, supporting the conclusion that NS3h is a relatively poor helicase in its monomeric form, as has been reported. In summary, NS3h translocates on ssDNA as a monomer, but the translocase activity does not correspond to comparable DNA unwinding activity or protein displacement activity under identical conditions.

Kinetic analysis of the genome packaging reaction in bacteriophage lambda

Bacteriophage lambda is a double-stranded DNA virus that infects the Escherichia coli bacterium. lambda genomic DNA is replicated via rolling circle replication, resulting in multiple genomes linked head to tail at the cos site. To insert a single lambda genome into the viral capsid, the lambda terminase enzyme introduces symmetric nicks, 12 bp apart, at the cos site, and then promotes a strand separation reaction, releasing the tail end of the previous genome and leaving a binary complex consisting of lambda terminase bound to the head end of the adjacent genome. Next, the genome is translocated into the interior of the capsid particle, in a process that requires ATP hydrolysis by lambda terminase. Even though DNA packaging has been studied extensively, currently no bulk assays are available that have been optimized to report directly on DNA translocation. Rather, these assays are sensitive to assembly steps reflecting formation of the active, DNA packaging machine. In this work, we have modified the DNase protection assay commonly used to study DNA packaging in several bacteriophage systems, such that it reports directly on the kinetics of the DNA packaging reaction. We have analyzed our DNA packaging data according to an N-step sequential minimal kinetic model and have estimated an overall packaging rate of 119 +/- 8 bp/s, at 4 degrees C and 1 mM ATP. Furthermore, we have measured an apparent step size for the this reaction (m(obs)) of 410 +/- 150 bp. The magnitude of this value indicates that our assay is most likely sensitive to both mechanical steps associated with DNA insertion as well as occasional slow steps that are repeated every >410 bp. These slow steps may be reflective of the pausing events observed in recent single-molecule studies of DNA packaging in bacteriophage lambda [Fuller, D. N., et al. (2007) J. Mol. Biol. 373, 1113-1122]. Finally, we show that either ATP or ADP is required for terminase cutting at cos, to generate the active, DNA packaging complex.

Influence of DNA end structure on the mechanism of initiation of DNA unwinding by the Escherichia coli RecBCD and RecBC helicases

Escherichia coli RecBCD is a bipolar DNA helicase possessing two motor subunits (RecB, a 3'-to-5' translocase, and RecD, a 5'-to-3' translocase) that is involved in the major pathway of recombinational repair. Previous studies indicated that the minimal kinetic mechanism needed to describe the ATP-dependent unwinding of blunt-ended DNA by RecBCD in vitro is a sequential n-step mechanism with two to three additional kinetic steps prior to initiating DNA unwinding. Since RecBCD can "melt out" approximately 6 bp upon binding to the end of a blunt-ended DNA duplex in a Mg(2+)-dependent but ATP-independent reaction, we investigated the effects of noncomplementary single-stranded (ss) DNA tails [3'-(dT)(6) and 5'-(dT)(6) or 5'-(dT)(10)] on the mechanism of RecBCD and RecBC unwinding of duplex DNA using rapid kinetic methods. As with blunt-ended DNA, RecBCD unwinding of DNA possessing 3'-(dT)(6) and 5'-(dT)(6) noncomplementary ssDNA tails is well described by a sequential n-step mechanism with the same unwinding rate (mk(U)=774+/-16 bp s(-1)) and kinetic step size (m=3.3+/-1.3 bp), yet two to three additional kinetic steps are still required prior to initiation of DNA unwinding (k(C)=45+/-2 s(-1)). However, when the noncomplementary 5' ssDNA tail is extended to 10 nt [5'-(dT)(10) and 3'-(dT)(6)], the DNA end structure for which RecBCD displays optimal binding affinity, the additional kinetic steps are no longer needed, although a slightly slower unwinding rate (mk(U)=538+/-24 bp s(-1)) is observed with a similar kinetic step size (m=3.9+/-0.5 bp). The RecBC DNA helicase (without the RecD subunit) does not initiate unwinding efficiently from a blunt DNA end. However, RecBC does initiate well from a DNA end possessing noncomplementary twin 5'-(dT)(6) and 3'-(dT)(6) tails, and unwinding can be described by a simple uniform n-step sequential scheme, without the need for the additional k(C) initiation steps, with a similar kinetic step size (m=4.4+/-1.7 bp) and unwinding rate (mk(obs)=396+/-15 bp s(-1)). These results suggest that the additional kinetic steps with rate constant k(C) required for RecBCD to initiate unwinding of blunt-ended and twin (dT)(6)-tailed DNA reflect processes needed to engage the RecD motor with the 5' ssDNA.

Bacillus stearothermophilus PcrA monomer is a single-stranded DNA translocase but not a processive helicase in vitro

Structural studies of the Bacillus stearothermophilus PcrA protein along with biochemical studies of the single-stranded (ss) DNA translocation activity of PcrA monomers have led to the suggestion that a PcrA monomer possesses processive helicase activity in vitro. Yet definitive studies testing whether the PcrA monomer possesses processive helicase activity have not been performed. Here we show, using single turnover kinetic methods, that monomers of PcrA are able to translocate along ssDNA, in the 3' to 5' direction, rapidly and processively, whereas these same monomers display no detectable helicase activity under the same solution conditions in vitro. The PcrA monomer ssDNA translocation activity, although necessary, is not sufficient for processive helicase activity, and thus the translocase and helicase activities of PcrA are separable. These results also suggest that the helicase activity of PcrA needs to be activated either by self-assembly or through interactions with accessory proteins. This same behavior is displayed by both the Escherichia coli Rep and UvrD monomers. Hence, all three of these SF1 enzymes are ssDNA translocases as monomers but do not display processive helicase activity in vitro unless activated. The fact that the translocase and helicase activities are separable suggests that each activity may be used for different functions in vivo.

Mechanism of ATP-dependent translocation of E.coli UvrD monomers along single-stranded DNA

Escherichia coli UvrD protein is a 3' to 5' SF1 DNA helicase involved in methyl-directed mismatch repair and nucleotide excision repair of DNA. Using stopped-flow methods we have examined the kinetic mechanism of translocation of UvrD monomers along single-stranded DNA (ssDNA) in vitro by monitoring the transient kinetics of arrival of protein at the 5'-end of the ssDNA. Arrival at the 5'-end was monitored by the effect of protein on the fluorescence intensity of fluorophores (Cy3 or fluorescein) attached to the 5'-end of a series of oligodeoxythymidylates varying in length from 16 to 124 nt. We find that UvrD monomers are capable of ATP-dependent translocation along ssDNA with a biased 3' to 5' directionality. Global non-linear least-squares analysis of the full kinetic time-courses in the presence of a protein trap to prevent rebinding of free protein to the DNA using the methods described in the accompanying paper enabled us to obtain quantitative estimates of the kinetic parameters for translocation. We find that UvrD monomers translocate in discrete steps with an average kinetic step-size, m=3.68(+/-0.03) nt step(-1), a translocation rate constant, kt=51.3(+/-0.6) steps s(-1), (macroscopic translocation rate, mkt=189.0(+/-0.7) nt s(-1)), with a processivity corresponding to an average translocation distance of 2400(+/-600) nt before dissociation (10 mM Tris-HCl (pH 8.3), 20 mM NaCl, 20% (v/v) glycerol, 25 degrees C). However, in spite of its ability to translocate rapidly and efficiently along ssDNA, a UvrD monomer is unable to unwind even an 18 bp duplex in vitro. DNA helicase activity in vitro requires a UvrD dimer that unwinds DNA with a similar kinetic step-size of 4-5 bp step(-1), but an approximately threefold slower unwinding rate of 68(+/-9) bp s(-1) under the same solution conditions, indicating that DNA unwinding activity requires more than the ability to simply translocate directionally along ss-DNA.

ATP-dependent translocation of proteins along single-stranded DNA: models and methods of analysis of pre-steady state kinetics

Processive DNA helicases are able to translocate along single-stranded DNA (ssDNA) with biased directionality in a nucleoside triphosphate-dependent reaction, although translocation is not generally sufficient for helicase activity. An understanding of the mechanism of protein translocation along ssDNA requires pre-steady state transient kinetic experiments. Although ensemble experimental approaches have been developed recently for the study of translocation of proteins along DNA, quantitative analysis of the complete time-courses from these experiments, which is needed to obtain quantitative estimates of translocation kinetic parameters (rate constants, processivity, step sizes and ATP coupling) has been lacking. We discuss three ensemble transient kinetic experiments that can be used to study protein translocation along ssDNA, along with the advantages and limitations of each approach. We further describe methods to analyze the complete kinetic time-courses obtained from such experiments performed with a series of ssDNA lengths under "single-round" conditions (i.e. in the absence of re-binding of dissociated protein to DNA). These analysis methods utilize a sequential "n-step" model for protein translocation along ssDNA and enable quantitative determinations of the rate constant, processivity and step size for translocation through global non-linear least-squares fitting of the full time-courses.

Effects of temperature and ATP on the kinetic mechanism and kinetic step-size for E.coli RecBCD helicase-catalyzed DNA unwinding

The kinetic mechanism by which Escherichia coli RecBCD helicase unwinds duplex DNA was studied using a fluorescence stopped-flow method. Single turnover DNA unwinding experiments were performed using a series of fluorescently labeled DNA substrates containing duplex DNA regions ranging from 24 bp to 60 bp. All or no DNA unwinding time courses were obtained by monitoring the changes in fluorescence resonance energy transfer between a Cy3 donor and Cy5 acceptor fluorescent pair placed on opposite sides of a nick in the duplex DNA. From these experiments one can determine the average rates of DNA unwinding as well as a kinetic step-size, defined as the average number of base-pairs unwound between two successive rate-limiting steps repeated during DNA unwinding. In order to probe how the kinetic step-size might relate to a mechanical step-size, we performed single turnover experiments as a function of [ATP] and temperature. The apparent unwinding rate constant, kUapp, decreases with decreasing [ATP], exhibiting a hyperbolic dependence on [ATP] (K1/2=176(+/-30) microM) and a maximum rate of kUapp=204(+/-4) steps s(-1) (mkUapp=709(+/-14) bp s(-1)) (10 mM MgCl2, 30 mM NaCl (pH 7.0), 5% (v/v) glycerol, 25.0 degrees C). kUapp also increases with increasing temperature (10-25 degrees C), with Ea=19(+/-1) kcal mol(-1). However, the average kinetic step-size, m=3.9(+/-0.5) bp step(-1), remains independent of [ATP] and temperature. This indicates that even though the values of the rate constants change, the same elementary kinetic step in the unwinding cycle remains rate-limiting over this range of conditions and this kinetic step remains coupled to ATP binding. The implications of the constancy of the measured kinetic step-size for the mechanism of RecBCD-catalyzed DNA unwinding are discussed.

General methods for analysis of sequential "n-step" kinetic mechanisms: application to single turnover kinetics of helicase-catalyzed DNA unwinding

Helicase-catalyzed DNA unwinding is often studied using "all or none" assays that detect only the final product of fully unwound DNA. Even using these assays, quantitative analysis of DNA unwinding time courses for DNA duplexes of different lengths, L, using "n-step" sequential mechanisms, can reveal information about the number of intermediates in the unwinding reaction and the "kinetic step size", m, defined as the average number of basepairs unwound between two successive rate limiting steps in the unwinding cycle. Simultaneous nonlinear least-squares analysis using "n-step" sequential mechanisms has previously been limited by an inability to float the number of "unwinding steps", n, and m, in the fitting algorithm. Here we discuss the behavior of single turnover DNA unwinding time courses and describe novel methods for nonlinear least-squares analysis that overcome these problems. Analytic expressions for the time courses, f(ss)(t), when obtainable, can be written using gamma and incomplete gamma functions. When analytic expressions are not obtainable, the numerical solution of the inverse Laplace transform can be used to obtain f(ss)(t). Both methods allow n and m to be continuous fitting parameters. These approaches are generally applicable to enzymes that translocate along a lattice or require repetition of a series of steps before product formation.

Fluorescence Stopped-flow Studies of Single Turnover Kinetics of E.coli RecBCD Helicase-catalyzed DNA Unwinding

We have developed and optimized a stopped-flow fluorescence assay for use in studying DNA unwinding catalyzed by Escherichia coli RecBCD helicase. This assay monitors changes in fluorescence resonance energy transfer (FRET) between a pair of fluorescent probes (Cy3 donor and Cy5 acceptor) placed on opposite sides of a nick in duplex DNA. As such, this is an "all-or-none" DNA unwinding assay. Single turnover DNA unwinding experiments were performed using a series of eight fluorescent DNA substrates containing duplex DNA regions ranging from 24 bp to 60 bp. The time-courses obtained by monitoring Cy3 fluorescence display a distinct lag phase that increases with increasing duplex DNA length, reflecting the transient formation of partially unwound DNA intermediates. These Cy3 FRET time-courses are identical with those obtained using a chemical quenched-flow kinetic assay developed previously. The signal from the Cy5 fluorescence probe shows additional effects that appear to specifically monitor the RecD helicase subunit. The continuous nature of this fluorescence assay enabled us to acquire more precise time-courses for many more duplex DNA lengths in a significantly reduced amount of time, compared to quenched-flow methods. Global analysis of the Cy3 and Cy5 FRET time-courses, using an n-step sequential DNA unwinding model, indicates that RecBCD unwinds duplex DNA with an average unwinding rate constant of kU = 200(+/-40) steps s(-1) (mkU = 680(+/-12)bp s(-1)) and an average kinetic step size, m = 3.4 (+/-0.6) bp step(-1) (5 mM ATP, 10 mM MgCl(2), 30 mM NaCl, pH 7.0, 5% (v/v) glycerol, 25.0 degrees C), in excellent agreement with the kinetic parameters determined using quenched-flow techniques. Under these same conditions, the RecBC enzyme unwinds DNA with a very similar rate. These methods will facilitate detailed studies of the mechanisms of DNA unwinding and translocation of the RecBCD and RecBC helicases.

Dynamics of the IRE RNA hairpin loop probed by 2-aminopurine fluorescence and stochastic dynamics simulations

The iron responsive element (IRE) RNA hairpin loop contains six phylogenetically conserved nucleotides, which constitute part of the sequence-specific binding site of the IRE-binding protein. The NMR structure of the loop has been solved, showing that 3 of the 6 nt are poorly constrained. Here, two purine nucleotides in the IRE loop are individually replaced with the fluorescent purine analog 2-aminopurine (2AP). Steady-state and time-resolved fluorescence methods are used to describe the structure and dynamics of 2AP in the IRE loop. The data indicate that 2AP at the position of the adenosine in the loop moves between stacked and unstacked positions, whereas 2AP at the adjacent guanosine is predominantly solvent exposed. Stochastic dynamics simulations are used to provide a physical description of how those nucleotides might move.

Self-association equilibria of Escherichia coli UvrD helicase studied by analytical ultracentrifugation

The Escherichia coli UvrD protein (helicase II) is an SF1 superfamily helicase required for methyl-directed mismatch repair and nucleotide excision repair of DNA. We have characterized quantitatively the self-assembly equilibria of the UvrD protein as a function of [NaCl], [glycerol], and temperature (5-35 degrees C; pH 8.3) using analytical sedimentation velocity and equilibrium techniques, and find that UvrD self-associates into dimeric and tetrameric species over a range of solution conditions (t</=25 degrees C). Increasing [NaCl] from 20mM to 200 mM decreases the dimerization equilibrium constant (L(20)) from 2.33(+/-0.30) microM(-1) to 0.297(+/-0.006) microM(-1) (pH 8.3, 20% (v/v) glycerol, 25 degrees C). The overall tetramerization equilibrium constant (L(40)) is 5.11(+/-0.80) microM(-3) at 20mM NaCl, but decreases so that it is not measurable at 200 mM NaCl. At 500 mM NaCl, only UvrD monomers are detectable. Increasing [glycerol] over the range from 20% to 40% (v/v) decreases both L(20) and L(40). We find no evidence for hexamer formation, although a species consistent in size with an octamer is detected at 35 degrees C. Inclusion of either ADP or ATPgammaS does not affect either L(20) or L(40) significantly, and does not induce the formation of additional assembly states. We also investigated the stoichiometry of UvrD binding to a 3'-(dT)(20)-18 bp DNA substrate by sedimentation equilibrium. At saturating concentrations of UvrD, three UvrD monomers can bind to the DNA substrate, although only two UvrD monomers are required to form a processive helicase complex. When the total DNA substrate concentration is about twofold greater than the total UvrD concentration, the vast majority of the DNA is bound by a single UvrD monomer.

DNA unwinding step-size of E. coli RecBCD helicase determined from single turnover chemical quenched-flow kinetic studies

The mechanism by which Escherichia coli RecBCD DNA helicase unwinds duplex DNA was examined in vitro using pre-steady-state chemical quenched-flow kinetic methods. Single turnover DNA unwinding experiments were performed by addition of ATP to RecBCD that was pre-bound to a series of DNA substrates containing duplex DNA regions ranging from 24 bp to 60 bp. In each case, the time-course for formation of completely unwound DNA displayed a distinct lag phase that increased with duplex length, reflecting the transient formation of partially unwound DNA intermediates during unwinding catalyzed by RecBCD. Quantitative analysis of five independent sets of DNA unwinding time courses indicates that RecBCD unwinds duplex DNA in discrete steps, with an average unwinding "step-size", m=3.9(+/-1.3)bp step(-1), with an average unwinding rate of k(U)=196(+/-77)steps s(-1) (mk(U)=790(+/-23)bps(-1)) at 25.0 degrees C (10mM MgCl(2), 30 mM NaCl (pH 7.0), 5% (v/v) glycerol). However, additional steps, not linked directly to DNA unwinding are also detected. This kinetic DNA unwinding step-size is similar to that determined for the E.coli UvrD helicase, suggesting that these two SF1 superfamily helicases may share similar mechanisms of DNA unwinding.

A functional role for correlated motion in the N-terminal RNA-binding domain of human U1A protein

The N-terminal RNA-binding domain of the human U1A protein (RBD1) undergoes local conformational changes upon binding to its target RNA. Here, the wild-type RBD1 and two mutants are examined with molecular dynamics simulations that are analyzed using the reorientational eigenmode dynamics (RED) formalism. The results reveal changes in the magnitude and extent of coupled intra-domain motions resulting from single amino acid substitutions. Interpretation of the novel RED results and corresponding NMR relaxation data suggests that the loss of collective motions in the mutants could account for their weak RNA-binding.

Experimental and computational studies of the G[UUCG]C RNA tetraloop

In prokaryotic ribosomal RNAs, most UUCG tetraloops are closed by a C-G base-pair. However, this preference is greatly reduced in eukaryotic rRNA species where many UUCG tetraloops are closed by G-C base-pairs. Here, biophysical properties of the C[UUCG]G and G[UUCG]C tetraloops are compared, using experimental and computational methods. Thermal denaturation experiments are used to derive thermodynamic parameters for the wild-type G[UUCG]C tetraloop and variants containing single deoxy substitutions in the loop. A comparison with analogous experiments on the C[UUCG]G motif shows that the two RNA species exhibit similar patterns in response to the substitutions, suggesting that their loop structures are similar. This conclusion is supported by NMR data that suggest that the essential UUCG loop structure is maintained in both tetraloops. However, NMR results show that the G[UUCG]C loop structure is disrupted prior to melting of the stem; this behavior is in contrast to the two-state behavior of the C[UUCG]G molecule. Stochastic dynamics simulations using the GB/SA continuum solvation model, run as a function of temperature, show rare conformational transitions in several G[UUCG]C simulations. These results lead to the conclusion that substitution of a G-C for a C-G closing base-pair increases the intrinsic flexibility of the UUCG loop.