Oxyanion-mediated inhibition of serine proteases

Crystal structure of the Michaelis complex of trypsin with N‐α‐benzoyl‐l‐arginine ethyl ester

The crystal structure of Michaelis complex of the bovine trypsin with N‐α‐benzoyl‐l‐arginine ethyl ester (BAEE) was determined at 2.30 Å resolution. The structure of enzyme‐substrate complex was solved in space group H32. The active site of trypsin was found to be occupied with the N‐α‐benzoyl‐l‐arginine ethyl ester. The hydrolyzed product of substrate molecules was also crystallized with trypsin. The substrate was embedded within the S1 binding site of the enzyme, and showed contacts with Gly‐195, Ser‐216, and Ser‐192. Some water molecules were also found within the vicinity of catalytic residues of enzyme. Interestingly, the hydrolyzed product present within the crystal lattice was found to be with inverted configuration. The crystal structure of trypsin in‐complex with N‐α‐benzoyl‐l‐arginine ethyl ester has never been reported previously.

Bisbenzamidine and Bisbenzguanidine Ureas Act as Antibacterial Agents against Pseudomonas aeruginosa

Due to the global rise in the number of antibiotic resistant bacterial infections over the past 20 years, there is a dire need for the development of small molecule antibiotics capable of overcoming resistance mechanisms in pathogenic bacteria. Antibiotic development against Gram-negative pathogens, such as Pseudomonas aeruginosa, is especially challenging due to their additional outer membrane which reduces antibiotic entry. Recently, it has been shown that a broad range of nitrogen functionality, including guanidines, amidines, primary amines, imidazolines, and imidazoles, promote antibiotic and adjuvant activity in Gram-negative bacteria, but few of these have been targeted towards Pseudomonas aeruginosa specifically despite this pathogen being deemed a critical threat by the United States Centers for Disease Control and Prevention. Herein, we examined a small series of known and unknown nitrogenous dimers, with guanidine, amidine, dimethyl amine, and pyridine functionality, for antibacterial activity against multidrug resistant Pseudomonas aeruginosa. We found that two, with bisbenzguanidine and bisbenzamidine functionality, are potent against clinical isolates of multidrug resistant and biofilm forming Pseudomonas aeruginosa.

Factor XIa Inhibitors in Anticoagulation Therapy: Recent Advances and Perspectives

Factor XIa (FXIa) in the intrinsic pathway of the coagulation process has been proven to be an effective and safe target for anticoagulant discovery with limited or no bleeding. Numerous small-molecule FXIa inhibitors (SMFIs) with various scaffolds have been identified in the early stages of drug discovery. They have served as the foundation for the recent discovery of additional promising SMFIs with improved potency, selectivity, and pharmacokinetic profiles, some of which have entered clinical trials for the treatment of thrombosis. After reviewing the coagulation process and structure of FXIa, this perspective discusses the rational or structure-based design, discovery, structure-activity relationships, and development of SMFIs disclosed in recent years. Strategies for identifying more selective and druggable SMFIs are provided, paving the way for the design and discovery of more useful SMFIs for anticoagulation therapy.

Synthesis of 2-guanidinyl pyridines and their trypsin inhibition and docking

A range of guanidine-based pyridines, and related compounds, have been prepared (19 examples). These compounds were evaluated in relation to their competitive inhibition of bovine pancreatic trypsin. Results demonstrate that compounds in which the guanidinyl substituent can form an intramolecular hydrogen bond (IMHB) with the pyridinyl nitrogen atom (6a-p) are better trypsin inhibitors than their counterparts (10-13) that are unable to form an IMHB. Among the compounds 6a-p, examples containing a 5-halo substituent were, generally, found to be better trypsin inhibitors. This trend was inversely related to electronegativity, thus, 1-(5-iodopyridin-2-yl)guanidinium ion 6e (K_i = 0.0151 mM) was the optimum inhibitor in the 5-halo series. Amongst the isomeric methyl substituted compounds, 1-(3-methylpyridin-2-yl)guanidinium ion 6h demonstrated optimum levels of trypsin inhibition (K_i = 0.0140 mM). In order to rationalise the measured enzyme inhibition, selected compounds were docked with bovine and human trypsin with a view to understanding active site occupancy and taken together with the K_i values the order of inhibitory ability suggests that the 5-halo 2-guanidinyl pyridine inhibitors form a halogen bond with the catalytically active serine hydroxy group.

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

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

FLIP: An assisting software in structure based drug design using fingerprint of protein-ligand interaction profiles

With the growing number of labor-intensive data in the pharmaceutical industries and public domain for protein-ligand complexes, a significant challenge is still remaining in managing and leveraging this vast information. Here, a standalone application is presented for analysis, organization, and illustration of structural data and molecular interactions for exploiting 3D-structures into simple 1D fingerprints. The utility of the approach was shown in unraveling a feasible solution for post-processing of docking results in parallel with providing fruitful analysis for users in order to investigate molecular interactions. Remarkably, all interaction possibilities including (hydrogen bond, water-bridged, electrostatic, and hydrophobic as well as π- π and cation-π interactions) are supported both in the form of fingerprints and compelling reports. These investigations are mainly considered based on right orientation, location, and geometry of the interacting pairs rather than the acquisition of the energy terms. The reasonable efficiency of our application in different models was comparable to recent methods It is clearly presented that FLIP provides a faster way to generate usable fingerprints for ligand and protein binding modes. FLIP is free for academic use and is available at: http://zistrayan.com/development/download/flip/package.zip.

PharmDock: a pharmacophore-based docking program

Background

Protein-based pharmacophore models are enriched with the information of potential interactions between ligands and the protein target. We have shown in a previous study that protein-based pharmacophore models can be applied for ligand pose prediction and pose ranking. In this publication, we present a new pharmacophore-based docking program PharmDock that combines pose sampling and ranking based on optimized protein-based pharmacophore models with local optimization using an empirical scoring function.

Results

Tests of PharmDock on ligand pose prediction, binding affinity estimation, compound ranking and virtual screening yielded comparable or better performance to existing and widely used docking programs. The docking program comes with an easy-to-use GUI within PyMOL. Two features have been incorporated in the program suite that allow for user-defined guidance of the docking process based on previous experimental data. Docking with those features demonstrated superior performance compared to unbiased docking.

Conclusion

A protein pharmacophore-based docking program, PharmDock, has been made available with a PyMOL plugin. PharmDock and the PyMOL plugin are freely available from http://people.pharmacy.purdue.edu/~mlill/software/pharmdock.

Ureasil-poly(ethylene oxide) hybrid matrix for selective adsorption and separation of dyes from water

Herein, we present a cross-linked ureasil-polyether-siloxane hybrid (labeled PEO500) that can function as a stimuli-sensitive material; it swells or shrinks in response to changes in the environmental conditions and it can also, effectively and selectively, remove dyes from water solution. We also developed a methodology to separate a mixture of cationic and anionic dyes present in water. Addition of PEO500 to an aqueous solution of the anionic orange II (OII) or the ponceau S (PS) dye rendered the solution colorless, but an aqueous solution of cationic methylene blue (MB) remained unchanged after 2 h of contact with the insoluble matrix. In situ small-angle X-ray scattering (SAXS) showed that the distance of siloxane nanodomains are strongly affected by the swelling or shriking. By in situ UV-vis adsorption experiments, we found that the kinetics of OII and PS removal followed a pseudo-first-order rate equation. We accomplished B3LYP calculations, to establish which sites on the matrix interacted with the dyes and to investigate the nature of the matrix-dye chemical bonds. On the basis of the experimental and theoretical investigations, we proposed some mechanisms to explain how PEO500 adsorbs anionic dyes efficiently. This "smart" matrix is potentially applicable as an efficient, fast, selective, and convenient device in water treatment and stimuli-sensitive response materials.

Implementation of the Hungarian algorithm to account for ligand symmetry and similarity in structure-based design

False negative docking outcomes for highly symmetric molecules are a barrier to the accurate evaluation of docking programs, scoring functions, and protocols. This work describes an implementation of a symmetry-corrected root-mean-square deviation (RMSD) method into the program DOCK based on the Hungarian algorithm for solving the minimum assignment problem, which dynamically assigns atom correspondence in molecules with symmetry. The algorithm adds only a trivial amount of computation time to the RMSD calculations and is shown to increase the reported overall docking success rate by approximately 5% when tested over 1043 receptor–ligand systems. For some families of protein systems the results are even more dramatic, with success rate increases up to 16.7%. Several additional applications of the method are also presented including as a pairwise similarity metric to compare molecules during de novo design, as a scoring function to rank-order virtual screening results, and for the analysis of trajectories from molecular dynamics simulation. The new method, including source code, is available to registered users of DOCK6 (http://dock.compbio.ucsf.edu).

A novel and efficient ligand-based virtual screening approach using the HWZ scoring function and an enhanced shape-density model

In this work, we extend our previous ligand shape-based virtual screening approach by using the scoring function Hamza-Wei-Zhan (HWZ) score and an enhanced molecular shape-density model for the ligands. The performance of the method has been tested against the 40 targets in the Database of Useful Decoys and compared with the performance of our previous HWZ score method. The virtual screening results using the novel ligand shape-based approach demonstrated a favorable improvement (area under the receiver operator characteristics curve AUC = .89 ± .02) and effectiveness (hit rate HR(1%) = 53.0% ± 6.3 and HR(10%) = 71.1% ± 4.9). The comparison of the overall performance of our ligand shape-based method with the highest ligand shape-based virtual screening approach using the data fusion of multi queries showed that our strategy takes into account deeper the chemical information of the set of active ligands. Furthermore, the results indicated that our method are suitable for virtual screening and yields superior prediction accuracy than the other study derived from the data fusion using five queries. Therefore, our novel ligand shape-based screening method constitutes a robust and efficient approach to the 3D similarity screening of small compounds and open the door to a whole new approach to drug design by implementing the method in the structure-based virtual screening.

Contact-based ligand-clustering approach for the identification of active compounds in virtual screening

Evaluation of docking results is one of the most important problems for virtual screening and in silico drug design. Modern approaches for the identification of active compounds in a large data set of docked molecules use energy scoring functions. One of the general and most significant limitations of these methods relates to inaccurate binding energy estimation, which results in false scoring of docked compounds. Automatic analysis of poses using self-organizing maps (AuPosSOM) represents an alternative approach for the evaluation of docking results based on the clustering of compounds by the similarity of their contacts with the receptor. A scoring function was developed for the identification of the active compounds in the AuPosSOM clustered dataset. In addition, the AuPosSOM efficiency for the clustering of compounds and the identification of key contacts considered as important for its activity, were also improved. Benchmark tests for several targets revealed that together with the developed scoring function, AuPosSOM represents a good alternative to the energy-based scoring functions for the evaluation of docking results.

Prediction of trypsin/molecular fragment binding affinities by free energy decomposition and empirical scores

Two families of binding affinity estimation methodologies are described which were utilized in the SAMPL3 trypsin/fragment binding affinity challenge. The first is a free energy decomposition scheme based on a thermodynamic cycle, which included separate contributions from enthalpy and entropy of binding as well as a solvent contribution. Enthalpic contributions were estimated with PM6-DH2 semiempirical quantum mechanical interaction energies, which were modified with a statistical error correction procedure. Entropic contributions were estimated with the rigid-rotor harmonic approximation, and solvent contributions to the free energy were estimated with several different methods. The second general methodology is the empirical score LISA, which contains several physics-based terms trained with the large PDBBind database of protein/ligand complexes. Here we also introduce LISA+, an updated version of LISA which, prior to scoring, classifies systems into one of four classes based on a ligand's hydrophobicity and molecular weight. Each version of the two methodologies (a total of 11 methods) was trained against a compiled set of known trypsin binders available in the Protein Data Bank to yield scaling parameters for linear regression models. Both raw and scaled scores were submitted to SAMPL3. Variants of LISA showed relatively low absolute errors but also low correlation with experiment, while the free energy decomposition methods had modest success when scaling factors were included. Nonetheless, re-scaled LISA yielded the best predictions in the challenge in terms of RMS error, and six of these models placed in the top ten best predictions by RMS error. This work highlights some of the difficulties of predicting binding affinities of small molecular fragments to protein receptors as well as the benefit of using training data.

Implementation and evaluation of a docking-rescoring method using molecular footprint comparisons

A docking-rescoring method, based on per-residue van der Waals (VDW), electrostatic (ES), or hydrogen bond (HB) energies has been developed to aid discovery of ligands that have interaction signatures with a target (footprints) similar to that of a reference. Biologically useful references could include known drugs, inhibitors, substrates, transition states, or side-chains that mediate protein-protein interactions. Termed footprint similarity (FPS) score, the method, as implemented in the program DOCK, was validated and characterized using: (1) pose identification, (2) crossdocking, (3) enrichment, and (4) virtual screening. Improvements in pose identification (6–12%) were obtained using footprint-based (FPS(VDW+ES)) vs. standard DOCK (DCE(VDW+ES)) scoring as evaluated on three large datasets (680–775 systems) from the SB2010 database. Enhanced pose identification was also observed using FPS (45.4% or 70.9%) compared with DCE (17.8%) methods to rank challenging crossdocking ensembles from carbonic anhydrase. Enrichment tests, for three representative systems, revealed FPSVDW+ES scoring yields significant early fold enrichment in the top 10% of ranked databases. For EGFR, top FPS poses are nicely accommodated in the molecular envelope defined by the reference in comparison with DCE, which yields distinct molecular weight bias toward larger molecules. Results from a representative virtual screen of ca. 1 million compounds additionally illustrate how ligands with footprints similar to a known inhibitor can readily be identified from within large commercially available databases. By providing an alternative way to rank ligand poses in a simple yet directed manner we anticipate that FPS scoring will be a useful tool for docking and structure-based design.

Structural binding evidence of the trypanocidal drugs berenil and pentacarinate active principles to a serine protease model

Bovine trypsin is a model system for the serine protease class of enzymes, which is an important target for contemporary medicinal chemistry. Some structural and thermodynamic reports are available on its interaction with benzamidine-based compounds but no structural information is available so far on its binding modes to the active principles of the trypanocidal drugs Pentacarinate (pentamidine) and Berenil (diminazene). The crystallographic structures of bovine beta-trypsin in complex with the ligands were determined to a resolution of 1.57 A (diminazene) and 1.70 A (diminazene and pentamidine). The second benzamidine moieties in these inhibitors are bound to the enzyme in different hot spots and only few hydrogen bonds mediate these interactions. Thermodynamic parameters for the association of pentamidine with beta-trypsin reveal that this inhibitor has about 1.3-fold lower affinity than diminazene. Moreover its binding mode resembles other benzamidine-based compounds that assess the aryl binding pocket of the enzyme; however, with almost 2.5-fold higher affinity. This is the first structural evidence of the binding of Berenil and Pentacarinate active principles trypanocidal drugs to serine proteases.

Hydration of acetone in the gas phase and in water solvent

In water solvent, the hydration of acetone proceeds by a cyclic (cooperative) process in which concurrent C–O bond formation and proton transfer to oxygen take place through a solvent and (or) catalyst bridge. Reactivity is determined primarily by the concentration of a reactant complex and not the barrier from this complex. This situation is reversed in the gas phase; although the concentrations of reactive complexes are much higher than in solution, the barriers are also higher and dominant in determining reactivity. Calculations of isotope effects suggest that multiple hydron transfers are synchronous in the gas phase to avoid zwitterionic transition states. In solution, such transition states are stabilized by solvation and hydron transfers can be asynchronous.

Crystal structure of a glycosylated Fab from an IgM cryoglobulin with properties of a natural proteolytic antibody

The 2.6 A (1 A=0.1 nm) resolution structure has been determined for the glycosylated Fab (fragment antigen binding) of an IgM (Yvo) obtained from a subject with Waldenström's macroglobulinaemia. Dynamic light scattering was used to estimate the gel point and monitor the formation of an ordered hydroscopic gel of Yvo IgM upon cooling. If a cryoglobulin forms gels in peripheral tissues and organs, the associated swelling and damage to microvasculature can result in considerable morbidity and mortality. The three-dimensional structure of the branched N-linked oligosaccharide associated with the CH1 domain (first constant domain of heavy chain) is reported. The carbohydrate may act to shield part of the lateral surface of the CH1 domain and crowd the junction between the CH1 and CH2 domains, thereby limiting the segmental flexibility of the Fab arms in intact Yvo IgM, especially at low temperatures. Recently, Yvo IgM was shown to have the properties of a naturally occurring proteolytic antibody [Paul, Karle, Planque, Taguchi, Salas, Nishiyama, Handy, Hunter, Edmundson and Hanson (2004) J. Biol. Chem. 279, 39611-39619; Planque, Bangale, Song, Karle, Taguchi, Poindexter, Bick, Edmundson, Nishiyama and Paul (2004) J. Biol Chem. 279, 14024-14032]. The Yvo protein displayed the ability to cleave, by a nucleophilic mechanism, the amide bonds of a variety of serine protease substrates and the gp120 coat protein of HIV. An atypical serine, arginine and glutamate motif is located in the middle of the Yvo antigen-binding site and displays an overall geometry that mimics the classical serine, histidine and aspartate catalytic triad of serine proteases. Our present findings indicate that pre-existing or natural antibodies can utilize at least one novel strategy for the cleavage of peptide bonds.

Noncovalent inhibition of the serine proteases, alpha-chymotrypsin and trypsin by trifluoro(organo)borates

A series of potassium organotrifluoroborates were synthesized. Their stability to hydrolysis was determined in D2O, TRIS and phosphate buffer. It was found that in both D2O and TRIS buffers, these compounds are quite stable, whereas in phosphate buffer rapid hydrolysis occurs. Based on these results, a study was undertaken to determine whether potassium organotrifluoroborates can serve as protease inhibitors. It was found that potassium organotrifluoroborates increased inhibition by at least an order of magnitude over the corresponding boronates. Dixon plots showed that these compounds are reversible competitive inhibitors of alpha-chymotrypsin and trypsin. Based on 19F NMR, we speculate that they inactivate the enzymes as a result of the formation of hydrogen-bonds between fluorine atoms of the inhibitors and the serine protease.

Enzyme Binding Selectivity Prediction: α-Thrombin vs Trypsin Inhibition

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

Inhibition of arginine gingipains (RgpB and HRgpA) with benzamidine inhibitors: zinc increases inhibitory potency

We assayed several benzamidine derivatives for inhibition potency with HRgpA and RgpB gingipains, enzymes which are involved in the pathogenesis of gingivitis and periodontal disease. The benzamidine derivatives proved to be effective inhibitors of HRgpA and RgpB, with the best inhibitor being a bis-benzamidine with a urea linker (Ki=30 microM). The inhibition potency was increased 2-3 fold in the presence of low concentrations of zinc with the benzamidines containing a urea moiety linking the two aromatic rings. We propose an inhibition model involving a tetrahedral zinc atom coordinated with the active site Cys and His of gingipain and the urea linker in the benzamidine inhibitor. In summary, we have discovered a new series of effective inhibitors for the gingipains and found a novel way to increase inhibitor potency with the HRgpA and RgpB gingipains using zinc.

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

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

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.

Recruiting Zn2+ to mediate potent, specific inhibition of serine proteases

As regulators of ubiquitous biological processes, serine proteases can cause disease states when inappropriately expressed or regulated, and are thus rational targets for inhibition by drugs. Recently we described a new inhibition mechanism applicable for the development of potent, selective small molecule serine protease inhibitors that recruit physiological Zn2+ to mediate high affinity (sub-nanomolar) binding. To demonstrate some of the structural principles by which the selectivity of Zn2+-mediated serine protease inhibitors can be developed toward or against a particular target, here we determine and describe the structures of thrombin-BABIM-Zn2+, -keto-BABIM-Zn2+, and -hemi-BABIM-Zn2+ (where BABIM is bis(5-amidino-2-benzimidazolyl)methane, keto-BABIM is bis(5-amidino-2-benzimidazolyl)methane ketone, and hemi-BABIM is (5-amidino-2-benzimidazolyl)(2-benzimidazolyl)methane), and compare them with the corresponding trypsin-inhibitor-Zn2+ complexes. Inhibitor binding is mediated by a Zn ion tetrahedrally coordinated by two benzimidazole nitrogen atoms of the inhibitor, by N(epsilon2)His57, and by O(gamma)Ser195. The structures of Zn2+-free trypsin-BABIM and -hemi-BABIM were also determined at selected pH values for comparison with the corresponding Zn2+-mediated complexes. To assess some of the physiological parameters important for harnessing Zn2+ as a co-inhibitor, crystal structures at multiple pH and [Zn2+] values were determined for trypsin-keto-BABIM. The Kdvalue of Zn2+ for the binary trypsin-keto-BABIM complex was estimated to be <12 nM at pH 7.06 by crystallographic determination of the occupancy of bound Zn2+ in trypsin-keto-BABIM crystals soaked at this pH in synthetic mother liquor containing inhibitor and 100 nM Zn2+. In synthetic mother liquor saturated in Zn2+, trypsin-bound keto-BABIM is unhydrated at pH 9.00 and 9.93, and has an sp2 hybridized ketone carbon bridging the 5-amidinobenzimidazoles, whereas at pH 7.00 and 8.00 it undergoes hydration and a change in geometry upon addition of water to the bridging carbonyl group. To show how Zn2+ could be recruited as a co-inhibitor of other enzymes, a method was developed for locating in protein crystals Zn2+ binding sites where design of Zn2+-mediated ligands can be attempted. Thus, by soaking trypsin crystals in high concentrations of Zn2+ in the absence of a molecular inhibitor, the site where Zn2+ mediates binding of BABIM and analogs was identified, as well as another Zn2+ binding site.