Comparison of dipeptidyl carboxypeptidase and endopeptidase activities in the three enkephalin-hydrolysing metallopeptidases: "angiotensin-converting enzyme", thermolysin and "enkephalinase"

Angiotensin-converting enzyme (ACE), thermolysin and "enkephalinase", three metallopeptidases cleaving the Gly3-Phe4 amide bond of enkephalins, were compared regarding substrate specificity and effects of butanedione, an arginyl-directed reagent. The hydrolysis of enkephalins and analogues was more affected by the nature of P1 and P2 residues in the case of thermolysin than in those of ACE or "enkephalinase"; amidation of the C-terminal carboxylate decreased drastically the hydrolysis by ACE but only marginally by thermolysin and the effect was intermediate for "enkephalinase". With adequate model substrates, the ratio of dipeptidylcarboxypeptidase to tripeptidylcaroxypeptidase (endopeptidase) activities were of 25 for ACE, 3 for "enkephalinase" and only 0.3 for thermolysin. Finally a butanedione treatment increased thermolysin activity, but abolished ACE activity; it reduced "enkephalinase" activity by 80% when measured with a free C-terminal carboxylate enkephalin analogue but only slightly with the corresponding amidated derivative. A critical role of an Arg residue in ACE and, to a lesser extent, in "enkephalinase" (but not in thermolysin) is suggested to be responsible for the preferential dipeptidylcarboxypeptidase activity of these two enzymes.     

Analysis and predication of structural motifs in the glycolytic enzymes

Protein crystallography has determined the three-dimensional structures of 10 of the 13 enzymes of the glycolytic pathway. Diagrams and details of these enzyme structures are given in the paper. Most of the enzyme domains are variations and extensions of a many (4--9)-stranded, predominantly or totally parallel, beta-sheet that is shielded from solvent by alpha-helices (i.e. alpha/beta structures). There are strong structural similarities between the domains of some, but not all, of the enzymes. In particular the dinucleotide binding fold of lactate dehydrogenase and the beta-barrel of triose phosphate isomerase are found in other domains. General rules governing the topology and packing of alpha-helices against a beta-sheet provide a basis for the combinatorial prediction of the tertiary fold of glycolytic domains from their amino acid sequence and observed secondary structure. The predication algorithm demonstrates that there are severe restrictions on the number of possible structures. However, these restrictions do not fully explain some of the remarkable structural similarities between different enzymes that probably result from evolution from a common ancestor.     

Domain association in immunoglobulin molecules. The packing of variable domains

We have analyzed the structure of the interface between VL and VH domains in three immunoglobulin fragments: Fab KOL, Fab NEW and Fab MCPC 603. About 1800 A2 of protein surface is buried between the domains. Approximately three quarters of this interface is formed by the packing of the VL and VH beta-sheets in the conserved "framework" and one quarter from contacts between the hypervariable regions. The beta-sheets that form the interface have edge strands that are strongly twisted (coiled) by beta-bulges. As a result, the edge strands fold back over their own beta-sheet at two diagonally opposite corners. When the VL and VH domains pack together, residues from these edge strands form the central part of the interface and give what we call a three-layer packing; i.e. there is a third layer composed of side-chains inserted between the two backbone side-chain layers that are usually in contact. This three-layer packing is different from previously described beta-sheet packings. The 12 residues that form the central part of the three observed VL-VH packings are absolutely or very strongly conserved in all immunoglobulin sequences. This strongly suggests that the structure described here is a general model for the association of VL and VH domains and that the three-layer packing plays a central role in forming the antibody combining site.     

Revised 2.3 A structure of porcine pepsin: evidence for a flexible subdomain

A revised three-dimensional crystal structure of ethanol-inhibited porcine pepsin refined to an R-factor of 0.171 at 2.3 A resolution is presented and compared to the refined structures of the fungal aspartic proteinases: penicillopepsin, rhizopuspepsin, and endothiapepsin. Pepsin is composed of two nearly equal N and C domains related by an intra dyad. The overall polypeptide fold and active site structures are homologous for pepsin and the fungal enzymes. The weak inhibition of pepsin by ethanol can be explained by the presence of one or more ethanol molecules, in the vicinity of the active site carboxylates, which slightly alter the hydrogen-bonding network and which may compete with substrate binding in the active site. Structural superposition analysis showed that the N domains aligned better than the C-domains for pepsin and the fungal aspartic proteinases: 107-140 C alpha pairs aligned to 0.72-0.85 A rms for the N domains; 64-95 C alpha pairs aligned to 0.78-1.03 A rms for the C domains. The major structural difference between pepsin and the fungal enzymes concerns a newly described subdomain whose conformation varies markedly among these enzyme structures. The subdomain in pepsin comprises nearly 100 residues and is composed of two contiguous segments within the C domain (residues 192-212 and 223-299). the subdomain is connected, or "hinged," to a mixed beta-sheet that forms one of the structurally invariant, active site psi-loops. Relative subdomain displacements as large as a 21.0 degrees rotation and a 5.9 A translation were observed among the different enzymes. There is some suggestion in pepsin that the subdomain may be flexible and perhaps plays a structural role in mediating substrate binding, determining the substrate specificity, or in the activation of the zymogen.     

Interactions between an alpha-helix and a beta-sheet. Energetics of alpha/beta packing in proteins

Conformational energy computations have been carried out to determine the favorable ways of packing a right-handed alpha-helix on a right-twisted antiparallel or parallel beta-sheet. Co-ordinate transformations have been developed to relate the position and orientation of the alpha-helix to the beta-sheet. The packing was investigated for a CH3CO-(L-Ala)16-NHCH3 alpha-helix interacting with five-stranded beta-sheets composed of CH3CO-(L-Val)6-NHCH3 chains. All internal and external variables for both the alpha-helix and the beta-sheet were allowed to change during energy minimization. Four distinct classes of low-energy packing arrangements were found for the alpha-helix interacting with both the parallel and the anti-parallel beta-sheet. The classes differ in the orientation of the axis of the alpha-helix relative to the direction of the strands of the right-twisted beta-sheet. In the class with the most favorable arrangement, the alpha-helix is oriented along the strands of the beta-sheet, as a result of attractive non-bonded side-chain-side-chain interactions along the entire length of the alpha-helix. A class with nearly perpendicular orientation of the helix axis to the strands is also of low energy, because it allows similarly extensive attractive interactions. In the other two classes, the helix is oriented diagonally relative to the strands of the beta-sheet. In one of them, it interacts with the convex surface near the middle of the saddle-shaped twisted beta-sheet. In the other, it is oriented along the concave diagonal of the beta-sheet and, therefore, it interacts only with the corner regions of the sheet, so that this packing is energetically less favorable. The packing arrangements involving an antiparallel and a parallel beta-sheet are generally similar, although the antiparallel beta-sheet has been found to be more flexible. The major features of 163 observed alpha/beta packing arrangements in 37 proteins are accounted for in terms of the computed structural preferences. The energetically most favored packing arrangement is similar to the right-handed beta alpha beta crossover structure that is observed in proteins; thus, the preference for this connectivity arises in large measure from this energetically favorable interaction.     

Conformational dynamics of free and catalytically active thermolysin are indistinguishable by hydrogen/deuterium exchange mass spectrometry

Conformational dynamics are thought to be a prerequisite for the catalytic activity of enzymes. However, the exact relationship between structural fluctuations and function is not well understood. In this work hydrogen/deuterium exchange (HDX) and electrospray ionization mass spectrometry (ESI-MS) are used for exploring the conformational dynamics of thermolysin. Amide HDX reflects the internal mobility of proteins; regions that undergo frequent unfolding-refolding show faster exchange than segments that are highly stable. Thermolysin is a zinc protease with an active site that is located between two lobes. Substrate turnover is associated with hinge bending that leads to a closed conformation. Product release regenerates the open form, such that steady-state catalysis involves a continuous closing/opening cycle. HDX/ESI-MS with proteolytic peptide mapping in the absence of substrate shows that elements in the periphery of the two lobes are most mobile. A comparison with previous X-ray data suggests that these peripheral regions undergo quite pronounced structural changes during the catalytic cycle. In contrast, active site residues exhibit only a moderate degree of backbone flexibility, and the central zinc appears to be in a fairly rigid environment. The presence of both rigid and moderately flexible elements in the active site may reflect a carefully tuned balance that is required for function. Interestingly, the HDX behavior of catalytically active thermolysin is indistinguishable from that of the free enzyme. This result is consistent with the view that catalytically relevant motions preexist in the resting state and that enzyme function can only be performed within the limitations given by the intrinsic dynamics of the protein. The data presented in this work indicate the prevalence of stochastic elements in the function of thermolysin, rather than supporting a deterministic mechanism.     

Substrate and inhibitor studies of thermolysin-like neutral metalloendopeptidase from kidney membrane fractions. Comparison with bacterial thermolysin

The inhibitory constants of a series of synthetic N-carboxymethyl peptide inhibitors and the kinetic parameters (Km, kcat, and kcat/Km) of a series of model synthetic substrates were determined for the membrane-bound kidney metalloendopeptidase isolated from rabbit kidney and compared with those of bacterial thermolysin. The two enzymes show striking similarities with respect to structural requirements for substrate binding to the hydrophobic pocket at the S1' subsite of the active site. Both enzymes showed the highest reaction rates with substrates having leucine residues in this position while phenylalanine residues gave the lowest Km. The two enzymes were also inhibited by the same N-carboxymethyl peptide inhibitors. Although the mammalian enzyme was more susceptible to inhibition than its bacterial counterpart, structural variations in the inhibitor molecules affected the inhibitory constants for both enzymes in a similar manner. The two enzymes differed significantly, however, with respect to the effect of structural changes in the P1 and P2' positions of the substrate on the kinetic parameters of the reaction. The mammalian enzyme showed the highest reaction rates and specificity constants with substrates having the sequence -Phe-Gly-Phe- or -Phe-Ala-Phe- in positions P2, P1, and P1', respectively, while the sequence -Ala-Phe-Phe- was the most favored by the bacterial enzyme. The sequence -Gly-Gly-Phe- as found in enkephalins was not favored by either of the enzymes. Of the substrates having an aminobenzoate group in the P2' position, the mammalian enzyme favored those with the carboxyl group in the meta position while the bacterial enzyme favored those with the carboxyl group in the para position.(ABSTRACT TRUNCATED AT 250 WORDS)     

Subtle functional collective motions in pancreatic-like ribonucleases: from ribonuclease A to angiogenin

The analysis of the dynamic behavior of enzymes is fundamental to structural biology. A direct relationship between protein flexibility and biological function has been shown for bovine pancreatic ribonuclease (RNase A) (Rasmussen et al., Nature 1992;357:423-424). More recently, crystallographic studies have shown that functional motions in RNase A involve the enzyme beta-sheet regions that move concertedly on substrate binding and release (Vitagliano et al., Proteins 2002;46:97-104). These motions have been shown to correspond to intrinsic dynamic properties of the native enzyme by molecular dynamics (MD) simulations. To unveil the occurrence of these collective motions in other members of pancreatic-like superfamily, we carried out MD simulations on human angiogenin (Ang). Essential dynamics (ED) analyses performed on the trajectories reveal that Ang exhibits collective motions similar to RNase A, despite the limited sequence identity (33%) of the two proteins. Furthermore, we show that these collective motions are also present in ensembles of experimentally determined structures of both Ang and RNase A. Finally, these subtle concerted beta-sheet motions were also observed for other two members of the pancreatic-like superfamily by comparing the ligand-bound and ligand-free structures of these enzymes. Taken together, these findings suggest that pancreatic-like ribonucleases share an evolutionary conserved dynamic behavior consisting of subtle beta-sheet motions, which are essential for substrate binding and release.     

Interactions between two beta-sheets. Energetics of beta/beta packing in proteins

The analysis of the interactions between regularly folded segments of the polypeptide chain contributes to an understanding of the energetics of protein folding. Conformational energy-minimization calculations have been carried out to determine the favorable ways of packing two right-twisted beta-sheets. The packing of two five-stranded beta-sheets was investigated, with the strands having the composition CH3CO-(L-Ile)6-NHCH3 in one beta-sheet and CH3CO-(L-Val)6-NHCH3 in the other. Two distinct classes of low-energy packing arrangements were found. In the class with lowest energies, the strands of the two beta-sheets are aligned nearly parallel (or antiparallel) with each other, with a preference for a negative orientation angle, because this arrangement corresponds to the best complementary packing of the two twisted saddle-shaped beta-sheets. In the second class, with higher interaction energies, the strands of the two beta-sheets are oriented nearly perpendicular to each other. While the surfaces of the two beta-sheets are not complementary in this arrangement, there is good packing between the corner of one beta-sheet and the interior part of the surface of the other, resulting in a favorable energy of packing. Both classes correspond to frequently observed orientations of beta-sheets in proteins. In proteins, the second class of packing is usually observed when the two beta-sheets are covalently linked, i.e. when a polypeptide strand passes from one beta-sheet to the other, but we have shown here that a large contribution to the stabilization of this packing arrangement arises from noncovalent interactions.     

Effect of N-Methyl-4-Phenylpyridinium Ion on Monoamine Oxidase in a Clonal Rat Pheochromocytoma Cell Line, PC 12h

The effects of the neurotoxin N-methyl-4-phenylpyridinium ion (MPP+) on the enzymes involved in synthesis and catabolism of catecholamines were examined using a clonal rat pheochromocytoma cell line, PC12h, as a model of dopaminergic neurons. MPP+ added in the culture medium was found to be accumulated in PC12h cells after 30-min incubation. Monoamine oxidase (MAO) activity in PC12h cells was inhibited by MPP+ in a dose-dependent way from 10 nM to 10 microM, but concentrations of MPP+ higher than 100 microM were found to increase the MAO activity. At the lower concentrations MPP+ inhibited MAO noncompetitively with respect to the substrate, kynuramine, and at the higher concentrations it increased both the Km and the Vmax values of MAO toward the substrate. On the other hand, tyrosine hydroxylase activity and the dopamine concentrations in PC12 cells were not changed by incubation with MPP+ for 30 min, 60 min, or 24 h.     

Structure and mechanism of metallocarboxypeptidases

Metallocarboxpeptidases cleave C-terminal residues from peptide substrates and participate in a wide range of physiological processes, but they also contribute to human pathology. On the basis of structural information, we can distinguish between two groups of such metallopeptidases: cowrins and funnelins. Cowrins comprise protozoan, prokaryotic, and mammalian enzymes related to both neurolysin and angiotensin-converting enzyme and their catalytic domains contain 500-700 residues. They are ellipsoidal and traversed horizontally by a long, deep, narrow active-site cleft, in which the C-terminal residues are cut from oligopeptides and unstructured protein tails. The consensus cowrin structure contains a common core of 17 helices and a three-stranded beta-sheet, which participates in substrate binding. This protease family is characterized by a set of spatially conserved amino acids involved in catalysis, HEXXH+EXXS/G+H+Y/R+Y. Funnelins comprise structural relatives of the archetypal bovine carboxypeptidase A1 and feature mammalian, insect and bacterial proteins with strict carboxypeptidase activity. Their approximately 300-residue catalytic domains evince a consensus central eight-stranded beta-sheet flanked on either side by a total of eight helices. They also contain a characteristic set of conserved residues, HXXE+R+NR+H+Y+E, and their active-site clefts are rather shallow and lie at the bottom of a funnel-like cavity. Therefore, these enzymes act on a large variety of well-folded proteins. In both cowrins and funnelins, substrate hydrolysis follows a common general base/acid mechanism. A metal-bound solvent molecule ultimately performs the attack on the scissile peptide bond with the assistance of a strictly conserved glutamate residue.     

Crystal structure of botulinum neurotoxin type G light chain: serotype divergence in substrate recognition

The seven serotypes (A-G) of botulinum neurotoxins (BoNTs) block neurotransmitter release through their specific proteolysis of one of the three proteins of the soluble N-ethylmaleimide-sensitive-factor attachment protein receptor (SNARE) complex. BoNTs have stringent substrate specificities that are unique for metalloprotease in that they require exceptionally long substrates (1). To understand the molecular reasons for the unique specificities of the BoNTs, we determined the crystal structure of the catalytic light chain (LC) of Clostridium botulinum neurotoxin type G (BoNT/G-LC) at 2.35 A resolution. The structure of BoNT/G-LC reveals a C-terminal beta-sheet that is critical for LC oligomerization and is unlike that seen in the other LC structures. Its structural comparison with thermolysin and the available pool of LC structures reveals important serotype differences that are likely to be involved in substrate recognition of the P1' residue. In addition, structural and sequence analyses have identified a potential exosite of BoNT/G-LC that recognizes a SNARE recognition motif of VAMP.     

Atomic-resolution crystal structure of Borrelia burgdorferi outer surface protein A via surface engineering

Outer surface protein A (OspA) from Borrelia burgdorferi has an unusual dumbbell-shaped structure in which two globular domains are connected with a "single-layer" beta-sheet (SLB). The protein is highly soluble, and it has been recalcitrant to crystallization. Only OspA complexes with Fab fragments have been successfully crystallized. OspA contains a large number of Lys and Glu residues, and these "high entropy" residues may disfavor crystal packing because some of them would need to be immobilized in forming a crystal lattice. We rationally designed a total of 13 surface mutations in which Lys and Glu residues were replaced with Ala or Ser. We successfully crystallized the mutant OspA without a bound Fab fragment and extended structure analysis to a 1.15 Angstroms resolution. The new high-resolution structure revealed a unique backbone hydration pattern of the SLB segment in which water molecules fill the "weak spots" on both faces of the antiparallel beta-sheet. These well-defined water molecules provide additional structural links between adjacent beta-strands, and thus they may be important for maintaining the rigidity of the SLB that inherently lacks tight packing afforded by a hydrophobic core. The structure also revealed new information on the side-chain dynamics and on a solvent-accessible cavity in the core of the C-terminal globular domain. This work demonstrates the utility of extensive surface mutation in crystallizing recalcitrant proteins and dramatically improving the resolution of crystal structures, and provides new insights into the stabilization mechanism of OspA.     

Determination of the cleavage site of the presequence by mitochondrial processing peptidase on the substrate binding scaffold and the multiple subsites inside a molecular cavity

Mitochondrial processing peptidase (MPP) recognizes a large variety of basic presequences of mitochondrial preproteins and cleaves the single site, often including arginine, at the -2 position (P(2)). To elucidate the recognition and specific processing of the preproteins by MPP, we mutated to alanines at acidic residues conserved in a large internal cavity formed by the MPP subunits, alpha-MPP and beta-MPP, and analyzed the processing efficiencies for various preproteins. We report here that alanine mutations at a subsite in rat beta-MPP interacting with the P(2) arginine cause a shift in the processing site to the C-terminal side of the preprotein. Because of reduced interactions with the P(2) arginine, the mutated enzymes recognize not only the N-terminal authentic cleavage site with P(2) arginine but also the potential C-terminal cleavage site without a P(2) arginine. In fact, it competitively cleaves the two sites of the preprotein. Moreover, the acidified site of alpha-MPP, which binds to the distal basic site in the long presequence, recognized the authentic P(2) arginine as the distal site in compensation for ionic interaction at the proximal site in the mutant MPP. Thus, MPP seems to scan the presequence from beta- to alpha-MPP on the substrate binding scaffold inside the MPP cavity and finds the distal and P(2) arginines on the multiple subsites on both MPP subunits. A possible mechanism for substrate recognition and cleavage is discussed here based on the notable character of a subsite-deficient mutant of MPP in which the substrate specificity is altered.     

Need for Aromatic Residue at Position 115 for Proteolytic Activity Found by Site-directed Mutagenesis of Tryptophan 115 in Thermolysin

In thermolysin, tryptophan 115 seems to be at the S2 subsite. Trp-115 was replaced with tyrosine, phenylalanine, leucine, and valine during site-directed mutagenesis in order to evaluate the role of Trp-115 in the proteolytic activity of thermolysin. The mutant enzymes with Tyr-115 or Phe-115 had as much proteolytic activity as the wild-type enzyme, but the other two mutant enzymes had no activity. We found earlier that the substitution of Trp-115 with alanine, glutamic acid, lysine, and glutamine causes the enzyme to lose all activity, so an aromatic amino acid at position 115 seems to be essential for thermolysin.     

Structural features of neutral protease from Bacillus subtilis deduced from model-building and limited proteolysis experiments

The overall folding of neutral protease from Bacillus subtilis has been predicted by computer-aided modelling, taking as a basis the known three-dimensional structure of thermolysin. As expected from the 50% similarity of sequence between the two proteins, the structure of B. subtilis protease is similar to that of thermolysin, including the two-domain topology and location of elements of regular secondary structure (helices and strands), whereas specific differences were predicted in loop regions. A protruding and loose loop predicted in B. subtilis has been detected also experimentally by a limited proteolysis approach. Incubation of B. subtilis protease at pH 9.0 for 24 h at room temperature with trypsin at 20:1 ratio (by mass) leads to a specific and almost quantitative fission of the Arg214-Asn215 peptide bond located in a highly exposed, and thus probably flexible, loop of the protease. On the other hand, thermolysin was completely resistant to tryptic hydrolysis when reacted under identical conditions. The 'nicked' B. subtilis protease can be isolated by gel filtration chromatography at neutral pH, whereas the two constituting fragments 1-214 and 215-300 are separated under protein-denaturing conditions. Overall, these results indicate that the limited proteolysis approach can pinpoint a peculiar difference in surface structure between the two similar protein molecules of B. subtilis neutral protease and thermolysin and emphasize the potential use of proteolytic enzymes as structural probes of globular proteins.     

Comparative sequence and structure analysis reveal features of cold adaptation of an enzyme in the thermolysin family

Knowledge about the structural features underlying cold adaptation is important for designing enzymes of different industrial relevance. Vibriolysin from Antarctic bacterium strain 643 (VAB) is at present the only enzyme of the thermolysin family from an organism that thrive in extremely cold climate. In this study comparative sequence-structure analysis and molecular dynamics (MD) simulations were used to reveal the molecular features of cold adaptation of VAB. Amino acid sequence analysis of 44 thermolysin enzymes showed that VAB compared to the other enzymes has: (1) fewer arginines, (2) a lower Arg/(Lys + Arg) ratio, (3) a lower fraction of large aliphatic side chains, expressed by the (Ile + Leu)/(Ile + Leu + Val) ratio, (4) more methionines, (5) more serines, and (6) more of the thermolabile amino acid asparagine. A model of the catalytic domain of VAB was constructed based on homology with pseudolysin. MD simulations for 3 ns of VAB, pseudolysin, and thermolysin supported the assumption that cold-adapted enzymes have a more flexible three-dimensional (3D) structure than their thermophilic and mesophilic counterparts, especially in some loop regions. The structural analysis indicated that VAB has fewer intramolecular cation-pi electron interactions and fewer hydrogen bonds than its mesophilic (pseudolysin) and thermophilic (thermolysin) counterparts. Lysine is the dominating cationic amino acids involved in salt bridges in VAB, while arginine is dominating in thermolysin and pseudolysin. VAB has a greater volume of inaccessible cavities than pseudolysin and thermolysin. The electrostatic potentials on the surface of the catalytic domain were also more negative for VAB than for thermolysin and pseudolysin. Thus, the MD simulations, the structural patterns, and the amino acid composition of VAB relative to other enzymes of the thermolysin family suggest that VAB possesses the biophysical properties generally following adaptation to cold climate.     

A convenient automated method for the determination of proteolytic enzymes

An automated method for the assay of proteolytic enzymes is described. The insoluble dye-protein complex, Remazolbrilliant Blue-hide powder, is used as the substrate and is readily digested by trypsin, elastase, subtilisin, thermolysin, chymotrypsin, and to a lesser extent pronase, pepsin, and bacterial collagenase. The proteolytic activity of crude microbial culture preparations is expressed in terms of an equivalent concentration of crystalline trypsin, which itself can be readily determined within the concentration range 0.25–5.0 μg/ml.     

The crystal structure of pectate lyase Pel9A from Erwinia chrysanthemi

The "family 9 polysaccharide lyase" pectate lyase L (Pel9A) from Erwinia chrysanthemi comprises a 10-coil parallel beta-helix domain with distinct structural features including an asparagine ladder and aromatic stack at novel positions within the superhelical structure. Pel9A has a single high affinity calcium-binding site strikingly similar to the "primary" calcium-binding site described previously for the family Pel1A pectate lyases, and there is strong evidence for a common second calcium ion that binds between enzyme and substrate in the "Michaelis" complex. Although the primary calcium ion binds substrate in subsite -1, it is the second calcium ion, whose binding site is formed by the coming together of enzyme and substrate, that facilitates abstraction of the C5 proton from the sacharride in subsite +1. The role of the second calcium is to withdraw electrons from the C6 carboxylate of the substrate, thereby acidifying the C5 proton facilitating its abstraction and resulting in an E1cb-like anti-beta-elimination mechanism. The active site geometries and mechanism of Pel1A and Pel9A are closely similar, but the catalytic base is a lysine in the Pel9A enzymes as opposed to an arginine in the Pel1A enzymes.     

A novel method for the modelling of peptide ligands to their receptors

A knowledge-based approach to the modelling of enzyme-peptide inhibitor complexes is described. Given the structure of an enzyme, and knowledge of its binding site, the method seeks to predict the binding geometry of a peptide ligand. This novel method involves using examples of side-chain packing derived from proteins of known three-dimensional structure to define possible packing arrangements of a peptide inhibitor group to its binding site. A suite of programs, GEMINI, was written and used to predict the packing of pairs of amino acid groups from three inhibitors complexed to their enzymes for which the X-ray structures were available. These included the Phe group of the inhibitor H142 bound to endothiapepsin, the Leu group of CLT complexed to thermolysin and the C-terminus of Gly-L-Tyr bound to carboxypeptidase A. A detailed comparison of the modelled and observed inhibitor coordinates was made. This approach may be extended to modelling other types of protein interactions.