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.     

Protein kinase C-alpha and -delta are required for NADPH oxidase activation in WKYMVm-stimulated IMR90 human fibroblasts

Gene duplications in rodents have given rise to a family of proteases that are expressed exclusively in placenta. To define the biological role of these enzymes specific inhibitors are needed to differentiate their activities from other more ubiquitously expressed proteases, such as cathepsins B and L. Libraries of peptidyl inhibitors based upon a 4-cyclohexanone pharmacophore were screened for inhibition of cathepsins P, L, and B. The tightest binding dipeptidyl inhibitor for cathepsin P contained Tyr in P(2) and Trp in P(2)('), consistent with the specificity of this enzyme for hydrophobic amino acids at these sites in synthetic substrates. An inhibitor containing Trp in both P(2) and P(2)(') provided better discrimination between cathepsin P and cathepsins B and L. Extension of the inhibitors to include P(3), and P(3)(') amino acids identified an inhibitor with Trp in P(2), P(2)('), and P(3), and Phe in P(3)(') that bound to cathepsin P with a K(i) of 32 nM. This specificity for inhibitors with hydrophobic aromatic amino acids in these four positions is unique among the lysosomal cysteine proteases. This inhibitor bound to cathepsin P an order of magnitude tighter than to mouse and human cathepsin L and two orders of magnitude tighter than to human cathepsin B. Cbz-Trp-Trp-4-cyclohexanone-Trp-Phe-OMe can discriminate cathepsin P from cathepsins B and L and consequently can be used to specifically inhibit and identify cathepsin P in cellular systems.     

Comparative study on specificities of rat cathepsin L and papain: amino acid differences at substrate-binding sites are involved in their specificities

Sixty-nine rat cathepsin L-susceptible peptide bonds were analyzed employing various peptide substrates. The proteolytic specificities of rat cathepsin L and papain were compared and the results are discussed in relation to differences in amino acid residues around their binding sites. The specificity of cathepsin L, which is characterized by a remarkable preference for hydrophobic amino acids at the P2 site of the scissile peptide bonds, was analogous to that of papain as a whole. This analogous specificity suggests that the binding sites of the two proteases are analogous, as expected from their homologous amino acid sequences. However, there is a slight difference in the preference for S3 site between them. That is, cathepsin L showed a greater preference for bulky and hydrophobic amino acids at the S3 site than did papain. Based on the computer-graphically deduced structure of the binding sites of cathepsin L, the preferences for hydrophobic amino acids at the S2 site and for bulky and hydrophobic amino acids at the S3 site of the protease are supposed to be related to the compensating amino acid substitutions at the S2 site (V133A and V157L) and the reduction in size at the S3 site (Y61Q and Y67L), respectively. The discussion of the effect of the amino acid substitutions on the proteolytic activities of cathepsin L and papain in this paper provides a basis for more advanced studies of the relationship between structure and function of proteases belonging to the papain superfamily by means of protein engineering.     

Characterization of amino acid variation at strategic positions in parasite and human proteases for selective inhibition of falcipains in Plasmodium falciparum

Falcipains (FP) of Plasmodium falciparum are important virulence factors marked as potential targets for antimalarial drug discovery. In this study, the previously uncharacterized fp2B (PF11_0161) was shown to be highly expressed as an active enzyme during the erythrocytic stage. With three related proteases in the FP family and the existence of human homologues, it is prudent to identify clusters of residues unique to the parasite proteases that can be targeted selectively for drug design. Using bioinformatic tools, we have carefully mapped out a highly conserved and unique region constituted by I85, S149, and A151 in the plasmodial proteases that can influence the development of compounds capable of inhibiting the entire FP family. Taking drug interactions with the human homologues into consideration, these residues in FP2B were replaced with the cognate residues found in human cathepsin L (catL) for evaluation. Despite the high sequence similarity between the FP2 isozymes (97.5%), FP2B is found to be more tolerant to amino acid substitution at position 149 than FP2A. This structural disparity implied that residues mediating peptide substrate interactions are not fully conserved across the FP family and warrant attention in the design and evaluation of protease inhibitors focused on the FPs. The simultaneous substitution of the neighboring residues (I85 or A151) rendered the double mutants (S149A/I85M and S149A/A151D) completely inactive. Significantly, the mutations did not result in 'catL-like' specificity, suggesting that substrate-based inhibitors could be rationally designed against these important parasite-specific structural determinants.     

Selective inhibition of the collagenolytic activity of human cathepsin K by altering its S2 subsite specificity

The primary specificity of papain-like cysteine proteases (family C1, clan CA) is determined by S2-P2 interactions. Despite the high amino acid sequence identities and structural similarities between cathepsins K and L, only cathepsin K is capable of cleaving interstitial collagens in their triple helical domains. To investigate this specificity, we have engineered the S2 pocket of human cathepsin K into a cathepsin L-like subsite. Using combinatorial fluorogenic substrate libraries, the P1-P4 substrate specificity of the cathepsin K variant, Tyr67Leu/Leu205Ala, was determined and compared with those of cathepsins K and L. The introduction of the double mutation into the S2 subsite of cathepsin K rendered the unique S2 binding preference of the protease for proline and leucine residues into a cathepsin L-like preference for bulky aromatic residues. Homology modeling and docking calculations supported the experimental findings. The cathepsin L-like S2 specificity of the mutant protein and the integrity of its catalytic site were confirmed by kinetic analysis of synthetic di- and tripeptide substrates as well as pH stability and pH activity profile studies. The loss of the ability to accept proline in the S2 binding pocket by the mutant protease completely abolished the collagenolytic activity of cathepsin K whereas its overall gelatinolytic activity remained unaffected. These results indicate that Tyr67 and Leu205 play a key role in the binding of proline residues in the S2 pocket of cathepsin K and are required for its unique collagenase activity.     

Human cathepsin V functional expression, tissue distribution, electrostatic surface potential, enzymatic characterization, and chromosomal localization

Cathepsin V, a thymus and testis-specific human cysteine protease, was expressed in Pichia pastoris, and its physicokinetic properties were determined. Recombinant procathepsin V is autocatalytically activated at acidic pH and is effectively inhibited by various cysteine protease class-specific inhibitors. The S2P2 subsite specificity of cathepsin V was found to be intermediate between those of cathepsins S and L. The substrate binding pocket, S2, accepted both aromatic and nonaromatic hydrophobic residues, whereas cathepsins L and S preferred either an aromatic or nonaromatic hydrophobic residue, respectively. In contrast to cathepsin L, but similar to cathepsin S, cathepsin V exhibited only a very weak collagenolytic activity. Furthermore, cathepsin V was determined to be significantly more stable at mildly acidic and neutral pH than cathepsin L, but distinctly less stable than cathepsin S. A homology structure model of cathepsin V revealed completely different electrostatic potentials on the molecular surface when compared with human cathepsin L. The model-based electrostatic potential of human cathepsin V was neutral to weakly positive at and in the vicinity of the active site cleft, whereas that of cathepsin L was negative over extended regions of the surface. Surprisingly, the electrostatic potential of the human cathepsin V model structure resembled that of the model structure of mouse cathepsin L. These differences in the electrostatic potential at the molecular surfaces provide a reactivity determinant that may be the source of differences in substrate selectivity and pH stability. Cathepsin V was mapped to the chromosomal region 9q22.2, a site adjacent to the cathepsin L locus. The high sequence identity and the overlapping chromosomal gene loci suggest that both proteases evolved from an ancestral cathepsin L-like precursor by gene duplication.     

Crystal structure of the parasite protease inhibitor chagasin in complex with a host target cysteine protease

Chagasin is a protein produced by Trypanosoma cruzi, the parasite that causes Chagas' disease. This small protein belongs to a recently defined family of cysteine protease inhibitors. Although resembling well-known inhibitors like the cystatins in size (110 amino acid residues) and function (they all inhibit papain-like (C1 family) proteases), it has a unique amino acid sequence and structure. We have crystallized and solved the structure of chagasin in complex with the host cysteine protease, cathepsin L, at 1.75 A resolution. An inhibitory wedge composed of three loops (L2, L4, and L6) forms a number of contacts responsible for high-affinity binding (K(i), 39 pM) to the enzyme. All three loops interact with the catalytic groove, with the central loop L2 inserted directly into the catalytic center. Loops L4 and L6 embrace the enzyme molecule from both sides and exhibit distinctly different patterns of protein-protein recognition. Comparison with a 1.7 A structure of uncomplexed chagasin, also determined in this study, demonstrates that a conformational change of the first binding loop (L4) allows extended binding to the non-primed substrate pockets of the enzyme active site cleft, thereby providing a substantial part of the inhibitory surface. The mode of chagasin binding is generally similar, albeit distinctly different in detail, when compared to those displayed by cystatins and the cysteine protease inhibitory p41 fragment of the invariant chain. The chagasin-cathepsin L complex structure provides details of how the parasite protein inhibits a host enzyme of possible importance in host defense. The high level of structural and functional similarity between cathepsin L and the T. cruzi enzyme cruzipain gives clues to how the cysteine protease activity of the parasite can be targeted. This information will aid in the development of synthetic inhibitors for use as potential drugs for the treatment of Chagas disease.     

Cathepsin D inactivates cysteine proteinase inhibitors, cystatins

The formation of inactive complexes in excess molar amounts of human cathepsins H and L with their protein inhibitors human stefin A, human stefin B and chicken cystatin at pH 5.6 has been shown by measurement of enzyme activity coupled with reverse-phase HPLC not to involve covalent cleavage of the inhibitors. Inhibition must be the direct result of binding. On the contrary the interaction of cystatins with aspartic proteinase cathepsin D at pH 3.5 for 60 min followed by HPLC resulted in their inactivation accompanied by peptide bond cleavage at several sites, preferentially those involving hydrophobic amino acid residues. The released peptides do not inhibit papain and cathepsin L. These results explain reported elevated levels of cysteine proteinases and lead to the proposal that cathepsin D exerts an important function, through inactivation of cystatins, in the increased activities of cysteine proteinases in human diseases including muscular distrophy.     

Cathepsin H functions as an aminopeptidase in secretory vesicles for production of enkephalin and galanin peptide neurotransmitters

Peptide neurotransmitters function as key intercellular signaling molecules in the nervous system. These peptides are generated in secretory vesicles from proneuropeptides by proteolytic processing at dibasic residues, followed by removal of N- and/or C-terminal basic residues to form active peptides. Enkephalin biosynthesis from proenkephalin utilizes the cysteine protease cathepsin L and the subtilisin-like prohormone convertase 2 (PC2). Cathepsin L generates peptide intermediates with N-terminal basic residue extensions, which must be removed by an aminopeptidase. In this study, we identified cathepsin H as an aminopeptidase in secretory vesicles that produces (Met)enkephalin (ME) by sequential removal of basic residues from KR-ME and KK-ME, supported by in vivo knockout of the cathepsin H gene. Localization of cathepsin H in secretory vesicles was demonstrated by immunoelectron microscopy and immunofluorescence deconvolution microscopy. Purified human cathepsin H sequentially removes N-terminal basic residues to generate ME, with peptide products characterized by nano-LC-MS/MS tandem mass spectrometry. Cathepsin H shows highest activities for cleaving N-terminal basic residues (Arg and Lys) among amino acid fluorogenic substrates. Notably, knockout of the cathepsin H gene results in reduction of ME in mouse brain. Cathepsin H deficient mice also show a substantial decrease in galanin peptide neurotransmitter levels in brain. These results illustrate a role for cathepsin H as an aminopeptidase for enkephalin and galanin peptide neurotransmitter production.     

Inhibition of cathepsin G by 2-amino-3,1-benzoxazin-4-ones: kinetic investigations and docking studies

A series of benzoxazinones was used to investigate the interaction of human cathepsin G with acyl-enzyme inhibitors. With respect to the primary specificity of cathepsin G, inhibitors with hydrophobic or basic residues at position 2 were included in the study. Parameters of the enzyme acylation and deacylation were determined by slow-binding kinetics in the presence of a chromogenic substrate. For selected inhibitors, the time course of the enzyme-catalyzed conversion of the inhibitors was followed. This approach was suitable to elucidate a rate-determining deacylation step. Docking simulations of the noncovalent enzyme-inhibitor complexes were performed and several clusters were analyzed for each inhibitor. The amino acids of the active site that participate in the binding of the inhibitors were determined. The arrangements in several clusters of an inhibitor were not uniform with respect to the orientation by which the inhibitor was bound in the S(1) pocket. Docking of the basic piperazino derivatives 6 and 10 indicated an interaction with Glu 226 at the bottom of the S(1) specificity pocket. The (N-methyl)benzylamino derivative 1 showed the strongest acylation rate (k(on)=1200 M(-1) s(-1)), which was attributed to a high extent of pseudo-productive orientations of the noncovalent preassociation complex.     

Identification and discrimination of extracellularly active cathepsins B and L in high-invasive melanoma cells

We established a novel protocol for lithium dodecyl sulfate (LDS) gelatin zymography, which operates under reducing conditions and at a slightly acidic pH value (6.5). This zymographic assay is based on polyacrylamide gel electrophoresis and facilitates the electrophoretic separation of human cathepsins in an active state. By this technique, activity of purified human liver cathepsin B was detected at a concentration as low as 50 ng and was blocked only in the presence of the cysteine protease inhibitor E-64 and the specific cathepsin B inhibitor CA-074 but not by aspartate, serine, or matrix metalloprotease inhibitors. The method was applied to analyze cathepsin activities in cell culture supernatants of the high-invasive melanoma cell line MV3. Interestingly, LDS zymography of MV3 cell supernatants in combination with specific inhibitors of cathepsins B and L identified three forms of extracellularly active cathepsin B and two forms of proteolytically active cathepsin L. We herein describe the generation and biochemical significance of acidic LDS zymography. This novel method permits not only the enzymatic analysis of purified cysteine proteases but also the identification and discrimination of different cathepsin activities in biological fluids, cell lysates, or supernatants, especially of cathepsins B and L, which are closely linked to major inflammatory and malignant processes.     

Elastin Degradation by Cathepsin V Requires Two Exosites

Cathepsin V is a highly effective elastase and has been implicated in physiological and pathological extracellular matrix degradation. However, its mechanism of action remains elusive. Whereas human cathepsin V exhibits a potent elastolytic activity, the structurally homologous cathepsin L, which shares a 78% amino acid sequence, has only a minimal proteolytic activity toward insoluble elastin. This suggests that there are distinct structural domains that play an important role in elastinolysis. In this study, a total of 11 chimeras of cathepsins V and L were generated to identify elastin-binding domains in cathepsin V. Evaluation of these chimeras revealed two exosites contributing to the elastolytic activity of cathepsin V that are distant from the active cleft of the protease and are located in surface loop regions. Replacement of exosite 1 or 2 with analogous residues from cathepsin L led to a 75 and 43% loss in the elastolytic activity, respectively. Replacement of both exosites yielded a non-elastase variant similar to that of cathepsin L. Identification of these exosites may contribute to the design of inhibitors that will only affect the elastolytic activity of cysteine cathepsins without interfering with other physiological protease functions.     

Novel cathepsin D inhibitors block the formation of hyperphosphorylated tau fragments in hippocampus

Lysosomal disturbances may be a contributing factor to Alzheimer's disease. We used novel compounds to test if suppression of the lysosomal protease cathepsin D blocks production of known precursors to neurofibrillary tangles. Partial lysosomal dysfunction was induced in cultured hippocampal slices with a selective inhibitor of cathepsins B and L. This led within 48 h to hyperphosphorylated tau protein fragments recognized by antibodies against human tangles. Potent nonpeptidic cathepsin D inhibitors developed using combinatorial chemistry and structure-based design blocked production of the fragments in a dose-dependent fashion. Threshold was in the submicromolar range, with higher concentrations producing complete suppression. The effects were selective and not accompanied by pathophysiology. Comparable results were obtained with three structurally distinct inhibitors. These results support the hypothesis that cathepsin D links lysosomal dysfunction to the etiology of Alzheimer's disease and suggest a new approach to treating the disease.     

Site-directed mutagenesis, proteolytic cleavage, and activation of human proheparanase

Heparanase is an endo-beta-D-glucuronidase that degrades heparan sulfate in the extracellular matrix and cell surfaces. Human proheparanase is produced as a latent 65-kDa polypeptide undergoing processing at two potential proteolytic cleavage sites, located at Glu109-Ser110 (site 1) and Gln157-Lys158 (site 2). Cleavage of proheparanase yields 8- and 50-kDa subunits that heterodimerize to form the active enzyme. The fate of the linker segment (Ser110-Gln157) residing between the two subunits, the mode of processing, and the protease(s) engaged in proheparanase processing are currently unknown. We applied multiple site-directed mutagenesis and deletions to study the nature of the potential cleavage sites and amino acids essential for processing of proheparanase in transfected human choriocarcinoma cells devoid of endogenous heparanase but possessing the enzymatic machinery for proper processing and activation of the proenzyme. Although mutagenesis at site 1 and its flanking sequences failed to identify critical residues for proteolytic cleavage, processing at site 2 required a bulky hydrophobic amino acid at position 156 (i.e. P2 of the cleavage site). Substitution of Tyr156 by Ala or Glu, but not Val, resulted in cleavage at an upstream site in the linker segment, yielding an improperly processed inactive enzyme. Processing of the latent 65-kDa proheparanase in transfected Jar cells was inhibited by a cell-permeable inhibitor of cathepsin L. Moreover, recombinant 65-kDa proheparanase was processed and activated by cathepsin L in a cell-free system. Altogether, these results suggest that proheparanase processing at site 2 is brought about by cathepsin L-like proteases. The involvement of other members of the cathepsin family with specificity to bulky hydrophobic residues cannot be excluded. Our results and a three-dimensional model of the enzyme are expected to accelerate the design of inhibitory molecules capable of suppressing heparanase-mediated enhancement of tumor angiogenesis and metastasis.     

Differences in local conformation in human apolipoprotein B-100 of plasma low density and very low density lipoproteins as identified by cathepsin D

The conformational changes of human apolipoprotein (apo) B-100 which accompany the conversion of plasma very low density lipoproteins (VLDL) to low density lipoproteins (LDL) were investigated by studying the accessibility of apoB-100 in LDL and VLDL to limited proteolysis with cathepsin D, an aspartyl proteinase involved in intracellular protein degradation. We characterized the proteolytic products of apoB-100 by sodium dodecyl sulfate-polyacrylamide gel electrophoresis followed by NH2-terminal sequence analysis to locate cleavage sites. The results identified at least 10 cleavage products generated from apoB-100 and showed differential accessibility of cleavage sites for cathepsin D in apoB-100 between LDL and VLDL. We identified a specific peptide region (residues 2660-2710), which is preferentially accessible to limited proteolysis by cathepsin D but inaccessible to limited proteolysis by 12 other enzymes tested. Within this peptide region, cathepsin D cleaved apoB-100 of LDL and VLDL preferentially at different sites, separated by 33-36 amino acids (2665-2666 or 2668-2669 (LDL) and 2701-2702 (VLDL]. In addition, we identified a cleavage site, located at residues 3272-3273, specific for cathepsin D, which is contained within the COOH-terminal enzyme-accessible peptide region (residues 3180-3280), which we have demonstrated using 12 endoproteases with various specificities. The previously identified NH2-terminal region (residues 1280-1320) appears to be resistant to limited cleavage by cathepsin D. However, a new site was revealed only approximately 66 kDA from the NH2 terminus. We conclude that differential accessibility and the shift of the novel scission site for cathepsin D by 33-36 amino acids indicate significant differences in local conformation at these sites in apoB-100 as VLDL are converted to LDL.     

Arginine-based structures are specific inhibitors of cathepsin C. Application of peptide combinatorial libraries.

Novel synthetic peptide inhibitors of lysosomal cysteine proteinase cathepsin C have been designed through the use of soluble peptide combinatorial libraries. The uncovered structural inhibitory module consists of the N-terminal cluster of L-arginine residues. Its modification with D-amino acids or arginine derivatives did not increase the inhibition strength. Inhibitory potency of oligoarginines improves with the elongation of peptide chain reaching a maximum for octa-L-arginine. The oligoarginines specifically interact with the cathepsin C active site as shown by competitive-type inhibition kinetics (Ki approximately 10-5 M) and intrinsic fluorescence measurements. The inhibitory interaction of oligoarginines is established through the specific spatial contact of a net of guanidino groups in the arginine side-chains, as indicated by comparison with inhibitory action of low molecular mass guanidine derivatives (Ki approximately 10-3 M). Nonarginine polyionic compounds cannot mimic the inhibitory effect of oligoarginines. The arginine-based peptide inhibitors were selective towards cathepsin C among other cysteine proteinases tested.     

Sequence analysis, tissue distribution, and expression of rat cathepsin S.

Cysteine proteases are involved in many diverse cellular processes ranging from processing of precursor proteins to intracellular degradation. In an effort to identify novel cysteine proteases, we used the polymerase chain reaction and primers directed against the catalytic sites of previously cloned cysteine proteases. From rat brain mRNA, a 600-base pair band was amplified; cloning and partial sequence analysis of this band resulted in the identification of cathepsins B and L and five novel sequences. The novel cDNAs contained a number of residues conserved in lysosomal cysteine proteases, including the active site residue His159 (papain numbering). In addition, the amino acid homology between the novel sequences and either cathepsins B, L, or H, ranged from 63 to 32%. The insert with highest homology was used to screen a rat brain cDNA library; a 1334-base pair cDNA was isolated and the nucleotide sequence determined. This sequence encodes an open reading frame of 330 amino acids which is 82% homologous to human cathepsin S, suggesting that this sequence represents rat cathepsin S. Northern blot analysis for rat cathepsin S revealed tissue-specific expression distinct from the distribution of cathepsin B and L. The regulation of expression of rat cathepsin S mRNA in response to thyroid-stimulating hormone was studied in a rat thyroid cell line FRTL-5. The level of cathepsin S mRNA was substantially increased in response to thyroid-stimulating hormone, whereas cathepsin B and cathepsin L mRNA levels were not altered by this treatment. A portion of cDNA encoding the predicted mature protein of rat cathepsin S was expressed as a glutathione S-transferase-fusion protein. The affinity-purified protein exhibited proteolytic activity with properties similar to bovine cathepsin S. Taken together, these results imply highly specific functions for cathepsin S.     

DNA binding affinity of pentapeptides selected from combinatorial library.

The combinatorial method has been applied to determine peptide ligands to the duplex DNA by using the solid-state pentapeptide library and the target-DNA conjugated magnetic beads. Seventy-one sequences were determined as ligands for AT duplex. Interestingly, hydrophobic amino acids such as Phe, Ile and Gly were most frequently determined. Relative binding affinity of the selected pentapeptides with the various DNA sequences was estimated by ethidium displacement assay in 10 mM SHE buffer. FQGII constituted of amino acids that were most frequently determined in the random screening showed highest binding affinity to the duplex DNA.     

An exploration of the primary specificity site of cathepsin B

Peptidyl diazomethyl ketones inactivate cathepsin B apparently by alkylation of the active center thiol following complex formation as in the case of benzyloxycarbonyl (Cbz)-Phe-AlaCHN2. The phenylalanine contributes considerably to binding in the secondary specificity site. In order to define the topography of the active center region comprising the primary specificity site of beef spleen cathepsin B, a series of peptidyl diazomethyl ketones having the general structure Cbz-Phe-X-CHN2 has now been synthesized. The amino acid, X, has been varied in size to include rather large side chains which might reveal available binding potential or limitations. Some of the reagents, in fact, were not inhibitory even at 10(-4) M. Others, however, that did measurably inactivate cathepsin B provided a range of reactivities that extended over 5 orders of magnitude and correlated with affinity in the reversible phase of inactivation. Some large side chains, for example, that of tryptophan, were very poorly tolerated in this region of the active center, whereas others, such as O-benzyl threonine, provided remarkably active inhibitors. A topographical rationalization of the results is offered.     

Structural and functional conservation profiles of novel cathepsin L-like proteins identified in the Drosophila melanogaster genome

Cathepsin L is a cysteine protease which degrades connective tissue proteins including collagen, elastin, and fibronectin. In this study, five well-characterized cathepsin L proteins from different arthropods were used as query sequences for the Drosophila genome database. The search yielded 10 cathepsin L-like sequences, of which eight putatively represent novel cathepsin L-like proteins. To understand the phylogenetic relationship among these cathepsin L-like proteins, a phylogenetic tree was constructed based on their sequences. In addition, models of the tertiary structures of cathepsin L were constructed using homology modeling methods and subjected to molecular dynamics simulations to obtain reasonable structure to understand its dynamical behavior. Our findings demonstrate that all of the potential Drosophila cathepsin L-like proteins contain at least one cathepsin propeptide inhibitor domain. Multiple sequence alignment and homology models clearly highlight the conservation of active site residues, disulfide bonds, and amino acid residues critical for inhibitor binding. Furthermore, comparative modeling indicates that the sequence/structure/function profiles and active site architectures are conserved.