Enable selecting to show driv force, m3h plot in notebook.

Author: pgleeson

Tutorial/Source/HodgkinHuxley.py

@@ -15,8 +15,8 @@ def __init__(self, C_m=1, gmax_Na=120, gmax_K=36, gmax_L=0.3, E_Na=50,
                  I_inj_amplitude=0, I_inj_duration=0, I_inj_delay=0,
                  vc_delay=10, vc_duration=30, vc_condVoltage=-65,
                  vc_testVoltage=10, vc_returnVoltage=-65, runMode='iclamp',
-                 injected_current_plot=True, gating_plot=True, scaling_plot=True,
-                 cond_dens_plot=True, driving_force_plot=True,
+                 injected_current_plot=True, gating_plot=True, cond_scaling_plot=False,
+                 cond_dens_plot=True, driving_force_plot=False,
                  current_plot=True, memb_pot_plot=True):
 
         self.C_m  = C_m
@@ -79,14 +79,14 @@ def __init__(self, C_m=1, gmax_Na=120, gmax_K=36, gmax_L=0.3, E_Na=50,
         # plotting conditionals
         self.injected_current_plot = injected_current_plot
         self.gating_plot = gating_plot
-        self.scaling_plot = scaling_plot
+        self.cond_scaling_plot = cond_scaling_plot
         self.cond_dens_plot = cond_dens_plot
         self.driving_force_plot = driving_force_plot
         self.current_plot = current_plot
         self.memb_pot_plot = memb_pot_plot
 
         self.num_plots = (int(self.injected_current_plot) +
-                          int(self.gating_plot)+ int(self.scaling_plot) +
+                          int(self.gating_plot)+ int(self.cond_scaling_plot) +
                           int(self.cond_dens_plot) + int(self.driving_force_plot) +
                           int(self.current_plot) + int(self.memb_pot_plot))
 
@@ -321,7 +321,7 @@ def Main(self, init_values=[-64.99584, 0.05296, 0.59590, 0.31773]):
                 plt.legend()
                 self.plot_count += 1
 
-            if self.scaling_plot:
+            if self.cond_scaling_plot:
                 try:
                     plt.subplot(self.num_plots,1,self.plot_count+1, sharex = ax1)
                 except NameError:
@@ -364,8 +364,8 @@ def Main(self, init_values=[-64.99584, 0.05296, 0.59590, 0.31773]):
                 ax_here.fill_between(self.t, dna, color='c', alpha=0.5)
                 ax_here.fill_between(self.t, dk, color='y', alpha=0.5)
 
-                plt.plot(self.t, dna, 'c', label='$V - E_{Na}$', linewidth=0.8)
-                plt.plot(self.t, dk, 'y', label='$V - E_{K}$', linewidth=0.8)
+                plt.plot(self.t, dna, 'c', label='$V_{m} - E_{Na}$', linewidth=0.8)
+                plt.plot(self.t, dk, 'y', label='$V_{m} - E_{K}$', linewidth=0.8)
                 plt.plot(self.t, zero, 'k', linestyle='dashed', linewidth=0.5)
                 plt.ylabel('Driving force (mV)')
                 plt.legend()

notebooks/Python_HH_version/Python_Notebook_HH.ipynb

@@ -37,7 +37,7 @@
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "012cf4db99d64e63b2aba86949411ced",
+       "model_id": "c131ca61e49a441e95dff9c74a0b29f9",
        "version_major": 2,
        "version_minor": 0
       },
@@ -51,7 +51,7 @@
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "76a0b6d59696460dbbf15f755c87eee9",
+       "model_id": "3080e0df7a9645299d1e6f2f7992ef53",
        "version_major": 2,
        "version_minor": 0
       },
@@ -88,11 +88,15 @@
     "\n",
     "2) The second plot shows the **activation/inactivation variables** of the ion channels in the neuron: activation variable **m** and inactivation variable **h** for *Na*, and activation variable **n** for *K*. Activation variables increase as the cell becomes depolarised (**V<sub>m</sub>** increases) and the inactivation variable decreases. **m** responds to these changes quicker than **h** or **n**.\n",
     "\n",
-    "3) The third plot shows the **conductance densities** for Sodium (**g<sub>Na</sub> = gmax<sub>Na</sub> m<sup>3</sup> h**) and Potassium (**g<sub>K</sub> = gmax<sub>K</sub> n<sup>4</sup>**). Note, the leak conductance density value, **g<sub>L</sub>** is not voltage dependent, and so does not vary with time. The maximum conductance densities for the 3 conductances can be changed or set to zero.\n",
+    "3) The third (optional) plot shows the **conductance scaling factors** for Sodium (**m<sup>3</sup> * h**) and Potassium (**n<sup>4</sup>**). Note: these only ever have a value between 0 and 1. \n",
     "\n",
-    "4) The fourth plot shows the influx (negative y-axis) and outflux (positive y-axis) of **ionic currents (I<sub>Na</sub> = g<sub>Na</sub> (V<sub>m</sub> - E<sub>Na</sub>), I<sub>K</sub> = g<sub>K</sub> (V<sub>m</sub> - E<sub>K</sub>) and I<sub>L</sub> = g<sub>L</sub> (V<sub>m</sub> - E<sub>L</sub>**) passing through each type of ion channel being modeled. The reversal potentials of each of these (**E<sub>X</sub>**) can be altered above. \n",
+    "4) The fourth plot shows the **conductance densities** for Sodium (**g<sub>Na</sub> = gmax<sub>Na</sub> m<sup>3</sup> h**) and Potassium (**g<sub>K</sub> = gmax<sub>K</sub> n<sup>4</sup>**). Note, the leak conductance density value, **g<sub>L</sub>** is not voltage dependent, and so does not vary with time. The maximum conductance densities for the 3 conductances can be changed or set to zero.\n",
     "\n",
-    "5) The bottom plot shows the **neural membrane potential activity, V<sub>m</sub>**. The spikes here are called \"action potentials\". This variable is governed by the expression: **dV<sub>m</sub>/dt = (I<sub>inj</sub> - I<sub>Na</sub> - I<sub>K</sub> - I<sub>L</sub>)/C<sub>m</sub>**\n",
+    "5) The fifth (optional) plot shows the **driving force** for the Sodium (**V<sub>m</sub> - E<sub>Na</sub>**) and Potassium (**V<sub>m</sub> - E<sub>K</sub>**) currents. Note, the important value here is how different these are from the zero line, i.e. the magnitude of the force driving ions in or out of the cell. \n",
+    "\n",
+    "6) The sixth plot shows the influx (negative y-axis) and outflux (positive y-axis) of **ionic currents** (**I<sub>Na</sub> = g<sub>Na</sub> (V<sub>m</sub> - E<sub>Na</sub>), I<sub>K</sub> = g<sub>K</sub> (V<sub>m</sub> - E<sub>K</sub>)** and **I<sub>L</sub> = g<sub>L</sub> (V<sub>m</sub> - E<sub>L</sub>)**) passing through each type of ion channel being modeled. The reversal potentials of each of these (**E<sub>X</sub>**) can be altered above. \n",
+    "\n",
+    "7) The bottom plot shows the **neural membrane potential activity, V<sub>m</sub>**. The spikes here are called \"action potentials\". This variable is governed by the expression: **dV<sub>m</sub>/dt = (I<sub>inj</sub> - I<sub>Na</sub> - I<sub>K</sub> - I<sub>L</sub>)/C<sub>m</sub>**\n",
     "\n",
     "\n",
     "![hh](../../Tutorial/_media/equivalentCircuit.png)\n"
@@ -104,7 +108,9 @@
    "source": [
     "<font size=\"5\">Notes</font>\n",
     "\n",
-    ">**Note: Sign of channel current $I_{Na}$, $I_{K}$ and $I_{L}$ are opposite to NeuroML tutorial<br />**\n",
+    ">**Note 1: Python source of the HH model<br />** The source of the HH model in Python as used in this notebook can be found [here](https://github.com/openworm/hodgkin_huxley_tutorial/blob/master/Tutorial/Source/HodgkinHuxley.py).\n",
+    "\n",
+    ">**Note 2: Sign of channel current $I_{Na}$, $I_{K}$ and $I_{L}$ are opposite to NeuroML tutorial<br />**\n",
     "The Python code is interested in the current coming through the channel (so positive if flowing \"out\" like K, negative if flowing \"in\", like Na), whereas for the NeuroML channels (more specifically the channel density element) is what's flowing into the cell (so positive for Na)."
    ]
   },

notebooks/Python_HH_version/ui_widget.py

@@ -206,13 +206,15 @@ def highlight_slider():
 button_row=ipywidgets.HBox([reset_button,showValue_togglebtn])
 
 #plot selectors
-header_plotting = ipywidgets.HTMLMath(value=r"<b> Select plots</b>")
-injected_current_plot_value = ipywidgets.Checkbox(value=True, description="Current injection", disabled=False,)
-gating_plot_value = ipywidgets.Checkbox(value=True, description="Gating variables", disabled=False,)
-cond_dens_plot_value = ipywidgets.Checkbox(value=True, description="Conductance densities", disabled=False,)
-current_plot_value = ipywidgets.Checkbox(value=True, description="Current densities", disabled=False,)
-memb_pot_plot_value = ipywidgets.Checkbox(value=True, description="Membrane potential", disabled=False,)
-plot_selection_row=ipywidgets.VBox([header_plotting, ipywidgets.HBox([injected_current_plot_value, gating_plot_value, cond_dens_plot_value]), ipywidgets.HBox([current_plot_value, memb_pot_plot_value])])
+header_plotting = ipywidgets.HTMLMath(value=r"<b> Select plots to show</b>")
+injected_current_plot_value = ipywidgets.Checkbox(value=True, description="1) Current injection", disabled=False,)
+gating_plot_value = ipywidgets.Checkbox(value=True, description="2) Gating variables", disabled=False,)
+cond_scaling_value = ipywidgets.Checkbox(value=False, description="3) Conductance scaling", disabled=False,)
+cond_dens_plot_value = ipywidgets.Checkbox(value=True, description="4) Conductance densities", disabled=False,)
+driving_force_value = ipywidgets.Checkbox(value=False, description="5) Driving force", disabled=False,)
+current_plot_value = ipywidgets.Checkbox(value=True, description="6) Current densities", disabled=False,)
+memb_pot_plot_value = ipywidgets.Checkbox(value=True, description="7) Membrane potential", disabled=False,)
+plot_selection_row=ipywidgets.VBox([header_plotting, ipywidgets.HBox([injected_current_plot_value, gating_plot_value, cond_scaling_value, cond_dens_plot_value]), ipywidgets.HBox([driving_force_value, current_plot_value, memb_pot_plot_value])])
 
 #layout vertically all the widgets defined above
 modelInputs=ipywidgets.VBox([h1,h2,h3,h4,h5,h6,h7,h8,h9,h10,runMode_iclamp,runMode_vclamp,button_row,defalultValues,plot_selection_row])
@@ -230,7 +232,7 @@ def launch_interactive_widget():
     def runHH(C_m, g_Na, g_K, g_L, E_Na, E_K, E_L, t_0, t_n, delta_t,
               I_inj_max, I_inj_width, I_inj_trans, vc_delay, vc_duration,
               vc_condVoltage, vc_testVoltage, vc_returnVoltage, runMode,
-              injected_current_plot, gating_plot, cond_dens_plot, current_plot, memb_pot_plot):
+              injected_current_plot, gating_plot, cond_scaling_plot, cond_dens_plot, driving_force_plot, current_plot, memb_pot_plot):
 
         highlight_slider()
         runner = HHmodel.HodgkinHuxley(C_m, g_Na, g_K, g_L, E_Na, E_K, E_L,
@@ -241,7 +243,9 @@ def runHH(C_m, g_Na, g_K, g_L, E_Na, E_K, E_L, t_0, t_n, delta_t,
                                        runMode,
                                        injected_current_plot=injected_current_plot,
                                        gating_plot=gating_plot,
+                                       cond_scaling_plot=cond_scaling_plot,
                                        cond_dens_plot=cond_dens_plot,
+                                       driving_force_plot=driving_force_plot,
                                        current_plot=current_plot,
                                        memb_pot_plot=memb_pot_plot)
         # init_values are the steady state values for v,m,h,n at zero current injection
@@ -257,7 +261,9 @@ def runHH(C_m, g_Na, g_K, g_L, E_Na, E_K, E_L, t_0, t_n, delta_t,
                                             'vc_testVoltage':textBox_testVoltage,'vc_returnVoltage':textBox_returnVoltage,
                                             'runMode':runMode_togglebtns, 'injected_current_plot': injected_current_plot_value,
                                                       'gating_plot':gating_plot_value,
+                                                      'cond_scaling_plot':cond_scaling_value,
                                                       'cond_dens_plot':cond_dens_plot_value,
+                                                      'driving_force_plot':driving_force_value,
                                                       'current_plot':current_plot_value,
                                                       'memb_pot_plot':memb_pot_plot_value})