I'm curious if it is possible to use the rxd module to model active transport of intracellular solutes. For a simple example, suppose I have two sections:
Code: Select all
from neuron import h, rxd
soma = h.Section()
dend = h.Section()
dend.connect(soma , 1, 0)
Lets say that some mRNA transcript of ('m') interest is produced in the soma at some constant rate and degraded at a rate proportional to m. Below is my best guess as to how to implement this. (Disclaimer: I've only just started reading into the rxd module, my apologies for any errors/misunderstanding):
Code: Select all
m = rxd.Species([soma, dend], d=1, name='m')
k1 = 1.0 # rate at which m is produced in the soma
k2 = 1.0 # degradation rate of m
soma_reaction = rxd.Rate(m,k1-(k2*m),regions=soma)
Within the dendritic compartment, m is only degraded (not produced). For simplicity, lets assume the rate at which m is degraded in the dendrites is the same as in the soma.
Code: Select all
dend_reaction = rxd.Rate(m,-(k2*m),regions=dend)
The final component I'm looking for (assuming the above code is basically correct) is to have the mRNA actively transported into dend from the soma. I presume that this could be achieved by adding a drift term to the diffusion equation... Any ideas on how this could be implemented? Thanks in advance for any help.