Corrigendum: Sparse system identification of leptin dynamics in women with obesity

^{1}Department of Electrical and Computer Engineering, University of Houston, Houston, TX, United States^{2}Department of Nutrition and Health, Louis Bolk Instituut, Bunnik, Netherlands^{3}Department of Psychiatry, Weill Institute for Neurosciences, University of California, San Francisco, San Francisco, CA, United States

The prevalence of obesity is increasing around the world at an alarming rate. The interplay of the hormone leptin with the hypothalamus-pituitary-adrenal axis plays an important role in regulating energy balance, thereby contributing to obesity. This study presents a mathematical model, which describes hormonal behavior leading to an energy abnormal equilibrium that contributes to obesity. To this end, we analyze the behavior of two neuroendocrine hormones, leptin and cortisol, in a cohort of women with obesity, with simplified minimal state-space modeling. Using a system theoretic approach, coordinate descent method, and sparse recovery, we deconvolved the serum leptin-cortisol levels. Accordingly, we estimate the secretion patterns, timings, amplitudes, number of underlying pulses, infusion, and clearance rates of hormones in eighteen premenopausal women with obesity. Our results show that minimal state-space model was able to successfully capture the leptin and cortisol sparse dynamics with the multiple correlation coefficients greater than 0.83 and 0.87, respectively. Furthermore, the Granger causality test demonstrated a negative prospective predictive relationship between leptin and cortisol, 14 of 18 women. These results indicate that increases in cortisol are prospectively associated with reductions in leptin and vice versa, suggesting a bidirectional negative inhibitory relationship. As dysregulation of leptin may result in an abnormality in satiety and thereby associated to obesity, the investigation of leptin-cortisol sparse dynamics may offer a better diagnostic methodology to improve better treatments plans for individuals with obesity.

## 1 Introduction

Obesity, dubbed the “Global Epidemic” by the World Health Organization, is said to cause or aggravate various other health problems, worsening one’s life expectancy (1, 2). The prevalence of obesity is increasing worldwide at an alarming rate. The estimate shows that about 38% of the adult American population suffer from obesity (3). Obesity is associated with a reduced life expectancy (4). It also increases the risk of other chronic disorders such as diabetes and cardiovascular disease (5). Available studies suggest different lifestyle interventions for treating obesity such as a low-calorie diet, increased physical activity, and bariatric surgery (6). Nevertheless, many studies have pointed the determinants of failure of such approaches including adherence to lifestyle interventions (7), genetic background, adaptive changes in basal metabolic rate, hunger and satiety hormones that occur with weight loss (8). While there exists pathogenic approaches to treat the condition through studying its development, a salutogenic model focused on holistic wellbeing might be more effective as both a preventive and remedial measure (9, 10). For example, obesity can be conceptualized as a deficiency of energy regulation - i.e., the brain may fail to respond to negative feedback hormonal signals from adipose tissue, thereby perpetuating non-homeostatic eating behaviors that drive obesity. Thus, identifying the role of neuroendocrine hormones and adipokines in facilitating brain-adipose communication and energy regulation may inform novel strategies for the treatment and prevention of obesity. In this research, we propose models to understand leptin and cortisol behavior, which, in the future, can further be generalized for other relevant hormones, such as growth hormone, insulin, and ghrelin. Leptin is a signaling hormone and adipokine that is essential to activating central nervous system (CNS) networks involved in the suppression of appetite. It regulates food intake, metabolism, energy expenditure, and body weight (11). Although yet to be proven, Jacquier et al. (12) mathematically demonstrated how leptin resistance can be a result of persistent leptin infusion. They further pointed out that temporal alteration in some parameter values related to food intake can shift the equilibrium. The shift in equilibrium can lead to leptin resistance and obesity state. Blood leptin levels also correlate with changes in fat mass (13). Leptin is a hormone produced and secreted primarily by adipocytes. Hence, excessive amounts of adipose tissue in the body correspond to an increase in leptin production rate, and thus, higher serum leptin levels (11). Leptin interacts with the hypothalamus-pituitary-adrenal-axis (HPA-axis), thereby constituting an important biological mechanism whereby experiences of psychological stress may influence or be influenced by body weight and energy regulation (14, 15).

Another frequently studied hormone in obesity research is the glucocorticoid hormone cortisol (16). In response to different physiological demands such as blood glucose regulation, inflammation suppression, physiological stress, significant HPA-axis responses are induced across all ages and genders (17). Björntorp et al. (16) compared several studies and concluded that the association between obesity and cortisol is complicated. They concluded that there are several factors such regulatory systems, including the HPA, gonadal, growth-hormone, leptin axes, the sympathetic nervous system, the central adrenergic, serotoninergic, and dopaminergic systems that are related to obesity. Despite cortisol’s short half-life and circadian secretion pattern, most studies use cortisol measures such as serum, urine or salivary cortisol that do not reflect long-term cortisol exposure (18). Björntorp et al. (16) suggest that it is necessary to examine the cortisol secretion process to further understand the roots of obesity.

In an attempt to study obesity, various models based on ordinary differential equations considering regulations in energy and metabolism (19–23), as well as on the effects of different hormones such as ghrelin, cholecystokinin, and leptin have been hypothesized. Traditional studies perform statistical analysis on the measured serum hormone levels (24, 25). Some models include leptin, leptin resistance, and leptin receptors to understand obesity (12, 26, 27). Other studies show the association of glucocorticoids and weight gain (28, 29). Since obesity is associated with the polypeptide hormone, leptin (11, 12, 14), and the glucocorticoid hormone, cortisol (16, 30), studying them will prove vital in designing an approach to treating the condition. For example, we cannot measure leptin resistance in a living human, nevertheless, investigation of the association between cortisol and leptin through mathematical modeling may lead to some insights into the presence of leptin resistance.

Aschbacher et al. (23) propose a mathematical model based on the systems-level understanding of the HPA-leptin axis. Incorporating leptin in the model helped them to obtain three parameters: (I) inhibitory cortisol feedback signal, (II) Adrenocorticotropic hormone (ACTH)-adrenal interaction, and (III) leptin-cortisol proportionality. The authors suggested that the leptin-cortisol relationship may be important for understanding the neuroendocrine starvation response, and that poor system control may ultimately contribute to obesity. The mathematical model presented in the research by Aschbacher et al. (23) explains the rate of changes in cortisol concentration, ACTH, and leptin dynamics. The model used in this paper is an extension of the HPA system dynamics model used previously in (31) incorporating leptin’s impact on cortisol. Though leptin regulation is dependent on its relation with different hormones, understanding its behavior requires the study of leptin secretion and regulation. Investigating cortisol and leptin separately, and based on their physiology, will be a vital approach in interpreting how they behave concerning each other and other systems of the body. Moreover, cortisol and leptin each have circadian and ultradian rhythms, which our approach may better characterize. Jacquier et al. (12) introduce a mathematical model to study leptin resistance based on leptin synthesis and receptor activity. They study the changes in plasma leptin concentration, the density of leptin receptors, and food intake by varying leptin infusion and food consumption. This model studies leptin resistance dynamics based on many parameters such as food intake, fat mass, leptin receptors, and leptin. However, it is challenging to simultaneously collect this kind of clinical data. Therefore, we propose a simplified model considering only plasma leptin levels as the observation. As leptin secretion is a pulsatile process, we aim to exploit this characteristic by inferring underlying impulses.

In this research, we propose a simplified system to understand leptin regulation by exploiting the sparse nature of leptin secretion. If plasma leptin levels are sampled at a 10-minute sampling rate over 24-hours, we obtain a maximum of 40 leptin secretion events out of a total of 1440 secretion events, therefore, the nature being sparse. The sparse system identification of leptin dynamics over time makes it possible to recover the number of leptin pulses. It further allows us to recover information about the pulse amplitudes and intervals of occurrence from collected data. For an extensive understanding of leptin secretion, it is necessary to unveil the pulses originating from the hypothalamus as an outcome to leptin signaling. The pulses resulting in leptin regulation and secretion are the outcome of several contributing factors such as effects of other hormones on leptin regulation, density of leptin receptors, and food intake. Similar to leptin, cortisol secretion is a sparse process. Sparse system identification of cortisol will also lead us to recover the abstraction of secretion events coming from the HPA-axis. It also provides us important insight about the way it is infused by the adrenal glands and cleared by the liver. The proposed architecture has been previously verified for cortisol in healthy subjects and also on patients suffering from chronic fatigue and fibromyalgia syndromes (32–34). Here, we propose to study leptin-cortisol behavior by first observing the underlying pulses and estimating the physiological infusion and clearance rates in the system using two separate regulatory minimal models, one each for cortisol and leptin, and then determining the relationship between them. Similar to the results of previous study by Aschbacher et al. (23), we observe a negative relationship between cortisol and leptin. This relationship is observed based on a negative correlation. To further strengthen this observation, we perform the Granger causality test.

Figure 1A shows a pictorial representation of the leptin regulation and secretion model. Figure 1B shows a pictorial over all of the approach used in this research. This approach yields three important parameters: the timing and amplitudes of the leptin secretion events, the infusion rate of leptin by the adipose tissue, and the clearance rate by the renal system. We further conduct statistical analysis on the measured and estimated hormone levels. Finally, we propose that this method can be utilized to identify difference between patient and matched control participants for future characterization related to leptin deregulation similar to our previous approach (32).

**Figure 1** Leptin Regulation Model and an Overview of System-Theoretic Approach. **(A)** Shows the leptin secretion & regulation model. Leptin is the hormone to signal the hypothalamus regarding energy storage and use. It is produced and secreted primarily by adipocytes. The renal system is responsible for clearing it. **(B)** Shows the overall approach used in this study. A state-space model with physiological constraints is designed. A deconvolution approach is then used to estimate the secretion events and the physiological parameters. A statistical analysis is done on the estimations to understand the relationship of leptin-cortisol rhythms and dynamics in detail.

## 2 Results

We utilize hormone assay from eighteen premenopausal women with mean body mass index (BMI) 33 with range 30–41kgm^{-2}; mean age 37.5 with range 22–51 years; at morning (at 9 AM) mean cortisol level 11.02µg/dl and the corresponding standard deviation 3.65 *μg*/*dl*. Th distribution of average leptin level for all subject is 3.3583 ± 1.0633 *μg*/*dl*. The corresponding fat distribution is 39.61 ± 3.08 *kg* and lean body mass distribution is 55.094 ± 6.55 *kg*. Total cholesterol 4.7 ± 0.2 mmol/l (range 3.7–5.8 mmol/l), LDL cholesterol 2.99 ± 1.57 mmol/l (range 2.03– 4.00 mmol/l), and HDL cholesterol 1.54 ± 0.08 mmol/l (range 1.03–2.32 mmol/l). Homeostasis model assessment for insulin resistance (HOMA-IR) was estimated to be 3.25 ± 0.48 with placebo and 2.32 ± 0.19 after treating with Bromocriptine (35). The statistical test between two distribution showed significant difference with *p* = 0.01.

### 2.1 Pulsatile Nature of Leptin

Figure 2 shows the example deconvolution results of measured and reconstructed blood leptin levels of women with obesity for both experimental and simulated data for one subject. Figure 2A exhibits the estimated amplitudes and timings of hormonal secretory events. It also presents the experimental leptin data and estimates predicted by the model. The amplitude variations of the pulses are due to the circadian rhythm of underlying leptin release, and the variations in the timings are because of the ultradian rhythm (36). The number of recovered secretory events for all subjects were within the physiologically plausible range (i.e., 20-40 events per day). This experimental data includes 18 subjects with obesity. The black diamonds in Figure 2A represent the measured leptin levels obtained from blood samples. By performing deconvolution, we obtain the hormonal secretory pulses, which are used to obtain the reconstructed signal. The square of the multiple correlation coefficient (*R*^{2}) is between 0.8327 and 0.98603. *Y*_{1} and *Y*_{2} are the infusion rate of leptin by adipose tissue and the clearance rate of leptin by the renal system, respectively.

**Figure 2** Deconvolved Twenty-Four Hours Leptin Levels in Patient with Obesity. **(A)** Sub-plot shows the measured 24-hour leptin time series (black diamonds), the reconstructed leptin levels (cyan curve), the estimated pulse timings, and amplitudes (magenta vertical lines) for experimental data. **(B)** Sub-plot shows the simulated 24-hour leptin time series (black diamonds), the reconstructed leptin levels (cyan curve), the simulated pulse timings and amplitudes (cyan vertical lines), the recovered pulse timings, and amplitudes (magenta vertical lines) for simulated data.

Moreover, we simulated the data to further validate the proposed model. We simulated 18 leptin datasets, each corresponding to an experimental profile. An example of the result is presented in Figure 2B. The datasets are simulated from the estimated pulses of the experimental measurements shown in Figure 2A. We added Gaussian noise with standard deviation based on interassay co-efficient of variability (σ) provided in (23), defined in the Supplementary Information. Figure 2B shows an example of the ground truth of the sparse input, the estimated input, and the simulated leptin data. The blue stars show the estimated 24-hour leptin data, the black curve shows the estimated leptin levels, and the vertical blue lines show the amplitudes and timings of the simulated data. The vertical red lines in Figure 2B show the amplitudes and timings of the estimated hormone secretory events. Figure 3 shows the sample distribution for infusion and clearance rates for simulated and experimental leptin levels.

**Figure 3** Box-plot of Physiological Parameters of Leptin Levels. Subplots, respectively illustrate the sample distribution of **(A)** the infusion rate for simulated leptin data, **(B)** the clearance rate for simulated leptin data, **(C)** the infusion rate for experimental leptin data, and **(D)** the clearance rate for experimental leptin data, depicting the median (red line), the lower (Q1) to upper (Q3) quartile range (blue rectangle), and 9 to 91 percentile range (black line and black dashed line).

### 2.2 Pulsatile Nature of Cortisol

Figure 4 shows the measured and reconstructed blood cortisol levels of women with obesity. It shows the measured cortisol levels, reconstructed cortisol levels, and estimated amplitudes and timings of hormonal secretory events. The amplitude and timing variations of the pulses are due to the circadian rhythm and ultradian rhythm of the underlying cortisol release (37). The black diamonds in Figure 4 represent the measured cortisol levels obtained from blood samples. By performing deconvolution, we obtain the abstraction of hormonal secretory events (blue vertical lines in Figure 4), which are used to obtain the reconstructed signal (red curve in Figure 4). The square of the multiple correlation coefficient (*R*^{2}) is between 0.87557 and 0.98023. As this approach has been previously examined for cortisol data in different studies, we do not validate it on simulated data (33, 34, 38).

**Figure 4** Deconvolved Experimental Twenty-Four Hour Cortisol Levels in Patient with Obesity. Plot shows the measured 24-hour cortisol time series (black diamonds), the reconstructed cortisol levels (cyan curve), the estimated pulse timings and amplitudes (magenta vertical lines).

### 2.3 Leptin-Cortisol Antagonism

We seek to test how the leptin hormone samples and cortisol hormone samples are related. In this study, we look at the direct correlation between the sample regardless of other factors. Then, we also perform Granger causality analysis. This way we obtain an overall understanding of the apparent relationship between the leptin and cortisol hormone regardless of the other factors for the current dataset. Consideration of more flexible models that accounts for different factors such as age, gender, BMI, etc., and complexities related to nonlinearity is kept for future work.

#### 2.3.1 Spearman Correlation

We obtain Spearman correlation coefficient to measure the strength of a monotonic association between two variables based on the nonparametric measure of rank correlation. Table 1 shows the Spearman correlation coefficient between measured leptin levels and measured cortisol levels for 18 subjects. It also shows correlation coefficients for estimated as well as measured leptin and cortisol levels. The outcome for both these cases is very similar and shows that the proposed model retains the properties previously known.

Figure 5 shows the correlation plots for leptin and cortisol levels in subjects 5 and 12 who showed the highest negative correlation for both measured and estimated cases. We plot for both the measured and estimated levels to show the similarities. Except for in subjects 3, 15, and 18, we have observed significant correlation (i.e. *p* ≤ 0.05). All other correlation plots are provided in the Supplementary Information. This analysis shows, a negative association exists between leptin and cortisol levels across all time points in most, but not all, patients. The negative correlation on the same experimental dataset has been first identified by Aschbacher et al. (23). In this study, we see that the outcome also hold on the reconstructed data from the proposed model and sparse deconvolution.

**Figure 5** Correlation Plot between Leptin and Cortisol Levels in Patients with Obesity. The correlation plots are as follows: **(A)** measured leptin and cortisol levels for subject 5, **(B)** estimated leptin and cortisol levels for subject 5, **(C)** measured leptin and cortisol levels for subject 12, **(D)** estimated leptin and cortisol levels for subject 5. Each correlation plot incorporates a histogram and a scatter plot: the top left plot shows a histogram depicting leptin distributions and the bottom right plot shows a histogram depicting cortisol distributions, the bottom left plot and top right plot shows the scatter plot depicting correlation between leptin and cortisol.

#### 2.3.2 Granger Causality

We perform Granger causality analysis to infer the prospective prediction between leptin and cortisol. Figure 6 shows the results from the Granger causality analysis. Fourteen out of eighteen subjects showed statistically significant causal relationship from cortisol to leptin (i.e., cortisol samples from past can predict the present leptin samples). The corresponding lags are 36.64 ± 21.68 samples. Out of these fourteen subjects, ten showed the negative causal relationship. Seven out of eighteen showed statistically significant causal relationships from leptin to cortisol and four showed statistically significant negative relationships. The corresponding lags are 10.43 ± 17.95 samples. The evidence of possible negative association is in accordance with the previous findings in previous correlation based study by Aschbacher et al. (23).

**Figure 6** Granger Causality Map Estimated for Leptin and Cortisol for 18 Women with Obesity. Red color shows a negative causal influence and green color shows a positive causal influence, respectively. The direction of each arrow indicates that the variable at the start of the arrow is causing the variable at the end of the arrow. The thickness of the arrows represent the magnitude of the interaction. Starts next to the arrows denote that the corresponding prospective predictions are statistically significant with *p* < 0.05. * denotes statistical significance.

#### 2.3.3 Secretion Characteristics Based on Between Subject Variability

In order to investigate the between subject variability, we perform investigation of the leptin secretion characteristics such as the *l*_{2} -norm of the recovered leptin pulse vector (i.e. ||*u** _{l}*||

_{2}), infusions rates, and clearance rates against the factors such as BMI, lean body mass (LBM), HOMA-IR, and age with regression analysis. Similarly, we investigated the cortisol secretion characteristics such as the

*l*

_{2}-norm of the recovered leptin pulse vector (i.e. ||

*u**||*

_{c}_{2}), infusions rates, and clearance rates against the same factors factors. All the investigation were statistically insignificant (

*p*> 0) except for the two cases. The regression analysis results between the

*l*

_{2}-norm of leptin secretion (

*u*

_{l}) vs. LBM and the

*l*

_{2}-norm of leptin secretion (

*u*

_{l}) vs. BMI showed statistically significant relationships (

*p*< 0.05). Figure 7 show the scatter plot and the corresponding regression lines.

**Figure 7** Leptin Secretion Characteristics vs. BMI/LBM. Each panel shows: **(A)** the scatter plot (blue cross) of BMI vs l_{2}-norm of leptin secretion u_{2}, fitted regression line (red solid line), corresponding 95% confidence interval (dotted red curve), **(B)** the scatter plot (blue cross) of LBM vs l_{2}-norm of leptin secretion u_{1}, fitted regression line (red solid line), corresponding 95% confidence interval (dotted red curve).

## 3 Discussion

A complete model for the representation of hormonal variations must include all important intrinsic parameters such as forward and backward linkages between the hypothalamus, adipose tissue, anterior pituitary, adrenal gland, renal system, and liver as well as extrinsic parameters like stress, sleep, light, and food (37). It is a challenge to include all these factors while collecting human data, therefore, we propose this simplified minimal model considering the leptin measurements only. With the proposed minimal system theoratic model we formulate an optimization problem to infer unknown underlying the leptin secretion events as well as the corresponding system parameters. Furthermore, we performed the correlation analysis and the Granger causality analysis to investigate some possible directions to combine both cortisol and leptin in a single model in future. Investigation of secretion dynamics my offer a deeper understanding of the diseased conditions.

Understanding the secretion dynamics of leptin and cortisol in obese subjects and designing a model to understand their irregularities is difficult and challenging because the range of the number of leptin secretion events over a 24-hour period has not been accurately established. Licinio et al. (39) compared the leptin levels of women with obesity against healthy women and found that the concentration of independent pulsatile parameters of leptin, such as pulse duration and frequency, remained consistent between the groups, and that the excessive leptin levels in the persons with obesity were due to only an increased pulse height during secretion. The average number of leptin secretory events was found to be 32.0 ± 1.5 per day (39), within a range of 29 to 39 pulses among the subjects. In another study (40), they reported the average number of leptin pulses to be 30.0 ± 1 with the range of pulses between 21 to 39 over a 24-hour sampling period. These studies give a likely average of around 30 pulses per day of leptin, with a range between 20 to 40 being a possibility. We assumed there to be between 20-40 secretory events during a 24-hour period. Since there is no prior knowledge about the exact range of secretion events in women with obesity, we relax the constraints on the pulse range. We relaxed the upper and lower limits of this problem while taking care of the overfitting problem using Generalized Cross Validation - FOCal Under-determined System Solver (GCV-FOCUSS) to find *λ*. Although the upper limit was set to 50, we obtained less than or equal to 40 pulses for all patients. The proposed approach is applicable for a shorter length of hormone assay than a 24-hour assay, as long as it contains at least one secretion event. Therefore, the samples collected over an hour or more can be utilized to obtain the corresponding secretion event and the system parameters. As it can be seen in Figure 3, there is a difference in the distribution of the estimated parameters from the experimental data and the simulated data. This could be an indication of the consideration of a complex model in the optimization problem relative to the number of data-point available. Therefore proposed optimization problem might suffer from overfitting as pointed out in (41). As a future work of this, appropriate probabilistic priors can be applied on the system parameters to prevent such overfitting (42).

Unlike leptin, the mathematical modeling of glucocorticoids such as cortisol is a more widely studied problem. Using the coordinate descent approach, we recover the underlying pulses of cortisol along with the infusion of cortisol into plasma by the adrenal glands and clearance by the liver. After recovering the underlying leptin and cortisol pulses, we check the relation between them. In this research, to find the correlation between cortisol and leptin secretory patterns in women with obesity, we keep the units for both hormones consistent. Both cortisol and leptin measurements were converted to *μg*/*dL*. As compared to cortisol measurements, the leptin measurements have lower levels, making it more challenging to deconvolve. Aschbacher et al. observed a leptin-cortisol antagonism (23). This research further suggests that leptin-cortisol dynamics might provide an indirect or functional biomarker of the neuroendocrine starvation response (23). Based on our statistical analysis, we observe a negative correlation for 14 out of 18 subjects between measured leptin and cortisol levels, as well as a negative correlation between the estimated leptin and cortisol levels. However, as no lag has been considered during the correlation analysis which may explain the inhibition action, the negative correlation in this case may suggest an existence of an underlying factor which is causing both of these hormones to change in opposite direction. Therefore, the proposed simplified minimal model to study cortisol and leptin, and then finding the relationship between them, retains the previous properties. Consequently, exploiting the sparse and pulsatile nature of these hormones provides us with a more coherent insight.

As the distributions of the measured and the estimated samples are not Gaussian, we consider the rank-based Spearman correlation measure. A negative relationship has been observed between leptin and cortisol based on Spearman correlation for most of the subjects. The negative correlation between the reconstructed leptin and cortisol data also demonstrate that the proposed model and the corresponding sparse deconvolution also preserves the negative association. To further verify the interdependence of leptin and cortisol dynamics, we investigate the causal relationship between them. We perform the Granger causality analysis on leptin and cortisol to obtain the directional relationship. We observe a negative relationship, i.e., leptin-cortisol antagonism in 14 out of 18 subjects. leptin inhibits the HPA axis, which, *via* neuroendocrine hormones like ACTH, results in a reduction in cortisol (43, 44). Cortisol increases or stimulates production of leptin in adipose tissue (44). It is possible that this relationship may help explain behavioral contributors to obesity such as “stress-eating.” In the remaining 4 subjects, we do not observe an antagonism. It is possible that there exists a potential leptin resistance (45–48) in these subjects. Leptin resistance is a condition in which, even though the body has enough high leptin levels, receptors are desensitized leading to a “low signal.” Consequently, the brain may respond as if it were starving leading to a decrease in postprandial satiety. It is also possible that, due to a faulty leptin signaling, the brain is unable to realize the leptin levels in the body and produce high leptin levels during high cortisol levels. On the other hand, persistent higher concentration of leptin level plasma can also lead to leptin resistance (45, 47). Izquierdo et al. (47), also stated that obese people have high leptin levels and that the treatment of leptin resistance is still a great challenge in the treatment of obesity. Therefore, negative causation further provides a strong reason to study the negative relationship between leptin-cortisol. Furthermore, for some participants, the negative relationship is statistically not significant, while for other participants the relationship is positive. The result shows an indication of individual differences in leptin-cortisol regulation. Moreover, there are studies that did not find any relationship between cortisol and leptin in different groups such as children and adolescents investigating the relation in a single morning measurement of serum cortisol (49). Therefore, it would be very interesting to study these individual and group differences which will be crucial to develop more personalized obesity treatment programs. In other words, if there are differences in the hormone regulation dynamics in presence of obesity, these differences should be studied more thoroughly to before any intervention. This study has limitations in addition to observed variability including lack of measurement of insulin and ACTH, no control group without obesity, no data on the participants’ eating behaviors/physical activity/metabolism, use of Granger causality test, etc. All these factors should be also considered for understanding before considering any treatment program. Personalize treatment programs can be based on many parameters and preferences from patients, which can include possible existence of leptin resistance. That way personalized programs might benefit from taking the leptin cortisol dynamics and possible leptin resistance into account. Although the results obtained from Granger causality test provide some lights towards possible explanation, a conclusive understanding is still elusive which require further future investigation. The challenges includes limited number of subjects and observed between-subject variability. Further investigation is required with consideration of more related factors with larger population.

The leptin regulation model used in this research may be further be useful to understand the regulation dynamics in different symptoms such as leptin congenital deficiency, Alzheimer’s disease, and metabolic syndrome (50, 51). The advantage of using a model based on human physiology is that it can be used to further isolate the organs or tissues responsible for causing a particular deficiency. Determining a system model for the relationship between leptin and cortisol concentrations provides us with a more comprehensive understanding of the biological system’s behavior. This model can provide insight into the effects of the relationship of hormones on weight gain. There exist several bioscience and biosystems technologies used presently to simplify laboratory procedures and ease human life (52, 53). Similar to those available to obtain insulin levels in the human body (54), there will be bioscience technologies that would be able to obtain leptin levels in the human body in the near future.

Similar to leptin and cortisol, skin conductance responses are sparse. Wickramasuriya et al. presents a method to relate internal arousal state to the changes in skin conductance responses (55). Wickramasuriya et al. employ a deconvolution algorithm along with the Bayesian estimation to estimate the hidden arousal state. Azgomi et al. use control technique methods to close the loop and regulate the estimated arousal state within a desired range (56). Furthermore, Wickramasuriya et al. utilize the cortisol observations to track energy regulation in Cushing’s patients (57). Azgomi et al. employ the model presented in (57) to design a control system for closing the loop and regulating energy state in patients with hypercortisolism (58). Similarly, by deconvolving leptin levels and obtaining the hormonal pulses, one may develop a framework to obtain hunger states and regulate it. In another study, researchers in (57) have revealed the fact that cortisol secretion in healthy humans has a circadian nature (57). They derive multi-day data and observe differences between healthy humans and patients. A similar study can be done for leptin since leptin secretion is also influenced by the circadian rhythm. Similarly, another possible outcome of this study is to understand the leptin regulation behavior in women with obesity against their matched healthy subjects. Pednekar et al. (32) uses deconvolution to obtain the underlying hormonal pulses from cortisol data for fibromyalgia, chronic fatigue syndrome, and their matched healthy subjects. The healthy subjects are matched to the patients based on several factors such as age, sex, and weight. The comparison between the parameters and pulses obtained from the patients against their matched healthy subjects showed significant differences in clearance rates and hormonal pulses. A similar study can be done in the future, using leptin data in women with obesity. As another future direction, the model could be further extended to include ACTH and study the differences in hormonal secretory events in women with obesity with respect to healthy subjects (38). These hormonal secretory events are an abstraction of their pulsatile nature and provide a good estimation of the pulses originating in the CNS. A thorough investigation of these events could further help us to understand the role of the CNS as well as kidneys and liver in causing obesity.

In this study, we have proposed a state-space model for describing sparse leptin hormone dynamics. Furthermore, we have proposed a coordinate descent approach for the sparse deconvolution of the sampled leptin data to obtain the pulsatile secretion events. The deconvolution framework also provides estimates of the state-space model parameters. The successful modeling of leptin is the first step to a more comprehensive understanding of dynamics and its interaction with cortisol. Our analysis with Spearman correlation and Granger causality has provided some evidence of negative relation between leptin and cortisol. However, some variability has also been observed which prevents us from concluding a concrete explanation. A more thorough investigation is required as the study has limitation such as lack of including insulin and ACTH measurement, no control group without obesity, no data on the participants’ eating behaviors/physical activity/metabolism, use of Granger causality test that only consider two factors, a small number of participants, low number of samples, etc.

## 4 Methods

### 4.1 Experimental Dataset

In this research, we utilize blood sample data collected from eighteen premenopausal women with obesity for our analysis of leptin and cortisol pulse recovery (59). The participants have given informed consent for participation in the study. They collected data for 18 premenopausal women with obesity over a 24-hour period with 10-minute sampling frequency. All participants were required to have regular menstrual cycles. All studies were done in the early follicular phase of the menstrual cycle. The subjects were within the age group of 22-51 and with BMI between 30-41 kg/m^{2}. Subjects were admitted to the research center at 7:00 after an overnight fast. After the subjects rested for 45 min, indirect calorimetry was performed using a ventilated hood for 30 min. Thereafter, a cannula for blood sampling was inserted into an antecubital vein, which was attached to a three-way stopcock and kept patent by a continuous 0.9% NaCl and heparin (1 U/ml) infusion (500 ml/24 h). Each subject followed a consistent schedule with breakfast, lunch, and dinner provided, and had an 8.5 hour sleeping period between 23:00 and 7:30. The blood samples were collected without waking the subjects during this period. Mean fasting glucose concentration was 5.0 ± 0.1 mmol/l (range 4.2– 6.3 mmol/l), insulin 15.3 ± 1.7 mU/l (range 7–28 mU/l), glycosylated hemoglobin (HbA1C) 4.7 ± 0.1% (range 3.9 –5.3%). The above study was approved by the Medical Ethics Committee of Leiden University (23). Subjects were required to be free of chronic disease, and exclusion criteria were fixed, such as shift work, depression, alcohol abuse, and oral contraceptives. The blood samples were assayed every 10 minutes for cortisol and every 20 minutes for leptin and intermediate points were interpolated. Plasma leptin concentrations were estimated by radioimmunoassay with a detection limit of 0.5 *ng*/*liter*, and interassay coefficient of variation of 3.6% to 6.8%. Plasma cortisol was measured with a radioimmunoassay with a detection limit is 25 *nmol*^{–1} and the intra-assay coefficient of variation ranges between 2% and 4%. A detailed description of the experiment is provided in (23, 59). In this research, we analyze the leptin and cortisol abstraction pulses for 18 subjects obtained from this dataset.

### 4.2 Leptin Model Formulation

The state-space model introduced in this research considers the first-order differential system of equations for leptin synthesis in the adipose tissue and clearance by the renal system. In a 24-hour period, the blood samples were collected at a 10-minute interval, i.e., 144samples in a day. However, leptin was assayed every 20 min and concentrations was obtained for every 10 min by interpolation. If we assume that leptin secretion can happen in every minutes then there are 1440 possible events for 24-h duration. In a discrete space of 1440 events, we observe between 20 to 40 pulses making it sparse in nature. In this model, we intend to utilize the sparse nature of the hormonal secretory events along with other physiological constraints in a state-space model to estimate the amplitude and frequency of hormonal secretory events. The rate of change of leptin concentration in the adipose tissue is equal to the difference between leptin synthesis rate and the leptin infusion rate from adipose tissue into the blood. Similarly, the rate of change of leptin concentration in the serum is equal to the difference between the leptin infusion rate from adipose tissue into the blood and the leptin clearance rate by the renal system (60). The leptin secretion dynamics are represented as follows:

where *x*_{1}(*t*) and *x*_{2}(*t*) represent the effective leptin concentrations inside adipose tissue and leptin concentration in the blood serum, respectively. The model parameters *Y*_{1} and *Y*_{2} represent the infusion rate of leptin by the adipose tissue and the clearance rate of leptin by the kidneys, respectively. Input *u _{l}*(

*t*) represents the abstraction of effective pulses mainly produced by the adipose tissue.

*u*(

_{l}*t*) can be modeled as a summation of delta functions with ${u}_{l}(t)={\Sigma}_{i=1}^{N}\text{}{q}_{i}\text{}\delta (t-{\tau}_{i})$ where

*q*is the magnitude of the pulse initiated at time

_{i}*τ*. If no pulse occurs at

_{i}*τ*,

_{i}*q*will equal zero. The pulses are assumed to occur at integer minutes, i.e., in a 24-hour period, there are 1440 distinct locations (

_{i}*N*= 1440). The blood samples are collected every 10 minutes for

*M*samples (

*M*= 144). Blood samples were collected beginning at ${y}_{{l}_{o}}$ and then assayed for leptin. Let ${y}_{{t}_{10}},{y}_{{t}_{20}},\cdots ,{y}_{{t}_{10M}},{t}_{k}:k=10,20,\dots ,10M$. Furthermore, ${y}_{{l}_{o}}={y}_{{t}_{10}}$.

where ${y}_{{t}_{k}}$ and ${v}_{{t}_{k}}$ represent the observed leptin level and hormone measurement error at a time *t _{k}*, respectively. We model ${v}_{{t}_{k}}$ random variables as Gaussian random variable. Then for the system output, we assume that the observed output value at time

*t*will be equal to the previous value, y

_{k}_{0}, multiplied by a decay term and added to a secretion value and the error term. We can represent the solution for every time point

*t*with the following equation:

_{k}where ${\mathbf{\text{b}}}_{{\mathbf{\text{t}}}_{\mathbf{\text{k}}}}={\left[\frac{{\gamma}_{1}}{{\gamma}_{1}+{\gamma}_{2}}({e}^{-{\gamma}_{2}k}-{e}^{-{\gamma}_{1}k})\frac{{\gamma}_{1}}{{\gamma}_{1}-{\gamma}_{2}}({e}^{-{\gamma}_{2}(k-1)}-{e}^{-{\gamma}_{1}(k-1)})\cdots \frac{{\gamma}_{1}}{{\gamma}_{1}-{\gamma}_{2}}({e}^{-{\gamma}_{2}}-{e}^{-{\gamma}_{1}})\underset{N-k}{\underset{\backslash underbrace\backslash}{\begin{array}{ccc}0& \cdots & 0\end{array}}}\right]}^{\text{'}}$ and ${a}_{{t}_{k}}={e}^{-{\gamma}_{2}k}$. *u** _{l}* represents the input over the entire 24-hour sampling period, with values

*q*over

_{i}*i*= 1,…,1440. ${v}_{{t}_{k}}$ stands for the vector representation of the random variable ${v}_{{t}_{k}}$. We then solve for ${\text{a}}_{{t}_{k}}$ and ${b}_{{t}_{k}}$ using a forced solution approach, multiplying each side of the equation by

*e*

^{–}

*and using mathematical methods to obtain the solutions.*

^{at}From these matrices, we can form a combined representation for the system at any time, using

We then represent the system output as:

Where output vector *y _{l}* is dependent on the initial signal value, the input pulses, and the error term. The values of matrices

*A*

_{γ}and

*B*

_{γ}are dependent on the values of

*y*for the given subject.

### 4.3 Cortisol Model Formulation

In this research, we use the cortisol secretion model provided by Faghih et al. (33). It exploits the sparse nature of hormonal secretory events and other physiological constraints to estimate the amplitude and timings of the secretory events using a state-space model. The rate at which cortisol concentration changes in the adrenal glands is equal to the difference between the rate at which it is infused in the blood and the rate at which it is secreted. The rate of change of serum cortisol concentration is equal to the difference between the rate at which it is infused by the adrenal gland into the blood and the rate at which it is cleared by the liver from the blood. The cortisol secretion dynamics are represented as follows:

where *x*_{3}(*t*) and *x*_{4}(*t*) represent the cortisol concentrations in the adrenal glands and the blood serum, respectively. The physiological model parameters *ψ*_{1} and *ψ*_{2} represent the infusion rate of cortisol by the adrenal glands and the cortisol clearance rate by the liver respectively. Input *u _{c}*(

*t*) represent the abstraction of pulses coming from the hypothalamus responsible for the production of cortisol.

*u*(

_{c}*t*) can be modeled as a summation of delta functions with ${u}_{c}(t)={\Sigma}_{i=1}^{N}\text{}{p}_{i}\text{}\delta (t-{\tau}_{i})$.

*p*is the magnitude of the pulse initiated at time

_{i}*τ*. If no pulse occurs at

_{i}*τ*,

_{i}*q*will equal zero. The pulses are assumed to occur at integer minutes, i.e. in a 24-hour period, there are 1440 distinct locations (

_{i}*N*= 1440). The blood samples are collected every 10 minutes, for

*M*samples (M = 144). As explained in Section 4.2 the measurement output can be represented as:

where ${y}_{{t}_{k}}$ and ${v}_{{t}_{k}}$ are the observed cortisol concentration and the measurement error, respectively. The initial concentration of cortisol in adrenal glands and serum is assumed to be zero and ${y}_{{c}_{o}}$. The system is further expressed as:

where

#### 4.3.1 Granger Causality Analysis

We utilize Granger Causality to determine the prospective prediction between leptin and cortisol. Let ${y}_{{c}_{{t}_{k}}}$ and ${y}_{{l}_{{t}_{k}}}$ are two stationary processes sampled every 10 minutes. A simple causal model for the interaction between leptin and cortisol is given as follows (61):

According to the definition of the causality, ${y}_{{c}_{{t}_{k}}}$ is causing ${y}_{{l}_{{t}_{k}}}$ provided some *β _{l,j}* is not zero and ${y}_{{l}_{{t}_{k}}}$ is causing ${y}_{{c}_{{t}_{k}}}$ if some

*β*is not zero. If both

_{c,j}*β*and

_{c,j}*β*are not zero, there is exist a feedback relationship between ${y}_{{c}_{{t}_{k}}}$ and ${y}_{{l}_{{t}_{k}}}$. Here,

_{l,j}*ha*and

_{j,j}*α*represents the coefficients related to the self interaction. ${\in}_{{l}_{t}}$ and ${\in}_{{c}_{t}}$ represent the model errors. We model them as i.i.d Gaussian random variables.

_{c,j}We use the Bayesian information criteria (BIC) to find ${m}_{{l}_{1}}$ and ${m}_{{c}_{1}}$ considering all *β _{l,j}* =

*β*= 0. In contrast, we find the ${m}_{{l}_{2}}$ and ${m}_{{c}_{2}}$ using BIC considering ${m}_{{l}_{1}}$ and ${m}_{{c}_{1}}$ numbers of nonzero

_{c,j}*α*and

_{l,j}*α*in (7) and (8), respectively. We first fit the models in (7) and (8) to the data with

_{c,j}*β*=

_{l,j}*β*= 0. We obtain the likelihoods as ${L}_{l}({y}_{{l}_{{t}_{k}}})$ and ${L}_{c}({y}_{{c}_{{t}_{k}}})$. Then, we assume

_{c,j}*β*,

_{l,j}*β*≠ 0 and fit the model to the data. We obtain the corresponding likelihood ${L}_{l}({y}_{{l}_{{t}_{k}}}|{y}_{{c}_{{t}_{k}}})$ and $I{L}_{c}({y}_{{c}_{{t}_{k}}}|{y}_{{l}_{{t}_{k}}})$. Therefore, the Granger causality depicts how cortisol is causing leptin,

_{c,j}and how leptin is causing cortisol,

where negative sign denotes inhibitory, positive sign denotes excitatory, and zero means no prospective prediction. The magnitude gives the relative strength of the prospective prediction. Furthermore, we calculate the F-statistics for testing whether the prospective prediction is statistically significant or not (61).

## 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 human participants were reviewed and approved by Medical Ethics Committee of Leiden University. The patients/participants provided their written informed consent to participate in this study.

## Author Contributions

RF conceived and designed the study. MA and DP performed the research, analyzed data and wrote the manuscript. RF and MA developed the algorithms and analysis tools. All authors revised the manuscript. All authors contributed to the article and approved the submitted version.

## Funding

This work was supported in part by NSF grant 1942585 - CAREER: MINDWATCH: Multimodal Intelligent Noninvasive brain state Decoder for Wearable AdapTive Closed-loop arcHitectures and 1755780 - CRII: CPS: Wearable-Machine Interface Architectures. RF as the senior author. NYU faculty start-up funds covered the article processing charges.

## 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.

## 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.2022.769951/full#supplementary-material

## References

1. Ng M, Fleming T, Robinson M, Thomson B, Graetz N, Margono C, et al. Global, Regional, and National Prevalence of Overweight and Obesity in Children and Adults During 1980–2013: A Systematic Analysis for the Global Burden of Disease Study 2013. *Lancet* (2014) 384:766–81. doi: 10.1016/S0140-6736(14)60460-8

2. Wang Y, Beydoun MA, Min J, Xue H, Kaminsky LA, Cheskin LJ. Has the Prevalence of Overweight, Obesity and Central Obesity Levelled Off in the United States? Trends, Patterns, Disparities, and Future Projections for the Obesity Epidemic. *Int J Epidemiol* (2020) 49(3):810–23. doi: 10.1093/ije/dyz273

3. Liou D, Kulik L. Self-Efficacy and Psychosocial Considerations of Obesity Risk Reduction Behaviors in Young Adult White Americans. *PLoS One* (2020) 15:e0235219. doi: 10.1371/journal.pone.0235219

4. Blüher M. Metabolically Healthy Obesity. *Endocrine Rev* (2020) 41:405–20. doi: 10.1038/s41366-020-00667-4

5. Jastreboff AM, Kotz CM, Kahan S, Kelly AS, Heymsfield SB. Obesity as a Disease: The Obesity Society 2018 Position Statement. *Obesity* (2019) 27:7–9. doi: 10.1002/oby.22378

6. Cardel MI, Atkinson MA, Taveras EM, Holm JC, Kelly AS. Obesity Treatment Among Adolescents: A Review of Current Evidence and Future Directions. *JAMA Pediatr* (2020) 174:609–17. doi: 10.1001/jamapediatrics.2020.0085

7. Burgess E, Hassmén P, Pumpa KL. Determinants of Adherence to Lifestyle Intervention in Adults With Obesity: A Systematic Review. *Clin Obes* (2017) 7:123–35. doi: 10.1111/cob.12183

8. Reinehr T. Lifestyle Intervention in Childhood Obesity: Changes and Challenges. *Nat Rev Endocrinol* (2013) 9:607–14. doi: 10.1038/nrendo.2013.149

9. WHO, World Health Organization. *Obesity: Preventing and Managing the Global Epidemic*, Vol. 894. Geneva, Switzerland: Geneva WHO (2000).

10. Seo YG, Lim H, Kim Y, Ju YS, Lee HJ, Jang HB, et al. The Effect of a Multidisciplinary Lifestyle Intervention on Obesity Status, Body Composition, Physical Fitness, and Cardiometabolic Risk Markers in Children and Adolescents With Obesity. *Nutrients* (2019) 11:137. doi: 10.3390/nu11010137

11. Mantzoros CS, Magkos F, Brinkoetter M, Sienkiewicz E, Dardeno TA, Kim SY, et al. Leptin in Human Physiology and Pathophysiology. *Am J Physiology-Endocrinology Metab* (2011) 301:E567–84. doi: 10.1152/ajpendo.00315.2011

12. Jacquier M, Soula HA, Crauste F. A Mathematical Model of Leptin Resistance. *Math Biosci* (2015) 267:10–23. doi: 10.1016/j.mbs.2015.06.008

13. Ahima RS, Flier JS. Leptin. *Annu Rev Physiol* (2000) 62:413–37. doi: 10.1146/annurev.physiol.62.1.413

14. Buckley TM, Schatzberg AF. On the Interactions of the Hypothalamic-Pituitary-Adrenal (Hpa) Axis and Sleep: Normal Hpa Axis Activity and Circadian Rhythm, Exemplary Sleep Disorders. *J Clin Endocrinol Metab* (2005) 90:3106–14. doi: 10.1210/jc.2004-1056

15. Adam TC, Epel ES. Stress, Eating and the Reward System. *Physiol Behav* (2007) 91:449–58. doi: 10.1016/j.physbeh.2007.04.011

16. Björntorp P, Rosmond R. Obesity and Cortisol. *Nutrition* (2000) 16:924–36. doi: 10.1016/S0899-9007(00)00422-6

17. Kudielka B, Buske-Kirschbaum A, Hellhammer D, Kirschbaum C. Hpa Axis Responses to Laboratory Psychosocial Stress in Healthy Elderly Adults, Younger Adults, and Children: Impact of Age and Gender. *Psychoneuroendocrinology* (2004) 29:83–98. doi: 10.1016/S0306-4530(02)00146-4

18. Rodriguez ACI, Epel ES, White ML, Standen EC, Seckl JR, Tomiyama AJ. Hypothalamic-Pituitary-Adrenal Axis Dysregulation and Cortisol Activity in Obesity: A Systematic Review. *Psychoneuroendocrinology* (2015) 62:301–18. doi: 10.1016/j.psyneuen.2015.08.014

19. de Graaf AA, Freidig AP, De Roos B, Jamshidi N, Heinemann M, Rullmann JA, et al. Nutritional Systems Biology Modeling: From Molecular Mechanisms to Physiology. *PLoS Comput Biol* (2009) 5:e1000554. doi: 10.1371/journal.pcbi.1000554

20. Chow CC, Hall KD. Short and Long-Term Energy Intake Patterns and Their Implications for Human Body Weight Regulation. *Physiol Behav* (2014) 134:60–5. doi: 10.1016/j.physbeh.2014.02.044

21. Horgan GW. The Behaviour of a Neutral Model of Weight Regulated Only by Body Mass. *J Theor Biol* (2011) 270:1–6. doi: 10.1016/j.jtbi.2010.11.001

22. Hall KD. Predicting Metabolic Adaptation, Body Weight Change, and Energy Intake in Humans. *Am J Physiology-Endocrinology Metab* (2009) 298:E449–66. doi: 10.1152/ajpendo.00559.2009

23. Aschbacher K, Rodriguez-Fernandez M, van Wietmarschen H, Tomiyama AJ, Jain S, Epel E, et al. The Hypothalamic–Pituitary–Adrenal–Leptin Axis and Metabolic Health: A Systems Approach to Resilience, Robustness and Control. *Interface Focus* (2014) 4:20140020. doi: 10.1098/rsfs.2014.0020

24. Kanaley JA, Weltman JY, Pieper KS, Weltman A, Hartman ML. Cortisol and Growth Hormone Responses to Exercise at Different Times of Day. *J Clin Endocrinol Metab* (2001) 86:2881–9. doi: 10.1210/jc.86.6.2881

25. Antunes H, Santos C, Carvalho S. Serum Leptin Levels in Overweight Children and Adolescents. *Br J Nutr* (2008) 101:1262–6. doi: 10.1017/S0007114508055682

26. Tam J, Fukumura D, Jain RK. A Mathematical Model of Murine Metabolic Regulation by Leptin: Energy Balance and Defense of a Stable Body Weight. *Cell Metab* (2009) 9:52–63. doi: 10.1016/j.cmet.2008.11.005

27. De Gaetano A, Hardy T, Beck B, Abu-Raddad E, Palumbo P, Bue-Valleskey J, et al. Mathematical Models of Diabetes Progression. *Am J Physiology-Endocrinology Metab* (2008) 295:E1462–79. doi: 10.1152/ajpendo.90444.2008

28. Udden J, Björntorp P, Arner P, Barkeling B, Meurling L, Rössner S. Effects of Glucocorticoids on Leptin Levels and Eating Behaviour in Women. *J Internal Med* (2003) 253:225–31. doi: 10.1046/j.1365-2796.2003.01099.x

29. Sominsky L, Spencer SJ. Eating Behavior and Stress: A Pathway to Obesity. *Front Psychol* (2014) 5:434. doi: 10.3389/fpsyg.2014.00434

30. Myers MG Jr., Leibel RL, Seeley RJ, Schwartz MW. Obesity and Leptin Resistance: Distinguishing Cause From Effect. *Trends Endocrinol Metab* (2010) 21:643–51. doi: 10.1016/j.tem.2010.08.002

31. Aschbacher K, Adam EK, Crofford LJ, Kemeny ME, Demitrack MA, Ben-Zvi A. Linking Disease Symptoms and Subtypes With Personalized Systems-Based Phenotypes: A Proof of Concept Study. *Brain Behav Immun* (2012) 26:1047–56. doi: 10.1016/j.bbi.2012.06.002

32. Pednekar DD, Amin MR, Azgomi HF, Aschbacher K, Crofford LJ, Faghih RT. Characterization of Cortisol Dysregulation in Fibromyalgia and Chronic Fatigue Syndromes: A State-Space Approach. *IEEE Trans Biomed Eng* (2020) 67:3163–72. doi: 10.1109/TBME.2020.2978801

33. Faghih RT, Dahleh MA, Adler GK, Klerman EB, Brown EN. Deconvolution of Serum Cortisol Levels by Using Compressed Sensing. *PLoS One* (2014) 9:e85204. doi: 10.1371/journal.pone.0085204

34. Pednekar DD, Amin MR, Fekri H, Aschbacher CLJK, Faghih RT. A System Theoretic Investigation of Cortisol Dysregulation in Fibromyalgia Patients With Chronic Fatigue. In: 2019 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC). Berlin: IEEE (2019).

35. Kok P, Roelfsema F, Frolich M, van Pelt J, Stokkel MP, Meinders AE, et 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

36. Simon C, Gronfier C, Schlienger J, Brandenberger G. 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

37. Brown EN, Meehan PM, Dempster AP. A Stochastic Differential Equation Model of Diurnal Cortisol Patterns. *Am J Physiology-Endocrinology Metab* (2001) 280:E450–61. doi: 10.1152/ajpendo.2001.280.3.E450

38. Faghih RT, Dahleh MA, Adler GK, Klerman EB, Brown EN. Quantifying Pituitary-Adrenal Dynamics and Deconvolution of Concurrent Cortisol and Adrenocorticotropic Hormone Data by Compressed Sensing. *IEEE Trans Biomed Eng* (2015) 62:2379–88. doi: 10.1109/TBME.2015.2427745

39. Licinio J, Mantzoros C, Negrão AB, Cizza G, Wong ML, Bongiorno PB, et al. Human Leptin Levels Are Pulsatile and Inversely Related to Pituitary–Ardenal Function. *Nat Med* (1997) 3:575. doi: 10.1038/nm0597-575

40. Licinio J, Negrão AB, Mantzoros C, Kaklamani V, Wong ML, Bongiorno PB, et al. Sex Differences in Circulating Human Leptin Pulse Amplitude: Clinical Implications. *J Clin Endocrinol Metab* (1998) 83:4140–7. doi: 10.1210/jc.83.11.4140

41. Ghasemi Z, Lee JC, Kim CS, Cheng HM, Sung SH, Chen CH, et al. Estimation of Cardiovascular Risk Predictors From Non-Invasively Measured Diametric Pulse Volume Waveforms *via* Multiple Measurement Information Fusion. *Sci Rep* (2018) 8:1–11. doi: 10.1038/s41598-018-28604-6

42. Amin MR, Faghih RT. Identification of Sympathetic Nervous System Activation From Skin Conductance: A Sparse Decomposition Approach With Physiological Priors. *IEEE Trans Biomed Eng* (2020) 68:1726–36. doi: 10.1109/TBME.2020.3034632

43. Bornstein SR, Uhlmann K, Haidan A, Ehrhart-Bornstein M, Scherbaum WA. Evidence for a Novel Peripheral Action of Leptin as a Metabolic Signal to the Adrenal Gland: Leptin Inhibits Cortisol Release Directly. *Diabetes* (1997) 46:1235–8. doi: 10.2337/diab.46.7.1235

44. Pralong FP, Roduit R, Waeber G, Castillo E, Mosimann F, Thorens B, et al. Leptin Inhibits Directly Glucocorticoid Secretion by Normal Human and Rat Adrenal Gland. *Endocrinology* (1998) 139:4264–8. doi: 10.1210/endo.139.10.6254

45. Enriori PJ, Evans AE, Sinnayah P, Cowley MA. Leptin Resistance and Obesity. *Obesity* (2006) 14:254S–8S. doi: 10.1038/oby.2006.319

46. Engin A. Diet-Induced Obesity and the Mechanism of Leptin Resistance. *Obes Lipotoxicity* (2017) 381–97. doi: 10.1007/978-3-319-48382-5_16

47. Izquierdo AG, Crujeiras AB, Casanueva FF, Carreira MC. Leptin, Obesity, and Leptin Resistance: Where Are We 25 Years Later? *Nutrients* (2019) 11:2704. doi: 10.3390/nu11112704

48. Gruzdeva O, Borodkina D, Uchasova E, Dyleva Y, Barbarash O. Leptin Resistance: Underlying Mechanisms and Diagnosis. *Diabetes Metab Syndrome Obesity: Targets Ther* (2019) 12:191. doi: 10.2147/DMSO.S182406

49. Sudi K, Gallistl S, Weinhandl G, Payer C, Cartellieri M, Borkenstein MH. No Relationship Between Leptin and Cortisol in Obese Children and Adolescents. *J Pediatr Endocrinol Metab* (2000) 13:913–22. doi: 10.1515/JPEM.2000.13.7.913

50. Yupanqui-Lozno H, Bastarrachea RA, Yupanqui-Velazco ME, Alvarez-Jaramillo M, Medina-Méndez E, Giraldo-Peña AP, et al. Congenital Leptin Deficiency and Leptin Gene Missense Mutation Found in Two Colombian Sisters With Severe Obesity. *Genes* (2019) 10:342. doi: 10.3390/genes10050342

51. Maioli S, Lodeiro M, Merino-Serrais P, Falahati F, Khan W, Puerta E, et al. Alterations in Brain Leptin Signalling in Spite of Unchanged Csf Leptin Levels in Alzheimer’s Disease. *Aging Cell* (2015) 14:122–9. doi: 10.1111/acel.12281

52. Salehi SM. Engineering Cell-Free Protein Synthesis Technology for Codon Reassignment, Biotherapeutics Production Using Just-Add-Water System, and Biosensing Endocrine Disrupting Compounds. (2017). doi: 10.1021/acs.analchem.6b04034

53. Ghosh M, Sharma N, Singh AK, Gera M, Pulicherla KK, Jeong DK. Transformation of Animal Genomics by Next-Generation Sequencing Technologies: A Decade of Challenges and Their Impact on Genetic Architecture. *Crit Rev Biotechnol* (2018) 38:1157–75. doi: 10.1080/07388551.2018.1451819

54. Dias D, Paulo Silva Cunha J. Wearable Health Devices—Vital Sign Monitoring, Systems and Technologies. *Sensors* (2018) 18:2414. doi: 10.3390/s18082414

55. Wickramasuriya DS, Amin M, Faghih RT. Skin Conductance as a Viable Alternative for Closing the Deep Brain Stimulation Loop in Neuropsychiatric Disorders. *Front Neurosci* (2019) 13:780. doi: 10.3389/fnins.2019.00780

56. Azgomi HF, Wickramasuriya DS, Faghih RT. State-Space Modeling and Fuzzy Feedback Control of Cognitive Stress, In: 2019 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC). Berlin, Germany: IEEE (2019). pp. 6327–30.

57. Wickramasuriya DS, Faghih RT. A Cortisol-Based Energy Decoder for Investigation of Fatigue in Hypercortisolism, In: 2019 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC). Berlin, Germany: IEEE (2019). pp. 11–4.

58. Azgomi HF, Faghih RT. A Wearable Brain Machine Interface Architecture for Regulation of Energy in Hypercortisolism, In: 2019 53rd Asilomar Conference on Signals, Systems, and Computers. Pacific Grove, CA, USA: IEEE (2019). pp. 254–8.

59. Kok P, Roelfsema F, Frolich M, van Pelt J, Meinders AE, Pijl H. 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

60. Wolf G, Chen S, Han DC, Ziyadeh FN. Leptin and Renal Disease. *Am J Kidney Dis* (2002) 39:1–11. doi: 10.1053/ajkd.2002.29865

Keywords: obesity, leptin, cortisol, state-space, sparse recovery, deconvolution, correlation, grange-causality

Citation: Amin MR, Pednekar DD, Azgomi HF, van Wietmarschen H, Aschbacher K and Faghih RT (2022) Sparse System Identification of Leptin Dynamics in Women With Obesity. *Front. Endocrinol.* 13:769951. doi: 10.3389/fendo.2022.769951

Received: 02 September 2021; Accepted: 24 January 2022;

Published: 05 April 2022.

Edited by:

Hanneke M. Van Santen, Wilhelmina Children’s Hospital, NetherlandsReviewed by:

Ana Paula Santos-Silva, Federal University of Rio de Janeiro, BrazilOzair Abawi, Erasmus Medical Center, Netherlands

Copyright © 2022 Amin, Pednekar, Azgomi, van Wietmarschen, Aschbacher 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

^{†}Present address: Biomedical Engineering, New York University, New York City, NY, United States