Quantitative evaluation of each catalytic subsite of cathepsin B for inhibitory activity based on inhibitory activity-binding mode relationship of epoxysuccinyl inhibitors by X-ray crystal structure analyses of complexes

Nature-Inspired Gallinamides Are Potent Antischistosomal Agents: Inhibition of the Cathepsin B1 Protease Target and Binding Mode Analysis

Schistosomiasis, caused by a parasitic blood fluke of the genus Schistosoma, is a global health problem for which new chemotherapeutic options are needed. We explored the scaffold of gallinamide A, a natural peptidic metabolite of marine cyanobacteria that has previously been shown to inhibit cathepsin L-type proteases. We screened a library of 19 synthetic gallinamide A analogs and identified nanomolar inhibitors of the cathepsin B-type protease SmCB1, which is a drug target for the treatment of schistosomiasis mansoni. Against cultured S. mansoni schistosomula and adult worms, many of the gallinamides generated a range of deleterious phenotypic responses. Imaging with a fluorescent-activity-based probe derived from gallinamide A demonstrated that SmCB1 is the primary target for gallinamides in the parasite. Furthermore, we solved the high-resolution crystal structures of SmCB1 in complex with gallinamide A and its two analogs and describe the acrylamide covalent warhead and binding mode in the active site. Quantum chemical calculations evaluated the contribution of individual positions in the peptidomimetic scaffold to the inhibition of the target and demonstrated the importance of the P1' and P2 positions. Our study introduces gallinamides as a powerful chemotype that can be exploited for the development of novel antischistosomal chemotherapeutics.

Cathepsin B: Active site mapping with peptidic substrates and inhibitors

The potential of papain-like cysteine proteases, such as cathepsin B, as drug discovery targets for systemic human diseases has prevailed over the past years. The development of potent and selective low-molecular cathepsin B inhibitors relies on the detailed expertise on preferred amino acid and inhibitor residues interacting with the corresponding specificity pockets of cathepsin B. Such knowledge might be obtained by mapping the active site of the protease with combinatorial libraries of peptidic substrates and peptidomimetic inhibitors. This review, for the first time, summarizes a wide spectrum of active site mapping approaches. It considers relevant X-ray crystallographic data and discloses propensities towards favorable protein-ligand interactions in case of the therapeutically relevant protease cathepsin B.

Bioinformatics characterization of a cathepsin B transcript from the giant river prawn, Macrobrachium rosenbergii: Homology modeling and expression analysis after Aeromonas hydrophila infection

Cathepsin B is a lysosomal proteolytic enzyme that has been suggested to play a role in pathological processes of immune system. In this study, the full-length cDNA sequence of cathepsin B transcript in the giant river prawn Macrobrachium rosenbergii (MrCTSB) was obtained from 454 pyrosequencing of cDNAs from hepatopancreas and muscle. It was 1158 bp in length, containing an open reading frame (ORF) of 987 bp corresponding to 328 amino acids. The predicted molecular mass and pI of MrCTSB protein was 36.04 kDa and 4.73. The major characteristics of MrCTSB protein consisted of a propeptide of C1 peptidase family at the N-terminus and a cysteine protease (Pept_C1) domain at the C-terminus. The 3-dimentional structure of MrCTSB was constructed by computer-assisted homology modeling. The folding of MrCTSB was highly conserved to human CTSB structure and the modeled MrCTSB displayed characteristics of cysteine proteinases superfamily. The docking study was performed to investigate binding interactions between known inhibitors against MrCTSB. Known inhibitors were oriented in the groove of catalytic site cleft. They bound to subsites from S2, S1, S1', and S2', respectively, with key residues in each subsite. Challenge of juvenile prawns with Aeromonas hydrophila revealed that the MrCTSB transcript in hepatopancreas significantly increased at 60-96 h post injection (hpi). This suggested that MrCTSB may play roles in innate immunity of M. rosenbergii. Our results provide useful information for a more comprehensive study in immune-related functions of MrCTSB.

Crystallographic and docking (Cathepsins B, K, L and S) studies on bioactive halotelluroxetanes

The molecular structures of the halotelluroxetanes p-MeOC6H4Te(X)[C(=C(H)X′)C(CH2)nO], X=X′=Cl and n=6 (1) and X=Cl, X′=Br and n=5 (4), show similar binuclear aggregates sustained by {· · ·Te–O}2 cores comprising covalent Te–O and secondary Te· · ·O interactions. The resulting C2ClO2(lone-pair) sets define pseudo-octahedral geometries. In each structure, C–X· · ·π(arene) interactions lead to supramolecular layers. Literature studies have shown these and related compounds (i.e. 2: X=X′=Cl and n=5; 3: X=X′=Br and n=5) to inhibit Cathepsins B, K, L and S to varying extents. Molecular docking calculations have been conducted on ligands (i.e. cations derived by removal of the tellurium-bound X atoms) 1′–3′ (note 3′=4′) enabling correlations between affinity for sub-sites and inhibition. The common feature of all docked complexes was the formation of a Te–S covalent bond with cysteine residues, the relative stability of the ligands with an E-configuration and the formation of a C–O· · ·π interaction with the phenyl ring; for 1′ the Te–S covalent bond was weak, a result correlating with its low inhibition profile. At the next level differences are apparent, especially with respect to the interactions formed by the organic-ligand-bound halides. While these atoms do not form specific interactions in Cathepsins B and K, in Cathepsin L, these halides are involved in C–O· · ·X halogen bonds.

Design, Synthesis, and Structure-Activity Relationship Study of Epoxysuccinyl-Peptide Derivatives as Cathepsin B Inhibitors

Cathepsin B is a lysosomal cysteine protease involved in many diseases. The present research demonstrates that derivatives of epoxysuccinyl-peptide are effective and selective cathepsin B inhibitors. We synthesized a series of epoxysuccinyl-peptide derivatives based on the well-known cathepsin B inhibitor E64d. Specifically, we substituted the 2-methylpropane group at the R¹ position of E64d with a sulfane, such as ethyl(methyl) sulfane or benzyl(methyl) sulfane. We also designed and synthesized a library of molecules with various substituents at the R² position of E64d to replace 2-methylbutane. By studying the structure-activity relationships of these newly synthesized molecules as cathepsin B inhibitors, we demonstrated that substituting ethyl(methyl) sulfane for 2-methylbutane (R²) of E64d improves the inhibitory activity and selectivity for cathepsin B inhibition. Our new cathepsin B inhibitors were highly effective and selective.

Crystallographic, DFT and docking (cathepsin B) studies on an organotellurium(IV) compound

Some biologically active organotellurium compounds exhibit inhibitory potency against cathepsin B. In this study, an alkyl derivative, viz. [CH3(CH2)2C(I)=C(H)](nBu)TeI2, 1, has been structurally characterised by X-ray crystallography and shown to be coordinated within a C2I2 donor set. When the stereochemically active lone pair of electrons is taken into account, a distorted trigonal bipyramidal geometry results with the iodide atoms in axial positions. Both intra- and inter-molecular Te···I interactions are also noted. If all interactions are considered, the coordination geometry is based on a Ψ-pentagonal bipyramidal geometry. An unusual feature of the structure is the curving of the functionalised C5 chain. This feature has been explored by DFT methods and shown to arise as a result of close C–H···I interactions. A docking study (cathepsin B) was performed to understand the inhibition mechanism and to compare the new results with previous observations. Notably, 1 has the same pose exhibited by analogous biologically active compounds with aryl groups. Thus, the present study suggests that (alkyl)2TeX2 compounds should also be evaluated for biological activity.

Development of cell-active non-peptidyl inhibitors of cysteine cathepsins

Cysteine cathepsins are an important class of enzymes that coordinate a variety of important cellular processes, and are implicated in various types of human diseases. However, small molecule inhibitors that are cell-permeable and non-peptidyl in nature are scarcely available. Herein the synthesis and development of sulfonyloxiranes as covalent inhibitors of cysteine cathepsins are reported. From a library of compounds, compound 5 is identified as a selective inhibitor of cysteine cathepsins. Live cell imaging and immunocytochemistry of metastatic human breast carcinoma MDA-MB-231 cells document the efficacy of compound 5 in inhibiting cysteine cathepsin activity in living cells. A cell-motility assay demonstrates that compound 5 is effective in mitigating the cell-migratory potential of highly metastatic breast carcinoma MDA-MB-231 cells.

Development of new cathepsin B inhibitors: combining bioisosteric replacements and structure-based design to explore the structure-activity relationships of nitroxoline derivatives

Human cathepsin B has many house-keeping functions, such as protein turnover in lysosomes. However, dysregulation of its activity is associated with numerous diseases, including cancers. We present here the structure-based design and synthesis of new cathepsin B inhibitors using the cocrystal structure of 5-nitro-8-hydroxyquinoline in the cathepsin B active site. A focused library of over 50 compounds was prepared by modifying positions 5, 7, and 8 of the parent compound nitroxoline. The kinetic parameters and modes of inhibition were characterized, and the selectivities of the most promising inhibitors were determined. The best performing inhibitor 17 was effective in cell-based in vitro models of tumor invasion, where it significantly abrogated invasion of MCF-10A neoT cells. These data show that we have successfully explored the structure-activity relationships of nitroxoline derivatives to provide new inhibitors that could eventually lead to compounds with clinical usefulness against the deleterious effects of cathepsin B in cancer progression.

Vitamin K2 covalently binds to Bak and induces Bak-mediated apoptosis

Vitamin K2 (VK2, menaquinone) is known to have anticancer activity in vitro and in vivo. Although its effect is thought to be mediated, at least in part, by the induction of apoptosis, the underlying molecular mechanism remains elusive. Here, we identified Bcl-2 antagonist killer 1 (Bak) as a molecular target of VK2-induced apoptosis. VK2 directly interacts with Bak and induces mitochondrial-mediated apoptosis. Although Bak and Bcl-2-associated X protein (Bax), another member of the Bcl-2 family, are generally thought to be functionally redundant, only Bak is necessary and sufficient for VK2-induced cytochrome c (cyt c) release and cell death. Moreover, VK2-2,3 epoxide, an intracellular metabolite of VK2, was shown to covalently bind to the cysteine-166 residue of Bak. Several lines of evidence suggested that the covalent attachment of VK2 is critical for apoptosis induction. Thus this study reveals a specific role for Bak in mitochondria-mediated apoptosis. This study also provides insight into the anticancer effects of VK2 and suggests that Bak may be a potential target of cancer therapy.

[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.

Structural basis for inhibition of cathepsin B drug target from the human blood fluke, Schistosoma mansoni

Schistosomiasis caused by a parasitic blood fluke of the genus Schistosoma afflicts over 200 million people worldwide. Schistosoma mansoni cathepsin B1 (SmCB1) is a gut-associated peptidase that digests host blood proteins as a source of nutrients. It is under investigation as a drug target. To further this goal, we report three crystal structures of SmCB1 complexed with peptidomimetic inhibitors as follows: the epoxide CA074 at 1.3 Å resolution and the vinyl sulfones K11017 and K11777 at 1.8 and 2.5 Å resolutions, respectively. Interactions of the inhibitors with the subsites of the active-site cleft were evaluated by quantum chemical calculations. These data and inhibition profiling with a panel of vinyl sulfone derivatives identify key binding interactions and provide insight into the specificity of SmCB1 inhibition. Furthermore, hydrolysis profiling of SmCB1 using synthetic peptides and the natural substrate hemoglobin revealed that carboxydipeptidase activity predominates over endopeptidolysis, thereby demonstrating the contribution of the occluding loop that restricts access to the active-site cleft. Critically, the severity of phenotypes induced in the parasite by vinyl sulfone inhibitors correlated with enzyme inhibition, providing support that SmCB1 is a valuable drug target. The present structure and inhibitor interaction data provide a footing for the rational design of anti-schistosomal inhibitors.

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.

Multi-conformation 3D QSAR study of benzenesulfonyl-pyrazol-ester compounds and their analogs as cathepsin B inhibitors

Cathepsin B has been found being responsible for many human diseases. Inhibitors of cathepsin B, a ubiquitous lysosomal cysteine protease, have been developed as a promising treatment for human diseases resulting from malfunction and over-expression of this enzyme. Through a high throughput screening assay, a set of compounds were found able to inhibit the enzymatic activity of cathepsin B. The binding structures of these active compounds were modeled through docking simulation. Three-dimensional (3D) quantitative structure-activity relationship (QSAR) models were constructed using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) based on the docked structures of the compounds. Strong correlations were obtained for both CoMFA and CoMSIA models with cross-validated correlation coefficients (q²) of 0.605 and 0.605 and the regression correlation coefficients (r²) of 0.999 and 0.997, respectively. The robustness of these models was further validated using leave-one-out (LOO) method and training-test set method. The activities of eight (8) randomly selected compounds were predicted using models built from training set of compounds with prediction errors of less than 1 unit for most compounds in CoMFA and CoMSIA models. Structural features for compounds with improved activity are suggested based on the analysis of the CoMFA and CoMSIA contour maps and the property map of the protein ligand binding site. These results may help to provide better understanding of the structure-activity relationship of cathepsin B inhibitors and to facilitate lead optimization and novel inhibitor design. The multi-conformation method to build 3D QSAR is very effective approach to obtain satisfactory models with high correlation with experimental results and high prediction power for unknown compounds.

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.

QSAR models for predicting cathepsin B inhibition by small molecules--continuous and binary QSAR models to classify cathepsin B inhibition activities of small molecules

Cathepsin B is a potential target for the development of drugs to treat several important human diseases. A number of inhibitors targeting this protein have been developed in the past several years. Recently, a group of small molecules were identified to have inhibitory activity against cathepsin B through high throughput screening (HTS) tests. In this study, traditional continuous and binary QSAR models were built to classify the biological activities of previously identified compounds and to distinguish active compounds from inactive compounds for drug development based on the calculated molecular and physicochemical properties. Strong correlations were obtained for the continuous QSAR models with regression correlation coefficients (r²) and cross-validated correlation coefficients (q²) of 0.77 and 0.61 for all compounds, and 0.82 and 0.68 for the compound set excluding 3 outliers, respectively. The models were further validated through the leave-one-out (LOO) method and the training-test set method. The binary models demonstrated a strong level of predictability in distinguishing the active compounds from inactive compounds with accuracies of 0.89 and 0.94 for active and inactive compounds, respectively, in non-cross-validated models. Similar results were obtained for the cross-validated models. Collectively, these results demonstrate the models’ ability to discriminate between active and inactive compounds, suggesting that the models may be used to pre-screen compounds to facilitate compound optimization and to design novel inhibitors for drug development.

Computational analysis of the cathepsin B inhibitors activities through LR-MMPBSA binding affinity calculation based on docked complex

Cathepsin B, a ubiquitous lysosomal cysteine protease, is involved in many biological processes related to several human diseases. Inhibitors targeting the enzyme have been investigated as possible diseases treatments. A set of 37 compounds were recently found active in a high throughput screening assay to inhibit the catalytic activity of Cathepsin B, with chemical structures and biological test results available to the public in the PubChem BioAssay Database (AID 820). In this study, we compare these experimental activities to the results of theoretical predictions from binding affinity calculation with a LR-MM-PNSA approach based on docked complexes. Strong correlations (r(2) = 0.919 and q(2) = 0.887 for the best) are observed between the theoretical predictions and experimental biological activity. The models are cross-validated by four independent predictive experiments with randomly split compounds into training and test sets. Our results also show that the results based on protein dimer show better correlations with experimental activity when compared to results based on monomer in the in silico calculations.

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.

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.