L5 PFC microcircuit used to study persistent activity (Papoutsi et al. 2014, 2013)

 Download zip file   Auto-launch 
Help downloading and running models
Using a heavily constrained biophysical model of a L5 PFC microcircuit we investigate the mechanisms that underlie persistent activity emergence (ON) and termination (OFF) and search for the minimum network size required for expressing these states within physiological regimes.
1 . Papoutsi A, Sidiropoulou K, Cutsuridis V, Poirazi P (2013) Induction and modulation of persistent activity in a layer V PFC microcircuit model. Front Neural Circuits 7:161 [PubMed]
2 . Papoutsi A, Sidiropoulou K, Poirazi P (2014) Dendritic nonlinearities reduce network size requirements and mediate ON and OFF states of persistent activity in a PFC microcircuit model. PLoS Comput Biol 10:e1003764 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Dendrite; Connectionist Network;
Brain Region(s)/Organism: Neocortex;
Cell Type(s): Neocortex L5/6 pyramidal GLU cell;
Channel(s): I Na,p; I Na,t; I L high threshold; I A; I CAN; I Potassium; I R; I_AHP;
Gap Junctions:
Receptor(s): GabaA; GabaB; AMPA; NMDA;
Simulation Environment: NEURON;
Model Concept(s): Active Dendrites; Working memory;
Implementer(s): Papoutsi, Athanasia [athpapoutsi at gmail.com];
Search NeuronDB for information about:  Neocortex L5/6 pyramidal GLU cell; GabaA; GabaB; AMPA; NMDA; I Na,p; I Na,t; I L high threshold; I A; I CAN; I Potassium; I R; I_AHP;
// This function is used to balance the membrane potential to v_init 
// written by Yiota Poirazi, July 2001, poirazi@LNC.usc.edu

proc current_balance() {

	printf("Balancing each compartment to %d mV\n", $1)

	forsec pcell_list {
		for (x) {
			if (ismembrane("na_ion")) {e_pas(x)=v(x)+ina(x)/g_pas(x)} 
			if (ismembrane("k_ion")) {e_pas(x)=e_pas(x)+ik(x)/g_pas(x)} 
			if (ismembrane("ca_ion")) {e_pas(x)=e_pas(x)+ica(x)/g_pas(x)}
			if (ismembrane("in_ion")) {e_pas(x)=e_pas(x)+in(x)/g_pas(x)} 
			if (ismembrane("h")) {e_pas(x)=e_pas(x)+ihi(x)/g_pas(x)}