Complexation of single stranded RNA with an ionizable lipid: an all-atom molecular dynamics simulation study

Abstract

Complexation of a lipid-based ionizable cationic molecule (referred to as DML: see main text) with RNA in an aqueous medium was examined in detail by means of fully atomistic molecular dynamics simulations. The different stages of the DML–RNA association process were explored, while the structural characteristics of the final complex were described. The self-assembly process of the DML molecules was examined in the absence and in the presence of nucleotide sequences of different lengths. The formed DML clusters were described in detail in terms of their size and composition and were found to share common features in all the examined systems. Different timescales related to their self-assembly and their association with RNA were identified. It was found that beyond a time period of a few tens of ns, a conformationally stable DML–RNA complex was formed, characterized by DML clusters covering the entire contour of RNA. In a system with a 642-nucleotide sequence, the average size of the complex in the longest dimension was found to be close to 40 nm. The DML clusters were characterized by a rather low surface charge, while a propensity for the formation of larger size clusters close to RNA was noted. Apart from hydrophobic and electrostatic interactions, hydrogen bonding was found to play a key-role in the DML–DML and in the DML–RNA association. The information obtained regarding the structural features of the final complex, the timescales and the driving forces associated with the complexation and the self-assembly processes provide new insight towards a rational design of optimized lipid-based ionizable cationic gene delivery vectors.