Mechanistic Study of Icaritin-Induced Inactivation of Cytochrome P450 2C9.

Icaritin (ICT) is a prenylflavonoid derivative that has been approved by National Medical Products Administration for the treatment of hepatocellular carcinoma. This study aims to evaluate the potential inhibitory effect of ICT against cytochrome P450 (CYP) enzymes and to elucidate the inactivation mechanisms. Results showed that ICT inactivated CYP2C9 in a time-, concentration-, and NADPH-dependent manner with Ki = 1.896 µM, Kinact = 0.02298 minutes-1, and Kinact/Ki = 12 minutes-1 mM-1, whereas the activities of other CYP isozymes was minimally affected. Additionally, the presence of CYP2C9 competitive inhibitor, sulfaphenazole, superoxide dismutase/catalase system, and GSH all protected CYP2C9 from ICT-induced activity loss. Moreover, the activity loss was neither recovered by washing the ICT-CYP2C9 preincubation mixture nor the addition of potassium ferricyanide. These results, collectively, implied the underlying inactivation mechanism involved the covalent binding of ICT to the apoprotein and/or the prosthetic heme of CYP2C9. Furthermore, an ICT-quinone methide (QM)-derived GSH adduct was identified, and human glutathione S-transferases (GST) isozymes GSTA1-1, GSTM1-1, and GSTP1-1 were shown to be substantially involved in the detoxification of ICT-QM. Interestingly, our systematic molecular modeling work predicted that ICT-QM was covalently bound to C216, a cysteine residue located in the F-G loop downstream of substrate recognition site (SRS) 2 in CYP2C9. The sequential molecular dynamics simulation confirmed the binding to C216 induced a conformational change in the active catalytic center of CYP2C9. Lastly, the potential risks of clinical drug-drug interactions triggered by ICT as a perpetrator were extrapolated. In summary, this work confirmed that ICT was an inactivator of CYP2C9. SIGNIFICANCE STATEMENT: This study is the first to report the time-dependent inhibition of CYP2C9 by icaritin (ICT) and the intrinsic molecular mechanism behind it. Experimental data indicated that the inactivation was via irreversible covalent binding of ICT-quinone methide to CYP2C9, while molecular modeling analysis provided additional evidence by predicting C216 as the key binding site which influenced the structural confirmation of CYP2C9's catalytic center. These findings suggest the potential of drug-drug interactions when ICT is co-administered with CYP2C9 substrates clinically.

An efficient multistep ligand-based virtual screening approach for GPR40 agonists.

G protein-coupled receptor 40/free fatty acid receptor 1 (GPR40/FFAR1) is a member of the GPCR superfamily, and GPR40 agonists have therapeutic potential for type 2 diabetes. With the crystal structure of GPR40 currently unavailable, various ligand-based virtual screening approaches can be applied to identify novel agonists of GPR40. It is known that each ligand-based method has its own advantages and limitations. To improve the efficiency of individual ligand-based methods, an efficient multistep ligand-based virtual screening approach is presented in this study, including the pharmacophore-based screening, physicochemical property filtering, protein-ligand interaction fingerprint similarity analysis, and 2D-fingerprint structural similarity search. A focused decoy library was generated and used to evaluate the efficiency of this virtual screening protocol. This multistep workflow not only significantly improved the hit rate compared with each individual ligand-based method, but also identified diverse known actives from decoys. This protocol may serve as an efficient virtual screening tool for the targets without crystal structures available to discover novel active compounds.

An integrated virtual screening approach for VEGFR-2 inhibitors.

In recent years, various virtual screening (VS) tools have been developed, and many successful screening campaigns have been showcased. However, whether by conventional molecular docking or pharmacophore screening, the selection of virtual hits is based on the ranking of compounds by scoring functions or fit values, which remains the bottleneck of VS due to insufficient accuracy. As the limitations of individual methods persist, a comprehensive comparison and integration of different methods may provide insights into selecting suitable methods for VS. Here, we evaluated the performance of molecular docking, fingerprint-based 2D similarity and multicomplex pharmacophore in an individual and a combined manner, through a retrospective VS study on VEGFR-2 inhibitors. An integrated two-layer workflow was developed and validated through VS of VEGFR-2 inhibitors against the DUD-E database, which demonstrated improved VS performance through a ligand-based method ECFP_4, followed by molecular docking, and then a strict multicomplex pharmacophore. Through a retrospective comparison with six published papers, this integrated approach outperformed 43 out of 45 methods, indicating a great effectiveness. This kind of integrated VS approach can be extended to other targets for the screening and discovery of inhibitors.

Fragment-based strategy for structural optimization in combination with 3D-QSAR.

Fragment-based drug design has emerged as an important methodology for lead discovery and drug design. Different with other studies focused on fragment library design and active fragment identification, a fragment-based strategy was developed in combination with three-dimensional quantitative structure-activity relationship (3D-QSAR) for structural optimization in this study. Based on a validated scaffold or fragment hit, a series of structural optimization was conducted to convert it to lead compounds, including 3D-QSAR modelling, active site analysis, fragment-based structural optimization and evaluation of new molecules. 3D-QSAR models and active site analysis provided sufficient information for confirming the SAR and pharmacophoric features for fragments. This strategy was evaluated through the structural optimization on a c-Met inhibitor scaffold 5H-benzo[4,5]cyclohepta[1,2-b]pyridin-5-one, which resulted in an c-Met inhibitor with high inhibitory activity. Our study suggested the effectiveness of this fragment-based strategy and the druggability of our newly explored active region. The reliability of this strategy indicated it could also be applied to facilitate lead optimization of other targets.

De novo design of N-(pyridin-4-ylmethyl)aniline derivatives as KDR inhibitors: 3D-QSAR, molecular fragment replacement, protein-ligand interaction fingerprint, and ADMET prediction.

Vascular endothelial growth factor (VEGF) and its receptor tyrosine kinase VEGFR-2 or kinase insert domain receptor (KDR) have been identified as promising targets for novel anticancer agents. To achieve new potent inhibitors of KDR, we conducted molecular fragment replacement (MFR) studies for the understanding of 3D-QSAR modeling and the docking investigation of arylphthalazines and 2-((1H-Azol-1-yl)methyl)-N-arylbenzamides-based KDR inhibitors. Two favorable 3D-QSAR models (CoMFA with q(2), 0.671; r(2), 0.969; CoMSIA with q(2), 0.608; r(2), 0.936) have been developed to predict the biological activity of new compounds. The new molecular database generated by MFR was virtually screened using Glide (docking) and further evaluated with CoMFA prediction, protein-ligand interaction fingerprint (PLIF) and ADMET analysis. 44 N-(pyridin-4-ylmethyl)aniline derivatives as novel potential KDR inhibitors were finally obtained. In this paper, the work flow developed could be applied to de novo drug design and virtual screening potential KDR inhibitors, and use hit compounds to further optimize and design new potential KDR inhibitors.