Ventromedial Thalamocortical Neuron (Bichler et al 2021)

 Download zip file 
Help downloading and running models
"Biophysical computer modeling of a thalamic neuron demonstrated that an increase in rebound spiking can also be accounted for by a decrease in the M-type potassium current. Modeling also showed that an increase in sag with hyperpolarizing steps found after 6-OHDA treatment could in part but not fully be accounted for by the decrease in M-type current. These findings support the hypothesis that homeostatic changes in BGMT neural properties following 6-OHDA treatment likely influence the signal processing taking place in the BG thalamocortical network in Parkinson's disease."
1 . Bichler EK, Cavarretta F, Jaeger D (2021) Changes in Excitability Properties of Ventromedial Motor Thalamic Neurons in 6-OHDA Lesioned Mice. eNeuro [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism: Thalamus;
Cell Type(s):
Channel(s): I A; I h; I L high threshold; I T low threshold; I K,Ca; I Sodium; I Potassium; I Na,p; I K;
Gap Junctions:
Simulation Environment: NEURON; Python;
Model Concept(s): Pathophysiology;
Search NeuronDB for information about:  I Na,p; I L high threshold; I T low threshold; I A; I K; I h; I K,Ca; I Sodium; I Potassium;
from plots import Plot, VPlot
from pynwb import NWBHDF5IO
import fi_curve_analysis as FI
import numpy as np

class FIplot(Plot):
  def __init__(self, ax, markersize=7.5, linewidth=2, markeredgecolor='black', markeredgewidth=2, labelsize=10, fontweight='normal'):
    super(FIplot, self).__init__(ax,
      xlabel='Current (pA)',
      ylabel='Frequency (Hz)',
      legend=True)[-2.5, 30])[0,400]), 31, 10).tolist()), 401, 100).tolist())







def plot_fi_curve(ax, label=''):
  io = NWBHDF5IO('fi_curve.nwb', 'r')
  nwbfile =
  _fi_plot = FIplot(ax)
  ii, ff = FI.get_fi_curve(nwbfile, "control_%g", np.arange(0.03, 0.36, 0.01))
  _fi_plot.plot(ii*1000, ff, color='black', label='Control')
  ii, ff = FI.get_fi_curve(nwbfile, "ih2_%g", np.arange(0.03, 0.36, 0.01))
  _fi_plot.plot(ii*1000, ff, color='pink', label='+200% ih')
  ii, ff = FI.get_fi_curve(nwbfile, "ih3_%g", np.arange(0.03, 0.36, 0.01))
  _fi_plot.plot(ii*1000, ff, color='red', label='+300% ih')


def plot_vm_curve(ax, shift=0.05, stim_start=5000, stim_end=7000, xlim=[5000-50,5760], key='control', color='black', title=''):
  io = NWBHDF5IO('fi_curve.nwb', 'r')
  nwbfile =
  _fi_plot = VPlot(ax, title=title)
  _ii_rheobase_control = shift+FI.get_rheobase(
    *FI.get_fi_curve(nwbfile, "control_%g", np.arange(0.01, 0.36, 0.01))
  trace = FI.analysis.read_trace(nwbfile, "%s_%g" % (key, _ii_rheobase_control), stim_start=stim_start, stim_end=stim_end)
  _fi_plot.plot(trace['T'], trace['V'], xlim=xlim, color=color)



Loading data, please wait...