This study aimed to construct a physiologically based pharmacokinetic (PBPK) model of rifampicin that can accurately and quantitatively predict complex drug-drug interactions (DDIs) involving its saturable hepatic uptake and auto-induction. Using in silico and in vitro parameters, and reported clinical pharmacokinetic data, rifampicin PBPK model was built and relevant parameters for saturable hepatic uptake and UDP-glucuronosyltransferase auto-induction were optimized by fitting. The parameters for cytochrome P450 (CYP) 3A and CYP2C9 induction by rifampicin were similarly optimized using clinical DDI data with midazolam and tolbutamide as probe substrates, respectively. For validation, our current PBPK model was applied to simulate complex DDIs with glibenclamide (a substrate of CYP3A/2C9 and hepatic organic anion transporting polypeptide (OATP)). Simulated results were in quite good accordance with the observed data. Altogether, our constructed PBPK model of rifampicin demonstrates the robustness and utility in quantitatively predicting CYP3A/2C9 induction- and/or OATP inhibition-mediated DDIs with victim drugs. This article is protected by copyright. All rights reserved.