A systematic density functional theory exploration of various reactive steps together with benchmark coupled cluster results are used to propose an accurate model of the mechanism of the Morita–Baylis–Hillman (MBH) reaction in organic chemistry. This reaction has attracted considerable interest from the synthetic and mechanistic points of view in recent years, with both computational and experimental mechanistic studies. It has recently (R. E. Plata and D. A. Singleton, J. Am. Chem. Soc., 2015, 137, 3811–3826) been correctly pointed out that previous computational studies failed to reproduce known mechanistic features of the reaction. The same study argued that computation is simply not able at the present time to provide accurate models for such reactions. This second claim is shown by our present work to overstate the problem: by using current ‘state of the art’ methodology, our results are fully consistent with observed behavior within the expected error bars of 1–5 kcal mol−1, far smaller than the errors reported in Plata and Singleton's study. On the basis of exhaustive calculations reported here, we suggest that our proposed approaches for modeling electronic structure, solvation, and entropy should be able to provide accurate predictions for many more reactions. We also suggest that reactions like the MBH reaction, where solvation and entropy effects are particularly large, are among the hardest for computational mechanistic studies.
21 August 2009 in The Journal of Organic Chemistry
Plasmodium falciparumplasmepsin V (PfPMV) is an essential aspartic protease required for parasite survival, thus, considered as a potential drug target. This study reports the first detailed structural analysis and molecular dynamics simulation of PfPMV as an apoenzyme and its complexes with the substrate PEXEL as well as with the inhibitor saquinavir. The presence of pro-peptide in PfPMV may not structurally hinder the formation of a functionally competent catalytic active site. The structure of PfPMV-PEXEL complex shows that the unique positions of Glu179 and Gln222 are responsible for providing the specificity of PEXEL substrate with arginine at P3 position. The structural analysis also reveals that the S4 binding pocket in PfPMV is occupied by Ile94, Ala98, Phe370 and Tyr472 and therefore, does not allow binding of pepstatin, a potent inhibitor of most pepsin-like aspartic proteases. Among the screened inhibitors, the HIV-1 protease inhibitors and KNI compounds have higher binding affinities for PfPMV with saquinavir having the highest value. The presence of a flexible group at P2 and a bulky hydrophobic group at P3 position of the inhibitor is preferred in the PfPMV substrate binding pocket. Results from the present study will aid in the design of potent inhibitors of PMV.