The cytochrome P450 (Cyp) proteins Cyp1A1 and Cyp1A2 are strongly induced in the mouse liver by the potent environmental toxicant 2, 3, 7, 8‐tetrachlorodibenzo‐p‐dioxin (TCDD), acting through the aryl hydrocarbon receptor (AHR). The induction of Cyp1A1 is localized within the centrilobular regions of the mouse liver at low doses of TCDD, progressing to pan‐lobular induction at higher doses. Even without chemical perturbation, metabolic functions and associated genes are basally zonated in the liver lobule along the central‐to‐portal axis. To investigate the mechanistic basis of spatially restricted gene induction by TCDD, we have developed a multiscale computational model of the mouse liver lobule with single‐cell resolution. The spatial location of individual hepatocytes in the model was calibrated from previously published high‐resolution images. A systems biology model of the network of biochemical signaling pathways underlying Cyp1A1 and Cyp1A2 induction was then incorporated into each hepatocyte in the model. Model simulations showed that a negative feedback loop formed by binding of the induced Cyp1A2 protein to TCDD, together with cooperative gene induction by the β‐catenin/AHR/TCDD transcription factor complex and β‐catenin, help produce the spatially localized induction pattern of Cyp1A1. Although endogenous WNT regulates the metabolic zonation of many genes, it was not a driver of zonal Cyp1A1 induction in our model. Conclusion: In this work, we used data‐driven computational modeling to identify the mechanistic basis of zonally restricted gene expression induced by the potent and persistent environmental pollutant TCDD. The multiscale model and derived results clarify the mechanisms of dose‐dependent hepatic gene induction responses to TCDD. Additionally, this work contributes to our broader understanding of spatial gene regulation along the liver lobule.
Read full abstract