most biochemical reactions in living cells are open system interacting with environment through chemostats. at a mesoscopic scale, the number of  each species in those biochemical reactions can be modeled by the random time-changed poisson processes. to characterize the macroscopic behaviors in the large volume limit, the law of large number in path space determines a mean-field limit nonlinear kurtz ode, while the wkb expansion yields a hamilton-jacobi equation and the corresponding lagrangian gives the good rate function in the large deviation principle. a parametric variation principle can be formulated to compute the reaction paths. we propose a gauge-symmetry criteria for a class of non-equilibrium chemical reactions including  enzyme reactions, which identifies a new concept of balance within the same reaction vector due to flux grouping degeneracy. with this criteria,  we (i) formulate an onsager-type gradient flow structure in terms of the energy landscape given by a steady solution to the hamilton-jacobi equation;  (ii) find transition paths between multiple non-equilibrium steady states (rare events in biochemical reactions).  we illustrate this idea through a bistable catalysis  reaction. in contrast to the standard diffusion approximations via kramers-moyal expansion, a new drift-diffusion approximation sharing the same gauge-symmetry is constructed based on the onsager-type gradient flow formulation to compute the correct energy barrier.