Statistical Analysis and Compound Selection of Combinatorial Libraries for Soluble Epoxide Hydrolase

A Random Forest Model to Predict the Activity of a Large Set of Soluble Epoxide Hydrolase Inhibitors Solely Based on a Set of Simple Fragmental Descriptors

BACKGROUND

The Soluble Epoxide Hydrolase (sEH) is a ubiquitously expressed enzyme in various tissues. The inhibition of the sEH has shown promising results to treat hypertension, alleviate pain and inflammation.

OBJECTIVE

In this study, the power of machine learning has been employed to develop a predictive QSAR model for a large set of sEH inhibitors.

METHODS

In this study, the random forest method was employed to make a valid model for the prediction of sEH inhibition. Besides, two new methods (Treeinterpreter python package and LIME, Local Interpretable Model-agnostic Explanations) have been exploited to explain and interpret the model.

RESULTS

The performance metrics of the model were as follows: R2=0.831, Q2=0.565, RMSE=0.552 and R2 pred=0.595. The model also demonstrated good predictability on the two extra external test sets at least in terms of ranking. The Spearman's rank correlation coefficients for external test set 1 and 2 were 0.872 and 0.673, respectively. The external test set 2 was a diverse one compared to the training set. Therefore, the model could be used for virtual screening to enrich potential sEH inhibitors among a diverse compound library.

CONCLUSION

As the model was solely developed based on a set of simple fragmental descriptors, the model was explained by two local interpretation algorithms, and this could guide medicinal chemists to design new sEH inhibitors. Moreover, the most important general descriptors (fragments) suggested by the model were consistent with the available crystallographic data. The model is available as an executable binary at http://www.pharm-sbg.com and https://github.com/shamsaraj.

Predicting the binding modes and sites of metabolism of xenobiotics

Metabolism studies are an essential integral part of ADMET profiling of drug candidates to evaluate their safety and efficacy. Cytochrome P-450 (CYP) metabolizes a wide variety of xenobiotics/drugs. The binding modes of these compounds with CYP and their intrinsic reactivities decide the metabolic products. We report here a novel computational protocol, which comprises docking of ligands to heme-containing CYPs and prediction of binding energies through a newly developed scoring function, followed by analyses of the docked structures and molecular orbitals of the ligand molecules, for predicting the sites of metabolism (SOM) of ligands. The calculated binding free energies of 121 heme-containing protein-ligand docked complexes yielded a correlation coefficient of 0.84 against experiment. Molecular orbital analyses of the resultant top three unique poses of the docked complexes provided a success rate of 87% in identifying the experimentally known sites of metabolism of the xenobiotics. The SOM prediction methodology is freely accessible at .

Predicting the DPP-IV inhibitory activity pIC₅₀ based on their physicochemical properties

The second development program developed in this work was introduced to obtain physicochemical properties of DPP-IV inhibitors. Based on the computation of molecular descriptors, a two-stage feature selection method called mRMR-BFS (minimum redundancy maximum relevance-backward feature selection) was adopted. Then, the support vector regression (SVR) was used in the establishment of the model to map DPP-IV inhibitors to their corresponding inhibitory activity possible. The squared correlation coefficient for the training set of LOOCV and the test set are 0.815 and 0.884, respectively. An online server for predicting inhibitory activity pIC50 of the DPP-IV inhibitors as described in this paper has been given in the introduction.