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

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.

Recent progresses in the exploration of machine learning methods as in-silico ADME prediction tools

In-silico methods have been explored as potential tools for assessing ADME and ADME regulatory properties particularly in early drug discovery stages. Machine learning methods, with their ability in classifying diverse structures and complex mechanisms, are well suited for predicting ADME and ADME regulatory properties. Recent efforts have been directed at the broadening of application scopes and the improvement of predictive performance with particular focuses on the coverage of ADME properties, and exploration of more diversified training data, appropriate molecular features, and consensus modeling. Moreover, several online machine learning ADME prediction servers have emerged. Here we review these progresses and discuss the performances, application prospects and challenges of exploring machine learning methods as useful tools in predicting ADME and ADME regulatory properties.

Mass spectrometry based identification of geometric isomers during metabolic stability study of a new cytotoxic sulfonamide derivatives supported by quantitative structure-retention relationships

A set of 15 new sulphonamide derivatives, presenting antitumor activity have been subjected to a metabolic stability study. The results showed that besides products of biotransformation, some additional peaks occurred in chromatograms. Tandem mass spectrometry revealed the same mass and fragmentation pathway, suggesting that geometric isomerization occurred. Thus, to support this hypothesis, quantitative structure-retention relationships were applied. Human liver microsomes were used as an in vitro model of metabolism. The biotransformation reactions were tracked by liquid chromatography assay and additionally, fragmentation mass spectra were recorded. In silico molecular modeling at a semi-empirical level was conducted as a starting point for molecular descriptor calculations. A quantitative structure-retention relationship model was built applying multiple linear regression based on selected three-dimensional descriptors. The studied compounds revealed high metabolic stability, with a tendency to form hydroxylated biotransformation products. However, significant chemical instability in conditions simulating human body fluids was noticed. According to literature and MS data geometrical isomerization was suggested. The developed in sillico model was able to describe the relationship between the geometry of isomer pairs and their chromatographic retention properties, thus it supported the hypothesis that the observed pairs of peaks are most likely geometric isomers. However, extensive structural investigations are needed to fully identify isomers' geometry. An effort to describe MS fragmentation pathways of novel chemical structures is often not enough to propose structures of potent metabolites and products of other chemical reactions that can be observed in compound solutions at early drug discovery studies. The results indicate that the relatively non-expensive and not time- and labor-consuming in sillico approach could be a good supportive tool assisting the identification of cis-trans isomers based on retention data. This methodology can be helpful during the structural identification of biotransformation and degradation products of new chemical entities--potential new drugs.

Evaluation of the HC-04 cell line as an in vitro model for mechanistic assessment of changes in hepatic cytochrome P450 3A during adenovirus infection

HC-04 cells were evaluated as an in vitro model for mechanistic study of changes in the function of hepatic CYP3A during virus infection. Similar to in vivo observations, infection with a first generation recombinant adenovirus significantly inhibited CYP3A4 catalytic activity in an isoform-specific manner. Virus (MOI 100) significantly reduced expression of the retinoid X receptor (RXR) by 30% 96 hours after infection. Cytoplasmic concentrations of the pregnane X receptor (PXR) were reduced by 50%, whereas the amount of the constitutive androstane receptor (CAR) in the nuclear fraction doubled with respect to uninfected controls. Hepatocyte nuclear factor 4α (HNF-4α) and peroxisome proliferator-activated receptor γ coactivator 1α (PGC-1α) were also reduced by ∼70% during infection. Virus suppressed CYP3A4 activity in the presence of the PXR agonist rifampicin and did not affect CYP3A4 activity in the presence of the CAR agonist CITCO [6-(4-chlorophenyl) imidazo[2,1-b][1,3]thiazole-5-carbaldehyde-O-(3,4-dichlorobenzyl)oxime], suggesting that virus-induced modification of PXR may be responsible for observed changes in hepatic CYP3A4. The HC-04 cell line is easy to maintain, and CYP3A4 in these cells was responsive to known inducers and suppressors. Dexamethasone (200 μM) and phenobarbital (500 μM) increased activity by 230 and 124%, whereas ketoconazole (10 μM) and lipopolysaccharide (LPS) (10 μg/ml) reduced activity by 90 and 92%, respectively. This suggests that HC-04 cells can be a valuable tool for mechanistic study of drug metabolism during infection and for routine toxicological screening of novel compounds prior to use in the clinic.

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.

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.

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.

Structure based discovery of small molecules to regulate the activity of human insulin degrading enzyme

BACKGROUND

Insulin-degrading enzyme (IDE) is an allosteric Zn(+2) metalloprotease involved in the degradation of many peptides including amyloid-β, and insulin that play key roles in Alzheimer's disease (AD) and type 2 diabetes mellitus (T2DM), respectively. Therefore, the use of therapeutic agents that regulate the activity of IDE would be a viable approach towards generating pharmaceutical treatments for these diseases. Crystal structure of IDE revealed that N-terminal has an exosite which is ∼30 Å away from the catalytic region and serves as a regulation site by orientation of the substrates of IDE to the catalytic site. It is possible to find small molecules that bind to the exosite of IDE and enhance its proteolytic activity towards different substrates.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we applied structure based drug design method combined with experimental methods to discover four novel molecules that enhance the activity of human IDE. The novel compounds, designated as D3, D4, D6, and D10 enhanced IDE mediated proteolysis of substrate V, insulin and amyloid-β, while enhanced degradation profiles were obtained towards substrate V and insulin in the presence of D10 only.

CONCLUSION/SIGNIFICANCE

This paper describes the first examples of a computer-aided discovery of IDE regulators, showing that in vitro and in vivo activation of this important enzyme with small molecules is possible.

Optimization based tumor classification from microarray gene expression data

BACKGROUND

An important use of data obtained from microarray measurements is the classification of tumor types with respect to genes that are either up or down regulated in specific cancer types. A number of algorithms have been proposed to obtain such classifications. These algorithms usually require parameter optimization to obtain accurate results depending on the type of data. Additionally, it is highly critical to find an optimal set of markers among those up or down regulated genes that can be clinically utilized to build assays for the diagnosis or to follow progression of specific cancer types. In this paper, we employ a mixed integer programming based classification algorithm named hyper-box enclosure method (HBE) for the classification of some cancer types with a minimal set of predictor genes. This optimization based method which is a user friendly and efficient classifier may allow the clinicians to diagnose and follow progression of certain cancer types.

METHODOLOGY/PRINCIPAL FINDINGS

We apply HBE algorithm to some well known data sets such as leukemia, prostate cancer, diffuse large B-cell lymphoma (DLBCL), small round blue cell tumors (SRBCT) to find some predictor genes that can be utilized for diagnosis and prognosis in a robust manner with a high accuracy. Our approach does not require any modification or parameter optimization for each data set. Additionally, information gain attribute evaluator, relief attribute evaluator and correlation-based feature selection methods are employed for the gene selection. The results are compared with those from other studies and biological roles of selected genes in corresponding cancer type are described.

CONCLUSIONS/SIGNIFICANCE

The performance of our algorithm overall was better than the other algorithms reported in the literature and classifiers found in WEKA data-mining package. Since it does not require a parameter optimization and it performs consistently very high prediction rate on different type of data sets, HBE method is an effective and consistent tool for cancer type prediction with a small number of gene markers.