We have carried out a large scale computational investigation to assess the utility of common small‐molecule force fields for computational screening of low energy conformers of typical organic molecules. Using statistical analyses on the energies and relative rankings of up to 250 diverse conformers of 700 different molecular structures, we find that energies from widely used classical force fields (MMFF94, UFF, and GAFF) show unconditionally poor energy and rank correlation with semiempirical (PM7) and Kohn–Sham density functional theory (DFT) energies calculated at PM7 and DFT optimized geometries. In contrast, semiempirical PM7 calculations show significantly better correlation with DFT calculations and generally better geometries. With these results, we make recommendations to more reliably carry out conformer screening.
- PAR ID:
- 10448893
- Date Published:
- Journal Name:
- Molecules
- Volume:
- 28
- Issue:
- 3
- ISSN:
- 1420-3049
- Page Range / eLocation ID:
- 1277
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
null (Ed.)The basicity constant, or p K b , is an intrinsic physical property of bases that gives a measure of its proton affinity in macroscopic environments. While the p K b is typically defined in reference to the bulk aqueous phase, several studies have suggested that this value can differ significantly at the air–water interface (which can have significant ramifications for particle surface chemistry and aerosol growth modeling). To provide mechanistic insight into surface proton affinity, we carried out ab initio metadynamics calculations to (1) explore the free-energy profile of dimethylamine and (2) provide reasonable estimates of the p K b value in different solvent environments. We find that the free-energy profiles obtained with our metadynamics calculations show a dramatic variation, with interfacial aqueous dimethylamine p K b values being significantly lower than in the bulk aqueous environment. Furthermore, our metadynamics calculations indicate that these variations are due to reduced hydrogen bonding at the air–water surface. Taken together, our quantum mechanical metadynamics calculations show that the reactivity of dimethylamine is surprisingly complex, leading to p K b variations that critically depend on the different atomic interactions occurring at the microscopic molecular level.more » « less
-
more » « less
Abstract Computational modeling of chemical and biological systems at atomic resolution is a crucial tool in the chemist’s toolset. The use of computer simulations requires a balance between cost and accuracy: quantum-mechanical methods provide high accuracy but are computationally expensive and scale poorly to large systems, while classical force fields are cheap and scalable, but lack transferability to new systems. Machine learning can be used to achieve the best of both approaches. Here we train a general-purpose neural network potential (ANI-1ccx) that approaches CCSD(T)/CBS accuracy on benchmarks for reaction thermochemistry, isomerization, and drug-like molecular torsions. This is achieved by training a network to DFT data then using transfer learning techniques to retrain on a dataset of gold standard QM calculations (CCSD(T)/CBS) that optimally spans chemical space. The resulting potential is broadly applicable to materials science, biology, and chemistry, and billions of times faster than CCSD(T)/CBS calculations.
-
Performing alchemical transformations, in which one molecular system is nonphysically changed to another system, is a popular approach adopted in performing free energy calculations associated with various biophysical processes, such as protein–ligand binding or the transfer of a molecule between environments. While the sampling of alchemical intermediate states in either parallel (e.g., Hamiltonian replica exchange) or serial manner (e.g., expanded ensemble) can bridge the high-probability regions in the configurational space between two end states of interest, alchemical methods can fail in scenarios where the most important slow degrees of freedom in the configurational space are, in large part, orthogonal to the alchemical variable, or if the system gets trapped in a deep basin extending in both the configurational and alchemical space. To alleviate these issues, we propose to use alchemical variables as an additional dimension in metadynamics, making it possible to both sample collective variables and to enhance sampling in free energy calculations simultaneously. In this study, we validate our implementation of “alchemical metadynamics” in PLUMED with test systems and alchemical processes with varying complexities and dimensionalities of collective variable space, including the interconversion between the torsional metastable states of a toy system and the methylation of a nucleoside both in the isolated form and in a duplex. We show that multidimensional alchemical metadynamics can address the challenges mentioned above and further accelerate sampling by introducing configurational collective variables. The method can trivially be combined with other metadynamics-based algorithms implemented in PLUMED. The necessary PLUMED code changes have already been released for general use in PLUMED 2.8.more » « less
-
more » « less
The capabilities of the polarizable force fields for alchemical free energy calculations have been limited by the high computational cost and complexity of the underlying potential energy functions. In this work, we present a GPU‐based general alchemical free energy simulation platform for polarizable potential AMOEBA. Tinker‐OpenMM, the OpenMM implementation of the AMOEBA simulation engine has been modified to enable both absolute and relative alchemical simulations on GPUs, which leads to a ∼200‐fold improvement in simulation speed over a single CPU core. We show that free energy values calculated using this platform agree with the results of Tinker simulations for the hydration of organic compounds and binding of host–guest systems within the statistical errors. In addition to absolute binding, we designed a relative alchemical approach for computing relative binding affinities of ligands to the same host, where a special path was applied to avoid numerical instability due to polarization between the different ligands that bind to the same site. This scheme is general and does not require ligands to have similar scaffolds. We show that relative hydration and binding free energy calculated using this approach match those computed from the absolute free energy approach. © 2017 Wiley Periodicals, Inc.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/molecules28031277
