Structural factors controlling ligand binding to myoglobin: A kinetic hole-burning study

Using temperature-derivative spectroscopy in the temperature range below 100 K, we have studied the dependence of the Soret band on the recombination barrier in sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation at 12 K. The spectra were separated into contributions from the photodissociated species, Mb*CO, and CO-bound myoglobin. The line shapes of the Soret bands of both photolyzed and liganded myoglobin were analyzed with a model that takes into account the homogeneous bandwidth, coupling of the electronic transition to vibrational modes, and static conformational heterogeneity. The analysis yields correlations between the activation enthalpy for rebinding and the model parameters that characterize the homogeneous subensembles within the conformationally heterogeneous ensemble. Such couplings between spectral and functional parameters arise when they both originate from a common structural coordinate. This effect is frequently denoted as “kinetic hole burning.” The study of these correlations gives direct insights into the structure–function relationship in proteins. On the basis of earlier work that assigned spectral parameters to geometric properties of the heme, the connections with the heme geometry are discussed. We show that two separate structural coordinates influence the Soret line shape, but only one of the two is coupled to the enthalpy barrier for rebinding. We give evidence that this coordinate, contrary to widespread belief, is not the iron displacement from the mean heme plane.

Stability and dynamics in a hyperthermophilic protein with melting temperature close to 200°C

The rubredoxin protein from the hyperthermophilic archaebacterium Pyrococcus furiosus was examined by a hydrogen exchange method. Even though the protein does not exhibit reversible thermal unfolding, one can determine its stability parameters—free energy, enthalpy, entropy, and melting temperature—and also the distribution of stability throughout the protein, by using hydrogen exchange to measure the reversible cycling of the protein between native and unfolded states that occurs even under native conditions.

Energetics of the interaction between water and the helical peptide group and its role in determining helix propensities

The alanine helix provides a model system for studying the energetics of interaction between water and the helical peptide group, a possible major factor in the energetics of protein folding. Helix formation is enthalpy-driven (−1.0 kcal/mol per residue). Experimental transfer data (vapor phase to aqueous) for amides give the enthalpy of interaction with water of the amide group as ≈−11.5 kcal/mol. The enthalpy of the helical peptide hydrogen bond, computed for the gas phase by quantum mechanics, is −4.9 kcal/mol. These numbers give an enthalpy deficit for helix formation of −7.6 kcal/mol. To study this problem, we calculate the electrostatic solvation free energy (ESF) of the peptide groups in the helical and β-strand conformations, by using the delphi program and parse parameter set. Experimental data show that the ESF values of amides are almost entirely enthalpic. Two key results are: in the β-strand conformation, the ESF value of an interior alanine peptide group is −7.9 kcal/mol, substantially less than that of N-methylacetamide (−12.2 kcal/mol), and the helical peptide group is solvated with an ESF of −2.5 kcal/mol. These results reduce the enthalpy deficit to −1.5 kcal/mol, and desolvation of peptide groups through partial burial in the random coil may account for the remainder. Mutant peptides in the helical conformation show ESF differences among nonpolar amino acids that are comparable to observed helix propensity differences, but the ESF differences in the random coil conformation still must be subtracted.

Binding of the cellulose-binding domain of exoglucanase Cex from Cellulomonas fimi to insoluble microcrystalline cellulose is entropically driven.

Isothermal titration microcalorimetry is combined with solution-depletion isotherm data to analyze the thermodynamics of binding of the cellulose-binding domain (CBD) from the beta-1,4-(exo)glucanase Cex of Cellulomonas fimi to insoluble bacterial microcrystalline cellulose. Analysis of isothermal titration microcalorimetry data against two putative binding models indicates that the bacterial microcrystalline cellulose surface presents two independent classes of binding sites, with the predominant high-affinity site being characterized by a Langmuir-type Ka of 6.3 (+/-1.4) x 10(7) M-1 and the low-affinity site by a Ka of 1.1 (+/-0.6) x 10(6) M-1. CBDCex binding to either site is exothermic, but is mainly driven by a large positive change in entropy. This differs from protein binding to soluble carbohydrates, which is usually driven by a relatively large exothermic standard enthalpy change for binding. Differential heat capacity changes are large and negative, indicating that sorbent and protein dehydration effects make a dominant contribution to the driving force for binding.

Mimicry of the calcium-induced conformational state of troponin C by low temperature under pressure.

Calcium binding to the N-domain of troponin C initiates a series of conformational changes that lead to muscle contraction. Calcium binding provides the free energy for a hydrophobic region in the core of N-domain to assume a more open configuration. Fluorescence measurements on a tryptophan mutant (F29W) show that a similar conformational change occurs in the absence of Ca2+ when the temperature is lowered under pressure. The conformation induced by subzero temperatures binds the hydrophobic probe bis-aminonaphthalene sulfonate, and the tryptophan has the same fluorescence lifetime (7 ns) as in the Ca2+-bound form. The decrease in volume (delta V = -25.4 ml/mol) corresponds to an increase in surface area. Thermodynamic measurements suggest an enthalpy-driven conformational change that leads to an intermediate with an exposed N-domain core and a high affinity for Ca2+.

Binding of a hairpin polyamide in the minor groove of DNA: sequence-specific enthalpic discrimination.

Hairpin polyamides are synthetic ligands for sequence-specific recognition in the minor groove of double-helical DNA. A thermodynamic characterization of the DNA-binding properties exhibited by a six-ring hairpin polyamide, ImPyPy-gamma-PyPyPy-beta-Dp (where Im = imidazole, Py = pyrrole, gamma = gamma-aminobutyric acid, beta = beta-alanine, and Dp = dimethylaminopropylamide), reveals an approximately 1-2 kcal/mol greater affinity for the designated match site, 5'-TGTTA-3', relative to the single base pair mismatch sites, 5'-TGGTA-3' and 5'-TATTA-3'. The enthalpy and entropy data at 20 degrees C reveal this sequence specificity to be entirely enthalpic in origin. Correlations between the thermodynamic driving forces underlying the sequence specificity exhibited by ImPyPy-gamma-PyPyPy-beta-Dp and the structural properties of the heterodimeric complex of PyPyPy and ImPyPy bound to the minor groove of DNA provide insight into the molecular forces that govern the affinity and specificity of pyrrole-imidazole polyamides.

Persistent confusion of total entropy and chemical system entropy in chemical thermodynamics.

The change in free energy with temperature at constant pressure of a chemical reaction is determined by the sum (dS) of changes in entropy of the system of reagents, dS(i), and the additional entropy change of the surroundings, dS(H), that results from the enthalpy change, W. A faulty identification of the total entropy change on reaction with dS(i) has been responsible for the attribution of general validity to the expressions (d deltaG/dT)p = -deltaS(i) and d(deltaG/T)/d(1/T)= deltaH, which are found in most textbooks and in innumerable papers.

Sequence-specific DNA-binding dominated by dehydration.

Fluorescence spectroscopy and isothermal titration calorimetry were used to study the thermodynamics of binding of the glucocorticoid receptor DNA-binding domain to four different, but similar, DNA-binding sites. The binding sites are two naturally occurring sites that differ in the composition of one base pair, i.e., an A-T to G-C mutation, and two sites containing chemical intermediates of these base pairs. The calorimetrically determined heat capacity change (Delta C(p)o(obs)) for glucocorticoid receptor DNA-binding domain binding agrees with that calculated for dehydration of solvent-accessible surface areas. A dominating effect of dehydration or solvent reorganization on the thermodynamics is also consistent with an observed linear relationship between observed enthalpy change (Delta Ho(obs)) and observed entropy change (Delta So(obs)) with a slope close to the experimental temperature. Comparisons with structural data allow us to rationalize individual differences between Delta Ho(obs) (and Delta So(obs)) for the four complexes. For instance, we find that the removal of a methyl group at the DNA-protein interface is enthalpically favorable but entropically unfavorable, which is consistent with a replacement by an ordered water molecule.

Prediction of the stability of DNA triplexes.

We present rules that allow one to predict the stability of DNA pyrimidine.purine.pyrimidine (Y.R.Y) triple helices on the basis of the sequence. The rules were derived from van't Hoff analysis of 23 oligonucleotide triplexes tested at a variety of pH values. To predict the enthalpy of triplex formation (delta H degrees), a simple nearest-neighbor model was found to be sufficient. However, to accurately predict the free energy of the triplex (delta G degrees), a combination model consisting of five parameters was needed. These parameters were (i) the delta G degrees for helix initiation, (ii) the delta G degrees for adding a T-A.T triple, (iii) the delta G degrees for adding a C(+)-G.C triple, (iv) the penalty for adjacent C bases, and (v) the pH dependence of the C(+)-G.C triple's stability. The fitted parameters are highly consistent with thermodynamic data from the basis set, generally predicting both delta H degrees and delta G degrees to within the experimental error. Examination of the parameters points out several interesting features. The combination model predicts that C(+) -G.C. triples are much more stabilizing than T-A.T triples below pH 7.0 and that the stability of the former increases approximately equal to 1 kcal/mol per pH unit as the pH is decreased. Surprisingly though, the most stable sequence is predicted to be a CT repeat, as adjacent C bases partially cancel the stability of one another. The parameters successfully predict tm values from other laboratories, with some interesting exceptions.

Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

Thermal unfolding of human high-density apolipoprotein A-1: implications for a lipid-free molten globular state.

Apolipoprotein A-1 (apoA-1) in complex with high-density lipoprotein is critically involved in the transport and metabolism of cholesterol and in the pathogenesis of atherosclerosis. We reexamined the thermal unfolding of lipid-free apoA-1 in low-salt solution at pH approximately 7, by using differential scanning calorimetry and circular dichroism. At protein concentrations <5 mg/ml, thermal unfolding of apoA-1 is resolved as an extended peak (25 degrees C-90 degrees C) that can be largely accounted for by a single reversible non-two-state transition with midpoint Tm 57 +/- 1 degree C, calorimetric enthalpy deltaH(Tm)= 200 +/- 20 kcal/mol (1 kcal = 4.18 kJ), van't Hoff enthalpy deltaHv(Tm) approximately 32.5 kcal/mol, and cooperativity deltaHv(Tm)/deltaH(Tm) approximately 0.16. The enthalpy deltaH(Tm) can be accounted for by melting of the alpha-helical structure that is inferred by CD to constitute approximately 60% of apoA-1 amino acids. Farand near-UV CD spectra reveal noncoincident melting of the secondary and tertiary structural elements and indicate a well-defined secondary structure but a largely melted tertiary structure for apoA-1 at approximately 37 degrees C and pH 7. This suggests a molten globular-like state for lipid-free apoA-1 under near-physiological conditions. Our results suggest that in vivo lipid binding by apoA-1 may be mediated via the molten globular apolipoprotein state in plasma.

Temperature effects in hydrophobic interaction chromatography.

The effect of temperature from 5 degrees C to 50 degrees C on the retention of dansyl derivatives of amino acids in hydrophobic interaction chromatography (HIC) was investigated by HPLC on three stationary phases. Plots of the logarithmic retention factor against the reciprocal temperature in a wide range were nonlinear, indicative of a large negative heat capacity change associated with retention. By using Kirchoff's relations, the enthalpy, entropy, and heat capacity changes were evaluated from the logarithmic retention factor at various temperatures by fitting the data to a logarithmic equation and a quadratic equation that are based on the invariance and on an inverse square dependence of the heat capacity on temperature, respectively. In the experimental temperature interval, the heat capacity change was found to increase with temperature and could be approximated by the arithmetic average. For HIC retention of a set of dansylamino acids, both enthalpy and entropy changes were positive at low temperatures but negative at high temperatures as described in the literature for other processes based on the hydrophobic effect. The approach presented here shows that chromatographic measurements can be not only a useful adjunct to calorimetry but also an alternative means for the evaluation of thermodynamic parameters.

Least activation path for protein folding: investigation of staphylococcal nuclease folding by stopped-flow circular dichroism.

Is the pathway of protein folding determined by the relative stability of folding intermediates, or by the relative height of the activation barriers leading to these intermediates? This is a fundamental question for resolving the Levinthal paradox, which stated that protein folding by a random search mechanism would require a time too long to be plausible. To answer this question, we have studied the guanidinium chloride (GdmCl)-induced folding/unfolding of staphylococcal nuclease [(SNase, formerly EC 3.1.4.7; now called microbial nuclease or endonuclease, EC 3.1.31.1] by stopped-flow circular dichroism (CD) and differential scanning microcalorimetry (DSC). The data show that while the equilibrium transition is a quasi-two-state process, kinetics in the 2-ms to 500-s time range are triphasic. Data support the sequential mechanism for SNase folding: U3 <--> U2 <--> U1 <--> N0, where U1, U2, and U3 are substates of the unfolded protein and N0 is the native state. Analysis of the relative population of the U1, U2, and U3 species in 2.0 M GdmCl gives delta-G values for the U3 --> U2 reaction of +0.1 kcal/mol and for the U2 --> U1 reaction of -0.49 kcal/mol. The delta-G value for the U1 --> N0 reaction is calculated to be -4.5 kcal/mol from DSC data. The activation energy, enthalpy, and entropy for each kinetic step are also determined. These results allow us to make the following four conclusions. (i) Although the U1, U2, and U3 states are nearly isoenergetic, no random walk occurs among them during the folding. The pathway of folding is unique and sequential. In other words, the relative stability of the folding intermediates does not dictate the folding pathway. Instead, the folding is a descent toward the global free-energy minimum of the native state via the least activation path in the vast energy landscape. Barrier avoidance leads the way, and barrier height limits the rate. Thus, the Levinthal paradox is not applicable to the protein-folding problem. (ii) The main folding reaction (U1 --> N0), in which the peptide chain acquires most of its free energy (via van der Waals' contacts, hydrogen bonding, and electrostatic interactions), is a highly concerted process. These energy-acquiring events take place in a single kinetic phase. (iii) U1 appears to be a compact unfolded species; the rate of conversion of U2 to U1 depends on the viscosity of solution. (iv) All four relaxation times reported here depend on GdmCl concentrations: it is likely that none involve the cis/trans isomerization of prolines. Finally, a mechanism is presented in which formation of sheet-like chain conformations and a hydrophobic condensation event precede the main-chain folding reaction.

Negative activation enthalpies in the kinetics of protein folding.

Although the rates of chemical reactions become faster with increasing temperature, the converse may be observed with protein-folding reactions. The rate constant for folding initially increases with temperature, goes through a maximum, and then decreases. The activation enthalpy is thus highly temperature dependent because of a large change in specific heat (delta Cp). Such a delta Cp term is usually presumed to be a consequence of a large decrease in exposure of hydrophobic surfaces to water as the reaction proceeds from the denatured state to the transition state for folding: the hydrophobic side chains are surrounded by "icebergs" of water that melt with increasing temperature, thus making a large contribution to the Cp of the denatured state and a smaller one to the more compact transition state. The rate could also be affected by temperature-induced changes in the conformational population of the ground state: the heat required for the progressive melting of residual structure in the denatured state will contribute to delta Cp. By examining two proteins with different refolding mechanisms, we are able to find both of these two processes; barley chymotrypsin inhibitor 2, which refolds from a highly unfolded state, fits well to a hydrophobic interaction model with a constant delta Cp of activation, whereas barnase, which refolds from a more structured denatured state, deviates from this ideal behavior.

Significant discrepancies between van't Hoff and calorimetric enthalpies.

In this paper we show that the usual assumption in studies of the temperature variation of equilibrium constants for equilibria of the form A+B <-->AB that a plot of ln K vs. 1/T (K = equilibrium constant, T = temperature in degrees kelvin) is a straight line with slope equal to -delta HvH/R (delta HvH = van't Hoff or apparent enthalpy, R = gas constant) is not valid in many cases. In all the cases considered here, delta HvH is temperature dependent and is significantly different from the true or calorimetrically measured enthalpy, and the respective values for delta Cp are also significantly different.