The effects of nitrogen-heme-iron coordination on substrate affinities for cytochrome P450 2E1

A descriptor based computational model was developed for cytochrome P450 2E1 (CYP2E1) based on inhibition constants determined for inhibition of chlorzoxazone, or 4-nitrophenol, metabolism. An empirical descriptor for type II binding was developed and tested for a series of CYP2E1 inhibitors. Inhibition constants where measured for 51 different compounds. A fast 2-dimensional predictive model was developed based on 40 compounds, and tested on 8 compounds of diverse structure. The trained model (n = 40) had an r² value of 0.76 and an RMSE of 0.48. The correlation between the predicted and actual pK_i values of the test set of compounds not included in the model gives an r² value of 0.78. The features that described binding include heme coordination (type II binding), molecular volume, octanol/water partition coefficient, solvent accessible surface area, and the sum of the atomic polarizabilities. The heme coordination parameter assigns an integer between 0 and 6 depending on structure, and is a new descriptor, based on simple quantum chemical calculations with correction for steric effects. The type II binding parameter was found to be important in obtaining a good correlation between predicted and experimental inhibition constants increasing the r ² value from 0.38 to 0.77.