A physiologically-based flow network model for hepatic drug elimination III: 2D/3D DLA lobule models

BACKGROUND: One of the major issues in current pharmaceutical development is potential hepatotoxicity and drug-induced liver damage. This is due to the unique metabolic processes performed in the liver to prevent accumulation of a wide range of chemicals in the blood. Recently, we developed a physiologically-based lattice model to address the transport and metabolism of drugs in the liver lobule (liver functional unit). METHOD: In this paper, we extend our idealized model to consider structural and spatial variability in two and three dimensions. We introduce a hexagonal-based model with one input (portal vein) and six outputs (hepatic veins) to represent a typical liver lobule. To capture even more realistic structures, we implement a novel sequential diffusion-limited aggregation (DLA) method to construct a morphological sinusoid network in the lobule. A 3D model constructed with stacks of multiple 2D sinusoid realizations is explored to study the effects of 3D structural variations. The role of liver zonation on drug metabolism in the lobule is also addressed, based on flow-based predicted steady-state O(2) profiles used as a zonation indicator. RESULTS: With this model, we analyze predicted drug concentration levels observed exiting the lobule with their detailed distribution inside the lobule, and compare with our earlier idealized models. In 2D, due to randomness of the sinusoidal structure, individual hepatic veins respond differently (i.e. at different times) to injected drug. In 3D, however, the variation of response to the injected drug is observed to be less extreme. Also, the production curves show more diffusive behavior in 3D than in 2D. CONCLUSION: Although, the individual producing ports respond differently, the average lobule production summed over all hepatic veins is more diffuse. Thus the net effect of all these variations makes the overall response smoother. We also show that, in 3D, the effect of zonation on drug production characteristics appears quite small. Our new biophysical structural analysis of a physiologically-based 3D lobule can therefore form the basis for a quantitative assessment of liver function and performance both in health and disease ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12976-016-0034-5) contains supplementary material, which is available to authorized users.