Molecular Dynamics Investigation of Clustering in Aqueous Glycine Solutions

Uncovering the Pathway of Serine Octamer Magic Number Cluster Formation during Electrospray Ionization: Experiments and Simulations

Electrospray ionization (ESI) of serine (Ser) solution generates Ser8H+ as an abundant magic number cluster. ESI clustering of most other solutes yields nonspecific stoichiometries. It is unclear why Ser8H+ dominates in the case of Ser, and how Ser8H+ forms during ESI. Even the location of Ser8H+ formation is contentious (in solution, in ESI droplets, or elsewhere). Here we unravel key aspects of the l-Ser8H+ formation pathway. Harsh ion sampling conditions promote the collision-induced dissociation (CID) of regular ESI analytes. Unexpectedly, Ser8H+ was seemingly resistant against CID during ion sampling, despite its extremely low tandem mass spectrometry (MS/MS) stability. This unusual behavior reveals that Ser8H+ forms during ion sampling. We propose the following pathway: (1) Nonspecific Ser clusters are released when ESI droplets evaporate to dryness. These initial clusters cover a wide size range, from a few Ser to hundreds or thousands of monomers. (2) The clusters undergo dissociation during ion sampling, mostly via successive loss of neutral monomers. For any source activation voltage, there is a subpopulation of clusters for which this CID cascade tends to terminate at the octamer level, culminating in Ser8H+-dominated product distributions. Mobile proton molecular dynamics simulations were used to model the entire pathway. Ser8H+ structures formed in these simulations were consistent with ion mobility experiments. The most compact structures resembled the model of [Scutelnic, V. J. Am. Chem. Soc. 2018, 140, 7554-7560], with numerous intermolecular salt bridges and H-bonds. Our findings illustrate how the interplay of association and dissociation reactions across phase boundaries can culminate in magic number clusters.

Orientational Behavior and Vibrational Response of Glycine at Aqueous Interfaces

Aqueous glycine plays many different roles in living systems, from being a building block for proteins to being a neurotransmitter. To better understand its fundamental behavior, we study glycine's orientational behavior near model aqueous interfaces, in the absence and presence of electric fields and biorelevant ions. To this purpose, we use a surface-specific technique called heterodyne-detected vibrational sum-frequency generation spectroscopy (HD-VSFG). Using HD-VSFG, we directly probe the symmetric and antisymmetric stretching vibrations of the carboxylate group of zwitterionic glycine. From their relative amplitudes, we infer the zwitterion's orientation near surfactant-covered interfaces and find that it is governed by both electrostatic and surfactant-specific interactions. By introducing additional ions, we observe that the net orientation is altered by the enhanced ionic strength, indicating a change in the balance of the electrostatic and surfactant-specific interactions.

Structure, energetics and dynamics in crowded amino acid solutions: a molecular dynamics study

A comprehensive understanding of crowding effects on biomolecular processes necessitates investigating the bulk thermodynamic and kinetic properties of the solutions with an accurate molecular representation of the crowded milieu. Recent studies have reparameterized the non-bonded dispersion interaction of solutes to precisely model intermolecular interactions, which would circumvent artificial aggregation as shown by the original force-fields. However, the performance of this reparameterization is yet to be assessed for concentrated crowded solutions in terms of investigating the hydration shell structure, energetics and dynamics. In this study, we perform molecular dynamics simulations of crowded aqueous solutions of five zwitterionic neutral amino acids (Gly, Ala, Thr, Pro, and Ser), mimicking the molecular crowding environment, using a modified AMBER ff99SB-ILDN force-field. We systematically examine and show that the reproducibility of the osmotic coefficients, density, viscosity and self-diffusivity of amino acids improves using the modified force-field in crowded concentrations. The modified force-field also improves the structuring of the solute solvation shells, solute interaction energy and convergence of tails of radial distribution functions, indicating reduction in the artificial aggregation. Our results also indicate that the hydrogen bonding network of water weakens and water molecules anomalously diffuse at small time scales in the crowded solutions. These results underscore the significance of examining the solution properties and anomalous hydration behaviour of water in crowded solutions, which have implications in shaping the structure and dynamics of biomolecules. The findings also illustrate the improvement in predicting bulk solution properties using the modified force-field, thereby providing an approach towards accurate modeling of crowded molecular solutions.

Molecular Dynamics Investigation of Giant Clustering in Small-Molecule Solutions: The Case of Aqueous PEHA

The importance of the formation of giant clusters in solution, in nature and industry, is increasingly recognized. However, relatively little attention has been paid to the formation of giant clusters in solutions of small, relatively soluble but nonamphiphilic molecules. In this work, we present a general methodology based on molecular dynamics that can be used to investigate such systems. As a case study, we focus on the formation of apparently stable clusters of pentaethylenehexamine (PEHA) in water. These clusters have been used as templates for the construction of bioinspired silica nanoparticles. To better understand clustering in this system, we study the effect of PEHA protonation state (neutral, +1, and +2) and counterion type (chloride or acetate) on PEHA clustering in dilute aqueous solutions (200 and 400 mM) using large-scale classical molecular dynamics. We find that large stable clusters are formed by singly charged PEHA with chloride or acetate as the counterion, although it is not clear for the case with acetate whether bulk phase separation, that might lead to precipitation, would eventually occur. Large clusters also appear to be stable for doubly charged PEHA with acetate, the less soluble counterion. We attribute this behavior to a form of complex coacervation, observed here for relatively small and highly soluble molecules (PEHA + counterion) rather than the large polyions usually found to form such coacervates. We discuss whether this behavior might also be described by an effective SALR (short-range attraction, long-range repulsion) interaction. This work might help future studies of additives for the design of novel bioinspired templated nanomaterials and of giant clustering in small-molecule solutions more generally.