Metallothionein-Like Multinuclear Clusters of Mercury(II) and Sulfur in

Aug 2, 2011 - Department of Civil, Environmental, and Architectural Engineering, University of Colorado Boulder, 428 UCB, Boulder,. Colorado 80309-042...
4 downloads 9 Views 1MB Size
ARTICLE pubs.acs.org/est

Metallothionein-Like Multinuclear Clusters of Mercury(II) and Sulfur in Peat Kathryn L. Nagy,*,† Alain Manceau,‡ Jarrod D. Gasper,§,^ Joseph N. Ryan,|| and George R. Aiken§ †

)

Department of Earth and Environmental Sciences, MC-186, 845 West Taylor Street, University of Illinois at Chicago, Chicago, Illinois 60607, United States ‡ ISTerre-Maison des Geosciences, CNRS and Universite J. Fourier, BP 53, 38041 Grenoble Cedex 9, France § United States Geological Survey, Water Resources Discipline, 3215 Marine Street, Boulder, Colorado 80303, United States Department of Civil, Environmental, and Architectural Engineering, University of Colorado Boulder, 428 UCB, Boulder, Colorado 80309-0428, United States

bS Supporting Information ABSTRACT: Strong mercury(II) sulfur (Hg-SR) bonds in natural organic matter, which influence mercury bioavailability, are difficult to characterize. We report evidence for two new Hg-SR structures using X-ray absorption spectroscopy in peats from the Florida Everglades with added Hg. The first, observed at a mole ratio of organic reduced S to Hg (Sred/Hg) between 220 and 1140, is a Hg4Sx type of cluster with each Hg atom bonded to two S atoms at 2.34 Å and one S at 2.53 Å, and all Hg atoms 4.12 Å apart. This model structure matches those of metal thiolate clusters in metallothioneins, but not those of HgS minerals. The second, with one S atom at 2.34 Å and about six C atoms at 2.97 to 3.28 Å, occurred at Sred/Hg between 0.80 and 4.3 and suggests Hg binding to a thiolated aromatic unit. The multinuclear Hg cluster indicates a strong binding environment to cysteinyl sulfur that might impede methylation. Along with a linear Hg(SR)2 unit with Hg—S bond lengths of 2.34 Å at Sred/Hg of about 10 to 20, the new structures support a continuum in Hg-SR binding strength in natural organic matter.

’ INTRODUCTION Reduced sulfur (Sred) ligands are the main pool of strong binding sites for mercury(II) (Hg) in natural organic matter (OM) that may influence methylation rates, distribution and transport, and ecotoxicity of mercury. The Sred derives from biomolecules produced by organisms, free sulfide species formed during microbial sulfate reduction, and recycled products of OM degradation and diagenesis. Total amounts of Sred in OM range typically from about 0.1 2 wt %, sufficient to bind even unusually high concentrations of Hg. Small organic and inorganic Hg Sred species can enable methylation,1,2 but effects of Hg bound to larger OM molecules on microbial methylation and direct toxicity to fish and filterfeeders are ambiguous or unknown.3 5 Unlike small biomolecules, which are labile and can rapidly exchange thiol-bonded Hg,6 macromolecular OM may have binding environments that are influential in decreasing bioavailability of Hg. Because peat and soil OM contain more macromolecular plant-derived compounds, as well as more reduced and total S, compared to dissolved OM (DOM), they may sequester Hg from the aqueous phase. The binding strength of Hg to OM, expressed as conditional HgLR+ and HgL2R formation constants, spans from about 1010 to 1038,7 14 corresponding to mole ratios of Sred to bound Hg (Sred/Hg) from less than 1 to about 77000 in the solid phase at equilibrium with dissolved mercury. In the form Hg(LR)2 the r 2011 American Chemical Society

largest constants are g1043, matching values for thiol-binding in low molecular weight organic molecules.15 The strongest binding sites containing Sred are fewer in number than the markedly weaker sites containing O and/or N.8 Binding constants must be linked to molecular structures to quantify reactions between OM and Hg species. Such structures have been characterized in peat, soil OM, and humic acid with added metal using extended X-ray absorption fine-structure (EXAFS) spectroscopy. At mole ratios of Sred/Hg = 0.3 to 1.2 a structure with one S and one O or N bond (i.e., R(O/N)HgSR) was observed16,17 whereas at Sred/Hg = 1.6 to 4, 60 to 90% of added Hg bonded to 1.3 to 1.7 S atoms, suggesting up to 70% Hg(SR)2 in a linear configuration.17 At Sred/Hg = 20 to 100 a Hg(SR)2 3 complex with two S ligands at 2.34 Å and a tentative third near 3 Å in a T-shape geometry was reported.18,19 Possible three-coordinate binding is supported by the average release of 2.7 protons for each adsorbed Hg atom in peat humic acid at Sred/Hg = 120.13 Peat typically contains about 0.1 ppm Hg (0.1 mg Hg/kg OM, equivalent to about 106 mol Sred/mol Hg). This amount of Received: March 27, 2011 Accepted: July 14, 2011 Revised: July 12, 2011 Published: August 02, 2011 7298

dx.doi.org/10.1021/es201025v | Environ. Sci. Technol. 2011, 45, 7298–7306

Environmental Science & Technology mercury is about 1% of the adsorbed amount (10 ppm Hg) at which the strongest binding constants in OM have been measured,8,9,12 and about 0.1% of the lowest value at which EXAFS spectra are measurable today. Because only about 0.02% of Sred is needed to bind 0.1 ppm Hg, molecular structures identified by EXAFS at around 100 ppm Hg, or about 20% of the Sred, may reasonably approximate the strongest sites. Our goal was to distinguish binding environments that involved Sred by obtaining and interpreting EXAFS spectra of Hg-containing peats, including spectra obtained at liquid He temperature. The peats were collected from two locations (F1 and 2BS) in the Everglades that had different inorganic sulfide (HS ) concentrations and Hg methylation rates.20 They were shown to have dominantly strong to weak binding at adsorbed Hg concentrations over the range of 1 400000 ppm.9 We observed two new molecular binding environments out of a total of four between 59 and 99100 ppm added Hg (Sred/Hg mole ratio from 1140 to 0.50) that expand our understanding of how the continuum of sulfur binding sites in OM may form and influence bioavailability of mercury.

’ EXPERIMENTAL SECTION We adsorbed mercury to the peats using the experimental procedure described previously9 except that the concentration of peat was 10 times higher (0.1 g peat in 0.25 L of solution) to obtain enough material for EXAFS spectroscopic analysis. Experiments were conducted at initial solution concentrations of mercury from 3  10 7 to 2  10 4 mol/L. Mercury was analyzed using cold vapor atomic fluorescence spectroscopy (PS Analytical Ltd., Kent, UK). EXAFS spectra at the Hg-L3 edge were recorded at room temperature (RT) and nine months later at liquid helium temperature (He-T; 8 16 K) at beamline BM30B (FAME) of the European Synchrotron Radiation Facility (ESRF) and analyzed by standard techniques. X-ray absorption near-edge structure (XANES) spectra at the S-K edge were recorded on unreacted peat at RT at beamline 10.3.2, Advanced Light Source (ALS). Fractional abundances of reduced and oxidized sulfur species were obtained by linear combination fitting of reference spectra. Binding structures for Hg were identified by comparing spectral features of samples to those of reference compounds, fitting sample spectra with linear combinations of spectra, and modeling Fourier-transformed k3-weighted spectra using amplitude and phase shift functions calculated from first principles. Details of the experiments and analyses are in the Supporting Information (SI). ’ RESULTS AND DISCUSSION In contrast to the generally linear trend in the adsorption isotherms, which may be explained by at most two binding sites, the EXAFS spectral analysis revealed four structural binding environments, three of which involved sulfur ligands. Two of the binding environments identified are reported for the first time. One is a multinuclear Hg-SR coordination structure observed at low amounts of adsorbed Hg and the other is indicative of Hg bonded to thiolated aromatic subunits observed at high amounts of adsorbed Hg. We first present results for the adsorption isotherms and XANES spectra of the sulfur speciation in the peats, followed by results for the individual binding environments derived from the EXAFS spectra.

ARTICLE

Figure 1. S-XANES spectra for F1 and 2BS peats. The peats were fit with approximately the same amounts of cystine (29 31%), sulfoxide (9 11%), and sulfate (6 7%), but F1 had more cysteine (41 vs 34%) and 2BS more sulfonate (23 vs 16%). The first difference is seen in the ratio of the peak amplitudes of reduced (Sred) to oxidized sulfur (Sox), which is higher for F1 than 2BS, and the second in the asymmetrical shape of the oxidized sulfur peak for 2BS which crests at the energy of the sulfonate component (2480.5 eV).

S-XANES Spectra and Adsorption Isotherms. Initial values of the mole ratio Sred/Hg in the peats were 3  105 for 2BS and 8  105 for F1, where Sred was determined by XANES spectroscopy (Figure 1 and SI). The peats adsorbed 1 to 99100 ppm Hg (Sred/Hg = 78300 to 0.50, respectively) as aqueous Hg concentration increased, similarly to earlier observations9 (SI). At low aqueous concentrations of Hg the F1 peat, which had more total and reduced S (1.79 wt % STotal; 1.25 wt % Sred), adsorbed slightly more Hg than the 2BS peat (1.27 wt % STotal; 0.83 wt % Sred). One inflection point between 10 8 and 10 7 mol Hg/L indicates the data might be fit with two Langmuir isotherms representing two binding sites. EXAFS Spectra at Room Temperature. The k3-weighted spectra of peats with 59 99100 ppm Hg (Sred/Hg = 1140 to 0.50) obtained at RT have features that differentiate three binding environments (Figure 2a). Two are evident by isosbestic points for samples having 4170 to 99100 ppm Hg, resulting from a rightward shift in oscillation frequency and a decrease in amplitude with increasing Hg concentration (Figure 2b). These features show that, with decreasing Sred/Hg mole ratio, Hg binding transitions from dominant longer bonds (short wave frequency) with S atoms to shorter bonds (long wave frequency) with atoms of smaller ionic radii (O, N).17,18 At any amount of Hg, the F1 spectrum is shifted slightly left of the 2BS spectrum to the point that F1-99 (99100 ppm) is more like 2BS-19 (19000 ppm) than 2BS-99 (99100 ppm). The systematic shift is related to the 50% more Sred in the F1 peat and confirms the significance of Hg S bonds in the EXAFS signal even at high amounts of adsorbed Hg (i.e., low values of Sred/Hg).16,17 A third binding environment is indicated at lower Hg concentrations by slight leftward shifts of the oscillations, which imply even longer bonds (Figure 2c). This observation motivated measurements at He-T of six samples representing the variation in RT spectra and the spread of Sred/Hg values. Freezing slows the thermal motion of atoms, and for an 7299

dx.doi.org/10.1021/es201025v |Environ. Sci. Technol. 2011, 45, 7298–7306

Environmental Science & Technology

ARTICLE

Figure 2. EXAFS spectra for F1 and 2BS peat samples. Sample code suffixes represent adsorbed Hg in ppm (mg Hg/kg peat): 006 and 007 (59 and 69 ppm); 03 (∼350 ppm); 4 (∼4100 ppm); 18 and 19 (18400 and 19000 ppm); 99 (99100 ppm) (SI). (a) Room temperature (RT) spectra cut at k = 10 Å 1 for clarity except for 006 and 007 spectra which were stopped at k = 9.8 Å 1 because of interfering absorbance at the Se K-edge. (b) RT spectra at k = 6 to 8 Å 1 showing rightward shift to lower frequency at higher Hg concentration, signifying shorter Hg-ligand distances. Spectra at 59 to 350 ppm are omitted for clarity. (c) RT spectra at k = 4.2 to 5.6 Å 1 showing leftward shift to higher frequency at lower Hg concentration, signifying longer Hg-ligand distances. (d) Spectra at liquid helium temperature (He-T) cut to k = 10 Å 1.

individual bond increases the amplitude of the EXAFS signal. Thus, we expected an improved signal-to-noise ratio and better resolution of interatomic distances at liquid He temperature. EXAFS Spectra at Liquid He Temperature. At He-T, longer Hg—S bonds were more evident by the farther leftward shift of spectra at low Hg concentrations and sinistral asymmetry of the first two oscillations (Figures 2d and 3a). Other bonds, also more apparent at He-T, caused extra small oscillations at about k = 7.8 and 9.2 Å 1 in F1-007 (69 ppm Hg; Sred/Hg = 1140) and F1-03 (353 ppm Hg; Sred/Hg = 220) (Figure 2d) and also at 10.6 and 12 Å 1 in F1-03 (SI). Detailed analyses of the spectra and modeled coordination structures from low to high adsorbed Hg concentrations (i.e., high to low mole ratios of Sred/Hg) are presented below. Multinuclear Hg—S Binding at Low Hg Concentration. The spectra of F1-03 and F1-007 changed the most between RT and He-T, implicating dominant soft metal-soft ligand Hg—S bonds. Fourier transforms of the F1-03 spectra show that (1) the peak at R + δ ≈ 2 Å (distance uncorrected for phase shift) with a small shoulder on the right at RT split into two peaks at about 1.8 and 2.3 Å and (2) the peak at about 4 Å increased in amplitude at He-T (Figure 3a). These changes are understood by comparison

to the changes in cinnabar (R-HgS) spectra between RT and He-T (Figure 3b). At RT the spectrum of R-HgS is dominated by the signal from Hg S bonds of length 2.38 Å (R + δ ≈ 2 Å); the signal from Hg Hg pairs at R + δ ≈ 4 Å is weak. At He-T the signal amplitude is higher and the local structure of R-HgS fully agrees with crystallographic data. Longer Hg—S bonds at 3.07 and 3.26 Å are detected, modeling shows two distinct Hg Hg pairs at 3.79 Å and 4.11 Å, and longer Hg Hg pairs at 5.67 and 6.82 Å appear (SI). Similar spectral features are seen in cinnabar spectra obtained at 77 K21 (SI). Changes for the F1 peats with temperature, although similar, are distinct from those for cinnabar. Differences include the outof-phase patterns between k ≈ 10.2 and 12.8 Å 1 and the peak shapes at R + δ ≈ 4 Å at He-T. Using the cinnabar spectrum at He-T to guide bond assignments, two Hg S bonds of length 2.34 Å and one of length 2.53 Å were derived from the F1-03 peat spectrum at He-T along with one Hg—C or Hg—S bond at 3.25 Å, and three Hg Hg pairs at 4.12 Å (Figure 4a, SI). The counterintuitive weakening of spectral amplitude at low k values in F1-03 at He-T compared to RT is a result of the relative difference in enhancement of amplitude of the two Hg—S bonds at low temperature (SI). 7300

dx.doi.org/10.1021/es201025v |Environ. Sci. Technol. 2011, 45, 7298–7306

Environmental Science & Technology

ARTICLE

Figure 3. EXAFS spectra (left) and Fourier transforms (right) for F1 03 peat (353 ppm Hg) and cinnabar (R-HgS). Imaginary parts of the Fourier transforms show that the differences in the EXAFS frequencies observed in k-space at RT and He-T come from the magnification at low temperature of distant Hg—S and Hg—Hg pairs. (a) F1-03 at RT and He-T. (b) R-HgS at RT and He-T. Amplitudes for the spectra and Fourier transforms of R-HgS were reduced by multiplicative factors to aid comparison of peak shapes and positions.

Comparison of the Hg-SR structure in the peat and metacinnabar (β-HgS) at He-T is also warranted, but is not straightforward. Although the spectra for metacinnabar and cinnabar (R-HgS) are distinct at RT, the spectrum of β-HgS is similar to that of cinnabar (R-HgS) at He-T, indicating that the structure of the metastable mercuric sulfide (β-HgS) approaches that of cinnabar (R-HgS) near absolute zero (SI). Such a structural transition is consistent with the relative thermodynamic stability of the two phases: metacinnabar is the stable phase at higher temperature.22 The transition is also consistent with a reduction in crystallographic symmetry from cubic (β-HgS) to trigonal (R-HgS) as temperature is lowered. A similar transition was reported recently for magnetite which changed structure from cubic to monoclinic at about 125 K.23 The spectrum of the structurally transformed metacinnabar does not match that of F1-03 at He-T (SI). We used the bond distances derived for F1-03 at He-T as initial conditions for modeling the RT F1-03 spectrum. A good fit was obtained with two S atoms at 2.35 Å, 0.6 S at 2.57 Å, and 1.6 Hg at 4.14 Å (SI); the 3.25 Å distance was not detected. Like for cinnabar, fewer longer Hg—S and Hg—Hg distances were revealed at RT than at He-T, but none changed significantly (SI). In addition, F1-03 does not match either metacinnabar, which has all Hg—S bonds at 2.51 Å, or cinnabar at RT (SI). Freezing did not change the unique binding between Hg and S in the peat: the multinuclear Hg-SR structure, which is nanosized and surrounded by a relatively flexible organic matrix, is stable.

Organic materials are commonly examined in the frozen state by EXAFS spectroscopy to avoid formation of free radicals that cause photoreduction.24,25 However, there is virtually no evidence for degraded RT structures during synchrotron measurements of Hg binding to OM,16 18 likely in part because of the strong bonds that form between Hg and Sred. In fact, the small oscillation near 7.8 Å 1 at He-T in F1-007 and F1-03 (Sred/Hg = 1140 and 220) caused by Hg Hg pairs is also perceptible in RT spectra of soil OM at Sred/Hg = 100 and 20.18 This is not only further evidence against photoreduction of strongly bound Hg at RT but also supports the conclusion that freezing did not induce the multinuclear structure. The features of Hg Hg pairs observed in the He-T spectra of the peats may be visible in the spectra of soil OM at RT18 because the experimental soil OM Hg solutions equilibrated for 1 wk18 compared to 15 h in the peat OM-Hg solutions. Mercury atoms could have had time to move to higher affinity protein-like Sred sites even though the soil OM had 75% less Sred than the Everglades peats. A similar kinetic effect may have contributed to visualizing the longer bond distances in our He-T analyses given the nine-month interval between the measurements at the two temperatures. In addition, the greater amount of Sred in the peats may have accelerated the change in binding enough to detect the multinuclear cluster in the modeled RT spectra despite the short equilibration time. The final pH of our peat experiments was 5 6 vs a pH near 4 for the soil OM.18 If apparent pKa values of thiols in peats are in the range of those for plant 7301

dx.doi.org/10.1021/es201025v |Environ. Sci. Technol. 2011, 45, 7298–7306

Environmental Science & Technology

ARTICLE

Figure 4. EXAFS spectra, Fourier transforms, and model structures of dominant Hg binding structures (left to right) and correspondence with Sred/Hg in the peats (far right). (a) F1 03 (353 ppm Hg) at RT and He-T with model fits. Schematic multinuclear Hg4Sx structure is at right. Data analysis shows that the Hg—S distance at 2.53 2.57 Å contributes to the right tail of the first peak (SI). Imaginary parts of the Fourier transforms were omitted for clarity. (b) F1-4 (~4100 ppm Hg) and Hg(Cys)2 reference solution at He-T. The modeled Hg(SR)2 structure is shown at the upper right. (c) F1-99 (99100 ppm) and 2BS-99 (99100 ppm) at He-T. A schematic of the five-membered (O,N) chelate ring structure in 2BS-99 is shown at the lower right. Modeled bond distances and coordination numbers in F1-99 suggest the hypothetical structure linking Hg to aromatic rings at the upper right. Atoms coded as oxygen are indistinguishable from nitrogen.

metallothioneins,26 then accelerated binding of Hg in the peat also could be explained by more deprotonated sites at higher pH. Similarity of Multinuclear Hg S Structure to Metal Thiolate Clusters in Metallothionein Proteins. The average Hg S structure can be described as a tetrameric cluster (Hg4Sx) in which three-coordinate molecular units of one Hg atom with two short bonds (2.34 Å) and one longer bond (2.53 Å) to S atoms are arranged so that any two Hg atoms are about 4.12 Å apart (Figure 4a). The bonds around any Hg atom form a Y-shape,19 precluding successful modeling of the peak at R + δ ≈ 4 Å in the Fourier transforms at both temperatures as multiple-scattering from a linear Hg(SR)2 complex18 (SI). At least 0.5 mol % of Sred must be in the cluster, similar to the 0.15 1% deduced for strong binding constants of Hg in other OM samples.8,9,15 The Hg4Sx model cluster resembles polynuclear-thiolate clusters in Cu(I) metalloproteins and metal-peptide compounds.27

Cu(I)4S6-type structures with three-coordinate metal centers occur in metallochaperone and metallothionein (MT) proteins with CXC (CysXaaCys) amino acid consensus sequences from mammals and yeast.28 32 These structures mimic organometallic [Cu4(SPh)6]2-, which contains thiolated aromatics.33 Structures of Hg metalloproteins are few and varied. One three-coordinate unit with Hg—S bond length of 2.43(2) Å was observed in the merR-regulatory protein,34 and one with either two 2.34 Å and one 2.47 Å or three 2.39 Å Hg—S bonds in a de novo Hg(II)-metalloprotein.35 In Hg7-MT, Hg is bonded to either two S atoms at 2.33(2) Å (and possibly two other S at 3.4 Å)36 or to three S at 2.42(2) Å,37 whereas in Hg18-MT Hg is bonded to two S atoms at 2.42(3) Å and one Cl at 2.57 Å;36 no Hg Hg distances were reported for either Hg-MT. The crystal structure of the metallochaperone Hah1 in the presence of Hg has a distorted tetrahedral coordination with one Hg bonded to 7302

dx.doi.org/10.1021/es201025v |Environ. Sci. Technol. 2011, 45, 7298–7306

Environmental Science & Technology one S at 2.3 Å, two S at 2.5 Å, and one S at 2.8 Å.38 One Hg atom is thought to bind to two Hah1 molecules each with CXXC sequences. Molecules with CXXC sequences usually form twocoordinate structures as in the MerP Hg(II) transport protein.39 Three-coordinate binding in a multinuclear cluster was inferred for Hg-loaded MT from Cicer arietinum (chickpea) using UV-vis spectroscopy.40 Hg(SR)2, Chelate Ring, and Thiolated Aromatic Structures at Higher Hg Concentrations. A linear Hg(SR)2 structure was dominant in the F1-4 peat with 4180 ppm Hg (Sred/Hg = 18.7). This structure has the same shorter Hg—S bonds (2.33 Å at HeT and 2.35 Å at RT) as the cluster, but not the longer Hg—S bond (2.53 Å) and Hg Hg distance (4.12 Å): there is neither a spectral shift to lower k nor extra oscillations at higher k (Figure 4b). The spectrum and Fourier transform at He-T match those for our 1:2 Hg-cysteine solution (Hg(Cys)2), and similarly frozen Hg-cysteine solids,41 and Hg-cysteine and Hg-glutathione solutions.41,42 Linearity is confirmed by the peak at R + δ ≈ 4 Å in the Fourier transform at He-T, caused by multiple scattering paths (SI). Only 65% of the Hg(Cys)2 spectrum is needed to match the F1 4 spectrum indicating that at least 20 25% of the Hg is bonded to O and/or N (Figure 4b). The other sample having a dominant binding environment was 2BS-99 (99100 ppm Hg and Sred/Hg = 0.50). The spectra indicate shorter bonds between Hg and lower Z atoms (Figures 2 and 4c). In the best model, Hg bonds to one O atom at 2.02 Å, one O or N atom at 2.45 Å, one C atom at 2.87 Å, and one C atom at 3.28 Å (SI). The structure is likely a five-membered chelate ring with one carboxylate O and amino N (R-HN2) or alcohol O (R-OH) from an R-substituted carbon, similar to structures for monomethylmercury binding to amino acids43,44 and for Hg and Cu(II) complexed to soil OM and humic acid.18,24 The second new Hg S structure appeared in F1-99 (also having 99100 ppm Hg, but a higher Sred/Hg of 0.80 compared to 2BS-99). A small oscillation at about 6.7 Å 1 and the shift of the third oscillation at 8 Å 1 in the He-T spectrum are distinctive (Figure 4c). The best model gives an average structure around Hg of one O or N atom at 2.02 Å, 0.4 S atoms at 2.34 Å, two C atoms at 2.97 Å, and four C atoms at 3.28 Å (SI). The number of C atoms along with S suggests the presence of thiolated aromatic units.45 However, the low amount of S indicates that this structure is mixed mainly with the (O,N) chelate ring seen in 2BS-99. The spectra of F1-18 (18400 ppm Hg; Sred/Hg = 4.3) could not be modeled successfully from first principles because there was no dominant binding environment. Instead, RT and He-T spectra were analyzed by linear combination fitting using spectra from the other peat samples. Sample F1-18 has about 67% F1-4 type bonds and 33% of the mixture of bonds modeled in F1-99 (SI), further supporting a continuum in binding environments. Environmental Implications. As the concentration of adsorbed Hg increased, the Hg atoms filled increasing amounts of progressively weaker binding sites. Binding of highest to lowest affinity corresponds in order to the Hg4Sx cluster, linear Hg(SR)2, thiolated aromatic, and five-membered chelate ring sites (Figure 4). The first three reinforce the importance of thiolate ligands even at relatively low values of Sred/Hg, where OH and N ligands are in great excess. Other relatively weak binding sites with O and N ligands may exist at low values of Sred/ Hg, but could not be distinguished in this study. Also, the EXAFS spectroscopic data obtained here cannot be used to differentiate O from N ligands such as the imidazole component of histidine which can bind metals in protein structures.46 A four-coordinate

ARTICLE

Hg-SR structure like that observed, for example, in the HgHah1 metallochaperone protein38 or the Hg(Cys)4 species in Hg-cysteine solutions with H2Cys/Hg(II) > 541 might form at higher values of Sred/Hg than investigated here. However, the longer bond-length expected for such a structure19 was never dominant in the peat samples. Thiols may be incorporated into OM as components of biomolecules such as metallothioneins and phytochelatins. Metallothioneins, which have high contents of cysteinyl thiols, are produced by bacteria and eukaryoa and possibly by archaea.26,47,48 Photosynthetic organisms also produce phytochelatins, which are glutathione-based peptides with (γ-Glu-Cys)nGly (PCn, n = 2 11). In addition, thiols can be formed when hydrogen sulfide, typically generated by bacterial sulfate reduction,49 51 adds to unsubstituted double bonds in alkenes and aromatics from tannins, lignins, and terpenoids.52,53 Reduced sulfur in the peats has XANES character like that of thiols in cysteine and disulfide bridges in cystine (SI); that is, the reduced sulfur that forms metal thiolate clusters in metalloproteins. Such clusters typically have 2 11 metal centers54with evidence in bacterial, animal, and plant MTs that histidine residues stabilize the negative charge of the cluster.55 Histidine residues can also be part of the coordination structure.46 Values equal to 2 2.7 observed for the ratios of protons exchanged for low amounts of adsorbed Hg in DOM and soil OM at pH e 710,13 imply x = 8 11 in the Hg4Sx structure, in the range of typical numbers of S atoms in MT metal thiolate clusters. Living organisms can control the binding of metals to stable clusters in MTs by cellular redox reactions 56,57 and rapid metal exchange can occur through protein folding and unfolding.58 Metal thiolate clusters in plant MTs tend to be somewhat less stable than those in most vertebrate MTs at comparable pH,26,59 which may be related to fewer cysteine residues and longer linker chains between domains in plant MTs. Upon cellular death, metal thiolate clusters might become incorporated in OM whereas amino acid linker segments likely decompose. Mercury(II) binds more strongly than Cu(I), Cd(II), Pb(II), and Zn(II) to MTs and model peptides in vitro60 62 and in phytochelatins the strength of metal binding is higher with length of the peptides (i.e., with more cysteine).63,64 Phytochelatin metal complexes can be stabilized if plants or fungi produce excess sulfide.65 67 Phytochelatins, MTs, and intracellular sulfide production may all explain Hg-rich granules in roots of plants grown in vivo68 70 and in a contaminated saltmarsh.71 In the Everglades peats, the Sred in the thiolate ligands that formed the experimental Hg4Sx type of clusters may originate as cysteine in MTs or phytochelatins in sawgrass (2BS), cattail (F1) (SI), and periphyton72 or in MTs produced by aquatic invertebrates,73 fungi, and bacteria. Thiols also could have been added to the peats during bacterial reduction of aqueous sulfate.20 The two new Hg-SR structures we identified may elucidate other interactions between OM and Hg species. Formation of mercury-thiolate clusters might explain why metacinnabar was inhibited from precipitating in DOM solutions to which Hg was added first (