Highly selective mechanism-based thrombin inhibitors: structures of thrombin and trypsin inhibited with rigid peptidyl aldehydes

Novel inhibitors and activity-based probes targeting serine proteases

Serine proteases play varied and manifold roles in important biological, physiological, and pathological processes. These include viral, bacterial, and parasitic infection, allergic sensitization, tumor invasion, and metastasis. The use of activity-based profiling has been foundational in pinpointing the precise roles of serine proteases across this myriad of processes. A broad range of serine protease-targeted activity-based probe (ABP) chemotypes have been developed and we have recently introduced biotinylated and "clickable" peptides containing P1 N-alkyl glycine arginine N-hydroxy succinimidyl (NHS) carbamates as ABPs for detection/profiling of trypsin-like serine proteases. This present study provides synthetic details for the preparation of additional examples of this ABP chemotype, which function as potent irreversible inhibitors of their respective target serine protease. We describe their use for the activity-based profiling of a broad range of serine proteases including trypsin, the trypsin-like protease plasmin, chymotrypsin, cathepsin G, and neutrophil elastase (NE), including the profiling of the latter protease in clinical samples obtained from patients with cystic fibrosis.

Assessing the Performance of Mixed Strategies To Combine Lipophilic Molecular Similarity and Docking in Virtual Screening

The accuracy of structure-based (SB) virtual screening (VS) is heavily affected by the scoring function used to rank a library of screened compounds. Even in cases where the docked pose agrees with the experimental binding mode of the ligand, the limitations of current scoring functions may lead to sensible inaccuracies in the ability to discriminate between actives and inactives. In this context, the combination of SB and ligand-based (LB) molecular similarity may be a promising strategy to increase the hit rates in VS. This study explores different strategies that aim to exploit the synergy between LB and SB methods in order to mitigate the limitations of these techniques, and to enhance the performance of VS studies by means of a balanced combination between docking scores and three-dimensional (3D) similarity. Particularly, attention is focused to the use of measurements of molecular similarity with PharmScreen, which exploits the 3D distribution of atomic lipophilicity determined from quantum mechanical-based continuum solvation calculations performed with the MST model, in conjunction with three docking programs: Glide, rDock, and GOLD. Different strategies have been explored to combine the information provided by docking and similarity measurements for re-ranking the screened ligands. For a benchmarking of 44 datasets, including 41 targets, the hybrid methods increase the identification of active compounds, according to the early (ROCe%) and total (AUC) enrichment metrics of VS, compared to pure LB and SB methods. Finally, the hybrid approaches are also more effective in enhancing the chemical diversity of active compounds. The datasets employed in this work are available in https://github.com/Pharmacelera/Molecular-Similarity-and-Docking.

Crystal structure of highly glycosylated human leukocyte elastase in complex with an S2' site binding inhibitor

Glycosylated human leukocyte elastase (HLE) was crystallized and structurally analysed in complex with a 1,3-thiazolidine-2,4-dione derivative that had been identified as an HLE inhibitor in preliminary studies. In contrast to previously described HLE structures with small-molecule inhibitors, in this structure the inhibitor does not bind to the S1 and S2 substrate-recognition sites; rather, this is the first HLE structure with a synthetic inhibitor in which the S2' site is blocked that normally binds the second side chain at the C-terminal side of the scissile peptide bond in a substrate protein. The inhibitor also induces the formation of crystalline HLE dimers that block access to the active sites and that are also predicted to be stable in solution. Neither such HLE dimers nor the corresponding crystal packing have been observed in previous HLE crystal structures. This novel crystalline environment contributes to the observation that comparatively large parts of the N-glycan chains of HLE are defined by electron density. The final HLE structure contains the largest structurally defined carbohydrate trees among currently available HLE structures.

A potent, proteolysis-resistant inhibitor of kallikrein-related peptidase 6 (KLK6) for cancer therapy, developed by combinatorial engineering

Human tissue kallikrein (KLK) proteases are hormone-like signaling molecules with important functions in cancer pathophysiology. KLK-related peptidase 6 (KLK6), specifically, is highly up-regulated in several types of cancer, where its increased activity promotes cancer invasion and metastasis. This characteristic suggests KLK6 as an attractive target for therapeutic interventions. However, inhibitors that specifically target KLK6 have not yet been reported, possibly because KLK6 shares a high sequence homology and structural similarity with other serine proteases and resists inhibition by many polypeptide inhibitors. Here, we present an innovative combinatorial approach to engineering KLK6 inhibitors via flow cytometry-based screening of a yeast-displayed mutant library of the human amyloid precursor protein Kunitz protease inhibitor domain (APPI), an inhibitor of other serine proteases, such as anionic and cationic trypsins. On the basis of this screening, we generated APPI_{M17L,I18F,S19F,F34V} (APPI-4M), an APPI variant with a KLK6 inhibition constant (K_i ) of 160 pm and a turnover time of 10 days. To the best of our knowledge, APPI-4M is the most potent KLK6 inhibitor reported to date, displaying 146-fold improved affinity and 13-fold improved proteolytic stability compared with WT APPI (APPI_WT). We further demonstrate that APPI-4M acts as a functional inhibitor in a cell-based model of KLK6-dependent breast cancer invasion. Finally, the crystal structures of the APPI_WT/KLK6 and APPI-4M/KLK6 complexes revealed the structural and mechanistic bases for the improved KLK6 binding and proteolytic resistance of APPI-4M. We anticipate that APPI-4M will have substantial translational potential as both imaging agent and therapeutic.

Factor XIa Inhibitors as New Anticoagulants

With the introduction of thrombin and factor Xa inhibitors to the oral anticoagulant market, significant improvements in both efficacy and safety have been achieved. Early clinical and preclinical data suggest that inhibitors of factor XIa can provide a still safer alternative, with expanded efficacy for arterial indications. This Perspective provides an overview of target rationale and details of the discovery and development of inhibitors of factor XIa as next generation antithrombotic agents.

Statistical potential for modeling and ranking of protein-ligand interactions

Applications in structural biology and medicinal chemistry require protein-ligand scoring functions for two distinct tasks: (i) ranking different poses of a small molecule in a protein binding site and (ii) ranking different small molecules by their complementarity to a protein site. Using probability theory, we developed two atomic distance-dependent statistical scoring functions: PoseScore was optimized for recognizing native binding geometries of ligands from other poses and RankScore was optimized for distinguishing ligands from nonbinding molecules. Both scores are based on a set of 8,885 crystallographic structures of protein-ligand complexes but differ in the values of three key parameters. Factors influencing the accuracy of scoring were investigated, including the maximal atomic distance and non-native ligand geometries used for scoring, as well as the use of protein models instead of crystallographic structures for training and testing the scoring function. For the test set of 19 targets, RankScore improved the ligand enrichment (logAUC) and early enrichment (EF(1)) scores computed by DOCK 3.6 for 13 and 14 targets, respectively. In addition, RankScore performed better at rescoring than each of seven other scoring functions tested. Accepting both the crystal structure and decoy geometries with all-atom root-mean-square errors of up to 2 Å from the crystal structure as correct binding poses, PoseScore gave the best score to a correct binding pose among 100 decoys for 88% of all cases in a benchmark set containing 100 protein-ligand complexes. PoseScore accuracy is comparable to that of DrugScore(CSD) and ITScore/SE and superior to 12 other tested scoring functions. Therefore, RankScore can facilitate ligand discovery, by ranking complexes of the target with different small molecules; PoseScore can be used for protein-ligand complex structure prediction, by ranking different conformations of a given protein-ligand pair. The statistical potentials are available through the Integrative Modeling Platform (IMP) software package (http://salilab.org/imp) and the LigScore Web server (http://salilab.org/ligscore/).

Molecular insights into substrate specificity of prostate specific antigen through structural modeling

Prostate Specific Antigen's (PSA) role as a biomarker for prostate cancer is well established but the physiological role of its serine protease activity in the pathobiology of normal prostate and prostate carcinogenesis remains largely unknown. In light of recent studies that implicate PSA's enzymatic activity in the initiation and/or progression of prostate cancer, we performed a molecular modeling study of substrate binding at the catalytic site of PSA wherein a PSA-selective substrate (HSSKLQ) was docked in an acyl-enzyme conformation to a three-dimensional homology model of PSA. Additionally, virtual positional scanning studies were conducted to gain mechanistic insights into substrate recognition of PSA. Subsequently, 13 novel peptide substrates of 6-aa length and four peptide substrates with varying length were synthesized and assayed for PSA hydrolysis to evaluate the experimental validity of docking insights. Additionally, six novel aldehyde-containing transition state analog inhibitors were synthesized and tested for their inhibitory potencies. The experimental data on the hydrolysis rates of the newly synthesized substrates and inhibitory potencies of the aldehyde peptides agreed with the docking predictions, providing validation of the docking methodology and demonstrating its utility towards the design of substrate-mimetic inhibitors that can be used to explore PSA's role in the pathobiology of prostate cancer.

Using peptidyl aldehydes in activity-based proteomics

The broad inhibitory spectrum of aldehydes and the possibility that amino acid residues modulate their specificity point to the potential of using peptidyl aldehydes as activity-based probes. Here, we establish the potential of peptidyl aldehydes in activity-based proteomics by synthesizing different probes and using them to specifically label a well-known serine protease in an activity-dependent manner. From our results, peptidyl aldehydes emerge as promising activity-based probes that enable multiple enzymatic-class detection by substrate recognition and can be used in diverse activity-based proteomics applications like protein identification and activity profiling.

Active and Exo-site Inhibition of Human Factor Xa: Structure of des-Gla Factor Xa Inhibited by NAP5, a Potent Nematode Anticoagulant Protein from Ancylostoma caninum

Hookworms are hematophagous nematodes capable of growth, development and subsistence in living host systems such as humans and other mammals. Approximately one billion, or one in six, people worldwide are infected by hookworms causing gastrointestinal blood loss and iron deficiency anemia. The hematophagous hookworm Ancylostoma caninum produces a family of small, disulfide-linked protein anticoagulants (75-84 amino acid residues). One of these nematode anticoagulant proteins, NAP5, inhibits the amidolytic activity of factor Xa (fXa) with K(i)=43 pM, and is the most potent natural fXa inhibitor identified thus far. The crystal structure of NAP5 bound at the active site of gamma-carboxyglutamic acid domainless factor Xa (des-fXa) has been determined at 3.1 A resolution, which indicates that Asp189 (fXa, S1 subsite) binds to Arg40 (NAP5, P1 site) in a mode similar to that of the BPTI/trypsin interaction. However, the hydroxyl group of Ser39 of NAP5 additionally forms a hydrogen bond (2.5 A) with His57 NE2 of the catalytic triad, replacing the hydrogen bond of Ser195 OG to the latter in the native structure, resulting in an interaction that has not been observed before. Furthermore, the C-terminal extension of NAP5 surprisingly interacts with the fXa exosite of a symmetry-equivalent molecule forming a short intermolecular beta-strand as observed in the structure of the NAPc2/fXa complex. This indicates that NAP5 can bind to fXa at the active site, or the exosite, and to fX at the exosite. However, unlike NAPc2, NAP5 does not inhibit fVIIa of the fVIIa/TF complex.

Classification and comparison of ligand-binding sites derived from grid-mapped knowledge-based potentials

We describe the application of knowledge-based potentials implemented in the MOE program to compare the ligand-binding sites of several proteins. The binding probabilities for a polar and a hydrophobic probe are calculated on a grid to allow easy comparison of binding sites of superimposed related proteins. The method is fast and simple enough to simultaneously use structural information of multiple proteins of a target family. The method can be used to rapidly cluster proteins into subfamilies according to the similarity of hydrophobic and polar fields of their ligand-binding sites. Regions of the binding site which are common within a protein family can be identified and analysed for the design of family-targeted libraries or those which differ for improvement of ligand selectivity. The field-based hierarchical clustering is demonstrated for three protein families: the ligand-binding domains of nuclear receptors, the ATP-binding sites of protein kinases and the substrate binding sites of proteases. More detailed comparisons are presented for serine proteases of the chymotrypsin family, for the peroxisome proliferator-activated receptor subfamily of nuclear receptors and for progesterone and androgen receptor. The results are in good accordance with structure-based analysis and highlight important differences of the binding sites, which have been also described in the literature.

Phenolic P2/P3 core motif as thrombin inhibitors—Design, synthesis, and X-ray co-crystal structure

Prototypical thrombin inhibitors were synthesized based on a trisubstituted phenol as a core motif. A naphthylsulfonamide analogue showed excellent antithrombin activity. An X-ray co-crystal structure showed the expected interactions.

Physics-Based Scoring of Protein−Ligand Complexes: Enrichment of Known Inhibitors in Large-Scale Virtual Screening

We demonstrate that using an all-atom molecular mechanics force field combined with an implicit solvent model for scoring protein-ligand complexes is a promising approach for improving inhibitor enrichment in the virtual screening of large compound databases. The rescoring method is evaluated by the extent to which known binders for nine diverse, therapeutically relevant enzymes are enriched against a background of approximately 100,000 drug-like decoys. The improvement in enrichment is most robust and dramatic within the top 1% of the ranked database, that is, the first thousand compounds; below the first few percent of the ranked database, there is little overall improvement. The improved early enrichment is likely due to the more realistic treatment of ligand and receptor desolvation in the rescoring procedure. We also present anecdotal but encouraging results assessing the ability of the rescoring method to predict specificity of inhibitors for structurally related proteins.

Synthesis and SAR of thrombin inhibitors incorporating a novel 4-amino-morpholinone sscaffold: analysis of X-ray crystal structure of enzyme inhibitor complex

A 4-amino-2-carboxymethyl-3-morpholinone structural motif derived from malic acid has been used to mimic d-Phe-Pro in the thrombin inhibiting tripeptide d-Phe-Pro-Arg. The arginine in D-Phe-Pro-Arg was replaced by the more rigid P1 truncated p-amidinobenzylamine (Pab). These new thrombin inhibitors were used to probe the inhibitor binding site of alpha-thrombin. The best candidate in this series of thrombin inhibitors exhibits an in vitro IC50 of 0.130 microM. Interestingly, the stereochemistry of the 4-amino-2-carboxymethyl-3-morpholinone motif is reversed for the most active compounds compared to that of a previously reported 2-carboxymethyl-3-morpholinone series. The X-ray crystal structure of the lead inhibitor cocrystallized with alpha-thrombin is discussed.

Non-covalent thrombin inhibitors featuring P3-Heterocycles with P1-Bicyclic arginine surrogates

Novel, potent, and highly selective classes of thrombin inhibitors were identified, which resulted from judicious combination of P4-aromatics and P2-P3-heterocyclic dipeptide surrogates with weakly basic (calcd pKa approximately non-basic-8.6) bicyclic P1-arginine mimics. The design, synthesis, and biological activity of achiral, non-covalent, orally bioavailable inhibitors NC1-NC44 featuring P1-indazoles, benzimidazoles, indoles, benzotriazoles, and aminobenzisoxazoles is disclosed.

A comparison of the pharmacophore identification programs: Catalyst, DISCO and GASP

Three commercially available pharmacophore generation programs, Catalyst/HipHop, DISCO and GASP, were compared on their ability to generate known pharmacophores deduced from protein-ligand complexes extracted from the Protein Data Bank. Five different protein families were included Thrombin, Cyclin Dependent Kinase 2, Dihydrofolate Reductase, HIV Reverse Transcriptase and Thermolysin. Target pharmacophores were defined through visual analysis of the data sets. The pharmacophore models produced were evaluated qualitatively through visual inspection and according to their ability to generate the target pharmacophores. Our results show that GASP and Catalyst outperformed DISCO at reproducing the five target pharmacophores.

Non-covalent thrombin inhibitors featuring P(3)-heterocycles with P(1)-monocyclic arginine surrogates

Investigations on P(2)-P(3)-heterocyclic dipeptide surrogates directed towards identification of an orally bioavailable thrombin inhibitor led us to pursue novel classes of achiral, non-covalent P(1)-arginine derivatives. The design, synthesis, and biological activity of inhibitors NC1-NC30 that feature three classes of monocyclic P(1)-arginine surrogates will be disclosed: (1) (hetero)aromatic amidines, amines and hydroxyamidines, (2) 2-aminopyrazines, and (3) 2-aminopyrimidines and 2-aminotetrahydropyrimidines.

Novel, potent non-covalent thrombin inhibitors incorporating p(3)-lactam scaffolds

Evolution of P(1)-argininal inhibitor prototypes led to a series of non-covalent P(3)-7-membered lactam inhibitors 1a-w, featuring novel peptidomimetic units that probe each of the S(1), S(2), and S(3) specificity pockets of thrombin. Rigid P(1)-arginine surrogates possessing a wide range of basicity (calcd pK(a)'s approximately neutral-14) were surveyed. The design, synthesis, and biological activity of these targets are presented.

An oxyanion-hole selective serine protease inhibitor in complex with trypsin

p-amidinophenylmethylphosphinic acid (AMPA) was designed, synthesized and crystallized in complex with trypsin to study interactions with the oxyanion hole at the S1 site. In comparison to benzamidine, AMPA shows improved activity, which the crystal structure demonstrates to result from hydrogen bonds between the negatively charged phosphinic acid group and the catalytic residues at the oxyanion hole.

New targets for antithrombotic drugs

The normal hemostatic process is initiated by disruption in the vascular continuity and exposure of the subendothelial components. Platelets adhere to subendothelium-bound von Willebrand factor via glycoprotein (GP) Ib complex. This initial interaction per se and the release of platelet agonists transduce signals that lead to the rise in intracellular Ca(2+). The rise in Ca(2+) induces shape change, prostaglandin synthesis, release of granular contents and conformational changes in platelet Gp IIb-IIIa. Gp IIb-IIIa in activated platelets becomes competent to bind fibrinogen and other adhesive proteins and mediates platelet cohesion (primary hemostatic plug). Furthermore, the activated platelet surface provides an efficient catalytic surface for the coagulation reactions, ultimately resulting in the formation of fibrin (secondary hemostasis). Normally the hemostatic process plays a delicate balance between keeping the blood in the fluid state to maintain flow and rapidly forming an occluding plug following vessel injury. Thrombosis occurs because of alteration in this delicate balance. Consequences of thrombosis are a major cause of morbidity and mortality in industrialized countries. Arterial thrombosis occurs in the setting of previous vessel wall injury mostly because of atherosclerosis, while venous thrombosis occurs in areas of stasis. The recent advances in our understanding of the hemostatic process have led to a better elucidation of the mechanism of action of many antithrombotic drugs and identification of new targets for drug development. The molecular target of the well known antiplatelet drug ticlopidine has been identified. Large numbers of IIb-IIIa inhibitors have been developed based on the crystal structure of a potent antagonist echistatin. The mechanism of action of heparin has been defined at the molecular level. As a result a synthetic pentasaccharide, based on antithrombin-binding domain of heparin, has been developed and tested successfully in clinical trials. New generation direct thrombin inhibitors are being developed based on the crystal structure of thrombin. Factor Xa has a critical position at the convergence of intrinsic and extrinsic pathway ways. The clinical tolerability and the efficacy of low molecular weight heparins led to the concept that inhibition of further thrombin generation, by blocking factor Xa alone, can be an effective way of preventing thrombus growth without inactivating thrombin. A large number of specific factor Xa inhibitors are under development. Some of these drugs have already undergone preliminary clinical trials and appear to be promising. Future clinical trials will determine whether these new drugs will provide better risk-benefit ratio in treatment of thrombotic disorders.

Prediction of enzyme binding: human thrombin inhibition study by quantum chemical and artificial intelligence methods based on X-ray structures

Thrombin is a serine protease which plays important roles in the human body, the key one being the control of thrombus formation. The inhibition of thrombin has become a target for new antithrombotics. The aim of our work was to (i) construct a model which would enable us to predict Ki values for the binding of an inhibitor into the active site of thrombin based on a database of known X-ray structures of inhibitor-enzyme complexes and (ii) to identify the structural and electrostatic characteristics of inhibitor molecules crucially important to their effective binding. To retain as much of the 3D structural information of the bound inhibitor as possible, we implemented the quantum mechanical/molecular mechanical (QM/MM) procedure for calculating the molecular electrostatic potential (MEP) at the van der Waals surfaces of atoms in the protein's active site. The inhibitor was treated quantum mechanically, while the rest of the complex was treated by classical means. The obtained MEP values served as inputs into the counter-propagation artificial neural network (CP-ANN), and a genetic algorithm was subsequently used to search for the combination of atoms that predominantly influences the binding. The constructed CP-ANN model yielded Ki values predictions with a correlation coefficient of 0.96, with Ki values extended over 7 orders of magnitude. Our approach also shows the relative importance of the various amino acid residues present in the active site of the enzyme for inhibitor binding. The list of residues selected by our automatic procedure is in good correlation with the current consensus regarding the importance of certain crucial residues in thrombin's active site.

Protease inhibitors: Part 4. Synthesis of weakly basic thrombin inhibitors incorporating pyridinium-sulfanilylaminoguanidine moieties

Three series of derivatives have been prepared by reaction of sulfanilylaminoguanidine with pyrylium salts, with the pyridinium derivatives of glycine and with the pyridinium derivatives of beta-alanine, respectively. The new compounds were assayed as inhibitors of two serine proteases, thrombin and trypsin. The study showed that in contrast to the leads, possessing KI's around 100-300 nM against thrombin, and 450-1420 nM against trypsin, respectively, the new derivatives showed inhibition constants in the range of 15-50 nM against thrombin, whereas their affinity for trypsin remained relatively low. Derivatives of beta-alanine were more active than the corresponding glycine derivatives, which in turn were more inhibitory than the pyridinium derivatives of sulfanilylaminoguanidine possessing the same substitution pattern at the pyridinium ring. Thus, the present study proposes two novel approaches for the preparation of high affinity, specific thrombin inhibitors: a novel S1 anchoring moiety in the already large family of arginine/amidine-based inhibitors, i.e., the SO2NHNHC(=NH)NH2 group, and novel non-peptidomimetic scaffolds obtained by incorporating alkyl-/aryl-substituted-pyridinium moieties in the hydrophobic binding site(s). The first one is important for obtaining bioavailable thrombin inhibitors, devoid of the high basicity of the commonly used arginine/amidine-based inhibitors, whereas the second one may lead to improved water solubility of such compounds.

Interaction of the factor XIII activation peptide with alpha -thrombin. Crystal structure of its enzyme-substrate analog complex

The serine protease thrombin proteolytically activates blood coagulation factor XIII by cleavage at residue Arg(37); factor XIII in turn cross-links fibrin molecules and gives mechanical stability to the blood clot. The 2.0-A resolution x-ray crystal structure of human alpha-thrombin bound to the factor XIII-(28-37) decapeptide has been determined. This structure reveals the detailed atomic level interactions between the factor XIII activation peptide and thrombin and provides the first high resolution view of this functionally important part of the factor XIII molecule. A comparison of this structure with the crystal structure of fibrinopeptide A complexed with thrombin highlights several important determinants of thrombin substrate interaction. First, the P1 and P2 residues must be compatible with the geometry and chemistry of the S1 and S2 specificity sites in thrombin. Second, a glycine in the P5 position is necessary for the conserved substrate conformation seen in both factor XIII-(28-37) and fibrinopeptide A. Finally, the hydrophobic residues, which occupy the aryl binding site of thrombin determine the substrate conformation further away from the catalytic residues. In the case of factor XIII-(28-37), the aryl binding site is shared by hydrophobic residues P4 (Val(34)) and P9 (Val(29)). A bulkier residue in either of these sites may alter the substrate peptide conformation.

Novel, potent and selective chimeric FXa inhibitors featuring hydrophobic P1-ketoamide moieties

Judicious combination of P-region sequences of highly potent anticoagulant proteins including NAP5, NAP6, Ecotin, and Antistasin with SAR from small molecule FXa inhibitors led to a series of chimeric inhibitors of formula 1a-j. We report herein the design, synthesis, and biological activity of this novel family of FXa inhibitors that express both high in vitro potency and superb selectivity against related serine proteases.

GRID/CPCA: a new computational tool to design selective ligands

We present a computational procedure aimed at understanding enzyme selectivity and guiding the design of drugs with respect to selectivity. It starts from a set of 3D structures of the target proteins characterized by the program GRID. In the multivariate description proposed, the variables are organized and scaled in a different way than previously published methodologies. Then, consensus principal component analysis (CPCA) is used to analyze the GRID descriptors, allowing the straightforward identification of possible modifications in the ligand to improve its selectivity toward a chosen target. As an important new feature the computational method is able to work with more than two target proteins and with several 3D structures for each protein. Additionally, the use of a 'cutout tool' allows to focus on the important regions around the active site. The method is validated for a total number of nine structures of the three homologous serine proteases thrombin, trypsin, and factor Xa. The regions identified by the method as being important for selectivity are in excellent agreement with available experimental data and inhibitor structure-activity relationships.

Protease inhibitors. Part 2. Weakly basic thrombin inhibitors incorporating sulfonyl-aminoguanidine moieties as S1 anchoring groups: synthesis and structure-activity correlations

Two series of derivatives have been prepared and assayed as inhibitors of two physiologically relevant serine proteases, human thrombin and human trypsin. The first series includes alkyl-/ aralkyl-/aryl- and hetarylsulfonyl-aminoguanidines. It was thus observed that sulfanilyl-aminoguanidine possesses moderate but intrinsically selective thrombin inhibitory properties, with KI values around 90 and 1400 nM against thrombin and trypsin respectively. Further elaboration of this molecule afforded compounds that inhibited thrombin with KI values in the range 10-50 nM, whereas affinity for trypsin remained relatively low. Such compounds were obtained either by attaching benzyloxycarbonyl- or 4-toluenesulfonylureido-protected amino acids (such as D-Phe, L-Pro) or dipeptides (such as Phe-Pro, Gly His, beta-Ala-His or Pro-Gly) to the N-4 atom of the lead molecule, sulfanilyl-aminoguanidine, or by attaching substituted-pyridinium propylcarboxamido moieties to this lead. Thus, this study brings novel insights regarding a novel non-basic S1 anchoring moiety (i.e., SO2NHNHC(=NH)NH2), and new types of peptidomimetic scaffolds obtained by incorporating tosylureido-amino acids/pyridinium-substituted-GABA moieties in the hydrophobic binding site(s). Structure-activity correlations of the new serine protease inhibitors are also discussed based on a QSAR model described previously for a large series of structurally-related derivatives (Supuran et al. (1999) J. Med. Chem., in press).

Protease inhibitors: synthesis and QSAR study of novel classes of nonbasic thrombin inhibitors incorporating sulfonylguanidine and O-methylsulfonylisourea moieties at P1

Using benzamidine as a lead molecule, two series of alkyl/aralkyl/arylsulfonylguanidines/sulfonyl-O-methylisoureas+ ++ have been prepared and assayed as inhibitors of two serine proteases, thrombin and trypsin. The study showed that sulfaguanidine and its corresponding O-methylisourea derivative possess moderate but intrinsically selective thrombin inhibitory properties, with K(I)'s around 100 nM against thrombin and 1350-1500 nM against trypsin. Further elaboration of these two molecules afforded compounds that inhibited thrombin with K(I)'s in the range of 12-50 nM, whereas affinity for trypsin remained relatively low. Such compounds were obtained by attaching benzyloxycarbonyl- or 4-toluenesulfonylureido-protected amino acids (such as L- and D-Phe or L-Pro) or dipeptides (such as Phe-Pro, Gly-His, beta-Ala-His, or Pro-Gly) to the two leads mentioned above, sulfaguanidine and 4-aminobenzenesulfonyl-O-methylisourea. Thus, the present study proposes two novel approaches for the preparation of high-affinity, specific thrombin inhibitors: two novel S1 anchoring moieties in the already large family of arginine/amidine-based inhibitors and novel peptidomimetic scaffolds obtained by incorporating tosylureido amino acids in the hydrophobic binding site(s). The first one is important for obtaining bioavailable thrombin inhibitors, devoid of the high basicity of the commonly used arginine/amidine-based inhibitors, whereas the second one may lead to improved water solubility of such compounds due to facilitated metal (sodium) salts formation (at the relatively acidic SO(2)NHCO protons) as well as increased stability at hydrolysis (in vivo). A QSAR study also explained the activity in terms of global properties of the molecules, electronic properties of the sulfonylguanidine/sulfonylisourea moiety, and novel descriptors, the frontier orbital phase angles (FOPA), that account for the directions of the nodes in the pi orbitals in the aromatic portion of those of the drugs in which the sulfonyl group was bound to a benzene ring. For thrombin inhibition, the size of the molecule was the dominant influence, while for trypsin inhibition the FOPA was the principal determinant of activity. The dependence of activity on the FOPA variables is perhaps the clearest example of a quantum effect in pharmacology and suggests a promising new tool for drug design.

Structural basis of the thrombin selectivity of a ligand that contains the constrained arginine mimic (2S)-2-amino-(3S)-3-(1-carbamimidoyl- piperidin-3-yl)-propanoic acid at P1

We have studied the thrombin and trypsin complexed structures of a pair of peptidomimetic thrombin inhibitors, containing different P1 fragments. The first has arginine as its P1 fragment, and the second contains the constrained arginine mimic (2S)-2-amino-(3S)-3-(1-carbamimidoyl-piperidin-3-yl)-propano ic acid (SAPA), a fragment known to enhance thrombin/trypsin selectivity of inhibitors. On the basis of an analysis of the nonbonded interactions present in the structures of the trypsin and thrombin complexes of the two inhibitors, the calculated accessible surfaces of the enzymes and inhibitors in the four complexes, data on known structures of trypsin complexes of inhibitors, and factor Xa inhibitory potency of these compounds, we conclude that the ability of this arginine mimic to increase thrombin selectivity of an inhibitor is mediated by its differential interaction with the residue at position 192 (chymotrypsinogen numbering). Thrombin has a glutamic acid at residue 192, and trypsin has a glutamine. The analysis also suggests that this constrained arginine mimic, when present in an inhibitor, might enhance selectivity against other trypsin-like enzymes that have a glutamine at residue position 192.

Design, synthesis and structure-activity relationship of a series of arginine aldehyde factor Xa inhibitors. Part 1: structures based on the (D)-Arg-Gly-Arg tripeptide sequence

A series of arginine aldehyde inhibitors was designed as transition state (TS) analogues based on the known factor Xa specific substrate Cbz-D-Arg-Gly-Arg-pNA. BnSO2-(D)Arg-Gly-Arg-H (20) was found to be the most potent and selective inhibitor of factor Xa and prothrombinase activity in this series.

Exploratory solid-phase synthesis of factor Xa inhibitors: discovery and application of p3-heterocyclic amides as novel types of non-basic arginine surrogates

A series of novel FXa inhibitors 2a-m and 3a-f was discovered that feature heterocyclic carboxamides tethered to a d-diaminobutyric acid sidechain. These neutral amide derivatives serve as novel P3 d-arginine mimics. Pyrazine carboxamide scaffolds afforded the most potent FXa inhibitors (e.g., 2b IC50 = 4.6 nM). The synthesis and biological activity of two focused libraries are reported.

Protease inhibitors - Part 3. Synthesis of non-basic thrombin inhibitors incorporating pyridinium-sulfanilylguanidine moieties at the P1 site

Using benzamidine and sulfaguanidine as lead molecules, three series of derivatives have been prepared by reaction of sulfaguanidine with pyrylium salts, with the pyridinium derivatives of glycine and with the pyridinium derivatives of beta-alanine, respectively. The new compounds were assayed as inhibitors of two serine proteases, thrombin and trypsin. The study showed that in contrast to the leads, possessing K(I)'s around 100-300 nM against thrombin, and 1200-1500 nM against trypsin, respectively, the new derivatives showed inhibition constants in the range of 15-50 nM against thrombin, whereas their affinity for trypsin remained relatively low. Derivatives of beta-alanine were more active than the corresponding Gly derivatives, which in turn were more inhibitory than the pyridinium derivatives of sulfaguanidine possessing the same substitution pattern at the pyridinium ring. Thus, the present study proposes two novel approaches for the preparation of high affinity, specific thrombin inhibitors: a novel S1 anchoring moiety in the already large family of arginine/amidine-based inhibitors, i.e., the SO(2)N=C(NH(2))(2) group, and novel non-peptidomimetic scaffolds obtained by incorporating alkyl-/aryl-substituted-pyridinium moieties in the hydrophobic binding site(s). The first one is important for obtaining bioavailable thrombin inhibitors, devoid of the high basicity of the commonly used arginine/amidine-based inhibitors, whereas the second one may lead to improved water solubility of such compounds due to facilitated salt formation as well as increased stability at hydrolysis (in vivo).

Asymmetric synthesis of novel quaternary alpha-hydroxy-delta-lactam dipeptide surrogates

Application of the Sharpless AD protocol to a series of alpha-(E)-benzylidene-delta-lactam precursors followed by selective deoxygenation provided efficient synthetic routes to the chiral quaternary alpha-hydroxy-gammalactam derivatives 4 and 5. These functionalized intermediates and the diol precursors 3 are regarded as novel types of D-Phe-Pro dipeptide surrogates that are useful as enzyme active site probes.

Novel benzo-fused lactam scaffolds as factor Xa inhibitors

Rigid benzolactam P3-P2 dipeptide mimics were designed and prepared as potential inhibitors of blood coagulation factor Xa. Methoxy substitution of the tetrahydrobenzazepinone scaffold led to potent and selective inhibitors. The synthesis and biological activities of these derivatives are reported herein.

Investigation of the S3 site of thrombin: design, synthesis and biological activity of 4-substituted 3-amino-2-pyridones incorporating P1-argininals

A novel scaffold for P4-P2 dipeptide mimics containing a rigid pyridone spacer was designed based on a virtual library strategy. Several selected nonpeptidic 4-aralkyl or 4-alkylpyridones incorporating a P1-argininal sequence were prepared. The modeling studies, synthesis and biological activities of these unique pyridone derivatives are reported herein.

Structure of extracellular tissue factor complexed with factor VIIa inhibited with a BPTI mutant

The event that initiates the extrinsic pathway of blood coagulation is the association of coagulation factor VIIa (VIIa) with its cell-bound receptor, tissue factor (TF), exposed to blood circulation following tissue injury and/or vascular damage. The natural inhibitor of the TF.VIIa complex is the first Kunitz domain of tissue factor pathway inhibitor (TFPI-K1). The structure of TF. VIIa reversibly inhibited with a potent (Ki=0.4 nM) bovine pancreatic trypsin inhibitor (BPTI) mutant (5L15), a homolog of TFPI-K1, has been determined at 2.1 A resolution. When bound to TF, the four domain VIIa molecule assumes an extended conformation with its light chain wrapping around the framework of the two domain TF cofactor. The 5L15 inhibitor associates with the active site of VIIa similar to trypsin-bound BPTI, but makes several unique interactions near the perimeter of the site that are not observed in the latter. Most of the interactions are polar and involve mutated positions of 5L15. Of the eight rationally engineered mutations distinguishing 5L15 from BPTI, seven are involved in productive interactions stabilizing the enzyme-inhibitor association with four contributing contacts unique to the VIIa.5L15 complex. Two additional unique interactions are due to distinguishing residues in the VIIa sequence: a salt bridge between Arg20 of 5L15 and Asp60 of an insertion loop of VIIa, and a hydrogen bond between Tyr34O of the inhibitor and Lys192NZ of the enzyme. These interactions were used further to model binding of TFPI-K1 to VIIa and TFPI-K2 to factor Xa, the principal activation product of TF.VIIa. The structure of the ternary protein complex identifies the determinants important for binding within and near the active site of VIIa, and provides cogent information for addressing the manner in which substrates of VIIa are bound and hydrolyzed in blood coagulation. It should also provide guidance in structure-aided drug design for the discovery of potent and selective small molecule VIIa inhibitors.

Alkoxide-catalyzed ring-opening of a novel homosaccharin derivative: synthesis of potent, selective P3-lactam thrombin inhibitors containing P4-o-alkoxycarbonylbenzylsulfonamide residues

A series of lactam derivatives 1b-g featuring P4-o-alkoxycarbonylbenzylsulfonamide residues along with the potential P4-homosaccharin prodrug candidate 1h was prepared in order to probe the thrombin S3 specificity pocket. The synthesis and alkoxide-catalyzed ring opening of the novel homosaccharin intermediate 7 followed by subsequent elaboration delivered the targets 1b-h which were potent and selective thrombin inhibitors. The design, synthesis, and biological activity of these targets will be presented.

Potent and selective thrombin inhibitors featuring hydrophobic, basic P3-P4-aminoalkyllactam moieties

Crystal structure and evolving SAR considerations of potent, selective benzylsulfonamide lactam thrombin inhibitors and related serine protease inhibitors have led to the design of novel thrombin inhibitors 1a-g, featuring hydrophobic, basic, P4-alkylaminolactam scaffolds that serve as novel types of P3-P4 dipeptide mimics. The design, synthesis, and biological activity of these targets is presented.

Design and construction of novel thrombin inhibitors featuring P3-P4 quaternary lactam dipeptide surrogates

Potent serine protease inhibitor 1a featuring a hybrid P3-P4 quaternary lactam dipeptide surrogate was prepared based upon SAR and molecular modeling investigations and in order to further probe the S2/S3 thrombin and FXa subsites. An efficient and concise synthetic route to the key aminolactam intermediate 4 was developed. The design, synthesis, and biological activity of this target and its P3-P4 diastereomer 1b is presented.