Zhan_2005_Biophys.J_89_3863

The geometries of the transition states, intermediates, and prereactive enzyme-substrate complex and the corresponding energy barriers have been determined by performing hybrid quantum mechanical/molecular mechanical (QM/MM) calculations on butyrylcholinesterase (BChE)-catalyzed hydrolysis of (-)- and (+)-cocaine. The energy barriers were evaluated by performing QM/MM calculations with the QM method at the MP2/6-31+G* level and the MM method using the AMBER force field. These calculations allow us to account for the protein environmental effects on the transition states and energy barriers of these enzymatic reactions, showing remarkable effects of the protein environment on intermolecular hydrogen bonding (with an oxyanion hole), which is crucial for the transition state stabilization and, therefore, on the energy barriers. The calculated energy barriers are consistent with available experimental kinetic data. The highest barrier calculated for BChE-catalyzed hydrolysis of (-)- and (+)-cocaine is associated with the third reaction step, but the energy barrier calculated for the first step is close to the highest and is so sensitive to the protein environment that the first reaction step can be rate determining for (-)-cocaine hydrolysis catalyzed by a BChE mutant. The computational results provide valuable insights into future design of BChE mutants with a higher catalytic activity for (-)-cocaine.