The basic enzyme reaction boundary value problem is described and approximate expressions for substrate and product concentrations are given. This model was initially designed using the classical differential equations and it is extended to the Caputo fractional derivative (FDE’s) of order μ. Non-linear reaction equations with a non-linear term related to enzymatic reaction can be approximated and analytically solved using the Homotopy Perturbation method. Dimensionless reaction diffusion parameters ε^μ,k^μ,λ^μ are used to discuss the relevant analytical solutions for the substrate, enzyme, substrate-enzyme, and product concentration profiles.