Explicit solvent molecular dynamics simulation of duplex formed by the modified oligonucleotide with alternating phosphate/phosphonate internucleoside linkages and its natural counterpart.

Impact of the internucleoside linkage modification by inserting a methylene group on the ability of the modified oligonucleotide to hybridize with a natural DNA strand was studied by fully solvated molecular dynamics (MD) simulations. Three undecamer complexes were analyzed: natural dT(11).dA(11) duplex as a reference and two its analogs with alternating modified and natural linkages in the deoxyadenosine chain. The isopolar, non-isosteric modified linkages were of 5'-O-PO(2)-CH(2)-O-3' (5'PC3') or 5'-O-CH(2)-PO(2)-O-3' (5'CP3') type. Simulations were performed by using the AMBER 5.0 software package with the force field completed by a set of parameters needed to model the modified segments. Both modifications were found to lead to double helical complexes, in which the thymidine strand as well as deoxyriboses and unmodified linkages in the adenosine strand adopted conformations typical for the B-type structure. For each of the two conformational richer modified linkages two stable conformations were found at 300 K: the -ggt and ggt for the 5'PC3' and ggg, tgg for the 5'CP3', respectively. Both modified chains adopted helical conformations with heightened values of the inclination parameter but without affecting the Watson-Crick hydrogen bonds.