We developed a kinetic model that describes a heterogeneous reaction system for the production of D-p-hydroxyphenylglycine from D,L-p-hydroxyphenylhydantoin using D-hydantoinase of Bacillus stearothermophilus SD1 and N-carbamoylase of Agrobacterium tumefaciens NRRL B11291. As a biocatalyst, whole cells with separately or co-expressed enzymes were used. The reaction system involves dissolution of substrate particles, enzymatic conversion, racemization of the L-form substrate, and transfer of the dissolved substrate, intermediate, and product through the cell membrane. Because the two enzymes have, different pH optimum, kinetic parameters were evaluated at different pH for the reaction systems. The model was simulated using the kinetic parameters and compared with experimental data, and it was found that the kinetic model-well describes the behavior of the reaction systems using whole cells with separately and co-expressed enzymes. Factors affecting the kinetics of the reaction systems were analyzed on the basis of the kinetic model. In the reaction system with separately expressed enzymes, racemization rate and transport of the reaction intermediate (N-carbamoyl-D-p-hydroxyphenylglycine) were revealed to be the limiting factors at neutral pH, resulting in accumulation of intermediate in the reaction medium. At alkaline condition, on the other hand, inhibition of N-carbamoylase by ammonia was severe and thereby the reaction rate significantly reduced. In the co-expressed enzyme system, accumulation of intermediate was negligible in the reaction medium, and the improved performance was observed compared to that with separately expressed enzymes. The present model might be applied for the optimization and development of the reaction system using two sequential enzymes.