Structural, dynamic, and electrostatic properties of fully hydrated DMPC bilayers from molecular dynamics simulations...

Area of Science:

• Biophysics
• Computational Chemistry
• Materials Science

Background:

• Lipid bilayers are crucial biological environments.
• Accurate molecular dynamics simulations are essential for understanding membrane properties.
• Advancements in computational power are needed for detailed atomic-resolution modeling.

Purpose of the Study:

• To present molecular dynamics simulations of hydrated DMPC bilayers using GPUs.
• To validate the FEN ZI code and particle mesh Ewald method (PME) implementation for electrostatic interactions.
• To assess the impact of system size on simulated bilayer properties.

Main Methods:

• Utilized state-of-the-art non-polarizable force fields and a GPU-enabled molecular dynamics code (FEN ZI).
• Employed the particle mesh Ewald method (PME) for exact treatment of electrostatic interactions.
• Performed simulations on DMPC bilayer systems with varying atom counts (17,004 to 273,936 atoms) in explicit solvent.

Main Results:

• Simulated bilayer structural properties (densities, order parameters) and electrostatic properties (dipole potential) showed excellent agreement with experimental data and previous simulations.
• No statistically significant differences in structural or electrostatic properties were observed between small ( < 100 lipids) and large ( > 1000 lipids) bilayer systems.
• The study highlights the suitability of GPU technology for large-scale molecular dynamics of lipid bilayers.

Conclusions:

• Graphics processing unit (GPU) technology and the FEN ZI code provide a powerful platform for accurate molecular dynamics simulations of lipid bilayers.
• Small-scale simulations can effectively capture key structural and electrostatic properties, offering computational efficiency.
• This approach enables enhanced exploration of biological membrane systems and integral membrane proteins at atomic resolution.