Protein molecular dynamics with electrostatic force entirely determined by a single Poisson-Boltzmann calculation

Area of Science:

• Computational chemistry
• Molecular dynamics simulations
• Biophysics

Background:

• Accurate simulation of electrostatic forces is crucial for understanding protein dynamics.
• Traditional methods like molecular dynamics (MD) can be computationally intensive.
• Integrating electrostatic calculations into stochastic dynamics offers a potential improvement.

Purpose of the Study:

• To introduce and validate a new finite difference stochastic dynamics (FDSD) procedure.
• To incorporate electrostatic interactions, including intramolecular Coulombic forces and solvation effects, into protein dynamics simulations.
• To determine the optimal relative dielectric constant for protein interior simulations.

Main Methods:

• The finite difference Poisson-Boltzmann (FDPB) method was used to calculate electrostatic forces.
• FDPB was integrated with the GROMOS96 force field to create the FDSD procedure.
• Simulations were performed on an insulin dimer, testing various dielectric constants (optimal found to be 17).
• Results were compared against explicit water MD and conventional stochastic dynamics (SD) simulations.

Main Results:

• The FDSD method, with an optimal protein dielectric constant of 17, successfully simulated an insulin dimer.
• The simulation data showed good agreement with established MD and SD methods.
• The study highlights the importance of selecting an appropriate dielectric constant for accurate results.

Conclusions:

• The FDSD method, coupled with the GROMOS96 force field and an optimal dielectric constant, is a suitable approach for studying protein dynamics and structure in solution.
• This method provides a computationally efficient alternative for molecular simulations.
• Further studies can explore its application to a wider range of protein systems.