Curvature Footprints of Transmembrane Proteins in Simulations with the Martini Force Field.

Membranes play essential roles in biological systems and are tremendously diverse in the topologies and chemical and elastic properties that define their functions. In many cases, a given membrane may display considerable heterogeneity, with localized clusters of lipids and proteins exhibiting distinct characteristics compared to adjoining regions. These lipid-protein assemblies can span nanometers to micrometers and are associated with cellular processes such as transport and signaling. While lipid-protein assemblages are dynamic, they can be stabilized by coupling between local membrane composition and shape. Due to the inherent difficulty in resolving atomistic details of membrane proteins in their native lipid environments, these complexes are notoriously challenging to study experimentally; however, molecular dynamics (MD) simulations might be a viable alternative. Here, we aim to assess the utility of coarse-grained (CG) MD simulations with the Martini force field for studying membrane curvature induced by transmembrane (TM) proteins that are reported to generate local curvature. The direction and magnitude of curvature induced by five different TM proteins, as well as certain lipid-protein and protein-protein interactions, were found to be in good agreement with available reference data.