This is a one-site modification cycle with two competing kinases.
i1 : oneSiteModificationC() o1 = S_0+E_k-->X_k X_k-->S_0+E_k X_k-->E_k+S_1 S_1+F-->Y Y-->S_1+F Y-->S_0+F o1 : ReactionNetwork
To rename species, use substitute(ReactionNetwork,List).
The object oneSiteModificationC is a method function.