Role of the hydrophobic effect in stability of site-specific protein-DNA complexes

Large contributions of coupled protonation equilibria to the observed enthalpy and heat capacity changes for ssDNA binding to Escherichia coli SSB protein

Many macromolecular interactions, including protein-nucleic acid interactions, are accompanied by a substantial negative heat capacity change, the molecular origins of which have generated substantial interest. We have shown previously that temperature-dependent unstacking of the bases within oligo(dA) upon binding to the Escherichia coli SSB tetramer dominates the binding enthalpy, DeltaH(obs), and accounts for as much as a half of the observed heat capacity change, DeltaC(p). However, there is still a substantial DeltaC(p) associated with SSB binding to ssDNA, such as oligo(dT), that does not undergo substantial base stacking. In an attempt to determine the origins of this heat capacity change, we have examined by isothermal titration calorimetry (ITC) the equilibrium binding of dT(pT)(34) to SSB over a broad pH range (pH 5. 0-10.0) at 0.02 M, 0.2 M NaCl and 1 M NaCl (25 degrees C), and as a function of temperature at pH 8.1. A net protonation of the SSB protein occurs upon dT(pT)(34) binding over this entire pH range, with contributions from at least three sets of protonation sites (pK(a1) = 5.9-6.6, pK(a2) = 8.2-8.4, and pK(a3) = 10.2-10.3) and these protonation equilibria contribute substantially to the observed DeltaH and DeltaC(p) for the SSB-dT(pT)(34) interaction. The contribution of this coupled protonation ( approximately -260 to -320 cal mol(-1) K(-1)) accounts for as much as half of the total DeltaC(p). The values of the "intrinsic" DeltaC(p,0) range from -210 +/- 33 cal mol(-1) degrees K(-1) to -237 +/- 36 cal mol(-1)K(-1), independent of [NaCl]. These results indicate that the coupling of a temperature-dependent protonation equilibria to a macromolecular interaction can result in a large negative DeltaC(p), and this finding needs to be considered in interpretations of the molecular origins of heat capacity changes associated with ligand-macromolecular interactions, as well as protein folding.

Thermostable and site-specific DNA binding of the gene product ORF56 from the Sulfolobus islandicus plasmid pRN1, a putative archael plasmid copy control protein

There is still a lack of information on the specific characteristics of DNA-binding proteins from hyperthermophiles. Here we report on the product of the gene orf56 from plasmid pRN1 of the acidophilic and thermophilic archaeon Sulfolobus islandicus. orf56 has not been characterised yet but low sequence similarily to several eubacterial plasmid-encoded genes suggests that this 6.5 kDa protein is a sequence-specific DNA-binding protein. The DNA-binding properties of ORF56, expressed in Escherichia coli, have been investigated by EMSA experiments and by fluorescence anisotropy measurements. Recombinant ORF56 binds to double-stranded DNA, specifically to an inverted repeat located within the promoter of orf56. Binding to this site could down-regulate transcription of the orf56 gene and also of the overlapping orf904 gene, encoding the putative initiator protein of plasmid replication. By gel filtration and chemical crosslinking we have shown that ORF56 is a dimeric protein. Stoichiometric fluorescence anisotropy titrations further indicate that ORF56 binds as a tetramer to the inverted repeat of its target binding site. CD spectroscopy points to a significant increase in ordered secondary structure of ORF56 upon binding DNA. ORF56 binds without apparent cooperativity to its target DNA with a dissociation constant in the nanomolar range. Quantitative analysis of binding isotherms performed at various salt concentrations and at different temperatures indicates that approximately seven ions are released upon complex formation and that complex formation is accompanied by a change in heat capacity of -6.2 kJ/mol. Furthermore, recombinant ORF56 proved to be highly thermostable and is able to bind DNA up to 85 degrees C.

Structural and thermodynamic strategies for site-specific DNA binding proteins

BACKGROUND

Site-specific protein-DNA complexes vary greatly in structural properties and in the thermodynamic strategy for achieving an appropriate binding free energy. A better understanding of the structural and energetic engineering principles might lead to rational methods for modification or design of such proteins.

RESULTS

A novel analysis of ten site-specific protein-DNA complexes reveals a striking correspondence between the degree of imposed DNA distortion and the thermodynamic parameters of each system. For complexes with relatively undistorted DNA, favorable enthalpy change drives unfavorable entropy change, whereas for complexes with highly distorted DNA, unfavorable DeltaH degrees is driven by favorable DeltaS degrees. We show for the first time that protein-DNA associations have isothermal enthalpy-entropy compensation, distinct from temperature-dependent compensation, so DeltaH degrees and DeltaS degrees do not vary independently. All complexes have favorable DeltaH degrees from direct protein-DNA recognition interactions and favorable DeltaS degrees from water release. Systems that strongly distort the DNA nevertheless have net unfavorable DeltaH degrees as the result of molecular strain, primarily associated with the base pair destacking. These systems have little coupled protein folding and the strained interface suffers less immobilization, so DeltaS degrees is net favorable. By contrast, systems with little DNA distortion have net favorable DeltaH degrees, which must be counterbalanced by net unfavorable DeltaS degrees, derived from loss of vibrational entropy (a result of isothermal enthalpy-entropy compensation) and from coupling between DNA binding and protein folding.

CONCLUSIONS

Isothermal enthalpy-entropy compensation implies that a structurally optimal, unstrained fit is achieved only at the cost of entropically unfavorable immobilization, whereas an enthalpically weaker, strained interface entails smaller entropic penalties.

Macromolecular hydration changes associated with BamHI binding and catalysis

In this report, the effects of osmotic pressure on BamHI cognate binding and catalysis were investigated and compared with a previous study on EcoRI (Robinson, C. R. and Sligar, S. G. (1998) Proc. Natl. Acad. Sci. U.S.A. 95, 2186-2191). Our observation of the dependence of binding and catalytic parameters on osmotic pressure has allowed for the comparison of hydration changes associated with site-specific DNA recognition for both endonucleases. Over a large range of osmotic pressures (pi), the dependence of BamHI on osmotic stress during cognate binding and catalysis was very different from that of the related endonuclease EcoRI. The binding of EcoRI to cognate DNA was dominated by a dehydration of the endonuclease-DNA complex, whereas binding by BamHI to its cognate sequence was accompanied by a solvent release corresponding to some 125 fewer waters. Catalytic analysis at elevated osmotic pressures indicated that both endonucleases had undergone a net hydration of the complex with BamHI displaying a much greater dependence on osmotic stress than EcoRI. Although the enzymes shared core structural motifs, comparisons of high resolution x-ray structures revealed many different secondary structural features of the complexed endonucleases. The large difference in hydration changes by both BamHI and EcoRI could be attributed to these dissimilar secondary structural features, as well as the functional differences of the two endonucleases during site-specific DNA recognition.

Speeding molecular recognition by using the folding funnel: the fly-casting mechanism

Protein folding and binding are kindred processes. Many proteins in the cell are unfolded, so folding and function are coupled. This paper investigates how binding kinetics is influenced by the folding of a protein. We find that a relatively unstructured protein molecule can have a greater capture radius for a specific binding site than the folded state with its restricted conformational freedom. In this scenario of binding, the unfolded state binds weakly at a relatively large distance followed by folding as the protein approaches the binding site: the "fly-casting mechanism." We illustrate this scenario with the hypothetical kinetics of binding a single repressor molecule to a DNA site and find that the binding rate can be significantly enhanced over the rate of binding of a fully folded protein.

A thermodynamic and structural analysis of DNA minor-groove complex formation

As part of an effort to develop a better understanding of the structural and thermodynamic principles of DNA minor groove recognition, we have investigated complexes of three diphenylfuran dications with the d(CGCGAATTCGCG)(2) duplex. The parent compound, furamidine (DB75), has two amidine substituents while DB244 has cyclopentyl amidine substituents and DB226 has 3-pentyl amidines. The structure for the DB244-DNA complex is reported here and is compared to the structure of the DB75 complex. Crystals were not obtained with DB226 but information from the DB75 and DB244 structures as well as previous NMR results on DB226 indicate that all three compounds bind in the minor groove at the AATT site of the duplex. DB244 and DB75 penetrate to the floor of the groove and form hydrogen bonds with T8 on one strand and T20 on the opposite strand while DB226 forms a complex with fewer interactions. Binding studies by surface plasmon resonance (SPR) yield -delta G degrees values in the order DB244>DB75>DB226 that are relatively constant with temperature. The equilibrium binding constants for DB244 are 10-20 times greater than that for DB226. Isothermal titration calorimetric (ITC) experiments indicate that, in contrast to delta G degrees, delta H degrees varies considerably with temperature to yield large negative delta Cp degrees values. The thermodynamic results, analyzed in terms of structures of the DNA complexes, provide an explanation of why DB244 binds more strongly to DNA than DB75, while DB266 binds more weakly. All three compounds have a major contribution to binding from hydrophobic interactions but the hydrophobic term is most favorable for DB244. DB244 also has strong contributions from molecular interactions in its DNA complex and all of these factors combine to give it the largest-delta G degrees for binding. Although the factors that influence the energetics of minor groove interactions are varied and complex, results from the literature coupled with those on the furan derivatives indicate that there are some common characteristics for minor groove recognition by unfused heterocyclic cations that can be used in molecular design.

Large-scale expression and thermodynamic characterization of a glutamate receptor agonist-binding domain

The ionotropic glutamate receptors (GluR) are the primary mediators of excitatory synaptic transmission in the brain. GluR agonist binding has been localized to an extracellular domain whose core is homologous to the bacterial periplasmic binding proteins (PBP). We have established routine, baculovirus-mediated expression of a complete ligand-binding domain construct at the 10-L scale, yielding 10-40 milligrams of purified protein. This construct contains peptides that lie outside the PBP-homologous core and that connect the domain core to the transmembrane domains of the channel and to the N-terminal 'X'-domain. These linker peptides have been implicated in modulating channel physiology. Such extended constructs have proven difficult to express in bacteria, but the protein described here is stable and monomeric. Isothermal titration calorimetry reveals that glutamate binding to the domain involves a substantial heat capacity change and that at physiological temperatures, the reaction is both entropically and enthalpically favorable.

Energetics of drug-DNA interactions

Understanding the thermodynamics of drug binding to DNA is of both practical and fundamental interest. The practical interest lies in the contribution that thermodynamics can make to the rational design process for the development of new DNA targeted drugs. Thermodynamics offer key insights into the molecular forces that drive complex formation that cannot be obtained by structural or computational studies alone. The fundamental interest in these interactions lies in what they can reveal about the general problems of parsing and predicting ligand binding free energies. For these problems, drug-DNA interactions offer several distinct advantages, among them being that the structures of many drug-DNA complexes are known at high resolution and that such structures reveal that in many cases the drug acts as a rigid body, with little conformational change upon binding. Complete thermodynamic profiles (delta G, delta H, delta S, delta Cp) for numerous drug-DNA interactions have been obtained, with the help of high-sensitivity microcalorimetry. The purpose of this article is to offer a perspective on the interpretation of these thermodynamics parameters, and in particular how they might be correlated with known structural features. Obligatory conformational changes in the DNA to accommodate intercalators and the loss of translational and rotational freedom upon complex formation both present unfavorable free energy barriers for binding. Such barriers must be overcome by favorable free energy contributions from the hydrophobic transfer of ligand from solution into the binding site, polyelectrolyte contributions from coupled ion release, and molecular interactions (hydrogen and ionic bonds, van der Waals interactions) that form within the binding site. Theoretical and semiempirical tools that allow estimates of these contributions to be made will be discussed, and their use in dissecting experimental data illustrated. This process, even at the current level of approximation, can shed considerable light on the drug-DNA binding process.

Protein-DNA recognition complexes: conservation of structure and binding energy in the transition state

This paper considers how enzymes that catalyze reactions at specific DNA sites have been engineered to overcome the problem of competitive inhibition by excess nonspecific binding sites on DNA. The formation of a specific protein-DNA recognition complex is discussed from both structural and thermodynamic perspectives, and contrasted with formation of nonspecific complexes. Evidence (from EcoRI and BamHI endonucleases) is presented that a wide variety of perturbations of the DNA substrate alter binding free energy but do not affect the free energy of activation for the chemical step; that is, many energetic factors contribute equally to the recognition complex and the transition-state complex. This implies that the specific recognition complex bears a close resemblance to the transition-state complex, such that very tight binding to the recognition site on the DNA substrate does not inhibit catalysis, but instead provides energy that is efficiently utilized along the path to the transition state. It is suggested that this view can be usefully extended to "noncatalytic" site-specific DNA-binding proteins like transcriptional activators and general transcription factors.

Thermodynamics of DNA binding of MM17, a 'single chain dimer' of transcription factor MASH-1

MASH-1, a member of the basic helix-loop-helix (bHLH) family of transcriptional regulators, is a central factor for the regulation of the differentiation of committed neuronal precursor cells of the peripheral nervous system. We have previously produced MM17, a single chain version of this dimeric protein, by linking the C-terminal end of the first subunit to the N-terminal residue of the second subunit through a flexible peptide linker. We have now determined by isothermal titration calorimetry the thermodynamic parameters characterising the DNA binding reactions of MM17. The DNA binding specificity was relatively low and comparable to that observed for wild-type MASH bHLH. At 32 degrees C and pH 7, the concentration of MM17 at which 50% DNA binding occurred was determined as 22.8 and 152 nM for binding to MCK-S and the heterologous SP-1, respectively. Similarly to MASH bHLH the free energy of the association was only slightly temperature dependent, while both the entropy and the enthalpy change were strong functions of temperature. The free energy of DNA binding was independent of the pH for the pH range between 6 and 8. Dissection of the entropy change of the association reaction suggested that the two basic domains and the linker region between the subunits underwent a folding transition from a mainly unfolded to a predominantly ordered conformation. Therefore, like wild-type MASH bHLH, the DNA binding reaction of MM17 follows an induced fit mechanism.

Temperature dependence of DNA affinity chromatography of transcription factors

Oligonucleotides bound by the CAAT enhancer binding protein (C/EBP), the lactose repressor, and Gal4 were chemically coupled to cyanogen bromide-activated Sepharose and the temperature dependence of transcription factor chromatography was characterized. Each transcription factor was applied to the appropriate column and eluted using a salt gradient at several temperatures. Each transcription factor showed a unique behavior. As temperature was increases, less salt was required to elute C/EBP, more salt was required to elute lac repressor, while Gal4 showed a biphasic dependency with the amount of salt first decreasing between 4 and 19 degrees C and then increasing above 19 degrees C. This temperature dependence is not due to protein or DNA unfolding but rather is a property of complex formation. By loading a column, washing it at a permissive temperature, and then rapidly changing the column temperature, highly selective elution can be obtained. The thermodynamics of this temperature effect are different for the binding of specific and nonspecific DNA sequences, making chromatography at different temperatures a potentially important way of purifying transcription factors.

The energetics of HMG box interactions with DNA: thermodynamics of the DNA binding of the HMG box from mouse sox-5

The energetics of the Sox-5 HMG box interaction with DNA duplexes, containing the recognition sequence AACAAT, were studied by fluorescence spectroscopy, isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC). Fluorescence titration showed that the association constant of this HMG box with the duplexes is of the order 4x10(7) M(-1), increasing somewhat with temperature rise, i.e. the Gibbs energy is -40 kJ mol(-1) at 5 degrees C, decreasing to -48 kJ mol(-1) at 32 degrees C. ITC measurements of the enthalpy of association over this temperature range showed an endothermic effect below 17 degrees C and an exothermic effect above, suggesting a heat capacity change on binding of about -4 kJ K(-1) mol(-1), a value twice larger than expected from structural considerations. A straightforward interpretation of ITC data in heat capacity terms assumes, however, that the heat capacities of all participants in the association reaction do not change over the considered temperature range. Our previous studies showed that over the temperature range of the ITC experiments the HMG box of Sox-5 starts to unfold, absorbing heat and the heat capacities of the DNA duplexes also increase significantly. These heat capacity effects differ from that of the DNA/Sox-5 complex. Correcting the ITC measured binding enthalpies for the heat capacity changes of the components and complex yielded the net enthalpies which exhibit a temperature dependence of about -2 kJ K(-1) mol(-1), in good agreement with that predicted on the basis of dehydration of the protein-DNA interface. Using the derived heat capacity change and the enthalpy and Gibbs energy of association measured at 5 degrees C, the net enthalpy and entropy of association of the fully folded HMG box with the target DNA duplexes was determined over a broad temperature range. These functions were compared with those for other known cases of sequence specific DNA/protein association. It appears that the enthalpy and entropy of association of minor groove binding proteins are more positive than for proteins binding in the major groove. The observed thermodynamic characteristics of protein binding to the A+T-rich minor groove of DNA might result from dehydration of both polar and non-polar groups at the interface and release of counterions. The expected entropy of dehydration was calculated and found to be too large to be compensated by the negative entropy of reduction of translational/rotational freedom. This implies that DNA/HMG box association proceeds with significant decrease of conformational entropy, i.e. reduction in conformational mobility.

Thermodynamics of T cell receptor binding to peptide-MHC: evidence for a general mechanism of molecular scanning

Antigen-dependent activation of T lymphocytes requires T cell receptor (TCR)-mediated recognition of specific peptides, together with the MHC molecules to which they are bound. To achieve this recognition in a reasonable time frame, the TCR must scan and discriminate rapidly between thousands of MHC molecules differing from each other only in their bound peptides. Kinetic analysis of the interaction between a TCR and its cognate peptide-MHC complex indicates that both association and dissociation depend heavily on the temperature, indicating the presence of large energy barriers in both phases. Thermodynamic analysis reveals changes in heat capacity and entropy that are characteristic of protein-ligand associations in which local folding is coupled to binding. Such an "induced-fit" mechanism is characteristic of sequence-specific DNA-binding proteins that must also recognize specific ligands in the presence of a high background of competing elements. Here, we propose that induced fit may endow the TCR with its requisite discriminatory capacity and suggest a model whereby the loosely structured antigen-binding loops of the TCR rapidly explore peptide-MHC complexes on the cell surface until some critical structural complementarity is achieved through localized folding transitions. We further suggest that conformational changes, implicit in this model, may also propagate beyond the TCR antigen-binding site and directly affect self-association of ligated TCRs or TCR-CD3 interactions required for signaling.

Thermodynamics of nucleotide binding to NBS-I of the Bacillus subtilis preprotein translocase subunit SecA

SecA is the dissociatable nucleotide and preprotein binding subunit of the bacterial translocase. The thermodynamics of nucleotide binding to soluble SecA at nucleotide binding site I were determined by isothermal titration calorimetry. Binding of ADP and non-hydrolyzable ATPgammaS is enthalpy-driven (DeltaH(0) of -14.44 and -5.56 kcal/mol, respectively), but is accompanied by opposite entropic contributions (DeltaS(0) of -18.25 and 9.55 cal/mol/K, respectively). ADP binding results in a large change in the heat capacity of SecA (DeltaC(p)=-780 cal/mol/K). It is suggested that ADP binding promotes the interaction between the two thermodynamically discernible domains of SecA which is accompanied by a shielding of hydrophobic surface from solvent.

Single strand DNA specificity analysis of human nucleoside diphosphate kinase B

Nucleoside diphosphate kinases (NDP kinases) form a family of oligomeric enzymes present in all organisms. Eukaryotic NDP kinases are hexamers composed of identical subunits (approximately 17 kDa). A distinctive property of human NDPK-B encoded by the gene nm23-H2 is its ability to stimulate the gene transcription. This property is independent of its catalytic activity and is possibly related to the role of this protein in cellular events including differentiation and tumor metastasis. In this paper, we report the first characterization of human NDPK-B.DNA complex formation using a filter-binding assay and fluorescence spectroscopy. We analyzed the binding of several oligonucleotides mimicking the promoter region of the c-myc oncogene including variants in sequence, structure, and length of both strands. We show that NDPK-B binds to single-stranded oligonucleotides in a nonsequence specific manner, but that it exhibits a poor binding activity to double-stranded oligonucleotides. This indicates that the specificity of recognition to DNA is a function of the structural conformation of DNA rather than of its specific sequence. Moreover, competition experiments performed with all nucleotides provide evidence for the contribution of the six active sites in the DNA.protein complex formation. We propose a mechanism through which human NDPK-B could stimulate transcription of c-myc or possibly other genes involved in cellular differentiation.

Retention in reversed-phase chromatography: partition or adsorption?

A unified framework within the hermeneutics of the solvophobic theory is employed for the treatment of experimental data with nonpolar and weakly polar substances in reversed-phase chromatography (RPC), oil-water partitioning and adsorption on activated charcoal from dilute aqueous solution. This approach sheds light on the energetic similarities between such processes driven by the hydrophobic effect. Among several stationary phase models that have been proposed in the literature for the physical representation of alkyl-silica bonded phases, the isolated solvated hydrocarbon chains model is adopted for the retention in RPC since it represents most closely the stationary phase configuration and is not based a priori on a partition or adsorption mechanism as some other models are for the retention in RPC. Using the fundamental framework of the solvophobic theory, the free energy change per unit nonpolar surface area for octanol-water and hexadecane-water partitioning, retention in RPC as well as adsorption on activated charcoal from dilute aqueous solution at 25 degrees C are evaluated and they are found to be in good agreement with the corresponding experimental data. Furthermore, such quantities are very similar for all the above mentioned processes involving aqueous solution, in contradistinction to the predictions by the lattice theory. From the results it follows that these apparently disparate processes are subject to the same physicochemical principle. The present study demonstrates the capability of the solvophobic theory in describing the energetics of processes involving hydrophobic interactions, and exposes the difficulties in distinguishing between partition and adsorption mechanisms in RPC by using partition models based on the lattice approach. It is concluded that a clear distinction between partition and adsorption in RPC of nonpolar elites is not apparent from thermodynamic analysis.

DNA footprints of the two kinetically significant intermediates in formation of an RNA polymerase-promoter open complex: evidence that interactions with start site and downstream DNA induce sequential conformational changes in polymerase and DNA

Kinetic studies of formation and dissociation of open-promoter complexes (RPo) involving Esigma70 RNA polymerase (R) and the lambdaPR promoter (P) demonstrate the existence of two kinetically significant intermediates, designated I1 and I2, and facilitate the choice of conditions under which each accumulates. For such conditions, we report the results of equilibrium and transient DNase I and KMnO4 footprinting studies which characterize I1 and I2. At 0 degreesC, where extrapolation of equilibrium data indicates I1 is the dominant complex, DNA bases in the vicinity of the transcription start site (+1) do not react with KMnO4, indicating that this region is closed in I1. However, the DNA backbone in I1 is extensively protected from DNase I cleavage; the DNase I footprint extends approximately 30 bases downstream and at least approximately 40 bases upstream from the start site. I1 has a short lifetime (</=15 seconds), based on its sensitivity to competition with heparin. Shortly after a temperature downshift from 37 degreesC to 0 degreesC, in the time-range where we conclude that the dominant, transiently accumulated complex is I2, DNase I and KMnO4 footprinting reveal a complex with a closed-start site and an extended DNase I footprint like that of I1. However, unlike I1, I2 is insensitive to heparin competition and has a much longer dissociation lifetime at 0 degreesC. Based on footprinting, kinetic and thermodynamic studies, we conclude that in the short-lived intermediate I1 the promoter start site and downstream region are bound in a cleft defined by the open clamp-like jaws of Esigma70. We propose that binding of the start site and downstream DNA in this cleft triggers massive, relatively slow conformational changes which likely include RNA polymerase jaw closing with coupled folding. These proposed conformational changes occur prior to opening of the promoter start site region, and are responsible for the much longer lifetime of I2. Closing of the jaws of polymerase around the downstream region of promoter DNA appears to trigger opening of the start site region. From a quantitative analysis of the biphasic decay of KMnO4 reactivity of RPo at 0 degreesC, we obtain the equilibrium constant K3 for the conversion of I2 to RPo and the rate constant k-2 for the conversion of I2 to I1 (i.e. jaw opening). These quantitative results were previously unavailable at any temperature, and are necessary for the dissection of dissociation kinetic data at higher temperatures.

Electrostatic contributions to the binding free energy of the lambdacI repressor to DNA

A model based on the nonlinear Poisson-Boltzmann (NLPB) equation is used to study the electrostatic contribution to the binding free energy of the lambdacI repressor to its operator DNA. In particular, we use the Poisson-Boltzmann model to calculate the pKa shift of individual ionizable amino acids upon binding. We find that three residues on each monomer, Glu34, Glu83, and the amino terminus, have significant changes in their pKa and titrate between pH 4 and 9. This information is then used to calculate the pH dependence of the binding free energy. We find that the calculated pH dependence of binding accurately reproduces the available experimental data over a range of physiological pH values. The NLPB equation is then used to develop an overall picture of the electrostatics of the lambdacI repressor-operator interaction. We find that long-range Coulombic forces associated with the highly charged nucleic acid provide a strong driving force for the interaction of the protein with the DNA. These favorable electrostatic interactions are opposed, however, by unfavorable changes in the solvation of both the protein and the DNA upon binding. Specifically, the formation of a protein-DNA complex removes both charged and polar groups at the binding interface from solvent while it displaces salt from around the nucleic acid. As a result, the electrostatic contribution to the lambdacI repressor-operator interaction opposes binding by approximately 73 kcal/mol at physiological salt concentrations and neutral pH. A variety of entropic terms also oppose binding. The major force driving the binding process appears to be release of interfacial water from the protein and DNA surfaces upon complexation and, possibly, enhanced packing interactions between the protein and DNA in the interface. When the various nonelectrostatic terms are described with simple models that have been applied previously to other binding processes, a general picture of protein/DNA association emerges in which binding is driven by the nonpolar interactions, whereas specificity results from electrostatic interactions that weaken binding but are necessary components of any protein/DNA complex.

Stopped-flow kinetic analysis of the interaction of Escherichia coli RNA polymerase with the bacteriophage T7 A1 promoter

We have conducted a detailed kinetic and thermodynamic analysis of open complex formation between Escherichia coli RNA polymerase and the A1 promoter from bacteriophage T7 by monitoring alterations in the intrinsic protein fluorescence of RNA polymerase in stopped-flow kinetic studies. The stopped-flow kinetic data are consistent with a minimal model involving four steps for the formation of the open complex. Arrhenius plots for both the association and dissociation reactions for the equilibrium binding step leading to the formation of the closed complex were linear. With a positive van't Hoff enthalpy (DeltaHobs=18(+/-3) kcal mol-1) and a positive entropy (DeltaSobs=94(+/-15) e.u.) change for the equilibrium binding process, formation of the closed complex is entropy driven. The value of the apparent association rate constant for this binding step was approximately three orders of magnitude less than that expected for facilitated binding. Thus, a minimum of two steps is required to describe the formation of the closed complex. A fast facilitated binding step appears to be followed by a conformational change in RNA polymerase which leads to the formation of the closed complex. A non-linear Arrhenius plot obtained for the isomerization step in the conversion of the closed complex to an open one indicates that there are at least two steps in the conversion of the closed complex to an open one. We have assigned the apparent activation energy of 9.1(+/-1.9) kcal mol-1 to the step involving a conformational change in the protein and nucleation of strand separation and the apparent activation energy of 46(+/-12) kcal mol-1 to the step involving strand separation. At 37 degreesC, the value of the macroscopic isomerization rate constant (0.26(+/-0.02) s-1) in the conversion of the closed complex to an open one was an order of magnitude greater than the value reported in abortive initiation assays. This suggests that open complex formation is not the rate-determining step in the initiation of transcription in the case of the A1 promoter. To gain greater insight into the mechanism of initiation at the A1 promoter, we investigated the process of abortive product formation (pppApU) under conditions of non-saturating concentrations of the initiating nucleotide. A comparison of the lag times in the approach to the steady-state rate of abortive product formation when the reaction was initiated by the addition of UTP, ATP, the enzyme and the A1 promoter, respectively, indicates that the initiating nucleotide plays a key regulatory role in the initiation of transcription in the case of the A1 promoter.

Hydrophobic forces dominate the thermodynamic characteristics of UvrA-DNA damage interactions

The Escherichia coli DNA repair proteins UvrA, UvrB and UvrC work together to recognize and incise DNA damage during the process of nucleotide excision repair (NER). To gain an understanding of the damage recognition properties of UvrA, we have used fluorescence spectroscopy to study the thermodynamics of its interaction with a defined DNA substrate containing a benzo[a]pyrene diol epoxide (BPDE) adduct. Oligonucleotides containing a single site-specifically modified N2-guanine (+)-trans-, (-)-trans-, (+)-cis-, or (-)-cis-BPDE adducts were ligated into 50-base-pair DNA fragments. All four stereoisomers of DNA-BPDE adducts show an excitation maximum at 350 nm and an emission maximum around 380 to 385 nm. Binding of UvrA to the BPDE-DNA adducts results in a five to sevenfold fluorescence enhancement. Titration of the BPDE-adducted DNA with UvrA was used to generate binding isotherms. The equilibrium dissociation constants for UvrA binding to (+)-trans-, (-)-trans-, (+)-cis-, and (-)-cis- BPDE adduct were: 7.4+/-1.9, 15. 8+/-5.4, 11.3+/-2.7 and 22.4+/-2.0 nM, respectively. There was a large negative change in heat capacity DeltaCpo,obs, (-3.3 kcal mol-1 K-1) accompanied by a relatively unchanged DeltaGoobs with temperature. Furthermore, varying the concentration of KCl showed that the number of ions released upon formation of UvrA-DNA complex is about 3.4, a relatively small value compared to the contact size of UvrA with the substrate. These data suggest that hydrophobic interactions are an important driving force for UvrA binding to BPDE-damaged DNA.

The effects of salt on the TATA binding protein-DNA interaction from a hyperthermophilic archaeon

This study investigates the thermodynamics of the interaction of the TATA box binding protein (TBP) from Pyrococcus woesei (Pw) with an oligonucleotide containing a specific binding site. Pw is a hyperthermophilic archeal organism which exists under conditions of high salt and high temperature. A measurable protein-DNA interaction only occurs at high salt concentrations. Isothermal titration calorimetric binding studies were performed under a range of salts (potassium chloride, potassium phosphate, potassium acetate and sodium acetate) at varying concentrations (0.8 to 1.6 M). At the high salt concentrations used the observed equilibrium binding constant increases with increasing salt concentration. This is very different to the effect reported for all other protein-DNA interactions which have been studied at lower salt concentrations. Thermodynamic data suggest that the protein-DNA interaction at high salt concentration is accompanied by the removal of large numbers of water molecules from the buried hydrophobic surface area. In addition, the involvement of ions appears to influence the binding which can be explained by binding of cations in the interface between the electrostatically negative lateral lobes on the protein and the negatively charged DNA.

Calorimetric studies of E. coli SSB protein-single-stranded DNA interactions. Effects of monovalent salts on binding enthalpy

Isothermal titration calorimetry (ITC) was used to examine the effects of monovalent salts (NaCl, NaBr, NaF and ChCl) on the binding enthalpy (DeltaHobs) for E. coli SSB tetramer binding to the single-stranded oligodeoxythymidylates, dT(pT)69 and dT(pT)34 over a wide range of salt concentrations from 10 mM to 2.0 M (25 degrees C, pH 8.1), and when possible, the binding free energy and entropy (DeltaG degrees obs, DeltaS degrees obs). At low monovalent salt concentrations (<0.1 M), the total DeltaHobs for saturating all sites on the SSB tetramer with ssDNA shows little dependence on salt concentration, but is extremely large and exothermic (DeltaHobs=-150(+/-5) kcal/mol). This is much larger than any DeltaHobs previously reported for a protein-nucleic acid interaction. However, at salt concentrations above 0.1 M, DeltaHobs is quite sensitive to NaCl and NaBr concentration, becoming less negative with increasing salt concentration (DeltaHobs=-70(+/-1)-kcal/mol in 2 M NaBr). These salt effects on DeltaHobs were mainly a function of anion type and concentration, with the largest effects observed in NaBr, and then NaCl, with little effect of [NaF]. These large effects of salt on DeltaHobs appear to be coupled to a net release of weakly bound anions (Br- and Cl-) from the SSB protein upon DNA binding. However, at lower salt concentrations (</=0.1 M), specific cation effects on DeltaHobs also are observed. Under conditions where we can determine DeltaG degrees obs, DeltaS degrees obs, and DeltaHobs (25 degrees C, pH 8.1, 0.17 to 2 M NaBr), SSB binding to dT(pT)69 is enthalpically driven with a large unfavorable entropic contribution, both of which are dependent upon [NaBr]. These studies show that weak anion binding to a protein can result in large effects of salt concentration on DeltaHobs (as well as DeltaG degrees obs and DeltaS degrees obs) for a protein-ssDNA interaction. The possibility of such effects needs to be considered in any interpretation of the thermodynamics of this and other protein-nucleic acid interactions.

Isothermal titration calorimetric studies on the associations of putidaredoxin to NADH-putidaredoxin reductase and P450cam

Putidaredoxin (Pdx), an iron-sulfur protein containing a 2Fe-2S cluster, serves as a physiological electron mediator from NADH-putidaredoxin reductase (PdR) to P450cam in the P450cam monooxygenation reaction cycle. Previous studies have revealed that the associations of Pdx with P450cam and PdR are not strongly dominated by electrostatic interactions, although such interactions stabilize most electron-transfer complexes [A.R. De Pascalis, I. Jelesarov, F. Ackermann, W.H. Koppenol, M. Hiroasawa, D.B. Knaff, H.R. Bosshard, Protein Sci. 2 (1993) 1126-1135]. In the present study, to elucidate the interactions dominating the specific associations in the electron-transfer reaction mediated by Pdx, the thermodynamic properties--entropy (delta S), enthalpy (delta H), and heat capacity changes (delta Cp)--for PdR/Pdx and P450cam/Pdx association reactions have been examined by isothermal titration calorimetry (ITC). Although the binding enthalpy change, delta Hbind, for the PdR/Pdx association is positive at 10 degrees C, it declines linearly with temperature in the range 10-22 degrees C and becomes negative above 11 degrees C. On the other hand, the binding entropy change, delta Sbind, is positive at all temperatures examined in this study, indicating that the association of Pdx to PdR is entropically driven. On the basis of the temperature dependence of delta Hbind, delta Cpbind for the association of Pdx to PdR was estimated as -1.24 kJ mol-1 K-1. This value is larger than those reported for other electron-transfer protein systems (e.g., -0.68 kJ mol-1 K-1 for ferredoxin/ferredoxin: NADP+ reductase), suggesting that the PdR/Pdx association may be dominated by hydrophobic rather than electrostatic components. For the P450cam/Pdx association, the negative delta Sbind and highly favorable delta Hbind were observed, behavior that stands in sharp contrast to the association reactions in other electron-transfer proteins. The energetics of the P450cam/Pdx association are similar to those of binding reaction of antibody to antigen in which van der Waals and hydrogen bonding interactions are dominant, resulting in high specificity in the association of Pdx with P450cam.

Transcription factor-green fluorescent protein chimeric fusion proteins and their use in studies of DNA affinity chromatography

A new plasmid, pJ22, was produced by introducing the enhanced green fluorescent protein (GFP) coding sequence into the pET28 plasmid while retaining much of the multiple cloning site. This new plasmid was then used to produce a chimeric fusion protein containing the DNA-binding region of the rat liver CAAT enhancer binding protein (C/EBP) fused to the COOH-terminus of GFP. This new GFP-C/EBP fusion protein also contains (His)6 to facilitate purification by Ni(2+)-agarose and several other useful features. The plasmid and protein were developed to allow us to more rapidly investigate the DNA-Sepharose affinity chromatography of transcription factors. The GFP-C/EBP protein is virtually identical in its DNA-binding properties to a well-characterized, bacterially expressed protein called C/EBP 62 which has been shown to mimic rat wild-type C/EBP DNA-binding. GFP-C/EBP also binds to DNA-Sepharose which contains the CAAT element and is eluted by a salt gradient. Salt-dependent elution was highly temperature-dependent over the range of 4-19 degrees C. Since temperature-dependent DNA-binding has also been reported for other DNA-binding proteins, this may also occur with other transcription factors. DNA-affinity chromatography gave higher purity than that obtained by Ni(2+)-agarose chromatography and chromatography on the same DNA-Sepharose column at two different temperatures resulted in the greatest purification, to near homogeneity. This temperature-dependent affinity chromatography provides an important new approach to transcription factor purification.

Thermodynamic characterization of non-sequence-specific DNA-binding by the Sso7d protein from Sulfolobus solfataricus

We used isothermal titration calorimetry and fluorescence spectroscopy to investigate the thermodynamics of non-sequence-specific DNA-binding by the Sso7d protein from the archaeon Sulfolobus solfataricus. We report the Sso7d-poly(dGdC) binding thermodynamics as a function of buffer composition (Tris-HCl or phosphate), temperature (15 to 45 degrees C), pH (7.1 to 8.0), osmotic stress and solvent (H2O/2H2O), and compare it to poly (dAdT) binding; and we have previously also reported the salt concentration dependence. Binding isotherms can be represented by the McGhee-von Hippel model for non-cooperative binding, with a binding site size of four to five DNA base-pairs and binding free energies in the range DeltaG degrees approximately -7 to DeltaG degrees approximately -10 kcal mol-1, depending on experimental conditions. The non-specific nature of the binding is reflected in similar thermodynamics for binding to poly(dAdT) and poly(dGdC). The native lysine methylation of Sso7d has only minor effects on the binding thermodynamics. Sso7d binding to poly(dGdC) is endothermic at 25 degrees C with a binding enthalpy DeltaH degrees approximately 10 kcal mol-1 in both phosphate and Tris-HCl buffers at pH 7.6, indicating that DeltaH degrees does not include large contributions from coupled buffer ionization equilibria at this pH. The binding enthalpy is temperature dependent with a measured heat capacity change DeltaCp degrees=-0.25(+/-0.01) kcal mol-1 K-1 and extrapolations of thermodynamic data indicate that the complex is heat stable with exothermic binding close to the growth temperature (75 to 80 degreesC) of S. solfataricus. Addition of neutral solutes (osmotic stress) has minor effects on DeltaG degrees and the exchange of H2O for 2H2O has only a small effect on DeltaH degrees, consistent with the inference that complex formation is not accompanied by net changes in surface hydration. Thus, other mechanisms for the heat capacity change must be found. The observed thermodynamics is discussed in relation to the nature of non-sequence-specific DNA-binding by proteins.

Changes in solvation during DNA binding and cleavage are critical to altered specificity of the EcoRI endonuclease

Restriction endonucleases such as EcoRI bind and cleave DNA with great specificity and represent a paradigm for protein-DNA interactions and molecular recognition. Using osmotic pressure to induce water release, we demonstrate the participation of bound waters in the sequence discrimination of substrate DNA by EcoRI. Changes in solvation can play a critical role in directing sequence-specific DNA binding by EcoRI and are also crucial in assisting site discrimination during catalysis. By measuring the volume change for complex formation, we show that at the cognate sequence (GAATTC) EcoRI binding releases about 70 fewer water molecules than binding at an alternate DNA sequence (TAATTC), which differs by a single base pair. EcoRI complexation with nonspecific DNA releases substantially less water than either of these specific complexes. In cognate substrates (GAATTC) kcat decreases as osmotic pressure is increased, indicating the binding of about 30 water molecules accompanies the cleavage reaction. For the alternate substrate (TAATTC), release of about 40 water molecules accompanies the reaction, indicated by a dramatic acceleration of the rate when osmotic pressure is raised. These large differences in solvation effects demonstrate that water molecules can be key players in the molecular recognition process during both association and catalytic phases of the EcoRI reaction, acting to change the specificity of the enzyme. For both the protein-DNA complex and the transition state, there may be substantial conformational differences between cognate and alternate sites, accompanied by significant alterations in hydration and solvent accessibility.

Thermodynamics of specific and non-specific DNA binding by the c-Myb DNA-binding domain

The thermodynamics of the c-Myb DNA-binding domain (R2R3) interaction with its target DNA have been analyzed using isothermal titration calorimetry and amino acid mutagenesis. The enthalpy of association between the standard R2R3, the Cys130 mutant substituted with Ile, and the cognate DNA is -12.5 (+/- 0.1) kcal mol-1 at pH 7.5 and at 20 degrees C, and this interaction is enthalpically driven throughout the physiological temperature range. In order to understand the DNA recognition mechanism, several pairs of interactions were investigated using single and multiple-base alterations with single and multiple-amino acid substituted mutants. The interactions between the standard R2R3 and many non-cognate DNAs were accompanied by binding enthalpy changes and heat capacity changes, although their affinities were reduced. The roles of the electrostatic interactions in binding to the cognate and the non-cognate DNAs were also analyzed from the dependency of the thermodynamic parameters on the salt concentration. The heat capacity change was found to be significantly dependent upon the salt concentration. Several mutant proteins bound to the multiple-base altered DNA with very small enthalpy changes, although they bound to the cognate and the single-base altered DNAs with detectable enthalpy and heat capacity changes. From the thermodynamic cycles derived from the DNA binding of the amino acid substituted R2R3 to the base substituted DNA duplexes, the individual thermodynamic mechanisms of the specific DNA recognition of R2R3 were dissected. The local folding mechanism was highlighted by the substitution of Pro with either Gly or Ala at the linker between R2 and R3. The characteristic thermodynamic features of specific and non-specific DNA binding are discussed.

DNA-binding characteristics of the Escherichia coli CytR regulator: a relaxed spacing requirement between operator half-sites is provided by a flexible, unstructured interdomain linker

The Escherichia coli CytR regulator belongs to the LacI family of sequence-specific DNA-binding proteins and prevents CRP-mediated transcription in the CytR regulon. Unlike the other members of this protein family, CytR binds with only modest affinity to its operators and transcription repression thus relies on the formation of nucleoprotein complexes with the cAMP-CRP complex. Moreover, CytR exhibits a rotational and translational flexibility in operator binding that is unprecedented in the LacI family. In this report we examined the effect of changing the spacing between CytR half-operators on CytR regulation in vivo and on CytR binding in vitro. Maximum repression was seen with the short spacing variants: repression peaks when the half-operators lie on the same face of the DNA helix. Repression was retained for most spacing variants with centre separations of half-operators < or = 3 helical turns. Our data confirm and extend the view that CytR is a highly flexible DNA binder that can adapt many different conformations for co-operative binding with CRP. Furthermore, limited proteolysis of radiolabelled CytR protein showed that the interdomain linker connecting the DNA binding domains and the core part of CytR does not become structured upon DNA binding. We conclude that CytR does not use hinge alpha-helices for minor groove recognition. Rather, CytR possesses a highly flexible interdomain linker that allows it to form complexes with CRP at promoters with quite different architecture.

Comparison of simulated and experimentally determined dynamics for a variant of the Lacl DNA-binding domain, Nlac-P

Recent advances in the experimentally determined structures and dynamics of the domains within LacI provide a rare context for evaluating dynamics calculations. A 1500-ps trajectory was simulated for a variant of the LacI DNA-binding domain, which consists of the first three helices in LacI and the hinge helix of the homologous PurR. Order parameters derived from dynamics simulations are compared to those obtained for the LacI DNA-binding domain with 15N relaxation NMR spectroscopy (Slijper et al., 1997. Biochemistry. 36:249-254). The MD simulations suggest that the unstructured loop between helices II and III does not exist in a discrete state under the conditions of no salt and neutral pH, but occupies a continuum of states between the DNA-bound and free structures. Simulations also indicate that the unstructured region between helix III and the hinge helix is very mobile, rendering motions of the hinge helix essentially independent of the rest of the protein. Finally, the alpha-helical hydrogen bonds in the hinge helix are broken after 1250 ps, perhaps as a prelude to helix unfolding.

Calorimetric observation of a GroEL-protein binding reaction with little contribution of hydrophobic interaction

Binding of Escherichia coli chaperonin, GroEL, to substrate proteins with non-native structure, reduced alpha-lactalbumin (rLA) and denatured pepsin, were analyzed by isothermal titration calorimetry at various temperatures in the presence of salt (0.2 M KCl). Both proteins bound to GroEL with 1:1 stoichiometry and micromolar affinity at all temperatures tested. However, thermodynamic properties of their binding to GroEL are remarkably different from each other. While heat capacity changes (DeltaCp) of rLA-GroEL binding showed large negative values, -4.19 kJ mol-1 K-1, that of denatured pepsin-GroEL binding was only -0.2 kJ mol-1 K-1. These values strongly indicate that the hydrophobic interaction is a major force of rLA-GroEL binding but not so for denatured pepsin-GroEL binding. When salt was omitted from the solution, the affinity and DeltaCp of the rLA-GroEL binding reaction were not significantly changed whereas denatured pepsin lost affinity to GroEL. Thus, in the non-native protein-GroEL binding reaction, thermodynamic properties, as well as the effect of salt, differ from protein to protein and hydrophobic interaction may not always be a major driving force.

Electrostatic and non-electrostatic contributions to the binding free energies of anthracycline antibiotics to DNA

The knowledge about molecular factors driving simple ligand-DNA interactions is still limited. The aim of the present study was to investigate the electrostatic and non-electrostatic contributions to the binding free energies of anthracycline compounds with DNA. Theoretical calculations based on continuum methods (Poisson-Boltzmann and solvent accessible surface area) were performed to estimate the binding free energies of five selected anthracycline ligands (daunomycin, adriamycin, 9-deoxyadriamycin, hydroxyrubicin, and adriamycinone) to DNA. The free energy calculations also took into account the conformational change that DNA undergoes upon ligand binding. This conformational change appeared to be very important for estimating absolute free energies of binding. Our studies revealed that the absolute values of all computed contributions to the binding free energy were quite large compared to the total free energy of binding. However, the sum of these large positive and negative values produced a small negative value of the free energy around -10 kcal/mol. This value is in good agreement with experimental data. Experimental values for relative binding free energies were also reproduced for charged ligands by our calculations. Together, it was found that the driving force for ligand-DNA complex formation is the non-polar interaction between the ligand and DNA even if the ligand is positively charged.

Lactose repressor protein: functional properties and structure

The lactose repressor protein (LacI), the prototype for genetic regulatory proteins, controls expression of lactose metabolic genes by binding to its cognate operator sequences in E. coli DNA. Inducer binding elicits a conformational change that diminishes affinity for operator sequences with no effect on nonspecific binding. The release of operator is followed by synthesis of mRNA encoding the enzymes for lactose utilization. Genetic, chemical and physical studies provided detailed insight into the function of this protein prior to the recent completion of X-ray crystallographic structures. The structural information can now be correlated with the phenotypic data for numerous mutants. These structures also provide the opportunity for physical and chemical studies on mutants designed to examine various aspects of lac repressor structure and function. In addition to providing insight into protein structure-function correlations, LacI has been utilized in a wide variety of applications both in prokaryotic gene expression and in eukaryotic gene regulation and studies of mutagenesis.

The trp RNA-binding attenuation protein (TRAP) from Bacillus subtilis binds to unstacked trp leader RNA

TRAP (trp RNA-binding attenuation protein) is a tryptophan-activated RNA-binding protein that regulates expression of the trp biosynthetic genes by binding to a series of GAG and UAG trinucleotide repeats generally separated by two or three spacer nucleotides. Previously, we showed that TRAP contains 11 identical subunits arranged in a symmetrical ring. Based on this structure, we proposed a model for the TRAP.RNA interaction where the RNA wraps around the protein with each repeat of the RNA contacting one or a combination of two adjacent subunits of the TRAP oligomer. Here, we have shown that RNAs selected in vitro based on their ability to bind tryptophan-activated TRAP contain multiple G/UAG repeats and show a strong bias for pyrimidines as the spacer nucleotides between these repeats. The affinity of the TRAP.RNA interaction displays a nonlinear temperature dependence, increasing between 5 degrees C and 47 degrees C and then decreasing from 47 degrees C to 67 degrees C. Differential scanning calorimetry and circular dichroism spectroscopy demonstrate that TRAP is highly thermostable with few detectable changes in the structure between 25 degrees C and 70 degrees C, suggesting that the temperature dependence of this interaction reflects changes in the RNA. Results from circular dichroism and UV absorbance spectroscopy support this hypothesis, demonstrating that trp leader RNA becomes unstacked upon binding TRAP. We propose that the bias toward pyrimidines in the spacer nucleotides of the in vitro selected RNAs represents the inability of Us and Cs to form stable base stacking interactions, which allows the flexibility needed for the RNA to wrap around the TRAP oligomer.

Binding of the estrogen receptor to DNA. The role of waters

Molecular dynamics simulations are carried out to investigate the binding of the estrogen receptor, a member of the nuclear hormone receptor family, to specific and non-specific DNA. Two systems have been simulated, each based on the crystallographic structure of a complex of a dimer of the estrogen receptor DNA binding domain with DNA. One structure includes the dimer and a consensus segment of DNA, ds(CCAGGTCACAGTGACCTGG); the other structure includes the dimer and a nonconsensus segment of DNA, ds(CCAGAACACAGTGACCTGG). The simulations involve an atomic model of the protein-DNA complex, counterions, and a sphere of explicit water with a radius of 45 A. The molecular dynamics package NAMD was used to obtain 100 ps of dynamics for each system with complete long-range electrostatic interactions. Analysis of the simulations revealed differences in the protein-DNA interactions for consensus and nonconsensus sequences, a bending and unwinding of the DNA, a slight rearrangement of several amino acid side chains, and inclusion of water molecules at the protein-DNA interface region. Our results indicate that binding specificity and stability is conferred by a network of direct and water mediated protein-DNA hydrogen bonds. For the consensus sequence, the network involves three water molecules, residues Glu-25, Lys-28, Lys-32, Arg-33, and bases of the DNA. The binding differs for the nonconsensus DNA sequence in which case the fluctuating network of hydrogen bonds allows water molecules to enter the protein-DNA interface. We conclude that water plays a role in furnishing DNA binding specificity to nuclear hormone receptors.

Electrostatics and hydration at the homeodomain-DNA interface: chemical probes of an interfacial water cavity

Electrostatics and hydration of a homeodomain-DNA complex are dissected by chemical modification. Selective neutralization of phosphate charges by methylphosphonate substitution demonstrates the differential importance of short- and long-range electrostatic interactions. Whereas the footprint of direct contacts is in accord with crystal structures, interference is also observed at non-contacted sites. Such sites adjoin a novel interfacial water cavity in the major groove. Non-contacted phosphodiester groups in the cavity are proposed to contribute to long-range ordering of an extended protein-water-DNA interface. Use of isolated S(p) and R(p) methylphosphonate diastereomers demonstrates that interference at this extended interface is stereoselective and charge-independent. Attenuation of protein binding presumably reflects groove-specific reorganization of bound water. Surprisingly, such attenuation can exceed that due to neutralization of a direct phosphate-side-chain salt bridge. These results support the hypothesis that hydration of an interfacial cavity functions as a non-covalent extension of the DNA surface. Stereo-specific interrogation of bound water by chemical synthesis provides a general method to assess the coupling between solvation and DNA recognition.

Thermodynamics of the interactions of lac repressor with variants of the symmetric lac operator: effects of converting a consensus site to a non-specific site

What are the thermodynamic consequences of the stepwise conversion of a highly specific (consensus) protein-DNA interface to one that is nonspecific? How do the magnitudes of key favorable contributions to complex stability (burial of hydrophobic surfaces and reduction of DNA phosphate charge density) change as the DNA sequence of the specific site is detuned? To address these questions we investigated the binding of lac repressor (LacI) to a series of 40 bp fragments carrying symmetric (consensus) and variant operator sequences over a range of temperatures and salt concentrations. Variant DNA sites contained symmetrical single and double base-pair substitutions at positions 4 and/or 5 [sequence: see text] in each 10 bp half site of the symmetric lac operator (Osym). Non-specific interactions were examined using a 40 bp non-operator DNA fragment. Disruption of the consensus interface by a single symmetrical substitution reduces the observed equilibrium association constant (K(obs)) for Osym by three to four orders of magnitude; double symmetrical substitutions approach the six orders in magnitude difference between specific and non-specific binding to a 40 bp fragment. At these adjacent positions in the consensus site, the free energy effects of multiple substitutions are non-additive: the first reduces /deltaG(obs)o/ by 3 to 5 kcal mol(-1), approximately halfway to the non-specific level, whereas the second is less deleterious, reducing /deltaG(obs)o/ by less than 3 kcal mol(-1). Variant-specific dependences of K(obs) on temperature and salt concentration characterize these LacI-operator interactions. In general, binding constants and standard free energies of binding both exhibit characteristic extrema near 290 K. As a consequence, both the enthalpic and entropic contributions to stability of Osym and variant complexes change from positive (i.e. entropy driven) at lower temperatures to negative (i.e. enthalpy driven) at higher temperatures, indicating that the heat capacity change upon binding, deltaC(obs)o, is large and negative. In general, /deltaC(obs)o/ decreases as the specificity and stability of the variant complex decreases. Stabilities of complexes of LacI with Osym and all variant operators are strongly [salt]-dependent. Binding constants for the variant complexes exhibit a power-dependence on [salt] that is larger in magnitude (i.e. more negative) than for Osym, but no obvious trend relates changes in contributions from the polyelectrolyte effect and the observed reductions in stability (delta deltaG(obs)o). These variant-specific thermodynamic signatures provide novel insights into the consequences of converting a consensus interface to a less specific one; such insights are not obtained from comparisons at the level of delta deltaG(obs)o. We propose that this variant-specific behavior arises from a strong effect of operator sequence on the extent of induced conformational changes in the protein (and possibly also in the DNA site) which accompany binding.

Thermodynamic analysis of monoclonal antibody binding to duplex DNA

A technique based on fluorescence polarization (anisotropy) was used to measure the binding of antibodies to DNA under a variety of conditions. Fluorescein-labeled duplexes of 20 bp in length were employed as the standard because they are stable even at low ionic strength yet sufficiently short so that both arms of an IgG cannot bind to the same duplex. IgG Jel 274 binds duplexes in preference to single-stranded DNA; in 80 mM NaCl Kobs for (dG)20.(dC)20 is 4.1x10(7) M-1 compared with 6.4x10(5) M-1 for d(A5C10A5). There is little sequence specificity, but the interaction is very dependent on ionic strength. From plots of log Kobs against log[Na+] it was deduced that five or six ion pairs are involved in complex formation. At low ionic strength,Kobs is independent of temperature and complex formation is entropy driven with DeltaH degrees obs and DeltaC degrees p,obs both zero. In contrast, in 80 mM NaCl DeltaC degrees p,obs is -630 and -580 cal mol-1K-1 for [d(TG)]10.[d(CA)]10 and (dG)20.(dC)20 respectively. IgG Jel 241 also binds more tightly to duplexes than single-stranded DNA, but sequence preferences were apparent. The values for Kobs to [d(AT)]20 and [d(GC)]20 are 2.7x10(8) and 1.3x10(8) M-1 respectively compared with 5.7x10(6) M-1 for both (dA)20. (dT)20 and (dG)20.(dC)20. As with Jel 274, the binding of Jel 241 is very dependent on ionic strength and four or five ionic bonds are involved in complex formation with all the duplex DNAs which were tested. DeltaC degrees p,obs for Jel 241 binding to [d(AT)]20 was negative (-87 cal mol-1K-1) in 80 mM NaCl but was zero at high ionic strength (130 mM NaCl). Therefore, for duplex-specific DNA binding antibodies DeltaC degrees p,obs is dependent on [Na+] and a large negative value does not correlate with sequence-specific interactions.

Electrostatic effects in homeodomain-DNA interactions

We report here an investigation of the role of electrostatics in homeodomain-DNA interactions using techniques based around the use of the Poisson-Boltzmann equation. In the present case such a study is of particular interest, since in contrast to other proteins previously studied with this method, the homeodomain is a small, highly charged protein that forms extensive ion pairs upon binding DNA. We have investigated the salt dependence of the binding constant for specific association and for a variety of models for non-specific association. The results indicate that, in line with the models proposed by Manning and Record, the entropy of counterion release accounts for a significant fraction of the salt dependence of the binding free energy, though this is perhaps due to fortuitous cancellation of other contributing terms. The thermodynamic effects of a number of specific homeodomain mutants were also investigated, and partly rationalized in terms of favorable electrostatic interactions in the major goove of DNA. Investigation of the temperature-dependence of the free energy of association indicates that the electrostatic contributions become increasingly favorable as the temperature rises. For this particular system, however, there appears to be no significant electrostatic contribution to the delta(delta C(p)) of association. Finally, an analysis of the free energy of interaction when the homeodomain is moved ca one Debye length from the DNA suggests that pure electrostatic forces are able to steer the homeodomain into a partially correct orientation for binding to the DNA.

Thermodynamic and structural analysis of microtubule assembly: the role of GTP hydrolysis

Different models have been proposed that link the tubulin heterodimer nucleotide content and the role of GTP hydrolysis with microtubule assembly and dynamics. Here we compare the thermodynamics of microtubule assembly as a function of nucleotide content by van't Hoff analysis. The thermodynamic parameters of tubulin assembly in 30-100 mM piperazine-N,N'-bis(2-ethanesulfonic acid), 1 mM MgSO4, 2 mM EGTA, pH 6.9, in the presence of a weakly hydrolyzable analog, GMPCPP, the dinucleotide analog GMPCP plus 2 M glycerol, and GTP plus 2 M glycerol were obtained together with data for taxol-GTP/GDP tubulin assembly (GMPCPP and GMPCP are the GTP and GDP nucleotide analogs where the alpha beta oxygen has been replaced by a methylene, -CH2-). All of the processes studied are characterized by a positive enthalpy, a positive entropy, and a large, negative heat capacity change. GMPCP-induced assembly has the largest negative heat capacity change and GMPCPP has the second largest, whereas GTP/2 M glycerol- and taxol-induced assembly have more positive values, respectively. A large, negative heat capacity is most consistent with the burial of water-accessible hydrophobic surface area, which gives rise to the release of bound water. The heat capacity changes observed with GTP/2 M glycerol-induced and with taxol-induced assembly are very similar, -790 +/- 190 cal/mol/k, and correspond to the burial of 3330 +/- 820 A2 of nonpolar surface area. This value is shown to be very similar to an estimate of the buried nonpolar surface in a reconstructed microtubule lattice. Polymerization data from GMPCP- and GMPCPP-induced assembly are consistent with buried nonpolar surface areas that are 3 and 6 times larger. A linear enthalpy-entropy and enthalpy-free energy plot for tubulin polymerization reactions verifies that enthalpy-entropy compensation for this system is based upon true biochemical correlation, most likely corresponding to a dominant hydrophobic effect. Entropy analysis suggests that assembly with GTP/2 M glycerol and with taxol is consistent with conformational rearrangements in 3-6% of the total amino acids in the heterodimer. In addition, taxol binding contributes to the thermodynamics of the overall process by reducing the delta H degree and delta S degree for microtubule assembly. In the presence of GMPCPP or GMPCP, tubulin subunits associate with extensive conformational rearrangement, corresponding to 10% and 26% of the total amino acids in the heterodimer, respectively, which gives rise to a large loss of configurational entropy. An alternative, and probably preferable, interpretation of these data is that, especially with GMPCP-tubulin, additional isomerization or protonation events are induced by the presence of the methylene moiety and linked to microtubule assembly. Structural analysis shows that GTP hydrolysis is not required for sheet closure into a microtubule cylinder, but only increases the probability of this event occurring. Sheet extensions and sheet polymers appear to have a similar average length under various conditions, suggesting that the minimum cooperative unit for closure of sheets into a microtubule cylinder is approximately 400 nm long. Because of their low level of occurrence, sheets are not expected to significantly affect the thermodynamics of assembly.

The role of water in protein-DNA interactions

It is becoming well accepted that water plays an important role in both the specificity and affinity of protein-DNA interactions. Recently, a combination of structural, biochemical and thermodynamic techniques has particularly enhanced our understanding of the role of water in complexes between DNA and three different proteins: the trp repressor; the homeodomain; and the glucocorticoid receptor DNA-binding domain.

Eukaryotic transcription factor-DNA complexes

Eukaryotes have three distinct RNA polymerases that catalyze transcription of nuclear genes. RNA polymerase II is responsible for transcribing nuclear genes encoding the messenger RNAs and several small nuclear RNAs. Like RNA polymerases I and III, polymerase II cannot recognize its target promoter directly and initiate transcription without accessory factors. Instead, this large multisubunit enzyme relies on general transcription factors and transcriptional activators and coactivators to regulate transcription from class II promoters. X-ray crystallography and nuclear magnetic resonance spectroscopy have been used to study complexes of general transcription factors and transcriptional activators with their specific DNA targets. This work has provided important structural insights into transcription initiation by polymerase II and the more general problem of DNA sequence recognition.

Charge density-dependent strength of hydration and biological structure

Small ions of high charge density (kosmotropes) bind water molecules strongly, whereas large monovalent ions of low charge density (chaotropes) bind water molecules weakly relative to the strength of water-water interactions in bulk solution. The standard heat of solution of a crystalline alkali halide is shown here to be negative (exothermic) only when one ion is a kosmotrope and the ion of opposite charge is a chaotrope; this standard heat of solution is known to become proportionally more positive as the difference between the absolute heats of hydration of the corresponding gaseous anion and cation decreases. This suggests that inner sphere ion pairs are preferentially formed between oppositely charged ions with matching absolute enthalpies of hydration, and that biological organization arises from the noncovalent association of moieties with matching absolute free energies of solution, except where free energy is expended to keep them apart. The major intracellular anions (phosphates and carboxylates) are kosmotropes, whereas the major intracellular monovalent cations (K+; arg, his, and lys side chains) are chaotropes; together they form highly soluble, solvent-separated ion pairs that keep the contents of the cell in solution.

Characterization of a cooperativity domain mutant Lys3 --> Ala (K3A) T4 gene 32 protein

The N-terminal "B" domain of T4 gene 32 protein contains a Lys3-Arg4-Lys5 sequence that has been postulated to provide a major determinant for cooperative binding. In this report, the equilibrium binding properties of a Lys3 --> Ala substitution mutant of gp32 (K3A gp32) and described and compared to a set of substitution mutants of Arg4 previously described (Villemain, J. L., and Giedroc, D. P. (1993) Biochemistry 32, 11235-11246) and further characterized here. K3A gp32 exhibits binding behavior which mirrors that of R4Q gp32. Despite an 6-8-fold decrease in overall binding affinity (Kapp = Kint x omega) at pH 8.1, 0.20 M NaCl, 20 degrees C, the magnitude of the cooperativity parameter is at most 2-3-fold smaller than that of the wild-type protein. The magnitude of omega is independent of salt concentration and type over the range in [NaCl] from 0.125 to 0. 225 M and [NaF] between 0.20 and 0.32 M (log omega = 2.86 +/- 0.19). For comparison, log omega for wild-type gp32 is 2.91 (+/- 0.21) resolved at 0.275 M NaCl and 3.39 +/- 0.11 in [NaF] between 0.40 and 0.45 M. In contrast to omega, the [NaCl] dependence of Kapp is large and markedly nonlinear for both wild-type and K3A gp32s over a [NaCl] range extending from 0.05 M to 0.40 M NaCl. Modeling of the complete salt dependence of Kapp for wild-type, K3A, and R4T gp32s in NaCl and NaF with a simple ion-exchange model suggests that substitutions within the basic Lys3-Arg4-Lys5 sequence do not strongly modulate the net displacement of cations and anions upon poly(A) complex formation by gp32.