OSCR

The Role of the Glutamate-Glutamine Cycle in Synaptic Transmission During Ischemia and Recovery.

Code ↔ Paper

6 matches between paragraphs of the paper and lines of its authors' code, computed by the harvester (lexical-v1). Click a colored paragraph or line to see its counterpart.

The 6 matches
  1. [1] § Methods › The Model Equations ↔ Dictionary.m, lines 31–45 · score 0.77 · membrane water permeability, universal gas constant, Faraday, moles
  2. [2] § Methods ↔ Dictionary.m, lines 71–125 · score 0.74 · oxygen diffusion, glutamate recycling, oxygen bath, cotransporter, NKCC1, NCX
  3. [3] § Methods ↔ Dictionary.m, lines 71–125 · score 0.67 · oxygen diffusion, glutamate transporter, oxygen bath, cotransporter, NKCC1, NCX
  4. [4] § Results › Glutamate Dynamics During Ischemia › Moderate Ischemia ↔ Plots.m, lines 35–73 · score 0.57 · moderate ischemia, ion concentrations, membrane potentials, potassium, calcium, oxygen
  5. [5] § Results › Glutamate Dynamics During Ischemia › Severe Ischemia ↔ Plots.m, lines 35–73 · score 0.56 · moderate ischemia, ion concentrations, glutamine concentrations, GG, calcium, oxygen
  6. [6] § Methods ↔ Model.m, lines 142–227 · score 0.54 · glutamate cycle, square, duration, KCC, SAT, NCX

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

MATLAB · 197 lines · 8.3 KB · no license · 3 matches

  1. % Creates a dictionary for all constants,parameters
  2. % Written by Hannah van Susteren
  3. function [M] = Dictionary(varargin)
  4. if length(varargin) ==1
  5. varargin = varargin{1};
  6. end
  7. %Put all the default values in a dictionary M.
  8. M = dictionary;
  9. %%
  10. M('tstart') = 0;
  11. M('tend') = 0;
  12. M('time1_OD') = 0;
  13. M('time2_OD') = 0;
  14. M('time1_SB') = 0;
  15. M('time2_SB') = 0;
  16. M('time1_GS') = 0;
  17. M('time2_GS') = 0;
  18. M('time1_EAAT') = 0;
  19. M('time2_EAAT') = 0;
  20. M('alpha_block') = 0;
  21. M('plot') = 0;
  22. M('excite') = 0; % enter a vector of 1x5 of [TSTART TFINAL CURRENT WAVELENGTH PULSEDUTY] to excite the neuron with (CURRENT)pA current starting at TSTART and ending at TEND with wavelength WAVELENGTH and duty DUTY. Values are saved as "excite 1" upto "excite 5"
  23. %% Universal constants
  24. M('F') = 96485.333; % C/mol %Faraday's constant
  25. M('R') = 8314.4598; % (C mV)/(mol K)%Universal gas constant
  26. M('T') = 310; % K %Absolute temperature
  27. M('Zna') = 1; % valence Na
  28. M('Zk') = 1; % valence K
  29. M('ZB') = 1; % valence
  30. M('ZCa') = 2; % valence Ca
  31. M('ZGlu') = -1; % valence Glu
  32. M('Zcl') = -1; % valence Cl
  33. M('ZA') = -1; % valence
  34. M('C') = 20; % pF %membrane constant
  35. M('HeOHa') = 2/3; % Proton ratio (ex:in) in neurons (fixed)
  36. M('Lh2o_n') = 2e-14; % m^3/(ms bar) %Effective membrane water permeability
  37. M('Lh2o_a') = 2e-14; % m^3/(ms bar) %Effective membrane water permeability
  38. %% Permeabilities
  39. %Neuronal channels
  40. M('Pt_Na_n') = 800e-6; %mm^3/ms %Maximal transient Na permeability
  41. M('Pl_Na_n') = 2e-6; %mm^3/ms %Leak Na permeability
  42. M('Pd_K_n') = 400e-6; %mm^3/ms %Maximal delayed rectifier K permeability
  43. M('Pl_K_n') = 2e-5; %mm^3/ms %Leak k permeability
  44. M('Pg_Cl_n') = 19.5e-6; %mm^3/ms %Maximal voltage-gated Cl permeability
  45. M('Pl_Cl_n') = 2.5e-6; %mm^3/ms %Leak Cl permeability
  46. M('Pg_Ca_n') = 1.6e-5; %mm^3/ms %gated Ca permeability
  47. M('Pl_Ca_n') = 4e-10; %mm^3/ms %Leak Ca permeability
  48. M('Pl_Glu_n') = 1e-7; %mm^3/ms %Leak Glu permeability
  49. %Astrocytic channels
  50. M('Pl_Ca_a') = 4e-11; %mm^3/ms %Leak Ca permeability
  51. M('Pl_Glu_a') = 1e-10; %mm^3/ms %Leak Glu permeability
  52. M('Pl_Na_a') = 10e-7;
  53. M('Pl_K_a') = 12e-5;
  54. M('Pl_Cl_a') = 2e-6;
  55. M('Pt_Na_a') = 0.01;
  56. M('g_KDR')=3566.4;
  57. M('P_KDR')=50;
  58. %% Transporters
  59. %KCC transporter
  60. M('Ukcc') = 1.3e-6; % f mol/(ms mV) %KCL cotransporter strength
  61. % %NKA pump from Manu
  62. M('NKA_na') = 13.0; % Half saturation concentration for intracellular Na+
  63. M('NKA_k') = 0.2; % Half saturation concentration for extracellular K+
  64. M('P_NKA') = 110; % Baseline NKA pump strength
  65. M('P_NKA_a') = 90;
  66. %Glutamate recycling
  67. M('k1max') = 1.4e-4; % Glu recycling: Max forward reaction rate
  68. M('KM') = 0.0023; % Glu recycling: Ca half-saturation concentration
  69. M('KDV') = 0.1; % Glu recycling: Half saturation for forward reaction rate
  70. M('k20') = 0.021*1e-3; % Glu recycling: Uncatalysed forward reaction rate
  71. M('k2cat') = 20*1e-3; % Glu recycling: Catalysed forward reaction rate
  72. M('kmin20') = 0.017*1e-3; % Glu recycling: Uncatalysed backward reaction rate
  73. M('kmin1') = 0.05*1e-3; % Glu recycling: Backward reaction rate
  74. M('k3') = 4.4; % Glu recycling: Forward reaction rate
  75. M('kmin3') = 0.001; % Glu recycling: Backward reaction rate
  76. M('k4') = 1.45; % Glu recycling: Fusion rate
  77. M('trec') = 50; % Glu recycling: Vesicle fusion factor
  78. %EAAT neuron
  79. M('alpha_EAAT_n') = 0.0032; % Glutamate transport fitting parameter
  80. M('beta_EAAT_n') = 0.0288; % Glutamate transport fitting parameter
  81. M('rg') = 5e-4;
  82. M('sg') = 9e-2;
  83. %EAAT astrocyte
  84. M('alpha_EAAT_a') = 0.0032; % Glutamate transport fitting parameter
  85. M('beta_EAAT_a') = 0.0288; % Glutamate transport fitting parameter
  86. M('P_EAAT_a') = 120; % EAAT conductance
  87. %NCX
  88. M('P_NCX_n') = 35; % NCX: conductance (1/15 of NKA strength)
  89. M('alphaNaNCX') = 87.5; % NCX: Na half saturation concentration
  90. M('alphaCaNCX') = 1.38; % NCX: Ca half saturation concentration
  91. M('eNCX') = 0.35; % NCX: position of energy barrier
  92. M('ksatNCX') = 0.1; % NCX: saturation factor
  93. M('INCXi0') = -1.4959e-05; % Copied from Manu's computation
  94. %oxygen dynamics parameters
  95. M('alpha_O2') = 5.3/32; % Conversion factor, no unit
  96. M('eps_O2') = 25e-6; % Oxygen diffusion rate, unit: 1/ms
  97. M('O2_e0') = 30;
  98. M('O2_bath') = 40;
  99. M('O2_bath_OD') = 40; % Lower oxygen bath during oxygen deprivation
  100. %Transporters
  101. M('P_Kir') = 0.1;
  102. M('KCe_thres') = 13;
  103. M('P_NKCC1_a') = 100e-7;
  104. M('P_NCX_a') = 1; % NCX conductance
  105. %% Glutamate-Glutamine cycle
  106. %Glutaminase
  107. M('phi_GM') = 1e-7;
  108. M('Km_GM_Gln') = 0.6; % Reference: Kvamme, JNR, 2001
  109. %Glutamine synthetase
  110. M('phi_GS') = 0.3e-6;
  111. M('Km_GS_Glu') = 2.5;
  112. M('Km_GS_NH4') = 0.2;
  113. M('Km_GS_ATP') = 2.3;
  114. M('NH4_a') = 0.15;
  115. %SN transporter
  116. M('Km_SN_Gln') = 1.57;
  117. M('Km_SN_Na') = 31;
  118. M('Km_SN_H') = 100e-6;
  119. M('phi_SN') = 6e-5;
  120. M('P_SN') = 1.755e-4; % fmol/100 mu m^2 ms
  121. M('H_a') = 6.3096e-5;
  122. M('H_e') = 6.3069e-5;
  123. %SAT glutamine transporter
  124. M('Km_SAT_Gln') = 0.3;
  125. M('Km_SAT_Na') = 10;
  126. M('phi_SAT') = 0.00017;
  127. %% Unitial volumes, areas
  128. M('VolPreSyn_n') = 1e-3; % 1000 mu m^3 Presynaptic terminal volume (fixed)
  129. M('VolCleft') = 1e-3; % 1000 mu m^3 Cleft volume (fixed)
  130. M('VolPreSyn_a') = 1e-3;
  131. M('RadiusPreSyn_n') = (M('VolPreSyn_n')*3/4/pi)^(1/3);
  132. M('RadiusPreSyn_a') = (M('VolPreSyn_a')*3/4/pi)^(1/3);
  133. M('AreaPreSyn_n') = 4*pi*M('RadiusPreSyn_n')^2;
  134. M('AreaPreSyn_a') = 4*pi*M('RadiusPreSyn_a')^2;
  135. %% Total molar amounts
  136. %Total sodium, potassium, chloride and volume
  137. M('Tot_N_Na') = 180;
  138. M('Tot_N_K') = 478;
  139. M('Tot_N_Cl') = 1.841874545811657e+02;
  140. M('Tot_Vol') = 4.625000000000000;
  141. %Total calcium,glutamate in the system
  142. M('N_Ca_tot') = 0.001000183047559;
  143. M('N_Gln_tot') = 0.007800000000000;
  144. M('N_Glu_tot') = 0.007010408693616;
  145. M('N_B_a') = 82.485416968700662;
  146. M('N_O2_e0') = 27.750000000000000;
  147. M('N_A_e') = 45;
  148. M('N_A_n') = 2.966822780028186e+02;
  149. M('N_A_a') = 210;
  150. M('P_EAAT_n')=M('P_EAAT_a')/9;
  151. %%
  152. if length(varargin)>1
  153. for i=1:2:length(varargin) % Assign input values to parameters
  154. M(varargin{i}); % If you get an error, one of the labels is misspelled.
  155. if length(varargin{i+1})>1
  156. M(varargin{i})= 1; %the specific key is turned on
  157. for j=1:length(varargin{i+1})
  158. key = join([varargin{i}, j],'_');
  159. M(key)= varargin{i+1}(j); %The different parts of value are saved under 'key 1', 'key 2', ..., 'key n'
  160. end
  161. else
  162. M(varargin{i})= varargin{i+1}; % Assign new values (varargin{i+1}) to labels varargin{i}
  163. end
  164. end
  165. end

Dictionary.m at commit 37d2787, no license · at the source

Overview

  1. Department of Applied Mathematics, University of Twente, Enschede, the Netherlands
  2. Institute of Neurobiology, Heinrich Heine University, Düsseldorf, Germany
  3. Clinical Neurophysiology Group, Department of Science and Technology, University of Twente, Enschede, the Netherlands
  4. Medisch Spectrum Twente, Enschede, the Netherlands
Institutions: University of Twente (Netherlands); Heinrich Heine University Düsseldorf (Germany)
Journal: The European journal of neuroscience, volume 64, issue 1, article e70604
Dates: received 10 November 2025; accepted 12 June 2026; published online 10 July 2026; in print July 2026
Type: Research article · Language: English
License: CC BY
Identifiers: DOI 10.1111/ejn.70604 · PMID 42427222 · PMCID PMC13351822 · OpenAlex W4416186287
Open access: hybrid, a free copy (OpenAlex)
Status: code verified
Categories: computational modeling (no new data) (modality), human (organism), stroke (population), cellular / molecular (subfield)
Keywords: cerebral ischemia, computational model, glutamate, glutamine, neurotransmitter, synaptic transmission
MeSH: Brain Ischemia*, Glutamic Acid*, Glutamine*, Models, Neurological*, Synaptic Transmission*, Animals, Astrocytes, Humans, Neurons (* major topic)
Topic: Neuroscience and Neuropharmacology Research (Cellular and Molecular Neuroscience, Neuroscience), according to OpenAlex
Funding: Deutsche Forschungsgemeinschaft (DFG) (Ro2327/14-2, Ro2327/13-2)
Citations: cited by 1 paper (Europe PMC); 67 references in the paper

Abstract

Cerebral ischemia impairs neuronal and glial function, ranging from transient synaptic failure to irreversible damage. The effects of ischemia on excitatory synaptic transmission remain incompletely understood. Here, we present a detailed biophysical model, including the first full implementation of the glutamate–glutamine cycle (GG‐cycle), which is essential for proper functioning of glutamatergic synapses. We model a presynaptic neuron and an astrocyte in a finite extracellular space (ECS), surrounded by an oxygen bath as a proxy for energy supply. The model includes ionic currents with corresponding channels and transporters such as the sodium‐potassium ATPase. To model synaptic transmission, we combine calcium‐dependent glutamate release, its uptake by the sodium‐dependent excitatory amino acid transporters (EAATs), and the GG‐cycle, including glutamine synthesis. We simulate ischemia by blocking energy supply completely. This drives the neuron into depolarization block, with pathological ion concentrations and extracellular glutamate accumulation despite disrupted synaptic release. Synaptic transmission failure is not primarily caused by excessive glutamate release or by failure of glutamine synthetase, but mainly results from EAAT dysfunction, driven by the collapse of the sodium gradient. Restoring synaptic transmission is not possible by solely targeting glutamate dynamics but is possible by restoring ion gradients by inhibition of the voltage‐gated Na+‐channel. Our study highlights the critical role of ion homeostasis, in particular the sodium gradient, in failure and recovery of synaptic function and the EAAT during metabolic stress.

Reproduced under the paper's license (CC BY), from the paper cited above.

Repository

Its files are read in the Code ↔ Paper reader above, with 6 matches between paragraphs and lines of code.

HannahvanSusteren3/GGsyntrans

License: none: the authors keep all their rights
State: the link answers, verified on 27 September 2026
Evidence: files inventoried
Commit: 37d2787612916a1ee16bcabff10586d08bed1aec, 6 November 2025
Languages: MATLAB (6)
Size: 13 files, 6 scripts
Software Heritage: not archived
Found in: “Data Availability Statement”
Holds: README
Not found: license file, CITATION.cff, environment file, tests, continuous integration, documentation
Availability: 1 check, the latest on 27 September 2026: the link answers
  • 27 September 2026: the link answers
7 files

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:

  • 1 repository of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
  • 6 scripts, each with its path and the digest of its content;
  • 6 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 Statement

The code for all the simulations performed is available at https://github.com/HannahvanSusteren3/GGsyntrans.

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 2, 28 September 2026

  • Publisher: n/a → Wiley

Version 1, 27 September 2026: the first record

Recorded: type, language, journal, volume, issue, pages, dates, 4 authors, 6 keywords, 9 MeSH terms, 1 funder, 64 references.

Cite

This paper

van Susteren, H., Rose, C. R., van Putten, M. J. A. M., & Meijer, H. G. E. (2026). The Role of the Glutamate-Glutamine Cycle in Synaptic Transmission During Ischemia and Recovery. The European journal of neuroscience, 64(1), e70604. https://doi.org/10.1111/ejn.70604

BibTeX

@article{vansusteren2026role,
author = {van Susteren, Hannah and Rose, Christine R and van Putten, Michel J A M and Meijer, Hil G E},
title = {{The Role of the Glutamate-Glutamine Cycle in Synaptic Transmission During Ischemia and Recovery}},
journal = {The European journal of neuroscience},
year = {2026},
month = jul,
volume = {64},
number = {1},
pages = {e70604},
publisher = {Wiley},
issn = {0953-816X},
doi = {10.1111/ejn.70604},
url = {https://doi.org/10.1111/ejn.70604},
pmid = {42427222},
pmcid = {PMC13351822}
}

RIS

TY - JOUR
AU - van Susteren, Hannah
AU - Rose, Christine R
AU - van Putten, Michel J A M
AU - Meijer, Hil G E
TI - The Role of the Glutamate-Glutamine Cycle in Synaptic Transmission During Ischemia and Recovery
T2 - The European journal of neuroscience
J2 - Eur J Neurosci
PY - 2026
DA - 2026/07/01
VL - 64
IS - 1
SP - e70604
SN - 0953-816X
PB - Wiley
DO - 10.1111/ejn.70604
UR - https://doi.org/10.1111/ejn.70604
LA - en
ER -

CSL-JSON

{
"id": "10.1111/ejn.70604",
"type": "article-journal",
"title": "The Role of the Glutamate-Glutamine Cycle in Synaptic Transmission During Ischemia and Recovery",
"container-title": "The European journal of neuroscience",
"author": [
{
"family": "van Susteren",
"given": "Hannah"
},
{
"family": "Rose",
"given": "Christine R"
},
{
"family": "van Putten",
"given": "Michel J A M"
},
{
"family": "Meijer",
"given": "Hil G E"
}
],
"container-title-short": "Eur J Neurosci",
"volume": "64",
"issue": "1",
"page": "e70604",
"DOI": "10.1111/ejn.70604",
"PMID": "42427222",
"PMCID": "PMC13351822",
"ISSN": "0953-816X",
"publisher": "Wiley",
"URL": "https://doi.org/10.1111/ejn.70604",
"language": "en",
"issued": {
"date-parts": [
[
2026,
7,
1
]
]
}
}

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.1093/braincomms/fcag144 [code]
Corticothalamic dynamics during postictal recovery of self-orientation after electroconvulsive therapy.
Journal: Brain communications
In common: 1 reference, author Michel J A M van Putten
[2] doi:10.1038/s41467-026-73435-z [code]
Cellular and subcellular heterogeneity of astrocytic Na⁺ homeostasis tuning astrocytes into functionally distinct subgroups in the mouse brain.
Journal: Nature communications
In common: cellular / molecular, 1 reference, author Christine R Rose
[3] doi:10.1371/journal.pcbi.1014701 [code]
Computer models predict differential dendritic vulnerability with ischemia and spreading depression.
Journal: PLoS computational biology
In common: computational modeling (no new data), 3 references
[4] doi:10.1186/s40478-026-02399-4
Overactive EAAT1 Cl<sup>-</sup> channels impair GABAergic tonic inhibition in SLC1A3-associated episodic ataxia.
Journal: Acta neuropathologica communications
In common: cellular / molecular, author Christine R Rose
[5] doi:10.1038/s41467-026-71391-2 [code]
Accelerating Leigh syndrome drug discovery through deep learning screening in brain organoids.
Journal: Nature communications
In common: author Christine R Rose
[6] doi:10.1038/s42003-026-10044-y [code]
Astrocyte-mediated higher-order control of synaptic plasticity.
Journal: Communications biology
In common: computational modeling (no new data), cellular / molecular, 2 references
[7] doi:10.1016/j.jbc.2026.113091
Structure and enzymology of glutaminase S482C and H461L variants associated with excess brain glutamate and neurological disease.
Journal: The Journal of biological chemistry
In common: cellular / molecular, 2 references
[8] doi:10.1016/j.isci.2026.117487
Endothelial SLC38A2 controls brain D-serine homeostasis and synapse development.
Journal: iScience
In common: cellular / molecular, 2 references
[9] doi:10.1007/s10827-026-00935-8 [code]
A biophysically grounded model of glutamatergic synaptic transmission integrating glutamate transport, receptor kinetics, and electrotonic effects.
Journal: Journal of computational neuroscience
In common: computational modeling (no new data), cellular / molecular, 1 reference
[10] doi:10.1038/s41467-026-76812-w [code]
Assessing molecular, cellular and transcriptomic bases of laminar perfusion and cytoarchitecture coupling in the human cortex.
Journal: Nature communications
In common: stroke, cellular / molecular, 1 reference

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.

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.