Imputation of sensory properties using deep learning.

Predicting the sensory properties of compounds is challenging due to the subjective nature of the experimental measurements. This testing relies on a panel of human participants and is therefore also expensive and time-consuming. We describe the application of a state-of-the-art deep learning method, Alchemite™, to the imputation of sparse physicochemical and sensory data and compare the results with conventional quantitative structure-activity relationship methods and a multi-target graph convolutional neural network. The imputation model achieved a substantially higher accuracy of prediction, with improvements in R2 between 0.26 and 0.45 over the next best method for each sensory property. We also demonstrate that robust uncertainty estimates generated by the imputation model enable the most accurate predictions to be identified and that imputation also more accurately predicts activity cliffs, where small changes in compound structure result in large changes in sensory properties. In combination, these results demonstrate that the use of imputation, based on data from less expensive, early experiments, enables better selection of compounds for more costly studies, saving experimental time and resources.

Online chemical modeling environment (OCHEM): web platform for data storage, model development and publishing of chemical information.

The Online Chemical Modeling Environment is a web-based platform that aims to automate and simplify the typical steps required for QSAR modeling. The platform consists of two major subsystems: the database of experimental measurements and the modeling framework. A user-contributed database contains a set of tools for easy input, search and modification of thousands of records. The OCHEM database is based on the wiki principle and focuses primarily on the quality and verifiability of the data. The database is tightly integrated with the modeling framework, which supports all the steps required to create a predictive model: data search, calculation and selection of a vast variety of molecular descriptors, application of machine learning methods, validation, analysis of the model and assessment of the applicability domain. As compared to other similar systems, OCHEM is not intended to re-implement the existing tools or models but rather to invite the original authors to contribute their results, make them publicly available, share them with other users and to become members of the growing research community. Our intention is to make OCHEM a widely used platform to perform the QSPR/QSAR studies online and share it with other users on the Web. The ultimate goal of OCHEM is collecting all possible chemoinformatics tools within one simple, reliable and user-friendly resource. The OCHEM is free for web users and it is available online at http://www.ochem.eu.

Predicting inhibitors of acetylcholinesterase by regression and classification machine learning approaches with combinations of molecular descriptors.

PURPOSE: Acetylcholinesterase (AChE) is both a therapeutic target for Alzheimer's disease and a target for organophosphorus, carbamates and chemical warfare agents. Prediction of the likelihood of compounds interacting with this enzyme is therefore important from both therapeutic and toxicological perspectives. MATERIALS AND METHODS: Support vector machine classification and regression models with molecular descriptors derived from Shape Signatures and the Molecular Operating Environment (MOE) application software were built and tested using a set of piperidine AChE inhibitors (N = 110). RESULTS: The combination of the alignment free Shape Signatures and 2D MOE descriptors with the Support Vector Regression method outperforms the models based solely on 2D and internal 3D (i3D) MOE descriptors, and is comparable with the best previously reported PLS model based on CoMFA molecular descriptors (r(2)(test,SVR) = 0.48 vs. r(2)(test,PLS) = 0.47 from Sutherland et al. J Med Chem 47:5541-5554, 2004). Support Vector Classification algorithms proved superior to a classifier based on scores from the molecular docking program GOLD, with the overall prediction accuracies being Q(SVC(10CV)) = 74% and Q(SVC(LNO)) = 67% vs. Q(GOLD) = 56%. CONCLUSIONS: These new machine learning models with combined descriptor schemes may find utility for predicting novel AChE inhibitors.

Shape signatures: new descriptors for predicting cardiotoxicity in silico.

Shape Signatures is a new computational tool that is being evaluated for applications in computational toxicology and drug discovery. The method employs a customized ray-tracing algorithm to explore the volume enclosed by the surface of a molecule and then uses the output to construct compact histograms (i.e., signatures) that encode for molecular shape and polarity. In the present study, we extend the application of the Shape Signatures methodology to the domain of computational models for cardiotoxicity. The Shape Signatures method is used to generate molecular descriptors that are then utilized with widely used classification techniques such as k nearest neighbors ( k-NN), support vector machines (SVM), and Kohonen self-organizing maps (SOM). The performances of these approaches were assessed by applying them to a data set of compounds with varying affinity toward the 5-HT(2B) receptor as well as a set of human ether-a-go-go-related gene (hERG) potassium channel inhibitors. Our classification models for 5-HT(2B) represented the first attempt at global computational models for this receptor and exhibited average accuracies in the range of 73-83%. This level of performance is comparable to using commercially available molecular descriptors. The overall accuracy of the hERG Shape Signatures-SVM models was 69-73%, in line with other computational models published to date. Our data indicate that Shape Signatures descriptors can be used with SVM and Kohonen SOM and perform better in classification problems related to the analysis of highly clustered and heterogeneous property spaces. Such models may have utility for predicting the potential for cardiotoxicity in drug discovery mediated by the 5-HT(2B) receptor and hERG.

Hybrid scoring and classification approaches to predict human pregnane X receptor activators.

PURPOSE: The human pregnane X receptor (PXR) is a transcriptional regulator of many genes involved in xenobiotic metabolism and excretion. Reliable prediction of high affinity binders with this receptor would be valuable for pharmaceutical drug discovery to predict potential toxicological responses MATERIALS AND METHODS: Computational models were developed and validated for a dataset consisting of human PXR (PXR) activators and non-activators. We used support vector machine (SVM) algorithms with molecular descriptors derived from two sources, Shape Signatures and the Molecular Operating Environment (MOE) application software. We also employed the molecular docking program GOLD in which the GoldScore method was supplemented with other scoring functions to improve docking results. RESULTS: The overall test set prediction accuracy for PXR activators with SVM was 72% to 81%. This indicates that molecular shape descriptors are useful in classification of compounds binding to this receptor. The best docking prediction accuracy (61%) was obtained using 1D Shape Signature descriptors as a weighting factor to the GoldScore. By pooling the available human PXR data sets we revealed those molecular features that are associated with human PXR activators. CONCLUSIONS: These combined computational approaches using molecular shape information may assist scientists to more confidently identify PXR activators.

New predictive models for blood-brain barrier permeability of drug-like molecules.

PURPOSE: The goals of the present study were to apply a generalized regression model and support vector machine (SVM) models with Shape Signatures descriptors, to the domain of blood-brain barrier (BBB) modeling. MATERIALS AND METHODS: The Shape Signatures method is a novel computational tool that was used to generate molecular descriptors utilized with the SVM classification technique with various BBB datasets. For comparison purposes we have created a generalized linear regression model with eight MOE descriptors and these same descriptors were also used to create SVM models. RESULTS: The generalized regression model was tested on 100 molecules not in the model and resulted in a correlation r2 = 0.65. SVM models with MOE descriptors were superior to regression models, while Shape Signatures SVM models were comparable or better than those with MOE descriptors. The best 2D shape signature models had 10-fold cross validation prediction accuracy between 80-83% and leave-20%-out testing prediction accuracy between 80-82% as well as correctly predicting 84% of BBB+ compounds (n = 95) in an external database of drugs. CONCLUSIONS: Our data indicate that Shape Signatures descriptors can be used with SVM and these models may have utility for predicting blood-brain barrier permeation in drug discovery.

Prediction of Protein Loop Conformations using the AGBNP Implicit Solvent Model and Torsion Angle Sampling.

The OPLS-AA all-atom force field and the Analytical Generalized Born plus Non-Polar (AGBNP) implicit solvent model, in conjunction with torsion angle conformational search protocols based on the Protein Local Optimization Program (PLOP), are shown to be effective in predicting the native conformations of 57 9-residue and 35 13-residue loops of a diverse series of proteins with low sequence identity. The novel nonpolar solvation free energy estimator implemented in AGBNP augmented by correction terms aimed at reducing the occurrence of ion pairing are important to achieve the best prediction accuracy. Extended versions of the previously developed PLOP-based conformational search schemes based on calculations in the crystal environment are reported that are suitable for application to loop homology modeling without the crystal environment. Our results suggest that in general the loop backbone conformation is not strongly influenced by crystal packing. The application of the temperature Replica Exchange Molecular Dynamics (T-REMD) sampling method for a few examples where PLOP sampling is insufficient are also reported. The results reported indicate that the OPLS-AA/AGBNP effective potential is suitable for high-resolution modeling of proteins in the final stages of homology modeling and/or protein crystallographic refinement.

Polyfunctional compounds containing the 4,6-dialkoxy-7-arylthioheptene moiety as synthetically useful intermediates. The course of lewis acid-induced transformations.

Data on the selectivity of the Lewis acids induced transformations of the title compounds are presented, and the routes leading to formation of products containing either cyclohexane or 1,3-diene units are described.