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

Extracting binding energies and binding modes from biomolecular simulations of fragment binding to endothiapepsin

Fragment-based drug discovery (FBDD) aims to discover a set of small binding fragments that may be subsequently linked together. Therefore, in-depth knowledge of the individual fragments' structural and energetic binding properties is essential. In addition to experimental techniques, the direct simulation of fragment binding by molecular dynamics (MD) simulations became popular to characterize fragment binding. However, former studies showed that long simulation times and high computational demands per fragment are needed, which limits applicability in FBDD. Here, we performed short, unbiased MD simulations of direct fragment binding to endothiapepsin, a well-characterized model system of pepsin-like aspartic proteases. To evaluate the strengths and limitations of short MD simulations for the structural and energetic characterization of fragment binding, we predicted the fragments' absolute free energies and binding poses based on the direct simulations of fragment binding and compared the predictions to experimental data. The predicted absolute free energies are in fair agreement with the experiment. Combining the MD data with binding mode predictions from molecular docking approaches helped to correctly identify the most promising fragments for further chemical optimization. Importantly, all computations and predictions were done within 5 days, suggesting that MD simulations may become a viable tool in FBDD projects.

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.

An in silico approach for identification of novel inhibitors as potential therapeutics targeting COVID-19 main protease

Respiratory disease caused by a novel coronavirus, COVID-19, has been labeled a pandemic by the World Health Organization. Very little is known about the infection mechanism for this virus. More importantly, there are no drugs or vaccines that can cure or prevent a person from getting COVID-19. In this study, the binding affinity of 2692 protease inhibitor compounds that are known in the protein data bank, are calculated against the main protease of the novel coronavirus with docking and molecular dynamics (MD). Both the docking and MD methods predict the macrocyclic tissue factor-factor VIIa (PubChem ID: 118098670) inhibitor to bind strongly with the main protease with a binding affinity of -10.6 and -10.0 kcal/mol, respectively. The TF-FVIIa inhibitors are known to prevent the coagulation of blood and have antiviral activity as shown in the case of SARS coronavirus. Two more inhibitors, phenyltriazolinones (PubChem ID: 104161460) and allosteric HCV NS5B polymerase thumb pocket 2 (PubChem ID: 163632044) have shown antiviral activity and also have high affinity towards the main protease of COVID-19. Furthermore, these inhibitors interact with the catalytic dyad in the active site of the COVID-19 main protease that is especially important in viral replication. The calculated theoretical dissociation constants of the proposed COVID-19 inhibitors are found to be very similar to the experimental dissociation constant values of similar protease-inhibitor systems.Communicated by Ramaswamy H. Sarma.

Advancements in macromolecular crystallography: from past to present

Protein Crystallography or Macromolecular Crystallography (MX) started as a new discipline of science with the pioneering work on the determination of the protein crystal structures by John Kendrew in 1958 and Max Perutz in 1960. The incredible achievements in MX are attributed to the development of advanced tools, methodologies, and automation in every aspect of the structure determination process, which have reduced the time required for solving protein structures from years to a few days, as evident from the tens of thousands of crystal structures of macromolecules available in PDB. The advent of brilliant synchrotron sources, fast detectors, and novel sample delivery methods has shifted the paradigm from static structures to understanding the dynamic picture of macromolecules; further propelled by X-ray Free Electron Lasers (XFELs) that explore the femtosecond regime. The revival of the Laue diffraction has also enabled the understanding of macromolecules through time-resolved crystallography. In this review, we present some of the astonishing method-related and technological advancements that have contributed to the progress of MX. Even with the rapid evolution of several methods for structure determination, the developments in MX will keep this technique relevant and it will continue to play a pivotal role in gaining unprecedented atomic-level details as well as revealing the dynamics of biological macromolecules. With many exciting developments awaiting in the upcoming years, MX has the potential to contribute significantly to the growth of modern biology by unraveling the mechanisms of complex biological processes as well as impacting the area of drug designing.

Low barrier hydrogen bonds in protein structure and function

Low-barrier hydrogen bonds (LBHBs) are a special type of short hydrogen bond (HB) that is characterized by the equal sharing of a hydrogen atom. The existence and catalytic role of LBHBs in proteins has been intensely contested. Advancements in X-ray and neutron diffraction methods has revealed delocalized hydrogen atoms involved in potential LBHBs in a number of proteins, while also demonstrating that short HBs are not necessarily LBHBs. More importantly, a series of experiments on ketosteroid isomerase (KSI) have suggested that LBHBs are significantly stronger than standard HBs in the protein microenvironment in terms of enthalpy, but not free energy. The discrepancy between the enthalpy and free energy of LBHBs offers clues to the challenges, and potential solutions, of the LBHB debate, where the unique strength of LBHBs plays a special role in the kinetic processes of enzyme function and structure, together with other molecular forces in a pre-organized environment.

Characterization of the endothiapepsin-like protein in the entomopathogenic fungus Beauveria bassiana and its virulence effect on the silkworm, Bombyx mori

Endothiapepsin is an aspartic proteinase that was first isolated from the plant pathogenic fungus Endothia parasitica. In previous studies, we reported on three endothiapepsin-like proteins in the entomopathogenic fungus Beauveria bassiana; the genes were up-regulated in B. bassiana hyper-virulent strain GXsk1011 at early stage infection in the silkworm. However, whether these proteins play a role in pathogenicity or not remains unknown. In this study, we cloned one protein, BbepnL-1 gene (BBA-07766), that has 98% homology with B. bassiana strain Bb2860, and expressed it in the yeast Pichia pastoris to investigate its function. The endothiapepsin-like protein is a secreted proteinase of molecular weight approximately 40 kDa. It has an N-glycosylation site and a mutation in the C-terminal conserved domain- a Thr was mutated to Gly in B. bassiana GXsk1011 and is different than the endothiapepsin of Endothia parasitica. The recombinant endothiapepsin-like protein showed enzyme activity and degraded the protein components of the silkworm cuticle. To further investigate the activity of the endothiapepsin-like protein, we knocked out the gene BbepnL-1 and showed that the loss of BbepnL-1 reduced the virulence in the silkworm. These results demonstrated that the endothiapepsin-like protein of B. bassiana is a virulence factor.

Deciphering the mechanism of potent peptidomimetic inhibitors targeting plasmepsins - biochemical and structural insights

Malaria is a deadly disease killing worldwide hundreds of thousands people each year and the responsible parasite has acquired resistance to the available drug combinations. The four vacuolar plasmepsins (PMs) in Plasmodium falciparum involved in hemoglobin (Hb) catabolism represent promising targets to combat drug resistance. High antimalarial activities can be achieved by developing a single drug that would simultaneously target all the vacuolar PMs. We have demonstrated for the first time the use of soluble recombinant plasmepsin II (PMII) for structure-guided drug discovery with KNI inhibitors. Compounds used in this study (KNI-10742, 10743, 10395, 10333, and 10343) exhibit nanomolar inhibition against PMII and are also effective in blocking the activities of PMI and PMIV with the low nanomolar K_i values. The high-resolution crystal structures of PMII-KNI inhibitor complexes reveal interesting features modulating their differential potency. Important individual characteristics of the inhibitors and their importance for potency have been established. The alkylamino analog, KNI-10743, shows intrinsic flexibility at the P2 position that potentiates its interactions with Asp132, Leu133, and Ser134. The phenylacetyl tripeptides, KNI-10333 and KNI-10343, accommodate different ρ-substituents at the P3 phenylacetyl ring that determine the orientation of the ring, thus creating novel hydrogen-bonding contacts. KNI-10743 and KNI-10333 possess significant antimalarial activity, block Hb degradation inside the food vacuole, and show no cytotoxicity on human cells; thus, they can be considered as promising candidates for further optimization. Based on our structural data, novel KNI derivatives with improved antimalarial activity could be designed for potential clinical use. DATABASE: Structural data are available in the PDB under the accession numbers 5YIE, 5YIB, 5YID, 5YIC, and 5YIA.

Short Carboxylic Acid-Carboxylate Hydrogen Bonds Can Have Fully Localized Protons

Short hydrogen bonds (H-bonds) have been proposed to play key functional roles in several proteins. The location of the proton in short H-bonds is of central importance, as proton delocalization is a defining feature of low-barrier hydrogen bonds (LBHBs). Experimentally determining proton location in H-bonds is challenging. Here, bond length analysis of atomic (1.15-0.98 Å) resolution X-ray crystal structures of the human protein DJ-1 and its bacterial homologue, YajL, was used to determine the protonation states of H-bonded carboxylic acids. DJ-1 contains a buried, dimer-spanning 2.49 Å H-bond between Glu15 and Asp24 that satisfies standard donor-acceptor distance criteria for a LBHB. Bond length analysis indicates that the proton is localized on Asp24, excluding a LBHB at this location. However, similar analysis of the Escherichia coli homologue YajL shows both residues may be protonated at the H-bonded oxygen atoms, potentially consistent with a LBHB. A Protein Data Bank-wide screen identifies candidate carboxylic acid H-bonds in approximately 14% of proteins, which are typically short [⟨d_O-O⟩ = 2.542(2) Å]. Chemically similar H-bonds between hydroxylated residues (Ser/Thr/Tyr) and carboxylates show a trend of lengthening O-O distance with increasing H-bond donor pK_a. This trend suggests that conventional electronic effects provide an adequate explanation for short, charge-assisted carboxylic acid-carboxylate H-bonds in proteins, without the need to invoke LBHBs in general. This study demonstrates that bond length analysis of atomic resolution X-ray crystal structures provides a useful experimental test of certain candidate LBHBs.

Atomic resolution crystal structure of Sapp2p, a secreted aspartic protease from Candida parapsilosis

The virulence of the Candida pathogens is enhanced by the production of secreted aspartic proteases, which therefore represent possible targets for drug design. Here, the crystal structure of the secreted aspartic protease Sapp2p from Candida parapsilosis was determined. Sapp2p was isolated from its natural source and crystallized in complex with pepstatin A, a classical aspartic protease inhibitor. The atomic resolution of 0.83 Å allowed the protonation states of the active-site residues to be inferred. A detailed comparison of the structure of Sapp2p with the structure of Sapp1p, the most abundant C. parapsilosis secreted aspartic protease, was performed. The analysis, which included advanced quantum-chemical interaction-energy calculations, uncovered molecular details that allowed the experimentally observed equipotent inhibition of both isoenzymes by pepstatin A to be rationalized.

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.

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 small nonrule of 3 compatible fragment library provides high hit rate of endothiapepsin crystal structures with various fragment chemotypes

Druglike molecules are defined by Lipinski's rule of 5, to characterize fragment thresholds, they have been reduced from 5 to 3 (Astex's rule of 3). They are applied to assemble fragment libraries, and providers use them to select fragments for commercial offer. We question whether these rules are too stringent to compose fragment libraries with candidates exhibiting sufficient room for chemical subsequent growing and merging modifications as appropriate functional groups for chemical transformations are required. Usually these groups exhibit properties as hydrogen bond donors/acceptors and provide entry points for optimization chemistry. We therefore designed a fragment library (364 entries) without strictly applying the rule of 3. For initial screening for endothiapepsin binding, we performed a biochemical cleavage assay of a fluorogenic substrate at 1 mM. "Hits" were defined to inhibit the enzyme by at least 40%. Fifty-five hits were suggested and subsequently soaked into endothiapepsin crystals. Eleven crystal structures could be determined covering fragments with diverse binding modes: (i) direct binding to the catalytic dyad aspartates, (ii) water-mediated binding to the aspartates, (iii) no direct interaction with the dyad. They occupy different specificity pockets. Only 4 of the 11 fragments are consistent with the rule of 3. Restriction to this rule would have limited the fragment hits to a strongly reduced variety of chemotypes.

Internal motion in protein crystal structures

The binding states of the substrates and the environment have significant influence on protein motion. We present the analysis of such motion derived from anisotropic atomic displacement parameters (ADPs) in a set of atomic resolution protein structures. Local structural motion caused by ligand binding as well as functional loops showing cooperative patterns of motion could be inferred. The results are in line with proposed protonation states, hydrogen bonding patterns and the location of distinctly flexible regions: we could locate the mobile active site loop in a virus integrase, distinguish the subdomains in RNAse A and hydroxynitrile lyase, and reconstruct the molecular architecture in a xylanase. We demonstrate that the ADP-based motion analysis provides information at high level of detail and that the structural changes needed for substrate attachment or release may be derived from single X-ray structures.

Synthesis of α-carboxyphosphinopeptides derived from norleucine

In the present study, we describe in detail the synthesis of a relatively rare class of phosphorus compounds, α-carboxyphosphinopeptides. We prepared several norleucine-derived α-carboxyphosphinic pseudopeptides of the general formula Nle-Ψ[PO(OH)]-Gly. These compounds could have important applications as transition state-mimicking inhibitors for methionine or leucine aminopeptidases or other enzymes. For the preparation of the key α-carboxyphosphinate protected precursors, we investigated, compared and improved two different synthetic methods described in literature: the Arbuzov reaction of a silylated N-protected phosphinic acid with a bromoacetate ester and the nucleophilic addition of a mixed O-methyl S-phenyl N-protected phosphonic acid or a methyl N-protected phosphonochloridate with tert-butyl lithioacetate. We also prepared two N-Fmoc protected synthons, Fmoc-Nle-Ψ[PO(OH)]-Gly-COOH and Fmoc-Nle-Ψ[PO(OAd)]-Gly-COOH, and demonstrated that these precursors are suitable building blocks for the solid-phase synthesis of α-carboxyphosphinopeptides.

Cysteine pKa depression by a protonated glutamic acid in human DJ-1

Human DJ-1, a disease-associated protein that protects cells from oxidative stress, contains an oxidation-sensitive cysteine (C106) that is essential for its cytoprotective activity. The origin of C106 reactivity is obscure, due in part to the absence of an experimentally determined p K a value for this residue. We have used atomic-resolution X-ray crystallography and UV spectroscopy to show that C106 has a depressed p K a of 5.4 +/- 0.1 and that the C106 thiolate accepts a hydrogen bond from a protonated glutamic acid side chain (E18). X-ray crystal structures and cysteine p K a analysis of several site-directed substitutions at residue 18 demonstrate that the protonated carboxylic acid side chain of E18 is required for the maximal stabilization of the C106 thiolate. A nearby arginine residue (R48) participates in a guanidinium stacking interaction with R28 from the other monomer in the DJ-1 dimer and elevates the p K a of C106 by binding an anion that electrostatically suppresses thiol ionization. Our results show that the ionizable residues (E18, R48, and R28) surrounding C106 affect its p K a in a way that is contrary to expectations based on the typical ionization behavior of glutamic acid and arginine. Lastly, a search of the Protein Data Bank (PDB) produces several candidate hydrogen-bonded aspartic/glutamic acid-cysteine interactions, which we propose are particularly common in the DJ-1 superfamily.

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.

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.

NMR Study of the Inhibition of Pepsin by Glyoxal Inhibitors: Mechanism of Tetrahedral Intermediate Stabilization by the Aspartyl Proteases

Z-Ala-Ala-Phe-glyoxal (where Z is benzyloxycarbonyl) has been shown to be a competitive inhibitor of pepsin with a Ki = 89 +/- 24 nM at pH 2.0 and 25 degrees C. Both the ketone carbon (R13COCHO) and the aldehyde carbon (RCO13CHO) of the glyoxal group of Z-Ala-Ala-Phe-glyoxal have been 13C-enriched. Using 13C NMR, it has been shown that when the inhibitor is bound to pepsin, the glyoxal keto and aldehyde carbons give signals at 98.8 and 90.9 ppm, respectively. This demonstrates that pepsin binds and preferentially stabilizes the fully hydrated form of the glyoxal inhibitor Z-Ala-Ala-Phe-glyoxal. From 13C NMR pH studies with glyoxal inhibitor, we obtain no evidence for its hemiketal or hemiacetal hydroxyl groups ionizing to give oxyanions. We conclude that if an oxyanion is formed its pKa must be >8.0. Using 1H NMR, we observe four hydrogen bonds in free pepsin and in pepsin/Z-Ala-Ala-Phe-glyoxal complexes. In the pepsin/pepstatin complex an additional hydrogen bond is formed. We examine the effect of pH on hydrogen bond formation, but we do not find any evidence for low-barrier hydrogen bond formation in the inhibitor complexes. We conclude that the primary role of hydrogen bonding to catalytic tetrahedral intermediates in the aspartyl proteases is to correctly orientate the tetrahedral intermediate for catalysis.

Catalysis and Linear Free Energy Relationships in Aspartic Proteases

Aspartic proteases are receiving considerable attention as potential drug targets in several serious diseases, such as AIDS, malaria, and Alzheimer's disease. These enzymes cleave polypeptide chains, often between specific amino acid residues, but despite the common reaction mechanism, they exhibit large structural differences. Here, the catalytic mechanism of aspartic proteases plasmepsin II, cathepsin D, and HIV-1 protease is examined by computer simulations utilizing the empirical valence bond approach in combination with molecular dynamics and free energy perturbation calculations. Free energy profiles are established for four different substrates, each six amino acids long and containing hydrophobic side chains in the P1 and P1' positions. Our simulations reproduce the catalytic effect of these enzymes, which accelerate the reaction rate by a factor of approximately 10(10) compared to that of the corresponding uncatalyzed reaction in water. The calculations elucidate the origin of the catalytic effect and allow a rationalization of the fact that, despite large structural differences between plasmepsin II/cathepsin D and HIV-1 protease, the magnitude of their rate enhancement is very similar. Amino acid residues surrounding the active site together with structurally conserved water molecules are found to play an important role in catalysis, mainly through dipolar (electrostatic) stabilization. A linear free energy relationship for the reactions in the different enzymes is established that also demonstrates the reduced reorganization energy in the enzymes compared to that in the uncatalyzed water reaction.

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.

Binding of Phosphinate and Phosphonate Inhibitors to Aspartic Proteases: A First-Principles Study

Phosphinate and phosphonate derivatives are potent inhibitors of aspartic proteases (APs). The affinity for the enzyme might be caused by the presence of low barrier hydrogen bonds between the ligand and the catalytic Asp dyad in the cleavage site. We have used density functional theory calculations along with hybrid quantum mechanics/molecular mechanics Car-Parrinello molecular dynamics simulations to investigate the hydrogen-bonding pattern at the binding site of the complexes of human immunodeficiency virus type-1 AP and the eukaryotic endothiapepsin and penicillopepsin. Our calculations are in fair agreement with the NMR data available for endothiapepsin (Coates et al. J. Mol. Biol. 2002, 318, 1405-1415) and show that the most stable active site configuration is the diprotonated, negatively charged form. In the viral complex both protons are located at the catalytic Asp dyad, while in the eukaryotic complexes the proton shared by the closest oxygen atoms is located at the phosphinic/phosphonic group.

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.

A phenylnorstatine inhibitor binding to HIV-1 protease: geometry, protonation, and subsite-pocket interactions analyzed at atomic resolution

The X-ray structure of a complex of HIV-1 protease (PR) with a phenylnorstatine inhibitor Z-Pns-Phe-Glu-Glu-NH(2) has been determined at 1.03 A, the highest resolution so far reported for any HIV PR complex. The inhibitor shows subnanomolar K(i) values for both the wild-type PR and the variant representing one of the most common mutations linked to resistance development. The structure comprising the phenylnorstatine moiety of (2R,3S)-chirality displays a unique pattern of hydrogen bonding to the two catalytic aspartate residues. This high resolution makes it possible to assess the donor and acceptor relations of this hydrogen bonding and to indicate a proton shared by the two catalytic residues. A structural mechanism for the unimpaired inhibition of the protease Val82Ala mutant is also suggested, based on energy calculations and analyses.

Characterization of Fluxional Hydrogen-Bonded Complexes of Acetic Acid and Acetate by NMR: Geometries and Isotope and Solvent Effects

1H, (2)H, and (13)C NMR spectra of enriched CH(3)(13)COOH acid without and in the presence of tetra-n-butylammonium acetate have been measured around 110 K using a liquefied Freon mixture CDF(3)/CDF(2)Cl as a solvent, as a function of the deuterium fraction in the mobile proton sites. For comparison, spectra were also taken of the adduct CH(3)(13)COOH.SbCl(5) 1 and of CH(2)Cl(13)COOH under similar conditions, as well as of CH(3)(13)COOH and CH(3)(13)COO(-) dissolved in H(2)O and D(2)O at low and high pH at 298 K. The low temperatures employed allowed us to detect several well-known and novel hydrogen-bonded complexes in the slow hydrogen bond exchange regime and to determine chemical shifts and coupling constants as well as H/D isotope effects on chemical shifts from the fine structure of the corresponding signals. The measurements show that self-association of both carboxylic acids in Freon solution gives rise exclusively to the formation of cyclic dimers 2 and 3 exhibiting a rapid degenerate double proton transfer. For the first time, a two-bond coupling of the type (2)J(CH(3)COOH) between a hydrogen-bonded proton and the carboxylic carbon has been observed, which is slightly smaller than half of the value observed for 1. In addition, the (1)H and (2)H chemical shifts of the HH, HD, and the DD isotopologues of 2 and 3 have been determined as well as the corresponding HH/HD/DD isotope effects on the (13)C chemical shifts. Similar "primary", "vicinal", and "secondary" isotope effects were observed for the novel 2:1 complex "dihydrogen triacetate" 5 between acetic acid and acetate. Another novel species is the 3:1 complex "trihydrogen tetraacetate" 6, which was also characterized by a complex degenerate combined hydrogen bond- and proton-transfer process. For comparison, the results obtained previously for hydrogen diacetate 4 and hydrogen maleate 7 are discussed. Using an improved (1)H chemical shift-hydrogen bond geometry correlation, the chemical shift data are converted into hydrogen bond geometries. They indicate cooperative hydrogen bonds in the cyclic dimers; i.e., widening of a given hydrogen bond by H/D substitution also widens the other coupled hydrogen bond. By contrast, the hydrogen bonds in 5 are anticooperative. The measurements show that ionization shifts the (13)C signal of the carboxyl group to low field when the group is immersed in water, but to high field when it is embedded in a polar aprotic environment. This finding allows us to understand the unusual ionization shift of aspartate groups in the HIV-pepstatin complex observed by Smith, R.; Brereton, I. M.; Chai, R. Y.; Kent, S. B. H. Nature Struct. Biol. 1996, 3, 946. It is demonstrated that the Freon solvents used in this study are better environments for model studies of amino acid interactions than aqueous or protic environments. Finally, a novel correlation of the hydrogen bond geometries with the H/D isotope effects on the (13)C chemical shifts of carboxylic acid groups is proposed, which allows one to estimate the hydrogen bond geometries and protonation states of these groups. It is shown that absence of such an isotope effect is not only compatible with an isolated carboxylate group but also with the presence of a short and strong hydrogen bond.

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.

Kinetic studies on beta-site amyloid precursor protein-cleaving enzyme (BACE). Confirmation of an iso mechanism

The steady-state kinetic mechanism of beta-amyloid precursor protein-cleaving enzyme (BACE)-catalyzed proteolytic cleavage was evaluated using product and statine- (Stat(V)) or hydroxyethylene-containing (OM99-2) peptide inhibition data, solvent kinetic isotope effects, and proton NMR spectroscopy. The noncompetitive inhibition pattern observed for both cleavage products, together with the independence of Stat(V) inhibition on substrate concentration, suggests a uni-bi-iso kinetic mechanism. According to this mechanism, the enzyme undergoes multiple conformation changes during the catalytic cycle. If any of these steps are rate-limiting to turnover, an enzyme form preceding the rate-limiting conformational change should accumulate. An insignificant solvent kinetic isotope effect (SKIE) on k(cat)/K(m), a large inverse solvent kinetic isotope effect on k(cat), and the absence of any SKIE on the inhibition onset by Stat(V) during catalysis together indicate that the rate-limiting iso-step occurs after formation of a tetrahedral intermediate. A moderately short and strong hydrogen bond (at delta 13.0 ppm and phi of 0.6) has been observed by NMR spectroscopy in the enzyme-hydroxyethylene peptide (OM99-2) complex that presumably mimics the tetrahedral intermediate of catalysis. Collapse of this intermediate, involving multiple steps and interconversion of enzyme forms, has been suggested to impose a rate limitation, which is manifested in a significant SKIE on k(cat). Multiple enzyme forms and their distribution during catalysis were evaluated by measuring the SKIE on the noncompetitive (mixed) inhibition constants for the C-terminal reaction product. Large, normal SKIE values were observed for these inhibition constants, suggesting that both kinetic and thermodynamic components contribute to the K(ii) and K(is) expressions, as has been suggested for other iso-mechanism featuring enzymes. We propose that a conformational change related to the reprotonation of aspartates during or after the bond-breaking event is the rate-limiting segment in the catalytic reaction of beta-amyloid precursor protein-cleaving enzyme, and ligands binding to other than the ground-state forms of the enzyme might provide inhibitors of greater pharmacological relevance.