DNA backbone interactions impact the sequence specificity of DNA sulfur-binding domains: revelations from structural analyses

Harnessing sulfur-binding domains to separate Sp and Rp isomers of phosphorothioate oligonucleotides

Chemical synthesis of phosphoromonothioate oligonucleotides (PS-ONs) is not stereo-specific and produces a mixture of Rp and Sp diastereomers, whose disparate reactivity can complicate applications. Although the current methods to separate these diastereomers which rely on chromatography are constantly improving, many Rp and Sp diastereomers are still co-eluted. Here, based on sulfur-binding domains that specifically recognize phosphorothioated DNA and RNA in Rp configuration, we developed a universal separation system for phosphorothioate oligonucleotide isomers using immobilized SBD (SPOIS). With the scalable SPOIS, His-tagged SBD is immobilized onto Ni-nitrilotriacetic acid-coated magnetic beads to form a beads/SBD complex, Rp isomers of the mixture can be completely bound by SBD and separated from Sp isomers unbound in liquid phase, then recovered through suitable elution approach. Using the phosphoromonothioate single-stranded DNA as a model, SPOIS separated PS-ON diastereomers of 4 nt to 50 nt in length at yields of 60-90% of the starting Rp isomers, with PS linkage not locating at 5' or 3' end. Within this length range, PS-ON diastereomers that co-eluted in HPLC could be separated by SPOIS at yields of 84% and 89% for Rp and Sp stereoisomers, respectively. Furthermore, as each Rp phosphorothioate linkage can be bound by SBD, SPOIS allowed the separation of stereoisomers with multiple uniform Sp configurations for multiple phosphorothioate modifications. A second generation of SPOIS was developed using the thermolabile and non-sequence-specific SBDPed, enabling fast and high-yield recovery of PS substrate stereoisomers for the DNAzyme Cd16 and further demonstrating the efficiency of this method. KEY POINTS: • SPOIS allows isomer separations of the Rp and Sp isomers co-eluted on HPLC. • SPOIS can obtain Sp isomers with 5 min and Rp in 20 min from PS-ON diastereomers. • SPOIS was successfully applied to separate isomers of PS substrates of DNAzyme.

The binding affinity-dependent inhibition of cell growth and viability by DNA sulfur-binding domains

Increasing evidence suggests that DNA phosphorothioate (PT) modification serves several purposes in the bacterial host, and some restriction enzymes specifically target PT-DNA. PT-dependent restriction enzymes (PDREs) bind PT-DNA through their DNA sulfur binding domain (SBD) with dissociation constants (KD) of 5 nM~1 μM. Here, we report that SprMcrA, a PDRE, failed to dissociate from PT-DNA after cleavage due to high binding affinity, resulting in low DNA cleavage efficiency. Expression of SBDs in Escherichia coli cells with PT modification induced a drastic loss of cell viability at 25°C when both DNA strands of a PT site were bound, with one SBD on each DNA strand. However, at this temperature, SBD binding to only one PT DNA strand elicited a severe growth lag rather than lethality. This cell growth inhibition phenotype was alleviated by raising the growth temperature. An in vitro assay mimicking DNA replication and RNA transcription demonstrated that the bound SBD hindered the synthesis of new DNA and RNA when using PT-DNA as the template. Our findings suggest that DNA modification-targeting proteins might regulate cellular processes involved in DNA metabolism in addition to being components of restriction-modification systems and epigenetic readers.

Understanding base and backbone contributions of phosphorothioate DNA for molecular recognition with SBD proteins

Bacterial DNA phosphorothioate (PT) modification provides a specific anchoring site for sulfur-binding proteins (SBDs). Besides, their recognition patterns include phosphate links and bases neighboring the PT-modified site, thereby bringing about genome sequence-dependent properties in PT-related epigenetics. Here, we analyze the contributions of the DNA backbone (phosphates and deoxyribose) and bases bound with two SBD proteins in Streptomyces pristinaespiralis and coelicolor (SBDSco and SBDSpr). The chalcogen-hydrophobic interactions remained constantly at the anchoring site while the adjacent bases formed conditional and distinctive non-covalent interactions. More importantly, SBD/PT-DNA interactions were not limited within the traditional "4-bp core" range from 5'-I to 3'-III but extended to upstream 5'-II and 5'-III bases and even 5''-I to 5''-III at the non-PT-modified complementary strand. From the epigenetic viewpoint, bases 3'-II, 5''-I, and 5''-III of SBDSpr and 3'-II, 5''-II, and 5''-III of SBDSco present remarkable differentiations in the molecular recognitions. From the protein viewpoint, H102 in SBDSpr and R191 in SBDSco contribute significantly while proline residues at the PT-bound site are strictly conserved for the PT-chalcogen bond. The mutual and make-up mutations are proposed to alter the SBD/PT-DNA recognition pattern, besides additional chiral phosphorothioate modifications on phosphates 5'-II, 5'-II, 3'-I, and 3'-II.

Characterization of a promiscuous DNA sulfur binding domain and application in site-directed RNA base editing

Phosphorothioate (PT)-modification was discovered in prokaryotes and is involved in many biological functions such as restriction-modification systems. PT-modification can be recognized by the sulfur binding domains (SBDs) of PT-dependent restriction endonucleases, through coordination with the sulfur atom, accompanied by interactions with the DNA backbone and bases. The unique characteristics of PT recognition endow SBDs with the potential to be developed into gene-targeting tools, but previously reported SBDs display sequence-specificity for PT-DNA, which limits their applications. In this work, we identified a novel sequence-promiscuous SBDHga from Hahella ganghwensis. We solved the crystal structure of SBDHga complexed with PT-DNA substrate to 1.8 Å resolution and revealed the recognition mechanism. A shorter L4 loop of SBDHga interacts with the DNA backbone, in contrast with previously reported SBDs, which interact with DNA bases. Furthermore, we explored the feasibility of using SBDHga and a PT-oligonucleotide as targeting tools for site-directed adenosine-to-inosine (A-to-I) RNA editing. A GFP non-sense mutant RNA was repaired at about 60% by harnessing a chimeric SBD-hADAR2DD (deaminase domain of human adenosine deaminase acting on RNA), comparable with currently available RNA editing techniques. This work provides insights into understanding the mechanism of sequence-specificity for SBDs and for developing new tools for gene therapy.

A rapid nucleic acid detection platform based on phosphorothioate-DNA and sulfur binding domain

Nucleic acid detection plays a key role in diverse diagnosis and disease control. Currently available nucleic acid detection techniques are challenged by trade-offs among speed, simplicity, precision and cost. Here, we described a novel method, designated SENSOR (Sulfur DNA mediated nucleic acid sensing platform), for rapid nucleic acid detection. SENSOR was developed from phosphorothioate (PT)-DNA and sulfur binding domain (SBD) which specifically binds double-stranded PT-modified DNA. SENSOR utilizes PT-DNA oligo and SBD as targeting module, which is linked with split luciferase reporter to generate luminescence signal within 10 min. We tested detection on synthesized nucleic acid and COVID-19 pseudovirus, achieving attomolar sensitivity combined with an amplification procedure. Single nucleotide polymorphisms (SNP) could also be discriminated. Indicating SENSOR a new promising nucleic acid detection technique.

High-throughput sequencing of EcoWI restriction fragments maps the genome-wide landscape of phosphorothioate modification at base resolution

Phosphorothioation (PT), in which a non-bridging oxygen is replaced by a sulfur, is one of the rare modifications discovered in bacteria and archaea that occurs on the sugar-phosphate backbone as opposed to the nucleobase moiety of DNA. While PT modification is widespread in the prokaryotic kingdom, how PT modifications are distributed in the genomes and their exact roles in the cell remain to be defined. In this study, we developed a simple and convenient technique called EcoWI-seq based on a modification-dependent restriction endonuclease to identify genomic positions of PT modifications. EcoWI-seq shows similar performance than other PT modification detection techniques and additionally, is easily scalable while requiring little starting material. As a proof of principle, we applied EcoWI-seq to map the PT modifications at base resolution in the genomes of both the Salmonella enterica cerro 87 and E. coli expressing the dnd+ gene cluster. Specifically, we address whether the partial establishment of modified PT positions is a stochastic or deterministic process. EcoWI-seq reveals a systematic usage of the same subset of target sites in clones for which the PT modification has been independently established.

Large number of bacteria have modified their DNA mainly as part of a strategy to resist virus infection. Most of the modifications are chemical variations on the canonical bases A, T, C or G with phosphorothioate (PT) being a rare exception of a modification that happens on the backbone of the DNA. Interestingly, this PT modification was first chemically synthesized for specific biotechnological processes before scientists discovered that bacteria and archaea naturally perform this modification using their enzymes. The exact roles of phosphorothioation in bacteria and archaea is still under investigation. To enable further investigation of PT modifications, we designed EcoWI-seq, a method to identify the exact positions of these modifications in bacterial genomes. Notably, we applied the EcoWI-seq to several strains of E. coli for which PT modification has been induced by cloning into these strains, the necessary genes for making such modification. We found that these strains, despite being independently made, followed a precise pattern of PT modification with always the same sites being modified. This result indicates a deterministic process for the establishment of PT modification.

Origin of iodine preferential attack at sulfur in phosphorothioate and subsequent P-O or P-S bond dissociation

Iodine-induced cleavage at phosphorothioate DNA (PT-DNA) is characterized by extremely high sensitivity (∼1 phosphorothioate link per 106 nucleotides), which has been used for detecting and sequencing PT-DNA in bacteria. Despite its foreseeable potential for wide applications, the cleavage mechanism at the PT-modified site has not been well established, and it remains unknown as to whether or not cleavage of the bridging P-O occurs at every PT-modified site. In this work, we conducted accurate ωB97X-D calculations and high-performance liquid chromatography-mass spectrometry to investigate the process of PT-DNA cleavage at the atomic and molecular levels. We have found that iodine chemoselectively binds to the sulfur atom of the phosphorothioate link via a strong halogen-chalcogen interaction (a type of halogen bond, with binding affinity as high as 14.9 kcal/mol) and thus triggers P-O bond cleavage via phosphotriester-like hydrolysis. Additionally, aside from cleavage of the bridging P-O bond, the downstream hydrolyses lead to unwanted P-S/P-O conversions and a loss of the phosphorothioate handle. The mechanism we outline helps to explain specific selectivity at the PT-modified site but also predicts the dynamic stoichiometry of P-S and P-O bond breaking. For instance, Tris is involved in the cascade derivation of S-iodo-phosphorothioate to S-amino-phosphorothioate, suppressing the S-iodo-phosphorothioate hydrolysis to a phosphate diester. However, hydrolysis of one-third of the Tris-O-grafting phosphotriester results in unwanted P-S/P-O conversions. Our study suggests that bacterial DNA phosphorothioation may more frequently occur than previous bioinformatic estimations have predicted from iodine-induced deep sequencing data.

Molecular recognition between bacterial phosphorothioate DNA and sulfur-binding domain (SBD): competition between the water cage and chalcogen-hydrophobic packet

Bacterial DNA phosphorothioation (PT) physiologically and stereo-specifically replaces a non-bridging oxygen in a phosphate link with a sulfur atom, which can be recognized by a highly conserved sulfur-binding domain (SBD). Here we conducted thermodynamic integration (TI), molecular dynamics simulation, and quantum chemical calculations to decipher the specific molecular interactions between PT-DNA and SBD in Streptomyces coelicolor type IV restriction enzyme ScoMcrA. The TI-calculated binding affinity of (5'-CCG_Rp-PSGCCGG-3')₂ is larger than that of (5'-CCGGCCGG-3')₂ by about 7.4-7.7 kcal mol^-1. The binding difference dominantly stems from hydration energy of non-phosphorothioate DNA (9.8-10.6 kcal mol^-1) in aqueous solution, despite the persistent preference of 2.6-3.2 kcal mol^-1 in the DNA-SBD MD simulations. Furthermore, the quantum chemical calculations reveal an unusual non-covalent interaction in the phosphorothioate-binding scenario, where the ^PS⋯N^P165 chalcogen bond prevails the ^PS⋯HC_β vdW interactions from the adjacent residues H116-R117-Y164-P165-A168. Thus, the chalcogen-hydrophobic interaction pulls PT-DNA into the SBD binding pocket while the water cage pulls a normal DNA molecule out. The synergetic mechanism suggests the special roles of the proline pyrrolidine group in the SBD proteins, consistent with the experimental observations in the X-ray crystallography and structural bioinformatics analysis.

Phosphorothioate-DNA bacterial diet reduces the ROS levels in C. elegans while improving locomotion and longevity

DNA phosphorothioation (PT) is widely distributed in the human gut microbiome. In this work, PT-diet effect on nematodes was studied with PT-bioengineering bacteria. We found that the ROS level decreased by about 20–50% and the age-related lipofuscin accumulation was reduced by 15–25%. Moreover, the PT-feeding worms were more active at all life periods, and more resistant to acute stressors. Intriguingly, their lifespans were prolonged by ~21.7%. Comparative RNA-seq analysis indicated that many gene expressions were dramatically regulated by PT-diet, such as cysteine-rich protein (scl-11/12/13), sulfur-related enzyme (cpr-2), longevity gene (jnk-1) and stress response (sod-3/5, gps-5/6, gst-18/20, hsp-12.6). Both the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis suggested that neuroactivity pathways were upregulated, while phosphoryl transfer and DNA-repair pathways were down-regulated in good-appetite young worms. The findings pave the way for pro-longevity of multicellular organisms by PT-bacterial interference.

Qiang Huang et al. fed C. elegans with E. coli containing phosphorothioate (PT) DNA or a control strain and evaluated the impact on animal physiology. They observed that worms fed PT( + ) diets exhibited low reactive oxygen species, more active movement, and a longer lifespan compared to controls, suggesting that PT-DNA may have a positive effect on animal health.

Bacterial YedK represses plasmid DNA replication and transformation through its DNA single-strand binding activity

The SOS response-associated peptidase (SRAP) is an ancient protein superfamily in all domains of life. The mammalian SRAP was recently reported to covalently bind to the abasic sites (AP) in single stranded (ss) DNA to shield the chromosome integrity. YedK, the Escherichia coli SRAP, is not functionally characterized. Here we report the fortuitous pull-down of YedK from bacterial cell lysates by short (<20 bp) double stranded (ds) DNAs, further enrichment of YedK was observed when single stranded (ss) DNA was added. YedK can bind multiple DNA substrates, particularly with a high affinity to DNA duplex with single strand segment. As a SRAP protein, the involvement of YedK in SOS response was extensively examined, however yedK mutant of Escherichia coli showed no difference from the wild type strain upon the treatments with UV and various DNA damaging reagents, indicating its non-essentiality or redundancy in E. coli. Surprisingly, yedK mutants derived from Escherichia coli and Samonella enterica both showed an increased plasmid DNA transformation efficiency compared to the wild types. In accordance with this, induction of YedK effectively decreased the copy number of plasmid DNA. Site-directed mutagenesis of YedK demonstrated that residues involved in single strand DNA binding and cysteine residue at position 2 from N-terminus can discharge the repression of the plasmid transformation efficiency.

Characterization of BisI Homologs

BisI is a sequence-specific and 5-methylcytosine (m5C)-dependent restriction endonuclease (REase), that cleaves the modified DNA sequence Gm5CNGC (G indicates that the cytosine opposite to G is modified). We expressed and purified a number of BisI homologs from sequenced bacterial genomes and used Illumina sequencing to determine the Pam7902I (Esp638I-like) cleavage sites in phage Xp12 DNA. One BisI homolog KpnW2I is EcoBLMcrX-like, cleaving GCNGC/RCNGY/RCNRC sites with m5C. We also cloned and expressed three BisI homologs from metagenome sequences derived from thermophilic sources. One enzyme EsaTMI is active at 37 to 65°C. EsaHLI cleaves GCNGC sites with three to four m5C and is active up to 50°C. In addition, we determined the number and position of m5C in BisI sites for efficient cleavage. BisI cleavage efficiency of GCNGC site is as following: Gm5CNGC (two internal m5C) > Gm5CNGC (one internal m5C) > GCNGm5C (one external m5C) > > GCNGC (unmodified). Three or four m5C in GCNGC site also supports BisI cleavage although partial inhibition was observed on duplex oligos with four m5C. BisI can be used to partially cleave a desired GCNGC site targeted with a complementary oligonucleotide (hemi-methylated). The m5C-dependent BisI variants will be useful for epigenetic research.

Ultrasensitive and easily reproducible biosensor based on novel doped MoS2 nanowires field-effect transistor in label-free approach for detection of hepatitis B virus in blood serum

An ultrasensitive field-effect transistor (FET) for hepatitis B virus deoxyribonucleic acid (HBV DNA) detection in label free approach and easily reproducible setup was reported. The fabricated FET biosensor was materialized by ZnO doped MoS₂ nanowires (NWs). This report introduced a novel structure of the MoS₂ in bio-sensing approach. Because of unique electrical and structural properties of MoS₂, HBV biosensor could demonstrate the high sensitivity and showed the detection limit of 1 fM. The MoS₂ NWs fabrication was materialized through ZnO based vapor-liquid-solid (VLS) technique. The fabricated device could measure the DNA targets in a linear concentration range from 0.5 pM to 50 μM. The dynamic response time of FET biosensor was 25 s. The functionality of the NWs biosensor for label-free measurements could be repeated for several times without any significant malfunction and biosensor could retain 96% of its initial response after eight weeks maintenance. The HBV biosensor showed high selectivity by discrimination the complementary DNA oligonucleotides from non-complementary and the mismatch (1, 2 and 3 bases) oligonucleotides. The materialized platform was desirably reproduced for HBV concentrations in human serum. The specificity of the biosensor was evaluated against several different types of DNAs and the fabricated device showed the outstanding performance. In order to optimize the device functionality, the biosensor was checked for two different human samples and device could distinguish the samples from each other in the same manner.

Protein Domain Guided Screen for Sequence Specific and Phosphorothioate-Dependent Restriction Endonucleases

Modification dependent restriction endonucleases (MDREs) restrict modified DNA, typically with limited sequence specificity (∼2-4 bp). Here, we focus on MDREs that have an SRA and/or SBD (sulfur binding domain) fused to an HNH endonuclease domain, cleaving cytosine modified or phosphorothioated (PT) DNA. We independently characterized the SBD-SRA-HNH endonuclease ScoMcrA, which preferentially cleaves 5hmC modified DNA. We report five SBD-HNH endonucleases, all recognizing GpsAAC/GpsTTC sequence and cleaving outside with a single nucleotide 3' stagger: EcoWI (N7/N6), Ksp11411I (N5/N4), Bsp305I (N6/N4-5), Mae9806I [N(8-10)/N(8-9)], and Sau43800I [N(8-9)/N(7-8)]. EcoWI and Bsp305I are more specific for PT modified DNA in Mg2+ buffer, and promiscuous with Mn2+. Ksp11411I is more PT specific with Ni2+. EcoWI and Ksp11411I cleave fully- and hemi-PT modified oligos, while Bsp305I cleaves only fully modified ones. EcoWI forms a dimer in solution and cleaves more efficiently in the presence of two modified sites. In addition, we demonstrate that EcoWI PT-dependent activity has biological function: EcoWI expressing cells restrict dnd+ GpsAAC modified plasmid strongly, and GpsGCC DNA weakly. This work establishes a framework for biotechnology applications of PT-dependent restriction endonucleases (PTDRs).