Atomic-Scale Molecular Dynamics Simulations of DNA-Polycation Complexes: Two Distinct Binding Patterns

Synthetic cationic polymers represent a promising class of delivery vectors for gene therapy. Here, we employ atomistic molecular dynamics simulations to gain insight into the structure and properties of complexes of DNA with four linear polycations: polyethylenimine (PEI), poly-L-lysine (PLL), polyvinylamine (PVA), and polyallylamine (PAA). These polycations differ in their polymer geometries, protonation states, and hydrophobicities of their backbone chains. Overall, our results demonstrate for the first time the existence of two distinct patterns of binding of DNA with polycations. For PEI, PLL, and PAA, the complex is stabilized by the electrostatic attraction between protonated amine groups of the polycation and phosphate groups of DNA. In contrast, PVA demonstrates an alternative binding pattern as it gets embedded into the DNA major groove. It is likely that both the polymer topology and affinity of the backbone chain of PVA to the DNA groove are responsible for such behavior. The differences in binding patterns can have important biomedical implications: embedding PVA into a DNA groove makes it less sensitive to changes in the aqueous environment (pH level, ionic strength, etc.) and could therefore hinder the intracellular release of genetic material from a delivery vector, leading to lower transfection activity.