SQUEEZE-E: The Optimal Solution for Molecular Simulations with Periodic Boundary Conditions

Area of Science:

• Computational chemistry
• Structural biology
• Biophysics

Background:

• Molecular simulations require efficient solvent management to minimize computational costs.
• Periodic boundary conditions are standard, necessitating small simulation boxes and minimal image distances.
• Conformational flexibility can lead to artifacts if not accounted for in simulation cell setup.

Purpose of the Study:

• To develop a method for setting up simulation cells for macromolecular ensembles that accounts for conformational flexibility.
• To avoid periodicity artifacts in molecular dynamics simulations.
• To reduce computational resources spent on solvent-solvent interactions.

Main Methods:

• Derivation of an optimal packing framework for molecular ensembles.
• Development of a fast approximation method for α-hull and contact bodies.
• Implementation of a routine for finding reduced lattice bases for minimal cell volumes.

Main Results:

• New algorithms significantly reduce the time for calculating single-body packings from minutes/hours to seconds.
• The method ensures conformational ensembles fit within simulation cells, preventing artifacts.
• Demonstrated efficacy using ensembles from NMR, MD simulations, and elastic network modeling.

Conclusions:

• The developed computational geometry methods enable efficient and artifact-free molecular simulations of flexible macromolecules.
• The implementation is available online and as an option for GRID MD servers.
• This approach optimizes resource utilization in computational biophysics.