Abstract:Natural organic matter (NOM, or humic substance) has a known tendency to form colloidal aggregates in aqueous environments, with the composition and concentration of cationic species in solution, pH, temperature, and the composition of the NOM itself playing important roles. Strong interaction of carboxylic groups of NOM with dissolved metal cations is thought to be the leading chemical interaction in NOM supramolecular aggregation. Computational molecular dynamics (MD) study of the interactions of Na(+), Mg(2+), and Ca(2+) with the carboxylic groups of a model NOM fragment and acetate anions in aqueous solutions provides new quantitative insight into the structure, energetics, and dynamics of the interactions of carboxylic groups with metal cations, their association, and the effects of cations on the colloidal aggregation of NOM molecules. Potentials of mean force and the equilibrium constants describing overall ion association and the distribution of metal cations between contact ion pairs and solvent-separated ions pairs were computed from free MD simulations and restrained umbrella sampling calculations. The results provide insight into the local structural environments of metal-carboxylate association and the dynamics of exchange among these sites. All three cations prefer contact ion pair to solvent-separated ion pair coordination, and Na(+) and Ca(2+) show a strong preference for bidentate contact ion pair formation. The average residence time of a Ca(2+) ion in a contact ion pair with the carboxylic groups is of the order of 0.5 ns, whereas the corresponding residence time of a Na(+) ion is only between 0.02 and 0.05 ns. The average residence times of a Ca(2+) ion in a bidentate coordinated contact ion pair vs a monodentate coordinated contact ion pair are about 0.5 and 0.08 ns, respectively. On the 10 ns time scale of our simulations, aggregation of the NOM molecules occurs in the presence of Ca(2+) but not Na(+) or Mg(2+). These results agree with previous experimental observations and are explained by both Ca(2+) ion bridging between NOM molecules and decreased repulsion between the NOM molecules due to the reduced net charge of the NOM-metal complexes. Simulations on a larger scale are needed to further explore the relative importance of the different aggregation mechanisms and the stability of NOM aggregates.

Binding of Biologically Relevant Divalent Cations to Aqueous Carboxylates: Molecular Simulations Guided by Raman Spectroscopy

Binding of Divalent Cations to Aqueous Acetate: Molecular Simulations Guided by Raman Spectroscopy

Association Free Energies of Metal Cations with Mesylate and Acetate in Brine Calculated via Molecular Dynamics Simulation

Theoretical Binding Affinities And Spectra Of Complexes Formed By A Cyclic Beta-Peptoid With Amino Acids

Thermodynamics of Metal-Acetate Interactions

Ab Initio Molecular Dynamics Study on the Interactions between Carboxylate Ions and Metal Ions in Water

Metal cation complexation with natural organic matter in aqueous solutions: molecular dynamics simulations and potentials of mean force.

Quantum effects in cation interactions with first and second coordination shell ligands in metalloproteins

Interaction of divalent cations and amino acids in bulk water: Molecular simulations with neural network potentials

Study of interactions between metal ions and protein model compounds by energy decomposition analyses and the AMOEBA force field.

Molecular dynamics simulation of cationic complexation with natural organic matter

Simulating Metal-Imidazole Complexes

Binding of Phosphate Species to Ca 2+ and Mg 2+ in Aqueous Solution

Binding thermodynamics of divalent metal ions to several biological buffers

Metal cation dependence of interactions with amino acids: bond dissociation energies of Rb(+) and Cs(+) to the acidic amino acids and their amide derivatives

A Qicar Approach for Quantifying Binding Constants for Metal-Ligand Complexes

Conformations and Metal Ion Affinities of Glutamine Binding with Alkali and Alkaline Earth Metal Cations: an Ab Initio Study

[Measurement of Binding Constants of Metal Complex by Indirect Ultraviolet Affinity Capillary Electrophoresis].

Ambivalent Metal-Ion Binding-Properties of Cytidine in Aqueous-Solution

Coarse-Grained Parameters for Divalent Cations within the SIRAH Force Field

Interaction Structure and Affinity of Zwitterionic Amino Acids with Important Metal Cations (Cd<SUP>2+</SUP>, Cu<SUP>2+</SUP>, Fe<SUP>3+</SUP>, Hg<SUP>2+</SUP>, Mn<SUP>2+</SUP>, Ni<SUP>2+</SUP> and Zn<SUP>2+</SUP>) in Aqueous Solution: A Theoretical Study