The coarse-grained Martini force field is widely used in biomolecular simulations. Here, we present the refined model, Martini 3 (http://cgmartini.nl), with an improved interaction balance, new bead types, and expanded ability to include specific interactions representing, e.g. hydrogen bonding and electronic polarizability. The new model allows more accurate predictions of molecular packing and interactions in general, which is exemplified with a vast and diverse set of applications, ranging from oil/water partitioning and miscibility data to complex molecular systems, involving protein-protein and protein-lipid interactions and material science applications as ionic liquids and aedamers.

Triacylglycerols (TG) are synthesized at the endoplasmic reticulum (ER) bilayer and packaged into organelles called lipid droplets (LDs). LDs are covered by a single phospholipid monolayer contiguous with the ER bilayer. This connection is used by several monotopic integral membrane proteins, with hydrophobic membrane association domains (HDs), to diffuse between the organelles. However, how proteins partition between ER and LDs is not understood. Here, we employed synthetic model systems and found that HD-containing proteins strongly prefer monolayers and returning to the bilayer is unfavorable. This preference for monolayers is due to a higher affinity of HDs for TG over membrane phospholipids. Protein distribution is regulated by PC/PE ratio via alterations in monolayer packing and HD-TG interaction. Thus, HD-containing proteins appear to non-specifically accumulate to the LD surface. In cells, protein editing mechanisms at the ER membrane would be necessary to prevent unspecific relocation of HD-containing proteins to LDs.

No abstract