(533g) Modeling the Spatio-Temporal Dynamics of Calcium Signal Propagation in Liver Lobules
Although the behavior of Ca2+ signal propagation at the lobular scale is understood fairly well, the contribution of hepatocyte-specific Ca2+ signaling components to the lobular scale Ca2+ wave propagation remains largely unknown. In the present study, we performed a causality analysis of time series data on Ca2+ signals in single hepatocytes to identify a cell-cell influence network underlying the Ca2+ wave propagation in liver lobules. We also developed a computational model of spatio-temporal Ca2+signaling and cell-cell interactions to predict the coordinated changes in the intracellular Ca2+ signaling parameters sufficient to yield the lobular scale Ca2+ propagation dynamics identified in the time series data analysis.
We acquired time-resolved images of Ca2+ spikes caused by different vasopressindoses in the isolated perfusedrat liver at the lobular level. We segmented the acquired images to obtain hepatocyte-specific cytosolic Ca2+ profiles. For a given time lag, we estimated pairwise transfer entropies (TEs) for Ca2+ signal information received by a hepatocyte from all the other hepatocytes within the two-dimensional lobular visual field to obtain hepatocyte- and lag-specific empirical TE distributions. For each hepatocyte, we identified the set of statistically significant causal influences by comparing the TE values from its immediate neighbors with the hepatocyte-specific empirical TE distribution. Our analysis revealed that at the lobular scale, for the vasopressin levels employed in the experiment, Ca2+ waves propagate from the pericentral to the periportal region along the hepatocyte cords. However, the transfer entropy analysis revealed that the direction of causal influence between successive neighboring hepatocytes may not be consistently in the pericentral-to-periportal direction, i.e., in some instances, Ca2+ signal information may â??flowâ? in the direction opposite to lobular scale Ca2+ wave propagation between neighboring hepatocytes lying within a hepatocyte cord.
We utilized a computational modeling approach to explore how such local switches in the direction of Ca2+ signal influence could arise, and evaluated their influence on the overall lobular scale wave propagation. We developed an ODE based computational model of cytosolic Ca2+ spiking in hepatocytes to identify the spatial patterns of cell-intrinsic Ca2+ signaling parameters that can reproduce the experimentally observed lobular Ca2+ wave propagation. Our model consists of individual hepatocytes in which the intracellular Ca2+ signals are elicited via a G-protein coupled receptor mediated pathway. The hepatocytes are oriented along a cord and exchange IP3 with the adjacent hepatocytes via gap junctions. Model simulations revealed that 1) gap junction coupling results in coordination of Ca2+ spikes in hepatocytes along a cord, despite variability in intracellular Ca2+ signaling parameters across cells, and 2) correlated spatial gradients of intracellular Ca2+ signaling parameters are required for yielding unidirectional flow of Ca2+ waves at the lobular scale. In the absence of gap junction coupling, Ca2+ spikes did not occur sequentially along a hepatocyte cord in response to stimuli that lasted over extended periods of time. Disruption of spatial gradients of intracellular signaling parameters resulted in local deviations from the overall direction of Ca2+ wave propagation along hepatocyte cords.
Next, we tested the effect of local disruptions in cell-specific Ca2+ signaling ability of hepatocytes on the lobular scale Ca2+ wave propagation. Inability of hepatocytes to respond to Ca2+ agonists may contribute to aberrant bile trafficking and promote cholestasis. We initialized our model of a hepatocyte cord in which a randomly distributed subset of hepatocytes was insensitive to Ca2+ agonists. Simulations revealed that Ca2+ waves do not propagate throughout the lobule when Ca2+ signals are not elicited in a subset of hepatocytes. Furthermore, increase in the strength of gap junction coupling was sufficient to restore the lobular scale Ca2+ signaling pattern even in the presence of local disruptions in Ca2+ response, suggesting that efficient intercellular communication can yield robust lobular scale Ca2+ response despite local insufficiencies in Ca2+ signaling.
In summary, we have developed and illustrated a computational model based approach to analyze and interpret time series data on Ca2+ signals in single cells in intact liver tissue. Our results led to new predictions on the role of gap junctions in overriding the cell-intrinsic variability towards a robust and spatially-organized lobular scale Ca2+ response to external stimuli.