Prediction of phenotypic consequences of mutations constitutes an important aspect of precision medicine. Current computational tools mostly rely on evolutionary conservation and have been calibrated on variants associated with disease, which poses conceptual problems for assessment of variants in poorly conserved pharmacogenes. Here, we evaluated the performance of 18 current functionality prediction methods leveraging experimental high-quality activity data from 337 variants in genes involved in drug metabolism and transport and found that these models only achieved probabilities of 0.1–50.6% to make informed conclusions. We therefore developed a functionality prediction framework optimized for pharmacogenetic assessments that significantly outperformed current algorithms. Our model achieved 93% for both sensitivity and specificity for both loss-of-function and functionally neutral variants, and we confirmed its superior performance using cross validation analyses. This novel model holds promise to improve the translation of personal genetic information into biological conclusions and pharmacogenetic recommendations, thereby facilitating the implementation of Next-Generation Sequencing data into clinical diagnostics.

Cytochrome P450 2D6 (CYP2D6) is an enzyme of potential importance for the metabolism of drugs used clinically, and it exhibits genetic polymorphism with interindividual differences in metabolic activity. To date, 21 CYP2D6 allelic variants have been identified in the Japanese population. The aim of this study was to investigate the functional characterization of CYP2D6 variants identified in Japanese subjects. Wild-type CYP2D6 and its variants, namely, CYP2D6.2, CYP2D6.10, CYP2D6.14A, CYP2D6.14B, CYP2D6.18, CYP2D6.27, CYP2D6.36, CYP2D6.39, CYP2D6.47, CYP2D6.48, CYP2D6.49, CYP2D6.50, CYP2D6.51, CYP2D6.53, CYP2D6.54, CYP2D6.55, and CYP2D6.57 were transiently expressed in COS-7 cells, and enzymatic activities of the CYP2D6 variant proteins were characterized using bufuralol and dextromethorphan. Functional characterization of 17 CYP2D6 variants revealed an absence of enzyme activity in four (CYP2D6.14A, CYP2D6.36, CYP2D6.47, and CYP2D6.57), low activity in eight (CYP2D6.10, CYP2D6.14B, CYP2D6.18, CYP2D6.49, CYP2D6.50, CYP2D6.51, CYP2D6.54, and CYP2D6.55), and high activity in one (CYP2D6.53) compared with the wild type. Analysis of CYP2D6 variant proteins can be useful for predicting CYP2D6 phenotypes and could be applied to personalized drug therapy.

The cytochrome P450 (CYP) superfamily is one of the most important groups of enzymes involved in drug metabolism. It is responsible for the metabolism of a large number of drugs. Many CYP isoforms are expressed polymorphically, and catalytic alterations of allelic variant proteins can affect the metabolic activities of many drugs. The CYP2D6, CYP2C9, CYP2C19, and CYP2B6 genes are particularly polymorphic, whereas CYP1A1, CYP1A2, CYP2E1, and CYP3A4 are relatively well conserved without common functional polymorphisms. In vitro studies using cDNA expression systems are useful tools for evaluating functional alterations of the allelic variants of CYP, particularly for low-frequency alleles. Recombinant CYPs have been successfully expressed in bacteria, yeast, baculoviruses, and several mammalian cells. Determination of CYP variant-mediated kinetic parameters (Km and Vmax) in vitro can be useful for predicting drug dosing and clearance in humans. This review focuses on the advantages and disadvantages of the various cDNA-expression systems used to determine the kinetic parameters for CYP allelic variants, the methods for determining the kinetic parameters, and the findings of in vitro studies on highly polymorphic CYPs, including CYP2D6, CYP2C9, CYP2C19, and CYP2B6.