Stereochemical analysis of peptide bond hydrolysis catalyzed by the aspartic proteinase penicillopepsin

The X-ray crystal structures of native penicillopepsin and of its complex with a synthetic analogue of the inhibitor pepstatin have been refined recently at 1.8-A resolution. These highly refined structures permit a detailed examination of peptide hydrolysis in the aspartic proteinases. Complexes of penicillopepsin with substrate and catalytic intermediates were modeled, by using computer graphics, with minimal perturbation of the observed inhibitor complex. A thallium ion binding experiment shows that the position of solvent molecule O39, between Asp-33(32) and Asp-213(215) in the native structure, is favorable for cations, a fact that places constraints on possible mechanisms. A mechanism for hydrolysis is proposed in which Asp-213(215) acts as an electrophile by protonating the carbonyl oxygen of the substrate, thereby polarizing the carbon-oxygen bond, a water molecule bound to Asp-33(32) (O284 in the native structure) attacks the carbonyl carbon as the nucleophile in a general-base mechanism, the newly pyramidal peptide nitrogen is protonated, either from the solvent after nitrogen inversion or by an internal proton transfer via Asp-213(215) from a hydroxyl of the tetrahedral carbon, and the tetrahedral intermediate breaks down in a manner consistent with the stereoelectronic hypothesis. The models permit the rationalization of observed subsite preferences for substrates and may be useful in predicting subsite preferences of other aspartic proteinases.     

Five atomic resolution structures of endothiapepsin inhibitor complexes: implications for the aspartic proteinase mechanism

Endothiapepsin is derived from the fungus Endothia parasitica and is a member of the aspartic proteinase class of enzymes. This class of enzyme is comprised of two structurally similar lobes, each lobe contributing an aspartic acid residue to form a catalytic dyad that acts to cleave the substrate peptide bond. The three-dimensional structures of endothiapepsin bound to five transition state analogue inhibitors (H189, H256, CP-80,794, PD-129,541 and PD-130,328) have been solved at atomic resolution allowing full anisotropic modelling of each complex. The active sites of the five structures have been studied with a view to studying the catalytic mechanism of the aspartic proteinases by locating the active site protons by carboxyl bond length differences and electron density analysis. In the CP-80,794 structure there is excellent electron density for the hydrogen on the inhibitory statine hydroxyl group which forms a hydrogen bond with the inner oxygen of Asp32. The location of this proton has implications for the catalytic mechanism of the aspartic proteinases as it is consistent with the proposed mechanism in which Asp32 is the negatively charged aspartate. A number of short hydrogen bonds (approximately 2.6 A) with ESD values of around 0.01 A that may have a role in catalysis have been identified within the active site of each structure; the lengths of these bonds have been confirmed using NMR techniques. The possibility and implications of low barrier hydrogen bonds in the active site are considered.     

Structural insights into the activation and inhibition of histo-aspartic protease from Plasmodium falciparum

Histo-aspartic protease (HAP) from Plasmodium falciparum is a promising target for the development of novel antimalarial drugs. The sequence of HAP is highly similar to those of pepsin-like aspartic proteases, but one of the two catalytic aspartates, Asp32, is replaced with histidine. Crystal structures of the truncated zymogen of HAP and of the complex of the mature enzyme with inhibitor KNI-10395 have been determined at 2.1 and 2.5 ? resolution, respectively. As in other proplasmepsins, the propeptide of the zymogen interacts with the C-terminal domain of the enzyme, forcing the N- and C-terminal domains apart, thereby separating His32 and Asp215 and preventing formation of the mature active site. In the inhibitor complex, the enzyme forms a tight domain-swapped dimer, not previously seen in any aspartic proteases. The inhibitor is found in an unprecedented conformation resembling the letter U, stabilized by two intramolecular hydrogen bonds. Surprisingly, the location and conformation of the inhibitor are similar to those of the fragment of helix 2 comprising residues 34p-38p in the prosegments of the zymogens of gastric aspartic proteases; a corresponding helix assumes a vastly different orientation in proplasmepsins. Each inhibitor molecule is in contact with two molecules of HAP, interacting with the carboxylate group of the catalytic Asp215 of one HAP protomer through a water molecule, while also making a direct hydrogen bond to Glu278A' of the other protomer. A comparison of the shifts in the positions of the catalytic residues in the inhibitor complex presented here with those published previously gives further hints regarding the enzymatic mechanism of HAP.     

Strongly hydrogen-bonded water molecules in the Schiff base region of rhodopsins.

In many rhodopsins, a positively charged retinal chromophore is stabilized by a negatively charged carboxylate, and the presence of bound water molecules has been found in the Schiff base region by X-ray crystallography of various rhodopsins. Low-temperature Fourier-transform infrared (FTIR) spectroscopy can directly monitor hydrogen-bonding alterations of internal water molecules of rhodopsins. In particular, we found that a bridged water molecule between the Schiff base and Asp 85 in bacteriorhodopsin (BR), a light-driven proton-pump protein, forms an extremely strong hydrogen bond. It is likely that a hydration switch of the water from Asp 85 to Asp 212 plays an important role in the proton transfer in the Schiff base region of BR. Comprehensive studies of archaeal and visual rhodopsins have revealed that strongly hydrogen-bonded water molecules are only found in the proteins exhibiting proton-pump activities. Strongly hydrogen-bonded water molecules and its transient weakening may be essential for the proton-pump function of rhodopsins.     

The catalytic mechanism of aspartic proteinases

The highly symmetric active site of an aspartic proteinase, endothiapepsin, binds a water molecule ideally situated for nucleophilic attack on a substrate peptide bond whose distortion from planarity is stabilised by interactions of the substrate with the extended binding cleft. The apparent electrophilicity of the catalysis results from this distortion. The scissile peptide bond is orientated with the carbonyl oxygen hydrogen bonding to the tip of the beta-hairpin 'flap' which lies over the cleft. Nucleophilic attack by the bound water leads to a tetrahedral intermediate similar to observed complexes with hydroxyl inhibitors and stabilised by hydrogen bonds with the flap.     

Caught in the Act: the 1.5 A resolution crystal structures of the HIV-1 protease and the I54V mutant reveal a tetrahedral reaction intermediate

HIV-1 protease (PR) is the target for several important antiviral drugs used in AIDS therapy. The drugs bind inside the active site cavity of PR where normally the viral polyprotein substrate is bound and hydrolyzed. We report two high-resolution crystal structures of wild-type PR (PRWT) and the multi-drug-resistant variant with the I54V mutation (PRI54V) in complex with a peptide at 1.46 and 1.50 A resolution, respectively. The peptide forms a gem-diol tetrahedral reaction intermediate (TI) in the crystal structures. Distinctive interactions are observed for the TI binding in the active site cavity of PRWT and PRI54V. The mutant PRI54V/TI complex has lost water-mediated hydrogen bond interactions with the amides of Ile50 and Ile50' in the flap. Hence, the structures provide insight into the mechanism of drug resistance arising from this mutation. The structures also illustrate an intermediate state in the hydrolysis reaction. One of the gem-diol hydroxide groups in the PRWT complex forms a very short (2.3 A) hydrogen bond with the outer carboxylate oxygen of Asp25. Quantum chemical calculations based on this TI structure are consistent with protonation of the inner carboxylate oxygen of Asp25', in contrast to several theoretical studies. These TI complexes and quantum calculations are discussed in relation to the chemical mechanism of the peptide bond hydrolysis catalyzed by PR.     

A neutron Laue diffraction study of endothiapepsin: implications for the aspartic proteinase mechanism.

Current proposals for the catalytic mechanism of aspartic proteinases are largely based on X-ray structures of bound oligopeptide inhibitors possessing nonhydrolyzable analogues of the scissile peptide bond. However, the positions of protons on the catalytic aspartates and the ligand in these complexes have not been determined with certainty. Thus, our objective was to locate crucial protons at the active site of an inhibitor complex since this will have major implications for a detailed understanding of the mechanism of action. We have demonstrated that high-resolution neutron diffraction data can be collected from crystals of the fungal aspartic proteinase endothiapepsin bound to a transition state analogue (H261). The neutron structure of the complex has been refined at a resolution of 2.1 A to an R-factor of 23.5% and an R(free) of 27.4%. This work represents the largest protein structure studied to date by neutron crystallography at high resolution. The neutron data demonstrate that 49% of the main chain nitrogens have exchanged their hydrogen atoms with D2O in the mother liquor. The majority of residues resisting exchange are buried within core beta-sheet regions of the molecule. The neutron maps confirm that the protein has a number of buried ionized carboxylate groups which are likely to give the molecule a net negative charge even at very low pH, thereby accounting for its low pI. The functional groups at the catalytic center have clearly undergone H-D exchange despite being buried by the inhibitor occupying the active site cleft. Most importantly, the data provide convincing evidence that Asp 215 is protonated and that Asp 32 is the negatively charged residue in the transition state complex. This has an important bearing on mechanistic proposals for this class of proteinase.     

Crystallographic analysis of transition state mimics bound to penicillopepsin: difluorostatine- and difluorostatone-containing peptides.

Difluorostatine- and difluorostatone-containing peptides have been evaluated as potent inhibitors of penicillopepsin, a member of the aspartic proteinase family of enzymes. Isovaleryl-Val-Val-StaF2NHCH3 [StaF2 = (S)-4-amino-2,2-difluoro-(R)-3-hydroxy-6-methylheptanoic acid] and isovaleryl-Val-Val-StoF2NHCH3 [StoF2 = (S)-4-amino-2,2-difluoro-3-oxo-6-methylheptanoic acid] have measured Ki's of 10 x 10(-9) and 1 x 10(-9) M, respectively, with this fungal proteinase. The StoF2-containing peptide binds 32-fold more tightly to the enzyme than the analogous peptide containing the non-fluorinated statine ethyl ester. Each compound was cocrystallized with penicillopepsin, intensity data were collected to 1.8-A resolution, and the atomic coordinates were refined to an R factor [formula: see text] of 0.131 for both complexes. The inhibitors bind in the active site of penicillopepsin in much the same fashion as do other statine-containing inhibitors of penicillopepsin analyzed earlier [James, M. N. G., Sielecki, A. Salituro, F., Rich, D. H., & Hofmann, T. (1982) Proc. Natl. Acad. Sci. U.S.A. 79, 6137-6141; James, M.N.G., Sielecki, A., & Hofmann, T. (1985) in Aspartic Proteinases and their Inhibitors (Kosta, V., Ed.) pp 163-177, Walter deGruyter, Berlin]. The (R)-3-hydroxyl group in StaF2 binds between the active site carboxyl groups of Asp33 and Asp213, making hydrogen-bonding contacts to each one. The ketone functional group of the StoF2 inhibitor is bound as a hydrated species, with the gem-diol situated between the two aspartic acid carboxyl groups in a manner similar to that predicted for the tetrahedral intermediate expected during the catalytic hydrolysis of a peptide bond [James, M. N. G., & Sielecki, A. (1985) Biochemistry 24, 3701-3713]. One hydrogen-bonding interaction from the "outer" hydroxyl group is made to O delta 1 of Asp33, and the "inner" hydroxyl group forms two hydrogen-bonding contacts, one to each of the carboxyl groups of Asp33 (O delta 2) and Asp213 (O delta 2). The only structural difference between the StaF2 and StoF2 inhibitors that accounts for the factor of 10 in their Ki's is the additional (R)-3-OH group on the tetrahedral sp3 carbon atom of the hydrated StoF2 inhibitor. The intermolecular interactions involving the fluorine atoms of each inhibitor are normal van der Waals contacts to one of the carboxyl oxygen atoms of Asp213 (F2-O delta 2 Asp213, 2.9 A). The observed stereochemistry of the bound StoF2 group in the active site of penicillopepsin has stimulated our reappraisal of the catalytic pathway for the aspartic proteinases.(ABSTRACT TRUNCATED AT 400 WORDS)     

A Combined Quantum/Classical Molecular Dynamics Study of the Catalytic Mechanism of HIV Protease

Based on available three-dimensional structures of enzyme-inhibitor complexes, the mechanism of the reaction catalysed by HIV protease is studied using molecular dynamics simulations with molecular mechanics and combined quantum-mechanics/molecular-mechanics potential energy functions. The results support the general acid/general base catalysis mechanism, with Asp25′ protonated in the enzyme-substrate complex. In the enzyme-substrate complex, the lytic water molecule binds at a position different from the positions of the hydroxyl groups in various aspartic protease-inhibitor complexes. The carboxyl groups at the active site also adopt a different orientation. However, when the lytic water molecule approaches the scissile peptide, the reaction centre changes gradually to a conformation close to that derived from X-ray diffraction studies of various enzyme-inhibitor complexes. The proton transfer processes can take place only after the lytic water molecule has approached the scissile peptide bond to a certain degree. Qualitatively, the free-energy barrier associated with the nucleophilic attack step, which takes place at physiological pH, is comparable with the acid or base-catalysed reactions of model systems. The structure of the tetrahedral intermediate resulting from the nucleophilic attack step also indicates a straightforward pathway of the next reaction step, i.e. the breaking of the C-N bond.     

NMR study of human mutant hemoglobins synthesized in Escherichia coli. Consequences of tyrosine alpha 42 substitutions

The hydroxyl group of Tyr alpha 42 in human hemoglobin forms a hydrogen bond with the carboxylate of Asp beta 99 which is considered to be one of the most important hydrogen bonds for stabilizing the "T-state." However, no spontaneous mutation at position 42 of the alpha subunit has been reported, and the role of the tyrosine has not been tested experimentally. Two artificial human mutant hemoglobins in which Tyr alpha 42 was replaced by phenylalanine or histidine were synthesized in Escherichia coli, and their proton NMR spectra were studied with particular attention to the hyperfine-shifted and hydrogen-bonded proton resonances. The site-directed mutagenesis of the Tyr alpha 42----Phe removes the hydrogen bond described above and prevents transition to the T-state so that the mutant Hb is rather similar to the "R-state" even when deoxygenated. On the other hand, the mutation from tyrosine to histidine causes less drastic structural changes, and its quaternary and tertiary structures are almost the same as native deoxy-Hb A. This may be attributed to the formation of a new hydrogen bond between His alpha 1(42) and Asp beta 2(99). These observations indicate that the hydrogen bond formed between Tyr alpha 42 and Asp beta 99 is required to convert unliganded Hb to the T-state.     

Functional plasticity in the substrate binding site of beta-secretase

The aspartic protease beta-secretase (BACE) cleaves the amyloid precursor protein into a 42 residue beta-peptide, which is the principal biochemical marker of Alzheimer's disease. Multiple explicit-water molecular dynamics simulations of the apo and inhibitor bound structures of BACE indicate that both open- and closed-flap conformations are accessible at room temperature and should be taken into account for inhibitor design. Correlated motion is observed within each of the two lobes of BACE, as well as for the interfacial region. A self-inhibited conformation with the side chain of Tyr71 occupying the S(1) pocket is present in some of the unbound simulations. The reversible loss of the side chain hydrogen bond between the catalytic Asp32 and Ser35, due to the concomitant reorientation of the Ser35 hydroxyl group and a water molecule conserved in pepsin-like enzymes, provides further evidence for the suggestion that Ser35 assists in proton acceptance and release by Asp32 during catalysis.     

The role of hydrogen bonding in the enzymatic reaction catalyzed by HIV-1 protease

The hydrogen-bond network in various stages of the enzymatic reaction catalyzed by HIV-1 protease was studied through quantum-classical molecular dynamics simulations. The approximate valence bond method was applied to the active site atoms participating directly in the rearrangement of chemical bonds. The rest of the protein with explicit solvent was treated with a classical molecular mechanics model. Two possible mechanisms were studied, general-acid/general-base (GA/GB) with Asp 25 protonated at the inner oxygen, and a direct nucleophilic attack by Asp 25. Strong hydrogen bonds leading to spontaneous proton transfers were observed in both reaction paths. A single-well hydrogen bond was formed between the peptide nitrogen and outer oxygen of Asp 125. The proton was diffusely distributed with an average central position and transferred back and forth on a picosecond scale. In both mechanisms, this interaction helped change the peptide-bond hybridization, increased the partial charge on peptidyl carbon, and in the GA/GB mechanism, helped deprotonate the water molecule. The inner oxygens of the aspartic dyad formed a low-barrier, but asymmetric hydrogen bond; the proton was not positioned midway and made a slightly elongated covalent bond, transferring from one to the other aspartate. In the GA/GB mechanism both aspartates may help deprotonate the water molecule. We observed the breakage of the peptide bond and found that the protonation of the peptidyl amine group was essential for the peptide-bond cleavage. In studies of the direct nucleophilic mechanism, the peptide carbon of the substrate and oxygen of Asp 25 approached as close as 2.3 A.     

Aspartic proteinases: Fourier transform infrared spectroscopic studies of a model of the active side.

We synthesized and studied by Fourier transform infrared spectroscopy nine monosalts of diamides as models for the active side of aspartic proteinases. One compound, the monosalt of meta-aminobenzoic acid diamide of fumaric acid (m-FUM), shows the same biological activity as pepsin with regard to the splitting of peptide bonds of the Pro-Thi-Glu-Phe-Phe(4-NO2)-Arg-Leu heptapeptide. The monosalt of m-FUM forms with oxindole a complex in which the carboxylic acid group of the monosalt of m-FUM is strongly hydrogen bonded with the O atom of the peptide bond of oxindole. When one water molecule is added to this complex, the strong field of the carboxylate group destabilizes an O-H bond of the water molecule. The distorted water molecule attacks the carbon atom of the peptide group, and the water proton transfers to the peptide N atom. Simultaneously, the C-N bond of the amide group is broken. Hence it is demonstrated that the catalytic mechanism of aspartic acid proteinases is a base catalysis. The results show that for this catalytic mechanism there are sufficient carboxylic and carboxylate groups, as well as a water molecule in the correct arrangement. It was also demonstrated with other monosalts of dicarboxylic acids that well-defined steric conditions of the carboxylic acid and the carboxylate group must be fulfilled to show hydrolytic activity with regard to oxindole molecules.     

Binding of alkylurea inhibitors to epoxide hydrolase implicates active site tyrosines in substrate activation

The structures of two alkylurea inhibitors complexed with murine soluble epoxide hydrolase have been determined by x-ray crystallographic methods. The alkyl substituents of each inhibitor make extensive hydrophobic contacts in the soluble epoxide hydrolase active site, and each urea carbonyl oxygen accepts hydrogen bonds from the phenolic hydroxyl groups of Tyr(381) and Tyr(465). These hydrogen bond interactions suggest that Tyr(381) and/or Tyr(465) are general acid catalysts that facilitate epoxide ring opening in the first step of the hydrolysis reaction; Tyr(465) is highly conserved among all epoxide hydrolases, and Tyr(381) is conserved among the soluble epoxide hydrolases. In one enzyme-inhibitor complex, the urea carbonyl oxygen additionally interacts with Gln(382). If a comparable interaction occurs in catalysis, then Gln(382) may provide electrostatic stabilization of partial negative charge on the epoxide oxygen. The carboxylate side chain of Asp(333) accepts a hydrogen bond from one of the urea NH groups in each enzyme-inhibitor complex. Because Asp(333) is the catalytic nucleophile, its interaction with the partial positive charge on the urea NH group mimics its approach toward the partial positive charge on the electrophilic carbon of an epoxide substrate. Accordingly, alkylurea inhibitors mimic features encountered in the reaction coordinate of epoxide ring opening, and a structure-based mechanism is proposed for leukotoxin epoxide hydrolysis.     

Conformation and intramolecular hydrogen-bonding in the crystal structure of potassium D-gluconate monohydrate

The D-gluconate ion is found to have the planar, extended carbon-chain conformation in the crystal structure of potassium D-gluconate monohydrate, with an intramolecular hydrogen-bond between 0-2 and 0-4. The D-gluconate ions and water molecules are linked in puckered sheets by a series of intermolecular hydrogen-bonds that involve the water molecules, the carboxylate groups, and pairs of hydroxyl groups. One hydroxyl group in the ion does not form a hydrogen bond. The potassium ions lie between the puckered sheets, with an eight-fold coordination of six D-gluconate groups and two water oxygen atoms. The crystal structure was determined from three-dimensional, CuKα, X-ray diffraction data taken on an automatic diffractometer.     

Active-site mobility in human immunodeficiency virus, type 1, protease as demonstrated by crystal structure of A28S mutant.

The mutation Ala28 to serine in human immunodeficiency virus, type 1, (HIV-1) protease introduces putative hydrogen bonds to each active-site carboxyl group. These hydrogen bonds are ubiquitous in pepsin-like eukaryotic aspartic proteases. In order to understand the significance of this difference between HIV-1 protease and homologous, eukaryotic aspartic proteases, we solved the three-dimensional structure of A28S mutant HIV-1 protease in complex with a peptidic inhibitor U-89360E. The structure has been determined to 2.0 A resolution with an R factor of 0.194. Comparison of the mutant enzyme structure with that of the wild-type HIV-1 protease bound to the same inhibitor (Hong L, Treharne A, Hartsuck JA, Foundling S, Tang J, 1996, Biochemistry 35:10627-10633) revealed double occupancy for the Ser28 hydroxyl group, which forms a hydrogen bond either to one of the oxygen atoms of the active-site carboxyl or to the carbonyl oxygen of Asp30. We also observed marked changes in orientation of the Asp25 catalytic carboxyl groups, presumably caused by the new hydrogen bonds. These observations suggest that catalytic aspartyl groups of HIV-1 protease have significant conformational flexibility unseen in eukaryotic aspartic proteases. This difference may provide an explanation for some unique catalytic properties of HIV-1 protease.     

Quantum chemical study of the "catalytic triad" of serine proteases

Using the semi-empirical MNDO/H method several systems simulating the reaction of tetrahedral intermediate formation in the active site of serine proteases have been studied. The role played by elements of the "catalytic triad" in increasing the reactivity of serine hydroxyl has been discussed. The formation of a strong hydrogen bond between His and Asp was shown to be important in lowering the activation energy in the reaction of Ser with substrate. The change in position of the proton located between Ser and His and between His and Asp was analysed. The influence of substrate distortion on the energy of intermediate formation has been considered.     

X-ray structure of HIV-1 protease in situ product complex

HIV-1 protease is an effective target for design of different types of drugs against AIDS. HIV-1 protease is also one of the few enzymes that can cleave substrates containing both proline and nonproline residues at the cleavage site. We report here the first structure of HIV-1 protease complexed with the product peptides SQNY and PIV derived by in situ cleavage of the oligopeptide substrate SQNYPIV, within the crystals. In the structure, refined against 2.0-A resolution synchrotron data, a carboxyl oxygen of SQNY is hydrogen-bonded with the N-terminal nitrogen atom of PIV. At the same time, this proline nitrogen atom does not form any hydrogen bond with catalytic aspartates. These two observations suggest that the protonation of scissile nitrogen, during peptide bond cleavage, is by a gem-hydroxyl of the tetrahedral intermediate rather than by a catalytic aspartic acid.     

Structure of cobalt carbonic anhydrase complexed with bicarbonate.

The three-dimensional structure of a complex between catalytically active cobalt(II) substituted human carbonic anhydrase II and its substrate bicarbonate was determined by X-ray crystallography (1.9 A). One water molecule and two bicarbonate oxygen atoms are found at distances between 2.3 and 2.5 A from the cobalt ion in addition to the three histidyl ligands contributed by the peptide chain. The tetrahedral geometry around the metal ion in the native enzyme with a single water molecule 2.0 A from the metal is therefore lost. The geometry is difficult to classify but might best be described as distorted octahedral. The structure is suggested to represent a water-bicarbonate exchange state relevant also for native carbonic anhydrase, where the two unprotonized oxygen atoms of the substrate are bound in a carboxylate binding site and the hydroxyl group is free to move closer to the metal thereby replacing the metal-bound water molecule. A reaction mechanism based on crystallographically determined enzyme-ligand complexes is represented.     

Protein engineering of chymosin; modification of the optimum pH of enzyme catalysis

The aspartic proteinase chymosin exhibits a local network of hydrogen bonds involving the active site aspartates and surrounding residues which may have an influence on the rate and optimal pH of substrate cleavage. We have introduced into chymosin B the following substitutions: Asp304 to Ala (D304A), Thr218 to Ala (T218A) and Gly244 to Asp (G244D, chymosin A), using oligonucleotide-directed mutagenesis. Kinetic analysis of these active mutants shows shifts in their pH optima to 4.4 D304A, 4.2 T218A and 4.0 G244D compared with 3.8 for chymosin B using a synthetic octapeptide substrate. The upward shift of the D304A and T218A may be due to the loss of hydrogen bond interactions indirectly affecting the catalytic aspartates 32 and 215. The G244D mutation which is in a flexible loop on the surface of the enzyme may alter the conformation of the specificity pockets on the prime side of the scissile bond.