Subsite mapping of enzymes. Application of the depolymerase computer model to two alpha-amylases.

In the preceding paper (Allen and Thoma, 1976) we developed a depolymerase computer model, which uses a minimization routine to establish a subsite map for a depolymerase. In the present paper we show how the model is applied to experimental data for two alpha-amylases. Michaelis parameters and bond-cleavage frequencies for substrates of chain lengths up to twelve glucosyl units have been reported for Bacillus amyloliquefaciens, and a subsite map has been proposed for this enzyme [Thoma et al. (1971) J. Biol. Chem. 246, 5621-5635]. By applying the computer model to the experimental data, we have arrived at a ten-subsite map. We find that a significant improvement in this map is achieved by allowing the hydrolytic rate coefficient to vary as a function of the number of occupied subsites comprising the enzyme-binding region. The bond-cleavage frequencies, the enzyme is found to have eight subsites. A partial subsite map is arrived at, but the entire binding region cannot be mapped because Michaelis parameters are complicated by transglycosylation reactions. The hydrolytic rate coefficients for this enzyme are not constant.     

Subsite mapping of enzymes. Double inhibition studies

Our earlier subsite mapping studies led us to believe that ground state distortion occurred when a glucopyranoside ring filled the site which held the substrate monomer unit transferred to water during hydrolysis. We tested this hypothesis by performing double inhibitor studies on two amylases (EC 3.2.1.1) of bacterial origin. A general theory for multiple inhibition of this type is developed and applied to these two enzymes. Our data are consistent with the hypothesis that ground state strain occurs when substrates are bound to carbohydrates. An explanation is offered to account for the fact that monomers give strictly competitive inhibition patterns. The subsite model predicts that noncompetitive or mixed inhibition patterns can occur.     

Subsite mapping of the binding region of alpha-amylases with a computer program.

A computer program has been evaluated for subsite map calculations of depolymerases. The program runs in windows and uses the experimentally determined bond cleavage frequencies (BCFs) for determination of the number of subsites, the position of the catalytic site and for calculation of subsite binding energies. The apparent free energy values were optimized by minimization of the differences of the measured and calculated BCF data. The program called suma (SUbsite Mapping of alpha-Amylases) is freely available for research and educational purposes via the Internet (E-mail: gyemant@tigris.klte.hu). The advantages of this program are demonstrated through alpha-amylases of different origin, e.g. porcine pancreatic alpha-amylase (PPA) studied in our laboratory, in addition to barley and rice alpha-amylases published in the literature. Results confirm the popular 'five subsite model' for PPA with three glycone and two aglycone binding sites. Calculations for barley alpha-amylase justify the '6 + 2 + (1) model' prediction. The binding area of barley alpha-amylase is composed of six glycone, two aglycone binding sites followed by a barrier subsite at the reducing end of the binding site. Calculations for rice alpha-amylase represent an entirely new map with a '(1) + 2 + 5 model', where '(1)' is a barrier subsite at the nonreducing end of the binding site and there are two glycone and five aglycone binding sites. The rice model may be reminiscent of the action of the bacterial maltogenic amylase, that is, suggesting an exo-mechanism for this enzyme.     

Subsite mapping of enzymes. Use of subsite map to simulate complete time course of hydrolysis of a polymeric substrate.

In this paper we extend our earlier work on subsite mapping and show that our model for depolymerase action can be used to accurately predict product ratios vs the extent of reaction when a polymer is hydrolyzed. The experimental product ratios for Bacillus amyloliquefaciens α-amylase acting on reducing end-labeled ¹⁴C-maltodextrins ranging in chain length 3 to 10 are reported. These data and Michaelis parameters are used with a depolymerase computer model (J. D. Allen, 1977, Ph.D. thesis, University of Arkansas; J. D. Allen and J. A. Thoma, 1976, Biochem. J.159, 105) to compute an optimized subsite map. The depolymerase computer model generates a 10-subsite map for B. amyloliquefaciens α-amylase with the catalytic site located to the left of subsite 7. The binding affinities of the subsites are then used as the sole input in another computer program to quantitatively predict the mole fraction of products vs the extent of hydrolysis for substrates of varying chain length. Excellent agreement is obtained between the computed and experimental data for seven maltodextrins examined.     

Subsite mapping of Aspergillus niger endopolygalacturonase II by site-directed mutagenesis

To assess the subsites involved in substrate binding in Aspergillus niger endopolygalacturonase II, residues located in the potential substrate binding cleft stretching along the enzyme from the N to the C terminus were subjected to site-directed mutagenesis. Mutant enzymes were characterized with respect to their kinetic parameters using polygalacturonate as a substrate and with respect to their mode of action using oligogalacturonates of defined length (n = 3-6). In addition, the effect of the mutations on the hydrolysis of pectins with various degrees of esterification was studied. Based on the results obtained with enzymes N186E and D282K it was established that the substrate binds with the nonreducing end toward the N terminus of the enzyme. Asn(186) is located at subsite -4, and Asp(282) is located at subsite +2. The mutations D183N and M150Q, both located at subsite -2, affected catalysis, probably mediated via the sugar residue bound at subsite -1. Tyr(291), located at subsite +1 and strictly conserved among endopolygalacturonases appeared indispensable for effective catalysis. The mutations E252A and Q288E, both located at subsite +2, showed only slight effects on catalysis and mode of action. Tyr(326) is probably located at the imaginary subsite +3. The mutation Y326L affected the stability of the enzyme. For mutant E252A, an increased affinity for partially methylesterified substrates was recorded. Enzyme N186E displayed the opposite behavior; the specificity for completely demethylesterified regions of substrate, already high for the native enzyme, was increased. The origin of the effects of the mutations is discussed.     

Subsite mapping of human salivary alpha-amylase and the mutant Y151M

This study characterizes the substrate-binding sites of human salivary alpha-amylase (HSA) and its Y151M mutant. It describes the first subsite maps, namely, the number of subsites, the position of cleavage sites and apparent subsite energies. The product pattern and cleavage frequencies were determined by high-performance liquid chromatography, utilizing a homologous series of chromophore-substituted maltooligosaccharides of degree of polymerization 3-10 as model substrates. The binding region of HSA is composed of four glycone and three aglycone-binding sites, while that of Tyr151Met is composed of four glycone and two aglycone-binding sites. The subsite maps show that Y151M has strikingly decreased binding energy at subsite (+2), where the mutation has occurred (-2.6 kJ/mol), compared to the binding energy at subsite (+2) of HSA (-12.0 kJ/mol).     

Subsite mapping of Aspergillus niger glucoamylases I and II with malto- and isomaltooligosaccharides

Glucoamylase, industrially derived from Aspergillus niger, was chromatographically separated into forms I and II and purified to near homogeneity. Preparations were proved to be free of D-glucosyltransferase by electrophoretic and differential inhibition tests. Maximum rates and Michaelis constants were obtained for both glucoamylases I and II with maltooligosaccharides from maltose to maltoheptaose and with isomaltooligosac-charides from isomaltose to isomaltohexaose. Subsite maps were calculated from these kinetic data and were not significantly different for the two forms. Subsites in both forms had lower affinities for D-glucosyl residues contained in isomaltooligosaccharides than for D-glucosyl residues in maltooligosaccharides.     

Towards 3-D modelling of epithelia by computer simulation.

The present work is a preliminary step towards dynamic 3-D modelling by computer graphics simulation of the structure of normal and pathological epithelia, using an expert system. In its present state, Esexsy (Epithelium Simulation by EXpert SYstem) allows the construction, through iterative steps, of a simple 3-D representation of the nasal epithelium, based on the positions, sizes and shapes of nuclei. The iterative process is based on statistical comparisons between distributions of parameter values calculated from real (2-D) histological sections and those issued from an equivalent computer 'section' through the simulated 3-D image. We show the results of attempts at simulating normal, metaplastic and dysplastic states of the nasal epithelium, the latter two being characterized by a progressive architectural disorganization, accompanied by nuclear size/shape alterations. The representation takes into account the size, shape, orientation and spatial arrangement of nuclei, with one or several layers from the basal lamina to the lumen. A modified Poisson point process is used at present to position the nuclei, which are modelled by bi-axial spheroids (from prolate to oblate through spherical), with random orientation and size/shape deviations. It should be possible to use the same computer program to simulate other types of epithelia and to achieve increasingly realistic representations by incorporating, notably, nuclear deformations and chromatin texture.     

Human leukocyte cathepsin G. Subsite mapping with 4-nitroanilides, chemical modification, and effect of possible cofactors

The extended substrate binding site of cathepsin G from human leukocytes has been mapped by using a series of peptide 4-nitroanilide substrates. The enzyme has a significant preference for substrates with a P1 Phe over those with the other aromatic amino acids Tyr and Trp. The S2 subsite was mapped with the substrates Suc-Phe-AA-Phe-NA where AA was 13 of the 20 amino acid residues commonly found in proteins. The best residues were Pro and Met. The S3 subsite was mapped with the sequence Suc-AA-Pro-Phe-NA by using 14 different amino acid residues for AA. The two best residues were the isosteric Val and Thr. No significant improvement in reactivity was obtained by extending the substrate to include seven different P4 residues. The kinetic parameters for cathepsin G are significantly slower than those for many other serine proteases. Changes in the reaction conditions and addition of possible cofactors or ligands were in general found to have little effect on the enzymatic activity, while chemical modifications and proteolysis destroyed the activity of cathepsin G. Cathepsin G hydrolyzed peptides containing model desmosine residues and prefers the hydrophobic picolinoyllysine derivative over lysine by substantial margins at both the S4 and S2 subsites but will not tolerate it at S3. Substrates with sequences related to the cathepsin G cleavage site in angiotensin I and angiotensinogen, and the reactive site of alpha 1-antichymotrypsin, were hydrolyzed effectively by enzyme, but with unexceptional rates. Our results indicate that the natural substrate(s) and function(s) of cathepsin G still remain to be discovered.     

Subsite mapping of an acidic amino acid-specific endopeptidase from Streptomyces griseus, GluSGP, and protease V8.

The substrate specificities of an acidic amino acid-specific endopeptidase of Streptomyces griseus, GluSGP, and protease V8 [EC 3.4.21.19] were investigated with peptide p-nitroanilide substrates which have a Glu residue at the P1 position. GluSGP and protease V8 favored Pro and Leu residues at S2, respectively, while the S3 subsite of GluSGP preferred Phe over either Ala or Leu. The S3 subsite of protease V8 preferred Leu over either Ala or Phe. The best substrates for GluSGP and for protease V8 were Boc-Ala-Phe-Pro-Glu-pNA with a Km value of 0.41 mM (0.1 M Tris-HCl, pH 8.8) and Boc-Ala-Leu-Leu-Glu-pNA with a Km value of 0.25 mM (0.1 M phosphate, pH 7.8), respectively. The kcat/Km values for these substrates obtained with GluSGP were about one hundred to twenty thousand times larger than those obtained with protease V8. Protease V8 exhibited a single optimal pH of around 8 for the hydrolysis of Boc-Ala-Ala-Leu-Glu-pNA and Boc-Ala-Leu-Leu-Asp-pNA.     

Subsite mapping of the human pancreatic alpha-amylase active site through structural, kinetic, and mutagenesis techniques

We report a multifaceted study of the active site region of human pancreatic alpha-amylase. Through a series of novel kinetic analyses using malto-oligosaccharides and malto-oligosaccharyl fluorides, an overall cleavage action pattern for this enzyme has been developed. The preferred binding/cleavage mode occurs when a maltose residue serves as the leaving group (aglycone sites +1 and +2) and there are three sugars in the glycon (-1, -2, -3) sites. Overall it appears that five binding subsites span the active site, although an additional glycon subsite appears to be a significant factor in the binding of longer substrates. Kinetic parameters for the cleavage of substrates modified at the 2 and 4' ' positions also highlight the importance of these hydroxyl groups for catalysis and identify the rate-determining step. Further kinetic and structural studies pinpoint Asp197 as being the likely nucleophile in catalysis, with substitution of this residue leading to an approximately 10(6)-fold drop in catalytic activity. Structural studies show that the original pseudo-tetrasaccharide structure of acarbose is modified upon binding, presumably through a series of hydrolysis and transglycosylation reactions. The end result is a pseudo-pentasaccharide moiety that spans the active site region with its N-linked "glycosidic" bond positioned at the normal site of cleavage. Interestingly, the side chains of Glu233 and Asp300, along with a water molecule, are aligned about the inhibitor N-linked glycosidic bond in a manner suggesting that these might act individually or collectively in the role of acid/base catalyst in the reaction mechanism. Indeed, kinetic analyses show that substitution of the side chains of either Glu233 or Asp300 leads to as much as a approximately 10(3)-fold decrease in catalytic activity. Structural analyses of the Asp300Asn variant of human pancreatic alpha-amylase and its complex with acarbose clearly demonstrate the importance of Asp300 to the mode of inhibitor binding.     

Subsite mapping of purified glucoamylase I, II, III produced by Arthrobotrys amerospora ATCC 34468

Arthrobotrys amerospora ATCC 34468 produced glucoamylase in a medium containing maize starch as carbon source. On native PAGE, crude glucoamylase showed three isoenzymes which were designated as Glu I, Glu II, Glu III according to their electrophoretic mobility. These were purified by column chromatography techniques. The energy of binding for each glucoamylase was calculated using Hiromi's kinetic based calculation. At subsite 1, the binding energies for Glu I, II and III were found to be negative.     

Cognate ligand domain mapping for enzymes

Here, we present an automatic assignment of potential cognate ligands to domains of enzymes in the CATH and SCOP protein domain classifications on the basis of structural data available in the wwPDB. This procedure involves two steps; firstly, we assign the binding of particular ligands to particular domains; secondly, we compare the chemical similarity of the PDB ligands to ligands in KEGG in order to assign cognate ligands. We find that use of the Enzyme Commission (EC) numbers is necessary to enable efficient and accurate cognate ligand assignment. The PROCOGNATE database currently has cognate ligand mapping for 3277 (4118) protein structures and 351 (302) superfamilies, as described by the CATH and (SCOP) databases, respectively. We find that just under half of all ligands are only and always bound by a single domain, with 16% bound by more than one domain and the remainder of the ligands showing a variety of binding modes. This finding has implications for domain recombination and the evolution of new protein functions. Domain architecture or context is also found to affect substrate specificity of particular domains, and we discuss example cases. The most popular PDB ligands are all found to be generic components of crystallisation buffers, highlighting the non-cognate ligand problem inherent in the PDB. In contrast, the most popular cognate ligands are all found to be universal cellular currencies of reducing power and energy such as NADH, FADH2 and ATP, respectively, reflecting the fact that the vast majority of enzymatic reactions utilise one of these popular co-factors. These ligands all share a common adenine ribonucleotide moiety, suggesting that many different domain superfamilies have converged to bind this chemical framework.     

A Novel PHB Depolymerase from a Thermophilic Streptomyces Sp.

A novel PHB depolymerase from a thermophilic Streptomyces sp. MG was purified to homogeneity by hydrophobic interaction chromatography and gel filtration. The molecular mass of the purified enzyme was 43 kDa as determined by size exclusion chromatography and 41 kDa by SDS-PAGE. The optimum pH and temperature were 8.5 and 60 °C respectively. The enzyme was stable at 50 °C and from pH 6.5–8.5. The enzyme hydrolyzed not only bacterial polyesters, i.e. poly(3-hydroxybutyric acid and poly(3-hydroxybutyrate-co-3-hydroxyvalerate), but also synthetic, aliphatic polyesters such as polypropiolactone, poly(ethylene adipate) and poly(ethylene succinate). Revisions requested 9 November 2005; Revisions received 12 December 2005     

Structure of donor side components in photosystem II predicted by computer modelling.

Thirty-one and eleven sequences for the photosystem II reaction centre proteins D1 and D2 respectively, were compared to identify conserved single amino acid residues and regions in the sequences. Both proteins are highly conserved. One important difference is that the lumenal parts of the D1 protein are more conserved than the corresponding parts in the D2 protein. The three-dimensional structures around the electron donors tyrosineZ and tyrosineD on the oxidizing side of photosystem II have been predicted by computer modelling using the photosynthetic reaction centre from purple bacteria as a framework. In the model the tyrosines occupy two cavities close to the lumenal surface of the membrane. They are symmetrically arranged around the primary donor P680 and the distances between the centre of the tyrosines and the closest Mg ion in P680 are around 14 A. Both tyrosineZ and tyrosineD are suggested to form a hydrogen bond with histidine 190 from the loop connecting helices C and D in the D1 and D2 proteins, respectively. The Mn cluster in the oxygen evolving complex has been localized by using known and estimated distances from the tyrosine radicals. It is suggested that a binding region for the Mn cluster is constituted by the lumenal ends of helices A and B and the loop connecting them in the D1 protein. This part of the D1 protein contains a large number of strictly conserved carboxylic acid residues and histidines which could participate in the Mn binding. There is little probability that the Mn cluster binds on the lumenal surface of the D2 protein.     

Origin of Polysaccharide Depolymerase Associated with Bacteriophage Infection

Analyses, by construction of phage growth curves, indicated that the polysaccharide depolymerase was synthesized by Pseudomonas aeruginosa strains B and BI after infection with phage 2. The kinetics of biosynthesis of the depolymerase were found to parallel closely the rate of formation of phage-directed virions, and alterations in the experimental conditions of infection were reflected by alterations in the production of enzyme. Infection with other Pseudomonas phages, 84 and 1197, did not result in the synthesis of depolymerase. The enzyme was not detectable in uninfected cultures, and no evidence was obtained for the existence of inhibitors or activators of enzyme activity in extracts of uninfected or infected cells. The results of experiments employing chloramphenicol or an auxotorphic mutant (BI arg(-)) suggested that protein synthesis de novo was essential for production of the enzyme. Various mutants of phage 2 (pdp(1), pdp(2)), which alter the synthesis of the polysaccharide depolymerase, have been isolated. These experimental results strongly support the role of the phage genome in the synthesis of this enzyme.