Skip to main content

Transitioning from acute to chronic pain: a simulation study of trajectories of low back pain

Abstract

Background

Identifying how pain transitions from acute to chronic is critical in designing effective prevention and management techniques for patients’ well-being, physically, psychosocially, and financially. There is an increasingly pressing need for a quantitative and predictive method to evaluate how low back pain trajectories are classified and, subsequently, how we can more effectively intervene during these progression stages.

Methods

In order to better understand pain mechanisms, we investigated, using computational modeling, how best to describe pain trajectories by developing a platform by which we studied the transition of acute chronic pain.

Results

The present study uses a computational neuroscience-based method to conduct such trajectory research, motivated by the use of hypothalamic–pituitary–adrenal (HPA) axis activity-history over a time-period as a way to mimic pain trajectories. A numerical simulation study is presented as a “proof of concept” for this modeling approach.

Conclusions

This model and its simulation results have highlighted the feasibility and the potential of developing such a broader model for patient evaluations.

Background

This paper is a continuation of an earlier study of the same subject [16], where we focused on the transition of acute pain to chronic pain. As we mentioned, the need to distinguish common and reproducible pain trajectories in the population is of great importance. The 2011 Institute of Medicine (IOM) Report was released in response to the rising pain costs and prevalence, and sparked formation of a number of interdisciplinary teams and strategies to address the pain epidemic and determine areas of concern. Understanding and defining pain trajectories was one such area in which the knowledge was still currently lacking. As previously discussed in Gatchel et al. [16], identifying how pain transitions from acute to chronic is critical in designing effective prevention and management techniques for patients’ well–being, physically, psychosocially, and financially. Also, in our previous paper [16], we reviewed the various trajectories that have been delineated by different clinical research groups, which suggested the need for a comprehensive model for understanding them. We proposed that focusing this research on low back pain (LBP) is advantageous because LBP is the most prevalent form of musculoskeletal pain, and amasses billions of dollars in associated costs each year [16, 22]. Indeed, LBP presents an opportunity to observe all stages of pain (acute, subacute, chronic), as well as distinguishing different trajectories within this group.

Kongsted et al. [27] summarized ten (10) studies of LBP trajactories over a ten-year period from 2006 to 2015 [2, 3, 7, 8, 10, 13, 14, 23, 26, 36, 51] over research on adult patients of 10 cohorts by LBP trajectory methodology. In these studies, participants with a main complaint of LBP were followed from 3 to 12 months with data collection at four (4) to fifty-two (52) time-points. Outcome measures were LBP intensity, LBP frequency (number of LBP days per week) and activity limitation. Trajectory patterns were identified using either hierarchical cluster analysis, latent class analysis, or latent class growth analysis. From two (2) to twelve (12) discrete LBP trajectory patterns have been identified in these published studies. They suggested that trajectory differentiation between acute and chronic LBP is overly simplistic, and the next step is to shift from this paradigm to one that focuses on trajectories over time. Our proposed approach to modeling the time trajectory is well aligned with this new goal to build upon these earlier pain trajectory studies.

As is well known, pain pathways include multiple brain regions, and the connection among the anterior cingulate cortex, the parieto-insular cortex, the thalamus, as well as the amygdala. Several neuronal-processing mechanisms of pain signals have been proposed, although their temporal behavior, especially during an extended period of several weeks or months, is still not fully understood. It has been known that different brain areas are involved in neuronal electric and chemical activities in response to pain, and form a distributed pain-processing network, mostly centered on the somatosensory cortex and the thalamic axis, closely associated with the pain signals and/or pain-induced stress. While a full study of pain and its transition between general acute and chronic pain is beyond our reach at this point, due to the complexity of the brain network and brain activities, LBP involves relatively isolated areas of the brain, and has been attributed to a correlation to HPA-axis activity [18, 20]. Although the sequence of brain-activity events leading to LBP can still be complex or even convoluted, it is generally agreed that the principle regions of the pain response are localized in the paraventricular nucleus of the hypothalamus, the anterior lobe of the pituitary gland, and the adrenal gland, commonly referred to as the HPA axis [45]. The HPA axis plays an important role in balancing hormonal levels for the brain, and generates high concentrations of hormones in response to pain (considered as a form of stress), which leads to many “downstream” changes [4]. A number of measurement data, such as patient self-reported trajectory data, EEG data, and cortisol level changes, can be used for the quantitative study of the HPA axis, which is the primary neuronal responding mechanism for stress and pain. From the quantitative measurements, we can utilize mechanism-based computational models to predict the temporal behavior of hormone concentrations and infer trends in pain trajectories.

Along with cortisol, adrenocorticotropin (ACTH) is one of the major hormones secreted by the HPA in the anterior pituitary region in response to severe pain and other stressors [31, 63]. Tennant et al. [53] measured ACTH serum levels in fifty-five (55) severe, chronic pain patients. This study supports other reports that pituitary–adrenal function may be altered during uncontrolled pain period and after return to normal when pain control is achieved [20, 50]. Corticotropin-releasing factor (CRF) is also released from the hypothalamus and in widespread areas of the brain following the stress or pain episode. Lariviere and Melzack [30] presented evidence that CRF can act at all levels of the neuraxis to produce analgesia; inflammation must be present for local CRF to evoke analgesia and the analgesic effects of CRF presented significant specificity for prolonged pain. The similar mechanism is also studied in human subjects [6]. Overall, pain results in a hyperarousal of the hypothalamic–pituitary–adrenal system which results in elevated serum hormone levels such as adrenocorticotropin, cortisol, and pregnenolone [52].

The current methodology of measuring hormones in the hypothalamic pituitary adrenal (HPA) axis was reviewed recently by Yeo et al. [62]. Different types of dexamethasone suppression testing are compared and described in detail in [61, 62]. Common serum cortisol measurements are performed by Mass Spectrometry, Beckman Coulter, Roche Diagnostics, Siemens ADVIA Centaur, Siemens Immulite, Abbott Architect, Vitros, Tosoh Bioscience AIA-PACK Test Cups, and urinary cortisol measurements are by Mass Spectrometry, Beckman Coulter, Roche Diagnostics, Siemens ADVIA Centaur, Abbott Architect, Vitros. Salivary cortisol is measured by Mass Spectrometry, Roche. ACTH measurements however have been less common, and they can be influenced by patients with cortisol producing adrenal adenomas or patients taking exogenous steroids. The preferred specimen for ACTH is plasma, and the measurements are done by Siemens Immulite and Roche. Urine measurement can also be done [20]. CRF (or Corticotropin-releasing Hormone, CRH)—measurement is not common. Although it is known that high levels CRH are associated with high levels of CRH binding protein, few data is available in setting up a reliable measurement protocol in this area. The evaluations of the HPA axis functions are usually through biochemical measurements, imaging studies can complement the hormonal evaluations, providing valuable information for prognosis and management. Pituitary T2 MRI, Pituitary CT scan have been used in combination with biochemical measurement [54], and Adrenal CT scan, MRI have been used to adrenal measurements in Cushing syndrome [42]. In summary, cortisol measurements remain a superior method for HPA function evaluations for convenience and reliability. In particular, the modeling work proposed here is based on continuous time data and cortisol measurements using electrochemical impedance sensor [56] provides a platform for data acquisition.

Computational neuroscience-modeling techniques (using systems of ordinary differential equations) have been utilized to develop a reliable predictive outcome model for HPA-related issues, such as ultradian and circadian patterns in normal, depressed, post-traumatic stress disorder states, and their comorbid pain [1, 47, 49, 59]. In Prince et al. [40], the authors modeled acute pain from a gating mechanism point of view, which assumes the flow of inputs through the spinal cord to the brain as a major contributor for the pain experience, and constructed a biologically plausible mathematical model. Research efforts in the past have been also focused on developing mechanistic models for stress and pain, using differential equations and computational simulation of HPA-axis activities [43]. They found similar quantitative trajectories as LBP. Therefore, one can then use computational simulation and bifurcation analysis to study the complex biological processes that are involved, complementary to experimental work. A bifurcation study is a computational investigation of trajectory patterns in a system’s parametric space, and is a powerful tool useful for investigating possible pathways of the HPA network accounting for variations in neural receptors, synaptic plasticity, as well as other factors such as neuron-conductivity degeneration. Indeed the HPA-axis process and its abnormality are highly associated with the pain and stress of muscular movement involved in LBP. The present proposed predictive-outcome model is a first step in our efforts to understand the acute/chronic pain transition mechanism and a proof of concept to show the feasibility of applying the computational model to pain transition, which is a high-impact medical research issue.

Methods

Queries

In a clinical study of 131 persons with chronic widespread pain, 267 ‘at risk’ and 56 controls [34] with abnormal HPA-axis function were associated with a high rate of fibromyalgia, a syndrome characterized by chronic widespread body pain. These results confirmed the hypothesis that the psychosocial-stress caused by pain, as well as the pain, are linked to altered HPA-axis functioning. Increased cortisol levels due to alterations in the regulation of the HPA axis at young age are also believed to be one mediating mechanism for several adult conditions, including metabolic syndrome, cardiovascular disease, and psychiatric disorders [34]. Infants born after significant exposure to stressful conditions are often small for gestational age (SGA), based on standardized growth norms. A study of 37 participants, including infants of gestational age ranging from 34 to 41 weeks, suggests SGA neonates have blunted HPA axis responses to stressors [37]. These findings are consistent with animal models showing that adverse intrauterine conditions can result in blunted cortisol responses to acute stressors, and may provide a mechanism for adult susceptibility to disease through HPA axis abnormality. The importance of the hormonal response system to stress through the HPA axis is also seen in other clinical studies, particularly stress hormones known as glucocorticoids (and primarily cortisol) that are key factors in patient alcohol dependence [48].

There is growing evidence that relative hypocortisolism, as a marker of stress-induced HPA axis dysfunction, may increase vulnerability to pain and chronic pain disorders [17,18,19, 28, 29, 35, 58]. HPA activities, especially cortisol trajectory, are also a useful tool in longitudinal cohort studies of pain. In a study by Paananen et al. [38], for example, their sample group included 805 participants from the Western Australian Pregnancy Cohort (Raine) Study, who participated in the Trier Social Stress Test (TSST) at age 18-years. The number of pain sites, pain duration, pain intensity and pain frequency were assessed at age 22 in order to measure severity of musculoskeletal (MS) pain. An abnormal HPA axis response to psychosocial stress at 18-years of age was shown to be associated with MS pain alone, and MS pain combined with increased pain sensitivity at 22-years of age.

More recent evidence also suggests that in a disordered stress-system pathway, the HPA axis activity may be responsible for the genesis and maintenance of long-term sensory and emotional problems that lead to post-traumatic pain and disability. In a study by [57], the authors used hair-cortisol and hair-normalized salivary cortisol as biomarkers of distress following traumatic injuries of whiplash or distal radius fractures. Small sample results indicated that the cortisol-waking response may become a useful biomarker of current distress as measured using the Pain Catastrophizing Scale [57], especially when it is normalized to three-month hair cortisol. The hair-normalized cortisol waking-response also had predictive capacity, correlating with three-month self-reported disability. Tomas et al. [55] also provided a review of HPA axis dysfunction in chronic fatigue syndrome patients. The study included evidence of enhanced corticosteroid-induced negative feedback, basal hypocortisolism, attenuated diurnal variation, and a reduced responsivity to challenge. A putative causal role for a genetic profile, childhood trauma, and oxidative stress were considered. In addition, gender was also determined as a factor, in addition to an increased frequency of HPA axis dysregulation in females. The mechanisms of the HPA appears to affect LBP by following similar pathways as chronic fatigue syndrome. Abnormal cortisol concentrations reflect differences in the biological mediation of the stress response, or may be consequent on the differential nature/magnitude of the stressor engendered by the LBP. For example, in the chronic fatigue study, an attenuated diurnal variation has also been shown [39, 55], particularly with a loss of the morning-peak of ACTH or cortisol, while challenge studies often, but not invariably, show a diminished HPA axis responsivity. This has been assessed using ACTH, cortisol, and/or 11-deoxycortisol response to pharmacological challenge employing, for example, dexamethasone combined with corticotropin-releasing hormone (CRH), insulin, inflammatory cytokines, metyrapone, a psychosocial-challenge (e.g., using the Trier Social Stress Test), and to a physiological-challenge (such as a wakening; Tomas et al. [55]).

In the proposed model, we focus on the cortisol-level as a biomarker for pain trajectory for the continuation study over the time. There are a few current in vivo techniques that allow us to monitor the cortisol-level continuously for a long duration of time. For example, [56] used an electrochemical impedance sensor to measure cortisol concentration in the interstitial fluid of a human subject. The interstitial fluid is extracted by means of vacuum pressure from micropores created on the stratum corneum layer of the skin. Other measurements of the cortisol such as the saliva sampling protocol [46], and electrochemical immunosensing platforms that are well-known for their sensitive and selective detection of cortisol in biofluids [44], have also developed as tools for continuous in vivo measurements. They will guide predictive modeling development. The recent nanosheet technology [24] also allows monitoring cortisol in low-volume perspired human sweat (cortisol dynamic range from 1 to 500 ng/mL with a limit of detection of 1 ng/mL) using electrochemical impedance spectroscopy.

In this study, we will focus on the longitudinal study of LBP, in particular the transition from acute pain to chronic pain. The HPA activity history over time is of high interest. Our modeling strategy calls for a coupled system of ordinary differential equations to represent the network of brain regions along the HPA axis, and its cortisol and adrenaline production. Specifically relating to a trajectory modeling study in LBP, we focus on outcome measures such as LBP intensity, LBP frequency (number of LBP days per week), but will ignore activity-limitation which will be a subject of a future study. Trajectory patterns in a group of experimental subjects were identified using Hierarchical Cluster Analysis, Latent Class Analysis, or Latent Class Growth Analysis [27].

Modeling consideration

Following Sriram et al. [47], in developing the mathematical model of the cortisol molecular network for the HPA, we have made two assumptions: (i) The first order dilution-rate due to the transport of hormones and autonomous degradation are considered together. Apart from dilution/autonomous degradation, Michaelis–Menten kinetics are separately considered for the degradation of the hormones, and hormone complexes, within each specific region of the brain (hypothalamus, pituitary and adrenal); and (ii) A sufficient number of molecules is present for the reactions to take place, using continuum kinetics so that stochastic fluctuations (internal noise) are minimal.

Hypothalamic–pituitary–adrenal (HPA) axis

The anatomical structure that mediates the stress response is mainly in the central nervous system (CNS) and its peripheral organs. The principle-effectors of the stress-response are localized in the paraventricular nucleus of the hypothalamus, the anterior lobe of the pituitary gland, and the adrenal gland. This collection of brain areas is commonly referred to as the Hypothalamic–pituitary–adrenal (HPA) axis [45]. The HPA axis plays an important role in balancing hormonal levels for the brain, and generate high concentrations of hormones in response to stress, which leads to “downstream changes” [43]. In response to stress over a period of time, the paraventricular nucleus of hypothalamus, which contains neuroendocrine neurons, releases corticotropin-releasing hormone (CRH). The anterior lobe of the pituitary gland is stimulated by CRH to secrete adrenal corticotrophic hormone (ACTH). The adrenal cortex then produces cortisol hormones in response to stimulation by ACTH. Cortisol is a major stress-related hormone, and has effects on many tissues in the body, particularly in the brain. In the brain, cortisol acts on two types of receptors: mineralocorticoid receptors and glucocorticoid receptors [4, 9, 41]. For “down regulating” CRH, it is known that cortisol inhibits the secretion of CRH through the glucocorticoid-receptors complex [11]. With a strong affinity to the glucocorticoid-receptors complex, cortisol, in turn, acts on the hypothalamus and pituitary in a negative feedback cycle to “down regulate” the production of cortisol [47]. Figure 1 illustrates the conceptual connection of the different brain areas in this process. Based on the biophysical relation, we established the underlying nonlinear ordinary differential equations for the cortisol network, as presented in Tables 1, 2, 3. Table 1 provides the mathematical model as in a system of ordinary differential equations. Table 2 provides the biological meaning of the parameters. A set of HPA parameters for LBP-modeling is provided in Table 3 for LBP modeling study. The model described in Table 1 is constructed by the chemical kinetics relations shown in Fig. 1. The four brain areas have excitatory or inhibitory inputs to other areas through mostly hormonal chemical kinetics relations and synaptic couplings [47]. The coefficients [47] are either from: the scientific literature; obtained by other researchers through direct measurement; or are best estimates to match experimental data through global optimization [60].

Fig. 1
figure 1

A Diagram of Hypothalamic–Pituitary–Adrenal (HPA) axis. In response to stress, the hypothalamus releases CRH, which activates the pituitary which secretes ACTH. ACTH stimulates the adrenal to secrete glucocorticoid. Glucocorticoid performs a negative feedback on the secretion of CRH and ACTH

Table 1 Hypothalamic pituitary adrenocortical (HPA) axis model
Table 2 A list of model variables
Table 3 The model-kinetic parameters

Computer simulation technique

We simulated the time-series data of cortisol level during a period of 100 days, based on a computational HPA-axis model, which is described in Tables 1 and 3. MATLAB solver ode45 (https://www.mathworks.com/help/matlab/ref/ode45.html) was used for its accuracy and relatively fast speed in computing this model. It should be pointed out that, in our preliminary simulation study to demonstrate how mathematical modeling techniques can be used to develop a better understanding of these trajectory results, it was not our intention to explore all possible scenarios.

Results

Preliminary simulation study

In this section, we present some examples based on the cortisol-dynamics model reported by [47], with parameters modified to fit into LBP problems. We then conducted a series of modeling pain-trajectories, based on cortisol that revealed patterns that resemble the pain trajectory patterns reported by Kongsted et al. [26]. Because cortisol secretion can be monitored in real-time, its time trajectory provides insight for studying LBP due to the inherent structure of the HPA in the brain hormone system under stress. The cortisol-level values are oscillatory as the time goes on, in response to the pain-process during LBP.

The time trajectory of cortisol-level is presented in Fig. 2, during the first 100 days for the HPA computational model. The cortisol values are given at each hour. Cortisol data can be simulated for any instance, but we only present the data at each hour, and each day contains 24-cortisol data points. Due to the inherent structure of the HPA brain hormone system under stress, the cortisol level values are oscillatory as the time goes forward, in response to the pain-process during LBP. The values are between 0.13440 and 0.13465 μ/DL for the entire times, and the up and down cycle is about 24 h. This simulated cortisol data can be compared with clinical data in the future for validation purpose, whereas cortisol secretion of patients can be monitored in real time and provide further insight for LBP and pain trajectory patterns.

Fig. 2
figure 2

Simulated data of cortisol level during a period of 100 days, based on a HPA computational model that is described in Tables 1 and 3. The x-axis represents the time (unit: day) and the y-axis is the cortisol level. MATLAB solver ode45 (https://www.mathworks.com/help/matlab/ref/ode45.html) was used for its accuracy and relatively fast speed in computing this model. Due to the inherent structure of the HPA brain hormone system under stress, the cortisol level values are oscillatory as the time goes forward, in response to the pain-process during LBP. Cortisol data can be simulated for each hour, and each day contains 24-cortisol data points. This can be compared with clinical data in the future, where cortisol secretion of patients can be monitored in real time and provide further insight for LBP and pain trajectory patterns. The mathematical HPA model, and its biological meaning of the parameters, are provided in Tables 1, 2. A typical set of HPA parameters for LBP is given in Table 3

In Fig. 3, a simulated pain trajectory, based on HPA-modeling analysis, is presented.

Fig. 3
figure 3

A simulated 2-state pain trajectory, based on the HPA modeling. The cortisol level values obtained from HPA model are classified into two different groups, and they present high-state and low-state patterns. One can use the higher cortisol level to model the high intensity pain, and the lower level to represent low-pain or no-pain. We first take a daily average of cortisol level (24 hourly data points). A subjective threshold (0.1345 in this model; around the average of extreme cortisol values) of cortisol level is used to separate the first 100 days into two groups: high-pain day (labeled in the high group as “cort = 0.13463”) for the days when cort 0.1345, or low-pain day when cort 0.1345 (labeled in the low group as “cort = 0.13442”) for trajectory study. The simulated pain trajectory is representative for a LBP episode that has chronic and intermittent pain

The cortisol-level values were classified into two different groups, and they present high-state and low-state patterns. The higher-cortisol level was used to model the high-intensity pain, and the lower-level to represent the low-pain or no-pain. We took a daily average of cortisol values for the study period. Then, for the averaged cortisol values, a subjective threshold of cortisol level was taken to separate the daily pain states into two groups (high-pain or low-pain) for the trajectory study. The simulated pain trajectory in two states (high-pain and low-pain) is presented in Fig. 3. We first took a daily average of cortisol level (24 hourly cortisol data points). A subjective threshold (0.1345 in this model; around the average of extreme cortisol values) of cortisol level was used to separate the first 100 days into two groups: high-pain day for these days when cort  0.1345, or low-pain day when cort < 0.1345, for pain trajectory study. The simulated pain trajectory depicted in Fig. 3 is representative for a LBP episode that has chronic and intermittent pain.

The HPA-system parameters that reflect the synaptic connectivity are crucial in setting pain-states. One can also analyze the dynamics-pattern changes induced by any parameter variation, in order to understand LBP transition from acute to chronic pain. The HPA-system parameters reflect the synaptic coupling strength or neuron degradation, and they are crucial in setting pain-states. One can systematically analyze the dynamic pattern-changes induced by any parameter variation, in order to understand LBP transition from acute to chronic pain.

In Fig. 4a, b, a case of cortisol dynamics of acute pain was simulated over a 100-day period, with acute pain that quickly subsides in 1–2 weeks, whereas a HPA with a reduced adreno-corticotrophin hormone degradation was considered (Vs4 was reduced from 15 µg d to 0.907 µg d, while others remain at normal level). Figure 4a depicts the cortisol level in the 100 days. The cortisol level quickly returns to a flat level in Fig. 4a, after a few initial oscillations. In Fig. 4b, we presented a 3-state pain trajectory calculated from simulated cortisol data in Fig. 4a. Similarly, we first took a daily average of cortisol level (24 hourly data points). Subjective thresholds (0.23986 and 0.2985 in this model) of cortisol level are used to separate the days into three groups: high-pain day (labeled in the high-pain group on the top) for the days when cort  0. 23986, or healed day when 0.23986 > cort  0. 23985 (labeled in the healed or no-pain group in the middle), and low-pain for cort < 0. 23985, for trajectory studies. The simulated pain trajectory is representative for an acute LBP episode, and the corresponding pain trajectory converges quickly to the healed or no-pain state, as shown in Fig. 4b.

Fig. 4
figure 4

A simulated cortisol level data and pain trajectory for acute pain. a A simulated cortisol level of acute pain over the 100-day period. Here, we present a case of simulated cortisol dynamics that is with acute pain that quickly subsides in 1–2 weeks. We modify a HPA model with a reduced adreno-corticotrophin hormone degradation rate (Vs4 is reduced from 15 to 0.907 µg days), while others remain at a normal level. The cortisol level quickly returns to a flat level in a, after a few initial oscillations. b Simulated 3-state pain trajectory calculated from simulated cortisol data in a. The cortisol level values obtained from HPA model are classified into three different groups, and they present high pain state, low pain state and no pain (healed) state patterns. Again, we first take a daily average of cortisol level (24 hourly data points). Subjective thresholds (0.23986 and 0.2985 in this model) of cortisol level are used to separate the first 100 days into three groups: high-pain day (labeled in the high-pain group) for the days when cort  0. 23986, or healed day when 0.23986 > cort  0. 23985 (labeled in the healed or no-pain group), and low-pain for cort < 0. 23985 for trajectory studies. The simulated pain trajectory is representative for an acute LBP episode and the corresponding pain trajectory will converge quickly to the healed or no-pain state

Then, in Fig. 5, an HPA-system of “manic” pain was simulated, when the stress is elevated (Kstress is increased from 10.1 to 30 µg days), with other parameters identical as in Fig. 4; the cortisol trajectory transforms into a “manic” oscillatory pattern, as shown in Fig. 5a. The nearly periodic patterns have a period around 10, and alternate between high-amplitude oscillations and low-amplitude ones during the 100 day time period. In Fig. 5b, pain trajectory calculated from data in Fig. 5a shows a similar manic pattern. The cortisol level values obtained from HPA model are classified into three different groups, and they present high pain state, low pain state and no pain (healed) state patterns. We first took a daily average of cortisol level (24 hourly data points). Subjective thresholds (0.23990 and 0.29980 in this model) of cortisol level were used to separate the first 100 days into three groups: high-pain day (labeled in the high-pain group in the top) for the days when cort  0. 23990; or healed day when 0.23990 > cort  0. 23980 (labeled in the healed or no-pain group in the middle); and low-pain for cort < 0. 23980 (in the bottom) for trajectory studies. The simulated pain trajectory is representative for a manic-style pain LBP.

Fig. 5
figure 5

The cortisol trajectory data for a manic-style pain. a A simulated cortisol level of manic-style pain over the 100-day period. This is simulated with a patient who can quickly recover from acute pain with identical parameters as in Fig. 3, when the stress is elevated (K-stress is increased from 10.1 to 30 µg days). The cortisol trajectory can transition from acute pain to a “manic” oscillatory pattern, based on simulated cortisol level data for the HPA model with the new parameter set. The cortisol level is presented here. The nearly periodic patterns have a period around 10, and alternate between high-amplitude oscillations and low-amplitude ones during the 100 days time period. b Pain trajectory calculated from data in a for 100 days of the HPA model simulation. The cortisol level values obtained from HPA model are classified into three different groups, and they present high pain state, low pain state and no pain (healed) state patterns. We take a daily average of cortisol level (24 hourly data points). Subjective thresholds (0.23990 and 0.29980 in this model) of cortisol level are used to separate the first 100 days into three groups: high-pain day (labeled in the high-pain group in the top) for the days when cort  0. 23990, or healed day when 0.23990 > cort  0.23980 (labeled in the healed or no-pain group in the middle), and low-pain for cort < 0.23980 (in the bottom) for trajectory studies. The simulated pain trajectory is representative for a manic-style pain LBP

Pain trajectory of patient groups

We now a present a computer-simulated study of pain trajectories in a group of patients. Unlike an individual patient whose hormonal levels are uniquely decided, we expect a distribution of cortisol levels among different patients. We assume the cortisol dynamics values are normalized (as in [57], so the discrepancy is due to random fluctuation or variation among individual patients. We chose a sample of 100 random values of a reduced adreno-corticotrophin hormone degradation rate (Vs4) and the stress coefficient (Kstress). We first evaluated the cortisol dynamics, and took the averaged values of 100 patients. Then, we developed the pain trajectory for this group of patients. We also calculated the variance of the cortisol-level and pain-trajectories. Figure 6a and b depict the average dynamics and their variances in cortisol level during the first 100 days, We kept all the same parameters in the HPA model as in Fig. 4, except Kstress, and we choose 100 randomized values for Kstress. The 100 values were randomly generated by the computer; the average Kstress is 0.9922, and the variance is 0.2907. Then, we simulated the HPA model for each Kstress. The cortisol trajectories vary among patients, from high-amplitude oscillations and low-amplitude ones during the 100-day time period. Each patient’s cortisol data is presented in a, with different color as shown in Fig. 6a. We also calculated the average daily cortisol-levels to be used to calculate the group-pain trajectory. The level variance is also marked in the group. Both daily values are presented in b. For the group, the average and variance become stabilized after 60 days. The average cortisol level values were then used to obtain a 3-state pain trajectory. We took a daily average of cortisol level (24 hourly data points). Subjective thresholds (0.556 and 0.552 in this model) of cortisol level were used to separate the first 100 days into three groups: high-pain day (labeled in the high-pain group in the top) for the days when cort  0.556; or healed day when 0.556 > cort  0.552 (labeled in the healed or no-pain group in the middle); and low-pain for cort < 0.552 (in the bottom) for trajectory studies. Shown in Fig. 6c, the simulated pain trajectory is representative for a group of LBP patients who recovered from acute pain in various severities.

Fig. 6
figure 6figure 6

The cortisol trajectory and pain trajectory modeling for a group of 100 patients. a Simulated cortisol level for LBP by the HPA model in a group of 100 patients. We keep all same parameters in the HPA model as in Fig. 4 except Kstress, and we choose 100 randomized values for Kstress. The twenty values are randomly generated by computer, the average Kstress is 0.9922, and the variance is 0.2907. Then we simulated the HPA model for each Kstress. The cortisol trajectories vary among patients, from high-amplitude oscillations and low-amplitude ones during the 100-day time period. Each patient’s cortisol data is presented in a with different color. b Average cortisol level for LBP by the HPA model in a group of 100 patients. We also calculated the average daily cortisol-levels to be used to calculate the group-pain trajectory. The level variance is also marked in the group. For the group, the average and variance become stabilized after 60 days. c Pain trajectory calculated from data in Fig. 4a for 100 days of the patients group. The average cortisol level values obtained from HPA model are classified into three different groups, and they present high pain state, low pain state and no pain (healed) state patterns respectively. We take a daily average of cortisol level (24 hourly data points). Subjective thresholds (0.556 and 0.552 in this model) of cortisol level are used to separate the first 100 days into three groups: high-pain day (labeled in the high-pain group in the top) for the days when cort  0.556, or healed day when 0.556 > cort  0.552 (labeled in the healed or no-pain group in the middle), and low-pain for cort < 0.552 (in the bottom) for trajectory studies. The simulated pain trajectory is representative for a group of LBP patients who recovered from acute pain in various severities

Discussion

The simulations revealed in the present study demonstrate the feasibility of studying pain trajectory and pain transition, based on cortisol dynamics of a HPA-computational model. From a representative cortisol-value change, we can construct a fairly good pain trajectory for an individual patient, as well as a group of patients. With variation of parameters representing synaptic connectivity and neural degrading, a cortisol dynamics and pain trajectory study can provide various scenarios of acute, manic-style pain and intermittent/chronical pain. However, one will need to “fine-tune” the parameters with other pain-related measures and biomarkers in the future. Our ultimate goal is to use computational-bifurcation analysis in order to predict the outcome of the initial pain trajectory, and then the transition to different types of trajectories (i.e., acute pain and chronic pain in LBP), based on stable-patterns in different parameter-sets. This mathematical tool has been very successful in predicting the pattern-transition during epileptic seizures in an extended Taylor model (Fang et al. in press), and breathing irregularity induced by the changes in Pre-Bőtzinger Complex (a pacemaker in deep-brain regions; [12, 15].

It has been established that interdisciplinary treatment strategies are generally most effective for pain management, allowing for the customization of a treatment plan for the individual patient [5]. Understanding how to identify pain trajectories would be of great importance in tailoring treatment-strategies to patients experiencing pain, both at the acute and chronic stages. Kongsted et al. [26] iterate that other than known factors, such as activity-limitation, work participation, history of back or leg pain, anxiety, and catastrophizing, there were also observed relationships between pain patients and mixed-evidence concerning sleep-disturbances and patient-education level. They further note that these variables were not used to group patients into trajectory patterns. It would be useful to determine if these associations will affect the trajectory a patient follows long-term. Our computational modeling is currently not sophisticated enough to include all these factors. Fortunately, there are currently available cutting-edge methods, such as daily electronic diaries using smart phones (e.g., [21, 23, 25, 32, 33] to systematically track patients from the acute LBP injury and over the following one-year period, using an inception-cohort study design. Finally, these methods will supplement the HPA-axis predictive data to further provide guidance on biomarkers that may be used to better understand the underlying neuroscientific-transition mechanisms involved in order to help prevent CLBP.

Conclusions

Computational-modeling based on HPA-axis can also provide insight into the mechanisms of pain, and the pathways from acute injury to chronic conditions, and collaborate with pain trajectory patterns. There are, though, currently limitations of predictive modeling, as it primarily depends on HPA-hormonal dynamics only. However, illuminating pain trajectories and subsequently developing effective and customizable prevention and intervention techniques for the growing adult population is critical [16]. Many trajectory patterns have been proposed in the CLBP population, and we must develop a cohesive understanding of how and why these patterns are emerging in order to provide the best possible medical care for those suffering with pain. Promoting ideal outcomes, particularly in the aging population, requires all-inclusive comprehension of risk factors and triggers.

Availability of data and materials

All data generated or analyzed during this study are included in this published article.

References

  1. Andersen M, Vinther F, Ottesen J. Mathematical modeling of the hypothalamic-pituitary-adrenal gland (HPA) axis, including hippocampal mechanisms. Math Biosci. 2013;246(1):122–38. https://doi.org/10.1016/j.mbs.2013.08.010 (Epub 2013 Sep 4).

    Article  CAS  PubMed  Google Scholar 

  2. Axen I, Bodin L, Bergstrom G, Halasz L, Lange F, Lovgren PW, Rosenbaum A, Leboeuf-Yde C, Jensen I. Clustering patients on the basis of their individual course of low back pain over a six month period. BMC Musculoskelet Disord. 2011;12:99.

    Article  Google Scholar 

  3. Axen I, Leboeuf-Yde C. Trajectories of low back pain. Best Pract Res Clin Rheumatol. 2013;27(5):601–12. https://doi.org/10.1016/j.berh.2013.10.004.

    Article  PubMed  Google Scholar 

  4. Besedovsky H, Chrousos G, Del Rey A. The hypothalamus-pituitary-adrenal axis, 1st edition. 2008; Amsterdam: Elsevier. ISBN 9780444530400.

    Chapter  Google Scholar 

  5. Bevers K, Watts L, Kishino ND, Gatchel RJ. The biopsychosocial model of the assessment, prevention, and treatment of chronic pain. US Neurol. 2016;12(2):98–104. https://doi.org/10.17925/usn.2016.12.02.98.

    Article  Google Scholar 

  6. Carpenter L, Tyrka AR, McDougle CJ, Malison RT, Owens MJ, Nemeroff CB, Price LH. Cerebrospinal fluid corticotropin-releasing factor and perceived early-life stress in depressed patients and healthy control subjects. Neuropsychopharmacology. 2004;29:777–84.

    Article  CAS  Google Scholar 

  7. Chen C, Hogg-Johnson S, Smith P. The recovery patterns of back pain among workers with compensated occupational back injuries. Occup Environ Med. 2007;64(8):534–40.

    Article  Google Scholar 

  8. Deyo RA, Bryan M, Comstock BA, Turner JA, Heagerty P, Friedly J, Avins AL, Nedeljkovic SS, Nerenz DR, Jarvik JG. Trajectories of symptoms and function in older adults with low back disorders. Spine (Phila Pa 1972). 2015;40(17):1352–62.

    Article  Google Scholar 

  9. Douglas AJ. Central noradrenergic mechanisms underlying acute stress responses of the Hypothalamopituitaryadrenal axis: adaptations through pregnancy and lactation. Stress. 2005;8(1):5–18. https://doi.org/10.1080/10253890500044380.

    Article  CAS  PubMed  Google Scholar 

  10. Downie AS, Hancock MJ, Rzewuska M, Williams CM, Lin CC, Maher CG. Trajectories of acute low back pain: a latent class growth analysis. Pain. 2015;157:225–34.

    Article  Google Scholar 

  11. Drouin J, Sun YL, Tremblay S, Lavender P, Schmidt TJ, de Lan A, Nemer M. Homodimer formation is rate-limiting for high-affinity DNA binding by glucocorticoid receptor. Mol Endocrinol. 1992;6(8):1299–309.

    CAS  PubMed  Google Scholar 

  12. Duan L, Chen Z, Tang X, Su J. Bursting transition dynamics within the pre-bőtzinger complex. Int J Bifurcat Chaos. 2016;26(10):1650162. https://doi.org/10.1142/S0218127416501625.

    Article  Google Scholar 

  13. Dunn KM, Jordan K, Croft PR. Characterizing the course of low back pain: a latent class analysis. Am J Epidemiol. 2006;163(8):754–61.

    Article  Google Scholar 

  14. Dunn KM, Campbell P, Jordan KP. Long-term trajectories of back pain: cohort study with 7-year follow-up. BMJ Open. 2013;3(12):e003838.

    Article  Google Scholar 

  15. Fan D, Wang Q, Su J, Xi H. Stimulus-induced transitions between spike-wave discharges and spindles with the modulation of thalamic reticular nucleus. J Comput Neurosci. 2017;43:203. https://doi.org/10.1007/s10827-017-0658-4.

    Article  PubMed  Google Scholar 

  16. Gatchel RJ, Bevers K, Licciardone JC, Su J, Du Y, Brotto M. Transitioning from acute to chronic pain: an examination of different trajectories of low-back pain. Healthcare. 2018;6(48):1–12. https://doi.org/10.3390/healthcare6020048.

    Article  Google Scholar 

  17. Geiss A, Rohleder N, Kirschbaum C, Steinbach K, Bauer HW, Anton F. Predicting the failure of disc surgery by a hypofunctional HPA axis: evidence from a prospective study on patients undergoing disc surgery. Pain. 2005;114:104–17. https://doi.org/10.1016/j.pain.2004.12.007.

    Article  CAS  PubMed  Google Scholar 

  18. Generaal E, Vogelzangs N, Macfarlane GJ, Geenen R, Smit JH, Penninx BW, Dekker J. Reduced hypothalamic–pituitary–adrenal axis activity in chronic multi-site musculoskeletal pain: partly masked by depressive and anxiety disorders. BMC Musculoskeletal Disord. 2014;15:227. https://doi.org/10.1186/1471-2474-15-227.

    Article  Google Scholar 

  19. Godfrey KM, Strachan E, Dansie E, Crofford LJ, Buchwald D, Goldberg J, Poeschla B, Succop A, Noonan C, Afari N. Salivary cortisol and cold pain sensitivity in female twins. Ann Behav Med. 2014;47:180–8. https://doi.org/10.1007/s12160-013-9532-4.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Griep EN, Boersma JW, Lentjes EG, Prins AP, van der Korst JK, de Kloet ER. Function of the hypothalamic–pituitary–adrenal axis inpatients with fibromyalgia and low back pain. J Rheumatol. 1998;25:1374–81.

    CAS  PubMed  Google Scholar 

  21. Jamison RN, Raymond SA, Slawsby EA, McHugo GJ, Baird JC. Pain assessment in patients with low back pain: comparison of weekly recall and momentary electronic data. J Pain. 2006;7(3):192–9. https://doi.org/10.1016/j.pain.2005.10.006.

    Article  PubMed  Google Scholar 

  22. Katz JN. Lumbar disc disorders and low back pain: socioeconomic factors and consequences. J Bone Joint Surg Am. 2006;88(suppl 2):21–4.

    PubMed  Google Scholar 

  23. Kent P, Kongsted A. Identifying clinical course patterns in SMS data using cluster analysis. Chiropr Man Therap. 2012;20(1):20.

    Article  Google Scholar 

  24. Kinnamon D, Ghanta R, Muthukumar K, Lin S, Prasad S. Portable biosensor for monitoring cortisol in low-volume perspired human sweat. Sci Rep. 2017;7:13312. https://doi.org/10.1038/s41598-017-13684-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Kongsted A, Leboeuf-Yde C. The Nordic back pain subpopulation program—individual patterns of low back pain established by means of text messaging: a longitudinal pilot study. Chiropractic Osteopathy. 2009. https://doi.org/10.1186/1746-1340-17-11.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Kongsted A, Kent P, Hestbaek L, Vach W. Patients with low back pain has distinct clinical course patterns that were typically neither complete recovery nor constant pain. A latent class analysis of longitudinal data. Spine J. 2015;15(5):885–94.

    Article  Google Scholar 

  27. Kongsted A, Kent P, Axen I, Downie AS, Dunn KM. What have we learned from ten years of trajectory research in low back pain? BMC Musculoskeletal Disord. 2016;17(220):1–11. https://doi.org/10.1186/s12891-016-1071-2.

    Article  Google Scholar 

  28. Kuehl LK, Michaux GP, Richter S, Schanchinger H, Anton F. Increased basal mechanical pain sensitivity but decreased perceptual wind-up in a human model of relative hypocortisolism. Pain. 2010;149:539–46. https://doi.org/10.1016/j.pain.2010.03.026.

    Article  CAS  PubMed  Google Scholar 

  29. Janssens KA, Oldehinkel AJ, Verhulst FC, Hunfeld JA, Ormel J, Rosmalen JG. Symptom-specific associations between low cortisol responses and functional somatic symptoms: the TRAILS study. Psychoneuroendocrinology. 2012;37:332–40. https://doi.org/10.1016/j.psyneuen.2011.06.016.

    Article  CAS  PubMed  Google Scholar 

  30. Lariviere W, Melzack M. The role of corticotropin-releasing factor in pain and analgesia. Pain. 2000;84(1):1–12.

    Article  CAS  Google Scholar 

  31. Lattin CR, Bauer CM, de Bruijn R, Michael Romero L. Hypotha-lamus–pituitary–adrenal axis activity and the subsequent response to chronicstress differ depending upon life history stage. Gen Comp Endocrinol. 2012;178:494–501.

    Article  CAS  Google Scholar 

  32. Leboeuf-Yde C, Lemeunier N, Wedderkopp N, Kjaer P. Evidence-based classification of low back pain in the general population: one-year data collected with SMS track. Chiropractic Manual Ther. 2013;21:30. https://doi.org/10.1186/2045-709x-21-30.

    Article  Google Scholar 

  33. Lemeunier N, Leboeuf-Yde C, Gagey O. The natural course of low back pain: a systematic critical literature review. Chiropractic Manual Ther. 2012;20(1):33.

    Article  Google Scholar 

  34. McBeth J, Chiu YH, Silman AJ, Ray D, Morriss R, Dickens C, Gupta A, Macfarlane GJ. Hypothalamic–pituitary–adrenal stress axis function and the relationship with chronic widespread pain and its antecedents. Arthrit Res Ther. 2005;7(5):R992–1000. https://doi.org/10.1186/ar1772.

    Article  CAS  Google Scholar 

  35. McBeth J, Silman AJ, Gupta A, Chiu YH, Ray D, Morriss R, Dickens C, King Y, Macfarlane GJ. Moderation of psychosocial risk factors through dysfunction of the hypothalamic–pituitary–adrenal stress axis in the onset of chronic widespread musculoskeletal pain: findings of a population-based prospective cohort study. Arthritis Rheum. 2007;56:360–71. https://doi.org/10.1002/art.22336.

    Article  CAS  PubMed  Google Scholar 

  36. Macedo LG, Maher CG, Latimer J, McAuley JH, Hodges PW, Rogers WT. Nature and determinants of the course of chronic low back pain over a 12-month period: a cluster analysis. Phys Ther. 2014;94(2):210–21.

    Article  Google Scholar 

  37. Osterholm MT, Kelley NS, Sommer A, Belongia EA. Efficacy and effectiveness of influenza vaccines: a systematic review and meta-analysis. Lancet Infect Dis. 2012;12(1):36–44.

    Article  Google Scholar 

  38. Paananen M, O’Sullivan P, Straker L, Beales D, Coenen P, Karppinen J, Pennell C, Smith A. A low cortisol response to stress is associated with musculoskeletal pain combined with increased pain sensitivity in young adults: a longitudinal cohort study. Arthritis Res Ther. 2015;17:355. https://doi.org/10.1186/s13075-015-0875-z.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Papadopoulos A, Cleare A. Hypothalamic–pituitary–adrenal axis dysfunction in chronic fatigue syndrome. Nat Rev Endocrinol. 2012;8:22–32.

    Article  CAS  Google Scholar 

  40. Prince K, Campbell J, Picton P, Turner S. A computational model of acute pain. Int J Simult Syst Sci Technol. 2005;6:1.

    Google Scholar 

  41. Pruessner JC, Hellhammer DH, Kirschbaum C. Burnout, perceived stress, and cortisol responses to awakening. Psychosom Med. 1999;61(2):197–204.

    Article  CAS  Google Scholar 

  42. Rockall AG, et al. CT and MR imaging of the adrenal glands in ACTH-independent cushing syndrome. Radiographics. 2004;24(2):435–52.

    Article  Google Scholar 

  43. Savic D, Jelic S, Buric N. Stability of a general delay differential model of the hypothalamo–pituitary–adrenocortical system. Int J Bifurcati Chaos. 2006;16(10):3079–85. https://doi.org/10.1142/s0218127406016665.

    Article  Google Scholar 

  44. Singh A, Kaushik A, Kumar R, Nair M, Bhansali S. Electrochemical sensing of cortisol: a recent update. Appl Biochem Biotechnol. 2014;174(3):1115–26. https://doi.org/10.1007/s12010-014-0894-2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Smith SM, Vale W. The role of the hypothalamic–pituitary–adrenal axis in neuroendocrine responses to stress. Dialog Clin Neurosci. 2006;8(4):383–95.

    Google Scholar 

  46. Smyth AR, Geraldi NR, Piehler MF. Oyster-mediated benthic-pelagic coupling modifies nitrogen pools and processes. MEPS. 2013;493:23–30.

    Article  CAS  Google Scholar 

  47. Sriram K, Rodriguez-Fernandez M, Doyle FJ III. Modeling cortisol dynamics in the neuro-endocrine axis distinguishes normal, depression, and post-traumatic stress disorder (PTSD) in humans. PLoS Comput Biol. 2012;8(2):e1002379. https://doi.org/10.1371/journal.pcbi.1002379.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Stephens MAC, Wand G. Stress and the HPA axis: role of glucocorticoids in alcohol dependence. Alcohol Res. 2012;34(4):468–83.

    PubMed  PubMed Central  Google Scholar 

  49. Stephenson ER, Kojouharov HV. A mathematical model of skeletal muscle regeneration. J Math Methods Appl Sci Biomath. 2018. https://doi.org/10.1002/mma.4908 (in Press).

    Article  Google Scholar 

  50. Strittmatter M, Bianchi O, Ostertag D, Grauer M, Paulus C, Fischer C, et al. Altered function of the hypothalamic–pituitary–adrenal axis in patients with acute, chronic and episodic pain. Schmerz. 2005;19:109–16 (Article in German).

    Article  CAS  Google Scholar 

  51. Tamcan O, Mannion AF, Eisenring C, Horisberger B, Elfering A, Muller U. The course of chronic and recurrent low back pain in the general population. Pain. 2010;150(3):451–7.

    Article  Google Scholar 

  52. Tennant F. The physiologic effects of pain on the endocrine system. Pain Ther. 2013;2(2):75–86. https://doi.org/10.1007/s40122-013-0015-x.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Tennant F. Adrenocorticotropin (ACTH) in chronic pain. J Appl Biobehav Res. 2017;22:e12063.

    Article  Google Scholar 

  54. Testa RM, et al. The usefulness of combined biochemical tests in the diagnosis of Cushing’s disease with negative pituitary magnetic resonance imaging. Eur J Endocrinol. 2007;156(2):241–8.

    Article  CAS  Google Scholar 

  55. Tomas C, Newton J, Watson S. A review of hypothalamic-pituitary-adrenal axis function in chronic fatigue syndrome. ISRN Neurosci. 2013;2013:784520.

    Article  Google Scholar 

  56. Venugopal M, Arya SK, Chornokur G, Bhansali S. A realtime and continuous assessment of cortisol in ISF using electrochemical impedance spectroscopy. Sens Actuat A Phys. 2011;172(1):154–60. https://doi.org/10.1016/j.sna.2011.04.028.

    Article  CAS  Google Scholar 

  57. Walton D, MacDermid J, Russell J, Koren G, Van Uum S. Hair-normalized cortisol waking response as a novel biomarker of hypothalamic–pituitary–adrenal axis activity following acute trauma: a proof-of-concept study with pilot results. Pain Res Treatment. 2013;2013:9. https://doi.org/10.1155/2013/876871.

    Article  CAS  Google Scholar 

  58. Wingenfeld K, Heim C, Schmidt I, Wagner D, Meinlschmidt G, Hellhammer DH. HPA axis reactivity and lymphocyte glucocorticoid sensitivity in fibromyalgia syndrome and chronic pelvic pain. Psychosom Med. 2008;70:65–72. https://doi.org/10.1097/PSY.0b013e31815ff3ce.

    Article  CAS  PubMed  Google Scholar 

  59. Xiao P, Duan L, Su J. The synaptic plasticity variability in a post-traumatic stress disorder model. In: Wang R, Pan X, editors. Advances in cognitive neurodynamics (V). Advances in cognitive neurodynamics. Singapore: Springer; 2016. p. 603–10. https://doi.org/10.1007/978-981-10-0207-6_81.

    Chapter  Google Scholar 

  60. Yehuda R, Teicher M, Trestman R, Levengood R, Siever L. Cortisol regulation in posttraumatic stress disorder and major depression: a chronobiological analysis. Biol Psychiatry. 1996;40:79–88.

    Article  CAS  Google Scholar 

  61. Yanovski JA, et al. Corticotropin-releasing hormone stimulation following low-dose dexamethasone administration. A new test to distinguish Cushing’s syndrome from pseudo-Cushing’s states. JAMA. 1993;269(17):2232–8.

    Article  CAS  Google Scholar 

  62. Yeo JKT, Babic N, Hannoush ZC, et al. Endocrine testing protocols: hypothalamic pituitary adrenal axis. In: Feingold KR, Anawalt B, Boyce A, et al., editors. Endotext. South Dartmouth: MDText.com, Inc.; 2000.

    Google Scholar 

  63. Yoshihara T, Shigeta K, Hasegawa H, Ishitani N, Masumoto Y, Yamasaki Y. Neuroendocrine responses to psychological stress in patients with myofascial pain. J Orofacial Pain. 2005;19:202–8.

    PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

Ying Du is supported by the National Natural Science Foundation of China (No.11672107).

Author information

Authors and Affiliations

Authors

Contributions

GR designed the overall project. GR, KB, JL and MB analyzed and interpreted the patient data regarding the low back pain. JS, YD and PX discussed and designed the mathematical model, and YD implemented the simulation study of the trajectories, JS was a major contributor in writing the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Robert J. Gatchel.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Su, J., Du, Y., Bevers, K. et al. Transitioning from acute to chronic pain: a simulation study of trajectories of low back pain. J Transl Med 17, 306 (2019). https://doi.org/10.1186/s12967-019-2030-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-019-2030-0

Keywords