Ionizable lipids penetrate phospholipid bilayers with high phase transition temperatures: perspectives from free energy calculations

The efficacies of modern gene-therapies strongly depend on their contents. At the same time the most potent formulations might not contain the best compounds. In this work we investigated the effect of phospholipids and their saturation on the binding ability of (6Z,9Z,28Z,31Z)-heptatriacont-6,9,28,31-tetraene-19-yl 4-(dimethylamino) butanoate (DLin-MC3-DMA) to model membranes at the neutral pH. We discovered that DLin-MC3-DMA has affinity to the most saturated monocomponent lipid bilayer 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and an aversion to the unsaturated one 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). The preference to a certain membrane was also well-correlated to the phase transition temperatures of phospholipid bilayers, and to their structural and dynamical properties. Additionally, in the case of the presence of DLin-MC3-DMA in the membrane with DOPC the ionizable lipid penetrated it, which indicates possible synergistic effects. Comparisons with other ionizable lipids were performed using a model lipid bilayer of 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC). Particularly, the lipids heptadecan-9-yl 8-[2-hydroxyethyl-(6-oxo-6-undecoxyhexyl)amino]octanoate (SM-102) and [(4-hydroxybutyl) azanediyl] di(hexane-6,1-diyl) bis(2-hexyldecanoate) (ALC-0315) from modern mRNA-vaccines against COVID-19 were investigated and force fields parameters were derived for those new lipids. It was discovered that ALC-0315 binds strongest to the membrane, while DLin-MC3-DMA is not able to reside in the bilayer center. The ability to penetrate the membrane POPC by SM-102 and ALC-0315 can be related to their saturation, comparing to DLin-MC3-DMA.