Molecular dynamics simulation of ssDNA and cationic polythiophene

Abstract

In this work, molecular dynamics simulations of complexes composed of single strand DNA (ssDNA) sequences and cationic oligothiophenes are performed to understand experimental findings and the sensing ability of polythiophene electrolytes toward ssDNA. The simulation results exhibit no significant structural effect for replacing the cationic amine moiety with imidazole derivative on the side group of the oligomer. Adding a homopurine strand elongates the oligomer backbone; on the contrary, mixing up the homopyrimidine strand causes compression. On the other hand, these ssDNAs do not notably affect the compactness of the oligomer backbones. The anion-cation interactions play an essential role in the structural and spectroscopic change of cationic polythiophenes (CPTs) upon complexation with ssDNAs. The red shift of CPTs in the UV-VIS spectra with the addition of homopurine strands might be explained by the strong anion-cation, weak π -cation interactions, and high binding affinities. Nonpolar interactions (vdW and SA) and complex solvation energies dominate binding free energies. Hydrogen interaction analyses show that oligomers most likely approach the ssDNAs from their backbone upon complexation except for the duplex containing homopyrimidine strand and oligothiophene possessing imidazole derivative side chain.

Keywords: Cationic polythiophene; MM-PBSA; Molecular dynamics simulation; NAMD; Polyelectrolyte; ssDNA.