An engineered PET depolymerase to break down and recycle plastic bottles.

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.

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.

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.

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

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.

An assessment of COMPOSER: a rule-based approach to modelling protein structure.

We describe COMPOSER, a computer program for modelling proteins which uses three-dimensional structures defined by X-ray analysis together with rules defined by their analysis and comparison. We discuss the modelling of phospholipases and serine proteinases of known structure to assess the errors associated with such models.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.