Understanding liver exposure of hepatic transporter substrates in clinical studies is often critical, as it typically governs pharmacodynamics, drug-drug interactions, and toxicity for certain drugs. However, this is a challenging task since there is currently no easy method to directly measure drug concentration in the human liver. Using bosentan as an example, we demonstrate a new approach to estimate liver exposure based on observed systemic pharmacokinetics from clinical studies using physiologically based pharmacokinetic modeling. The prediction was verified to be both accurate and precise using sensitivity analysis. For bosentan, the predicted pseudo steady-state unbound liver-to-unbound systemic plasma concentration ratio was 34.9 (95% confidence interval: 4.2, 50). Drug-drug interaction (i.e., CYP3A and CYP2B6 induction) and inhibition of hepatic transporters (i.e., bile salt export pump, multidrug resistance-associated proteins, and sodium-taurocholate cotransporting polypeptide) were predicted based on the estimated unbound liver tissue or plasma concentrations. With further validation and refinement, we conclude that this approach may serve to predict human liver exposure and complement other methods involving tissue biopsy and imaging.