In this paper the complex dynamics of a smallest biochemical system model in three-dimensional systems with the reaction scheme. This model is described by a system of three nonlinear ordinary differential equations with five positive real parameters, are analyzed and studied. We present a thorough analysis of their invariant algebraic surfaces and exponential factors and investigate the integrability and nonintegrabilty of this model. Particularly, we show the non-existence of polynomial, rational, Darboux and local analytic first integrals in a neighborhood of the equilibrium. Moreover, we prove that, the model is not integrable in the sense of Bogoyavlensky in the class of rational functions.