A combined quantum mechanical (QM) and molecular mechanical (MM) method was used to study the nucleophilic addition reaction of methanethiolate to N-methylacetamide (NMA) in the gas phase and aqueous solution. At the B3LYP/aug-cc-pVDZ//HF/6-31+G(d) level, the ion-dipole complex was found to be the global minimum on the potential energy surface in the gas phase with a binding energy of 21.2 kcal/mol. The complex has a C-S distance of 4.33 Å, and no stabilized tetrahedral intermediate was located. The computed potential of mean force in water shows that solvent effects stabilize the reactants over the tetrahedral adduct by 36.5 kcal/mol, and that the tetrahedral intermediate does not exist for the present reaction in water. The present study provides an initial step for modeling the cysteine protease hydrolysis reactions in enzymes.
- Combined QM/MM simulations
- Monte Carlo simulation
- Nucleophilic addition reaction
- Potential of mean force
- Solvent effects