In silico self-assembly and complexation dynamics of cationic lipids with DNA nanocages to enhance lipofection

DNA nanostructures are promising materials for drug delivery due to their unique topology, shape, size control, biocompatibility, structural stability and blood-brain-barrier (BBB) penetration capability. However, their cellular permeability is hindered by strong electrostatic repulsion from negatively charged cellular membranes, posing a significant obstacle to the use of DNA nanostructures as a drug delivery vehicle. Recent experimental studies have shown enhanced cellular uptake for the conjugate binary mixtures of DNA Tetrahedron (TDN) with cationic lipid N-[1-(2,3-dioleyloxy)propyl]-N,N,N-trimethylammonium chloride (DOTMA) compared to TDN alone. However, the cationic DOTMA lipids binding mechanism with the TDN nucleotides is still elusive. Using fully atomistic MD simulations, we aim to understand the molecular interactions that drive the formation and stability of the TDN-DOTMA binary complexes in a physiological environment. Our results uncovered that lipid concentration plays a crucial role in the energetics of the TDN-DOTMA association. We also report that distinct timescales are associated with the self-assembly of cationic DOTMA lipids first, followed by the complexation of self-assembled DOTMA lipid clusters with the TDN nucleotides, where electrostatics, hydrophobicity, and hydrogen bonding are the key interactions that drive the formation and stability of these complexes. These results provide molecular insights into TDN-DOTMA interactions, highlighting the lipid self-assembly dynamics, complex stability, and morphology, paving the way for the better rational design of cationic lipid functionalized DNA nanostructures for efficient drug delivery and transfection.