skip to main content

Title: Automation of AMOEBA polarizable force field for small molecules: Poltype 2

A next‐generation protocol (Poltype 2) has been developed which automatically generates AMOEBA polarizable force field parameters for small molecules. Both features and computational efficiency have been drastically improved. Notable advances include improved database transferability using SMILES, robust torsion fitting, non‐aromatic ring torsion parameterization, coupled torsion‐torsion parameterization, Van der Waals parameter refinement using ab initio dimer data and an intelligent fragmentation scheme that produces parameters with dramatically reduced ab initio computational cost. Additional improvements include better local frame assignment for atomic multipoles, automated formal charge assignment, Zwitterion detection, smart memory resource defaults, parallelized fragment job submission, incorporation of Psi4 quantum package, ab initio error handling, ionization state enumeration, hydration free energy prediction and binding free energy prediction. For validation, we have applied Poltype 2 to ~1000 FDA approved drug molecules from DrugBank. The ab initio molecular dipole moments and electrostatic potential values were compared with Poltype 2 derived AMOEBA counterparts. Parameters were further substantiated by calculating hydration free energy (HFE) on 40 small organic molecules and were compared with experimental data, resulting in an RMSE error of 0.59 kcal/mol. The torsion database has expanded to include 3543 fragments derived from FDA approved drugs. Poltype 2 provides a convenient utility for applications including binding free energy prediction for computational drug discovery. Further improvement will focus on automated parameter refinement by experimental liquid properties, expansion of the Van der Waals parameter database and automated parametrization of modified bio‐fragments such as amino and nucleic acids.

more » « less
Award ID(s):
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Journal of Computational Chemistry
Page Range / eLocation ID:
p. 1530-1542
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We describe a strategy of integrating quantum mechanical (QM), hybrid quantum mechanical/molecular mechanical (QM/MM) and MM simulations to analyze the physical properties of a solid/water interface. This protocol involves using a correlated ab initio (CCSD(T)) method to first calibrate Density Functional Theory (DFT) as the QM approach, which is then used in QM/MM simulations to compute relevant free energy quantities at the solid/water interface using a mean-field approximation of Yang et al. that decouples QM and MM thermal fluctuations; gas-phase QM/MM and periodic DFT calculations are used to determine the proper QM size in the QM/MM simulations. Finally, the QM/MM free energy results are compared with those obtained from MM simulations to directly calibrate the force field model for the solid/water interface. This protocol is illustrated by examining the orientations of an alkyl amine ligand at the gold/water interface, since the ligand conformation is expected to impact the chemical properties ( e.g. , charge) of the solid surface. DFT/MM and MM simulations using the INTERFACE force field lead to consistent results, suggesting that the effective gold/ligand interactions can be adequately described by a van der Waals model, while electrostatic and induction effects are largely quenched by solvation. The observed differences among periodic DFT, QM/MM and MM simulations, nevertheless, suggest that explicitly including electronic polarization and potentially charge transfer in the MM model can be important to the quantitative accuracy. The strategy of integrating multiple computational methods to cross-validate each other for complex interfaces is applicable to many problems that involve both inorganic/metallic and organic/biomolecular components, such as functionalized nanoparticles. 
    more » « less
  2. Computational simulation of biomolecules can provide important insights into protein design, protein-ligand binding interactions, and ab initio biomolecular folding, among other applications. Accurate treatment of the solvent environment is essential in such applications, but the use of explicit solvents can add considerable cost. Implicit treatment of solvent effects using a dielectric continuum model is an attractive alternative to explicit solvation since it is able to describe solvation effects without the inclusion of solvent degrees of freedom. Previously, we described the development and parameterization of implicit solvent models for small molecules. Here, we extend the parameterization of the generalized Kirkwood (GK) implicit solvent model for use with biomolecules described by the AMOEBA force field via the addition of corrections to the calculation of effective radii that account for interstitial spaces that arise within biomolecules. These include element-specific pairwise descreening scale factors, a short-range neck contribution to describe the solvent-excluded space between pairs of nearby atoms, and finally tanh-based rescaling of the overall descreening integral. We then apply the AMOEBA/GK implicit solvent to a set of ten proteins and achieve an average coordinate root mean square deviation for the experimental structures of 2.0 Å across 500 ns simulations. Overall, the continued development of implicit solvent models will help facilitate the simulation of biomolecules on mechanistically relevant timescales. 
    more » « less
  3. The structure of zymonic acid (systematic name: 4-hydroxy-2-methyl-5-oxo-2,5-dihydrofuran-2-carboxylic acid), C 6 H 6 O 5 , which had previously eluded crystallographic determination, is presented here for the first time. It forms by intramolecular condensation of parapyruvic acid, which is the product of aldol condensation of pyruvic acid. A redetermination of the crystal structure of pyruvic acid (systematic name: 2-oxopropanoic acid), C 3 H 4 O 3 , at low temperature (90 K) and with increased precision, is also presented [for the previous structure, see: Harata et al. (1977). Acta Cryst. B 33 , 210–212]. In zymonic acid, the hydroxylactone ring is close to planar (r.m.s. deviation = 0.0108 Å) and the dihedral angle between the ring and the plane formed by the bonds of the methyl and carboxylic acid carbon atoms to the ring is 88.68 (7)°. The torsion angle of the carboxylic acid group relative to the ring is 12.04 (16)°. The pyruvic acid molecule is almost planar, having a dihedral angle between the carboxylic acid and methyl-ketone groups of 3.95 (6)°. Intermolecular interactions in both crystal structures are dominated by hydrogen bonding. The common R 2 2 (8) hydrogen-bonding motif links carboxylic acid groups on adjacent molecules in both structures. In zymonic acid, this results in dimers about a crystallographic twofold of space group C 2/ c , which forces the carboxylic acid group to be disordered exactly 50:50, which scrambles the carbonyl and hydroxyl groups and gives an apparent equalization of the C—O bond lengths [1.2568 (16) and 1.2602 (16) Å]. The other hydrogen bonds in zymonic acid (O—H...O and weak C—H...O), link molecules across a 2 1 -screw axis, and generate an R 2 2 (9) motif. These hydrogen-bonding interactions propagate to form extended pleated sheets in the ab plane. Stacking of these zigzag sheets along c involves only van der Waals contacts. In pyruvic acid, inversion-related molecules are linked into R 2 2 (8) dimers, with van der Waals interactions between dimers as the only other intermolecular contacts. 
    more » « less
  4. null (Ed.)
    The ability to locally tune solute–water interactions and thus control the hydrophilic/hydrophobic character of a solute is key to control molecular self-assembly and to develop new drugs and biocatalysts; it has been a holy grail in synthetic chemistry and biology. To date, the connection between (i) the hydrophobicity of a functional group; (ii) the local structure and thermodynamics of its hydration shell; and (iii) the relative influence of van der Waals (dispersion) and electrostatic interactions on hydration remains unclear. We investigate this connection using spectroscopic, classical simulation and ab initio methods by following the transition from hydrophile to hydrophobe induced by the step-wise fluorination of methyl groups. Along the transition, we find that water–solute hydrogen bonds are progressively transformed into dangling hydroxy groups. Each structure has a distinct thermodynamic, spectroscopic and quantum-mechanical signature connected to the associated local solute hydrophobicity and correlating with the relative contribution of electrostatics and dispersion to the solute–water interactions. 
    more » « less
  5. Recent technology development of logic devices based on 2-D semiconductors such as MoS2, WS2, and WSe2 has triggered great excitement, paving the way to practical applications. Making low-resistance p-type contacts to 2-D semiconductors remains a critical challenge. The key to addressing this challenge is to find high-work function metallic materials which also introduce minimal metal-induced gap states (MIGSs) at the metal/semiconductor interface. In this work, we perform a systematic computational screening of novel metallic materials and their heterojunctions with monolayer WSe2 based on ab initio density functional theory and quantum device simulations. Two contact strategies, van der Waals (vdW) metallic contact and bulk semimetallic contact, are identified as promising solutions to achieving Schottky-barrier-free and low-contact-resistance p-type contacts for WSe2 p-type field-effect transistor (pFETs). Good candidates of p-type contact materials are found based on our screening criteria, including 1H-NbS2, 1H-TaS2, and 1T-TiS2 in the vdW metal category, as well as Co3Sn2S2 and TaP in the bulk semimetal category. Simulations of these new p-type contact materials suggest reduced MIGS, less Fermi-level pinning effect, negligible Schottky barrier height and small contact resistance (down to 20 Ωμm ) 
    more » « less