Article
Subscriber access provided by SUNY DOWNSTATE
Direct Identification of Diverse Alicyclic Terpenoids in Suwannee River Fulvic Acid Neal Arakawa, and Lihini I. Aluwihare Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/es5055176 • Publication Date (Web): 13 Mar 2015 Downloaded from http://pubs.acs.org on March 25, 2015
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 28
Environmental Science & Technology
1
Direct Identification of Diverse Alicyclic Terpenoids
2
in Suwannee River Fulvic Acid
3
Neal Arakawa*. Lihini Aluwihare
4
Scripps Institution of Oceanography, University of California San Diego, 9500 Gilman Drive, La
5
Jolla, CA 92093-0244.
6
The chemical complexity of natural dissolved organic matter (NOMDOM) obstructs our ability
7
to definitively recover source compounds from within NOMDOM, an objective which has the
8
capacity to alter our understanding of carbon sequestration on a global scale. To advance
9
compositional studies of NOM DOM we have applied a previously published reduction method
10
to an environmental standard, Suwannee River Fulvic Acid (SRFA). The reduction products,
11
comprising 12% of the pre-reduced carbon, were then separated by comprehensive two-
12
dimensional gas chromatography time-of-flight mass spectrometry (GCxGC-TOF-MS). Results
13
indicate that the majority of observed reduced compounds corresponded to alicyclic
14
hydrocarbons in the size range C10 to C17. Cyclic terpenoids are the only biomolecule class with
15
contiguous, alicyclic carbon backbones of this size. These terpenoid reduction products contain
16
series offset by CH2 and exhibit great isomeric diversity, features previously inferred from
17
ultrahigh resolution mass spectrometry and NMR studies of unreduced SRFA. Reduction of
18
Taxoddinium leaf litter as a source material to SRFA confirmed the prevalence of terpenoids in
19
SRFA and provided insight into the parent compounds that must be diagenetically modified on
1
ACS Paragon Plus Environment
Environmental Science & Technology
Page 2 of 28
20
relatively short timescales. These data corroborate several recent studies that suggest alicyclic
21
hydrocarbons to be important components of longer-lived DNOM.
22
Abstract Art
23 24
Introduction
25
Dissolved organic carbon (DOC) from aquatic environments is most easily divided into two
26
classes of material; compounds of identifiable structure and biochemical origin (e.g., lignin
27
phenols, sugars, lipids), and compounds with unknown structure (e.g., typically termed humic
28
substances [humic acid, fulvic acid, humin], protokerogen). The majority of Dissolved organic
29
matter (DOM) represents DOC falls into the latter group, which makes it a quantitatively
30
significant reservoir of the global carbon cycle. 1, 2 Yet, its provenance remains elusive, as much
31
of it has been immune to degradative compositional interrogation. Recent studies employing
32
nuclear magnetic resonance spectroscopy3-5 (NMR) and Fourier transform ion cyclotron
33
resonance spectrometry6,7 (FT-ICR-MS) have implicated linear and alicyclic terpenoids as major
34
precursors of terrestrial and marine DOM. Sources of terpenoids and their biosynthetic pathways
35
are well known, and so, direct identification of the relevant structures and their diversity could
36
hold significant promise for linking DOM composition to source.8
2
ACS Paragon Plus Environment
Formatted: Indent: First line: 0"
Page 3 of 28
Environmental Science & Technology
37
Some of the earliest reports of terpenoid precursors in DOM are from studies of Suwannee River
38
Fulvic Acid (SRFA), a widely used environmental DOM standard. For this reason, numerous
39
attempts have been made to compositionally deconstruct humic substances (used loosely here to
40
refer to DOM that is poorly characterized). For such studies, Suwannee River Fulvic Acid
41
(SRFA) has been heavily used as an environmental natural organic matter (NOM) standard in the
42
development of analytical techniques.
43
The parent material of SRFA is organic matter leached from leaves or soil.3 However,
44
biological and abiotic processes have transformed these compounds into an extremely
45
chemically diverse mixture with a broad molecular size range. Humic substances and SRFA
46
specifically, have received decades of attention.4,
47
application of ultra-high resolution Fourier transform ion cyclotron resonance mass spectrometry
48
(FT-ICR-MS) to detect thousands of compounds, many of which are assigned a specific
49
elemental formula.6, 7 In addition, 13C and multi-dimensional NMR techniques have been used to
50
provide a broad compositional context for SRFA including the identification of individual
51
chemical structures. The sum of the body of current work indicates that SRFA is derived from a
52
mixture of lignin, tannin, and terpenoid (including sterol) precursors. 8-113-5, 9 Yet, quantitatively
53
little lignin, and no definitive tannin or terpenoid structures, have been isolated from humic
54
substances12 substances10
5
More recent studies have focused on the
55
The recovery of specific structures from SRFA and other natural organic matter (NOM) would
56
not only confirm source material identity, but also has the potential to inform our understanding
57
of the processes that alter DOM on the molecular level. As with previous studies we use SRFA
58
as a model for DOM and test a method that is specifically designed to access hydrocarbon
59
backbone structures.
3
ACS Paragon Plus Environment
Environmental Science & Technology
60
Page 4 of 28
Although fractionation of SRFA via liquid chromatography is improving, isolation of
61
individual compounds is currently unattainable.13
11
62
separating capabilitieson, but much of SRFA is likely too large and polar to be analyzed
63
natively.6, 7 However, chemical degradation, such as chemical reduction, is a viable method for
64
making SRFA more GC-amenable. Previous studies employing an n-Butylsilane reduction
65
method reported the separation of various products from SRFA, but much of the mixture was left
66
unresolved.14, 1512,13 In this study, the reduction method has been adapted and combined with
67
comprehensive gas chromatography (GCxGC) coupled to time-of-flight mass spectrometry
68
(TOF-MS), to provide additional molecular-level data on the composition of SRFA.
69
GCxGC technique uses an additional column to provide unsurpassed chromatographic
70
resolution, especially when applied to hydrocarbons.16
71
presence of cyclic terpenoid derived carbon backbones compounds in SRFA, and demonstrates
72
both the unique capability of this approach and its potential for further use.. The unique
73
capability of this approach complements the existing suite of degradation techniques and can
74
provide a more comprehensive molecular level understanding of DOM source and fate.
75
Materials and Methods
76
Materials. All solvents (ACS grade or better), model compounds and hydrocarbon standards
77
were purchased commercially, and used as supplied. Suwannee River Fulvic Acid (2S101F) was
78
purchased from the International Humic Substances Society.
79
Tris(pentafluorophenyl)borane (B(C6F5)3) were purchased from Sigma Aldrich, and stored in an
80
Argon filled dessiccator when not in use. Dichloromethane (CH2Cl2) was stored over activated
81
molecular sieves (3 Å, 10% wv). Freshly fallen Lleaf litter (based on intactness of leaves) was
Gas chromatography (GC) offers better
14
This
Data presented here confirms the
n-Butylsilane (n-BS) and
4
ACS Paragon Plus Environment
Page 5 of 28
Environmental Science & Technology
82
collected from Tifton, GA, and verified to be of the genus Taxodinium. Leaves were freeze-
83
dried for 1 week, and ground using a mortar and pestle before reduction.
84
Reduction Procedure. The reduction procedure was modeled after previous studies.14
85
model compound reductions, compounds (~1 mg of each) were weighed (balance accurate to 0.1
86
mg) into a single, flame dried 2 mL vial equipped with a stir bar. Under argon atmosphere, 10
87
mg. of B(C6F5)3 was added (100 μL of 100 mg B(C6F5)3/1 mL CH2Cl2 solution). Immediately
88
after addition, 100 μL of n-BS was added, also under argon. The mixture was then allowed to
89
stir-overnight. Following reduction, the samples were treated by the following procedure to
90
remove excess catalyst and siloxane by-products. Subsequently 50 μL of the pot was syringed
91
into 250 μL of conc. potassium hydroxide in methanol (KOH/MeOH). After 1 minute, the
92
mixture was extracted with pentane (3x100 µuL). The pentane fractions were collected, washed
93
with H2O, and dried over sodium sulfate. The organic fraction was then dried under N2 to 100
94
μL for GC-FID/MS and GCxGC-TOF-MS analysis. The same procedure was used for the
95
reduction of SRFA and leaf litter, with 10 mg. of starting material for each. All reductions were
96
carried out in duplicate. While SRFA was reduced to dissolution, particulates were still observed
97
in the leaf litter reduction were still observed.
12
For
98
Reaction conditions were selected based on results from model compound reductions. The
99
amount of B(C6F5)3 used was found to be the most critical parameter affecting the reduction
100
yield of model compounds.
Chart 1 lists model compounds and their observed reduction
101
products. In a previous work with reduction of 1a (Chart 1, 10 mg), which was the most
102
chemically resistant of the model compounds tested, it was observed that increasing the amount
103
of B(C6F5)3 to 10 mg also increased reduction yields, whereas increasing n-BS did not. The
104
ratio of 10 mg B(C6F5)3 to 10 mg (or ~5 mg organic carbon) starting material was therefore
5
ACS Paragon Plus Environment
Environmental Science & Technology
Page 6 of 28
105
selected to ensure complete reduction of oxygen containing groups in model compounds as well
106
as in environmental samples.
107
B(C6F5)3 addition, as environmental samples are extremely polydisperse and therefore make the
108
use of molar equivalents impractical. This quantity of catalyst gave the most predictable and
109
consistent results across all model compounds and SRFA.
110
Gas Chromatography-Flame Ionization Detection/Mass Spectrometry (GC-MS/FID). An
111
Agilent 7890A gas chromatograph coupled to an Agilent 5975C quadrupole mass spectrometer
112
was used for chromatographic separation and compounds were detected simultaneously with
113
Electron Ionization MS and FID, enabled by a splitter. Splitless injection with 1 μL of analyte
114
was performed onto a 5% phenyl poly(dimethylsiloxane) column (J&W 123-5731DB-ht, 30 m,
115
320 μm i.d., 0.1 μm film) and separated using a temperature program from 50 °C ( hold time 0.5
116
min) to 140 °C (at 10 °C/min, hold time 0 min) to 320 °C (at 15 °C/min, hold time 10 min).
117
Helium was the carrier gas with a constant flow of 1.8 mL/min, the FID was operated at 310°C
118
and the mass spectrometer conditions included 70 eV ionization with scanning from 50-750 m/z.
119
The transfer line was set at 280 °C. Quantification of reduced model compounds and reduced
120
SRFA was accomplished using a hydrocarbon standard (decahydronapthalene) and assuming a
121
constant FID response factor across hydrocarbons. When available, reduced model compound
122
structures were identifiedconfirmed by comparison to the NIST Mass Spectral Standard
123
Reference Database (2005). For SRFA, manual integration was used to remove peaks associated
124
with catalyst degradation and residual siloxanes, and sample MS are located in Supporting
125
Information.
126
Comprehensive Two-Dimensional Gas Chromatography (GCxGC-TOF-MS).
127
Pegasus 4D GCxGC-Time-of-Flight Mass Spectrometer (GCxGC-TOF-MS) was used for
Quantity is the most practical measure for both n-BS and
A LECO
6
ACS Paragon Plus Environment
Page 7 of 28
Environmental Science & Technology
128
comprehensive GC analysis.15 The term GCxGC refers to the use of two orthogonal columns,
129
that is, each column has a different chemical selectivity which were selected to increase the
130
separation of aliphatic from aromatic hydrocarbons. Compounds are separated primarily by
131
volatility in the first column and polarity in the second. The analysis is comprehensive because
132
all of the effluent from the first column is cryofocused and transferred onto the second column.
133
Effluent from the second column is then analyzed by the TOF-MS, which benefits from high
134
spectral acquisition speed (50 Hz-500 Hz). Finally, the data is compiled into a two-dimensional
135
chromatogram that is visualized and processed by ChromaTOF® software. Both columns are
136
housed within an Agilent 7890A gas chromatograph. The splitless inlet temperature was set at
137
300 °C.
138
polysiloxane column (Restek Rxi-17Si, 30 m length, 0.2550mm i.d., 0.25 μm film thickness).
139
The column was programmed to remain isothermal at 40°C for 1 min, and ramped to 315 °C at 3
140
°C per minute. The modulator temperature was offset by +30 °C to the primary oven. The
141
secondary oven (within the GC) housed the second dimension non polar Crossbond® dimethyl
142
polysiloxane column (Restek Rxi-1, 1.58m length, 0.250mm i.d., 0.25μm film thickness). The
143
secondary oven temperature was offset by +25°C to the primary oven. The modulation period
144
was 2.5 s, with a hot pulse time of 1.05 s and a cool time of 0.2 s. The carrier gas was Helium at
145
a constant flow of 1.5 mL/min. The acquisition delay on the TOF-MS was set to 160 s. The
146
acquisition rate was set to 50 Hz, with a range of 5-1000 m/z. The transfer line was set at 280
147
°C. Electron Ionization Impact was run at 70 eV. Again, reduced model compound structures
148
were confirmed identified by comparison to the NIST Mass Spectral Standard Reference
149
Database (2011).
150
Results and Discussion
The first dimension column was a semi-polar Crossbond® diphenyl dimethyl
7
ACS Paragon Plus Environment
Environmental Science & Technology
Page 8 of 28
151
Model Compound Analysis. The reduction process was applied to a collection of chemically
152
distinct model compounds with functional groups that are expected to also be present in SRFA.
153
This was done to verify that the products were both consistent with previously published
154
literature and consistently reproducible. We observed that the catalyst (B(C 6F5)3) and by-product
155
siloxane polymers often coelute in high concentrations with analytes of interest, impacting
156
accurate GC-FID quantification and complicating GCxGC analysis.
157
KOH/MeOH to remove these by products, which resulted in cleaner chromatograms. However,
158
this method also removes reduction products which are polysilylated, such as phenols (the
159
reduction products of lignins and tannins) and some hydrosilylated (e.g. 2g, 2h, Chart 1)
160
compounds.
161
the topic of this study. Below, we systematically describe the effect of chemical reduction on
162
different functional groups, as supported by selected substrates and their reduced products.
163
Chart 1. Reduction of model compounds with n-butylsilane/ B(C6F5)3 Entry
We used concentrated
Recovery of phenols can be achieved by acid hydrolysis of silyl ethers, but are not
Substrate (1)
Product(s) (2)
GC Yield
1
58%
2
78%
3
65%
4
65%
5
72%
8
ACS Paragon Plus Environment
Page 9 of 28
Environmental Science & Technology
6
No observed product
N.A.
7
*25%
8
*25%
9
51%
10
34%
164 165
*Products 2g and 2h coelute. GC yield calculated as average
166 167
The reduction process was adapted from previous work,14 , 1512,13 and the reaction mechanism is
168
reproduced in Supporting Information. All of the carboxylic acids in model compounds (1a, 1e,
169
1i, and 1j) reduced completely to methyl groups, as predicted. The structure of 1a is likely to
170
rival the most sterically-hindered arrangement of oxygen containing functional groups in SRFA,
171
and is potentially a good chemical proxy. Ketones such as in 1b are labile, reducing completely
172
to methylene groups. Aldehydes such as in 1c are also labile, reducing completely to methyl
173
groups. Ester linkages, such as in 1d, cleaved and reduced completely. This product was not
174
expected based on previous studies14
175
preservation of the newly created ether linkage was observed. Aliphatic alcohols are reduced,
176
but as is shown in 1e they can also be eliminated in particular cases, resulting in a double bond.
studies12
in which reduction of the carbonyl but
9
ACS Paragon Plus Environment
Environmental Science & Technology
Page 10 of 28
177
In contrast to previous studies, we have found that with a few exceptions double bonds are
178
labile, and this can be attributed to the higher catalyst concentrations employed here. In the
179
reduction of linear diene 1f no products are observed. This is believed to be the result of a
180
hydrosilylation reaction at the double bond(s).
181
reacted with reducing agent by-products, forming compounds which are lost during the post-
182
reduction workup. This behavior is observed in branched diene 1g and cyclic alkene 1h. The
183
hydrosilylation products 2g and 2h are in fact isolated but in low yields, which indicates that
184
material is lost to further reactions. Hydrosilylation of double bonds has been observed in
185
previous literature under similar conditions.17
186
easily identified by distinctive fragmentation patterns in mass spectra.
16
The resulting organosilane product further
The incorporation of n-BS into products is
187
In addition to hydrosilylation, double bonds can also be reduced under reaction conditions. In
188
the reduction of 1i, the isolated product 2i is fully reduced. However, the reduction product 2j is
189
only partially reduced. This is because double bonds present in a ring system are more resistant
190
to reduction, as is observed for products 2c and 2e. That double bonds are most stable between
191
tertiary carbons is clearly seen in the difference between the reduction products 2c and 2h, where
192
the double bond is preserved in 2c but hydrosilylated in 2h. The mechanism by which double
193
bonds are reduced vs. hydrosilylated is not yet fully understood. There does appear to be a
194
correlation of reduction with conjugated compounds (1i, 1j), however in a later section we
195
observe that non-conjugated terpenes are also reduced.
196
The reduction of 1e and 1j produces multiple isomers of products with identical molecular ions
197
and similar fragmentation patterns. Previous studies have shown that Lewis acid catalysts enable
198
backbone rearrangement (double bond migration and methyl shifts) in unsaturated sterol
199
hydrocarbons.1178 In those studies and our current work, rearrangements give rise to multiple
10
ACS Paragon Plus Environment
Page 11 of 28
Environmental Science & Technology
200
isomers from the reduction of a single compound, and the relative amounts of these products can
201
be altered by the amount of catalyst used. A combination of GC-MS, 13C NMR, and 1H NMR of
202
the reduction products was used to determine the location of the double bond in 2e. A detailed
203
description of this derivation is provided in Supporting Information. Thus far, isomer production
204
has only been observed in fused cyclic structures, which contained either double bonds or
205
hydroxyl groups directly bonded to the ring system prior to reduction. It is further proposed,
206
though not tested here, that other functional groups incorporated into the ring backbone, such as
207
ketones, would also behave in a similar manner to the hydroxyl reduction.
208
Finally, we have observed that the reduction of 1e produces both fully (2e’) and partially
209
reduced products (2e), with a molecular ion offset of 2 m/z. It appears that the partially reduced
210
products are not intermediates to the fully reduced product 2e’, as an increase in catalyst only
211
results in the loss of 2e (likely due to hydrosilylation) without an increase in 2e’. As such, it
212
appears that both are produced simultaneously by different pathwaysdifferent pathways generate
213
both reduction products simultaneously.
214
behavior. The production of multiple reduction products with different numbers of double bonds
215
from the same reactant ultimately prevents assignment of unsaturation state to a reactant of
216
unknown composition.
Of the compounds tested, only 1e exhibited this
217
In summary, it was found that all oxygen-containing groups, with the exception of phenols, are
218
reduced, that aromatic double bonds are stable but that aliphatic double bonds are only preserved
219
in ring systems, and that carbon-carbon bonds are not made or broken. The more complete
220
reduction compared to previous work12,13 is attributed to the elevated catalyst concentration that
221
is employed in the current method.
222
standards and samples.
Excess catalyst has been consistently used across all
11
ACS Paragon Plus Environment
Environmental Science & Technology
Page 12 of 28
223
Comprehensive GCxGC-TOF-MS. The distribution of hydrocarbon standards and reduced
224
model compounds in GCxGC space (Figure 1) mapped a chemical landscape that enables
225
preliminary assignment of broad functionality to unknown hydrocarbons based on 1st and 2nd
226
dimension retention times (1RT and 2RT respectively). Colored areas indicate where saturated
227
(alkane), unsaturated, and aromatic hydrocarbons are expected to elute. A thorough description
228
of the assignment of these regions is provided as Supporting Information.
229 230
Figure 1. A GCxGC selected-ion chromatogram (71 m/z) identifying the elution geography of
231
n-alkane standards. Circled regions overlay areas where depicted reduced model compounds and
232
hydrocarbons elute. Boxed compounds are isomers of aromatic C12H18 compounds. Note that
233
the reduction products of 1e and 1j elute as a series of peaks over a large region.
234
Terpenoids in SRFA: As previously mentioned, we were unable to separate the UCM into
235
discrete peaks using GC-FID/MS. However, GC-MS data still proved to be useful for structural
236
characterization. In a time averaged, baseline subtracted spectrum (Figure 2b) taken along the
237
UCM we observe clusters of prominent fragmentation ions which are spaced 14 m/z apart
238
(Figure 2c). Based on previous knowledge of EI-MS hydrocarbon fragmentation we suspected
239
that many closely related compounds contributed these fragments and used this knowledge to
240
further interrogate the sample by GCxGC-TOF-MS.
241
Three series of compounds are largely responsible for the diversity of fragment ions observed
242
in reduced SRFA by GCxGC-TOF-MS (Figures 2b, 2c). Each series contains many individual
12
ACS Paragon Plus Environment
Page 13 of 28
Environmental Science & Technology
243
compounds which occupy a two dimensional GCxGC space. These three series are displayed in
244
Figures 3a, 4a and 5a in 97 m/z, 95 m/z, and 107 m/z single ion chromatograms. Within each
245
series individual compounds yield very similar mass spectra, as certain fragment ions are highly
246
conserved across the series. For example, in the hydrocarbon series identified with 97 m/z
247
(Figure 3a) we observe 3 prominent fragmentation ions: 83 m/z, 97 m/z and 111 m/z (Figure 3c).
248
These ions are spaced at 14 m/z, the same spacing observed in GC-MS of the UCM, and are
249
present in almost all observed peaks from this series. Compounds identified in Figure 4a have
250
prominent fragment ions of 81 m/z, 95 m/z and 109 m/z (Figure 4c), and compounds in Figure
251
5a have ions of 79 m/z, 93 m/z and 107 m/z (Figure 5c). All together, these three series
252
represent the prominent fragment ions observed in GC-MS analysis. Additionally, the large
253
molecular size range (1RT range) and abundance distribution of each series are similar to the
254
features observed in the GC-FID chromatogram (Fig. 2a), and confirms the GCxGC observation
255
that the UCM is largely represented by the three series of compounds presented.
256
13
ACS Paragon Plus Environment
Environmental Science & Technology
Page 14 of 28
257 258
Figure 2-5. GC-FID chromatogram (2a) of reduced SRFA. GC-MS mass spectrum (2b) of
259
material eluting from 10.0-17.6 minutes.
260
fragmentation ions.
261
Selected- ion chromatograms of 97 m/z (3a), 95 m/z (4a) and 107 m/z (5a) display three distinct
262
series within a single GCxGC-TOF-MS run of reduced SRFA. Each series is comprised of
Inset in 2b shows (2c) clusters of prominent
Ions offset by 14 m/z are colored to indicate potential connectivity.
14
ACS Paragon Plus Environment
Page 15 of 28
Environmental Science & Technology
263
alicyclic hydrocarbons, with circled regions indicating isomers of the presented molecular
264
formula. Insets show possible structures for regions. Mass spectra (3b, 4b, 5b) of compounds
265
representative of the series are presented, and insets (3c, 4c, 5c) display prominent fragmentation
266
ions. * denotes peak location in chromatograms, molecular ion in spectra, prominent ions in 2c
267
It was further determined that the three series in reduced SRFA are comprised of alicyclic
268
hydrocarbons.
269
of the spectra from any of the series are consistent with hydrosilylated compounds. Model
270
compound work clearly demonstrated that all oxygen containing groups besides phenols are
271
reduced to hydrocarbons. Thus, chemical formulas for these compounds were easily calculated.
272
In the series depicted in Figure 3a, the peaks enclosed by a particular region all have the same
273
chemical formula. From these chemical formulae, it was determined that all three of the series
274
are comprised of hydrocarbons with varying degrees of unsaturation (n, where n = number of
275
rings plus double bonds). This assignment is also supported by elution of the three series in the
276
same unsaturated hydrocarbon region established in Figure 1. In model compound work it was
277
established that double bonds are only preserved on alicyclic compounds. Therefore, under these
278
conditions, reduction products that are unsaturated hydrocarbons must are bemost likely alicyclic
279
(with or without double bonds).
280
Many of the identified peaks have MS with detectable molecular ions, and none
These alicyclic reduction products are notable for a number of reasons. First, they are
281
broadly distributed over a wide molecular size range.
The most abundant compounds are
282
centered between C13 to C14, but some material appears to be larger than C 20 (Figure 4a).
283
Unfortunately, the analytical method did not preserve molecular ions for compounds >C17. The
284
size-range of compounds observed in the reduced sample has important implications for source
285
material structure, in that terpenes and terpenoids are the only biomolecules with alicyclic
15
ACS Paragon Plus Environment
Environmental Science & Technology
Page 16 of 28
286
backbones of this size.
Terpenoids are biomolecules synthesized from isoprene units that
287
contain five carbons. Thus terpenoid carbon backbones often contain 10, 15, 20, 25 and 30
288
carbons (although, larger terpenoids can be decarboxylated). By convention, terpenes are strictly
289
hydrocarbons, while terpenoids are oxygenated.918
290
Within a given series, alicyclic reduction products are offset by a CH2 group (14 m/z).
291
This methyl offset is not an artifact of the reduction because saturated C-C bonds are stable to
292
the reduction, and must be a feature of native, un-reduced SRFA. In fact, offsets of carboxylic
293
acid groups between individual compounds in Suwannee River Fulvic Acid have been observed
294
using FTICR-MS.7 When reduced, these offsets would appear as a difference of 14 m/z, or the
295
equivalent of a CH2 group in the molecular formula.
296
All compounds within a series are unsaturated to the same degree, and each series varies
297
by one degree of unsaturation, or 2 m/z. All compounds in Figure 3a have one degree of
298
unsaturation, meaning that they contain one alicyclic ring. Compounds displayed in Figure 4a
299
have two degrees of unsaturation. These compounds may contain a double bond on an alicyclic
300
structure, a fused cyclic structure, or two non-fused rings. Finally, compounds in Figure 5a have
301
three degrees of unsaturation. Possible structures that satisfy each chemical formula and the
302
requirements outlined above are depicted in insets within Figures 3a, 4a and 5a. Unfortunately,
303
exact structures are unavailable because the product MS are too ambiguous. The recovery of
304
compounds offset by 2 m/z (Figures 2b, 3b and 4b) would corroborate other studies which have
305
identified this feature in SRFA.7 However, as established in model compound work, the
306
unsaturation state of products is not necessarily reflective of the unsaturation state of reactants.
307
Thus, while the appearance of series offset by degrees of unsaturation is congruent with other
16
ACS Paragon Plus Environment
Page 17 of 28
Environmental Science & Technology
308
previous studies of SRFA, the potential that reduction artifacts contribute to this feature in the
309
current dataset must be considered a caveat.
310
There are a large number of peaks within an enclosed region that share the same
311
molecular ion. Mass spectra of these compounds contain the same fragment ions, but the
312
relative abundances of these ions vary. This confirms that each deconvoluted mass spectrum
313
belongs to a unique compound, but that all compounds within a region are structurally related.
314
Therefore, peaks within a region are very likely to be stereoisomers and/or positional isomers.
315
Based on this body of evidence, that compounds in reduced SRFA elute within the
316
chromatographic region occupied by unsaturated, alicyclic standard compounds, that double
317
bonds are only stable to theare unlikely to be stable in this reduction system reduction if unless
318
they are present within a ring system, that contiguous (no ether or ester bonds) carbon skeletons
319
between C10 and C17 dominate the chromatogram and are not created during the reduction, it can
320
be concluded that the majority of SRFA carbon that is detectable by GC following this reduction
321
(i.e., 12% of the SRFA carbon) resembles cyclic terpenes. This observation is consistent with
322
previous NMR-based studies that have suggested that terpenoids must be present in SRFA19
323
SRFA9 but these compounds have previously not been chromatographically isolated. To provide
324
some context to the reduction products observed in SRFA, a potential source material was
325
examined.
326
Terpenoids in bald cypress leaves. In the Suwannee River, a major DOM source is organic
327
matter leached from leaf litter.3
328
ran the chemical reduction on freshly fallen leaf litter gathered from Tifton, GA. The leaves are
329
of the genus Ttaxodinium, primarily from the bald-cypress (Ttaxodinium distichum). Bald-
330
cypress is present throughout the entire Suwannee River watershed, and is certain to be a major
19
To confirm the presence of terpenoids in source material we
17
ACS Paragon Plus Environment
Environmental Science & Technology
Page 18 of 28
331
contributor of DOM to SRFA. Although the potential of (leaf colonizing) bacterial and fungal
332
terpenoids must be considered, the intactness of our sample implies that any observed
333
compounds are ultimately plant derived.
334
separated using GCxGC and displayed using the 95 m/z chromatogram (Figure. 6).
The reduced leaf standard was chromatographically
335
The chemical reduction of these source terpenoids to hydrocarbons, and their distribution in
336
the unsaturated region of GCxGC space, confirms that compounds detected in reduced SRFA are
337
indeed derived from terpenes. For example, there are at least two isomers of pinene, a type of
338
bicyclic monoterpene, which are found in the leaves of conifers20 conifers20 and which are
339
present in the extracts of the leaf litter. Under the reaction conditions used, pinene would be
340
expected to reduce to a chemical formula of C10H18. We find that such products are present, and
341
elute near the classification region of reduced SRFA of the same chemical formula. Conifer
342
leaves are also known to contain bicyclic sesquiterpenes such as α-Cadinene2121 (found in leaf
343
extracts). The reduction product of such bicyclic terpenes would have the chemical formula
344
C15H28. Again, we find reduced hydrocarbons of that chemical formula elute in the same
345
GCxGC space as reduced SRFA products. This correspondence does not mean that intact
346
terpenes are present in SRFA. Rather, these smaller terpenoids function as additional model
347
compounds by further establishing the domain of reduced alicyclic compounds in GCxGC space.
348
Such a correspondence is especially reassuring for C15 alicyclic compounds, which are
349
prominent in reduced SRFA, but for which we have no purchased standards.
350
Consistent with previous studies of Taxodium wood,22 Tthe leaf mixture also contained
351
abundant reduction products of the diterpenoid resin acids and triterpenoid phytosterols. It was
352
previously mentioned that the reduction products of 1c are isomers, which elute over a
353
chromatographic region. This region is overlaid in Figure 6. There are at least four (and
18
ACS Paragon Plus Environment
Page 19 of 28
Environmental Science & Technology
354
potentially eight or greater) diterpenoid2s020, 2322 that are naturally produced by conifers. These
355
other compounds would also be expected to reduce to isomers, and this is confirmed by the
356
numerous isomers of diterpenoid reduction products that appear in the chromatogram of the
357
reduced leaf mixture (Figure 6).
358 359
360 361
Figure 6. 95 m/z selected- ion chromatogram of reduced bald-cypress leaf litter. Dark markings
362
indicate compound peaks. Regions corresponding to reduced terpenes/terpenoids are colored in
363
orange and labeled with size class. Also overlaid are series from Fig. 4a (black) and selected
364
reduced compounds and hydrocarbon standard from Fig. 1 (red).
365
In contrast to the resin acids, the reduction products of triterpenoid phytosterols appear as
366
discrete, but highly abundant, individual peaks. The small number of peaks indicates that
367
phytosterols are less diverse and/or not (as easily) isomerized by reaction conditions. Reduction
368
products of intact resin acids and phytosterols are not identified in the reduction products of
369
SRFA, which suggests that these compounds are rapidly removed during diagenesis or are
370
transformed into material that escapes detection. Reduced sesterterpenes, which are considered
371
to be rare in nature,243 are also present in the leaf litter sample. The reduced compounds elute in
372
the same GCxGC space as the reduction products of model compound 1e.
19
ACS Paragon Plus Environment
Environmental Science & Technology
Page 20 of 28
373
When considering the possible contribution of terpenoids as source material to SRFA, resin
374
acids are an attractive source, especially considering that the distribution of reduced isomers is
375
similar in shape to that of reduced SRFA. However, there are clearly SRFA reduction products
376
larger than C20 (Figure 4a), indicating that phytosterols are also likely to be a significant
377
contributor to SRFA. The cyclic portion of larger terpenoids such as carotenoids could also be a
378
source of the alicyclic compounds observed in SRFA. .
379
material confirms that cyclic terpenoids larger than C15 appear to be a major source for the
380
alicyclic compounds observed in reduced SRFA.
In short, analysis of potential source
381
Finally, in Scheme 1 we present a plausible pathway by which microbial degradation of parent
382
terpenoid compounds could generate the cyclic hydrocarbons observed in reduced SRFA. We
383
used the resin acid Abietic acid (1j) as our precursor model compound. Partially degraded
384
abietanes a and b have been observed as bacterial metabolites of 1j.254 These metabolites only
385
differ by one carbon number, which is the result of decarboxylation at the C-4 position. Further
386
degradation is hypothesized to occur via a c-ring cleavage, catalyzed by a known meta-cleavage
387
dioxygenase,2254 producing 2-isopropyl malic acid (c) in the process. However, no degradation
388
products beyond a and b have been reported on in the literature.
20
ACS Paragon Plus Environment
Page 21 of 28
Environmental Science & Technology
389 390
Scheme 1. Hypothetical degradation pathway from cyclic terpenoid to proposed metabolites
391
(components of native SRFA), which then could be further degraded to produce the compounds
392
observed in reduced SRFA.
393
We suggest that further metabolites (d, e) could exist that are either not amenable to recovery,
394
or are extremely short lived. Such bicyclic metabolites could degrade further via a biological
395
Baeyer-Villiger oxidation, forming f and g. The enzyme which catalyzes this reaction is non-
396
specific and capable of oxidizing a large range of cyclic substrates.265 These could then be
397
environmentally hydrolyzed, forming h and i, and/or potentially ester linked into larger
398
compounds.
Reduction of these hypothetical metabolites (d, e, h, i) could produce carbon
21
ACS Paragon Plus Environment
Environmental Science & Technology
Page 22 of 28
399
backbones of the size and unsaturation observed in reduced SRFA. Although speculative, it
400
appears that degradation via decarboxylation (and then reduction) is the simplest explanation for
401
the observed CH2 offset of products in reduced SRFA.
402
Molecular Diversity and Future DirectionsMethod contributions to DOM analysis.
403
In a greatly simplified framework, the extreme molecular diversity of NOM can be divided into
404
compositional and structural diversity.We have effectively shown that chemical reduction
405
followed by GCxGC analysis of SRFA allows for moderate interrogation of both structural and
406
compositional diversity, where C compositional diversity refers to the diversity of molecular
407
formulae within a mixture of compounds, and s. Structural diversity refers to the arrangement of
408
elements within a single compound. Both are indicators of source and diagenetic history. We
409
have effectively shown that chemical reduction followed by GCxGC analysis of SRFA allows
410
for moderate interrogation of both structural and compositional diversity.
411
Although capable of revealing compositional diversity in the reaction products, compositional
412
diversity in native SRFA is not directly supported. As the reduction method is comprehensive, it
413
is not possible to determine how many oxygenated functional groups were reduced to form the
414
observed compounds, and therefore it is impossible to extrapolate the exact original composition
415
of reduced compounds. However, we can be certain that direct carbon-carbon bonds were not
416
broken or created during the reduction.
417
compositional diversity, in the context of carbon skeletons, in unreduced SRFA. The distribution
418
of reduced compounds has a maximum in intensity at a molecular formula of C14H26. Unreduced
419
SRFA has a peak maximum at around 350 m/z, which would correspond to a carbon number per
420
molecule at ~15-20.0.2626 It is therefore unlikely that individual SRFA compounds would contain
This allows us to make a reasonable argument for
22
ACS Paragon Plus Environment
Formatted: Indent: First line: 0"
Page 23 of 28
Environmental Science & Technology
421
more than one degraded terpenoid backbone. If true, then the observed compositional diversity
422
in reduced SRFA is in part a result of the compositional diversity of unreduced SRFA.
423
There is a great deal of structural diversity present within the UCM, supported by the many
424
isomers identified for each region of a given chemical formula. As a caveat, an unknown
425
amount of diversity could be an artifact of catalyst induced isomerization. However this effect
426
has only been observed in fused cyclic model compounds. The series displayed in Figure 2a are
427
alicyclic, but contain only one ring and no double bonds. Thus, the many compounds observed
428
within each region are certainly not artifacts of the reduction. This demonstrates structural
429
variation within compositionally identical carbon backbones of reduced SRFA, and therefore
430
native SRFA. Broad signals in NMR spectra have been attributed to structural diversity, but to
431
the best of our knowledge this is the first time that structurally distinct isomers have been
432
recovered from SRFA.
433
The isolation of a significant fraction of terpenoid-derived material indicates that all classes of
434
biomolecules are not preserved equally within humic substances. This observation is striking,
435
considering that terpenoids are likely to account for only a small percentage of total organic
436
matter input.
437
preferentially preserved within SRFA.
438
sediments28 is well known, but the current findings support long-term accumulation of these
439
biomarkers in the dissolved phase. Further studies are necessary to confirm if this trend is also
440
observed across different ecological settings, including marine environments.
The results presented here indicate that terpenoid derived compounds are The preferential accumulation of terpenoids in
441 442
ASSOCIATED CONTENT
23
ACS Paragon Plus Environment
Environmental Science & Technology
Page 24 of 28
443
Supporting Information. Previously published mechanism of chemical reduction; rationale for
444
GCxGC geography; structure rational for produced isomers. This material is available free of
445
charge via the internet at http://pubs.acs.org.
446
AUTHOR INFORMATION
447
Corresponding Author
448
* email:
[email protected].
449
ACKNOWLEDGMENTS
450
The authors thank the National Science Foundation for funding (NSF OCE-1155269 to LIA), as
451
well as the following researchers for inspiration, chemical insights, chromatographic support
452
and foliage: Cameron McIntyre, Vladimir Gevorgyan, Rama Nimmagadda, Brad Allpike, Jack
453
Cochran, Christoph Aeppli, Bob Nelson, Andrew Mehring and Casey Harris.
454 455 456 457 458
REFERENCES 1. In Humic Substances in Soil, Sediment and Water; Aiken, G. R.; Mcknight, D. M.; Wershaw, R. L. MacCarthy, P., Eds.; Wiley: New York, 1985 2. Benner, R.
Chemical Composition and Reactivity.
In Biogeochemistry of Marine
459
Dissolved Organic Matter; Hansell, D. A., Carlson, C. A., Eds.; Academic Press: San
460
Diego, 2002; pp 59
461
3. Leenheer, J. A.; Nanny, M. A.; McIntyre, C. Terpenoids as major precursors of dissolved
462
organic matter in landfill leachates, surface water, and groundwater. Environ. Sci.
463
Technol. 2003, 37(11), 2323-2331
24
ACS Paragon Plus Environment
Page 25 of 28
Environmental Science & Technology
464
4. Lam, B.; Baer, A.; Alaee, M.; Lefebvre, B.; Moser, A.; Williams, A.; Simpson, A.J.
465
Major structural components in freshwater dissolved organic matter. Environ. Sci.
466
Technol. 2007, 41(24), 8240-8247.
467
5. Woods, G. C.; Simpson, M. J.; Simpson, A. J. Oxidized sterols as a significant
468
component of dissolved organic matter: Evidence from 2D HPLC in combination with
469
2D and 3D NMR spectroscopy. Water Research. 2012, 46, 3398-3408.
470
6. Kujawinski, E. B.; Hatcher, P. G.; Freitas, M. A. High-resolution Fourier transform ion
471
cyclotron resonance mass spectrometry of humic and fulvic acids: Improvements and
472
comparisons. Anal. Chem. 2002, 74, 413-419
473
7. Stenson, A. C.; Marshall, A. G.; Cooper, W. T. Exact masses and chemical formulas of
474
individual Suwannnee River fulvic acids from ultrahigh resolution electrospray ionization
475
Fourier transform ion cyclortron resonance mass spectra. Anal. Chem. 2003, 75, 1275-
476
1284
477 478
7.8.Ourisson, G. The General Role of Terpenes and their global significance. Pure & Applied Chemistry, 1990, 62, 1401-1404
479
8.9.Leenheer J.A.; Rostad, C.E. Tannins and terpendoids as major precursors of Suwannee
480
River fulvic acid. U.S. Geological Survey investigations Report 2004-5276, Denver, CO;
481
U.S. Geological Survey, 2004.
482 483
9.10.
Ertel, J. R.; Hedges, J. I., Perdue, E. M., Lignin Signature of Aquatic Humic
Substances. Science, 1984, 223(4635), 485-487.
25
ACS Paragon Plus Environment
Environmental Science & Technology
484
10.11.
Page 26 of 28
Woods, G. C; Simpson, M. J.; Koerner, P. J.; Napoli, A.; Simpson, A. J. HILIC-
485
NMR: Toward the Identification of individual molecular components in dissolved
486
organic matter. Environ. Sci. Technol. 2011, 45(9), 3880-3886
487
11.12.
Nimmagadda, R. D.; McRae, C. R.; Characterisation of the backbone structures of
488
several fulvic acids using a novel selective chemical reduction method. Org. Geochem.,
489
2007, 38, 1061-1072.
490
12.13.
491 492
Advanced characterization of natural organic matter (NOM) in Australian water
supplies; Research report 80; Water quality research Australia: Adelaide, SA, 2010 13.14.
Frysinger, G. S.; Caines, R. B.; Xu, L.; Reddy, C. M. Resolving the unresolved
493
complex mixture in petroleum-contaminated sediments. Environ. Sci. Technol. 2003, 37,
494
1653-1662
495
15. Ball, G.I., Aluwihare, L.I. CuO-oxidized dissolved organic matter (DOM) elucidated by
496
comprehensive two dimensional gas chromatography-time of flight-mass spectrometry
497
(GCxGC-TOF-MS) Organic Geochemistry. 2014, 87-98
498
14.16.
499
Rubin, M.; Schwier, T.; Gevorgyan, V. Highly efficient B(C6F5)3-Catalyzed
Hydrosilylation of Olefins. J. Org. Chem. 2002, 67, 1936-1940
500
17. Peakman, T. M.; Maxwell, J.R. Acid-catalysed rearrangements of steroid alkenes. Part 1.
501
Rearrangement of 5α-Cholest-7-ene. J. Chem. Soc. Perkin. Trans. 1. 1988, 1065-1070
502
15.
503
16.18.
504
Zwenger, S.; Basu, C. Plant terpenoids: applications and future potentials.
Biotechnology and molecular biology reviews. 2008, 1-7
26
ACS Paragon Plus Environment
Page 27 of 28
505
Environmental Science & Technology
17.19.
Averett, R. C.; Leenheer, J. A.; McKnight, D. M.; Thorn, K. A. Humic Substances
506
in the Suwanneee River, Georgia: Interactions, Properties, and Proposed Structures:
507
U.S. Geological Survey water-Supply Paper 2373; U.S. Government Printing Office:
508
Washington, DC, 1994
509 510 511 512
18.20.
Gershenzon, G. J.; Dudareva, N. The function of terpene natural products in the
natural world. Nat. Chem. Bio.2007, 3(7), 408-414 19.21.
Yatagai, M.; Sato, Y.; Takahashi, T. Terpenes of leaf oils from Cupressaceae.
Biochm. System. A. Ecol. 1985, 13, 377-385
513
22. Zhong, A.; Sleighter, R.L.; Salmon, E.; McKee, G.E.; Hatcher, P.G. Combining advanced
514
NMR techniques with ultrahigh resolution mass spectrometry: A new strategy for
515
molecular scale characterization of macromolecular components of soil and sedimentary
516
organic matter. Organic Geochemistry. 2011, 903-916
517
20.23.
Himejima, M.; Hobson, K. R.; Otsuka, T.; Wood, D. L.; Kubo, I. Antimicrobial
518
terpenes from oleoresin of ponderosa pine tree Pinus ponderosa: A defense mechanism
519
against microbial invasion. Journal of Chem. Ecol. 1992, 18, 1809-1818
520 521 522
21.24.
Cordell, G. A. The occurrence, structure and elucidation and biosynthesis of the
sesterterpenes. Phytochemistry.1974, 13, 2343-2364 22.25.
Martin, V. J. J.; Yu, Z.; Mohn, W. W.; Recent advances in understanding resin
523
acid biodegradation: microbial diversity and metabolism. Arch Microbiol. 1999, 172,
524
131-138
27
ACS Paragon Plus Environment
Environmental Science & Technology
525
Page 28 of 28
26. Mihovilovic, M. D.; Muller, B.; Stanetty, P. Monooxygenase-mediated Baeyer-Villiger
526
Oxidations. Eur. J. Org. Chem. 2002, 3711-3730
527
23.
528
27. Witt, M.; Fuchser, J.; Koch, B.P. Fragmentation studies of fulvic acids using collision
529
induced dissociation fourier transform ion cyclotron resonance mass spectrometr. Anal.
530
Chem.2009, 2688-2694
531 532
24. Brassell, S.C.; Eglinton, G.; Maxwell, J.R. The geochemistry of terpenoids and steroids.
533 534
Biochem Soc Trans. 1983, 5,575-86 25. Stuermer, D.H.; Harvey, G.R. Structural studies on marine humus: a new reduction
535
sequence for carbon skeleton determination. Marine Chemistry 6.1 1978, 55-70
536
26. Leenheer, J.A.; Wershaw, R.L.; Reddy, M.M. Strong-acid, carboxyl-group structures in
537
fulvic acid from the Suwannee River, Georgia. 2. Major structures. Environ. Sci. Technol. 1995,
538
399-405
539 540 541
Formatted: Indent: First line: 0.13", No bullets or numbering
27. Advanced characterization of natural organic matter (NOM) in Australian water supplies; Research report 80; Water quality research Australia: Adelaide, SA, 2010 28.
Formatted: TF_References_Section, Space After: 0 pt, Numbered + Level: 1 + Numbering Style: 1, 2, 3, … + Start at: 1 + Alignment: Left + Aligned at: 0.25" + Indent at: 0.5"
28
ACS Paragon Plus Environment