Sequence,Structural Analysis and Metrics to Define the Unique Dynamic Features of the Flap Regions Among Aspartic Proteases

Aspartic proteases are a class of hydrolytic enzymes that have been implicated in a number of diseases such as HIV, malaria, cancer and Alzheimer’s. The flap region of aspartic proteases is a characteristic unique structural feature of these enzymes; and found to have a profound impact on protein overall structure, function and dynamics. Flap dynamics also plays a crucial role in drug binding and drug resistance. Therefore, understanding the structure and dynamic behavior of this flap regions is crucial in the design of potent and selective inhibitors against aspartic proteases. Defining metrics that can describe the flap motion/dynamics has been a challenging topic in literature. This review is the first attempt to compile comprehensive information on sequence, structure, motion and metrics used to assess the dynamics of the flap region of different aspartic proteases in “one pot”. We believe that this review would be of critical importance to the researchers from different scientific domains.     

The primary structure of staphylococcal protease.

The amino acid sequence of staphylococcal protease has been determined by analysis of tryptic peptides obtained from cyanogen bromide fragments. Selected peptides obtained from digests with staphylococcal protease, thermolysin, and chymotrypsin provided the information necessary to align the tryptic peptides and the cyanogen bromide fragments. The protease is a single polypeptide chain of some 250 amino acids and is devoid of sulfhydryl groups. The COOH-terminal tryptic peptide of of the protease molecule contains some 43 residues, most of which are aspartic acids, asparagines, and prolines. The amino acid sequence of this peptide was not determined. The primary structure near the active serine residue indicates that staphylococcal protease is related to the pancreatic serine proteases. However, it has little or no additional sequence homologies with these enzymes except for the regions near histidine-50 and aspartic acid - 91. These regions have striking similarities with the corresponding regions of protease B and the trypsin-like enzyme of Streptomyces griseus.     

Evolution in the structure and function of aspartic proteases

Aspartic proteases (EC3.4.23) are a group of proteolytic enzymes of the pepsin family that share the same catalytic apparatus and usually function in acid solutions. This latter aspect limits the function of aspartic proteases to some specific locations in different organisms; thus the occurrence of aspartic proteases is less abundant than other groups of proteases, such as serine proteases. The best known sources of aspartic proteases are stomach (for pepsin, gastricsin, and chymosin), lysosomes (for cathepsins D and E), kidney (for renin), yeast granules, and fungi (for secreted proteases such as rhizopuspepsin, penicillopepsin, and endothiapepsin). These aspartic proteases have been extensively studied for their structure and function relationships and have been the topics of several reviews or monographs (Tang: Acid Proteases, Structure, Function and Biology. New York: Plenum Press, 1977; Tang: J Mol Cell Biochem 26:93-109, 1979; Kostka: Aspartic Proteinases and Their Inhibitors. Berlin: Walter de Gruyter, 1985). All mammalian aspartic proteases are synthesized as zymogens and are subsequently activated to active proteases. Although a zymogen for a fungal aspartic protease has not been found, the cDNA structure of rhizopuspepsin suggests the presence of a "pro" enzyme (Wong et al: Fed Proc 44:2725, 1985). It is probable that other fungal aspartic proteases are also synthesized as zymogens. It is the aim of this article to summarize the major models of structure-function relationships of aspartic proteases and their zymogens with emphasis on more recent findings. Attempts will also be made to relate these models to other aspartic proteases.     

Conserved interactions of the active carboxyls in pepsin-like enzymes and retroviral proteases

In addition to previous studies, 30 crystal structures of retroviral proteases corresponding to the highest resolution were inspected to analyze the interactions of the active carboxyl with surroundings groups. The outer oxygen of the active carboxyl in retroviral enzymes form contacts only with the water molecule participating in catalysis. This is an important difference between retroviral proteases and pepsin-like enzymes, which form a net of hydrogen bonds of these outer oxygen with residues neighboring the catalytic site in 3D structures. At the same time, it was found that in all aspartic proteases the inner oxygen of the active carboxyl are also involved in the chain of interactions through peptide groups Thr-Gly adjacent to the active residues. Polarization of these peptide groups influences the donor-acceptor properties of the active carboxyl. The found chain of interactions is more extensive in retroviral than in pepsin-like proteases; however, its main part is conserved for the whole class of these enzymes. Some implications of the role of these interactions are discussed.     

A dual target of Plasmepsin IX and X: Unveiling the atomistic superiority of a core chemical scaffold in malaria therapy

Plasmepsin IX and X, members of the prominent aspartic family of proteases whose function were hitherto unknown have only recently been established as key mediators of erythrocyte invasion and egress of the virulent malarial parasite. Inhibitor 49c, a potent antimalarial peptidomimetic inhibitor initially developed to target Plasmepsin II has lately been proven to exhibit potent inhibitory activity against Plasmepsin IX and X. However, the molecular and structural dynamics supporting its inhibitory activity remain inconclusive. Hindering the motion of the flap and hinge region of an aspartic protease remains essential for disabling the catalytic activity of the enzyme. Integrating molecular dynamic simulations coupled with other advanced biocomputational tools, we reveal the enhanced structural mechanistic competence of 49c in complex with Plasmepsin IX and X relative to Pepstatin. Pepstatin, a known aspartic protease inhibitor which actively hinders the opening and closing of the flap tip and flexible loop and consequently limits access to the catalytic aspartic residues, however, its administration has been related to elevated levels of toxicity. Thermodynamic calculations reveal a higher relative binding free energy associated with Plasmepsin IX and X in complex with 49c as opposed to Pepstatin. A relatively compact and structurally rigid 49c bound complexes sequel into the restriction of the flap and hinge residues by restraining cohesive movement, consequently hindering their “twisting motion” from transpiring. Findings unveil an atomistic perspective into the structural superiority of 49c in complex with Plasmepsin IX and X.     

Conserved Interactions of the Active Carboxyls in Pepsin-like Enzymes and Retroviral Proteases

In addition to previous studies, 30 crystal structures of retroviral proteases corresponding to the highest resolution were inspected to analyze the interactions of the active carboxyls with surroundings groups. The outer oxygens of the active carboxyls in retroviral enzymes form contacts only with the water molecule participating in catalysis. This is an important difference between retroviral proteases and pepsin-like enzymes, which form a net of hydrogen bonds of these outer oxygens with residues neighboring the catalytic site in 3D structures. At the same time, it was found that in all aspartic proteases the inner oxygens of the active carboxyls are also involved in the chain of interactions through peptide groups Thr–Gly adjacent to the active residues. Polarization of these peptide groups influences the donor–acceptor properties of the active carboxyls. The found chain of interactions is more extensive in retroviral than in pepsin-like proteases; however, its main part is conserved for the whole class of these enzymes. Some implications of the role of these interactions are discussed.     

Dynamic properties of the native free antithrombin from molecular dynamics simulations: computational evidence for solvent- exposed Arg393 side chain

While antithrombin (AT) has small basal inhibitory activity, it reaches its full inhibitory potential against activated blood coagulation factors, FXa, FIXa, and FIIa (thrombin), via an allosteric and/or template (bridging) mechanism by the action of heparin, heparan sulfate, or heparin-mimetic pentasaccharides (PS). From the numerous X-ray structures available for different conformational states of AT, only indirect and incomplete conclusions can be drawn on the inherently dynamic properties of AT. As a typical example, the basal inhibitory activity of AT cannot be interpreted on the basis of “non-activated” free antithrombin X-ray structures since the Arg393 side chain, playing crucial role in antithrombin-proteinase interaction, is not exposed. In order to reveal the intrinsic dynamic properties and the reason of basal inhibitory activity of antithrombin, 2 μs molecular dynamics simulations were carried out on its native free-forms. It was shown from the simulation trajectories that the reactive center loop which is functioning as “bait” for proteases, even without any biasing potential can populate conformational state in which the Arg393 side chain is solvent exposed. It is revealed from the trajectory analysis that the peptide sequences correspond to the helix D extension, and new helix P formation can be featured with especially large root-mean-square fluctuations. Mutual information analyses of the trajectory showed remarkable (generalized) correlation between those regions of antithrombin which changed their conformations as the consequence of AT–PS complex formation. This suggests that allosteric information propagation pathways are present even in the non-activated native form of AT.     

The protonation state of the catalytic aspartates in plasmepsin II

Assigning the correct protonation state to the catalytic residues is essential for a realistic modelling of an enzyme's active site. Plasmepsins are pharmaceutically relevant aspartic proteases involved in haemoglobin degradation by Plasmodium spp. In aspartic proteases, one of the two catalytic aspartates is protonated, while the other is negatively charged. Here, multiple explicit-water molecular dynamics simulations of plasmepsin II, uncomplexed and with a hydroxypropylamine peptidomimetic inhibitor, indicate that protonation of Asp214 favours a stable active site structure. Moreover, the protonation state of the catalytic aspartate has a strong influence on a linear chain of hydrogen bonds with the adjacent side chains.     

Homology cloning of aspartic proteases from an endocrine cell line using the polymerase chain reaction

Oligonucleotides directed towards the active site regions of aspartic proteases were used as primers for the polymerase chain reaction to identify a unique sequence (asppcr1) from the AtT-20 anterior pituitary corticotrope cell line. Asppcr1 showed the greatest similarity (85% identity) to human cathepsin E [(1989) J. Biol. Chem. 264, 16748-16753]. Northern blot analysis of AtT-20 RNA revealed a single 1.9 kB message. Nuclease protection experiments indicated that asppcr1 mRNA was present in pancreas, spleen, testis and liver at low levels and undetectable in heart and brain. This contrasted with the lysosomal aspartic protease, cathepsin D whose mRNA showed a broader tissue distribution. The restricted message distribution of asppcr1 supports a more specific role for this aspartic protease in aspect(s) of cellular physiology.     

Widespread tissue expression of nepenthesin-like aspartic protease genes in Arabidopsis thaliana.

There are approximately 69 genes encoding aspartyl protease homologues in Arabidopsis thaliana, and most of the gene products constitute a novel subfamily of aspartic proteases. However, their physiological roles are largely unknown. As an initial step to shed light on the roles of these nepenthesin-like aspartic proteases (NAPs), a phylogenetic tree was constructed, which indicated that these proteases are classified into several distinct sub-sub-groups. Based on these results, specific primers were designed for genes selected from several of these groups and their tissue expression was investigated using RT-PCR. The results indicated that these genes are widely expressed in several tissues, such as leaves, stems, seeds and pods, suggesting ubiquitous occurrence and multiple functions of the corresponding proteases in the tissues of A. thaliana.     

Crystal structure of an activation intermediate of cathepsin E

Cathepsin E is an intracellular, non-lysosomal aspartic protease expressed in a variety of cells and tissues. The protease has proposed physiological roles in antigen presentation by the MHC class II system, in the biogenesis of the vasoconstrictor peptide endothelin, and in neurodegeneration associated with brain ischemia and aging. Cathepsin E is the only A1 aspartic protease that exists as a homodimer with a disulfide bridge linking the two monomers. Like many other aspartic proteases, it is synthesized as a zymogen which is catalytically inactive towards its natural substrates at neutral pH and which auto-activates in an acidic environment. Here we report the crystal structure of an activation intermediate of human cathepsin E at 2.35A resolution. The overall structure follows the general fold of aspartic proteases of the A1 family, and the intermediate shares many features with the intermediate 2 on the proposed activation pathway of aspartic proteases like pepsin C and cathepsin D. The pro-sequence is cleaved from the protease and remains stably associated with the mature enzyme by forming the outermost sixth strand of the interdomain beta-sheet. However, different from these other aspartic proteases the pro-sequence of cathepsin E remains intact after cleavage from the mature enzyme. In addition, the active site of cathepsin E in the crystal is occupied by N-terminal amino acid residues of the mature protease in the non-primed binding site and by an artificial N-terminal extension of the pro-sequence from a neighboring molecule in the primed site. The crystal structure of the cathepsin E/pro-sequence complex, therefore, provides further insight into the activation mechanism of aspartic proteases.     

Aspartic proteases of Plasmodium falciparum and other parasitic protozoa as drug targets

All parasitic protozoa contain multiple proteases, some of which are attracting attention as drug targets. Aspartic proteases are already the targets of some clinically useful drugs (e.g. chemotherapy of HIV infection) and a variety of factors make these enzymes appealing to those seeking novel antiparasite therapies. This review provides a critical analysis of the current knowledge on Plasmodium aspartic proteases termed plasmepsins, proposes a definitive nomenclature for this group of enzymes, and compares these enzymes with aspartic proteases of humans and other parasitic protozoa. The present status of attempts to obtain specific inhibitors of the parasite enzymes that will be useful as drugs is outlined and suggestions for future research priorities are proposed.     

Aspartic proteases of Plasmodium vivax are highly conserved in wild isolates

The plasmepsins are the aspartic proteases of malaria parasites. Treatment of aspartic protease inhibitor inhibits hemoglobin hydrolysis and blocks the parasite development in vitro suggesting that these proteases might be exploited their potentials as antimalarial drug targets. In this study, we determined the genetic variations of the aspartic proteases of Plasmodium vivax (PvPMs) of wild isolates. Two plasmepsins (PvPM4 and PvPM5) were cloned and sequenced from 20 P. vivax Korean isolates and two imported isolates. The sequences of the enzymes were highly conserved except a small number of amino acid substitutions did not modify key residues for the function or the structure of the enzymes. The high sequence conservations between the plasmepsins from the isolates support the notion that the enzymes could be reliable targets for new antimalarial chemotherapeutics.     

Microbial aspartic proteases: current and potential applications in industry

Aspartic proteases are a relatively small group of proteolytic enzymes that are active in acidic environments and are found across all forms of life. Certain microorganisms secrete such proteases as virulence agents and/or in order to break down proteins thereby liberating assimilable sources of nitrogen. Some of the earlier applications of these proteolytic enzymes are found in the manufacturing of cheese where they are used as milk-clotting agents. Over the last decade, they have received tremendous research interest because of their involvement in human diseases. Furthermore, there has also been a growing interest on these enzymes for their applications in several other industries. Recent research suggests in particular that they could be used in the wine industry to prevent the formation of protein haze while preserving the wines’ organoleptic properties. In this mini-review, the properties and mechanisms of action of aspartic proteases are summarized. Thereafter, a brief overview of the industrial applications of this specific class of proteases is provided. The use of aspartic proteases as alternatives to clarifying agents in various beverage industries is mentioned, and the potential applications in the wine industry are thoroughly discussed.     

Evaluation of the presence of aspartic proteases from Centaurea calcitrapa during seed germination

Aspartic proteinases are present in a variety of organisms including plants. Common features of aspartic proteases include an active site cleft that contains two catalytic aspartic residues, acid pH optima for enzymatic activity, inhibition by pepstatin A. Plant aspartic proteinases occur in seeds and may be involved in the processing of storage proteins. Many of them have been purified and characterized. The presence of aspartic proteases in seeds of Centaurea calcitrapa during germination was investigated by measuring the activity on enzyme extracts. The aspartic proteases are present mainly in the beginning of seed germination suggesting that they could initiate the degradation of protein reserves in germinating seeds.

These proteases were purified by salt precipitation followed by anion-exchange chromatography. Purified aspartic proteases have an optimal pH between 3.5 and 4.5, using FTC-hemoglobin as substrate and an optimal temperature at 52 °C. The ability of seed extracts for milk clotting was tested and the clotting time that was achieved is in the same range found for flower extracts appropriated for special cheeses in which weak clotting agents are required.     

Conformational homogeneity in molecular recognition by proteolytic enzymes

Crystal structures for several hundred protease-inhibitor complexes have been analysed and their superimpositions have been used to demonstrate a universal relationship between inhibitor/substrate conformation and molecular recognition by all aspartic, serine, cysteine and metallo proteases. Proteases universally recognize an extended beta strand conformation in all their peptidic (and non-peptidic) inhibitors and substrate analogues without significant exceptions. This conformational homogeneity is illustrated here for a subset of 180 protease-inhibitor structures which are displayed as (a) structural overlays of multiple inhibitors for each of eight aspartic, eight serine, six metallo and five cysteine proteases; (b) single inhibitors each bound to different proteases; and (c) Ramachandran plots of peptide or pseudo-peptide dihedral angle pairs which demonstrate beta strands (Phi -54 degrees to -173 degrees, Psi 24 degrees to 174 degrees ) like those normally found paired in proteins as beta sheets. However, unlike beta sheets, alpha and 3(10) helices, beta and gamma turns, where the folded main chain amide components are intramolecularly hydrogen bonded and thus unavailable for interaction with proteins, an inhibitor/substrate in an isolated beta strand conformation provides maximum exposure of its hydrogen bonding donors/acceptors and side chain components to a putative protease receptor. This analysis highlights the advantages of a strand conformation over other elements of secondary structure for protease recognition and may lead to generic strategies for inhibitor design.     

The S8 serine, C1A cysteine and A1 aspartic protease families in Arabidopsis

The Arabidopsis thaliana genome has over 550 protease sequences representing all five catalytic types: serine, cysteine, aspartic acid, metallo and threonine (MEROPS peptidase database, http://merops.sanger.ac.uk/), which probably reflect a wide variety of as yet unidentified functions performed by plant proteases. Recent indications that the 26S proteasome, a T1 family-threonine protease, is a regulator of light and hormone responsive signal transduction highlight the potential of proteases to participate in many aspects of plant growth and development. Recent discoveries that proteases are required for stomatal distribution, embryo development and disease resistance point to wider roles for four additional multigene families that include some of the most frequently studied (yet poorly understood) plant proteases: the subtilisin-like, serine proteases (family S8), the papain-like, cysteine proteases (family C1A), the pepsin-like, aspartic proteases (family A1) and the plant matrixin, metalloproteases (family M10A). In this report, 54 subtilisin-like, 30 papain-like and 59 pepsin-like proteases from Arabidopsis, are compared with S8, C1A and A1 proteases known from other plant species at the functional, phylogenetic and gene structure levels. Examples of structural conservation between S8, C1A and A1 genes from rice, barley, tomato and soybean and those from Arabidopsis are noted, indicating that some common, essential plant protease roles were established before the divergence of monocots and eudicots. Numerous examples of tandem duplications of protease genes and evidence for a variety of restricted expression patterns suggest that a high degree of specialization exists among proteases within each family. We propose that comprehensive analysis of the functions of these genes in Arabidopsis will firmly establish serine, cysteine and aspartic proteases as regulators and effectors of a wide range of plant processes.     

Building self-avoiding lattice models of proteins using a self-consistent field optimization

We present an algorithm to build self-avoiding lattice models of chain molecules with low RMS deviation from their actual 3D structures. To find the optimal coordinates for the lattice chain model, we minimize a function that consists of three terms: (1) the sum of squared deviations of link coordinates on a lattice from their off-lattice values, (2) the sum of “short-range” terms, penalizing violation of chain connectivity, and (3) the sum of “long-range” repulsive terms, penalizing chain self-intersections. We treat this function as a chain molecule “energy” and minimize it using self-consistent field (SCF) theory to represent the pairwise link repulsions as 3D fields acting on the links. The statistical mechanics of chain molecules enables computation of the chain distribution in this field on the lattice. The field is refined by iteration to become self-consistent with the chain distribution, then dynamic programming is used to find the optimal lattice model as the “lowest-energy” chain pathway in this SCF. We have tested the method on one of the coarsest (and most difficult) lattices used for model building on proteins of all structural types and show that the method is adequate for building self-avoiding models of proteins with low RMS deviations from the actual structures. © 1996 Wiley-Liss, Inc.     

Biosynthesis of lysosomal endopeptidases

Despite the clear differences between the amino acid sequence and enzymatic specificity of aspartic and cysteine endopeptidases, the biosynthetic processing of lysosomal members of these two families is very similar. With in vitro translation and pulse-chase analysis in tissue culture cells, the biosynthesis of cathepsin D, a aspartic protease, and cathepsins B, H and L, cysteine proteases, are compared. Both aspartic and cysteine endopeptidases undergo cotranslational cleavage of an amino-terminal signal peptide that mediates transport across the endoplasmic reticulum (ER) membrane. Addition of high-mannose carbohydrate also occurs cotranslationally in the lumen of the ER. Proteases of both enzyme classes are initially synthesized as inactive proenzymes possessing amino-terminal activation peptides. Removal of the propeptide generates an active single-chain enzyme. Whether the single-chain enzyme undergoes asymmetric cleavage into a light and a heavy chain appears to be cell type specific. Finally, late during their biosynthesis both classes of enzymes undergo amino acid trimming, losing a few amino acid residues at the cleavage site between the light and heavy chains and/or at their carboxyltermini. During biosynthesis these enzymes are also secreted to some extent. In most cells the secreted enzyme is the proenzyme bearing some complex carbohydrate. Under certain physiological conditions the inactive secreted enzymes may become activated as a result of a conformational change that may or may not result in autolysis. Analysis of the biochemical nature of the various processing steps helps define the cellular pathway followed by newly synthesized proteases targeted to the lysosome.