Multiscale molecular dynamics (MD) simulations using Martini coarse-grained (CG) and all-atom (AA) forcer fields are commonly used in membrane protein studies. In particular, reverse-mapping an equilibrated CG model to an AA model offers an efficient way of preparing large membrane protein systems with complex protein shapes and lipid compositions. Here, we report that this hybrid CG-equilibrium-AA-production protocol may artificially increase lipid density and decrease hydration in ion channel pores walled with transmembrane gaps. To understand the origin of this conundrum, we conducted replicas of CG, AA, and CG reverse-mapped AA simulations of the pore domain of the mechanosensitive Piezo1 channel in a non-conducting conformation. Lipid/water density analysis and free energy calculations reveal that the lack of initial pore hydration allows adjacent lipids to enter the pore lumen through gaps between pore helices during CG simulation. Due to the mismatch between CG and AA lipid kinetics, these pore lipids remain trapped in the subsequent AA simulations, despite unfavorable binding free energy. We tested several CG equilibrium protocols and found that a protocol restraining the whole lipid produces pore hydration consistent with AA results, thus eliminating this artifact for further studies of lipid gating and protein-lipid interactions.

Importance of molecular dynamics equilibrium protocol on protein-lipid Interactions near channel pore

Where thousands of chemistry professionals meet to share ideas and advance scientific and technical knowledge. Spring 2023 abstracts, posters, and presentations will launch at the end of March.

ACS Spring 2023

The official archive and repository for all ACS Meetings & Expositions abstracts and select posters or presentations.