Abstract
Background/Objectives: Understanding the stabilization mechanisms of amino acid conformations in different solvent environments is crucial for elucidating biomolecular interactions and crystallization processes. This study presents a comprehensive computational investigation of glycine, the simplest amino acid, in both its canonical and zwitterionic forms when interacting with dimethyl sulfoxide (DMSO) molecules. Methods: Using density functional theory (DFT) calculations at the B3LYP/6-311++G(d,p) level with empirical dispersion corrections, we examined the conformational landscape of glycine-DMSO clusters with one and two DMSO molecules, as well as implicit solvent calculations, and compared them with analogous water clusters. Results: Our results demonstrate that while a single water molecule is insufficient to stabilize the zwitterionic form of glycine, one DMSO molecule successfully stabilizes this form through specific interactions between the S=O and the methyl groups of DMSO and the NH(3)(+) and the oxoanion group of zwitterionic glycine, respectively. Topological analysis of the electron density using QTAIM and NCI methods reveals the nature of these interactions. When comparing the relative stability between canonical and zwitterionic forms, we found that two DMSO molecules significantly reduce the energy gap to approximately 12 kJ mol(-1), suggesting that increasing DMSO coordination could potentially invert this stability. Implicit solvent calculations indicate that in pure DMSO medium, the zwitterionic form becomes more stable below 150 K, while remaining less stable at room temperature, contrasting with aqueous environments where the zwitterionic form predominates. Conclusions: These findings provide valuable insights into DMSO's unique role in biomolecular stabilization and have implications for protein crystallization protocols where DMSO is commonly used as a co-solvent.