(eng): Our study presents an application of one special technique, further called the Bohl- Marek formulation, related to the mathematical modeling of biochemical networks with mass conservation properties. Nonlinear cooperative systems possessing certain conservation laws, arising fre- quently in biology and chemistry, were studied by Erich Bohl and Ivo Marek in different papers [1], [2]. The conservation law guarantees a lot of mathematical properties, e.g. the existence and uniqueness results, the positivity of the solutions (for positive inputs), periodicity and controllability, or proposing a stability theorem. These properties are consequence of the fact that the underlying nonlinear dynamical systems describing a class of biochemical networks are built up of linear evolutions with negative M-matrices whose entries depend on the dynamical variables of the other subsystems involved. Thus, the nonlinearity of the whole system is created via this dependence, i.e., the matrix of the complete system is blockwise diagonal. On a case study, namely the Michaelis-Menten enzyme-catalyzed reaction with a substrate transport chain [3], we show how to transform the system of nonlinear ODEs into a set of smaller, quasi-linear subsystems of ODEs with negative M- matrices. For the two model formulations, the classical nonlinear formulation and the quasi-linear Bohl-Marek formulation, we determine and compare the results and show the computational advantages of the latter formulation.