Cyclic (hydroxyphosphinyl)acyl dipeptides: a new class of angiotensin converting enzyme inhibitors

Conformationally constrained phenylbutyl(hydroxyphosphinyl)acyl dipeptides are potent inhibitors of angiotensin converting enzyme. The activity enhancement obtained by introducing conformational constraint into these molecules is greater than for related sulfhydryl and carboxyl analogs. The results are interpreted in terms of a binding model which optimally positions both zinc binding and hydrophobic groups for active site binding.     

Think Twice: Understanding the High Potency of Bis(phenyl)methane Inhibitors of Thrombin

Successful design of potent and selective protein inhibitors, in terms of structure-based drug design, strongly relies on the correct understanding of the molecular features determining the ligand binding to the target protein. We present a case study of serine protease inhibitors with a bis(phenyl)methane moiety binding into the S3 pocket. These inhibitors bind with remarkable potency to the active site of thrombin, the blood coagulation factor IIa. A combination of X-ray crystallography and isothermal titration calorimetry provides conclusive insights into the driving forces responsible for the surprisingly high potency of these inhibitors. Analysis of six well-resolved crystal structures (resolution 1.58-2.25 Å) along with the thermodynamic data allows an explanation of the tight binding of the bis(phenyl)methane inhibitors. Interestingly, the two phenyl rings contribute to binding affinity for very different reasons — a fact that can only be elucidated by a structure-based approach. The first phenyl moiety occupies the hydrophobic S3 pocket, resulting in a mainly entropic advantage of binding. This observation is based on the displacement of structural water molecules from the S3 pocket that are observed in complexes with inhibitors that do not bind in the S3 pocket. The same classic hydrophobic effect cannot explain the enhanced binding affinity resulting from the attachment of the second, more solvent-exposed phenyl ring. For the bis(phenyl)methane inhibitors, an observed adaptive rotation of a glutamate residue adjacent to the S3 binding pocket attracted our attention. The rotation of this glutamate into salt-bridging distance with a lysine moiety correlates with an enhanced enthalpic contribution to binding for these highly potent thrombin binders. This explanation for the magnitude of the attractive force is confirmed by data retrieved by a Relibase search of several thrombin-inhibitor complexes deposited in the Protein Data Bank exhibiting similar molecular features.Special attention was attributed to putative changes in the protonation states of the interaction partners. For this purpose, two analogous inhibitors differing mainly in their potential to change the protonation state of a hydrogen-bond donor functionality were compared. Buffer dependencies of the binding enthalpy associated with complex formation could be traced by isothermal titration calorimetry, which revealed, along with analysis of the crystal structures (resolution 1.60 and 1.75 Å), that a virtually compensating proton interchange between enzyme, inhibitor and buffer is responsible for the observed buffer-independent thermodynamic signatures.     

New chloro,fluorobenzylindole derivatives as integrase strand-transfer inhibitors (INSTIs) and their mode of action

The life cycle of HIV-1 requires extensive assistance from the integrase (IN) enzyme which therefore constitutes an attractive therapeutic target for the development of anti-AIDS agents. We herein report the synthesis and biological evaluation of new HIV integrase strand-transfer inhibitors (INSTIs) which proved to be also potent anti-HIV agents. The binding mode of the most representative molecules were also studied by induced-fit docking (IFD). The obtained IFD results were consistent with the mechanism of action proposed for this class of IN inhibitors, that is metal chelating/binding agents.     

Boron-based phosphodiesterase inhibitors show novel binding of boron to PDE4 bimetal center

We have used boron-based molecules to create novel, competitive, reversible inhibitors of phosphodiesterase 4 (PDE4). The co-crystal structure reveals a binding configuration which is unique compared to classical catechol PDE4 inhibitors, with boron binding to the activated water in the bimetal center. These phenoxybenzoxaboroles can be optimized to generate submicromolar potency enzyme inhibitors, which inhibit TNF-α, IL-2, IFN-γ, IL-5 and IL-10 activities in vitro and show safety and efficacy for topical treatment of human psoriasis. They provide a valuable new route for creating novel potent anti-PDE4 inhibitors.     

Identification of Novel Human Dipeptidyl Peptidase-IV Inhibitors of Natural Origin (Part I): Virtual Screening and Activity Assays

Background

There has been great interest in determining whether natural products show biological activity toward protein targets of pharmacological relevance. One target of particular interest is DPP-IV whose most important substrates are incretins that, among other beneficial effects, stimulates insulin biosynthesis and secretion. Incretins have very short half-lives because of their rapid degradation by DPP-IV and, therefore, inhibiting this enzyme improves glucose homeostasis. As a result, DPP-IV inhibitors are of considerable interest to the pharmaceutical industry. The main goals of this study were (a) to develop a virtual screening process to identify potential DPP-IV inhibitors of natural origin; (b) to evaluate the reliability of our virtual-screening protocol by experimentally testing the in vitro activity of selected natural-product hits; and (c) to use the most active hit for predicting derivatives with higher binding affinities for the DPP-IV binding site.

Methodology/Principal Findings

We predicted that 446 out of the 89,165 molecules present in the natural products subset of the ZINC database would inhibit DPP-IV with good ADMET properties. Notably, when these 446 molecules were merged with 2,342 known DPP-IV inhibitors and the resulting set was classified into 50 clusters according to chemical similarity, there were 12 clusters that contained only natural products for which no DPP-IV inhibitory activity has been previously reported. Nine molecules from 7 of these 12 clusters were then selected for in vitro activity testing and 7 out of the 9 molecules were shown to inhibit DPP-IV (where the remaining two molecules could not be solubilized, preventing the evaluation of their DPP-IV inhibitory activity). Then, the hit with the highest activity was used as a lead compound in the prediction of more potent derivatives.

Conclusions/Significance

We have demonstrated that our virtual-screening protocol was successful in identifying novel lead compounds for developing more potent DPP-IV inhibitors.     

Structural and dynamics analysis of matrix metalloproteinases MMP-2 complexed with chemically modified tetracyclines (CMTs)

Matrix metalloproteinases (MMPs) play a critical role in physiological processes and pathological conditions such tumor invasion and metastasis. In recent years, a number of MMP inhibitors have been proposed, including the chemically modified tetracyclines (CMTs), which have been evaluated in preclinical cancer models showing promising results. This work provides insights into the structure and dynamics of the MMP-2 catalytic domain complexed with seven CMT (CMT-n), based on the analysis of molecular dynamics trajectories in solution. The comparative analysis of various relevant molecular aspects of the different complexes of MMP-2 and CMT-n derivatives was performed aiming to elucidate the effect of ligands on the enzyme structure. These include the radial distribution function of the water molecules around the catalytic zinc, the solvent accessible surface area for the inhibitors and the root-mean-square fluctuation for all amino acid residues. The results help to understand the differences in the binding modes of related compounds and, therefore, add to further design of novel tetracycline-based inhibitors for MMP enzymes.     

Molecular docking study and development of an empirical binding free energy model for phosphodiesterase 4 inhibitors

In the present work, several computational methodologies were combined to develop a model for the prediction of PDE4B inhibitors' activity. The adequacy of applying the ligand docking approach, keeping the enzyme rigid, to the study of a series of PDE4 inhibitors was confirmed by a previous molecular dynamics analysis of the complete enzyme. An exhaustive docking procedure was performed to identify the most probable binding modes of the ligands to the enzyme, including the active site metal ions and the surrounding structural water molecules. The enzyme-inhibitor interaction enthalpies, refined by using the semiempirical molecular orbital approach, were combined with calculated solvation free energies and entropy considerations in an empirical free energy model that enabled the calculation of binding free energies that correlated very well with experimentally derived binding free energies. Our results indicate that both the inclusion of the structural water molecules close to the ions in the binding site and the use of a free energy model with a quadratic dependency on the ligand free energy of solvation are important aspects to be considered for molecular docking investigations involving the PDE4 enzyme family.     

Structures of complexes of rhizopuspepsin with pepstatin and other statine-containing inhibitors.

The three-dimensional structures of the complexes of the aspartic proteinase from Rhizopus chinensis (Rhizopuspepsin, EC 3.4.23.6) with pepstatin and two pepstatin-like peptide inhibitors of renin have been determined by X-ray diffraction methods and refined by restrained least-squares procedures. The inhibitors adopt an extended conformation and lie in the deep groove located between the two domains of the enzyme. Inhibitor binding is accompanied by a conformational change at the "flap," a beta-hairpin loop region, that projects over the binding cleft and closes down over the inhibitor, excluding water molecules from the vicinity of the scissile bond. The hydroxyl group of the central statyl residue of the inhibitors replaces the water molecule found between the two active aspartates, Asp-35 and Asp-218, in the native structure. The refined structures provide additional data to define the specific subsites of the enzyme and also show a system of hydrogen bonding to the inhibitor backbone similar to that observed for a reduced inhibitor.     

Synthesis and evaluation of unsaturated caprolactams as interleukin-1beta converting enzyme (ICE) inhibitors

Peptidomimetic compounds possessing a caprolactam ring constraint were prepared and evaluated as interleukin-1beta converting enzyme (ICE) inhibitors. The caprolactam ring was used to constrain the P3 region of our inhibitors. This strategy proved to be effective for the synthesis of ICE inhibitors, maintaining key hydrogen bond interactions with the enzyme and invoking a preferred conformation for binding. Several compounds exhibited IC(50) values less than 10nM in a caspase-1 enzyme assay and less than 100nM in a THP-1 whole cell assay measuring IL-1beta production. Two compounds, 13c and 13j, were found to have good oral bioavailability (>50%) in rats when administered as prodrugs.     

Establishing the principles of recognition in the adenine-binding region of an aminoglycoside antibiotic kinase [APH(3')-IIIa

The protein-based molecular recognition of the adenine ring has implications throughout biological systems. In this paper, we discuss the adenine-binding region of an aminoglycoside antibiotic kinase [APH(3')-IIIa], which serves as an excellent model system for proteins that bind the adenine ring. This enzyme employs a hydrogen-bonding network involving water molecules along with enzyme backbone/side-chain atoms and a pi-pi stacking interaction to recognize the adenine ring. Our approach utilized site-directed mutagenesis, adenosine analogues and a variety of biophysical methods to probe the contacts in the adenine-binding region of APH(3')-IIIa. The results point to the polar nature of an adenine-Met90 contact in this binding pocket and the important role that Met90, the "gatekeeper" residue in structurally similar Ser/Thr protein kinases, plays in adenine binding. The results also suggest that small changes in the structure of the adenine ring can lead to significant changes in the ability of these analogues to occupy the adenine-binding region of the enzyme. Additional computational experiments indicate that both size and electronic factors are important in the binding of aromatic systems in this interaction-rich pocket. The principles governing adenine recognition established in this study may be applied to other protein-ligand complexes and used to navigate future studies directed at discovering potent and selective inhibitors of APH-type enzymes.     

Comparative study of glutamine synthetase bound lanthanide(III) ions using NMR relaxation and lanthanide(III) luminescence techniques

Changes in the intrinsic fluorescence intensity of glutamine synthetase induced by lanthanide(III) ion binding demonstrate the existence of three types of sites for these ions. The sites are populated sequentially during titrations of the enzyme, and the first two have a stoichiometry of 1 per enzyme subunit. The number of water molecules coordinated to Eu(III) bound to the first site was determined by luminescence lifetime techniques to be 4.1 +/- 0.5. The hydration of Gd(III) bound to the same site was studied by magnetic field dependent water proton longitudinal relaxation rate measurements, and by water proton and deuteron relaxation measurements of one sample at single magnetic fields. The magnetic resonance techniques also yield a value of 4 for the hydration number.     

Active-site gorge and buried water molecules in crystal structures of acetylcholinesterase from Torpedo californica

Buried water molecules and the water molecules in the active-site gorge are analyzed for five crystal structures of acetylcholinesterase from Torpedo californica in the resolution range 2.2-2.5 A (native enzyme, and four inhibitor complexes). A total of 45 buried hydration sites are identified, which are populated with between 36 and 41 water molecules. About half of the buried water is located in a distinct region neighboring the active-site gorge. Most of the buried water molecules are very well conserved among the five structures, and have low displacement parameters, B, of magnitudes similar to those of the main-chain atoms of the central beta-sheet structure. The active-site gorge of the native enzyme is filled with over 20 water molecules, which have poor hydrogen-bond coordination with an average of 2.9 polar contacts per water molecule. Upon ligand binding, distinct groups of these water molecules are displaced, whereas the others remain in positions similar to those that they occupy in the native enzyme. Possible roles of the buried water molecules are discussed, including their possible action as a lubricant to allow large-amplitude fluctuations of the loop structures forming the gorge wall. Such fluctuations are required to facilitate traffic of substrate, products and water molecules to and from the active-site. Because of their poor coordination, the gorge water molecules can be considered as "activated" as compared to bulk water. This should allow their easy displacement by incoming substrate. The relatively loose packing of the gorge water molecules leaves numerous small voids, and more efficient space-filling by substrates and inhibitors may be a major driving force of ligand binding.     

Crystal structures of tRNA-guanine transglycosylase (TGT) in complex with novel and potent inhibitors unravel pronounced induced-fit adaptations and suggest dimer formation upon substrate binding

The bacterial tRNA-guanine transglycosylase (TGT) is a tRNA modifying enzyme catalyzing the exchange of guanine 34 by the modified base preQ1. The enzyme is involved in the infection pathway of Shigella, causing bacterial dysentery. As no crystal structure of the Shigella enzyme is available the homologous Zymomonas mobilis TGT was used for structure-based drug design resulting in new, potent, lin-benzoguanine-based inhibitors. Thorough kinetic studies show size-dependent inhibition of these compounds resulting in either a competitive or non-competitive blocking of the base exchange reaction in the low micromolar range. Four crystal structures of TGT-inhibitor complexes were determined with a resolution of 1.58-2.1 A. These structures give insight into the structural flexibility of TGT necessary to perform catalysis. In three of the structures molecular rearrangements are observed that match with conformational changes also noticed upon tRNA substrate binding. Several water molecules are involved in these rearrangement processes. Two of them demonstrate the structural and catalytic importance of water molecules during TGT base exchange reaction. In the fourth crystal structure the inhibitor unexpectedly interferes with protein contact formation and crystal packing. In all presently known TGT crystal structures the enzyme forms tightly associated homodimers internally related by crystallographic symmetry. Upon binding of the fourth inhibitor the dimer interface, however, becomes partially disordered. This result prompted further analyses to investigate the relevance of dimer formation for the functional protein. Consultation of the available TGT structures and sequences from different species revealed structural and functional conservation across the contacting residues. This suggests that bacterial and eukaryotic TGT could possibly act as homodimers in catalysis. It is hypothesized that one unit of the dimer performs the catalytic reaction whereas the second is required to recognize and properly orient the bound tRNA for the catalytic reaction.     

Mechanism of inhibition of yeast squalene synthase by substrate analog inhibitors.

Squalene synthase catalyzes the reductive condensation of two identical substrate molecules, farnesyl diphosphate, to the hydrocarbon squalene via an obligatory intermediate, presqualene pyrophosphate. Since the kinetic mechanism of the transformation is sequential, two substrate binding pockets that recognize the same molecule must exist in the enzyme active site. This raises the possibility of a choice of binding pockets for inhibitors that are designed as substrate or reaction intermediate analogs and thus may provide some information on the mechanism of differentiation of the two identical molecules. In this report, we have investigated the mechanism of inhibition of a series of farnesyl diphosphate analog inhibitors. The inhibitors fall into two categories. One class of compounds binds to free enzyme as well as the enzyme substrate complex, and the binding is refractory to the concentration of the substrate. The second class binds only to the free enzyme, and its binding is significantly modulated by the substrate concentration. Very modest structural changes in the compounds appear to dictate which class of inhibitor any compound may fall into. The significance of these observations with respect to the mechanism of the enzyme are discussed.     

1,3-Phenylene bis(ketoacid) derivatives as inhibitors of Escherichia coli dihydrodipicolinate synthase

Dihydrodipicolinate synthase is a key enzyme in the lysine biosynthesis pathway that catalyzes the condensation of pyruvate and aspartate semi-aldehyde. A series of phenolic ketoacid derivatives that mimic the proposed enzymatic intermediate were designed as potential inhibitors of this enzyme and were synthesized from simple precursors. The ketoacid derivatives were shown to act as slow and slow-tight binding inhibitors. Mass spectrometric experiments provided further evidence to support the proposed model of inhibition, demonstrating either an encounter complex or a condensation product for the slow and slow-tight binding inhibitors, respectively.     

Crystal structure of a barnase-d(GpC) complex at 1.9 A resolution

The ribonuclease excreted by Bacillus amyloliquefaciens, Barnase, was co-crystallized with the deoxy-dinucleotide d(GpC). The crystal structure was determined by molecular replacement from a model of free Barnase previously derived by Mauguen et al. Refinement was carried out using data to 1.9 A resolution. The final model, which has a crystallographic R factor of 22%, includes 869 protein atoms, 38 atoms from d(GpC), a sulfate ion and 73 water molecules. Only minor differences from free Barnase are seen in the protein moiety, the root-mean-square C alpha movement being 0.45 A. The dinucleotide has a folded conformation. It is located near the active site of the enzyme, but outside the protein molecule and making crystal packing contacts with neighboring molecules. The guanine base is stacked on the imidazole ring of active site His102, rather than binding to the so-called recognition loop as it does in other complexes of guanine nucleotides with microbial nucleases. The deoxyguanosine is syn, with the sugar ring in C-2'-endo conformation; the deoxycytidine is anti and C-4'-exo. In addition to the stacking interaction, His102 hydrogen bonds to the free 5' hydroxyl, which is located near the position where the 3' phosphate group is found in other inhibitors of microbial ribonucleases. While the mode of binding observed with d(GpC) and Barnase would be non-productive for a dinucleotide substrate, it may define a site for the nucleotide product on the 3' side of the hydrolyzed bond.     

Discovering new inhibitors of bacterial glucosamine-6P synthase (GlmS) by docking simulations

Results of an in silico screening of a freely accessible database encompassing 50,000 commercial compounds on bacterial glucosamine-6P synthase (Glms) are described. Each product was docked with the GOLD software in a region of 20A surrounding the sugar binding site and ranked according to its score. Among the 14 best-scored molecules, three molecules exhibited good experimental inhibition properties (IC(50)=70 microM) giving a high hit rate (H.R.: 0.23). Interestingly, these molecules are predicted to interact with a protein region that forms a pocket at the interface between the two enzyme monomers, opening the route to dimerization inhibitors.     

Exploring interfacial water trapping in protein-ligand complexes with multithermal titration calorimetry

In this work, we examine the hypothesis about how trapped water molecules at the interface between triosephosphate isomerase (TIM) and either of two phosphorylated inhibitors, 2-phosphoglycolate (2PG) or phosphoglycolohydroxamate (PGH), can explain the anomalous highly negative binding heat capacities (ΔC_p,b) of both complexes, TIM–2PG and TIM–PGH. We performed fluorimetric titrations of the enzyme with PGH inhibitor under osmotic stress conditions, using various concentrations of either osmolyte: sucrose, ethylene glycol or glycine betaine. We also analyze the binding processes under various stressor concentrations using a novel calorimetric methodology that allows ΔC_p,b determinations in single experiments: Multithermal Titration Calorimetry. The binding constant of the TIM–PGH complex decreased gradually with the concentration of all osmolytes, but at diverse extents depending on the osmolyte nature. According to the osmotic stress theory, this decrease indicates that the number of water molecules associated with the enzyme increases with inhibitor binding, i.e. some solvent molecules became trapped. Additionally, the binding heat capacities became less negative at higher osmolyte concentrations, their final values depending on the osmolyte. These effects were also observed in the TIM–2PG complex using sucrose as stressor. Our results strongly suggest that some water molecules became immobilized when the TIM-inhibitor complexes were formed. A computational analysis of the hydration state of the binding site of TIM in both its free state and its complexed form with 2PG or PGH, based on molecular dynamics (MD) simulations in explicit solvent, showed that the binding site effectively immobilized additional water molecules after binding these inhibitors.     

2-(oxalylamino)-benzoic acid is a general, competitive inhibitor of protein-tyrosine phosphatases

Protein-tyrosine phosphatases (PTPs) are critically involved in regulation of signal transduction processes. Members of this class of enzymes are considered attractive therapeutic targets in several disease states, e.g. diabetes, cancer, and inflammation. However, most reported PTP inhibitors have been phosphorus-containing compounds, tight binding inhibitors, and/or inhibitors that covalently modify the enzymes. We therefore embarked on identifying a general, reversible, competitive PTP inhibitor that could be used as a common scaffold for lead optimization for specific PTPs. We here report the identification of 2-(oxalylamino)-benzoic acid (OBA) as a classical competitive inhibitor of several PTPs. X-ray crystallography of PTP1B complexed with OBA and related non-phosphate low molecular weight derivatives reveals that the binding mode of these molecules to a large extent mimics that of the natural substrate including hydrogen bonding to the PTP signature motif. In addition, binding of OBA to the active site of PTP1B creates a unique arrangement involving Asp(181), Lys(120), and Tyr(46). PTP inhibitors are essential tools in elucidating the biological function of specific PTPs and they may eventually be developed into selective drug candidates. The unique enzyme kinetic features and the low molecular weight of OBA makes it an ideal starting point for further optimization.     

Dimerization of matrix metalloproteinase-2 (MMP-2): functional implication in MMP-2 activation

Matrix metalloproteinase-2 (MMP-2) functions in diverse biological processes through the degradation of extracellular and non-extracellular matrix molecules. Because of its potential for tissue damage, there are several ways to regulate MMP-2 activity, including gene expression, compartmentalization, zymogen activation, and enzyme inactivation by extracellular inhibitors. Enzyme regulation through zymogen activation is important for the regulation of MMP-2 activity. In our previous studies, we showed that thrombin directly cleaved the propeptide of MMP-2 at specific sites for enzyme activation. We also demonstrated that heparan sulfate was required for thrombin-mediated activation of pro-MMP-2 by binding to thrombin, presumably through conformational changes at the active site of the enzyme. This suggests a regulatory mechanism for thrombin-mediated activation of pro-MMP-2. In this study, we found that MMP-2 formed a reduction-sensitive homodimer in a controlled manner and that Ca(2+) ion was essential for homodimerization of MMP-2. Homodimerization was not associated with protein kinase C-mediated phosphorylation of MMP-2. MMP-2 formed a homodimer through an intermolecular disulfide bond between Cys(102) and the neighboring Cys(102). Homodimerization of MMP-2 enhanced thrombin-mediated activation of pro-MMP-2. Moreover, the MMP-2 homodimer could cleave a small peptide substrate without removal of the propeptide. Taken together, our experimental data suggest a novel regulatory mechanism for pro-MMP-2 activation that is modulated through homodimerization of MMP-2.