An official website of the United States government
Two new flexible-monomer two-body ab initio potential energy surfaces (PESs) for the neon and krypton van der Waals complexes with carbon dioxide were developed, extending our previous work on the Ar–CO2 molecule. The accuracy of the PESs was validated by their agreement with the vibrational spectrum of the rare-gas complexes. The intermolecular and intramolecular vibrational excitation energies were computed at the vibrational self-consistent field and vibrational configuration interaction levels of theory. Overall, the agreement between theory and experiment is excellent throughout the vibrational spectra. The observed slight splitting of the bending modes, resulting from their nondegeneracy in the complexes, is confirmed by our computations, and the results qualitatively agree with the experiment. The splitting increases with increasing polarizability of the rare-gas atom. Additionally, we explain a discrepancy in the mode assignment in the intermolecular region of the neon complex with our VCI character assignment.
more »
« less
Automation of AMOEBA polarizable force field for small molecules: Poltype 2
Abstract 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):
- 1856173
- PAR ID:
- 10375629
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Journal of Computational Chemistry
- Volume:
- 43
- Issue:
- 23
- ISSN:
- 0192-8651
- Page Range / eLocation ID:
- p. 1530-1542
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
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
-
Hydration free energies of small molecules are commonly used as benchmarks for solvation models. However, errors in predicting hydration free energies are partially due to the force fields used and not just the solvation model. To address this, we have used the 3D reference interaction site model (3D-RISM) of molecular solvation and existing benchmark explicit solvent calculations with a simple element count correction (ECC) to identify problems with the non-bond parameters in the general AMBER force field (GAFF). 3D-RISM was used to calculate hydration free energies of all 642 molecules in the FreeSolv database, and a partial molar volume correction (PMVC), ECC, and their combination (PMVECC) were applied to the results. The PMVECC produced a mean unsigned error of 1.01±0.04kcal/mol and root mean squared error of 1.44±0.07kcal/mol, better than the benchmark explicit solvent calculations from FreeSolv, and required less than 15 s of computing time per molecule on a single CPU core. Importantly, parameters for PMVECC showed systematic errors for molecules containing Cl, Br, I, and P. Applying ECC to the explicit solvent hydration free energies found the same systematic errors. The results strongly suggest that some small adjustments to the Lennard–Jones parameters for GAFF will lead to improved hydration free energy calculations for all solvent models.more » « less
-
We measure and model monolayers of concentrated diffusing colloidal probes interacting with polymerized liquid crystal (PLC) planar surfaces. At topological defects in local nematic director profiles at PLC surfaces, we observe time-averaged two-dimensional particle density profiles of diffusing colloidal probes that closely correlate with spatial variations in PLC optical properties. An inverse Monte Carlo analysis of particle concentration profiles yields two-dimensional PLC interfacial energy landscapes on the kT -scale, which is the inherent scale of many interfacial phenomena ( e.g. , self-assembly, adsorption, diffusion). Energy landscapes are modelled as the superposition of macromolecular repulsion and van der Waals attraction based on an anisotropic dielectric function obtained from the liquid crystal birefringence. Modelled van der Waals landscapes capture most net energy landscape variations and correlate well with experimental PLC director profiles around defects. Some energy landscape variations near PLC defects indicate either additional local repulsive interactions or possibly the need for more rigorous van der Waals models with complete spectral data. These findings demonstrate direct, sensitive measurements of kT -scale van der Waals energy landscapes at PLC interfacial defects and suggest the ability to design interfacial anisotropic materials and van der Waals energy landscapes for colloidal assembly.more » « less
