1
|
Draper MR, Waterman A, Dannatt JE, Patel P. Integrating multiscale and machine learning approaches towards the SAMPL9 log P challenge. Phys Chem Chem Phys 2024; 26:7907-7919. [PMID: 38376855 PMCID: PMC10938873 DOI: 10.1039/d3cp04140a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/21/2024]
Abstract
The partition coefficient (log P) is an important physicochemical property that provides information regarding a molecule's pharmacokinetics, toxicity, and bioavailability. Methods to accurately predict the partition coefficient have the potential to accelerate drug design. In an effort to test current methods and explore new computational techniques, the statistical assessment of the modeling of proteins and ligands (SAMPL) has established a blind prediction challenge. The ninth iteration challenge was to predict the toluene-water partition coefficient (log Ptol/w) of sixteen drug molecules. Herein, three approaches are reported broadly under the categories of quantum mechanics (QM), molecular mechanics (MM), and data-driven machine learning (ML). The three blind submissions yield mean unsigned errors (MUE) ranging from 1.53-2.93 log Ptol/w units. The MUEs were reduced to 1.00 log Ptol/w for the QM methods. While MM and ML methods outperformed DFT approaches for challenge molecules with fewer rotational degrees of freedom, they suffered for the larger molecules in this dataset. Overall, DFT functionals paired with a triple-ζ basis set were the simplest and most effective tool to obtain quantitatively accurate partition coefficients.
Collapse
Affiliation(s)
- Michael R Draper
- Chemistry Department, University of Dallas, Irving, Texas, 75062, USA.
| | - Asa Waterman
- Chemistry Department, University of Dallas, Irving, Texas, 75062, USA.
| | | | - Prajay Patel
- Chemistry Department, University of Dallas, Irving, Texas, 75062, USA.
| |
Collapse
|
2
|
Aceves-Luna H, Glossman-Mitnik D, Flores-Holguín N. Permeability of antioxidants through a lipid bilayer model with coarse-grained simulations. J Biomol Struct Dyn 2023:1-19. [PMID: 37768552 DOI: 10.1080/07391102.2023.2262044] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Accepted: 09/13/2023] [Indexed: 09/29/2023]
Abstract
Oxidative stress caused by pollution and lifestyle changes causes an excess of free radicals that react chemically with cell constituents leading to irreversible damage. There are molecules known as antioxidants that reduce the levels of free radicals. Some pigments of fruits and vegetables known as anthocyanins have antioxidant properties. Their interaction with the cell membrane becomes a crucial step in studying these substances. In this research, molecular dynamics simulations, particularly, coarse-grained molecular dynamics (CGMD) were used. This technique aims to replace functional groups with corresponding beads that represent their level of polarity and affinities to other chemical groups. Also, umbrella sampling was carried out to obtain free energy profiles that describe well the orientation and location of antioxidants in a membrane considering Trolox, Cyanidin, Gallic Acid, and Resveratrol molecules to study the structural effects they cause on it. It was concluded in this study that an antioxidant when crossing the membrane does not cause either damage to the structural properties or the loss of packing and stratification of phospholipids. it was also observed that the most reactive part of the molecules could easily approach area A prone to lipid oxidation, which can describe the antioxidant capacity of these molecules.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Hugo Aceves-Luna
- Laboratorio Virtual NANOCOSMOS, Departamento de Medio Ambiente y Energía, Centro de Investigación en Materiales Avanzados, Chihuahua, Chih, Mexico
| | - Daniel Glossman-Mitnik
- Laboratorio Virtual NANOCOSMOS, Departamento de Medio Ambiente y Energía, Centro de Investigación en Materiales Avanzados, Chihuahua, Chih, Mexico
| | - Norma Flores-Holguín
- Laboratorio Virtual NANOCOSMOS, Departamento de Medio Ambiente y Energía, Centro de Investigación en Materiales Avanzados, Chihuahua, Chih, Mexico
| |
Collapse
|
3
|
Jia Q, Ni Y, Liu Z, Gu X, Cui Z, Fan M, Zhu Q, Wang Y, Ma J. Fast Prediction of Lipophilicity of Organofluorine Molecules: Deep Learning-Derived Polarity Characters and Experimental Tests. J Chem Inf Model 2022; 62:4928-4936. [PMID: 36223527 DOI: 10.1021/acs.jcim.2c01201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Fast and accurate estimation of lipophilicity for organofluorine molecules is in great demand for accelerating drug and materials discovery. A lipophilicity data set of organofluorine molecules (OFL data set), containing 1907 samples, is constructed through density functional theory (DFT) calculations and experimental measurements. An efficient and interpretable model, called PoLogP, is developed to predict the n-octanol/water partition coefficient, log Po/w, of organofluorine molecules on the basis of the descriptors of polarization, which is a combination of polarity descriptors, including the molecular polarity index and molecular polarizability (α), and hydrogen bond (HBs) index, consisting of the number of donors (NHBD) and acceptors (NHBA and NHB-FA). The present PoLogP with a combination of polarity descriptors is demonstrated to perform better than the dipole moment (μ) alone for the F-contained molecules. With the aid of a multilevel attention graph convolutional neural network model, the fast generation of polarity descriptors of organofluorine molecules could be achieved with the DFT accuracy based only on a topological molecular graph structure. The performance of PoLogP is further validated on synthesized organofluorine molecules and 2626 non-fluorinated molecules with satisfactory accuracy, highlighting the potential usage of PoLogP in high-throughput screening of the functional molecules with the desired solubility in various solvent media.
Collapse
Affiliation(s)
- Qingqing Jia
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China
| | - Yifan Ni
- Jiangsu Key Laboratory of Advanced Organic Materials, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China
| | - Ziteng Liu
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China
| | - Xu Gu
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China
| | - Ziyi Cui
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China
| | - Mengting Fan
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China
| | - Qiang Zhu
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China
| | - Yi Wang
- Jiangsu Key Laboratory of Advanced Organic Materials, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China
| | - Jing Ma
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China.,Jiangsu Key Laboratory of Advanced Organic Materials, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China
| |
Collapse
|
4
|
Zhu Q, Jia Q, Liu Z, Ge Y, Gu X, Cui Z, Fan M, Ma J. Molecular partition coefficient from machine learning with polarization and entropy embedded atom-centered symmetry functions. Phys Chem Chem Phys 2022; 24:23082-23088. [PMID: 36134471 DOI: 10.1039/d2cp02648a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Efficient prediction of the partition coefficient (log P) between polar and non-polar phases could shorten the cycle of drug and materials design. In this work, a descriptor, named 〈q - ACSFs〉conf, is proposed to take the explicit polarization effects in the polar phase and the conformation ensemble of energetic and entropic significance in the non-polar phase into consideration. The polarization effects are involved by embedding the partial charge directly derived from force fields or quantum chemistry calculations into the atom-centered symmetry functions (ACSFs), together with the entropy effects, which are averaged according to the Boltzmann distribution of different conformations taken from the similarity matrix. The model was trained with high-dimensional neural networks (HDNNs) on a public dataset PhysProp (with 41 039 samples). Satisfactory log P prediction performance was achieved on three other datasets, namely, Martel (707 molecules), Star & Non-Star (266) and Huuskonen (1870). The present 〈q - ACSFs〉conf model was also applicable to n-carboxylic acids with the number of carbons ranging from 2 to 14 and 54 kinds of organic solvent. It is easy to apply the present method to arbitrary sized systems and give a transferable atom-based partition coefficient.
Collapse
Affiliation(s)
- Qiang Zhu
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education Institute of Theoretical and Computational Chemistry School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, 210023, P. R. China.
| | - Qingqing Jia
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education Institute of Theoretical and Computational Chemistry School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, 210023, P. R. China.
| | - Ziteng Liu
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education Institute of Theoretical and Computational Chemistry School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, 210023, P. R. China.
| | - Yang Ge
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education Institute of Theoretical and Computational Chemistry School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, 210023, P. R. China.
| | - Xu Gu
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education Institute of Theoretical and Computational Chemistry School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, 210023, P. R. China.
| | - Ziyi Cui
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education Institute of Theoretical and Computational Chemistry School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, 210023, P. R. China.
| | - Mengting Fan
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education Institute of Theoretical and Computational Chemistry School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, 210023, P. R. China.
| | - Jing Ma
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education Institute of Theoretical and Computational Chemistry School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, 210023, P. R. China.
| |
Collapse
|
5
|
Prediction of Partition Coefficients in SDS Micelles by DFT Calculations. Symmetry (Basel) 2021. [DOI: 10.3390/sym13091750] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
A computational methodology using Density-Functional Theory (DFT) calculations was developed to determine the partition coefficient of a compound in a solution of Sodium Dodecyl Sulfate (SDS) micelles. Different sets of DFT calculations were used to predict the free energy of a set of 63 molecules in 15 different solvents with the purpose of identifying the solvents with similar physicochemical characteristics to the studied micelles. Experimental partition coefficients were obtained from Micellar Electrokinetic Chromatography (MEKC) measurements. The experimental partition coefficient of these molecules was compared with the predicted partition coefficient in heptane/water, cyclohexane/water, N-dodecane/water, pyridine/water, acetic acid/water, decan-1-ol/water, octanol/water, propan-2-ol/water, acetone/water, propan-1-ol/water, methanol/water, 1,2-ethane diol/water, dimethyl sulfoxide/water, formic acid/water, and diethyl sulphide/water systems. It is observed that the combination of pronan-1-ol/water solvent was the most appropriated to estimate the partition coefficient for SDS micelles. This approach allowed us to estimate the partition coefficient orders of magnitude faster than the classical molecular dynamics simulations. The DFT calculations were carried out with the well-known exchange correlation functional B3LYP and with the global hybrid functional M06-2X from the Minnesota functionals with 6-31++G ** basis set. The effect of solvation was considered by the continuum model based on density (SMD). The proposed workflow for the prediction rate of the participation coefficient unveiled the symmetric balance between the experimental data and the computational methods.
Collapse
|
6
|
SAMPL7 blind challenge: quantum-mechanical prediction of partition coefficients and acid dissociation constants for small drug-like molecules. J Comput Aided Mol Des 2021; 35:841-851. [PMID: 34164769 DOI: 10.1007/s10822-021-00402-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Accepted: 06/17/2021] [Indexed: 02/02/2023]
Abstract
The physicochemical properties of a drug molecule determine the therapeutic effectiveness of the drug. Thus, the development of fast and accurate theoretical approaches for the prediction of such properties is inevitable. The participation to the SAMPL7 challenge is based on the estimation of logP coefficients and pKa values of small drug-like sulfonamide derivatives. Thereby, quantum mechanical calculations were carried out in order to calculate the free energy of solvation and the transfer energy of 22 drug-like compounds in different environments (water and n-octanol) by employing the SMD solvation model. For logP calculations, we studied eleven different methodologies to calculate the transfer free energies, the lowest RMSE value was obtained for the M06L/def2-TZVP//M06L/def2-SVP level of theory. On the other hand, we employed an isodesmic reaction scheme within the macro pKa framework; this was based on selecting reference molecules similar to the SAMPL7 challenge molecules. Consequently, highly well correlated pKa values were obtained with the M062X/6-311+G(2df,2p)//M052X/6-31+G(d,p) level of theory.
Collapse
|
7
|
Patel P, Kuntz DM, Jones MR, Brooks BR, Wilson AK. SAMPL6 logP challenge: machine learning and quantum mechanical approaches. J Comput Aided Mol Des 2020; 34:495-510. [PMID: 32002780 PMCID: PMC10817701 DOI: 10.1007/s10822-020-00287-0] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2019] [Accepted: 01/08/2020] [Indexed: 10/25/2022]
Abstract
Two different types of approaches: (a) approaches that combine quantitative structure activity relationships, quantum mechanical electronic structure methods, and machine-learning and, (b) electronic structure vertical solvation approaches, were used to predict the logP coefficients of 11 molecules as part of the SAMPL6 logP blind prediction challenge. Using electronic structures optimized with density functional theory (DFT), several molecular descriptors were calculated for each molecule, including van der Waals areas and volumes, HOMO/LUMO energies, dipole moments, polarizabilities, and electrophilic and nucleophilic superdelocalizabilities. A multilinear regression model and a partial least squares model were used to train a set of 97 molecules. As well, descriptors were generated using the molecular operating environment and used to create additional machine learning models. Electronic structure vertical solvation approaches considered include DFT and the domain-based local pair natural orbital methods combined with the solvated variant of the correlation consistent composite approach.
Collapse
Affiliation(s)
- Prajay Patel
- Department of Chemistry, Michigan State University, East Lansing, MI, 48824-1322, USA
| | - David M Kuntz
- Department of Chemistry and Center for Advanced Scientific Computing and Modeling (CASCaM), University of North Texas, Denton, TX, 76203-5070, USA
| | - Michael R Jones
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD, 20852-5690, USA
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD, 20852-5690, USA
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, MI, 48824-1322, USA.
- Department of Chemistry and Center for Advanced Scientific Computing and Modeling (CASCaM), University of North Texas, Denton, TX, 76203-5070, USA.
| |
Collapse
|
8
|
Jones MR, Brooks BR. Quantum chemical predictions of water-octanol partition coefficients applied to the SAMPL6 logP blind challenge. J Comput Aided Mol Des 2020; 34:485-493. [PMID: 32002778 DOI: 10.1007/s10822-020-00286-1] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2019] [Accepted: 01/08/2020] [Indexed: 11/30/2022]
Abstract
Theoretical approaches for predicting physicochemical properties are valuable tools for accelerating the drug discovery process. In this work, quantum chemical methods are used to predict water-octanol partition coefficients as a part of the SAMPL6 blind challenge. The SMD continuum solvent model was employed with MP2 and eight DFT functionals in conjunction with correlation consistent basis sets to determine the water-octanol transfer free energy. Several tactics towards improving the predictions of the partition coefficient were examined, including increasing the quality of basis sets, considering tautomerization, and accounting for inhomogeneities in the water and n-octanol phases. Evaluation of these various schemes highlights the impact of modeling approaches across different methods. With the inclusion of tautomers and adjustments to the permittivity constants, the best predictions were obtained with smaller basis sets and the O3LYP functional, which yielded an RMSE of 0.79 logP units. The results presented correspond to the SAMPL6 logP submission IDs: DYXBT, O7DJK, and AHMTF.
Collapse
Affiliation(s)
- Michael R Jones
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD, 20892-5690, USA.
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD, 20892-5690, USA
| |
Collapse
|
9
|
Ding Y, Xu Y, Qian C, Chen J, Zhu J, Huang H, Shi Y, Huang J. Predicting partition coefficients of drug-like molecules in the SAMPL6 challenge with Drude polarizable force fields. J Comput Aided Mol Des 2020; 34:421-435. [PMID: 31960252 DOI: 10.1007/s10822-020-00282-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2019] [Accepted: 01/08/2020] [Indexed: 12/14/2022]
Abstract
The water-octanol partition coefficient is an important physicochemical property for small molecule drug design. Here, we report our participation in the SAMPL6 logP prediction challenge with free energy perturbation (FEP) calculations in the water phase and in the 1-octanol phase using Drude polarizable force fields. Root mean square error (RMSE) and mean absolute error (MAE) of our prediction are equal to 1.85 and 1.25 logP units. The errors are not evenly distributed. Out of eleven SAMPL6 solutes, FEP/Drude performed very badly on three molecules (deviations all larger than 2 logP units) but good on the remaining eight (deviations all less than 1 logP unit). We find while FEP converges well within one nanosecond in water, simulations in 1-octanol need much longer simulation time and possibly more independent runs for sampling. We also find out that 1-octanol, albeit being a non-polar solvent, still polarizes solute molecules and forms stable hydrogen bonds with them. At the end, we attempt to reweight FEP trajectories with QM/Drude calculations and discuss possible caveats in our simulation setup.
Collapse
Affiliation(s)
- Ye Ding
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - You Xu
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Cheng Qian
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Jinfeng Chen
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Jian Zhu
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Houhou Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Yi Shi
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China
| | - Jing Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China. .,Institute of Biology, Westlake Institute for Advanced Study, 18 Shilongshan Road, Hangzhou, 310024, Zhejiang Province, China.
| |
Collapse
|
10
|
Guan D, Lui R, Matthews S. LogP prediction performance with the SMD solvation model and the M06 density functional family for SAMPL6 blind prediction challenge molecules. J Comput Aided Mol Des 2020; 34:511-522. [PMID: 31939103 DOI: 10.1007/s10822-020-00278-1] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2019] [Accepted: 01/08/2020] [Indexed: 12/13/2022]
Abstract
This work presents a quantum mechanical model for predicting octanol-water partition coefficients of small protein-kinase inhibitor fragments as part of the SAMPL6 LogP Prediction Challenge. The model calculates solvation free energy differences using the M06-2X functional with SMD implicit solvation and the def2-SVP basis set. This model was identified as dqxk4 in the SAMPL6 Challenge and was the third highest performing model in the physical methods category with 0.49 log Root Mean Squared Error (RMSE) for predicting the 11 compounds in SAMPL6 blind prediction set. We also collaboratively investigated the use of empirical models to address model deficiencies for halogenated compounds at minimal additional computational cost. A mixed model consisting of the dqxk4 physical and hdpuj empirical models found improved performance at 0.34 log RMSE on the SAMPL6 dataset. This collaborative mixed model approach shows how empirical models can be leveraged to expediently improve performance in chemical spaces that are difficult for ab initio methods to simulate.
Collapse
Affiliation(s)
- Davy Guan
- Pharmacoinformatics Laboratory, Discipline of Pharmacology, School of Medical Sciences, Faculty of Medicine and Health, The University of Sydney, Sydney, NSW, 2006, Australia
| | - Raymond Lui
- Pharmacoinformatics Laboratory, Discipline of Pharmacology, School of Medical Sciences, Faculty of Medicine and Health, The University of Sydney, Sydney, NSW, 2006, Australia
| | - Slade Matthews
- Pharmacoinformatics Laboratory, Discipline of Pharmacology, School of Medical Sciences, Faculty of Medicine and Health, The University of Sydney, Sydney, NSW, 2006, Australia.
| |
Collapse
|
11
|
Nedyalkova MA, Madurga S, Tobiszewski M, Simeonov V. Calculating the Partition Coefficients of Organic Solvents in Octanol/Water and Octanol/Air. J Chem Inf Model 2019; 59:2257-2263. [PMID: 31042037 DOI: 10.1021/acs.jcim.9b00212] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Partition coefficients define how a solute is distributed between two immiscible phases at equilibrium. The experimental estimation of partition coefficients in a complex system can be an expensive, difficult, and time-consuming process. Here a computational strategy to predict the distributions of a set of solutes in two relevant phase equilibria is presented. The octanol/water and octanol/air partition coefficients are predicted for a group of polar solvents using density functional theory (DFT) calculations in combination with a solvation model based on density (SMD) and are in excellent agreement with experimental data. Thus, the use of quantum-chemical calculations to predict partition coefficients from free energies should be a valuable alternative for unknown solvents. The obtained results indicate that the SMD continuum model in conjunction with any of the three DFT functionals (B3LYP, M06-2X, and M11) agrees with the observed experimental values. The highest correlation to experimental data for the octanol/water partition coefficients was reached by the M11 functional; for the octanol/air partition coefficient, the M06-2X functional yielded the best performance. To the best of our knowledge, this is the first computational approach for the prediction of octanol/air partition coefficients by DFT calculations, which has remarkable accuracy and precision.
Collapse
Affiliation(s)
- Miroslava A Nedyalkova
- Inorganic Chemistry Department, Faculty of Chemistry and Pharmacy , University of Sofia , Sofia 1164 , Bulgaria
| | - Sergio Madurga
- Departament de Ciència de Materials i Química Física and Institut de Química Teòrica i Computacional (IQTCUB) , Universitat de Barcelona , 08028 Barcelona , Catalonia , Spain
| | - Marek Tobiszewski
- Department of Analytical Chemistry, Faculty of Chemistry , Gdańsk University of Technology (GUT) , 80-233 Gdańsk , Poland
| | - Vasil Simeonov
- Analytical Chemistry Department, Faculty of Chemistry and Pharmacy , University of Sofia , Sofia 1164 , Bulgaria
| |
Collapse
|
12
|
Chen X, He S, Liang Z, Li QX, Yan H, Hu J, Liu X. Biodegradation of pyraclostrobin by two microbial communities from Hawaiian soils and metabolic mechanism. JOURNAL OF HAZARDOUS MATERIALS 2018; 354:225-230. [PMID: 29753191 DOI: 10.1016/j.jhazmat.2018.04.067] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/03/2018] [Revised: 04/08/2018] [Accepted: 04/25/2018] [Indexed: 06/08/2023]
Abstract
Pyraclostrobin has been widely and long-termly applicated to agricultural fields. The removal of pyraclostrobin from ecological environment has received wide attention. In this study, using sequential enrichments with pyraclostrobin as a sole carbon source, two microbial communities (HI2 and HI6) capable of catabolizing pyraclostrobin were obtained from Hawaiian soils. The microfloras analysis indicated that only Proteobacteria and Bacteroides could survive in HI2-soil after acclimatization, whereas the number of Proteobacteria in HI6-soil accounted for more than 99%. The percentages of Pseudomonas in the HI2 and HI6 microfloras were 69.3% and 59.3%, respectively. More than 99% of pyraclostrobin (C0 = 100 mg L-1) was degraded by the HI2 and HI6 microorganisms within five days. A unique metabolite was identified by high performance liquid chromatography tandem quadrupole time-of-flight mass spectrometry (HPLC-QTOF-MS/MS). A metabolic pathway involving carbamate hydrolysis was proposed. The tertiary amine group of pyraclostrobin was hydrolyzed to primary amine group with the decarboxylation, which facilitated pyraclostrobin detoxification because carboxylester was an important functional group. The metabolic mechanism suggested that Pseudomonas expressing carboxylesterase might be able to degrade carbamate chemicals. Therefore, Pseudomonas might be an ideal candidate for expression and cloning of carbamate-degrading gene in genomics studies. The current study would have important implications in detoxification and bioremediation of carbamates through the CN bond cleavage of methyl carbamate.
Collapse
Affiliation(s)
- Xiaoxin Chen
- College of Chemistry and Environmental Science, Hebei University, Baoding City, Hebei Province, 071002, PR China.
| | - Sheng He
- College of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing, 100083, PR China.
| | - Zhibin Liang
- Department of Molecular Biosciences and Bioengineering, University of Hawaii at Manoa, Honolulu, HI, 96822, USA.
| | - Qing X Li
- Department of Molecular Biosciences and Bioengineering, University of Hawaii at Manoa, Honolulu, HI, 96822, USA.
| | - Hai Yan
- College of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing, 100083, PR China.
| | - Jiye Hu
- College of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing, 100083, PR China.
| | - Xiaolu Liu
- College of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing, 100083, PR China.
| |
Collapse
|
13
|
König G, Reetz MT, Thiel W. 1-Butanol as a Solvent for Efficient Extraction of Polar Compounds from Aqueous Medium: Theoretical and Practical Aspects. J Phys Chem B 2018; 122:6975-6988. [PMID: 29897756 DOI: 10.1021/acs.jpcb.8b02877] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The extraction of polar molecules from aqueous solution is a challenging task in organic synthesis. 1-Butanol has been used sporadically as an eluent for polar molecules, but it is unclear which molecular features drive its efficiency. Here, we employ free energy simulations to study the partitioning of 15 solutes between water and 1-butanol. The simulations demonstrate that the high affinity of polar molecules to the wet 1-butanol phase is associated with its nanostructure. Small inverse micelles of water are able to accommodate polar solutes and locally mimic an aqueous environment. We verify the simulations based on partition coefficients between water and 1-octanol, and include a blind prediction of the water/1-butanol partition coefficient of cyclohexane-1,2-diol. The calculations are in excellent agreement with experiment, reaching root-mean-square deviations below 0.7 kcal/mol. Actual extractions of cyclohexane-1,2-diol from buffer solutions that mimic cell lysates and suspensions in biocatalytic reactions further exemplify our findings. The yields highlight that extractions with 1-butanol can be significantly more efficient than the conventional protocol based on ethyl acetate.
Collapse
Affiliation(s)
- Gerhard König
- Max-Planck-Institut für Kohlenforschung , 45470 Mülheim an der Ruhr , Germany.,Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research, and Department of Chemistry and Chemical Biology , Rutgers University , Piscataway , New Jersey 08854 , United States
| | - Manfred T Reetz
- Max-Planck-Institut für Kohlenforschung , 45470 Mülheim an der Ruhr , Germany.,Department of Chemistry , Philipps-University Marburg , 35032 Marburg , Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung , 45470 Mülheim an der Ruhr , Germany
| |
Collapse
|
14
|
Fizer O, Fizer M, Sidey V, Studenyak Y, Mariychuk R. Benchmark of different charges for prediction of the partitioning coefficient through the hydrophilic/lipophilic index. J Mol Model 2018; 24:141. [PMID: 29855716 DOI: 10.1007/s00894-018-3692-x] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2018] [Accepted: 05/22/2018] [Indexed: 10/14/2022]
Abstract
A few different theoretical methods for assigning the partial atomic charges were benchmarked for calculation of the hydrophilic/lipophilic index (HLI). The coefficients were selected to produce the best correlation of the HLI values with the experimental octanol-water partition. Different parameters were checked in calculations of partial charges to get the best performance of the HLI values obtained. Thus, four partitioning schemes (Coulson, Mulliken, Merz-Kollman, Ford-Wang) were benchmarked for calculations of atomic charges with six semiempirical methods (AM1, PM3, RM1, PM6, PM6-D3H4, PM7). Moreover, five distinct types of partial atomic charges (Mulliken, Hirshfeld, Löwdin, CHELPG, NPA), obtained at the Hartree-Fock and DFT levels of theory with three basis sets, were tested for their ability to produce the HLI values with the best correlation to experimental logP coefficients of 50 mono-charged organic anions. In the case of the semiempirical methods, the best correlation between the HLI and logP values (the correlation coefficient r = 0.9216) was obtained with the AM1 Ford-Wang parametric electrostatic potential charges. The Mulliken and Coulson charges calculated with the PM7 method can be used as an alternative to AM1, with the r values of 0.9107 and 0.8984, respectively. In the case of the DFT, the PBE/def2-TZVP natural population analysis charges produce the best correlation (r = 0.9220). Nevertheless, in spite of a marginally lower performance (r = 0.9159), the NPA charges computed at the PBE/def2-SVP level are more robust and can be regarded as the optimum choice for calculating the HLI values. Graphical abstract The hydrophilic/lipophilic index (HLI).
Collapse
Affiliation(s)
- Oksana Fizer
- Faculty of Chemistry, Uzhhorod National University, Pidhirna Str. 46, Uzhhorod, 88000, Ukraine.,Faculty of Humanity and Natural Sciences, University of Prešov in Prešov, 17th November 1, Prešov, 08116, Slovak Republic
| | - Maksym Fizer
- Faculty of Chemistry, Uzhhorod National University, Pidhirna Str. 46, Uzhhorod, 88000, Ukraine.
| | - Vasyl Sidey
- Research Institute for Physics and Chemistry of Solid State, Uzhhorod National University, Pidhirna Str. 46, Uzhhorod, 88000, Ukraine
| | - Yaroslav Studenyak
- Faculty of Chemistry, Uzhhorod National University, Pidhirna Str. 46, Uzhhorod, 88000, Ukraine
| | - Ruslan Mariychuk
- Faculty of Humanity and Natural Sciences, University of Prešov in Prešov, 17th November 1, Prešov, 08116, Slovak Republic
| |
Collapse
|
15
|
Bannan CC, Burley KH, Chiu M, Shirts MR, Gilson MK, Mobley DL. Blind prediction of cyclohexane-water distribution coefficients from the SAMPL5 challenge. J Comput Aided Mol Des 2016; 30:927-944. [PMID: 27677750 PMCID: PMC5209301 DOI: 10.1007/s10822-016-9954-8] [Citation(s) in RCA: 87] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2016] [Accepted: 08/22/2016] [Indexed: 12/14/2022]
Abstract
In the recent SAMPL5 challenge, participants submitted predictions for cyclohexane/water distribution coefficients for a set of 53 small molecules. Distribution coefficients (log D) replace the hydration free energies that were a central part of the past five SAMPL challenges. A wide variety of computational methods were represented by the 76 submissions from 18 participating groups. Here, we analyze submissions by a variety of error metrics and provide details for a number of reference calculations we performed. As in the SAMPL4 challenge, we assessed the ability of participants to evaluate not just their statistical uncertainty, but their model uncertainty-how well they can predict the magnitude of their model or force field error for specific predictions. Unfortunately, this remains an area where prediction and analysis need improvement. In SAMPL4 the top performing submissions achieved a root-mean-squared error (RMSE) around 1.5 kcal/mol. If we anticipate accuracy in log D predictions to be similar to the hydration free energy predictions in SAMPL4, the expected error here would be around 1.54 log units. Only a few submissions had an RMSE below 2.5 log units in their predicted log D values. However, distribution coefficients introduced complexities not present in past SAMPL challenges, including tautomer enumeration, that are likely to be important in predicting biomolecular properties of interest to drug discovery, therefore some decrease in accuracy would be expected. Overall, the SAMPL5 distribution coefficient challenge provided great insight into the importance of modeling a variety of physical effects. We believe these types of measurements will be a promising source of data for future blind challenges, especially in view of the relatively straightforward nature of the experiments and the level of insight provided.
Collapse
Affiliation(s)
- Caitlin C Bannan
- Department of Chemistry, University of California, 147 Bison Modular, Irvine, CA, 92697, USA
| | - Kalistyn H Burley
- Department of Pharmaceutical Sciences, University of California, 147 Bison Modular, Irvine, CA, 92697, USA
| | - Michael Chiu
- Qualcomm Institute, University of California, San Diego, CA, USA
| | - Michael R Shirts
- Chemical and Biological Engineering, University of Colorado, Boulder, CO, USA
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, CA, USA
| | - David L Mobley
- Department of Chemistry, University of California, 147 Bison Modular, Irvine, CA, 92697, USA.
- Department of Pharmaceutical Sciences, University of California, 147 Bison Modular, Irvine, CA, 92697, USA.
| |
Collapse
|