Diffusion constant of K+ inside Gramicidin A: a comparative study of four computational methods

Area of Science:

• Biophysics
• Computational Chemistry
• Ion Channel Function

Background:

• Gramicidin A (GA) is a channel protein crucial for ion transport across cell membranes.
• Understanding ion diffusion within GA is vital for cellular electrophysiology and drug development.
• Accurate calculation of local diffusion constants is essential for modeling ion transport.

Purpose of the Study:

• To determine the local diffusion constant of potassium ions (K+) within the Gramicidin A channel.
• To compare the reliability of four distinct molecular dynamics (MD) simulation methods for calculating ion diffusion.
• To identify the most accurate computational approaches for studying ion transport in confined environments.

Main Methods:

• Four computational methods based on molecular dynamics (MD) simulations were employed: Mean Square Displacement (MSD), Velocity Autocorrelation Function (VACF), Second Fluctuation Dissipation Theorem (SFDT), and Generalized Langevin Equation for a Harmonic Oscillator (GLE-HO).
• Methods were validated using K+ diffusion in bulk water.
• Systematic forces from the membrane-channel environment were analyzed for their influence on diffusion calculations.

Main Results:

• MSD and VACF methods were found to be unreliable for K+ diffusion within GA due to bias from systematic forces.
• SFDT and GLE-HO techniques successfully unbias the influence of systematic forces.
• SFDT and GLE-HO predicted a K+ diffusion constant within GA approximately 10 times smaller than in bulk water.
• Reliable K+ diffusion constants within GA using SFDT and GLE-HO require extensive MD sampling (tens of nanoseconds).

Conclusions:

• SFDT and GLE-HO are reliable methods for calculating ion diffusion constants within biological channels like Gramicidin A.
• Ion diffusion within the Gramicidin A channel is significantly hindered compared to bulk water.
• Extensive molecular dynamics simulation sampling is critical for obtaining accurate diffusion properties in complex biological systems.