Crystal structure of papain-E64-c complex. Binding diversity of E64-c to papain S2 and S3 subsites

Structural Characterization and Physiological Role of Bombyx mori Fibroinase in the Silk Gland Development

Silkworm is a highly valuable insect that produces silk through secretion by a silk gland. Within this gland, a type of cathepsin L protease called Fibroinase was identified as an enzyme for hydrolyzing the primary components of silk, including fibroin and sericin. Here, we determined the crystal structure of Fibroinase fromBombyx mori at a resolution of 1.56 Å. Comparative structural analysis revealed that Fibroinase adopted a similar structural pattern with papain-type cathepsin, consisting of an N-terminal domain and a C-terminal domain. The interface between the domains forms a substrate-binding cleft, where the E64 inhibitor noncovalently binds in a novel manner. Additionally, computational simulations combined with biochemical analysis allowed us to define the binding mode and inhibition mechanism of physiological inhibitor Bombyx cysteine protease inhibitor (BCPI) with Fibroinase. Moreover, the expression profiles and RNA interference of Fibroinase indicated its critical role in removing silk proteins in the silk gland lumen and the destruction of silk gland tissue during the larval-pupal metamorphosis. These findings enhance our understanding of the structural and biochemical features of Fibroinase and its inhibitors, while also providing evidence for the physiological role of Fibroinase in silk gland development.

Unveiling the Mechanistic Singularities of Caspases: A Computational Analysis of the Reaction Mechanism in Human Caspase-1

Caspases are cysteine proteases in charge of breaking a peptide bond next to an aspartate residue. Caspases constitute an important family of enzymes involved in cell death and inflammatory processes. A plethora of diseases, including neurological and metabolic diseases and cancer, are associated with the poor regulation of caspase-mediated cell death and inflammation. Human caspase-1 in particular carries out the transformation of the pro-inflammatory cytokine pro-interleukin-1β into its active form, a key process in the inflammatory response and then in many diseases, such as Alzheimer's disease. Despite its importance, the reaction mechanism of caspases has remained elusive. The standard mechanistic proposal valid for other cysteine proteases and that involves the formation of an ion pair in the catalytic dyad is not supported by experimental evidence. Using a combination of classical and hybrid DFT/MM simulations, we propose a reaction mechanism for the human caspase-1 that explains experimental observations, including mutagenesis, kinetic, and structural data. In our mechanistic proposal, the catalytic cysteine, Cys285, is activated after a proton transfer to the amide group of the scissile peptide bond, a process facilitated by hydrogen-bond interactions with Ser339 and His237. The catalytic histidine does not directly participate in any proton transfer during the reaction. After formation of the acylenzyme intermediate, the deacylation step takes place through the activation of a water molecule by the terminal amino group of the peptide fragment formed during the acylation step. The overall activation free energy obtained from our DFT/MM simulations is in excellent agreement with the value derived from the experimental rate constant, 18.7 vs 17.9 kcal·mol^-1, respectively. Simulations of the H237A mutant support our conclusions and agree with the reported reduced activity observed for this caspase-1 variant. We propose that this mechanism can explain the reactivity of all cysteine proteases belonging to the CD clan and that differences with respect to other clans could be related to the larger preference showed by enzymes of the CD clan for charged residues at position P1. This mechanism would avoid the free energy penalty associated with the formation of an ion pair. Finally, our structural description of the reaction process can be useful to assist in the design of inhibitors of caspase-1, a target in the treatment of several human diseases.

The Structural Basis of African Swine Fever Virus pS273R Protease Binding to E64 through Molecular Dynamics Simulations

Identification of novel drugs for anti-African swine fever (ASF) applications is of utmost urgency, as it negatively affects pig farming and no effective vaccine or treatment is currently available. African swine fever virus (ASFV) encoded pS273R is a cysteine protease that plays an important role in virus replication. E64, acting as an inhibitor of cysteine protease, has been established as exerting an inhibitory effect on pS273R. In order to obtain a better understanding of the interaction between E64 and pS273R, common docking, restriction docking, and covalent docking were employed to analyze the optimal bonding position between pS273R-E64 and its bonding strength. Additionally, three sets of 100 ns molecular dynamics simulations were conducted to examine the conformational dynamics of pS273R and the dynamic interaction of pS273R-E64, based on a variety of analytical methods including root mean square deviation (RMSD), root mean square fluctuation (RMSF), free energy of ligand (FEL), principal component analysis (PCA), and molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) analysis. The results show that E64 and pS273R exhibited close binding degrees at the activity center of ASFV pS273R protease. The data of these simulations indicate that binding of E64 to pS273R results in a reduction in flexibility, particularly in the ARM region, and a change in the conformational space of pS273R. Additionally, the ability of E64 to interact with polar amino acids such as ASN158, SER192, and GLN229, as well as charged amino acids such as LYS167 and HIS168, seems to be an important factor in its inhibitory effect. Finally, Octet biostratigraphy confirmed the binding of E64 and pS273R with a KD value of 903 uM. Overall, these findings could potentially be utilized in the development of novel inhibitors of pS273R to address the challenges posed by ASFV.

Giant magnetoresistive biosensors for real-time quantitative detection of protease activity

Proteases are enzymes that cleave proteins and are crucial to physiological processes such as digestion, blood clotting, and wound healing. Unregulated protease activity is a biomarker of several human diseases. Synthetic peptides that are selectively hydrolyzed by a protease of interest can be used as reporter substrates of unregulated protease activity. We developed an activity-based protease sensor by immobilizing magnetic nanoparticles (MNPs) to the surface of a giant magnetoresistive spin-valve (GMR SV) sensor using peptides. Cleavage of these peptides by a protease releases the magnetic nanoparticles resulting in a time-dependent change in the local magnetic field. Using this approach, we detected a significant release of MNPs after 3.5 minutes incubation using just 4 nM of the cysteine protease, papain. In addition, we show that proteases in healthy human urine do not release the MNPs, however addition of 20 nM of papain to the urine samples resulted in a time-dependent change in magnetoresistance. This study lays the foundation for using GMR SV sensors as a platform for real-time, quantitative detection of protease activity in biological fluids.

Fluorescent Probes for Studying Thioamide Positional Effects on Proteolysis Reveal Insight into Resistance to Cysteine Proteases

Thioamide substitutions of the peptide backbone have been shown to reduce proteolytic degradation, and this property can be used to generate competitive protease inhibitors and to stabilize peptides toward degradation in vivo. Here, we present a straightforward sensor design that allows a systematic study of the positional effects of thioamide substitution by using real-time fluorescence. Thioamide scanning in peptide substrates of five papain family cysteine proteases demonstrates that a thioamide at or near the scissile bond can slow proteolysis in all cases, but that the magnitude of the effects varies with position and protease in spite of high sequence homology. Mechanistic investigation of papain proteolysis reveals that the thioamide effects derive from reductions in both affinity (KM ) and turnover number (kcat ). Computational modeling allows these effects to be understood based on disruption of key enzyme-substrate hydrogen bonds, providing a model for future rational use of thioamides to confer cysteine protease resistance.

Dual binding mode of "bitter sugars" to their human bitter taste receptor target

The 25 human bitter taste receptors (hTAS2Rs) are responsible for detecting bitter molecules present in food, and they also play several physiological and pathological roles in extraoral compartments. Therefore, understanding their ligand specificity is important both for food research and for pharmacological applications. Here we provide a molecular insight into the exquisite molecular recognition of bitter β-glycopyranosides by one of the members of this receptor subclass, hTAS2R16. Most of its agonists have in common the presence of a β-glycopyranose unit along with an extremely structurally diverse aglycon moiety. This poses the question of how hTAS2R16 can recognize such a large number of "bitter sugars". By means of hybrid molecular mechanics/coarse grained molecular dynamics simulations, here we show that the three hTAS2R16 agonists salicin, arbutin and phenyl-β-D-glucopyranoside interact with the receptor through a previously unrecognized dual binding mode. Such mechanism may offer a seamless way to fit different aglycons inside the binding cavity, while maintaining the sugar bound, similar to the strategy used by several carbohydrate-binding lectins. Our prediction is validated a posteriori by comparison with mutagenesis data and also rationalizes a wealth of structure-activity relationship data. Therefore, our findings not only provide a deeper molecular characterization of the binding determinants for the three ligands studied here, but also give insights applicable to other hTAS2R16 agonists. Together with our results for other hTAS2Rs, this study paves the way to improve our overall understanding of the structural determinants of ligand specificity in bitter taste receptors.

Plant Proteases Involved in Regulated Cell Death

Each plant genome encodes hundreds of proteolytic enzymes. These enzymes can be divided into five distinct classes: cysteine-, serine-, aspartic-, threonine-, and metalloproteinases. Despite the differences in their structural properties and activities, members of all of these classes in plants are involved in the processes of regulated cell death - a basic feature of eukaryotic organisms. Regulated cell death in plants is an indispensable mechanism supporting plant development, survival, stress responses, and defense against pathogens. This review summarizes recent advances in studies of plant proteolytic enzymes functioning in the initiation and execution of distinct types of regulated cell death.

Electrostatic complementarity in pseudoreceptor modeling based on drug molecule crystal structures: the case of loxistatin acid (E64c)

A combination of pseudoreceptor modeling and electrostatic complementarity maps properties of a native pocket for an enzyme ligand.

A Mitofusin-2-dependent inactivating cleavage of Opa1 links changes in mitochondria cristae and ER contacts in the postprandial liver

Hepatic metabolism requires mitochondria to adapt their bioenergetic and biosynthetic output to accompany the ever-changing anabolic/catabolic state of the liver cell, but the wiring of this process is still largely unknown. Using a postprandial mouse liver model and quantitative cryo-EM analysis, we show that when the hepatic mammalian target of rapamycin complex 1 (mTORC1) signaling pathway disengages, the mitochondria network fragments, cristae density drops by 30%, and mitochondrial respiratory capacity decreases by 20%. Instead, mitochondria-ER contacts (MERCs), which mediate calcium and phospholipid fluxes between these organelles, double in length. These events are associated with the transient expression of two previously unidentified C-terminal fragments (CTFs) of Optic atrophy 1 (Opa1), a mitochondrial GTPase that regulates cristae biogenesis and mitochondria dynamics. Expression of Opa1 CTFs in the intermembrane space has no effect on mitochondria morphology, supporting a model in which they are intermediates of an Opa1 degradation program. Using an in vitro assay, we show that these CTFs indeed originate from the cleavage of Opa1 at two evolutionarily conserved consensus sites that map within critical folds of the GTPase. This processing of Opa1, termed C-cleavage, is mediated by the activity of a cysteine protease whose activity is independent from that of Oma1 and presenilin-associated rhomboid-like (PARL), two known Opa1 regulators. However, C-cleavage requires Mitofusin-2 (Mfn2), a key factor in mitochondria-ER tethering, thereby linking cristae remodeling to MERC assembly. Thus, in vivo, mitochondria adapt to metabolic shifts through the parallel remodeling of the cristae and of the MERCs via a mechanism that degrades Opa1 in an Mfn2-dependent pathway.

Influence of azide incorporation on binding affinity by small papain inhibitors

In order to develop affinity-based biosensor platforms, appropriate ligands with a functional handle for immobilization onto a biosensor surface are required. To this end, a library of papain inhibitors was designed and synthesized, containing different azide linkers for subsequent immobilization by 'click' chemistry, in this particular case by copper-free, strain-promoted azide-alkyne cycloaddition (SPAAC). Furthermore, a molecular docking study was performed to obtain a better insight as to at which position such azide handles could be tolerated without affecting binding affinity. Although the azide moiety is small, in some cases its introduction strongly influenced the binding affinity. For one class of inhibitors a swapped binding mode was proposed to explain the results. In addition, a specific site for linker introduction was identified, which did not significantly affect the binding affinity.

Quantum mechanics/molecular mechanics studies of the mechanism of falcipain-2 inhibition by the epoxysuccinate E64

Because of the increasing resistance of malaria parasites to antimalarial drugs, the lack of highly effective vaccines, and an inadequate control of mosquito vectors, the problem is growing, especially in the developing world. New approaches to drug development are consequently required. One of the proteases involved in the degradation of human hemoglobin is named falcipain-2 (FP2), which has emerged as a promising target for the development of novel antimalarial drugs. However, very little is known about the inhibition of FP2. In this paper, the inhibition of FP2 by the epoxysuccinate E64 has been studied by molecular dynamics (MD) simulations using hybrid AM1d/MM and M06-2X/MM potentials to obtain a complete picture of the possible free energy reaction paths. A thorough analysis of the reaction mechanism has been conducted to understand the inhibition of FP2 by E64. According to our results, the irreversible attack of Cys42 on E64 can take place on both carbon atoms of the epoxy ring because both processes present similar barriers. While the attack on the C2 atom presents a slightly smaller barrier (12.3 vs 13.6 kcal mol(-1)), the inhibitor-protein complex derived from the attack on C3 appears to be much more stabilized. In contrast to previous hypotheses, our results suggest that residues such as Gln171, Asp170, Gln36, Trp43, Asn81, and even His174 would be anchoring the inhibitor in a proper orientation for the reaction to take place. These results may be useful for the rational design of new compounds with higher inhibitory activity.

Chemical biology of antigen presentation by MHC molecules

MHC class I and MHC class II molecules present peptides to the immune system to drive proper T cell responses. Pharmacological modulation of T-cell responses can offer treatment options for a range of immune-related diseases. Pharmacological downregulation of MHC molecules may find application in treatment of auto-immunity and transplantation rejection while pharmacological activation of antigen presentation would support immune responses to infection and cancer. Since the cell biology of MHC class I and MHC class II antigen presentation is understood in great detail, many potential targets for manipulation have been defined over the years. Here, we discuss how antigen presentation by MHC molecules can be modulated by pharmacological agents and how chemistry can further support the study of antigen presentation in general. The chemical biology of antigen presentation by MHC molecules shows surprising options for immune modulation and the development of future therapies.

Dissecting the active site of the collagenolytic cathepsin L3 protease of the invasive stage of Fasciola hepatica

Background

A family of secreted cathepsin L proteases with differential activities is essential for host colonization and survival in the parasitic flatworm Fasciola hepatica. While the blood feeding adult secretes predominantly FheCL1, an enzyme with a strong preference for Leu at the S₂ pocket of the active site, the infective stage produces FheCL3, a unique enzyme with collagenolytic activity that favours Pro at P₂.

Methodology/Principal Findings

Using a novel unbiased multiplex substrate profiling and mass spectrometry methodology (MSP-MS), we compared the preferences of FheCL1 and FheCL3 along the complete active site cleft and confirm that while the S₂ imposes the greatest influence on substrate selectivity, preferences can be indicated on other active site subsites. Notably, we discovered that the activity of FheCL1 and FheCL3 enzymes is very different, sharing only 50% of the cleavage sites, supporting the idea of functional specialization. We generated variants of FheCL1 and FheCL3 with S₂ and S₃ residues by mutagenesis and evaluated their substrate specificity using positional scanning synthetic combinatorial libraries (PS-SCL). Besides the rare P₂ Pro preference, FheCL3 showed a distinctive specificity at the S₃ pocket, accommodating preferentially the small Gly residue. Both P₂ Pro and P₃ Gly preferences were strongly reduced when Trp67 of FheCL3 was replaced by Leu, rendering the enzyme incapable of digesting collagen. In contrast, the inverse Leu67Trp substitution in FheCL1 only slightly reduced its Leu preference and improved Pro acceptance in P₂, but greatly increased accommodation of Gly at S₃.

Conclusions/Significance

These data reveal the significance of S₂ and S₃ interactions in substrate binding emphasizing the role for residue 67 in modulating both sites, providing a plausible explanation for the FheCL3 collagenolytic activity essential to host invasion. The unique specificity of FheCL3 could be exploited in the design of specific inhibitors selectively directed to specific infective stage parasite proteinases.

The flatworm Fasciola hepatica is responsible for fasciolosis, one of the most common parasitic diseases of livestock worldwide, with increased incidence of human cases. When contaminated plants are ingested, infective larvae are released and transverse the gut wall before migrating to the bile ducts within the liver. Migrating liver flukes erode host tissue while adults feed on blood and they mature and release thousands of eggs. Several developmentally-regulated cathepsin L like proteolytic enzymes (FheCLs) are essential to the migrating and feeding processes. Despite being similar in structure and sequence these enzymes show specialization attacking preferentially different substrates and taking part in the diverse process of invasion, immune evasion and feeding. Our analyses reveal unique differences in activity between the major infective juvenile (FheCL3) and adult (FheCL1) enzymes, and demonstrate that the juvenile enzyme has a particular active site that allows it to degrade collagen, the main component of connective tissues. We demonstrate that a single position on the active site, residue 67, is essential to this collagenolytic activity critical for parasite invasion.

[Overview of structural study on conformations and intermolecular interactions of biomolecules]

Information on the conformational feature and specific intermolecular interaction of biomolecules is important to understand the biological function and to develop device for treating disorder caused by the abnormal function. Thus the 3D structures of the biologically active molecules and the specific interactions with their target molecules at the atomic level have been investigated by various physicochemical approaches. Herein, the following five subjects are reviewed: (1) function-linked conformations of biomolecules including natural annular products, opioid peptides and neuropeptides; (2) π-π stacking interactions of tryptophan derivatives with coenzymes and nucleic acid bases; (3) mRNA cap recognition of eukaryotic initiation factor 4E and its regulation by 4E-binding protein; (4) conformational feature of histamine H2 receptor antagonists and design of cathepsin B inhibitors; (5) self-aggregation mechanism of tau protein and its inhibition.

Subclassification and biochemical analysis of plant papain-like cysteine proteases displays subfamily-specific characteristics

Papain-like cysteine proteases (PLCPs) are a large class of proteolytic enzymes associated with development, immunity, and senescence. Although many properties have been described for individual proteases, the distribution of these characteristics has not been studied collectively. Here, we analyzed 723 plant PLCPs and classify them into nine subfamilies that are present throughout the plant kingdom. Analysis of these subfamilies revealed previously unreported distinct subfamily-specific functional and structural characteristics. For example, the NPIR and KDEL localization signals are distinctive for subfamilies, and the carboxyl-terminal granulin domain occurs in two PLCP subfamilies, in which some individual members probably evolved by deletion of the granulin domains. We also discovered a conserved double cysteine in the catalytic site of SAG12-like proteases and two subfamily-specific disulfides in RD19A-like proteases. Protease activity profiling of representatives of the PLCP subfamilies using novel fluorescent probes revealed striking polymorphic labeling profiles and remarkably distinct pH dependency. Competition assays with peptide-epoxide scanning libraries revealed common and unique inhibitory fingerprints. Finally, we expand the detection of PLCPs by identifying common and organ-specific protease activities and identify previously undetected proteases upon labeling with cell-penetrating probes in vivo. This study provides the plant protease research community with tools for further functional annotation of plant PLCPs.

Foot and mouth disease leader protease (Lbpro): Investigation of prime side specificity allows the synthesis of a potent inhibitor

Foot and mouth disease virus expresses its genetic information as a single polyprotein that is translated from the single-stranded RNA genome. Proteinases contained within the polyprotein then generate the mature viral proteins. The leader protease (Lb(pro)) performs the initial cleavage by freeing itself from the growing polypeptide chain; subsequently, Lb(pro) cleaves the two homologues of the host cell protein eukaryotic initiation factor 4G (eIF4G). We showed that Lb(pro) possesses specific binding sites at the non prime side from S(1) down to S(7) [Santos et al. (2009) Biochemistry, 48, 7948-7958]. Here, we demonstrate that Lb(pro) has high prime side specificity at least down to the S'(5) site. Lb(pro) is thus not only one of the smallest papain-like cysteine peptidases but also one of the most specific. It can still however cleave between both K↓G and G↓R pairs. We further determined the two-step irreversible inhibition (E + I ↔ EI→ E - I) kinetic parameters of two known irreversible epoxide-based inhibitors of cysteine proteinases, E64 and CA074 on Lb(pro) that show for the reversible step (E + I ↔ EI) K(i) = 3.4 μM and 11.6 μM, and for the irreversible step (EI→E-I) k(4) = 0.16 and 0.06 min(-1), respectively. Knowledge of the Lb(pro) specificity led us to extend E64 by addition of the dipeptide R-P. This compound, termed E64-R-P-NH(2), irreversibly inhibited Lb(pro) with a K(i) = 30 nM and k(4) = 0.01 min(-1) and can serve as the basis for design of specific inhibitors of FMDV replication.

Cathepsin cleavage potentiates the Ebola virus glycoprotein to undergo a subsequent fusion-relevant conformational change

Cellular entry of Ebola virus (EBOV), a deadly hemorrhagic fever virus, is mediated by the viral glycoprotein (GP). The receptor-binding subunit of GP must be cleaved (by endosomal cathepsins) in order for entry and infection to proceed. Cleavage appears to proceed through 50-kDa and 20-kDa intermediates, ultimately generating a key 19-kDa core. How 19-kDa GP is subsequently triggered to bind membranes and induce fusion remains a mystery. Here we show that 50-kDa GP cannot be triggered to bind to liposomes in response to elevated temperature but that 20-kDa and 19-kDa GP can. Importantly, 19-kDa GP can be triggered at temperatures ∼10°C lower than 20-kDa GP, suggesting that it is the most fusion ready form. Triggering by heat (or urea) occurs only at pH 5, not pH 7.5, and involves the fusion loop, as a fusion loop mutant is defective in liposome binding. We further show that mild reduction (preferentially at low pH) triggers 19-kDa GP to bind to liposomes, with the wild-type protein being triggered to a greater extent than the fusion loop mutant. Moreover, mild reduction inactivates pseudovirion infection, suggesting that reduction can also trigger 19-kDa GP on virus particles. Our results support the hypothesis that priming of EBOV GP, specifically to the 19-kDa core, potentiates GP to undergo subsequent fusion-relevant conformational changes. Our findings also indicate that low pH and an additional endosomal factor (possibly reduction or possibly a process mimicked by reduction) act as fusion triggers.

Selective and reversible thiol-pegylation, an effective approach for purification and characterization of five fully active ficin (iso)forms from Ficus carica latex

The latex of Ficus carica constitutes an important source of many proteolytic components known under the general term of ficin (EC 3.4.22.3) which belongs to the cysteine proteases of the papain family. So far, no data on the purification and characterization of individual forms of these proteases are available. An effective strategy was used to fractionate and purify to homogeneity five ficin forms, designated A, B, C, D1 and D2 according to their sequence of elution from a cation-exchange chromatographic support. Following rapid fractionation on a SP-Sepharose Fast Flow column, the different ficin forms were chemically modified by a specific and reversible monomethoxypolyethylene glycol (mPEG) reagent. In comparison with their un-derivatized counterparts, the mPEG-protein derivatives behaved differently on the ion-exchanger, allowing us for the first time to obtain five highly purified ficin molecular species titrating 1mol of thiol group per mole of enzyme. The purified ficins were characterized by de novo peptide sequencing and peptide mass fingerprinting analyzes, using mass spectrometry. Circular dichroism measurements indicated that all five ficins were highly structured, both in term of secondary and tertiary structure. Furthermore, analysis of far-UV CD spectra allowed calculation of their secondary structural content. Both these data and the molecular masses determined by MS reinforce the view that the enzymes belong to the family of papain-like proteases. The five ficin forms also displayed different specific amidase activities against small synthetic substrates like dl-BAPNA and Boc-Ala-Ala-Gly-pNA, suggesting some differences in their active site organization. Enzymatic activity of the five ficin forms was completely inhibited by specific cysteine and cysteine/serine proteases inhibitors but was unaffected by specific serine, aspartic and metallo proteases inhibitors.

Influenza virus hemagglutinin spike neck architectures and interaction with model enzymes evaluated by MALDI-TOF mass spectrometry and bioinformatics tools

Interactions between model enzymes and the influenza virus hemagglutinin (HA) homotrimeric spike were addressed. We digested influenza virions (naturally occurring strains and laboratory reassortants) with bromelain or subtilisin Carlsberg and analyzed by MALDI-TOF mass spectrometry the resulting HA2 C-terminal segments. All cleavage sites, together with (minor) sites detected in undigested HAs, were situated in the linker region that connects the transmembrane domain to the ectodomain. In addition to cleavage at highly favorable amino acids, various alternative enzyme preferences were found that strongly depended on the HA subtype/type. We also evaluated the surface electrostatic potentials, binding cleft topographies and spatial dimensions of stem bromelain (homologically modeled) and subtilisin Carlsberg (X-ray resolved). The results show that the enzymes (∼45Å(3)) would hardly fit into the small (∼18-20Å) linker region of the HA-spike. However, the HA membrane proximal ectodomain region was predicted to be intrinsically disordered. We propose that its motions allow steric adjustment of the enzymes' active sites to the neck of the HA spike. The subtype/type-specific architectures in this region also influenced significantly the cleavage preferences of the enzymes.

Peptidase inhibitors in the MEROPS database

The MEROPS website (http://merops.sanger.ac.uk) includes information on peptidase inhibitors as well as on peptidases and their substrates. Displays have been put in place to link peptidases and inhibitors together. The classification of protein peptidase inhibitors is continually being revised, and currently inhibitors are grouped into 67 families based on comparisons of protein sequences. These families can be further grouped into 38 clans based on comparisons of tertiary structure. Small molecule inhibitors are important reagents for peptidase characterization and, with the increasing importance of peptidases as drug targets, they are also important to the pharmaceutical industry. Small molecule inhibitors are now included in MEROPS and over 160 summaries have been written.

Structural studies of specific intermolecular interactions and self-aggregation of biomolecules and their application to drug design

Information on the structural basis of intermolecular recognition or self-aggregation of biomolecules at the atomic level is important to understand biological functions and to develop devices for treating disorders caused by abnormal functions. Thus structural analysis of specific intermolecular or intramolecular interactions of biomolecules has been performed using various physicochemical approaches. Herein, the following three subjects are reviewed: (1) structural analyses of mRNA cap structure recognition by eukaryotic initiation factor 4E and its functional regulation by endogenous 4E-binding protein; (2) structural studies of self-aggregation mechanism of microtubule-binding domain in tau protein and aggregation inhibitor; and (3) molecular design of cathepsin B-specific inhibitor.

Atomistic insights into the inhibition of cysteine proteases: first QM/MM calculations clarifying the stereoselectivity of epoxide-based inhibitors

Due to their important role in many diseases, cysteine proteases represent new promising drug targets. An important class of cysteine-protease inhibitors is derived from the naturally occurring compound E64, possessing an epoxysuccinyl moiety as warhead. Experimental studies show stereoselectivity concerning the inhibition potency, e.g., a trans-configured epoxide ring is essential for inhibition, and furthermore, in most cases, the ( S, S)-configured inhibitors have a higher inhibition potency than their ( R, R)-counterparts. However, the underlying effects are not fully understood. In this work, such effects are investigated by classical molecular dynamics simulations and combined quantum mechanics/molecular modeling (QM/MM) calculations for the E64c-cathepsin B complex. Our computations reveal that the hydrogen bonding network between the enzyme and the E64c (or its derivatives) determines the stereoselectivity of the subsequent ring opening reaction by governing the distance between the attacking thiolate and the attacked C2 atom of the epoxide ring. For the ( S, S)-configuration, a strong network can be realized which enables a close contact between the reacting centers, so that the irreversible step becomes very efficient. The ( R, S)-configuration ( cis-configuration) can only form networks in which the two reacting centers are so far away from each other that the irreversible step can hardly happen. The ( R, R)-configuration is in between, less optimal than the ( S, S)-configuration but much better than the ( R, S)-configuration. Exceptions where the ( R, R)-configurations shows higher potency than the ( S, S) ones are also explained.

Molecular docking of cathepsin L inhibitors in the binding site of papain

The papain/CLIK-148 coordinate system was employed as a model to study the interactions of a nonpeptide thiocarbazate inhibitor of cathepsin L ( 1). This small molecule inhibitor, a thiol ester containing a diacyl hydrazine functionality and one stereogenic center, was most active as the S-enantiomer, with an IC 50 of 56 nM; the R-enantiomer ( 2) displayed only weak activity (33 microM). Correspondingly, molecular docking studies with Extra Precision Glide revealed a correlation between score and biological activity for the two thiocarbazate enantiomers when a structural water was preserved. The molecular interactions between 1 and papain were very similar to the interactions observed for CLIK-148 ( 3a and 3b) with papain, especially with regard to the hydrogen-bonding and lipophilic interactions of the ligands with conserved residues in the catalytic binding site. Subsequent docking of virtual compounds in the binding site led to the identification of a more potent inhibitor ( 5), with an IC 50 of 7.0 nM. These docking studies revealed that favorable energy scores and correspondingly favorable biological activities could be realized when the virtual compound design included occupation of the S2, S3, and S1' subsites by hydrophobic and aromatic functionalities of the ligand, and at least three hydrogen bonding contacts between the ligand and the conserved binding site residues of the protein.

Atomistic insights into the inhibition of cysteine proteases: first QM/MM calculations clarifying the regiospecificity and the inhibition potency of epoxide- and aziridine-based inhibitors

Epoxides and aziridines are important building blocks for inhibitors of cysteine proteases which are promising drug targets for many diseases. In spite of the large amount of experimental data concerning inhibition potency, structure-activity relationships, and structural arrangements of enzyme-inhibitor complexes, little is known about the basic principles which connect the substitution pattern with the resulting activities. To shed some light on this issue which is essential for the rational design of improved compounds, we have studied the inhibition processes theoretically for various inhibitors using quantum mechanical/molecular mechanical hybrid approaches and classical molecular dynamics simulations. The careful analysis of the computational results allows insight into the interactions which govern the regio- and stereospecificity of the interactions. Known structure-activity relationships are rationalized in terms of the same interactions that determine the measured pH dependencies. Inconsistencies in existing X-ray structures are resolved through comparison with the computed structures, which leads to a reassessment of the factors that control the inhibition potency. Similarities and differences in the mode of action of epoxide- and aziridine-based inhibitors are elucidated. Finally the small reaction barriers computed for the irreversible step in E64 analogues call into question the commonly accepted two-step model of inhibition since the second, irreversible step is predicted to be so fast that suitably oriented enzyme-inhibitor complexes will react rather than dissociate and equilibrate.

Fungal effector protein AVR2 targets diversifying defense-related cys proteases of tomato

The interaction between the fungal pathogen Cladosporium fulvum and its host tomato (Solanum lycopersicum) is an ideal model to study suppression of extracellular host defenses by pathogens. Secretion of protease inhibitor AVR2 by C. fulvum during infection suggests that tomato papain-like cysteine proteases (PLCPs) are part of the tomato defense response. We show that the tomato apoplast contains a remarkable diversity of PLCP activities with seven PLCPs that fall into four different subfamilies. Of these PLCPs, transcription of only PIP1 and RCR3 is induced by treatment with benzothiadiazole, which triggers the salicylic acid-regulated defense pathway. Sequencing of PLCP alleles of tomato relatives revealed that only PIP1 and RCR3 are under strong diversifying selection, resulting in variant residues around the substrate binding groove. The doubled number of variant residues in RCR3 suggests that RCR3 is under additional adaptive selection, probably to prevent autoimmune responses. AVR2 selectively inhibits only PIP1 and RCR3, and one of the naturally occurring variant residues in RCR3 affects AVR2 inhibition. The higher accumulation of PIP1 protein levels compared with RCR3 indicates that PIP1 might be the real virulence target of AVR2 and that RCR3 acts as a decoy for AVR2 perception in plants carrying the Cf-2 resistance gene.

Identification and characterization of 3-substituted pyrazolyl esters as alternate substrates for cathepsin B: the confounding effects of DTT and cysteine in biological assays

Substituted pyrazole esters were identified as hits in a high throughput screen (HTS) of the NIH Molecular Libraries Small Molecule Repository (MLSMR) to identify inhibitors of the enzyme cathepsin B. Members of this class, along with functional group analogs, were synthesized in an effort to define the structural requirements for activity. Analog characterization was hampered by the need to include a reducing agent such as dithiothreitol (DTT) or cysteine in the assay, highlighting the caution required in interpreting biological data gathered in the presence of such nucleophiles. Despite the confounding effects of DTT and cysteine, our studies demonstrate that the pyrazole 1 acts as alternate substrate for cathepsin B, rather than as an inhibitor.

Crystal structure of NS-134 in complex with bovine cathepsin B: a two-headed epoxysuccinyl inhibitor extends along the entire active-site cleft

The crystal structure of the inhibitor NS-134 in complex with bovine cathepsin B reveals that functional groups attached to both sides of the epoxysuccinyl reactive group bind to the part of active-site cleft as predicted. The -Leu-Pro-OH side binds to the primed binding sites interacting with the His110 and His111 residues with its C-terminal carboxy group, whereas the -Leu-Gly-Meu (-Leu-Gly-Gly-OMe) part (Meu, methoxycarbonylmethyl) binds along the non-primed binding sites. Comparison with the propeptide structures of cathepsins revealed that the binding of the latter part is least similar to the procathepsin B structure; this result, together with the two-residue shift in positioning of the Leu-Gly-Gly part, suggests that the propeptide structures of the cognate enzymes may not be the best starting point for the design of reverse binding inhibitors.

Structural basis for development of cathepsin B-specific noncovalent-type inhibitor: crystal structure of cathepsin B-E64c complex

In order to elucidate the substrate specificity of the Sn subsites (n=1-3) of cathepsin B, its crystal structure inhibited by E64c [(+)-(2S,3S)-3-(1-[N-(3-methylbutyl)amino]-leucylcarbonyl)oxirane-2-carboxylic acid] was analyzed by the X-ray diffraction method. Iterative manual rebuilding and convenient conjugate refinement of structure decreased R- and free R-factors to 19.7% and to 23.9%, respectively, where 130 water molecules were included for the refinement using 14,759 independent reflections from 10 to 2.3 A resolution. The epoxy carbonyl carbon of E64c was covalently bonded to the Cys(29) S(gamma) atom and the remaining parts were located at Sn subsites (n=1-3). The substrate specificity of these subsites was characterized based on their interactions with the inhibitor. Base on these structural data, we developed a novel cathepsin B-specific noncovalent-type inhibitor, which may bind to S2'-S3. The molecular design of possessing structural elements of both CA074 and E64c, assisted by energy minimization and molecular dynamics (MD) simulation, may lead to a new lead noncovalent-type inhibitor.

Insight to structural subsite recognition in plant thiol protease-inhibitor complexes : understanding the basis of differential inhibition and the role of water

BACKGROUND

This work represents an extensive MD simulation / water-dynamics studies on a series of complexes of inhibitors (leupeptin, E-64, E-64-C, ZPACK) and plant cysteine proteases (actinidin, caricain, chymopapain, calotropin DI) of papain family to understand the various interactions, water binding mode, factors influencing it and the structural basis of differential inhibition.

RESULTS

The tertiary structure of the enzyme-inhibitor complexes were built by visual interactive modeling and energy minimization followed by dynamic simulation of 120 ps in water environment. DASA study with and without the inhibitor revealed the potential subsite residues involved in inhibition. Though the interaction involving main chain atoms are similar, critical inspection of the complexes reveal significant differences in the side chain interactions in S2-P2 and S3-P3 pairs due to sequence differences in the equivalent positions of respective subsites leading to differential inhibition.

CONCLUSION

The key finding of the study is a conserved site of a water molecule near oxyanion hole of the enzyme active site, which is found in all the modeled complexes and in most crystal structures of papain family either native or complexed. Conserved water molecules at the ligand binding sites of these homologous proteins suggest the structural importance of the water, which changes the conventional definition of chemical geometry of inhibitor binding domain, its shape and complimentarity. The water mediated recognition of inhibitor to enzyme subsites (Pn.H2O.Sn) of leupeptin acetyl oxygen to caricain, chymopapain and calotropinDI is an additional information and offer valuable insight to potent inhibitor design.

Structural basis of inhibition of cysteine proteases by E-64 and its derivatives

This paper focuses on the inhibitory mechanism of E-64 and its derivatives (epoxysuccinyl-based inhibitors) with some cysteine proteases, based on the binding modes observed in the x-ray crystal structures of their enzyme-inhibitor complexes. E-64 is a potent irreversible inhibitor against general cysteine proteases, and its binding modes with papain, actinidin, cathepsin L, and cathepsin K have been reviewed at the atomic level. E-64 interacts with the Sn subsites of cysteine proteases. Although the Sn-Pn (n = 1-3) interactions of the inhibitor with the main chains of the active site residues are similar in respective complexes, the significant difference is observed in the side-chain interactions of S2-P2 and S3-P3 pairs because of different residues constituting the respective subsites. E-64-c and CA074 are representative derivatives developed from E-64 as a clinical usable and a cathepsin B-specific inhibitors, respectively. In contrast with similar binding/inhibitory modes of E-64-c and E-64 for cysteine proteases, the inhibitory mechanism of cathepsin B-specific CA074 results from the binding to the Sn' subsite.

Probing the specificity of cysteine proteinases at subsites remote from the active site: analysis of P4, P3, P2' and P3' variations in extended substrates

We have determined the kinetic parameters for the hydrolysis by papain, cathepsin B and cathepsin L of internally quenched fluorescent peptides derived from the lead peptides Abz-AAFRSAQ-EDDnp [in which Abz and EDDnp stand for o-aminobenzoic acid and N-(2,4-dinitrophenyl)ethylenediamine respectively], to map the specificity of S(4) and S(3) subsites, and Abz-AFRSAAQ-EDDnp, to identify the specificity of S(2)' and S(3)'. Abz and EDDnp were the fluorescent quencher pair. These two series of peptides were cleaved at the Arg-Ser bond and systematic modifications at P(4), P(3), P(2)' and P(3)' were made. The S(4) to S(2)' subsites had a significant influence on the hydrolytic efficiencies of the three enzymes. Only papain activity was observed to be dependent on S(3)', indicating that its binding site is larger than those of cathepsins B and L. Hydrophobic amino acids were accepted at S(4), S(3), S(2)' and S(3)' of the three enzymes. The best substrates for cathepsins L and B had Trp and Asn at P(2)' respectively; variations at this position were less accepted by these enzymes. The best substrates for papain were peptides containing Trp, Tyr or Asn at P(3)'. Basic residues at P(3) and P(4) were well accepted by cathepsin L and papain. We also explored the susceptibility of substrates Abz-AFRSXAQ-EDDnp, modified at P(2)' (X), to human cathepsin B mutants from which one or two occluding loop contacts had been removed. The modifications at His(111) (H111A) and His(110) (H110A) of cathepsin B led to an increase in k(cat) values of one or two orders of magnitude. The hydrolytic efficiencies of these cathepsin B mutants became closer to those of papain or cathepsin L.

Inhibition mechanism of cathepsin L-specific inhibitors based on the crystal structure of papain-CLIK148 complex

Papain was used as an experimental model structure to understand the inhibition mechanism of newly developed specific inhibitors of cathepsin L, the papain superfamily. Recently, we developed a series of cathepsin L-specific inhibitors which are called the CLIK series [(1999) FEBS Lett. 458, 6-10]. Here, we report the complex structure of papain with CLIK148, which is a representative inhibitor from the CLIK series. The inhibitor complex structure was solved at 1.7 A resolution with conventional R 0.177. Unlike other epoxisuccinate inhibitors (E64, CA030, and CA074), CLIK148 uses both prime and nonprime sites, which are important for the specific inhibitory effect on cathepsin L. Also, the specificity for cathepsin L could be explained by the existence of Phe in the P2 site and hydrophobic interaction of N-terminal pyridine ring.

Binding modes of a new epoxysuccinyl-peptide inhibitor of cysteine proteases. Where and how do cysteine proteases express their selectivity?

Papain from Carica papaya, an easily available cysteine protease, is the best-studied representative of this family of enzymes. The three dimensional structure of papain is very similar to that of other cysteine proteases of either plant (actinidin, caricain, papaya protease IV) or animal (cathepsins B, K, L, H) origin. As abnormalities in the activities of mammalian cysteine proteases accompany a variety of diseases, there has been a long-lasting interest in the development of potent and selective inhibitors for these enzymes. A covalent inhibitor of cysteine proteases, designed as a combination of epoxysuccinyl and peptide moieties, has been modeled in the catalytic pocket of papain. A number of its configurations have been generated and relaxed by constrained simulated annealing-molecular dynamics in water. A clear conformational variability of this inhibitor is discussed in the context of a conspicuous conformational diversity observed earlier in several solid-state structures of other complexes between cysteine proteases and covalent inhibitors. The catalytic pockets S2 and even more so S3, as defined by the pioneering studies on the papain-ZPACK, papain-E64c and papain-leupeptin complexes, appear elusive in view of the evident flexibility of the present inhibitor and in confrontation with the obvious conformational scatter seen in other examples. This predicts limited chances for the development of selective structure-based inhibitors of thiol proteases, designed to exploit the minute differences in the catalytic pockets of various members of this family. A simultaneous comparison of the three published proenzyme structures suggests the enzyme's prosegment binding loop-prosegment interface as a new potential target for selective inhibitors of papain-related thiol proteases.

Design and synthesis of dipeptidyl alpha',beta'-epoxy ketones, potent irreversible inhibitors of the cysteine protease cruzain

The dipeptidyl epoxy ketone 4 is a potent, irreversible inhibitor of cruzain, a cysteine protease isolated from Trypanosoma cruzi, with an apparent second order inhibition rate constant of 330,000 sec-1M-1.

X-ray crystal structure of papain complexed with cathepsin B-specific covalent-type inhibitor: substrate specificity and inhibitory activity

The Ile-Pro sequence of CA074, potent covalent-type inhibitor, is necessary to exhibit the specificity for cathepsin B, but not for papain. In order to elucidate how its sequence binds to papain and why such binding does not exhibit the specificity for papain at the atomic level, two CA074-related compounds, 1 (N-(L-3-carboxyloxirane-2-carbonyl)-L-isoleucyl-L-proline) and 2 (N-(L-3-carboxyloxirane-2-carbonyl)-L-isoleucyl-diethylamide), were designed and their structure--inhibitory activity relationship was investigated by the X-ray crystal analyses of the complexes with papain. The Ile-Pro moiety of 1 was located at the S2 and S3 subsites consisting of Val-133, Val-157, and Asp-158 and of Tyr-61, Gly-66, and Tyr-67 residues of papain, respectively, which is in contrast with the binding of CA074 to S'n (n = 1 approximately 2) subsites in the complex with cathepsin B. Although 2 in the complex with papain showed the similar binding pattern to 1, its inhibitory activity was about two-fold higher than of 1, suggesting the importance of tight S3-P3 hydrophobic interaction for the activity. The difference of the substrate specificity between papain and cathepsin B has also been discussed based on the X-ray results of the present and cathepsin B-inhibitor complexes.

Revised definition of substrate binding sites of papain-like cysteine proteases

A review of kinetic and structural data has enabled us to reconsider the definition of substrate binding sites in papain-like cysteine proteases. Only three substrate binding sites, S2, S1 and S1', involve main as well as side chain contacts between substrate and enzyme residues. Interactions between the enzymes and the substrate P3 and P2' residues are based on side chains (an exception is cathepsin B which is a carboxydipeptidase), so their interaction surface spreads over a relatively wide area. The location and definition of substrate binding sites beyond S3 and S2' is even more questionable.

Substrate/propeptide-derived endo-epoxysuccinyl peptides as highly potent and selective cathepsin B inhibitors

Based on recent information about the anti-substrate binding mode of the propeptide portion of procathepsin B and the well established substrate-like binding of epoxysuccinyl-dipeptide carboxylates to the S' subsites of cathepsin B a new endo-trans-epoxysuccinyl peptide was synthesized that contains the dipeptide moiety Leu-Pro-OH for the P1'-P2' substrate positions and the tripeptide moiety Leu-Gly-Gly-OMe (sequence portion 46-48 of the propeptide) for the P2-P4 positions in anti-substrate orientation. With an unequivocal (2S,3S) configuration this new trans-epoxysuccinyl peptide derivative was found to inhibit cathepsin B with an apparent second-order rate constant of 1,520,000 M(-1) s(-1) which represents so far the most potent inhibitor among E-64-derived compounds. Conversely, the (2R,3R) diastereomer exhibited a significantly lower inhibition potency. This observation fully agrees with our previous findings that inhibitor/enzyme interactions at the S subsites are favored by the (2S,3S) and reverse interactions at the S' subsites by the (2R,3R) configuration of the trans-epoxysuccinyl moiety.

E-64 analogues as inhibitors of cathepsin B. On the role of the absolute configuration of the epoxysuccinyl group

A series of trans-epoxysuccinyl-peptide derivatives based on the natural inhibitor E-64 were synthesized in the (2R,3R) and (2S,3S) configuration in order to analyze the role of the stereochemistry of this residue in dictating inhibitory potency and selectivity for cysteine proteases. We confirmed that binding of E-64 like trans-epoxysuccinyl compounds is remarkably favored by the (2S,3S) configuration, but we also found that CA030-type compounds are stronger inhibitors in the (2R,3R) configuration than the related diastereomers. Consequently, the structural requirements for exploiting both the S and S' subsites are not additive and a structure-based design of bis-peptidyl derivatives of trans-epoxysuccinic acid to increase selective inhibition becomes even more difficult. Additional contrasting effects were observed for the pH optima required in the electrostatic interactions at the S and S' subsites.

The crystal structure of human cathepsin L complexed with E-64

We have determined the three dimensional structure of the complex of human cathepsin L and E-64, an irreversible inhibitor of cysteine proteases, at 2.5 A resolution. The overall structure was similar to that of other known cysteine proteases and apparently identical to the mature region of procathepsin L. The electron density for E-64 is clearly visible except for the guanidinobutane moiety. From comparison of the active sites of cathepsin L and B, we found the following: (1) The S' subsites of cathepsin L and B are totally different because of the 'occluding loop' lying on the end of the S' subsites of cathepsin B. (2) The S2 pocket of cathepsin L is shallow and narrow compared to that of cathepsin B. (3) The S3 subsites of the two enzymes are more similar than the other subsites, but cathepsin L may accommodate a more bulky group at this site. Knowledge of the active site structure of cathepsin L should be helpful for the structure-based design of potent and specific inhibitors which are of therapeutic importance.

A structure of the complex between concanavalin A and methyl-3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside reveals two binding modes

The structure of concanavalin A in complex with the trimannoside methyl-3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside has been determined in a novel space group. In three of the four subunits of the concanavalin A tetramer, the interactions between the protein and the bound saccharide are essentially identical to those reported previously by other authors (Naismith, J. H., and Field, R. A. (1996) J. Biol. Chem. 271, 972-976). In the fourth subunit, however, the alpha1-->3 linkage has a different conformation, resulting in a different part of the alpha1-->3-linked mannose interacting with essentially the same surface of the protein. Furthermore, significant differences are observed in the quaternary associations of the subunits compared with the saccharide-free structures and other carbohydrate complexes, suggesting that the concanavalin A tetramer is a rather flexible entity.

Inhibition of cathepsin B by its propeptide: use of overlapping peptides to identify a critical segment

Ten overlapping 15-mer peptides (peptidyl amides) spanning the proregion of rat cathepsin B (residues 1p-60p) were constructed to identify minimal segments having inhibitory activity towards the mature enzyme, that could be used to develop a new generation of peptide-derived inhibitors specifically targeting the active site of the corresponding proteinase. Three synthetic peptides, containing the pentapeptide Leu-Cys-Gly-Thr-Val (residues 41p-45p) in their sequence, inhibited cathepsin B with Ki values in the micromolar range. Alkylation of the thiol group of Cys-42p of peptide PB8 (36p-50p) resulted in its rapid proteolytic degradation, suggesting that this residue is essential for inhibition. The inhibition constant was slightly improved (Ki = 2 microM) using a longer peptide (26p-50p) which was completely resistant to cleavage even after a prolonged incubation. Alkylation of its cysteinyl residue also resulted in rapid cleavage of the peptide chain. Peptides derived from the rat cathepsin B prosequence also inhibited human cathepsin B with similar Ki values. Unlike rat cathepsin B, which cleaves peptide PB8 at the G47p-G48p bond after prolonged incubation, the human enzyme cleaved both PB8 and PB11 at the Lys-40p-Leu-41p bond, in agreement with the different kinetic properties of these two proteinases. New probes with improved specificity for cysteine proteinases may therefore be designed based on the sequences of their propeptides.