Molecular Dynamics Simulations Provide Insight into the Loading Efficiency of Proresolving Lipid Mediators Resolvin D1 and D2 in Cell Membrane-Derived Nanovesicles

Resolvins D1 and D2 (RvDs) are structural isomers and metabolites of docosahexaenoic acid, an omega-3 fatty acid, enzymatically produced in our body in response to acute inflammation or microbial invasion. Resolvins have been shown to play an essential role in the resolution of inflammation, tissue repair, and return to homeostasis and thus are actively pursued as potential therapeutics in treating inflammatory disorders and infectious diseases. However, effective in vivo delivery of RvDs continues to be a challenging task. Recent studies demonstrated that RvD1 or RvD2 loaded in cell membrane-derived nanovesides significantly increased therapeutic efficacy in treating murine peritonitis and ischemic stroke, respectively. The mechanistic details of how the subtle structural difference between RvD1 and RvD2 alters their molecular interactions with the membrane lipids of the nanovesides and thus affects the loading efficiency remain unknown. Here, we report the encapsulation profiles of the neutral and ionized species of both RvD1 and RvD2 determined with the cell membrane-derived nanovesides at pH values 5.4 and 7.4, respectively. Also, we performed microsecond time-scale all-atom molecular dynamics (MD) simulations in explicit water to elucidate the molecular interactions of both neutral and ionized species of RvD1 and RvD2 with the lipid bilayer using a model membrane system, containing 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and cholesterol. We found that the differences in the position and chirality of hydroxyl groups in RvD1 and RvD2 affected their location, orientation, and conformations within the bilayer. Surprisingly, the deprotonation of their carboxyl group caused their orientation and conformation to change from a fully extended one that is oriented in parallel to the membrane plane to a J-shaped bent conformation that is oriented perpendicular to the bilayer plane. Our studies offer valuable insight into the molecular interactions of RvD1/D2 with the lipid bilayer in atomistic details and provide a mechanistic explanation for the observed differences in the encapsulation profiles of RvD1 and RvD2, which may facilitate the rational design of nanovesicle-based therapeutics for treating inflammatory diseases.