High-throughput screening of enzyme inhibitors: automatic determination of tight-binding inhibition constants

A community effort in SARS-CoV-2 drug discovery

The COVID-19 pandemic continues to pose a substantial threat to human lives and is likely to do so for years to come. Despite the availability of vaccines, searching for efficient small-molecule drugs that are widely available, including in low- and middle-income countries, is an ongoing challenge. In this work, we report the results of an open science community effort, the "Billion molecules against COVID-19 challenge", to identify small-molecule inhibitors against SARS-CoV-2 or relevant human receptors. Participating teams used a wide variety of computational methods to screen a minimum of 1 billion virtual molecules against 6 protein targets. Overall, 31 teams participated, and they suggested a total of 639,024 molecules, which were subsequently ranked to find 'consensus compounds'. The organizing team coordinated with various contract research organizations (CROs) and collaborating institutions to synthesize and test 878 compounds for biological activity against proteases (Nsp5, Nsp3, TMPRSS2), nucleocapsid N, RdRP (only the Nsp12 domain), and (alpha) spike protein S. Overall, 27 compounds with weak inhibition/binding were experimentally identified by binding-, cleavage-, and/or viral suppression assays and are presented here. Open science approaches such as the one presented here contribute to the knowledge base of future drug discovery efforts in finding better SARS-CoV-2 treatments.

DNA-encoded chemical libraries yield non-covalent and non-peptidic SARS-CoV-2 main protease inhibitors

The development of SARS-CoV-2 main protease (Mpro) inhibitors for the treatment of COVID-19 has mostly benefitted from X-ray structures and preexisting knowledge of inhibitors; however, an efficient method to generate Mpro inhibitors, which circumvents such information would be advantageous. As an alternative approach, we show here that DNA-encoded chemistry technology (DEC-Tec) can be used to discover inhibitors of Mpro. An affinity selection of a 4-billion-membered DNA-encoded chemical library (DECL) using Mpro as bait produces novel non-covalent and non-peptide-based small molecule inhibitors of Mpro with low nanomolar Ki values. Furthermore, these compounds demonstrate efficacy against mutant forms of Mpro that have shown resistance to the standard-of-care drug nirmatrelvir. Overall, this work demonstrates that DEC-Tec can efficiently generate novel and potent inhibitors without preliminary chemical or structural information.

Broad-spectrum coronavirus 3C-like protease peptidomimetic inhibitors effectively block SARS-CoV-2 replication in cells: Design, synthesis, biological evaluation, and X-ray structure determination

Despite the approval of vaccines, monoclonal antibodies and restrictions during the pandemic, the demand for new efficacious and safe antivirals is compelling to boost the therapeutic arsenal against the COVID-19. The viral 3-chymotrypsin-like protease (3CL^pro) is an essential enzyme for replication with high homology in the active site across CoVs and variants showing an almost unique specificity for Leu-Gln as P2-P1 residues, allowing the development of broad-spectrum inhibitors. The design, synthesis, biological activity, and cocrystal structural information of newly conceived peptidomimetic covalent reversible inhibitors are herein described. The inhibitors display an aldehyde warhead, a Gln mimetic at P1 and modified P2-P3 residues. Particularly, functionalized proline residues were inserted at P2 to stabilize the β-turn like bioactive conformation, modulating the affinity. The most potent compounds displayed low/sub-nM potency against the 3CL^pro of SARS-CoV-2 and MERS-CoV and inhibited viral replication of three human CoVs, i.e. SARS-CoV-2, MERS-CoV, and HCoV 229 in different cell lines. Particularly, derivative 12 exhibited nM-low μM antiviral activity depending on the virus, and the highest selectivity index. Some compounds were co-crystallized with SARS-CoV-2 3CL^pro validating our design. Altogether, these results foster future work toward broad-spectrum 3CL^pro inhibitors to challenge CoVs related pandemics.

Highly potent inhibitors of cathepsin K with a differently positioned cyanohydrazide warhead: structural analysis of binding mode to mature and zymogen-like enzymes

Cathepsin K (CatK) is a target for the treatment of osteoporosis, arthritis, and bone metastasis. Peptidomimetics with a cyanohydrazide warhead represent a new class of highly potent CatK inhibitors; however, their binding mechanism is unknown. We investigated two model cyanohydrazide inhibitors with differently positioned warheads: an azadipeptide nitrile Gü1303 and a 3-cyano-3-aza-β-amino acid Gü2602. Crystal structures of their covalent complexes were determined with mature CatK as well as a zymogen-like activation intermediate of CatK. Binding mode analysis, together with quantum chemical calculations, revealed that the extraordinary picomolar potency of Gü2602 is entropically favoured by its conformational flexibility at the nonprimed-primed subsites boundary. Furthermore, we demonstrated by live cell imaging that cyanohydrazides effectively target mature CatK in osteosarcoma cells. Cyanohydrazides also suppressed the maturation of CatK by inhibiting the autoactivation of the CatK zymogen. Our results provide structural insights for the rational design of cyanohydrazide inhibitors of CatK as potential drugs.

Potent macrocycle inhibitors of the human SAGA deubiquitinating module

The Spt-Ada-Gcn5 acetyltransferase (SAGA) transcriptional coactivator contains a four-protein subcomplex called the deubiquitinating enzyme (DUB) module that removes ubiquitin from histone H2B-K120. The human DUB module contains the catalytic subunit ubiquitin-specific protease 22 (USP22), which is overexpressed in a number of cancers that are resistant to available therapies. We screened a massive combinatorial library of cyclic peptides and identified potent inhibitors of USP22. The top hit was highly specific for USP22 compared with a panel of 44 other human DUBs. Cells treated with peptide had increased levels of H2B monoubiquitination, demonstrating the ability of the cyclic peptides to enter human cells and inhibit H2B deubiquitination. These macrocycle inhibitors are, to our knowledge, the first reported inhibitors of USP22/SAGA DUB module and show promise for development.

DNA-encoded chemistry technology yields expedient access to SARS-CoV-2 Mpro inhibitors

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has killed more than 4 million humans globally, but there is no bona fide Food and Drug Administration-approved drug-like molecule to impede the COVID-19 pandemic. The sluggish pace of traditional therapeutic discovery is poorly suited to producing targeted treatments against rapidly evolving viruses. Here, we used an affinity-based screen of 4 billion DNA-encoded molecules en masse to identify a potent class of virus-specific inhibitors of the SARS-CoV-2 main protease (M^pro) without extensive and time-consuming medicinal chemistry. CDD-1714, the initial three-building-block screening hit (molecular weight [MW] = 542.5 g/mol), was a potent inhibitor (inhibition constant [K _i] = 20 nM). CDD-1713, a smaller two-building-block analog (MW = 353.3 g/mol) of CDD-1714, is a reversible covalent inhibitor of M^pro (K _i = 45 nM) that binds in the protease pocket, has specificity over human proteases, and shows in vitro efficacy in a SARS-CoV-2 infectivity model. Subsequently, key regions of CDD-1713 that were necessary for inhibitory activity were identified and a potent (K _i = 37 nM), smaller (MW = 323.4 g/mol), and metabolically more stable analog (CDD-1976) was generated. Thus, screening of DNA-encoded chemical libraries can accelerate the discovery of efficacious drug-like inhibitors of emerging viral disease targets.

Accurate quantitative determination of affinity and binding kinetics for tight binding inhibition of xanthine oxidase

The accurate quantitative determination of affinity and binding kinetics (BK) for tight binding inhibition is extraordinary important from both the continuous optimization of compounds, particularly in developing structure-activity relationships (SAR), and the prediction of in vivo target occupancy (TO). Due to the unique properties for tight binding inhibition that the inhibitors are characterized by the ultrahigh-affinity, relatively fast association to the target enzyme combined with extremely slow dissociation of the inhibitor-enzyme binary complex, the classical steady state equilibrium methods are no longer valid. Here, we made several recommendations of how to design the optimal experiments and apply special mathematical calculation approaches to quantitatively evaluate the accurate affinity and BK as the examples of two tight binding inhibitors against the xanthine oxidase (XO), as well as compared the differences in the results calculated from the different data analytical methods and analyzed the influence of these differences on the XO engagement in human. Analysis of the results displayed that the accurate apparent dissociation constant (K_i^*,app) was 0.2 ± 0.06 nM for topiroxotstat and was 0.45 ± 0.2 nM for febuxostat; that on-rate (k_on) was (4.3 ± 1.1) × 10⁶ M^-1s^-1 for topiroxotstat and was(133.3 ± 3.5) × 10⁶ M^-1s^-1 for febuxostat, and off-rate (k_off) was (1.0±0.2) × 10^-5 s^-1 for topiroxotstat and was ≤ 0.16 × 10^-5 s^-1for febuxostat. Moreover, there were significant differences in the K_i^*,app and k_off values estimated using the appropriate specialized methods for tight binding inhibition versus classical steady state equilibrium methods, with the substantial differences of 14-fold and 32-fold reduction for topiroxostat, respectively, and of 9.6-fold and ≥ 213-fold reduction for febuxostat, while the k_on values remain the moderate differences for the two inhibitors. The obvious greater AUC of XO engagement time courses and longer durations of above 70% engagement by the appropriate specialized methods for tight binding inhibition were observed that the results display the differences of 70.1% and 88%, respectively for topiroxostat and of 38.1% and 35.0%, respectively for febuxostat in human liver cell than by classical steady state equilibrium methods. Again, our studies provide several valuable recommendations of the optimal experiment protocols and appropriate analytical approaches for accurately quantitatively assessing the affinity and BK parameters as well as demonstrate the ability of our recommended methods to generate reliable data for tight binding inhibitors against XO.

Reviewing the experimental and mathematical factors involved in tight binding inhibitors Ki values determination: The bi-functional protease inhibitor SmCI as a test model

Different methodologies for determining the dissociation equilibrium constant (K_i) of protein tight binding inhibitors are frequently found in the scientific literature. Taking into account that the K_i value is the main parameter characterizing the inhibition strength, its determination often represents the first step during the characterization of a potential drug. The purpose of this review is to summarize the current information related to tight binding inhibitors K_i values determination and discuss about the importance of different factors as the enzyme concentration, the inhibitor concentration dilution series, the enzyme-inhibitor incubation time and the dose-response data mathematical fitting. For this aim, the bi-functional SmCI protease inhibitor is used as a tool for exemplifying the experimental and mathematical steps performed during tight binding inhibitors K_i values determination. In addition, the natural and the different recombinant forms of SmCI were used to go deeply into the comparison of some mathematic approaches that are frequently used in the literature. Finally, other biochemical techniques that could be potentially used for tight binding inhibitors K_i values determination are also commented.

I36T↑T mutation in South African subtype C (C-SA) HIV-1 protease significantly alters protease-drug interactions

The efficacy of HIV-1 protease (PR) inhibition therapies is often compromised by the emergence of mutations in the PR molecule that reduces the binding affinity of inhibitors while maintaining viable catalytic activity and affinity for natural substrates. In the present study, we used a recombinant HIV-1 C-SA PR and a recently reported variant for inhibition (Ki, IC50) and thermodynamic studies against nine clinically used inhibitors. This is the first time that binding free energies for C-SA PR and the mutant are reported. This variant PR harbours a mutation and insertion (I36T↑T) at position 36 of the C-SA HIV-1 PR, and did not show a significant difference in the catalytic effect of the HIV-1 PR. However, the nine clinically approved HIV PR drugs used in this study demonstrated weaker inhibition and lower binding affinities toward the variant when compared to the wild type HIV-1 PR. All the protease inhibitors (PIs), except Amprenavir and Ritonavir exhibited a significant decrease in binding affinity (p<0.0001). Darunavir and Nelfinavir exhibited the weakest binding affinity, 155- and 95-fold decreases respectively, toward the variant. Vitality values for the variant PR, against the seven selected PIs, confirm the impact of the mutation and insertion on the South African HIV-1 subtype C PR. This information has important clinical implications for thousands of patients in Sub-Saharan Africa.

A fluorimetric readout reporting the kinetics of nucleotide-induced human ribonucleotide reductase oligomerization

Human ribonucleotide reductase (hRNR) is a target of nucleotide chemotherapeutics in clinical use. The nucleotide-induced oligomeric regulation of hRNR subunit α is increasingly being recognized as an innate and drug-relevant mechanism for enzyme activity modulation. In the presence of negative feedback inhibitor dATP and leukemia drug clofarabine nucleotides, hRNR-α assembles into catalytically inert hexameric complexes, whereas nucleotide effectors that govern substrate specificity typically trigger α-dimerization. Currently, both knowledge of and tools to interrogate the oligomeric assembly pathway of RNR in any species in real time are lacking. We therefore developed a fluorimetric assay that reliably reports on oligomeric state changes of α with high sensitivity. The oligomerization-directed fluorescence quenching of hRNR-α, covalently labeled with two fluorophores, allows for direct readout of hRNR dimeric and hexameric states. We applied the newly developed platform to reveal the timescales of α self-assembly, driven by the feedback regulator dATP. This information is currently unavailable, despite the pharmaceutical relevance of hRNR oligomeric regulation.

Quantitative characterization of the activation steps of mannan-binding lectin (MBL)-associated serine proteases (MASPs) points to the central role of MASP-1 in the initiation of the complement lectin pathway

Mannan-binding lectin (MBL)-associated serine proteases, MASP-1 and MASP-2, have been thought to autoactivate when MBL/ficolin·MASP complexes bind to pathogens triggering the complement lectin pathway. Autoactivation of MASPs occurs in two steps: 1) zymogen autoactivation, when one proenzyme cleaves another proenzyme molecule of the same protease, and 2) autocatalytic activation, when the activated protease cleaves its own zymogen. Using recombinant catalytic fragments, we demonstrated that a stable proenzyme MASP-1 variant (R448Q) cleaved the inactive, catalytic site Ser-to-Ala variant (S646A). The autoactivation steps of MASP-1 were separately quantified using these mutants and the wild type enzyme. Analogous mutants were made for MASP-2, and rate constants of the autoactivation steps as well as the possible cross-activation steps between MASP-1 and MASP-2 were determined. Based on the rate constants, a kinetic model of lectin pathway activation was outlined. The zymogen autoactivation rate of MASP-1 is ∼3000-fold higher, and the autocatalytic activation of MASP-1 is about 140-fold faster than those of MASP-2. Moreover, both activated and proenzyme MASP-1 can effectively cleave proenzyme MASP-2. MASP-3, which does not autoactivate, is also cleaved by MASP-1 quite efficiently. The structure of the catalytic region of proenzyme MASP-1 R448Q was solved at 2.5 Å. Proenzyme MASP-1 R448Q readily cleaves synthetic substrates, and it is inhibited by a specific canonical inhibitor developed against active MASP-1, indicating that zymogen MASP-1 fluctuates between an inactive and an active-like conformation. The determined structure provides a feasible explanation for this phenomenon. In summary, autoactivation of MASP-1 is crucial for the activation of MBL/ficolin·MASP complexes, and in the proenzymic phase zymogen MASP-1 controls the process.

Adenosine analogue-oligo-arginine conjugates (ARCs) serve as high-affinity inhibitors and fluorescence probes of type I cGMP-dependent protein kinase (PKGIalpha)

INTRODUCTION

Type I cGMP-dependent protein kinase (PKGIalpha) belongs to the family of cyclic nucleotide-dependent protein kinases and is one of the main effectors of cGMP. PKGIalpha is involved in regulation of cardiac contractility, vasorelaxation, and blood pressure; hence, the development of potent modulators of PKGIalpha would lead to advances in the treatment of a variety of cardiovascular diseases.

AIM

Representatives of ARC-type compounds previously characterized as potent inhibitors and high-affinity fluorescent probes of PKA catalytic subunit (PKAc) were tested towards PKGIalpha to determine that ARCs could serve as activity regulators and sensors for the latter protein kinase both in vitro and in complex biological systems.

RESULTS

Structure-activity profiling of ARCs with PKGIalpha in vitro demonstrated both similarities as well as differences to corresponding profiling with PKAc, whereas ARC-903 and ARC-668 revealed low nanomolar displacement constants and inhibition IC(50) values with both cyclic nucleotide-dependent kinases. The ability of ARC-based fluorescent probes to penetrate cell plasma membrane was demonstrated in the smooth muscle tissue of rat cerebellum isolated arteries, and the compound with the highest affinity in vitro (ARC-903) showed also potential for in vivo applications, fully abolishing the PKG1alpha-induced vasodilation.

Development and validation of a cytochrome c-coupled assay for pteridine reductase 1 and dihydrofolate reductase

Activity of the pterin- and folate-salvaging enzymes pteridine reductase 1 (PTR1) and dihydrofolate reductase-thymidylate synthetase (DHFR-TS) is commonly measured as a decrease in absorbance at 340 nm, corresponding to oxidation of nicotinamide adenine dinucleotide phosphate (NADPH). Although this assay has been adequate to study the biology of these enzymes, it is not amenable to support any degree of routine inhibitor assessment because its restricted linearity is incompatible with enhanced throughput microtiter plate screening. In this article, we report the development and validation of a nonenzymatically coupled screening assay in which the product of the enzymatic reaction reduces cytochrome c, causing an increase in absorbance at 550 nm. We demonstrate this assay to be robust and accurate, and we describe its utility in supporting a structure-based design, small-molecule inhibitor campaign against Trypanosoma brucei PTR1 and DHFR-TS.

One scaffold, three binding modes: novel and selective pteridine reductase 1 inhibitors derived from fragment hits discovered by virtual screening

The enzyme pteridine reductase 1 (PTR1) is a potential target for new compounds to treat human African trypanosomiasis. A virtual screening campaign for fragments inhibiting PTR1 was carried out. Two novel chemical series were identified containing aminobenzothiazole and aminobenzimidazole scaffolds, respectively. One of the hits (2-amino-6-chloro-benzimidazole) was subjected to crystal structure analysis and a high resolution crystal structure in complex with PTR1 was obtained, confirming the predicted binding mode. However, the crystal structures of two analogues (2-amino-benzimidazole and 1-(3,4-dichloro-benzyl)-2-amino-benzimidazole) in complex with PTR1 revealed two alternative binding modes. In these complexes, previously unobserved protein movements and water-mediated protein-ligand contacts occurred, which prohibited a correct prediction of the binding modes. On the basis of the alternative binding mode of 1-(3,4-dichloro-benzyl)-2-amino-benzimidazole, derivatives were designed and selective PTR1 inhibitors with low nanomolar potency and favorable physicochemical properties were obtained.

High-affinity bisubstrate probe for fluorescence anisotropy binding/displacement assays with protein kinases PKA and ROCK

The bisubstrate fluorescent probe ARC-583 (Adc-Ahx-(D-Arg)(6)-d-Lys(5-TAMRA)-NH2) and its application for the characterization of both ATP- and protein/peptide substrate-competitive inhibitors of protein kinases PKA (cyclic AMP-dependent protein kinase) and ROCK (rho kinase) in fluorescence polarization-based assay are described. High affinity of the probe (K(D)=0.48 nM toward PKA) enables its application for the characterization of inhibitors with nanomolar and micromolar potency and determination of the active concentration of the kinase in individual experiments as well as in the high-throughput screening format. The probe can be used for the assessment of protein-protein interactions (e.g., between regulatory and catalytic subunits of PKA) and as a cyclic AMP biosensor.

Bifunctional molecules evade cytochrome P(450) metabolism by forming protective complexes with FK506-binding protein

Despite their large size and complexity, the macrolide natural products rapamycin and FK506 have excellent pharmacological characteristics. We hypothesize that these unexpected properties may arise from protective, high affinity interactions with the cellular FK506-binding protein, FKBP. In this model, the drug-FKBP complex might sequester the small molecule and limit its degradation by restricting access to metabolic enzymes. In support of this idea, we found that adding FKBP blocks binding of FK506 to the common cytochrome P(450) enzyme CYP3A4 in vitro. To further test this idea, we have systematically modified a small collection of otherwise unrelated compounds, such that they acquire affinity for FKBP. Strikingly, we found that many of these synthetic derivatives, but not the unmodified parent compounds, are also protected from CYP3A4-mediated metabolism. Depending on the properties of the linker, the bifunctional molecules exhibited up to a 3.5-fold weaker binding to CYP3A4, and this protective effect was observed in the presence of either purified FKBP or FKBP-expressing cells. Together, these results suggest that the surprising pharmacology of rapamycin and FK506 might arise, in part, from binding to their abundant, intracellular target, FKBP. Furthermore, these findings provide a framework by which other small molecules might be systematically modified to impart this protective effect.

Expression and characterization of recombinant γ-tryptase

Tryptases are trypsin-like serine proteases whose expression is restricted to cells of hematopoietic origin, notably mast cells. gamma-Tryptase, a recently described member of the family also known as transmembrane tryptase (TMT), is a membrane-bound serine protease found in the secretory granules or on the surface of degranulated mast cells. The 321 amino acid protein contains an 18 amino acid propeptide linked to the catalytic domain (cd), followed by a single-span transmembrane domain. gamma-Tryptase is distinguished from other human mast cell tryptases by the presence of two unique cysteine residues, Cys(26) and Cys(145), that are predicted to form an intra-molecular disulfide bond linking the propeptide to the catalytic domain to form the mature, membrane-anchored two-chain enzyme. We expressed gamma-tryptase as either a soluble, single-chain enzyme with a C-terminal His tag (cd gamma-tryptase) or as a soluble pseudozymogen activated by enterokinase cleavage to form a two-chain protein with an N-terminal His tag (tc gamma-tryptase). Both recombinant proteins were expressed at high levels in Pichia pastoris and purified by affinity chromatography. The two forms of gamma-tryptase exhibit comparable kinetic parameters, indicating the propeptide does not contribute significantly to the substrate affinity or activity of the protease. Substrate and inhibitor library screening indicate that gamma-tryptase possesses a substrate preference and inhibitor profile distinct from that of beta-tryptase. Although the role of gamma-tryptase in mast cell function is unknown, our results suggest that it is likely to be distinct from that of beta-tryptase.

Mixed-type noncompetitive inhibition of anthrax lethal factor protease by aminoglycosides

We report a detailed kinetic investigation of the aminoglycosides neomycin B and neamine as inhibitors of the lethal factor protease from Bacillus anthracis. Both inhibitors display a mixed-type, noncompetitive kinetic pattern, which suggests the existence of multiple enzyme-inhibitor binding sites or the involvement of multiple structural binding modes at the same site. Quantitative analysis of the ionic strength effects by using the Debye-Hückel model revealed that the average interionic distance at the point of enzyme-inhibitor attachment is likely to be extremely short, which suggests specific, rather than nonspecific, binding. Only one ion pair seems to be involved in the binding process, which suggests the presence of a single binding site. Combining the results of our substrate competition studies with the ionic strength effects on the apparent inhibition constant, we propose that aminoglycoside inhibitors, such as neomycin B, bind to the lethal factor protease from B. anthracis in two different structural orientations. These results have important implications for the rational design of lethal factor protease inhibitors as possible therapeutic agents against anthrax. The strategies and methods we describe are general and can be employed to investigate in depth the mechanism of inhibition by other bioactive compounds.

Novel 5-azaindole factor VIIa inhibitors

The discovery and development of 5-azaindole factor VIIa inhibitors will be described.

Optimization of Subsite Binding to the β5 Subunit of the Human 20S Proteasome Using Vinyl Sulfones and 2-Keto-1,3,4-oxadiazoles: Syntheses and Cellular Properties of Potent, Selective Proteasome Inhibitors

Beginning with the peptide sequence Cbz-Ile-Glu(OtBu)-Ala-Leu found in PSI (3), a series of vinyl sulfones (VS) were synthesized for evaluation as inhibitors of the chymotrypsin-like activity of the 20S proteasome. Variations at the key P3 position confirmed the importance of a long side chain capped with a hydrophobic group for optimal potency, consistent with a model of binding to the S3 subsite. The tert-butyl glutamic ester initially used at P3 gave plasma unstable, insoluble compounds and was replaced with the better isostere, N-beta-neopentyl asparagine. The inhibitors were shortened by replacing the N-terminal Cbz-isoleucine with a p-tosyl group without loss of potency. Small l-amino acids were used at P2, where d-substitution was not tolerated. The resulting optimized P4-P3-P2 sequence was grafted onto a novel proteasome inhibitor warhead, 2-keto-1,3,4-oxadiazoles (KOD), to produce reversible, subnanomolar proteasome inhibitors that were 1000-fold selective versus cathepsin B (CatB), cathepsin S (CatS), and trypsin-like as well as PGPH-like proteasome activity. A number of compounds in both the VS and the KOD series exhibited growth inhibitory effects against the human prostate cancer cell line PC3 at submicromolar concentrations.

Design and synthesis of tri-ring P3 benzamide-containing aminonitriles as potent, selective, orally effective inhibitors of cathepsin K

We have prepared a series of achiral aminoacetonitriles, bearing tri-ring benzamide moieties and an aminocyclohexanecarboxylate residue at P2. This combination of binding elements resulted in sub-250 pM, reversible, selective, and orally bioavailable cathepsin K inhibitors. Lead compounds displayed single digit nanomolar inhibition in vitro (of rabbit osteoclast-mediated degradation of bovine bone). The best compound in this series, 39n (CRA-013783/L-006235), was orally bioavailable in rats, with a terminal half-life of over 3 h. 39n was dosed orally in ovariectomized rhesus monkeys once per day for 7 days. Collagen breakdown products were reduced by up to 76% dose-dependently. Plasma concentrations of 39n above the bone resorption IC50 after 24 h indicated a correlation between functional cellular and in vivo assays. Inhibition of collagen breakdown by cathepsin K inhibitors suggests this mechanism of action may be useful in osteoporosis and other indications involving bone resorption.

Mutations conferring resistance to SCH6, a novel hepatitis C virus NS3/4A protease inhibitor. Reduced RNA replication fitness and partial rescue by second-site mutations

Drug resistance is a major issue in the development and use of specific antiviral therapies. Here we report the isolation and characterization of hepatitis C virus RNA replicons resistant to a novel ketoamide inhibitor of the NS3/4A protease, SCH6 (originally SCH446211). Resistant replicon RNAs were generated by G418 selection in the presence of SCH6 in a dose-dependent fashion, with the emergence of resistance reduced at higher SCH6 concentrations. Sequencing demonstrated remarkable consistency in the mutations conferring SCH6 resistance in genotype 1b replicons derived from two different strains of hepatitis C virus, A156T/A156V and R109K. R109K, a novel mutation not reported previously to cause resistance to NS3/4A inhibitors, conferred moderate resistance only to SCH6. Structural analysis indicated that this reflects unique interactions of SCH6 with P'-side residues in the protease active site. In contrast, A156T conferred high level resistance to SCH6 and a related ketoamide, SCH503034, as well as BILN 2061 and VX-950. Unlike R109K, which had minimal impact on NS3/4A enzymatic function, A156T significantly reduced NS3/4A catalytic efficiency, polyprotein processing, and replicon fitness. However, three separate second-site mutations, P89L, Q86R, and G162R, were capable of partially reversing A156T-associated defects in polyprotein processing and/or replicon fitness, without significantly reducing resistance to the protease inhibitor.

Expression, crystallization, and three-dimensional structure of the catalytic domain of human plasma kallikrein

Plasma kallikrein is a serine protease that has many important functions, including modulation of blood pressure, complement activation, and mediation and maintenance of inflammatory responses. Although plasma kallikrein has been purified for 40 years, its structure has not been elucidated. In this report, we described two systems (Pichia pastoris and baculovirus/Sf9 cells) for expression of the protease domain of plasma kallikrein, along with the purification and high resolution crystal structures of the two recombinant forms. In the Pichia pastoris system, the protease domain was expressed as a heterogeneously glycosylated zymogen that was activated by limited trypsin digestion and treated with endoglycosidase H deglycosidase to reduce heterogeneity from the glycosylation. The resulting protein was chromatographically resolved into four components, one of which was crystallized. In the baculovirus/Sf9 system, homogeneous, crystallizable, and nonglycosylated protein was expressed after mutagenizing three asparagines (the glycosylation sites) to glutamates. When assayed against the peptide substrates, pefachrome-PK and oxidized insulin B chain, both forms of the protease domain were found to have catalytic activity similar to that of the full-length protein. Crystallization and x-ray crystal structure determination of both forms have yielded the first three-dimensional views of the catalytic domain of plasma kallikrein. The structures, determined at 1.85 A for the endoglycosidase H-deglycosylated protease domain produced from P. pastoris and at 1.40 A for the mutagenically deglycosylated form produced from Sf9 cells, show that the protease domain adopts a typical chymotrypsin-like serine protease conformation. The structural information provides insights into the biochemical and enzymatic properties of plasma kallikrein and paves the way for structure-based design of protease inhibitors that are selective either for or against plasma kallikrein.

Functional profiling of recombinant NS3 proteases from all four serotypes of dengue virus using tetrapeptide and octapeptide substrate libraries

Regulated proteolysis by the two-component NS2B/NS3 protease of dengue virus is essential for virus replication and the maturation of infectious virions. The functional similarity between the NS2B/NS3 proteases from the four genetically and antigenically distinct serotypes was addressed by characterizing the differences in their substrate specificity using tetrapeptide and octapeptide libraries in a positional scanning format, each containing 130,321 substrates. The proteases from different serotypes were shown to be functionally homologous based on the similarity of their substrate cleavage preferences. A strong preference for basic amino acid residues (Arg/Lys) at the P1 positions was observed, whereas the preferences for the P2-4 sites were in the order of Arg > Thr > Gln/Asn/Lys for P2, Lys > Arg > Asn for P3, and Nle > Leu > Lys > Xaa for P4. The prime site substrate specificity was for small and polar amino acids in P1' and P3'. In contrast, the P2' and P4' substrate positions showed minimal activity. The influence of the P2 and P3 amino acids on ground state binding and the P4 position for transition state stabilization was identified through single substrate kinetics with optimal and suboptimal substrate sequences. The specificities observed for dengue NS2B/NS3 have features in common with the physiological cleavage sites in the dengue polyprotein; however, all sites reveal previously unrecognized suboptimal sequences.

Enzyme family-specific and activity-based screening of chemical libraries using enzyme microarrays

The potential of protein microarrays in high-throughput screening (HTS) still remains largely unfulfilled, essentially because of the difficulty of extracting meaningful, quantitative data from such experiments. In the particular case of enzyme microarrays, low-molecular-weight fluorescent affinity labels (FALs) can function as ideally suited activity probes of the microarrayed enzymes. FALs form covalent bonds with enzymes in an activity-dependent manner and therefore can be used to characterize enzyme activity at each enzyme's address, as predetermined by the microarraying process. Relying on this principle, we introduce herein thematic enzyme microarrays (TEMA). In a kinetic setup we used TEMAs to determine the full set of kinetic constants and the reaction mechanism between the microarrayed enzymes (the theme of the microarray) and a family-wide FAL. Based on this kinetic understanding, in an HTS setup we established the practical and theoretical methodology for quantitative, multiplexed determination of the inhibition profile of compounds from a chemical library against each microarrayed enzyme. Finally, in a validation setup, K(i)(app) values and inhibitor profiles were confirmed and refined.

3,4-Disubstituted azetidinones as selective inhibitors of the cysteine protease cathepsin K. Exploring P3 elements for potency and selectivity

The synthesis of a series of highly potent and selective inhibitors of cathepsin K based on the 3,4-disubstituted azetidin-2-one warhead is reported. A high degree of potency and selectivity was achieved by introducing a basic nitrogen into the distal part of the P3 element of the molecule. Data from kinetic and mass spectrometry experiments are consistent with the interpretation that compounds of this series transiently acylate the sulfhydrile of cathepsin K.

Design and synthesis of beta-amino-alpha-hydroxy amide derivatives as inhibitors of MetAP2 and HUVEC growth

The rational design and synthesis of beta-amino-alpha-hydroxy amide derivatives as reversible inhibitors of methionine aminopeptidase-2 (MetAP2) with anti-proliferative activity against human umbilical vein endothelial cells (HUVECs) is described.

Biochemical characterization of prostasin, a channel activating protease

Human prostasin was recently identified as a potential regulator of epithelial sodium channel (ENaC) function. Through the use of positional scanning combinatorial substrate libraries, prostasin was shown to have a preference for poly-basic substrates: in position P4 preference was for arginine or lysine; in P3 preference was for histidine, lysine or arginine; in P2 preference was for basic or large hydrophobic amino acids; and in P1 preference was for arginine and lysine. P1', P2', and P3' displayed broad selectivity with the exception of a lack of activity for isoleucine, and P4' had a preference for small, unbranched, amino acids such as alanine and serine. A prostasin-preferred poly-basic cleavage site was found in the extracellular domains of the ENaC alpha- and beta-subunits, and may present a mechanism for prostasin activation. The absence of activity seen with substrates containing isoleucine in position P1' explains the inability of prostasin to autoactivate and suggests that prostasin proteolytic activity is regulated by an upstream protease. Prostasin activity was highly influenced by mono- and divalent metal ions which were potent inhibitors and substrate specific modulators of enzymatic activity. In the presence of sub-inhibitory concentrations of zinc, the activity of prostasin increased several-fold and its substrate specificity was significantly altered in favor of a strong preference for histidine in positions P3 or P4 of the substrate.

Dissecting and Designing Inhibitor Selectivity Determinants at the S1 Site Using an Artificial Ala190 Protease (Ala190 uPA)

A site-directed mutant of the serine protease urokinase-type plasminogen activator (uPA), was produced to assess the contribution of the Ser190 side-chain to the affinity and selectivity of lead uPA inhibitors in the absence of other differences present in comparisons of natural proteases. Crystallography and enzymology involving WT and Ala190 uPA were used to calculate free energy binding contributions of hydrogen bonds involving the Ser190 hydroxyl group (O(gamma)(Ser190)) responsible for the remarkable selectivity of 6-halo-5-amidinoindole and 6-halo-5-amidinobenzimidazole inhibitors toward uPA and against natural Ala190 protease anti-targets. Crystal structures of uPA complexes of novel, active site-directed arylguanidine and 2-aminobenzimidazole inhibitors of WT uPA, together with associated K(i) values for WT and Ala190 uPA, also indicate a significant role of Ser190 in the binding of these classes of uPA inhibitors. Structures and associated K(i) values for a lead inhibitor (CA-11) bound to uPA and to five other proteases, as well as for other leads bound to multiple proteases, help reveal the features responsible for the potency (K(i)=11nM) and selectivity of the remarkably small inhibitor, CA-11. The 6-fluoro-5-amidinobenzimidzole, CA-11, is more than 1000-fold selective against natural Ala190 protease anti-targets, and more than 100-fold selective against other Ser190 anti-targets.

Potent bivalent inhibition of human tryptase-beta by a synthetic inhibitor

Human tryptase-beta (HTbeta) is a unique serine protease exhibiting a frame-like tetramer structure with four active sites directed toward a central pore. Potent inhibition of HTbeta has been attained using CRA-2059. This compound has two phenylguanidinium head groups connected via a linker capable of spanning between two active sites. The properties of the CRA-2059:HTbeta interaction were defined in this study. Tight-binding reversible inhibition was observed with an inhibition constant (Ki) of 620 pM, an association rate constant of 7x10(7) M(-1) s(-1) and a relatively slow dissociation rate constant of 0.04 s(-1). Bivalent inhibition was demonstrated by displacement of p-aminobenzamidine from the primary specificity pocket with a stoichiometry, [CRA-2059]0/[HTbeta]0, of 0.5. The potency of the bivalent interaction was illustrated by CRA-2059 inhibition of HTbeta, 24% or 53% inhibited by pre-incubation with an irreversible inhibitor. Two interactions were observed consistent with mono- and bi-valent binding; the Ki value for bivalent inhibition was at least 10(4)-fold lower than that for monovalent inhibition. Comparison of the affinities of CRA-2059 and phenylguanidine for HTbeta finds an approximate doubling of the free energy change upon bivalent binding. This doubling suggests that the linker portion minimally hinders the binding of CRA-2059 to HTbeta. The potency of CRA-2059 is thus attributable to effective bivalent binding.

A new member of the bacterial ribonuclease inhibitor family from Saccharopolyspora erythraea

We have identified Sti, the gene of a ribonuclease inhibitor from Saccharopolyspora erythraea, by using a T7 phage display system. A specific phage has been isolated from a genome library by a biopanning procedure, using RNase Sa3, a ribonuclease from Streptomyces aureofaciens, as bait. Sti, a protein of 121 amino acid residues, with molecular mass 13059 Da, is a homolog of barstar and other microbial ribonuclease inhibitors. To overexpress its gene in Escherichia coli, we optimized the secondary structure of its mRNA by introducing a series of silent mutations. Soluble protein was isolated and purified to homogeneity. Inhibition constants of complex of Sti and RNase Sa3 or barnase were determined at pH 7 as 5 x 10(-12) or 7 x 10(-7), respectively.

Determination of accurate KI values for tight-binding enzyme inhibitors: an in silico study of experimental error and assay design

Determination of accurate K(I) values for tight-binding enzyme inhibitors is important from both a basic biochemistry point of view (understanding the differences in affinity of related molecules) and a medicinal chemistry vantage (developing structure-activity relationships (SAR)). It is advantageous to directly fit the quadratic equation describing tight-binding behavior, known commonly as the Morrison equation, to obtain these K(I) values. The results of simulated experiments that examine the effect of assay design and experimental error on the ability to accurately determine K(I) values at several [E]0/K(I-app) ratios are described. Input ("true") values of the uninhibited velocity, inhibition constant, and total enzyme concentration were used to calculate the velocity at various inhibitor concentrations. Gaussian error was introduced into the velocities and the simulated reactions were fit to estimate upsilon0, K(I), and [E]0. Recommendations for optimizing the inhibitor dilutions within the context of a 96-well-plate format and simple serial dilution steps are made. These include using three points to determine the enzyme concentration ([I]=0, 0.5[E]0, and [E]0), using a narrow dilution series with only two or three points to determine the asymptote at high inhibitor concentration, and avoiding fixing [E]0 to a constant value in the fitting if at all possible. The risks and rewards of fixing [E]0 to a constant value, especially the effect on SAR, are also examined.

Peptide Ketobenzoxazole Inhibitors Bound to Cathepsin K

Potent inhibitors of human cysteine proteases of the papain family have been made and assayed versus a number of relevant family members. We describe the synthesis of peptide alpha-ketoheterocyclic inhibitors that occupy binding subsites S1'-S3 of the cysteine protease substrate recognition cleft and that form a reversible covalent bond with the Cys 25 nucleophile. X-ray crystal structures of cathepsin K both unbound and complexed with inhibitors provide detailed information on protease/inhibitor interactions and suggestions for the design of tight-binding, selective molecules.

Kinetic determination of tight-binding impurities in enzyme inhibitors

A novel rate equation to characterize the dose-response behavior of a moderately potent ("classical") enzyme inhibitor contaminated with a very potent ("tight-binding") impurity is derived. Mathematical properties of the new rate equation show that, for such contaminated materials, experimentally observed I(50) values are ambiguous. The four-parameter logistic equation, conventionally used to determine I(50) values, cannot be used to detect the presence of tight-binding impurities in inhibitor samples. In contrast, fitting the newly derived rate equation to inhibitor dose- response curves can, in favorable cases, reveal whether the unknown material is chemically homogeneous or whether it is contaminated with a tight-binding impurity. The limitations of our method with respect to the detectable range of inhibition constants (both classical and tight-binding) were examined by using Monte-Carlo simulations. To test the new analytical procedure experimentally, we added a small amount (0.02 mole%) of a tight-binding impurity (K(i)=0.065 nM) to an otherwise weak inhibitor of human mast-cell tryptase (K(i)=50.4 microM). The resulting material was treated as "unknown." Our kinetic equation predicts that such adulterated material should show I(50)=0.40 microM, which was identical to the experimentally observed value. The best-fit value of the apparent inhibition constants for the tight-binding inhibitor was K(i)=(0.107+/-0.035)nM, close to the true value of 0.065 nM.

Elaborate manifold of short hydrogen bond arrays mediating binding of active site-directed serine protease inhibitors

An extensive structural manifold of short hydrogen bond-mediated, active site-directed, serine protease inhibition motifs is revealed in a set of over 300 crystal structures involving a large suite of small molecule inhibitors (2-(2-phenol)-indoles and 2-(2-phenol)-benzimidazoles) determined over a wide range of pH (3.5-11.4). The active site hydrogen-bonding mode was found to vary markedly with pH, with the steric and electronic properties of the inhibitor, and with the type of protease (trypsin, thrombin or urokinase type plasminogen activator (uPA)). The pH dependence of the active site hydrogen-bonding motif is often intricate, constituting a distinct fingerprint of each complex. Isosteric replacements or minor substitutions within the inhibitor that modulate the pK(a) of the phenol hydroxyl involved in short hydrogen bonding, or that affect steric interactions distal to the active site, can significantly shift the pH-dependent structural profile characteristic of the parent scaffold, or produce active site-binding motifs unique to the bound analog. Ionization equilibria at the active site associated with inhibitor binding are probed in a series of the protease-inhibitor complexes through analysis of the pH dependence of the structure and environment of the active site-binding groups involved in short hydrogen bond arrays. Structures determined at high pH (>11), suggest that the pK(a) of His57 is dramatically elevated, to a value as high as approximately 11 in certain complexes. K(i) values involving uPA and trypsin determined as a function of pH for a set of inhibitors show pronounced parabolic pH dependence, the pH for optimal inhibition governed by the pK(a) of the inhibitor phenol involved in short hydrogen bonds. Comparison of structures of trypsin, thrombin and uPA, each bound by the same inhibitor, highlights important structural variations in the S1 and active sites accessible for engineering notable selectivity into remarkably small molecules with low nanomolar K(i) values.

Novel route to the synthesis of peptides containing 2-amino-1′-hydroxymethyl ketones and their application as cathepsin K inhibitors

Cathepsin K is highly expressed in human osteoclasts, and is implicated in bone resorption. This makes it an attractive target for the treatment of osteoporosis. Peptides containing 2-amino-1'-hydroxymethyl ketones and 2-amino-1'-alkoxymethyl ketones were discovered as potent inhibitors of cathepsin K. A novel synthetic route was devised to facilitate rapid elucidation of the SAR of these inhibitors. The synthesis and SAR of hydroxymethyl ketones are presented.

Contribution of Multicentered Short Hydrogen Bond Arrays to Potency of Active Site-Directed Serine Protease Inhibitors

We describe and compare the pH dependencies of the potencies and of the bound structures of two inhibitor isosteres that form multicentered short hydrogen bond arrays at the active sites of trypsin, thrombin, and urokinase type plasminogen activator (urokinase or uPA) over certain ranges of pH. Depending on the pH, short hydrogen bond arrays at the active site are mediated by two waters, one in the oxyanion hole (H(2)O(oxy)) and one on the other (S2) side of the inhibitor (H(2)O(S2)), by one water (H(2)O(oxy)), or by no water. The dramatic variation in the length of the active site hydrogen bonds as a function of pH, of inhibitor, and of enzyme, along with the involvement or absence of ordered water, produces a large structural manifold of active site hydrogen bond motifs. Diverse examples of multicentered and two-centered short hydrogen bond arrays, both at and away from the active site, recently discovered in several protein crystal systems, suggest that short hydrogen bonds in proteins may be more common than has been recognized. The short hydrogen bond arrays resemble one another with respect to ionic nature, highly polar environment, multitude of associated ordinary hydrogen bonds, and disparate pK(a) values of participating groups. Comparison of structures and K(i) values of trypsin complexes at pH values where the multicentered short hydrogen bond arrays mediating inhibitor binding are present or absent indicate that these arrays have a minor effect on inhibitor potency. These features suggest little covalent nature within the short hydrogen bonds, despite their extraordinary shortness (as short as 2.0 A).

Exploiting subsite S1 of trypsin-like serine proteases for selectivity: potent and selective inhibitors of urokinase-type plasminogen activator

A nonselective inhibitor of trypsin-like serine proteases, 2-(2-hydroxybiphenyl-3-yl)-1H-indole-5-carboxamidine (1) (Verner, E.; Katz, B. A.; Spencer, J.; Allen, D.; Hataye, J.; Hruzewicz, W.; Hui, H. C.; Kolesnikov, A.; Li, Y.; Luong, C.; Martelli, A.; Radika. K.; Rai, R.; She, M.; Shrader, W.; Sprengeler, P. A.; Trapp, S.; Wang, J.; Young, W. B.; Mackman, R. L. J. Med. Chem. 2001, 44, 2753-2771) has been optimized through minor structural changes on the S1 binding group to afford remarkably selective and potent inhibitors of urokinase-type plasminogen activator (uPA). The trypsin-like serine proteases(1) that comprise drug targets can be broadly categorized into two subfamilies, those with Ser190 and those with Ala190. A single-atom modification, for example, replacement of hydrogen for chlorine at the 6-position of the 5-amidinoindole P1 group on 1, generated up to 6700-fold selectivity toward the Ser190 enzymes and against the Ala190 enzymes. The larger chlorine atom displaces a water molecule (H(2)O1(S1)) that binds near residue 190 in all the complexes of 1, and related inhibitors, in uPA, thrombin, and trypsin. The water molecule, H(2)O1(S1), in both the Ser190 or Ala190 enzymes, hydrogen bonds with the amidine N1 nitrogen of the inhibitor. When it is displaced, a reduction in affinity toward the Ala190 enzymes is observed due to the amidine N1 nitrogen of the bound inhibitor being deprived of a key hydrogen-bonding partner. In the Ser190 enzymes the affinity is maintained since the serine hydroxyl oxygen O gamma(Ser190) compensates for the displaced water molecule. High-resolution crystallography provided evidence for the displacement of the water molecule and validated the design rationale. In summation, a novel and powerful method for engineering selectivity toward Ser190 proteases and against Ala190 proteases without substantially increasing molecular weight is described.

Engineering inhibitors highly selective for the S1 sites of Ser190 trypsin-like serine protease drug targets

BACKGROUND

Involved or implicated in a wide spectrum of diseases, trypsin-like serine proteases comprise well studied drug targets and anti-targets that can be subdivided into two major classes. In one class there is a serine at position 190 at the S1 site, as in urokinase type plasminogen activator (urokinase or uPA) and factor VIIa, and in the other there is an alanine at 190, as in tissue type plasminogen activator (tPA) and factor Xa. A hydrogen bond unique to Ser190 protease-arylamidine complexes between O gamma(Ser190) and the inhibitor amidine confers an intrinsic preference for such inhibitors toward Ser190 proteases over Ala190 counterparts.

RESULTS

Based on the structural differences between the S1 sites of Ser190 and Ala190 protease-arylamidine complexes, we amplified the selectivity of amidine inhibitors toward uPA and against tPA, by factors as high as 220-fold, by incorporating a halo group ortho to the amidine of a lead inhibitor scaffold. Comparison of K(i) values of such halo-substituted and parent inhibitors toward a panel of Ser190 and Ala190 proteases demonstrates pronounced selectivity of the halo analogs for Ser190 proteases over Ala190 counterparts. Crystal structures of Ser190 proteases, uPA and trypsin, and of an Ala190 counterpart, thrombin, bound by a set of ortho (halo, amidino) aryl inhibitors and of non-halo parents reveal the structural basis of the exquisite selectivity and validate the design principle.

CONCLUSIONS

Remarkable selectivity enhancements of exceptionally small inhibitors are achieved toward the uPA target over the highly similar tPA anti-target through a single atom substitution on an otherwise relatively non-selective scaffold. Overall selectivities for uPA over tPA as high as 980-fold at physiological pH were realized. The increase in selectivity results from the displacement of a single bound water molecule common to the S1 site of both the uPA target and the tPA anti-target because of the ensuing deficit in hydrogen bonding of the arylamidine inhibitor when bound in the Ala190 protease anti-target.

Development of serine protease inhibitors displaying a multicentered short (<2.3 A) hydrogen bond binding mode: inhibitors of urokinase-type plasminogen activator and factor Xa

Novel scaffolds that bind to serine proteases through a unique network of short hydrogen bonds to the catalytic Ser195 have been developed. The resulting potent serine protease inhibitors were designed from lead molecule 2-(2-hydroxyphenyl)1H-benzoimidazole-5-carboxamidine, 6b, which is known to display several modes of binding. For instance, 6b can recruit zinc and bind in a manner similar to that reported by bis(5-amidino-2-benzimidazolyl)methane (BABIM) (Nature 1998, 391, 608-612).(1) Alternatively, 6b can bind in the absence of zinc through a multicentered network of short (<2.3 A) hydrogen bonds. The lead structure was optimized in the zinc-independent binding mode toward a panel of six human serine proteases to yield optimized inhibitors such as 2-(3-bromo-2-hydroxy-5-methylphenyl)-1H-indole-5-carboxamidine, 22a, and 2-(2-hydroxybiphenyl-3-yl)-1H-indole-5-carboxamidine, 22f. Structure-activity relationships determined that, apart from the amidine function, an indole or benzimidazole and an ortho substituted phenol group were also essential components for optimal potency. The affinities (K(i)) of 22a and 22f, for example, bearing these groups ranged from 8 to 600 nM toward a panel of six human serine proteases. High-resolution crystal structures revealed that the binding mode of these molecules in several of the enzymes was identical to that of 6b and involved short (<2.3 A) hydrogen bonds among the inhibitor hydroxyl oxygen, Ser195, and a water molecule trapped in the oxyanion hole. In summation, novel and potent trypsin-like serine protease inhibitors possessing a unique mode of binding have been discovered.

Neutral inhibitors of the serine protease factor Xa

A neutral inhibitor of the serine protease factor Xa was identified via a high-throughput screen of a commercial library. The initial lead 1 demonstrated reversible and competitive inhibition kinetics for factor Xa and possessed a high degree of selectivity versus other related serine proteases. Initial modeling efforts and the generation of a series of analogues of 1 are described.

Mechanistic implications of methylglyoxal synthase complexed with phosphoglycolohydroxamic acid as observed by X-ray crystallography and NMR spectroscopy

Methylglyoxal synthase (MGS) and triosephosphate isomerase (TIM) share neither sequence nor structural similarities, yet the reactions catalyzed by both enzymes are similar, in that both initially convert dihydroxyacetone phosphate to a cis-enediolic intermediate. This enediolic intermediate is formed from the abstraction of the pro-S C3 proton of DHAP by Asp-71 of MGS or the pro-R C3 proton of DHAP by Glu-165 of TIM. MGS then catalyzes the elimination of phosphate from this enediolic intermediate to form the enol of methylglyoxal, while TIM catalyzes proton donation to C2 to form D-glyceraldehyde phosphate. A competitive inhibitor of TIM, phosphoglycolohydroxamic acid (PGH) is found to be a tight binding competitive inhibitor of MGS with a K(i) of 39 nM. PGH's high affinity for MGS may be due in part to a short, strong hydrogen bond (SSHB) from the NOH of PGH to the carboxylate of Asp-71. Evidence for this SSHB is found in X-ray, 1H NMR, and fractionation factor data. The X-ray structure of the MGS homohexamer complexed with PGH at 2.0 A resolution shows this distance to be 2.30-2.37 +/- 0.24 A. 1H NMR shows a PGH-dependent 18.1 ppm signal that is consistent with a hydrogen bond length of 2.49 +/- 0.02 A. The D/H fractionation factor (phi = 0.43 +/- 0.02) is consistent with a hydrogen bond length of 2.53 +/- 0.01 A. Further, 15N NMR suggests a significant partial positive charge on the nitrogen atom of bound PGH, which could strengthen hydrogen bond donation to Asp-71. Both His-98 and His-19 are uncharged in the MGS-PGH complex on the basis of the chemical shifts of their Cdelta and C(epsilon) protons. The crystal structure reveals that Asp-71, on the re face of PGH, and His-19, on the si face of PGH, both approach the NO group of the analogue, while His-98, in the plane of PGH, approaches the carbonyl oxygen of the analogue. The phosphate group of PGH accepts nine hydrogen bonds from seven residues and is tilted out of the imidate plane of PGH toward the re face. Asp-71 and phosphate are thus positioned to function as the base and leaving group, respectively, in a concerted suprafacial 1,4-elimination of phosphate from the enediolic intermediate in the second step of the MGS reaction. Combined, these data suggest that Asp-71 is the one base that initially abstracts the C3 pro-S proton from DHAP and subsequently the 3-OH proton from the enediolic intermediate. This mechanism is compared to an alternative TIM-like mechanism for MGS, and the relative merits of both mechanisms are discussed.

A novel serine protease inhibition motif involving a multi-centered short hydrogen bonding network at the active site

We describe a new serine protease inhibition motif in which binding is mediated by a cluster of very short hydrogen bonds (<2.3 A) at the active site. This protease-inhibitor binding paradigm is observed at high resolution in a large set of crystal structures of trypsin, thrombin, and urokinase-type plasminogen activator (uPA) bound with a series of small molecule inhibitors (2-(2-phenol)indoles and 2-(2-phenol)benzimidazoles). In each complex there are eight enzyme-inhibitor or enzyme-water-inhibitor hydrogen bonds at the active site, three of which are very short. These short hydrogen bonds connect a triangle of oxygen atoms comprising O(gamma)(Ser195), a water molecule co-bound in the oxyanion hole (H(2)O(oxy)), and the phenolate oxygen atom of the inhibitor (O6'). Two of the other hydrogen bonds between the inhibitor and active site of the trypsin and uPA complexes become short in the thrombin counterparts, extending the three-centered short hydrogen-bonding array into a tetrahedral array of atoms (three oxygen and one nitrogen) involved in short hydrogen bonds. In the uPA complexes, the extensive hydrogen-bonding interactions at the active site prevent the inhibitor S1 amidine from forming direct hydrogen bonds with Asp189 because the S1 site is deeper in uPA than in trypsin or thrombin. Ionization equilibria at the active site associated with inhibitor binding are probed through determination and comparison of structures over a wide range of pH (3.5 to 11.4) of thrombin complexes and of trypsin complexes in three different crystal forms. The high-pH trypsin-inhibitor structures suggest that His57 is protonated at pH values as high as 9.5. The pH-dependent inhibition of trypsin, thrombin, uPA and factor Xa by 2-(2-phenol)benzimidazole analogs in which the pK(a) of the phenol group is modulated is shown to be consistent with a binding process involving ionization of both the inhibitor and the enzyme. These data further suggest that the pK(a) of His57 of each protease in the unbound state in solution is about the same, approximately 6.8. By comparing inhibition constants (K(i) values), inhibitor solubilities, inhibitor conformational energies and corresponding structures of short and normal hydrogen bond-mediated complexes, we have estimated the contribution of the short hydrogen bond networks to inhibitor affinity ( approximately 1.7 kcal/mol). The structures and K(i) values associated with the short hydrogen-bonding motif are compared with those corresponding to an alternate, Zn(2+)-mediated inhibition motif at the active site. Structural differences among apo-enzymes, enzyme-inhibitor and enzyme-inhibitor-Zn(2+) complexes are discussed in the context of affinity determinants, selectivity development, and structure-based inhibitor design.

A broad spectrum Kunitz type serine protease inhibitor secreted by the hookworm Ancylostoma ceylanicum

Although blood-feeding hookworms infect over a billion people worldwide, little is known about the molecular mechanisms through which these parasitic nematodes cause gastrointestinal hemorrhage and iron deficiency anemia. A cDNA corresponding to a secreted Kunitz type serine protease inhibitor has been cloned from adult Ancylostoma ceylanicum hookworm RNA. The translated sequence of the A. ceylanicum Kunitz type inhibitor 1 (AceKI-1) cDNA predicts a 16-amino acid secretory signal sequence, followed by a 68-amino acid mature protein with a molecular mass of 7889 daltons. Recombinant protein (rAceKI-1) was purified from induced lysates of Escherichia coli transformed with the rAceKI-1/pET 28a plasmid, and in vitro studies demonstrate that rAceKI-1 is a tight binding inhibitor of the serine proteases chymotrypsin, pancreatic elastase, neutrophil elastase, and trypsin. AceKI-1 inhibitory activity is present in soluble protein extracts and excretory/secretory products of adult hookworms but not the infective third stage larvae. The native AceKI-1 inhibitor has been purified to homogeneity from soluble extracts of adult A. ceylanicum using size exclusion and reverse-phase high pressure liquid chromatography. As a potent inhibitor of mammalian intestinal proteases, AceKI-1 may play a role in parasite survival and the pathogenesis of hookworm anemia.