Protocol for MM/PBSA molecular dynamics simulations of proteins

Area of Science:

• Computational chemistry
• Biophysics
• Molecular modeling

Background:

• Continuum solvent models are crucial for studying biomolecular processes like protein folding.
• Integrating atomic detail molecular dynamics with continuum solvent models is an active research area.
• The Poisson-Boltzmann model is fundamental for biomolecular electrostatics but rarely used for solvation forces in dynamics.

Purpose of the Study:

• To develop a fast and reliable methodology for implementing continuum solvent forces within molecular dynamics.
• To enable the use of continuum solvent models for calculating solvation forces in molecular dynamics simulations.
• To bridge the gap between atomic detail simulations and continuum solvent approaches.

Main Methods:

• Implementation of continuum forces into standard molecular mechanics and dynamics algorithms.
• Development of a novel methodology for fast and reliable continuum force calculations.
• Conducting a 1-nanosecond unrestrained molecular dynamics simulation of a small protein.

Main Results:

• The proposed methodology efficiently incorporates continuum forces into molecular dynamics.
• Solvation forces were successfully calculated using the Poisson-Boltzmann solvent accessible surface area model.
• Simulation results quantitatively matched those obtained from explicit solvent molecular dynamics.

Conclusions:

• The developed method provides a fast and reliable way to implement continuum forces in molecular dynamics.
• This approach offers a viable alternative to computationally expensive explicit solvent simulations for certain applications.
• The findings support the use of continuum solvent models for accurate biomolecular simulations.