Abstract

Equilibrium study on complex formation of L-histidine with Pb(II), Cd(II), and Hg(II) has been investigated pH metrically in DMSO-water mixtures (0–60% v/v) at 303 K and 0.16 mol L−1 ionic strength. The predominant species detected for Pb(II) and Cd(II) are ML2H4, ML2H3, ML2H2, ML2H, and ML2 and those for Hg(II) are ML2H4, ML2H3, ML2, and ML. The appropriateness of experimental conditions is verified by introducing errors intentionally in the concentrations of ingredients. The models containing different numbers of species were refined by using the computer program MINIQUAD75. The best-fit chemical models were arrived at based on statistical parameters. The trend in variation of stability constants of the complexes with dielectric constant of the medium is attributed to the electrostatic and nonelectrostatic forces. The species distribution and the plausible equilibria for the formation of the species are also presented.

1. Introduction

The speciation study of toxic metal ion complexes is useful to understand the role played by the active site cavities in biological molecules and the bonding behavior of protein residues with the metal ion. The species refined and their relative concentrations under the experimental conditions represent the possible forms of aminoacids in biofluids.

Due to its numerous uses and high persistence, lead is a major environmental contaminant [1]. Lead is toxic even at low concentrations for living organisms, who can absorb it in various ways [2]. Lead intake by humans can be due to the consumption of crop plants grown on soils with high plant-available metal concentrations [3]. It can, however, migrate through the soil with dissolved organic matter [4] or mobilized by certain plants [5]. Moreover, carried from the air to the soils as fine particles, lead could be released more easily in soil solution [6].

Cadmium causes iron deficiency by binding to cysteine, glutamate, aspartate, and histidine ligands [7]. Cadmium inhibits enzymes that participate in bilirubine conjunction [8]. It increases urine Ca2+ excretion which can cause severe bone pathology [9]. The possible effects of long term low-level exposure to cadmium are of concern because it is readily distributed to tissues of liver and kidney, which are the main target organs in acute and chronic cadmium exposure [10, 11]. Other tissues involved in cadmium toxicity include the testis, heart, bone, eye, and brain [12].

Mercury is one of the most toxic elements and has negative health effects in human populations, highly dependent on fish consumption [13, 14]. Recent research concluded that neither vitamin B12 nor the acetyl-CoA pathways are required for bacterial methylation of mercury [15]. Cysteine enhances Hg methylation by facilitating uptake of Hg2+ and promoting enzymatic formation of monomethylmercury [1618]. Humans are exposed to Hg primarily as monomethylmercury [19], which is the form of Hg that accumulates readily in organisms and biomagnifies in food webs [20].

L-Histidine (His) is one of the strongest metal coordinating ligands among the aminoacids and plays an important role in the binding of metal ions by proteins. Most of the active sites of the biomolecules have histidine as one of the amino acid residues. It has three potential metal-binding sites, namely, carboxylate oxygen (Ocaboxyl), imidazole imidonitrogen (Nim) and aminonitrogen (Nam) [21] (Figure 1). Due to the high reactivity of its imidazole group, histidine residue is often found at the active site of enzymes and involved directly in catalysis [22]. It controls the transmission of metals in biological bases [23] and has been reported to act as a neurotransmitter or neuromodulator in mammalian central nervous system. His is involved in a large number of biochemical processes such as biosynthesis of histamine, secretion of prolactin and antidiuretic hormone, and production of red and white blood cells. It possesses vasodilating and hypotensive actions and may boost the activity of soothing alpha waves in the brain. His is used in the treatment of anemia, allergies, rheumatoid arthritis, and other inflammatory reactions [2428].

Dimethyl sulfoxide (DMSO) has the largest dielectric constant of the common dipolar aprotic solvents and it has been termed the super solvent. Another important phenomenon about the structure of the DMSO molecule is its trigonal pyramidal geometry. There is a highly directional lone pair of electrons at the apex of the pyramid, which helps in complexation.

His is an important component of active sites of biomolecules and there is a very likely hood for toxic metals to interact with His. Hence the speciation of the title systems has been studied in DMSO-water mixtures. DMSO has been chosen as a medium to decrease the dielectric constant of the medium, since the permittivity of active site cavities is very low [29]. Protonation equilibria of L-His in DMSO-water mixtures have already been studied in this laboratory [30].

2. Experimental

2.1. Materials

L-Histidine.HCl (E-Merck, Germany) solution (0.05 mol L−1) was prepared in triple-distilled deionised water by maintaining 0.05 mol L−1 nitric acid concentration to increase the solubility. DMSO (Qualigens, India) was used as received. 2 mol L−1 sodium nitrate (Qualigens, India) was prepared to maintain the ionic strength in the titrand. Solutions of Pb(II), Cd(II), and Hg(II) nitrates (0.05 mol L−1) were prepared by dissolving G.R. Grade (E-Merck, Germany) salts in triple-distilled water maintaining 0.05 mol L−1 nitric acid to suppress the hydrolysis of metal salts. All the solutions were standardized by standard methods. To assess the errors that might have crept into the determination of the concentrations, the data were subjected to analysis of variance of one way classification [31]. The strengths of alkali and mineral acid were determined using the Gran plot method [32, 33].

2.2. Apparatus

The titrimetric data were obtained using ELICO (Model LI-120) pH meter (readability 0.01), which was calibrated with 0.05 mol L−1 potassium hydrogen phthalate in acidic region and 0.01 mol L−1 borax solution in basic region. The glass electrode was equilibrated in a well stirred DMSO-water mixture containing the inert electrolyte. All the titrations were carried out in the medium containing varying concentrations of DMSO-water mixtures (0–60% v/v) by maintaining an ionic strength of 0.16 mol L−1 with sodium chloride at 303.0±0.1 K. The effect of variation in asymmetry potential, liquid junction potential, activity coefficient, sodium ion error, and dissolved carbon dioxide on the response of glass electrode was accounted for in the form of correction factor [34, 35].

2.3. Procedure

For the determination of stability constants of metal-ligand binary species, initially titrations of strong acid with alkali were carried out at regular intervals to check whether complete equilibration was achieved. Then the calomel electrode was refilled with DMSO-water mixture of equivalent composition as that of titrand. In each of the titrations, the titrand consisted of approximately 1 mmol mineral acid in a total volume of 50 mL. Titrations with different ratios (1 : 2.5, 1 : 3.75 and 1 : 5.0 in the case of Pb(II) and Cd(II) and 1 : 7.5, 1 : 8.5 and 1 : 10.0 in the case of Hg(II)) of metal-to-ligand were carried out with 0.4 mol L−1 sodium hydroxide. Other experimental details are given elsewhere [36].

2.4. Modeling Strategy

The computer program SCPHD [37] was used to calculate the correction factor. By using the pH-metric titration data, the binary stability constants were calculated with the computer program MINIQUAD75 [38], which exploit the advantage of the constrained least-squares method in the initial refinement and reliable convergence of Marquardt algorithm. During the refinement of binary systems, the correction factor and the protonation constants of histidine were fixed. The variation of stability constants with the dielectric constant of the medium was analyzed on electrostatic grounds on the basis of solute-solute and solute-solvent interactions.

3. Result and Discussion

The results of exhaustive modeling for the Cd(II)-His system in 40% DMSO-water mixture are given in Table 1. The models gave better statistics as the number of species was increased, confirming better fit. There was no further improvement in the fit on inclusion of some more species in the model containing CdL2H4, CdL2H3, CdL2H2, CdL2H, and CdL2. This indicates that the final model appropriately fits the experimental data. Such exhaustive modeling was performed for all the systems. The results of the final best-fit models that contain the stoichiometry of the complex species and their overall formation constants along with some of the important statistical parameters are given in Table 2. Very low-standard deviation in overall stability constants (log β) signifies the precision of these constants. The small values of 𝑈corr (sum of squares of deviations in concentrations of ingredients at all experimental points) corrected for degrees of freedom, small values of mean, standard deviation, and mean deviation for the systems are validated by the residual analysis [39].

3.1. Residual Analysis

In data analysis with least squares methods, the residuals (the differences between the experimental data and the data simulated based on model parameters) are assumed to follow Gaussian or normal distribution. When the data are fit into the models, the residuals should ideally be equal to zero. If statistical measures of the residuals and the errors assumed in the models are not significantly different from each other, the model is said to be adequate. Further, a model is considered adequate only if the residuals do not show any trend. Respecting the hypothesis that the errors are random, the residuals are tested for normal distribution. Such tests are χ2, Skewness, Kurtosis, and R-factor. These statistical parameters show that the best-fit models portray the metal-ligand species in DMSO-water mixtures, as discussed below.

In the present study, the χ2 values are less than the table values, and so the models are accepted. The kurtosis values in this study indicate that the residuals form leptokurtic pattern. The values of skewness recorded in Table 2 are between 0.06 and 1.39 for Pb(II), −0.93 and 2.45 for Cd(II) and −0.78, and 0.19 for Hg(II). These data evince that the residuals form part of a normal distribution. Hence, least square method can be applied to the present data. The sufficiency of the model is further evident from crystallographic 𝑅-values. These statistical parameters thus show that the best-fit models portray the metal-ligand species in DMSO media.

3.2. Effect of Systematic Errors on Best Fit Model

In order to rely upon the best-fit chemical model for critical evaluation and application under varied experimental conditions with different accuracies of data acquisition, an investigation was undertaken by introducing pessimistic errors in the influential parameters like concentrations of alkali, mineral acid, ligand, metal, log𝐹, and volume (Table 3). The order of the ingredients that influence the magnitudes of stability constants due to incorporation of errors is alkali > acid > metal > ligand > volume > log F. Some species were even rejected when errors were introduced in the concentrations. The rejection of some species and increased standard deviations in the stability constants on introduction of errors confirm the suitability of the experimental conditions (concentrations of ingredients) and choice of the best-fit models.

3.3. Effect of Solvent

The linear variation of stability constants (log β) of his complexes with Pb(II), Cd(II), and Hg(II) with variation of 1/𝐷 (𝐷 is the dielectric constant) of DMSO-water mixtures are given in Figure 2. DMSO is a polar aprotic and coordinating solvent. It is a structure former and it enhances the water structure in DMSO-water mixtures; hence, it removes water from coordination sphere of metal ions, making it more reactive towards the ligands. As a result, the stability of the complexes is expected to increase. At the same time, it is a coordinating solvent and competes with the ligands for coordinating the metals. This decreases the stability of the complexes. Hence, the stability of the complexes is expected to either increase or decrease. The linear variation indicates that electrostatic forces dominate the equilibrium process under the present experimental conditions.

3.4. Distribution Diagrams

Histidine is a tridentate ligand that has one dissociable (carboxylate group) and two associable (amino, imidazole) protons. The different forms of His are LH3+2  LH2+ LH, and L in the pH range below 4.0, 2.0–6.0, 6.0–10.0, and above 9.0, respectively. Hence, the plausible binary metal-ligand complexes can be predicted from these data. The present investigation reveals the existence of ML2H4, ML2H3, ML2H2, ML2H, and ML2 for Pb(II), and Cd(II). Hg(II) forms ML2H4, ML2H3, ML2, and ML species. The formation of various his complex species is shown in the following equilibria. M(II)+2LH32+ML2H44++2H+,(1)M(II)+2LH32+ML2H33++3H+,(2)ML2H44+ML2H33++H+,(3)M(II)+2LH32+ML2H22++H+,(4)M(II)+2LH2+ML2H22++2H+,(5)ML2H44+ML2H22++2H+,(6)ML2H33+ML2H22++H+,(7)M(II)+2LH32+ML2H++5H+,(8)M(II)+2LH2+ML2H++3H+,(9)ML2H44+ML2H++3H+,(10)ML2H33+ML2H++2H+,(11)ML2H22+ML2H++H+,(12)M(II)+LH32+ML++3H+,(13)M(II)+LH2+ML++2H+,(14)M(II)+2LH2+ML2+4H+,(15)M(II)+2LHML2+2H+,(16)ML++LHML2+H+,(17)ML2H44+ML2+4H+,(18)ML2H33+ML2+3H+,(19)ML2H22+ML2+2H+,(20)ML2H+ML2+H+.(21) The species distribution diagrams in DMSO-water mixtures are shown in Figure 3. They indicate that the his complexes of Pb(II), Cd(II), and Hg(II) are formed in the pH range 2.0–9.0. At lower pH, ML2H44+ species are formed from free-metal ion that interacts with LH32+ form of the ligand [Equilibrium (1)]. The species ML2H32+ may be formed from M(II) and LH32+ or by the deprotonation of ML2H44+ (Equilibria (2) & (3)). In the same manner, the species ML2H2 may be formed either from free-metal ion and LH3+/LH2+ (Equilibria (4) & (5)) or deprotonation of ML2H44+ and ML2H32+ (Equilibria (6) & (7)). At lower pH, the interaction of free-metal ion with LH32+ or LH2+ (Equilibria (4) & (5)) and at higher pH, deprotonation of ML2H44+, and ML2H32+ (Equilibria (6) & (7)) results in the formation of ML2H2 species. Species ML2H is formed from M(II) and LH32+ or LH2+ or at lower pH, deprotonation of ML2H44+, ML2H33+, and ML2H22+ species at higher pH (Equilibria (8)–(12)). The ML2 species are formed from interaction of free-metal ion, with LH2+ or LH (Equilibria 15 & 16) or ML+ and LH (Equilibrium (17)) or deprotonation of ML2H44+, ML2H33+, ML2H22+, and ML2H+ species (Equilibria (18)–(21)). In the case of Hg(II) (Figure 2(c)) the ML species are formed from interaction of free-metal ion and LH32+ or LH2+ (Equilibria (13) & (14)). It is formed at near neutral pH.

4. Conclusions

The following conclusions have been drawn from the modeling studies of the L-Histidine complexes of Pb(II), Cd(II) and Hg(II) in DMSO-water mixture. (1) L-Histidine forms both protonated and unprotonated complexes under pH range 2.0–9.0. The species detected for the toxic metals Pb(II) and Cd(II) are ML2H4, ML2H3, ML2H2, ML2H, and ML2. ML2H4, ML2H3, ML2, and ML species are detected the case of Hg(II). (2) The linear variation of stability constants as a function of 1/𝐷 of the medium indicates the dominance of electrostatic forces over nonelectrostatic forces. The linear increasing trend with DMSO content supports the dominance of the structure forming nature of DMSO over its complexing ability. (3) The order of ingredients in influencing the magnitudes of stability constants due to incorporation of errors in their concentrations is alkali > acid > metal > ligand > volume > log F. (4) Some species are stabilized due to electrostatic interactions and some are destabilized due to the decreased dielectric constant.

Acknowledgment

The authors (K. B. K. Naik, B. A. Kumar, S. Raju) thank The University Grants Commission, New Delhi, India for financial support under Faculty Development Programme.