Homology-Based Variant-Effect Predictors Break Down on Cytochrome P450 Pharmacogenes

Abstract

Cytochrome P450 (CYP) enzymes metabolize roughly three-quarters of clinically used drugs; genetic variation in these enzymes is a leading source of interindividual differences in drug response. Predicting a variant’s functional effect is therefore critical, yet the consequences of most CYP variants remain unknown. Many state-of-the-art variant-effect predictors rest on a homology-based paradigm that scores variants by evolutionary conservation — an assumption that pharmacogenes including CYPs violate. Indeed, focusing on human CYPs, we show that homology-based models fail systematically. AlphaMissense (AM) assigns variants to its “ambiguous” class at nearly twice the proteome-wide rate across six CYPs, and within that class the scores are essentially uncorrelated with CYP2C9 DMS activity (Spearman’s ρ = 0.069); Evolutionary Scale Modeling 2 (ESM-2) shows the same pattern. We further hypothesized that non-homology-based features (sequence position, substitution chemistry, binding-site distance, and secondary structure) might help resolve the ambiguous calls, but their explanatory power is weak. Using CYP2C9 DMS activity as ground truth, we built a k-nearest-neighbors model over ESM-2 embeddings and ensembled it with AM and ESM-2 masked marginal probability, improving the ambiguous-class correlation roughly ten-fold, from ρ = 0.069 to 0.715 (overall ρ from 0.638 to 0.825). However, this markedly improved accuracy does not translate into agreement with clinical annotations. Drawing on evidence that a variant’s effect can depend on the drug, we hypothesize that substrate identity is the key missing feature in current models, and that predicting function for these multi-substrate enzymes may require redefining function as substrate-conditioned.