Cognitive flexibility is essential to modify our behavior in a non-stationary environment and is often explored by reversal learning tasks. The basal ganglia (BG) dopaminergic system, under a top-down control of the pre-frontal cortex, is known to be involved in flexible action selection through reinforcement learning. However, how adaptive dopamine changes regulate this process and learning mechanisms for training the striatal synapses remain open questions. The current study uses a neurocomputational model of the BG, based on dopamine-dependent direct (Go) and indirect (NoGo) pathways, to investigate reinforcement learning in a probabilistic environment through a task that associates different stimuli to different actions. Here, we investigated: the efficacy of several versions of the Hebb rule, based on covariance between pre- and post-synaptic neurons, as well as the required control in phasic dopamine changes crucial to achieving a proper reversal learning. Furthermore, an original mechanism for modulating the phasic dopamine changes is proposed, assuming that the expected reward probability is coded by the activity of the winner Go neuron before a reward/punishment takes place. Simulations show that this original formulation for an automatic phasic dopamine control allows the achievement of a good flexible reversal even in difficult conditions. The current outcomes may contribute to understanding the mechanisms for active control of dopamine changes during flexible behavior. In perspective, it may be applied in neuropsychiatric or neurological disorders, such as Parkinson’s or schizophrenia, in which reinforcement learning is impaired.

Reinforcement learning, crucial for behavior in dynamic environments, is driven by rewards and punishments, modulated by dopamine (DA) changes. This study explores the dopaminergic system's influence on learning, particularly in Parkinson's disease (PD), where medication leads to impaired adaptability. Highlighting the role of tonic DA in signaling the valence of actions, this research investigates how DA affects response vigor and decision-making in PD. DA not only influences reward and punishment learning but also indicates the cognitive effort level and risk propensity in actions, which are essential for understanding and managing PD symptoms. In this work, we adapt our existing neurocomputational model of basal ganglia (BG) to simulate two reversal learning tasks proposed by Cools et al. We first optimized a Hebb rule for both probabilistic and deterministic reversal learning, conducted a sensitivity analysis (SA) on parameters related to DA effect, and compared performances between three groups: PD-ON, PD-OFF, and control subjects. In our deterministic task simulation, we explored switch error rates after unexpected task switches and found a U-shaped relationship between tonic DA levels and switch error frequency. Through SA, we classify these three groups. Then, assuming that the valence of the stimulus affects the tonic levels of DA, we were able to reproduce the results by Cools et al. As for the probabilistic task simulation, our results are in line with clinical data, showing similar trends with PD-ON, characterized by higher tonic DA levels that are correlated with increased difficulty in both acquisition and reversal tasks. Our study proposes a new hypothesis: valence, signaled by tonic DA levels, influences learning in PD, confirming the uncorrelation between phasic and tonic DA changes. This hypothesis challenges existing paradigms and opens new avenues for understanding cognitive processes in PD, particularly in reversal learning tasks.

Virtual population (VPop) generation is a central component of quantitative systems pharmacology (QSP), involving the sampling of parameter sets that represent physiologically plausible patients and capture observed inter-individual variability in clinical outcomes. This approach poses challenges due to the high dimensionality and often non-identifiability nature of many QSP models. In this study, we evaluate the performance of the DREAM(ZS) algorithm, a multi-chain adaptive Markov chain Monte Carlo (MCMC) method for generating VPop. Using the Van De Pas model of cholesterol metabolism as a case study, we compare DREAM(ZS) to the single-chain Metropolis-Hastings (MH) algorithm introduced by Rieger et al. Our comparison focuses on convergence behavior, parametric diversity, and posterior coverage, in relation to the ability of each method to explore complex parameter distributions and maintain correlations. DREAM(ZS) demonstrates superior exploration of the parameter space, reducing boundary accumulation effects common in traditional MH sampling, and restoring parameter correlation structures. These advantages are attributed in part to its adaptive proposal mechanism and the use of a bias-corrected likelihood formulation, which together contribute to a better parameters space sampling without compromising model fit. Our findings contribute to the ongoing development of efficient sampling methodologies for high-dimensional biological models, introducing a promising and easy to use alternative for VPop generation in QSP, expanding the methodological approaches for in silico trial simulation.