(eng): This study presents an application of one special technique, further called as Bohl-Marek decomposition, related to the mathematical modeling of biochemical networks with mass conservation properties. We continue in direction of papers devoted to inverse problems of parameter estimation for mathematical models describing the drug-induced enzyme production networks [3]. However, being aware of the complexity of general physiologically based pharmacokinetic (PBPK) models, here we focus on the case of enzyme-catalyzed reactions with a substrate transport chain [5]. Although our ultimate goal is to develop a reliable method for fitting the model parameters to given experimental data, here we study certain numerical issues within the framework of optimal experimental design [6]. Before starting an experiment on a real biochemical network, we formulate an optimization problem aiming to maximize the information content of the corresponding experiment. For the above-sketched optimization problem, the computational costs related to the two formulations of the same biochemical network, being (i) the classical formulation x˙(t) = Ax(t) + b(t) and (ii) the 'quasi-linear' Bohl-Marek formulation x˙M(t) = M(x(t)) xM(t), can be determined and compared.