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

Urokinase plasminogen activator as an anti-metastasis target: inhibitor design principles, recent amiloride derivatives, and issues with human/mouse species selectivity

The urokinase plasminogen activator (uPA) is a widely studied anticancer drug target with multiple classes of inhibitors reported to date. Many of these inhibitors contain amidine or guanidine groups, while others lacking these groups show improved oral bioavailability. Most of the X-ray co-crystal structures of small molecule uPA inhibitors show a key salt bridge with the side chain carboxylate of Asp189 in the S1 pocket of uPA. This review summarises the different classes of uPA inhibitors, their binding interactions and experimentally measured inhibitory potencies and highlights species selectivity issues with attention to recently described 6-substituted amiloride and 5‑N,N-(hexamethylene)amiloride (HMA) derivatives.

Disruption of Water Networks is the Cause of Human/Mouse Species Selectivity in Urokinase Plasminogen Activator (uPA) Inhibitors Derived from Hexamethylene Amiloride (HMA)

The urokinase plasminogen activator (uPA) plays a critical role in tumor cell invasion and migration and is a promising antimetastasis target. 6-Substituted analogues of 5-N,N-(hexamethylene)amiloride (HMA) are potent and selective uPA inhibitors that lack the diuretic and antikaliuretic properties of the parent drug amiloride. However, the compounds display pronounced selectivity for human over mouse uPA, thus confounding interpretation of data from human xenograft mouse models of cancer. Here, computational and experimental findings reveal that residue 99 is a key contributor to the observed species selectivity, whereby enthalpically unfavorable expulsion of a water molecule by the 5-N,N-hexamethylene ring occurs when residue 99 is Tyr (as in mouse uPA). Analogue 7 lacking the 5-N,N-hexamethylene ring maintained similar water networks when bound to human and mouse uPA and displayed reduced selectivity, thus supporting this conclusion. The study will guide further optimization of dual-potent human/mouse uPA inhibitors from the amiloride class as antimetastasis drugs.

How a Fragment Draws Attention to Selectivity Discriminating Features between the Related Proteases Trypsin and Thrombin

In the S₁ pocket, the serine proteases thrombin and trypsin commonly feature Asp189 and a Ala190Ser and Glu192Gln exchange. Nevertheless, thrombin cleaves peptide chains solely after Arg, and trypsin after Lys and Arg. Thrombin exhibits a Na+-binding site next to Asp189, which is missing in trypsin. The fragment benzylamine shows direct H-bonding to Asp189 in trypsin, while in thrombin, it forms an H-bond to Glu192. A series of fragments and expanded ligands were studied against both enzymes and mutated variants by crystallography and ITC. The selectivity-determining features of both S₁ pockets are difficult to assign to one dominating factor. The Ala190Ser and Glu192Gln replacements may be regarded as highly conserved as no structural and affinity changes are observed between both proteases. With respect to charge distribution, Glu192, together with the thrombin-specific sodium ion, helps in creating an electrostatic gradient across the S₁ pocket. This feature is definitely absent in trypsin but important for selectivity along with solvation-pattern differences in the S₁ pocket.

Potent DNA gyrase inhibitors bind asymmetrically to their target using symmetrical bifurcated halogen bonds

Novel bacterial type II topoisomerase inhibitors (NBTIs) stabilize single-strand DNA cleavage breaks by DNA gyrase but their exact mechanism of action has remained hypothetical until now. We have designed a small library of NBTIs with an improved DNA gyrase-binding moiety resulting in low nanomolar inhibition and very potent antibacterial activity. They stabilize single-stranded cleavage complexes and, importantly, we have obtained the crystal structure where an NBTI binds gyrase–DNA in a single conformation lacking apparent static disorder. This directly proves the previously postulated NBTI mechanism of action and shows that they stabilize single-strand cleavage through asymmetric intercalation with a shift of the scissile phosphate. This crystal stucture shows that the chlorine forms a halogen bond with the backbone carbonyls of the two symmetry-related Ala68 residues. To the best of our knowledge, such a so-called symmetrical bifurcated halogen bond has not been identified in a biological system until now.

The mechanism of DNA gyrase inhibitor stabilization of single-strand DNA cleavage breaks by DNA gyrase has been hypothetical. Here, the authors show experimental evidence of the mechanism using a library of inhibitors with improved binding and employ crystal analysis to show bifurcated halogen bonding.

Structural analysis of experimental drugs binding to the SARS-CoV-2 target TMPRSS2

The emergence of SARS-CoV-2 has prompted a worldwide health emergency. There is an urgent need for therapeutics, both through the repurposing of approved drugs and the development of new treatments. In addition to the viral drug targets, a number of human drug targets have been suggested. In theory, targeting human proteins should provide an advantage over targeting viral proteins in terms of drug resistance, which is commonly a problem in treating RNA viruses. This paper focuses on the human protein TMPRSS2, which supports coronavirus life cycles by cleaving viral spike proteins. The three-dimensional structure of TMPRSS2 is not known and so we have generated models of the TMPRSS2 in the apo state as well as in complex with a peptide substrate and putative inhibitors to aid future work. Importantly, many related human proteases have 80% or higher identity with TMPRSS2 in the S1-S1' subsites, with plasminogen and urokinase-type plasminogen activator (uPA) having 95% identity. We highlight 376 approved, investigational or experimental drugs targeting S1A serine proteases that may also inhibit TMPRSS2. Whilst the presence of a relatively uncommon lysine residue in the S2/S3 subsites means that some serine protease inhibitors will not inhibit TMPRSS2, this residue is likely to provide a handle for selective targeting in a focused drug discovery project. We discuss how experimental drugs targeting related serine proteases might be repurposed as TMPRSS2 inhibitors to treat coronaviruses.

Enhancing Water Sampling in Free Energy Calculations with Grand Canonical Monte Carlo

The prediction of protein-ligand binding affinities using free energy perturbation (FEP) is becoming increasingly routine in structure-based drug discovery. Most FEP packages use molecular dynamics (MD) to sample the configurations of proteins and ligands, as MD is well-suited to capturing coupled motion. However, MD can be prohibitively inefficient at sampling water molecules that are buried within binding sites, which has severely limited the domain of applicability of FEP and its prospective usage in drug discovery. In this paper, we present an advancement of FEP that augments MD with grand canonical Monte Carlo (GCMC), an enhanced sampling method, to overcome the problem of sampling water. We accomplished this without degrading computational performance. On both old and newly assembled data sets of protein-ligand complexes, we show that the use of GCMC in FEP is essential for accurate and robust predictions for ligand perturbations that disrupt buried water.

Synthesis and Structure-Activity Relationships of N-(4-Benzamidino)-Oxazolidinones: Potent and Selective Inhibitors of Kallikrein-Related Peptidase 6

Kallikrein-related peptidase 6 (KLK6) is a secreted serine protease that belongs to the family of tissue kallikreins. Aberrant expression of KLK6 has been found in different cancers and neurodegenerative diseases, and KLK6 is currently studied as a potential target in these pathologies. We report a novel series of KLK6 inhibitors discovered in a high-throughput screen within the European Lead Factory program. Structure-guided design based on docking studies enabled rapid progression of a hit cluster to inhibitors with improved potency, selectivity and pharmacokinetic properties. In particular, inhibitors 32 ((5R)-3-(4-carbamimidoylphenyl)-N-((S)-1-(naphthalen-1-yl)propyl)-2-oxooxazolidine-5-carboxamide) and 34 ((5R)-3-(6-carbamimidoylpyridin-3-yl)-N-((1S)-1-(naphthalen-1-yl)propyl)-2-oxooxazolidine-5-carboxamide) have single-digit nanomolar potency against KLK6, with over 25-fold and 100-fold selectivities against the closely related enzyme trypsin, respectively. The most potent compound, 32, effectively reduces KLK6-dependent invasion of HCT116 cells. The high potency in combination with good solubility and low clearance of 32 make it a good chemical probe for KLK6 target validation in vitro and potentially in vivo.

Proteolysis is the most fundamental property of malignancy and its inhibition may be used therapeutically (Review)

The mortality rates of cancer patients decreased by ~1.5% per year between 2001 and 2015, although the decrease depends on patient sex, ethnic group and type of malignancy. Cancer remains a significant global health problem, requiring a search for novel treatments. The most common property of malignant tumors is their capacity to invade adjacent tissue and to metastasize, and this cancer aggressiveness is contingent on overexpression of proteolytic enzymes. The components of the plasminogen activation system (PAS) and the metal-loproteinase family [mainly matrix metalloproteinases (MMPs)] are overexpressed in malignant tumors, driving the local invasion, metastasis and angiogenesis. This is the case for numerous types of cancer, such as breast, colon, prostate and oral carcinoma, among others. Present chemotherapeutics agents typically attack all dividing cells; however, for future therapeutic agents to be clinically successful, they need to be highly selective for a specific protein(s) and act on the cancerous tissues without adverse systemic effects. Inhibition of proteolysis in cancerous tissue has the ability to attenuate tumor invasion, angiogenesis and migration. For that purpose, inhibiting both PAS and MMPs may be another approach, since the two groups of enzymes are overexpressed in cancer. In the present review, the roles and new findings on PAS and MMP families in cancer formation, growth and possible treatments are discussed.

Can components of the plasminogen activation system predict the outcome of kidney transplants?

Proteolytic and antiproteolytic enzymes play a critical role in the physiology and pathology of different stages of human life. One of the important members of the proteolytic family is the plasminogen activation system (PAS), which includes several elements crucial for this review: the 50 kDa glycoprotein plasminogen activator inhibitor 1 (PAI-1) that inhibits tissue-type (tPA) and urokinase-type plasminogen activator (uPA). These two convert plasminogen into its active form named plasmin that can lyse a broad spectrum of proteins. Urokinase receptor (uPAR) is the binding site of uPA. This glycoprotein on the cell surface facilitates urokinase activation of plasminogen, creating high proteolytic activity close to the cell surface. PAS activities have been reported to predict the outcome of kidney transplants. However, reports on expression of PAS in kidney transplants seem to be controversial. On the one hand there are reports that impaired proteolytic activity leads to induction of chronic allograft nephropathy, while on the other hand treatment with uPA and tPA can restore function of acute renal transplants. In this comprehensive review we describe the complexity of the PAS as well as biological effects of the PAS on renal allografts, and provide a possible explanation of the reported controversy.

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.

3D proteochemometrics: using three-dimensional information of proteins and ligands to address aspects of the selectivity of serine proteases

The high similarity between certain sub-pockets of serine proteases may lead to low selectivity of protease inhibitors. Therefore the application of proteochemometrics (PCM), which quantifies the relationship between protein/ligand descriptors and affinity for multiple ligands and targets simultaneously, is useful to understand and improve the selectivity profiles of potential inhibitors. In this study, protein field-based PCM that uses knowledge-based and WaterMap derived fields to describe proteins in combination with 2D (RDKit and MOE fingerprints) and 3D (4 point pharmacophoric fingerprints and GRIND) ligand descriptors was used to model the bioactivities of 24 homologous serine proteases and 5863 inhibitors in an integrated fashion. Of the multiple field-based PCM models generated based on different ligand descriptors, RDKit fingerprints showed the best performance in terms of external prediction with Rtest2 of 0.72 and RMSEP of 0.81. Further, visual interpretation of the models highlights sub-pocket specific regions that influence affinity and selectivity of serine protease inhibitors.

Relationships between hydrogen bonds and halogen bonds in biological systems

The recent recognition that halogen bonding (XB) plays important roles in the recognition and assembly of biological molecules has led to new approaches in medicinal chemistry and biomolecular engineering. When designing XBs into strategies for rational drug design or into a biomolecule to affect its structure and function, we must consider the relationship between this interaction and the more ubiquitous hydrogen bond (HB). In this review, we explore these relationships by asking whether and how XBs can replace, compete against or behave independently of HBs in various biological systems. The complex relationships between the two interactions inform us of the challenges we face in fully utilizing XBs to control the affinity and recognition of inhibitors against their therapeutic targets, and to control the structure and function of proteins, nucleic acids and other biomolecular scaffolds.

Lessons from Hot Spot Analysis for Fragment-Based Drug Discovery

Analysis of binding energy hot spots at protein surfaces can provide crucial insights into the prospects for successful application of fragment-based drug discovery (FBDD), and whether a fragment hit can be advanced into a high-affinity, drug-like ligand. The key factor is the strength of the top ranking hot spot, and how well a given fragment complements it. We show that published data are sufficient to provide a sophisticated and quantitative understanding of how hot spots derive from a protein 3D structure, and how their strength, number, and spatial arrangement govern the potential for a surface site to bind to fragment-sized and larger ligands. This improved understanding provides important guidance for the effective application of FBDD in drug discovery.

Discovery of novel P1 groups for coagulation factor VIIa inhibition using fragment-based screening

A multidisciplinary, fragment-based screening approach involving protein ensemble docking and biochemical and NMR assays is described. This approach led to the discovery of several structurally diverse, neutral surrogates for cationic factor VIIa P1 groups, which are generally associated with poor pharmacokinetic (PK) properties. Among the novel factor VIIa inhibitory fragments identified were aryl halides, lactams, and heterocycles. Crystallographic structures for several bound fragments were obtained, leading to the successful design of a potent factor VIIa inhibitor with a neutral lactam P1 and improved permeability.

Application of molecular modeling to urokinase inhibitors development

Urokinase-type plasminogen activator (uPA) plays an important role in the regulation of diverse physiologic and pathologic processes. Experimental research has shown that elevated uPA expression is associated with cancer progression, metastasis, and shortened survival in patients, whereas suppression of proteolytic activity of uPA leads to evident decrease of metastasis. Therefore, uPA has been considered as a promising molecular target for development of anticancer drugs. The present study sets out to develop the new selective uPA inhibitors using computer-aided structural based drug design methods. Investigation involves the following stages: computer modeling of the protein active site, development and validation of computer molecular modeling methods: docking (SOL program), postprocessing (DISCORE program), direct generalized docking (FLM program), and the application of the quantum chemical calculations (MOPAC package), search of uPA inhibitors among molecules from databases of ready-made compounds to find new uPA inhibitors, and design of new chemical structures and their optimization and experimental examination. On the basis of known uPA inhibitors and modeling results, 18 new compounds have been designed, calculated using programs mentioned above, synthesized, and tested in vitro. Eight of them display inhibitory activity and two of them display activity about 10 μM.

Identification of the first synthetic inhibitors of the type II transmembrane serine protease TMPRSS2 suitable for inhibition of influenza virus activation

TMPRSS2 (transmembrane serine proteinase 2) is a multidomain type II transmembrane serine protease that cleaves the surface glycoprotein HA (haemagglutinin) of influenza viruses with a monobasic cleavage site, which is a prerequisite for virus fusion and propagation. Furthermore, it activates the fusion protein F of the human metapneumovirus and the spike protein S of the SARS-CoV (severe acute respiratory syndrome coronavirus). Increased TMPRSS2 expression was also described in several tumour entities. Therefore TMPRSS2 emerged as a potential target for drug design. The catalytic domain of TMPRSS2 was expressed in Escherichia coli and used for an inhibitor screen with previously synthesized inhibitors of various trypsin-like serine proteases. Two inhibitor types were identified which inhibit TMPRSS2 in the nanomolar range. The first series comprises substrate analogue inhibitors containing a 4-amidinobenzylamide moiety at the P1 position, whereby some of these analogues possess inhibition constants of approximately 20 nM. An improved potency was found for a second type derived from sulfonylated 3-amindinophenylalanylamide derivatives. The most potent derivative of this series inhibits TMPRSS2 with a K(i) value of 0.9 nM and showed an efficient blockage of influenza virus propagation in human airway epithelial cells. On the basis of the inhibitor studies, a series of new fluorogenic substrates containing a D-arginine residue at the P3 position was synthesized, some of them were efficiently cleaved by TMPRSS2.

Structural model of haptoglobin and its complex with the anticoagulant ecotin variants: structure-activity relationship study and analysis of interactions

Recently the literature described the binding of Haptoglobin (HP) with ecotin, a fold-specific serine-proteases inhibitor with an anticoagulant profile and produced by Escherichia coli. In this work, we used some in silico and in vitro techniques to evaluate HP 3D-fold and its interaction with wild-type ecotin and two variants. Our data showed HP models conserved trypsin fold, in agreement to the in vitro immunological recognition of HP by trypsin antibodies. The analysis of the three ecotin-HP complexes using the mutants RR and TSRR/R besides the wild type revealed several hydrogen bonds between HP and ecotin secondary site. These data are in agreement with the in vitro PAGE assays that showed the HP-RR complex in native gel conditions. Interestingly, the ternary complex interactions varied depending on the inhibitor structure and site-directed mutation. The interaction of HP with TSRR/R involved new residues compared to wild type, which infers a binding energy increase caused by the mutation.

Natural and engineered plasmin inhibitors: applications and design strategies

The serine protease plasmin is ubiquitously expressed throughout the human body in the form of the zymogen plasminogen. Conversion to active plasmin occurs through enzymatic cleavage by plasminogen activators. The plasminogen activator/plasmin system has a well-established function in the removal of intravascular fibrin deposition through fibrinolysis and the inhibition of plasmin activity; this has found widespread clinical use in reducing perioperative bleeding. Increasing evidence also suggests diverse, although currently less defined, roles for plasmin in a number of physiological and pathological processes relating to extracellular matrix degradation, cell migration and tissue remodelling. In particular, dysregulation of plasmin has been linked to cancer invasion/metastasis and various chronic inflammatory conditions; this has prompted efforts to develop inhibitors of this protease. Although a number of plasmin inhibitors exist, they commonly suffer from poor potency and/or specificity of inhibition that either results in reduced efficacy or prevents clinical use. Consequently, there is a need for further development of high-affinity plasmin inhibitors that maintain selectivity over other serine proteases. This review summarises clearly defined and potential applications for plasmin inhibition. The properties of naturally occurring and engineered plasmin inhibitors are discussed in the context of current knowledge regarding plasmin structure, specificity and function. This includes design strategies to obtain the potency and specificity of inhibition in addition to controlled temporal and spatial distribution tailored for the intended use.

Hot spot analysis for driving the development of hits into leads in fragment-based drug discovery

Fragment-based drug design (FBDD) starts with finding fragment-sized compounds that are highly ligand efficient and can serve as a core moiety for developing high-affinity leads. Although the core-bound structure of a protein facilitates the construction of leads, effective design is far from straightforward. We show that protein mapping, a computational method developed to find binding hot spots and implemented as the FTMap server, provides information that complements the fragment screening results and can drive the evolution of core fragments into larger leads with a minimal loss or, in some cases, even a gain in ligand efficiency. The method places small molecular probes, the size of organic solvents, on a dense grid around the protein and identifies the hot spots as consensus clusters formed by clusters of several probes. The hot spots are ranked based on the number of probe clusters, which predicts the binding propensity of the subsites and hence their importance for drug design. Accordingly, with a single exception the main hot spot identified by FTMap binds the core compound found by fragment screening. The most useful information is provided by the neighboring secondary hot spots, indicating the regions where the core can be extended to increase its affinity. To quantify this information, we calculate the density of probes from mapping, which describes the binding propensity at each point, and show that the change in the correlation between a ligand position and the probe density upon extending or repositioning the core moiety predicts the expected change in ligand efficiency.

Contribution of explicit solvent effects to the binding affinity of small-molecule inhibitors in blood coagulation factor serine proteases

The prevention of blood coagulation is important in treating thromboembolic disorders, and several serine proteases involved in the coagulation cascade have been classified as pharmaceutically relevant. Whereas structure-based drug design has contributed to the development of some serine protease inhibitors, traditional computational methods have not been able to fully describe structure-activity relationships (SAR). Here, we study the SAR for a number of serine proteases by using a method that calculates the thermodynamic properties (enthalpy and entropy) of the water that solvates the active site. We show that the displacement of water from specific subpockets (such as S1-4 and the ester binding pocket) of the active site by the ligand can govern potency, especially for cases in which small chemical changes (i.e., a methyl group or halogen) result in a substantial increase in potency. Furthermore, we describe how relative binding free energies can be estimated by combining the water displacement energy with complementary terms from an implicit solvent molecular mechanics description binding.

The plasmin-antiplasmin system: structural and functional aspects

The plasmin-antiplasmin system plays a key role in blood coagulation and fibrinolysis. Plasmin and α(2)-antiplasmin are primarily responsible for a controlled and regulated dissolution of the fibrin polymers into soluble fragments. However, besides plasmin(ogen) and α(2)-antiplasmin the system contains a series of specific activators and inhibitors. The main physiological activators of plasminogen are tissue-type plasminogen activator, which is mainly involved in the dissolution of the fibrin polymers by plasmin, and urokinase-type plasminogen activator, which is primarily responsible for the generation of plasmin activity in the intercellular space. Both activators are multidomain serine proteases. Besides the main physiological inhibitor α(2)-antiplasmin, the plasmin-antiplasmin system is also regulated by the general protease inhibitor α(2)-macroglobulin, a member of the protease inhibitor I39 family. The activity of the plasminogen activators is primarily regulated by the plasminogen activator inhibitors 1 and 2, members of the serine protease inhibitor superfamily.

Application of a polarizable force field to calculations of relative protein-ligand binding affinities

An explicitly polarizable force field based exclusively on quantum data is applied to calculations of relative binding affinities of ligands to proteins. Five ligands, differing by replacement of an atom or functional group, in complexes with three serine proteases-trypsin, thrombin, and urokinase-type plasminogen activator-with available experimental binding data are used as test systems. A special protocol of thermodynamic integration was developed and used to provide sufficiently low levels of systematic error along with high numerical efficiency and statistical stability. The calculated results are in excellent quantitative (rmsd = 1.0 kcal/mol) and qualitative (R(2) = 0.90) agreement with experimental data. The potential of the methodology to explain the observed differences in the ligand affinities is also demonstrated.

Identification of a novel inhibitor of urokinase-type plasminogen activator

Urokinase-type plasminogen activator (uPA), a highly restricted serine protease, plays an important role in the regulation of diverse physiologic and pathologic processes. Strong clinical and experimental evidence has shown that elevated uPA expression is associated with cancer progression, metastasis, and shortened survival in patients. uPA has been considered as a promising molecular target for development of anticancer drugs. Here, we report the identification of several new uPA inhibitors using a high-throughput screen from a chemical library. From these uPA inhibitors, molecular modeling and docking studies identified 4-oxazolidinone as a novel lead pharmacophore. Optimization of the 4-oxazolidinone pharmacophore resulted in a series of structurally modified compounds with improved potency and selectivity. One of the 4-oxazolidinone analogues, UK122, showed the highest inhibition of uPA activity. The IC(50) of UK122 in a cell-free indirect uPA assay is 0.2 micromol/L. This compound also showed no or little inhibition of other serine proteases such as thrombin, trypsin, plasmin, and the tissue-type plasminogen activator, indicating its high specificity against uPA. Moreover, UK122 showed little cytotoxicity against CFPAC-1 cells (IC(50) >100 micromol/L) but significantly inhibited the migration and invasion of this pancreatic cancer cell line. Our data show that UK122 could potentially be developed as a new anticancer agent that prevents the invasion and metastasis of pancreatic cancer.

Substituted thiophene-anthranilamides as potent inhibitors of human factor Xa

A series of thiophene-containing non-amidine factor Xa inhibitors is described. Simple methyl-substituted thiophene analogs were relatively weak inhibitors. However, introduction of hydrophilic substituents at C-4 or C-5 of the thiophene afforded inhibitors with low nanomolar potency. Optimization of the thiophene substituent at C-4 afforded subnanomolar inhibitors with improved in vitro anticoagulant activity. Incorporating basic amine substituents on the thiophene increased hydrophilicity and improved anticoagulant activity. The pharmacokinetic profile of one inhibitor was evaluated in dogs, and the X-ray crystal structure of this compound bound to factor Xa provides insight into the observed SAR for binding to factor Xa.

A fluorine scan at the catalytic center of thrombin: C--F, C--OH, and C--OMe bioisosterism and fluorine effects on pKa and log D values

A series of 16 tricyclic thrombin inhibitors was prepared by using the 1,3-dipolar cycloaddition of azomethine ylides derived from 3- or 4-hydroxyproline and 4-bromobenzaldehyde, with N-(4-fluorobenzyl)maleimide as the key step. The terminal pyrrolidine ring of the inhibitors was systematically substituted to explore the potential bioisosteric behavior of C-F, C-OH, and C-OMe residues pointing into the environment of the catalytic center of a serine protease. X-ray crystal structure analyses revealed a distinct puckering preference of this ring. Substitution by F, HO, and MeO has a strong effect on the basicity of the adjacent pyrrolidine nitrogen center which originates from two sigma-inductive pathways between this center and the electronegative O and F atoms. gem-Difluorination decreases the pKa value of this tertiary amine center to <2, making the conjugated ammonium ion a moderately strong acid. Unexpectedly, F substitution next to the nitrogen center reduced the lipophilicity of the ligands, as revealed by measurements of the logarithmic partition coefficient log D. The biological assays showed that all compounds are thrombin inhibitors with activities between Ki=0.08 and 2.17 microM. Bioisosteric behavior of F, HO, and MeO substituents was observed. Their electronegative F and O atoms undergo energetically similar polar interactions with positively polarized centers, such as the N atom of His 57 which is hydrogen bonded to the catalytic Ser 195. However, for energetically similar polar interactions of C-F, C-OH, and C-OMe to occur, sufficient space is necessary for the accommodation of the Me group of the C-OMe residue, and a H-bond acceptor must be present to prevent unfavorable desolvation of the C-OH residue.

Mapping the Fluorophilicity of a Hydrophobic Pocket: Synthesis and Biological Evaluation of Tricyclic Thrombin Inhibitors Directing Fluorinated Alkyl Groups into the P Pocket

In the completion of our fluorine scan of tricyclic inhibitors to map the fluorophilicity/fluorophobicity of the thrombin active site, a series of 11 new ligands featuring alkyl, alkenyl, and fluoroalkyl groups was prepared to explore fluorine effects on binding into the hydrophobic proximal (P) pocket, lined by Tyr 60A and Trp 60D, His 57, and Leu 99. The synthesis of the tricyclic scaffolds was based on the 1,3-dipolar cycloaddition of azomethine ylides, derived from L-proline and 4-bromobenzaldehyde, with N-(4-fluorobenzyl)maleimide. Introduction of alkyl, alkenyl, and partially fluorinated alkyl residues was achieved upon substitution of a sulfonyl group by mixed Mg/Zn organometallics followed by oxidation/deoxyfluorination, as well as oxidation/reduction/deoxyfluorination sequences. In contrast, the incorporation of perfluoroalkyl groups required a stereoselective nucleophilic addition reaction at the "upper" carbonyl group of the tricycles, thereby yielding scaffolds with an additional OH, F, or OMe group, respectively. All newly prepared inhibitors showed potent biological activity, with inhibitory constants (K(i) values) in the range of 0.008-0.163 microM. The X-ray crystal structure of a protein-ligand complex revealed the exact positioning of a difluoromethyl substituent in the tight P pocket. Fluorophilic characteristics are attributed to this hydrophobic pocket, although the potency of the inhibitors was found to be modulated by steric rather than electronic factors.

Exploring structurally conserved solvent sites in protein families

Protein-bound water molecules are important components of protein structure, and therefore, protein function and energetics. Although structural conservation of solvent has been studied in a few protein families, a lack of suitable computational tools has hindered more comprehensive analyses. Herein we present a semiautomated computational approach for identifying solvent sites that are conserved among proteins sharing a common three-dimensional structure. This method is tested on six protein families: (1) monodomain cytochrome c, (2) fatty-acid binding protein, (3) lactate/malate dehydrogenase, (4) parvalbumin, (5) phospholipase A2, and (6) serine protease. For each family, the method successfully identified previously known conserved solvent sites. Moreover, the method discovered 22 novel conserved solvent sites, some of which have higher degrees of conservation than the previously known sites. All six families studied had solvent sites with more than 90% conservation and these sites were invariably located in regions of the protein with very high sequence conservation. These results suggest that highly conserved solvent sites, by virtue of their proximity to conserved residues, should be considered as one of the defining three-dimensional structural characteristics of protein families and folds.

A Urokinase-type Plasminogen Activator-inhibiting Cyclic Peptide with an Unusual P2 Residue and an Extended Protease Binding Surface Demonstrates New Modalities for Enzyme Inhibition

To find new principles for inhibiting serine proteases, we screened phage-displayed random peptide repertoires with urokinase-type plasminogen activator (uPA) as the target. The most frequent of the isolated phage clones contained the disulfide bridge-constrained sequence CSWRGLENHRMC, which we designated upain-1. When expressed recombinantly with a protein fusion partner, upain-1 inhibited the enzymatic activity of uPA competitively with a temperature and pH-dependent K(i), which at 25 degrees C and pH 7.4 was approximately 500 nm. At the same conditions, the equilibrium dissociation constant K(D), monitored by displacement of p-aminobenzamidine from the specificity pocket of uPA, was approximately 400 nm. By an inhibitory screen against other serine proteases, including trypsin, upain-1 was found to be highly selective for uPA. The cyclical structure of upain-1 was indispensable for uPA binding. Alanine-scanning mutagenesis identified Arg(4) of upain-1 as the P(1) residue and indicated an extended binding interaction including the specificity pocket and the 37-, 60-, and 97-loops of uPA and the P(1), P(2), P(3)', P(4)', and the P(5)' residues of upain-1. Substitution with alanine of the P(2) residue, Trp(3), converted upain-1 into a distinct, although poor, uPA substrate. Upain-1 represents a new type of uPA inhibitor that achieves selectivity by targeting uPA-specific surface loops. Most likely, the inhibitory activity depends on its cyclical structure and the unusual P(2) residue preventing the scissile bond from assuming a tetrahedral geometry and thus from undergoing hydrolysis. Peptide-derived inhibitors such as upain-1 may provide novel mechanistic information about enzyme-inhibitor interactions and alternative methodologies for designing effective protease inhibitors.

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.

Synthesis of a Hydroxyethylene Isostere of the Tripeptide Arg-Gly-Leu via a Convergent Acyl-like Radical Addition Strategy

[reaction: see text] A hydroxyethylene isostere of the tripeptide Arg-Gly-Leu, representing an important fragment of a novel cyclic-peptide-based uPA inhibitor, was synthesized in few steps employing as the key step a samarium diiodide promoted coupling of either the 4-thiopyridyl ester of N(alpha)-Fmoc- or N(alpha)-Cbz-protected L-ornithine with the N-acryloyl derivative of L-leucine methyl ester. Epimerization under the coupling conditions at the chiral center in the alpha-position to the ketone was demonstrated not to take place. A stereoselective reduction of the Cbz-protected aminoketone obtained from this radical reaction was promoted by the same single-electron reducing agent in the presence of methanol providing the syn-amino alcohol with a diastereoselectivity of 85:15. With the use of lithium tri-tert-butoxyaluminum hydride in methanol, the corresponding anti-isomer was obtained almost exclusively. Subsequent elaboration of the ornithine moiety in the anti-isomer by introduction of the guanidine group followed by hydrolysis of the C-terminal ester bond and protection of the alcohol as its tert-butyldimethylsilyl ether provided the desired tripeptide mimic. The long reaction times required for the radical addition reactions with N(delta)-Boc-L-ornithine (up to 5 days) led to a short study where a series of 4-thiopyridyl esters of Cbz-protected amino acids were reacted with two acrylates. Whereas N(delta)-Boc-L-ornithine, alanine, phenylalanine, proline, and leucine all provided the aminoketone in 43-79% yield, valine only afforded traces of the coupling product.

Crystal Structures of the FXIa Catalytic Domain in Complex with Ecotin Mutants Reveal Substrate-like Interactions

Thrombosis can lead to life-threatening conditions such as acute myocardial infarction, pulmonary embolism, and stroke. Although commonly used anti-coagulant drugs, such as low molecular weight heparin and warfarin, are effective, they carry a significant risk of inducing severe bleeding complications, and there is a need for safer drugs. Activated Factor XI (FXIa) is a key enzyme in the amplification phase of the coagulation cascade. Anti-human FXI antibody significantly reduces thrombus growth in a baboon thrombosis model without bleeding problems (Gruber, A., and Hanson, S. R. (2003) Blood 102, 953-955). Therefore, FXIa is a potential target for anti-thrombosis therapy. To determine the structure of FXIa, we derived a recombinant catalytic domain of FXI, consisting of residues 370-607 (rhFXI370-607). Here we report the first crystal structure of rhFXI370-607 in complex with a substitution mutant of ecotin, a panserine protease protein inhibitor secreted by Escherichia coli, to 2.2 A resolution. The presence of ecotin not only assisted in the crystallization of the enzyme but also revealed unique structural features in the active site of FXIa. Subsequently, the sequence from P5 to P2' in ecotin was mutated to the FXIa substrate sequence, and the structures of the rhFXI370-607-ecotin mutant complexes were determined. These structures provide us with an understanding of substrate binding interactions of FXIa, the structural information essential for the structure-based design of FXIa-selective inhibitors.

S1 subsite in snake venom thrombin-like enzymes: can S1 subsite lipophilicity be used to sort binding affinities of trypsin-like enzymes to small-molecule inhibitors?

Thrombin-like enzymes isolated from snake venoms comprise a group of serine proteinases responsible for many important coagulation disorders in the envenomed victims. Besides, these proteinases have great biotechnological interest as antithrombotic agents and as diagnostic tools. However, in spite of the recent overflow of snake venom thrombin-like enzymes (SVTLEs) on protein sequence databases, there is a lack of three-dimensional (3D) structural information on this family. Without such 3D structures available many aspects of the biological function and biochemical properties of these enzymes still remain obscure. Therefore, we have gone through a series of computational techniques, which enabled us to identify the set of residues involved in molecular recognition of inhibitors bound to the S1 subsite of snake venom thrombin-like enzymes (SVTLEs) and ultimately conclude that nonpolar (van der Waals) intermolecular interactions and ligand's hydrophobicity are the most important factors affecting binding affinities to the S1 subsite of a SVTLE isolated from the venom of Lachesis muta muta (Lmm-TLE). Consequently, we have proposed that S1 subsite lipophilicity may be used to sort binding affinities of trypsin-like enzymes to small molecules by showing that the inhibitory potency of several S1-directed compounds follows subsite lipophilicity among Lmm-TLE and other three homologous proteases. Noteworthy, in the course of our analyses we determined that thrombin's S1 subsite should, in fact, be considered less lipophilic than that of trypsin if we account for the presence of the sodium-controlled water channel communicating with the S1 subsite in the coagulant enzyme.

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.

A fluorine scan of the phenylamidinium needle of tricyclic thrombin inhibitors: effects of fluorine substitution on pKa and binding affinity and evidence for intermolecular C-F...CN interactions

The H-atoms of the phenylamidinium needle of tricyclic thrombin inhibitors, which interacts with Asp189 at the bottom of the selectivity pocket S1 of the enzyme, were systematically exchanged with F-atoms in an attempt to improve the pharmacokinetic properties by lowering the pK(a) value. Both the pK(a) values and the inhibitory constants K(i) against thrombin and trypsin were decreased upon F-substitution. Interestingly, linear free energy relationships (LFERs) revealed that binding affinity against thrombin is much more affected by a decrease in pK(a) than the affinity against trypsin. Surprising effects of F-substitutions in the phenylamidinium needle on the pK(a) value of the tertiary amine centre in the tricyclic scaffold of the inhibitors were observed and subsequently rationalised by X-ray crystallographic analysis and ab initio calculations. Evidence for highly directional intermolecular C-F...CN interactions was obtained by analysis of small-molecule X-ray crystal structures and investigations in the Cambridge Structural Database (CSD).

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.

Azaindoles: moderately basic P1 groups for enhancing the selectivity of thrombin inhibitors

Starting from a 2-amino-6-methylpyridine P1 group and following a strategy of enlarging it whilst reducing its polarity, we have developed a series of potent, moderately basic azaindoles which are intrinsically much more selective for thrombin versus trypsin. Certain pyrazinone acetamide azaindole derivatives have pharmacokinetic parameters after oral administration to dogs, or efficacy in vitro, comparable to an optimized pyrazinone acetamide 2-amino-6-methylpyridine derivative.

Reconstructing the binding site of factor Xa in trypsin reveals ligand-induced structural plasticity

In order to investigate issues of selectivity and specificity in protein-ligand interactions, we have undertaken the reconstruction of the binding pocket of human factor Xa in the structurally related rat trypsin by site-directed mutagenesis. Three sequential regions (the "99"-, the "175"- and the "190"- loops) were selected as representing the major structural differences between the ligand binding sites of the two enzymes. Wild-type rat trypsin and variants X99rT and X(99/175/190)rT were expressed in yeast, and analysed for their interaction with factor Xa and trypsin inhibitors. For most of the inhibitors studied, progressive loop replacement at the trypsin surface resulted in inhibitory profiles akin to factor Xa. Crystals of the variants were obtained in the presence of benzamidine (3), and could be soaked with the highly specific factor Xa inhibitor (1). Binding of the latter to X99rT results in a series of structural adaptations to the ligand, including the establishment of an "aromatic box" characteristic of factor Xa. In X(99/175/190)rT, introduction of the 175-loop results in a surprising re-orientation of the "intermediate helix", otherwise common to trypsin and factor Xa. The re-orientation is accompanied by an isomerisation of the Cys168-Cys182 disulphide bond, and burial of the critical Phe174 side-chain. In the presence of (1), a major re-organisation of the binding site takes place to yield a geometry identical to that of factor Xa. In all, binding of (1) to trypsin and its variants results in significant structural rearrangements, inducing a binding surface strongly reminiscent of factor Xa, against which the inhibitor was optimised. The structural data reveal a plasticity of the intermediate helix, which has been implicated in the functional cofactor dependency of many trypsin-like serine proteinases. This approach of grafting loops onto scaffolds of known related structures may serve to bridge the gap between structural genomics and drug design.

Function first: a powerful approach to post-genomic drug discovery

In the post-genomic era, pharmaceutical researchers must evaluate vast numbers of protein sequences and formulate novel, intelligent strategies for identifying valid targets and discovering leads against them. The identification of small molecules that selectively target proteins or protein families will be aided by knowing the function and/or the structure of the target(s). By identifying protein function first, efficiencies are gained that allow subsequent focus of resources on particular protein families of interest. This article reviews current proteomic-scale approaches to identifying function as a way of accelerating lead discovery.

2-(2-Hydroxy-3-alkoxyphenyl)-1H-benzimidazole-5-carboxamidine derivatives as potent and selective urokinase-type plasminogen activator inhibitors

The development of potent and selective urokinase-type plasminogen activator (uPA) inhibitors based on the lead molecule 2-(2-hydroxy-3-ethoxyphenyl)-1H-benzimidazole-5-carboxamidine (3a) is described.

4-Aminoarylguanidine and 4-aminobenzamidine derivatives as potent and selective urokinase-type plasminogen activator inhibitors

The structure-based design of potent and selective urokinase-type plasminogen activator (uPA) inhibitors with 4-aminoarylamidine or 4-aminoarylguanidine S1 binding groups, is described.