Contribution to the thermodynamics of protein folding from the reduction in water-accessible nonpolar surface area

Protein folding and the transfer of hydrocarbons from a dilute aqueous solution to the pure liquid phase are thermodynamically similar in that both processes remove nonpolar surface from water and both are accompanied by anomalously large negative heat capacity changes. On the basis of a limited set of published surface areas, we previously proposed that heat capacity changes (delta C degrees p) for the transfer of hydrocarbons from water to the pure liquid phase and for the folding of globular proteins exhibit the same proportionality to the reduction in water-accessible nonpolar surface area (delta Anp) [Spolar, R.S., Ha, J.H., & Record, M.T., Jr. (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 8382-8385]. The consequence of this proposal is that the experimental delta C degrees p for protein folding can be used to obtain estimates of delta Anp and of the contribution to the stability of the folded state from removal of a nonpolar surface from water. In this paper, a rigorous molecular surface area algorithm [Richmond, T.J. (1984) J. Mol. Biol. 178, 63-89] is applied to obtain self-consistent values of the water-accessible nonpolar surface areas of the native and completely denatured states of the entire set of globular proteins for which both crystal structures and delta C degrees p of folding have been determined and for the set of liquid and liquefiable hydrocarbons for which delta C degrees p of transfer are known. Both processes (hydrocarbon transfer and protein folding) exhibit the same direct proportionality between delta C degrees p and delta Anp. We conclude that the large negative heat capacity changes observed in protein folding and other self-assembly processes involving proteins provide a quantitative measure of the reduction in the water-accessible nonpolar surface area and of the contribution of the hydrophobic effect to the stability of the native state and to protein assembly.     

Equilibrium constants under physiological conditions for the reactions of choline kinase and the hydrolysis of phosphorylcholine to choline and inorganic phosphate.

The observed equilibrium constants (Kobs) of the P-choline hydrolysis reaction have been determined under physiological conditions of temperature (38 degrees) and ionic strength (0.25 M) and physiological ranges of pH and free [Mg2+]. Using sigma and square brackets to indicate total concentrations: (see article.) The value of Kobs has been found to be relatively insensitive to variations in pH and free [Mg2+]. At pH 7.0 and taking the standard state of liquid water to have unit activity ([H2O] = 1), Kobs = 26.6 M at free [Mg2+] = 0 [epsilon G0obs = -2.03 kcal/mol(-8.48 kJ/mol)], 26.8 M at free [Mg2+] = 10(-3) M, and 28.4 M at free [Mg2+] = 10(-2) M. At pH 8.0, Kobs = 18.8 M at free [Mg2+] = 0, 19.2 M at free [Mg2+] = 10(-3), and 22.2 M at free [Mg2+] = 10(-2) M. These values apply only to situations where choline and Pi concentrations are both relatively low (such as the conditions found in most tissues). At higher concentrations of phosphate and choline, the value of Kobs becomes significantly increased since HPO42- complexes choline weakly (association constant = 3.3 M-1). The value of K at 38 degrees and I = 0.25 M is calculated to be 16.4 +/- 0.3 M [epsilonG0 = 1.73 kcal/mol (-7.23 kJ/mol)]. The K for the P-choline hydrolysis reaction has been combined with the K for the ATP hydrolysis reaction determined previously under physiological conditions to calculate a value of 4.95 X 10(-3 M [deltaG0 j.28 kcal/mol (13.7 kJ/mol] for the K of the choline kinase reaction (EC 2.7.1.32), an important step in phospholipid metabolism: (see article.) Likewise, values for Kobs for the choline kinase reaction at 38 degrees, pH 7.0, and I = 0.25 M have been calculated to be 5.76 X 10(4) [deltaG0OBS = -6.77 KCAL/MOL (-28.3 KJ/mol)] at [Mg2+] = 0; 1.24 X 10(4) [deltaG0obs = -5.82 kcal/mol (-24.4 kJ/mol)] at [Mg2+] = 10(-3) M and 8.05 X 10(3) [delta G0obs = -5.56 kcal/mol (-23.3 kJ/mol)] at [Mg2+ = 10(-2) M. Attempts to determine the Kobs of the choline kinase reaction directly were unsuccessful because of the high value of the constant. The results indicate that in contrast to the high deltaG0obs for the hydrolysis of the ester bond of acetylcholine, the deltaG0obs for the hydrolysis of the ester bond of P-choline is quite low, among the lowest known for phosphate ester bonds of biological interest.     

Heat capacity changes for protein-peptide interactions in the ribonuclease S system.

Two fragments of pancreatic ribonuclease A, a truncated version of S-peptide (residues 1-15) and S-protein (residues 21-124), combine to give a catalytically active complex designated ribonuclease S. We have substituted the wild-type residue Met-13 with six other hydrophobic residues ranging in size from alanine to phenylalanine and have determined the thermodynamic parameters associated with binding of these analogues to S-protein by titration calorimetry in the temperature range 5-25 degrees C. The heat capacity change (delta Cp) associated with binding was obtained from a global analysis of the temperature dependences of the free energies and enthalpies of binding. The delta Cp's were not correlated in any simple fashion with the nonpolar surface area (delta Anp) buried upon binding.     

Interactions of cationic ligands and proteins with small nucleic acids: analytic treatment of the large coulombic end effect on binding free energy as a function of salt concentration

For nonspecific binding of oligopeptides and other cationic ligands, including proteins, to nucleic acid oligomers, we develop a model capable of quantifying and predicting the salt concentration dependence of the binding free energy (deltaG(o)obs) by way of an analytic treatment of the Coulombic end effect (CEE). Ligands, nucleic acids, and their complexes (species j of valence Zj) are modeled as finite lattices with absolute value(Zj) charged residues; the CEE is quantified by its characteristic length Ne (specified in charged residues) and its consequences for the free energy and ion association of the oligomer. Expressions are developed for the individual site binding constants Ki as a function of position (site number i) of a bound ligand on a nucleic acid and for the observed binding constant Kobs as an ensemble average of Ki. Analysis of deltaG(o)obs = -RT ln Kobs and Sa Kobs identical with (partial differential ln Kobs)/(partial differential ln a(+/-)) for binding of the oligopeptide KWK6 (ZL = +8) to single-stranded (ss) dT(pdT)(absolute value(ZD) oligomers (dT-mers) where ZD = {-6, -10, -11, -14, -15} in the range 0.1-0.25 M Na+ yields Ne = 9.0 +/- 0.8 residues at each end, demonstrating that both KWK6 and the above dT-mers are sufficiently short so that the CEE extends over the entire molecule. The dependences of Kobs and of Sa Kobs on absolute value(ZD) for a given ZL are determined by the difference between 2Ne and the net number of charged residues Q in the complex (Q identical with absolute value(ZD) - ZL). For Q < 2Ne, characteristic of complexes of KWK6 with this set of dT-mers, the distribution of binding free energies deltaG(o)obs = -RT ln Ki for sites along the DNA oligomer is parabolic, and Kobs and Sa Kobs are strongly dependent on absolute value(ZD). For Q > or = 2Ne, the distribution of binding free energies deltaG(o)obs is trapezoidal, and the dependence of Kobs and Sa Kobs on absolute value(ZD) is weaker. Application of the model to nonspecific binding of human DNA polymerase beta to ssDNA demonstrates the significance of the CEE in determining Kobs and Sa Kobs of binding of a cationic site on a protein to a DNA oligomer.     

Thermodynamics of single-stranded RNA binding to oligolysines containing tryptophan.

The equilibrium binding to the synthetic RNA poly(U) of a series of oligolysines containing one, two, or three tryptophans has been examined as a function of pH, monovalent salt concentration (MX), temperature, and Mg2+. Oligopeptides containing lysine (K) and tryptophan (W) of the type KWKp-NH2 and KWKp-CO2 (p = 1-8), as well as peptides containing additional tryptophans or glycines, were studied by monitoring the quenching of the peptide tryptophan fluorescence upon binding poly(U). Equilibrium association constants, K(obs), and the thermodynamic quantities delta G(o)obs, delta H(o)obs, and delta S(o)obs describing peptide-poly(U) binding were measured as well as their dependences on monovalent salt concentration, temperature, and pH. In all cases, K(obs) decreases significantly with increasing monovalent salt concentration, with (delta log K(obs)/delta log [K+]) = -0.74 (+/- 0.04)z, independent of temperature and salt concentration, where z is the net positive charge on the peptide. The origin of these salt effects is entropic, consistent with the release of counterions from the poly(U) upon formation of the complex. Upon extrapolation to 1 M K+, the value of delta G(o)obs is observed to be near zero for all oligolysines binding to poly(U), supporting the conclusion that these complexes are stabilized at lower salt concentrations due to the increase in entropy accompanying the release of monovalent counterions from the poly(U). Only the net peptide charge appears to influence the thermodynamics of these interactions, since no effects of peptide charge distribution were observed. The binding of poly(U) to the monotryptophan peptides displays interesting behavior as a function of the peptide charge. The extent of tryptophan fluorescence quenching, Qmax, is dependent upon the peptide charge for z less than or equal to +4, and the value of Qmax correlates with z-dependent changes in delta H(o)obs and delta S(o)obs(1 M K+), whereas for z greater than or equal to +4, Qmax, delta H(o)obs, and delta S(o)obs (1 M K+) are constant. The correlation between Qmax and delta H(o)obs and delta S(o)obs(1 M K+) suggests a context (peptide charge)-dependence of the interaction of the peptide tryptophan with poly(U). However the interaction of the peptide tryptophan does not contribute substantially to delta G(o)obs for any of the peptides, independent of z, due to enthalpy-entropy compensations. Each of the tryptophans in multiple Trp-containing peptides appear to bind to poly(U) independently, with delta H(o)Trp = -2.9 +/- 0.7, although delta G(o)Trp is near zero due to enthalpy-entropy compensations.(ABSTRACT TRUNCATED AT 400 WORDS)     

Clustering of large hydrophobes in the hydrophobic core of two-stranded alpha-helical coiled-coils controls protein folding and stability

The de novo design and biophysical characterization of two 60-residue peptides that dimerize to fold as parallel coiled-coils with different hydrophobic core clustering is described. Our goal was to investigate whether designing coiled-coils with identical hydrophobicity but with different hydrophobic clustering of non-polar core residues (each contained 6 Leu, 3 Ile, and 7 Ala residues in the hydrophobic core) would affect helical content and protein stability. The disulfide-bridged P3 and P2 differed dramatically in alpha-helical structure in benign conditions. P3 with three hydrophobic clusters was 98% alpha-helical, whereas P2 was only 65% alpha-helical. The stability profiles of these two analogs were compared, and the enthalpy and heat capacity changes upon denaturation were determined by measuring the temperature dependence by circular dichroism spectroscopy and confirmed by differential scanning calorimetry. The results showed that P3 assembled into a stable alpha-helical two-stranded coiled-coil and exhibited a native protein-like cooperative two-state transition in thermal melting, chemical denaturation, and calorimetry experiments. Although both peptides have identical inherent hydrophobicity (the hydrophobic burial of identical non-polar residues in equivalent heptad coiled-coil positions), we found that the context dependence of an additional hydrophobic cluster dramatically increased stability of P3 (Delta Tm approximately equal to 18 degrees C and Delta[urea](1/2) approximately equal to 1.5 M) as compared with P2. These results suggested that hydrophobic clustering significantly stabilized the coiled-coil structure and may explain how long fibrous proteins like tropomyosin maintain chain integrity while accommodating polar or charged residues in regions of the protein hydrophobic core.     

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.     

Use of liquid hydrocarbon and amide transfer data to estimate contributions to thermodynamic functions of protein folding from the removal of nonpolar and polar surface from water.

This extension of the liquid hydrocarbon model seeks to quantify the thermodynamic contributions to protein stability from the removal of nonpolar and polar surface from water. Thermodynamic data for the transfer of hydrocarbons and organic amides from water to the pure liquid phase are analyzed to obtain contributions to the thermodynamics of folding from the reduction in water-accessible surface area. Although the removal of nonpolar surface makes the dominant contribution to the standard heat capacity change of folding (delta C0fold), here we show that inclusion of the contribution from removal of polar surface allows a quantitative prediction of delta C0fold within the uncertainty of the calorimetrically determined value. Moreover, analysis of the contribution of polar surface area to the enthalpy of transfer of liquid amides provides a means of estimating the contributions from changes in nonpolar and polar surface area as well as other factors to the enthalpy of folding (delta H0fold). In addition to estimates of delta H0fold, this extension of the liquid hydrocarbon model provides a thermodynamic explanation for the observation [Privalov, P. L., & Khechinashvili, N. N. (1974) J. Mol. Biol. 86, 665-684] that the specific enthalpy of folding (cal g-1) of a number of globular proteins converges to a common value at approximately 383 K. Because amounts of nonpolar and polar surface area buried by these proteins upon folding are found to be linear functions of molar mass, estimates of both delta C0fold and delta H0fold may be obtained given only the molar mass of the protein of interest.(ABSTRACT TRUNCATED AT 250 WORDS)     

Protein chemistry at membrane interfaces: non-additivity of electrostatic and hydrophobic interactions

Non-specific binding of proteins and peptides to charged membrane interfaces depends upon the combined contributions of hydrophobic (DeltaG(HPhi)) and electrostatic (DeltaG(ES)) free energies. If these are simply additive, then the observed free energy of binding (DeltaG(obs)) will be given by DeltaG(obs)=DeltaG(HPhi)+DeltaG(ES), where DeltaG(HPhi)=-sigma(NP)A(NP) and DeltaG(ES)=zFphi. In these expressions, A(NP) is the non-polar accessible area, sigma(NP) the non-polar solvation parameter, z the formal peptide valence, F the Faraday constant, and phi the membrane surface potential. But several lines of evidence suggest that hydrophobic and electrostatic binding free energies of proteins at membrane interfaces, such as those associated with cell signaling, are not simply additive. In order to explore this issue systematically, we have determined the interfacial partitioning free energies of variants of indolicidin, a cationic proline-rich antimicrobial peptide. The synthesized variants of the 13 residue peptide covered a wide range of hydrophobic free energies, which allowed us to examine the effect of hydrophobicity on electrostatic binding to membranes formed from mixtures of neutral and anionic lipids. Although DeltaG(obs) was always a linear function of DeltaG(HPhi), the slope depended upon anionic lipid content: the slope was 1.0 for pure, zwitterionic phosphocholine bilayers and 0.3 for pure phosphoglycerol membranes. DeltaG(obs) also varied linearly with surface potential, but the slope was smaller than the expected value, zF. As observed by others, this suggests an effective peptide valence z(eff) that is smaller than the formal valence z. Because of our systematic approach, we were able to establish a useful rule-of-thumb: z(eff) is reduced relative to z by about 20 % for each 3 kcal mol(-1) (1 kcal=4.184 kJ) favorable increase in DeltaG(HPhi). For neutral phosphocholine interfaces, we found that DeltaG(obs) could be predicted with remarkable accuracy using the Wimley-White experiment-based interfacial hydrophobicity scale.     

Thermal and urea-induced unfolding of the marginally stable lac repressor DNA-binding domain: a model system for analysis of solute effects on protein processes

Thermodynamic and structural evidence indicates that the DNA binding domains of lac repressor (lacI) exhibit significant conformational adaptability in operator binding, and that the marginally stable helix-turn-helix (HTH) recognition element is greatly stabilized by operator binding. Here we use circular dichroism at 222 nm to quantify the thermodynamics of the urea- and thermally induced unfolding of the marginally stable lacI HTH. Van't Hoff analysis of the two-state unfolding data, highly accurate because of the large transition breadth and experimental access to the temperature of maximum stability (T(S); 6-10 degrees C), yields standard-state thermodynamic functions (deltaG(o)(obs), deltaH(o)(obs), deltaS(o)(obs), deltaC(o)(P,obs)) over the temperature range 4-40 degrees C and urea concentration range 0 相似文献   

An early transition state for folding of the P4-P6 RNA domain

Tertiary folding of the 160-nt P4-P6 domain of the Tetrahymena group I intron RNA involves burying of substantial surface area, providing a model for the folding of other large RNA domains involved in catalysis. Stopped-flow fluorescence was used to monitor the Mg2+-induced tertiary folding of pyrene-labeled P4-P6. At 35 degrees C with [Mg2+] approximately 10 mM, P4-P6 folds on the tens of milliseconds timescale with k(obs) = 15-31 s(-1). From these values, an activation free energy deltaG(double dagger) of approximately 8-16 kcal/mol is calculated, where the large range for deltaG(double dagger) arises from uncertainty in the pre-exponential factor relating k(obs) and delta G(double dagger). The folding rates of six mutant P4-P6 RNAs were measured and found to be similar to that of the wild-type RNA, in spite of significant thermodynamic destabilization or stabilization. The ratios of the kinetic and thermodynamic free energy changes phi = delta deltaG(double dagger)/delta deltaG(o') are approximately 0, implying a folding transition state in which most of the native-state tertiary contacts are not yet formed (an early folding transition state). The k(obs) depends on the Mg2+ concentration, and the initial slope of k(obs) versus [Mg2+] suggests that only approximately 1 Mg2+ ion is bound in the rate-limiting folding step. This is consistent with an early folding transition state, because folded P4-P6 binds many Mg2+ ions. The observation of a substantial deltaG(double dagger) despite an early folding transition state suggests that a simple two-state folding diagram for Mg2+-induced P4-P6 folding is incomplete. Our kinetic data are some of the first to provide quantitative values for an activation barrier and location of a transition state for tertiary folding of an RNA domain.     

Equilibrium constants of the reactions of choline acetyltransferase, carnitine acetyltransferase, and acetylcholinesterase under physiological conditions.

The observed equilibrium constant (Kobs) for the reaction of choline acetyltransferase (EC 2.3.1.6) has been determined under physiological conditions. Using sigma and square brackets to indicate total concentrations of all ionic species present: (see article). The value of Kobs has been determined to be 12.3 plus or minus 0.6 at 38 degrees, pH 7.0 and ionic strength 0.25 M. The value at 25 degrees is not significantly different, and the constant has been found to be insensitive to variations in ionic strength (0.03 to 0.375 M), pH (6.5 TO 7.5) OR FREE [Mg-2+] (0 to 5 mM). The Kobs of this reaction reflects the difference between the observed standard free energy change (delta G-oobs) for the hydrolysis of acetylcholine and the delta G-oobs for the hydrolysis of acetyl-CoA. Since the delta G-oobs for the hydrolysis of acetyl-CoA has been previously determined to be minus 8.54 kcal/mol (minus 35.75 kJ/mol under the same physiological conditions, the delta G-oobs for the reaction of acetylcholinesterase (EC 3.1.1.7): (SEE ARTICLE). Can be calculated to be minus 6.99 kcal/mol (minus 29.26 kJ/mol) at pH ionic strength 0.25 M and 38 degrees, taking the standard state of liquid water to have unit activity ([H2O] equals 1). The pKa for acetic acid under the same conditions, has been determined to be 4.60 plus or minus 0.01, allowing the Kobs for the pH-independent reaction (see article). To be calculated to be 3.28 times 10-2 M. Choline and carnitine are chemical analogues. The Kobs for the corresponding reaction of carnitine acetyltransferase (EC 2.3.1.7). (SEE ARTICLE). Under the same physiological conditions of pH (7.0), ionic strength (0.25 M), and temperature (38 degrees) has been determined to be 1.73 plus or minus 0.05, making the delta G-oobs for the hydrolysis of acetylcholine only 1.21 kcal/mol (5.06 kJ) less negative than that for the hydrolysis of acetylcarnitine.     

Effects of pH and free Mg2+ on the Keq of the creatine kinase reaction and other phosphate hydrolyses and phosphate transfer reactions.

The observed equilibrium constants (Kobs) of the creatine kinase (EC 2.7.3.2), myokinase (EC 2.7.4.3), glucose-6-phosphatase (EC 3.1.3.9), and fructose-1,6-diphosphatase (EC 3.1.3.11) reactions have been determined at 38 degrees C, pH 7.0, ionic strength 0.25, and varying free magnesium concentrations. The equilibrium constant (KCK) for the creatine kinase reaction defined as: KCK = [sigma ATP] [sigma creatine] divided by ([sigma ADP] [sigma creatine-P] [H+]) was measured at 0.25 ionic strength and 38 degrees C and was shown to vary with free [Mg2+]. The value was found to be 3.78 x 10(8) M-1 at free [Mg2+] = 0 and 1.66 x 10(9) M-1 at free [Mg2+] = 10(-3) M. Therefore, at pH 7.0, the value of Kobs, defined as Kobs = KCK[H+] = [sigma ATP] [sigma creatine] divided by ([sigma ADP] [sigma creatine-P] was 37.8 at free [Mg2+] = 0 and 166 at free [Mg2+] = 10(-3) M. The Kobs value for the myokinase reaction, 2 sigma ADP equilibrium sigma AMP + sigma ATP, was found to vary with free [Mg2+], being 0.391 at free [Mg2+] = 0 and 1.05 at free [Mg2+] = 10(-3) M. Taking the standard state of water to have activity equal to 1, the Kobs of glucose-6-P hydrolysis, sigma glucose-6-P + H2O equilibrium sigma glucose + sigma Pi, was found not to vary with free [Mg2+], being 110 M at both free [Mg2+] = 0 and free [Mg2+] = 10(-3) M. The Kobs of fructose-1,6-P2 hydrolysis, sigma fructose-1,6-P2 equilibrium sigma fructose-6-P + sigma Pi, was found to vary with free [Mg2+], being 272 M at free [Mg2+] = 0 and 174 M at free [Mg2+] = 0.89 x 10(-3) M.     

Surface hydrophobic amino acid residues in cellulase molecules as a structural factor responsible for their high denim-washing performance

The denim-washing performance of six purified fungal cellulases (four endo-1,4-beta-D-glucanases and two cellobiohydrolases) was compared using a model microassay. The performance of cellobiohydrolases per mg of protein was much lower than that of endoglucanases. For endoglucanases, it varied up to 5 times between the best and the worst enzyme. Experiments with amino acids immobilized on cross-linked agarose showed that their side chains may bind indigo owing to hydrophobic interactions and formation of hydrogen bonds. The best binding effects provided Tyr and Phe. Analysis of three-dimensional structures of cellulase molecules showed that a certain correlation exists between the washing performance of enzyme and (i) quantity (percentage) of aromatic residues exposed to solvent on the surface of protein globule or (ii) overall percentage of the surface hydrophobic residues. Data presented provide an evidence that the molecules of certain cellulases, which have hydrophobic domains (clusters of closely located non-polar residues) on their surface, may bind indigo and thus act as emulsifiers helping the dye to float out of cellulose fibers to the bulk solution.     

A comparison of dogfish and bovine chymotrypsins

Bovine and dogfish chymotrypsins were compared to determine if chymotrypsin from a poikilothermic organism (spiny dogfish (Squalus acanthias] adapted to low temperatures possessed catalytic properties different from those of the same enzyme from a warm-blooded animal. An improved procedure was developed for isolating dogfish pancreatic chymotrypsin. The least hydrophobic and smallest substrate used, p-nitrophenyl acetate, had similar enthalpies of association (delta Ha) with both enzymes, whereas larger, more hydrophobic substrates had delta Ha values that were of opposite sign for the two enzymes. As the temperature increased, the association constants (1/Ks) for p-nitrophenyl valerate and p-nitrophenyltrimethyl acetate increased for dogfish chymotrypsin and decreased for bovine chymotrypsin, while the free energies of association (delta Ga) remained relatively constant. Acylation of chymotrypsin was 1.5-2.5 times slower in the dogfish enzyme than in the bovine enzyme except below 15 degrees C with p-nitrophenyltrimethyl acetate. delta H++ for acylation by p-nitrophenyltrimethyl acetate were 2.0 kcal/mol for the dogfish enzyme and 10.2 kcal/mol for the bovine, whereas delta H++ values were only slightly lower in the dogfish enzyme for the other two substrates. For all substrates, the deacylation rate constant (kcat) was greater with dogfish chymotrypsin than bovine. However, the free energies of activation (delta G++) for deacylation were nearly equal between the two enzymes for each of the substrates.     

Regulatory behavior of Escherichia coli aspartate transcarbamylase altered by site-specific mutation of Tyr240----Phe in the catalytic chain

Isotopic exchange kinetics at chemical equilibrium have been used to identify changes in the regulatory properties of aspartate transcarbamylase (ATCase) caused by site-specific mutation of Tyr240----Phe (Y240F) in the catalytic chain. With both wild-type and the mutant enzymes, ATP activates both [14C]Asp in equilibrium N-carbamyl-L-aspartate (C-Asp) and the [32P]carbamyl phosphate (C-P) in equilibrium Pi exchanges. In contrast, with wild-type enzyme, CTP inhibits both exchanges, but with Y240F mutant enzyme CTP inhibits Asp in equilibrium C-Asp exchange and activates C-P in equilibrium Pi exchange. The bisubstrate analog N-(phosphonacetyl-L-aspartate), PALA, activates Asp in equilibrium C-Asp at a lower concentration with the Y240F enzyme, but the extent of activation is decreased, relative to wild-type enzyme. PALA activation of C-P in equilibrium Pi observed with wild-type enzyme disappears completely with the Y240F mutant enzyme. Analysis of perturbations of exchange rates by ATP and CTP were carried out by systematic methods plus computer-based simulations with the ISOBI program. These analyses indicate that (a) ATP increases the rates of association and dissociation for both C-P and Asp, but (b) CTP differentially increases the rate of C-P association to a greater degree than dissociation, but also decreases the rates for Asp association and dissociation in equal proportion. In addition, Arrhenius plots for Y240F ATCase suggest that ATP and CTP act by different mechanisms: ATP increases Vmax (decreases delta G not equal to) uniformly at all temperatures, whereas CTP does not alter either Vmax (delta G not equal to) or the Arrhenius slope (delta H not equal to).     

Proton nuclear magnetic resonance studies on the kinetics of tryptic fragments of calmodulin upon calcium binding

The kinetics of the Ca2+-dependent conformational change of the tryptic fragments F12 (residues 1-75) and F34 (residues 78-148) of calmodulin were studied by 1H-NMR. Resonances of two phenylalanines, 16 (or 19) and 65 (or 68), N epsilon, N epsilon, N epsilon-trimethyllysine-115 and tyrosine-138 were examined by the saturation-transfer technique or computer-aided line-shape simulation to obtain the rate of the conformational exchange between the Ca2+-free form and the Ca2+-bound form. The rates for F12 and F34 in the presence of 0.2 M KCl at 22 degrees C were 300-500 s-1 and 3-10 s-1, respectively. Activation parameters are as follows: Delta H not equal to = 11(+/- 2) kcal X M-1 and delta S not equal to = -9(+/- 5) cal X K-1 X M-1 for F12, and delta H not equal to = 16(+/- 2) kcal X M-1 and delta S not equal to = -2(+/- 5) cal X K-1 X M-1 for F34. These kinetic data for the conformational exchange are in agreement with those of Ca2+ dissociation from the binding sites obtained by 43Ca-NMR and stopped-flow fluorescence studies.     

Grand canonical Monte Carlo molecular and thermodynamic predictions of ion effects on binding of an oligocation (L8+) to the center of DNA oligomers.

Grand canonical Monte Carlo (GCMC) simulations are reported for aqueous solutions containing excess univalent salt (activities a +/- = 1.76-12.3 mM) and one of the following species: an octacationic rod-like ligand, L8+; a B-DNA oligomer with N phosphate charges (8 < or = N < or = 100); or a complex resulting from the binding of L8+ at the center of an N-mer (24 < or = N < or = 250). Simplified models of these multiply charged species are used in the GCMC simulations to predict the fundamental coulombic contributions to the following experimentally relevant properties: 1) the axial distance over which ligand binding affects local counterion concentrations at the surface of the N-mer; 2) the dependence on N of GCMC preferential interaction coefficients, gamma 32MC identical to delta C3/delta C2l a +/-, T, where C3 and C2 are, respectively, the molar concentrations of salt and the multiply charged species (ligand, N-mer or complex); and 3) the dependence on N of SaKobs identical to d in Kobs/d in a +/- = delta (magnitude of ZJ + 2 gamma 32J), where Kobs is the equilibrium concentration quotient for the binding of L8+ to the center of an N-mer and delta denotes the stoichiometric combination of terms, each of which pertains to a reactant or product J having magnitude of ZJ charges. The participation of electrolyte ions in the ligand binding interaction is quantified by the magnitude of SaKobs, which reflects the net (stoichiometrically weighted) difference in the extent of thermodynamic binding of salt ions to the products and reactants. Results obtained here from GCMC simulations yield a picture of the salient molecular consequences of binding a cationic ligand, as well as thermodynamic predictions whose applicability can be tested experimentally. Formation of the central complex is predicted to cause a dramatic reduction in the surface counterion (e.g., Na+) concentration over a region including but extending well beyond the location of the ligand binding site. For binding a cationic ligand, SaKobs is predicted to be negative, indicating net electrolyte ion release in the binding process. At small enough N, -SaKobs is predicted to decrease strongly toward zero with decreasing N. At intermediate N, -SaKobs appears to exceed its limiting value as N-->infinity.     

Contribution of the hydrophobic effect to globular protein stability.

The decrease in conformational stability, delta(delta G), has been measured for 72 aliphatic side-chain mutants from four proteins in which a larger side-chain is replaced by a smaller side-chain so that steric effects are minimal. When these delta(delta G) values are corrected to the same accessibility, namely 100% buried, then the following -delta(delta G) values per -CH2- group (in kcal/mol) are obtained: Ile----Val (1.26), Ala (1.26), Gly (1.26); Leu----Ala (1.16), Gly (1.21); Val----Ala (1.23), Gly (1.53). The average of these values is 1.27(+/- 0.07) kcal/mol. The 72 individual values range from 0 to 2.4 kcal/mol with an average value of 1.27(+/- 0.51) (standard deviation) kcal/mol. When the delta Gtr values from n-octanol to water are corrected for the difference in volume between the solutes and the solvents, the average value for the same substitutions is 1.25(+/- 0.05) kcal/mol. This suggests that proteins gain 1.3(+/- 0.5) kcal/mol in stability for each -CH2- group buried in folding, and, furthermore, that the volume corrected delta Gtr values for n-octanol for the amino acid side-chains provide good estimates of the contribution of the hydrophobic effect to globular protein stability.     

Thermodynamic stoichiometries of participation of water, cations and anions in specific and non-specific binding of lac repressor to DNA. Possible thermodynamic origins of the "glutamate effect" on protein-DNA interactions.

The objective of this study is to quantify the contributions of cations, anions and water to stability and specificity of the interaction of lac repressor (lac R) protein with the strong-binding symmetric lac operator (Osym) DNA site. To this end, binding constants Kobs and their power dependences on univalent salt (MX) concentration (SKobs = d log Kobs/d log[MX]) have been determined for the interactions of lac R with Osym operator and with non-operator DNA using filter binding and DNA cellulose chromatography, respectively. For both specific and non-specific binding of lac R, Kobs at fixed salt concentration [KX] increases when chloride (Cl-) is replaced by the physiological anion glutamate (Glu-). At 0.25 M-KX, the increase in Kobs for Osym is observed to be approximately 40-fold, whereas for non-operator DNA the increase in Kobs is estimated by extrapolation to be approximately 300-fold. For non-operator DNA, SKobsRD is independent of salt concentration within experimental uncertainty, and is similar in KCl (SKobs,RDKCl = -9.8(+/- 1.0) between 0.13 M and 0.18 M-KCl) and KGlu (SKobs,RDKGlu = -9.3(+/- 0.7) between 0.23 M and 0.36 M-KGlu). For Osym DNA, SKobsRO varies significantly with the nature of the anion, and, at least in KGlu appears to decrease in magnitude with increasing [KGlu]. Average magnitudes of SKobsRO are less than SKobsRD, and, for specific binding decrease in the order [SKobsRO,KCl[>[SKobsRO,KAc[>[SKobsRO,KGlu[ . Neither KobsRO nor SKobsRO is affected by the choice of univalent cation M+ (Na+, K+, NH4+, or mixtures thereof, all as the chloride salt), and SKobsRO is independent of [MCl] in the range examined (0.125 to 0.3 M). This behavior of SKobsRO is consistent with that expected for a binding process with a large contribution from the polyelectrolyte effect. However, the lack of an effect of the nature of the cation on the magnitude of KobsRO at a fixed [MX] is somewhat unexpected, in view of the order of preference of cations for the immediate vicinity of DNA (NH4+ > K+ > Na+) observed by 23Na nuclear magnetic resonance. For both specific and non-specific binding, the large stoichiometry of cation release from the DNA polyelectrolyte is the dominant contribution to SKobs. To interpret these data, we propose that Glu- is an inert anion, whereas Ac- and Cl- compete with DNA phosphate groups in binding to lac repressor. A thermodynamic estimate of the minimum stoichiometry of water release from lac repressor and Osym operator (210(+/- 30) H2O) is determined from analysis of the apparently significant reduction in [SKobsRO,KGlu[ with increasing [KGlu] in the range 0.25 to 0.9 M. According to this analysis, SKobs values of specific and non-specific binding in KGlu differ primarily because of the release of water in specific binding. In KAc and KCl, we deduce that anion competition affects Kobs and SKobs to an extent which differs for different anions and for the different binding modes.