Glutathione limits aquacopper(I) to sub-femtomolar concentrations through cooperative assembly of a tetranuclear cluster

The tripeptide glutathione (GSH) is a crucial intracellular reductant and radical scavenger, but it may also coordinate the soft Cu(I) cation and thereby yield pro-oxidant species. The GSH-Cu(I) interaction is thus a key consideration for both redox and copper homeostasis in cells. However, even after nearly four decades of investigation, the nature and stability of the GSH-Cu(I) complexes formed under biologically relevant conditions remain controversial. Here, we revealed the unexpected predominance of a tetranuclear [Cu4(GS)6] cluster that is sufficiently stable to limit the effective free aquacopper(I) concentration to the sub-femtomolar regime. Combined spectrophotometric–potentiometric titrations at biologically realistic GSH/Cu(I) ratios, enabled by our recently developed Cu(I) affinity standards and corroborated by low-temperature phosphorescence studies, established cooperative assembly of [Cu4(GS)6] as the dominant species over a wide pH range, from 5.5 to 7.5. Our robust model for the glutathione-Cu(I) equilibrium system sets a firm upper limit on the thermodynamic availability of intracellular copper that is 3 orders of magnitude lower than previously estimated. Taking into account their ability to catalyze the production of deleterious superoxide, the formation of Cu(I)-glutathione complexes might be avoided under normal physiological conditions. The actual intracellular Cu(I) availability may thus be regulated a further 3 orders of magnitude below the GSH/Cu(I) affinity limit, consistent with the most recent affinity determinations of Cu(I) chaperones.

Glutathione (GSH) 2 is a ubiquitous thiol-containing non-ribosomal tripeptide that plays a central role in cellular redox homeostasis (1)(2)(3), detoxification (4), and signaling pathways (5). Present at millimolar concentrations, GSH serves as a cytosolic redox buffer (6) and is essential for combating oxidative stress. By virtue of its sulfhydryl group, GSH functions not only as an intracellular reductant and radical scavenger but may also serve as a ligand for soft transition metal cations. The formation of complexes with copper ions is of particular interest as GSH plays a critical role in cellular copper homeostasis (7). For example, depletion of GSH in HEK293 cells resulted in a 50% decrease in the initial rate of copper uptake (8). At the same time, recent in vitro experiments demonstrated that such complexes may catalyze the production of superoxide from dioxygen (9), thus challenging the presumed role of GSH as an immediate copper ligand. Because Cu(II) is rapidly converted to Cu(I) within the reducing cellular environment, the coordination chemistry of GSH with Cu(I) is key toward understanding its potential role as a ligand in cellular copper homeostasis. Despite the importance of the GSH-Cu(I) interaction, studies spanning the past four decades have yet to reach a consensus on the nature and stability of the complexes produced under physiological conditions (10 -16). Previously reported stability constants, derived from equilibrium titrations using stoichiometric quantities of Cu(I) and GSH, imply predominance of species with 1:1 and 1:2 stoichiometry (10 -14), which is contradicted by X-ray absorption data showing a trigonal CuS 3 coordination environment (15,16). Although these discrepancies might be attributed to differences in experimental conditions, the usefulness of complexation studies is contingent on their applicability to biological systems, where GSH is present in large excess compared with copper. Standard potentiometric techniques are not applicable at biologically realistic GSH/copper ratios, and the redox instability of free aquacopper(I) at typical analytical concentrations poses an additional challenge for determining reliable formation constants.
We recently reported a set of three water-soluble monovalent copper ligands, MCL-1, MCL-2, and MCL-3 (Scheme 1), that form discrete air-stable Cu(I) complexes with well-defined 1:1 metal-ligand binding stoichiometry (17). By offering complementary buffer windows stretching from 10 Ϫ10 to 10 Ϫ17 M, this ensemble of ligands is well suited to determine Cu(I)-complex dissociation constants in the pico-to sub-femtomolar range by means of ligand competition titrations. Taking advantage of these properties, we employed the monovalent copper ligand reagents as affinity standards to elucidate the nature of the Cu(I)-GSH equilibrium system at biologically relevant concentrations, where GSH is present at large excess compared with Cu(I).

Competition titration at constant pH
The three MCL-affinity standards form complexes with apparent dissociation constants ranging from 10 pM to 100 aM at pH 7 (17,18), thus covering the full range of Cu(I) affinities previously reported for GSH. Both the free ligands and their Cu(I) complexes are transparent across the visible and near-UV range, permitting spectroscopic observation of Cu(I) complexation equilibria in their presence. Monitored by UV absorption spectroscopy under anaerobic conditions (Fig. 1A, red trace), addition of 100 M aqueous Cu(II) aq to 4 mM GSH in pH 7 buffer produced a strong charge-transfer band with a shoulder near 300 nm. Consistent with quantitative reduction to Cu(I), an identical spectrum was observed after addition of [Cu(I)-(MCL-2)]PF 6 . (Fig. 1B). Back titration with the high-affinity ligand MCL-1 (17), which offers a logK of 16.0 at pH 7, yielded a steady decrease in absorbance (Fig. 1A), indicating competitive chelation of GSH-bound Cu(I). Although evolving factor analysis of these data supported only formation of a single new absorbing species, non-linear least-squares fitting using equilibrium models with either 1:1 or 1:2 Cu(I)-GSH stoichiometries resulted in trending residuals (Fig. 1C).
Consistent with the large excess of unbound GSH, changes to the GSH stoichiometry or inclusion of multiple mononuclear species had no discernible effect on the fit quality. In contrast, varying the nuclearity revealed a sharp minimum in fitting error for a tetranuclear model (Fig. 1D). This behavior could be explained by cooperative assembly of an adamantane-like {Cu 4 S 6 } cluster ( Fig. 2A), which is the most common structural motif for crystalline Cu(I)-thiolate complexes (19). With macroscopic pK a values of 1.98, 3.49, 8.64, and 9.36 (20), corresponding to sequential deprotonation of the Glu-␣-COOH, Gly-COOH, -SH, and Glu-NH 3 ϩ groups, GSH exists primarily as a monoanion at neutral pH. A species distribution diagram for the stepwise protonation of glutathione as a function of pH is provided in Fig. S1. As Cu(I)-thiolate cluster formation entails loss of the -SH proton, the corresponding pK a should be taken into account when extracting the stability constant of the cluster from titration data collected at pH 7. Employing a thiol pK a of 8.75, derived by correcting the macroscopic pK a of 8.64 from a proton concentration to a proton activity scale (17,21), a Cu 4 (GS) 6 complexation model indeed yielded small, random residuals (Fig. 1C) with a reproducible stability constant of log␤ 46 ϭ 85.0.

Low-temperature phosphorescence
Because many Cu(I)-thiolate clusters are phosphorescent, we turned to luminescence spectroscopy to further probe the nature of the multinuclear Cu(I)-glutathione complex. When employing conditions identical to above complexation studies, no emission was observed at ambient temperature; however, at 77 K a vitrified sample revealed a strong and unusually narrow phosphorescence band at 423 nm (Fig. 2B) with a luminescence lifetime of ϭ 42 s (Fig. S2). Because propylene glycol (33% v/v) was added to achieve vitrification, we also acquired the phosphorescence spectrum of a glycol-free frozen solution as control. These measurements revealed the same band, albeit slightly broadened. Furthermore, glycol addition had no significant effect on the UV spectrum at room temperature.
The distinctive signal at 423 nm can be attributed to a cluster-centered emission previously reported for complexes containing an adamantane-type {Cu 4 S 6 } core, including crystalline [Et 4 N ϩ ] 2 [Cu 4 (SPh) 6 2Ϫ ] (22) and certain copper-loaded metallothioneins (23,24), whereas no other Cu(I)-thiolate cluster arrangement is known to emit below 550 nm (22,23,25). Although the 423-nm band remains the dominant spectral feature over total Cu(I) concentrations spanning at least 3 orders of magnitude from 1 M to 1 mM (Fig. 2B), close examination of the baseline revealed a second, much weaker emission at ϳ590 nm (Fig. 2B, inset).
To determine whether the long-wavelength band is due to a change in protonation state or a different cluster stoichiometry, we explored the effect of GSH concentration and pH on the phosphorescence spectrum. As noted above, free GSH exists primarily in monoanionic form at pH 7, whereas cluster formation presumably entails S-deprotonation to the corresponding dianion, which we represent as GS Ϫ . An equilibrium between clusters differing in Cu(I)/glutathione ratio would therefore depend on [GS Ϫ ], which is influenced by both [GSH] and pH. Once the latter is significantly below the thiol pK a of 8.75, each further unit drop in pH will decrease the GS Ϫ /GSH ratio by 10-fold. We therefore measured pairs of solutions containing either 1 mM GSH or 10 mM GSH at 1 unit lower pH.
As evident from Fig. 2C, the long-wavelength band steadily increases relative to the 423-nm emission upon acidification from pH 8 to 5; however, each spectrum is essentially superimposable with its matched [GS Ϫ ] counterpart. These data unequivocally demonstrate that the long-wavelength band originates from a distinct Cu(I)-thiolate cluster with a lower GS/Cu(I) ratio than the predominant {Cu 4 S 6 } species. Each order of magnitude decrease in [GS Ϫ ] increased the abundance of the minor species by only 3-fold, or 0.5 orders of magnitude, implying that the equilibrium between {Cu 4 S 6 } and the minor complex yields only half an equivalent of GS Ϫ per eq of the latter. Of all previously observed Cu(I)-thiolate cluster stoichiometries, only {Cu 5 S 7 } is consistent with this result (see below).
Although the emission maximum of the long-wavelength band appears unchanged at the scale presented in Fig. 2C, nor-

Scheme 1
Glutathione-copper(I) equilibrium system malization of the 1 mM GSH spectra at 590 nm revealed a subtle shift from 588 nm at pH 7 and above to 594 nm at pH 6 ( Fig. 2, D and E). Increasing the excitation wavelength to 350 nm had no effect at pH 7 but further shifted the long wavelength band to 611 nm at pH 6, indicating that the onset of the shift corresponds to the appearance of a third coordination species. This was verified by decreasing the total glutathione concentration under otherwise identical conditions (100 M Cu(I), pH 5, ex 310 nm), which shifted the long-wavelength emission from 593 nm at 10 mM GSH to 605 and 617 nm at 4 and 1 mM, respectively ( Fig. 2F). At pH 5 and 1 mM GSH, a fourth emission is apparent as a shoulder on the 423-nm band but is already greatly diminished at 4 mM GSH. Precipitation occurs upon further acidifi-cation as the pH approaches 4, suggesting formation of a coordination polymer. Taking the 423 nm {Cu 4 S 6 } emission as a reference point, the band in the longest wavelength region (Ͼ630 nm) increases 21-fold in response to a 10-fold reduction in GSH concentration. This implies the presence of a third species whose net formation from Cu 4 (GS) 6 yields more than 1 molar eq of GS Ϫ . The most likely candidate is a {Cu 5 S 6 } cluster, given the previously reported long-wavelength phosphorescence at 618 nm of Et 4 N ϩ [Cu 5 (SAd) 6 Ϫ ] (Ad ϭ 1-adamantyl) (25). The symmetrical {Cu 5 S 6 } cage is the second most frequently crystallized Cu(I)-thiolate motif after {Cu 4 S 6 } and is structurally related to the {Cu 5 S 7 } by a missing bridging thiolate ( Fig. 2A).

Glutathione-copper(I) equilibrium system Combined potentiometric and spectrophotometric titrations
To fully characterize the GSH-Cu(I) solution equilibrium system, including the two additional clusters detected by phosphorescence, we performed further ligand competition titrations using pH as the independent variable. Because the thiol group deprotonates upon Cu(I) coordination, the position of the GSH-Cu(I) equilibrium system is sensitive toward pH changes. Increasing the pH in the presence of a non-protic competitor ligand should therefore produce a distinct binding isotherm, reflective of the complex stoichiometry, for each coordination species. Simulated titrations predicted that MCL-3, a non-protic Cu(I) chelator with a logK of 13.8 (17), should yield titration midpoints in the acidic range ( Fig.  3A), where all three Cu(I)-glutathione species can be detected according to the phosphorescence data. With variable MCL-3 concentrations as a second titration dimension, discrimination between species of similar Cu/GS ratio such as Cu 4 (GS) 6 and Cu 5 (GS) 7 can be further improved (Fig. 3B).
We thus performed a series of spectrophotometric pH titrations using 4 mM GSH, 100 M Cu(I), and 16 different MCL-3 concentrations ranging from 0.1 to 1.7 mM (Fig. 4). Each titration was conducted in duplicate and entailed 50 absorption traces for which the pH was varied between 4.5 and 7.5 by addi-tion of KOH in the presence of a non-coordinating universal buffer ( Fig. 4A; see "Experimental procedures" for details). An additional pair of titrations was conducted in the absence of MCL-3 to aid in spectral discrimination of the minor species (Fig. 4B). The cumulative data set, which entails a total of 900 spectra or 86,400 individual data points, was subjected to a global least-squares fitting procedure implemented in BeerOz (26). As evident from the absorption profiles at 308 nm, the previously proposed models with 1:1 or 1:2 Cu(I)-GSH stoichiometries yielded either a poor fit (Fig. 4C) or a close fit only at intermediate MCL-3 concentrations (Fig. 4D). In contrast, an equilibrium system comprising Cu 4 (GS) 6 , Cu 5 (GS) 7 , and Cu 5 (GS) 6 produces an excellent fit across the entire data set, including the titrations devoid of MCL-3 (Fig. 4E). The corresponding stability constants for each of the three species are compiled in Table 1, and detailed residual maps for the global fitting results are provided in Fig. S3. It is noteworthy that the stability constant of log␤ 46 ϭ 84.9 for Cu 4 (GS) 6 closely mirrors the value of 85.0 determined using MCL-1 under an orthogonal titration mode with a single-component fit at pH 7 (Fig. 1A).

Speciation of the Cu(I)-GSH equilibrium system
Consistent with this remarkable agreement, the fitted stability constants for all three species (Table 1) indicate a distribution dominated by Cu 4 (GS) 6 above pH 5 (Fig. 5A). Accordingly, fitting the entire titration dataset with this component alone yields an essentially unchanged log␤ 46 of 85.0, albeit with a reduced fit quality most apparent in the titrations without MCL-3 (Fig. 4F).
Because such direct, independent corroborating evidence is not available for the minor species apparent at low pH, we verified our equilibrium model by comparing the species distribution calculated from the titration-derived stability constants to the relative intensities of the phosphorescence bands over conditions selected to cover a broad range of different speciations (Fig. 5B). The close agreement between the two independent techniques supports the assignment of all major components of the Cu(I)-GSH equilibrium system within the concentration and pH ranges investigated and reaffirms the predominance of Cu 4 (GS) 6 under biologically relevant conditions. It should be noted that the observed 1:1 correspondence is expected only if the product of molar absorptivity and phosphorescence quantum yield is comparable for all three clusters. This condition is apparently met given the similar molar absorptivities at 310 nm indicated by the deconvoluted absorption spectra and the suppression of most non-radiative excited-state deactivation pathways at 77 K. The actual phosphorescence quantum yields could not be measured due to strong and variable scattering of the excitation beam by the samples.

Cluster speciation
Both the nuclearity analysis of the MCL-1 titration dataset (Fig. 1D) and the predominance of the 423-nm phosphorescence emission (Fig. 2B) suggest that the glutathione-Cu(I) solution equilibrium is dominated by formation of an adaman-

Glutathione-copper(I) equilibrium system
tane-like {Cu 4 S 6 } cluster, which remains the predominant coordination species over a Cu(I) concentration range spanning at least 3 orders of magnitude. This result is also fully consistent with previous X-ray absorption studies, which indicate a trigonal coordination environment (15,16). Each face of the {Cu 4 S 6 } cluster consists of four trigonal CuS 3 faces ( Fig. 2A), and the Cu-S and Cu-Cu bond distances of 2.26 and 2.69 Å, which were determined for the Cu(I)-GSH system at pH 8 (16) ] (27). Although the literature extended X-ray absorption fine structure data and our titration results at pH 7 can be fully accounted for by {Cu 4 S 6 } alone, the phosphorescence data reveal the presence of a minor coordination species characterized by a broad emission at 590 nm. Although several Cu(I) clusters with various stoichiometries have been reported to exhibit low-temperature phosphorescence within this spectral region, we deduced that the GSH-derived species is {Cu 5 S 7 } based on the observed dependence of the phosphorescence spectrum on the concentration of S-deprotonated glutathione (GS Ϫ , Fig. 2C). The fact that the spectrum does not respond directly to pH but rather to [GS Ϫ ] indicates both that -SH is the only protic group linked to the equilibrium between the two clusters and that GS Ϫ appears directly in the equilibrium expression. Within the experimental pH window, GSH can therefore be modeled as a mono-protic thiol RSH and the tetranuclear cluster as Cu 4 (RS) 6 2Ϫ . Thus, the relative ratio of Cu 4 (RS) 6 2Ϫ and a minor species with stoichiometry Cu x (RS) y zϪ is related to pH-dependent [GS Ϫ ] according to Equation 1, with nominal cluster charge z ϭ y Ϫ x and stoichiometric coefficients A ϭ 4/x and B ϭ 6 Ϫ Ay. The corresponding equilibrium constant for this reaction can be expressed as in Equation 2, and solving for the concentration of the minor cluster yields Equation 3, Given the large excess of glutathione over Cu(I), the concentration of [GS Ϫ ] remains essentially unaltered for a given pH regardless of the cluster stoichiometry. Furthermore, over the range of RS Ϫ concentrations where Cu 4 (RS) 6 2Ϫ remains the predominant Cu(I) species, its concentration can be approximated as constant, making the concentration of the minor species Cu x (RS) y zϪ inversely proportional to [RS Ϫ ] B/A . Thus, an

Table 1 Experimental formation constants for Cu(I)-glutathione clusters at 25°C (0.1 M ionic strength)
Data were determined based on global non-linear least-squares fitting of spectrophotometric titrations data covering the pH range from 4.5 to 7.5 (see text for details).

Glutathione-copper(I) equilibrium system
order of magnitude reduction in [RS Ϫ ] will increase the abundance of the minor species by a factor of 10 B/A . Concluding from the B/A ratios calculated for all stoichiometric combinations up to 12 Cu(I) and 12 thiolates (Table S1), only formation of species with 1:1, 3:4, 5:7, or 7:10 Cu(I)-glutathione stoichiometry could account for the observed 0.5 order of magnitude increase in relative abundance per 10-fold reduction in [RS Ϫ ]. Within these possibilities, only Cu 5 (RS) 7 2Ϫ has been experimentally observed and structurally characterized. Furthermore, {Cu 5 S 7 } cages have been synthesized under similar conditions or even crystallized from the same solution as the corresponding {Cu 4 S 6 } species, implying that the two species coexist in equilibrium (28,29).
Because the accuracy of the preceding approximate treatment is expected to diminish once the minor cluster represents a significant fraction of total copper, we confirmed the results by calculating the actual change in species ratio versus Cu 4 (GS) 6 upon lowering the pH from 7 to 6 for all crystallographically observed Cu(I)-thiolate cluster ratios where B/A is in the range of 0 -2 (Table S2), as well as hypothetical species with B/A ϭ 0.5 for comparison. Using the program HySS (30) and an activitycorrected thiol pK a of 8.75 for GSH, we adjusted the formation constant of each minor species such that its abundance is 10% relative to Cu 4 (GS) 6 (log␤ 46 ϭ 84.9) at pH 7, 100 M Cu(I), and 1 mM GSH, consistent with the relative areas of the integrated long-wavelength and short-wavelength phosphorescence bands observed under these conditions. The results in Table S2 confirm that Cu 5 (GS) 7 is the only stoichiometry that is both consistent with the phosphorescence behavior and corresponds to a crystallographically observed Cu(I)-thiolate cluster structure. Furthermore, the predicted 22-fold increase in [Cu 5 (GS) 6 ] per unit reduction in pH is consistent with the 21-fold increase in long-wavelength phosphorescence observed upon reduction of [GSH] from 10 to 1 mM at pH 5, where {Cu 5 S 6 } is expected to overtake {Cu 5 S 7 } as the most abundant pentanuclear species. The excellent correlation between the phosphorescence spectra and the calculated species distributions derived from the independent competition titration dataset (Fig. 5B) suggests that both pentanuclear species have been correctly assigned.

Role of GSH in copper homeostasis
The equilibria unraveled in this work have important implications for intracellular copper transport. Remarkably, GSH alone is sufficient to hold the thermodynamically available Cu(I) aq concentration below 1 fM under cytosolic conditions, even in a scenario where the entire cellular copper content would be available for GSH complexation (Fig. 5C). Furthermore, the cooperative assembly of a tetranuclear cluster provides a clamping effect, insulating Ϫlog[Cu(I) aq ] (pCu) 4-fold more effectively from increases in total Cu(I) than a 1:1 binding stoichiometry would imply. This is best illustrated by the simplified pCu Equation 4, which was derived from the formation constant expression of Cu 4 (GS) 6  Although our data indicate that GSH has the ability to clamp cytosolic Cu(I) aq to the sub-femtomolar range, it remains debatable whether this mechanism serves a physiological role under ordinary conditions. Because the GSH-Cu(I)-O 2 system can produce deleterious superoxide (9), it would appear reasonable that intracellular copper is regulated below the range where assembly of Cu(I)-glutathione clusters can occur. This would be consistent with the most recent affinity determinations of Cu(I) chaperones, which yielded dissociation constants in the attomolar range (31)(32)(33)(34). Perhaps the remarkable Cu(I) cluster-forming tendency of glutathione, and potentially other cysteine-derived monothiols as well, is an underlying reason why these chaperones have evolved such large Cu(I) affinities.

Implications for cellular copper stores
It is important to note that the high affinity of proteins involved in cellular copper transport does not necessarily imply a slow metal-exchange kinetics. Below the picomolar range, metal exchange may still proceed rapidly through associative  Table 1). B, linear correlation between the calculated mole fraction of Cu 4 (GS) 6

Glutathione-copper(I) equilibrium system
mechanisms without intermediacy of the free aquated Cu(I) cation. In fact, single molecule FRET studies of the copper metallochaperone Atox1 and its binding partners imply a highly dynamic exchange mechanism involving formation of multiple interaction complexes characterized by short lifetimes (35)(36)(37). Although the transfer of Cu(I) could not be observed directly based on these measurements, the short lifetime of the underlying protein-protein interaction complexes is consistent with rapid Cu(I) transfer through an associative exchange as originally proposed by O'Halloran and co-workers (38).
An associative metal-transfer mechanism is likely also responsible for shuttling Cu(I) in and out of proteins maintaining cellular copper stores. Although specific proteins that contribute to this labile pool may vary between the cytosol, different compartments, and organelles, or across various organisms, there is increasing evidence that, as a unifying theme, the underlying Cu(I)-binding sites are rich in thiols and characterized by sub-femtomolar dissociation constants. For example, Pushkar et al. (39) reported copper accumulation in the glia of the sub-ventricular zone in rodent brain. Quantified by X-ray fluorescence microscopy, the copper levels were thousands of times higher compared with cells in the surrounding tissue. Microprobe K-edge X-ray absorption spectra of copper-rich locations were consistent with reduced Cu(I) and resembled the spectrum of copper-metallothionein, suggesting that metallothionein might in fact serve as a ligand for copper storage. Given the high Cu(I) affinity of metallothionein with reported dissociation constants between 10 Ϫ19 and 10 Ϫ20 (40, 41), GSH would not be able to compete for Cu(I) binding, even when present at tens of millimolar concentrations. A similar accumulation of copper in the form of intracellular deposits was also observed in the brain tissue of Ctr2 knock-out mice (42), although the speciation of this pool was not further investigated. The utilization of dedicated high-affinity ligands for maintaining safe copper stores appears to apply also to prokaryotes as underscored by the recent discovery of a novel Cu(I) storage protein (Csp1) isolated from a methane-oxidizing bacteria requiring large quantities of copper (43). Composed of a cysteine-rich four-helix bundle, the tetrameric protein is capable of binding up to 52 Cu(I) ions with an average dissociation constant of 10 Ϫ17 .
In summary, based on the presented data we established a robust model for the GSH-Cu(I) equilibrium system under physiologically relevant conditions. Anchored by low-temperature phosphorescence studies, our model is consistent with cooperative assembly of a tetranuclear cluster that sets a firm upper limit of ϳ1 fM on the thermodynamic availability of intracellular copper. Conversely, the existence of intracellular Cu(I) pools that are buffered to the upper femtomolar or picomolar range, as implied from studies with synthetic fluorescent probes (44 -46), is precluded by normal physiological levels of GSH, both in the cytosol as well as in GSH-rich compartments and organelles. Further studies of intracellular Cu(I) fluxes thus will require analytical tools capable of detecting Cu(I) at subfemtomolar buffered concentrations. Although GSH occupies a prominent role as cellular redox buffer, the interaction of Cu(I) with less abundant low-molecular weight thiols, including cysteine, coenzyme A, or lipoic acid, might prove important toward a comprehensive understanding of cellular copper homeostasis and copper-induced stress response. Likewise, the interaction of Cu(I) with prokaryotic counterparts to GSH such as mycothiol (47) or bacillithiol (48), might have important implications for organisms that do not utilize GSH for maintaining the cellular redox balance. Although we can only speculate regarding the speciation of these systems, the formation of high-affinity tetranuclear clusters might take shape as a unifying concept, especially given the prevalence of tetranuclear clusters in inorganic Cu(I)-thiolate coordination chemistry (19).

Instrumentation
Automated spectrophotometric titrations were performed using a Varian Cary Bio50 spectrophotometer equipped with a temperature-controlled cuvette holder (Quantum Northwest), a computer-controlled syringe pump (J-KEM Scientific), and a pH meter (Orion 3 Star) equipped with a double junction glass electrode (Orion 9110DJWP). The electrode was calibrated with standard pH buffers (VWR) at 25°C (water bath). Titrations were conducted under a continuous stream of argon (passed through a 5-m filter and pre-humidified by bubbling through deionized water).
Steady-state phosphorescence spectra were acquired on a PTI fluorimeter with cold finger Dewar accessory. Interference from second-order diffraction of the scattered excitation beam was eliminated using a 400-nm long-pass filter (Thorlabs FEL0400). Spectra were corrected for the intensity of the excitation source, the spectral response of the detector, and the transmittance of the filter element. Phosphorescence decay profiles were recorded on a Horiba Fluorolog-3 instrument equipped with a 75-watt xenon flash lamp (3-s full width at half-maximum) employing the same cold finger Dewar and long-pass filter as described above.

Preparation and handling of Cu(I)-GSH solutions
All experiments were conducted under exclusion of oxygen using the Schlenk technique. To avoid contamination with UVabsorbing impurities, rubber septa were washed with two portions of acetonitrile followed by two portions of water prior to use. For molar ratio titrations and phosphorescence measurements, the appropriate amount of reduced glutathione was dissolved in 10 ml of aqueous buffer and purged with argon for 15 Glutathione-copper(I) equilibrium system min. Addition of CuSO 4 or the specified Cu(I) complex followed by pH adjustment was carried out under argon. The precise copper concentrations of the prepared stock solutions were determined by inductively coupled plasma atomic emission spectroscopy.

Low-temperature phosphorescence measurements
To account for dilution by propylene glycol, GSH-Cu(I) solutions were prepared using 1.5 times the target concentrations. Solutions with final concentrations of 1-100 M Cu(I) were obtained by means of in situ reduction of CuSO 4 supplied from 3 to 60 mM stock solutions. The initial GSH concentration was adjusted to compensate for stoichiometric oxidation to GSSG upon Cu(II) reduction. For each sample, 1.0 ml of aqueous GSH-Cu(I) solution and 500 l of deoxygenated of propylene glycol were combined, and a 200-l aliquot of the resulting solution was transferred into a quartz NMR tube and vitrified by immersion in liquid nitrogen. For variable pH experiments, a single bulk GSH-Cu(I) solution was prepared in universal buffer at pH 5 for each GSH concentration, and the pH was adjusted progressively by addition of 1 M aq. KOH. Each solution was allowed to equilibrate for at least 10 min following pH adjustment before vitrification.

Universal buffer
Full-range buffering was achieved using only non-coordinating reagents with the exception of biologically relevant chloride. The buffer consisted of 10 mM MES, 10 mM PIPPS, 5 mM DEPP dihydrochloride, and 5 mM KOH in 90 mM aq. KCl for an initial pH near 3.7 and a total chloride concentration of 100 mM. The calculated pH response to added alkali is nearly linear (R 2 ϭ 0.9986) up to a pH of 9 with a slope of 0.17 pH-units per mM of added hydroxide. Commercial-grade DEPP was purified as the dihydrochloride salt. To this end, 12 ml of aq. HCl (37%) were added to a solution of 10.1 g (70.7 mmol) of DEPP in 200 ml of ethanol, and the slurry was concentrated under reduced pressure. After taking up the thick paste in boiling ethanol (200 ml), methanol was added to the point of dissolution, and the mixture was boiled down to the onset of crystallization. The mixture was diluted with ethanol, allowed to cool, and filtered under argon. After washing with ethanol, the colorless crystalline product was dried under vacuum. Yield was 13.2 g (87%).

Spectrophotometric competition titrations at constant pH
A 1-cm path length screw-top quartz cuvette was fitted with a magnetic stir bar, sealed with a PTFE-faced silicone septum, placed in a thermostatted holder at 25°C, and purged with humidified argon for 15 min. Glutathione was dissolved in pH 7 buffer (10 mM PIPES, 0.1 M KCl) at a concentration of 4.0 mM, and the pH was adjusted to 7.00 with 1 M KOH. A 3.0-ml aliquot of the resulting solution was delivered to the sealed cuvette via syringe, and CuSO 4 was added from a stock solution (5 l, 60 mM). A solution of MCL-1 was prepared in 0.1 M KCl, adjusted to pH 7.0 with HCl, deoxygenated, and loaded into the syringe pump. Argon was introduced through a pipette tip (0.3-mm inner diameter) and exited through a 25-gauge needle to a water-filled bubbler. A total of 400 l of titrant was injected into 10-l aliquots with a 10-min equilibration time between spectral acquisitions.
The UV-absorption traces were analyzed by evolving factor analysis and non-linear least-squares fitting with the SPECFIT software package (49) to yield the corresponding apparent stability constant for pH 7.0 at 0.1 M ionic strength (KCl) and 25°C. The published protonation constant of GSH (pK a ϭ 8.64) (20) was adjusted upward by 0.11 to correct for the proton activity at 0.1 M ionic strength (17,21).

Spectrophotometric pH titrations
Titrations were conducted in a modified quartz cuvette with a 3-cm extended neck, a sidearm fitted with a rubber septum for introduction of argon and titrant, and an internally threaded top port with an O-ring compression seal holding the pH electrode. Humidified argon was introduced at the bottom of the cuvette neck through a polyethylene tube (1.5-mm outer diameter) and vented at the top of the neck to a water bubbler.
For variable pH titrations in the presence of MCL-3, stock solutions of 100 M [Cu(I)MCL-3]PF 6 were prepared in 100-ml batches using the above universal buffer. Each batch was purged with argon for 15 min and stored at 4°C under 1.3 atm of argon in a flask sealed with a rubber septum. For each set of titrations, 40 mol (12.29 mg) of solid glutathione were placed in a 25-ml round-bottom flask, which was sealed with a rubber septum, evacuated, and back-filled with argon. Using a syringe, 10.0 ml of [Cu(I)MCL-3]PF 6 stock solution were added followed by the corresponding amount of free MCL-3 supplied from a 100 mM aqueous stock solution. A 3.5-ml aliquot of the resulting working solution was transferred to the cuvette under a blanket of argon. The solution was titrated at 25°C with 4-l aliquots of 0.3 M KOH up to a total volume of 400 l with 10-min equilibration time between aliquots.
For variable pH titrations in the absence of MCL-3, a 4.1 mM glutathione solution was prepared in universal buffer and deoxygenated by bubbling with argon. A 3.5-ml aliquot was transferred by syringe to the cuvette under a blanket of argon and adjusted to approximately pH 4.5 with KOH (28 l, 1.0 M) before addition of CuSO 4 stock solution (5.83 l, 60 mM). The titration was then carried out as described above.
The proton titration data were analyzed by fitting the cumulative UV absorption traces to a single equilibrium model using the global non-linear least-squares algorithm implemented in BeerOz (26). The data set included a total of 18 independent titrations with MCL-3 concentrations varying between 0 and 1.7 mM, whereas the initial GSH and Cu(I) concentrations were kept constant at 4 mM and 100 M, respectively. Specifically, duplicate spectrophotometric pH titrations were carried out at total MCL-3 concentrations of 0, 100, 120, 150, 200, 300, 500, 1000, and 1700 M. For each titration, a total of 50 absorption traces were acquired over the range of pH 4.5 to 7.5. For the global fitting procedure, spectral data between 285 and 380 nm were utilized, corresponding to a total of 86,400 independent data points. To account for proton activity, the protonation constant of GSH (pK a ϭ 8.64) (20) was corrected upward by 0.11 pH units, whereas all other species were expressed in terms of their concentrations (17). Species distributions and free