Chemical Mimics of Aspartate-Directed Proteases: Predictive and Strictly Specific Hydrolysis of a Globular Protein at Asp-X Sequence Promoted by Polyoxometalate Complexes Rationalized by a Combined Experimental and Theoretical Approach.

Creating efficient and residue-directed artificial proteases is a challenging task due to the extreme inertness of the peptide bond, combined with the difficulty of achieving specific interactions between the catalysts and the protein side chains. Herein we report strictly site-selective hydrolysis of a multi-subunit globular protein, hemoglobin (Hb) from bovine blood, by a range of ZrIV -substituted polyoxometalates (Zr-POMs) in mildly acidic and physiological pH solutions. Among 570 peptide bonds in Hb, selective cleavage was observed at only eleven sites, each occurring at Asp-X peptide bonds located in the positive patches on the protein surface. The molecular origins of the observed Asp-X selectivity were rationalized by means of molecular docking, DFT-based binding, and mechanistic studies on model peptides. The proposed mechanism of hydrolysis involves coordination of the amide oxygen to ZrIV followed by a direct nucleophilic attack of the side chain carboxylate group on the C-terminal amide carbon atom with formation of a cyclic anhydride, which is further hydrolyzed to give the reaction products. The activation energy for the cleavage of the structurally related Glu-X sequence compared to Asp-X was calculated to be higher by 1.4 kcal mol-1 , which corresponds to a difference of about one order of magnitude in the rates of hydrolysis. The higher activation energy is attributed to the higher strain present in the six-membered ring of glutaric anhydride (Glu-X), as compared to the five-membered ring of the succinic anhydride (Asp-X) intermediate. Similarly, the cleavage at X-Asp and X-Glu bonds are predicted to be kinetically less likely as the corresponding activation energies were 6 kcal mol-1 higher, explaining the experimentally observed selectivity. The synergy between the negatively charged polyoxometalate cluster, which binds at positive patches on protein surfaces, and selective activation of Asp-X peptide bonds located in these regions by ZrIV ions, results in a novel class of artificial proteases with aspartate-directed reactivity, which is very rare among naturally occurring proteases.

Phosphate Ester Bond Hydrolysis Promoted by Lanthanide-Substituted Keggin-type Polyoxometalates Studied by a Combined Experimental and Density Functional Theory Approach.

Hydrolytic cleavage of 4-nitrophenyl phosphate (NPP), a commonly used DNA model substrate, was examined in the presence of series of lanthanide-substituted Keggin-type polyoxometalates (POMs) [Me2NH2]11[CeIII(PW11O39)2], [Me2NH2]10[CeIV(PW11O39)2] (abbreviated as (CeIV(PW11)2), and K4[EuPW11O39] by means of NMR and luminescence spectroscopies and density functional theory (DFT) calculations. Among the examined complexes, the Ce(IV)-substituted Keggin POM (CeIV(PW11)2) showed the highest reactivity, and its aqueous speciation was fully determined under different conditions of pD, temperature, concentration, and ionic strength by means of 31P and 31P diffusion-ordered NMR spectroscopy. The cleavage of the phosphoester bond of NPP in the presence of (CeIV(PW11)2) proceeded with an observed rate constant kobs = (5.31 ± 0.06) × 10-6 s-1 at pD 6.4 and 50 °C. The pD dependence of NPP hydrolysis exhibits a bell-shaped profile, with the fastest rate observed at pD 6.4. The formation constant (Kf = 127 M-1) and catalytic rate constant (kc = 19.41 × 10-5 s-1) for the NPP-Ce(IV)-Keggin POM complex were calculated, and binding between CeIV(PW11)2 and the phosphate group of NPP was also evidenced by the change of the chemical shift of the 31P nucleus in NPP upon addition of the POM complex. DFT calculations revealed that binding of NPP to the parent catalyst CeIV(PW11)2 is thermodynamically unlikely. On the contrary, formation of complexes with the monomeric 1:1 species, CeIVPW11, is considered to be more favorable, and the most stable complex, [CeIVPW11(H2O)2(NPP-κO)2]7-, was found to involve two NPP ligands coordinated to the CeIVcenter of CeIVPW11 in the monodentate fashion. The formation of such species is considered to be responsible for the hydrolytic activity of CeIV(PW11)2 toward phosphomonoesters. On the basis of these findings a principle mechanism for the hydrolysis of NPP by the POM is proposed.

Molecular Insight from DFT Computations and Kinetic Measurements into the Steric Factors Influencing Peptide Bond Hydrolysis Catalyzed by a Dimeric Zr(IV)-Substituted Keggin Type Polyoxometalate.

Peptide bond hydrolysis of several peptides with a Gly-X sequence (X = Gly, Ala, Val, Leu, Ile, Phe) catalyzed by a dimeric Zr(IV)-substituted Keggin type polyoxometalate (POM), (Et2NH2)8[{α-PW11O39Zr(µ-OH)(H2O)}2]·7H2O (1), was studied by means of kinetic experiments and (1)H NMR spectroscopy. The observed rate of peptide bond hydrolysis was found to decrease with increase of the side chain bulkiness, from 4.44 × 10(-7) s(-1) for Gly-Gly to 0.81 × 10(-7) s(-1) for Gly-Ile. A thorough DFT investigation was performed to elucidate (a) the nature of the hydrolytically active species in solution, (b) the mechanism of peptide bond hydrolysis, and (c) the influence of the aliphatic residues on the rate of hydrolysis. Formation of substrate-catalyst complexes of the dimeric POM 1 was predicted as thermodynamically unlikely. Instead, the substrates prefer to bind to the monomerization product of 1, [α-PW11O39Zr(OH)(H2O)](4-) (2), which is also present in solution. In the hydrolytically active complex two dipeptide ligands are coordinated to the Zr(IV) center of 2. The first ligand is bidentate-bound through its amino nitrogen and amide oxygen atoms, while the second ligand is monodentate-bound through a carboxylic oxygen atom. The mechanism of hydrolysis involves nucleophilic attack by a solvent water molecule on the amide carbon atom of the bidentate-bound ligand. In this process the uncoordinated carboxylic group of the same ligand acts as a general base to abstract a proton from the attacking water molecule. The decrease of the hydrolysis rate with an increase in the side chain bulkiness is mostly due to the increased ligand conformational strain in the rate-limiting transition state, which elevates the reaction activation energy. The conformational strain increases first upon substitution of Hα in Gly-Gly with the aliphatic α substituent and second with the ß branching of the α substituent.

Multinuclear diffusion NMR spectroscopy and DFT modeling: a powerful combination for unraveling the mechanism of phosphoester bond hydrolysis catalyzed by metal-substituted polyoxometalates.

A detailed reaction mechanism is proposed for the hydrolysis of the phosphoester bonds in the DNA model substrate bis(4-nitrophenyl) phosphate (BNPP) in the presence of the Zr(IV)-substituted Keggin type polyoxometalate (Et2NH2)8[{α-PW11O39Zr(µ-OH)(H2O)}2]⋅7 H2O (ZrK 2:2) at pD 6.4. Low-temperature (31)P DOSY spectra at pD 6.4 gave the first experimental evidence for the presence of ZrK 1:1 in fast equilibrium with ZrK 2:2 in purely aqueous solution. Moreover, theoretical calculations identified the ZrK 1:1 form as the potentially active species in solution. The reaction intermediates involved in the hydrolysis were identified by means of (1)H/(31)P NMR studies, including EXSY and DOSY NMR spectroscopy, which were supported by DFT calculations. This experimental/theoretical approach enabled the determination of the structures of four intermediate species in which the starting compound BNPP, nitrophenyl phosphate (NPP), or the end product phosphate (P) is coordinated to ZrK 1:1. In the proposed reaction mechanism, BNPP initially coordinates to ZrK 1:1 in a monodentate fashion, which results in hydrolysis of the first phosphoester bond in BNPP and formation of NPP. EXSY NMR studies showed that the bidentate complex between NPP and ZrK 1:1 is in equilibrium with monobound and free NPP. Subsequently, hydrolysis of NPP results in P, which is in equilibrium with its monobound form.

A computational study of the glycylserine hydrolysis at physiological pH: a zwitterionic versus anionic mechanism.

The hydrolysis of GlySer at physiological pH was investigated by modeling the most feasible reaction mechanisms in aqueous phase at the MP2/6-311+(2df,2p)//SMD-M06/6-311+(2df,2p) level of the theory. To refine the energies of the most relevant transition states along the reaction paths the cluster-continuum concept was adopted. The hydrolytic process could proceed through two competitive mechanisms involving either the zwitterionic or the anionic form of GlySer. The calculations suggest that at physiological pH the actual mechanism is most probably mixed, anionic-zwitterionic. In this reaction scheme the first stage of N→O acyl transfer involves the anionic form whereas the second stage, during which the resultant ester is hydrolyzed, most likely involves the zwitterionic ester form of GlySer. The energy requirement for the first reaction stage is estimated to be slightly lower than for the second one. The calculated activation parameters (e.g. ΔG(#) = 27.8 kcal mol(-1)) for the nucleophilic addition of a water molecule to the ester carbonyl group of the zwitterionic ester are in good agreement with the experimentally determined values at pD 7.4 (ΔG(#) = 28.7 kcal mol(-1)).

A mechanistic study of the spontaneous hydrolysis of glycylserine as the simplest model for protein self-cleavage.

A common feature of several classes of intrinsically reactive proteins with diverse biological functions is that they undergo self-catalyzed reactions initiated by an N→O or N→S acyl shift of a peptide bond adjacent to a serine, threonine, or cysteine residue. In this study, we examine the N→O acyl shift initiated peptide-bond hydrolysis at the serine residue on a model compound, glycylserine (GlySer), by means of DFT and ab initio methods. In the most favorable rate-determining transition state, the serine COO(-) group acts as a general base to accept a proton from the attacking OH function, which results in oxyoxazolidine ring closure. The calculated activation energy (29.4 kcal mol(-1) ) is in excellent agreement with the experimental value, 29.4 kcal mol(-1) , determined by (1) H NMR measurements. A reaction mechanism for the entire process of GlySer dipeptide hydrolysis is also proposed. In the case of proteins, we found that when no other groups that may act as a general base are available, the N→O acyl shift mechanism might instead involve a water-assisted proton transfer from the attacking serine OH group to the amide oxygen. However, the calculated energy barrier for this process is relatively high (33.6 kcal mol(-1) ), thus indicating that in absence of catalytic factors the peptide bond adjacent to serine is no longer a weak point in the protein backbone. An analogous rearrangement involving the amide N-protonated form, rather than the principle zwitterion form of GlySer, was also considered as a model for the previously proposed mechanism of sea-urchin sperm protein, enterokinase, and agrin (SEA) domain autoproteolysis. The calculated activation energy (14.3 kcal mol(-1) ) is significantly lower than the experimental value reported for SEA (≈21 kcal mol(-1) ), but is still in better agreement as compared to earlier theoretical attempts.