Molecular dynamics simulation and binding energy calculation for estimation of oligonucleotide duplex thermostability in RNA-based therapeutics.

For oligonucleotide-based therapeutics, a thorough understanding of the thermodynamic properties of duplex formation is critical to developing stable and potent drugs. For unmodified small interfering RNA (siRNA), DNA antisense oligonucleotide (AON) and locked nucleic acid (LNA), DNA/LNA modified oligonucleotides, nearest neighbor (NN) methods can be effectively used to quickly and accurately predict duplex thermodynamic properties such as melting point. Unfortunately, for chemically modified olignonucleotides, there has been no accurate prediction method available.