Electrotonic transform and EPSCs for WT and Q175+/- spiny projection neurons (Goodliffe et al 2018)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:236310
This model achieves electrotonic transform and computes mean inward and outward attenuation from 0 to 500 Hz input; and randomly activates synapses along dendrites to simulate AMPAR mediated EPSCs. For electrotonic analysis, in Elec folder, the entry file is MSNelec_transform.hoc. For EPSC simulation, in Syn folder, the entry file is randomepsc.hoc. Run read_EPSCsims_mdb_alone.m next with the simulated parameter values specified to compute the mean EPSC.
Reference:
1 . Goodliffe JW, Song H, Rubakovic A, Chang W, Medalla M, Weaver CM, Luebke JI (2018) Differential changes to D1 and D2 medium spiny neurons in the 12-month-old Q175+/- mouse model of Huntington's Disease. PLoS One 13:e0200626 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Synapse;
Brain Region(s)/Organism: Striatum;
Cell Type(s): Neostriatum spiny neuron;
Channel(s):
Gap Junctions:
Receptor(s): AMPA;
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Detailed Neuronal Models; Membrane Properties; Electrotonus; Synaptic-input statistic;
Implementer(s):
Search NeuronDB for information about:  AMPA;
/
GoodliffeEtAl2018
Elec
tau_tables
kir.mod *
actionPotentialPlayer.hoc *
all_tau_vecs.hoc
analyticFunctions.hoc *
aux_procs.hoc
baseline_values.txt
basic_procs.hoc
colorDendrites.hoc
electro_procs.hoc *
fixnseg.hoc *
load_scripts.hoc *
measureMeanAtten.hoc
MSN_fixDiams.hoc
MSNelect.hoc
MSNelect_transform.hoc
Nov3IR3a.hoc
Nov9IR2a_spine.hoc
readcell.hoc
                            
//First, run nrnivmodl *.mod to compile the KIR channel mod file
//Then, specify the morphology model in MSNname and mBase as the same format as the Nov3IR3a.hoc
//Run nrngui MSNelect_transform.hoc to compute inward/outward attenuation.

//load_file("nrngui.hoc")

// now load morph;  recall that 'basic_shape()' resets the morphology
strdef MSNname, mBase


MSNname = "Nov3IR3a.hoc"
mBase="Nov3IR3a"






load_file(MSNname)


objref dendritic, apical
    dendritic = new SectionList()    
//    forsec "dend"  {
	forsec proximal {
        dendritic.append()
    }





xopen("baseline_values.txt")
xopen("all_tau_vecs.hoc")
xopen("basic_procs.hoc")
xopen("load_scripts.hoc")

load_file("aux_procs.hoc")

HALF_diam = soma.diam/2
HALF_L = soma.L/2
//xopen("electro_procs.hoc")
set_pas(G_PAS)
set_kir(G_KIR)



soma printf("%s {L = %g  diam = %g}\n",secname(), L, diam)
//forall print secname()," ",L,diam


geom_nseg(500,0.1)



xopen("MSN_fixDiams.hoc")
//fix the morphology, let each dendrite with diameter<0.5 to 0.5µm

  xopen("measureMeanAtten.hoc")
  soma {
    nseg = 1
    soma_ref=new SectionRef()
  }


// {printf("File:  %s\n",MSNname)}
//{printf("File:  %s\n",MSNname)}
//{printf("%s_out = [",MSNname)}
//{printf("\nOutward Atten, dendritic\n")}
  meanOutwardAttenuationAllFrequencies(soma_ref,dendritic)
// compute mean outward attenuation

//{printf("];\n\n%s_in = [",MSNname)}
//{printf("\n\nInward Atten, dendritic\n")}
  meanInwardAttenuationAllFrequencies(soma_ref,dendritic)
// compute mean inward attenuation

{printf("];\n")}

/***
objref Mylist
forsec proximal { 
	x { 
	print secname()
	Mylist = new SectionList()
	Mylist.subtree()
	meanOutAttenAllFreqs_SecList(soma_ref,Mylist)
	dendriticLen = 0
	
	forsec Mylist {
		dendriticLen=dendriticLen + L
	}
	printf(" %g\n", dendriticLen)
}
}
***/
// This is a subroutine to compute the outward attenuation along each primary dendritic tree. 

Loading data, please wait...