Substrate Binding Subsites Research Articles

Enzymatic Properties of Rhea Lysozyme

Rhea lysozyme was analyzed for its enzymatic properties both lytic and oligomer activities to reveal the structural and functional relationships of goose type lysozyme. Rhea lysozyme had the highest lytic activity at pH 6, followed by ostrich and goose at pH 5.5-6, whereas the optimum of cassowary was at pH 5. pH profile was correlated to the net charge of each molecule surface. On the other hand, the pH optimum for oligomer substrate was found to be pH 4, indicating the mechanism of rhea catalysis as a general acid. The time-course of the reaction was studied using beta-1,4-linked oligosaccharide of N-acetylglucosamine (GlcNAc) with a polymerization degree of n ((GlcNAc)n) (n=4, 5, and 6) as the substrate. This enzyme hydrolyzed (GlcNAc)6 in an endo-splitting manner, which produced (GlcNAc)3+(GlcNAc)3 predominating over that to (GlcNAc)2+ (GlcNAc)4. This indicates that the lysozyme hydrolyzed preferentially the third glycosidic linkage from the nonreducing end. Theoretical analysis has shown the highest rate constant value at 1.5 s-1 with (GlcNAc)6. This confirmed six substrate binding subsites as goose lysozyme (Honda, Y., and Fukamizo, T., Biochim. Biophys. Acta, 1388, 53-65 (1998)). The different binding free energy values for subsites B, C, F, and G from goose lysozyme might responsible for the amino acid substitutions, Asn122Ser and Phe123Met, located at the subsite B.

Comparison of the substrate specificity of two potyvirus proteases

The substrate specificity of the nuclear inclusion protein a (NIa) proteolytic enzymes from two potyviruses, the tobacco etch virus (TEV) and tobacco vein mottling virus (TVMV), was compared using oligopeptide substrates. Mutations were introduced into TEV protease in an effort to identify key determinants of substrate specificity. The specificity of the mutant enzymes was assessed by using peptides with complementary substitutions. The crystal structure of TEV protease and a homology model of TVMV protease were used to interpret the kinetic data. A comparison of the two structures and the experimental data suggested that the differences in the specificity of the two enzymes may be mainly due to the variation in their S4 and S3 binding subsites. Two key residues predicted to be important for these differences were replaced in TEV protease with the corresponding residues of TVMV protease. Kinetic analyses of the mutants confirmed that these residues play a role in the specificity of the two enzymes. Additional residues in the substrate-binding subsites of TEV protease were also mutated in an effort to alter the specificity of the enzyme.

Homology models of main proteinase from coronavirus associated with SARS

In this study, two homology models of the main proteinase (M pro) from the novel coronavirus associated with severe acute respiratory syndrome (SARS-CoV) were constructed. These models reveal three distinct functional domains, in which an intervening loop connecting domains II and III as well as a catalytic cleft containing the substrate binding subsites S1 and S2 between domains I and II are observed. S2 exhibits structural variations more significantly than S1 during the 200 ps molecular dynamics simulations because it is located at the open mouth of the catalytic cleft and the amino acid residues lining up this subsite are least conserved. In addition, the higher structural variation of S2 makes it flexible enough to accommodate a bulky hydrophobic residue from the substrate.

Homology Models and Molecular Dynamics Simulations of Main Proteinase from Coronavirus Associated with Severe Acute Respiratory Syndrome (SARS).

In this study, two structural models (denoted as MproST and MproSH) of the main proteinase (Mpro) from the novel coronavirus associated with severe acute respiratory syndrome (SARS‐CoV) were constructed based on the crystallographic structures of Mpro from transmissible gastroenteritis coronavirus (TGEV) (MproT) and human coronavirus HcoV‐229E (MproH), respectively. Various 200 ps molecular dynamics simulations were subsequently performed to investigate the dynamics behaviors of several structural features. Both MproST and MproSH exhibit similar folds as their respective template proteins. These structural models reveal three distinct functional domains as well as an intervening loop connecting domains II and III as found in both template proteins. In addition, domain III of these structures exhibits the least secondary structural conservation. A catalytic cleft containing the substrate binding subsites S1 and the S2 between domains I and II are also observed in these structural models. Although these structures share many common features, the most significant difference occurs at the S2 subsite, where the amino acid residues lining up this subsite are least conserved. It may be a critical challenge for designing anti‐SARS drugs by simply screening the known database of proteinase inhibitors.

Characterization of trans- and cis-cleavage activity of the SARS coronavirus 3CL pro protease: basis for the in vitro screening of anti-SARS drugs

Severe acute respiratory syndrome (SARS) has been globally reported. A novel coronavirus (CoV), SARS-CoV, was identified as the etiological agent of the disease. SARS-CoV 3C-like protease (3CL pro) mediates the proteolytic processing of replicase polypeptides 1a and 1ab into functional proteins, playing an important role in viral replication. In this study, we demonstrated the expression of the SARS-CoV 3CL pro in Escherichia coli and Vero cells, and then characterized the in vitro trans-cleavage and the cell-based cis-cleavage by the 3CL pro. Mutational analysis of the 3CL pro demonstrated the importance of His41, Cys145, and Glu166 in the substrate-binding subsite S1 for keeping the proteolytic activity. In addition, alanine substitution of the cleavage substrates indicated that Gln- P1 in the substrates mainly determined the cleavage efficiency. Therefore, this study not only established the quantifiable and reliable assay for the in vitro and cell-based measurement of the 3CL pro activity, but also characterized the molecular interaction of the SARS-CoV 3CL pro with the substrates. The results will be useful for the rational development of the anti-SARS drugs.

Tandem repeat of a seven-bladed beta-propeller domain in oligoxyloglucan reducing-end-specific cellobiohydrolase.

Oligoxyloglucan reducing-end-specific cellobiohydrolase (OXG-RCBH; EC 3.2.1.150) is an exoglucanase that recognizes the reducing end of oligoxyloglucan and releases two glucosyl residue segments from the main chain. The X-ray crystal structure of OXG-RCBH determined at 2.2 A resolution reveals a unique feature of this enzyme; OXG-RCBH consists of a tandem repeat of two similar domains, which are both folded into seven-bladed beta-propeller structures. The sequence alignment of the propeller blades, based on the structure, indicates that a weak repeat of the amino acid sequence occurred seven times to construct each domain. There is a cleft that can accommodate the substrate oligosaccharide between the two domains, which is a putative substrate binding subsite. Mutation of either Asp35 or Asp465, located in the putative catalytic center, to Asn resulted in a protein with no detectable catalytic activity, indicating the critical role of these amino acids in catalysis.

The structure of endo-beta-1,4-galactanase from Bacillus licheniformis in complex with two oligosaccharide products.

The beta-1,4-galactanase from Bacillus licheniformis (BLGAL) is a plant cell-wall-degrading enzyme involved in the hydrolysis of beta-1,4-galactan in the hairy regions of pectin. The crystal structure of BLGAL was determined by molecular replacement both alone and in complex with the products galactobiose and galactotriose, catching a first crystallographic glimpse of fragments of beta-1,4-galactan. As expected for an enzyme belonging to GH-53, the BLGAL structure reveals a (betaalpha)(8)-barrel architecture. However, BLGAL betaalpha-loops 2, 7 and 8 are long in contrast to the corresponding loops in structures of fungal galactanases determined previously. The structure of BLGAL additionally shows a calcium ion linking the long betaalpha-loops 7 and 8, which replaces a disulphide bridge in the fungal galactanases. Compared to the substrate-binding subsites predicted for Aspergillus aculeatus galactanase (AAGAL), two additional subsites for substrate binding are found in BLGAL, -3 and -4. A comparison of the pattern of galactan and galactooligosaccharides degradation by AAGAL and BLGAL shows that, although both are most active on substrates with a high degree of polymerization, AAGAL can degrade galactotriose and galactotetraose efficiently, whereas BLGAL prefers longer oligosaccharides and cannot hydrolyze galactotriose to any appreciable extent. This difference in substrate preference can be explained structurally by the presence of the extra subsites -3 and -4 in BLGAL.

Gene amplification and cold adaptation of pepsin in Antarctic fish. A possible strategy for food digestion at low temperature

Cold-adapted organisms have developed a number of adjustments at the molecular level to maintain metabolic functions at low temperatures. Among other features, they can produce enzymes characterized by a high turnover number or a high catalytic efficiency. The present work is aimed at investigating the process of food digestion at low temperature through the study of pepsins in Antarctic notothenioids. For such a purpose, we have cloned and sequenced three forms of pepsin A and a single form of gastricsin from the gastric mucosa of Trematomus bernacchii (rock cod). Phylogenetic analysis has suggested that the three pepsin A isotypes arose from two gene duplication events leading to the most ancestral pepsin A3 and to the most recent forms represented by pepsin A1 and pepsin A2. Molecular modeling has unraveled significant structural differences in these enzymes with respect to their mesophilic counterparts. Hydropathy and flexibility determined on the substrate-binding subsites of Antarctic and mesophilic pepsins have shown for pepsin A2 reduced hydropathy and increased flexibility at the level of the substrate cleft, features typical of cold-adapted enzymes. Northern blot analysis of RNA from rock cod gastric mucosa hybridized with molecular probes designed on specific regions of different pepsin forms has shown that rock cod pepsin genes are expressed at comparable levels. The present results suggest that the Antarctic rock cod adopted two different strategies to accomplish efficient protein digestion at low temperature. One mechanism is the gene duplication that increases enzyme production to compensate for the reduced kinetic efficiency, the other is the expression of a new enzyme provided with features typical of cold-adapted enzymes.

Evolutionary relationships of the prolyl oligopeptidase family enzymes

The prolyl oligopeptidase (POP) family of serine proteases includes prolyl oligopeptidase, dipeptidyl peptidase IV, acylaminoacyl peptidase and oligopeptidase B. The enzymes of this family specifically hydrolyze oligopeptides with less than 30 amino acids. Many of the POP family enzymes have evoked pharmaceutical interest as they have roles in the regulation of peptide hormones and are involved in a variety of diseases such as dementia, trypanosomiasis and type 2 diabetes. In this study we have clarified the evolutionary relationships of these four POP family enzymes and analyzed POP sequences from different sources. The phylogenetic trees indicate that the four enzymes were present in the last common ancestor of all life forms and that the beta-propeller domain has been part of the family for billions of years. There are striking differences in the mutation rates between the enzymes and POP was found to be the most conserved enzyme of this family. However, the localization of this enzyme has changed throughout evolution, as three archaeal POPs seem to be membrane bound and one third of the bacterial as well as two eukaryotic POPs were found to be secreted out of the cell. There are also considerable distinctions between the mutation rates of the different substrate binding subsites of POP. This information may help in the development of species-specific POP inhibitors.

ACE2 X-Ray Structures Reveal a Large Hinge-bending Motion Important for Inhibitor Binding and Catalysis

The angiotensin-converting enzyme (ACE)-related carboxypeptidase, ACE2, is a type I integral membrane protein of 805 amino acids that contains one HEXXH + E zinc-binding consensus sequence. ACE2 has been implicated in the regulation of heart function and also as a functional receptor for the coronavirus that causes the severe acute respiratory syndrome (SARS). To gain further insights into this enzyme, the first crystal structures of the native and inhibitor-bound forms of the ACE2 extracellular domains were solved to 2.2- and 3.0-Å resolution, respectively. Comparison of these structures revealed a large inhibitor-dependent hinge-bending movement of one catalytic subdomain relative to the other (∼16°) that brings important residues into position for catalysis. The potent inhibitor MLN-4760 ((S,S)-2-{1-carboxy-2-[3-(3,5-dichlorobenzyl)-3H-imidazol4-yl]-ethylamino}-4-methylpentanoic acid) makes key binding interactions within the active site and offers insights regarding the action of residues involved in catalysis and substrate specificity. A few active site residue substitutions in ACE2 relative to ACE appear to eliminate the S2′ substrate-binding subsite and account for the observed reactivity change from the peptidyl dipeptidase activity of ACE to the carboxypeptidase activity of ACE2.

Tyrosine 105 and Threonine 212 at Outermost Substrate Binding Subsites –6 and +4 Control Substrate Specificity, Oligosaccharide Cleavage Patterns, and Multiple Binding Modes of Barley α-Amylase 1

The role in activity of outer regions in the substrate binding cleft in alpha-amylases is illustrated by mutational analysis of Tyr(105) and Thr(212) localized at subsites -6 and +4 (substrate cleavage occurs between subsites -1 and +1) in barley alpha-amylase 1 (AMY1). Tyr(105) is conserved in plant alpha-amylases whereas Thr(212) varies in these and related enzymes. Compared with wild-type AMY1, the subsite -6 mutant Y105A has 140, 15, and <1% activity (k(cat)/K(m)) on starch, amylose DP17, and 2-chloro-4-nitrophenyl beta-d-maltoheptaoside, whereas T212Y at subsite +4 has 32, 370, and 90% activity, respectively. Thus engineering of aromatic stacking interactions at the ends of the 10-subsite long binding cleft affects activity very differently, dependent on the substrate. Y105A dominates in dual subsite -6/+4 [Y105A/T212(Y/W)]AMY1 mutants having almost retained and low activity on starch and oligosaccharides, respectively. Bond cleavage analysis of oligosaccharide degradation by wild-type and mutant AMY1 supports that Tyr(105) is critical for binding at subsite -6. Substrate binding is improved by T212(Y/W) introduced at subsite +4 and the [Y105A/T212(Y/W)]AMY1 double mutants synergistically enhanced productive binding of the substrate aglycone. The enzymatic properties of the series of AMY1 mutants suggest that longer substrates adopt several binding modes. This is in excellent agreement with computed distinct multiple docking solutions observed for maltododecaose at outer binding areas of AMY1 beyond subsites -3 and +3.

Coevolution of insect trypsins and inhibitors.

Many plant proteinase inhibitors have lysine at the P1 position, presumably to avoid hydrolysis by insect trypsins. Lepidopteran trypsins appear to have adapted to resist proteinase inhibitors through increased inhibitor hydrolysis and decreased binding to inhibitor hydrophilic reactive sites. Lepidopteran digestive trypsins prefer lysine at the P1 position and have substrate binding subsites more hydrophobic than trypsins from insects in other orders. All available sequences of sensitive and inhibitor-insensitive insect trypsins were aligned with porcine trypsin, for which interactions with Kunitz and Bowman-Birk inhibitor are known. After discounting conserved positions and positions not typical of sensitive or insensitive trypsins, the following residues were considered important to insect trypsin-PI interactions (chymotrypsin numbering): 60, 94, 97, 98, 99, 188, 190, 213, 215, 217, 219, 228. These residues support the Neighbor Joining analysis tree branches separating sensitive and insensitive trypsin sequences. Primary sequences interacting with PIs are around the active site, with some forming part of the S1 (188, 217, 219 and 228) or S4 (99, 215) pockets.

Papain-like lysosomal cysteine proteases and their inhibitors: drug discovery targets?

Papain-like lysosomal cysteine proteases are processive and digestive enzymes that are expressed in organisms from bacteria to humans. Increasing knowledge about the physiological and pathological roles of cysteine proteases is bringing them into the focus of drug discovery research. These proteases have rather short active-site clefts, comprising three well defined substrate-binding subsites (S2, S1 and S1') and additional broad binding areas (S4, S3, S2' and S3'). The geometry of the active site distinguishes cysteine proteases from other protease classes, such as serine and aspartic proteases, which have six and eight substrate-binding sites respectively. Exopeptidases (cathepsins B, C, H and X), in contrast with endopeptidases (such as cathepsins L, S, V and F), possess structural features that facilitate the binding of N- and C-terminal groups of substrates into the active-site cleft. Other than a clear preference for free chain termini in the case of exopeptidases, the substrate-binding sites exhibit no strict specificities. Instead, their subsite preferences arise more from the specific exclusion of substrate types. This presents a challenge for the design of inhibitors to target a specific cathepsin: only the cumulative effect of an assembly of inhibitor fragments will bring the desired result.

Engineering cyclodextrin glycosyltransferase into a starch hydrolase with a high exo-specificity

Cyclodextrin glycosyltransferase (CGTase) enzymes from various bacteria catalyze the formation of cyclodextrins from starch. The Bacillus stearothermophilus maltogenic α-amylase (G2-amylase is structurally very similar to CGTases, but converts starch into maltose. Comparison of the three-dimensional structures revealed two large differences in the substrate binding clefts. (i) The loop forming acceptor subsite +3 had a different conformation, providing the G2-amylase with more space at acceptor subsite +3, and (ii) the G2-amylase contained a five-residue amino acid insertion that hampers substrate binding at the donor subsites −3/−4 (Biochemistry, 38 (1999) 8385). In an attempt to change CGTase into an enzyme with the reaction and product specificity of the G2-amylase, which is used in the bakery industry, these differences were introduced into Thermoanerobacterium thermosulfurigenes CGTase. The loop forming acceptor subsite +3 was exchanged, which strongly reduced the cyclization activity, however, the product specificity was hardly altered. The five-residue insertion at the donor subsites drastically decreased the cyclization activity of CGTase to the extent that hydrolysis had become the main activity of enzyme. Moreover, this mutant produces linear products of variable sizes with a preference for maltose and had a strongly increased exo-specificity. Thus, CGTase can be changed into a starch hydrolase with a high exo-specificity by hampering substrate binding at the remote donor substrate binding subsites.

Specificity determinants of human cathepsin s revealed by crystal structures of complexes.

Cathepsin S, a lysosomal cysteine protease of the papain superfamily, has been implicated in the preparation of MHC class II alphabeta-heterodimers for antigen presentation to CD4+ T lymphocytes and is considered a potential target for autoimmune-disease therapy. Selective inhibition of this enzyme may be therapeutically useful for attenuating the hyperimmune responses in a number of disorders. We determined the three-dimensional crystal structures of human cathepsin S in complex with potent covalent inhibitors, the aldehyde inhibitor 4-morpholinecarbonyl-Phe-(S-benzyl)Cys-Psi(CH=O), and the vinyl sulfone irreversible inhibitor 4-morpholinecarbonyl-Leu-Hph-Psi(CH=CH-SO(2)-phenyl) at resolutions of 1.8 and 2.0 A, respectively. In the structure of the cathepsin S-aldehyde complex, the tetrahedral thiohemiacetal adduct favors the S-configuration, in which the oxygen atom interacts with the imidazole group of the active site His164 rather than with the oxyanion hole. The present structures provide a detailed map of noncovalent intermolecular interactions established in the substrate-binding subsites S3 to S1' of cathepsin S. In the S2 pocket, which is the binding affinity hot spot of cathepsin S, the Phe211 side chain can assume two stable conformations that accommodate either the P2-Leu or a bulkier P2-Phe side chain. This structural plasticity of the S2 pocket in cathepsin S explains the selective inhibition of cathepsin S over cathepsin K afforded by inhibitors with the P2-Phe side chain. Comparison with the structures of cathepsins K, V, and L allows delineation of local intermolecular contacts that are unique to cathepsin S.

Lysosomal cysteine proteases (cathepsins): promising drug targets.

Papain-like lysosomal cysteine proteases are processive and digestive enzymes expressed in organisms from bacteria to humans. Their ubiquity alone makes them potential drug targets, with the assumption that appropriate specificities may be achieved. These enzymes have rather short active-site clefts, comprising three well defined substrate-binding subsites (S2, S1 and S1') and additionally have comparatively broad binding areas (S4, S3, S2', S3'). This geometry distinguishes them from other protease classes, such as serine and aspartic proteases, with six and eight substrate-binding sites, respectively. Exopeptidases (cathepsins B, C, H and X), in contrast to endopeptidases (such as cathepsins L, S, V and F), possess structural features that facilitate binding of N- and C-terminal groups of substrates in the active-site cleft. Other than a clear preference for free chain termini in the case of exopeptidases, the substrate-binding sites exhibit no strict specificities. Instead, their subsite preferences arise more from specific exclusions of substrate type. This presents a challenge for the design of inhibitors to target a specific cathepsin: only the cumulative effect of an assembly of inhibitor fragments can produce the desired result. The small number of papain-like lysosomal cysteine proteases (11 human enzymes are known) and the small number of substrate-binding sites calls for a innovative and empirical approach.

Family 19 chitinase from rice (Oryza sativa L.): substrate-binding subsites demonstrated by kinetic and molecular modeling studies.

A family 19 chitinase (OsChia1c, class I) from rice, Oryza sativa L., and its chitin-binding domain-truncated mutant (OsChia1c deltaCBD, class II) were produced by the Pichia expression system, and the hydrolytic mechanism toward N-acetylglucosamine hexasaccharide [(GlcNAc)6] was investigated by HPLC analysis of the reaction products. The profile of the time-course of (GlcNAc)6 degradation obtained by OsChia1c was identical to that obtained by OsChia1c deltaCBD, indicating that the chitin-binding domain does not significantly participate in oligosaccharide hydrolysis. From the theoretical analysis of the reaction time-course of OsChia1c deltaCBD, the free energy changes of sugar residue binding were estimated to be -0.4, -4.7, +3.4, -0.5, -2.3, and -1.0 kcal/mol for the individual subsites of (-3), (-2), (-1), (+1), (+2), and (+3), respectively. The hexasaccharide substrate appears to bind to the enzyme through interactions at the high-affinity sites, (-2) and (+2), and the sugar residues at both ends more loosely bind to the corresponding subsites, (-3) and (+3). The docking study of (GlcNAc)6 with the modeled structure of OsChia1c deltaCBD supported the subsite structure estimated from the experimental time-course of hexasaccharide degradation. Since the class II chitinase from barley seeds was reported to possess a similar subsite structure from (-3) to (+3) and a similar free energy distribution, substrate-binding mode of plant chitinases of this class would be similar to each other.

Aspergillus aculeatus beta-1,4-galactanase: substrate recognition and relations to other glycoside hydrolases in clan GH-A.

The three-dimensional structure of Aspergillus aculeatus beta-1,4-galactanase (AAGAL), an enzyme involved in pectin degradation, has been determined by multiple isomorphous replacement to 2.3 and 1.8 A resolution at 293 and 100 K, respectively. It represents the first known structure for a polysaccharidase with this specificity and for a member of glycoside hydrolase family 53 (GH-53). The enzyme folds into a (beta/alpha)(8) barrel with the active site cleft located at the C-terminal side of the barrel consistent with the classification of GH-53 in clan GH-A, a superfamily of enzymes with common fold and catalytic machinery but diverse specificities. Putative substrate-enzyme interactions were elucidated by modeling of beta-1,4-linked galactobioses into the possible substrate binding subsites. The structure and modeling studies identified five potential subsites for the binding of galactans, of which one is a pocket suited for accommodating the arabinan side chain in arabinogalactan, one of the natural substrates. A comparison with the substrate binding grooves of other Clan GH-A enzymes suggests that shape complementarity is crucial in determining the specificity of polysaccharidases.

Monomeric Structures of the Zymogen and Active Catalytic Domain of Complement Protease C1r: Further Insights into the C1 Activation Mechanism

C1r is the serine protease (SP) that mediates autoactivation of C1, the complex that triggers the classical complement pathway. We have determined the crystal structure of two fragments from the human C1r catalytic domain, each encompassing the second complement control protein (CCP2) module and the SP domain. The wild-type species has an active structure, whereas the S637A mutant is a zymogen. The structures reveal a restricted hinge flexibility of the CCP2-SP interface, and both are characterized by the unique α-helical conformation of loop E. The zymogen activation domain exhibits high mobility, and the active structure shows a restricted access to most substrate binding subsites. Further implications relevant to the C1r self-activation process are derived from protein-protein interactions in the crystals.

The 1.62 Å structure of Thermoascus aurantiacus endoglucanase: completing the structural picture of subfamilies in glycoside hydrolase family 5

The crystal structure of Thermoascus aurantiacus endoglucanase (Cel5A), a family 5 glycoside hydrolase, has been determined to 1.62 Å resolution by multiple isomorphous replacement with anomalous scattering. It is the first report of a structure in the subfamily to which Cel5A belongs. Cel5A consists solely of a catalytic module with compact eight-fold β/α barrel architecture. The length of the tryptophan-rich substrate binding groove suggests the presence of substrate binding subsites −4 to +3. Structural comparison shows that two glycines are completely conserved in the family, in addition to the two catalytic glutamates and six other conserved residues previously identified. Gly 44 in particular is part of a type IV C-terminal helix capping motif, whose disruption is likely to affect the position of an essential conserved arginine. One aromatic residue (Trp 170 in Cel5A), not conserved in term of sequence, is nonetheless spatially conserved in the substrate binding groove. Its role might be to force the bend that occurs in the polysaccharide chain on binding, thus favoring substrate distortion at subsite −1.