Package: currentscape
CurrentscapeLatest ReleaseDocumentationLicenseBuild StatusCoverageGitterCitationIntroductionCurrentscape is a Python tool enabling scientists to easily plot the currents in electrical neuron models.
The code is based on the paperAlonso 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.CitationWhen 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}
}SupportWe are providing support onGitter. We suggest you create tickets on theGithub issue trackerin case you encounter problems while using the software or if you have some suggestions.Main dependenciesPython 3.8+Numpy(automatically installed by pip)Palettable(automatically installed by pip)InstallationCurrentscape can be pip installed with the following command:pipinstallcurrentscapeIf you want to be able to run the Currentscape examples, you will need to also install the example dependencies:pipinstallcurrentscape[example]Quick StartBelow 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:pipinstallneuronWhen you then execute the following python code, a window should open with the currentscape plot:importnumpyasnpfromneuronimporthfromneuron.unitsimportms,mVimportcurrentscapedefmain():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()defrun_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.6157dend.L=200dend.diam=1forsecinh.allsec():sec.Ra=100# Axial resistance in Ohm * cmsec.cm=1# Membrane capacitance in micro Farads / cm^2# Insert active Hodgkin-Huxley current in the somasoma.insert('hh')forseginsoma:seg.hh.gnabar=0.12# Sodium conductance in S/cm2seg.hh.gkbar=0.036# Potassium conductance in S/cm2seg.hh.gl=0.0003# Leak conductance in S/cm2seg.hh.el=-54.3# Reversal potential in mV# Insert passive current in the dendritedend.insert('pas')forsegindend:seg.pas.g=0.001# Passive conductance in S/cm2seg.pas.e=-65# Leak reversal potential mVstim=h.IClamp(dend(1))stim.delay=5stim.dur=10stim.amp=0.1t_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 pAvoltage=np.asarray(v_vec)potassium=np.asarray(ik_vec)*to_pAsodium=np.asarray(ina_vec)*to_pAleak=np.asarray(il_vec)*to_pAreturnvoltage,potassium,sodium,leakif__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):TutorialA more detailed explanation on how to use Currentscape, as well as other examples can be found on thetutorial page.API DocumentationThe API documentation can be found onReadTheDocs.Funding & AcknowledgementsWe wish to thank the authors ofAlonso and Marder, 2019to let us integrate a part of theircodeinto 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.