# Interfacial Ca2+ environments in nanocrystalline apatites revealed by dynamic nuclear polarization enhanced 43Ca NMR spectroscopy

Daniel Lee, César Leroy, Charlene Crevant, Laure Bonhomme-Coury, Florence Babonneau, Danielle Laurencin, Christian Bonhomme, Gaël de Paëpe

► To cite this version:

Daniel Lee, César Leroy, Charlene Crevant, Laure Bonhomme-Coury, Florence Babonneau, et al. Interfacial Ca2+ environments in nanocrystalline apatites revealed by dynamic nuclear polarization enhanced 43Ca NMR spectroscopy. Nature Communications, 2017, 8, pp.14104. 10.1038/ncomms14104. hal-01448516

HAL Id: hal-01448516

https://hal.science/hal-01448516

Submitted on 13 Feb 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Interfacial Ca²⁺ environments in nanocrystalline apatites revealed by dynamic nuclear polarization enhanced ⁴³Ca NMR spectroscopy

Daniel Lee<fnref n="1" />,², César Leroy³, Charlène Crevant³, Laure Bonhomme-Coury³, Florence Babonneau³, Danielle Laurencin⁴, Christian Bonhomme³ & Gaël De Paëpe<fnref n="1" />,²

The interfaces within bones, teeth and other hybrid biomaterials are of paramount importance but remain particularly difficult to characterize at the molecular level because both sensitive and selective techniques are mandatory. Here, it is demonstrated that unprecedented insights into calcium environments, for example the differentiation of surface and core species of hydroxyapatite nanoparticles, can be obtained using solid-state NMR, when combined with dynamic nuclear polarization. Although calcium represents an ideal NMR target here (and *de facto* for a large variety of calcium-derived materials), its stable NMR-active isotope, calcium-43, is a highly unreceptive probe. Using the sensitivity gains from dynamic nuclear polarization, not only could calcium-43 NMR spectra be obtained easily, but natural isotopic abundance 2D correlation experiments could be recorded for calcium-43 in short experimental time. This opens perspectives for the detailed study of interfaces in nanostructured materials of the highest biological interest as well as calcium-based nanosystems in general.

<footnote n="1">¹Univ. Grenoble Alpes, INAC, F-38000 Grenoble, France. ²CEA, INAC, F-38000 Grenoble, France. ³Sorbonne Universités, UPMC Univ Paris 06, Collège de France, UMR CNRS 7574, Laboratoire de Chimie de la Matière Condensée de Paris, 4 place Jussieu, 75252 Paris cedex 05, France. ⁴Institut Charles Gerhardt de Montpellier, UMR5253, CNRS UM ENSCM, CC1701, Place Eugène Bataillon, 34095 Montpellier cedex 05, France. Correspondence and requests for materials should be addressed to G.D.P. (email: gael.depaepe@cea.fr).</footnote>
Calcium is an important element in many natural and synthetic materials. A key calcium-containing material is apatite, a calcium phosphate that is naturally present in sedimentary rocks in the form of fluorinated francolite, as well as in bones and teeth of vertebrae in the form of carbonated hydroxyapatite (Ca₁₀(PO₄)₆(OH)₂). The broad diversity of structural roles played by calcium compounds undoubtedly comes from the ability of Ca²⁺ cations to adopt a large variety of coordination environments, with coordination numbers varying from 4 to 9, and many different coordination geometries. Thus, determining the local environment of calcium is of fundamental importance to fully understand the structure and properties of these materials. Unfortunately, very few techniques are available to probe the local structure around calcium atoms, especially in disordered or amorphous materials, or in materials with several different calcium environments, as is the case in organic-inorganic biomaterials like bone, where calcium is present not only in the mineral phase but also at the organic-mineral interface. In a first rough approximation, nanosized synthetic carbonated hydroxyapatite (C-HAp), obtained by precipitation in aqueous medium, can act as a structural analogue for the inorganic component of bone tissues. Such nano-objects present intrinsic distributions of both cationic and anionic sites in the core as well as at the surface of the particles, leading to challenging analysis and interpretation.

Although solid-state nuclear magnetic resonance (NMR) has been used to study Ca²⁺ environments¹⁻⁵, it has only rarely been applied to complex materials because of severe sensitivity drawbacks. Indeed, the natural abundance of ⁴³Ca, the NMR-active stable isotope of calcium, is only 0.14%. Moreover, it has a small gyromagnetic ratio (|γ<sub>Ca-43</sub>/γ<sub>H-1</sub> ~ 0.07) and spin-7/2, rendering it a highly insensitive nucleus. Thus, natural abundance experiments require large sample volumes (using 7 mm rotors or larger, limiting magic angle spinning (MAS) rates to ~5 kHz) in combination with the highest magnetic field available and signal-enhancing NMR pulse sequences²,⁵⁻¹⁰, and even in these favourable conditions, experimental times for recording ¹D ⁴³Ca MAS NMR spectra can reach ~7–10 h. It follows that attempting to gain further information from 2D experiments, recorded at natural abundance, remains unrealistic. In the case of ⁴³Ca-labelled samples, ⁴³Ca 2D and triple-resonance NMR experiments have been reported, showing that much greater insight into Ca²⁺ local environments can be accessed, but so far only a limited set of such experiments have been performed because the commercially-available ⁴³Ca-labelled starting chemical, ⁴³Ca-labelled calcite (CaCO₃), is highly expensive⁵,¹¹⁻¹⁵.

Recently, high-field dynamic nuclear polarization combined with magic angle spinning (MAS-DNP)¹⁶ has been shown to alleviate sensitivity limitations of NMR for both biological¹⁷ and material¹⁸ analytes and, accordingly, some of the difficulties when trying to observe low-gamma (|γ| ≤ |γ<sub>N-15</sub>|) nuclei¹⁹⁻²². This is because in MAS-DNP the large polarization of electron spins can be transferred to nuclei, greatly increasing their sensitivity in solid-state NMR experiments. Moreover, MAS-DNP experiments are generally performed with low sample temperatures of ~100 K, which has been shown to be beneficial not only for the efficiency of the DNP process but also for cross-polarization (CP) of quadrupolar²³ and low-gamma²⁰ nuclei. Recent work using MAS-DNP has also shown its potential for studies of quadrupolar nuclei with low natural abundance, such as ¹⁷O (refs 24,25). Nevertheless, fast ⁴³Ca NMR at natural isotopic abundance still seemed a pipe dream owing to the intrinsic combination of NMR obstacles.

Herein, it is not only demonstrated that MAS-DNP can be used to facilitate the fast acquisition of ⁴³Ca NMR spectra, but it is also shown that this can be done for samples at their natural isotopic

abundance, even allowing the successful acquisition of 2D correlation experiments, which provide here the unambiguous discrimination of core and surface calcium sites in hydroxyapatite nanoparticles.

## Results

**The C-HAp nanoparticles.** The synthetic protocol to obtain the C-HAp nanoparticles is given in the ‘Methods’ section and the routine characterization (transmission electron microscopy (TEM), X-ray powder diffraction (XRD), Fourier transform infrared spectroscopy (FTIR), thermal gravimetric analysis (TGA)) of the sample is presented in Supplementary Figs 1–4. Figure 1 displays a schematic representation of stoichiometric HAp. Carbonate groups can partially substitute OH⁻ (A-type) and/or PO₄³⁻ moieties (B-type)²⁶. The OH⁻ groups form columns (parallel to the c axis of the structure) that can be considered as rather isolated (the distance between two consecutive columns is ~0.9 nm) and can act as paths for quasi one-dimensional ¹H-¹H spin-diffusion. It is worth noting that Klimavicius *et al*.²⁷ demonstrated very recently that moderate MAS rates (up to 9 kHz) had little impact on the spin-diffusion in nanostructured HAp. Following the synthetic protocol, some Ca²⁺ cations were substituted by Na⁺ (see the chemical analysis in the Supplementary Methods).

**Preparation for DNP experiments.** Approximately 30 mg of the white nanoparticulate (diameter ~30 nm, see Supplementary Fig. 1) C-HAp powder was wetted with a 30 µl aliquot of a glass-forming solvent matrix of glycerol-d₈/D₂O/H₂O (60/30/10; v/v/v) containing 10 mM AMUPol²⁸ biradical polarizing agent to create a DNP-ready sample²⁹. The resulting slightly damp powder was packed into a 3.2 mm outer-diameter thin-wall zirconia MAS rotor and this was then inserted into the pre-cooled (100 K) MAS-DNP triple-resonance (¹H/¹³C/⁴³Ca) NMR probe and pneumatically spun to the desired rate using cold nitrogen gas. Under these MAS conditions, a few watts of microwave (µw) irradiation at approximately the Larmor frequency of the unpaired electron spins of AMUPol (~263 GHz) are applied to the sample, inducing a transfer of the large spin polarization of the electrons to nearby nuclei³⁰,³¹. For ¹H nuclei, this is then equilibrated throughout the system via efficient ¹H-¹H spin-diffusion, resulting in a large net polarization gain among the ¹H spins, as demonstrated in Fig. 2. Comparison with the

<figure><img src="image_3.png" /><figcaption>**Figure 1** | Representation of the unit cell of hexagonal HAp, Ca₁₀(PO₄)₆(OH)₂. Columns of OH⁻ (A-sites) are represented vertically (c axis) (Ca in green, O in red, H in white, PO₄ tetrahedra in brown). ICSD-26204: Hollysprings hydroxyapatite, showing the orientational disorder along the OH columns.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 2** | ¹H Hahn-echo MAS NMR spectra of C-HAp. All were recorded at ~100 K and 9.4 T. The pure powder (orange, MAS rate = 8.5 kHz) and the DNP-ready sample with (purple) and without (green) microwave irradiation suitable for DNP (both recorded using a MAS rate of ~14 kHz) are shown and scaled according to the factors given.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 3** | <sup>43</sup>Ca MAS NMR spectra of C-HAp. Both were recorded at ~100 K, 9.4 T and using a MAS rate of 8.5 kHz. Direct CT excitation of <sup>43</sup>Ca after DFS of the pure powder (orange) is compared with the DNP-enhanced {¹H-}<sup>43</sup>Ca CP experiment (purple), recorded using a CP contact time of 5 ms. Note that the time required to obtain much greater S/N with the DNP-enhanced experiment was <1/5 compared with the conventional method; the experimental time-saving with the DNP-enhanced experiment is a factor of 15.</figcaption></figure>

signal returned in the absence of µw irradiation and also with the signal for a full rotor of the pure powder only (that is, no solvent matrix or polarizing agent), highlights the significant gain in signal intensity through the use of MAS-DNP. The absolute sensitivity ratio (ASR)³² is the best method to compare gains with MAS-DNP to conventional MAS NMR and gives a factor of ~35 here for the OH⁻ groups, corresponding to an experimental time-saving of three orders of magnitude (35²). The OH⁻ groups are characterized by a sharp peak centred at δ<sub>iso</sub>(¹H) ~0 p.p.m., whereas a broad deshielded component assigned to hydrogen-bonded water/glycerol molecules is observed at δ<sub>iso</sub>(¹H) ~6.5 p.p.m. (independent experiments, not shown here, have demonstrated that HPO<sub>4</sub><sup>2-</sup> entities were almost entirely absent in the C-HAp sample). Note also the presence of a small, broad contribution at δ<sub>iso</sub>(¹H) ~6.8 p.p.m. for the pure powder, which arises from surface-absorbed water³³.

To better understand the process by which the sample becomes hyperpolarized from the polarizing agents that are contained in the solvent, ¹H NMR can be used here. The ¹H spectra of the DNP-ready sample in Fig. 2 were obtained at a MAS rate of ~14 kHz. The sharpening of the OH⁻ resonance compared with the spectrum recorded using a MAS rate of 8.5 kHz corroborates that ¹H-¹H spin-diffusion is still efficient under these moderate MAS rates as it implies the continued presence of residual ¹H-¹H dipolar couplings. This experimental fact is very important as it demonstrates that the ¹H nuclei are still communicating along the OH⁻ columns. The communication is established via the residual ¹H dipolar homonuclear interaction. The efficiency of this spin-diffusion process can be estimated more quantitatively: indeed, the ratio of signal intensities measured in the presence and absence of µw irradiation, ¹H<sub>εon/off</sub>, is 120 for water/glycerol and 110 for C-HAp OH⁻ groups. Notably, these discrete but similar values can be measured since the OH⁻ resonance remains clearly distinguished from the solvent contribution. Such a similarity tends to prove that almost all protons belonging to the C-HAp nanoparticles are involved in the DNP process. Interestingly though, after saturation, the monoexponential magnetization build-ups for OH⁻ and the solvent (in the presence of biradicals) have time constants, T<sub>B</sub>(¹H), that are quite different: 5 and 3 s, respectively. It follows that both ¹H spin baths can be considered as quite disconnected. However, connections are observed between the OH⁻ and solvent resonances in a DNP-enhanced ¹H-¹H spin-diffusion correlation experiment, albeit only at relatively long spin-diffusion mixing times

(Supplementary Fig. 5), indicating that the ¹H spin baths do share some small connection. Hence, the hyperpolarization of ¹H spins near to the polarizing agent can propagate from the solvent matrix into the nanoparticles via the OH⁻ channels (which are not large enough to accommodate the AMUPol molecules). In a second step, the ¹H hyperpolarization can then be transferred to other nuclei of interest (here <sup>43</sup>Ca (and <sup>13</sup>C)), for example through a CP³⁴ process.

**Studying Ca environments at natural abundance.** Many sensitivity enhancing experiments for solid-state NMR of quadrupolar nuclei have been devised³⁵ and one of the most efficient involving population transfer, besides DNP, is a so-called double-frequency sweep (DFS)³⁶ that can result in an increase of the population difference across the central transition (CT) by saturating the satellite transitions under MAS, giving NMR signal increases for the more readily observed CT of spin-7/2 nuclei of around a factor of 4 (for full satellite transition inversions a factor of 7 is attainable, at least theoretically)¹⁵. Figure 3 shows the <sup>43</sup>Ca MAS NMR CT spectrum of C-HAp, recorded on the pure powder, after DFS, at a sample temperature of 100 K. Already by only using the low temperatures, the NMR signal-to-noise ratio (S/N) will be ~3.6 times greater compared with conventional room temperature measurements³⁷ as a result of the increased thermal nuclear polarization and decreased Johnson-Nyquist noise in the radio-frequency circuit of the NMR probe. This, combined with the DFS, will result in experiments requiring only one two-hundredth of the conventional acquisition time (time-saving = (3.6 × 4)²). Nevertheless, even with the low temperatures and the DFS, the aforementioned difficulties of <sup>43</sup>Ca NMR cannot be surmounted and the S/N from the natural isotopic abundance sample is poor, even after 5.5 h of signal averaging. The long experimental times required to obtain sufficient S/N then render the acquisition of more complex experiments (such as 2D schemes) completely infeasible.

Figure 3 also shows the DNP-enhanced <sup>43</sup>Ca MAS NMR spectrum of C-HAp, recorded using a ¹H-<sup>43</sup>Ca CP step (that is, indirect DNP of <sup>43</sup>Ca). Although CP to half-integer spin quadrupolar nuclei is usually very inefficient owing to the
complex spin dynamics under MAS³⁸, particularly for the case of ⁴³Ca⁵, it is evident that the S/N is significantly better compared with the spectrum of the pure powder recorded with DFS, even with the much reduced experimental time (<1 h). It is worth highlighting that no data was recorded for the same experiment but in the absence of µw irradiation since this would have required months to obtain even a small signal. Therefore, the low efficiency of ¹H-⁴³Ca CP is completely negated by the large sensitivity enhancement from the MAS-DNP at low temperature, allowing fast acquisition of ⁴³Ca NMR spectra at natural isotopic abundance.

There are two non-equivalent calcium sites, Ca(I) and Ca(II) (hexagonal symmetry, Ca(II)/Ca(I) = 1.5, present in the structure of stoichiometric HAp (tentative fits of ⁴³Ca NMR spectra based on previously published data⁵,³⁹ are shown in Supplementary Fig. 6), with the Ca(II) sites being closer than Ca(I) sites to OH⁻ groups. The experimental spectra are compatible with the previously published data⁵,³⁹, though slightly broadened. Such broadening has already been described experimentally at moderate field (8.45 T) and is assigned to distributions of the quadrupolar parameters⁵,³⁹. At moderate fields, such distributions can play a significant role, especially in the case of nanoparticles for which surface and sub-surface defects are necessarily present in relatively high proportions. Furthermore, the substitution of some OH⁻/PO₄³⁻ groups by carbonates and some Ca²⁺ by Na⁺ (see the chemical analysis in the Supplementary Methods) can contribute as another source of local disorder and thus also result in distributions of quadrupolar parameters and associated spectral broadening. Indeed, such broadening due to extended substitutions has been observed in natural samples such as teeth and bones⁴⁰. It should also be noted that using CP can result in distorted quadrupolar lineshapes under MAS⁴¹.

**Interfaces in C-HAp nanoparticles.** As 1D DNP-enhanced {¹H – }⁴³Ca CP MAS NMR spectra could be recorded in a short experimental time for C-HAp, this experiment was easily extended to its 2D heteronuclear correlation (HETCOR) version⁴², which is used to probe proximities between the

heteronuclei. A ¹H-⁴³Ca HETCOR spectrum of C-HAp is presented in Fig. 4b. To the best of our knowledge, this is the first ever obtained without any ⁴³Ca labelling. Moreover, this was acquired in 15 h on a limited mass of sample (~30 mg). For comparison, a 2D ¹H-detected ¹H-⁴³Ca dipolar correlation spectrum of 30 mg 60% ⁴³Ca-enriched HAp required ~8 h of total experimental time at 14 T¹⁴, meaning that ~150 days would be needed for a natural abundance (0.14%) sample. Most interestingly, two correlation peaks are clearly distinguished using a short CP contact time of 3 ms (in contrast with experiments at longer contact times of 7 ms, in which only the one correlation peak at δᵢₐ(¹H) ~0 p.p.m. is present—see Supplementary Fig. 7). As expected, ⁴³Ca sites are correlated to the columnar OH⁻ groups at δᵢₐₕ(¹H) ~0 p.p.m. A strong correlation is also evidenced at δᵢₐ₅(¹H) ~3 p.p.m., although a ¹H resonance related to this chemical shift can only be observed as a shoulder in Fig. 2. The global linewidth of the ¹H projection centred at δᵢₐ₅(¹H) ~3 p.p.m. is estimated to ~2.7 p.p.m. (from 1.7 to 4.4 p.p.m.). Such a range of ¹H chemical shifts is compatible with those of hydrogen-bonded hydroxyls⁴³, indicating an assignment to surface OH⁻ groups that are hydrogen-bonded to surface-bound water or surface-adsorbed solvent molecules. Contributions from pure solvent protons can be ruled out because these are much more deshielded, at δᵢₐ₅(¹H) ~6.5 p.p.m. as shown in Fig. 2. This is corroborated by the ¹H-¹³C HETCOR spectrum of C-HAp, which shows correlations for the ¹³C resonances only with ¹H species with δᵢₐ₅(¹H) = 6.3 p.p.m. (FWHM ~3 p.p.m.) (Supplementary Fig. 7a). Since the DNP solvent mixture (glycerol-d₈/D₂O/H₂O (60/30/10; v/v/v)) was prepared at room temperature there will most likely be deuterium/proton equilibration of the exchangeable groups, resulting in 1.8 times more OH groups from water compared with glycerol, with both of these proton environments then contributing to this correlation at δᵢₐ₅(¹H) = 6.3 p.p.m. Furthermore, it should be noted that polyols such as glycerol (and carbohydrates as well) can chelate Ca²⁺ cations⁴⁴⁻⁴⁷. So it cannot be ruled out that there is a ¹H-⁴³Ca dipolar interaction between residual (non-hydrogen-bonded) glycerol molecules from the DNP solvent that are chelated to surface Ca²⁺ sites,

<figure><img src="image_6.png" /><figcaption>**Figure 4** | MAS-DNP HETCOR spectra of C-HAp. Both were recorded at ~100 K, 9.4 T, and using a MAS rate of 8.5 kHz. (a) ¹H-¹³C HETCOR spectrum with associated ¹³C skyline projection (top), recorded using a CP contact time of 7 ms. (b) ¹H-⁴³Ca HETCOR spectrum with associated ⁴³Ca (top) and ¹H (right) skyline projections, recorded using a CP contact time of 3 ms. Both spectra were recorded using FSLG homonuclear decoupling⁴⁸ during the indirect acquisition and SPINAL-64 heteronuclear decoupling⁴⁹ during the direct acquisition. Also shown (far right, purple) is the DNP-enhanced ¹H Hahn-echo MAS NMR spectrum, for comparison.</figcaption></figure>
leading to a correlation at $\delta_{iso}({}^1\text{H}) \sim 3$ p.p.m. Therefore, it follows that this correlation is then fully attributable to ${}^1\text{H}$ nuclei at the surface (that is, not in the columnar channels) of the C-HAp nanoparticles, either from $\text{OH}^-$ groups or chelating glycerol molecules, that are coupled to surface or near-surface $\text{Ca}^{2+}$ cations.

At longer ${}^1\text{H} - {}^{43}\text{Ca}$ CP contact times spin decoherence effects will become more important, reducing the efficiency of CP for local spin systems with larger ${}^1\text{H} - {}^1\text{H}$ dipolar couplings and those that are in closer proximity to the electron spins from the polarizing agent. Hence, CP for the isolated $\text{OH}^-$ groups of the columns (vide supra) should be less affected than for the surface $\text{OH}^-$ groups (or solvent molecules), which is indeed the case here (see the ${}^1\text{H} - {}^{43}\text{Ca}$ HETCOR spectrum of Supplementary Fig. 7b where a CP contact time of 7 ms was used, showing only a dominant correlation peak at $\delta_{iso}({}^1\text{H}) \sim 0$ p.p.m.). The correlation centred at $\delta_{iso}({}^1\text{H}) \sim 0$ p.p.m. in Fig. 4b is assigned to ${}^{43}\text{Ca}$ nuclei located near the hydroxyl columns, that is, in the core of the nanoparticles. A comparison of the ${}^{43}\text{Ca}$ cross-sections taken through the surface and core correlations shows that the surface related ${}^{43}\text{Ca}$ spectrum is slightly broader (Supplementary Fig. 8). This could be expected because of structural disorder at the surface, leading to additional distributions of chemical shift and/or quadrupolar parameters.

The ${}^1\text{H}-{}^{13}\text{C}$ HETCOR spectrum of C-HAp is presented in Fig. 4a. It is clearly demonstrated that even in the case of ${}^{13}\text{C}$ at natural abundance (1.1%) and with only $\sim 1$ wt% in carbon, 2D DNP-enhanced ${}^1\text{H}-{}^{13}\text{C}$ HETCOR spectra can be obtained in short experimental time (here 5.5 h). One notes also that ${}^{13}\text{C}$ resolution remains almost identical to that observed at room temperature. Following Beshah et al.$^{48}$, the most shielded ${}^{13}\text{C}$ resonance ($\delta_{iso}({}^{13}\text{C}) = 166$ p.p.m.) can be assigned to A-site carbonates whereas all deshielded resonances correspond to B-site carbonates. Looking carefully at the correlations in the 2D ${}^1\text{H}-{}^{13}\text{C}$ HETCOR experiment provides more detailed information. (i) As expected, ${}^{13}\text{C}$ resonances are coupled to channel $\text{OH}^-$ protons, corresponding then to core carbonates. (ii) Some B-type resonances are correlated to the most deshielded part of the ${}^1\text{H}$ spectrum, corresponding then to surface carbonates. There seems to be two discrete ${}^1\text{H}$ environments associated to these surface carbonates: at $\delta_{iso}({}^1\text{H}) \sim 7$ p.p.m. and $\delta_{iso}({}^1\text{H}) \sim 5$ p.p.m., most likely from distinct surface OH groups. Notably, the proton correlations related to these surface carbonates are clearly distinguished from the one related to the surface $\text{Ca}^{2+}$ cations. In other words, the multinuclear MAS-DNP experiments are selective towards the different protonated species present at the surface of the C-HAp nanoparticles. (iii) The A-site resonance ($\delta_{iso}({}^{13}\text{C}) = 166$ p.p.m.) is mainly correlated to the $\text{OH}^-$ groups in the c axis columns. (iv) The lack of any ${}^{13}\text{C}$ correlations with a peak at $\delta_{iso}({}^1\text{H}) \sim 3$ p.p.m. (using CP contact times of 0.25, 1 and 7 ms (only the latter shown here)) allows the refinement of the previous assignment (vide supra). This peak must then be specifically from (hydrogen-bonded) calcium-bound hydroxyls on the surface of the nanoparticles, or from calcium-chelated glycerol molecules (which have a lesser affinity to carbonates).

The ability to selectively probe either the surface or the core of materials based on HAp at natural isotopic abundance is important, particularly so for the former. Indeed, the structure of tissues like bone and teeth is still the matter of much debate$^{49-52}$. Many recent studies show that there is a real need to develop new ways of analysing these materials at different length-scales, and notably at the molecular level$^{53}$, to determine how the different individual components (proteins, mineral, water molecules and so on) are associated together, and how this can have an impact on the functional properties of the tissue. In this context, the

organic-mineral interface has drawn much attention, revealing the importance of water molecules, small molecular anions like citrate, and non-collagenous proteins like osteopontin and osteocalcin, though detailed insight on their mode of binding to bone mineral is still missing. Given that these species are expected to coordinate to the surface $\text{Ca}^{2+}$ sites (through ionic interactions and/or coordination bonds), directly analysing the surface $\text{Ca}^{2+}$ sites of HAp in native tissue would be a way to directly investigate this mode of binding. The natural isotopic abundance MAS-DNP experiments proposed here illustrate that such experiments are now feasible, which clearly opens avenues for understanding the structure of the mineral interface in these important materials. A preliminary attempt was made to record a DNP-enhanced $\{ {}^1\text{H} - {}^{43}\text{Ca} \} \text{CP MAS NMR spectrum of mice teeth}$ (Supplementary Fig. 9). The observation of a ${}^{43}\text{Ca}$ NMR signal shows the applicability of the experiment to natural samples, although further investigation into the sample preparation and acquisition conditions is required to optimize the MAS-DNP experiments.

It has been demonstrated that C-HAp nanoparticles can be studied efficiently and effectively using MAS-DNP. Although ${}^{43}\text{Ca}$ should be a fundamental NMR target in the field of biomaterials and biominerals, it suffers from severe drawbacks in terms of intrinsic sensitivity (low $\gamma$, high spin number and very low natural abundance). Nevertheless, ${}^1\text{H}-{}^{43}\text{Ca}$ 2D experiments were recorded at natural isotopic abundance in short experimental time, notably for a restricted mass of sample ($\sim 30$ mg). To date, such experiments would have been impossible to implement using conventional NMR. An analysis of the ${}^1\text{H}-{}^{43}\text{Ca}$ and ${}^1\text{H}-{}^{13}\text{C}$ correlation data showed that surface species ($\text{Ca}^{2+}$ and carbonates) could be distinguished from the core entities. Such results open avenues for the detailed atomic-scale description of interfaces in bones and teeth, which is of primary importance to understand their structure and properties$^{53}$, as well as for calcium-derived nanomaterials in general.

## Methods

**Carbonate substituted hydroxyapatite (C-HAp).** The protocol used corresponded to the standard precipitation of HAp at room temperature starting from aqueous solutions (initially decarbonated by boiling the solutions and using an argon flux during the syntheses). A solution of ammonium hydrogenphosphate (NH$_4$)$_2$HPO$_4$, 0.30 $\text{mol l}^{-1}$, pH $\sim 10$) was added to a solution of calcium nitrate tetrahydrate (Ca(NO$_3$)$_2$·4H$_2$O, 0.50 $\text{mol l}^{-1}$, pH $\sim 5$) using a titration apparatus (808 Titrando, Methrom, 3 ml min $^{-1}$). Carbonates were introduced by adding sodium hydrogencarbonate (NaHCO$_3$) to the phosphate solution. All experiments were performed under argon flux and magnetic stirring. The initial molar ratios were: Ca:10, P:6, $CO_3^{2-}$:3. After the addition of the phosphate solution, the obtained precipitates were further stirred for 24 h. After centrifugation, the precipitates were rinsed four times with distilled water (final pH of water $\sim 7$). All samples were then heat treated (400 °C) for 48 h to eliminate water and ammonia molecules. Routine characterization (TEM, XRD, FTIR, TGA) of the sample is presented in Supplementary Figs 1-4.

**Solid-state NMR experiments.** The glycerol-d$_8$ (98% in D atoms) and D$_2$O (> 99% in D atoms) used for the glass-forming solvent matrix were standard commercial reagents and the AMUPol$^{28}$ biradical polarizing agent was purchased through SATT Sud Est, Marseille, France. Solid-state NMR experiments were recorded using a Bruker DNP-NMR AVANCE III 400 MHz spectrometer equipped with a gyrotron and associated transmission line capable of delivering $\sim 5$ W of $\sim 264$ GHz µw irradiation at the sample. All experiments were recorded with a 3.2 mm HXY triple-resonance MAS probe (in $\lambda/4$ ${}^1\text{H}$ transmission mode) at $\nu_0({}^1\text{H}) = 400.33$ MHz corresponding to the maximum ${}^1\text{H}$ enhancement field position for AMUPol at $\nu_0(\text{e}^-) = 263.7$ GHz, with the X-channel tuned to $\nu_0({}^{13}\text{C}) = 100.66$ MHz and the Y-channel tuned to $\nu_0({}^{43}\text{Ca}) = 26.94$ MHz. A home-adapted triple-resonance ${}^1\text{H}/{}^{13}\text{C}/{}^{43}\text{Ca}$ circuit was constructed using a 15 pF capacitor in parallel to an insert coil of seven turns, used to bridge the X and Y channels, and also a 64 pF shunt capacitance in parallel to the variable Y capacitor. The Y-channel tuning range at low temperature ($\sim 100$ K) using this setup is $26.9 \pm 0.8$ MHz.

${}^1\text{H}$ ($\pi/2$-) and $\pi$-pulses, as well as the (SPINAL-64 (ref. 54)) heteronuclear decoupling, were applied to induce a nutation frequency of 100 kHz. ${}^{13}\text{C}$ CP spin-locking was at $\sim 50$ kHz with corresponding ramped (50-100%) ${}^1\text{H}$ CP
spin-locking at ~76 kHz (100%) to meet a CP-matching condition and facilitate CP transfer. The ¹H spin-lock during ¹H-⁴³Ca CP transfer was applied in a low nutation frequency regime (~12 kHz) to permit relatively efficient CT-selective CP of ⁴³Ca nuclei, for which a ($n=1$, zero quantum) CP-matching condition was met using ⁴³Ca radio-frequency irradiation at ~0.9 kHz. A long dead time of 120 µs was used to avoid the acquisition of ⁴³Ca signals resulting from probe ringing. For ¹³C acquisition, this dead time was only 7 µs. The {¹H-}⁴³Ca CP and ⁴³Ca direct excitation after DFS preparation were calibrated using a (60%) ⁴³Ca-labelled sample of Mg-substituted HAP, which was synthesized using a similar protocol to that described previously¹¹. For the ⁴³Ca NMR spectra shown in Fig. 3, the recycle delays were 6.5 s for the MAS-DNP {¹H-}⁴³Ca CP NMR experiment (value directly optimized on the sample) and 1 s for the DFS-⁴³Ca NMR experiment (the ⁴³Ca relaxation of precipitated hydroxyapatite is known to be fast compared with other Ca-containing materials⁵⁵).

**Data availability.** The data that support the findings of this study are available from the corresponding author upon reasonable request.

## References

1. MacKenzie, K. J. D. & Smith, M. E. *Multinuclear Solid-State Nuclear Magnetic Resonance of Inorganic Materials* (Pergamon Press, 2002).

2. Bryce, D. L. Calcium binding environments probed by ⁴³Ca NMR spectroscopy. *Dalton Trans.* **39**, 8593–8602 (2010).

3. Moudrakovski, I. L. in *Annual Reports on NMR Spectroscopy* (ed. Webb, G. A.) 129–240 (Elsevier Academic Press Inc, 2013).

4. Schurko, R. W. & Jaroszewicz, M. J. in *The Lightest Metals: Science and Technology from Lithium to Calcium* (ed. Hanusa, T. P.) 117–172 (John Wiley & Sons, Ltd, 2015).

5. Laurencin, D. & Smith, M. E. Development of ⁴³Ca solid state NMR spectroscopy as a probe of local structure in inorganic and molecular materials. *Prog. Nucl. Magn. Reson. Spectrosc.* **68**, 1–40 (2013).

6. Dupree, R., Howes, A. P. & Kohn, S. C. Natural abundance solid state ⁴³Ca NMR. *Chem. Phys. Lett.* **276**, 399–404 (1997).

7. Bryce, D. L., Bultz, E. B. & Aebi, D. Calcium-43 chemical shift tensors as probes of calcium binding environments. Insight into the structure of the vaterite CaCO₃ polymorph by ⁴³Ca solid-state NMR spectroscopy. *J. Am. Chem. Soc.* **130**, 9282–9292 (2008).

8. Xu, J. *et al.* Natural-abundance ⁴³Ca solid-state NMR spectroscopy of bone. *J. Am. Chem. Soc.* **132**, 11504–11509 (2010).

9. Colas, H. *et al.* Whewellite, CaC₂O₄·H₂O: structural study by a combined NMR, crystallography and modelling approach. *CrystEngComm.* **15**, 8840–8847 (2013).

10. Widdifield, C. M., Moudrakovski, I. & Bryce, D. L. Calcium-43 chemical shift and electric field gradient tensor interplay: a sensitive probe of structure, polymorphism, and hydration. *Phys. Chem. Chem. Phys.* **16**, 13340–13359 (2014).

11. Laurencin, D., Wong, A., Hanna, J. V., Dupree, R. & Smith, M. E. A high-resolution ⁴³Ca solid-state NMR study of the calcium sites of hydroxyapatite. *J. Am. Chem. Soc.* **130**, 2412–2413 (2008).

12. Laurencin, D. *et al.* Implementation of high resolution ⁴³Ca solid state NMR spectroscopy: toward the elucidation of calcium sites in biological materials. *J. Am. Chem. Soc.* **131**, 13430–13440 (2009).

13. Wong, A., Aguiar, P. M., Charpentier, T. & Sakellariou, D. A low-cost strategy for ⁴³Ca solid-state NMR spectroscopy. *Chem. Sci.* **2**, 815–818 (2011).

14. Wong, A., Laurencin, D., Dupree, R. & Smith, M. E. Two-dimensional ⁴³Ca-¹H correlation solid-state NMR spectroscopy. *Solid State Nucl. Magn. Reson.* **35**, 32–36 (2009).

15. Burgess, K. M. N., Perras, F. A., Moudrakovski, I. L., Xu, Y. & Bryce, D. L. High sensitivity and resolution in ⁴³Ca solid-state NMR experiments. *Can. J. Chem.* **93**, 799–807 (2015).

16. Hall, D. A. *et al.* Polarization-enhanced NMR spectroscopy of biomolecules in frozen solution. *Science* **276**, 930–932 (1997).

17. Ni, Q. Z. *et al.* High frequency dynamic nuclear polarization. *Acc. Chem. Res.* **46**, 1933–1941 (2013).

18. Rossini, A. J. *et al.* Dynamic nuclear polarization surface enhanced NMR spectroscopy. *Acc. Chem. Res.* **46**, 1942–1951 (2013).

19. Zagdoun, A. *et al.* A slowly relaxing rigid biradical for efficient dynamic nuclear polarization surface-enhanced NMR spectroscopy: expeditious characterization of functional group manipulation in hybrid materials. *J. Am. Chem. Soc.* **134**, 2284–2291 (2012).

20. Blanc, F. *et al.* Dynamic nuclear polarization NMR of low-γ nuclei: structural insights into hydrated yttrium-doped BaZrO₃. *J. Phys. Chem. Lett.* **5**, 2431–2436 (2014).

21. Guo, Z. *et al.* Selective host-guest interaction between metal ions and metal-organic frameworks using dynamic nuclear polarization enhanced solid-state NMR spectroscopy. *Chem. Eur. J.* **20**, 16308–16313 (2014).

22. Gutmann, T. *et al.* Natural abundance ¹⁵N NMR by dynamic nuclear polarization: fast analysis of binding sites of a novel amine-carboxyl-linked immobilized dirhodium catalyst. *Chem. Eur. J.* **21**, 3798–3805 (2015).

23. Lee, D., Duong, N. T., Lafon, O. & De Paëpe, G. Primostrato solid-state NMR enhanced by dynamic nuclear polarization: penta-coordinated Al are only located at the surface of hydrated γ-alumina. *J. Phys. Chem. C* **118**, 25065–25076 (2014).

24. Blanc, F. *et al.* Dynamic nuclear polarization enhanced natural abundance ¹⁷O spectroscopy. *J. Am. Chem. Soc.* **135**, 2975–2978 (2013).

25. Perras, F. A., Kobayashi, T. & Pruski, M. Natural abundance ¹⁷O DNP two-dimensional and surface-enhanced NMR spectroscopy. *J. Am. Chem. Soc.* **137**, 8336–8339 (2015).

26. Fleet, M. E. *Carbonated Hydroxyapatite: Materials, Synthesis, and Applications* (Pan Stanford Publishing Pte Ltd, 2015).

27. Klimavicius, V., Dagys, L. & Balevicius, V. Subnanoscale order and spin diffusion in complex solids through the processing of cross-polarization kinetics. *J. Phys. Chem. C* **120**, 3542–3549 (2016).

28. Sauvée, C. *et al.* Highly efficient, water-soluble polarizing agents for dynamic nuclear polarization at high frequency. *Angew. Chem. Int. Ed.* **52**, 10858–10861 (2013).

29. Lesage, A. *et al.* Surface enhanced NMR spectroscopy by dynamic nuclear polarization. *J. Am. Chem. Soc.* **132**, 15459–15461 (2010).

30. Mentink-Vigier, F. *et al.* Fast passage dynamic nuclear polarization on rotating solids. *J. Magn. Reson.* **224**, 13–21 (2012).

31. Thurber, K. R. & Tycko, R. Theory for cross effect dynamic nuclear polarization under magic-angle spinning in solid state nuclear magnetic resonance: the importance of level crossings. *J. Chem. Phys.* **137**, 084508 (2012).

32. Takahashi, H. *et al.* Rapid natural-abundance 2D ¹³C-¹³C correlation spectroscopy using dynamic nuclear polarization enhanced solid-state NMR and matrix-free sample preparation. *Angew. Chem. Int. Ed.* **51**, 11766–11769 (2012).

33. Yesinowski, J. P. & Eckert, H. Hydrogen environments in calcium phosphates: proton MAS NMR at high spinning speeds. *J. Am. Chem. Soc.* **109**, 6274–6282 (1987).

34. Pines, A., Gibby, M. G. & Waugh, J. S. Proton-enhanced nuclear induction spectroscopy. A method for high resolution NMR of dilute spins in solids. *J. Chem. Phys.* **56**, 1776–1777 (1972).

35. Perras, F. A., Viger-Gravel, J., Burgess, K. M. N. & Bryce, D. L. Signal enhancement in solid-state NMR of quadrupolar nuclei. *Solid State Nucl. Magn. Reson.* **51-52**, 1–15 (2013).

36. Iuga, D., Schäfer, H., Verhagen, R. & Kentgens, A. P. Population and coherence transfer induced by double frequency sweeps in half-integer quadrupolar spin systems. *J. Magn. Reson.* **147**, 192–209 (2000).

37. Takahashi, H. *et al.* Optimization of an absolute sensitivity in a glassy matrix during DNP-enhanced multidimensional solid-state NMR experiments. *J. Magn. Reson.* **239**, 91–99 (2014).

38. Ashbrook, S. E. & Wimperis, S. Spin-locking of half-integer quadrupolar nuclei in nuclear magnetic resonance of solids: second-order quadrupolar and resonance offset effects. *J. Chem. Phys.* **131**, 194509 (2009).

39. Laurencin, D., Wong, A., Dupree, R. & Smith, M. E. Natural abundance ⁴³Ca solid-state NMR characterisation of hydroxyapatite: identification of the two calcium sites. *Magn. Reson. Chem.* **46**, 347–350 (2008).

40. Laurencin, D. *et al.* Probing the calcium and sodium local environment in bones and teeth using multinuclear solid state NMR and X-ray absorption spectroscopy. *Phys. Chem. Chem. Phys.* **12**, 1081–1091 (2010).

41. Hayashi, S. & Hayamizu, K. Line shapes in CP/MAS NMR spectra of half-integer quadrupolar nuclei. *Chem. Phys. Lett.* **203**, 319–324 (1993).

42. Roberts, J. E., Vega, S. & Griffin, R. G. Two-dimensional heteronuclear chemical shift correlation spectroscopy in rotating solids. *J. Am. Chem. Soc.* **106**, 2506–2512 (1984).

43. Grünberg, B. *et al.* Hydrogen bonding of water confined in mesoporous silica MCM-41 and SBA-15 studied by ¹H solid-state NMR. *Chem. Eur. J.* **10**, 5689–5696 (2004).

44. Charley, P. & Saltman, P. Chelation of calcium by lactose: its role in transport mechanisms. *Science* **139**, 1205–1206 (1963).

45. Fujii, K. & Kondo, W. Calcium glycerglycerides formed in the system of calcium oxide-glycerol. With 4 figures. *Z. Anorg. Allg. Chem.* **359**, 296–304 (1968).

46. Bugg, C. E. Calcium binding to carbohydrates: crystal structure of a hydrated calcium bromide complex of lactose. *J. Am. Chem. Soc.* **95**, 908–913 (1973).

47. Dheu-Andries, M. L. & Pérez, S. Geometrical features of calcium—carbohydrate interactions. *Carbohydr. Res.* **124**, 324–332 (1983).

48. Beshah, K., Rey, C., Glimcher, M. J., Schimizu, M. & Griffin, R. G. Solid state carbon-13 and proton NMR studies of carbonate-containing calcium phosphates and enamel. *J. Solid State Chem.* **84**, 71–81 (1990).

49. Hu, Y.-Y., Rawal, A. & Schmidt-Rohr, K. Strongly bound citrate stabilizes the apatite nanocrystals in bone. *Proc. Natl Acad. Sci. USA* **107**, 22425–22429 (2010).

50. Wang, Y. *et al.* Water-mediated structuring of bone apatite. *Nat. Mater.* **12**, 1144–1153 (2013).
51. Davies, E. et al. Citrate bridges between mineral platelets in bone. *Proc. Natl Acad. Sci. USA* **111**, E1354–E1363 (2014).

52. Stock, S. R. The mineral–collagen interface in bone. *Calcif. Tissue Int.* **97**, 262–280 (2015).

53. Duer, M. J. The contribution of solid-state NMR spectroscopy to understanding biomineralization: atomic and molecular structure of bone. *J. Magn. Reson.* **253**, 98–110 (2015).

54. Fung, B. M., Khitrin, A. K. & Ermolaev, K. An improved broadband decoupling sequence for liquid crystals and solids. *J. Magn. Reson.* **142**, 97–101 (2000).

55. Gervais, C. et al. New perspectives on calcium environments in inorganic materials containing calcium-oxygen bonds: a combined computational-experimental $^{43}$Ca NMR approach. *Chem. Phys. Lett.* **464**, 42–48 (2008).

56. Bielecki, A., Kolbert, A. C. & Levitt, M. H. Frequency-switched pulse sequences: homonuclear decoupling and dilute spin NMR in solids. *Chem. Phys. Lett.* **155**, 341–346 (1989).

## Acknowledgements

This work was supported by the French National Research Agency through the Labex ARCANE (ANR-11-LABX-0003-01), the 'programme blanc' (ANR-12-BS08-0016-01) and the European Research Council (ERC-CoG-2015, No. 682895). The RTB is acknowledged for financial support. Katharina Märker is thanked for her help in setting up the experiments and contributions from Sabine Hediger, Frédéric Mentink-Vigier, Dimitris Sakellariou and Jochen Lehmann-Horn are greatly appreciated. Dr Sara Laurencin-Dalicieux (Centre de Physiopathologie Toulouse-Purpan, INSERM UMR 1043, France) is gratefully acknowledged for providing the mice teeth samples.

## Author contributions

D.Lee designed and conducted the solid-state NMR characterization, and drafted the manuscript with C.B. and G.D.P. C.L. and C.C. prepared the C–HAp sample. L.B.-C.

and F.B. contributed to the analysis of the results. D.Laurencin contributed to the $^{43}$Ca NMR experiments, and also prepared the $^{43}$Ca-labelled sample of Mg-substituted HAp used for calibration. C.B. and G.D.P. designed the project. All authors contributed to the writing of the manuscript.

## Additional information

**Supplementary Information** accompanies this paper at http://www.nature.com/naturecommunications

**Competing financial interests:** The authors declare no competing financial interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**How to cite this article:** Lee, D. et al. Interfacial $Ca^{2+}$ environments in nanocrystalline apatites revealed by dynamic nuclear polarization enhanced $^{43}$Ca NMR spectroscopy. *Nat. Commun.* **8**, 14104 doi: 10.1038/ncomms14104 (2017).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2017