Structure and mechanism of the pepsin-like family of aspartic peptidases

Characterization of a novel aspartic protease from Trichoderma asperellum for the preparation of duck blood peptides

A novel aspartic protease gene (TaproA1) from Trichoderma asperellum was successfully expressed in Komagataella phaffii (Pichia pastoris). TaproA1 showed 52.8% amino acid sequence identity with the aspartic protease PEP3 from Coccidioides posadasii C735. TaproA1 was efficiently produced in a 5 L fermenter with a protease activity of 4092 U/mL. It exhibited optimal reaction conditions at pH 3.0 and 50 °C and was stable within pH 3.0-6.0 and at temperatures up to 45 °C. The protease exhibited broad substrate specificity with high hydrolysis activity towards myoglobin and hemoglobin. Furthermore, duck blood proteins (hemoglobin and plasma protein) were hydrolyzed by TaproA1 to prepare bioactive peptides with high ACE inhibitory activity. The IC50 values of hemoglobin and plasma protein hydrolysates from duck blood proteins were 0.105 mg/mL and 0.091 mg/mL, respectively. Thus, the high yield and excellent biochemical characterization of TaproA1 presented here make it a potential candidate for the preparation of duck blood peptides. KEY POINTS: • An aspartic protease (TaproA1) from Trichoderma asperellum was expressed in Komagataella phaffii. • TaproA1 exhibited broad substrate specificity and the highest activity towards myoglobin and hemoglobin. • TaproA1 has great potential for the preparation of bioactive peptides from duck blood proteins.

Selective Protonation of Catalytic Dyad for γ-Secretase-Mediated Hydrolysis Revealed by Multiscale Simulations

γ-Secretase plays a crucial role in producing disease-related amyloid-β proteins by cleaving the amyloid precursor protein (APP). The enzyme employs its catalytic dyad containing two aspartates (Asp257 and Asp385) to hydrolyze the substrate by a general acid-base catalytic mechanism, necessitating monoprotonation of the two aspartates for efficient hydrolysis. However, the precise protonation states of the aspartates remain uncertain. In this study, we employed a multiscale computational approach to investigate the dependence of the catalytic efficiency of γ-secretase on the protonation states of its catalytic dyad. Over 200 ms unbiased atomistic simulations of the substrate-enzyme complex reveal diverse orientations of the scissile bond of the bound substrate and accessible structural ensembles of the catalytic dyad with Asp257-Asp385 distances fluctuating between 4 and 10 Å. With a quantum mechanics/molecular mechanics (QM/MM) approach accelerated by enhanced sampling techniques, we find that the first step of the hydrolysis reaction, i.e., the formation of a gem-diol intermediate, experiences a higher reaction barrier by ∼2 kcal/mol when Asp385 is protonated. Furthermore, we find that Arg269 of the enzyme is most likely responsible for this preference of the protonation state: its basic side chain is spatially close to that of Asp257 and specifically stabilizes the transition state electrostatically when Asp257 is protonated. Collectively, our study suggests that Asp257 is likely the favored protonation site for APP cleavage by γ-secretase.

"Blade of Polarized Water Molecule" Is the Key to Hydrolase Catalysis Regulation

Hydrolysis catalyzed by aspartic proteases is a crucial reaction in many biological processes. However, anchoring water molecules and unifying multiple catalytic pathways remain significant challenges. Consequently, molecular design often compromises by focusing on enhancing substrate specificity. Using our self-developed polarizable point charge (PPC) force field, we determined the significant role of polarization in the hydrolase of pepsin for the first time. To be stably anchored in the active site, the water should be intensely polarized with a charge higher than -0.94e. Induced by this polarization, the pepsin was shown to support three general base/general acid pathways, with a preference for the gemdiol-intermediate-based pathway. Consequently, we proposed the "Blade of Polarized Water Molecule" model for rational enzyme design, highlighting that the polarization of both the attacking water and the attacked carbonyl is crucial for enhancing hydrolysis. Mutants D290Q and S172P showed activity enhancements of 191.23% and 324.70%, respectively. The improved polarization of water, carbonyl, and relevant nucleophilic attack distances in the mutants reaffirmed the crucial role of polarization in improving hydrolysis. This study provides a new perspective on hydrolase analysis and modification.

First Insight into the Degradome of Aspergillus ochraceus: Novel Secreted Peptidases and Their Inhibitors

Aspergillus fungi constitute a pivotal element within ecosystems, serving as both contributors of biologically active compounds and harboring the potential to cause various diseases across living organisms. The organism's proteolytic enzyme complex, termed the degradome, acts as an intermediary in its dynamic interaction with the surrounding environment. Using techniques such as genome and transcriptome sequencing, alongside protein prediction methodologies, we identified putative extracellular peptidases within Aspergillus ochraceus VKM-F4104D. Following manual annotation procedures, a total of 11 aspartic, 2 cysteine, 2 glutamic, 21 serine, 1 threonine, and 21 metallopeptidases were attributed to the extracellular degradome of A. ochraceus VKM-F4104D. Among them are enzymes with promising applications in biotechnology, potential targets and agents for antifungal therapy, and microbial antagonism factors. Thus, additional functionalities of the extracellular degradome, extending beyond mere protein substrate digestion for nutritional purposes, were demonstrated.

Discovery of Borosin Catalytic Strategies and Function through Bioinformatic Profiling

Borosins are ribosomally synthesized and post-translationally modified peptides (RiPPs) containing backbone α-N-methylations. These modifications confer favorable pharmacokinetic properties including increased membrane permeability and resistance to proteolytic degradation. Previous studies have biochemically and bioinformatically explored several borosins, revealing (1) numerous domain architectures and (2) diverse core regions lacking conserved sequence elements. Due to these characteristics, large-scale computational identification of borosin biosynthetic genes remains challenging and often requires additional, time-intensive manual inspection. This work builds upon previous findings and updates the genome-mining tool RODEO to automatically evaluate borosin biosynthetic gene clusters (BGCs) and identify putative precursor peptides. Using the new RODEO module, we provide an updated analysis of borosin BGCs identified in the NCBI database. From our data set, we bioinformatically predict and experimentally characterize a new fused borosin domain architecture, in which the modified natural product core is encoded N-terminal to the methyltransferase domain. Additionally, we demonstrate that a borosin precursor peptide is a native substrate of shewasin A, a reported aspartyl peptidase with no previously identified substrates. Shewasin A requires post-translational modification of the leader peptide for proteolytic maturation, a feature not previously observed in RiPPs. Overall, this work provides a user-friendly and open-access tool for the analysis of borosin BGCs and we demonstrate its utility to uncover additional biosynthetic strategies within the borosin class of RiPPs.

Extracting binding energies and binding modes from biomolecular simulations of fragment binding to endothiapepsin

Fragment-based drug discovery (FBDD) aims to discover a set of small binding fragments that may be subsequently linked together. Therefore, in-depth knowledge of the individual fragments' structural and energetic binding properties is essential. In addition to experimental techniques, the direct simulation of fragment binding by molecular dynamics (MD) simulations became popular to characterize fragment binding. However, former studies showed that long simulation times and high computational demands per fragment are needed, which limits applicability in FBDD. Here, we performed short, unbiased MD simulations of direct fragment binding to endothiapepsin, a well-characterized model system of pepsin-like aspartic proteases. To evaluate the strengths and limitations of short MD simulations for the structural and energetic characterization of fragment binding, we predicted the fragments' absolute free energies and binding poses based on the direct simulations of fragment binding and compared the predictions to experimental data. The predicted absolute free energies are in fair agreement with the experiment. Combining the MD data with binding mode predictions from molecular docking approaches helped to correctly identify the most promising fragments for further chemical optimization. Importantly, all computations and predictions were done within 5 days, suggesting that MD simulations may become a viable tool in FBDD projects.

Tuning Active Site Electron Density for Enhanced Molecular Recognition and Catalysis

Enzymes have an extraordinary ability to utilize aromatic interactions for molecular recognition and catalysis. We here report molecularly imprinted nanoparticle receptors. The aromatic "wall" material in the imprinted binding site is used to enhance the molecular recognition of aromatic guests that have similar charges, shapes, and sizes but differ in π-electron density. Additionally, aromatic interactions are employed to activate an electron-rich aryl leaving group on a glycoside, mimicking the nucleoside hydrolase of the parasite Trypanosoma vivax.

Aspartyl proteases identified as candidate genes of a fiber length QTL, qFLD05, that regulates fiber length in cotton (Gossypium hirsutum L.)

KEY MESSAGE

GhAP genes were identified as the candidates involved in cotton fiber length under the scope of fine mapping a stable fiber length QTL, qFLD05. Moreover, the transcription factor GhWRKY40 positively regulated GhAP3 to decrease fiber length. Fiber length (FL) is an economically important fiber quality trait. Although several genes controlling cotton fiber development have been identified, our understanding of this process remains limited. In this study, an FL QTL (qFLD05) was fine-mapped to a 216.9-kb interval using a secondary F2:3 population derived from the upland hybrid cultivar Ji1518. This mapped genomic segment included 15 coding genes, four of which were annotated as aspartyl proteases (GhAP1-GhAP4). GhAPs were identified as candidates for qFLD05 as the sequence variations in GhAPs were associated with FL deviations in the mapping population, and functional validation of GhAP3 and GhAP4 indicated a longer FL following decreases in their expression levels through virus-induced gene silencing (VIGS). Subsequently, the potential involvement of GhWRKY40 in the regulatory network was revealed: GhWRKY40 positively regulated GhAP3's expression according to transcriptional profiling, VIGS, yeast one-hybrid assays and dual-luciferase experiments. Furthermore, alterations in the expression of the eight previously reported cotton FL-responsive genes from the above three VIGS lines (GhAP3, GhAP4 and GhWRKY40) implied that MYB5_A12 was involved in the GhWRKY40-GhAP network. In short, we unveiled the unprecedented FL regulation roles of GhAPs in cotton, which was possibly further regulated by GhWRKY40. These findings will reveal the genetic basis of FL development associated with qFLD05 and be beneficial for the marker-assisted selection of long-staple cotton.

Advances in screening, synthesis, modification, and biomedical applications of peptides and peptide aptamers

Peptides and peptide aptamers have emerged as promising molecules for a wide range of biomedical applications due to their unique properties and versatile functionalities. The screening strategies for identifying peptides and peptide aptamers with desired properties are discussed, including high-throughput screening, display screening technology, and in silico design approaches. The synthesis methods for the efficient production of peptides and peptide aptamers, such as solid-phase peptide synthesis and biosynthesis technology, are described, along with their advantages and limitations. Moreover, various modification techniques are explored to enhance the stability, specificity, and pharmacokinetic properties of peptides and peptide aptamers. This includes chemical modifications, enzymatic modifications, biomodifications, genetic engineering modifications, and physical modifications. Furthermore, the review highlights the diverse biomedical applications of peptides and peptide aptamers, including targeted drug delivery, diagnostics, and therapeutic. This review provides valuable insights into the advancements in screening, synthesis, modification, and biomedical applications of peptides and peptide aptamers. A comprehensive understanding of these aspects will aid researchers in the development of novel peptide-based therapeutics and diagnostic tools for various biomedical challenges.

Lipid environment modulates processivity and kinetics of a presenilin homolog acting on multiple substrates in vitro

Intramembrane proteases (IPs) hydrolyze peptides in the lipid membrane. IPs participate in a number of cellular pathways including immune response and surveillance, and cholesterol biosynthesis, and they are exploited by viruses for replication. Despite their broad importance across biology, how activity is regulated in the cell to control protein maturation and release of specific bioactive peptides at the right place and right time remains largely unanswered, particularly for the intramembrane aspartyl protease (IAP) subtype. At a molecular biochemical level, different IAP homologs can cleave non-biological substrates, and there is no sequence recognition motif among the nearly 150 substrates identified for just one IAP, presenilin-1, the catalytic component of γ-secretase known for its involvement in the production of amyloid-β plaques associated with Alzheimer disease. Here we used gel-based assays combined with quantitative mass spectrometry and FRET-based kinetics assays to probe the cleavage profile of the presenilin homolog from the methanogen Methanoculleus marisnigri JR1 as a function of the surrounding lipid-mimicking environment, either detergent micelles or bicelles. We selected four biological IAP substrates that have not undergone extensive cleavage profiling previously, namely, the viral core protein of Hepatitis C virus, the viral core protein of Classical Swine Fever virus, the transmembrane segment of Notch-1, and the tyrosine receptor kinase ErbB4. Our study demonstrates a proclivity toward cleavage of substrates at positions of low average hydrophobicity and a consistent role for the lipid environment in modulating kinetic properties.

Charge Distribution Patterns of IA3 Impact Conformational Expansion and Hydration Diffusivity of the Disordered Ensemble

IA3 is a 68 amino acid natural peptide/protein inhibitor of yeast aspartic proteinase A (YPRA) that is intrinsically disordered in solution with induced N-terminal helicity when in the protein complex with YPRA. Based on the intrinsically disordered protein (IDP) parameters of fractional net charge (FNC), net charge density per residue (NCPR), and charge patterning (κ), the two domains of IA3 are defined to occupy different domains within conformationally based subclasses of IDPs, thus making IA3 a bimodal domain IDP. Site-directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy and low-field Overhauser dynamic nuclear polarization (ODNP) spectroscopy results show that these two domains possess different degrees of compaction and hydration diffusivity behavior. This work suggests that SDSL EPR line shapes, analyzed in terms of their local tumbling volume (VL), provide insights into the compaction of the unstructured IDP ensemble in solution and that protein sequence and net charge distribution patterns within a conformational subclass can impact bound water hydration dynamics, thus possibly offering an alternative thermodynamic property that can encode conformational binding and behavior of IDPs and liquid-liquid phase separations.

Hydrolysis of Gluten-Derived Celiac Disease-Triggering Peptides across a Broad pH Range by RmuAP1: A Novel Aspartic Peptidase Isolated from Rhodotorula mucilaginosa

An aspartate peptidase with proteolytic activity toward gluten was identified from an isolated red yeast Rhodotorula mucilaginosa strain. This peptidase consists of 425 amino acids, comprising an N-terminal signal peptide, a propeptide, and a C-terminal catalytic domain. The catalytic domain, termed RmuAP1CD, could be secreted by the recombinant oleaginous yeast Yarrowia lipolytica, whose genome contains the expression cassette for RmuAP1CD. RmuAP1CD exhibited optimum activity at pH 2.5 when acting on bovine serum albumin. Moreover, it facilitated the hydrolysis of gluten-derived immunogenic peptides (GIPs), which are responsible for triggering celiac disease symptoms, across a pH range of 3.0-6.0. The preferred cleavage sites are P-Q-Q-↓-P-Q in the 26-mer and P-Q-L-↓-P-Y in the 33-mer GIPs. Conversely, porcine pepsin cannot hydrolyze these two GIPs. The ability of RmuAP1CD to degrade GIPs under acidic conditions of the stomach indicates its potential as a viable oral enzyme therapy for celiac disease.

Genome-Wide Analysis of Haemonchus contortus Proteases and Protease Inhibitors Using Advanced Informatics Provides Insights into Parasite Biology and Host-Parasite Interactions

Biodiversity within the animal kingdom is associated with extensive molecular diversity. The expansion of genomic, transcriptomic and proteomic data sets for invertebrate groups and species with unique biological traits necessitates reliable in silico tools for the accurate identification and annotation of molecules and molecular groups. However, conventional tools are inadequate for lesser-known organismal groups, such as eukaryotic pathogens (parasites), so that improved approaches are urgently needed. Here, we established a combined sequence- and structure-based workflow system to harness well-curated publicly available data sets and resources to identify, classify and annotate proteases and protease inhibitors of a highly pathogenic parasitic roundworm (nematode) of global relevance, called Haemonchus contortus (barber's pole worm). This workflow performed markedly better than conventional, sequence-based classification and annotation alone and allowed the first genome-wide characterisation of protease and protease inhibitor genes and gene products in this worm. In total, we identified 790 genes encoding 860 proteases and protease inhibitors representing 83 gene families. The proteins inferred included 280 metallo-, 145 cysteine, 142 serine, 121 aspartic and 81 "mixed" proteases as well as 91 protease inhibitors, all of which had marked physicochemical diversity and inferred involvements in >400 biological processes or pathways. A detailed investigation revealed a remarkable expansion of some protease or inhibitor gene families, which are likely linked to parasitism (e.g., host-parasite interactions, immunomodulation and blood-feeding) and exhibit stage- or sex-specific transcription profiles. This investigation provides a solid foundation for detailed explorations of the structures and functions of proteases and protease inhibitors of H. contortus and related nematodes, and it could assist in the discovery of new drug or vaccine targets against infections or diseases.

Potential of Human Hemoglobin as a Source of Bioactive Peptides: Comparative Study of Enzymatic Hydrolysis with Bovine Hemoglobin and the Production of Active Peptide α137-141

Cruor, the main component responsible for the red color of mammalian blood, contains 90% haemoglobin, a protein considered to be a rich source of bioactive peptides. The aim of the present study is to assess the potential of human hemoglobin as a source of bioactive peptides, compared with bovine hemoglobin, which has been extensively studied in recent years. More specifically, the study focused on the α137-141 fragment of bovine haemoglobin (TSKYR), a small (653 Da) hydrophilic antimicrobial peptide. In this work, the potential of human hemoglobin to contain bioactive peptides was first investigated in silico in comparison with bovine hemoglobin-derived peptides using bioinformatics tools. The blast results showed a high identity, 88% and 85% respectively, indicating a high similarity between the α and β chains. Peptide Cutter software was used to predict cleavage sites during peptide hydrolysis, revealing major conservation in the number and location of cleavage sites between the two species, while highlighting some differences. Some peptides were conserved, notably our target peptide (TSKYR), while others were specific to each species. Secondly, the two types of hemoglobin were subjected to similar enzymatic hydrolysis conditions (23 °C, pH 3.5), which showed that the hydrolysis of human hemoglobin followed the same reaction mechanism as the hydrolysis of bovine hemoglobin, the 'zipper' mechanism. Concerning the peptide of interest, α137-141, the RP-UPLC analyses showed that its identification was not affected by the increase in the initial substrate concentration. Its production was rapid, with more than 60% of the total α137-141 peptide production achieved in just 30 min of hydrolysis, reaching peak production at 3 h. Furthermore, increasing the substrate concentration from 1% to 10% (w/v) resulted in a proportional increase in α137-141 production, with a maximum concentration reaching 687.98 ± 75.77 mg·L-1, approximately ten-fold higher than that obtained with a 1% (w/v) concentration. Finally, the results of the UPLC-MS/MS analysis revealed the identification of 217 unique peptides in bovine hemoglobin hydrolysate and 189 unique peptides in human hemoglobin hydrolysate. Of these, 57 peptides were strictly common to both species. This revealed the presence of several bioactive peptides in both cattle and humans. Although some had been known previously, new bioactive peptides were discovered in human hemoglobin, such as four antibacterial peptides (α37-46 PTTKTYFPHF, α36-45 FPTTKTYFPH, α137-141 TSKYR, and α133-141 STVLTSKYR), three opioid peptides (α137-141 TSKYR,β31-40 LVVYPWTQRF,β32-40, VVYPWTQRF), an ACE inhibitor (β129-135 KVVAGVA), an anticancer agent (β33-39 VVYPWTQ), and an antioxidant (α137-141 TSKYR). To the best of our knowledge, these peptides have never been found in human hemoglobin before.

Hydrogen Bonding Compensation on the Convex Solvent-Exposed Helical Face of IA3, an Intrinsically Disordered Protein

Saccharomyces cerevisiae IA₃ is a 68 amino acid peptide inhibitor of yeast proteinase A (YPRA) characterized as a random coil when in solution, folding into an N-terminal amphipathic alpha helix for residues 2-32 when bound to YPRA, with residues 33-68 unresolved in the crystal complex. Circular dichroism (CD) spectroscopy results show that amino acid substitutions that remove hydrogen-bonding interactions observed within the hydrophilic face of the N-terminal domain (NTD) of IA₃-YPRA crystal complex reduce the 2,2,2-trifluoroethanol (TFE)-induced helical transition in solution. Although nearly all substitutions decreased TFE-induced helicity compared to wild-type (WT), each construct did retain helical character in the presence of 30% (v/v) TFE and retained disorder in the absence of TFE. The NTDs of 8 different Saccharomyces species have nearly identical amino acid sequences, indicating that the NTD of IA₃ may be highly evolved to adopt a helical fold when bound to YPRA and in the presence of TFE but remain unstructured in solution. Only one natural amino acid substitution explored within the solvent-exposed face of the NTD of IA₃ induced TFE-helicity greater than the WT sequence. However, chemical modification of a cysteine by a nitroxide spin label that contains an acetamide side chain did enhance TFE-induced helicity. This finding suggests that non-natural amino acids that can increase hydrogen bonding or alter hydration through side-chain interactions may be important to consider when rationally designing intrinsically disordered proteins (IDPs) with varied biotechnological applications.

Overview of the Properties of Glutamic Peptidases That Are Present in Plant and Bacterial Pathogens and Play a Role in Celiac Disease and Cancer

Seven peptidase (proteinase) families─aspartic, cysteine, metallo, serine, glutamic, threonine, and asparagine─are in the peptidase database MEROPS, version 12.4 (https://www.ebi.ac.uk/merops/). The glutamic peptidase family is assigned two clans, GA and GB, and comprises six subfamilies. This perspective summarizes the unique features of their representatives. (1) G1, scytalidoglutamic peptidase, has a β-sandwich structure containing catalytic residues glutamic acid (E) and glutamine (Q), thus the name eqolisin. Most family members are pepstatin-insensitive and act as plant pathogens. (2) G2, preneck appendage protein, originates in phages, is a transmembrane protein, and its catalytic residues consist of glutamic and aspartic acids. (3) G3, strawberry mottle virus glutamic peptidase, originates in viruses and has a β-sandwich structure with catalytic residues E and Q. Neprosin has propyl endopeptidase activity, is associated with celiac disease, has a β-sandwich structure, and contains catalytic residues E-E and Q-tryptophan. (4) G4, Tiki peptidase, of the erythromycin esterase family, is a transmembrane protein, and its catalytic residues are E-histidine pairs. (5) G5, RCE1 peptidase, is associated with cancer, is a transmembrane protein, and its catalytic residues are E-histidine and asparagine-histidine. Microcystinase, a bacterial toxin, is a transmembrane protein with catalytic residues E-histidine and asparagine-histidine. (6) G6, Ras/Rap1-specific peptidase, is a bacterial pathogen, a transmembrane protein, and its catalytic residues are E-histidine pairs. This family's common features are that their catalytic residues consist of a glutamic acid and another (variable) amino acid and that they exhibit a diversity of biological functions─plant and bacterial pathogens and involvement in celiac disease and cancer─that suggests they are viable drug targets.

Golden 2-like transcription factor contributes to the major QTL against rice black-streaked dwarf virus disease

A major resistance QTL was identified on chromosome 6 in rice variety Wuke; both overexpression and knockdown experiments confirmed that OsGLK1 is the candidate gene for association with Rice black-streaked dwarf virus disease. Rice black-streaked dwarf virus disease is one of the most destructive rice viral diseases in China and East Asia. Progress has been limited in RBSDVD resistance breeding due to inadequate knowledge on the underlying functional genes. In this study, a major QTL for RBSDV (rice black-streaked dwarf virus) independent of SBPH (small brown planthopper) resistance was mapped in a 1.8 Mb interval on chromosome 6 by using an F_2:3 population originated from resistant rice variety Wuke. Representative transcripts within this region were analysed and three genes showing amino acid sequence variation in functional domains were selected for transformation. Overexpression experiments showed that one gene exhibited significant enhanced resistance compared to control lines, encoding protein involving Myb domain and probable transcription factor Golden 2-like1 (GLK1). Furthermore, OsGLK1 knockdown rice lines were investigated and the resistance ability was significantly declined without this gene compared to the wild type. Taken together, both overexpression and knockdown experiments strongly suggested that OsGLK1 plays an important role for RBSDV resistance and contributes to the major QTL. The study paves the way for elucidating the molecular mechanism underlying RBSDVD resistance and the molecular markers associated with OsGLK1 may be used for marker-assisted selection.

Serine-Carboxyl Peptidases, Sedolisins: From Discovery to Evolution

Sedolisin is a proteolytic enzyme, listed in the peptidase database MEROPS as a founding member of clan SB, family S53. This enzyme, although active at low pH, was originally shown not to be inhibited by an aspartic peptidase specific inhibitor, S-PI (pepstatin Ac). In this Perspective, the S53 family is described from the moment of original identification to evolution. The representative enzymes of the family are sedolisin, kumamolisin, and TPP-1. They exhibit the following unique features. (1) The fold of the molecule is similar to that of subtilisin, but the catalytic residues consist of a triad, Ser/Glu/Asp, that is unlike the Ser/His/Asp triad of subtilisin. (2) The molecule is expressed as a pro-form composed of the amino-terminal prosegment and the active domain. Additionally, some members of this family have an additional, carboxy-terminal prosegment. (3) Their optimum pH for activity is in the acidic region, not in the neutral to alkaline region where subtilisin is active. (4) Their distribution in nature is very broad across the three kingdoms of life. (5) Some of these enzymes from fungi and bacteria are pathogens to plants. (6) Some of them have significant potential applications for industry. (7) The lack of a TPP-1 gene in human brain is the cause of incurable juvenile neuronal ceroid lipofuscinosis (Batten's disease).

Transcriptional and post-transcriptional regulation of ethylene biosynthesis by exogenous acetylsalicylic acid in kiwifruit

Levels of ethylene, implicated in the induction of fruit ripening in a diverse array of plants, are influenced by genetic and environmental factors, such as other plant hormones. Among these, salicylic acid (SA) and its derivative, acetylsalicylic acid (ASA), have been demonstrated to inhibit ethylene biosynthesis in fruit, yet the underlying regulatory mechanisms remain elusive. Here, we showed that treatment with exogenous ASA dramatically reduced ethylene production, as well as activities of ACC synthase (ACS) and ACC oxidase (ACO), in kiwifruit tissues. Comparative transcriptome analysis indicated the differential expression of ethylene biosynthetic genes (AdACS1/2 and AdACO5). A screen of transcription factors indicated that AdERF105L and AdWRKY29 were ASA-responsive regulators of AdACS1/2 and AdACO5, respectively. In addition to these genes, AdACS3 and AdACO3 were abundantly expressed in both ASA-treated and control tissues. AdACS3 protein was phosphorylated and stabilized by AdMPK16, a mitogen-activated protein kinase, while AdACO3 activity was enhanced by AdAP, an aspartic peptidase. Exogenous ASA downregulated AdMPK16 and AdAP, thereby influencing ethylene biosynthesis at a post-transcriptional level. These findings led us to propose a multidimensional system for inhibition of ethylene biosynthesis by ASA, inducing differential expression of some ethylene biosynthesis genes, as well as differential effects on protein activity on other targets.

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.

Structures of plasmepsin X from Plasmodium falciparum reveal a novel inactivation mechanism of the zymogen and molecular basis for binding of inhibitors in mature enzyme

Plasmodium falciparum plasmepsin X (PfPMX), involved in the invasion and egress of this deadliest malarial parasite, is essential for its survival and hence considered as an important drug target. We report the first crystal structure of PfPMX zymogen containing a novel fold of its prosegment. A unique twisted loop from the prosegment and arginine 244 from the mature enzyme is involved in zymogen inactivation; such mechanism, not previously reported, might be common for apicomplexan proteases similar to PfPMX. The maturation of PfPMX zymogen occurs through cleavage of its prosegment at multiple sites. Our data provide thorough insights into the mode of binding of a substrate and a potent inhibitor 49c to PfPMX. We present molecular details of inactivation, maturation, and inhibition of PfPMX that should aid in the development of potent inhibitors against pepsin-like aspartic proteases from apicomplexan parasites.

An aspartic protease 47 causes quantitative recessive resistance to rice black-streaked dwarf virus disease and southern rice black-streaked dwarf virus disease

Rice black-streaked dwarf virus disease (RBSDVD) and southern rice black-streaked dwarf virus disease (SRBSDVD) are the most destructive viral diseases in rice. Progress is limited in breeding due to lack of resistance resource and inadequate knowledge on the underlying functional gene. Using genome-wide association study (GWAS), linkage disequilibrium (LD) decay analyses, RNA-sequencing, and genome editing, we identified a highly RBSDVD-resistant variety and its first functional gene. A highly RBSDVD-resistant variety W44 was identified through extensive evaluation of a diverse international rice panel. Seventeen quantitative trait loci (QTLs) were identified among which qRBSDV6-1 had the largest phenotypic effect. It was finely mapped to a 0.8-1.2 Mb region on chromosome 6, with 62 annotated genes. Analysis of the candidate genes underlying qRBSDV6-1 showed high expression of aspartic proteinase 47 (OsAP47) in a susceptible variety, W122, and a low resistance variety, W44. OsAP47 overexpressing lines exhibited significantly reduced resistance, while the knockout mutants exhibited significantly reduced SRBSDVD and RBSDVD severity. Furthermore, the resistant allele Hap1 of OsAP47 is almost exclusive to Indica, but rare in Japonica. Results suggest that OsAP47 knockout by editing is effective for improving RBSDVD and SRBSDVD resistance. This study provides genetic information for breeding resistant cultivars.

Flap Dynamics in Pepsin-Like Aspartic Proteases: A Computational Perspective Using Plasmepsin-II and BACE-1 as Model Systems

The flexibility of β hairpin structure known as the flap plays a key role in catalytic activity and substrate intake in pepsin-like aspartic proteases. Most of these enzymes share structural and sequential similarity. In this study, we have used apo Plm-II and BACE-1 as model systems. In the apo form of the proteases, a conserved tyrosine residue in the flap region remains in a dynamic equilibrium between the normal and flipped states through rotation of the χ1 and χ2 angles. Independent MD simulations of Plm-II and BACE-1 remained stuck either in the normal or flipped state. Metadynamics simulations using side-chain torsion angles (χ1 and χ2 of tyrosine) as collective variables sampled the transition between the normal and flipped states. Qualitatively, the two states were predicted to be equally populated. The normal and flipped states were stabilized by H-bond interactions to a tryptophan residue and to the catalytic aspartate, respectively. Further, mutation of tyrosine to an amino-acid with smaller side-chain, such as alanine, reduced the flexibility of the flap and resulted in a flap collapse (flap loses flexibility and remains stuck in a particular state). This is in accordance with previous experimental studies, which showed that mutation to alanine resulted in loss of activity in pepsin-like aspartic proteases. Our results suggest that the ring flipping associated with the tyrosine side-chain is the key order parameter that governs flap dynamics and opening of the binding pocket in most pepsin-like aspartic proteases.

Genome-Wide Analyses of Aspartic Proteases on Potato Genome (Solanum tuberosum): Generating New Tools to Improve the Resistance of Plants to Abiotic Stress

Aspartic proteases are proteolytic enzymes widely distributed in living organisms and viruses. Although they have been extensively studied in many plant species, they are poorly described in potatoes. The present study aimed to identify and characterize S. tuberosum aspartic proteases. Gene structure, chromosome and protein domain organization, phylogeny, and subcellular predicted localization were analyzed and integrated with RNAseq data from different tissues, organs, and conditions focused on abiotic stress. Sixty-two aspartic protease genes were retrieved from the potato genome, distributed in 12 chromosomes. A high number of intronless genes and segmental and tandem duplications were detected. Phylogenetic analysis revealed eight StAP groups, named from StAPI to StAPVIII, that were differentiated into typical (StAPI), nucellin-like (StAPIIIa), and atypical aspartic proteases (StAPII, StAPIIIb to StAPVIII). RNAseq data analyses showed that gene expression was consistent with the presence of cis-acting regulatory elements on StAP promoter regions related to water deficit. The study presents the first identification and characterization of 62 aspartic protease genes and proteins on the potato genome and provides the baseline material for functional gene determinations and potato breeding programs, including gene editing mediated by CRISPR.

Rapid and Simple Analysis of the Human Pepsin Secondary Structure Using a Portable Raman Spectrometer

Human pepsin is a digestive protease that plays an important role in the human digestive system. The secondary structure of human pepsin determines its bioactivity. Therefore, an in-depth understanding of human pepsin secondary structure changes is particularly important for the further improvement of the efficiency of human pepsin biological function. However, the complexity and diversity of the human pepsin secondary structure make its analysis difficult. Herein, a convenient method has been developed to quickly detect the secondary structure of human pepsin using a portable Raman spectrometer. According to the change of surface-enhanced Raman spectroscopy (SERS) signal intensity and activity of human pepsin at different pH values, we analyze the change of the human pepsin secondary structure. The results show that the content of the β-sheet gradually increased with the increase in the pH in the active range, which is in good agreement with circular dichroism (CD) measurements. The change of the secondary structure improves the sensitivity of human pepsin SERS detection. Meanwhile, human pepsin is a commonly used disease marker for the noninvasive diagnosis of gastroesophageal reflux disease (GERD); the detection limit of human pepsin we obtained is 2 μg/mL by the abovementioned method. The real clinical detection scenario is also simulated by spiking pepsin solution in saliva, and the standard recovery rate is 80.7-92.3%. These results show the great prospect of our method in studying the protein secondary structure and furthermore promote the application of SERS in clinical diagnosis.

Effect of simulated in vitro upper gut digestion of processed cowpea beans on phenolic composition, antioxidant properties and cellular protection

The effect of simulated in vitro upper gut digestion on the phenolic composition and antioxidant properties of processed cowpea beans was studied. The samples comprised four cowpea cultivars: a cream, brownish-cream and two reddish-brown cultivars. Dry cowpea seeds were soaked in water, blended into paste and deep-fried in vegetable oil. The fried samples were taken through in vitro upper gut digestion followed by freeze-drying of the supernatant. Phenolic composition of extracts from the supernatants were determined using HPLC-MS. Radical scavenging activities were documented using the TEAC, ORAC and nitric oxide (NO) assays. In vitro digestion of the processed cowpeas resulted in phenolic-peptide complexes that were identified for the first time, and decreased extractable phenolic compounds. However, the radical scavenging activities increased. The processed cowpeas and their digests inhibited cellular NO production, and oxidative DNA and cellular damage. In conclusion, deep-fried cowpeas when consumed, could potentially help alleviate oxidative stress-related conditions.

Negatively charged phospholipids accelerate the membrane fusion activity of the plant-specific insert domain of an aspartic protease

Various plants use antimicrobial proteins/peptides to resist phytopathogens. In the potato, Solanum tuberosum, the plant-specific insert (PSI) domain of an aspartic protease performs this role by disrupting phytopathogen plasma membranes. However, the mechanism by which PSI selects target membranes has not been elucidated. Here, we studied PSI-induced membrane fusion, focusing on the effects of lipid composition on fusion efficiency. Membrane fusion by the PSI involves an intermediate state whereby adjacent liposomes share their bilayers. We found that increasing the concentration of negatively charged phosphatidylserine (PS) phospholipids substantially accelerated PSI-mediated membrane fusion. NMR data demonstrated that PS did not affect the binding between the PSI and liposomes but had seminal effects on the dynamics of PSI interaction with liposomes. In PS-free liposomes, the PSI underwent significant motion, which was suppressed on PS-contained liposomes. Molecular dynamics simulations showed that the PSI binds to PS-containing membranes with a dominant angle ranging from −31° to 30°, with respect to the bilayer, and is closer to the membrane surfaces. In contrast, PSI is mobile and exhibits multiple topological states on the surface of PS-free membranes. Taken together, our data suggested that PS lipids limit the motion of the anchored PSI, bringing it closer to the membrane surface and efficiently bridging different liposomes to accelerate fusion. As most phytopathogens have a higher content of negatively charged lipids as compared with host cells, these results indicate that the PSI selectively targets negatively charged lipids, which likely represents a way of distinguishing the pathogen from the host.

Toxoplasma gondii UBL-UBA shuttle proteins regulate several important cellular processes

Toxoplasma gondii is an obligate intracellular apicomplexan parasite causing lethal diseases in immunocompromised patients. UBL-UBA shuttle proteins (DDI1, RAD23, and DSK2) are important components of the ubiquitin-proteasome system. By degrading ubiquitinated proteins, UBL-UBA shuttle proteins regulate many cellular processes. However, the specific processes regulated by UBL-UBA shuttle proteins remain elusive. Here, we revealed that the deletion of shuttle proteins results in a selective accumulation of ubiquitinated proteins in the nucleus and aberrant DNA replication. ROP18 was mistargeted and accumulated in the shuttle protein mutant strain, resulting in the recruitment of immunity-related GTPases to the parasitophorous vacuole membrane (PVM). Furthermore, the mistargeting of ROP18 and the recruitment of Irgb6 to the PVM were also observed in the DDI1 mutant strain. DDI1 is a nonclassical UBL-UBA shuttle protein homologous to the HIV-1 protease. Molecular docking showed that DDI1 was a potential target of HIV-1 protease inhibitors. However, these inhibitors blocked the growth of T gondii in vitro but not in vivo. In conclusion, the Toxoplasma UBL-UBA shuttle protein regulates several important cellular processes and the mistargeting of ROP18 may be a representative of the abnormal homeostasis caused by shuttle protein mutation.

Improving the alkaline stability of pepsin through rational protein design using renin, an alkaline-stable aspartic protease, as a structural and functional reference

The present study sought to identify the structural determinants of aspartic protease structural stability and activity at elevated pH. Various hypotheses have been published regarding the features responsible for the unusual alkaline structural stability of renin, however, few structure-function studies have verified these claims. Using pepsin as a model system, and renin as a template for functional and structural alkaline stability, a rational re-design of pepsin was undertaken to identify residues contributing to the alkaline instability of pepsin-like aspartic proteases in regards to both structure and function. We constructed 13 mutants based on this strategy. Among them, mutants D159 L and D60A led to an increase in activity at elevated pH levels (p ≤ 0.05) and E4V and H53F were shown to retain native-like structure at elevated pH (p ≤ 0.05). Previously suggested carboxyl groups Asp¹¹, Asp¹¹⁸, and Glu¹³ were individually shown not to be responsible for the structural instability or lack of activity at neutral pH in pepsin. The importance of the β-barrel to structural stability was highlighted as the majority of the stabilizing residues identified, and 39% of the weakly conserved residues in the N-terminal lobe, were located in β-sheet strands of the barrel. The results of the present study indicate that alkaline stabilization of pepsin will require reduction of electrostatic repulsions and an improved understanding of the role of the hydrogen bonding network of the characteristic β-barrel.

Ribosomal Synthesis of Macrocyclic Peptides with Linear γ4- and β-Hydroxy-γ4-amino Acids

Herein, we report the ribosomal elongation of linear γ⁴- and β-hydroxy-γ⁴-amino acids (statines) to expand the nonproteinogenic amino acid repertoire of natural product-like combinatorial peptide libraries. First, we demonstrated the successful ribosomal incorporation of four linear γ⁴-amino acids (γ⁴Gly, (S)-γ⁴Ala, (S)-γ⁴Nva, and (R)-γ⁴Leu) into a 10-mer macrocyclic peptide scaffold. Given the promising effects reported for statines on the cell permeability of macrocyclic peptides, we also designed and tested the ribosomal incorporation of six statines derived from Ala and d-val. Four Ala-derived statines were successfully incorporated into peptides, and γ^4SAla^3R-OH (GP2) showed a similar efficiency of incorporation to that of (S)-β²hAla and l-Ala. These new building blocks might confer the important pharmacological properties of protease resistance and membrane permeability to macrocyclic peptides and expand the diversity of future combinatorial peptide libraries that can be translated by the ribosome.

Structural determinants for subnanomolar inhibition of the secreted aspartic protease Sapp1p from Candida parapsilosis

Pathogenic Candida albicans yeasts frequently cause infections in hospitals. Antifungal drugs lose effectiveness due to other Candida species and resistance. New medications are thus required. Secreted aspartic protease of C. parapsilosis (Sapp1p) is a promising target. We have thus solved the crystal structures of Sapp1p complexed to four peptidomimetic inhibitors. Three potent inhibitors (K_i: 0.1, 0.4, 6.6 nM) resembled pepstatin A (K_i: 0.3 nM), a general aspartic protease inhibitor, in terms of their interactions with Sapp1p. However, the weaker inhibitor (K_i: 14.6 nM) formed fewer nonpolar contacts with Sapp1p, similarly to the smaller HIV protease inhibitor ritonavir (K_i: 1.9 µM), which, moreover, formed fewer H-bonds. The analyses have revealed the structural determinants of the subnanomolar inhibition of C. parapsilosis aspartic protease. Because of the high similarity between Saps from different Candida species, these results can further be used for the design of potent and specific Sap inhibitor-based antimycotic drugs.

Re-emerging Aspartic Protease Targets: Examining Cryptococcus neoformans Major Aspartyl Peptidase 1 as a Target for Antifungal Drug Discovery

Cryptococcosis is an invasive infection that accounts for 15% of AIDS-related fatalities. Still, treating cryptococcosis remains a significant challenge due to the poor availability of effective antifungal therapies and emergence of drug resistance. Interestingly, protease inhibitor components of antiretroviral therapy regimens have shown some clinical benefits in these opportunistic infections. We investigated Major aspartyl peptidase 1 (May1), a secreted Cryptococcus neoformans protease, as a possible target for the development of drugs that act against both fungal and retroviral aspartyl proteases. Here, we describe the biochemical characterization of May1, present its high-resolution X-ray structure, and provide its substrate specificity analysis. Through combinatorial screening of 11,520 compounds, we identified a potent inhibitor of May1 and HIV protease. This dual-specificity inhibitor exhibits antifungal activity in yeast culture, low cytotoxicity, and low off-target activity against host proteases and could thus serve as a lead compound for further development of May1 and HIV protease inhibitors.

Introduction of a Glycine Linker Connecting the Heavy and Light Chains in Synthetic Cardosin B-Derived Rennet Changes the Specificity of Subpocket S3'

The development of plant-based synthetic rennets is of high commercial interest, due to the current great consumer demand for animal product alternatives. A previously developed recombinant form of the aspartic protease cardosin B with a three-glycine linker showed great potential due to its good performance in milk coagulation. This enzyme was found to be more specific and less proteolytically active than the native form for milk clotting, but the underlying structural causes for these activity changes were not completely clear. Here, we have performed molecular dynamics simulations with the recombinant enzyme with and without the linker. Our results showed that the introduction of the linker changes the subpocket S3', which is located more than 4 nm away. These results showcase how small modifications in proteins can have significant effects in distant regions in the protein structure that affect their biotechnological applications.

Selective Hydrolysis of Aryl Esters under Acidic and Neutral Conditions by a Synthetic Aspartic Protease Mimic

Aspartic proteases use a pair of carboxylic acids to activate water molecules for nucleophilic attack. Here we report a nanoparticle catalyst with a similar catalytic motif capable of generating a hydroxide ion in its active site even under acidic reaction conditions. The synthetic enzyme accelerated the hydrolysis of para-nitrophenyl acetate (PNPA) by 91,000 times and could also hydrolyze nonactivated aryl esters at pH 7. The distance between the two acids and, in particular, the flexibility of the catalytic groups in the active site controlled the catalytic efficiency. The synthetic enzyme readily detected the addition of a single methyl on the acyl group of the substrate, as well as the substitution pattern on the phenyl ring.

Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies

Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.

In vivo tumor growth inhibition by Solanum tuberosum aspartic protease 3 (StAP3) treatment

Solanum tuberosum aspartic Proteases (StAPs) show selective plasma membrane permeabilization, inducing cytotoxicity of cancer cells versus normal cells in vitro. Herein, we aimed to evaluate both StAP3 systemic toxicity and antitumoral activity against human melanoma in vivo. The toxicity of a single high dose of StAP3 (10 µg/g body weight, intraperitoneally) was assessed in a Balb/c mice model. Subcutaneous A375 human melanoma xenografts in athymic nude (nu/nu) mice were induced. Once tumors developed (mean larger dimension = 3.8 ± 0.09 mm), mice were StAP3-treated (6 µg/g body weight, subcutaneously under the tumor at a single dose). For both models, controls were treated with physiologic saline solution. StAP3-treated mice showed a significant inhibition of tumor growth (p < 0.05) compared with controls. No signs of toxicity were detected in StAP3-treated mice in both models. These results suggest the potential of these plant proteases as anticancer agents.

Pepsin-like aspartic proteases (PAPs) as model systems for combining biomolecular simulation with biophysical experiments

Pepsin-like aspartic proteases (PAPs) are a class of aspartic proteases which shares tremendous structural similarity with human pepsin. One of the key structural features of PAPs is the presence of a β-hairpin motif otherwise known as flap. The biological function of the PAPs is highly dependent on the conformational dynamics of the flap region. In apo PAPs, the conformational dynamics of the flap is dominated by the rotational degrees of freedom associated with χ1 and χ2 angles of conserved Tyr (or Phe in some cases). However it is plausible that dihedral order parameters associated with several other residues might play crucial roles in the conformational dynamics of apo PAPs. Due to their size, complexities associated with conformational dynamics and clinical significance (drug targets for malaria, Alzheimer's disease etc.), PAPs provide a challenging testing ground for computational and experimental methods focusing on understanding conformational dynamics and molecular recognition in biomolecules. The opening of the flap region is necessary to accommodate substrate/ligand in the active site of the PAPs. The BIG challenge is to gain atomistic details into how reversible ligand binding/unbinding (molecular recognition) affects the conformational dynamics. Recent reports of kinetics (K i, K d) and thermodynamic parameters (ΔH, TΔS, and ΔG) associated with macro-cyclic ligands bound to BACE1 (belongs to PAP family) provide a perfect challenge (how to deal with big ligands with multiple torsional angles and select optimum order parameters to study reversible ligand binding/unbinding) for computational methods to predict binding free energies and kinetics beyond typical test systems e.g. benzamide-trypsin. In this work, i reviewed several order parameters which were proposed to capture the conformational dynamics and molecular recognition in PAPs. I further highlighted how machine learning methods can be used as order parameters in the context of PAPs. I then proposed some open ideas and challenges in the context of molecular simulation and put forward my case on how biophysical experiments e.g. NMR, time-resolved FRET etc. can be used in conjunction with biomolecular simulation to gain complete atomistic insights into the conformational dynamics of PAPs.

Degradation of a Main Plastic Pollutant Polyethylene Terephthalate by Two Distinct Proteases (Neprilysin and Cutinase-like Enzyme)

In this DFT study, hydrolysis of polyethylene terephthalate (PET), a major cause of plastic pollution, by two distinct enzymes, neprilysin (NEP, a mononuclear metalloprotease) and cutinase-like enzyme (CLE, a serine protease), has been investigated. These enzymes utilize different mechanisms for the degradation of PET. NEP uses either the metal-bound hydroxide attack (MH) mechanism or reverse protonation (RP) mechanism, while CLE utilizes a general acid/base mechanism that includes acylation and deacylation processes. Additionally, the RP mechanism of NEP can proceed through three pathways, RP0, RP1, and RP2. The DFT calculations predict that, among all these mechanisms, the MH mechanism is the energetically most favorable one for the NEP enzyme. In comparison, CLE catalyzes this reaction with a significantly higher barrier. These results suggest that the Lewis acid and nucleophile activations provided by the Zn metal center of NEP are more effective than the hydrogen bonding interactions afforded by the catalytic Ser85-His180-Asp165 triad of CLE. They have provided intrinsic details regarding PET degradation and will pave the way for the design of efficient metal-based catalysts for this critical reaction.

Defense and Offense Strategies: The Role of Aspartic Proteases in Plant-Pathogen Interactions

Plant aspartic proteases (APs; E.C.3.4.23) are a group of proteolytic enzymes widely distributed among different species characterized by the conserved sequence Asp-Gly-Thr at the active site. With a broad spectrum of biological roles, plant APs are suggested to undergo functional specialization and to be crucial in developmental processes, such as in both biotic and abiotic stress responses. Over the last decade, an increasing number of publications highlighted the APs' involvement in plant defense responses against a diversity of stresses. In contrast, few studies regarding pathogen-secreted APs and AP inhibitors have been published so far. In this review, we provide a comprehensive picture of aspartic proteases from plant and pathogenic origins, focusing on their relevance and participation in defense and offense strategies in plant-pathogen interactions.

Tunable Artificial Enzyme-Cofactor Complex for Selective Hydrolysis of Acetals

Enzymes frequently use unimpressive functional groups such as weak carboxylic acids for efficient, highly selective catalysis including hydrolysis of acetals and even amides. Much stronger acids generally have to be used for such purposes in synthetic systems. We report here a method to position an acidic group near the acetal oxygen of 2-(4-nitrophenyl)-1,3-dioxolane bound by an artificial enzyme. The hydrolytic activity of the resulting artificial enzyme-cofactor complex was tuned by the number and depth of the active site as well as the hydrophobicity and acidity of the cofactor. The selectivity of the complex was controlled by the size and shape of the active site and enabled less reactive acetals to be hydrolyzed over more reactive ones.

Genome-wide identification, evolution and expression analysis of the aspartic protease gene family during rapid growth of moso bamboo (Phyllostachys edulis) shoots

BACKGROUND

Aspartic proteases (APs) are a class of aspartic peptidases belonging to nine proteolytic enzyme families whose members are widely distributed in biological organisms. APs play essential functions during plant development and environmental adaptation. However, there are few reports about APs in fast-growing moso bamboo.

RESULT

In this study, we identified a total of 129 AP proteins (PhAPs) encoded by the moso bamboo genome. Phylogenetic and gene structure analyses showed that these 129 PhAPs could be divided into three categories (categories A, B and C). The PhAP gene family in moso bamboo may have undergone gene expansion, especially the members of categories A and B, although homologs of some members in category C have been lost. The chromosomal location of PhAPs suggested that segmental and tandem duplication events were critical for PhAP gene expansion. Promoter analysis revealed that PhAPs in moso bamboo may be involved in plant development and responses to environmental stress. Furthermore, PhAPs showed tissue-specific expression patterns and may play important roles in rapid growth, including programmed cell death, cell division and elongation, by integrating environmental signals such as light and gibberellin signals.

CONCLUSION

Comprehensive analysis of the AP gene family in moso bamboo suggests that PhAPs have experienced gene expansion that is distinct from that in rice and may play an important role in moso bamboo organ development and rapid growth. Our results provide a direction and lay a foundation for further analysis of plant AP genes to clarify their function during rapid growth.

Molecular characterization of a novel aspartyl protease-1 from Trichinella spiralis

The aim of this study was to characterize the biological properties of a novel aspartic protease-1 from Trichinella spiralis (TsASP1) and evaluate its potential in inducing immune response. TsASP1 gene was cloned and expressed in Escherichia coli BL21 (DE3). On Western blotting analysis with anti-rTsASP1 serum, native TsASP1 was detected in various T. spiralis phases other than newborn larvae (NBL). qPCR results showed that TsASP1 transcription was the highest in intestinal infective larvae (IIL) and the lowest in the NBL stage. Immunofluorescence test result shows that native TsASP1 was principally localized in stichosome, muscle cells of muscle larvae (ML) and IIL, and surrounded intrauterine embryos in female adult worms (AW). After silencing TsASP1 gene of the ML by siRNA, the worm development was significantly inhibited, showed by shorter AW and more wrinkles and longitudinal crack on epicuticle of AW on scanning electron microscopy; the AW and ML burdens were reduced by 41.82 and 56.36% respectively, compared with the control siRNA or PBS group (P < 0.001). Immunization of mice with rTsASP1 elicited an evident antibody response (serum IgG, IgG1/IgG2a and enteral sIgA), and systemic (spleen) and intestinal local mucosal (mesenteric lymph node) cellular immune response, demonstrated by a prominent elevation of IFN-γ and IL-4. The results suggested TsASP1 participated in T. spiralis development and survival in host, and immunization of mice with rTsASP1 induced systemic/intestinal local mucosal humoral and cellular immune response against Trichinella.

The crystal structure of the tetrameric human vasohibin-1-SVBP complex reveals a variable arm region within the structural core

Vasohibins regulate angiogenesis, tumor growth, metastasis and neuronal differentiation. They form a complex with small vasohibin-binding protein (SVBP) and show tubulin tyrosine carboxypeptidase activity. Recent crystal structure determinations of vasohibin-SVBP complexes have provided a molecular basis for complex formation, substrate binding and catalytic activity. However, the regulatory mechanism and dynamics of the complex remain elusive. Here, the crystal structure of the VASH1-SVBP complex and a molecular-dynamics simulation study are reported. The overall structure of the complex was similar to previously reported structures. Importantly, however, the structure revealed a domain-swapped heterotetramer that was formed between twofold symmetry-related molecules. This heterotetramerization was stabilized by the mutual exchange of ten conserved N-terminal residues from the VASH1 structural core, which was intramolecular in other structures. Interestingly, a comparison of this region with previously reported structures revealed that the patterns of hydrogen bonding and hydrophobic interactions vary. In the molecular-dynamics simulations, differences were found between the heterotetramer and heterodimer, where the fluctuation of the N-terminal region in the heterotetramer was suppressed. Thus, heterotetramer formation and flexibility of the N-terminal region may be important for enzyme activity and regulation.

Characterization of the cathepsin D in Procambarus clarkii and its biological role in innate immune responses

Cathepsin D belongs to aspartic protease family, produced in the rough endoplasmic reticulum, and then transported to lysosomes, where it participates in various physiological processes. Despite its importance, only a few reports available on the functional role of cathepsin D in crustaceans. Herein, we cloned a cDNA fragment of cathepsin D from the hepatopancreas of the red swamp crayfish, Procambarus clarkii (Pc-cathepsin D) for the first time. It included 1158 base pairs open reading frame, encoding a protein of 385 amino acids. Multiple alignment analysis confirmed the presence of aspartic proteinase active sites and N glycosylation sites. Pc-cathepsin D mRNA expression was high in the gills followed by gut, heart, hepatopancreas of P. clarkii. At different time points post-infection with lipopolysaccharides, peptidoglycan, or polyinosinic polycytidylic acid, Pc-cathepsin D mRNA expression significantly enhanced compared with the control group. Knockdown of the Pc-cathepsin D by double-stranded RNA, strikingly, changed the expression of all the tested P. clarkii immune-associated genes, including Pc-Toll, Pc-lectin, Pc-cactus, Pc-anti-lipopolysaccharide factor, Pc-phospholipase, and Pc-sptzale. Altogether, these results suggest that Pc-cathepsin D is needed to confer innate immunity against microbial pathogens by modulating the expression of crucial transcripts that encode immune-associated genes.