  The estimation of the kinetic parameters requires the careful design of
experiments under a constrained set of conditions. Many estimates reported in
the literature incorporate protocols that leverage simplified mathematical
models known as quasi-steady-state reductions. Such reductions often - but not
always - emerge as the result of a singular perturbation scenario. However, the
utilization of the singular perturbation reduction method requires knowledge of
a dimensionless parameter, $\varepsilon$, that is proportional to the ratio of
the reaction's fast and slow timescales. Using techniques from differential
equations, Fenichel theory, and center manifold theory, we derive the
appropriate $\varepsilon$ whose magnitude regulates the validity of the
quasi-steady-state reduction employed in the reported experimental procedures
for intermolecular autocatalytic zymogen activation reaction. Although the
model equations are two-dimensional, the fast/slow dynamics are rich. The phase
plane exhibits a dynamic transcritical bifurcation point in a particular
singular limit. The existence of such a bifurcation is relevant, because the
critical manifold losses normal hyperbolicity and classical Fenichel theory is
inapplicable. Furthermore, we show that in some cases chemical reversibility
can be interpreted dynamically as an imperfection, since the presence of
reversibility can destroy the bifurcation structure present in the singular
limit. We show that the reduction method by which QSS reductions are justified
can depend on the path taken in parameter space. Specifically, we show that the
standard quasi-steady-state reduction for this reaction is justifiable by
center manifold theory in one limit, and via Fenichel theory in a different
limit.
