Abstract
The complexation equilibria of folic acid (FA) and glycine (Gly) were studied in aqueous solutions at room temperature (ca. 298.15 K) and in fixed ionic strength (0.15 mol.dm−3 NaNO3), by means of potentiometry technique. The ferric (Fe3+), ferrous (Fe2+), calcium (Ca2+), and zinc (Zn2+) complex-forming capacities of FA and Gly, and their overall stability constants in aqueous solutions were obtained by applying the potentiometric data in Hyperquad 2008 computer program. From the determined stability constants of the metal complex species, the concentration distribution of the various metal ion complex species involving FA and Gly in solutions was estimated using HySS 2009 software. The complex species distribution diagrams were explained and discussed. Supplementary, the Gibbs free energies and the molecular structures of the formed complexes were evaluated and predicted using Gaussian 09 software molecular modeling and calculations.
GRAPHICAL ABSTRACT
Introduction
Vitamin B9 (folic acid, pteroyl-l-monoglutamic acid, FA), a water-soluble compound of the vitamin B9 group [Citation1], is added to many food products to prevent folate deficiency in individuals [Citation2]. Supplementation with FA is particularly important in pregnant women, as insufficient FA can cause neural tube defects in the developing fetus [Citation3]. Moreover, folate deficiency is the most common cause of anemia after iron deficiency and is thought to increase the likelihood of heart attacks and strokes [Citation4].
Glycine (Gly) is a non-essential amino acid that can be produced by body at a rate of 125 mg kg−1 body weight with adequate total protein intake. Gly is generated from serine by glycine hydroxymethyl transferase, and threonine by L-threonine dehydrogenase and 2-amino-3-ketobutyrate coenzyme A ligase. For dietary source, Gly can be found in all foods, particularly beef (55 mg g−1 protein) and pork (57 mg g−1). Daily oxidation of Gly is about 90 mg kg−1 body weight with an energy yield of 2.011 kcal g−1 [Citation5]. As an amino acid, Gly is the building blocks of protein and nucleotide (purine). As an antioxidant, Gly is one of the amino acids that make vital antioxidant peptide glutathione. Moreover, Gly is a major inhibitory neurotransmitter in brain [Citation6].
Vitamin B9 (FA) is one of the most biologically important compounds, and information about their metal ions complexation properties are rare in the literature. However, the interaction between some transition metal ions and vitamin B9 resulted in the formation of a number of binary transition metals–FA complexes in which number of drugs can interfere FA supplement [Citation7–14]. The present study was done to determine the global stoichiometric stability constants of calcium/zinc/iron metal ions, FA, and Gly binary and mixed ligand complexes using potentiometric technique, and by applying the row data on Hyperquad 2008 computer program. The elucidation of the plausible molecular structures of the vitamin B9–Gly complexes in solutions and the possible binding sites of FA to bind calcium, zinc, and iron metal ions, and the Gibbs free energy calculations of the complex species were done using Gaussian 09 molecular modeling software. In addition, calcium, zinc, and iron metal speciation studies were performed using Hyperquad Simulation and Speciation (HySS 2008) program, in order to evaluate the suitability of supplementation of vitamin B9 in pregnant women by studying the interaction of FA with iron, zinc, and calcium metal ions in the presence of Gly, the simplest amino acid.
Experimental
Materials and solutions
All the chemical materials used are of analytical reagent grade and were used without further purifications. FA and Gly were purchased from Sigma Aldrich (USA). A nitric acid (Panreac, Spain) solution was prepared and used after being standardized. Carbonate-free sodium hydroxide (NaOH) from Across Organics (USA) was standardized with potassium hydrogen phthalate (Sigma Aldrich, USA). All metal salts [calcium chloride dihydrate (CaCl2.2H2O) from Sigma Aldrich (USA), ferric chloride hexahydrate (FeCl3.6H2O), ferrous sulfate (Fe(SO4)2), and zinc nitrate hexahydrate (Zn(NO3)2.6H2O) from Across Organics (USA)] were weighted accurately before preparing the solutions. To maintain the ionic strength in the solution, sodium nitrate (NaNO3) from Across Organics (USA) was used. All solutions used throughout the experiments were prepared freshly in ultrapure water obtained from a NANO-pure ultrapure water system in which water was distilled and deionized with a resistance of 18.3 MΩ cm−1.
Apparatus and procedure
The pH-potentiometric titrations were performed using a Metrohm 836 Titrando with a model 801 magnetic stirrer, coupled with a Dosino burette model 803. The electrode response can be read to the third decimal place in terms of pH units with a precision of ±0.001 and the potential with a precision of ±0.1 mV. The titroprocessor was coupled to a personal computer and the titration software TINET version 2.4 was used to control the titration and data acquisition [Citation15]. The pH titrations were carried out in a 150 cm3 commercial glass vessel. The ionic strength of the solutions was maintained at a constant level by using the desired concentration of NaNO3 solution as supporting electrolyte. The pH meter was calibrated with standard buffer solutions (pH 4.0 and 7.0) before and after each series of pH measurements. The results of strong acid versus alkali titrations were analyzed using a computer program (GLEE, glass electrode evaluation) that has been used for the calibration of glass electrodes by means of a strong acid–strong base titration [Citation16]. The program GLEE has been developed as part of the Hyperquad suite of programs for stability constant determination [Citation16]. This program provides an estimate of the carbonate concentration of the base, the pseudo-Nernstian standard potential and slope of the electrode, and optionally the concentration of the base and the values of pKw (p at 25°C).
For the protonation constants estimation of FA and Gly, the following solutions were prepared (total volume 50 cm3) and titrated potentiometrically against a standard carbonate-free NaOH (0.1 mol dm−3) solution: (a) 0.003 mol dm−3 HNO3+0.15 mol dm−3 NaNO3; (b) Solution (a) +0.001 mol dm−3 FA; and (c) Solution (a) +0.001 mol dm−3 Gly. For the determination of the binary metal complexes, solutions containing FA, Gly, and metal ions (calcium, zinc, and iron) were titrated at 1:1, 1:2, and 1:3 metal ion to ligand mole ratios to fulfill the maximum coordination number of the metal ion, and for ternary systems, the 1:1:1 metal ion/FA/Gly ratio was used (). The concentration of metal ion and ligand solutions in the titrated samples was always the same and varied in the range 0.00004–0.0001 mol·dm−3 (). Each solution was thermostated at 298.15 K, where the solutions were left to stand for several minutes before titration. A magnetic stirrer was used during all titrations. Each titration was repeated at least 3 times under carefully controlled experimental conditions. Typically, more than 80 pH readings (points of potentiometric measurements) were collected and taken into account for each titration.
Table 1. Metal and ligand concentrations that were used to determine the overall formation constants of the complexes at 298.15 K and in 0.15 mol·dm−3 NaNO3 solutions.
Data analysis
The Hyperquad2008 program was used to obtain protonation and complex stability constants from potentiometric data [Citation18,Citation19]. For this purpose, a fitting criterion based on the minimization of the nonlinear least-squares sum defined by the difference between the calculated and the experimental data of the titration curves was used (equation (1)):
Hyperquad 2008 program permits the determination of formation constants of different complex species that can be formed in the aqueous solution simultaneously [Citation17,Citation18]. For each experimental titration point, the concentrations of the species in solutions are defined by the following mass balance equations:
In this work, Gaussian 09 program was used to predict the binding site of the ligands that contributed to complex formation and the structure of the complexes [Citation19]. The model structure of these complexes were optimized using density functional theory (DFT) with the Becke's three-parameter hybrid method [Citation20] in correlation with the Lee–Yang–Parr (LYP) (B3LYP) method, and the 6-31G basis set was used with diffuse and polarization +d for hydrogen and heavy atoms [Citation21–25]. Along with geometry optimization, frequency analysis was done to obtain the thermo-chemical properties of the complexes. The B3LYP/6−31+G(d) basis set was chosen since it has been shown to be one of the suitable basis sets to approach the molecular orbital of metal complexes [Citation22–26]. The validity of the structural optimization was checked by using normal-mode frequency analysis, in which the real minimum structure must exhibit positive value for all frequencies. The modeling was performed by employing a number of assumptions as simplification, such as the use of fully deprotonated ligands, only considering the complex formation of one metal ion and one molecule of ligand, not considering the usage of salt to maintain the solution's ionic strength [Citation22–25].
Results and discussions
From the refinement of the pH-potentiometric titration data by Hyperquad 2008 program, the protonation constants as overall formation constants (logβ) were obtained according to the following equations:
Four protonation constants as overall equilibrium constants [log β1, log β2, log β3, and log (related to α–COOH), 12.1547 (related to β–COOH), 16.3104 (related to NH2 amine group – site 3), and 19.9647 (related NH/CO (N-3 site) fragment of the pteridine ring)] of FA were determined by fitting the pH-potentiometric experimental data of FA on Hyperquad 2008 program ( and ). Another two low protonation constant values (pKa<−1.5, 0.20) of FA investigated can be associated with (N-5 site) and (N-10 site), respectively [Citation26]. As shown, these two values are very low and not significant in the physiological pH region. Thus, these values are not used in the pH-potentiometric data analysis measured in the range of 2.00≤ pH ≤ 12.5. The graphical distribution species diagrams in showed that the protonated FA species FAH4 started to form in very low acidic conditions (pH <2.5). While the FA species FAH3 was formed in acidic conditions (pH ∼2.5–5), conversely, the species FAH2 started to form at normal acidic conditions and extremely formed at neutral conditions (pH ∼3–7). The folate species FAH started to form in normal acidic conditions (pH ∼4) and completely formed in neutral and alkaline conditions (pH ∼5.5–9.5).
The global stability constants (log βpqrs) for the binary and mixed ligand systems involving metal ions (Fe3+, Fe2+, Ca+, and Zn+2) with FA and Gly as well as the pH intervals in which the data were collected are presented in and . The complexation between metal ion (M) and FA and Gly can be described by the following general equilibrium equations:
Table 2. Global stability constants (logβpqrs) of metal ion complexes involving FA and Gly at 298.15°K and in 0.15 mol·dm−3 NaNO3 solutions.
Table 3. Overall formation constants of metal ions (M)–FA–Gly mixed ligand complexes at 298.15°K and in 0.15 mol·dm−3 NaNO3 solutions.
In all titrations, the experimental pH titration curves recorded for the binary and ternary metal ion–FA–Gly systems, when compared with the pH titration curves of FA and Gly ligands alone, prove that complex formation between the metal ions studied (Fe3+, Fe2+, Ca+, and Zn+2) and the two ligands (FA and Gly) did occur. The ternary complex formation may proceed through either a stepwise or a simultaneous mechanism depending on the chelating potentials of FA and Gly. No precipitation occurred in all solutions during the titrations. However, slight turbidity was observed in the case of the Ca+–FA–Gly ternary system at pH ∼6.71–9.37.
The bioavailability of metals and their physiological and toxicological effects depend on the actual species present. The equilibrium distributions of various complex species formed in water solutions are shown as functions of pH (), and showed the following:
For the iron (III) systems involving FA and Gly, the protonated and deprotonated binary complexes, [, [], [FeGlyH−1], and [FeGly]+ start to form in very low acidic media and reach maximum of 92%, 95%, 5%, and 18% at pH values 3.7, 5.7, 4, and 7.3, respectively. The formed complexes [, [FeFA2H2], [FeGlyH−1], and [FeGly]+ disappeared at pH values 6.2, 7.9, 5.8, and 10, respectively. While the binary and mixed ligand complexes [FeFAH−1], [, [, [FeFAGlyH4], and [ start to form in slightly acidic media and the two complex species [FeFAH−1] and [ achieved maximum formation (88% and 84%) at pH values 6.8 and 9, respectively. The two mixed ligand complex species [FeFAGlyH4] and [ achieved maximum complex formation (18% and 13%) in neutral media and disappeared in alkali media.
For the divalent iron complexes involving FA and Gly, the protonated and deprotonated binary complexes [, [], [, [, and [FeGly]+ also start to form in very low acidic media and reach maximum of 64%, 72%, 31%, 51%, and 67% at pH values 4.8, 4.9, 4.2, 7, and 7.7, respectively. While the binary and mixed ligand complex species [FeFAH−1], [, and [FeFAGlyH4] start to form in slightly acidic media and achieved maximum formation (47%, 26%, and 18%) in neutral media, respectively. These complexes disappeared in slightly alkaline media.
The calcium binary and mixed ligand complex species [, [CaGly]+, [CaFAH−1], [, [CaGlyH−1], and [CaFAGlyH4] started to form in very low acidic media. The complex species [, [CaGly]+, and [CaFAH−1] achieved maximum complex formation (53%, 52%, and 42%) in slightly acidic media, and ended in slightly neutral and alkaline media. The mixed ligand complex formation [ started to form at pH ∼6.7 and reached maximum (54%) in neutral solution and disappeared in slightly alkaline media.
The binary and mixed ligand complexes involving zinc metal ion [, [ZnFA2H2], [ZnFAH−1], [, [, and [ZnFAGlyH4] started to formed in very low acidic media, and both [ and [ZnFA2H2] complex species achieved complex formation maximum (60% and 53%) in low acidic media, and disappeared at pH values ∼4 and 7, respectively. While the complexes [ZnFAH−1] and [ completely formed (47% and 51%) in slightly neutral media, the complex species [ZnGly]+, [, and [ZnGly]− started to form in neutral and slightly alkaline media, and achieved their maximum complex formation in alkaline media.
The logarithms of the overall stability constants (log βpqrs) of the binary and mixed ligand complexes investigated are given in and . From these data, the following was observed:
The global stability constants data of the binary complexes involving FA, and Fe3+, Fe2+, Ca+, and Zn2+ metal ions are decreased in the order of Fe complexes.
The stability constants of ferric (III), ferrous (II), calcium (I), and zinc (II) binary complexes of Gly agreed well with my previous data [Citation26].
By comparing the stability constants of binary and mixed ligand complexes, it was demonstrated that mixed ligand complexation is responsible for the stabilization of the mixed complexes.
The negative values indicate the absence of undeprotonated amide nitrogen contribution in the complexation of metal ions with Gly. Gly, the smallest and the simplest ligand, can form five-membered ring complexes with metal ions via deprotonated amine and carboxyl group of the amino acid during the ternary complex formation (scheme 2). The phenomena of amide nitrogen contribution in the complex stability were further supported by the optimization and frequency analysis using Gaussian 09. From the geometry optimization, the structure of the complexes that formed can be predicted (scheme 2), and while the thermochemistry properties such as Gibbs free energy (G) can be obtained from the frequency analysis [Citation27]. In this program, the complex equilibrium was presented as the Gibbs free energy of reaction (ΔrG) that was calculated using the program. According to the aforementioned correlation, the more negative the ΔrG, the larger the equilibrium constant (Kc) value will be, and the complex was found to be more stable ().
Table 4. Complex species Gibbs free energy of reaction obtained from optimization and frequency calculation by using Gaussian 09 program.
Funding
This work was financially supported by the King Khalid University, Abha, Kingdom of Saudi Arabia, through the project number [KKU_S112_85].
References
- T. Tamura, M.F. Picciano. Am. J. Clin. Nutr., 83, 993 (2006).
- O. Stanger. Curr. Drug. Metab., 3, 211 (2002).
- M. McCarthy. The Lancet, 347, 682 (1996).
- M. Lucock. Mol. Genet. Metab., 71, 121 (2000).
- J.C. Bulinski. Int. Rev. Cytol., 103, 281 (1986).
- R.Y. Gundersen, P. Vaagenes, T. Breivik, F. Fonnum, P.K. Opstad. Acta Anaesthesiol Scand, 49, 1108 (2005).
- M. Poe. J. Biol. Chem., 252, 3724 (1977).
- J. Ghasemi, S. Ghobadi, B. Abbasi, M. Kubista. J. Korean Chem. Soc., 49, 269 (2005).
- Z. Szakács, B. Noszál. Electrophoresis, 27, 3399 (2006).
- N. Nombona, W. Chidawanyika, T. Nyokong. Polyhedron, 30, 654 (2011).
- R.C. Mehrotra, R. Bohra. Metal Carboxylates, Academic Press, London (1983).
- R.C. Mehrotra, A. Singh. Prog. Inorg. Chem., 46, 239 (1997).
- E. Hamed, M.S. Attia, K. Bassiouny. Bioinorg. Chem. Appl., 6 (2009). Article ID 979680.
- M.G. Abd El-Wahed, M.S. Refat, S.M. El-Megharbel. Spectrochim Acta, Part A, 70, 916 (2008).
- A.G. Metrohm. Instructions for Use for 6.6012.X40 Software TiNet 2.4 CH-9101, pp. 1–148, Metrohm AG, Herisau (1999).
- P. Gans, O. Sullivan, B. Glee. Talanta, 51, 33 (2000).
- P. Gans, A. Sabatini, A. Vacca. Talanta, 43, 1739 (1996).
- L. Alderighi, P. Gans, A. Ienco, D. Peters, A. Sabatini, A. Vacca. Coord. Chem. Rev., 184, 311 (1999).
- M.J. Frisch, G.W. Trucks, H.B. Schlegel, G.E. Scuseria, M.A. Robb, J.R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci, G.A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H.P. Hratchian, A.F. Izmaylov, J. Bloino, G. Zheng, J.L. Sonnenberg, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, J.A.J. Montgomery, J.E. Peralta, F. Ogliaro, M. Bearpark, J.J. Heyd, E. Brothers, K.N. Kudin, V.N. Staroverov, R. Kobayashi, J. Normand, K. Raghavachari, A. Rendell, J.C. Burant, S.S. Iyengar, J. Tomasi, M. Cossi, N. rega, J.M. Millam, M. Klene, J.E. Knox, J.B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R.E. Stratmann, O. Yazyev, A.J. Austin, R. Cammi, C. Pomelli, J.W. Ochterski, E.L. Martin, K. Morokuma, V.G. Zakrzewski, G.A. Voth, P. Salvador, J.J. Dannenberg, S. Dapprich, A.D. Daniels, O. Farkas, J.B. Foresman, J.V. Ortiz, J. Cioslowski, D.J. Fox. Gaussian 09, Gaussian, Inc., Wallingford, CT (2009).
- A.D. Becke. J. Chem. Phys., 98, 5648 (1993).
- C. Lee, W. Yang, R.G. Parr. Phys. Rev. B, 37, 785 (1998).
- V.A. Rassolov, J.A. Pople, M.A. Ratner, T.L. Windus. J. Chem. Phys., 109, 1223 (1998).
- J.M. Ramos, O.J.F. Versiane, C.A.T. Soto. Spectrochim. Acta Part A, 68, 1370 (2007).
- J.M. Ramos, O. Versiane, J. Felcman, C.A.T. Soto. Spectrochim. Acta Part A, 72, 182 (2009).
- D.V. Chachkov, O.V. Mikhailov. Russ. J. Inorg. Chem., 54, 1952 (2009).
- A.Y. Rajhi, Y.-H. Ju, A.E. Angkawijaya, A.E. Fazary. J. Solution Chem., 42, 2409 (2013).
- A.E. Fazary. J. Chem. Eng. Data, 58, 2219 (2013).