The motions of water molecules, the acyl moiety, the catalytic triad, and the oxyanion binding site of acyl-chymotrypsin were studied by means of a stochastic boundary molecular dynamics simulation. A water molecule that could provide the nucleophilic OH- for the deacylation stage of the catalysis was found to be trapped between the imidazole ring of His-57 and the carbonyl carbon of the acyl group. It makes a hydrogen bond with the N epsilon 2 of His-57 and is held in place through a network of hydrogen-bonded water molecules in the active site. The water molecule was found as close as 2.8 A to the carbonyl carbon. This appears to be due to the constraints imposed by nonbonded interaction in the active site. Configurations were found in which one hydrogen of the trapped water shared a bifurcated hydrogen bond with His-57-N epsilon 2 and Ser-195-O gamma, with the water oxygen very close to the carbonyl carbon. The existence of such a water molecule suggests that large movement of the His-57 imidazole ring between positions suitable for providing general-base catalyzed assistance and for providing general-acid catalyzed assistance may not be required during the reaction. The simulation indicates that the side chains of residues involved in catalysis (i.e., His-57, Ser-195, and Asp-102) are significantly less flexible than other side chains in the protein. The 40% reduction in rms fluctuations is consistent with a comparable reduction calculated from the temperature factors obtained in the X-ray crystallographic data of gamma-chymotrypsin. The greater rigidity of active site residues seems to result from interconnected hydrogen bonding networks among the residues and between the residues and the solvent water in the active site.