CYP Inhibition
Inhibition of CYPs is important for understanding drug toxicity and drug-drug interactions. The Metabolism Module includes generic inhibition models for nine recombinant CYPs: 1A2, 2A6, 2B6, 2C8, 2C9, 2C19, 2D6, 2E1, and 3A4.
The structures and biological activity data of corresponding CYP isoforms were extracted from the National Center for Biotechnology Information (NCBI) PubChem database AID 1851. According to the Pubchem Bioassay database, a compound was classified as a CYP Inhibitor if the AC50 value (compound concentration that leads to 50% of the activity of an inhibition control) is ≤10 μM and as a Non-inhibitor if the AC50 value is >57 μM. If the AC50 value falls between 10 and 57 μM, the compound is defined as inconclusive. Only the compounds defined as actives and inactives were considered for further analysis. The inconclusive compounds were excluded to prevent any uncertainty during model building.
These datasets were appended to our in-house collection of CYP inhibition data. All the datasets that were curated for regression models were used for building the classification models. Compounds with Ki less than 10 μM were labeled as “Inhibitors” and those with Ki greater than 30 μM were labeled as “Non-inhibitors”.
The combined datasets were pre-processed to eliminate entries containing non-covalent complexes, inorganic compounds and mixtures. Moreover, the entries containing salts were converted to their corresponding acids or bases, and water molecules were removed from hydrates. The duplicate entries were removed using the following criteria. If the duplicates had identical endpoints, then one entry was retained, and others eliminated.
For entries with conflicting end-points:
-
If one of the duplicate entries contained an in-house (hand-curated) entry, the in-house data point was retained while that from the PubChem dataset was eliminated.
-
If the duplicate entries originated from the PubChem dataset, both entries (with conflicting end-points) were removed from the final dataset.
The table below lists the statistics for each model. Note that concordance, sensitivity and specificity are above 80% for the training and test sets of all models. The CYP inhibition classification models are used as a filter for the CYP Ki models: i.e., if the classification model indicates that the compound does not inhibit that isoform, then the Ki model's prediction will be "NonInhibitor".
|
Model |
Dataset |
Negatives |
Positives |
Total |
Correct |
Concordance |
Sensitivity |
Specificity |
|---|---|---|---|---|---|---|---|---|
|
CYP1A2_Inh |
Training |
5445 |
4999 |
10444 |
9150 |
87.6% |
87.5% |
87.7% |
|
Test |
1386 |
1226 |
2612 |
2237 |
85.6% |
86.1% |
85.3% |
|
|
CYP2A6_Inh |
Training |
776 |
215 |
991 |
939 |
94.8% |
94.0% |
95.0% |
|
Test |
194 |
54 |
248 |
233 |
94.0% |
92.6% |
94.3% |
|
|
CYP2B6_Inh |
Training |
267 |
159 |
426 |
396 |
93.0% |
90.6% |
94.4% |
|
Test |
67 |
40 |
107 |
94 |
87.9% |
87.5% |
88.1% |
|
|
CYP2C8_Inh |
Training |
202 |
450 |
652 |
597 |
91.6% |
92.0% |
90.6% |
|
Test |
45 |
119 |
164 |
146 |
89.0% |
89.9% |
86.7% |
|
|
CYP2C9_Inh |
Training |
6739 |
3712 |
10451 |
8750 |
83.7% |
81.0% |
85.2% |
|
Test |
1705 |
908 |
2613 |
2175 |
83.2% |
80.3% |
84.8% |
|
|
CYP2C19_Inh |
Training |
6127 |
4881 |
11008 |
9235 |
83.9% |
83.5% |
84.2% |
|
Test |
1584 |
1168 |
2752 |
2283 |
83.0% |
82.0% |
83.6% |
|
|
CYP2D6_Inh |
Training |
8810 |
2996 |
11806 |
9978 |
84.5% |
81.3% |
85.6% |
|
Test |
2171 |
781 |
2952 |
2435 |
82.5% |
80.0% |
83.4% |
|
|
CYP2E1_Inh |
Training |
305 |
93 |
398 |
349 |
87.7% |
88.2% |
87.5% |
|
Test |
69 |
31 |
100 |
87 |
87.0% |
87.1% |
87.0% |
|
|
CYP3A4_Inh |
Training |
6225 |
5263 |
11488 |
9468 |
82.4% |
82.4% |
82.4% |
|
Test |
1569 |
1304 |
2873 |
2323 |
80.9% |
80.1% |
81.5% |