
Looking for a SteadyState Flux DistributionFrom MathBioIntroductionAs seen in our diffusion approximation, a biologically realistic steadystate flux distribution does not exist for our C5a production model. We are interested studying cases where a steadystate flux distribution arises when diffusion is considered in addition to reaction stoichiometry, so we now turn to the canonical MichaelisMenten model and add diffusion using our "flux" method. Classical MichaelisMenten Model With DiffusionWe will begin by looking at the classical model considered earlier, in which it is assumed that substrate flows into some domain containing enzyme to form the substrateenzyme complex. The product formed leaves the domain, while the enzyme remains internal to the spatial domain. We will continue to consider diffusion along a line as in our C5a production model and reaction occurring at the surface of the enzyme (bound by substrate). We approximate the reaction space as a single point and simplify the system so that it effectively involves the substrate flowing into this point from the left, binding with enzyme at our single point, and then undergoing a conformational change to produce a product, which flows out of the domain to the right. We label the reaction space as point 0, the points to the left of the reaction space where the substrate diffuses in as ..., 3, 2, 1, and the points to the right of the reaction space where the product leaves as +1, +2, +3, ... However, we will begin by approximating this system at three discrete points: 1, some distance to the left of the reaction space, 0, at the reaction space, and +1, some distance to the right of the reaction space. Then our system of differential equations is:
Where is the diffusion rate of the substrate into the reaction and is the diffusion rate of the product out of the reaction. Note that we assume some nonzero initial concentration at point 1 that flows into the reaction space, which has zero initial concentration of the substrate. We also assume that the enzyme and complex remain internal to the reaction space. We now construct the stoichiometricdiffusion matrix. As seen in the C5a production case, note that the diffusion is considered as a gradient, such that the diffusion flux terms seen in S, 0 and S, 1 are equal but opposite. Thus, we have created an matrix such that
,
where
We now conduct our SVD analysis to identify the steadystate flux distribution. As seen in our C5a production example, our results for a finite number of points at wellchosen locations can be generalized to the real infinite case. The rank , so the dimension of the right null space is 1. By SVD, a basis for this space is given by: This corresponds to a steadystate arising when , or when the substrateenzyme binding reaction is equal and opposite to the complex decomposition with all other reaction fluxes zero. Of course, this is not a biologically realistic scenario. A similar result had been found for the closed case previously. Considering the Reaction Space Only However, note that this SVD method takes into account steadystates for the complex at point 1 and for the product at point +1, which we are not necessarily concerned with. We will see what happens when these two differential equations are removed from our stoichiometricdiffusion matrix:
We now have a matrix with rank , so the dimension of the right null space is 2, given by the following orthonormal basis: Any linear combination of these vectors will give a steadystate flux distribution; our goal is now to find biological relevance from these results. We begin by trying , as this will give us all positive entries. And a mathematical check shows that this flux distribution indeed gives rise to a steadystate flux distribution. The biological significance of this particular reaction stoichiometry is uncertain, but would indeed be biologically realistic. It is worth noting, however, that as part of this steadystate we have substrate flowing into the reaction at the same rate that product flows out of it, as expected. 