Team:UNSW Australia/Model/EKD

Enzyme Kinetics and Diffusion


Why Make a Mathematical Model?

Previous experiments involving enzyme scaffolding have shown various yields, ranging from fractional improvement to 77-fold increases in product titre 2,7. The exact improvements on reaction efficiency made via scaffolding enzymes vary across experiments. Notably, different scaffolds, enzymes, and substrates may result in significantly different yields 2. We are employing the enzymes used in Indole Acetic Acid synthesis in our system, utilizing their kinetic parameters to dictate the model. Additionally we model the taxol synthesis pathway as it is a potential commercial application of our scaffold. We consider the diffusion of the intermediary substrate when the enzymes are placed various distances apart. We believe the model was able to assist in our understanding of how scaffolding enzymes can alter reaction efficiency.

Aim of Modeling

We aimed to model the effect of distance between enzymes on reaction efficiency, hypothesizing that co-localising enzymes would increase reaction efficiency by decreasing intermediary loss to diffusion. Under this model, we assumed that intermediary loss to diffusion is significant enough to limit rates of product formation. We further aimed to incorporate the specific kinetics of enzymes in our system, and create an adaptable model that could generate data for other enzyme reaction pathways.

Boundary Conditions

To define the scope of the model, some initial boundary conditions were required. The boundary conditions are as follows, with justification following:

  1. The first enzyme produces the second substrate at a constant rate.
  2. The enzymes are some distance 0 < R < ∞ apart (they are not touching, nor are they infinitely far apart).
  3. The rate of production of the final product is dependent on the concentration of the second substrate that diffuses to enzyme B.
  4. The enzymes are not points, but rather small spheres given the size of the enzymes relative to the distance between them.
  5. The reactant can escape to infinity.

Condition (1) creates a linear diffusion gradient for the second substrate’s diffusion to enzyme B. This follows Fick’s Law of diffusion 4. In order to meet this condition the concentration of the first substrate must be held constant. Since the first substrate is being consumed by enzyme A, it must be kept in flow to the system.

Condition (2) describes the physical nature of the system; the enzymes attached to prefoldin are a finite distance apart.

For Condition (3), the second enzyme functions as a sink, removing S2 particles from the system and creating the final product IAA. Similarly the rate of reaction for enzyme B is dependent on the concentration of substrate 2 present; as with any enzyme reaction, thus it is a first order sink 6.

Condition (4) is important to define the diffusion gradient and was verified by the molecular dynamics simulations.

The final condition determines how the model treats the concentration of substrate at very large distances from the source enzyme. Allowing the substrate to escape to infinity implies the concentration tends to zero as distance tends to infinity. This avoids the complexities of determining how long it might take for the reaction to complete by considering the probability of the prefoldin-enzyme complex and escaped substrate colliding due to random walk 3.

Enzyme Kinetics

The IAA reaction is a 2-step sequential reaction:

Figure 1: The enzymatic patway for the production of IAA.

Each enzyme interacts with a single substrate to produce a single product, and each enzyme has its own set of kinetic parameters. We employed the indole acetic acid pathway as a model enzyme reaction, as the enzymes are readily available with well-defined kinetic parameters.

The Michaelis-Menten (MM) equation (equation 1) describes the kinetics of an enzyme-substrate interaction under a given set of environmental conditions. MM kinetics were used to describe the effective rates of product formation and consumption at each enzyme in the sequential reaction.

$$\frac{d[P]}{dt} = \frac{Vmax[S]}{Km + [S]} \tag{1} $$

The Km or Michaelis constant is an affinity constant for the particular enzyme-substrate interaction, whilst the Vmax is the maximum speed of the enzyme when substrate concentration is saturated. The equation is pseudo first order in this sense, as substrate concentration is high, enzyme speed (V) approaches Vmax, but whilst it is low, the MM equation holds true.

To simplify the model, we assume that the first substrate (tryptophan) is maintained in the system at a constant, saturating concentration, thus the first enzyme produces the second substrate (IAM) at a constant rate. This leads to a simple, linear diffusion of IAM from enzyme A. Thus, the production rate at enzyme A of substrate 2 is simply the Vmax of enzyme A (equation 2).

$$\frac{d[S2]}{dt} = Vmax \tag{2} $$

During consultations with industry about our model, the idea of keeping the first substrate in constant supply to the system arose. Dr Warren King suggested that it would be preferable in a commercial setting, as continuous systems generally make more efficient use of resources. Dr King also confirmed it would simplify the remainder of our analysis in developing the model.

Fick’s law is used to describe the diffusion of IAM from enzyme A into the system in two dimensions (equation 3).

$$J = -D \nabla \varphi \tag{3} $$

Where J is the diffusion flux in amount per unit area per unit time, D is the diffusivity in area per unit time, and ϕ denotes the concentration in amount per unit volume. The driving force of diffusion using Fick’s law is a concentration gradient. The concentration gradient in this model is determined by the difference in production and consumption of substrate 2 between the first and second enzyme. Thus the rate of substrate consumption at the second enzyme will determine the diffusive flux toward second enzyme.

We define the second enzyme in the system as the sink for diffusion, with consumption rate a 1st order reaction defined by the Michaelis-Menten equation. At the edge of the second enzyme (represented by E4 on the diagram), the consuming boundary condition is defined by equation 3.

$$\frac{d[P]}{dt} = \frac{Vmax[S2]}{Km + [S2]} \tag{4}$$

Vmax and Km are now specific to the second enzyme; indole acetamide hydrolase. As the substrate needs to diffuse from enzyme A, the concentration of substrate reaching enzyme B is expected to be low, therefore the reaction rate at the second enzyme is dependent on the amount of substrate reaching enzyme B through diffusion. Thus the 1st order form of the Michaelis-Menten equation applies.