Code to Replicate the Computational Results from Romano et al. Nature Communications 2025

The computational model and corresponding code have been designed and implemented by Jorge Mejias, from the University of Amsterdam.
Contact: j.f.mejias@uva.nl

---

Instructions

Follow these steps to run the code and obtain the simulation results:

1. Run main1.m:
   Modify any parameters within this file to reproduce the different results shown in the manuscript.

2. Replicating Figures:

   • Figure 1:
     Run main.m for the following conditions:

     • Premovement: PFCsig = 100
     • Movement: PFCsig = 0
       Compare the results by running the script figure1.m.

   • Figures 2 and 3:
     Inactivate different populations (e.g., PN or VL) using the variable optopulse.

     • For inactivation, use negative input values.
       Analyze the results by building and using similar scripts.

   • Figure 4:
     To simulate the stimulation of PC, simply run figure4.m.