The liver is the primary metabolic processor of nutrients, and hence, it is deeply involved in digestive physiology. Feeding behavior is a complex set of events controlled by neural and hormonal communication between structures of the central nervous system (mainly hypothalamic structures such as the arcuate, lateral, ventromedial and paraventricular nuclei) and a variety of organs related to nutrient handling in the digestive tract . Given that in nature the food sometimes is scarce, organisms have developed strategies to optimize the finding of food, the processing of nutrients, and the assimilation of biomolecules precisely when mealtime result a predictable event. Part of these adaptations is a timing system that underlies the events involved in the circadian rhythmicity. It is in this context that, in conditions of restricted access to food, the expression of an oscillator synchronized by food (FEO) is observed, which is an circadian alternative to the role played by the suprachiasmatic nuclei [2,3]. When animals are under a restricted feeding schedule over a period of several weeks, they display a behavior known as food anticipatory activity (FAA), an arousal behavior that precedes the availability of food. This FAA behavior is associated with the FEO .
The anatomical location of the FEO is still undetermined. However, the liver is likely playing a role in FEO physiology since it acts as a time-driven metabolic integrator of nutrients and has been associated with the control of the hunger-satiety cycle [5 - 7]. For example, a decrease in liver ATP triggers feeding behavior as well as an increase in the hepatic cytosolic Ca2+ concentration after administration of the metabolic inhibitor fructose analogue 2,5-anhydro-D-mannitol (2,5-AM) . Upon re-feeding, rats display compensatory hyperphagia with a concomitant increase in ATP synthesis . Many parameters related to hepatic physiology display circadian rhythmicity driven by the so-called clock genes, but they are also entrained by food availability . In this context, the liver acts as a peripheral oscillator with the capacity to exhibit circadian rhythms in metabolic and physiological activities, eliciting an extremely rapid response when a restricted feeding (RF) protocol is applied .
During FAA, changes occur in: 1) the expression of ~80 hepatic genes related to biochemical reactions , 2) the cytoplasmic and mitochondrial redox state, which becomes oxidized , 3) ATP levels and mitochondrial respiration, which show an elevation , and 5) glycogenolytic activity, which is reduced in comparison with 24-h-fasted rats . In addition, RF promotes lipid mobilization from adipose tissue followed by an increased catabolism of fatty acids within the liver [16 - 18]. The metabolic pattern found in rats entrained by daily restricted-food access/expressing the FEO is different from the one displayed by rats fed ad-libitum and under a 24 h of fasting, suggesting that during the RF protocol, hepatic physiology adopts a novel regulatory condition known as rheostasis or allostasis, a term meaning stability through regulated changes .
In the liver, intracellular calcium modulates glucose metabolism (glycogenolytic and glyconeogenic activities), protein folding, mitochondrial function, gene transcription, apoptosis, cell proliferation, and bile secretion . Calcium modulates all these processes based on the temporal and spatial transients that function as a metabolic and transcriptional control code . Intracellular calcium dynamics involves the coordinated action of a variety of proteins responsible for calcium mobilization outside of and within the cytosolic space. The liver expresses the 2 principal intracellular, calcium-release channels: the inositol 1,4,5-trisphosphate receptor (IP3R) (types 1 and 2) [22,23] and the ryanodine receptor type 1 (RyR), detected as a truncated but functional channel-protein . The hepatic metabolic pumps that extrude cytosolic calcium are the sarco/endoplasmatic reticulum calcium ATPase (SERCA) splicing isoform 2b and the plasmatic membrane calcium ATPase (PMCA) type 1 and 4, but isoform 2 is also expressed at a low level .
Liver acinus shows two distinctive zones in which hepatocytes exhibit biochemical heterogeneity depending on the arrival of oxygen and nutrients. Hepatocytes next to portal vein are named periportal (PP) and show higher rates of glyconeogenesis, urea synthesis, bile formation, and lipid catabolic activity. Hepatocytes near the central vein are called pericentral (PC), and show mainly glycolysis, glycogenolysis, a high content of cytochromes P-450, and detoxification activities . It was reported that IP3R is heterogeneously distributed along the hepatic acinus, being more abundant in the PP than in the PC zone . No reports exist regarding liver zonational distribution for other calcium-handling proteins.
Liver calcium signaling responds to the energy status and hence to the feeding condition. Hepatocytes from fasting rats show significantly higher cytosolic calcium levels . In addition, as mentioned previously , rats treated with 2,5-AM showed a significant reduction of hepatic ATP levels and a concomitant elevated intracellular calcium response. Our hypothesis is that the biochemical properties and zonal location of receptor-channels and ATPases handling intracellular calcium in the liver will be regulated by the timing system and will be responsive to a feeding protocol of food access restricted to 2 h during daytime. Hence, the aim of this report was to explore the influence of restricted feeding and the associated FEO expression on the daily variations of the main hepatic intracellular calcium-handling proteins: IP3Rs, RyR, SERCA, and PMCAs.
Antibodies against PER1 (sc-7724), IP3R1(sc-6093) and IP3R2 (sc-7278), SERCA2 (sc-8094), PMCA1 (sc-16488) and PMCA4 (sc-22080), and Actin as well as alkaline phosphatase (AP)-conjugated rabbit anti-goat and goat anti-mouse secondary antibodies were obtained from Santa Cruz Biotechnology (Santa Cruz, CA, USA), and the Ryanodine Receptor (ab9078) was from Millipore (MA, USA). [3H]-IP3 and [3H]-ryanodine were purchased from New England Nuclear (NEN, MA, USA). (1,4,5)-Inositol trisphosphate and ryanodine were from Calbiochem (CA, USA). Protease inhibitors, Stains all reactive and all other chemicals were obtained from Sigma (MO, USA), and Western blot equipment and reagents were from Bio-Rad (CA, USA).
Animals and housing
Adult male Wistar rats weighing 200 ± 20 g (11–12 weeks old) at the beginning of the experiment were maintained under a 12 h:12 h light–dark cycle (light on at 08:00 h) at constant temperature (22 ± 1°C). Rats were kept in separate groups of 4 in transparent acrylic cages (40 × 50 × 20 cm), with free access to water and balanced Purina Chow meal except during food restriction, fasting, or re-feeding conditions. Illumination during the light period was obtained from 40 W fluorescent bulbs that generated 120 lux at the cage lid. All experimental procedures were conducted in accordance with our Institutional Guide for Care and Use of Animal Experimentation (Universidad Nacional Autónoma de México) and in agreement to international ethical standards .
Rats were randomly assigned to 4 groups: 1) rats fed ad libitum (AL) for 3 weeks; 2) rats exposed to a daily restricted feeding schedule (RF) with access to food only between ZT4 to ZT6 (ZT0 is the time of lights on) for 3 consecutive weeks; 3) rats under 21 h of food deprivation (Fasted) starting at ZT6, and 4) rats that were fasted for 22 h (starting at ZT6 on the first day) and re-fed (Refed) for 2 h (from ZT4 to ZT 6 on the second day). At the end of the third week, rats from the AL and RF groups were sacrificed at 3-h intervals to complete a 24-h, day-night cycle (from ZT0 to ZT21). The Fasted and Refed group were sacrificed at ZT3 and ZT6, respectively. RF group was used to characterize the effect of temporal restricted feeding on calcium-handling proteins 24 h profiles, meanwhile Fasted and Refed groups were used as acute feeding condition controls for the RF groups at ZT3 (before food access and during FAA) and ZT6 (after feeding) respectively .
Samples were collected from trunk blood and centrifuged at 4,000 g for 15 min at 4°C. Plasma was then stored at −80°C for subsequent measurements of corticosterone. Corticosterone concentrations were measured in duplicate using a commercial ELISA kit (Assay Designs, MI, USA).
Hepatic tissue was fractionated as reported by . The liver was removed (≈5 g), immediately placed in ice-cold homogenization buffer (HB, 1:10 w/v), and disrupted with a Potter-Elvehjam teflon-glass homogenizer (40 rpm for 15–20 s). The HB contained: 225 mM sucrose, 0.3 mM EGTA, 10 mM Tris/HCl (pH 7.4), and 1 mM DTT, supplemented with a mixture of protease inhibitors (0.1 mM PMSF, 0.1 mM benzamidine, 10 μM pepstatin A, 1 μg/ml aprotinin, 1 μg/ml o-phenanthroline, and 10 μg/ml leupeptin). The liver homogenate was centrifuged at 1,000 g for 15 min (in a Sorvall SS34 centrifuge), and the resulting supernatant was decanted. The supernatant was centrifuged 2 times at 7,700 g for 15 min to precipitate the mitochondrial fraction. The resultant supernatant was ultracentrifuged (Beckman 70Ti rotor) at 100,000 g for 60 min. The pellet, corresponding to the endoplasmic reticulum fraction (ER), and the supernatant, corresponding to the cytosol fraction, were collected, aliquoted, and kept at −70°C. The plasma membrane fraction (PM) was obtained from the first 7700 × g pellet, as described by . The pellet was resuspended in 20 ml of 250 mM sucrose - 10 mM Tris/HCl (pH 7.5) solution, then mixed with 2.6 ml of Percoll (1.13 g/ml density) and 0.4 ml of 2 M sucrose, placed at the top of a Percoll gradient, and centrifuged at 35,000 g for 20 min to yield the crude PM fraction in the supernatant. This fraction was carefully recovered (~15 ml), layered on 5 ml of the same buffer previously described, with the addition of sucrose/Tris buffer containing 1.3 M CaCl2, and centrifuged at 37,000 g for 25 min; the PM fraction was obtained from the middle of the gradient, separated from other cellular components. The total membrane fraction from skeletal muscle was processed as described by  and was used as a tissue control as refereed forward. Subcellular fractions were resuspended in HB and stored in aliquots at −70°C until further use. Protein concentration was determined by the Lowry method  using bovine serum albumin (BSA) as standard.
Glucose-6-phosphatase activity (EC 126.96.36.199) was assayed as described by  and used as an endoplasmic reticulum (ER) marker; 5´nucleotidase activity (EC 188.8.131.52) was determined following the method of , and was used to evaluate the purity of the plasma membrane fraction (PM). Both enzymes were determined first in the total homogenate and later in the isolated fractions.
The [3H]-IP3 binding assay was performed as reported by Furiuchi et al.  with some modifications: 100–300 μg of microsomal and PM fractions were incubated in triplicate for 30 min in 120 μl of a solution containing 25 mM Tris/HCl (pH 8.0), 5 mM NaHCO3, 1 mM EDTA, 0.25 mM DTT, and 0.5 to 100 nM of [3H]-IP3. Non-specific binding was measured as the radioactivity not displaced by non-radioactive 10 μM IP3. Each sample was washed 5 times with 5 ml of a cold buffer containing 25 mM Tris/HCl (pH 8.0), 5 mM NaHCO3, and 1 mM EDTA. Filters were counted in 10 ml tritosol  in a LS6500 Beckman multi-purpose scintillation counting system.
[3H]-Ryanodine binding was evaluated as described by Hamilton et al. . Briefly, 400 μg of microsomal protein was incubated in triplicate for 14–16 h with [3H]-ryanodine at concentrations ranging from 0.5 to 500 nM, in a final volume of 250 μl of incubating buffer: 300 mM KCl, 100 μM CaCl2, 100 μg/ml BSA, and 20 mM MOPS (pH 7.4) at room temperature. The assay was terminated by addition of 5 ml ice-cold 0.3 mM KCl and filtration on Whatman GF/F filters, followed by another 4 washes. Non-specific binding was defined as that not displaced by the addition of 10 μM of non-radioactive ryanodine. Radioactivity bound and binding parameters (Bmax and Kd) were determined as described previously for IP3R binding assay.
Ca2+ATPase activity assays
ATPase activities were measured by a standard coupled enzymatic assay in which the rate of ATP hydrolysis was linked to NADH oxidation, and the optical density was recorded at 340 nm (ϵ = 6.22 mM-1.cm-1). ER and PM samples were used for SERCA and PMCA assays as described by Saborido et al. . Aliquots (50 to 100 μg) of each fraction were incubated in buffer containing 25 mM MOPS (pH 7.4), 0.2 mM EGTA, 5 mM MgCl2, 100 mM KCl, 0.6 mM phosphoenolpyruvate, 2.4 unit/ml pyruvate kinase, 10 unit/ml lactate dehydrogenase, 4 μM ionophore A23187, 0.27 mM NADH, and 1 mM CaCl2 for total activity (21 mM CaCl2 for basal activity) in a final volume of 1 ml. After pre-incubation of the assay mixture for 5 min at 37°C, the reaction was started by adding 1 mM ATP (final concentration). SERCA and PUMP Ca2+-ATPase activities were calculated as total activity minus the basal activity. Thapsigargin (1 μM) and Eosin (2 μM) were used as SERCA and PMCA inhibitors, respectively, in order to test the specificity of the measured activity of the two pumps .
Western blot analysis
Cellular fractions were resuspended in SDS sample buffer in reducing conditions, as described by Laemli . Proteins were separated by SDS-PAGE using 5% acrylamide for IP3R (type 1 and 2) and RyR, and 7.5% acrylamide for PER1, SERCA2, and PMCA (types 1 and 4). After electrophoresis, proteins were transferred to Protean nitro-cellulose membranes and blocked with PBS containing 0.1% Tween and 5% defatted milk for 1 h at room temperature. Membranes were incubated with the primary antibody overnight. All the antibodies were used at 1:500 dilutions, except that the antibody used to detect the skeletal muscle RyR protein as positive control not detected the hepatic RyR protein (even using other RyR antibodies from the companies Santa Cruz and Abcam). Then stains-all protocol was applied in order to detect the hepatic RyR protein using its characteristic very high molecular weight (more than 500 kDa) as a defining parameter . After repeated washes with PBS-Tween buffer, membranes were incubated for 2 h with the appropriate alkaline phosphatase (AP)-conjugated secondary antibody at 1:5000 dilution, and the bands were visualized using the AP conjugate substrate kit (Bio-Rad, CA, USA) according to the manufacturer’s instructions. Membrane samples from cerebellum (IP3R1), kidney (IP3R2), skeletal muscle (RyR) and brain (SERCA2b and PMCA1-4 proteins) were used as positive controls. Blots were digitalized and analyzed with the Image J® software (version 1.38, USA).
To determine the acinus (PP and PC) distribution of RyR, IP3R1 and 2, SERCA2, and PMCA1 and 4, slices of rat liver (10 μm thick) were incubated with specific antibodies. Immunohistochemistry was performed as reported by Clair et al.  using freshly isolated liver slices fixed in 4% (w/v) paraformaldehyde-PBS, pH 7.5 overnight, and cryo-preserved in 30% (w/v) sucrose for another overnight cycle. After permeabilization with 0.05% Triton-PBS for 10 min and blocking with PBS containing 5% de-fatted milk plus 0.1% BSA, slices were incubated overnight at 4°C with primary antibody at concentration suggested by the manufacturer. After rinsing in PBS 3 times, the liver tissue slices were incubated at room temperature for 1 h with a secondary antibody conjugated to FITC (fluorescein 5-isothiocyanate) supplied by Sigma (MO, USA). Actin related to apical membrane was stained with Rhodamine-conjugated Phalloidine in order to outline hepatic cells .
Visualization and quantification of fluorescence
Liver slices were visualized using a CX31 Olympus microscope. Images were collected with a DP71 Olympus camera and visualized using Image-Pro Plus software (version 6.0, MD, USA). Throughout the study, standardized fluorescence was set in all proteins studied against their respective negative control (without specific primary antibody). Fluorescent intensity in PP and PC hepatocytes was quantified as described by Lahm et al. .
Calculations and statistics
The results are expressed as mean ± SEM of at least 4 individual experimental observations. Statistical analysis was done using the Prism version 5.0 program (GraphPad software, USA). To detect possible time-condition differences, RF and Al groups were compared with a two-way ANOVA. Significance was estimated by the Tukey test with an α level set at 0.05. For the feeding condition controls, a Student´s t-test was used to detect significant differences between fasted vs re-fed, fasted vs RF ZT3, re-fed vs RF ZT6, and RF ZT3 vs RF ZT6. For chronobiological analysis, first a one-way ANOVA was performed in each group and then a 24-h period single-cosinor method was used as previously described . For rhythmic interpretation of the results, the following parameters were considered: acrophase (time of peak value), MESOR (midline estimating statistic of rhythm), amplitude (half of the total variation of the rhythm), and rhythmicity, which corresponds to a p value (< 0.01) of an F test of fitting the original results to an expected sinusoidal curve with a 24-h period.
Establishment of the FEO expression / restricted feeding protocol
To confirm the metabolic and physiological adaptations associated with the protocol for RF (food access only during 2 daytime hours) and the FEO expression, serum corticosterone levels and the presence of liver PER1 protein were determined over a 24-h cycle (Figure 1). It is well accepted that RF promotes a significant shift in the peaks of these 2 parameters during the time of the FAA, prior to food access [45,46]. Indeed, it is observed in Figure 1A that AL and RF rats showed an elevation of circulating corticosterone previous to the end of the light period when the animals are in the transition between the sleep and the awake periods (around 50% from the trough; AL group, p < 0.009, one-way ANOVA). But, in addition, the RF group showed another, even larger peak at ZT3 before feeding (100% increase with respect to the trough; RF group, p < 0.004, one-way ANOVA). A similar pattern was observed in the 24-h rhythm of PER1 (Figure 1B): the AL group showed a clear PER1 peak during the dark period (ZT15–Z21; p < 0.009, one-way ANOVA) when the animals are active, confirming previous reports . In contrast, RF rats showed a shift of the PER1 peak towards the light period (ZT3-ZT6), which corresponds to the period immediately before and after feeding. According to the two-way ANOVA test, levels of corticosterone and PER1 protein from AL and RF groups (Figure 1) differed significantly, indicating a rhythmicity that was characterized by a Cosinor test. Corticosterone and PER1 daily rhythms from the AL group have significant 24-h rhythmicity (Table 1), but in the RF group these parameters exhibited a phase advance (of ~7 h and ~12 h, respectively). PER1, but not corticosterone, showed a significant 24-h rhythm (Table 1) in the RF group. This result is due to the biphasic pattern of daily corticosterone levels (one corresponding to light entrained before dark period and the other related to meal cue before food access time). Both results demonstrated that our protocol successfully induces hepatic entrainment due to daytime restricted feeding.
Acrophase ZT (h)
Rhythmicity ( P < 0.05 )
||09:30 ± 00:26
||02:15 ± 00:35*
||3.13 ± 0.23
||3.39 ± 0.30
||17.6 ± 3.9
||20.7 ± 4.3
|PER1 (WB)||17:30 ± 00:19||05:06 ± 00:24*||0.27 ± 0.015||0.22 ± 0.011||0.55 ± 0.07||0.64 ± 0.06||0.002||0.04|
IP3R And RyR: ligand binding properties and protein expression during restricted feeding / FEO expression
Diurnal fluctuations of liver IP3R and RyR were determined in subcellular fractions. To characterize the hepatic endoplasmic reticulum (ER) and plasma membrane (PM) fractions, the specific activities of glucose-6-phosphatase (ER marker) and 5´nucleotidase (PM marker) were measured (data not shown). The recovery yield was similar in AL and RF groups at the different tested time (data not shown). The maximum contamination of PM in the ER fraction was 13%, whereas the ER contamination in the PM fraction was 18%. According to the one-way ANOVA test, the [3H]-IP3 binding in ER membranes of the AL group showed significant differences over the 24-h period (Figure 2A). Although the RF group showed significant differences at certain time points, there was a clear decrease in the amplitude of [3H]-IP3 binding rhythm. Both groups showed 24-h rhythmicity. In addition, the RF group showed an almost 50% reduction in amplitude and MESOR as well as a phase advance of 8 h with respect to the AL group (Table 2). Interestingly, [3H]-IP3 binding in the ER fraction of the Fasted group increase in comparison to the Refed group but also reach the highest value of all experimental groups. The same pattern occur in AL group since [3H]-IP3 binding decrease after rats start it activity period (Figure 2A). In contrast, the RF group did not show any difference between ZT3 (before feeding) and ZT6 (just after feeding), but both values were significantly lower than the corresponding control groups of feeding condition (Figure 2B). In the hepatic ER fraction, the IP3R type 1 is the main isoform, but type 2 is also present at lower levels . Hence, the [3H]-IP3 binding was also measured in the PM fraction which has been suggested as an IP3R type 2-enriched subcellular site ; the data indicated significant differences among time points for both groups (p < 0.001, one-way ANOVA, Figure 2C) as well as between the AL and RF groups (p < 0.01, one-way ANOVA). These results are consistent with the differences shown in Table 2 regarding 24-h rhythmicity: RF promoted a 5-h advance, a 50% reduction of amplitude (relative to the AL group) and a ~25% decrease of MESOR. Regarding feeding condition control groups, a decrease in [3H]-IP3 binding was observed in the PM fraction of the Refed group (Figure 2D). However, the RF group again displayed a dissimilar pattern: At ZT6, RF rats showed a significant PM [3H]-IP3 binding elevation in comparison to the Refed group (Figure 2D).
Acrophase ZT (h)
Rhythmicity ( P < 0.05 )
Ca2+ releasing proteins
IP3R (ER Binding assay)
||12:20 + 00:21
||04:02 + 00:20*
||1.3 + 0.22
||0.6 + 0.10*
||2.6 + 0.44
||1.5 + 0.21*
||12:05 + 00:35
||05:37 + 00:47*
||0.28 + 0.04
||0.16 + 0.03*
||0.56 + 0.08
||0.37 + 0.06*
IP3R (PM Binding assay)
||09:57 + 00:48
||04:22 + 00:28*
||42.8 + 10.3
||36.3 + 6.1
||82.9 + 21.7
||61.9 + 12.3*
||11:26 + 00:29
||04:03 + 00:31*
||0.29 + 0.034
||0.25 + 0.036
||0.54 + 0.07
||0.31 + 0.07*
RyR (Binding assay)
||17:43 + 00:41
||17:10 + 00:32
||3.2 + 0.46
||4.3 + 0.78
||9.7 + 0.9
||10.2 + 1.5
||17:50 + 00:50
||17:30 + 00:58
||0.29 + 0.033
||0.22 + 0.052
||0.59 + 0.06
||0.62 + 0.10
|Ca2+ extruding proteins||
||18:35 + 00:31
||03:51 + 00:39*
||5.36 + 0.9
||20.9 + 2.4*
||23.2 + 2.1
||47.4 + 5.8*
||16:48 + 00:29
||02:55 + 00:24*
||0.24 + 0.02
||0.34 + 0.03*
||0.49 + 0.04
||0.66 + 0.07*
||21:48 + 01:09
||23:44 + 01:07
||6.2 + 1.2.
||2.7 + 0.8*
||19.8 + 2.3
||21.6 + 1.4
||01:35 + 0:59
||09:15 + 00:45*
||0.1 + 0.02
||0.2 + 0.03
||0.8 + 0.04
||0.67 + 0.06
|PMCA4 (WB)||12:48 + 00:32||17:37 + 00:24*||0.24 + 0.04||0.3 + 0.05||0.66 + 0.08||0.53 + 1.0||0.06NS||0.04|
The RyR, the other important intracellular calcium-releasing channel in hepatic tissue, was also studied for ligand binding. According to a one-way ANOVA, [3H]-ryanodine binding showed differences among time points in both AL and RF groups (Figure 2E). Interestingly, there were no differences between these groups, as tested by a two-way ANOVA test, suggesting that the RF schedule did not affect the diurnal rhythmicity of the RyR. Clear differences in [3H]-ryanodine binding were observed in the RF group between ZT3 (before feeding) and ZT6 (after feeding) but not in the control groups (Fasted vs Refed) (Figure 2F). Hence, the RF group exhibited discrete modifications in the daily rhythm of [3H]-ryanodine binding (a phase advance of ~3 h, Table 2), but a significant effect at ZT3 in comparison to the Fasted group.
To investigate if these changes were associated with fluctuating levels of receptor expression, IP3R and RyR were analyzed by Western blot. IP3R1 (the main isoform in the hepatic ER fraction) showed a pattern similar to the [3H]-IP3 binding assay: AL and RF displayed 24-h rhythmicity (one-way ANOVA, Figure 3A), and significant differences were found due to time and feeding condition in both groups (tested by two-way ANOVA). However, very different characteristics were observed in their daily rhythm: the AL group showed a clear peak at the beginning of the dark period (ZT12), whereas in the RF group the peak shifted 8 h toward the middle part of the light period (before and after RF schedule) and had a 50% lower amplitude (Table 2). Acute fasting promoted a significant increase in the IP3R1 expression that was reverted upon re-feeding (Figure 3B). This effect was not observed in the rats under RF, since the level of IP3R1 was similar before and after food access (RF at ZT3 and ZT6, Figure 3B). The PM fraction was examined for the presence of hepatic IP3R type 2. The level of liver IP3R2 protein showed daily rhythmicity in the AL and RF groups, but the amplitude was ~50% lower in the latter (Figure 3C); in addition, RF promoted a shift in the IP3R2 peak, from the transition between light and dark periods (ZT12) observed in AL rats, to the time of daytime feeding (ZT6) (Table 2), a phase advance of 6 h. Fasted and Refed controls also showed an evident decrease in comparison to the RF group at ZT3 and ZT6, respectively (Figure 3D). These results strongly suggest that the diurnal rhythms of hepatic IP3R1 and 2 are profoundly modified by RF (resulted from 3 weeks) but in a distinct way than in the Fasted and Refed groups, which resulted from an acute condition (1 day).
The other calcium-release channel located in liver microsomal membranes, the RyR, also showed a clear diurnal rhythm in the AL group, with a significant elevation of the RyR protein (≈40%) during the dark period (Figure 3E). RF promoted 2 main changes in this rhythm: 1) after feeding (ZT6), an increase (≈80%) was observed; 2) at the beginning of the dark period (ZT15) a second, even larger, peak was detected. In addition, after food access (ZT6), the RF group showed an increase (≈30%) in comparison to the group before meal time (ZT3) that resulted significant among the Fasted versus Refed group (Figure 3F).
SERCA and PMCA: activity and protein expression during restricted feeding / FEO expression
SERCA activity in the AL group did not change over 24 h (Figure 4A). In contrast, RF rats showed a clear peak and an evident increase in the SERCA activity during the light phase (P < 0.001; one-way ANOVA, Figure 4A). Hence, significant rhythmicity was detected in the RF, but not in the AL group (Table 2). This clear difference associated with the restricted feeding protocol/FEO expression was even more remarkable since no changes were detected in feeding control rats (Fasted vs Refed). PMCA activity was reduced (≈50%) in the AL group at the end of the light period, whereas RF rats showed no daily pattern of PMCA activity (except for a reduction observed at ZT9-12) (Figure 4C); no daily rhythmicity was detected for these 2 groups. In addition, no change was observed in the feeding control groups (Figure 4D).
To complement the enzymatic measurements of SERCA and PMCA activity, the subcellular fractions were assayed for these proteins by Western blot. The results shown in Figure 5A confirmed that SERCA2 protein (the isoform expressed in liver) did not change significantly during the 24-h period in AL rats, even though ZT18 resulted higher than ZT3. The Cosinor test failed to detect a daily rhythm, but he RF group resulted rhythmic as in the case of the SERCA activity. This result was confirmed by the Cosinor test (Table 2). A significant difference (P < 0.001; Student t-test) between the Fasted and RF-ZT3 group was also observed (Figure 5B).
Liver expresses the PMCA isoforms 1 and 4 . Since the enzymatic assay did not distinguish among PMCA isoforms, a Western blot for isoforms 1 and 4 was done. The PMCA1 protein did not oscillate in the AL group, whereas the RF group showed a significant daily rhythmicity (Figure 5C). The PMCA1 acrophase was at the end of the light period (Table 2). A similar pattern was seen for PMCA4: AL did not present a daily rhythmicity and RF groups showed significant differences with time (Figure 5E). Interestingly, the Acrophase of PMCA4 in the RF group was delayed by ~ 5 h (Table 2). However, the protein level of neither PMCA1 nor PMCA4 was affected by alimentary condition (Fasted vs Refed) (Figure 5D and F).
Rhythmical analysis of the hepatic calcium handling channels and ATPases during restricted feeding / FEO expression
The cosinor test helps to determine a potential 24-h rhythm that is consistent with a sinusoidal curve. The rhythms for all the proteins studied in this project were confirmed first by a one-way ANOVA prior to the 24-h cosinor test (Table 1). Consistent with the corticosterone and Per1 levels measured in rats expressing the FEO (Figure 1), their acrophases were clearly shifted around the time of food access (ZT4 and ZT6) (Table 1). A 24-h rhythm was detected for IP3R types 1 and 2, RyR, and PMCA types 1 and 4, either with AL or RF protocols or with both AL and RF protocols (P < 0.05, Table 2).
Effect of restricted feeding/FEO expression on the hepatic zonal distribution of calcium-handling proteins
To compare the changes in the binding and activities of the hepatic calcium channels and ATPases during FEO expression in the specialized hepatocyte populations, we analyzed for the presence of these proteins within the hepatic acinus, distinguishing among the PP and PC zones (Figure 6). Since the main metabolic and physiological changes in the parameters studied were detected at ZT 3 (before food access in RF protocol) and ZT6 (after feeding in RF protocol), this analysis was done only at these times. The immunohistochemistry of IP3R1 (Figure 6A) showed a homogeneous distribution along the hepatic acinus except in the AL group at ZT3, which showed an increase in the periportal zone, and the Refed control group in which this protein was predominantly in the pericentral zone. The IP3R2 distribution in control rats did not show zonal heterogeneity; in contrast, for RF rats, it was higher in the periportal zone at ZT3 and ZT6 (Figure 6B). Again, the control Refed group showed an increase in IP3R2 in the pericentral zone in comparison to the RF rats at ZT6. The ryanodine receptor was found primarily in the periportal zone in almost all conditions (Figure 6C). RyR was detected in the pericentral zone only at ZT6 in RF rats, as well as in the Refed control group. SERCA2 showed a periportal distribution in the AL group, but not in the RF or control feeding groups (Figure 6D). PMCA1 occurred in the perportal zone in the RF group at ZT6, but was uniformly distributed throughout the hepatic acinus in all other groups (Figure 6E). PMCA4 distribution also showed no zonal preference, unlike the other calcium-handling proteins (Figure 6F).
Many reports indicate that food entrainment profoundly affects liver physiology, from gene expression studied by microarrays  to a variety of metabolic adaptations [11,13,14,17,18]. In this context, a report from Stokkan et al.,  indicates that the peak of Per1, measured by bioluminescent luciferase activity, was modified to coincide with the time of food access. A similar result was reproduced in Figure 1, indicating a successful entrainment to mealtime in the RF group. A similar entrainment was also evident in the daily corticosterone daily profile, since the RF group showed 2 peaks (one corresponding to the FFA and the other to the light–dark transition), whereas the AL group showed only the second peak. This result suggests that even though the daytime food cue strongly modifies the timing system, the influence of the SCN is still present during the FEO expression.
A role for intracellular calcium dynamics in regulating the circadian clock has been reported in several models . This secondary messenger is involved in the entrainment process  and the 24-h rhythmicity of clock genes  as well as in the components of the circadian output . For example, RyR2 seems to participate in the clock machinery function of the SCN by modulating the membrane potential . However, so far no reports exist regarding the daily rhythm characterization of the liver calcium-handling proteins in the protocol of daytime restricted feeding/FEO expression. A previous report, using a pharmacological approach in hepatic explants, showed that IP3R, RyR, and SERCA modulated the clock gene Per1 rhythmicity .
This project demonstrated that the activity and protein levels of the hepatic calcium releasing channels IP3R and RyR and the calcium pump SERCA have a daily rhythm that, in addition, is modified by changes in meal access time. It is known that IP3R types 1 and 3 as well as the RyR type 2 oscillate with a circadian periodicity in the SCN [49,53]. In addition, IP3R has been suggested as an element that contributes to the SCN entrainment mechanism . In liver tissue, the IP3R has an important role regulating the calcium oscillations involved in metabolic processes such as bile production, mitochondrial activity, and the gluconeogenic pathway [54,55]. It was reported that the period of the clock gene Per1 was lengthened by inhibiting the IP3R using 2-APB (2-Aminoethoxydiphenyl borate) in liver explants . These data strongly suggest that calcium handled by the IP3R participates in regulating the hepatic molecular clock rhythmicity. The present project showed that the levels of hepatic IP3R protein and binding activity displayed robust daily rhythms in the AL group, but as a consequence of daytime food restriction, their daily fluctuations are modified. The changes in the properties of the IP3Rs were also observed in the control groups of feeding conditions (Fasted and Refed), which suggests that the IP3Rs are under circadian and metabolic control. In food restriction, the circadian control could be exerted by the SCN as well as by the FEO. The IP3R has been postulated to be a positive modulator of glyconeogenisis in the liver as an element responsive to glucagon stimulus . Our group has previously reported daily fluctuations of glucagon in RF conditions that show good correspondence to the daily variations of IP3Rs observed in this report (both showed peaks in the activity and protein levels of IP3R at the time when animals are expecting their meal) . In addition, our data of daily variations of IP3R properties during restricted feeding showed good coincidence with the 24-h changes observed in gluconeogenic enzymes (manuscript in preparation).
The hepatic RyR has been related to mitochondrial respiratory activity, glycogen catabolism, liver regeneration, and protection against hypoxic stress [27,56,57]. The isoform of RyR within the liver is truncated, but it still shows the pharmacological profile of the better-studied skeletal muscle RyR isoform . At 1 mM, ryanodine inhibited RyR and lengthened the circadian period of Per1-luc expression in liver explants of rats feeding ad libitum. In contrast, liver explants from animals under restricted meal schedule did not show modification in the PER1::LUC rhythmicity. Our data indicated a clear daily rhythm in RyR protein levels and ligand-binding properties in the AL group. Interestingly, the RF group showed similar patterns, which supports the finding in liver explants , and suggests that RyR is not affected by the RF protocol. These data seem to indicate that the role played by the RyR in the timing system of the liver differs from the one in the SCN .
It is well accepted that the calcium ATPases SERCA and PMCA function as enzymes that maintain calcium homeostasis by intruding this cation to the ER or extruding it to the extracellular space, respectively. Both calcium pumps modulate the frequency of [Ca2+] oscillation in response to hepatic endocrine stimuli . When SERCA or PMCA are overexpressed, expression of the other decreases, suggesting a finely tuned communication between the two enzymes . SERCA is also involved in the ER stress response as well as in regulating lipid metabolism. Overexpression of hepatic SERCA2 is associated with a significant rise in lipogenic activity, whereas at the same time, it mitigates ER stress in obese mice . However, there are no reports studying daily fluctuations of these ATPases. One of the most remarkable results of this study is the large increases in SERCA protein and activity observed in the RF group before and after the feeding time (from ZT0 to ZT9), which is coincident with the onset of FAA (before meal) and the intense hyperphagic event (after mealtime). A possible interpretation of elevated SERCA activity during the restricted feeding/FEO expression is that high calcium within the ER lumen and the low cytoplasmic calcium could be needed in the reticular response in preparation for food processing and in the lipolytic response during nutrient handling .
Regarding the daily patterns characterized in the hepatic calcium-release channels and pumps, almost all the proteins that had 24-h rhythmicity, also showed a shift in their Acrophase (Figure 7). Only RyR peaks (activity and protein level) were not totally shifted to meal time in the RF protocol. The RyR binding assay demonstrated a modest 3 h phase advance in the RF group with respect to the AL group (Figure 2E, Table 2), whereas the other calcium-handling proteins showed phase advances with an average of 9 h (Figure 7, Table 2). Since amplitude is another rhythmic parameter that must be taken into account, it is noteworthy that SERCA activity amplitude in RF rats increased almost 4 fold in comparison to the AL group (Figure 4A, Table 2). MESOR was no exception, since changes in both IP3R isoforms (a 50% decrease) and SERCA activity (~2-fold increase) during RF protocol/FEO expression were detected. The temporal profile of PMCA activity did not show a 24-h rhythm according to cosinor for both the AL and RF group. However, an ultradian 8-h rhythm was detected for this calcium pump activity in the AL group (data not shown); meanwhile a 24-h period rhythm was observed for both PMCA isoforms protein levels in the case of the RF group (Table 2).
Although IP3R2 zonal distribution within the hepatic acinus has been described , those for RyR, SERCA, and PMCA were studied for the first time, including the possible influence by feeding or chronobiological factors. Using a protocol of liver regeneration, IP3R type 2 as well as hormone receptors responsible for IP3R activation (glucagon and adrenaline receptors) were detected mainly in PP hepatocytes [20,27]. A significant PC distribution was seen for adrenergic receptors . In addition, it was tested that hepatic calcium waves induced by the activation of the vassopresin receptor start in the PP and then disperse towards the PC zone . Circadian variation of glycogen deposition was reported by histological studies for both zones. The proteins studied in our project differed in amount and activity between fasted and ZT3 or between re-fed and ZT6, even though they are under similar feeding conditions; this can be explained in the context in which FEO expression (induced by the RF condition) produces characteristic rheostatic adaptations in liver physiology . This could be related with glycogen breakdown, we found that glycogen is more abundant at ZT3 than after a single 21-h fast, which confirms that circadian meal entrainment establishes a new physiological state in hepatic tissue .
Even though our data did not show a direct link between modifications in calcium-handling proteins and intracellular calcium dynamics associated with food entrainment, our results indicate a clear chronostatic adaptation for the liver calcium channels and ATPases considered in this study. Some authors have suggested a possible role of CamKII (Calmodulin Kinase II) as mediator of the coupling between intracellular calcium and circadian oscillations . More experiments are needed to define if the calcium-handling proteins have an influence on: 1) the molecular clock of the liver, 2) the synchronization process, or 3) the output response of the liver oscillator.
As summarized in Figure 7, the activity, presence, and zonal distribution of liver IP3Rs, RyR, SERCA, and PMCA were modified differentially in ad libitum and daytime restricted feeding protocols. Whereas IP3R and SERCA changed during preprandial time, RyR and PMCA modifications were largely postprandial. The phase of the peaks in the daily rhythms of IP3Rs, RyR, and PMCA1 in the AL group was coincident with the active time of the rats, whereas the temporal pattern of activity and protein levels of IP3Rs, SERCA, and PMCA1-4 changed toward the meal access schedule in the animals under RF. These data suggest that during the reosthatic or allostatic adaptations shown by the liver associated with FEO expression, there exists a differential coordination of these calcium-handling elements which could have a direct impact on hepatic physiology and metabolism.
Our results indicate that a daily rhythmic regulation occurs in the calcium-handling proteins, and very likely in hepatic calcium signaling. The RF condition promotes an adjustment in the activity, protein level, and zonal distribution of these calcium channels and ATPases (Figure 7). Hence, these results demonstrate that important elements in the intracellular calcium dynamics of the liver exhibit daily variations in the control condition of ad libitum feeding. This rhythmicity is further modulated during a protocol of daytime RF and the concomitant FEO expression.
The authors declare that they have no conflicts of interest, financial or otherwise.
AB-R, KC-G and OV-M performed experiments; AB-R, RA-R and MD-M analyzed data; AB-R prepared figures; AB-R, RA-R and MD-M conceived and designed the research; AB-R, RA-R and MD-M interpreted the results; AB-R and MD-M edited and revised manuscript. All authors approved read and approved the final manuscript.