Model Prediction and Validation of an Order Mechanism Controlling the Spatiotemporal Phenotype of Early Hepatocellular Carcinoma

Recently, hepatocyte–sinusoid alignment (HSA) has been identified as a mechanism that supports the coordination of hepatocytes during liver regeneration to reestablish a functional micro-architecture (Hoehme et al. in Proc Natl Acad Sci 107(23):10371–10376, 2010). HSA means that hepatocytes preferentially align along the closest micro-vessels. Here, we studied whether this mechanism is still active in early hepatocellular tumors. The same agent-based spatiotemporal model that previously correctly predicted HSA in liver regeneration was further developed to simulate scenarios in early tumor development, when individual initiated hepatocytes gain increased proliferation capacity. The model simulations were performed under conditions of realistic liver micro-architectures obtained from 3D reconstructions of confocal laser scanning micrographs. Interestingly, the established model predicted that initiated hepatocytes at first arrange in elongated patterns. Only when the tumor progresses to cell numbers of approximately 4000, does it adopt spherical structures. This prediction may have relevant consequences, since elongated tumors may reach critical structures faster, such as larger vessels, compared to a spherical tumor of similar cell number. Interestingly, this model prediction was confirmed by analysis of the spatial organization of initiated hepatocytes in a rat liver tumor initiation study using single doses of 250 mg/kg of the genotoxic carcinogen N-nitrosomorpholine (NNM). Indeed, small clusters of GST-P positive cells induced by NNM were elongated, almost columnar, while larger GDT-P positive foci of approximately the size of liver lobuli adopted spherical shapes. From simulations testing numerous possible mechanisms, only HSA could explain the experimentally observed initial deviation from spherical shape. The present study demonstrates that the architecture of small cell clusters of hepatocytes early after initiation is still controlled by physiological mechanisms. However, this coordinating influence is lost when the tumor grows to approximately 4000 cells, leading to further growth in spherical shape. Our findings stress the potential importance of organ micro-architecture in understanding tumor phenotypes.

[1]  S. Giordano,et al.  MicroRNA/gene profiling unveils early molecular changes and nuclear factor erythroid related factor 2 (NRF2) activation in a rat model recapitulating human hepatocellular carcinoma (HCC) , 2014, Hepatology.

[2]  A. Zell,et al.  Tumor promotion and inhibition by phenobarbital in livers of conditional Apc-deficient mice , 2016, Archives of Toxicology.

[3]  Curtis D. Klaassen,et al.  Casarett and Doull's Toxicology. The Basic Science of Poisons , 1981 .

[4]  D. Drasdo,et al.  Buckling instabilities of one-layered growing tissues. , 2000, Physical review letters.

[5]  Jean Paul Thiery,et al.  Johnson-Kendall-Roberts theory applied to living cells. , 2005, Physical review letters.

[6]  Nick Jagiella,et al.  Simulating tissue mechanics with agent-based models: concepts, perspectives and some novel results , 2015 .

[7]  S. Moolgavkar,et al.  Quantitative analysis of tumor initiation in rat liver: role of cell replication and cell death (apoptosis). , 2000, Carcinogenesis.

[8]  Alissa M. Weaver,et al.  Tumor Morphology and Phenotypic Evolution Driven by Selective Pressure from the Microenvironment , 2006, Cell.

[9]  Xiaoye S. Li,et al.  An overview of SuperLU: Algorithms, implementation, and user interface , 2003, TOMS.

[10]  Jan Böttger,et al.  Protocols for staining of bile canalicular and sinusoidal networks of human, mouse and pig livers, three-dimensional reconstruction and quantification of tissue microarchitecture by image processing and analysis , 2014, Archives of Toxicology.

[11]  M. Loeffler,et al.  Individual-based models to growth and folding in one-layered tissues: intestinal crypts and early development , 2001 .

[12]  A. Zell,et al.  Proinflammatory mesenchymal effects of the non-genotoxic hepatocarcinogen phenobarbital: a novel mechanism of antiapoptosis and tumor promotion. , 2015, Carcinogenesis.

[13]  J. Wagg,et al.  Vascular phenotype identification and anti-angiogenic treatment recommendation: A pseudo-multiscale mathematical model of angiogenesis. , 2016, Journal of theoretical biology.

[14]  Gruler,et al.  Random walk and directed movement: Comparison between inert particles and self-organized molecular machines. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[15]  Jan G. Hengstler,et al.  TiQuant: Software for tissue analysis, quantification and surface reconstruction , 2014, Bioinform..

[16]  Robert J Gillies,et al.  Impact of metabolic heterogeneity on tumor growth, invasion, and treatment outcomes. , 2015, Cancer research.

[17]  C. Zhu,et al.  Determining force dependence of two-dimensional receptor-ligand binding affinity by centrifugation. , 1998, Biophysical journal.

[18]  Kurt Binder,et al.  Monte Carlo Simulations in Statistical Physics , 2009, Encyclopedia of Complexity and Systems Science.

[19]  Easl EASL-EORTC clinical practice guidelines: management of hepatocellular carcinoma , 2012 .

[20]  C. Zhu,et al.  Measuring two-dimensional receptor-ligand binding kinetics by micropipette. , 1998, Biophysical journal.

[21]  U. Hofmann,et al.  Model-guided identification of a therapeutic strategy to reduce hyperammonemia in liver diseases. , 2016, Journal of hepatology.

[22]  Z. Stachura,et al.  Elasticity of normal and cancerous human bladder cells studied by scanning force microscopy , 1999, European Biophysics Journal.

[23]  S. Hoehme,et al.  On the Role of Physics in the Growth and Pattern Formation of Multi-Cellular Systems: What can we Learn from Individual-Cell Based Models? , 2007 .

[24]  R. Guthke,et al.  Integrated metabolic spatial‐temporal model for the prediction of ammonia detoxification during liver damage and regeneration , 2014, Hepatology.

[25]  D. Beysens,et al.  Cell sorting is analogous to phase ordering in fluids. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[26]  O. Dirsch,et al.  Cholestasis‐induced adaptive remodeling of interlobular bile ducts , 2016, Hepatology.

[27]  D. Drasdo,et al.  A single-cell-based model of tumor growth in vitro: monolayers and spheroids , 2005, Physical biology.

[28]  Alexander R. A. Anderson,et al.  Single-Cell-Based Models in Biology and Medicine , 2007 .

[29]  Stefan Hoehme,et al.  A cell-based simulation software for multi-cellular systems , 2010, Bioinform..

[30]  Vittorio Cristini,et al.  Patient-calibrated agent-based modelling of ductal carcinoma in situ (DCIS): from microscopic measurements to macroscopic predictions of clinical progression. , 2012, Journal of theoretical biology.

[31]  Nick Jagiella,et al.  Inferring Growth Control Mechanisms in Growing Multi-cellular Spheroids of NSCLC Cells from Spatial-Temporal Image Data , 2016, PLoS Comput. Biol..

[32]  M. Hayakari,et al.  Nonclonal growth of preneoplastic cells positive for glutathione S‐transferase P‐form in the rat liver , 2012, Cancer science.

[33]  Ursula Klingmüller,et al.  The virtual liver: state of the art and future perspectives , 2014, Archives of Toxicology.

[34]  R. A. Malmgren Observations on a liver mitotic stimulant present in tumor tissue. , 1956, Cancer research.

[35]  Ignacio Ramis-Conde,et al.  From genotypes to phenotypes: classification of the tumour profiles for different variants of the cadherin adhesion pathway , 2012, Physical biology.

[36]  Ignacio Ramis-Conde,et al.  Modeling the influence of the E-cadherin-beta-catenin pathway in cancer cell invasion: a multiscale approach. , 2008, Biophysical journal.

[37]  Stefan Hoehme,et al.  How predictive quantitative modelling of tissue organisation can inform liver disease pathogenesis. , 2014, Journal of hepatology.

[38]  H. Holzhütter,et al.  Modeling function–perfusion behavior in liver lobules including tissue, blood, glucose, lactate and glycogen by use of a coupled two-scale PDE–ODE approach , 2014, Biomechanics and Modeling in Mechanobiology.

[39]  F. MacKintosh,et al.  Scanning probe-based frequency-dependent microrheology of polymer gels and biological cells. , 2000, Physical review letters.

[40]  Matthias Hermes,et al.  Prediction and validation of cell alignment along microvessels as order principle to restore tissue architecture in liver regeneration , 2010, Proceedings of the National Academy of Sciences.

[41]  Michael Thomas,et al.  The novel histologic International Association for the Study of Lung Cancer/American Thoracic Society/European Respiratory Society classification system of lung adenocarcinoma is a stage-independent predictor of survival. , 2012, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[42]  S. V. Sotirchos,et al.  Glucose diffusivity in multicellular tumor spheroids. , 1988, Cancer research.

[43]  A. Henney,et al.  The virtual liver: a multidisciplinary, multilevel challenge for systems biology , 2012, Wiley interdisciplinary reviews. Systems biology and medicine.

[44]  M. Roncalli,et al.  YAP activation is an early event and a potential therapeutic target in liver cancer development. , 2014, Journal of hepatology.

[45]  Dirk Drasdo,et al.  Coarse Graining in Simulated Cell Populations , 2005, Adv. Complex Syst..

[46]  A. Henney,et al.  Unraveling liver complexity from molecular to organ level: challenges and perspectives. , 2015, Progress in biophysics and molecular biology.

[47]  R. Schulte‐Hermann,et al.  Antagonistic effects of selenium and lipid peroxides on growth control in early hepatocellular carcinoma , 2012, Hepatology.

[48]  S. Døskeland,et al.  Cell cycle parameters of adult rat hepatocytes in a defined medium. A note on the timing of nucleolar DNA replication , 1987, Journal of cellular physiology.

[49]  O. Dirsch,et al.  A biphasic model for sinusoidal liver perfusion remodeling after outflow obstruction , 2010, Biomechanics and modeling in mechanobiology.

[50]  S. Giordano,et al.  MicroRNA/gene profiling unveils early molecular changes and NRF2 activation in a rat model recapitulating human HCC. , 2014 .

[51]  G. Dai,et al.  Four Waves of Hepatocyte Proliferation Linked with Three Waves of Hepatic Fat Accumulation during Partial Hepatectomy-Induced Liver Regeneration , 2012, PloS one.

[52]  G. Oster,et al.  How do sea urchins invaginate? Using biomechanics to distinguish between mechanisms of primary invagination. , 1995, Development.

[53]  Ben Fabry,et al.  Microrheology of human lung epithelial cells measured by atomic force microscopy. , 2003, Biophysical journal.

[54]  Adam S. Hayward,et al.  Recent advances in 2D and 3D in vitro systems using primary hepatocytes, alternative hepatocyte sources and non-parenchymal liver cells and their use in investigating mechanisms of hepatotoxicity, cell signaling and ADME , 2013, Archives of Toxicology.

[55]  M. Loeffler,et al.  Modeling the effect of deregulated proliferation and apoptosis on the growth dynamics of epithelial cell populations in vitro. , 2005, Biophysical journal.

[56]  Sylvain V Costes,et al.  Phenotypic transition maps of 3D breast acini obtained by imaging-guided agent-based modeling. , 2011, Integrative biology : quantitative biosciences from nano to macro.

[57]  G. Ippolito,et al.  Induction of autophagy promotes the growth of early preneoplastic rat liver nodules , 2015, Oncotarget.