Effects of neuronal magnetic fields on MRI: numerical analysis with axon and dendrite models

Can MRI Be Used as a Sensor to Record Neural Activity?

Magnetic resonance provides exquisite anatomical images and functional MRI monitors physiological activity by recording blood oxygenation. This review attempts to answer the following question: Can MRI be used as a sensor to directly record neural behavior? It considers MRI sensing of electrical activity in the heart and in peripheral nerves before turning to the central topic: recording of brain activity. The primary hypothesis is that bioelectric current produced by a nerve or muscle creates a magnetic field that influences the magnetic resonance signal, although other mechanisms for detection are also considered. Recent studies have provided evidence that using MRI to sense neural activity is possible under ideal conditions. Whether it can be used routinely to provide functional information about brain processes in people remains an open question. The review concludes with a survey of artificial intelligence techniques that have been applied to functional MRI and may be appropriate for MRI sensing of neural activity.

Neural magnetic field dependent fMRI toward direct functional connectivity measurements: A phantom study

Recently, the main issue in neuroscience has been the imaging of the functional connectivity in the brain. No modality that can measure functional connectivity directly, however, has been developed yet. Here, we show the novel MRI sequence, called the partial spinlock sequence toward direct measurements of functional connectivity. This study investigates a probable measurement of phase differences directly associated with functional connectivity. By employing partial spinlock imaging, the neural magnetic field might influence the magnetic resonance signals. Using simulation and phantom studies to model the neural magnetic fields, we showed that magnetic resonance signals vary depending on the phase of an externally applied oscillating magnetic field with non-right flip angles. These results suggest that the partial spinlock sequence is a promising modality for functional connectivity measurements.

Estimation of phase signal change in neuronal current MRI for evoke response of tactile detection with realistic somatosensory laminar network model

Magnetic field generated by neuronal activity could alter magnetic resonance imaging (MRI) signals but detection of such signal is under debate. Previous researches proposed that magnitude signal change is below current detectable level, but phase signal change (PSC) may be measurable with current MRI systems. Optimal imaging parameters like echo time, voxel size and external field direction, could increase the probability of detection of this small signal change. We simulate a voxel of cortical column to determine effect of such parameters on PSC signal. We extended a laminar network model for somatosensory cortex to find neuronal current in each segment of pyramidal neurons (PN). 60,000 PNs of simulated network were positioned randomly in a voxel. Biot-savart law applied to calculate neuronal magnetic field and additional phase. The procedure repeated for eleven neuronal arrangements in the voxel. PSC signal variation with the echo time and voxel size was assessed. The simulated results show that PSC signal increases with echo time, especially 100/80 ms after stimulus for gradient echo/spin echo sequence. It can be up to 0.1 mrad for echo time = 175 ms and voxel size = 1.48 × 1.48 × 2.18 mm(3). With echo time less than 25 ms after stimulus, it was just acquired effects of physiological noise on PSC signal. The absolute value of the signal increased with decrease of voxel size, but its components had complex variation. External field orthogonal to local surface of cortex maximizes the signal. Expected PSC signal for tactile detection in the somatosensory cortex increase with echo time and have no oscillation.

Direct neural current imaging in an intact cerebellum with magnetic resonance imaging

The ability to detect neuronal currents with high spatiotemporal resolution using magnetic resonance imaging (MRI) is important for studying human brain function in both health and disease. While significant progress has been made, we still lack evidence showing that it is possible to measure an MR signal time-locked to neuronal currents with a temporal waveform matching concurrently recorded local field potentials (LFPs). Also lacking is evidence that such MR data can be used to image current distribution in active tissue. Since these two results are lacking even in vitro, we obtained these data in an intact isolated whole cerebellum of turtle during slow neuronal activity mediated by metabotropic glutamate receptors using a gradient-echo EPI sequence (TR=100ms) at 4.7T. Our results show that it is possible (1) to reliably detect an MR phase shift time course matching that of the concurrently measured LFP evoked by stimulation of a cerebellar peduncle, (2) to detect the signal in single voxels of 0.1mm(3), (3) to determine the spatial phase map matching the magnetic field distribution predicted by the LFP map, (4) to estimate the distribution of neuronal current in the active tissue from a group-average phase map, and (5) to provide a quantitatively accurate theoretical account of the measured phase shifts. The peak values of the detected MR phase shifts were 0.27-0.37°, corresponding to local magnetic field changes of 0.67-0.93nT (for TE=26ms). Our work provides an empirical basis for future extensions to in vivo imaging of neuronal currents.

Similar Spectral Power Densities Within the Schumann Resonance and a Large Population of Quantitative Electroencephalographic Profiles: Supportive Evidence for Koenig and Pobachenko

In 1954 and 1960 Koenig and his colleagues described the remarkable similarities of spectral power density profiles and patterns between the earth-ionosphere resonance and human brain activity which also share magnitudes for both electric field (mV/m) and magnetic field (pT) components. In 2006 Pobachenko and colleagues reported real time coherence between variations in the Schumann and brain activity spectra within the 6-16 Hz band for a small sample. We examined the ratios of the average potential differences (~3 μV) obtained by whole brain quantitative electroencephalography (QEEG) between rostral-caudal and left-right (hemispheric) comparisons of 238 measurements from 184 individuals over a 3.5 year period. Spectral densities for the rostral-caudal axis revealed a powerful peak at 10.25 Hz while the left-right peak was 1.95 Hz with beat-differences of ~7.5 to 8 Hz. When global cerebral measures were employed, the first (7-8 Hz), second (13-14 Hz) and third (19-20 Hz) harmonics of the Schumann resonances were discernable in averaged QEEG profiles in some but not all participants. The intensity of the endogenous Schumann resonance was related to the 'best-of-fitness' of the traditional 4-class microstate model. Additional measurements demonstrated real-time coherence for durations approximating microstates in spectral power density variations between Schumann frequencies measured in Sudbury, Canada and Cumiana, Italy with the QEEGs of local subjects. Our results confirm the measurements reported by earlier researchers that demonstrated unexpected similarities in the spectral patterns and strengths of electromagnetic fields generated by the human brain and the earth-ionospheric cavity.

Computational Modeling of Neuronal Current MRI Signals with Rat Somatosensory Cortical Neurons

Magnetic field generated by active neurons has recently been considered to determine location of neuronal activity directly with magnetic resonance imaging (MRI), but controversial results have been reported about detection of such small magnetic fields. In this study, multiple neuronal morphologies of rat tissue were modeled to investigate better estimation of MRI signal change produced by neuronal magnetic field (NMF). Ten pyramidal neurons from layer II to VI of rat somatosensory area with realistic morphology, biophysics, and neuronal density were modeled to simulate NMF of neuronal tissue, from which effects of NMF on MRI signals were obtained. Neuronal current MRI signals, which consist of relative magnitude signal change (RMSC) and phase signal change (PSC), were at least three and one orders of magnitude less than a tissue with single neuron type, respectively. Also, a reduction in voxel size could increase signal alterations. Furthermore, with selection of zenith angle of external main magnetic field related to tissue surface near to 90°, RMSC could be maximized. This value for PSC would be 90° for small voxel size and zero degree for large ones.

Magnetic resonance imaging of ionic currents in solution: the effect of magnetohydrodynamic flow

PURPOSE

Reliably detecting MRI signals in the brain that are more tightly coupled to neural activity than blood-oxygen-level-dependent fMRI signals could not only prove valuable for basic scientific research but could also enhance clinical applications such as epilepsy presurgical mapping. This endeavor will likely benefit from an improved understanding of the behavior of ionic currents, the mediators of neural activity, in the presence of the strong magnetic fields that are typical of modern-day MRI scanners.

THEORY

Of the various mechanisms that have been proposed to explain the behavior of ionic volume currents in a magnetic field, only one-magnetohydrodynamic flow-predicts a slow evolution of signals, on the order of a minute for normal saline in a typical MRI scanner.

METHODS

This prediction was tested by scanning a volume-current phantom containing normal saline with gradient-echo-planar imaging at 3 T.

RESULTS

Greater signal changes were observed in the phase of the images than in the magnitude, with the changes evolving on the order of a minute.

CONCLUSION

These results provide experimental support for the MHD flow hypothesis. Furthermore, MHD-driven cerebrospinal fluid flow could provide a novel fMRI contrast mechanism.

Bloch simulations towards direct detection of oscillating magnetic fields using MRI with spin-lock sequence

A new MRI method using the spin-lock sequence has attracted wide attention because of its potential for detecting small oscillating magnetic fields. However, as the mechanism involved is complicated, we visualized the magnetization performance during the spin-lock sequence in order to better understand interaction of the spin-lock pulse and the externally applied oscillating magnetic fields by means of a fast-and-simple method using matrix operations to solve a time-dependent Bloch equation. To improve spin-lock imaging in the detection of small magnetic fields (in an fMRI experiment that modeled neural magnetic fields), we observed that the phenomenon decreases MR signals, which led us to investigate how spin-lock parameters cause the MR signal to decrease; based on this, we determined that MR signals decrease in oscillating magnetic fields that are resonant with the spin-lock pulse. We also determined that MR signals decrease is directly proportional to spin-lock duration. Our results suggest that MRI can feasibly detect oscillating magnetic fields directly by using of the spin-lock sequence.

Is it possible to detect dendrite currents using presently available magnetic resonance imaging techniques?

The action currents of a dendrite, peripheral nerve or skeletal muscle create their own magnetic field. Many investigators have attempted to detect neural and dendritic currents directly using magnetic resonance imaging that can cause the phase of the spins to change. Our goal in this paper is to use the calculated magnetic field of a dendrite to estimate the resulting phase shift in the magnetic resonance signal. The field produced by a dense collection of simultaneously active dendrites may be just detectable under the most ideal circumstances, but in almost every realistic case the field cannot be detected using current MRI technology.

Direct detection of a single evoked action potential with MRS in Lumbricus terrestris

Functional MRI (fMRI) measures neural activity indirectly by detecting the signal change associated with the hemodynamic response following brain activation. In order to alleviate the temporal and spatial specificity problems associated with fMRI, a number of attempts have been made to detect neural magnetic fields (NMFs) with MRI directly, but have thus far provided conflicting results. In this study, we used MR to detect axonal NMFs in the median giant fiber of the earthworm, Lumbricus terrestris, by examining the free induction decay (FID) with a sampling interval of 0.32 ms. The earthworm nerve cords were isolated from the vasculature and stimulated at the threshold of action potential generation. FIDs were acquired shortly after the stimulation, and simultaneous field potential recordings identified the presence or absence of single evoked action potentials. FIDs acquired when the stimulus did not evoke an action potential were summed as background. The phase of the background-subtracted FID exhibited a systematic change, with a peak phase difference of (-1.2 ± 0.3) × 10(-5) radians occurring at a time corresponding to the timing of the action potential. In addition, we calculated the possible changes in the FID magnitude and phase caused by a simulated action potential using a volume conductor model. The measured phase difference matched the theoretical prediction well in both amplitude and temporal characteristics. This study provides the first evidence for the direct detection of a magnetic field from an evoked action potential using MR.

Modeling magnitude and phase neuronal current MRI signal dependence on echo time

To enhance sensitivity in measuring neuronal current MRI (ncMRI) signal using T(2)*-weighted sequences, appropriate selection of echo time (TE) is vital for optimizing data acquisition strategy. The purpose of this study is to establish the contrast-to-noise ratio of neuronal current MRI signal dependence on TE and determine the optimum TE (TE(opt)) in achieving its highest detection power. The TE(opt) in human brain and tissue preparation at 1.5, 3, and 7 T are estimated with different voxel sizes. Our results show that TE(opt) values are different between magnitude and phase images, and TE(opt) is larger in magnitude than phase imaging. This suggests that a dual-echo data acquisition strategy would provide the best efficiency in detecting magnitude and phase neuronal current MRI signals simultaneously. Our results also indicated that the detection sensitivity will be stronger at lower magnetic fields for human brain, whereas the sensitivity will be enhanced/reduced as field strength increases for phase/magnitude imaging on tissue preparation.

Direct MRI detection of the neuronal magnetic field: the effect of the dendrite branch

In recent years, neuronal current MRI (nc-MRI) was proposed as a new imaging method to directly map the magnetic field change caused by neuronal activity. Nc-MRI could offer improved spatial and temporal resolution compared to blood hemodynamics-based functional magnetic resonance imaging (fMRI). In this paper, with a finite current dipole as the model of dendrite or dendrite branch, we investigated the spatial distribution of the magnetic field generated by synchronously activated neurons to evaluate the possibility of nc-MRI. Our simulations imply that the existence of a dendrite branch may not only increase the strength of the neuronal magnetic field (NMF), but also raise the non-uniform and unsymmetry of the NMF; therefore, it can enhance the detectability of the neuronal current magnetic field by MRI directly. The results show that the signal phase shift is enlarged, but it is unstable and is still very small, <<1 radian, while the magnitude signal may be strong enough for a typical MRI voxel to be detected. We suggest making further efforts to measure the magnitude signal which may induce a large effect in an nc-MRI experiment.

Magnetic resonance imaging of oscillating electrical currents

Functional MRI has become an important tool of researchers and clinicians who seek to understand patterns of neuronal activation that accompany sensory and cognitive processes. However, the interpretation of fMRI images rests on assumptions about the relationship between neuronal firing and hemodynamic response that are not firmly grounded in rigorous theory or experimental evidence. Further, the blood-oxygen-level-dependent effect, which correlates an MRI observable to neuronal firing, evolves over a period that is 2 orders of magnitude longer than the underlying processes that are thought to cause it. Here, we instead demonstrate experiments to directly image oscillating currents by MRI. The approach rests on a resonant interaction between an applied rf field and an oscillating magnetic field in the sample and, as such, permits quantitative, frequency-selective measurements of current density without spatial or temporal cancellation. We apply this method in a current loop phantom, mapping its magnetic field and achieving a detection sensitivity near the threshold required for the detection of neuronal currents. Because the contrast mechanism is under spectroscopic control, we are able to demonstrate how ramped and phase-modulated spin-lock radiation can enhance the sensitivity and robustness of the experiment. We further demonstrate the combination of these methods with remote detection, a technique in which the encoding and detection of an MRI experiment are separated by sample flow or translation. We illustrate that remotely detected MRI permits the measurement of currents in small volumes of flowing water with high sensitivity and spatial resolution.

Can high-field MREIT be used to directly detect neural activity? Theoretical considerations

We sought to determine the feasibility of directly studying neural tissue activity by analysis of differential phase shifts in MRI signals that occurred when trickle currents were applied to a bath containing active or resting neural tissue. We developed a finite element bidomain model of an aplysia abdominal ganglion in order to estimate the sensitivity of this contrast mechanism to changes in cell membrane conductance occurring during a gill-withdrawal reflex. We used our model to determine both current density and magnetic potential distributions within a sample chamber containing an isolated ganglion when it was illuminated with current injected synchronously with the MR imaging sequence and predicted the resulting changes in MRI phase images. This study provides the groundwork for attempts to image neural function using Magnetic Resonance Electrical Impedance Tomography (MREIT). We found that phase noise in a candidate 17.6 T MRI system should be sufficiently low to detect phase signal differences between active and resting membrane states at resolutions around 1 mm(3). We further delineate the broad dependencies of signal-to-noise ratio on activity frequency, current application time and active tissue fractions and outline strategies that can be used to lower phase noise below that presently observed in conventional MREIT techniques. We also propose the idea of using MREIT as an alternative means of studying neuromodulation.

Looking for neuronal currents using MRI: an EEG-fMRI investigation of fast MR signal changes time-locked to frequent focal epileptic discharges

RATIONALE

Reproducible direct measurement of neuronal electrical activity using MRI signal changes due to local magnetic field perturbations would represent a step change in neuroimaging methods. While some previous studies using experiments based on evoked and spontaneous activity provided encouraging results no clear demonstration of neuronal current-related MR changes in the human brain has emerged to date. The availability of simultaneously acquired EEG-fMRI in patients with frequent interictal epileptic discharges (IED), which have significantly greater amplitude than evoked potentials, offers the opportunity to further investigate the phenomenon.

METHODS

We re-analysed simultaneously acquired EEG-fMRI data in 6 epilepsy patients with very frequent focal IED and a well-localised generator. A model of MRI signal changes due to fast activity and BOLD signal changes was used to identify fast MR signal changes, potentially directly reflecting neuronal activity. Simultaneously-acquired EEG allowed the comparison of electrical source localisation (ESI), clinical epilepsy localisation and BOLD signal changes with the fast MR signal changes.

RESULTS

Clusters of IED-related fast MR signal change were observed in all cases. Spatial correspondence between the IED-related fast MR, BOLD, ESI clusters and irritative zone (IZ) was observed in one slice of a single dataset. The other IED-related fast MR clusters were remote from electro-clinically determined generators of interictal activity. The sign and magnitude of the fast MR signal changes varied across regions and subjects.

CONCLUSION

The observed fast MR changes cannot be confidently attributed to the direct effect of neuronal currents due to lack of spatial concordance with generators of interictal activity, IED-related BOLD clusters and ESI estimates.

Physiologically evoked neuronal current MRI in a bloodless turtle brain: detectable or not?

Contradictory reports regarding the detection of neuronal currents have left the feasibility of neuronal current MRI (ncMRI) an open question. Most previous ncMRI studies in human subjects are suspect due to their inability to separate or eliminate hemodynamic effects. In this study, we used a bloodless turtle brain preparation that eliminates hemodynamic effects, to explore the feasibility of detecting visually-evoked ncMRI signals at 9.4 T. Intact turtle brains, with eyes attached, were dissected from the cranium and placed in artificial cerebral spinal fluid. Light flashes were delivered to the eyes to evoke neuronal activity. Local field potential (LFP) and MRI signals were measured in an interleaved fashion. Robust visually-evoked LFP signals were observed in turtle brains, but no significant signal changes synchronized with neuronal currents were found in the ncMRI images. In this study, detection thresholds of 0.1% and 0.1 degrees were set for MRI magnitude and phase signal changes, respectively. The absence of significant signal changes in the MRI images suggests that visually-evoked ncMRI signals in the turtle brain are below these detectable levels.

Comparison of BOLD and direct-MR neuronal detection (DND) in the human visual cortex at 3T

Direct-MR neuronal detection (DND) of transient magnetic fields has recently been investigated as a novel imaging alternative to the conventional BOLD functional MRI (fMRI) technique. However, there remain controversial issues and debate surrounding this methodology, and this study attempts clarification by comparing BOLD responses in the human visual system with those of DND. BOLD relies on indirectly measuring blood oxygenation and flow changes as a result of neuronal activity, whereas the putative DND method is based on the hypothesis that the components of the in vivo neuronal magnetic fields, which lie parallel to the B(0) field, can potentially modulate the MR signal, thus providing a means of direct detection of nerve impulses. Block paradigms of checkerboard patterns were used for visual stimulation in both DND and BOLD experiments, allowing detection based on different frequency responses. This study shows colocalization of some voxels with slow BOLD responses and putative fast DND responses using General Linear Model (GLM) analysis. Frequency spectra for the activated voxel cluster are also shown for both stimulated and control data. The mean percentage signal change for the DND responses is 0.2%, corresponding to a predicted neuronal field of 0.14 nT.

The application of magnets directs the orientation of neurite outgrowth in cultured human neuronal cells

Electric and magnetic fields have been known to influence cellular behavior. In the present study, we hypothesized that the application of static magnetic fields to neurons will cause neurites to grow in a specific direction. In cultured human neuronal SH-SY5Y cells or PC12 cells, neurite outgrowth was induced by forskolin, retinoic acid, or nerve growth factor (NGF). We applied static magnetic fields to the neurons and analyzed the direction and morphology of newly formed neuronal processes. In the presence of the magnetic field, neurites grew in a direction perpendicular to the direction of the magnetic field, as revealed by the higher orientation index of neurites grown under the magnetic field compared to that of the neurites grown in the absence of the magnetic field. The neurites parallel to the magnetic field appeared to be dystrophic, beaded or thickened, suggesting that they would hinder further elongation processes. The co-localized areas of microtubules and actin filaments were arranged into the vertical axis to the magnetic field, while the levels of neurofilament and synaptotagmin were not altered. Our results suggest that the application of magnetic field can be used to modulate the orientation and direction of neurite formation in cultured human neuronal cells.

Stimulus-induced Rotary Saturation (SIRS): a potential method for the detection of neuronal currents with MRI

Neuronal currents produce local transient and oscillatory magnetic fields that can be readily detected by MEG. Previous work attempting to detect these magnetic fields with MR focused on detecting local phase shifts and dephasing in T(2) or T(2)-weighted images. For temporally biphasic and multi-phasic local currents the sensitivity of these methods can be reduced through the cancellation of the accrued phase induced by positive and negative episodes of the neuronal current. The magnitude of the phase shift is also dependent on the distribution of the current within the voxel. Since spins on one side of a current source develop an opposite phase shift relative to those on the other side, there is likely to be significant cancellation within the voxel. We introduce a potential method for detecting neuronal currents though their resonant T(1rho) saturation during a spin-lock preparation period. The method is insensitive to the temporal and spatial cancellation effects since it utilizes the multi-phasic nature of the neuronal currents and thus is not sensitive to the sign of the local field. To produce a T(1)(rho) reduction, the Larmor frequency in the rotating frame, which is set by gammaB(1lock) (typically 20 Hz-5 kHz), must match the major frequency components of the stimulus-induced neuronal currents. We validate the method in MRI phantom studies. The rotary saturation spectra showed a sharp resonance when a current dipole within the phantom was driven at the Larmor frequency in the rotating frame. A 7 min block-design experiment was found to be sensitive to a current dipole strength of 56 nAm, an approximate magnetic field of 1 nT at 1.5 mm from the dipole. This dipole moment is similar to that seen using the phase shift method in a similar experimental setup by Konn et al. [Konn, D., Gowland, P., Bowtell, R., 2003. MRI detection of weak magnetic fields due to an extended current dipole in a conducting sphere: a model for direct detection of neuronal currents in the brain. Magn. Reson. Med. 50, 40-49], but is potentially less encumbered by temporal and spatial cancellation effects.

Microscopic investigation of the resonant mechanism for the implementation of nc-MRI at ultra-low field MRI

In this paper the possible use of the resonant mechanism between some spectral components of the neuronal activity and the spin dynamics in ultra-low field MRI experiments--for the implementation of the nc-MRI techniques and proposed by Kraus et al., 2008--is investigated by means of "realistic" simulations of the neuronal activity of a modelled neuronal network. Previously characterized digital neurons are used to reproduce neuronal currents based on biophysical details and the distribution of the local magnetic field inside a MRI cubic voxel (having a dimension of 1.2 mm) is evaluated. The properties of the water proton spin dynamics as a consequence of the neuronal field and of external applied fields are extrapolated integrating the Bloch equations. The characteristics of the expected MR signals are discussed in relation to the specifics of the NMR sequence used and to the properties of the neuronal activity. The great potentialities of the technique are provided by: a) the possible easy implementation of the technique, b) the possible cheap instrumentation required; c) the flexibility of the ultra-low field systems.

Lorentz effect imaging of ionic currents in solution

Current functional MRI techniques relying on hemodynamic modulations are inherently limited in their ability to accurately localize neural activity in space and time. To address these limitations, we previously proposed a novel technique based on the Lorentz effect and demonstrated its ability to directly image minute electrical activity with a millisecond temporal resolution in gel phantoms containing conductive wires as well as in the human median nerve in vivo. To better characterize its contrast mechanism and ultimately further improve its sensitivity for in vivo applications, we now apply this technique to image ionic currents in solution, which serve as a better model for neural conduction in biological systems than the electronic currents in conductive wires used in previous phantom studies. Our results demonstrate that ionic currents with durations and current densities on the same order of magnitude as those induced by neuroelectric activity in nerve fibers and in the brain can be detected.

Realistic simulations of neuronal activity: a contribution to the debate on direct detection of neuronal currents by MRI

Many efforts have been done in order to preview the properties of the magnetic resonance (MR) signals produced by the neuronal currents using simulations. In this paper, starting with a detailed calculation of the magnetic field produced by the neuronal currents propagating over single hippocampal CA1 pyramidal neurons placed inside a cubic MR voxel of length 1.2 mm, we proceeded on the estimation of the phase and magnitude MR signals. We then extended the results to layers of parallel and synchronous similar neurons and to ensembles of layers, considering different echo times, voxel volumes and neuronal densities. The descriptions of the neurons and of their electrical activity took into account the real neuronal morphologies and the physiology of the neuronal events. Our results concern: (a) the expected time course of the MR signals produced by the neuronal currents in the brain, based on physiological and anatomical properties; (b) the different contributions of post-synaptic potentials and of action potentials to the MR signals; (c) the estimation of the equivalent current dipole and the influence of its orientation with respect to the external magnetic field on the observable MR signal variations; (d) the size of the estimated neuronal current induced phase and magnitude MR signal changes with respect to the echo time, voxel-size and neuronal density. The inclusion of realistic neuronal properties into the simulation introduces new information that can be helpful for the design of MR sequences for the direct detection of neuronal current effects and the testing of bio-electromagnetic models.