A Physical Framework to Study the Effect of Magnetic Fields on the Spike-Time Coding

Abstract

A temporal neural code reliant on the pattern of spike times rather than spike rates offers a feasible mechanism for encoding information from weak periodic external stimuli, such as static or extremely low-frequency electromagnetic fields. Our model focuses on the influence of magnetic fields on neurotransmitter dynamics near the neuron membrane. Neurotransmitter binding to specific receptor sites on membrane proteins can regulate biochemical reactions. The duration a neurotransmitter spends in the bonded state serves as a metric for the magnetic field’s capacity as a chemical regulator. By initiating a physical analysis of ligand-receptor binding, utilizing the alpha function for synaptic conductance, and employing a modified version of Bell’s law, we quantified the impact of magnetic fields on the bond half-life time and, consequently, on postsynaptic spike timing.

Keywords

Magnetic fields bond half-life time bell model alpha function synaptic conductance spike-time coding

Introduction

The existence of biological effects resulting from non-thermal exposure to static or extremely low frequency (ELF) magnetic fields is undeniable; however, it remains a topic of ongoing debate. Investigating the impact of electromagnetic (EM) fields on neuronal activity poses significant challenges. The biological effects triggered by electric or magnetic fields can be examined independently or in conjunction by applying an EM field. Numerous variables come into play, including the broad spectrum of frequencies, varying field strengths, signal waveform, signal duration, field application duration and the shape of the signal. Additionally, the development of experimental models presents considerable diversity, encompassing in vitro or in vivo studies involving animals or human volunteers. This diversity significantly complicates efforts to consolidate results into one or a few unified models.

Experimental data have progressively converged to identify the cell membrane as the primary target structure for the interaction of ELF electromagnetic fields with living organisms.^1
-4 The cell membrane’s capacity to amplify biological signals has been elucidated through various physical models incorporating non-linear mechanisms.^5
-8

Moggia et al⁹ present experimental studies aimed at establishing that among the various explanatory mechanisms involved in the interaction between biological membranes and external EM fields, ligand-receptor binding stands out as the pivotal process. While the Langevin-Lorentz model remains prevalent in the study of cell receptors,^10,11 there has been ongoing development and refinement of this model over the last two decades.

In recent years, numerous studies have utilized magnetic fields to modulate cell activity by targeting cell-surface receptors. An applied magnetic field has been shown to modify Notch receptors, initiating receptor signaling.¹² For E-cadherin, low magnetic forces have been observed to cluster the molecule, thus initiating F-actin assembly.¹³ A similar mechanism has been employed to mechanically deflect and activate inner ear hair cells by targeting membrane glycoproteins.¹⁴ Magnetic field displacement of the hair bundle has been found to increase intracellular calcium, indicating the activation of the hair cells.¹⁵ Magnetic activation of endogenous receptors has also been utilized to enhance stem cell differentiation.¹⁶ Other studies have employed magnetic fields to investigate the role of mechanical force on neurons. Tay et al¹⁷ cultured cortical neurons and discovered that, in the presence of a magnetic field, a mechanical force in the piconewton range was exerted on the neurons, increasing intracellular calcium by targeting membrane calcium channels. Magnetic fields have also been employed to regulate apoptosis. Cancer cells, expressing high levels of the DR4 death receptor, have shown clustering upon application of magnetic fields, subsequently activating caspase-8 and caspase-3, ultimately leading to cell death.¹⁷

The synapse, also known as the synaptic cleft, refers to the space containing the extracellular fluid situated between the pre- and postsynaptic neurons. Its diminutive size, approximately 20 nm wide, facilitates effective neurotransmitter-receptor binding. Fast responses at these synapses are mediated by ligand-gated ion channels (LGICs).¹⁸

The duration, amplitude, and frequency of synaptic responses are determined by the kinetics of channel opening, closing, and desensitization. Binding to LGICs entails a specific and reversible reaction, as it typically involves weak interactions (such as ionic bonds, Van der Waals forces, and hydrogen bonds) between the receptor and ligand.^2,19

Naundorf et al discovered that key features of spike initiation dynamics -such as the spike-time pattern- cannot be explained by the original Hodgkin-Huxley (HH) model. They analyzed the action potential dynamics in cortical neurons by in vivo, in vitro, and computational models²⁰; and found that the spike onset was 10 times greater than the one predicted by the original HH model. Unlike the original HH model, in which the channel opening and closing are independent processes, they proposed a model based on the cooperative action of sodium channels.

While the Naundorf et al model elucidates spike onset characteristics not accounted for in the original HH model, more comprehensive explanations have emerged in recent years. In 2009, Colwell et al introduced a HH-type model incorporating a stochastic term attributed to background synaptic activity. They utilized the path integral formalism to derive an analytical formula that correlates onset speed with onset potential.²¹

In the Colwell HH model, the dynamics of neuron membrane potential $(V$ ) are described as follows,

C_{m} \frac{d V}{d t} = I_{s y n} - I_{N a} - I_{K} - I_{L}

(1)

I_{s y n} = g_{e} (t) (V - E_{e}) + g_{I} (t) (V - E_{I})

(2)

where the input current is implicitly included in equation (1) as an initial condition of the differential equation.

Background synaptic activity is modeled by assuming that synaptic conductance is stochastic, comprising excitatory conductance ( $g_{e}$ ) with a reversal potential $E_{e}$ and inhibitory conductance ( $g_{I}$ ) with a reversal potential $E_{I}$ . The conductances $g_{e} (t)$ and $g_{I} (t)$ are represented by white Gaussian noises, with noise diffusion coefficients $D_{e}$ and $D_{I}$ respectively. Colwell et al examined the correlation between voltage threshold variation and onset speed, enabling them to elucidate the findings of Naundorf et al.

The voltage threshold is defined as the membrane potential value at which $\frac{d V}{d t} = 10 \frac{m V}{m s}$ . As a result of stochastic synaptic background activity, a voltage distribution arises where the voltage threshold is attained. The action potential is triggered when $V$ reaches $V^{*}$ , where $V^{*}$ represents an unstable equilibrium of the aforementioned equation in the absence of noise. Below $V^{*}$ , the membrane potential returns to its resting state, while above $V^{*}$ , a spike can be initiated. Anastassiou et al have studied the problem of spike-timing control through external electromagnetic fields, we provide a more detailed explanation in the Supplemental Information section. Hence, in his model, the enhancement in coding properties stems from the stimulus input not arising from stochastic fluctuations but from a deterministic external field. Nonetheless, the model proposed by Anastassiou et al is purely phenomenological. Here, our aim is to establish a mechanistic model aligned with the one we will detail in subsequent sections, potentially enabling us to broaden the findings for studying temporal coding within local neural networks.

In this paper, we propose a modified HH model based on Colwell’s advancements. In our model, the noise term arises not only from a stochastic component but also from a deterministic one, which encompasses the synapse-induced magnetic field. We adjust the Bell equation, expressing the bond half-life time, to account for the dynamics of neurotransmitters. The paper is structured as follows: firstly, we introduce a physical framework to examine the influence of magnetic fields on biological dipoles under physiological conditions. Secondly, we establish a connection between ligand-receptor binding and the synapse-induced magnetic field, deriving the explicit solution for synaptic conductance. Finally, we analyze the variation in postsynaptic spike timing as a function of magnetic fields and the stability analysis of the solutions.

Methods

Ionic conductance properties and parameters

The ionic conductances, excluding the leakage current $g_{L}$ , which was assumed constant, are described by the following equations (3) to (4) for potassium conductance and (5) to (7) for sodium conductance,

g_{K} = {\bar{g}}_{K} n^{4},

(3)

\frac{d n}{d t} = α_{n} (1 - n) - β_{n} n,

(4)

g_{N a} = {\bar{g}}_{N a} m^{3} h,

(5)

\frac{d m}{d t} = α_{m} (1 - m) - β_{m} m,

(6)

\frac{d h}{d t} = α_{h} (1 - h) - β_{h} h,

(7)

where $α_{n}, β_{n}, α_{m}, β_{m}, α_{h}$ and $β_{h}$ are free parameters and ${\bar{g}}_{K}$ and ${\bar{g}}_{N a}$ represent the maximum ionic conductances. The dimensionless variables $m, n$ and $h$ describe the open probability of ion channels (and thus the expected fraction of the full ionic conductance across the membrane at a given time and voltage), while the $α$ and $β$ terms denote rate constants.

After selecting rate parameters $α$ and $β$ for various clamped voltages, Hodgkin and Huxley empirically proposed functions capable of explaining the voltage dependency of these parameters, resulting in the following equations (8) to (9) for potassium conductance and (10) to (13) for sodium conductance,

{\hat{α}}_{n} (V) = \frac{0.01 (V + 10)}{e x p (\frac{V + 10}{10}) - 1},

(8)

{\hat{β}}_{n} (V) = 0.125 \exp (\frac{V}{80}),

(9)

{\hat{α}}_{m} (V) = \frac{0.1 (V + 25)}{e x p (\frac{V + 25}{10}) - 1},

(10)

{\hat{β}}_{m} (V) = 4 \exp (\frac{V}{18}),

(11)

{\hat{α}}_{h} (V) = 0.07 \exp (\frac{V}{20}),

(12)

{\hat{β}}_{h} (V) = \frac{1}{e x p (\frac{V + 30}{10}) - 1},

(13)

All HH’s conductance data were sourced from Marom and Abbott.²² Given the absence of computers in 1952 for conducting numerical solutions, HH’s numerical algorithm was designed for manual calculations using the tools available at the time. The numerical method employed in this study is known as Euler’s method.²³ All simulations using Euler’s method were performed using a contemporary implementation in MATLAB,²⁴ equipped with metadata tags to facilitate adjustment of the $α$ and $β$ parameters (equations (4), (6), and (7)) through a fitting algorithm. Similarly, the constants in equations (8) to (13) were replaced with free, externally adjustable variables as outlined in equations (14) to (19), thereby introducing 5 free parameters for the voltage-dependent potassium conductance $g_{K}$ and 9 free parameters for the voltage-dependent sodium conductance $g_{N a}$ .

When fitting the limited six-parameter form of the model (equations (4), (6), and (7)), the rate parameters $α_{n}, β_{n}, α_{m}, β_{m}, α_{h}$ , and $β_{h}$ are defined by the protocol to be constants, which are set by the fitting algorithm. When fitting the full 14-parameter voltage-dependent model, these six rate parameters are in turn parametrized according to the following equations (14) to (15) for potassium conductance type and (16) to (19) for sodium conductance, with the 14 free parameters set by the fitting algorithm,

α_{n} (V) = \frac{k_{α_{n} 1} (V + k_{α_{n} 2})}{e x p (\frac{V + k_{α_{n} 2}}{k_{α_{n} 3}}) - 1},

(14)

β_{n} (V) = k_{β_{n} 1} \exp (\frac{V}{k_{β_{n} 2}}),

(15)

α_{m} (V) = \frac{k_{α_{m} 1} (V + k_{α_{m} 2})}{e x p (\frac{V + k_{α_{m} 2}}{k_{α_{m} 3}}) - 1},

(16)

β_{m} (V) = k_{β_{m} 1} \exp (\frac{V}{k_{β_{m} 2}}),

(17)

α_{h} (V) = k_{α_{h} 1} \exp (\frac{V}{k_{α_{h} 2}}),

(18)

β_{h} (V) = \frac{1}{e x p (\frac{V + k_{β_{h} 1}}{k_{β_{h} 2}}) + 1},

(19)

$k_{α_{n} 1}, k_{α_{n} 2}, k_{α_{n} 3}, k_{β_{n} 1}, k_{β_{n} 2}, k_{α_{m} 1}, k_{α_{m} 2}, k_{α_{m} 3}, k_{β_{m} 1}, k_{β_{m} 2}, k_{α_{h} 1}, k_{α_{h} 2}, k_{β_{h} 1}$ , and $k_{β_{h} 2}$ are free parameters.

Membrane potential properties and parameters

All membrane potential data were sourced from Marom and Abbott.²² The voltage differential between the cell’s interior and the synapse influences both the excitatory synaptic response (EPSP) and spike generation. Under resting conditions, a neuron typically maintains a membrane potential of approximately $- 65 mV$ (ranging from $- 60$ to $- 70 mV$ ), with the interior of the cell being negatively charged relative to the external environment (ie, $V_{t = 0} = V_{R} = V_{i n} - V_{o u t} = - 65 mV$ ).

The membrane potential of the neuron is contingent upon the disparate concentrations of $K^{+}$ and $N a^{+}$ ions on either side of the membrane, as well as the varying membrane permeability to these ions, with $K^{+}$ exhibiting significantly higher permeability.

Potassium ions ( $K^{+}$ ) are predominantly concentrated within the cell’s interior, with a concentration reaching approximately 150 mM ( ${[K^{+}]}_{i}$ ), while outside the cell ( ${[K^{+}]}_{e}$ ), it typically remains around 4 to 5 mM. By employing the Nernst equation, we can calculate the equilibrium potential of the equivalent electrochemical cell,

E_{K} = \frac{k_{B} T}{z F} L n \frac{{[K^{+}]}_{i}}{{[K^{+}]}_{e}}

(20)

Here, $k_{B}$ represents the Boltzmann constant, $T$ denotes the absolute temperature, $F$ stands for the Faraday constant, and $z$ signifies the number of electrons transferred in the electrochemical reduction reaction (which is $+ 1$ ). At the mentioned physiological concentrations, $E_{k}$ is approximately $- 90 mV$ .

Conversely, $N a^{+}$ ions have an internal concentration of $5 - 10 mM$ and an external concentration of $140 mM$ , resulting in an equilibrium potential $E_{N a}$ of approximately $+ 50 mV$ . $K^{+}$ and $N a^{+}$ ions possess opposing electrochemical gradients, with the gradient for $K^{+}$ ions driving outward currents (positive current) from the cell, while that for $N a^{+}$ ions pushes them from the external to the internal cell space (negative current). The significantly higher membrane permeability to $K^{+}$ ions maintains the membrane potential at a depolarization level equivalent to the equilibrium potential for $K^{+}$ ( $- 90 mV$ ). Sodium-potassium channels act as electrical pumps for the neuron membrane.

The neuron expends metabolic energy to uphold the membrane voltage at a less depolarized level of $- 65 mV$ . This function is carried out by biochemical electrogenic pumps that transport $K^{+}$ ions from the external to the internal cell space and $N a^{+}$ ions in the opposite direction, against their electrochemical gradients. These pumps, known as $N a^{+}$ / $K^{+}$ pumps, rely on the ATP biochemical conversion of energy and establish the resting membrane potential (a potential where the net current crossing the membrane is $0 A$ ) at $V (t) = V_{R} = - 65 mV$ . At this level, the voltage-dependent conductance for $N a^{+}$ and $K^{+}$ ions ( $G_{N a}$ and $G_{K}$ , respectively) is inactivated ( $0 S$ ), and the resistance to these ions can be considered infinite (factoring in the pump action), resulting in respective currents of are 0 A ( $I_{N a} = I_{K} = 0$ in equation (1)). $G_{N a}$ and $G_{K}$ undergo rapid increases for membrane potentials depolarized to almost $+ 10 mV$ , which can be deemed an activation threshold value ( $V_{t h}$ = $V_{r} + 10 mV = - 65 + 10 = - 55 mV$ ).

At the threshold level ( $V_{t h} = - 55 mV$ ), there is a rapid membrane depolarization triggered by a negative net current resulting from the swift inward movement of $N a^{+}$ ions. This is attributed to the significantly higher driving force at this potential level for this ion. Specifically, the ionic current can be represented as $I_{i o n}$ = $G_{i o n} (V (t) - E_{i o n}$ , where the driving forces generating the respective $N a^{+}$ and $K^{+}$ currents are calculated as follows: $V (t) - E_{N a} = - 55 mV - (+ 50 m V) = - 105 mV$ ; $V (t) - E_{K} = - 55 mV - (- 90 m V) = + 35 mV$ .

The negative, inward total current not only governs the rapid rise phase of the spike but also contributes to the synaptic current. This current initiates depolarization of the membrane potential, reaching and often surpassing $0 mV$ , leading to overshooting potentials of $+ 20$ , $+ 30 mV$ at the spike peak. As the membrane potential $(V (t))$ enters positive territory, the driving force on $N a^{+}$ ions diminishes while that on $K^{+}$ ions increases. The net current definitively reverses direction at the spike peak ( $+ 20$ , $+ 30 mV$ ), and the membrane begins repolarization as the $K^{+}$ ion battery approaches the ion equilibrium potential ( $- 90 mV$ ). This repolarization is facilitated by the rapid inactivation of $N a^{+}$ conductance ( $G_{N a}$ returning swiftly to $0 S$ ). Subsequently, the energy-dependent work of the $N a^{+}$ / $K^{+}$ pump restores the resting potential value ( $V_{r} = - 65 mV$ ), completing the spike cycle. The membrane dynamics, encompassing both spike and synaptic activity, can be likened to an electrical circuit operating in accordance with the aforementioned parameters and following the dynamics outlined in the HH model.

The sole distinction, from an electrical standpoint, between synaptic response and spike generation lies in the driving force behind the synaptic current. Unlike spike generation, where the driving force stems from the disparity between membrane voltage and the ion-specific equilibrium potential, excitatory receptors (primarily AMPA receptors) constitute nonspecific ionic channels. Consequently, the net current is determined by the synaptic conductance (dependent on neurotransmitter-receptor binding) and a nonspecific synaptic equilibrium potential, typically around $- 10 mV$ . This signifies that membrane potentials exceeding $- 10 mV$ result in inward currents (negative and depolarizing), while those below produce outward currents (repolarizing), primarily due to outward flow of $K^{+}$ ions from the cell. This framework offers a straightforward explanation for the dynamics of neuron membrane voltage.

Physical Derivation of the Model

Electric dipole moment of neurotransmitters under physiological conditions

Finite electric dipoles typically consist of two-point charges of equal magnitude: a negative charge ( $- q$ ) and a positive charge ( $q$ ), with the relative position vector $\vec{d}$ extending from the first to the second. Therefore, we define the electric dipole moment as follows,²⁵

\vec{p} = q \cdot \vec{d}

(21)

The point dipole moment is obtained in the limit as $q \to \infty$ and $d \to 0$ . Thus, the expressions in equation (24) can be derived directly from Tsang,²⁶

\vec{p} = \hat{k} \int_{0}^{π} \int_{0}^{2 π} a^{2} \sin θ d θ d φ ρ_{0} c o s^{2} θ = \hat{k} \frac{4 π a^{3} ρ_{0}}{3}

(22)

here, $a$ represents the radius of the electric dipole, assuming the dipole as a spherical particle, and $ρ_{0}$ denotes the surface charge density.

Sivasankar et al²⁷ determined the value of $ρ_{0}$ for neurotransmitters under physiological conditions and obtained an average value of $12.5 \frac{m C}{m^{2}}$ (range $1 - 20$ ). Assuming a neurotransmitter size of $a$ = $3 \overset{\circ}{A}$ (and $ρ_{0} = 12.5 \frac{m C}{m^{2}}$ ),²⁸ we obtain a dipole moment ( $p$ ) value of $1.979 \cdot 10^{- 30} m \cdot C$ . Antosiewicz and Porschke²⁹ determined the dipole moment ( $p$ ) value for neurotransmitters in aqueous solution at physiological pH and obtained a range of variation of $(1 - 10) \cdot 10^{- 30} m \cdot C$ .³⁰

Dynamics of electric dipoles in magnetic fields

A current element $I d l$ immersed in a magnetic field $B$ will experience a force $d F$ given by,³¹

d F = I \times B d l

(23)

The current $I$ in equation (26) is attributed to electrons moving within the medium and is referred to as conduction current. We can adapt equation (23) to explicitly address charges moving in free space, constituting what is known as a convection current. A charge $q$ moving with velocity $v$ is equivalent to a current element $I d l = q v$ , thus we obtain,

F = q v \times B (24)

(24)

When a charged particle encounters a static field at right angles, it undergoes a force known as the Lorentz force, which alters the particle’s motion into uniform circular motion around the magnetic field lines. According to Newton’s second law, the radial force acting on the particle is $F_{r} = m \cdot a_{r} = q \cdot v_{⊥} \cdot B$ , where $v_{⊥} = v \cdot s i n θ$ represents the velocity component perpendicular to the magnetic field. The radial or centripetal acceleration is given by $a_{r} = \frac{{v_{⊥}}^{2}}{r} = w^{2} \cdot r$ , where $r$ is the radius of the circle and $w$ is the angular velocity of the charged particle orbiting with mass $m$ . The parallel component of velocity, $v_{∥} = v \cdot c o s θ$ , remains unaffected by the presence of $B$ . Therefore, when a particle enters a magnetic field at an angle to the magnetic field lines, it follows a helical path around the magnetic field lines, with $v_{∥}$ unchanged and $v_{⊥}$ causing the circular motion (Figure 1).³²

Figure 1.

(a) Trajectories of charged particles in a magnetic field and (b) a charged particle with a velocity component parallel to a magnetic field and another perpendicular to it, moves in a helical path around the field lines.

The energy principle can be applied to derive force distributions. To develop a simple model for force density distribution, we require the expression for the force on an electric dipole in polarizable media, which can be derived from the Lorentz force law. The force on a stationary electric charge is governed by the Lorentz law with a frequency of zero. Considering two charges $q$ of opposite polarity separated by a vector distance $d$ , the total force on the dipole is the sum of the forces acting on the individual charges.

f = q E

(25)

where $E$ is the electric field. The total force on the dipole is the sum of the forces acting on the individual charges, given by,

f = q [E (r + d) - E (r)]

(26)

Unless the electric field at the location $r + d$ of the positive charge differs from that at the location $r$ of the negative charge, the separate contributions cancel each other out. To derive an expression for the force on the dipole in the limit where the spacing $d$ of the charges is small compared to distances over which the field varies appreciably, equation (26) is expressed in Cartesian coordinates. The field at the positive charge is then expanded about the position of the negative charge. Thus, the $x$ component is,

\begin{matrix} f_{x} = q [E_{x} (x + d_{x}, y + d_{y}, z + d_{z}) - E_{x} (x, y, z)] \\ = q [\begin{array}{l} E_{x} (x, y, z) + d_{x} \frac{\partial E_{x}}{\partial x} + d_{y} \frac{\partial E_{y}}{\partial y} + d_{z} \frac{\partial E_{z}}{\partial z} \\ + \dots - E_{x} (x, y, z) \end{array}] \end{matrix}

(27)

The first and last terms cancel each other out. In a more concise notation, this expression is thus represented as,

f_{x} = p \cdot \nabla E_{x}

(28)

where we have identified the dipole moment $p = q d$ . The other force components are derived similarly, with $y$ and $z$ taking the place of $x$ . These 3 components are then summarized in the vector expression,

f = p \cdot \nabla E

(29)

Time-scale hierarchy in complex systems

In dynamics involving far-from-equilibrium irreversible processes, the relax integrability function tends to converge toward stable asymptotic configurations. The relaxation time is contingent upon the dissipative process parameters operating at each level. Neuronal time scales exhibit dynamism, with a succession of relaxation times discernible, each succeeding level characterized by an increase in magnitude^33,34:

Microsecond scale, associated with the rearrangement of ionic concentrations and reactive phenomena in tissue electrolytes.

Hundreds of microseconds scale, associated with the activation of sodium channels and subsequently with calcium channels involved in membrane excitation.

Millisecond scale, associated with the membrane time constant, pulse and stimulation duration, and electrochemical double layer polarization at the membrane-electrolyte interface.

Hundreds of milliseconds scale, associated with the mechanisms of activation and inactivation of membrane channels related to recovery phenomena.

Minute scale, involved in membrane adaptation phenomena (slow excitability variations dependent on ionic accumulation in regions where transport processes are disturbed).

The differences in the numerical order of characteristic relaxation times among hierarchical scales in dissipative systems often enable simplification of their description by applying 2 connecting scale principles. First, on a particular time scale, the variables of the higher time scale can be treated as constant parameters (as their relaxation times are significantly greater). Second, the variables of the lower time scale can be disregarded (as their relaxation times are much smaller).³⁵

Suppose we have an electric dipole initially at rest within a static vertical magnetic field. The movement of the dipole depends on the EM field configuration, specifically on the electric potential $(V_{e})$ generating the field, rather than the magnetic field itself, which remains constant. Zemped et al investigated the effect of small extracellular fields on spike-time variation. Although they did not propose an explanatory mechanism, they provided an experimental model to refute a potential explanation based on the resonant properties of neurons.³⁶ In their studies, they applied magnetic fields in the millitesla (mT) range, and the voltage measurement was approximately 2 mV, with a range of 1 to 5 mV. We will use this value in the following calculation.

While the magnetic field may not be conservative, the EM field as a whole is conservative. Therefore, applying the principle of conservation of energy to a dipole positioned in a region of potential $V_{e}$ , the potential energy $\frac{p}{d} \cdot V_{e}$ is converted into kinetic energy.

\frac{1}{2} m {v_{z}}^{2} = q V_{e} \Rightarrow v_{z} = \sqrt{\frac{2 q V_{e}}{m}} = \sqrt{\frac{2 \frac{p}{d} V_{e}}{m}}

(30)

Taking $p = 1.979 \cdot 10^{- 30} m \cdot C$ ; $d = 2 \cdot a = 2 \cdot 3 \cdot 10^{- 10} = 6 \cdot 10^{- 10} m e t e r s; m = 6 \cdot 023 \cdot 10^{- 26} k g$ , $V_{e} = 2 m V$ and a synaptic cleft width of $20 n m,$ we obtain a $v_{z}$ value of $14.491 \frac{m e t e r s}{s}$ ³⁷ and a transit time of $1.38 n s (1.38 \cdot 10^{- 9} s)$ .³⁸

Therefore, the electrical potential associated with extremely weak static magnetic fields induces temporal variations in neurotransmitter dynamics at a much lower time scale, from a hierarchical perspective, compared to what is typical in classical studies of chemical synapses. This suggests that the very weak synapse-induced magnetic fields may serve as fine control mechanisms in shaping the overall process of postsynaptic spike-time patterns.

Dynamics of hybrid multiscale systems: Combining deterministic and stochastic elements

Deterministic and stochastic hybrid systems no longer evolve with a deterministic trajectory in the parameter space, but rather with a certain probability distribution defined within that same space. A general formulation for these systems is based on nonlinear differential equations incorporating both a stochastic noise term, which models fluctuations, and a deterministic term.^39,40 The statistical properties of the noise completely define the system through its stochastic differential equation (SDE). This additive noise, resulting from fluctuation timescales and system response time, is typically assumed to have zero correlation time (white noise). Additionally, there may be an additional source of fluctuations arising from environmental influences, known as external noise.

The most common assumption regarding the stochastic term present in the SDE is that it follows a normal distribution, meaning it is Gaussian. In many systems, noise represents the influence of random factors whose cumulative effect, according to the central limit theorem, converges to a Gaussian distribution. White noise intensity is constant in the frequency-domain. Its correlation function is defined as follows,⁴¹

ξ (t) ξ (t^{'}) = 2 D δ (t - t^{'})

(31)

Here, the parameter $D$ represents the intensity of the noise. White noise represents the extreme case of a stochastic force evolving at an exceedingly rapid pace. From a physical standpoint, this indicates that the fluctuation time scale is significantly smaller compared to typical system response times. This heuristic approximation method is used to account for internal fluctuations and serves as a starting point for expanding the analysis to include external fluctuations in oscillating and excitable systems.

In a region of local instabilities in phase space induced by external stochastic influences, such as thermal motions or fluctuations in boundary conditions, a particle can be significantly affected by a deterministic force field. The existence of deterministic chaos is not a prerequisite for this phenomenon. If a deterministic force field is sufficiently strong and persistent, deviations from a stochastic trajectory can be profound, leading to predictable, even deterministic trajectories that are asymptotically stable. The outcome is noisy behavior that would not be observable without the influence of deterministic dynamics but cannot be sustained by deterministic dynamics alone.

The probability distribution $P (x, t)$ for systems with white noise follows a Fokker-Planck equation. Technically, white noise is a random signal where the values of the signal at 2 different times are not statistically correlated. However, white noise serves as an approximation of true noise, as real systems subjected to external fluctuations typically do not have a zero-correlation time. The most effective approach to address this issue is to model the random signal as an Ornstein-Uhlenbeck noise, which is a Gaussian noise characterized by a correlation function that decreases exponentially, expressed as,⁴²

ξ (t) ξ (t^{'}) = \frac{D t}{τ} e^{- \frac{| t - t^{'} |}{τ}}

(32)

The parameters $D$ and $τ$ represent the intensity and the noise correlation time, respectively. White noise corresponds to the limit as $τ$ approaches zero.

Extended kinetics of the ligand-receptor bell model incorporating synaptic magnetic fields

The temporal evolution of the reversible chemical reaction between 2 reactants A and $B$ , forming the product $A B$ ( $A + B \Leftrightarrow A B$ ), can be described by the equation,

\frac{d [A B]}{d t} = k_{o n} [A] \cdot [B] - k_{o f f} [A B]

(33)

Here, $[A]$ , $[B]$ , and $[A B]$ represent the concentrations of the reactants and product, respectively. $k_{o n}$ and $k_{o f f}$ denote the association and dissociation rate constants, which govern the kinetics of the interaction. $k_{o n}$ is influenced by the diffusive properties of the molecules and depends on factors such as distance and orientation. On the other hand, $k_{o f f}$ is linked to the bond half-life time, denoted as $τ = {(k_{o f f})}^{- 1}$ ,⁴³ providing insights into the specificity and spontaneity of the reaction,

k_{o f f} = A \cdot e^{\frac{- Δ G^{0}}{k_{B} T}}

(34)

In this equation, $A$ represents the Arrhenius constant, $k_{B}$ stands for the Boltzmann constant, $T$ denotes the absolute temperature, and $Δ G^{0}$ represents the activation free energy of the reaction under standard conditions.

In our study, we introduce a theoretical model rooted in mechanics force fields theory to simulate the impact of magnetic forces on molecular interactions, particularly chemical bonds. When an external force is exerted, it alters the activation energy of the chemical reaction. The Bell model establishes a relationship between the magnitude of the force and the dissociation constant of the reaction,⁴⁴

k_{o f f} (F) = {(k_{o f f})}_{0} \cdot e^{- (\frac{F \cdot x_{β}}{k_{B} \cdot T})}

(35)

In the equation, $F$ represents the external force, ${(k_{o f f})}_{0}$ denotes the equilibrium dissociation rate constant when the external force on the system is zero, and $x_{β}$ represents the activation energy of the chemical reaction. Rewriting the equation in terms of the bond half-life time,

τ = τ_{0} \cdot e^{\frac{F \cdot x_{β}}{k_{B} \cdot T}}

(36)

where $τ_{0} = {(k_{o f f})}_{0}^{- 1}$ . If we incorporate the typical parameters of LGIC-type binding into equation (39),⁴⁵ we derive $τ_{0} = 10^{- 8}$ seg. In this context, we suggest employing a slightly adjusted version of the Bell equation that incorporates the synapse-induced magnetic field $B$ , expressed as,

τ = τ_{0} \cdot e^{\frac{B \cdot \frac{p}{d}}{k_{B} \cdot T}}

(37)

Here, $p$ and $d$ represent the dipole moment and the diameter of the neurotransmitter, respectively.

Neuron models that integrate alpha function conductance-based synapses

Synaptic noise acts as the primary origin of membrane potential fluctuations in neurons and plays a significant role in shaping their integrative behaviors, including synaptic conductance. Destexche and Rudolph-Lilith⁴⁶ conducted a review of experimental synaptic noise measurements and formulated a model grounded in stochastic noise processes.

A common approach to estimating the alteration in conductance within the postsynaptic neuron is through the alpha function model. In numerous synapses, the progression of synaptic conductance undergoes a finite-duration phase. The alpha function delineates this phase of conductance progression as not infinitely rapid, implying it possesses a finite rise time,

g_{s} (t) = \frac{t}{τ} e^{- \frac{t}{τ}}

(38)

It is worth noting the similarity between equations (32) and (38), which validates the application of stochastic-deterministic noise theory to conductance-based synapses. Consequently, we can express the kinetics of synaptic conductance as follows,

g_{s y n} = {\bar{g}}_{s y n} \cdot g_{s} (t) = {\bar{g}}_{s y n} \cdot \frac{t}{τ} e^{- \frac{t}{τ}}

(39)

Hence, we eventually deduce a membrane potential dynamics equation rooted in the Colwell HH model, integrating the synapse-induced magnetic field,

I_{s y n} = g_{s y n} (t) [V (t) - E_{s y n}]

(40)

Here, $E_{s y n}$ denotes the constant voltage source.

Figure 2 depicts the temporal progression of synaptic conductance ( $g_{s y n} (t)$ ) for different magnitudes of $B$ .

Figure 2.

Synaptic conductance analytic solutions using the alpha-function model at different synapse magnetic field values ( $B$ ).

Results and Discussion

Cell membrane proteins can be arranged as a kind of electric amplifiers, such as is used in electronics. Every periodic current can be analyzed using Fourier series as a sum of a constant component (DC) and alternating components. These are entirely transformed by the transformer, all in the same ratio, while the DC component induces no voltage across the secondary. Consequently, the transformed voltage is, except for the transformation factor, identical to the initial voltage, minus its mean value—(Figure 3a)—(the transformed voltage has a zero-mean value).

Figure 3.

In (a), we can observe the decomposition of a rectified alternating current (primary current) into a direct current (untransformed) and a rectified alternating current with zero average offset. In (b), the phase shift of the intensity of a sinusoidal alternating current due to an external disturbance is presented, in this case, a magnetic field.

The phase shift $(Δ ϕ)$ is contingent upon both the phase ( $ϕ$ ) and the magnitude of the disturbance. For extremely brief perturbations, $Δ ϕ$ remains zero across the entire phase spectrum. This domain corresponds to the oscillation’s refractory zone, wherein no effect is observed regardless of the phase value of the applied disturbance. Beyond this zone, the disturbance consistently induces a phase shift, the magnitude of which hinges on the phase value. At the boundary between the refractory and sensitive zones, there exists a discontinuity in amplitude directly influenced by the intensity of the disturbance. This discontinuity underscores the excitability nature of the system.

From the perspective of alternating current circuit theory, it’s feasible to manipulate the phase of the resulting electric current (the initiation timing) by adjusting the amplitude of the sinusoidal potential at the moment $t = 0$ when the circuit closes. Let’s consider a capacitor (cell membrane), initially in an open circuit, perturbed by a sinusoidal voltage $U = U_{0} s i n w t$ at $t = 0$ (corresponding to the periodic extracellular EM disturbance associated with the induced magnetic field). Initially, the current remains zero until the circuit is closed, and begins to rise as the capacitor charges. However, this represents a transient regime primarily due to the non-zero resistance of the capacitor (the duration of this transient regime is approximately of the order of RC, where R and C denote the resistance and capacitance of the capacitor). Subsequently, the steady state is achieved, where the current leads the voltage. In other words, the EM signal doesn’t manifest instantaneously but undergoes a transient phase (Figure 3b)).

A single excitatory synapse, such as a glutamatergic synapse, which accounts for nearly 90% of synapses in the cortex, generates a current ranging from 5 to 110 pA, with an average of 25 pA, and a maximal conductance of 700 to 1000 pS (ranging from 0 to 1000 pS). The reliability of a synapse to elicit or influence spikes depends not only on its conductance but also on its proximity to the spike initiation point, typically located near the soma where the axon emerges, and often distant from the synapse itself. At most, a single synapse produces 100 pA at the origin, with a local membrane potential variation of less than 3 to 5 mV, insufficient to surpass the threshold for firing. As the potential wave travels from the origin to the soma, where the spike is generated, these values decrease due to dendritic cable properties. Typically, a synaptic event occurring in the dendritic tree induces less than 0.5 mV in membrane potential, which, when added to the resting level, fails to reach the threshold. These observations lead to a significant conclusion: firing a spike necessitates the temporal integration of currents arriving from multiple synapses.

A post-synaptic current intensity of $10^{- 10} A$ generates a magnetic field near the dendrite, at the external cranial base, of approximately $10^{- 4} fT (1 fT = 1 f e m t o t e s l a = 10^{- 15} T e s l a)$ . Consequently, this field is undetectable, as current sensors typically necessitate a minimum field strength of $100 fT$ . To register this field with a sensor, a simultaneous contribution of 1 million post-synaptic currents would be required.

In a recent study,⁴⁷ we explored the induced magnetic fields within living tissues. Our findings suggest that if we acknowledge the presence of off-diagonal terms in the synapse conductivity tensor, an exogenous EM field can indeed trigger a synapse magnetic field, thereby altering neurotransmitter dynamics.

Hence, the following question emerges: Could an exogenous EM field, generating a weak synapse magnetic field ( $10 - 100 fT)$ , potentially modify the temporal patterns of postsynaptic action potentials? Presently, there is limited research exploring the connection between magnetic fields and synapses. Nevertheless, a few studies have investigated the impact of magnetic fields within the mT range on action potential patterns.^48
-50

When applying a continuous input of 1 $μ A$ , while maintaining the remaining chemical-physical variables constant, but altering the magnetic field within the range of $0.001 - 5 mT$ , we observe a nonlinear variation in spike timing with increasing $B$ (Figure 4a)). Changes in spike timing directly correlate with fluctuations in bond half-life and, consequently, neurotransmitter dynamics. As a result, magnetic fields at the synapse induce a time-varying spike pattern (Figure 4b)). Additionally, we have modeled the scenario with a continuous input of 5 $μ A$ , considering magnetic field values ranging from $0.001 - 10 mT$ in Figure 5. As we may observe in Figure 4b the postsynaptic spike-time results are highly non-linear, resulting from a number of non-linear voltage dependent currents acting in context. Instead of instantiating these complexities, we focused on making the analysis of the results as simple as possible and relegated ourselves to establish the qualitative conclusions.

Figure 4.

Numerical simulation of the model with input driving current of 1 $μ$ A and induced magnetic fields of 0.001, 0.1, 1, 2.5, and 5 mT. In (a) membrane voltaje time-variation. In (b) time-varying spike pattern.

Figure 5.

Numerical simulation of the model with input driving current of 5 $μ$ A and induced magnetic fields of 0.001, 1, and 10 mT.

Anastassiou et al^51,52 investigated the impact of spatially inhomogeneous extracellular electromagnetic fields on neurons, specifically focusing on alterations in spike timing. They observed that extracellular magnetic fields, within the same order of magnitude as those utilized in our study, induced changes in somatic membrane potential of less than 0.5 mV under subthreshold conditions. Furthermore, they conducted simultaneous recordings from up to 4 patched neurons located proximally to each other. Their findings suggest that despite their small magnitude, these magnetic fields could significantly influence spike timing, aligning with the theoretical outcomes of our study.

We can analyze the stability of disturbances within the framework of stable oscillator theory. Figure 6 illustrates the phase shift among identical oscillators with matching periods, initially synchronized, upon application of a continuous current input of 100 $μ A$ , albeit with varying $B$ values. The spike trains could be represented like a sine wave, its amplitude is represented on vertical axis and time on horizontal axis. As we have previously discussed (Figure 3b), this process has an initial phase with highly unstable behavior (transient phase, varying in duration depending on the magnitude of the applied magnetic field) with many artifacts due to the large number of nonlinear equations involved. So, in Figure 6 we focused on the difference between the first spikes of the spike trains corresponding to the different applied magnetic fields $(Δ ϕ)$ .

Figure 6.

Numerical simulation of the model with input driving current of 100 $μ$ A and induced magnetic fields of 0.001 and 1 mT.

Stability analysis of the spike-timing model based on bi-dimensional dynamic systems

Our modified HH model is a dynamic system consisting of the state variables $V, n, m y h$ , governed by a 4-dimensional system of ordinary differential equations. This HH model can be reduced to a 2-dimensional model and still produce the same action potentials.

The system state is a 2-dimensional vector $(V, n)$ in the phase plane $R^{2}$ because the current exhibits instantaneous activation kinetics, such that its conductance can be considered maximal, $m_{\infty}$ , over most of the time interval. Activation kinetics for $n$ are much slower and therefore need to be defined by its derivative.

First, we need to establish the equilibrium points to study the nullclines of the state variables. A nullcline is the set of locations in the phase plane where the state variable is at rest. In this system, it would correspond to the conditions $\dot{V} = 0$ and $\dot{n} = 0$ .

The intersections of these nullclines will be the points where none of the state variables are changing, thus indicating membrane equilibrium. In Figure 7, we observe the cubic-shaped nullcline for $V$ in red and the sigmoid-shaped nullcline for $n$ in green. Their intersection is the equilibrium point for this model.

Figure 7.

Representation of the nullclines of the model with $I = 1 p A$ and $B = 0 m T$ . The curves appear to have an equilibrium point approximately at $V = - 62 m V$ .

Unlike 1-dimensional systems, stability cannot be determined solely by the slope of the $\dot{V} v s V .$ curve. Often, a local linear analysis can determine the stability of an equilibrium $(u_{0}, w_{0})$ , in 2 dimensions. We know from linear algebra that 2 nullclines can be approximated by their linearizations near the equilibrium. We use the Jacobian matrix, $A$ , at the equilibrium to obtain the eigenvalues and eigenvectors.

Now we will analyze the stability of the equilibrium point of $I_{N a, p} + I_{k}$ (Figure 7), which was approximately at $V = - 62 m V$ . The partial derivatives required for the elements of the Jacobian matrix are as follows,

\frac{\partial \dot{V}}{\partial V} = \frac{1}{C_{m}} (\begin{array}{l} - g_{L} - g_{N a} m_{\infty} (V) - \\ g_{N a} (V - E_{N a}) (- 1) \\ {(1 + e x p (V_{\frac{1}{2}} - V) / k)}^{- 2} \cdot \\ (e x p (V_{\frac{1}{2}} - V) / k) \cdot (\frac{- 1}{k}) - g_{k} \cdot n \end{array}),

(41)

\frac{\partial \dot{V}}{\partial n} = - g_{k} (V - E_{K}),

(42)

\begin{matrix} \frac{\partial \dot{n}}{\partial V} = \frac{1}{τ_{V}} {(1 + e x p (V_{\frac{1}{2}} - V) / k)}^{- 2} \\ (e x p (V_{\frac{1}{2}} - V) / k) \cdot (\frac{1}{k}) \end{matrix}

(43)

\frac{\partial \dot{n}}{\partial n} = - \frac{1}{τ_{V}}

(44)

The corresponding value of $n$ can be computed using the equilibrium value of $V$ in equation $(41)$ . The equilibrium coordinates are $(- 62, 0.032)$ The Jacobian matrix is as follows,

A = (\begin{matrix} - 0, 679 & - 280 \\ 0, 00625 & - 1 \end{matrix})

(45)

The trace of the matrix is $T r (A) = - 1, 6779$ and the determinant $| A | = 2, 429$ . Finally, we can calculate the eigenvalues for this equilibrium point using $\frac{T r (A) \pm \sqrt{{(T r (A))}^{2} - 4 | A |}}{2}$ . The eigenvalues for $V = - 62$ are $- 0, 84 \pm 1, 31 i$ . Since they have a negative real part, this corresponds to an equilibrium point with asymptotic stability; more precisely, to a stable spiral point.

Next, in Figure 8, we can observe trajectories in the phase plane with different initial values of continuous current inputs (1 and 10 $μ A$ ) but with different $B$ values (0.1, 1, 2, 5, and 10 mT), corresponding to the external magnetic field, while maintaining the same initial value for the potassium activation variable, $n$ . We can appreciate that all trajectories follow a subthreshold path toward equilibrium. We can see that this model does not have a fixed threshold for membrane potential. All action potentials shown are transient; they all return to equilibrium values.

Figure 8.

Representation of various action potential trajectories with different continuous current inputs and different $B$ values: (a) $I = 1 μ A$ and (b) $I = 10 μ A$ .

Figure 9 shows the emergence of the limit cycle between the values of $I = 100 - 400$ μA, corresponding to an external magnetic field $B = 5 m T$ . Limits cycles exhibit repetitive patterns of stationary movement in contrast to critical points, which represent equilibrium states. By definition, a solution $γ (t)$ is called periodic with period $p$ when it satisfies $γ (t + p) = γ (t)$ , for all $t$ , but $γ (t + m) \neq γ (t)$ , when $m < p$ . In 2-dimensional linear systems, the only possible periodic trajectories are centers, where 2 purely imaginary eigenvalues appear. However, the center is not considered an ideal model for many cyclic processes due to its structurally unstable nature; any variation in the system parameters can alter its dynamic structure, changing its behavior. Nevertheless, nonlinearity can lead to stable cyclic behaviors that persist despite small changes in the system parameters. These are limit cycles, phenomena exclusively associated with nonlinear systems. Finally, note that in the caption of Figures 8 and 9, the lines with other colors correspond to the same lines as in Figure 7.

Figure 9.

Representation of various action potential trajectories with different continuous current inputs and $B$ = 5 mT: (a) $I = 100 μ A$ , (b) $I = 200 μ A$ , (c) $I = 300 μ A$ , and (d) $I = 400 μ A$ .

STDP, or Spike Timing-Dependent Plasticity, is a phenomenon whereby the precise timing of neuronal spikes influences both the direction and magnitude of synaptic strength alterations.^53,54 Often regarded as the fundamental learning rule governing synapses, STDP has undergone successive refinements to address predictive limitations. In proposing an alternative fundamental principle, we suggest that exogenous or endogenous EM fields induce changes in crucial biochemical intermediates, such as neurotransmitters, serving as more direct triggers of plasticity mechanisms. The modulation of chemical synapses by STDP is dictated by the temporal discrepancy between pre- and postsynaptic neuronal firing, essentially the synchronization between neurons.⁵⁵ Concurrently, the magnetic coupling among neurons can enhance signal exchange and foster global synchronization within the neural network.^56,57

We have observed that the mechanism underlying the temporal information processing of a neuron is determined by a non-linear current-voltage relationship associated with the interactions between the exogenous EM field and certain proteins within the excitable membrane. Therefore, the variability of the spike timing of a neuron would be strongly influenced by extracellular EM fields in close proximity to a bifurcation point. However, we have demonstrated that within the range of $I = 0 - 400 μ A$ and $B = 0 - 10 m T$ , there exist equilibrium situations that would ensure precise control of the STDP mechanism.

Conclusions

Cognitive neuroscience has revealed a direct correlation between structural connections throughout the brain’s temporal hierarchical levels. Considering that all cognitive functions exhibit temporal characteristics, it is conceivable that the accurate estimation of transient alterations within specific temporal hierarchies is crucial for interpreting information from weak periodic endogenous or exogenous stimuli, such as static or extremely low-frequency magnetic fields. Thus, investigating the effects of magnetic fields on spike-time dynamics and their influence on other core cognitive functions may shed light on fundamental processes in typical cognition and their dysregulation in neuropsychological and mental disorders.

The EM field theory enables the determination of fields at the mesoscopic scale based on the morphological properties and bioelectrical tissue parameters at the microscopic scale. In this study, structural and functional intricacies are simplified by delineating 2 domains (internal and external neuron spaces), characterized by a set of interconnected fields governing the movement of charged particles within a magnetic field.

Estimating synaptic conductance for any noise correlation time remains an unresolved issue. While we outline a physical framework with simplifying assumptions to aid in understanding or predicting phenomena, obtaining data from experimental or quasi-experimental studies would be valuable and insightful.

We establish a connection between the STPD phenomenon and equations that depict the movement of neurotransmitters under an induced magnetic field. Our findings uncover the pivotal role of EM fields in governing the precise timing of pre- and postsynaptic spikes. This timing intricately shapes the outcomes of the STDP rule, while the propagation delay of neural signals significantly influences the connectivity patterns and dynamics of neural networks.^58,59

The macroscopic study of complex phenomena such as STPD is based on the application of constitutive equations to a supposedly continuous model subject to given constraints. These equations arise from the incorporation into conservation relationships (mass, momentum, energy. . .) of phenomenological laws known as kinetics, which connect the various currents therein (diffusion, chemical reaction rates, stresses, heat. . .) with corresponding state variables and their spatial derivatives (Fick’s law, chemical kinetics, Fourier’s law. . .).

In the particular case of STPD, the mechanisms generating instability are the autocatalytic transport of ions across excitable membranes in nerve cells or positive or negative feedback between neurons. Metabolic oscillations result from the control of enzyme activity, while oscillations in protein synthesis result from genetic regulation. The common point of biological rhythms in our study is that they have an exogenous origin. They possess all the properties of dissipative structures, in the sense that they exhibit a phenomenon of temporal self-organization under conditions of non-equilibrium. It is the internal regulatory mechanisms of biological systems that, through the bias of nonlinear kinetics, lead to instability phenomena associated with the emergence of periodic behaviors beyond the bifurcation point. In our case, within the physiological parameter range of intensity ( $I = 0 - 400 μ A$ ) and magnetic field ( $B = 0 - 10 m T$ ), the temporary dissipative structures (spike-timing) of action potentials associated with the STPD mechanism manifest rhythmic characteristics in the form of limit cycles around a stable steady state, resulting from nonlinear effects due to catalytic actions in associated chemical, physical, and biological processes.

Our approach to predicting timing changes in spike induced by magnetic fields can be extended to local neural networks. In fact, some authors consider extracellular EM field oscillations as a marker of neural network oscillations.^60
-62 It is even suggested that information contained in the oscillations of endogenous extracellular potentials could provide particularly relevant insights into the temporal coding of network activity.^63,64 This offers promising prospects regarding potential therapies for pathological neural synchronization, such as in the case of Parkinson’s disease, for example. The results may suggest that, by appropriately choosing the delay time, a mechanism of control like STPD through locally applied external EM field to a neural subpopulation could suppress global neuron synchronization. The coupling mechanisms of network activity to exogenous EM field are not an alternative to classical approaches but complementary mechanisms that take into account the spatiotemporal characteristics of the involved EM field.

A fresh perspective on the macroscopic system, arguably one of the most debated topics in neuroscience research, emerges from examining the average synaptic activity time in neural networks through the lens of path integrals. This approach offers insights into the partial understanding of asymptotic behavior on a microscopic scale. In the long run, it becomes imperative to delve into the topological equations governing neural temporal hierarchical levels, an area where analytical findings are scarce. This poses an issue of exceptional interest, as theoretical conjectures supported by experimental evidence indicate that neural networks retain information over timescales corresponding to their anatomical hierarchy.

Supplemental Material

sj-docx-1-bec-10.1177_11795972241272380 – Supplemental material for A Physical Framework to Study the Effect of Magnetic Fields on the Spike-Time Coding

Supplemental material, sj-docx-1-bec-10.1177_11795972241272380 for A Physical Framework to Study the Effect of Magnetic Fields on the Spike-Time Coding by Manuel Rivas and Marina Martinez-Garcia in Biomedical Engineering and Computational Biology

Footnotes

Author Contributions

M.R. and M.M. wrote the main manuscript text and prepared all the figures. All listed authors have contributed significantly to the manuscript and consent to their names and order of authorship on the manuscript. Co-authors will be kept informed of editorial decisions and changes made. There is no conflict of interests involving the conduction or the report of this research.

Declaration Of Conflicting Interests:

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding:

The author disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Founding grant CIGE/2022/066 from Generalitat Valenciana, Spain.

Availability of Data and Materials

Hereby, we M. R. and M. M. consciously assure that for this manuscript the following is fulfilled: (1) This material is the authors’ own original work, which has not been previously published elsewhere; (2) The paper is not currently being considered for publication elsewhere; (3) The paper reflects the authors’ own research and analysis in a truthful and complete manner; (4) The paper properly credits the meaningful contributions of co-authors and co-researchers; (5) The results are appropriately placed in the context of prior and existing research; (6) All sources used are properly disclosed (correct citation). Literally copying of text must be indicated as such by using quotation marks and giving proper reference, and (7) All authors have been personally and actively involved in substantial work leading to the paper and will take public responsibility for its content.

Ethical Approval

Not applicable.

ORCID iD

Marina Martinez-Garcia

Supplemental Material

Supplemental material for this article is available online.

References

Adey

WR.

Tissue interactions with nonionizing electromagnetic fields. Physiol Rev. 1981;61:435-514.

Tenforde

TS.

Biological interactions of extremely-low-frequency electric and magnetic fields. Bioelectrochem Bioenerg. 1991;25:1-17.

Schwan

. Biophysics of the interaction of electromagnetic energy with cells and membranes. In: Polk

Posto

eds. CRC Handbook of Biological Effects of Electromagnetic Fields. CRC; 1986:213-231.

Barnes

FS.

The effects of time varying magnetic fields on biological materials. IEEE Trans Magn. 1990;26:2092-2097.

Lawrence

Adey

WR.

Nonlinear wave mechanisms in interactions between excitable tissue and electromagnetic fields. Neurol Res. 1982;4:115-153.

Blackman

CF.

The biogical influences of low-frequency sinusoidal electromagnetic signals alone and superimposed on RF carrier waves. In: Chiabrera

Nicolini

Scham

eds. Interactions Between Electromagnetic Fields and Cells NATO Series. Plenum Press; 1985;521-536.

Liboff

McLeod

BR.

Kinetics of channelized membrane ions in magnetic fields. Bioelectromagnetics. 1988;9:39-51.

del Moral

Azanza

. Model for the effect of static magnetic fields on neurons. J Magn Mat. 1992;114:220-242.

Moggia

Chiabrera

Bianco

Fokker–Planck analysis of the Langevin–Lorentz equation: application to ligand-receptor binding under electromagnetic exposure. J Appl Phys. 1997;82:4669-4677.

10.

Hill TL, H.

Effect of rotation on the diffusion-controlled rate of ligand-protein association. Proc Natl Acad Sci USA. 1975;72:4918-4922.

11.

Rubinow

SI.

Introduction to Mathematical Biology. Wiley; 1975.

12.

Stanley

Gagner

Damanpour

, et al Radio-wave heating of iron oxide nanoparticles can regulate plasma glucose in mice. Science. 2012;336:604-608.

13.

Seo

Southard

Kim

, et al A mechanogenetic toolkit for interrogating cell signaling in space and time. Cell. 2017;169:1357-1518.

14.

Munshi

Qadri

Zhang

, et al Magnetothermal genetic deep brain stimulation of motor behaviors in awake, freely moving mice. eLife. 2017;6:1–26.

15.

Lee

Kim

Levy

, et al Magnetic nanoparticles for ultrafast mechanical control of inner ear hair cells. ACS Nano. 2014;8:6590-6598.

16.

Henstock

Rotherham

Rashidi

Shakesheff

El Haj

AJ.

Remotely activated mechanotransduction via magnetic nanoparticles promotes mineralization synergistically with bone morphogenetic protein 2: applications for injectable cell therapy. Stem Cells Transl Med. 2014;3:1363-1374.

17.

Tay

Kunze

Murray

Di Carlo

Induction of calcium influx in cortical neural networks by nanomagnetic forces. ACS Nano. 2016;10:2331-2341.

18.

Aarts

Liu

, et al Treatment of ischemic brain damage by perturbing NMDA receptor- PSD-95 protein interactions. Science. 2002;298:846-850.

19.

Galvanovskis

J.m

Sandblom

Berqvist

Hamnerius

Cytoplasmic ELF magnetic field intracellular Ca²⁺ oscillations. Proceedings of the Sixteenth Annual Meeting of the Bioelectrochemistry Society. 1999:90.

20.

Naundorf

Wolf

Volgushev

Unique features of action potential initiation in cortical neurons. Nature. 2006;440:1060-1063.

21.

Colwell

Brenner

MP.

Action potential initiation in the Hodgkin-Huxley model. PLoS Comput Biol. 2009;5:E1000265.

22.

Marom

Abbott

LF.

Modeling state-dependent inactivation of membrane currents. Biophys J. 1994;67:515-520.

23.

Mascagni

The backward Euler method for numerical solution of the Hodgkin–Huxley equations of nerve conduction. SIAM J Numer Anal. 1990;27:941-962.

24.

Clay

Paydarfar

Forger

DB.

A simple modification of the Hodgkin and Huxley equations explains type 3 excitability in squid giant axons. J R Soc Interface. 2008;5:1421-1428.

25.

Purcell

EM.

Electricity and Magnetism. Mc Graw-Hill; 1965.

26.

Tsang

Classical Electrodynamics. World Scientific; 1997:59.

27.

Sivasankar

Subramaniam

Leckband

Direct molecular level measurements of the electrostatic properties of a protein surface. Proc Natl Acad Sci USA. 1998;95:12961-12966.

28.

Zheng

Jensen

Savtchenko

, et al Nanoscale diffusion in the synaptic cleft and beyond measured with time-resolved fluorescence anisotropy imaging. Sci Rep. 2017;7:42022.

29.

Antosiewicz

Porschke

. The nature of protein dipole moments: experimental and calculated permanent dipole of .alpha.-chymotrypsin. Biochemistry. 1989;28:10072-10078.

30.

Mellor

Cortés

Busath

Mazzeo

BA.

Method for estimating the internal permittivity of proteins using dielectric spectroscopy. J Phys Chem B. 2011;115:2205-2013.

31.

Panofsky

WKH

Phillips

eds. Classical Electricity and Magnetism. 2nd ed. Dover Publications; 1962:333-338.

32.

Kholmetskii

Missevitch

Yarman

Electric/magnetic dipolein an electromagnetic field: force, torque and energy. Eur Phys J Plus. 2014;129:215-222.

33.

Friedman

Ito

Brinkman

, et al Universal critical dynamics in high resolution neuronal avalanche data. Phys Rev Lett. 2012;108:208102.

34.

Newell

Liu

Mayer-Kress

Time scales in motor learning and development. Psychol Rev. 2001;108:57-82.

35.

Grigolini

Rocco

West

BJ.

Fractional calculus as a macroscopic manifestation of randomness. Phys Rev E Stat Nonlin Soft Matter Phys. 1999;59:2603-2613.

36.

Zempel

Snyder

Raichle

ME.

The temporal structures and functional significance of scale-free brain activity. Neuron. 2010;66:353-369.

37.

Zuo

A geometric method to determine the electric field due to a uniformly charged line segment. Am J Phys. 2015;83:567-569.

38.

Woosley

Roth

Wikswo

JP.

The magnetic field of a single axon: a volume conductor model. Math Biosci. 1985;76:1-36.

39.

Grove

Ladas

Periodicities in Nonlinear Difference Equations. Chapman & Hall/CRC; 2004.

40.

Quintero-Quiroz

Pigolotti

Torrent

Masoller

Numerical and experimental study of the effects of noise on the permutation entropy. New J Phys. 2015;17:093002.

41.

Bandt

Pompe

Permutation entropy: a natural complexity measure for time series. Phys Rev Lett. 2002;88:174102.

42.

Rosso

Masoller

Detecting and quantifying temporal correlations in stochastic resonance via information theory measures. Eur Phys J B. 2009;69:37-43.

43.

Bizzarri

Cannistraro

Atomic force spectroscopy in biological complex formation: strategies and Perspectives. J Phys Chem B. 2009;113:16449-16464.

44.

Bell

GI.

Models for the specific adhesion of cell to cells: a theoretical framework for adhesion mediated by reversible bonds between cell surface molecules. Science. 1978;200:618-627.

45.

Pecht

Chemical Relaxation in Molecular Biology. Eds Springer-Verlag; 1976.

46.

Destexche

Rudolph-Lilith

Chapter 9. Synaptic noise experiments, computational consequences and methodsto analyseexperiments. In: Laing

Lord

eds. Stochastic Method in Neuroscience. Oxford University Press; 2009:242–271.

47.

Rivas

Martinez-García

The effect if exogenously induced magnetic fields on neurotransmitter dynamics. J Neurol Disord. 2022;10:2-9.

48.

Azanza

Calvo

del Moral

Evidence of synchronization of neuronal activity of molluscan brain ganglia induced by alternating 50 Hz applied magnetic field. Electromagn Biol Med. 2002;21:209-220.

49.

Azanza

del Moral

Cell membrane biochemistry and neurobiological approach to biomagnetism. Prog Neurobiol. 1994;44:517-601.

50.

del Moral

Azanza

Calvo

Pérez Bruzón

. Cooperative diamagnetism and Ca²⁺ liberation of plasma membrane molecules explains the neuron response to applied static and extremely low frequency magnetic fields. Proceedings of the 2nd International Workshop: Biological Effects of EMFs, Rhodas. 2002:651-659.

51.

Anastassiou

Montgomery

Barahona

Buzsáki

Koch

The effect of spatially inhomogeneous extracellular electric fields on neurons. J Neurosci. 2010;30:1925-1936.

52.

Anastassiou

Perin

Markram

Koch

Ephaptic coupling of cortical neurons. Nat Neurosci. 2011;14:217-223.

53.

Bell

Han

Sugawara

Grant

Synaptic plasticity in a cerebellum-like structure depends on temporal order. Nature. 1997;387:278-281.

54.

Markram

Lübke

Frotscher

Sakmann

Regulation of synaptic efficacy by coincidence of postsynaptic APs and EPSPs. Science. 1997;275:213-215.

55.

Takeuchi

Duszkiewicz

Morris

RG.

The synaptic plasticity and memory hypothesis: encoding, storage and persistence. Phil Trans R Soc A. 2014;369:20130288.

56.

Jia

Hayat

Alsaedi

Collective responses in electrical activities of neurons under field coupling. Sci Rep. 2018;8:1349.

57.

Guo

Wang

, et al Collective response, synapse coupling and field coupling in neuronal network. Chaos Solitons Fractals. 2017;105:120-127.

58.

Madadi Asl

Valizadeh

Tass

. Dendritic and axonal propagation delays determine emergent structures of neuronal networks with plastic synapses. Sci Rep. 2017;7:39682.

59.

Madadi Asl

Valizadeh

Tass

. Delay-induced multistability and loop formation in neuronal networks with spike-timing-dependent plasticity. Sci Rep. 2018;8:12068.

60.

Csicsvari

Jamieson

Wise

Buzsáki

Mechanisms of gamma oscillations in the hippocampus of the behaving rat. Neuron. 2003;37:311-322.

61.

Nunez

PL.

Electric Fields of the Brain: The Neurophysics of EEG. Oxford University Press; 2006.

62.

Sarnthein

Jeanmonod

High thalamocortical theta coherence in patients with Parkinson's disease. J Neurosci. 2007;27:124-131.

63.

Scherberger

Jarvis

Andersen

RA.

Cortical local field potential encodes movement intentions in the posterior parietal cortex. Neuron. 2005;46:347-354.

64.

Wang

Local field potential spectral tuning in motor cortex during reaching. IEEE Trans Rehabil Eng. 2006;14:180–183.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.02 MB