Metal Complexation in Xylem Fluid : II. THEORETICAL EQUILIBRIUM MODEL AND COMPUTATIONAL COMPUTER PROGRAM

Theoretical considerations of metal complex formation in aqueous solutions were used to develop a computer program (CHELATE) to calculate all equilibrium species (free metal ions, metal complexes, etc.) in any user-defined system, such as xylem fluid. Mass-balance equations were established to describe each free metal ion and each free ligand concentration as a function of solution pH, total metal or total ligand, hydrogen-association constants, and the stability constants of known metal complexes. A default data base can be altered by the user to define any desired system covered by the stored equilibrium data. The program can currently handle nine metal ions, 35 ligands, and 500 complex species. The validity of the program was confirmed by using experimental test systems in which free-metal ion activity measurements were made with ion-selective electrodes.     

A simple method for calculation of free metal ion concentrations in the presence of a chelating agent and multiple metal ions by successive approximation

Calculation of free metal ion concentrations in the presence of a chelating agent (ligand) and multiple metal ions is complicated. In this paper we describe a simple method for calculation of free ion concentrations from given total ion concentrations. The outline of this method is as follows: i) Using an arbitrarily chosen provisional value of free ligand concentration (p-Lf), calculate a total ligand concentration (Lt). ii) Divide the p-Lf by the ratio, the calculated Lt (c-Lt)/the specified Lt (s-Lt). Take the resulting value as the next p-Lf and repeat the calculation until c-Lt is close enough to s-Lt. At this point, p-Lf is supposed to be set to a good approximation of true free ligand concentration. iii) Finally, calculate free metal concentrations from the above approximation.     

CHELATOR: an improved method for computing metal ion concentrations in physiological solutions.

An algorithm is presented for the calculation of metal ion concentrations from given total metal concentrations (and vice versa) in physiological media containing metal-chelating compounds. In such media, conditions differ from those used for stability constant determination of metal-chelator equilibria; therefore calculated metal ion concentrations are incorrect. We recompute stability constants to reflect the effects of ionic strength and temperature of physiological solutions. Twelve different equilibria can be considered per metal-chelator pair. The computer program also calculates the contribution of ionized species of metals, chelator, complexes and pH buffers to ionic strength. Measurements with a Ca-selective electrode and with fura-2 show that calculated ionic Ca2+ concentrations are correct from 10 nM up to the millimolar range. The importance of the correct calculation of metal ion concentrations in physiological experiments is demonstrated by data, and derived kinetic parameters, on Na+/Ca2+ exchange and the ATP-dependent Ca2+ pump of enterocyte plasma membrane vesicles. The program is written in Turbo Pascal and will run on IBM-compatible computers. It is menu-driven and supports the use of a Microsoft mouse.     

Transporters of ligands for essential metal ions in plants

Essential metals are required for healthy plant growth but can be toxic when present in excess. Therefore plants have mechanisms of metal homeostasis which involve coordination of metal ion transporters for uptake, translocation and compartmentalization. However, very little metal in plants is thought to exist as free ions. A number of small, organic molecules have been implicated in metal ion homeostasis as metal ion ligands to facilitate uptake and transport of metal ions with low solubility and also as chelators implicated in sequestration for metal tolerance and storage. Ligands for a number of essential metals have been identified and proteins involved in the transport of these ligands and of metal-ligand complexes have been characterized. Here we review recent advances in understanding the role of mugineic acid, nicotianamine, organic acids (citrate and malate), histidine and phytate as ligands for iron (Fe), zinc (Zn), copper (Cu), manganese (Mn) and nickel (Ni) in plants, and the proteins identified as their transporters.     

Inhibition of yeast growth by molybdenum-hydroxylamido complexes correlates with their presence in media at differing pH values

The effects of Mo-hydroxylamido complexes on cell growth were determined in Saccharomyces cerevisiae to investigate the biological effects of four different Mo complexes as a function of pH. Studies with yeast, an eukaryotic cell, are particularly suited to examine growth at different pH values because this organism grows well from pH 3 to 6.5. Studies can therefore be performed both in the presence of intact complexes and when the complexes have hydrolyzed to ligand and free metal ion. One of the complexes we examined was structurally characterized by X-ray crystallography. Yeast growth was inhibited in media solutions containing added Mo-dialkylhydroxylamido complexes at pH 3-7. When combining the yeast growth studies with a systematic study of the Mo-hydroxylamido complexes' stability as a function of pH and an examination of their speciation in yeast media, the effects of intact complexes can be distinguished from that of ligand and metal. This is possible because different effects are observed with complex present than when ligand or metal alone is present. At pH 3, the growth inhibition is attributed to the forms of molybdate ion that exist in solution because most of the complexes have hydrolyzed to oxomolybdate and ligand. The monoalkylhydroxylamine ligand inhibited yeast growth at pH 5, 6 and 7, while the dialkylhydroxylamine ligands had little effect on yeast growth. Growth inhibition of the Mo-dialkylhydroxylamido complexes is observed when a complex exists in the media. A complex that is inert to ligand exchange is not effective even at pH 3 where other Mo-hydroxylamido complexes show growth inhibition as molybdate. These results show that the formation of some Mo complexes can protect yeast from the growth inhibition observed when either the ligand or Mo salt alone are present.     

Calculation of the concentrations of free cations and cation-ligand complexes in solutions containing multiple divalent cations and ligands.

The method described permits the computation of the concentrations of free ions and ion-ligand complexes in a solution containing arbitrary numbers of divalent cations and ligands. It is required that the pH be known, along with appropriate sets of ligand-hydrogen and ligand-divalent cation concentration binding constants. It is assumed that these sets of constants are chosen to be consistent with the ionic strength of the complete solution which contains the divalent cations and ligands. The technique is an iterative one which provides upper and lower bounds for the values of the unknowns. The method does not require initial guesses at the values of the unknowns, and it gives correct answers even when the concentrations involved are many orders of magnitude apart. The present formulation of the problem is restricted to the case where only one cation can bind to a given ligand at any one time. The method is applicable to large molecules with multiple "sub-ligands" provided these sub-ligands are independent in their function as ion-binding sites. These sub-ligands need not all have the same properties. It is also shown that a simple modification of the method permits the determination of the subset of total ion concentrations that are required in order to produce a specified subset of free ion concentrations. The modifications required to include monovalent cation binding are presented in outline form.     

The determination of the binding constant of metalloenzymes for their active site metal ion from ligand inhibition data: Theoretical analysis and application to the inhibition of thermolysin by 1,10-phenanthroline

An equation is found relating the fractional activity, (v/v₀), of an enzyme assay mixture to the total concentrations of metalloenzyme, active site metal ion, metal-binding ligand and substrate and the stability constants of the complexes present. When (v/v₀) is measured as a function of the total ligand concentration, this equation offers a way of data-plotting which yields straight lines and permits the calculation of the metal-binding constant K_ME from either the slope or the intercept, provided that mixed complexes (enzyme-metal ion-ligand) do not contribute significantly to the change in (v/v₀). Since deviations from linearity occur in the latter case, the proposed inhibition plot serves as a diagnostic tool for the recognition of such complexes. Application to the inhibition of thermolysin by 1,10-phenanthroline gives a value of 2.1 × 10¹¹m⁻¹ for K_ZnE, the binding constant of the active site zinc ion, at pH 7.50, 25°C and ionic strength 0.1. The equation also allows the rapid calculation of the ligand concentration necessary to attain a desired degree of inhibition when the total enzyme and active site metal ion concentrations of the solution are known.     

Conformational studies of A23187 with mono-, di- and trivalent metal ions by circular dichroism spectroscopy

CD studies carried out on A23187 indicate a solvent-dependent conformation for the free acid. Alkali metal ions were found to bind to the ionophore weakly. Divalent metal ions such as Mg2+, Ca2+, Sr2+, Ba2+ and Co2+ and trivalent lanthanide metal ions like La3+ were found to form predominantly 2:1 (ionophore-metal ion) complexes at low concentrations of metal ions, but both 2:1 and 1:1 complexes were formed with increasing salt concentration. Mg2+ and Co2+ exhibit similar CD behaviour that differs from that observed for the other divalent and lanthanide metal ions. The structure of 2:1 complexes involves two ligand molecules coordinated to the metal ion through the carboxylate oxygen, benzoxazole nitrogen and keto-pyrrole oxygen from each ligand molecule along with one or more solvent molecules. Values of the binding constant were determined for 2:1 complexes of the ionophore with divalent and lanthanide metal ions.     

UPTAKE OF CADMIUM BY FRESHWATER GREEN ALGAE: EFFECTS OF PH AND AQUATIC HUMIC SUBSTANCES1

The effects of humic substances and low pH on short‐term Cd uptake by Pseudokirchneriella subcapitata (Korshikov) Hindak and Chlamydomonas reinhardtii Dang were investigated under defined exposure conditions. The uptake experiments were run in the presence of either a synthetic organic ligand (nitrilotriacetate) or natural organic ligands (Suwannee River fulvic or humic acid). An ion‐exchange method was used to measure the free Cd²⁺ concentrations in the exposure solutions. At pH 5, measured free Cd²⁺ concentrations agreed with estimations made using the geochemical equilibrium model WHAM, but at pH 7 the model overestimated complexation by both Suwannee River fulvic and humic acids compared with the ion‐exchange measurements. Consistent with the metal internalization step being rate limiting for overall short‐term uptake, intracellular Cd uptake was linear for exposure times less than 20 min at pH 5 or pH 7 for both algal species. After taking into account complexation of Cd in solution, Suwannee River humic substances had no additional effects on cadmium uptake at pH 7, as would be predicted by the free ion model. This absence of effects other than complexation persisted at pH 5, where the tendency of humic substances to adsorb to the algal cell surface is favored. Changes in pH strongly influenced Cd uptake, with the intracellular flux of Cd being at least 20 times lower at pH 5 than at pH 7 for P. subcapitata. Our results support models such as the free ion model or the biotic ligand model, in which humic substances act indirectly on Cd uptake by reducing the bioavailability of Cd by complexation in solution.     

Chromatographic separations of serum proteins on immobilized metal ion stationary phases

The separation of proteins on stationary phases consisting of a bound organic chelator and a chelated divalent transition metal has been studied as a function of (A) metal ion species; (B) mobile phase composition and pH; and (C) anion and cation concentration. Optimum separation was observed at alkaline pH on chelated nickel stationary phases. Ammonium and Tris salts reduced the affinity of the metal chelate packing for serum proteins. Halide ions caused the proteins to be more strongly bound to the stationary phase. High salt concentrations had only a small effect on the binding of serum proteins in the absence of amine containing buffers or salts. It was also observed that the ease of elution and the recovery of protein were dependent on pH and upon the presence of halides. The general order of elution of serum proteins, based on isoelectric focusing, was independent of metal ion species and elution conditions, suggesting that a single mechanism or a unique sequence of mechanisms was operative. The results suggest that ligand exchange is the major mechanism of separation under basic conditions and that hydrophobic effects are the result of the competition of nonnitrogen ions with ammonium ions or amines for ligand binding sites modifying or participating in protein binding. Protein binding studies under weak acidic conditions are also presented although the mechanism responsible for protein binding is unclear.     

Expanding the range of free calcium regulation in biological solutions

Many biological systems use ethylene glycol bis (beta-aminoethylether)-N,N,N',N'-tetraacetic acid (EGTA) to regulate the free calcium concentration ([Ca(2+)](free)) in the presence of physiological levels of free Mg(2+) ([Mg(2+)](free)). Frequently, it is necessary to work at [Ca(2+)](free) beyond EGTA's buffering capabilities. Therefore, we have developed methods to extend the buffering range by adding nitrilotriacetic acid (NTA) to solutions containing EGTA. This extension results from NTA having a lower K'(dCa) than EGTA. Such equilibria are solved by pCa Calculator, a computer program designed to aid in the study of Ca(2+)-dependent physiological processes while accounting for the effects of pH, temperature, and ionic strength. With multiple chelators and pH buffers from which to choose, pCa Calculator calculates the total concentration of each species required to achieve specified free concentrations of Ca(2+), ATP, and Mg(2+). The program is intuitive, user-friendly, and flexible enough to fix or vary the [Mg-ATP(2-)] and ionic strength. Moreover, it can account for increases in experimental volume from calcium addition. A comparative analysis is reported for testing solutions in the presence and absence of NTA by measuring the calcium binding affinity of fluorescent cardiac troponin C. These findings demonstrate that EGTA, when used in conjunction with NTA, improves and expands the regulation of free calcium in solution.     

A method of computing drug distribution in plasma using stepwise association constants: clofibrate acid as an illustrative example.

A program is described to compute the distribution of ligand among various complexes of ligand and protein using stepwise association constants. The key input data consists of stepwise association constants and the concentration of macromolecule and various concentrations of ligand. The key output data consists of the concentration of free drug and protein and the concentration of the various complexes.     

Metal binding to ligands: cadmium complexes with glutathione revisited

We studied the interaction of gamma-L-glutamyl-L-cysteinyl-glycine (glutathione, GSH) with cadmium ions (Cd(2+)) by first performing classical potentiometric pH titration measurements and then turning to additional spectroscopic methods. To estimate the residual concentrations of free cadmium, we studied the competition of glutathione with a Cd(2+)-sensitive dye, either an absorbing dye (murexide) or a fluorescent one (FluoZin-1), and consistent results were obtained with the two dyes. In KCl-containing Tes, Mops, or Tris buffer at pH 7.0 to 7.1 and 37 degrees C (and at a total Cd(2+) concentration of 0.01 mM), results suggest that free cadmium concentration is halved when the concentration of glutathione is approximately 0.05 mM; this mainly reflects the combined apparent dissociation constant for the Cd(glutathione) 1:1 complex under these conditions. To identify the other complexes formed, we used far-UV spectroscopy of the ligand-to-metal charge transfer absorption bands. The Cd(glutathione)(2) 1:2 complex predominated over the 1:1 complex only at high millimolar concentrations of total glutathione and not at low submillimolar concentrations of total glutathione. The apparent conditional constants derived from these spectroscopy results made it possible to discriminate between sets of absolute constants that would otherwise have simulated the pH titration data similarly well in this complicated system. Related experiments showed that although the Cl(-) ions in our media competed (modestly) with glutathione for binding to Cd(2+), the buffers we had chosen did not bind Cd(2+) significantly under our conditions. Our experiments also revealed that Cd(2+) may be adsorbed onto quartz or glass vessel walls, reducing the accuracy of theoretical predictions of the concentrations of species in solution. Lastly, the experiments confirmed the rapid kinetics of formation and dissociation of the UV-absorbing Cd(glutathione)(2) 1:2 complexes. The methods described here may be useful for biochemists needing to determine conditional binding constants for charge transfer metal-ligand complexes under their own conditions.     

Phosphodiesterolytic activity of lanthanide (III) complexes with α-amino acids

Kinetics of the hydrolysis of BNPP (bis(4-nitrophenyl)phosphate) mediated by lanthanide - samarium (III) and ytterbium (III) - alone and in the presence of various alfa amino acids has been systematically studied at 37.0 °C and I = 0.15 M in NaClO₄, in the pH interval of 7-9. The rate of BNPP cleavage is sensitive to metal ion concentration, pH, and ligand to metal molar ratio. Hydrolysis follows Michaelis-Menten-type saturation kinetics. For both metals, high pH values markedly increase the observed activity. Besides, potentiometric titrations of all these systems under identical conditions allowed us to identify the active coordination compounds towards hydrolysis. The results show that complexes with phosphodiesterolytic activity are monomeric cationic species such as [Ln(aa)₃(OH)]²⁺ or [Ln(aa)₂(OH)₂]⁺. Since phosphodiesterolytic activity is evident above pH 7 and it is increased with increasing pH, hydrolytic reactions of the metals are competitive processes that could lead to their precipitation as Ln(OH)₃(s). In this sense, ligand excess (for example, ligand to metal molar ratio equal to 30) was employed. Furthermore, due to its more extended hydrolysis, ytterbium shows, in general, less activity than samarium under the studied conditions. In general, a good phosphodiesterolytic activity is observed for these complexes under similar conditions to the physiological ones. Amino acids could be easily derivatized without changing their coordinating ability, leading to lanthanide complexes possibly capable of efficiently hydrolyzing the phosphodiester linkages of nucleic acids.     

The concentrations of free Mg2+ and free Zn2+ in equine blood plasma

The enzyme phosphoglucomutase can be used as a metal ion indicator to measure the concentrations of free Mg2+ and free Zn2+ in physiological fluids. In horse plasma, the concentration of free Mg2+ is close to 0.5 mM, whereas that of free Zn2+ is about 2 X 10(-10) M, although numerous physiological roles for Zn2+ have been postulated that would require free Zn2+ concentration orders of magnitude higher than this. A titration of plasma with Zn2+ shows that the fractional increase in free Zn2+ is essentially the same as the fractional increase in total exchangeable Zn2+, and the results are consistent with a model in which essentially all of the Zn2+ in plasma is bound to albumin. Regardless of the model, the buffering capacity of plasma for free Zn2+ is intrinsically low; however, its capacity relative to the total (exchangeable) Zn2+ present is maximal. The implications of this type of buffering for homeostasis of plasma Zn2+ are considered. Treatment of plasma with a strong reducing agent such as dithiothreitol (0.1 mM) substantially increases the apparent binding of Zn2+ and thus reduces the free Zn2+ concentration. However, the concentration of free Zn2+ appears to be insensitive to decreases in the physiological concentrations of reduced glutathione and cysteine. The concentrations of free Zn2+ and free Mg2+ in plasma are similar to those that have been reported for muscle tissue (rabbit). Their ratio is about 4 X 10(-7). The physiological implications of these concentrations are considered. In some cases, if the Zn2+ and Mg2+ complexes of an uncharacterized vertebrate protein exhibit significantly different properties, their relative importance under physiological conditions can be approximated by evaluating those of the mixed complexes present in a solution that contains the physiological concentration of free Mg2+, plus Zn2+ buffered with histidine, at the appropriate pH and ionic strength. Other metal ion/chelon systems that come close to reproducing the concentrations of free Mg2+ and free Zn2+ in horse plasma also are considered.     

Metal speciation in xylem and phloem exudates

Summary Two computer programs based on simultaneous chemical equilibria were compared for calculation of chemical species in xylem exudates. The first program, CHELATE, was developed to calculate the chemical species in xylem exudates while GEOCHEM was developed to calculate the speciation of natural aquatic systems. The output of the two programs should be similar since they are based on similar calculations. Data input to the programs consisted of concentration data for Ca, Cu, Fe, Mg, Mn, Zn, NH₄, PO₄, pH and 28 organic ligands reported for xylem exudates from soybean (Glycine max (L.) Merr.) and tomato (Lycopersicon esculentum Mill.) plants grown in nutrient solution¹⁷. The organic ligands included amino acids and low molecular weight organic acids (e.g., citric and malic). With the exception of Fe, there were large differences between CHELATE and GEOCHEM in the calculated speciation of nearly all metals in the xylem exudates. In general, there was better agreement between the programs for the speciation of alkaline earth metals than for transition metals. Discrepancies between the two programs were attributed to differences in 1) species considered and 2) stability constants. GEOCHEM considered a greater number of possible complexes. In addition, stability constants for some complexes differed by as much as 10 fold between the two programs. When the data base for GEOCHEM and CHELATE were the same, the output from CHELATE and GEOCHEM was almost identical. Thus, computations performed by the two programs are equally valid, but it is essential that the data base used in chemical models be verified before interpreting the output. Average concentration data for Al, Au, Ca, Cu, Fe, K, La, Mg, Mn, Na, Rb, Zn, Cl, MoO₄, PO₄, SO₄, HVO₄, pH and 18 organic ligands in phloem exudates from Yucca (Yucca flaccida Haw.) were complied from available literature and analyzed by GEOCHEM. Amino acids were the predominant organic ligand analyzed. Calculations revealed that alkali metals existed almost totally as the free ionic species (≥99%) whereas alkaline earth metals were transported as complexes with organic acids (oxalic, malic, and asparagine). Aluminum and Fe were present as hydroxyl species while <1% of micronutrients were transported as the free ion. Major micronutrient species were Cu-glutamine, Mn-asparagine and Zn-alanine. Information on calculated species present in phloem exudates could be useful to guide studies for isolation of metal-ligand complexes in phloem exudates.     

Vanadium(IV/V) speciation of pyridine-2,6-dicarboxylic acid and 4-hydroxy-pyridine-2,6-dicarboxylic acid complexes: potentiometry,EPR spectroscopy and comparison across oxidation states

Evaluation of stability of vanadium(IV) and (V) complexes under similar conditions is critical for the interpretation and assessment of bioactivity of various vanadium species. Detailed understanding of the chemical properties of these complexes is necessary to explain differences observed their activity in biological systems. These studies are carried out to link the chemistry of both vanadium(IV) and (V) complexes of two ligands, 2,6-pyridinedicarboxylic acid (dipicolinic acid, H(2)dipic) and 4-hydroxy-2,6-pyridinedicarboxylic acid (H(2)dipic-OH). Solution speciation of the two 2,6-pyridinedicarboxylic acids with vanadium(IV) and vanadium(V) ions was determined by pH-potentiometry at I=0.2 M (KCl) ionic strength and at T=298 K. The stability and the metal affinities of the ligands were compared. Vanadium(V) complexes were found to form only tridentate coordinated 1:1 complexes, while vanadium(IV) formed complexes with both 1:1 and 1:2 stoichiometries. The formation constant reflects hindered coordination of a second ligand molecule, presumably because of the relatively small size of the metal ion. The most probable binding mode of the complexes was further explored using ambient and low temperature EPR spectroscopy for vanadium(IV) and 51V NMR spectroscopy for vanadium(V) systems. Upon complex formation the pyridinol-OH in position 4 deprotonates with pK approximately 3.7-4.1, which is approximately 6 orders of magnitude lower than that of the free ligand. The deprotonation enhances the ligand metal ion affinity compared to the parent ligand dipicolinic acid. In the light of the speciation and stability data of the metal complexes, the efficiency of the two ligands in transporting the metal ion in the two different oxidation states are assessed and discussed.     

Determination of the total concentration and speciation of Al(III) in tea infusions

The aluminium species in different tea infusions were investigated, by determining their stability constants and concentration. This was done for some particular samples using a simple experimental method based on the sorption of aluminium on the strongly sorbing resin Chelex 100, by a batch procedure. From the thermodynamic information obtained it is possible to calculate the concentration of the different species, and in particular that of the free metal ion, which is very important for evaluating the adsorption of aluminium on biological membranes. It was found that aluminium in the tea infusions here considered is present at high total concentration, approximately 0.1 mM, but mainly linked to strong complexes, for instance with side reaction coefficient higher than 10(5.11) at pH 3.95 in one case (tea 1). This could be the reason for the low toxicity of aluminium in tea. These strong complexes were not dissociated even in the presence of Chelex 100. In this case only a limiting value of the reaction coefficient could be evaluated. The presence of the very strong complexes was found in all the tea sample here considered. In two of the considered samples (one black and one green tea) a part of Al(III) was linked to less strong complexes, for example with a reaction coefficient 10(4.14) (tea 2, pH 4.20). The presence in the considered tea infusions of other substances able to complex aluminium was also detected, by the well known ligand titration procedure, at concentration ranging from 0.65 to 3.37 mM in three tea infusions, and at somewhat higher concentration in the case of the ready drink, which was also considered for comparison.     

Metal binding to angiotensin converting enzyme: implications for the metal binding site

Angiotensin converting enzyme interacts with the chelator, 1,10-phenanthroline (OP) to form an OP-Zn-ACE ternary complex, which subsequently dissociates to OP-Zn and apoenzyme. The association and dissociation rate constants for the reaction OP + Zn-ACE in equilibrium OP-Zn-ACE have been determined and compared with those of known OP-metal complexes. Such constants were also used to calculate the rate constant for formation of the OP-Zn complex from OP-Zn-ACE. The rate of dissociation of zinc from ACE has been measured in the presence of EDTA (which acts only as a metal scavenger) as a function of chelator concentration, at different pH values, and with different buffers. The stability constant for the binding of zinc to apoACE log Kc = 8.2, determined by equilibrium dialysis using atomic absorption spectroscopy to assess metal concentration, is much smaller than that for Zn-carboxypeptidase A. Zn-thermolysin, or Zn-carbonic anhydrase. This weak binding is attributable to the zinc dissociation rate constant of ACE, 7.5 X 10(-3) sec-1 at pH 7.0, which is much greater than that of the other zinc metalloenzymes. These results lead to inferences regarding the metal binding site of ACE.