Roles of metal ions in nucleases

Elucidation of the catalytic mechanism of a single-metal dependent homing endonuclease using QM and QM/MM approaches: the case study of I-PpoI

Homing endonucleases (HEs) are highly specific DNA cleaving enzymes, with I-PpoI having been suggested to use a single metal to accelerate phosphodiester bond cleavage. Although an I-PpoI mechanism has been proposed based on experimental structural data, no consensus has been reached regarding the roles of the metal or key active site amino acids. This study uses QM cluster and QM/MM calculations to provide atomic-level details of the I-PpoI catalytic mechanism. Minimal QM cluster and large-scale QM/MM models demonstrate that the experimentally-proposed pathway involving direct Mg2+ coordination to the substrate coupled with leaving group protonation through a metal-activated water is not feasible due to an inconducive I-PpoI active site alignment. Despite QM cluster models of varying size uncovering a pathway involving leaving group protonation by a metal-activated water, indirect (water-mediated) metal coordination to the substrate is required to afford this pathway, which renders this mechanism energetically infeasible. Instead, QM cluster models reveal that the preferred pathway involves direct Mg2+-O3' coordination to stabilize the charged substrate and assist leaving group departure, while H98 activates the water nucleophile. These calculations also underscore that both catalytic residues that directly interact with the substrate and secondary amino acids that position or stabilize these residues are required for efficient catalysis. QM/MM calculations on the solvated enzyme-DNA complex verify the preferred mechanism, which is fully consistent with experimental kinetic, structural, and mutational data. The fundamental understanding of the I-PpoI mechanism of action, gained from the present work can be used to further explore potential uses of this enzyme in biotechnology and medicine, and direct future computational investigations of other members of the understudied HE family.

Is Metal Stabilization of the Leaving Group Required or Can Lysine Facilitate Phosphodiester Bond Cleavage in Nucleic Acids? A Computational Study of EndoV

Endonuclease V (EndoV) is a single-metal-dependent enzyme that repairs deaminated DNA nucleobases in cells by cleaving the phosphodiester bond, and this enzyme has proven to be a powerful tool in biotechnology and medicine. The catalytic mechanism used by EndoV must be understood to design new disease detection and therapeutic solutions and further exploit the enzyme in interdisciplinary applications. This study has used a mixed molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) approach to compare eight distinct catalytic pathways and provides the first proposed mechanism for bacterial EndoV. The calculations demonstrate that mechanisms involving either direct or indirect metal coordination to the leaving group of the substrate previously proposed for other nucleases are unlikely for EndoV, regardless of the general base (histidine, aspartate, and substrate phosphate moiety). Instead, distinct catalytic pathways are characterized for EndoV that involve K139 stabilizing the leaving group, a metal-coordinated water stabilizing the transition structure, and either H214 or a substrate phosphate group activating the water nucleophile. In silico K139A and H214A mutational results support the newly proposed roles of these residues. Although this is a previously unseen combination of general base, general acid, and metal-binding architecture for a one-metal-dependent endonuclease, our proposed catalytic mechanisms are fully consistent with experimental kinetic, structural, and mutational data. In addition to substantiating a growing body of literature, suggesting that one metal is enough to catalyze P-O bond cleavage in nucleic acids, this new fundamental understanding of the catalytic function will promote the exploration of new and improved applications of EndoV.

An intracellular, non-oxidative factor activates in vitro chromatin fragmentation in pig sperm

BACKGROUND

In vitro incubation of epididymal and vas deferens sperm with Mn2+ induces Sperm Chromatin Fragmentation (SCF), a mechanism that causes double-stranded breaks in toroid-linker regions (TLRs). Whether this mechanism, thought to require the participation of topoisomerases and/or DNAses and thus far only described in epididymal mouse sperm, can be triggered in ejaculated sperm is yet to be elucidated. The current study aimed to determine if exposure of pig ejaculated sperm to divalent ions (Mn2+ and Mg2+) activates SCF, and whether this has any impact on sperm function and survival. For this purpose, sperm DNA integrity was evaluated through the Comet assay and Pulsed Field Gel Electrophoresis (PFGE); sperm motility and agglutination were assessed with computer assisted sperm analysis (CASA); and sperm viability and levels of total reactive oxygen species (ROS) and superoxides were determined through flow cytometry.

RESULTS

Incubation with Mn2+/Ca2+ activated SCF in a dose-dependent (P < 0.05) albeit not time-dependent manner (P > 0.05); in contrast, Mg2+/Ca2+ only triggered SCF at high concentrations (50 mM). The PFGE revealed that, when activated by Mn2+/Ca2+ or Mg2+/Ca2+, SCF generated DNA fragments of 33-194 Kb, compatible with the size of one or multiple toroids. Besides, Mn2+/Ca2+ affected sperm motility in a dose-dependent manner (P < 0.05), whereas Mg2+/Ca2+ only impaired this variable at high concentrations (P < 0.05). While this effect on motility was concomitant with an increase of agglutination, neither viability nor ROS levels were affected by Mn2+/Ca2+ or Mg2+/Ca2+ treatments.

CONCLUSION

Mn2+/Ca2+ and Mn2+/Ca2+ were observed to induce SCF in ejaculated sperm, resulting in DNA cleavage at TLRs. The activation of this mechanism by an intracellular, non-oxidative factor sheds light on the events taking place during sperm cell death.

Secreting recombinant barnase by Lactococcus lactis and its application in reducing RNA from forages

Barnase is a ribonuclease used for plasmid purification, targeted gene therapy and studies of protein interactions. To make the use of barnase easier, the barnase gene from Bacillus amyloliquefaciens BH072 was cloned into Lactococcus lactis under the control of the P_P5 or P_nisA promoters. Four recombinant expression vectors were constructed with one or two signal peptides to control the enzyme secretion. 310 mg/L barnase was obtained in the presence of its inhibitor barstar after 36 h induction. The properties of barnase were investigated, showing that the optimal reaction temperature and pH were 50 °C and 5.0, respectively, and the highest enzyme activity reached 16.5 kU/mL. Barnase stored at 40 °C for 72 h retained 90 % of its initial activity, and maintained more than 80 % of its initial activity after 72 h of storage at pH 5.0-9.0. Furthermore, the optimal conditions for enzymatic reduction of nucleic acids in single-cell proteins (SCP) forages was investigated. 1 % salt solution with an SCP-enzyme ratio of 1000:1, pH 5.0 and incubated at 50 °C for 1 h, allowed 82 % RNA content reduction. Finally, homology modeling of barnase demonstrates its three-dimensional structure, and substrate simulation docking predicts key active residues as well as bonding patterns.

Selectivity in Enzymatic Phosphorus Recycling from Biopolymers: Isotope Effect, Reactivity Kinetics, and Molecular Docking with Fungal and Plant Phosphatases

Among ubiquitous phosphorus (P) reserves in environmental matrices are ribonucleic acid (RNA) and polyphosphate (polyP), which are, respectively, organic and inorganic P-containing biopolymers. Relevant to P recycling from these biopolymers, much remains unknown about the kinetics and mechanisms of different acid phosphatases (APs) secreted by plants and soil microorganisms. Here we investigated RNA and polyP dephosphorylation by two common APs, a plant purple AP (PAP) from sweet potato and a fungal phytase from Aspergillus niger. Trends of δ¹⁸O values in released orthophosphate during each enzyme-catalyzed reaction in ¹⁸O-water implied a different extent of reactivity. Subsequent enzyme kinetics experiments revealed that A. niger phytase had 10-fold higher maximum rate for polyP dephosphorylation than the sweet potato PAP, whereas the sweet potato PAP dephosphorylated RNA at a 6-fold faster rate than A. niger phytase. Both enzymes had up to 3 orders of magnitude lower reactivity for RNA than for polyP. We determined a combined phosphodiesterase-monoesterase mechanism for RNA and terminal phosphatase mechanism for polyP using high-resolution mass spectrometry and ³¹P nuclear magnetic resonance, respectively. Molecular modeling with eight plant and fungal AP structures predicted substrate binding interactions consistent with the relative reactivity kinetics. Our findings implied a hierarchy in enzymatic P recycling from P-polymers by phosphatases from different biological origins, thereby influencing the relatively longer residence time of RNA versus polyP in environmental matrices. This research further sheds light on engineering strategies to enhance enzymatic recycling of biopolymer-derived P, in addition to advancing environmental predictions of this P recycling by plants and microorganisms.

Generating a Small Shuttle Vector for Effective Genetic Engineering of Methanosarcina mazei Allowed First Insights in Plasmid Replication Mechanism in the Methanoarchaeon

Due to their role in methane production, methanoarchaea are of high ecological relevance and genetic systems have been ever more established in the last two decades. The system for protein expression in Methanosarcina using a comprehensive shuttle vector is established; however, details about its replication mechanism in methanoarchaea remain unknown. Here, we report on a significant optimisation of the rather large shuttle vector pWM321 (8.9 kbp) generated by Metcalf through a decrease in its size by about 35% by means of the deletion of several non-coding regions and the ssrA gene. The resulting plasmid (pRS1595) still stably replicates in M. mazei and—most likely due to its reduced size—shows a significantly higher transformation efficiency compared to pWM321. In addition, we investigate the essential gene repA, coding for a rep type protein. RepA was heterologously expressed in Escherichia coli, purified and characterised, demonstrating the significant binding and nicking activity of supercoiled plasmid DNA. Based on our findings we propose that the optimised shuttle vector replicates via a rolling circle mechanism with RepA as the initial replication protein in Methanosarcina. On the basis of bioinformatic comparisons, we propose the presence and location of a double-strand and a single-strand origin, which need to be further verified.

Imprinted Polymeric Nanoparticles as Artificial Enzymes for Ester Hydrolysis at Room Temperature and pH 7

Natural esterases hydrolyze esters under physiological pHs but chemists often have to use strongly acidic or basic conditions for the same hydrolysis. We report synthetic nanoparticle catalysts that hydrolyze nonactivated alkyl esters at room temperature and neutral pH, with enzyme-like catalytic mechanisms and exquisite substrate selectivity. Unlike natural enzymes that denature easily at elevated temperatures, the synthetic catalysts become more active at higher temperatures.

Nucleic acid amplification-integrated single-molecule fluorescence imaging for in vitro and in vivo biosensing

Single-molecule fluorescence imaging is among the most advanced analytical technologies and has been widely adopted for biosensing due to its distinct advantages of simplicity, rapidity, high sensitivity, low sample consumption, and visualization capability. Recently, a variety of nucleic acid amplification approaches have been developed to provide a straightforward and highly efficient way for amplifying low abundance target signals. The integration of single-molecule fluorescence imaging with nucleic acid amplification has greatly facilitated the construction of various fluorescent biosensors for in vitro and in vivo detection of DNAs, RNAs, enzymes, and live cells with high sensitivity and good selectivity. Herein, we review the advances in the development of fluorescent biosensors by integrating single-molecule fluorescence imaging with nucleic acid amplification based on enzyme (e.g., DNA polymerase, RNA polymerase, exonuclease, and endonuclease)-assisted and enzyme-free (e.g., catalytic hairpin assembly, entropy-driven DNA amplification, ligation chain reaction, and hybridization chain reaction) strategies, and summarize the principles, features, and in vitro and in vivo applications of the emerging biosensors. Moreover, we discuss the remaining challenges and future directions in this area. This review may inspire the development of new signal-amplified single-molecule biosensors and promote their practical applications in fundamental and clinical research.

The core Cas1 protein of CRISPR-Cas I-B in Leptospira shows metal-tunable nuclease activity

Leptospira interrogans serovar Copenhageni strain Fiocruz L1-130 is the causative agent of leptospirosis in animals and humans. This organism carries a functional cas1 gene classified under CRISPR-Cas I-B. In this study, using various nuclease assays and bioinformatics analysis, we report that the recombinant Cas1 (LinCas1) possesses metal-ion dependent DNase activity, which is inhibited upon substitution or chelation of metal-ion and/or interaction with recombinant Cas2 (LinCas2) of L. interrogans. Model of LinCas1 structure shows a shorter N-terminal domain unlike other Cas1 orthologs reported to date. The C-terminal domain of LinCas1 contains conserved divalent-metal binding residues (Glu108, His176, and Glu191) and the mutation of these residues leads to abolition in DNase activity. Immunoassay using anti-LinCas2 demonstrates that LinCas1 interacts with LinCas2 and attains a saturation point. Moreover, the nuclease activity of the LinCas1-Cas2 mixture on ds-DNA displayed a reduction in activity compared to the pure core LinCas proteins under in vitro condition. The DNase activity for LinCas1 is consistent with a role for this protein in the recognition/cleavage of foreign DNA and integration of foreign DNA as spacer into the CRISPR array.

Toward Catalytic Antibiotics: Redesign of Fluoroquinolones to Catalytically Fragment Chromosomal DNA

A library of ciprofloxacin-nuclease conjugates was designed and synthesized to investigate their potential as catalytic antibiotics. The Cu(II) complexes of the new designer compounds (i) showed excellent in vitro hydrolytic and oxidative DNase activity, (ii) showed good antibacterial activity against both Gram-negative and Gram-positive bacteria, and (iii) proved to be highly potent bacterial DNA gyrase inhibitors via a mechanism that involves stabilization of the fluoroquinolone-topoisomerase-DNA ternary complex. Furthermore, the Cu(II) complexes of two of the new designer compounds were shown to fragment supercoiled plasmid DNA into linear DNA in the presence of DNA gyrase, demonstrating a "proof of concept" in vitro. These ciprofloxacin-nuclease conjugates can therefore serve as models with which to develop next-generation, in vivo functioning catalytic antimicrobials.

R, Sissi C, Pecina A, Riccardi L, De Vivo M, Mancin F, Scrimin P. A Gold Nanoparticle Nanonuclease Relying on a Zn(II) Mononuclear Complex

Similarly to enzymes, functionalized gold nanoparticles efficiently catalyze chemical reactions, hence the term nanozymes. Herein, we present our results showing how surface-passivated gold nanoparticles behave as synthetic nanonucleases, able to cleave pBR322 plasmid DNA with the highest efficiency reported so far for catalysts based on a single metal ion mechanism. Experimental and computational data indicate that we have been successful in creating a catalytic site precisely mimicking that suggested for natural metallonucleases relying on a single metal ion for their activity. It comprises one Zn(II) ion to which a phosphate diester of DNA is coordinated. Importantly, as in nucleic acids-processing enzymes, a positively charged arginine plays a key role by assisting with transition state stabilization and by reducing the pKa of the nucleophilic alcohol of a serine. Our results also show how designing a catalyst for a model substrate (bis-p-nitrophenylphosphate) may provide wrong indications as for its efficiency when it is tested against the real target (plasmid DNA).

Metal-Peptide Complexes as Promising Antibiotics to Fight Emerging Drug Resistance: New Perspectives in Tuberculosis

In metal-peptide interactions, cations form stable complexes through bonds with coordinating groups as side chains of amino acids. These compounds, among other things, exert a wide variety of antimicrobial activities through structural changes of peptides upon metal binding and redox chemistry. They exhibit different mechanisms of action (MOA), including the modification of DNA/RNA, protein and cell wall synthesis, permeabilization and modulation of gradients of cellular membranes. Nowadays, the large increase in antibiotic resistance represents a crucial problem to limit progression at the pandemic level of the diseases that seemed nearly eradicated, such as tuberculosis (Tb). Mycobacterium tuberculosis (Mtb) is intrinsically resistant to many antibiotics due to chromosomal mutations which can lead to the onset of novel strains. Consequently, the maximum pharmaceutical effort should be focused on the development of new therapeutic agents and antimicrobial peptides can represent a valuable option as a copious source of potential bioactive compounds. The introduction of a metal center can improve chemical diversity and hence specificity and bioavailability while, in turn, the coordination to peptides of metal complexes can protect them and enhance their poor water solubility and air stability: the optimization of these parameters is strictly required for drug prioritization and to obtain potent inhibitors of Mtb infections with novel MOAs. Here, we present a panoramic review of the most recent findings in the field of metal complex-peptide conjugates and their delivery systems with the potential pharmaceutical application as novel antibiotics in Mtb infections.

Factors Influencing the Activity of Nanozymes in the Cleavage of an RNA Model Substrate

A series of 2-nm gold nanoparticles passivated with different thiols all featuring at least one triazacyclonanone-Zn(II) complex and different flanking units (a second Zn(II) complex, a triethyleneoxymethyl derivative or a guanidinium of arginine of a peptide) were prepared and studied for their efficiency in the cleavage of the RNA-model substrate 2-hydroxypropyl-p-nitrophenyl phosphate. The source of catalysis for each of them was elucidated from the kinetic analysis (Michaelis–Menten profiles, pH dependence and kinetic isotope effect). The data indicated that two different mechanisms were operative: One involving two Zn(II) complexes and the other one involving a single Zn(II) complex and a flanking guanidinium cation. The mechanism based on a dinuclear catalytic site appeared more efficient than the one based on the cooperativity between a metal complex and a guanidinium.

Peroxidase from proso millet exhibits endonuclease-like activity

In this study, the mechanism of DNA cleavage by cationic peroxidase from proso millet (PmPOD) was investigated. PmPOD cleaved supercoiled circular DNA into both nicked circular and linear forms via a cleavage mechanism that resembles those of native endonucleases. Inhibition and ligation studies demonstrated that reactive oxygen species and the ferriprotoporphyrin IX moiety in PmPOD are not involved in PmPOD-mediated DNA cleavage. Similar to other endonucleases, Mg ions considerably enhance the DNA cleavage activity of PmPOD. Further studies suggested that PmPOD can disrupt phosphodiester bonds in DNA and mononucleotides, indicating that it is a phosphatase. The phosphatase activity of PmPOD is higher than that of horseradish peroxidase (HRP), but the peroxidase activity of PmPOD was lower than that of HRP. PmPOD-mediated hydrolytic cleavage of DNA observed in this study is different from those reported for heme proteins. This study provides valuable insights into the distinct mechanisms underlying DNA cleavage by heme proteins.

Metal complexes promoting catalytic cleavage of nucleic acids-biochemical tools and therapeutics

The development of metal complexes that promote degradation of nucleic acids has garnered significant interest as a result of their broad range of potential application. This review focuses on recent progress in the design and synthesis of metal complexes as artificial nucleases that promote either hydrolytic or oxidative cleavage of nucleic acids. Illustrative examples demonstrate the versatility of artificial nucleases for in vitro applications as molecular tools to address biochemical problems, as well as their potential use as therapeutic agents. We also address future challenges for improvement and avenues for further investigation.

Structure and dynamics of mesophilic variants from the homing endonuclease I-DmoI

I-DmoI, from the hyperthermophilic archaeon Desulfurococcus mobilis, belongs to the LAGLIDADG homing endonuclease protein family. Its members are highly specific enzymes capable of recognizing long DNA target sequences, thus providing potential tools for genome manipulation. Working towards this particular application, many efforts have been made to generate mesophilic variants of I-DmoI that function at lower temperatures than the wild-type. Here, we report a structural and computational analysis of two I-DmoI mesophilic mutants. Despite very limited structural variations between the crystal structures of these variants and the wild-type, a different dynamical behaviour near the cleavage sites is observed. In particular, both the dynamics of the water molecules and the protein perturbation effect on the cleavage site correlate well with the changes observed in the experimental enzymatic activity.

Structure-based domain assignment in Leishmania infantum EndoG: characterization of a pH-dependent regulatory switch and a C-terminal extension that largely dictates DNA substrate preferences

Mitochondrial endonuclease G from Leishmania infantum (LiEndoG) participates in the degradation of double-stranded DNA (dsDNA) during parasite cell death and is catalytically inactive at a pH of 8.0 or above. The presence, in the primary sequence, of an acidic amino acid-rich insertion exclusive to trypanosomatids and its spatial position in a homology-built model of LiEndoG led us to postulate that this peptide stretch might act as a pH sensor for self-inhibition. We found that a LiEndoG variant lacking residues 145–180 is indeed far more active than its wild-type counterpart at pH values >7.0. In addition, we discovered that (i) LiEndoG exists as a homodimer, (ii) replacement of Ser211 in the active-site SRGH motif with the canonical aspartate from the DRGH motif of other nucleases leads to a catalytically deficient enzyme, (iii) the activity of the S211D variant can be restored upon the concomitant replacement of Ala247 with Arg and (iv) a C-terminal extension is responsible for the observed preferential cleavage of single-stranded DNA (ssDNA) and ssDNA–dsDNA junctions. Taken together, our results support the view that LiEndoG is a multidomain molecular machine whose nuclease activity can be subtly modulated or even abrogated through architectural changes brought about by environmental conditions and interaction with other binding partners.

Comprehensive classification of the PIN domain-like superfamily

PIN-like domains constitute a widespread superfamily of nucleases, diverse in terms of the reaction mechanism, substrate specificity, biological function and taxonomic distribution. Proteins with PIN-like domains are involved in central cellular processes, such as DNA replication and repair, mRNA degradation, transcription regulation and ncRNA maturation. In this work, we identify and classify the most complete set of PIN-like domains to provide the first comprehensive analysis of sequence–structure–function relationships within the whole PIN domain-like superfamily. Transitive sequence searches using highly sensitive methods for remote homology detection led to the identification of several new families, including representatives of Pfam (DUF1308, DUF4935) and CDD (COG2454), and 23 other families not classified in the public domain databases. Further sequence clustering revealed relationships between individual sequence clusters and showed heterogeneity within some families, suggesting a possible functional divergence. With five structural groups, 70 defined clusters, over 100,000 proteins, and broad biological functions, the PIN domain-like superfamily constitutes one of the largest and most diverse nuclease superfamilies. Detailed analyses of sequences and structures, domain architectures, and genomic contexts allowed us to predict biological function of several new families, including new toxin-antitoxin components, proteins involved in tRNA/rRNA maturation and transcription/translation regulation.

Unraveling the essential role of CysK in CDI toxin activation

Contact-dependent growth inhibition (CDI) is a widespread mechanism of bacterial competition. CDI(+) bacteria deliver the toxic C-terminal region of contact-dependent inhibition A proteins (CdiA-CT) into neighboring target bacteria and produce CDI immunity proteins (CdiI) to protect against self-inhibition. The CdiA-CT(EC536) deployed by uropathogenic Escherichia coli 536 (EC536) is a bacterial toxin 28 (Ntox28) domain that only exhibits ribonuclease activity when bound to the cysteine biosynthetic enzyme O-acetylserine sulfhydrylase A (CysK). Here, we present crystal structures of the CysK/CdiA-CT(EC536) binary complex and the neutralized ternary complex of CysK/CdiA-CT/CdiI(EC536) CdiA-CT(EC536) inserts its C-terminal Gly-Tyr-Gly-Ile peptide tail into the active-site cleft of CysK to anchor the interaction. Remarkably, E. coli serine O-acetyltransferase uses a similar Gly-Asp-Gly-Ile motif to form the "cysteine synthase" complex with CysK. The cysteine synthase complex is found throughout bacteria, protozoa, and plants, indicating that CdiA-CT(EC536) exploits a highly conserved protein-protein interaction to promote its toxicity. CysK significantly increases CdiA-CT(EC536) thermostability and is required for toxin interaction with tRNA substrates. These observations suggest that CysK stabilizes the toxin fold, thereby organizing the nuclease active site for substrate recognition and catalysis. By contrast, Ntox28 domains from Gram-positive bacteria lack C-terminal Gly-Tyr-Gly-Ile motifs, suggesting that they do not interact with CysK. We show that the Ntox28 domain from Ruminococcus lactaris is significantly more thermostable than CdiA-CT(EC536), and its intrinsic tRNA-binding properties support CysK-independent nuclease activity. The striking differences between related Ntox28 domains suggest that CDI toxins may be under evolutionary pressure to maintain low global stability.

Role of magnesium ions in the reaction mechanism at the interface between Tm1631 protein and its DNA ligand

A protein, Tm1631 from the hyperthermophilic organism Thermotoga maritima belongs to a domain of unknown function protein family. It was predicted that Tm1631 binds with the DNA and that the Tm1631–DNA complex is an endonuclease repair system with a DNA repair function (Konc et al. PLoS Comput Biol 9(11): e1003341, 2013). We observed that the severely bent, strained DNA binds to the protein for the entire 90 ns of classical molecular dynamics (MD) performed; we could observe no significant changes in the most distorted region of the DNA, where the cleavage of phosphodiester bond occurs. In this article, we modeled the reaction mechanism at the interface between Tm1631 and its proposed ligand, the DNA molecule, focusing on cleavage of the phosphodiester bond. After addition of two Mg²⁺ ions to the reaction center and extension of classical MD by 50 ns (totaling 140 ns), the DNA ligand stayed bolted to the protein. Results from density functional theory quantum mechanics/molecular mechanics (QM/MM) calculations suggest that the reaction is analogous to known endonuclease mechanisms: an enzyme reaction mechanism with two Mg²⁺ ions in the reaction center and a pentacovalent intermediate. The minimum energy pathway profile shows that the phosphodiester bond cleavage step of the reaction is kinetically controlled and not thermodynamically because of a lack of any energy barrier above the accuracy of the energy profile calculation. The role of ions is shown by comparing the results with the reaction mechanisms in the absence of the Mg²⁺ ions where there is a significantly higher reaction barrier than in the presence of the Mg²⁺ ions.Graphical abstract

A family of metal-dependent phosphatases implicated in metabolite damage-control

DUF89 family proteins occur widely in both prokaryotes and eukaryotes, but their functions are unknown. Here we define three DUF89 subfamilies (I, II, and III), with subfamily II being split into stand-alone proteins and proteins fused to pantothenate kinase (PanK). We demonstrated that DUF89 proteins have metal-dependent phosphatase activity against reactive phosphoesters or their damaged forms, notably sugar phosphates (subfamilies II and III), phosphopantetheine and its S-sulfonate or sulfonate (subfamily II-PanK fusions), and nucleotides (subfamily I). Genetic and comparative genomic data strongly associated DUF89 genes with phosphoester metabolism. The crystal structure of the yeast (Saccharomyces cerevisiae) subfamily III protein YMR027W revealed a novel phosphatase active site with fructose 6-phosphate and Mg(2+) bound near conserved signature residues Asp254 and Asn255 that are critical for activity. These findings indicate that DUF89 proteins are previously unrecognized hydrolases whose characteristic in vivo function is to limit potentially harmful buildups of normal or damaged phosphometabolites.

Reversal of DNA damage induced Topoisomerase 2 DNA-protein crosslinks by Tdp2

Mammalian Tyrosyl-DNA phosphodiesterase 2 (Tdp2) reverses Topoisomerase 2 (Top2) DNA-protein crosslinks triggered by Top2 engagement of DNA damage or poisoning by anticancer drugs. Tdp2 deficiencies are linked to neurological disease and cellular sensitivity to Top2 poisons. Herein, we report X-ray crystal structures of ligand-free Tdp2 and Tdp2-DNA complexes with alkylated and abasic DNA that unveil a dynamic Tdp2 active site lid and deep substrate binding trench well-suited for engaging the diverse DNA damage triggers of abortive Top2 reactions. Modeling of a proposed Tdp2 reaction coordinate, combined with mutagenesis and biochemical studies support a single Mg(2+)-ion mechanism assisted by a phosphotyrosyl-arginine cation-π interface. We further identify a Tdp2 active site SNP that ablates Tdp2 Mg(2+) binding and catalytic activity, impairs Tdp2 mediated NHEJ of tyrosine blocked termini, and renders cells sensitive to the anticancer agent etoposide. Collectively, our results provide a structural mechanism for Tdp2 engagement of heterogeneous DNA damage that causes Top2 poisoning, and indicate that evaluation of Tdp2 status may be an important personalized medicine biomarker informing on individual sensitivities to chemotherapeutic Top2 poisons.

Fine tuning of the catalytic activity of colicin E7 nuclease domain by systematic N-terminal mutations

The nuclease domain of colicin E7 (NColE7) promotes the nonspecific cleavage of nucleic acids at its C-terminal HNH motif. Interestingly, the deletion of four N-terminal residues (446-449 NColE7 = KRNK) resulted in complete loss of the enzyme activity. R447A mutation was reported to decrease the nuclease activity, but a detailed analysis of the role of the highly positive and flexible N-terminus is still missing. Here, we present the study of four mutants, with a decreased activity in the following order: NColE7  >> KGNK > KGNG ∼ GGNK > GGNG. At the same time, the folding, the metal-ion, and the DNA-binding affinity were unaffected by the mutations as revealed by linear and circular dichroism spectroscopy, isothermal calorimetric titrations, and gel mobility shift experiments. Semiempirical quantum chemical calculations and molecular dynamics simulations revealed that K446, K449, and/or the N-terminal amino group are able to approach the active centre in the absence of the other positively charged residues. The results suggested a complex role of the N-terminus in the catalytic process that could be exploited in the design of a controlled nuclease.

The cutting edges in DNA repair, licensing, and fidelity: DNA and RNA repair nucleases sculpt DNA to measure twice, cut once

To avoid genome instability, DNA repair nucleases must precisely target the correct damaged substrate before they are licensed to incise. Damage identification is a challenge for all DNA damage response proteins, but especially for nucleases that cut the DNA and necessarily create a cleaved DNA repair intermediate, likely more toxic than the initial damage. How do these enzymes achieve exquisite specificity without specific sequence recognition or, in some cases, without a non-canonical DNA nucleotide? Combined structural, biochemical, and biological analyses of repair nucleases are revealing their molecular tools for damage verification and safeguarding against inadvertent incision. Surprisingly, these enzymes also often act on RNA, which deserves more attention. Here, we review protein-DNA structures for nucleases involved in replication, base excision repair, mismatch repair, double strand break repair (DSBR), and telomere maintenance: apurinic/apyrimidinic endonuclease 1 (APE1), Endonuclease IV (Nfo), tyrosyl DNA phosphodiesterase (TDP2), UV Damage endonuclease (UVDE), very short patch repair endonuclease (Vsr), Endonuclease V (Nfi), Flap endonuclease 1 (FEN1), exonuclease 1 (Exo1), RNase T and Meiotic recombination 11 (Mre11). DNA and RNA structure-sensing nucleases are essential to life with roles in DNA replication, repair, and transcription. Increasingly these enzymes are employed as advanced tools for synthetic biology and as targets for cancer prognosis and interventions. Currently their structural biology is most fully illuminated for DNA repair, which is also essential to life. How DNA repair enzymes maintain genome fidelity is one of the DNA double helix secrets missed by James Watson and Francis Crick, that is only now being illuminated though structural biology and mutational analyses. Structures reveal motifs for repair nucleases and mechanisms whereby these enzymes follow the old carpenter adage: measure twice, cut once. Furthermore, to measure twice these nucleases act as molecular level transformers that typically reshape the DNA and sometimes themselves to achieve extraordinary specificity and efficiency.

Using single-turnover kinetics with osmotic stress to characterize the EcoRV cleavage reaction

Type II restriction endonucleases require metal ions to specifically cleave DNA at canonical sites. Despite the wealth of structural and biochemical information, the number of Mg(2+) ions used for cleavage by EcoRV, in particular, at physiological divalent ion concentrations has not been established. In this work, we employ a single-turnover technique that uses osmotic stress to probe reaction kinetics between an initial specific EcoRV-DNA complex formed in the absence of Mg(2+) and the final cleavage step. With osmotic stress, complex dissociation before cleavage is minimized and the reaction rates are slowed to a convenient time scale of minutes to hours. We find that cleavage occurs by a two-step mechanism that can be characterized by two rate constants. The dependence of these rate constants on Mg(2+) concentration and osmotic pressure gives the number of Mg(2+) ions and water molecules coupled to each kinetic step of the EcoRV cleavage reaction. Each kinetic step is coupled to the binding 1.5-2.5 Mg(2+) ions, the uptake of ∼30 water molecules, and the cleavage of a DNA single strand. We suggest that each kinetic step reflects an independent, rate-limiting conformational change of each monomer of the dimeric enzyme that allows Mg(2+) ion binding. This modified single-turnover protocol has general applicability for metalloenzymes.

Taking into Account the Ion-induced Dipole Interaction in the Nonbonded Model of Ions

Metal ions exist in almost half of the proteins in the protein databank and they serve as structural, electron-transfer and catalytic elements in the metabolic processes of organisms. Molecular Dynamics (MD) simulation is a powerful tool that provides information about biomolecular systems at the atomic level. Coupled with the growth in computing power, algorithms like the Particle Mesh Ewald (PME) method have become the accepted standard when dealing with long-range interactions in MD simulations. The nonbonded model of metal ions consists of an electrostatic plus 12-6 Lennard Jones (LJ) potential and is used largely because of its speed relative to more accurate models. In previous work we found that ideal parameters do not exist that reproduce several experimental properties for M(II) ions simultaneously using the nonbonded model coupled with the PME method due to the underestimation of metal ion-ligand interactions. Via a consideration of the nature of the nonbonded model, we proposed that the observed error largely arises from overlooking charge-induced dipole interactions. The electrostatic plus 12-6 LJ potential model works reasonably well for neutral systems but does struggle with more highly charged systems. In the present work we designed and parameterized a new nonbonded model for metal ions by adding a 1/r4 term to the 12-6 model. We call it the 12-6-4 LJ-type nonbonded model due to its mathematical construction. Parameters were determined for 16 +2 metal ions for the TIP3P, SPC/E and TIP4PEW water models. The final parameters reproduce the experimental hydration free energies (HFE), ion-oxygen distances (IOD) in the first solvation shell and coordination numbers (CN) accurately for the metal ions investigated. Preliminary tests on MgCl2 at different concentrations in aqueous solution and Mg2+--nucleic acid systems show reasonable results suggesting that the present parameters can work in mixed systems. The 12-6-4 LJ-type nonbonded model is readily adopted into standard force fields like AMBER, CHARMM and OPLS-AA with only a modest computational overhead. The new nonbonded model doesn't consider charge-transfer effects explicitly and, hence, may not suitable for the simulation of systems where charge-transfer effects play a decisive role.

mtsslSuite: In silico spin labelling, trilateration and distance-constrained rigid body docking in PyMOL

Nanometer distance measurements based on electron paramagnetic resonance methods in combination with site-directed spin labelling are powerful tools for the structural analysis of macromolecules. The software package mtsslSuite provides scientists with a set of tools for the translation of experimental distance distributions into structural information. The package is based on the previously published mtsslWizard software for in silico spin labelling. The mtsslSuite includes a new version of MtsslWizard that has improved performance and now includes additional types of spin labels. Moreover, it contains applications for the trilateration of paramagnetic centres in biomolecules and for rigid-body docking of subdomains of macromolecular complexes. The mtsslSuite is tested on a number of challenging test cases and its strengths and weaknesses are evaluated.

Rational Design of Particle Mesh Ewald Compatible Lennard-Jones Parameters for +2 Metal Cations in Explicit Solvent

Metal ions play significant roles in biological systems. Accurate molecular dynamics (MD) simulations on these systems require a validated set of parameters. Although there are more detailed ways to model metal ions, the nonbonded model, which employs a 12-6 Lennard-Jones (LJ) term plus an electrostatic potential is still widely used in MD simulations today due to its simple form. However, LJ parameters have limited transferability due to different combining rules, various water models and diverse simulation methods. Recently, simulations employing a Particle Mesh Ewald (PME) treatment for long-range electrostatics have become more and more popular owing to their speed and accuracy. In the present work we have systematically designed LJ parameters for 24 +2 metal (M(II)) cations to reproduce different experimental properties appropriate for the Lorentz-Berthelot combining rules and PME simulations. We began by testing the transferability of currently available M(II) ion LJ parameters. The results showed that there are differences between simulations employing Ewald summation with other simulation methods and that it was necessary to design new parameters specific for PME based simulations. Employing the thermodynamic integration (TI) method and performing periodic boundary MD simulations employing PME, allowed for the systematic investigation of the LJ parameter space. Hydration free energies (HFEs), the ion-oxygen distance in the first solvation shell (IOD) and coordination numbers (CNs) were obtained for various combinations of the parameters of the LJ potential for four widely used water models (TIP3P, SPC/E, TIP4P and TIP4PEW). Results showed that the three simulated properties were highly correlated. Meanwhile, M(II) ions with the same parameters in different water models produce remarkably different HFEs but similar structural properties. It is difficult to reproduce various experimental values simultaneously because the nonbonded model underestimates the interaction between the metal ions and water molecules at short range. Moreover, the extent of underestimation increases successively for the TIP3P, SPC/E, TIP4PEW and TIP4P water models. Nonetheless, we fitted a curve to describe the relationship between ε (the well depth) and radius (Rmin/2) from experimental data on noble gases to facilitate the generation of the best possible compromise models. Hence, by targeting different experimental values, we developed three sets of parameters for M(II) cations for three different water models (TIP3P, SPC/E and TIP4PEW). These parameters we feel represent the best possible compromise that can be achieved using the nonbonded model for the ions in combination with simple water models. From a computational uncertainty analysis we estimate that the uncertainty in our computed HFEs is on the order of ±1kcal/mol. Further improvements will require more advanced non-bonded models likely with inclusion of polarization.

Effect of common metal ions on the rate of degradation of 4-nitrophenol by a laccase-Cu2+ synergistic system

Various metal ions are present in industrial wastewater. In this study, the role of Cu(2+) in pollutant degradation by laccase and the effect of common metal ions such as Na(+), K(+), Ca(2+), Mg(2+), Zn(2+), Hg(2+), Mn(2+), Fe(2+), Co(2+), Ni(2+), Fe(3+), and Al(3+) on the degradation velocity of the laccase-Cu(2+) synergistic system using 4-nitrophenol as the target pollutant were investigated. The results show that the laccase-Cu(2+) system achieved a higher degradation velocity than the system without Cu(2+). The activity of the laccase-Cu(2+) synergistic system was inhibited when a monovalent metal ion (Na(+) or K(+)) was added into the system regardless of the concentration. The addition of relatively low concentrations of any divalent metal ions also resulted in inhibition of the activity. However, increasing concentrations of Ca(2+) or Fe(2+) increased the rate of degradation of 4-nitrophenol by the laccase-Cu(2+) system, whereas increasing concentrations of other divalent metal ions suppressed the system. The inhibition effect of the added trivalent metal ions (Al(3+) or Fe(3+)) was significant during the entire process, indicating that these trivalent metal ions can be a serious obstacle to the activity of the laccase-Cu(2+) synergistic system.

Ca(2+) binding to the ExDxD motif regulates the DNA cleavage specificity of a promiscuous endonuclease

Most of the restriction endonucleases (REases) are dependent on Mg(2+) for DNA cleavage, and in general, Ca(2+) inhibits their activity. R.KpnI, an HNH active site containing ββα-Me finger nuclease, is an exception. In presence of Ca(2+), the enzyme exhibits high-fidelity DNA cleavage and complete suppression of Mg(2+)-induced promiscuous activity. To elucidate the mechanism of unusual Ca(2+)-mediated activity, we generated alanine variants in the putative Ca(2+) binding motif, E(132)xD(134)xD(136), of the enzyme. Mutants showed decreased levels of DNA cleavage in the presence of Ca(2+). We demonstrate that ExDxD residues are involved in Ca(2+) coordination; however, the invariant His of the catalytic HNH motif acts as a general base for nucleophile activation, and the other two active site residues, D148 and Q175, also participate in Ca(2+)-mediated cleavage. Insertion of a 10-amino acid linker to disrupt the spatial organization of the ExDxD and HNH motifs impairs Ca(2+) binding and affects DNA cleavage by the enzyme. Although ExDxD mutant enzymes retained efficient cleavage at the canonical sites in the presence of Mg(2+), the promiscuous activity was greatly reduced, indicating that the carboxyl residues of the acidic triad play an important role in sequence recognition by the enzyme. Thus, the distinct Ca(2+) binding motif that confers site specific cleavage upon Ca(2+) binding is also critical for the promiscuous activity of the Mg(2+)-bound enzyme, revealing its role in metal ion-mediated modulation of DNA cleavage.

Electrophilic assistance to the cleavage of an RNA model phopshodiester via specific and general base-catalyzed mechanisms

Kinetics of transesterification of the RNA model substrate 2-hydroxypropyl 4-nitrophenyl phosphate promoted by Mg(2+) and Ca(2+), the most common biological metals acting as cofactors for nuclease enzymes and ribozymes, as well as by Co(NH(3))(6)(3+), Co(en)(3)(3+), Li(+), and Na(+) cations, often employed as mechanistic probes, was studied in 80% v/v (50 mol %) aqueous DMSO, a medium that allows one to discriminate easily specific base (OH(-)-catalyzed) and general base (buffer-catalyzed) reaction paths. All cations assist the specific base reaction, but only Mg(2+) and Na(+) assist the general base reaction. For Mg(2+)-assisted reactions, the solvent deuterium isotope effects are 1.23 and 0.25 for general base and specific base mechanisms, respectively. Rate constants for Mg(2+)-assisted general base reactions measured with different bases fit the Brønsted correlation with a slope of 0.38, significantly lower than the slope for the unassisted general base reaction (0.77). Transition state binding constants for catalysts in the specific base reaction (K(‡)(OH)) both in aqueous DMSO and pure water correlate with their binding constants to 4-nitrophenyl phosphate dianion (K(NPP)) used as a minimalist transition state model. It was found that K(‡)(OH) ≈ K(NPP) for "protic" catalysts (Co(NH(3))(6)(3+), Co(en)(3)(3+), guanidinium), but K(‡)(OH) ≫ K(NPP) for Mg(2+) and Ca(2+) acting as Lewis acids. It appears from results of this study that Mg(2+) is unique in its ability to assist efficiently the general base-catalyzed transesterification often occurring in active sites of nuclease enzymes and ribozymes.