Search

Present estimates suggest that of the 359 million tons of plastics produced annually worldwide1, 150–200 million tons accumulate in landfill or in the natural environment2. Poly(ethylene terephthalate) (PET) is the most abundant polyester plastic, with almost 70 million tons manufactured annually worldwide for use in textiles and packaging3. The main recycling process for PET, via thermomechanical means, results in a loss of mechanical properties4. Consequently, de novo synthesis is preferred and PET waste continues to accumulate. With a high ratio of aromatic terephthalate units—which reduce chain mobility—PET is a polyester that is extremely difficult to hydrolyse5. Several PET hydrolase enzymes have been reported, but show limited productivity6,7. Here we describe an improved PET hydrolase that ultimately achieves, over 10 hours, a minimum of 90 per cent PET depolymerization into monomers, with a productivity of 16.7 grams of terephthalate per litre per hour (200 grams per kilogram of PET suspension, with an enzyme concentration of 3 milligrams per gram of PET). This highly efficient, optimized enzyme outperforms all PET hydrolases reported so far, including an enzyme8,9 from the bacterium Ideonella sakaiensis strain 201-F6 (even assisted by a secondary enzyme10) and related improved variants11–14 that have attracted recent interest. We also show that biologically recycled PET exhibiting the same properties as petrochemical PET can be produced from enzymatically depolymerized PET waste, before being processed into bottles, thereby contributing towards the concept of a circular PET economy. Computer-aided engineering produces improvements to an enzyme that breaks down poly(ethylene terephthalate) (PET) into its constituent monomers, which are used to synthesize PET of near-petrochemical grade that can be further processed into bottles.

Sequence-specific targeting of double-stranded DNA and non-coding RNA via triple-helix-forming peptide nucleic acids (PNAs) has attracted considerable attention in therapeutic, diagnostic and nanotechnological fields. An E-base (3-oxo-2,3-dihydropyridazine), attached to the polyamide backbone of a PNA Hoogsteen strand by a side-chain linker molecule, is typically used in the hydrogen bond recognition of the 4-oxo group of thymine and uracil nucleic acid bases in the major groove. We report on the application of quantum chemical computational methods, in conjunction with spatial constraints derived from the experimental structure of a homopyrimidine PNA·DNA-PNA hetero-triplex, to investigate the influence of linker flexibility on binding interactions of the E-base with thymine and uracil bases in geometry-optimised model systems. Hydrogen bond formation between the N2 E-base atom and target pyrimidine base 4-oxo groups in model systems containing a β-alanine linker (J Am Chem Soc 119:11116, 1997) was found to incur significant internal strain energy and the potential disruption of intra-stand aromatic base stacking interactions in an oligomeric context. In geometry-optimised model systems containing a 3-trans olefin linker (Bioorg Med Chem Lett 14:1551, 2004) the E-base swung out away from the target pyrimidine bases into the solvent. These findings are in qualitative agreement with calorimetric measurements in hybridisation experiments at T-A and U-A inversion sites. In contrast, calculations on a novel 2-cis olefin linker design indicate that it could permit simultaneous E-base hydrogen bonding with the thymine 4-oxo group, circumvention and solvent screening of the thymine 5-methyl group, and maintenance of triplex intra-stand base stacking interactions.

Shortcomings in the definition of effective free-energy surfaces of proteins are recognized to be a major contributory factor responsible for the low success rates of existing automated methods for computational protein design (CPD). The formulation of an atomistic statistically effective energy function (SEEF) suitable for a wide range of CPD applications and its derivation from structural data extracted from protein domains and protein-ligand complexes are described here. The proposed energy function comprises nonlocal atom-based and local residue based SEEFs, which are coupled using a novel atom connectivity number factor to scale short-range, pairwise, nonbonded atomic interaction energies and a surface-area-dependent cavity energy term. This energy function was used to derive additional SEEFs describing the unfolded-state ensemble of any given residue sequence based on computed average energies for partially or fully solvent-exposed fragments in regions of irregular structure in native proteins. Relative thermal stabilities of 97 T4 bacteriophage lysozyme mutants were predicted from calculated energy differences for folded and unfolded states with an average unsigned error (AUE) of 0.84 kcal mol(-1) when compared to experiment. To demonstrate the utility of the energy function for CPD, further validation was carried out in tests of its capacity to recover cognate protein sequences and to discriminate native and near-native protein folds, loop conformers, and small-molecule ligand binding poses from non-native benchmark decoys. Experimental ligand binding free energies for a diverse set of 80 protein complexes could be predicted with an AUE of 2.4 kcal mol(-1) using an additional energy term to account for the loss in ligand configurational entropy upon binding. The atomistic SEEF is expected to improve the accuracy of residue-based coarse-grained SEEFs currently used in CPD and to extend the range of applications of extant atom-based protein statistical potentials.

A synthetic pathway for the production of 2,4-dihydroxybutyric acid from homoserine, composed of two consecutive enzymatic reaction steps has been recently reported. An important step in this pathway consists in the reduction of 2-keto-4-hydroxybutyrate (OHB) into (L)-dihydroxybutyrate (DHB), by an enzyme with OHB reductase activity. In this study, we used a rational approach to engineer an OHB reductase by using the cytosolic (L)-malate dehydrogenase from Escherichia coli (Ec-Mdh) as the template enzyme. Structural analysis of (L)-malate dehydrogenase and (L)-lactate dehydrogenase enzymes acting on sterically cognate substrates revealed key residues in the substrate and co-substrate binding sites responsible for substrate discrimination. Accordingly, amino acid changes were introduced in a step-wise manner into these regions of the protein. This rational engineering led to the production of a Ec-Mdh-5E variant (I12V/R81A/M85E/G179D/D86S) with a turnover number (k cat ) on OHB that was increased by more than 2,000 fold (from 0.03 up to 65.0 s -1 ), which turned out to be 7 fold higher than that on its natural substrate oxaloacetate. Further kinetic analysis revealed the engineered enzyme to possess comparable catalytic efficiencies (k cat /K m ) between natural and synthetic OHB substrates (84 and 31 s -1 mM -1 , respectively). Shake-flask cultivation of an homoserine-overproducing E. coli strain expressing this improved OHB reductase together with a transaminase encoded by aspC able to convert homoserine to OHB resulted in 89 % increased DHB production as compared to our previous report using a E. coli host strain expressing an OHB reductase derived from the lactate dehydrogenase A of Lactococcus lactis

An end-point ADP/NAD(+) acid/alkali assay procedure, directly applicable to library screening of any type of ATP-utilising/ADP producing enzyme activity, was implemented. Typically, ADP production is coupled to NAD(+) co-enzyme formation by the conventional addition of pyruvate kinase and lactate dehydrogenase. Transformation of enzymatically generated NAD(+) into a photometrically active alkali derivative product is then achieved through the successive application of acidic/alkali treatment steps. The assay was successfully miniaturized to search for malate kinase activity in a structurally-guided library of LysC aspartate kinase variants comprising 6,700 clones. The screening procedure enabled the isolation of nine positive variants showing novel kinase activity on (L)-malate, the best mutant, LysC V115A: E119S:E434V exhibited strong substrate selectivity for (L)-malate compared to (L)-aspartate with a (k(cat)/K-m)(malate)/(k(cat)/K-m)(aspartate) ratio of 86. Double mutants V115A:E119S, V115A:E119C and E119S:E434V were constructed to further probe the origins of stabilising substrate binding energy gains for (L)-malate due to mutation. The introduction of less sterically hindering side-chains in engineered enzymes carrying E119S and V115A mutations increases the effective volume available for substrate binding in the catalytic pocket. Improved binding of the (L)-malate substrate may be assisted by less hindered movement of the Phe184 aromatic side-chain. Additional favourable long-range electostatic effects on binding arising from the E434V surface mutation are conditionally dependent upon the presence of the V115A mutation close to Phe184 in the active-site.

Relative amino acid residue solvent accessibility values allow the quantitative comparison of atomic solvent-accessible surface areas in different residue types and physical environments in proteins and in protein structural alignments. Geometry-optimised tri-peptide structures in extended solvent-exposed reference conformations have been obtained for 43 amino acid residue types at a high level of quantum chemical theory. Significant increases in side-chain solvent accessibility, offset by reductions in main-chain atom solvent exposure, were observed for standard residue types in partially geometry-optimised structures when compared to non-minimised models built from identical sets of proper dihedral angles abstracted from the literature. Optimisation of proper dihedral angles led most notably to marked increases of up to 54% in proline main-chain atom solvent accessibility compared to literature values. Similar effects were observed for fully-optimised tri-peptides in implicit solvent. The relief of internal strain energy was associated with systematic variation in N, Cα and Cβ atom solvent accessibility across all standard residue types. The results underline the importance of optimisation of ‘hard’ degrees of freedom (bond lengths and valence bond angles) and improper dihedral angle values from force field or other context-independent reference values, and impact on the use of standardised fixed internal co-ordinate geometry in sampling approaches to the determination of absolute values of protein amino acid residue solvent accessibility. Quantum chemical methods provide a useful and accurate alternative to molecular mechanics methods to perform energy minimisation of peptides containing non-standard (chemically modified) amino acid residues frequently present in experimental protein structure data sets, for which force field parameters may not be available. Reference tri-peptide atomic co-ordinate sets including hydrogen atoms are made freely available.

2,4-Dihydroxybutyric acid (DHB) is a molecule with considerable potential as a versatile chemical synthon. Notably, it may serve as a precursor for chemical synthesis of the methionine analogue 2-hydroxy-4-(methylthio)butyrate, thus, targeting a considerable market in animal nutrition. However, no natural metabolic pathway exists for the biosynthesis of DHB. Here we have therefore conceived a three-step metabolic pathway for the synthesis of DHB starting from the natural metabolite malate. The pathway employs previously unreported malate kinase, malate semialdehyde dehydrogenase and malate semialdehyde reductase activities. The kinase and semialdehyde dehydrogenase activities were obtained by rational design based on structural and mechanistic knowledge of candidate enzymes acting on sterically cognate substrates. Malate semialdehyde reductase activity was identified from an initial screening of several natural enzymes, and was further improved by rational design. The pathway was expressed in a minimally engineered Escherichia coli strain and produces 1.8 g l−1 DHB with a molar yield of 0.15., 2,4-Dihydroxybutyric acid has potential to be a precursor to a range of industrially important products, however a natural metabolic pathway for its synthesis does not exist. Here the authors rationally design a synthetic pathway in E. coli by engineering enzymes from malate metabolism.

The amylosucrase from Neisseria polysaccharea is a transglucosidase from the GH13 family of glycoside-hydrolases that naturally catalyzes the synthesis of α-glucans from the widely available donor sucrose. Interestingly, natural molecular evolution has modeled a dense hydrogen bond network at subsite −1 responsible for the specific recognition of sucrose and conversely, it has loosened interactions at the subsite +1 creating a highly promiscuous subsite +1. The residues forming these subsites are considered to be likely involved in the activity as well as the overall stability of the enzyme. To assess their role, a structure-based approach was followed to reshape the subsite −1. A strategy based on stability change predictions, using the FoldX algorithm, was considered to identify the best candidates for site-directed mutagenesis and guide the construction of a small targeted library. A miniaturized purification protocol was developed and both mutant stability and substrate promiscuity were explored. A range of 8°C between extreme melting temperature values was observed and some variants were able to synthesize series of oligosaccharides with distributions differing from that of the parental enzyme. The crucial role of subsite −1 was thus highlighted and the biocatalysts generated can now be considered as starting points for further engineering purposes.

The POLYFIT rigid-body algorithm for automated global pairwise and multiple protein structural alignment is presented. Smith-Waterman local alignment is used to establish a set of seed equivalences that are extended using Needleman-Wunsch dynamic programming techniques. Structural and functional interaction constraints provided by evolution are encoded as one-dimensional residue physical environment strings for alignment of highly structurally overlapped protein pairs. Local structure alignment of more distantly related pairs is carried out using rigid-body conformational matching of 15-residue fragments, with allowance made for less stringent conformational matching of metal-ion and small molecule ligand-contact, disulphide bridge, and cis-peptide correspondences. Protein structural plasticity is accommodated through the stepped adjustment of a single empirical distance parameter value in the calculation of the Smith-Waterman dynamic programming matrix. Structural overlap is used both as a measure of similarity and to assess alignment quality. Pairwise alignment accuracy has been benchmarked against that of 10 widely used aligners on the Sippl and Wiederstein set of difficult pairwise structure alignment problems, and more extensively against that of Matt, SALIGN, and MUSTANG in pairwise and multiple structural alignments of protein domains with low shared sequence identity in the SCOP-ASTRAL 40% compendium. The results demonstrate the advantages of POLYFIT over other aligners in the efficient and robust identification of matching seed residue positions in distantly related protein targets and in the generation of longer structurally overlapped alignment lengths. Superposition-based application areas include comparative modeling and protein and ligand design. POLYFIT is available on the Web server at http://polyfit.insa-toulouse.fr. Proteins 2013; 81:1823-1839. (c) 2013 Wiley Periodicals, Inc.

Nociceptin and the ORL1 receptor share high sequence similarity with opioid peptides, particularly dynorphin A, and their receptors. However, nociceptin and dynorphin A may use distinct molecular pathways to bind and activate their cognate receptors. Activation of the kappa-opioid receptor by dynorphin A is thought to require interactions of its N-terminal hydrophobic domain (Y(1)GGF) with the receptor opioid binding pocket, located within the transmembrane helix bundle, while activation of the ORL1 receptor appears to require interactions of the positively charged core (R(8)KSARK) of nociceptin with the negatively charged second extracellular receptor loop.

The theoretical foundation has been laid for the investigation of catalytic systems using first-order kinetics and for a general kinetic method of investigation of the active site content, E(a), of enzymes, catalytic antibodies, and other enzyme-like catalysts. The method involves a combination of steady-state and single-turnover kinetics to provide Vmax and Km and k(lim)(obs) and K(app)(m), respectively. The validity of the method is shown to remain valid for two extensions of the simple two-step enzyme catalysis model (a) when the catalyst preparation contains molecules (Eb) that bind substrate but fail to catalyse product formation and (b) when the catalyst itself binds substrate non-productively as well as productively. The former is a particularly serious complication for polyclonal catalytic antibodies and the latter a potential complication for all catalysts. For the simple model and for (b) Vmax/k(lim)(obs) provides the value of [Ea]T and for (a) its upper limit. This can be refined by consideration of the relative values of Km and the equilibrium dissociation constant of EbS. For the polyclonal catalytic antibody preparation investigated, the fact that K(app/m)Km demonstrates for the first time the presence of a substrate-binding but non-catalytic component in a polyclonal preparation. First-order behaviour in catalytic systems occurs not only with a large excess of catalyst over substrate but also with lower catalyst/substrate ratios, including the equimolar condition, when K(app)(m)[S]0, a phenomenon that is not widely appreciated.

A systematic analysis of peptide nucleic acid (PNA) complexes deposited in the Protein Data Bank has been carried out using a set of contiguous atom torsion angle definitions. The analysis is complemented by molecular mechanics adiabatic potential energy calculations on hybrid PNA-nucleic acid model systems. Hitherto unobserved correlations in the values of the (α and ϵ) dihedral angles flanking the backbone secondary amide bond are found. This dihedral coupling forms the basis of a PNA backbone conformation classification scheme. Six conformations are thus characterised in experimental structures. Helix morphology is found to exert a significant influence on backbone conformation and flexibility: Watson-Crick PNA strands in complexes with DNA and RNA, that possess A -like base-pair stacking, adopt backbone conformations distinct from those in PNA·DNA-PNA triplex and PNA-PNA duplex P-helix forms. Solvation effects on Watson-Crick PNA backbone conformation in heterotriplexes are discussed and the possible involvement of inter-conformational transitions and dihedral angle uncoupling in asymmetric heteroduplex base-pair breathing is suggested.

Les agents therapeutiques qui empechent la replication et la transcription en se liant a l'ADN double-brin nous donnent la possibilite de combattre le cancer en minimisant les effets secondaires. Nous decrivons ici la parametrisation d'une fonction d'energie pour des molecules qui se lient au sillon mineur de l'ADN. Les parametres etaient ajustes pour reproduire les resultats de calculs de chimie quantique ab initio. Le champs de force derive devrait etre utile dans la comprehension d'aspects fondamentaux dans la conceptualisation de medicaments qui ciblent l'ADN.

The POLYFIT rigid-body algorithm for automated global pairwise and multiple protein structural alignment is presented. Smith-Waterman local alignment is used to establish a set of seed equivalences that are extended using Needleman-Wunsch dynamic programming techniques. Structural and functional interaction constraints provided by evolution are encoded as one-dimensional residue physical environment strings for alignment of highly structurally overlapped protein pairs. Local structure alignment of more distantly related pairs is carried out using rigid-body conformational matching of 15-residue fragments, with allowance made for less stringent conformational matching of metal-ion and small molecule ligand-contact, disulphide bridge, and cis-peptide correspondences. Protein structural plasticity is accommodated through the stepped adjustment of a single empirical distance parameter value in the calculation of the Smith-Waterman dynamic programming matrix. Structural overlap is used both as a measure of similarity and to assess alignment quality. Pairwise alignment accuracy has been benchmarked against that of 10 widely used aligners on the Sippl and Wiederstein set of difficult pairwise structure alignment problems, and more extensively against that of Matt, SALIGN, and MUSTANG in pairwise and multiple structural alignments of protein domains with low shared sequence identity in the SCOP-ASTRAL 40% compendium. The results demonstrate the advantages of POLYFIT over other aligners in the efficient and robust identification of matching seed residue positions in distantly related protein targets and in the generation of longer structurally overlapped alignment lengths. Superposition-based application areas include comparative modeling and protein and ligand design. POLYFIT is available on the Web server at http://polyfit.insa-toulouse.fr.

International audience; To metabolize both dietary fiber constituent carbohydrates and host glycans lining the intestinal epithelium, gut bacteria produce a wide range of carbohydrate-active enzymes, of which glycoside hydrolases are the main components. In this study, we describe the ability of phosphorylases to participate in the breakdown of human N-glycans, from an analysis of the substrate specificity of UhgbMP, a mannoside phosphorylase of the GH130 protein family discovered by functional metagenomics. UhgbMP is found to phosphorolyze -d-Manp-1,4--d-GlcpNAc-1,4-d-GlcpNAc and is also a highly efficient enzyme to catalyze the synthesis of this precious N-glycan core oligosaccharide by reverse phosphorolysis. Analysis of sequence conservation within family GH130, mapped on a three-dimensional model of UhgbMP and supported by site-directed mutagenesis results, revealed two GH130 subfamilies and allowed the identification of key residues responsible for catalysis and substrate specificity. The analysis of the genomic context of 65 known GH130 sequences belonging to human gut bacteria indicates that the enzymes of the GH130_1 subfamily would be involved in mannan catabolism, whereas the enzymes belonging to the GH130_2 subfamily would rather work in synergy with glycoside hydrolases of the GH92 and GH18 families in the breakdown of N-glycans. The use of GH130 inhibitors as therapeutic agents or functional foods could thus be considered as an innovative strategy to inhibit N-glycan degradation, with the ultimate goal of protecting, or restoring, the epithelial barrier.

The amylosucrase from Neisseria polysaccharea is a transglucosidase from the GH13 family of glycoside-hydrolases that naturally catalyzes the synthesis of -glucans from the widely available donor sucrose. Interestingly, natural molecular evolution has modeled a dense hydrogen bond network at subsite -1 responsible for the specific recognition of sucrose and conversely, it has loosened interactions at the subsite +1 creating a highly promiscuous subsite +1. The residues forming these subsites are considered to be likely involved in the activity as well as the overall stability of the enzyme. To assess their role, a structure-based approach was followed to reshape the subsite -1. A strategy based on stability change predictions, using the FoldX algorithm, was considered to identify the best candidates for site-directed mutagenesis and guide the construction of a small targeted library. A miniaturized purification protocol was developed and both mutant stability and substrate promiscuity were explored. A range of 8 degrees C between extreme melting temperature values was observed and some variants were able to synthesize series of oligosaccharides with distributions differing from that of the parental enzyme. The crucial role of subsite -1 was thus highlighted and the biocatalysts generated can now be considered as starting points for further engineering purposes.

Chymopapain M, the monothiol cysteine proteinase component of the chymopapain band eluted after chymopapains A and B in cation-exchange chromatography, was isolated from the dried latex of Carica papaya and characterized by kinetic and chromatographic analysis. This late-eluted chymopapain is probably a component of the cysteine proteinase fraction of papaya latex discovered by Schack [(1967) Compt. Rend. Trav. Lab. Carlsberg 36, 67-83], named papaya peptidase B by Lynn [(1979) Biochim. Biophys. Acta 569, 193-201] and partially characterized by Polgár [(1981) Biochim. Biophys. Acta 658, 262-269] and is the enzyme with unusual specificity characteristics (papaya proteinase IV) that Buttle, Kembhavi, Sharp, Shute, Rich and Barrett [Biochem. J. (1989) 261, 469-476] claimed to be a previously undetected cysteine proteinase eluted from a cation-exchange column near to the early-eluted chymopapains. A study of the time-dependent chromatographic consequences of thiol-dependent proteolysis of the components of papaya latex is reported. Chymopapain M was isolated by (i) affinity chromatography followed by separation from papain using cation-exchange f.p.l.c. on a Mono S HR5/5 column and (ii) cation-exchange chromatography followed by an unusual variant of covalent chromatography by thiol-disulphide interchange. The existence in chymopapain M of a nucleophilic interactive Cys/His catalytic-site system analogous to those in papain (EC 3.4.22.2) and other cysteine proteinases was deduced from the characteristics shape of the pH-second-order rate constant (k) profiles for its reactions with 2,2′-dipyridyl disulphide and ethyl 2-pyridyl disulphide. Analysis of the pH-k data for the reactions of chymopapain M with the 2-pyridyl disulphides and with 4,4′-dipyridyl disulphide permits the assignment of molecular pKa values of 3.4 and 8.7 to the formation and subsequent dehydronation of the Cys-S-/His-Im+H state of the catalytic site and reveals three other kinetically influential ionizations with pKa values 3.4, 4.3 and 5.6. The pH-dependences of kcat./Km for the hydrolysis of N-acetyl-L-Phe-Gly-4-nitroanilide at 25.0 degrees C and I0.1 M catalysed by chymopapain M and papain were determined. For both enzymes, little catalytic activity (5-7% of the maximal) develops consequent on formation of the catalytic site Cys-S-/His-Im+H ion-pair state (across pKa 3.4 for both enzymes). For papain, full expression of Kcat./Km for the uncharged substrate requires only the additional hydronic dissociation with pKa 4.2. By contrast, full expression of kcat./Km for chymopapain M requires additional hydronic dissociation with pKa values of 4.3 and 5.6.(ABSTRACT TRUNCATED AT 400 WORDS)

A three-dimensional model of human renin has been constructed using COMPOSER, a rule-based approach to comparative modelling. The model was constructed from the three-dimensional structures of two homologous aspartic proteinases, pepsin and chymosin, which have been defined by X-ray analysis. The variable regions were obtained from aspartic proteinase structures for most of the loops, but also from unrelated structures where appropriate conformers were not available within the family. Using a 3.5A cut-off for pairs of equivalent atoms, 280 C-(Y atom airs of the model and the experimentally defined structure superpose with an rmsd of 0.841. Larger deviations occur in some of the loop regions, especially around residues 290 where there is a proline-rich sequence with CIS-peptides. The analysis shows that comparative modelling can give a reasonably accurate estimate of homologous protein structures especially in the active-site region and that these provide a useful basis for structure-based design in the absence of experimental data.

The family of aspartic proteinases includes several human enzymes that may play roles in both physiological and pathophysiological processes. The human lysosomal aspartic proteinase cathepsin D is thought to function in the normal degradation of intracellular and endocytosed proteins but has also emerged as a prognostic indicator of breast tumor invasiveness. Presented here are results from a continuing effort to elucidate the factors that contribute to specificity of ligand binding at individual subsites within the cathepsin D active site. The synthetic peptide Lys-Pro-Ile-Glu-Phe*Nph-Arg-Leu has proven to be an excellent chromogenic substrate for cathepsin D yielding a value of kcat/Km = 0.92 x 10(-6) s-1 M-1 for enzyme isolated from human placenta. In contrast, the peptide Lys-Pro-Ala-Lys-Phe*Nph-Arg-Leu and all derivatives with Ala-Lys in the P3-P2 positions are either not cleaved at all or cleaved with extremely poor efficiency. To explore the binding requirements of the S3 and S2 subsites of cathepsin D, a series of synthetic peptides was prepared with systematic replacements at the P2 position fixing either Ile or Ala in P3. Kinetic parameters were determined using both human placenta cathepsin D and recombinant human fibroblast cathepsin D expressed in Escherichia coli. A rule-based structural model of human cathepsin D, constructed on the basis of known three-dimensional structures of other aspartic proteinases, was utilized in an effort to rationalize the observed substrate selectivity.

Conformationally constrained environment-dependent amino acid residue substitution tables have been constructed from a database comprising 33 homologous families of protein sequences aligned on the basis of their three-dimensional structures. Residues are allotted to one of 216 (or 54) classes of combinations of structural features. These include nine main-chain conformation classes, three classes of side-chain accessibility and eight (or two) classes of side-chain involvement in three types of hydrogen bond. Seven different main-chain conformational classes outside of regions of regular structure were identified in an analysis of the distributions of φ-ψ torsion angles in 84 high-resolution crystallographic structures. Residue substitutions at equivalent positions in the structural alignments are included where the main-chain conformational class is conserved. Frequency data in the form of 216 (or 54) environment specific (20 × 20 residue type) matrices are then converted to probabilities. Two smoothing regimes incorporating entropy-driven weights were applied to the set of 54 tables. Predicted residue substitutions have been generated for individual residue positions in β-hairpins and the hypervariable regions of the immunoglobulins. These have been compared with the observed sequence variation at the same positions using rank correlation methods. Measurements of χ2 distances demonstrate the considerable improvement in predictive power at key residue positions identified from interactive graphics studies when compared to the Dayhoff MDM250 mutation matrix. An illustrative example is given of an application of the method in the ranking of loop fragments in model building studies of structurally variable regions in two subtilisins. A combined template scoring procedure is found to be 26-fold more discriminatory than the Dayhoff matrix. The success rate is approximately 85%.

1. 4-Nitrophenyl 4′-(3-aza-2-oxoheptyl)phenyl carbonate (I), an amide conjugate (XI) involving the carboxy group of 4-nitrophenyl 4′-carboxymethylphenyl phosphate and an amino group of keyhole-limpet haemocyanin, and a fluorescein derivative (XVII) were synthesized. 2. The conjugate (XI) was used as an immunogen with which to raise polyclonal antibodies in multigeneration cross-bred sheep; the fluorescent derivative (XVII) was used for the initial assessment of the antisera via binding assays monitored by fluorescence polarization; the carbonate ester (I) was used as a chromogenic substrate for the investigation of catalytic activity. 3. The IgG from the antiserum of sheep no. 270 was isolated by Na2SO4 precipitation and chromatography on Protein G-Sepharose. 4. This preparation of IgG catalysed the hydrolysis of the carbonate ester (I); the catalysis at pH 8.0 and 25 degrees C obeyed Michaelis-Menten kinetics with at least 25 turnovers, Km = 3.34 microM, and lower limits for kcat. of 0.029 s-1 and for kcat./Km of 8.77 x 10(3) M-1.S-1, on the unlikely assumption that the concentration of catalytic antibody is provided by twice the total IgG concentration (two sites per molecule); probable estimates of the fraction of the total IgG that is anti-haptenic IgG and of the fraction of this that is catalytically active suggest that the values of kcat./Km are actually very much larger than these lower limits. 5. The failure of the antibody preparation to catalyse the hydrolysis of the isomeric 2-nitrophenyl carbonate (II), which differs from compound (I) only in the position of the nitro substituent in the leaving group, compels the view that catalytic activity is due to antibody rather than contaminant enzyme; this conclusion is supported by (a) the failure of the following to discriminate effectively between the isomeric substrates (I) and (II): pig liver carboxylesterase, rabbit liver carboxylesterase (collectively EC 3.1.1.1), whole serum from a non-immunized sheep and whole serum from a sheep immunized with a derivative of 3-O-methylnoradrenaline and (b) the lack of catalytic activity in IgG preparations from sheep immunized with sulphoxide or sulphone analogues of immunogen (XI). 6. The various parameters used for the comparison of the kinetic characteristics of hydrolytic catalytic antibodies are discussed. 7. The characteristics of hydrolysis of compound (I) catalysed by the present polyclonal antibody preparation are shown to be substantially better in most respects than those of analogous reactions of two other carbonate esters catalysed by monoclonal antibodies.

A generalized theoretical treatment of the kinetics of an enzyme-catalysed reaction in the presence of an unstable irreversible inhibitor (or activator) is presented. Analytical expressions describing the time-dependence of product formation have been derived in coefficient form amenable to non-linear regression analysis for two operationally distinct types of reaction mechanism dependent on whether the reaction of the unstable modifier (X) with either or both the free enzyme (E) and enzyme-substrate complex (ES) occurs as a simple bimolecular process, or proceeds through the intermediacy of either or both adsorptive enzyme-modifier (EX) and enzyme-modifier-substrate (EXS) complexes in what may be considered as an extension of the Botts-Morales general modifier mechanism for (stable) reversible enzyme inhibitors and activators. Special cases of both models are classified in an analogous way to the traditional naming of reversible enzyme modifications, and guidelines concerning tests of mechanism and determination of kinetic parameters are given. In particular, it has been shown that kinetic constants describing enzyme inactivation by an unstable site-specific inhibitor forming a reversible EX complex prior to covalent modification step may be determined from a single progress curve. Kinetic analysis of the extended Botts-Morales mechanism describing irreversible enzyme inactivation has demonstrated that analytical expressions describing the time-course of product formation may be derived for a stable modifier by retaining the usual steady-state assumptions regarding the fluxes around ES and EXS provided quasi-equilibrium modifier binding to E and ES is assumed, but for unstable modifiers all of the binding steps must be assumed to be at quasi-equilibrium in the steady-state, except under restrictive circumstances.

Kinetic characterization of enzymes and analogous catalysts such as catalytic antibodies requires knowledge of the molarity of functional sites. Various stoichiometric titration methods are available for the determination of active-site concentrations of some enzymes and these are exemplified in the second part of this article. Most of these are not general in that they require the existence of certain types of either intermediate or active-site residues that are susceptible to specific covalent modification. Thus they are not readily applicable to many enzymes and they are rarely available currently for titration of catalytic antibody active sites. In the first part of the article we discuss a general kinetic method for the investigation of active-site availability in preparations of macromolecular catalysts. The method involves steady-state kinetics to provide Vmax and Km and single-turnover first-order kinetics using excess of catalyst over substrate to provide the analogous parameters k(obs)lim and K(m)app. The active-site contents of preparations that contain only active catalyst (Ea) and inert material (Ei) may be calculated as [Ea](T) = Vmax)/k(obs)lim. This is true even if nonproductive binding to E(a) occurs. For polyclonal catalytic antibody preparations, which may contain binding but noncatalytic material (Eb) in addition to Ea and Ei, the significance of Vmax/k(obs)lim is more complex but provides an upper limit to E(a). This can be refined by consideration of the relative values of Km and the equilibrium dissociation constant of EbS. Analysis of the Ea, Eb, Ei system requires the separate determination of Ei. For catalytic antibodies this may be achieved by analytical affinity chromatography using an immobilized hapten or hapten analog and an ELISA procedure to ensure the clean separation of Ei from the Ea + Eb mixture.

Evidence from statistical cluster analyses of a multiple sequence alignment of G protein-coupled receptor seven-helix folds supports the existence of structurally conserved transmembrane (TM) ligand binding sites in the opioid/opioid receptor-like (ORL1) and amine receptor families. Based on the expectation that functionally conserved regions in homologous proteins will display locally higher levels of sequence identity compared with global sequence similarities that pertain to the overall fold, this approach may have wider applications in functional genomics to annotate sequence data. Binding sites in models of the kappa-opioid receptor seven-helix bundle built from the rhodopsin templates of Baldwin et al. (1997) [J. Mol. Biol., 272, 144-164] and Herzyk and Hubbard (1998) [J. Mol. Biol., 281, 742-751] are compared. The Herzyk and Hubbard template is found to be in better accord with experimental studies of amine, opioid and rhodopsin receptors owing to the reduced physical separation of the extracellular parts of TM helices V and VI and differences in the rotational orientation of the N-terminal of helix V that reveal side chain accessibilities in the Baldwin et al. structure to be out of phase with relative alkylation rates of engineered cysteine residues in the TM binding site of the alpha(2A)-adrenergic receptor. TM helix V in the Baldwin et al. template has been remodelled with a different proline kink to satisfy experimental constraints. A recent proposal that rotation of helix V is associated with receptor activation is critically discussed.

The heptadecapeptide nociceptin, also known as orphanin FQ, is the endogenous agonist of the opioid receptor-like 1 (ORL1) G protein-coupled receptor. An affinity labeling approach has been implemented to probe the interactions of the neuropeptide with the receptor using the photolabile nociceptin derivative, [p-benzoyl-l-Phe(10),Tyr(14)]nociceptin ([Bpa(10),Tyr(14)]noc). In recombinant Chinese hamster ovary cells expressing the human ORL1 receptor, [Bpa(10),Tyr(14)]noc binds the receptor with high affinity (K(i) approximately 0.7 nm) and is as potent as nociceptin in the inhibition of forskolin-induced cAMP synthesis (EC(50) approximately 0.5 nm). UV irradiation at 365 nm of the complex formed by the ORL1 receptor and radioiodinated [Bpa(10),Tyr(14)]noc results in the irreversible labeling of a glycoprotein of approximately 65 kDa, determined by SDS-polyacrylamide gel electrophoresis. Complete digestion of the partially purified 65-kDa complex with kallikrein generates a single labeled fragment (approximately 6.5 kDa) that is readily cleaved by endoproteinase Glu-C to yield a labeled fragment of approximately 3.2 kDa. Kallikrein treatment of the photoaffinity cross-linked Glu(295) --Asp mutant receptor also yields a single labeled fragment of approximately 6.5 kDa but is resistant to further cleavage by endoproteinase Glu-C. Based upon the expected proteolytic fingerprint of the labeled receptor, the photoreactive region can be identified as ORL1-(296-302; residues Thr-Ala-Val-Ala-Ile-Leu-Arg) spanning the C terminus of extracellular loop 3 and the N terminus of transmembrane helix VII. Molecular modeling of the ORL1 receptor complex with [Bpa(10)]noc suggests that reaction of the Bpa carbonyl group may occur with the side chain of Ile(300) within the experimentally identified photoreactive region.

A potentially general kinetic method for the investigation of active-site availability in preparations of macromolecular catalysts was developed. Three kinetic models were considered: (a) the conventional two-step model of enzyme catalysis, where the preparation contains only active catalyst (E(a)) and inert (i.e. non-binding, non-catalytic) material (E(i)); (b) an extension of the conventional model (a) involving only E(a) and E(i), but with non-productive binding to E(a) (in addition to productive binding); (c) a model in which the preparation contains also binding but non-catalytic material (E(b)), predicted to be present in polyclonal catalytic antibody preparations. The method involves comparing the parameters V(max) and K(m) obtained under catalytic conditions where substrate concentrations greatly exceed catalyst concentration with those (klim/obs, the limiting value of the first-order rate constant, k(obs), at saturating concentrations of catalyst; and Kapp/m) for single-turnover kinetics, in which the reverse situation obtains. The active-site contents of systems that adhere to model (a) or extensions that also lack E(b), such as the non-productive binding model (b), may be calculated using [E(a)](T)=V(max)/klim/obs. This was validated by showing that, for alpha-chymotrypsin, identical values of [E(a)](T) were obtained by the kinetic method using Suc-Ala-Ala-Pro-Phe-4-nitroanilide as substrate and the well-known 'all-or-none' spectroscopic assay using N-trans-cinnamoylimidazole as titrant. For systems that contain E(b), such as polyclonal catalytic antibody preparations, V(max)/klim/obs is more complex, but provides an upper limit to [E(a)](T). Use of the kinetic method to investigate PCA 271-22, a polyclonal catalytic antibody preparation obtained from the antiserum of sheep 271 in week 22 of the immunization protocol, established that [E(a)](T) is less than approx. 8% of [IgG], and probably less than approx. 1% of [IgG].

An approach to the prediction of mutant stability is described using knowledge of amino acid replacements that are tolerated within the families of homologous proteins of known 3-D structure. Amino acid variations in families of homologous proteins are converted to propensity and substitution tables; these provide quantitative information about the existence of an amino acid in a structural environment and the probability of replacement by any other amino acid. The tables are used to calculate a 'stability difference score', analogous to the difference in free energy between a mutant and the wild type. The method has been developed and tested using the high-resolution structures for T4 lysozyme and 159 site-specific mutants. We show that differences in stability scores are correlated with experimentally observed free energy differences and differences in melting temperature. Blind tests, using only structural information derived from the parent wild-type crystal structures, on a combined set of 83 staphylococcal nuclease and 68 barnase mutants showed a correlation of 0.80 in the predicted stability changes with experimental thermodynamic data. Approximately 86% of the predictions were correctly classified as destabilizing or stabilizing.

We have used two structurally well-characterized serine proteinase variants, subtilisins Carlsberg and BPN', to produce (Cys)-S-/(His)-Im+H ion-pairs by chemical mutation in well defined, different, electrostatic microenvironments. These ion-pairs have been characterized by pH-dependent rapid reaction kinetics using, as reactivity probes, thiol-specific time dependent inhibitors, 2,2'-dipyridyl disulfide and 4,4'-dipyrimidyl disulfide, that differ in the protonation states of their leaving groups in acidic media, computer modelling and electrostatic potential calculations. Both ion-pairs possess nucleophilic character, identified by the striking rate maxima in their reactions with 2,2'-dipyridyl disulfide in acid media. In the Carlsberg enzyme, the (Cys220)-S-/(His63)-Im+H ion-pair is produced by protonic dissociation with pKa 4.1 and its reactivity is not perturbed by any detectable electrostatic influence other than the deprotonation of His63 (pKa 10.2). In the BPN' enzyme, the analogous, (Cys221)-S-/(His64)-Im+H ion-pair is produced by protonic dissociation with pKa 5.1 and its reactivity is affected by an ionization with pKa 3.5 in addition to the deprotonation of His64 (pKa > or = 10.35). It is a striking result that calculations using finite difference solutions of the Poisson-Boltzmann equation provide a value of the pKa difference between the two enzyme catalytic sites (0.97) in close agreement with the value (1.0) determined by reactivity probe kinetics when a protein dielectric constant of 2 is assumed and water molecules within 5 A of the catalytic site His residue are included. The pKa difference is calculated to be 0.84 when the water molecules are not included and a protein dielectric constant of 20 is assumed. The calculations also identify Glu156 in the BPN' enzyme (which is Ser in the Carlsberg enzyme) as the main individual source of the pKa shift. The additional kinetically influential pKa of 3.5 is assigned to Glu156 by examining the non-covalent interactions between the 2-pyridyl disulfide reactivity probe and the enzyme active centre region.

A model of the 3-D structure of a major house dust mite allergen Der p I associated with hypersensitivity reactions in humans was built from its amino acid sequence and its homology to three known structures, papain, actinidin and papaya proteinase omega of the cysteine proteinase family. Comparative modelling using COMPOSER was used to arrive at an initial model. This was refined using interactive graphics and energy minimization with the AMBER force field incorporated in SYBYL (Tripos Associates). Compatibility of the Der p I amino acid sequence with the cysteine proteinase fold was checked using an environment-dependent amino acid propensity table incorporated into a new program HARMONY with a variable length windowing facility. A five-residue window was used to probe local conformational integrity. Propensities were derived from a structural alignment database of homologous proteins using a robust entropy-driven smoothing procedure. Der p I shares essential structural and mechanistic features with other papain-like cysteine proteinases, including cathepsin B. The active-site thiolate-imidazolium ion pair comprises the side chains of Cys34 and His170. A cystine disulfide not present in other known structures bridges residue 4 of an N-terminal extension and the core residue 117. Two conserved disulfide bridges are formed by residues 31 and 71 and residues 65 and 103. Model building of peptide substrate analogue complexes suggests a preference for phenylalanyl or basic residues at the P2 position, whilst selectivity may be of minor importance at the S1 subsite. The electrostatic influences on the Der p I active-site ion pair and extended peptide binding region are markedly different from those in known structures. A highly immunogenic surface exposed region (residues 107-131), comprising several overlapping T cell epitope sites, has no shared sequence identity with human liver cathepsin B and contains three insertion-deletion sites. The structure provides a basis for testing the substrate specificity of Der p I and the potential role of proteinase activity in hypersensitivity reactions. These studies may offer a new treatment strategy by hyposensitization with inactive mutants or mutants with significantly altered proteinase activity, either alone or complexed with antibody.

1. A new thiol-specific reactivity probe 4,4′-dipyrimidyl disulphide [compound (VII), m.p. 110 degrees C, pKa of its monohydronated form 0.91] was synthesized and used to resolve the ambiguity of interpretation of the behaviour of papain (EC 3.4.22.2) in alkaline media known to depend to varying extents on two ionizations with pKa values approx. 8.0-8.5 and > or = 9.5 respectively. 2. A new extensive pH-second-order rate constant (k) data set for the reaction of papain with 2-(acetamido)-ethyl 2′-pyridyl disulphide (IV) demonstrated the existence of a striking rate maximum at pH approx. 4, the independence of k around pH 8 and the increase in k with increase in pH across a pKa value of 10.0, behaviour similar to that of other 2-pyridyl disulphides (R-S-S-2-Py) that lack key substrate-like binding sites in R. 3. Although the simplest interpretation of the pKa value of 10.0 assigns it to the formation of (Cys-25)-S-/(His-159)-Im from the ion-pair state of the papain catalytic site, another interpretation may be conceived in which this pKa value is assigned to another group remote from the catalytic site, the state of ionization of which modulates catalytic-site behaviour. This alternative assignment is shown to require compensating effects in the pH region around 8 such that the formation of (Cys-25)-S-/(His-159)-Im across pKa 8.0-8.5 is without net kinetic effect in the reactions of simple 2-pyridyl disulphides such as compound (IV) and 2,2′-dipyridyl disulphide (II). 4. The lower basicity of compound (VII) relative to that of compound (II) (pKa 2.45) was predicted to diminish or abolish the compensation postulated as a possibility in reactions of 2-pyridyl disulphides because of the decreased effectiveness of reaction via a (His-159)-Im+H-assisted transition state. The characteristics of the pH-dependence of the reaction of papain with compound (VII) which are quite different from those for its reaction with compound (II) support both this prediction and the alternative assignment with a value of 8.3 for the pKa of the formation of (Cys-25)-S-/(His-159)-Im. 5. Evidence that the behaviour of papain towards both substrates and some substrate-derived time-dependent inhibitors is determined not only by the loss of the (Cys-25)-S-/(His-159)-Im+H ion-pair state by dehydronation with pKa 8.3 but also by another ionization of pKa approx. 10.0 is briefly discussed.

1. A model of the three-dimensional structure of papaya proteinase omega, the most basic cysteine proteinase component of the latex of papaya (Carica papaya), was built from its amino acid sequence and the two currently known high-resolution crystal structures of the homologous enzymes papain (EC 3.4.22.2) and actinidin (EC 3.4.22.14). The method used a knowledge-based approach incorporated in the COMPOSER suite of programs and refinement by using the interactive graphics program FRODO on an Evans and Sutherland PS 390 and by energy minimization using the GROMOS program library. 2. Functional similarities and differences between the three cysteine proteinases revealed by analysis of pH-dependent kinetics of the acylation process of the catalytic act and of the reactions of the enzyme catalytic sites with substrate-derived 2-pyridyl disulphides as two-hydronic-state reactivity probes are reported and discussed in terms of the knowledge-based model. 3. To facilitate analysis of complex pH-dependent kinetic data, a multitasking application program (SKETCHER) for parameter estimation by interactive manipulation of calculated curves and a simple method of writing down pH-dependent kinetic equations for reactions involving any number of reactive hydronic states by using information matrices were developed. 4. Papaya proteinase omega differs from the other two enzymes in the ionization characteristics of the common (Cys)-SH/(His)-Im+H catalytic-site system and of the other acid/base groups that modulate thiol reactivity towards substrate-derived inhibitors and the acylation process of the catalytic act. The most marked difference in the Cys/His system is that the pKa for the loss of the ion-pair state to form -S-/-Im is 8.1-8.3 for papaya proteinase omega, whereas it is 9.5 for both actinidin and papain. Papaya proteinase omega is similar to actinidin in that it lacks the second catalytically influential group with pKa approx. 4 present in papain and possesses a catalytically influential group with pKa 5.5-6.0. 5. Papaya proteinase omega occupies an intermediate position between actinidin and papain in the sensitivity with which hydrophobic interaction in the S2 subsite is transmitted to produce changes in transition-state geometry in the catalytic site, a fact that may be linked with differences in specificity in P2-S2 interaction exhibited by the three enzymes.(ABSTRACT TRUNCATED AT 400 WORDS)

Ct iRIS7'OPHER M. TOPHAM.* JOHN OVERINGTON.* MAIREAD O'DRISCOLL,t# ERDJAN SALIH,tII MARK THOMAS.tll EMRYS W. THOMAS$ and KEITH BROCKLEH URSTt 11 * Liiborii ton of'A4olecirlcrr Biology, Ilep(irttnetii of ~ ' t ~ . s t i i l l o g r u ~ ~ i i y , Hirk beck College, U t i i vmi t y of Lotidoti, M i i l c ~ Strcw, Lotidoti W C I E 7/IX, U. K . , t Depirttnetit of Hioc~hctiiistty S I B i v t l i o h i e w :v ilospiiul Mediciil C 'ollege, Utiiiwsity of' Lotidoti, Cliurterhoirse Sqiriire. Lotidoti EC'IM OHQ, U. K . utid $ Dtyurttmwt ofHiologicu1 Scierices. Utii iwvity of Siilford. Srilfortl A45 4 W I ; U. K .

1. 2-(N'-Acetyl-D-phenylalanylamino)ethyl 2′-pyridyl disulphide (compound I) [m.p. 123-124 degrees C; [alpha]20D -7.1 degrees (c 0.042 in methanol)] was synthesized, and the results of a study of the pH-dependence of the second-order rate constant (k) for its reaction with the catalytic-site thiol group of papain (EC 3.4.22.2), together with existing kinetic data for the analogous reaction of the L-enantiomer (compound II), were used to evaluate the consequences for transition-state geometry of the difference in chirality at the P2 position of the probe molecule. 2. The kinetic data suggest that the D-enantiomer binds approx. 40-fold less tightly to papain than the L-enantiomer but that the binding-site-catalytic-site signalling that results in a (His-159)-Im(+)-H-assisted transition state occurs equally effectively in the interaction of the former probe as in that of the latter. This results in pH-k profiles for the reactions of both enantiomers each characterized by four macroscopic pKa values (3.7-3.9, 4.1-4.3, 7.9-8.3 and 9.4-9.5) in which k is maximal at pH approx. 6 where the -Im(+)-H-assisted transition state is most fully developed. 3. Model building indicates that both enantiomers can bind to papain such that the phenyl ring of the N-acetylphenylalanyl group makes hydrophobic contacts in the binding pocket of the S2 subsite with preservation of the three hydrogen-bonding interactions involving the substrate analogue reagent and (Asp-158) C = O, (Gly-66) C = O, and (Gly-66)-N-H of papain. Earlier predictions that binding of N-acyl-D-phenylalanine derivatives to papain would be prevented on steric grounds [Berger & Schechter (1970) Philos. Trans. R. Soc. London B 257, 249-264; Lowe & Yuthavong (1971) Biochem. J. 124, 107-115; Lowe (1976) Tetrahedron 32, 291-302] were based on assumed models that are not consistent with the X-ray-diffraction data for papain inhibited by alkylation of Cys-25 with N-benzyloxycarbonyl-Phe-Ala-chloromethane [Drenth, Kalk & Swen (1976) Biochemistry 15, 3731-3738]. 4. The possibility that the kinetic expression of P2-S2 stereospecificity may depend on the nature of the chemistry occurring in the catalytic site of papain is discussed.(ABSTRACT TRUNCATED AT 400 WORDS)

1. N-Acetyl-L-phenylalanylglycine 4-nitroanilide and its D-enantiomer were synthesized and characterized and used as substrates with which to evaluate stereochemical selectivity in papain (EC 3.4.22.2)-catalysed hydrolysis. 2. Kinetic analysis at pH 6.0, I 0.1, 8.3% (v/v) NN-dimethylformamide and 25 degrees C by using initial-rate data with [S] much less than Km and weighted non-linear regression provided values of kcat./Km for the catalysed hydrolysis of both enantiomers as (kcat./Km)L = 2040 +/- 48 M-1.S-1 and (kcat./Km)D = 5.9 +/- 0.07 M-1.S-1. These data, taken together with individual values of kcat. and Km for the hydrolysis of the L-enantiomer (a) estimated in the present work as kcat. = 3.2 +/- 1.2 S-1 and Km = 1.5 +/- 0.6 mM and (b) reported by Lowe & Yuthavong [(1971) Biochem. J. 124, 107-115] for the reaction at pH 6.0 in 10% (v/v) NN-dimethylformamide and 35 degrees C, as kcat. = 1.3 +/- 0.2 S-1 and Km = 0.88 +/- 0.1 mM, suggest that (kcat./Km)L congruent to 2000 M-1.S-1 and thus that (kcat./Km)L/(kcat./Km)D congruent to 330.3. Model building indicates that both enantiomeric 4-nitroanilides can bind to papain such that the phenyl ring of the N-acetylphenylalanyl group makes hydrophobic contacts in the S2 subsite with preservation of mechanistically relevant hydrogen-bonding interactions and that the main difference is in the positioning of the beta-methylene group. 4. The dependence of P2-S2 stereochemical selectivity of papain on the nature of the catalytic-site chemistry for reactions involving derivatives of N-acetylphenylalanine is discussed. The variation in the index of stereochemical selectivity (ratio of the appropriate kinetic or thermodynamic parameter for a given pair of enantiomeric ligands), from 330 for the overall acylation process of the catalytic act, through 40 and 31 for the reaction at electrophilic sulphur in 2-pyridyl disulphides respectively without and with assistance by (His-159)-Im(+)-H, to 5 for the formation of thiohemiacetal adducts by reaction at aldehydic carbon, is interpreted in terms of the extent to which conformational variation of the bound ligand in the catalytic-site region permits the binding mode of the -CH2-Ph group of the D-enantiomer to approach that of the L-enantiomer.

1. The pH-dependence of the second-order rate constant (k) for the reaction of actinidin (EC 3.4.22.14) with 2-(N'-acetyl-L-phenylalanylamino)ethyl 2'-pyridyl disulphide was determined and the contributions to k of various hydronic states were evaluated. 2. The data were used to assess the consequences for transition-state geometry of providing P2/S2 hydrophobic contacts in addition to hydrogen-bonding opportunities in the S1-S2 intersubsite region. 3. The P2/S2 contacts (a) substantially improve enzyme-ligand binding, (b) greatly enhance the contribution to reactivity of the hydronic state bounded by pKa 3 (the pKa characteristic of the formation of catalytic-site-S-/-ImH+ state) and pKa 5 (a relatively minor contributor in reactions that lack the P2/S2 contacts), such that the major rate optimum occurs at pH 4 instead of at pH 2.8-2.9, and (c) reveal the kinetic influence of a pKa approx. 6.3 not hitherto observed in reactions of actinidin. 4. Possibilities for the interplay of electrostatic effects and binding interactions in both actinidin and papain (EC 3.4.22.2) are discussed.

1. 2-(N'-Acetyl-L-phenylalanyl)hydroxyethyl 2′-pyridyl disulphide (compound V) was synthesized, and a study of the pH-dependence of the second-order rate constant (k) for its reaction with the catalytic-site thiol group of papain (EC 3.4.22.2) was used to evaluate the consequences for transition-state geometry of the presence of a hydrophobic occupant for the S2 subsite of the enzyme in the absence of the N-H component of the P1-P2 amide bond. 2. Comparison of the pH-dependences of K for reactions of compound (V), 2-(acetamido)ethyl 2′-pyridyl disulphide (compound I) and 2-(acetoxy)ethyl 2′-pyridyl disulphide (compound III) with the cysteine-proteinase minimal catalytic-site model, benzimidazol-2-ylmethanethiol, established the activation of all of these pyridyl disulphides by hydronation and that their reactivities are relatively insensitive to structural change in the non-pyridyl part of the molecule. The marked differences in their reactivities towards papain therefore derive from binding, either directly, or indirectly via signalling mechanisms. 3. Comparison of the kinetic data for the reaction of papain with compound (V) with those for analogous reactions with reactivity probes that provide opportunities for a variety of binding interactions in the S1-S2 intersubsite region and in the S2 subsite itself lead to the following conclusions: (a) the (Gly-66) N-H...O = C less than (P1-P2 ester) interaction of papain with compound (III) provides for better binding relative to that for a probe with a simple hydrocarbon side chain, but no signalling to the catalytic site to provide a (His-159)-ImH+-assisted transition state; (b) when this interaction is augmented either by a (P1-P2 amide) N-H...O = C less than (Asp-158) interaction (compound I) or hydrophobic P2/S2 contacts (compound V), signalling to the catalytic region occurs to provide the assisted transition state; (c) when both the P2/S2 contacts and the interaction involving Gly-66 exist, provision additionally of the (P1-P2 amide) N-H...O = C less than (Asp-158) interaction [as in 2-(N'-acetyl-L-phenylalanylamino)ethyl 2′-pyridyl disulphide] serves only to assist the binding without an additional signalling effect. 4. Such studies promise to allow binding interactions that merely locate substrates in appropriate enzyme loci to be distinguished from those that transmit signals with a chemical consequence to catalytic sites.

Summary: Sporocytophaga myxococcoides NCIB 8639 utilized a number of cellulosic substrates and produced extracellular carboxymethylcellulase and activity towards Avicel (a milled micro-crystalline cellulose powder). Both types of enzyme were synthesized in media containing glucose but the activity was influenced by the nature and concentration of the carbon source. The ratio of carboxymethylcellulase activities measured by reducing sugar production and by viscosity decrease was also influenced by the nature of the cellulosic substrate supplied to cultures. The greatest extracellular enzyme production occurred in cultures grown on microcrystalline cellulose powders such as Avicel and Whatman CC31 at concentrations up to 4% (w/v). The greatest degradation rate observed was in cultures supplied with Avicel where 60 % of the substrate was degraded in 4 d.

[2-18O]Ribulose 5-phosphate was prepared and shown to be converted enzymically by 6-phosphogluconate dehydrogenase from sheep liver into 6-phosphogluconate with complete retention of the heavy isotope. This finding unequivocally excludes the possibility of a Schiff-base mechanism for the enzyme. The involvement of metal ions has already been excluded, and other possible mechanisms are discussed. The enzyme was purified by an improved large-scale procedure, which is briefly described.

Computer simulations of the irreversible inhibition of an enzyme by an unstable inhibitor are presented. Data obtained at the end point of reaction are shown to conform poorly in many situations with relationships derived from integrated rate equations by setting t = infinity, and the implications concerning the experimental use of this method to determine kinetic constants describing inactivation are considered. The alternative approach of conducting experiments under conditions of inhibitor excess over enzyme is further discussed, and a graphical procedure is suggested for the description of time courses of reaction of enzyme with unstable inhibitor when an enzyme-inhibitor adsorptive complex is involved.

The half-time method for the determination of Michaelis parameters from enzyme progress-curve data (Wharton, C.W. and Szawelski, R.J. (1982) Biochem. J. 203, 351-360) has been adapted for analysis of the kinetics of irreversible enzyme inhibition by an unstable site-specific inhibitor. The method is applicable to a model in which a product (R) of the decomposition of the site-specific reagent, retaining the chemical moiety responsible for inhibitor specificity, binds reversibly to the enzyme with dissociation constant Kr: (formula; see text). Half-time plots of simulated enzyme inactivation time-course data are shown to be unbiased, and excellent estimates of the apparent second-order rate constant for inactivation (k +2/Ki) and Kr can be obtained from a series of experiments with varying initial concentrations of inhibitor. Reliable estimates of k +2 and Ki individually are dependent upon the relative magnitudes of the kinetic parameters describing inactivation. The special case, Kr = Ki, is considered in some detail, and the integrated rate equation describing enzyme inactivation shown to be analogous to that for a simple bimolecular reaction between enzyme and an unstable irreversible inhibitor without the formation of a reversible enzyme-inhibitor complex. The half-time method can be directly extended to the kinetics of enzyme inactivation by an unstable mechanism-based (suicide) inhibitor, provided that the inhibitor is not also a substrate for the enzyme.