Molecular Dynamics Simulation of the Structures, Dynamics, and Aggregation of Dissolved Organic Matter

Deepa Devarajan,Liyuan Liang,Baohua Gu,Scott C. Brooks,Jerry M. Parks,Jeremy C. Smith
DOI: https://doi.org/10.1021/acs.est.0c01176
2020-09-28
Abstract:Dissolved organic matter (DOM) plays a significant role in the transport and transformation of pollutants in the aquatic environment. However, the experimental characterization of DOM has been limited mainly to bulk properties, and the molecular-level interactions among various components of DOM remain to be fully characterized. Here, we use molecular dynamics (MD) simulations to probe the structural properties of model DOM systems at atomic detail. The 200 ns simulations, validated by available experimental data, reveal processes and mechanisms by which chemical species (cations, peptides, lipids, lignin, carbohydrates, and some low-molecular-weight aliphatic and aromatic compounds) aggregate to form complex DOM. The DOM aggregates are dynamic, consisting of a hydrophobic core and amphiphilic exterior. The lipid tails and other hydrophobic fragments form the core, with hydrophilic and amphiphilic groups exposed to water, making DOM accessible to both polar and nonpolar species. Thus, the lipid component acts as a nucleator, whereas cations (especially Ca<sup>2+</sup>) connect the molecular fragments on the surface by coordinating with the O-containing functional groups of DOM. The structural details revealed here provide new insights including surface accessible atoms, overall assemblage, and interactions among the molecules of DOM for understanding the kinetics and mechanisms through which DOM interacts with metal and other contaminants.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.est.0c01176?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.est.0c01176</a>.Component molecules; van Krevelen diagram; construction of DOM models; p<i>K</i><sub>a</sub> values; DOM1–7 compositions; percentage distribution of O and N; calculated <i>R</i><sub>g</sub>; method details; distributions of the standard deviations of contacts; RMSD of contacts; radial distribution function; representative snapshots; PDF; interactions of component molecules in DOM6; contact probabilities; modeling results of clusters; calculated <i>D</i>; calculated MSD; linear correlation of the calculated <i>D</i>; molecular and atomic contributions to SASA; effect of other cations on DOM bulk properties; contact maps for DOM6a and DOM7a; number of hydrogen bonds per heteroatom; and calculated bulk properties of the DOM models (<a class="ext-link" href="/doi/suppl/10.1021/acs.est.0c01176/suppl_file/es0c01176_si_001.pdf">PDF</a>)Movies of the MD trajectories for DOM6 (<a class="ext-link" href="/doi/suppl/10.1021/acs.est.0c01176/suppl_file/es0c01176_si_002.mp4">MP4</a>)Movies of the MD trajectories for DOM7 (<a class="ext-link" href="/doi/suppl/10.1021/acs.est.0c01176/suppl_file/es0c01176_si_003.mp4">MP4</a>)Coordinates of representative snapshots of the equilibrated model systems (DOM6) (<a class="ext-link" href="/doi/suppl/10.1021/acs.est.0c01176/suppl_file/es0c01176_si_004.pdb">PDB</a>)Coordinates of representative snapshots of the equilibrated model systems (DOM7) (<a class="ext-link" href="/doi/suppl/10.1021/acs.est.0c01176/suppl_file/es0c01176_si_005.pdb">PDB</a>)This article has not yet been cited by other publications.
environmental sciences,engineering, environmental
What problem does this paper attempt to address?