Optimizing Coarse-Grained Models for Large-Scale Membrane Protein Simulation

Coarse-grained (CG) models have been developed for studying membrane proteins at physiologically relevant scales. Such methods, including popular CG lipid models, exhibit stability and efficiency at moderate scales, but they can become impractical or even unusable beyond a critical size due to various technical issues. Here, we report that these scale-dependent issues can arise from progressively slower relaxation dynamics and become confounded by unforeseen instabilities observed only at larger scales. To address these issues, we systemically optimized a 4-site solvent-free CG lipid model that is suitable for conducting micron-scale molecular dynamics simulations of membrane proteins under various membrane properties. We applied this lipid model to explore the long-range membrane deformation induced by a large mechanosensitive ion channel, PIEZO. We show that the optimized CG models are powerful in elucidating the structural and dynamic interplay between PIEZO and the membrane. Furthermore, we anticipate that our methodological insights can prove useful for resolving issues stemming from scale-dependent limitations of similar CG methodologies.