Dynamic mass balances on each metabolite
Time derivatives of metabolite concentrations are linear combination of the reaction rates. The reaction rates are non-linear functions of the metabolite concentrations (typically from in vitro kinetics).
Where vj is the jth reaction rate, b is the transport rate vector,
Sij is the “Stoichiometric matrix” = moles of metabolite i produced in reaction j