A Multiscale, Systems-level, Neuropharmacological Model of Cortico-Basal Ganglia System for Arm Reaching under Normal, Parkinsonian and Levodopa Medication Conditions

Abstract

The root cause of Parkinson's disease (PD) is the death of dopaminergic neurons in Substantia Nigra pars compacta (SNc). The exact cause of this cell death is still not known. Loss of SNc cells manifest as the cardinal symptoms of PD, including tremor, rigidity, bradykinesia, and postural imbalance. To investigate the PD condition in detail and understand the link between loss of cells in SNc and PD symptoms, it is important to have an integrated multiscale computational model that can replicate the symptoms at the behavioural level by evoking the key cellular and molecular underlying mechanisms that contribute to the pathology. In line with this objective, we present a multiscale integrated model of cortico-basal ganglia motor circuitry for arm reaching task, incorporating a detailed biophysical model of SNc dopaminergic neuron. Earlier researchers have shown that fluctuations in dopamine (DA) signals are analogous to reward/punishment signals, thereby prompting application of concepts from reinforcement learning (RL) to modelling the basal ganglia system. In our model, we replace the abstract representations of reward with the realistic variable of extracellular DA released by a network of SNc cells and incorporate it with the RL-based behavioural model, which simulates the arm reaching task. Our results showed that as SNc cell loss increases, the percentage success rate to reach the target decreases, and average time to reach the target increases. With levodopa (L-DOPA) medication, both the success rate and the average time to reach the target improved significantly. The proposed model also exhibits how differential dopaminergic axonal degeneration in basal ganglia results in various cardinal symptoms of PD as manifest in reaching movements. From the model results, we were able to show the side effects of L-DOPA mediation, such as wearing off and peak dosage dyskinesias. Moreover, from the results, we were able to predict the optimum dosage for varying degrees of cell loss and L-DOPA medication. The proposed model has a potential clinical application where drug dosage can be optimized as per patient characteristics. We conclude that our model presents a realistic and efficient way of simulating the PD pathology conditions and the effect of levodopa medication, thereby giving a reliable indicator towards the optimization of the drug dosage.