Generalized Proteochemometric Model of Multiple Cytochrome P450 Enzymes and Their Inhibitors

A big data approach with artificial neural network and molecular similarity for chemical data mining and endocrine disruption prediction

CONTEXT:

Chemical toxicity prediction at early stage drug discovery phase has been researched for years, and newest methods are always investigated. Research data comprising chemical physicochemical properties, toxicity, assay, and activity details create massive data which are becoming difficult to manage. Identifying the desired featured chemical with the desired biological activity from millions of chemicals is a challenging task.

AIMS:

In this study, we investigate and explore big data technologies and machine learning approaches to do an efficient chemical data mining for endocrine receptor disruption prediction and virtual compound screening. The power of artificial neural network (ANN) in predicting chemicals' activity toward androgen receptor (AR) and estrogen receptor (ER) and thereby classifying into human endocrine disruptor or nondisruptor is investigated.

SUBJECTS AND METHODS:

Molecules are collected along with their Inhibitory Concentration (IC ₅₀) values toward AR and ER. Training and test datasets are created with active and inactive classes of molecules. Molecular fingerprints of Electro Topological State (E-State) are generated for describing every compound. ANN machine learning model is created using Apache Spark and implemented in Hadoop big data environment. Test chemical's structural similarity toward active class of training compounds is estimated and combined with ANN model for improving prediction accuracy.

RESULTS:

AR and ER predictive models applied on corresponding test datasets gave 86.31% and 89.57% accuracies, respectively, in correctly classifying molecules as disruptor or nondisruptor. Molecular fragments and functional groups are ranked based on their importance in forming ANN model and influence toward the AR and ER disruption behavior. Training molecules that are specific to the test molecules' endocrine disruption prediction are retrieved based on the structural similarity values.

CONCLUSIONS:

The current study demonstrates a new approach of chemical endocrine receptor disruption prediction combining ANN machine learning method and molecular similarity in a big data environment. This method of predictive modeling can be further tested with more receptors and hormones and predictive power can be examined.

Computational Analysis of LOX1 Inhibition Identifies Descriptors Responsible for Binding Selectivity

Lipoxygenases are a family of cytosolic, peripheral membrane enzymes, which catalyze the hydroperoxidation of polyunsaturated fatty acids and are implicated in the pathogenesis of major human diseases. Over the years, a substantial number of scientific reports have introduced inhibitors active against one or another subtype of the enzyme, but the selectivity issue has proved to be a major challenge for drug design. In the present work, we assembled a dataset of 317 structurally diverse molecules hitherto reported as active against 15S-LOX1, 12S-LOX1, and 15S-LOX2 and identified, using supervised machine learning, a set of structural descriptors responsible for the binding selectivity toward the enzyme 15S-LOX1. We subsequently incorporated these descriptors in the training of QSAR models for LOX1 activity and selectivity. The best performing classifiers are two stacked models that include an ensemble of support vector machine, random forest, and k-nearest neighbor algorithms. These models not only can predict LOX1 activity/inactivity but also can discriminate with high accuracy between molecules that exhibit selective activity toward either one of the isozymes 15S-LOX1 and 12S-LOX1.

Exploring the origin of phosphodiesterase inhibition via proteochemometric modeling

A proteochemometric study of a set of phosphodiesterase 4B and 4D inhibitors sheds light on the origin of their inhibition and selectivities.

Proteochemometric Modeling of the Interaction Space of Carbonic Anhydrase and its Inhibitors: An Assessment of Structure-based and Sequence-based Descriptors

Due to its physiological and clinical roles, carbonic anhydrase (CA) is one of the most interesting case studies. There are different classes of CAinhibitors including sulfonamides, polyamines, coumarins and dithiocarbamates (DTCs). However, many of them hardly act as a selective inhibitor against a specific isoform. Therefore, finding highly selective inhibitors for different isoforms of CA is still an ongoing project. Proteochemometrics modeling (PCM) is able to model the bioactivity of multiple compounds against different isoforms of a protein. Therefore, it would be extremely applicable when investigating the selectivity of different ligands towards different receptors. Given the facts, we applied PCM to investigate the interaction space and structural properties that lead to the selective inhibition of CA isoforms by some dithiocarbamates. Our models have provided interesting structural information that can be considered to design compounds capable of inhibiting different isoforms of CA in an improved selective manner. Validity and predictivity of the models were confirmed by both internal and external validation methods; while Y-scrambling approach was applied to assess the robustness of the models. To prove the reliability and the applicability of our findings, we showed how ligands-receptors selectivity can be affected by removing any of these critical findings from the modeling process.

The recent progress in proteochemometric modelling: focusing on target descriptors, cross-term descriptors and application scope

As an extension of the conventional quantitative structure activity relationship models, proteochemometric (PCM) modelling is a computational method that can predict the bioactivity relations between multiple ligands and multiple targets. Traditional PCM modelling includes three essential elements: descriptors (including target descriptors, ligand descriptors and cross-term descriptors), bioactivity data and appropriate learning functions that link the descriptors to the bioactivity data. Since its appearance, PCM modelling has developed rapidly over the past decade by taking advantage of the progress of different descriptors and machine learning techniques, along with the increasing amounts of available bioactivity data. Specifically, the new emerging target descriptors and cross-term descriptors not only significantly increased the performance of PCM modelling but also expanded its application scope from traditional protein-ligand interaction to more abundant interactions, including protein-peptide, protein-DNA and even protein-protein interactions. In this review, target descriptors and cross-term descriptors, as well as the corresponding application scope, are intensively summarized. Additionally, we look forward to seeing PCM modelling extend into new application scopes, such as Target-Catalyst-Ligand systems, with the further development of descriptors, machine learning techniques and increasing amounts of available bioactivity data.

Computational chemogenomics: is it more than inductive transfer?

High-throughput assays challenge us to extract knowledge from multi-ligand, multi-target activity data. In QSAR, weights are statically fitted to each ligand descriptor with respect to a single endpoint or target. However, computational chemogenomics (CG) has demonstrated benefits of learning from entire grids of data at once, rather than building target-specific QSARs. A possible reason for this is the emergence of inductive knowledge transfer (IT) between targets, providing statistical robustness to the model, with no assumption about the structure of the targets. Relevant protein descriptors in CG should allow one to learn how to dynamically adjust ligand attribute weights with respect to protein structure. Hence, models built through explicit learning (EL) by including protein information, while benefitting from IT enhancement, should provide additional predictive capability, notably for protein deorphanization. This interplay between IT and EL in CG modeling is not sufficiently studied. While IT is likely to occur irrespective of the injected target information, it is not clear whether and when boosting due to EL may occur. EL is only possible if protein description is appropriate to the target set under investigation. The key issue here is the search for evidence of genuine EL exceeding expectations based on pure IT. We explore the problem in the context of Support Vector Regression, using more than 9,400 pKi values of 31 GPCRs, where compound-protein interactions are represented by the concatenation of vectorial descriptions of compounds and proteins. This provides a unified framework to generate both IT-enhanced and potentially EL-enabled models, where the difference is toggled by supplied protein information. For EL-enabled models, protein information includes genuine protein descriptors such as typical sequence-based terms, but also the experimentally determined affinity cross-correlation fingerprints. These latter benchmark the expected behavior of a quasi-ideal descriptor capturing the actual functional protein-protein relatedness, and therefore thought to be the most likely to enable EL. EL- and IT-based methods were benchmarked alongside classical QSAR, with respect to cross-validation and deorphanization challenges. A rational method for projecting benchmarked methodologies into a strategy space is given, in the aims that the projection will provide directions for the types of molecule designs possible using a given methodology. While EL-enabled strategies outperform classical QSARs and favorably compare to similar published results, they are, in all respects evaluated herein, not strongly distinguished from IT-enhanced models. Moreover, EL-enabled strategies failed to prove superior in deorphanization challenges. Therefore, this paper raises caution that, contrary to common belief and intuitive expectation, the benefits of chemogenomics models over classical QSAR are quite possibly due less to the injection of protein-related information, and rather impacted more by the effect of inductive transfer, due to simultaneous learning from all of the modeled endpoints. These results show that the field of protein descriptor research needs further improvements to truly realize the expected benefit of EL.

Visually interpretable models of kinase selectivity related features derived from field-based proteochemometrics

Achieving selectivity for small organic molecules toward biological targets is a main focus of drug discovery but has been proven difficult, for example, for kinases because of the high similarity of their ATP binding pockets. To support the design of more selective inhibitors with fewer side effects or with altered target profiles for improved efficacy, we developed a method combining ligand- and receptor-based information. Conventional QSAR models enable one to study the interactions of multiple ligands toward a single protein target, but in order to understand the interactions between multiple ligands and multiple proteins, we have used proteochemometrics, a multivariate statistics method that aims to combine and correlate both ligand and protein descriptions with affinity to receptors. The superimposed binding sites of 50 unique kinases were described by molecular interaction fields derived from knowledge-based potentials and Schrödinger's WaterMap software. Eighty ligands were described by Mold(2), Open Babel, and Volsurf descriptors. Partial least-squares regression including cross-terms, which describe the selectivity, was used for model building. This combination of methods allows interpretation and easy visualization of the models within the context of ligand binding pockets, which can be translated readily into the design of novel inhibitors.

A unified proteochemometric model for prediction of inhibition of cytochrome p450 isoforms

A unified proteochemometric (PCM) model for the prediction of the ability of drug-like chemicals to inhibit five major drug metabolizing CYP isoforms (i.e. CYP1A2, CYP2C9, CYP2C19, CYP2D6 and CYP3A4) was created and made publicly available under the Bioclipse Decision Support open source system at www.cyp450model.org. In regards to the proteochemometric modeling we represented the chemical compounds by molecular signature descriptors and the CYP-isoforms by alignment-independent description of composition and transition of amino acid properties of their protein primary sequences. The entire training dataset contained 63 391 interactions and the best PCM model was obtained using signature descriptors of height 1, 2 and 3 and inducing the model with a support vector machine. The model showed excellent predictive ability with internal AUC = 0.923 and an external AUC = 0.940, as evaluated on a large external dataset. The advantage of PCM models is their extensibility making it possible to extend our model for new CYP isoforms and polymorphic CYP forms. A key benefit of PCM is that all proteins are confined in one single model, which makes it generally more stable and predictive as compared with single target models. The inclusion of the model in Bioclipse Decision Support makes it possible to make virtual instantaneous predictions (∼100 ms per prediction) while interactively drawing or modifying chemical structures in the Bioclipse chemical structure editor.

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.

Identification of cytochrome P450 2D6 and 2C9 substrates and inhibitors by QSAR analysis

This paper presents four new QSAR models for CYP2C9 and CYP2D6 substrate recognition and inhibitor identification based on human clinical data. The models were used to screen a large data set of environmental chemicals for CYP activity, and to analyze the frequency of CYP activity among these compounds. A large fraction of these chemicals were found to be CYP active, and thus potentially capable of affecting human physiology. 20% of the compounds within applicability domain of the models were predicted to be CYP2C9 substrates, and 17% to be inhibitors. The corresponding numbers for CYP2D6 were 9% and 21%. Where the majority of CYP2C9 active compounds were predicted to be both a substrate and an inhibitor at the same time, the CYP2D6 active compounds were primarily predicted to be only inhibitors. It was demonstrated that the models could identify compound classes with a high occurrence of specific CYP activity. An overrepresentation was seen for poly-aromatic hydrocarbons (group of procarcinogens) among CYP2C9 active and mutagenic compounds compared to CYP2C9 inactive and mutagenic compounds. The mutagenicity was predicted with a QSAR model based on Ames in vitro test data.

Toward in silico structure-based ADMET prediction in drug discovery

Quantitative structure-activity relationship (QSAR) methods and related approaches have been used to investigate the molecular features that influence the absorption, distribution, metabolism, excretion and toxicity (ADMET) of drugs. As the three-dimensional structures of several major ADMET proteins become available, structure-based (docking-scoring) computations can be carried out to complement or to go beyond QSAR studies. Applying docking-scoring methods to ADMET proteins is a challenging process because they usually have a large and flexible binding cavity; however, promising results relating to metabolizing enzymes have been reported. After reviewing current trends in the field we applied structure-based methods in the context of receptor flexibility in a case study involving the phase II metabolizing sulfotransferases. Overall, the explored concepts and results suggested that structure-based ADMET profiling will probably join the mainstream during the coming years.

RS-predictor: a new tool for predicting sites of cytochrome P450-mediated metabolism applied to CYP 3A4

This article describes RegioSelectivity-Predictor (RS-Predictor), a new in silico method for generating predictive models of P450-mediated metabolism for drug-like compounds. Within this method, potential sites of metabolism (SOMs) are represented as "metabolophores": A concept that describes the hierarchical combination of topological and quantum chemical descriptors needed to represent the reactivity of potential metabolic reaction sites. RS-Predictor modeling involves the use of metabolophore descriptors together with multiple-instance ranking (MIRank) to generate an optimized descriptor weight vector that encodes regioselectivity trends across all cases in a training set. The resulting pathway-independent (O-dealkylation vs N-oxidation vs Csp(3) hydroxylation, etc.), isozyme-specific regioselectivity model may be used to predict potential metabolic liabilities. In the present work, cross-validated RS-Predictor models were generated for a set of 394 substrates of CYP 3A4 as a proof-of-principle for the method. Rank aggregation was then employed to merge independently generated predictions for each substrate into a single consensus prediction. The resulting consensus RS-Predictor models were shown to reliably identify at least one observed site of metabolism in the top two rank-positions on 78% of the substrates. Comparisons between RS-Predictor and previously described regioselectivity prediction methods reveal new insights into how in silico metabolite prediction methods should be compared.

Services for prediction of drug susceptibility for HIV proteases and reverse transcriptases at the HIV drug research centre

SUMMARY

The HIV Drug Research Centre (HIVDRC) has established Web services for prediction of drug susceptibility for HIV proteases and reverse transcriptases. The services are based on two proteochemometric models which accepts a protease or reverse transcriptase sequence in amino acid form, and outputs the predicted drug susceptibility values. The predictions are based on a comprehensive analysis where all the relevant inhibitors are included, resulting in models with excellent predictive capabilities.

AVAILABILITY AND IMPLEMENTATION

The services are implemented as interoperable Web services (REST and XMPP), with supporting web pages to allow for individual analyses. A set of plugins were also developed which make the services available from the Bioclipse workbench for life science. Services are available at http://www.hivdrc.org/services.

Proteochemometric modeling as a tool to design selective compounds and for extrapolating to novel targets

Proteochemometric modeling is founded on the principles of QSAR but is able to benefit from additional information in model training due to the inclusion of target information.

Classification of cytochrome P450 inhibitors with respect to binding free energy and pIC50 using common molecular descriptors

Virtual screening of chemical libraries following experimental assays of drug candidates is a common procedure in structure based drug discovery. However, the relationship between binding free energies and biological activities (pIC50) of drug candidates is still an unsolved issue that limits the efficiency and speed of drug development processes. In this study, the relationship between them is investigated based on a common molecular descriptor set for human cytochrome P450 enzymes (CYPs). CYPs play an important role in drug-drug interactions, drug metabolism, and toxicity. Therefore, in silico prediction of CYP inhibition by drug candidates is one of the major considerations in drug discovery. The combination of partial least-squares regression (PLSR) and a variety of classification algorithms were employed by considering this relationship as a classification problem. Our results indicate that PLSR with classification is a powerful tool to predict more than one output such as binding free energy and pIC50 simultaneously. PLSR with mixed-integer linear programming based hyperboxes predicts the binding free energy and pIC50 with a mean accuracy of 87.18% (min: 81.67% max: 97.05%) and 88.09% (min: 79.83% max: 92.90%), respectively, for the cytochrome p450 superfamily using the common 6 molecular descriptors with a 10-fold cross-validation.

Organizing bioactive compound discovery in target families

The sequencing of genomes gave access to the complete set of building blocks for organisms of various species. A plethora of "-omics"-technologies has been developed to investigate the dynamic interactions of the building blocks in order to understand the functioning of living organisms. This has given rise to the clustering of proteins into target families based on the phylogenetic and structural similarities. In this chapter we will discuss how the concept of target families enables to investigate and modulate biochemical function in the quest to chart Chemical and Biological Spaces.