Linkage of EcoRI dissociation from its specific DNA recognition site to water activity, salt concentration, and pH: separating their roles in specific and non-specific binding

Cosolutes Modify Alkaline Phosphatase Catalysis through Osmotic Stress and Crowding Mechanisms

Examining the effects of different cosolutes on in vitro enzyme kinetics yielded glimpses into their potential behavior when functioning in their natural, complex, in vivo milieu. Viewing cosolute in vitro influences on a model enzyme, calf intestinal alkaline phosphatase, as a combination of competitive and uncompetitive behaviors provided quantitative insights into their effects on catalysis. Observed decreases in the apparent specificity constant, K _asp, caused by the presence of polyethylene glycols or betaine in the reaction solution, indicated interference with enzyme-substrate complex formation. This competitive inhibition appeared to be driven by osmotic stress. Dextran 6 K and sucrose strongly impeded the subsequent conversion of the bound substrate into a free product, which was marked by sharp reductions in V _max, uncompetitive inhibition. For the same step, smaller noncarbohydrate cosolutes, triethylene glycol, polyethylene glycol 400, and betaine, also behaved as uncompetitive inhibitors but to a lesser extent. However, polyethylene glycol 8000 and 20,000 were uncompetitive activators, increasing V _max. Polyethylene glycol of molecular weight 1000 displayed intermediate effects between these two groups of noncarbohydrate cosolutes. These results suggested that crowding has a strong influence on free product formation. The combination of competitive and uncompetitive effects and mixed behaviors, caused by the cosolutes on calf intestinal alkaline phosphatase kinetics, was consistent with the trends seen in similar enzyme-cosolute studies. It is proposed that the double-displacement mechanism of alkaline phosphatases, shared by many other enzymes, could be the root of this general observation.

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

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

Mechanisms of Protein Translocation on DNA Are Differentially Responsive to Water Activity

Water plays important but poorly understood roles in the functions of most biomolecules. We are interested in understanding how proteins use diverse search mechanisms to locate specific sites on DNA; here we present a study of the role of closely associated waters in diverse translocation mechanisms. The bacterial DNA adenine methyltransferase, Dam, moves across large segments of DNA using an intersegmental hopping mechanism, relying in part on movement through bulk water. In contrast, other proteins, such as the bacterial restriction endonuclease EcoRI, rely on a sliding mechanism, requiring the protein to stay closely associated with DNA. Here we probed how these two mechanistically distinct proteins respond to well-characterized osmolytes, dimethyl sulfoxide (DMSO), and glycerol. The ability of Dam to move over large segments of DNA is not impacted by either osmolyte, consistent with its minimal reliance on a sliding mechanism. In contrast, EcoRI endonuclease translocation is significantly enhanced by DMSO and inhibited by glycerol, providing further corroboration that these proteins rely on distinct translocation mechanisms. The well-established similar effects of these osmolytes on bulk water, and their differential effects on macromolecule-associated waters, support our results and provide further evidence of the importance of water in interactions between macromolecules and their ligands.

Effect of (2)H and (18)O water isotopes in kinesin-1 gliding assay

We show for the first time the effects of heavy-hydrogen water (²H₂O) and heavy-oxygen water (H₂¹⁸O) on the gliding speed of microtubules on kinesin-1 coated surfaces. Increased fractions of isotopic waters used in the motility solution decreased the gliding speed of microtubules by a maximum of 21% for heavy-hydrogen and 5% for heavy-oxygen water. We also show that gliding microtubule speed returns to its original speed after being treated with heavy-hydrogen water. We discuss possible interpretations of these results and the importance for future studies of water effects on kinesin and microtubules. We also discuss the implication for using heavy waters in biomolecular devices incorporating molecular motors.

Using single-turnover kinetics with osmotic stress to characterize the EcoRV cleavage reaction

Type II restriction endonucleases require metal ions to specifically cleave DNA at canonical sites. Despite the wealth of structural and biochemical information, the number of Mg(2+) ions used for cleavage by EcoRV, in particular, at physiological divalent ion concentrations has not been established. In this work, we employ a single-turnover technique that uses osmotic stress to probe reaction kinetics between an initial specific EcoRV-DNA complex formed in the absence of Mg(2+) and the final cleavage step. With osmotic stress, complex dissociation before cleavage is minimized and the reaction rates are slowed to a convenient time scale of minutes to hours. We find that cleavage occurs by a two-step mechanism that can be characterized by two rate constants. The dependence of these rate constants on Mg(2+) concentration and osmotic pressure gives the number of Mg(2+) ions and water molecules coupled to each kinetic step of the EcoRV cleavage reaction. Each kinetic step is coupled to the binding 1.5-2.5 Mg(2+) ions, the uptake of ∼30 water molecules, and the cleavage of a DNA single strand. We suggest that each kinetic step reflects an independent, rate-limiting conformational change of each monomer of the dimeric enzyme that allows Mg(2+) ion binding. This modified single-turnover protocol has general applicability for metalloenzymes.

DNA concentration-dependent dissociation of EcoRI: direct transfer or reaction during hopping

Direct transfer of proteins between DNA helices is a recognized important feature of the recognition site search process. Direct transfer is characterized by a dissociation rate that depends on total DNA concentration. This is taken as evidence for the formation of an intermediate DNA-protein-DNA ternary complex. We find that the dissociation rate of EcoRI-DNA-specific complexes at 80 mM NaCl depends on the concentration of competitor oligonucleotide suggesting that direct transfer contributes to EcoRI dissociation. This dependence on competitor DNA concentration is not seen at 180 mM salt. A careful examination of the salt concentration dependence of the dissociation rate, however, shows that the predictions for the formation of a ternary complex are not observed experimentally. The findings can be rationalized by considering that just after dissociating from a DNA fragment the protein remains in close proximity to that fragment, can reassociate with it, and diffuse back to the recognition site rather than bind to an oligonucleotide in solution, a hopping excursion. The probability that a protein will bind to an oligonucleotide during a hop can be approximately calculated and shown to explain the data. A dependence of the dissociation rate of a DNA-protein complex on competitor DNA concentration does not necessarily mean direct transfer.

A model of sequence-dependent protein diffusion along DNA

We introduce a probabilistic model for protein sliding motion along DNA during the search of a target sequence. The model accounts for possible effects due to sequence-dependent interaction between the nonspecific DNA and the protein. Hydrogen bonds formed at the target site are used as the main sequence-dependent interaction between protein and DNA. The resulting dynamical properties and the possibility of an experimental verification are discussed in details. We show that, while at large times the process reaches a linear diffusion regime, it initially displays a sub-diffusive behavior. The sub-diffusive regime can last sufficiently long to be of biological interest.

Nanopore translocation dynamics of a single DNA-bound protein

We study the translocation dynamics of a single protein molecule attached to a double-stranded DNA that is threaded through a solid-state nanopore by optical tweezers and an electric field (nanopore force spectroscopy). We find distinct asymmetric and retarded force signals that depend on the protein charge, the DNA elasticity and its counterionic screening in the buffer. A theoretical model where an isolated charge on an elastic, polyelectrolyte strand is experiencing an anharmonic nanopore potential was developed. Its results compare very well with the measured force curves and explain the experimental findings that the force depends linearly on the applied electric field and exhibits a small hysteresis during back and forth translocation cycles. Moreover, the translocation dynamics reflects the stochastic nature of the thermally activated hopping between two adjacent states in the nanopore that can be adequately described by Kramers rate theory.

Solution parameters modulating DNA binding specificity of the restriction endonuclease EcoRV

The DNA binding stringency of restriction endonucleases is crucial for their proper function. The X-ray structures of the specific and non-cognate complexes of the restriction nuclease EcoRV are considerably different suggesting significant differences in the hydration and binding free energies. Nonetheless, the majority of studies performed at pH 7.5, optimal for enzymatic activity, have found a < 10-fold difference between EcoRV binding constants to the specific and nonspecific sequences in the absence of divalent ions. We used a recently developed self-cleavage assay to measure EcoRV-DNA competitive binding and to evaluate the influence of water activity, pH and salt concentration on the binding stringency of the enzyme in the absence of divalent ions. We find the enzyme can readily distinguish specific and nonspecific sequences. The relative specific-nonspecific binding constant increases strongly with increasing neutral solute concentration and with decreasing pH. The difference in number of associated waters between specific and nonspecific DNA-EcoRV complexes is consistent with the differences in the crystal structures. Despite the large pH dependence of the sequence specificity, the osmotic pressure dependence indicates little change in structure with pH. The large osmotic pressure dependence means that measurement of protein-DNA specificity in dilute solution cannot be directly applied to binding in the crowded environment of the cell. In addition to divalent ions, water activity and pH are key parameters that strongly modulate binding specificity of EcoRV.

A high throughput molecular force assay for protein-DNA interactions

An accurate and genome-wide characterization of protein-DNA interactions such as transcription factor binding is of utmost importance for modern biology. Powerful screening methods emerged. But the vast majority of these techniques depend on special labels or markers against the ligand of interest and moreover most of them are not suitable for detecting low-affinity binders. In this article a molecular force assay is described based on measuring comparative unbinding forces of biomolecules for the detection of protein-DNA interactions. The measurement of binding or unbinding forces has several unique advantages in biological applications since the interaction between certain molecules and not the mere presence of one of them is detected. No label or marker against the protein is needed and only specifically bound ligands are detected. In addition the force-based assay permits the detection of ligands over a broad range of affinities in a crowded and opaque ambient environment. We demonstrate that the molecular force assay allows highly sensitive and fast detection of protein-DNA interactions. As a proof of principle, binding of the protein EcoRI to its DNA recognition sequence is measured and the corresponding dissociation constant in the sub-nanomolar range is determined. Furthermore, we introduce a new, simplified setup employing FRET pairs on the molecular level and standard epi-fluorescence for readout. Due to these advancements we can now demonstrate that a feature size of a few microns is sufficient for the measurement process. This will open a new paradigm in high-throughput screening with all the advantages of force-based ligand detection.

Osmotic stress regulates the strength and kinetics of sugar binding to the maltoporin channel

We study the effect of osmotic stress, exerted by salts, on carbohydrate binding to the sugar-specific bacterial channel maltoporin. When the channel is reconstituted into planar lipid bilayers, single events of its occlusion by sugar are seen as transient interruptions in the flow of small ions. We find that, for most salts, changes in the free energy of maltoporin-sugar binding vary linearly with solution osmotic pressure. Such a change in binding with solution osmolarity indicates that for each salt a constant number of salt-excluding water molecules is released upon sugar-maltoporin association at all salt concentrations. We find that larger numbers of water molecules are released upon binding of the cyclic carbohydrate β-cyclodextrin (CD) than upon binding of the corresponding linear homologue maltoheptaose (m7). Remarkably, the extent to which salts affect the binding constants and rates depends sensitively on the type of salt; dehydration in solutions of different anions corresponds to the Hofmeister series. In sodium sulfate solutions, CD and m7 respectively release about 120 and 35 salt-excluding water molecules; in sodium chloride solutions, 35 and 15 waters. No water release is observed with sodium bromide. Finally, by adding adamantane, known to form an inclusion complex with CD, we can infer that CD not only dehydrates but also undergoes a conformational change upon binding to the channel. As a practical outcome, our results also demonstrate how osmotic stress can improve single-molecule detection of different solutes using protein-based nanopores.

Molecular diagnostics for personal medicine using a nanopore

Semiconductor nanotechnology has created the ultimate analytical tool: a nanopore with single molecule sensitivity. This tool offers the intriguing possibility of high-throughput, low cost sequencing of DNA with the absolute minimum of material and preprocessing. The exquisite single molecule sensitivity obviates the need for costly and error-prone procedures like polymerase chain reaction amplification. Instead, nanopore sequencing relies on the electric signal that develops when a DNA molecule translocates through a pore in a membrane. If each base pair has a characteristic electrical signature, then ostensibly a pore could be used to analyze the sequence by reporting all of the signatures in a single read without resorting to multiple DNA copies. The potential for a long read length combined with high translocation velocity should make resequencing inexpensive and allow for haplotyping and methylation profiling in a chromosome.

Stabilizing labile DNA-protein complexes in polyacrylamide gels

The electrophoretic mobility-shift assay (EMSA) is one of the most popular tools in molecular biology for measuring DNA-protein interactions. EMSA, as standardly practiced today, works well for complexes with association binding constants K(a)>10(9) M(-1) under normal conditions of salt and pH. Many DNA-protein complexes are not stable enough so that they dissociate while moving through the gel matrix giving smeared bands that are difficult to quantitate reliably. In this work we demonstrate that the addition of the osmolyte triethylene glycol to polyacrylamide gels dramatically stabilizes labile restriction endonuclease EcoRI complexes with nonspecific DNA sequences enabling quantitation of binding using EMSA. The significant improvement of the technique resulting from the addition of osmolytes to the gel matrix greatly extends the range of binding constants of protein-DNA complexes that can be investigated using this widely used assay. Extension of this approach to other techniques used for separating bound and free components such as gel chromatography and CE is straightforward.

Diffusion of the restriction nuclease EcoRI along DNA

Many specific sequence DNA binding proteins locate their target sequence by first binding to DNA nonspecifically, then by linearly diffusing or hopping along DNA until either the protein dissociates from the DNA or it finds the recognition sequence. We have devised a method for measuring one-dimensional diffusion along DNA based on the ratio of the dissociation rate of protein from DNA fragments containing one specific binding site to the dissociation rate from DNA fragments containing two specific binding sites. Our extensive measurements of dissociation rates and specific-nonspecific relative binding constants of the restriction nuclease EcoRI enable us to determine the diffusion rate of nonspecifically bound protein along the DNA. By varying the distance between the two binding sites, we confirm a linear diffusion mechanism. The sliding rate is relatively insensitive to salt concentration and osmotic pressure, indicating that the protein moves smoothly along the DNA probably following the helical phosphate-sugar backbone of DNA. We calculate a diffusion coefficient for EcoRI of 3 x 10(4) bp(2) s(-)(1) EcoRI is able to diffuse approximately 150 bp, on average, along the DNA in 1 s. This diffusion rate is about 2000-fold slower than the diffusion of free protein in solution. A factor of 40-50 can be accounted for by rotational friction resulting from following the helical path of the DNA backbone. Two possibilities could account for the remaining activation energy: salt bridges between the DNA and the protein are transiently broken, or the water structure at the protein-DNA interface is disrupted as the two surfaces move past each other.

The dynamic side of the Hofmeister effect: a single-molecule nanopore study of specific complex formation

Beyond measurements of equilibria: An alpha-hemolysin nanopore is used as a single-molecule sensor to follow the effects of different salts on the complexation reaction of gamma-cyclodextrin and adamantane carboxylate (see picture). The kinetics underlying the dynamic equilibrium are studied to reveal qualitatively different dynamic actions of various cosolute salts.

Analyzing the forces binding a restriction endonuclease to DNA using a synthetic nanopore

Restriction endonucleases are used prevalently in recombinant DNA technology because they bind so stably to a specific target sequence and, in the presence of cofactors, cleave double-helical DNA specifically at a target sequence at a high rate. Using synthetic nanopores along with molecular dynamics (MD), we have analyzed with atomic resolution how a prototypical restriction endonuclease, EcoRI, binds to the DNA target sequence--GAATTC--in the absence of a Mg(2+) ion cofactor. We have previously shown that there is a voltage threshold for permeation of DNA bound to restriction enzymes through a nanopore that is associated with a nanonewton force required to rupture the complex. By introducing mutations in the DNA, we now show that this threshold depends on the recognition sequence and scales linearly with the dissociation energy, independent of the pore geometry. To predict the effect of mutation in a base pair on the free energy of dissociation, MD is used to qualitatively rank the stability of bonds in the EcoRI-DNA complex. We find that the second base in the target sequence exhibits the strongest binding to the protein, followed by the third and first bases, with even the flanking sequence affecting the binding, corroborating our experiments.

Energy-based prediction of amino acid-nucleotide base recognition

Despite decades of investigations, it is not yet clear whether there are rules dictating the specificity of the interaction between amino acids and nucleotide bases. This issue was addressed by determining, in a dataset consisting of 100 high-resolution protein-DNA structures, the frequency and energy of interaction between each amino acid and base, and the energetics of water-mediated interactions. The analysis was carried out using HINT, a non-Newtonian force field encoding both enthalpic and entropic contributions, and Rank, a geometry-based tool for evaluating hydrogen bond interactions. A frequency- and energy-based preferential interaction of Arg and Lys with G, Asp and Glu with C, and Asn and Gln with A was found. Not only favorable, but also unfavorable contacts were found to be conserved. Water-mediated interactions strongly increase the probability of Thr-A, Lys-A, and Lys-C contacts. The frequency, interaction energy, and water enhancement factors associated with each amino acid-base pair were used to predict the base triplet recognized by the helix motif in 45 zinc fingers, which represents an ideal case study for the analysis of one-to-one amino acid-base pair contacts. The model correctly predicted 70.4% of 135 amino acid-base pairs, and, by weighting the energetic relevance of each amino acid-base pair to the overall recognition energy, it yielded a prediction rate of 89.7%.

Enthalpies of DNA melting in the presence of osmolytes

The melting of DNA in the presence of osmolytes has been studied with the intention of obtaining information about how base pair stability is affected by changes in solution conditions. In previous investigations, the melting enthalpies were assumed to be constant as osmolalities change, but no systematic evaluation of whether this condition is true has been offered. This paper presents calorimetric data on the melting of two synthetic DNA samples in the presence of a number of common osmolytes. Poly(dAdT)*poly(dTdA) and poly(dGdC)*poly(dCdG) melting have been examined by differential scanning calorimetry in solutions containing ethylene glycol, glycerol, sucrose, urea, betaine, PEG 200 and PEG 1450 at increasing osmolalities. The results show small, but significant changes in the enthalpy of melting of the two polynucleotides that are different, depending on the structure of the cosolvent. The polyols, ethylene glycol, glycerol, PEG 200 and also urea all show decreases in melting enthalpy, while betaine and sucrose display increases with increasing concentration of cosolvent. The large stabilizing PEG 1450 shows no change within the experimental errors. Using concepts relating to preferential interactions of the cosolvents with the DNA base pairs, it is possible to interpret some of the observed changes in the thermodynamic properties of melting. The results indicate that there is strong entropy-enthalpy compensation upon melting base pairs, but entropy increases dominate to cause the decreases in stability with increased cosolvent concentration. Excess hydration parameters are evaluated and their magnitudes discussed in terms of changes in cosolvent interactions with the DNA base pairs.

The EcoRI-DNA complex as a model for investigating protein-DNA interactions by atomic force microscopy

Atomic force microscopy (AFM) is a technique widely used to image protein-DNA complexes, and its application has now been extended to the measurements of protein-DNA binding constants and specificities. However, the spreading of the protein-DNA complexes on a flat substrate, generally mica, is required prior to AFM imaging. The influence of the surface on protein-DNA interactions is therefore an issue which needs to be addressed. For that purpose, the extensively studied EcoRI-DNA complex was investigated with the aim of providing quantitative information about the surface influence. The equilibrium binding constant of the complex was determined by AFM at both low and high ionic strengths and compared to electrophoretic mobility shift assay measurements (EMSA). In addition, the effect of the DNA length on dissociation of the protein from its specific site was analyzed. It turned out that the AFM measurements are similar to those obtained by EMSA at high ionic strengths. We then advance the idea that this effect is due to the high counterion concentration near the highly negatively charged mica surface. In addition, a dissociation of the complexes once they are adsorbed onto the surface was observed, which is weakly dependent on the ionic strength contrary to what occurs in solution. Finally, a two-step mechanism, which describes the adsorption of the EcoRI-DNA complexes on the surface, is proposed. This model could also be extended to other protein-DNA complexes.

Energetics of the protein-DNA-water interaction

Background

To understand the energetics of the interaction between protein and DNA we analyzed 39 crystallographically characterized complexes with the HINT (Hydropathic INTeractions) computational model. HINT is an empirical free energy force field based on solvent partitioning of small molecules between water and 1-octanol. Our previous studies on protein-ligand complexes demonstrated that free energy predictions were significantly improved by taking into account the energetic contribution of water molecules that form at least one hydrogen bond with each interacting species.

Results

An initial correlation between the calculated HINT scores and the experimentally determined binding free energies in the protein-DNA system exhibited a relatively poor r² of 0.21 and standard error of ± 1.71 kcal mol^-1. However, the inclusion of 261 waters that bridge protein and DNA improved the HINT score-free energy correlation to an r² of 0.56 and standard error of ± 1.28 kcal mol^-1. Analysis of the water role and energy contributions indicate that 46% of the bridging waters act as linkers between amino acids and nucleotide bases at the protein-DNA interface, while the remaining 54% are largely involved in screening unfavorable electrostatic contacts.

Conclusion

This study quantifies the key energetic role of bridging waters in protein-DNA associations. In addition, the relevant role of hydrophobic interactions and entropy in driving protein-DNA association is indicated by analyses of interaction character showing that, together, the favorable polar and unfavorable polar/hydrophobic-polar interactions (i.e., desolvation) mostly cancel.

Protein phase diagrams II: nonideal behavior of biochemical reactions in the presence of osmolytes

In the age of biochemical systems biology, proteomics, and high throughput methods, the thermodynamic quantification of cytoplasmatic reaction networks comes into reach of the current generation of scientists. What is needed to efficiently extract the relevant information from the raw data is a robust tool for evaluating the number and stoichiometry of all observed reactions while providing a good estimate of the thermodynamic parameters that determine the molecular behavior. The recently developed phase-diagram method, strictly speaking a graphical representation of linkage or Maxwell Relations, offers such capabilities. Here, we extend the phase diagram method to nonideal conditions. For the sake of simplicity, we choose as an example a reaction system involving the protein RNase A, its inhibitor CMP, the osmolyte urea, and water. We investigate this system as a function of the concentrations of inhibitor and osmolyte at different temperatures ranging from 280 K to 340 K. The most interesting finding is that the protein-inhibitor binding equilibrium depends strongly on the urea concentration--by orders-of-magnitude more than expected from urea-protein interaction alone. Moreover, the m-value of ligand binding is strongly concentration-dependent, which is highly unusual. It is concluded that the interaction between small molecules like urea and CMP can significantly contribute to cytoplasmic nonideality. Such a finding is highly significant because of its impact on renal tissue where high concentrations of cosolutes occur regularly.

Differences in hydration coupled to specific and nonspecific competitive binding and to specific DNA Binding of the restriction endonuclease BamHI

Using the osmotic stress technique together with a self-cleavage assay we measure directly differences in sequestered water between specific and nonspecific DNA-BamHI complexes as well as the numbers of water molecules released coupled to specific complex formation. The difference between specific and nonspecific binding free energy of the BamHI scales linearly with solute osmolal concentration for seven neutral solutes used to set water activity. The observed osmotic dependence indicates that the nonspecific DNA-BamHI complex sequesters some 120-150 more water molecules than the specific complex. The weak sensitivity of the difference in number of waters to the solute identity suggests that these waters are sterically inaccessible to solutes. This result is in close agreement with differences in the structures determined by x-ray crystallography. We demonstrate additionally that when the same solutes that were used in competition experiments are used to probe changes accompanying the binding of free BamHI to its specific DNA sequence, the measured number of water molecules released in the binding process is strikingly solute-dependent (with up to 10-fold difference between solutes). This result is expected for reactions resulting in a large change in a surface exposed area.

Sequestered water and binding energy are coupled in complexes of lambda Cro repressor with non-consensus binding sequences

We use the osmotic pressure dependence of dissociation rates and relative binding constants to infer differences in sequestered water among complexes of lambda Cro repressor with varied DNA recognition sequences. For over a 1000-fold change in association constant, the number of water molecules sequestered by non-cognate complexes varies linearly with binding free energy. One extra bound water molecule is coupled with the loss of approximately 150 cal/mol complex in binding free energy. Equivalently, every tenfold decrease in binding constant at constant salt and temperature is associated with eight to nine additional water molecules sequestered in the non-cognate complex. The relative insensitivity of the difference in water molecules to the nature of the osmolyte used to probe the reaction suggests that the water is sterically sequestered. If the previously measured changes in heat capacity for lambda Cro binding to different non-cognate sequences are attributed solely to this change in water, then the heat capacity change per incorporated water is almost the same as the difference between ice and water. The associated changes in enthalpies and entropies, however, indicate that the change in complex structure involves more than a simple incorporation of fixed water molecules that act as adaptors between non-complementary surfaces.

Trapping DNA-protein binding reactions with neutral osmolytes for the analysis by gel mobility shift and self-cleavage assays

We take advantage of our previous observation that neutral osmolytes can strongly slow down the rate of DNA–protein complex dissociation to develop a method that uses osmotic stress to ‘freeze’ mixtures of DNA–protein complexes and prevent further reaction enabling analysis of the products. We apply this approach to the gel mobility shift assay and use it to modify a self-cleavage assay that uses the nuclease activity of the restriction endonucleases to measure sensitively their specific binding to DNA. At sufficiently high concentrations of neutral osmolytes the cleavage reaction can be triggered at only those DNA fragments with initially bound enzyme. The self-cleavage assay allows measurement of binding equilibrium and kinetics directly in solution avoiding the intrinsic problems of gel mobility shift and filter binding assays while providing the same sensitivity level. Here we compare the self-cleavage and gel mobility shift assays applied to the DNA binding of EcoRI and BamHI restriction endonucleases. Initial results indicate that BamHI dissociation from its specific DNA sequence is strongly linked to water activity with the half-life time of the specific complex increasing ∼20-fold from 0 to 1 osmolal betaine.

Differences between EcoRI nonspecific and "star" sequence complexes revealed by osmotic stress

The binding of the restriction endonuclease EcoRI to DNA is exceptionally specific. Even a single basepair change ("star" sequence) from the recognition sequence, GAATTC, decreases the binding free energy of EcoRI to values nearly indistinguishable from nonspecific binding. The difference in the number of waters sequestered by the protein-DNA complexes of the "star" sequences TAATTC and CAATTC and by the specific sequence complex determined from the dependence of binding free energy on water activity is also practically indistinguishable at low osmotic pressures from the 110 water molecules sequestered by nonspecific sequence complexes. Novel measurements of the dissociation rates of noncognate sequence complexes and competition equilibrium show that sequestered water can be removed from "star" sequence complexes by high osmotic pressure, but not from a nonspecific complex. By 5 Osm, the TAATTC "star" sequence complex has lost almost 90 of the approximately 110 waters initially present. It is more difficult to remove water from the CAATTC "star" sequence complex. The sequence dependence of water loss correlates with the known sequence dependence of "star" cleavage activity.

Comparison of the effect of water release on the interaction of the Saccharomyces cerevisiae TATA binding protein (TBP) with "TATA Box" sequences composed of adenosine or inosine

The formation of sequence-specific complexes of TATA binding protein (TBP) with the minor groove of DNA results in the burial of large nonpolar surfaces and the exclusion of water from these interfaces. The release of water is thus expected to provide a significant entropic driving force for formation of the transcription-preinitiated complexes mediated by the binding of TBP to specific sequences. In this article are described equilibrium-binding studies of Saccharomyces cerevisiae TBP to 14 bp oligonucleotides bearing either the tightly bound and efficiently transcribed adenovirus major late promoter (TATAAAAG) or its inosine-substituted derivative (TITIIIIG) as a function of neutral osmolyte concentration. These two DNA sequences present the same pattern of minor groove hydrogen-bond donors and acceptors to the protein. TBP-DNA complex formation was monitored by steady-state fluorescence resonance energy transfer measurements of the oligonucleotides end-labeled with fluorescein (donor) and TAMRA (acceptor). Correct interpretation of the results obtained with the inosine-substituted sequence required careful consideration of the optical properties of the dyes as a function of osmolyte concentration to demonstrate that the relative change in the end-to-end distances for TATAAAAG- and TITIIIIG-bearing oligonucleotides is the same upon TBP binding. Although the affinity of TBP is slightly greater for the adenosine compared with the inosine-substituted TATA sequence in the absence of osmolyte, the end-to-end distances of the bound DNA in complex with TBP, the enthalpic and electrostatic components of binding, are identical within experimental precision. However, approximately 18 additional molecules of water are released upon TBP binding the TATAAAAG as compared with the TITIIIIG sequence resulting in an entropic advantage to the binding of the natural promoter sequence. These results are considered with regard to differences in the flexibility and hydration of the two DNA sequences.

Energetics of the sequence-specific binding of single-stranded DNA by the F factor relaxase domain

Transfer of conjugative plasmids between bacteria requires the activity of relaxases or mobilization proteins. These proteins nick the plasmid in a site- and strand-specific manner prior to transfer of the cut strand from donor to recipient. TraI36, the relaxase domain of TraI from plasmid F factor, binds a single-stranded DNA (ssDNA) oligonucleotide containing an F factor sequence with high affinity and sequence specificity. To better understand the energetics of this interaction, we examined the temperature, salt, and pH dependence of TraI36 recognition. Binding is entropically driven below 25 degrees C and enthalpically driven at higher temperatures. van't Hoff analysis yields an estimated deltaC(P)(0) of binding (-3300 cal x mol(-1) x K(-1)) that is larger and more negative than that observed for most double-stranded DNA (dsDNA)-binding proteins. Based on analyses of circular dichroism data and the crystal structure of the unliganded protein, we attribute the deltaC(P)(0) to both burial of hydrophobic surface area and coupled folding and binding of the protein. The salt dependence of the binding indicates that several ssDNA phosphates are buried in the complex, and the pH dependence of the binding suggests that some of these ssDNA phosphates form ionic interactions with basic residues of the protein. Although data are available for relatively few sequence-specific ssDNA-binding proteins, sufficient differences exist between TraI36 and other proteins to indicate that, like dsDNA-binding proteins, ssDNA-binding proteins use different motifs and combinations of forces to achieve specific recognition.

Investigation of restriction enzyme cofactor requirements: a relationship between metal ion properties and sequence specificity

Restriction enzymes are important model systems for understanding the mechanistic contributions of metal ions to nuclease activity. These systems are unique in that they combine distinct functions which have been shown to depend on metal ions: high-affinity DNA binding, sequence-specific recognition of DNA, and Mg(II)-dependent phosphodiester cleavage. While Ca(II) and Mn(II) are commonly used to promote DNA binding and cleavage, respectively, the metal ion properties that are critical to the support of these functions are not clear. To address this question, we assessed the abilities of a series of metal ions to promote DNA binding, sequence specificity, and cleavage in the representative PvuII endonuclease. Among the metal ions tested [Ca(II), Sr(II), Ba(II), Eu(III), Tb(III), Cd(II), Mn(II), Co(II), and Zn(II)], only Mn(II) and Co(II) were similar enough to Mg(II) to support detectable cleavage activity. Interestingly, cofactor requirements for the support of DNA binding are much more permissive; the survey of DNA binding cofactors indicated that Cd(II) and the heavier and larger alkaline earth metal ions Sr(II) and Ba(II) were effective cofactors, stimulating DNA binding affinity 20-200-fold. Impressively, the trivalent lanthanides Tb(III) and Eu(III) promoted DNA binding as efficiently as Ca(II), corresponding to an increase in affinity over 1000-fold higher than that observed under metal-free conditions. The trend for DNA binding affinity supported by these ions suggests that ionic radius and charge are not critical to the promotion of DNA binding. To examine the role of metal ions in sequence discrimination, we determined specificity factors [K(a)(specific)/K(a)(nonspecific)] in the presence of Cd(II), Ba(II), and Tb(III). Most interestingly, all of these ions compromised sequence specificity to some degree compared to Ca(II), by either increased affinity for a noncognate sequence, decreased affinity for the cognate sequence, or both. These results suggest that while amino acid-base contacts are important for specificity, the properties of metal ion cofactors at the catalytic site are also critical for sequence discrimination. This insight is invaluable to our efforts to understand and subsequently design sequence-specific nucleases.

Volumetric properties of proteins

Structural and thermodynamic characterizations of a variety of intra- and intermolecular interactions stabilizing/destabilizing protein systems represent a major part of multidisciplinary efforts aimed at solving the problems of protein folding and binding. To this end, volumetric techniques have been successfully used to gain insights into protein hydration and intraglobular packing. Despite the fact that the use of volumetric measurements in protein-related studies dates back to the 1950s, such measurements still represent a relatively untapped yet potentially informative means for tackling the problems of protein folding and binding. This notion has been further emphasized by recent advances in the development of highly sensitive volumetric instrumentation that has led to intensifying volumetric investigations of protein systems. This paper reviews the volumetric properties of proteins and their low-molecular-weight analogs, in particular, discussing the recent progress in the use of volumetric data for studying conformational transitions of proteins as well as protein-ligand, protein-protein, and protein-nucleic acid interactions.

A protein-DNA binding mechanism proceeds through multi-state or two-state parallel pathways

The DNA-binding mechanism of the dimeric C-terminal domain of the papillomavirus E2 protein with its specific DNA target was investigated and shown to proceed through two parallel pathways. A sequential multi-step reaction is initiated by the diffusion-controlled formation of an encounter complex, with no evidence of base sequence discrimination capacity. Following a substantial conformational rearrangement of the protein, a solvent exclusion step leading to the formation of a final protein-DNA complex was identified. This last step involves the largest burial of surface area from the interface and involves the consolidation of the direct readout of the DNA bases. Double-jump stopped-flow experiments allowed us to characterize the sequence of events and demonstrated that a fast-formed consolidated complex can take place through a parallel route. We present the simplest model for the overall mechanism with a description of all the intermediate species in energetic terms.

The thermodynamics of protein-protein recognition as characterized by a combination of volumetric and calorimetric techniques: the binding of turkey ovomucoid third domain to alpha-chymotrypsin

We have used ultrasonic velocimetry, high-precision densimetry, and fluorescence spectroscopy, in conjunction with isothermal titration and differential scanning calorimetry, to characterize the binding of turkey ovomucoid third domain (OMTKY3) to alpha-chymotrypsin. We report the changes in volume and adiabatic compressibility that accompany the association of these proteins at 25 degrees C and pH 4.5. In addition, we report the changes in free energy, enthalpy, entropy, and heat capacity upon the binding of OMTKY3 to alpha-chymotrypsin over a temperature range of 20-40 degrees C. Our volume and compressibility data, in conjunction with X-ray crytsallographic data on the OMTKY3-alpha-chymotrypsin complex, suggest that 454(+/-22) water molecules are released to the bulk state upon the binding of OMTKY3 to alpha-chymotrypsin. Furthermore, these volumetric data suggest that the intrinsic compressibility of the two proteins decreases by 7%. At each temperature studied, OMTKY3 association with alpha-chymotrypsin is entropy driven with a large, unfavorable enthalpy contribution. The observed entropy of the binding reflects interplay between two very large favorable and unfavorable terms. The favorable term reflects an increase in the hydrational entropy resulting from release to the bulk of 454 water molecules. The unfavorable term is related to a decrease in the configurational entropy and, consequently, a decrease in the conformational dynamics of the two proteins. In general, we discuss the relationship between macroscopic and microscopic properties, in particular, identifying and quantifying the role of hydration in determining the thermodynamics of protein recognition as reflected in volumetric and calorimetric parameters.