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. Here we describe the potential of estimating melting temperature (T(m)) for nonstandard oligonucleotides by using the correlation of the experimental T(m) with the calculated duplex binding energy (BE) for oligonucleotides of a given length. This method has been automated into a standardized molecular dynamics (MD) protocol through Pipeline Pilot (PP) using the CHARMm component in Discovery Studio (DS). Results will be presented showing the correlation of the predicted data with experiment for both standard and chemically modified siRNA and AON.