Molecular Dynamics Simulations of the Catalytic Pathway of a Cysteine Protease: A Combined QM/MM Study of Human Cathepsin K

Citations Over TimeTop 15% of 2007 papers

Abstract

Molecular dynamics simulations using a combined QM/MM potential have been performed to study the catalytic mechanism of human cathepsin K, a member of the papain family of cysteine proteases. We have determined the two-dimensional free energy surfaces of both acylation and deacylation steps to characterize the reaction mechanism. These free energy profiles show that the acylation step is rate limiting with a barrier height of 19.8 kcal/mol in human cathepsin K and of 29.3 kcal/mol in aqueous solution. The free energy of activation for the deacylation step is 16.7 kcal/mol in cathepsin K and 17.8 kcal/mol in aqueous solution. The reduction of free energy barrier is achieved by stabilization of the oxyanion in the transition state. Interestingly, although the "oxyanion hole" has been formed in the Michaelis complex, the amide units do not donate hydrogen bonds directly to the carbonyl oxygen of the substrate, but they stabilize the thiolate anion nucleophile. Hydrogen-bonding interactions are induced as the substrate amide group approaches the nucleophile, moving more than 2 A and placing the oxyanion in contact with Gln19 and the backbone amide of Cys25. The hydrolysis of peptide substrate shares a common mechanism both for the catalyzed reaction in human cathepsin K and for the uncatalyzed reaction in water. Overall, the nucleophilic attack by Cys25 thiolate and the proton-transfer reaction from His162 to the amide nitrogen are highly coupled, whereas a tetrahedral intermediate is formed along the nucleophilic reaction pathway.

Related Papers

• → Ab initio and density functional theory studies of the catalytic mechanism for ester hydrolysis in serine hydrolases(1998)56 cited
• → Mechanism of action of cysteine proteinases: oxyanion binding site is not essential in the hydrolysis of specific substrates(1985)47 cited
• → Role of catalytic residues in the formation of a tetrahedral adduct in the acylation reaction of bovine β-trypsin(1984)27 cited
• → Quantum mechanical/molecular mechanical study of three stationary points along the deacylation step of the catalytic mechanism of elastase(2001)12 cited
• Energetics of Catalytic Reaction of Acetylcholinesterase (AChE) with Acetylcholine (ACh): Role of the Oxyanion Hole #(2003)