472778 Modeling the Spatio-Temporal Dynamics of Calcium Signal Propagation in Liver Lobules

Wednesday, November 16, 2016: 2:18 PM
Continental 8 (Hilton San Francisco Union Square)
Aalap Verma1,2, Anil Antony3, Hirenkumar Makadia4, Jan Hoek3, Babatunde A. Ogunnaike5 and Rajanikanth Vadigepalli6, (1)Department of Pathology, Anatomy and Cell Biology, Thomas Jefferson University, Daniel Baugh Institute for Functional Genomics and Computational Biology, Philadelphia, PA, (2)Biomedical Engineering, University of Delaware, Newark, DE, (3)Department of Pathology, Anatomy and Cell Biology, Thomas Jefferson University, Philadelphia, PA, (4)Daniel Baugh Institute for Functional Genomics and Computational Biology, Department of Pathology Anatomy and Cell Biology, Thomas Jefferson University, Philadelphia, PA, (5)Department of Chemical Engineering, University of Delaware, Newark, DE, (6)Daniel Baugh Institute for Functional Genomics and Computational Biology, Department of Pathology, Anatomy and Cell Biology, Thomas Jefferson University, Philadelphia, PA

The liver performs a wide variety of physiological functions such as regulation of intermediary metabolism, lipid synthesis, bile production and xenobiotic detoxification. Normal liver function requires tight regulation of intracellular processes as well as intercellular coordination. Free Ca2+ in the intracellular domain participates in the regulation of hepatic functions such as energy metabolism, bile secretion, proliferation and apoptosis. Regulation of cytosolic Ca2+ is particularly important in hepatocytes, the cells responsible for the bulk of metabolic and detoxification activities in the liver. Consequently, disruption of Ca2+ dynamics can potentially lead to pathological conditions. Extracellular stimulants such as hormones elicit a signaling cascade involving phospholipase C (PLC) cleavage of PI-4,5-P2, IP3 generation, and subsequent Ca2+ release from intracellular stores. Gap-junction mediated transfer of IP3 molecules between adjacent hepatocytes induces, in neighboring hepatocytes, Ca2+ oscillations that spread through the tissue. This cell-cell interaction leads to a coordinated response at the tissue level. Measuring lobular Ca2+ responses upon treatment of perfused rodent livers with PLC-activating agonists has revealed a remarkable spatial organization. For low stimulus levels, Ca2+ signals in liver lobules start near a specific lobular locus, in the pericentral region, and propagate towards the periportal region in a wave-like fashion.

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.

Extended Abstract: File Not Uploaded