I am interested in using Neuron software to model miniature excitatory post synaptic currents (mEPSCs) built from single channel kinetic schemes. These mEPSCs would contain noise fluctuations, like those used in peak scaled non-stationary fluctuation analysis. I then want to be able to place these mEPSCs at different locations along the dendritic tree of a multi-compartmental neuron model to look at the effects of dentritic filtering on estimates of single channel conductance. Has this been done? Advice on related downloadable Neuron demos or modelling papers would be helpful.
Hello,
I think that this point has already been partly addressed .
Your NMDA model should be a point process in order to be able to locate it anywhere. A lot of NMDA mod files are available on the Neuron site (model DB). (https://senselab.med.yale.edu/modeldb/)
In addition, you should register to the Neuron forum and post any question. It's quite reactive and very helpful.
http://www.neuron.yale.edu/phpBB/
Regards
Obviously this is possible. Like everything in NEURON, you can make it as simple, or as complex as you would like.
For a vague idea, neuron comes with demo.hoc, where you can insert a stochastic single channel HH process.
But if it were me, I'd be tempted to start with a very very simple model. I'd start by coding a point process that was just one single channel. I would open syn.mod and use that as a template. I would have a state parameter of P, that would be the open probability. When t > onset, set P to some value and have it decay via P = pmax * exp(-(t-onset)/tau). It should also have a state parameter O (I think it should be state), which is whether the channel is open or not.
You can generate a random number with scop_random() (uniform random between 0 and 1). If then random number is less than P, then the channel should open (O = 1). I would then say that the channel should be open for some period of time t (presumably selected again from a random process, but could also be fixed for your purposes I believe). Once a period of time greater than that had passed, the channel should close (O = 0) and you're back to square one.
Then your BREAKPOINT block simply is i = O*g*(v-e)
Obviously, this is incredibly simple, and has lots of places where it would fall down (like non fixed time step, and it is very far from efficient) but it will give you some idea, and should in principle work for your purposes. Then you can just insert as many single channels at any given place you want.
As I've stressed, this is very very simple. Have a look in model DB at Marteen Kole's Ih model ( 64195 ) for a more complex approach to single channel dynamics.
But the approach I mention should be able to produce something like the attached, which is the sum of 8 single channels, coded with the above approach.
- Badges
- Science topic
- Similar topics
- Biological Science
- Neuroscience