| Latest Release | 
|
| Documentation | |
| License | 
|
| Build Status |  |
| Coverage |  |
| Gitter | 
|
| Citation | 
|
Currentscape is a Python tool enabling scientists to easily plot the currents in electrical neuron models. The code is based on the paper Alonso and Marder, 2019.
Currentscape figures plot the percentage of inward and outward ionic membrane currents, the total inward and outward currents, as well as the voltage in function of time. It allows modellers to see which currents play a role at any given time during a simulation, and check in depth the current dynamics.
When you use this Currentscape software for your research, we ask you to cite the following publications (this includes poster presentations):
@article {alonsomarder2019,
article_type = {journal},
title = {Visualization of currents in neural models with similar behavior and different conductance densities},
author = {Alonso, Leandro M and Marder, Eve},
editor = {Westbrook, Gary L and Skinner, Frances K and Lankarany, Milad and Britton, Oliver},
volume = 8,
year = 2019,
month = {jan},
pub_date = {2019-01-31},
pages = {e42722},
citation = {eLife 2019;8:e42722},
doi = {10.7554/eLife.42722},
url = {https://doi.org/10.7554/eLife.42722},
abstract = {Conductance-based models of neural activity produce large amounts of data that can be hard to visualize and interpret. We introduce visualization methods to display the dynamics of the ionic currents and to display the models’ response to perturbations. To visualize the currents’ dynamics, we compute the percent contribution of each current and display them over time using stacked-area plots. The waveform of the membrane potential and the contribution of each current change as the models are perturbed. To represent these changes over a range of the perturbation control parameter, we compute and display the distributions of these waveforms. We illustrate these procedures in six examples of bursting model neurons with similar activity but that differ as much as threefold in their conductance densities. These visualization methods provide heuristic insight into why individual neurons or networks with similar behavior can respond widely differently to perturbations.},
keywords = {neuronal oscillators, Na+ channels, Ca++ channels, K+ channels, conductance-based, ionic channels},
journal = {eLife},
issn = {2050-084X},
publisher = {eLife Sciences Publications, Ltd},
}
@article{currentscape,
title={Currentscape},
DOI={10.5281/zenodo.8046373},
abstractNote={Currentscape is a Python tool enabling scientists to easily plot the currents in electrical neuron models. The code is based on the paper Alonso and Marder, 2019. Currentscape figures plot the percentage of inward and outward ionic membrane currents, the total inward and outward currents, as well as the voltage in function of time. It allows modellers to see which currents play a role at any given time during a simulation, and check in depth the current dynamics.},
publisher={Zenodo},
author={Jaquier, Aurélien and Tuncel, Anil and Van Geit, Werner and Alonso, Leandro M and Marder, Eve},
year={2023},
month={Jun}
}
We are providing support on Gitter. We suggest you create tickets on the Github issue tracker in case you encounter problems while using the software or if you have some suggestions.
- Python 3.8+
- Numpy (automatically installed by pip)
- Palettable (automatically installed by pip)
Currentscape can be pip installed with the following command:
pip install currentscapeIf you want to be able to run the Currentscape examples, you will need to also install the example dependencies:
pip install currentscape[example]Below is an example of a ball and stick model in NEURON with simple Hodgkin-Huxley mechanisms, to which a step stimulus is applied.
The voltage and ionic currents are recorded and fed to Currentscape, along with a configuration dictionary containing the current names to be displayed in the legend.
To run the code you will first have to install NEURON package:
pip install neuronWhen you then execute the following python code, a window should open with the currentscape plot:
import numpy as np
from neuron import h
from neuron.units import ms, mV
import currentscape
def main():
current_names = ["Potassium", "Sodium", "Leak"]
voltage, potassium, sodium, leak = run_sim()
config = {
"output": {
"savefig": True,
"dir": ".",
"fname": "quickstart_plot",
"extension": "png",
"dpi": 300,
"transparent": False
},
"current": {"names": current_names},
"voltage": {"ylim": [-90, 50]},
"legendtextsize": 5,
"adjust": {
"left": 0.15,
"right": 0.8,
"top": 1.0,
"bottom": 0.0
}
}
fig = currentscape.plot(voltage, [potassium, sodium, leak], config)
fig.show()
def run_sim():
h.load_file('stdrun.hoc')
soma = h.Section(name='soma')
dend = h.Section(name='dend')
dend.connect(soma(1))
soma.L = soma.diam = 12.6157
dend.L = 200
dend.diam = 1
for sec in h.allsec():
sec.Ra = 100 # Axial resistance in Ohm * cm
sec.cm = 1 # Membrane capacitance in micro Farads / cm^2
# Insert active Hodgkin-Huxley current in the soma
soma.insert('hh')
for seg in soma:
seg.hh.gnabar = 0.12 # Sodium conductance in S/cm2
seg.hh.gkbar = 0.036 # Potassium conductance in S/cm2
seg.hh.gl = 0.0003 # Leak conductance in S/cm2
seg.hh.el = -54.3 # Reversal potential in mV
# Insert passive current in the dendrite
dend.insert('pas')
for seg in dend:
seg.pas.g = 0.001 # Passive conductance in S/cm2
seg.pas.e = -65 # Leak reversal potential mV
stim = h.IClamp(dend(1))
stim.delay = 5
stim.dur = 10
stim.amp = 0.1
t_vec = h.Vector()
v_vec = h.Vector()
ik_vec = h.Vector()
ina_vec = h.Vector()
il_vec = h.Vector()
t_vec.record(h._ref_t)
v_vec.record(soma(0.5)._ref_v)
ik_vec.record(soma(0.5)._ref_ik)
ina_vec.record(soma(0.5)._ref_ina)
il_vec.record(soma(0.5)._ref_il_hh)
h.finitialize(-65 * mV)
h.continuerun(25 * ms)
to_pA = 10 * soma(0.5).area() # turn mA/cm2 (*um2) into pA
voltage = np.asarray(v_vec)
potassium = np.asarray(ik_vec) * to_pA
sodium = np.asarray(ina_vec) * to_pA
leak = np.asarray(il_vec) * to_pA
return voltage, potassium, sodium, leak
if __name__ == "__main__":
main()When you run this code in Python, it will generate the following currentscape plot (in a window, and on disk as quickstart_plot.png):
A more detailed explanation on how to use Currentscape, as well as other examples can be found on the tutorial page.
The API documentation can be found on ReadTheDocs.
We wish to thank the authors of Alonso and Marder, 2019 to let us integrate a part of their code into this repository.
The part of the code in this repository developed by the EPFL Blue Brain Project was supported by funding to the Blue Brain Project, a research center of the École polytechnique fédérale de Lausanne (EPFL), from the Swiss government's ETH Board of the Swiss Federal Institutes of Technology.
Copyright (c) 2023-2024 Blue Brain Project/EPFL