Detection of nucleic acids and other low abundance components in native bone and osteosarcoma extracellular matrix by isotope enrichment and DNP-enhanced NMR

Abstract

Sensitivity enhancement by isotope enrichment and DNP NMR enables detection of minor but biologically relevant species in native intact bone, including nucleic acids, choline from phospholipid headgroups, and histidinyl and hydroxylysyl groups. Labelled matrix from the aggressive osteosarcoma K7M2 cell line confirms the assignments of nucleic acid signals arising from purine, pyrimidine, ribose, and deoxyribose species. Detection of these species is an important and necessary step in elucidating the atomic level structural basis of their functions in intact tissue.

---

The extracellular matrix (ECM) is a complex, heterogeneous, and active, environment, playing structural roles, but also strongly influencing cell behaviour in both homeostasis and pathology. The difficulty of understanding the composition and regulation of the ECM parallels the situation within cells, but with the added complication that many molecular components are insoluble. Solid-state NMR (ssNMR) can contribute to a better understanding of the ECM through direct detection of sparse species, made possible by a strategy of isotopic enrichment of native biomaterials¹ combined with rapidly evolving dynamic nuclear polarization (DNP) enhancement techniques.² Briefly these are based on the transfer of the much greater electron spin polarization, excited by a microwave laser at the electron resonance frequency of an organic free radical dopant, to the (NMR active) nuclei of the molecules of interest, at low temperatures to maximize polarization. The conversion of electron to nuclear polarization results in NMR sensitivity enhancements of up to several orders of magnitude. Recently we have reported direct observation of important low abundance collagen post-translational modifications (PTMs) in intact native, and model, ECM.³ While we focussed on species derived from one specific amino acid (Lys) component of the ECM, it is also of great potential importance to investigate whether other sparse ECM species become observable with DNP enhancement. Non-protein minor species, besides low abundance amino acids, are of particular interest since they represent a realm in which NMR can provide atomic length scale molecular structural information complementing other current biophysical techniques. Incorporating non-protein molecules into the overall structural biology picture will broaden conceptions of which ECM components to investigate in further detail towards an integrated model of the biological and material properties of tissues.

Here, we report DNP-enhanced NMR on mouse bone¹ enriched in a wide range of ¹³C, ¹⁵N-amino acids. Where necessary to distinguish between low abundance protein, and non-protein, signals, we confirmed amino acid type assignments by comparison with ECM produced by cells cultured on media supplemented with one or a few specific amino acids, using methods already described.^4,5 Our previous analysis of bone focussed on the relatively abundant species observable with conventional ssNMR. The enhancement provided by DNP now reveals non-protein biomolecules biosynthesized from isotopically enriched amino acids and/or sugars, and their incorporation into bone. The preparation of stable isotope enriched biomaterials, and our DNP NMR methodology, is described in detail in the ESI,† Materials and methods sections.

We begin our investigation on bone with a suite of ssNMR characterizations under DNP signal enhancement conditions. A full 1D ¹⁵N DNP NMR spectrum is shown in ESI Fig. S1.† Apart from the predictable major signals (from backbone and sidechain amides, and abundant Lys and Arg) the sensitivity enhancement of enrichment-DNP reveals a number of more minor resonances (Fig. 1A) from lower abundance species; the shifts of many of those at high spectral frequency coincide closely with the common (deoxy)nucleotide base ring nitrogens (Biological Magnetic Resonance Database, BMRB, as summarized in ESI Schemes S1† (DNA) and S2 (RNA)). Signals from the nitrogen atoms of protonated His imidazolium (ca. 180 ppm) and unprotonated imidazole (ca. 250 ppm) ring forms are also observed.^6,7 The His assignments are reinforced by ¹⁵N data from model osteoblast (bone precursor) ECM selectively enriched in (U-¹³C₆, ¹⁵N₃)-His (His*) (Fig. 1B; see ESI† for preparation and comparison of ¹⁵N–¹³C correlation spectra, Fig. S2†). The low frequency region of the bone DARR-assisted ¹⁵N–¹³C correlation spectrum (Fig. 1C) displays a clear cross peak (65 ppm ¹³C, 52 ppm ¹⁵N) which we ascribe to the (CH₃)₃N̲⁺–C̲H₂–CH₂–O– correlation in the head groups of phosphatidylcholine (PC) lipids (¹³C assignments based on Kroon et al.;⁸ see ESI Fig. S3† for overlays with model compound spectra). The PC ethanolamine substructure is biosynthesized from a single serine unit so, in our bone material, all atoms will be concurrently labelled. Finally Fig. 1C also shows the distinctive Hyl N_ζ–C_δ fingerprint correlations; the close correspondence with (U-¹³C₆, ¹⁵N₂)-Lys enriched ECM from vascular smooth muscle cells³ is shown in ESI Fig. S4,† proving that our labelling – DNP NMR combination readily detects the Hyl PTM even in this complex organic–inorganic composite native biomaterial, as well as in simpler native (skin) and model ECMs.³ Chemical shifts are summarized in Table 1.

Fig. 1 (A) High frequency region of the ¹⁵N DNP-NMR spectrum showing assignments of resolved base signals of nucleic acids, and ring His signals; (B) ¹⁵N spectrum of His* FSOb ECM highlighting correspondence with the signal from the imidazole form of the His ring; (C) DARR assisted ¹⁵N–¹³C correlation spectrum of labelled bone, transmitting magnetization along the Hyl (and Lys) sidechain to emphasize correlations between Hyl (as well as the intense Lys) N_ζ and C_δ. Also highlighted is the signal from the (CH₃)₃N̲⁺–C̲H₂–CH₂–O– correlation from phosphatidylcholine headgroups. Atom numbering of the nucleic acid purine (top) and pyrimidine skeletons is also shown.

Table 1 Chemical shifts of resolved signals from low abundance species detected in ¹³C, ¹⁵N enriched bone by DNP NMR

a Protonated imidazolium form.b Unprotonated imidazole form.c Single unresolved signal.
Amino acids/PTMs
Hyl	Cδ	Cε	Nζ
	69	44	31
His	Cγ	Cδ2	Cε1	Nδ1^a	Nε2^a	Nδ1/ε2^b
	129	117	135	180	177	250
(Deoxy)nucleotide bases
		N1	N3	N7	N9
Ade		224	212	235^c	174^c
Gua		148		235^c	174^c
Cyt		148	197		197
Thy/Ura		162	148
Lipids
PC	–N̲(CH3)3^+	–C̲H2–N^+
	52	65

---

Detection of His, lipid components, and Hyl, in native biomaterials is not unexpected. In contrast nucleic acids in bone have been little discussed, although this tissue is a favourable source of forensic genetic material.⁹ Given this underrepresentation of bone in nucleic acid literature, we turned to ECM from the aggressively metastasizing K7M2 osteosarcoma cell line¹⁰ to substantiate these signals. A ¹³C DNP NMR spectrum of K7M2 ECM supplemented with (U-¹³C₂, ¹⁵N)-Gly (Gly*) and (U-¹³C₆)-glucose (Glc*) is shown in Fig. 2A; it shows a number of clear signals with shifts consistent with the common (deoxy)nucleotide bases¹¹ (ESI Schemes S1 and S2†). Abundant nucleic acid in the K7M2 ECM is confirmed by conventional in situ propidium iodide staining and confocal fluorescence microscopy (ESI Fig. S5†). The biosynthetic routes to the bases (ESI Scheme S3†) dictates that all the base carbon atoms will be at least partially ¹³C enriched by Gly* and Glc*, although only the purine C₄–C₅–N₇ substructure (from intact Gly*), and to some extent pyrimidine C₄–C₅–C₆ from Glc* via Ser,¹² will necessarily be concurrently labelled. (Deoxy)riboses are biosynthesized from Glc* and so all (deoxy)sugar rings will be concurrently ¹³C enriched. Fig. 2B details results of a through-bond (J-transmitted) INADEQUATE,¹³ which clearly shows the predicted purine C₄–C₅ and pyrimidine C₄–C₅ cross peaks. DARR¹⁴ (Fig. 2C) confirms these and indicates pyrimidine C₆ signals in Cyt and Thy, and Ura from RNA. The presence of both DNA dRib, and RNA Rib, is confirmed by the low frequency part of the INADEQUATE spectrum (Fig. 2D). ¹⁵N DNP NMR shows a clear signal from purine N₇'s, corresponding closely to that of the putative purine N₇'s in bone (Fig. 2E). Also a ¹³C{¹⁵N} TEDOR¹⁵ (Fig. 2F) mapping short carbon-nitrogen distances, shows the expected purine N₇–C₅ and weaker N₇–C₄ correlations, but also purine N₇–C₈ cross peaks (the latter atom acquired from Glc* via N-formyl tetrahydrofolate). All DNA and RNA shifts measured in K7M2 ECM are summarized in Table 2.

Fig. 2 ¹³C DNP NMR spectra of Gly*, Glc* supplemented K7M2 ECM. (A) ¹³C CP-MAS; nucleotide base signals assigned by 2D experiments and database shift values (ESI Schemes S1 and S2†) are labelled; (B) high frequency region of INADEQUATE J-transmitted SQ-DQ correlation spectrum (265 ppm ¹³C DQ sweep width to avoid aliasing; ^#, * – correlations between lipid olefinic carbons); (C) low frequency region of the same spectrum tracing out correlations within ribose (Rib), and deoxyribose (dRib), rings; (D) high frequency region of the ¹³C–¹³C DARR correlation spectrum base; (E) ¹⁵N DNP NMR of Gly*Glc* K7M2 ECM compared with that of labelled bone, supporting the assignment of the 235 ppm signal as purine N₇; (F) ¹³C{¹⁵N} TEDOR spectrum showing the expected purine N₇–C₅ and N₇–C₄, as well as the putative purine ring N₇–C₈, through-space correlations (¹⁵N sweep width 438 ppm to cover full ¹⁵N range).

Table 2 Chemical shifts of resolved (deoxy)nucleotide signals in the K7M2 rat osteosarcoma cell line

	C2	C4	C5	C6	C8	N7
a Signals observed in 1D spectrum although no SQ-DQ or DARR correlations observable; assignments are inferred from average shifts reported in BMRB.b Assignments to Ade and Gua are based on rank order of BMRB average shifts, and may be reversed.c From J-transmitted INADEQUATE SQ-DQ correlation.d From ^13C–^13C DARR correlation.e From ^15N–^13C TEDOR.
Ade	153^a	147^b^,^c	116^b^,^c	153^a	139^b^,^e	236^b^,^e
Gua	153^a	150^b^,^c	114^b^,^c	167^a	136^b^,^e	239^b^,^e
Cyt	153^a	165^c	96^c	140^d
Thy	153^a	167^c	110^c	137^d
Ura			105^d	145^d

---

	C1′	C2′	C3′	C4′	C5′
dRib	82	37	78	84	65
Rib	90	72	69	79	61

---

There are large number of catalogued nucleic acid shifts in the BMRB but as yet no robust equivalent of the chemical shift index¹⁶ of protein secondary structural elements, especially with respect to ¹⁵N. Nucleotide base ¹³C and ¹⁵N shifts are particularly sensitive to secondary structural features^17–20 but linewidths and overlap in our ¹³C and ¹⁵N spectra currently preclude structural conclusions about K7M2 ECM. In principle though this work suggests that the potential of ¹⁵N shifts as indices of nucleic acid structure^20,21 in whole tissues can be realized through several conformational “reporter” nitrogen atoms in the well resolved high frequency NMR window above 150 ppm. Nucleic acid in bone may be predominantly intracellular, unlike that of K7M2 where a higher extracellular proportion is expected, but it is interesting that these distinctive ¹⁵N signals can be observed in both systems. We envisage that detection by DNP-NMR in native tissue is an essential step in fulfilling this potential to characterize nucleic acid structure in the intact tissue and ECM environment.

NMR detection of nucleic acids in bone raises fresh questions about their biological role, and avenues to new insights. NMR detection is an essential first step in elucidating these roles at the molecular and atomic level in native materials. Extracellular roles e.g. signalling, of nucleic acids are likely to operate in K7M2 ECM, and potentially also in bone and bone diseases e.g. osteoporosis.²² Bone undergoes constant remodelling so a high degree of transcriptional, translational, and apoptotic activities are expected,²³ leading to significant amounts of nucleic acid being produced in bone cells and released in the ECM. Once in the extracellular environment, some nucleic acids may have a different function in addition to their intracellular roles. For instance, nucleic acids have pronounced affinity for apatite mineral,²⁴ and may participate in biomineralization,²⁵ as calcium sequestering and mineral templating scaffolds, a role that is already proposed for structurally related polynucleotides like poly-ADP-ribose.⁴

Phospholipids²⁶ are abundant in mineralized bone.²⁷ They may represent the presence of matrix vesicles – extracellular structures ca. 100 nm in size, “pinched off” from the outer membranes of osteoblasts – which act as a nidus for initial calcium phosphate precipitation in normal²⁸ and pathological²⁹ biomineralization. The molecular mechanisms underlying the initial bioprecipitation events and processes remain to be elucidated, in particular the role of matrix vesicle phospholipids. Direct observation of the PC headgroup by NMR in healthy bone is a significant step in atomic level understanding of the relevant intermolecular and interionic interactions. In addition to matrix vesicles, of 40–100 nm sized exosomes are the basis of a pathway for cancer metastasis³⁰ and therapy resistance³¹ through delivery of proteins, lipids, RNA and DNA to remote cells.³² The cargo carried by these vesicles may play unexpected roles if released directly in the extracellular space and allowed to interact with other components of the ECM.³³ DNP-enhanced NMR is potential non-targeted method of better understanding the role played by these extracellular vesicles.

Hyl and its subsequent modifications are essential to the structure, properties, and functions, of collagen.³⁴ The type I collagen heterotrimer contains 12 His residues. Both of those in each α1 (I) chain are two residues from those Lys residues which are sequence specifically δ-hydroxylated to Hyl, in a strongly conserved Lys-Gly-His-Arg tetrad.³⁵ These His residues have been proposed to play a part in targeting lysyl hydroxylases and thereby the formation of the specific Lys and Hyl crosslinks³⁶ essential to maintaining collagen and tissue structure. Some of the same His residues are not only hydroxylated to Hyl, but of these some are subsequently also enzymatically glycosylated to (Glc)GalHyl, an essential event in directing the as-yet poorly understood process of collagen fibril assembly. It has recently been proposed that these self same residues are also preferentially glycated in non-enzymatic, usually harmful, reactions with metabolic sugars.³⁷ Finally there is evidence for a trivalent crosslink, histidinohydroxylysinonorleucine, between these same Lys/Hyl and His residues.³⁸ Thus His and Lys/Hyl residues co-exist in crosslinking and glycosylation/glycation “hotspots” which may be key to better understanding the spatial relationships underlying orderly, and pathological, collagen assembly and structure. Detecting His, Lys and Hyl simultaneously, as enabled by isotopic enrichment and DNP-enhanced NMR, in native and model ECM, opens fertile new possibilities for studying these processes.

This communication demonstrates again³ the ability of ssNMR, combined with signal enhancement by DNP and targeted isotope labelling, to detect low abundance components in intact tissue that may have surprising biological roles, thus complementing other biophysical and molecular biology techniques in exploring diverse areas of matrix biology.

Conflicts of interest

There are no conflicts of interest. Animal procedures were unregulated in terms of the U. K. Animals (Scientific Procedures) Act 1986, but nevertheless were reviewed for compliance with institutional ethical guidelines by the Animal Welfare Ethical Review Board, University of Cambridge.

Acknowledgements

U. K. Medical Research Council (DGR, RR, RG75828), the Engineering and Physical Sciences Research Council PhD studentship (IG, AP), CSC Cambridge Scholarship for PhD studentship (RL), postdoctoral fellowship from the DAAD and Leibniz Association (WYC), Raymond and Beverly Sackler Fund for Physics of Medicine, University of Cambridge (AP), and support by iNEXT, grant number 653706, funded by the Horizon 2020 programme of the European Commission.

Notes and references

1. V. W. C. Wong, D. G. Reid, W. Y. Chow, R. Rajan, M. Green, R. A. Brooks and M. J. Duer, J. Biomol. NMR, 2015, 63, 119.
2. A. S. L. Thankamony, J. J. Wittmann, M. Kaushik and B. Corzilius, Prog. Nucl. Magn. Reson. Spectrosc., 2017, 102–3, 120.
3. W. Y. Chow, R. Li, I. Goldberga, D. G. Reid, R. Rajan, J. Clark, H. Oschkinat, M. J. Duer, R. Hayward and C. M. Shanahan, Chem. Commun., 2018, 54, 12570.
4. W. Y. Chow, R. Rajan, K. H. Muller, D. G. Reid, J. N. Skepper, W. C. Wong, R. A. Brooks, M. Green, D. Bihan, R. W. Farndale, D. A. Slatter, C. M. Shanahan and M. J. Duer, Science, 2014, 344, 742.
5. R. Li, R. Rajan, W. C. V. Wong, D. G. Reid, M. J. Duer, V. J. Somovilla, N. Martinez-Saez, G. J. L. Bernardes, R. Hayward and C. M. Shanahan, Chem. Commun., 2017, 53, 13316.
6. S. H. Li and M. Hong, J. Am. Chem. Soc., 2011, 133, 1534.
7. F. Hu, K. Schmidt-Rohr and M. Hong, J. Am. Chem. Soc., 2012, 134, 3703.
8. P. A. Kroon, D. M. Quinn and E. H. Cordes, Biochemistry, 1982, 21, 2745.
9. T. Kitayama, Y. Ogawa, K. Fujii, H. Nakahara, N. Mizuno, K. Sekiguchi, K. Kasai, N. Yurino, T. Yokoi, Y. Fukuma, K. Yamamoto, T. Oki, H. Asamura and H. Fukushima, Leg. Med., 2010, 12, 84.
10. C. Khanna, J. Khan, P. Nguyen, J. Prehn, J. Caylor, C. Yeung, J. Trepel, P. Meltzer and L. Helman, Cancer Res., 2001, 61, 3750.
11. I. V. Sergeyev, L. A. Day, A. Goldbourt and A. E. McDermott, J. Am. Chem. Soc., 2011, 133, 20208.
12. D. L. Nelson and H. M. Cox, Lehninger: Biochemistry, W. H. Freeman & Co., New York, 6th edn, 2013.
13. A. Lesage, C. Auger, S. Caldarelli and L. Emsley, J. Am. Chem. Soc., 1997, 119, 7867.
14. K. Takegoshi, S. Nakamura and T. Terao, Chem. Phys. Lett., 2001, 344, 631.
15. C. P. Jaroniec, C. Filip and R. G. Griffin, J. Am. Chem. Soc., 2002, 124, 10728.
16. D. S. Wishart, Prog. Nucl. Magn. Reson. Spectrosc., 2011, 58, 62.
17. P. N. Borer, S. R. LaPlante, N. Zanatta and G. C. Levy, Nucleic Acids Res., 1988, 16, 2323.
18. S. R. LaPlante, E. A. Boudreau, N. Zanatta, G. C. Levy, P. N. Borer, J. Ashcroft and D. Cowburn, Biochemistry, 1988, 27, 7902.
19. B. L. Gaffney, B. Goswami and R. A. Jones, J. Am. Chem. Soc., 1993, 115, 12607.
20. X. P. Xu and S. C. F. Au-Yeung, J. Phys. Chem. B, 2000, 104, 5641.
21. X. Gao and R. A. Jones, J. Am. Chem. Soc., 1987, 109, 3169.
22. M. P. Yavropoulou and J. G. Yovos, J. Musculoskeletal Neuronal Interact., 2018, 18, 18.
23. R. L. Jilka, B. Noble and R. S. Weinstein, Bone, 2013, 54, 264.
24. O. M. Loreille, T. M. Diegoli, J. A. Irwin, M. D. Coble and T. J. Parsons, Forensic Sci. Int.: Genet., 2007, 1, 191.
25. R. Coscas, M. Bensussan, M. P. Jacob, L. Louedec, Z. Massy, J. Sadoine, M. Daudon, C. Chaussain, D. Bazin and J. B. Michel, Atherosclerosis, 2017, 259, 60.
26. A. L. Boskey and D. M. Timchak, Metab. Bone Dis. Relat. Res., 1983, 5, 81.
27. R. A. Haraszti, M. C. Didiot, E. Sapp, J. Leszyk, S. A. Shaffer, H. E. Rockwell, F. Gao, N. R. Narain, M. DiFiglia, M. A. Kiebish, N. Aronin and A. Khvorova, J. Extracell. Vesicles, 2016, 5, 32570.
28. E. E. Golub, Biochim. Biophys. Acta, 2009, 1790, 1592.
29. H. H. T. Hsu, N. Camacho and H. Aono, Atherosclerosis, 1999, 144, 94.
30. C. Williams, R. Rodriguez-Barrueco, J. M. Silva, W. J. Zhang, S. Hearn, O. Elemento, N. Paknejad, K. Manova-Todorova, K. Welte, J. Bromberg, H. Peinado and D. Lyden, Cell Res., 2014, 24, 766.
31. I. Li and B. Y. Nabet, Mol. Cancer, 2019, 18, 32.
32. T. Lagerweij, M. Perez-Lanzon and S. R. Baglio, J. Vis. Exp., 2018, 135, e56932.
33. L. Min, J. Shen, C. Q. Tu, F. Hornicek and Z. F. Duan, Cancer Metastasis Rev., 2016, 35, 377.
34. J. Myllyharju, Top. Curr. Chem., 2005, 247, 115.
35. F. Rodriguez-Pascual and D. A. Slatter, Sci. Rep., 2016, 6, 37374,  DOI:10.1038/srep37374.
36. D. L. Helseth, J. H. Lechner and A. Veis, Biopolymers, 1979, 18, 3005.
37. D. M. Hudson, M. Archer, K. B. King and D. R. Eyre, J. Biol. Chem., 2018, 293, 15620.
38. G. L. Mechanic, E. P. Katz, M. Henmi, C. Noyes and M. Yamauchi, Biochemistry, 1987, 26, 3500.

---

Footnote

† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9ra03198g

---