MF-LOGP, a new method for determining a single component octanol–water partition coefficients (
The logarithm of
- Award ID(s):
- 1955260
- NSF-PAR ID:
- 10398582
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Journal of Computational Chemistry
- Volume:
- 44
- Issue:
- 13
- ISSN:
- 0192-8651
- Page Range / eLocation ID:
- p. 1300-1311
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract $$LogP$$ $$LogP$$ $$RMSE$$ $$MAE$$ $${R}^{2}$$ $$RMSE$$ $$MAE$$ $${R}^{2}$$ Graphical Abstract -
more » « less
Abstract Accurate estimation of solvation free energy (SFE) lays the foundation for accurate prediction of binding free energy. The Poisson‐Boltzmann (PB) or generalized Born (GB) combined with surface area (SA) continuum solvation method (PBSA and GBSA) have been widely used in SFE calculations because they can achieve good balance between accuracy and efficiency. However, the accuracy of these methods can be affected by several factors such as the charge models, polar and nonpolar SFE calculation methods and the atom radii used in the calculation. In this work, the performance of the ABCG2 (AM1‐BCC‐GAFF2) charge model as well as other two charge models, that is, RESP (Restrained Electrostatic Potential) and AM1‐BCC (Austin Model 1‐bond charge corrections), on the SFE prediction of 544 small molecules in water by PBSA/GBSA was evaluated. In order to improve the performance of the PBSA prediction based on the ABCG2 charge, we further explored the influence of atom radii on the prediction accuracy and yielded a set of atom radius parameters for more accurate SFE prediction using PBSA based on the ABCG2/GAFF2 by reproducing the thermodynamic integration (TI) calculation results. The PB radius parameters of carbon, oxygen, sulfur, phosphorus, chloride, bromide and iodine, were adjusted. New atom types,
on oi hn1 hn2 hn3 -
more » « less
Abstract This work examines methods for predicting the partition coefficient (log
P ) for a dataset of small molecules. Here, we use atomic attributes such as radius and partial charge, which are typically used as force field parameters in classical molecular dynamics simulations. These atomic attributes are transformed into index‐invariant molecular features using a recently developed method called geometric scattering for graphs (GSG). We call this approach “ClassicalGSG” and examine its performance under a broad range of conditions and hyperparameters. We train ClassicalGSG logP predictors with neural networks using 10,722 molecules from the OpenChem dataset and apply them to predict the logP values from four independent test sets. The ClassicalGSG method's performance is compared to a baseline model that employs graph convolutional networks. Our results show that the best prediction accuracies are obtained using atomic attributes generated with the CHARMM generalized force field and 2D molecular structures. -
more » « less
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
Abstract Soil thermal conductivity (λ) is an important thermal property for environmental, agricultural, and engineering heat transfer applications. Existing λ models for frozen soils are complicated to use because they require estimates of both liquid water content and ice content. This study introduces a new approach to estimate λ of partially frozen soils from air‐filled porosity (
n a), which can be determined by using an oven‐drying method. A λ andn arelationship was established based on measurements for 28 partially frozen soils. A strong exponential relationship between λ andn awas found (withR 2of 0.82). Independent tests on 10 partially frozen soils showed that the exponential λ–n amodel produced reliable λ estimates with a RMSE of 0.319 W m−1K−1, which was smaller than those of two widely used λ models for partially frozen soils. The λ–n amodel is easier to use than existing models, because it requires fewer parameters. Note that the λ‐n amodel ignores the effect of temperature on λ of frozen soils and is most applicable to soil at temperatures of at least −4 °C.
