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

Abstract

In molecular simulations of macromolecules, it is desirable to limit the amount of solvent in the system to avoid spending computational resources on uninteresting solvent−solvent interactions. As a consequence, periodic boundary conditions are commonly used, with a simulation box chosen as small as possible, for a given minimal distance between images. Here, we describe how such a simulation cell can be set up for ensembles, taking into account a priori available or estimable information regarding conformational flexibility. Doing so ensures that any conformation present in the input ensemble will satisfy the distance criterion during the simulation. This helps avoid periodicity artifacts due to conformational changes. The method introduces three new approaches in computational geometry: (1) The first is the derivation of an optimal packing of ensembles, for which the mathematical framework is described. (2) A new method for approximating the α-hull and the contact body for single bodies and ensembles is presented, which is orders of magnitude faster than existing routines, allowing the calculation of packings of large ensembles and/or large bodies. 3. A routine is described for searching a combination of three vectors on a discretized contact body forming a reduced base for a lattice with minimal cell volume. The new algorithms reduce the time required to calculate packings of single bodies from minutes or hours to seconds. The use and efficacy of the method is demonstrated for ensembles obtained from NMR, MD simulations, and elastic network modeling. An implementation of the method has been made available online at http://haddock.chem.uu.nl/services/SQUEEZE/ and has been made available as an option for running simulations through the weNMR GRID MD server at

Glycan–protein interactions determine kinetics of N-glycan remodeling

“…The glycoprotein conformations were pre-sampled with the REST2 Hamiltonian replica exchange method 42 as described in Section S1 (ESI † ). The enhanced conformational sampling was used in two ways: (i) construction of an optimized triclinic box size 43 as described in Section S2 (ESI † ), and (ii) the sampled structures served as starting configurations for the production runs.…”

Molecular architecture and dynamics of SARS-CoV-2 envelope by integrative modeling

“…Furthermore, to measure convergence of the simulation we estimated the auto-correlation and equilibration times as previously proposed (69, 70). In addition to this global shape analysis, a more detailed view was determined following a procedure derived from an efficient method to calculate a molecular hull and contact bodies (71), which consists of extruding a triangulated sphere to the surface of the virion, in this case marked by the M dimer TMD centers of mass. The extrusion of each vertex was determined from a kernel density estimate of TMD distances from the center in angular coordinates.…”

N-Glycosylation Enhances Conformational Flexibility of Protein Disulfide Isomerase Revealed by Microsecond Molecular Dynamics and Markov State Modeling

“…Before the production runs, the box size was optimized using the near-densest-lattice-packing (NDLP) algorithm. , This procedure utilized REST2 sampled structures to render a triclinic box such that the minimum distance between the periodic images of every input conformation is at least 3.5 nm. The conservative choice of 3.5 nm minimum distance is far beyond twice the cutoff range and presents a compromise between the number of solvent molecules and hydrodynamic interactions across periodic boundaries.…”