Episelection: novel Ki approximately nanomolar inhibitors of serine proteases selected by binding or chemistry on an enzyme surface

Advances in computational methods for ligand binding kinetics

Binding kinetic parameters can be correlated with drug efficacy, which in recent years led to the development of various computational methods for predicting binding kinetic rates and gaining insight into protein-drug binding paths and mechanisms. In this review, we introduce and compare computational methods recently developed and applied to two systems, trypsin-benzamidine and kinase-inhibitor complexes. Methods involving enhanced sampling in molecular dynamics simulations or machine learning can be used not only to predict kinetic rates, but also to reveal factors modulating the duration of residence times, selectivity, and drug resistance to mutations. Methods which require less computational time to make predictions are highlighted, and suggestions to reduce the error of computed kinetic rates are presented.

The β-link motif in protein architecture

The β-link is a composite protein motif consisting of a G1β β-bulge and a type II β-turn, and is generally found at the end of two adjacent strands of antiparallel β-sheet. The 1,2-positions of the β-bulge are also the 3,4-positions of the β-turn, with the result that the N-terminal portion of the polypeptide chain is orientated at right angles to the β-sheet. Here, it is reported that the β-link is frequently found in certain protein folds of the SCOPe structural classification at specific locations where it connects a β-sheet to another area of a protein. It is found at locations where it connects one β-sheet to another in the β-sandwich and related structures, and in small (four-, five- or six-stranded) β-barrels, where it connects two β-strands through the polypeptide chain that crosses an open end of the barrel. It is not found in larger (eight-stranded or more) β-barrels that are straightforward β-meanders. In some cases it initiates a connection between a single β-sheet and an α-helix. The β-link also provides a framework for catalysis in serine proteases, where the catalytic serine is part of a conserved β-link, and in cysteine proteases, including M^pro of human SARS-CoV-2, in which two residues of the active site are located in a conserved β-link.

Characterization of the trypsin-III from Monterey sardine (Sardinops caeruleus): Insights on the cold-adaptation from the A236N mutant

Trypsins (E.C. 3.4.21.4) are digestive enzymes that catalyze the hydrolysis of peptide bonds containing arginine and lysine residues. Some trypsins from fish species are active at temperatures just above freezing, and for that are called cold-adapted enzymes, having many biotechnological applications. In this work, we characterized a recombinant trypsin-III from Monterey sardine (Sardinops caeruleus) and studied the role of a single residue on its cold-adapted features. The A236N mutant from sardine trypsin-III showed higher activation energy for the enzyme-catalyzed reaction, it was more active at higher temperatures, and exhibited a higher thermal stability than the wild-type enzyme, suggesting a key role of this residue. The thermodynamic activation parameters revealed an increase in the activation enthalpy for the A236N mutant, suggesting the existence of more intramolecular contacts during the activation step. Molecular models for both enzymes suggest that a hydrogen-bond involving N236 may contact the C-terminal α-helix to the vicinity of the active site, thus affecting the biochemical and thermodynamic properties of the enzyme.

Insight to the residue in P2 position prevents the peptide inhibitor from being hydrolyzed by serine proteases

Peptidic inhibitors of proteases are attracting increasing interest not only as drug candidates but also for studying the function and regulation mechanisms of these enzymes. Previously, we screened out a cyclic peptide inhibitor of human uPA [Formula: see text] and found that Ala substitution of P2 residue turns upain-1 to a substrate. To further investigate the effect of P2 residue on the peptide behavior transformation, we constructed upain-1-W3F, which has Phe replacement in the P2 position. We determined K_D and K_i of upain-1-W3F and found that upain-1-W3F might still exist as an inhibitor. Furthermore, the high-resolution crystal structure of upain-1-W3F·uPA reveals that upain-1-W3F indeed stays as an intact inhibitor bind to uPA. We thus propose that the P2 residue plays a nonnegligible role in the conversion of upain-1 to a substrate. These results also proposed a strategy to optimize the pharmacological properties of peptide-based drug candidates by hydrophobicity and steric hindrance.Abbreviations : uPA: urokinase-type plasminogen activator; SPD: serine protease domain; S1 pocket: specific substrate-binding pocket.

Target-Directed Self-Assembly of Homodimeric Drugs Against β-Tryptase

Tryptase, a serine protease released from mast cells, is implicated in many allergic and inflammatory disorders. Human tryptase is a donut-shaped tetramer with the active sites facing inward forming a central pore. Bivalent ligands spanning two active sites potently inhibit this configuration, but these large compounds have poor drug-like properties. To overcome some of these challenges, we developed self-assembling molecules, called coferons, which deliver a larger compound in two parts. Using a pharmacophoric core and reversibly binding linkers to span two active sites, we have successfully produced three novel homodimeric tryptase inhibitors. Upon binding to tryptase, compounds reassembled into flexible homodimers, with significant improvements in IC₅₀ (0.19 ± 0.08 μM) over controls (5.50 ± 0.09 μM), and demonstrate good activity in mast cell lines. These studies provide validation for this innovative technology that is especially well-suited for the delivery of dimeric drugs to modulate intracellular macromolecular targets.

Target-directed Dynamic Combinatorial Chemistry: A Study on Potentials and Pitfalls as Exemplified on a Bacterial Target

Target-directed dynamic combinatorial chemistry (DCC) is an emerging technique for the efficient identification of inhibitors of pharmacologically relevant targets. In this contribution, we present an application for a bacterial target, the lectin FimH, a crucial virulence factor of uropathogenic E. coli being the main cause of urinary tract infections. A small dynamic library of acylhydrazones was formed from aldehydes and hydrazides and equilibrated at neutral pH in presence of aniline as nucleophilic catalyst. The major success factors turned out to be an accordingly adjusted ratio of scaffolds and fragments, an adequate sample preparation prior to HPLC analysis, and the data processing. Only then did the ranking of the dynamic library constituents correlate well with affinity data. Furthermore, as a support of DCC applications especially to larger libraries, a new protocol for improved hit identification was established.

Synthesis of α-aminoboronic acids

This review describes available methods for the preparation of α-aminoboronic acids in their racemic or in their enantiopure form. Both, highly stereoselective syntheses and asymmetric procedures leading to the stereocontrolled generation of α-aminoboronic acid derivatives are included. The preparation of acyclic, carbocyclic and azacyclic α-aminoboronic acid derivatives is covered. Within each section, the different synthetic approaches have been classified according to the key bond which is formed to complete the α-aminoboronic acid skeleton.

Insights into the serine protease mechanism based on structural observations of the conversion of a peptidyl serine protease inhibitor to a substrate

BACKGROUND

Serine proteases are one of the most studied group of enzymes. Despite the extensive mechanistic studies, some crucial details remain controversial, for example, how the cleaved product is released in the catalysis reaction. A cyclic peptidyl inhibitor (CSWRGLENHRMC, upain-1) of a serine protease, urokinase-type plasminogen activator (uPA), was found to become a slow substrate and cleaved slowly upon the replacement of single residue (W3A).

METHODS

By taking advantage of the unique property of this peptide, we report the high-resolution structures of uPA in complex with upain-1-W3A peptide at four different pH values by X-ray crystallography.

RESULTS

In the structures obtained at low pH (pH4.6 and 5.5), the cyclic peptide upain-1-W3A was found to be intact and remained in the active site of uPA. At 7.4, the scissile bond of the peptide was found cleaved, showing that the peptide became a uPA substrate. At pH9.0, the C-terminal part of the substrate was no longer visible, and only the P1 residue occupying the S1 pocket was identified.

CONCLUSIONS

The analysis of these structures provides explanations why the upain-1-W3A is a slow substrate. In addition, we clearly identified the cleaved fragments of the peptide at both sides of the scissile bond in the active site of the enzyme, showing a slow release of the cleaved peptide.

GENERAL SIGNIFICANCE

This work indicates that the quick release of the cleaved P' fragment after the first step of hydrolysis may not always be needed for the second hydrolysis.

Dynamic combinatorial/covalent chemistry: a tool to read, generate and modulate the bioactivity of compounds and compound mixtures

Reversible covalent bond formation under thermodynamic control adds reactivity to self-assembled supramolecular systems, and is therefore an ideal tool to assess complexity of chemical and biological systems. Dynamic combinatorial/covalent chemistry (DCC) has been used to read structural information by selectively assembling receptors with the optimum molecular fit around a given template from a mixture of reversibly reacting building blocks. This technique allows access to efficient sensing devices and the generation of new biomolecules, such as small molecule receptor binders for drug discovery, but also larger biomimetic polymers and macromolecules with particular three-dimensional structural architectures. Adding a kinetic factor to a thermodynamically controlled equilibrium results in dynamic resolution and in self-sorting and self-replicating systems, all of which are of major importance in biological systems. Furthermore, the temporary modification of bioactive compounds by reversible combinatorial/covalent derivatisation allows control of their release and facilitates their transport across amphiphilic self-assembled systems such as artificial membranes or cell walls. The goal of this review is to give a conceptual overview of how the impact of DCC on supramolecular assemblies at different levels can allow us to understand, predict and modulate the complexity of biological systems.

Effect of ligand binding on the dynamics of trypsin. Comparison of different approaches

The intramolecular signal transduction induced by the binding of ligands to trypsin was investigated by molecular dynamics simulations. Ligand binding changes the residue-residue interaction energies and suppresses the mobility of loops that are in direct contact with the ligand. The reduced mobility of these loops results in the altered flexibility of the nearby loops and thereby transmits the information from ligand binding site to the remote sites. The analysis of the flexibility of all residues confirmed the coupling between loops L1 (185-188) and L2 (221-224) and the residues in the active center. The significance of S1 pocket residues for the signal transduction from the active center to the substrate-binding site was confirmed by the dynamical network and covariance matrix analyses. Gaussian network model and principal component analysis demonstrated that the active center residues had zero amplitude in the slowest fluctuations acting as hinges or anchors. Overall, our results provide a new insight into protein-ligand interactions and show how the allosteric signaling may occur.

Rational design of a transition state analogue with picomolar affinity for Pseudomonas aeruginosa PvdQ, a siderophore biosynthetic enzyme

The Pseudomonas aeruginosa enzyme PvdQ can process different substrates involved in quorum-sensing or in siderophore biosynthesis. Substrate selectivity was evaluated using steady-state kinetic constants for hydrolysis of N-acyl-homoserine lactones (HSLs) and p-nitrophenyl fatty acid esters. PvdQ prefers substrates with alkyl chains between 12 and 14 carbons long that do not bear a 3-oxo substitution and is revealed here to have a relatively high specificity constant for selected N-acyl-HSLs (kcat/KM = 10(5) to 10(6) M(-1) s(-1)). However, endogenous P. aeruginosa N-acyl-HSLs are ≥100-fold disfavored, supporting the conclusion that PvdQ was not primarily evolved to regulate endogenous quorum-sensing. PvdQ plays an essential biosynthetic role for the siderophore pyoverdine, on which P. aeruginosa depends for growth in iron-limited environments. A series of alkylboronate inhibitors was found to be reversible, competitive, and extremely potent (Ki ≥ 190 pM). A 1.8 Å X-ray structure shows that 1-tridecylboronic acid forms a monocovalent bond with the N-terminal β-chain Ser residue in the PvdQ heterodimer, mimicking a reaction transition state. This boronic acid inhibits growth of P. aeruginosa in iron-limited media, reproducing the phenotype of a genetic pvdQ disruption, although co-administration of an efflux pump inhibitor is required to maintain growth inhibition. These findings support the strategy of designing boron-based inhibitors of siderophore biosynthetic enzymes to control P. aeruginosa infections.

Screening and in situ synthesis using crystals of a NAD kinase lead to a potent antistaphylococcal compound

Making new ligands for a given protein by in situ ligation of building blocks (or fragments) is an attractive method. However, it suffers from inherent limitations, such as the limited number of available chemical reactions and the low information content of usual chemical library deconvolution. Here, we describe a focused screening of adenosine derivatives using X-ray crystallography. We discovered an unexpected and biocompatible chemical reactivity and have simultaneously identified the mode of binding of the resulting products. We observed that the NAD kinase from Listeria monocytogenes (LmNADK1) can promote amide formation between 5'-amino-5'-deoxyadenosine and carboxylic acid groups. This unexpected reactivity allowed us to bridge in situ two adenosine derivatives to fully occupy the active NAD site. This guided the design of a close analog showing micromolar inhibition of two human pathogenic NAD kinases and potent bactericidal activity against Staphylococcus aureus in vitro.

Design, synthesis and biological study of pinacolyl boronate-substituted stilbenes as novel lipogenic inhibitors

A pilot library of novel 4,4,5,5-tetramethyl-2-(4-substitutedstyrylphenyl)-1,3,2 dioxaborolane derivatives has been synthesized. 4-(4,4,5,5-Tetramethyl-1,3,2-dioxaboratophenyl)-methyl triphenylphosphonium bromide 3 was treated with various aldehydes in the presence of 3 equiv of (t)BuONa in DMF, and stirred at room temperature for 4-6h to yield the corresponding boron-containing stilbene derivatives in 71-94% yields. Several of them, including BF102 and BF175, have the lipogenesis inhibitory effect by suppressing lipogenic gene expression in mammalian hepatocytes. Further, BF102 also inhibits cholesterol biosynthesis by suppressing HMG-CoA reductase gene expression in hepatocytes. Interestingly, our preliminary in vivo data suggests that BF102 has no significant toxicity in mice at the highest possible dose we can administered. Thus, BF102 is a potential lead for the next generation of lipid-lowering drugs.

Rosetta in CAPRI rounds 13-19

Modeling the conformational changes that occur on binding of macromolecules is an unsolved challenge. In previous rounds of the Critical Assessment of PRediction of Interactions (CAPRI), it was demonstrated that the Rosetta approach to macromolecular modeling could capture side chain conformational changes on binding with high accuracy. In rounds 13-19 we tested the ability of various backbone remodeling strategies to capture the main-chain conformational changes observed during binding events. These approaches span a wide range of backbone motions, from limited refinement of loops to relieve clashes in homologous docking, through extensive remodeling of loop segments, to large-scale remodeling of RNA. Although the results are encouraging, major improvements in sampling and energy evaluation are clearly required for consistent high accuracy modeling. Analysis of our failures in the CAPRI challenges suggest that conformational sampling at the termini of exposed beta strands is a particularly pressing area for improvement.

MDockPP: A hierarchical approach for protein-protein docking and its application to CAPRI rounds 15-19

A hierarchical approach has been developed for protein-protein docking. In the first step, a Fast Fourier Transform (FFT)-based docking algorithm is used to globally sample all putative binding modes, in which the protein is represented by a reduced model, that is, each side chain on the protein surface is represented by its center of mass. Compared to conventional FFT docking with all-atom models, the FFT docking method with a reduced model is expected to generate more hits because it allows larger side-chain flexibility. Next, the filtered binding modes (normally several thousands) are refined by an iteratively derived knowledge-based scoring function ITScorePP and by considering backbone/loop flexibility using an ensemble docking algorithm. The distance-dependent potentials of ITScorePP were extracted by a physics-based iterative method, which circumvents the long-standing reference state problem in the knowledge-based approaches. With this hierarchical protocol, we have participated in the CAPRI experiments for Rounds 15-19 of 11 targets (T32-T42). In the predictor experiments, we achieved correct binding modes for six targets: three are with high accuracy (T40 for both distinct binding modes, T41, and T42), two are with medium accuracy (T34 and T37), and one is acceptable (T32). In the scorer experiments, of the seven target complexes that contain at least one acceptable mode submitted by the CAPRI predictor groups, we obtained correct binding modes for four targets: three are with high accuracy (T37, T40, and T41) and one is with medium accuracy (T34), suggesting good accuracy and robustness of ITScorePP.

CAPRI targets T29-T42: proving ground for new docking procedures

The critical assessment of protein interactions (CAPRI) experiment provides a unique opportunity for unbiased assessment of docking procedures. The recent CAPRI targets T29-T42 entailed docking of bound, unbound, and modeled structures, presenting a wide range of prediction difficulty. We submitted accurate predictions for targets T40, T41, and T42, a good prediction for T32 and acceptable predictions for T29 and T34. The accuracy of our docking results generally matched the prediction difficulty; hence, docking of modeled proteins produced less accurate results. However, there were interesting exceptions: an accurate prediction was submitted for the dimer of modeled tetratricopeptide repeat (T42) and only an acceptable prediction for the bound/unbound case T29. The ensembles of docking models produced in the scans included an acceptable or better prediction for every target. We show here that our recently developed postscan reevaluation procedure, which tests propensity and solvation measures of the whole interface and the interface core, successfully distinguished these predictions from false docking models. For enzyme-inhibitor targets, we show that the distance of the interface from the enzyme's centroid ranked high native like docking models. Also, for one case we demonstrate that docking of an ensemble of conformers produced by normal modes analysis can improve the accuracy of the prediction.

Strengths and weaknesses of data-driven docking in critical assessment of prediction of interactions

The recent CAPRI rounds have introduced new docking challenges in the form of protein-RNA complexes, multiple alternative interfaces, and an unprecedented number of targets for which homology modeling was required. We present here the performance of HADDOCK and its web server in the CAPRI experiment and discuss the strengths and weaknesses of data-driven docking. HADDOCK was successful for 6 out of 9 complexes (6 out of 11 targets) and accurately predicted the individual interfaces for two more complexes. The HADDOCK server, which is the first allowing the simultaneous docking of generic multi-body complexes, was successful in 4 out of 7 complexes for which it participated. In the scoring experiment, we predicted the highest number of targets of any group. The main weakness of data-driven docking revealed from these last CAPRI results is its vulnerability for incorrect experimental data related to the interface or the stoichiometry of the complex. At the same time, the use of experimental and/or predicted information is also the strength of our approach as evidenced for those targets for which accurate experimental information was available (e.g., the 10 three-stars predictions for T40!). Even when the models show a wrong orientation, the individual interfaces are generally well predicted with an average coverage of 60% ± 26% over all targets. This makes data-driven docking particularly valuable in a biological context to guide experimental studies like, for example, targeted mutagenesis.

A generalized approach to sampling backbone conformations with RosettaDock for CAPRI rounds 13-19

In CAPRI rounds 13-19, the most native-like structure predicted by RosettaDock resulted in two high, one medium, and one acceptable accuracy model out of 13 targets. The current rounds of CAPRI were especially challenging with many unbound and homology modeled starting structures. Novel docking methods, including EnsembleDock and SnugDock, allowed backbone conformational sampling during docking and enabled the creation of more accurate models. For Target 32, α-amylase/subtilisin inhibitor-subtilisin savinase, we sampled different backbone conformations at an interfacial loop to produce five high-quality models including the most accurate structure submitted in the challenge (2.1 Å ligand rmsd, 0.52 Å interface rmsd). For Target 41, colicin-immunity protein, we used EnsembleDock to sample the ensemble of nuclear magnetic resonance (NMR) models of the immunity protein to generate a medium accuracy structure. Experimental data identifying the catalytic residues at the binding interface for Target 40 (trypsin-inhibitor) were used to filter RosettaDock global rigid body docking decoys to determine high accuracy predictions for the two distinct binding sites in which the inhibitor interacts with trypsin. We discuss our generalized approach to selecting appropriate methods for different types of docking problems. The current toolset provides some robustness to errors in homology models, but significant challenges remain in accommodating larger backbone uncertainties and in sampling adequately for global searches.

SwarmDock and the use of normal modes in protein-protein docking

Here is presented an investigation of the use of normal modes in protein-protein docking, both in theory and in practice. Upper limits of the ability of normal modes to capture the unbound to bound conformational change are calculated on a large test set, with particular focus on the binding interface, the subset of residues from which the binding energy is calculated. Further, the SwarmDock algorithm is presented, to demonstrate that the modelling of conformational change as a linear combination of normal modes is an effective method of modelling flexibility in protein-protein docking.

Trypsin-ligand binding free energy calculation with AMOEBA

The binding free energies of several benzamidine -like inhibitors to trypsin were examined using a polarizable molecular mechanics potential. All the computed binding free energies are in good agreement with the experimental data. From free energy decomposition, electrostatic interaction was indicated to be the driving force for the binding. MD simulations show that the ligands form hydrogen bonds with trypsin and water molecules nearby in a competitive fashion. While the binding free energy is independent of the ligand dipole moment, it shows a strong correlation with the ligand molecular polarizability.

Current status of the AMOEBA polarizable force field

Molecular force fields have been approaching a generational transition over the past several years, moving away from well-established and well-tuned, but intrinsically limited, fixed point charge models toward more intricate and expensive polarizable models that should allow more accurate description of molecular properties. The recently introduced AMOEBA force field is a leading publicly available example of this next generation of theoretical model, but to date, it has only received relatively limited validation, which we address here. We show that the AMOEBA force field is in fact a significant improvement over fixed charge models for small molecule structural and thermodynamic observables in particular, although further fine-tuning is necessary to describe solvation free energies of drug-like small molecules, dynamical properties away from ambient conditions, and possible improvements in aromatic interactions. State of the art electronic structure calculations reveal generally very good agreement with AMOEBA for demanding problems such as relative conformational energies of the alanine tetrapeptide and isomers of water sulfate complexes. AMOEBA is shown to be especially successful on protein-ligand binding and computational X-ray crystallography where polarization and accurate electrostatics are critical.

Trypsin-ligand binding free energies from explicit and implicit solvent simulations with polarizable potential

We have calculated the binding free energies of a series of benzamidine-like inhibitors to trypsin with a polarizable force field using both explicit and implicit solvent approaches. Free energy perturbation has been performed for the ligands in bulk water and in protein complex with molecular dynamics simulations. The binding free energies calculated from explicit solvent simulations are well within the accuracy of experimental measurement and the direction of change is predicted correctly in all cases. We analyzed the molecular dipole moments of the ligands in gas, water and protein environments. Neither binding affinity nor ligand solvation free energy in bulk water shows much dependence on the molecular dipole moments of the ligands. Substitution of the aromatic or the charged group in the ligand results in considerable change in the solvation energy in bulk water and protein whereas the binding affinity varies insignificantly due to cancellation. The effect of chemical modification on ligand charge distribution is mostly local. Replacing benzene with diazine has minimal impact on the atomic multipoles at the amidinium group. We have also utilized an implicit solvent based end-state approach to evaluate the binding free energies of these inhibitors. In this approach, the polarizable multipole model combined with Poisson-Boltzmann/surface area (PMPB/SA) provides the electrostatic interaction energy and the polar solvation free energy. Overall the relative binding free energies obtained from the MM-PMPB/SA model are in good agreement with the experimental data.

Structure of a serine protease poised to resynthesize a peptide bond

The serine proteases are among the most thoroughly studied enzymes, and numerous crystal structures representing the enzyme-substrate complex and intermediates in the hydrolysis reactions have been reported. Some aspects of the catalytic mechanism remain controversial, however, especially the role of conformational changes in the reaction. We describe here a high-resolution (1.46 A) crystal structure of a complex formed between a cleaved form of bovine pancreatic trypsin inhibitor (BPTI) and a catalytically inactive trypsin variant with the BPTI cleavage site ideally positioned in the active site for resynthesis of the peptide bond. This structure defines the positions of the newly generated amino and carboxyl groups following the 2 steps in the hydrolytic reaction. Comparison of this structure with those representing other intermediates in the reaction demonstrates that the residues of the catalytic triad are positioned to promote each step of both the forward and reverse reaction with remarkably little motion and with conservation of hydrogen-bonding interactions. The results also provide insights into the mechanism by which inhibitors like BPTI normally resist hydrolysis when bound to their target proteases.

Protein/ligand binding free energies calculated with quantum mechanics/molecular mechanics

The calculation of binding affinities for flexible ligands has hitherto required the availability of reliable molecular mechanics parameters for the ligands, a restriction that can in principle be lifted by using a mixed quantum mechanics/molecular mechanics (QM/MM) representation in which the ligand is treated quantum mechanically. The feasibility of this approach is evaluated here, combining QM/MM with the Poisson-Boltzmann/surface area model of continuum solvation and testing the method on a set of 47 benzamidine derivatives binding to trypsin. The experimental range of the absolute binding energy (DeltaG = -3.9 to -7.6 kcal/mol) is reproduced well, with a root-mean-square (RMS) error of 1.2 kcal/mol. When QM/MM is applied without reoptimization to the very different ligands of FK506 binding protein the RMS error is only 0.7 kcal/mol. The results show that QM/MM is a promising new avenue for automated docking and scoring of flexible ligands. Suggestions are made for further improvements in accuracy.

Inhibition of chymotrypsin by a complex of ortho-vanadate and benzohydroxamic acid: structure of the inert complex and its mechanistic interpretation

Serine proteases, like serine beta-lactamases, are rapidly and covalently inhibited by suitably designed phosph(on)ates. The active sites of these enzymes must, therefore, be able to stabilize the pentacoordinated transition states of phosphyl transfer reactions as well as the tetrahedral transition states of acyl transfers. It follows that these enzymes should also be inhibited by molecules capable of generating inert pentacoordinated species. We (J.H.B. and R.F.P.) have previously shown that these enzymes are, in fact, rapidly and reversibly inhibited by 1:1 complexes of vanadate and hydroxamic acids. In this paper, we present the first crystal structure of an acyl transferase inhibited by vanadate. The complex of vanadate and benzohydroxamic acid is a competitive inhibitor of alpha-chymotrypsin with a KI value of 16 muM. In the structure, obtained at a resolution of 1.5 A, the protein is conformationally little different from the apoenzyme. The vanadium, in a distorted octahedral ligand field, is covalently bound to the active site serine oxygen group. One oxgen ligand, presumably anionic, is located in the oxyanion hole. Another is directed roughly in the direction of the acyl transfer leaving group, and a third in the direction of the S2 site. The hydroxamate is bound to vanadium through the hydroxyl oxygen and also, more weakly, through the carbonyl group, to form a five-membered chelate ring. The effect of this chelation is to place the phenyl group of the inhibitor into the important S1 specificity site. The hydroxamate oxygen is directed in line away from the Ser195 Ogamma, approximating the direction of departure of a leaving group in phosphyl transfer. The entire complex can be seen as a reasonable mimic of a phosphyl transfer transition state where the leaving group is extended into the S1 site.

Synthesis and biological evaluation of boronic acid containing cis-stilbenes as apoptotic tubulin polymerization inhibitors

A series of boronic acid containing cis-stilbenes as potent inhibitors of tubulin polymerization was synthesized by the introduction of boronic acid as an acceptor-type functional group into the aromatic ring B of the combretastatin framework. High cell-growth inhibition was observed with boron compounds 13 c and 13 d, in which a hydroxy group on the aromatic ring B of combretastatin A-4 was replaced with boronic acid; IC50 values toward B-16 and 1-87 cell lines are 0.48-2.1 microM. Compounds 13 c and 13 d exhibited significant inhibitory activity toward tubulin polymerization (IC50=21-22 microM). The carboxylic acid derivative 17, which can be considered as a mimic of boronic acid 13 c, did not show significant inhibition of cell growth or tubulin polymerization. According to the FACScan analysis using Jurkat cells, apoptosis was induced after incubation for 8 h with 13 c at a concentration of >10(-8) M. Growth inhibitory experiments against a panel of 39 human cancer cell lines revealed 13 c to inhibit growth differently than combretastatin A-4; the correlation coefficient (r) between the two compounds was 0.553 in the COMPARE analysis.

Utilization of Self-Sorting Processes To Generate Dynamic Combinatorial Libraries with New Network Topologies

The synthesis of water-soluble, organometallic macrocycles is described. They were obtained by self-assembly in reactions of the half-sandwich complexes [[Ru(C6H5Me)Cl2]2], [[Ru(p-cymene)Cl2]2], [[Rh(Cp)Cl2]2], and [[Ir(Cp*)Cl2]2] with the ligand 5-dimethylaminomethyl-3-hydroxy-2-methyl-4-(1H)-pyridone in buffered aqueous solution at pH 8. The structure of the Ru-(p-cymene) complex was determined by single-crystal X-ray crystallography. Upon mixing, these complexes undergo scrambling reactions to give dynamic combinatorial libraries. In combination with structurally related complexes based on amino-methylated 3-hydroxy-2-(1H)-pyridone ligands, an exchange of metal fragments but no mixing of ligands was observed. This self-sorting behavior was used to construct dynamic combinatorial libraries of macrocycles, in which two four-component sub-libraries are connected by two common building blocks. This type of network topology influences the adaptive behavior of the library as demonstrated in selection experiments with lithium ions as the target.

Rigidity and Flexibility of Dipeptidyl Peptidase IV: Crystal Structures of and Docking Experiments with DPIV

Dipeptidyl peptidase IV (DPIV) is an alpha,beta-hydrolase-like serine exopeptidase, which removes dipeptides, preferentially with a C-terminal l-Pro residue, from the N terminus of longer peptide substrates. Previously, we determined the tetrameric 1.8A crystal structure of native porcine DPIV. Each monomer is composed of a beta-propeller and a catalytic domain, which together embrace an internal cavity housing the active centre. This cavity is connected to the bulk solvent by a "propeller opening" and a "side opening". Here, we analyse DPIV complexes with a t-butyl-Gly-Pro-Ile tripeptide, Pro-boroPro, a piperazine purine compound, and aminoethyl phenyl sulfonylfluoride. The latter two compounds bind to the active-site groove in a compact and a quite bulky manner, respectively, causing considerable shifts of the catalytic Ser630 side-chain and of the Tyr547 phenolic group, which forms the oxyanion hole. The tripeptide, mimicking a peptide substrate, is clamped to the active site through tight interactions via its N-terminal alpha-ammonium group, the P2 carbonyl group, the P1-l-Pro side-chain, the C-terminal carboxylate group, and the stable orthoacid ester amide formed between the scissile peptide carbonyl group and Ser630 O(gamma). This stable trapping of the tripeptide could be due to stabilization of the protonated His740 imidazolium cation by the adjacent negatively charged C-terminal carboxylate group, preventing proton transfer to the leaving group nitrogen atom. Docking experiments with the compact rigid 58 residue protein aprotinin, which had been shown to be processed by DPIV, indicate that the Arg1-Pro2 N terminus can access the DPIV active site only upon widening of its side openings, probably by separation of the first and the last propeller blades, and/or of the catalytic and the propeller domain.

Design, synthesis, and biological evaluation of aminoboronic acids as growth-factor receptor inhibitors of EGFR and VEGFR-1 tyrosine kinases

A series of aminoboronic acids was synthesized based on the structure of lavendustin pharmacophore 1. Their inhibitory activities against the epidermal growth-factor receptor (EGFR) and vascular endothelial growth-factor receptor-1 (VEGFR-1, Flt-1) protein tyrosine kinases, and various protein kinases, PKA, PKC, PTK, and eEF2K were evaluated. Selective inhibition activities were observed in a series of aminoboronic acids. 4-Methoxy-3-((2- methoxyphenylamino)methyl)phenylboronic acid 10 inhibited EGFR tyrosine kinase, whereas 4-(2,5-dihydroxybenzylamino)phenylboronic acid 12 inhibited Flt-1 protein kinase, although lavendustin pharmacophore 1 inhibited both EGFR and Flt-1 kinases at a compound concentration of 1.0 microg mL(-1). The selective inhibition of EGFR by 10 is considered to be due to the substitution of the dihydroxy groups on the benzyl moiety for a boronic acid group at the para position, whereas the selective inhibition of Flt-1 by 12 is due to the substitution of the carboxyl group on the aniline moiety in the lavendustin pharmacophore 1 for a boronic acid group.

Synthesis and evaluation of macrocyclic transition state analogue inhibitors for alpha-chymotrypsin

Lactams 1 and 2 are readily formed from acyclic precursors in the presence of trypsin and chymotrypsin, respectively, identifying the macrocyclic ring system as a potential motif for constrained transition state analogue inhibitors of the serine peptidases. Ketone 3 was synthesized and shown to be a modest inhibitor of chymotrypsin (Ki = 220 microM), albeit 4-fold more potent than the acyclic hydroxy acid 25 (Ki = 1.5 mM as a mixture of epimers). A precursor (31) to the amino boronic acid 4 was also prepared; although this derivative was a potent inhibitor of chymotrypsin (Ki = 130 nM) by virtue of the boronic acid moiety, it showed no advantage over the des-amino analogue 32 (Ki = 120 nM), which is not capable of cyclizing.

Method for computing protein binding affinity

A Monte Carlo method is given to compute the binding affinity of a ligand to a protein. The method involves extending configuration space by a discrete variable indicating whether the ligand is bound to the protein and a special Monte Carlo move, which allows transitions between the unbound and bound states. Provided that an accurate protein structure is given, that the protein-ligand binding site is known, and that an accurate chemical force field together with a continuum solvation model is used, this method provides a quantitative estimate of the free energy of binding.

On the transferability of hydration-parametrized continuum electrostatics models to solvated binding calculations

Using molecular mechanics force field partial atomic charges, we show the nonuniqueness of the parametrization of continuum electrostatics models with respect to solute atomic radii and interior dielectric constant based on hydration (vacuum-to-water transfer) free energy data available for small molecules. Moreover, parameter sets that are optimal and equivalent for hydration free energy calculations lead to large variations of calculated absolute and relative electrostatic binding free energies. Hence, parametrization of solvation effects based on hydration data, although a necessary condition, is not sufficient to guarantee its transferability to the calculation of binding free energies in solution.

Profiling charge complementarity and selectivity for binding at the protein surface

A novel analysis and representation of the protein surface in terms of electrostatic binding complementarity and selectivity is presented. The charge optimization methodology is applied in a probe-based approach that simulates the binding process to the target protein. The molecular surface is color coded according to calculated optimal charge or according to charge selectivity, i.e., the binding cost of deviating from the optimal charge. The optimal charge profile depends on both the protein shape and charge distribution whereas the charge selectivity profile depends only on protein shape. High selectivity is concentrated in well-shaped concave pockets, whereas solvent-exposed convex regions are not charge selective. This suggests the synergy of charge and shape selectivity hot spots toward molecular selection and recognition, as well as the asymmetry of charge selectivity at the binding interface of biomolecular systems. The charge complementarity and selectivity profiles map relevant electrostatic properties in a readily interpretable way and encode information that is quite different from that visualized in the standard electrostatic potential map of unbound proteins.

Optimizing the reversibility of hydrazone formation for dynamic combinatorial chemistry

Hydrazones from hydrazines bearing electron withdrawing groups, and aromatic or aliphatic aldehydes form and hydrolyse rapidly in water at neutral pH.

Crystal structure and nucleotide sequence of an anionic trypsin from chum salmon (Oncorhynchus keta) in comparison with Atlantic salmon (Salmo salar) and bovine trypsin

The nucleotide sequence and crystal structure of chum salmon trypsin (CST) are now reported. The cDNA isolated from the pyloric caeca of chum salmon encodes 222 amino acid residues, the same number of residues as the anionic Atlantic salmon trypsin (AST), but one residue less than bovine beta-trypsin (BT). The net charge on CST determined from the sum of all charged amino acid side-chains is -3. There are 79 sequence differences between CST and BT, but only seven sequence differences between CST and AST. Anionic CST isolated from pyloric caeca has also been purified and crystallized; the structure of the CST-benzamidine complex has been determined to 1.8A resolution. The overall tertiary structure of CST is similar to that of AST and BT, but some differences are observed among the three trypsins. The most striking difference is at the C terminus of CST, where the expected last two residues are absent. The absence of these residues likely increases the flexibility of CST by the loss of important interactions between the N and C-terminal domains. Similarly, the lack of Tyr151 in CST (when compared with BT) allows more space for Gln192 in the active site thereby increasing substrate accessibility to the binding pocket. Lys152 in CST also adopts the important role of stabilizing the loop from residue 142 to 153. These observations on CST provide a complementary view of a second cold-adapted trypsin, which in comparison with the structures of AST and BT, suggest a structural basis for differences in enzymatic activity between enzymes from cold-adapted species and mammals.

Can the calculation of ligand binding free energies be improved with continuum solvent electrostatics and an ideal-gas entropy correction?

The prediction of a ligand binding constant requires generating three-dimensional structures of the complex concerned and reliably scoring these structures. Here, the scoring problem is investigated by examining benzamidine-like inhibitors of trypsin, a system for which errors in the structures are small. Precise and consistent binding free energies for the inhibitors are determined experimentally for this test system. To examine possible improvement of scoring methods, we test the suitability of continuum electrostatics to account for solvation effects and use an ideal-gas entropy correction to account for the changes in the degrees of freedom of the ligand. The small observed root-mean-square deviation of 0.55 kcal/mol of the calculated relative to the experimental values indicates that the essentials of the binding process have been captured. Even though all six ligands make the same salt bridge and H-bonds to the protein, the electrostatic contribution varies among the ligands by as much as 2 kcal/mol. Moreover, although the ligands are rigid and similar in size, the entropic terms also significantly affect the relative binding affinities (by up to 2.7 kcal/mol). The present approach to solvation and entropy may allow the ranking of the ligands to be considerably improved at a cost that makes the method applicable to the optimization of lead compounds or to the screening of small collections of ligands.