Molecular dynamics simulations of membrane proteins under asymmetric ionic concentrations

Area of Science:

• Computational biophysics
• Membrane biophysics
• Ion channel electrophysiology

Background:

• Simulating biomembrane systems with realistic ionic gradients is computationally challenging.
• Maintaining periodic boundary conditions while imposing asymmetric salt concentrations requires novel approaches.
• Accurate modeling of ion transport is crucial for understanding cellular function.

Purpose of the Study:

• To develop a computational method for molecular dynamics simulations of biomembranes under realistic ionic gradients.
• To validate the method's accuracy by comparing simulation results with experimental data.
• To enable the simulation of complex ion channel systems under physiological conditions.

Main Methods:

• Introduction of a nonperiodic energy step for ionic species at the simulation cell edge.
• Application to membrane slab and phospholipid bilayer models.
• Simulation of current-voltage (I-V) curves for bacterial porin OmpF and Kv1.2 channel.

Main Results:

• The method successfully simulates biomembranes with asymmetric salt concentrations using periodic boundary conditions.
• Calculated reversal potential for OmpF (28.6 mV) closely matches experimental values (26-27 mV).
• First-time simulation of Kv1.2 channel under physiological K+ and Na+ gradients and electrostatic potential.

Conclusions:

• The developed nonperiodic energy-step method allows for accurate, quantitative simulations of nonequilibrium membrane transport.
• This computational approach provides a powerful tool for studying ion channel function and biomembrane dynamics.
• The method opens new possibilities for simulating complex biological systems under physiologically relevant conditions.