"Flexible hinge" dynamics in mismatched DNA revealed by fluorescence correlation spectroscopy.

Altered unwinding/bending fluctuations at DNA lesion sites are implicated as plausible mechanisms for damage sensing by DNA-repair proteins. These dynamics are expected to occur on similar timescales as one-dimensional (1D) diffusion of proteins on DNA if effective in stalling these proteins as they scan DNA. We examined the flexibility and dynamics of DNA oligomers containing 3 base pair (bp) mismatched sites specifically recognized in vitro by nucleotide excision repair protein Rad4 (yeast ortholog of mammalian XPC). A previous Forster resonance energy transfer (FRET) study mapped DNA conformational distributions with cytosine analog FRET pair primarily sensitive to DNA twisting/unwinding deformations (Chakraborty et al. Nucleic Acids Res. 46: 1240-1255 (2018)). These studies revealed B-DNA conformations for nonspecific (matched) constructs but significant unwinding for mismatched constructs specifically recognized by Rad4, even in the absence of Rad4. The timescales of these unwinding fluctuations, however, remained elusive. Here, we labeled DNA with Atto550/Atto647N FRET dyes suitable for fluorescence correlation spectroscopy (FCS). With these probes, we detected higher FRET in specific, mismatched DNA compared with matched DNA, reaffirming unwinding/bending deformations in mismatched DNA. FCS unveiled the dynamics of these spontaneous deformations at ~ 300 µs with no fluctuations detected for matched DNA within the ~ 600 ns-10 ms FCS time window. These studies are the first to visualize anomalous unwinding/bending fluctuations in mismatched DNA on timescales that overlap with the < 500 µs "stepping" times of repair proteins on DNA. Such "flexible hinge" dynamics at lesion sites could arrest a diffusing protein to facilitate damage interrogation and recognition.

Two-step interrogation then recognition of DNA binding site by Integration Host Factor: an architectural DNA-bending protein.

The dynamics and mechanism of how site-specific DNA-bending proteins initially interrogate potential binding sites prior to recognition have remained elusive for most systems. Here we present these dynamics for Integration Host factor (IHF), a nucleoid-associated architectural protein, using a µs-resolved T-jump approach. Our studies show two distinct DNA-bending steps during site recognition by IHF. While the faster (∼100 µs) step is unaffected by changes in DNA or protein sequence that alter affinity by >100-fold, the slower (1-10 ms) step is accelerated ∼5-fold when mismatches are introduced at DNA sites that are sharply kinked in the specific complex. The amplitudes of the fast phase increase when the specific complex is destabilized and decrease with increasing [salt], which increases specificity. Taken together, these results indicate that the fast phase is non-specific DNA bending while the slow phase, which responds only to changes in DNA flexibility at the kink sites, is specific DNA kinking during site recognition. Notably, the timescales for the fast phase overlap with one-dimensional diffusion times measured for several proteins on DNA, suggesting that these dynamics reflect partial DNA bending during interrogation of potential binding sites by IHF as it scans DNA.

Global analysis of ion dependence unveils hidden steps in DNA binding and bending by integration host factor.

Proteins that recognize and bind to specific sites on DNA often distort the DNA at these sites. The rates at which these DNA distortions occur are considered to be important in the ability of these proteins to discriminate between specific and nonspecific sites. These rates have proven difficult to measure for most protein-DNA complexes in part because of the difficulty in separating the kinetics of unimolecular conformational rearrangements (DNA bending and kinking) from the kinetics of bimolecular complex association and dissociation. A notable exception is the Integration Host Factor (IHF), a eubacterial architectural protein involved in chromosomal compaction and DNA recombination, which binds with subnanomolar affinity to specific DNA sites and bends them into sharp U-turns. The unimolecular DNA bending kinetics has been resolved using both stopped-flow and laser temperature-jump perturbation. Here we expand our investigation by presenting a global analysis of the ionic strength dependence of specific binding affinity and relaxation kinetics of an IHF-DNA complex. This analysis enables us to obtain each of the underlying elementary rates (DNA bending/unbending and protein-DNA association/dissociation), and their ionic strength dependence, even under conditions where the two processes are coupled. Our analysis indicates interesting differences in the ionic strength dependence of the bi- versus unimolecular steps. At moderate [KCl] (100-500 mM), nearly all the ionic strength dependence to the overall equilibrium binding affinity appears in the bimolecular association/dissociation of an initial, presumably weakly bent, encounter complex, with a slope SK(bi) ≈ 8 describing the loglog-dependence of the equilibrium constant to form this complex on [KCl]. In contrast, the unimolecular equilibrium constant to form the fully wrapped specific complex from the initial complex is nearly independent of [KCl], with SK(uni) < 0.5. This result is counterintuitive because there are at least twice as many ionic protein-DNA contacts in the fully wrapped complex than in the weakly bent intermediate. The following picture emerges from this analysis: in the bimolecular step, the observed [KCl]-dependence is consistent with the number of DNA counterions expected to be released when IHF binds nonspecifically to DNA whereas in the unimolecular reorganization step, the weak [KCl]-dependence suggests that two effects cancel one another. On one hand, formation of additional protein-DNA contacts in the fully wrapped complex releases bound counterions into bulk solution, which is entropically favored by decreasing [salt]. On the other hand, formation of the fully wrapped complex also releases tightly bound water molecules, which is osmotically favored by increasing [salt]. More generally, our global analysis strategy is applicable to other protein-DNA complexes, and opens up the possibility of measuring DNA bending rates in complexes where the unimolecular and bimolecular steps are not easily separable.

A kinetic zipper model with intrachain interactions applied to nucleic acid hairpin folding kinetics.

Single-stranded DNA and RNA hairpin structures with 4-10 nucleotides (nt) in the loop and 5-8 basepairs (bp) in the stem fold on 10-100 µs timescale. In contrast, theoretical estimate of first contact time of two ends of an ideal semiflexible polymer of similar lengths (with persistence length ~2-nt) is 10-100 ns. We propose that this three-orders-of-magnitude difference between these two timescales is a result of roughness in the folding free energy surface arising from intrachain interactions. We present a statistical mechanical model that explicitly includes all misfolded microstates with nonnative Watson-Crick (WC) and non-WC contacts. Rates of interconversion between different microstates are described in terms of two adjustable parameters: the strength of the non-WC interactions (ΔG(nWC)) and the rate at which a basepair is formed adjacent to an existing basepair (k(bp)(+)). The model accurately reproduces the temperature and loop-length dependence of the measured relaxation rates in temperature-jump studies of a 7-bp stem, single-stranded DNA hairpin with 4-20-nt-long poly(dT) loops, with ΔG(nWC) ≈ -2.4 kcal/mol and k(bp)(+) ≥ (1 ns)(-1), in 100 mM NaCl. Thus, our model provides a microscopic interpretation of the slow hairpin folding times as well as an estimate of the strength of intrachain interactions.

Exploring the energy landscape of nucleic acid hairpins using laser temperature-jump and microfluidic mixing.

We have investigated the multidimensionality of the free energy landscape accessible to a nucleic acid hairpin by measuring the relaxation kinetics in response to two very different perturbations of the folding/unfolding equilibrium, either a laser temperature-jump or ion-jump (from rapid mixing with counterions). The two sets of measurements carried out on DNA hairpins (4 or 5 base pairs in the stem and 21-nucleotide polythymine loop), using FRET between end labels or fluorescence of 2-aminopurine in the stem as conformational probes, yield distinctly different relaxation kinetics in the temperature range 10-30 °C and salt range 100-500 mM NaCl, with rapid mixing exhibiting slower relaxation kinetics after an initial collapse of the chain within 8 µs of the counterion mixing time. The discrepancy in the relaxation times increases with increasing temperatures, with rapid mixing times nearly 10-fold slower than T-jump times at 30 °C. These results rule out a simple two-state scenario with the folded and unfolded ensemble separated by a significant free energy barrier, even at temperatures close to the thermal melting temperature T(m). Instead, our results point to the scenario in which the conformational ensemble accessed after counterion condensation and collapse of the chain is distinctly different from the unfolded ensemble accessed with T-jump perturbation. Our data suggest that, even at temperatures in the vicinity of T(m) or higher, the relaxation kinetics obtained from the ion-jump measurements are dominated by the escape from the collapsed state accessed after counterion condensation.

Fast folding of RNA pseudoknots initiated by laser temperature-jump.

RNA pseudoknots are examples of minimal structural motifs in RNA with tertiary interactions that stabilize the structures of many ribozymes. They also play an essential role in a variety of biological functions that are modulated by their structure, stability, and dynamics. Therefore, understanding the global principles that determine the thermodynamics and folding pathways of RNA pseudoknots is an important problem in biology, both for elucidating the folding mechanisms of larger ribozymes as well as addressing issues of possible kinetic control of the biological functions of pseudoknots. We report on the folding/unfolding kinetics of a hairpin-type pseudoknot obtained with microsecond time-resolution in response to a laser temperature-jump perturbation. The kinetics are monitored using UV absorbance as well as fluorescence of extrinsically attached labels as spectroscopic probes of the transiently populated RNA conformations. We measure folding times of 1-6 ms at 37 °C, which are at least 100-fold faster than previous observations of very slow folding pseudoknots that were trapped in misfolded conformations. The measured relaxation times are remarkably similar to predictions of a computational study by Thirumalai and co-workers (Cho, S. S.; Pincus, D.L.; Thirumalai, D. Proc. Natl. Acad. Sci. U. S. A. 2009, 106, 17349-17354). Thus, these studies provide the first observation of a fast-folding pseudoknot and present a benchmark against which computational models can be refined.

Loop dependence of the stability and dynamics of nucleic acid hairpins.

Hairpin loops are critical to the formation of nucleic acid secondary structure, and to their function. Previous studies revealed a steep dependence of single-stranded DNA (ssDNA) hairpin stability with length of the loop (L) as approximately L(8.5 +/- 0.5), in 100 mM NaCl, which was attributed to intraloop stacking interactions. In this article, the loop-size dependence of RNA hairpin stabilities and their folding/unfolding kinetics were monitored with laser temperature-jump spectroscopy. Our results suggest that similar mechanisms stabilize small ssDNA and RNA loops, and show that salt contributes significantly to the dependence of hairpin stability on loop size. In 2.5 mM MgCl2, the stabilities of both ssDNA and RNA hairpins scale as approximately L(4 +/- 0.5), indicating that the intraloop interactions are weaker in the presence of Mg2+. Interestingly, the folding times for ssDNA hairpins (in 100 mM NaCl) and RNA hairpins (in 2.5 mM MgCl2) are similar despite differences in the salt conditions and the stem sequence, and increase similarly with loop size, approximately L(2.2 +/- 0.5) and approximately L(2.6 +/- 0.5), respectively. These results suggest that hairpins with small loops may be specifically stabilized by interactions of the Na+ ions with the loops. The results also reinforce the idea that folding times are dominated by an entropic search for the correct nucleating conformation.

New insights into the transition pathway from nonspecific to specific complex of DNA with Escherichia coli integration host factor.

To elucidate the nature of the transition-state ensemble along the reaction pathway from a nonspecific protein-DNA complex to the specific complex, we have carried out measurements of DNA bending/unbending dynamics on a cognate DNA substrate in complex with integration host factor (IHF), an architectural protein from E. coli that bends its cognate site by approximately 180 degrees . We use a laser temperature jump to perturb the IHF-DNA complex and monitor the relaxation kinetics with time-resolved FRET measurements on DNA substrates end-labeled with a FRET pair. Previously, we showed that spontaneous bending/kinking of DNA, from thermal disruption of base-pairing/-stacking interactions, may be the rate-limiting step in the formation of the specific complex (Kuznetsov, S. V.; Sugimura, S.; Vivas, P.; Crothers, D. M.; Ansari, A. Proc. Natl. Acad. Sci. USA 2006, 103, 18515). Here, we probe the effect of varying [KCl], which affects the stability of the complex, on this rate-limiting step. We find that below approximately 250 mM KCl, the observed relaxation kinetics are from the unimolecular bending/unbending of DNA, and the relaxation rate kr is independent of [KCl]. Above approximately 300 mM KCl, dissociation of the IHF-DNA complex becomes significant, and the observed relaxation process includes contributions from the association/dissociation step, with kr decreasing with increasing [KCl]. The DNA bending step occurs with a positive activation enthalpy, despite the large negative enthalpy change reported for the specific IHF-DNA complex (Holbrook, J. A.; Tsodikov, O. V.; Saecker, R. M.; Record, M. T., Jr. J. Mol. Biol. 2001, 310, 379). Our conclusion from these studies is that in the uphill climb to the transition state, the DNA is kinked, but with no release of ions, as indicated by the salt-independent behavior of k(r) at low [KCl]. Any release of ions in the unimolecular process, together with conformational changes in the protein-DNA complex that facilitate favorable interactions and that contribute to the negative enthalpy change, must occur as the system leaves the transition state, downhill to the final complex.

Microsecond dynamics of protein-DNA interactions: direct observation of the wrapping/unwrapping kinetics of single-stranded DNA around the E. coli SSB tetramer.

The Escherichia coli single-stranded DNA binding protein (SSB) binds selectively to single-stranded (ss) DNA intermediates during DNA replication, recombination and repair. Each subunit of the homo-tetrameric protein contains a potential ssDNA binding site, thus the protein can bind to ssDNA in multiple binding modes, one of which is the (SSB)(65) mode, in which a 65 nucleotide stretch of ssDNA interacts with and wraps around all four subunits of the tetramer. Previous stopped-flow kinetic studies of (SSB)(65) complex formation using the oligodeoxynucleotide, (dT)70, were unable to resolve the initial binding step from the rapid wrapping of ssDNA around the tetramer. Here we report a laser temperature-jump study with resolution in the approximately 500 ns to 4 ms time range, which directly detects these ssDNA wrapping/unwrapping steps. Biphasic time courses are observed with a fast phase that is concentration-independent and which occurs on a time-scale of tens of microseconds, reflecting the wrapping/unwrapping of ssDNA around the SSB tetramer. Analysis of the slower binding phase, in combination with equilibrium binding and stopped-flow kinetic studies, also provides evidence for a previously undetected intermediate along the pathway to forming the (SSB)(65) complex.

Is hairpin formation in single-stranded polynucleotide diffusion-controlled?

An intriguing puzzle in biopolymer science is the observation that single-stranded DNA and RNA oligomers form hairpin structures on time scales of tens of microseconds, considerably slower than the estimated time for loop formation for a semiflexible polymer of similar length. To address the origin of the slow kinetics and to determine whether hairpin dynamics are diffusion-controlled, the effect of solvent viscosity (eta) on hairpin kinetics was investigated using laser temperature-jump techniques. The viscosity was varied by addition of glycerol, which significantly destabilizes hairpins. A previous study on the viscosity dependence of hairpin dynamics, in which all the changes in the measured rates were attributed to a change in solvent viscosity, reported an apparent scaling of relaxation times (tau(r)) on eta as tau(r) approximately eta(0.8). In this study, we demonstrate that if the effect of viscosity on the measured rates is not deconvoluted from the inevitable effect of change in stability, then separation of tau(r) into opening (tau(o)) and closing (tau(c)) times yields erroneous behavior, with different values (and opposite signs) of the apparent scaling exponents, tau(o) approximately eta(-0.4) and tau(c) approximately eta(1.5). Under isostability conditions, obtained by varying the temperature to compensate for the destabilizing effect of glycerol, both tau(o) and tau(c) scale as approximately eta(1.1+/-0.1). Thus, hairpin dynamics are strongly coupled to solvent viscosity, indicating that diffusion of the polynucleotide chain through the solvent is involved in the rate-determining step.

Mapping the transition state for DNA bending by IHF.

How DNA-bending proteins recognize their specific sites on DNA remains elusive, particularly for proteins that use indirect readout, which relies on sequence-dependent variations in DNA flexibility/bendability. The question remains as to whether the protein bends the DNA (protein-induced bending) or, alternatively, "prebent" DNA conformations are thermally accessible, which the protein captures to form the specific complex (conformational capture). To distinguish between these mechanisms requires characterization of reaction intermediates and, in particular, snapshots of the transition state along the recognition pathway. We present such a snapshot, from measurements of DNA bending dynamics in complex with Escherichia coli integration host factor (IHF), an architectural protein that bends specific sites on λ-DNA in a U-turn by creating two sharp kinks in DNA. Fluorescence resonance energy transfer measurements in response to laser temperature-jump perturbation monitor DNA bending. We find that nicks or mismatches that enhance DNA flexibility at the site of the kinks show 3- to 4-fold increase in DNA bending rates that reflect a 4- to 11-fold increase in binding affinities, while sequence modifications away from the kink sites, as well as mutations in IHF designed to destabilize the complex, have negligible effect on DNA bending rates despite >250-fold decrease in binding affinities. These results support the scenario that the bottleneck in the recognition step for IHF is spontaneous kinking of cognate DNA to adopt a partially prebent conformation and point to conformational capture as the underlying mechanism of initial recognition, with additional protein-induced bending occurring after the transition state.

Direct observation of DNA bending/unbending kinetics in complex with DNA-bending protein IHF.

Regulation of gene expression involves formation of specific protein-DNA complexes in which the DNA is often bent or sharply kinked. Kinetics measurements of DNA bending when in complex with the protein are essential for understanding the molecular mechanism that leads to precise recognition of specific DNA-binding sites. Previous kinetics measurements on several DNA-bending proteins used stopped-flow techniques that have limited time resolution of few milliseconds. Here we use a nanosecond laser temperature-jump apparatus to probe, with submillisecond time resolution, the kinetics of bending/unbending of a DNA substrate bound to integration host factor (IHF), an architectural protein from Escherichia coli. The kinetics are monitored with time-resolved FRET, with the DNA substrates end-labeled with a FRET pair. The temperature-jump measurements, in combination with stopped-flow measurements, demonstrate that the binding of IHF to its cognate DNA site involves an intermediate state with straight or, possibly, partially bent DNA. The DNA bending rates range from approximately 2 ms(-1) at approximately 37 degrees C to approximately 40 ms(-1) at approximately 10 degrees C and correspond to an activation energy of approximately 14 +/- 3 kcal/mol. These rates and activation energy are similar to those of a single A:T base pair opening inside duplex DNA. Thus, our results suggest that spontaneous thermal disruption in base-paring, nucleated at an A:T site, may be sufficient to overcome the free energy barrier needed to partially bend/kink DNA before forming a tight complex with IHF.

Spectroscopic studies of structural changes in two beta-sheet-forming peptides show an ensemble of structures that unfold noncooperatively.

A characterization of the conformation and stability of model peptide systems that form beta-sheets in aqueous solutions is considerably important in gaining insights into the mechanism of beta-sheet formation in proteins. We have characterized the conformation and equilibrium folding and unfolding of two 20-residue peptides whose NMR spectra suggest a three-stranded beta-sheet topology in aqueous solution: Betanova [Kortemme, T., Ramirez-Alvarado, M., and Serrano, L. (1998) Science 281, 253-256] and (D)P(D)P with d-Pro-Gly segments at the turns [Schenck, H. L., and Gellman, S. H. (1998) J. Am. Chem. Soc. 120, 4869-4870]. Both circular dichroism (CD) and infrared measurements indicate only 20-26% beta-sheet-like structure at 5 degrees C for Betanova and 42-59% beta-sheet for (D)P(D)P. For both peptides, the CD and infrared spectra change nearly linearly with increasing temperatures (or urea concentrations) and lack a sigmoidal signature characteristic of cooperative unfolding. Fluorescence resonance energy transfer (FRET) measurements between donor and acceptor molecules attached to the two ends confirm that Betanova is largely unstructured even at 10 degrees C; the average end-to-end distance estimated from FRET is closer to that of a random coil than a structured beta-sheet. In (D)P(D)P, the FRET results indicate a more compact structure that remains compact even at high temperatures (approximately 80 degrees C) or high urea concentrations (approximately 8 M). These results indicate that both these peptides access an ensemble of conformations at all temperatures or denaturant concentrations, with no significant free energy barrier separating the "folded" and "unfolded" conformations.