  This paper studies a system of Ordinary Differential Equations modeling a
chemical reaction network and derives from it a simulation tool mimicking Loss
of Function and Gain of Function mutations found in cancer cells. More
specifically, from a theoretical perspective, our approach focuses on the
determination of moiety conservation laws for the system and their relation
with the corresponding stoichiometric surfaces. Then we show that Loss of
Function mutations can be implemented in the model via modification of the
initial conditions in the system, while Gain of Function mutations can be
implemented by eliminating specific reactions. Finally, the model is utilized
to examine in detail the G1-S phase of a colorectal cancer cell.
