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

Extracellular proteolysis in cancer: Proteases, substrates, and mechanisms in tumor progression and metastasis

A vast ensemble of extracellular proteins influences the development and progression of cancer, shaped and reshaped by a complex network of extracellular proteases. These proteases, belonging to the distinct classes of metalloproteases, serine proteases, cysteine proteases, and aspartic proteases, play a critical role in cancer. They often become dysregulated in cancer, with increases in pathological protease activity frequently driven by the loss of normal latency controls, diminished regulation by endogenous protease inhibitors, and changes in localization. Dysregulated proteases accelerate tumor progression and metastasis by degrading protein barriers within the extracellular matrix (ECM), stimulating tumor growth, reactivating dormant tumor cells, facilitating tumor cell escape from immune surveillance, and shifting stromal cells toward cancer-promoting behaviors through the precise proteolysis of specific substrates to alter their functions. These crucial substrates include ECM proteins and proteoglycans, soluble proteins secreted by tumor and stromal cells, and extracellular domains of cell surface proteins, including membrane receptors and adhesion proteins. The complexity of the extracellular protease web presents a significant challenge to untangle. Nevertheless, technological strides in proteomics, chemical biology, and the development of new probes and reagents are enabling progress and advancing our understanding of the pivotal importance of extracellular proteolysis in cancer.

Discovery of new Schistosoma mansoni aspartyl protease inhibitors by structure-based virtual screening

BACKGROUND

Schistosomiasis is a neglected tropical disease caused by trematodes of the genus Schistosoma, with a limited treatment, mainly based on the use of praziquantel (PZQ). Currently, several aspartic proteases genes have already been identified within the genome of Schistosoma species. At least one enzyme encoded from this gene family (SmAP), named SmCD1, has been validated for the development of schistosomicidal drugs, since it has a key role in haemoglobin digestion by worms.

OBJECTIVE

In this work, we integrated a structure-based virtual screening campaign, enzymatic assays and adult worms ex vivo experiments aiming to discover the first classes of SmCD1 inhibitors.

METHODS

Initially, the 3D-structures of SmCD1, SmCD2 and SmCD3 were generated using homology modelling approach. Using these models, we prioritised 50 compounds from 20,000 compounds from ChemBridge database for further testing in adult worm aqueous extract (AWAE) and recombinant SmCD1 using enzymatic assays.

FINDINGS

Seven compounds were confirmed as hits and among them, two compounds representing new chemical scaffolds, named 5 and 19, had IC50 values against SmCD1 close to 100 μM while presenting binding efficiency indexes comparable to or even higher than pepstatin, a classical tight-binding peptide inhibitor of aspartyl proteases. Upon activity comparison against mammalian enzymes, compound 50 was selective and the most potent against the AWAE aspartic protease activity (IC50 = 77.7 μM). Combination of computational and experimental results indicate that compound 50 is a selective inhibitor of SmCD2. Compounds 5, 19 and 50 tested at low concentrations (10 uM) were neither cytotoxic against WSS-1 cells (48 h) nor could kill adult worms ex-vivo, although compounds 5 and 50 presented a slight decrease on female worms motility on late incubations times (48 or 72 h).

MAIN CONCLUSION

Overall, the inhibitors identified in this work represent promising hits for further hit-to-lead optimisation.

Elucidating the multiple roles of hydration for accurate protein-ligand binding prediction via deep learning

Accurate and efficient prediction of protein-ligand interactions has been a long-lasting dream of practitioners in drug discovery. The insufficient treatment of hydration is widely recognized to be a major limitation for accurate protein-ligand scoring. Using an integration of molecular dynamics simulations on thousands of protein structures with novel big-data analytics based on convolutional neural networks and deep Taylor decomposition, we consistently identify here three different patterns of hydration to be essential for protein-ligand interactions. In addition to desolvation and water-mediated interactions, the formation of enthalpically favorable networks of first-shell water molecules around solvent-exposed ligand moieties is identified to be essential for protein-ligand binding. Despite being currently neglected in drug discovery, this hydration phenomenon could lead to new avenues in optimizing the free energy of ligand binding. Application of deep neural networks incorporating hydration to docking provides 89% accuracy in binding pose ranking, an essential step for rational structure-based drug design.

Active Site Mapping of an Aspartic Protease by Multiple Fragment Crystal Structures: Versatile Warheads To Address a Catalytic Dyad

Crystallography is frequently used as follow-up method to validate hits identified by biophysical screening cascades. The capacity of crystallography to directly screen fragment libraries is often underestimated, due to its supposed low-throughput and need for high-quality crystals. We applied crystallographic fragment screening to map the protein-binding site of the aspartic protease endothiapepsin by individual soaking experiments. Here, we report on 41 fragments binding to the catalytic dyad and adjacent specificity pockets. The analysis identifies already known warheads but also reveals hydrazide, pyrazole, or carboxylic acid fragments as novel functional groups binding to the dyad. A remarkable swapping of the S1 and S1' pocket between structurally related fragments is explained by either steric demand, required displacement of a well-bound water molecule, or changes of trigonal-planar to tetrahedral geometry of an oxygen functional group in a side chain. Some warheads simultaneously occupying both S1 and S1' are promising starting points for fragment-growing strategies.

Allosteric-Activation Mechanism of Bovine Chymosin Revealed by Bias-Exchange Metadynamics and Molecular Dynamics Simulations

The aspartic protease, bovine chymosin, catalyzes the proteolysis of κ-casein proteins in milk. The bovine chymosin-κ-casein complex is of industrial interest as the enzyme is used extensively in the manufacturing of processed dairy products. The apo form of the enzyme adopts a self-inhibited conformation in which the side chain of Tyr77 occludes the binding site. On the basis of kinetic, mutagenesis, and crystallographic data, it has been widely reported that a HPHPH sequence in the P8-P4 residues of the natural substrate κ-casein acts as the allosteric activator, but the mechanism by which this occurs has not previously been elucidated due to the challenges associated with studying this process by experimental methods. Here we have employed two computational techniques, molecular dynamics and bias-exchange metadynamics simulations, to study the mechanism of allosteric activation and to compute the free energy surface for the process. The simulations reveal that allosteric activation is initiated by interactions between the HPHPH sequence of κ-casein and a small α-helical region of chymosin (residues 112-116). A small conformational change in the α-helix causes the side chain of Phe114 to vacate a pocket that may then be occupied by the side chain of Tyr77. The free energy surface for the self-inhibited to open transition is significantly altered by the presence of the HPHPH sequence of κ-casein.

Secreted fungal aspartic proteases: A review

The aspartic proteases, also called aspartyl and aspartate proteases or acid proteases (E.C.3.4.23), belong to the endopeptidase family and are characterized by the conserved sequence Asp-Gly-Thr at the active site. These enzymes are found in a wide variety of microorganisms in which they perform important functions related to nutrition and pathogenesis. In addition, their high activity and stability at acid pH make them attractive for industrial application in the food industry; specifically, they are used as milk-coagulating agents in cheese production or serve to improve the taste of some foods. This review presents an analysis of the characteristics and properties of secreted microbial aspartic proteases and their potential for commercial application.

Crystallization of nepenthesin I using a low-pH crystallization screen

Nepenthesins are aspartic proteases secreted by carnivorous pitcher plants of the genus Nepenthes. They significantly differ in sequence from other plant aspartic proteases. This difference, which provides more cysteine residues in the structure of nepenthesins, may contribute to their unique stability profile. Recombinantly produced nepenthesin 1 (rNep1) from N. gracilis in complex with pepstatin A was crystallized under two different crystallization conditions using a newly formulated low-pH crystallization screen. The diffraction data were processed to 2.9 and 2.8 Å resolution, respectively. The crystals belonged to space group P212121, with unit-cell parameters a = 86.63, b = 95.90, c = 105.40 Å, α = β = γ = 90° and a = 86.28, b = 97.22, c = 103.78 Å, α = β = γ = 90°, respectively. Matthews coefficient and solvent-content calculations suggest the presence of two molecules of rNep1 in the asymmetric unit. Here, the details of the crystallization experiment and analysis of the X-ray data are reported.

Aspartic protease inhibitors containing tertiary alcohol transition-state mimics

Aspartic proteases (APs) are a class of enzymes engaged in the proteolytic digestion of peptide substrates. APs play important roles in physiological and infectious pathways, making them plausible drug targets. For instance in the treatment of HIV infections, access to an efficient combination of protease and reverse transcriptase inhibitors have changed a terminal illness to a chronic but manageable disease. However, the benefits have been limited due to the emergence of drug resistant viral strains, poor pharmacokinetic properties of peptidomimetic inhibitors and adverse effects associated with the treatment. In the 1980s, D. Rich and co-workers proposed a novel strategy for the development of AP inhibitors by replacing the secondary hydroxyl group with a tertiary alcohol as part of the transition state (TS) mimicking moiety. This strategy has been extensively explored over the last decade with a common belief that masking of the polar group, e.g. by intramolecular hydrogen bonding, has the potential to enhance transcellular transport. This is the first review presenting the advances of AP inhibitors comprising a tertiary hydroxyl group. The inhibitors have been classified into different tert-hydroxy TS mimics and their design strategies, synthesis, biological activities, structure-activity-relationships and X-ray structures are discussed.

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.

Interactions of a low molecular weight inhibitor from Streptomyces sp. MBR04 with human cathepsin D: implications in mechanism of inactivation

Cathepsin D, a lysosomal aspartic protease, is of potential interest as a target for drug design due to its implication in breast and ovarian cancer. The article reports a low molecular weight cathepsin D inhibitor from Streptomyces sp. MBR04. The Mr of the inhibitor was 1,078 Da as determined by MALDI-TOF, and the amino acid analysis showed the presence of Asp, Asp, Gly, Ala, Lys, Leu, Tyr, Trp residues. The steady-state kinetic interactions revealed reversible, competitive, slow-tight-binding nature of the inhibitor with an IC50 and K i values of 3.2 and 2.5 nM, respectively. The binding of the inhibitor with the enzyme and the subsequent conformational changes were monitored by exploiting the intrinsic fluorescence of the surface exposed Trp-54 residue. Based on the fluorescence and circular dichroism studies, we demonstrate that the inhibitor binds to the active site of cathepsin D and causes inactivation. All these kinetic, thermodynamic, and quenching studies suggest that the newly isolated peptidic inhibitor could be a potential scaffold to study and can be used to develop new potent therapeutic lead molecule for the development of drugs. The inhibitor will be significant as a potential lead molecule to target cathepsin D.

Aspartic peptidases of human pathogenic trypanosomatids: perspectives and trends for chemotherapy

Aspartic peptidases are proteolytic enzymes present in many organisms like vertebrates, plants, fungi, protozoa and in some retroviruses such as human immunodeficiency virus (HIV). These enzymes are involved in important metabolic processes in microorganisms/virus and play major roles in infectious diseases. Although few studies have been performed in order to identify and characterize aspartic peptidase in trypanosomatids, which include the etiologic agents of leishmaniasis, Chagas’ disease and sleeping sickness, some beneficial properties of aspartic peptidase inhibitors have been described on fundamental biological events of these pathogenic agents. In this context, aspartic peptidase inhibitors (PIs) used in the current chemotherapy against HIV (e.g., amprenavir, indinavir, lopinavir, nelfinavir, ritonavir and saquinavir) were able to inhibit the aspartic peptidase activity produced by different species of Leishmania. Moreover, the treatment of Leishmania promastigotes with HIV PIs induced several perturbations on the parasite homeostasis, including loss of the motility and arrest of proliferation/growth. The HIV PIs also induced an increase in the level of reactive oxygen species and the appearance of irreversible morphological alterations, triggering parasite death pathways such as programed cell death (apoptosis) and uncontrolled autophagy. The blockage of physiological parasite events as well as the induction of death pathways culminated in its incapacity to adhere, survive and escape of phagocytic cells. Collectively, these results support the data showing that parasites treated with HIV PIs have a significant reduction in the ability to cause in vivo infection. Similarly, the treatment of Trypanosoma cruzi cells with pepstatin A showed a significant inhibition on both aspartic peptidase activity and growth as well as promoted several and irreversible morphological changes. These studies indicate that aspartic peptidases can be promising targets in trypanosomatid cells and aspartic proteolytic inhibitors can be benefic chemotherapeutic agents against these human pathogenic microorganisms.

The structure of endothiapepsin complexed with a Phe-Tyr reduced-bond inhibitor at 1.35 Å resolution

Endothiapepsin is a typical member of the aspartic proteinase family. The catalytic mechanism of this family is attributed to two conserved catalytic aspartate residues, which coordinate the hydrolysis of a peptide bond. An oligopeptide inhibitor (IC50 = 0.62 µM) based on a reduced-bond transition-state inhibitor of mucorpepsin was co-crystallized with endothiapepsin and the crystal structure of the enzyme-inhibitor complex was determined at 1.35 Å resolution. A total of 12 hydrogen bonds between the inhibitor and the active-site residues were identified. The resulting structure demonstrates a number of novel subsite interactions in the active-site cleft.

A high performance cloud-based protein-ligand docking prediction algorithm

The potential of predicting druggability for a particular disease by integrating biological and computer science technologies has witnessed success in recent years. Although the computer science technologies can be used to reduce the costs of the pharmaceutical research, the computation time of the structure-based protein-ligand docking prediction is still unsatisfied until now. Hence, in this paper, a novel docking prediction algorithm, named fast cloud-based protein-ligand docking prediction algorithm (FCPLDPA), is presented to accelerate the docking prediction algorithm. The proposed algorithm works by leveraging two high-performance operators: (1) the novel migration (information exchange) operator is designed specially for cloud-based environments to reduce the computation time; (2) the efficient operator is aimed at filtering out the worst search directions. Our simulation results illustrate that the proposed method outperforms the other docking algorithms compared in this paper in terms of both the computation time and the quality of the end result.

Capturing the reaction pathway in near-atomic-resolution crystal structures of HIV-1 protease

Snapshots of three consecutive steps in the proteolytic reaction of HIV-1 protease (PR) were obtained in crystal structures at resolutions of 1.2-1.4 Å. Structures of wild-type protease and two mutants (PR(V32I) and PR(I47V)) with V32I and I47V substitutions, which are common in drug resistance, reveal the gem-diol tetrahedral intermediate, the separating N- and C-terminal products, and the C-terminal product of an autoproteolytic peptide. These structures represent three stages in the reaction pathway and shed light on the reaction mechanism. The near-atomic-resolution geometric details include a short hydrogen bond between the intermediate and the outer carboxylate oxygen of one catalytic Asp25 that is conserved in all three structures. The two products in the complex with mutant PR(I47V) have a 2.2 Å separation of the amide and carboxyl carbon of the adjacent ends, suggesting partial cleavage prior to product release. The complex of mutant PR(V32I) with a single C-terminal product shows density for water molecules in the other half of the binding site, including a partial occupancy water molecule interacting with the product carboxylate end and the carbonyl oxygen of one conformation of Gly27, which suggests a potential role of Gly27 in recycling from the product complex to the ligand-free enzyme. These structural details at near-atomic resolution enhance our understanding of the reaction pathway and will assist in the design of mechanism-based inhibitors as antiviral agents.

Covalent inhibitors of fatty acid amide hydrolase: a rationale for the activity of piperidine and piperazine aryl ureas

Recently, covalent drugs have attracted great interest in the drug discovery community, with successful examples that have demonstrated their therapeutic effects. Here, we focus on the covalent inhibition of the fatty acid amide hydrolase (FAAH), which is a promising strategy in the treatment of pain and inflammation. Among the most recent and potent FAAH inhibitors (FAAHi), there are the cyclic piperidine and piperazine aryl ureas. FAAH hydrolyzes efficiently the amide bond of these compounds, forming a covalent enzyme-inhibitor adduct. To rationalize this experimental evidence, we performed an extensive computational analysis centered on piperidine-based PF750 (1) and piperazine-based JNJ1661010 (2), two potent lead compounds used to generate covalent inhibitors as clinical candidates. We found that FAAH induces a distortion of the amide bond of the piperidine and piperazine aryl ureas. Quantum mechanics/molecular mechanics ΔE(LUMO-HOMO) energies indicate that the observed enzyme-induced distortion of the amide bond favors the formation of a covalent FAAH-inhibitor adduct. These findings could help in the rational structure-based design of novel covalent FAAHi.

Initial stage of cheese production: a molecular modeling study of bovine and camel chymosin complexed with peptides from the chymosin-sensitive region of κ-casein

Bovine chymosin has long been the preferred enzyme used to coagulate cow's milk, in the initial stage of cheese production, during which it cleaves a specific bond in the milk protein κ-casein. Recently, camel chymosin has been shown to have a 70% higher clotting activity toward cow's milk and, moreover, to cleave κ-casein more selectively. Bovine chymosin, on the other hand, is a poor clotting agent toward camel's milk. This paper reports a molecular modeling study aimed at understanding this disparity, based on homology modeling and molecular dynamics simulations using peptide fragments of κ-casein from cow and camel in both bovine and camel chymosin. The results show that the complex between bovine chymosin and the fragment of camel κ-casein is indeed less stable in the binding pocket. The results also indicate that this in part may be due to charge repulsion between a lysine residue in bovine chymosin and an arginine residue in the P4 position of camel κ-casein.

Bovine chymosin: a computational study of recognition and binding of bovine kappa-casein

Bovine chymosin is an aspartic protease that selectively cleaves the milk protein kappa-casein. The enzyme is widely used to promote milk clotting in cheese manufacturing. We have developed models of residues 97-112 of bovine kappa-casein complexed with bovine chymosin, using ligand docking, conformational search algorithms, and molecular dynamics simulations. In agreement with limited experimental evidence, the model suggests that the substrate binds in an extended conformation with charged residues on either side of the scissile bond playing an important role in stabilizing the binding pose. Lys111 and Lys112 are observed to bind to the N-terminal domain of chymosin displacing a conserved water molecule. A cluster of histidine and proline residues (His98-Pro99-His100-Pro101-His102) in kappa-casein binds to the C-terminal domain of the protein, where a neighboring conserved arginine residue (Arg97) is found to be important for stabilizing the binding pose. The catalytic site (including the catalytic water molecule) is stable in the starting conformation of the previously proposed general acid/base catalytic mechanism for 18 ns of molecular dynamics simulations.

Statistical coupling analysis of aspartic proteinases based on crystal structures of the Trichoderma reesei enzyme and its complex with pepstatin A

The crystal structures of an aspartic proteinase from Trichoderma reesei (TrAsP) and of its complex with a competitive inhibitor, pepstatin A, were solved and refined to crystallographic R-factors of 17.9% (R(free)=21.2%) at 1.70 A resolution and 15.8% (R(free)=19.2%) at 1.85 A resolution, respectively. The three-dimensional structure of TrAsP is similar to structures of other members of the pepsin-like family of aspartic proteinases. Each molecule is folded in a predominantly beta-sheet bilobal structure with the N-terminal and C-terminal domains of about the same size. Structural comparison of the native structure and the TrAsP-pepstatin complex reveals that the enzyme undergoes an induced-fit, rigid-body movement upon inhibitor binding, with the N-terminal and C-terminal lobes tightly enclosing the inhibitor. Upon recognition and binding of pepstatin A, amino acid residues of the enzyme active site form a number of short hydrogen bonds to the inhibitor that may play an important role in the mechanism of catalysis and inhibition. The structures of TrAsP were used as a template for performing statistical coupling analysis of the aspartic protease family. This approach permitted, for the first time, the identification of a network of structurally linked residues putatively mediating conformational changes relevant to the function of this family of enzymes. Statistical coupling analysis reveals coevolved continuous clusters of amino acid residues that extend from the active site into the hydrophobic cores of each of the two domains and include amino acid residues from the flap regions, highlighting the importance of these parts of the protein for its enzymatic activity.

The catalytic mechanism of an aspartic proteinase explored with neutron and X-ray diffraction

Hydrogen atoms play key roles in enzyme mechanism, but as this study shows, even high-quality X-ray data to a resolution of 1 A cannot directly visualize them. Neutron diffraction, however, can locate deuterium atoms even at resolutions around 2 A. Both neutron and X-ray diffraction data have been used to investigate the transition state of the aspartic proteinase endothiapepsin. The different techniques reveal a different part of the story, revealing the clearest picture yet of the catalytic mechanism by which the enzyme operates. Room temperature neutron and X-ray diffraction data were used in a newly developed joint refinement software package to visualize deuterium atoms within the active site of the enzyme when a gem-diol transition state analogue inhibitor is bound at the active site. These data were also used to estimate their individual occupancy, while analysis of the differences between the bond lengths of the catalytic aspartates was performed using atomic resolution X-ray data. The two methods are in agreement on the protonation state of the active site with a transition state analogue inhibitor bound confirming the catalytic mechanism at which the enzyme operates.

Human cathepsin D

A literature survey was performed of human cathepsin D gene, cathepsin D biosynthesis, posttranslatory modifications, transport within the cell, substrate specificity and catalytic effect. Methods used to determine the activity and level of this proteinase as well as its role in the biochemistry and pathobiochemistry of cells, tissues and organs were considered.

Serine protease mechanism-based mimics. Direct evidence for a transition state bridge proton in stable potentials

We have synthesized 1-(2-hydroxyacetyl)piperidine-2-one (2) and 1-(2-hydroxyacetyl)azepan-2-one (3). Equilibrium (K(f)) between the free alcohol (open form) and the tetrahedral intermediate (cyclol) is readily established, and both forms are observed in the D(2)O (1)H NMR spectra of 2 and 3. Therefore, their interconversion can be considered as an almost thermoneutral non-identical one. Pseudo-first-order rate constants (k(obs)) were obtained by simulating the AB (1)H NMR system observed for the cyclol. By best fitting the experimental points of a k(obs) versus pD profile to the equation k(obs) = 0.5k(0r) + 0.5k(r) K(ac)/(K(ac) + [D(+)]) + 0.5k(f)K(ao)/(K(ao)+ [D(+)]), the parameters involved were obtained: rate constants of rupture and formation (k(0r) and k(0f) = K(f)k(0r)) catalyzed by water, rate constants of rupture (k(r)) and formation (k(f)) from the conjugated bases of the cyclol form and the open form, and their acidity equilibrium constants K(ac) and K(ao). The system studied mimics the serine alcohol attack on the peptide bond and its reverse reaction in serine protease enzymes. In fact, the reaction rates are similar or perhaps even faster than the ones obtained for enzymatic reactions. The results also show the participation of water molecules forming catalytic proton bridges in stable potentials with the two interconverted forms. The position change of the bridged proton is sensitive to lactam ring size, and it is manifested by considerable change in the pKa values of both cyclol and open forms. Other evidence such as kinetics, DeltaS degrees , DeltaS, and proton inventory experiments and semiempirical molecular calculations support this proposal.

Preliminary neutron and ultrahigh-resolution X-ray diffraction studies of the aspartic proteinase endothiapepsin cocrystallized with a gem-diol inhibitor

Endothiapepsin has been cocrystallized with the gem-diol inhibitor PD-135,040 in a low solvent-content (39%) unit cell, which is unprecedented for this enzyme-inhibitor complex and enables ultrahigh-resolution (1.0 A) X-ray diffraction data to be collected. This atomic resolution X-ray data set will be used to deduce the protonation states of the catalytic aspartate residues. A room-temperature neutron data set has also been collected for joint refinement with a room-temperature X-ray data set in order to locate the H/D atoms at the active site.

Assigning the protonation states of the key aspartates in β-Secretase using QM/MM X-ray structure refinement

β-Secretase, a.k.a. β-APP cleaving enzyme (BACE), is an aspartyl protease that has been implicated as a key target in the pathogenesis of Alzheimer's disease (AD). The identification of the protonation states of the key aspartates in β-secretase is of great interest both in understanding the reaction mechanism and in guiding the design of drugs against AD. However, the resolutions of currently available crystal structures for BACE are not sufficient to determine the hydrogen atom locations. We have assigned the protonation states of the key aspartates using a novel method, QM/MM X-ray refinement. In our approach, an energy function is introduced to the refinement where the atoms in the active site are modeled by quantum mechanics (QM) and the other atoms are represented by molecular mechanics (MM). The gradients derived from the QM/MM energy function are combined with those from the X-ray target to refine the crystal structure of a complex containing BACE and an inhibitor. A total number of 8 protonation configurations of the aspartyl dyad were considered and QM/MM X-ray refinements were performed for all of them. The relative stability of the refined structures was scored by constructing the thermodynamic cycle using the energetics calculated by fully quantum mechanical self-consistent reaction field (QM/SCRF) calculations. While all 8 refined structures fit the observed electron density about equally well, we find the mono-protonated configurations to be strongly favored energetically, especially the configuration with the inner oxygen of Asp32 protonated and the hydroxyl of the inhibitor pointing towards Asp228. It was also found that these results depend on the constraints imposed by the X-ray data. We suggest that one of the strengths of this approach is that the resulting structures are a consensus of theoretical and experimental data and remark on the significance of our results in structure based drug design and mechanistic studies.

An unusual functional group interaction and its potential to reproduce steric and electrostatic features of the transition states of peptidolysis

The donor-acceptor interaction between a tertiary amine and an aldehyde, first observed among a select class of alkaloids, was deliberately established in a peptidomimetic (1a-c) to mimic features of the two principal transition states of peptide hydrolysis. Compounds 1a-c show preferential adoption in methanol and water of a 'folded' conformation displaying the interaction. Proportions of the folded form in MeOH range from 45% to 70% and can reach 84% in buffer. Significantly, three tendencies for the folded/unfolded equilibrium are observed: increasing solubility and polarity of the medium and decreasing temperature results in a higher extent of folding. In the absence of any parameter set available for this weak bond, no modeling studies were conducted to aid in the design of 1a-c. The successful straightforward synthesis of 1 and its folding and inhibition results with HIV-1 peptidase using FRET technology encourage studies to further pre-organize candidate molecules and to screen the structure space by modeling and parallel combinatorial chemistry.

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.

Molecular complexes at a glance: automated generation of two-dimensional complex diagrams

MOTIVATION

In this paper a new algorithmic approach is presented, which automatically generates structure diagrams of molecular complexes. A complex diagram contains the ligand, the amino acids of the protein interacting with the ligand and the hydrophilic interactions schematized as dashed lines between the corresponding atoms. The algorithm is based on a combinatorial optimization strategy which solves parts of the layout problem non-heuristically. The depicted molecules are represented as structure diagrams according to the chemical nomenclature. Due to the frequent usage of complex diagrams in the scientific literature as well as in text books dealing with structural biology, biochemistry and medicinal chemistry, the new algorithm is a key element for computer applications in these areas.

RESULTS

The method was implemented in the new software tool PoseView. It was tested on a representative dataset containing 305 protein-ligand complexes in total from the Brookhaven Protein Data Bank. PoseView was able to find collision-free layouts for more than three quarters of all complexes. In the following the layout generation algorithm is presented and, additional to the statistical results, representative test cases demonstrating the challenges of the layout generation will be discussed.

AVAILABILITY

The method is available as a webservice at http://www.zbh.uni-hamburg.de/poseview.

Evolutionarily conserved functional mechanics across pepsin-like and retroviral aspartic proteases

The biological function of the aspartic protease from HIV-1 has recently been related to the conformational flexibility of its structural scaffold. Here, we use a multistep strategy to investigate whether the same mechanism affects the functionality in the pepsin-like fold. (i) We identify the set of conserved residues by using sequence-alignment techniques. These residues cluster in three distinct regions: near the cleavage-site cavity, in the four beta-sheets cross-linking the two lobes, and in a solvent-exposed region below the long beta-hairpin in the N-terminal lobe. (ii) We elucidate the role played by the conserved residues for the enzymatic functionality of one representative member of the fold family, the human beta-secretase, by means of classical molecular dynamics (MD). The conserved regions exhibit little overall mobility and yet are involved into the most important modes of structural fluctuations. These modes influence the substrate-catalytic aspartates distance through a relative rotation of the N- and C-terminal lobes. (iii) We investigate the effects of this modulation by estimating the reaction free energy at different representative substrate/enzyme conformations. The activation free energy is strongly affected by large-scale protein motions, similarly to what has been observed in the HIV-1 enzyme. (iv) We extend our findings to all other members of the two eukaryotic and retroviral fold families by recurring to a simple, topology-based, energy functional. This analysis reveals a sophisticated mechanism of enzymatic activity modulation common to all aspartic proteases. We suggest that aspartic proteases have been evolutionarily selected to possess similar functional motions despite the observed fold variations.

Switching of turn conformation in an aspartate anion peptide fragment by NH · · · O−hydrogen bonds

Aspartic acid protease model peptides Z-Phe-Asp(COOH)-Thr-Gly-Ser-Ala-NHCy (1) and AdCO-Asp(COOH)-Val-Gly-NHBzl (3), and their aspartate anions (NEt4)[Z-Phe-Asp(COO-)-Thr-Gly-Ser-Ala-NHCy] (2) and (NEt4)[AdCO-Asp(COO-)-Val-Gly-NHBzl] (4), having an invariant primary sequence of the Asp-X(Thr,Ser)-Gly fragment, were synthesized and characterized by 1H-NMR, CD, and infrared (IR) spectroscopies. NMR structure analyses indicate that the Asp O(delta) atoms of the aspartate peptide 2 are intramolecularly hydrogen-bonded with Gly, Ser, Ala NH, and Ser OH, supporting the rigid beta-turn-like conformation in acetonitrile solution. The tripeptide in the aspartic acid 3 forms an inverse gamma-turn structure, which is converted to a beta-turn-like conformation because of the formation of the intramolecular NH . . . O- hydrogen bonds with the Asp O(delta) in 4. Such a conformational change is not detected between dipeptides AdCO-Asp(COOH)-Va-NHAd (5) and (NEt4)[AdCO-Asp(COO-)-Val-NHAd] (6). The pK(a) value of side-chain carboxylic acid (5.0) for 3 exhibits a lower shift (0.3 unit) from that of 5 in aqueous polyethyleneglycol lauryl ether micellar solution. NMR structure analyses for 3 in an aqueous micellar solution indicate that the preorganized turn structure, which readily forms the NH . . . O- hydrogen bonds, lowers the pK(a) value and that resulting hydrogen bonds stabilize the rigid conformation in the aspartate anion state. We found that the formation of the NH . . . O- hydrogen bonds involved in the hairpin turn is correlated with the protonation and deprotonation state of the Asp side chain in the conserved amino acid fragments.

The structure of the C-C bond hydrolase MhpC provides insights into its catalytic mechanism

2-Hydroxy-6-ketonona-2,4-diene-1,9-dioic acid 5,6-hydrolase (MhpC) is a 62 kDa homodimeric enzyme of the phenylpropionate degradation pathway of Escherichia coli. The 2.1 A resolution X-ray structure of the native enzyme determined from orthorhombic crystals confirms that it is a member of the alpha/beta hydrolase fold family, comprising eight beta-strands interconnected by loops and helices. The 2.8 A resolution structure of the enzyme co-crystallised with the non-hydrolysable substrate analogue 2,6-diketo-nona-1,9-dioic acid (DKNDA) confirms the location of the active site in a buried channel including Ser110, His263 and Asp235, postulated contributors to a serine protease-like catalytic triad in homologous enzymes. It appears that the ligand binds in two separate orientations. In the first, the C6 keto group of the inhibitor forms a hemi-ketal adduct with the Ser110 side-chain, the C9 carboxylate group interacts, via the intermediacy of a water molecule, with Arg188 at one end of the active site, while the C1 carboxylate group of the inhibitor comes close to His114 at the other end. In the second orientation, the C1 carboxylate group binds at the Arg188 end of the active site and the C9 carboxylate group at the His114 end. These arrangements implicated His114 or His263 as plausible contributors to catalysis of the initial enol/keto tautomerisation of the substrate but lack of conservation of His114 amongst related enzymes and mutagenesis results suggest that His263 is the residue involved. Variability in the quality of the electron density for the inhibitor amongst the eight molecules of the crystal asymmetric unit appears to correlate with alternative positions for the side-chain of His114. This might arise from half-site occupation of the dimeric enzyme and reflect the apparent dissociation of approximately 50% of the keto intermediate from the enzyme during the catalytic cycle.

Modeling the protonation states of the catalytic aspartates in beta-secretase

Beta-secretase (BACE) is a critical enzyme in the production of beta-amyloid, a protein that has been implicated as a potential cause of Alzheimer's disease (AD). There are two aspartic acid residues (Asp 32 and Asp 228) present in the catalytic region of BACE that can adopt multiple protonation states. The protonation state and precise location of the protons for these two residues, particularly in the presence of an inhibitor, are subjects of great interest since they have a direct bearing on the mechanism of aspartyl proteases and efforts to model beta-secretase. We have carried out full liner-scaling quantum mechanical (QM) calculations that include Poisson-Boltzmann solvation in order to identify the preferred protonation state and proton location in the presence and absence of an inhibitor. These calculations favor the monoprotonated state in the presence of ligand, and di-deprotonated state in the absence of ligand. Further the proton in the monoprotonated state is located on the inner oxygen of Asp 228. These results have implications for the catalytic mechanism of BACE and related aspartyl proteases. They also provide a reference state for the protein in structure-based modeling studies of this therapeutically important target.

Crystal structure of aspartic proteinase from Irpex lacteus in complex with inhibitor pepstatin

The crystal structure of Irpex lacteus aspartic proteinase (ILAP) in complex with pepstatin (a six amino acid residue peptide-like inhibitor) was determined at 1.3A resolution. ILAP is a pepsin-like enzyme, widely distributed in nature, with high milk-clotting activity relative to proteolytic activity. The overall structure was in good topological agreement with pepsin and other aspartic proteases. The structure and interaction pattern around the catalytic site were conserved, in agreement with the other aspartic proteinase/inhibitor complex structures reported previously. The high-resolution data also supported the transition state model, as proposed previously for the catalytic mechanism of aspartic proteinase. Unlike the other aspartic proteinases, ILAP was found to require hydrophobic residues either in the P(1) or P(1') site, and also in the P(4) and/or P(3) site(s) for secondary interactions. The inhibitor complex structure also revealed the substrate binding mechanism of ILAP at the P(3) and P(4) site of the substrate, where the inserted loop built up the unique hydrophobic pocket at the P(4) site.

Atomic resolution analysis of the catalytic site of an aspartic proteinase and an unexpected mode of binding by short peptides

The X-ray structures of native endothiapepsin and a complex with a hydroxyethylene transition state analog inhibitor (H261) have been determined at atomic resolution. Unrestrained refinement of the carboxyl groups of the enzyme by using the atomic resolution data indicates that both catalytic aspartates in the native enzyme share a single negative charge equally; that is, in the crystal, one half of the active sites have Asp 32 ionized and the other half have Asp 215 ionized. The electron density map of the native enzyme refined at 0.9 A resolution demonstrates that there is a short peptide (probably Ser-Thr) bound noncovalently in the active site cleft. The N-terminal nitrogen of the dipeptide interacts with the aspartate diad of the enzyme by hydrogen bonds involving the carboxyl of Asp 215 and the catalytic water molecule. This is consistent with classical findings that the aspartic proteinases can be inhibited weakly by short peptides and that these enzymes can catalyze transpeptidation reactions. The dipeptide may originate from autolysis of the N-terminal Ser-Thr sequence of the enzyme during crystallization.

Crystal structures of Aspergillus oryzae aspartic proteinase and its complex with an inhibitor pepstatin at 1.9A resolution

The X-ray structures of Aspergillus oryzae aspartic proteinase (AOAP) and its complex with inhibitor pepstatin have been determined at 1.9A resolution. AOAP was crystallized in an orthorhombic system with the space group P2(1)2(1)2(1) and cell dimensions of a=49.4A, b=79.4A, and c=93.6A. By the soaking of pepstatin, crystals are transformed into a monoclinic system with the space group C2 and cell dimensions of a=106.8A, b=38.6A, c=78.7A, and beta=120.3 degrees. The structures of AOAP and AOAP/pepstatin complex were refined to an R-factor of 0.177 (R(free)=0.213) and of 0.185 (0.221), respectively. AOAP has a crescent-shaped structure with two lobes (N-lobe and C-lobe) and the deep active site cleft is constructed between them. At the center of the active site cleft, two Asp residues (Asp33 and Asp214) form the active dyad with a hydrogen bonding solvent molecule between them. Pepstatin binds to the active site cleft via hydrogen bonds and hydrophobic interactions with the enzyme. The structures of AOAP and AOAP/pepstatin complex including interactions between the enzyme and pepstatin are very similar to those of other structure-solved aspartic proteinases and their complexes with pepstatin. Generally, aspartic proteinases cleave a peptide bond between hydrophobic amino acid residues, but AOAP can also recognize the Lys/Arg residue as well as hydrophobic amino acid residues, leading to the activation of trypsinogen and chymotrypsinogen. The X-ray structure of AOAP/pepstatin complex and preliminary modeling show two possible sites of recognition for the positively charged groups of Lys/Arg residues around the active site of AOAP.

A history of pepsin and related enzymes

Studies on gastric digestion during 1820-1840 led to the discovery of pepsin as the agent which, in the presence of stomach acid, causes the dissolution of nutrients such as meat or coagulated egg white. Soon afterward it was shown that these protein nutrients were cleaved by pepsin to diffusible products named peptones. Efforts to isolate and purify pepsin were spurred by its widespread adoption for the treatment of digestive disorders, and highly active preparations were available by the end of the nineteenth century. There was uncertainty, however, as to the chemical nature of pepsin, for some preparations exhibited the properties of proteins while other preparations failed to do so. The question was not settled until after 1930, when Northrop crystallized swine pepsin and provided convincing evidence for its identity as a protein. The availability of this purified pepsin during the 1930s also led to the discovery of the first synthetic peptide substrates for pepsin, thus providing needed evidence for the peptide structure of native proteins, a matter of debate at that time. After 1945, with the introduction of new separation methods, notably chromatography and electrophoresis, and the availability of specific proteinases, the amino acid sequences of many proteins, including pepsin and its precursor pepsinogen, were determined. Moreover, treatment of pepsin with chemical reagents indicated the participation in the catalytic mechanism of two aspartyl units widely separated in the linear sequence. Studies on the kinetics of pepsin action on long chain synthetic peptides suggested that the catalytic site was an extended structure. Similar properties were found for other "aspartyl proteinases," such as chymosin (used in cheese making), some intracellular proteinases (cathepsins), and plant proteinases. After 1975, the three-dimensional structures of pepsin and many of its relatives were determined by means of x-ray diffraction techniques, greatly extending our insight into the mechanism of the catalytic action of these enzymes. That knowledge has led to the design of new inhibitors of aspartyl proteinases, which are participants in the maturation of human immunodeficiency virus and in the generation of Alzheimer's disease.