Dear Forum,
I'd be grateful to hear people's general thoughts on the following.
I would like to use a cable model to simulate action potentials along human CNS axons for specific longrange white matter fibre bundles. The main objective is to obtain an estimate of the expected conduction delay for a given fibre bundle given three pieces of information that can now be obtained from noninvasive imaging:
1. Axon length (from diffusionweighted MRI tractography, not particularly new)
2. Axon diameter distribution (from a slightly newer diffusionweighted MRI technique based on MR signature of restricted diffusion in cylinders)
3. Gratio (from another fairly new technique  magnetization transfer imaging with bound pool fractions).
Ultimately I would like to plug these three quantities into a series of cable models (one for each bin in the diameter distribution) to get a predicted conduction delay distribution, which will then be used as an anatomicallyinformed prior on the conduction delay between the pair regions connected by the fibre bundle in question in an EEG/MEG connectivity analysis.
I have found the discussions on the details of cable models in this forum extremely helpfuly, but I am still left with the general question: is this approach (simulating human CNS conduction delays from length, diameter, and gratio) a feasible thing to try and do with a cable model? Are there any specific and/or principled reasons why this cannot / should not be done? If not, could anyone point me to the most appropriate published model (ideally on the modelDB, but also in the general literature) for this kind of analysis?
The main motivation for pursuing this is because it potentially provides a physiologically wellmotivated link between structural connectivity information and dynamic (functional/effective) connectivity models that are used to analyze EEG/MEG data. The other motivation is that I want to simulate the effects of agerelated demyelination on action potential propagation, particularly with respect to conduction delays, conduction block, and spike frequency limits. I'm fully aware of the extensive literature looking at demyelination in MS, which focuses mostly on peripheral motor pathways. One of the ways that my application differs from this I want to simulate demyelination on specific fibre bundles in individual subjects, where I know the 3D shape, trajectory, length, and (in theory) spatial distribution of demyelination. So my second question is a general one: would it be useful to construct realistic simulations with this rich geometry and microstructure information? Or would it be just as informative to look at the effects of different spatial distributions of demyelination on a set of generic 30 node axons, and try to generalize the results to specific fibre bundles with appropriate scaling etc. as a second step? Relatedly, aside from whether or not it would be useful to go for the more complex simulation, and assuming that computation time isn't an issue, would there be any actual numerical problems with using a largescale model like this, and with trying to simulate delays directly rather than calculating them from length/(simulated conduction velocity) ?
(Just FYI: So far I am using the Hursh/Rushton etc. approximation of conduction velocity = (k/g)*d, where g is the gratio, k is a scaling constant (~5.6; c.f. Caminiti et al. 2009), and d is the axon diameter. Conduction delay is then axon length/CV. Does anyone have any passionate objections to using this expression to estimate delays for fibre bundles in the human CNS when length, g, and d are known? )
Answers to my (three, I think there were) questions and general comments are greatly appreciated.
Many thanks,
JG
simulating conduction delays in human CNS

 Posts: 3
 Joined: Thu Jul 07, 2011 12:58 pm

 Site Admin
 Posts: 5471
 Joined: Wed May 18, 2005 4:50 pm
 Location: Yale University School of Medicine
 Contact:
Re: simulating conduction delays in human CNS
It is certainly "feasible" to use computational modeling to study the relationship between demyelination and conduction delay or velocity. There are plenty of models of spike initiation and propagation in myelinated axons, and at least a few of them might be a suitable starting point for such a study. You'll need to know something about how the anatomical and biophysical properties of axons are affected by demyelination, and how they change as lesions evolve; presumably there are clinical and experimental observations that can help with this aspect, and you are almost certainly far more familiar with the relevant literature than I am.
Length and diameter, yes. gratio? What's that? You're not going to find anything called gratio in any of the published axon models that I have seen.simulating human CNS conduction delays from length, diameter, and gratio