
Multiple BinsFrom MathBioIntroductionBiological systems of varying complexities are frequently simplified by division into bins, generally based on spatial domains. These bins are frequently connected by diffusion reactions, as products from one reaction (in one bin) flow to a different spatial domain (a different bin) and are taken up as a reactant in a separate reaction. We are particularly interested in the composition of a stoichiometricdiffusion matrix for multiple bins. As an example, we consider our open MichaelisMenten example in the context of bins. Points 1, 0, and +1 can be in a sense taken as representative of three separate bins, where the 1 bin "holds" the substrate, which flows into the 0 bin. The 0 bin "holds" the main MichaelisMenten reactions, involving the substrate and enzyme combining into a complex and subsequently forming a product. Finally, the +1 bin "holds" the product leaving the reaction domain to fulfill whatever nefarious purposes it might have. It is worth noting that our method of incorporating diffusion into our stoichiometric matrix using the second derivative and Fick's Laws involves the diffusing substance being "held" at a fixed distance h, as opposed to the biologically realistic case of various distances radially around the fixed enzyme. In that sense, the 1 bin "holding" the substrate and the +1 bin "holding" the product might be better described as existing around the 0 reaction bin, as opposed to existing to the left and to the right, respectively. Additionally, note that the reaction gradient method of approximating diffusion as a flux automatically accounts for "net" diffusion between bins, including in the stoichiometric matrix. For example, the substrate diffusion term represents diffusion into the reaction stoichiometry when the amount of substrate in the 1 bin is greater, and represents diffusion out of the reaction stoichiometry when the amount of substrate in the 0 bin is greater.
Two Connected MichaelisMenten ReactionsIn order to develop an understanding of the structure of a stoichiometricdiffusion matrix in the case of multiple bins, we begin by examining the simple case of two connected MichaelisMenten reactions, where the product of the first reaction forms the substrate for a second substrateenzyme reaction. Such a cascade is relevant to our interest in the complement cascade found in the human innate immune system. We can then consider the system as being composed of four distinct bins: the 1 bin again holds the substrate which reacts in the first MichaelisMenten reaction, which we denote . The enzyme and the complex for the first reaction are denoted and , respectively, and will be assumed to remain native to the 0 bin, which will hold the first MichaelisMenten reaction. The product produced in this scenario is denoted and flows into the +1 bin, where the second MichaelisMenten reaction occurs. The enzyme and the complex for the second reaction are denoted and , respectively, and will be assumed to remain native to the +1 bin. Finally, the +2 bin will hold the outflow of the product of the second MichaelisMenten reaction, denoted . We keep the same reaction names for the first MichaelisMenten reaction and assign as the corresponding reactions for the second MichaelisMenten reaction. We formulate our system of differential equations:
Using the same diffusion approximation as mentioned previously, we have generated a size matrix such that
,
where
Considering the Reaction Space Only
Right null space, rank r = 6 