The initial and boundary value problem in enzyme reactions mechanism for inhabitation process is discussed. Approximate analytical expressions for the concentrations of substrate and product are presented. Approximate analytical solutions of non-linear reaction equations containing non-linear terms related to enzymatic reaction mechanism are solved using Homotopy perturbation method. The relevant analytical expression for the substrate and product concentration profiles is discussed in terms of dimensionless reaction diffusion parameters α, β, γE, and γS. Numerical solution is also obtained using Matlab program. Our analytical expression compared with numerical estimation and good agreement is noted