Isotope-filtered nD NMR spectroscopy of complex mixtures to unravel the molecular structures of phenolic compounds in tagged soil organic matter

Unravelling structures of molecules contained in complex, chromatographically inseparable mixtures is a challenging task. Due to the number of overlapping resonances in NMR spectra of these mixtures, unambiguous chemical shift correlations attributable to individual molecules cannot be achieved and thus their structure determination is elusive by this technique. Placing a tag carrying an NMR active nucleus onto a subset of molecules enables (i) to eliminate signals from the non-tagged molecules, and (ii) to obtain a set of correlated chemical shifts and coupling constants belonging to a single molecular type. This approach provides an opportunity for structure determination without the need for compound separation. Focusing on the most abundant functional groups of natural organic matter molecules, the carboxyl and hydroxyl groups were converted into esters and ethers, respectively by introducing (13)CH3O groups. A set of (13)C-filtered nD NMR experiments was designed yielding structures/structural motives of tagged molecules. The relative sensitivity of these experiments was compared and a step-by-step guide how to use these experiments to analyse the structures of methylated phenolics is provided. The methods are illustrated using an operational fraction of soil organic matter, fulvic acid isolated from a Scottish peat bog. Analysis of 33 structures identified in this sample revealed a correlation between the position of the methoxy cross-peaks in the (1)H, (13)C HSQC spectra and the compound type. This information enables profiling of phenolic compounds in natural organic matter without the need to acquire a full set of experiments described here or access to high field cryoprobe NMR spectrometers.


Introduction
Natural organic matter (NOM), an assembly of thousands of organic molecules, is regarded as the most complex mixture on Earth. 1 Formed by biotic and abiotic degradation of plant and animal residues, 2 NOM represents an immense continuum of organic structures found in soils, natural waters and air.4][5] Despite its significance, the composition and the transformations of NOM on a molecular level are poorly understood.The plethora of compounds contained within NOM have resisted scientists' use of modern analytical techniques to unravel their molecular structures.This is mostly because of our inability to isolate individual molecular con-stituents of NOM. 6 Nevertheless, knowledge of molecular structure is crucial to the comprehension of the various roles of NOM.It will allow us to study many important processes on a molecular level e.g. the ability of phenolics 7 or Sphagnum moss molecules 8 to hinder peat bog degradation.To advance and fully develop the structure-function relationships of NOM molecules, we need to tackle the issue of how to determine the molecular composition of complex mixtures.
2][13][14][15] However, at present, users of both techniques face significant challenges.The sheer number of compounds constituting NOM prevents the interpretation of standard multidimensional NMR spectra, making it impossible to determine the structure of individual molecules.FT-ICR MS offers incredible resolution and high mass accuracy; nevertheless, methodologies are yet to be developed that will take this technique beyond molecular formulae determination.Imagina- † Electronic supplementary information (ESI) available.See DOI: 10.1039/ c6an00999a a EastCHEM School of Chemistry, University of Edinburgh, King's Buildings, David Brewster Road, Edinburgh, EH9 3FJ, UK.E-mail: dusan.uhrin@ed.ac.uk tive approaches are required that will push the boundaries of these powerful techniques to structurally characterise molecules buried in complex mixtures.
There are few examples of recent developments, such as the application of 2D HPLC-NMR 16 or MS n techniques 17 that provide some directions.We have proposed a combination of chemical modifications and isotope-filtered nD NMR spectroscopy as a way of simplifying NMR spectra of NOM to determine the structure of individual molecules. 18,19In this approach, NMR active isotopes contained within a molecular tag, are attached to a subset of NOM molecules.The structures of the tagged molecules are then investigated by using the tags as molecular spies that report on their immediate environment.We have proposed the use of 13 CH 3tags to convert the most abundant functional groups of NOM molecules, the carboxyls and hydroxyls into esters and ethers, respectively.Focusing on the phenolic compounds we have identified 33 structures/structural fragments in an operational fraction of soil organic matter, fulvic acid isolated from a Scottish peat bog 19 that was prepared according to the International Humic Substance Society (IHSS) protocols. 20ere we extend our original reports by describing a set of 13 C-filtered nD NMR experiments which lead to the structure elucidation of these molecules.Even though the application of this methodology requires access to cryoprobe instruments, preferably ≥600 MHz, we suggest that its end result can be used for profiling of phenolic compounds using routine NMR instruments.

Methods
All NMR spectra were acquired in CDCl 3 (550 µl) on a Bruker Avance III 800 MHz spectrometer equipped with a 5 mm TCI cryoprobe at 288 K to minimise gradient induced losses of signals caused by diffusion during long evolution intervals in the 3D INEPT-INADEQUATE-HSQC pulse sequence.The pulse sequences of proposed NMR experiments are shown in Fig. 2 and 4 and the general acquisition parameters are summarised in Table 1.
The methylation of the model mixture and the peat soil FA fraction has been described previously. 18The model mixture of nine hydroxyl benzoic acids contained the following molecules: 2-hydroxybenzoic acid, 3-hydroxybenzoic acid, 2,3dihydroxybenzoic acid, 2,4-dihydroxybenzoic acid, 3,5-dihydroxybenzoic acid, 2,5-dihydroxybenzoic acid, 3,4-dihydroxybenzoic acid, 3,4,5-trihydroxybenzoic acid and 4-hydroxybenzoic acid was 13 C-methylated yielding the final concentration of 1.4 mM of each compound.The NMR tubes were sealed to avoid evaporation of CDCl 3 over time.The parameters used to acquire the spectra of the model mixture were similar to those used for the methylated FA sample.The latter are detailed below.
The 2D 1 H, 13 C HSQC spectrum was recorded using the standard Bruker hsqcetgp sequence.The delays were optimised for a 1 J CH = 145 Hz.The 1 H and 13 C carrier frequencies were set to 3.8 and 55 ppm.XY-32 decoupling 21 was used during t 2 at γB 1 /2π = 3571 Hz.
The 3D INEPT-INADEQUATE-HSQC spectrum was acquired using the pulse sequence published previously. 18The experiment was optimised for n J CC = 6 Hz and 1 J CH = 153 Hz.The 1 H and 13 C carrier frequencies were set to 5.5 and 110 ppm, respectively.2 ms composite adiabatic CHIRP 13 C pulse 22 was used to refocus carbon resonances, while another 500 μs CHIRP adiabatic pulse was applied to invert carbon resonances.The acquisition in t 1 and t 2 was delayed by half of the dwell time and the ϕ 1 phase was incremented by 60°simultaneously with the t 1 incrementation.The corresponding zero and first order corrections in F 1 were 24°, −180°, respectively.Bi-level adiabatic decoupling was used during t 3 .
The 3D HMQC-HMBC spectrum was acquired using the pulse sequence shown in Fig. 2A.The delays were optimised for a 1 J CH = 145 Hz and n J CH = 6 Hz.The 1 H and 13 C carrier frequencies were set to 3.8 and 115 ppm, respectively.The acquisition in t 1 and t 2 was delayed by half of the dwell time and the ϕ 1 phase was incremented by −50°simultaneously with t 1 incrementation.The corresponding zero and first order corrections in F 1 were 153°, −180°, respectively.A 2 ms composite adiabatic CHIRP 13 C pulse was utilised for refocusing of carbon resonances while a 150 μs r-SNOB 180°pulse was used to refocus the 1 J CH 3 evolution during the HMBC step.XY-32 decoupling was used during t 3 at γB 1 /2π = 3571 Hz.
The 4D HCCH 3 spectrum was acquired using the pulse sequence published previously. 19The delays were optimised for a 1 J CHarom = 155 Hz, 3 J CC = 6 Hz and 1 J CH 3 = 145 Hz.To reduce off-resonance effects the carrier frequencies were changed at specific places.The 13 C frequency offsets were changed from 56 to 110 to 96 then back to 56 ppm, while the 1 H offsets were changed from 7.5 to 5.5 and finally to 3.8 ppm during the pulse sequence.The acquisition in t 1 , t 2 and t 3 was delayed by half of the dwell time and the ϕ 5 phase was incremented by 135°, while the ϕ 6 phase was incremented by −20°s imultaneously with t 2 and t 1 incrementation, respectively.This required the zero and first order phase corrections to be set to 90°, −180°in F 1 , 188°, −180°in F 2 and 154°, −180°in F 3 .2 ms composite adiabatic CHIRP 13 C pulses were utilised for refocusing carbon resonances.A 1 ms trim pulse was applied after the Double Pulsed Field Gradient Spin Echo (DPFGSE) with full power, while 1 ms r-SNOB 180°pulses were used for the DPFGSE.XY-32 decoupling was used during t 4 at γB 1 /2π = 3571 Hz.The 3D HcCH 3 spectrum was acquired using the pulse sequence shown in Fig. 2B.The acquisition in t 1 and t 2 was delayed by half of the dwell time and phases ϕ 1 to ϕ 4 were incremented by −20°simultaneously with t 1 incrementation.This required setting the zero and first order phase corrections to 83°, −180°in F 1 and 90°, −180°in F 2 .Other parameters were the same as for the 4D HCCH 3 experiment.
The 3D hCCH 3 spectrum was acquired using the pulse sequence shown in Fig. 2C.The acquisition in t 1 and t 2 was delayed by half of the dwell time and phase ϕ 5 was incremented by +135°, while ϕ 6 was incremented by −90°simultaneously with t 1 and t 2 incrementation, respectively.This required to set the zero and first order phase corrections to 28°, −180°in F 1 and 129°, −180°in F 2 .Other parameters were the same as for the 4D HCCH 3 experiment.
The 3D HMQC-NOESY spectrum was acquired using the pulse sequence shown in Fig. 4A.The NOESY mixing time was set to 600 ms.The carrier frequencies in 1 H and 13 C dimensions were 3.9 and 56.6 ppm, respectively.The acquisition in t 1 and t 2 was delayed by half of the dwell time, with the corresponding phases set to 90°in the DPFGSE.Two 2.8 ms 180°I BURP2 pulses were applied at 3.9 and 7.4 ppm two times at specified intervals during the NOESY mixing interval.The second pulse was created by time inverting the IBURP2 profile and phase shifting it to invert the aromatic protons.
The 3D HMQC-NOESY-TOCSY spectrum was acquired using the pulse sequence shown in Fig. 4B.The NOESY and TOCSY mixing times were 600 and 80 ms, respectively.Two hard 180°p ulses were used during the mixing time. 23In order to remove zero quantum coherences 24 a 20 ms adiabatic CHIRP pulse was applied simultaneously with 4% gradient after the DIPSI2 spin-lock using γB 1 /2π = 10 kHz.Other parameters were the same as for the 3D HMQC-NOESY experiment.
Pulse sequences of the proposed NMR experiments are given in the ESI.†

Description of 13 CH 3 O-filtered nD NMR experiments
Isotope tagging places an NMR active nucleus onto the target molecules by modifying their functional groups.Such a nucleus serves a dual purpose.It eliminates NMR signals from non-tagged molecules in the NMR spectra and reports on the chemical shifts of nearby nuclei.In this work the 13 C tagging was achieved using 13 CH 3 I as described previously. 18,19Methylation by 13 CH 3 I is non-selective and modifies both phenols and aliphatic hydroxyl groups.In addition, it also converts aromatic and aliphatic carboxyl groups into methyl esters.In order to extract the chemical shifts of the tagged molecules nuclei, several nD 13 CH 3 O-filtered NMR experiments are required.These combine multiple polarisation transfer steps, initiating with passage through the nuclei of the 13 CH 3 O group enabled by numerous J couplings of this tag (Fig. 1).For phenolic compounds, the relevant coupling constants (Table 2) were measured on a mixture of nine model compounds. 18In the following, 13 C-filtered 3D experiments utilising these tags to report on their environment are described.Their relative sensitivity and individual merits are discussed.An example is given to illustrate how best to analyse the obtained spectra to achieve molecular level characterisation of phenolic compounds.

3D INEPT-INADEQUATE-HSQC
This experiment, applied to a mixture of nine 13 CH 3 O-tagged model compounds, has been described previously, 18 therefore only a brief outline is presented here.INADEQUATE 25 is the ultimate NMR experiment for the structure determination of organic molecules, as it allows their carbon skeletons to be traced.However, its reliance on the presence of two 13 C nuclei in the same molecule lowers the sensitivity and consequently the application scope of this experiment.However, this is not an issue here, as only the n J ( 13 C, 13 C) couplings involving 100% enriched 13 CH 3 O-groups are utilised.To boost the sensitivity of the experiment further proton detection along the lines of the 2D INEPT-INADEQUATE experiment was employed. 26This experiment provides double-quantum (DQ) 13 C frequencies of coupled carbon nuclei and the 1 H chemical shifts of protons attached to carbons forming the DQ coherences.The lost 13 C chemical shift information is recovered by incorporating constant-time single quantum (SQ) chemical shift labelling 27,28 into one of the long evolution intervals of the experiment.The resulting 3D INEPT-INADEQUATE-HSQC is a very powerful NMR experiment that provides numerous chemical shifts of tagged molecules. 18Specifically, carbon chemical shifts of ipso, ortho and occasionally meta carbons relative to the methoxy groups, proton chemical shifts of the ortho protons and the proton and carbon chemical shifts of the methoxy groups.A partial 3D spectrum of 13 C-methylated peat soil fulvic acid obtained from this experiment is shown in Fig. 3, together with the spectra of the rest of the experiments discussed below.

3D HMQC-HMBC
The 3D HMQC-HMBC experiment (Fig. 2A) correlates the chemical shifts of the quaternary carbon carrying a 13 CH 3 Ogroup with the chemical shifts of the methoxy group nuclei.The pulse sequence of the 3D HMQC-HMBC starts with a HMQC part, where the 13 C chemical shifts of the methoxy carbons are labelled (t 1 ), followed by the creation of an antiphase methyl proton magnetisation with respect to the quaternary ipso carbons.Once this is converted into MQ coherences, the chemical shifts of the quaternary carbons are labelled in the HMBC part of the sequence, followed by a back conversion to proton magnetisation prior to the detection under 13 C decoupling.The evolution of 1 J CH 3 couplings of fully labelled methoxy groups during the long evolution intervals of the HMBC part of the pulse sequence need to be suppressed.This is achieved by (i) setting the HMQC Δ 1 interval to the average 1/2 1 J CH 3 value followed by conversion by a 90°1 3 C pulse of any imperfectly refocused signals at the end of the HMQC sequence into MQ coherences; (ii) selective inversion of methoxy carbons mid-way through the defocusing and refocusing periods of the HMBC part of the pulse sequence.This is a far more effective way to eliminate the overall 1 J CH 3 evolution than setting the long-range evolution intervals to an even multiple of 1/2 1 J CH 3 using an average value of 1 J CH 3 .

3D/4D HCCH 3
The HCCH 3 experiments provide chemical shifts of the CH moieties ortho to a methoxy group as well as the proton and carbon chemical shifts of this methoxy group.Although the 3D INEPT-INADEQUATE-HSQC experiment 18 also provides in principle the chemical shifts of all HCCH 3 nuclei in a single experiment, the S/N of aromatic protons in the INADEQUATE spectra is noticeably lower than that of the methoxy protons.This is because the three protons of the methoxy group show as singlets in NMR spectra, while in general, aromatic protons are presented as multiplets in the INADEQUATE spectra.The HCCH 3 experiments record the aromatic proton chemical shifts in the indirectly detected dimension, thus alleviating this problem.Chemical shift labelling of all nuclei along the HCCH 3 polarisation transfer pathway naturally leads to a 4D experiment.The 4D HCCH 3 experiment provides an excellent insight into the substitution patterns of phenolic molecules by separating signals into four dimensions, as demonstrated previously. 19A limiting factor, common to all 4D NMR experiments, is that adding a 4 th dimension reduces the sensitivity and compromises the achievable resolution relative to 3D experiments.This issue is addressed here by two complementary 3D experiments, 3D HcCH 3 and 3D hCCH 3 , which in addition to the nuclei of the methoxy group, only label chemical shifts of either the aromatic protons or carbons.
The 3D HCCH 3 experiments (Fig. 2B and C) start with the excitation of aromatic protons and end with the detection of methyl protons.The initial selection of the aromatic protons is achieved by a DPFGSE 29 containing r-SNOB pulses, 30 which were chosen due to their near 'top-hat' profile.The DPFGSE eliminates phase errors across the inverted frequency range and is therefore preferable to a shorter SPFGSE sequence.Both 3D HCCH 3 experiments follow the magnetisation pathway of the 4D HCCH 3 , but differ in whether the chemical shits of the aromatic protons or carbons are labelled.In the HcCH 3 experiment, the chemical shifts of the aromatic protons are labelled during the t 1 period before an INEPT transfer to the aromatic carbons, while in the hCCH 3 experiment the chemical shifts of the aromatic carbons are labelled after the transfer.Once the magnetisation is on the aromatic carbons, the spin system evolves to create an antiphase magnetisation with regard to the methoxy carbon while refocusing the 1 J CHarom couplings.The magnetisation is then passed to the methoxy carbons for chemical shift labelling during t 2 , before transfer to the methyl protons for detection during t 3 .The carbon decoupling over a narrow range of methyl carbon resonances with reduced power allows longer acquisition and thus very good signal resolution.

3D HMQC-NOESY
The 3D HMQC-NOESY experiment is an alternative for obtaining the chemical shifts of protons ortho to the methoxy group.Unlike in the experiments described above, here the magnetisation does not pass through the natural abundance 13 C atoms thus avoiding the inevitable sensitivity loss of heterocorrelated experiments.The observed efficiency of the NOE transfer between methoxy and aromatic ring protons investigated by a 1D NOESY experiments on the mixture of nine methylated compounds was around 4% (data not shown).This exceeds the natural abundance of 13 C (1%), thus providing higher sensitivity, albeit with reduced information content by not including the 13 C chemical shifts.Nevertheless, the 3D HMQC-NOESY experiment proved very useful, in part because it also yields the proton-proton couplings, thus providing additional structural information.
The 3D HMQC-NOESY experiment (Fig. 4A) starts with a DPFGSE, selecting methoxy protons, followed by chemical shift labelling of methoxy carbons during the t 2 period of the refocused HMQC pulse sequence.Chemical shifts of the methoxy protons are labelled during the t 1 period and their magnetisation is then positioned along the z-axis, by the action of a 90°1H pulse, for NOE transfer.An efficient NOE transfer for small molecules requires long NOE build-up times (>0.5 s), during which considerable auto relaxation occurs.This leads to the appearance of axial peaks at the zero F 1 frequency and an overall increase in the t 1 noise.These effects can be reduced by inverting resonances during the NOE mixing time. 29In addition, using selective rather than Fig. 2 Pulse sequences of 13 C-filtered 3D NMR experiments.The nuclei whose chemical shift are sampled are highlighted on the structures shown on the right.The green, black, purple, red and cyan boldfaced text represent the methyl protons, methyl carbons, aromatic protons, aromatic carbons and carboxyl carbons, respectively.The thin and thick rectangles represent 90°and 180°pulses, respectively.Unless stated otherwise the pulses were applied from the x-axis.(A) 3D HMQC-HMBC.Δ 1 = 1/2 1 J CH 3 , δ = t 180 ( 1 H), Δ 2 = 1/4 n J CH − t 180 ( 13 C SNOB )/2, δ 1 = δ G + 0.5t 180 ( 13 C adia ) + 0.5t 180 ( 1 H), where δ G is the gradient duration plus recovery time, while t 180 is the duration of 180°pulse.The following phase cycle was used φ 1 = 2x, 2(−x); φ 2 = x, −x; φ 3 = 4(x), 4(−x), ψ = x, −x, −x, x, −x, x, x, −x. 1 ms gradient pulses were applied at the following strengths: G 1 = 10%, G 2 = 80% and G 3 = 47.85%.Open and filled gradients were applied in odd and even t 2 increments, respectively.Phase φ 1 was incremented according to the States-TPPI protocol.Quadrature detection in t 2 was achieved using the echo, anti-echo protocol with phases φ 2 and ψ incremented by 180°simultaneously with t 2 incrementation.The carrier frequencies were changed at specified places, indicated by the balloon topped arrows: ① = middle of methyl carbons (55 ppm) and ② = middle of aromatic and methyl carbons (115 ppm).A CHIRP refocusing pulse was applied at the end of the t 2 period and two r-SNOB pulses were used to invert the methoxy carbon resonances in the middle of the 2Δ 2 intervals.The 90°1 3 C pulse labelled with an asterisk converts any non-refocused methyl proton magnetisation into MQ coherences.(B) 3D HcCH 3 .1 ms r-SNOB pulses are applied during the initial DPFGSE to select aromatic 1 H resonances.The following parameters were used: Δ 1 = 1/4 1 JCH arom , Δ 2 = 1/4 3 J CC , Δ 3 = 0.0915/ 1 J CH 3 , Δ 4 = 1/4 1 J CH 3 .The following phase cycle was used: φ 1 = x, φ 2 = 4(x), 4(y), φ 3 = 2(x), 2(y), φ 4 = y, φ 5 = x, −x, ψ = x, −x, −x, x, −x, x, x, −x.Phases φ 1 , φ 4 (and φ 5 ) were incremented according to the States-TPPI protocol.When phase shifting was applied in t 1 , (and t 2 ), the phases φ 1 , φ 2 , φ 3 and φ 4 (and φ 5 ) were phase shifted as stated in the Methods section. 1 (open) and 0.6 ms (dashed) gradient pulses were applied at the following strengths: nonselective pulses to achieve this, the NOE transfer is restricted between the groups of inverted protons. 31When only aromatic and methoxy protons are inverted, leakage of magnetisation from methoxy to other aliphatic protons is therefore eliminated.As shown by Stott et al. placing one 180°pulse in the middle of the mixing time interval is less efficient, compared to two pulses spaced appropriately.It is best to use two 180°p ulses, one at ∼0.25τ mix and the other at ∼0.75τ mix to get efficient artefact suppression.Band selective IBURP2 pulses 32 were used here because of their 'top hat' inversion profile.In order to avoid the so-called 'close encounter' effects between soft pulses when two neighbouring regions are inverted simulteneously, 33 the pulses inverting the methoxy and aromatic resonances were applied consecutively.There is no need for carbon decoupling in this pulse sequence, which allowed longer acquisition times to be used, providing accurate readouts of the proton-proton coupling constants of aromatic protons.

3D HMQC-NOESY-TOCSY
It is straightforward to extent the 3D HMQC-NOESY pulse sequence to construct a 3D HMQC-NOESY-TOCSY experiment (Fig. 4B).This involves appending the NOESY pulse sequence with a DIPSI-2 spin-lock 34 containing a ZQ interference removal element 24 applied after the TOCSY spinlock.In this experiment the band selective inversion pulses applied during the NOESY mixing time were replaced by non-selective pulses, as the molecules in the studied sample have no prospect of NOE transfer from methoxy to other aliphatic protons.This modification can also be applied to the 3D HMQC-NOESY experiment.

Relative sensitivity of the isotope-filtered nD experiments
The discussed pulse sequences produce high quality 3D spectra, nevertheless, as there is a partial redundancy in the obtained chemical shifts, it is useful to compare their performance.The S/N ratios achievable by individual experiments were therefore investigated by acquiring 3D/4D NMR spectra of a model mixture of nine methylated aromatic compounds.As much as possible, similar experimental conditions were used for this comparison.In particular, the digital resolution in the indirectly detected dimensions was kept the same.In order to achieve this, while at the same time accommodating different spectral widths required by individual experiments, different overall acquisitions times had to be used.These differences were eliminated by scaling the measured S/N ratios relative to the S/N of a 3D INEPT-INADEQUATE-HSQC spectrum obtained in a 14 hour long experiment. 18As all experiments provide 1 H, 13 C correlations of the methoxy groups, 1D traces were extracted from the HSQC/HMQC planes at individual methoxy carbon chemical shifts.The S/N in these traces was measured for methoxy protons except for the 3D HMQC-NOESY and 3D HMQC-NOESY-TOCSY, where the S/N was determined for aromatic protons.The NOESY peaks were not included when evaluating the TOCSY spectra.The S/N ratio of the 3D INEPT-INADEQUATE-HSQC spectrum was measured in both aromatic and methoxy regions.The S/N of the 4D HCCH 3 spectrum was obtained by extracting 3D cubes at individual methoxy 13 C chemical shifts, followed by the extraction of 1D slices from the 2D HSQC planes, as done for the 3D experiments.The obtained data are presented in two graphs (Fig. 5A  and B) for clarity in the form of bar charts, where the count represents the number of cross-peaks with a particular S/N ratio within a 0-3000 range, increasing in steps of 300.From the comparison of the S/N ratios obtained, it is clear that the discussed experiments are not vastly different in terms of their sensitivity.As expected, the 3D HCCH 3 experiments show better S/N compared to the 4D HCCH 3 experiment.The INEPT-INADEQUATE-HSQC has the lowest sensitivity, which again is to be expected due to the inclusion of DQ coherences in the polarisation transfer pathway of this experiment.The 3D HMQC-NOESY experiment is the most sensitive, followed by the 3D HMQC-HMBC experiment.
Practical guide to the interpretation of 13  thus be regarded as the most important experiment for the analysis of phenolic compounds, as it identifies all accessible carbon chemical shifts, including the ortho carbons (regardless if protonated or not).However, as can be seen from the S/N comparison presented in Fig. 5, the sensitivity of the 3D HMQC-HMBC and 3D HCCH 3 experiments is higher than that of the 3D INEPT-INADEQUATE-HSQC.The 3D HMQC-NOESY-TOCSY yields chemical shifts of additional aromatic protons and therefore provides valuable information extending the characterised structures.While the interpretation of 13 C-filtered nD spectra of 13 CH 3 O-tagged NOM can be approached in several different ways, a particular procedure that efficiently pieces together molecular structures of phenolic molecules is summarised next.This is illustrated on the structure determination of a compound carrying a methoxy group resonating at δ 13 C/ 1 H of 55.41/3.87ppm.The process can be divided into the following six steps: 1. Identify methoxy cross-peaks in the 2D 1 H, 13 C HSQC spectrum.
2. Using the 3D HMQC-HMBC spectrum, determine the chemical shifts of ipso carbons associated with individual methoxy cross-peaks of the 1 H, 13 C HSQC spectrum.
3. Using the HCCH 3 spectra, identify chemical shifts of protonated aromatic CH moieties ortho to the methoxy groups (if there is no corresponding cross-peaks in 4D or 3D HCCH 3 spectra, then the methoxy group is either surrounded by two quaternary carbons or it is not attached to an aromatic molecule).
4. Verify the assignments obtained in points 2 and 3 by inspecting the 3D INEPT-INADEQUATE-HSQC spectrum.Identify any additional cross-peaks, especially belonging to the quaternary carbons.Extra (weak) cross-peaks may belong to carbons meta to the methoxy group.
5. Extend the 1 H spin systems by inspecting the 3D HMQC-NOESY and HMQC-NOESY-TOCSY spectra.As a starting point, the 55.41 ppm F 1 F 3 2D plane of the 3D HMQC-HMBC spectrum is inspected.The chemical shift of the ipso carbon at 163.32 ppm (Fig. 6) confirms that this methoxy group is joined to an aromatic ring.
Next, any ortho CH moieties are identified by inspecting the 4D/3D HCCH 3 spectra.A 3D cube extracted from the 4D HCCH 3 spectrum through the methoxy carbon at 55.41 ppm is shown in Fig. 7A.The 55.41 ppm F 1 F 3 2D planes from the 3D hCCH 3 and HcCH 3 spectra (Fig. 7B and C) contain two crosspeaks; the one at the correct methoxy proton chemical shift of 3.87 ppm is indicated by the dotted lines.This analysis associates the HC atoms at 6.93/113.42ppm with the investigated methoxy group.As only one CH pair was identified, the molecule is either symmetrical around the methoxy group, or the other ortho position contains a quaternary carbon.
Fig. 8 shows the 55.41 ppm F 1 F 3 2D plane extracted from the 3D INEPT-INADEQUATE-HSQC spectrum.Its methoxy region (Fig. 8A) contains three cross-peaks at the methoxy proton chemical shift of 3.87 ppm.The two most intense signals with the 13 C chemical shifts of 113.42 and 163.32 ppm have already been identified by previous experiments as belonging to ortho and ipso carbons, respectively.An additional cross-peak at 131.68 ppm is much weaker than that at 113.42 ppm and therefore does not represent the second ortho proton.This analysis confirms that the corresponding fragment is likely symmetrical around the methoxy group.The aromatic region of the plane (Fig. 8B) shows only one aromatic proton resonating at 6.93 ppm, reflecting the symmetry of the molecule.Its signal appears as a doublet ( J = 8.5 Hz), indicating that the meta carbons relative to the methoxy group are protonated.
Fig. 9 shows the 55.41 ppm F 1 F 3 2D planes extracted from the 3D HMQC-NOESY and 3D HMQC-NOESY-TOCSY spectra.The signal at 6.93 ppm ( J = 8.5 Hz) belongs to protons ortho to the methoxy group.The only other signal identified in the 3D  HMQC-NOESY-TOCSY spectrum is a doublet at 7.99 ppm, confirming that this molecule is symmetrical.
In summary, information gathered thus far consistently indicates that the analysed molecule is a para di-substituted benzene with a methoxy group at one end of the molecule.
When the obtained chemical shifts of the investigated fragment were compared with the database chemical shifts of the methoxy benzene 35,36 it became clear that the substituent para to the methoxy group has a large ortho deshielding effect, as indicated by the increased proton chemical shift (7.99 ppm).This is typical for an electron-withdrawing group such as a carboxylic acid or an ester.To examine if the para position carries a methyl ester group, the carboxyl region of the 3D HMQC-HMBC was inspected (data not shown).Considering the database chemical shifts 35 of the ester group of methyl p-methoxy benzoate (166.82,51.82 and 3.86 ppm), the 51.89 ppm F 1 F 3 2D plane of the 3D INEPT-INADEQUATE-HSQC spectrum was inspected.It showed three cross-peaks at SQ 13 C chemical shifts of 166.91, 122.61 and 131.81 ppm, respectively (Fig. 8C).The first two intense cross-peaks belong to the carbonyl of the ester group and the ipso quaternary carbon of the aromatic ring, respectively.The remaining weak correlation at 131.84 ppm involves the aromatic carbons ortho to the methyl ester group.This correlation was also identified as a weak cross-peak in the methoxy ether region of the 3D INEPT-INADEQUATE-HSQC spectrum (Fig. 8A), which indicates that it belongs to the same molecule.A comparison of the database and the experimental chemical shifts (Fig. 10) showed an excellent agreement, indicating that the inspected molecule is methyl p-methoxybenzoate.It should be noted that the 4 J CH 3 ,C carbon-carbon couplings are small and correlations such as involving 13 C at 131.81 ppm mediated by these couplings will only be observed for the most abundant compounds.Nevertheless, even in the absence of this information, the structure of this compound would have been identified by analysing the chemical shifts associated with the methoxy group.
Analogous analysis starting with the methoxy cross-peaks present in the 1 H, 13 C HSQC spectrum of the methylated FA sample (Fig. 11) identified 33 structures/structural fragments as reported previously 19 (Fig. 1S of the ESI †).This spectrum shows a typical diagonal trend between the proton and carbon chemical shifts.Analysis of the identified compounds/structural fragments revealed that the cross-peaks can be divided into four groups labelled (a) to (d) as shown in Fig. 11.Region (a) contains resonances of methoxy groups of p-disubstituted benzene rings, while region (b) contains methoxy groups of 3,4-disubstituted rings, with an exception of cross-peak 25 that belongs to the one of the two 3,5 methoxy groups of a 3,5-disubstituted molecule 25/28.Region (c) is the most heterogeneous in terms of the structures of its molecules.Cross-peak 28 belongs to the other methoxy group of molecule 25/28, while cross-peaks 14 and 19 belong to compounds containing three methoxy groups.Molecule 19 only contains one 13 Clabelled and two unlabelled methoxy groups, while in molecule 14 all OH groups at positions 3, 4, 5 were 13 C-methylated (the C4 OMe group is not shown in Fig. 11, as its carbon resonates at 60.89 ppm).The remaining molecules of this group had only one methoxy group located next to an unsubstituted CH site on otherwise highly substituted aromatic rings containing 2 to 3 additional substituents.Region (d) contains signals of methoxy groups located next to an unsubstituted CH site on otherwise highly substituted rings containing 3 to 4 additional non-oxygen containing substituents, i.e. similar to some molecules of the previous group, albeit with a different character of substituents causing further deshielding of methoxy resonances.
This grouping of cross-peaks enables a quick assessment of compound types present in 13 C methylated mixtures without the need for a full analysis along the lines presented here.Based only on a 2D 1 H, 13 C HSQC spectra of 13 C methylated mixtures, which can be acquired on a standard NMR spectrometer equipped with a room temperature probe, this classification will improve when more mixtures are fully characterised.

Conclusions
A series of isotope-filtered nD NMR experiments designed for the structure elucidation of 13 CH 3 O-labelled mixture was presented.Their application to a peat FA sample illustrates how structures of aromatic compounds contained in a complex mixture can be obtained.Analysis of these experiments indicates that in an ideal case only two experiments are required for full molecular characterisation, 3D INEPT-INADEQUATE-HSQC and the HMQC-NOESY-TOCSY.However, for less abundant compounds yielding weaker signals, it is recommended to acquire a whole suite of experiments to overcome the sensitivity and potential overlap issues.It can be envisaged that similar experiments can be designed for different compound classes and/or different tags.Although illustrated here on a NOM sample, this approach is not limited to studies of the mixtures found in the environment and can applied other complex mixtures, such as food and beverages, bio fluids, plant extracts, or lignin degradation products.

Fig. 1
Fig. 1 Possible polarisation transfer pathways of 13 CH 3 O tagged phenolic molecules.Note: symmetrical pathways are depicted only once.The 4 J CC indicated by dashed lines are ∼1 Hz.
CH 3 O-filtered nD NMR spectra of tagged complex mixtures Only the 3D INEPT-INADEQUATE-HSQC and 3D HMQC-HMBC spectra provide chemical shifts of quaternary carbons.The 3D INEPT-INADEQUATE-HSQC spectra yield the chemical shifts of CH atoms ortho to the methoxy groups, as do the HCCH 3type experiments.The 3D INEPT-INADEQUATE-HSQC could

6 .
Examine cross-peaks in the 3D HMQC-HMBC and 3D INEPT-INADEQUATE-HSQC spectra associated with the methoxy ester groups.Cross reference with the weak crosspeaks associated with phenolic methoxy groups.

Fig. 6
Fig.6The 55.41 ppm F 1 F 3 2D plane extracted from the 3D HMQC-HMBC spectrum of13 C methylated FA.The cross-peak with the 3.87 ppm methoxy proton chemical shift is highlighted.

Fig. 7 (
Fig. 7 (A) A 3D cube extracted from the 4D HCCH 3 spectrum of 13 C methylated FA at the methoxy carbon chemical shift of 55.41 ppm.The 55.41 ppm F 1 F 3 2D planes extracted from (B) 3D hCCH 3 and (C) 3D HcCH 3 spectra.Dotted lines highlight the cross-peaks for the methoxy 1 H at 3.87 ppm.

Fig. 8
Fig. 8 The (A) methoxy and (B) aromatic regions of the 55.41 ppm F 1 F 3 2D plane and (C) methoxy region of the 51.89 ppm F 1 F 3 2D plane extracted from the 3D INEPT-INADEQUATE-HSQC spectrum of 13 C methylated FA.The cross-peaks are labeled using 13 C chemical shifts, not the DQ coherences.

Fig. 9
Fig. 9 The 55.41 ppm F 1 F 3 2D planes extracted from (A) the 3D HMQC-NOESY and (B) 3D HMQC-NOESY-TOCSY spectra of 13 C methylated FA.The dotted line indicates the methoxy proton chemical shift at 3.87 ppm.The position of protons relative to the methoxy group are indicated.

Fig. 10
Fig. 10 Completed structure of molecule 23 based on the chemical shifts of the highlighted nuclei obtained from the isotope-filtered nD NMR spectra.The experimental (black, top) and database (grey, bottom) chemical shifts are compared.

Fig. 11
Fig.11Methoxy region of 2D 1 H,13 C HSQC spectrum of13 C methylated FA.The cross-peak labels correspond to the compounds given in Fig.1Sof the ESI.† The dashed rectangles delineate regions discussed in the text.

Table 1
13quisition parameters of nD NMR spectra of13C methylated FA

Table 2
Typical coupling constants of phenolic compounds a