Theoretical Perspectives on the Reaction Mechanism of Serine Proteases: The Reaction Free Energy Profiles of the Acylation Process

Citations Over TimeTop 11% of 2003 papers

Abstract

The reaction mechanism of serine proteases (trypsin), which catalyze peptide hydrolysis, is studied theoretically by ab initio QM/MM electronic structure calculations combined with Molecular Dynamics-Free Energy Perturbation calculations. We have calculated the entire reaction free energy profiles of the first reaction step of this enzyme (acylation process). The present calculations show that the rate-determining step of the acylation is the formation of the tetrahedral intermediate, and the breakdown of this intermediate has a small energy barrier. The calculated activation free energy for the acylation is approximately 17.8 kcal/mol at QM/MM MP2/(aug)-cc-pVDZ//HF/6-31(+)G/AMBER level, and this reaction is an exothermic process. MD simulations of the enzyme-substrate (ES) complex and the free enzyme in aqueous phase show that the substrate binding induces slight conformational changes around the active site, which favor the alignment of the reactive fragments (His57, Asp102, and Ser195) together in a reactive orientation. It is also shown that the proton transfer from Ser195 to His57 and the nucleophilic attack of Ser195 to the carbonyl carbon of the scissile bond of the substrate occur in a concerted manner. In this reaction, protein environment plays a crucial role to lowering the activation free energy by stabilizing the tetrahedral intermediate compared to the ES complex. The polarization energy calculations show that the enzyme active site is in a very polar environment because of the polar main chain contributions of protein. Also, the ground-state destabilization effect (steric strain) is not a major catalytic factor. The most important catalytic factor of stabilizing the tetrahedral intermediate is the electrostatic interaction between the active site and particular regions of protein: the main chain NH groups in Gly193 and Ser195 (so-called oxyanion hole region) stabilize negative charge generated on the carbonyl oxygen of the scissile bond, and the main chain carbonyl groups in Ile212 approximately Ser214 stabilize a positive charge generated on the imidazole ring of His57.

Related Papers

• → Ab Initio QM/MM Study of Class A β-Lactamase Acylation: Dual Participation of Glu166 and Lys73 in a Concerted Base Promotion of Ser70(2005)156 cited
• → The Impact of Substituents on the Transition States of S_N2 and E2 Reactions in Aliphatic and Vinylic Systems: Remarkably Facile Vinylic Eliminations(2012)21 cited
• → The QM/MM Molecular Dynamics and Free Energy Simulations of the Acylation Reaction Catalyzed by the Serine-Carboxyl Peptidase Kumamolisin-As(2007)25 cited
• → Kinetic Study on the Acylation Step of α-Chymotrypsin-Catalyzed Hydrolysis of Acylimidazole(1980)5 cited
• → Mechanistic Investigation on the Reaction of O− with CH3CN Using Density Functional Theory(2010)1 cited