Model of nanoscale deformation mechanisms of cellular structures could render different results depending on the molecular dynamics (MD) simulator chosen. Also, the comparison of different MD simulators is typically an intricate task, requiring all configurations be converted appropriately with available parameter choices. This study aims to perform and compare MD simulations between two MD programs (GROMACS and LAMMPS), in which a phospholipid bilayer is deformed under different strain states. The two systems produced similar deformation behaviors and strain state effect on bilayer failure. However, GROMACS produced more pores at lower strains, lower stress, and higher damage values. Multiple setting options and algorithm variations have been considered as possible explanations for the differences. Overall, the study aids in the cross-check of parameter settings and simulation results in MD research, particularly on the mechanical damage of bilayer membranes. Besides, based on that, GROMACS and LAMMPS could be further exploited with better reproducibility.

No description available.

No description available.

No description available.

No description available.

Simulation of processes in cellular membranes Abstract Many important processes in cells involve ions, e.g., fusion of synaptic vesi- cles with neuronal cell membranes is controlled by a divalent cation Ca2+ ; and the exchange of Na+ and K+ drives the the fast electrical signal transmis- sion in neurons. We have investigated model phospholipid membranes and their interactions with these biologically relevant ions. Using state-of-the-art molecular dynamics simulations, we accurately quantified their respective affinites towards neutral and negatively charged phospholipid bilayers. In order to achieve that, we developed a new model of phospholipids termed ECC-lipids, which accounts for the electronic polarization via the electronic continuum correction implemented as charge rescaling. Our simulations with this new force field reach for the first time a quantitative agreement with the experimental lipid electrometer concept for POPC as well as for POPS with all the studied cations. We have also examined the effects of transmembrane voltage on phospholipid bilayers. The electric field induced by the voltage exists exclusively in the hydrophobic region of the membrane, where it has an almost constant strength. This field affects the structure of nearby water molecules highlighting its importance in electroporation. 1