Abstract
Introduction:
Leptin, primarily secreted by adipose tissue, is a critical hormone involved in regulating energy balance and food intake by inducing satiety. Although several hormones influence satiety, leptin plays a dominant role in long-term satiety regulation.
Methods:
We apply a state-space estimation framework using Bayesian filtering to infer continuous, long-term pro-satiety states from plasma leptin concentrations collected from premenopausal women with obesity. Our approach adopts methodologies previously applied to biosignals such as skin conductance and cortisol data to estimate latent states, leveraging the features in the leptin secretory pulses and plasma leptin levels. Additionally, we investigate the potential influence of meals, sleep, and bromocriptine treatment on the pro-satiety states. We introduce the High Satiety Index (HSI), a direct, long-term satiety measure based on leptin secretion dynamics, minimizing biases inherent in conventional assessment methods.
Results:
Comparisons of the estimated state in different time windows show that the pro-satiety state inferred by leptin secretion is significantly higher during sleep, aligning with a circadian rhythm. The estimated state does not show a significant variation in response to meal intake or bromocriptine treatment.
Discussion:
The leptin-based estimator reflects basal variations of satiety in women with obesity. This framework shows the feasibility of applying Bayesian filtering to track satiety and can be further developed to perform a multimodal estimation.
1 Introduction
Leptin is a hormone secreted primarily in the adipose tissues (fat cells) and is associated with body weight (, ). It regulates body energy and food intake by inducing satiety (, –). Although other hormones such as gut hormones also affect satiety, leptin is a major factor that affects long-term satiety (, ). Leptin provides information on the energy storage in the human body to the brain. As fat accumulates, more leptin is released and acts in the hypothalamus, leading to high satiety and behaviors that will reduce energy intake (). Previous studies have conducted and discussed several methods of satiety measurement, including bipolar and unipolar scales () and appetite ratings using the visual analogue scale (VAS) (, ). However, these measurements are subjective and may require consideration of various conditions under which the measurements are taken (). The measurement of satiety is often under the guideline of the satiety cascade, first introduced by Blundell et al. in (). The cascade describes short-term and long-term aspects that affect satiety. It also distinguishes short-term and long-term satiety by naming the former “satiation” and the latter “satiety” (). The former refers to the feeling that stops eating, meaning that it is mainly affected by sensory factors such as gastric distension (). Satiety, or more precisely long-term satiety, on the other hand, determines the beginning of meals and is often considered the opposite of hunger (). Long-term satiety is often measured between meals, but the measurement is usually complicated by many internal and external factors ().
Other hormones are also involved in the satiety cascade, such as ghrelin, a hunger-inducing hormone (, ). Ghrelin is secreted from the stomach and has a short-term effect of increasing appetite (). It does not affect leptin but acts as a counterpart to leptin and can be seen as an anti-satiety feature (), as they both act in the hypothalamus and affect downstream neural signaling that leads to food related behaviors (). Leptin’s primary role is to regulate energy storage in the body, but its regulation is asymmetric, as it is more effective in replenishing fat storage when it is low than resisting obesity (). Giving leptin’s characteristics, our goal is to decode a continuous pro-satiety state from plasma leptin concentration collected from a cohort of women with obesity, as an effort to measure the long-term basal satiety in an objective way.
Many latent state estimation schemes have been developed based on state-space models. Information extracted from biosignals and behavior measurements has been used as observations from which the state is estimated (). These state-space methods have been applied in various fields including behavioral learning (–), movement decoding (, ), latent state underlying neural spiking estimation (, ), and sleep onset process tracking (). The estimators developed in these studies often utilized Bayesian filtering for processing binary observations. Moreover, Wickramasuriya et al. () introduced a marked point process (MPP) plus continuous observation estimator that takes the binary sequence of impulse occurrences underlying sweat secretion as one of the observations to estimate the hidden sympathetic arousal state. In this research, we implement a similar state-space approach, leveraging information extracted from plasma leptin measurements and utilizing Bayesian filtering.
Similar to skin conductance signals that arise from sweat secretion by sweat glands near the measurement site, blood leptin comes from fat cells () and can be modeled by a two-compartment state-space model (). In addition, leptin pulse amplitudes and inter-pulse intervals have similar distributions to sweat secretory pulses (, ). The estimator applies Bayesian causal inference which is well-suited for decoding the latent arousal state from skin conductance and latent energy state from cortisol given the causal links (). But this estimator cannot directly decode long-term satiety from leptin data. Unlike the relationship between the arousal state and sweat secretion, where the hidden state causes the secretion, satiety is an effect of leptin (), which means that the direction of causality is the opposite. However, we assume a latent pro-satiety state that controls leptin secretion, and this pro-satiety state can be an indicator of a person’s overall satiety. Since the state estimation scheme assumes a causal link between the hidden state and leptin secretion, and leptin controls long-term satiety (, ), the estimated state based on leptin may display the dynamics of a person’s basal satiety over time. Figure 1 illustrates our proposed satiety inference scheme.
Figure 1
After state estimation, to investigate the latent state’s potential association with food intake, we compare the estimated state before, during, and after each meal. We also compare wake and sleep periods since it has been reported that leptin follows a circadian rhythm (, ). The original study of the dataset used in this study was designed partly to test the effect of bromocriptine on obesity (). With this satiety inference framework, we also examine the effect of bromocriptine treatment on the estimated pro-satiety state during these mentioned windows, providing a new perspective in testing the drug’s effect. We define a High Satiety Index (HSI) that provides a judgment of the participant’s basal satiety status based on the estimated state. With this objective method of satiety evaluation, continuous tracking and automated control systems can be implemented in portable devices to monitor and influence a person’s satiety. It is also a step toward a multimodal method of satiety inference and a monitor for dynamic latent health status in the human body.
2 Materials and methods
2.1 Dataset
We use leptin data collected in the previous clinical studies by Kok et al. (–). The original clinical study was approved by the Medical Ethics Committee of the Leiden University Medical Center and conducted based on the Helsinki Declaration (). The written acknowledgment of informed consent for participation was obtained from all participants. The participants were 18 healthy premenopausal women with obesity (BMI 30.1–40.5kg/m2, age 22–51 years with mean = 37.5 ± 1.7) who were not under the influence of any medication or drugs. Participants with acute or chronic disease, depression, head trauma, habits of smoking or alcohol consumption, recent transmeridian flights, nightshift work, weight change, blood donation, or participation in another clinical trial were excluded. The data were gathered during the early follicular phase of their menstrual cycles confirmed by plasma estradiol and progesterone levels, to minimize hormonal variability between individuals and across measurements. All, and all participants were served with the same meal and had a strict meal and sleep schedule. The meals consist of liquids (16% proteins, 49% carbohydrates, and 35% fat) for a total of 2100 kilocalories per day. The meals periods are from 09:30 to 10:00, 13:00 to 13:30, and 18:30 to 19:00. The sleep period is from 23:00 to 07:30 on the next day. The trial is 24 hours long, from 09:00 to 09:00 on the next day, following 7 days of placebo treatment. Plasma leptin concentration levels were measured from blood samples using radioimmunoassay, which has a detection limit of 0.5 ng/L, at 10-minute intervals throughout this period. The same protocol was repeated four weeks later with the same participants who received a 2.5 mg dose of bromocriptine twice per day for 7 days this time.
In their study, Reddy et al. () extracted the underlying secretory pulses as an MPP through a deconvolution method based on a state-space model. They also reconstructed the leptin level using the MPP. In this study, we utilize both the extracted MPP and the reconstructed leptin level as inputs to the estimator.
2.2 Leptin-related pro-satiety state estimation
To estimate the hidden pro-satiety states that give rise to leptin secretory events, we employ an MPP estimator with continuous observation that involves the inter-pulse intervals, pulse amplitudes, and the plasma leptin level altogether. Unlike a binary estimator that treats the secretory pulse sequence as a point process where pulse amplitudes are ignored () or an MPP estimator that estimates the hidden state based on the MPP only (), this approach assumes linear relationships between both the hidden state and the MPP and between the state and the plasma leptin level. In their previous study, Wickramasuriya et al. applied such an estimator to estimate arousal state based on skin conductance and energy state based on blood cortisol level (). The leptin pulses studied in this paper were recovered based on the same state-space model structure as applied in the estimation of skin conductance and cortisol pulses (, ). Since there are very few pulse events (29 to 45) compared to the total length of 1441 time instants where a pulse can occur, a substantial portion of the MPP is 0 which trivializes the nonzero amplitudes (). To solve this, we chose larger time bin sizes based on the minimum time interval between pulse events for each MPP and resampled the data. We use the same state dynamics and algorithm to decode leptin secretion. The state dynamics is:
where xj is the state at time instant j and is the process noise. This definition assumes a random walk process that controls the state’s evolution over time, which was chosen to model neural states in several previous works as well as the skin conductance study (–, , ). We denote the occurrence of a pulse in the MPP as nj∈ (0,1). njis a Bernoulli-distributed random variable with mass density . Let , then pj is related to the state variable xjvia a logit transformation according to the theory of generalized linear models ():
where β0 and β1 are coefficients to be determined.
We further define the linear relationship between the leptin secretory pulse amplitudes rj and the state variable xj as:
where γ0 and γ1 are coefficients to be determined, and is the measurement noise ().
Based on the binary variable njand the continuous variable rj the joint probability density function for an observation at a time instant given the hidden state is,
This is based on the assumption that the distribution of the amplitudes of the leptin secretory pulses can be modeled by a Gaussian model, which is proven to be effective in our previous study (). The Gaussian distribution is also assumed to be a proper model for the distribution of the skin conductance pulses (). When no pulse is observed, rj is not involved.
In addition to rj, we relate the reconstructed plasma leptin level sj to the state variable xjvia:
where δ0 and δ1 are coefficients to be estimated, and is the measurement or modeling noise (). The original continuous-valued leptin level from plasma measurements has a lower resolution (10-minute time bin size) than the extracted underlying pulses (1-minute time bin size). To address this problem, we reconstruct the continuous-valued leptin level using the extracted pulses. The reconstructed results according to the state-space model of leptin secretion defined in the previous chapter closely resemble the original measurement (). Now the state xj gives rise to both the MPP observation and the plasma leptin level observation with 1-minute time bin size. Although xj is linearly related to both the MPP and the reconstructed leptin level, the overall model is nonlinear due to xj’s relationship with pj ().
To estimate the hidden state xj and , we apply an expectation-maximization (EM) algorithm. The EM algorithm consists of two steps: E-step, where we predict the current state and the process noise limit based on the initial or previous state prediction; and M-step, where we estimate the parameters that maximize the likelihood of the predicted state variable. Based on Equations 1–6, the posterior density function for the E-step is ():
where yj denotes the collective observations of the MPP and the reconstructed leptin level.
2.2.1 E-Step
Based on Equation 7, we implement Bayesian filtering and utilize both a forward filter and a backward smoother in the E-step, as applied in the skin conductance study (, , ). The forward filter estimates using all previous observations , and the backward smoother estimates using all observations (). The following shows one iteration of the E-step which involves a Gaussian approximation (, ).
Predict:
Update:
if nj = 0
if nj= 1
Note that since pj|j is calculated using xj|j (), and () has xj|j on both sides of the equation which is solved numerically by the Newton-Raphson method. Besides, the algorithm switches between two ways of updating based on nj’s value. This switching is to account for the different cases in () for the MPP.
After predicting the state variable and the variance at each time bin, we improve the predictions by applying the backward smoother (, , , ):
2.2.2 M-Step
We derive by maximizing the expected log-likelihood L defined as
where denotes the sequence of locations where pulses occur. To estimate β0 and β1, we apply one of the two approaches mentioned in (). From Equation 2, we have
if xj = 0 at time j. We set β1 = 1 and assume that x0 ≈ 0, and the baseline probability pb approximates p0. pb is taken to be the number of pulses in the MPP divided by the length of the MPP:
The algorithm based on Equations 8–21 iterates between the two steps until convergence, which we consider is achieved when the mean distance between obtained variables in consecutive iterations is smaller than 10−8 (, , , ).
2.3 High satiety index
Satiety index has been given to food, suggesting how much satiety different types of food can induce in a person (). Inspired by the High Arousal Index defined in (), we define the High Satiety Index (HSI) based on the estimated pro-satiety state in humans as mentioned before. The HSI is expressed as the probability of the estimated state exceeding a certain threshold. It indicates the certainty of the participant having high satiety based on observation, similar to the ideal observer certainty defined in (). We set the threshold to be the median of the state value so that it may illustrate a potentially significant change from wakefulness to sleep period during a full day. Since the estimated state reflects long-term satiety, the HSI indicates the baseline change of the participant’s satiety. When the HSI is near 1, it suggests that the participant is very likely satiated at that time and has no appetite, whereas when it is near 0, the participant potentially wants to start eating. It is worth noting that the HSI might not reflect the direct feeling of hunger or fullness which are more related to short-term satiety or satiation ().
3 Results
We applied an MPP with a continuous observation estimator to uncover the leptin-related hidden state of the 18 participants. The estimator assumed causal links both between the hidden state and the leptin secretory pulses and between the hidden state and the reconstructed plasma leptin levels. Through the EM algorithm, we obtained the estimated pro-satiety state over 24 hours for the pre-treatment and post-treatment samples of each participant. For comparison, we also modified the estimator to apply two other versions. One only considers the MPP and the other only considers the continuous observation. The solely MPP-based estimator was implemented in (), where the authors linked the pulsatile sweat secretion modeled as an MPP to the sympathetic arousal state. The estimator design is the similar to the estimator described above where a logit transformation is applied to linearly link the pulse probability to the state variable. Similarly the other version that is solely based on the continuous observation considers the relationship in 6 alone, thus modifying the update equations and cost function in the EM algorithm. The overall structure remains the same in both versions and their results are compared to the main estimator, which is based on both MPP and continuous observation.
The results show that the estimated state based on MPP alone does not vary much over time compared to the other two versions while having relatively wide confidence intervals. An example is shown in Figure 2. From the solely MPP-based estimator, the estimated probability of pulse occurrence also has little oscillations through our the day, while the HSI stays close to the middle for the majority of the time. On the other hand, the estimator solely based on the continuous observation overfits the data. At the same time, the MPP plus continuous observation estimator balances the variances of the MPP and the leptin levels to output a meaningful result that reflects changes in both inputs and avoids overfitting.
Figure 2
Since no pulse sequence is involved in the continuous observation estimator, pulse occurrence probability is not meaningful; likewise, due to the lack of a relationship between the leptin levels and the hidden state, we cannot reconstruct leptin levels from the state estimates in the MPP estimator. A comparison of the results produced by these three versions of estimators for an example participant is shown in Figure 2. Results for all participants are available in the Supplementary Materials.
In the results of the MPP plus continuous estimator, we saw a dominant pattern of increased plasma leptin level during sleep in all samples, regardless of the bromocriptine treatment () as shown by Participant 6’s state estimation result in Figure 3. The estimated state shows a general trend that consists of an increase during the day and a rapid decrease before waking, with a peak around mid-sleep and a nadir in the morning.
Figure 3
This is also reflected by the HSI in the majority of cases. Example estimation results are shown in Figure 3. Estimation results of other participants can be found in the Supplementary Materials.
To evaluate potential meal influence on the pro-satiety state, we calculated the mean state value during different meal-related time periods. Figure 4 shows the distribution box plots of the mean state values of all participants. Specifically, in these time windows, a participant might be preparing to eat, eating, or digesting food. Comparing the estimated state during these time windows may inform pro-satiety state changes due to having meals. According to the trial design that includes a strict schedule for meals, we focus on 4 time windows for each meal:
29 minutes to 0 minutes before the start of the meal period (-30 window),
the 30-minute meal period (meal window),
1 minute to 30 minutes after the end of the meal period (+30 window), and
31 minutes to 60 minutes after the end of the meal period (+60 window).
Figure 4
We expect that the state values will not show significant change during each meal-related window than in all previous windows, since our estimator is solely leptin-based which is not sensitive to food intake in a short period (). To test if there is an increase or decrease and if the increase or decrease is consistent during and one hour after meals, we made six comparisons with the median state value for each meal-related time window. Specifically, we perform a two-tailed Wilcoxon signed rank test () on each pair of data to test if they are different with statistical significance.
According to Table 1 and Figure 4, there is a consistent increase for at least one hour since the beginning of lunch. The distribution of state values during the lunch related windows seem to spread out more across all participants after bromocriptine treatment, but the general increase remained statistically significant. Besides, both pre-treatment and post-treatment data show a much higher mean state value during sleep than during wakefulness (p = 0.0002 and p = 0.0033). Wilcoxon signed rank tests comparing pre-treatment and post-treatment data over the same time periods show no statistical significant results.
Table 1
| Meals | meal vs.-30 | +30 vs.-30 | +60 vs.-30 | +30 vs.meal | +60 vs.meal | +60 vs. +30 |
|---|---|---|---|---|---|---|
| Breakfast pre-treatment | 0.3491 | 0.3061 | 0.3720 | 0.6165 | 0.5566 | 0.3271 |
| Breakfast post-treatment | 0.1570 | 0.1701 | 0.1989 | 0.1701 | 0.2668 | 0.5566 |
| Lunch pre-treatment | 0.0096 | 0.0123 | 0.0043 | 0.0123 | 0.0005 | 0.0002 |
| Lunch post-treatment | 0.0249 | 0.0004 | 0.0002 | 0.0002 | 0.0002 | 0.0002 |
| Dinner pre-treatment | 0.0347 | 0.1119 | 0.1570 | 0.3958 | 0.3491 | 0.2311 |
| Dinner post-treatment | 0.1989 | 0.2145 | 0.0854 | 0.2145 | 0.0347 | 0.0429 |
P-values of two-tailed Wilcoxon signed rank test on estimated state from each two meal-related windows.
meal: meal window; −30: from 29 minutes before to the start of a meal; +30: 1 to 30 minutes following a meal; +60: 31 to 60 minutes following a meal. p-values less than 0.05 are in bold (indicating statistical significance).
4 Discussion
Using the MPP plus continuous observation estimator, we obtain the hidden pro-satiety state that leads to the observation of both the underlying leptin secretory pulses and the plasma leptin levels. We also applied an MPP-based estimator without linking the state to the reconstructed leptin levels. As demonstrated in Figure 2, this estimator does not produce state values that have large variations over time compared to the results from the other two estimators. This is likely due to the fact that leptin pulses are more spread out throughout the day and do not have a burst during any short time window. To strengthen the influence of the leptin secretory events, we downsampled the pulse sequences by dividing them into segments and using a single pulse to represent all pulses in each segment with an aggregated pulse amplitude. This helps concentrate the influence of each pulse. However, since almost no inter-pulse intervals are very small (), the estimated state is still relatively flat after downsampling the pulse sequence. For estimating the pro-satiety state, adding the continuous leptin concentration as a reference would make the estimated state more informative, as it would show additional features which might directly reflect the impact of leptin in the hypothalamus, leading to downstream satiety related neural activities ().
A simple state estimator based on continuous observation alone can extract important trends from the data. However, this estimator leaves out the influence of the timings and amplitudes of leptin secretory pulses, and thus cannot produce meaningful probabilities for the pulse occurrences. Algorithmically, since the estimator also updates and minimizes the variances of the input, it quickly overfits and presents a convergence issue. The formulation of the estimator involves a forward filter and a backward filter. Without the MPP, the forward filter resembles a Kalman filter that reduces noise based on the variance of the input (, ). However, since the estimator also minimizes the input variance, it overfits the input as the variance approaches 0. On the other hand, the output of the MPP plus continuous observation estimator outputs a state that reflects changes in the pulse sequences and the leptin concentration levels, balancing features from both inputs.
The estimated leptin-related hidden states capture the general trend of the reconstructed plasma leptin level which displays a peak around midnight and decreases rapidly until reaching a minimum in the morning. The HSI also shows a clear distinction between wake and sleep in most results. However, some cases such as the estimated state of the pre-treatment sample of participant 16 (shown in Figure 3) show more oscillation than the average results. And the HSI also shows a relatively chaotic pattern. A disrupted circadian rhythm of the biological system might have been a cause of the participant’s obesity (). Despite these cases, our comparison and test results show that the wake-sleep difference is statistically significant, and in most cases, the estimated states and the HSI have patterns similar to the results of participant 6 shown in Figure 3, suggesting a circadian rhythm (). As expected, this circadian rhythm resembles that of blood leptin levels (, ), but it is more apparent. The near-zero HSI during most of the day in most cases reflects low basal satiety, meaning that the participant may have a good appetite under the right conditions (). On the other hand, the high HSI during nighttime suggests high basal satiety in the participant who is unlikely to feel the need for food (). Moreover, leptin as a satiety hormone is not only regulated by wake and sleep periods but also active and resting periods (). The pro-satiety state also reflects this behavior. In many cases such as participant 6, both the pre-treatment and post-treatment samples indicate a high satiety after dinner and before sleep. This might suggest that the participants are relatively inactive though still awake (). The HSI derived from the pro-satiety state reveals basal satiety status that can hardly be accurately measured through self-reporting especially in a free-living setting due to underreporting and misreporting ().
By comparing the differences between time windows around each meal, we have found a consistent increase in the estimated state during the time windows after lunch across all participants. However, such a pattern is not apparent after dinner and does not show after breakfast. Further studies that focus on meal sizes and compositions might provide us with more details on the reason behind the discrepancy between the effects of the different meals (). Nonetheless, circadian rhythm is dominant and can overshadow the effect of meals on long-term satiety if there is any (). The decrease in the pro-satiety state before waking could extend to the morning and the breakfast period. The low or slowly rising HSI during these time windows further supports that food intake does not have an acute effect on the leptin-based pro-satiety state.
The original study on the dataset has found a decrease in overall mean leptin level after bromocriptine treatment (). We expanded on this comparison by comparing the pro-satiety states estimated from pre-treatment and post-treatment data different time windows for each participant. However, the results did not show a significant change during any of the investigated time windows. Although it is shown in the previous studies (, ) that the mean leptin level is reduced after the treatment, bromocriptine does not affect the diurnal pattern of the pro-satiety state. Future studies including real-world monitoring studies may also implement this method in testing the effects of other medication on satiety. The proposed dynamic inference of satiety provides more details than traditional tests on related hormones or surveys.
Feeling satiated involves an integration of hormonal signals in the brain, including ghrelin secreted from the stomach and insulin from the pancreas (, ). Insulin has a diurnal pattern exactly in phase with that of leptin while it is more sensitive to food intake, whereas ghrelin has a pattern reciprocal to that of insulin (, ). This study focused on leptin and its link to the long-term aspect of the satiety cascade (). Based on our satiety inference scheme, a multimodal estimator can be developed to incorporate more related biosignals such as insulin, and ghrelin, which have short reaction times to food intake (). Nonetheless, using an unimodal approach, we isolate the hormone and may have a clearer understanding of its specific effect on human homeostasis. To fully capture the satiety felt by a person, downstream neural links need to be studied. In other words, how leptin and other hormones signaling in the hypothalamus leads to behaviors such as food seeking, food consumption, and meal termination requires further investigation (). Both a unimodal analysis of each involved biosignals and a multimodal analysis of the biosignals combined might be needed for determining their exact roles in appetite control.
The presented unimodal algorithm may have better time efficiency when implemented in real-world applications like health monitoring devices. The leptin data collection required measuring plasma leptin concentrations from blood samples every 10 minutes for 24 hours (). The details are presented in the Methods section. This invasive way of data collection is not practical for long-term monitoring and the development of an automated wearable device. To increase temporal resolution, a higher frequency of sample collection might be preferred. Thus, non-invasive measurements are essential for a potential wearable device. Recently, studies have applied new technologies to measure and monitor glucose, glucocorticoids, and other biosignals in non-invasive or minimally invasive ways, such as electrochemical continuous glucose monitoring (CGM) () and sweat-sensor based CGM (). Other CGM schemes include using optics or microwaves (). In addition, U-RHYTHM, a portable device that collects interstitial fluid is capable of measuring glucocorticoids and other hormones for a long period (, ). These examples have shown that a better leptin measuring scheme is possible, and implementing our method in a wearable device is achievable. Since our long-term satiety estimator assumed relatively simple links between observations and the latent state, the estimation process is fast and can be done in real time. However, the deconvolution for recovering the underlying pulses may take a longer time, but a near real-time deconvolution algorithm is also possible (). In addition, the presented state estimation algorithm is flexible and can be customized depending on the individual by setting different hyperparameters (details on hyperparameters are available in the Supplementary Materials). This could be a crucial function for wearable devices that can adapt to users with different health conditions.
In conclusion, this study applied an estimator to infer the 24-hour basal satiety using leptin secretory pulses and the continuous leptin levels. The estimated hidden states display a circadian rhythm with a prominent increase during nighttime. This hidden state is an indicator of long-term basal satiety grounded in the assumed causal link between the hidden state and leptin secretion as well as between leptin secretion and the pro-satiety state (, ). Despite limited sample size, our method presents the feasibility of applying a Bayesian framework to track satiety, demonstrating its practicality in revealing person-specific hidden health status and providing a new direction in drug testing and patient monitoring for future studies with larger cohort sizes. Future studies may also incorporate other satiety-related hormones such as ghrelin and insulin and apply both a unimodal and a multimodal framework to reflect a more complete spectrum of satiety. Furthermore, the same method can be applied to investigate the differences in the pro-satiety state and HSI across different weight groups. The results might help design a more personalized treatment of satiety-related diseases and diets for general health improvement. larger groups of participants with different body weights should be studied under the same framework to consider its correspondence with body weight. Last but not least, with the rapid development of minimally invasive and portable devices, the proposed satiety inference method contributes to a real-world implementation of automated monitoring and medication testing for satiety and potentially other latent states.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by Medical Ethics Committee of the Leiden University Medical Center. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
QX: Writing – original draft, Formal Analysis, Software, Visualization, Writing – review & editing, Data curation. SK: Data curation, Validation, Software, Writing – review & editing, Methodology. RF: Methodology, Supervision, Conceptualization, Writing – review & editing, Project administration, Funding acquisition.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This work was supported in part by National Institutes of Health (NIH) grant R35GM151353 - MESH: Multimodal Estimators for Sensing Health and in part by New York University (NYU) start-up funds to R. T. Faghih.
Acknowledgments
The authors thank Dr. Daniel J. Vis for sharing the leptin data used in this article.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2025.1638568/full#supplementary-material
References
1
IzquierdoAGCrujeirasABCasanuevaFFCarreiraMC. Leptin, obesity, and leptin resistance: where are we 25 years later? Nutrients. (2019) 11:2704. doi: 10.3390/nu11112704
2
LandechoMFTueroCValentíVBilbaoIde la HigueraMFrühbeckG. Relevance of leptin and other adipokines in obesity-associated cardiovascular risk. Nutrients. (2019) 11:2664. doi: 10.3390/nu11112664
3
KolaczynskiJWConsidineRVOhannesianJMarcoCOpentanovaINyceMRet al. Responses of leptin to short-term fasting and refeeding in humans: a link with ketogenesis but not ketones themselves. Diabetes. (1996) 45:1511–5. doi: 10.2337/diab.45.11.1511
4
SchoellerDCellaLSinhaMCaroJ. Entrainment of the diurnal rhythm of plasma leptin to meal timing. J Clin Invest. (1997) 100:1882–7. doi: 10.1172/JCI119717
5
HeiniALara-CastroCKirkKConsidineRCaroJWeinsierR. Association of leptin and hungersatiety ratings in obese women. Int J Obes. (1998) 22:1084–7. doi: 10.1038/sj.ijo.0800731
6
BlundellJGoodsonSHalfordJ. Regulation of appetite: role of leptin in signalling systems for drive and satiety. Int J Obes. (2001) 25:S29–34. doi: 10.1038/sj.ijo.0801693
7
BenelamB. Satiation, satiety and their effects on eating behaviour. Nutr Bull. (2009) 34:126–73. doi: 10.1111/j.1467-3010.2009.01753.x
8
AndermannMLLowellBB. Toward a wiring diagram understanding of appetite control. Neuron. (2017) 95:757–78. doi: 10.1016/j.neuron.2017.06.014
9
MantzorosCS. The role of leptin in human obesity and disease: a review of current evidence. Ann Internal Med. (1999) 130:671–80. doi: 10.7326/0003-4819-130-8-199904200-00014
10
MerrillEKramerFCardelloASchutzH. A comparison of satiety measures. (2002) 39:181–3. doi: 10.1006/appe.2002.0496
11
FordeCG. Measuring satiation and satiety. In: Methods in Consumer Research, vol. 2. Amsterdam, Netherlands: Elsevier (2018). p. 151–82.
12
BlundellJEDe GraafKFinlaysonGHalfordJCHetheringtonMKingNet al. Measuring food intake, hunger, satiety and satiation in the laboratory. In: Handbook of assessment methods for eating behaviours and weight-related problems: Measures, theory and research, 2nd ed. Sage, Newbury Park, CA (2009). p. 283–325.
13
BlundellJRogersPHillA. Evaluating the satiating power of foods: implications for acceptance and consumption. In: Food acceptance and nutrition. London, UK: Academic Press (1987). p. 205–19.
14
De GraafCBlomWASmeetsPAStafleuAHendriksHF. Biomarkers of satiation and satiety. Am J Clin Nutr. (2004) 79:946–61. doi: 10.1093/ajcn/79.6.946
15
MüllerTDNogueirasRAndermannMLAndrewsZBAnkerSDArgenteJet al. Ghrelin. Mol Metab. (2015) 4:437–60. doi: 10.1016/j.molmet.2015.03.005
16
SchmidDAHeldKIsingMUhrMWeikelJCSteigerA. Ghrelin stimulates appetite, imagination of food, gh, acth, and cortisol, but does not affect leptin in normal controls. Neuropsychopharmacology. (2005) 30:1187–92. doi: 10.1038/sj.npp.1300670
17
WickramasuriyaDSFaghihRT. A marked point process filtering approach for tracking sympathetic arousal from skin conductance. IEEE Access. (2020) 8:68499–513. doi: 10.1109/Access.6287639
18
SmithACFrankLMWirthSYanikeMHuDKubotaYet al. Dynamic analysis of learning in behavioral experiments. J Neurosci. (2004) 24:447–61. doi: 10.1523/JNEUROSCI.2908-03.2004
19
PrerauMJSmithACEdenUTKubotaYYanikeMSuzukiWet al. Characterizing learning by simultaneous analysis of continuous and binary measures of performance. J Neurophysiol. (2009) 102:3060–72. doi: 10.1152/jn.91251.2008
20
ColemanTPYanikeMSuzukiWABrownEN. A mixed-filter algorithm for dynamically tracking learning from multiple behavioral and neurophysiological measures. In: The dynamic brain: An exploration of neuronal variability and its functional significance. Oxford, UK: Oxford University Press (2011). p. 3–28.
21
DengXFaghihRTBarbieriRPaulkACAsaadWFBrownENet al. (2015). Estimating a dynamic state to relate neural spiking activity to behavioral signals during cognitive tasks, in: 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC) (IEEE), . Milan, Italy: IEEE. pp. 7808–13.
22
ShanechiMMHuRCPowersMWornellGWBrownENWilliamsZM. Neural population partitioning and a concurrent brain-machine interface for sequential motor function. Nat Neurosci. (2012) 15:1715–22. doi: 10.1038/nn.3250
23
ShanechiMMHuRCWilliamsZM. A cortical–spinal prosthesis for targeted limb movement in paralysed primate avatars. Nat Commun. (2014) 5:3237. doi: 10.1038/ncomms4237
24
SmithACBrownEN. Estimating a state-space model from point process observations. Neural Comput. (2003) 15:965–91. doi: 10.1162/089976603765202622
25
SmithACScalonJDWirthSYanikeMSuzukiWABrownEN. State-space algorithms for estimating spike rate functions. Comput Intell Neurosci. (2010) 2010:426539. doi: 10.1155/2010/426539
26
PrerauMJHartnackKEObregon-HenaoGSampsonAMerlinoMGannonKet al. Tracking the sleep onset process: an empirical model of behavioral and physiological dynamics. PloS Comput Biol. (2014) 10:e1003866. doi: 10.1371/journal.pcbi.1003866
27
WickramasuriyaDSKhazaeiSKianiRFaghihRT. A bayesian filtering approach for tracking sympathetic arousal and cortisol-related energy from marked point process and continuous-valued observations. IEEE Access. (2023) 11:137204–47. doi: 10.1109/ACCESS.2023.3334974
28
ReddyRGuoYRajuVFaghihRT. (2023). Characterization of leptin secretion in premenopausal obese women treated with bromocriptine, in: 2023 45th Annual International Conference of the IEEE Engineering in Medicine & Biology Society (EMBC) (IEEE), . Sydney, Australia: IEEE. pp. 1–6.
29
SubramanianSBarbieriRBrownEN. Point process temporal structure characterizes electrodermal activity. Proc Natl Acad Sci. (2020) 117:26422–8. doi: 10.1073/pnas.2004403117
30
XiangQReddyRFaghihRT. Marked point process secretory events statistically characterize leptin pulsatile dynamics. J Endocrine Soc. (2024) 8:bvae149. doi: 10.1210/jendso/bvae149
31
SheaSAHiltonMFOrlovaCAyersRTMantzorosCS. Independent circadian and sleep/wake regulation of adipokines and glucose in humans. J Clin Endocrinol Metab. (2005) 90:2537–44. doi: 10.1210/jc.2004-2232
32
KokPRoelfsemaFFrolichMvan PeltJMeindersAEPijlH. Activation of dopamine d2 receptors lowers circadian leptin concentrations in obese women. J Clin Endocrinol Metab. (2006) 91:3236–40. doi: 10.1210/jc.2005-2529
33
KokPRoelfsemaFFrolichMvan PeltJStokkelMPMeindersAEet al. Activation of dopamine d2 receptors simultaneously ameliorates various metabolic features of obese women. Am J Physiology-Endocrinology Metab. (2006) 291:E1038–43. doi: 10.1152/ajpendo.00567.2005
34
KokPRoelfsemaFFrolichMvan PeltJMeindersAEPijlH. Short-term treatment with bromocriptine improves impaired circadian growth hormone secretion in obese premenopausal women. J Clin Endocrinol Metab. (2008) 93:3455–61. doi: 10.1210/jc.2008-0001
35
WickramasuriyaDSAminMRFaghihRT. Skin conductance as a viable alternative for closing the deep brain stimulation loop in neuropsychiatric disorders. Front Neurosci. (2019) 13:470066. doi: 10.3389/fnins.2019.00780
36
McCullaghPNelderJA. Generalized Linear Models Vol. 37. Boca Raton, FL: CRC Press (1989).
37
MendelJM. Lessons in estimation theory for signal processing, communications, and control. Singapore: Pearson Education Asia / Prentice Hall PTR (1995).
38
HoltSHBrand MillerJCPetoczPFarmakalidisE. A satiety index of common foods. Eur J Clin Nutr. (1995) 49:675–90.
39
WoolsonRF. Wilcoxon signed-rank test. In: Wiley encyclopedia of clinical trials. Hoboken, NJ, USA: John Wiley & Sons (2007). p. 1–3.
40
KimYBangHet al. Introduction to kalman filter and its applications. In: Introduction and Implementations of the Kalman Filter (London, UK: IntechOpen), vol. 1. (2018). p. 1–16.
41
WickramasuriyaDSFaghihRT. Bayesian Filter Design for Computational Medicine: A State-Space Estimation Framework. Cham, Switzerland: Springer Nature (2024).
42
HsuchouHWangYCornelissen-GuillaumeGGKastinAJJangEHalbergFet al. Diminished leptin signaling can alter circadian rhythm of metabolic activity and feeding. J Appl Physiol. (2013) 115:995–1003. doi: 10.1152/japplphysiol.00630.2013
43
PageAJChristieSSymondsELiH. Circadian regulation of appetite and time restricted feeding. Physiol Behav. (2020) 220:112873. doi: 10.1016/j.physbeh.2020.112873
44
SimonCGronfierCSchliengerJBrandenbergerG. Circadian and ultradian variations of leptin in normal man under continuous enteral nutrition: relationship to sleep and body temperature. J Clin Endocrinol Metab. (1998) 83:1893–9. doi: 10.1210/jcem.83.6.4864
45
LangendonkJGPijlHToornvlietACBurggraafJFrolichMSchoemakerRCet al. Circadian rhythm of plasma leptin levels in upper and lower body obese women: influence of body fat distribution and weight loss. J Clin Endocrinol Metab. (1998) 83:1706–12. doi: 10.1210/jcem.83.5.4717
46
DafneVJManuelMARocioCV. Chronobiotics, satiety signaling, and clock gene expression interplay. J Nutr Biochem. (2024) 126:109564. doi: 10.1016/j.jnutbio.2023.109564
47
CummingsDEPurnellJQFrayoRSSchmidovaKWisseBEWeigleDS. A preprandial rise in plasma ghrelin levels suggests a role in meal initiation in humans. Diabetes. (2001) 50:1714–9. doi: 10.2337/diabetes.50.8.1714
48
StoverPJFieldMSAndermannMLBaileyRLBatterhamRLCauffmanEet al. Neurobiology of eating behavior, nutrition, and health. J Internal Med. (2023) 294:582–604. doi: 10.1111/joim.13699
49
ZouYChuZGuoJLiuSMaXGuoJ. Minimally invasive electrochemical continuous glucose monitoring sensors: Recent progress and perspective. Biosensors Bioelectronics. (2023) 225:115103. doi: 10.1016/j.bios.2023.115103
50
ZafarHChannaAJeotiVStojanovićGM. Comprehensive review on wearable sweat-glucose sensors for continuous glucose monitoring. Sensors. (2022) 22:638. doi: 10.3390/s22020638
51
TangLChangSJChenCJLiuJT. Non-invasive blood glucose monitoring technology: a review. Sensors. (2020) 20:6925. doi: 10.3390/s20236925
52
MohanDRLerarioAM. Closing the loop on adrenal health, dysfunction, and disease. Sci Trans Med. (2023) 15:eadh4450. doi: 10.1126/scitranslmed.adh4450
53
UptonTJZavalaEMethliePKämpeOTsagarakisSØksnesMet al. High-resolution daily profiles of tissue adrenal steroids by portable automated collection. Sci Trans Med. (2023) 15:eadg8464. doi: 10.1126/scitranslmed.adg8464
54
AminRFaghihRT. Physiological characterization of electrodermal activity enables scalable near real-time autonomic nervous system activation inference. PLoS Comput Biol. (2022) 18:e1010275. doi: 10.1371/journal.pcbi.1010275
Summary
Keywords
satiety, leptin, obesity, marked point process, Bayesian filtering, circadian rhythm
Citation
Xiang Q, Khazaei S and Faghih RT (2025) A leptin-based Bayesian inference of a pro-satiety state reflects a basal circadian rhythm in women with obesity. Front. Endocrinol. 16:1638568. doi: 10.3389/fendo.2025.1638568
Received
30 May 2025
Accepted
25 August 2025
Published
16 September 2025
Volume
16 - 2025
Edited by
Rushendhiran Kesavan, University of Texas Southwestern Medical Center, United States
Reviewed by
Tatiana Romero-García, Universidad Autonoma de Baja California, Mexico
Satya Prakash Singh, Florida State University, United States
Updates
Copyright
© 2025 Xiang, Khazaei and Faghih.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rose T. Faghih, rfaghih@nyu.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.