Simulating the aggregation of DNA oligonucleotides

In this work we have studied, by means of Molecular Dynamics simulations, the process of denaturation and self-assembly of short oligonucleotides. Supramolecular ordering of DNA short strands is a promising field which is constantly enriched with new findings. Examples are provided by micellar and fibrils formations and due to the selectivity of DNA bindings, "intelligent" devices have been developed to perform simple logic operations. It is worth to notice that computer simulations of these DNA nanosystems would complement experiments with detailed insight into processes involved in self-assembly. In order to obtain an accurate description of the interactions involved in the complex structure of DNA we used oxDNA, a coarse-grained model developed by Ouldridge. We simulated the melting transition of 4, 6, and 8 base pair sequences. Sequence and length dependence were analyzed, specifically we compared thermodynamic parameters DeltaH, DeltaS and the melting temperature with literature results. Moreover, we have attempted to reproduce liquid crystal ordering of the ultrashort sequence GCCG at relatively high saline concentration, until now only experimentally observed in Bellini's works. We found that our simple model successfully reproduces the experimental phase sequence (isotropic, nematic, columnar) at T= 5 °C as a function of oligonucleotide concentration, and we fully characterized the microscopic structure of the three phases.