Structure-based ligand design to overcome CYP inhibition in drug discovery projects

Computational Binding Study Hints at Ecdysone 20-Mono-Oxygenase as the Hitherto Unknown Target for Ring C-Seco Limonoid-Type Insecticides

The insecticidal property of ring C-seco limonoids has been discovered empirically and the target protein identified, but, to date, the molecular mechanism of action has not been described at the atomic scale. We elucidate on computational grounds whether nine C-seco limonoids present sufficiently high affinity to bind specifically with the putative target enzyme of the insects (ecdysone 20-monooxygenase). To this end, 3D models of ligands and the receptor target were generated and their interaction energies estimated by docking simulations. As a proof of concept, the tetrahydro-isoquinolinyl propenamide derivative QHC is the reference ligand bound to aldosterone synthase in the complex with PDB entry 4ZGX. It served as the 3D template for target modeling via homology. QHC was successfully docked back to its crystal pose in a one-digit nanomolar range. The reported experimental binding affinities span over the nanomolar to lower micromolar range. All nine limonoids were found with strong affinities in the range of -9 < ΔG < -13 kcal/mol. The molt hormone ecdysone showed a comparable ΔG energy of -12 kcal/mol, whereas -11 kcal/mol was the back docking result for the liganded crystal 4ZGX. In conclusion, the nine C-seco limonoids were strong binders on theoretical grounds in an activity range between a ten-fold lower to a ten-fold higher concentration level than insecticide ecdysone with its known target receptor. The comparable or even stronger binding hints at ecdysone 20-monooxygenase as their target biomolecule. Our assumption, however, is in need of future experimental confirmation before conclusions with certainty can be drawn about the true molecular mechanism of action for the C-seco limonoids under scrutiny.

Cheminformatics identification of modulators of key carbohydrate-metabolizing enzymes from C. cujete for type-2 diabetes mellitus intervention

Purpose

The therapeutic use of oral hypoglycaemic agents in the management of type-2 diabetes mellitus (T2DM) is without adverse effects; thus, calls for alternative and novel candidates from natural products in medicinal plants.

Method

The study explored molecular docking and molecular dynamics (MD) simulation approaches to identify key antidiabetic metabolites from Crescentia cujete.

Results

Molecular docking results identified four and/or five best compounds against each target enzyme (alpha-glucosidase, dipeptidyl peptidase-IV, aldose reductase, and protein tyrosine phosphatase-1B (PTP-1B)) implicated in diabetes. The resulting complexes (except against PTP-1B) had higher docking scores above respective standards (acarbose, Diprotin A, ranirestat). The MD simulation results revealed compounds such as benzoic acid (-48.414 kcal/mol) and phytol (-45.112 kcal/mol) as well as chlorogenic acid (-42.978 kcal/mol) and naringenin (-31.292 kcal/mol) had higher binding affinities than the standards [acarbose (-28.248 kcal/mol), ranirestat (-21.042 kcal/mol)] against alpha-glucosidase and aldose reductase, respectively while Diprotin A (-45.112 kcal/mol) and ursolic acid (-18.740 kcal/mol) presented superior binding affinities than the compounds [luteolin (-41.957 kcal/mol and naringenin (-16.518 kcal/mol)] against DPP-IV and PTP-1B respectively.

Conclusion

While isoflavone (alpha-glucosidase), xylocaine (DPP-IV), luteolin (aldose reductase,) and chlorogenic acid (PTP-1B) were affirmed as the best inhibitors of respective enzyme targets, luteolin, and chlorogenic acid may be suggested and proposed as probable candidates against T2DM and related retinopathy complication based on their structural stability, compactness and affinity for three (DPP-IV, aldose reductase, and PTP-1B) of the four targets investigated. Further studies are warranted in vitro and in vivo on the antihyperglycaemic effects of these drug candidates.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40200-023-01249-7.

Molecular Modeling Identification of Key Secondary Metabolites from Xylopia aethiopica as Promising Therapeutics Targeting Essential Measles Viral Proteins

This study computationally screened three key compounds (vanillin (VAN), oxophoebine (OPB), and dihydrochalcone (DHC)) derived from Xylopia aethiopica (Guinea pepper), a medicinal plant with known antiviral activity, against key druggable measles virus (MV) proteins (fusion protein (FUP), haemagglutinin protein (HMG), and phosphoprotein (PSP)). Each molecular species was subjected to a 100 ns molecular dynamics (MD) simulation following docking, and a range of postdynamic parameters including free binding energy and pharmacokinetic properties were determined. The docking scores of the resulting OPB-FUP (-5.4 kcal/mol), OPB-HMG (-8.1 kcal/mol), and OPB-PSP (-8.0 kcal/mol) complexes were consistent with their respective binding energy values (-25.37, -28.74, and -40.68 kcal/mol), and higher than that of the reference standard, ribavirin (RBV) in each case. Furthermore, all the investigated compounds were thermodynamically compact and stable, especially HMG of MV, and this observation could be attributed to the resulting intermolecular interactions in each system. Overall, OPB may possess inhibitory properties against MV glycoproteins (FUP and HMG) and PSP that play important roles in the replication of MV and measles pathogenesis. While OPB could serve as a scaffold for the development of novel MV fusion and entry inhibitors, further in vitro and in vivo evaluation is highly recommended.

Network Pharmacology- and Molecular Dynamics Simulation-Based Bioprospection of Aspalathus linearis for Type-2 Diabetes Care

The medicinal herb Aspalathus linearis (rooibos) is globally recognized in type-2 diabetes mellitus (T2DM) treatment due to its known and distinctive compounds. This work utilized network pharmacology (NP) coupled with molecular dynamics simulation in gaining new insight into the anti-diabetic molecular mechanism of action of rooibos teas. It looked at the interactions between rooibos constituents with various relevant protein receptors and signaling routes associated with T2DM progression. The initial analysis revealed 197 intersecting gene targets and 13 bioactive rooibos constituents linked to T2DM. The interactions between proteins and compounds to the target matrix were generated with the Cystoscope platform and STRING database. These analyses revealed intersecting nodes active in T2DM and hypoxia-inducible factor 1 (HIF-1) as an integral receptors target. In addition, KEGG analysis identified 11 other pathways besides the hub HIF-1 signaling route which may also be targeted in T2DM progression. In final molecular docking and dynamics simulation analysis, a significant binding affinity was confirmed for key compound-protein matrices. As such, the identified rooibos moieties could serve as putative drug candidates for T2DM control and therapy. This study shows rooibos constituents' interaction with T2DM-linked signaling pathways and target receptors and proposes vitexin, esculin and isovitexin as well as apigenin and kaempferol as respective pharmacologically active rooibos compounds for the modulation of EGFR and IGF1R in the HIF-1 signaling pathway to maintain normal homeostasis and function of the pancreas and pancreatic β-cells in diabetics.

Generation of new inhibitors of selected cytochrome P450 subtypes- In silico study

Physicochemical and pharmacokinetic compound profile has crucial impact on compound potency to become a future drug. Ligands with desired activity profile cannot be used for treatment if they are characterized by unfavourable physicochemical or ADMET properties. In the study, we consider metabolic stability and focus on selected subtypes of cytochrome P450 - proteins, which take part in the first phase of compound transformations in the organism. We develop a protocol for generation of new potential inhibitors of selected cytochrome isoforms. Its subsequent stages are composed of generation and assessment of new derivatives of known cytochrome inhibitors, docking and evaluation of the compound possible inhibition on the basis of the obtained ligand-protein complexes. Besides the library of new potential agents inhibiting particular cytochrome subtypes, we also prepare a graph neural network that predicts the change in activity for all modifications of the starting molecule. In addition, we perform a systematic statistical study on the influence of particular substitutions on the potential inhibition properties of generated compounds (both mono- and di-substitutions are considered), provide explanations of the inhibitory predictions and prepare an on-line visualization platform enabling manual inspection of the results. The developed methodology can greatly support the design of new cytochrome P450 inhibitors with the overarching goal of generation of new metabolically stable compounds. It enables instant evaluation of possible compound-cytochrome interactions and selection of ligands with the highest potential of possessing desired biological activity.

Computer-Aided (In Silico) Modeling of Cytochrome P450-Mediated Food–Drug Interactions (FDI)

Modifications of the activity of Cytochrome 450 (CYP) enzymes by compounds in food might impair medical treatments. These CYP-mediated food–drug interactions (FDI) play a major role in drug clearance in the intestine and liver. Inter-individual variation in both CYP expression and structure is an important determinant of FDI. Traditional targeted approaches have highlighted a limited number of dietary inhibitors and single-nucleotide variations (SNVs), each determining personal CYP activity and inhibition. These approaches are costly in time, money and labor. Here, we review computational tools and databases that are already available and are relevant to predicting CYP-mediated FDIs. Computer-aided approaches such as protein–ligand interaction modeling and the virtual screening of big data narrow down hundreds of thousands of items in databanks to a few putative targets, to which the research resources could be further directed. Structure-based methods are used to explore the structural nature of the interaction between compounds and CYP enzymes. However, while collections of chemical, biochemical and genetic data are available today and call for the implementation of big-data approaches, ligand-based machine-learning approaches for virtual screening are still scarcely used for FDI studies. This review of CYP-mediated FDIs promises to attract scientists and the general public.

In Silico Design and SAR Study of Dibenzyl Trisulfide Analogues for Improved CYP1A1 Inhibition

Dibenzyl trisulfide (DTS) is a natural compound with potential cancer‐preventive properties occurring in Petiveria alliacea L., an ethnomedicinal plant native to the Americas. Previous studies revealed its inhibitory activity toward cytochrome P450 (CYP)1 enzymes, key in the activation of environmental pollutants. Accordingly, the aim of this study was to design novel DTS analogues, aimed at improving not only inhibitory activity, but also specificity toward CYP1A1. This was achieved by targeting interactions with CYP1A1 residues of identified importance. Three‐dimensional structures for the novel analogues were subjected to molecular docking with several CYP isoforms, before being ranked in terms of binding affinity to CYP1A1. With three hydrogen bond donors, two hydrogen bond acceptors, a molecular mass of 361 Da, and a log P of 3.72, the most promising DTS analogue obeys Lipinski's rule of five. Following synthesis and in vitro validation of its CYP1A1‐inhibitory properties, this compound may be useful in future cancer‐preventive approaches.

A Comprehensive Analysis of Human CYP3A4 Crystal Structures as a Potential Tool for Molecular Docking-Based Site of Metabolism and Enzyme Inhibition Studies

The notable ability of human liver cytochrome P450 3A4 (CYP3A4) to metabolize diverse xenobiotics encourages researchers to explore in-depth the mechanism of enzyme action. Numerous CYP3A4 protein crystal structures have been deposited in protein data bank (PDB) and are majorly used in molecular docking analysis. The quality of the molecular docking results depends on the three-dimensional CYP3A4 protein crystal structures from the PDB. Present review endeavors to provide a brief outline of some technical parameters of CYP3A4 PDB entries as valuable information for molecular docking research. PDB entries between 22 April 2004 and 2 June 2021 were compiled and the active sites were thoroughly observed. The present review identified 76 deposited PDB entries and described basic information that includes CYP3A4 from human genetic, Escherichia coli (E. coli) use for protein expression, crystal structure obtained from X-ray diffraction method, taxonomy ID 9606, Uniprot ID P08684, ligand–protein structure description, co-crystal ligand, protein site deposit and resolution ranges between 1.7[Formula: see text]Å and 2.95[Formula: see text]Å. The observation of protein–ligand interactions showed the various residues on the active site depending on the ligand. The residues Ala305, Ser119, Ala370, Phe304, Phe108, Phe213 and Phe215 have been found to frequently interact with ligands from CYP3A4 PDB. Literature surveys of 17 co-crystal ligands reveal multiple mechanisms that include competitive inhibition, noncompetitive inhibition, mixed-mode inhibition, mechanism-based inhibition, substrate with metabolite, inducer, or combination modes of action. This overview may help researchers choose a trustworthy CYP3A4 protein structure from the PDB database to apply the protein in molecular docking analysis for drug discovery.

Optimization of a Pyrimidinone Series for Selective Inhibition of Ca2+/Calmodulin-Stimulated Adenylyl Cyclase 1 Activity for the Treatment of Chronic Pain

Adenylyl cyclase type 1 (AC1) is involved in signaling for chronic pain sensitization in the central nervous system and is an emerging target for the treatment of chronic pain. AC1 and a closely related isoform AC8 are also implicated to have roles in learning and memory signaling processes. Our team has carried out cellular screening for inhibitors of AC1 yielding a pyrazolyl-pyrimidinone scaffold with low micromolar potency against AC1 and selectivity versus AC8. Structure-activity relationship (SAR) studies led to analogues with cellular IC50 values as low as 0.25 μM, selectivity versus AC8 and other AC isoforms as well as other common neurological targets. A representative analogue displayed modest antiallodynic effects in a mouse model of inflammatory pain. This series represents the most potent and selective inhibitors of Ca2+/calmodulin-stimulated AC1 activity to date with improved drug-like physicochemical properties making them potential lead compounds for the treatment of inflammatory pain.

Identification of Flavonoid C-Glycosides as Promising Antidiabetics Targeting Protein Tyrosine Phosphatase 1B

Protein tyrosine phosphatase 1B (PTP1B), a negative regulator of the insulin signaling pathway, has gained attention as a validated druggable target in the management of type 2 diabetes mellitus (T2DM). The lack of clinically approved PTP1B inhibitors has continued to prompt research in plant-derived therapeutics possibly due to their relatively lesser toxicity profiles. Flavonoid C-glycosides are one of the plant-derived metabolites gaining increased relevance as antidiabetic agents, but their possible mechanism of action remains largely unknown. This study investigates the antidiabetic potential of flavonoid C-glycosides against PTP1B in silico and in vitro. Of the seven flavonoid C-glycosides docked against the enzyme, three compounds (apigenin, vitexin, and orientin) had the best affinity for the enzyme with a binding score of -7.3 kcal/mol each, relative to -7.4 kcal/mol for the reference standard, ursolic acid. A further probe (in terms of stability, flexibility, and compactness) of the complexes over a molecular dynamics time study of 100 ns for the three compounds suggested orientin as the most outstanding inhibitor of PTP1B owing to its overall -34.47 kcal/mol binding energy score compared to ursolic acid (-19.24 kcal/mol). This observation was in accordance with the in vitro evaluation result, where orientin had a half maximal inhibitory concentration (IC₅₀) of 0.18 mg/ml relative to 0.13 mg/ml for the reference standard. The kinetics of inhibition of PTP1B by orientin was mixed-type with V _max and K _m values of 0.004 μM/s and 0.515 μM. Put together, the results suggest orientin as a potential PTP1B inhibitor and could therefore be further explored in the management T2DM as a promising therapeutic agent.

Machine learning-driven identification of drugs inhibiting cytochrome P450 2C9

Cytochrome P450 2C9 (CYP2C9) is a major drug-metabolizing enzyme that represents 20% of the hepatic CYPs and is responsible for the metabolism of 15% of drugs. A general concern in drug discovery is to avoid the inhibition of CYP leading to toxic drug accumulation and adverse drug-drug interactions. However, the prediction of CYP inhibition remains challenging due to its complexity. We developed an original machine learning approach for the prediction of drug-like molecules inhibiting CYP2C9. We created new predictive models by integrating CYP2C9 protein structure and dynamics knowledge, an original selection of physicochemical properties of CYP2C9 inhibitors, and machine learning modeling. We tested the machine learning models on publicly available data and demonstrated that our models successfully predicted CYP2C9 inhibitors with an accuracy, sensitivity and specificity of approximately 80%. We experimentally validated the developed approach and provided the first identification of the drugs vatalanib, piriqualone, ticagrelor and cloperidone as strong inhibitors of CYP2C9 with IC values <18 μM and sertindole, asapiprant, duvelisib and dasatinib as moderate inhibitors with IC50 values between 40 and 85 μM. Vatalanib was identified as the strongest inhibitor with an IC50 value of 0.067 μM. Metabolism assays allowed the characterization of specific metabolites of abemaciclib, cloperidone, vatalanib and tarafenacin produced by CYP2C9. The obtained results demonstrate that such a strategy could improve the prediction of drug-drug interactions in clinical practice and could be utilized to prioritize drug candidates in drug discovery pipelines.

Phenolics Profiling of Carpobrotus edulis (L.) N.E.Br. and Insights into Molecular Dynamics of Their Significance in Type 2 Diabetes Therapy and Its Retinopathy Complication

Adverse effects associated with synthetic drugs in diabetes therapy has prompted the search for novel natural lead compounds with little or no side effects. Effects of phenolic compounds from Carpobrotus edulis on carbohydrate-metabolizing enzymes through in vitro and in silico methods were assessed. Based on the half-maximal inhibitory concentrations (IC₅₀), the phenolic extract of the plant had significant (p < 0.05) in vitro inhibitory effect on the specific activity of alpha-amylase (0.51 mg/mL), alpha-glucosidase (0.062 mg/mL) and aldose reductase (0.75 mg/mL), compared with the reference standards (0.55, 0.72 and 7.05 mg/mL, respectively). Molecular interactions established between the 11 phenolic compounds identifiable from the HPLC chromatogram of the extract and active site residues of the enzymes revealed higher binding affinity and more structural compactness with procyanidin (−69.834 ± 6.574 kcal/mol) and 1,3-dicaffeoxyl quinic acid (−42.630 ± 4.076 kcal/mol) as potential inhibitors of alpha-amylase and alpha-glucosidase, respectively, while isorhamnetin-3-O-rutinoside (−45.398 ± 4.568 kcal/mol) and luteolin-7-O-beta-d-glucoside (−45.102 ± 4.024 kcal/mol) for aldose reductase relative to respective reference standards. Put together, the findings are suggestive of the compounds as potential constituents of C. edulis phenolic extract responsible for the significant hypoglycemic effect in vitro; hence, they could be exploited in the development of novel therapeutic agents for type-2 diabetes and its retinopathy complication.

An insight on the nature of biochemical interactions between glycyrrhizin, myricetin and CYP3A4 isoform

Drug interaction studies are imperative to gain insights into the beneficial or harmful effects of therapeutic and dietary agents. This study investigated the mechanism of modulatory roles of glycyrrhizin (GLH) and myricetin (MYC) on the human CYP3A4 isoform using in silico and in vitro methods. While MYC had concentration-dependent inhibitory effect on CYP3A4 (IC₅₀ : 10.5 ± 0.55 μM) with characteristic K_m and V_max values of 1.13 μM and 1.54 nM/min, respectively, GLH exhibited no inhibitory effect on CYP3A4 activity in vitro. These observations are consistent with the results of in silico evaluations where the effect of MYC compared well with that of ketoconazole (a known CYP3A4 inhibitor) against CYP3A4. Overall, the established interactions between the study compounds and CYP3A4 could potentiate clinically vital drug-drug interactions and has lent credence to the mechanism of modulatory effect of MYC and GLH on CYP3A4 that could guide their safe use as therapeutic agents. PRACTICAL IMPLICATIONS: Myricetin (MYR) and glycyrrhizin (GLH) occur freely in commonly ingested foods and their supplements are recommended for the treatment of several debilitating diseases such as diabetes, cancer, and cardiovascular complications. This study provided an insight on the possible interactions that could be established between these compounds (MYR and GLH) and CYP3A4 when ingested and metabolized by the liver. The results suggested possibilities of potential clinical drug-drug interactions and advocates for their cautious use within the therapeutic dose in food supplements or medications to avoid probable liver damage.

CYP2J2 Molecular Recognition: A New Axis for Therapeutic Design

Cytochrome P450 (CYP) epoxygenases are a special subset of heme-containing CYP enzymes capable of performing the epoxidation of polyunsaturated fatty acids (PUFA) and the metabolism of xenobiotics. This dual functionality positions epoxygenases along a metabolic crossroad. Therefore, structure-function studies are critical for understanding their role in bioactive oxy-lipid synthesis, drug-PUFA interactions, and for designing therapeutics that directly target the epoxygenases. To better exploit CYP epoxygenases as therapeutic targets, there is a need for improved understanding of epoxygenase structure-function. Of the characterized epoxygenases, human CYP2J2 stands out as a potential target because of its role in cardiovascular physiology. In this review, the early research on the discovery and activity of epoxygenases is contextualized to more recent advances in CYP epoxygenase enzymology with respect to PUFA and drug metabolism. Additionally, this review employs CYP2J2 epoxygenase as a model system to highlight both the seminal works and recent advances in epoxygenase enzymology. Herein we cover CYP2J2's interactions with PUFAs and xenobiotics, its tissue-specific physiological roles in diseased states, and its structural features that enable epoxygenase function. Additionally, the enumeration of research on CYP2J2 identifies the future needs for the molecular characterization of CYP2J2 to enable a new axis of therapeutic design.

Molecular dynamic mechanism(s) of inhibition of bioactive antiviral phytochemical compounds targeting cytochrome P450 3A4 and P-glycoprotein

P-glycoprotein (ABCB1) and cytochrome P450 3A4 (CYP3A4) metabolize almost all known human immunodeficiency virus' protease inhibitor drugs (PIs). Over induction of these proteins' activities has been linked to rapid metabolism of PIs which are then pumped out of the circulatory system, eventually leading to drug-resistance in HIV-positive patients. This study aims to determine, with the use of computational tools, the inhibitory potential of four phytochemical compounds (PCs) (epigallocatechin gallate (EGCG), kaempferol-7-glucoside (K7G), luteolin (LUT) and ellagic acid (EGA)) in inhibiting the activities of these drug-metabolizing proteins. The comparative analysis of the MM/GBSA results revealed that the binding affinity (ΔG_bind) of EGCG and K7G for CYP3A4 and ABCB1 are higher than LUT and EGA and fall between the ΔG_bind of the inhibitors of CYP3A4 and ABCB1 (Ritonavir (strong inhibitor) and Lopinavir (moderate inhibitor)). The structural analysis (RMSD, RMSF, RoG and protein-ligand interaction plots) also confirmed that EGCG and K7G showed similar inhibitory activities with the inhibitors. The study has shown that EGCG and K7G have inhibitory activities against the two proteins and assumes they could decrease intracellular efflux of PIs, consequently increasing the optimal concentration of PIs in the systemic circulation.Communicated by Ramaswamy H. Sarma.

Simultaneous Pharmacokinetics Estimation of Nateglinide and Pioglitazone by RP-HPLC: Computational Study to Unlock the Synergism

Nateglinide (NAT) and Pioglitazone (PIO) are an antidiabetic drugs combination and currently under clinical trial in countries like Japan. In this study, an alternative, a simple, sensitive high-performance liquid chromatography method has been developed (limit of detection: 15 ng/mL and limit of quantification: 50 ng/mL) for simultaneous estimation of this drug combination in rat plasma. Most remarkably, bioavailability of NAT has been increased markedly on coadministration with PIO, than when it was administered alone. Thus, PIO is assumed to retard the catabolism of NAT by inhibiting metabolic liver-microsomal enzyme, especially CYP2C9. Using a Waters Nova-Pak C 18 column (150 × 3.9 mm, 4 μm) and a mobile phase of acetonitrile: 10 mM KH2PO4 (60: 40, V/V (volume by volume)) pH 3.5, the analysis was performed at 210 nm with a flow rate of 1.5 mL/min. In silico docking via molecular dynamics simulation revealed that NAT-CYP2C9 binding affinity may be reduced after PIO attachment, presumably due to the binding site overlapping of the two drugs. Thus, it has been proposed that NAT and PIO may be an efficient synergistic fixed dose combination against diabetes mellitus, and the above method can foster a simple but highly sensitive bioanalytical estimation for routine analysis.

Molecular Docking and Preclinical Study of Five-Membered S,S-Palladaheterocycle as Hepatoprotective Agent

Purpose: In order to investigate mechanisms underlying the hepatoprotective action of S,Spalladaheterocycle, inhibition of cytochromes P450 has been modeled by molecular docking of four palladaheterocycle stereoisomers to the active sites of an enzymatic oxidase system. To obtain a deeper insight into biochemical aspects providing a basis for the therapeutic effects of five-membered palladacycles (as mixture of stereoisomers), a number of preclinical trials has been conducted

Methods: 2D and 3D structures of palladaheterocycle stereoisomers were obtained via converting into SDF files by means of software MarvinSketch. Binding of palladaheterocycle at the active sites of cytochromes P450 2E1 and P450 2C9 has been studied by molecular docking using LeadIT 2.3.2. Hepatoprotective activity of palladaheterocycle at 2.5, 25 and 250 mg/kg doses has been studied based on a model of acute intoxication by CCl₄ using in vivo methods.

Results: By molecular docking it was identify amino acid fragments responsible for binding with palladacyclic isomers. The tested compound is comparable, in terms of its activity to the hepatoprotective drug SAM according to the in vivo and in vitro experiments such as animal survival data, the efficiency of correction of the cytolytic syndrome, the liver excretory function, carbohydrate, protein and lipid metabolism, and the correction efficiency of the liver antitoxic function (the latter has been determined based on the results of a hexobarbital control experiment).

Conclusion: Taking into account results obtained in vivo, in vitro and in silico, it can be concluded that the five-membered S,S-palladaheterocycle effectively protect the liver against acute damage caused by CCl₄ , via activation of catalase and glucuronyltransferase, as well as via inhibition of the oxidative stress enzymes.

Differing Membrane Interactions of Two Highly Similar Drug-Metabolizing Cytochrome P450 Isoforms: CYP 2C9 and CYP 2C19

The human cytochrome P450 (CYP) 2C9 and 2C19 enzymes are two highly similar isoforms with key roles in drug metabolism. They are anchored to the endoplasmic reticulum membrane by their N-terminal transmembrane helix and interactions of their cytoplasmic globular domain with the membrane. However, their crystal structures were determined after N-terminal truncation and mutating residues in the globular domain that contact the membrane. Therefore, the CYP-membrane interactions are not structurally well-characterized and their dynamics and the influence of membrane interactions on CYP function are not well understood. We describe herein the modeling and simulation of CYP 2C9 and CYP 2C19 in a phospholipid bilayer. The simulations revealed that, despite high sequence conservation, the small sequence and structural differences between the two isoforms altered the interactions and orientations of the CYPs in the membrane bilayer. We identified residues (including K72, P73, and I99 in CYP 2C9 and E72, R73, and H99 in CYP 2C19) at the protein-membrane interface that contribute not only to the differing orientations adopted by the two isoforms in the membrane, but also to their differing substrate specificities by affecting the substrate access tunnels. Our findings provide a mechanistic interpretation of experimentally observed effects of mutagenesis on substrate selectivity.

Analysis of the Molecular Interactions between Cytochromes P450 3A4 and 1A2 and Aflatoxins: A Docking Study

Mycotoxins known as aflatoxins (AF) are produced as a secondary metabolite by some species of Aspergillus fungi. They are considered carcinogenic, hepatotoxic, teratogenic, and mutagenic. In this study, the molecular structure, chemical reactivity, and charge transfer values of AFB1, B2, G1, and G2 were analyzed using density functional theory. Different methodologies—B3LYP/6-311G(d,p) and M06-2X/6-311G(d,p)—were applied for geometrical calculations. Chemical reactivity parameters were used in the calculation of charge transfer values during the interaction between protein and ligand. The binding energy, the electrostatic interactions, and the amino acids of the active site were determined by molecular docking analysis between AF and cytochromes P450 (3A4 and 1A2), employing different PDB files (CYP3A4:1TQN, 2V0M, 4NY4 and 1W0E, and CYP1A2:2HI4). Molecular docking analysis indicated that the central rings of the AF are involved in the interaction with the HEM group of the active site. The differences in the molecular structure of the AF affect their position regarding the HEM group. The resulting configurations presented considerable variation in the amino acids and the position of the coupling. The charge transfer values showed that there is oxidative damage inside the active site and that the HEM group is responsible for the main charge transferences.

Active-site differences between substrate-free and ritonavir-bound cytochrome P450 (CYP) 3A5 reveal plasticity differences between CYP3A5 and CYP3A4

Cytochrome P450 (CYP) 3A4 is a major contributor to hepatic drug and xenobiotic metabolism in human adults. The related enzyme CYP3A5 is also expressed in adult liver and has broader age and tissue distributions. However, CYP3A5 expression is low in most Caucasians because of the prevalence of an allele that leads to an incorrectly spliced mRNA and premature termination of translation. When expressed, CYP3A5 expands metabolic capabilities and can augment CYP3A4-mediated drug metabolism, thereby reducing drug efficacy and potentially requiring dose adjustments. The extensive role of CYP3A4 in drug metabolism reflects in part the plasticity of the substrate-free enzyme to enlarge its active site and accommodate very large substrates. We have previously shown that the structure of the CYP3A5-ritonavir complex differs substantially from that of the CYP3A4-ritonavir complex. To better understand whether these differences are conserved in other CYP3A5 structures and how they relate to differential plasticity, we determined the X-ray crystallographic structure of the CYP3A5 substrate-free complex to 2.20 Å resolution. We observed that this structure exhibits a much larger active site than substrate-free CYP3A4 and displays an open substrate access channel. This reflected in part a lower trajectory of the helix F-F' connector in CYP3A4 and more extensive π-CH interactions between phenylalanine residues forming the roof of the active-site cavity than in CYP3A5. Comparison with the CYP3A5-ritonavir complex confirmed conserved CYP3A5 structural features and indicated differences in plasticity between CYP3A4 and CYP3A5 that favor alternative ritonavir conformations.

Amelioration of mechanism-based inactivation of CYP3A4 by a H-PGDS inhibitor

This work describes the rational amelioration of mechanism-based inactivation (MBI) of Cytochrome P450 (CYP) 3A4 in a human hematopoietic prostaglandin D synthase (hH-PGDS) inhibitor (cpd 1). We utilized metabolism reports in order to check if patterns in the metabolism of 1 and similar compounds by CYP3A4 could be deciphered. Then we used structure based design, first modifying the CYP3A4 crystal structure (pdb code: 4NY4) by adding an oxyferryl moiety to the heme, followed by validating the modified structure to obtain the 1' and 4 position oxidation products of midazolam and then recapitulating the metabolism patterns deciphered previously for 1 and analogs. We checked if the pattern deciphered could lead to a putative reactive moiety. Finally we used the docking pose of 1 into this model of the modified CYP3A4 crystal structure to guide transformation of 1 into MBI-free H-PGDS inhibitors.

Can the inhibition of cytochrome P450 in aquatic invertebrates due to azole fungicides be estimated with in silico and in vitro models and extrapolated between species?

Azole fungicides, designed to halt fungal growth by specific inhibition of fungal cytochrome P450 (CYP51), inhibit cytochrome P450s involved in the metabolism of xenobiotics in several non-target organisms thus raising environmental concern. The present study investigates the degree by which inhibition strengths of azoles toward cytochrome P450 in rat liver, the insect Chironomus riparius larvae and the snail Lymnaea stagnalis can be extrapolated from estimated in silico affinities. Azoles' affinities toward human cytochrome P450 isoforms involved in xenobiotic metabolism (CYP3A4, CYP2C9 and CYP2D6) as well as fungal CYP51 were estimated with a ligand-protein docking model based on the ChemScore scoring function. Estimated affinities toward the selected enzymatic structures correlated strongly with measured inhibition strengths in rat liver (ChemScore vs. logIC₅₀ among cytochrome P450 isoforms: -0.662 < r < -0.891, n = 17 azoles), while weaker correlations were found for C. riparius larvae (-0.167 < r < -0.733, n = 9) and L. stagnalis (-0.084 < r < -0.648, n = 8). Inhibition strengths toward C. riparius and rat liver activities were found to be highly correlated to each other (r: 0.857) while no significant relationship was found between either of the species and L. stagnalis. The inhibition of cytochrome P450 due to azole fungicides could be estimated in vitro and to a lesser extent in silico for C. riparius but not for L. stagnalis, possibly due to different enzymatic susceptibility toward azole inhibition among the species.

Steroid bioconjugation to a CYP3A4 allosteric site and its effect on substrate binding and coupling efficiency

Human cytochrome P450 3A4 (CYP3A4) is an important drug metabolizing enzyme involved in a number of drug-drug and food-drug interactions. As such, much effort has been devoted into investigating its mechanism of interaction with ligands. CYP3A4 has one of the highest levels of substrate promiscuity for an enzyme, and can even bind multiple ligands simultaneously. The location and orientation of these ligands depend on the chemical structure and stoichiometry, and are generally poorly understood. In the case of the steroid testosterone, up to three copies of the molecule can associate with the enzyme at once, likely two in the active site and one at a postulated allosteric site. Recently, we demonstrated that steroid bioconjugation at the allosteric site results in an increase in activity of CYP3A4 toward testosterone and 7-benzyloxy-4-trifluoromethylcoumarin oxidation. Here, using the established bioconjugation methodology, we show how steroid bioconjugation at the allosteric site affects the heme spin state, the binding affinity (K_S) of CYP3A4 for testosterone, as well as the enzyme coupling efficiency.

Predicting mammalian metabolism and toxicity of pesticides in silico

Pesticides must be effective to be commercially viable but they must also be reasonably safe for those who manufacture them, apply them, or consume the food they are used to produce. Animal testing is key to ensuring safety, but it comes late in the agrochemical development process, is expensive, and requires relatively large amounts of material. Surrogate assays used as in vitro models require less material and shift identification of potential mammalian toxicity back to earlier stages in development. Modern in silico methods are cost-effective complements to such in vitro models that make it possible to predict mammalian metabolism, toxicity and exposure for a pesticide, crop residue or other metabolite before it has been synthesized. Their broader use could substantially reduce the amount of time and effort wasted in pesticide development. This contribution reviews the kind of in silico models that are currently available for vetting ideas about what to synthesize and how to focus development efforts; the limitations of those models; and the practical considerations that have slowed development in the area. Detailed discussions are provided of how bacterial mutagenicity, human cytochrome P450 (CYP) metabolism, and bioavailability in humans and rats can be predicted. © 2018 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Insights into molecular mechanisms of drug metabolism dysfunction of human CYP2C9*30

Cytochrome P450 2C9 (CYP2C9) metabolizes about 15% of clinically administrated drugs. The allelic variant CYP2C9*30 (A477T) is associated to diminished response to the antihypertensive effects of the prodrug losartan and affected metabolism of other drugs. Here, we investigated molecular mechanisms involved in the functional consequences of this amino-acid substitution. Molecular dynamics (MD) simulations performed for the active species of the enzyme (heme in the Compound I state), in the apo or substrate-bound state, and binding energy analyses gave insights into altered protein structure and dynamics involved in the defective drug metabolism of human CYP2C9.30. Our data revealed an increased rigidity of the key Substrate Recognition Sites SRS1 and SRS5 and shifting of the β turn 4 of SRS6 toward the helix F in CYP2C9.30. Channel and binding substrate dynamics analyses showed altered substrate channel access and active site accommodation. These conformational and dynamic changes are believed to be involved in the governing mechanism of the reduced catalytic activity. An ensemble of representative conformations of the WT and A477T mutant properly accommodating drug substrates were identified, those structures can be used for prediction of new CYP2C9 and CYP2C9.30 substrates and drug-drug interactions.

Identification of karanjin isolated from the Indian beech tree as a potent CYP1 enzyme inhibitor with cellular efficacy via screening of a natural product repository

CYP1A1 is thought to mediate carcinogenesis in oral, lung and epithelial cancers. In order to identify a CYP1A1 inhibitor from an edible plant, 394 natural products in the IIIM's natural product repository were screened, at 10 μM concentration, using CYP1A1-Sacchrosomes™ (i.e. microsomal enzyme isolated from recombinant baker's yeast). Twenty-seven natural products were identified that inhibited 40-97% of CYP1A1's 7-ethoxyresorufin-O-deethylase activity. The IC50 values of the 'hits', belonging to different chemical scaffolds, were determined. Their selectivity was studied against a panel of 8 CYP-Sacchrosomes™. In order to assess cellular efficacy, the 'hits' were screened for their capability to inhibit CYP enzymes expressed within live recombinant human embryonic kidney (HEK293) cells from plasmids encoding specific CYP genes (1A2, 1B1, 2C9, 2C19, 2D6, 3A4). Isopimpinellin (IN-475; IC50, 20 nM) and karanjin (IN-195; IC50, 30 nM) showed the most potent inhibition of CYP1A1 in human cells. Isopimpinellin is found in celery, parsnip, fruits and in the rind and pulp of limes whereas different parts of the Indian beech tree, which contain karanjin, have been used in traditional medicine. Both isopimpinellin and karanjin negate the cellular toxicity of CYP1A1-mediated benzo[a]pyrene. Molecular docking and molecular dynamic simulations with CYP isoforms rationalize the observed trends in the potency and selectivity of isopimpinellin and karanjin.

Molecular dynamics investigations of regioselectivity of anionic/aromatic substrates by a family of enzymes: a case study of diclofenac binding in CYP2C isoforms

The CYP2C subfamily is of particular importance in the metabolism of drugs, food toxins, and procarcinogens. Like other P450 subfamilies, 2C enzymes share a high sequence identity, but significantly contribute in different ways to hepatic capacity to metabolize drugs. They often metabolize the same substrate to more than one product with different catalytic sites. Because it is challenging to characterize experimentally, much still remains unknown about the reason for why the substrate regioselectivity of these closely related subfamily members is different. Here, we have investigated the structural features of CYP2C8, CYP2C9, and CYP2C19 bound with their shared substrate diclofenac to elucidate the underlying molecular mechanism for the substrate regioselectivity of CYP2C subfamily enzymes. The obtained results demonstrate how a sequence divergence for the active site residues causes heterogeneous variations in the secondary structures and in major tunnel selections, and further affects the shape and chemical properties of the substrate-binding site. Structural analysis and free energy calculations showed that the most important determinants of regioselectivity among the CYP2C isoforms are the geometrical features of the active sites, as well as the hydrogen bonds and the hydrophobic interactions, mainly presenting as the various locations of Arg108 and substitutions of Phe205 for Ile205 in CYP2C8. The MM-GB/SA calculations combined with PMF results accord well with the experimental KM values, bridging the gap between the theory and the experimentally observed results of binding affinity differences. The present study provides important insights into the structure-function relationships of CYP2C subfamily enzymes, the knowledge of ligand binding characteristics and key residue contributions could guide future experimental and computational work on the synthesis of drugs with better pharmacokinetic properties so that CYP interactions could be avoided.

The effect of the intramolecular C–H⋯O interactions on the conformational preferences of bis-arylsulfones – 5-HT6 receptor antagonists and beyond

The development of compounds with enhanced activity and selectivity by a conserved spatial orientation of the pharmacophore elements has a long history in medicinal chemistry. Rigidified compounds are an example of this concept. However, the intramolecular interactions were seldom used as a basis for conformational restraints. Here, we show the weak intramolecular interactions that contribute to the relatively well-conserved geometry of N1-arylsulfonyl indole derivatives. The structure analysis along with quantum mechanics calculations revealed a crucial impact of the sulfonyl group on the compound geometry. The weak intramolecular C–H⋯O interaction stabilizes the mutual "facing" orientation of two aromatic fragments. These findings extend the pharmacological interpretation of the sulfonyl group role from the double hydrogen bond acceptor to the conformational scaffold based on intramolecular forces. This feature has, to date, been omitted in in silico drug discovery. Our results should increase the awareness of researchers to consider the conformational preference when designing new compounds or improving computational methods.

The impact of weak intramolecular C–H⋯O interactions on the conformational stability of bis-arylsulfones is discussed, suggesting different role of sulfonyl group in the ligand – 5HT₆ receptor interaction.

Pharmacogenomics of CYP2C9: Functional and Clinical Considerations

CYP2C9 is the most abundant CYP2C subfamily enzyme in human liver and the most important contributor from this subfamily to drug metabolism. Polymorphisms resulting in decreased enzyme activity are common in the CYP2C9 gene and this, combined with narrow therapeutic indices for several key drug substrates, results in some important issues relating to drug safety and efficacy. CYP2C9 substrate selectivity is detailed and, based on crystal structures for the enzyme, we describe how CYP2C9 catalyzes these reactions. Factors relevant to clinical response to CYP2C9 substrates including inhibition, induction and genetic polymorphism are discussed in detail. In particular, we consider the issue of ethnic variation in pattern and frequency of genetic polymorphisms and clinical implications. Warfarin is the most well studied CYP2C9 substrate; recent work on use of dosing algorithms that include CYP2C9 genotype to improve patient safety during initiation of warfarin dosing are reviewed and prospects for their clinical implementation considered. Finally, we discuss a novel approach to cataloging the functional capabilities of rare 'variants of uncertain significance', which are increasingly detected as more exome and genome sequencing of diverse populations is conducted.

Structural Basis of Single-Nucleotide Polymorphisms in Cytochrome P450 2C9

Single-nucleotide polymorphisms in drug-metabolizing cytochrome P450 (CYP) enzymes are important contributors to interindividual differences in drug metabolism leading to adverse drug reactions. Despite their extensive characterization and importance in pharmacogenetics of clinical drugs, the structural basis of CYP polymorphisms has remained scant. Here we report the crystal structures of human CYP2C9 and its polymorphic variants, *3 (I359L) and *30 (A477T), with an antihypertensive drug losartan. The structures show distinct interaction and occupation of losartan in the active site, the access channel, and the peripheral binding site. The I359L substitution located far from the active site remarkably altered the residue side chains near the active site and the access channel, whereas the T477 substitution illustrated hydrogen-bonding interaction with the reoriented side chain of Q214. The results yield structural insights into the reduced catalytic activity of the CYP2C9 variants and have important implications for understanding genetic polymorphisms in CYP-mediated drug metabolism.

The Cytochrome P450 3A4 has three Major Conformations: New Clues to Drug Recognition by this Promiscuous Enzyme

We computed the channels of the 3A4 isoform of the cytochrome P450 3A4 (CYP) on the basis of 24 crystal structures extracted from the Protein Data Bank (PDB). We identified three major conformations (denoted C, O1 and O2) using an enhanced version of the CCCPP software that we developed for the present work, while only two conformations (C and O² ) are considered in the literature. We established the flowchart of definition of these three conformations in function of the structural and physicochemical parameters of the ligand. The channels are characterized with qualitative and quantitative parameters, and not only with their surrounding secondary structures as it is usually done in the literature.

Identification of Potent and Selective CYP1A1 Inhibitors via Combined Ligand and Structure-Based Virtual Screening and Their in Vitro Validation in Sacchrosomes and Live Human Cells

Target structure-guided virtual screening (VS) is a versatile, powerful, and inexpensive alternative to experimental high-throughput screening (HTS). To discover potent CYP1A1 enzyme inhibitors for cancer chemoprevention, a commercial library of 50 000 small molecules was utilized for VS guided by both ligand and structure-based strategies. For experimental validation, 300 ligands were proposed based on combined analysis of fitness scores from ligand based e-pharmacophore screening and docking score, prime MMGB/SA binding affinity and interaction pattern analysis from structure-based VS. These 300 compounds were screened, at 10 μM concentration, for in vitro inhibition of CYP1A1-Sacchrosomes (yeast-derived microsomal enzyme) in the ethoxyresorufin-O-de-ethylase assay. Thirty-two compounds displayed >50% inhibition of CYP1A1 enzyme activity at 10 μM. 2-Phenylimidazo-[1,2-a]quinoline (5121780, 119) was found to be the most potent with 97% inhibition. It also inhibited ∼95% activity of CYP1B1 and CYP1A2, the other two CYP1 enzymes. The compound 5121780 (119) showed high selectivity toward inhibition of CYP1 enzymes with respect to CYP2 and CYP3 enzymes (i.e., there was no detectable inhibition of CYP2D6/CYP2C9/CYP2C19 and CYP3A4 at 10 μM). It was further investigated in live CYP-expressing human cell system, which confirmed that compound 5121780 (119) potently inhibited CYP1A1, CYP1A2, CYP1B1 enzymes with IC₅₀ values of 269, 30, and 56 nM, respectively. Like in Sacchrosomes, inhibition of CYP2D6/CYP2C9/CYP2C19 and CYP3A4 enzymes, expressed within live human cells, could hardly be detected at 10 μM. The compound 119 rescued CYP1A1 overexpressing HEK293 cells from CYP1A1 mediated benzo[a]pyrene (B[a]P) toxicity and also overcame cisplatin resistance in CYP1B1 overexpressing HEK293 cells. Molecular dynamics simulations of 5121780 (119) with CYP1 enzymes was performed to understand the interaction pattern to CYP isoforms. Results indicate that VS can successfully be used to identify promising CYP1A1 inhibitors, which may have potential in the development of novel cancer chemo-preventive agents.

Thermo-kinetic analysis space expansion for cyclophilin-ligand interactions - identification of a new nonpeptide inhibitor using Biacore™ T200

We have established a refined methodology for generating surface plasmon resonance sensor surfaces of recombinant his‐tagged human cyclophilin‐A. Our orientation‐specific stabilisation approach captures his‐tagged protein under ‘physiological conditions’ (150 mm NaCl, pH 7.5) and covalently stabilises it on Ni²⁺‐nitrilotriacetic acid surfaces, very briefly activated for primary amine‐coupling reactions, producing very stable and active surfaces (≥ 95% specific activity) of cyclophilin‐A. Variation in protein concentration with the same contact time allows straightforward generation of variable density surfaces, with essentially no loss of activity, making the protocol easily adaptable for studying numerous interactions; from very small fragments, ~ 100 Da, to large protein ligands. This new method results in an increased stability and activity of the immobilised protein and allowed us to expand the thermo‐kinetic analysis space, and to determine accurate and robust thermodynamic parameters for the cyclophilin‐A–cyclosporin‐A interaction. Furthermore, the increased sensitivity of the surface allowed identification of a new nonpeptide inhibitor of cyclophilin‐A, from a screen of a fragment library. This fragment, 2,3‐diaminopyridine, bound specifically with a mean affinity of 248 ± 60 μm. The X‐ray structure of this 109‐Da fragment bound in the active site of cyclophilin‐A was solved to a resolution of 1.25 Å (PDB: 5LUD), providing new insight into the molecular details for a potential new series of nonpeptide cyclophilin‐A inhibitors.

Quinazoline derivatives as selective CYP1B1 inhibitors

CYP1B1 is implicated to have a role in the development of breast, ovarian, renal, skin and lung carcinomas. It has been suggested that identification of potent and specific CYP1B1 inhibitors can lead to a novel treatment of cancer. Flavonoids have a compact rigid skeleton which fit precisely within the binding cavity of CYP1B1. Systematic isosteric replacement of flavonoid 'O' atom with 'N' atom led to the prediction that a 'quinazoline' scaffold could be the basis for designing potential CYP1B1 inhibitors. A total of 20 quinazoline analogs were synthesized and screened for CYP1B1 and CYP1A1 inhibition in Sacchrosomes™. IC₅₀ determinations of six compounds with capability of inhibiting CYP1B1 identified quinazolines 5c and 5h as the best candidates for CYP1B1 inhibition, with IC₅₀ values in the nM range. Further selectivity studies with homologous CYPs, belonging to the CYP1, CYP2 and CYP3 family of enzymes, showed that the compounds are likely to be free from critical drug-drug interaction liability. Molecular modelling studies were performed to rationalize the observed enzymatic inhibitions. Further biological studies in live yeast and human cells, harboring CYP1A1 and CYP1B1 enzymes, have illustrated the most potent compounds' cellular permeability and capability of potently inhibiting CYP1B1 enzyme expressed within live cells.

Combined Docking and Quantum Chemical Study on CYP-Mediated Metabolism of Estrogens in Man

Long-term exposure to estrogens seriously increases the incidence of various diseases including breast cancer. Experimental studies indicate that cytochrome P450 (CYP) enzymes catalyze the bioactivation of estrogens to catechols, which can exert their harmful effects via various routes. It has been shown that the 4-hydroxylation pathway of estrogens is the most malign, while 2-hydroxylation is considered a benign pathway. It is also known experimentally that with increasing unsaturation of ring B of estrogens the prevalence of the 4-hydroxylation pathway significantly increases. In this study, we used a combination of structural analysis, docking, and quantum chemical calculations at the B3LYP/6-311+G* level to investigate the factors that influence the regioselectivity of estrogen metabolism in man. We studied the structure of human estrogen metabolizing enzymes (CYP1A1, CYP1A2, CYP1B1, and CYP3A4) in complex with estrone using docking and investigated the susceptibility of estrone, equilin, and equilenin (which only differ in the unsaturation of ring B) to undergo 2- and 4-hydroxylation using several models of CYP enzymes (Compound I, methoxy, and phenoxy radical). We found that even the simplest models could account for the experimental difference between the 2- and 4- hydroxylation pathways and thus might be used for fast screening purposes. We also show that reactivity indices, specifically in this case the radical and nucleophilic condensed Fukui functions, also correctly predict the likeliness of estrogen derivatives to undergo 2- or 4-hydroxylation.

Functioning of drug-metabolizing microsomal cytochrome P450s: In silico probing of proteins suggests that the distal heme 'active site' pocket plays a relatively 'passive role' in some enzyme-substrate interactions

PURPOSE

The currently held mechanistic understanding of microsomal cytochrome P450s (CYPs) seeks that diverse drug molecules bind within the deep-seated distal heme pocket and subsequently react at the heme centre. To explain a bevy of experimental observations and meta-analyses, we indulge a hypothesis that involves a "diffusible radical mediated" mechanism. This new hypothesis posits that many substrates could also bind at alternate loci on/within the enzyme and be reacted without the pertinent moiety accessing a bonding proximity to the purported catalytic Fe-O enzyme intermediate.

METHODS

Through blind and heme-distal pocket centered dockings of various substrates and non-substrates (drug molecules of diverse sizes, classes, topographies etc.) of microsomal CYPs, we explored the possibility of access of substrates via the distal channels, its binding energies, docking orientations, distance of reactive moieties (or molecule per se) to/from the heme centre, etc. We investigated specific cases like- (a) large drug molecules as substrates, (b) classical marker drug substrates, (c) class of drugs as substrates (Sartans, Statins etc.), (d) substrate preferences between related and unrelated CYPs, (e) man-made site-directed mutants' and naturally occurring mutants' reactivity and metabolic disposition, (f) drug-drug interactions, (g) overall affinities of drug substrate versus oxidized product, (h) meta-analysis of in silico versus experimental binding constants and reaction/residence times etc.

RESULTS

It was found that heme-centered dockings of the substrate/modulator drug molecules with the available CYP crystal structures gave poor docking geometries and distances from Fe-heme centre. In conjunction with several other arguments, the findings discount the relevance of erstwhile hypothesis in many CYP systems. Consequently, the newly proposed hypothesis is deemed a viable alternate, as it satisfies Occam's razor.

CONCLUSIONS

The new proposal affords expanded scope for explaining the mechanism, kinetics and overall phenomenology of CYP mediated drug metabolism. It is now understood that the heme-iron and the hydrophobic distal pocket of CYPs serve primarily to stabilize the reactive intermediate (diffusible radical) and the surface or crypts of the apoprotein bind to the xenobiotic substrate (and in some cases, the heme distal pocket could also serve the latter function). Thus, CYPs enhance reaction rates and selectivity/specificity via a hitherto unrecognized modality.

Pharmacogenomics of the cytochrome P450 2C family: impacts of amino acid variations on drug metabolism

Pharmacogenomics investigates DNA and RNA variations in the human genome related to drug responses. Cytochrome P450 (CYP) is a supergene family of drug-metabolizing enzymes responsible for the metabolism of approximately 90% of human drugs. Among the major CYP isoforms, the CYP2C subfamily is of clinical significance because it metabolizes approximately 20% of clinically administrated drugs and represents several variant alleles leading to adverse drug reactions or altering drug efficacy. Here, we review recent progress on understanding the interindividual variability of the CYP2C members and the functional and clinical impact on drug metabolism. We summarize current advances in the molecular modeling of CYP2C polymorphisms and discuss the structural bases and molecular mechanisms of amino acid variants of CYP2C members that affect drug metabolism.

Novel free fatty acid receptor 1 (GPR40) agonists based on 1,3,4-thiadiazole-2-carboxamide scaffold

Free fatty acid receptor 1 (FFA1), previously known as GPR40 is a G protein-coupled receptor and a new target for treatment of type 2 diabetes. Two series of FFA1 agonists utilizing a 1,3,4-thiadiazole-2-caboxamide scaffold were synthetized. Both series offered significant improvement of the potency compared to the previously described 1,3,4-thiadiazole-based FFA1 agonists and high selectivity for FFA1. Molecular docking predicts new aromatic interactions with the receptor that improve agonist potency. The most potent compounds from both series were profiled for in vitro ADME properties (plasma and metabolic stability, LogD, plasma protein binding, hERG binding and CYP inhibition). One series suffered very rapid degradation in plasma and in presence of mouse liver microsomes. However, the other series delivered a lead compound that displayed a reasonable ADME profile together with the improved FFA1 potency.

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.