Abstract
It is widely accepted that synaptic transmissions and gap junctions are the major governing mechanisms for signal traveling in the neural system. Yet, a group of neural waves, either physiological or pathological, share the same speed of ∼0.1 m/s without synaptic transmission or gap junctions, and this speed is not consistent with axonal conduction or ionic diffusion. The only explanation left is an electrical field effect. We tested the hypothesis that endogenous electric fields are sufficient to explain the propagation with in silico and in vitro experiments. Simulation results show that field effects alone can indeed mediate propagation across layers of neurons with speeds of 0.12 ± 0.09 m/s with pathological kinetics, and 0.11 ± 0.03 m/s with physiologic kinetics, both generating weak field amplitudes of ∼2–6 mV/mm. Further, the model predicted that propagation speed values are inversely proportional to the cell-to-cell distances, but do not significantly change with extracellular resistivity, membrane capacitance, or membrane resistance. In vitro recordings in mice hippocampi produced similar speeds (0.10 ± 0.03 m/s) and field amplitudes (2.5–5 mV/mm), and by applying a blocking field, the propagation speed was greatly reduced. Finally, osmolarity experiments confirmed the model's prediction that cell-to-cell distance inversely affects propagation speed. Together, these results show that despite their weak amplitude, electric fields can be solely responsible for spike propagation at ∼0.1 m/s. This phenomenon could be important to explain the slow propagation of epileptic activity and other normal propagations at similar speeds.
SIGNIFICANCE STATEMENT Neural activity (waves or spikes) can propagate using well documented mechanisms such as synaptic transmission, gap junctions, or diffusion. However, the purpose of this paper is to provide an explanation for experimental data showing that neural signals can propagate by means other than synaptic transmission, gap junction, or diffusion. The results indicate that electric fields (ephaptic effects) are capable of mediating propagation of self-regenerating neural waves. This novel mechanism coupling cell-by-volume conduction could be involved in other types of propagating neural signals, such as slow-wave sleep, sharp hippocampal waves, theta waves, or seizures.
Introduction
Electrochemical signal propagation is one of the properties that enables communication among neurons. In the hippocampus, such communication is essential for basic neural function, such as memorization or spatial navigation, but can also cause pathological conditions, such as epilepsy. Interestingly, signal propagation speeds in various conditions are similar (∼0.1 m/s). Neural spikes generated by 4-aminopyridine (4-AP) travel with a longitudinal speed of 0.09 ± 0.03 m/s along the CA3 region (Kibler and Durand, 2011), whereas in the presence of picrotoxin, synchronous firing events propagate longitudinally at 0.14 ± 0.04 m/s (Miles et al., 1988). High K+-, low Mg2+-, and zero-Ca2+-triggered spikes again exhibit speeds of 0.07-0.1 m/s, 0.1–0.15 m/s, and 0.04–0.15 m/s, respectively (Haas and Jefferys, 1984; Quilichini et al., 2002; Liu et al., 2013). In normal tissue, theta oscillations travel with a speed of 0.08–0.107 m/s in the hippocampus of living rodent rats (Lubenov and Siapas, 2009), whereas carbachol-induced theta oscillations travels with a speed of 0.119 m/s along the CA1 cell layer and a 0.141 m/s along the CA3 cell layer (Cappaert et al., 2009). Together, it is clear that 0.1 m/s is a common propagation speed regardless of experimental models, therefore a shared fundamental mechanism may underlie these neural propagations.
Despite the fact that synaptic transmission and gap junctions are generally known to be responsible for neural communication, 4-AP-induced spikes still persisted when synapses and gap junctions were blocked in vitro, essentially suggesting a non-synaptic governing mechanism. These spikes propagate at the same or even higher speed in the presence of synaptic-blockers, such as the 4-AP/low-Ca2+ aCSF condition (Haas and Jefferys, 1984; Bikson et al., 1999; Zhang et al., 2014) and gap junction blockers, such as the 4-AP/mefloquine aCSF condition (Zhang et al., 2014). Other propagation mechanisms, such as extracellular ionic transients and axonal conduction mechanisms, have very different propagation speeds (0.0004–0.008 m/s for K+ diffusion and 0.3–0.5 m/s for axonal conduction, Miles et al., 1988; Lian et al., 2001; Francis et al., 2003; Meeks and Mennerick, 2007; Jensen, 2008; Kibler et al., 2012). Thus, it is unlikely that neural propagation at ∼0.1 m/s can be explained by synaptic transmission, gap junction, diffusion, or axonal propagation.
Another mechanism is electric field effects. Experiments have shown that extracellular fields could modulate the activity of single neurons and/or network activity. In hippocampal pyramidal cells, weak electrical fields can affect the excitability of cells and the synchronization of the network (Francis et al., 2003). In cortical neurons, endogenous weak fields could entrain action leading to phase-locking of spikes to the external field to reach a “spike-field synchrony” (Fröhlich and McCormick, 2010; Anastassiou et al., 2011) and its effect is enhanced by neural packing density (Vigmond et al., 1997). Using the Wilson and Cowan (1973) model for traveling waves, the speed of propagation is determined by excitability thresholds, spatial connectivity, and the refractory time constant (Richard et al., 1995). Locally applied electrical fields can also modulate the speed of propagation (Richard et al., 1995). Thus, we tested the hypothesis that an electrical field effect could underlie the propagation of 4-AP-induced spikes at a speed of 0.1 m/s in the hippocampus using computational modeling and experimental approaches. The propagation speed and network field amplitude were measured in two simulated conditions: epileptiform and normal, and the simulated results were compared with the propagation speeds and field strengths in in vitro experiments. Further, in vitro experiments were done to study the effect of applying a blocking field on propagation. Results from the model also predict that cell-to-cell distance (extracellular space) could inversely affect speed due to field effect, and this prediction was tested in vitro with osmolarity experiments.
Materials and Methods
Computational methods
Model structure.
The compartment model was built under NEURON 7.3 simulation environment (Yale School of Medicine, New Haven, CT). A CA1 hippocampal pyramidal network was simulated because it is a typical region where 4-AP-induced spike was observed experimentally. Two simulation models were tested: 4-AP/CA2+-free condition for epileptiform propagation and normal aCSF condition for physiological propagation. In the model, each cell contains 33 compartments, with one compartment for soma, 21 compartments for the apical dendrite, and 11 compartments for the basal dendrite. The lengths of the apical and basal dendrites were 735.3 and 490.2 μm, respectively, whereas both had a diameter of 5.2 μm (Shuai et al., 2003). Similar as to previous experiments and models, the diameter of each cell body was 10 μm (Warman et al., 1994; Shuai et al., 2003; Kurt et al., 2004; Anastassiou et al., 2011). The transmembrane potential Vi for each compartment is described by the relationship known as cable equation. where Cm,i, Vi, Im,i, are the membrane capacitance, transmembrane potential, and transmembrane current of the i,th compartment respectively; ri,i+1 is the cytoplasmic (axial) conductivity between the i,th and (i +1)th compartment. To simplify the model, only the soma compartment contains active conductance; for dendritic compartments, the transmembrane current only contains leak current Im,idend = Ileak,i = gi Vi, with gi representing the passive conductance of the i,th compartment.
All passive membrane properties for each pyramidal cell were set to the following values according to patch-clamp recordings and current-clamp data in Ca2+-free conditions by Shuai et al. (2003): somatic membrane resistivity Rmsoma = 680 Ωcm2, somatic membrane capacitance Cmsoma = 1 μF/cm, dendrite membrane resistivity Rmdend = 34,200 Ωcm2, dendritic capacitance Cmdend = 1 μF/cm, and cytoplasmic (axial) resistivity Ri 530Ωcm. For each cell in the model, the transmembrane current for somatic compartments Im,isoma are given by the following: where ILeak,i and IIons,i are the leak current and total active ionic currents, respectively.
For 4-AP/Ca2+-free kinetics, the calcium currents and calcium-dependent potassium currents are eliminated (Shuai et al., 2003). In vitro experiments and computational models for CA1 hippocampal pyramidal cells stated that there are generally five major active ionic conductances in a Ca2+-free extracellular environment: a fast sodium current INa, a delayed rectifier current IDR, an A-type transient current IA, a muscarinic current IM, and a persistent sodium current INaP (Warman et al., 1994; Shuai et al., 2003). Thus, the total ionic current in the somatic compartment consists of five active currents: The gating equations for each active current are implemented using the Hodgkin-Huxley formalism (Hodgkin and Huxley, 1952) as in Table 1. The maximum conductances are based on those of Shuai et al. (2003), designed to reproduce the electrophysiological behavior of CA1 pyramidal cells in low Ca2+ concentration environment. However, 4-AP is a blocker of the Kv1 voltage-activated K+ (A-type transient K) channels; thus, we greatly reduced the maximum conductance of the IA potassium channel, gA, to a near-zero value for 4-AP/Ca2+-free kinetics, Table 2 shows the maximum conductance values used in the model.
To test whether field effect could mediate propagation with physiologic channel kinetics (normal aCSF), we adapted an established computational model by Wimmer et al. (2010) and Royeck et al. (2008). With these kinetics, each cell contained six different distributed Ca2+ channels, three different Na+ channels, and seven K+ channels that essentially mimicked a normal experimental condition. The morphology of the cell, the network layout, and neural activity initiation were the same as the 4-AP/Ca2+-free model.
The neural network.
The layout of the network was designed in a way to capture the most important features of the physiological cellular layout for hippocampal activity propagation while minimizing the number of elements. The network is made up of three rows (Rows A–C) of neurons with 10 cell locations in each row (Fig. 1a,b). The distance between each two adjacent neurons (dc-c) in a physiologic pyramidal network ranges from 2 to 5 μm (Kurt et al., 2004; Anastassiou et al., 2011); in the model, we tested 10 trials for each of the mean dc-c equal to 2, 3, and 4 μm (n = 10 and variance = 0.1 μm with Gaussian distribution for each dc-c). To represent the cell stacking throughout the depth of a tissue slice, an arbitrary “stacking factor” (SF) was generated to take into account the actual number of cells around a certain location (Fig. 1c). The value of the SF was estimated from the volume cellular density reported from previous histology studies. In mice hippocampus, CA1 cellular volume density varies from ∼180,000 neurons/mm3 to 300,000 neurons/mm3 depending on specific locations in the CA1 region (Kurt et al., 2004; Jinno and Kosaka, 2010). Using a tissue thickness of ∼600 μm and the model network area equal to ∼4216–4788 μm2 (Fig. 1a; d̄ = 2 μm, area = 34 × 124 μm; d̄ = 4 μm, area = 38 × 136 μm), a SF between 15 and 28 would best represent the cellular layout in tissue. In the model, stacking factor values between 10 and 30 were tested.
Electrical field effect calculation.
To test the hypothesis that endogenous electric field alone could induce neural propagation, the synaptic connections, gap junctions, and diffusion effects were eliminated in the network. The cells in the first row were first initiated by an experimentally determined extracellular potential (Fig. 2a). The mean values of 10 in vitro signals acquired from each of the channels located at cell layer, apical dendrites, and basal dendrites by Michigan Shank were applied extracellularly at corresponding locations in each cell of Row A in the model (Fig. 2b). The electric field effect was calculated using the quasi-static formulation of the Maxwell equations assuming homogenous and linear volume conductors (He, 2005). According to Ohm's law for conductors in Maxwell's equations, the potential ϕ at the point P at a distance r relatively to the reference electrode in a medium of conductivity σ is as follows: Using superposition, Equation 1 can be generalized to j monopolar electrodes from a total of i cells from the previous row(s) with a current Ii,j→n located at a distance ri,j→n to each target node n of the target cell. The extracellular potential inserted at a target node n is given by the following: where Vex,n is the extracellular potential inserted at target node n in the target cell, Ii,j→n is the membrane current (assuming current going out of the cell to be positive direction) at a node j of cell i located ri,j→n distance from target node n, with a total j = 35 nodes in each cell (Fig. 2b); σ is the extracellular conductivity (typically the extracellular resistivity ρ = 1/σ is 250–380 Ωcm; 300 Ωcm was used in the model; Gold et al., 2007; Logothetis et al., 2007); SF is the stacking factor. The electric field effect could only propagate the action potential in the transverse direction to simplify the model; Row B cells were only affected by Row A cells, whereas Row C cells were affected by both Row A and Row B cells with the electric field effect.
Field amplitude and speed measurement.
To measure the resulting electrical field due to the network firing activity and compare it to the experimentally recorded waveforms, we placed three virtual electrodes (v1, v2, and v3) outside of the network (Fig. 3a) vertically to test the network field amplitude, similar as to the placement of Michigan Shank electrodes into the in vitro slice. The middle electrode was placed 30 μm away from the Row C middle cell to account for approximately three rows of dead cell around the Michigan Shank electrode, a situation observed experimentally. Each row's resulting extracellular voltage at the three virtual electrodes was calculated, and the field amplitude was calculated by finding the average of the spatial derivatives (Eq. 3). The network field amplitude was the summation of the field amplitudes generated by all three Row layers (Eq. 4). The propagation speed was measured based on the intracellular recording of each row's middle cell, where the first spike of the action potential peak time from each cell was recorded, and the delay from Row A to Row B and from Row B to Row C were calculated. The propagation speed was derived by taking the traveling distance divided by delay time.
Experimental methods
In vitro field amplitude measurements in unfolded hippocampus under the 4-AP/low Ca2+ condition.
In vitro measurements of extracellular activity and electric field in unfolded mice hippocampus were performed to measure the propagation speed and validate the results of the simulation. Experimental protocols were approved by the Institutional Animal Case and Usage Committee. CD1 mice from Charles River Laboratories (P10–P20; of either sex) were anesthetized and then decapitated. Subsequently, the whole hippocampus was separated from the brain and unfolded by following a previously developed protocol (Kibler and Durand, 2011). All unfolded hippocampal preparations were then transferred into a recovery solution of aCSF bubbled with 95% O2/5% CO2 and maintained at a room temperature (25°C) for at least 1 h before recording. Normal aCSF buffer consisted of the following (mm): 124 NaCl, 3.75 KCl, 1.25 KH2PO4, 2 MgSO4, 26 NaHCO2, 10 dextrose, and 2 CaCl2.
During recording, 4-AP was added into a low-calcium solution to increase the excitability of the tissue (Perreault and Avoli, 1989, 1992; Schechter, 1997) and to maximally block synapses. Low Ca2+/4-AP aCSF consisted of the following (mm): 0.1 4-AP, 124 NaCl, 5.25 KCl, 1.25 KH2PO4, 1.5 MgSO4, 26 NaHCO3, 10 dextrose, and 0.2 CaCl2. The extracellular potentials were acquired by a Michigan-type electrode array inserted into the unfolded hippocampus along the dendritic–soma axis. Electrical field amplitudes were calculated as the first difference of voltages measured by two adjacent channels on the shank, one at the depth of cell layer (200 μm from the stratum) and the other in apical dendrites. The Michigan shank was combined through a multichannel workstation with data acquisition system by Alpha ω (AlphaLab SNR).
In vitro propagation speed measurements in unfolded hippocampus under the 4-AP/low Ca2+ condition.
The propagation of the extracellular activity and its speed in the CA1 region of unfolded hippocampus was recorded by a custom-made 64-channel microelectrode array as previously described (Kibler et al., 2012; Zhang et al., 2014). The preparation for the tissue and the solution (4-AP/low-Ca2+) was the same as described above. Each microelectrode was coated with a biocompatible insulation material with only the tip being exposed. The electrodes were patterned as an 8 × 8 matrix with a 400 μm electrode-to-electrode space on the transverse direction and 300 μm space on the longitudinal direction. As the unfolded hippocampus was placed on the array, each microelectrode was inserted into the tissue and the extracellular activity at the location was recorded. The propagation speed was then calculated based on the distance between the two neighboring electrodes divided by the spike time delay recorded by the electrodes (Kibler et al., 2012; Zhang et al., 2014).
In vitro propagation speed reduction by applying a blocking field in longitudinal slices under the 4-AP condition.
Preliminary studies were conducted to observe the effect of applying an external blocking electric field on propagation of spontaneous activity in vitro. The 4-AP solution consisted of the following (mm): 0.1 4-AP, 124 NaCl, 3.75 KCl, 1.25 KH2PO4, 2 MgSO4, 26 NaHCO2, 10 dextrose, and 2.6 CaCl2. Longitudinal hippocampal slices were used to observe the effect of applying an external electric field in the opposite phase (“anti” field) to that of the propagating spike in between the two recording sites on the CA3 layer (see Fig. 6a). A stimulator and stimulus isolation unit was used to generate a field between 2 and 5 mV/mm. The electric field was generated across two 30 μm diameter tungsten-stimulating electrodes with (World Precision Instruments) that were placed across the tissue and between the two recording sites. 4-AP-induced spikes were generated when tissue was exposed to 4-AP/aCSF solution. Upon determining the propagation direction, the field strength was varied and the delay between the two recording channels was calculated.
In vitro osmolarity test on its effect of speed in longitudinal slices under the 4-AP condition.
To determine how extracellular space (osmolarity) affects the propagation speed and thus validating the role of the electrical field effect, we performed osmolarity experiments in vitro on longitudinal slices in 4-AP solution, as the propagation is usually very robust in that direction (Zhang et al., 2014). The extracellular osmolarity was both increased and decreased relative to the normal osmolarity condition. The control 4-AP solution was the same as described above. For the low-osmolarity experiment, 150 ml of deionized water and 15 μm of 4-AP were added into the control 4-AP solution, resulting in a 15% dilution and a 10% decrease in osmolarity. For high osmolarity experiment, 30 mm d-mannitol was added to the control 4-AP solution, causing an increase in osmolarity by 10%. Each osmolarity experiment (control, low, and high) started with the tissue being exposed to the control 4-AP solution for 30 min and field potential recorded from each trial. The high- or low-osmolarity solution was then introduced, and the field potential was recorded after 15 min of osmolarity change. The traveling speed of the spontaneous activity was measured by two recording glass micropipette electrodes (1.0 mm outer diameter and 0.5 mm inner diameter borosilicate glass filled with 150 mm NaCl solution) placed on the CA3 longitudinal hippocampal slices, and the speed was determined by the transit time between the electrode pair.
Results
Neural activity can propagate non-synaptically with weak field amplitudes (pathological propagation)
We first sought to constrain and validate the simulated pathological model by recording in vitro electric field amplitude during epileptiform activity (4-AP/low-Ca2+ solutions) in the hippocampus. Extracellular recordings obtained in the somatic layer (Va) and apical dendrites (Vb) by Michigan Shank showed trains of bursting events recorded under these conditions (Fig. 4a, Left). The electrical field amplitude near the cell layer was estimated by taking the spatial derivatives of the two voltage measurements. Among 65 spikes in five unfolded hippocampi, the field spiking amplitude was found to be between 1.60 and 7.13 mV/mm with a mean value of 3.5 ± 1.0 mV/mm; most of the amplitudes fell between 2.5 and 5 mV/mm (Fig. 4a, right). We then compared these signals with the field amplitude obtained in silico with 4-AP/Ca2+-free kinetics and observed similar characteristics (Fig. 4b, left). The median network field amplitude from the model was 4.58 mV/mm (Fig. 4b, left), and the field ranged between 3.26 and 5.80 mV/mm (Fig. 4b, right) among 180 simulated trials with randomized cell-to-cell distance dc-c (2–4 μm) and physiologically relevant SF values based on hippocampal cell density. These values fall within the known field amplitude range of ∼3–6 mV/mm (Ghai et al., 2000; Radman et al., 2007; Fröhlich and McCormick, 2010). It was also determined that the field amplitude increased with stacking factor SF (Fig. 4b, right). The extracellular activity generated in the model traveled from Row A to Row B and from Row B to Row C, showing that electrical field effect alone was responsible for action potential propagation in the CA1 hippocampal network (Fig. 4b, left, inset).
In addition, the intracellular waveform generated in silico was similar to those observed experimentally. Using the control parameters given in Table 2, the model reproduced the essential characteristics of the intracellular bursting activity seen in CA1 pyramidal cells. The resting potential, first spike amplitude, and first spike width at half-amplitude width were −61.5 mV, 93.9 mV, and 6.5 ms, respectively (Fig. 4c). These parameters fall within the experimentally measured ranges from previous CA1 pyramidal cell in vitro studies (resting potential: −75 mV∼ −52 mV; first spike amplitude: 70 ∼ 103.2 mV; and first spike width at half-amplitude: ∼1∼8 ms; Deuchars and Thomson, 1996; Su et al., 2001; Shuai et al., 2003; Golomb et al., 2006). The action potential fired at 12.09 mV above the resting potential when the stimulus was at 0.132 nA/10 ms, matching the in vivo measured threshold values (12.1 mV above the resting potential; Su et al., 2001; Fig. 4c).
Spiking events propagate by electrical field effect with speeds of ∼0.1 m/s in silico and in vitro (pathological propagation)
We then looked at the recorded epileptiform activity propagation speed in the CA1 region of unfolded hippocampus under 4-AP/low-Ca2+ environment. A high aspect ratio microelectrode array was used to record the activity and the spikes obtained in a solution containing 4-AP/low-Ca2+ at four adjacent spikes are shown in Figure 5a traveling at a speed of ∼0.1 m/s. Among 139 speed samples from the CA1 regions of six unfolded hippocampal preparations, the in vitro transverse speed ranged from 0.05 to 0.2 m/s, with a mean of 0.10 ± 0.03 m/s (Fig. 5b). These results clearly match the values obtained in the simulations, and are similar to those obtained in Zhang et al. (2014).
The experimental propagation speed in 4-AP/low-Ca2+ was compared with the simulated propagation speeds with the 4-AP/Ca2+-free kinetics in silico. The propagation speed across the network was evaluated based on the action potential peak time delay from each row's middle cell (namely, cell 5 of each row; Fig. 1a). This is because the middle cells experience the strongest field effect compared with other cells in the row according to Equations 1 and 2. We recorded action potential waveforms and calculated the propagation speeds by taking the distance traveled from Row A to Row C divided by the total action potential delay time (Δt1 + Δt2). Figure 5c shows an example of membrane potential signals with dc-c of 2.94 μm and SF set to 20 in the middle of the physiological range. Through only the electrical field effect in the network, the activity propagated from Row A to Row C with propagation speed of 0.093 m/s. Among 180 simulated trials with random cell-to-cell distance dc-c and physiological stacking factor SF values, the field effect driven activity in the model propagated at a mean speed of 0.23 ± 0.096 m/s, 0.10 ± 0.048 m/s, and 0.05 ± 0.020 m/s for trials with mean dc-c = 2, 3, and 4 μm, respectively, with a total mean speed of 0.12 ± 0.097 m/s; the propagation speed increased with higher values of SF (Fig. 5d).
Blocking the endogenous field (in vitro) and increasing cellular distance (in silico and in vitro) both result in a decrease in speed
To determine whether electric fields directly affect the speed of propagation, two electrodes were positioned to generate a DC electric field parallel to the dendritic tree main axis thereby opposing the endogenous field (Fig. 6b). Without stimulation, the average delay between two events obtained from the two recording electrodes was 66.01 ± 1.01 ms (42 spike-pair, 3 slices). With a field of 1 mV/mm applied across the cell layer the average delay increased to 74.86 ± 0.99 ms (40 spike-pair, 2 slices), resulting in a 13% decrease of the speed (Fig. 6c). Varying the field to 3 and 5 mV/mm resulted in 81.85 ± 1.05 ms (35 spike-pair, 2 slices) and 121.58 ± 15.67 ms (34 spike-pair, 2 slices) in delay time, respectively. This corresponds to 24% (3 mV/mm) and 84% (5 mV/mm) decreases in speed. These results confirm that the propagation mechanism of spontaneous activity in the hippocampus is consistent with electric fields.
The in silico model predicted that the speed of propagation should decrease with increasing distance between neurons. When setting dc-c = 3 μm as the control group, the propagation speed increased by 130% when dc-c decreased to ∼2 μm, and decreased by 50% when mean dc-c increased to ∼4 μm (Fig. 6d). To test this prediction, we subsequently performed osmolarity experiments on longitudinal slices in 4-AP solution in vitro and found a similar relationship between the change in extracellular space volume and the resulting 4-AP-induced spike propagation speed. By decreasing the osmolarity of 4-AP aCSF (cell swelling and dc-c decreases), the time delay between signals measured along the CA3 layer decreased from 10.12 ± 3.22 ms (n = 41 spike-pairs from 2 slices in normal osmolarity) to 4.53 ± 2.59 ms (n = 73 spike-pairs from the same 2 slices in low osmolarity >10 min after dilution). With a measured distance between the two glass pipettes equal to 0.8 mm, a 125% increase in speed was obtained. Conversely, by increasing the osmolarity of 4-AP aCSF (cell shrinkage and dc-c increases), the time delay increased from 5.95 ± 1.12 ms (26 spike-pair from 2 slices in normal osmolarity) to 9.76 ± 1.78 ms (50 spike-pair from the same 2 slices in high osmolarity), corresponding to a 46.8% decrease in speed (Fig. 6e). These results confirm the predictions of the model regarding the effect of cell-to-cell distance on the speed of propagation.
Electrical field propagation of 0.1 m/s also occurs in physiological conditions in silico
Previous experimental studies suggested that physiological propagations, such as theta wave also exhibits a traveling speed of 0.1 m/s (Cappaert et al., 2009; Lubenov and Siapas, 2009). With normal aCSF kinetics, the modeled cell produced an intracellular signal that captured physiological action potential characteristics upon initial stimulation (Fig. 7a). The speed and network field amplitude were collected with random dc-c centered at 2 μm and with physiological SF range (N = 40: n = 10 trials with random cell-to-cell distance for mean dc-c = 2 μm and variance = 0.1 μm and physiological SF = 15, 20, 25, and 30). Within the physiological range of SF, the mean propagation speed was 0.11 ± 0.03 m/s (Fig. 7c,e) and the field amplitude ranged from 2.33 to ∼3.06 mV/mm (Fig. 7b,d). These results imply that electrical field effect alone can produce neural propagation not only in pathologic epileptiform-inducing environments, but also in networks with normal channel kinetics.
Parameter sensitivity analysis indicates that the model is consistent with known applied field effects and robust to changes in intrinsic parameters (pathologic condition)
Previous studies related to the effect of an externally applied field on the speed of a traveling wave both in individual hippocampal CA1 pyramidal cells and a network of neocortical cells showed that an excitatory (or depolarizing) field can enhance the propagation speed (Richardson et al., 2005; Radman et al., 2007). Our results confirmed this relationship but with an endogenous electrical field. For both pathologic and physiologic conditions, as the field amplitude increased due to the increasing SF, the resulting speed also increased (Figs. 4b, right, 5d, for the pathologic condition; Fig. 7d,e for the physiological condition). Further, we found an inversely-linear relationship between the field amplitude and the propagation time delay for every cell-to-cell distance group (Fig. 8a), similar to the inversely-linear relationship between applied field and changes in action potential delay in a single cell shown by Radman et al. (2007). Changing the extracellular resistivity (Rex) also affected the field amplitude thus the propagation speed. We examined the effect on the propagation speed with ±40% change of the Rex value (n = 3 with dc-c = 2 μm, 3 μm, and 4 μm for each Rex). The speed increased with a larger value of Rex (Fig. 8b) because of the resulting higher field amplitude (Fig. 8c) within the physiologic range, but the effect was not significant (ANOVA one-way analysis and one-sample t test) for field amplitude or speed. These results indicate that the field effect does not depend on one important variable involving the field amplitude, namely the extracellular resistance within the physiological range of SF values. We also examined the dependence on membrane-specific resistance and capacitance (data not shown) and found that the effect is robust to changes in those variables as well.
Discussion
The observation that neural activity can propagate at a speed of ∼0.1 m/s in both the physiologic and pathologic conditions (Haas and Jefferys, 1984; Miles et al., 1988; Quilichini et al., 2002; Cappaert et al., 2009; Lubenov and Siapas, 2009; Kibler and Durand, 2011; Liu et al., 2013) suggests that a common mechanism for neural communication is involved. Despite the fact that chemical synapses and gap junctions usually maintain propagation, 4-AP-induced spikes still persisted in the presence or absence of chemical/electrical synapses (Zhang et al., 2014), thus, the mechanism could not be synaptic or gap junction mediated. Previous studies related to such type of propagation also eliminated the possibility for axonal conductions and ionic diffusion (Miles et al., 1988; Weissinger et al., 2000; Lian et al., 2001; Meeks and Mennerick, 2007; Jensen, 2008; Kibler and Durand, 2011) as their speeds are not consistent with a speed of 0.1 m/s. The only remaining explanation for the governing mechanism is endogenous electrical field effect. The field effect is said to be endogenous when no external current or voltage source generate the field; rather, it is the electrical field generated by the cells firing themselves that in turn influence the behavior of the adjacent cells in the same network. The electrical field effect is generally thought to be too small to give rise to action potentials (Chan and Nicholson, 1986; Anastassiou et al., 2011). However, several studies have shown that electrical fields applied externally can modulate neural activity at low amplitudes or with amplitudes that match the endogenous levels (Chan and Nicholson, 1986; Radman et al., 2007; Fröhlich and McCormick, 2010; Anastassiou et al., 2011). Francis et al. (2003) demonstrated that a neural network is much more sensitive to field modulation than single neurons, and that the highly structured and dense cellular packing in hippocampus makes this region of the brain highly sensitive to field effects. Moreover, endogenous electrical fields generated by antidromic activation can trigger after-discharges in the hippocampus in the absence of synaptic transmission, indicating that transient extracellular electrical fields (filed effect interaction) could be responsible for the observed synchrony (Taylor and Dudek, 1982). However, there is little or no information in the literature concerning the ability of electric field to sustain a regenerative wave propagating in densely packed neural network. Yet, the propagation of 4-AP-induced spikes (at 0.1 m/s) recently observed in the hippocampus in the absence of synapses and gap junctions was attributed to field effects by default (Zhang et al., 2014).
In this study, we used computer simulations in silico to test the feasibility of field effect as a propagation mechanism of action potentials within a hippocampal CA1 pyramidal network and performed experiments in vitro to confirm both speed and field amplitude results. Further, by applying a blocking “anti” field into the network in vitro, we observed speed reduction. The simulation in silico also provided information on how endogenous field strength (amplitude) relates to the propagation speed, which is consistent with previous experimental studies, as well as how robust the propagation speed is in regard to internal cell and network parameters. Finally, the model made predictions that were confirmed experimentally concerning the relation between cell distance and speed, and that the speed was not related to extracellular resistivity, membrane capacitance, or membrane resistance (in silico only). Together, this study confirms that the electric field effect with amplitude within the physiologic range can be sufficient to generate propagation of neural activity. Furthermore, we provided a demonstration that the speed of propagation is determined by the amplitude of the field.
The cellular network in the simulation consists of 30 hippocampal CA1 pyramidal cell locations (10 cells/row × 3 rows; Fig. 1), where cell-to-cell distances were set between ∼2 and ∼4 μm (Kurt et al., 2004; Anastassiou et al., 2011). An arbitrary SF was used to take into account the actual number of cells around each location. Its value was constrained by the physiological cell density in hippocampal CA1 region (physiological SF range = 15∼28 based on cell density). The model was validated as the field amplitude from the simulation with pathologic kinetics (∼3 to ∼6 mV/mm; Fig. 4b) and physiologic kinetics (∼2 to ∼3 mV/mm; Fig. 7d) matched the range obtained in our experimental recordings in hippocampus (2.5–5 mV/mm; Fig. 4a) and in the literature of recorded field amplitude in hippocampal slices (Ghai et al., 2000; Radman et al., 2007; ∼5 mV/mm) and in cortex (up to 4.52 mV/mm; Fröhlich and McCormick, 2010). Finally, the model successfully reproduced physiological intracellular spike characteristics both in a 4-AP/Ca2+-free (Fig. 4c) and in normal aCSF environment (Fig. 7a). These results indicate that the model is well validated to simulate a physiological system composed of layers of neurons.
Upon initial stimulation, simulated transverse propagation of the activity was observed both with 4-AP/Ca2+-free and normal channel kinetics in the simulation at an average speed of ∼0.1 m/s (Figs. 5c,d, 7c,e), matching the 0.1 m/s speed reported for 4-AP-induced spike propagation (Miles et al., 1988; Quilichini et al., 2002; Kibler and Durand, 2011; Liu et al., 2013), low Ca2+ field-bursting (Haas and Jefferys, 1984), and physiologic theta oscillations both in vitro (Cappaert et al., 2009) and in vivo (Lubenov and Siapas, 2009). Experiments in the CA1 region of unfolded hippocampus under 4-AP/low-Ca2+ condition in vitro also showed the same speed range (Fig. 5a,b). By applying an “anti” fields of matching amplitudes as the recorded field amplitudes in real time, a reduction of speed was seen (Fig. 6c). Together, both in silico and in vitro results indicate that field effect itself can be solely responsible for neural activity propagation with a speed of 0.1 m/s in a cellular network independently of channel kinetics.
Further, the electrical field effect is consistent with the results of the experiments on the effect of osmolarity (extracellular space) on neural propagation speed. In the study, the simulation predicted that a smaller dc-c could induce a higher speed (Fig. 6d). In vitro experiments confirmed this prediction, whereby 4-AP-induced epileptiform activity in hippocampal longitudinal slices traveled faster when extracellular osmolarity was lowered (cells swelling), and vice versa (Fig. 6e). This also matches with the relationship between osmolarity and propagation speed in low Mg2+-induced seizure activity in hippocampal CA1 region (Shahar et al., 2009). These results confirmed that the strength of electrical field plays a significant role in determining the neural activity propagation speed. For instance, a smaller dc-c increases the endogenous field amplitude according to Equations 1 and 2. Experiments from Francis et al. (2003) and Vigmond et al. (1997) also confirmed that a more tightly packed network increases the electrical impedance and field sensitivity between cells, thereby strengthening the field amplitude. On the other hand, parameter sensitivity analysis showed that the propagation speed is highly influenced by change in endogenous field amplitude but not by intrinsic parameters. It has been suggested that there exists a linear relationship between externally applied field amplitude and changes in action potential timing in a single cell such that the stronger the field strength, the earlier the cell fires (Radman et al., 2007). Interestingly, we found that the endogenous field amplitude and propagation delay share a similar inverse linear relationship (Fig. 8a), where the time delay decreases (speed increases) as field amplitude increases. Parameter sensitivity analysis showed that parameters that affect minimally the field amplitude, such as extracellular resistivity (Fig. 8b,c), membrane capacitance (data not shown), and passive membrane resistance (data not shown) in turn do not affect propagation speed. These results show that network endogenous field amplitude is the key factor in determining the propagation speed induced by field effect, and the speed is robust to extracellular resistivity, membrane capacitance, and membrane resistance. Therefore, the propagating mechanism is dependent on the network properties (eg, cell density) but not individual cell properties, suggesting that this phenomenon may be more general than originally anticipated.
Thus, the hypothesis that the propagation of neural activity can be carried out solely by electric fields is consistent with experimental data as well as computer simulations. The fact that the speed of the propagation remains constant under different experimental conditions can be explained by the presence of an electric field effect associated with neural firing generated by the network configuration and not by the properties of individual cells. This electrical field effect is revealed when synaptic transmission is blocked and both pathologic and normal propagation can exist simultaneously with field effect governing short range or local propagation, where synaptic transmission may govern long-range propagation and communication.
Conclusion
This study shows that in a simulated hippocampal pyramidal cell network, an electrical field of ∼2–5 mV/mm can underlie the transverse propagation of action potentials with a speed of ∼0.1 m/s, both with pathologic (4-AP/CA2+-free) and physiologic channel kinetics. This field-effect-induced propagation speed matches with the speed of 4-AP-induced activity, as well as theta waves recorded from both in vitro and previous in vivo experiments. Parameter sensitivity analysis revealed that the endogenous field amplitude is the main factor that plays an important role in influencing the speed (there exists an inverse linear relationship between field amplitude and propagation time delay), whereas intrinsic parameters, such as passive and active membrane properties, do not significantly affect the speed. Finally, it was shown both in experiments and simulations that a tightly packed network (either in silico or induced by low osmolarity) can produce an increase in propagation speed (and vice-versa), confirming that extracellular-space induced changes in field amplitude can influence speed. Moreover, applied blocking electric fields can significantly affect the speed of propagation. Together, these results reveal a novel propagation mechanism consistent with ephaptic field effects that could lead to better understanding of neural activity propagation in both the abnormal condition, such as epilepsy, or the normal condition, such as theta waves.
Footnotes
The experimental work was funded by National Institute of Health (National Institute of Neurological Disorders and Stroke) Grant 1R01NS060757-01. We thank Drs. Kenneth Gustafson, Hillel Chiel, Thomas Ladas, and Luis Gonzalez-Reyes for valuable comments and discussions.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dominique M. Durand, Department of Biomedical Engineering, Neural Engineering Center, Wickenden Building 112, 2071 Martin Luther King Jr Drive, Case Western Reserve University, Cleveland, OH 44106. dxd6{at}case.edu