Florida State University Libraries
Electronic Theses, Treatises and Dissertations The Graduate School
2012 Cellular Responses to Osmotic Perturbation: A High-Field ¹H and ²³Na Magnetic Resonance Microscopy Study John Walsh
Follow this and additional works at the FSU Digital Library. For more information, please contact [email protected] THE FLORIDA STATE UNIVERSITY
COLLEGE OF ENGINEERING
CELLULAR RESPONSES TO OSMOTIC PERTURBATION:
A HIGH-FIELD 1H AND 23Na MAGNETIC RESONANCE MICROSCOPY STUDY
By
JOHN WALSH
A Thesis submitted to the Department of Chemical and Biomedical Engineering in partial fulfillment of the requirements for the degree of Master of Science
Degree Awarded: Fall Semester, 2012 John Walsh defended this thesis on August 10, 2012.
The members of the supervisory committee were:
Samuel C. Grant Professor Directing Thesis
Teng Ma Committee Member
Jingjiao Guan Committee Member
The Graduate School has verified and approved the above-named committee members, and certifies that the thesis has been approved in accordance with university requirements.
ii
This thesis is dedicated to my family and friends who have continued to support me through this long process. I couldn’t have done it without you.
iii
ACKNOWLEDGEMENTS
I would first like to acknowledge the FAMU-FSU College of Engineering and the National High Magnetic Field Laboratory for use of their facilities to perform this work. Only with access to these high magnetic fields could a project of this nature have been undertaken. I would also like to acknowledge my supervising professor, Dr. Samuel C. Grant, as his knowledge of magnetic resonance imaging and RF coil design were critical in conducting the experiments and acquiring the high resolution images obtained. Additionally, my committee members Dr. Teng Ma and Dr. Jingjiao Guan provided much insight on the directions and results of this project. I would like to acknowledge Parastou Foroutan, a graduate student in Dr. Grant’s lab, who aided me with the initial sodium imaging and acquainted me with the techniques used throughout this project. Additionally, I would like to acknowledge Corey Falgas, another fellow graduate student, who aided in maintaining the animals used in this work. I would also like to thank my other colleagues Ihssan Masad, Jens Rosenberg, and Jose Muniz for their support, ideas, and help throughout the course of the project. Finally, funding for this project came from the American Heart Association (SE Division) and the FSU Department of Chemical and Biomedical Engineering. Additional funding came from the FSU Office of National Fellowships through an Undergraduate Research and Creative Activity Award. For all of these contributions, I am grateful.
iv
TABLE OF CONTENTS
List of Tables ...... vii List of Figures ...... ix Abstract ...... xiii 1. INTRODUCTION ...... 1 2. BACKGROUND AND THEORY ...... 4 2.1 Magnetic Resonance Imaging ...... 4 2.2 Cell Membranes and Osmosis ...... 12 2.3 Relation to Pathological Conditions ...... 14 2.4 Single Cell Models ...... 17 2.5 Tissue Models ...... 20 2.6 Diffusion Multi-Compartment Modeling...... 22 2.7 Sodium MRI ...... 24 3. OBJECTIVES ...... 28 3.1 High-Field MRI ...... 28 3.2 Radiofrequency (RF) Coil Design ...... 30 3.3 Neural Ganglia Model System ...... 32 3.4 Osmotic Perturbation ...... 33 3.5 Cell Viability ...... 34 4. METHODS ...... 36 4.1 Neural Ganglia Preparation ...... 36 4.2 Artificial Sea Water (ASW) Solutions ...... 38 4.3 MR Image Acquisition ...... 43 4.4 Statistical Analysis...... 50 5. RESULTS ...... 52 5.1 1H/23Na Imaging of Fixed Ganglia ...... 52 5.1.1 Isotonic Fixed Ganglia ...... 52 5.1.2 Hypertonic Fixed Ganglia ...... 60 5.1.3 Hypotonic Fixed Ganglia...... 67 5.1.4 Fixed Ganglia Summary ...... 74 5.2 1H/23Na Imaging of Viable and Nonviable Ganglia ...... 76 5.2.1 Isotonic Viable/Nonviable Ganglia ...... 76 5.2.2 Hypertonic Viable/Nonviable Ganglia ...... 85 5.2.3 Hypotonic Viable/Nonviable Ganglia ...... 93 5.2.4 Viable/Nonviable Ganglia Summary ...... 102 5.3 Comparison of Viable and Fixed Ganglia ...... 106
v
6. DISCUSSION ...... 111 7. CONCLUSIONS AND FUTURE DIRECTIONS ...... 118 References ...... 121 Biographical Sketch ...... 135
vi
LIST OF TABLES
1 Standard concentrations of key compounds in isotonic artificial sea water (ASW) ...... 38
2 Quantities of each component to make 1 L of ASW solutions of different tonicities ...... 42
1 3 Summary of MR acquisition parameters for the H spin echo T1 measurement ...... 45
1 4 Summary of MR acquisition parameters for the H spin echo T2 measurement ...... 46
1 5 Summary of MR acquisition parameters for the H gradient echo T2* measurement ...... 47
6 Summary of MR acquisition parameters for the 1H diffusion-weighted spin echo ADC measurement ...... 48
23 7 Summary of MR acquisition parameters for the Na 3D gradient echo T1 measurement ... 49
23 8 Summary of MR acquisition parameters for the Na 3D gradient echo T2* measurement ...... 50
9 Proton and sodium relaxation and diffusion measurements for the isotonic fixed ganglia .. 60
10 Proton and sodium relaxation and diffusion measurements for the hypertonic fixed ganglia ...... 67
11 Proton and sodium relaxation and diffusion measurements for the hypotonic fixed ganglia ...... 74
12 Summary of relaxation and diffusion trends in fixed ganglia at varying tonicities ...... 75
13 Proton relaxation and diffusion measurements for the isotonic viable ganglia ...... 84
14 Proton and sodium relaxation and diffusion measurements for the isotonic nonviable ganglia ...... 84
15 Proton relaxation and diffusion measurements for the hypertonic viable ganglia ...... 92
16 Proton and sodium relaxation and diffusion measurements for the hypertonic nonviable ganglia ...... 93
17 Proton relaxation and diffusion measurements for the hypotonic viable ganglia ...... 101
18 Proton and sodium relaxation and diffusion measurements for the hypertonic nonviable ganglia ...... 102
19 Summary of relaxation and diffusion trends in viable ganglia at varying tonicities ...... 104 vii
20 Summary of relaxation and diffusion trends in nonviable ganglia at varying tonicities .... 105
viii
LIST OF FIGURES
1 The orientation of magnetic moments when placed in an external magnetic field ...... 5
2 Diagram of a standard spin echo pulse sequence ...... 8
3 Representation of the net magnetization during a spin echo pulse sequence ...... 9
4 Diagram of a standard diffusion-weighted spin echo pulse sequence ...... 11
5 The fluid mosaic model of the cell membrane ...... 12
6 Mechanism of the sodium-potassium ATPase ...... 13
7 Mechanism of cell death during ischemia and stroke ...... 15
8 Diagram of the mechanism of osmotic redistributions that result due to ischemia and stroke ...... 16
9 Dissection process for isolating the abdominal ganglion from Aplysia ...... 37
10 Image of three viable isolated abdominal ganglia from Aplysia after dissection ...... 37
11 Cellular response of placing neuronal ganglia in extracellular environments of different tonicities ...... 43
12 Image of the solenoidal transmit-receive coil showing both resonant circuits for simultaneous tuning to both the proton and sodium resonant frequencies ...... 44
1 13 Representative T1-weighted H MR images of three fixed ganglia in isotonic ASW acquired at 11.7 T ...... 54
1 14 Representative T2-weighted H MR images of three fixed ganglia in isotonic ASW acquired at 11.7 T ...... 55
1 15 Representative T2*-weighted H MR images of three fixed ganglia in isotonic ASW acquired at 11.7 T ...... 56
16 Representative diffusion-weighted 1H MR images of three fixed ganglia in isotonic ASW acquired at 11.7 T ...... 57
17 Representative 23Na MR images of three fixed ganglia in isotonic ASW acquired at 11.7 T ...... 58
18 Representative regression curve fits for the isotonic fixed ganglia ...... 59
ix
1 19 Representative T1-weighted H MR images of three fixed ganglia in hypertonic ASW acquired at 11.7 T ...... 61
1 20 Representative T2-weighted H MR images of three fixed ganglia in hypertonic ASW acquired at 11.7 T ...... 62
1 21 Representative T2*-weighted H MR images of three fixed ganglia in hypertonic ASW acquired at 11.7 T ...... 63
22 Representative diffusion-weighted 1H MR images of three fixed ganglia in hypertonic ASW acquired at 11.7 T ...... 64
23 Representative 23Na MR images of three fixed ganglia in hypertonic ASW acquired at 11.7 T ...... 65
24 Representative regression curve fits for the hypertonic fixed ganglia ...... 66
1 25 Representative T1-weighted H MR images of three fixed ganglia in hypotonic ASW acquired at 11.7 T ...... 68
1 26 Representative T2-weighted H MR images of three fixed ganglia in hypotonic ASW acquired at 11.7 T ...... 69
1 27 Representative T2*-weighted H MR images of three fixed ganglia in hypotonic ASW acquired at 11.7 T ...... 70
28 Representative diffusion-weighted 1H MR images of three fixed ganglia in hypertonic ASW acquired at 11.7 T ...... 71
29 Representative 23Na MR images of three fixed ganglia in hypertonic ASW acquired at 11.7 T ...... 72
30 Representative regression curve fits for the hypotonic fixed ganglia ...... 73
1 31 Representative T1-weighted H MR images of three viable ganglia in isotonic ASW acquired at 11.7 T ...... 78
1 32 Representative T2-weighted H MR images of three viable ganglia in isotonic ASW acquired at 11.7 T ...... 79
1 33 Representative T2*-weighted H MR images of three viable ganglia in isotonic ASW acquired at 11.7 T ...... 80
34 Representative diffusion-weighted 1H MR images of three viable ganglia in isotonic ASW acquired at 11.7 T ...... 81
x
35 Representative 23Na MR images of three nonviable ganglia in isotonic ASW acquired at 11.7 T ...... 82
36 Representative regression curve fits for the isotonic viable ganglia...... 83
1 37 Representative T1-weighted H MR images of three viable ganglia in hypertonic ASW acquired at 11.7 T ...... 86
1 38 Representative T2-weighted H MR images of three viable ganglia in hypertonic ASW acquired at 11.7 T ...... 87
1 39 Representative T2*-weighted H MR images of three viable ganglia in hypertonic ASW acquired at 11.7 T ...... 88
40 Representative diffusion-weighted 1H MR images of three viable ganglia in hypertonic ASW acquired at 11.7 T ...... 89
41 Representative 23Na MR images of three nonviable ganglia in hypertonic ASW acquired at 11.7 T ...... 90
42 Representative regression curve fits for the hypertonic viable ganglia ...... 91
1 43 Representative T1-weighted H MR images of three viable ganglia in hypotonic ASW acquired at 11.7 T ...... 95
1 44 Representative T2-weighted H MR images of three viable ganglia in hypotonic ASW acquired at 11.7 T ...... 96
1 45 Representative T2*-weighted H MR images of three viable ganglia in hypotonic ASW acquired at 11.7 T ...... 97
46 Representative diffusion-weighted 1H MR images of three viable ganglia in hypotonic ASW acquired at 11.7 T ...... 98
47 Representative 23Na MR images of three nonviable ganglia in hypotonic ASW acquired at 11.7 T ...... 99
48 Representative regression curve fits for the hypotonic viable ganglia ...... 100
49 Comparisons of fixed, viable, and nonviable ganglia between hypotonic, isotonic, and 1 hypertonic extracellular tonicities for the H T1 measurement ...... 107
50 Comparisons of fixed, viable, and nonviable ganglia between hypotonic, isotonic, and 1 hypertonic extracellular tonicities for the H T2 measurement ...... 107
xi
51 Comparisons of fixed, viable, and nonviable ganglia between hypotonic, isotonic, and 1 hypertonic extracellular tonicities for the H T2* measurement ...... 108
52 Comparisons of fixed, viable, and nonviable ganglia between hypotonic, isotonic, and hypertonic extracellular tonicities for the 1H ADC measurement ...... 108
53 Comparisons of fixed, viable, and nonviable ganglia between hypotonic, isotonic, and 23 hypertonic extracellular tonicities for the Na T1 measurement ...... 109
54 Comparisons of fixed, viable, and nonviable ganglia between hypotonic, isotonic, and 23 hypertonic extracellular tonicities for the Na T2* measurement ...... 109
55 Comparisons of 23Na SNR in nonviable ganglia between hypotonic, isotonic, and hypertonic extracellular tonicities ...... 110
56 23Na images and signal intensity of a single viable ganglion in isotonic, hypertonic, and hypotonic ASW ...... 110
xii
ABSTRACT
Magnetic resonance imaging (MRI) relaxation and diffusion properties are sensitive to the physiological state of cells and tissues. In particular, cerebral ischemia and stroke is associated with a reduced apparent diffusion coefficient (ADC) and increased transverse 1 magnetization relaxation parameter (T2) in proton ( H) MRI making it useful for clinical diagnosis. Further, increases in sodium (23Na) signal intensity in 23Na MRI correlate with several neurodegenerative diseases as well as ischemia and stroke, which can be used to identify stroke lesions and predict stroke onset time. However, the contributing mechanisms underlying these 23 changes, including possible increases in intracellular sodium and/or lengthened Na T2 relaxation, are not well understood. It has been hypothesized that alterations in cell regulatory mechanisms result in cell swelling, which disrupts tissue microstructure and ionic distributions. Osmotic perturbations have been used on single neurons and neural tissue models to mimic these volume changes. Evaluated with 1H MRI, hypotonic perturbations mimic tissue ischemia and result in an increase in 1H signal intensity, a decrease in 1H ADC, and an increase 1 in H T2, while all trends are reversed for hypertonic perturbation. Single cell work has focused on the large L7 neuron from the abdominal ganglion in the sea hare Aplysia californica, which has been shown to have distinct nuclear and cytoplasmic compartments with differing relaxation and diffusion properties. Therefore, it would be useful to develop a tissue model comprised of large cells in which the intracellular contribution to the volume averaged signal could be determined. In this study, Aplysia abdominal ganglia were used due to their simple anatomy with a small collection of relatively large neurons up to 300 microns in diameter representing a simple neural tissue model system. The abdominal ganglia were dissected from the living animal, washed with isotonic artificial sea water (ASW; in mM: 460 NaCl, 10.4 KCl, 55 MgCl2, 11 CaCl2, 15 HEPES) and loaded into a 2.5-mm o.d. capillary containing isotonic, hypertonic, or hypotonic ASW. Hypertonic and hypotonic perturbations were introduced by changing the isotonic sodium chloride concentration (460 mM) to 545 mM and 345 mM, respectively. All MR imaging was performed at 11.75 T utilizing a homebuilt, double-tuned 1H/23Na solenoidal coil having a diameter of 3 mm. Imaging was performed in separate studies to quantify the proton T1,
T2, T2*, and ADC immediately following dissection on viable ganglia as well as sodium T1 and
T2* after the loss of cell viability. Changes in the MRI relaxation parameters (T1, T2 and T2*) as xiii well as the apparent diffusion coefficient (ADC) in isolated ganglia were used to assess changing cellular environments under osmotic perturbation in the context of the influence changing ionic distributions, notably sodium, have on cell swelling in disease states. With tonicity changes, the ganglionic sac remains intact although alterations in the signal intensity (particularly with respect to sodium) are evident. Proton T1 did not change significantly 1 1 with osmotic perturbation; however, H T2 and T2* decreased with increasing tonicity and H
T2* provided the most robust measurement for identifying cellular changes due to osmotic perturbation. All relaxation parameters increased with cell death and loss of viability and no changes in 1H ADC were observed. The trends observed in 1H relaxation are possibly due to an increased intracellular volume fraction while the lack of a clear trend in ADC values may be due to significant volume averaging between the nuclear, cytoplasmic, and interstitial environments and should be modeled utilizing a multi-compartmental approach. Sodium measurements reveal a general increasing trend in relaxation with osmotic perturbation but no clear conclusions could be drawn. Changes in sodium relaxation are evident in compromised ganglia; however, linearly increasing trends in sodium signal intensity are seen in all osmotic conditions and may represent changes occurring with loss of cell viability, providing evidence that the direct evaluation of ions may prove more sensitivity in assessing pathological osmotic changes. Future work with perfused ganglia as well as voxel-selective spectroscopy and contrast agents to isolate the intracellular sodium signal could be used to answer questions regarding these observed trends in sodium signal intensity and relaxation.
xiv
CHAPTER ONE
INTRODUCTION
All living cells maintain a delicate homeostatic balance of intracellular and extracellular osmolytes which drives transport across the cell membrane and is critical in regulating cell volume (Lang, et al. 1998). Ions, such as sodium, potassium, and chloride, as well as cellular proteins, maintain these osmotic gradients. In fact, a large percentage of the cell’s energy supply is expended in regulating cell volume by actively transporting sodium and other ions across the cell membrane to influence osmosis, the passive diffusion of water to regions of higher osmolyte concentration (Faller, 2008). Transport between the intracellular and extracellular environments is critical to the nourishment and health of cells. Additionally, ion distributions across the membrane create electrochemical gradients which act as the basis of signal transduction in neurons and myocytes. Often, external perturbations, both osmotic and neurotoxic, interfere with these transport systems resulting in cell volume and osmoregulation imbalances (Carpenter, et al. 1992). This loss of cellular homeostasis becomes apparent in pathological disease states such as stroke, cancer, and chronic neurodegeneration. Specifically, protein synthesis and degradation is altered in states of osmotic stress leading to protein distribution alterations that are found in many common neurodegenerative diseases such as Alzheimer’s, Parkinson’s and Huntington’s disease (Mao, et al. 2008). Further, the osmotic effects of water redistributions and changing ionic distributions are evident in states of tissue ischemia and acute stroke (Young, et al. 1987). Intricate osmoregulation pathways exist in organisms, including humans, that respond to states of osmotic changes in the extracellular environment through feedback mechanisms in an attempt to maintain cellular homeostasis (Bourque, 2008). Therefore, methods to quantify tissue microstructure, changes in cell volume, and the tissue sodium concentration in these pathological states are needed and have led to the incorporation of magnetic resonance microscopy as a tool to noninvasively evaluate changing tissue microstructure at sub-millimeter resolutions. Further, MR microscopy provides the added benefit of sensitivity to physiological function which, through different contrast mechanisms such as relaxation and diffusion weighting, can characterize the state and distribution of cellular water as well as relate restricted water diffusion to changing tissue microstructures in vivo (Beneviste & Blackband, 2002). Expanding the MR techniques to non-proton nuclei such as
1 sodium (23Na) allows the tissue sodium concentration and ionic distributions to be measured in abnormal tissues, providing another biomarker for many pathological conditions (Christensen, et al. 1996). By measuring intracellular and extracellular sodium and monitoring the rate of transport across the neural membrane through cell volume changes, a wealth of information about osmoregulation can be investigated in light of chemical and biological challenges to cell viability and function (Pasantes-Morales & Tuz, 2006). In addition, differences in microstructure and ion and protein concentration between the nucleus and cytoplasm are vital to intracellular equilibrium and provide a basis in understanding cellular responses to osmotic perturbation which may be expanded to the analysis of small collections of cells and tissues (Hsu, et al. 1996). Magnetic resonance imaging (MRI) provides a powerful tool for measuring sodium and water distributions in vivo. The correlation of MRI relaxation parameters and diffusion coefficients with the cellular response of changing intracellular volume fractions under osmotic (O’Shea, et al. 2000) and neurotoxic (Buckley, et al. 1998) perturbation allows the assessment of physiological conditions in disease states and is a direct indicator of changes in cell structure and water distribution involved in cell volume regulation. In the clinical setting, stroke and cerebral ischemia are associated with a reduced apparent diffusion coefficient (ADC) and increased 1 transverse magnetization relaxation parameter (T2) in proton ( H) MRI making it useful for diagnosis (Baird & Warach, 1998). Although usually applied to 1H imaging, the technique can be advanced to include 23Na MRI, which facilitates the direct measurement of sodium distributions in tissues, allowing the role of sodium concentration gradients and the direct maintenance of ionic sodium distributions in the regulation of cell volume to be examined in pathological conditions such as cerebral ischemia (Lin, et al. 2001). 23Na MRI can also be used to identify stroke lesions, but unlike 1H MRI, the sodium signal intensity increases in a linear fashion making it possible to predict stroke onset time which allows more appropriate treatment options to be pursued (Jones, et al. 2006). Further, increases in sodium signal intensity are evident in brain tumors and recently demonstrated in chronic neurodegeneration, opening the field of 23Na MRI into the investigation of a variety of pathologies (Thulborn, et al. 1999; Mellon, et al. 2009; Inglese, et al. 2010). MR diffusion and relaxation, as well as sodium signal intensity, measurements are aimed at examining differences in parameter values between normal and perturbed conditions in the context of understanding cells in diseased states.
2
The typical MRI relaxation parameters, T1, T2, and T2*, and the apparent diffusion coefficient (ADC) from diffusion-weighted magnetic resonance imaging are used to quantitatively assess changes in tissue microstructure and water distributions. Further, 23Na MRI provides physiological information on ionic concentrations and distributions within different tissue regions. In particular, changes in sodium signal and relaxation show a faster response to perturbation, prior to morphological change at the tissue level, indicating the benefit of this imaging modality in the early detection of pathological conditions (Schepkin, et al. 2011). Further insight can be gained through the incorporation of sodium MRI as a tissue viability marker and illustrating how these techniques can be expanded and incorporated into clinical image settings (Hussain, et al. 2009). 1 This study will be aimed at examining differences in H relaxation (T1, T2, and T2*) and 23 diffusion (ADC) as well as Na relaxation (T1 and T2*) in the abdominal neural ganglion from the sea hare Aplysia californica. As a small collection of relatively large cells, the abdominal ganglion will serve as a model system previously not characterized with MR imaging techniques. Further, as a majority of the work performed on single cells using MR involves the Aplysia L7 neuron from the abdominal ganglion, the characterization of the entire ganglion will provide insight between previous single cell studies and the more complex tissue models that have been explored. Although the measured relaxation and diffusion properties are an average of the response of the nucleus, cytoplasm, interstitial space, and extracellular medium, comparisons of cellular responses between different extracellular tonicities can provide insight into the collective ganglia response to osmotic stress and the nature of the single cell response within the tissue environment. Further, by imaging fixed, viable, and nonviable ganglia, MR changes that are characteristic of determining cell viability can be described. The results are presented in the context of understanding the mechanism of the MR changes observed in tissue ischemia and chronic neurodegeneration models.
3
CHAPTER TWO
BACKGROUND AND THEORY
2.1 Magnetic Resonance Imaging
Magnetic Resonance Imaging (MRI) is a non-ionizing medical imaging tool utilized for its tissue contrast and three-dimensional imaging capabilities. The use of MRI in the field of medical imaging was not introduced until the 1970s when it was shown that normal and abnormal tissue behaved differently when placed in a magnetic field. The increase in NMR relaxation rates (T1 and T2) in tumors as compared to normal tissues showed that this principle could be used to differentiate between tissues in disease states (Damadian, 1971). Since that time, MRI has revolutionized the medical community and become a prominent technique in the detection and diagnosis of disease. Although an insensitive imaging modality limited in spatial resolution on the order of 1-3 mm with a temporal resolution on the order of 10 minutes in clinical applications, the use of higher magnetic fields and advanced radiofrequency coil design has allowed the realm of magnetic resonance imaging to expand to resolutions on the order of tens of microns. Early from the conception of MRI, its use as a microscopic tool was recognized and currently resolutions less than 100 µm are feasible in research settings, allowing MRI to be termed magnetic resonance microscopy (Blackband, et al. 1999). The use of high-resolution imaging has allowed the cellular components of isolated tissues and single cells to be characterized through MRI relaxation and diffusion measurements. The principles of nuclear magnetism govern MRI, such that atomic nuclei with an odd atomic mass or atomic number, which results from the number and distribution of protons and neutrons in the nucleus, have a nonzero spin quantum number (denoted I). The nuclei are charged and, due to the nonzero spin, have a net angular momentum that causes them to act like magnets. Therefore, atomic nuclei which have a net nuclear spin are sensitive to magnetic resonance imaging techniques (Webb, 2003). Although clinically applied to hydrogen atoms in water molecules with a spin of 1/2 (I = 1/2), MRI can be expanded to include quadrupolar nuclei such as sodium which has a spin of 3/2 (I = 3/2). With a weaker signal emission due to its chemical nature and lower cellular concentration, sodium is much more difficult to image (Lin, et al. 2001).
4
The magnetic moment (denoted µ) of these nuclei can be related to the spin quantum number such that the magnetic moment is quantized and can only take certain discrete values as seen in equation (1) where h is Planck’s constant and γ is the gyromagnetic ratio, the ratio of the nuclear magnetic moment to the angular momentum.
(1) ⁄ | | [ ]
When placed in an external magnetic field, the nuclei align with the field; however, although the nuclei preferentially align in a low energy state parallel to the direction of the external magnetic field, there is a distribution between the low energy parallel state and the higher energy anti- parallel state as seen in Figure 1. Although illustrated as being directly aligned with the field, the magnetic moments are only partially aligned as they can only have discrete values. Also, the split into two energy levels only applies for spin 1/2 nuclei, corresponding to the nuclear magnetic quantum number which can have values of 1/2 and -1/2. For spin 3/2 nuclei, four energy levels are allowed and the magnetic moments will be distributed between the energy levels corresponding to magnetic quantum numbers of 3/2, 1/2, -1/2, and -3/2.
Figure 1. The orientation of magnetic moments when placed in an external magnetic field. Left: Without a magnetic field, the magnetic moments are random. Right: With the application of a magnetic field, the magnetic moments are aligned in a parallel or anti-parallel state resulting in a net magnetization (van Geuns, et al. 1999).
5
Each atomic nucleus has an energy related to its inherent magnet moment, gyromagnetic ratio, and the strength of the external magnetic field. For example, the proton nucleus has an energy difference between the parallel and antiparallel orientations shown in equation (2) where B0 is the strength of the externally applied magnetic field.
(2)
Based on the energy differences at a given temperature (T), there is a given distribution of spins between the parallel and antiparallel orientations as shown in equation (3) where k is the Boltzmann constant. The difference is small and represents the relative insensitivity of the MR imaging technique as it can only detect the differences between these states.
(3)
Upon being aligned with the magnetic field, each nucleus spins about the axis of the main magnetic field at a resonant frequency, often called the Larmor frequency, which is characteristic of each atomic nucleus and the external magnetic field strength as shown in equation (4). The gyromagnetic ratio (γ/2π) for the proton nucleus is 42.58 MHz/T and for the sodium nucleus is 11.26 MHz/T.
(4)
For NMR applications, the orientations of the nuclear magnetic moments are altered. Upon excitation with a radio frequency (RF) pulse at the resonant frequency, the spins absorb energy and are reversed into the high energy anti-parallel state. At the completion of the RF pulse, the net magnetization begins to return to the equilibrium state. After excitation, the nuclear spins relax to their previous thermal equilibrium state such that nuclear spins that were flipped into the anti-parallel state return to the parallel state, releasing energy that is detected as an induced voltage in a radio-frequency receiver coil. The Bloch equations represent the 6 mathematical formulation of the relaxation of the net magnetization after being altered from the equilibrium state in the x direction (5), y direction (6), and z direction (8). By definition, the z direction is aligned with the externally applied magnetic field. Before excitation, the net magnetization is only in the z direction; however, to obtain any MR signal, the net magnetization must be rotated into the x-y plane.
(5)
( )
(6)
( )
(7)
In particular, several characteristic relaxation parameters are used to assess the return of the net magnetization to equilibrium. The T1, or spin-lattice, relaxation parameter is the longitudinal relaxation in the direction of the externally applied magnetic field. The T1 is directly related to the correlation time and represents the thermal exchange of energy between the nuclear spins and the surrounding environment. The T1 relaxation parameter can be detected through a spin-echo pulse sequence that introduces a 90 degree RF pulse, rotating the net magnetization along the z axis into the x-y plane and observing the time response to the initial state prior to excitation (Kuhn, 1990). A typical spin echo pulse sequence is illustrated in Figure 2. The repetition times (TR) govern the amount of time the longitudinal magnetization has to relax back to equilibrium values before the next excitation. The repetition time, TR, or time between successive pulse sequences is varied and plotted against the mean signal detected in an exponential rise to max equation that takes the form of equation (8).
(8)
( )
7
The T1 relaxation parameter is also a measure of viscosity as more viscous substances have shorter T1 relaxation times as there is a more effective transfer of energy between nuclei and their environment.
Figure 2. Diagram of a standard spin echo pulse sequence. The echo time (TE) is the time between the 90˚ excitation pulse and the image acquisition. The repetition time (TR) is the time between successive 90˚ excitation pulses. The phase and frequency encoding gradients spatially encode the two-dimensional image within the slice encoded by the slice-selecting gradient (van Geuns, et al. 1999).
The T2, or spin-spin, relaxation parameter is the transverse relaxation in the direction perpendicular to the applied magnetic field and is due to a loss of phase coherence of oscillating nuclei. Also measured through a spin echo pulse sequence, inhomogeneities in the main magnetic field cause the transverse magnetization to disperse. Through the introduction of a 90 degree pulse followed by a 180 degree pulse, the net transverse magnetization is realigned into a characteristic echo response which can be detected as illustrated in Figure 3 (Kuhn, 1990). The echo time, TE, or time between initial excitation and signal measurement is varied and plotted against the mean signal in an exponentially decreasing function that takes the form of the equation (9) (Webb, 2003).
(9)
8
Figure 3. Representation of the net magnetization during a spin echo pulse sequence. The coordinate system shown is rotating at the resonant frequency of the nuclei. A) Application of the 90˚ pulse flips the magnetization into the transverse plane. B) The spins begin to diphase due to magnetic field inhomogeneities. C) Application of the 180˚ pulse flips the magnetization within the transverse plane and dephasing continues. D) At the echo time (TE) the spins rephase producing the spin echo (Kuhn, 1990).
Based on T2 relaxation, the T2* relaxation parameter is also the transverse relaxation but differs from T2 based on sample inhomogeneities and magnetic susceptibilities. Particularly in samples that are composed of regions that have different magnetic susceptibilities, such as in tissues, there are greater signal losses at the boundaries and is especially pronounced at interfaces with air, such as an air bubble located within the sample. Additionally, the T2* relaxation parameter is determined in a very different manner. Using a gradient echo pulse sequence, the characteristic echo response is generated not by using a 180 degree refocusing pulse, but rather through the application of a bipolar gradient that cancels phase differences that result from sample inhomogeneities and variations in the externally applied magnetic field. Again, the echo
9 time, TE, or time between initial excitation and signal measurement is varied and plotted against the mean signal in an exponentially decreasing function that takes the form of the equation (10).
(10)
The apparent diffusion coefficient (ADC) is a measure of nuclear movements within a sample. Based on the spin echo pulse sequence used for the determination of T1 and T2, an additional diffusion encoding gradient is introduced before and after the 180 degree pulse as illustrated in Figure 4. Nuclei that move within this time period are not seen in the resulting image which can be used to differentiate between diffusing and non-diffusing nuclei. The strength of the diffusion gradient, also termed the b values, is varied and an exponentially decreasing function results, as shown in equation (11) where the b values are defined in equation (12) and the ADC is the apparent diffusion coefficient (Minati &Weglarz, 2007). The strength (G), duration (δ) and delay (Δ) correspond to the pulse sequence in Figure 4.
(11)
(12)
( )
Diffusion measurements in inhomogeneous environments, such as tissues, are complicated by barriers to diffusion that impact the ADC. The theoretical models for how these factors (restrictions to due cellular membranes, diffusion anisotropy, diffusion times, and impact of specific experimental setup) contribute to diffusion measurements have been explored in some detail (Szafer et al. 1995; Yablonskiy & Sukstanskii, 2010). Factors such as a cell size distribution (Jespersen, et al. 2005) were shown to not significantly impact diffusion measurements while increases in cellularity were shown to decrease the ADC (Matsumoto, et al. 2009). Further, temperature impacts both relaxation (Kamman, et al. 1988) and diffusion (Matsumoto, et al. 2009) measurements. Frequently diffusion measurements result in biexponential signal decay resulting from two different populations of water nuclei (Minati &
10
Weglarz, 2007) such that exchange must be considered (Latt, et al. 2009). However, the biexponential behavior cannot be directly related to compartmentation found in tissues, and although the model fits experimental data well, it is not founded on physical barriers to diffusion (Kiselev & Il’yasov, 2007; Schwarz, et al. 2004; Sehy, et al. 2002). Further, the diffusion coefficient measured in MRI differs from the intrinsic dynamic ADC characterized from mean- square displacements (Grebenkov, 2010).
Figure 4. Diagram of a standard diffusion-weighted spin echo pulse sequence. The echo time (TE) remains the time between the end of the 90˚ excitation pulse and the image acquisition. The application of the first gradient encodes the spins with a phase. The second gradient rephases the spins. However, if diffusion has occurred the spins will only be partially rephased and image signal will be lost. The characteristic diffusion encoding parameters are the strength of the gradient pulse (G), duration of the pulse (δ) and diffusion time (Δ) (Sotak, 2004).
Through the variation of the given parameters (TR, TE and b) and linear regression curve fitting, the relaxation parameters and apparent diffusion coefficient can be determined through NMR techniques. In order to generate an MR image, the same principles are utilized, although applied in the presence of magnetic field gradients. Gradients are used in each of the three directions to isolate individual voxels with a given resonant frequency and phase. A slice is first selected by applying a slice-selection gradient during the initial excitation in the direction corresponding to the desired slice orientation. The remaining two directions are specified by a frequency-encoding gradient and phase-encoding gradient. The phase is encoded prior to data
11 acquisition while the frequency is encoded during data acquisition. By identifying each voxel with a specific frequency and phase, an image can be generated.
2.2 Cell Membranes and Osmosis
Cell membranes provide a selectively permeable boundary which regulates exchange between the intracellular and the extracellular environment. Modeled as a lipid bilayer, as illustrated in Figure 5, cell membranes do not facilitate the bulk transport of charged nuclei, such as sodium, or cellular proteins across the membrane. The composition and balance of intracellular and extracellular ions and osmolytes is critical in regulating osmosis (Lang, et al. 1998). Osmosis is the entropic process driven by water transport from regions where the water matrix structure is highly ordered owing to small concentrations of osmolytes to regions where this structure is more disordered due to higher concentrations of osmolytes. Gradients established by these concentration differences are the driving forces of transport across the cell membrane and the basis of signal transduction in excitable cells, such as neurons.
Figure 5. The fluid mosaic model of the cell membrane. The diagram illustrates the lipid bilayer integrated with transmembrane proteins for cellular transport (Singer, et al. 1972).
Neurons in particular are sensitive to the transport of sodium, potassium, and chloride ions across the membrane. These ions are the most abundant and play the largest role in establishing osmotic gradients within the cell. For excitable cells, the ionic distributions in the cell affect the electrochemical gradient established by the neuron, which may interfere with transport across the membrane leading to cell swelling or crenation (Pasantes-Morales & Tuz,
12
2006). Although most membranes are freely permeable to water through aquaporin channels, diffusion of sodium and other ions is limited to leakage channels found within the membrane. However, in living cells, the sodium-potassium pump actively transports sodium and potassium, maintaining osmotic gradients across the membrane and making possible the active processes of life. The sodium-potassium pump uses ATP as an energy source to actively transport three sodium ions out of the cell while simultaneously transporting two potassium ions into the intracellular environment using a mechanism illustrated in Figure 6. This transport process results in a low intracellular sodium concentration and high intracellular potassium concentration. Disruption of this mechanism through blocking the sodium potassium pump by removing all potassium from the extracellular environment (or using neurotoxins) results in drastic increases in neuronal volume in neurons from the sea slug Aplysia californica (Carpenter, et al. 1992). Similar increases in cell volume were found by changing the extracellular osmolarity, providing evidence that alterations in extracellular tonicity can be used to alter cell volume and mimic cellular changes seen in pathological conditions.
Figure 6. Mechanism of the sodium-potassium ATPase. Utilizing a majority of cellular energy in neurons through ATP and a process of dephosphorylation, the Na+/K+ pump actively transports three Na+ ions out of the cell and two K+ ions into the cell to maintain osmotic gradients, regulate cell volume, and maintain or alter the cell membrane potential to act as a + + signal transducer. The two configurations of the protein for binding Na or K are denoted E1 and E2 (Faller, 2008).
13
Although cells are usually surrounded by extracellular fluid of nearly constant osmolarity, changes in the extracellular osmolarity occur in diseased cells. The effect of differing extracellular osmolarities on mouse embryonic stem cells resulted in changes in confluence as well as cellular size (Mao et al, 2008). When isotonicity is lost, the confluence of the culture is greatly reduced as many cells lose viability. In surviving cells, cell volume changes occur as hypotonic osmotic perturbations result in cell swelling and hypertonic perturbation results in crenation. Further, changes in intracellular ion or protein concentrations can osmotically draw water into the intracellular environment, resulting in cell swelling. In molluscs, prolonged exposure to hypotonic perturbation results in decreased levels of organic osmolytes to balance the intracellular osmolarity (Silva & Wright, 1994). These regulatory mechanisms can result in either cell volume increases or decreases based on the tonicity of the extracellular environment to maintain cellular homeostasis (Strange, 2004).
2.3 Relation to Pathological Conditions
In particular, early studies have shown that water and sodium concentrations increase while potassium concentrations decrease in states of tissue ischemia with a mechanism is not fully understood (Young, et al. 1987). Assessed by atomic absorption spectroscopy, these changes are immediately evident, although continue to change with time, and could only be measured in nonliving tissues. In the mammalian brain, complex osmoregulatory circuits precisely control the intracellular and extracellular osmolarity and are maintained through feedback loops sensitive to the extracellular tonicity (Bourque, 2005). Although complex, these feedback mechanisms are in place to maintain cell volume and regulate ion and protein distributions in the intracellular and extracellular environments. However, disruptions of these transport processes are seen in pathological conditions such as tissue ischemia and stroke. In particular, ischemia leads to failure of all active transport processes as ATP becomes depleted leading to failure of the sodium-potassium pump and a corresponding increase in intracellular sodium (Lipton, 1999). Current understandings of the mechanisms responsible for cell death are summarized in Figure 7. Without ATP, the electron transport chain fails and an increase in glycolysis to produce ATP increases, leading to an increase in free radicals and changes in enzyme activity. This is further exacerbated by a decrease in protein synthesis which
14 only adds to the osmotic imbalances caused by ion transport failures. Focusing on ion transport failures with ischemia, Figure 8 illustrates the changes expected within the intracellular and extracellular environment due to ischemia. The increase in intracellular sodium depolarizes the neuron which leads to intracellular calcium accumulation (Pasantes-Morales & Tuz, 2006). As the interior of the cell becomes more positively charged and osmotic pressure increases, water and chloride ions flow into the cell causing the cells to swell. These increases in cell volume alter the tissue microstructure and can lead to cell death. Additionally, neurodegeneration is also thought to change protein synthesis and metabolic processes leading to osmotic changes in neurons and ultimately to neuronal death leading to brain atrophy.
Figure 7. Mechanism of cell death during ischemia and stroke. A lack of oxygen results in a lack of aerobic respiration leading to decreased ATP and a disruption of active processes. Notably, the sodium-potassium pump fails to maintain ionic distributions resulting in changing ionic distributions. Additional changes in protein synthesis, electron transport, and phospholipase activity disrupt many cellular physiological processes (Lipton, 1999).
15
Figure 8. Diagram of the mechanism of osmotic redistributions that result due to ischemia and stroke. Failure of the sodium-potassium pump results in increased intracellular sodium and depolarization of the cell leading to increased intracellular calcium. To preserve charge balance, chloride ions leak into the cell along with water resulting in cell swelling (Pasantes-Morales & Tuz, 2006).
Magnetic resonance imaging has been utilized to identify osmoregulation changes that occur in stroke (Baird & Warach, 1998; Merino & Warach, 2010). First studied in cats, MR images showed that ischemic stroke could be identified as early as six and a half hours after onset and had clinical implications (Buonanno, et al. 1982). Relaxation measurements, specifically increases in T2 which resulted in an increase in image signal intensity, were used as the basis for identifying ischemic stroke. However, the first diffusion-weighted MR images were shown to provide evidence of stroke at much earlier time points than conventional T2-weighted imaging in cats (Moseley, et al. 1990). One hour after the onset of ischemia, T2-weighted MR images appear normal, but significant increases in the signal intensity in diffusion-weighted MRI due to a reduced ADC were observed, indicating that diffusion-weighted imaging was more sensitive in identifying stroke. Through the combination of relaxation, diffusion, and perfusion MR measurements, tissue ischemia and stroke can be identified in rats with mismatches between the diffusion, perfusion, and relaxation infarct size (Hoehn-Berlage, et al. 1995; Calamante, et al. 1999; Beneviste, et al. 1992; van Dorsten, et al. 2002). The region of damage is underestimated
16 by relaxation measurements at early time points while perfusion imaging does not necessarily indicate salvageable tissue. Although first performed on a variety of animal models (Hoehn, et al. 2001), MR imaging is noninvasive and through relaxation and diffusion measurements, the onset time of stroke in the clinical environment can be estimated to allow better patient treatment (Petkova, et al. 2010; Jokivarsi, et al. 2010). Additionally, both cancerous tumors and neurodegeneration also exhibit drastic changes in cellular osmoregulation and cell viability. Similar MR relaxation and diffusion measurements are useful in identifying brain tumors (Naruse, et al. 1987) and identifying tumor changes with treatment (Chenevert, et al. 1997). In mouse models of Alzheimer’s disease, reductions in T2 were observed in several brain regions (Helpern, et al. 2004). Further, T2*-weighted MRI has been used to identify plaques in human tissue associated with Alzheimer’s disease (Beneviste, et al. 1999) and standard MR imaging protocols can detect reductions in hippocampal size that are indicative of neurodegeneration. Although the different pathological conditions can be identified, the cellular changes that result in different contrast in MR images are not fully understood. MRI typically looks at water in the human body and therefore changes in cellular water distributions are likely responsible; however, water distributions are heterogeneous and water diffusion is highly restricted leading to difficulties in understanding MR relaxation and diffusion measurements. However, MR diffusion decreases have been correlated with reductions in activity of the sodium-potassium pump resulting in increased intracellular sodium, decreased intracellular potassium and cell swelling (Mintorovitch, et al. 1992). Similar diffusion decreases also correlate with decreased extracellular volume fraction (increased intracellular volume fraction due to cell swelling) and increases in extracellular tortuosity (van der Toorn, et al. 1996).
2.4 Single Cell Models
As MRI hardware continues to improve and higher magnetic field strengths are achieved, higher resolution images become possible due to an improved signal-to-noise ratio. With resolutions below 100 µm, MR microscopy has been developed using small RF microcoils along with small samples to look at individual isolated cells (Ciobanu & Pennington, 2004). Initial efforts were put forth using the toad Xenopus laevis oocyte due to its large size (~1 mm) and unique regions including the nucleus, animal pole, and vegetal pole which could be clearly
17 differentiated (Aguayo, et al. 1986; Mallard, 1986). Increased signal intensity was observed in the animal pole indicating a greater concentration of free water which showed that physiological differences existed between different regions of the cell. Further, using water suppression techniques, it was shown that a majority of the lipids are isolated to the cytoplasm and the nucleus contains mostly mobile water (Posse & Aue, 1989). Since the cell is not differentiated, there are significant lipid fractions within the oocyte suitable for MR imaging (Sehy, et al. 2001) and spectroscopy (Lee, et al. 2006). Expanding on these findings, MR relaxation measurements
(T1 and T2) identified relaxation differences between these three regions. The nucleus had the longest T1 and T2 relaxation parameters at a resonant frequency of 300 MHz (T1: 1856 ± 100 ms;
T2: 29.45 ± 6.95 ms) followed by the animal pole (T1: 1207 ± 173 ms; T2: 18.81 ± 4.08 ms) and the vegetal pole (T1: 950 ± 191 ms; T2: 11.07 ± 3.05 ms) supporting the finding that the nucleus has a high proton spin density due to a large pool of free water (Pauser, et al. 1995). Relaxometry methods were used to differentiate the regions of the cell and identify the contributions of free and bound water to the MR signal. The identification of two distinct water pools was in agreement with NMR spectroscopy methods that were correlated with destructive methods to determine the water weight in entire oocytes and subcellular regions in response to osmotic perturbation using sucrose solutions (Cameron, et al. 1990). It was reported that no intracellular water pool was equivalent to bulk water or water in dilute solutions and that there were hydration layers in which the water molecules were bound and not in fast exchange. Further, the bound and free water fractions are dependent on the cell type and change in relaxation in protein concentration and other factors that may hinder water exchange. Diffusion measurements on the Xenopus oocyte were expanded to include osmotic perturbations of the extracellular environment and reveal a linear increase in ADC with increasing cell volume (Sehy, et al. 2004). Although these images were obtained on an oocyte and not a fully differentiated cell, it sparked interest in understanding the cell regulatory mechanisms that contribute to the volume- averaged tissue response seen in clinical MRI. In order to examine the effect of osmosis and ionic sodium distributions on a single neuron, many factors must be controlled in the experimental design. A common concern in large organisms is the heterogeneity of tissues and an uncontrolled cellular environment. To combat such concerns, a single isolated neuron was used for single cell studies in a controlled external environment (Grant, et al. 2000). An ideal animal model for such experimentation is the sea slug Aplysia californica which contains a large
18
L7 neuron capable of being isolated from the organism and remaining viable for extended periods of time after dissection (Schoeniger, et al. 1994; Bowtell, et al. 1995). With a diameter of 300-400 µm, the L7 neuron is much larger than typical cells and can be observed and manipulated using standard microscopic techniques. Also, neuronal cells are excitable in which sodium plays an important role in the propagation of electrical signals along the cell. Cell swelling of neurons is also of concern as damage to such cells is typically irreversible making an understanding of cell swelling due to osmosis of particular interest. In studying single neurons isolated from Aplysia californica, clear differences appear to exist between the nuclear and cytoplasmic compartments within the neuron. The nucleus and cytoplasm have a clear chemical composition of osmolytes including betaine and choline that are maintained in relatively high concentrations (hundreds of mmols) that are critical in maintaining osmotic balances in the neuron (Grant, et al. 2000). By imaging over extended time periods in upwards of eight hours, the detected signal of these osmolytes significantly decreased in both the nuclear and cytoplasmic compartments showing a loss of cell membrane integrity which leads to cell swelling and eventually cell death. The nucleus has longer T2 relaxation, resulting in an increase in signal intensity. Additionally, the ADC can be measured in both compartments (Hsu, et al. 1996), showing differences in diffusion between these intracellular regions where the biexponential nature of water diffusion is evident (Grant, et al. 2001). Since the membrane is semipermeable, it was thought that diffusion anisotropy may exist within cells as diffusion may become more prevalent in a single direction. However, upon encoding the ADC in two perpendicular directions, it was shown that tissue anisotropy is unlikely in cells as there are no clear differences in diffusion in the different direction (Hsu, et al. 1997). These results show that diffusion is likely isotropic and maintains this characteristic under osmotic perturbation. Under osmotic perturbation, diffusion in the nuclear and cytoplasmic compartments also change in different ways and can be attributed to differences in the permeability of the membranes as the nuclear membrane allows diffusion more readily. The T2 times increase under hypotonic perturbation, while the diffusion coefficient decreases in the nucleus. However, the ADC measured in the nucleus and cytoplasm remains relatively constant under osmotic perturbation (Hsu, et al. 1996). The response of osmotic perturbation of human brain slices illustrates this similar trend in an increased fast-diffusing fraction within the tissue as the tonicity is increased (Shepherd, et al. 2003). Additionally, diffusion in these human brain slices follow a
19 biexponential model, similar to that observed in the single cell analysis showing that the cellular response can be utilized in understanding how tissues respond to pathological insults and that animal models provide similar tissue environments to those found in humans.
2.5 Tissue Models
Osmotic, neurotoxic, and ischemic external perturbations of the extracellular environment interfere with cellular transport processes, resulting in cell volume and osmoregulation imbalances. Methods to quantify diffusion in these pathological states have led to the incorporation of magnetic resonance imaging (MRI), specifically diffusion-weighted imaging, as a tool to evaluate changing cellular microstructure and water distributions within cellular structures, as can be seen in an isolated neural tissue model. Using the rat hippocampus (van Pul, et al. 2005) and the turtle cerebellum (O'Shea, et al. 2000) as model tissue structures, diffusion, in the context of cell volume regulation under osmotic perturbation and oxygen and glucose deprivation, can be quantified. The apparent diffusion coefficient (ADC) is an indicator of molecular movement within a sample and can be determined using diffusion-weighted MRI. Of particular interest, the ADC decreases in states of tissue ischemia, particularly in acute stroke, with a biophysical mechanism that is not understood. O’Shea, et al. focused on the effect of osmotic perturbations on the intracellular and extracellular volume as assessed by the ADC and the MRI relaxation parameter,
T2. It was hypothesized that the ADC relates to both the intracellular and extracellular volume and that changing ADC values under osmotic perturbation leads to changing volume fractions. Similarly, van Pul, et al. tested the idea that osmotic perturbation and the resulting cell swelling only partially contributed to the decreased ADC seen in ischemic tissue, hypothesizing that oxygen and glucose deprivation also plays a role in the reduced ADC. Although, cell swelling leads to a decreased extracellular volume due to increased tortuosity, the ratio of restricted diffusion to diffusion in free water, changes in the intracellular environment and the tissue itself can lead to a reduced ADC. By separately investigating the effects due to osmotically-induced cell swelling and those due to ischemia, a new aspect of the cellular response can be quantified. As with any living in vitro tissue model, there are factors that cannot be controlled. Cellular
20 metabolism, as well as inherent repair and cell volume regulatory mechanisms, is not considered providing limitations to the applicability to in vivo model systems. To quantify these changes in the ADC, dissected tissues were perfused to maintain tissue viability during MRI analysis. Both the turtle cerebellum (O'Shea, et al. 2000) and rat cerebellum (van Pul, et al. 2005) were chosen because the characteristic intracellular and extracellular volume fractions are well defined. Additionally, comparisons of different tissue models validate the results as being applicable to different tissue types. The turtle cerebellum was held between a cotton cloth perfused with an oxygenated isotonic control solution and osmotic perturbations were introduced by changing the sodium chloride concentration. However, changing only this concentration could lead to neuronal activation which would alter the cellular response to perturbation. The rat hippocampus was held with a nylon mesh perfused with oxygenated artificial cerebrospinal fluid, however, the mesh could provide limitations to cell swelling. Osmotic perturbations were introduced by changing the amount of water added to solutions, which maintained the relative concentration of osmolytes; glucose and oxygen deprivation was introduced simply by removing glucose and bubbling the solution with nitrogen instead of oxygen. Similar diffusion-weighted MR images were acquired. The general conclusion of both studies shows that under hypertonic perturbation the ADC increased while under hypotonic perturbation the ADC decreased showing a nearly linear relationship between the ADC and the extracellular osmolarity. The average ADC was similar in both studies at 6.8 x 10-4 mm2/s compared to 4.69 x 10-4 mm2/s in the rat hippocampus and turtle cerebellum, respectively. Additionally, O’Shea et al. concluded that the T2 value decreased under hypertonic perturbation and increased under hypotonic perturbation. The diffusion coefficients are lower in the intracellular space by a factor of nearly 2, the response is much faster under hypertonic perturbation (time constant of 5 min for hypertonic and 11 min for hypotonic), and the ADC changes faster and more significantly than T2 (32-50% increase over baseline compared to 8-22% for T2). The results show comparable trends to both a fast and slow exchange model, which accounts for independent diffusion in the intracellular and extracellular compartments, showing the physical system likely obeys an intermediate combination of the two exchange models. However, when accounting for tortuosity, viscosity, and molecular interactions, these results show that the total tissue water does impact T2 and the ADC varies with the extracellular space fraction.
21
Similarly, van Pul, et al. showed little anisotropy of the diffusion within the tissue, as there is limited myelin enclosed structures and restrictive membranes. The ADC decreased by 8% under no perturbation, 83% under hypotonic osmotic perturbation, and 57% under glucose and oxygen deprivation. The decrease of the ADC under no perturbation shows there are additional cellular mechanisms accounting for ADC decrease that have not been accounted for. Additionally, oxygen deprivation is critical for the time scale of ADC changes as changes occurred much more rapidly when oxygen is removed, on the order of 3 hours to reach the final ADC deviation in comparison to 8 hours when oxygen is present. However, changes in the extracellular volume fraction are similar between perturbations. Of particular interest, osmotic perturbation resulted in little cell damage, while oxygen and glucose deprivation resulted in detrimental damage to the tissue. It follows that the ADC changes do not result exclusively from cell swelling as the intracellular ADC may be changing due to ischemic effects not related to osmotic perturbation. The ADC can be related to the intracellular and extracellular volume fractions. However, reductions in ADC are due not only to osmotically induced cell volume manipulations, but also to ischemia, which will have large impacts on future work related to the onset of ischemia and stroke in clinical situations. A method for imaging nervous tissue slices has been well developed for MR analysis (Blackband, et al. 1997). The applicability of these tissue model systems shows that both rat and human hippocampal slices exhibit similar MR relaxation and diffusion properties.
2.6 Diffusion Multi-Compartment Models
Together with other MRI studies, a multi-component diffusion model is used to quantitatively assess parameters that affect diffusion across the cell membrane through the use of diffusion weighted MRI. The model for diffusion in isolated neurons relates to the compartmentalization of water in the cell resulting in the existence of both bound and free water with slow transport properties between the intracellular and extracellular environments. A multi- compartment model of diffusion can adequately represent a tissue model of diffusion (Stanisz, 2003. Using human erythrocyte ghosts, water diffusion was found to be non-monoexponential with significant exchange between the intracellular and extracellular compartments (Thelwall, et al. 2002). This method can be used to observe how a cell responds under different osmotic
22 perturbations, the corresponding cell volume increase, and the rate of cellular osmosis. Cell membranes aid in the compartmentalization of the cell; however, the properties of the individual compartments within the cell must be known to be able to understand the heterogeneous system. A diffusion model must account for water exchange between the nucleus, the cytoplasm, and the extracellular environment, but the slow diffusion of such systems ensures that a large fraction of the water remains in the original compartment. Of concern in the diffusion weighted MRI studies is the SNR dependence on coil and voxel size becoming noticeably apparent in imaging when the resolution must be increased. In an effort to increase the resolution of such imaging sequences, the SNR must be increased significantly, often resulting in a tradeoff between SNR and resolution. The measurement of the diffusion coefficients under different osmotic perturbations, including both hypotonic and hypertonic solutions, will be a basis of understanding the mechanism of osmosis on cell volume regulation. Cellular modeling of diffusion becomes complicated by the fact that no tissue sample is homogeneous as cells, extracellular proteins, as well as the extracellular fluid can impact the rate of diffusion. Recently, efforts have been made to determine the parameters that most significantly affect the changes in ADC measured under states of perturbation by modeling the tissue as a collection of cuboidal cells. The ADC decreases by up to 50% following ischemia in tissues which can be accounted for by cell swelling under hypotonic perturbation (Harkins, et al. 2009). The intracellular diffusivity and heterogeneity of tissue structures decrease the ADC in tissues due to the presence of impermeable membranes. The intracellular diffusivity and membrane permeability most greatly effect diffusion as decreases in ADC can be accounted for by decreases in membrane permeability and diffusivity in the intracellular environment. The tortuous intracellular environment causes a faster signal decay leading to reduced ADC measurements, showing that the compartmentalization of tissues is required to accurately account for the measured ADC (Harkins, et al. 2009). It has also been shown that intracellular hydration increases under hypotonic perturbation and the diffusion coefficient decreases in states of increased protein concentration which can directly account for reductions in T2 relaxation values and ADC measurements (Kotek, et al. 2009). Changes in tissue microstructure directly impact both the intracellular and extracellular spaces and consequently alter diffusion (Norris, 2001).
23
2.7 Sodium MRI
Although sodium is much more difficult to image due to its quadrupolar nature, weak signal emission, and low cellular concentration, the use of high magnetic fields allows the imaging of sodium in tissues (Boada, et al. 2005). A quadrupolar nucleus is defined as having a quantum spin number greater than 1/2 (noting that the traditional spin 1/2 nucleus is 1H). Sodium (23Na) is a spin 3/2 nucleus and therefore has a quadrupole moment due to a higher degree of nuclear asymmetry. As opposed to spin 1/2 nuclei that split into two unique energy levels when an external magnetic field is applied, spin 3/2 nuclei split into four levels characterized by the four magnetic quantum numbers +3/2, +1/2, -1/2, and -3/2, representing four possible orientations in a static magnetic field. Therefore, the spin 3/2 nucleus has more than one allowed single quantum transition (3/2 → 1/2, 1/2 → -1/2, -1/2 → -3/2). When comparing sodium and proton MRI, it must be noted that sodium is not homogeneously distributed in tissues, but rather has an intracellular concentration of sodium of 10 mM and an extracellular concentration of 150 mM in humans. Due to poor spatial resolution, the sodium signal is volume averaged over both environments with a total sodium concentration near 50 mM in the brain, due to the large intracellular volume fraction. In comparison to proton imaging where the water concentration is frequently above 50 M, this presents a significant difference and the total sodium concentration limits the signal obtained from a given voxel. Further, several other important differences exist. 23Na is found in 100% natural abundance, similar to the near 100% natural abundance of 1H. However, protons have a gyromagnetic ratio of 42.58 MHz/T while sodium has a gyromagnetic ratio of 11.26 MHz/T, giving sodium a resonant frequency of 132 MHz at 11.75 T while protons have a resonant frequency of 500 MHz at 11.75 T. These differences give sodium a relative sensitivity of approximately 10% that of 1H. Although the sensitivity of 23Na is an order of magnitude below 1 H, its relaxation rates are two orders of magnitude faster with T1 ~ 30-40 ms and T2 ~ 30-40 ms. Factoring the low sodium concentration and reduced sensitivity, the total sensitivity of 23Na-MRI is 10,000 times lower than 1H-MRI, which limits spatial resolution and requires higher magnetic fields, small radiofrequency (RF) coils, and extensive averaging leading to longer imaging times. Due to these factors, the resolution of sodium is poor with the maximum resolution on the order of 100 μm, although maximum clinical resolutions are on the order of 1-2 mm.
24
Given the quadrupolar nature of the sodium nucleus with multiple quantum transitions, the signal obtained from each of these transitions is not uniform, but rather the relative intensities are 3:4:3. Static quadrupolar effects result due to the electric field gradient at the nucleus, while dynamic quadrupolar effects occur due to fluctuations in the electric field at the nucleus. Both longitudinal and transverse (T1 and T2) intrinsic relaxation mechanisms are affected by these interactions. Static quadrupolar effects can cause the two outer transitions (3/2 → 1/2 and -1/2 → -3/2) to experience different resonant frequencies and will be broader, distributed over the entirety of the spectral range (deGraff, 2007). In tissue samples heterogeneous broadening occurs when the tissue is not macroscopically oriented and the outer transitions become blurred together over the entire spectral frequency range. Dynamic quadrupolar effects occur due to fluctuations in the electric field at the nucleus. The extreme narrowing condition is related to the correlation time (ω0τc « 1) and is related to the resonant frequency and the rotational correlation time. If this condition is not fulfilled, biexponential behavior becomes evident. Further, the T1 and T2 relaxations of a quadrupolar nucleus are biexponential in nature. The first in vivo sodium images were produced in 1983 of a cat experiencing acute cerebral stroke (Hilal, et al. 1983). It became evident that sodium MRI was extremely sensitive to pathological conditions as one of the first changes in diseased tissue is the loss of ion homeostasis. By 1985, clinical sodium MRI of humans was feasible (Perman, et al. 1986). Although initial studies examined primarily the sodium signal changes without quantifying either the sodium relaxation parameters or the total sodium concentration, methods have been developed to address these concerns. It was observed by Christensen et al. that the signal intensity for high SNR images is expected to be linear with concentration, given that T2 weighting is minimized and the B0 and B1 fields are homogenous over the sample volume. Field homogeneity concerns can be minimized by shimming on the proton water signal using a double- tuned RF coil. The relationship between sodium signal intensity, I, and total sodium concentration, TSC, can be described in equation (13) when two sodium reference phantoms
(concentrations of C1 and C2 with corresponding image intensities of I1 and I2) are included in the imaging field of view.
25
(13)
It was shown that the total tissue sodium concentration in the rat brain using sodium MRI (45 ± 4 mM) was in excellent agreement with the 22Na radionuclide dilution technique (49 ± 6 mM), however the obvious advantage is that the sodium MRI can be used in vivo and does not require radiation (Christensen, et al. 1996). However, sodium distributions within ischemic tissue and in cancerous tumors have been shown to increase due to the perturbed states (Lin, et al. 2001; Thulburn, et al. 1999). Thulburn et al. focused on the use of sodium MRI to assess changing sodium distributions in tumors in the rat cerebellum with and without interventional treatments. MRI analysis showed a significant increase in sodium concentration in the tumor region as illustrated by an increase in signal intensity. The affected area correlated with morphological changes as seen in the proton image and with ex vivo histological staining confirming that sodium MRI was capable of identifying the affected region and monitoring changes in sodium distributions as a response to pathological conditions. The sodium concentration increased from around 40 mM to 65 mM, a significant increase, possibly illustrating that sodium has a key role in identifying diseased tissues. Although the specific treatments performed in the study did not alter the sodium concentration, this provides an initial starting point to assess therapeutic treatments for cancer. Sodium imaging was expanded by Lin et al. in comparing the sodium response to 1 ischemic conditions. Similar reductions in T2 and increase in ADC were reported for the H imaging studies which correlated with increases in 23Na signal intensity on the same animal subjects. The sodium signal intensity increased in the ischemic regions, possibly due to an increased spin density related to higher sodium concentrations. The exact relaxation mechanism, either T1 or T2, was undetermined and future studies must focus on relating these changes to the physiological conditions in perturbed tissue (Thulburn, et al. 1999). However, the sodium distributions change more rapidly and remained in an altered state for a longer period of time
26 than the corresponding proton images. Even at low resolution, quantitative analysis of diseased nervous tissue show changes in sodium distributions. Sodium imaging could appear normal, even after ischemia has set in; however, there is potential to observe the response to sodium in neurological conditions as a possible method for the early detection, prior to cell death, of neurological diseases. 23Na imaging has become more prominent due its promise of earlier identification of altered cell physiology in vivo. Similar to all MR imaging, 23Na MRI does not affect cognitive function or impose any serious health risks to patients (Atkinson, et al. 2010). Increases in 23Na signal intensity are seen in ischemia (Thulborn, et al. 2005) and tumors (Schepkin, et al. 2011). Although the increase in signal intensity may be due to increases in the total sodium concentration, it could also result from changes in the 23Na relaxation parameters. To provide 23 23 evidence of changing Na relaxation, the Na T2* relaxation parameter was mapped in the human brain and differences were seen between grey and white matter, as well as a much higher sodium content in cerebrospinal fluid (Fleysher, et al. 2009). Recently, the field of sodium MRI has expanded to look at sodium alterations that occur as a result of neurodegeneration in clinical settings. Increases in 23Na signal intensity were shown to correlate with Alzheimer’s disease (Mellon, et al. 2009) and multiple sclerosis (Inglese, et al. 2010); however contributions from 23Na relaxation were not included.
27
CHAPTER THREE
OBJECTIVES
The specific objectives of this work are based on characterizing the cellular response to osmotic perturbation, which results in cell swelling or crenation, through the MR relaxation and diffusion properties of the abdominal Aplysia californica ganglia. Both 1H and 23Na will be utilized to examine the benefit of using non-proton nuclei in assessing cell viability and cellular responses to insults of the extracellular environment. Changes in the MR relaxation parameters and diffusion properties can be used as direct indicators of change in the cell structure and water distribution involved in cell volume regulation. Therefore, MR analysis can show the impact altered extracellular tonicity has on cellular ion distributions as well as morphological change and cellular damage within the ganglia providing insight to the osmotic regulation of cell volume in isolated neural ganglia using a non-destructive technique. It is hypothesized that 23Na MRI at high magnetic field strengths will show that redistributions of ionic sodium occur prior to morphological changes in the cell or apoptotic events and are more sensitive to the physiological environment than conventional 1H MR imaging techniques. Although there are many aspects to cell volume regulation and redistributions of ionic sodium within the cell, this study will focus on altered cellular transport mechanisms due to osmotic perturbations, allowing quantitative characterizations of cell volume regulation in the model Aplysia abdominal ganglia. Specifically, 1 23 differences in H and Na magnetic resonance imaging contrast mechanisms (T1, T2, T2*, and ADC) in fixed, viable, and nonviable neural ganglia will be assessed. The project will be aimed at examining these differences in the context of understanding diseased cell states in pathological conditions.
3.1 High-Field MRI
Magnetic resonance imaging (MRI) will be used in the analysis of ion and water distributions in the isolated neural ganglia as well as in the measurement in the rate of osmosis, characterized by the apparent diffusion coefficient. MRI relies on magnetism to image atomic nuclei and both 1H and 23Na nuclei are detectable through magnetic resonance principles.
28
Although sodium is much more difficult to image, it nevertheless has critical functions in regulating cell volume and initiating both membrane depolarization and signal transduction in 1 excitable cells. Isolated Aplysia ganglia were characterized through H T1, T2, T2* and ADC 23 measurements, as well as Na T1 and T2* measurements. Changes in the relaxation and diffusion measurements were quantified and related to the extracellular osmolarity and viability of the ganglia. By imaging both 1H and 23Na, the role of sodium concentration gradients in the regulation of osmosis and the role of active transport processes on the maintenance of ionic sodium distributions in isolated neural ganglia becomes evident through examination of the MRI relaxation parameters. The MRI relaxation parameters (T1, T2 and T2*) and the apparent diffusion coefficient (ADC) from diffusion-weighted magnetic resonance imaging were used to quantitatively assess changes in tissue microstructure within the ganglia, allowing the effect of osmotic perturbations on cellular osmoregulation to be observed. Although averaged over the entirety of the ganglia, the ganglia is compartmentalized and the measurements represent the average tissue response of the nuclear and cytoplasmic regions of the neurons, as well as the interstitial environment. By quantifying changes in diffusion which is utilized as an indicator of the relative rate of osmosis, cell volume regulation through altered water distributions as well as changes in tissue microstructure due to altered barriers to diffusion are observed. The relaxation and diffusion measurements are aimed at examining changes in the context of understanding and reversing pathological insults in chronic neurodegenerative disease and tissue ischemia. It is hypothesized that 23Na MRI will permit the assessment of the osmotic regulation of cell volume in a prototypical isolated neural ganglia and that redistributions of ionic sodium will occur prior to morphological changes in the cell or apoptotic events. MR relaxation and diffusion measurements of isolated neurons from the Aplysia abdominal ganglia indicate that hypotonic perturbations result in increases in T2 and no apparent changes in ADC which serves as the foundation for work on the abdominal Aplysia ganglia (Hsu, et al. 1996). Both hypertonic and hypotonic perturbations of the extracellular environment were introduced and compared to isotonic conditions to illustrate the role changing ionic gradients have on the detection of cell swelling through MR relaxation and diffusion mechanisms.
29
3.2 Radiofrequency (RF) Coil Design
A radio frequency (RF) coil must be utilized in order to transmit RF pulses to the system and receive the MR signal. A solenoidal coil achieves the highest signal-to-noise ratio (SNR) of available coil configurations and therefore necessary to achieve the required signal to observe the cellular structure of neural ganglia samples (Glover & Mansfield, 2002). However, coil development continues to improve and a large number of coil configurations are possible (Neuberger & Webb, 2009). A tradeoff exists between spatial and temporal resolution as the resolution is directly related to the imaging time and signal obtained. For example, in imaging brain slices a resolution of 120x230x300 μm can be obtained in 4 minutes while a resolution of 15x15x300 μm necessitates a 14 hr image acquisition (Benveniste & Blackband, 2002). The resolution, however, can be improved without sacrificing imaging time through the improvement of the SNR, defined as the difference in mean signal and mean noise obtained divided by the standard deviation of the noise. A higher SNR corresponds to higher signal intensity in the acquired MR images (Minard & Wind, 2001). The SNR can be improved by utilizing higher magnetic field strengths as there is a strong dependence of SNR on the externally applied magnetic field strength, Bo, as shown in equation (14) (Peck, et al. 1995).
(14)
Further, the SNR can be improved through the use of small solenoidal coils. For a solenoidal coil of a specified length (L) and diameter (D), the SNR improves according to equation (15). It can be seen that the length to diameter ratio is important is maximizing the signal obtained (Peck, et al. 1995).
(15)
√ ⁄
√ ( ⁄ ) A radio frequency (RF) coil was constructed which could be used at multiple frequencies. The solenoidal coil was double tuned to both the resonant frequencies of 1H and 23Na. At field strength of 11.75 T, this corresponds to a coil tuned to 132 MHz and 500 MHz for 23Na and 1H, 30 respectively. The large range of required frequencies necessitates the use of a versatile coil which can be easily manipulated to tune to varying frequencies. Two separate channels were utilized on the coil and each was tuned to the operating frequency at the given magnetic field strength, which fluctuates slightly with time. Changing between the resonant frequency of the different nuclei required changing the tuning capacitor on the coil to a different capacitance, C, as the relation between the inductive and capacitive properties of the coil are shown in equation (16) where ω is the resonant frequency, L is the inductance, and C is the capacitance (Minard & Wind, 2001).
(16)
The coil was constructed by wrapping copper wire around a fused silica capillary. A capillary with a 3 mm diameter was used to facilitate the larger ganglia samples, although small microcoils (700 µm diameter) have been used to image single neurons (Grant, et al. 2001). SNR analysis will be useful in analyzing coil performance at varying field strengths (Peck, et al. 1995). Due to its low sensitivity, sodium is inherently difficult to image and requires the use of high magnetic field to achieve signal to noise ratios that yield signal strengths that can be clearly identified from the noise background. The SNR depends on the size of the RF coil, the size of the voxel which is being imaged, the static magnetic field strength, the spectral width and acquisition time (Webb, 2003). In order to increase the SNR, small coils must be used at high magnetic fields and the RF coil design is at the millimeter length scale to increase SNR (Minard & Wind, 2001). With an inherently small voxel size due to the microscopic nature of the neurons within the ganglia, these efforts must be taken to reduce the total acquisition time. This time reduction will reduce the number of averages in an effort to limit the imaging time to levels short enough to ensure cell viability. During imaging, the neural ganglia were placed in a capillary tube and centered in the RF solenoidal coil. Due to the coil configuration, 1H and 23Na imaging were performed in separate studies which did not allow for simultaneous data acquisition of both the 1H and 23Na signals. The experimental setup was simplified by not perfusing the ganglia (viability was expected to be maintained for up to seven hours after dissection) which allowed the ganglia samples to be directly inserted into the coil.
31
3.3 Neural Ganglia Model System
In understanding the role of the isolated neuron in the context of larger tissues and cellular interactions, isolated abdominal ganglia from the sea slug Aplysia californica were quantitatively described using both 1H and 23Na magnetic resonance imaging. The abdominal ganglion houses the L7 motor neuron which has been used extensively in single cell MRI studies. Expanding to the imaging of the abdominal ganglia allows comparisons of the entire ganglia to previous work using isolated neurons; cell to cell interactions and tissue homogeneity issues in the MR relaxation and diffusion properties of single neurons within a tissue environment can be explored. Further, the abdominal ganglion has been extensively characterized for cellular electrophysiology studies. It is a suitable model system due to its relatively simple anatomy, consisting of a small number of large neurons which can reach sizes more than 300 µm in diameter (Coggeshall, 1967). The general anatomy of the ganglion consists of two spherical hemispheres connected by a commissure and the entire ganglion is encapsulated in a fibrous connective tissue sheath (Frazier, et al. 1967). The sheath is not uniform, but is thickest where the nerves attach to the abdominal ganglia and thinner in other areas. Also, the ganglion is vascularized through the main dorsal aorta which branches into several arteries; however, the arteries end in the sheath surrounding the ganglia and do not penetrate into the interior of the ganglia, but rather empty the blood directly into the sheath indicating that the nervous tissue itself (consisting of the nerve and glial cells within the ganglion) are avascular and rely on diffusion for nourishment (Coggeshall, 1967). The connections of the large neurons within the abdominal ganglia have also been characterized (Kandel, et al. 1967). There is a 30-40% increase in the number of axons when large animals (200 g) are compared with juvenile animals (5 g). Further, increases in the size of Aplysia neurons can be related to temperature (Treistman & Grant, 1993) and osmoregulatory processes become slower as the animals age (Skinner &Peretz, 1989) making it important to account for factors that may alter the physiology of the ganglia. The neurons themselves are located along the periphery of the ganglion surrounding a central region containing neutrophil containing the presynaptic and postsynaptic processes of the neurons where synapses between neurons occur. Although the number of large neurons remain the same throughout development, as the animals age the total neurons increase by up to 40% mainly through additional small neurons at the boundary of the interior neutrophil. The large neurons are unipolar cells with a large nucleus that accounts for
32
30% of the cellular volume but as the cells become larger the nuclei become flatter and occupy less space. The abdominal ganglion contains between 1,000-2,000 cells and 15,000-25,000 axons making it a very simple nervous system model in comparison to higher organisms (Coggeshall, 1967). As the animals mature, several key changes occur within the neurons of the abdominal ganglion including an increase in the size of the nucleus and cell body with added irregularities in shape and granule distribution. The changes of developing neurons are evident in MR relaxation and diffusion measurements (Aiken, et al. 1996). In imaging Aplysia abdominal ganglia, temperature, age, and the general size and appearance of the ganglia were kept uniform to prevent unwanted changes in cellular physiology.
3.4 Osmotic Perturbation
By examining the cellular response to different osmotic insults of the extracellular environment, much can be learned about cell volume regulation and the influence of changing ionic concentration distributions on cell swelling and osmosis in diseased states. In states of reduced tonicity, there is a net influx of water into the cell and the cell swells, while in states of increased tonicity the cell shrinks. The morphological changes under state of osmotic stress lead to disordered cell shapes which impacts transport in cells and tissues (Chen, et al. 2000). Of the cellular ions, sodium is paramount in cell volume regulation due to the large variations of concentration in the intracellular and extracellular environments. Sodium ions are polar compounds incapable of diffusing through the highly non-polar interior region of the cell membrane in significant quantities. Relying on active transport processes, notably the sodium- potassium pump, sodium transport is directly controlled by each cell and is a regulatory mechanism which accounts for a great majority of cellular energy consumption (Lang, et al. 1998). The sodium-potassium pump is a key regulator of sodium transport across the cell 23 membrane and therefore in controlling osmosis. By quantifying the Na T1 and T2* relaxation parameters and the 23Na signal intensity, the sodium redistributions due to osmotic insult can be characterized. Through this understanding, further efforts may be undertaken to combat the effects of such insults which may be used in reversing the onset of cell death in neurological disease.
33
Osmosis is driven by osmotic pressure differences which can be determined using equation (17) where i is the dissociation factor, M is the concentration, R is the gas constant, T is the absolute temperature, and π is the osmotic pressure (Strange, 2004).
(17)
Following isotonic studies on neural tissue, the osmolarity of the perfusate will be varied by altering the ionic concentrations found in the extracellular medium. Similar MRI analysis will be performed and the parameters for the ganglia will be compared to the isotonic values initially obtained. In so doing, the effect of certain perturbations can be related to decreases in the intracellular volume fraction under states of hypertonic perturbation and increases under hypotonic perturbation. Therefore, it can be concluded that changes in the relaxation parameters are direct indicators of changes in the cell structure and water distribution involved in cell volume regulation (Hsu, et al. 1996). Although there are many aspects to cell volume regulation and redistributions of ionic sodium within the cell, studies have used MRI to examine the transport mechanisms within a cell under both isotonic conditions and with osmotic perturbations (Hsu, et al. 1996). Allowing quantitative characterizations of transport mechanisms to be developed in the model Aplysia californica isolated neuron, osmotic effects have shown to decrease the T2 relaxation time under hypotonic perturbation as is commonly encountered in tissue ischemia. Through the efforts of examining the effect of osmotic perturbation on isolated neural ganglia, the mechanism of cell volume regulation in neurons continues to be better developed to understand diseased cell states in neurodegeneration.
3.4 Cell Viability
Of concern in many single cell and tissue model systems is cell viability and how significantly physiological mechanisms are altered due to cell death. As cells undergo apoptosis, significant changes in water distributions alter MR relaxation and diffusion measurements (Kettunen & Brindle, 2005). By quantifying the proton relaxation and diffusion properties of the ganglia immediately after dissection and again after cell viability is lost, the changes associated with loss of viability can be identified. Sodium redistributions may present novel ways of
34 assessing cell viability since as cells become nonviable, active transport ceases and membrane permeability is altered resulting in significant ionic redistributions. Although the sodium relaxation parameters are only measured on nonviable ganglia, by observing the sodium signal intensity over time, insight into the changes in ionic distributions within the ganglia can be identified.
35
CHAPTER FOUR
METHODS
4.1 Neural Ganglia Preparation
Aplysia californica were acquired from the National Resource for Aplysia at the Rosenstiel School of Marine and Atmospheric Science at the University of Miami and were in the late juvenile stage of development with a weight of 46 – 75 g. The Aplysia were shipped in plastic bags filled with oxygen and sea water. Animals were maintained in a seawater aquarium at 18-20 ˚C, slightly higher than their natural water temperature, and approximate pH of 7.8-8.0. Aquarium salinity was adjusted by adding Instant Ocean sea salt to achieve a specific gravity of 1.020-1.025. The Aplysia could be maintained in this environment for several weeks by checking the salinity and feeding every three days with washed romaine lettuce. The abdominal, or visceral, neuronal ganglion from each Aplysia was obtained by first anesthetizing the slug with a magnesium chloride paralytic agent injected in three separate locations along the muscular foot of the animal. Approximately 15-20 mL of the anesthesia was injected into each slug in each of the three different injection sites in order to quickly disperse the anesthesia within the organism. Anesthesia was isotonic to artificial sea water (ASW) and contained 379 mM MgCl2 and 15 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES, a common buffering agent used to maintain pH). After verifying that there were no apparent muscular responses or gill reflexes, a ventral incision was made along the muscular foot and standardized dissecting techniques using scissors and tweezers were used to locate, isolate, and remove the ganglia (Figure 9). The ganglia were easily located in each animal as they were small orange masses located in the abdominal cavity whose connections to other ganglia in the organism could be seen. No digestive enzymes were used in the isolation of the ganglia as only the large portion of connective tissue surrounding the ganglia was removed. After cutting out around the abdominal ganglia and removing the connective tissue, the abdominal ganglia exhibit a characteristic “H” shape as seen in Figure 10. The ganglia were initially washed by being placed in a Petri dish containing isotonic ASW and allowed to adjust for several minutes. The ganglia were then moved to another Petri dish containing either isotonic, 20% hypertonic, or
36
20% hypotonic ASW. Three Aplysia were dissected for each repetition and all three ganglia were loaded into a 2.5-mm outer diameter capillary, acquired from Wilmad Lab Glass, approximately 2.5-3 cm in length which contained the ASW solution and was sealed at both ends using wax.
Figure 9. Dissection process for isolating the abdominal ganglion from Aplysia. Left: Anesthetized Aplysia prior to dissection. Middle: Aplysia after the ventral incision exposing the abdominal cavity. Right: Close-up of the abdominal ganglion which is removed.
Figure 10. Image of three viable isolated abdominal ganglia from Aplysia after dissection. All three ganglia were loaded into the same capillary and then imaged.
For fixed ganglia samples, the same dissection procedure as outlined above was used. The ganglia were initially washed by being placed in a Petri dish containing isotonic ASW and allowed to adjust. The ganglia were then placed in a 4 wt% paraformaldehyde ASW solution and refrigerated for several weeks. At the time of final sample preparation, the ganglia were then
37 washed several times with isotonic ASW to remove the as much of the fixative as possible. Finally the ganglia were moved to a Petri dish containing either isotonic, 20% hypertonic, or 20% hypotonic ASW and allowed to adjust. Again, three ganglia were loaded into a 2.5 mm outer diameter capillary containing the ASW solution and sealed with wax at both ends.
4.2 Artificial Sea Water (ASW) Solutions
Aplysia californica is a marine organism whose natural environment is sea water. Therefore, to assess osmotic changes in Aplysia neurons, artificial sea water (ASW) of similar composition is prepared. ASW has standard values of several key compounds, as shown in Table
1, including sodium chloride (NaCl), potassium chloride (KCl), magnesium chloride (MgCl2), calcium chloride (CaCl2), and HEPES, a buffering agent used to maintain pH.
Table 1. Standard concentrations of key compounds in isotonic artificial sea water (ASW). The concentration, C, is given in mmol/L.
Substance C / mM NaCl 460 KCl 10.4
MgCl2 55
CaCl2 11 HEPES 15
The osmotic pressure, defined as a measure of the tendency of water to move to a more concentrated solution, of the isotonic solution was calculated using the Morse equation, shown in (18), for osmotic pressure, assuming an ideal dilute solution and full dissociation of the salt compounds.
(18)
38
In this equation i is the van’t Hoff factor for dissociation, M is the concentration of the compound in mol/L, R is the universal gas constant (0.082 Lˑatm/molˑK), T is the absolute temperature in K, and π is the osmotic pressure in atm. When multiple compounds contribute to the osmotic pressure, this equation can be expanded as shown in (19) when n species are present.
(19)
∑
Assuming the solution to be at 25 °C (293 K) and letting species 1, 2, 3, 4, and 5 be sodium chloride, potassium chloride, magnesium chloride, calcium chloride, and HEPES, respectively, the osmotic pressure can be calculated. The van’t Hoff factors for each of these species are 2, 2, 3, 3, and 1, respectively. An example calculation for the isotonic ASW is shown in (20) below.