International Journal of Inorganic Chemistry

Volume 2012 (2012), Article ID 265249, 9 pages

http://dx.doi.org/10.1155/2012/265249

## Speciation Studies of L-Histidine Complexes of Pb(II), Cd(II), and Hg(II) in DMSO-Water Mixtures

Department of Inorganic & Analytical Chemistry, School of Chemistry, Andhra University, Visakhapatnam 530003, India

Received 5 November 2011; Revised 21 January 2012; Accepted 23 January 2012

Academic Editor: Wei-Yin Sun

Copyright © 2012 K. Bharath Kumar Naik et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### 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 ML_{2}H_{4}, ML_{2}H_{3}, ML_{2}H_{2}, ML_{2}H, and ML_{2} and those for Hg(II) are ML_{2}H_{4}, ML_{2}H_{3}, ML_{2}, 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 Ca^{2+} 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 Hg^{2+} and promoting enzymatic formation of monomethylmercury [16–18]. 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 (O_{caboxyl}), imidazole imidonitrogen (N_{im}) and aminonitrogen (N_{am}) [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 [24–28].

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 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 CdL_{2}H_{4}, CdL_{2}H_{3,} CdL_{2}H_{2}, CdL_{2}H, and CdL_{2}. 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 (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, , 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 ( 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 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 ML_{2}H_{4}, ML_{2}H_{3}, ML_{2}H_{2}, ML_{2}H, and ML_{2} for Pb(II), and Cd(II). Hg(II) forms ML_{2}H_{4}, ML_{2}H_{3}, ML_{2}, and ML species. The formation of various his complex species is shown in the following equilibria.
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, species are formed from free-metal ion that interacts with form of the ligand [Equilibrium (1)]. The species may be formed from M(II) and or by the deprotonation of (Equilibria (2) & (3)). In the same manner, the species ML_{2}H_{2} may be formed either from free-metal ion and / (Equilibria (4) & (5)) or deprotonation of and (Equilibria (6) & (7)). At lower pH, the interaction of free-metal ion with or (Equilibria (4) & (5)) and at higher pH, deprotonation of , and (Equilibria (6) & (7)) results in the formation of ML_{2}H_{2} species. Species ML_{2}H is formed from M(II) and or or at lower pH, deprotonation of , , and species at higher pH (Equilibria (8)–(12)). The ML_{2} species are formed from interaction of free-metal ion, with or LH (Equilibria 15 & 16) or ML^{+} and LH (Equilibrium (17)) or deprotonation of , , , and ML_{2}H^{+} 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 or (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 ML_{2}H_{4}, ML_{2}H_{3}, ML_{2}H_{2}, ML_{2}H, and ML_{2}. ML_{2}H_{4}, ML_{2}H_{3}, ML_{2}, and ML species are detected the case of Hg(II). (2) The linear variation of stability constants as a function of 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.

#### References

- Y. Chen, C. Wang, and Z. Wang, “Residues and source identification of persistent organic pollutants in farmland soils irrigated by effluents from biological treatment plants,”
*Environment International*, vol. 31, no. 6, pp. 778–783, 2005. View at Publisher · View at Google Scholar · View at Scopus - R. L. Canfield Jr., C. R. Henderson, D. A. Cory-Slechta, C. Cox, T. A. Jusko, and B. P. Lanphear, “Intellectual impairment in children with blood lead concentrations below 10
*μ*g per deciliter,”*New England Journal of Medicine*, vol. 348, no. 16, pp. 1517–1526, 2003. View at Publisher · View at Google Scholar · View at Scopus - P. D. Alexander, B. J. Alloway, and A. M. Dourado, “Genotypic variations in the accumulation of Cd, Cu, Pb and Zn exhibited by six commonly grown vegetables,”
*Environmental Pollution*, vol. 144, no. 3, pp. 736–745, 2006. View at Publisher · View at Google Scholar · View at Scopus - M. Cecchi, C. Dumat, A. Alric, B. Felix-Faure, P. Pradere, and M. Guiresse, “Multi-metal contamination of a calcic cambisol by fallout from a lead-recycling plant,”
*Geoderma*, vol. 144, no. 1-2, pp. 287–298, 2008. View at Publisher · View at Google Scholar · View at Scopus - M. Arshad, J. Silvestre, E. Pinelli et al., “A field study of lead phytoextraction by various scented
*Pelargonium*cultivars,”*Chemosphere*, vol. 71, no. 11, pp. 2187–2192, 2008. View at Publisher · View at Google Scholar · View at Scopus - G. J. K. Komarnicki, “Lead and cadmium in indoor air and the urban environment,”
*Environmental Pollution*, vol. 136, no. 1, pp. 47–61, 2005. View at Publisher · View at Google Scholar · View at Scopus - J. M. Castagnetto, S. W. Hennessy, V. A. Roberts, E. D. Getzoff, J. A. Tainer, and M. E. Pique, “MDB: the Metalloprotein Database and Browser at The Scripps Research Institute,”
*Nucleic Acids Research*, vol. 30, no. 1, pp. 379–382, 2002. View at Google Scholar · View at Scopus - L. Zeneli, H. Paçarizi, N. M. Daci, M. Daci-Ajvazi, and A. Prenaj, “The effects of air pollution and smoking on cadmium concentration in human blood and correlation with biochemical parameters,”
*American Journal of Biochemistry and Biotechnology*, vol. 5, no. 2, pp. 59–62, 2009. View at Publisher · View at Google Scholar · View at Scopus - S. Takebayashi, S. Jimi, M. Segawa, and Y. Kiyoshi, “Cadmium induces osteomalacia mediated by proximal tubular atrophy and disturbances of phosphate reabsorption. A study of 11 autopsies,”
*Pathology Research and Practice*, vol. 196, no. 9, pp. 653–663, 2000. View at Google Scholar · View at Scopus - S. Chater, T. Douki, A. Favier, M. Sakly, and H. Abdelmelek, “Changes in antioxidant status and biochemical parameters after orally cadmium administration in females rats,”
*Acta Biologica Hungarica*, vol. 60, no. 1, pp. 79–88, 2009. View at Publisher · View at Google Scholar · View at Scopus - S. O. Asagba and F. O. Obi, “Effect of cadmium on kidney and liver cell membrane integrity and antioxidant enzyme status: implications for Warri River cadmium level,”
*Tropical Journal of Environmental Science and Health*, vol. 3, no. 1, pp. 33–39, 2000. View at Google Scholar - S. O. Asagba, “Alterations in activities of tissue enzymes in oral cadmium Toxicity,”
*Nigerian Journal of Science and Environment*, vol. 6, pp. 91–102, 2007. View at Google Scholar - D. Mergler, H. A. Anderson, L. H. M. Chan et al., “Methylmercury exposure and health effects in humans: a worldwide concern,”
*Ambio*, vol. 36, no. 1, pp. 3–11, 2007. View at Publisher · View at Google Scholar · View at Scopus - U. S. Environmental Protection Agency, 2010.
- E. B. Ekstrom and F. M. M. Morel, “Cobalt limitation of growth and mercury methylation in sulfate-reducing bacteria,”
*Environmental Science and Technology*, vol. 42, no. 1, pp. 93–99, 2008. View at Publisher · View at Google Scholar · View at Scopus - L. Landner, “Biochemical model for the biological methylation of mercury suggested from methylation studies in vivo with neurospora crassa,”
*Nature*, vol. 230, no. 5294, pp. 452–454, 1971. View at Publisher · View at Google Scholar · View at Scopus - J. K. Schaefer and F. M. M. Morel, “High methylation rates of mercury bound to cysteine by Geobacter sulfurreducens,”
*Nature Geoscience*, vol. 2, no. 2, pp. 123–126, 2009. View at Publisher · View at Google Scholar · View at Scopus - R. Sparling, “Biogeochemistry: mercury methylation made easy,”
*Nature Geoscience*, vol. 2, no. 2, pp. 92–93, 2009. View at Publisher · View at Google Scholar · View at Scopus - W. F. Fitzgerald, C. H. Lamborg, and C. R. Hammerschmidt, “Marine biogeochemical cycling of mercury,”
*Chemical Reviews*, vol. 107, no. 2, pp. 641–662, 2007. View at Publisher · View at Google Scholar · View at Scopus - J. G. Wiener, D. P. Krabbenhoft, G. H. Heinz, and A. M. Scheuhammer, “Ecotoxicology of mercury,” in
*Handbook of Ecotoxicology*, Chapter 16, p. 407, CRC Press, 2nd edition, 461. View at Google Scholar - P. Deschamps, P. P. Kulkarni, M. Gautam-Basak, and B. Sarkar, “The saga of copper(II)-L-histidine,”
*Coordination Chemistry Reviews*, vol. 249, no. 9-10, pp. 895–909, 2005. View at Publisher · View at Google Scholar · View at Scopus - T. E. Creighton,
*Encyclopedia of Molecular Biology*, vol. 2, Wiley, New York, NY, USA, 1999. - C. Guo Nan, W. Xiao Ping, D. Jian Ping, and C. Hong Qing, “A study on electrochemistry of histidine and its metabolites based on the diazo coupling reaction,”
*Talanta*, vol. 49, no. 2, pp. 319–330, 1999. View at Publisher · View at Google Scholar · View at Scopus - A. K. Vinnikova, R. C. Kukreja, and M. L. Hess, “Singlet oxygen-induced inhibition of cardiac sarcolemmal Na
^{+}K^{+}-ATPase,”*Journal of Molecular and Cellular Cardiology*, vol. 24, no. 5, pp. 465–470, 1992. View at Publisher · View at Google Scholar · View at Scopus - N. G. Sitton, J. S. Dixon, C. Astbury, R. J. Francis, H. A. Bird, and V. Wright, “Kinetic investigations into the possible cause of low serum histidine in rheumatoid arthritis,”
*Annals of the Rheumatic Diseases*, vol. 47, no. 1, pp. 48–52, 1988. View at Google Scholar · View at Scopus - E. Hitomi-Ohmura, N. Amano, Y. Aoyama, and A. Yoshida, “The effect of a histidine-excess diet on cholesterol synthesis and degradation in rats,”
*Lipids*, vol. 27, no. 10, pp. 755–760, 1992. View at Publisher · View at Google Scholar · View at Scopus - Y. Aoyama, T. Tsuda, E. Hitomi-Ohmura, and A. Yoshida, “Effect of dietary excess-histidine on fructose 1,6-bisphosphatase and 6-phosphofructokinase activities, and activation of fructose 1,6-bisphosphatase by basic amino acids in rat liver,”
*International Journal of Biochemistry*, vol. 24, no. 6, pp. 981–985, 1992. View at Publisher · View at Google Scholar · View at Scopus - Q. Cai, G. Takemura, and M. Ashraf, “Antioxidative properties of histidine and its effect on myocardial injury during ischemia/reperfusion in isolated rat heart,”
*Journal of Cardiovascular Pharmacology*, vol. 25, no. 1, pp. 147–155, 1995. View at Google Scholar · View at Scopus - H. Sigel, R. B. Martin, and R. Tribolet, “An estimation of the equivalent solution dielectric constant in the active-site cavity of metalloenzymes. Dependence of carboxylate—metal-ion complex stabilities on the polarity of mixed aqueous/organic solvents,”
*European Journal of Biochemistry*, vol. 152, no. 1, pp. 187–193, 1985. View at Google Scholar · View at Scopus - S. Raju, B. Ananda Kumar, G. Pushpa Raju, K. Bharath Kumar Naik, and G. N. Rao, “Protonation equilibria of L-glutamic acid and L-histidine in Dmso-water mixtures,”
*Current Chemical Research*. In press. - R. S. Rao and G. N. Rao,
*Computer Applications in Chemistry*, Himalaya Publishing House, Mumbai, India, 2005. - G. Gran, “Determination of the equivalence point in potentiometric titrations. Part II,”
*The Analyst*, vol. 77, no. 920, pp. 661–670, 1952. View at Publisher · View at Google Scholar · View at Scopus - G. Gran, “Equivalence volumes in potentiometric titrations,”
*Analytica Chimica Acta*, vol. 206, pp. 111–123, 1988. View at Google Scholar · View at Scopus - M. P. Latha, V. M. Rao, T. S. Rao, and G. N. Rao, “Chemical speciation of Pb(II), Cd(II), Hg(II), Co(II), Ni(II), Cu(II) and Zn(II) binary complexes of L-methionine in 1,2-propanediol-water mixtures,”
*Bulletin of the Chemical Society of Ethiopia*, vol. 21, no. 3, pp. 363–372, 2007. View at Google Scholar · View at Scopus - B. B. V. Sailaja, T. Kebede, G. N. Rao, and M. S. P. Rao, “Effect of micelles on speciation of ternary complexes of Uranyl(VI) with oxalic and malonic acids,”
*Proceedings of National Academy of Sciences, India*, vol. 74, no. 4, pp. 399–412, 2004. View at Google Scholar - N. Padmaja, M. S. Babu, G. N. Rao, R. S. Rao, and K. V. Ramana, “Effect of urea on the speciation of 1,10-phenanthroline complexes of toxic metal ions-computer augmented modelling studies,”
*Polyhedron*, vol. 9, no. 20, pp. 2497–2506, 1990. View at Google Scholar · View at Scopus - G. N. Rao,
*Complex equilibria of biological importance in aquo organic media-Computer augmented modeling studies*, Ph.D. thesis, Andhra University, Visakhapatnam, India, 1989. - P. Gans, A. Sabatini, and A. Vacca, “An improved computer program for the computation of formation constants from potentiometric data,”
*Inorganica Chimica Acta*, vol. 18, pp. 237–239, 1976. View at Google Scholar · View at Scopus - G. N. Rao and A. Ramakrishna, “Speciation studies of nickel(II) complexes of L-glutamine and succinic acid in urea-water mixtures,”
*Proceedings of National Academy of Sciences, India*, vol. 75, no. 4, pp. 245–248, 2005. View at Google Scholar