The stability of Seeman JX DNA topoisomers of paranemic crossover (PX) molecules as a function of crossover number

We use molecular dynamics simulations in explicit water and salt (Na+) to determine the effect of varying the number of crossover points on the structure and stability of the PX65 paranemic crossover DNA molecule and its JXM topoisomers (M denotes the number of missing crossover points), recently synthesized by the Seeman group at New York University. We find that PX65, with six crossover points, is the most stable, and that the stability decreases monotonically with the number of crossover points PX65 > JX1 > JX2 > JX3 > JX4, with 6, 5, 4, 3 and 2 crossover points, respectively. Thus, for PX65/JX1, the strain energy is similar to3 kcal/mol/bp, while it is similar to13 kcal/mol/bp for JX2, JX3 and JX4. Another measure of the stability is the change in the structure from the minimum energy structure to the equilibrium structure at 300 K, denoted as root-mean-square deviation in coordinates (CRMSD). We find that CRMSD is similar to 3.5 Angstrom for PX65, increases to 6 Angstrom for JX1 and increases to 10 Angstrom for JX2/JX3/JX4. As the number of crossover points decreases, the distance between the two double helical domains of the PX/JX molecules increases from similar to 20 Angstrom for PX65 to 23 Angstrom for JX4. This indicates that JX2, JX3 and JX4 are less likely to form, at least in with Na+. However, in all the cases, the two double helical domains have average helicoidal parameters similar to a typical B-DNA of similar length and base sequence.