ADME evaluation in drug discovery. 2. Prediction of partition coefficient by atom-additive approach based on atom-weighted solvent accessible surface areas

A novel method for the calculations of 1-octanol/water partition coefficient (log P) of organic molecules has been presented here. The method, SLOGP v1.0, estimates the log P values by summing the contribution of atom-weighted solvent accessible surface areas (SASA) and correction factors. Altogether 100 atom/group types were used to classify atoms with different chemical environments, and two correlation factors were used to consider the intermolecular hydrophobic interactions and intramolecular hydrogen bonds. Coefficient values for 100 atom/group and two correction factors have been derived from a training set of 1850 compounds. The parametrization procedure for different kinds of atoms was performed as follows: first, the atoms in a molecule were defined to different atom/group types based on SMARTS language, and the correction factors were determined by substructure searching; then, SASA for each atom/group type was calculated and added; finally, multivariate linear regression analysis was applied to optimize the hydrophobic parameters for different atom/group types and correction factors in order to reproduce the experimental log P. The correlation based on the training set gives a model with the correlation coefficient (r) of 0.988, the standard deviation (SD) of 0.368 log units, and the absolute unsigned mean error of 0.261. Comparison of various procedures of log P calculations for the external test set of 138 organic compounds demonstrates that our method bears very good accuracy and is comparable or even better than the fragment-based approaches. Moreover, the atom-additive approach based on SASA was compared with the simple atom-additive approach based on the number of atoms. The calculated results show that the atom-additive approach based on SASA gives better predictions than the simple atom-additive one. Due to the connection between the molecular conformation and the molecular surface areas, the atom-additive model based on SASA may be a more universal model for log P estimation especially for large molecules.     

基于分子表面的水化自由能预测方法

报导了一种基于加权原子表面的水合自由能预测（SAWSA）.对于不同原子类型的溶剂化参数，其参数化分为三个步骤：首先用SMARTS 语言确定不同的原子类型；然后计算每个原子的溶剂可及化表面；最后用遗传算法来优化不同原子类型的溶剂化参数.采用该模型，计算了18个蛋白质分子的水合自由能，预测结果和PB/SA的计算结果呈现了很好的线性关系（r=0.99).计算表明，SAWSA模型对有机小分子和生物大分子都具有很好的预测能力.     

一种计算水溶解度的经验加合模型的适用范围与局限

我们发展了一种用于预测有机小分子化合物水溶解度(logS)的经验方法XLOGS. 它本质上是一种加合模型, 采用83种原子/基团类型和3个校正因子作为模型的描述符.该方法还可以根据一个合适的参照分子的logS实验值来计算未知化合物的logS值. 我们将XLOGS模型在由4171个化合物组成的训练集上进行了参数化, 多元线性回归获得的相关系数R²和标准偏差SD分别为0.82和0.96单位. 将该训练集进一步分为仅含液体化合物和仅含固体化合物的两个子集. XLOGS模型在这两个子集上的回归结果显示前者优于后者(标准偏差分别为0.65单位和0.94单位). 还利用log1/S和logP(脂水分配系数)之间的差值来研究XLOGS方法在液体和固体化合物数据集上的表现. 研究结果表明: XLOGS等加合法模型更适合应用于这两者差值接近于0的化合物. 我们还将XLOGS和其他三种流行的logS计算模型(包括Qikprop, MOE-logS和ALOGPS)在一个含有132个类药化合物的独立测试集上进行了比较. 总体而言, 我们的研究结果为加合法模型在水溶解度预测方面的合理应用提供了指导.     

Estimation of the aqueous solubility of organic molecules by the group contribution approach

Several group contribution methods to estimate the aqueous solubility of organic molecules are proposed and evaluated for their ability to predict the water solubility of new molecules. The learning set consisted of 1168 organic compounds with experimental data taken from the literature after critical evaluation. The best method, based on a new fragment atom scheme, leads to a squared correlation coefficient of 0.95 and an average absolute calculation error of 0.50 log unit, which is superior to other group contribution methods currently available. One of the advantages of this model is that it has upper and lower limits so that the predicted solubilities cannot be unrealistily high or low.     

Prediction of drug solubility from molecular structure using a drug-like training set

Using a training set of 191 drug-like compounds extracted from the AQUASOL database a quantitative structure-property relationship (QSPR) study was conducted employing a set of simple structural and physicochemical properties to predict aqueous solubility. The resultant regression model comprised five parameters (ClogP, molecular weight, indicator variable for aliphatic amine groups, number of rotatable bonds and number of aromatic rings) and demonstrated acceptable statistics (r2 = 0.87, s = 0.51, F = 243.6, n = 191). The model was applied to two test sets consisting of a drug-like set of compounds (r2 = 0.80, s = 0.68, n = 174) and a set of agrochemicals (r2 = 0.88, s = 0.65, n = 200). Using the established general solubility equation (GSE) on the training and drug-like test set gave poorer results than the current study. The agrochemical test set was predicted with equal accuracy using the GSE and the QSPR equation. The results of this study suggest that increasing molecular size, rigidity and lipophilicity decrease solubility whereas increasing conformational flexibility and the presence of a non-conjugated amine group increase the solubility of drug-like compounds. Indeed, the proposed structural parameters make physical sense and provide simple guidelines for modifying solubility during lead optimisation.     

Estimation of aqueous solubility of organic molecules by the group contribution approach. Application to the study of biodegradation.

A reliable and generally applicable aqueous solubility estimation method for organic compounds based on a group contribution approach has been developed. Two models have been established based on two different sets of parameters. One has a higher accuracy, while the other has a more general applicability. The prediction potentials of these two models have been evaluated through cross-validation experiments. For model I, the mean cross-validated r2 and SD for 10 such cross-validation experiments were 0.946 and 0.503 log units, respectively. While for model II, they were 0.953 and 0.546 log units, respectively. Applying our models to estimate the water solubility values for the compounds in an independent test set, we found that model I can be applied to 13 out of 21 compounds with a SD equal to 0.58 log unit and model II can be applied to all the 21 compounds with a SD equal to 1.25 log units. Our models compare favorably to all the current available water estimation methods. A program based on this approach has been written in FORTRAN77 and is currently running on a VAX/VMS system. The program can be applied to estimate the water solubility of the water solubility of any organic chemical with a good or fairly good accuracy except for except for electrolytes. Applying our aqueous solubility estimation models to biodegradation studies, we found that although the water solubility was not the sole factor controlling the rate of biodegradation, ring compounds with greater solubilities were more likely to biodegrade at a faster rate. The significance of the relationship between water solubility and biodegradation activity has been illustrated by predicting the biodegradation activity of 27 new chemicals based solely on their estimated solubility values.     

Anisotropic solvent model of the lipid bilayer. 1. Parameterization of long-range electrostatics and first solvation shell effects

A new implicit solvation model was developed for calculating free energies of transfer of molecules from water to any solvent with defined bulk properties. The transfer energy was calculated as a sum of the first solvation shell energy and the long-range electrostatic contribution. The first term was proportional to solvent accessible surface area and solvation parameters (σ(i)) for different atom types. The electrostatic term was computed as a product of group dipole moments and dipolar solvation parameter (η) for neutral molecules or using a modified Born equation for ions. The regression coefficients in linear dependencies of solvation parameters σ(i) and η on dielectric constant, solvatochromic polarizability parameter π*, and hydrogen-bonding donor and acceptor capacities of solvents were optimized using 1269 experimental transfer energies from 19 organic solvents to water. The root-mean-square errors for neutral compounds and ions were 0.82 and 1.61 kcal/mol, respectively. Quantification of energy components demonstrates the dominant roles of hydrophobic effect for nonpolar atoms and of hydrogen-bonding for polar atoms. The estimated first solvation shell energy outweighs the long-range electrostatics for most compounds including ions. The simplicity and computational efficiency of the model allows its application for modeling of macromolecules in anisotropic environments, such as biological membranes.     

离子液体中有机物溶解度的QSPR模型分析

采用偏最小二乘(PLS)分析方法对84个有机物在四种不同的离子液体中的溶解度进行基于VolSurf参数的定量构性关系(QSPR)研究,取得较好的结果.训练集模型对预测集具有良好的预测能力.参数分析表明有机物具有较大体积的亲水区域,对溶解度有利,且有机物与离子液体之间的相互作用能约为-0.84kJ·mol-1.一定的疏水性对溶解度也是有利因素,当离子液体具有小体积的疏水取代基,有机物具有不对称的局部疏水区域对溶解度有利,当离子液体具有大体积或多个疏水取代基,有机物较高的疏水体积对溶解度有利.多元线性回归(MLR)显示亲水参数W1最重要,表明分子的亲水性是影响有机物在离子液体中溶解的关键因素.     

Molecular dynamics study on hydrophobic effects in aqueous urea solutions.

Hydrophobic effects in aqueous urea were analyzed by molecular dynamics simulations. The contribution of solvents to the potential of mean force between two methane molecules was calculated by using molecular dynamics simulations and was compared with the solubility data of hydrocarbons in aqueous urea. Both the simulation results and the solubility data indicated that urea stabilizes methane-methane association. The stabilization was due to increasing the solvation free energies of small hydrocarbons such as methane by addition of urea. The solvation free energies of larger hydrocarbons, on the other hand, are decreased by addition of urea. This effect of the solute size on hydrophobic free energies in aqueous urea was also analyzed by using molecular dynamics simulations by means of division of the solvation process into two parts: the cavity formation and the introduction of the solute-solvent attractive interactions. In the cavity formation, urea increased hydrophobic free energies, and in the introduction of the solute-solvent attractive interactions, urea decreased hydrophobic free energies. The influence of urea on hydrophobic free energies was determined by the balance of effects of the two parts of the solvation process.     

Atom-type-based AI topological descriptors: application in structure-boiling point correlations of oxo organic compounds

Structure-boiling point relationships are studied for a series of oxo organic compounds by means of multiple linear regression (MLR) analysis. Excellent MLR models based on the recently introduced Xu index and the atom-type-based AI indices are obtained for the two subsets containing respectively 77 ethers and 107 carbonyl compounds and a combined set of 184 oxo compounds. The best models are tested using the leave-one-out cross-validation and an external test set, respectively. The MLR model produces a correlation coefficient of r = 0.9977 and a standard error of s = 3.99 degrees C for the training set of 184 compounds, and r(cv) = 0.9974 and s(cv) = 4.16 degrees C for the cross-validation set, and r(pred) = 0.9949 and s(pred) = 4.38 degrees C for the prediction set of 21 compounds. For the two subsets containing respectively 77 ethers and 107 carbonyl compounds, the quality of the models is further improved. The standard errors are reduced to 3.30 and 3.02 degrees C, respectively. Furthermore, the results obtained from this study indicate that the boiling points of the studied oxo compound dominantly depend on molecular size and also depend on individual atom types, especially oxygen heteroatoms in molecules due to strong polar interactions between molecules. These excellent structure-boiling point models not only provide profound insights into the role of structural features in a molecule but also illustrate the usefulness of these indices in QSPR/QSAR modeling of complex compounds.     

基于原子表面的蛋白质水合自由能预测模型

提出了一种计算蛋白质水合自由能的简化模型（SAWSA 2）.模型把蛋白质分子中的原子分为20种不同的原子类型,通过每类原子的溶剂可及化表面以及相应的溶剂化参数,就可以得到分子的水合自由能.不同原子类型的溶剂化参数通过110个蛋白质分子水合自由能拟合得到,水合自由能的标准值采用了基于求解Possion-Boltzmann方程（PB）以及分子表面计算（SA） 相结合的方法.采用得到的模型,预测了20个蛋白质分子的水合自由能,预测值的相对值和绝对值都能和PB/SA的计算值很好地吻合,大大优于两种已报导的水合自由能模型.