A mathematical model of biotransformation of D-methionine into L-methionine in the cascade of the enzymes such as, D-amino acid oxidase (D-AAO), L-phenylalanine dehydrogenase (L-PheDH) and formate dehydrogenase (FDH) is discussed. The model is based on a system of coupled nonlinear reaction equations under non steady-state conditions for biochemical reactions occurring in the batch reactor that describes the substrate and product concentration within the catalyst. Simple analytical expressions for the concentration of substrate and product have been derived for all values of reaction parameters using the new homotopy perturbation method (NHPM). Enzyme reaction rate in terms of concentration and kinetic parameters are also reported. The analytical results are also compared with experimental and numerical ones and a good agreement is obtained. The graphical procedure for estimating the kinetic parameters is also reported.