A Quantum Mechanics / Molecular Mechanics Study of the Protein-Ligand Interaction of two potent inhibitors of human O-GlcNAcase: PUGNAc and NAG-thiazoline

O-glycoprotein 2-acetamino-2-deoxy-β-D-glucopyranosidase (O-GlcNAcase) hydrolyzes 2-acetamido-2-deoxy-β-D-glucopyranose (O-GlcNAc) residues of serine/threonine residues of modified proteins. O-GlcNAc is present in many intracellular proteins and appears to have a role in the etiology of several diseases including cancer, Alzheimer’s disease, and type II diabetes. In this work, we have carried out molecular dynamics (MD) simulations using a hybrid Quantum Mechanics / Molecular Mechanics (QM/MM) approach, to determine the binding of two potent inhibitors, PUGNAc and NAG, with a bacterial O-GlcNAcase. The results of these simulations show that Asp-401, Asp-298 and Asp-297 residues play an important role in the protein-inhibitor interactions. These results might be useful to design compounds with more interesting inhibitory activity on the basis of its three-dimensional structur

QM/MM study of thymidylate synthase: enzymatic motions and the temperature dependence of the rate limiting step

Thymidylate synthase (TS) is an enzyme that catalyzes a complex cascade of reactions. A theoretical study of the reduction of an exocyclic methylene intermediate by hydride transfer from the 6S position of 5,6,7,8-tetrahydrofolate (H4folate), to form 2′-deoxyuridine 5′-monophosphate (dTMP) and 7,8-dihydrofolate (H2folate), has been carried out using hybrid quantum mechanics/molecular mechanics methods. This step is of special interest because it is the rate-limiting step of the reaction catalyzed by TS. The acceptor of this hydride is an intermediate that is covalently bound to the enzyme via a thioether bond to an overall conserved active site cysteine residue (Cys146 in Escherichia coli). Heretofore, whether the hydride transfer precedes the thiol abstraction that releases the product from the enzyme or whether these two processes are concerted has been an open question. We have examined this step in terms of free energy surfaces obtained at the same temperatures we previously used in experimental studies of this mechanistic step (273−313 K). Analysis of the results reveals that substantial features of the reaction and the nature of the H-transfer seem to be temperature independent, in agreement with our experimental data. The findings also indicate that the hydride transfer and the scission of Cys146 take place in a concerted but asynchronous fashion. This 1,3-SN2 substitution is assisted by arginine 166 and several other arginine residues in the active site that polarize the carbon−sulfur bond and stabilize the charge transferred from cofactor to substrate. Finally, the simulation elucidates the molecular details of the enzyme’s motion that brings the system to its transition state and, in accordance with the experimental data, indicates that this “tunneling ready” conformation is temperature independent

Theoretical study of the temperature dependence of dynamic effects in thymidylate synthase

A theoretical study of the temperature dependence of dynamic effects in the rate limiting step of the reaction catalyzed by thymidylate synthase is presented in this paper. From hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) optimizations of transition state structures within a fully flexible molecular model, free downhill molecular dynamics trajectories have been performed at four different temperatures. The analysis of the reactive and non-reactive trajectories in the enzyme environment has allowed us to study the geometric and electronic coupling between the substrate, the cofactor and the protein. The results show how the contribution of dynamic effects to the rate enhancement measured by the transmission coefficients is, at the four studied temperatures, negligible. Nevertheless, the rare event trajectories performed have shown how the hydride transfer and the scission of the conserved active site cysteine residue (Cys146 in E. coli) take place in a concerted but asynchronous way; the latter takes place once the transfer has occurred. The analysis of the dynamics of the protein reveals also how the relative movements of some amino acids, especially Arg166, and a water molecule, promotes the departure of the Cys146 from the dUMP. Finally, it seems that the protein environment creates an almost invariant electric field in the active site of the protein that stabilizes the transition state of the reaction, thus reducing the free energy barrier