Absolute pKa Values and Solvation Structure of Amino Acids from Density Functional Based Molecular Dynamics...

Area of Science:

• Computational Chemistry
• Biophysical Chemistry
• Molecular Dynamics

Background:

• Accurate pKa values are crucial for understanding protein function and biochemical reactions.
• Previous molecular dynamics methods faced challenges in precisely calculating pKa values for amino acid residues.

Purpose of the Study:

• To calculate absolute pKa values for key amino acid side chains and termini.
• To refine molecular dynamics simulations for improved pKa prediction accuracy.

Main Methods:

• Utilized a revised density functional-based molecular dynamics simulation technique.
• Treated acid deprotonation as a dissociation reaction, incorporating a hydronium ion correction.
• Calculated acidity constants from vertical energy gaps and solvation free energies.

Main Results:

• Achieved an unsigned mean error of 2.1 pKa units relative to experimental data.
• Reported a maximum error of 4.0 pKa units and a mean statistical uncertainty of ±1.1 pKa units.
• Analyzed solvation structures of protonated and deprotonated amino acids.

Conclusions:

• The revised simulation technique provides accurate absolute pKa calculations for amino acids.
• The findings offer valuable reference data for developing future molecular force fields.
• This method enhances the understanding of protonation states in biological systems.