Development of a Physiologically Based Model of Bilirubin Metabolism in Health and Disease and Its Comparison With Real‐World Data

Bilirubin is a breakdown product of erythrocytes and plays a crucial role in elimination of heme-containing proteins. After its synthesis in the reticuloendothelial system, unconjugated bilirubin is released into plasma and taken up into the liver. In hepatocytes, bilirubin is conjugated and excreted into the gastrointestinal tract via bile, where it is further converted to urobilinoids. There are various genetic factors causing abnormal bilirubin levels in plasma, such as Gilbert syndrome, Crigler-Najjar syndrome, Dubin-Johnson syndrome, and Rotor syndrome. To better understand bilirubin metabolism and its disorders, this study develops a physiologically based computational model incorporating published literature as well as real-world clinical data from the Explorys database. The model simulates bilirubin levels in both healthy individuals and patients with disorders of bilirubin metabolism. Population simulations show that Gilbert syndrome requires a substantial reduction in UDP-glucuronosyltransferase 1A1 activity, while Crigler-Najjar syndrome requires near-complete loss of its function. In contrast, Dubin-Johnson syndrome is characterized by a significant impairment of multidrug resistance-associated protein 2 activity. To also illustrate model behavior under targeted perturbations, we simulated administration of atazanavir in healthy individuals and patients with Gilbert syndrome to investigate its effect on bilirubin levels. Relative to baseline, unconjugated bilirubin maximum concentration (C_max) increased by 34% in healthy individuals but by 67% in Gilbert syndrome. Overall, this study provides a conceptual and mechanistically informed framework for studying bilirubin homeostasis and the functional consequences of drug administration in health and disease.