Classification analysis of P-glycoprotein substrate specificity

Prediction of P-glycoprotein substrate specificity (S-PGP) can be viewed as a constituent part of a compound's pharmaceutical profiling in drug design. This task is difficult to achieve due to several factors that raised many contradictory opinions: (i) the disparity between the S-PGP values obtained in different assays, (ii) the confusion between Pgp substrates and inhibitors, (iii) the confusion between lipophilicity and amphiphilicity of Pgp substrates, and (iv) the dilemma of describing class-specific relationships when Pgp has no binding sites of high ligand specificity. In this work, we compiled S-PGP data for 1000 compounds. All data were represented in a binary format, assigning S-PGP =1 for substrates and S-PGP =0 for non-substrates. Each value was ranked according to the reliability of experimental assay. Two data sets were considered. Set 1 included 220 compounds with S-PGP from polarized transport across MDR1 transfected cell monolayers. Set 2 included the entire list of 1000 compounds, with S-PGP values of generally lower reliability. Both sets were analysed using a stepwise classification structure-activity relationship (C-SAR) method, leading to derivation of simple rules for crude estimation of S-PGP values. The obtained rules are based on the following factors: (i) compound's size expressed through molar weight or volume, (ii) H-accepting given by the Abraham's beta (that can be crudely approximated by the sum of O and N atoms), and (iii) ionization given by the acid and base pK(a) values. Very roughly, S-PGP can be estimated by the rule of fours. Compounds with (N+O)greater than or equal to8, MW>400 and acid pK(a) > 4 are likely to be Pgp substrates, whereas compounds with (N+O)less than or equal to4, MW<400 and base pK(a) <8 are likely to be non-substrates. The obtained results support the view that Pgp functioning can be compared to a complex mini-pharmacokinetic system with fuzzy specificity. This system can be described by a probabilistic version of Abraham's solvation equation, suggesting a certain similarity between Pgp transport and chromatographic retention. The chromatographic model does not work in the case of marginal compounds with properties close to the global physicochemical cut-offs. In the latter case various class-specific rules must be considered. These can be associated with the amphiphilicity and biological similarity of compounds. The definition of class-specific effects entails construction of the knowledge base that can be very useful in ADME profiling of new drugs.