Detection of multiscale pockets on protein surfaces using mathematical morphology

Detection of pockets on protein surfaces is an important step toward finding the binding sites of small molecules. In a previous study, we defined a pocket as a space into which a small spherical probe can enter, but a large probe cannot. The radius of the large probes corresponds to the shallowness of pockets. We showed that each type of binding molecule has a characteristic shallowness distribution. In this study, we introduced fundamental changes to our previous algorithm by using a 3D grid representation of proteins and probes, and the theory of mathematical morphology. We invented an efficient algorithm for calculating deep and shallow pockets (multiscale pockets) simultaneously, using several different sizes of spherical probes (multiscale probes). We implemented our algorithm as a new program, ghecom (grid‐based HECOMi finder). The statistics of calculated pockets for the structural dataset showed that our program had a higher performance of detecting binding pockets, than four other popular pocket‐finding programs proposed previously. The ghecom also calculates the shallowness of binding ligands, R_inaccess (minimum radius of inaccessible spherical probes) that can be obtained from the multiscale molecular volume. We showed that each part of the binding molecule had a bias toward a specific range of shallowness. These findings will be useful for predicting the types of molecules that will be most likely to bind putative binding pockets, as well as the configurations of binding molecules. The program ghecom is available through the Web server ( http://biunit.naist.jp/ghecom ). Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Protein-RNA interactions: exploring binding patterns with a three-dimensional superposition analysis of high resolution structures

MOTIVATION: The recognition of specific RNA sequences and structures by proteins is critical to our understanding of RNA processing, gene expression and viral replication. The diversity of RNA structures suggests that RNA recognition is substantially different than that of DNA. RESULTS: The atomic coordinates of 41 protein-RNA complexes have been used to probe composite nucleoside binding pockets that form the structural and chemical underpinnings of base recognition. Composite nucleoside binding pockets were constructed using three-dimensional superpositions of each RNA nucleoside. Unlike protein-DNA interactions which are dominated by accessibility, RNA recognition frequently occurs in non-canonical and single-strand-like structures that allow interactions to occur from a much wider set of geometries and make fuller use of unique base shapes and hydrogen-bonding ability. By constructing composites that include all van der Waals, hydrogen-bonding, stacking and general non-polar interactions made to a particular nucleoside, the strategies employed are made readily visible. Protein-RNA interactions can result in the formation of a glove-like tight binding pocket around RNA bases, but the size, shape and non-polar binding patterns differ between specific RNA bases. We show that adenine can be distinguished from guanine based on the size and shape of the binding pocket and steric exclusion of the guanine N2 exocyclic amino group. The unique shape and hydrogen-bonding pattern for each RNA base allow proteins to make specific interactions through a very small number of contacts, as few as two in some cases. AVAILABILITY: The program ENTANGLE is available from http://www.bioc.rice.edu/~shamoo     

Anatomy of protein pockets and cavities: measurement of binding site geometry and implications for ligand design.

Identification and size characterization of surface pockets and occluded cavities are initial steps in protein structure-based ligand design. A new program, CAST, for automatically locating and measuring protein pockets and cavities, is based on precise computational geometry methods, including alpha shape and discrete flow theory. CAST identifies and measures pockets and pocket mouth openings, as well as cavities. The program specifies the atoms lining pockets, pocket openings, and buried cavities; the volume and area of pockets and cavities; and the area and circumference of mouth openings. CAST analysis of over 100 proteins has been carried out; proteins examined include a set of 51 monomeric enzyme-ligand structures, several elastase-inhibitor complexes, the FK506 binding protein, 30 HIV-1 protease-inhibitor complexes, and a number of small and large protein inhibitors. Medium-sized globular proteins typically have 10-20 pockets/cavities. Most often, binding sites are pockets with 1-2 mouth openings; much less frequently they are cavities. Ligand binding pockets vary widely in size, most within the range 10(2)-10(3)A3. Statistical analysis reveals that the number of pockets and cavities is correlated with protein size, but there is no correlation between the size of the protein and the size of binding sites. Most frequently, the largest pocket/cavity is the active site, but there are a number of instructive exceptions. Ligand volume and binding site volume are somewhat correlated when binding site volume is < or =700 A3, but the ligand seldom occupies the entire site. Auxiliary pockets near the active site have been suggested as additional binding surface for designed ligands (Mattos C et al., 1994, Nat Struct Biol 1:55-58). Analysis of elastase-inhibitor complexes suggests that CAST can identify ancillary pockets suitable for recruitment in ligand design strategies. Analysis of the FK506 binding protein, and of compounds developed in SAR by NMR (Shuker SB et al., 1996, Science 274:1531-1534), indicates that CAST pocket computation may provide a priori identification of target proteins for linked-fragment design. CAST analysis of 30 HIV-1 protease-inhibitor complexes shows that the flexible active site pocket can vary over a range of 853-1,566 A3, and that there are two pockets near or adjoining the active site that may be recruited for ligand design.     

Modeling the noncovalent interactions at the metabolite binding site in purine riboswitches

We present gas phase quantum chemical studies on the metabolite binding interactions in two important purine riboswitches, the adenine and guanine riboswitches, at the B3LYP/6-31G(d,p) level of theory. In order to gain insights into the strucutral basis of their discriminative abilities of regulating gene expression, the structural properties and binding energies for the gas phase optimized geometries of the metabolite bound binding pocket are analyzed and compared with their respective crystal geometries. Kitaura-Morokuma analysis has been carried out to calculate and decompose the interaction energy into various components. NBO and AIM analysis has been carried out to understand the strength and nature of binding of the individual aptamer bases with their respective purine metabolites. The Y74 base, U in case of adenine riboswitch and C in case of guanine riboswitch constitutes the only differentiating element between the two binding pockets. As expected, with W:W cis G:C74 interaction contributing more than 50% of the total binding energy, the interaction energy for metabolite binding as calculated for guanine (-46.43 Kcal/mol) is nearly double compared to the corresponding value for that of adenine (-24.73 Kcal/mol) in the crystal context. Variations in the optimized geometries for different models and comparison of relative contribution to metabolite binding involving four conserved bases reveal the possible role of U47:U51 W:H trans pair in the conformational transition of the riboswitch from the metabolite free to metabolite bound state. Our results are also indicative of significant contributions from stacking and magnesium ion interactions toward cooperativity effects in metabolite recognition.     

NADP-Dependent enzymes. I: Conserved stereochemistry of cofactor binding

The ubiquitous redox cofactors nicotinamide adenine dinucleotides [NAD and NADP] are very similar molecules, despite their participation in substantially different biochemical processes. NADP differs from NAD in only the presence of an additional phosphate group esterified to the 2′-hydroxyl group of the ribose at the adenine end and yet NADP is confined with few exceptions to the reactions of reductive biosynthesis, whereas NAD is used almost exclusively in oxidative degradations. The discrimination between NAD and NADP is therefore an impressive example of the power of molecular recognition by proteins. The many known tertiary structures of NADP complexes affords the possibility for an analysis of their discrimination. A systematic analysis of several crystal structures of NAD(P)-protein complexes show that: 1) the NADP coenzymes are more flexible in conformation than those of NAD; 2) although the protein-cofactor interactions are largely conserved in the NAD complexes, they are quite variable in those of NADP; and 3) in both cases the pocket around the nicotinamide moiety is substrate dependent. The conserved and variable interactions between protein and cofactors in the respective binding pockets are reported in detail. Discrimination between NAD and NADP is essentially a consequence of the overall pocket and not of a few residues. A clear fingerprint in NAD complexes is a carboxylate side chain that chelates the diol group at the ribose near the adenine, whereas in NADP complexes an arginine side chain faces the adenine plane and interacts with the phosphomonoester. The latter type of interaction might be a general feature of recognition of nucleotides by proteins. Other features such as strand-like hydrogen bonding between the NADP diphosphate moeties and the protein are also significant. The NADP binding pocket properties should prove useful in protein engineering and design. © 1997 Wiley-Liss Inc.     

X-ray crystallographic analysis of the structural basis for the interaction of pokeweed antiviral protein with guanine residues of ribosomal RNA

Pokeweed antiviral protein (PAP) is a ribosome-inactivating protein (RIP), which enzymatically removes a single adenine base from a conserved, surface exposed loop sequence of ribosomal rRNA. We now present unprecedented experimental evidence that PAP can release not only adenine but guanine as well from Escherichia coli rRNA, albeit at a rate 20 times slower than for adenine. We also report X-ray structure analysis and supporting modeling studies for the interactions of PAP with guanine. Our modeling studies indicated that PAP can accommodate a guanine base in the active site pocket without large conformational changes. This prediction was experimentally confirmed, since a guanine base was visible in the active site pocket of the crystal structure of the PAP-guanine complex.     

Establishing the principles of recognition in the adenine-binding region of an aminoglycoside antibiotic kinase [APH(3')-IIIa

The protein-based molecular recognition of the adenine ring has implications throughout biological systems. In this paper, we discuss the adenine-binding region of an aminoglycoside antibiotic kinase [APH(3')-IIIa], which serves as an excellent model system for proteins that bind the adenine ring. This enzyme employs a hydrogen-bonding network involving water molecules along with enzyme backbone/side-chain atoms and a pi-pi stacking interaction to recognize the adenine ring. Our approach utilized site-directed mutagenesis, adenosine analogues and a variety of biophysical methods to probe the contacts in the adenine-binding region of APH(3')-IIIa. The results point to the polar nature of an adenine-Met90 contact in this binding pocket and the important role that Met90, the "gatekeeper" residue in structurally similar Ser/Thr protein kinases, plays in adenine binding. The results also suggest that small changes in the structure of the adenine ring can lead to significant changes in the ability of these analogues to occupy the adenine-binding region of the enzyme. Additional computational experiments indicate that both size and electronic factors are important in the binding of aromatic systems in this interaction-rich pocket. The principles governing adenine recognition established in this study may be applied to other protein-ligand complexes and used to navigate future studies directed at discovering potent and selective inhibitors of APH-type enzymes.     

Shape variation in protein binding pockets and their ligands

A common assumption about the shape of protein binding pockets is that they are related to the shape of the small ligand molecules that can bind there. But to what extent is that assumption true? Here we use a recently developed shape matching method to compare the shapes of protein binding pockets to the shapes of their ligands. We find that pockets binding the same ligand show greater variation in their shapes than can be accounted for by the conformational variability of the ligand. This suggests that geometrical complementarity in general is not sufficient to drive molecular recognition. Nevertheless, we show when considering only shape and size that a significant proportion of the recognition power of a binding pocket for its ligand resides in its shape. Additionally, we observe a "buffer zone" or a region of free space between the ligand and protein, which results in binding pockets being on average three times larger than the ligand that they bind.     

Structural Identification and Systematic Comparison of Phorbol Ester,Dioleoylglycerol, Alcohol and Sevoflurane Binding Sites in PKCδ C1A Domain

Protein kinase C (PKC) is a family of signal transducing enzymes that have been implicated in anesthetic preconditioning signaling cascade. Evidences are emerging that certain exogenous neuromodulators such as n-alkanols and general anesthetics can stimulate PKC activity by binding to regulatory C1A domain of the enzyme. However, the accurate binding sites in C1A domain as well as the molecular mechanism underlying binding-stimulated PKC activation still remain unelucidated. Here, we report a systematic investigation of the intermolecular interaction of human PKCδ C1A domain with its natural activator phorbol ester (PE) and co-activator dioleoylglycerol (DOG) as well as exogenous stimulators butanol, octanol and sevoflurane. The domain is computationally identified to potentially have three spatially vicinal ligand-binding pockets 1, 2 and 3, in which the pockets 1 and 2 have previously been determined as the binding sites of PE and DOG, respectively. Systematic cross-binding analysis reveals that long-chain octanol and DOG are well compatible with the flat, nonpolar pocket 2, where the nonspecific hydrophobic contacts and van der Waals packing are primarily responsible for the binding, while the general anesthetic sevoflurane prefer to interact with the rugged, polar pocket 3 through specific hydrogen bonds and electrostatic forces. Short-chain butanol appears to bind effectively none of the three pockets. In addition, the pocket 1 consists of two angled arms 1 and 2 that are also involved in pockets 2 and 3, respectively. Dynamics characterization imparts that binding of long-chain octanol and DOG to pocket 2 or binding of sevoflurane to pocket 3 can induce a conformational displacement in arm 1 or 2, thus further opening the included angle and enlarging pocket 1, which can improve the pocket 1-PE affinity via an allosteric mechanism, consequently stimulating the PE-induced PKCδ activation.     

Identification of protein functional surfaces by the concept of a split pocket

The function of a protein is often fulfilled via molecular interactions on its surfaces, so identifying the functional surface(s) of a protein is helpful for understanding its function. Here, we introduce the concept of a split pocket, which is a pocket that is split by a cognate ligand. We use a geometric approach that is site‐specific. Specifically, we first compute a set of all pockets in the protein with its ligand(s) and a set of all pockets with the ligand(s) removed and then compare the two sets of pockets to identify the split pocket(s) of the protein. To reduce the search space and expedite the process of surface partitioning, we design probe radii according to the physicochemical textures of molecules. Our method achieves a success rate of 96% on a benchmark test set. We conduct a large‐scale computation to identify ～19,000 split pockets from 11,328 structures (1.16 million potential pockets); for each pocket, we obtain residue composition, solvent‐accessible area, and molecular volume. With this database of split pockets, our method can be used to predict the functional surfaces of unbound structures. Indeed, the functional surface of an unbound protein may often be found from its similarity to remotely related bound forms that belong to distinct folds. Finally, we apply our method to identify glucose‐binding proteins, including unbound structures. Our study demonstrates the power of geometric and evolutionary matching for studying protein functional evolution and provides a framework for classifying protein functions by local spatial patterns of functional surfaces. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Nucleoside binding site of herpes simplex type 1 thymidine kinase analyzed by X-ray crystallography

The crystal structures of the full-length Herpes simplex virus type 1 thymidine kinase in its unligated form and in a complex with an adenine analogue have been determined at 1.9 A resolution. The unligated enzyme contains four water molecules in the thymidine pocket and reveals a small induced fit on substrate binding. The structure of the ligated enzyme shows for the first time a bound adenine analogue after numerous complexes with thymine and guanine analogues have been reported. The adenine analogue constitutes a new lead compound for enzyme-prodrug gene therapy. In addition, the structure of mutant Q125N modifying the binding site of the natural substrate thymidine in complex with this substrate has been established at 2.5 A resolution. It reveals that neither the binding mode of thymidine nor the polypeptide backbone conformation is altered, except that the two major hydrogen bonds to thymidine are replaced by a single water-mediated hydrogen bond, which improves the relative acceptance of the prodrugs aciclovir and ganciclovir compared with the natural substrate. Accordingly, the mutant structure represents a first step toward improving the virus-directed enzyme-prodrug gene therapy by enzyme engineering.     

Electrostatic potential of nucleotide-free protein is sufficient for discrimination between adenine and guanine-specific binding sites

Despite sharing many common features, adenine-binding and guanine-binding sites in proteins often show a clear preference for the cognate over the non-cognate ligand. We have analyzed electrostatic potential (ESP) patterns at adenine and guanine-binding sites of a large number of non-redundant proteins where each binding site was first annotated as adenine/guanine-specific or non-specific from a survey of primary literature. We show that more than 90% of ESP variance at the binding sites is accounted for by only two principal component ESP vectors, each aligned to molecular dipoles of adenine and guanine. Projected on these principal component vectors, the adenine/guanine-specific and non-specific binding sites, including adenine-containing dinucleotides, show non-overlapping distributions. Adenine or guanine specificities of the binding sites also show high correlation with the corresponding electrostatic replacement (cognate by non-cognate ligand) energies. High correlation coefficients (0.94 for 35 adenine-binding sites and 1.0 for 20 guanine-binding sites) were obtained when adenine/guanine specificities were predicted using the replacement energies. Our results demonstrate that ligand-free protein ESP is an excellent indicator for discrimination between adenine and guanine-specific binding sites and that ESP of ligand-free protein can be used as a tool to annotate known and putative purine-binding sites in proteins as adenine or guanine-specific.     

Purine sensing by riboswitches

Structured mRNA elements called riboswitches control gene expression by binding to small metabolites. Over a dozen riboswitch classes have been characterized that target a broad range of molecules and vary widely in size and secondary structure. Four of the known riboswitch classes recognize purines or modified purines. Three of these classes are closely related in conserved sequence and secondary structure, but members of these classes selectively recognize guanine, adenine or 2'-deoxyguanosine. Members of the fourth riboswitch class adopt a distinct structure to form a selective binding pocket for the guanine analogue preQ(1) (7-aminomethyl-7-deazaguanine). All four classes of purine-sensing riboswitches are most likely to recognize their respective metabolites by utilizing a riboswitch residue to make a canonical Watson-Crick base-pair with the ligand. This review will provide a summary of the purine-sensing riboswitches, as well as discuss the complex functions and applications of these RNAs.     

Purine sensing by riboswitches.

Structured mRNA elements called riboswitches control gene expression by binding to small metabolites. Over a dozen riboswitch classes have been characterized that target a broad range of molecules and vary widely in size and secondary structure. Four of the known riboswitch classes recognize purines or modified purines. Three of these classes are closely related in conserved sequence and secondary structure, but members of these classes selectively recognize guanine, adenine or 2'-deoxyguanosine. Members of the fourth riboswitch class adopt a distinct structure to form a selective binding pocket for the guanine analogue preQ(1) (7-aminomethyl-7-deazaguanine). All four classes of purine-sensing riboswitches are most likely to recognize their respective metabolites by utilizing a riboswitch residue to make a canonical Watson-Crick base-pair with the ligand. This review will provide a summary of the purine-sensing riboswitches, as well as discuss the complex functions and applications of these RNAs.     

Molecular basis for tetanus toxin coreceptor interactions

Tetanus toxin (TeNT) elicits spastic paralysis through the cleavage of vesicle-associated membrane protein-2 (VAMP-2) in neurons at the interneuronal junction of the central nervous system. While TeNT retrograde traffics from peripheral nerve endings to the interneuronal junction, there is limited understanding of the neuronal receptors utilized by tetanus toxin for the initial entry into nerve cells. Earlier studies implicated a coreceptor for tetanus toxin entry into neurons: a ganglioside binding pocket and a sialic acid binding pocket and that GT1b bound to each pocket. In this study, a solid phase assay characterized the ganglioside binding specificity and functional properties of both carbohydrate binding pockets of TeNT. The ganglioside binding pocket recognized the ganglioside sugar backbone, Gal-GalNAc, independent of sialic acid-(5) and sialic acid-(7) and GM1a was an optimal substrate for this pocket, while the sialic acid binding pocket recognized sialic acid-(5) and sialic acid-(7) with "b"series of gangliosides preferred relative to "a" series gangliosides. The high-affinity binding of gangliosides to TeNT HCR required functional ganglioside and sialic acid binding pockets, supporting synergistic binding to coreceptors. This analysis provides a model for how tetanus toxin utilizes coreceptors for high-affinity binding to neurons.     

Structures of Escherichia coli NAD synthetase with substrates and products reveal mechanistic rearrangements

Nicotinamide adenine dinucleotide synthetases (NADS) catalyze the amidation of nicotinic acid adenine dinucleotide (NAAD) to yield the enzyme cofactor nicotinamide adenine dinucleotide (NAD). Here we describe the crystal structures of the ammonia-dependent homodimeric NADS from Escherichia coli alone and in complex with natural substrates and with the reaction product NAD. The structures disclosed two NAAD/NAD binding sites at the dimer interface and an adenosine triphosphate (ATP) binding site within each subunit. Comparison with the Bacillus subtilis NADS showed pronounced chemical differences in the NAAD/NAD binding sites and less prominent differences in the ATP binding pockets. In addition, the E. coli NADS structures revealed unexpected dynamical rearrangements in the NAAD/NAD binding pocket upon NAAD-to-NAD conversion, which define a catalysis state and a substrate/product exchange state. The two states are adopted by concerted movement of the nicotinysyl moieties of NAAD and NAD, Phe-170, and residues 224-228, which may be triggered by differential coordination of a magnesium ion to NAAD and NAD. Phylogenetic structure comparisons suggest that the present results are relevant for designing species-specific antibiotics.     

Side-chain flexibility in proteins upon ligand binding

Ligand binding may involve a wide range of structural changes in the receptor protein, from hinge movement of entire domains to small side-chain rearrangements in the binding pocket residues. The analysis of side chain flexibility gives insights valuable to improve docking algorithms and can provide an index of amino-acid side-chain flexibility potentially useful in molecular biology and protein engineering studies. In this study we analyzed side-chain rearrangements upon ligand binding. We constructed two non-redundant databases (980 and 353 entries) of "paired" protein structures in complexed (holo-protein) and uncomplexed (apo-protein) forms from the PDB macromolecular structural database. The number and identity of binding pocket residues that undergo side-chain conformational changes were determined. We show that, in general, only a small number of residues in the pocket undergo such changes (e.g., approximately 85% of cases show changes in three residues or less). The flexibility scale has the following order: Lys > Arg, Gln, Met > Glu, Ile, Leu > Asn, Thr, Val, Tyr, Ser, His, Asp > Cys, Trp, Phe; thus, Lys side chains in binding pockets flex 25 times more often then do the Phe side chains. Normalizing for the number of flexible dihedral bonds in each amino acid attenuates the scale somewhat, however, the clear trend of large, polar amino acids being more flexible in the pocket than aromatic ones remains. We found no correlation between backbone movement of a residue upon ligand binding and the flexibility of its side chain. These results are relevant to 1. Reduction of search space in docking algorithms by inclusion of side-chain flexibility for a limited number of binding pocket residues; and 2. Utilization of the amino acid flexibility scale in protein engineering studies to alter the flexibility of binding pockets.     

Development of rationally designed DNA N6 adenine methyltransferase inhibitors

A series of bisubstrate inhibitors for DNA N6 adenine methyltransferase (Dam) have been synthesized by linking an amine analogue of S-adenosylmethionine to an aryl moiety designed to probe the binding pocket of the DNA adenine base. An initial structure-activity relationship study has identified substituents that increase inhibitor potency to the ～10 μM range and improve selectivity against the human cytosine methyltransferase Dnmt1.     

A Structural Basis for the Recognition of 2′-Deoxyguanosine by the Purine Riboswitch

Riboswitches are noncoding RNA elements that are commonly found in the 5′-untranslated region of bacterial mRNA. Binding of a small-molecule metabolite to the riboswitch aptamer domain guides the folding of the downstream sequence into one of two mutually exclusive secondary structures that directs gene expression. The purine riboswitch family, which regulates aspects of purine biosynthesis and transport, contains three distinct classes that specifically recognize guanine/hypoxanthine, adenine, or 2′-deoxyguanosine (dG). Structural analysis of the guanine and adenine classes revealed a binding pocket that almost completely buries the nucleobase within the core of the folded RNA. Thus, it is somewhat surprising that this family of RNA elements also recognizes dG. We have used a combination of structural and biochemical techniques to understand how the guanine riboswitch could be converted into a dG binder and the structural basis for dG recognition. These studies reveal that a limited number of sequence changes to a guanine-sensing RNA are required to cause a specificity switch from guanine to 2′-deoxyguanosine, and to impart an altered structure for accommodating the additional deoxyribose sugar moiety.     

The 2.1 A structure of a cysteine protease with proline specificity from ginger rhizome, Zingiber officinale.

A cysteine protease from ginger rhizome (GP-II) cleaves peptides and proteins with proline at the P(2) position. The unusual specificity for proline makes GP-II an attractive tool for protein sequencing and identification of stably folded domains in proteins. The enzyme is a 221 amino acid glycoprotein possessing two N-linked oligosaccharide chains (8% glycosylated by weight) at Asn99 and Asn156. The availability of the sequence of these glycosyl chains afforded the opportunity to observe their structure and impact on protein conformation. The three-dimensional structure of GP-II has been determined by X-ray crystallography to a resolution of 2.1 A (overall R-factor = 0.214, free R = 0.248). The overall structure of GP-II is similar to that of the homologous cysteine proteases papain, actinidin, and glycyl endopeptidase, folding into two distinct domains of roughly equal size which are divided by a cleft. The observed N-linked glycosyl chains (half the total carbohydrate sequence) participate in both crystallographic and noncrystallographic contacts, tethering the proteins together via hydrogen bonds to the carbohydrate residues without intervening ordered water molecules. The putative S(2) binding pocket (the proline recognition site) was identified by superposition of the GP-II structure with structures of four previously determined papain-inhibitor complexes. The particular enzymic amino acids forming the S(2) pocket of GP-II (Trp, Met, and Ala) are similar to those found in the proline binding pockets of the unrelated enzymes alpha-lytic protease and cyclophilin. However, there is no conserved three-dimensional arrangement of these residues between the three enzymes (i.e., no proline binding motif). Thus, the particular amino acids found at S(2) are consistent with a binding pocket for a moiety with the steric characteristics and charge distribution of proline. Size exclusion is also a mechanism for selectivity compared to the S(2) binding pocket of papain. The S(2) binding pocket of GP-II greatly restricts the size of the side chain which could be bound because of the occurrence of a tryptophan in place of the corresponding tyrosine in papain. In light of the nature of the binding pocket, the specificity of GP-II for proline over other small nonpolar amino acids may be attributed to a direct effect of proline on the substrate peptide backbone conformation.