The world's population is estimated to rise to 11 billion by 2100,\[[@b1]\] putting unparalleled pressure on agricultural food production.\[[@b2]\] New sustainable means of increasing land fertility whilst avoiding the current overuse of NPK fertilizers are needed.\[[@b2c], [@b3]\] Humic substances, which exist in soil as complex mixtures of thousands of organic compounds derived from animal and plant remains, could form an integral part of green and sustainable solutions to the soil fertility problem.  Humic substances play a vital role in physical and biochemical soil functions,\[[@b4]\] stimulate plant growth,\[[@b5]\] and ameliorate the effect of contaminant metals and chemical residues, which may build up over time and hinder plant growth.\[[@b6]\] It is unclear, however, if all fractions of humic substances are beneficial to plant growth,\[[@b5c]\] and if certain humic molecules may indeed be harmful to humans.\[[@b7]\] Thus, to address the potential role of humic substances in restoring the fertility of intensively farmed agricultural soils, it is essential to determine their molecular composition.  Such an achievement will allow development of a molecular rather than a phenomenological description of soil---a fundamental step toward defining the structure--function relationships of humic substances.  Nevertheless, the two high-resolution analytical techniques, Fourier transform ion cyclotron resonance mass spectrometry (FT ICR MS)\[[@b8]\] and NMR spectroscopy,\[[@b6b], [@b9]\] which are essential in this endeavor, have serious shortcomings: MS only provides molecular formulae, and standard 2D or 3D NMR experiments cannot deliver the structures of individual compounds contained in complex mixtures. 

The development of more powerful multidimensional (*n*D) NMR experiments is essential in overcoming these limitations.\[[@b10]\] However, *n*D NMR alone cannot address the inherent complexity of chromatographically inseparable mixtures,\[[@b10]\] and some form of "spectroscopic separation" is still necessary.  To achieve this, we have recently developed an isotope-filtered *n*D NMR methodology\[[@b11]\]---a combination of isotopic tagging and *n*D NMR.  Unlike previous tagging approaches that focused on the signals of the tags,\[[@b6b], [@b12]\] we utilize them to probe directly the identity of the tagged molecules.  The power of this approach is in the ability to provide multiple correlated chemical shifts of individual molecules. 

Here we report how isotope-filtered *n*D NMR enables characterization of phenolic moieties of humic molecules and illustrate its principles using a 4D ^13^CH~3~O-filtered NMR experiment, 4D HCCH~3~, which correlates chemical shifts of four nuclei---the aromatic CH atoms *ortho* to methoxy groups and those of ^13^CH~3~O atoms.  The polarization transfer pathway of this experiment starts on aromatic protons and finishes on methyl protons utilizing ^1^*J*, ^3^*J*~CC~, and ^1^*J* couplings (Figure [1 a](#fig01){ref-type="fig"}).  The pulse sequence of 4D HCCH~3~ is shown in [Figure S1](#SD1){ref-type="supplementary-material"} in the Supporting Information. 

![ a) An exemplar aromatic compound highlighting nuclei accessible by ^13^CH~3~O-filtered *n*D NMR experiments.  The interactions mediating the polarization transfers in the 4D HCCH~3~ are shown as dashed double-headed arrows.  The color coding of nuclei is the same as used for the chemical shift axes of *n*D NMR spectra; b) the methoxy region of the 800 MHz 2D ^1^H,^13^C HSQC spectrum of ^13^CH~3~O-methylated FA.  Circled areas identify subregions as ester (orange), phenolic (cyan), carbohydrate (green), and aromatic/aliphatic sterically hindered methoxy groups (magenta); c) A stack plot of the phenolic methoxy resonances corresponding to the cyan-circled cross peaks in (b).  The region of the spectrum used to illustrate the 4D HCCH~3~ experiment below is labeled with an asterisk.  d) Exemplar 3D cuboids extracted from a 800 MHz 4D HCCH~3~ spectrum of ^13^C-methylated FA at ^13^C(H~3~) chemical shifts indicated by the dashed lines.  e) An F~1~F~4~ (H~ar~(C)H~3~) projection of a cuboid extracted at 55.36 ppm; f) An overlay of three 2D F~1~F~2~ (H~ar~C~ar~) planes extracted from this cuboid at proton methoxy chemical shifts indicated by the colored dashed lines in (e).  Insets show the identified structural fragments, which belong to compounds 23, 32, and 33 ([Figure S2](#SD1){ref-type="supplementary-material"}).](anie0054-8382-f1){#fig01}

The workings of this experiment are illustrated using a peat soil fulvic acid (FA) sample that was methylated using ^13^CH~3~I.\[[@b11]\] The sample was first characterized by acquiring a 2D ^1^H,^13^C HSQC spectrum (Figure [1 b](#fig01){ref-type="fig"}).  Its methoxy region shows a spread of ^13^CH~3~ cross peaks reflecting different chemical environments of the methoxy groups.  As the 4D HCCH~3~ experiment was designed to investigate phenolic compounds, an expansion of the phenolic methoxy region of the 2D ^1^H,^13^C HSQC spectrum is shown as a stacked plot in Figure [1 c](#fig01){ref-type="fig"}.  This presentation reveals the presence of several major and numerous minor phenolic compounds in this FA sample. 

A partial 4D HCCH~3~ spectrum of ^13^C-methylated FA shown in Figure [1 d](#fig01){ref-type="fig"} illustrates how isotope filtration combined with the dispersion of signals in a 4D space leads to comprehensive structural information.  A cuboid extracted at 55.36 ppm is inspected.  This region of the spectrum contains severely overlapping medium intensity signals labeled with an asterisk on Figure [1 c](#fig01){ref-type="fig"}.  An F~1~F~4~ (or H~ar~(C~ar~COC)H~3~) projection of the 55.36 ppm cuboid shows an adequate spread of signals in the directly detected F~4~ dimension (Figure [1 e](#fig01){ref-type="fig"}).  An overlay of three 2D F~1~F~2~ (or H~ar~C~ar~) planes taken from this cuboid at methyl proton chemical shifts indicated by the dashed lines in Figure [1 e](#fig01){ref-type="fig"} reveals similar, but non-identical, H~ar~/C~ar~ chemical shifts (Figure [1 f](#fig01){ref-type="fig"}).  The inspected region of the 4D HCCH~3~ spectrum thus provides unprecedented chemical shift separation of aromatic compounds which have identical aromatic ring substitution patterns, whereas their substituents show structural variations further away from the aromatic ring. 

Additional 3D ^13^CH~3~O-filtered experiments were designed to complement the 4D experiment.  3D variants of the 4D HCCH~3~, 3D H(C)CH~3~, and 3D (H)CCH~3~, do not label aromatic carbon or proton resonances, respectively, and hence provide similar, but more ambiguous information than the 4D spectra albeit with higher digital resolution.  Chemical shifts of nuclei other than those *ortho* to the ^13^CH~3~O groups ([Table S1](#SD1){ref-type="supplementary-material"}), were obtained by 3D INEPT-INADEQUATE-HSQC\[[@b11]\]/3D C~q~CH~3~ (quaternary aromatic carbons) and 3D CH~3~-NOESY/3D CH~3~-NOESY-TOCSY (*meta* and *para* protons relative to ^13^CH~3~O groups).  These NMR experiments utilize polarization transfer through various couplings shown in Figure [1 a](#fig01){ref-type="fig"}.  The NOESY- and TOCSY- based experiments also show the splitting of aromatic proton resonances due to ^3^*J*~HH~ or ^4^*J*~HH~ couplings, hence providing additional valuable structural information. 

In summary, multiple chemical shifts and coupling constants were obtained for individual molecules by analyzing the 3D/4D ^13^CH~3~O-filtered spectra.  When cross-referenced against database information,\[[@b13]\] this lead to the identification of nine major phenolic compounds of this FA sample (Figure [2](#fig02){ref-type="fig"}).  These comprise 1,3,4-, 1,3,5-trisubstituted as well as 1,4-disubstituted hydroxybenzenes that differ by the nature of the *para* substituent relative to the ^13^CH~3~O group. 

![ Major phenolic compounds identified in the methylated FA sample; their methoxy cross peaks are labeled on the 2D ^1^H,^13^C HSQC spectrum.](anie0054-8382-f2){#fig02}

Twenty-one additional structures/structural motifs were also identified ([Figure S2](#SD1){ref-type="supplementary-material"}).  Some of these contained highly substituted aromatic rings with four or five substituents.  These molecules mainly relate to plant-derived lignin precursors such as vanillin, coumaryl, syringyl, coniferyl monolignols, but also flavonoids.  The presence of several of these molecules in humic substances has been postulated based on the results of tetramethyl ammonium hydroxide (TMAH)-assisted pyrolysis GC-MS analysis.\[[@b14]\] Aromatic ring substitution patterns identified in our work were previously suggested based on the bulk matching of experimental and theoretical 2D ^1^H,^13^C HSQC spectra; the latter were calculated by considering the effects of -OR and -COOR groups on the ^1^H and ^13^C chemical shift of CH atoms.\[[@b15]\] These approaches relied on matching a limited number of descriptors to numerous possibilities.  In contrast, our methodology provides a multitude of correlated chemical shifts, which allows convergence to a single structure/structural motif, thus yielding unprecedented structural details for phenolic compounds in samples of humic substances. 

Further advances of this methodology are not limited to methylation.  In addition, other tags containing NMR-active nuclei such as ^15^N, ^19^F, and ^31^P, represent promising candidates for tagging various functional groups.  Once fully developed, this methodology will lay ground for the structure--function investigations of humic compounds thus enabling exploration of their roles in improving the soil fertility and sustainable food production.  Isotope-filtered *n*D NMR spectroscopy is also applicable to the investigation of natural organic matter in general, but also complex mixtures other than those found in the environment, for example, small-molecular metabolites, plant extracts, oil, food, and beverages. 

Experimental Section
====================

Peaty soil was collected from the Red Moss, an ombrotrophic peat bog near Balerno, central Scotland.  The FA sample (60 mg) was extracted using the International Humic Substance Society (IHSS) protocols (<http://www.humicsubstances.org/isolation.html>) and methylated as described previously\[[@b11]\] yielding ^13^C-methylated FA (10 mg).  The sample was dissolved in 550 μL of CDCl~3~ for NMR analysis at 15 °C.  The pulse sequence of the 4D HCCH~3~ is shown in [Figure S1](#SD1){ref-type="supplementary-material"}.  The parameters of NMR experiments are detailed in the Supporting Information. 

Supporting Information
======================

As a service to our authors and readers, this journal provides supporting information supplied by the authors.  Such materials are peer reviewed and may be re-organized for online delivery, but are not copy-edited or typeset.  Technical support issues arising from supporting information (other than missing files) should be addressed to the authors. 

###### 

miscellaneous_information

[^1]: This work was supported by the NERC grant NE/L00044X/1 to M.C.G. and D.U. N.G.A.B. would like to acknowledge the support of the University of Edinburgh Principal\'s Career Development Scholarship and NERC.  We thank J. Bella and Dr. L. Murray for maintaining NMR spectrometers.  The peat sample was collected with the permission of the Scottish Wildlife Trust. 

[^2]: Supporting information for this article is available on the WWW under <http://dx.doi.org/10.1002/anie.201503321>.