Central carbon metabolism plays an important role in the origin or progression of many diseases, including cancers, neurodegenerative diseases, diabetes, and nonalcoholic fatty liver disease (NAFLD) (Zilberter, Y. & Zilberter, M. J Neurosci Res. 2017 95:2217-2235; Pavlova, N. N. & Thompson, C. B. Cell Metabolism 2016 23:27-47; Samuel, V T. & Shulman, G I. Cell Metabolism 2018 27:22-41). Metabolic reprogramming is one of the key features of cancer cell metabolism. Many cancer cells consume tremendous amounts of glucose and generate ATP via glycolysis even under normoxic conditions (the Warburg effect) (Bensinger, S J. & Christofk, H R. Semin Cell Dev Biol. 2012 23:352-361; Warburg, O. Science 1956 124:269-270; Warburg, O. The J Cancer Res 1925 9:148-163). Many theories have been promulgated to explain this observation, ranging from impaired mitochondrial function as the primary causative agent of cancer (Warburg, O. Science 1956 123:309-314), to a trade-off of inefficient energy metabolism for faster incorporation of carbon into biomass (Vander Heiden, M G., et al. Science 2009 324:1029-1033), to a double selective advantage of hypoxia resistance and generation of local acidosis (Gatenby, R A. & Gillies, R J. Nat Rev Cancer. 2004 4:891-899). Molecular imaging techniques capable of identifying such differences have become very important for cancer diagnosis, disease progression, and treatment.
The kinetics of lactate and HDO (or 2HOH) production can serve as a quick and reliable metric for cancer diagnosis or staging. Glucose tracers have been widely utilized for accessing the metabolic conversion of glucose to lactate predominantly via glycolysis. Metabolic imaging with radioactive agents, which relies on physiologic hallmarks unique to cancer (Hanahan, D. & Weinberg, R A. Cell 2011 144:646-674), may be more sensitive but its widespread use is limited. In the clinic, the glucose analog [18F]2-fluoro-2-deoxy-glucose (18FDG) is used for positron emission tomography (PET) high-resolution maps of glucose uptake, but cannot report differences in metabolism and is associated with doses of ionizing radiation that carry an unknown, but real, risk of secondary cancer (Linet, M. S. et al. C A Cancer J Clin. 2012 62(2):75-100; Kaushik, A. et al. Indian J Med Res. 2015 142(6):721-31). These tests are therefore limited to patients with known or suspected malignancy. 18FDG is radioactive, which is also a limitation when repetitive scanning is required for evaluating disease progression (Kim, M M., et al. Nat Rev Clin Oncol 2016 13:725-739). CEST MRI was developed to overcome some of the limitations of PET for glucose imaging. CEST detects glucose uptake indirectly by saturation transfer from the exchangeable protons of glucose to water, which enhances sensitivity relative to direct detection, but the technique is hampered by water exchange rates at physiological pH (Walker-Samuel, S. et al. Nat Med 2013 19:1067-1072; van Zijl, P C., et al. Proc Natl Acad Sci USA. 2007 104(11):4359-64). Dissolution dynamic nuclear polarization (dDNP) was developed to enhance signal-to-noise ratio (SNR) in 13C MR via polarization transfer from an electron source to nuclei, resulting in >10,000 times signal enhancement in the MRI signal (Ardenkjaer-Larsen, J H. et al. Proc Natl Acad Sci USA. 2003 100:10158-10163; Kurhanewicz, J. et al. Neoplasia 2011 13:81-97). However, dDNP is limited to a small set of molecules having long relaxation times to measure the kinetics before polarized signal is lost. Many important metabolites, i.e. glucose, have short relaxation times, and are difficult to image using dDNP methods (Comment A, et al. Biochemistry. 2014 53(47):7333-57). Lactate pool size convolves with estimates of lactate production when estimated with hyperpolarized pyruvate, which in a strict biochemical sense is not measuring glycolysis. Analysis of lactate isotopomers has been carried out using 13C labeled glucose and spin editing, but the 13C enrichment is not easily determined in some cases due to overlap with fatty acid (FA) resonances (Jucker, B M., et al. J. Biol. Chem. 1997 272:10464-10473). Due to the multiplicities from nuclear spin-spin couplings of the 1H-13C in lactate 13C-isotopomers, 1H-NMR spectra become complex. Direct 13C acquisition provides a straightforward method for isotopomer analysis but requires long data acquisition time (Lloyd, S G., et al. Magn. Reson. Med. 2004 51:1279-1282). Therefore, improved methods for in vivo detection of glucose utilization by tumor cells are needed.
Non-alcoholic fatty liver disease (NAFLD) and non-alcoholic steatohepatitis (NASH) are public health risks that are rapidly growing (Younossi, Z, et al. Nature reviews Gastroenterology & Hepatology 2018 15(1):11). NASH, the more medically serious condition, is predicted to become the number one cause of liver transplant in the next decade (McPherson, S, et al. Journal of hepatology 2015 62(5):1148-1155). Liver biopsy has remained the primary criterion standard or “gold standard” in the evaluation of the etiology and extent of disease of the liver such as NAFLD and NASH. Percutaneous liver biopsy is the preferred method to determine NAFLD and to differentiate NASH from NAFLD (Nalbantoglu, I, et al. World journal of gastroenterology 2014 20(27):9026). Other biopsy methods are typically even more invasive and include transvenous and laparoscopic liver biopsy. Although liver biopsies are generally regarded as safe, they bare risks that are potentially lethal. Almost two thirds of complications of liver biopsy occur within two hours. Approximately 2% of patients undergoing liver biopsy require hospitalization for the management of an adverse event. Significant bleeding after a liver biopsy occurs in approximately 1% of patients who are biopsied (Perrault, J, et al. Gastroenterology 1978 74(1):103-106). If bleeding persists, a blood transfusion may be needed. Surgery or angiography, where the bleeding site is identified and treated, may be required if the bleeding is severe or does not stop on its own. Intraperitoneal hemorrhage is the most serious consequence of bleeding. Liver biopsy is also handicapped by sampling error, which can result in misdiagnosis (Ratziu, V, et al. Gastroenterology 2005 128(7):1898-1906). Thus, a robust minimally invasive test to reliably and efficiently diagnose NASH and differentiate the inflammatory indication NASH from the less harmful NAFLD, and which also provides an assessment of the entire liver, is needed.
As disclosed herein, deuterium (2H) magnetic resonance imaging (DMI), a new method for assessing metabolic flux in functioning tissues, cell culture, and in vivo, provides the requisite chemical, spatial and temporal resolution needed for detecting glucose metabolism in tissues and cancers (De Feyter, H M. et al. et al. Sci Adv 2018 4:eaat7314). Nonradioactive 2H-labeled substrates can be utilized to observe the downstream metabolites, due to very low natural abundance of 2H (0.0115%), requiring no natural abundance correction for analysis. Administration of deuterium-labeled substrates at tracer levels (<200 mg/kg 2H) is well established as safe (Klein, D. Stable Isotopes 1986 378-382; Jones, P J H. & Leatherdale, S T. Clin Sci (Lond). 1991 80(4):277-80; Koletzko, B., et al. Eur J Pediatr. 1997 156 Suppl 1:S12-7). In vivo 2H MR was investigated in the late 1980's (Ackerman, J J H., et al. Proc Natl Acad Sci USA. 1987 84:4099-4102; Ewy, C S., et al. Magnetic Resonance in Medicine 1988 8:35-44; Ewy, C S., et al. Magnetic Resonance Imaging 1986 4:407-411; Müller, S. & Seelig, J. Journal of Magnetic Resonance 1969 72:456-466), but has recently enjoyed a resurgence due to improved technology and the accumulation of decades of basic research on cancer metabolic phenotypes (Tee, S. & Keshari, K R. The Cancer Journal 2015 21:165-173). Previous studies using 2H MRI to image cancer metabolism relied on differential production of TCA cycle intermediates and lactate (De Feyter, H M. et al. Sci Adv. 2018 4(8):eaat7314). Deuterated tracers are inexpensive and straightforward to administer, but 2H MR sensitivity is low. Indeed, the experiment is only tractable because the quadrupole of the 2H nucleus makes its T1 short, facilitating very short recycle time (Tr) To cope with low sensitivity of 2H-NMR experiment for less abundant metabolites, more sensitive experiments and methodology were designed.
As disclosed herein, enrichment of 2HOH (i.e. HDO production) after administration of deuterium-labeled glucose will provide a sensitive marker of glucose utilization in tissues and tumor cells in vivo, with a possibility for a simplified imaging strategy. Likewise, administration of deuterium-labeled fatty acid will provide a sensitive marker of fatty acid (FA) beta-oxidation. All living tissues are oxidative in nature. Therefore, one could expect both glucose and any fatty acid to produce HDO by the disclosed methods in any functioning tissue including the brain, heart, muscle, skin, liver, pancreas, spleen, kidney, stomach, gut, or lungs. HDO generated from oxidation of glucose or a fatty acid could be used to detect either an increment or decrement of metabolic activity in any organ or tissue. These changes could be derived from inborn errors in metabolism.
Therefore, disclosed herein is a method for imaging a tissue in a subject, comprising administering to the subject a composition comprising deuterium-labeled glycolytic of FA substrate; and imaging the subject with deuterium magnetic resonance imaging (DMI) to detect hydrogen-deuterium oxide (HDO) in tissues of the subject, wherein detection of HDO in a location within the tissue of the subject is an indication of glucose utilization or FA oxidation in the tissue of the subject, respectively.
This method could therefore be used to detect any disease in the brain that results in changing metabolism as well as traumatic brain injury. This method could also be used to detect normal brain activation, in analogy to functional MRI (fMRI). As another example, HDO generation could be used to detect any increment/decrement in heart function, due to exercise, or due to heart hypertrophy or failure. Heart dysfunction would be detected by HDO generation from glucose and/or a fatty acid whether it was caused by over-pressure heart failure or via diabetic cardiomyopathy. HDO generation could be used to detect any increment/decrement in muscle function or in exercise physiology. Any disease of the muscle that changed its metabolism would be sensitive to HDO generation from glucose or a FA. HDO generated from glucose metabolism in the pancreas could also be used to detect glucose stimulated insulin secretion. HDO generated from glucose or FA metabolism could also be used to detect kidney function.
Lactate production from glucose is a hallmark of many cancers. The disclosed methods can therefore also be used to detect cancer based on glycolytic flux due to this Warburg Effect. If the glycolytic product is enriched with deuterium, HDO is produced along with the lactate, which can then be detected. The disclosed methods involve first administering to the subject a deuterium-labeled glycolytic substrate (e.g. glucose), which will be metabolized into lactate and hydrogen-deuterium oxide (HDO) in cancers. In addition, deuterium-induced perturbations of 1H-NMR chemical shifts (isotope shift) allows resolution of the lactate isotopomers.
Therefore, disclosed herein is a method for detecting a tumor in a subject, that involves administering to the subject a composition comprising deuterium-labeled glycolytic substrate; and imaging the subject with deuterium magnetic resonance imaging (DMI) to detect hydrogen-deuterium oxide (HDO) in tissues of the subject, wherein detection of elevated HDO in a location within the tissue of the subject is an indication of the presence of cancerous tissue in the subject.
In some embodiments, the deuterium-labeled glycolytic substrate comprises deuterated glucose. In preferred embodiments, the deuterated glucose comprises at least 3, 4, 5, 6, or 7 deuterium atoms. Therefore, in some embodiments, the deuterated glucose comprises [2H7]glucose.
In some embodiments, the DMI step is further used to detect deuterated lactate production in the tissue of the subject by detecting deuterium-induced perturbations of 1H-NMR chemical shifts. However, in preferred embodiments, the method does not involve detecting deuterated lactate production. This method is less viable, as the shifts are only easily detected using a two-stage experiment with acquisition with and without 2H spin decoupling. This will make the experiment twice as long, and also requires 2H decoupling, which is difficult to implement in humans due to issue with specific absorbance rate (SAR). The power associated with 2H decoupling can trigger SAR safety issues, as it can cause local heating in the subject and is therefore protected against.
Also disclosed herein is a metabolic imaging modality that can be used to screen and follow high-risk patients (i.e., obese patients) and diagnose early on the transition from “benign” NAFLD to NASH and allow prompt intervention to prevent severe liver disease.
Both NAFLD and NASH are characterized by aberrant fatty acid (FA) beta-oxidation in the liver. FA oxidation can be detected using 13C based magnetic resonance, but the protocol is cost prohibitive and relies upon a very long examination period for the subject that is clinically intractable. Fatty acid uptake can be estimated with 11C positron emission tomography, but this method cannot directly measure oxidative flux. However, when beta-oxidation proceeds in any living tissue, it generates hydrogen from the precursor FA. If the FA, of any chain length or specific chemistry in regards to double bond positioning, is perdeuterated, then the eliminated hydrogen is instead a deuterium (2H), which causes the generation of HDO instead of H2O. This HDO signal can be imaged via T1, weighted 1H MRI or by 2H based imaging. Deuterium magnetic resonance imaging (DMI) in particular provides the requisite chemical, spatial and temporal resolution needed for detecting HDO signal (De Feyter, H M. et al. et al. Sci Adv 2018 4:eaat7314).
Therefore, detection of HDO following administration of 2H labeled FA precursors, via any route, oral, intravenous, or intraperitoneal, should correlate with increasing or decreasing beta-oxidation. Administration of deuterium-labeled substrates at tracer levels (<200 mg/kg 2H) is well established as safe (Klein, D. Stable Isotopes 1986 378-382; Jones, P J H. & Leatherdale, S T. Clin Sci (Lond). 1991 80(4):277-80; Koletzko, B., et al. Eur J Pediatr. 1997 156 Suppl 1:S12-7). By developing 2H magnetic resonance imaging (MRI) as a sensitive detector of increased hepatic fat oxidation that accompanies the NAFLD/NASH transition, this method has the potential to radically improve patient care and have major clinical translational value for millions of affected individuals.
As disclosed herein, enrichment of 2HOH (i.e. HDO production) after administration of deuterium-labeled FA provides a sensitive marker of liver FA beta-oxidation in vivo, with a possibility for a simplified imaging strategy.
Therefore, disclosed herein is a method for detecting FA beta-oxidation in the liver of a subject, that involves administering to the subject a composition comprising deuterium-labeled FA substrate; and imaging the subject with deuterium magnetic resonance imaging (DMI) to detect hydrogen-deuterium oxide (HDO) in the liver of the subject, wherein detection of elevated HDO in a location within the liver of the subject is an indication of the presence of aberrant FA beta-oxidation in the subject, which is characteristic of NAFLD and NASH. In some embodiments, the DMI comprises 2H-NMR imaging. In some embodiments, the DMI comprises 1H-NMR imaging that uses T1 weighting.
As very long chain fatty acids can be oxidized in the peroxisome, which is not thought to share a similar role as mitochondria in NAFLD and NASH, the deuterated fatty acid substrate is preferably a small chain fatty acid (SCFA), a medium-chain fatty acid (MCFA), or a long-chain fatty acid (LCFA), but not a very long chain fatty acid (VLCFA). Therefore, in some embodiments, the fatty acid substrate is a C6-C19 fatty acid, such as a C6, C7, C8, C9, C10, C11, C12, C13, C14, C16, C17, C18, C19, CD20, or CD21 fatty acid.
In some embodiments, the fatty acid substrate is a saturated fatty acid. In some embodiments the fatty acid substrate is an unsaturated fatty acid. In some embodiments, the fatty acid substrate has 0, 1, 2, 3, 4, or 5 double bonds. In some embodiments, the fatty acid substrate is even-chained fatty acids (ECFA). In some embodiments, the fatty acid substrate is odd-chained fatty acids (OCFA). In some embodiments, the fatty acid substrate is a non-esterified fatty acids (NEFA) or free fatty acids (FFA).
In some embodiments, the substrate is composed of a partial or fully deuterated triglyceride (TRG). A TRG contains three fatty acid tails that can be the same length or different lengths. As disclosed herein, the TRG can have any combination of fatty acid tail lengths (e.g. C6-C21), wherein 1, 2, or 3 of the tails are partial or fully deuterated.
In some embodiments, the deuterium-labeled FA substrate comprises deuterated octanoate, which bypasses fatty acid transport mechanisms and is rapidly oxidized.
In some embodiments, the deuterated fatty acid is perdeuterated. In some embodiments, the deuterated fatty acid comprises at least 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, or 41 deuterium atoms (or up to 3×these numbers if a triglyceride).
One advantage of using HDO, e.g. as the surrogate of glucose utilization by cancers, or as a surrogate of beta-oxidation in the liver, is the reduced imaging time that is needed for imaging (DMI). In some embodiments, the imaging step(s) occurs and/or are completed within 60 minutes of the administration of the deuterated tracer, including 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 39, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, or 60 minutes of the administration of the deuterated tracer. As shown herein, the case of glucose, the imaging step(s) preferably occur in a window from 8 to 24 minutes. By choosing the correct imaging window temporally, it is possible to produce images that more accurately estimate metabolic flux without contamination from peripheral metabolism of the administered deuterated tracer. DMI signal to noise ratio (SNR) is defined by the highest peak intensity versus the noise. HDO is the strongest peak in the spectrum. By detecting it, e.g. instead of lactate, the peak with the highest SNR is being used, reducing acquisition time. Also, imaging using the HDO peaks allows other imaging methods than chemical shift imaging (CSI), an imaging sequence that is not easily accelerated by methods like sensitivity encoding for fast MRI (SENSE) (Pruessmann, Klaas P., et al. Magnetic Resonance in Medicine 42.5 1999 952-962.
In some embodiments, chemically selective imaging is accomplished using a suitable MRI pulse sequence. Chemically selective imaging allows the concentration of the starting substrate, e.g. glucose or a FA, to be ascertained apart from the HDO signal intensity or from other downstream metabolic products, e.g. of glucose or fatty acid metabolism.
Chemical shift imaging (CSI) is described in de Feyter, et. Al., Science advances. 2018 4(8):eaat7314, which is incorporated by reference herein for the teaching of this method. Further improvements upon CSI that increase the speed of imaging can be implemented, including spiral-CSI (Adalsteinsson E, et al. Magnetic resonance in medicine. 1998 39(6):889-98, which is incorporated by reference herein for the teaching of this method), CSI using centric sampling (Wilman A H, et al. Magnetic resonance in medicine. 1997 38(5):793-802, which is incorporated by reference herein for the teaching of this method), and CSI using compressed sensing for reconstruction (Lustig M, et al. Magnetic Resonance in Medicine. 2007 58(6):1182-95, which is incorporated by reference herein for the teaching of this method). Echo planar spectroscopic imaging (EPSI) is also possible (Mansfield P. Magnetic Resonance in Medicine. 1984 1(3):370-86, which is incorporated by reference herein for the teaching of this method). Newer implementations of EPSI also exist (Chen A P, et al. Magn Reson Med. 2007 58(6):1099-106; Jakob P M, et al. Magnetic resonance in medicine. 1995 33(4):573-8, which are incorporated by reference herein for the teaching of this method).
An alternative method for producing chemically selective images is to use accumulated phase and multiple sampling windows to produce images which can be combined to produce chemically selective images, commonly referred to as the Dixon method (Dixon W T. Radiology. 1984 153(1):189-94, which is incorporated by reference herein for the teaching of this method). This technique uses two images to produce in phase and out of phase combinations of two metabolites. It can be adapted for 3 resonances as well by using more echoes (Szumowski J, et al. Radiology. 1994 192(2):555-61, which is incorporated by reference herein for the teaching of this method). More sophisticated methods can use even more echoes (Wang J x, et al. Magnetic Resonance in Chemistry. 2016 54(8):665-73, which is incorporated by reference herein for the teaching of this method). This many-echo method is commonly referred to as Iterative decomposition of water and fat with echo asymmetry and least squares estimation (IDEAL) imaging.
An alternative method uses selective excitation of the desired resonance followed by imaging to produce chemically selective images. Long selective pulses can be used to excite the single resonance peak (Meyer C H, et al. Magnetic resonance in medicine. 1990 15(2):287-304, which is incorporated by reference herein for the teaching of this method), but these pulses are subject to RF homogeneity issues. Spectral-spatial pulses are widely used for this application and have superior properties for slice selection (Cunningham C H, et al. Journal of magnetic resonance. 2008 193(1):139-46; Sigfridsson A, et al. Magnetic Resonance in Medicine. 2015 73(5):1713-7, which are incorporated by reference herein for the teaching of this method).
Another method for imaging that could be suitable combines selective excitation with steady-state free precession (SSFP) to rapidly acquire images (Scheffler K, et al. Eur Radiol. 2003 13(11):2409-18, which is incorporated by reference herein for the teaching of this method). These sequences are rapid and produce excellent in plane resolution.
Deuterium MRI is handicapped by the short T2 of the nucleus. Methods that utilize very short time to echoes could also be applicable in certain circumstances (Hafner S. Magn Reson Imaging. 1994 12(7):1047-51; Robson M D, et al. Journal of computer assisted tomography. 2003 27(6):825-46, which are incorporated by reference herein for the teaching of this method). A similar technique is sweep imaging (ldiyatullin D, et al. Journal of Magnetic Resonance. 2006 181(2):342-9, which is incorporated by reference herein for the teaching of this method).
Each of these methods focus on producing images that separate HDO from any other resonances that would contaminate the image, reducing its diagnostic power.
In some embodiments, the method involves a single imaging step that detects elevated HDO in a tissue, e.g. by comparing it to other tissues. In other embodiments, the method involves imaging the subject at multiple time points in order to determine a rate of HDO enrichment in the tissue of the subject, to detect regional differences in glycolytic flux associated with cancer or other disease. Using kinetic modeling to estimate metabolic turnover has been used to study in vivo metabolism in many contexts (Cerdan S, et al. J Biol Chem. 1990 265(22):12916-26; Patel A B, et al. Proc Natl Acad Sci USA 2005 102(15):5588-93; Sibson N R, et al. Proc Natl Acad Sci USA 1997; 94(6):2699-704, which are incorporated by reference herein for the teaching of this method). These methods measure enrichment as a function of time to measure rates quantitatively. They are endemically more accurate than taking a single image, but require longer scan times for the patient. Patient scan time minimization is essential, as scans longer than 30 minutes are poorly tolerated by the ill.
The details of one or more embodiments of the invention are set forth in the accompanying drawings and the description below. Other features, objects, and advantages of the invention will be apparent from the description and drawings, and from the claims.
Before the present disclosure is described in greater detail, it is to be understood that this disclosure is not limited to particular embodiments described, and as such may, of course, vary. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting, since the scope of the present disclosure will be limited only by the appended claims.
Where a range of values is provided, it is understood that each intervening value, to the tenth of the unit of the lower limit unless the context clearly dictates otherwise, between the upper and lower limit of that range and any other stated or intervening value in that stated range, is encompassed within the disclosure. The upper and lower limits of these smaller ranges may independently be included in the smaller ranges and are also encompassed within the disclosure, subject to any specifically excluded limit in the stated range. Where the stated range includes one or both of the limits, ranges excluding either or both of those included limits are also included in the disclosure.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this disclosure belongs. Although any methods and materials similar or equivalent to those described herein can also be used in the practice or testing of the present disclosure, the preferred methods and materials are now described.
All publications and patents cited in this specification are herein incorporated by reference as if each individual publication or patent were specifically and individually indicated to be incorporated by reference and are incorporated herein by reference to disclose and describe the methods and/or materials in connection with which the publications are cited. The citation of any publication is for its disclosure prior to the filing date and should not be construed as an admission that the present disclosure is not entitled to antedate such publication by virtue of prior disclosure. Further, the dates of publication provided could be different from the actual publication dates that may need to be independently confirmed.
As will be apparent to those of skill in the art upon reading this disclosure, each of the individual embodiments described and illustrated herein has discrete components and features which may be readily separated from or combined with the features of any of the other several embodiments without departing from the scope or spirit of the present disclosure. Any recited method can be carried out in the order of events recited or in any other order that is logically possible.
Embodiments of the present disclosure will employ, unless otherwise indicated, techniques of chemistry, biology, and the like, which are within the skill of the art.
The following examples are put forth so as to provide those of ordinary skill in the art with a complete disclosure and description of how to perform the methods and use the probes disclosed and claimed herein. Efforts have been made to ensure accuracy with respect to numbers (e.g., amounts, temperature, etc.), but some errors and deviations should be accounted for. Unless indicated otherwise, parts are parts by weight, temperature is in ° C., and pressure is at or near atmospheric. Standard temperature and pressure are defined as 25° C. and 1 atmosphere.
Before the embodiments of the present disclosure are described in detail, it is to be understood that, unless otherwise indicated, the present disclosure is not limited to particular materials, reagents, reaction materials, manufacturing processes, or the like, as such can vary. It is also to be understood that the terminology used herein is for purposes of describing particular embodiments only, and is not intended to be limiting. It is also possible in the present disclosure that steps can be executed in different sequence where this is logically possible.
It must be noted that, as used in the specification and the appended claims, the singular forms “a,” “an,” and “the” include plural referents unless the context clearly dictates otherwise.
The term “subject” refers to any individual who is the target of administration or treatment. The subject can be a vertebrate, for example, a mammal. Thus, the subject can be a human or veterinary patient. The term “patient” refers to a subject under the treatment of a clinician, e.g., physician.
Disclosed herein is a method for detecting a tumor in a subject, that involves administering to the subject a composition comprising a deuterium-labeled glycolytic substrate; and imaging the subject with deuterium magnetic resonance imaging (DMI) to detect hydrogen-deuterium oxide (HDO) in tissues of the subject, wherein detection of elevated HDO in a location within the tissue of the subject is an indication of the presence of a tumor in the tissue of the subject. The likelihood that D2O can be formed by metabolism or exchange is very low, but its formation would not preclude use of this technology as stated.
Also disclosed herein is a method for detecting fatty acid beta-oxidation in the liver of a subject, that involves administering to the subject a composition comprising a deuterium-labeled fatty acid substrate; and imaging the subject with deuterium magnetic resonance imaging (DMI) to detect hydrogen-deuterium oxide (HDO) in liver of the subject, wherein detection of elevated HDO in a location within the tissue of the subject is an indication of the presence of fatty acid beta-oxidation of the subject.
In some embodiments, the deuterium-labeled substrate is administered to the subject as an oral composition. Compositions for oral administration include powders or granules, suspensions or solutions in water or non-aqueous media, capsules, sachets, or tablets. Thickeners, flavorings, diluents, emulsifiers, dispersing aids or binders may be desirable.
Formulations can be prepared using a pharmaceutically acceptable “carrier” composed of materials that are considered safe and effective and may be administered to an individual without causing undesirable biological side effects or unwanted interactions. The “carrier” is all components present in the pharmaceutical formulation other than the active ingredient or ingredients. The term “carrier” includes but is not limited to diluents, binders, lubricants, disintegrators, fillers, matrix-forming compositions and coating compositions.
The herein disclosed compositions, including pharmaceutical composition, may be administered in a number of ways depending on whether local or systemic treatment is desired, and on the area to be treated. For example, the disclosed compositions can be administered intravenously, intraperitoneally, intramuscularly, subcutaneously, intracavity, or transdermally. The compositions may be administered orally, parenterally (e.g., intravenously), by intramuscular injection, by intraperitoneal injection, transdermally, extracorporeally, ophthalmically, vaginally, rectally, intranasally, topically or the like, including topical intranasal administration or administration by inhalant.
In some embodiments, deuterium-labeled substrate is administered in a dose equivalent to parenteral administration of about 1 mg to about 10 g per kg of body weight, including about 1 mg to about 1 g, about 1 g to about 10 g, about 10 mg to about 1 g, and about 100 mg to about 10 g.
Parenteral administration of the composition, if used, is generally characterized by injection. Injectables can be prepared in conventional forms, either as liquid solutions or suspensions, solid forms suitable for solution of suspension in liquid prior to injection, or as emulsions. A revised approach for parenteral administration involves use of a slow release or sustained release system such that a constant dosage is maintained. See, e.g., U.S. Pat. No. 3,610,795, which is incorporated by reference herein.
Deuterium metabolic imaging (DMI) is a noninvasive approach that combines deuterium magnetic resonance spectroscopic imaging with oral intake or intravenous infusion of nonradioactive 2H-labeled substrates to generate three-dimensional metabolic maps. Deuterium NMR (2H-NMR) is NMR spectroscopy of deuterium (2H or D), an isotope of hydrogen. Deuterium is an isotope with spin=1, unlike hydrogen which is spin=½.
Deuterium NMR has a range of chemical shift similar to proton NMR but with poorer resolution.
Deuterium usually yields broad signals whose line-width typically varies between a few hertz and a few kilohertz. The spectrum has the same narrow chemical shift range as for H but with a lower resolution and lower sensitivity. Deuterium-deuterium couplings are about 40 times smaller than proton-proton couplings and are therefore not observed.
A number of embodiments of the invention have been described.
Nevertheless, it will be understood that various modifications may be made without departing from the spirit and scope of the invention. Accordingly, other embodiments are within the scope of the following claims.
Introduction
In this Example, highly sensitive 1H-NMR was used to measure the production of lactate isotopomers and 2H-NMR to measure HDO production and residual [2H7]glucose (eventually cellular glucose uptake) in the cell culture medium. The fate of the 2H tracer for lactate isotopomers, HDO production, and loss of 2H atoms from [2H7]glucose during glycolysis, sugar isomerization and pentose phosphate pathway (PPP) can be easily understood from
Methods
Cell Line, Chemicals, and Media
HUH-7 hepatocellular carcinoma cells were received. AML12 cell lines were purchased from the American Type Culture Collection (Manassas, Va.). Pyrazine, Pyrazine-D4 and [1,2,3,4,5,6,6-2H7]-D-glucose ([2H7]glucose) were purchased from Sigma Aldrich, (St. Louis, Mo., USA).
Cell Culture
HUH-7 or AML12 cell lines were cultured in a complete growth medium made up of Dulbecco's Modified Eagle Medium (DMEM) with 10% v/v FBS, 10,000 units/mL penicillin, 10,000 μg/mL streptomycin, 25 μg/mL amphotericin B and 5.5 mM [2H7]glucose. Cell lines were maintained at 37° C. in a 5% CO2 in air atmosphere in an air-jacketed incubator (Heracell Vios 160i, Thermo Scientific, Waltham, Mass.). Growth media was replaced every 3 days, and subcultured 1:5 at 80% confluence. 25 cm2 cell culture plates, containing 15 million HUH-7 or AML12 cells at 70-80% confluence, were washed once with warm PBS and incubated with 4 mL each of DMEM with deuterated tracer for 5h, withdrawing 100 μL aliquots at 20 min, 1h, 2h, and 5h.
Sample Preparation
No extraction was performed on the cell media samples. Pyrazine and pyrazine-D4 were mixed to achieve 5 mM concentration of each of the standards in each cell media samples of HUH-7 and AML12 cells. Pyrazine and pyrazine-D4 served as an internal standard to quantify the lactate isotopomers from 1H-NMR data and HDO as well as residual [2H7]glucose from 2H-NMR data respectively. 40 μL of cell medium was transferred into 1.7 mm NMR tubes for NMR analysis. 1H-NMR Spectroscopy
All NMR spectra were recorded on a 600.23 MHz NMR, equipped with a 1.7 mm TCI CryoProbe and Avance Neo Console (Bruker Biospin). All NMR experiments were set up in the lock off followed by sweep off mode. 2H-decoupling off and on 1H-NMR spectra were acquired with a sequence that alternated decoupling with each scan. Relaxation delay (d1) of 1.5s and acquisition time (AQ) of 2s (3.5s of recycling time, Tr) and a 90° pulse of 10.50 μs was used for acquisition of each of the spectra. The WET method containing selective pulses was applied on the strong water resonance (M o, H. & Raftery, D. J. Biomol. NMR 2008 41:105-111). Waltz16 2H-decoupling sequence consists of 90° pulse of 250 μs with decoupling power of 1.72watts, producing B, of 1 kHz during the acquisition period. 13,157 complex data points were acquired with the spectral width of 11 ppm. 128 scans were used to acquire each spectrum. 4 dummy scans were also used for equilibration of the spin states prior to acquisition. A modified T1 inversion recovery pulse sequence with 2H-decoupling on during the acquisition was used to measure the T1 of the lactate isotopomers and pyrazine standard. Relaxation delay (d1) of 50s and inversion recovery delays (T) of 0.001, 0.6, 0.9, 1.5, 2.0, 3.0, 6.0, 9.0, 14.0, 20.0, 30.0 and 40.0s were used to determine the T1 of lactate isotopomers.
2H-NMR spectroscopy
The deuterium lock channel was used to acquire the 1H-decoupled 2H-NMR spectrum at 92.13 MHz resonant frequency. Relaxation delay (d1) of 2s and acquisition time (AQ) of 1s (3s of recycling time, Tr) with a 90° pulse was used to acquire each of the 2H-NMR spectra. 1086 complex data points were digitized with the spectral width of 11 ppm using 1024 scans for each of the 5 FIDs (5120 scans) for each of the samples. T1 inversion recovery pulse sequence was used to calculate the T1 relaxation times of HDO and pyrazine-D4 standard. Relaxation delay (d1) of 5s, acquisition time (AQ) of 1s and inversion recovery delays (T) of 0.001, 0.1, 0.4, 0.8, 1.6, and 3.2s were used to acquire an array of FIDs to calculate the T1 of HDO and pyrazine-D4. All NMR experiments were carried out at the room temperature (25° C.). 1H-NMR Data Processing and Quantification of Lactate Isotopomers
NMR data analysis was performed using MestReNova v14.0.1-23284 (Mestrelab Research S. L.). For processing of the 1H-NMR spectrum, the spectra were zero-filled to 32,768 points with an exponential window function of 0.3 Hz before the Fourier Transform (FT). Manual phase and automatic spline baseline correction were performed on each of the spectra. 2H-decoupled 1H-NMR spectra of each time point were used for the extraction of the areas of lactate isotopomers. The MestReNova line fitting tool was used to extract the peak areas of the slightly overlapped signals for quantification of the lactate isotopomers. A T1 correction factor for each lactate isotopomer was used to calculate concentrations using the 1H-NMR data (Bharti, S. K. et al. Metabolomics 4:367-376 (2008)).
2H-NMR Data Processing and Quantification of HDO and Residual [2H7]Glucose
2H-NMR spectra were processed with zero-filling of FID to 4096 points followed by exponential window function of 1 Hz before the Fourier Transform (FT). Each spectrum was manually phase corrected and followed by automatic spline baseline correction. Five FIDs were acquired for each of the samples and all five 2H-NMR spectra were aligned and summed-up to compensate the peak shifting due to magnetic field drifting over the course of time. Concentrations of HDO and residual [2H7]glucose in the cell media were calculated using peak areas of the internal standard pyrazine-D4 (concentration 5 mM), HDO and residual [2H7]glucose, normalized with the number of deuterons responsible for corresponded resonances. Baseline 2H content of ultrapure water was also measured with laser cavity ring-down spectroscopy (CRDS) (Picarro Science Instruments).
Results 1H-NMR spectra of the cell culture media
1H-NMR spectra in 2H-decoupling off and on mode were recorded for all of the withdrawn cell media. NMR spectra from the 5 hour time point show pronounced isotope shifts for lactate-CH2D and lactate-CHD2 isotopomers (
2H-NMR Spectra of Cell Culture Media
Proton-decoupled 2H-NMR spectra were recorded for all of the withdrawn cell media samples for HUH-7 and AML12 cell lines. HDO signal intensity increases with incubation time period whereas [2H7]glucose decreases rapidly for HUH-7 cells as indicated by stacked plot of 2H-NMR spectra in
Longitudinal Relaxation Time (T1) Measurement of the Compounds
Once the lactate isotopomers were identified, it would be expected to observe differential changes in T1 due to a loss of 1H dipolar relaxation. T1 inversion recovery experiments were performed in the 2H-decoupling on mode to calculate the T1 relaxation time of the pyrazine and lactate isotopomers.
Quantification of the Lactate Isotopomers
Pyrazine was used as an internal standard to quantify each of the lactate isotopomers. Peak areas of the isotopomers were extracted from each of the {2H}1H—NMR spectra of the HUH-7 cells at different incubation times. After correction for T1 and for the number of Hs leading to each resonance, the concentration of the isotopomers was plotted (
Quantification of the HDO and Residual [2H7]Glucose
Pyrazine-D4 was used as an internal standard to quantify the HDO production and residual [2H7]glucose in each of the cell media samples of HUH-7 and AML12 cell lines. Peak areas of pyrazine-D4, HDO and residual [2H7]glucose was used to measure the concentrations of HDO and residual [2H7]glucose with respect to the 5 mM pyrazine-D4 internal standard. Cellular [2H7]glucose consumption is calculated from the residual [2H7]glucose in the cell media whereas HDO production was corrected for the natural abundance HDO in the cell media samples of each time points. Cavity ring-down spectroscopy (CRDS) measurements of starting water 2H content provided a secondary internal standard.
The 1H—NMR spectrum has a complicated pattern due to the very small geminal couplings (2JH-D) and vicinal couplings (3JH—H) of methyls of the partially deuterated lactate isotopomers (
Proton relaxation times of the lactate isotopomers increases incrementally with the number of deuterons (
Ben-Yoseph et. al. found a significant loss of 2H atoms from C1 and C6 of the [1-2H]glucose and [6,6-2H2]glucose respectively. Loss of 2H from C6 of glucose was mainly due to catalytic action of pyruvate kinase (PK) whereas loss of 2H from C1 of glucose was due to both PK and phosphomannose isomerase (PMI) (
Pentose phosphate pathway (PPP) is another possible source of diversion of 2H from C1 of [2H7]glucose derived intermediates (
Relative percent contributions of lactate isotopomers to the efflux pool size are shown in
The mechanism of each step of glycolysis, indicated that HDO would be produced by triose phosphate isomerase (TPI) from the C2 position of GA3P and DHAP, within the course of direct metabolism of one equivalent of [2H7]glucose to pyruvate. HDO production from the C1/C2 positions can occur due to the interconversion of F6P and M6P via phosphomannoseisomerase (PMI) (Chandramouli, V., et al. American Journal of Physiology-Endocrinology and Metabolism 277:E717-E723 (1999)). Enolase activity on phospho-enol pyruvate (PEP) precursor and keto-enol tautomerization in pyruvate are also prominent sites of HDO production during glycolysis (
The results indicate that metabolism of [2H7]glucose produces an increase in the HDO signal that correlates with cellular glucose uptake in the HUH7 cells, and a dramatic difference in HDO production has been observed between HUH-7 cancer cells and background precursor cells (AML12 cells) (
Warburg metabolism is phenotypic of many cancers. Increased lactate production is a consequence of increased glucose consumption and reduced pyruvate dehydrogenase flux in cancer cells. This study has successfully demonstrated the metabolic conversion of [2H7]glucose to lactate isotopomers and also quantitatively accounted for HDO production in these models. Conversion of DHAP and GA3P into lactate—CH2D and lactate-CHD2 can be measured efficiently by this method as well, which also indicates the 2H loss from [2H7]glucose during glycolysis. Quantitation of HDO production may provide a simplified means for metabolic imaging of the deuterium signal specific to cancer cell glucose uptake, analogous to [18F]-FDG PET, with superior signal-to-noise of other metabolite signals, as well as 2H MRS imaging.
Based on parameters for oral (Gunning, R. R. Jama 240:1611 (1978)) and intravenous administration of glucose in humans and estimates of tumor geometry and cellularity (Kondo, F. Hepatology International 3:283-293 (2009)), HDO enrichment in the peritumoral area may reach several percent, sufficient to observe lengthening of 1H T1 due to less efficient dipole-dipole relaxation (Sherry, A. D., et al. Analytical Biochemistry 52:415-420 (1973)). T1 weighted 1H imaging might therefore also be diagnostic, which would remove the need for purpose built 2H detection coils. A useful feature of our proposed tracer method versus methods based on TCA intermediate/lactate ratios is that it obviates the need to correct for 2H loss between the two species, reducing inter-subject variance and providing for easier extension to multiple tumor types.
From a practical point of view, {2H}1H—NMR experiments are easy to set up and require minimal time compared to direct 2H detection for 2H-lactate. Measuring the kinetics of the lactate isotopomer production from deuterated glucose by utilizing the more sensitive 1H—NMR experiment makes an easy method for measuring glycolytic rates. It is to be determined if this method could be translated for in vivo use. Lactate detection by 1H MRS is typically confounded by overlap of resonances derived from FAs. With a large FA signal the difference spectrum approach used here would have to be carefully implemented to prevent improper subtraction of the FA peaks and uncertainty in quantification of the 2H lactate. As compared to 13C based editing sequences, this method does not have to contend with a natural abundance correction. Also, the isotope shifts could potentially provide a more robust subtraction signal. The differential T1 of the lactate isotopomers is more difficult to take advantage of, as the T1 is lengthened upon 2H-enrichment. Standard T1 weighting would not provide contrast from the isotopomers. Specific choice of the delay time might allow the lactate—CH3 signal to be nulled which the 2H isotopomers remain using inversion recovery schemes.
In conclusion, this study demonstrates that administration of [2H7]glucose can trace Warburg metabolism in cancer cells using two readily applicable isotopomer detection methods: evolution of HDO, detected by 2H—MR, and of three different isotopomers of lactate, detected by {2H}1H—NMR with difference spectroscopy to remove interference from fatty acids. As compared to earlier results with [6,6-2H2]glucose, this method retains the selectivity of lactate enrichment, but with the increased likelihood of HDO generation, amplifies the sensitivity of the experiment to increased glycolytic flux. Either technique can be implemented on preclinical or clinical hardware and may provide the basis for a safe, zero-dose method of imaging cancer glucose uptake and Warburg metabolism.
Introduction
Glucose is the principal source of energy in the brain, and its dynamic uptake and utilization remain topical after decades of research. Metabolic imaging using MR based methods has traditionally been limited to 1H based spectroscopic detection of metabolites via chemical shift imaging (CSI) or single voxel spectroscopy (Frahm J, et al. Journal of Magnetic Resonance (1969) 1985 64:81-93; Bottomley P A. Ann N Y Acad Sci 1987 508:333-348). While these methods have been of immense use in clinical research, sensitivity issues have made their translation to the clinic problematic. More recently, hyperpolarized (HP)13C imaging has been used to study heart, liver, and brain metabolism, as well as changes in metabolism associated with various cancers (Merritt M E, et al. Proceedings of the National Academy of Sciences 2007 104:19773-19777; Salamanca-Cardona L, et al. Cancer Metab 2015 3:9). HP imaging uses the large gain in signal-to-noise ratio (SNR) made possible by dynamic nuclear polarization (DNP) to produce chemically selective images that reflect changes in metabolic flux (Comment A, et al. Biochemistry 2014 53:7333-7357). However, DNP is limited to molecules having long relaxation times amenable to sample transfer from the polarizer to the patient. Many important metabolites such as glucose have short relaxation times, making clinical translation more challenging (Kurhanewicz J, et al. Neoplasia 2019 21:1-16). 2-[18F]-fluorodeoxyglucose positron emission tomography (FDG-PET) is widely utilized for cancer diagnosis, staging, and treatment assessment (Kelloff G J. Clinical Cancer Research 2005 11:2785-2808). However, the specificity of FDG-PET is limited in the brain, where normally high glucose uptake and non-cancerous inflammation can give false positive results (Shreve P D, et al. Radiographics 1999 19:61-77). Additionally, FDG-PET is limited to glucose uptake and phosphorylation; glycolytic and TCA cycle metabolism are invisible to this technique.
Another approach that utilizes the chemical selectivity of MR utilizes 2H labelled substrates to image metabolic turnover. In vivo deuterium magnetic resonance spectroscopy (DMRS) has been demonstrated using [6,6-2H2]glucose to measure glucose uptake and glutamate/glutamine (glx) production over time, providing an alternative method for turnover studies in the brain (Lu M, et al. J Cereb Blood Flow Metab 2017 37:3518-3530). The same glucose isotopomer was used to image human glioblastoma (De Feyter H M, et al. Sci. Adv. 2018 4:eaat7314). However, the deuterated glucose isotopomer chosen for the study predetermined that downstream deuteration at lactate-C3 and glx-C4 would necessarily be the metabolic readouts of glucose metabolism. These peaks are slowly enriched as a function of time, and are most easily imaged using the CSI technique, as they are smaller peaks relative to the deuterated water (HDO) signal that arises both from natural abundance contributions and from a limited set of exchange reactions that eliminate the 2H label from the methyl position of lactate. As disclosed herein, [2H7]glucose is an alternative agent, and HDO evolution accurately reflects total glucose uptake (Mahar R, et al. Sci Rep 2020 10:8885), in analogy to FDG-PET (Vlassenko A G, et al. Neurolmage 2006 33:1036-1041), but with detection of glucose downstream metabolism and without the use of ionizing radiation.
This example demonstrates the application of this method in vivo to image glucose consumption in the rat brain. Administration of [2H7]glucose via the tail vein produced deuterated lactate and glx analogous to earlier reports, but shown here for the first time is that HDO generated from metabolism can be imaged in the rodent brain using gradient echo based methods as opposed to CSI. Using methods other than CSI to detect metabolic flux fundamentally transforms the technique, allowing excellent in plane resolution and faster acquisition times. By using HDO evolution as a metric of glucose utilization in the brain, a biomarker of cerebral glycolysis and oxidative flux is introduced that is imaged with superior sensitivity and resolution compared to previous approaches.
Methods
Chemicals and Tracer
[1,2,3,4,5,6,6-2H7]-D-glucose (or [2H7]glucose) was purchased from Sigma Aldrich, USA. Isoflurane anesthesia was purchased from Aspen Veterinary Resources, M O, USA. Sterile filtered saline solution was prepared to dissolve the [2H7]glucose solution. Heparin at 1000 USP/ml of glucose solution was also utilized to stop the clotting of blood during animal preparation.
Animal Experimental Setup
All experimental procedures were approved by the University of Florida Institutional Animal Care and Use Committee. Female Sprague Dawley rats (body weight: ˜220 g) were purchased from Jackson Laboratories and acclimatized in the University's animal care facilities. All rats were in the fed state. The rat was anesthetized with 2% isoflurane in a mixture of O2 and N2 gases. The tail vein was catheterized using a 23-gauge needle for the infusion of [2H7]glucose tracer. The body temperature of the rat was maintained at 37° C. using a heated circulating water pad. The respiratory rate was monitored during the entire period of the experiment. After shimming the magnet, three baseline 2H spectra were acquired. Subsequently, a 2 minute tail vain infusion of [2H7]glucose (1.95 g/kg body weight and dissolved in 2.3 mL saline) was administered. Blood was drawn after the end of spectroscopy and imaging experiments via saphenous vein and blood plasma was separated by centrifugation of the blood at 4° C. for 10 minutes. Circulating blood glucose was assessed by a blood draw and glucometer (Evancare, Medline Industries, Northfield, IL) prior to infusion and 3 hours after the glucose injection.
MRI System
Animal experiments were performed on an 11.1 T Magnex magnet (Magnex Scientific Ltd.) interfaced to a BrukerAvance Ill HD console running on ParaVision 6.0.1 (Bruker Instruments, Billerica, Mass.). The MRI system was equipped with RRI BFG-240/120-S6 gradient coil with the bore size of 120 mm and capable of producing 100 mT/m with 200 μs rise time. A shim map was generated for localized shimming on rat brain using an 85 mm home-built, actively-decoupled, linear volume 1H transmit-receive coil. The coil was tuned to the 1H NMR frequency at 470 MHz. A home-built 14 mm diameter, circular surface coil tuned to 72.26 MHz was used to acquire 2H spectra and images.
Deuterium Magnetic Resonance Spectroscopy
The 2H coil was placed over the rat head and a 3D printed phantom containing approximately 200 μl of natural abundance water was placed on the top of the coil to confirm positioning. The 2H RF pulse and flip angle were optimized to maximize the 2H signals from the rat brain. A partial signal contribution from the natural abundance 2H of the water phantom was unavoidable but was accounted for by subtraction of the signal areas from runs that followed the glucose injection. A single pulse acquire sequence was used to obtain 2H NMR spectra, with an observed linewidth of approximately 40 Hz. Spectra were acquired with a 1 kHz width and 256 data points for each free induction decay (FID), a repetition time of 300 ms, and a flip angle of 60° with 256 averages. The total acquisition time for each spectrum was 1.28 minutes. A total of 30 spectra were acquired for a 38-minute time period. Two additional 2H NMR spectra from the rat brain were acquired at 2 hours and 2.5 hours following [2H7]glucose infusion. Pyrazine-D4 (2.5 mM, internal standard) was added in each plasma sample to calculate the HDO concentration in plasma. 2H NMR spectra were collected on Bruker 14 T NMR magnet system equipped with a 1.7 mm TCl CryoProbe and Avance Neo console (Bruker BioSpin). 2H NMR spectra of plasma samples were acquired at 92.13 MHz resonant frequency with an acquisition time of 1.0 s, relaxation delay of 2 s, and 3072 scans.
Deuterium Magnetic Resonance Imaging
Axial 1H images of rat brain were acquired using a Bruker Fast Low Angle Shot (FLASH) gradient echo sequence. The two-dimensional (2D) acquisition matrix size of 128×128, field-of-view (FOV) of 50×50 mm2, and 8 averages were used to acquire 10 slices of rat brain with 2 mm thickness. Before [2H7]glucose infusion, a natural abundance deuterium 2H image of the rat brain was acquired using the FLASH sequence, with the 2H water signal on resonance. Two deuterium images were also acquired at 38 minutes and ˜2 hours post-[2H7]glucose infusion. The following parameters were used to acquire an axial 2H images of the rat brain: 2D acquisition matrix size 32×32, FOV of 50×50 mm2, slice thickness of 14 mm, repetition time of 100 ms, 256 signal averages, 1.416 ms echo time (TE), and 34.2° flip angle with an in-plane resolution of 1.56 mm2. Each image took 13 minutes to acquire. The Bruker multi-gradient echo (MGE) imaging experiment was utilized to acquire the echo images for separating HDO-only and [2H7]glucose-only images in the rat brain during the pseudo-steady state. Axial two-echo images were acquired with echo times of 11.1 ms (GE1) and 16.6 ms (GE2), acquisition matrix size of 32×32, FOV of 40×40 mm 2, and slice thickness of 14 mm. The in-plane resolution was 1.25 mm2. Proton intensity images of the rat brain was also acquired with the same FOV as deuterium images for co-registration purposes.
NMR and MRI Data Processing
2H NMR data analysis was performed using MestReNova v14.0.1-23284 (Mestrelab Research S. L., Santiago, Spain). For processing of the 2H NMR, each spectrum was Fourier transformed with an exponential window function of 2.0 Hz, zero-filling of 512 data points and automatic phase correction was also applied. Peak areas were extracted for HDO, [2H7]glucose including the anomeric signal, deuterated glx (glutamate/glutamine) and lactate from the 2H NMR spectra of the rat brain. 2D 1H and 2H images were processed in the ImageJ software (NIH, Bethesda, Md.). Since the acquired 2H image size was 32×32 and the 1H brain image was 128×128, the deuterium image was scaled using bicubic interpolation to 128×128 for co-registration and overlay of the images. This resulted in a final, interpolated, in plane resolution of 0.39 mm2. MGE echo images were reconstructed to 128×128 matrix size with gaussian windowing filter using Bruker ParaVision 6.0.1 software. Reconstructed echo images were imported to MATLAB for further processing to separate the contribution of HDO and [2H7]glucose to the rat brain DMRI. The Dixon reconstruction algorithm (Ragan D K, et al. J. Magn. Reson. Imaging 2010 31:510-514) was applied to reconstruct the HDO-only and [2H7]glucose-only images from the MGE images. T2 relaxation times of HDO and [2H7]glucose were also used as part of the Dixon reconstruction (De Feyter H M, et al. Sci Adv 2018 4:eaat7314).
Results
The initial blood glucose was measured as 8.0±1.4 mM (n=4) prior to infusion and the final blood glucose concentration at 3 hours after injection was 7.1±1.8 mM (n=4). The natural abundance HDO water signal (
Two gradient echo images were acquired at approximately 38 minutes and 2 hours post infusion. The initial image prior to injection shows a very faint background HDO signal, which increases significantly with time (
To distinguish the HDO signal from the [2H7]glucose substrate, an MGE based two-point Dixon method was implemented during the initial time period (approximately 10-24 minutes) of the post-[2H7]glucose infusion and metabolism (
Discussion
The brain image of
As detailed in the results, the first 2H image after glucose injection (
While the images here were acquired with a significant slice thickness that must cause significant partial volume effects, the overall SNR indicates that thinner slices will be readily obtainable once the protocol and hardware are further optimized. The fact that the HDO concentration in the brain remains below that of body water enrichment at 3 hours indicates that influx of HDO produced outside the brain is not as rapid as might be expected. This is confirmed by the HDO appearance at the longer time points (
Metabolic imaging of the brain has been a research target driven by both pathophysiological considerations, like stroke and cancer, and by studies of brain structure and cognition as a basic science. Glucose is the primary source of energy in the brain, but the dynamic response of glucose utilization remains topical. In cancer, aerobic glycolysis is a common metabolic phenotype of many brain malignancies. This change in metabolic programming, often termed Warburg metabolism (Warburg O. Science 1956 123:309-314), is effectively targeted using FDG-PET, which measures excessive glucose uptake but does not otherwise interrogate downstream glucose metabolism. FDG-PET has been used to study glucose metabolism in the brain for forty years (Phelps M E, et al. Ann Neurol. 1979 6:371-388; Villien M, et al. Neurolmage 2014 100:192-199), and has been used to measure functional responses to stimuli (Vlassenko A G, et al. Neurolmage 2006 33:1036-1041). However, total radiation burden makes PET unsuitable in pediatric populations and prevents its serial use in patients (Chawla S C, et al. Pediatr Radiol 2010 40:681-686). In cognitive psychology, tremendous effort has been expended to map the effects of brain activation using the BOLD methodology (fMRI) (Yamada H, t al. NeuroReport 1997 8:3775-3778; Ugurbil K, et al. Annu. Rev. Biomed. Eng. 2000 2:633-660). Similar techniques measuring cerebral metabolic rate of O2 consumption (CMRO2) have been used to estimate glucose oxidation as well (Rodgers Z B, et al. J. Cereb. Blood Flow Metab. 2016 36:1165-1185). All these methods indirectly infer glucose metabolism without measuring the downstream metabolites. Hyperpolarized 13C does offer the ability to directly detect lactate concentrations and lactate dehydrogenase flux when using HP pyruvate as the tracer (Kurhanewicz J, et al. Neoplasia 2019 21:1-16; Merritt M E, et al. Proceedings of the National Academy of Sciences 2007 104:19773-19777), but this technique demands a special pre-polarizer to produce the imaging agent. Furthermore, hyperpolarized glucose is obviously an analogous method, but the short carbon-13 T1 of glucose makes it difficult to translate to human studies (Kurhanewicz J, et al. Neoplasia 2019 21:1-16). In contrast, imaging of deuterated agents requires only the addition of a detection coil tuned to the appropriate frequency, and therefore could be easily disseminated to clinical sites. As shown in this Example, by imaging HDO, an end-product of [2H7]glucose metabolism, one can produce an image of glucose disposal. This modality should be complimentary to 1H based methods, and kinetic analysis should yield absolute rates of glycolysis, and TCA cycle turnover in healthy tissues.
Metabolic imaging of HDO can therefore be used as a biomarker of [2H7]glucose, lactate, and glx metabolism. Spectrally selective imaging of HDO and glucose using a two-point Dixon method was accomplished with a simple multi-gradient echo sequence. This modality facilitates higher in-plane resolution and faster acquisition of imaging data than possible with CSI methods.
Nonalcoholic steatohepatitis (NASH) affects about one-third of all patients with obesity and type 2 diabetes mellitus (T2DM) becoming a major risk factor for the development of cirrhosis and hepatocellular carcinoma (HCC) (Anstee, Q M, et al. Nature reviews Gastroenterology & hepatology 2019 16(7):411-428). NASH is the second leading cause of liver transplantation in the USA. The transition from nonalcoholic fatty liver disease (NAFLD) (considered a rather “benign” condition that affects 70% of obese insulin-resistant patients with T2DM, not leading to cirrhosis) to NASH (a stage of hepatocyte necrosis, inflammation, and often active fibrogenesis leading to cirrhosis) which involves induction of abnormal fatty acid p-oxidation, leading to chronic inflammation and liver fibrosis (Younossi, Z, et al. Nature reviews Gastroenterology & Hepatology 2018 15(1):11). Current plasma biomarkers or metabolic/other imaging techniques are needed that can identify NAFLD and end-stage disease (i.e., advanced fibrosis or cirrhosis) but not early disease (i.e., NASH) when the disease is at a stage that demands diagnosis and therapeutic intervention. Therefore, staging of the disease is central to the effective treatment.
NAFLD is becoming a nearly ubiquitous consequence of the western diet, with the most recent estimates suggesting a prevalence as high as 55-65% in patients with obesity or T2DM in the United States. As the disease progresses, NASH develops by mechanisms that remain unclear in about one-third of patients and in many it results in cirrhosis. Within the next 5 years, NASH will become the number one cause of liver transplant in the United States. NASH not only can lead to cirrhosis, but is also associated with a higher incidence of hepatocellular carcinoma (HCC). There has been tremendous effort expended on identification of biomarkers and imaging (FibroScan or MR elastography) for the diagnosis of fibrosis in patients with NASH, but both of these modalities are insensitive to fatty acid oxidation (increased or decreased) and inflammatory response, two primary characteristics of NASH onset. Ideally, one would like to make the diagnosis of NASH before the development of liver fibrosis but the available biomarkers and imaging have not been good enough for widespread use. The only definitive method for the diagnosis of NASH is liver biopsy, which is expensive, sampling dependent, poorly accepted by patients and clinicians, and in a small percentage associated with complications, rarely even death. Development of an imaging method which could directly assess global hepatic metabolism and monitor patients for the NAFLD to NASH transition would transform our ability to manage this increasing prevalent, costly disease.
In mice, the transition from NAFLD to NASH is accompanied by an upregulation of fatty acid (FA) oxidation in the tricarboxylic acid (TCA) cycle, in seeming contradiction to the visible deposition of fat in the liver. Upregulated fatty acid oxidation in NASH has also been implicated in human studies. Increased FA oxidation leads to the accumulation of lipotoxic intermediates as well as increased oxidative stress in the liver. Previously carbon-13 magnetic resonance (MR) has been used to interrogate TCA cycle turnover, but only at tremendous expense (˜$25,000/patient in [1-13C]acetate label alone). Recently, it was demonstrated that metabolic imaging with deuterated glucose and acetate, and 2H MRI is feasible in human brain (De Feyter, et.al., Science, 2018). The low sensitivity of the 2H nucleus can be overcome by fast repetition of the experiment, as the MR T1 of deuterium is very short. Octanoate is a medium chain fatty acid that does not use CPT-1 for transport into the mitochondria. Therefore it is avidly and rapidly p-oxidized by the liver. Therefore, 2H MRI in patients after a bolus of perdeuterated octanoate can produce a metabolically sensitive image suitable for identifying the NAFLD to NASH transition.
A robust imaging method for detecting NASH development could transform the field and radically modify current patient management. Early intervention in patients with NASH, has been shown to prevent or reverse hepatocyte necrosis and liver inflammation/fibrosis (lifestyle targeted to high-risk patients, bariatric surgery vitamin E and pioglitazone; AASLD 2018 guidelines). Early and accurate staging can allow the physician to determine the moment when patient care should be modified, reducing costs, and improving patient outcomes.
Unless defined otherwise, all technical and scientific terms used herein have the same meanings as commonly understood by one of skill in the art to which the disclosed invention belongs. Publications cited herein and the materials for which they are cited are specifically incorporated by reference.
Those skilled in the art will recognize, or be able to ascertain using no more than routine experimentation, many equivalents to the specific embodiments of the invention described herein. Such equivalents are intended to be encompassed by the following claims.
This application claims benefit of U.S. Provisional Application No. 62/992,726, filed Mar. 20, 2020, U.S. Provisional Application No. 63/004,145, filed Apr. 2, 2020, and U.S. Provisional Application No. 63/120,996, filed Dec. 3, 2020, which are hereby incorporated herein by reference in their entireties.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US2021/012400 | 1/7/2021 | WO |
Number | Date | Country | |
---|---|---|---|
62992726 | Mar 2020 | US | |
63004145 | Apr 2020 | US | |
63120996 | Dec 2020 | US |