Fitting Molecular Electrostatic Potentials from Quantum Mechanical Calculations

Mechanism of OMP decarboxylation in orotidine 5'-monophosphate decarboxylase

Despite extensive experimental and theoretical studies, the detailed catalytic mechanism of orotidine 5'-monophosphate decarboxylase (ODCase) remains controversial. In particular simulation studies using high level quantum mechanics have failed to reproduce experimental activation free energy. One common feature of many previous simulations is that there is a water molecule in the vicinity of the leaving CO2 group whose presence was only observed in the inhibitor bound complex of ODCase/BMP. Various roles have even been proposed for this water molecule from the perspective of stabilizing the transition state and/or intermediate state. We hypothesize that this water molecule is not present in the active ODCase/OMP complex. Based on QM/MM minimum free energy path simulations with accurate density functional methods, we show here that in the absence of this water molecule the enzyme functions through a simple direct decarboxylation mechanism. Analysis of the interactions in the active site indicates multiple factors contributing to the catalysis, including the fine-tuned electrostatic environment of the active site and multiple hydrogen-bonding interactions. To understand better the interactions between the enzyme and the inhibitor BMP molecule, simulations were also carried out to determine the binding free energy of this special water molecule in the ODCase/BMP complex. The results indicate that the water molecule in the active site plays a significant role in the binding of BMP by contributing approximately -3 kcal/mol to the binding free energy of the complex. Therefore, the complex of BMP plus a water molecule, instead of the BMP molecule alone, better represents the tight binding transition state analogue of ODCase. Our simulation results support the direct decarboxylation mechanism and highlight the importance of proper recognition of protein bound water molecules in the protein-ligand binding and the enzyme catalysis.

Free energies of chemical reactions in solution and in enzymes with ab initio quantum mechanics/molecular mechanics methods

Combined quantum mechanics/molecular mechanics (QM/MM) methods provide an accurate and efficient energetic description of complex chemical and biological systems, leading to significant advances in the understanding of chemical reactions in solution and in enzymes. Here we review progress in QM/MM methodology and applications, focusing on ab initio QM-based approaches. Ab initio QM/MM methods capitalize on the accuracy and reliability of the associated quantum-mechanical approaches, however, at a much higher computational cost compared with semiempirical quantum-mechanical approaches. Thus reaction-path and activation free-energy calculations based on ab initio QM/MM methods encounter unique challenges in simulation timescales and phase-space sampling. This review features recent developments overcoming these challenges and enabling accurate free-energy determination for reaction processes in solution and in enzymes, along with applications.

Numerical fitting of molecular properties to Hermite Gaussians

A procedure is presented to fit gridded molecular properties to auxiliary basis sets (ABSs) of Hermite Gaussians, analogous to the density fitting (DF) method (Dunlap; et al. J. Chem. Phys. 1979, 71, 4993). In this procedure, the ab initio calculated properties (density, electrostatic potential, and/or electric field) are fitted via a linear- or nonlinear-least-squares procedure to auxiliary basis sets (ABS). The calculated fitting coefficients from the numerical grids are shown to be more robust than analytic density fitting due to the neglect of the core contributions. The fitting coefficients are tested by calculating intermolecular Coulomb and exchange interactions for a set of dimers. It is shown that the numerical instabilities observed in DF are caused by the attempt of the ABS to fit the core contributions. In addition, this new approach allows us to reduce the number of functions required to obtain an accurate fit. This results in decreased computational cost, which is shown by calculating the Coulomb energy of a 4096 water box in periodic boundary conditions. Using atom centered Hermite Gaussians, this calculation is only 1 order of magnitude slower than conventional atom-centered point charges.