FTIR and Synchronous Fluorescence Heterospectral Two

Jan 22, 2015 - Dimensional Correlation Analyses on the Binding Characteristics of. Copper onto Dissolved Organic Matter. Wei Chen, Nuzahat Habibul, ...
0 downloads 0 Views 1MB Size
Subscriber access provided by SELCUK UNIV

Article

FTIR and Synchronous Fluorescence Heterospectral Two-Dimensional Correlation Analyses on the Binding Characteristics of Copper onto Dissolved Organic Matter Wei Chen, Nuzahat Habibul, Xiao-Yang Liu, Guo-Ping Sheng, and Han-Qing Yu Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/es5049495 • Publication Date (Web): 22 Jan 2015 Downloaded from http://pubs.acs.org on February 5, 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

FTIR and Synchronous Fluorescence Heterospectral Two-Dimensional Correlation Analyses on the Binding Characteristics of Copper onto Dissolved Organic Matter

Wei Chen, Nuzahat Habibul, Xiao-Yang Liu, Guo-Ping Sheng, Han-Qing Yu* CAS Key Laboratory of Urban Pollutant Conversion, Department of Chemistry, University of Science and Technology of China, Hefei, 230026, China

1 ACS Paragon Plus Environment

Environmental Science & Technology

1

Dissolved organic matter (DOM) is known to form strong complexes with heavy

2

metals and thus governs the distribution, toxicity, bioavailability, and the ultimate fate

3

of heavy metals in the environment. The relevant aspects of metal-organic interactions

4

remain unclear because the metal binding functionalities in DOM are substantially

5

non-uniform and the availability of the models is limited. In this work, two-

6

dimensional correlation spectroscopy (2DCOS) integrated with synchronous

7

fluorescence and infrared absorption spectroscopy was used to explore the binding

8

process of copper to DOM. A series of heterogeneous binding sites in humic acid

9

(HA), a representative DOM, and the subsequent subtle changes of these sites within

10

the molecular interactions were elucidated by 2DCOS method. The band assignments

11

and the correspondence between the results obtained by two spectral probes

12

(synchronous fluorescence and infrared absorption spectra) were verified by hetero-

13

2DCOS. Our results showed that, during the copper binding process, the carboxyl and

14

polysaccharide groups gave the fastest responses to copper binding. Then,

15

fluorescence quencing of fluorescent humic-like moieties occurred with vibrational

16

change of the related functionalities, i.e., phenolic and aryl carboxylic groups, which

17

further induces the fluorescence quenching of fulvic-like fractions. Finally, small

18

amounts of amide and aliphatic groups participated in the copper binding after the

19

fluorescence of the protein-like fraction decreased. With these promising results, a

20

comprehensive picture of structural changes of HA during copper binding process was

21

developed, highlighting the superior potential of 2D-heterospectral correlation

22

spectroscopy in studying complex interactions in environment.

2 ACS Paragon Plus Environment

Page 2 of 28

Page 3 of 28

Environmental Science & Technology

23

INTRODUCTION

24 25

Ubiquitously existing in aquatic and soil environments with various chemical

26

functional groups,1-3 dissolved organic matter (DOM) is known to form strong

27

complexes with heavy metals, thus determining the distribution, toxicity,

28

bioavailability, and the ultimate fate of heavy metals in the environment.4-6

29

Carboxylic and phenolic groups in DOM are considered as the predominant

30

coordination sites with metal ions.7, 8 Some recent reports, however, have suggested

31

the presence of other types of metal coordination arising from the heterogeneous

32

distribution of the metal binding sites in DOM.9, 10

33

Approaches studying the metal-binding characteristics of DOM are quite diverse,

34

including size exclusion chromatography, equilibrium dialysis, differential absorbance

35

titration, and fluorescence excitation-emission matrix coupled with PARAFAC

36

analysis.6, 11-13 Optical spectroscopy is the most frequently used tool for metal-DOM

37

interaction analysis owing to its characteristic responses to the abundant light

38

absorbing and emitting groups in metal-DOM complexes. Models like NICA-

39

Donnan14 or Stockholm Humic Model,15 and chemometric methods like PARAFAC

40

analysis,10 have been employed to enhance the spectral resolution to reveal the subtle

41

spectral changes and highly overlapping responses during the binding process.

42

However, since the metal binding functionalities in DOM are substantially non-

43

uniform, the availability of models is limited. Therefore, some aspects of metal-

44

organic interactions, such as the copper binding affinities of functional groups in HA,

45

are unclear yet.

46

Two-dimensional correlation spectra (2DCOS) is capable of resolving overlapped

47

peaks by extending spectra along the second dimension as well as providing

3 ACS Paragon Plus Environment

Environmental Science & Technology

48

information about the relative directions and sequential orders of structural variations.

49

Thus, it has been recently applied to probe the interaction mechanisms of some

50

environmental-related substances.9, 16-18 The major probes used in 2DCOS are UV-Vis,

51

fluorescence and IR, each of which gives particular molecular information about the

52

binding process. 2D hetero-spectral correlation analysis integrates two different types

53

of spectra for a system obtained by using multiple spectroscopic probes under the

54

same external perturbation.19, 20 Correlation between the two complementary spectral

55

signals can help to better understand the structural variation caused by the

56

perturbation. For these appealing features, hetero-2DCOS is among the most active

57

areas of 2DCOS research.21-23

58

In the present work, by using a novel hetero-2DCOS analysis which integrates

59

synchronous fluorescence and Fourier transform infrared (FTIR) absorption

60

spectroscopy, we aim to provide an in-depth understanding about the heterogeneous

61

binding characteristics of copper onto DOM at the molecular level. This might, to the

62

best of our knowledge, be the first time to use such a new hetero-2DCOS method to

63

study the metal-DOM interactions. Humic acid (HA) was used as the representative

64

DOM and copper concentration was used as the external perturbation. The 2D

65

fluorescence, 2D IR and fluorescence/IR hetero-2DCOS were compared and

66

combined to elucidate the copper-DOM interaction mechanisms.

67 68

EXPERIMENTAL SECTION

69 70

Sample Preparation. Commercial HAs (Sigma-Aldrich Co., USA) were purified by

71

dissolving in NaOH solution (pH = 13.0) with subsequent filtration and acidifying the

72

filtrate with HCl to pH 1.0. The precipitate was collected by filtration, washed

4 ACS Paragon Plus Environment

Page 4 of 28

Page 5 of 28

Environmental Science & Technology

73

extensively with 0.1 M HCl and freeze-dried. AR-grade CuCl2, NaOH and HCl were

74

purchased from Sinopharm Chemical Reagent Co., China, and were used as received.

75

The purified HAs were dissolved and diluted to a concentration of 12.5 mg/L (pH

76

7.0). Cu2+ titrations were carried out by adding appropriate volumes of 1 mM CuCl2

77

stock solution into 20 mL HA solution to generate a series of samples containing Cu2+

78

concentration in a range of 0-160 µM. The final volume was 25 mL and thus the HA

79

concentration was fixed at 10 mg/L. Background ionic strength was maintained by

80

adding 0.04 M KBr (IR grade). The pH was maintained constant by adjusting both

81

HA and Cu2+ stock solution to pH 7.0, and the mixed solution was fine tuned with 0.1

82

M NaOH or HCl solution. Each solution was shaken for 12 h at 25 oC to ensure

83

coordination equilibrium. Then, 5 mL of each solution was analyzed by synchronous

84

fluorescence spectroscopy and the rest was freeze-dried for further FTIR spectroscopy

85

analysis.

86

Procedures and Parameters of Measurements. The elemental compositions of

87

HA samples were analyzed using an elemental analyzer (vario EL cube, Elementar

88

Co., Germany). Synchronous fluorescence spectra of HA solutions with increasing

89

copper concentration were obtained by an average of 3 scans using a luminescence

90

spectrometer (Perkin–Elmer LS-55, USA). Excitation and emission slits were both

91

adjusted to 10 nm, and the excitation wavelengths ranging from 250 to 550 nm were

92

used in 0.5 nm increment at a scan rate of 1000 nm/min with a constant offset ∆λ of

93

60 nm. 60 nm was chosen as the offset to provide a higher fluorescence intensity and

94

a better resolution compared with other offsets (Figure S1).9 50 mg of the freeze-dried

95

HA samples with different copper concentration were ground, homogenized and

96

pressed under the irradiation of an infrared lamp to eliminate the influence of moisture

97

in samples. Their transmission IR spectra were recorded on a Vertex 70 spectrometer

5 ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 28

98

(Bruker Co., Germany) with a DTGS detector, and each spectrum was obtained after

99

32 scans with 4 cm-1 resolution. The IR spectra were smoothed, baseline-corrected

100

and finally transformed to absorbance spectra using OPUS 5.5 software for the

101

subsequent analysis.

102

2DCOS Analysis. To obtain the structural variation information of copper-bound

103

HA, 2DCOS was employed using synchronous fluorescence and FTIR spectra with

104

copper concentration as the external perturbation. The 2DCOS were produced

105

according to the method of Noda and Ozaki.24 The analytical spectral changes of y(x,

106

C), as a function of a spectral variable (x, here representing λ or v) and an external

107

variable (C, copper concentration) at m evenly spaced points, can be represented as:

108 109

110

  = ,  

(1)

A set of dynamic spectra  ,  is defined as follows:

   = 

   −   for 1 ≤  ≤  0 otherwise

(2)

111

where   = ! ∑! #   , denoting the reference spectrum, typically the m variable

112

averaged spectrum. A computational method based on the discrete Hilbert-Noda

113

transform is used to generate a pair of correlation intensity maps, i.e., synchronous (Ф)

114

and asynchronous (Ψ) correlation spectra: 0 if  = *

115

$̃   = ∑! '  , N' = ) '# N' ∙

116

Ф , .  = !, ∑!    ∙   .  #

(4)

117

Ψ , .  = !, ∑!    ∙ $̃  .  #

(5)

+',

otherwise

(3)

118

Spectral coordinates, intensities and signs of correlation peaks appearing on 2D

119

spectra can be interpreted by a set of well-established principles.24, 25 Synchronous

120

spectra, corresponding to the real part of the cross correlation function, consist of

6 ACS Paragon Plus Environment

Page 7 of 28

Environmental Science & Technology

121

auto-peaks located along the main diagonal of the map and cross-peaks located at the

122

off-diagonal positions. Higher susceptibility of intensity changes to the perturbation

123

generates greater positively signed auto-peak intensities in the synchronous spectra.

124

Cross-peaks suggest the coordinated changes of spectral intensities observed at two

125

different spectral variables (x1, x2). A positive cross peak indicates the same direction

126

of the intensities change at the corresponding spectral coordinates, while a negative

127

value suggests an opposite direction. Asynchronous spectra, the imaginary part of the

128

cross correlation function, show cross-peaks exclusively. Their signs reveal the

129

sequential order of the dynamics of spectral intensity variations induced by the

130

perturbation. Same signs of the spectral coordinate (x1, x2) in both the synchronous

131

and asynchronous maps indicate that the spectral intensity change at x1 occurs

132

predominantly before that of x2 along the perturbation variable axis. This order is

133

reversed when signs are opposite. The spectral changes occur simultaneously if only

134

the synchronous correlation intensity is observed, and the temporal relationship of

135

spectral intensity changes becomes indeterminate if only the asynchronous correlation

136

intensity is observed.

137

The spectra of each sample for 2DCOS analysis have subtracted the background

138

signal (without addition of copper) first. The data set was then transformed into a new

139

spectral matrix suitable for 2DCOS maps using 2D Shige software (Kwansei-Gakuin

140

University, Japan). Synchronous fluorescence 2DCOS, IR 2DCOS and synchronous

141

fluorescence/IR hetero-2DCOS maps were plotted using Origin 8.5 software.

142

Determination of Copper Binding Parameters. The fluorescence quenching

143

titration method was employed to quantify the copper binding characteristics of

144

fluorescent fractions of HA, and the parameters related to copper binding were

145

estimated by the following modified Stern-Volmer equation:26

7 ACS Paragon Plus Environment

Environmental Science & Technology

146

01

01 ,0

=

2∙34 ∙54

+

Page 8 of 28

(6)

2

147

where F, F0, f, and CM represent the measured fluorescence intensity, the initial

148

fluorescence intensity, the fraction of the initial fluorescence that corresponds to the

149

binding fluorophores, and the Cu2+ concentration, respectively. The conditional

150

stability constant KM can be obtained by plotting F0/ (F0-F) against 1/CM. The

151

fluorescence wavelengths were selected with a consideration of the pronounced peaks

152

from the 2DCOS maps.

153

154

RESULTS AND DISCUSSION

155 156

General Spectroscopic Observations. DOM is a complex mixture of organic

157

compounds formed by decomposition of plants and animal residues by microbial

158

activity, which contains carboxyl, phenol, quinonyl, ester, ketone, hydroxyl, amino

159

and other functional groups. Elemental analysis and structural characterization of the

160

purified HA are given in Table S1 and Figure S2. Figure 1 shows the change of

161

synchronous fluorescence of HA with the addition of Cu2+. Three fluorescence

162

regions, corresponding to the wavelength ranges of 250-300 nm, 300-380 nm and

163

380-550 nm, could be roughly assigned to protein-like, fulvic-like and humic-like

164

fluorescence fractions, respectively.27 The fluorescence intensity of protein-like

165

fractions was relatively low, compared with that of the sharp fulvic-like peaks and the

166

broad humic-like shoulders. An increase in copper concentration caused a higher

167

extent of fluorescence quenching for all the fractions, indicating the occurrence of

168

electronic structural changes in the fractions by forming complexes with copper.

8 ACS Paragon Plus Environment

Page 9 of 28

Environmental Science & Technology

169

The predominant infrared spectral characteristics of HA occurred in the range of

170

1750-750 cm-1, where almost all the vibrational information of HA backbones could

171

be identified. The spectral shift of HA upon addition of copper was shown in Figure

172

2(a), in which the characteristic bands changed to various degrees with variation of

173

copper concentration. The effect of copper to the structural variation of HA was

174

illustrated more clearly in differential IR absorbance spectra in Figure 2(b), in which

175

each of the spectra was subtracted by the spectra without copper addition and the

176

trend for the spectral change was amplified. Even so, some of the IR bands remained

177

strongly overlapped with uncertain assignments. It can be seen that the intensities of

178

these IR bands diminished with increasing copper concentration, except for the bands

179

centered at 800 cm-1. These observations indicated that during the binding process

180

with copper, the vibrational structural changes varied in the fractions of HA

181

backbones.

182

The one-dimensional synchronous fluorescence and FTIR spectra failed to give

183

detailed information on the exact binding characteristics of HA to copper, owing to

184

the strong overlapped absorbances caused by the diverse groups of HA. To achieve a

185

reliable determination of the conformational changes of HA induced by copper

186

binding, enhancement of the spectral resolution and understanding the correlation

187

between the fluorescence and IR responses are needed.

188

Correlation Analysis. 2DCOS helps to enhance the spectral resolution by

189

spreading the spectrum in a second dimension, thus facilitating the deconvolution of

190

overlapping peaks. Spectral features of different molecular origins that make similar

191

contributions to the spectra in the static case may behave differently (either delayed or

192

accelerated) in the case of a dynamic change. Such differences can be distinguished

193

by 2DCOS. The synchronous 2DCOS map for the synchronous fluorescence spectra

9 ACS Paragon Plus Environment

Environmental Science & Technology

194

of HA in the region 280-500 nm (Figure 3a) showed three predominant auto-peaks

195

centered at 445, 380 and 350 nm with a small peak at 292 nm identified from the

196

cross-peaks. Compared with one dimensional fluorescence spectra (Figure 1), the

197

2DCOS gave peaks responsible for fluorescent fractions of HA with a higher

198

resolution. Intensities of peaks decreased in the order of 445, 380, 350 and 292 nm,

199

suggesting that the fluorescence of humic-like fractions were more susceptible to

200

copper addition, while the fluorescence of protein-like fractions were less sensitive,

201

and the fulvic-like fractions could be divided into two components (380 and 350 nm).

202

All the cross-peaks were positive, indicating that the spectral changes proceeded in

203

the same direction as the variation of copper ion concentration.

204

The asynchronous map provided additional useful information about the

205

sequential relationship between two spectral origins during copper binding. The

206

fluorescence intensity changes at one wavelength raises or lags those at other

207

wavelengths. As shown in Figure 3b and Table S2, the cross-peaks located at the

208

upper-left corner of the asynchronous map all exhibited negative signs, which

209

suggested that, according to Noda’s rule, the change followed the order of 445 > 380,

210

350 > 292 nm. This suggested that copper bound to HA fractions in the following

211

sequence: humic-like fraction → fulvic-like fraction → protein-like fraction.

212

The various functional groups of HA always lead to overlapping IR peaks,

213

resulting in ambiguous band assignments confusing the conformational analysis of

214

HA. 2DCOS is a promising tool to solve this problem. Figure 4 showed the copper-

215

induced 2D-IR-COS of HA. The synchronous map gave seven predominant auto-

216

peaks centered at 1685, 1620, 1360, 1100, 1020, 850, 780 cm-1 and a small peak at

217

1260 cm-1 (Figure 4a). The band of 1380 cm-1 was shown to change most significantly.

218

In contrast, the smallest change was observed for bands at 1100 and 850 cm-1. It could

10 ACS Paragon Plus Environment

Page 10 of 28

Page 11 of 28

Environmental Science & Technology

219

be seen that the the signs of the cross-peaks in the synchronous map were all positive,

220

except for ones associated with 850 and 780 cm-1. This indicated that conformational

221

changes associated with most IR bands underwent the same direction, while those

222

associated with 850 and 780 cm-1 bands proceeded reversely.

223

The peak at 1685 cm-1 splitted into two parts at 1720 and 1660 cm-1 respectively

224

in the asynchronous map (Figure 4b), and the peak at 1360 cm-1 was overlapped with

225

peaks at 1400 and 1350 cm-1. The 1720 and 1260 cm-1 singals can be attributed to the

226

C=O streching and C-O streching (or OH deformation) of carboxylic acid groups,

227

while the bands at 1350 and 1020 cm-1 corresponded to the aliphatic C-H deformation

228

and C-OH streching, and the bands at 850 and 780 cm-1 corresponded to the aromatic

229

C-H deformation.2 The detailed assignments of bands and signs of their cross-peaks in

230

the asynchronous map were shown in Table 1. Using the sequential order rules, it

231

could be concluded that the structural change sequence of HA backbones by copper

232

binding followed the order: carboxyl → C-O of polysaccharides → phenolic groups

233

→ aryl → amide → aliphatic groups. The 2D-IR-COS results suggested that not only

234

the fluorescent groups, but also some non-fluorescent groups (like polysaccharides) of

235

HA were involved in the binding process of copper with HA, and the binding

236

affinities of hydrophilic sites were higher than that of hydrophobic sites.

237

Fluorescence/IR Hetero-2DCOS. To verify the band assignments as well as the

238

correspondence between the two complementary tehcniques, fluorescence and IR

239

spectra that represent different parts of HA structure, we employed a 2D

240

heterospectral fluorescence/IR correlation analysis. The hetero-2D spectra could serve

241

as a visual aid to help describe the structural change of HA during the copper binding.

242

Figure 5 shows the synchronous and asynchronous maps of the hetero-2DCOS,

243

where the IR frequency and the fluorescence shift were ploted on the x-axis and y-axis

11 ACS Paragon Plus Environment

Environmental Science & Technology

244

respectively. Remarkably, four narrow cross-peaks, located in the IR regions of 1650,

245

1350, 850 and 780 cm-1 and the corresponding fluorescence region of 450-350 nm,

246

appeared in the synchronous map (Figure 5a). The former two peaks showed positive

247

signs while the latter two were negative. These results indicated that the carboxyl,

248

phenolic and aryl groups of HA were the basic fluorescent units for humic-like and

249

fulvic-like fractions. In the asynchronous map (Figure 5b), the cross-peaks were

250

divided into different parts, and the humic-like fraction (445 nm) was related to the

251

phenolic (1380 cm-1) and aryl (1620, 850, 780 cm-1) groups, while the fulvic-like

252

fraction (380-350 nm) was related to the carboxyl (1720, 1260 cm-1), phenolic (1400

253

cm-1) and aryl groups (1620, 850 cm-1). Two new cross-peaks were also observed in

254

the asynchronous map at 1660 and 1350 cm-1 with the fluorescence peak at 290 nm,

255

corresponding to the fluorescence region of protein-like fraction. This suggested that

256

the protein-like fraction was related to the amide and aliphatic groups.

257

As a result of the difference in IR and fluorescence spectra probing the tinily

258

different aspects of molecular responses toward an identical perturbation (copper

259

concentration), asynchronicity was observed between two different types of

260

spectroscopic data of the same species. In the asynchronous map, there were two

261

negative signs at coordinates (1620 cm-1, 445 nm) and (1380 cm-1, 445 nm), and two

262

positive signs at (850 cm-1, 445 nm) and (780 cm-1, 445 nm), indicating that the

263

fluorescence response of humic-like fraction occurred before the spectral changes of

264

phenolic, aryl, and Ar-H groups detected by IR. On the other hand, the positive signs

265

of fluorophore coordinates at 380-350 nm and 290 nm suggested that the fluorescence

266

change of fulvic-like fraction occurred after the vibrational changes of carboxyl,

267

phenolic and aryl groups but before that of Ar-H groups, and the fluorescence change

268

of protein-like fractions occurred after the vibrational changes of amide and aliphatic

12 ACS Paragon Plus Environment

Page 12 of 28

Page 13 of 28

Environmental Science & Technology

269

groups. This order of intensity changes are in good agreement with that obtained by

270

2D-fluorescence and 2D-IR correlation studies.

271

To quantitatively confirm our interpretation of the 2DCOS results, copper binding

272

parameters were calculated by modified Stern–Volmer equation, which has been

273

widely used for characterizing the metal binding behavior of DOM, at the

274

wavelengths corresponding to the main peaks of the 2DCOS maps (Figure 6). The

275

calculated values of the conditional stability constants (Table 2) were comparable to

276

those reported in similar studies. And they decreased in the order 445 > 380 > 350

277

> 292 nm, supporting our aforementioned interpretation of the fluorescence 2DCOS

278

maps. In addition, the fraction of the initial fluorescence also decreased in this trend.

279

The degree of fluorescence quenching depends upon the hydrophobic/hydrophilic

280

nature and the origin of HA. The relative bigger Km of the fulvic-like and humic-like

281

fractions of HA might be attributed to the high contents of phenolic and aromatic

282

carboxylic groups, which could form highly stable ring structures with copper. These

283

results suggested that the 2DCOS was a promising tool to readily capture subtle

284

sequential differences in the spectral variations as well as to probe the structural

285

responses with external factors.

286

Significance of This Work. In the present work, a comprehensive picture about

287

the structural changes of DOM in copper binding process is established according to

288

the 2DCOS results at molecular scales. As copper ions are added to interact with HA,

289

the carboxyl and polysaccharide groups give the fastest responses, then the

290

fluorescent humic-like groups participate in the fluorescence queching, followed by

291

vibrational changes of related groups, i.e., phenolic and aryl groups. Subsequently, the

292

involvment of aromatic carboxyl groups further induces the fluorescence quenching

293

of fulvic-like fraction. Finally, the small amounts of amide and aliphatic groups

13 ACS Paragon Plus Environment

Environmental Science & Technology

294

become participated in the copper binding process, after which the fluorescence of

295

protein-like fraction decreases. Hetero-2DCOS offers an unique insight to understand

296

the transformation of HA in the presence of metals in aquatic systems and clarify the

297

structural relationship between the fluorescent groups and the corresponding IR

298

groups in HA in the copper binding process. This study might open a door to visit the

299

geochemical cycling process of metals in natural environments. On the other hand,

300

elucidation of the nature of interactions between HA and copper shown in this work

301

demonstrates the great potential of using the 2D-heterospectral correlation

302

spectroscopy for further applications in exploring complex interactions in natural and

303

engineered environments.

304 305 306 307

AUTHOR INFORMATION *Corresponding author: Prof. Han-Qing Yu, Fax: +86 551 63601592; E-mail: [email protected]

308 309

ACKNOWLEDGEMENTS

310

We thank the Natural Science Foundation of China (51129803), Hefei Center for

311

Physical Science and Technology (2012FXZY005) and the Program for Changjiang

312

Scholars and Innovative Research Team in University and the Collaborative

313

Innovation Center of Suzhou Nano Science and Technology of the Ministry of

314

Education of China for the support of this study.

315 316

ASSOCIATED CONTENT

317

Supporting Information Available. Elemental contents of the unpurified and

318

purified HA samples (Table S1), sign of each cross-peak in the synchronous (Ф) and

14 ACS Paragon Plus Environment

Page 14 of 28

Page 15 of 28

Environmental Science & Technology

319

asynchronous (Ψ, in the brackets) maps from the synchronous fluorescence spectra of

320

HA with copper binding (Table S2), excitation-emission matrix of HA samples and

321

synchronous spectra recorded with different ∆λ values (Figure S1), FTIR spectra for

322

HA samples (Figure S2) and possible assignment for each band. This information is

323

available free of charge via the Internet at http://pubs.acs.org/.

324 325

REFERENCES

326 327 328 329 330

(1) Nebbioso, A.; Piccolo, A. Molecular characterization of dissolved organic matter (DOM): A critical review. Anal. Bioanal. Chem. 2013, 405, (1), 109-124. (2) Stevenson, F. J., Humus Chemistry: Genesis, Composition, Reactions. John Wiley & Sons, Inc.: New York, 1994.

331

(3) Abdulla, H. A. N.; Hatcher, P. G. Dynamics of dissolved organic matter: A view

332

from two dimensional correlation spectroscopy techniques. J. Mol. Struct. 2014,

333

1069, (0), 313-317.

334 335

(4) Cabaniss, S. E. Forward modeling of metal complexation by NOM: II. Prediction of binding site properties. Environ. Sci. Technol. 2010, 45, (8), 3202-3209.

336

(5) Cabaniss, S. E. Forward modeling of metal complexation by NOM: I. A priori

337

prediction of conditional constants and speciation. Environ. Sci. Technol. 2009, 43,

338

(8), 2838-2844.

339

(6) Yan, M.; Korshin, G. V. Comparative examination of effects of binding of

340

different metals on chromophores of dissolved organic matter. Environ. Sci.

341

Technol. 2014, 48, (6), 3177-3185.

342

(7) Croué, J. P.; Benedetti, M. F.; Violleau, D.; Leenheer, J. A. Characterization and

343

copper binding of humic and nonhumic organic matter isolated from the south

15 ACS Paragon Plus Environment

Environmental Science & Technology

344

platte river: Evidence for the presence of nitrogenous binding site. Environ. Sci.

345

Technol. 2002, 37, (2), 328-336.

346

(8) Yan, M.; Benedetti, M. F.; Korshin, G. V. Study of iron and aluminum binding to

347

Suwannee River fulvic acid using absorbance and fluorescence spectroscopy:

348

Comparison of data interpretation based on NICA-Donnan and Stockholm humic

349

models. Water Res. 2013, 47, (14), 5439-5446.

350

(9) Hur, J.; Lee, B.-M. Characterization of binding site heterogeneity for copper

351

within dissolved organic matter fractions using two-dimensional correlation

352

fluorescence spectroscopy. Chemosphere 2011, 83, (11), 1603-1611.

353

(10) Yamashita, Y.; Jaffé, R. Characterizing the interactions between trace metals

354

and dissolved organic matter using excitation−emission matrix and parallel factor

355

analysis. Environ. Sci. Technol. 2008, 42, (19), 7374-7379.

356

(11) Park, J.-H. Spectroscopic characterization of dissolved organic matter and its

357

interactions with metals in surface waters using size exclusion chromatography.

358

Chemosphere 2009, 77, (4), 485-494.

359

(12) Bai, Y. C.; Wu, F. C.; Liu, C. Q.; Li, W.; Guo, J. Y.; Fu, P. Q.; Xing, B. S.;

360

Zheng, J. Ultraviolet absorbance titration for determining stability constants of

361

humic substances with Cu(II) and Hg(II). Anal. Chim. Acta 2008, 616, (1), 115-

362

121.

363

(13) Yan, M.; Dryer, D.; Korshin, G. V.; Benedetti, M. F. In situ study of binding of

364

copper by fulvic acid: Comparison of differential absorbance data and model

365

predictions. Water Res. 2013, 47, (2), 588-596.

366

(14) Xiong, J.; Koopal, L. K.; Tan, W.; Fang, L.; Wang, M.; Zhao, W.; Liu, F.;

367

Zhang, J.; Weng, L. Lead binding to soil fulvic and humic acids: NICA-Donnan

16 ACS Paragon Plus Environment

Page 16 of 28

Page 17 of 28

Environmental Science & Technology

368

modeling and XAFS spectroscopy. Environ. Sci. Technol. 2013, 47, (20), 11634-

369

11642.

370

(15) Gustafsson, J. P. Modeling the acid–base properties and metal complexation of

371

humic substances with the Stockholm Humic Model. J. Colloid Interface Sci.

372

2001, 244, (1), 102-112.

373

(16) Chen, W.; Qian, C.; Liu, X.-Y.; Yu, H.-Q. Two-dimensional correlation

374

spectroscopic analysis on the interaction between humic acids and TiO2

375

nanoparticles. Environ. Sci. Technol. 2014.

376

(17) Xu, H.; Yu, G.; Yang, L.; Jiang, H. Combination of two-dimensional correlation

377

spectroscopy and parallel factor analysis to characterize the binding of heavy

378

metals with DOM in lake sediments. J. Hazard. Mater. 2013, 263, Part 2, (0),

379

412-421.

380

(18) Nakashima, K.; Xing, S.; Gong, Y.; Miyajima, T. Characterization of humic

381

acids by two-dimensional correlation fluorescence spectroscopy. J. Mol. Struct.

382

2008, 883–884, (0), 155-159.

383

(19) Kim, H. J.; Kim, S. B.; Kim, J. K.; Jung, Y. M. Two-dimensional heterospectral

384

correlation analysis of wide-angle X-ray scattering and infrared spectroscopy for

385

specific chemical interactions in weakly interacting block copolymers. J. Phy.

386

Chem. B 2006, 110, (46), 23123-23129.

387

(20) Mee Jung, Y.; Noda, I. New approaches to generalized two ‐ dimensional

388

correlation spectroscopy and its applications. Appl. Spectrosc. Rev. 2006, 41, (5),

389

515-547.

390

(21) Zou, C.; Larisika, M.; Nagy, G.; Srajer, J.; Oostenbrink, C.; Chen, X.; Knoll, W.;

391

Liedberg, B.; Nowak, C. Two-dimensional heterospectral correlation analysis of

392

the redox-induced conformational transition in cytochrome c using surface-

17 ACS Paragon Plus Environment

Environmental Science & Technology

393

enhanced raman and infrared absorption spectroscopies on a two-layer gold

394

surface. J. Phy. Chem. B 2013, 117, (33), 9606-9614.

395

(22) Koeppe, B.; Tolstoy, P. M.; Guo, J.; Nibbering, E. T. J.; Elsaesser, T. Two-

396

dimensional UV–vis/NMR correlation spectroscopy: A heterospectral signal

397

assignment of hydrogen-bonded complexes. J. Phy. Chem. Lett. 2011, 2, (9),

398

1106-1110.

399

(23) Choi, H. C.; Ryu, S. R.; Ji, H.; Kim, S. B.; Noda, I.; Jung, Y. M. Two-

400

dimensional heterospectral correlation analysis of X-ray photoelectron spectra and

401

infrared spectra for spin-coated films of biodegradable poly(3-hydroxybutyrate-

402

co-3-hydroxyhexanoate) copolymers. J. Phy. Chem. B 2010, 114, (34), 10979-

403

10985.

404 405 406 407

(24) Noda, I.; Ozaki, Y., Two-Dimensional Correlation Spectroscopy - Applications in Vibrational and Optical Spectroscopy. John Wiley & Sons Ltd: London, 2004. (25) Noda, I. Close-up view on the inner workings of two-dimensional correlation spectroscopy. Vib. Spectrosc. 2012, 60, (0), 146-153.

408

(26) Berkovic, A. M.; Garcia Einschlag, F. S.; Gonzalez, M. C.; Pis Diez, R.; Martire,

409

D. O. Evaluation of the Hg2+ binding potential of fulvic acids from fluorescence

410

excitation-emission matrices. Photoch. Photobio. Sci. 2013, 12, (2), 384-392.

411

(27) Hur, J.; Jung, K.-Y.; Jung, Y. M. Characterization of spectral responses of humic

412

substances upon UV irradiation using two-dimensional correlation spectroscopy.

413

Water Res. 2011, 45, (9), 2965-2974.

18 ACS Paragon Plus Environment

Page 18 of 28

Page 19 of 28

Environmental Science & Technology

Table 1. 2D-IR-COS Results on the Assignment and Sign of Each Cross-Peak in Synchronous (Ф) and Asynchronous (Ψ, in the brackets) Maps of HA with Copper Binding position

sign

assignment

(cm-1) 1720

carboxyl vC=O

1660

amide vC=O

1620

aromatic vC=C

1400

phenolic vC-O, δO-H

1350

aliphatic δC-H

1100

polysaccharides vC-O

850

Ar-H, vC-H

1720 1660 1620

1400

1350

1100

850

+

+(+)

+(+)

+(+)

+(+)

+(+)

-(-)

+

+(-)

+(-)

+(+)

+(-)

-(-)

+

+(-)

+(+)

+(-)

-(-)

+

+(+)

+(-)

-(-)

+

+(-)

-(-)

+

-(-) +

(signs were obtained in the upper-left corner of the maps, + positive, - negative)

19 ACS Paragon Plus Environment

Environmental Science & Technology

Page 20 of 28

Table 2. Copper Binding Parameters Calculated by the Modified SternVolmer Equation fraction

peak (nm)

f

LogKM

R*

protein-like

292

0.542

4.61

0.987

350

0.769

4.99

0.994

380

0.771

5.09

0.992

445

1

5.26

0.999

fulvic-like

fumic-like *

Correlation coefficients (R).

20 ACS Paragon Plus Environment

Page 21 of 28

Environmental Science & Technology

Figure captions

Figure 1. Changes in the synchronous fluorescence spectra of 10 mg/L HA at pH 7.0 upon addition of copper with concentrations of 0, 1, 2, 5, 10, 20, 40, 60, 80, 100, 120, 160 µM. The dashed vertical lines shows the protein-like, the fulvic-like, and the humic-like regions. Each spectrum was averaged by three scans

Figure 2. FTIR absorbance spectra (a) and differential IR absorbance spectra (b) of the freeze-dried HA samples with varying copper concentrations. The spectra were an average of 32 scans

Figure 3. Synchronous (a) and asynchronous (b) 2D correlation maps generated from the synchronous fluorescence spectra of HA with increasing copper ion concentration in the region 280-500 nm

Figure 4. Synchronous (a) and asynchronous (b) 2D IR correlation maps from a FTIR analysis of HA with copper binding

Figure 5. 2D heterospectral correlation plots generated from the dynamic synchronous fluorescence and IR spectra of HA upon copper addition: (a) synchronous map; (b) asynchronous map

Figure 6. Modified Stern-Volmer plots for the fluorescence quenching of HA fractions by copper. Four wavelengths at 292, 350, 380 and 445 nm were selected corresponding to the peaks in the 2DCOS map in Figure 3a

21 ACS Paragon Plus Environment

Environmental Science & Technology

80

Proteinlike

Fulviclike

Humiclike

60

Intensity (a.u.)

Page 22 of 28

2+

[Cu ] increase

40

20

0 250

300

350

400

450

Wavelength (nm) Figure 1

22 ACS Paragon Plus Environment

500

550

Page 23 of 28

Environmental Science & Technology

2+

[Cu ] 160 µM 120 µM 100 µM 80 µM 60 µM 40 µM 20 µM 10 µM 5 µM 2 µM 1 µM 0 µM

IR Absorbance

(a)

1750

1500

1250

1000 -1

750

Differential Absorbance

Wavenumber (cm ) (b)

1750

160 µM 120 µM 100 µM 80 µM 60 µM 40 µM 20 µM 10 µM 5 µM 2 µM 1 µM

1500

1250

1000 -1

Wavenumber (cm )

Figure 2

23 ACS Paragon Plus Environment

750

Environmental Science & Technology

Figure 3

24 ACS Paragon Plus Environment

Page 24 of 28

Page 25 of 28

Environmental Science & Technology

Figure 4

25 ACS Paragon Plus Environment

Environmental Science & Technology

Figure 5

26 ACS Paragon Plus Environment

Page 26 of 28

Page 27 of 28

Environmental Science & Technology

15

F0/(F0-F)

10

292 nm 350 nm 380 nm 445 nm Modified Stern-Volmer model

5

0 0.0

0.2

0.4 2+

-1

1/[Cu ] (µM )

Figure 6

27 ACS Paragon Plus Environment

0.6

Environmental Science & Technology

Table of Contents (TOC) art

28 ACS Paragon Plus Environment

Page 28 of 28