Fibromyalgia (FM) is a central disorder characterized by chronic pain, fatigue, insomnia, depression, and other minor symptoms. Knowledge about pathogenesis is lacking, diagnosis difficult, clinical approach puzzling, and patient management disappointing. We conducted a theoretical study based on literature data and computational analysis, aimed at developing a comprehensive model of FM pathogenesis and addressing suitable therapeutic targets. We started from the evidence that FM must involve a dysregulation of central pain processing, is female prevalent, suggesting a role for the hypothalamus-pituitary-gonadal (HPG) axis, and is stress-related, suggesting a role for the HP-adrenocortical (HPA) axis. Central pathogenesis was supposed to involve a pain processing loop system including the thalamic ventroposterolateral nucleus (VPL), the primary somatosensory cortex (SSC), and the thalamic reticular nucleus (TRN). For decreasing GABAergic and/or increasing glutamatergic transmission, the loop system crosses a bifurcation point, switching from monostable to bistable, and converging on a high-firing-rate steady state supposed to be the pathogenic condition. Thereafter, we showed that GABAergic transmission is positively correlated with gonadal-hormone-derived neurosteroids, notably allopregnanolone, whereas glutamatergic transmission is positively correlated with stress-induced glucocorticoids, notably cortisol. Finally, we built a dynamic model describing a multistable, double-inhibitory loop between HPG and HPA axes. This system has a high-HPA/low-HPG steady state, allegedly reached in females under combined premenstrual/postpartum brain allopregnanolone withdrawal and stress condition, driving the thalamocortical loop to the high-firing-rate steady state, and explaining the connection between endocrine and neural mechanisms in FM pathogenesis. Our model accounts for FM female prevalence and stress correlation, suggesting the use of neurosteroid drugs as a possible solution to currently unsolved problems in the clinical treatment of the disease.

Mal de Debarquement Syndrome (MdDS) is a puzzling central vestibular disorder characterized by a long-lasting perception of oscillatory postural instability that may occur after sea travels or flights. We have postulated that MdDS originates from the post-disembarking persistence of an adaptive internal oscillator consisting of a loop system, involving the right and left vestibular nuclei, and the Purkinje cells of the right and left flocculonodular cerebellar cortex, connected by GABAergic and glutamatergic fibers. We have formulated here a mathematical model of the vestibulo-cerebellar loop system and carried out a computational analysis based on a set of differential equations describing the interactions among the loop elements and containing Hill functions that model input-output firing rates relationships among neurons. The analysis indicates that the system acquires a spontaneous and permanent oscillatory behavior for a decrease of threshold and an increase of sensitivity in neuronal input-output responses. These results suggest a role for synaptic plasticity in MdDS pathophysiology, thus reinforcing our previous hypothesis that MdDS may be the result of excessive synaptic plasticity acting on the vestibulo-cerebellar network during its entraining to an oscillatory environment. Hence, our study points to neuroendocrine pathways that lead to increased synaptic response as possible new therapeutic targets for the clinical treatment of the disorder.

Fibromyalgia (FM) is an unsolved central pain processing disturbance. We aim to provide a unifying model for FM pathogenesis based on a loop network involving thalamocortical regions, i.e., the ventroposterior lateral thalamus (VPL), the somatosensory cortex (SC), and the thalamic reticular nucleus (TRN). The dynamics of the loop have been described by three differential equations having neuron mean firing rates as variables and containing Hill functions to model mutual interactions among the loop elements. A computational analysis conducted with MATLAB has shown a transition from monostability to bistability of the loop behavior for a weakening of GABAergic transmission between TRN and VPL. This involves the appearance of a high-firing-rate steady state, which becomes dominant and is assumed to represent pathogenic pain processing giving rise to chronic pain. Our model is consistent with a bulk of literature evidence, such as neuroimaging and pharmacological data collected on FM patients, and with correlations between FM and immunoendocrine conditions, such as stress, perimenopause, chronic inflammation, obesity, and chronic dizziness. The model suggests that critical targets for FM treatment are to be found among immunoendocrine pathways leading to GABA/glutamate imbalance having an impact on the thalamocortical system.

Amyotrophic lateral sclerosis (ALS) is a poor-prognosis disease with puzzling pathogenesis and inconclusive treatments. We develop a mathematical model of ALS based on a system of interactive feedback loops, focusing on the mutant SOD1^G93A mouse. Misfolded mutant SOD1 aggregates in motor neuron (MN) mitochondria and triggers a first loop characterized by oxidative phosphorylation impairment, AMP kinase over-activation, 6-phosphofructo-2-kinase (PFK3) rise, glucose metabolism shift from pentose phosphate pathway (PPP) to glycolysis, cell redox unbalance, and further worsening of mitochondrial dysfunction. Oxidative stress then triggers a second loop, involving the excitotoxic glutamatergic cascade, with cytosolic Ca²⁺ overload, increase of PFK3 expression, and further metabolic shift from PPP to glycolysis. Finally, cytosolic Ca²⁺ rise is also detrimental to mitochondria and oxidative phosphorylation, thus closing a third loop. These three loops are overlapped and positive (including an even number of inhibitory steps), hence they form a candidate multistationary (bistable) system. To describe the system dynamics, we model the interactions among the functional agents with differential equations. The system turns out to admit two stable equilibria: the healthy state, with high oxidative phosphorylation and preferential PPP, and the pathological state, with AMP kinase activation, PFK3 over expression, oxidative stress, excitotoxicity and MN degeneration. We demonstrate that the loop system is monotone: all functional agents consistently act toward the healthy or pathological condition, depending on low or high mutant SOD1 input. We also highlight that molecular interactions involving PFK3 are crucial, as their deletion disrupts the system’s bistability leading to a single healthy equilibrium point. Hence, our mathematical model unveils that promising ALS management strategies should be targeted to mechanisms that keep low PFK3 expression and activity within MNs.

Introduction: Mal de Debarquement Syndrome (MdDS) is a poorly understood neurological disorder affecting mostly perimenopausal women. MdDS has been hypothesized to be a maladaptation of the vestibulo-ocular reflex, a neuroplasticity disorder, and a consequence of neurochemical imbalances and hormonal changes. Our hypothesis considers elements from these theories, but presents a novel approach based on the analysis of functional loops, according to Systems and Control Theory. Hypothesis: MdDS is characterized by a persistent sensation of self-motion, usually occurring after sea travels. We assume the existence of a neuronal mechanism acting as an oscillator, i.e., an adaptive internal model, that may be able to cancel a sinusoidal disturbance of posture experienced aboard, due to wave motion. Thereafter, we identify this mechanism as a multi-loop neural network that spans between vestibular nuclei and the flocculonodular lobe of the cerebellum. We demonstrate that this loop system has a tendency to oscillate, which increases with increasing strength of neuronal connections. Therefore, we hypothesize that synaptic plasticity, specifically long-term potentiation, may play a role in making these oscillations poorly damped. Finally, we assume that the neuromodulator Calcitonin Gene-Related Peptide, which is modulated in perimenopausal women, exacerbates this process thus rendering the transition irreversible and consequently leading to MdDS. Conclusion and Validation: The concept of an oscillator that becomes noxiously permanent can be used as a model for MdDS, given a high correlation between patients with MdDS and sea travels involving undulating passive motion, and an alleviation of symptoms when patients are re-exposed to similar passive motion. The mechanism could be further investigated utilizing posturography tests to evaluate if subjective perception of motion matches with objective postural instability. Neurochemical imbalances that would render individuals more susceptible to developing MdDS could be investigated through hormonal profile screening. Alterations in the connections between vestibular nuclei and cerebellum, notably GABAergic fibers, could be explored by neuroimaging techniques as well as transcranial magnetic stimulation. If our hypothesis were tested and verified, optimal targets for MdDS treatment could be found within both the neural networks and biochemical factors that are deemed to play a fundamental role in loop functioning and synaptic plasticity.

We performed a mathematical analysis of the dynamic control loops regulating the vasomotor tone of vascular smooth muscle, blood volume, and mean arterial pressure, which involve the arginine vasopressin (AVP) system, the atrial natriuretic peptide system (ANP), and the renin-angiotensin-aldosterone system (RAAS). Our loop analysis of the AVP-ANP-RAAS system revealed the concurrent presence of two different regulatory mechanisms, which perform the same qualitative function: one affects blood pressure by regulating vasoconstriction, the other by regulating blood volume. Both the systems are candidate oscillators consisting of the negative-feedback loop of a monotone system: they admit a single equilibrium that can either be stable or give rise to oscillatory instability. Also a subsystem, which includes ANP and AVP stimulation of vascular smooth muscle cells, turns out to be a candidate oscillator composed of a monotone system with multiple negative feedback loops, and we show that its oscillatory potential is higher when the delays along all feedback loops are comparable. Our results give insight into the physiological mechanisms ruling long-term homeostasis of blood hydraulic parameters, which operate based on dynamical loops of interactions.