Crystal structure of human osteoclast cathepsin K complex with E-64

Integrated machine learning-based virtual screening and biological evaluation for identification of potential inhibitors against cathepsin K

Cathepsin K is a type of cysteine proteinase that is primarily expressed in osteoclasts and has a key role in the breakdown of bone matrix protein during bone resorption. Many studies suggest that the deficiency of cathepsin K is concomitant with a suppression of osteoclast functioning, therefore rendering the resorptive properties of cathepsin K the most prominent target for osteoporosis. This innovative work has identified a novel anti-osteoporotic agent against Cathepsin K by using a comparison of machine learning and deep learning-based virtual screening followed by their biological evaluation. Out of ten shortlisted compounds, five of the compounds (JFD02945, JFD02944, RJC01981, KM08968 and SB01934) exhibit more than 50% inhibition of the Cathepsin K activity at 0.1 μM concentration and are considered to have a promising inhibitory effect against Cathepsin K. The comprehensive docking, MD simulation, and MM/PBSA investigations affirm the stable and effective interaction of these compounds with Cathepsin K to inhibit its function. Furthermore, the compounds RJC01981, KM08968 and SB01934 are represented to have promising anti-osteoporotic properties for the management of osteoporosis owing to their significantly well predicted ADMET properties.

Cathepsin H: molecular characteristics and clues to function and mechanism

Cathepsin H (CatH) is a lysosomal cysteine protease with a unique aminopeptidase activity that is extensively expressed in the lung, pancreas, thymus, kidney, liver, skin, and brain. Owing to its specific enzymatic activity, CatH has critical effects on the regulation of biological behaviours of cancer cells and pathological processes in brain diseases. Moreover, a neutral pH level is optimal for CatH activity, so it is expected to be active in the extra-lysosomal and extracellular space. In the present review, we describe the expression, maturation, and enzymatic properties of CatH, and summarize the available experimental evidence that mechanistically links CatH to various physiological and pathological processes. Finally, we discuss the challenges and potentials of CatH inhibitors in CatH-induced disease therapy.

Structural characterization of the hypothetical protein Lpg2622, a new member of the C1 family peptidases from Legionella pneumophila

The Legionella pneumophila type II secretion system can promote bacterial growth under a wide variety of conditions and mediates the secretion of more than 25 proteins, including the uncharacterized effector Lpg2622. Here, we determined the crystal structures of apo-Lpg2622 and Lpg2622 in complex with the cysteine protease inhibitor E64. Structural analysis suggests that Lpg2622 belongs to the C1 family peptidases. Interestingly, unlike the other structurally resolved papain-like cysteine proteases, the propeptide of Lpg2622 forms a novel super-secondary structural fold (hairpin-turn-helix) and can be categorized into a new group. In addition, the N-terminal β-sheet of the Lpg2622 propeptide plays a regulatory role on enzymatic activity. This study enhances our understanding of the classification and regulatory mechanisms of the C1 family peptidases.

Evaluation of variability in high-resolution protein structures by global distance scoring

Systematic analysis of the statistical and dynamical properties of proteins is critical to understanding cellular events. Extraction of biologically relevant information from a set of high-resolution structures is important because it can provide mechanistic details behind the functional properties of protein families, enabling rational comparison between families. Most of the current structural comparisons are pairwise-based, which hampers the global analysis of increasing contents in the Protein Data Bank. Additionally, pairing of protein structures introduces uncertainty with respect to reproducibility because it frequently accompanies other settings for superimposition. This study introduces intramolecular distance scoring for the global analysis of proteins, for each of which at least several high-resolution structures are available. As a pilot study, we have tested 300 human proteins and showed that the method is comprehensively used to overview advances in each protein and protein family at the atomic level. This method, together with the interpretation of the model calculations, provide new criteria for understanding specific structural variation in a protein, enabling global comparison of the variability in proteins from different species.

Analysis of non-peptidic compounds as potential malarial inhibitors against Plasmodial cysteine proteases via integrated virtual screening workflow

Falcipain-2 (FP-2) and falcipain-3 (FP-3), haemoglobin-degrading enzymes in Plasmodium falciparum, are validated drug targets for the development of effective inhibitors against malaria. However, no commercial drug-targeting falcipains has been developed despite their central role in the life cycle of the parasites. In this work, in silico approaches are used to identify key structural elements that control the binding and selectivity of a diverse set of non-peptidic compounds onto FP-2, FP-3 and homologues from other Plasmodium species as well as human cathepsins. Hotspot residues and the underlying non-covalent interactions, important for the binding of ligands, are identified by interaction fingerprint analysis between the proteases and 2-cyanopyridine derivatives (best hits). It is observed that the size and chemical type of substituent groups within 2-cyanopyridine derivatives determine the strength of protein-ligand interactions. This research presents novel results that can further be exploited in the structure-based molecular-guided design of more potent antimalarial drugs.

Affinity Crystallography: A New Approach to Extracting High-Affinity Enzyme Inhibitors from Natural Extracts

Natural products are an important source of novel drug scaffolds. The highly variable and unpredictable timelines associated with isolating novel compounds and elucidating their structures have led to the demise of exploring natural product extract libraries in drug discovery programs. Here we introduce affinity crystallography as a new methodology that significantly shortens the time of the hit to active structure cycle in bioactive natural product discovery research. This affinity crystallography approach is illustrated by using semipure fractions of an actinomycetes culture extract to isolate and identify a cathepsin K inhibitor and to compare the outcome with the traditional assay-guided purification/structural analysis approach. The traditional approach resulted in the identification of the known inhibitor antipain (1) and its new but lower potency dehydration product 2, while the affinity crystallography approach led to the identification of a new high-affinity inhibitor named lichostatinal (3). The structure and potency of lichostatinal (3) was verified by total synthesis and kinetic characterization. To the best of our knowledge, this is the first example of isolating and characterizing a potent enzyme inhibitor from a partially purified crude natural product extract using a protein crystallographic approach.

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.

Cathepsin K: a unique collagenolytic cysteine peptidase

Cathepsin K has emerged as a promising target for the treatment of osteoporosis in recent years. Initially identified as a papain-like cysteine peptidase expressed in high levels in osteoclasts, the important role of this enzyme in bone metabolism was highlighted by the finding that mutations in the CTSK gene cause the rare recessive disorder pycnodysostosis, which is characterized by severe bone anomalies. At the molecular level, the physiological role of cathepsin K is reflected by its unique cleavage pattern of type I collagen molecules, which is fundamentally different from that of other endogenous collagenases. Several cathepsin K inhibitors have been developed to reduce the excessive bone matrix degradation associated with osteoporosis, with the frontrunner odanacatib about to successfully conclude Phase 3 clinical trials. Apart from osteoclasts, cathepsin K is expressed in different cell types throughout the body and is involved in processes of adipogenesis, thyroxine liberation and peptide hormone regulation. Elevated activity of cathepsin K has been associated with arthritis, atherosclerosis, obesity, schizophrenia, and tumor metastasis. Accordingly, its activity is tightly regulated via multiple mechanisms, including competitive inhibition by endogenous macromolecular inhibitors and allosteric regulation by glycosaminoglycans. This review provides a state-of-the-art description of the activity of cathepsin K at the molecular level, its biological functions and the mechanisms involved in its regulation.

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.

Docking-based virtual screening of covalently binding ligands: an orthogonal lead discovery approach

In pharmaceutical industry, lead discovery strategies and screening collections have been predominantly tailored to discover compounds that modulate target proteins through noncovalent interactions. Conversely, covalent linkage formation is an important mechanism for a quantity of successful drugs in the market, which are discovered in most cases by hindsight instead of systematical design. In this article, the implementation of a docking-based virtual screening workflow for the retrieval of covalent binders is presented considering human cathepsin K as a test case. By use of the docking conditions that led to the best enrichment of known actives, 44 candidate compounds with unknown activity on cathepsin K were finally selected for experimental evaluation. The most potent inhibitor, 4-(N-phenylanilino)-6-pyrrolidin-1-yl-1,3,5-triazine-2-carbonitrile (CP243522), showed a K(i) of 21 nM and was confirmed to have a covalent reversible mechanism of inhibition. The presented approach will have great potential in cases where covalent inhibition is the desired drug discovery strategy.

Cysteine cathepsins: from structure, function and regulation to new frontiers

It is more than 50 years since the lysosome was discovered. Since then its hydrolytic machinery, including proteases and other hydrolases, has been fairly well identified and characterized. Among these are the cysteine cathepsins, members of the family of papain-like cysteine proteases. They have unique reactive-site properties and an uneven tissue-specific expression pattern. In living organisms their activity is a delicate balance of expression, targeting, zymogen activation, inhibition by protein inhibitors and degradation. The specificity of their substrate binding sites, small-molecule inhibitor repertoire and crystal structures are providing new tools for research and development. Their unique reactive-site properties have made it possible to confine the targets simply by the use of appropriate reactive groups. The epoxysuccinyls still dominate the field, but now nitriles seem to be the most appropriate “warhead”. The view of cysteine cathepsins as lysosomal proteases is changing as there is now clear evidence of their localization in other cellular compartments. Besides being involved in protein turnover, they build an important part of the endosomal antigen presentation. Together with the growing number of non-endosomal roles of cysteine cathepsins is growing also the knowledge of their involvement in diseases such as cancer and rheumatoid arthritis, among others. Finally, cysteine cathepsins are important regulators and signaling molecules of an unimaginable number of biological processes. The current challenge is to identify their endogenous substrates, in order to gain an insight into the mechanisms of substrate degradation and processing. In this review, some of the remarkable advances that have taken place in the past decade are presented. This article is part of a Special Issue entitled: Proteolysis 50 years after the discovery of lysosome.

Trypanosoma cruzi: insights into naphthoquinone effects on growth and proteinase activity

In this study we compared the effects of naphthoquinones (α-lapachone, β-lapachone, nor-β-lapachone and Epoxy-α-lap) on growth of Trypanosoma cruzi epimastigotes forms, and on viability of VERO cells. In addition we also experimentally analyzed the most active compounds inhibitory profile against T. cruzi serine- and cysteine-proteinases activity and theoretically evaluated them against cruzain, the major T. cruzi cysteine proteinase by using a molecular docking approach. Our results confirmed β-lapachone and Epoxy-α-lap with a high trypanocidal activity in contrast to α-lapachone and nor-β-lapachone whereas Epoxy-α-lap presented the safest toxicity profile against VERO cells. Interestingly the evaluation of the active compounds effects against T. cruzi cysteine- and serine-proteinases activities revealed different targets for these molecules. β-Lapachone is able to inhibit the cysteine-proteinase activity of T. cruzi proteic whole extract and of cruzain, similar to E-64, a classical cysteine-proteinase inhibitor. Differently, Epoxy-α-lap inhibited the T. cruzi serine-proteinase activity, similar to PMSF, a classical serine-proteinase inhibitor. In agreement to these biological profiles in the enzymatic assays, our theoretical analysis showed that E-64 and β-lapachone interact with the cruzain specific S2 pocket and active site whereas Epoxy-α-lap showed no important interactions. Overall, our results infer that β-lapachone and Epoxy-α-lap compounds may inhibit T. cruzi epimastigotes growth by affecting T. cruzi different proteinases. Thus the present data shows the potential of these compounds as prototype of protease inhibitors on drug design studies for developing new antichagasic compounds.

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.

Species differences between human and rat in the substrate specificity of cathepsin K

Cathepsin K is known to play an important role in bone resorption, and it has the P2 specificity for proline. Rat cathepsin K has 88% identity with the human enzyme. However, it has been reported that its enzymatic activity for a Cbz-Leu-Arg-MCA substrate is lower than that of human cathepsin K, and that the rat enzyme is not well inhibited by human cathepsin K inhibitors. For this study, we prepared recombinant enzyme to investigate the substrate specificity of rat cathepsin K. Cleavage experiments using the fragment of type I collagen and peptidic libraries demonstrated that rat cathepsin K preferentially hydrolyses the substrates at the P2 Hyp position. Comparison of the S2 site between rat and human cathepsin K sequences indicated that two S2 residues at Ser134 and Val160 in rat are varied to Ala and Leu, respectively, in the human enzyme. Cleavage experiments using two single mutants, S134A and V160L, and one double mutant, S134A/V160L, of rat cathepsin K showed that all the rat mutants lost the P2 Hyp specificity. The information obtained from our comparative studies on rat and human cathepsin K should make a significant impact on developing specific inhibitors of human cathepsin K since rat is usually used as test species.

The crystal and molecular structures of a cathepsin K:chondroitin sulfate complex

Cathepsin K is the major collagenolytic enzyme produced by bone-resorbing osteoclasts. We showed earlier that the unique triple-helical collagen-degrading activity of cathepsin K depends on the formation of complexes with bone-or cartilage-resident glycosaminoglycans, such as chondroitin 4-sulfate (C4-S). Here, we describe the crystal structure of a 1:n complex of cathepsin K:C4-S inhibited by E64 at a resolution of 1.8 A. The overall structure reveals an unusual "beads-on-a-string"-like organization. Multiple cathepsin K molecules bind specifically to a single cosine curve-shaped strand of C4-S with each cathepsin K molecule interacting with three disaccharide residues of C4-S. One of the more important sets of interactions comes from a single turn of helix close to the N terminus of the proteinase containing a basic amino acid triplet (Arg8-Lys9-Lys10) that forms multiple hydrogen bonds either to the caboxylate or to the 4-sulfate groups of C4-S. Altogether, the binding sites with C4-S are located in the R-domain of cathepsin K and are distant from its active site. This explains why the general proteolytic activity of cathepsin K is not affected by the binding of chondroitin sulfate. Biochemical analyses of cathepsin K and C4-S mixtures support the presence of a 1:n complex in solution; a dissociation constant, K(d), of about 10 nM was determined for the interaction between cathepsin K and C4-S.

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.

Molecular dynamics simulations of the catalytic pathway of a cysteine protease: a combined QM/MM study of human cathepsin K

Molecular dynamics simulations using a combined QM/MM potential have been performed to study the catalytic mechanism of human cathepsin K, a member of the papain family of cysteine proteases. We have determined the two-dimensional free energy surfaces of both acylation and deacylation steps to characterize the reaction mechanism. These free energy profiles show that the acylation step is rate limiting with a barrier height of 19.8 kcal/mol in human cathepsin K and of 29.3 kcal/mol in aqueous solution. The free energy of activation for the deacylation step is 16.7 kcal/mol in cathepsin K and 17.8 kcal/mol in aqueous solution. The reduction of free energy barrier is achieved by stabilization of the oxyanion in the transition state. Interestingly, although the "oxyanion hole" has been formed in the Michaelis complex, the amide units do not donate hydrogen bonds directly to the carbonyl oxygen of the substrate, but they stabilize the thiolate anion nucleophile. Hydrogen-bonding interactions are induced as the substrate amide group approaches the nucleophile, moving more than 2 A and placing the oxyanion in contact with Gln19 and the backbone amide of Cys25. The hydrolysis of peptide substrate shares a common mechanism both for the catalyzed reaction in human cathepsin K and for the uncatalyzed reaction in water. Overall, the nucleophilic attack by Cys25 thiolate and the proton-transfer reaction from His162 to the amide nitrogen are highly coupled, whereas a tetrahedral intermediate is formed along the nucleophilic reaction pathway.

Biochemical properties and regulation of cathepsin K activity

Cysteine cathepsins (11 in humans) are mostly located in the acidic compartments of cells. They have been known for decades to be involved in intracellular protein degradation as housekeeping proteases. However, the discovery of new cathepsins, including cathepsins K, V and F, has provided strong evidence that they also participate in specific biological events. This review focuses on the current knowledge of cathepsin K, the major bone cysteine protease, which is a drug target of clinical interest. Nevertheless, we will not discuss recent developments in cathepsin K inhibitor design since they have been extensively detailed elsewhere. We will cover features of cathepsin K structure, cellular and tissue distribution, substrate specificity, and regulation (pH, propeptide, glycosaminoglycans, oxidants), and its putative roles in physiological or pathophysiological processes. Finally, we will review the kinetic data of its inhibition by natural endogenous inhibitors (stefin B, cystatin C, H- and L-kininogens).

(2S,3S)-Oxirane-2,3-dicarboxylic acid: A privileged platform for probing human cysteine cathepsins

The notion that human cysteine cathepsins contribute only to general protein turnover within the lysosomes has changed in the last decade in a substantial manner. A continuously growing number of data accumulated in different fields of life sciences revealed that these enzymes are involved in a variety of pivotal physiological processes. To investigate these particular fraction of proteolytical activity of the human degradome even in a complex cellular environment, chemical probes that covalently label the corresponding proteases proved to be versatile tools. (2S,3S)-Oxirane-2,3-dicarboxylic acid provides an ideal platform for the design of such probing systems. Depending on the complexity of the attached recognition elements, either the activity of the entire group of human cysteine cathepsins or individual members can be detected.

Protein-reactive natural products

Researchers in the post-genome era are confronted with the daunting task of assigning structure and function to tens of thousands of encoded proteins. To realize this goal, new technologies are emerging for the analysis of protein function on a global scale, such as activity-based protein profiling (ABPP), which aims to develop active site-directed chemical probes for enzyme analysis in whole proteomes. For the pursuit of such chemical proteomic technologies, it is helpful to derive inspiration from protein-reactive natural products. Natural products use a remarkably diverse set of mechanisms to covalently modify enzymes from distinct mechanistic classes, thus providing a wellspring of chemical concepts that can be exploited for the design of active-site-directed proteomic probes. Herein, we highlight several examples of protein-reactive natural products and illustrate how their mechanisms of action have influenced and continue to shape the progression of chemical proteomic technologies like ABPP.

Crucial Commitment of Proteolytic Activity of a Purified Recombinant Major House Dust Mite Allergen Der p1 to Sensitization toward IgE and IgG Responses

The major proteolytic allergen derived from the house dust mite Dermatophagoides pteronyssinus, Der p1, is one of the most clinically relevant allergens worldwide. In the present study, we evaluate the contribution of the proteolytic activity and structure of a highly purified rDer p 1 to immune responses. Mice were i.p. immunized with three forms of rDer p 1 adsorbed to Alum: one enzymatically active, one treated with an irreversible cysteine protease-specific inhibitor, E-64, and one heat denatured. Immunization with E-64-treated or heat-denatured rDer p 1 elicited much less production of serum total IgE and not only rDer p 1-specific IgE but also IgGs compared with immunization with active rDer p 1. Assays for Ab-binding and its inhibition and structural analyses indicated that E-64-treated rDer p 1 retained its global structure and conformational B cell epitopes. A proliferative response and production of IL-5 by spleen cells restimulated with rDer p 1 were observed on immunization with the active rDer p 1 but not E-64-treated rDer p 1. The cells from mice immunized with heat-denatured rDer p 1 exhibited the highest levels of proliferation and production of IL-5 and IFN-gamma. The results indicate that the proteolytic activity of the highly purified rDer p 1 crucially commits to the sensitization process, including both IgE and IgG responses. Additionally, we demonstrated immunogenic differences by functional or structural manipulations of the rDer p 1. The findings have implications for sensitization to this relevant allergen in humans and for the design of modified allergen-vaccines for future allergen-specific immunotherapy.

Binding site-based classification of coronaviral papain-like proteases

The coronavirus replicase gene encodes one or two papain-like proteases (termed PL1pro and PL2pro) implicated in the N-terminal processing of the replicase polyprotein and thus contributing to the formation of the viral replicase complex that mediates genome replication. Using consensus fold recognition with the 3D-JURY meta-predictor followed by model building and refinement, we developed a structural model for the single PLpro present in the severe acute respiratory syndrome coronavirus (SCoV) genome, based on significant structural relationships to the catalytic core domain of HAUSP, a ubiquitin-specific protease (USP). By combining the SCoV PLpro model with comparative sequence analyses we show that all currently known coronaviral PLpros can be classified into two groups according to their binding site architectures. One group includes all PL2pros and some of the PL1pros, which are characterized by a restricted USP-like binding site. This group is designated the R-group. The remaining PL1pros from some of the coronaviruses form the other group, featuring a more open papain-like binding site, and is referred to as the O-group. This two-group, binding site-based classification is consistent with experimental data accumulated to date for the specificity of PLpro-mediated polyprotein processing and PLpro inhibition. It also provides an independent evaluation of the similarity-based annotation of PLpro-mediated cleavage sites, as well as a basis for comparison with previous groupings based on phylogenetic analyses.

Structure activity relationships of 5-, 6-, and 7-methyl-substituted azepan-3-one cathepsin K inhibitors

The syntheses, in vitro characterizations, and rat and monkey in vivo pharmacokinetic profiles of a series of 5-, 6-, and 7-methyl-substituted azepanone-based cathepsin K inhibitors are described. Depending on the particular regiochemical substitution and stereochemical configuration, methyl-substituted azepanones were identified that had widely varied cathepsin K inhibitory potency as well as pharmacokinetic properties compared to the 4S-parent azepanone analogue, 1 (human cathepsin K, K(i,app) = 0.16 nM, rat oral bioavailability = 42%, rat in vivo clearance = 49.2 mL/min/kg). Of particular note, the 4S-7-cis-methylazepanone analogue, 10, had a K(i,app) = 0.041 nM vs human cathepsin K and 89% oral bioavailability and an in vivo clearance rate of 19.5 mL/min/kg in the rat. Hypotheses that rationalize some of the observed characteristics of these closely related analogues have been made using X-ray crystallography and conformational analysis. These examples demonstrate the potential for modulation of pharmacological properties of cathepsin inhibitors by substituting the azepanone core. The high potency for inhibition of cathepsin K coupled with the favorable rat and monkey pharmacokinetic characteristics of compound 10, also known as SB-462795 or relacatib, has made it the subject of considerable in vivo evaluation for safety and efficacy as an inhibitor of excessive bone resorption in rat, monkey, and human studies, which will be reported elsewhere.

3D-QSAR and docking studies of aldehyde inhibitors of human cathepsin K

In order to better understand the structural and chemical features of human cathepsin K (CatK), which is an important cysteine protease in the pathogenesis of osteoporosis, the 3D-QSAR (CoMFA) studies were conducted on recently explored aldehyde compounds with known CatK inhibitory activities. The genetic algorithm of GOLD2.2 has been employed to position 59 aldehyde compounds into the active sites of CatK to determine the probable binding conformation. Good correlations between the predicted binding free energies and the experimental inhibitory activities suggested that the identified binding conformations of these potential inhibitors are reliable. The docking results also provided a reliable conformational alignment scheme for 3D-QSAR model. Based on the docking conformations, highly predictive comparative molecular field analysis (CoMFA) was performed with q2 value of 0.723. The predictive ability was validated by some compounds that were not included in the training set. Furthermore, the CoMFA model was mapped back to the binding sites of CatK, to get a better understanding of vital interactions between the aldehyde compounds and the protease. The CoMFA field distributions are in good agreement with the structural characteristics of the binding groove of the CatK, which suggested that the n-Bu in R4 position is the favor group substitute at P1 and moderate groups in R2 group are required on P2 substitute. In addition, 3D-QSAR results also demonstrated that aldehyde is an important pharmacophore because of electrostatic effect. These results, together with the good correlations between the inhibitory activities and the binding free energies predicted by GOLD2.2, demonstrated the power of combining docking/QSAR approach to explore the probable binding conformations of compounds at the active sites of the protein target, and further provided useful information in understanding the structural and chemical features of CatK in designing and finding new potential inhibitors.

Azepanone-Based Inhibitors of Human Cathepsin L

The extension of a previously reported cathepsin K azepanone-based inhibitor template to the design and synthesis of potent and selective inhibitors of the homologous cysteine protease cathepsin L is detailed. Structure-activity studies examining the effect of inhibitor selectivity as a function of the P3 and P2 binding elements of the potent cathepsin K inhibitor 1 revealed that incorporation of either a P3 quinoline-8-carboxamide or a naphthylene-1-carboxamide led to increased selectivity for cathepsin L over cathepsin K. Substitution of the P2 leucine of 1 with either a phenylalanine or a beta-naphthylalanine also resulted in an increased selectivity for cathepsin L over cathepsin K. Molecular modeling studies with the inhibitors docked within the active sites of both cathepsins L and K have rationalized the observed selectivities. Optimization of cathepsin L binding by the combination of the P3 naphthylene-1-carboxamide with the P2 beta-naphthylalanine provided 15, which is a potent, selective, and competitive inhibitor of human cathepsin L with a K(i) = 0.43 nM.

Inhibition of cathepsin K - a novel approach to antiresorptive therapy

Cathepsin K is a cysteine protease enzyme with high and specific expression in osteoclasts, the cells responsible for resorbing bone. Its discovery and the numerous lines of evidence that support a specialised role for cathepsin K in bone resorption are discussed. Inhibition of this enzyme would be anticipated to prevent bone resorption, and be of potential utility in diseases involving excess bone loss, such as osteoporosis. The pharmaceutical companies known to be involved in developing inhibitors of cathepsin K, and their progress to date, are described.

Comprehensive search for cysteine cathepsins in the human genome

Our study was aimed at examinating whether or not the human genome encodes for previously unreported cysteine cathepsins. To this end, we used analyses of the genome sequence and mRNA expression levels. The program TBLASTN was employed to scan the draft sequence of the human genome for the 11 known cysteine cathepsins. The cathepsin-like segments in the genome were inspected, filtered, and annotated. In addition to the known cysteine cathepsins, the scan identified three pseudogenes, closely related to cathepsin L, on chromosome 10, as well as two remote homologs, tubulointerstitial protein antigen and tubulointerstitial protein antigen-related protein. No new members of the family were identified. mRNA expression profiles for 10 known human cysteine cathepsins showed varying expression levels in 46 different human tissues and cell lines. No expression of any of the three cathepsin L-like pseudogenes was found. Based on these results, it is likely that to date all human cysteine cathepsins are known.

Dual concentration-dependent activity of thyroglobulin type-1 domain of testican: specific inhibitor and substrate of cathepsin L

The thyroglobulin type-1 (Tg-1) domain is a protein module that occurs in a variety of secreted and membrane proteins and is recognised as a potent inhibitor of cysteine peptidases. We present here some properties of the Tg-1 domain of human testican, a modularly organised proteoglycan secreted mainly by brain cells, the exactin vivofunction of which is not yet clear. The domain was prepared as a recombinant protein in aPichia pastorisexpression system and its activity was demonstrated by specific and selective inhibition of cathepsin L (Ki=0.14 nM). Interaction at high enzyme and inhibitor concentrations resulted in degradation of the domain by cathepsin L, which was not observed under conditions used for the determination of kinetic parameters. No inhibitory activity could be detected for cathepsin K, but it exhibited a very similar degradation pattern. Homology modelling provided a good explanation for the different behaviour observed with the two enzymes. Firstly, the steric fit between the interfaces of testican domain and cathepsin L is stabilised by numerous favourable forces, while no such interactions are evident in the complex with cathepsin K, and repulsive interactions even prevent access of the domain to the active site of papain. Secondly, the prolonged first loop of the domain occupies a position near the catalytic cysteine residue in a more substrate-like manner, enabling cleavage of the Gly22-Ala23bond.

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.

Crystal Structures of Calpain–E64 and –Leupeptin Inhibitor Complexes Reveal Mobile Loops Gating the Active Site

The endogenous calpain inhibitor, calpastatin, modulates some patho-physiological aspects of calpain signaling. Excess calpain can escape this inhibition and as well, many calpain isoforms and autolytically generated protease core fragments are not inhibited by calpastatin. There is a need, therefore, to develop specific, cell-permeable calpain inhibitors to block uncontrolled proteolysis and prevent tissue damage during brain and heart ischemia, spinal-cord injury and Alzheimer's diseases. Here, we report the first high-resolution crystal structures of rat mu-calpain protease core complexed with two traditional, low molecular mass inhibitors, leupeptin and E64. These structures show that access to a slightly deeper, but otherwise papain-like active site is gated by two flexible loops. These loops are divergent among the calpain isoforms giving a potential structural basis for substrate/inhibitor selectivity over other papain-like cysteine proteases and between members of the calpain family.

Structural Characterization of Vivapain-2 and Vivapain-3, Cysteine Proteases fromPlasmodium vivax: Comparative Protein Modeling and Docking Studies

Malaria remains one of the most important infectious diseases in the world. Plasmodial cysteine proteases are proposed to be promising targets for novel antimalarial drug design. Vivapain-2 and vivapain-3 are cysteine proteases from Plasmodium vivax and apparent orthologs of falcipain-2 and falcipain-3 from Plasmodium falciparum. Model structures of vivapain-2 and vivapain-3 have been derived using the comparative protein modeling approach and validated by various structure/geometry verification tools. Correlation between the interaction energies calculated based on the docking studies of the inhibitors and the corresponding association constants (k(ass)) provide additional validation for the structures. Moreover, some of the biochemical differences observed between the vivapains may be explained by the results of the docking studies. The overall structures of the two vivapains are similar to each other as well as to the falcipains with most of the catalytic residues conserved. At the same time, some important differences are observed between the sizes of the binding pockets as well as some of the residues involved in binding. The study suggests a likelihood of developing common inhibitors for these enzymes provided the interesting differences in the binding pockets of these enzymes are critically considered during such an attempt. The results of the current study can be utilized in de novo drug design for effective treatment of malaria.

Reaction mechanism of caspases: insights from QM/MM Car-Parrinello simulations

Caspases are fundamental targets for pharmaceutical interventions in a variety of diseases involving disregulated apoptosis. Here, we present a quantum mechanics/molecular mechanics Car-Parrinello study of key steps of the enzymatic reaction for a representative member of this family, caspase-3. The hydrolysis of the acyl-enzyme complex is described at the density functional (BLYP) level of theory while the protein frame and solvent are treated using the GROMOS96 force field. These calculations show that the attack of the hydrolytic water molecule implies an activation free energy of ca. DeltaF(A) approximately equal 19 +/- 4 kcal/mol in good agreement with experimental data and leads to a previously unrecognized gem-diol intermediate that can readily (DeltaF(A) approximately equal 5 +/- 3 kcal/mol) evolve to the enzyme products. Our findings assist in elucidating the striking difference in catalytic activity between caspases and other structurally well-characterized cysteine proteases (papains and cathepsins) and may help design novel transition-state analog inhibitors.

Molecular dynamics studies of caspase-3

Caspase-3 is a fundamental target for pharmaceutical interventions against a variety of diseases involving disregulated apoptosis. The enzyme is active as a dimer with two symmetry-related active sites, each featuring a Cys-His catalytic dyad and a selectivity loop, which recognizes the characteristic DEVD pattern of the substrate. Here, a molecular dynamics study of the enzyme in complex with two pentapeptide substrates DEVDG is presented, which provides a characterization of the dynamic properties of the active form in aqueous solution. The mobility of the substrate and that of the catalytic residues are rather low indicating a distinct preorganization effect of the Michaelis complex. An essential mode analysis permits us to identify coupled motions between the two monomers. In particular, it is found that the motions of the two active site loops are correlated and tend to steer the substrate toward the reactive center, suggesting that dimerization has a distinct effect on the dynamic properties of the active site regions. The selectivity loop of one monomer turns out to be correlated with the N-terminal region of the p12 subunit of the other monomer, an interaction that is also found to play a fundamental role in the electrostatic stabilization of the quaternary structure. To further characterize the specific influence of dimerization on the enzyme essential motions, a molecular dynamics analysis is also performed on the isolated monomer.

Probing the structure of falcipain-3, a cysteine protease from Plasmodium falciparum: comparative protein modeling and docking studies

Increasing resistance of malaria parasites to conventional antimalarial drugs is an important factor contributing to the persistence of the disease as a major health threat. The ongoing search for novel targets has resulted in identification and expression of several enzymes including cysteine proteases that are implicated in hemoglobin degradation. Falcipain-2 and falcipain-3 are considered to be the two principal cysteine proteases in this degradation, and hence, are potential drug targets. A homology model of falcipain-3 was built and validated by various structure/geometry verification tools as well as docking studies of known substrates. The correlation coefficient of 0.975 between interaction energies and K(m) values of these substrates provided additional support for the model. On comparison with the previously reported falcipain-2 homology model, the currently constructed falcipain-3 structure showed important differences between the S2 pockets that might explain the variations in the K(m) values of various substrates for these enzymes. Further, docking studies also provided insight into possible binding modes and interactions of ligands with falcipain-3. Results of the current study could be employed in de novo drug design leading to development of new antimalarial agents.

Functional sites in protein families uncovered via an objective and automated graph theoretic approach

We report a method for detection of recurring side-chain patterns (DRESPAT) using an unbiased and automated graph theoretic approach. We first list all structural patterns as sub-graphs where the protein is represented as a graph. The patterns from proteins are compared pair-wise to detect patterns common to a protein pair based on content and geometry criteria. The recurring pattern is then detected using an automated search algorithm from the all-against-all pair-wise comparison data of proteins. Intra-protein pattern comparison data are used to enable detection of patterns recurring within a protein. A method has been proposed for empirical calculation of statistical significance of recurring pattern. The method was tested on 17 protein sets of varying size, composed of non-redundant representatives from SCOP superfamilies. Recurring patterns in serine proteases, cysteine proteases, lipases, cupredoxin, ferredoxin, ferritin, cytochrome c, aspartoyl proteases, peroxidases, phospholipase A2, endonuclease, SH3 domain, EF-hand and lectins show additional residues conserved in the vicinity of the known functional sites. On the basis of the recurring patterns in ferritin, EF-hand and lectins, we could separate proteins or domains that are structurally similar yet different in metal ion-binding characteristics. In addition, novel recurring patterns were observed in glutathione-S-transferase, phospholipase A2 and ferredoxin with potential structural/functional roles. The results are discussed in relation to the known functional sites in each family. Between 2000 and 50,000 patterns were enumerated from each protein with between ten and 500 patterns detected as common to an evolutionarily related protein pair. Our results show that unbiased extraction of functional site pattern is not feasible from an evolutionarily related protein pair but is feasible from protein sets comprising five or more proteins. The DRESPAT method does not require a user-defined pattern, size or location of the pattern and therefore, has the potential to uncover new functional sites in protein families.

Small molecule affinity fingerprinting. A tool for enzyme family subclassification, target identification, and inhibitor design

Classifying proteins into functionally distinct families based only on primary sequence information remains a difficult task. We describe here a method to generate a large data set of small molecule affinity fingerprints for a group of closely related enzymes, the papain family of cysteine proteases. Binding data was generated for a library of inhibitors based on the ability of each compound to block active-site labeling of the target proteases by a covalent activity based probe (ABP). Clustering algorithms were used to automatically classify a reference group of proteases into subfamilies based on their small molecule affinity fingerprints. This approach was also used to identify cysteine protease targets modified by the ABP in complex proteomes by direct comparison of target affinity fingerprints with those of the reference library of proteases. Finally, experimental data were used to guide the development of a computational method that predicts small molecule inhibitors based on reported crystal structures. This method could ultimately be used with large enzyme families to aid in the design of selective inhibitors of targets based on limited structural/function information.

Preparation of active recombinant cathepsin K expressed in bacteria as inclusion body

Human cathepsin K (EC 3.4.22.38) is a member of the cysteine protease family with high primary sequence homology to cathepsins S, L, and B. It has been shown that cathepsin K plays a major role in the resorption of the bone matrix by osteoclasts. Cathepsin K has a potential as a drug target for the diseases related to bone matrix metabolism such as osteoporosis. We have expressed recombinant human procathepsin K in Escherichia coli as inclusion bodies. Purified procathepsin K had size of 38kDa which is in agreement with the predicted mass of the construct. Refolding was done by rapid dilution into 50mM Tris-HCl, pH 8.0 buffer containing 5mM EDTA, 10 mM GSH, 1mM GSSG, 0.7 M L-arginine, 0.5 M NaCl, and 1% CHAPS and further dialysis against 25 mM Tris-HCl, pH 8.0 containing 0.5 M NaCl. Mature active cathepsin K was prepared from refolded procathepsin K by incubating at 40 degrees C in pH 4.0 buffers with or without pepsin or cysteine. The presence of pepsin or cysteine in autocatalysis buffer did not have effect on the degree of conversion of nascent to mature cathepsin K, but reduced the autocatalysis time slightly. Proteolytic activity was confirmed using synthetic substrate, and Western blotting identified mature cathepsin K. Active cathepsin K had type I and II collagenolytic activities which could be inhibited by E-64.