Multivariate modeling of cytochrome P450 3A4 inhibition

An Uncertainty-Guided Deep Learning Method Facilitates Rapid Screening of CYP3A4 Inhibitors

Cytochrome P450 3A4 (CYP3A4), a prominent member of the P450 enzyme superfamily, plays a crucial role in metabolizing various xenobiotics, including over 50% of clinically significant drugs. Evaluating CYP3A4 inhibition before drug approval is essential to avoiding potentially harmful pharmacokinetic drug-drug interactions (DDIs) and adverse drug reactions (ADRs). Despite the development of several CYP inhibitor prediction models, the primary approach for screening CYP inhibitors still relies on experimental methods. This might stem from the limitations of existing models, which only provide deterministic classification outcomes instead of precise inhibition intensity (e.g., IC50) and often suffer from inadequate prediction reliability. To address this challenge, we propose an uncertainty-guided regression model to accurately predict the IC50 values of anti-CYP3A4 activities. First, a comprehensive data set of CYP3A4 inhibitors was compiled, consisting of 27,045 compounds with classification labels, including 4395 compounds with explicit IC50 values. Second, by integrating the predictions of the classification model trained on a larger data set and introducing an evidential uncertainty method to rank prediction confidence, we obtained a high-precision and reliable regression model. Finally, we use the evidential uncertainty values as a trustworthy indicator to perform a virtual screening of an in-house compound set. The in vitro experiment results revealed that this new indicator significantly improved the hit ratio and reduced false positives among the top-ranked compounds. Specifically, among the top 20 compounds ranked with uncertainty, 15 compounds were identified as novel CYP3A4 inhibitors, and three of them exhibited activities less than 1 μM. In summary, our findings highlight the effectiveness of incorporating uncertainty in compound screening, providing a promising strategy for drug discovery and development.

Computational prediction of cytochrome P450 inhibition and induction

Cytochrome P450 (CYP) enzymes play an important role in the phase I metabolism of many xenobiotics. Most drug-drug interactions (DDIs) associated with CYP are caused by either CYP inhibition or induction. The early detection of potential DDIs is highly desirable in the pharmaceutical industry because DDIs can cause serious adverse events, which can lead to poor patient health and drug development failures. Recently, many computational studies predicting CYP inhibition and induction have been reported. The current computational modeling approaches for CYP metabolism are classified as ligand- and structure-based; various techniques, such as quantitative structure-activity relationships, machine learning, docking, and molecular dynamic simulation, are involved in both the approaches. Recently, combining these two approaches have resulted in improvements in the prediction accuracy of DDIs. In this review, we present important, recent developments in the computational prediction of the inhibition of four clinically crucial CYP isoforms (CYP1A2, 2C9, 2D6, and 3A4) and three nuclear receptors (aryl hydrocarbon receptor, constitutive androstane receptor, and pregnane X receptor) involved in the induction of CYP1A2, 2B6, and 3A4, respectively.

Prediction of Drug-Drug Interactions Related to Inhibition or Induction of Drug-Metabolizing Enzymes

Drug-drug interaction (DDI) is the phenomenon of alteration of the pharmacological activity of a drug(s) when another drug(s) is co-administered in cases of so-called polypharmacy. There are three types of DDIs: pharmacokinetic (PK), pharmacodynamic, and pharmaceutical. PK is the most frequent type of DDI, which often appears as a result of the inhibition or induction of drug-metabolising enzymes (DME). In this review, we summarise in silico methods that may be applied for the prediction of the inhibition or induction of DMEs and describe appropriate computational methods for DDI prediction, showing the current situation and perspectives of these approaches in medicinal and pharmaceutical chemistry. We review sources of information on DDI, which can be used in pharmaceutical investigations and medicinal practice and/or for the creation of computational models. The problem of the inaccuracy and redundancy of these data are discussed. We provide information on the state-of-the-art physiologically- based pharmacokinetic modelling (PBPK) approaches and DME-based in silico methods. In the section on ligand-based methods, we describe pharmacophore models, molecular field analysis, quantitative structure-activity relationships (QSAR), and similarity analysis applied to the prediction of DDI related to the inhibition or induction of DME. In conclusion, we discuss the problems of DDI severity assessment, mention factors that influence severity, and highlight the issues, perspectives and practical using of in silico methods.

Ligand and structure-based classification models for prediction of P-glycoprotein inhibitors

The ABC transporter P-glycoprotein (P-gp) actively transports a wide range of drugs and toxins out of cells, and is therefore related to multidrug resistance and the ADME profile of therapeutics. Thus, development of predictive in silico models for the identification of P-gp inhibitors is of great interest in the field of drug discovery and development. So far in silico P-gp inhibitor prediction was dominated by ligand-based approaches because of the lack of high-quality structural information about P-gp. The present study aims at comparing the P-gp inhibitor/noninhibitor classification performance obtained by docking into a homology model of P-gp, to supervised machine learning methods, such as Kappa nearest neighbor, support vector machine (SVM), random fores,t and binary QSAR, by using a large, structurally diverse data set. In addition, the applicability domain of the models was assessed using an algorithm based on Euclidean distance. Results show that random forest and SVM performed best for classification of P-gp inhibitors and noninhibitors, correctly predicting 73/75% of the external test set compounds. Classification based on the docking experiments using the scoring function ChemScore resulted in the correct prediction of 61% of the external test set. This demonstrates that ligand-based models currently remain the methods of choice for accurately predicting P-gp inhibitors. However, structure-based classification offers information about possible drug/protein interactions, which helps in understanding the molecular basis of ligand-transporter interaction and could therefore also support lead optimization.

Three-dimensional quantitative structure-activity relationship analysis of inhibitors of human and rat cytochrome P4503A enzymes

Cytochrome P450 3A4 (CYP3A4) is a member of the CYP family and is an important enzyme in drug metabolism. A compound that inhibits CYP3A4 activity could also affect the pharmacokinetics of other substrates, resulting in drug-drug interactions (DDIs) that could cause side effects. Pharmacokinetic data from drug-development studies in rats often determine the dosage used in human clinical trials. It is therefore useful to understand differences in metabolism in different species at an early stage in drug development. Human and rat CYP3A enzymes show different inhibition profiles with different drugs, although the mechanisms involved are not yet clear. Here we built three-dimensional quantitative structure-activity relationship (3D-QSAR) models using structure-based comparative molecular field analysis (CoMFA), to predict the direct inhibitory activity of ligands for human CYP3A4 and rat CYP3A1, based on computer-ligand docking. The alignment of the ligand docking poses suggested that key amino acid-ligand interactions (e.g., Thr309 in CYP3A4 and Pro310 in CYP3A1) characterized the different potencies with which the ligands inhibited CYP3A4 and CYP3A1. The 3D-QSAR models for human and rat CYP3A family inhibitors predicted the potency of inhibitors and could be useful for assessing DDIs at an early stage in drug discovery.

Fingerprint-based in silico models for the prediction of P-glycoprotein substrates and inhibitors

P-Glycoprotein (P-gp, ABCB1) plays a significant role in determining the ADMET properties of drugs and drug candidates. Substrates of P-gp are not only subject to multidrug resistance (MDR) in tumor therapy, they are also associated with poor pharmacokinetic profiles. In contrast, inhibitors of P-gp have been advocated as modulators of MDR. However, due to the polyspecificity of P-gp, knowledge on the molecular basis of ligand-transporter interaction is still poor, which renders the prediction of whether a compound is a P-gp substrate/non-substrate or an inhibitor/non-inhibitor quite challenging. In the present investigation, we used a set of fingerprints representing the presence/absence of various functional groups for machine learning based classification of a set of 484 substrates/non-substrates and a set of 1935 inhibitors/non-inhibitors. Best models were obtained using a combination of a wrapper subset evaluator (WSE) with random forest (RF), kappa nearest neighbor (kNN) and support vector machine (SVM), showing accuracies >70%. Best P-gp substrate models were further validated with three sets of external P-gp substrate sources, which include Drug Bank (n = 134), TP Search (n = 90) and a set compiled from literature (n = 76). Association rule analysis explores the various structural feature requirements for P-gp substrates and inhibitors.

QSAR Modelling of CYP3A4 Inhibition as a Screening Tool in the Context of DrugDrug Interaction Studies

Drugdrug interaction potential (DDI), especially cytochrome P450 (CYP) 3A4 inhibition potential, is one of the most important parameters to be optimized before preclinical and clinical pharmaceutical development as regard to the number of marketed drug metabolized mainly by this CYP and potentially co-administered with the future drug. The present study aims to develop in silico models for CYP3A4 inhibition prediction to help medicinal chemists during the discovery phase and even before the synthesis of new chemical entities (NCEs), focusing on NCEs devoid of any inhibitory potential toward this CYP. In order to find a relevant relationship between CYP3A4 inhibition and chemical features of the screened compounds, we applied a genetic-algorithm-based QSAR exploratory tool SQS (Stochastic QSAR Sampler) in combination with different description approaches comprising alignment-independent Volsurf descriptors, ISIDA fragments and Topological Fuzzy Pharmacophore Triplets. The experimental data used to build models were extracted from an in-house database. We derived a model with good prediction ability that was confirmed on both newly synthesized compound and public dataset retrieved from Pubchem database. This model is a promising efficient tool for filtering out potentially problematic compounds.

Orthogonal PLS (OPLS) Modeling for Improved Analysis and Interpretation in Drug Design

Partial least squares (PLS) regression is a flexible data analytical approach, which can be made even more versatile and useful by various modifications. In this article we describe the extension into orthogonal PLS modeling, in terms of two new methods, called OPLS and O2PLS, with similar prediction capacity but improved model interpretation.

Modeling chemical interaction profiles: II. Molecular docking, spectral data-activity relationship, and structure-activity relationship models for potent and weak inhibitors of cytochrome P450 CYP3A4 isozyme

Polypharmacy increasingly has become a topic of public health concern, particularly as the U.S. population ages. Drug labels often contain insufficient information to enable the clinician to safely use multiple drugs. Because many of the drugs are bio-transformed by cytochrome P450 (CYP) enzymes, inhibition of CYP activity has long been associated with potentially adverse health effects. In an attempt to reduce the uncertainty pertaining to CYP-mediated drug-drug/chemical interactions, an interagency collaborative group developed a consensus approach to prioritizing information concerning CYP inhibition. The consensus involved computational molecular docking, spectral data-activity relationship (SDAR), and structure-activity relationship (SAR) models that addressed the clinical potency of CYP inhibition. The models were built upon chemicals that were categorized as either potent or weak inhibitors of the CYP3A4 isozyme. The categorization was carried out using information from clinical trials because currently available in vitro high-throughput screening data were not fully representative of the in vivo potency of inhibition. During categorization it was found that compounds, which break the Lipinski rule of five by molecular weight, were about twice more likely to be inhibitors of CYP3A4 compared to those, which obey the rule. Similarly, among inhibitors that break the rule, potent inhibitors were 2–3 times more frequent. The molecular docking classification relied on logistic regression, by which the docking scores from different docking algorithms, CYP3A4 three-dimensional structures, and binding sites on them were combined in a unified probabilistic model. The SDAR models employed a multiple linear regression approach applied to binned 1D ¹³C-NMR and 1D ¹⁵N-NMR spectral descriptors. Structure-based and physical-chemical descriptors were used as the basis for developing SAR models by the decision forest method. Thirty-three potent inhibitors and 88 weak inhibitors of CYP3A4 were used to train the models. Using these models, a synthetic majority rules consensus classifier was implemented, while the confidence of estimation was assigned following the percent agreement strategy. The classifier was applied to a testing set of 120 inhibitors not included in the development of the models. Five compounds of the test set, including known strong inhibitors dalfopristin and tioconazole, were classified as probable potent inhibitors of CYP3A4. Other known strong inhibitors, such as lopinavir, oltipraz, quercetin, raloxifene, and troglitazone, were among 18 compounds classified as plausible potent inhibitors of CYP3A4. The consensus estimation of inhibition potency is expected to aid in the nomination of pharmaceuticals, dietary supplements, environmental pollutants, and occupational and other chemicals for in-depth evaluation of the CYP3A4 inhibitory activity. It may serve also as an estimate of chemical interactions via CYP3A4 metabolic pharmacokinetic pathways occurring through polypharmacy and nutritional and environmental exposures to chemical mixtures.

Classification of anti-HIV compounds using counterpropagation artificial neural networks and decision trees

The main aim of the present work was to collect and categorize anti-HIV molecules in order to identify general structure-activity relationships. In this respect, a total of 5580 drugs and drug-like molecules was collected from 256 different articles published between 1992 and 2010. An algorithm called genetic algorithm-pattern search counterpropagation artificial neural networks (GPS-CPANN) was proposed for the classification of compounds. In addition, the CART (classification and regression trees) method was used for construction of decision trees and finding the best molecular descriptors. The results revealed that the developed CPANN models and decision tree can correctly classify the molecules according to their inhibition mechanisms and activities. Some general parameters such as molecular weight, average molecular weight, number of hydrogen atoms and number of hydroxyl groups were found to be important for describing the inhibition behaviour of anti-HIV agents. The developed classifier models in this work can be used to screen large libraries of compounds to identify those likely to display activity as anti-HIV agents.

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

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 IC₅₀ 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.

QSAR of cytochrome inhibitors

Cytochrome P450 (CYP450) enzymes are predominantly involved in the Phase I metabolism of xenobiotics. Metabolic inhibition and induction can give rise to clinically important drug-drug interactions. Metabolic stability is a prerequisite for sustaining the therapeutically relevant concentrations, and very often drug candidates are sacrificed due to poor metabolic profiles. Computational tools such as quantitative structure-activity relationships are widely used to study different metabolic end points successfully to accelerate the drug discovery process. There are a lot of computational studies on clinically important CYPs already reported in recent years. But other clinically significant families are to yet be explored computationally. Powerfulness of quantitative structure-activity relationship will drive computational chemists to develop new potent and selective inhibitors of different classes of CYPs for the treatment of different diseases with least drug-drug interactions. Furthermore, there is a need to enhance the accuracy, interpretability and confidence in the computational models in accelerating the drug discovery pathways.

Computational prediction of drug binding and rationalisation of selectivity towards cytochromes P450

BACKGROUND

Early in-vitro consideration of metabolism and inhibition of cytochrome P450 has proven its merits over the last 15 years. Simultaneously, many computational drug-design methods have been developed, and are being applied to study the interactions between drug candidates and cytochrome P450 enzymes (P450s).

OBJECTIVE

This review discusses the recent advances of these methods and the implications that are specific for P450s.

METHODS

Mainly focusing on the prediction of binding affinity and ligand selectivity, we outline the applicability of the different methods to answer specific questions. Special emphasis is put on the different levels of theory that are being used in recent computational descriptions of ligand-P450 interactions.

CONCLUSION

P450s offer an additional challenge for computational methods, considering the ambiguities of the catalytic cycle and the significant flexibility of the active site. Different computational methods display different limitations, which is crucial to take into account when choosing the method appropriate to each application.

Classification of highly unbalanced CYP450 data of drugs using cost sensitive machine learning techniques

In this paper, we study the classifications of unbalanced data sets of drugs. As an example we chose a data set of 2D6 inhibitors of cytochrome P450. The human cytochrome P450 2D6 isoform plays a key role in the metabolism of many drugs in the preclinical drug discovery process. We have collected a data set from annotated public data and calculated physicochemical properties with chemoinformatics methods. On top of this data, we have built classifiers based on machine learning methods. Data sets with different class distributions lead to the effect that conventional machine learning methods are biased toward the larger class. To overcome this problem and to obtain sensitive but also accurate classifiers we combine machine learning and feature selection methods with techniques addressing the problem of unbalanced classification, such as oversampling and threshold moving. We have used our own implementation of a support vector machine algorithm as well as the maximum entropy method. Our feature selection is based on the unsupervised McCabe method. The classification results from our test set are compared structurally with compounds from the training set. We show that the applied algorithms enable the effective high throughput in silico classification of potential drug candidates.

In Silico Prediction of Cytochrome P450 2D6 and 3A4 Inhibition Using Gaussian Kernel Weighted k-Nearest Neighbor and Extended Connectivity Fingerprints, Including Structural Fragment Analysis of Inhibitors versus Noninhibitors

Inhibition of cytochrome P450 (CYP) enzymes is unwanted because of the risk of severe side effects due to drug-drug interactions. We present two in silico Gaussian kernel weighted k-nearest neighbor models based on extended connectivity fingerprints that classify CYP2D6 and CYP3A4 inhibition. Data used for modeling consisted of diverse sets of 1153 and 1382 drug candidates tested for CYP2D6 and CYP3A4 inhibition in human liver microsomes. For CYP2D6, 82% of the classified test set compounds were predicted to the correct class. For CYP3A4, 88% of the classified compounds were correctly classified. CYP2D6 and CYP3A4 inhibition were additionally classified for an external test set on 14 drugs, and multidimensional scaling plots showed that the drugs in the external test set were in the periphery of the training sets. Furthermore, fragment analyses were performed and structural fragments frequent in CYP2D6 and CYP3A4 inhibitors and noninhibitors are presented.

QSAR modeling of in vitro inhibition of cytochrome P450 3A4

We report the QSAR modeling of cytochrome P450 3A4 (CYP3A4) enzyme inhibition using four large data sets of in vitro data. These data sets consist of marketed drugs and drug-like compounds all tested in four assays measuring the inhibition of the metabolism of four different substrates by the CYP3A4 enzyme. The four probe substrates are benzyloxycoumarin, testosterone, benzyloxyresorufin, and midazolam. We first show that using state-of-the-art QSAR modeling approaches applied to only one of these four data sets does not lead to predictive models that would be useful for in silico filtering of chemical libraries. We then present the development and the testing of a multiple pharmacophore hypothesis (MPH) that is formulated as a conceptual extension of the traditional QSAR approach to modeling the promiscuous binding of a large variety of drugs to CYP3A4. In the simplest form, the MPH approach takes advantage of the multiple substrate data sets and identifies the binding of test compounds as either proximal or distal relative to that of a given substrate. Application of the approach to the in silico filtering of test compounds for potential inhibitors of CYP3A4 is also presented. In addition to an improvement in the QSAR modeling for the inhibition of CYP3A4, the results from this modeling approach provide structural insights into the drug-enzyme interactions. The existence of multiple inhibition data sets in the BioPrint database motivates the original development of the concept of a multiple pharmacophore hypothesis and provides a unique opportunity for formulating alternative strategies of QSAR modeling of the inhibition of the in vitro metabolism of CYP3A4.

Comparison of methods for chemical-compound affinity prediction

The selection of effective features from various descriptors of chemical compounds and the exploitation of the most appropriate classifier is a momentous issue in improving overall accuracies of virtual screening of chemical compounds. In this article, the performance of various feature-selection methods and various classifiers of chemical compound-protein binding affinities are compared by using six series of compounds: cytochrome P450 2C9 inhibitors, multi-drug-resistance reversal compounds, estrogen receptor ligands, inhibitors of human ether-a-go-go-related genes, and ligands of serotonin receptor 5HT1A and 5HT2A. As a result, it was found that the genetic algorithm was superior to the other feature-selection methods, and its combination with Random Forests and Adaboosts or Baggings gave almost the same performance as support-vector machines and was superior to the other classifiers. The precision and recall of these methods were almost the same or ascendant to those of previous work. The automatically selected descriptors for each protein-compound affinity prediction were plausible and would be informative to interpret the resulting model.