Abstract

Baicalein, a typical flavonoid presented in Scutellariae radix, exhibits a unique metabolic profile during first-pass metabolism: parallel glucuronidation and sulfation pathways, with possible substrate inhibition in both pathways. In this project, we aimed to construct an agent-based model to dynamically represent baicalein pharmacokinetics and to verify the substrate inhibition hypothesis. The model consisted of three 3D spaces and two membranes: apical space (S1), intracellular space (S2), basolateral space (S3), apical membrane (M1), and basolateral membrane (M2). In silico enzymes (UDP-glucuronosyltransferases (UGTs) and sulfotransferases (SULTs)) and binder components were placed in S2. The model was then executed to simulate one-pass metabolism experiments of baicalein. With the implementation of a two-site enzyme design, the simulated results captured the preset qualitative and quantitative features of the wet-lab observations. The feasible parameter set showed that substrate inhibition happened in both conjugation pathways of baicalein. The simulation results suggested that the sulfation pathway was dominant at low concentrations and that SULT was more inclined to substrate inhibition than UGT. Cross-model validation was satisfactory. Our findings were consistent with a previously reported catenary model. We conclude that the mechanisms represented by our model are plausible. Our novel modeling approach could dynamically represent the metabolic pathways of baicalein in a Caco-2 system.