Analysis of crystal structures of aspartic proteinases: on the role of amino acid residues adjacent to the catalytic site of pepsin-like enzymes

Toxoplasma gondii aspartic protease 5: structural basis of substrate binding and inhibition mechanism

Toxoplasma gondii, a worldwide prevalent parasite is responsible for causing toxoplasmosis in almost all warm-blooded animals, including humans. Golgi-resident T. gondii aspartic protease 5 (TgASP5) plays an essential role in the maturation and export of the effector proteins those modulate the host immune system to establish a successful infection. Hence, inhibiting this enzyme can be a possible therapeutic strategy against toxoplasmosis. This is the first report of the detailed structural investigations of the TgASP5 mature enzyme using molecular modeling and an all-atom simulation approach which provide in-depth knowledge of the active site architecture of TgASP5. The analysis of the binding mode of the TEXEL (Toxoplasma EXport Element) substrate to TgASP5 highlighted the importance of the active site residues. Ser505, Ala776 and Tyr689 in the S2 binding pocket are responsible for the specificity towards Arg at the P2 position of TEXEL substrate. The molecular basis of inhibition by the only known inhibitor RRLStatine has been identified, and our results show that it blocks the active site by forming a hydrogen bond with a catalytic aspartate. Besides that, known aspartic protease inhibitors were screened against TgASP5 using docking, MD simulations and MM-PBSA binding energy calculations. The top-ranked inhibitors (SC6, ZY1, QBH) showed higher binding energy than RRLStatine. Understanding the structural basis of substrate recognition and the binding mode of these inhibitors will help to develop potent mechanistic inhibitors against TgASP5. This study will also provide insights into the structural features of pepsin-like aspartic proteases from other apicomplexan parasites for developing antiparasitic agents.Communicated by Ramaswamy H. Sarma.

Assessing protein digestibility in allergenicity risk assessment: A comparison of in silico and high throughput in vitro gastric digestion assays

The susceptibility of a novel food protein to digestion in the pepsin resistance test is widely used to inform the allergenicity risk assessment process. However, it does not model the variation in the intragastric environment found in vivo. Consequently a 96-well plate format in vitro gastric digestion protocol has been developed with a high and low pepsin activity test executed at pH 1.2, 2.5, 5.5 and 6.5. It was used to analyse seven allergens (from milk, egg, peach and peanut) and two non-allergens (cytochrome c and zein). Digestion was monitored using SDS-PAGE and densitometry. In silico predictions were not confirmed experimentally for most of the proteins studied. Proteins were ranked according to half-life and showed susceptibility to digestion was related to the stability of protein structure and protein solubility rather than allergenicity per se. Highly digestible proteins, such as β-casein and Ara h 1, generated abundant resistant fragments Mr > 3.5 kDa in the low pepsin activity test which could be immunologically significant within the context of allergenicity risk assessment for susceptible groups such as infants. The high- and low pepsin activity tests used in this study provided complementary data to support allergenicity risk assessment and used only 10 mg protein.

Exploring Aspartic Protease Inhibitor Binding to Design-Selective Antimalarials

Selectivity is a major issue in the development of drugs targeting pathogen aspartic proteases. Here, we explore the selectivity-determining factors by studying specifically designed malaria aspartic protease (plasmepsin) open-flap inhibitors. Metadynamics simulations are used to uncover the complex binding/unbinding pathways of these inhibitors and describe the critical transition states in atomistic resolution. The simulation results are compared with experimentally determined enzymatic activities. Our findings demonstrate that plasmepsin inhibitor selectivity can be achieved by targeting the flap loop with hydrophobic substituents that enable ligand binding under the flap loop, as such a behavior is not observed for several other aspartic proteases. The ability to estimate the selectivity of compounds before they are synthesized is of considerable importance in drug design; therefore, we expect that our approach will be useful in selective inhibitor designs against not only aspartic proteases but also other enzyme classes.

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.

Exploring the pH- and Ligand-Dependent Flap Dynamics of Malarial Plasmepsin II

Malaria remains a global health threat─over 400,000 deaths occurred in 2019. Plasmepsins are promising targets of antimalarial therapeutics; however, no inhibitors have reached the clinic. To fuel the progress, a detailed understanding of the pH- and ligand-dependent conformational dynamics of plasmepsins is needed. Here we present the continuous constant pH molecular dynamics study of the prototypical plasmepsin II and its complexed form with a substrate analogue. The simulations revealed that the catalytic dyads D34 and D214 are highly coupled in the apo protein and that the pepstatin binding enhances the difference in proton affinity, making D34 the general base and D214 the general acid. The simulations showed that the flap adopts an open state regardless of pH; however, upon pepstatin binding the flap can close or open depending on the protonation state of D214. These and other data are discussed and compared with the off-targets human cathepsin D and renin. This study lays the groundwork for a systematic investigation of pH- and ligand-modulated dynamics of the entire family of plasmepsins to help design more potent and selective 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.

Exploring the pH-Dependent Structure-Dynamics-Function Relationship of Human Renin

Renin is a pepsin-like aspartyl protease and an important drug target for the treatment of hypertension; despite three decades' research, its pH-dependent structure-function relationship remains poorly understood. Here, we employed continuous constant pH molecular dynamics (CpHMD) simulations to decipher the acid/base roles of renin's catalytic dyad and the conformational dynamics of the flap, which is a common structural feature among aspartyl proteases. The calculated pK_a's suggest that catalytic Asp38 and Asp226 serve as the general base and acid, respectively, in agreement with experiment and supporting the hypothesis that renin's neutral optimum pH is due to the substrate-induced pK_a shifts of the aspartic dyad. The CpHMD data confirmed our previous hypothesis that hydrogen bond formation is the major determinant of the dyad pK_a order. Additionally, our simulations showed that renin's flap remains open regardless of pH, although a Tyr-inhibited state is occasionally formed above pH 5. These findings are discussed in comparison to the related aspartyl proteases, including β-secretases 1 and 2, cathepsin D, and plasmepsin II. Our work represents a first step toward a systematic understanding of the pH-dependent structure-dynamics-function relationships of pepsin-like aspartyl proteases that play important roles in biology and human disease states.

Structural basis for the specificity of renin-mediated angiotensinogen cleavage

The renin–angiotensin cascade is a hormone system that regulates blood pressure and fluid balance. Renin-mediated cleavage of the angiotensin I peptide from the N terminus of angiotensinogen (AGT) is the rate-limiting step of this cascade; however, the detailed molecular mechanism underlying this step is unclear. Here, we solved the crystal structures of glycosylated human AGT (2.30 Å resolution), its encounter complex with renin (2.55 Å), AGT cleaved in its reactive center loop (RCL; 2.97 Å), and spent AGT from which the N-terminal angiotensin peptide was removed (2.63 Å). These structures revealed that AGT undergoes profound conformational changes and binds renin through a tail-into-mouth allosteric mechanism that inserts the N terminus into a pocket equivalent to a hormone-binding site on other serpins. These changes fully extended the N-terminal tail, with the scissile bond for angiotensin release docked in renin's active site. Insertion of the N terminus into this pocket accompanied a complete unwinding of helix H of AGT, which, in turn, formed key interactions with renin in the complementary binding interface. Mutagenesis and kinetic analyses confirmed that renin-mediated production of angiotensin I is controlled by interactions of amino acid residues and glycan components outside renin's active-site cleft. Our findings indicate that AGT adapts unique serpin features for hormone delivery and binds renin through concerted movements in the N-terminal tail and in its main body to modulate angiotensin release. These insights provide a structural basis for the development of agents that attenuate angiotensin release by targeting AGT's hormone binding pocket.

Primary structures of different isoforms of buffalo pregnancy-associated glycoproteins (BuPAGs) during early pregnancy and elucidation of the 3-dimensional structure of the most abundant isoform BuPAG 7

Pregnancy-associated glycoproteins (PAGs) are expressed during pregnancy by the trophoectodermal cells of fetus. Presence of PAGs in dam's circulation has been widely used in pregnancy diagnosis. The present study reports the identification and characterization of different PAG isoforms in buffalo during early stages of pregnancy. The PAG mRNAs isolated from fetal cotyledons (Pregnancy stages: 45, 75 and 90 days) were successfully cloned in pJET1.2 vector and transformed in E. coli. A total of 360 random clones were sequenced and correlated with their stages of expression. A total of 12 isoforms namely, BuPAG 1, 2, 4, 6, 7, 8, 9, 13, 15, 16, 18 and one new isoform were identified. BuPAG 7 was found as the most abundant isoform in all three stages followed by BuPAG 18. Further, a large number of variants were found for most of these isoforms. Phylogenetic relationship of identified BuPAGs showed that BuPAG 2 belonged to an ancient group while other members clustered with modern group. Three-dimensional (3D) structure of BuPAG 7 was determined by homology modeling and molecular dynamic (MD) simulations which displayed a typical fold represented by other aspartic proteinase (AP) family members. Molecular docking of Pepstatin inhibitor with BuPAG 7 revealed to interact through various hydrogen bonding and hydrophobic interactions. Various amino acid substitutions were observed in peptide-binding cleft of BuPAG 7. Superimposition of BuPAG 7 with homologous structures revealed the presence of a 35-41 amino acid long insertion (alpha helix connected by two loops) near the N- terminus which seems to be a unique feature of BuPAG 7 in AP family. This is the first report on identification and sequence characterization of PAG isoforms in buffalo with unique finding that these isoforms represent many transcript variants. We also report 3D structure of the most abundant isoform BuPAG 7 for the first time.

Modeling and resistant alleles explain the selectivity of antimalarial compound 49c towards apicomplexan aspartyl proteases

Toxoplasma gondii aspartyl protease 3 (TgASP3) phylogenetically clusters with Plasmodium falciparum Plasmepsins IX and X (PfPMIX, PfPMX). These proteases are essential for parasite survival, acting as key maturases for secreted proteins implicated in invasion and egress. A potent antimalarial peptidomimetic inhibitor (49c) originally developed against Plasmepsin II selectively targets TgASP3, PfPMIX, and PfPMX To unravel the molecular basis for the selectivity of 49c, we constructed homology models of PfPMIX, PfPMX, and TgASP3 that were first validated by identifying the determinants of microneme and rhoptry substrate recognition. The flap and flap-like structures of several reported Plasmepsins are highly flexible and critically modulate the access to the binding cavity. Molecular docking of 49c to TgASP3, PfPMIX, and PfPMX models predicted that the conserved phenylalanine residues in the flap, F344, F291, and F305, respectively, account for the sensitivity toward 49c. Concordantly, phenylalanine mutations in the flap of the three proteases increase twofold to 15-fold the IC₅₀ values of 49c. Compellingly the selection of mutagenized T. gondii resistant strains to 49c reproducibly converted F344 to a cysteine residue.

Proteolytic Cleavage-Mechanisms, Function, and "Omic" Approaches for a Near-Ubiquitous Posttranslational Modification

Proteases enzymatically hydrolyze peptide bonds in substrate proteins, resulting in a widespread, irreversible posttranslational modification of the protein's structure and biological function. Often regarded as a mere degradative mechanism in destruction of proteins or turnover in maintaining physiological homeostasis, recent research in the field of degradomics has led to the recognition of two main yet unexpected concepts. First, that targeted, limited proteolytic cleavage events by a wide repertoire of proteases are pivotal regulators of most, if not all, physiological and pathological processes. Second, an unexpected in vivo abundance of stable cleaved proteins revealed pervasive, functionally relevant protein processing in normal and diseased tissue-from 40 to 70% of proteins also occur in vivo as distinct stable proteoforms with undocumented N- or C-termini, meaning these proteoforms are stable functional cleavage products, most with unknown functional implications. In this Review, we discuss the structural biology aspects and mechanisms of catalysis by different protease classes. We also provide an overview of biological pathways that utilize specific proteolytic cleavage as a precision control mechanism in protein quality control, stability, localization, and maturation, as well as proteolytic cleavage as a mediator in signaling pathways. Lastly, we provide a comprehensive overview of analytical methods and approaches to study activity and substrates of proteolytic enzymes in relevant biological models, both historical and focusing on state of the art proteomics techniques in the field of degradomics research.

Fragment-based virtual screening approach and molecular dynamics simulation studies for identification of BACE1 inhibitor leads

Traditional structure-based virtual screening method to identify drug-like small molecules for BACE1 is so far unsuccessful. Location of BACE1, poor Blood Brain Barrier permeability and P-glycoprotein (Pgp) susceptibility of the inhibitors make it even more difficult. Fragment-based drug design method is suitable for efficient optimization of initial hit molecules for target like BACE1. We have developed a fragment-based virtual screening approach to identify/optimize the fragment molecules as a starting point. This method combines the shape, electrostatic, and pharmacophoric features of known fragment molecules, bound to protein conjugate crystal structure, and aims to identify both chemically and energetically feasible small fragment ligands that bind to BACE1 active site. The two top-ranked fragment hits were subjected for a 53 ns MD simulation. Principle component analysis and free energy landscape analysis reveal that the new ligands show the characteristic features of established BACE1 inhibitors. The potent method employed in this study may serve for the development of potential lead molecules for BACE1-directed Alzheimer's disease therapeutics.

pH Dependence of Chitosan Enzymolysis

As a means of making chitosan more useful in biotechnological applications, it was hydrolyzed using pepsin, chitosanase and α-amylase. The enzymolysis behavior of these enzymes was further systematically studied for its effectiveness in the production of low-molecular-weight chitosans (LMWCs) and other derivatives. The study showed that these enzymes depend on ion hydronium (H3O+), thus on pH with a pH dependence fitting R2 value of 0.99. In y = 1.484[H^+] + 0.114, the equation of pH dependence, when [H^+] increases by one, y (k_0/k_m) increases by 1.484. From the temperature dependence study, the activation energy (Ea) and pre-exponential factor (A) were almost identical for two of the enzymes, but a considerable difference was observed in comparison with the third enzyme. Chitosanase and pepsin had nearly identical Ea, but α-amylase was significantly lower. This serves as evidence that the hydrolysis reaction of α-amylase relies on low-barrier hydrogen bonds (LBHBs), which explains its low Ea in actual conditions. The confirmation of this phenomenon was further derived from a similarly considerable difference in the order magnitudes of A between α-amylase and the other two enzymes, which was more than five. Variation of the rate constants of the enzymatic hydrolysis of chitosan with temperature follows the Arrhenius equation.

Understanding the structural basis of substrate recognition by Plasmodium falciparum plasmepsin V to aid in the design of potent inhibitors

Plasmodium falciparum plasmepsin V (PfPMV) is an essential aspartic protease required for parasite survival, thus, considered as a potential drug target. This study reports the first detailed structural analysis and molecular dynamics simulation of PfPMV as an apoenzyme and its complexes with the substrate PEXEL as well as with the inhibitor saquinavir. The presence of pro-peptide in PfPMV may not structurally hinder the formation of a functionally competent catalytic active site. The structure of PfPMV-PEXEL complex shows that the unique positions of Glu179 and Gln222 are responsible for providing the specificity of PEXEL substrate with arginine at P3 position. The structural analysis also reveals that the S4 binding pocket in PfPMV is occupied by Ile94, Ala98, Phe370 and Tyr472, and therefore, does not allow binding of pepstatin, a potent inhibitor of most pepsin-like aspartic proteases. Among the screened inhibitors, the HIV-1 protease inhibitors and KNI compounds have higher binding affinities for PfPMV with saquinavir having the highest value. The presence of a flexible group at P2 and a bulky hydrophobic group at P3 position of the inhibitor is preferred in the PfPMV substrate binding pocket. Results from the present study will aid in the design of potent inhibitors of PMV.

N-Glycosylation Improves the Pepsin Resistance of Histidine Acid Phosphatase Phytases by Enhancing Their Stability at Acidic pHs and Reducing Pepsin's Accessibility to Its Cleavage Sites

N-Glycosylation can modulate enzyme structure and function. In this study, we identified two pepsin-resistant histidine acid phosphatase (HAP) phytases from Yersinia kristensenii (YkAPPA) and Yersinia rohdei (YrAPPA), each having an N-glycosylation motif, and one pepsin-sensitive HAP phytase from Yersinia enterocolitica (YeAPPA) that lacked an N-glycosylation site. Site-directed mutagenesis was employed to construct mutants by altering the N-glycosylation status of each enzyme, and the mutant and wild-type enzymes were expressed in Pichia pastoris for biochemical characterization. Compared with those of the N-glycosylation site deletion mutants and N-deglycosylated enzymes, all N-glycosylated counterparts exhibited enhanced pepsin resistance. Introduction of the N-glycosylation site into YeAPPA as YkAPPA and YrAPPA conferred pepsin resistance, shifted the pH optimum (0.5 and 1.5 pH units downward, respectively) and improved stability at acidic pH (83.2 and 98.8% residual activities at pH 2.0 for 1 h). Replacing the pepsin cleavage sites L197 and L396 in the immediate vicinity of the N-glycosylation motifs of YkAPPA and YrAPPA with V promoted their resistance to pepsin digestion when produced in Escherichia coli but had no effect on the pepsin resistance of N-glycosylated enzymes produced in P. pastoris. Thus, N-glycosylation may improve pepsin resistance by enhancing the stability at acidic pH and reducing pepsin's accessibility to peptic cleavage sites. This study provides a strategy, namely, the manipulation of N-glycosylation, for improvement of phytase properties for use in animal feed.

Conformational Dynamics and Binding Free Energies of Inhibitors of BACE-1: From the Perspective of Protonation Equilibria

BACE-1 is the β-secretase responsible for the initial amyloidogenesis in Alzheimer’s disease, catalyzing hydrolytic cleavage of substrate in a pH-sensitive manner. The catalytic mechanism of BACE-1 requires water-mediated proton transfer from aspartyl dyad to the substrate, as well as structural flexibility in the flap region. Thus, the coupling of protonation and conformational equilibria is essential to a full in silico characterization of BACE-1. In this work, we perform constant pH replica exchange molecular dynamics simulations on both apo BACE-1 and five BACE-1-inhibitor complexes to examine the effect of pH on dynamics and inhibitor binding properties of BACE-1. In our simulations, we find that solution pH controls the conformational flexibility of apo BACE-1, whereas bound inhibitors largely limit the motions of the holo enzyme at all levels of pH. The microscopic pK_a values of titratable residues in BACE-1 including its aspartyl dyad are computed and compared between apo and inhibitor-bound states. Changes in protonation between the apo and holo forms suggest a thermodynamic linkage between binding of inhibitors and protons localized at the dyad. Utilizing our recently developed computational protocol applying the binding polynomial formalism to the constant pH molecular dynamics (CpHMD) framework, we are able to obtain the pH-dependent binding free energy profiles for various BACE-1-inhibitor complexes. Our results highlight the importance of correctly addressing the binding-induced protonation changes in protein-ligand systems where binding accompanies a net proton transfer. This work comprises the first application of our CpHMD-based free energy computational method to protein-ligand complexes and illustrates the value of CpHMD as an all-purpose tool for obtaining pH-dependent dynamics and binding free energies of biological systems.

Formation of insoluble amyloid plaques in the vascular and hippocampal areas of the brain characterizes Alzheimer’s disease, a devastating neurodegenerative disorder causing dementia. Site-specific hydrolytic catalysis of β-secretase, or BACE-1, is responsible for production of oligomerative amyloid β-peptide. As the catalytic activity of BACE-1 is pH-dependent and its structural dynamics are intrinsic to the catalysis, we examine the dependence of dynamics of BACE-1 on solution pH and its implications on the catalytic mechanism of BACE-1. Also, we highlight the importance of accurate description of protonation states of the titratable groups in computer-aided drug discovery targeting BACE-1. We hope the understanding of pH dependence of the dynamics and inhibitor binding properties of BACE-1 will aid the structure-based inhibitor design efforts against Alzheimer’s disease.

Computational Insights into Substrate and Site Specificities, Catalytic Mechanism, and Protonation States of the Catalytic Asp Dyad of β -Secretase

In this review, information regarding substrate and site specificities, catalytic mechanism, and protonation states of the catalytic Asp dyad of β-secretase (BACE1) derived from computational studies has been discussed. BACE1 catalyzes the rate-limiting step in the generation of Alzheimer amyloid beta peptide through the proteolytic cleavage of the amyloid precursor protein. Due to its biological functioning, this enzyme has been considered as one of the most important targets for finding the cure for Alzheimer's disease. Molecular dynamics (MD) simulations suggested that structural differences in the key regions (inserts A, D, and F and the 10s loop) of the enzyme are responsible for the observed difference in its activities towards the WT- and SW-substrates. The modifications in the flap, third strand, and insert F regions were found to be involved in the alteration in the site specificity of the glycosylphosphatidylinositol bound form of BACE1. Our QM and QM/MM calculations suggested that BACE1 hydrolyzed the SW-substrate more efficiently than the WT-substrate and that cleavage of the peptide bond occurred in the rate-determining step. The results from molecular docking studies showed that the information concerning a single protonation state of the Asp dyad is not enough to run an in silico screening campaign.

Exploring the binding of BACE-1 inhibitors using comparative binding energy analysis (COMBINE)

BACKGROUND

The inhibition of the activity of β-secretase (BACE-1) is a potentially important approach for the treatment of Alzheimer disease. To explore the mechanism of inhibition, we describe the use of 46 X-ray crystallographic BACE-1/inhibitor complexes to derive quantitative structure-activity relationship (QSAR) models. The inhibitors were aligned by superimposing 46 X-ray crystallographic BACE-1/inhibitor complexes, and gCOMBINE software was used to perform COMparative BINding Energy (COMBINE) analysis on these 46 minimized BACE-1/inhibitor complexes. The major advantage of the COMBINE analysis is that it can quantitatively extract key residues involved in binding the ligand and identify the nature of the interactions between the ligand and receptor.

RESULTS

By considering the contributions of the protein residues to the electrostatic and van der Waals intermolecular interaction energies, two predictive and robust COMBINE models were developed: (i) the 3-PC distance-dependent dielectric constant model (built from a single X-ray crystal structure) with a q2 value of 0.74 and an SDEC value of 0.521; and (ii) the 5-PC sigmoidal electrostatic model (built from the actual complexes present in the Brookhaven Protein Data Bank) with a q2 value of 0.79 and an SDEC value of 0.41.

CONCLUSIONS

These QSAR models and the information describing the inhibition provide useful insights into the design of novel inhibitors via the optimization of the interactions between ligands and those key residues of BACE-1.

Structure and activity of the Pseudomonas aeruginosa hotdog-fold thioesterases PA5202 and PA2801

The hotdog fold is one of the basic protein folds widely present in bacteria, archaea and eukaryotes. Many of these proteins exhibit thioesterase activity against fatty acyl-CoAs and play important roles in lipid metabolism, cellular signalling and degradation of xenobiotics. The genome of the opportunistic pathogen Pseudomonas aeruginosa contains over 20 genes encoding predicted hotdog-fold proteins, none of which have been experimentally characterized. We have found that two P. aeruginosa hotdog proteins display high thioesterase activity against 3-hydroxy-3-methylglutaryl-CoA and glutaryl-CoA (PA5202), and octanoyl-CoA (PA2801). Crystal structures of these proteins were solved (at 1.70 and 1.75 Å for PA5202 and PA2801 respectively) and revealed a hotdog fold with a potential catalytic carboxylate residue located on the long α-helix (Asp(57) in PA5202 and Glu(35) in PA2801). Alanine residue replacement mutagenesis of PA5202 identified four residues (Asn(42), Arg(43), Asp(57) and Thr(76)) that are critical for its activity and are located in the active site. A P. aeruginosa PA5202 deletion strain showed an increased secretion of the antimicrobial pigment pyocyanine and an increased expression of genes involved in pyocyanin biosynthesis, suggesting a functional link between PA5202 activity and pyocyanin production. Thus the P. aeruginosa hotdog thioesterases PA5202 and PA2801 have similar structures, but exhibit different substrate preferences and functions.

Structural studies of vacuolar plasmepsins

Plasmepsins (PMs) are pepsin-like aspartic proteases present in different species of parasite Plasmodium. Four Plasmodium spp. (P. vivax, P. ovale, P. malariae, and the most lethal P. falciparum) are mainly responsible for causing human malaria that affects millions worldwide. Due to the complexity and rate of parasite mutation coupled with regional variations, and the emergence of P. falciparum strains which are resistant to antimalarial agents such as chloroquine and sulfadoxine/pyrimethamine, there is constant pressure to find new and lasting chemotherapeutic drug therapies. Since many proteases represent therapeutic targets and PMs have been shown to play an important role in the survival of parasite, these enzymes have recently been identified as promising targets for the development of novel antimalarial drugs. The genome of P. falciparum encodes 10 PMs (PMI, PMII, PMIV-X and histo-aspartic protease (HAP)), 4 of which (PMI, PMII, PMIV and HAP) reside within the food vacuole, are directly involved in degradation of human hemoglobin, and share 50-79% amino acid sequence identity. This review focuses on structural studies of only these four enzymes, including their orthologs in other Plasmodium spp.. Almost all original crystallographic studies were performed with PMII, but more recent work on PMIV, PMI, and HAP resulted in a more complete picture of the structure-function relationship of vacuolar PMs. Many structures of inhibitor complexes of vacuolar plasmepsins, as well as their zymogens, have been reported in the last 15 years. Information gained by such studies will be helpful for the development of better inhibitors that could become a new class of potent antimalarial drugs. This article is part of a Special Issue entitled: Proteolysis 50 years after the discovery of lysosome.

Structural insights into the activation and inhibition of histo-aspartic protease from Plasmodium falciparum

Histo-aspartic protease (HAP) from Plasmodium falciparum is a promising target for the development of novel antimalarial drugs. The sequence of HAP is highly similar to those of pepsin-like aspartic proteases, but one of the two catalytic aspartates, Asp32, is replaced with histidine. Crystal structures of the truncated zymogen of HAP and of the complex of the mature enzyme with inhibitor KNI-10395 have been determined at 2.1 and 2.5 Å resolution, respectively. As in other proplasmepsins, the propeptide of the zymogen interacts with the C-terminal domain of the enzyme, forcing the N- and C-terminal domains apart, thereby separating His32 and Asp215 and preventing formation of the mature active site. In the inhibitor complex, the enzyme forms a tight domain-swapped dimer, not previously seen in any aspartic proteases. The inhibitor is found in an unprecedented conformation resembling the letter U, stabilized by two intramolecular hydrogen bonds. Surprisingly, the location and conformation of the inhibitor are similar to those of the fragment of helix 2 comprising residues 34p-38p in the prosegments of the zymogens of gastric aspartic proteases; a corresponding helix assumes a vastly different orientation in proplasmepsins. Each inhibitor molecule is in contact with two molecules of HAP, interacting with the carboxylate group of the catalytic Asp215 of one HAP protomer through a water molecule, while also making a direct hydrogen bond to Glu278A' of the other protomer. A comparison of the shifts in the positions of the catalytic residues in the inhibitor complex presented here with those published previously gives further hints regarding the enzymatic mechanism of HAP.

The role of tyrosine 71 in modulating the flap conformations of BACE1

β-Site amyloid precursor protein cleaving enzyme 1 (BACE1) is a potential target for treating Alzheimer's disease. BACE1's binding site is partially covered by a flexible loop on its N-terminal domain, known as the "flap," which has been found in several conformations in crystal structures of BACE1 and other aspartyl proteases. The side chain of the invariant residue Tyr71 on the flap adopts several rotameric orientations, leading to our hypothesis that the orientation of this residue dictates the movement and conformations available to the flap. We investigated this hypothesis by performing 220 ns of molecular dynamics simulations of bound and unbound wild-type BACE1 as well as the unbound Y71A mutant. Our findings indicate that the flap exhibits various degrees of mobility and adopts different conformations depending on the Tyr71 orientation. Surprisingly, the "self-inhibited" form is stable in our simulations, making it a reasonable target for drug design. The alanine mutant, lacking a large side chain at position 71, displays significant differences in flap dynamics from wild type, freely sampling very open and closed conformations. Our simulations show that Tyr71, in addition to its previously determined functions in catalysis and substrate binding, has the important role of modulating flap conformations in BACE1.

Crystal structures of the histo-aspartic protease (HAP) from Plasmodium falciparum

The structures of recombinant histo-aspartic protease (HAP) from malaria-causing parasite Plasmodium falciparum as apoenzyme and in complex with two inhibitors, pepstatin A and KNI-10006, were solved at 2.5-, 3.3-, and 3.05-A resolutions, respectively. In the apoenzyme crystals, HAP forms a tight dimer not seen previously in any aspartic protease. The interactions between the monomers affect the conformation of two flexible loops, the functionally important "flap" (residues 70-83) and its structural equivalent in the C-terminal domain (residues 238-245), as well as the orientation of helix 225-235. The flap is found in an open conformation in the apoenzyme. Unexpectedly, the active site of the apoenzyme contains a zinc ion tightly bound to His32 and Asp215 from one monomer and to Glu278A from the other monomer, with the coordination of Zn resembling that seen in metalloproteases. The flap is closed in the structure of the pepstatin A complex, whereas it is open in the complex with KNI-10006. Although the binding mode of pepstatin A is significantly different from that in other pepsin-like aspartic proteases, its location in the active site makes unlikely the previously proposed hypothesis that HAP is a serine protease. The binding mode of KNI-10006 is unusual compared with the binding of other inhibitors from the KNI series to aspartic proteases. The novel features of the HAP active site could facilitate design of specific inhibitors used in the development of antimalarial drugs.

Statistical coupling analysis of aspartic proteinases based on crystal structures of the Trichoderma reesei enzyme and its complex with pepstatin A

The crystal structures of an aspartic proteinase from Trichoderma reesei (TrAsP) and of its complex with a competitive inhibitor, pepstatin A, were solved and refined to crystallographic R-factors of 17.9% (R(free)=21.2%) at 1.70 A resolution and 15.8% (R(free)=19.2%) at 1.85 A resolution, respectively. The three-dimensional structure of TrAsP is similar to structures of other members of the pepsin-like family of aspartic proteinases. Each molecule is folded in a predominantly beta-sheet bilobal structure with the N-terminal and C-terminal domains of about the same size. Structural comparison of the native structure and the TrAsP-pepstatin complex reveals that the enzyme undergoes an induced-fit, rigid-body movement upon inhibitor binding, with the N-terminal and C-terminal lobes tightly enclosing the inhibitor. Upon recognition and binding of pepstatin A, amino acid residues of the enzyme active site form a number of short hydrogen bonds to the inhibitor that may play an important role in the mechanism of catalysis and inhibition. The structures of TrAsP were used as a template for performing statistical coupling analysis of the aspartic protease family. This approach permitted, for the first time, the identification of a network of structurally linked residues putatively mediating conformational changes relevant to the function of this family of enzymes. Statistical coupling analysis reveals coevolved continuous clusters of amino acid residues that extend from the active site into the hydrophobic cores of each of the two domains and include amino acid residues from the flap regions, highlighting the importance of these parts of the protein for its enzymatic activity.

Crystal structure of an active form of BACE1, an enzyme responsible for amyloid beta protein production

BACE1 (beta-secretase) is a transmembrane aspartic protease that cleaves the beta-amyloid precursor protein and generates the amyloid beta peptide (Abeta). BACE1 cycles between the cell surface and the endosomal system many times and becomes activated interconvertibly during its cellular trafficking, leading to the production of Abeta. Here we report the crystal structure of the catalytically active form of BACE1. The active form has novel structural features involving the conformation of the flap and subsites that promote substrate binding. The functionally essential residues and water molecules are well defined and play a key role in the iterative activation of BACE1. We further describe the crystal structure of the dehydrated form of BACE1, showing that BACE1 activity is dependent on the dynamics of a catalytically required Asp-bound water molecule, which directly affects its catalytic properties. These findings provide insight into a novel regulation of BACE1 activity and elucidate how BACE1 modulates its activity during cellular trafficking.

Novel method for probing the specificity binding profile of ligands: applications to HIV protease

A detailed understanding of factors influencing the binding specificity of a ligand to a set of desirable targets and undesirable decoys is a key step in the design of potent and selective therapeutics. We have developed a general method for optimizing binding specificity in ligand-receptor complexes based on the theory of electrostatic charge optimization. This methodology can be used to tune the binding of a ligand to a panel of potential targets and decoys, along the continuum from narrow binding to only one partner to broad binding to the entire panel. Using HIV-1 protease as a model system, we probe specificity in three distinct ways. First, we probe interactions that could make the promiscuous protease inhibitor pepstatin more selective toward HIV-1 protease. Next, we study clinically approved HIV-1 protease inhibitors and probe ways to broaden the binding profiles toward both wild-type HIV-1 protease and drug-resistant mutants. Finally, we study a conformational ensemble of wild-type HIV-1 protease to 'design in' broad specificity to known drugs before resistance mutations arise. The results from this conformational ensemble were similar to those from the drug-resistant ensemble, suggesting the use of a conformational wild-type ensemble as a tool to develop escape-mutant-resistant inhibitors.

Assigning the protonation states of the key aspartates in β-Secretase using QM/MM X-ray structure refinement

β-Secretase, a.k.a. β-APP cleaving enzyme (BACE), is an aspartyl protease that has been implicated as a key target in the pathogenesis of Alzheimer's disease (AD). The identification of the protonation states of the key aspartates in β-secretase is of great interest both in understanding the reaction mechanism and in guiding the design of drugs against AD. However, the resolutions of currently available crystal structures for BACE are not sufficient to determine the hydrogen atom locations. We have assigned the protonation states of the key aspartates using a novel method, QM/MM X-ray refinement. In our approach, an energy function is introduced to the refinement where the atoms in the active site are modeled by quantum mechanics (QM) and the other atoms are represented by molecular mechanics (MM). The gradients derived from the QM/MM energy function are combined with those from the X-ray target to refine the crystal structure of a complex containing BACE and an inhibitor. A total number of 8 protonation configurations of the aspartyl dyad were considered and QM/MM X-ray refinements were performed for all of them. The relative stability of the refined structures was scored by constructing the thermodynamic cycle using the energetics calculated by fully quantum mechanical self-consistent reaction field (QM/SCRF) calculations. While all 8 refined structures fit the observed electron density about equally well, we find the mono-protonated configurations to be strongly favored energetically, especially the configuration with the inner oxygen of Asp32 protonated and the hydroxyl of the inhibitor pointing towards Asp228. It was also found that these results depend on the constraints imposed by the X-ray data. We suggest that one of the strengths of this approach is that the resulting structures are a consensus of theoretical and experimental data and remark on the significance of our results in structure based drug design and mechanistic studies.

Enzymic and structural characterization of nepenthesin, a unique member of a novel subfamily of aspartic proteinases

Carnivorous plants are known to secrete acid proteinases to digest prey, mainly insects, for nitrogen uptake. In the present study, we have purified, for the first time, to homogeneity two acid proteinases (nepenthesins I and II) from the pitcher fluid of Nepenthes distillatoria (a pitcher-plant known locally as badura) and investigated their enzymic and structural characteristics. Both enzymes were optimally active at pH approx. 2.6 towards acid-denatured haemoglobin; the specificity of nepenthesin I towards oxidized insulin B chain appears to be similar, but slightly wider than those of other APs (aspartic proteinases). Among the enzymic properties, however, the most notable is their unusual stability: both enzymes were remarkably stable at or below 50 degrees C, especially nepenthesin I was extremely stable over a wide range of pH from 3 to 10 for over 30 days. This suggests an evolutionary adaptation of the enzymes to their specific habitat. We have also cloned the cDNAs and deduced the complete amino acid sequences of the precursors of nepenthesins I and II (437 and 438 residues respectively) from the pitcher tissue of N. gracilis. Although the corresponding mature enzymes (each 359 residues) are homologous with ordinary pepsin-type APs, both enzymes had a high content of cysteine residues (12 residues/molecule), which are assumed to form six unique disulphide bonds as suggested by computer modelling and are supposed to contribute towards the remarkable stability of nepenthesins. Moreover, the amino acid sequence identity of nepenthesins with ordinary APs, including plant vacuolar APs, is remarkably low (approx. 20%), and phylogenetic comparison shows that nepenthesins are distantly related to them to form a novel subfamily of APs with a high content of cysteine residues and a characteristic insertion, named 'the nepenthesin-type AP-specific insertion', that includes a large number of novel, orthologous plant APs emerging in the gene/protein databases.