Is there a model that handles Signaling Pathways in rxd?

Extending NEURON to handle reaction-diffusion problems.

Moderators: hines, wwlytton, ramcdougal

Post Reply
Posts: 37
Joined: Sun Aug 13, 2017 9:51 am

Is there a model that handles Signaling Pathways in rxd?

Post by vogdb »

I'm trying to find an example that shows usage of RxD for signaling pathways like

Code: Select all

Glu + mGluR <-> Glu_mGluR
Glu_mGlur + G <-> Glu_G_mGluR
I found some good examples in NModl like here ... mod#tabs-2 but I'm curious is there an RxD implementation.

It would be great if the example triggers Ca wave but I am happy to have any cause I couldn't find anything.
Posts: 43
Joined: Tue Apr 18, 2017 10:05 am

Re: Is there a model that handles Signaling Pathways in rxd?

Post by adamjhn »

Sorry for the delay in responding. You can use NMODL for the reactions pathway and link the output to rxd using node.include_flux, as in the example you linked to; ... .py#tabs-2 (see the setmGLURflux function).

Alternatively you can model the pathway in rxd, but you still require NMODL to implement NET_RECEIVE. Here is the an example of the same part of metabotropic glutamate receptor pathway. The segment’s glutamate concentration is linked to the NMODL file as a pointer (G) by calling h.setpointer.

With rxd species are specified in terms of concentrations, so the volume where the reactions happens is relevant. In this case it is not unreasonable to assume the reactions occur in a spine head. The spine head must have nseg=1, so that there is only one concentration for each species. The connection weight is the concentration increase in glutamate following an action potential, it is determined by number of molecules of glutamate (I arbitrarily chose 1000 for this example) and the volume of the spine head ( ~0.4μm^3)

Here is the simple mod file;

Code: Select all


    CON = (1/micrometer3) (millimole/m3)  : count/volume to uM


NET_RECEIVE (weight)  {
  G = G + weight * CON
Here is Python code that uses rxd to implement the part of the metabotropic glutamate pathway from: ... mod#tabs-2

Code: Select all

from neuron import h, rxd
from matplotlib import pyplot

#Paramters from
# species concentrations are in uM like the mod file rather than the standard mM  

initmGluR =0.3e-3   #Bhalla & Iyenger Science  1999
K1 = 0.28	        # forward binding rate to receptor from Bhalla et al
K2 = 0.016          # backward (unbinding) rate of receptor from Bhalla et al
K_PLC = 5           # total concentration of PLC 
K_PIP2 = 160        # total concentration of PIP2 
K_G=25              #
     #kplc and Vmax describe aPLC catalyzing IP3 production from PIP2
kfplc = 0.83
kbplc = 0.68        #0.1/ms in the paper; added to Vmax1=0.58/ms in the paper
Vmax1 = 0.58      
     #D5 and D6 describe Glu_mGluR catalyzing G_alpha production, Km2=(D6f+D5B)/D5f
D5f = 15 
D5b = 7.2 
D6f = 1.8
    #G2 describe aG binding to PLC
G2f = 100 
G2b = 100 
    #degradation of aG (D7f) and IP3 (G9f)
D7f = 9
degGluRate = 1.0 # should lookup value

# Spine Head
head = h.Section()
head.nseg = 1
head_L = 0.8 #um
head_d = 0.8 #um
h.pt3dadd(0, 0, 0, head_d, sec=head)
h.pt3dadd(0, 0, head_L, head_d, sec=head) 
glu_per_spike = 1000   # number of molecules (arbitrarily choice for this example) 

# RxD
spine_head = rxd.Region([head],nrn_region='i')

Glu = rxd.Species(spine_head,name="Glu",initial=0)
mGluR = rxd.Species(spine_head,name="mGluR",initial=initmGluR)
Glu_mGluR = rxd.Species(spine_head,name="Glu_mGluR",initial=0)
react1 = rxd.Reaction(Glu + mGluR <> Glu_mGluR,K1,K2)

degGlu = rxd.Species(spine_head,name="degGlu",initial=0)
react2 = rxd.Reaction(Glu > degGlu,degGluRate)

G = rxd.Species(spine_head,name="G",initial=K_G)
GG_mGluR = rxd.Species(spine_head,name="GG_mGluR",initial=0)
react3 = rxd.Reaction(Glu_mGluR + G <> GG_mGluR,D5f,D5b)

aG = rxd.Species(spine_head,name="aG",initial=0)
react4 = rxd.Reaction(GG_mGluR > aG + mGluR,D6f)
react5 = rxd.Reaction(aG > G,D7f)

PLC = rxd.Species(spine_head,name="PLC",initial=K_PLC)
aPLC_aG = rxd.Species(spine_head,name="aPLC_aG",initial=0)
react6 = rxd.Reaction(aG + PLC <> aPLC_aG, G2f, G2b)

PIP2 = rxd.Species(spine_head,name="PIP2",initial=K_PIP2)
aPLC_PIP2 = rxd.Species(spine_head,name="aPLC_PIP2",initial=0)

react7 = rxd.Reaction(aPLC_aG+PIP2 <> aPLC_PIP2,kfplc,kbplc)

ip3 = rxd.Species(spine_head,name="ip3",d=1.415,initial=0)
react8 = rxd.Reaction(aPLC_PIP2 > ip3,Vmax1)

#ip3 degradation has been omitted in this example

# Add the synapse
syn = h.mGLURRxD(head(0.5))

# Connect the RxD concentration with the point process

# Noisy stimulation
stim = h.NetStim() 
stim.number = 100
stim.start = 0
stim.noise = 1
ncstim = h.NetCon(stim, syn)
ncstim.delay = 1
ncstim.weight[0] = glu_per_spike/(Glu.nodes[0].volume) # number per um^3

# Record
tvec = h.Vector()
Glu_vec = h.Vector()
ip3_vec = h.Vector()

# Run simulation

# Plot the result
pyplot.xlabel("t (ms)")
pyplot.ylabel("Glu (uM)")

pyplot.xlabel("t (ms)")
pyplot.ylabel("IP3 (nM)")
Here is the output;
Post Reply