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 (r ² = 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 (r ² = 0.80, s = 0.68, n = 174) and a set of agrochemicals (r ² = 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.     

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 for a diverse set of organic compounds based on molecular topology

An accurate and generally applicable method for estimating aqueous solubilities for a diverse set of 1297 organic compounds based on multilinear regression and artificial neural network modeling was developed. Molecular connectivity, shape, and atom-type electrotopological state (E-state) indices were used as structural parameters. The data set was divided into a training set of 884 compounds and a randomly chosen test set of 413 compounds. The structural parameters in a 30-12-1 artificial neural network included 24 atom-type E-state indices and six other topological indices, and for the test set, a predictive r2 = 0.92 and s = 0.60 were achieved. With the same parameters the statistics in the multilinear regression were r2 = 0.88 and s = 0.71, respectively.     

QSPR modelling of the soil sorption coefficient from training sets of different sizes

Quantitative structure–property relationship (QSPR) modelling has been used in many scientific fields. This approach has been extensively applied in environmental research to predict physicochemical properties of compounds with potential environmental impact. The soil sorption coefficient is an important parameter for the evaluation of environmental risks, and it helps to determine the final fate of substances in the environment. In the last few years, different QSPR models have been developed for the determination of the sorption coefficient. In this study, several QSPR models were generated and evaluated for the prediction of log K_oc from the relationship with log P. These models were obtained from an extensive and diverse training set (n = 639) and from subsets of this initial set (i.e. halves, fourths and eighths). The aim of this study was to investigate whether the size of the training set affects the statistical quality of the obtained models. Furthermore, statistical equivalence was verified between the models obtained from smaller sets and the model obtained from the total training set. The results confirmed the equivalence between the models, thus indicating the possibility of using smaller training sets without compromising the statistical quality and predictive capability, as long as most chemical classes in the test set are represented in the training set.     

Adsorption of a diverse set of organic vapors on the bulk water surface

The bulk water surface is of fundamental interest to physical as well as environmental chemistry. As there is a lack of wide-ranging adsorption data from the air to the bulk water surface, a large and diverse data set of adsorption coefficients of nonionic, organic compounds has been produced with inverse gas chromatography. The 61 compounds were chosen to cover a large range of properties, considering the intermolecular interactions between the compounds and the bulk water surface, i.e., van der Waals and electron-donor/acceptor interactions. The data set gained in this work was interpreted with a linear free energy relationship (LFER) based on these intermolecular interactions. From this LFER, a general adsorption model is derived, including compound (i) and surface (surf) properties: log K(i surf/air)(m(3)/m(2)) = 0.135(+/-0.003) log K(i hexadecane/air)(gamma(surf)(vdW))(0.5) + 5.11(+/-0.15)Sigma beta(i2)(H)EA(surf) + 3.60(+/-0.28)Sigma alpha(i2)(H)ED(surf) - 8.47. This adsorption model can be used for the characterization of adsorption to any other surface. The adsorption model for bulk water surface adsorption as well as the general adsorption model can be used as prediction tools in natural or technical systems.     

Comparative molecular field analysis (CoMFA) model using a large diverse set of natural, synthetic and environmental chemicals for binding to the androgen receptor

A large number of natural, synthetic and environmental chemicals are capable of disrupting the endocrine systems of experimental animals, wildlife and humans. These so-called endocrine disrupting chemicals (EDCs), some mimic the functions of the endogenous androgens, have become a concern to the public health. Androgens play an important role in many physiological processes, including the development and maintenance of male sexual characteristics. A common mechanism for androgen to produce both normal and adverse effects is binding to the androgen receptor (AR). In this study, we used Comparative Molecular Field Analysis (CoMFA), a three-dimensional quantitative structure-activity relationship (3D-QSAR) technique, to examine AR-ligand binding affinities. A CoMFA model with r2 = 0.902 and q2 = 0.571 was developed using a large training data set containing 146 structurally diverse natural, synthetic, and environmental chemicals with a 10(6)-fold range of relative binding affinity (RBA). By comparing the binding characteristics derived from the CoMFA contour map with these observed in a human AR crystal structure, we found that the steric and electrostatic properties encoded in this training data set are necessary and sufficient to describe the RBA of AR ligands. Finally, the CoMFA model was challenged with an external test data set; the predicted results were close to the actual values with average difference of 0.637 logRBA. This study demonstrates the utility of this CoMFA model for real-world use in predicting the AR binding affinities of structurally diverse chemicals over a wide RBA range.     

Chemical structure of soil organic matter

Journal of Thermal Analysis and Calorimetry - Selective preservation belongs among the important stabilization mechanisms of soil organic matter (SOM). Conceptually, it is based on non-covalent...     

Prediction of hydration free energies for the SAMPL4 diverse set of compounds using molecular dynamics simulations with the OPLS-AA force field

All-atom molecular dynamics computer simulations were used to blindly predict the hydration free energies of a range of small molecules as part of the SAMPL4 challenge. Compounds were parametrized on the basis of the OPLS-AA force field using three different protocols for deriving partial charges: (1) using existing OPLS-AA atom types and charges with minor adjustments of partial charges on equivalent connecting atoms and derivation of new parameters for a number of distinct chemical groups (N-alkyl imidazole, nitrate) that were not present in the published force field; (2) calculation of quantum mechanical charges via geometry optimization, followed by electrostatic potential (ESP) fitting, using Jaguar at the LMP2/cc-pVTZ(-F) level; and (3) via geometry optimization and CHelpG charges (Gaussian09 at the HF/6-31G* level), followed by two-stage RESP fitting. The absolute hydration free energy was computed by an established protocol including alchemical free energy perturbation with thermodynamic integration. The use of standard OPLS-AA charges (protocol 1) with a number of newly parametrized charges and the use of histidine derived parameters for imidazole yielded an overall root mean square deviation of the prediction from the experimental data of 1.75 kcal/mol. The precision of our results appears to be mainly limited by relatively poor reproducibility of the Lennard-Jones contribution towards the solvation free energy, for which we observed large variability that could be traced to a strong dependence on the initial system conditions.