Structure of actinidin: details of the polypeptide chain conformation and active site from an electron density map at 2-8 A resolution

A 2·8 Å resolution electron density map of the sulphydryl protease, actinidin, has been calculated. Two isomorphous heavy-atom derivatives, prepared with uranyl acetate and dichloroethylenediamineplatinum(II), were used to calculate phases by the method of isomorphous replacement, giving an overall figure of merit of 0·81. The polypeptide chain is well-defined in the present map and many side-chains can be identified from their appearance. The molecule consists of a single chain of 220 residues, the last two of which appear disordered in the map and contains at least two, and probably three, disulphide bridges. The conformation of the polypeptide chain is remarkably similar to that of papain. It is folded into two domains, domain I consisting of residues 19–115 and 214–218, and domain II residues 1–18 and 116–213. There are three significant stretches of α-helix, involving residues 25–42, 69–81 and 120–129, together with several shorter pieces, while much of domain II consists of a twisted β-sheet structure. When compared with papain, actinidin has two additional residues inserted between 59 and 60, one inserted between 78 and 79, and four between 168 and 169 (papain numbering) while one residue (194) has been deleted from the papain structure. All these changes are in external parts of the molecule and have little effect on the conformation. The positions and orientations of the catalytically-important side-chains in the active site are virtually identical with those in papain, but some of the side-chains lining the non-polar binding pocket are clearly different.     

Crystal and molecular structure of the sulfhydryl protease calotropin DI at 3·2 Å resolution

The three-dimensional structure of the sulfhydryl protease calotropin DI from the madar plant, Calotropis gigantea, has been determined at 3·2 Å resolution using the multiple isomorphous replacement method with five heavy atom derivatives. A Fourier synthesis based on protein phases with a mean figure of merit of 0·857 was used for model building. The polypeptide backbone of calotropin DI is folded to form two distinct lobes, one of which is comprised mainly of α-helices, while the other is characterized by a system of all antiparallel pleated sheets. The overall molecular architecture closely resembles those found in the sulfhydryl proteases papain and actinidin.Despite the unknown amino acid sequence of calotropin DI a number of residues around its active center could be identified. These amino acid side-chains were found in a similar arrangement as the corresponding ones in papain and actinidin. The polypeptide chain between residues 1 and 18 of calotropin DI folds in a unique manner, providing a possible explanation for the unusual inability of calotropin DI to hydrolyze those synthetic substrates that papain and actinidin act upon.     

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.     

Cis peptide bonds in proteins: residues involved, their conformations, interactions and locations.

An analysis of a non-redundant set of protein structures from the Brookhaven Protein Data Bank has been carried out to find out the residue preference, local conformation, hydrogen bonding and other stabilizing interactions involving cis peptide bonds. This has led to a reclassification of turns mediated by cis peptides, and their average geometrical parameters have been evaluated. The interdependence of the side and main-chain torsion angles of proline rings provided an explanation why such rings in cis peptides are found to have the DOWN puckering. A comparison of cis peptides containing proline and non-proline residues show differences in conformation, location in the secondary structure and in relation to the centre of the molecule, and relative accessibilities of residues. Relevance of the results in mutation studies and the cis-trans isomerization during protein folding is discussed.     

The interrelationships of side-chain and main-chain conformations in proteins

The accurate determination of a large number of protein structures by X-ray crystallography makes it possible to conduct a reliable statistical analysis of the distribution of the main-chain and side-chain conformational angles, how these are dependent on residue type, adjacent residue in the sequence, secondary structure, residue-residue interactions and location at the polypeptide chain termini. The interrelationship between the main-chain (phi, psi) and side-chain (chi 1) torsion angles leads to a classification of amino acid residues that simplify the folding alphabet considerably and can be a guide to the design of new proteins or mutational studies. Analyses of residues occurring with disallowed main-chain conformation or with multiple conformations shed some light on why some residues are less favoured in thermophiles.     

Role of amino acid residues at turns in the conformational stability and folding of human lysozyme

To clarify the role of amino acid residues at turns in the conformational stability and folding of a globular protein, six mutant human lysozymes deleted or substituted at turn structures were investigated by calorimetry, GuHCl denaturation experiments, and X-ray crystal analysis. The thermodynamic properties of the mutant and wild-type human lysozymes were compared and discussed on the basis of their three-dimensional structures. For the deletion mutants, Delta47-48 and Delta101, the deleted residues are in turns on the surface and are absent in human alpha-lactalbumin, which is homologous to human lysozyme in amino acid sequence and tertiary structure. The stability of both mutants would be expected to increase due to a decrease in conformational entropy in the denatured state; however, both proteins were destabilized. The destabilizations were mainly caused by the disappearance of intramolecular hydrogen bonds. Each part deleted was recovered by the turn region like the alpha-lactalbumin structure, but there were differences in the main-chain conformation of the turn between each deletion mutant and alpha-lactalbumin even if the loop length was the same. For the point mutants, R50G, Q58G, H78G, and G37Q, the main-chain conformations of these substitution residues located in turns adopt a left-handed helical region in the wild-type structure. It is thought that the left-handed non-Gly residue has unfavorable conformational energy compared to the left-handed Gly residue. Q58G was stabilized, but the others had little effect on the stability. The structural analysis revealed that the turns could rearrange the main-chain conformation to accommodate the left-handed non-Gly residues. The present results indicate that turn structures are able to change their main-chain conformations, depending upon the side-chain features of amino acid residues on the turns. Furthermore, stopped-flow GuHCl denaturation experiments on the six mutants were performed. The effects of mutations on unfolding-refolding kinetics were significantly different among the mutant proteins. The deletion/substitutions in turns located in the alpha-domain of human lysozyme affected the refolding rate, indicating the contribution of turn structures to the folding of a globular protein.     

Thiol proteases and aldehyde dehydrogenases: Evolution from a common thiolesterase precursor?

The C-terminal 222 residues of human liver aldehyde dehydrogenase can be aligned with the C-terminal 226 residues of a thiol protease from Dictyostelium discoideum to yield 47 residue identities, including matching active site cysteine residues. A multiple alignment with three more aldehyde dehydrogenases and three more thiol proteases yields three regions with clustered residue similarities. In the tertiary structure of papain, these three regions are in close proximity although widely separated in primary structure, and many conserved residues are located in the active site groove. The three-dimensional relationships, the common thiol ester mechanisms of the enzymes, the locations of exon boundaries in the dehydrogenase and protease genes, and the conservation of internal salt-bridging and disulfide-paired residues in papain, all appear compatible with the hypothesis of an ancestral relationship between thiol proteases and aldehyde dehydrogenases.     

High-resolution study of the three-dimensional structure of horse heart metmyoglobin

The three-dimensional structure of horse heart metmyoglobin has been refined to a final R-factor of 15.5% for all observed data in the 6.0 to 1.9 A resolution range. The final model consists of 1242 non-hydrogen protein atoms, 154 water molecules and one sulfate ion. This structure has nearly ideal bonding and bond angle geometry. A Luzzati plot of the variation in R-factor with resolution yields an estimated mean co-ordinate error of 0.18 A. An extensive analysis of the pattern of hydrogen bonds formed in horse heart metmyoglobin has been completed. Over 80% of the polypeptide chain is involved in eight helical segments, of which seven are composed mainly of alpha-helical (3.6(13))-type hydrogen bonds; the remaining helix is composed entirely of 3(10) hydrogen bonds. Altogether, of 102 hydrogen bonds between main-chain atoms only six are not involved in helical structures, and four of these six occur within beta-turns. The majority of water molecules in horse heart metmyoglobin are found in solvent networks that range in size from two to 35 members. The size of water molecule networks can be rationalized on the basis of three factors: the number of hydrogen bonds to the protein surface, the presence of charged side-chain atoms, and the ability to bridge to neighboring molecules in the crystal lattice. Bridging water networks form the dominant intermolecular interactions. The backbone conformation of horse heart metmyoglobin is very similar to sperm whale metmyoglobin, with significant differences in secondary structure occurring only near residues 119 and 120, where residues 120 to 123 in sperm whale form a distorted type I reverse turn and the horse heart protein has a type II turn at residues 119 to 122. Nearly all of the hydrogen bonds between main-chain atoms (occurring mainly in helical regions) are common to both proteins, and more than half of the hydrogen bonds involving side-chain atoms observed in horse heart are also found in sperm whale metmyoglobin. Unlike sperm whale metmyoglobin, the heme iron atom in horse heart metmyoglobin is not significantly displaced from the plane of the heme group.     

The thiol proteinases from the latex of Carica papaya L. II. The primary structure of proteinase omega

The complete primary structure of the proteinase omega isolated from the latex of the Carica papaya fruits is given. The polypeptide chain contains 216 amino-acid residues, the alignment of which was deduced from sequence analyses of the native enzyme, the tryptic, chymotryptic, peptic and thermolysinolytic peptides and facilitated due to the considerable degree of homology with papain and actinidin. The location of the three disulfide bridges could be established with the help of peptic and thermolysinolytic fragments. Proteinase omega shares 148 identical amino-acid residues (68.5%) with papain and 108 ones (50%) with actinidin, including the three disulfide bridges and the free cysteine residue required for activity, as well as most of the other amino-acid residues involved in the catalytic mechanism and two thirds of the glycine residues which are of structural significance. The homology with other cysteine proteinases of different origin is discussed.     

Comparative resonance Raman spectroscopic and kinetic studies of acyl-enzymes involving papain, actinidin and papaya peptidase II.

Resonance Raman spectra are reported for a series of dithioacyl-enzymes involving actinidin (EC 3.4.22.14) and papaya peptidase II (the more basic monothiol cysteine proteinase of Carica papaya). The acyl groups are N-benzoylglycine and N-(beta-phenylpropionyl)glycine containing C = S or 13C = S at the ester function. Comparison of the data with those for the corresponding papain (EC 3.4.22.2) analogues [Storer, Lee & Carey (1983) Biochemistry 22, 4789-4796] allows us to define the conformation of the dithioacyl group in the catalytic site. In each case the dithioacyl group is bound in a single conformation known as conformer B, in which the glycinic nitrogen atom comes into close contact with the sulphur atom of the catalytic-site cysteine residue. For the N-(beta-phenylpropionyl)glycine dithioacyl-enzymes the torsional angles of the NH-CH2-C(= S) bonds assume values typical of an essentially relaxed non-strained state. However, for the N-benzoylglycine dithioacyl-enzymes there is evidence for a slightly perturbed conformer B, and the perturbation is most pronounced for N-benzoylglycine dithioacyl-actinidin. Values of k+2/Ks and k+3 for the reactions of papain, actinidin and papaya peptidase II with N-benzoylglycine and N-(beta-phenylpropionyl)glycine methyl thionoesters were obtained by a pre-steady-state kinetic study. Wide variation was found in k+2/Ks, but the values of k+3 are all similar. This general picture is supported by the results from a steady-state kinetic study of the reactions of the three enzymes with N-benzoyl-L-arginine-p-nitroanilide and with N-benzyloxycarbonyl-L-lysine p-nitrophenyl ester. The similarity of the values of k+3, together with the invariance of conformer B geometry at the P1 site, suggests that the chemistry of the deacylation process is highly conserved among these three cysteine proteinases.     

Conservation of cis prolyl bonds in proteins during evolution

In proteins and peptides, the vast majority of peptide bonds occurs in trans conformation, but a considerable fraction (about 5%) of X-Pro bonds adopts the cis conformation. Here we study the conservation of cis prolyl residues in evolutionary related proteins. We find that overall, in contrast to local, protein sequence similarity is a clear indicator for the conformation of prolyl residues. We observe that cis prolyl residues are more often conserved than trans prolyl residues, and both are more conserved than the surrounding amino acids, which show the same extent of conservation as the whole protein. The pattern of amino acid exchanges differs between cis and trans prolyl residues. Also, the cis prolyl bond is maintained in proteins with sequence identity as low as 20%. This finding emphasizes the importance of cis peptide bonds in protein structure and function.     

Structural basis of the unusual stability and substrate specificity of ervatamin C, a plant cysteine protease from Ervatamia coronaria

Ervatamin C is an unusually stable cysteine protease from the medicinal plant Ervatamia coronaria belonging to the papain family. Though it cleaves denatured natural proteins with high specific activity, its activity toward some small synthetic substrates is found to be insignificant. The three-dimensional structure and amino acid sequence of the protein have been determined from X-ray diffraction data at 1.9 A (R = 17.7% and R(free) = 19.0%). The overall structure of ervatamin C is similar to those of other homologous cysteine proteases of the family, folding into two distinct left and right domains separated by an active site cleft. However, substitution of a few amino acid residues, which are conserved in the other members of the family, has been observed in both the domains and also at the region of the interdomain cleft. Consequently, the number of intra- and interdomain hydrogen-bonding interactions is enhanced in the structure of ervatamin C. Moreover, a unique disulfide bond has been identified in the right domain of the structure, in addition to the three conserved disulfide bridges present in the papain family. All these factors contribute to an increase in the stability of ervatamin C. In this enzyme, the nature of the S2 subsite, which is the primary determinant of specificity of these proteases, is similar to that of papain, but at the S3 subsite, Ala67 replaces an aromatic residue, and has the effect of eliminating sufficient hydrophobic interactions required for S3-P3 stabilization. This provides the possible explanation for the lower activity of ervatamin C toward the small substrate/inhibitor. This substitution, however, does not affect the binding of denatured natural protein substrates to the enzyme significantly, as there exist a number of additional interactions at the enzyme-substrate interface outside the active site cleft.     

Proposed amino acid sequence and the 1.63 A X-ray crystal structure of a plant cysteine protease,ervatamin B: some insights into the structural basis of its stability and substrate specificity

The crystal structure of a cysteine protease ervatamin B, isolated from the medicinal plant Ervatamia coronaria, has been determined at 1.63 A. The unknown primary structure of the enzyme could also be traced from the high-quality electron density map. The final refined model, consisting of 215 amino acid residues, 208 water molecules, and a thiosulfate ligand molecule, has a crystallographic R-factor of 15.9% and a free R-factor of 18.2% for F > 2sigma(F). The protein belongs to the papain superfamily of cysteine proteases and has some unique properties compared to other members of the family. Though the overall fold of the structure, comprising two domains, is similar to the others, a few natural substitutions of conserved amino acid residues at the interdomain cleft of ervatamin B are expected to increase the stability of the protein. The substitution of a lysine residue by an arginine (residue 177) in this region of the protein may be important, because Lys --> Arg substitution is reported to increase the stability of proteins. Another substitution in this cleft region that helps to hold the domains together through hydrogen bonds is Ser36, replacing a conserved glycine residue in the others. There are also some substitutions in and around the active site cleft. Residues Tyr67, Pro68, Val157, and Ser205 in papain are replaced by Trp67, Met68, Gln156, and Leu208, respectively, in ervatamin B, which reduces the volume of the S2 subsite to almost one-fourth that of papain, and this in turn alters the substrate specificity of the enzyme.     

Sequence homologies, hydrophobic profiles and secondary structures of cathepsins B, H and L: comparison with papain and actinidin

The comparison of the amino acid sequences of 5 cysteine proteinases: papain, actinidin, rat cathepsins B and H and chicken cathepsin L, demonstrates a striking homology among their sequences. The N-terminal region (residues 1-70 in papain) and C-terminal region (residues 118-212 in papain) display the highest sequence homologies, whereas the lowest sequence homologies are observed in the middle region (residues 71-117 in papain); a segment where most insertions/deletions are observed. The highest sequence homology is observed between rat cathepsin H and chicken cathepsin L. As shown by X-ray studies, papain and actinidin have a clearly defined double domain structure. Each domain contains a core of non-polar side chains, which are retained in cathepsins B, H and L, except for the non-polar residue 203 of the core which is replaced by glutamic acid in cathepsin B. The percentage and the location of alpha-helix and beta-sheets of cathepsins B, H and L, assessed using the methods of Garnier et al. (1978, J. Mol. Biol. 120, 97-120) and Chou and Fasman (1974, Biochemistry 13, 222-245), show that the main ordered structures in papain and actinidin are probably retained in cathepsins B, H and L. The differences observed occur essentially in the middle region, a place where sequences display the lowest homologies and which is far removed from the active site.     

NMR studies of structure, hydrogen exchange, and main-chain dynamics in a disrupted-core mutant of thioredoxin.

Core-packing mutants of proteins often approach molten globule states, and hence may have attributes of folding intermediates. We have studied a core-packing mutant of thioredoxin, L78K, in which a leucine residue is substituted by lysine, using 15N heteronuclear two- and three-dimensional NMR. Chemical shift differences between the mutant and wild-type main-chain resonances reveal that structural changes caused by the mutation are localized within 12 A of the altered side chain. The majority of resonances are unchanged, as are many 1H-1H NOEs indicative of the main-chain fold, suggesting that the structure of L78K is largely similar to wild type. Hydrogen exchange studies reveal that residues comprising the central beta-sheet of both mutant and wild-type proteins constitute a local unfolding unit, but with the unfolding/folding equilibrium approximately 12 times larger in L78K. The dynamics of main-chain NH bonds in L78K were studied by 15N spin relaxation and compared with a previous study of wild type. Order parameters for angular motion of NH bonds in the mutant are on average lower than in wild type, suggesting greater spatial freedom on a rapid time scale, but may also be related to different rotational correlation times in the two proteins. There is also evidence of greater conformational exchange in the mutant. Differences between mutant and wild type in hydrogen exchange and main-chain dynamics are not confined to the vicinity of the mutation. We infer that mispacking of the protein core in one location affects local dynamics and stability throughout.     

The refined 2.15 A X-ray crystal structure of human liver cathepsin B: the structural basis for its specificity.

From the lysosomal cysteine proteinase cathepsin B, isolated from human liver in its two-chain form, monoclinic crystals were obtained which contain two molecules per asymmetric unit. The molecular structure was solved by a combination of Patterson search and heavy atom replacement methods (simultaneously with rat cathepsin B) and refined to a crystallographic R value of 0.164 using X-ray data to 2.15 A resolution. The overall folding pattern of cathepsin B and the arrangement of the active site residues are similar to the related cysteine proteinases papain, actinidin and calotropin DI. 166 alpha-carbon atoms out of 248 defined cathepsin B residues are topologically equivalent (with an r.m.s. deviation of 1.04 A) with alpha-carbon atoms of papain. However, several large insertion loops are accommodated on the molecular surface and modify its properties. The disulphide connectivities recently determined for bovine cathepsin B by chemical means were shown to be correct. Some of the primed subsites are occluded by a novel insertion loop, which seems to favour binding of peptide substrates with two residues carboxy-terminal to the scissile peptide bond; two histidine residues (His110 and His111) in this "occluding loop' provide positively charged anchors for the C-terminal carboxylate group of such polypeptide substrates. These structural features explain the well-known dipeptidyl carboxypeptidase activity of cathepsin B. The other subsites adjacent to the reactive site Cys29 are relatively similar to papain; Glu245 in the S2 subsite favours basic P2-side chains. The above mentioned histidine residues, but also the buried Glu171 might represent the group with a pKa of approximately 5.5 near the active site, which governs endo- and exopeptidase activity. The "occluding loop' does not allow cystatin-like protein inhibitors to bind to cathepsin B as they do to papain, consistent with the reduced affinity of these protein inhibitors for cathepsin B compared with the related plant enzymes.     

The structure of actinidin at 5-5 A resolution.

The structure of actinidin, a sulphydryl protease obtained from the fruit of Actinidia chinensis, has been determined from X-ray crystallographic data to a resolution of 5.5 Å. Three isomorphous heavy atom derivatives, prepared with uranyl acetate, dichloroethylenediamine platinum(II) and potassium iodomercurate(II), were used in the phase calculation, giving a mean figure of merit of 0.88. The molecule can be described as an oblate ellipsoid with approximate dimensions 50 Å × 40 Å × 36 Å, and consists of two globular units separated by a shallow cleft. Binding studies with mercuric chloride reveal two sites of attachment, both within this cleft, and although both sites are of low occupancy it is probable that one or other marks the position of the active sulphydryl group. Although the folding of the polypeptide chain of actinidin cannot be followed with certainty, it appears to be closely similar to that of papain, suggesting that these are members of a family of homologous proteins.     

Homology modelling and protein engineering strategy of subtilases, the family of subtilisin-like serine proteinases.

Subtilases are members of the family of subtilisin-like serine proteases. Presently, greater than 50 subtilases are known, greater than 40 of which with their complete amino acid sequences. We have compared these sequences and the available three-dimensional structures (subtilisin BPN', subtilisin Carlsberg, thermitase and proteinase K). The mature enzymes contain up to 1775 residues, with N-terminal catalytic domains ranging from 268 to 511 residues, and signal and/or activation-peptides ranging from 27 to 280 residues. Several members contain C-terminal extensions, relative to the subtilisins, which display additional properties such as sequence repeats, processing sites and membrane anchor segments. Multiple sequence alignment of the N-terminal catalytic domains allows the definition of two main classes of subtilases. A structurally conserved framework of 191 core residues has been defined from a comparison of the four known three-dimensional structures. Eighteen of these core residues are highly conserved, nine of which are glycines. While the alpha-helix and beta-sheet secondary structure elements show considerable sequence homology, this is less so for peptide loops that connect the core secondary structure elements. These loops can vary in length by greater than 150 residues. While the core three-dimensional structure is conserved, insertions and deletions are preferentially confined to surface loops. From the known three-dimensional structures various predictions are made for the other subtilases concerning essential conserved residues, allowable amino acid substitutions, disulphide bonds, Ca(2+)-binding sites, substrate-binding site residues, ionic and aromatic interactions, proteolytically susceptible surface loops, etc. These predictions form a basis for protein engineering of members of the subtilase family, for which no three-dimensional structure is known.     

Modelling family 2 cystatins and their interaction with papain

Cystatins are extensively studied cysteine protease inhibitors, found in wide range of organisms with highly conserved structural folds. S-type of cystatins is well known for their abundance in saliva, high selectivity and poorer activity towards host cysteine proteases in comparison to their immediate ancestor cystatin C. Despite more than 90% sequence similarity, the members of this group show highly dissimilar binding affinity towards papain. Cystatin M/E is a potent inhibitor of legumain and papain like cysteine proteases and recognized for its involvement in skin barrier formation and potential role as a tumor suppressor gene. However, the structures of these proteins and their complexes with papain or legumain are still unknown. In the present study, we have employed computational methods to get insight into the interactions between papain and cystatins. Three-dimensional structures of the cystatins are generated by homology modelling, refined with molecular dynamics simulation, validated through numerous web servers and finally complexed with papain using ZDOCK algorithm in Discovery Studio. A high degree of shape complementarity is observed within the complexes, stabilized by numerous hydrogen bonds (HB) and hydrophobic interactions. Using interaction energy, HB and solvent accessible surface area analyses, we have identified a series of key residues that may be involved in papain–cystatin interaction. Differential approaches of cystatins towards papain are also noticed which are possibly responsible for diverse inhibitory activity within the group. These findings will improve our understanding of fundamental inhibitory mechanisms of cystatin and provide clues for further research.     

Identification of conserved residue patterns in small beta-barrel proteins

Our abilities to predict three-dimensional conformation of a polypeptide, given its amino acid sequence, remain limited despite advances in structure analysis. Analysis of structures and sequences of protein families with similar secondary structural elements, but varying topologies, might help in addressing this problem. We have studied the small beta-barrel class of proteins characterized by four strands (n = 4) and a shear number of 8 (S = 8) to understand the principles of barrel formation. Multiple alignments of the various protein sequences were generated for the analysis. Positional entropy, as a measure of residue conservation, indicated conservation of non-polar residues at the core positions. The presence of a type II beta-turn among the various barrel proteins considered was another strikingly invariant feature. A conserved glycyl-aspartyl dipeptide at the beta-turn appeared to be important in guiding the protein sequence into the barrel fold. Molecular dynamics simulations of the type II beta-turn peptide suggested that aspartate is a key residue in the folding of the protein sequence into the barrel. Our study suggests that the conserved type II beta-turn and the non-polar residues in the barrel core are crucial for the folding of the protein's primary sequence into the beta-barrel conformation.