Convex Arrhenius plots and their interpretation Donald G. Truhlar*† and Amnon Kohen†‡ *Department of Chemistry and Supercomputer Institute, University of Minnesota, Minneapolis, MN 55455-0431; and ‡Department of Chemistry, University of Iowa, Iowa City, IA 52242-1294 Edited by Rudolph A. Marcus, California Institute of Technology, Pasadena, CA, and approved November 16, 2000 (received for review September 12, 2000)
This paper draws attention to selected experiments on enzymecatalyzed reactions that show convex Arrhenius plots, which are very rare, and points out that Tolman’s interpretation of the activation energy places a fundamental model-independent constraint on any detailed explanation of these reactions. The analysis presented here shows that in such systems, the rate coefficient as a function of energy is not just increasing more slowly than expected, it is actually decreasing. This interpretation of the data provides a constraint on proposed microscopic models, i.e., it requires that any successful model of a reaction with a convex Arrhenius plot should be consistent with the microcanonical rate coefficient being a decreasing function of energy. The implications and limitations of this analysis to interpreting enzyme mechanisms are discussed. This model-independent conclusion has broad applicability to all fields of kinetics, and we also draw attention to an analogy with diffusion in metastable fluids and glasses.
I
n recent years, an increasing number of studies have measured the temperature dependence of enzyme reactions for both thermostable and mesostable enzymes (enzymes from organisms that grow optimally at high temperature and at normal body temperature, respectively). The temperature dependence is usually exhibited in an Arrhenius plot, which is the natural logarithm of the rate coefficient vs. reciprocal of the absolute temperature (Fig. 1). Some of the methodologies for measuring the temperature dependences of reaction rates have been applied to measure the temperature dependences of the kinetic isotope effects (KIEs) of the catalyzed reactions as well, and such experiments provide experimental data very relevant to characterizing the dynamical bottlenecks of the reactions (1–6). A large number of theoretical models have been advanced for the interpretation of this kind of data (7–19). Experiments on reactions catalyzed by a variety of thermostable and mesostable dehydrogenase and oxidase enzymes have shown that in some cases the Arrhenius plot for the chemical step is convex (20–23). In the present paper, we comment on the significance of this observation, and we provide an interpretation of this finding that places it in the broader perspective of nonbiological chemical as well as biochemical kinetics. This paper points out that the convex Arrhenius behavior places a severe constraint on the microcanonical unimolecular rate coefficient for the chemical step. In particular, we can go beyond the obvious interpretation that a convex Arrhenius plot can be interpreted as a lower than expected rate at high temperatures, and we can say, under certain conditions, that the microcanonical unimolecular rate coefficient for the chemical step must actually decrease with increasing energy. In general, it is desirable to separate the interpretation of experimental data into two steps. First, we infer general consequences that rest on solid foundations such as thermodynamics and general principles of statistical mechanics. Second, we proceed to more microscopically detailed models that rest on less sure footing. The new aspect of the present analysis is concerned entirely with the first step. The striking feature of the reactions in refs. 20–23 that we discuss here is the positive convexity C of their Arrhenius plots, where C ⬅ ⫺d2共ln k共T兲兲兾d2共1兾T兲, 848 – 851 兩 PNAS 兩 January 30, 2001 兩 vol. 98 兩 no. 3
[1]
k(T) is the catalytic rate coefficient, and T is temperature. We write T2 dEa R dT
[2]
d ln k共T兲 , d共1兾T兲
[3]
C⫽⫺ where E a ⬅ ⫺R
E a is the temperature-dependent phenomenological energy of activation, and R is the gas constant. Thus a positive convexity means that E a decreases with increasing temperature. In the broad field of kinetics, not restricting consideration to enzyme kinetics, when nonlinear Arrhenius plots can be observed, they are almost always concave, i.e., C is negative. This general result is explained by Tolman’s interpretation of the activation energy, which may be written as (25–27) ញ ⫺E , Ea ⫽ E
[4]
where a double overbar denotes an average over all reacting systems, and a single overbar denotes an average over all systems, whether or not they react. In words, Tolman’s Eq. 4 states that the negative of the local slope of an Arrhenius plot can be interpreted as equal to the average energy of molecules that react minus the average energy of all possible reactants. For unimolecular reactions of a thermally equilibrated species: ⫽ E
冕
⬁
dE
⫺⬁
dP E dE
[5]
and
冕 冕
⬁
ញ ⫽ E
dE
⫺⬁ ⬁
⫺⬁
dP k共E兲E dE
dP dE k共E兲 dE
,
[6]
where E is total energy, dP兾dE is the normalized probability density that a system has energy E in a canonical ensemble at temperature T, and k(E) is the unimolecular rate coefficient for ញ increases a microcanonical ensemble at energy E. § Usually, E with increasing T because k(E) is usually an more rapidly than E This paper was submitted directly (Track II) to the PNAS office. Abbreviation: KIE, Kinetic isotope effect. †E-mail:
[email protected] and
[email protected].
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be hereby marked “advertisement” in accordance with 18 U.S.C. §1734 solely to indicate this fact. §Recall
that a canonical ensemble is one characterized by a temperature, and a microcanonical ensemble is a slice out of a canonical ensemble such that the slice contains all
This result presents a constraint on any microscopic interpretation, and—equally importantly—it suggests new microscopic interpretations. For example, suppose that the reactant configuration space may be partitioned into two types of conformations, one of which (called type R) is reactive with rate coefficient k R(E) and the other of which (called type N) is nonreactive with zero rate coefficient. Then
increasing function of energy.¶ This explains why E a usually increases with T and hence why C is negative. A convex ញ decreases or that it increases less Arrhenius plot means that E as temperature increases, and in the next pararapidly than E graph we argue that this ordinarily implies that k(E) decreases with increasing energy over the region of energy that dominates the integrals in Eq. 6. One can make a more precise statement about the conseញ decreasing or increasing less rapidly than E as quence of E temperature increases. By noting that dP ⫽ e⫺E共E兲, dE
[7]
where  is 1兾RT, and (E) is the density of states, and by using Eqs. 4–6, one can easily show that 1 ញ2 2 dEa ញ 兲 ⫺ 共E ⫺ E 兲 兴. 关共E ⫺ E ⫽ dT RT2
[8]
Thus the statement that dE a兾dT is negative is the same as the statement that the distribution of reaction energies is narrower than the distribution of reactant energies. If dP兾dE were sym , then the sign of dE a兾dT would depend on the metric about E shape of the function k(E), and not just its sign. But dP兾dE is actually a highly skewed function. For typical molecules, its shape is well approximated by setting (E) in Eq. 7 equal to E n⫺1, where n is approximately equal to the number of degrees of freedom (31). If one assumes this form and also assumes that k(E) varies slowly and linearly over the important range of the integrals in Eq. 6, one can show that the right-hand side of Eq. 8 is proportional to dk兾dE with a positive constant of proportionality. For more complicated k(E) functions, one no longer has exact proportionality, but the direction of the effect is established by the positive skewedness of dP兾dE. Although the identity of the signs of dE a兾dT and dk兾dE is not a mathematical theorem, it is expected to hold for normal densities of states and smooth k(E) functions. systems with a given total energy. Note that Ea, Eញ , E , and dP兾dE are all associated with the canonical ensemble and therefore are functions of T, and dP兾dE equals the Boltzmannweighted density of states times a normalization constant (28). ¶Microcanonical
rate coefficients usually increase with energy in both over-barrier and tunneling regimes except at very low temperature where the reaction occurs entirely by tunneling in the reactant ground state (29, 30).
Truhlar and Kohen
where P R(E) is the conditional probability that a system having energy E is in a conformation of type R. This provides a mechanism for k(E) to decrease with E, i.e., P R(E) decreases with E faster than k R(E) increases. This could occur because higher-energy systems visit a wider region of phase space and hence spend a smaller fraction of their time in any single constrained region. This interpretation may be of special importance for enzymology. In recent years, several studies have been concerned with understanding the role of conformational sampling in enzyme catalysis (17, 18, 22, 32–37). In a recent review, Bruice and Benkovic (17) argued that the active configurations sampled by the enzymatic complex are not of higher energy than the general population of conformers and that populating these conformations may determine the reaction rate. In contrast to the sometimes detailed models in the papers cited, the present interpretation proceeds on very general grounds without invoking any particular microscopic molecular mechanism or model such as tunneling promoted by thermally excited enzyme fluctuations (11, 13, 38, 39), free energy of activation calculations including proper sample over conformations (40), or specific mechanisms of transitions between microstates of varying conformational flexibility (41–44). Nevertheless, it can clearly accommodate more specific interpretations such as a protein that is rigid at low T or a dynamical motion that requires population of specific low-energy vibrational states. However, nothing we have said here gives specific support to such models, whose validity should be examined by further experiments and detailed simulations. The present discussion does not specifically invoke effects associated with the breakdown of transition state theory. Quantum effects like tunneling can be accommodated by simply including them in k(E) or k R(E). In the context of our analysis, the results mentioned above for thermophilic alcohol dehydrogenase ADH-hT (22, 32) can be inញ and E when terpreted in terms of the relative trends in E temperature is changed. A widely advanced hypothesis is that thermophilic enzymes have evolved to function in their physiological niche, that is, at elevated temperature, and are less active at room temperature (even less active, in some cases, than anticipated by their enthalpy of activation at physiological temperature). Our analysis may offer a connection between that notion and calculations concerned with the identification of active and nonactive conformations (16, 17). For example, Wrba et al. (21) reported an empirical observation that, for enzymes with different melting temperatures, the difference in the temperatures of maximum curvature in Arrhenius plots is similar to the difference in melting temperature. This may be related to shifts in P R(E) and k R(E) dependencies on E in the ranges that are important at the transition temperatures for denaturation and activity. As an additional example, we note that it was recently suggested that it is ‘‘perhaps not surprising that, for COH-bond cleavage reactions, a partitioning of thermal energy into low-frequency protein modes is a major factor. These modes can become significantly excited at ambient temperature, in contrast to the COH stretching mode, which is expected to be largely confined to its zero-point energy level’’ (22). These kinds of considerations can be related to the temperature dependence PNAS 兩 January 30, 2001 兩 vol. 98 兩 no. 3 兩 849
BIOCHEMISTRY
Fig. 1. Concave, linear, and convex Arrhenius plots where k is rate coefficient, and T is the absolute temperature.
[9]
CHEMISTRY
k共E兲 ⫽ P R共E兲kR共E兲,
ញ and E and hence to the concavity or convexity of Arrhenius of E plots. The general conclusion about the energy dependence of the microcanonical rate coefficient is also applicable to convex Arrhenius plots that are occasionally observed in other studies. If one is to apply these arguments, one must be careful to ensure that the convex Arrhenius plot represents a true elementary reaction with internally equilibrated reactants. For example, in the prototype steady-state problem k1 k2 A| L; B O ¡ C, k ⫺1
[10]
d关C兴 ⫽ k关A兴 dt
[11]
where
with k⫽
k 2k 1 , k 2 ⫹ k ⫺1
[12]
it can be applied to k 1, k ⫺1, or k 2, but not generally to k. However, if k ⫺1 and k 1 are much larger than k 2, then k becomes an elementary rate coefficient for A 3 C, and no conclusions are changed by the pre-equilibrium before the rate-determining step. This further applies to the case in which k 2 is much greater than k 1 and k ⫺1, where k becomes k 1. A special case where the steady-state approach may be applied is that in which A corresponds to unbound reactant and enzyme, state B is the Michaelis complex, and state C is the product (45). Note that in this case k 1 (in Eq. 12) should be a second order rate constant leading to the Michaelis equation: rate ⫽
Vmax关S兴 , KM ⫹ 关S兴
[13]
where [S] is the reactant (substrate) concentration, V max is the reaction rate at [S] much higher than K M, which is the Michaelis constant, (k 2 ⫹ k ⫺1)兾k 1, for Eq. 10. The discussion does not apply to effective rate coefficients resulting from a change in the rate-determining step or to changes in effective enzyme concentration [irreversible denaturation or a reversible inhibition equilibrium of the reactants (46)], nor does it apply to rate coefficients with hysteresis attributable to glass formation. This point should be emphasized for enzymatic reactions: do the rates measured actually represent a meaningful consistent rate constant over the whole temperature range? For example, in an early report of convex Arrhenius plot for the thermophilic alcohol dehydrogenase from Bacillus stearothermophilus (ADH-hT), all measurements were conducted with the same substrate concentration (that was saturating the enzyme at room temperature), and V max was presented (47). A later report (22) found that K M increases with temperature so at high temperature [S] was not sufficiently larger than K M, and the measured rate was found to be slower than V max, presenting an apparently convex Arrhenius plot. Even for the simplest realistic enzyme kinetics, mechanism 10 should be augmented by a product release step, and the data should be analyzed with due regard for the commitment to catalysis and the ratio of catalysis (48, 49); otherwise, the measured k cat may show a different temperature dependence than the intrinsic k 2. Refs. 20–24 have been selected as examples where particular care was taken to ensure and demonstrate that the Arrhenius plot corresponds to a single rate-limiting elementary chemical reaction step as function of temperature. Rubach and Plapp (personal communication of unpublished results), for 850 兩 www.pnas.org
example, used stopped-flow pre-steady state techniques to measure KIEs on microscopic rate constants to study the reaction catalyzed by horse liver alcohol dehydrogenase. They found convex Arrhenius plots for both hydride and deuteride transfer, but no temperature dependence of the KIEs, and they demonstrated that no change in the rate-determining step was likely to affect their results. We further excluded cases for which the temperature range allows for irreversibly melted or partially melted protein. In contrast to the enzymatic examples showing convex Arrhenius behavior involving a single elementary step, we mention the case of skin permeation of drugs, where convex Arrhenius behavior has been observed, but where the elementary steps have not been or cannot be sorted out (50, 51). For comparison with permeation through skin, we also mention that the somewhat better defined case of non-Arrhenius temperature dependence of diffusion in amorphous solids has been modeled by a multistep site-hopping model (52, 53). An example from organic chemistry provides a useful illustration of these concepts in conjunction with a temperaturedependent structural reorganization. Limbach has observed convex Arrhenius plots (personal communication of unpublished results) for the base-catalyzed intramolecular proton transfer in 1,3-bis-(4-fluorophenyl)-triazene dissolved in ethyl methyl ether [this can be thought of as a model system for enzymatic processes (55)]. The rate coefficient represents an average over different hydrogen-bonded forms of the reactant, especially the free molecule and the complex with base. The proton transfer takes place in the complex. The convex Arrhenius curvature arises from the fact that the molecule changes at high temperature into the noncomplexed form where the barrier to proton transfer is larger (ref. 55 and H. H. Limbach, personal communication). Thus the Tolman interpretation does not strictly apply because the rate coefficient does not represent a single elementary step. The Arrhenius plot of the reaction from either form separately would presumably be quasilinear or concave. If, however, the different hydrogen-bonded forms are considered to be conformers of a single reactant, the Tolman interpretation may be applied, and it provides an illustration of a molecule that reacts more slowly at higher energy because it explores a wider region of phase space that includes essentially nonreactive conformers. Similarly, Massey et al. (23) suggest a possible explanation of their convex Arrhenius plot involving two competing enzymatic forms, each dominating in a different temperature range. This model is one particular example of a model that could be consistent with the model-independent result presented here. For nonbiological reactions we expect that particularly reactive configurations will often be high-energy species. It is possible that convex Arrhenius plots will be more common in biological systems because an enzyme may have specifically evolved such that the lowest-energy states (i.e., the most accessible ones at ambient T) are the reactive ones, thereby ញ even as T is increased. leading to low E Convex Arrhenius plots have also been observed for diffusion in supercooled liquids and glasses. Such cases, like the skin permeation problem mentioned above, may be complicated by multiple parallel pathways, but nevertheless have yielded to detailed analysis. One widely quoted explanation is that molecules in low-temperature liquids diffuse by crossing appreciable potential barriers, whereas at high-temperature, diffusive motions are almost free (24, 56). This model recognizes that as temperature is raised, thermal energies will become larger compared to barrier heights; what is required for a convex Arrhenius plot is that the average energy of molecules making transitions (in this case, diffusive hops from one local minimum to another) rises more slowly than the average energy of all molecules. There are many detailed theories of non-Arrhenius behavior in metastable fluids and glasses that do not explicitly address the temperature dependence of the activation energy in Truhlar and Kohen
We appreciate helpful discussion with Dan Geselter, Jan Jensen, Judith Klinman, and Hans Limbach. This work was supported in part by the National Science Foundation.
Kohen, A. & Klinman, J. P. (1998) Acc. Chem. Res. 31, 397–404. Kohen, A. & Klinman, J. P. (1999) Chem. Biol. 6, R191–R198. Rickert, K. W. & Klinman, J. P. (1999) Biochemistry 38, 12218–12228. Basran, J., Sutcliffe, M. J. & Scrutton, N. S. (1999) Biochemistry 38, 3218–3222. Scrutton, N. S. (1999) Biochem. Soc. Trans. 27, 767–779. Harris, R. J., Meskys, R., Sutcliffe, M. J. & Scrutton, N. S. (2000) Biochemistry 39, 1189–1198. Warshel, A. (1984) Proc. Natl. Acad. Sci. USA 81, 444–448. Warshel, A., Sussman, F. & Hwang, J.-K. (1988) J. Mol. Biol. 201, 139–159. Warshel, A. (1991) Computer Modeling of Chemical Reactions in Enzymes and Solutions (Plenum, New York). Hwang, J.-K. & Warshel, A. (1996) J. Am. Chem. Soc. 118, 11745–11751. Borgis, D. & Hynes, J. T. (1989) in The Enzyme Catalysis Process, eds. Cooper, A., Houben, J. & Chien, L. (Plenum, New York), pp. 293–303. Bruno, W. J. & Bialek, W. (1992) Biophys. J. 63, 689–699. Antoniou, D. & Schwartz, S. D. (1997) Proc. Natl. Acad. Sci. USA 94, 12360–12365. Kuznetsov, A. M. & Ulstrup, J. (1999) Can. J. Chem. 77, 1085–1096. Truhlar, D. G. & Morokuma, K. (1999) Transition State Modeling for Catalysis (Am. Chem. Soc., Washington, DC). Alhambra, C., Corchado, J. C., Sa´nchez, M. L., Gao, J. & Truhlar, D. J. (2000) J. Am. Chem. Soc. 122, 8197–8203. Bruice, T. C. & Benkovic, S. J. (2000) Biochemistry 39, 6267–6273. Radkiewicz, J. L. & Brooks, C. L. (2000) J. Am. Chem. Soc. 122, 225–231. Wilson, E. K. (July 17, 2000) Chem. Eng. News 78, 42–45. Lakatos, S., Halasz, G. & Zavodszky, P. (1978) Biochem. Soc. Trans. 6, 1195–1197. Wrba, A., Schweiger, A., Schultes, V., Jaenicke, R. & Zavodszky, P. (1990) Biochemistry 29, 7584–7592. Kohen, A., Cannio, R., Bartolucci, S. & Klinman, J. P. (1999) Nature (London) 399, 496–499. Massey, V., Curti, B. & Ganther, H. J. (1966) Biol. Chem. 241, 2347–2357. Ullo, J. & Yip, S. (1989) Phys. Rev. A 39, 5877–5886. Tolman, R. C. (1920) J. Am. Chem. Soc. 42, 2506–2528. Tolman, R. C. (1927) Statistical Mechanics with Applications to Physics and Chemistry (Chemical Catalog Co., New York), pp. 260–270. Truhlar, D. G. (1978) J. Chem. Educ. 55, 309–311. Denbigh, K. (1971) The Principles of Chemical Equilibrium (Cambridge Univ. Press, Cambridge, U.K.), p. 342. Skinner, J. L. & Trommsdorff, H. P. (1988) J. Chem. Phys. 89, 897–907.
30. Wonchoba, S. E., Hu, W. P. & Truhlar, D. G. (1995) Phys. Rev. B 51, 9985–10002. 31. Rice, O. K. (1967) Statistical Mechanics, Thermodynamics, and Kinetics (Freeman, San Francisco), p. 565. 32. Kohen, A. & Klinman, J. P. (2000) J. Am. Chem. Soc. 122, 10738–10739. 33. Gregory, R. B. (1995) Protein-Solvent Interactions (Dekker, New York). 34. Scrutton, N. S., Basran, J. & Sutcliffe, M. J. (1999) Eur. J. Biochem. 264, 666–671. 35. Northrop, D. B. (1999) J. Am. Chem. Soc. 121, 3521–3524. 36. Ostermann, A., Waschipkyo, R., Parak, F. G. & Nienhaus, G. U. (2000) Nature (London) 404, 205–208. 37. Liebl, U., Lipowski, G., Ne´grerie, M., Lambry, J. C., Martin, J. L. & Vos, M. H. (1999) Nature (London) 401, 181–184. 38. Antoniou, D. & Schwartz, S. D. (1998) J. Chem. Phys. 108, 3620–3625. 39. Borgis, D. & Hynes, J. T. (1993) Chem. Phys. 170, 315–346. 40. Warshel, A. (1992) Curr. Opin. Struct. Biol. 2, 230–236. 41. Vihinen, M. (1987) Prot. Eng. 1, 477–480. 42. Varley, P. G. & Pain, R. H. (1991) J. Mol. Biol. 220, 531–538. 43. Rasmussen, B. F., Stock, A. M., Ringe, D. & Petsko, G. A. (1992) Nature (London) 357, 423–424. 44. Zavodszky, P., Kardos, J., Svingor, A. & Petsko, G. A. (1998) Proc. Natl. Acad. Sci. USA 95, 7406–7411. 45. Fersht, A. (1998) Structure and Mechanism in Protein Sciences: A Guide to Enzyme Catalysis and Protein Folding (Freeman, New York). 46. Kistiakowsky, G. B. & Lumry, R. (1949) J. Am. Chem. Soc. 71, 2006–2013. 47. Guagliardi, A., Martino, M., Iaccarino, I., Derosa, M., Rossi, M. & Bartolucci, S. (1996) Int. J. Biochem. Cell Biol. 28, 239–246. 48. Northrop, D. B. (1991) in Enzyme Mechanism from Isotope Effects, ed. Cook, P. F. (CRC, Boca Raton, FL), pp. 181–202. 49. Malany, S., Sawai, M., Sikorski, R. S., Seravalli, J., Quinn, D. M., Radi, Z., Taylor, P., Kronman, C., Velan, B. & Shafferman, A. (2000) J. Am. Chem. Soc. 122, 2981–2987. 50. Ohara, N., Takayama, K., Machida, Y. & Nagai, K. (1994) Int. J. Pharm. 105, 31–38. 51. Ohara, N., Takayama, K. & Nagai, T. (1995) Biol. Pharm. Bull. 18, 439–442. 52. Chang, W. & Xu, G. (1993) J. Chem. Phys. 99, 2001–2003. 53. Xu, G. J. (1994) Non-Cryst. Sol. 170, 68–72. 54. Ediger, M. D., Angell, C. A. & Nagel, S. R. (1996) J. Phys. Chem. 100, 13200–13212. 55. Mannle, F. & Limbach, H. H. (1996) Agnew. Chem. Int. Ed. 35, 441–442. 56. Goldstein, M. (1969) J. Chem. Phys. 51, 3728–3739.
1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20. 21. 22. 23. 24. 25. 26. 27. 28.
CHEMISTRY
29.
BIOCHEMISTRY
chemistry, in cluster reactions, etc.) where convex Arrhenius plots may be observed, and we have therefore kept the presentation as general as possible, while citing recent experiments with thermostable and mesostable dehydrogenase and oxidase catalysts as specific examples exhibiting the phenomenon. We have also made a connection to the theory of diffusion. We hope that the point of view espoused in this paper will stimulate new ways of thinking about reaction rates, catalysis, and diffusion.
terms of Tolman’s interpretation (58), but any consistent explanation of a convex Arrhenius plot in terms of local equilibrium elementary rate coefficients must be consistent with Tolman’s result at some underlying level, even if this is not explicit or not even recognized. In summary, the present article has pointed out the fundamental microcanonical information that is provided by an observation of Arrhenius plot convexity, but which has apparently not previously been interpreted. The analysis has widespread applicability and makes a definite model-independent prediction about average energies of systems making transitions. The result is relevant to all fields of kinetics (in enzymes, in atmospheric
Truhlar and Kohen
PNAS 兩 January 30, 2001 兩 vol. 98 兩 no. 3 兩 851