Human alpha-thrombin inhibition by the highly selective compounds N-ethoxycarbonyl-D-Phe-Pro-alpha-azaLys p-nitrophenyl ester and N-carbobenzoxy-Pro-alpha-azaLys p-nitrophenyl ester: a kinetic, thermodynamic and X-ray crystallographic study

Functional Characterization of Primordial Protein Repair Enzyme M38 Metallo-Peptidase From Fervidobacterium islandicum AW-1

The NA23_RS08100 gene of Fervidobacterium islandicum AW-1 encodes a keratin-degrading β-aspartyl peptidase (FiBAP) that is highly expressed under starvation conditions. Herein, we expressed the gene in Escherichia coli, purified the recombinant enzyme to homogeneity, and investigated its function. The 318 kDa recombinant FiBAP enzyme exhibited maximal activity at 80°C and pH 7.0 in the presence of Zn²⁺. Size-exclusion chromatography revealed that the native enzyme is an octamer comprising a tetramer of dimers; this was further supported by determination of its crystal structure at 2.6 Å resolution. Consistently, the structure of FiBAP revealed three additional salt bridges in each dimer, involving 12 ionic interactions that might contribute to its high thermostability. In addition, the co-crystal structure containing the substrate analog N-carbobenzoxy-β-Asp-Leu at 2.7 Å resolution revealed binuclear Zn²⁺-mediated substrate binding, suggesting that FiBAP is a hyperthermophilic type-I IadA, in accordance with sequence-based phylogenetic analysis. Indeed, complementation of a Leu auxotrophic E. coli mutant strain (ΔiadA and ΔleuB) with FiBAP enabled the mutant strain to grow on isoAsp-Leu peptides. Remarkably, LC-MS/MS analysis of soluble keratin hydrolysates revealed that FiBAP not only cleaves the C-terminus of isoAsp residues but also has a relatively broad substrate specificity toward α-peptide bonds. Moreover, heat shock-induced protein aggregates retarded bacterial growth, but expression of BAP alleviated the growth defect by degrading damaged proteins. Taken together, these results suggest that the viability of hyperthermophiles under stressful conditions may rely on the activity of BAP within cellular protein repair systems.

Enzyme-Inhibitor Interactions and a Simple, Rapid Method for Determining Inhibition Modality

Contemporary chemical biology and drug discovery are increasingly focused on the discovery of inhibitory molecules that interact with enzyme targets in specific ways, such as allosteric or orthosteric binding. Hence, there is increasing interest in evaluating hit compounds from high-throughput diversity screening to determine their mode of interaction with the target. In this work, the common inhibition modalities are reviewed and clarified. The impact of substrate concentration, relative to substrate K_M, for each common inhibition modality is also reviewed. The pattern of changes in IC₅₀ that accompany increasing substrate concentration are shown to be diagnostic of specific inhibition modalities. Thus, replots of IC₅₀ as a function of the ratio [S]/K_M are recommended as a simple and rapid means of assessing inhibition modality. Finally, specific recommendations are offered for ideal experimental conditions for the determination of inhibition modality through the use of IC₅₀ replots.

Rational design and characterization of D-Phe-Pro-D-Arg-derived direct thrombin inhibitors

The tremendous social and economic impact of thrombotic disorders, together with the considerable risks associated to the currently available therapies, prompt for the development of more efficient and safer anticoagulants. Novel peptide-based thrombin inhibitors were identified using in silico structure-based design and further validated in vitro. The best candidate compounds contained both l- and d-amino acids, with the general sequence d-Phe(P3)-Pro(P2)-d-Arg(P1)-P1′-CONH₂. The P1′ position was scanned with l- and d-isomers of natural or unnatural amino acids, covering the major chemical classes. The most potent non-covalent and proteolysis-resistant inhibitors contain small hydrophobic or polar amino acids (Gly, Ala, Ser, Cys, Thr) at the P1′ position. The lead tetrapeptide, d-Phe-Pro-d-Arg-d-Thr-CONH₂, competitively inhibits α-thrombin's cleavage of the S2238 chromogenic substrate with a K_i of 0.92 µM. In order to understand the molecular details of their inhibitory action, the three-dimensional structure of three peptides (with P1′ l-isoleucine (fPrI), l-cysteine (fPrC) or d-threonine (fPrt)) in complex with human α-thrombin were determined by X-ray crystallography. All the inhibitors bind in a substrate-like orientation to the active site of the enzyme. The contacts established between the d-Arg residue in position P1 and thrombin are similar to those observed for the l-isomer in other substrates and inhibitors. However, fPrC and fPrt disrupt the active site His57-Ser195 hydrogen bond, while the combination of a P1 d-Arg and a bulkier P1′ residue in fPrI induce an unfavorable geometry for the nucleophilic attack of the scissile bond by the catalytic serine. The experimental models explain the observed relative potency of the inhibitors, as well as their stability to proteolysis. Moreover, the newly identified direct thrombin inhibitors provide a novel pharmacophore platform for developing antithrombotic agents by exploring the conformational constrains imposed by the d-stereochemistry of the residues at positions P1 and P1′.

Reproducing crystal binding modes of ligand functional groups using Site-Identification by Ligand Competitive Saturation (SILCS) simulations

The applicability of a computational method, Site Identification by Ligand Competitive Saturation (SILCS), to identify regions on a protein surface with which different types of functional groups on low-molecular weight inhibitors interact is demonstrated. The method involves molecular dynamics (MD) simulations of a protein in an aqueous solution of chemically diverse small molecules from which probability distributions of fragments types, termed FragMaps, are obtained. In the present application, SILCS simulations are performed with an aqueous solution of 1 M benzene and propane to map the affinity pattern of the protein for aromatic and aliphatic functional groups. In addition, water hydrogen and oxygen atoms serve as probes for hydrogen-bond donor and acceptor affinity, respectively. The method is tested using a set of 7 proteins for which crystal structures of complexes with several high affinity inhibitors are known. Good agreement is obtained between FragMaps and the positions of chemically similar functional groups in inhibitors as observed in the X-ray crystallographic structures. Quantitative capabilities of the SILCS approach are demonstrated by converting FragMaps to free energies, termed Grid Free Energies (GFE), and showing correlation between the GFE values and experimental binding affinities. For proteins for which ligand decoy sets are available, GFE values are shown to typically score the crystal conformation and conformations similar to it more favorable than decoys. Additionally, SILCS is tested for its ability to capture the subtle differences in ligand affinity across homologous proteins, information which may be of utility toward specificity-guided drug design. Taken together, our results show that SILCS can recapitulate the known location of functional groups of bound inhibitors for a number of proteins, suggesting that the method may be of utility for rational drug design.

Mapping the energetics of water-protein and water-ligand interactions with the "natural" HINT forcefield: predictive tools for characterizing the roles of water in biomolecules

The energetics and hydrogen bonding pattern of water molecules bound to proteins were mapped by analyzing structural data (resolution better than 2.3A) for sets of uncomplexed and ligand-complexed proteins. Water-protein and water-ligand interactions were evaluated using hydropatic interactions (HINT), a non-Newtonian forcefield based on experimentally determined logP(octanol/water) values. Potential water hydrogen bonding ability was assessed by a new Rank algorithm. The HINT-derived binding energies and Ranks for second shell water molecules were -0.04 kcal mol(-1) and 0.0, respectively, for first shell water molecules -0.38 kcal mol(-1) and 1.6, for active site water molecules -0.45 kcal mol(-1) and 2.3, for cavity water molecules -0.55 kcal mol(-1) and 3.3, and for buried water molecules -0.56 kcal mol(-1) and 4.4. For the last four classes, similar energies indicate that internal and external water molecules interact with protein almost equally, despite different degrees of hydrogen bonding. The binding energies and Ranks for water molecules bridging ligand-protein were -1.13 kcal mol(-1) and 4.5, respectively. This energetic contribution is shared equally between protein and ligand, whereas Rank favors the protein. Lastly, by comparing the uncomplexed and complexed forms of proteins, guidelines were developed for prediction of the roles played by active site water molecules in ligand binding. A water molecule with high Rank and HINT score is unlikely to make further interactions with the ligand and is largely irrelevant to the binding process, while a water molecule with moderate Rank and high HINT score is available for ligand interaction. Water molecule displaced for steric reasons were characterized by lower Rank and HINT score. These guidelines, tested by calculating HINT score and Rank for 50 water molecules bound in the active site of four uncomplexed proteins (for which the structures of the liganded forms were also available), correctly predicted the ultimate roles (in the complex) for 76% of water molecules. Some failures were likely due to ambiguities in the structural data.

Virtual Screening and Scaffold Hopping Based on GRID Molecular Interaction Fields

In this study, a set of strategies for structure-based design using GRID molecular interaction fields (MIFs) to derive a pharmacophoric representation of a protein is reported. Thrombin, one of the key enzymes involved in the blood coagulation cascade, was chosen as the model system since abundant published experimental data are available related to both crystal structures and structurally diverse sets of inhibitors. First, a virtual screening methodology was developed either using a pharmacophore representation of the protein based on GRID MIFs or using GRID MIFs from the 3D structure of a set of chosen thrombin inhibitors. The search was done in a 3D multiconformation version of the Available Chemical Directory (ACD) database, which had been spiked with 262 known thrombin inhibitors (multiple conformers available per compound). The model managed to find 80% of the known thrombin inhibitors among the 74,291 conformers in the ACD by only searching 5% of the database; hence, a 15-fold enrichment of the library was achieved. Second, a scaffold hopping methodology was developed using GRID MIFs, giving the scaffold interaction pattern and the shape of the scaffold, together with the distance between the anchor points. The scaffolds reported by Dolle in the Journal of Combinatorial Chemistry summaries (2000 and 2001) and scaffolds built or derived from ligands cocomplexed with the thrombin enzyme were parameterized using a new set of descriptors and saved into a searchable database. The scaffold representation from the database was then compared to a template scaffold (from a thrombin crystal structure), and the thrombin-derived scaffolds included in the database were found among the top solutions. To validate the usefulness of the methodology to replace the template scaffold, the entire molecule was built (scaffold and side chains) and the resulting compounds were docked into the active site of thrombin. The docking solutions showed the same binding pattern as the cocomplexed compound, hence, showing that this method can be a valuable tool for medicinal chemists to select interchangeable core structures (scaffolds) in an easy manner and retaining the binding properties from the original ligand.

Enzyme Binding Selectivity Prediction: α-Thrombin vs Trypsin Inhibition

In the present work we explore the possibility of an in-depth computational analysis of available experimental X-ray structures in the specific case of a series of alpha-thrombin and trypsin complexes with their respective inhibitors for the development of a novel scoring function based on molecular electrostatic potential computed at the contact surface in the enzyme-inhibitor molecular complex. We subsequently employ the chemometrical approach to determine which are the interactions in the large volume of data that determine the resulting experimental binding constant between ligand and receptor. The results of the model evaluated with molecules in the independent validation set show that a reasonable average error of 1.30 log units of the difference between experimental and calculated binding constants was achieved in the system thrombin-trypsin, which is comparable with those of methods from the literature. Furthermore, by a careful preparation of the Kohonen top layer in the artificial neural network approach that is normally perceived as a "black box device", we have been able to follow the implications of the structure of the inhibitor-enzyme complex for the inhibitor's binding constant. The method appears to be suitable for evaluation of selectivity in structurally similar enzymatic systems, which is currently an important problem in drug design.

Design of weakly basic thrombin inhibitors incorporating novel P1 binding functions: molecular and X-ray crystallographic studies

To prepare weakly basic thrombin inhibitors with modified S1 anchoring groups, two series of compounds were synthesized by reaction of guanidine or aminoguanidine with acyl halides and N,N-disubstituted carbamoyl chlorides. pK(a) measurements of these acylated guanidines/aminoguanidines showed a reduced basicity, with pK(a) values in the range of 8.4-8.7. These molecules typically showed inhibition constants in the range of 150-425 nM against thrombin and 360-965 nM against trypsin, even though some bulky derivatives, such as N,N-diphenylcarbamoylguanidine/aminoguanidine and their congeners, showed much stronger thrombin inhibitory activity, with inhibition constants in the range of 24-42 nM. Unexpectedly, very long incubation times with both proteases revealed that aminoguanidine derivatives behaved as irreversible inhibitors. To assess the molecular basis responsible for the high affinity observed for these molecules toward thrombin, the crystal structure of the thrombin-hirugen-N,N-diphenylcarbamoylaminoguanidine complex has been solved at 1.90 A resolution. The structural analysis of the complex revealed an unexpected interaction mode with the protease, resulting in an N,N-diphenylcarbamoyl intermediate covalently bound to the catalytic serine as a consequence of its hydrolysis together with the release of the aminoguanidine moiety. Surprisingly, in this covalent adduct a phenyl group was found in the S1 specificity pocket, which usually recognizes positively charged residues. These findings provide new insights in the design of low basicity serine protease inhibitors.

A theoretical study of conformational properties of N-methyl azapeptide derivatives

The conformational properties of azapeptide derivatives, Ac-azaGly-NHMe (1), Ac-azaAla-NHMe (2), Ac-NMe-azaGly-NHMe (3), Ac-NMe-azaAla-NHMe (4), Ac-azaGly-NMe(2) (5), Ac-azaAla-NMe(2) (6), Ac-NMe-azaGly-NMe(2) (7), and Ac-NMe-azaAla-NMe(2) (8), were systematically examined by using ab initio MO and DFT methods. Structural perturbations in azapeptides resulting from cyclic substitution of a methyl group at three N-positions of an azaamino acid were studied on the basis of the structure of the simplest model azapeptide, 1. Potential energy surfaces were generated at the HF/6-31G level for 1-4 and at the HF/6-31G//HF/3-21G level for 5-8 by rotating two key dihedral angles (phi, psi) in increments of 30 degrees. The backbone (phi, psi) angles of the minima for 1-4 are observed at the i + 2 position to form the betaI(I')-, betaII(II')-, betaVI-turns or the polyproline II structure according to the orientation of the acetyl group and the positions of the N-methyl groups. Compounds 5-8 coupled to a secondary amine were found to preferentially adopt polyproline II, betaI(III)-turn, or alpha-helical structure or even extended conformations depending on the orientation of the acetyl group and the positions of the N-methyl groups. Furthermore, N-methyl groups, depending on their positions, were found to affect the orientation of the amide group in the lowest energy conformations, the pyramidality of the N2 atom, and the bond length in azapeptide derivatives. These unique theoretical conformations of N-methyl azapeptide derivatives could be utilized in the definite design of secondary structure for peptides and proteins, and in the development of new drugs and molecular machines.

Serine proteinase inhibition by the active site titrant N alpha-(N, N-dimethylcarbamoyl)-alpha-azaornithine p-nitrophenyl ester. A comparative study

Kinetics for the hydrolysis of the chromogenic active-site titrant N alpha-(N,N-dimethylcarbamoyl)-alpha-azaornithine p-nitrophenyl ester (Dmc-azaOrn-ONp) catalysed by bovine beta-trypsin, bovine alpha-thrombin, bovine Factor Xa, human alpha-thrombin, human Factor Xa, human Lys77-plasmin, human urinary kallikrein, Mr 33 000 and Mr 54 000 species of human urokinase, porcine pancreatic beta-kallikrein-A and -B and Ancrod (the coagulating serine proteinase from the Malayan pit viper Agkistrodon rhodostoma venom) have been obtained between pH 6.0 and 8.0, at 21.0 degrees C, and analysed in parallel with those for the enzymatic cleavage of N alpha-(N,N-dimethylcarbamoyl)-alpha-azalysine p-nitrophenyl ester (Dmc-azaLys-ONp). The enzyme kinetics are consistent with the minimum three-step catalytic mechanism of serine proteinases, the rate-limiting step being represented by the deacylation process. Bovine beta-trypsin kinetics are modulated by the acid-base equilibrium of the His57 catalytic residue (pKa approximately 6.9). Dmc-azaOrn-ONp and Dmc-azaLys-ONp bind stoichiometrically to the serine proteinase active site, and allow the reliable determination of the active enzyme concentration between 1.0 x 10-6 M and 3.0 x 10-4 M. The affinity and the reactivity for Dmc-azaOrn-ONp (expressed by Ks and k+2/Ks, respectively) of the serine proteinases considered are much lower than those for Dmc-azaLys-ONp. The very different affinity and reactivity properties for Dmc-azaOrn-ONp and Dmc-azaLys-ONp have been related to the different size of the ornithine/lysine side chains, and to the ensuing different positioning of the active-site titrants upon binding to the enzyme catalytic centre (i.e. to P1-S1 recognition). These data represent the first detailed comparative investigation on the catalytic properties of serine proteinases towards an ornithine derivative (i. e. Dmc-azaOrn-ONp).

Novel natural product 5,5-trans-lactone inhibitors of human alpha-thrombin: mechanism of action and structural studies

High-throughput screening of methanolic extracts from the leaves of the plant Lantana camara identified potent inhibitors of human alpha-thrombin, which were shown to be 5,5-trans-fused cyclic lactone euphane triterpenes [O'Neill et al. (1998) J. Nat. Prod. (submitted for publication)]. Proflavin displacement studies showed the inhibitors to bind at the active site of alpha-thrombin and alpha-chymotrypsin. Kinetic analysis of alpha-thrombin showed tight-binding reversible competitive inhibition by both compounds, named GR133487 and GR133686, with respective kon values at pH 8.4 of 1.7 x 10(6) s-1 M-1 and 4.6 x 10(6) s-1 M-1. Electrospray ionization mass spectrometry of thrombin/inhibitor complexes showed the tight-bound species to be covalently attached, suggesting acyl-enzyme formation by reaction of the active-site Ser195 with the trans-lactone carbonyl. X-ray crystal structures of alpha-thrombin/GR133686 (3.0 A resolution) and alpha-thrombin/GR133487 (2.2 A resolution) complexes showed continuous electron density between Ser195 and the ring-opened lactone carbonyl, demonstrating acyl-enzyme formation. Turnover of inhibitor by alpha-thrombin was negligible and mass spectrometry of isolated complexes showed that reversal of inhibition occurs by reformation of the trans-lactone from the acyl-enzyme. The catalytic triad appears undisrupted and the inhibitor carbonyl occupies the oxyanion hole, suggesting the observed lack of turnover is due to exclusion of water for deacylation. The acyl-enzyme inhibitor hydroxyl is properly positioned for nucleophilic attack on the ester carbonyl and therefore relactonization; furthermore, the higher resolution structure of alpha-thrombin/GR133487 shows this hydroxyl to be effectively superimposable with the recently proposed deacylating water for peptide substrate hydrolysis [Wilmouth, R. C., et al. (1997) Nat. Struct.Biol. 4, 456-462], suggesting the alpha-thrombin/GR133487 complex may be a good model for this reaction.