A computational model of a small DRG neuron to explore pain (Verma et al. 2019, 2020)

 Download zip file 
Help downloading and running models
This is a Hodgkin-Huxley type model for a small DRG neuron consisting of four voltage-gated ion channels: sodium channels 1.7 and 1.8, delayed rectifier potassium, and A-type transient potassium channels. This model was used to explore the dynamics of this neuron using bifurcation theory, with the motive to investigate pain since small DRG neuron is a pain-sensing neuron.
1 . Verma P, Kienle A, Flockerzi D, Ramkrishna D (2020) Using Bifurcation Theory for Exploring Pain Industrial & Engineering Chemistry Research 59:2524-2535
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism:
Cell Type(s): Dorsal Root Ganglion (DRG) cell;
Channel(s): I Sodium; I Potassium;
Gap Junctions:
Gene(s): Nav1.7 SCN9A; Nav1.8 SCN10A;
Simulation Environment: MATLAB; XPP;
Model Concept(s): Action Potentials; Bifurcation;
Implementer(s): Verma, Parul [parulv.iitb at gmail.com];
Search NeuronDB for information about:  I Sodium; I Potassium;
Verma P, Kienle A, Flockerzi D, Ramkrishna D. Using Bifurcation Theory for Exploring Pain. Industrial & Engineering Chemistry Research. 2019.

Verma P, Kienle A, Flockerzi D, Ramkrishna D. Computational analysis of a 9D model for a small DRG neuron. arXiv preprint arXiv:2001.04915. 2020 Jan 14.

Verma P, Eaton M, Kienle A, Flockerzi D, Yang Y, Ramkrishna D. A mathematical investigation of chemotherapy-induced peripheral neuropathy. bioRxiv. 2020 Jan 1.

Both XPP (DRG.ode) and MATCONT (*.m) files are provided in the respective folders.

In MATCONT, analytical Jacobian was used. You can run DRG_Jacobian.m to generate the Jacobian matrix functions (J_D.m, JhKA_D.m, JN_D.m, JP_D.m). DRG.m contains the equations. DRG_Steady_Periodic_I.m can be used to perform continuation of equilibrium and periodic solutions for the set of equations. Note that an extra parameter V12 is introduced in the MATCONT equations. This is kept as zero, but can be varied as was done for the first aforementioned paper.

The MATCONT equations are scaled, while XPP equations are not.

20200521 An update from Parul Verma supplied another XPP file and the
bioRxiv citation was added to this Readme

Loading data, please wait...