Achieving High-Accuracy Intermolecular Interactions by Combining Coulomb-Attenuated Second-Order Møller–Plesset Perturbation Theory with Coupled Kohn–Sham Dispersion

Improved Description of Intra- and Intermolecular Interactions through Dispersion-Corrected Second-Order Møller-Plesset Perturbation Theory

ConspectusThe quantum chemical modeling of organic crystals and other molecular condensed-phase problems requires computationally affordable electronic structure methods which can simultaneously describe intramolecular conformational energies and intermolecular interactions accurately. To achieve this, we have developed a spin-component-scaled, dispersion-corrected second-order Møller-Plesset perturbation theory (SCS-MP2D) model. SCS-MP2D augments canonical MP2 with a dispersion correction which removes the uncoupled Hartree-Fock dispersion energy present in canonical MP2 and replaces it with a more reliable coupled Kohn-Sham treatment, all evaluated within the framework of Grimme's D3 dispersion model. The spin-component scaling is then used to improve the description of the residual (nondispersion) portion of the correlation energy.The SCS-MP2D model improves upon earlier corrected MP2 models in a few ways. Compared to the highly successful dispersion-corrected MP2C model, which is based solely on intermolecular perturbation theory, the SCS-MP2D dispersion correction improves the description of both inter- and intramolecular interactions. The dispersion correction can also be evaluated with trivial computational cost, and nuclear analytic gradients are computed readily to enable geometry optimizations. In contrast to earlier spin-component scaling MP2 models, the optimal spin-component scaling coefficients are only mildly sensitive to the choice of training data, and a single global parametrization of the model can describe both thermochemistry and noncovalent interactions.The resulting dispersion-corrected, spin-component-scaled MP2 (SCS-MP2D) model predicts conformational energies and intermolecular interactions with accuracy comparable to or better than those of many range-separated and double-hybrid density functionals, as is demonstrated on a variety of benchmark tests. Among the functionals considered here, only the revDSD-PBEP86-D3(BJ) functional gives consistently smaller errors in benchmark tests. The results presented also hint that further improvements of SCS-MP2D may be possible through a more robust fitting procedure for the seven empirical parameters.To demonstrate the performance of SCS-MP2D further, several applications to molecular crystal problems are presented. The three chosen examples all represent cases where density-driven delocalization error causes GGA or hybrid density functionals to artificially stabilize crystals exhibiting more extended π-conjugation. Our pragmatic strategy addresses the delocalization error by combining a periodic density functional theory (DFT) treatment of the infinite lattice with intramolecular/conformational energy corrections computed with SCS-MP2D. For the anticancer drug axitinib, applying the SCS-MP2D conformational energy correction produces crystal polymorph stabilities that are consistent with experiment, in contrast to earlier studies. For the crystal structure prediction of the ROY molecule, so named for its colorful red, orange, and yellow crystals, this approach leads to the first plausible crystal energy landscape, and it reveals that the lowest-energy polymorphs have already been found experimentally. Finally, in the context of photomechanical crystals, which transform light into mechanical work, these techniques are used to predict the structural transformations and extract design principles for maximizing the work performed.

Computational Study of the Addition of Methanethiol to 40+ Michael Acceptors as a Model for the Bioconjugation of Cysteines

A long series of Michael acceptors are studied computationally as potential alternatives to the maleimides that are used in most antibody-drug conjugates to link Cys of mAbs with cytotoxic drugs. The products of the reaction of methanethiol (CH3SH/MeSH, as a simple model of Cys) with N-methylated ethynesulfonamide, 2-ethynylpyridinium ion, propynamide, and methyl ethynephosphonamidate (that is, with HC≡C-EWG) are predicted by the M06-2X/6-311+G(d,p) method to be thermodynamically more stable, in relation to their precursors, than that of MeSH with N-methylmaleimide and, in general, with H2C═CH-EWG; calculations with AcCysOMe and tBuSH are also included. However, for the addition of the anion (MeS-), which is the reactive species, the order changes and N-methylated 2-vinylpyridinium ion, 2,3-butadienamide, and maleimide may give more easily the anionic adducts than several activated triple bonds; moreover, the calculated ΔG⧧ values increase following the order HC≡C-SO2NHMe, N-methylmaleimide, HC≡C-PO(OMe)NHMe, and HC≡C-CONHMe. In other words, MeS- is predicted to react more rapidly with maleimides than with ethynephosphonamidates and with propynamides, in agreement with the experimental results. New mechanistic details are disclosed regarding the advantageous use of some amides, especially of ethynesulfonamides, which, however, are more prone to double additions and exchange reactions.

The Good, the Bad, and the Ugly: "HiPen", a New Dataset for Validating (S)QM/MM Free Energy Simulations

Indirect (S)QM/MM free energy simulations (FES) are vital to efficiently incorporating sufficient sampling and accurate (QM) energetic evaluations when estimating free energies of practical/experimental interest. Connecting between levels of theory, i.e., calculating Δ A l o w → h i g h , remains to be the most challenging step within an indirect FES protocol. To improve calculations of Δ A l o w → h i g h , we must: (1) compare the performance of all FES methods currently available; and (2) compile and maintain datasets of Δ A l o w → h i g h calculated for a wide-variety of molecules so that future practitioners may replicate or improve upon the current state-of-the-art. Towards these two aims, we introduce a new dataset, "HiPen", which tabulates Δ A g a s M M → 3 o b (the free energy associated with switching from an M M to an S C C - D F T B molecular description using the 3ob parameter set in gas phase), calculated for 22 drug-like small molecules. We compare the calculation of this value using free energy perturbation, Bennett's acceptance ratio, Jarzynski's equation, and Crooks' equation. We also predict the reliability of each calculated Δ A g a s M M → 3 o b by evaluating several convergence criteria including sample size hysteresis, overlap statistics, and bias metric ( Π ). Within the total dataset, three distinct categories of molecules emerge: the "good" molecules, for which we can obtain converged Δ A g a s M M → 3 o b using Jarzynski's equation; "bad" molecules which require Crooks' equation to obtain a converged Δ A g a s M M → 3 o b ; and "ugly" molecules for which we cannot obtain reliably converged Δ A g a s M M → 3 o b with either Jarzynski's or Crooks' equations. We discuss, in depth, results from several example molecules in each of these categories and describe how dihedral discrepancies between levels of theory cause convergence failures even for these gas phase free energy simulations.

Modeling Polymorphic Molecular Crystals with Electronic Structure Theory

Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.

Importance of the Electron Correlation and Dispersion Corrections in Calculations Involving Enamines, Hemiaminals, and Aminals. Comparison of B3LYP, M06-2X, MP2, and CCSD Results with Experimental Data

While B3LYP, M06-2X, and MP2 calculations predict the ΔG° values for exchange equilibria between enamines and ketones with similar acceptable accuracy, the M06-2X/6-311+G(d,p) and MP2/6-311+G(d,p) methods are required for enamine formation reactions (for example, for enamine 5a, arising from 3-methylbutanal and pyrrolidine). Stronger disagreement was observed when calculated energies of hemiaminals (N,O-acetals) and aminals (N,N-acetals) were compared with experimental equilibrium constants, which are reported here for the first time. Although it is known that the B3LYP method does not provide a good description of the London dispersion forces, while M06-2X and MP2 may overestimate them, it is shown here how large the gaps are and that at least single-point calculations at the CCSD(T)/6-31+G(d) level should be used for these reaction intermediates; CCSD(T)/6-31+G(d) and CCSD(T)/6-311+G(d,p) calculations afford ΔG° values in some cases quite close to MP2/6-311+G(d,p) while in others closer to M06-2X/6-311+G(d,p). The effect of solvents is similarly predicted by the SMD, CPCM, and IEFPCM approaches (with energy differences below 1 kcal/mol).