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. 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:
- The first enzyme produces the second substrate at a constant rate.
- The enzymes are some distance 0 < R < ∞ apart (they are not touching, nor are they infinitely far apart).
- The rate of production of the final product is dependent on the concentration of the second substrate that diffuses to enzyme B.
- The enzymes are not points, but rather small spheres given the size of the enzymes relative to the distance between them.
- The reactant can escape to infinity.
Condition (1) creates a linear diffusion gradient for the second substrate’s diffusion to enzyme 2. 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 1, 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 2 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.
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).
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 1 into the system in two dimensions (equation 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.
Vmax and Km are now specific to the second enzyme; indole acetimide 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.
Our System
Figure 2: Diffusion of substrate from Enzyme A to Enzyme B.
Assumptions of This Model
At each enzyme, employing the Michaelis-Menten equation assumes the law of mass action, meaning all substrate that reaches the boundary at each enzyme is completely converted to product 5. In reality, product levels are expected to be slightly lower than predicted due to imperfect substrate-enzyme interactions.
Employing Fick’s law of diffusion assumes that a concentration gradient is the only force driving substrate motion in the system. In reality, charge interactions may create a substrate channelling effect, or substrate-enzyme repulsion, leading to increased or decreased true product yields when compared to the predicted data 8. Fick’s law also does not account for the shape of each molecule, and how angular momentum of each of the substrate may affect motion of the substrate between enzymes 9.
The distance between enzymes is assumed to be relatively constant as the reaction takes place. In reality, the enzymes may be moving rapidly and non-constantly, creating variation in the amount of product produced over a given time 9. By employing a constant distance model we are effectively taking an average of the distance between enzymes.
Simplifying the geometry to 2 dimensions, and a half sphere further assumes that diffusion would not be biased in any directional plane in the real system. This assumes an equal and opposite flux in the plane perpendicular to the plane of distance between enzymes.
Figure 2: Geometry of the model developed in MATLAB looking at diffusion.
Approximating the enzymes as spheres is the last assumption employed to simplify the system, the enzymes are approximated as spheres as opposed to points due to the size of the enzymes relative to the distance between them (either reference MD or reference trigonometric analysis). Additionally the exact position of the active site of the enzyme is not considered in this case. This may have an affect on the accuracy of the model in cases where the active site of the enzyme is pointed away from the next enzyme in the reaction pathway 2. A more realistic model would involve assigning the site of reaction (production or consumption) to be a point on each of the half-spheres. The point could then move over the surface of the half sphere, mimicking the changing orientation of the active site in a real enzyme reaction. This was not possible due to an inability to assign point-node boundary conditions on element meshes in MATLAB, however other more specialised software for solving finite element problems may be available to solve this problem.
Analysis and Results
The concentration of substrate reaching the sink over a 0.2 second time-frame was measured at multiple distances. The largest distance between enzymes was 1500 angstroms, with the smallest distance measured being 50 angstroms. The distance between enzymes was reduced by 50-angstrom increments to measure the effect of distance on final product yield. The model predicts a 5-fold improvement in IAA production when the enzymes are clustered together using prefoldin, compared to when they are free in solution (Figure 3).
Figure 3: From left to right, each graph represents a smaller distance between the two enzymes. The uppermost graph presents the concentration-time curve when the enzymes are separated by a distance of 50 angstroms, with the lowermost graph displaying the concentration-time curve for a distance of 1500 angstroms.
The graphs show a faster achievement of steady-state IAA production when the two enzymes are clustered together, and higher overall concentration of the second substrate at enzyme B.
By integrating under each curve, we calculated the amount of product formed at each distance over the same time period (Figure 4). As stated above, this calculation assumed perfect substrate to product conversion.
Figure 4: Each point represents the amount of IAA produced in micromoles when the enzymes are separated by a certain distance in angstroms.
We were less concerned with the exact values of product generated by the model, and more with the relationship between distance and amount of product formed. From this result the relationship between distance and product formation was made clear; product formation increases non-linearly as the separation between enzymes is reduced.
MATLAB
The MATLAB PDE toolbox was used to solve the system. The files for generating the geometry and solving the system are available on GitHub, with a readme.txt file supplied to describe how the parameters can be manipulated.
Reflections on the Model
Overall, the model served to effectively support our initial hypothesis that reduced distance between enzymes would increase product yield. We believed the boundary conditions and assumptions did not limit the possibility of validating the model against experimental data. However, we did not reach the stage in building our prefoldin scaffold where such validation was possible. Thus the model served mainly to verify our design process, and could not be modified based on experimental results. Furthermore we believe there is better-suited software specialized for designing molecular meshes, and solving numerical problems, such as the FeNics toolkit. However due to inexperience in computational mathematics we were not able to build our model in a more complex environment.