1. Introduction
Gout is the most common form of inflammatory arthritis and is caused by a chronic elevation of the serum uric acid concentration above the saturation level for monosodium urate crystal formation [1]. The deposition of monosodium urate crystals, which occurs predominantly in peripheral joints and long-term deposition of monosodium urate crystals, can result in joint damage. Consequently, gout has a substantial effect on the physical functioning, productivity and quality of life of those affected, not to mention high health care costs. Gout is also associated with metabolic syndrome, cardiovascular diseases and renal diseases [2] [3] [4]. Furthermore, Gout has no cure and it can lead to a higher risk of death [5].
The proportion of individuals, within a population, with gout is highly variable across various regions of the world. However, gout is common in most countries in North America, Western Europe and Asia [1] [6]. More than eight million people in the U.S., almost 4% of the population, are estimated to suffer from gout [7]. Although no formal survey has been undertaken in Canada, gout is generally thought to affect 3% of adults [8]. Greece has the highest reported prevalence of gout in Europe, at 4.75% of the adult population [9]. Using 2009 data from the Aotearoa New Zealand Health Tracker, in which the case definition of gout was consistent across ethnic groups, Pacific islanders and Maori had a threefold greater risk of gout than those of European descent after adjustment for age and sex. The crude prevalence estimates in adults being 7.63%, 6.06% and 3.24% in these three groups respectively [10].
Recently, mathematical models have been proven to be an effective research methodology and have become a widely accepted tool in biology and the medical sciences. For instance, Tolic et al. [11], analyzed a mathematical model of the insulin-glucose feedback regulation in men and the results suggest that interactions between the oscillatory insulin supply and the receptor dynamics can be of minute significance only. Li and Kuang [12] analyzed a mathematical model of glucose-insulin with time delay that provides qualitatively robust dynamics for clinical applications. Song and Thomas [13] developed a differential equation model describing the dynamics of stored energy in the form of fat mass, lean body mass, and ketone body mass during prolonged starvation. Pearson et al. [14] derived a system of differential equations that describes the transport between and storage in different tissues of the human body and their results are in agreement with experimental data. Similarly, Song et al. [15] analyzed a model of insulin delivery with an open-loop control and time delays in an insulin pump system, their results suggest that a smaller dose with higher delivery frequency has a better effect on continuous subcutaneous insulin injection administration. Furthermore, Jacquier et al. [16] proposed a mathematical model of the leptin-leptin receptor system based on the assumption that leptin is a regulator of its own receptor activity. Similarly, while the interaction between insulin, leptin, uric acid and the metabolism is complex, the dynamics could be accurately captured using a mathematical model similar to the ones previously mentioned.
The rest of this paper is organized as follows. Methodology and underlying mechanisms for our model are presented in Section 2. In Section 3, we present a mathematical model to study the dynamics of uric acid. In Section 4, we run simulations and present some results directly linking fat mass, diet and uric acid concentration and finally in Section 5, we discuss our models and provide possible directions for future study to gout symptoms.
2. Methodology
Gout is associated with increased blood levels of uric acid, called hyperuricemia. Uric acid is the end product of purine metabolism in humans due to the loss of uricase activity by the evolution of hominids, which leads to higher uric acid levels [17]. Although there is no universally accepted definition of hyperuricemia, it is often defined as an in-blood uric acid concentration greater than 7 mg/dL for men and greater than 6 mg/dL for women [18]. In women, high levels of uric acid are mainly found in postmenopausal women, African American, patients with renal disease and alcohol intake [19]. Physiologically, uric acid plasma concentrations increases with age; they are lower in women of childbearing age, but post menopause women reach similar values to those found in males [2].
Uric acid is formed by the liver and excreted mainly by the kidneys and intestines [20]. Several factors are associated as cause and consequences of high uric acid concentration, e.g. diet, obesity [17] [21] [22]. Observations show that uric acid concentration increases with insulin resistance, probably because hyperinsulinemia would cause lower renal uric acid excretion [23]. In addition, some studies found that uric acid serum concentrations are independently related to leptin concentration [24], thus suggesting it as a pathogenic factor responsible for uric acid increase in obese patients [25]. Furthermore, low muscle mass is negatively associated with the concentration of uric acid [26], meanwhile the oxidative stress produced by excessive uric acid can influence muscle mass reduction. The role played by a diet on high uric acid concentration has not yet been fully clarified, but high intake of fructose-rich industrialized food and high beer intake seem to influence uricemia [27]. The main mechanism to reduce the concentration of uric acid occurs by means of renal excretion, hence, glomerular function markers are positively associated with hyperuricemia. As a result, in this paper, we would like to study the relationship between Uric acid (U), Fat mass (F), Muscle mass (M), Leptin (L) and Insulin (I) using mathematical models.
3. Model Description
Our model is primarily based on the biological process proposed by Erick Prado de Oliveira and Roberto Carlos Burini [26]. Since uric acid, the end product of purine metabolism in humans, is mainly known for its harmful effects such as gout and uric lithiasis, as well as its association with renal disease, metabolic syndrome and hypertension [1] [28], in this paper we use uric acid as a biomarker for the development of gout symptoms. Generally, nutrients are assimilated and ingested in the bloodstream after eating. Once visceral adiposity reaches a threshold of lipid storage, the synthesis and secretion of leptin and insulin begins. In turn, leptin production and insulin resistance are able to increase the uric acid serum concentration since an increase in their synthesis inhibits uric acid excretion. Furthermore, fat mass and muscle mass reduction are associated with low-intensity chronic inflammation, which in turn promotes an increase in uric acid levels in order to protect the organism against the moderate oxidative stress.
Based on these mechanisms, we suppose that the main factors affecting the uric acid concentration are the following: fat mass (F), muscle mass (M), leptin (L) and insulin (I) (see Table 1 and Figure 1).
After incorporating the aforementioned mechanisms, we are able to derive the following system of differential equations and represent the specified interactions in Figure 1.
(1)
Figure 1. Schematic diagram of the mathematical model.
Table 1. Definition of state variables.
with initial conditions
The total mass is
, where B is the mass of bone. In this research we assume that B is a constant.
All parameter definition can be found in Table 2. For model (1), the equilibria and stability analysis are complex, hence we will be running simulations and presenting the results in the following section.
4. Simulation and Results
Based on simulations of our model combined with the knowledge that fructose is strongly associated with uric acid production, we now proceed to analyze the relationships between uric acid and energy intake. In particular, we want to explore the specific uric acid growth rates (
) stimulated by the proportion of fructose in the overall energy intake. Using the following parameters and initial conditions (see Table 3).
Table 2. Definition of the parameters.
Table 3. Average body composition by gender.
we were able to reach high levels of uric acid before losing all muscle mass.
Using the parameters corresponding to males between 40 and 59, we were able to determine that the amount of time it will take to surpass the 7 μg/dL threshold when
is about 10,000 days. Figure 2 shows the amount of time (in days) it will take to develop a high uric acid condition when
(green, blue and red respectively). Similarly, Figure 3, shows the same threshold for American women (20 - 39 years old) when
(green, blue and red respectively), suggesting that it would only take 5000 days and a production rate of
.
Figure 2. The relationship of uric acid concentration with the uric acid production rate as a result of energy intake in American males with ages 40 - 59.
Figure 3. The relationship of uric acid concentration with the uric acid production rate as a result of energy intake in American females with ages 20 - 39.
From the simulations, we are able to identify a strong correlation with initial fat mass in combination with uric acid production rate and the amount of time required to surpass the threshold concentration for both males and females. Nonetheless, studies show that only post-menopause women have comparable uric acid dynamics to men. This suggest that the energy intake stimulates lower levels of uric acid production when dealing with women in the first two age categories in Table 3 and thus an inhibition or an extra parameter directly impacting uric acid production or excretion needs to be incorporated.
5. Discussion and Future Work
Developed countries tend to have a higher burden of gout and with continuous improvements in living standards they seem to have an increasing prevalence and incidence of the disease [1]. This suggests that the implementation of mathematical models to study gout will become more and more important. The dynamic model we proposed in this paper, describes the concentration of uric acid in the blood as a result of fat mass and energy intake (diet); since a high uric acid condition is the first stage in the development of gout. In fact, due to the complexity of the biological mechanisms of the uric acid action pathway and energy balance in an individual, we had to use a simplified form of the biological process for the concentration of uric acid. Thus, there are still many ways to improve and better calibrate this process.
In addition to fat mass, muscle mass, insulin, leptin, and diet, there are many factors, such as blood pressure, triglycerides and inflammation, known to affect the concentration and production of uric acid. While the model in this paper ignores these factors that would likely alter the stability of the system, in the future, we are going to incorporate the effect these factors have on the production and excretion of uric acid.
Furthermore, it is known that muscle mass reduction is associated with increasing uric acid concentration levels. This is believed to be a defensive mechanism in order to protect the organism against the moderate oxidative stress associated with muscle mass decrease. In this paper, we only think about the bilinear functional effect between the interaction of muscle mass and uric acid, however, it has not yet been clarified what the cause or effect is and whether there are differences in these reactions when we have individuals of both genders.
In addition, uric acid concentrations are independently related to leptin concentration [25] and insulin resistance. Elevated levels of both would lower renal uric acid excretion rates [23]. While our model incorporates that effect of insulin and leptin on uric acid concentration using the function
, there are different functions that show this inhibiting influence that need to be considered, for example,
.
Finally, data collection could have a tremendous impact in the actual calibration of the model parameters, something that at this point is difficult to do since the dynamics of this process are not measured in gout patients. Nonetheless, further studies via mathematical models could help develop important hypotheses regarding this biological process, hence motivating in depth studies at the molecular level and thus the collection of appropriate data from individuals at different stages in the development and progression of gout.
Acknowledgements
We would like to thank the Mathematical and Theoretical Biology Institute (MTBI) at Arizona State University (ASU) for giving us the opportunity to conduct this project. This project was partially supported by grants from the National Science Foundation (DMS1263374), the National Security Agency (H98230-15-1-0021), the Office of the President of ASU, and the Office of the Provost at ASU. Yan’s work was partially supported by the National Natural Science Foundation of China (11671260) and Jiangxi Education Foundation (GJJ170816) of China.