Phosphate mono- and diesterase activities of the trinuclear zinc enzyme nuclease P1--insights from quantum chemical calculations

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

Variations in the enzymatic activity of S1-type nucleases results from differences in their active site structures

BACKGROUND

S1-like nucleases are widespread enzymes commonly used in biotechnology and molecular biology. Although it is commonly believed that they are mainly Zn2+-dependent acidic enzymes, we have found that numerous members of this family deviate from this rule. Therefore, in this work, we decided to check how broad is the range of non‑zinc-dependent S1-like nucleases and what is the molecular basis of their activities.

METHODS

S1-like nucleases chosen for analysis were achieved through heterologous expression in appropriate eukaryotic hosts. To characterize nucleases' active-site properties, point mutations were introduced in selected positions. The enzymatic activities of wild-type and mutant nucleases were tested by in-gel nuclease activity assay.

RESULTS

We discovered that S1-like nucleases encoded by non-vascular plants and single-celled protozoa, like their higher plant homologues, exhibit a large variety of catalytic properties. We have shown that these individual properties are determined by specific non-conserved active site residues.

CONCLUSIONS

Our findings demonstrate that mutations that occur during evolution can significantly alter the catalytic properties of S1-like nucleases. As a result, different ions can compete for particular S1-type nucleases' active sites. This phenomenon undermines the existing classification of S1-like nucleases.

GENERAL SIGNIFICANCE

Our findings have numerous implications for applications and understanding the S1-like nucleases' biological functions. For example, new biotechnological applications should take into account their unexpected catalytic properties. Moreover, these results demonstrate that the trinuclear zinc-based model commonly used to characterize the catalytic activities of S1-like nucleases is insufficient to explain the actions of non‑zinc-dependent members of this family.

Expression and function of an S1-type nuclease in the digestive fluid of a sundew, Drosera adelae

BACKGROUND AND AIMS

Carnivorous plants trap and digest insects and similar-sized animals. Many studies have examined enzymes in the digestive fluids of these plants and have gradually unveiled the origins and gene expression of these enzymes. However, only a few attempts have been made at characterization of nucleases. This study aimed to reveal gene expression and the structural, functional and evolutionary characteristics of an S1-type nuclease (DAN1) in the digestive fluid of an Australian sundew, Drosera adelae, whose trap organ shows unique gene expression and related epigenetic regulation.

METHODS

Organ-specificity in Dan1 expression was examined using glandular tentacles, laminas, roots and inflorescences, and real-time PCR. The methylation status of the Dan1 promoter in each organ was clarified by bisulphite sequencing. The structural characteristics of DAN1 were studied by a comparison of primary structures of S1-type nucleases of three carnivorous and seven non-carnivorous plants. DAN1 was prepared using a cell-free protein synthesis system. Requirements for metal ions, optimum pH and temperature, and substrate preference were examined using conventional methods.

KEY RESULTS

Dan1 is exclusively expressed in the glandular tentacles and its promoter is almost completely unmethylated in all organs. This is in contrast to the S-like RNase gene da-I of Dr. adelae, which shows similar organ-specific expression, but is controlled by a promoter that is specifically unmethylated in the glandular tentacles. Comparison of amino acid sequences of S1-type nucleases identifies seven and three positions where amino acid residues are conserved only among the carnivorous plants and only among the non-carnivorous plants, respectively. DAN1 prefers a substrate RNA over DNA in the presence of Zn2+, Mn2+ or Ca2+ at an optimum pH of 4.0.

CONCLUSIONS

Uptake of phosphates from prey is suggested to be the main function of DAN1, which is very different from the known functions of S1-type nucleases. Evolution has modified the structure and expression of Dan1 to specifically function in the digestive fluid.

The metal dependence of single-metal mediated phosphodiester bond cleavage: a QM/MM study of a multifaceted human enzyme

Nucleases catalyze the cleavage of phosphodiester bonds in nucleic acids using a range of metal cofactors. Although it is well accepted that many nucleases rely on two metal ions, the one-metal mediated pathway is debated. Furthermore, one-metal mediated nucleases maintain activity in the presence of many different metals, but the underlying reasons for this broad metal specificity are unknown. The human apurinic/apyrimidinic endonuclease (APE1), which plays a key role in DNA repair, transcription regulation, and gene expression, is a prototypical example of a one-metal dependent nuclease. Although Mg²⁺ is the native metal cofactor, APE1 remains catalytically active in the presence of several metals, with the rate decreasing as Mg²⁺ > Mn²⁺ > Ni²⁺ > Zn²⁺, while Ca²⁺ completely abolished the activity. The present work uses quantum mechanics-molecular mechanics techniques to map APE1-facilitated phosphodiester bond hydrolysis in the presence of these metals. The structural differences in stationary points along the reaction pathway shed light on the interplay between several factors that allow APE1 to remain catalytically active for various metals, with the trend in the barrier heights correlating with the experimentally reported APE1 catalytic activity. In contrast, Ca²⁺ significantly changes the metal coordination and active site geometry, and thus completely inhibits catalysis. Our work thereby provides support for the controversial single-metal mediated phosphodiester bond cleavage and clarifies uncertainties regarding the role of the metal and metal identity in this important reaction. This information is key for future medicinal and biotechnological applications including disease diagnosis and treatment, and protein engineering.

Frontiers of metal-coordinating drug design

Introduction: The occurrence of metal ions in biomolecules is required to exert vital cellular functions. Metal-containing biomolecules can be modulated by small-molecule inhibitors targeting their metal-moiety. As well, the discovery of cisplatin ushered the rational discovery of metal-containing-drugs. The use of both drug types exploiting metal-ligand interactions is well established to treat distinct pathologies. Therefore, characterizing and leveraging metal-coordinating drugs is a pivotal, yet challenging, part of medicinal chemistry.Area covered: Atomic-level simulations are increasingly employed to overcome the challenges met by traditional drug-discovery approaches and to complement wet-lab experiments in elucidating the mechanisms of drugs' action. Multiscale simulations, allow deciphering the mechanism of metal-binding inhibitors and metallo-containing-drugs, enabling a reliable description of metal-complexes in their biological environment. In this compendium, the authors review selected applications exploiting the metal-ligand interactions by focusing on understanding the mechanism and design of (i) inhibitors targeting iron and zinc-enzymes, and (ii) ruthenium and gold-based anticancer agents targeting the nucleosome and aquaporin protein, respectively.Expert opinion: The showcased applications exemplify the current role and the potential of atomic-level simulations and reveal how their synergic use with experiments can contribute to uncover fundamental mechanistic facets and exploit metal-ligand interactions in medicinal chemistry.

A Transient and Flexible Cation-π Interaction Promotes Hydrolysis of Nucleic Acids in DNA and RNA Nucleases

Metal-dependent DNA and RNA nucleases are enzymes that cleave nucleic acids with great efficiency and precision. These enzyme-mediated hydrolytic reactions are fundamental for the replication, repair, and storage of genetic information within the cell. Here, extensive classical and quantum-based free-energy molecular simulations show that a cation-π interaction is transiently formed in situ at the metal core of Bacteriophage-λ Exonuclease (Exo-λ), during catalysis. This noncovalent interaction (Lys131-Tyr154) triggers nucleophile activation for nucleotide excision. Then, our simulations also show the oscillatory dynamics and swinging of the newly formed cation-π dyad, whose conformational change may favor proton release from the cationic Lys131 to the bulk solution, thus restoring the precatalytic protonation state in Exo-λ. Altogether, we report on the novel mechanistic character of cation-π interactions for catalysis. Structural and bioinformatic analyses support that flexible orientation and transient formation of mobile cation-π interactions may represent a common catalytic strategy to promote nucleic acid hydrolysis in DNA and RNA nucleases.

Dimerisation induced formation of the active site and the identification of three metal sites in EAL-phosphodiesterases

The bacterial second messenger cyclic di-3',5'-guanosine monophosphate (c-di-GMP) is a key regulator of bacterial motility and virulence. As high levels of c-di-GMP are associated with the biofilm lifestyle, c-di-GMP hydrolysing phosphodiesterases (PDEs) have been identified as key targets to aid development of novel strategies to treat chronic infection by exploiting biofilm dispersal. We have studied the EAL signature motif-containing phosphodiesterase domains from the Pseudomonas aeruginosa proteins PA3825 (PA3825^EAL) and PA1727 (MucR^EAL). Different dimerisation interfaces allow us to identify interface independent principles of enzyme regulation. Unlike previously characterised two-metal binding EAL-phosphodiesterases, PA3825^EAL in complex with pGpG provides a model for a third metal site. The third metal is positioned to stabilise the negative charge of the 5'-phosphate, and thus three metals could be required for catalysis in analogy to other nucleases. This newly uncovered variation in metal coordination may provide a further level of bacterial PDE regulation.

Structural and Catalytic Properties of S1 Nuclease from Aspergillus oryzae Responsible for Substrate Recognition, Cleavage, Non-Specificity, and Inhibition

The single-strand-specific S1 nuclease from Aspergillus oryzae is an archetypal enzyme of the S1-P1 family of nucleases with a widespread use for biochemical analyses of nucleic acids. We present the first X-ray structure of this nuclease along with a thorough analysis of the reaction and inhibition mechanisms and of its properties responsible for identification and binding of ligands. Seven structures of S1 nuclease, six of which are complexes with products and inhibitors, and characterization of catalytic properties of a wild type and mutants reveal unknown attributes of the S1-P1 family. The active site can bind phosphate, nucleosides, and nucleotides in several distinguished ways. The nucleoside binding site accepts bases in two binding modes-shallow and deep. It can also undergo remodeling and so adapt to different ligands. The amino acid residue Asp65 is critical for activity while Asn154 secures interaction with the sugar moiety, and Lys68 is involved in interactions with the phosphate and sugar moieties of ligands. An additional nucleobase binding site was identified on the surface, which explains the absence of the Tyr site known from P1 nuclease. For the first time ternary complexes with ligands enable modeling of ssDNA binding in the active site cleft. Interpretation of the results in the context of the whole S1-P1 nuclease family significantly broadens our knowledge regarding ligand interaction modes and the strategies of adjustment of the enzyme surface and binding sites to achieve particular specificity.

μ3-Oxo stabilized by three metal cations is a sufficient nucleophile for enzymatic hydrolysis of phosphate monoesters

Diverse species have previously been proposed to be effective nucleophiles in the enzymatic hydrolysis of phosphate esters. A novel penta-metal cluster (two Fe(3+) and three Ca(2+)) was recently discovered in the active site of PhoX alkaline phosphatase, with the revelation of the architecture of μ3-oxo bridging one Ca(2+) and two antiferromagnetically coupled Fe(3+). In this work, using density functional theory calculations, the μ3-oxo stabilized by three cations has been demonstrated to be a new type of effective nucleophile. The calculations give strong support to the "ping-pong" mechanism involving the nucleophilic attack of the μ3-oxo on the substrate phosphor and the subsequent hydrolysis of the covalent phospho-enzyme intermediate. A base mechanism with the μ3-oxo acting as a general base to activate an additional water molecule has further been demonstrated to be inaccessible. The results advance the understanding of the enzymatic hydrolysis of phosphate esters and may give important inspiration for the exploration of multinuclear biomimetic catalysts.

Theoretically exploring the key role of the Lys412 residue in the conversion of N2O to N2by nitrous oxide reductase from Achromobacter cycloclastes

The active Cu_Zcluster of NOR provides strong back-donation to coordinated N₂O and activates the O atom of the N₂O group facilitating H-bonding and protonationviathe Lys412 residue.

Visualizing phosphodiester-bond hydrolysis by an endonuclease

The enzymatic hydrolysis of DNA phosphodiester bonds has been widely studied, but the chemical reaction has not yet been observed. Here we follow the generation of a DNA double-strand break (DSB) by the Desulfurococcus mobilis homing endonuclease I-DmoI, trapping sequential stages of a two-metal-ion cleavage mechanism. We captured intermediates of the different catalytic steps, and this allowed us to watch the reaction by 'freezing' multiple states. We observed the successive entry of two metals involved in the reaction and the arrival of a third cation in a central position of the active site. This third metal ion has a crucial role, triggering the consecutive hydrolysis of the targeted phosphodiester bonds in the DNA strands and leaving its position once the DSB is generated. The multiple structures show the orchestrated conformational changes in the protein residues, nucleotides and metals during catalysis.

Phosphate monoester hydrolysis by trinuclear alkaline phosphatase; DFT study of transition States and reaction mechanism

Alkaline phosphatase (AP) is a trinuclear metalloenzyme that catalyzes the hydrolysis of a broad range of phosphate monoesters to form inorganic phosphate and alcohol (or phenol). In this paper, by using density functional theory with a model based on a crystal structure, the AP-catalyzed hydrolysis of phosphate monoesters is investigated by calculating two substrates, that is, methyl and p-nitrophenyl phosphates, which represent alkyl and aryl phosphates, respectively. The calculations confirm that the AP reaction employs a "ping-pong" mechanism involving two chemical displacement steps, that is, the displacement of the substrate leaving group by a Ser102 alkoxide and the hydrolysis of the phosphoseryl intermediate by a Zn2-bound hydroxide. Both displacement steps proceed via a concerted associative pathway no matter which substrate is used. Other mechanistic aspects are also studied. Comparison of our calculations with linear free energy relationships experiments shows good agreement.

A dinuclear zinc(II) complex of a new unsymmetric ligand with an N(5)O(2) donor set: a structural and functional model for the active site of zinc phosphoesterases

The dinuclear complex [Zn(2)(DPCPMP)(pivalate)](ClO4), where DPCPMP is the new unsymmetrical ligand [2-(N-(3-((bis((pyridin-2-yl)methyl)amino)methyl)-2-hydroxy-5-methylbenzyl)-N-((pyridin-2-yl)methyl)amino)acetic acid], has been synthesized and characterized. The complex is a functional model for zinc phosphoesterases with dinuclear active sites. The hydrolytic efficacy of the complex has been investigated using bis-(2,4-dinitrophenyl)phosphate (BDNPP), a DNA analog, as substrate. Speciation studies using potentiometric titrations have been performed for both the ligand and the corresponding dizinc complex to elucidate the formation of the active hydrolysis catalyst; they reveals that the dinuclear zinc(II) complexes, [Zn(2)(DPCPMP)](2+) and [Zn(2)(DPCPMP)(OH)](+) predominate the solution above pH4. The relatively high pK(a) of 8.38 for water deprotonation suggests that a terminal hydroxide complex is formed. Kinetic investigations of BDNPP hydrolysis over the pH range 5.5-11.0 and with varying metal to ligand ratio (metal salt:ligand=0.5:1 to 3:1) have been performed. Variable temperature studies gave the activation parameters ΔH(‡)=95.6kJmol(-1), ΔS(‡)=-44.8Jmol(-1)K(-1), and ΔG(‡)=108.0 kJmol(-1). The cumulative results indicate the hydroxido-bridged dinuclear Zn(II) complex [Zn(2)(DPCPMP)(μ-OH)](+) as the effective catalyst. The mechanism of hydrolysis has been probed by computational modeling using density functional theory (DFT). Calculations show that the reaction goes through one concerted step (S(N)2 type) in which the bridging hydroxide in the transition state becomes terminal and performs a nucleophilic attack on the BDNPP phosphorus; the leaving group dissociates simultaneously in an overall inner sphere type activation. The calculated free energy barrier is in good agreement with the experimentally determined activation parameters.

The plant s1-like nuclease family has evolved a highly diverse range of catalytic capabilities

Plant S1-like nucleases, often referred to as nuclease I enzymes, are the main class of enzymes involved in nucleic acid degradation during plant programmed cell death. The catalytically active site of these enzymes shows a significant similarity to the well-described P1 nuclease from Penicillium citrinum. Previously published studies reported that plant S1-like nucleases possess catalytic activities similar to their fungal orthologs, i.e. they hydrolyze single-stranded DNA and RNA, and less efficiently double-stranded DNA, in the presence of zinc ions. Here we describe a comprehensive study of the nucleolytic activities of all Arabidopsis S1-like paralogs. Our results revealed that different members of this family are characterized by a surprisingly large variety of catalytic properties. We found that, in addition to Zn(2+)-dependent enzymes, this family also comprises nucleases activated by Ca(2+) and Mn(2+), which implies that the apparently well-known S1 nuclease active site in plant nucleases is able to cooperate with different activatory ions. Moreover, particular members of this class differ in their optimum pH value and substrate specificity. These results shed new light on the widely accepted classification of plant nucleases which is based on the assumption that the catalytic requirements of plant nucleases reflect their phylogenetic origin. Our results imply the need to redefine the understanding of the term 'nuclease I'. Analysis of the phylogenetic relationships between S1-like enzymes shows that plant representatives of this family evolve toward an increase in catalytic diversity. The importance of this process for the biological functions of plant S1-type enzymes is discussed.