**This page is under construction.**

## Introduction

The closed Michaelis-Menten model can easily be modified to allow for a more biologically realistic open system, by introducing a substrate input rate and a product outflow flux , such that:

This simple Michaelis-Menten system is being studied because its behavior in response to certain inputs offers a great deal of insight into reactant behavior in much more complex pathways, including the complement system of the human immune system, which I will study later as a biological application.

My focus will be on the open Michaelis-Menten system, in which traditional analysis of the closed Michaelis-Menten system will prove important. In this case, the stoichiometric matrix for the open reaction (as depicted in the introduction) is given by:

Using the Law of Mass Action, this corresponds to the following set of differential equations:

## Maximum Allowable Reaction Rate

One particular question we would like to answer is the maximum allowable substrate input rate such that the substrate does not grow out of control over time, but instead approaches a steady state. Implicit in this requirement is that the complex and product both approach a steady-state as well.

By taking the (right) null space of the stoichiometric matrix, we can find spanning vectors for the steady state flux vector . The rank of the matrix is 3, so the dimension of must be 2. The two biologically significant spanning vectors are and , which correspond to a pathway through the system and the reversible reaction, respectively. In essence, these vectors tell us that the steady-state flux balance for the system is as follows:

This makes sense: the system will have reached equilibrium if the fluxes are balanced as such. Thus, it follows that the maximum substrate input rate *below but not equal to which* the substrate will approach a steady-state and not grow out of control is:

.

The flux is the most logical to use, as we can make use of the conservation relationship to note that cannot approach a steady state if the system has all complex and no enzyme in the conservation relationship at any time; that is, if . In this case, the enzyme will not be able to convert substrate to product quickly enough, so the substrate input rate must be set at a value *lower* than for the system to approach a steady-state.

**Steady-State Concentrations**

Algebraic expressions for the steady-state concentrations may be useful. By setting all four differential equations in the open system to 0 to simulate the steady-state, we can algebraically solve for all four steady-state concentrations.