Theoretical Study of Phosphodiester Hydrolysis and Transesterification Catalyzed by an Unsymmetric Biomimetic Dizinc Complex

Unveiling the Cleavage Mechanism of an RNA Model Compound on the whole pH Scale: Computations Meet Experiments in the Determination of Reaction Rates

The knowledge of the mechanism of reactions occurring in solution is a primary research line both in the context of theoretical-computational chemistry and in the field of organic and bio-organic chemistry. Given the importance of the hydrolysis of nucleic acids in life-related phenomena, here we present a combined experimental and computational study on the cleavage of an RNA model compound. This phosphodiester features a cleavage rate strictly dependent on the pH with three different dependence domains. Such experimental evidence, highlighted by an in-depth kinetic investigation, unequivocally suggests a change in the reaction mechanism along the pH scale. In order to interpret the data and to explain the experimental behavior, we have applied a theoretical-computational procedure, involving a hybrid quantum/classical approach, able to model chemical reactions in complex environments, i. e. in solution. This study turns out to quantitatively reproduce the experimental data with accuracy and, in addition, provides useful mechanistic insight into the transesterification process of the investigated compound. The study indicates that the cleavage can occur through anA N D N ${A_N D_N }$ , anA N + D N ${A_N + D_N }$ , and aD N A N ${D_N A_N }$ mechanism depending on the pH values.

Distinct chemical factors in hydrolytic reactions catalyzed by metalloenzymes and metal complexes

The selective hydrolysis of the extremely stable phosphoester, peptide and ester bonds of molecules by bio-inspired metal-based catalysts (metallohydrolases) is required in a wide range of biological, biotechnological and industrial applications. Despite the impressive advances made in the field, the ultimate goal of designing efficient enzyme mimics for these reactions is still elusive. Its realization will require a deeper understanding of the diverse chemical factors that influence the activities of both natural and synthetic catalysts. They include catalyst-substrate complexation, non-covalent interactions and the electronic nature of the metal ion, ligand environment and nucleophile. Based on our computational studies, their roles are discussed for several mono- and binuclear metallohydrolases and their synthetic analogues. Hydrolysis by natural metallohydrolases is found to be promoted by a ligand environment with low basicity, a metal bound water and a heterobinuclear metal center (in binuclear enzymes). Additionally, peptide and phosphoester hydrolysis is dominated by two competing effects, i.e. nucleophilicity and Lewis acid activation, respectively. In synthetic analogues, hydrolysis is facilitated by the inclusion of a second metal center, hydrophobic effects, a biological metal (Zn, Cu and Co) and a terminal hydroxyl nucleophile. Due to the absence of the protein environment, hydrolysis by these small molecules is exclusively influenced by nucleophile activation. The results gleaned from these studies will enhance the understanding of fundamental principles of multiple hydrolytic reactions. They will also advance the development of computational methods as a predictive tool to design more efficient catalysts for hydrolysis, Diels-Alder reaction, Michael addition, epoxide opening and aldol condensation.

Transesterification of RNA model induced by novel dinuclear copper (II) complexes with bis-tridentate imidazole derivatives

Two novel bis-tridentate imidazole derivatives were conveniently synthesized using a 'one-pot' method. Their dinuclear (Cu₂L¹Cl₄, Cu₂L²Cl₄) and mononuclear (CuL¹Cl₂, CuL²Cl₂∙H₂O) copper (II) complexes were synthesized to comparably evaluate their reactivities in the hydrolytic cleavage of 2-hydroxypropyl p-nitrophenyl phosphate (HPNP) as a classic RNA model. Single crystals of Cu₂L¹Cl₄ and Cu₂L²Cl₄ indicate that both of them are centrosymmetric, and each central copper ion is penta-coordinated. Regarding the transesterification of HPNP, both of dinuclear ones exhibited excess one order of magnitude rate enhancement in contrast with auto-hydrolysis reaction. Under comparable conditions, dinuclear complexes displayed no more than twofold increase in activity over their mononuclear analogues, which verifies the lack of binuclear cooperation effect due to long Cu-to-Cu space.

Intra- and Intermolecular Cooperativity in the Catalytic Activity of Phosphodiester Cleavage by Self-Assembled Systems Based on Guanidinylated Calix[4]arenes

The calix[4]arene scaffold, blocked in the cone conformation through alkylation with long alkyl chains, and decorated at the upper rim with four guanidine or arginine units, effectively catalyzes the cleavage of the phosphodiester bond of DNA and RNA model compounds in water. An exhaustive kinetic investigation unequivocally points to the existence of spontaneous aggregation phenomena, driven by hydrophobic effect, occurring at different critical concentrations that depend on the identity of the compound. A pronounced superiority of the assembled structures compared with the monomers in solution was observed. Moreover, the catalytically active units, clustered on the macrocyclic tetrafunctional scaffold, were proved to efficiently cooperate in the catalytic mechanism and result in improved reaction rates compared to those of the monofunctional model compounds. The kinetic analysis is also integrated and corroborated with further experiments based on fluorescence spectroscopy and light scattering. The advantage of the supramolecular assemblies based on tetrafunctional calixarenes leads to believe that the active units can cooperate not only intramolecularly but also intermolecularly. The molecules in the aggregates can probably mold, flex and rearrange but, at the same time, keep an ordered structure that favors phosphodiester bond cleavage. This dynamic preorganization can allow the catalytic units to reach a better fitting with the substrates and perform a superior catalytic activity.

Nuclease-like metalloscissors: Biomimetic candidates for cancer and bacterial and viral infections therapy

Despite the extensive and rapid discovery of modern drugs for treatment of cancer, microbial infections, and viral illnesses; these diseases are still among major global health concerns. To take inspiration from natural nucleases and also the therapeutic potential of metallopeptide antibiotics such as the bleomycin family, artificial metallonucleases with the ability of promoting DNA/RNA cleavage and eventually affecting cellular biological processes can be introduced as a new class of therapeutic candidates. Metal complexes can be considered as one of the main categories of artificial metalloscissors, which can prompt nucleic acid strand scission. Accordingly, biologists, inorganic chemists, and medicinal inorganic chemists worldwide have been designing, synthesizing and evaluating the biological properties of metal complexes as artificial metalloscissors. In this review, we try to highlight the recent studies conducted on the nuclease-like metalloscissors and their potential therapeutic applications. Under the light of the concurrent Covid-19 pandemic, the human need for new therapeutics was highlighted much more than ever before. The nuclease-like metalloscissors with the potential of RNA cleavage of invading viral pathogens hence deserve prime attention.

Binuclear metal complexes with a novel hexadentate imidazole derivative for the cleavage of phosphate diesters and biomolecules: distinguishable mechanisms

Oxidative cleavage of phosphate diesters (HPNP, BNPP) is highly faster than the hydrolytic one by binuclear metal complexes with novel imidazole derivative, producing a non-lactone phosphate monoester due to the direct attack of free radicals.

Electronic Effect on Bimetallic Catalysts: Cleavage of Phosphodiester Mediated by Fe(III)-Zn(II) Purple Acid Phosphatase Mimics

The bimetallic system is an important strategy for the catalytic hydrolysis of phosphodiester. The purple acid phosphatase (PAPs) enzyme is a typical bimetallic catalyst in this field. Mechanistic details for the hydrolysis cleavage of the DNA dinucleotide analogue BNPP^- (BNPP^- = bis(p-nitrophenyl) phosphate) by hetero-binuclear [Fe^III(μ-OH)Zn^IIL]²⁺ complexes (L = 2-[N-bis(2-pyridylmethyl)-aminomethyl]-4-methyl-6-[N'-(2-pyridylmethyl)(2-hydroxybenzyl) aminomethyl] phenol) were investigated using density functional theory calculations. The catalysts with single-bridged hydroxyl and double-bridged hydroxyl groups were compared. The calculation results show that the doubly hydroxide-bridged complex could better bind to substrates. For the BNPP^- hydrolysis, the doubly hydroxide-bridged reactant isomerizes into a single hydroxide-bridged complex, and then the attack is initiated by the hydroxyl group on the iron center. In addition, the catalyst with the electron-donating group (Me) was determined to take precedence over electron-withdrawing groups (Br and NO₂ groups) in the hydrolysis reaction. This is because the substituents affect the high-lying occupied molecular orbitals, tuning the Lewis acidity of iron and pK_a values of the metal-bonded water. These factors influence the hydroxyl nucleophilicity, leading to changes in catalytic activity. To further examine substituent effects, the occupied orbital energies were calculated with several different substituent groups (-CF₃, -OMe, -OH, -NH₂, and -N(Me)₂). It was found that the HOMO or HOMO-1 energy decreases with the increase of the σ_p value. Further, the catalyst activity of the [Fe^III(μ-OH)Zn^IIL]²⁺ complexes was found to be mainly affected by the phenolate ligand (B) coordinated to the iron and zinc centers. These fundamental aspects of the hydrolysis reactions of BNPP^- catalyzed by [Fe^III(μ-OH)Zn^IIL]²⁺ complexes should contribute to improved understanding of the mechanism and to catalyst design involving hetero-binuclear metals complexes.

Executing a Series of Zinc(II) Complexes of Homologous Schiff Base Ligands for a Comparative Analysis on Hydrolytic, Antioxidant, and Antibacterial Activities

Six zinc(II) complexes, namely, [Zn(HL¹H)Cl₂] (1), [Zn(HL¹H)Br₂] (2), [Zn₂(HL¹H)₂(OH)I₂]·I (3), [Zn(HL²)Cl] (4), [Zn₂(HL²)Br₃] (5), and [Zn(HL²)I] (6) have been manufactured by using two homologous Schiff base ligands H₂L¹ and H₂L² for the purpose of perlustrating their phosphatase-like activity, antioxidant activity, and antibacterial activity. Complexes 1, 2, 4, and 5 have been reported earlier by us, whereas complexes 3 and 6 have been synthesized and structurally characterized by regular physicochemical methods The hydrolytic property of the six complexes has been evaluated by checking the hydrolysis of the P-O bond of a widely used substrate, namely, disodium salt of (para-nitrophenyl)phosphate (PNPP) in 97.5% (v/v) mixture of N,N-dimethylformamide and water (DMF-water). Complexes 2-5 have profound efficiency toward hydrolysis of phosphate ester bonds, and complexes 1 and 6 were noted to be inactive toward hydrolysis. Complex 3 displayed the highest efficacy among the six complexes. Additionally, antioxidant and antibacterial activities of the complexes were studied thoroughly. A detailed study of their antioxidant property revealed that complex 3 manifested superior radical scavenging activity, thus exhibiting the highest antioxidant property. The antibacterial activity was tested using four investigating bacteria, specifically Listeria monocytogenes ATCC19111, Staphylococcus aureus ATCC 700699, Salmonella typhimurium ATCC 23564, and Escherichia coli ATCC 25922 by determining minimum inhibitory concentration (MIC) values using the microdilution method. Here as well, complex 3 exhibited the highest activity to both Gram positive and Gram negative bacteria. The chemistry behind these experimental findings has been manifested by shedding light upon the structural features of the complexes. The suitable choice of ligand H₂L¹ where one methylene group is less than its homologous ligand and metal precursor (ZnI₂) imparts a unique hydroxo-bridged molecular geometry and 2D hydrogen bonding network which in turn probably enhances the hydrolytic and biological activities of complex 3.

Effects of the Metal Ion on the Mechanism of Phosphodiester Hydrolysis Catalyzed by Metal-Cyclen Complexes

In this study, mechanisms of phosphodiester hydrolysis catalyzed by six di- and tetravalent metal-cyclen (M-C) complexes (Zn-C, Cu-C, Co-C, Ce-C, Zr-C and Ti-C) have been investigated using DFT calculations. The activities of these complexes were studied using three distinct mechanisms: (1) direct attack ( DA ), (2) catalyst-assisted ( CA ), and (3) water-assisted ( WA ). All divalent metal complexes (Zn-C, Cu-C and Co-C) coordinated to the BNPP substrate in a monodentate fashion and activated its scissile phosphoester bond. However, all tetravalent metal complexes (Ce-C, Zr-C, and Ti-C) interacted with BNPP in a bidentate manner and strengthened this bond. The DA mechanism was energetically the most feasible for all divalent M-C complexes, while the WA mechanism was favored by the tetravalent complexes, except Ce-C. The divalent complexes were found to be more reactive than their tetravalent counterparts. Zn-C catalyzed the hydrolysis with the lowest barrier among all M-C complexes, while Ti-C was the most reactive tetravalent complex. The activities of Ce-C and Zr-C, except Ti-C, were improved with an increase in the coordination number of the metal ion. The structural and mechanistic information provided in this study will be very helpful in the development of more efficient metal complexes for this critical reaction.

Mechanistic Studies of Homo- and Heterodinuclear Zinc Phosphoesterase Mimics: What Has Been Learned?

Phosphoesterases hydrolyze the phosphorus oxygen bond of phosphomono-, di- or triesters and are involved in various important biological processes. Carboxylate and/or hydroxido-bridged dizinc(II) sites are a widespread structural motif in this enzyme class. Much effort has been invested to unravel the mechanistic features that provide the enormous rate accelerations observed for enzymatic phosphate ester hydrolysis and much has been learned by using simple low-molecular-weight model systems for the biological dizinc(II) sites. This review summarizes the knowledge and mechanistic understanding of phosphoesterases that has been gained from biomimetic dizinc(II) complexes, showing the power as well as the limitations of model studies.

A calix[4]arene with acylguanidine units as an efficient catalyst for phosphodiester bond cleavage in RNA and DNA model compounds

Conjugated carbonyl units in a calixarene scaffold provide the right amount of flexibility for catalysis with a minimum entropic cost.

Mapping the working route of phosphate monoester hydrolysis catalyzed by copper based models with special emphasis on the role of oxoanions by experimental and theoretical studies

The mechanistic pathway of phosphate-ester bond hydrolysis with special emphasis on the role of oxoanions was explored by experimental and theoretical study.

Mixed Explicit-Implicit Solvation Approach for Modeling of Alkane Complexation in Water-Soluble Self-Assembled Capsules

The host-guest binding properties of a water-soluble resorcinarene-based cavitand are examined using density functional theory methodology. Experimentally, the cavitand has been observed to self-assemble in aqueous solution into both 1:1 and 2:1 host/guest complexes with hydrophobic guests such as n-alkanes. For n-decane, equilibrium was observed between the 1:1 and 2:1 complexes, while 1:1 complexes are formed with shorter n-alkanes and 2:1 complexes are formed with longer ones. These findings are used to assess the standard quantum chemical methodology. It is first shown that a rather advanced computational protocol (B3LYP-D3(BJ)/6-311+G(2d,2p) with COSMO-RS and quasi-rigid-rotor-harmonic-oscillator) gives very large errors. Systematic examination of the various elements of the methodology shows that the error stems from the implicit solvation model. A mixed explicit-implicit solvation protocol is developed that involves a parametrization of the hydration free energy of water such that water cluster formation in water is predicted to be thermoneutral. This new approach is demonstrated to lead to a major improvement in the calculated binding free energies of n-alkanes, reproducing very well the 1:1 versus 2:1 host/guest binding trends.

Catalytic degradation of organophosphorous nerve agent simulants by polymer beads@graphene oxide with organophosphorus hydrolase-like activity based on rational design of functional bimetallic nuclear ligand

The degradation of organophosphorous nerve agents is of primary concern due to the severe toxicity of these agents. Based on the active center of organophosphorus hydrolase (OPH), a bimetallic nuclear ligand, (5-vinyl-1,3-phenylene)bis(di(1H-imidazol-2-yl) methanol) (VPIM), was designed and synthesized, which contains four imidazole groups to mimic the four histidines at OPH active center. By grafting VPIM on graphene oxide (GO) surface via polymerization, the VPIM-polymer beads@GO was produced. The obtained OPH mimics has an impressive activity in dephosphorylation reactions (turnover frequency (TOF) towards paraoxon: 2.3 s^-1). The synergistic catalytic effect of the bimetallic Zn²⁺ nuclear center and carboxyl groups on surface of GO possibly contributes to the high hydrolysis on organophosphate substrate. Thus, a biomimetic catalyst for efficient degradation of some organophosphorous nerve agent simulants, such as paraoxon and chlorpyrifos, was prepared by constructing catalytic active sites. The proposed mechanism and general synthetic strategy open new avenues for the engineering of functional GOs for biomimetic catalysts.

Higher Flexibility of Glu-172 Explains the Unusual Stereospecificity of Glyoxalase I

Despite many studies during the latest two decades, the reason for the unusual stereospecificity of glyoxalase I (GlxI) is still unknown. This metalloenzyme converts both enantiomers of its natural substrate to only one enantiomer of its product. In addition, GlxI catalyzes reactions involving some substrate and product analogues with a stereospecificity similar to that of its natural substrate reaction. For example, the enzyme exchanges the pro- S, but not the pro- R, hydroxymethyl proton of glutathiohydroxyacetone (HOC-SG) with a deuterium from D₂O. To find some clues to the unusual stereospecificity of GlxI, we have studied the stereospecific proton exchange of the hydroxymethyl proton of HOC-SG by this enzyme. We employed density functional theory and molecular dynamics (MD) simulations to study the proton exchange mechanism and origin of the stereospecificity. The results show that a rigid cluster model with the same flexibility for the two active-site glutamate residues cannot explain the unusual stereospecificity of GlxI. However, using a cluster model with full flexibility of Glu-172 or a larger model with the entire glutamates, extending the backbone into the neighboring residues, the results showed that there is no way for HOC-SG to exchange its protons if the alcoholic proton is directed toward Glu-99. However, if the hydroxymethyl proton instead is directed toward the more flexible Glu-172, we find a catalytic reaction mechanism for the exchange of the H_S proton by a deuterium, in accordance with experimental findings. Thus, our results indicate that the special stereospecificity of GlxI is caused by the more flexible environment of Glu-172 in comparison to that of Glu-99. This higher flexibility of Glu-172 is also confirmed by MD simulations. We propose a reaction mechanism for the stereospecific proton exchange of the hydroxymethyl proton of HOC-SG by GlxI with an overall energy barrier of 15 kcal/mol.

Hydrolysis of Phosphate Esters Catalyzed by Inorganic Iron Oxide Nanoparticles Acting as Biocatalysts

Phosphorus ester hydrolysis is one of the key chemical processes in biological systems, including signaling, free-energy transaction, protein synthesis, and maintaining the integrity of genetic material. Hydrolysis of this otherwise kinetically stable phosphoester and/or phosphoanhydride bond is induced by enzymes such as purple acid phosphatase. Here, I report that, as in previously reported aged inorganic iron ion solutions, the iron oxide nanoparticles in the solution, which are trapped in a dialysis membrane tube filled with the various iron oxides, significantly promote the hydrolysis of the various phosphate esters, including the inorganic polyphosphates, with enzyme-like kinetics. This observation, along with those of recent studies of iron oxide, vanadium pentoxide, and molybdenum trioxide nanoparticles that behave as mimics of peroxidase, bromoperoxidase, and sulfite oxidase, respectively, indicates that the oxo-metal bond in the oxide nanoparticles is critical for the function of these corresponding natural metalloproteins. These inorganic biocatalysts challenge the traditional concept of replicator-first scenarios and support the metabolism-first hypothesis. As biocatalysts, these inorganic nanoparticles with enzyme-like activity may work in natural terrestrial environments and likely were at work in early Earth environments as well. They may have played an important role in the C, H, O, S, and P metabolic pathway with regard to the emergence and early evolution of life. Key Words: Enzyme-Hydrolysis-Iron oxide-Nanoparticles-Origin of life-Phosphate ester. Astrobiology 18, 294-310.

A phenoxo-bridged dicopper(ii) complex as a model for phosphatase activity: mechanistic insights from a combined experimental and computational study

A μ-phenoxo-bis(μ₂-1,3-acetato)-bridged dicopper(ii) complex [Cu(L¹)(μ-O₂CMe)₂][NO₃] (1) has been synthesized from the perspective of modeling phosphodiesterase activity. Structural characterization was done initially with 1·3Et₂O (vapour diffusion of Et₂O into MeOH solution of 1; poor crystal quality) and finally with its perchlorate salt [Cu(L¹)(μ-O₂CMe)₂][ClO₄]·1.375MeCN·0.25H₂O, crystallized from vapour diffusion of n-pentane into a MeCN-MeOH mixture (comparatively better crystal quality). An asymmetric unit of such a crystal contains two independent molecules of compositions [Cu(L¹)(μ-O₂CMe)₂][ClO₄] and [Cu(L¹)(μ-O₂CMe)₂(MeCN)][ClO₄] (coordinated MeCN with 0.75 occupancy), and two molecules of MeCN and H₂O (each H₂O molecule with 0.25 occupancy) as the solvent of crystallization. These two cations, each having five-coordinate (μ-phenoxo)bis(μ-acetato)-bridged Cu^II ions, differ by only the coordination environment of only one Cu^II ion, which has a weakly coordinated acetonitrile molecule in its sixth position. Temperature-dependent magnetic studies on 1 reveal that the copper(ii) centres are antiferromagnetically coupled with the exchange-coupling constant J = -124(1) cm^-1. Theoretically calculated J = -126.51 cm^-1, employing a broken-symmetry DFT approach, is in excellent agreement with the experimental value. The dicopper(ii) complex has been found to be catalytically efficient in the hydrolysis of 2-hydroxypropyl-p-nitrophenylphosphate (HPNP). Detailed kinetic experiments and solution studies (potentiometry, species distribution and ESI-MS) were performed to elucidate the reaction mechanism. DFT calculations were performed to discriminate between different possible mechanistic pathways. The free-energy barrier for HPNP hydrolysis catalyzed by 1 is comparable to that obtained from the experimentally-determined value. The involvement of non-covalent (hydrogen-bonding) interaction has also been probed by DFT calculations. The activity of 1 is found to be the highest, compared to the structurally-characterized Mn, Co, Ni and Zn complexes of L¹(-) reported earlier, under identical experimental conditions, in which each metal centre is six-coordinate.

Magnetic Properties of a Dinuclear Nickel(II) Complex with 2,6-Bis[(2-hydroxyethyl)methylaminomethyl]-4-methylphenolate

Magnetic properties of dinuclear nickel(II) complex [Ni₂(sym-hmp)₂](BPh₄)₂·3.5DMF·0.5(2-PrOH) (1), where (sym-hmp)^- is 2,6-bis[(2-hydroxyethyl)methylaminomethyl]-4-methylphenolate anion and DMF indicates dimethylformamide, were investigated using high-frequency and -field electron paramagnetic resonance (HFEPR). To magnetically characterize the mononuclear nickel(II) species forming the dimer, its two dinuclear zinc(II) analogues, [Zn₂(sym-hmp)₂](BPh₄)₂·3.5DMF·0.5(2-PrOH) (2) and [Zn₂(sym-hmp)₂](BPh₄)₂·2acetone·2H₂O (2'), were prepared. One of them (2') was structurally characterized by X-ray diffractometry and doped with 5% mol nickel(II) ions to prepare a mixed crystal 3. From the HFEPR results on complex 1 obtained at 40 K, the spin Hamiltonian parameters of the first excited spin state (S = 1) of the dimer were accurately determined as |D₁| = 9.99(2) cm^-1, |E₁| = 1.62(1) cm^-1, and g₁ = [2.25(1), 2.19(2), 2.27(2)], and for the second excited spin state (S = 2) at 150 K estimated as |D₂| ≈ 3.5 cm^-1. From these numbers, the single-ion zero-field splitting (ZFS) parameter of the Ni(II) ions forming the dimer was estimated as |D_Ni| ≈ 10-10.5 cm^-1. The HFEPR spectra of 3 yielded directly the single-ion parameters for D_Ni = +10.1 cm^-1, |E_Ni| = 3.1 cm^-1, and g_iso = 2.2. On the basis of the HFEPR results, the previously obtained magnetic data (Sakiyama, H.; Tone, K.; Yamasaki, M.; Mikuriya, M. Inorg. Chim. Acta 2011, 365, 183) were reanalyzed, and the isotropic interaction parameter between the Ni(II) ions was determined as J = -70 cm^-1 (H_ex = -J S_A·S_B). Finally, density functional theory calculations yielded the J value of -90 cm^-1 in a qualitative agreement with the experiment.