Structural investigation into physiological DNA phosphorothioate modification

Crystal structures of monomeric BsmI restriction endonuclease reveal coordinated sequential cleavage of two DNA strands

BsmI, a thermophilic Type IIS restriction endonuclease from Bacillus stearothermophilus, presents a unique structural composition, housing two distinct active sites within a single monomer. Recognition of the non-symmetrical 5'-GAATGC-3' sequence enables precise cleavage of the top and bottom DNA strands. Synthetic biology interventions have led to the transformation of BsmI into Nb.BsmI, a nicking endonuclease. Here we introduce Nt*.BsmI, tailored for top-strand cleavage, which is inactive on standard double-stranded DNA, but active on bottom-strand nicked DNA, suggesting a sequential cleavage mechanism. Crystallographic structures of pre- and post-reactive complexes with cognate DNA show one major conformational change, a retractable loop possibly governing sequential active site accessibility. The x-ray structures reveal the position of the divalent metal ions in the active sites and the DNA:protein interactions, while the models predicted by Alphafold3 are incorrect. This comprehensive structural and functional study lays a foundation for rational enzyme redesign and potential applications in biotechnology.

Unravelling the Link between Oligonucleotide Structure and Diastereomer Separation in Hydrophilic Interaction Chromatography

Therapeutic oligonucleotides (ONs) commonly incorporate phosphorothioate (PS) modifications. These introduce chiral centers and generate ON diastereomers. The increasing number of ONs undergoing clinical trials and reaching the market has led to a growing interest to better characterize the ON diastereomer composition, especially for small interfering ribonucleic acids (siRNAs). In this study, and for the first time, we identify higher-order structures as the major cause of ON diastereomer separation in hydrophilic interaction chromatography (HILIC). We have used conformational predictions and melting profiles of several representative full-length ONs to first analyze ON folding and then run mass spectrometry and HILIC to underpin the link between their folding and diastereomer separation. On top, we show how one can either enhance or suppress diastereomer separation depending on chromatographic settings, such as column temperature, pore size, stationary phase, mobile-phase ionic strength, and organic modifier. This work will significantly facilitate future HILIC-based characterization of PS-containing ONs; e.g., enabling monitoring of batch-to-batch diastereomer distributions in full-length siRNAs, a complex task that is now for the first time shown as possible on this delicate class of therapeutic double-stranded ONs.

Analytical techniques for characterizing diastereomers of phosphorothioated oligonucleotides

Oligonucleotides have emerged as powerful therapeutics for treating diverse diseases. To fully unlock the therapeutic potential of oligonucleotides, there is still a great need to further improve their drug-like properties. Numerous chemical modifications have been explored to achieve this goal, with phosphorothioation being one of the most widely used strategies. However, phosphorothioate modification produces diastereomers that are reported to have different properties and performances, demanding detailed characterization of these diastereomers. Here we provide an overview of phosphorothioated oligonucleotide diastereomers, covering their origin and configurations, physicochemical and pharmacological properties, and stereo-selective chemical synthesis, followed by a summary of currently available analytical techniques for characterizing these diastereomers, with a focus on liquid chromatography-based approaches, including ion-pair reversed-phase liquid chromatography, anion exchange chromatography, mixed-mode chromatography, and hybrid approaches. Non-chromatographic techniques, such as capillary electrophoresis, spectroscopy and other methods, are also being reviewed.

Antisense oligonucleotides: absorption, distribution, metabolism, and excretion

INTRODUCTION

Antisense oligonucleotides (ASOs) have emerged as a promising novel drug modality that aims to address unmet medical needs. A record of six ASO drugs have been approved since 2016, and more candidates are in clinical development. ASOs are the most advanced class within the RNA-based therapeutics field.

AREAS COVERED

This review highlights the two major backbones that are currently used to build the most advanced ASO platforms - the phosphorodiamidate morpholino oligomers (PMOs) and the phosphorothioates (PSs). The absorption, distribution, metabolism, and excretion (ADME) properties of the PMO and PS platforms are discussed in detail.

EXPERT OPINION

Understanding the ADME properties of existing ASOs can foster further improvement of this cutting-edge therapy, thereby enabling researchers to safely develop ASO drugs and enhancing their ability to innovate.

ABBREVIATIONS

2'-MOE, 2'-O-methoxyethyl; 2'PS, 2 modified PS; ADME, absorption, distribution, metabolism, and excretion; ASO, antisense oligonucleotide; AUC, area under the curve; BNA, bridged nucleic acid; CPP, cell-penetrating peptide; CMV, cytomegalovirus; CNS, central nervous system; CYP, cytochrome P; DDI, drug-drug interaction; DMD, Duchenne muscular dystrophy; FDA, Food and Drug Administration; GalNAc3, triantennary N-acetyl galactosamine; IT, intrathecal; IV, intravenous; LNA, locked nucleic acid; mRNA, messenger RNA; NA, not applicable; PBPK, physiologically based pharmacokinetics; PD, pharmacodynamic; PK, pharmacokinetic; PMO, phosphorodiamidate morpholino oligomer; PMOplus, PMOs with positionally specific positive molecular charges; PPMO, peptide-conjugated PMO; PS, phosphorothioate; SC, subcutaneous; siRNA, small-interfering RNA; SMA, spinal muscular atrophy.

Multipolar Description of Atom-Atom Electrostatic Interaction Energies in Single/Double-Stranded DNAs

Molecular force field simulation is an effective method to explore the properties of DNA molecules in depth. Almost all current popular force fields calculate atom-atom electrostatic interaction energies for DNAs based on the atomic charge and dipole or quadrupole moments, without considering high-rank atomic multipole moments for more accurate electrostatics. Actually, the distribution of electrons around atomic nuclei is not spherically symmetric but is geometry dependent. In this work, a multipole expansion method that allows us to combine polarizability and anisotropy was applied. One single-stranded DNA and one double-stranded DNA were selected as pilot systems. Deoxynucleotides were cut out from pilot systems and capped by mimicking the original DNA environment. Atomic multipole moments were integrated instead of fixed-point charges to calculate atom-atom electrostatic energies to improve the accuracy of force fields for DNA simulations. Also, the applicability of modeling the behavior of both single-stranded and double-stranded DNAs was investigated. The calculation results indicated that the models can be transferred from pilot systems to test systems, which is of great significance for the development of future DNA force fields.

Different reactivity of Sp and Rp isomers of phosphorothioate-modified oligonucleotides in a duplex structure

The presence of a stereoisomeric center at the phosphorus atom in phosphorothioate-modified oligonucleotides (PS-ONs) has been recognized as an important feature since the early stages of their development. Therefore, several studies have been conducted on the chirality of PS-ONs. In this study, we evaluated the stereo-biased chemistry of PS-ON duplexes. Depending on their absolute configurations, PS-ON duplexes were found to have significantly different and stereospecific reactivities towards simple alkylating reagent.

Molecular Dynamics Study of the Hybridization between RNA and Modified Oligonucleotides

MicroRNAs (miRNAs) are attractive drug candidates for many diseases as they can modulate the expression of gene networks. Recently, we discovered that DNAs targeting microRNA-22-3p (miR-22-3p) hold the potential for treating obesity and related metabolic disorders (type 2 diabetes mellitus, hyperlipidemia, and nonalcoholic fatty liver disease (NAFLD)) by turning fat-storing white adipocytes into fat-burning adipocytes. In this work, we explored the effects of chemical modifications, including phosphorothioate (PS), locked nucleic acid (LNA), and peptide nucleic acid (PNA), on the structure and energy of DNA analogs by using molecular dynamics (MD) simulations. To achieve a reliable prediction of the hybridization free energy, the AMOEBA polarizable force field and the free energy perturbation technique were employed. The calculated hybridization free energies are generally compatible with previous experiments. For LNA and PNA, the enhanced duplex stability can be explained by the preorganization mechanism, i.e., the single strands adopt stable helical structures similar to those in the duplex. For PS, the S and R isomers (Sp and Rp) have preferences for C2'-endo and C3'-endo sugar puckering conformations, respectively, and therefore Sp is less stable than Rp in DNA/RNA hybrids. In addition, the solvation penalty of Rp accounts for its destabilization effect. PS-LNA is similar to LNA as the sugar puckering is dominated by the locked sugar ring. This work demonstrated that MD simulations with polarizable force fields are useful for the understanding and design of modified nucleic acids.

Investigating the interactions between DNA and DndE during DNA phosphorothioation

The DNA phosphorothioate modification is a novel physiological variation in bacteria. DndE controls this modification by binding to dsDNA via a mechanism that remains unclear. Structural analysis of the wild-type DndE tetramer suggests that a positively charged region in its center is important for DNA binding. In the present study, we replaced residues G21 and G24 in this region with lysines, which increases the DNA binding affinity but does not affect the DNA degradation phenotype. Structural analysis of the mutant indicates that it forms a new tetrameric conformation and that DndE interacts with DNA as a monomer rather than as a tetramer. A structural model of the DndE-DNA complex, based on its structural homolog P22 Arc repressor, indicates that two flexible loops in DndE are determinants of DNA binding.

DNA phosphorothioate modification-a new multi-functional epigenetic system in bacteria

Synthetic phosphorothioate (PT) internucleotide linkages, in which a nonbridging oxygen is replaced by a sulphur atom, share similar physical and chemical properties with phosphodiesters but confer enhanced nuclease tolerance on DNA/RNA, making PTs a valuable biochemical and pharmacological tool. Interestingly, PT modification was recently found to occur naturally in bacteria in a sequence-selective and RP configuration-specific manner. This oxygen-sulphur swap is catalysed by the gene products of dndABCDE, which constitute a defence barrier with DndFGH in some bacterial strains that can distinguish and attack non-PT-modified foreign DNA, resembling DNA methylation-based restriction-modification (R-M) systems. Despite their similar defensive mechanisms, PT- and methylation-based R-M systems have evolved to target different consensus contexts in the host cell because when they share the same recognition sequences, the protective function of each can be impeded. The redox and nucleophilic properties of PT sulphur render PT modification a versatile player in the maintenance of cellular redox homeostasis, epigenetic regulation and environmental fitness. The widespread presence of dnd systems is considered a consequence of extensive horizontal gene transfer, whereas the lability of PT during oxidative stress and the susceptibility of PT to PT-dependent endonucleases provide possible explanations for the ubiquitous but sporadic distribution of PT modification in the bacterial world.

Investigation of factors influencing the separation of diastereomers of phosphorothioated oligonucleotides

This study presents a systematic investigation of factors influencing the chromatographic separation of diastereomers of phosphorothioated pentameric oligonucleotides as model solutes. Separation was carried out under ion-pairing conditions using an XBridge C₁₈ column. For oligonucleotides with a single sulfur substitution, the diastereomer selectivity was found to increase with decreasing carbon chain length of the tertiary alkylamine used as an ion-pair reagent. Using an ion-pair reagent with high selectivity for diastereomers, triethylammonium, it was found the selectivity increased with decreased ion-pair concentration and shallower gradient slope. Selectivity was also demonstrated to be dependent on the position of the modified linkage. Substitutions at the center of the pentamer resulted in higher diastereomer selectivity compared to substitutions at either end. For mono-substituted oligonucleotides, the retention order and stereo configuration were consistently found to be correlated, with Rp followed by Sp, regardless of which linkage was modified. The type of nucleobase greatly affects the observed selectivity. A pentamer of cytosine has about twice the diastereomer selectivity of that of thymine. When investigating the retention of various oligonucleotides eluted using tributylammonium as the ion-pairing reagent, no diastereomer selectivity could be observed. However, retention was found to be dependent on both the degree and position of sulfur substitution as well as on the nucleobase. When analyzing fractions collected in the front and tail of overloaded injections, a significant difference was found in the ratio between Rp and Sp diastereomers, indicating that the peak broadening observed when using tributylammonium could be explained by partial diastereomer separation.

Diastereomers of a mono-substituted phosphoryl guanidine trideoxyribonucleotide: Isolation and properties

Recently, a new type of nucleic acid analogues with modified phosphate group, namely, phosphoryl guanidine oligonucleotides, has been described. In the present work, we assess the difference between diastereomers of a mono-substituted phosphoryl guanidine oligonucleotide and analyze their resistance to nuclease digestion. Individual diastereomers ('fast' and 'slow') of a trideoxynucleotide d (TpCp*A) were isolated by reverse-phase HPLC. Snake venom phosphodiesterase digestion showed that the native trideoxynucleotide was fully degraded after 30 min, whereas both 'fast' and 'slow' diastereomers of d (TpCp*A) were not completely digested even after 7 days. UV and CD spectra revealed similarities in the structure of the diastereomers. Structural analysis by 1D and 2D NMR spectroscopy also uncovered significant similarity in the properties of Rp and Sp diastereomers. Structural analysis of nuclear Overhauser effect spectroscopy (NOESY) data and restrained molecular dynamics methods showed very flexible single-stranded oligonucleotide structures. Detailed computational analysis of restraint penalty energies via restrained molecular dynamics simulations with the 2D NMR interproton distance data allowed us to conclude that most likely, the 'fast' isomer is the Sp diastereomer, and the 'slow' isomer is the Rp diastereomer.

Structural basis for the recognition of sulfur in phosphorothioated DNA

There have been very few reports on protein domains that specifically recognize sulfur. Here we present the crystal structure of the sulfur-binding domain (SBD) from the DNA phosphorothioation (PT)-dependent restriction endonuclease ScoMcrA. SBD contains a hydrophobic surface cavity that is formed by the aromatic ring of Y164, the pyrolidine ring of P165, and the non-polar side chains of four other residues that serve as lid, base, and wall of the cavity. The SBD and PT-DNA undergo conformational changes upon binding. The S¹⁸⁷RGRR¹⁹¹ loop inserts into the DNA major groove to make contacts with the bases of the G_PSGCC core sequence. Mutating key residues of SBD impairs PT-DNA association. More than 1000 sequenced microbial species from fourteen phyla contain SBD homologs. We show that three of these homologs bind PT-DNA in vitro and restrict PT-DNA gene transfer in vivo. These results show that SBD-like PT-DNA readers exist widely in prokaryotes.

P-Stereodefined phosphorothioate analogs of glycol nucleic acids-synthesis and structural properties

Enantiomerically pure, protected acyclic nucleosides of the GNA type (glycol nucleic acids) (GN'), obtained from (R)-(+)- and (S)-(-)-glycidols and the four canonical DNA nucleobases (Ade, Cyt, Gua and Thy), were transformed into 3'-O-DMT-protected 2-thio-4,4-pentamethylene-1,3,2-oxathiaphospholane derivatives (OTP-GN') containing a second stereogenic center at the phosphorus atom. These monomers were chromatographically separated into P-diastereoisomers, which were further used for the synthesis of P-stereodefined "dinucleoside" phosphorothioates GNPST (GN = GA, GC, GG, GT). The absolute configuration at the phosphorus atom for all eight GNPST was established enzymatically and verified chemically, and correlated with chromatographic mobility of the OTP-GN' monomers on silica gel. The GNPS units (derived from (R)-(+)-glycidol) were introduced into self-complementary PS-(DNA/GNA) octamers of preselected, uniform absolute configuration at P-atoms. Thermal dissociation experiments showed that the thermodynamic stability of the duplexes depends on the stereochemistry of the phosphorus centers and relative arrangement of the GN units in the oligonucleotide strands. These results correlate with the changes of conformation assessed from circular dichroism spectra.

HELM Software for Biopolymers

Hierarchical Editing Language for Macromolecules (HELM version 2.0) is a molecular line notation similar to SMILEs but specifically for communicating and managing biopolymer structures. The HELM project, part of the Pistoia Alliance nonprofit organization, has been tasked to develop and promote HELM as a global exchange format and recently released version 2.0 of the specification. Here we will describe the specifics of the HELM v2.0 notation along with the large ecosystem of software to support HELM-based structure management. We will highlight a recent open-source software and database for HELM monomers and a new, simpler approach to deploying a large complicated molecular management system.

Mechanistic Investigation on ROS Resistance of Phosphorothioated DNA

Phosphorothioated DNA (PT-DNA) exhibits a mild anti-oxidant property both in vivo and in vitro. It was found that 8-OHdG and ROS levels were significantly lower in dnd+ (i.e. S⁺) E. coli., compared to a dnd− (i.e. S⁻) strain. Furthermore, different from traditional antioxidants, phosphorothioate compound presents an unexpectedly high capacity to quench hydroxyl radical. Oxidative product analysis by liquid chromatography-mass spectrometry and quantum mechanistic computation supported its unique anti-oxidant characteristic of the hydroxyl selectivity: phosphorothioate donates an electron to either hydroxyl radical or guanine radical derived from hydroxyl radical, leading to a PS^• radical; a complex of PS^• radical and OH⁻ (i.e. the reductive product of hydroxyl radical) releases a highly reductive HS^• radical, which scavenges more equivalents of oxidants in the way to high-covalent sulphur compounds such as sulphur, sulphite and sulphate. The PS-PO conversion (PS and PO denote phosphorus-sulphur and phosphorus-oxygen compounds, respectively) made a switch of extremely oxidative OH^• to highly reductive HS^• species, endowing PT-DNA with the observed high capacity in hydroxyl-radical neutralization. This plausible mechanism provides partial rationale as to why bacteria develop the resource-demanding PT modification on guanine-neighboring phosphates in genome.

DNA Phosphorothioate Modification Plays a Role in Peroxides Resistance in Streptomyces lividans

DNA phosphorothioation, conferred by dnd genes, was originally discovered in the soil-dwelling bacterium Streptomyces lividans, and thereafter found to exist in various bacterial genera. However, the physiological significance of this sulfur modification of the DNA backbone remains unknown in S. lividans. Our studies indicate that DNA phosphorothioation has a major role in resistance to oxidative stress in the strain. Although Streptomyces species express multiple catalase/peroxidase and organic hydroperoxide resistance genes to protect them against peroxide damage, a wild type strain of S. lividans exhibited two-fold to 10-fold higher survival, compared to a dnd⁻ mutant, following treatment with peroxides. RNA-seq experiments revealed that, catalase and organic hydroperoxide resistance gene expression were not up-regulated in the wild type strain, suggesting that the resistance to oxidative stress was not due to the up-regulation of these genes by DNA phosphorothioation. Quantitative RT-PCR analysis was conducted to trace the expression of the catalase and the organic hydroperoxide resistance genes after peroxides treatments. A bunch of these genes were activated in the dnd⁻ mutant rather than the wild type strain in response to peroxides. Moreover, the organic hydroperoxide peracetic acid was scavenged more rapidly in the presence than in the absence of phosphorothioate modification, both in vivo and in vitro. The dnd gene cluster can be up-regulated by the disulfide stressor diamide. Overall, our observations suggest that DNA phosphorothioate modification functions as a peroxide resistance system in S. lividans.