Comparative Analysis of Structural and Dynamical Properties of Lipid Membranes Simulated with the AMBER Lipid21 ForceField Using SPC/E, TIP3P, TIP3P-FB, TIP4P-FB, TIP4P-Ew, TIP4P/2005, TIP4P-D, and OPC Water Models

We have conducted all atom molecular dynamics simulations of POPC and DPPC lipid bilayers using AMBER Lipid21 force field with eight different water models, including SPC/E, TIP3P, TIP3P-FB, TIP4P-FB, TIP4P-Ew, TIP4P/2005, TIP4P-D, and OPC, to identify the most compatible one without any modification. A number of parameters have been computed in order to understand the structure of the lipid bilayer: Area per lipid, Isothermal compressibility modulus, average Volume per lipid, electron density profile, bilayer thickness, X-ray and neutron scattering form factors, deuterium order parameter, and radial distribution function. The estimated Area per lipid, Isothermal compressibility factor, volume per lipid and bilayer thickness are highly consistent with experimental results for the SPC/E water model, indicating its suitability with the AMBER Lipid21 force field, insted of any modification. The bilayer electron density profiles of both the lipid bilayers demonstrate a little augmentation of water penetration with respect to the membrane surface for TIP4P-D water model. However, the experimental X-ray and neutron scattering form factors are aligning well with the simulated results for all studied water models, and TIP4P-D shows better for X-ray data. The deuterium order parameter for lipid acyl chains value less than 0.25 for all observed water models, depicting their disorderness for both the lipid bilayers. The lateral diffusion and reorientation autocorrelation function of the lipid molecules in both the bilayers are computed to reveal their dynamics across all water models. In comparison to other water models, the simulated trajectories predict better structure and reasonably fair dynamic properties for the SPC/E water model. The TIP4P-Ew water model reproduces the lateral diffusion co-efficient in close agreement with experiment. Reorientational dynamics for both the lipids in the bilayers for eight different water models are observed; the presence of slow and slowest time components corresponds to the lipid axial motion (wobble motion) and Twist/Splay motions. So, in view of the overall performance of the different water models with the AMBER Lipid21 all atom force field in reproducing membrane physical properties, the SPC/E water model appears to be an optimal choice.