- Open Access
Transmembrane potential induced on the internal organelle by a time-varying magnetic field: a model study
Journal of NeuroEngineering and Rehabilitation volume 7, Article number: 12 (2010)
When a cell is exposed to a time-varying magnetic field, this leads to an induced voltage on the cytoplasmic membrane, as well as on the membranes of the internal organelles, such as mitochondria. These potential changes in the organelles could have a significant impact on their functionality. However, a quantitative analysis on the magnetically-induced membrane potential on the internal organelles has not been performed.
Using a two-shell model, we provided the first analytical solution for the transmembrane potential in the organelle membrane induced by a time-varying magnetic field. We then analyzed factors that impact on the polarization of the organelle, including the frequency of the magnetic field, the presence of the outer cytoplasmic membrane, and electrical and geometrical parameters of the cytoplasmic membrane and the organelle membrane.
The amount of polarization in the organelle was less than its counterpart in the cytoplasmic membrane. This was largely due to the presence of the cell membrane, which "shielded" the internal organelle from excessive polarization by the field. Organelle polarization was largely dependent on the frequency of the magnetic field, and its polarization was not significant under the low frequency band used for transcranial magnetic stimulation (TMS). Both the properties of the cytoplasmic and the organelle membranes affect the polarization of the internal organelle in a frequency-dependent manner.
The work provided a theoretical framework and insights into factors affecting mitochondrial function under time-varying magnetic stimulation, and provided evidence that TMS does not affect normal mitochondrial functionality by altering its membrane potential.
Time-varying magnetic fields have been used to stimulate neural tissues since the start of 20th century . Today, pulsed magnetic fields are used in stimulating the central nervous system, via a technique named transcranial magnetic stimulation (TMS). TMS is being explored in the treatment of depression , seizures [3, 4], Parkinson's disease , and Alzheimer's disease [6, 7]. It also facilitates long-lasting plastic changes induced by motor practice, leading to more remarkable and outlasting clinical gains during recovery from stroke or traumatic brain injury .
When exposed to a time-varying magnetic field, the neural tissue is stimulated by an electric current via electromagnetic induction , which induces an electrical potential that is superimposed on the resting membrane potential of the cell. The polarization could be controlled by appropriate geometrical positioning of the magnetic coil [10–12]. To investigate the effects of stimulation, theoretical studies have been performed to compute the magnetically induced electric field and potentials in the neuronal tissue, using models that represent nerve fibers [13–18] or cell bodies .
Mitochondria are involved in a large range of physiological processes such as supplying cellular energy, signaling, cellular differentiation, cell death, as well as the control of cell cycle and growth . Their large negative membrane potential (-180 mV) in the mitochondrial inner membrane, which is generated by the electron-transport chain, is the main driving force in these regulatory processes [21–23]. Alteration of this large negative membrane potential has been associated with disruption in cellular homeostasis that leads to cell death in aging and many neurological disorders [24–27]. Thus, mitochondria can be a therapeutic target in many neurodegenerative diseases such as Alzheimer's disease and Parkinson's disease.
Two lines of evidences suggest that the physiology of mitochondria could be affected by the magnetic field via its induced transmembrane potential. First, magnetic fields can induce electric fields in the neural tissue, and it has been shown that exposure of a cell to an electrical field could introduce a voltage on the mitochondrial membrane . This induced potential has led to many physiological/pathological changes, such as opening of the mitochondrial permeability transition pore complex . Nanosecond pulsed electric fields (nsPEFs) can affect mitochondrial membrane [30, 31], cause calcium release from internal stores , and induce mitochondria-dependent apoptosis under severe stress [33, 34]. Secondly, there is evidence that magnetic fields could alter several important physiological processes that are related to the mitochondrial membrane potential, including ATP synthesis [35, 36], metabolic activities [37, 38] and Ca2+ handling [39, 40]. An analysis of the mitochondrial membrane potential is of experimental significance in understanding its physiology/pathology under magnetic stimulation.
In this theoretical work, we have provided the first analytical solution for the transmembrane potential in an internal organelle (i.e., mitochondrion) that is induced by a time-varying magnetic field. The model was a two-shell cell structure, with the outer shell representing the cell membrane and the inner shell representing the membrane of an internal organelle. Factors that affect the amount of organelle polarization were investigated by parametric analysis, including field frequency, and properties of the cytoplasmic and organelle membranes. We also estimated to what degree magnetic fields used in TMS practice affect organelle polarization.
Spherical cell and internal organelle model in a time-varying magnetic field
Figure 1 shows the model representation of the cell membrane and the internal organelle, and their orientation to the coil that generates the magnetic field. Two coordinate systems were utilized to represent the cell and the coil, respectively.
The co-centric spherical cell and the organelle were represented in a spherical coordinate system (r, θ, ϕ) centered at point O. The cell membrane was represented as a very thin shell with inner radius R-, outer radius R+ and thickness D. The organelle membrane was represented as a very thin shell with inner radius r-, outer radius r+ and thickness d. The two membrane shells divided the cellular environment into five homogenous, isotropic regions: extracellular medium (0#), cytoplasm membrane (1#), intracellular cytoplasm (2#), organelle membrane (3#) and the organelle internal (4#). The dielectric permittivities and the conductivities in the five regions were ε i and σ i , respectively, where i represents the region number.
The low-frequency magnetic field was represented in a cylindrical coordinate system (r', ϕ', z'). The distance between the center of the cell (O) and the axis of the coil (O') was C. The externally applied, sinusoidally alternating magnetic field was symmetric about the O' Z' axis. The magnetic field was represented as , where was the unit vector in the direction of O' Z', ω was the angular frequency of the magnetic field, and was the imaginary unit.
Table 1 lists the parameters used for the model. To quantitatively investigate the amount of polarization on both the cytoplasmic and organelle membranes, we chose their geometrical and electrical parameters (standard values, the lower and upper limits) from the literature . The frequency range of interest was determined to be between 2 - 200 kHz. The upper limit was determined by calculating the reciprocal value of the rising phase of a current pulse during peripheral nerve stimulation [42, 43]. Most frequencies used in the experimental practices were lower than this value . The intensity of the magnetic field was 2 Tesla from TMS practice. The standard frequency of the magnetic field was estimated to be 10 kHz, as the rising time of single pulses was ~100 μs during TMS. This yielded the peak value of dB/dt = 2 × 104T/s .
Governing equations for potentials and electric fields induced by the time-varying magnetic field
The electric field induced by the time varying magnetic field in the biological media was
where is the magnetic vector potential induced by the current in the coil. The potential V was the electric scalar potential due to charge accumulation that appears from the application of a time-varying magnetic field . In spherical coordinates (r, θ, ϕ), . Using quasi-static approximations, in charge-free regions, V was obtained by solving Laplace's equation
Four boundary conditions were considered in the derivation of the potentials induced by the time-varying magnetic field.
(A). The potential was continuous across the boundary of two different media. In this paper, this refers to the extracellular media/membrane interface (0#1#), the cell membrane/intracellular cytoplasm interface (1#2#), the intracellular cytoplasm/organelle membrane interface (2#3#), and the organelle membrane/organelle interior interface (3#4#).
(B). The normal component of the current density was continuous across two different media. For materials such as pure conductors, it was equal to the product of the electric field and the conductivity of the media. During time-varying field stimulation, the complex conductivity, defined as S = σ +jωε, was used to account for the dielectric permittivity of the material . Here, σ was the conductivity, ε was the dielectric permittivity of the tissue, ω was the angular frequency of the field. Therefore, on the extracellular media/membrane interface (0#1#),
On the cell membrane/intracellular cytoplasm interface (1#2#),
On the intracellular cytoplasm/organelle membrane interface (2#3#),
On the organelle membrane/organelle interior interface (3#4#),
where S0 = σ0+jωε0, S1 = σ1+jωε1, S2 = σ2+jωε2, S3 = σ3+jωε3 and S4 = σ4+jωε4 were the complex conductivities of the five media, respectively.
(C). The electric field at an infinite distance from the cell was not perturbed by the presence of the cell.
(D). The potential inside the organelle (r = 0) was finite.
Magnetic vector potential
When the center of the magnetic field was at point O', was in the direction of since
In order to calculate the potential distribution in the model cell, one needs to have an expression for in spherical coordinates(r, θ, ϕ). By coordinate transformation (Appendix B in ), we obtained the magnetic vector potential in spherical coordinates (r, θ, ϕ):
The vector potential components in the , , directions were:
Derivations of the equations were done with Mathematica 6.0 (Wolfram Research, Inc. Champaign, IL). Numerical simulations were done with Matlab 7.4.0 (The MathWorks, Inc. Natick, MA).
Transmembrane potentials induced by a time-varying magnetic field
In spherical coordinates (r, θ, ϕ), the solution for Laplace's equation (2) can be written in the form
where C n , D n were unknown coefficients (n = 0,1,2,3,4,5). We solved for those coefficients (Appendix) and substituted them into equation (13) to obtain the potential terms in the five model regions. Next, the transmembrane potential in a membrane can be obtained by subtracting the membrane potential at the inner surface from that at the outer surface.
In the cell membrane, the induced transmembrane potential was
In the organelle membrane, the induced transmembrane potential was
Similar regional polarization patterns were observed between the cell membrane and the organelle membrane, since they both depended on a sinθ cosϕ term. Since θ and ϕ were determined by the relative orientation of the coil to the cell, the patterns of polarization in the target cell and the organelle both depended on their orientations to the stimulation coil.
ψ cell and ψ org at one instant moment were plotted for 10 KHz and 100 KHz, respectively (Figure 2). The locations for the maximal polarization were on the equators of the cell and of the organelle membranes (θ = 90° or z = 0 plane). The two membranes were maximally depolarized at ϕ = 180° (deep red) and maximally hyperpolarized at ϕ = 0 (deep blue) on the equator, respectively. The cell and the organelle membranes were not polarized on the two poles corresponding to θ = 0° and θ = 180°, respectively. The full cycle of polarization by the time-varying magnetic field was also illustrated (see Additional file 1).
Both ψ cell and ψ org depended on the geometrical parameters of the cell (R+, R-, C) and the organelle (r+, r-), and the electrical properties of the five media considered in the model (S0, S1, S2, S3, S4). These parameters did not affect the polarization pattern. Therefore, we chose maximal polarizations (corresponding to the point that is defined by θ = 90°, ϕ = 270°) on the cell and organelle membranes (Figures 1 and 2) for the further analysis of their dependency on the field frequency.
Two factors contribute to the frequency-dependency of the polarizations (magnitude and phase) in the two membranes. First, the magnitude of the electrical field is proportional to the frequency of the externally applied magnetic field, as required by Faraday's law. Second, the dielectric properties of the material considered in the model are frequency-dependent.
With the standard values, ψ cell was always greater than and ψ org (Figure 3A). At 10 kHz, the maximal polarization on the cell membrane was 9.397 mV, and the maximum polarization on the internal organelle was only 0.08 mV. Figure 3B plots the ratio of the two polarizations. As the frequency increased, ψ org became quantitatively comparable to ψ cell . At extremely high frequency (~100 MHz), the ratio reached a plateau of 1 (not shown).
The phase was defined as the phase difference between the externally applied magnetic field and membrane polarization, which was computed as the phase angle of the complex transmembrane potentials. Phase in the cell membrane was insensitive to the frequency change below 10 KHz. At 10 KHz, the phase in the cell membrane is -91.23°, which meant that an extra -1.23° was added to the membrane phase, due to frequency-dependent capacitive features of the tissue. On the other hand, phase response in the organelle membrane was more sensitive to the frequency change than the cell membrane, showing the dependence as low as 50 Hz. At 10 KHz, the phase in the organelle was -5.69°. Above 10 KHz, phases in both membranes increased with frequency. At 200 KHz, the phase in the cell membrane was -113.1°, and in the organelle membrane was -33.07°. Figure 3D plots the difference between the two phases as a function of frequency. At very low frequency (< 50 Hz), the two membranes demonstrated an in-phase polarization. At 10 KHz, their polarizations were nearly 90° out-of-phase.
"Interaction" between the cell membrane and the organelle membrane
Previous studies have shown that the cell membrane "shields" the internal cytoplasm and prevent the external field from penetrating inside the cell in electric stimulation [48, 49]. Will similar phenomenon occur under magnetic stimulation? To estimate the impact of cell membrane on organelle polarization, we compared ψ org with and without the presence of the cell membrane. The later was done by letting S1 = S0 and S2 = S0 in equation (15), which removed the cell membrane,
Removal of the cell membrane allowed greater organelle polarization (Figure 4A). At 10 KHz, ψ org was 2.82 mV in the absence of the cell membrane, which was considerably greater than 0.08 mV for the case with the cell membrane. This screening effect was more prominent at 200 KHz, where ψ org was only 28.78 mV in the intact cell, and 55.87 mV without the cell membrane.
The phase response for the isolated organelle was similar to a cell membrane that was directly exposed in the field (Figure 4B). Therefore, presence of the cell membrane not only" shielded" the internal mitochondria from excessive polarization by the external field, but also provides an extra phase term that reduce the phase delay between the field and the organelle response.
Alteration in the organelle polarization by removing the cell membrane suggested an "interactive" effect between the two membranes via electric fields. We next asked if the presence of the internal organelle might have the reciprocal effects on ψ cell . To test this possibility, we removed the internal organelle and investigated its effect on ψ cell . This was done by letting S3 = S2 and S4 = S2 in equation (14). Removal of the internal organelle did not introduce significant changes on ψ cell (Figure 5). Removal of the organelle led to a 0.001 mV increase in ψ cell at 10 KHz, and a 1.3 mV increase at 200 KHz, respectively. The phase change caused by organelle removal was only 0.7 degrees at 200 KHz. These results suggest that the presence of the internal organelle only had trivial effects on the cytoplasmic membrane.
Dependency of ψ org on the cell membrane parameters
To further investigate the shielding effects of the cell membrane on ψ org , we systemically varied the cell membrane parameters within their physiological ranges, and studied their individual impacts on the organelle polarization. These parameters included the geometrical properties (radius and membrane thickness) and the electrical properties (cell membrane conductivity and dielectric permittivity) of the cell membrane. This was done by varying one parameter through its given range but maintaining the others at their standard values. Since the dielectric properties of the tissues were frequency dependent, the parameter sweep was done within a frequency range (2 - 200 KHz). This generated a set of data that could be depicted in a color plot of ψ org (amplitude or phase) as a function of frequency and the studied parameters (Figures 6).
At a low frequency band (< 10 KHz), ψ org was trivial, since the magnitude of the induced electric field was small. ψ org became considerably large beyond 10 KHz. Increase in the cell radius facilitates this polarization (Figure 6A left). Increase in the cell radius did not significantly change the phase-frequency relation in the organelle. However, it increased the phase at relatively high frequency (~100 KHz, Figure 6A right). Increase in the cell membrane thickness compromised ψ org , so that higher frequency was needed to induce considerable polarization in the organelle (Figure 6B left). Variation in membrane thickness did not significantly alter the phase of the organelle polarization (Figure 6B right). Since removal of the low-conductive cell membrane enhanced organelle polarization (Figure 4A), one might expect that an increase in the membrane conductivity could have a similar effect. However, within the physiological range considered in this paper, ψ org was insensitive to the cell membrane conductivity (Figure 6C left). The cell membrane conductivity did have a significant impact on the phase of mitochondria polarization. At extremely low values (<10-7S/m), ψ org demonstrated a phase advance at frequency lower than 1 KHz (Figure 6C right), rather than a phase delay, as was the case for the standard values (Figure 3C). The cell membrane dielectric permittivity represents the capacitive property of the membrane. Increase in this parameter facilitated ψ org , so that ψ org became noticeable at relatively lower frequency range (Figure 6D left). An increase in this parameter also led to a decrease in the phase delay in the organelle polarization, which was most prominent at the frequency above 100 Hz (Figure 6D right).
Dependency of ψ org on its own biophysics
Previous studies have shown that polarization of a neuronal structure depends on its own membrane properties under both electrical , and magnetic stimulations . How do the membrane properties of the organelle membrane affect its own polarization?
An increase in the organelle radius led to a greater ψ org (Figure 7A, left). The phase-frequency relationship differentiated at a radius value around 1.1 um. Above this value, the phase response followed a pattern depicted in Figure 3C, i.e., the phase delay was -90 degree for low frequency and decreased to 0 at around 10 K Hz. Below this value, the phase showed a 90-degree advance instead of a lag in the low frequency range < 10 K Hz (Figure 7A, right). The membrane thickness has been generally agreed to be least significant to membrane polarization . Varying membrane thickness in the organelle did not cause significant change in the magnitude (Figure 7B, left) nor the phase (Figure 7B, right) of ψ org . ψ org was also insensitive to its own electrical properties. Varying membrane conductivity (Figure 7C) or dielectricity (Figure 7D) in the organelle did not alter the frequency-dependent polarization in this structure.
Similarities and differences to electrical stimulation
Analysis of ψ org under magnetic stimulation reveals several commonalities and differences to that under electric stimulation. The build up of ψ org requires the electric field to penetrate through the cytoplasmic membrane. In electric stimulation, this is achieved by directly applied electric current via electrodes. In magnetic stimulation, electric field is produced by electromagnetic induction.
Analysis on ψ org under electric field has been performed in two recent publications. Vajrala et al.  developed a three-membrane model that included the inner and our membranes of a mitochondrion, and have analytically solved ψ cell and ψ org under oscillatory electric fields. Another study  has modeled the internal membrane response to the time-varying electric field, and has investigated the condition under which ψ org can temporarily exceed ψ cell under nanosecond duration pulsed electric fields.
Results obtained from this magnetic study share several commonalities with those from AC electric stimulation. Under both stimulation conditions, ψ org can never exceed ψ cell . The ratio between the (organelle/cell) increases with frequency, and this ratio can reach 1 at very high frequency (108 Hz, data not shown). The phase responses of the organelle within a cell have not been analyzed previously under electric stimulation, which prevent direct comparison with this work. For an isolated mitochondrion, its response is similar to a single cell membrane under AC electric field stimulation , except that an extra -90° phase is introduced by electromagnetic induction (Figure 4B).
Stimulation on the internal organelle by time-varying magnetic field, though, has its own uniqueness. First, as a non-invasive method, magnetic stimulation is achieved by current induction inside the tissue, which prevents direct contact with the electrodes and introduces minimal discomfort. Second, the frequency responses of the internal organelle are different under the two stimulation protocols. In electric stimulation, magnitude of the field is independent of its frequency. In magnetic stimulation, however, the magnitude of the induced electric field is proportional to the frequency of the magnetic field (Faraday's law). Consequently, alteration in the field frequency could also contribute to ψ org . Low frequency field (< 1 KHz) is insufficient in building up noticeable ψ org and ψ cell (Figure 3A). Both ψ org and ψ cell increase with field frequency (Figure 3A). Therefore, it is unlikely possible to use high-frequency magnetic field to specifically target internal organelles, such as been done under AC electric stimulation with nanosecond pulses, for mitochondria electroporation and for the induction of mitochondria-dependent apoptosis .
Cellular factors that influence ψ cell
When a neuron is exposed to an electric field, a transmembrane potential is induced on its membrane. Attempts to analytically solve ψ cell began as early as the 1950s [51, 52]. Later works added more complexity to the modeled cell and provided insights into the factors affecting ψ cell . These include electrical properties [49, 50, 53, 54] of the cell, such as its membrane conductivity. Geometrical properties of the cell could also affect ψ cell , such as its shape [55, 56] and orientation to the field [57, 58].
Presence of neighboring cells affect ψ cell in a tissue with high-density cells, For example, isthmo-optic cells in pigeons can be excited by electrical field effect through ephaptic interaction produced by the nearby cells whose axons were activated by electric stimulation, suggesting that electrical field effect may play important roles in interneuronal communications . In infinite cell suspensions, ψ cell depended on cell volume fraction and cell arrangement . Theoretical studies have proved that presence of a single cell affected ψ cell in its neighboring cells, without direct physical contact between the two cells .
This work investigates another important factor that might affect ψ cell , i.e., presence of the internal organelle. We have previously solved ψ cell for a spherical cell model under magnetic field stimulation, without considering the presence of the internal organelle . This work extends the previous study by including an internal organelle in the cell model. Here, adding an organelle to the cell internal did not significantly change the magnitude and phase of ψ cell (Figure 5).
Factors that influence ψ org during magnetic stimulation
Biological tissue is composed of many non-homogenous, anisotropic components, such as the cellular/axonal membrane, the internal organelles and the extracellular medium. The electrical properties (i.e., conductivities) of the tissue may vary with location in the tissue, even at a microscopic level. Under magnetic stimulation, several studies have provided insights into the impact of tissue properties on field distribution and tissue polarization [42, 61].
This work further illustrates that the effects of magnetic stimulation are a function of tissue properties, by providing evidence that both the geometrical and electrical parameters of the cell/organelle membranes affect ψ org . Both the radius of the cell and the organelle strongly affect ψ org , which is in agree with previous studies [48, 62]. Radius of the neuronal structure is important in determining the threshold for its own membrane polarization, as proved by in vitro studies on eukaryotic  and bacterial cells . This model prediction is potentially testable with voltage-sensitive dyes that can provide both high temporal and high spatial resolutions [23, 65]. Another model prediction is that the amount of ψ org is insensitive to the change in cell membrane conductivity. Evidence has shown that electric field can cause long-lasting increase in passive electrical conductance of the cell membrane, probably by opening of stable conductance pores . The opening and closing of ion channels can also alter the membrane conductance. This model prediction can be tested by varying membrane conductivity, using ion-channel blockers applied to the cell membrane.
Implications for transcranial magnetic stimulation (TMS)
Another important finding in this study that within the frequency band used TMS, ψ org is insignificant comparing with ψ cell . At 10 KHz, a frequency that corresponds to the rising time of the electric pulses used in clinical TMS, the field causes considerable amount of change in ψ cell , but only 0.08 mV change in ψ org (Figure 3A). It is worth noting that even this value was probably a consequence of overestimation in the magnetic field intensity (B0). To simplify the calculation, B0 was a constant (2 Tesla) everywhere in the modeled region. In reality, the intensity of the magnetic field generated by a coil could decay quickly in the tissue far away from the coil [67, 68]. The duration of the stimulation time was also likely overestimated. During TMS, neuronal responses are induced by pulses, as opposed to the mathematically more tractable sinusoidal stimulus used in this model. Under this scenario, the magnetically-induced electric field in the tissue (essentially the change in the transmembrane potential) is determined by , which means the transmembrane potential can only be induced during the rise time (and decay time) during a step in the B field. Indeed, rise times of the field affect stimulation in clinic practice, and a faster rise time pulse is more efficient . Therefore, ψ org is unlikely significant enough in TMS to have physiological implications, and internal organelles such as mitochondria are not likely be the target in TMS practice. This conclusion is made after extensive analysis on model parameters with the values in broad physiological ranges (Table 1). To our knowledge and based on a Medline search, there have been no reports on mitochondria-related effects in TMS practice.
This paper provides two mechanisms to account for the ineffectiveness of magnetically-induced polarization in internal organelles under TMS parameters. First, the cell membrane, which is made up of lipids and proteins, provides a dominant "shielding effect" on the organelles and prevents certain amount of electric fields to penetrate into the cell membrane and polarize the organelle membrane (Figure 4). Second, the radius of the organelle is always much smaller than that of the cell, which render them relatively insensitive to the magnetic field.
Several simplifying assumptions were proposed in this model to facilitate the derivation of the analytical solutions. The model assumed that the cell was located in an electrically homogenous extracellular medium, which was an over-simplification of the true electrically anisotropic extracellular environment. Both the extracellular medium and cytoplasmic environment are not truly homogenous [69, 70]. We found that neither parameter significantly affects the organelle or cytoplasmic membrane polarization (not shown).
Both the cell membrane and the mitochondria membranes were modeled as a single spherical shell. In reality, however, cellular structures have irregular shape, which may play an important role in the dynamics of membrane polarization [71, 72]. The interior sphere was centered inside the cell to allow for mathematical simplicity of the model. However, as organelle locations vary spatially in a cell, we hypothesize that organelles located off-center of the cell or closer to the exterior cell membrane may be more sensitive to the applied field. Also, we believe the "shielding effect" of the cell membrane persists even when the separation distance between the two membranes is small (data not shown). The membrane of the organelle was modeled as a single internal shell as in a previous study , rather than a two-shell structure, representative of the inner and our membranes of a mitochondrion, respectively . The highly curved projections of the cell body and the organelle membrane may provide focal points for even greater changes in the induced transmembrane potential . Future study will use numerical methods with multi-compartment modeling or finite element meshes to represent these structure complexities.
All the dielectric permittivities in the model were assumed to be frequency-independent, which was valid for the low frequencies considered (10-200 kHz). When field frequency exceeds several hundreds of megahertz, the finite mobility of molecular dipoles starts to weaken the polarization processes . This phenomenon, known as dielectric relaxation, is characterized with decrease in the permittivities and increase in the conductivity. When this happens, the complex conductivity should be defined as S = σ (ω) + jωε (ω), where σ (ω) and ε (ω) are frequency-dependent conductivity and permittivity, respectively. By implementing this term in equations (14) and (15), one can estimate the transmembrane potentials in the cell and in the organelle when dielectric relaxation occurs.
This work provides the first analytical solution for the transmembrane potentials in an internal organelle (ψ org ) in response to time-varying magnetic stimulation. The frequency response of the membrane under magnetic stimulation is different from that under electric field stimulation. This work provides evidence that the presence of the internal organelle does not significantly affect polarization of the cell membrane (ψ cell ). Moreover, ψ org is always smaller than ψ cell under low frequency range (< 200 KHz), largely due to the "shielding effect" imposed by the presence of the cell membrane. Both the geometrical and electrical properties of the cell membrane affect ψ org in a frequency-dependent manner. The properties of the organelle membrane also affect ψ org in a frequency-dependent manner. Finally, the present study provides evidence that normal mitochondrial functionality is not likely affected by transcranial magnetic stimulation, via altering its membrane potential.
Determining unknown coefficients C n , D n in equation (13) using boundary conditions
Since V was bounded at r = 0 and r → ∞, from equation (13) we had
Therefore, expressions for the potential distribution in the extracellular media, the cell membrane, the cytoplasm, the organelle membrane, and organelle interior are:
We substituted A0r(equation 10) and the components of ∇V in the five regions into (1) to yield the expressions of the normal components of the electric fields in the five regions:
Following boundary condition (A), V was continuous at the extracellular media/membrane (r = R+), the membrane/intracellular cytoplasm interfaces (r = R-), the cytoplasm/organelle interface and the organelle membrane/organelle interior interface.
We then used the boundary condition (B), that the normal components of the current densities were continuous between two different media (equations 3-6), to obtain the following equations:
We solved (A-11) to (A-18) the last eight unknown coefficients D0-D3, C1-C4. (see Additional file 2).
Thompson SP: A Physiological Effect of an Alternating Magnetic Field. Proc R Soc 1910, 82: 396-398. 10.1098/rspb.1910.0032
Epstein CM, Davey KR: Iron-core coils for transcranial magnetic stimulation. J Clin Neurophysiol 2002, 19: 376-381. 10.1097/00004691-200208000-00010
Anninos PA, Tsagas N, Sandyk R, Derpapas K: Magnetic stimulation in the treatment of partial seizures. Int J Neurosci 1991, 60: 141-171.
Anninos PA, Tsagas N, Jacobson JI, Kotini A: The biological effects of magnetic stimulation in epileptic patients. Panminerva Med 1999, 41: 207-215.
Sandyk R, Anninos PA, Tsagas N, Derpapas K: Magnetic fields in the treatment of Parkinson's disease. Int J Neurosci 1992, 63: 141-150. 10.3109/00207459208986664
Cotelli M, Manenti R, Cappa SF, Zanetti O, Miniussi C: Transcranial magnetic stimulation improves naming in Alzheimer disease patients at different stages of cognitive decline. Eur J Neurol 2008, 15: 1286-1292. 10.1111/j.1468-1331.2008.02202.x
Julkunen P, Jauhiainen AM, Westeren-Punnonen S, Pirinen E, Soininen H, Kononen M, Paakkonen A, Maatta S, Karhu J: Navigated TMS combined with EEG in mild cognitive impairment and Alzheimer's disease: a pilot study. J Neurosci Methods 2008, 172: 270-276. 10.1016/j.jneumeth.2008.04.021
Bolognini N, Pascual-Leone A, Fregni F: Using non-invasive brain stimulation to augment motor training-induced plasticity. J Neuroeng Rehabil 2009, 6: 8. 10.1186/1743-0003-6-8
Kobayashi M, Pascual-Leone A: Transcranial magnetic stimulation in neurology. Lancet Neurol 2003, 2: 145-156. 10.1016/S1474-4422(03)00321-1
Cohen D: Feasibility of a Magnetic Stimulator For the Brain. In Biomagnetism: Application and Theory. Edited by: Weinberg H, Stroink G, Katila T. Pergamon Press; 1984:466-470.
Nagarajan SS, Durand DM, Hsuing-Hsu K: Mapping location of excitation during magnetic stimulation: effects of coil position. Ann Biomed Eng 1997, 25: 112-125. 10.1007/BF02738543
Roth Y, Zangen A, Hallett M: A coil design for transcranial magnetic stimulation of deep brain regions. J Clin Neurophysiol 2002, 19: 361-370. 10.1097/00004691-200208000-00008
Roth BJ, Basser PJ: A model of the stimulation of a nerve fiber by electromagnetic induction. IEEE Trans Biomed Eng 1990, 37: 588-597. 10.1109/10.55662
Ravazzani P, Ruohonen J, Grandori F, Tognola G: Magnetic stimulation of the nervous system: induced electric field in unbounded, semi-infinite, spherical, and cylindrical media. Ann Biomed Eng 1996, 24: 606-616. 10.1007/BF02684229
Esselle KP, Stuchly MA: Quasi-static electric field in a cylindrical volume conductor induced by external coils. IEEE Trans Biomed Eng 1994, 41: 151-158. 10.1109/10.284926
Esselle KP, Stuchly MA: Cylindrical tissue model for magnetic field stimulation of neurons: effects of coil geometry. IEEE Trans Biomed Eng 1995, 42: 934-941. 10.1109/10.412660
Roth BJ, Cohen LG, Hallett M, Friauf W, Basser PJ: A theoretical calculation of the electric field induced by magnetic stimulation of a peripheral nerve. Muscle Nerve 1990, 13: 734-741. 10.1002/mus.880130812
Nagarajan SS, Durand DM: Analysis of magnetic stimulation of a concentric axon in a nerve bundle. IEEE Trans Biomed Eng 1995, 42: 926-933. 10.1109/10.412659
Ye H, Cotic M, Carlen PL: Transmembrane potential induced in a spherical cell model under low-frequency magnetic stimulation. J Neural Eng 2007, 4: 283-293. 10.1088/1741-2560/4/3/014
McBride HM, Neuspiel M, Wasiak S: Mitochondria: more than just a powerhouse. Curr Biol 2006, 16: R551-560. 10.1016/j.cub.2006.06.054
Berg J, Tymoczko J, Stryer L: Biochemistry. 5th edition. WH Freeman and Company; 2002.
Gunter KK, Gunter TE: Transport of calcium by mitochondria. J Bioenerg Biomembr 1994, 26: 471-485. 10.1007/BF00762732
Nicholls DG, Budd SL: Mitochondria and neuronal survival. Physiol Rev 2000, 80: 315-360.
Rizzuto R, Bastianutto C, Brini M, Murgia M, Pozzan T: Mitochondrial Ca2+ homeostasis in intact cells. J Cell Biol 1994, 126: 1183-1194. 10.1083/jcb.126.5.1183
Beal MF: Mitochondrial dysfunction in neurodegenerative diseases. Biochim Biophys Acta 1998, 1366: 211-223. 10.1016/S0005-2728(98)00114-5
Miller RJ: Mitochondria - the Kraken wakes! Trends Neurosci 1998, 21: 95-97. 10.1016/S0166-2236(97)01206-X
Schapira AH, Gu M, Taanman JW, Tabrizi SJ, Seaton T, Cleeter M, Cooper JM: Mitochondria in the etiology and pathogenesis of Parkinson's disease. Ann Neurol 1998, 44: S89-98.
Vajrala V, Claycomb JR, Sanabria H, Miller JH Jr: Effects of oscillatory electric fields on internal membranes: an analytical model. Biophys J 2008, 94: 2043-2052. 10.1529/biophysj.107.114611
Weaver J: Electroporation of Biological Membranes From Multicellular to Nano Scales. IEEE Trans Dielectr Electr Insul 2003, 10: 754-768. 10.1109/TDEI.2003.1237325
White JA, Blackmore PF, Schoenbach KH, Beebe SJ: Stimulation of capacitative calcium entry in HL-60 cells by nanosecond pulsed electric fields. J Biol Chem 2004, 279: 22964-22972. 10.1074/jbc.M311135200
Tekle E, Oubrahim H, Dzekunov SM, Kolb JF, Schoenbach KH, Chock PB: Selective field effects on intracellular vacuoles and vesicle membranes with nanosecond electric pulses. Biophys J 2005, 89: 274-284. 10.1529/biophysj.104.054494
Scarlett SS, White JA, Blackmore PF, Schoenbach KH, Kolb JF: Regulation of intracellular calcium concentration by nanosecond pulsed electric fields. Biochim Biophys Acta 2009, 1788: 1168-1175. 10.1016/j.bbamem.2009.02.006
Beebe SJ, Fox PM, Rec LJ, Willis EL, Schoenbach KH: Nanosecond, high-intensity pulsed electric fields induce apoptosis in human cells. FASEB J 2003, 17: 1493-1495.
Schoenbach KH, Beebe SJ, Buescher ES: Intracellular effect of ultrashort electrical pulses. Bioelectromagnetics 2001, 22: 440-448. 10.1002/bem.71
Feng HL, Yan L, Cui LY: Effects of repetitive transcranial magnetic stimulation on adenosine triphosphate content and microtubule associated protein-2 expression after cerebral ischemia-reperfusion injury in rat brain. Chin Med J (Engl) 2008, 121: 1307-1312.
Dimroth P, Kaim G, Matthey U: Crucial role of the membrane potential for ATP synthesis by F(1)F(o) ATP synthases. J Exp Biol 2000, 203: 51-59.
Yamashita K, Saito M: Effects of Middle-Level Static Magnetic Field on Metabolic Activity of Mitochondria. Electr Eng Jpn 2001, 137: 36-41. 10.1002/eej.1078
Valdez LB, Zaobornyj T, Boveris A: Mitochondrial metabolic states and membrane potential modulate mtNOS activity. Biochim Biophys Acta 2006, 1757: 166-172. 10.1016/j.bbabio.2006.02.013
Belyavskaya NA: Biological effects due to weak magnetic field on plants. Adv Space Res 2004, 34: 1566-1574. 10.1016/j.asr.2004.01.021
Kovacs R, Kardos J, Heinemann U, Kann O: Mitochondrial calcium ion and membrane potential transients follow the pattern of epileptiform discharges in hippocampal slice cultures. J Neurosci 2005, 25: 4260-4269. 10.1523/JNEUROSCI.4000-04.2005
Kotnik T, Miklavcic D: Theoretical evaluation of voltage inducement on internal membranes of biological cells exposed to electric fields. Biophys J 2006, 90: 480-491. 10.1529/biophysj.105.070771
Krasteva VT, Papazov SP, Daskalov IK: Peripheral nerve magnetic stimulation: influence of tissue non-homogeneity. Biomed Eng Online 2003, 2: 19. 10.1186/1475-925X-2-19
Ruohonen J, Panizza M, Nilsson J, Ravazzani P, Grandori F, Tognola G: Transverse-field activation mechanism in magnetic stimulation of peripheral nerves. Electroencephalogr Clin Neurophysiol 1996, 101: 167-174. 10.1016/0924-980X(95)00237-F
Mansfield P, Harvey PR: Limits to neural stimulation in echo-planar imaging. Magn Reson Med 1993, 29: 746-758. 10.1002/mrm.1910290606
Barker A, Freeston IL, Garnham CW: Measurement of Cortical and Peripheral Neural Membrane Time Constant in Man Using Nerve Stimulation. J Physiol (Lond) 1990, 423: 66.
Stratton J: Electromagnetic Theory. New York: McGraw-Hill; 1941.
Kotnik T, Miklavcic D: Second-order model of membrane electric field induced by alternating external electric fields. IEEE Trans Biomed Eng 2000, 47: 1074-1081. 10.1109/10.855935
Mossop BJ, Barr RC, Zaharoff DA, Yuan F: Electric fields within cells as a function of membrane resistivity--a model study. IEEE Trans Nanobioscience 2004, 3: 225-231. 10.1109/TNB.2004.833703
Lee DC, Grill WM: Polarization of a spherical cell in a nonuniform extracellular electric field. Ann Biomed Eng 2005, 33: 603-615. 10.1007/s10439-005-2397-3
Kotnik T, Bobanovic F, Miklavcic D: Sensitivity of Transmembrane Voltage Induced by Applied Electric Fields--A Theoretical Analysis. Bioelectrochem Bioenerg 1997, 43: 285-291. 10.1016/S0302-4598(97)00023-8
Fricke H: The Electric Permittivity of a Dilute Suspension of Membrane-Covered Ellipsoids. J Appl Phys 1953, 24: 644-646. 10.1063/1.1721343
Schwan HP: Electrical properties of tissue and cell suspensions. Adv Biol Med Phys 1957, 5: 147-209.
DeBruin KA, Krassowska W: Modeling electroporation in a single cell. II. Effects Of ionic concentrations. Biophys J 1999, 77: 1225-1233. 10.1016/S0006-3495(99)76974-2
DeBruin KA, Krassowska W: Modeling electroporation in a single cell. I. Effects Of field strength and rest potential. Biophys J 1999, 77: 1213-1224. 10.1016/S0006-3495(99)76973-0
Gimsa J, Wachner D: Analytical description of the transmembrane voltage induced on arbitrarily oriented ellipsoidal and cylindrical cells. Biophys J 2001, 81: 1888-1896. 10.1016/S0006-3495(01)75840-7
Kotnik T, Miklavcic D: Analytical description of transmembrane voltage induced by electric fields on spheroidal cells. Biophys J 2000, 79: 670-679. 10.1016/S0006-3495(00)76325-9
Pavlin M, Pavselj N, Miklavcic D: Dependence of induced transmembrane potential on cell density, arrangement, and cell position inside a cell system. IEEE Trans Biomed Eng 2002, 49: 605-612. 10.1109/TBME.2002.1001975
Valic B, Golzio M, Pavlin M, Schatz A, Faurie C, Gabriel B, Teissie J, Rols MP, Miklavcic D: Effect of electric field induced transmembrane potential on spheroidal cells: theory and experiment. Eur Biophys J 2003, 32: 519-528. 10.1007/s00249-003-0296-9
Hu J, Li WC, Xiao Q, Wang SR: Electrical interaction between neurons in the pigeon isthmo-optic nucleus. Brain Res Bull 2000, 51: 159-163. 10.1016/S0361-9230(99)00211-7
Pucihar G, Kotnik T, Valic B, Miklavcic D: Numerical determination of transmembrane voltage induced on irregularly shaped cells. Ann Biomed Eng 2006, 34: 642-652. 10.1007/s10439-005-9076-2
Miranda PC, Hallett M, Basser PJ: The electric field induced in the brain by magnetic stimulation: a 3-D finite-element analysis of the effect of tissue heterogeneity and anisotropy. IEEE Trans Biomed Eng 2003, 50: 1074-1085. 10.1109/TBME.2003.816079
Farkas DL, Korenstein R, Malkin S: Electrophotoluminescence and the electrical properties of the photosynthetic membrane. I. Initial kinetics and the charging capacitance of the membrane. Biophys J 1984, 45: 363-373. 10.1016/S0006-3495(84)84160-0
Rols MP, Delteil C, Serin G, Teissie J: Temperature effects on electrotransfection of mammalian cells. Nucleic Acids Res 1994, 22: 540. 10.1093/nar/22.3.540
Calvin NM, Hanawalt PC: High-efficiency transformation of bacterial cells by electroporation. J Bacteriol 1988, 170: 2796-2801.
Perez Velazquez JL, Frantseva MV, Huzar DV, Carlen PL: Mitochondrial porin required for ischemia-induced mitochondrial dysfunction and neuronal damage. Neuroscience 2000, 97: 363-369. 10.1016/S0306-4522(99)00569-2
Pakhomov AG, Shevin R, White JA, Kolb JF, Pakhomova ON, Joshi RP, Schoenbach KH: Membrane permeabilization and cell damage by ultrashort electric field shocks. Arch Biochem Biophys 2007, 465: 109-118. 10.1016/j.abb.2007.05.003
Tofts PS, Branston NM: The measurement of electric field, and the influence of surface charge, in magnetic stimulation. Electroencephalogr Clin Neurophysiol 1991, 81: 238-239. 10.1016/0168-5597(91)90077-B
Eaton H: Electric field induced in a spherical volume conductor from arbitrary coils: application to magnetic stimulation and MEG. Med Biol Eng Comput 1992, 30: 433-440. 10.1007/BF02446182
Holsheimer J: Electrical conductivity of the hippocampal CA1 layers and application to current-source-density analysis. Exp Brain Res 1987, 67: 402-410. 10.1007/BF00248560
Tyner KM, Kopelman R, Philbert MA: "Nanosized voltmeter" enables cellular-wide electric field mapping. Biophys J 2007, 93: 1163-1174. 10.1529/biophysj.106.092452
Durand DM: Electric field effects in hyperexcitable neural tissue: a review. Radiat Prot Dosimetry 2003, 106: 325-331.
Teruel MN, Meyer T: Electroporation-induced formation of individual calcium entry sites in the cell body and processes of adherent cells. Biophys J 1997, 73: 1785-1796. 10.1016/S0006-3495(97)78209-2
Neumann E, Kakorin S, Toensing K: Fundamentals of electroporative delivery of drugs and genes. Bioelectrochem Bioenerg 1999, 48: 3-16. 10.1016/S0302-4598(99)00008-2
This work was supported by CIHR and a Canadian Heart and Stroke Foundation postdoctoral fellowship to Hui Ye. The authors thank Joe Hayek for valuable comments to the paper.
The authors declare that they have no competing interests.
HY was involved with model equation derivation, data analysis, and drafting of the manuscript. MC was involved in generating figures. MGF and PLC supervised and coordinated the study. In addition, MC, EEK, MGF and PLC helped in drafting of the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.