<<xlr8r.m;

model={{ES→Enz+S1,Flux[0,v1,1,1,0]},{S1→S2,Flux[-∞,v2,∞,1,0]},{Enz+S2→ES,Flux[-∞,v3,∞,1,1]}}

interpret[model]

StoichiometryMatrix[model]//TableForm

OptimizeFluxModel[model]