Next Article in Journal
Reservoir Control Operations and Water Resources Management
Previous Article in Journal
Optimizing Local Materials in Green Roofs Through Citizen Science Activities at a Primary School in Azores
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

An Evaluation of Metal Binding Constants to Cell Surface Receptors in Freshwater Organisms, and Their Application in Biotic Ligand Models to Predict Metal Toxicity

1
Geochem Australia, Kiama, NSW 2533, Australia
2
Aquatic Solutions International, North Arm Cove, NSW 2324, Australia
3
School of Natural Sciences, Macquarie University, Macquarie Park, NSW 2109, Australia
*
Author to whom correspondence should be addressed.
Water 2024, 16(20), 2999; https://doi.org/10.3390/w16202999
Submission received: 20 August 2024 / Revised: 27 September 2024 / Accepted: 8 October 2024 / Published: 21 October 2024
(This article belongs to the Special Issue Ecotoxicity of Pollutants on Aquatic Species)

Abstract

:
Biotic ligand models (BLMs) predict the toxicity of metals in aquatic environments by accounting for metal interactions with cell surface receptors (biotic ligands) in organisms, including water chemistry (metal speciation) and competing cations. Metal binding constants (log KMBL values), which indicate the affinity of metals for cell surface receptors, are fundamental to BLMs, but have only been reported for a few commonly investigated metals and freshwater species. This review evaluated literature toxicity and uptake data for seven key metals (cadmium (Cd), cobalt (Co), copper (Cu), nickel (Ni), lead (Pb), uranium (U), and zinc (Zn)) and four key competing cations (protons (H), calcium (Ca), magnesium (Mg), and sodium (Na)), to derive average metal binding constants for freshwater organisms/taxa. These constants will improve current BLMs for Cd, Cu, Ni, Pb, and Zn, and aid in developing new BLMs for Co and U. The derived metal binding constants accurately predicted metal toxicity for a wide range of freshwater organisms (75–88% of data were within a factor of two and 88–98% of data were within a factor of three of the ideal 1:1 agreement line), when considering metal speciation, competing cations and the fraction of cell receptors ((fC)M50%) occupied by the metal at the median (50%) effect concentration (EC50). For many organisms, toxicity occurs when 50% of cell surface receptors are occupied by the metal, though this threshold can vary. Some organisms exhibit toxicity with less than 50% receptor occupancy, while others with protective mechanisms show reduced toxicity, even with similar log KMBL values. For Cu, U, and Pb, the toxic effect of the metal hydroxide (as MOH+) must be considered in addition to the free metal ion (M2+), as these metals hydrolyse in circumneutral freshwaters (pH 5.5 to 8.5), contributing to toxicity.

1. Introduction

Metal toxicity in aquatic organisms occurs when metals bind to active functional groups (ligands) on the cell surface, such as cell membrane transporters (receptors), thereby impairing their vital functions (Figure 1). For instance, copper (Cu+) interferes with sodium (Na+) uptake, while cadmium (Cd2+), zinc (Zn2+), and lead (Pb2+) disrupt calcium (Ca2+) uptake [1]. Additionally, copper (Cu2+) can also disrupt Ca2+ uptake [2].
Metal toxicity to aquatic organisms is influenced largely by the following factors, which stem from three key conceptual models for metal–organism interaction—the free ion activity model (FIAM) [3,4], the biotic ligand model (BLM) [5,6], and the electrostatic toxicity model (ETM) [7,8]:
(i)
The metal’s binding affinity to ligands at cell surface receptors: a higher affinity means toxicity will occur at a lower metal concentration;
(ii)
Metal speciation in water: Free metal ions (Mz+) are typically toxic, but other positively charged species (MOH(z−1)+, MF(z−1)+) can also cause toxicity (Figure 1), especially for metals that hydrolyse within the pH range of natural freshwaters (pH 5.5–8.5), such as uranium (as uranyl, UO22+), copper (Cu2+), and lead (Pb2+). Predicting toxicity is challenging since total or dissolved metal concentrations are not reliable indicators of metal bioavailability [9]. Binding with anions, like bicarbonate or dissolved organic matter, reduces the concentration of free metal ions and reduces toxicity;
(iii)
Metals disrupt essential element uptake at cell surface receptors: however, the presence of essential elements and similar cations, such as Ca2+, magnesium (Mg2+), sodium (Na+), and protons (H+), in freshwater can ameliorate metal toxicity (Figure 1);
(iv)
Toxicity depends on the proportion of cell surface receptor sites bound by metals: If a toxic response occurs when 50% of sites are occupied, the TR50 can be calculated from the inverse of the metal concentration causing the effect. If less than 50% of sites are occupied at the TR50, then the metal’s binding affinity depends on both the metal concentration and the fraction of occupied sites.
The existing conceptual models (FIAM, BLM and ETM) that describe metal binding at cell surface receptors (biotic ligands) are widely acknowledged to share a common framework, with metal binding (stability) constants being a central component. However, metal binding constants have generally been reported for only a limited number of freshwater taxa (fish, crustaceans and micro-organisms), and often show considerable variability within and among species [10,11]. The lack of binding constants for metals such as cobalt (Co2+), UO22+ and Pb2+, and limited data for nickel (Ni2+), Cd2+, Cu2+ and Zn2+, constrains the wider use of these models in predicting metal toxicity across freshwater taxa.
The binding constants of metals with cell surface receptors (represented as log K values) can be measured directly through uptake (accumulation) studies or derived from toxicity data. When using toxicity data, the apparent binding constant (which is the inverse of the free metal ion concentration at the TR50) requires adjustment for ameliorative cation concentrations (H+, Na+, Mg2+, Ca2+) and receptor occupancy to obtain the true binding constant. If less than 50% of cell surface receptors are occupied at the TR50, increased toxicity is observed, but the binding constant decreases. Organisms with the ability to bind metals at other high-affinity, low-capacity, receptors can reduce metal toxicity.
This study evaluated the literature on the toxicity and uptake of seven key metals (Cd2+, Co2+, Cu2+, Ni2+, Pb2+, UO22+ and Zn2+) in freshwater systems. The primary objective was to compile a comprehensive database of metal binding (stability) constants at cell surface receptors (biotic ligands) for a diverse range of freshwater taxa, including fish (e.g., Oncorhynchus mykiss and Pimephales promelas), crustaceans (e.g., Daphnia magna and Ceriodaphnia dubia), microalgae (e.g., Raphidocelis subcapitata) and molluscs (e.g., Lymnaea stagnalis). Binding constants for the ameliorative cations Ca2+, Mg2+, Na+ and H+, were also included. Potassium (K+) was excluded due to its relatively minor ameliorative effect. A secondary objective was to compare the measured and predicted toxicity (EC/LC50 values) of the seven metals, as a key measure of the reliability of the derived metal binding constants.
The following section outlines the development and unification of the existing conceptual models, providing the foundation for deriving metal binding constants and predicting metal toxicity to aquatic organisms.

2. Conceptual Models of Metal Binding to Biotic Ligands

To cause a toxic response, a metal must first interact with transporters (receptors) on the cell membrane surface, comprising macromolecules dominated by oxygen-containing functional (donor) groups, such as carboxylate (COOH) [12]. These functional groups have an average proton dissociation constant between 4 and 6 [11] and, as such, in most natural waters (circumneutral pH), the cell surface receptors will be deprotonated and consist of negatively charged surface sites, {X–cell}. The transport and binding of soluble metals at cell surface receptors is considered to occur rapidly [12]. The complexation of the metal with the cell surface results in the formation of {M–X–cell}. The seven divalent metals investigated are more likely to bind with oxygen-rich donor groups, like carboxylate (see Section 3.13). In contrast, toxic metals like Ag+ and Hg2+ show a preference for nitrogen or sulphur-rich donor groups over oxygen.
In the development of the extended FIAM, Brown and Markich [4] demonstrated that the toxic response (TR) of an organism to a metal could be described by Equation (1)
TR = TR max X cell T n M n X cell E n ( K 1 + [ M ] ) n + X cell T n [ M ] n
where TRmax is the maximum toxic response, {X–cell}T is the total number of receptor sites on the cell membrane, {X–cell}E is the concentration of receptor sites bound by metal ions which causes half of the maximal toxic response (TR50), K−1 is the inverse of the binding (or stability) constant (which will be defined as KMBL, so K1 = 1/K−1 = KMBL) between the metal and the receptor sites on the cell membrane surface and n determines the slope of the concentration-response curve.
Equation (1) is the product of two processes: (1) the reaction of the metal with receptor sites on the cell membrane surface and (2) the response of the organism to the presence of the toxic metal. In the first process, an explicit equation can be derived for the concentration of the total number of surface receptor sites bound by the metal (i.e., {M–X–cell}), using the law of mass action [4], as described by Equation (2).
{ M X cell } = { X cell } T [ M ] K 1 + [ M ]
Similarly, the response of the organism to the interaction of the metal with receptor sites on the cell membrane surface is described by Equation (3).
TR = TR max M X cell n X cell E n + { M X cell } n
The substitution of Equation (2) into (3) leads to Equation (1). Equation (3), through the exponent n, accounts for the fact that concentration-response curves are generally sigmoidal in shape and have large slopes (i.e., n ≥ 5) [4,13].
A closer inspection of Equation (1) is useful as it contains important information. Moving TRmax to the left-hand side of Equation (1), inverting followed by simplification leads to Equation (4):
TR max TR = 1 f TR = X cell E n K 1 + M n X cell T n M n + 1
where fTR is the ratio of the toxic response to the maximum response (e.g., TR50/TRmax = f50 = 0.5). The rearrangement of Equation (4) gives
1 f TR 1 X cell T n [ M ] n = X cell E n ( K 1 + [ M ] ) n
Equation (4) can be simplified, by taking the n-th root, and following rearrangement gives
1 f TR 1 1 / n { X cell } T X cell E = ( K 1 + [ M ] ) [ M ] = 1 + K 1 M K 1 M
Further rearrangement leads to Equation (7),
[ M ] = K 1 1 / f TR 1 1 / n ( 1 / f C M 50 % ) 1
where (fC)M50% (= {Xcell}E/{Xcell}T) is the fraction of receptor sites occupied at the TR50. If the TR50 (i.e., fTR = 0.5) is considered, then Equation (7) can be reduced to
[ M ] = f C M 50 % K 1 1 f C M 50 %
In the absence of competing ions, such as protons (H+), Na+ and the hardness cations, Ca2+ and Mg2+, Equation (8) is identical to that given by De Schamphelaere and Janssen [9] in the development of a BLM to predict the effects of such ions (hardness and protons) on acute copper (Cu2+) toxicity. Additionally, if the TR50 occurs when half of the receptor sites are occupied, then (fC)M50% is equivalent to (1 − (fC)M50%) [i.e., (fC)M50% = 0.5] and the metal concentration at the TR50 is given by
[M] = K−1 = 1/K1 = 1/KMBL
Alternatively, if less than half of the receptor sites are occupied at the TR50, then the fraction ((fC)M50%/(1 − (fC)M50%) in Equation (8) is less than unity since (fC)M50% < 0.5 and, as such, K−1 is larger than the metal concentration which causes the TR50. It is not possible for more than 50% of the receptor cells to be occupied at the TR50. To demonstrate these behaviours, Figure S1 shows the calculated TR when the following parameter values are used in Equation (1): TRmax = 100; K−1 = 0.5; n = 10; {Xcell}T = 1 and (fC)M50% = 0.2, 0.5 and 0.8 (units are not provided for brevity and the values used are for simplicity and the purpose of the demonstration). It is clear from Figure S1 that it is not possible to achieve the true TRmax if more than half the receptor sites are occupied at the TR50, thereby invalidating such scenarios. Consequently, the minimum value K−1 can have is the metal concentration that occurs at the TR50 and this would happen when the biological response causing the TR50 occurs when 50% of the receptor sites are occupied. Conversely, when (fC)M50% ≤ 0.5, the true TRmax is attained; however, as (fC)M50% continues to decrease below 0.5, the TR50 moves progressively to lower metal concentrations, as indicated by Equation (8).
Brown and Markich [4] showed that the FIAM could be adapted to systems where metal ion binding to surface sites is ameliorated by protons (H+). In this scenario, the expression for {M–X–cell} is modified to
M X cell = K MBL X cell T M 1 + K MBL M + K HBL H
(i.e., an extension of Equation (2)), where KMBL has replaced K1 and KHBL is the binding constant of H+ with receptor sites on the cell membrane surface. The extension of the sequence of Equation (2) to Equation (10) above, in the case of amelioration by H+, leads to the metal concentration (free ion) at the TR50:
[ M ] = f C M 50 % 1 f C M 50 % 1 K MBL 1 + K HBL H
The expansion of Equation (11) to include amelioration by Na+, and the hardness ions, Ca2+ and Mg2+, gives (the metal concentration at the TR50)
[ M ] = f C M 50 % 1 f C M 50 % 1 K MBL 1 + K HBL H + K NaBL Na + K CaBL Ca + K MgBL Mg
where KNaBL, KCaBL and KMgBL are the binding constants of Na, Ca and Mg cations, respectively, with receptor sites on the cell membrane surface (Figure 1). Equation (12) is identical to that given by De Schamphelaere and Janssen [9] in the development of a BLM to predict the effects of such ions on acute metal toxicity. Most studies have found that the potassium ion (K+) does not provide an ameliorative affect to metal toxicity and, as such, is not considered further.
It has frequently been found that toxicity is not only caused by the free metal ion, but other cationic metal–ligand species, ML (such as metal hydroxide (MOH) or fluoride (MF)), can also elicit a toxic response. Brown and Markich [4] extended the FIAM to account for the toxicity caused by such species. Again, the model equations extend Equation (2):
M X cell + LM X cell = { X cell } T K 1 + K 3 K 4 L [ M ] 1 + K 1 + K 3 K 4 L [ M ]
where K1 is again used for KMBL, K3 is the binding constant of ML with receptor sites on the cell membrane surface and K4 is the binding constant of the metal (M) and the ligand (L). The rearrangement of Equation (13) and combination of it with a version of Equation (3) that contains both {M–X–cell} and {LM–X–cell} gives the following:
K MBL M + K LMBL ML = f C M 50 % 1 f C M 50 % 1 + K HBL H + K NaBL Na + K CaBL Ca + K MgBL Mg
Rearrangement of Equation (14), leads to:
M + K LMBL K MBL ML = f C M 50 % 1 f C M 50 % 1 K MBL 1 + K HBL H + K NaBL Na + K CaBL Ca + K MgBL Mg
Markich et al. [14,15] demonstrated that the toxic response of bivalves to dioxouranium(II) or copper(II) was best described by combining the activities of [M2+] and [MOH+], where [M2+] had twice the effect of [MOH+] (i.e., [M2+] + [MOH+]/2). Reanalysis of Cu(II) toxicity data for D. magna from De Schamphalaere and Janssen [9] revealed a similar relationship (i.e., [Cu2+] + [CuOH+]/2) (refer to Supplementary Information). Crémazy et al. [16] reported that scandium (Sc(III)) uptake by the microalga Chlamydomonas reinhardtii was best explained by a combination of Sc3+ and its hydroxide complexes (such as ScOH2+ and Sc(OH)2+) above pH 6.5. Wilkinson et al. [17] found that both [M3+] and [MF2+] were needed to describe the toxicity of aluminium(III) to Atlantic salmon (O. mykiss), with the activity of [Al3+] being 1.5 times that of [AlF2+] [4]. Similarly, Crémazy et al. [18] reported that both [Sc3+] and [ScF2+] were necessary for explaining Sc(III) uptake in C. reinhardtii, with binding constants having a ratio of 1.9 (see Equation (15)). It can be seen from these examples that the toxic response induced by [ML] is related to that of [M] by dividing [ML] by the ratio of the charges of the two species. This behaviour then allows the metal concentration at the TR50 to be described by the following:
M + Σ z ML z M ML = f C M 50 % 1 f C M 50 % 1 K MBL 1 + K HBL H + K NaBL Na + K CaBL Ca + K MgBL Mg
where zM is the charge of the metal ion, and zML is the charge of each ML species that also causes a toxic response (the summation is used to define the number of ML species that may cause a toxic response). A comparison of Equations (15) and (16) indicates that KLMBL must be less than KMBL because zML/zM is less than one (i.e., for M2+ and MOH+, KHOMBL/KMBL = 0.5). The use of charge alone in Equation (16) might be overly simplistic. Binding constants generally consist of two components: an intrinsic term and an electrostatic term. This approach has been successfully used to explain metal binding to natural dissolved organic matter (DOM), such as humic and fulvic acids [19]. The surfaces of natural DOM resemble those of cell membrane surfaces. Equation (16) may imply that the intrinsic terms for M and ML are similar (see Section 7 for more detail).
The model development outlined above demonstrates that the formulation of the BLM can be derived directly from rearrangement of the formulation of the FIAM, and both can be equally utilised to describe the toxic response of organisms to metals. The models show that toxic response is a function of the strength of binding that metals have at organism cell surface receptors, a response that can be ameliorated by the presence of several cations in aqueous systems, such as, H+, Na+, Ca2+ and Mg2+. The concentration at which a toxic response occurs is also related to the fraction of sites on the cell surface that are bound by the metal at the TR50; the lower the fraction of sites on the cell surface that are bound, the lower the metal concentration that elicits the TR50. It is also important to understand that the toxic response is caused by the free metal ion as well as other cationic species of the metal when bound to various ligands (e.g., MOH and MF).
Metals bound to other ligands (such as sulphate, chloride, carbonate, phosphate or organic macromolecules) in solution do not typically cause a toxic response and, as such, the total metal concentration at the TR50 will generally be greater than the concentration of metal species that cause the toxic response (e.g., M, MOH, MF, etc.). The combination of the two conceptual models also demonstrates how concentration-response curves are frequently sigmoidal in shape and have large slopes (i.e., n ≥ 5). Importantly, the extended FIAM demonstrates a cascade of effects that includes (a) the reaction of a metal with receptor sites on the cell membrane surface, (b) the response of an organism to the presence of the toxic metal and (c) the magnification of toxicity depending on the fraction of cell surface sites occupied at the TR50 (Figure 2).
In addition to the BLM and FIAM, ETMs have been used to describe metal and metalloid uptake and toxicity to plants [7,8]. These models describe an electrical double layer above the cell membrane surface with three key electrical properties of the cell membrane [20]. There is a transmembrane electrical potential (Em) that describes the electrical potential difference between the cell interior and the bulk phase medium. The second is a negative electrical potential (ψ0) on the cell membrane surface, due to acidic amino acids in membrane proteins and the phosphate groups of membrane phospholipids [20]. Finally, there is the electrical potential difference (Em,surf) from surface to surface, which is the driving force of metal transport across the cell membrane. Classical electrostatic theory is used together with ion binding through definition of the cell membrane surface charge density (σ) and the electrical potential of the exterior surface of the cell membrane (ψ0). The model provides an elegant depiction of binding at cell surfaces that can be utilised to describe the toxicity and uptake of metals in plants and that could equally be utilised to describe the toxicity of metals to aquatic organisms.
An important difference between ETMs and BLMs (including the FIAM) is that the former describes two distinct binding sites, with metal binding defined by the following two reactions:
M z + + X cell   M X cell ( z 1 ) +
M z + + Y cell   M Y cell z +
with Equation (17) representing a reaction with a negatively charged site (as described above) and Equation (18) defining the reaction of a metal with an uncharged site. These sites likely describe carboxylic acid (-COOH) binding and amine (-NH2) binding, respectively. As such, the binding of an amine group will occur at a greater pH. In some studies of aquatic organisms, the toxicity of a metal has been found to have an inflexion point as the pH is increased. Such an inflexion point may be the result of the binding of metals to a second cell surface site at the increased pH (see Section 8).

3. Metal Binding Constants

This section provides the binding (or stability) constants (log KMBL values) for seven environmentally important metals (Zn2+, Ni2+, Cd2+, Co2+, Cu2+, UO22+ and Pb2+; Table 1) and four key ameliorative cations (H+, Na+, Ca2+ and Mg2+; Table 2) with cell surface receptors (biotic ligands) in freshwater organisms. For Zn, Ni, Cd and Co, it is assumed that only the free metal ion (M2+) induces a toxic response in circumneutral (pH 5.5 to 8.5) oxic freshwaters. However, for Cu, UO2 and Pb, it is assumed that both the free metal ion (M2+) and the cationic hydrolysed species, MOH+, induces a toxic response. The Supplementary Materials include a summary of each study used to derive the log KMBL values for the seven divalent metals (Table S1) and four ameliorative cations (Table S2). A summary of the relative binding affinities of these metals/cations to biotic ligands, along with comparisons between different freshwater taxa and key organisms, is provided in Section 3.12.

3.1. Methodology

A comprehensive literature review on the toxicity and uptake of Cd2+, Co2+, Cu2+, Ni2+, Pb2+, UO22+ and Zn2+ (and Ca2+, Mg2+, Na+ and H+) in freshwater organisms was conducted using the Web of Science, Scopus, Google Scholar and ECOTOX databases (up to July 2024). For toxicity studies, the median effect concentrations (TR50) were assessed using ecologically relevant endpoints such as survival (LC50), reproduction and growth (EC50), since TR50 values are necessary for deriving log KMBL values (see Section 2). Both acute and chronic exposures, as defined by Warne et al. [23], were considered. Only studies reporting the concentrations of (i) key complexing ligands (bicarbonate (alkalinity), sulphate, chloride, phosphate (microalgae/macrophytes) and/or DOC), and (ii) ameliorative cations (protons (pH), Ca2+, Mg2+ (or hardness) and Na+) in test waters (natural or reconstituted) were considered, as these are necessary for determining metal speciation (e.g., free metal ions) and competition with test metals at cell surface biotic ligands, respectively. These studies also needed to demonstrate a robust experimental design (e.g., sufficient controls, test concentrations, replicates and control acceptability criteria) and data analysis (established statistical tests to derive EC/LC50 values), with quality scores ≥ 75% (see Warne et al. [23]). Preference was given to studies that measured test metal concentrations, and data were excluded if the solubility of the test metal was exceeded (and concentrations were not measured).
The speciation, or chemical form, of Cd, Co, Cu, Ni, Pb, UO2 and Zn (and Na, Ca and Mg) in test waters was predicted using the WHAM (v7.05) chemical speciation code [24]. Inorganic equilibrium constants were taken from Markich and Brown [25] and corrected for ionic strength using the Davies equation. Metal binding to natural DOM was estimated using the Humic Ion Binding Model (HIBM) VII (a component of WHAM), assuming DOM contained 50% carbon with 65% active cation binding sites [26], represented as fulvic acid (FA). For example, a DOC concentration of 1 mg/L was modelled as 1.3 mg/L FA. The model also assumed that iron(III) and aluminium(III) activities were controlled by hydroxide solubility. Adjustments to the default FA binding constants (log KMA and ∆LK2) for Cd, Co, Cu, Ni, Pb, UO2 and Zn were made to improve the accuracy of the free metal ion (M2+) activity predictions in natural freshwaters, which HIBM VII tends to underestimate for Cd, Cu, Pb and UO2, or overestimate for Ni, Co and Zn [27,28,29,30,31,32,33]. These adjustments, based on previous research, aimed to improve the accuracy of metal binding with natural DOM and provide better estimates of M2+ and/or MOH+ needed to derive log KMBL values. The modelling incorporated physicochemical parameters from the test waters as input data.
For toxicity experiments, log KMBL values were derived by adjusting TR50 values according to the percentages of M2+ and/or MOH+ (from speciation modelling) and accounting for the competitive (or ameliorative) effects of H+, Ca2+, Mg2+ and Na+, as outlined in Section 2. For metal uptake experiments, log KMBL values were derived using the inverse of the Michaelis–Menton constant (1/Km), as outlined in Section 2). Metal uptake is usually described by Michaelis–Menton kinetics (i.e., fast Langmuirian adsorption of M2+ and/or MOH+ at cell surface receptor sites), followed by a first-order rate-limiting internalisation step [34]. It is essential to include only intracellular metal (i.e., metal taken up via the cell surface receptor sites/transporters; Figure 1), ensuring that surface-bound (extracellular) metal is excluded (often using EDTA washes; see Hassler et al. [35] for more information on chemical extractions).
When log KMBL values were derived from both uptake and toxicity experiments under the same test conditions, (fC)M50% values (the fraction of receptor sites occupied at the TR50; Section 2) were directly determined (see Section 6). The (fC)M50% values were also derived using the methodology described by De Schamphalaere and Janssen [9]. Frequently, it was necessary to use a derived (fC)M50% value across multiple studies for the same metal/organism, or for multiple metals with the same organism.
A one-way analysis of variance (ANOVA) was used to determine if average log KMBL values differed significantly (p ≤ 0.05) among (i) taxa (e.g., fish vs. crustaceans vs. molluscs), (ii) key individual species (e.g., C. dubia vs. D. magna vs. O. mykiss vs. P. promelas vs. R. subcapitata vs. Chlorella sp.) and (iii) metals (e.g., Cd2+ vs. Co2+ vs. Cu2+ vs. Ni2+ vs. Pb2+ vs. UO22+ vs. Zn2+). The Tukey Honest Significant Difference test was used to compare average log KMBL values where ANOVA indicated significant (p ≤ 0.05) differences. ANOVA assumptions [36] were tested and model adequacy was confirmed using either raw or transformed (log10) data. The overall consistency of log KMBL values across uptake and toxicity studies allowed for the combination of values for each metal (Table S3). Both acute and chronic toxicity data (as defined by Warne et al. [23]) were employed to derive metal binding constants (log KMBL values)—see Section 5 for further discussion.
To simplify terminology, the unified conceptual model of metal–organism interaction (Section 2) will hereafter be referred to as the BLM. The capacity of the BLM to predict metal toxicity (EC/LC50) was assessed using the factor of two rule, which allows for a range of four around the ideal 1:1 agreement line in comparisons between measured and predicted toxicity. Additionally, a factor of three assessment was considered, in line with analyses by Peters et al. [37] and Price et al. [38]. This assessment method has shown broad applicability to organisms across different freshwater taxa, including fish, crustaceans and microalgae [38,39].
For trace metals, binding constants (log KMBL) can be determined directly, but for ameliorative cations, these constants are generally derived through competition experiments. In these experiments, the metal binding to cell surface receptors varies as a function of the ameliorative cation concentration. As the concentration of an ameliorative cation increases (or pH decreases), the metal concentration required to reach the TR50 rises. Binding constants for ameliorative cations are derived from the change in TR50 with respect to ameliorative cation concentrations. For Mg, binding constants were also derived directly from toxicity tests, corroborating the findings from competition experiments.
It might be expected that the binding constants for ameliorative cations would differ depending on the metal being examined (e.g., the binding constant for Ca2+ may differ between experiments with Cd2+ and Pb2+). However, the binding constants for each of the four ameliorative cations (Ca2+, Mg2+, H+ and Na+) are consistent across all metals (see Section 3.9, Section 3.10, Section 3.11 and Section 3.12), indicating that the biotic ligand binding these metals is either the same or very similar for all seven metals examined.

3.2. Zinc

Table 1 provides the Zn2+ binding constants for 42 freshwater organisms from eight taxa, with an average log KZnBL value of 6.0 ± 0.1 (95% confidence limit). Additional information is provided in Table S1 and the Supplementary Materials. The log KZnBL values were uniformly distributed across molluscs (12 species), fish (10 species) and crustaceans (10 species), followed by the rotifers (five species) and microalgae (two species). The average log KZnBL values were not significantly (p > 0.05) different across molluscs (6.0 ± 0.2), fish (5.9 ± 0.2), crustaceans (6.0 ± 0.2), rotifers (6.0 ± 0.4) and microalgae (6.0 ± 0.1) (Table 1). The average log KZnBL values for the cladocerans C. dubia (5.9 ± 0.2) and D. magna (5.8 ± 0.3), the two most studied crustaceans, were comparable. For the two most studied fish species, the average log KZnBL value for rainbow trout O. mykiss (5.9 ± 0.1) was consistent with that of the fathead minnow P. promelas (5.7 ± 0.4). Likewise, the average log KZnBL values were similar for the two most-studied microalgae, R. subcapitata (6.1 ± 0.1) and Chlorella sp. (5.8 ± 0.1) (Table 1). The average log KZnBL value of 6.0 ± 0.1 (all organisms) is higher than the reported range of 5.3–5.6 used for fish and crustaceans in Zn BLMs [10], but consistent with the average log KZnBL value of 5.8 derived directly from uptake experiments (Table S3).

3.3. Nickel

Table 1 provides the Ni2+ binding constants for 31 freshwater organisms from seven taxa, with an average log KNiBL value of 4.8 ± 0.1 (95% confidence limit). Additional information is provided in Table S1 and the Supplementary Materials. Most log KNiBL values were determined for crustaceans (16 species), followed by molluscs (seven species), fish (three species) and microalgae (three species). The average log KNiBL values were not significantly (p > 0.05) different across crustaceans (4.9 ± 0.1), molluscs (4.9 ± 0.1), microalgae (4.8 ± 0.2) and fish (4.7 ± 0.2) (Table 1). Among molluscs, the average log KNiBL values for bivalves (4.9 ± 0.1) and gastropods (4.9 ± 0.5) were identical. Similarly, among fish, the average log KNiBL value for O. mykiss (4.7 ± 0.5) was identical to P. promelas (4.7 ± 0.7) and Neogobius melanostomus (4.7). Although the average log KNiBL value for C. dubia (5.2 ± 0.5) was not significantly (p > 0.05) different to D. magna (4.9 ± 0.3) (Table 1), the 0.3 log unit difference between the two cladoceran species may stem from an overestimation in two of the four derived log KNiBL values for C. dubia, with values of 5.4 [40] and 5.5 [37] being higher than 5.0 [41] and 4.9 [42] (Table S1). This is likely due to the use of (fC)Ni50% values that were too high (see Supplementary Materials for further details).
The average log KNiBL value of 4.8 (all organisms) exceeds the value of 4.0 used by Santore et al. [43] in a recent analysis of the Ni BLM for a range of freshwater taxa. The authors also report log KMBL values for Ca (4.3), Mg (3.6), Na (1.0) and H (4.7), which are inconsistent with those provided in the present study. They also report a log KNiOHBL value of −5.5, which is irrelevant for circumneutral freshwaters (pH 5.5–8.5), as Ni does not hydrolyse until pH 9.5 (~10%). It seems that the toxicity data were used to “re-calibrate” the log KMBL values. However, the average log KNiBL value of 4.8 from the present study is comparable to the value of 4.9 used for crustaceans in the Ni BLM of Kozlova et al. [44], and in close agreement with the average log KNiBL value of 4.7 derived directly from uptake experiments (Table S3).

3.4. Cadmium

Table 1 provides the Cd2+ binding constants for 37 freshwater organisms from seven taxa, with an average log KCdBL value of 7.0 ± 0.1 (95% confidence limit). Additional information is provided in Table S1 and the Supplementary Materials. The majority of log KCdBL values were determined for crustaceans (13 species), followed by molluscs (11 species) and fish (9 species). The average log KCdBL value for fish (7.3 ± 0.1) was significantly (p ≤ 0.05) higher than that of crustaceans (6.9 ± 0.2), molluscs (6.8 ± 0.2) and microalgae (7.0 ± 0.4) (Table 1). Among fish, the average log KCdBL value for O. mykiss (7.5 ± 0.1) was higher than that of P. promelas (7.1 ± 0.2). In contrast, the average log KCdBL values for the cladocerans C. dubia (7.0 ± 0.9) and D. magna (6.9 ± 0.3) were comparable (Table 1). Likewise, the average log KCdBL values for bivalves (6.9 ± 0.3) and gastropods (6.8 ± 0.6) were similar. Although fish had a higher average log KCdBL value compared to other taxa, their ranges still overlapped (fish: 7.0–7.6, crustaceans: 6.5–7.4, molluscs: 6.5–7.5). The average log KCdBL value of 7.0 is lower than the range of 7.3–7.6 used for fish in Cd BLMs [10], but comparable to the value of 7.0 used for crustaceans in the Cd BLM of Clifford and McGeer [45], and consistent with the average log KCdBL value of 7.2 derived directly from uptake studies (Table S3).

3.5. Cobalt

Table 1 provides the Co2+ binding constants for 15 freshwater organisms from six taxa, with an average log KCoBL value of 5.8 ± 0.2 (95% confidence limit). Additional information is provided in Table S1 and the Supplementary Materials. Most log KCoBL values were determined for molluscs (seven species), followed by crustaceans (three species) and fish (two species). The average log KCoBL values for crustaceans (5.6 ± 0.5), fish (5.5 ± 0.4) and microalgae (5.7 ± 0.1) were similar, but all were significantly (p ≤ 0.05) lower than the average log KCoBL value for molluscs (6.0 ± 0.1) (Table 1). Among molluscs, the average log KCoBL values for bivalves (6.0 ± 0.1) and gastropods (5.9) were comparable. Cobalt has the least amount of data among the seven metals, with most of it derived from only two studies [46,47] with single values per organism, making taxa comparisons less reliable. At present, no Co BLM exists for aquatic organisms.

3.6. Copper

Table 1 provides the Cu2+ binding constants for 87 freshwater organisms from nine taxa, with an average log KCuBL value of 7.9 ± 0.1 (95% confidence limit). Additional information is provided in Table S1 and the Supplementary Materials. Most log KCuBL values were determined for crustaceans (32 species) and molluscs (32 species), followed by fish (ten species), microalgae (three species) and rotifers (two species). The average log KCuBL values were not significantly (p > 0.05) different across crustaceans (7.9 ± 0.1), molluscs (7.9 ± 0.1), fish (7.8 ± 0.2), microalgae (7.9 ± 0.1) and rotifers (8.0 ± 0.1) (Table 1). The average log KCuBL value for C. dubia (8.0 ± 0.2) was similar to that of D. magna (7.8 ± 0.1), and the cladoceran group (7.9 ± 0.1) (Table 1). Among molluscs, the average log KCuBL values for bivalves (7.8 ± 0.1) and gastropods (8.0 ± 0.1) were comparable. The average log KCuBL values were identical for the two microalgae, Chlorella sp. (7.8 ± 0.1) and R. subcapitata (7.8 ± 0.2). The average log KCuBL value for fish (7.8 ± 0.2) largely reflects the difference between O. mykiss (7.9 ± 0.1) and P. promelas (7.6 ± 0.1). The average log KCuBL value of 7.9 is at the upper end of the reported range of 7.4–8.0 used for fish and crustaceans in Cu BLMs [10], but is in excellent agreement with the average log KCuBL value of 7.9 derived directly from uptake studies (Table S3).

3.7. Uranium

Table 1 provides the UO22+ (uranyl) binding constants for 20 freshwater organisms from eight taxa, with an average log KUO2BL value of 7.7 ± 0.1 (95% confidence limit). Additional information is provided in Table S1 and the Supplementary Materials. Most log KUO2BL values were determined for crustaceans (six species), followed by fish, microalgae and macrophytes (three species) and molluscs (two species). The average log KUO2BL values were not significantly (p > 0.05) different across crustaceans (7.7 ± 0.4), fish (7.4 ± 0.1), microalgae (7.9 ± 0.1), macrophytes (7.6 ± 0.2) and molluscs (7.9 ± 0.9) (Table 1). Most organisms are represented by a single log KUO2BL value, hindering sensible comparisons. The average log KUO2BL value of 7.7 is in excellent agreement with the average log KUO2BL value of 7.7 derived directly from uptake studies (Table S3). At present, no UO2 BLM exists for aquatic organisms.

3.8. Lead

Table 1 provides the Pb2+ binding constants for 21 freshwater organisms from seven taxa, with an average log KPbBL value of 6.8 ± 0.1 (95% confidence limit). Additional information is provided in Table S1 and the Supplementary Materials. Most log KPbBL values were determined for molluscs (eight species), followed by crustaceans (five species), fish (three species), rotifers (two species) and microalgae (one species). The average log KPbBL value for molluscs (6.6 ± 0.1) was significantly (p ≤ 0.05) lower than that of crustaceans (6.9 ± 0.2), fish (7.1 ± 0.2) and microalgae (7.0 ± 0.2) (Table 1).
It should be noted that six of the eight data points for mollucs were from one study (average log KPbBL of 6.5) with unionid bivalves [46], which skewed the overall average downward, particularly when compared to the average (6.8) of the other two other unionid species. The average log KPbBL value for O. mykiss (7.1 ± 0.6) was similar than that of P. promelas (6.8 ± 0.2) (Table 1). The average log KPbBL value for C. dubia (7.0 ± 0.4) was comparable with that of D. magna (6.7 ± 0.3). More toxicity data/log KPbBL values are required for molluscs to re-assess comparisons with other taxa. The average log KPbBL value of 6.8 is generally consistent with the value of 6.7 used for a range of freshwater taxa in the Pb BLM of DeForest et al. [48], and the average log KPbBL value of 7.0 derived directly from uptake studies (Table S3).

3.9. Calcium

Calcium binding constants (log KCaBL) can be derived from experimental studies where freshwater organisms exposed to metals at varying Ca2+ concentrations show reduced metal toxicity as the Ca2+ concentration increases. Table 2 provides the Ca2+ binding constants for 14 freshwater organisms from six taxa, with an average log KCaBL value of 3.4 ± 0.1 (95% confidence limit). Additional information is provided in Table S2 and the Supplementary Materials. Most log KCaBL values were determined for crustaceans (five species), followed by fish (four species) and microalgae (two species). The average log KCaBL values were not significantly (p > 0.05) different among crustaceans (3.3 ± 0.1), fish (3.5 ± 0.1) and microalgae (3.4 ± 1.4) (Table 2). These values align with those currently used in BLMs for invertebrates (3.5) and fish (3.4) [49].
The average log KCaBL values from metal-specific studies were found to be 3.3 ± 0.2 for Zn2+, 3.3 ± 0.2 for Ni2+, 3.5 ± 0.1 for Cd2+, 3.3 ± 0.2 for Cu2+, 3.2 for UO22+ and 3.2 for Pb2+, with no data available for Co2+ (uncertainties apply where studies reported five or more values). These averages align with the overall value of 3.4 ± 0.1 (Table 2), with no significant differences observed (p > 0.05) across metals, indicating metal binding occurs at the same or similar receptor sites on the cell surface.

3.10. Magnesium

Magnesium binding constants (log KMgBL) can be derived from experimental studies where freshwater organisms exposed to metals at varying Mg2+ concentrations show reduced metal toxicity as the Mg2+ concentration increases. Table 2 provides the Mg2+ binding constants for 12 freshwater organisms from six taxa, with an average log KMgBL value of 2.8 ± 0.1 (95% confidence limit). Additional information is provided in Table S2 and the Supplementary Materials. Most log KMgBL values were determined for crustaceans (five species), followed by fish and microalgae (two species). Two of the Mg2+ binding constants were determined directly from toxicity studies [50,51] rather than from reduced metal toxicity with increasing Mg2+ concentration. The EC/LC50 values from the two studies were corrected for Mg speciation (e.g., fraction of Mg occurring as Mg2+) and further adjusted for the concentrations of the other ameliorative cations (Ca2+, Na+ and H+). The log KMgBL values obtained were consistent with those derived from the ameliorative effect of Mg (see Table S2). The average log KMgBL values were not significantly (p > 0.05) different among crustaceans (2.8 ± 0.1), fish (2.9 ± 0.3) and microalgae (2.8 ± 0.3) (Table 2). These values are lower than those currently used in BLMs for invertebrates (3.5) and fish (4.0) [49].
The average log KMgBL values from metal-specific studies were found to be 2.9 ± 0.2 (Zn2+), 2.8 ± 0.2 (Ni2+), 2.8 ± 0.2 (Cd2+), 2.8 (Cu2+) and 2.7 (UO22+), with no data available for Co2+ or Pb2+ (uncertainties apply where studies reported five or more values). These averages align with the overall value of 2.8 ± 0.1 (Table 2), with no significant differences observed (p > 0.05) across metals, indicating metal binding occurs at the same or similar receptor sites on the cell surface. The two log KMgBL values obtained from toxicity studies averaged 2.7, which matches the group average.

3.11. Sodium

Sodium binding constants (log KNaBL) can be derived from experimental studies where freshwater organisms exposed to metals at varying Na+ concentrations show reduced metal toxicity as the Na+ concentration increases. Table 2 provides the Na+ binding constants for 11 freshwater organisms from three taxa, with an average log KNaBL value of 2.2 ± 0.1 (95% confidence limit). Additional information is provided in Table S2 and the Supplementary Materials. Most log KNaBL values were determined for crustaceans (nine species), followed by fish and microlagae (one species). The average log KNaBL values were not significantly (p > 0.05) different among crustaceans (2.2 ± 0.1), fish (2.4) and microalgae (2.2) (Table 2). These values are lower than those currently used in BLMs for invertebrates (2.7) and fish (3.0) [49]. The binding capacities of Ca2+, Mg2+ and Na+ to organic carboxylate ligands [25] suggest that these metals would bind to biotic ligands in the order Ca2+ > Mg2+ > Na+, a sequence that is corroborated by the present study.
The average log KNaBL values from metal-specific studies were found to be 2.2 ± 0.3 for Zn2+, 2.1 for Cd2+ and 2.3 ± 0.1 for Cu2+, with no data available for the four other divalent metals (uncertainties apply where studies reported five or more values). These averages align with the overall value of 2.2 ± 0.1 (Table 2), with no significant differences observed (p > 0.05) across metals, indicating metal binding occurs at the same or similar receptor sites on the cell surface.
Sodium transporters on cell surfaces are widely recognised as essential for Na+ uptake by freshwater organisms [52], much like Ca and Mg transporters, providing a mechanistic basis for using Na+ binding constants. Moreover, Na transporters, including electrogenic Na+ channels and Na+/H+ exchangers found on the apical cell surface, are also capable of transporting non-essential metals of a similar size to Na, such as lithium (Li+). Of the four ameliorative cations, Na+ has the weakest binding affinity (or lowest average log KML value), being 3–4 orders of magnitude lower than H+, 16 times lower than Ca, and four times lower than Mg2+, at cell surface receptors. Some researchers have suggested that the “Na-effect” may be due to changes in activity coefficients and/or a charge screening effect on cell surface receptors, rather than from Na binding at Na receptors (transporters). This hypothesis is still largely untested, but it could offer a valuable complementary perspective.

3.12. Protons (H+)

Proton binding constants (log KHBL) can be derived from experimental studies where freshwater organisms exposed to metals at varying pH levels show reduced metal toxicity as pH decreases. Table 2 provides the proton binding constants for ten freshwater organisms from three taxa, with an average log KHBL value of 5.8 ± 0.1 (95% confidence limit). Additional information is provided in Table S2 and the Supplementary Materials. The log KHBL values were evenly distributed across the three major taxa: four crustacean species, three fish species and three microalgal species. The average log KHBL values were not significantly (p ≤ 0.05) different among crustaceans (5.7 ± 0.1), fish (5.7 ± 0.2) and microalgae (5.9 ± 0.2) (Table 2). These values are within the range previously reported for invertebrates (5.4–6.1) and fish (5.0–6.7) [10].
The average log KHBL values from metal-specific studies were found to be 5.9 ± 0.1 for Zn2+, 5.7 for Ni2+, 5.8 ± 0.1 for Cd2+, 5.6 ± 0.2 for Cu2+, 5.9 for UO22+ and 5.8 for Pb2+, with no data available for Co2+ (uncertainties apply where studies reported five or more values). These averages align with the overall value of 5.8 ± 0.1 (Table 2), with no significant differences observed (p > 0.05) across metals, indicating metal binding occurs at the same or similar receptor sites on the cell surface. The magnitude of the proton (H+) binding constant indicates that the cell surface receptor is a carboxylate group [10].

3.13. Overview

The average log KMBL values for all freshwater taxa indicate the following ranking of metal binding affinity to biotic ligands: Cu2+ > UO22+ > Cd2+ ≈ Pb2+ > Zn2+ ≈ Co2+ > Ni2+ (see Table 1). Copper had the highest binding affinity (log KCuBL = 7.9), while Ni had the lowest (log KNiBL = 4.8), with a ~1000-fold difference. This relative ranking largely reflects the four major taxa, which are highly consistent: fish (Cu2+ > UO22+ ≈ Cd2+ ≈ Pb2+ > Zn2+ ≈ Co2+ > Ni2+), crustaceans (Cu2+ ≈ UO22+ > Cd2+ ≈ Pb2+ > Zn2+ ≈ Co2+ > Ni2+), molluscs (Cu2+ ≈ UO22+ > Cd2+ ≈ Pb2+ > Zn2+ ≈ Co2+ > Ni2+), and microalgae (Cu2+ ≈ UO22+ > Cd2+ ≈ Pb2+ > Zn2+ ≈ Co2+ > Ni2+). For the ameliorative cations, the ranking of the average log KMBL values across all freshwater taxa was: H+ >> Ca2+ > Mg2+ > Na+. Protons (H+) had the highest binding affinity (log KHBL = 5.8), while Na+ had the lowest (log KNaBL = 2.2), with a ~4000-fold difference.
An average binding constant (log KML) for a specific metal can be applied across all freshwater taxa, though taxa- or species-specific log KML values may be applied if available. The uncertainty (error), or variability, in the average log KML values is relatively small (typically ± 0.1 log units; Table 1 and Table 2) and unlikely to significantly affect the overall accuracy or reliability of the BLM for predicting metal toxicity. For Cu, data from 42 species across nine phyla, resulted in an average log KML value with a relatively small uncertainty (8.0 ± 0.1), while for Co, data from only 15 species across six phyla led to a higher uncertainty for the average log KML value (5.8 ± 0.2). For new freshwater organisms (not yet part of the database) that are highly sensitive to certain metals, the use of average metal binding constants may be constrained due to unknown (fC)M50% factors. The accuracy or reliability of the BLM is most likely to be influenced by metal speciation (calculated via speciation models), specifically the percentage of free metal ions (M2+) for Zn, Ni, Cd and Co, or the percentages of M2+ and MOH+ for Cu, UO2 and Pb. Metal speciation is largely governed by key water chemistry parameters such as alkalinity, dissolved organic carbon (DOC) concentration and pH.
A key challenge in applying the BLM to less-studied metals like Co and UO2 is the limited toxicity and uptake data (compared to the other five metals studied) across freshwater species, which increases uncertainty (or error) in the binding constants (log KML values). This can lead to greater variability between measured and predicted toxicity, resulting in fewer data points within a factor of two of the 1:1 agreement line (as shown for Co and UO2 in Section 4). Both Co and UO2 tend to interact with the same cell surface receptors or transporters as other divalent metals like Cd, Pb, Ni and Zn, but with different binding affinities, affecting toxicity. Section 4 evaluates the reliability of the log KMBL values in predicting metal toxicity, and Section 5 explores the link between metal binding constants/affinity and toxicity.
The ameliorative cations (H+, Ca2+, Mg2+ and Na+) compete with trace metals (like Zn2+, Ni2+, Cd2+ and Co2+) for binding sites on cell surface receptors (biotic ligands) in aquatic organisms. The ability to influence (or modulate) metal binding is determined by their binding affinity at these biotic ligands (i.e., H+ >> Ca2+ > Mg2+ > Na+; Table 2). The binding constants (log KMBL values) for each ameliorative cation were shown to be independent of the trace metal examined. However, a key distinction is that H+, Ca2+, Mg2+ and Na+ compete with both M2+ and MOH+ for Cu2+, UO22+ and Pb2+, while they only compete with M2+ for Zn2+, Ni2+, Cd2+ and Co2+. In practical terms, as water hardness increases (i.e., the concentrations of Ca and/or Mg increase), Ca2+, and to a lesser extent Mg2+, will more effectively compete with trace metals like Cd2+, Zn2+ and Pb2+ at cell surface receptors. Likewise, as conductivity increases, Na+ will compete more effectively with trace metals, including Cu2+. When pH decreases, the higher proton (H+) concentration (a ten-fold increase for every unit drop in pH) will increasingly compete with trace metals at cell surface receptors.
The seven divalent metals investigated bind to an oxygen-containing donor ligand, such as carboxylate (-COOH). However, their binding affinities vary considerably, spanning three orders of magnitude, with Ni exhibiting the lowest affinity and Cu the highest (Table 1). Each metal’s binding affinity for cell surface receptors is influenced by its affinity for carboxylate. This is illustrated in Figure S2, which shows a strong positive linear regression (r2 = 0.97; p ≤ 0.001) between the binding affinity of six of the divalent metals (log KMBL) with cell surface receptors and acetate, a common monocarboxylic acid (carboxylate ligand). Uranium (UO22+) is excluded from this analysis because its two oxygen atoms cause steric hinderance, which reduces its binding to cell surface receptors and results in a lower log KML value than expected based on its affinity for acetate. The linear free energy relationship (LFER) shown in Figure S2 indicates that the binding constants of metals to biotic ligands are proportional to their binding with simple carboxylate ligands. Carbonaro and Di Toro [53] further explored LFERs between metal and proton (H+) binding constants of carboxylate, phenolic and hydroxo-ligands, deriving Irving-Rossotti slopes, which aligned well with metal binding constants for natural DOM (humic and fulvic acids). Given the similar metal binding behaviour of humic and fulvic acids to that of cell membrane surfaces, LFERs are a useful tool for determining metal binding constants to carboxylate-based receptors on cell surfaces.
To determine if the log KMBL values derived for freshwater organisms could extend to other organisms, a selection of metal toxicity studies was evaluated for marine and terrestrial organisms (see Supplementary Materials). The log KMBL values were consistent with those from freshwater organisms, suggesting a wider relevance of the data, although a more detailed comparison is necessary.

4. Comparison of Measured and Predicted Metal Toxicity

Figure 3 illustrates the comparison between the measured and predicted toxicity values (EC/LC50) for Zn2+, Ni2+, Cd2+ and Co2+, while Figure 4 shows the comparison for Cu2+, UO22+ and Pb2+. Predicted EC/LC50 values were derived using Equation (12), incorporating the average log KMBL values for Ca2+, Mg2+, Na+ and H+ (Table 2) and adjustments for M2+ for all seven metals, plus 0.5 × MOH+ for Cu2+, UO22+ and Pb2+, in each water quality scenario. The comparisons assess the accuracy of the derived log KMBL values, as well as those for the ameliorative cations (Ca2+, Mg2+, Na+ and H+) and the derived (fC)M50% values. Although previous studies have compared Cd2+, Cu2+ and Zn2+, this study uses 2–3 times more toxicity data, whilst maintaining or improving predictive accuracy. The Ni2+ and Pb2+ datasets are comparable to those from Santore et al. [43] and DeForest et al. [48], respectively. Additionally, this study provides the first comparisons for Co2+ and UO22+. The accuracy of the predictions is reflected by how closely the measured and predicted toxicity values match. This study found that 75 to 88% of toxicity values fell within a factor of two of the ideal 1:1 agreement line (Table 3 and Figure 3 and Figure 4). The best agreement was observed for Cd2+ (88%), while the worst was for Co2+ and UO22+ (75%), which had the smallest datasets and the largest uncertainties in their log KMBL values (Table 1, Table 2 and Table 3). Furthermore, 89 to 98% of toxicity values were within a factor of three of the ideal 1:1 agreement line (Table 3 and Figure 3 and Figure 4). Overall, the strong alignment between measured and predicted toxicity for all seven metals supports the reliability of the derived log KMBL and (fC)M50% values. An overview for each metal is provided in Section 4.1, Section 4.2, Section 4.3, Section 4.4, Section 4.5, Section 4.6 and Section 4.7.
The comparisons presented in Figure 4 for Cu, UO2, and Pb are based on predicted toxicity that considers both M2+ and (0.5×) MOH+. Among the seven metals, MOH+ occurs in test solutions only for these three metals. The significance of including MOH+ in toxicity assessments is most evident with uranium (UO2), as UO2OH+ has a higher concentration than UO22+ in nearly all studies (since it hydrolyses at about pH 5). A reanalysis of the predicted toxicity of uranium using only UO22+ (rather than UO22+ plus 0.5 × UO2OH+) in Figure 5 shows that just one out of the 52 data points (2% of data) falls within a factor of two of the ideal 1:1 line, with eight data points within a factor of three (15%), but the majority (85%) fall well outside the factor of three line. This suggests that uranium would be much more toxic than observed.
The differences for Cu and Pb are less pronounced because they hydrolyse around pH 7, resulting in substantial pH regions (pH 5.5–7.5) where the M2+ ion dominates. However, toxicity is still predicted to increase significantly for both metals, leading to notable discrepancies between measured and predicted toxicity. For Pb, the average log KPbBL value increased significantly (p ≤ 0.05) from 6.8 ± 0.1, when both Pb2+ and PbOH+ are considered, to 7.2 ± 0.2, when only Pb2+ is considered. Consequently, the average predicted Pb toxicity would increase (EC/LC50 values would decrease), by a factor of 2.5, indicating that fewer data points would lie above the 1:1 line, with many falling below the lower factor of three line. Similarly for Cu, the average log KCuBL value increased significantly (p ≤ 0.05) from 7.9 ± 0.1, when both Cu2+ and CuOH+ are considered, to 8.2 ± 0.1, when only Cu2+ is considered. Consequently, the average predicted Cu toxicity would increase (EC/LC50 values would decrease), by a factor of two, again indicating that a significant portion of the data would fall below the lower factor of three line when comparing measured and predicted data. Therefore, it is evident that in environmental pH ranges where MOH+ forms, both M2+ and MOH+ need to be accounted for as contributors to toxicity. Consulting with BLM practitioners is necessary to identify the best approach for integrating this into current software. The simplest solution would be to use the sum of the concentrations of M2+ and 0.5 × MOH+ from speciation modelling calculations for Cu and Pb (at pH ≥ 7.0) or UO2 (at pH ≥ 5.0) to improve, or optimise, the BLM.
The conceptual model demonstrates how a cascade of effects result in the concentration at which metals become toxic. Understanding these effects is vital for quantifying the varying toxicity metals impose on aquatic organisms. In many aquatic organisms, toxicity occurs when 50% of biotic ligands are occupied ((fC)M50%). However, in some cases, toxicity can be higher (i.e., lower TR50) because of lower ligand occupancy at the TR50, while other organisms experience reduced toxicity (due to protective mechanisms), despite a constant metal binding constant. These organisms might have additional low-affinity, high-capacity receptors that effectively reduce the effect of the metal. Further research is needed on the underlying physiological mechanisms governing (fC)M50% (i.e., why is organism sensitivity enhanced?), as well as the relationship between toxicity and (fC)M50%.
The (fC)M50% factor may be applied to multiple metals for a single species or to the same metal across related species. Further research is needed to quantify (fC)M50% for a wider range of freshwater species. Collaboration between BLM practitioners is required to determine the most effective way to integrate an adjustable (fC)M50% factor into existing software.

4.1. Zinc

Figure 3a compares the measured and predicted EC/LC50 values for 142 water quality scenarios across seven organisms from 13 studies. The predicted EC/LC50 values incorporate the average log KMBL values for Zn2+ (Table 1), Ca2+, Mg2+, Na+ and H+ (Table 2), and adjustments for Zn2+ activity in each scenario. Figure 3a shows that the predicted EC/LC50 values align closely with the measured values over nearly three orders of magnitude, with 81% of data points within a factor of two, and 94% within a factor of three, of the ideal 1:1 line. The data points appear to be uniformly distributed around the 1:1 line.
The sensitivity factor ((fC)Zn50%) was applied to eight organisms (two fish species, four crustacean species and a microalgae and insect), ranging from 0.0521 to 0.45 (Table S4)—with the default maximum of 0.5 being applied to all other species. According to Equation (12), as (fC)Zn50% decreases, the Zn EC/LC50 should also decrease. Figure 6 illustrates the relationship between (fC)Zn50% and the corresponding EC/LC50 values. Despite the small dataset (n = 8) leading to an insignificant (p > 0.05) linear relationship, the data still suggest that a decrease in (fC)Zn50% is associated with a decrease in EC/LC50.

4.2. Nickel

Figure 3b compares the measured and predicted EC/LC50 values for 256 water quality scenarios across ten organisms from 14 studies. The predicted EC/LC50 values incorporate the average log KMBL values for Ni2+ (Table 1), Ca2+, Mg2+, Na+ and H+ (Table 2), and adjustments for Ni2+ activity in each scenario. Figure 3b shows that the predicted EC/LC50 values align closely with the measured values over nearly five orders of magnitude, with 83% of data points within a factor of two, and 97% within a factor of three, of the ideal 1:1 line. The data points appear to be uniformly distributed around the 1:1 line.
The sensitivity factor ((fC)Ni50%) was applied to 17 organisms (two fish species, eight crustacean species, two mollusc species, two microalgal species, two macrophyte species and a rotifer), ranging from 0.0033 to 0.499 (Table S4)—with the default maximum of 0.5 being applied to all other species. According to Equation (12), as (fC)Ni50% decreases, the Ni EC/LC50 should also decrease. Figure 6 illustrates the positive linear relationship (r2 = 0.76; p < 0.001) between (fC)Ni50% and the corresponding EC/LC50 values. Therefore, as (fC)Ni50% decreases, freshwater organisms become more sensitive to Ni.

4.3. Cadmium

Figure 3c compares the measured and predicted EC/LC50 values for 83 water quality scenarios across six organisms from seven studies. The predicted EC/LC50 values incorporate the average log KMBL values for Cd2+ (Table 1), Ca2+, Mg2+, Na+ and H+ (Table 2), and adjustments for Cd2+ activity in each scenario. Figure 3c shows that the predicted EC/LC50 values align closely with the measured values over two orders of magnitude, with 88% of data points within a factor of two, and 98% within a factor of three, of the ideal 1:1 line. The data points appear to be uniformly distributed around the 1:1 line.
The sensitivity factor ((fC)Cd50%) was applied to two organisms (O. mykiss and Hyalella azteca), ranging from 0.00676 to 0.216 (Table S4)—with the default maximum of 0.5 being applied to all other species. Two of the three (fC)Cd50% values that were less than 0.5 were based on studies investigating the Cd toxicity and/or uptake in O. mykiss [77,109]. Liao et al. [77] reported a (fC)Cd50% value of 0.0128, with a range of 0 to 0.0634. A re-analysis of their data produced a (fC)Cd50% of 0.00676. Niyogi et al. [77] determined a log KCdBL of 7.5 from uptake measurements, but a log KCdBL of 8.0 from toxicity measurements (see Section 6). To harmonise the binding constant between these measurements, a (fC)Cd50% value of 0.216 was attributed to the toxicity data. Additionally, a (fC)Cd50% value of 0.184 was derived using data from Schroeder [72] for the toxicity of Cd2+ to H. azteca. The (fC)Cd50% value derived from Liao et al. [109] is substantially lower than corresponding values derived from Niyogi et al. [77] or Schroeder [72].

4.4. Cobalt

Figure 3d compares the measured and predicted EC/LC50 values for 12 water quality scenarios across nine organisms from four studies. The predicted EC/LC50 incorporate the average log KMBL values for Co2+ (Table 1), Ca2+, Mg2+, Na+ and H+ (Table 2), and adjustments for Co2+ activity in each scenario. This calculation accounts for both the competitive/ameliorative effects of Ca2+, Mg2+, Na+ and H+ and the effect of speciation (complexation). Figure 3d shows that the predicted EC/LC50 values align closely with the measured values over two orders of magnitude, with 75% of data points within a factor of two, and 92% within a factor of three, of the ideal 1:1 line. There appears to be a slight bias of data points below the 1:1 line, but the dataset is limited.
The sensitivity factor ((fC)Co50%) was applied to only one organism (the crustacean, C. dubia), with a value of 0.0521 (Table S4)—with the default maximum of 0.5 being applied to all other species.

4.5. Copper

Figure 4a compares the measured and predicted EC/LC50 values for 383 water quality scenarios across 12 organisms from 15 studies. The predicted EC/LC50 values incorporated the average log KMBL values for Cu2+ (Table 1), Ca2+, Mg2+, Na+ and H+ (Table 2), and adjustments for Cu2+ and (0.5×) CuOH+ activities in each scenario. Figure 4a shows that the predicted EC/LC50 values align closely with the measured values over three orders of magnitude, with 79% of data points within a factor of two, and 91% within a factor of three, of the ideal 1:1 line. The data points appear to be uniformly distributed around the 1:1 line.
The sensitivity factor ((fC)Cu50%) was applied to 12 organisms (nine crustacean species, two microalgal species and an insect), ranging from 0.002 to 0.499 (Table S4)—with the default maximum of 0.5 being applied to all other species. According to Equation (12), as (fC)Cu50% decreases, the Cu EC/LC50 should also decrease. Figure 6 illustrates the positive linear relationship (r2 = 0.73; p ≤ 0.01) between (fC)Cu50% and the corresponding EC/LC50 values. Therefore, as (fC)Cu50% decreases, freshwater organisms become more sensitive to Cu. Interestingly, the slope obtained for the linear regression for Cu (4.2 ± 0.8) is similar to that found for nickel (3.0 ± 0.4), within uncertainty. Both metals also have an extremely low (fC)M50% value in their dataset, which seems to depart from the rest of the linear relationship obtained for each metal (see Figure 6).

4.6. Uranium

Figure 4b compares the measured and predicted EC/LC50 values for 52 water quality scenarios across ten organisms from six studies. The predicted EC/LC50 values incorporate the average log KMBL values for UO22+ (Table 1), Ca2+, Mg2+, Na+ and H+ (Table 2), and adjustments for UO22+ and (0.5×) UO2OH+ activities in each scenario. Figure 4b shows that the predicted EC/LC50 values align closely with the measured values over two orders of magnitude, with 75% of data points within a factor of two, and 88% within a factor of three, of the ideal 1:1 line. The data points appear to be uniformly distributed around the 1:1 line.
The sensitivity factor ((fC)UO250%) was applied to only two organisms (the crustacean, C. dubia, and the macrophyte, Ceratophyllum demersum), with values of 0.0521 and 0.47 (Table S4), respectively—with the default maximum of 0.5 being applied to all other species. The (fC)UO250% value of 0.052 was used for the toxicity data of Goulet et al. [95] where C. dubia was studied (the same value as used for other metals, except Cd).

4.7. Lead

Figure 4c compares the measured and predicted EC/LC50 values for 110 water quality scenarios across five organisms from eight studies. The predicted EC/LC50 values incorporate the average log KMBL values for Pb2+ (Table 1), Ca2+, Mg2+, Na+ and H+ (Table 2), and adjustments for the Pb2+ and (0.5×) PbOH+ activities in each scenario. Figure 4c shows that the predicted EC/LC50 values align closely with the measured values over three orders of magnitude, with 80% of data points within a factor of two, and 95% within a factor of three, of the ideal 1:1 line. The data points are evenly distributed around the 1: 1 line.
The sensitivity factor ((fC)Pb50%) was applied to only three organisms (the crustaceans, C. dubia, D. magna and Daphnia similis), with values between 0.0521 and 0.24 (Table S4)—with the default maximum of 0.5 being applied to all other species.

5. Relationships between Metal Binding Constants and Toxicity

Since both acute and chronic toxicity endpoints are commonly used in aquatic ecotoxicology, it was important to incorporate all available metal toxicity data to derive metal binding constants (log KML values) for the BLM. The toxicity endpoints, which ranged from acute lethal effects (e.g., 2-day survival for crustaceans) to chronic sublethal effects (e.g., 28-day growth for gastropod molluscs), exhibited a bimodal distribution across taxa. Metal binding constants were largely derived from acute toxicity data for fish (84%), crustaceans (85%) and rotifers (75%), while chronic toxicity data dominated for molluscs (90%), microalgae (100%) and macrophytes (100%). No significant (p > 0.05) differences in Zn, Cu or Pb binding constants were observed between acute and chronic exposures in the cladoceran C. dubia (the only species where sufficient data were available for a direct comparison), with similar results for Zn across eight cladoceran species and for Cu in gastropods. Although more data are needed for other species/taxa comparisons, the current findings support combining acute and chronic toxicity (LC/EC50) data to derive metal binding constants. Moreover, metal binding to cell surface receptors should be consistent irrespective of the exposure duration or type (lethal or sublethal effects).
The log KMBL values for seven freshwater organisms (O. mykiss, P. promelas, C. dubia, D. magna, H. azteca, R. subcapitata and Lemna minor) were derived for all seven metals (Cd2+, Co2+, Cu2+, Ni2+, Pb2+, UO22+ and Zn2+). For these organisms, the measured toxicity of each metal is a function of the derived log KMBL value, the (fC)M50% value, and the concentrations of ameliorative cations. An assessment of metal toxicity can be made through the linear regression of the log KMBL values among the organisms. Table S5 presents the coefficients of determination (r2 values) and slopes of the linear regressions fitted to the log KMBL values across the seven organisms. All regressions showed strong positive linear relationships (p < 0.01), with r2 values between 0.82 and 0.99 (average of 0.93), and slopes ranging from 0.79 to 1.23 (average of 1.00 ± 0.12). The average slope met the expected slope of unity, within uncertainty.
Given that (fC)M50% values and the concentrations of ameliorative cations can strongly affect log KMBL values, it is important to apply a consistent (fC)M50% value when comparing toxicity across various organisms for all seven metals. Linear regressions of metal toxicity (expressed as 1/[M], where [M] is the toxicity concentration in µmol/L) and metal binding constants (log KMBL values) were fitted for C. dubia (with a (fC)M50% value of 0.0521), D. magna and O. mykiss (both with a (fC)M50% value of 0.5). Figure 7 shows the positive linear regressions (p ≤ 0.01) for these species, with r2 values ranging from 0.91 for C. dubia to 0.97 for O. mykiss, indicating that log KMBL values explain between 91% and 97% of the variability in metal toxicity. As predicted by Equation (19), the slope of the relationship between log KMBL and log (1/[M]) should be one. The regression slopes for C. dubia, D. magna and O. mykiss ranged from 0.88 to 1.10, with an average slope of 1.02, which is in close agreement with the expected slope of 1.00. Overall, the results indicate that higher log KMBL values correspond to increased metal toxicity in these freshwater organisms.
The linear regression for C. dubia (Figure 7a) omits a Cd2+ value because its (fC)M50% value was 0.5, while the other six metals had a value of 0.052. Similarly, O. mykiss lacks a Cd2+ value in the regression (Figure 7b) due to its (fC)M50% value of 0.216, compared to 0.5 for the other six metals. To illustrate the effect of (fC)M50% on metal toxicity, Cd2+ values were added as an overlay on Figure 7b using both the uncorrected (0.216; open red circle) and corrected (i.e., × (1 − 0.216)/0.216); open purple circle) (fC)M50% value. The correction brings Cd2+ in line with the other six metals.

6. Organism Sensitivity and Its Relationship to the Fraction of Cells Bound at the TR50 ((fC)M50%)

The toxic response of aquatic organisms to metals is driven by a sequence, or cascade, of effects, as described in Equation (2). This cascade is a key feature of the extended FIAM, but is not explicit in the BLM. The extended FIAM highlights two distinct and sequential processes in the toxic response to metals. Additionally, a third mechanism is revealed in the cascade when the TR50 is reached with less than 50% of cell surface receptor sites occupied (Equation (8)). Figure 2 illustrates this cascade, showing the key terms (or parameters) at each stage.
Altszyler et al. [110] demonstrated that a cascade of interactions can result in an ultrasensitive response. In the context of the metal toxicity in aquatic organisms, the (fC)M50% parameter can serve as an indicator of an organism’s sensitivity to metals. Altszyler et al. [110] showed that a cascade of effects can amplify toxicity, making organisms more sensitive. Changes in the sensitivity could lead to significant variation in an organism’s response to metal toxicants. Sensitive organisms may transport metals into cells more rapidly than less sensitive organisms, even though the metal binding affinity is equivalent. This behaviour suggests that organisms can employ defensive mechanisms, such as mucous coatings to reduce metal binding at cell surface receptors, providing protection from toxicity. Such effects can be further elucidated from studies that have examined both metal uptake by, and toxicity to, freshwater organisms.
Niyogi et al. [77] investigated the uptake and toxicity of Cd in rainbow trout (O. mykiss). From their uptake data, they calculated a log KCdBL value of 7.5 (Table S1), and also determined a 96 h LC50 value of 10.3 nmol/L (corrected for the effects of Ca and DOC) from toxicity tests. Using Equation (9), this LC50 value results in a log KCd-tox of 8.0, which the authors referred to as the “true estimate” of Cd binding affinity at the biotic ligand. However, metal uptake measurements reflect the true binding affinity of metals for biotic ligands. The concentration at which LC50 or EC50 values are observed is related to the binding constant through the sensitivity factor, (fC)M50%, as defined in Equation (8), and re-expressed below in terms of KMBL and KM-tox:
M = f C M 50 % 1 f C M 50 % 1 K MBL = 1 K M - tox
The log KCdML (termed log KCd-gill by the authors]) and log KCd-tox values from Niyogi et al. [52] can be equated using Equation (19) and a (fC)M50% of 0.249 was calculated (i.e., almost a quarter of the cell surface binding sites are occupied at the LC50).
Niyogi et al. [111] previously investigated Cd uptake and toxicity in O. mykiss, deriving a log KCdML value of 7.3 from uptake data (listed in Table S1), which aligned with their later findings. The toxicity tests produced an LC50 of 169 nmol/L for total Cd, which when corrected for Cd speciation and the concentrations of Ca and DOM resulted in an adjusted value of 27.1 nmol/L. This led to a log KCd-tox of 7.6, larger than the log KCdML value. From these results, a (fC)M50% of 0.372 was calculated. The log KCdML, log KCd-tox and (fC)Cd50% values are in good agreement across both studies. The study also investigated Cd uptake and toxicity in yellow perch (Perca flavescens), obtaining a log KCdML value of 7.2 from the uptake data (listed in Table S1), which is similar to the value for O. mykiss (7.5). Based on Cd toxicity data, P. flavescens had an LC50 of 72,420 nmol/L, more than 400 times higher than for O. mykiss. After accounting for Cd speciation and concentrations of Ca and DOM, the corrected LC50 value for P. flavescens was 6140 nmol/L, resulting in a log KCd-tox value of 5.2, two orders of magnitude lower than for O. mykiss.
Although both fish species share a similar Cd binding capacity, P. flavescens appears to have developed a defence mechanism that protects against Cd exposure. This is reflected in the substantially higher 3 h and 24 h LA50 values (17.80 and 31.67 nmol per g wet gill tissue, respectively) compared to O. mykiss (0.34 and 0.55 nmol per g wet gill tissue, respectively). For O. mykiss, both LA50 values were well within the concentration of saturable sites at the gill (0.85 nmol per g wet gill tissue) whereas for P. flavescens, both LA50 values were considerably larger than the concentration of saturable sites (0.67 nmol per g gill wet tissue). Consequently, in P. flavescens, in addition to filling up to 50% of the saturable binding sites (at the LA50), a substantial amount of other low-affinity, high-capacity sites were filled [111], which enables the species to protect itself from Cd exposure, even though the binding affinity for the saturable sites is very similar to that of O. mykiss. It appears that P. flavescens has another layer in the cascade of interactions that enables it to reduce the effect of exposure to Cd. Conversely, for O. mykiss, the two LA50 values occur around the concentration at which at least half of the receptor (saturable) sites are occupied, with only a small fraction of other sites also being occupied. Similar behaviour for Cu was previously reported by Taylor et al. [112].
Leonard et al. [113] investigated Ni uptake and toxicity in two species of fish, O. mykiss and the round goby (Neogobius melanostomus). They found that the log KNiBL values derived from uptake data were up to an order of magnitude (one log unit) higher than those derived from toxicity data for both species. However, these results need to be reassessed. For N. melanostomus, the authors found a Kd value of 17.8 μmol/L for Ni uptake, corresponding to a log KNiBL value of 4.8, which is in good agreement with the values listed in Table 1. Using toxicity data, the authors measured an LC50 of 104.1 μmol/L, which reduced to 80.2 μmol/L after adjusting for Ni speciation. After further adjustments for pH and the concentrations of Ca2+, Mg2+ and Na+, applying average binding constants from this study, the LC50 reduced to 21.1 μmol/L. This corresponds to a log KNiBL value of 4.7, which aligns with their uptake data and the results of this study (Table 1).
For O. mykiss, Leonard et al. [113] found a Kd value of 46.5 μmol/L from whole fish uptake data, which was unusual as the value for N. melanostomus was derived only from the gills. The Kd corresponds to a log KNiBL of 4.3, lower than measurements from other studies on O. mykiss (4.5–4.9; Table S1) and other freshwater taxa (4.6–5.2; Table 1). Based on toxicity data for O. mykiss, the LC50 concentration of 228.1 μmol/L was reduced to 175.6 μmol/L after adjustment for Ni speciation (Ni2+), and further to 46.9 μmol/L after accounting for the ameliorative cations. This final concentration corresponds to a log KNiBL value of 4.3, matching the uptake data. This analysis underscores the importance of accounting for the free metal ion in toxicity tests (but not uptake tests), and the ameliorative effects of cations (Ca2+, Mg2+, Na+, H+) when comparing log KMBL values derived from metal uptake and toxicity data.

7. The Binding Mechanism of MOH+

While free metal ions are often considered the best predictors of metal uptake and toxicity in freshwater organisms, several studies have shown that metal hydroxide complexes (e.g., UO2OH+, CuOH+ and Sc(OH)2+/Sc(OH)2+) improve the prediction of free metal ion behaviour [9,14,15,16]. Hydroxide complexes become important in freshwater systems after metal hydrolysis, typically starting around pH 5.0 for UO22+ and pH 7.0 for Cu2+ and Pb2+. Although PbOH+ has yet to be empirically confirmed as complementary to Pb2+, it is expected to behave similarly to Cu, based on its first hydrolysis constant [25].
Equations (15) and (16) were formulated based on empirical evidence that KLMBL/KMBL = 0.5 [9,14,15], indicating that the relative effects of M and ML at cell surface receptors correlate to their charges. The binding of M and ML at cell surface receptors (biotic ligands) can be described by the following equations:
M + BL   MBL
ML + BL   LMBL
where BL (biotic ligand) is the binding site, M is the free metal ion and ML is the complexed metal ion with ligand L (with OH as L in this study). Chemical binding reactions are often divided into intrinsic (chemical affinity) and coulombic (electrostatic interactions) components [19]. For reactions (20) and (21), the intrinsic binding constants (Kint) are likely similar, given that both reactions involve the same ligand (cell surface receptor) and metal. However, the coulombic binding constants (Kcoul), or electrostatic interactions, will vary due to charge differences between M and ML.
Brown and Sylva [114] developed a theoretical model to describe metal–ligand interactions, including those at cell surface receptors. Their model predicts that when ML binds to a cell surface receptor (biotic ligand), it will have a similar intrinsic binding constant to that of M binding to the same receptor, but with different electrostatic interactions. This difference is determined by the ratio of two factors: (i) the number of binding groups to the central metal ion (M interacts solely with the receptor (one group), while the metal ion in ML interacts with both L (such as OH) and the receptor (two groups), creating a 2:1 ratio), and (ii) a proportionality constant (k), which applies specifically to ML binding. Thus, the ratio between these two factors is 2/k, and for this ratio to reflect the observed two-fold binding difference between ML and M at cell surface receptors, the value of k must be one.

8. A Second Cell Surface Binding Site at a Moderately Alkaline pH?

According to ETMs, a second neutral site on the cell surface could bind metals in moderately alkaline (pH 7.0 to 8.5) environments [7,8]. This interaction is described by Equation (18). At these pH values, metals are more likely to bind with an amine group (-NH2) on the cell surface, while at a lower pH (5.5–7.0), metal binding with a carboxyl group (-COOH) is most likely. Analysis of the data from these studies shows both a second (higher) proton binding constant, consistent with binding at a higher pH, and higher metal binding constants to the second cell binding site—consistent with the correlation between binding stability and proton dissociation [25].
Deleebeeck et al. [69] observed a pH breakpoint in their study on Ni toxicity to O. mykiss (Figure S3). At or below a pH of 6.8, the log KHBL value was 5.7, consistent with other data. Above this pH, the log KHBL value increased to 6.8. This is the only study on this species that extensively examines the role of proton (H+) concentration in Ni toxicity, and one of three studies that explore the effect of pH (all three studies are reviewed here). Deleebeeck et al. [68] later observed similar behaviour in Ni toxicity to R. subcapitata (Figure S4), with a pH breakpoint at 7.2. Below this pH, a log KHBL value of 5.7 was derived, and above it, a value of 6.7, consistent with the value from their previous study. Similarly, Schroeder [72] reported comparable behaviour for H. azteca (Figure S5), although the pH breakpoint occurred at a higher pH of 8.0. Below this pH, a log KHBL value of 5.5 was derived, and above it, a value of 6.9, consistent with earlier studies by Deleebeeck et al. [68,69]. The trend in pH breakpoints corresponds with the (fC)Ni50% values across the studies: as the pH breakpoint increases, the (fC)Ni50% value decreases. The average log KHBL value from the three studies, possibly related to an amine site, is 6.8.
Heijerick et al. [61] investigated the toxicity of Zn to R. subcapitata and observed a similar pH-dependent behaviour as that found by Deleebeeck et al. [68] for Ni. Below and at a pH of 6.8, a log KHBL value of 5.8 was derived from their data, while above this pH, a larger log KHBL value of 8.3 was obtained (Figure S6). Wilde et al. [115] found a less distinct pH breakpoint (Figure S7), likely occurring around the same pH as Heijerick et al. [61]. Analysis of their data above pH 6.8 yielded a log KHBL value of 8.3, matching the earlier study of Heijerick et al. [61]. Van Regenmortel et al. [59], using a different model to assess Zn toxicity to R. subcapitata, could not identify a clear pH breakpoint but derived a high-pH log KHBL value of 7.9, reasonably consistent with the other two values. The average log KHBL value from these Zn studies is 8.2, significantly higher than the value obtained for Ni.
Few studies have identified a pH-dependent breakpoint in metal toxicity to freshwater organisms, and only for Ni and Zn, the metals with two of the lowest log KMBL values. This could imply that amine group binding is possible only when carboxyl group binding is relatively weak. However, given the limited research showing a pH-related toxicity breakpoint, additional studies are necessary to confirm this phenomenon. Although there is good agreement in log KHBL values for each metal individually (Ni or Zn), there is less agreement between these two metals. Consistent log KHBL values have been found for carboxyl binding sites across all seven metals, but the lack of consistency of log KHBL values reported for amine binding sites raises doubts about the reliability of the data. This reinforces the need for further research on this effect with Ni and Zn.

9. Conclusions and Recommendations

This study presents the integration of the extended free ion activity model (FIAM), biotic ligand model (BLM) and electrostatic toxicity model (ETM) into a unified conceptual framework to describe metal interactions at cell surface receptor sites (biotic ligands). The integrated model shows that (i) the free metal ion (Mz+) is primarily responsible for toxicity in aquatic organisms, although other cationic species (e.g., MOH(z−1)+, MF(z−1)+) may also contribute, (ii) toxicity can be ameliorated by major cations, including Ca2+, Mg2+, Na+ and H+ and (iii) toxicity is strongly influenced by the proportion of receptor sites (biotic ligands) occupied at the TR50 ((fC)M50%), where toxicity increases as (fC)M50% decreases. The key factor determining toxicity is the binding affinity of metals for biotic ligands; the stronger the affinity (or higher the binding constant, log KMBL), the greater the toxicity to aquatic organisms.
A comprehensive literature evaluation on the toxicity and uptake of seven key metals (Ni2+, Zn2+, Co2+, Cd2+, Cu2+, UO22+ and Pb2+) and four ameliorative cations (H+, Na+, Ca2+ and Mg2+) has produced an extensive database of binding constants (log KMBL values) for freshwater organisms. This database is a key tool for predicting metal toxicity to aquatic organisms. The binding constants incorporate the competitive effects of the ameliorative cations, the fraction of cells occupied by a metal at the TR50 ((fC)M50%) and metal speciation. For metals like Cu2+, UO22+ and Pb2+ that hydrolyze in circumneutral freshwater (pH 5.5–8.5), MOH+ was included by adding half of its activity at the TR50 to that of Mz+. The average binding constants (log KMBL values) for Ca2+ (3.4), Mg2+ (2.8) and Na+ (2.2) are notably lower than those of the seven selected metals (4.8–7.9). The proton (H+) binding constant (5.8) is similar to, or higher than, those of Ni2+ (4.8), Zn2+ (6.0) and Co2+ (5.8), but lower than those for Cd2+ (7.0), Cu2+ (7.9), Pb2+ (6.8) or UO22+ (7.7). There is a 1000-fold difference in binding affinity between Ni2+ (lowest) and Cu2+ (highest) for biotic ligands. In general, the average binding constants (log KMBL values) derived from toxicity studies were consistent with those derived directly from uptake studies. The derived binding constants will improve current BLMs for Cd2+, Cu2+, Ni2+, Pb2+ and Zn2+ and aid in developing new BLMs for Co2+ and UO22+.
The derived metal binding constants accurately predicted metal toxicity for a wide range of freshwater organisms, where 75–88% of data were within a factor of two, and 88–98% of data were within a factor of three, of the ideal 1:1 agreement line. The best agreement was observed for Cd2+ (88%) while the worst was for Co2+ and UO22+ (75%), which had the smallest datasets and the largest uncertainties in their log KMBL values (Table 1, Table 2 and Table 3). Overall, the strong alignment between measured and predicted toxicity for all seven metals supports the reliability of the derived log KMBL and (fC)M50% values. Additional toxicity and uptake experiments are required for Co2+, UO22+ and Pb2+, under a range of water chemistry conditions, to provide further data on binding constants for comparing measured versus predicted toxicity. Experiments designed to evaluate the fraction of cells occupied by a metal at the TR50 ((fC)M50%) for key organisms across freshwater taxa are much needed, as some organisms show great variability in derived (fC)M50% values. It was generally concluded that higher log KMBL values correspond to increased metal toxicity in freshwater.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/w16202999/s1, Figure S1: Calculated concentration-response curves when (fC)M50% = 0.2 (less than 50% of receptor sites are occupied at the TR50), 0.5 (50% of receptor sites are occupied at the TR50) and 0.8 (more than 50% of receptor sites are occupied at the TR50). Since the true TRmax is not reached when (fC)M50% > 0.5, such scenarios are invalid. Figure S2: Linear regression between metal binding constants (log KMBL) with cell surface receptors and metal stability constants (log KML) with acetate (a simple carboxylate ligand). Figure S3: Linear regression between nickel toxicity (EC50) in Oncorhynchus mykiss (fish) and pH (as proton (H+) concentration) from Deleebeeck et al. (2007a). Figure S4: Linear regression between nickel toxicity (EC50) in Raphidocelis subcapitata (microalgae) and pH (as proton (H+) concentration) from Deleebeeck et al. (2009). Figure S5: Linear regression between nickel toxicity (EC50) in Hyallela azteca (crustacean) and pH (as proton (H+) concentration) from Schroeder (2008). Figure S6: Linear regression between zinc toxicity (EC50) in Raphidocelis subcapitata (microalgae) and pH (as proton (H+) concentration) from Heijerick et al. (2002b). Figure S7: Linear regression between zinc toxicity (EC50) in Chlorella sp. (microalgae) and pH (as proton (H+) concentration) from Wilde et al. (2006). Figure S8: Linear regression between sodium (log KNaBL) and copper (log KNaBL) binding constants for cladocerans (crustaceans) from De Schamphalaere et al. (2007). The red data point is for A. elongatus and was not used in the linear regression because it required a (fC)M50% factor greater than 0.5 (nevertheless, it is consistent with the linear relationship). Table S1: Binding constants (log KMBL) of zinc (Zn), nickel (Ni), cadmium (Cd), cobalt (Co), copper (Cu), uranium (UO2) and lead (Pb) at cell surface receptor sites (biotic ligands) in freshwater organisms (including the toxic endpoint for each study). Table S2: Binding constants (log KMBL) of the ameliorative cations (calcium (Ca), magnesium (Mg), sodium (Na) and protons (H+)) at cell surface receptors (biotic ligands) in freshwater organisms (including the toxic endpoint for each study). Table S3: Comparison of metal binding constants (log KMBL) at cell surface receptors (biotic ligands) from uptake and toxicity experiments. Table S4: Values derived for (fC)M50% (for values < 0.5). Table S5: Coefficients of determination (r2 values—above the diagonal) and slopes (below the diagonal) for linear regressions of log KMBL values among freshwater organisms (see Table S1 for full species names).

Author Contributions

Conceptualization, P.L.B. and S.J.M.; methodology, P.L.B. and S.J.M.; formal analysis, P.L.B. and S.J.M.; investigation, P.L.B. and S.J.M.; writing—original draft preparation, P.L.B. and S.J.M.; writing—review and editing, P.L.B. and S.J.M.; visualisation, P.L.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in this study are included in the article/Supplementary Materials. Further enquiries can be directed to the corresponding author.

Acknowledgments

The authors thank two anonymous reviewers for their constructive feedback on an earlier manuscript.

Conflicts of Interest

Paul. L. Brown was employed by the company Geochem Australia and Scott J. Markich by the company Aquatic Solutions International. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Martinez-Finley, E.J.; Chakraborty, S.; Fretham, S.J.B.; Aschner, M. Cellular transport and homeostasis of essential and nonessential metals. Metallomics 2012, 4, 593–605. [Google Scholar] [CrossRef] [PubMed]
  2. Weishaupt, A.K.; Lamann, K.; Tallarek, E.; Pezacki, A.T.; Matier, C.D.; Schwerdtle, T.; Aschner, M.; Chang, C.J.; Stürzenbaum, S.R.; Bornhorst, J. Dysfunction in atox-1 and ceruloplasmin alters labile Cu levels and consequently Cu homeostasis in C. elegans. Front. Mol. Biosci. 2024, 11, 1354627. [Google Scholar] [CrossRef] [PubMed]
  3. Campbell, P.G. Interactions between trace metals and aquatic organisms: A critique of the free-ion activity model in metal speciation and bioavailability. In Aquatic Systems; Tessier, A., Turner, D.R., Eds.; John Wiley and Sons: New York, NY, USA, 1995; pp. 45–102. [Google Scholar]
  4. Brown, P.L.; Markich, S.J. Evaluation of the free ion activity model of metal-organism interaction: Extension of the conceptual model. Aquat. Toxicol. 2000, 51, 177–194. [Google Scholar] [CrossRef] [PubMed]
  5. Di Toro, D.M.; Allen, H.E.; Bergman, H.L.; Meyer, J.S.; Paquin, P.R.; Santore, R.C. Biotic ligand model of the acute toxicity of metals. I. Technical basis. Environ. Toxicol. Chem. 2001, 20, 2383–2396. [Google Scholar] [CrossRef]
  6. Campbell, P.G.; Fortin, C. Biotic ligand model. In Encyclopedia of Aquatic Ecotoxicology; Férard, J.F., Blaise, C., Eds.; Springer: Dordrecht, The Netherlands, 2013; pp. 237–245. [Google Scholar]
  7. Wang, P.; Kinraide, T.B.; Zhou, D.; Kopittke, P.M.; Peijnenburg, W.J.G.M. Plasma membrane surface potential: Dual effects upon ion uptake and toxicity. Plant Physiol. 2011, 155, 808–820. [Google Scholar] [CrossRef]
  8. Zhou, D.M.; Wang, P. A novel approach for predicting the uptake and toxicity of metallic and metalloid ions. Plant Signal. Behav. 2011, 6, 461–465. [Google Scholar] [CrossRef]
  9. De Schamphelaere, K.A.C.; Janssen, C.R. A biotic ligand model predicting acute copper toxicity for Daphnia magna: The effects of calcium, magnesium, sodium, potassium, and pH. Environ. Sci. Technol. 2002, 36, 48–54. [Google Scholar] [CrossRef]
  10. Ardestani, M.M.; van Straalen, N.M.; van Gestel, C.A.M. The relationship between metal toxicity and biotic ligand binding affinities in aquatic and soil organisms: A review. Environ. Pollut. 2014, 195, 133–147. [Google Scholar] [CrossRef]
  11. Slaveykova, V.; Wilkinson, K.J. Predicting the bioavailability of metals and metal complexes: Critical review of the biotic ligand model. Environ. Chem. 2005, 2, 9–24. [Google Scholar] [CrossRef]
  12. Campbell, P.G.C.; Errécalde, O.; Fortin, C.; Hiriart-Baer, V.P.; Vigneault, B. Metal bioavailability to phytoplanktonApplicability of the biotic ligand model. Comp. Biochem. Physiol. Part C Toxicol. Pharmacol. 2002, 133, 189–206. [Google Scholar] [CrossRef]
  13. Rand, G.M.; Wells, P.G.; McCarty, L.S. Introduction to aquatic toxicology. In Fundamentals of Aquatic Toxicology, 2nd ed.; Rand, G.M., Ed.; Taylor and Francis: Washington, DC, USA, 1995; pp. 3–66. [Google Scholar]
  14. Markich, S.J.; Brown, P.L.; Jeffree, R.A.; Lim, R.P. Valve movement responses of Velesunio angasi (Bivalvia: Hyriidae) to manganese and uranium: An exception to the free ion activity model. Aquat. Toxicol. 2000, 51, 155–175. [Google Scholar] [CrossRef] [PubMed]
  15. Markich, S.J.; Brown, P.L.; Jeffree, R.A.; Lim, R.P. The effects of pH and dissolved organic carbon on the toxicity of cadmium and copper to a freshwater bivalve: Further support for the extended free ion activity model. Arch. Environ. Contam. Toxicol. 2003, 45, 479–491. [Google Scholar] [CrossRef] [PubMed]
  16. Crémazy, A.; Campbell, P.G.; Fortin, C. The biotic ligand model can successfully predict the uptake of a trivalent ion by a unicellular alga below pH 6.5 but not above: Possible role of hydroxo-species. Environ. Sci. Technol. 2013, 47, 2408–2415. [Google Scholar] [CrossRef] [PubMed]
  17. Wilkinson, K.J.; Campbell, P.G.C.; Couture, P. Effect of fluoride complexation on aluminum toxicity towards juvenile Atlantic salmon (Salmo salar). Can. J. Fish. Aquat. Sci. 1990, 47, 1446–1452. [Google Scholar] [CrossRef]
  18. Crémazy, A.; Campbell, P.G.; Fortin, C. In the presence of fluoride, free Sc3+ is not a good predictor of Sc bioaccumulation by two unicellular algae: Possible role of fluoro-complexes. Environ. Sci. Technol. 2014, 48, 9754–9761. [Google Scholar] [CrossRef]
  19. Koopal, L.K.; van Riemsdijk, W.H.; Kinniburgh, D.G. Humic matter and contaminants. General aspects of modeling metal ion binding. Pure Appl. Chem. 2001, 72, 2005–2016. [Google Scholar] [CrossRef]
  20. Wang, P.; Zhou, D.-M.; Li, L.-Z.; Luo, X.-S. Evaluating the biotic ligand model for toxicity and the alleviation of toxicity in terms of cell membrane surface potential. Environ. Toxicol. Chem. 2010, 29, 1503–1511. [Google Scholar] [CrossRef]
  21. Bánki, O.; Roskov, Y.; Döring, M.; Ower, G.; Hernández Robles, D.R.; Plata Corredor, C.A.; Stjernegaard Jeppesen, T.; Örn, A.; Vandepitte, L.; Pape, T.; et al. Catalogue of Life (Version 2024-08-29); Catalogue of Life: Amsterdam, The Netherlands, 2024. [Google Scholar]
  22. Deleebeeck, N.M.; Muyssen, B.T.; De Laender, F.; Janssen, C.R.; De Schamphelaere, K.A. Comparison of nickel toxicity to cladocerans in soft versus hard surface waters. Aquat. Toxicol. 2007, 84, 223–235. [Google Scholar] [CrossRef]
  23. Warne, M.; Batley, G.E.; van Dam, R.A.; Chapman, J.C.; Fox, D.R.; Hickey, C.W.; Stauber, J.L. Revised Method for Deriving Australian and New Zealand Water Quality Guideline Values for Toxicants; Australian Government Department of Agriculture and Water Resources: Canberra, Australia, 2018. [Google Scholar]
  24. Tipping, E.; Lofts, S.; Sonke, J.E. Humic ion-binding model VII: A revised parameterisation of cation-binding by humic substances. Environ. Chem. 2011, 8, 225–235. [Google Scholar] [CrossRef]
  25. Markich, S.J.; Brown, P.L. Thermochemical Data (log K) for Metals and Metalloids, ASI C4/24; Aquatic Solutions International: Sydney, Australia, 2024. [Google Scholar]
  26. Bryan, S.E.; Tipping, E.; Hamilton-Taylor, J. Comparison and measured and modelled copper binding by natural organic matter in freshwaters. Comp. Biochem. Physiol. Part C Toxicol. Pharmacol. 2002, 133, 37–49. [Google Scholar] [CrossRef]
  27. Unsworth, E.R.; Wanken, K.W.; Zhang, H.; Davison, W.; Black, F.; Buffle, J.; Cao, J.; Cleven, R.; Galceran, J.; Gunkel, P.; et al. Model predictions of metal speciation in freshwaters compared to measurements made by in situ techniques. Environ. Sci. Technol. 2006, 40, 1942–1949. [Google Scholar] [CrossRef] [PubMed]
  28. Van Laer, L.; Smolders, E.; Degryse, F.; Janssen, C.; De Schamphelaere, K.A. Speciation of nickel in surface waters measured with the Donnan membrane technique. Anal. Chim. Acta 2006, 578, 195–202. [Google Scholar] [CrossRef] [PubMed]
  29. Mueller, K.K.; Lofts, S.; Fortin, C.; Campbell, P.G. Trace metal speciation predictions in natural aquatic systems: Incorporation of dissolved organic matter (DOM) spectroscopic quality. Environ. Chem. 2012, 9, 356–368. [Google Scholar] [CrossRef]
  30. Mueller, K.K.; Couture, R.M.; Fortin, C.; Campbell, P.G. Nickel and copper complexation by natural dissolved organic matterTitration of two contrasting lake waters and comparison of measured and modelled free metal ion concentrations. Environ. Chem. 2024, 21, EN23021. [Google Scholar] [CrossRef]
  31. Chen, W.; Guéguen, C.; Smith, D.S.; Galceran, J.; Puy, J.; Companys, E. Metal (Pb, Cd and Zn) binding to diverse organic matter samples and implications for speciation modeling. Environ. Sci. Technol. 2018, 52, 4163–4172. [Google Scholar] [CrossRef]
  32. Wang, P.; Ding, Y.; Liang, Y.; Liu, M.; Lin, X.; Ye, Q.; Shi, Z. Linking molecular composition to proton and copper binding ability of fulvic acid: A theoretical modeling approach based on FT-ICR-MS analysis. Geochim. Cosmochim. Acta 2021, 312, 279–298. [Google Scholar] [CrossRef]
  33. Mebane, C.A. Bioavailability and toxicity models of copper to freshwater life: The state of regulatory science. Environ. Toxicol. Chem. 2023, 42, 2529–2563. [Google Scholar] [CrossRef]
  34. Town, R.M.; van Leeuwen, H.P. Chemodynamic features of nickel(II) and its complexes: Implications for bioavailability in freshwaters. Ecotoxicol. Environ. Saf. 2022, 241, 113840. [Google Scholar] [CrossRef]
  35. Hassler, C.S.; Slaveykova, V.I.; Wilkinson, K.J. Discriminating between intra- and extracellular metals using chemical extractions. Limnol. Oceanogr. Methods 2004, 2, 237–247. [Google Scholar] [CrossRef]
  36. Warner, R.M. Applied Statistics I: Basic Bivariate Techniques, 3rd ed.; Sage Publications: Los Angeles, CA, USA, 2020. [Google Scholar]
  37. Peters, A.; Merrington, G.; Schlekat, C.; De Schampelaere, K.; Stauber, J.; Batley, G.; Harford, A.; van Dam, R.; Pease, C.; Mooney, T.; et al. Validation of the nickel biotic ligand model for locally relevant species in Australian freshwaters. Environ. Toxicol. Chem. 2018, 37, 2566–2574. [Google Scholar] [CrossRef]
  38. Price, G.A.; Stauber, J.L.; Stone, S.; Koppel, D.J.; Holland, A.; Jolley, D. Does toxicity test variability support bioavailability model predictions being within a factor of 2? Environ. Chem. 2022, 19, 177–182. [Google Scholar] [CrossRef]
  39. Meyer, J.S.; Traudt, E.M.; Ranville, J.F. Is the factor-of-2-rule broadly applicable for evaluating the prediction accuracy of metal-toxicity models? Bull. Environ. Contam. Toxicol. 2018, 100, 64–68. [Google Scholar] [CrossRef] [PubMed]
  40. Nys, C.; Janssen, C.R.; Blust, R.; Smolders, E.; De Schamphalaere, K.A. Reproductive toxicity of binary and ternary mixture combinations of nickel, zinc, and lead to Ceriodaphnia dubia is best predicted with the independent action model. Environ. Toxicol. Chem. 2016, 35, 1796–1805. [Google Scholar] [CrossRef]
  41. Schubauer-Berigan, M.K.; Dierkes, J.R.; Monson, P.D.; Ankley, G.T. pH-dependent toxicity of Cd, Cu, Ni, Pb and Zn to Ceriodaphnia dubia, Pimephales promelas, Hyalella azteca and Lumbriculus variegates. Environ. Toxicol. Chem. 1993, 12, 1261–1266. [Google Scholar]
  42. Keithly, J.; Brooker, J.A.; DeForest, D.K.; Wu, B.K.; Brix, K.V. Acute and chronic toxicity of nickel to a cladoceran (Ceriodaphnia dubia) and an amphipod (Hyalella azteca). Environ. Toxicol. Chem. 2004, 23, 691–696. [Google Scholar] [CrossRef] [PubMed]
  43. Santore, R.C.; Croteau, K.; Ryan, A.C.; Schlekat, C.; Middleton, E.; Garman, E.; Hoang, T. A review of water quality factors that affect nickel bioavailability to aquatic organisms: Refinement of the biotic ligand model for nickel in acute and chronic exposures. Environ. Toxicol. Chem. 2021, 40, 2121–2134. [Google Scholar] [CrossRef]
  44. Kozlova, T.; Wood, C.M.; McGeer, J.C. The effect of water chemistry on the acute toxicity of nickel to the cladoceran Daphnia pulex and the development of a biotic ligand model. Aquat. Toxicol. 2009, 91, 221–228. [Google Scholar] [CrossRef]
  45. Clifford, M.; McGeer, J.C. Development of a biotic ligand model to predict the acute toxicity of cadmium to Daphnia pulex. Aquat. Toxicol. 2010, 98, 1–7. [Google Scholar] [CrossRef]
  46. Markich, S.J. Sensitivity of the glochidia (larvae) of freshwater mussels (Bivalvia: Unionida: Hyriidae) to cadmium, cobalt, copper, lead, nickel and zinc: Differences between metals, species and exposure time. Sci. Total Environ. 2017, 601–602, 1427–1436. [Google Scholar] [CrossRef]
  47. Stubblefield, W.A.; Van Genderen, E.; Cardwell, A.S.; Heijerick, D.G.; Janssen, C.R.; De Schamphelaere, K.A.C. Acute and chronic toxicity of cobalt to freshwater organisms: Using a species sensitivity distribution approach to establish international water quality standards. Environ. Toxicol. Chem. 2020, 39, 799–811. [Google Scholar] [CrossRef]
  48. DeForest, D.K.; Santore, R.C.; Ryan, A.C.; Church, B.; Chowdhury, M.J.; Brix, K.V. Development of biotic ligand model–based freshwater aquatic life criteria for lead following US environmental protection agency guidelines. Environ. Toxicol. Chem. 2017, 36, 2695–2973. [Google Scholar] [CrossRef] [PubMed]
  49. Nys, C.; Van Sprang, P.; Lofts, S.; Baken, S.; Delbeke, K.; De Schamphelaere, K. Updated chronic copper biovailability models for invertebrates and algae. Environ. Toxicol. Chem. 2024, 43, 450–467. [Google Scholar] [CrossRef] [PubMed]
  50. Okamoto, A.; Yamamuro, M.; Tatarazako, N. Acute toxicity of 50 metals to Daphnia magna. J. Appl. Toxicol. 2015, 35, 824–830. [Google Scholar] [CrossRef] [PubMed]
  51. Kleinhenz, L.S.; Nugegoda, D.; Trenfield, M.A.; van Dam, R.A.; Humphrey, C.L.; Mooney, T.J.; Harford, A.J. Acute and chronic toxicity of magnesium to the early life stages of two tropical freshwater mussel species. Ecotoxicol. Environ. Saf. 2019, 184, 109638. [Google Scholar] [CrossRef]
  52. Lee, C.E.; Charmantier, G.; Lorin-Nebel, C. Mechanisms of Na+ uptake from freshwater habitats in animals. Front. Physiol. 2022, 13, 1006113. [Google Scholar] [CrossRef]
  53. Carbonaro, R.F.; Di Toro, D.M. Linear free energy relationships for metal-ligand complexation: Monodentate binding to negatively-charged oxygen donor atoms. Geochim. Cosmochim. Acta 2007, 71, 3958–3968. [Google Scholar] [CrossRef]
  54. Bringolf, R.B.; Morris, B.A.; Boese, C.J.; Santore, R.C.; Allen, H.E.; Meyer, J.S. Influence of dissolved organic matter on acute toxicity of zinc to larval fathead minnows (Pimephales promelas). Arch. Environ. Contam. Toxicol. 2006, 51, 438–444. [Google Scholar] [CrossRef]
  55. De Schamphelaere, K.A.; Janssen, C.R. Bioavailability and chronic toxicity of zinc to juvenile rainbow trout (Oncorhynchus mykiss): Comparison with other fish species and development of a biotic ligand model. Environ. Sci. Technol. 2004, 38, 6201–6209. [Google Scholar] [CrossRef]
  56. De Schamphelaere, K.A.; Janssen, C.R. Development and Validation of Biotic Ligand Models for Predicting Chronic Zinc Toxicity to Fish, Daphnids and Algae; Report ZEB-WA-01; Ghent University: Ghent, Belgium, 2004. [Google Scholar]
  57. Heijerick, D.G.; De Schamphelaere, K.A.C.; Van Sprang, P.A.; Janssen, C.R. Development of a chronic zinc biotic ligand model for Daphnia magna. Ecotoxicol. Environ. Saf. 2005, 62, 1–10. [Google Scholar] [CrossRef]
  58. Cooper, N.L.; Bidwell, J.R.; Kumar, A. Toxicity of copper, lead, and zinc mixtures to Ceriodaphnia dubia and Daphnia carinata. Ecotoxicol. Environ. Saf. 2009, 72, 1523–1528. [Google Scholar] [CrossRef]
  59. Van Regenmortel, T.; Berteloot, O.; Janssen, C.R.; De Schamphelaere, K.A.C. Analysing the capacity of the Daphnia magna and Pseudokirchneriella subcapitata bioavailability models to predict chronic zinc toxicity at high pH and low calcium concentrations and formation of a generalised bioavailability model for D. magna. Environ. Toxicol. Chem. 2017, 36, 2781–2798. [Google Scholar] [CrossRef] [PubMed]
  60. Clifford, M.; McGeer, J.C. Development of a biotic ligand model for the acute toxicity of zinc to Daphnia pulex in soft waters. Aquat. Toxicol. 2009, 91, 26–32. [Google Scholar] [CrossRef] [PubMed]
  61. Heijerick, D.G.; De Schamphelaere, K.A.C.; Janssen, C.R. Biotic ligand model development predicting Zn toxicity to the alga Pseudokirchneriella subcapitata: Possibilities and limitations. Comp. Biochem. Physiol. Part C Toxicol. Pharmacol. 2002, 133, 207–218. [Google Scholar] [CrossRef] [PubMed]
  62. Hyne, R.V.; Pablo, F.; Julli, M.; Markich, S.J. Influence of water chemistry on the acute toxicity of copper and zinc to the cladoceran Ceriodaphnia dubia. Environ. Toxicol. Chem. 2005, 24, 1667–1675. [Google Scholar] [CrossRef] [PubMed]
  63. Paulauskis, J.D.; Winner, R.W. Effects of water hardness and humic acid on zinc toxicity to Daphnia magna Straus. Aquat. Toxicol. 1988, 12, 273–290. [Google Scholar] [CrossRef]
  64. Cusimano, R.F.; Brakke, D.F.; Chapman, G.A. Effects of pH on the toxicities of cadmium, copper and zinc to steelhead trout (Salmo gairdneri). Can. J. Fish. Aquat. Sci. 1986, 43, 1497–1503. [Google Scholar] [CrossRef]
  65. Stauber, J.L.; Gadd, J.; Price, G.A.V.; Evans, A.; Holland, A.; Albert, A.; Batley, G.E.; Binet, M.T.; Golding, L.A.; Hickey, C.; et al. Applicability of chronic multiple linear regression models for predicting zinc toxicity in Australian and New Zealand freshwaters. Environ. Toxicol. Chem. 2023, 42, 2614–2629. [Google Scholar] [CrossRef]
  66. Deleebeeck, N.M.E.; De Schamphelaere, K.A.C.; Heijerick, D.G.; Bossuyt, B.T.A.; Janssen, C.R. The acute toxicity of nickel to Daphnia magna: Predictive capability of bioavailability models in artificial and natural waters. Ecotoxicol. Environ. Saf. 2008, 70, 67–78. [Google Scholar] [CrossRef]
  67. Schlekat, C.E.; Van Genderen, E.; De Schamphelaere, K.A.C.; Artunes, P.M.C.; Rogevich, E.C.; Stubblefield, W.A. Cross-species extrapolation of chronic nickel biotic ligand models. Sci. Total Environ. 2010, 408, 6148–6157. [Google Scholar] [CrossRef]
  68. Deleebeeck, N.M.E.; De Schamphelaere, K.A.C.; Janssen, C.R. Effects of Mg2+ and H+ on the toxicity of Ni2+ to the unicellular green alga Pseudokirchneriella subcapitata: Model development and validation with surface waters. Sci. Total Environ. 2009, 407, 1901–1914. [Google Scholar] [CrossRef]
  69. Deleebeeck, N.M.E.; De Schamphelaere, K.A.C.; Janssen, C.R. A bioavailability model predicting the toxicity of nickel to rainbow trout (Oncorhynchus mykiss) and fathead minnow (Pimephales promelas) in synthetic and natural waters. Ecotoxicol. Environ. Saf. 2007, 67, 1–13. [Google Scholar] [CrossRef] [PubMed]
  70. Mano, H.; Shinohara, N. Acute toxicity of nickel to Daphnia magna: Validation of bioavailability models in Japanese rivers. Water Air Soil Pollut. 2020, 231, 459. [Google Scholar] [CrossRef]
  71. Meyer, J.S.; Santore, R.C.; Bobbitt, J.P.; Debrey, L.D.; Boese, C.J.; Paquin, P.R.; Allen, H.E.; Bergman, H.L.; DiToro, D.M. Binding of nickel and copper to fish gills predicts toxicity when water hardness varies, but free-ion activity does not. Environ. Sci. Technol. 1999, 33, 913–916. [Google Scholar] [CrossRef]
  72. Schroeder, J.E. Development of Models for the Prediction of Short-term and Long-term Toxicity to Hyalella azteca from Separate Exposures to Nickel and Cadmium. Ph.D. Thesis, University of Waterloo, Waterloo, ON, Canada, 2008. [Google Scholar]
  73. Nys, C.; Janssen, C.R.; Van Sprang, P.; De Schamphalaere, K.A. The effect of pH on chronic aquatic nickel toxicity is dependent on the pH itself: Extending the chronic nickel bioavailability models. Environ. Toxicol. Chem. 2016, 35, 1097–1106. [Google Scholar] [CrossRef]
  74. Hoang, T.C.; Tomasso, J.R.; Klaine, S.J. Influence of water quality and age on nickel toxicity to fathead minnows (Pimephales promelas). Environ. Toxicol. Chem. 2004, 23, 86–92. [Google Scholar] [CrossRef]
  75. He, J.; Wang, C.; Schlekat, C.E.; Wu, F.; Middleton, E.; Garman, E.; Peters, A. Validation of nickel bioavailability models for algae, invertebrates, and fish in Chinese surface waters. Environ. Toxicol. Chem. 2023, 42, 1257–1265. [Google Scholar] [CrossRef]
  76. Chan, K. The Influence of Calcium and Dissolved Organic Matter on the Acute and Chronic Toxicity of Nickel to Hyalella azteca. Master’s Thesis, Wilfrid Laurier University, Waterloo, ON, Canada, 2013. [Google Scholar]
  77. Niyogi, S.; Kent, R.; Wood, C.M. Effects of water chemistry variables on gill binding and acute toxicity of cadmium in rainbow trout (Oncorhynchus mykiss): A biotic ligand model (BLM) approach. Comp. Biochem. Physiol. Part C Toxicol. Pharmacol. 2008, 148, 305–314. [Google Scholar] [CrossRef]
  78. Tan, Q.-G.; Wang, W.-X. Acute toxicity of cadmium in Daphnia magna under different calcium and pH conditions: Importance of influx rate. Environ. Sci. Technol. 2011, 45, 1970–1976. [Google Scholar] [CrossRef]
  79. Källqvist, T. Effect of water hardness on the toxicity of cadmium to the green alga Pseudokirchneriella subcapitata in an artificial growth medium and nutrient-spiked natural lake waters. J. Toxicol. Environ. Health 2009, 72, 277–283. [Google Scholar] [CrossRef]
  80. Clifford, M.S. A Study of the Waterborne and Dietary Toxicity of Cadmium to Hydra attenuata and Daphnia pulex in Soft Waters and the Development of Biotic Ligand Models to Predict such Toxicity. Ph.D. Thesis, Wilfrid Laurier University, Waterloo, ON, Canada, 2009. [Google Scholar]
  81. Jackson, B.P.; Lasier, P.J.; Miller, W.P.; Winger, P.W. Effects of calcium, magnesium and sodium on alleviating cadmium toxicity to Hyalella azteca. Bull. Environ. Contam. Toxicol. 2000, 64, 279–286. [Google Scholar] [CrossRef]
  82. Marr, J.C.A.; Hansen, J.A.; Meyer, J.S.; Cacela, D.; Podrabsky, T.; Lipton, J.; Bergman, H.L. Toxicity of cobalt and copper to rainbow trout: Application of a mechanistic model for predicting survival. Aquat. Toxicol. 1998, 43, 225–238. [Google Scholar] [CrossRef]
  83. dos Reis, L.L.; de Abreu, C.B.; Gebara, R.C.; Rocha, G.S.; Longo, E.; da Silva Mansano, A.; da Graça Gama Melão, M. Isolated and combined effects of cobalt and nickel on the microalga Raphidocelis subcapitata. Ecotoxicology 2024, 33, 104–118. [Google Scholar] [CrossRef] [PubMed]
  84. Alsop, D.H.; Wood, C.M. Kinetic analysis of zinc accumulation in the gills of juvenile rainbow trout: Effects of zinc acclimation and implications for biotic ligand modelling. Environ. Toxicol. Chem. 2000, 19, 1911–1918. [Google Scholar] [CrossRef]
  85. Long, K.E.; Van Genderen, E.J.; Klaine, S.J. The effects of low hardness and pH on copper toxicity to Daphnia magna. Environ. Toxicol. Chem. 2004, 23, 72–75. [Google Scholar] [CrossRef]
  86. Crémazy, A.; Wood, C.M.; Ng, T.Y.-T.; Smith, D.S.; Chowdhury, M.J. Experimentally derived acute and chronic copper biotic ligand models for rainbow trout. Aquat. Toxicol. 2017, 192, 224–240. [Google Scholar] [CrossRef]
  87. Erickson, R.J.; Benoit, D.A.; Mattson, V.R.; Nelson, H.P.; Leonard, E.N. The effects of water chemistry on the toxicity of copper to fathead minnows. Environ. Toxicol. Chem. 1996, 15, 181–193. [Google Scholar] [CrossRef]
  88. De Schamphelaere, K.A.C.; Bossuyt, B.T.A.; Janssen, C.R. Variability of the protective effect of sodium on the acute toxicity of copper to freshwater cladocerans. Environ. Toxicol. Chem. 2007, 26, 535–542. [Google Scholar] [CrossRef]
  89. De Schamphelaere, K.A.C.; Heijerick, D.G.; Janssen, C.R. Refinement and field validation of a biotic ligand model predicting acute copper toxicity to Daphnia magna. Comp. Biochem. Physiol. Part C Toxicol. Pharmacol. 2002, 133, 243–258. [Google Scholar] [CrossRef]
  90. Ryan, A.C.; Tomasso, J.R.; Klaine, S.J. Influence of pH, hardness, dissolved organic carbon concentration, and dissolved organic matter source on the acute toxicity of copper to Daphnia magna in soft waters: Implications for the biotic ligand model. Environ. Toxicol. Chem. 2009, 28, 1663–1670. [Google Scholar] [CrossRef]
  91. Villavicencio, G.; Urrestarazu, P.; Carvajal, C.; De Schamphelaere, K.A.C.; Janssen, C.R.; Torres, J.C.; Rodriguez, P.H. Biotic ligand model prediction of copper toxicity to Daphnids in a range of natural waters in Chile. Environ. Toxicol. Chem. 2005, 24, 1287–1299. [Google Scholar] [CrossRef]
  92. Welsh, P.G.; Parrott, J.L.; Dixon, D.G.; Hodson, P.V.; Spry, D.J.; Mierle, G. Estimating acute copper toxicity to larval fathead minnow (Pimephales promelas) in soft water from measurements of dissolved organic carbon, calcium and pH. Can. J. Fish. Aquat. Sci. 1996, 53, 1263–1271. [Google Scholar] [CrossRef]
  93. Sciera, K.L.; Isely, J.J.; Tomasso, J.R.; Klaine, S.J. Influence of multiple water-quality characteristics on copper toxicity to fathead minnows (Pimephales promelas). Environ. Toxicol. Chem. 2004, 23, 2900–2905. [Google Scholar] [CrossRef] [PubMed]
  94. Kramer, K.J.M.; Jak, R.G.; van Hattum, B.; Hooftman, R.N.; Zwolsman, J.J.G. Copper toxicity in relation to surface water-dissolved organic matter: Biological effects to Daphnia magna. Environ. Toxicol. Chem. 2004, 23, 2971–2980. [Google Scholar] [CrossRef] [PubMed]
  95. Goulet, R.R.; Thompson, P.A.; Serben, K.C.; Eickhoff, C.V. Impact of environmentally based chemical hardness on uranium speciation and toxicity in six aquatic species. Environ. Toxicol. Chem. 2015, 34, 562–574. [Google Scholar] [CrossRef]
  96. Semaan, M.; Holdway, D.A.; van Dam, R.A. Comparative sensitivity of three populations of the cladoceran Moinodaphnia macleayi to acute and chronic uranium exposure. Environ. Toxicol. 2001, 16, 365–376. [Google Scholar] [CrossRef]
  97. Markich, S.J. Water hardness reduces the accumulation and toxicity of uranium in a freshwater macrophyte (Ceratophyllum demersum). Sci. Total Environ. 2013, 443, 582–589. [Google Scholar] [CrossRef]
  98. Trenfield, M.A.; Ng, J.C.; Noller, B.N.; Markich, S.J.; van Dam, R.A. Dissolved organic carbon reduces uranium bioavailability and toxicity. 2. Uranium(VI) speciation and toxicity to three tropical freshwater organisms. Environ. Sci. Technol. 2011, 45, 3082–3089. [Google Scholar] [CrossRef]
  99. Charles, A.L.; Markich, S.J.; Stauber, J.L.; De Filippis, L.F. The effect of water hardness on the toxicity of uranium to a tropical freshwater alga (Chlorella sp.). Aquat. Toxicol. 2002, 60, 61–73. [Google Scholar] [CrossRef]
  100. Franklin, N.M.; Stauber, J.L.; Markich, S.J.; Lim, R.P. pH-dependent toxicity of copper and uranium to a tropical freshwater alga (Chlorella sp.). Aquat. Toxicol. 2000, 48, 275–289. [Google Scholar] [CrossRef]
  101. De Schamphelaere, K.A.; Nys, C.; Janssen, C.R. Toxicity of lead (Pb) to freshwater green algae: Development and validation of a bioavailability model and inter-species sensitivity comparison. Aquat. Toxicol. 2014, 155, 348–359. [Google Scholar] [CrossRef]
  102. Esbaugh, A.J.; Brix, K.V.; Mager, E.M.; Grosell, M. Multi-linear regression models predict the effects of water chemistry on acute lead toxicity to Ceriodaphnia dubia and Pimephales promelas. Comp. Biochem. Physiol. Part C Toxicol. Pharmacol. 2011, 154, 137–145. [Google Scholar] [CrossRef] [PubMed]
  103. Mager, E.M.; Esbaugh, A.J.; Brix, K.V.; Ryan, A.C.; Grosell, M. Influences of water chemistry on the acute toxicity of lead to Pimephales promelas and Ceriodaphnia dubia. Comp. Biochem. Physiol. Part C Toxicol. Pharmacol. 2011, 153, 82–90. [Google Scholar] [CrossRef] [PubMed]
  104. Antunes, P.M.C.; Kreager, N.J. Lead toxicity to Lemna minor predicted using a metal speciation chemistry approach. Environ. Toxicol. Chem. 2014, 33, 2225–2233. [Google Scholar] [CrossRef] [PubMed]
  105. Nys, C.; Janssen, C.R.; Mager, E.M.; Esbaugh, A.J.; Brix, K.V.; Grosell, M.; Stubblefield, W.A.; Holtze, K.; De Schamphelaere, K.A.C. Development and validation of a biotic ligand model for predicting chronic toxicity of lead to Ceriodaphnia dubia. Environ. Toxicol. Chem. 2014, 33, 394–403. [Google Scholar] [CrossRef]
  106. Grosell, M.; Gerdes, R.; Brix, K.V. Influence of Ca, humic acid and pH on lead accumulation and toxicity in the fathead minnow during prolonged water-borne lead exposure. Comp. Biochem. Physiol. Part C Toxicol. Pharmacol. 2006, 143, 473–483. [Google Scholar] [CrossRef]
  107. Bircneau, O.; Chowdhury, M.J.; Gillis, P.L.; McGeer, J.C.; Wood, C.M.; Wilkie, M.P. Modes of metal toxicity and impaired branchial ionregulation in rainbow trout exposed to mixtures of Pb and Cd in soft water. Aquat. Toxicol. 2008, 89, 222–231. [Google Scholar] [CrossRef]
  108. Macdonald, A.; Silk, L.; Schwartz, M.; Playle, R.C. A lead–gill binding model to predict acute lead toxicity to rainbow trout (Oncorhynchus mykiss). Comp. Biochem. Physiol. Part C Toxicol. Pharmacol. 2002, 133, 227–242. [Google Scholar] [CrossRef]
  109. Liao, C.-M.; Ju, Y.-R.; Chen, W.Y. Subcelllular partitioning links BLM-based toxicokinetics for assessing cadmium toxicity to rainbow trout. Environ. Toxicol. 2010, 26, 600–609. [Google Scholar] [CrossRef]
  110. Altszyler, E.; Ventura, A.C.; Colman-Lerner, A.; Chernomoretz, A. Ultrasensitivity in signaling cascades revisited: Linking local and global ultrasensitivity estimations. PLoS ONE 2017, 12, e0180083. [Google Scholar]
  111. Niyogi, S.; Couture, P.; Pyle, G.; McDonald, D.G.; Wood, C.M. Acute cadmium biotic ligand model characteristics of laboratory-reared and wild yellow perch (Perca flavescens) relative to rainbow trout (Oncorhynchus mykiss). Can. J. Fish. Aquat. Sci. 2004, 61, 942–953. [Google Scholar] [CrossRef]
  112. Taylor, L.N.; Wood, C.M.; McDonald, D.G. An evaluation of sodium loss and gill metal binding properties to explain species differences in copper tolerance. Environ. Toxicol. Chem. 2003, 22, 2159–2166. [Google Scholar] [CrossRef] [PubMed]
  113. Leonard, E.M.; Marentette, J.R.; Balshine, S.; Wood, C.M. Critical body residues, Michaelis-Menten analysis of bioaccumulation, lethality and behaviour as endpoints of waterborne Ni toxicity in two teleosts. Ecotoxicology 2014, 23, 147–162. [Google Scholar] [CrossRef] [PubMed]
  114. Brown, P.L.; Sylva, R.N. Unified Theory of Metal-Ion Complex Formation Constants; Australian Atomic Energy Commission: Sydney, Australia, 1987. [Google Scholar]
  115. Wilde, K.L.; Stauber, J.L.; Markich, S.J.; Franklin, N.M.; Brown, P.L. The effect of pH on the uptake and toxicity of copper and zinc in a tropical freshwater alga (Chlorella sp.). Arch. Environ. Contam. Toxicol. 2006, 51, 174–185. [Google Scholar] [CrossRef] [PubMed]
Figure 1. A schematic overview of the biotic ligand model (BLM) for divalent metals, demonstrating the binding of the free metal ion (M2+) and metal hydroxide (MOH+) to receptors (transporters) on the cell membrane surface of aquatic organisms. These metal ions may enter the cell, potentially inducing toxicity. The metal binding constant (log KMBL) represents the binding affinity of metals with cell surface receptors. The M2+ and MOH+ attached to purple spheres represent extracellular (or surface-bound) metal (not taken up by cells), typically removed by chemical extraction prior to measuring intracellular (internal) metal. The key cations (H+, Ca2+, Mg2+, and Na+) that may ameliorate the binding of M2+ and MOH+ at the cell surface receptor are also shown.
Figure 1. A schematic overview of the biotic ligand model (BLM) for divalent metals, demonstrating the binding of the free metal ion (M2+) and metal hydroxide (MOH+) to receptors (transporters) on the cell membrane surface of aquatic organisms. These metal ions may enter the cell, potentially inducing toxicity. The metal binding constant (log KMBL) represents the binding affinity of metals with cell surface receptors. The M2+ and MOH+ attached to purple spheres represent extracellular (or surface-bound) metal (not taken up by cells), typically removed by chemical extraction prior to measuring intracellular (internal) metal. The key cations (H+, Ca2+, Mg2+, and Na+) that may ameliorate the binding of M2+ and MOH+ at the cell surface receptor are also shown.
Water 16 02999 g001
Figure 2. An illustration of the “cascade of effects” leading to metal toxicity in aquatic organisms. Each stage of the cascade highlights the appearance of the terms {M–X–cell} (metal bound to a receptor site on the cell surface), [M] (free metal ion) and (fC)M50% (the fraction of receptor sites bound by a metal at its median toxic effect concentration (TR50)), in both the numerator and denominator on the right-hand side of the equation.
Figure 2. An illustration of the “cascade of effects” leading to metal toxicity in aquatic organisms. Each stage of the cascade highlights the appearance of the terms {M–X–cell} (metal bound to a receptor site on the cell surface), [M] (free metal ion) and (fC)M50% (the fraction of receptor sites bound by a metal at its median toxic effect concentration (TR50)), in both the numerator and denominator on the right-hand side of the equation.
Water 16 02999 g002
Figure 3. The predictive capacity of the biotic ligand model for (a) zinc (Zn) with seven freshwater species (13 studies) and 142 water quality scenarios (all values expressed as dissolved Zn in mg/L), (b) nickel (Ni) with ten freshwater species (14 studies) and 256 water quality scenarios (all values expressed as dissolved Ni in mg/L), (c) cadmium (Cd) with six freshwater species (seven studies) and 83 water quality scenarios (all values expressed as dissolved Cd in mg/L) and (d) cobalt (Co) with nine freshwater species (four studies) and 12 water quality scenarios (all values expressed as dissolved Co in mg/L). The solid lines indicate ideal (1:1) agreement between the measured and predicted EC/LC50 values, the dashed lines represent ratios of ½ and 2 between the measured and predicted EC/LC50 values and the dotted and dashed lines represent ratios of ⅓ and 3. Within the legend, the values in parentheses denote the number of data points from each study. References: Nys et al. (2016) [40], Keithly et al. (2004) [42], Kozlova et al. (2009) [44], Clifford and McGeer (2010) [45], Stubblefield et al. (2020) [47], Bringolf et al. (2006) [54], De Schamphelaere and Janssen (2004) [55], De Schamphelaere and Janssen (2004) [56], Heijerick et al. (2005) [57], Cooper et al. (2009) [58], Van Regenmortal et al. (2017) [59], Clifford and McGeer (2009) [60], Heijerick et al. (2002) [61], Hyne et al. (2005) [62], Paulauskis and Winner (1988) [63], Cusimano et al. (1986) [64], Stauber et al. (2023) [65], Deleebeek et al. (2008) [66], Schlekat et al. (2010) [67], Deleebeek et al. (2009) [68], Deleebeek et al. (2007) [69], Mano and Shinohara (2020) [70], Meyer et al. (1999) [71], Schroeder (2008) [72], Nys et al. (2016) [73], Hoang et al. (2004) [74], He et al. (2023) [75], Chan (2013) [76], Niyogi et al. (2008) [77], Tan and Wang (2011) [78], Kallqvist (2009) [79], Clifford (2009) [80], Jackson et al. (2000) [81], Marr et al. (1998) [82], dos Reis et al. (2024) [83] and Alsop and Wood (2000) [84].
Figure 3. The predictive capacity of the biotic ligand model for (a) zinc (Zn) with seven freshwater species (13 studies) and 142 water quality scenarios (all values expressed as dissolved Zn in mg/L), (b) nickel (Ni) with ten freshwater species (14 studies) and 256 water quality scenarios (all values expressed as dissolved Ni in mg/L), (c) cadmium (Cd) with six freshwater species (seven studies) and 83 water quality scenarios (all values expressed as dissolved Cd in mg/L) and (d) cobalt (Co) with nine freshwater species (four studies) and 12 water quality scenarios (all values expressed as dissolved Co in mg/L). The solid lines indicate ideal (1:1) agreement between the measured and predicted EC/LC50 values, the dashed lines represent ratios of ½ and 2 between the measured and predicted EC/LC50 values and the dotted and dashed lines represent ratios of ⅓ and 3. Within the legend, the values in parentheses denote the number of data points from each study. References: Nys et al. (2016) [40], Keithly et al. (2004) [42], Kozlova et al. (2009) [44], Clifford and McGeer (2010) [45], Stubblefield et al. (2020) [47], Bringolf et al. (2006) [54], De Schamphelaere and Janssen (2004) [55], De Schamphelaere and Janssen (2004) [56], Heijerick et al. (2005) [57], Cooper et al. (2009) [58], Van Regenmortal et al. (2017) [59], Clifford and McGeer (2009) [60], Heijerick et al. (2002) [61], Hyne et al. (2005) [62], Paulauskis and Winner (1988) [63], Cusimano et al. (1986) [64], Stauber et al. (2023) [65], Deleebeek et al. (2008) [66], Schlekat et al. (2010) [67], Deleebeek et al. (2009) [68], Deleebeek et al. (2007) [69], Mano and Shinohara (2020) [70], Meyer et al. (1999) [71], Schroeder (2008) [72], Nys et al. (2016) [73], Hoang et al. (2004) [74], He et al. (2023) [75], Chan (2013) [76], Niyogi et al. (2008) [77], Tan and Wang (2011) [78], Kallqvist (2009) [79], Clifford (2009) [80], Jackson et al. (2000) [81], Marr et al. (1998) [82], dos Reis et al. (2024) [83] and Alsop and Wood (2000) [84].
Water 16 02999 g003aWater 16 02999 g003bWater 16 02999 g003c
Figure 4. The predictive capacity of the biotic ligand model for (a) copper (Cu) with 12 freshwater species (15 studies) and 383 water quality scenarios (all values expressed as dissolved Cu in μg/L), (b) uranium (UO2) with ten freshwater species (six studies) and 52 water quality scenarios (all values expressed as dissolved U in μg/L) and (c) lead (Pb) with five freshwater species (eight studies) and 110 water quality scenarios (all values expressed as dissolved Pb in μg/L). The solid line indicates ideal (1:1) agreement between the measured and predicted EC/LC50 values, the dashed lines represent ratios of ½ and 2 between the measured and predicted EC/LC50 values and the dotted and dashed lines represent ratios of ⅓ and 3. Within the legend, the values in parentheses denote the number of data points from each study. References: Di Toro et al. (2001) [5], De Schampheleare and Janssen (2002) [9], Meyer et al. (1999) [39], Cooper et al. (2009) [58], Cusimano et al. (1986) [64], Long et al. (2004) [85], Crémazy et al. (2017) [86], Erickson et al. (1996) [87], De Schampheleare et al. (2007) [88], De Schampheleare et al. (2002) [89], Ryan et al. (2009) [90], Villavicencio et al. (2005) [91], Welsh et al. (1996) [92], Sciera et al. (2004) [93], Kramer et al. (2004) [94], Goulet et al. (2015) [95], Semaan et al. (2001) [96], Markich (2013) [97], Trenfield et al. (2011) [98], Charles et al. (2002) [99], Franklin et al. (2001) [100], De Schampheleare et al. (2014) [101], Esbaugh et al. (2011) [102], Mager et al. (2011) [103], Antunes and Kreager (2014) [104], Nys et al. (2014) [105], Grosell et al. (2006) [106], Bircneau et al. (2008) [107] and McDonald et al. (2002) [108].
Figure 4. The predictive capacity of the biotic ligand model for (a) copper (Cu) with 12 freshwater species (15 studies) and 383 water quality scenarios (all values expressed as dissolved Cu in μg/L), (b) uranium (UO2) with ten freshwater species (six studies) and 52 water quality scenarios (all values expressed as dissolved U in μg/L) and (c) lead (Pb) with five freshwater species (eight studies) and 110 water quality scenarios (all values expressed as dissolved Pb in μg/L). The solid line indicates ideal (1:1) agreement between the measured and predicted EC/LC50 values, the dashed lines represent ratios of ½ and 2 between the measured and predicted EC/LC50 values and the dotted and dashed lines represent ratios of ⅓ and 3. Within the legend, the values in parentheses denote the number of data points from each study. References: Di Toro et al. (2001) [5], De Schampheleare and Janssen (2002) [9], Meyer et al. (1999) [39], Cooper et al. (2009) [58], Cusimano et al. (1986) [64], Long et al. (2004) [85], Crémazy et al. (2017) [86], Erickson et al. (1996) [87], De Schampheleare et al. (2007) [88], De Schampheleare et al. (2002) [89], Ryan et al. (2009) [90], Villavicencio et al. (2005) [91], Welsh et al. (1996) [92], Sciera et al. (2004) [93], Kramer et al. (2004) [94], Goulet et al. (2015) [95], Semaan et al. (2001) [96], Markich (2013) [97], Trenfield et al. (2011) [98], Charles et al. (2002) [99], Franklin et al. (2001) [100], De Schampheleare et al. (2014) [101], Esbaugh et al. (2011) [102], Mager et al. (2011) [103], Antunes and Kreager (2014) [104], Nys et al. (2014) [105], Grosell et al. (2006) [106], Bircneau et al. (2008) [107] and McDonald et al. (2002) [108].
Water 16 02999 g004aWater 16 02999 g004b
Figure 5. The predictive capacity of the biotic ligand model for uranium (UO2) with ten freshwater species (six studies) and 52 water quality scenarios (all values expressed as dissolved U in μg/L) when only UO22+ is considered responsible for toxicity. The solid line indicates ideal (1:1) agreement between the measured and predicted EC/LC50 values, the dashed lines represent ratios of ½ and 2 between the measured and predicted EC/LC50 values and the dotted and dashed lines represent ratios of ⅓ and 3. Within the legend, the values in parentheses denote the number of data points from each study. References: Goulet et al. (2015) [95], Semaan et al. (2001) [96], Markich (2013) [97], Trenfield et al. (2011) [98], Charles et al. (2002) [99] and Franklin et al. (2001) [100].
Figure 5. The predictive capacity of the biotic ligand model for uranium (UO2) with ten freshwater species (six studies) and 52 water quality scenarios (all values expressed as dissolved U in μg/L) when only UO22+ is considered responsible for toxicity. The solid line indicates ideal (1:1) agreement between the measured and predicted EC/LC50 values, the dashed lines represent ratios of ½ and 2 between the measured and predicted EC/LC50 values and the dotted and dashed lines represent ratios of ⅓ and 3. Within the legend, the values in parentheses denote the number of data points from each study. References: Goulet et al. (2015) [95], Semaan et al. (2001) [96], Markich (2013) [97], Trenfield et al. (2011) [98], Charles et al. (2002) [99] and Franklin et al. (2001) [100].
Water 16 02999 g005
Figure 6. Linear regressions between (fC)M50% and their corresponding EC/LC50 values for nickel (Ni–green) and copper (Cu–blue). The red dashed lines represent the expected behaviour of Zn (red stars), despite no significant (p > 0.05) linear relationship, and are roughly parallel to the fitted regressions for Ni and Cu (aligning with both the upper and lower observed Zn values).
Figure 6. Linear regressions between (fC)M50% and their corresponding EC/LC50 values for nickel (Ni–green) and copper (Cu–blue). The red dashed lines represent the expected behaviour of Zn (red stars), despite no significant (p > 0.05) linear relationship, and are roughly parallel to the fitted regressions for Ni and Cu (aligning with both the upper and lower observed Zn values).
Water 16 02999 g006
Figure 7. Linear regressions between log KMBL (metal binding affinity) and log (1/[M]) (a measure of metal toxicity) for (a) Ceriodaphnia dubia and Daphnia magna (crustaceans), and (b) Oncorhynchus mykiss (fish). For O. mykiss, Cd2+ values were added as an overlay using both the uncorrected (0.216; open red circle) and corrected (i.e., × (1 − 0.216)/0.216); open purple circle) (fC)M50% value. The correction brings Cd2+ in line with the other six metals (see Section 5).
Figure 7. Linear regressions between log KMBL (metal binding affinity) and log (1/[M]) (a measure of metal toxicity) for (a) Ceriodaphnia dubia and Daphnia magna (crustaceans), and (b) Oncorhynchus mykiss (fish). For O. mykiss, Cd2+ values were added as an overlay using both the uncorrected (0.216; open red circle) and corrected (i.e., × (1 − 0.216)/0.216); open purple circle) (fC)M50% value. The correction brings Cd2+ in line with the other six metals (see Section 5).
Water 16 02999 g007
Table 1. Binding constants (log KMBL) for zinc (Zn2+), nickel (Ni2+), cadmium (Cd2+), cobalt (Co2+), copper (Cu2+), uranium (UO22+) and lead (Pb2+) at cell surface receptor sites (biotic ligands) in freshwater organisms.
Table 1. Binding constants (log KMBL) for zinc (Zn2+), nickel (Ni2+), cadmium (Cd2+), cobalt (Co2+), copper (Cu2+), uranium (UO22+) and lead (Pb2+) at cell surface receptor sites (biotic ligands) in freshwater organisms.
TaxaOrganism alog KMBL
Zn2+Ni2+Cd2+Co2+Cu2+UO22+Pb2+
Fish cAcipenser transmontanus6.1 7.3 8.2 7.3
Cottus bairdii6.1 (3) b 7.2 8.1
Ctenopharyngodon idella5.7
Danio rerio 7.7
Mogurnda mogurnda6.1 7.37.4
Neogobius melanostomus 4.7
Oncorhynchus clarkii5.8 7.8
Oncorhynchus kisutch 7.8
Oncorhynchus mykiss5.9 (13)4.7 (3)7.5 (13)5.3 (3)7.9 (12)7.37.1 (3)
Oncorhynchus nerka 7.7
Oncorhynchus tshawytscha6.1 7.6 7.7 (2)
Oreochromis mossambicus 7.0
Perca flavescens 7.2
Pimephales promelas5.7 (3)4.7 (3)7.1 (3)5.77.6 (11)7.46.8 (5)
Prosopium williamsoni5.5
Salmo trutta 7.6
Salvelinus confluentus5.7 7.2
Average (fish)5.94.77.35.57.87.47.1
Crustaceans
(cladocerans)
Acantholeberis curvirostris 7.8
Acroperus elongatus 8.2
Acroperus harpae 8.3
Alona sp. 8.1
Alona affinis 4.7
Alona quadrangularis 8.0
Bosmina coregoni 5.0
Bosmina longirostris 7.9 (2)
Camptocercus lilljeborgi 4.9
Ceriodaphnia cornuta 7.9
Ceriodaphnia dubia5.9 (8)5.2 (4)7.0 (3)5.58.0 (15)7.47.0 (7)
Ceriodaphnia pulchella 4.8 7.7 (2)
Ceriodaphnia quadrangula 4.8
Crustaceans
(cladocerans)
Ceriodaphnia reticulata6.2 6.5 8.1 (3)
Ceriodaphnia rigaudi 7.0
Ceriodaphnia silvestri 6.9
Chydorus ovalis 4.7 7.9
Chydorus sphearicus 8.1 (2)
Daphnia ambigua 7.9
Daphnia carinata5.7 7.4 6.9 (2)
Daphnia exilis 7.8
Daphnia galeata 7.6 (3)
Daphnia longispina 4.9 7.8 (2)
Daphnia magna5.8 (9)4.9 (6)6.9 (7)5.4 (2)7.8 (13)8.06.7 (4)
Daphnia obtusa 7.8
Daphnia pulex5.9 (2)4.7 (2)6.6 (2) 7.9 (4)
Daphnia similis 6.8
Daphnia thomsoni6.3
Disparalona rostrata 7.7
Eurycercus lamellatus 8.1
Moina affinis6.2 7.0 8.0
Moinodaphnia macleayi6.3 7.78.2
Peracantha truncata 4.6
Pleuroxus truncatus 7.7
Pseudosida ramosa 7.2
Pseudosida variabilis 8.3
Scapholeberis microcephala 8.2
Scapholeberis mucronata 8.1
Simocephalus expinosus 7.8 (2)
Simocephalus serrulatus 5.06.9
Simocephalus vetulus 4.86.9 7.9 (3)
Average (cladocerans)6.04.86.95.57.97.96.8
Crustaceans
(amphipods)
Gammarus sp. 7.5
Gammarus pseudolimnaeus 6.5
Gammarus pulex 4.7
Hyalella azteca5.9 (2)5.1 (4)7.3 (2)5.87.87.46.9
Crustaceans
(copepod)
Notodiaptomus iheringi 7.4 7.5
Crustaceans
(ostracod)
Cypridopsis vidua 7.8
Crustaceans
(decapod)
Macrobrachium rosenbergi5.6
Average (crustaceans)6.04.96.95.67.97.76.9
Molluscs d
(bivalves)
Alathyria profuga5.94.86.65.97.4 6.4
Cambarunio iris 8.1 (3)
Corbicula fluminea 7.77.6
Cucumerunio novaehollandiae6.05.06.86.17.7 6.6
Echyridella menziesii5.7 8.1
Epioblasma rangiana 7.9
Epioblasma triquetra 8.2
Hyridella australis6.05.06.86.17.6 6.6
Hyridella depressa6.04.96.76.07.5 6.5
Molluscs
(bivalves)
Hyridella drapeta5.94.96.66.07.5 6.5
Lampsilis abrupta 7.9
Lampsilis fasciola 8.0 (2)
Lampsilis rafinesqueana6.3 7.4 7.8 6.8
Lampsilis siliquoidea6.3 7.5 7.8 (3) 6.7
Obovaria subrotunda 7.8
Ortmanniana ligamentina 7.6 (2)
Paetulunio fabalis 8.2
Potamilus ohiensis 7.9
Velesunio sp. 8.2
Velesunio anbiguus5.84.86.65.97.4 6.4
Velesunio angasi6.3
Venustaconcha ellipsiformis 8.1
Average (bivalves)6.04.96.96.07.87.66.6
Molluscs
(gastropods)
Amerianna cumingi6.2 7.98.1
Filopaludina bengalensis 7.9
Fluminicola sp. 8.2
Fontigens aldrichi 7.9
Lymnaea stagnalis 4.9 (3)6.6 (2)5.98.0 (4)
Physella acuta 8.3
Physella gyrina 7.7
Planorbella pilsbryi 8.2
Pomacea paludosa5.8
Potamopyrgus antipodarum 7.3
Pyrgulopsis robusta 8.0
Racesina luteola 6.5 8.2
Taylorconcha serpenticola 8.1
Average (gastropods)6.04.96.85.98.08.1
Average (molluscs)6.04.96.86.07.97.96.6
RotifersAnuraeopsis fissa5.8
Brachionus calicyflorus 4.6 7.4
Brachionus rubens5.5
Euchlanis dilatata6.1 7.9
Lecane inermis6.3 8.0
Lecane quadridentata6.2 7.5
Average (rotifers)6.04.68.07.4
Microalgae
(Chlorophyceae)
Ankistrodesmus arcuatus 4.6 8.0
Chlamydomonas reinhardtii 7.7 (2)
Chlorella sp.5.8 (4)5.0 (2) 7.8 (11)7.8 (4)
Raphidocelis subcapitata6.1 (9)4.7 (6)7.0 (7)5.7 (2)7.8 (8)8.17.0 (2)
Average (microalgae)6.04.87.05.77.97.97.0
MacrophytesCeratophyllum demersum 7.27.7
Lemna aequinoctialis 7.77.4 (3)
Lemna minor6.24.7 (3)7.16.1 (2)7.47.8 (2)7.2
HydraHydra circumcincta 7.1
Hydra viridissima6.15.2 7.7 (2)7.3 (3)
InsectsChironomus tentans 7.6
Tanytarsus dissimilis6.1 7.3 7.5 6.6
EuglenidEuglena gracilis 7.7
AnnelidAeolosoma sp. 5.4
BacteriaErwinia sp. 7.7
Average (all)6.04.87.05.87.97.76.8
95% confidence limit0.10.10.10.20.10.10.1
Number of organisms (n)42313715842021
Note(s): a Taxonomic classification as per the Catalogue of Life [21]. b Values with parentheses (n) are averages (individual organism log KMBL values and references are provided in Table S1). c Fish are all teleostei, except for Acipenser transmontanus (chondrostei). d Bivalves are all unionids, except for Corbicula fluminea (venerid).
Table 2. Binding constants (log KMBL) for the ameliorative cations (calcium (Ca2+), magnesium (Mg2+), sodium (Na+) and protons (H+)) at cell surface receptor sites (biotic ligands) in freshwater organisms.
Table 2. Binding constants (log KMBL) for the ameliorative cations (calcium (Ca2+), magnesium (Mg2+), sodium (Na+) and protons (H+)) at cell surface receptor sites (biotic ligands) in freshwater organisms.
TaxaOrganism alog KMBL
Ca2+Mg2+Na+H+
FishCtenopharyngodon idella 2.8
Mogurnda mogurnda 5.8
Oncorhynchus mykiss3.5 (7) b3.02.45.8 (5)
Oreochromis mossambicus3.5
Perca flavescens3.7
Pimephales promelas3.2 (4) 5.6 (2)
Average (fish)3.52.92.45.7
Crustaceans
(cladocerans)
Ceriodaphnia dubia3.2 (2)2.62.4 (2)5.6
Ceriodaphnia pulchella 2.4
Ceriodaphnia reticulata 2.1
Daphnia galeata 2.2
Daphnia longispina 2.1
Daphnia magna3.2 (5)2.8 (6)2.2 (5)5.8 (3)
Daphnia pulex3.5 (2)3.0 (2)2.15.7
Simocephalus expinosus 2.0
Simocephalus vetulus 2.4
Ten species of cladocerans c3.22.6
Average (cladocerans)3.32.72.25.7
Crustaceans
(amphipod)
Hyalella azteca3.4 (3)2.8 (2) 5.8
Average (crustaceans)3.32.82.25.7
Molluscs
(bivalves)
Corbicula fluminea3.5
Velesunio angasi 2.7
Average (molluscs)3.52.7
MicoalgaeChlamydomonas reinhardtii3.72.8 6.2
Chlorella sp. 5.8 (2)
Raphidocelis subcapitata3.12.9 (2)2.25.7 (3)
Average (microalgae)3.42.82.25.9
MacrophyteCeratophyllum demersum3.22.7
HydraHydra circumcincta3.72.8
Average3.42.82.25.8
95% confidence limit0.10.10.10.1
Number of organisms (n)14121110
Note(s): a Taxonomic classification as per Catalogue of Life [21]. b Values with parentheses (n) are averages (individual organism log KMBL values and references are provided in Table S2). c Group average calculated by Deleebeeck et al. [22].
Table 3. A comparison of measured versus predicted metal toxicity using the conceptual model.
Table 3. A comparison of measured versus predicted metal toxicity using the conceptual model.
Metal% of Data within a
Factor of Two of the 1:1 Agreement Line
% of Data within a
Factor of Three of the
1:1 Agreement Line
Number of Water Quality ScenariosNumber of
Organisms
(Species)
Number of
Studies
Zinc8194142713
Nickel83972561014
Cadmium88988367
Cobalt75921294
Copper79913831215
Uranium758852106
Lead809511058
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Brown, P.L.; Markich, S.J. An Evaluation of Metal Binding Constants to Cell Surface Receptors in Freshwater Organisms, and Their Application in Biotic Ligand Models to Predict Metal Toxicity. Water 2024, 16, 2999. https://doi.org/10.3390/w16202999

AMA Style

Brown PL, Markich SJ. An Evaluation of Metal Binding Constants to Cell Surface Receptors in Freshwater Organisms, and Their Application in Biotic Ligand Models to Predict Metal Toxicity. Water. 2024; 16(20):2999. https://doi.org/10.3390/w16202999

Chicago/Turabian Style

Brown, Paul L., and Scott J. Markich. 2024. "An Evaluation of Metal Binding Constants to Cell Surface Receptors in Freshwater Organisms, and Their Application in Biotic Ligand Models to Predict Metal Toxicity" Water 16, no. 20: 2999. https://doi.org/10.3390/w16202999

APA Style

Brown, P. L., & Markich, S. J. (2024). An Evaluation of Metal Binding Constants to Cell Surface Receptors in Freshwater Organisms, and Their Application in Biotic Ligand Models to Predict Metal Toxicity. Water, 16(20), 2999. https://doi.org/10.3390/w16202999

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