Search

Human aromatase (CYP19A1) aromatizes the androgens to form estrogens via a three-step oxidative process. The estrogens are necessary in humans, mainly in women, because of the role they play in sexual and reproductive development. However, these also are involved in the development and growth of hormone-dependent breast cancer. Therefore, inhibition of the enzyme aromatase, by means of drugs known as aromatase inhibitors, is the frontline therapy for these types of cancers. Exemestane is a suicidal third-generation inhibitor of aromatase, currently used in breast cancer treatment. In this study, the hydroxylation of exemestane catalyzed by aromatase has been studied by means of hybrid QM/MM methods. The Free Energy Perturbation calculations provided a free energy of activation for the hydrogen abstraction step (rate-limiting step) of 17 kcal/mol. The results reveal that the hydroxylation of exemestane is not the inhibition stage, suggesting a possible competitive mechanism between the inhibitor and the natural substrate androstenedione in the first catalytic subcycle of the enzyme. Furthermore, the analysis of the interaction energy for the substrate and the cofactor in the active site shows that the role of the enzymatic environment during this reaction consists of a transition state stabilization by means of electrostatic effects. [ABSTRACT FROM AUTHOR]

A theoretical study of dynamic effects on the rate-limiting step of the thymidylate synthase catalyzed reaction has been carried out by means of Grote-Hynes theory, successfully predicting the values of the recrossing effects for a chemical reaction that involves the transfer of a classical light particle. The transmission coefficients, obtained at 278, 293, 303, and 313 K, are almost invariant and in all cases far from unity, revealing a significant coupling of the environment motions and the reaction coordinate. Nevertheless, their energetic contribution to the activation free energy represents less than 0.50 kcal/mol for each of the four tested temperatures. Calculation of the transmission coefficient for the isotopically labeled hydride transfer has rendered almost the same values, in agreement with the experimentally observed temperature-independent KIEs. Fourier transform of the time-dependent friction kernel at these four temperatures has allowed obtaining the transition-state friction spectra, which present very small dependence with temperature. Their analysis has led to the identification of some key vibrational modes governing the coupling between the reaction coordinate and the protein environment, thus identifying the relevant motions in the active site and obtaining a full picture of the role of each amino acid.

The reaction mechanism of phosphate monoester hydrolysis in alkaline phosphatase is analyzed by means of hybrid QM/MM simulations. A recently developed semiempirical Hamiltonian, AM1/d-PhoT, which takes into account the d orbitals on the phosphorus atom, has been employed. The reaction mechanism obtained is either associative or dissociative, depending on the size of the QM subsystem. The results are rationalized on the basis of the degree of charge transfer from the reacting fragments to the two zinc ions present in the active site, which has been observed to be dependent on whether or not metal atoms and their coordination spheres are included in the QM region. The description obtained using the largest QM region agrees with the picture obtained from experimental data.

O-glycoprotein 2-acetamino-2-deoxy-beta- d-glucopyranosidase ( O-GlcNAcase) hydrolyzes 2-acetamido-2-deoxy-beta- 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 simulations using a hybrid quantum mechanics/molecular mechanics 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 structure.

Hybrid quantum mechanics/molecular mechanics (QM/MM) techniques are widely used to study chemical reactions in large systems. Because of the computational cost associated with the high dimensionality of these systems, the quantum description is usually restricted to low-level methods, such as semiempirical Hamiltonians. In some cases, the description obtained at this computational level is quite poor and corrections must be considered. We here propose a simple but efficient way to include higher-level corrections to be used in potential energy surface explorations and in the calculation of potentials of mean force. We evaluate a correction energy term as the difference between a high-level and a low-level calculation on the QM subsystem, employing either the polarized or the gas-phase wave function, obtained as a function of two geometrical coordinates relevant in the process considered. Through the use of two-dimensional bicubic splines this correction energy is included in the simulations, ensuring the continuity and derivability of the energy function. We have tested the proposed scheme with two prototypical examples: the chorismate to prephenate rearrangement catalyzed by Bacillus subtilis chorismate mutase and the catechol methylation catalyzed by catechol O-methyltransferase. In both cases the use of interpolated corrections clearly improves the energetic and geometric descriptions of the reaction.