Models of the hydrogenoxalate (bioxalate, charge -1) and oxalate (charge -2) anions were developed for classical molecular dynamics (CMD) simulations and parametrized against ab initio molecular dynamics (AIMD) data from our previous study (Kroutil et al. (2016) J Mol Model 22:210). The interactions of the anions with water were described using charges scaled according to the electronic continuum correction approach with rescaling of nonbonded parameters (ECCR), and those descriptions of anion interactions were found to agree well with relevant AIMD and experimental results. The models with full RESP charges showed excessively strong electrostatic interactions between the solute and water molecules, leading to an overstructured solvation shell around the anions and thus to a diffusion coefficient that was much too low. The effect of charge scaling was more evident for the oxalate dianion than for the hydrogenoxalate anion. Our work provides CMD models for ions of oxalic acid and extends previous studies that showed the importance of ECCR for modeling divalent ions and ions of organic compounds. Graphical abstract The radial distribution function between a water oxygen (Ow) and an oxygen of the oxalate dianion (Ox) significantly improved when scaled charges were applied to the anion.

• Klíčová slova
• AIMD, Ab initio molecular dynamics, CMD, Classical molecular dynamics, ECCR, Electronic continuum correction, Hydrogenoxalate, Oxalate, Oxalic acid anions,
• Publikační typ
• časopisecké články MeSH

A combination of Fourier transform infrared and phase transition measurements as well as molecular computer simulations, and thermodynamic modeling were performed to probe the mechanisms by which guanidinium (Gnd+) salts influence the stability of the collapsed versus uncollapsed state of an elastin-like polypeptide (ELP), an uncharged thermoresponsive polymer. We found that the cation's action was highly dependent upon the counteranion with which it was paired. Specifically, Gnd+ was depleted from the ELP/water interface and was found to stabilize the collapsed state of the macromolecule when paired with well-hydrated anions such as SO42-. Stabilization in this case occurred via an excluded volume (or depletion) effect, whereby SO42- was strongly partitioned away from the ELP/water interface. Intriguingly, at low salt concentrations, Gnd+ was also found to stabilize the collapsed state of the ELP when paired with SCN-, which is a strong binder for the ELP. In this case, the anion and cation were both found to be enriched in the collapsed state of the polymer. The collapsed state was favored because the Gnd+ cross-linked the polymer chains together. Moreover, the anion helped partition Gnd+ to the polymer surface. At higher salt concentrations (>1.5 M), GndSCN switched to stabilizing the uncollapsed state because a sufficient amount of Gnd+ and SCN- partitioned to the polymer surface to prevent cross-linking from occurring. Finally, in a third case, it was found that salts which interacted in an intermediate fashion with the polymer (e.g., GndCl) favored the uncollapsed conformation at all salt concentrations. These results provide a detailed, molecular-level, mechanistic picture of how Gnd+ influences the stability of polypeptides in three distinct physical regimes by varying the anion. It also helps explain the circumstances under which guanidinium salts can act as powerful and versatile protein denaturants.

• MeSH
• guanidin chemie   MeSH
• hydrofobní a hydrofilní interakce MeSH
• kationty MeSH
• peptidy chemie   MeSH
• spektroskopie infračervená s Fourierovou transformací MeSH
• termodynamika MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Názvy látek
• guanidin MeSH
• kationty MeSH
• peptidy MeSH