Computational Electrophysiology from a Single Molecular Dynamics Simulation and the Electrodiffusion Model

Area of Science:

• Biophysics
• Computational Biology
• Molecular Modeling

Background:

• High-resolution ion channel structures enable accurate electrophysiological property calculations.
• Traditional methods for computing voltage-dependent properties are often computationally intensive.

Purpose of the Study:

• To develop novel theoretical approaches for calculating electrophysiological properties from molecular dynamics (MD) simulations.
• To enable accurate predictions of ionic currents and selectivities using MD data at a single voltage or even without applied voltage.

Main Methods:

• Two theoretical approaches combining MD simulations with the electrodiffusion model were developed.
• Method 1: Utilizes one-sided ion fluxes and density profiles to determine the free energy profile.
• Method 2: Employs committor probabilities for ion transport at a selected voltage to determine the free energy profile.

Main Results:

• Both developed approaches accurately reproduced the potential of mean force (PMF) compared to traditional methods.
• The calculated current-voltage (I-V) dependence showed excellent agreement with computationally demanding methods.
• The computationally challenging reversal potential was readily calculated.
• Key electrodiffusion model assumptions were validated, and membrane voltage distribution was characterized.

Conclusions:

• The developed methods provide a computationally efficient and accurate way to determine ion channel electrophysiological properties from MD simulations.
• These approaches simplify the calculation of crucial properties like I-V dependence and reversal potential.
• The findings validate the electrodiffusion model and offer insights into voltage drop across membranes.