Combined Pharmacophore Modeling, 3D-QSAR, Homology Modeling and Docking Studies on CYP11B1 Inhibitors

Molecules

Rui Yu¹  Juan Wang³  Rui Wang¹  Yong Lin²  Yong Hu¹  Yuanqiang Wang¹  Mao Shu¹  Zhihua Lin¹ 

[1] School of Pharmacy and Bioengineering, Chongqing University of Technology, Chongqing 400054, China; E-Mails:;School of Chemistry and Chemical Engineering, Chongqing University of Technology, Chongqing 400054, China; E-Mail:;College of Bioengineering, Chongqing University, Chongqing 400044, China; E-Mail:

关键词: CYP11B1 inhibitors;    3D-QSAR;    pharmacophore model;    Cushing’s syndrome;   

DOI  :  10.3390/molecules20011014

来源: mdpi

【 摘 要 】

The mitochondrial cytochrome P450 enzymes inhibitor steroid 11β-hydroxylase (CYP11B1) can decrease the production of cortisol. Therefore, these inhibitors have an effect in the treatment of Cushing’s syndrome. A pharmacophore model generated by Genetic Algorithm with Linear Assignment for Hypermolecular Alignment of Datasets (GALAHAD) was used to align the compounds and perform comparative molecular field analysis (CoMFA) with Q² = 0.658, R² = 0.959. The pharmacophore model contained six hydrophobic regions and one acceptor atom, and electropositive and bulky substituents would be tolerated at the A and B sites, respectively. A three-dimensional quantitative structure-activity relationship (3D-QSAR) study based on the alignment with the atom root mean square (RMS) was applied using comparative molecular field analysis (CoMFA) with Q² = 0.666, R² = 0.978, and comparative molecular similarity indices analysis (CoMSIA) with Q² = 0.721, R² = 0.972. These results proved that all the models have good predictability of the bioactivities of inhibitors. Furthermore, the QSAR models indicated that a hydrogen bond acceptor substituent would be disfavored at the A and B groups, while hydrophobic groups would be favored at the B site. The three-dimensional (3D) model of the CYP11B1 was generated based on the crystal structure of the CYP11B2 (PDB code 4DVQ). In order to probe the ligand-binding modes, Surflex-dock was employed to dock CYP11B1 inhibitory compounds into the active site of the receptor. The docking result showed that the imidazolidine ring of CYP11B1 inhibitors form H bonds with the amino group of residue Arg155 and Arg519, which suggested that an electronegative substituent at these positions could enhance the activities of compounds. All the models generated by GALAHAD QSAR and Docking methods provide guidance about how to design novel and potential drugs for Cushing’s syndrome treatment.

【 授权许可】

CC BY   
© 2015 by the authors; licensee MDPI, Basel, Switzerland.