Development of a New Predictive Model for Interactions with Human Cytochrome P450 2A6 Using Pharmacophore Ensemble/Support Vector Machine (PhE/SVM) Approach

Insight Analysis of Promiscuous Estrogen Receptor α-Ligand Binding by a Novel Machine Learning Scheme

Estrogen receptor α (ERα) plays a significant role in occurrence of breast cancer and may cause various adverse side-effects when ERα is an off-target protein. A theoretical model was derived to predict the binding affinity of ERα using the pharmacophore ensemble/support vector machine (PhE/SVM) scheme to consider the promiscuous characteristic of ERα. The estimations by PhE/SVM were discovered to be in good agreement with the observed values for those training molecules ( n = 31, r² = 0.80, q_CV² = 0.77, RMSE = 0.57, s = 0.58), test molecules ( n = 179, q² = 0.91-0.96, RMSE = 0.33, s = 0.26) and outliers ( n = 15, q² = 0.80-0.86, RMSE = 0.56, s = 0.49). When subjected to various statistical validations, the PhE/SVM model consistently fulfilled the strictest criteria. A mock test also asserted its predictivity. When compared with crystal structures, the calculated results are consistent with the reported ERα-ligand co-complex structure, and the plasticity nature of ERα is also disclosed. Consequently, this precise, fast, and robust model can be adopted to predict ERα-ligand binding affinities and to design safer non-ERα-targeted pharmaceuticals in the process of drug discovery and development.

Discovery of a selective, safe and novel anti-malarial compound with activity against chloroquine resistant strain of Plasmodium falciparum

In recent years the DNA minor groove has attracted much attention for the development of anti-malarial agents. In view of this we have attempted to discover novel DNA minor groove binders through in-silico and in-vitro workflow. A rigorously validated pharmacophore model comprising of two positive ionizable (PI), one hydrophobic (HY) and one ring aromatic (RA) features was used to mine NCI chemical compound database. This led to retrieval of many hits which were screened on the basis of estimated activity, fit value and Lipinski's violation. Finally two compounds NSC639017 and NSC371488 were evaluated for their in-vitro anti-malarial activities against Plasmodium falciparum 3D7 (CQ sensitive) and K1 (CQ resistant) strains by SYBR green-I based fluorescence assay. The results revealed that out of two, NSC639017 posses excellent anti-malarial activity particularly against chloroquine resistant strain and moreover NSC639017 also appeared to be safe (CC50 126.04 μg/ml) and selective during cytotoxicity evaluation.

In silico methods for predicting drug-drug interactions with cytochrome P-450s, transporters and beyond

Drug-drug interactions (DDIs) are associated with severe adverse effects that may lead to the patient requiring alternative therapeutics and could ultimately lead to drug withdrawal from the market if they are severe. To prevent the occurrence of DDI in the clinic, experimental systems to evaluate drug interaction have been integrated into the various stages of the drug discovery and development process. A large body of knowledge about DDI has also accumulated through these studies and pharmacovigillence systems. Much of this work to date has focused on the drug metabolizing enzymes such as cytochrome P-450s as well as drug transporters, ion channels and occasionally other proteins. This combined knowledge provides a foundation for a hypothesis-driven in silico approach, using either cheminformatics or physiologically based pharmacokinetics (PK) modeling methods to assess DDI potential. Here we review recent advances in these approaches with emphasis on hypothesis-driven mechanistic models for important protein targets involved in PK-based DDI. Recent efforts with other informatics approaches to detect DDI are highlighted. Besides DDI, we also briefly introduce drug interactions with other substances, such as Traditional Chinese Medicines to illustrate how in silico modeling can be useful in this domain. We also summarize valuable data sources and web-based tools that are available for DDI prediction. We finally explore the challenges we see faced by in silico approaches for predicting DDI and propose future directions to make these computational models more reliable, accurate, and publically accessible.

Directed evolution of cytochrome P450 enzymes for biocatalysis: exploiting the catalytic versatility of enzymes with relaxed substrate specificity

Cytochrome P450 enzymes are renowned for their ability to insert oxygen into an enormous variety of compounds with a high degree of chemo- and regio-selectivity under mild conditions. This property has been exploited in Nature for an enormous variety of physiological functions, and representatives of this ancient enzyme family have been identified in all kingdoms of life. The catalytic versatility of P450s makes them well suited for repurposing for the synthesis of fine chemicals such as drugs. Although these enzymes have not evolved in Nature to perform the reactions required for modern chemical industries, many P450s show relaxed substrate specificity and exhibit some degree of activity towards non-natural substrates of relevance to applications such as drug development. Directed evolution and other protein engineering methods can be used to improve upon this low level of activity and convert these promiscuous generalist enzymes into specialists capable of mediating reactions of interest with exquisite regio- and stereo-selectivity. Although there are some notable successes in exploiting P450s from natural sources in metabolic engineering, and P450s have been proven repeatedly to be excellent material for engineering, there are few examples to date of practical application of engineered P450s. The purpose of the present review is to illustrate the progress that has been made in altering properties of P450s such as substrate range, cofactor preference and stability, and outline some of the remaining challenges that must be overcome for industrial application of these powerful biocatalysts.

In silico prediction of inhibition of promiscuous breast cancer resistance protein (BCRP/ABCG2)

BACKGROUND

Breast cancer resistant protein has an essential role in active transport of endogenous substances and xenobiotics across extracellular and intracellular membranes along with P-glycoprotein. It also plays a major role in multiple drug resistance and permeation of blood-brain barrier. Therefore, it is of great importance to derive theoretical models to predict the inhibition of both transporters in the process of drug discovery and development. Hitherto, very limited BCRP inhibition predictive models have been proposed as compared with its P-gp counterpart.

METHODOLOGY/PRINCIPAL FINDINGS

An in silico BCRP inhibition model was developed in this study using the pharmacophore ensemble/support vector machine scheme to take into account the promiscuous nature of BCRP. The predictions by the PhE/SVM model were found to be in good agreement with the observed values for those molecules in the training set (n= 22, r2 =0.82, qCV2=0.73, RMSE= 0.40, s = 0.24), test set (n =97, q2=0.75-0.89, RMSE= 0.31, s= 0.21), and outlier set (n= 16, q2 =0.72-0.91, RMSE= 0.29, s=0.17). When subjected to a variety of statistical validations, the developed PhE/SVM model consistently met the most stringent criteria. A mock test by HIV protease inhibitors also asserted its predictivity.

CONCLUSIONS/SIGNIFICANCE

It was found that this accurate, fast, and robust PhE/SVM model can be employed to predict the BCRP inhibition of structurally diverse molecules that otherwise cannot be carried out by any other methods in a high-throughput fashion to design therapeutic agents with insignificant drug toxicity and unfavorable drug-drug interactions mediated by BCRP to enhance clinical efficacy and/or circumvent drug resistance.

Classification models for safe drug molecules

Frequent failure of drug candidates during development stages remains the major deterrent for an early introduction of new drug molecules. The drug toxicity is the major cause of expensive late-stage development failures. An early identification/optimization of the most favorable molecule will naturally save considerable cost, time, human efforts and minimize animal sacrifice. (Quantitative) Structure Activity Relationships [(Q)SARs] represent statistically derived predictive models correlating biological activity (including desirable therapeutic effect and undesirable side effects) of chemicals (drugs/toxicants/environmental pollutants) with molecular descriptors and/or properties. (Q)SAR models which categorize the available data into two or more groups/classes are known as classification models. Numerous techniques of diverse nature are being presently employed for development of classification models. Though there is an increasing use of classification models for prediction of either biological activity or toxicity, the future trend will naturally be towards the development of classification models capable of simultaneous prediction of biological activity, toxicity, and pharmacokinetic parameters so as to accelerate development of bioavailable safe drug molecules.

The promiscuous binding of pharmaceutical drugs and their transporter-mediated uptake into cells: what we (need to) know and how we can do so

A recent paper in this journal sought to counter evidence for the role of transport proteins in effecting drug uptake into cells, and questions that transporters can recognize drug molecules in addition to their endogenous substrates. However, there is abundant evidence that both drugs and proteins are highly promiscuous. Most proteins bind to many drugs and most drugs bind to multiple proteins (on average more than six), including transporters (mutations in these can determine resistance); most drugs are known to recognise at least one transporter. In this response, we alert readers to the relevant evidence that exists or is required. This needs to be acquired in cells that contain the relevant proteins, and we highlight an experimental system for simultaneous genome-wide assessment of carrier-mediated uptake in a eukaryotic cell (yeast).

Protein flexibility in docking and surface mapping

Structure-based drug design has become an essential tool for rapid lead discovery and optimization. As available structural information has increased, researchers have become increasingly aware of the importance of protein flexibility for accurate description of the native state. Typical protein-ligand docking efforts still rely on a single rigid receptor, which is an incomplete representation of potential binding conformations of the protein. These rigid docking efforts typically show the best performance rates between 50 and 75%, while fully flexible docking methods can enhance pose prediction up to 80-95%. This review examines the current toolbox for flexible protein-ligand docking and receptor surface mapping. Present limitations and possibilities for future development are discussed.

Prediction of promiscuous p-glycoprotein inhibition using a novel machine learning scheme

BACKGROUND

P-glycoprotein (P-gp) is an ATP-dependent membrane transporter that plays a pivotal role in eliminating xenobiotics by active extrusion of xenobiotics from the cell. Multidrug resistance (MDR) is highly associated with the over-expression of P-gp by cells, resulting in increased efflux of chemotherapeutical agents and reduction of intracellular drug accumulation. It is of clinical importance to develop a P-gp inhibition predictive model in the process of drug discovery and development.

METHODOLOGY/PRINCIPAL FINDINGS

An in silico model was derived to predict the inhibition of P-gp using the newly invented pharmacophore ensemble/support vector machine (PhE/SVM) scheme based on the data compiled from the literature. The predictions by the PhE/SVM model were found to be in good agreement with the observed values for those structurally diverse molecules in the training set (n = 31, r(2) = 0.89, q(2) = 0.86, RMSE = 0.40, s = 0.28), the test set (n = 88, r(2) = 0.87, RMSE = 0.39, s = 0.25) and the outlier set (n = 11, r(2) = 0.96, RMSE = 0.10, s = 0.05). The generated PhE/SVM model also showed high accuracy when subjected to those validation criteria generally adopted to gauge the predictivity of a theoretical model.

CONCLUSIONS/SIGNIFICANCE

This accurate, fast and robust PhE/SVM model that can take into account the promiscuous nature of P-gp can be applied to predict the P-gp inhibition of structurally diverse compounds that otherwise cannot be done by any other methods in a high-throughput fashion to facilitate drug discovery and development by designing drug candidates with better metabolism profile.

Computational prediction of metabolism: sites, products, SAR, P450 enzyme dynamics, and mechanisms

Metabolism of xenobiotics remains a central challenge for the discovery and development of drugs, cosmetics, nutritional supplements, and agrochemicals. Metabolic transformations are frequently related to the incidence of toxic effects that may result from the emergence of reactive species, the systemic accumulation of metabolites, or by induction of metabolic pathways. Experimental investigation of the metabolism of small organic molecules is particularly resource demanding; hence, computational methods are of considerable interest to complement experimental approaches. This review provides a broad overview of structure- and ligand-based computational methods for the prediction of xenobiotic metabolism. Current computational approaches to address xenobiotic metabolism are discussed from three major perspectives: (i) prediction of sites of metabolism (SOMs), (ii) elucidation of potential metabolites and their chemical structures, and (iii) prediction of direct and indirect effects of xenobiotics on metabolizing enzymes, where the focus is on the cytochrome P450 (CYP) superfamily of enzymes, the cardinal xenobiotics metabolizing enzymes. For each of these domains, a variety of approaches and their applications are systematically reviewed, including expert systems, data mining approaches, quantitative structure–activity relationships (QSARs), and machine learning-based methods, pharmacophore-based algorithms, shape-focused techniques, molecular interaction fields (MIFs), reactivity-focused techniques, protein–ligand docking, molecular dynamics (MD) simulations, and combinations of methods. Predictive metabolism is a developing area, and there is still enormous potential for improvement. However, it is clear that the combination of rapidly increasing amounts of available ligand- and structure-related experimental data (in particular, quantitative data) with novel and diverse simulation and modeling approaches is accelerating the development of effective tools for prediction of in vivo metabolism, which is reflected by the diverse and comprehensive data sources and methods for metabolism prediction reviewed here. This review attempts to survey the range and scope of computational methods applied to metabolism prediction and also to compare and contrast their applicability and performance.

Substrate selectivity of drug-metabolizing cytochrome P450s predicted from crystal structures and in silico modeling

Enormous efforts toward predicting the metabolic fate of a drug have been driven by the high attrition rate in drug development. To accelerate such efforts, it is critical to elucidate the molecular mechanisms of drug recognition by drug-metabolizing enzymes. Therefore, it is not surprising that an increasing number of crystal structures have been determined (by X-ray crystallography) and numerous insightful in silico (computational) models have been established for the most important metabolic enzymes, cytochrome P450s (CYPs). In this review, we provide a detailed analysis of the available crystal structures for CYPs to reveal the structural features and protein flexibility determining substrate selectivity. The ligand-based in silico models (including pharmacophore and molecular field analysis models) are also discussed, with a focus on their ability to characterize the structural features of the substrates for various CYP isoforms.

Identification of novel CYP2A6 inhibitors by virtual screening

The human CYP2A6 enzyme metabolises several xenobiotics including nicotine, the addictive component in tobacco. Reduced activity of CYP2A6, either for genetic reasons or by administering inhibitors of CYP2A6, reduces tobacco smoking. The aim was to design novel inhibitors of CYP2A6 using 3D-QSAR analysis combined with virtual screening. A 3D-QSAR model was utilised to identify the most important features of the inhibitors, and this knowledge was used to design inhibitors for CYP2A6. Chemical database screening yielded several potent inhibitor candidates such as alkylamine derivatives (compound no. 5, IC(50)=0.1 μM) and 1-benzothiophene-3-carbaldehyde that can be used as lead compounds in the development of drugs for smoking reduction therapy.

Pharmaceutical drug transport: the issues and the implications that it is essentially carrier-mediated only

All cells necessarily contain tens, if not hundreds, of carriers for nutrients and intermediary metabolites, and the human genome codes for more than 1000 carriers of various kinds. Here, we illustrate using a typical literature example the widespread but erroneous nature of the assumption that the 'background' or 'passive' permeability to drugs occurs in the absence of carriers. Comparison of the rate of drug transport in natural versus artificial membranes shows discrepancies in absolute magnitudes of 100-fold or more, with the carrier-containing cells showing the greater permeability. Expression profiling data show exactly which carriers are expressed in which tissues. The recognition that drugs necessarily require carriers for uptake into cells provides many opportunities for improving the effectiveness of the drug discovery process.

Machine learning algorithms for the prediction of hERG and CYP450 binding in drug development

IMPORTANCE OF THE FIELD

The cost of developing new drugs is estimated at approximately $1 billion; the withdrawal of a marketed compound due to toxicity can result in serious financial loss for a pharmaceutical company. There has been a greater interest in the development of in silico tools that can identify compounds with metabolic liabilities before they are brought to market.

AREAS COVERED IN THIS REVIEW

The two largest classes of machine learning (ML) models, which will be discussed in this review, have been developed to predict binding to the human ether-a-go-go related gene (hERG) ion channel protein and the various CYP isoforms. Being able to identify potentially toxic compounds before they are made would greatly reduce the number of compound failures and the costs associated with drug development.

WHAT THE READER WILL GAIN

This review summarizes the state of modeling hERG and CYP binding towards this goal since 2003 using ML algorithms.

TAKE HOME MESSAGE

A wide variety of ML algorithms that are comparable in their overall performance are available. These ML methods may be applied regularly in discovery projects to flag compounds with potential metabolic liabilities.

Classification of cytochrome p(450) activities using machine learning methods

The cytochrome P(450) (CYP) system plays an integral part in the metabolism of drugs and other xenobiotics. Knowledge of the structural features required for interaction with any of the different isoforms of the CYP system is therefore immensely valuable in early drug discovery. In this paper, we focus on three major isoforms (CYP 1A2, CYP 2D6, and CYP 3A4) and present a data set of 335 structurally diverse drug compounds classified for their interaction (as substrate, inhibitor, or any interaction) with these isoforms. We also present machine learning models using a variety of commonly used methods (k-nearest neighbors, decision tree induction using the CHAID and CRT algorithms, random forests, artificial neural networks, and support vector machines using the radial basis function (RBF) and homogeneous polynomials as kernel functions). We discuss the physicochemical features relevant for each end point and compare it to similar studies. Many of these models perform exceptionally well, even with 10-fold cross-validation, yielding corrected classification rates of 81.7 to 91.9% for CYP 1A2, 89.2 to 92.9% for CYP 2D6, and 87.4 to 89.9% for CYP3A4. Our models help in understanding the structural requirements for CYP interactions and can serve as sensitive tools in virtual screenings and lead optimization for toxicological profiles in drug discovery.

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.