Abstract:The highly polymorphic human cytochrome P450 2D6 (CYP2D6) metabolizes about 25% of currently used drugs. In this study, we have explored the interaction of a large number of substrates (n = 120) with wild-type and mutated CYP2D6 by molecular docking using the CDOCKER module. Before we conducted the molecular docking and virtual mutations, the pharmacophore and QSAR models of CYP2D6 substrates were developed and validated. Finally, we explored the interaction of a traditional Chinese herbal formula, Fangjifuling decoction, with CYP2D6 by virtual screening. The optimized pharmacophore model derived from 20 substrates of CYP2D6 contained two hydrophobic features and one hydrogen bond acceptor feature, giving a relevance ratio of 76% when a validation set of substrates were tested. However, our QSAR models gave poor prediction of the binding affinity of substrates. Our docking study demonstrated that 117 out of 120 substrates could be docked into the active site of CYP2D6. Forty one out of 117 substrates (35.04%) formed hydrogen bonds with various active site residues of CYP2D6 and 53 (45.30%) substrates formed a strong pi-pi interaction with Phe120 (53/54), with only carvedilol showing pi-pi interaction with Phe483. The active site residues involving hydrogen bond formation with substrates included Leu213, Lys214, Glu216, Ser217, Gln244, Asp301, Ser304, Ala305, Phe483, and Phe484. Furthermore, the CDOCKER algorithm was further applied to study the impact of mutations of 28 active site residues (mostly non-conserved) of CYP2D6 on substrate binding modes using five probe substrates including bufuralol, debrisoquine, dextromethorphan, sparteine, and tramadol. All mutations of the residues examined altered the hydrogen bond formation and/or aromatic interactions, depending on the probe used in molecular docking. Apparent changes of the binding modes have been observed with the Glu216Asp and Asp301Glu mutants. Overall, 60 compounds out of 130 from Fangjifuling decoction matched our pharmacophore model for CYP2D6 substrates. Fifty four out of these 60 compounds could be docked into the active site of CYP2D6 and 24 of 54 compounds formed hydrogen bonds with Glu216, Asp301, Ser304, and Ala305 in CYP2D6. These results have provided further insights into the factors that determining the binding modes of substrates to CYP2D6. Screening of high-affinity ligands for CYP2D6 from herbal formula using computational models is a useful approach to identify potential herb-drug interactions.

Exploring the Chemical Space of CYP17A1 Inhibitors Using Cheminformatics and Machine Learning

QSAR Research of Novel Tetrandrine Derivatives Against Human Hepatocellular Carcinoma.

Three dimensional quantitative structure activity relationship of P45017α inhibitors of 17-substituted steroids

Molecular dynamic investigations of the mutational effects on structural characteristics and tunnel geometry in CYP17A1.

Classification of Cytochrome P450 Inhibitors and Noninhibitors Using Combined Classifiers.

Structural Features and Dynamic Investigations of the Membrane-Bound Cytochrome P450 17A1

Three-Dimensional Quantitative Structure Activity Relationship for Cyp2d6 Substrates

Multiple machine learning, molecular docking, and ADMET screening approach for identification of selective inhibitors of CYP1B1

3D structure modeling of cytochrome P450 2C19 and its implication for personalized drug design.

3d Structure Modeling Of Cytochrome P4502c19 And Its Implication For Personalized Drug Design

Discovery of novel 1,2,3,4-tetrahydrobenzo[4, 5]thieno[2, 3-c]pyridine derivatives as potent and selective CYP17 inhibitors.

Prediction of Cytochrome P450 Inhibition Using a Deep Learning Approach and Substructure Pattern Recognition

Pharmacophore, Qsar, and Binding Mode Studies of Substrates of Human Cytochrome P450 2d6 (cyp2d6) Using Molecular Docking and Virtual Mutations and an Application to Chinese Herbal Medicine Screening

3D QSAR Studies, Pharmacophore Modeling and Virtual Screening on a Series of Steroidal Aromatase Inhibitors.

Combined Pharmacophore Modeling, 3D-Qsar, Homology Modeling and Docking Studies on Cyp11b1 Inhibitors

Comparative Dynamics and Functional Mechanisms of the CYP17A1 Tunnels Regulated by Ligand Binding.

An improved large-scale prediction model of CYP1A2 inhibitors by using combined fragment descriptors

Finding the Molecular Scaffold of Nuclear Receptor Inhibitors Through High-Throughput Screening Based on Proteochemometric Modelling

In Silico Prediction of Chemicals Binding to Aromatase with Machine Learning Methods.

Prediction of cytochrome P4502C9 inhibitors based on molecular indices

Initial Development of Automated Machine Learning-Assisted Prediction Tools for Aryl Hydrocarbon Receptor Activators