Peptide Hydrolysis by the Binuclear Zinc Enzyme Aminopeptidase from Aeromonas proteolytica: A Density Functional Theory Study

Nomogram for preoperative estimation of microvascular invasion risk in hepatocellular carcinoma

Microvascular invasion (MVI) is an adverse prognostic indicator of tumor recurrence after surgery for hepatocellular carcinoma (HCC). Therefore, developing a nomogram for estimating the presence of MVI before liver resection is necessary. We retrospectively included 260 patients with pathologically confirmed HCC at the Fifth Medical Center of Chinese PLA General Hospital between January 2021 and April 2024. The patients were randomly divided into a training cohort (n = 182) for nomogram development, and a validation cohort (n = 78) to confirm the performance of the model (7:3 ratio). Significant clinical variables associated with MVI were then incorporated into the predictive nomogram using both univariate and multivariate logistic analyses. The predictive performance of the nomogram was assessed based on its discrimination, calibration, and clinical utility. Serum carnosine dipeptidase 1 ([CNDP1] OR 2.973; 95 % CI 1.167-7.575; p = 0.022), cirrhosis (OR 8.911; 95 % CI 1.922-41.318; p = 0.005), multiple tumors (OR 4.095; 95 % CI 1.374-12.205; p = 0.011), and tumor diameter ≥3 cm (OR 4.408; 95 % CI 1.780-10.919; p = 0.001) were independent predictors of MVI. Performance of the nomogram based on serum CNDP1, cirrhosis, number of tumors and tumor diameter was achieved with a concordance index of 0.833 (95 % CI 0.771-0.894) and 0.821 (95 % CI 0.720-0.922) in the training and validation cohorts, respectively. It fitted well in the calibration curves, and the decision curve analysis further confirmed its clinical usefulness. The nomogram, incorporating significant clinical variables and imaging features, successfully predicted the personalized risk of MVI in HCC preoperatively.

Comprehensive pan-cancer investigation of carnosine dipeptidase 1 and its prospective prognostic significance in hepatocellular carcinoma

Carnosine dipeptidase 1 (CNDP1), an enzyme integral to the hydrolysis of dipeptides containing histidine, plays an indispensable role in myriad physiological processes, including hydrolysis of proteins, maturation of specific biochemical functionalities within proteins, tissue regeneration, and regulation of cell cycle. However, the implications of CNDP1 in oncogenesis and its prognostic value are not yet fully elucidated. Initially, we procured the GSE40367 dataset from the Gene Expression Omnibus and established a protein-protein interaction network. Thereafter, we conducted functional and pathway enrichment analyses utilizing GO, KEGG, and GSEA. Moreover, we undertook an association analysis concerning the expression of CNDP1 with immune infiltration, along with survival analysis across various cancers and specifically in hepatocellular carcinoma (HCC). Our study uncovered a total of 2,248 differentially expressed genes, with a down-regulation of CNDP1 in HCC and other cancers. Our explorations into the relationship between CNDP1 and immune infiltration disclosed a negative correlation between CNDP1 expression and the presence of immune cells in HCC. Survival analyses revealed that diminished expression of CNDP1 correlates with an adverse prognosis in HCC and several other types of cancer. These observations intimate that CNDP1 holds promise as a novel prognostic biomarker for both pan-cancer and HCC.

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.

Benchmarking Density Functional Theory Methods for Metalloenzyme Reactions: The Introduction of the MME55 Set

We present a new benchmark set of metalloenzyme model reaction energies and barrier heights that we call MME55. The set contains 10 different enzymes, representing eight transition metals, both open and closed shell systems, and system sizes of up to 116 atoms. We use four DLPNO-CCSD(T)-based approaches to calculate reference values against which we then benchmark the performance of a range of density functional approximations with and without dispersion corrections. Dispersion corrections improve the results across the board, and triple-ζ basis sets provide the best balance of efficiency and accuracy. Jacob's ladder is reproduced for the whole set based on averaged mean absolute (percent) deviations, with the double hybrids SOS0-PBE0-2-D3(BJ) and revDOD-PBEP86-D4 standing out as the most accurate methods for the MME55 set. The range-separated hybrids ωB97M-V and ωB97X-V also perform well here and can be recommended as a reliable compromise between accuracy and efficiency; they have already been shown to be robust across many other types of chemical problems, as well. Despite the popularity of B3LYP in computational enzymology, it is not a strong performer on our benchmark set, and we discourage its use for enzyme energetics.

Binding Modes of Carnostatine, Homocarnosine, and Ophidine to Human Carnosinase 1

Carnosine (CAR), anserine (ANS), homocarnosine (H-CAR), and ophidine (OPH) are histidine-containing dipeptides that show a wide range of therapeutic properties. With their potential physiological effects, these bioactive dipeptides are considered as bioactive food components. However, such dipeptides display low stability due to their rapid degradation by human serum carnosinase 1 (CN1). A dimeric CN1 hydrolyzes such histidine-containing compounds with different degrees of reactivities. A selective CN inhibitor, carnostatine (CARN), was reported to effectively inhibit CN's activity. To date, the binding mechanisms of CAR and ANS have been recently reported, while no clear information about H-CAR, OPH, and CARN binding is available. Thus, in this work, molecular dynamics simulations were employed to elucidate the binding mechanism of H-CAR, OPH, and CARN. Among all, the amine end and imidazole ring are the main players for trapping all of the ligands in a pocket. OPH shows the poorest binding affinity, while CARN displays the tightest binding. Such firm binding is due to the longer amine chain and the additional hydroxyl (-OH) group of CARN. H-CAR and CARN are analogous, but the absence of the -OH moiety in H-CAR significantly enhances its mobility, resulting in the reduction in binding affinity. For OPH which is an ANS analogue, the methylated imidazole ring destroys the OPH-CN1 interaction network at this region, consequentially leading to the poor binding ability. An insight into how CN recognizes and binds its substrates obtained here will be useful for designing an effective strategy to prolong the lifetime of CAR and its analogues after ingestion.

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.

Exploring the structural and dynamic differences between human carnosinase I (CN1) and II (CN2)

Human carnosinases (CNs) are dimeric dipeptidases in the metallopeptidase M20 family. Two isoforms of carnosinases (Zn²⁺ -containing carnosinase 1 (CN1) found in serum and Mn²⁺ -carnosinase 2 (CN2) in tissue) were identified. Both CNs cleave histidine-containing (Xaa-His) dipeptides such as carnosine where CN2 was found to accept a broader spectrum of substrates. A loss of CN function, resulting in a high carnosine concentration, reduces risk for diabetes and neurological disorders. Although several studies on CN activities and its Michaelis complex were conducted, all shed the light on CN1 activity where the CN2 data is limited. Also, the molecular details on CN1 and CN2 similarity and dissimilarity in structure and function remain unclear. Thus, in this work, molecular dynamics (MD) simulations were employed to study structure and dynamics of human CN1 and CN2 in comparison. The results show that the different catalytic ability of both CNs is due to their pocket size and environment. CN2 can accept a wider range of substrate due to the wider mouth of a binding pocket. The L1 loop seems to play a role in gating activity. Comparing to CN2, CN1 provides more electronegative entrance, more wettability, and higher stability of catalytic metal ion-pair in the active site which allow more efficient water-mediated catalysis. The microscopic understanding obtained here can serve as a basis for CN inhibition strategies resulting in higher carnosine levels and consequently mitigating complications associated with diseases such as diabetes and neurological disorder.

Molecular insights into the binding of carnosine and anserine to human serum carnosinase 1 (CN1)

Carnosine (CAR) and anserine (ANS) are histidine-containing dipeptides that show the therapeutic properties and protective abilities against diabetes and cognitive deficit. Both dipeptides are rich in meat products and have been used as a supplement. However, in humans, both compounds have a short half-life due to the rapid degradation by dizinc carnosinase 1 (CN1) which is a hurdle for its therapeutic application. To date, a comparative study of carnosine- and anserine-CN1 complexes is limited. Thus, in this work, molecular dynamics (MD) simulations were performed to explore the binding of carnosine and anserine to CN1. CN1 comprises 2 chains (Chains A and B). Both monomers are found to work independently and alternatingly. The displacement of Zn2+ pair is found to disrupt the substrate binding. CN1 employs residues from the neighbour chain (H235, T335, and T337) to form the active site. This highlights the importance of a dimer for enzymatic activity. Anserine is more resistant to CN 1 than carnosine because of its bulky and dehydrated imidazole moiety. Although both dipeptides can direct the peptide oxygen to the active Zn2+ which can facilitate the catalytic reaction, the bulky methylated imidazole on anserine promotes various poses that can retard the hydrolytic activity in contrast to carnosine. Anserine is likely to be the temporary competitive inhibitor by retarding the carnosine catabolism.

Engineered ACE2-Fc counters murine lethal SARS-CoV-2 infection through direct neutralization and Fc-effector activities

Soluble angiotensin-converting enzyme 2 (ACE2) constitutes an attractive antiviral capable of targeting a wide range of coronaviruses using ACE2 as their receptor. Using structure-guided approaches, we developed a series of bivalent ACE2-Fcs harboring functionally and structurally validated mutations that enhance severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) receptor binding domain recognition by up to ~12-fold and remove angiotensin enzymatic activity. The lead variant M81 potently cross-neutralized SARS-CoV-2 variants of concern (VOCs), including Omicron, at subnanomolar half-maximal inhibitory concentration and was capable of robust Fc-effector functions, including antibody-dependent cellular cytotoxicity, phagocytosis, and complement deposition. When tested in a stringent K18-hACE2 mouse model, Fc-enhanced ACE2-Fc delayed death by 3 to 5 days or effectively resolved lethal SARS-CoV-2 infection in both prophylactic and therapeutic settings via the combined effects of neutralization and Fc-effector functions. These data add to the demonstrated utility of soluble ACE2 as a valuable SARS-CoV-2 antiviral and indicate that Fc-effector functions may constitute an important component of ACE2-Fc therapeutic activity.

Promiscuous Catalytic Activity of a Binuclear Metallohydrolase: Peptide and Phosphoester Hydrolyses

In this study, chemical promiscuity of a binuclear metallohydrolase Streptomyces griseus aminopeptidase (SgAP) has been investigated using DFT calculations. SgAP catalyzes two diverse reactions, peptide and phosphoester hydrolyses, using its binuclear (Zn-Zn) core. On the basis of the experimental information, mechanisms of these reactions have been investigated utilizing leucine p-nitro aniline (Leu-pNA) and bis(4-nitrophenyl) phosphate (BNPP) as the substrates. The computed barriers of 16.5 and 16.8 kcal/mol for the most plausible mechanisms proposed by the DFT calculations are in good agreement with the measured values of 13.9 and 18.3 kcal/mol for the Leu-pNA and BNPP hydrolyses, respectively. The former was found to occur through the transfer of two protons, while the latter with only one proton transfer. They are in line with the experimental observations. The cleavage of the peptide bond was the rate-determining process for the Leu-pNA hydrolysis. However, the creation of the nucleophile and its attack on the electrophile phosphorus atom was the rate-determining step for the BNPP hydrolysis. These calculations showed that the chemical nature of the substrate and its binding mode influence the nucleophilicity of the metal bound hydroxyl nucleophile. Additionally, the nucleophilicity was found to be critical for the Leu-pNA hydrolysis, whereas double Lewis acid activation was needed for the BNPP hydrolysis. That could be one of the reasons why peptide hydrolysis can be catalyzed by both mononuclear and binuclear metal cofactors containing hydrolases, while phosphoester hydrolysis is almost exclusively by binuclear metallohydrolases. These results will be helpful in the development of versatile catalysts for chemically distinct hydrolytic reactions.

Predicting Effects of Site-Directed Mutagenesis on Enzyme Kinetics by QM/MM and QM Calculations: A Case of Glutamate Carboxypeptidase II

Quantum and molecular mechanics (QM/MM) and QM-only (cluster model) modeling techniques represent the two workhorses in mechanistic understanding of enzyme catalysis. One of the stringent tests for QM/MM and/or QM approaches is to provide quantitative answers to real-world biochemical questions, such as the effect of single-point mutations on enzyme kinetics. This translates into predicting the relative activation energies to 1-2 kcal·mol^-1 accuracy; such predictions can be used for the rational design of novel enzyme variants with desired/improved characteristics. Herein, we employ glutamate carboxypeptidase II (GCPII), a dizinc metallopeptidase, also known as the prostate specific membrane antigen, as a model system. The structure and activity of this major cancer antigen have been thoroughly studied, both experimentally and computationally, which makes it an ideal model system for method development. Its reaction mechanism is quite well understood: the reaction coordinate comprises a "tetrahedral intermediate" and two transition states and experimental activation Gibbs free energy of ∼17.5 kcal·mol^-1 can be inferred for the known k_cat ≈ 1 s^-1. We correlate experimental kinetic data (including the E424H variant, newly characterized in this work) for various GCPII mutants (k_cat = 8.6 × 10^-5 s^-1 to 2.7 s^-1) with the energy profiles calculated by QM/MM and QM-only (cluster model) approaches. We show that the near-quantitative agreement between the experimental values and the calculated activation energies (ΔH^⧧) can be obtained and recommend the combination of the two protocols: QM/MM optimized structures and cluster model (QM) energetics. The trend in relative activation energies is mostly independent of the QM method (DFT functional) used. Last but not least, a satisfactory correlation between experimental and theoretical data allows us to provide qualitative and fairly simple explanations of the observed kinetic effects which are thus based on a rigorous footing.

Engineered ACE2-Fc counters murine lethal SARS-CoV-2 infection through direct neutralization and Fc-effector activities

Soluble Angiotensin-Converting Enzyme 2 (ACE2) constitutes an attractive antiviral capable of targeting a wide range of coronaviruses utilizing ACE2 as their receptor. Here, using structure-guided approaches, we developed divalent ACE2 molecules by grafting the extracellular ACE2-domain onto a human IgG1 or IgG3 (ACE2-Fc). These ACE2-Fcs harbor structurally validated mutations that enhance spike (S) binding and remove angiotensin enzymatic activity. The lead variant bound tightly to S, mediated in vitro neutralization of SARS-CoV-2 variants of concern (VOCs) with sub-nanomolar IC 50 and was capable of robust Fc-effector functions, including antibody-dependent-cellular cytotoxicity, phagocytosis and complement deposition. When tested in a stringent K18-hACE2 mouse model, it delayed death or effectively resolved lethal SARS-CoV-2 infection in a prophylactic or therapeutic setting utilizing the combined effect of neutralization and Fc-effector functions. These data confirm the utility of ACE2-Fcs as valuable agents in preventing and eliminating SARS-CoV-2 infection and demonstrate that ACE2-Fc therapeutic activity require Fc-effector functions.

A novel binuclear Pd(ii) complex displaying synergic peptide cleavage behaviour

Mononuclear Pd(ii) complexes with two leaving groups are able to promote His-, Cys- and Met-orientated peptide hydrolysis, and exploring the peptide cleavage behavior of a novel Pd(ii) complex may provide "Omics" studies a promising artificial protease. In this study, a novel binuclear Pd(ii) complex [Pd₂(μ-O-L-H)(μ-Cl)](ClO₄)₂ (L = 2,6-bis(N-2'-aminoethylaminomethyl)-p-cresol) was constructed to promote peptide hydrolysis. Although each Pd(ii) center has only one leaving group (Cl) in this complex, electrophoresis and LC-MS-MS determination discloses that this complex enables myoglobin cleavage on the second upstream peptide bond from His and Met. A study on peptide cleavage also confirms the His- and Met-orientated peptide hydrolysis, yet no Cys-orientated hydrolysis was observed, although the cysteine-induced peptide/complex binding is distinct. Cysteine in the peptide even prevents the complex from promoting His-orientated hydrolysis, whereas the oxidized cysteine residue recovers the His-orientated hydrolysis. This peptide cleavage behavior is quite different from the simultaneous His-, Cys-, and Met-orientated hydrolysis promoted by the mononuclear Pd(ii) complexes. A theoretical study suggests that the two Pd(ii) centers of this complex might promote His- and Met-orientated hydrolysis in a synergic manner: one Pd(ii) center binds selectively on peptides or proteins and the other coordinates with the amide bond and water favoring nucleophilic attack on the peptide bond. The thiol group of cysteine is inclined to bridge the two Pd(ii) centers to form a "closed" sulphur-bridged structure, disfavoring the Cys-orientated hydrolysis. This study not only demonstrates the peptide cleavage behavior of this binuclear Pd(ii) complex, but also provides a polynuclear strategy to regulate the peptide cleavage behavior of Pd(ii) complexes.

Using Atomic Confining Potentials for Geometry Optimization and Vibrational Frequency Calculations in Quantum-Chemical Models of Enzyme Active Sites

Quantum-chemical studies of enzymatic reaction mechanisms sometimes use truncated active-site models as simplified alternatives to mixed quantum mechanics molecular mechanics (QM/MM) procedures. Eliminating the MM degrees of freedom reduces the complexity of the sampling problem, but the trade-off is the need to introduce geometric constraints in order to prevent structural collapse of the model system during geometry optimizations that do not contain a full protein backbone. These constraints may impair the efficiency of the optimization, and care must be taken to avoid artifacts such as imaginary vibrational frequencies. We introduce a simple alternative in which terminal atoms of the model system are placed in soft harmonic confining potentials rather than being rigidly constrained. This modification is simple to implement and straightforward to use in vibrational frequency calculations, unlike iterative constraint-satisfaction algorithms, and allows the optimization to proceed without constraint even though the practical result is to fix the anchor atoms in space. The new approach is more efficient for optimizing minima and transition states, as compared to the use of fixed-atom constraints, and also more robust against unwanted imaginary frequencies. We illustrate the method by application to several enzymatic reaction pathways where entropy makes a significant contribution to the relevant reaction barriers. The use of confining potentials correctly describes reaction paths and facilitates calculation of both vibrational zero-point and finite-temperature entropic corrections to barrier heights.

Mechanism and Inhibitor Exploration with Binuclear Mg Ketol-Acid Reductoisomerase: Targeting the Biosynthetic Pathway of Branched-Chain Amino Acids

Binuclear Mg ketol-acid reductoisomerase (KARI), which converts (S)-2-acetolactate into (R)-2,3-dihydroxyisovalerate, is responsible for the second step of the biosynthesis of branched-chain amino acids in plants and microorganisms and thus serves as a key inhibition target potentially without effects on mammals. Here, through the use of density functional calculations and a chemical model, the KARI-catalyzed reaction has been demonstrated to include the initial deprotonation of the substrate C2 hydroxy group, bridged by the two Mg ions, alkyl migration from the C2-alkoxide carbon atom to the C3-carbonyl carbon atom, and hydride transfer from a nicotinamide adenine dinucleotide phosphate [NAD(P)H] cofactor to C2. A dead-end mechanism with a hydride transferred to the C3 carbonyl group has been ruled out. The nucleophilicity (migratory aptitude) of the migrating carbon atom and the provision of additional negative charge to the di-Mg coordination sphere have significant effects on the steps of alkyl migration and hydride transfer, respectively. Other important mechanistic characteristics are also revealed. Inspired by the mechanism, an inhibitor (2-carboxylate-lactic acid) was designed and predicted by barrier analysis to be effective in inactivating KARI, hence probably enriching the antifungal and antibacterial library. Two types of slow substrate analogues (2-trihalomethyl acetolactic acids and 2-glutaryl lactic acid) were also found.

How does Mo-dependent perchlorate reductase work in the decomposition of oxyanions?

The mechanisms of Mo-dependent perchlorate reductase (PcrAB)-catalyzed decomposition of perchlorate, bromate, iodate, and nitrate were revealed by density functional calculations.

How does binuclear zinc amidohydrolase FwdA work in the initial step of methanogenesis: From formate to formyl-methanofuran

The initial step of methanogenesis is the fixation of CO₂ to formyl-methanofuran (formyl-MFR) catalyzed by formyl-MFR dehydrogenase, which can be divided into two half reactions. Herein, the second half reaction catalyzed by FwdA (formyl-methanofuran dehydrogenase subunit A), i.e., from formate to formyl-methanofuran, has been investigated using density functional theory and a chemical model based on the X-ray crystal structure. The calculations indicate that, compared with other well-known di-zinc hydrolases, the FwdA reaction employs a reverse mechanism, including the nucleophilic attack of MFR amine on formate carbon leading to a tetrahedral gem-diolate intermediate, two steps of proton transfer from amine to formate moieties assisted by the Asp385, and the CO bond dissociation to form the formyl-MFR product. The second step of proton transfer from the amine moiety to the Asp385 is rate-limiting with an overall barrier of 21.2 kcal/mol. The two zinc ions play an important role in stabilizing the transition states and intermediates, in particular the negative charge at the formate moiety originated from the nucleophilic attack of the MFR amine. The work here appends a crucial piece in the methanogenic mechanistics and advances the understanding of the global carbon cycle.

Evidence for Cleavage of the Metalloprotease Vsm from Vibrio splendidus Strain JZ6 by an M20 Peptidase (PepT-like Protein) at Low Temperature

Metalloprotease Vsm is a major extracellular virulence factor of Vibrio splendidus. The toxicity of Vsm from V. splendidus strain JZ6 has been characterized, and production of this virulence factor proved to be temperature-regulated. The present study provides evidence that two forms (JZE1 and JZE2) of Vsm protein exist in extracellular products (ECPs) of strain JZ6, and a significant conversion of these two forms was detected by SDS-PAGE and immunoblotting analyses of samples obtained from cells grown at 4, 10, 16, 20, 24, and 28°C. Mass spectroscopy confirmed that JZE1 was composed only of the peptidase_M4 domain of Vsm, and JZE2 contained both the PepSY domain and the peptidase_M4 domain. An M20 peptidase T-like protein (PepTL) was screened from the transcriptome data of strain JZ6, which was considered as a crucial molecule to produce the active Vsm (JZE1) by cleavage of the propeptide. Similar to that of Vsm, PepTL mRNA accumulation was highest at 4°C (836.82-fold of that at 28°C), decreased with increasing of temperature and reached its lowest level at 28°C. Deletion of the gene encoding the PepTL resulted in a mutant strain that did not produce the JZE1 cleavage product. The peptidase activity of PepTL recombinant protein (rPepTL) was confirmed by cleaving the Vsm in ECPs with an in vitro degradation reaction. These results demonstrate that PepTL participates in activating Vsm in strain JZ6 by proteolytic cleavage at low temperature.

Mechanisms of peptide hydrolysis by aspartyl and metalloproteases

Peptide hydrolysis has been involved in a wide range of biological, biotechnological, and industrial applications.

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.

Enhanced thermal stability and hydrolytic ability of Bacillus subtilis aminopeptidase by removing the thermal sensitive domain in the non-catalytic region

Besides the catalytic ability, many enzymes contain conserved domains to perform some other physiological functions. However, sometimes these conserved domains were unnecessary or even detrimental to the catalytic process for industrial application of the enzymes. In this study, based on homology modeling and molecular dynamics simulations, we found that Bacillus subtilis aminopeptidase contained a thermal sensitive domain (protease-associated domain) in the non-catalytic region, and predicted that deletion of this flexible domain can enhance the structure stability. This prediction was then verified by the deletion of protease-associated domain from the wild-type enzyme. The thermal stability analysis showed that deletion of this domain improved the T50 (the temperature required to reduce initial activity by 50% in 30 min) of the enzyme from 71 °C to 77 °C. The melting temperature (Tm) of the enzyme also increased, which was measured by thermal denaturation experiments using circular dichroism spectroscopy. Further studies indicated that this deletion did not affect the activity and specificity of the enzyme toward aminoacyl-p-nitroanilines, but improved its hydrolytic ability toward a 12-aa-long peptide (LKRLKRFLKRLK) and soybean protein. These findings suggested the possibility of a simple technique for enzyme modification and the artificial enzyme obtained here was more suitable for the protein hydrolysis in food industry than the wild-type enzyme.

Carnosinases, their substrates and diseases

Carnosinases are Xaa-His dipeptidases that play diverse functions throughout all kingdoms of life. Human isoforms of carnosinase (CN1 and CN2) under appropriate conditions catalyze the hydrolysis of the dipeptides carnosine (β-alanyl-l-histidine) and homocarnosine (γ-aminobutyryl-l-histidine). Alterations of serum carnosinase (CN1) activity has been associated with several pathological conditions, such as neurological disorders, chronic diseases and cancer. For this reason the use of carnosinase levels as a biomarker in cerebrospinal fluid (CSF) has been questioned. The hydrolysis of imidazole-related dipeptides in prokaryotes and eukaryotes is also catalyzed by aminoacyl-histidine dipeptidases like PepD (EC 3.4.13.3), PepV (EC 3.4.13.19) and anserinase (EC 3.4.13.5). The review deals with the structure and function of this class of enzymes in physiological and pathological conditions. The main substrates of these enzymes, i.e., carnosine, homocarnosine and anserine (β-alanyl-3-methyl-l-histidine) will also be described.

A dominant homolytic O-Cl bond cleavage with low-spin triplet-state Fe(IV)=O formed is revealed in the mechanism of heme-dependent chlorite dismutase

Chlorite dismutase (Cld) is a heme-dependent enzyme that catalyzes the decomposition of toxic chlorite (ClO2(-)) into innocuous chloride and O2. In this paper, using the hybrid B3LYP density functional theory (DFT) method including dispersion interactions, the Cld reaction mechanism has been studied with a chemical model constructed on the X-ray crystal structure. The calculations indicate that the reaction proceeds along a stepwise pathway in the doublet state, i.e. a homolytic O-Cl bond cleavage of the substrate leading to an O-Fe(heme) species and a ClO˙ radical, followed by a rebinding O-O bond formation between them. The O-Fe(heme) species is demonstrated to be a low-spin triplet-state Fe(IV)=O diradicaloid. A low-spin singlet-state Fe(IV)=O is much less stable than the former, with an energy difference of 9.2 kcal mol(-1). The O-Cl bond cleavage is rate-limiting with a barrier of 10.6 kcal mol(-1), in good agreement with the experimental reaction rate of 2.0 × 10(5) s(-1). Furthermore, a heterolytic O-Cl bond dissociation in the initial step is shown to be unreachable, which ensures the high efficiency of the Cld enzyme by avoiding the generation of chlorate byproduct observed in the reactions of synthetic Fe porphyrins. Also, the pathways in the quartet and sextet states are unfavorable for the Cld reaction. The present results reveal a detailed mechanism III (defined in the text) including an interesting di-radical intermediate composed of a low-spin triplet-state Fe(IV)=O and a ClO˙ radical. Compared to a competitive heterolytic Cl-O cleavage in synthetic Fe porphyrins, the revelation of the domination of homolysis in Cld indicates not only the high efficiency of enzyme, but also the sensitivity of a heme and the significance of the enzymatic active-site surroundings (the His170 and Arg183 residues in the present case), which gives more insights into heme chemistry.

Structure-based approach to alter the substrate specificity of Bacillus subtilis aminopeptidase

Aminopeptidases can selectively catalyze the cleavage of the N-terminal amino acid residues from peptides and proteins. Bacillus subtilis aminopeptidase (BSAP) is most active toward p-nitroanilides (pNAs) derivatives of Leu, Arg, and Lys. The BSAP with broad substrate specificity is expected to improve its application. Based on an analysis of the predicted structure of BSAP, four residues (Leu 370, Asn 385, Ile 387, and Val 396) located in the substrate binding region were selected for saturation mutagenesis. The hydrolytic activity toward different aminoacyl-pNAs of each mutant BSAP in the culture supernatant was measured. Although the mutations resulted in a decrease of hydrolytic activity toward Leu-pNA, N385L BSAP exhibited higher hydrolytic activities toward Lys-pNA (2.2-fold) and Ile-pNA (9.1-fold) than wild-type BSAP. Three mutant enzymes (I387A, I387C and I387S BSAPs) specially hydrolyzed Phe-pNA, which was undetectable in wild-type BSAP. Among these mutant BSAPs, N385L and I387A BSAPs were selected for further characterized and used for protein hydrolysis application. Both of N385L and I387A BSAPs showed higher hydrolysis efficiency than the wild-type BASP and a combination of the wild-type and N385L and I387A BSAPs exhibited the highest hydrolysis efficiency for protein hydrolysis. This study will greatly facilitate studies aimed on change the substrate specificity and our results obtained here should be useful for BSAP application in food industry.

What occurs by replacing Mn2+ with Co2+ in human arginase I: first-principles computational analysis

The reaction mechanism of the dinuclear cobalt enzyme arginase is investigated using density functional theory. As an arginase-containing binuclear Mn(2)(2+) cluster, it catalyzes the hydrolysis of L-arginine in L-ornithine and urea. The bridging hydroxide is capable of performing nucleophilic attack on the iminium carbon ion from its bridging position, resulting in the formation of a tetrahedral intermediate, as was already obtained in a previous theoretical study on the manganese enzyme. Our theoretical investigation allows us to obtain an accurate potential energy profile and confirms that the coordination mode of the substrate to the dimetallic center is quite similar to that present in the manganese enzyme. In agreement with the experimental observations, our results show that both Mn- and Co-containing enzymes catalyze the same reaction with quite comparable energy barriers.

A quantum-chemical study of the binding ability of βXaaHisGlyHis towards copper(II) ion

The present study analyzed binding of Cu²⁺ to tetrapeptides in water solution at several levels of theoretical approximation. The methods used to study the energetic and structural properties of the complexes in question include semiempirical hamiltonians, density functional theory as well as ab initio approaches including electron correlation effects. In order to shed light on the character of interactions between Cu²⁺ and peptides, which are expected to be mainly electrostatic in nature, decomposition of interaction energy into physically meaningful components was applied.

Mechanism of peptide hydrolysis by co-catalytic metal centers containing leucine aminopeptidase enzyme: a DFT approach

In this density functional theory study, reaction mechanisms of a co-catalytic binuclear metal center (Zn1-Zn2) containing enzyme leucine aminopeptidase for two different metal bridging nucleophiles (H(2)O and -OH) have been investigated. In addition, the effects of the substrate (L-leucine-p-nitroanilide → L-leucyl-p-anisidine) and metal (Zn1 → Mg and Zn2 → Co, i.e., Mg1-Zn2 and Mg1-Co2 variants) substitutions on the energetics of the mechanism have been investigated. The general acid/base mechanism utilizing a bicarbonate ion followed by this enzyme is divided into two steps: (1) the formation of the gem-diolate intermediate, and (2) the cleavage of the peptide bond. With the computed barrier of 17.8 kcal/mol, the mechanism utilizing a hydroxyl nucleophile was found to be in excellent agreement with the experimentally measured barrier of 18.7 kcal/mol. The rate-limiting step for reaction with L-leucine-p-nitroanilide is the cleavage of the peptide bond with a barrier of 17.8 kcal/mol. However, for L-leucyl-p-anisidine all steps of the mechanism were found to occur with similar barriers (18.0-19.0 kcal/mol). For the metallovariants, cleavage of the peptide bond occurs in the rate-limiting step with barriers of 17.8, 18.0, and 24.2 kcal/mol for the Zn1-Zn2, Mg1-Zn2, and Mg1-Co2 enzymes, respectively. The nature of the metal ion was found to affect only the creation of the gem-diolate intermediate, and after that all three enzymes follow essentially the same energetics. The results reported in this study have elucidated specific roles of both metal centers, the nucleophile, indirect ligands, and substrates in the catalytic functioning of this important class of binuclear metallopeptidases.

Protonation of a hydroxide anion bridging two divalent magnesium cations in water probed by first-principles metadynamics simulation

The protonation of a hydroxide anion (OH(-)) located between two magnesium cations (Mg(2+)) in aqueous solution has been investigated by first-principles metadynamics simulation. We observe that the complex Mg(2+)-OH(-)-Mg(2+) is stabilized by the coparticipation of the hydroxide anion to the first hydration shells of both the Mg(2+) cations. Contrary to the cases of OH(-) in pure water, the transfer of protons in the presence of the divalent metal ions turns out to be a slow chemical event. This can be ascribed to the decreased proton affinity of the bridging OH(-). Metadynamics simulation, used to overcome the difficulty of the long time scale required by the protonation of the bridging OH(-), has shown that the system possesses a great stability on the reactant state, characterized by a bioctahedral (6,6) solvation structure around the two Mg(2+) cations. The exploration of the free energy landscape shows that this stable bioctahedral configuration converts into a lower coordinated (5,6) structure, leading to a proton transfer from a water molecule belonging to the first solvation shell of the Mg(2+) ion having the lower coordination to the bridging OH(-); the free energy barrier for the protonation reaction is 11 kcal/mol, meaning that the bridging hydroxide is a weak base. During the proton transfer, the bridging OH(-) reverts to an H(2)O molecule, and this breaks the electrostatic coupling of the two Mg(2+) ions, which depart independently with their own hydration shells, one of which is entirely formed by water molecules. The second one carries the newly created OH(-). Our results show that the flexibility in the metal coordination plays a crucial role in both the protonation process of the bridging OH(-) and the separation of the metal cations, providing useful insight into the nature of proton transfer in binuclear divalent metal ions, with several biological implications, such as, for instance, transesterification of catalytic RNA.