Polarizable AMOEBA Model for Simulating Mg2+·Protein·Nucleotide Complexes.

Molecular mechanics (MM) simulations have the potential to provide detailed insights into the mechanisms of enzymes that utilize nucleotides as cofactors. In most cases, the activities of these enzymes also require the binding of divalent cations to catalytic sites. However, modeling divalent cations in MM simulations has been challenging. The inclusion of explicit polarization was considered promising, but despite improvements over nonpolarizable force fields and despite the inclusion of "Nonbonded-fix (NB-fix)" corrections, errors in interaction energies of divalent cations with proteins remain large. Importantly, the application of these models fails to reproduce the experimental structural data on Mg2+·Protein·ATP complexes. Focusing on these complexes, here we provide a systematic assessment of the polarizable AMOEBA model and recommend critical changes that substantially improve its predictive performance. Our key results are as follows. We first show that our recent revision of the AMOEBA protein model (AMOEBABIO18-HFC), which contains high field corrections (HFCs) to induced dipoles, dramatically improves Mg2+-protein interaction energies, reducing the mean absolute error (MAE) from 17 to 10 kcal/mol. This further supports the general applicability of AMOEBABIO18-HFC. The inclusion of many-body NB-fix corrections further reduces MAE to 6 kcal/mol, which amounts to less than 2% error. The errors are estimated with respect to vdW-inclusive density functional theory that we benchmark against CCSD(T) calculations and experiments. We also present a new model of ATP with revised polarization parameters to better capture its high field response, as well as new vdW and dihedral parameters. The ATP model accurately predicts experimental Mg2+-ATP binding free energy in the aqueous phase and provides new insights into how Mg2+ associates with ATP. Finally, we show that molecular dynamics (MD) simulations of Mg2+·Kinase·ATP complexes carried out with these improvements lead to a better agreement in global and local catalytic site structures between MD and X-ray crystallography.

Stereocontrol in Diphenylprolinol Silyl Ether Catalyzed Michael Additions: Steric Shielding or Curtin-Hammett Scenario?

The enantioselectivity of amine-catalyzed reactions of aldehydes with electrophiles is often explained by simple steric arguments emphasizing the role of the bulky group of the catalyst that prevents the approach of the electrophile from the more hindered side. This standard steric shielding model has recently been challenged by the discovery of stable downstream intermediates, which appear to be involved in the rate-determining step of the catalytic cycle. The alternative model, referred to as the Curtin-Hammett scenario of stereocontrol, assumes that the enantioselectivity is related to the stability and reactivity of downstream intermediates. In our present computational study, we examine the two key processes of the catalytic Michael reaction between propanal and ß-nitrostyrene that are relevant to the proposed stereoselectivity models, namely the C-C bond formation and the protonation steps. The free energy profiles obtained for the pathways leading to the enantiomeric products suggest that the rate- and stereodetermining steps are not identical as implied by the previous models. The stereoselectivity can be primarily controlled by C-C bond formation even though the reaction rate is dictated by the protonation step. This kinetic scheme is consistent with all observations of experimental mechanistic studies including those of mass spectrometric back reaction screening experiments, which reveal a mismatch between the stereoselectivity of the back and the forward reactions.