Abstract
The stereoselective stability of oligodeoxyribonucleotide methylphosphonates is examined using free-energy computer simulations. These modified phosphate linkages have the potential to be important antisense therapeutics that can be targeted on specific sequences of single- and double-stranded DNA, as well as crucial RNA messages. The stability of hybrid duplexes that contain these modified linkages is known experimentally to depend on the configuration of the chiral phosphonate center. Free-energy decomposition calculations were performed on three DNA oligomers to determine the origin of the structural interactions and physical properties that influence the relative stability of R and S methylphosphonate diastereomers. The strategy applied used free-energy decomposition methods to evaluate the free-energy contributions from selected groups. The results indicated that only three groups have a steric effect on the stability: the C2' and C3' substituents on the S diastereomer (5' side) and the C5' substituents on the R diastereomer (3' side). The balance considerably favors the R configuration in all the isomers studied and is not sequence dependent. The electrostatic effects were much more variable and were shown to be dependent on the conformation of duplex. The solvent interactions, however, were consistent and contributed favorably to the stability of the R over the S diastereomer. This favorable solvation energy for the R diastereomer was surprising (since the methyl group is more solvent exposed in this configuration) and was further supported by ab initio and associated free-energy calculations. This study concludes that oligonucleotides containing R-methylphosphonate linkages will normally form more stable duplexes than the corresponding S diastereomer irrespective of sequence, but also points that conformational changes may allow for a reversal in stability.
Keywords
Affiliated Institutions
Related Publications
Finite representation of an infinite bulk system: Solvent boundary potential for computer simulations
An approach is developed to obtain statistical properties similar to those of an infinite bulk system from computer simulations of a finite cluster. A rigorous theoretical formu...
A fast adaptive multigrid boundary element method for macromolecular electrostatic computations in a solvent
A fast multigrid boundary element (MBE) method for solving the Poisson equation for macromolecular electrostatic calculations in a solvent is developed. To convert the integral ...
gmx_MMPBSA: A New Tool to Perform End-State Free Energy Calculations with GROMACS
Molecular mechanics/Poisson-Boltzmann (Generalized-Born) surface area is one of the most popular methods to estimate binding free energies. This method has been proven to balanc...
Universal solvation model based on conductor-like screening model
Atomic surface tensions are parameterized for use with solvation models in which the electrostatic part of the calculation is based on the conductor-like screening model (COSMO)...
Ewald artifacts in computer simulations of ionic solvation and ion–ion interaction: A continuum electrostatics study
The use of Ewald and related methods to handle electrostatic interactions in explicit-solvent simulations of solutions imposes an artificial periodicity on systems which are inh...
Publication Info
- Year
- 1991
- Type
- article
- Volume
- 1
- Issue
- 3
- Pages
- 243-254
- Citations
- 10
- Access
- Closed
External Links
Social Impact
Social media, news, blog, policy document mentions
Citation Metrics
Cite This
Identifiers
- DOI
- 10.1089/ard.1991.1.243