The method for QSAR modelling of rat acute toxicity based on the combination of QNA (Quantitative Neighbourhoods of Atoms) descriptors, PASS (Prediction of Activity Spectra for Substances) predictions and self-consistent regression (SCR) is presented. PASS predicted biological activity profiles are used as independent input variables for QSAR modelling with SCR. QSAR models were developed using LD50 values for compounds tested on rats with four types of administration (oral, intravenous, intraperitoneal, subcutaneous). The proposed method was evaluated on the set of compounds tested for acute rat toxicity with oral administration (7286 compounds) used for testing the known QSAR methods in T.E.S.T. 3.0 program (U.S. EPA). The several other sets of compounds tested for acute rat toxicity by different routes of administration selected from SYMYX MDL Toxicity Database were used too. The method was compared with the results of prediction of acute rodent toxicity for noncongeneric sets obtained by ACD/Labs Inc. The test sets were predicted with regards to the applicability domain. Comparison of accuracy for QSAR models obtained separately using QNA descriptors, PASS predictions, nearest neighbours' assessment with consensus models clearly demonstrated the benefits of consensus prediction. Free available web-service for prediction of LD50 values of rat acute toxicity was developed: http://www.pharmaexpert.ru/GUSAR/AcuToxPredict/.

QSAR Modelling of Rat Acute Toxicity on the Basis of PASS Prediction

Adverse side effects of drug-drug interactions induced by human cytochrome P450 (CYP) inhibition is an important consideration, especially, during the research phase of drug discovery. It is highly desirable to develop computational models that can predict the inhibitive effect of a compound against a specific CYP isoform. In this study, inhibitor predicting models were developed for five major CYP isoforms, namely 1A2, 2C9, 2C19, 2D6, and 3A4, using a combined classifier algorithm on a large data set containing more than 24,700 unique compounds, extracted from PubChem. The combined classifiers algorithm is an ensemble of different independent machine learning classifiers including support vector machine, C4.5 decision tree, k-nearest neighbor, and naive Bayes, fused by a back-propagation artificial neural network (BP-ANN). All developed models were validated by 5-fold cross-validation and a diverse validation set composed of about 9000 diverse unique compounds. The range of the area under the receiver operating characteristic curve (AUC) for the validation sets was 0.764 to 0.815 for CYP1A2, 0.837 to 0.861 for CYP2C9, 0.793 to 0.842 for CYP2C19, 0.839 to 0.886 for CYP2D6, and 0.754 to 0.790 for CYP3A4, respectively, using the new developed combined classifiers. The overall performance of the combined classifiers fused by BP-ANN was superior to that of three classic fusion techniques (Mean, Maximum, and Multiply). The chemical spaces of data sets were explored by multidimensional scaling plots, and the use of applicability domain improved the prediction accuracies of models. In addition, some representative substructure fragments differentiating CYP inhibitors and noninhibitors were characterized by the substructure fragment analysis. These classification models are applicable for virtual screening of the five major CYP isoforms inhibitors or can be used as simple filters of potential chemicals in drug discovery.

Classification of cytochrome P450 inhibitors and noninhibitors using combined classifiers.

The human cytochrome P450 enzymes (P450s) catalyze oxidative reactions of a broad spectrum of substrates and play a critical role in the metabolism of xenobiotics, such as drugs and dietary compounds. CYP3A4 is known to be the main enzyme involved in the metabolism of drugs and most other xenobiotics. Dietary compounds, of which polyphenolics are the most studied, have been shown to interact with CYP3A4 and alter its expression and activity. Traditionally, the liver was considered the prime site of CYP3A-mediated first-pass metabolic extraction, but in vitro and in vivo studies now suggest that the small intestine can be of equal or even greater importance for the metabolism of polyphenolics and drugs. Recent studies have pointed to the role of gut microbiota in the metabolic fate of polyphenolics in human, suggesting their involvement in the complex interactions between dietary polyphenols and CYP3A4. Last but not least, all the above suggests that coadministration of drugs and foods that are rich in polyphenols is expected to stimulate undesirable clinical consequences. This review focuses on interactions between dietary polyphenols and CYP3A4 as they relate to structural considerations, food-drug interactions, and potential negative consequences of interactions between CYP3A4 and polyphenols.

/pdf/interactions-between-cyp3a4-and-dietary-polyphenols-4zdmc3cg0u.pdf

Interactions between CYP3A4 and Dietary Polyphenols

Introduction: With the emergence of the ‘big data’ era, the biomedical research community has great interest in exploiting publicly available chemical information for drug discovery. PubChem is an example of public databases that provide a large amount of chemical information free of charge.Areas covered: This article provides an overview of how PubChem’s data, tools, and services can be used for virtual screening and reviews recent publications that discuss important aspects of exploiting PubChem for drug discovery.Expert opinion: PubChem offers comprehensive chemical information useful for drug discovery. It also provides multiple programmatic access routes, which are essential to build automated virtual screening pipelines that exploit PubChem data. In addition, PubChemRDF allows users to download PubChem data and load them into a local computing facility, facilitating data integration between PubChem and other resources. PubChem resources have been used in many studies for developing bioactivi...

/pdf/getting-the-most-out-of-pubchem-for-virtual-screening-vxg7yw4csj.pdf

Getting the most out of PubChem for virtual screening.

Quantitative structure-activity relationship (QSAR) models have long been used for making predictions and data gap filling in diverse fields including medicinal chemistry, predictive toxicology, environmental fate modeling, materials science, agricultural science, nanoscience, food science, and so forth. Usually a QSAR model is developed based on chemical information of a properly designed training set and corresponding experimental response data while the model is validated using one or more test set(s) for which the experimental response data are available. However, it is interesting to estimate the reliability of predictions when the model is applied to a completely new data set (true external set) even when the new data points are within applicability domain (AD) of the developed model. In the present study, we have categorized the quality of predictions for the test set or true external set into three groups (good, moderate, and bad) based on absolute prediction errors. Then, we have used three criteria [(a) mean absolute error of leave-one-out predictions for 10 most close training compounds for each query molecule; (b) AD in terms of similarity based on the standardization approach; and (c) proximity of the predicted value of the query compound to the mean training response] in different weighting schemes for making a composite score of predictions. It was found that using the most frequently appearing weighting scheme 0.5-0-0.5, the composite score-based categorization showed concordance with absolute prediction error-based categorization for more than 80% test data points while working with 5 different datasets with 15 models for each set derived in three different splitting techniques. These observations were also confirmed with true external sets for another four endpoints suggesting applicability of the scheme to judge the reliability of predictions for new datasets. The scheme has been implemented in a tool "Prediction Reliability Indicator" available at http://dtclab.webs.com/software-tools and http://teqip.jdvu.ac.in/QSAR_Tools/DTCLab/, and the tool is presently valid for multiple linear regression models only.

How Precise Are Our Quantitative Structure–Activity Relationship Derived Predictions for New Query Chemicals?

This study presents a new type of acute toxicity (LD(50)) prediction that enables automated assessment of the reliability of predictions (which is synonymous with the assessment of the Model Applicability Domain as defined by the Organization for Economic Cooperation and Development). Analysis involved nearly 75,000 compounds from six animal systems (acute rat toxicity after oral and intraperitoneal administration; acute mouse toxicity after oral, intraperitoneal, intravenous, and subcutaneous administration). Fragmental Partial Least Squares (PLS) with 100 bootstraps yielded baseline predictions that were automatically corrected for non-linear effects in local chemical spaces--a combination called Global, Adjusted Locally According to Similarity (GALAS) modelling methodology. Each prediction obtained in this manner is provided with a reliability index value that depends on both compound's similarity to the training set (that accounts for similar trends in LD(50) variations within multiple bootstraps) and consistency of experimental results with regard to the baseline model in the local chemical environment. The actual performance of the Reliability Index (RI) was proven by its good (and uniform) correlations with Root Mean Square Error (RMSE) in all validation sets, thus providing quantitative assessment of the Model Applicability Domain. The obtained models can be used for compound screening in the early stages of drug development and prioritization for experimental in vitro testing or later in vivo animal acute toxicity studies.

Estimation of reliability of predictions and model applicability domain evaluation in the analysis of acute toxicity (LD50).

A new structure–activity relationship model predicting the probability for a compound to inhibit human cytochrome P450 3A4 has been developed using data for >800 compounds from various literature sources and tested on PubChem screening data. Novel GALAS (Global, Adjusted Locally According to Similarity) modeling methodology has been used, which is a combination of baseline global QSAR model and local similarity based corrections. GALAS modeling method allows forecasting the reliability of prediction thus defining the model applicability domain. For compounds within this domain the statistical results of the final model approach the data consistency between experimental data from literature and PubChem datasets with the overall accuracy of 89%. However, the original model is applicable only for less than a half of PubChem database. Since the similarity correction procedure of GALAS modeling method allows straightforward model training, the possibility to expand the applicability domain has been investigated. Experimental data from PubChem dataset served as an example of in-house high-throughput screening data. The model successfully adapted itself to both data classified using the same and different IC50 threshold compared with the training set. In addition, adjustment of the CYP3A4 inhibition model to compounds with a novel chemical scaffold has been demonstrated. The reported GALAS model is proposed as a useful tool for virtual screening of compounds for possible drug-drug interactions even prior to the actual synthesis.

Trainable structure-activity relationship model for virtual screening of CYP3A4 inhibition.

This article briefly introduces the results of in silico prediction of the most probable metabolism sites for the human cytochrome P450 3A4 and 2D6 isoforms. Ligand-based QSAR models have been developed using a novel GALAS modeling approach, and provide probabilities of being a target of CYP3A4 or CYP2D6 for any atom in a molecule. The GALAS-model development methodology allows evaluation of the reliability of predictions in the form of estimated prediction Reliability Indices (RIs). For all the models considered in this study, the number of misclassifications and inconclusive results was reduced significantly when only predictions of high quality (RI>0.5) were taken into account, demonstrating that RI reflects accuracy of prediction. The applicability domain of regioselectivity models is shown to be easily expandable to cover compound classes of interest to the user. The results obtained so far show promising perspectives for the utilization of the GALAS modeling in the analysis of regioselectivity for other important biotransformation enzymes – a work currently in progress.

Andrius Sazonovas

Papers

Estimation of reliability of predictions and model applicability domain evaluation in the analysis of acute toxicity (LD50).

Trainable structure-activity relationship model for virtual screening of CYP3A4 inhibition.

Probabilistic prediction of the human CYP3A4 and CYP2D6 metabolism sites.

Acute toxicity (LD50) prediction involving fragmental QSAR model, similarity analysis and reliability of predictions

Physicochemical QSAR analysis of hERG inhibition revisited: towards a quantitative potency prediction