Mathematical modelling of liver regeneration after intoxication with CCl(4)

Liver regeneration is a complex process, having evolved to protect animals from the consequences of liver loss caused by food toxins. In this study, we established a mathematical spatial-temporal model of the liver lobule regenerating after CCl(4) intoxication. The aim of modelling the regeneration process by matching experimental observations with those from a mathematical model is to gain a better understanding of the process and to recognize which parameters are relevant for specific phenomena. In order to set up a realistic minimal model, we first reconstructed a schematised liver lobule after determination of: (i) the mean number of hepatocytes between the central vein and the periphery of the lobule, (ii) the mean size of the hepatocytes and (iii) the mean number of hepatocyte columns in the inner, midzonal and peripheral ring of the lobule. In a next step, we determined the time course of cell death and BrdU incorporation after intoxication of male Sprague Dawley rats with CCl(4), thereby differentiating between inner, midzonal and peripheral hepatocytes. These parameters were used to construct a model. The basic unit of this model is the individual cell. The detailed behaviour of the cells is studied, controlled by the model parameters: (1) probability of cell division at defined positions of the lobule at a given time, (2) "coordinated cell orientation", i.e., the ability of the cells to align during the regeneration process into columns towards the central vein of a liver lobule, (3) cell cycle duration, (4) the migration activity and (5) the polarity of the hepatocytes resulting in polar cell-cell adhesion between them. In a schematised lobule, the model shows that CCl(4) initially induced cell death of a pericentral ring of hepatocytes, followed by a wave of proliferation that starts in the surviving hepatocytes next to the inner ring of dead cells and continues to the peripheral hepatocytes, finally restoring the characteristic micro-architecture of the lobule in a 7-day process. This model was used to systematically analyze the influence of parameters 1-5. Interestingly, coordinated cell orientation and cell polarity were identified to be the most critical parameters. Elimination led to destruction of the characteristic micro-architecture of the lobule and to a high degree of disorder characterized by hexagonal cell structures. Our model suggests that the ability of hepatocytes to realign after cell division by a process of coordinated cell orientation (model parameter 2) in combination with cell polarity (model parameter 5) may be at least as critical as hepatocyte proliferation (model parameter 1) itself.