A comparative theoretical study of dipeptide solvation in water

Area of Science:

• Computational Chemistry
• Biophysics
• Theoretical Chemistry

Background:

• Understanding solvation and electrostatic properties of peptides is crucial in biophysical studies.
• Molecular dynamics (MD) simulations are widely used to investigate these properties.

Purpose of the Study:

• To compare different theoretical methods for simulating the solvation and electrostatic properties of the zwitterionic dipeptide glycine-alanine in water.
• To assess the accuracy of classical force fields and hybrid quantum mechanical/molecular mechanical (QM/MM) methods against fully quantum mechanical (QM) calculations.

Main Methods:

• Molecular dynamics simulations were conducted on the zwitterionic form of glycine-alanine in water.
• Methods ranged from classical force fields to mixed QM/MM simulations and fully QM Car-Parrinello calculations.

Main Results:

• Solvation patterns were largely consistent across all employed methods for most atoms.
• Significant differences in solvation were observed for the carboxy terminus, amino terminus, and the backbone amide NH group.
• Hybrid QM/MM simulations effectively reproduced the solvation patterns obtained from fully QM simulations.

Conclusions:

• The choice of theoretical method can influence the depiction of solvation for specific charged groups in dipeptides.
• These findings have implications for theoretical studies involving peptides and proteins with charged side chains in aqueous environments.
• Hybrid QM/MM methods offer a reliable approach for capturing essential solvation characteristics comparable to more computationally intensive QM methods.