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.     

Theoretical calculations on the acidity of the active site in aspartic proteinases

Semiempirical minimal neglect of differential overlap-self-consistent field calculations, corrected and modified for multiple hydrogen-bonding interactions, were applied to models of the active site of aspartic proteinases (AP). The propensities of the two active-site aspartates to ionize were compared under the influence of various neighboring residues and of water molecules. Asp-32 and Asp-215 in three aspartic proteinases (endothiapepsin, Rhizopus pepsin, and penicillopepsin) are found to be basically asymmetric, Asp-32 being preferentially (by 2-3 kcal) ionized with respect to Asp-215. In penicillopepsin, this asymmetry is compensated by effects of surrounding residues. In our largest model for the active site, which includes such other residues, near equality is found for the ionizing tendency of Asp-32 and Asp-215. The pK difference is rationalized in terms of first and second ionizations of the full active-site model. Its ionization enthalpies correlate well with those of other small organic diacids. This "gas-phase" approach to AP active-site interactions represents the main possible contributions to the acidity of the active site.     

X-ray analyses of aspartic proteinases. The three-dimensional structure at 2.1 A resolution of endothiapepsin

The molecular structure of endothiapepsin (EC 3.4.23.6), the aspartic proteinase from Endothia parasitica, has been refined to a crystallographic R-factor of 0.178 at 2.1 A resolution. The positions of 2389 protein non-hydrogen atoms have been determined and the present model contains 333 solvent molecules. The structure is bilobal, consisting of two predominantly beta-sheet domains that are related by an approximate 2-fold axis. Of approximately 170 residues, 65 are topologically equivalent when one lobe is superimposed on the other. Twenty beta-strands are arranged as five beta-sheets and are connected by regions involving 29 turns and four helices. A central sheet involves three antiparallel strands from each lobe organized around the dyad axis. Each lobe contains a further local dyad that passes through two sheets arranged as a sandwich and relates two equivalent motifs of four antiparallel strands (a, b, c, d) followed by a helix or an irregular helical region. Sheets 1N and 1C, each contain two interpenetrating psi structures contributed by strands c,d,d' and c',d',d, which are related by the intralobe dyad. A further sheet, 2N or 2C, is formed from two extended beta-hairpins from strands b,c and b',c' that fold above the sheets 1N and 1C, respectively, and are hydrogen-bonded around the local intralobe dyad. Asp32 and Asp215 are related by the interlobe dyad and form an intricate hydrogen-bonded network with the neighbouring residues and comprise the most symmetrical part of the structure. The side-chains of the active site aspartate residues are held coplanar and the nearby main chain makes a "fireman's grip" hydrogen-bonding network. Residues 74 to 83 from strands a'N and b'N in the N-terminal lobe form a beta-hairpin loop with high thermal parameters. This "flap" projects over the active site cleft and shields the active site from the solvent region. Shells of water molecules are found on the surface of the protein molecule and large solvent channels are observed within the crystal. There are only three regions of intermolecular contacts and the crystal packing is stabilized by many solvent molecules forming a network of hydrogen bonds. The three-dimensional structure of endothiapepsin is found to be similar to two other fungal aspartic proteinases, penicillopepsin and rhizopuspepsin. Even though sequence identities of endothiapepsin with rhizopuspepsin and penicillopepsin are only 41% and 51%, respectively, a superposition of the three-dimensional structures of these three enzymes shows that 237 residues (72%) are within a root-mean-square distance of 1.0 A.     

Quantum mechanical modeling of aspartic proteinase interactions: difference in binding of diastereomeric statine models

Quantum mechanical calculations were carried out for the interaction of two diastereomeric model inhibitors of aspartic proteinases with a model for the active site, based on crystallographic coordinates of endothiapepsin. The model inhibitor is formamide N-(2-hydroxy 3-methyl propane) and the active site is represented by the full backbone and most of the side chains of the two partial sequences D32-T33-G34-S35 and D215-T216-G217-T218. Those calculations demonstrate that the best binding mode for this short inhibitor is consistent with the X-ray experiments and somewhat stronger with the inhibitor in a 2(S) configuration, compared to 2(R). Another binding mode is possible for this model inhibitor only in the 2(S)-configuration, and is weaker than the first.     

Aspartic proteinases in Antarctic fish

The present review surveys several recent studies of the aspartic proteinases from Antarctic Notothenioidei, a dominating fish group that has developed a number of adjustments at the molecular level to maintain metabolic function at low temperatures. Given the unique peculiarities of the Antarctic environment, studying the features of Antarctic aspartic proteinases could provide new insights into the role of these proteins in fish physiology. We describe here: (1) the biochemical properties of a cathepsin D purified from the liver of the hemoglobinless icefish Chionodraco hamatus; (2) the biochemical characterization of Trematomus bernacchii pepsins variants A1 and A2 obtained by heterologous expression in bacteria; and (3) the identification of two closely related, novel aspartic proteinases from the liver of the two Antarctic fish species mentioned above. Overall, the results show that Notothenioidei aspartic proteinases display a number of characteristics that are remarkably different from those of mammalian aspartic proteinases, including high turnover number or high catalytic efficiency. We have named the newly identified aspartic proteinases "Nothepsins" and classified them relative to aspartic proteinases from other species.     

Revised 2.3 A structure of porcine pepsin: evidence for a flexible subdomain

A revised three-dimensional crystal structure of ethanol-inhibited porcine pepsin refined to an R-factor of 0.171 at 2.3 A resolution is presented and compared to the refined structures of the fungal aspartic proteinases: penicillopepsin, rhizopuspepsin, and endothiapepsin. Pepsin is composed of two nearly equal N and C domains related by an intra dyad. The overall polypeptide fold and active site structures are homologous for pepsin and the fungal enzymes. The weak inhibition of pepsin by ethanol can be explained by the presence of one or more ethanol molecules, in the vicinity of the active site carboxylates, which slightly alter the hydrogen-bonding network and which may compete with substrate binding in the active site. Structural superposition analysis showed that the N domains aligned better than the C-domains for pepsin and the fungal aspartic proteinases: 107-140 C alpha pairs aligned to 0.72-0.85 A rms for the N domains; 64-95 C alpha pairs aligned to 0.78-1.03 A rms for the C domains. The major structural difference between pepsin and the fungal enzymes concerns a newly described subdomain whose conformation varies markedly among these enzyme structures. The subdomain in pepsin comprises nearly 100 residues and is composed of two contiguous segments within the C domain (residues 192-212 and 223-299). the subdomain is connected, or "hinged," to a mixed beta-sheet that forms one of the structurally invariant, active site psi-loops. Relative subdomain displacements as large as a 21.0 degrees rotation and a 5.9 A translation were observed among the different enzymes. There is some suggestion in pepsin that the subdomain may be flexible and perhaps plays a structural role in mediating substrate binding, determining the substrate specificity, or in the activation of the zymogen.     

A theoretical study of torsional flexibility in the active site of aspartic proteinases: Implications for catalysis

We have performed ab initio Hartree-Fock self-consistent field calculations on the active site of endothiapepsin. The active site was modeled as a formic acid/formate anion moiety (representing the catalytic aspartates, Asp-32 and -215) and a bound water molecule. Residues Gly-34, Ser-35, Gly-217, and Thr-218, which all form hydrogen bonds to the active site, were modeled using formamide and methanol molecules. The water molecule, which is generally believed to function as the attacking nucleophile in catalysis, was allowed to bind to the active site in four distinct configurations. The geometry of each configuration was optimized using two basis sets (4-31G and 4-31G*). The results indicate that in the native enzyme the nucleophilic water is bound in a catalytically inert configuration. However, by rotating the carboxyl group of Asp-32 by about 90° the water molecule can be reorientated to attack the scissile bond of the substrate. A model of the bound enzyme-substrate complex was constructed from the crystal structure of a difluorostatone inhibitor complexed with endothiapepsin. This model suggests that the substrate itself initiates the reorientation of the nucleophilic water immediately prior to catalysis by forcing the carboxyl group of Asp-32 to rotate. The theoretical results predict that the active site of endothiapepsin undergoes a large distortion during substrate binding and this observation has been used to explain some of the kinetics results which have been reported for mutant aspartic proteinases.     

Differential localization of two acid proteinases in germinating barley (Hordeum vulgare) seed

A cathepsin D-like aspartic proteinase (EC 3.4.23) is abundant in ungerminated barley ( Hordeum vulgare ) seed while a 30 kDa cysteine endoproteinase (EC 3.4.22) is one of the proteinases synthesized de novo in the germinating seed. In this work, the localization of these two acid proteinases was studied at both the tissue and subcellular levels by immunomicroscopy. The results confirm that they have completely different functions. The aspartic proteinase was present in the ungerminated seed and, during germination, it appeared in all the living tissues of the grain, including the shoot and root. Contrary to previous suggestions, it was not observed in the starchy endosperm. By immunoblotting, the high molecular mass form of the enzyme (32 + 16 kDa) was found in all the living tissues, whereas the low molecular mass form (29 + 11 kDa) was not present in the shoot or root, indicating that the two enzyme forms have different physiological roles. The aspartic proteinase was localized first in the scutellar protein bodies of germinating seed, and later in the vacuoles which are formed by fusion of the protein bodies. In contrast to the aspartic proteinase, the expression of the 30 kDa cysteine proteinase began during the first germination day, and it was secreted into the starchy endosperm; first from the scutellum and later from the aleurone layer. It was not found in either shoots or roots. The 30 kDa cysteine proteinase was detected in the Golgi apparatus and in the putative secretory vesicles of the scutellar epithelium. These results suggest that the aspartic proteinase functions only in the living tissues of the grain, as opposed to the 30 kDa cysteine proteinase which is apparently one of the proteases initiating the hydrolysis of storage proteins in the starchy endosperm.     

Aspartic proteinase genes in the Brassicaceae Arabidopsis thaliana and Brassica napus

Active aspartic proteinase is isolated from Brassica napus seeds and the peptide sequence is used to generate primers for PCR. We present here cDNA and genomic clones for aspartic proteinases from the closely related Brassicaceae Arabidopsis thaliana and Brassica napus. The Arabidopsis cDNA represents a single gene, while Brassica has at least 4 genes. Like other plant aspartic proteases, the two Brassicaceae enzymes contain an extra protein domain of about 100 amino acids relative to the mammalian forms. The intron/exon arrangement in the Brassica genomic clone is significantly different from that in mammalian genes. As the proteinase is isolated from seeds, the same tissue where 2S albumins are processed, this implies expression of one of the aspartic proteinase genes there.     

X-ray crystallographic analysis of inhibition of endothiapepsin by cyclohexyl renin inhibitors.

The crystal structures of endothiapepsin, a fungal aspartic proteinase (EC 3.4.23.6), cocrystallized with two oligopeptide renin inhibitors, PD125967 and PD125754, have been determined at 2.0-A resolution and refined to R-factors of 0.143 and 0.153, respectively. These inhibitors, which are of the hydroxyethylene and statine types, respectively, possess a cyclohexylalanine side chain at P1 and have interesting functionalities at the P3 position which, until now, have not been subjected to crystallographic analysis. PD125967 has a bis(1-naphthylmethyl)acetyl residue at P3, and PD125754 possesses a hydroxyethylene analogue of the P3-P2 peptide bond for proteolytic stability. The structures reveal that the S3 pocket accommodates one naphthyl ring with conformational changes of the Asp 77 and Asp 114 side chains, the other naphthyl group residing in the S4 region. The P3-P2 hydroxyethylene analogue of PD125754 forms a hydrogen bond with the NH of Thr 219, thereby making the same interaction with the enzyme as the equivalent peptide groups of all inhibitors studied so far. The absence of side chains at the P2 and P1' positions of this inhibitor allows water molecules to occupy the respective pockets in the complex. The relative potencies of PD125967 and PD125754 for endothiapepsin are consistent with the changes in solvent-accessible area which take place on inhibitor binding.     

X-ray analyses of aspartic proteinases. IV. Structure and refinement at 2.2 A resolution of bovine chymosin

The structure of calf chymosin (EC 3.4.23.3), the aspartic proteinase from the gastric mucosa, was solved using the technique of molecular replacement. We describe the use of different search models based on distantly related fungal aspartic proteinases and investigate the effect of using only structurally conserved regions. The structure has been refined to a crystallographic R-factor of 17% at 2.2 A resolution with an estimated co-ordinate error of 0.21 A. In all, 136 water molecules have been located of which eight are internal. The structure of chymosin resembles that of pepsin and other aspartic proteinases. However, there is a considerable rearrangement of the active-site "flap" and, in particular, Tyr75 (pepsin numbering), which forms part of the specificity pockets S1 and S1'. This is probably a consequence of crystal packing. Electrostatic interactions on the edge of the substrate binding cleft appear to account for the restricted proteolysis of the natural substrate kappa-casein by chymosin. The local environment of invariant residues is examined, showing that structural constraints and side-chain hydrogen bonding can play an important role in the conservation of particular amino acids.     

High-resolution X-ray diffraction study of the complex between endothiapepsin and an oligopeptide inhibitor: the analysis of the inhibitor binding and description of the rigid body shift in the enzyme.

The conformation of the synthetic renin inhibitor CP-69,799, bound to the active site of the fungal aspartic proteinase endothiapepsin (EC 3.4.23.6), has been determined by X-ray diffraction at 1.8 A resolution and refined to the crystallographic R factor of 16%. CP-69,799 is an oligopeptide transition--state analogue inhibitor that contains a new dipeptide isostere at the P1-P1' position. This dipeptide isostere is a nitrogen analogue of the well-explored hydroxyethylene dipeptide isostere, wherein the tetrahedral P1' C alpha atom has been replaced by trigonal nitrogen. The inhibitor binds in the extended conformation, filling S4 to S3' pockets, with hydroxyl group of the P1 residue positioned symmetrically between the two catalytic aspartates of the enzyme. Interactions between the inhibitor and the enzyme include 12 hydrogen bonds and extensive van der Waals contacts in all the pockets, except for S3'. The crystal structure reveals a bifurcated orientation of the P2 histidine side chain and an interesting relative rotation of the P3 phenyl ring to accommodate the cyclohexyl side chain at P1. The binding of the inhibitor to the enzyme, while producing no large distortions in the enzyme active site cleft, results in small but significant change in the relative orientation of the two endothiapepsin domains. This structural change may represent the action effected by the proteinase as it distorts its substrate towards the transition state for proteolytic cleavage.     

Purification and properties of an aspartic protease from potato tuber that is inhibited by a basic chitinase

A protease was isolated from potato ( Solanum tuberosum L. cv. Huinkul) tuber disks after 24 h of aeration when proteolysis is markedly increased. Purification was performed by ammonium sulfate precipitation, ion exchange chromatography, and affinity chromatography. A size of 40 kDa was estimated by sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE) and gel filtration, it is monomeric and its properties are consistent with those of aspartic proteinases (EC 3.4.23): it had a pH optimum between 4 and 5 and it was inhibited by pepstatin. Partial homology with other plant aspartic proteinases was observed in two sequenced tryptic fragments. It binds to Sepharose-concanavalin A and can be eluted with α -methyl mannoside, indicating that it is possibly glycosylated. Unlike other aspartic proteinases from Solanaceae that degrade pathogenesis-related proteins, it is unable to cleave a basic chitinase from potato. Moreover, this aspartic protease is strongly inhibited by the basic chitinase; the 50% inhibition is obtained when the molar ratio approaches 1, the same as with pepstatin. The interaction between this aspartic protease and a new type of endogenous inhibitor may be an interesting starting point to study the regulation of these aspartic proteases during stress.     

Localization and characterization of hemoglobin-degrading aspartic proteinases from the malarial parasite Plasmodium falciparum.

Three hemoglobin-degrading proteinases were partially purified from food vacuoles isolated from trophozoite-stage forms of the malarial parasite Plasmodium falciparum. Two of the proteinases (M1 and M2) were solubilized by repeated sonication. The remaining proteinase (M3) was solubilized by treatment of the particulate fraction with taurocholic acid, suggesting that proteinase M3 is a membrane-bound proteinase whereas proteinases M1 and M2 are weakly associated with parasite membrane. The location of these proteinases suggests that they may participate in the digestion of host cytosolic protein. After partial purification, but not before, proteinases M1, M2 and M3 are highly sensitive to pepstatin, supporting their designation as aspartic proteinases. These aspartic proteinases show broad specificity for protein substrates. Native hemoglobin, acid denatured hemoglobin and oxidatively damaged hemoglobin are comparable substrates. Hemoglobin within the food vacuole was shown to be primarily native hemoglobin. Chemical modification studies indicate that these three aspartic proteinases have similar properties. The peptide maps from degradation of hemoglobin, however, suggest that aspartic proteinases M1, M2 and M3 are distinct proteinases.     

Structure and function of plant aspartic proteinases.

Aspartic proteinases of the A1 family are widely distributed among plant species and have been purified from a variety of tissues. They are most active at acidic pH, are specifically inhibited by pepstatin A and contain two aspartic residues indispensible for catalytic activity. The three-dimensional structure of two plant aspartic proteinases has been determined, sharing significant structural similarity with other known structures of mammalian aspartic proteinases. With a few exceptions, the majority of plant aspartic proteinases identified so far are synthesized with a prepro-domain and subsequently converted to mature two-chain enzymes. A characteristic feature of the majority of plant aspartic proteinase precursors is the presence of an extra protein domain of about 100 amino acids known as the plant-specific insert, which is highly similar both in sequence and structure to saposin-like proteins. This insert is usually removed during processing and is absent from the mature form of the enzyme. Its functions are still unclear but a role in the vacuolar targeting of the precursors has been proposed. The biological role of plant aspartic proteinases is also not completely established. Nevertheless, their involvement in protein processing or degradation under different conditions and in different stages of plant development suggests some functional specialization. Based on the recent findings on the diversity of A1 family members in Arabidopsis thaliana, new questions concerning novel structure-function relationships among plant aspartic proteinases are now starting to be addressed.     

X-ray, neutron and NMR studies of the catalytic mechanism of aspartic proteinases

Current proposals for the catalytic mechanism of aspartic proteinases are largely based on X-ray structures of bound oligopeptide inhibitors possessing non-hydrolysable analogues of the scissile peptide bond. Until recent years, the positions of protons on the catalytic aspartates and the ligand in these complexes had not been determined with certainty due to the inadequate resolution of these analyses. There has been much interest in locating the catalytic protons at the active site of aspartic proteinases since this has major implications for detailed understanding of the mechanism of action and the design of improved transition state mimics for therapeutic applications. In this review we discuss the results of studies which have shed light on the locations of protons at the catalytic centre. The first direct determination of the proton positions stemmed from neutron diffraction data collected from crystals of the fungal aspartic proteinase endothiapepsin bound to a transition state analogue (H261). The neutron structure of the complex at a resolution of 2.1 A provided evidence that Asp 215 is protonated and that Asp 32 is the negatively charged residue in the transition state complex. Atomic resolution X-ray studies of inhibitor complexes have corroborated this finding. A similar study of the native enzyme established that it, unexpectedly, has a dipeptide bound at the catalytic site which is consistent with classical reports of inhibition by short peptides and the ability of pepsins to catalyse transpeptidation reactions. Studies by NMR have confirmed the findings of low-barrier and single-well hydrogen bonds in the complexes with transition state analogues.     

Direct observation by X-ray analysis of the tetrahedral "intermediate" of aspartic proteinases.

We report the X-ray analysis at 2.0 A resolution for crystals of the aspartic proteinase endothiapepsin (EC 3.4.23.6) complexed with a potent difluorostatone-containing tripeptide renin inhibitor (CP-81,282). The scissile bond surrogate, an electrophilic ketone, is hydrated in the complex. The pro-(R) (statine-like) hydroxyl of the tetrahedral carbonyl hydrate is hydrogen-bonded to both active-site aspartates 32 and 215 in the position occupied by a water in the native enzyme. The second hydroxyl oxygen of the hydrate is hydrogen-bonded only to the outer oxygen of Asp 32. These experimental data provide a basis for a model of the tetrahedral intermediate in aspartic proteinase-mediated cleavage of the amide bond. This indicates a mechanism in which Asp 32 is the proton donor and Asp 215 carboxylate polarizes a bound water for nucleophilic attack. The mechanism involves a carboxylate (Asp 32) that is stabilized by extensive hydrogen bonding, rather than an oxyanion derivative of the peptide as in serine proteinase catalysis.     

Aspartic proteinase from barley grains is related to mammalian lysosomal cathepsin D

Resting barley (Hordeum vulgare L.) grains contain acid-proteinase activity. The corresponding enzyme was purified from grain extracts by affinity chromatography on a pepstatin-Sepharose column. The pH optimum of the affinity-purified enzyme was between 3.5 and 3.9 as measured by hemoglobin hydrolysis and the enzymatic activity was completely inhibited by pepstatin a specific inhibitor of aspartic proteinases (EC 3.4.23). Further purification on a Mono S column followed by activity measurements and sodium dodecyl sulfate-polyacrylamide gel electrophoresis revealed that the affinity-purified enzyme preparation contained two active heterodimeric aspartic proteinases: a larger 48k Da enzyme, consisting of 32-kDa and 16-kDa subunits and a smaller one of 40 kDa, consisting of 29-kDa and 11-kDa subunits. Separation and partial amino acid sequence analysis of each subunit indicate that the 40-kDa enzyme is formed by proteolytic processing of the 48k Da form. Amino-acid sequence alignment and inhibition studies showed that the barley aspartic proteinase resembles mammalian lysosomal cathepsin D (EC 3.4.23.5).