# Biomod/2013/NanoUANL/Enzyme

(Difference between revisions)
 Revision as of 23:51, 11 October 2013 (view source)← Previous diff Current revision (22:56, 26 October 2013) (view source) (16 intermediate revisions not shown.) Line 1: Line 1: - ==What is an enzyme?== + {{Template:Nanouanl}} - In biological systems, chemical transformations are typically accelerated by enzymes, macromolecules capable of turning one or more compounds into others (substrates and products). The activity is determined greatly by their three-dimensional structure. Most enzymes are proteins, although several catalytic RNA molecules have been identified. They may also need to employ organic and inorganic cofactors for the reaction to occur. The process is based upon the diminishment of the activation energy needed for a reaction, greatly increasing its rate of reaction. The rate enhancement provided by these proteins can be as high as 10^19, while maintaining high substrate specificity. + - Because of this, reaction rates are millions of times faster than un-catalyzed reactions. Enzymes are not consumed by the reactions that they take part in, and they do not alter the equilibrium. Enzyme activity can be affected by a wide variety of factors. Inhibitors and activators intervene directly in the reaction rate, environmental factors like temperature, pressure, pH and substrate concentration also play a part in these kinetics. For temperature and pH, usually exist a range of values for which the enzyme works better (optimal conditions). The enzyme activity lowers dramatically as you get farther away from this range of values. As for concentration, other kind of relationship is observed. With increasing concentration, enzyme activity increases, until we reach the most optimal performance. Further increase of concentration won’t have an impact on the enzyme activity. + + + - Being able to determine these conditions allow us to manipulate the enzyme activity, thus achieving greater control over the reaction. + + + + +

+    ENZYME

+

+

+

+

+

+ Enzyme

+

+ In biological systems, chemical transformations are typically accelerated by enzymes, macromolecules capable of turning one or more compounds into others (substrates and products). The activity is determined greatly by their three-dimensional structure. Most enzymes are proteins, although several catalytic RNA molecules have been identified. They may also need to employ organic and inorganic cofactors for the reaction to occur. The process is based upon the diminishment of the activation energy needed for a reaction, greatly increasing its rate of reaction. The rate enhancement provided by these proteins can be as high as 10^19, while maintaining high substrate specificity.

+

+ Because of this, reaction rates are millions of times faster than un-catalyzed reactions. Enzymes are not consumed by the reactions that they take part in, and they do not alter the equilibrium. Enzyme activity can be affected by a wide variety of factors. Inhibitors and activators intervene directly in the reaction rate, environmental factors like temperature, pressure, pH and substrate concentration also play a part in these kinetics. For temperature and pH, usually exist a range of values for which the enzyme works better (optimal conditions). The enzyme activity lowers dramatically as you get farther away from this range of values. As for concentration, other kind of relationship is observed. With increasing concentration, enzyme activity increases, until we reach the most optimal performance. Further increase of concentration generally won’t have an impact on the enzyme activity.

+

+ Being able to determine these conditions allow us to manipulate the enzyme activity, thus achieving greater control over the reaction.

+

+ Michaelis-Menten kinetics

+

+ Michaelis-Menten kinetics is one of the oldest models for describing the catalytic activity of enzymes. The reaction cycle is divided into two basic steps: the reversible binding between the enzyme and substrate to form an intermediate complex, and the irreversible catalytic step to generate the product and release the enzyme; in which the first step is affected by the constants k1 and k-1, whereas the irreversible step only takes into account  k2.

- ===Horseradish peroxidase=== + + E + S \leftrightarrow ES \rightarrow E^0 + P + +

+ The rate of consumption can be expressed by the formation of the ES complex in the following equation:

- Enzyme information was gathered from research papers: + + \frac{d[ES]}{dt}=k_1[E][S]-k_{-1}[ES]-k_2[ES] + - * Isoelectric Point: Isozymes range from 3.0 – 9.0 pH +

- * pH Dependence: range of 5.0 to 9.0, optimum (6.0 to 6.5) + Using a steady-approximation and rearranging eq. 2 we obtain:

- * k1 = 1 x 107 + - * k-1 = between 20 and 250 + - * k2 = 200 at pH 6 + - ===Enzyme kinetics=== + - Michaelis-Menten kinetics is one of the oldest models for describing the catalytic activity of enzymes. The reaction cycle is divided into two basic steps: the reversible binding between the enzyme and substrate to form an intermediate complex, and the irreversible catalytic step to generate the product and release the enzyme; in which the first step is affected by the constants ''k1;'' and ''k-1'', whereas the irreversible step only takes into account ''k2''. + [ES]= \frac{[E][S]}{K_M+[S]} + - '''1.1''' +

+ where ''KM'' is the Michaelis constant defined as

- $+ - E + S \leftrightarrow ES \rightarrow E^0 + P + K_M= \frac{k_{-1} + k_2}{k_1} -$ + - The rate of consumption can be expressed by the formation of the ES complex in the following equation: +

+ As it was mentioned in the introduction, single-enzyme studies have proven that the "traditional" enzyme kinetics do not apply, and a new approach is needed. Enzyme concentration is meaningless in a single-molecule level, so it is more appropriate to consider the probability ''PE(t)'' for the enzyme to find a catalytically active enzyme in a time ''t'' in the process. This is because the reaction is a stochastic event.

+

+ Therefore, the rate equations of each species are:

- '''1.2''' + - $+ \frac{d[E]}{dt}=-k_1[E][S]+k_{-1}[ES] - \frac{d[ES]}{dt}=k_1[E][S]-k_{-1}[ES]-k_2[ES] + -$ + - Using a steady-approximation and rearranging 1.2 we obtain: + + \frac{d[ES]}{dt}=k_1[E][S]-(k_{-1}+k_2)[ES] + - '''1.3''' + - $+ \frac{d[E^0]}{dt}=\frac{d[P]}{dt}=k_2[ES] - [ES]= \frac{[E][S]}{K_M+[S]} + \end{equation} - + - where ''KM'' is the Michaelis constant defined as + + where ''t'' is the elapsed time, the initial conditions are [ES]=0 and [E<sup>0</sup>]=0 at ''t''=0. To derive the rate equations that describe the corresponding single-molecule Michaelis-Menten kinetics, the concentrations in equations 5-7 are replaced by the probabilities P of finding the single enzyme molecule in the states E, ES, and E<sup>0</sup> , leading to the equations: - '''1.4''' + + \frac{dP_E(t)}{dt}=-k_1^0P_E(t)+k_{-1}P_{ES}(t) + - [itex] + \begin{equation} - \operatorname{K_M}= + \frac{dP_{ES}(t)}{dt}=k_1^0P_E(t)-(k_{-1}+k_2)P_{ES}(t) - \tfrac{k_{-1} +k_2}{k_1} + -$ + - As it was mentioned in the introduction, single-enzyme studies have proven that the "traditional" enzyme kinetics do not apply, and a new approach is needed. Enzyme concentration is meaningless in a single-molecule level, so it is more appropriate to consider the probability ''PE(t)'' for the enzyme to find a catalytically active enzyme in a time ''t'' in the process. This is because the reaction is a stochastic event. + + \frac{dP_E^0(t)}{dt}=k_2P_{ES}(t) + - Therefore, the rate equations of each species are: +

+ These equations must satisfy the conditions ''PE''(0)=1, ''PES''(0)=0 and ''PE0''=0 at ''t''=0 (start of the reaction). Also, ''PE''(t) + ''PES''(t) + ''PE0''(t)=1. The rate constant ''k''10 can be taken as ''k''10=''k''1[S], assuming [S] is time-independent. Given that a single enzyme is unlikely to deplete all the substrate presence, [S] can be considered constant, virtually being unaffected.

- '''1.5''' +

- \frac{d[E]}{dt}=-k_1[E][S]+k_{-1}[ES] + Equations 8-10 become a system of linear first-order differential equations that can be solved exactly for ''PE''(t), ''PES''(t) and ''PE0''(t).

- '''1.6''' +

- \frac{d[ES]}{dt}=k_1[E][S]-(k_{-1}+k_2)[ES] + Knowing ''PE0''(t), the probability that a turnover occurs between ''t'' and ''t +'' Δ''t'' is ''f(t)Δt, the same as Δ''PE0''(t). Taking this into account, in the limit of infinitesimal Δ''t'':

- '''1.7''' + - $\frac{d[E^0]}{dt}=\tfrac{d[P]}{dt}=k_2[ES]$ + f(t)=\frac{dP_E^0(t)}{dt}=k_2P_{ES}(t) + \end{equation} - where ''t'' is the elapsed time, the initial conditions are [ES]=0 and [E0]=0 at ''t''=0. To derive the rate equations that describe the corresponding single-molecule Michaelis-Menten kinetics, the concentrations in equations 5-7 are replaced by the probabilities P of finding the single enzyme molecule in the states E, ES, and E0 , leading to the equations: +

+ Solving equations 8-10, and using equation 11:

- '''1.8''' + \begin{equation} - $\frac{dP_E(t)}{dt}=-k_1^0P_E(t)+k_{-1}P_{ES}(t)$ + f(t)=\frac{k_1k_2[S]}{2A}[exp(A+B)t-exp(B-A)t] + - '''1.9''' +

- \frac{dP_{ES}(t)]}{dt}=k_1^0P_E(t)-(k_{-1}+k_2)P_{ES}(t) + in which:

- '''1.10''' + \begin{equation} - $\frac{dP_E^0(t)}{dt}=k_2P_{ES}(t) + A=\sqrt{(k_1[S]+k_{-1}+k_2)^2/4-k_1k_2[S]} + - These equations must satisfy the conditions ''PE''(0)=1, ''PES''(0)=0 and ''PE0''=0 at ''t''=0 (start of the reaction). Also, ''PE''(t) + ''PES''(t) + ''PE0''(t)=1. The rate constant ''k''10 can be taken as ''k''10=''k''1[S], assuming [S] is time-independent. Given that a single enzyme is unlikely to deplete all the substrate presence, [S] can be considered constant, virtually being unaffected. + + B=- \frac{k_1[S]+k_{-1}+k_2}{2} + - Equations 8-10 become a system of linear first-order differential equations that can be solved exactly for ''PE''(t), ''PES''(t) and ''PE0''(t). + + when the substrate concentration dependence [S] has been shown throught the relation ''k10 = k1''[S]. - Knowing ''PE0''(t), the probability that a turnover occurs between ''t'' and ''t +'' Δ''t'' is ''f(t)Δt, the same as Δ''PE0''(t). Taking this into account, in the limit of infinitesimal Δ''t'': + + With these equations, we have used different values for ''k-1'' and [S], using values found in literature. - '''1.11''' + - f(t)=\frac{dP_E^0(t)}{dt}=k_2P_{ES}(t) + + At the first moment of ''f(t)'', which would be the mean waiting time for the reaction, ''<t>'', and its reciprocal can be taken as the average reaction rate. Starting eq. 12: - Solving equations 8-10, and using equation 11: + + \frac{1}{\langle t \rangle)}=-\frac{(A^2-B^2)^2}{2Bk_1k_2[S]} + + + + \frac{1}{\langle t \rangle}=\frac{k_2[S]}{[S]+K_M} + + + + Randomness parameter + The probability density ''f(t)'' completely characterizes single-enzyme kinetics, with the ''n''<sup>th</sup> moment being given by: + + \begin{equation} + f(t)= \int_0^\infty dt f(t)t^n + + + + Although the first moment of ''f(t)'' can be described eq. 11, higher moments of ''f(t)'' are usually calculated along with a "randomness parameter". Implying no dynamic disorder, ''r'' is given by: + + + r=\frac{(k_1[S]+k_2+k_{-1})^2-2k_1k_2[S]}{(k1[S]+k_2+k_{-1})^2} + + + + As substrate concentration increases, ''r'' decreases, indicating the formation of an intermediate enzyme-substrate complex ES. At even higher concentrations, the catalytic step limits the reaction, as is the same when concentration is low and the substrate binding limits the rate. + + + + Enzyme Reaction in a Container + + Some reactions occur in spaces which have such a small volume that the normal approach to enzyme kinetics does not work. Even when taking into account probability instead of concentration, reaction cannot be modeled by aninfinite volume with a constant density. Examples of these systems include nanoparticles, zeolites or silica gels. In systems containing a small number of interacting species a relatively large fluctuationin the number of reactants is inherent. Therefore, the partitioning of the reagentspecies in small volume and the small number of the species present in eachindividual volume make it difficult to apply classical chemical kinetics to describe individual reactions in the container. + + For a spherical system, in which molecules A and B react, the survival probability p(r,t) at the time t when A,B are at a d distance: + + Enzyme Reaction in a Container + + Some reactions occur in spaces which have such a small volume that the normal approach to enzyme kinetics does not work. Even when taking into account probability instead of concentration, reaction cannot be modeled by aninfinite volume with a constant density. Examples of these systems include nanoparticles, zeolites or silica gels. In systems containing a small number of interacting species a relatively large fluctuationin the number of reactants is inherent. Therefore, the partitioning of the reagentspecies in small volume and the small number of the species present in eachindividual volume make it difficult to apply classical chemical kinetics to describe individual reactions in the container. + + For a spherical system, in which molecules A and B react, the survival probability p(r,t) at the time t when A,B are at a d distance: - '''1.12''' - [itex] - f(t)=\frac{k_1k_2[S]}{2A}[exp(A+B)t-exp(B-A)t] -$ - in which: + + \frac{δp(r,t)}{δt}=D\bigtriangledown^2p(r,t) + - '''1.13''' + - $+ p(r,0)=1 - A=\sqrt((k_1[S]+k_{-1}+k_2)^2/4-k_1k_2[S]) + \end{equation} -$ + - '''1.14''' + - $+ p(d,t)=0 - B=\frac{-(k_1[S]+k_{-1}+k_2}{2} + \end{equation} -$ + - when the substrate concentration dependence [S] has been shown throught the relation ''k10 = k1''[S]. + + \bigg[ \frac{δp(r,t)}{δr} \bigg]_r=R =0 + - With these equations, we have used different values for ''k-1'' and [S], using values found in literature. +

+ These equations describe the random motion of A with respect to B, the system’s initial condition, distance needed to interact, and finally a boundary condition for the sphere.

+

+ In a general case where a sphere contains one B molecule and a number N of Amolecules that move mutually independently, the survival probability P(t) can bewritten as a product of the survival probabilities of all (A…B) pairs:

+ + P(t) = \prod_{i=1}^N p(r,t) + +

+ In this equation, riis the initial position of the ithA molecule. Averaging this last equation over therandom initial distribution of A molecules inside the sphere results in:

- [Descripcion y link del codigo de simulink] + + P(r,t)= \bigg[\sum_{i=1}^\infty H_i•exp(-α_i^2 \frac{Dt}{R^2}\bigg]^N + - [GRAFICAS Y ESO] + + H_i=\frac{6\bigg(\frac{d}{R}\bigg)^2}{\alpha_i^2\bigg(1-\frac{d}{R}^2\bigg)\bigg[\frac{\alpha_i^2}{1+\alpha_i^2}-\frac{d}{R}^2\bigg]} + +

+ Thus, P(r,t) decays exponentially, indicating that the contained reaction obeys first-order kinetics.

- At the first moment of ''f(t)'', = \int_0^\infty dt tf(t), which would be the mean waiting time for the reaction, '''', and its reciprocal can be taken as the average reaction rate. Starting eqtn 12: +

+ RESULTS

+

+

+

+ Single-enzyme Michalis-Menten

+

+ After the solution of the differential equations system (eqs. 8-10), the analytical function of waiting time distribution f(t) was obtained. Then we introduced our parameters to understand the behavior of our specific system. The following constants obtained from literature:

+ + $$k_1=10^7 M^{-1} s^{-1}$$ +

+

+ $$k_2=201 s^{-1}$$ +

+

+ $$k_{-1}=50-250 s^{-1}$$ +

+

+ $$[S]=0.0001 - 0.1 mM$$ - $\frac{1}{}=-\frac{(A^2-B^2)^2}{2Bk_1k_2[S]}$ - $\frac{1}{}=\frac{k_2[S]}{[S]+K_M}$ - ===Randomness parameter=== +

- The probability density ''f(t)'' completely characterizes single-enzyme kinetics, with the ''n''th moment being given by: +

+

+ Using Matlab for calculating and drawing the plots for this data, we obtained the following figures, varying the values of [S]:

+

+

+

+ Fig.1 Waiting time distribution, as a function of time and [S]

+

+ At concentrations lower than 0.01 M, we can observe the lowest values of waiting time distributions. With increasing concentrations, the maximum contributions of waiting time distribution were observed around the first 0.02 seconds. After this time, their values drop dramatically, as we expected, because it is less probable to have a reaction after a longer elapsed time (because it must have already happened).

+

+ The same was done varying the values of k-1:

+

+

+

+ Fig.2 Waiting time distribution, as a function of time and k-1

+

+ We can observe a change in the height of the curve´s peak aroung 0 and 0.02 seconds, as the values of k-1 increase. However, the same general behavior is observed for the waiting time distribution.

+

+ Then, with this function, we obtained the probability of survival (no reaction between the enzyme and the substrate) by the following equation:

+ + $$P(t)= 1- \int f(t) dt$$ - '''1.15''' +

- f(t)= \int_0^\infty dt f(t)t^n +

+

+ Fig.3 Probability of survival, as a function of time with varying [S] for the single-enzyme Michaelis-Menten model

+

+ As expected the probability drops dramatically as time increases, and drops faster for higher [S].

+

+

+

+ Confined space model

+

+ To prove the feasibility of the confined space model, we tested the special case of a container with an infinite radius, which has to behave as the previous model. For this, the probability of survival simplifies as follows:

+

+

+ + P(t)=exp(-4\pi cdDt)\end{equation} +

+ where c is the concentration of the reactive species. During the reaction model, it was found that d = 1.6 x 1012 A2s-1. Again, using Matlab to compute these equations, we obtained the following figure:

+

+

+

+ Fig.4 Probability of survival, as a function of time with varying [S] for the confined space model

+

+ The general behavior of this model is very similar to that of the single-enzyme Michaelis-Menten model, although the timescale of the former is in a different order of magnitude. We speculate that this is because the confined space model does is only taking into account geometric considerations (the reaction occurs when the molecules are close enough to each other). This model does not consider many other variables of the phenomenon, such as the location of active sites and substrate specificity, which are taken into account experimentaly by the rate of reaction coefficients in the single-enzyme Michaelis-Menten model. As such, it is apparent that the confined space model is a good model for this phenomenon, but needs to take into account more specific environmental conditions to fit with experimental data before applying it to confined spaces.

- Although the first moment of ''f(t)'' can be described eqtn 11, higher moments of ''f(t)'' are usually calculated along with a "randomness parameter".  Implying no dynamic disorder, ''r'' is given by: - '''1.16''' - $r=\frac{(k_1[S]+k_2+k_{-1})^2-2k_1k_2[S]}{(k1[S]+k_2+k_{-1})^2}$ - As substrate concentration increases, ''r'' decreases, indicating the formation of an intermediate enzyme-substrate complex ES. At even higher concentrations, the catalytic step limits the reaction, as is the same when concentration is low and the substrate binding limits the rate. - ===Conclusions=== + - We calculated the probability density of the waiting times of the enzymatic reaction, taking into account that the process is a mainly stochastic event at a single-molecule level. The reciprocal of the first moment of ''f(t)'' appears to have the same hyperbolic dependence on [S] as is described in conventional Michaelis-Menten kinetics (when assuming steady-state). +

## Current revision

ENZYME

Enzyme

In biological systems, chemical transformations are typically accelerated by enzymes, macromolecules capable of turning one or more compounds into others (substrates and products). The activity is determined greatly by their three-dimensional structure. Most enzymes are proteins, although several catalytic RNA molecules have been identified. They may also need to employ organic and inorganic cofactors for the reaction to occur. The process is based upon the diminishment of the activation energy needed for a reaction, greatly increasing its rate of reaction. The rate enhancement provided by these proteins can be as high as 10^19, while maintaining high substrate specificity.

Because of this, reaction rates are millions of times faster than un-catalyzed reactions. Enzymes are not consumed by the reactions that they take part in, and they do not alter the equilibrium. Enzyme activity can be affected by a wide variety of factors. Inhibitors and activators intervene directly in the reaction rate, environmental factors like temperature, pressure, pH and substrate concentration also play a part in these kinetics. For temperature and pH, usually exist a range of values for which the enzyme works better (optimal conditions). The enzyme activity lowers dramatically as you get farther away from this range of values. As for concentration, other kind of relationship is observed. With increasing concentration, enzyme activity increases, until we reach the most optimal performance. Further increase of concentration generally won’t have an impact on the enzyme activity.

Being able to determine these conditions allow us to manipulate the enzyme activity, thus achieving greater control over the reaction.

Michaelis-Menten kinetics

Michaelis-Menten kinetics is one of the oldest models for describing the catalytic activity of enzymes. The reaction cycle is divided into two basic steps: the reversible binding between the enzyme and substrate to form an intermediate complex, and the irreversible catalytic step to generate the product and release the enzyme; in which the first step is affected by the constants k1 and k-1, whereas the irreversible step only takes into account  k2.

$$E + S \leftrightarrow ES \rightarrow E^0 + P$$

The rate of consumption can be expressed by the formation of the ES complex in the following equation:

$$\frac{d[ES]}{dt}=k_1[E][S]-k_{-1}[ES]-k_2[ES]$$

Using a steady-approximation and rearranging eq. 2 we obtain:

$$[ES]= \frac{[E][S]}{K_M+[S]}$$

where ''KM'' is the Michaelis constant defined as

$$K_M= \frac{k_{-1} + k_2}{k_1}$$

As it was mentioned in the introduction, single-enzyme studies have proven that the "traditional" enzyme kinetics do not apply, and a new approach is needed. Enzyme concentration is meaningless in a single-molecule level, so it is more appropriate to consider the probability ''PE(t)'' for the enzyme to find a catalytically active enzyme in a time ''t'' in the process. This is because the reaction is a stochastic event.

Therefore, the rate equations of each species are:

$$\frac{d[E]}{dt}=-k_1[E][S]+k_{-1}[ES]$$ $$\frac{d[ES]}{dt}=k_1[E][S]-(k_{-1}+k_2)[ES]$$ $$\frac{d[E^0]}{dt}=\frac{d[P]}{dt}=k_2[ES]$$

where ''t'' is the elapsed time, the initial conditions are [ES]=0 and [E<sup>0</sup>]=0 at ''t''=0. To derive the rate equations that describe the corresponding single-molecule Michaelis-Menten kinetics, the concentrations in equations 5-7 are replaced by the probabilities P of finding the single enzyme molecule in the states E, ES, and E<sup>0</sup> , leading to the equations:

$$\frac{dP_E(t)}{dt}=-k_1^0P_E(t)+k_{-1}P_{ES}(t)$$ $$\frac{dP_{ES}(t)}{dt}=k_1^0P_E(t)-(k_{-1}+k_2)P_{ES}(t)$$ $$\frac{dP_E^0(t)}{dt}=k_2P_{ES}(t)$$

These equations must satisfy the conditions ''PE''(0)=1, ''PES''(0)=0 and ''PE0''=0 at ''t''=0 (start of the reaction). Also, ''PE''(t) + ''PES''(t) + ''PE0''(t)=1. The rate constant ''k''10 can be taken as ''k''10=''k''1[S], assuming [S] is time-independent. Given that a single enzyme is unlikely to deplete all the substrate presence, [S] can be considered constant, virtually being unaffected.

Equations 8-10 become a system of linear first-order differential equations that can be solved exactly for ''PE''(t), ''PES''(t) and ''PE0''(t).

Knowing ''PE0''(t), the probability that a turnover occurs between ''t'' and ''t +'' Δ''t'' is ''f(t)Δt, the same as Δ''PE0''(t). Taking this into account, in the limit of infinitesimal Δ''t'':

$$f(t)=\frac{dP_E^0(t)}{dt}=k_2P_{ES}(t)$$

Solving equations 8-10, and using equation 11:

$$f(t)=\frac{k_1k_2[S]}{2A}[exp(A+B)t-exp(B-A)t]$$

in which:

$$A=\sqrt{(k_1[S]+k_{-1}+k_2)^2/4-k_1k_2[S]}$$ $$B=- \frac{k_1[S]+k_{-1}+k_2}{2}$$

when the substrate concentration dependence [S] has been shown throught the relation ''k10 = k1''[S].

With these equations, we have used different values for ''k-1'' and [S], using values found in literature.

At the first moment of ''f(t)'', which would be the mean waiting time for the reaction, ''<t>'', and its reciprocal can be taken as the average reaction rate.  Starting eq. 12:

$$\frac{1}{\langle t \rangle)}=-\frac{(A^2-B^2)^2}{2Bk_1k_2[S]}$$ $$\frac{1}{\langle t \rangle}=\frac{k_2[S]}{[S]+K_M}$$

Randomness parameter
The probability density ''f(t)'' completely characterizes single-enzyme kinetics, with the ''n''<sup>th</sup> moment being given by:

$$f(t)= \int_0^\infty dt f(t)t^n$$

Although the first moment of ''f(t)'' can be described eq. 11, higher moments of ''f(t)'' are usually calculated along with a "randomness parameter". Implying no dynamic disorder, ''r'' is given by:

$$r=\frac{(k_1[S]+k_2+k_{-1})^2-2k_1k_2[S]}{(k1[S]+k_2+k_{-1})^2}$$

As substrate concentration increases, ''r'' decreases, indicating the formation of an intermediate enzyme-substrate complex ES. At even higher concentrations, the catalytic step limits the reaction, as is the same when concentration is low and the substrate binding limits the rate.

Enzyme Reaction in a Container

Some reactions occur in spaces which have such a small volume that the normal approach to enzyme kinetics does not work. Even when taking into account probability instead of concentration, reaction cannot be modeled by aninfinite volume with a constant density. Examples of these systems include nanoparticles, zeolites or silica gels. In systems containing a small number of interacting species a relatively large fluctuationin the number of reactants is inherent. Therefore, the partitioning of the reagentspecies in small volume and the small number of the species present in eachindividual volume make it difficult to apply classical chemical kinetics to describe individual reactions in the container.

For a spherical system, in which molecules A and B react, the survival probability p(r,t) at the time t when A,B are at a d distance:

Enzyme Reaction in a Container

Some reactions occur in spaces which have such a small volume that the normal approach to enzyme kinetics does not work. Even when taking into account probability instead of concentration, reaction cannot be modeled by aninfinite volume with a constant density. Examples of these systems include nanoparticles, zeolites or silica gels. In systems containing a small number of interacting species a relatively large fluctuationin the number of reactants is inherent. Therefore, the partitioning of the reagentspecies in small volume and the small number of the species present in eachindividual volume make it difficult to apply classical chemical kinetics to describe individual reactions in the container.

For a spherical system, in which molecules A and B react, the survival probability p(r,t) at the time t when A,B are at a d distance:

$$\frac{δp(r,t)}{δt}=D\bigtriangledown^2p(r,t)$$ $$p(r,0)=1$$ $$p(d,t)=0$$ $$\bigg[ \frac{δp(r,t)}{δr} \bigg]_r=R =0$$

These equations describe the random motion of A with respect to B, the system’s initial condition, distance needed to interact, and finally a boundary condition for the sphere.

In a general case where a sphere contains one B molecule and a number N of Amolecules that move mutually independently, the survival probability P(t) can bewritten as a product of the survival probabilities of all (A…B) pairs:

$$P(t) = \prod_{i=1}^N p(r,t)$$

In this equation, riis the initial position of the ithA molecule. Averaging this last equation over therandom initial distribution of A molecules inside the sphere results in:

$$P(r,t)= \bigg[\sum_{i=1}^\infty H_i•exp(-α_i^2 \frac{Dt}{R^2}\bigg]^N$$ $$H_i=\frac{6\bigg(\frac{d}{R}\bigg)^2}{\alpha_i^2\bigg(1-\frac{d}{R}^2\bigg)\bigg[\frac{\alpha_i^2}{1+\alpha_i^2}-\frac{d}{R}^2\bigg]}$$

Thus, P(r,t) decays exponentially, indicating that the contained reaction obeys first-order kinetics.

RESULTS

Single-enzyme Michalis-Menten

After the solution of the differential equations system (eqs. 8-10), the analytical function of waiting time distribution f(t) was obtained. Then we introduced our parameters to understand the behavior of our specific system. The following constants obtained from literature:

$$k_1=10^7 M^{-1} s^{-1}$$

$$k_2=201 s^{-1}$$

$$k_{-1}=50-250 s^{-1}$$

$$[S]=0.0001 - 0.1 mM$$

Using Matlab for calculating and drawing the plots for this data, we obtained the following figures, varying the values of [S]:

Fig.1 Waiting time distribution, as a function of time and [S]

At concentrations lower than 0.01 M, we can observe the lowest values of waiting time distributions. With increasing concentrations, the maximum contributions of waiting time distribution were observed around the first 0.02 seconds. After this time, their values drop dramatically, as we expected, because it is less probable to have a reaction after a longer elapsed time (because it must have already happened).

The same was done varying the values of k-1:

Fig.2 Waiting time distribution, as a function of time and k-1

We can observe a change in the height of the curve´s peak aroung 0 and 0.02 seconds, as the values of k-1 increase. However, the same general behavior is observed for the waiting time distribution.

Then, with this function, we obtained the probability of survival (no reaction between the enzyme and the substrate) by the following equation:

$$P(t)= 1- \int f(t) dt$$

Fig.3 Probability of survival, as a function of time with varying [S] for the single-enzyme Michaelis-Menten model

As expected the probability drops dramatically as time increases, and drops faster for higher [S].

Confined space model

To prove the feasibility of the confined space model, we tested the special case of a container with an infinite radius, which has to behave as the previous model. For this, the probability of survival simplifies as follows:

$$P(t)=exp(-4\pi cdDt)$$

where c is the concentration of the reactive species. During the reaction model, it was found that d = 1.6 x 1012 A2s-1. Again, using Matlab to compute these equations, we obtained the following figure:

Fig.4 Probability of survival, as a function of time with varying [S] for the confined space model

The general behavior of this model is very similar to that of the single-enzyme Michaelis-Menten model, although the timescale of the former is in a different order of magnitude. We speculate that this is because the confined space model does is only taking into account geometric considerations (the reaction occurs when the molecules are close enough to each other). This model does not consider many other variables of the phenomenon, such as the location of active sites and substrate specificity, which are taken into account experimentaly by the rate of reaction coefficients in the single-enzyme Michaelis-Menten model. As such, it is apparent that the confined space model is a good model for this phenomenon, but needs to take into account more specific environmental conditions to fit with experimental data before applying it to confined spaces.