Improved homology model of the human all-trans retinoic acid metabolizing enzyme CYP26A1

A new CYP26A1 homology model was built based on the crystal structure of cyanobacterial CYP120A1. The model quality was examined for stereochemical accuracy, folding reliability, and absolute quality using a variety of different bioinformatics tools. Furthermore, the docking capabilities of the model were assessed by docking of the natural substrate all-trans-retinoic acid (atRA), and a group of known azole- and tetralone-based CYP26A1 inhibitors. The preferred binding pose of atRA suggests the (4S)-OH-atRA metabolite production, in agreement with recently available experimental data. The distances between the ligands and the heme group iron of the enzyme are in agreement with corresponding distances obtained for substrates and azole inhibitors for other cytochrome systems. The calculated theoretical binding energies agree with recently reported experimental data and show that the model is capable of discriminating between natural substrate, strong inhibitors (R116010 and R115866), and weak inhibitors (liarozole, fluconazole, tetralone derivatives).