The beta-blocker metoprolol (the sixth most commonly prescribed drug in the USA in 2017) is subject to considerable drug-gene interaction (DGI) effects caused by genetic variations of the CYP2D6 gene. CYP2D6 poor metabolizers (5.7% of US population) show approximately five-fold higher metoprolol exposure compared to CYP2D6 normal metabolizers. This study aimed to develop a whole-body physiologically based pharmacokinetic (PBPK) model to predict CYP2D6 DGIs with metoprolol. The metoprolol (R)- and (S)-enantiomers as well as the active metabolite α-hydroxymetoprolol were implemented as model compounds, employing data of 48 different clinical studies (dosing range 5-200 mg). To mechanistically describe the effect of CYP2D6 polymorphisms, two separate metabolic CYP2D6 pathways (α-hydroxylation and O-demethylation) were incorporated for both metoprolol enantiomers. The good model performance is demonstrated in predicted plasma concentration-time profiles compared to observed data, goodness-of-fit plots, and low geometric mean fold errors of the predicted AUClast (1.27) and Cmax values (1.23) over all studies. For DGI predictions, 18 out of 18 DGI AUClast ratios and 18 out of 18 DGI Cmax ratios were within two-fold of the observed ratios. The newly developed and carefully validated model was applied to calculate dose recommendations for CYP2D6 polymorphic patients and will be freely available in the Open Systems Pharmacology repository.
Keywords: cytochrome P450 2D6 (CYP2D6); dose adaptation; drug-gene interactions (DGIs); metoprolol; metoprolol enantiomers; model-informed precision dosing; physiologically based pharmacokinetic (PBPK) modeling; α-hydroxymetoprolol.