Janjic P, Solev D, Kocarev L. Non-trivial dynamics in a model of glial membrane voltage driven by open potassium pores.
Biophys J 2023;
122:1470-1490. [PMID:
36919241 PMCID:
PMC10147837 DOI:
10.1016/j.bpj.2023.03.013]
[Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2022] [Revised: 02/01/2023] [Accepted: 03/07/2023] [Indexed: 03/16/2023] Open
Abstract
Despite the molecular evidence that a nearly linear steady-state current-voltage relationship in mammalian astrocytes reflects a total current resulting from more than one differentially regulated K+ conductance, detailed ordinary differential equation (ODE) models of membrane voltage Vm are still lacking. Various experimental results reporting altered rectification of the major Kir currents in glia, dominated by Kir4.1, have motivated us to develop a detailed model of Vm dynamics incorporating the weaker potassium K2P-TREK1 current in addition to Kir4.1, and study the stability of the resting state Vr. The main question is whether, with the loss of monotonicity in glial I-V curve resulting from altered Kir rectification, the nominal resting state Vr remains stable, and the cell retains the trivial, potassium electrode behavior with Vm after EK. The minimal two-dimensional model of Vm near Vr showed that an N-shape deformed Kir I-V curve induces multistability of Vm in a model that incorporates K2P activation kinetics, and nonspecific K+ leak currents. More specifically, an asymmetrical, nonlinear decrease of outward Kir4.1 conductance, turning the channels into inward rectifiers, introduces instability of Vr. That happens through a robust bifurcation giving birth to a second, more depolarized stable resting state Vdr > -10 mV. Realistic recordings from electrographic seizures were used to perturb the model. Simulations of the model perturbed by constant current through gap junctions and seizure-like discharges as local field potentials led to depolarization and switching of Vm between the two stable states, in a downstate-upstate manner. In the event of prolonged depolarizations near Vdr, such catastrophic instability would affect all aspects of the glial function, from metabolic support to membrane transport, and practically all neuromodulatory roles assigned to glia.
Collapse