Simply conceiving the Arrhenius law and absolute kinetic constants using the geometric distribution Denis Michel Universite de Rennes1-IRSET. Campus de Beaulieu Bat. 13. 35042 Rennes cedex. [email protected] Although first-order rate constants are mately from a strong initial impulsion spreading basic ingredients of physical chemistry, in the universe. Dynamic modeling of all molec- biochemistry and systems modeling, their ular systems, from enzymology to systems biol- innermost nature is derived from com- ogy, is grounded on the generalized use of kinetic plex physical chemistry mechanisms. The rates, usually written "k" in traditional biochem- present study suggests that equivalent istry. Molecular events cover a wide range of conclusions can be more straightly ob- phenomena including motion, covalent modifica- tained from simple statistics. The dif- tions or interactions. Although they are associ- ferent facets of kinetic constants are first ated to all these events, kinetic rates are subtle classified and clarified with respect to time enough. Their physical meaning is generally ig- and energy and the equivalences between nored by biochemists, who need only using the fa- traditional flux rate and modern proba- mous equation of Arrhenius [1] as a general rela- bilistic modeling are summarized. Then, tionship for the temperature dependence of rates a naive but rigorous approach is proposed of reaction. According to this equation, derived to concretely perceive how the Arrhenius from the thermodynamic study of van't Hoff [2], law naturally emerges from the geometric the logarithm of the rate constant is negatively distribution. It appears that (1) the dis- proportional to the inverse of temperature. Be- tribution in time of chemical events as well side this macroscopic thermodynamic approach, as (2) their mean frequency, are both dic- microscopic rate theories have then been devel- tated by randomness only and as such, are oped, for example based on the probability of accurately described by time-based and occurrence of a transition state theory [3] or the spatial exponential processes respectively. diffusion out of an energy well [4]. For complete reviews on the subject, see [5, 6, 7]. Very atom- Keywords: Arrhenius law; Rate constant; istic models will not be considered here but the Bose-Einstein distribution; Geometric law. simplified version of the transition state theory most widespread in biochemistry textbooks will be shown to suffer from several drawbacks. The arXiv:1307.4583v2 [physics.data-an] 19 Dec 2013 present study is aimed at recovering formally the Arrhenius law and establishing a simple defini- 1 Introduction tion of kinetic constants, expurgated from mech- anistic considerations, using minimalist hypothe- Over zero kelvin, matter particles move, collide ses and simple statistical tools. and react with a certain probability. All these events result from a transfer of energy provided by either radiation or previous motions and ulti- 1 2 Modeling the events nonlinear mechanisms responsible for the refined behaviors of living systems [9]. The basic condi- 2.1 The classical "mass action" ap- tion necessary to understand molecular events is proach to molecular events. that they occur randomly, that is to say without The very founder definition of chemical tran- memory. If a given event is memoryless and has sitions is based on rate fluxes, which are as- a mean waiting time hT i, then, the only possible sumed to be proportional to the amounts of law governing its probability of occurrence is the transformable reactants. In this view, which long exponential distribution. At every time point t, proved very efficient, the constants k are coeffi- this probability is simply cients of proportionality whose dimensions allow t − hT i to equalize the units in kinetic equations (time−1 P (X > t) = e (2a) for first order reactions, or M−1time−1 for second − t order reactions). For example, the elementary P (X ≤ t) = 1 − e hT i (2b) transition S !k P can be modeled as follows This result can be easily demonstrated. The d[S] − = k[S] (1a) absence of memory can be explicitly transcribed dt or into d[S] 8t ; t 2 = −kdt (1b) 1 2 R+ [S] P (X > t + t \ X > t ) = P (X > t ) (3) yielding after integration, 1 2 1 2 [S]t Indeed, P (X > t1 + t2jX > t1) = ln = −k(t − t0) (1c) [S]0 P (X > t1 + t2 \ X > t1)=P (X > t1) = P (X > or t1 + t2)=P (X > t1) −kt [S]t = [S]0e (1d) Reciprocally, Hence, Eq.(3) becomes −kt [P ]t = [S]0 − [S]t = [S]0(1 − e ) (1e) P (X > t1 + t2) = P (X > t1)P (X > t2) (4a) The probabilistic view described below can be anticipated from this classical mass action ap- This relationship imposes to handle t as an proach if identifying the proportions of molecules exponent, such that with probabilities, but the probabilistic view can be introduced even more straightly. e−k(t1+t2) = e−kt1 e−kt2 (4b) 2.2 Minimalist mathematical ap- proaches more suitable for single in which k is introduced to cancel the dimen- events in biochemistry sion of the exponent and turns to correspond to a frequency (time−1) and the average transition A living cell is not a dish as those used in chem- time is Z 1 1 istry labs, because certain subcellular compart- hT i = kt e−ktdt = (5) ments contain very few macromolecules of the t=0 k same kind. For example, the notion of concen- tration is obviously meaningless for a gene from 2.3 The ergodic correspondence the X chromosome, present in a single copy in The recent outset of single molecule-based per- each cell. To circumvent this problem, the single spectives, particularly in biochemistry, radically molecule mathematical approach is preferable, il- changed the traditional concepts of chemistry luminating for understanding biochemical mech- (Table 1). anisms [8], and now widely used for modeling the 2 Table 1: Correspondences between the traditional bulk approach of chemistry and single molecule- based perspective Mass action approach Single molecule approach Concentration State probability Rate flux Poisson parameter Total concentration of the leading macromolecule 1 (so that Vmax = kcat in enzymology) Space Time Of particular interest is the replacement of space (in which the concentration of reactants is defined), by time. This switch is fundamen- tally rooted in the principle of ergodicity, ac- cording to which the collection of states taken at a single time by an infinite number of parti- cles, corresponds to those taken by a single par- ticle during an infinite time window. Experi- Fig.1. The traditional energy barrier responsible for mental studies suggest that this assumption is reaction delays. acceptable, as illustrated by the single enzyme activity monitored in [10]. At equilibrium, the It is intriguing that when using this reasoning, classical Michaelis-Menten fraction of occupied the elementary rate k to be defined, is now re- z macromolecules can be defined as well as a frac- placed by 3 such rates (k1, k−1 and k2). The S tional occupation time of a single macromolecule intermediate is considered as short-lived. Ma- [11]. Although the principle of ergodicity is gen- nipulating rates as parameters of exponential erally presented as a non-demonstrated axiom, functions implies that the global transition from the identical results obtained from the mass ac- S to P takes the form k = k1k2=(k−1 + k2). New tion (section 2.1) and mathematical (section 2.2) hypotheses should be introduced to simplify this approaches allow to understand its basis. value 1st hypothesis, k2 >> k−1. Once gener- ated, Sz rapidly converts into P because of the 3 Modeling the frequencies of absence of energetic barrier between Sz and P. the events One haves k = k1 so that the definition of k is simply replaced by that of k1, no way solving the 3.1 Inconsistensies in the classical de- initial question: what is a k? scription of the transition state 2d hypothesis, k2 << k−1. This is the theory hypothesis classically retained in the activated complex theory. In this case, k = k2(k1=k−1), or z The version of the transition state theory most k = k2K . Now, the rate constant is defined by widespread in textbooks suffers from several an equilibrium (as initially proposed by Arrhe- drawbacks. Transitions from a starting point S nius himself, between normal and hypothetical to a final point P are supposed to follow a certain reaction-prone molecules). delay because they are restricted by an energy barrier. The transition is conditioned by prelim- z ∆Ez z [S ] − inary passage through an unstable reaction in- K = = e kB T (6) [S0] termediate named activated complex and usually written with a double-dagger label z. The result- from which the Arrhenius law can be recov- ing two-step reaction is represented in Fig.1. ered 3 black beads represent the separations between − Ea k = A e kB T (7) the particles. Adjacent white beads clustered be- tween two black beads correspond to the number In this expression, the mysterious constant of energy quanta (written ") in a single particle. "A", corresponding to a constant of integration If the beads added to the string are drawn ran- in the thermodynamic approach, is called "pre- domly from a bag containing a huge number of exponential coefficient" and is extremely impor- well-mixed black and white pearls (equivalent to tant since it gives to k the dimension of a rate a drawing with replacement), then the probabil- (time−1) [12]. In the present description, it cor- ity that any given bead on the string is white is responds to k . According to physical principles, 2 W=(W + B). Hence, the probability that an ar- this rate (written below kz) includes an universal bitrarily chosen stretch of white beads contains component k T=h time−1, corresponding to the B more than "z beads is mean vibrational or translational energy (kBT ) where T is the temperature and k is Boltz- z B W " mann's constant, divided by Planck's constant h.