Assessment of Reversibility for Covalent Cysteine Protease Inhibitors Using Quantum Mechanics/Molecular Mechanics Free Energy Surfaces.

We have used molecular dynamics (MD) simulations with hybrid quantum mechanics/molecular mechanics (QM/MM) potentials to investigate the reaction mechanism for covalent inhibition of cathepsin K and assess the reversibility of inhibition. The computed free energy profiles suggest that a nucleophilic attack by the catalytic cysteine on the inhibitor warhead and proton transfer from the catalytic histidine occur in a concerted manner. The results indicate that the reaction is more strongly exergonic for the alkyne-based inhibitors, which bind irreversibly to cathepsin K, than for the nitrile-based inhibitor odanacatib, which binds reversibly. Gas-phase energies were also calculated for the addition of methanethiol to structural prototypes for a number of warheads of interest in cysteine protease inhibitor design in order to assess electrophilicity. The approaches presented in this study are particularly applicable to assessment of novel warheads, and computed transition state geometries can be incorporated into molecular models for covalent docking.

Computational study of conformational changes in human 3-hydroxy-3-methylglutaryl coenzyme reductase induced by substrate binding.

The enzyme 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) is mainly involved in the regulation of cholesterol biosynthesis. HMGR catalyses the reduction of 3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) to mevalonate at the expense of two NADPH molecules in a two-step reversible reaction. In the present study, we constructed a model of human HMGR (hHMGR) to explore the conformational changes of HMGR in complex with HMG-CoA and NADPH. In addition, we analysed the complete sequence of the Flap domain using molecular dynamics (MD) simulations and principal component analysis (PCA). The simulations revealed that the Flap domain plays an important role in catalytic site activation and substrate binding. The apo form of hHMGR remained in an open state, while a substrate-induced closure of the Flap domain was observed for holo hHMGR. Our study also demonstrated that the phosphorylation of Ser872 induces significant conformational changes in the Flap domain that lead to a complete closure of the active site, suggesting three principal conformations for the first stage of hHMGR catalysis. Our results were consistent with previous proposed models for the catalytic mechanism of hHMGR. Communicated by Ramaswamy H. Sarma.