Abstract
The allostatic theory of drug abuse describes the brain’s reward system alterations as substance misuse progresses. Neural adaptations arising from the reward system itself and from the antireward system provide the subject with functional stability, while affecting the person’s mood. We propose a computational hypothesis describing how a virtual subject’s drug consumption, cognitive substrate, and mood interface with reward and antireward systems. Reward system adaptations are assumed interrelated with the ongoing neural activity defining behavior toward drug intake, including activity in the nucleus accumbens, ventral tegmental area, and prefrontal cortex (PFC). Antireward system adaptations are assumed to mutually connect with higher-order cognitive processes occurring within PFC, orbitofrontal cortex, and anterior cingulate cortex. The subject’s mood estimation is a provisional function of reward components. The presented knowledge repository model incorporates pharmacokinetic, pharmacodynamic, neuropsychological, cognitive, and behavioral components. Patterns of tobacco smoking exemplify the framework’s predictive properties: escalation of cigarette consumption, conventional treatments similar to nicotine patches, and alternative medical practices comparable to meditation. The primary outcomes include an estimate of the virtual subject’s mood and the daily account of drug intakes. The main limitation of this study resides in the 21 time-dependent processes which partially describe the complex phenomena of drug addiction and involve a large number of parameters which may underconstrain the framework. Our model predicts that reward system adaptations account for mood stabilization, whereas antireward system adaptations delineate mood improvement and reduction in drug consumption. This investigation provides formal arguments encouraging current rehabilitation therapies to include meditation-like practices along with pharmaceutical drugs and behavioral counseling.
Introduction
The principle of allostasis was established to enhance the homeostatic model whereby the well-balanced functional state of a living being is sustained by the constant conservation of the organism’s inner environment. Each divergence from the normal state of the organism is counterbalanced by negative feedback mechanisms which support the reinstatement of original setpoints. Instead, the allostatic model advances that the internal state of the organism continuously adapts to the surrounding natural world, attaining functional stability through the adaptation of physiological thresholds (). As defined by Sterling and Eyer, “allostasis provides for continuous re-evaluation of need and for continuous readjustment of all parameters toward new setpoints” (). In humans, this continuous adaptation to the environment is reached by means of neural and endocrine processes that are able to take priority over homeostatic regulations (). The present paper is concerned with understanding and modeling the role of allostasis in drug addiction.
According to the hypothesis put forward by Koob and colleagues, a drug addict attains the allostatic state through the chronic deviation of their hedonic baseline. The addict’s physiological state is maintained operative by means of this affective adaptation, rather than by reinstatement of the original homeostatic balance. Symptoms of allostasis are manifested as changes in the addict’s mood or state of mind (). The concept of allostasis accounts for the transitions in the gradual development of a dependency whereby drugs are first experienced to “feel high” but subsequently to “feel normal.” The allostatic framework of addiction relies on changes observed in the subject’s nervous and endocrinal systems which occur as addiction perpetuates, causing continuous and progressive distortions of the subject’s affective state (). The raison d’être of these distortions is to guarantee the functional stability of the organism while its hedonic homeostatic state is corrupted.
The hedonic effect of an addictive substance on the brain’s reward system is orchestrated by within-system neuroadaptations and between-system neuroadaptations (). Experimental observations show that rats undergo a continuous degradation of hedonic valence during extended periods of cocaine consumption (). Similarly, the negative hedonic valence of the individual is increased by within-system and between-system neuroadaptations, and chronically impacts the person’s mood (). Initial drug consumption disrupts the normal synaptic physiology of the reward system () which aims to reinstate its equilibrium by means of within-system neuroadaptations (). Within-system adaptations act at the molecular or cellular levels that define the brain reward circuitry () and increase the magnitude of the ideal threshold of the reward (). For instance, if the effect of the consumed drug of abuse relies on the availability of a particular neurotransmitter, within-system adaptations will diminish its amount within the reward system (). With repeated drug administration the reward system becomes accustomed to extended activations of the within-system component, which eventually causes withdrawal symptoms during periods of abstinence (). Dopamine in the nucleus accumbens (NAc) and extended amygdala plays an important role in within-system adaptations (). As consumption further advances, the brain’s expectation for future rewards increases and within-system neuroadaptations become progressively inadequate and eventually fail to provide the individual with a well-adjusted functional state.
Due to deficiencies of within-system adaptations, brain structures different than the ones defining the reward system are recruited through the deployment of between-system neuroadaptations to further counterbalance the effect of the drug (). These brain structures, delineated by Koob and Le Moal, embody the “antireward systems” (). Between-system adaptations increase the baseline reward threshold () and originate in the brain stress system (). For example, if a drug has a particular effect on the reward circuit, between-system activations could promote a hormonal response leading to the opposite effect (). The corticotropin-releasing factor operating in the amygdala, stria terminalis, and ventral tegmental area (VTA) plays an important role in between-system adaptations ().
The current investigation was undertaken to further explore the allostatic theory of addiction and combines two existing computational models: (i) a pharmacological model which describes how decreased reward function induces compulsive cocaine self-administration in rats () and (ii) a dynamical system model that represents hypotheses in terms of behavior, cognition, and neuropsychology, occurrences of spontaneous remissions from drug use in humans (). More specifically, these models are connected by two computational hypotheses and one provisional assumption. The first hypothesis outlines within-system neuroadaptations, the second hypothesis delineates between-system neuroadaptations, and the provisional assumption relates to the mood of a virtual subject. In particular, within-system adaptations reflect changes in the reward system, whereas between-system adaptations embody modification in the antireward system. Mood estimation relies on the direct effect of the drug and its influence on these two neuroadaptations.
Methods: Computational Framework for Allostasis
Figure 1 represents the model at the center of this investigation. The color code denotes different scales of observation: behavior (blue), neuropsychology (green), cognition (red), healing (orange), and pharmacology (light blue). The outputs of this discrete-time model reside in the pharmacological scale which includes computational predictions of drug intakes, Z(t), and mood, M(t). The processes composing the model are estimated over time-scales of minutes (t) and hours (t*). With the exception of Z, which is a binary variable, all other variables have continuous values with units of measurements that are left unspecified since corresponding human measurements are not available at this time. The current section describes the two hypotheses, the provisional assumption, and the dynamical system model at the center of the present investigation. Formal definitions of the processes in Figure 1 are included in the Supplementary Material.
Figure 1
The first hypothesis identifies reward system adaptations with ongoing neural activities within NAc, VTA, and prefrontal cortex (PFC); the second hypothesis associates antireward system adaptations with higher-order cognitive processes within PFC, orbitofrontal cortex (OFC), and anterior cingulate cortex (ACC); and the provisional assumption considers the virtual subject’s mood as a combination of reward components. Hypotheses and the provisional assumption reside in the pharmacological scale of the framework which extends the pharmacokinetic/pharmacodynamic (PK/PD) model of allostasis for laboratory rats of Ahmed and Koob (
Within-system alterations are represented by changes in the drug potency index, T50, and between-system adjustments are described by variations in the baseline reward threshold, T0. Ahmed and Koob (
Hypothesis 1: Within-system neuroadaptations (pharmacological scale)
Ahmed and Koob (
Along those lines, our model is based on the following hypothesis: within-system adaptations are expressed by the reward set point, TS, which deteriorates for maladaptive behaviors and improves for healthy behaviors. More specifically, within-system neuroadaptations depend upon the virtual subject’s current attitude toward drug use, G(t*), which emulates the translation of ongoing neural activities within the NAc, VTA, and PFC into behavior. Computationally, after the first drug intake, TS is assumed to monotonically increase for healthy behavior and exponentially decrease for maladaptive behavior: where λ, β, and γ, are constants; tc corresponds to the last time t* when G changed its sign; d is a temporal unit-step counter reset to 0 when G changes its sign, ensuring mathematical continuity to TS; and ΣZ ≥ 1 denotes that at least one drug intake occurred up to time t*.
Hypothesis 2: Between-system neuroadaptations (pharmacological scale)
Goldstein and Volkow (
Along those lines, our model considers the following hypothesis: between-system adaptations are expressed by the baseline reward threshold, T0, which increases when maladaptive behavior occurs and decreases during healthy behavior. After some drug intakes, cognitive weights which adapt to maladaptive behavior will cause T0 to increase. In case the virtual subject regains a healthy behavior, these cognitive weights cause T0 to decrease. More specifically, between-system alterations are defined as a function of cognitive time-dependent weights, ωS(t*), ωP(t*), and ωD(t*), and healing intervention, H(t*). Cognitive weights mimic associative learning between the drug and its pleasurable effect, whereas H emulates positive cognitive changes toward remission from drug intake. Computationally, after a number of drug intakes, the cognitive substrate of the virtual subject starts to stochastically influence T0 by means of ωS, ωP, ωD and the healing intervention H. The modulating influence of ωS, ωP, and ωD on T0 has opposite valence when a healing intervention occurs (H = 1) than when it does not (H = 0): where δT0 and α are constants, and Σ ≥ α denotes the period subsequent to at least α drug intakes.
Provisional assumption: Mood (behavioral/pharmacological scale)
A classical validation schema, where simulations are compared to laboratory data is not suitable for the proposed model, as human measures of the considered processes are not currently available. Instead, for qualitative evaluation purposes, an additional provisional assumption is formulated to predict the subject’s mood, M, as an aggregate of neural and psychological components. The former component relies on the direct rush/comedown effect of the drug, rc(t), and the latter upon cognitive distortions, cd(t), emulated as a function of current and previous hedonic adaptations.
The evolution of rc is defined as the summation of piece-wise sinusoidal functions each of one period with slightly exponentially decaying tails that initiate when a drug intake occurs. Cognitive distortions related to addiction are assumed to depend on the overall current reward state of the virtual subject which includes T, TS, and T0. Healthy individuals should not suffer from cognitive distortions: no current drug’s effect on reward threshold (T) should arise, nor should any negative hedonic valences from within- and between-system neuroadaptations (TS and T0). The speculative cd combines the current reward state of the individual, the current activation of within- and between-system neuroadaptations, alongside of previously experienced hedonic adaptations:
where γM is constant; ΣZ ≥ 1 denotes that at least one drug intake occurred up to t; and ΔTSO stands for the arithmetic difference between TS and T0. The formulation of cd is suggested by the temporal difference component employed in the first model of learning mechanism associated with dopaminergic neurons in the basal ganglia (
Behavioral, cognitive, and neuropsychological scales of the model
The tendency for drug-seeking behavior, G(t*), defines the behavioral scale of the model (in blue on Figure 1). Healthy behavior, i.e., avoidance of drug use, corresponds to positive values of G, and maladaptive behavior, i.e., a tendency toward drug-seeking behavior, corresponds to negative values. G assesses how neural activity of brain regions sensitive to addictive drugs affects the ratio of compulsion and inhibition (
The cognitive apparatus of the virtual subject (in red on Figure 1) transforms the ongoing neuropsychological activities into information accessible at the behavioral scale. The ratio of compulsion and inhibition is driven by the cognitive state, cs(t*), which is a mathematical transformation into the real interval [0, 1] of the virtual subject’s rationality density, rd(t*). A more compulsive behavior is expressed when cs is equal to 0, whereas a stronger inhibitory behavior is expressed when cs is equal to 1. In terms of economic theories of addiction, cs and rd delineate whether the virtual subject is irrational, imperfectly rational, or rational (
The neuropsychological activities (in green on Figure 1) encompass internal processes (S, P, D, and Q) and external triggers (AS, AP, AD, and AQ). The negative affective state of nervousness, anxiety, or stress, S(t*), of an addict expands during withdrawal phases (
Severe stressors, AS(t*), such as electric foot-shocks for laboratory rats (
Healing interventions (in orange on Figure 1), H(t*), can cause immediate abstinence by direct alteration of the virtual subject’s rationality rd and by adapting cognitive weights of internal and external processes. These modifications can endure over time and reflect real-life occurrences of “maturing out” of addiction (
As discussed above, the pharmacological scale (in light blue on Figure 1) includes the binary decision toward Z which depends on the arithmetic difference between the lowering effect on reward threshold, T(t), and the reward set point, TS(t*), similarly to the model in Ref. (
Illustrative Results
The experimental set up includes mathematical definitions of the processes described above and detailed in the Supplementary Material, as well as their correspondent implementation in MATLAB®. The initial conditions of the presented simulations are defined by 71 parameters whose values, reported in Table S1 in Supplementary Material, are chosen according to Ref. (
Three Case Studies based on the same experimental set up, and narrating tobacco smoking, are presented herein to exemplify the methods described in the previous section: (i) escalation from early to heavy smoking, (ii) conventional therapeutic process (e.g., nicotine patches), and (iii) alternative medical treatment (e.g., meditation). Case Study 1 stands for the evaluation baseline; Case Study 2 is identical but includes two long-lasting healing interventions, H; whereas Case Study 3 encompasses eight short-lasting healing interventions. The total amount of time for which H is active in these Case Studies is respectively 0, 10, and 5 days.
Each Case Study includes three Evaluations: in the first both TS and T0 are time-dependent processes according to Eqs 1 and 2; in the second TS is constant and T0 is time-dependent; and in the third TS is time-dependent and T0 is constant. The simulation results include means for 100 runs and 95% simulation envelopes corresponding to a period of 160 days with the drug becoming available on the fifth day. Changes in the allostatic state of the virtual subject are estimated through variations of the mood M. The daily abstinence index is the ratio (as a percentage) between the runs in which the virtual subject don’t use the addictive drug and the total number of simulation runs. This index assesses the virtual subject’s stage according to the classification by DiFranza et al. (
Case study 1: Allostatic state trajectory during escalation of drug consumption
Figures 2 and 3 depict the behavior of a virtual subject who engages in cigarette-smoking 5 days after the simulation begins. In Evaluation 1, changes in ωS, ωP, and ωD are gradual; T is at first weaker than TS but eventually surpasses it, and T0 continually increases; M increasingly oscillates around its downslope; the average drug consumption increases and the abstinence index diminishes. Upon completion of the simulation, this virtual subject is characterized by an average consumption of ~42 intakes/day (~2 packs) and an abstinence index of ~11%, comparable to a severe stage 4 (
Figure 2

Case Study 1 (cognitive weights and reward components). Simulations of virtual behavior for cigarette consumption over a period of 160 days. Cigarettes are available on the fifth day. Results are the average of 100 runs. Evaluation 1 simulates both TS and T0 as time-dependent processes. TS is constant and T0 time-dependent in Evaluation 2, TS is time-dependent and T0 constant in Evaluation 3. (A1–A3) show the evolution of cognitive weights ωS, ωP, and ωD; and (B1–B3) the progression of T, T0, and TS. The shades correspond to 95% simulation envelopes. The time-scales are hours for (A1–A3), and minutes for (B1–B3). Further details of these simulations are reported in Figures S1 and S2 in Supplementary Material.
Figure 3

Case Study 1 (mood and health state assessments). Simulations of virtual behavior for cigarette consumption over a period of 160 days. Cigarettes are available on the fifth day. Results are the average of 100 runs. Evaluation 1 simulates both TS and T0 as time-dependent processes. TS is constant and T0 time-dependent in Evaluation 2, TS is time-dependent and T0 constant in Evaluation 3. (A) Shows the evolution of the subject’s mood M; (B) the average number of drug intakes; and (C) the abstinence index. The shades correspond to 95% simulation envelopes. The time-scales are minutes for (A), and days for (B,C). Further details of these simulations are reported in Figures S1 and S2 in Supplementary Material.
In Evaluation 2, the cognitive adaptations occur significantly faster than in Evaluation 1. T and TS are approximately equal at first, but eventually T becomes larger than TS as T0 constantly increases. M abruptly decreases during the loading phase and subsequently manifests a negative trend enclosed by minor oscillations. The average consumption quickly increases and the virtual subject reaches a satiety state of ~48 cigarettes/day. Note that the number of intakes defining satiety is not an explicit constraint defined in the model. The abstinence index is consistently at zero.
In Evaluation 3, the evolution of the virtual subject’s processes is very similar to the predictions for Evaluation 1, but T decreases and M has a less negative downslope. This evaluation ends with ~34 intakes/day and ~25% abstinence.
Additional details can be found in the Supplementary Material. Figures S1 and S2 in Supplementary Material include the fine details for Case Study 1, and Figures S3–6 in Supplementary Material show how different probabilities defining changes in ωS, ωP, and ωD influence the predicted consumption rates and mood downslope. If the first smoked cigarette within the simulations is considered as the first ever in the life of the virtual subject, then Figures S3–6 in Supplementary Material can be considered to relate to different rates of progression from recreational smoking toward heavy smoking.
Case study 1 indicates that a virtual subject consuming drugs for the first time, or relapsing after a period of abstinence, undergoes a continuous negative shift in mood baseline which directly correlates with the strength of cognitive learning facilitating drug consumption. In addition, the virtual subject suffers growing mood swings during protracted consumption. When TS is constant, the mood substantially decreases during the first number of drug intakes, and the virtual subject rapidly reaches the satiety consumption rate. With T0 constant, the simulated mood has a weaker negative tendency and oscillates less.
This Case Study shows that escalation in drug consumption occurs together with chronic depression of mood. When just between-system adaptations are operative, these simulations predict that the virtual subject’s mood strongly drops, whereas it moderately decreases when just within-system adaptations are operative.
Case study 2: Allostatic state trajectory during conventional therapies
The profile presented in Figures 4 and 5 is similar to Case Study 1 but also includes healing interventions (H). Five-day long H events are activated at t = 1920 and t = 2280 [hours]. This is intended to emulate 25 days of a replacement therapy using a nicotine transdermal system for two 5 days periods separated from each other by 15 days. All evaluations for Case Studies 1 and 2 are similar until H is activated.
Figure 4

Case Study 2 (cognitive weights and reward components). Simulations of virtual behavior for cigarette consumption over a period of 160 days. Cigarettes are available on the fifth day. Results are the average of 100 runs. Evaluation 1 simulates both TS and T0 as time-dependent processes. TS is constant and T0 time-dependent in Evaluation 2, TS is time-dependent and T0 constant in Evaluation 3. In all evaluations, the recovery process H is activated at t = 1920 and t = 2280 [hours] (in light pink) and stays active for 120 h (dark pink). (A1–A3) show the evolution of cognitive weights ωS, ωP, and ωD; and (B1–B3) the progression of T, T0, and TS. The shades correspond to 95% simulation envelopes. The time-scales are hours for (A1–A3), minutes for (B1–B3). Further details of these simulations are reported in Figures S7 and S8 in Supplementary Material.
Figure 5

Case Study 2 (mood and health state assessments). Simulations of virtual behavior for cigarette consumption over a period of 160 days. Cigarettes are available on the fifth day. Results are the average of 100 runs. Evaluation 1 simulates both TS and T0 as time-dependent processes. TS is constant and T0 time-dependent in Evaluation 2, TS is time-dependent and T0 constant in Evaluation 3. In all evaluations, the recovery process H is activated at t = 1920 and t = 2280 [hours] (in light pink) and stays active for 120 h (dark pink). (A) shows the evolution of the subject’s mood M; (B) the average number of drug intakes; and (C) the abstinence index. The shades correspond to 95% simulation envelopes. The time-scales are minutes for (A), and days for (B,C). Further details of these simulations are reported in Figures S7 and S8 in Supplementary Material.
During the first therapeutic period in Figures 4 and 5, ωS, ωP, and ωD are influenced by H to promote healthier behavior. This positive effect partially persists after the first period of therapy and is further strengthened by the second. In Evaluation 1, the activation of H causes a small upswing in T, a strong upswing in TS, and a decrease in the upslope of T0. The degradation of M becomes less accentuated after the treatment. The average drug consumption drops and the abstinence index rises while H is active. This experiment’s endpoint is comparable to an advanced stage 3 or an intermediate stage 4 (
Additional fine detail for Case Study 2 can be found in Figures S7 and S8 in Supplementary Material. Figures S9–12 in Supplementary Material show how different probabilities defining the influence of H affect the permanent predicted consumption rates. Higher probabilities lead the virtual subject to stage 1 or intermediate stage 2, whereas lower probabilities to advanced stage 4. The same sets of probabilities are tested when TS is constant (Figures S13–16 in Supplementary Material) and when T0 is constant (Figures S17–20 in Supplementary Material). For TS constant, the virtual subject always gets to a satiety consumption rate and reveals a slight shy positive trend in M for the highest probabilities along with a decrease in average consumption. For T0 constant, M becomes constant after the therapy and its variations become smaller as the tested probabilities become higher.
Case study 2 shows how a few, though long, healing interventions diminish the negative trend in the virtual addict’s mood. Early indications of mood increase appear when healing signals are highly effective. When TS is constant, curative effects on the mood are negligible, unless cognitive learning is exceptionally successful. Even though positive, the influences on mood for this extreme case are quite limited. With T0 constant, the mood stabilizes after healing interventions and remains roughly constant.
This Case Study exemplifies the effect of prolonged healing interventions. When just between-system adaptations are operative, these simulations predict that the virtual subject’s mood continues to worsen in spite of therapeutic events. When just within-system adaptations are operative, the predicted mood stabilizes as a consequence of healing periods.
Case study 3: Allostatic state trajectory during alternative medical treatments
The profile in Figures 6 and 7 present a different type of healing intervention than in Case Study 2. At each of the simulated times t ∈ {1920, 1960, 2000, 2040, 2280, 2320, 2360, 2400} [hours], a 15 h H event is activated. This is intended to emulate two 5 day healing periods during which the virtual subject undergoes four meditation practices. The benefit of each practice lasts for 15 h, and the two healing periods are 10 days apart. Other than H event durations and activation times, all the parameters defining this case study are the same as in Case Studies 1 and 2.
Figure 6

Case Study 3 (cognitive weights and reward components). Simulations of virtual behavior for cigarette consumption over a period of 160 days. Cigarettes are available on the fifth day. Results are the average of 100 runs. Evaluation 1 simulates both TS and T0 as time-dependent processes. TS is constant and T0 time-dependent in Evaluation 2, TS is time-dependent and T0 constant in Evaluation 3. In all evaluations, the recovery process H is activated at t ∈ {1920, 1960, 2000, 2040, 2280, 2320, 2360, 2400} [hours] (in light pink) and stays active for 15 h (dark pink). (A1–A3) show the evolution of cognitive weights ωS, ωP, and ωD; and (B1–B3) the progression of T, T0, and TS. The shades correspond to 95% simulation envelopes. The time-scales are hours for (A1–A3), and minutes for (B1–B3). Further details of these simulations are reported in Figures S21 and S22 in Supplementary Material.
Figure 7

Case Study 3 (mood and health state assessments). Simulations of virtual behavior for cigarette consumption over a period of 160 days. Cigarettes are available on the fifth day. Results are the average of 100 runs. Evaluation 1 simulates both TS and T0 as time-dependent processes. TS is constant and T0 time-dependent in Evaluation 2, TS is time-dependent and T0 constant in Evaluation 3. In all evaluations, the recovery process H is activated at t ∈ {1920, 1960, 2000, 2040, 2280, 2320, 2360, 2400} [hours] (in light pink) and stays active for 15 h (dark pink). (A) shows the evolution of the subject’s mood M; (B) the average number of drug intakes; and (C) the abstinence index. The gray shades correspond to 95% simulation envelopes. The time-scales are minutes for (A), and days for (B,C). Further details of these simulations are reported in Figures S21 and S22 in Supplementary Material.
During the first period of meditation, ωS, ωP, and ωD strongly adjust. The enduring positive changes are additionally expanded after the second period of meditation. In Evaluation 1, T decreases before meditation, increases during meditation, and finally decreases again afterward. There is an upswing of TS when H is active. T0 initially increases, becomes constant after the first healing session, and decreases after the second. After the first period of meditation, M stops declining, and after it increases with reduced oscillation. The average drug consumption significantly decreases after the first healing period and further decreases after the second. The abstinence index robustly increases during the treatment. The endpoint of this evaluation corresponds to ~1 intake/day and ~97% abstinence, placing the virtual subject in stage 1 or early stage 2 (
Additional details for Case Study 3 can be found in the Supplementary Material. Figures S21 and S22 in Supplementary Material show the fine details. Figures S23–26 in Supplementary Material illustrate how different probabilities defining the influence of H permanently impact the predicted consumption rates. Higher probabilities lead the virtual subject to cease using the drug, whereas lower probabilities lead the subject to advanced stage 1 or intermediate stage 2. The same sets of probabilities are tested when TS is constant (Figures S27–30 in Supplementary Material) and when T0 is constant (Figures S31–34 in Supplementary Material). For TS constant, the virtual subject always gets to its satiety consumption rate, and after the treatment M increases. For T0 constant, M tends to become constant after the therapy and its variations become smaller as the probability becomes higher.
Case study 3 displays how phasic healing interventions increase the mood of the virtual addict. This increase becomes very conspicuous while healing signals have high effectiveness. Also, when TS is constant, mood increases as a direct correlation of healing success. With T0 constant, after healing interventions the mood stabilizes and becomes approximately constant.
This case study exemplifies the effect of brief healing interventions that follow one another at short intervals. When just between-system adaptations are operative, these simulations predict that the virtual subject’s mood increases as a result of curative events. When just within-system adaptations are operative, the predicted mood stabilizes but does not improve after healing periods.
Discussion
The computational framework for allostasis presented above unites and elaborates two earlier formal models (
The present article is an exploratory instance of knowledge repository (KR) modeling for addiction (
The computational description of a complex phenomenon such as drug use and abuse requires finding a compromise between two desirable but incompatible objectives. On the one hand, the biological components defining the model embrace a simplified ontology of addiction, and on the other hand, the mathematical features of the framework include a sizable number of elements and parameters making the model underconstrained. Moreover, a useful formal system should suggest testable hypotheses to further advance the investigated field. These perspectives are considered herein.
Biological conjectures and limitations
The simulations discussed in this investigation represent archetypal patterns of drug-seeking including transitions from recreational to heavy use, and rehabilitation. They express the comorbidity between addiction and mood disorders for an addict vulnerable to both reward and antireward system neuroadaptations. As drug intake proceeds, the model estimates a steady decrease in the addict’s mood that increasingly oscillates until the virtual subject reaches the satiety rate of drug consumption. This computational model also suggests a possible remission of the individual’s healthy state as a consequence of cognitive adjustments induced by conventional or alternative treatments. The model predicts a contraction in the addict’s negative mood tendency and fluctuations while solely reward system neuroadaptations influence the hedonic valence of the individual. This rigidity endures during healing interventions as the model predicts stabilization of the subject’s mood, rather than its increase. When the unique source of neuroadaptations affecting the brain’s reward system relies on the antireward system, the model predicts a noteworthy negative deflection of the subject’s mood during the first number of drug intakes. The model also predicts that neuroadaptations occurring during healing periods, and that are uniquely induced by the antireward system, empower the individual with the possibility to regain a healthy mood state. These simulations also suggest that the satiety rate of drug consumption is reached more rapidly when the individual expresses only antireward system adaptations.
An important biological limitation of this model resides in the omission of mechanisms responsible for the increase of pharmacodynamic tolerance, which arises when receptors or second messengers are blunted by drugs such as alcohol or opiates (
The discussed model could be further elaborated to consider alcohol dependence and treatment (
The model’s high-dimensionality
A KR model is inclined toward a high-dimensionality due to a large number of descriptive variables since its objective is to describe the studied phenomenon as comprehensively as possible. The predictions presented in this investigation rely on the high-dimensional dynamical model shown in Figure 1 that comprises 21 time-dependent biological processes translated into the same number of mathematical expressions. In addition, there are 71 parameters, the setting of which can dramatically affect the model’s behavior. Such high-dimensionality could simulate an assortment of dynamics larger than the ones expressed by living creatures, consequently limiting the model’s predictive power. The natural processes included in the presented computational framework are sizable, yet their descriptions are determined conservatively. For instance, the processes comprising the neuropsychological scale of the model have limited domains of definition which facilitate their tractability. These restricted domains reflect biological plausibility and ease the sensitivity analysis (
Another attempt at mathematical moderation resides in the definition of healing interventions. The same mathematical definition used with different calibrations to mimic conventional and alternative cures was intentionally deployed as a lower-bound estimation of real-life cleansing episodes. In fact, the minimal duration of nicotine replacement therapies ranges from 3 weeks to 3 months (
Implications for treatments
This article provides formal arguments, conditional upon the validity of the hypotheses defining the computational framework, in favor of a stronger consideration of the addict’s cognitive state evolution throughout the treatment process. In particular, this investigation suggests that higher rates of rehabilitation from drug addiction in humans can be reached by combining medical therapies that employ pharmaceutical drugs and counseling along with non-conventional treatments.
Several pharmacotherapies are available for smoking cessation, as for example nicotine in various forms (gums, patches, inhalers, tablets) and antidepressant drugs (
Behavioral counseling and personal support are instances of behavioral-cognitive therapies: non-pharmacological interventions which have flourished since the 1960s for the treatment of depression and anxiety (67–69). A complementary category of non-pharmacological interventions resides in mind-body practices, which include mindfulness meditation, guided imagery, and relaxation (70). Both the definitions of behavioral-cognitive (69) and mind-body (70) practices rely on the beneficial impact that healthy cognitive states exert on the overall person’s well being. A patient undergoing behavioral-cognitive therapy learns how to recognize and manage real-life situations that are negatively evaluated because of cognitive distortions, whereas mind-body practices provide the patient with a more realistic awareness that decreases irrational thoughts. In both cases the objective is to lead the patient to healthier physical and psychological states.
Behavioral-cognitive practices have demonstrated their positive impact on smoking cessation (66, 71), and mind-body techniques related to the treatment of nicotine addiction are restrained by a shortfall of related investigations (72), even though recent studies demonstrate their great potential. Preliminary experimental support in favor of mindfulness meditation as a practice decreasing relapse rates for post-rehabilitation patients was provided in a study including 168 participants who ceased the use of substances including alcohol, cocaine, and methamphetamines (73).
The rational speculation that arises while considering pharmacological and non-pharmacological healing practices suggests that current therapies deploying one or multiple pharmacological means along with counseling will raise their success rates by uniting with alternative medical practices. The computational framework presented in this investigation provides formal arguments to endorse this conjecture as the allostatic state of an addict, assessed through mood variations, is shown to improve because of cognitive interventions provided by practices comparable to those of conventional and alternative medicine.
If the predictions delivered by the computational model discussed in this investigation constitute a fair approximation to describe how cigarette-smoking influences the allostatic state of a human addicted to nicotine, then it is expected that an integrative medicine approach to drug rehabilitation will provide higher cessation rates and lower relapse rates than current therapies. Given that the allostatic theory of addiction is not limited to the description of a particular substance of abuse, this prediction should apply by extension also to addictive substances including heroin, alcohol, and others.
Concluding remarks
The computational model presented in this article considers drug addiction as a disease (74). Such viewpoint may lead to “increasing alienation, stigmatization, and social distance” (75) of human beings abusing drugs. Multi-leveled overlooks of addiction that include biological/psychological/social (76), and even spiritual (77) elements are suggested as potential candidates to restrain such an undesirable possibility (78). Along these lines, the multiscale standpoint of the framework shown in Figure 1 aims to promote a comprehensive understanding of addiction that provides prospect for recovery, which seems to occur more often than commonly believed (79). The present investigation develops a multiscale computational model to further explore the allostatic theory of addiction (
Supplementary Material
The Supplementary Material for this article can be found online at http://www.frontiersin.org/Journal/10.3389/fpsyt.2013.00167/abstract
Statements
Acknowledgments
The authors acknowledge valuable discussions with Serge H. Ahmed, Philip S. Thomas, Bettina Schaefli, Dimitri V. Nowicki, Jérémie Cabessa, Lorenzo Cantelli, and Bertrand Steiner. We thank the reviewers for their helpful and insightful comments. This work was supported in part by the National Science Foundation under NSF grant #CNS-0619337. Any opinions, findings conclusions or recommendations expressed here are the authors and do not necessarily reflect those of the sponsors.
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.
References
1
SterlingPEyerJ. Allostasis: a new paradigm to explain arousal pathology. In: FisherSReasonJ editors. Handbook of Life Stress Cognition and Health. Oxford: Wiley & Sons (1998). p. 629–49.
2
KoobGFLe MoalM. Drug addiction, dysregulation of reward, and allostasis. Neuropsychopharmacology (2001) 24(2):97–129.10.1016/S0893-133X(00)00195-0
3
KoobGFLe MoalM. Drug abuse: hedonic homeostatic dysregulation. Science (1997) 278(5335):52–8.10.1126/science.278.5335.52
4
KoobGFBloomFE. Cellular and molecular mechanisms of drug dependence. Science (1998) 242(4879):715–23.10.1126/science.2903550
5
AhmedSHKennyPJKoobGFMarkouA. Neurobiological evidence for hedonic allostasis associated with escalating cocaine use. Nat Neurosci (2002) 5(7):625–6.10.1038/nn872
6
CamíJFarréM. Mechanisms of disease: drug addiction. N Engl J Med (2003) 349(10):975–86.10.1056/NEJMra023160
7
KoobGFLe MoalM. Addiction and the brain antireward system. Annu Rev Psychol (2008) 59:29–53.10.1146/annurev.psych.59.103006.093548
8
AhmedSHKoobGF. Transition to drug addiction: a negative reinforcement model based on an allostatic decrease in reward function. Psychopharmacology (2005) 180(3):473–90.10.1007/s00213-005-2180-z
9
KoobGF. A role for brain stress systems in addiction. Neuron (2008) 59(1):11–34.10.1016/j.neuron.2008.06.012
10
GeorgeOLe MoalMKoobGF. Allostasis and addiction: role of the dopamine and corticotropin-releasing factor systems. Physiol Behav (2012) 106(1):58–64.10.1016/j.physbeh.2011.11.004
11
KoobGFLe MoalM. Plasticity of reward neurocircuitry and the ‘dark side’ of drug addiction. Nat Neurosci (2005) 8(11):1442–4.10.1038/nn1105-1442
12
LevyYZLevyDMeyerJS. Computational hypothesis for maturing out of addiction and mindfulness-based cognitive techniques. In: SamsonovichAVJóhannsdóttirKR editors. Biologically Inspired Cognitive Architectures 2011 – Proceedings of the Second Annual Meeting of the BICA Society. Amsterdam: IOS Press (2011). p. 222–31.
13
WiseRA. Addictive drugs and brain stimulation reward. Annu Rev Neurosci (1996) 19:319–40.10.1146/annurev.ne.19.030196.001535
14
AhmedSHKoobGF. Transition from moderate to excessive drug intake: change in hedonic set point. Science (1998) 282(5387):298–300.10.1126/science.282.5387.298
15
GoldsteinRZVolkowND. Drug addiction and its underlying neurobiological basis: neuroimaging evidence for the involvement of the frontal cortex. Am J Psychiatry (2002) 159(10):1642–52.10.1176/appi.ajp.159.10.1642
16
YuricaCLDiTomassoRA. Cognitive distortions. In: FreemanAFelgoiseSHNezuCMNezuAMReineckeMA editors. Encyclopedia of Cognitive Behavior Therapy. New York: Springer (2005). p. 117–22.
17
MorgensternJMcCradyBS. Curative factors in alcohol and drug treatment: behavioral and disease model perspectives. Br J Addict (1992) 87(6):901–12.10.1111/j.1360-0443.1992.tb01985.x
18
ChaseHWClarkL. Gambling severity predicts midbrain response to near-miss outcomes. J Neurosci (2010) 30(18):6180–7.10.1523/JNEUROSCI.5758-09.2010
19
FortuneEEGoodieAS. Cognitive distortions as a component and treatment focus of pathological gambling: a review. Psychol Addict Behav (2012) 26(2):298–310.10.1037/a0026422
20
PowersPSPerezABoydFRosemurgyA. Eating pathology before and after bariatric surgery: a prospective study. Int J Eat Disord (1999) 25(3):293–300.10.1002/(SICI)1098-108X(199904)25:3<293::AID-EAT7>3.0.CO;2-G
21
DavisRA. A cognitive-behavioral model of pathological Internet use. Comput Hum Behav (2001) 17(2):187–95.10.1080/00207594.2010.541256
22
BartoAG. Adaptive critics and the basal ganglia. In: HoukJCDavisJLBeiserDG editors. Models of Information Processing in the Basal Ganglia. Cambridge, MA: MIT Press (1995). p. 215–32.
23
HoukJCAdamsJLBartoAG. A model of how the basal ganglia generate and use neural signals that predict reinforcement. In: HoukJCDavisJLBeiserDG editors. Models of Information Processing in the Basal Ganglia. Cambridge, MA: MIT Press (1995). p. 249–70.
24
LevyYZLevyDMeyerJSSiegelmannHT. Drug addiction: a computational multiscale model combining neuropsychology, cognition, and behavior. In: EncarnaçãoPVelosoA editors. Proceedings of the International Conference on Bio-inspired Systems and Signal Processing. Portugal: INSTICC Press (2009). p. 87–94.
25
RobinsonTBerridgeK. The neural basis of drug craving: an incentive-sensitization theory of addiction. Brain Res Rev (1993) 18(3):247–291.10.1016/0165-0173(93)90013-P
26
SloanFAWangY. Economic theory and evidence on smoking behavior of adults. Addiction (2008) 103(11):1777–85.10.1111/j.1360-0443.2008.02329.x
27
GoldsteinRZVolkowND. Dysfunction of the prefrontal cortex in addiction: neuroimaging findings and clinical implications. Nat Rev Neurosci (2011) 12(11):652–69.10.1038/nrn3119
28
VolkowNDWangG-JFowlerJSTomasiD. Addiction circuitry in the human brain. Annu Rev Pharmacol Toxicol (2012) 52:321–36.10.1146/annurev-pharmtox-010611-134625
29
HodginsDCEl-GuebalyNArmstrongS. Prospective and retrospective reports of mood states before relapse to substance use. J Consult Clin Psych (1995) 63(3):400–7.10.1037/0022-006X.63.3.400
30
Aston-JonesGHarrisGC. Brain substrates for increased drug seeking during protracted withdrawal. Neuropharmacology (2004) 47(Suppl 1):167–79.10.1016/j.neuropharm.2004.06.020
31
De AlbaISametJHSaitzR. Burden of medical illness in drug- and alcohol-dependent persons without primary care. Am J Addict (2004) 13(1):33–45.10.1080/10550490490265307
32
KalivasPWDuffyP. Effect of acute and daily cocaine treatment on extracellular dopamine in the nucleus accumbens. Synapse (1990) 5(1):48–58.10.1002/syn.890050104
33
HeidbrederCAThompsonACShippenbergTS. Role of extracellular dopamine in the initiation and long-term expression of behavioral sensitization to cocaine. J Pharmacol Exp Ther (1996) 278(2):490–502.
34
VolkowNDWangGJFowlerJSHitzemannRAngristBGatleySJet alAssociation of methylphenidate-induced craving with changes in right striato-orbitofrontal metabolism in cocaine abusers: implications in addiction. Am J Psychiatry (1999) 156(1):19–26.
35
ErbSShahamYStewartJ. Stress reinstates cocaine-seeking behavior after prolonged extinction and a drug-free period. Psychopharmacology (1996) 128(4):408–12.10.1007/s002130050150
36
SinhaRFuseTAubinLRO’MalleySS. Psychological stress, drug-related cues and cocaine craving. Psychopharmacology (2000) 152(2):140–8.10.1007/s002130000499
37
PolliniRAMcCallLMehtaSHVlahovDStrathdeeSA. Non-fatal overdose and subsequent drug treatment among injection drug users. Drug Alcohol Depend (2006) 83(2):104–10.10.1016/j.drugalcdep.2005.10.015
38
BarthJCritchleyJBengelJ. Efficacy of psychosocial interventions for smoking cessation in patients with coronary heart disease: a systematic review and meta-analysis. Ann Behav Med (2006) 32(1):10–20.10.1207/s15324796abm3201_2
39
de WitHStewartJ. Drug reinstatement of heroin-reinforced responding in the rat. Psychopharmacology (1983) 79(1):29–31.10.1007/BF00433012
40
SpealmanRDBarrett-LarimoreRLRowlettJKPlattDMKhroyanTV. Pharmacological and environmental determinants of relapse to cocaine-seeking behavior. Pharmacol Biochem Behav (1999) 64(2):327–36.10.1016/S0091-3057(99)00049-0
41
SeeRE. Neural substrates of conditioned-cued relapse to drug-seeking behavior. Pharmacol Biochem Behav (2002) 71(3):517–29.10.1016/S0091-3057(01)00682-7
42
HymanSE. Addiction: a disease of learning and memory. Am J Psychiatry (2005) 162(8):1414–22.10.1176/appi.ajp.162.8.1414
43
RobinsonTEBerridgeKC. Addiction. Annu Rev Psychol (2003) 54:25–53.10.1146/annurev.psych.54.101601.145237
44
SteadLFPereraRBullenCMantDLancasterT. Nicotine replacement therapy for smoking cessation. Cochrane Database Syst Rev (2008) 1:CD000146.10.1002/14651858.CD000146.pub3
45
BrewerJAMallikSBabuscioTANichCJohnsonHEDeleoneCMet alMindfulness training for smoking cessation: results from a randomized controlled trial. Drug Alcohol Depend (2011) 119(1–2):72–80.10.1016/j.drugalcdep.2011.05.027
46
MuzicRFBerridgeMSFriedlandRPZhuNNelsonAD. PET quantification of specific binding of carbon-11-nicotine in human brain. J Nucl Med (1998) 39(12):2048–54.
47
DiFranzaJRWellmanRJSavageauJA. Does progression through the stages of physical addiction indicate increasing overall addiction to tobacco?Psychopharmacology (2012) 219(3):815–22.10.1007/s00213-011-2411-4
48
BobashevGCostenbaderEGutkinBS. Comprehensive mathematical modeling in drug addiction sciences. Drug Alcohol Depend (2007) 89(1):102–6.10.1016/j.drugalcdep.2006.12.029
49
KitanoH. Systems biology: a brief overview. Science (2002) 295(5560):1662–4.10.1126/science.1069492
50
MeyerJSQuenzerLF. Psychopharmacology: Drugs, the Brain, and Behavior. Sunderland, MA: Sinauer (2005).
51
GutkinBSDehaeneSChangeuxJ-P. A neurocomputational hypothesis for nicotine addiction. Proc Natl Acad Sci U S A (2006) 103(4):1106–11.10.1073/pnas.0510220103
52
LevyDLevyYZSiegelmannHT. Modeling addiction in terms of fluctuating neurotransmitters within the reward system. Tech Rep UMass Amherst (2011).
53
KovatchevBBretonMJohnsonB. In silico models of alcohol dependence and treatment. Front Psychiatry (2012) 3:4.10.3389/fpsyt.2012.00004
54
BradyKTSinhaR. Co-occurring mental and substance use disorders: the neurobiological effects of chronic stress. Am J Psychiatry (2005) 162(8):1483–93.10.1176/appi.ajp.162.8.1483
55
AgrawalAVerweijKJGillespieNAHeathACLessov-SchlaggarCNMartinNGet alThe genetics of addiction-a translational perspective. Transl Psychiatry (2012) 2:e140.10.1038/tp.2012.54
56
LevyYZLevyDMeyerJSSiegelmannHT. Identification and control of intrinsic bias in a multiscale computational model of drug addiction. In: ShinSYOssowskiSSchumacherMPalakalMJHungC-C editors. Proceedings of the 2010 ACM Symposium on Applied Computing. New York: ACM (2010). p. 2389–93.
57
TangY-YLuQFanMYangYPosnerMI. Mechanisms of white matter changes induced by meditation. Proc Natl Acad Sci U S A (2012) 109(26):10570–4.10.1073/pnas.1207817109
58
XueYXLuoYXWuPShiHSXueLFChenCet alA memory retrieval-extinction procedure to prevent drug craving and relapse. Science (2012) 336(6078):241–5.10.1126/science.1215070
59
MiltonALEverittBJ. Wiping drug memories. Science (2012) 336(6078):167–8.10.1126/science.1221691
60
AlpertHRConnollyGNBienerL. A prospective cohort study challenging the effectiveness of population-based medical intervention for smoking cessation. Tob Control (2013) 22(1):32–7.10.1136/tobaccocontrol-2011-050129
61
BoltDMPiperMETheobaldWEBakerTB. Why two smoking cessation agents work better than one: role of craving suppression. J Consult Clin Psychol (2012) 80(1):54–65.10.1037/a0026366
62
ThompsonGHHunterDA. Nicotine replacement therapy. Ann Pharmacother (1998) 32(10):1067–75.10.1345/aph.17382
63
FioreMCSmithSSJorenbyDEBakerTB. The effectiveness of the nicotine patch for smoking cessation. A meta-analysis. JAMA (1994) 271(24):1940–7.
64
U.S. Department of Health and Human Services, Public Health Service. Treating Tobacco Use and Dependence: 2008 Update – Clinical Practice Guideline. (2008). Available from: http://www.ahrq.gov/clinic/tobacco/treating_tobacco_use08.pdf
65
LohWYPiperMESchlamTRFioreMCSmithSSJorenbyDEet alShould all smokers use combination smoking cessation pharmacotherapy? Using novel analytic methods to detect differential treatment effects over 8 weeks of pharmacotherapy. Nicotine Tob Res (2012) 14(2):131–41.10.1093/ntr/ntr147
66
McCarthyDEPiaseckiTMJorenbyDELawrenceDLShiffmanSBakerTB. A multi-level analysis of non-significant counseling effects in a randomized smoking cessation trial. Addiction (2010) 105(12):2195–208.10.1111/j.1360-0443.2010.03089.x
67
BeckAT. Cognitive therapy. A 30-year retrospective. Am Psychol (1991) 46(4):368–75.10.1037/0003-066X.46.4.368
68
ChildressABurnsD. The basics of cognitive therapy. Psychosomatics (1981) 22(12):1017–27.10.1016/S0033-3182(81)73443-1
69
KendallPCHollonSD. Cognitive-behavioral interventions: overview and current status. In: KendallPCHollonSD editors. Cognitive-Behavioral Interventions: Theory, Research and Procedures. New York: Academic Press. (1979). p. 1–13.
70
U.S. Department of Health and Human Services, National Institutes of Health, National Center for Complementary and Alternative Medicine. What is Complementary and Alternative Medicine? (2008). Available from: http://nccam.nih.gov/sites/nccam.nih.gov/files/D347_05-25-2012.pdf
71
HendricksPSDelucchiKLHallSM. Mechanisms of change in extended cognitive behavioral treatment for tobacco dependence. Drug Alcohol Depend (2010) 109(1-3):114–9.10.1016/j.drugalcdep.2009.12.021
72
HathawayKMErasmusY. Integrative therapies for tobacco cessation. In: ChlanLHertzMI editors. Integrative Therapies in Lung Health and Sleep. New York: Springer (2012). p. 223–46.
73
BowenSChawlaNCollinsSEWitkiewitzKHsuSGrowJet alMindfulness-based relapse prevention for substance use disorders: a pilot efficacy trial. Subst Abus (2009) 30(4):295–305.10.1080/08897070903250084
74
LeshnerAI. Addiction is a brain disease, and it matters. Science (1997) 278(5335):45–7.10.1126/science.278.5335.45
75
BuchmanDZIllesJReinerPB. The paradox of addiction neuroscience. Neuroethics (2011) 4(2):65–77.10.1007/s12152-010-9079-z
76
EngelGL. The need for a new medical model: a challenge for biomedicine. Science (1977) 196(4286):129–36.10.1126/science.847460
77
LeukefeldCGLeukefeldS. Primary socialization theory and a bio/psycho/social/spiritual practice model for substance use. Subst Use Misuse (1999) 34(7):983–91.10.3109/10826089909039390
78
BuchmanDZSkinnerWIllesJ. Negotiating the relationship between addiction, ethics, and brain science. AJOB Neurosci (2010) 1(1):36–45.10.1080/21507740903508609
79
HeymanGM. Addiction: a Disorder of Choice. Cambridge, MA: Harvard UP (2010).
80
TretterFGebicke-HaerterPJAlbusMan der HeidenUSchweglerH. Systems biology and addiction. Pharmacopsychiatry (2009) 42:S11–33.10.1055/s-0029-1220699
Summary
Keywords
drug addition, allostasis, reward system adaptations, antireward system adaptations, mood, drug intake prediction, multiscale computational model, knowledge repository model
Citation
Levy YZ, Levy DJ, Barto AG and Meyer JS (2013) A Computational Hypothesis for Allostasis: Delineation of Substance Dependence, Conventional Therapies, and Alternative Treatments. Front. Psychiatry 4:167. doi: 10.3389/fpsyt.2013.00167
Received
26 June 2013
Accepted
30 November 2013
Published
19 December 2013
Volume
4 - 2013
Edited by
Lorenzo Leggio, National Institutes of Health, USA
Reviewed by
Primavera Alessandra Spagnolo, National Institutes of Health, USA; Georgiy Bobashev, RTI International, USA
Copyright
© 2013 Levy, Levy, Barto and Meyer.
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) or licensor 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: Yariv Z. Levy, School of Computer Science, University of Massachusetts Amherst, 140 Governors Drive, Amherst, MA 01003-9264, USA e-mail: ylevy@cs.umass.edu
This article was submitted to Addictive Disorders and Behavioral Dyscontrol, a section of the journal Frontiers in Psychiatry.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.