-
Notifications
You must be signed in to change notification settings - Fork 1
HH-type model of fast-spiking parvalbumin interneurons in spinal dorsal horn (Ma et al, 2023)
ModelDBRepository/267614
Folders and files
| Name | Name | Last commit message | Last commit date | |
|---|---|---|---|---|
Repository files navigation
<html>
<head>
<style>
th, td {
border-color: #FFFFFF;
}
th, tfoot td {
font-weight: normal;
text-align: left
}
</style>
</head>
<body>
<p>This is the readme for the models associated with the manuscript in preparation: </p>
<p> "Ma, X., Miraucourt, L., Qiu, H., Sharif-Naeini, R., Khadra, A. (2023). Calcium buffering tunes intrinsic excitability of spinal dorsal horn parvalbumin-expressing interneurons: A computational model." </p>
<p>The code can be used to produce the time series simulations and bifurcation analysis of PVIN model.</p>
<p>
---------------------------------------------<br>
Author: Xinyue Ma <br>
Email: xinyue.ma@mail.mcgill.ca <br>
Integrated Program in Neuroscience <br>
McGill University <br>
Montreal, QC, H3A 1A1 <br>
Canada <br>
---------------------------------------------
</p>
<h2>I. Description</h2>
<p>
The PVIN model is a Hodgkin-Huxley type model of parvalbumin-expressing interneurons (PVINs) in the dorsal horn of the spinal cord. It consists ion channels including fast-transient sodium channel (Nav), slow delayed-rectifier potassium channel (Kv1), fast delayed-rectifier potassium channel (Kv3), high voltage-gated calcium channel (Cav), small conductance calcium-activated potassium channel (SK), and a leak current. It also uses the flux-balanced equation describing the intracellular calcium dynamics. The PVIN model is simulated using a stable step current protocol and a randomized synaptic current protocol. The code of PVIN model written in Matlab and XPP are both provided.
</p>
<p>
The model is adopted from Bischop et al. 2012 and reparametrized to fit the electrical activity of spinal dorsal horn PVINs from naive mice. It reproduces the firing pattern change of PVINs from tonic to transient following nerve injury, which is realized by reducing cytosolic calcium buffer concentration. The bifurcation analysis of PVIN model further explains how the firing pattern transits as the injection current increases, in a manner similar to that is seen in our transient firing PVIN recordings.
</p>
<p>
The code also includes a “in vivo-like” neural circuit model of Aβ fiber-mediate nociceptive neural circuit. The circuit model is stimulated by Poisson-distributed excitatory synaptic currents representing the presynaptic inputs from the Aβ fibers. It includes the PVIN model above, as well as another HH type model describing the excitability of a PVIN post-synaptic target: the excitatory interneuron expressing protein kinase C gamma (PKCγIN). The PKCγIN model is adapted from Medlock et al. 2022. The Aβ fiber-like presynaptic current was applied on both the inhibitory PVIN model and the excitatory PKCγIN model, the latter of which also received inhibitory synaptic input from the PVIN model.
</p>
<p>
<b>Reference:</b> <br>
<p> Bischop, D. P., Orduz, D., Lambot, L., Schiffmann, S. N., & Gall, D. (2012). Control of neuronal excitability by calcium binding proteins: a new mathematical model for striatal fast-spiking interneurons. Frontiers in molecular neuroscience, 5, 78.</p>
<p> Medlock L, Sekiguchi K, Hong S, Dura-Bernal S, Lytton WW, Prescott SA (2022) Multiscale Computer Model of the Spinal Dorsal Horn Reveals Changes in Network Processing Associated with Chronic Pain. J Neurosci 42:3133–3149.
</p>
</p>
<h2>II. Simulators</h2>
MATLAB Version: 9.12.0 (R2022a) <br>
XPPAUT Version: 8.0 <br>
Time-stamp: 2023-March-05 <br>
<h2>III. File list</h2>
<h3>XPP</h3>
<p>
<table style="width:100%">
<tr>
<th>PVIN_2023.ode</th>
<th>PVIN model for time series simulations</th>
</tr>
</table>
<!--
 PVIN_2023.ode   PVIN model for time series simulations and one-parameter bifurcation analysis
//-->
</p>
<h3>MATLAB</h3>
<p>
<table style="width:100%">
<tr>
<th>PVIN_STEP_Simu.m</th>
<th>Example code to generate step current-stimulated voltage response of PVINs, as in Figs. 2,3 & 5 in Ma et al. 2023</th>
</tr>
<tr>
<th>PVIN_Syn_Simu.m</th>
<th>Example code to generate Poisson-distributed synaptic current-stimulated voltage response of PVINs, as in Fig. 9 in Ma et al. 2023</th>
</tr>
<tr>
<th>PVIN_TwoParBifur.m</th>
<th>Example code to reproduce the two-parameter bifurcation diagram shown in Fig. 8 in Ma et al. 2023</th>
</tr>
<tr>
<th>NeuralCircuit.mlx</th>
<th>Example code to simulate the spinal dorsal neural circuit model, as in Fig. 9 in Ma et al. 2023</th>
</tr>
<tr>
<th>PVIN_HH.m</th>
<th>ODEs of PVIN model</th>
</tr>
<tr>
<th>ePKCmodel.m</th>
<th>ODEs of ePKCγ interneuron model</th>
</tr>
<tr>
<th>PVIN_Cai.m</th>
<th>Function: MATCONT equations of PVIN model</th>
</tr>
<tr>
<th>PVIN_Cai.mat</th>
<th>Function: MATCONT variables of PVIN model</th>
</tr>
<tr>
<th>runHHmodel_STEP.m</th>
<th>Function: Step current stimulation protocol of PVIN model</th>
</tr>
<tr>
<th>runHHmodel_AbetaPoisson.m</th>
<th>Function: Synaptic current stimulation protocol of PVIN model</th>
</tr>
<tr>
<th>mySynInput.m</th>
<th>Function: Generate the time series of AMPA and NMDA conductances at customized firing rate</th>
</tr>
<tr>
<th>poissonSpikeGen.m</th>
<th>Function: Generate Poisson-distributed timing of spike trains</th>
</tr>
<tr>
<th>factor_syn.m</th>
<th>Function: Generate a factor to correct the exponential equation peaks at 1</th>
</tr>
<tr>
<th>genSyn.m</th>
<th>Function: Generate the time series of a customized synaptic conductance</th>
</tr>
</table>
<!--
 PVIN_STEP_Simu.m   Example code to generate step current-stimulated voltage response of PVINs, as shown in Fig. 2-4 <br>
 PVIN_Syn_Simu.m   Example code to generate Poisson-distributed synaptic current-stimulated voltage response of PVINs, as shown in Fig. 7 <br>
 PVIN_TwoParBifur.m   Example code to reproduce the two-parameter bifurcation diagram in Fig. 6 <br>
 NeuralCircuit.mlx   Example code to simulate the spinal dorsal neural circuit model, as shown in Fig. 7 <br>
 PVIN_HH.m   ODEs of PVIN model <br>
 ePKCmodel.m   ODEs of ePKCγ interneuron model <br>
 PVIN_Cai.m   Function: MATCONT equations of PVIN model <br>
 PVIN_Cai.mat   Function: MATCONT variables of PVIN model <br>
 runHHmodel_STEP.m   Function: Step current stimulation protocol of PVINs <br>
 runHHmodel_AbetaPoisson.m   Function: Synaptic current stimulation protocol of PVINs <br>
 mySynInput.m   Function: Generate the time series of AMPA and NMDA conductances at customized firing rate <br>
 poissonSpikeGen.m   Function: Generate Poisson-distributed timing of spike trains <br>
 factor_syn.m   Function: Generate a factor to correct the exponential equation peaks at 1
 genSyn.m   Function: Generate the time series of a customized synaptic conductance
//-->
</p>
<h2>IV. Instruction</h2>
<h3>XPPAUT<h3>
<h4> Simulate the voltage response of PVIN model to step current stimulation</h4>
<ol>
<li>Open <code>PVIN_2023.ode</code> with XPPAUT</li>
<li>Parameter specification: injection current value <br>
click 'Parameters', input 'iapp', input the value (e.g., 200), press 'Enter'
</li>
<li>Parameter specification: calcium buffer concentration <br>
click 'Parameters', input 'bt', input the value (e.g., 90), press 'Enter'
</li>
<li>Run the simulation <br>
click 'Initialconds' + '(G)o', and the result is like below <br>
<img src="./Figures/example_xpp.png" alt="PVIN model step current (xpp)" style="width:560px;height:400px;">
</li>
</ol>
<h3>Matlab<h3>
<h4> Simulate the voltage response of PVIN model to step current stimulation</h4>
<ol>
<li>Run <code>PVIN_STEP_Simu.m</code></li>
<li>In the command window, input the condition: 1 - naive PVIN | 2 - CCI PVIN </li>
<li>In the command window, input the value of applied step current value (pA) </li>
<li>A figure like below appears. The top panel shows voltage simulation and the bottom one is the applied current.</li>
</ol>
<p><img src="./Figures/example_step.png" alt="PVIN model step current stimulation" style="width:560px;height:400px;"></p>
<h4> Simulate the voltage response of PVINs to synaptic current stimulation</h4>
<ol>
<li>Run <code>PVIN_Syn_Simu.m</code></li>
<li>In the command window, input the condition: 1 - naive PVIN | 2 - CCI PVIN </li>
<li>In the command window, input PVIN presynaptic current firing rate (Hz) </li>
<li>A figure like below appears. The top panel shows voltage simulation and the bottom one is the presyaptic conductances.</li>
</ol>
<p><img src="./Figures/example_syn.png" alt="PVIN model synaptic current stimulation" style="width:560px;height:400px;"></p>
<h4> Reproduce the two-parameter bifurcation diagram shown in Fig. 8</h4>
<ol>
<li>Run <code>PVIN_TwoParBifur.m</code></li>
<li>After a few seconds, the two-parameter bifurcation diagram of PVIN model appears</li>
</ol>
<p><img src="./Figures/example_twopar.png" alt="PVIN model two-parameter bifurcation analysis" style="width:600px;height:600px;"></p>
<h4> Simulate the PVIN embedded neural circuit model</h4>
<p>Open and run <Code>NeuralCircuit.mlx</Code>. This is a livescript presenting the simulation steps of the neural circuit model associated with Fig. 9 in Ma et al. 2023. The steps include: <br>
<ol>
<li>Create Poisson-distributed synaptic conductances representing the presynaptic input from Aβ fibers </li>
<img src="./Figures/example_nc_abeta.png" alt="neural circuit model: Abeta fibersr" style="width:600px;height:a00px;"></p>
<li>Simulate the PVIN model using synaptic currents from Aβ fibers</li>
<img src="./Figures/example_nc_pvin.png" alt="neural circuit model: PVIN" style="width:600px;height:100px;"></p>
<li>Simulate the PKCγ model using synaptic currernts from Aβ fibers and PVIN model</li>
<img src="./Figures/example_nc_pkcg.png" alt="neural circuit model: PKCgammaIN" style="width:600px;height:100px;"></p>
</ol>
<p>
</p>
**************************************************************************************************
</body></html>
About
HH-type model of fast-spiking parvalbumin interneurons in spinal dorsal horn (Ma et al, 2023)