Next Article in Journal
Monitoring Arsenic Species Content in Seaweeds Produced off the Southern Coast of Korea and Its Risk Assessment
Previous Article in Journal
Assessing the Social Sustainability Indicators in Vernacular Architecture—Application of a Green Building Assessment Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of the Interaction of Gadolinium with Several Organic Ligands and Humic Acid by Ligand Competition Using 4-(2-Pyridylazo)-Resorcinol (PAR)

by
Spencer Steinberg
*,
Vernon Hodge
and
Luis Becerra-Hernandez
Department of Chemistry, University of Nevada, Las Vegas, NV 89154, USA
*
Author to whom correspondence should be addressed.
Submission received: 22 July 2020 / Revised: 27 August 2020 / Accepted: 28 August 2020 / Published: 3 September 2020

Abstract

:
Gd3+ forms a strongly colored complex with 4-(2-pyridylazo)-resorcinol (PAR) in aqueous solutions. We characterized the Gd3+-PAR complex in order to use it as a probe of Gd3+ speciation in the presence of environmentally relevant ligands. The formation of the Gd3+-PAR complex was investigated from pH 5 to 8 in the presence of excess PAR. The absorbance of the Gd3+-PAR complex dramatically increased from pH 5 to 8 and application of the method of continuous variation indicates that the complex was primarily 1:2 Gd(PAR)2 at pH 8. Stability constants for Gd3+ with other ligands can be quantified by competitive displacement of the PAR ligand. To establish the viability of this approach, we measured the stability constants between Gd3+ and several organic acids and carbonate. Our measurements show reasonable agreement with the literature values. We used the competitive displacement approach to establish that humic acids can competitively displace PAR from the Gd(PAR)2 complex.

1. Introduction

Due to the use of Gd3+ in magnetic resonance imaging and manufacturing, large quantities of gadolinium have been released into the environment [1,2,3]. Gd3+ can be extremely toxic under some circumstances, and an investigation of its chemistry in the environment is warranted. The goal of this work was to develop a method that can probe the interaction of Gd3+ with environmentally important ligands including organic acids and humic materials. In order to effectively utilize PAR as a speciation probe, it is important to establish the stoichiometry of the complexes formed and a value for the stability constants under the pH, ionic strength, and temperature conditions utilized. In addition, it is important to gauge the impact of potentially competitive side reactions such as hydrolysis on the formation of the PAR complex. These experiments were conducted with Gd3+ at micro-molar (uM) concentrations. While these concentrations are considerably higher than the observed concentrations for Gd3+ in the environment, they are considerably lower than the concentrations used in most previous studies.
Gadolinium is a rare earth that occurs primarily in the +3 oxidation state. It has a complex environmental chemistry that includes the potential for complex formation with a variety of organic and inorganic ligands, surface sorption and hydrolysis, co-precipitation and precipitation [4]. The hydrolysis behavior of Gd3+ is important for accessing the speciation of Gd3+ in the environment. The hydrolysis reactions are shown below. While more complex hydrolysis products are often evoked to rationalize potentiometric titration results, our measurements were all performed at low Gd3+ concentrations and in the presence of excess PAR, thus we are including only the first two hydrolysis reactions in our modeling for which data are available.
Gd + 3 + H 2 O β h 2   = GdOH + 2 + H + Gd 3 + + 2 H 2 O β h 2   = Gd ( OH ) 2 + + 2 H +
The hydrolysis of Gd3+ has been studied by a number of authors [5,6,7]. For the calculations of speciation in the presence of 4-(2-pyridylazo) resorcinol (PAR) performed for this investigation, both the first and second hydrolysis constants from Garcia-Rodenas and Liberman [5] were incorporated, after these values were adjusted to 0.100 ionic strength using the Davies equation. However, including the hydrolysis reactions did not make a significant difference in our calculations because of the relatively high concentrations of PAR.
PAR is a colorimetric complexing agent that has been utilized as a titration indicator and colorimetric reagent for various metals. The stoichiometry of the metal–PAR complexes has been reported as 1:1 or 1:2 by various authors. Iwamoto (1961) [8] investigated the acid–base properties of PAR and chelate formation with a several transition metals and lead. Iwamoto [8] indicated Fe3+, Co2+, and Ni2+ form PAR complexes in 1:2 at pH 4 while Cu2+, Zn2+ Cd, Hg2+, Pb2+ form 1:1 complexes at pH 10. Kocyła et al. used PAR as a metallochromic indicator in a competitive ligand (spectrophotometry) approach to measure the affinity of Zn2+ in Zn protein finger domains [9]. Kocyła et al. [9] indicated that several different protonated species of the PAR complexes may exist in solution. Under the conditions utilized in Kocłya et al. [9], the Zn2+ complex appeared to be primarily 1:2. Kocyła et al. [9] emphasized the conditional aspect of all of these measurements. PAR is not a simple ligand and several species of varying spectral properties may be present. Therefore, the observed extinction coefficient is highly pH dependent. Holm and Smothers (2006) [10] characterized the lead-complexing properties of polyphosphate water treatment products using competing-ligand spectrophotometry with PAR. These experiments employed high Pb2+: PAR ratios, and the authors indicated that 1:1 Pb2+: PAR complexes were forming under these conditions, with the possibility of multiple adsorbing species present in solution. They considered the possibility of mixed complexes between PAR and the polyphosphates but concluded that these complexes were not necessary to explain their data. Degnall et al. (1965) [11] developed an analytical method for the determination of lead with PAR. This method incorporated an elaborate solvent extraction procedure to remove interfering metal ions before forming the PAR complex. Degnall et al. [11] indicated that a 1:1 Pb2+-PAR complex was formed. Neas and Guyon (1972) [12] developed a method for the indirect determination of oxalate using the competitive disruption of the highly colored uranium and PAR complex by oxalate at pH 4.8. Florence and Farrar (1963) [13] used PAR for a spectrophotometric determination of uranium. Corsini et al. (1962) [14] conducted a potentiometric investigation of PAR with Mn2+, Co2+, Ni2+, and Zn2+ in 50% dioxane–water where their results indicated that 1:2 metal: PAR complexes were formed. Pollard et al. (1959) [15] suggested PAR as a possible analytical reagent for the colorimetric estimation of cobalt, lead, and uranium. These authors demonstrated excellent sensitivity for these metals and recommended the use of the reagent after separation using methods such as paper chromatography. Ghasemi et al. (2007) [16] studied complex formation between PAR and Al3+, Fe3+, Zn2+, and Cd2+ in an aqueous solution at 0.1 M ionic strength. McCall and Fierke (2000) [17] utilized PAR to investigate the binding of several transition metals to the enzyme carbonic anhydrase and measured the metal released from carbonic anhydrase that was previously incubated with various metal ion buffers. Their method incorporated 4 M guanidine to denature the carbonic anhydrase and release the metals into solution where they would be complexed by PAR. Similarly, Hunt et al. (1985) [18] used PAR to measure zinc release from Escherichia coli aspartate transcarbamoylase. PAR was used to measure the instantaneous concentration of Zn2+ at pH 7.0 after the aspartate transcarbamoylase was treated with the organo-mercury reagent (PMPS).
Mori et al. (1988) [19] developed an assay for mercury(II) using PAR. In addition, Mori et al. used the disruption of the Hg2+–PAR complex by cyanide in order to quantify cyanide in solution. Munshi and Dey (1964, 1971) [20,21] reported the spectrophotometric determination of lanthanides using PAR. They established 1:2 metal: PAR stoichiometry by several methods. Sommer and Novotná (1967) [22] investigated the complexation of aluminum, yttrium, lanthanum, and lanthanides with PAR and their results indicate that Al3+ forms 1:1 complexes at low pH while yttrium forms a 1:2 complex.
Thompson and Byrne (1988) [23] used PAR as an indicator complex in the study of europium carbonate equilibrium at an ionic strength of ~0.7 (similar to seawater). They used a large excess of PAR and observed competitive complexation by both the carbonate and oxalate. They assumed that the Gd3+–oxalate constant was well established and were able to estimate the Gd3+–carbonate equilibrium relative to the Gd3+–oxalate equilibrium. The mathematical development presented in the Thompson and Byrne [23] manuscript implies that both mono- and bis-PAR complexes exist.
Furthermore their data analysis including the possibility of mixed complexes, however their results indicated that they were quantitatively unimportant in their contribution to absorbance.
The pKa values for PAR have been reported by a number of groups [8,9,24]. PAR when fully protonated has a charge of +1. At circumneutral pH conditions, mono-protonated HPAR- is the dominant species. Munshi and Dey (1971) [21] measured stability constants for PAR chelates of the lanthanides and reported that the lanthanides formed 1:2 complexes with PAR. Munshi and Dey (1971) [21] reported a stability constant with Gd3+ of 109.9 for the formation of the 1:2 complex. Ohyoshi (1984) [25] reported a stability constant for the formation of Gd(HPAR)2+ and Gd(PAR)+ of 104.28 and 1010.25, respectively. Ohyoshi [25] significantly implied a 1:1 stoichiometry for the complex as these measurements were conducted with metal ions far in excess over PAR. The differences in experimental conditions makes these two results difficult to reconcile.
It is clear that PAR is a complex ligand that reacts with metal ions to produce a number of light absorbing species that are simultaneously present in solution. The speciation of PAR, and most likely the chelates, are highly concentration, pH, and ionic strength dependent. To utilize PAR in a competitive ligand study, it is important that pH, ionic strength, and concentration of PAR be well controlled during the measurements as stability constants and complex absorptivity are clearly conditional in nature. In addition, it is important that the stoichiometry of the complex be well established in order to successfully model competitive ligand interactions.

2. Materials and Methods

2.1. Materials

Water was prepared by redistilling house deionized water using an all glass still. All reagents used were of American Chemical Society (ACS) grade or better and employed without further purification. Gadolinium standards were prepared by dilution of 1000 ppm atomic absorption standards purchased from VWR Scientific (Visalia, CA). Two commercial humic acids were investigated. A leonardite humic acid was obtained from TeraVita (Lancaster, PA. This commercially-available humic material is marketed as a soil amendment (TeraVita product information). A commercial humic acid was obtained from Fluka Chemicals. A second leonardite humic acid and Pahokee peat humic acid were obtained from the International Humic Substances Society (IHSS). HEPES (4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid), MES (2-(N-morpholino)-ethanesulfonic acid), and TEA triethanolamine hydrochloride were purchased from Sigma Aldrich. Acetic acid was obtained from Baker. Buffer concentrations were 0.1 M in all experiments. The buffer pH values were adjusted with 6 M HCl or 10 M NaOH purchased from Baker. The pH of the solution was monitored and verified to within 0.02 pH units of the prescribed pH with a pH meter calibrated with National Institute of Standards and Technology (NIST) standard buffers.
Low molecular weight organic acids and sodium bicarbonate were of analytical quality and used without any further purification. PAR (4-(2-pyidylAzo) resorcinol) was purchased from Sigma Aldrich. PAR solutions of 1.20 mM were prepared gravimetrically in the various 0.10 M buffers.
An Ocean Optics SD2000 fiber optic spectrophotometer with a DT-Mini 26s deuterium tungsten light source was used to perform all measurements. A 1 cm path-length cuvette, containing a Teflon coated stir bar, was used. The cuvette holder for the spectrophotometer sits on a magnetic stirrer that operated continuously during all measurements. The temperature of the cuvette was maintained at 25 °C with a circulating water bath (VWR Scientific, Radnor, PA, USA). The instrument response was confirmed to be linear to an absorbance value of 1.8.

2.2. Methods

Measurements were generally performed as follows. A 2.4 mL volume of 0.100 M buffer was added to the cuvette followed by 0.600 mL of the PAR solution. The cuvette was allowed to equilibrate at 25 °C for 10 min and then Gd3+ (6.35 mM) was added in 2.0 uL aliquots with a micropipette. The solution was allowed to equilibrate for 5 min in between Gd3+ additions, although a steady reading was generally achieved in less than one minute. For ligand competition experiments, a concentrated solution of the ligand (varied with ligand) was prepared in buffer. The ligand was then added to the cuvette (21.1 uM in Gd3+) in microliter volumes. Concentrations of the ligands were adjusted to limit the volume addition to less than 25 uL in order to minimize dilution. Again, the solution was allowed to equilibrate for 5 min between each ligand addition. For the tannic acid and humic acid experiments, a stock solution containing ~100 mg/L of the polyelectrolyte was prepared in 0.100 M buffer. These stock solutions were then diluted with 0.100 M buffer to obtain various concentrations in solution. A 2.400 mL volume of the humic (or tannic) acid solution was pipetted into the cuvette and mixed with 0.600 mL of the 1.2 mM PAR. After temperature equilibration, the mixture was titrated with the Gd3+ stock solution (0 to 10 uL).

3. Results and Discussion

3.1. Gd3+-PAR Complexation Color Development as a Function of pH

For this study, most of the measurements were conducted at pH 8, where the response was favorable and hydrolysis was not extensive due to the presence of PAR concentrations that were in significant excess over the total Gd3+ concentrations. The objective was that PAR concentrations should be sufficient to assure that Gd(PAR)2 formation was essentially quantitative. The linearity of the absorbance at 510 nm with total Gd3+ concentration indicates that these conditions were achieved. With 0.6 mL of the 1.2 mM PAR reagent added to the cuvette, the resulting concentration of PAR was 0.240 mM. As a result, the PAR: Gd3+ ratio was greater than 10 for the majority of these experiments.
The maximum absorbance occurred at 510 nm for the Gd3+-PAR complex at all pH values investigated. Calibration plots for four buffers are shown in Figure 1, and demonstrate that the response at 510 nm increased linearly with Gd3+ concentration in all four buffers. In addition, we found that the Gd3+ response was significantly higher using the pH 8 TEA buffer than in the pH 8 HEPES buffer.

3.2. Gd3+: PAR Complex Ratio Investigation

In order to establish the stoichiometry of the Gd3+–PAR complex, the method of continuous variation was performed in 0.100 M TEA buffer adjusted to pH 8. A constant volume of 3.00 mL was maintained in the cuvette during these measurements. The total PAR plus Gd3+ concentration was maintained at 1.00 × 10−4 M, while the mole ratio of Gd3+: PAR was varied from 0:1 to 1:0. The results are shown in Figure 2.
These results are clear evidence that the 1:2 complex predominates under our experimental conditions. Experiments with the HEPES buffer also support the 1:2 complex. As noted above, the response (absorbance intensities) with the TEA buffer is higher than with the HEPES buffer. Although HEPES is generally considered to have a negligible interaction with metal cations, while TEA buffer is known to form metal complexes, it is clear that in this application, the formation of the PAR complex with Gd3+ was evidently more favorable in TEA [26]. There is no clear explanation for this, other than possible reduced interference of hydrolysis and adsorption in TEA. As a result of these observations, most of our experiments were conducted using the TEA buffer.
In order to determine the absorptivity of the Gd(PAR)2 complex, the following experiment was performed. A stock solution of 1.2 mM PAR in 100 mM TEA buffer (pH 8) was prepared. A series of calibration experiments was performed at pH 8 (TEA buffer) using variable PAR concentrations. The response factor for these calibrations was found to be constant with a PAR concentration above 0.200 mM. Therefore a PAR concentration of 0.240 mM was selected for these studies. With a PAR concentration of 0.240 mM, the above experiment indicates that Gd3+ is essentially quantitatively converted to the Gd(PAR)2 complex and the slope of the absorbance at 510 nm against gadolinium(III) concentrations was taken as the extinction coefficient for this complex. An average molar absorptivity (extinction coefficient) of 6.54(±0.15) × 104 M−1 was measured.
The response factor at 0.240 mM PAR was re-measured before each experiment and did vary ±5%. The variation may be instrumental in origin or due to small variations in the buffer pH or temperature. Using this extinction coefficient estimate, the MCV results were utilized to estimate the stability constant for the PAR complex by optimizing the fit of a simple mathematical model to these measurements. The absorbance at 510 nm was converted to complex concentration and the concentration of unbound PAR was calculated by mass balance. These measurements were fit to a speciation model as outlined below. In this exercise, only the first and second hydrolysis of Gd3+ were included in the analysis [5]. The hydrolysis constants have been estimated in previous studies using potentiometric methods and estimates for the second hydrolysis, constant, are not as well established. In order to apply potentiometric methods to determine the hydrolysis constants, the concentration of Gd3+ has to be at least on the order of mM. At these concentrations, the precipitation of Gd3+ interferes with the determination of the second and third hydrolysis constant. In addition, formation of dimeric and trimeric Gd3+ species is also a possibility [7,27]. Under the high PAR concentration conditions used in this investigation, hydrolysis made a very small contribution to Gd3+ speciation. We did not see any evidence of precipitation during our measurements.
Our simple model that considers “free Gd3+, the first and second hydrolysis product and the PAR complex, to be the only significant Gd3+ species in solution. The MCV results were fit to Equation (1). In Equation (1), Gdt is the total Gd3+ concentration and βPAR is the stability constant for the formation of Gd(PAR)2. The curve was fit using CurveExpert Professional Software 2.6.5 (Hyams Development). Similar results were obtained with using the Excel Solver to minimize the difference in the sum of square differences between the concentration of Gd(PAR)2 and the calculated concentration of Gd(PAR)2.
[ Gd ( PAR ) 2 ] = β PAR Gd t [ PAR ] 2 1 + β h 1 [ H + ] + β h 2 [ H + ] 2 + β PAR [ PAR ] 2
The MCV data can be fit to Equation (1) by optimizing βpar. A best fit value of log βPAR = 9.51 ± 0.06 was determined by this method. This value was utilized to assign stability constants to other ligands. We are cognizant that this estimate for log βpar likely only applies to the conditions of this experiment. We have not made any attempt to determine the protonation state of this complex. However, our results did demonstrate that absorbance increased strongly with the increase in pH.

3.3. Gadolinium Interaction with Various Ligands

In order to demonstrate the feasibility of this experimental approach to estimate the stability constant of Gd3+ with various ligands, the following procedure was followed. The determination of the stability constant for Gd3+ with oxalate will serve as an example. A 21.1 uM solution of Gd3+ in 0.1 M TEA (pH 8) and 0.240 mM PAR was prepared. This was done with 2400 uL of 100 mM TEA buffer and 600 uL of 1.2 mM PAR that was dissolved in 100 mM TEA (total volume 3.00 mL). Using the value of βPAR for the formation of Gd(PAR)2, a “reasonable” preliminary assumption is that all of the Gd3+ in this solution is essentially Gd(PAR)2. A 0.200 M solution of sodium oxalate was prepared and added to the solution (2–5 uL at a time). The objective was to minimize any dilution during the addition of oxalate as a competitive ligand. The results show that oxalate (at the concentrations introduced) can compete for Gd3+.
In order to estimate the stability constant of Gd+3 and oxalate, it was assumed that in this solution, Gd3+ exists as the Gd(PAR)2 complex or as an oxalate complex. The hydrolysis of Gd3+ was ignored in this preliminary calculation. It was assumed that free Gd3+ and Gd3+ hydrolysis products were of minor importance. It was also assumed that GdOx+ was the dominant oxalate complex. If this assumption is valid, the concentration of Gd(PAR)2 is given by the following equation (to a good approximation).
[ Gd ( PAR ) 2 ] = β PAR Gd t [ PAR ] 2 β OX [ O x 2 ] + β PAR [ PAR ] 2
Gdt is the total amount of gadolinium in solution. [Gd(PAR)2] is estimated from the absorbance at 510 nm. [PAR] in this equation is estimated from the mass balance and is given by subtracting the Gd(PAR)2 concentration from the total concentration of PAR in the solution. As the βox << βPAR for most of the plot [Ox2−] is equal to the total concentration of oxalate in the solution. This equation can be rearranged to:
Gd t [ Gd ( PAR ) 2 ] = 1 + β OX [ O x 2 ] β PAR [ PAR ] 2
The data from the above experiment were plotted according to the above equation in Figure 3. The fact that the plot is linear indicates that the previous speciation assumptions were reasonable and the slope of this plot was used to determine a preliminary estimate for βox. This approach was utilized for obtaining preliminary estimates for all of the low molecular weight ligands that we examined.
From the slope of the plot, βoxal was estimated as 105.02. This estimated value was refined by fitting the data to Equation (4) by using the software package CurveExpert Professional.
[ Gd ( PAR ) 2 ] = β PAR Gd t [ PAR ] 2 1 + β h 1 [ H + ] + β h 2 [ H + ] 2 + β PAR [ PAR ] 2 + β OX [ O x 2 ]
The decrease in the concentration of Gd(PAR)2 resulting from the addition of oxalate to the solution is illustrated in the Figure 4. Including a bis–oxalato complex did not improve the fit. If Gd(Ox)2− exists, it is clearly a minor species under these conditions. The value for the GdOx+ stability constant (log βoxal = 4.96) obtained is of course conditional. This estimate procedure was simplified by the rather low pKa values for oxalate, enabling us to consider all of the oxalate completely ionized at pH 8. This same approach could be utilized for the other low molecular weight organic acids examined in this study. We cannot rule out the occurrence of mixed complexes between oxalate and PAR. However, it is clear that our data reasonably fits this simplified model.
The results for other organic acids examined by this method are shown in Table 1. The stability constants from Martel and Smith (1977) and Gustafsson [28,29] are shown for comparison. The carbonate stability constant was extrapolated from the stability constant incorporated into the Visual MINTEQ software [29]. Extrapolation of the zero ionic strength value from Lee and Byrne (1993) [30] resulted in a log K of 6.63. All of the results from this study were of the same order of magnitude and were well correlated (R = 0.94) with the literature results. Differences are likely a function of the experimental methods utilized.
Wood (1993) [31] tabulated critical stability constants for rare earth element complexes with simple carboxylic acids at 25 °C and 1 bar. He reported a log β stability constant for oxalate and malonate of 6.68 (estimated by correlation) and 5.87, respectively at infinite dilution. Using the Davies equation to correct these values to 0.1 M ionic strength yielded values for oxalate and malonate of 5.39 and 4.59, respectively. The oxalate number was ~0.1 log units lower than this report, while the malonate value was considerably higher than our measurement based on competition with PAR.
The overall reasonable correlation between our results and previous measurements demonstrates that the competitive ligand approach can offer valuable insight into Gd3+ speciation in an aqueous solution in the presence of low molecular weight ligands. We next report the application of this method quantifying Gd3+ interaction with polyelectrolytes of environmental interest.

3.4. Tannic Acid

We tested tannic acid as a potential ligand for Gd3+ and as a humic acid analog. Tannic acid is a biopolymer comprised of polygallate esters with a carbohydrate core. It is a natural product that probably exists in soils and is likely incorporated into the humic fraction. The molecular weight of the material we utilized was reported by the manufacturer as 1701.2 g/mol.
Metal interaction with tannic acid should be primarily through the phenolic hydroxyl groups. Phenolic hydroxyl groups have been hypothesized by many workers to be important functionalities for the interaction of humic acids with various metals. We examined the influence of tannic acid at various concentrations on the Gd(PAR)2 response. The solutions were prepared by mixing different proportions of a tannic acid stock solution (in TEA) with TEA pH 8 buffer. The total volume of these mixtures again amounted to 2400 µL. A 600 µL aliquot of 1.2 mM PAR in 0.1 M TEA pH 8 buffer was added and the mixture, which was placed in the cuvette. Gd3+ was then titrated into the solution. The results for different concentrations of tannic acid at pH 8 are presented in Figure 5.
For these experiments, we titrated considerably more Gd3+ into solution than in some of the other experiments reported. The shape of the titration curves were fairly consistent and appear to indicate a finite number of binding sites for Gd3+ on the tannic acid molecules. Note that the slope of the curve (response) increased as additional Gd3+ was added and as the complexing capacity of the dissolved tannic acid was apparently being exceeded. From the molar concentration of tannic acid, it is clear that more than one Gd3+ binds to each tannic acid molecule. This observation is an important consideration for modeling the interaction of Gd3+ with polyelectrolytes.

3.5. Gadolinium Interaction with Humic Acids

3.5.1. International Humic Acid Society (IHSS) Leonardite HA

Leonardites are oxidation products (weathering products) of lignite coals. The humic acids are extracted from the leonardite with alkali. A leonardite humic acid was purchased from the IHSS. A 116 mg/L humic acid stock solution were prepared in 0.1 M TEA pH 8 buffer. Various portions of TEA pH 8 and the humic stock solution were mixed, keeping the total volume 2400 μL. Then, 600 μL of 1.2 mM PAR in TEA pH 8 was added. The solution was titrated with Gd3+. The results of this titration are presented in Figure 6.
Clearly, the humic acid is capable of sequestering Gd3+, rendering it unavailable to the PAR. However, it is not a simple task to assign a stability constant to the Gd3+ humic complex. IHSS does provide estimates of carboxylic acid and phenolic group content for samples. However, assuming a molecular weight for the humic acid of ~1000 g/mol, it is apparent that like tannic acid, it is likely that more than one Gd+3 is bound per molecule of humic acid.

3.5.2. TeraVita LeonarditeIHSS

In our previous study (Steinberg and Hodge, 2018) [32], we demonstrated the ability of TeraVita humic acid to complex copper as Cu2+. In this study, experiments with this humic acid were done in TEA pH 8 buffer with a humic acid stock solution of 122 mg/L. Various portions of TEA pH 8 and the humic stock solution were mixed, keeping the total volume 2400 μL. Then, 600 μL of 1.2 mM PAR in TEA pH 8 was added. The solution was titrated with Gd3+.
The addition of Gd3+ was done at four different humic acid concentrations and the results are presented in Figure 7. The suppression of the formation of Gd(PAR)2 is apparent and increases with increasing humic acid concentration.

3.6. IHSS Pahokee Peat Humic Acid

The ability of the IHSS Pahokee peat humic acid to sequester gadolinium at pH 8 was tested. These experiments were done in TEA pH 8 buffer using 600 µL of 1.2 mM PAR. The data for the Pahokee peat humic acid is reported in Figure 8. As the concentration of humic acid increases the Gd3+ response is increasingly suppressed demonstrating the competitive complexation of Gd3+ by the Pahokee peat humic acid.

3.7. Fluka Humic Acid

The ability of the Fluka humic acid to sequester gadolinium was tested in TEA pH 8 buffer. The humic acid was prepared in 0.1 M TEA pH 8 buffer at a concentration 110 mg/L. Various portions of TEA pH 8 and the humic stock solution were mixed, keeping the total volume 2400 μL. Then, 600 μL of 1.2 mM PAR in TEA pH 8 was added. The solution was titrated with Gd3+. The Gd3+ stock solution was a commercial 1000 mg/L solution. The data for all of the various concentrations tested are summarized in Figure 9.

3.8. Estimated Stability Constants for Gd3+ and Humic Acids

In the cases of the Pahokee peat humic acid and the IHSS Leonardite humic acid, the IHSS has provided measurements of the carboxylic and phenolic functional group contents. For the Pahokee peat humic acid, the carboxylic acid and phenolic acid groups occurred at 9.01 and 1.91 meq/g, respectively. The Leonardite humic acid is reported to contain 7.46 and 2.31 meq/g for carboxylic acids and phenolic acids, respectively. As these functionalities are generally thought to be involved in humic metal binding, we examined if competition between PAR and these two humic substances could be related to the concentrations of these two functional groups. Using the Gd3+ binding data at 21 uM total Gd3+ and variable concentrations of PAR, we adapted Equation (3) to analyze the data and found that there was a correlation between Gdt/Gd(PAR)2 and the ratio of the square of the sum of the total carboxylic acid and phenolic functional group content to PAR2 ratio. This correlation yielded r2 values of 0.993 and 0.996 for the Pahokee peat and Leonardite humic acids. We assumed for the Pahokee peat and the IHSS Leonardite results that Gd(PAR)2 and a Gd3+ humic acid complex were the predominant solution species. The calculated log β values for the interaction of the Gd3+ and humic acid were 9.5 and 9.4 for the Pahokee peat and the IHSS Leonardite, respectively. These calculated constants were very close to the conditional stability constant for PAR. It is tempting to interpret this simple correlation as an indication of bidentate interaction with the humic acid. However, this correlation does not consider the heterogeneous nature of the humic acid or the potential contribution of electrostatic potentials to metal humic interaction, so this is at best a very preliminary result. We believe that our observations can be incorporated a humic acid-metal binding model that considers heterogeneity and electrostatic effects. If this can be accomplished, it should lead to results that may have a predictive value.

4. Conclusions

We have demonstrated that through a competitive ligand methodology using PAR, it is possible to measure stability constants for various ligands of potential environmental importance. The sensitivity of the method is sufficiently high so that low (uM) concentrations of Gd3+ can be utilized. With excess PAR present in solution, the complications of hydrolysis are largely avoided. What is clear from this data is that various organic acids and humic acids can influence the speciation of Gd3+ under circumneutral pH conditions. On a mass basis, the Gd3+ binding to all of the humic materials was very similar and we anticipate that similar stability constants would be measured.
A large fraction of the anthropogenic Gd3+ released into the environment is likely from medical applications and is initially bound by a chelating agent. Apparently, these chelating agents, which are meant to make Gd3+ biologically inert during MRI testing, survive wastewater treatment [33]. Port et al. (2008) reviewed the thermodynamic and kinetic stability of various gadolinium chelates [34]. It is not clear how long these complexes will survive long term exposure to environmental chemical, photo-chemical, and biological degradation.
It has been demonstrated that Gd3+ may be partially displaced from the chelates by competition with other metals such as Ca2+ and Mg2+ when waste streams are mixed with seawater [35]. In addition, there is some evidence that zinc, calcium, and iron may be able to partially displace Gd3+ from medical chelating agents in vivo [36] despite the very high stability constants. If this occurs, then the environmental chemistry of Gd3+ will be determined by the abundance of alternative inorganic and organic ligands. These ligands including humic substances will then determine the bio-availability and toxicity of Gd3+.
This study demonstrated that various humic materials are capable of disrupting the Gd3+-PAR complex. Our estimates of stability of the stability constants for Gd3+ and humic acids are very preliminary. In our future work, will plan to improve the quantitative interpretation of these results using a humic acid metal binding model. We are presently exploring the incorporation of our results using the software available for one of the currently available humic acid metal binding models.

Author Contributions

S.S.: Conceptualization, methodology, experimental measurements, writing—original draft preparation. V.H.: Writing—review and editing. L.B.-H.: experimental measurements. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rogowska, J.; Olkowska, E.; Ratajczyk, W.; Wolska, L. Gadolinium as a New Emerging Contaminant of Aquatic Environments. Environ. Toxicol. Chem. 2018, 37, 1523–1534. [Google Scholar] [CrossRef]
  2. Ebrahimi, P.; Barbieri, M. Gadolinium as an Emerging Microcontaminant in Water Resources: Threats and Opportunities. Geosciences 2019, 9, 93. [Google Scholar] [CrossRef] [Green Version]
  3. Hatje, V.; Bruland, K.; Flegal, A. Increases in Anthropogenic Gadolinium Anomalies and Rare Earth Element Concentrations in San Francisco Bay over a 20 Year Record. Environ. Sci. Technol. 2016, 50, 4159–4168. [Google Scholar] [CrossRef]
  4. Migaszewski, Z.; Gałuszka, A. The Characteristics, Occurrence, and Geochemical Behavior of Rare Earth Elements in the Environment: A Review. Crit. Rev. Environ. Sci. Technol. 2015, 45, 429–471. [Google Scholar] [CrossRef]
  5. Rodenas, L.G.; Lieberman, S.J. Hydrolysis of Gadolinium (III) in light and heavy water. Talanta 1991, 38, 313–318. [Google Scholar] [CrossRef]
  6. Klungness, G.D.; Byrne, R.H. Comparative hydrolysis behavior of the rare earths and yttrium: The influence of temperature and ionic strength. Polyhedron 2000, 19, 99–107. [Google Scholar] [CrossRef]
  7. Yakubovich, Y.Y.; Alekseev, V.G. Hydrolysis Constants of Tervalent Lanthanum and Lanthanide Ions in 0.1M KNO3 Solution. Russ. J. Inorg. Chem. 2012, 57, 911–915. [Google Scholar] [CrossRef]
  8. Iwamoto, T. Acid-Base Property and Metal Chelate Formation of 4-(2-Pyridylazo)-Resorcinol. Bull. Chem. Soc. Jpn. 1961, 34, 605–610. [Google Scholar] [CrossRef] [Green Version]
  9. Kocyła, A.; Pomorski, A.; Krężel, A. Molar absorption coefficients and stability constants of metal complexes of 4-(2-pyridylazo)-resorcinol (PAR): Revisiting common chelating probe for the study of metalloproteins. J. Inorg. Biochem. 2015, 152, 82–92. [Google Scholar] [CrossRef]
  10. Holm, T.R.; Smothers, S.H. Characterizing the Lead-Complexing Properties of Polyphosphate Water Treatment Products by Competing-Ligand Spectrophotometry Using 4-(2-Pyridylazo)-Resorcinol. Intern. J. Environ. Anal. Chem. 2006, 41, 71–82. [Google Scholar] [CrossRef]
  11. Degnall, R.M.; West, T.S.; Young, P. Determination of Lead with 4-(2-Pyridylazo)-Resorcinol and solvent extraction. Talanta 1965, 12, 583–588. [Google Scholar] [CrossRef]
  12. Neas, R.E.; Guyon, J.C. Indirect Spectrophotometric Determination of Oxalate Using Uranium and 4-(2-pyridylzao)-resorcinol. Anal. Chem. 1972, 44, 799–805. [Google Scholar] [CrossRef]
  13. Florence, T.M.; Farrar, Y. Spectrophotometric Determination of Uranium with 4-(2-Pyridylazo)-resorcinol. Anal. Chem. 1963, 35, 1613–1616. [Google Scholar] [CrossRef]
  14. Corsini, A.; Yih, I.M.-L.; Fernand, Q.; Freiser, H. Potentiometric Investigation of the Metal Complexes of 1-(2-pyridylazo)-2-napthol and 4-(pyridylazo)-resorcinol. Anal. Chem. 1962, 34, 1090–1093. [Google Scholar] [CrossRef]
  15. Pollard, F.H.; Hanson, P.; Geary, W.J. 4-(2-Pyridylazo)-Resorcinol as a possible analytical reagent for the colorimetric estimation of cobalt, lead and uranium. Anal. Chim. Acta 1959, 20, 26–31. [Google Scholar] [CrossRef]
  16. Ghasemi, J.; Peyman, H.; Meloun, M. Study of Complex Formation between 4-(Pyridylazo) Resorcinol and Al3+, Fe3+, Zn2+, Cd2+ Ions in an Aqueous Solution at 0.1M Ionic Strength. J. Chem. Eng. Data 2007, 52, 1171–1178. [Google Scholar] [CrossRef]
  17. McCall, K.A.; Fierke, C.A. Colorimetric and Fluorimetric Assays to Quantitate Micromolar Concentrations of Transition Metals. Anal. Biochem. 2000, 284, 307–315. [Google Scholar] [CrossRef]
  18. Hunt, J.B.; Neece, S.H.; Ginsburg, A. The Use of 4-(2-Pyridylazo) resorcinol in Studies of Zinc Release from Escherichia coli Aspartate Transcarbamoylase. Anal. Biochem. 1985, 146, 150–157. [Google Scholar] [CrossRef]
  19. Mori, I.; Fujita, Y.; Fujita, K.; Nakahashi, Y.; Kato, K.; Nakamura, T. Color reaction between 4-(2-pyridylazo)resorcinol and mercury(II) in the presence of surfactant, and improved spectrophotometric determination of mercury (II) and cyanide with its color reaction. Anal. Let. 1988, 21, 2359–2379. [Google Scholar] [CrossRef]
  20. Munshi, K.N.; Dey, A.K. Spectrophotometric determination of Rare earth metals with 4-(2-Pyridylazo) Resorcinol. Anal. Chem. 1964, 36, 2003–2004. [Google Scholar] [CrossRef]
  21. Munshi, K.N.; Dey, A.K. Spectrophotometric Determination of Lanthanides Using 4-(2-Pyridylazo)-Resorcinol. Microchim. Acta 1971, 59, 751–756. [Google Scholar] [CrossRef]
  22. Sommer, L.; Novotná, H. Complexation of Aluminum, Yttrium, Lanthanum and Lanthanides with 4-(2-Pyridylazo)-Resorcinol (PAR). Talanta 1967, 14, 457–471. [Google Scholar] [CrossRef]
  23. Thompson, S.W.; Byrne, R.H. Indicator Ligands in Metal Complexation Studies: Role of 4-(2-Pyriylazo)-resorcinol in Europium Carbonate Equilibrium. Anal. Chem. 1988, 60, 19–22. [Google Scholar] [CrossRef]
  24. Schwarzenbach, G.; Biedermann, W.; Komplexone, X. Erdalkalikomplexe von o, o’-Dioxyazofarbstoffen. Helv. Chim. Acta 1948, 31, 678–687. [Google Scholar] [CrossRef] [PubMed]
  25. Ohyoshi, E. Spectrophotometric determination of formation constants of 1:1 complexes of Lanthanides with 4-(2-pyridylazo) resorcinol (PAR). Talanta 1984, 31, 1129–1132. [Google Scholar] [CrossRef]
  26. Ferreira, C.M.H.; Pinto, I.S.S.; Soares, E.V.; Soares, H.M.V.M. (Un)suitability of the use of pH buffers in biological, biochemical and environmental studies and their interaction with metal ions—A review. RSC Adv. 2015, 5, 30989–31003. [Google Scholar] [CrossRef] [Green Version]
  27. Djurdjević, P.; Jelic, R.; Joksović, L.; Lazarević, I.; Jelikić-Stankov, M. Study of Solution Equilibria between Gadolinium (III) Ion and Moxifloxacin. Acta Chim. Slov. 2010, 57, 386–397. [Google Scholar]
  28. Martell, A.E.; Smith, R.M. Critical Stability Constants Vol 3: Other Ligands; Springer Science + Business Media LLC: New York, NY, USA, 1977; pp. 329–344. ISBN 978-1-4757-1568-2. [Google Scholar]
  29. Gustafsson, J.P. Visual MINTEQ. Available online: https://vminteq.lwr.kth.se (accessed on 1 April 2019).
  30. Lee, J.H.; Byrne, R.H. Complexation of trivalent rare earth elements (Ce, Eu, Gd, Tb, Yb) by carbonate ions. Geochim. Cosmochim. Acta 1993, 57, 295–302. [Google Scholar]
  31. Wood, S.A. The aqueous geochemistry of rare-earth elements: Critical stability constants for complexes with simple carboxylic acids at 25 °C and 1 bar and their application to nuclear waste management. Eng. Geol. 1993, 34, 229–259. [Google Scholar] [CrossRef]
  32. Steinberg, S.; Hodge, V. Copper Complexation by Dissolved Organic Matter in arid Soils: A Voltammetric Study. Environments 2018, 5, 125. [Google Scholar] [CrossRef] [Green Version]
  33. Verplank, P.L.; Taylor, H.E.; Nordstron, D.K.; Barber, L.B. Aqueous Stability of Gadolinium in Surface Waters Receiving Sewage Treatment Plant Effluent, Boulder Creek, Colorado. Environ. Sci. Technol. 2005, 39, 6923–6929. [Google Scholar] [CrossRef] [PubMed]
  34. Port, M.; Idee’, J.-M.; Medina, C.; Robic, C.; Sabatou, M.; Corot, C. Efficiency, thermodynamic and kinetic stability of marketed gadolinium chelates and their possible clinical consequences: A critical review. Biometals 2008, 21, 469–490. [Google Scholar] [CrossRef] [PubMed]
  35. Schijf, J.; Christy, I.J. Effect of Mg and Ca on the Stability of the MRI Contrast Agent Gd–DTPA in Seawater. Front. Mar. Sci. 2018, 5, 111. [Google Scholar] [CrossRef] [Green Version]
  36. Idée, J.-M.; Port, M.; Raynal, I.; Schaefer, M.; Le Greneur, S.; Corot, C. Clinical and biological consequences of transmetallation induced by contrast agents for magnetic resonance imaging: A review. Fundam. Clin. Pharmacol. 2006, 20, 563–576. [Google Scholar] [CrossRef]
Figure 1. Absorbance at 510 nm as a function of the total Gd3+ concentration was linear in all buffers investigated. All of the buffers were 0.100 M.
Figure 1. Absorbance at 510 nm as a function of the total Gd3+ concentration was linear in all buffers investigated. All of the buffers were 0.100 M.
Environments 07 00069 g001
Figure 2. Results for the method of continuous variation in 0.100 M TEA buffer pH 8. The total concentration of PAR and Gd3+ was fixed at 1 × 10−4 M. The position of maximum absorbance (~35 ± 3%) is consistent with a 1:2 Gd3+: PAR ratio for the complex.
Figure 2. Results for the method of continuous variation in 0.100 M TEA buffer pH 8. The total concentration of PAR and Gd3+ was fixed at 1 × 10−4 M. The position of maximum absorbance (~35 ± 3%) is consistent with a 1:2 Gd3+: PAR ratio for the complex.
Environments 07 00069 g002
Figure 3. Graphical estimate of the Gd3+–oxalate stability constant. The linearity of this plot indicates that GdOx+ and Gd(PAR)2 are the predominant species under these experimental conditions.
Figure 3. Graphical estimate of the Gd3+–oxalate stability constant. The linearity of this plot indicates that GdOx+ and Gd(PAR)2 are the predominant species under these experimental conditions.
Environments 07 00069 g003
Figure 4. The reduction of Gd(PAR)2 concentration with increasing oxalate concentration at pH 8 at 25 °C is presented. The experimental data have been fit to Equation (4).
Figure 4. The reduction of Gd(PAR)2 concentration with increasing oxalate concentration at pH 8 at 25 °C is presented. The experimental data have been fit to Equation (4).
Environments 07 00069 g004
Figure 5. Summary of the Gd3+ titration data, at pH 8, for the various concentrations of tannic acid that were investigated using this technique.
Figure 5. Summary of the Gd3+ titration data, at pH 8, for the various concentrations of tannic acid that were investigated using this technique.
Environments 07 00069 g005
Figure 6. The Gd3+ titration data for the IHSS Leonardite humic acid at pH 8 and at different humic acid concentrations. The curvature is interpreted as a finite complexing capacity for the humic acid.
Figure 6. The Gd3+ titration data for the IHSS Leonardite humic acid at pH 8 and at different humic acid concentrations. The curvature is interpreted as a finite complexing capacity for the humic acid.
Environments 07 00069 g006
Figure 7. The titration data for the various concentrations of the TeraVita humic acid at pH 8 and at various humic acid concentrations are summarized. The curvature in the plot indicates finite complexing capacity.
Figure 7. The titration data for the various concentrations of the TeraVita humic acid at pH 8 and at various humic acid concentrations are summarized. The curvature in the plot indicates finite complexing capacity.
Environments 07 00069 g007
Figure 8. The data for the titration with Gd3+ of Pahokee Peat humic acid at pH 8 and at various concentrations is summarized in this figure. The curvature in the plot indicates finite complexing capacity.
Figure 8. The data for the titration with Gd3+ of Pahokee Peat humic acid at pH 8 and at various concentrations is summarized in this figure. The curvature in the plot indicates finite complexing capacity.
Environments 07 00069 g008
Figure 9. Summary of the titration of Fluka humic acid at pH 8 and various concentrations with Gd3+. The finite complexing capacity of the humic acid is apparent.
Figure 9. Summary of the titration of Fluka humic acid at pH 8 and various concentrations with Gd3+. The finite complexing capacity of the humic acid is apparent.
Environments 07 00069 g009
Table 1. Stability constants for the low molecular weight ligands determined in this study. Literature values for the stability constant are shown for comparison.
Table 1. Stability constants for the low molecular weight ligands determined in this study. Literature values for the stability constant are shown for comparison.
Ligandlog β (This Work)log β (Literature)
Oxalate4.96 ± 0.035.72 b,*
Citrate6.38 ± 0.057.83 a
Tartrate4.81 ± 0.024.15 a
Malonate3.92 ± 0.044.32 a
Lactate3.28 ± 0.092.89 a
Carbonate5.62 ± 0.046.10 b,*
Caffeate4.08 ± 0.08-
* Corrected to 0.100M ionic strength; a from Martell and Smith [28]; b from Gustafsson [29].

Share and Cite

MDPI and ACS Style

Steinberg, S.; Hodge, V.; Becerra-Hernandez, L. Investigation of the Interaction of Gadolinium with Several Organic Ligands and Humic Acid by Ligand Competition Using 4-(2-Pyridylazo)-Resorcinol (PAR). Environments 2020, 7, 69. https://0-doi-org.brum.beds.ac.uk/10.3390/environments7090069

AMA Style

Steinberg S, Hodge V, Becerra-Hernandez L. Investigation of the Interaction of Gadolinium with Several Organic Ligands and Humic Acid by Ligand Competition Using 4-(2-Pyridylazo)-Resorcinol (PAR). Environments. 2020; 7(9):69. https://0-doi-org.brum.beds.ac.uk/10.3390/environments7090069

Chicago/Turabian Style

Steinberg, Spencer, Vernon Hodge, and Luis Becerra-Hernandez. 2020. "Investigation of the Interaction of Gadolinium with Several Organic Ligands and Humic Acid by Ligand Competition Using 4-(2-Pyridylazo)-Resorcinol (PAR)" Environments 7, no. 9: 69. https://0-doi-org.brum.beds.ac.uk/10.3390/environments7090069

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop