A two-compartment mathematical model of neuroglial metabolism using [1-(11)C] acetate.

The purpose of this study was to develop a two-compartment metabolic model of brain metabolism to assess oxidative metabolism from [1-(11)C] acetate radiotracer experiments, using an approach previously applied in (13)C magnetic resonance spectroscopy (MRS), and compared with an one-tissue compartment model previously used in brain [1-(11)C] acetate studies. Compared with (13)C MRS studies, (11)C radiotracer measurements provide a single uptake curve representing the sum of all labeled metabolites, without chemical differentiation, but with higher temporal resolution. The reliability of the adjusted metabolic fluxes was analyzed with Monte-Carlo simulations using synthetic (11)C uptake curves, based on a typical arterial input function and previously published values of the neuroglial fluxes V(tca)(g), V(x), V(nt), and V(tca)(n) measured in dynamic (13)C MRS experiments. Assuming V(x)(g)=10 × V(tca)(g) and V(x)(n)=V(tca)(n), it was possible to assess the composite glial tricarboxylic acid (TCA) cycle flux V(gt)(g) (V(gt)(g)=V(x)(g) × V(tca)(g)/(V(x)(g)+V(tca)(g))) and the neurotransmission flux V(nt) from (11)C tissue-activity curves obtained within 30 minutes in the rat cortex with a beta-probe after a bolus infusion of [1-(11)C] acetate (n=9), resulting in V(gt)(g)=0.136±0.042 and V(nt)=0.170±0.103 µmol/g per minute (mean±s.d. of the group), in good agreement with (13)C MRS measurements.

Feasibility of in vivo 15N MRS detection of hyperpolarized 15N labeled choline in rats.

The increase of total choline in tumors has become an important biomarker in cancer diagnosis. Choline and choline metabolites can be measured in vivo and in vitro using multinuclear MRS. Recent in vivo(13)C MRS studies using labeled substrates enhanced via dynamic nuclear polarization demonstrated the tremendous potential of hyperpolarization for real-time metabolic studies. The present study demonstrates the feasibility of detecting hyperpolarized (15)N labeled choline in vivo in a rat head at 9.4 T. We furthermore report the in vitro (172 +/- 16 s) and in vivo (126 +/- 15 s) longitudinal relaxation times. We conclude that with appropriate infusion protocols it is feasible to detect hyperpolarized (15)N labeled choline in live animals.

Hyperpolarized lithium-6 as a sensor of nanomolar contrast agents.

Lithium is widely used in psychotherapy. The (6)Li isotope has a long intrinsic longitudinal relaxation time T(1) on the order of minutes, making it an ideal candidate for hyperpolarization experiments. In the present study we demonstrated that lithium-6 can be readily hyperpolarized within 30 min, while retaining a long polarization decay time on the order of a minute. We used the intrinsically long relaxation time for the detection of 500 nM contrast agent in vitro. Hyperpolarized lithium-6 was administered to the rat and its signal retained a decay time on the order of 70 sec in vivo. Localization experiments imply that the lithium signal originated from within the brain and that it was detectable up to 5 min after administration. We conclude that the detection of submicromolar contrast agents using hyperpolarized NMR nuclei such as (6)Li may provide a novel avenue for molecular imaging.

Mathematical modeling of (13)C label incorporation of the TCA cycle: the concept of composite precursor function.

A novel approach for the mathematical modeling of (13)C label incorporation into amino acids via the TCA cycle that eliminates the explicit calculation of the labeling of the TCA cycle intermediates is described, resulting in one differential equation per measurable time course of labeled amino acid. The equations demonstrate that both glutamate C4 and C3 labeling depend in a predictable manner on both transmitochondrial exchange rate, V(X), and TCA cycle rate, V(TCA). For example, glutamate C4 labeling alone does not provide any information on either V(X) or V(TCA) but rather a composite "flux". Interestingly, glutamate C3 simultaneously receives label not only from pyruvate C3 but also from glutamate C4, described by composite precursor functions that depend in a probabilistic way on the ratio of V(X) to V(TCA): An initial rate of labeling of glutamate C3 (or C2) being close to zero is indicative of a high V(X)/V(TCA). The derived analytical solution of these equations shows that, when the labeling of the precursor pool pyruvate reaches steady state quickly compared with the turnover rate of the measured amino acids, instantaneous labeling can be assumed for pyruvate. The derived analytical solution has acceptable errors compared with experimental uncertainty, thus obviating precise knowledge on the labeling kinetics of the precursor. In conclusion, a substantial reformulation of the modeling of label flow via the TCA cycle turnover into the amino acids is presented in the current study. This approach allows one to determine metabolic rates by fitting explicit mathematical functions to measured time courses.

Accelerating MR elastography: a multiecho phase-contrast gradient-echo sequence.

PURPOSE: To demonstrate the feasibility of using a multiecho phase-contrast (PC) gradient-echo sequence with motion-sensitizing gradient (MSG) to accelerate MR elastography (MRE) acquisitions in comparison to single-echo PC sequences. MATERIALS AND METHODS: The sequence was implemented and compared with a conventional single-echo sequence as the standard of reference in both agarose phantoms and in vivo in the biceps of three healthy volunteers. For reconstruction of the elasticity modulus, a local frequency estimation (LFE) algorithm was used. ETL factors of 1-16 were evaluated. RESULTS: Phantom experiments demonstrated excellent consistency between single-echo and multiecho measurements in terms of wave equivalency, SNR, and reconstructed shear modulus. Additionally, the in vivo MRE examinations showed an excellent correspondence to the single-echo results. Minor loss of wave amplitude was observed at higher ETL factors. CONCLUSION: The results demonstrate that a multiecho sequence is suitable for accelerating MRE in nearly homogeneous tissue, such as muscle. It provides equivalent elasticity values in a significantly reduced scan time compared to a single-echo sequence. The maximum achievable ETL factor must be individually determined for the target tissue.

A diffusion tensor imaging analysis of gender differences in water diffusivity within human skeletal muscle.

The diffusive properties of adjacent muscles at rest were evaluated in male (n = 12) and female (n = 12) subjects using diffusion tensor imaging (DTI). The principle, second and third eigenvalues, trace of the diffusion tensor [Tr(D)], and two anisotropic parameters, ellipsoid eccentricity (e) and fractional anisotropy (FA), of various muscles in the human calf were calculated from the diffusion tensor. Seven muscles were investigated in this study from images acquired of the left calf: the soleus, lateral gastrocnemius, medial gastrocnemius, posterior tibialis, anterior tibialis, extensor digitorum longus and peroneus longus. A mathematical model was also derived that relates the eigenvalues of the diffusion tensor to the muscle fiber volume fraction, which is defined as the volume of muscle fibers within a well-defined arbitrary muscle volume. Females on average had higher eigenvalues and Tr(D) compared with males, with the majority of muscles being statistically different between the sexes. In contrast, males on average had higher e and FA than females, with the large plantar flexors--soleus, lateral gastrocnemius, and medial gastrocnemius--producing statistically different results. The behavior of the mathematical model for variations in fiber volume fraction produced similar trends to those seen when the experimental data were fit to the model. The model predicts that a larger volume fraction of skeletal muscle in males is devoted to fibers than in females, but the true underlying source of the gender discrepancy remains unclear. Although the model does not fully account for other transport processes, it does provide some insight into the limiting factors that affect the diffusion of water in skeletal muscle measured by DTI.

Cerebral activation using a MR-compatible piezoelectric actuator with adjustable vibration frequencies and in vivo wave propagation control.

Functional magnetic resonance imaging (fMRI) studies are increasingly used in patients with brain tumors near the sensory motor cortex for planning of therapy. Passive stimuli can be helpful for reproducible results. The purpose of our study was to investigate frequency and amplitude dependencies of cerebral activation patterns using a vibratory stimulus that involves sensory and motor function and allows exact adjustment of vibratory frequencies and direct control of penetration depth into the tissue. Fifteen volunteers were studied with fMRI during vibratory stimulation of the right biceps muscle utilizing a block design (frequencies: 150 and 300 Hz, amplitudes: 400, 600, and 800 microm). In addition, visualization of the wave propagation into the biceps tissue itself was performed with a modified phase contrast sequence. A specially developed MR-compatible mechanical oscillator was used to apply the vibrotactile sensations. fMRI revealed activation of the left primary somatosensory cortex during application of both vibratory frequencies. Additionally, activity of the primary and supplementary motor cortex was revealed using 150-Hz stimuli, while only minimal at 300 Hz. The activity strength correlated with increasing stimulus amplitudes and the visualized penetration depth. Activation of motor areas was more pronounced at the beginning of the rest period. In conclusion, sensory motor areas can be activated using a piezoelectric actuator, with less pronounced motor area activation at higher frequencies. Our setup allowed local control of stimulus penetration through the tissue correlated to central activation, providing objective stimulus control. The pronounced activation of the motor cortex during the rest condition may reflect the subjective feeling of arm movement after the end of the stimulus.

Diffusive sensitivity to muscle architecture: a magnetic resonance diffusion tensor imaging study of the human calf.

The aim of this study was to examine the diffusive properties of adjacent muscles at rest, and to determine the relationship between diffusive and architectural properties, which are task-specific to muscles. The principle, second, and third eigenvalues, trace of the diffusion tensor, and two anisotropic parameters, ellipsoid eccentricity (e) and fractional anisotropy (FA), of various muscles in the human calf were calculated by diffusion tensor imaging (DTI). Linear correlations of the calculated parameters to the muscle physiological cross-sectional area (PCSA), which is proportional to maximum muscle force, were performed to ascertain any linear relation between muscle architecture and diffusivity. Images of the left calf were acquired from six healthy male volunteers. Seven muscles were investigated in this study. These comprised the soleus, lateral gastrocnemius, medial gastrocnemius, posterior tibialis, anterior tibialis, extensor digitorum longus, and peroneus longus. All data were presented as the mean and standard error of the mean (SEM). In general, differences in diffusive parameter values occurred primarily between functionally different muscles. A strong correlation was also found between PCSA and the third eigenvalue, e, and FA. A mathematical derivation revealed a linear relationship between PCSA and the third eigenvalue as a result of their dependence on the average radius of all fibers within a single muscle. These findings demonstrated the ability of DTI to differentiate between functionally different muscles in the same region of the body on the basis of their diffusive properties.

In vivo elasticity measurements of extremity skeletal muscle with MR elastography.

MR elastography (MRE) has been shown to be capable of non-invasively measuring tissue elasticity even in deep-lying regions. Although limited studies have already been published examining in vivo muscle elasticity, it is still not clear over what range the in vivo elasticity values vary. The present study intends to produce further information by examining four different skeletal muscles in a group of 12 healthy volunteers in the age range of 27-38 years. The examinations were performed in the biceps brachii, the flexor digitorum profundus, the soleus and the gastrocnemius. The average shear modulus was determined to be 17.9 (+/- 5.5), 8.7 (+/- 2.8), 12.5 (+/- 7.3) and 9.9 (+/- 6.8) kPa for each muscle, respectively. To ascertain the reproducibility of the examination, the stiffness measurements in two volunteers were repeated seven times for the biceps brachii. These examinations yielded a mean shear modulus of 11.3 +/-.7 and 13.3 +/- 4.7 kPa for the two subjects. For elasticity reconstruction, an automated reconstruction algorithm is introduced which eliminates variation due to subjective manual image analysis. This study yields new information regarding the expected variation in muscle elasticity in a healthy population, and also reveals the expected variability of the MRE technique in skeletal muscle.