Activity-dependent neuromodulation and calcium homeostasis cooperate to produce robust and modulable neuronal function.
The 5 matches
- [1] § Results › Sharp neuromodulation interferes with calcium homeostasis whereas controlled neuromodulation naturally leads to robust and modulable neuronal function ↔ fig2_controlled/STG_utils.jl, lines 45–128 · score 0.73 · interspike interval, interburst frequency, intraburst frequency, spike bursts, burstiness, firing
- [2] § Materials and methods › The neuromodulation controller ↔ fig6/network_STG_DIC.jl, lines 323–386 · score 0.72 · sensitivity matrix, linear system, threshold voltage, gating variables, DICs, Vth
- [3] § Materials and methods › The neuromodulation controller ↔ fig2_controlled/DA_ODE.jl, lines 1–51 · score 0.62 · sensitivity matrix, threshold voltage, gating variables, Vth, model
- [4] § Materials and methods › Conductance-based model ↔ fig2_controlled/STG_models.jl, lines 1–54 · score 0.58 · delayed rectifier potassium, STG model, slow calcium, activated, sodium
- [5] § Materials and methods › Conductance-based model ↔ fig6/network_STG_animation.jl, lines 1–39 · score 0.53 · reversal potentials, membrane capacitance, synaptic, gated, variables, STG
Paper
Loaded from Europe PMC by your browser, not stored by OSCR: doi.org · Europe PMC
The paper is loaded when this pane is shown.
The authors' code
Julia · 128 lines · 3.5 KB · MIT · 1 match
- #=
- This file contains functions to extract characteristics of the firing pattern
- as well as some functions to plot complicated graphs
- =#
- using Statistics, Plots, StatsPlots, LaTeXStrings, Printf
- ## Functions extracting characteristics of the firing pattern
- # This function extracts the spiking frequency of a spiking firing pattern
- function extract_frequency(V, t)
- # Defining thresholds
- spike_up_threshold = 10.
- spike_down_threshold = 0.
- # Detecting spikes
- spike_detected = 0
- spike_times = []
- for i in 1:length(V)
- if V[i] > spike_up_threshold && spike_detected == 0 # Start of spike
- append!(spike_times, t[i])
- spike_detected = 1
- end
- if V[i] < spike_down_threshold && spike_detected == 1 # End of spike
- spike_detected = 0
- end
- end
- # If the neuron is silent
- if length(spike_times) < 2
- return NaN
- end
- # Calculating all interspike intervals
- ISI=[]
- for i in 2 : length(spike_times)
- append!(ISI, spike_times[i] - spike_times[i-1])
- end
- # If the neuron is silent
- if length(ISI) < 2
- return NaN
- end
- # Computing the spiking frequency
- T = mean(ISI) / 1000 # in seconds
- f = 1 / T # in Hz
- return f
- end
- # This function extracts characteristics of a bursting firing pattern
- function extract_burstiness(V, t)
- # Defining thresholds
- spike_up_threshold = 10.
- spike_down_threshold = 0.
- # Detecting spikes
- spike_detected = 0
- spike_times = []
- for i in 1 : length(V)
- if V[i] > spike_up_threshold && spike_detected == 0 # Start of spike
- append!(spike_times, t[i])
- spike_detected = 1
- end
- if V[i] < spike_down_threshold && spike_detected == 1 # End of spike
- spike_detected = 0
- end
- end
- # If the neuron is silent
- if length(spike_times) < 3
- return NaN, NaN, NaN, NaN
- end
- # Calculating all interspike intervals
- ISI = []
- for i in 2 : length(spike_times)
- append!(ISI, spike_times[i] - spike_times[i-1])
- end
- # Defining a threshold to separate intraburst from interburst ISI
- max_ISI = maximum(ISI)
- min_ISI = minimum(ISI)
- half_ISI = (max_ISI+min_ISI)/2
- # If ISI too constant, neuron is spiking
- if max_ISI - min_ISI < 25
- return NaN, NaN, NaN, NaN
- end
- # Detecting the first spike of a burst
- first_spike_burst = findall(x -> x > half_ISI, ISI)
- # Computing the interburst frequency
- Ts = ISI[first_spike_burst]
- interburst_T = mean(Ts) / 1000 # in seconds
- interburst_f = 1 / interburst_T # in Hz
- # Computing the number of spikes per burst
- nb_spike_burst = []
- for i in 2 : length(first_spike_burst)
- append!(nb_spike_burst, first_spike_burst[i] - first_spike_burst[i-1])
- end
- # If spiking
- if length(nb_spike_burst) < 2
- return NaN, NaN, NaN, NaN
- end
- nb_spike_per_burst = round(mean(nb_spike_burst))
- # If no bursting
- if nb_spike_per_burst < 1.5 || nb_spike_per_burst > 500
- burstiness = NaN
- intraburst_f = NaN
- nb_spike_per_burst = NaN
- interburst_f = NaN
- else # Else, bursting: computing the intraburst frequency
- intra_spike_burst = findall(x -> x < half_ISI, ISI)
- Ts_intraburst = ISI[intra_spike_burst]
- T_intraburst = mean(Ts_intraburst) / 1000 # in seconds
- intraburst_f = 1 / T_intraburst # in Hz
- burstiness = (nb_spike_per_burst * intraburst_f) / interburst_T
- end
- return burstiness, nb_spike_per_burst, intraburst_f, interburst_f
- end
STG_utils.jl at commit 4b83c2d, under MIT · at the source
Overview
Abstract
Neurons rely on two interdependent mechanisms — homeostasis and neuromodulation — to maintain robust and adaptable functionality. Calcium homeostasis stabilizes neuronal activity by adjusting ionic conductances, whereas neuromodulation dynamically modifies ionic properties in response to external signals carried by neuromodulators. Combining these mechanisms in conductance-based models often produces unreliable outcomes, particularly when sharp neuromodulation interferes with calcium-homeostatic tuning. This study explores how a biologically inspired neuromodulation controller can harmonize with calcium homeostasis to ensure reliable neuronal function. Using computational models of stomatogastric ganglion and dopaminergic neurons, we demonstrate that controlled neuromodulation preserves neuronal firing patterns while calcium homeostasis simultaneously maintains target intracellular calcium levels. Unlike sharp neuromodulation, the neuromodulation controller integrates activity-dependent feedback through mechanisms mimicking G-protein-coupled receptor cascades. The interaction between these controllers critically depends on the existence of an intersection in conductance space, representing a balance between target calcium levels and neuromodulated firing patterns. Maximizing neuronal degeneracy enhances the likelihood of such intersections, enabling robust modulation and compensation for channel blockades. We further show that this controller pairing extends to network-level activity, reliably modulating the rhythmic activity of central pattern generators. This study highlights the complementary roles of calcium homeostasis and neuromodulation, proposing a unified control framework for maintaining robust and adaptive neural activity under physiological and pathological conditions.
Reproduced under the paper's license (CC BY), from the paper cited above.
Repositories
Its files are read in the Code ↔ Paper reader above, with 5 matches between paragraphs and lines of code.
arthur-fyon/BIOCONTROL_2025
4b83c2ddce05461bdadca66c1bb6352527ebf8fe, 18 May 2026Availability: 1 check, the latest on 29 September 2026: the link answers
- 29 September 2026: the link answers
69 files
- dependencies.jl — Julia, 16 lines
- fig2_controlled/
DA_ODE.jl — Julia, 163 lines, 1 match - fig2_controlled/
DA_kinetics.jl — Julia, 42 lines - fig2_controlled/
DA_large_simulations.jl — Julia, 87 lines - fig2_controlled/
DA_neuromod_CaLCaN_large — Jupyter, 686 linessimu.ipynb - fig2_controlled/
STG_DIC.jl — Julia, 302 lines - fig2_controlled/
STG_gs_derivatives.jl — Julia, 40 lines - fig2_controlled/
STG_kinetics.jl — Julia, 43 lines - fig2_controlled/
STG_large_simulations.jl — Julia, 175 lines - fig2_controlled/
STG_models.jl — Julia, 302 lines, 1 match - fig2_controlled/
STG_neuromod_CaSA_larges — Jupyter, 333 linesimu.ipynb - fig2_controlled/
STG_neuromod_CaSA_larges — Jupyter, 594 linesimu_washout.ipynb - fig2_controlled/
STG_neuromod_CaSH_larges — Jupyter, 436 linesimu.ipynb - fig2_controlled/
STG_neuromodulation.jl — Julia, 73 lines - fig2_controlled/
STG_utils.jl — Julia, 128 lines, 1 match - fig2_sharp/
STG_DIC.jl — Julia, 302 lines - fig2_sharp/
STG_gs_derivatives.jl — Julia, 40 lines - fig2_sharp/
STG_kinetics.jl — Julia, 43 lines - fig2_sharp/
STG_large_simulations.jl — Julia, 469 lines - fig2_sharp/
STG_models.jl — Julia, 647 lines - fig2_sharp/
STG_neuromod_CaSA_larges — Jupyter, 318 linesimu.ipynb - fig2_sharp/
STG_neuromod_CaSA_larges — Jupyter, 516 linesimu_H.ipynb - fig2_sharp/
STG_neuromod_CaSA_larges — Jupyter, 334 linesimu_single_traces.ipynb - fig2_sharp/
STG_neuromod_CaSA_larges — Jupyter, 523 linesimu_washout.ipynb - fig2_sharp/
STG_neuromod_CaSA_other_ — Jupyter, 552 linestraces.ipynb - fig2_sharp/
STG_neuromod_CaSA_voltag — Jupyter, 410 linese_traces.ipynb - fig2_sharp/
STG_neuromodulation.jl — Julia, 105 lines - fig2_sharp/
STG_utils.jl — Julia, 128 lines - fig3_controlled/
STG_DIC.jl — Julia, 247 lines - fig3_controlled/
STG_gs_derivatives.jl — Julia, 40 lines - fig3_controlled/
STG_kinetics.jl — Julia, 43 lines - fig3_controlled/
STG_large_simulations.jl — Julia, 95 lines - fig3_controlled/
STG_models.jl — Julia, 131 lines - fig3_controlled/
STG_neuromod_CaSA_larges — Jupyter, 558 linesimu.ipynb - fig3_controlled/
STG_neuromod_CaSA_larges — Jupyter, 352 linesimu_single_traces.ipynb - fig3_controlled/
STG_neuromodulation.jl — Julia, 73 lines - fig3_controlled/
STG_utils.jl — Julia, 128 lines - fig3_sharp/
STG_DIC.jl — Julia, 247 lines - fig3_sharp/
STG_gs_derivatives.jl — Julia, 40 lines - fig3_sharp/
STG_kinetics.jl — Julia, 43 lines - fig3_sharp/
STG_large_simulations.jl — Julia, 144 lines - fig3_sharp/
STG_models.jl — Julia, 95 lines - fig3_sharp/
STG_neuromod_CaSA_larges — Jupyter, 344 linesimu.ipynb - fig3_sharp/
STG_neuromodulation.jl — Julia, 73 lines - fig3_sharp/
STG_utils.jl — Julia, 128 lines - fig5/
STG_DIC.jl — Julia, 247 lines - fig5/
STG_gs_derivatives.jl — Julia, 40 lines - fig5/
STG_kinetics.jl — Julia, 43 lines - fig5/
STG_large_simulations.jl — Julia, 420 lines - fig5/
STG_models.jl — Julia, 607 lines - fig5/
STG_neuromod_CaSA_larges — Jupyter, 248 linesimu.ipynb - fig5/
STG_neuromod_CaSA_larges — Jupyter, 413 linesimu_plots.ipynb - fig5/
STG_neuromod_CaSA_neuron — Jupyter, 549 lines1.ipynb - fig5/
STG_neuromod_CaSA_neuron — Jupyter, 521 lines2.ipynb - fig5/
STG_neuromod_CaSA_neuron — Jupyter, 514 lines3.ipynb - fig5/
STG_neuromodulation.jl — Julia, 73 lines - fig5/
STG_utils.jl — Julia, 128 lines - fig6/
baseline_activity.ipynb — Jupyter, 312 lines - fig6/
network_STG_DIC.jl — Julia, 386 lines, 1 match - fig6/
network_STG_animation.jl — Julia, 347 lines, 1 match - fig6/
network_STG_gs_derivativ — Julia, 40 lineses.jl - fig6/
network_STG_kinetics.jl — Julia, 54 lines - fig6/
network_STG_models.jl — Julia, 368 lines - fig6/
network_STG_neuromodulat — Julia, 73 linesion.jl - fig6/
network_STG_utils.jl — Julia, 128 lines - fig6/
network_simulation_contr — Jupyter, 286 linesolled.ipynb - fig6/
network_simulation_sharp — Jupyter, 326 lines.ipynb - LICENSE — License, 21 lines
- README.md — Text, 49 lines
Zenodo 17827990
Availability: 1 check, the latest on 29 September 2026: the link answers (HTTP 200)
- 29 September 2026: the link answers (HTTP 200)
The paper's code and data availability statement is in the Data section.
Tracing map
Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.
What the map holds:
- 2 repositories of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
- 67 scripts, each with its path and the digest of its content;
- 5 matches between paragraphs of the paper and lines of the code (method lexical-v1);
- neither the text of the paper nor the code itself.
Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.
Data
No dataset and no data link were found in the paper.
Data Availability
There are no primary data in the paper; all materials are available at https://
Reproduced under the paper's license (CC BY), from the paper cited above.
Versions
The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.
Version 1, 29 September 2026: the first record
Recorded: type, language, journal, volume, issue, pages, dates, 2 authors, 10 MeSH terms, 2 funders, 67 references.
Cite
This paper
Fyon, A., & Drion, G. (2026). Activity-dependent neuromodulation and calcium homeostasis cooperate to produce robust and modulable neuronal function. PLoS computational biology, 22(4), e1014177. https://
BibTeX
@article{fyon2026activit
author = {Fyon, Arthur and Drion, Guillaume},
title = {{Activity-dependent neuromodulation and calcium homeostasis cooperate to produce robust and modulable neuronal function}},
journal = {PLoS computational biology},
year = {2026},
month = apr,
volume = {22},
number = {4},
pages = {e1014177},
publisher = {PLOS},
issn = {1553-734X},
doi = {10.1371/
url = {https://
pmid = {41984991},
pmcid = {PMC13102308}
}
RIS
TY - JOUR
AU - Fyon, Arthur
AU - Drion, Guillaume
TI - Activity-dependent neuromodulation and calcium homeostasis cooperate to produce robust and modulable neuronal function
T2 - PLoS computational biology
J2 - PLoS Comput Biol
PY - 2026
DA - 2026/
VL - 22
IS - 4
SP - e1014177
SN - 1553-734X
PB - PLOS
DO - 10.1371/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1371/
"type": "article-journal",
"title": "Activity-dependent neuromodulation and calcium homeostasis cooperate to produce robust and modulable neuronal function",
"container-title": "PLoS computational biology",
"author": [
{
"family": "Fyon",
"given": "Arthur"
},
{
"family": "Drion",
"given": "Guillaume"
}
],
"container-title-short":
"volume": "22",
"issue": "4",
"page": "e1014177",
"DOI": "10.1371/
"PMID": "41984991",
"PMCID": "PMC13102308",
"ISSN": "1553-734X",
"publisher": "PLOS",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
4,
15
]
]
}
}
The tracing map gets a citation of its own once an author has validated it and it has a DOI.
Similar papers
The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.
- [1] doi:10.1371/journal.pcbi.1014337 [code]
- Fast reconstruction of degenerate populations of conductance-based neuron models from spike times.Journal: PLoS computational biologyIn common: computational modeling (no new data), 13 references, author Arthur Fyon
- [2] doi:10.1007/s10827-026-00936-7 [code]
- When can neuronal activity-dependent homeostatic plasticity maintain circuit-level properties?Journal: Journal of computational neuroscienceIn common: none (in silico), computational modeling (no new data), 13 references
- [3] doi:10.1038/s41540-026-00749-5 [code]
- A novel approach to quantify out-of-distribution uncertainty in Neural and Universal Differential Equations.Journal: NPJ systems biology and applicationsIn common: Makie, DifferentialEquations.jl, Plots.jl, computational modeling (no new data), 1 reference
- [4] doi:10.1093/pnasnexus/pgag213 [code]
- Two-factor synaptic plasticity enables memory consolidation during neuronal burst firing.Journal: PNAS nexusIn common: DifferentialEquations.jl, Plots.jl, 3 references
- [5] doi:10.1371/journal.pcbi.1014458 [code]
- Neuronal excitability and parameter variability in the Hodgkin-Huxley model.Journal: PLoS computational biologyIn common: none (in silico), computational modeling (no new data), 5 references
- [6] doi:10.7554/elife.89629 [code]
- Active dendrites enable robust spiking computations despite timing jitter.Journal: eLifeIn common: Makie, DifferentialEquations.jl, Plots.jl
- [7] doi:10.1523/eneuro.0029-26.2026 [code]
- Electrical Coupling within Thalamocortical Networks Cumulatively Reduces Cortical Correlation to Sensory Inputs.Journal: eNeuroIn common: Makie, DifferentialEquations.jl, Plots.jl
- [8] doi:10.1162/imag.a.1249 [code]
- Global search metaheuristics for neural mass model calibration.Journal: Imaging neuroscience (Cambridge, Mass.)In common: DifferentialEquations.jl, Plots.jl, none (in silico), 1 reference
- [9] doi:10.1371/journal.pcbi.1014549 [code]
- Inward rectifier potassium channels interact with calcium channels to promote robust and physiological bistability.Journal: PLoS computational biologyIn common: DifferentialEquations.jl, Plots.jl, 1 reference
- [10] doi:10.7554/elife.108352 [code]
- Analysis of dendritic input currents during place field dynamics.Journal: eLifeIn common: 4 references
Contribute
The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.
Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.
Claim this paper
Correct its record
Say what each link of this record is, remove the ones that are not the paper's, add the ones that are missing. The correction becomes a new version of the record, in its Versions section.
Validate its tracing map
You validate the map as this page shows it: 2 repositories of the authors' code, each at its verified commit and with its license, 67 scripts, and 5 matches between paragraphs and code (see the Code and Map sections). It then receives a DOI on Zenodo, with you (your ORCID iD) and OSCR as its creators; the code itself is not deposited.
The map's fingerprint: sha256:dba8ec2bc304d6c6…
Add the badge to its README
The badge links the code to this page. Copy one of these into the README of the paper's code: only you decide where it goes, and nothing is changed for you.
Markdown
[, paste the snippet at the top, then “Commit changes…” and, to review it first, “Create a new branch and start a pull request”. You open the pull request; OSCR asks for no permission.
Request its removal
To ask OSCR to remove this record, the copies of its authors' scripts or its tracing map, use the removal request page: signed in, you say who you are, what to remove and why, then review and confirm the request. Published rules decide every request (how).
Discussion, reproductions, activity
Discussion: questions and error reports about this paper and its code, from signed-in readers and its authors. It opens with sign-in.
Reproductions: reports from readers who ran the authors' code: what they reproduced, with which environment, commit and data. It opens with sign-in.
Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.
