Exploration of Enzyme−Ligand Interactions in CYP2D6 & 3A4 Homology Models and Crystal Structures Using a Novel Computational Approach

Computational strategies to explore antimalarial thiazine alkaloid lead compounds based on an Australian marine sponge Plakortis Lita

In this work, an attempt was made to propose new leads based on the natural scaffold Thiaplakortone-A active against malaria. The 2D QSAR studies suggested that three descriptors correlate with the anti-malarial activity with an R² value of 0.814. Robustness, reliability, and predictive power of the model were tested by internal validation, external validation, Y-scrambling, and applicability domain analysis. HQSAR studies were carried out as an additional tool to find the sub-structural fingerprints. The CoMFA and CoMSIA models gave Q² values of 0.813 and 0.647, and [Formula: see text] values of 0.994 and 0.984, respectively. Using the 2D-QSAR equation, the activity values of the seven modified compounds were calculated and it was found that three molecules showed good anti-malarial activity. Molecular docking of the 42 Thiaplakortone-A derivatives with Plasmodium falciparum calcium-dependent protein kinase 1 (PfCDPK1) was carried out to find out protein-ligand interactions. Data mining of the bioassay data-set AID: 504850 using the classifier based on Random Forest of Weka suggested that all of the eight molecules selected and three out of the seven virtual molecules were anti-malarial active. Both the virtual molecules and drug molecules were docked with CYP3A4, indicating that the virtual molecules could metabolize easily. Toxicity studies using Osiris shows that three molecules showed no toxic characters.

Integrated structure- and ligand-based in silico approach to predict inhibition of cytochrome P450 2D6

MOTIVATION

Cytochrome P450 (CYP) is a superfamily of enzymes responsible for the metabolism of drugs, xenobiotics and endogenous compounds. CYP2D6 metabolizes about 30% of drugs and predicting potential CYP2D6 inhibition is important in early-stage drug discovery.

RESULTS

We developed an original in silico approach for the prediction of CYP2D6 inhibition combining the knowledge of the protein structure and its dynamic behavior in response to the binding of various ligands and machine learning modeling. This approach includes structural information for CYP2D6 based on the available crystal structures and molecular dynamic simulations (MD) that we performed to take into account conformational changes of the binding site. We performed modeling using three learning algorithms--support vector machine, RandomForest and NaiveBayesian--and we constructed combined models based on topological information of known CYP2D6 inhibitors and predicted binding energies computed by docking on both X-ray and MD protein conformations. In addition, we identified three MD-derived structures that are capable all together to better discriminate inhibitors and non-inhibitors compared with individual CYP2D6 conformations, thus ensuring complementary ligand profiles. Inhibition models based on classical molecular descriptors and predicted binding energies were able to predict CYP2D6 inhibition with an accuracy of 78% on the training set and 75% on the external validation set.

Rationalization of stereospecific binding of propranolol to cytochrome P450 2D6 by free energy calculations

Abstract

Cytochrome P450 2D6 is a major drug-metabolising enzyme with a wide substrate range. A single-point mutation introduced in this enzyme induces stereoselective binding of R and S-propranolol whereas the wild type has no preference. The system has previously been studied both experimentally and computationally (de Graaf et al. in Eur Biophys J 36:589–599, 2007a). The in silico study reported hysteresis and significant deviations from closure of thermodynamic cycles, probably because of lack of sampling. Here, we focus on the effect of prolonged simulation time and enhanced sampling methods, such as Hamiltonian replica exchange, to reduce these problems and to improve the precision of free energy calculations. Finally we rationalize the results at a molecular level and compare data with experimental findings and previously estimated free energies.

Graphical Abstract

Propranolol (PPD) binds stereospecifically to cytochrome P450 2D6, if the mutation F483A is introduced, by reducing binding affinity to R-PPD. This observation has previously been studied by molecular dynamics, by use of several different thermodynamic cycles. Previous computational results suffered from internal inconsistencies, probably because of insufficient sampling. Here we use prolonged simulation time and Hamiltonian replica exchange to overcome the inconsistencies, and give a possible rationalization of the phenomenon at a molecular level.

Development of a computational tool to rival experts in the prediction of sites of metabolism of xenobiotics by p450s

The metabolism of xenobiotics--and more specifically drugs--in the liver is a critical process controlling their half-life. Although there exist experimental methods, which measure the metabolic stability of xenobiotics and identify their metabolites, developing higher throughput predictive methods is an avenue of research. It is expected that predicting the chemical nature of the metabolites would be an asset for designing safer drugs and/or drugs with modulated half-lives. We have developed IMPACTS (In-silico Metabolism Prediction by Activated Cytochromes and Transition States), a computational tool combining docking to metabolic enzymes, transition state modeling, and rule-based substrate reactivity prediction to predict the site of metabolism (SoM) of xenobiotics. Its application to sets of CYP1A2, CYP2C9, CYP2D6, and CYP3A4 substrates and comparison to experts' predictions demonstrates its accuracy and significance. IMPACTS identified an experimentally observed SoM in the top 2 predicted sites for 77% of the substrates, while the accuracy of biotransformation experts' prediction was 65%. Application of IMPACTS to external sets and comparison of its accuracy to those of eleven other methods further validated the method implemented in IMPACTS.

High-throughput screening technologies for drug glucuronidation profiling

A significant number of endogenous and exogenous compounds, including many therapeutic agents, are metabolized in humans via glucuronidation, catalysed by uridine diphosphoglucurono-syltransferases (UGTs). The study of the UGTs is a growing field of research, with constantly accumulated and updated information regarding UGT structure, purification, substrate specificity and inhibition, including clinically relevant drug interactions. Development of reliable UGT assays for the assessment of individual isoform substrate specificity and for the discovery of novel isoform-specific substrates and inhibitors is crucial for understanding the function and regulation of the UGT enzyme family and its clinical and pharmacological relevance. High-throughput screening (HTS) is a powerful technology used to search for novel substrates and inhibitors for a wide variety of targets. However, application of HTS in the context of UGTs is complicated because of the poor stability, low levels of expression, low affinity and broad substrate specificity of the enzymes, combined with difficulties in obtaining individual UGT isoforms in purified format, and insufficient information regarding isoform-specific substrates and inhibitors. This review examines the current status of HTS assays used in the search for novel UGT substrates and inhibitors, emphasizing advancements and challenges in HTS technologies for drug glucuronidation profiling, and discusses possible avenues for future advancement of the field.

New insights into the structural characteristics and functional relevance of the human cytochrome P450 2D6 enzyme

To date, the crystal structures of at least 12 human CYPs (1A2, 2A6, 2A13, 2C8, 2C9, 2D6, 2E1, 2R1, 3A4, 7A1, 8A1, and 46A1) have been determined. CYP2D6 accounts for only a small percentage of all hepatic CYPs (< 2%), but it metabolizes approximately 25% of clinically used drugs with significant polymorphisms. CYP2D6 also metabolizes procarcinogens and neurotoxins, such as 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine, 1,2,3,4-tetrahydroquinoline, and indolealkylamines. Moreover, the enzyme utilizes hydroxytryptamines and neurosteroids as endogenous substrates. Typical CYP2D6 substrates are usually lipophilic bases with an aromatic ring and a nitrogen atom, which can be protonated at physiological pH. Substrate binding is generally followed by oxidation (5-7 A) from the proposed nitrogen-Asp301 interaction. A number of homology models have been constructed to explore the structural features of CYP2D6, while antibody studies also provide useful structural information. Site-directed mutagenesis studies have demonstrated that Glu216, Asp301, Phe120, Phe481, and Phe483 play important roles in determining the binding of ligands to CYP2D6. The structure of human CYP2D6 has been recently determined and shows the characteristic CYP fold observed for other members of the CYP superfamily. The lengths and orientations of the individual secondary structural elements in the CYP2D6 structure are similar to those seen in other human CYP2 members, such as CYP2C9 and 2C8. The 2D6 structure has a well-defined active-site cavity located above the heme group with a volume of approximately 540 A(3), which is larger than equivalent cavities in CYP2A6 (260 A(3)), 1A2 (375 A(3)), and 2E1 (190 A(3)), but smaller than those in CYP3A4 (1385 A(3)) and 2C8 (1438 A(3)). Further studies are required to delineate the molecular mechanisms involved in CYP2D6 ligand interactions and their implications for drug development and clinical practice.

SHOP: receptor-based scaffold HOPping by GRID-based similarity searches

A new field-derived 3D method for receptor-based scaffold hopping, implemented in the software SHOP, is presented. Information from a protein-ligand complex is utilized to substitute a fragment of the ligand with another fragment from a database of synthetically accessible scaffolds. A GRID-based interaction profile of the receptor and geometrical descriptions of a ligand scaffold are used to obtain new scaffolds with different structural features and are able to replace the original scaffold in the protein-ligand complex. An enrichment study was successfully performed verifying the ability of SHOP to find known active CDK2 scaffolds in a database. Additionally, SHOP was used for suggesting new inhibitors of p38 MAP kinase. Four p38 complexes were used to perform six scaffold searches. Several new scaffolds were suggested, and the resulting compounds were successfully docked into the query proteins.

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.