Direct Identification of Diverse Alicyclic Terpenoids in Suwannee

To advance compositional studies of DOM we have applied a previously published reduction method to an environmental standard, Suwannee River Fulvic Ac...
0 downloads 0 Views 647KB Size
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