Exploring QSAR for CYP11B2 binding affinity and CYP11B2/CYP11B1 selectivity of diverse functional compounds using GFA and G/PLS techniques

In silico selectivity modeling of pyridine and pyrimidine based CYP11B1 and CYP11B2 inhibitors: A case study

DESIGN

of selective drug candidates for highly structural similar targets is a challenging task for researchers. The main objective of this study was to explore the selectivity modeling of pyridine and pyrimidine scaffold towards the highly homologous targets CYP11B1 and CYP11B2 enzymes by in silico (Molecular docking and QSAR) approaches. In this regard, a big dataset (n = 228) of CYP11B1 and CYP11B2 inhibitors were gathered and classified based on heterocyclic ring and the exhaustive analysis was carried out for pyridine and pyrimidinescaffolds. The LibDock algorithm was used to explore the binding pattern, screening, and identify the structural feature responsible for the selectivity of the ligands towards the studied targets. Finally, QSAR analysis was done to explore the correlation between various binding parameters and structural features responsible for the inhibitory activity and selectivity of the ligands in a quantitative way. The docking and QSAR analysis clearly revealed and distinguished the importance of structural features, functional groups attached for CYP11B2 and CYP11B1 selectivity for pyridine and pyrimidine analogs. Additionally, the docking analysis highlighted the differentiating amino acids residues for selectivity for ligands for each of the enzymes. The results obtained from this research work will be helpful in designing the selective CYP11B1/CYP11B2 inhibitors.

Combined pharmacophore modeling, 3D-QSAR, homology modeling and docking studies on CYP11B1 inhibitors

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.