Molecular Dynamics Simulations of Liposomes: Structure, Dynamics, and Applications

Liposomes are nanoscale, spherical vesicles composed of phospholipid bilayers, typically ranging from 50 to 200 nm in diameter. Their unique ability to encapsulate both hydrophilic and hydrophobic molecules makes them powerful nanocarriers for drug delivery, diagnostics, and vaccine formulations. Several FDA-approved formulations such as Doxil^® (Baxter Healthcare Corporation, Deerfield, IL, USA), AmBisome^® (Gilead Sciences, Inc., Foster City, CA, USA), and Onivyde^® (Ipsen Biopharmaceuticals, Inc., Basking Ridge, NJ, USA) highlight their clinical significance. This review provides a comprehensive synthesis of how molecular dynamics (MD) simulations, particularly coarse-grained (CG) and atomistic approaches, advance our understanding of liposomal membranes. We explore key membrane biophysical properties, including area per lipid (APL), bilayer thickness, segmental order parameter (SCD), radial distribution functions (RDFs), bending modulus, and flip-flop dynamics, and examine how these are modulated by cholesterol concentration, PEGylation, and curvature. Special attention is given to curvature-induced effects in spherical vesicles, such as lipid asymmetry, interleaflet coupling, and stress gradients across the leaflets. We discuss recent developments in vesicle modeling using tools such as TS2CG, CHARMM-GUI Martini Maker, and Packmol, which have enabled the simulation of large-scale, compositionally heterogeneous systems. The review also highlights simulation-guided strategies for designing stealth liposomes, tuning membrane permeability, and enhancing structural stability under physiological conditions. A range of CG force fields, MARTINI, SPICA, SIRAH, ELBA, SDK, as well as emerging machine learning (ML)-based models, are critically assessed for their strengths and limitations. Despite the efficiency of CG models, challenges remain in capturing long-timescale events and atomistic-level interactions, driving the development of hybrid multiscale frameworks and AI-integrated techniques. By bridging experimental findings with in silico insights, MD simulations continue to play a pivotal role in the rational design of next-generation liposomal therapeutics.