Mechanistic insight into the catalytic activity of ββα-metallonucleases from computer simulations: Vibrio vulnificus periplasmic nuclease as a test case

The Impact of DFT Functional, Cluster Model Size, and Implicit Solvation on the Structural Description of Single-Metal-Mediated DNA Phosphodiester Bond Cleavage: The Case Study of APE1

Phosphodiester bond hydrolysis in nucleic acids is a ubiquitous reaction that can be facilitated by enzymes called nucleases, which often use metal ions to achieve catalytic function. While a two-metal-mediated pathway has been well established for many enzymes, there is growing support that some enzymes require only one metal for the catalytic step. Using human apurinic/apyrimidinic endonuclease (APE1) as a prototypical example and cluster models, this study clarifies the impact of DFT functional, cluster model size, and implicit solvation on single-metal-mediated phosphodiester bond cleavage and provides insight into how to efficiently model this chemistry. Initially, a model containing 69 atoms built from a high-resolution X-ray crystal structure is used to explore the reaction pathway mapped by a range of DFT functionals and basis sets, which provides support for the use of standard functionals (M06-2X and B3LYP-D3) to study this reaction. Subsequently, systematically increasing the model size to 185 atoms by including additional amino acids and altering residue truncation points highlights that small models containing only a few amino acids or β carbon truncation points introduce model strains and lead to incorrect metal coordination. Indeed, a model that contains all key residues (general base and acid, residues that stabilize the substrate, and amino acids that maintain the metal coordination) is required for an accurate structural depiction of the one-metal-mediated phosphodiester bond hydrolysis by APE1, which results in 185 atoms. The additional inclusion of the broader enzyme environment through continuum solvation models has negligible effects. The insights gained in the present work can be used to direct future computational studies of other one-metal-dependent nucleases to provide a greater understanding of how nature achieves this difficult chemistry.

Characterization of Vibrio cholerae's Extracellular Nuclease Xds

The Gram-negative bacterium Vibrio cholerae encodes two nucleases, Dns and Xds, which play a major role during the human pathogen's lifecycle. Dns and Xds control three-dimensional biofilm formation and bacterial detachment from biofilms via degradation of extracellular DNA and thus contribute to the environmental, inter-epidemic persistence of the pathogen. During intestinal colonization the enzymes help evade the innate immune response, and therefore promote survival by mediating escape from neutrophil extracellular traps. Xds has the additional function of degrading extracellular DNA down to nucleotides, which are an important nutrient source for V. cholerae. Thus, Xds is a key enzyme for survival fitness during distinct stages of the V. cholerae lifecycle and could be a potential therapeutic target. This study provides detailed information about the enzymatic properties of Xds using purified protein in combination with a real time nuclease activity assay. The data define an optimal buffer composition for Xds activity as 50 mM Tris/HCl pH 7, 100 mM NaCl, 10 mM MgCl2, and 20 mM CaCl2. Moreover, maximal activity was observed using substrate DNA with low GC content and ambient temperatures of 20-25°C. In silico analysis and homology modeling predicted an exonuclease domain in the C-terminal part of the protein. Biochemical analyses with truncated variants and point mutants of Xds confirm that the C-terminal region is sufficient for nuclease activity. We also find that residues D787 and H837 within the predicted exonuclease domain are key to formation of the catalytic center.

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.

A comparative study of cold- and warm-adapted Endonucleases A using sequence analyses and molecular dynamics simulations

The psychrophilic and mesophilic endonucleases A (EndA) from Aliivibrio salmonicida (VsEndA) and Vibrio cholera (VcEndA) have been studied experimentally in terms of the biophysical properties related to thermal adaptation. The analyses of their static X-ray structures was no sufficient to rationalize the determinants of their adaptive traits at the molecular level. Thus, we used Molecular Dynamics (MD) simulations to compare the two proteins and unveil their structural and dynamical differences. Our simulations did not show a substantial increase in flexibility in the cold-adapted variant on the nanosecond time scale. The only exception is a more rigid C-terminal region in VcEndA, which is ascribable to a cluster of electrostatic interactions and hydrogen bonds, as also supported by MD simulations of the VsEndA mutant variant where the cluster of interactions was introduced. Moreover, we identified three additional amino acidic substitutions through multiple sequence alignment and the analyses of MD-based protein structure networks. In particular, T120V occurs in the proximity of the catalytic residue H80 and alters the interaction with the residue Y43, which belongs to the second coordination sphere of the Mg2+ ion. This makes T120V an amenable candidate for future experimental mutagenesis.

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.

A new insight into the zinc-dependent DNA-cleavage by the colicin E7 nuclease: a crystallographic and computational study

The crystal structure of a colicin E7 metallonuclease mutant complemented by QM/MM calculations suggests an alternative catalytic mechanism of Zn²⁺-containing HNH nucleases.