On Discovering Low Order Models in Biochemical Reaction Kinetics