The stability of peptide bonds is a critical aspect of biological chemistry and therapeutic protein applications. Recent studies found elevated nonenzymatic hydrolysis in the hinge region of antibody molecules, but no mechanism was identified. As a first step in providing a mechanistic interpretation, this computational study examines the rate-determining step of the hydrolytic reaction of a peptide bond under acidic pH by a path sampling technique using a model compound N-MAA. Most previous computational studies did not include explicit water molecules, whose effects are significant in solution chemistry, nor did they provide a dynamic picture for the reaction process in aqueous conditions. Because no single trajectory can be used to describe the reaction dynamics due to fluctuations at finite temperatures, a variant version of the transition path sampling technique, the aimless shooting algorithm, was used to sample dynamic trajectories and to generate an ensemble of transition trajectories according to their statistical weights in the trajectory space. Each trajectory was computed as the time evolution of the molecular system using the Car Parrinello molecular dynamics technique. The likelihood maximization procedure and its modification were used in extracting dynamically relevant degrees of freedom in the system, and approximations of the reaction coordinate were compared. Its low log-likelihood score and poor p(B) histogram indicate that the C-O distance previously assumed as the reaction coordinate for the rate-determining step is inadequate in describing the dynamics of the reaction. More than one order parameter in a candidate set including millions of geometric quantities was required to produce a convergent reaction coordinate model; its involvement of many degrees of freedom suggests that this hydrolytic reaction step is very complex. In addition to affecting atoms directly involved in bond-making and -breaking processes, the water network also has determining effects on the hydrolytic reaction, a fact that is manifest in the expression of the one-dimensional best-ranked reaction coordinate, which includes three geometric quantities. The p(B) histograms were computed to verify the results of the likelihood maximization and to verify the accuracy of approximation to the "true" reaction coordinate.