Controlled sulfur-based engineering confers mouldability to phosphorothioate antisense oligonucleotides

Abstract Phosphorothioates (PS) have proven their effectiveness in the area of therapeutic oligonucleotides with applications spanning from cancer treatment to neurodegenerative disorders. Initially, PS substitution was introduced for the antisense oligonucleotides (PS ASOs) because it confers an increased nuclease resistance meanwhile ameliorates cellular uptake and in-vivo bioavailability. Thus, PS oligonucleotides have been elevated to a fundamental asset in the realm of gene silencing therapeutic methodologies. But, despite their wide use, little is known on the possibly different structural changes PS-substitutions may provoke in DNA·RNA hybrids. Additionally, scarce information and significant controversy exists on the role of phosphorothioate chirality in modulating PS properties. Here, through comprehensive computational investigations and experimental measurements, we shed light on the impact of PS chirality in DNA-based antisense oligonucleotides; how the different phosphorothioate diastereomers impact DNA topology, stability and flexibility to ultimately disclose pro-Sp S and pro-Rp S roles at the catalytic core of DNA Exonuclease and Human Ribonuclease H; two major obstacles in ASOs-based therapies. Altogether, our results provide full-atom and mechanistic insights on the structural aberrations PS-substitutions provoke and explain the origin of nuclease resistance PS-linkages confer to DNA·RNA hybrids; crucial information to improve current ASOs-based therapies.

[1]  Serdal Kirmizialtin,et al.  Refining the RNA Force Field with Small-Angle X-ray Scattering of Helix–Junction–Helix RNA , 2022, The journal of physical chemistry letters.

[2]  R. Bridges,et al.  Open reading frame correction using splice-switching antisense oligonucleotides for the treatment of cystic fibrosis , 2022, Proceedings of the National Academy of Sciences.

[3]  E. Rosta,et al.  The Role of Conserved Residues in the DEDDh Motif: the Proton-Transfer Mechanism of HIV-1 RNase H , 2021, ACS Catalysis.

[4]  D. Gaudet,et al.  Efficacy and safety of volanesorsen in patients with multifactorial chylomicronaemia (COMPASS): a multicentre, double-blind, randomised, placebo-controlled, phase 3 trial. , 2021, The lancet. Diabetes & endocrinology.

[5]  Pengfei Li,et al.  Alchemical Binding Free Energy Calculations in AMBER20: Advances and Best Practices for Drug Discovery , 2020, J. Chem. Inf. Model..

[6]  M. Orozco,et al.  DFFR: A NEW METHOD FOR HIGH-THROUGHPUT RECALIBRATION OF AUTOMATIC FORCE-FIELDS FOR DRUGS. , 2020, Journal of chemical theory and computation.

[7]  Michael R. Shirts,et al.  Best Practices for Alchemical Free Energy Calculations , 2020, 2008.03067.

[8]  S. Crooke,et al.  The Interaction of Phosphorothioate Containing RNA Targeted Drugs with Proteins is a Critical Determinant of The Therapeutic Effects of These Agents. , 2020, Journal of the American Chemical Society.

[9]  A. Khvorova,et al.  The chemical structure and phosphorothioate content of hydrophobically modified siRNAs impact extrahepatic distribution and efficacy , 2020, Nucleic acids research.

[10]  S. Crooke,et al.  Gapmer Antisense Oligonucleotides Targeting 5S Ribosomal RNA Can Reduce Mature 5S Ribosomal RNA by Two Mechanisms. , 2020, Nucleic acid therapeutics.

[11]  M. Orozco,et al.  Surviving the deluge of biosimulation data , 2019, WIREs Computational Molecular Science.

[12]  E. Mercuri,et al.  Increased dystrophin production with golodirsen in patients with Duchenne muscular dystrophy , 2020, Neurology.

[13]  S. Crooke,et al.  Understanding the effect of controlling phosphorothioate chirality in the DNA gap on the potency and safety of gapmer antisense oligonucleotides , 2020, Nucleic acids research.

[14]  M. De Vivo,et al.  Recruiting Mechanism and Functional Role of a Third Metal Ion in the Enzymatic Activity of 5′ Structure-Specific Nucleases , 2020, Journal of the American Chemical Society.

[15]  Jinkuk Kim,et al.  Patient-Customized Oligonucleotide Therapy for a Rare Genetic Disease. , 2019, The New England journal of medicine.

[16]  T. Yokota,et al.  Viltolarsen for the treatment of Duchenne muscular dystrophy. , 2019, Drugs of today.

[17]  M. Marcia,et al.  A Transient and Flexible Cation-π Interaction Promotes Hydrolysis of Nucleic Acids in DNA and RNA Nucleases. , 2019, Journal of the American Chemical Society.

[18]  R. Guigó,et al.  The Origins and the Biological Consequences of the Pur/Pyr DNA·RNA Asymmetry , 2019, Chem.

[19]  C. Bennett Therapeutic Antisense Oligonucleotides Are Coming of Age. , 2019, Annual review of medicine.

[20]  M. De Vivo,et al.  The Catalytic Mechanism of DNA and RNA Polymerases , 2018, ACS Catalysis.

[21]  S. Solomon,et al.  Inotersen Treatment for Patients with Hereditary Transthyretin Amyloidosis , 2018, The New England journal of medicine.

[22]  H. Lönnberg,et al.  Phosphodiester models for cleavage of nucleic acids , 2018, Beilstein journal of organic chemistry.

[23]  P. Carloni,et al.  A Strategically Located Arg/Lys Residue Promotes Correct Base Paring During Nucleic Acid Biosynthesis in Polymerases. , 2018, Journal of the American Chemical Society.

[24]  Robert M. Dirks,et al.  RNA force field with accuracy comparable to state-of-the-art protein force fields , 2018, Proceedings of the National Academy of Sciences.

[25]  M. Marcia,et al.  Second-Shell Basic Residues Expand the Two-Metal-Ion Architecture of DNA and RNA Processing Enzymes , 2018, Structure.

[26]  R. J. Ramamurthi,et al.  Nusinersen versus Sham Control in Infantile‐Onset Spinal Muscular Atrophy , 2017, The New England journal of medicine.

[27]  Naoki Iwamoto,et al.  Control of phosphorothioate stereochemistry substantially increases the efficacy of antisense oligonucleotides , 2017, Nature Biotechnology.

[28]  S. Crooke Molecular Mechanisms of Antisense Oligonucleotides , 2017, Nucleic acid therapeutics.

[29]  Modesto Orozco,et al.  How accurate are accurate force-fields for B-DNA? , 2017, Nucleic acids research.

[30]  P. Carloni,et al.  A Self-Activated Mechanism for Nucleic Acid Polymerization Catalyzed by DNA/RNA Polymerases. , 2016, Journal of the American Chemical Society.

[31]  Mallory A. Havens,et al.  Splice-switching antisense oligonucleotides as therapeutic drugs , 2016, Nucleic acids research.

[32]  Modesto Orozco,et al.  Multiscale simulation of DNA. , 2016, Current opinion in structural biology.

[33]  M. dal Peraro,et al.  Cooperative motion of a key positively charged residue and metal ions for DNA replication catalyzed by human DNA Polymerase-η , 2016, Nucleic acids research.

[34]  Jordi Torres,et al.  BIGNASim: a NoSQL database structure and analysis portal for nucleic acids simulation data , 2015, Nucleic Acids Res..

[35]  D. Case,et al.  PARMBSC1: A REFINED FORCE-FIELD FOR DNA SIMULATIONS , 2015, Nature Methods.

[36]  Thomas E Cheatham,et al.  Improved Force Field Parameters Lead to a Better Description of RNA Structure. , 2015, Journal of chemical theory and computation.

[37]  C. Simmerling,et al.  ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. , 2015, Journal of chemical theory and computation.

[38]  Carsten Kutzner,et al.  Tackling Exascale Software Challenges in Molecular Dynamics Simulations with GROMACS , 2015, EASC.

[39]  A. Cavalli,et al.  Catalytic metal ions and enzymatic processing of DNA and RNA. , 2015, Accounts of chemical research.

[40]  Bert L. de Groot,et al.  pmx: Automated protein structure and topology generation for alchemical perturbations , 2014, J. Comput. Chem..

[41]  Fritz Eckstein,et al.  Phosphorothioates, essential components of therapeutic oligonucleotides. , 2014, Nucleic acid therapeutics.

[42]  Eric E. Swayze,et al.  Synthesis, biophysical properties and biological activity of second generation antisense oligonucleotides containing chiral phosphorothioate linkages , 2014, Nucleic acids research.

[43]  G. Hummer,et al.  Calcium Inhibition of Ribonuclease H1 Two-Metal Ion Catalysis , 2014, Journal of the American Chemical Society.

[44]  Pengfei Li,et al.  Taking into Account the Ion-induced Dipole Interaction in the Nonbonded Model of Ions. , 2014, Journal of chemical theory and computation.

[45]  J. Mendell,et al.  Eteplirsen for the treatment of Duchenne muscular dystrophy , 2013, Annals of neurology.

[46]  Carlos González,et al.  NAFlex: a web server for the study of nucleic acid flexibility , 2013, Nucleic Acids Res..

[47]  G. Watts,et al.  Mipomersen and other therapies for the treatment of severe familial hypercholesterolemia , 2012, Vascular health and risk management.

[48]  B. L. de Groot,et al.  Towards computional specificity screening of DNA-binding proteins , 2011, Nucleic acids research.

[49]  Christophe Blanchet,et al.  CURVES+ web server for analyzing and visualizing the helical, backbone and groove parameters of nucleic acid structures , 2011, Nucleic Acids Res..

[50]  Brigitta Elsässer,et al.  Atomistic details of the associative phosphodiester cleavage in human ribonuclease H. , 2010, Physical chemistry chemical physics : PCCP.

[51]  T. Wada,et al.  Stereocontrolled solid-phase synthesis of oligonucleoside H-phosphonates by an oxazaphospholidine approach. , 2009, Angewandte Chemie.

[52]  T. Wada,et al.  Solid-phase synthesis of stereoregular oligodeoxyribonucleoside phosphorothioates using bicyclic oxazaphospholidine derivatives as monomer units. , 2008, Journal of the American Chemical Society.

[53]  Matteo Dal Peraro,et al.  Phosphodiester cleavage in ribonuclease H occurs via an associative two-metal-aided catalytic mechanism. , 2008, Journal of the American Chemical Society.

[54]  T. Cheatham,et al.  Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations , 2008, The journal of physical chemistry. B.

[55]  F. J. Luque,et al.  Theoretical analysis of antisense duplexes: determinants of the RNase H susceptibility. , 2008, Journal of the American Chemical Society.

[56]  F. Javier Luque,et al.  Towards a molecular dynamics consensus view of B-DNA flexibility , 2008, Nucleic acids research.

[57]  Ivan S Ufimtsev,et al.  Quantum Chemistry on Graphical Processing Units. 1. Strategies for Two-Electron Integral Evaluation. , 2008, Journal of chemical theory and computation.

[58]  Wei Yang,et al.  Structure of human RNase H1 complexed with an RNA/DNA hybrid: insight into HIV reverse transcription. , 2007, Molecular cell.

[59]  P. Chacón,et al.  Thorough validation of protein normal mode analysis: a comparative study with essential dynamics. , 2007, Structure.

[60]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[61]  Adam W Van Wynsberghe,et al.  Comparison of mode analyses at different resolutions applied to nucleic acid systems. , 2005, Biophysical journal.

[62]  Modesto Orozco,et al.  Exploring the Essential Dynamics of B-DNA. , 2005, Journal of chemical theory and computation.

[63]  F. J. Luque,et al.  Structure, recognition properties, and flexibility of the DNA.RNA hybrid. , 2005, Journal of the American Chemical Society.

[64]  Modesto Orozco,et al.  Relative flexibility of DNA and RNA: a molecular dynamics study. , 2004, Journal of molecular biology.

[65]  G. B. Tafoya,et al.  Evaluation of 15 Polymerases and Phosphorothioate Primer Modification for Detection of UV-induced C:G to T:A Mutations by Allele-specific PCR¶ , 2004 .

[66]  W. Stec,et al.  Dynamic NMR structures of [Rp]- and [Sp]-phosphorothioated DNA-RNA hybrids: is flexibility required for RNase H recognition? , 2003, Biophysical journal.

[67]  W. Stec,et al.  Stereodefined phosphorothioate analogues of DNA: relative thermodynamic stability of the model PS-DNA/DNA and PS-DNA/RNA complexes. , 2002, Biochemistry.

[68]  S. Melchionna,et al.  Origin of the low-frequency modes of globular proteins. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[69]  T. James,et al.  Structural effect of complete [Rp]-phosphorothioate and phosphorodithioate substitutions in the DNA strand of a model antisense inhibitor-target RNA complex. , 1999, Journal of molecular biology.

[70]  T. Steitz,et al.  Structures of normal single-stranded DNA and deoxyribo-3'-S-phosphorothiolates bound to the 3'-5' exonucleolytic active site of DNA polymerase I from Escherichia coli. , 1999, Biochemistry.

[71]  G. Hessler,et al.  Structure of a stereoregular phosphorothioate DNA/RNA duplex , 1998, Nature Structural Biology.

[72]  T. Steitz,et al.  Structural principles for the inhibition of the 3'-5' exonuclease activity of Escherichia coli DNA polymerase I by phosphorothioates. , 1998, Journal of molecular biology.

[73]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997, J. Comput. Chem..

[74]  W Roush,et al.  Antisense Aims for a Renaissance , 1997, Science.

[75]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[76]  N. Go,et al.  Harmonicity and anharmonicity in protein dynamics: A normal mode analysis and principal component analysis , 1995, Proteins.

[77]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[78]  T. James,et al.  Structure and dynamics of a DNA.RNA hybrid duplex with a chiral phosphorothioate moiety: NMR and molecular dynamics with conventional and time-averaged restraints. , 1995, Biochemistry.

[79]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[80]  Kremer,et al.  Molecular dynamics simulation for polymers in the presence of a heat bath. , 1986, Physical review. A, General physics.

[81]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[82]  H R Drew,et al.  Structure of a B-DNA dodecamer: conformation and dynamics. , 1981, Proceedings of the National Academy of Sciences of the United States of America.

[83]  Levi N. Naden,et al.  Best Practices for Alchemical Free Energy Calculations [Article v1.0] , 2019, Living Journal of Computational Molecular Science.

[84]  M. Wood,et al.  Antisense oligonucleotides: the next frontier for treatment of neurological disorders , 2018, Nature Reviews Neurology.

[85]  F. J. Luque,et al.  The relative flexibility of B-DNA and A-RNA duplexes: database analysis. , 2004, Nucleic acids research.

[86]  W. Stec,et al.  Stereodifferentiation--the effect of P chirality of oligo(nucleoside phosphorothioates) on the activity of bacterial RNase H. , 1995, Nucleic acids research.