The kinetic behavior of the ␣-chymotrypsincatalyzed hydrolysis of the two p-nitroanilide substrates succinyl-L-Ala-L-Ala-L-Pro-L-Phe-p-nitroanilide (Suc-Ala-Ala-Pro-Phe-pNA) and benzoyl-L-Tyr-p-nitroanilide (Bz-Tyr-pNA) was modeled and simulated for two different systems, namely for an aqueous solution and for a vesicle system, which was composed of phospholipid vesicles containing entrapped ␣-chymotrypsin. In the case of the vesicles, the substrate was added to the bulk, exovesicular aqueous phase. The experimentally determined time-dependence of product (p-nitroaniline) formation was modeled by considering the kinetic behavior of the enzyme and-in the case of vesicles-the substrate permeability across the bilayer membrane. In aqueous solution-without vesicles-the kinetic constants k cat and K S (respectively K M ) were determined from fitting the model to experimental data of batch product concentration-time curves. The results were in good agreement with the corresponding values obtained from initial velocity measurements. For the vesicle system, using the phospholipid 1-palmitoyl-2-oleoyl-snglycero-3-phosphocholine (POPC), simulation showed that the substrate permeation across the bilayer was rate limiting. Using experimental data, we could obtain the substrate permeability coefficient for Bz-Tyr-pNA by parametric fitting as 2.45 × 10 -7 cm/s.