Peptide nucleic acid Hoogsteen strand linker design for major groove recognition of DNA thymine bases

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.

Construction of a synthetic pathway for the production of 1,3-propanediol from glucose

In this work, we describe the construction of a synthetic metabolic pathway enabling direct biosynthesis of 1,3-propanediol (PDO) from glucose via the Krebs cycle intermediate malate. This non-natural pathway extends a previously published synthetic pathway for the synthesis of (L)-2,4-dihydroxybutyrate (L-DHB) from malate by three additional reaction steps catalyzed respectively, by a DHB dehydrogenase, a 2-keto-4-hydroxybutyrate (OHB) dehydrogenase and a PDO oxidoreductase. Screening and structure-guided protein engineering provided a (L)-DHB dehydrogenase from the membrane-associated (L)-lactate dehydrogenase of E. coli and OHB decarboxylase variants derived from the branched-chain keto-acid decarboxylase encoded by kdcA from Lactococcus lactis or pyruvate decarboxylase from Zymomonas mobilis. The simultaneous overexpression of the genes encoding these enzymes together with the endogenous ydhD-encoded aldehyde reductase enabled PDO biosynthesis from (L)-DHB. While the simultaneous expression of the six enzymatic activities in a single engineered E. coli strain resulted in a low production of 0.1 mM PDO from 110 mM glucose, a 40-fold increased PDO titer was obtained by co-cultivation of an E. coli strain expressing the malate-DHB pathway with another strain harboring the DHB-to-PDO pathway.

A generic HTS assay for kinase screening: Validation for the isolation of an engineered malate kinase

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 (kcat/Km)malate/(kcat/Km)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.

An Atomistic Statistically Effective Energy Function for Computational Protein Design

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.

Tri-peptide reference structures for the calculation of relative solvent accessible surface area in protein amino acid residues

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.

Probing impact of active site residue mutations on stability and activity of Neisseria polysaccharea amylosucrase

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.

Role of glycoside phosphorylases in mannose foraging by human gut bacteria

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.

Orientation preferences of backbone secondary amide functional groups in peptide nucleic acid complexes: quantum chemical calculations reveal an intrinsic preference of cationic D-amino acid-based chiral PNA analogues for the P-form

Geometric descriptions of nonideal interresidue hydrogen bonding and backbone-base water bridging in the minor groove are established in terms of polyamide backbone carbonyl group orientation from analyses of residue junction conformers in experimentally determined peptide nucleic acid (PNA) complexes. Two types of interresidue hydrogen bonding are identified in PNA conformers in heteroduplexes with nucleic acids that adopt A-like basepair stacking. Quantum chemical calculations on the binding of a water molecule to an O2 base atom in glycine-based PNA thymine dimers indicate that junctions modeled with P-form backbone conformations are lower in energy than a dimer comprising the predominant conformation observed in A-like helices. It is further shown in model systems that PNA analogs based on D-lysine are better able to preorganize in a conformation exclusive to P-form helices than is glycine-based PNA. An intrinsic preference for this conformation is also exhibited by positively charged chiral PNA dimers carrying 3-amino-D-alanine or 4-aza-D-leucine residue units that provide for additional rigidity by side-chain hydrogen bonding to the backbone carbonyl oxygen. Structural modifications stabilizing P-form helices may obviate the need for large heterocycles to target DNA pyrimidine bases via PNA.DNA-PNA triplex formation. Quantum chemical modeling methods are used to propose candidate PNA Hoogsteen strand designs.

Kinetic and titration methods for determination of active site contents of enzyme and catalytic antibody preparations

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.

Direct identification of a peptide binding region in the opioid receptor-like 1 receptor by photoaffinity labeling with [Bpa(10),Tyr(14)]nociceptin

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.

The nociceptin (ORL1) receptor: molecular cloning and functional architecture

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.

On the spatial disposition of the fifth transmembrane helix and the structural integrity of the transmembrane binding site in the opioid and ORL1 G protein-coupled receptor family

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 kinetic basis of a general method for the investigation of active site content of enzymes and catalytic antibodies: first-order behaviour under single-turnover and cycling conditions

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.

Functional inactivation of the nociceptin receptor by alanine substitution of glutamine 286 at the C terminus of transmembrane segment VI: evidence from a site-directed mutagenesis study of the ORL1 receptor transmembrane-binding domain

A site-directed mutagenesis approach has been used to gain insight into the molecular events whereby the heptadecapeptide nociceptin binds and activates the opioid receptor-like 1 (ORL1) receptor, a G protein-coupled receptor. Alanine mutation, in the human ORL1 receptor, of transmembrane amino acid residues that are conserved in opioid receptors, Asp(130) and Tyr(131) in transmembrane segment (TM) III, Phe(220) and Phe(224) in TM V, and Trp(276) in TM VI, yields mutant receptors with reduced affinity, and proportionally decreased reactivity, toward nociceptin. Least to most deleterious in this respect are Ala substitutions of Phe(220) approximately W276A < Tyr(131) << Phe(224) </= Asp(130). The dramatic effects of the D130A mutation on nociceptin binding and activity are not reversed in the D130N mutant, whereas those of the Y131A mutation are totally suppressed in Y131F. This suggests that a negative charge at position 130, and a phenyl ring at position 131 in TM III, are critical for occupancy and/or activation of the receptor by nociceptin. Alanine replacement of glutamine 286, located at the C terminus of TM VI, yields a mutant receptor that binds nociceptin with nearly the same affinity as does the wild-type receptor (K(d) values of 0.13 and 0.22 nM, respectively) but, unlike the latter, is unable to mediate nociceptin inhibition of forskolin-induced cAMP synthesis in recombinant Chinese hamster ovary cells (ED(50) > 10,000 nM compared with 0.8 nM at the wild-type receptor). In all respects, this mutant receptor appears to be functionally inactive, indicating that residue Gln(286) may play a pivotal role in ORL1 receptor-mediated transduction of the nociceptin signal.

A general kinetic approach to investigation of active-site availability in macromolecular catalysts

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].

The influence of helix morphology on co-operative polyamide backbone conformational flexibility in peptide nucleic acid complexes

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 (alpha and epsilon) 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.

Molecular modelling of the ORL1 receptor and its complex with nociceptin

The opioid receptor like (ORL1) receptor is a G-protein coupled receptor superfamily, and regulates a plethora of neurophysiological functions. The structural requirements for receptor activation by its endogenous agonist, nociceptin (FGGFTGARKSARKLANQ), differ markedly from those of the kappa-opioid receptor and its putative peptide agonist, dynorphin A (YGGFLRRIRPKLKWDNQ). In order to probe the functional architecture of the ORL1 receptor, a molecular model of the receptor has been built, including the TM domain and the extra- and intracellular loops. An extended binding site able to accommodate nociceptin-(1-13), the shortest fully active analogue of nociceptin, has been characterized. The N-terminal FGGF tetrapeptide is proposed to bind in a highly conserved region, comprising two distinct hydrophobic pockets in a cavity formed by TM helices 3, 5, 6 and 7, capped by the acidic second extracellular (EL2) loop controlling access to the TM elements of the peptide binding site. The nociceptin conformation provides for the selective preference of the ORL1 receptor for nociceptin over dynorphin A, conferred by residue positions 5 and 6 (TG versus LR), and the favourable interaction of its highly positively charged core (residues 8-13) with the EL2 loop, thought to mediate receptor activation. The functional roles of the EL2 loop and the conserved N-terminal tetrapeptide opioid 'message' binding site are discussed in the context of the different structural requirements of the ORL1 and kappa-opioid receptors for activation.

Prediction of the stability of protein mutants based on structural environment-dependent amino acid substitution and propensity tables

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.

Characterization of the electrostatic perturbation of a catalytic site (Cys)-S-/(His)-Im+H ion-pair in one type of serine proteinase architecture by kinetic and computational studies on chemically mutated subtilisin variants

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.

Comparative modelling of major house dust mite allergen Der p I: structure validation using an extended environmental amino acid propensity table

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.

Structure of chymopapain M the late-eluted chymopapain deduced by comparative modelling techniques and active-centre characteristics determined by pH-dependent kinetics of catalysis and reactions with time-dependent inhibitors: the Cys-25/His-159 ion-pair is insufficient for catalytic competence in both chymopapain M and papain

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)

Ionization characteristics of the Cys-25/His-159 interactive system and of the modulatory group of papain: resolution of ambiguity by electronic perturbation of the quasi-2-mercaptopyridine leaving group in a new pyrimidyl disulphide reactivity probe

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.

Fragment ranking in modelling of protein structure. Conformationally constrained environmental amino acid substitution tables

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 phi-psi 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 x 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 beta-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 chi 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%.

In defence of the general validity of the Cha method of deriving rate equations. The importance of explicit recognition of the thermodynamic box in enzyme kinetics

1. The suggestion by Segel & Martin [(1988) J. Theor. Biol. 135, 445-453] that the well-known schematic method for the derivation of rate equations for combined quasi-equilibrium-steady-state models proposed by Cha [(1968) J. Biol. Chem. 243, 820-825] may not be generally applicable was shown to be incorrect. By contrast, Cha's method was shown (a) to yield correct initial-rate equations that are exact and (b) not to require any constraints on the relative magnitudes of rate constants for slow steps outside the quasi-equilibrium segments of the kinetic model, including those suggested by Segel & Martin. 2. Examination of residual King-Altman patterns for the general modifier model of Botts & Morales [(1953) Trans. Faraday Soc. 49, 696-707] revealed the reasons for the erroneous conclusions reached by Segel & Martin. The errors arise from the failure to take account of fluxes in parallel pathways that connect the two isolated groups of enzyme species existing in quasi-equilibrium with the modifier. 3. A similar failure to take explicit account of parallel pathways in a thermodynamic box that delayed proper appreciation of the form of pH-dependence of kcat/Km is briefly discussed.

Structure-function relationships in the cysteine proteinases actinidin, papain and papaya proteinase omega. Three-dimensional structure of papaya proteinase omega deduced by knowledge-based modelling and active-centre characteristics determined by two-hydronic-state reactivity probe kinetics and kinetics of catalysis

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)

A polyclonal antibody preparation with Michaelian catalytic properties

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 modifier

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.

A re-appraisal of the structural basis of stereochemical recognition in papain. Insensitivity of binding-site-catalytic-site signalling to P2-chirality in a time-dependent inhibition

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)

Dependence of the P2-S2 stereochemical selectivity of papain on the nature of the catalytic-site chemistry. Quantification of selectivity in the catalysed hydrolysis of the enantiomeric N-acetylphenylalanylglycine 4-nitroanilides

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.

The interplay of electrostatic fields and binding interactions determining catalytic-site reactivity in actinidin. A possible origin of differences in the behaviour of actinidin and papain

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.

Identification of signalling and non-signalling binding contributions to enzyme reactivity. Alternative combinations of binding interactions provide for change in transition-state geometry in reactions of papain

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.

Supracrystallographic resolution of interactions contributing to enzyme catalysis by use of natural structural variants and reactivity-probe kinetics

1. The influence on the reactivities of the catalytic sites of papain (EC 3.4.22.2) and actinidin (3.4.22.14) of providing for interactions involving the S1-S2 intersubsite regions of the enzymes was evaluated by using as a series of thiol-specific two-hydronic-state reactivity probes: n-propyl 2-pyridyl disulphide (I) (a 'featureless' probe), 2-(acetamido)ethyl 2'-pyridyl disulphide (II) (containing a P1-P2 amide bond), 2-(acetoxy)ethyl 2'-pyridyl disulphide (III) [the ester analogue of probe (II)] and 2-carboxyethyl 2'-pyridyl disulphide N-methylamide (IV) [the retroamide analogue of probe (II)]. Syntheses of compounds (I), (III) and (IV) are reported. 2. The reactivities of the two enzymes towards the four reactivity probes (I)-(IV) and also that of papain towards 2-(N'-acetyl-L-phenylalanylamino)ethyl 2'-pyridyl disulphide (VII) (containing both a P1-P2 amide bond and an L-phenylalanyl side chain as an occupant for the S2 subsite), in up to four hydronic (previously called protonic) states, were evaluated by analysis of pH-dependent stopped-flow kinetic data (for the release of pyridine-2-thione) by using an eight-parameter rate equation [described in the Appendix: Brocklehurst & Brocklehurst (1988) Biochem. J. 256, 556-558] to provide pH-independent rate constants and macroscopic pKa values. The analysis reveals the various ways in which the two enzymes respond very differently to the binding of ligands in the S1-S2 intersubsite regions despite the virtually superimposable crystal structures in these regions of the molecules. 3. Particularly striking differences between the behaviour of papain and that of actinidin are that (a) only papain responds to the presence of a P1-P2 amide bond in the probe such that a rate maximum at pH 6-7 is produced in the pH-k profile in place of the rate minimum, (b) only in the papain reactions does the pKa value of the alkaline limb of the pH-k profile change from 9.5 to approx. 8.2 [the value characteristic of a pH-(kcat./Km) profile] when the probe contains a P1-P2 amide bond, (c) only papain reactivity is affected by two positively co-operative hydronic dissociations with pKI congruent to pKII congruent to 4 and (d) modulation of the reactivity of the common -S(-)-ImH+ catalytic-site ion-pair (Cys-25/His-159 in papain and Cys-25/His-162 in actinidin) by hydronic dissociation with pKa approx. 5 is more marked and occurs more generally in reactions of actinidin than is the case for papain reactions.(ABSTRACT TRUNCATED AT 400 WORDS)

Ill-conditioning associated with the "end-point" method for the determination of kinetic parameters describing irreversible enzyme inactivation by an unstable inhibitor

Study of the complete time-course of irreversible enzyme inhibition by an unstable inhibitor yields more information than can be obtained by recording data only at the end point of reaction. Time-course analysis of co-operative irreversible enzyme inhibition by an unstable inhibitor has been shown to be considerably less susceptible to ill-conditioning than the "end-point" method for the determination of kinetic parameters describing inactivation. As a result, mechanisms that cannot be distinguished by the "end-point" method are readily differentiated by time-course analysis without the need to isolate intermediate species.

Half-time analysis of the kinetics of irreversible enzyme inhibition by an unstable site-specific reagent

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.

Computer simulations of the kinetics of irreversible enzyme inhibition by an unstable inhibitor

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.

Kinetic studies of 6-phosphogluconate dehydrogenase from sheep liver

The steady-state kinetics of the oxidative decarboxylation of 6-phosphogluconate catalysed by 6-phosphogluconate dehydrogenase from sheep liver in triethanolamine and phosphate buffers (pH 7.0) have been reinvestigated. In triethanolamine buffer the enzyme is inhibited by high NADP+ concentrations in the presence of low fixed concentrations of 6-phosphogluconate. Data are consistent with an asymmetric sequential mechanism in which NADP+ and 6-phosphogluconate bind randomly and product release is ordered. The pathway through the enzyme--6-phosphogluconate complex appears to be preferred in triethanolamine buffer. Pre-steady-state studies of the oxidative decarboxylation reaction at pH 6.0-8.0 show that hydride transfer is greater than 900 s-1. After the fast formation of NADPH in amounts equivalent to about half of the enzyme-active-centre concentration, the rate of NADPH formation is equal to the steady-state rate. Two possible interpretations are considered. Rapid fluorescence measurements of the displacement of NADPH from its complex with the enzyme at pH 6.0 and 7.0 indicate that the dissociation of NADPH is fast (greater than 800 s-1) and cannot be the rate-limiting step in oxidative decarboxylation. Coenzyme binding studies at equilibrium have been extended to include the determination of the dissociation constants for the binary complexes of enzyme with NADPH and NADP+ at pH 6.0-8.0 and the dissociation constant for NADPH in the ternary enzyme--6-phosphogluconate--NADPH complex in triethanolamine buffer, pH 7.0.

The chemical mechanism of sheep liver 6-phosphogluconate dehydrogenase. A Schiff-base intermediate is not involved

[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.

Chemical modification of sheep-liver 6-phosphogluconate dehydrogenase by diethylpyrocarbonate. Evidence for an essential histidine residue

Sheep liver 6-phosphogluconate dehydrogenase is shown to be inactivated by diethylpyrocarbonate in a biphasic manner at pH 6.0, 25 degrees C. After allowing for the hydrolysis of the reagent, rate constants of 56 M-1 s-1 and 11.0 M-1 s-1 were estimated for the two processes. The complete reactivation of partially inactivated enzyme by neutral hydroxylamine, the elimination of the possibility that modification of cysteine or tyrosine residues are responsible for inactivation, and the magnitudes of the rate constants for inactivation relative to the experimentally determined value for the reaction of diethylpyrocarbonate with N alpha-acetylhistidine (2.2 M-1 s-1), all suggested that enzyme inactivation occurs solely by modification of histidine residues. Comparison of the experimental plot of residual fractional activity versus the number of modified histidine residues per subunit with simulated plots for three hypothetical models, each predicting biphasic kinetics, indicated that inactivation results from the modification of at most one essential histidine residue per subunit, although it appears that other (non-essential) histidines react independently. This histidine is thought to be His-242 and is present in the active site. Evidence in support of its role in catalysis is briefly discussed. Both 6-phosphogluconate and organic phosphate protect against inactivation, and a kinetic analysis of the protection indicated a dissociation constant of 2.1 X 10(-6) M for the enzyme--6-phosphogluconate complex. NADP+ also protected, but this might be due, at least in part, to a reduction in the effective concentration of diethylpyrocarbonate.