Three and four dimensional-quantitative structure activity relationship (3D/4D-QSAR) analyses of CYP2D6 inhibitors

Three- and four-dimensional quantitative structure activity relationship (3D/4D-QSAR) pharmacophore models of competitive inhibitors of CYP2D6 were constructed using data from our laboratory or the literature. The 3D-QSAR pharmacophore models of the common structural features of CYP2D6 inhibitors were built using the program Catalyst (Molecular Simulations, San Diego, CA, USA). These 3D-QSAR models were compared with 3D and 4D-QSAR partial least squares (PLS) models which were constructed using molecular surface-weighted holistic invariant molecular (MS-WHIM) descriptors of size and shape of inhibitors. The first Catalyst model was generated from multiple conformers of competitive inhibitors (n = 20) of CYP2D6 mediated bufurolol 1'-hydroxylation. This model demonstrated a correlation of observed and predicted Ki (apparent) values of r = 0.75. A second Catalyst model was constructed from literature derived Ki (apparent) values (n = 31) for the inhibition of CYP2D6. This model provided a correlation of observed and predicted inhibition for CYP2D6 of r = 0.91. Both Catalyst Ki pharmacophores were then validated by predicting the Ki (apparent) of a test set of known CYP2D6 inhibitors (n = 15). Ten out of 15 of these Ki (apparent) values were predicted to be within one log residual of the observed value using our CYP2D6 inhibitor model, while the literature model predicted nine out of 15 values. Similarly, 3D- and 4D-QSARs derived from PLS MS-WHIM for our dataset yielded predictable models as assessed using cross-validation. The corresponding cross-validated PLS MS-WHIM model for the literature dataset yielded a comparable 3D-QSAR and improved 4D-QSAR value. Such computational models will aid in future prediction of drug-drug interactions.