A Phenylboronic Acid-Based Transition State Analogue Yields Nanomolar Inhibition of Mandelate Racemase

Mandelate racemase (MR) catalyzes the Mg²⁺-dependent interconversion of (R)- and (S)-mandelate by stabilizing the altered substrate in the transition state (TS) by ∼26 kcal/mol. The enzyme has been employed as a model to explore the limits to which the free energy of TS stabilization may be captured by TS analogues to effect strong binding. Herein, we determined the thermodynamic parameters accompanying binding of a series of bromo-, chloro-, and fluoro-substituted phenylboronic acids (PBAs) by MR and found that binding was predominately driven by favorable entropy changes. 3,4-Dichloro-PBA was discovered to be the most potent inhibitor yet identified for MR, binding with a K_d^app value of 11 ± 2 nM and exceeding the binding of the substrate by ∼72,000-fold. The ΔC_p value accompanying binding (-488 ± 18 cal·mol^-1 K^-1) suggested that dispersion forces contribute significantly to the binding. The pH-dependence of the inhibition revealed that MR preferentially binds the anionic, tetrahedral form of 3,4-dichloro-PBA with a pH-independent K_i value of 5.7 ± 0.5 nM, which was consistent with the observed upfield shift of the ¹¹B NMR signal. The linear free energy relationship between log(k_cat/K_m) and log(1/K_i) for wild-type and 11 MR variants binding 3,4-dichloro-PBA had a slope of 0.8 ± 0.2, indicating that MR recognizes the inhibitor as an analogue of the TS. Hence, halogen substitution may be utilized to capture additional free energy of TS stabilization arising from dispersion forces to enhance the binding of boronic acid inhibitors by MR.

Ligand binding free-energy calculations with funnel metadynamics

The accurate resolution of the binding mechanism of a ligand to its molecular target is fundamental to develop a successful drug design campaign. Free-energy calculations, which provide the energy value of the ligand-protein binding complex, are essential for resolving the binding mode of the ligand. The accuracy of free-energy calculation methods is counteracted by their poor user-friendliness, which hampers their broad application. Here we present the Funnel-Metadynamics Advanced Protocol (FMAP), which is a flexible and user-friendly graphical user interface (GUI)-based protocol to perform funnel metadynamics, a binding free-energy method that employs a funnel-shape restraint potential to reveal the ligand binding mode and accurately calculate the absolute ligand-protein binding free energy. FMAP guides the user through all phases of the free-energy calculation process, from preparation of the input files, to production simulation, to analysis of the results. FMAP delivers the ligand binding mode and the absolute protein-ligand binding free energy as outputs. Alternative binding modes and the role of waters are also elucidated, providing a detailed description of the ligand binding mechanism. The entire protocol on the paradigmatic system benzamidine-trypsin, composed of ~105 k atoms, took ~2.8 d using the Cray XC50 piz Daint cluster at the Swiss National Supercomputing Centre.

Binding Constants of Substituted Benzoic Acids with Bovine Serum Albumin

Experimental data on the affinity of various substances to albumin are essential for the development of empirical models to predict plasma binding of drug candidates. Binding of 24 substituted benzoic acid anions to bovine serum albumin was studied using spectrofluorimetric titration. The equilibrium constants of binding at 298 K were determined according to 1:1 complex formation model. The relationships between the ligand structure and albumin affinity are analyzed. The binding constant values for m- and p-monosubstituted acids show a good correlation with the Hammett constants of substituents. Two- and three-parameter quantitative structure–activity relationship (QSAR) models with theoretical molecular descriptors are able to satisfactorily describe the obtained values for the whole set of acids. It is shown that the electron-density distribution in the aromatic ring exerts crucial influence on the albumin affinity.

Thermodynamic profiling for fragment-based lead discovery and optimization

Introduction: The enthalpic and entropic components of the ligand-protein binding free energy reflect the type and quality of the interactions and relate to the physicochemical properties of the ligands. These findings have significance in medicinal chemistry optimizations since they suggest that the thermodynamic profiling of the binding may help monitor and control the unfavorable size and hydrophobicity increase typically accompanying affinity improvements and leading to suboptimal pharmacokinetic properties.Areas covered: This review describes the ligand-protein binding event in terms of elementary steps, their associated interactions, and their enthalpic and entropic consequences. The relationships among the breaking and forming interactions, the binding thermodynamic profile, and the physicochemical properties of the ligands are also discussed.Expert opinion: Analysis of the size dependence of available affinity and favorable enthalpy highlights the limitation of the simultaneous optimization of these quantities. Indeed, moderate, rather than very high affinities can be conciliated with favorable physicochemical and pharmacokinetic profiles as it is supported by the affinity range of historical oral drugs. Although thermodynamic quantities are not suitable endpoints for medicinal chemistry optimizations owing to the complexity of the binding thermodynamics, thermodynamic profiling together with structural studies can be advantageously used to understand the details of the binding process and to optimize it.

A Fragmenting Protocol with Explicit Hydration for Calculation of Binding Enthalpies of Target-Ligand Complexes at a Quantum Mechanical Level

Optimization of the enthalpy component of binding thermodynamics of drug candidates is a successful pathway of rational molecular design. However, the large size and missing hydration structure of target-ligand complexes often hinder such optimizations with quantum mechanical (QM) methods. At the same time, QM calculations are often necessitated for proper handling of electronic effects. To overcome the above problems, and help the QM design of new drugs, a protocol is introduced for atomic level determination of hydration structure and extraction of structures of target-ligand complex interfaces. The protocol is a combination of a previously published program MobyWat, an engine for assigning explicit water positions, and Fragmenter, a new tool for optimal fragmentation of protein targets. The protocol fostered a series of fast calculations of ligand binding enthalpies at the semi-empirical QM level. Ligands of diverse chemistry ranging from small aromatic compounds up to a large peptide helix of a molecular weight of 3000 targeting a leukemia protein were selected for systematic investigations. Comparison of various combinations of implicit and explicit water models demonstrated that the presence of accurately predicted explicit water molecules in the complex interface considerably improved the agreement with experimental results. A single scaling factor was derived for conversion of QM reaction heats into binding enthalpy values. The factor links molecular structure with binding thermodynamics via QM calculations. The new protocol and scaling factor will help automated optimization of binding enthalpy in future molecular design projects.

Broad-scale analysis of thermodynamic signatures in medicinal chemistry: are enthalpy-favored binders the better development option?

Thermodynamic profiles of ligand binding, particularly enthalpically favored binding signatures, have been suggested as a criterion to support the decision-making process around which compounds to select for further optimization in drug development. The concept was enthusiastically taken up, but turned out to be too superficial, either because many aspects determining thermodynamic profiles are insufficiently appreciated or because it is difficult to compare such data on a global scale. The impact of water, changes in protonation states, along with buffer dependencies and incompatible measurement conditions that are far from standard conditions hamper such broad-scale comparisons. However, thermodynamic signatures can make us aware of the impact of these aspects and provide important hints for improving our understanding of the binding process and defining criteria for drug optimization.

Thiophene-Conjugated Ligand Probe for Nonenzymatic Turn-On Electrochemical Protein Detection

A new type of turn-on electrochemical protein detection is developed using an electropolymerizable molecular probe. To detect trypsin, a benzamidine ligand is conjugated with a thiophene moiety. Encapsulation of the probe in the trypsin pocket prevents electropolymerization, leading to efficient electron transfer from the electrolyte to the electrode. In contrast, unbound probes can become electropolymerized, yielding a polythiophene layer on the electrode. The polythiophene formed this way suppressed electron transfer. The detection limit of trypsin using this electrochemical strategy is 50 nM. The method is shown to be useful for nonenzymatic turn-on electrochemical detection.

Influence of Alkylammonium Acetate Buffers on Protein-Ligand Noncovalent Interactions Using Native Mass Spectrometry

We investigate the influence of three volatile alkylammonium acetate buffers on binding affinities for protein-ligand interactions determined by native electrospray ionization-mass spectrometry (ESI-MS). Four different types of proteins were chosen for this study. A charge-reduction effect was observed for all the cases studied, in comparison to the ions formed in ammonium acetate solution. When increasing the collision energy, the complexes of trypsin and the ligand were found to be more stable when sprayed from alkylammonium acetate buffers than from ammonium acetate. The determined dissociation constant (K_d) also exhibited a drop (up to 40%) when ammonium acetate was replaced by alkylammonium acetate buffers for the case of lysozyme and the ligand. The prospective uses of these ammonium acetate analogs in native ESI-MS are discussed in this paper as well. Graphical Abstract ᅟ.

Solvent effects on ligand binding to a serine protease

ITC experiments and MD simulations reveal the mechanism behind enthalpy/entropy compensation upon trypsin-benzamidine binding at different solvation conditions.

Discovery of a Highly Potent, Selective, and Orally Bioavailable Macrocyclic Inhibitor of Blood Coagulation Factor VIIa-Tissue Factor Complex

Inhibitors of the tissue factor (TF)/factor VIIa complex (TF-FVIIa) are promising novel anticoagulants which show excellent efficacy and minimal bleeding in preclinical models. Starting with an aminoisoquinoline P1-based macrocyclic inhibitor, optimization of the P' groups led to a series of highly potent and selective TF-FVIIa inhibitors which displayed poor permeability. Fluorination of the aminoisoquinoline reduced the basicity of the P1 group and significantly improved permeability. The resulting lead compound was highly potent, selective, and achieved good pharmacokinetics in dogs with oral dosing. Moreover, it demonstrated robust antithrombotic activity in a rabbit model of arterial thrombosis.

PyEMMA 2: A Software Package for Estimation, Validation, and Analysis of Markov Models

Markov (state) models (MSMs) and related models of molecular kinetics have recently received a surge of interest as they can systematically reconcile simulation data from either a few long or many short simulations and allow us to analyze the essential metastable structures, thermodynamics, and kinetics of the molecular system under investigation. However, the estimation, validation, and analysis of such models is far from trivial and involves sophisticated and often numerically sensitive methods. In this work we present the open-source Python package PyEMMA ( http://pyemma.org ) that provides accurate and efficient algorithms for kinetic model construction. PyEMMA can read all common molecular dynamics data formats, helps in the selection of input features, provides easy access to dimension reduction algorithms such as principal component analysis (PCA) and time-lagged independent component analysis (TICA) and clustering algorithms such as k-means, and contains estimators for MSMs, hidden Markov models, and several other models. Systematic model validation and error calculation methods are provided. PyEMMA offers a wealth of analysis functions such that the user can conveniently compute molecular observables of interest. We have derived a systematic and accurate way to coarse-grain MSMs to few states and to illustrate the structures of the metastable states of the system. Plotting functions to produce a manuscript-ready presentation of the results are available. In this work, we demonstrate the features of the software and show new methodological concepts and results produced by PyEMMA.

A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation

Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.

Design and Synthesis of Phenylpyrrolidine Phenylglycinamides As Highly Potent and Selective TF-FVIIa Inhibitors

Inhibitors of the Tissue Factor/Factor VIIa (TF-FVIIa) complex are promising novel anticoagulants that show excellent efficacy and minimal bleeding in preclinical models. On the basis of a zwitterionic phenylglycine acylsulfonamide 1, a phenylglycine benzylamide 2 was shown to possess improved permeability and oral bioavailability. Optimization of the benzylamide, guided by X-ray crystallography, led to a potent TF-FVIIa inhibitor 18i with promising oral bioavailability, but promiscuous activity in an in vitro safety panel of receptors and enzymes. Introducing an acid on the pyrrolidine ring, guided by molecular modeling, resulted in highly potent, selective, and efficacious TF-FVIIa inhibitors with clean in vitro safety profile. The pyrrolidine acid 20 showed a moderate clearance, low volume of distribution, and a short t 1/2 in dog PK studies.

Monte Carlo Free Ligand Diffusion with Markov State Model Analysis and Absolute Binding Free Energy Calculations

Obtaining absolute binding free energies from unbiased ligand diffusion has attracted a significant amount of attention due to its implications in drug design. Several studies have used special purpose computers and software to achieve microsecond molecular dynamics which, combined with a Markov state model analysis, are capable of providing absolute binding free energies. We have recently developed a Monte Carlo based technique, PELE, capable of performing a dynamical exploration of the protein-ligand energy landscape including free ligand diffusion into the active site, at a fraction of the computational cost of molecular dynamics techniques. We demonstrate here the capabilities of our Monte Carlo technique in obtaining absolute binding free energies for a series of benzamidine like inhibitors into trypsin. Our results are in good agreement with experimental data and other molecular dynamics simulations, indicating that PELE can be a useful tool for quick estimates of binding free energies and mechanisms.

Funnel metadynamics as accurate binding free-energy method

A detailed description of the events ruling ligand/protein interaction and an accurate estimation of the drug affinity to its target is of great help in speeding drug discovery strategies. We have developed a metadynamics-based approach, named funnel metadynamics, that allows the ligand to enhance the sampling of the target binding sites and its solvated states. This method leads to an efficient characterization of the binding free-energy surface and an accurate calculation of the absolute protein-ligand binding free energy. We illustrate our protocol in two systems, benzamidine/trypsin and SC-558/cyclooxygenase 2. In both cases, the X-ray conformation has been found as the lowest free-energy pose, and the computed protein-ligand binding free energy in good agreement with experiments. Furthermore, funnel metadynamics unveils important information about the binding process, such as the presence of alternative binding modes and the role of waters. The results achieved at an affordable computational cost make funnel metadynamics a valuable method for drug discovery and for dealing with a variety of problems in chemistry, physics, and material science.

Influence of dimehylsulfoxide on protein-ligand binding affinities

Because of its favorable physicochemical properties, DMSO is the standard solvent for sample storage and handling of compounds in drug discovery. To date, little attention was given to how DMSO influences protein-ligand binding strengths. In this study we investigated the effects of DMSO on different noncovalent protein-ligand complexes, in particular in terms of the binding affinities, which we determined using nanoESI-MS. For the investigation, three different protein-ligand complexes were chosen: trypsin-Pefabloc, lysozyme-tri-N-acetylchitotriose (NAG3), and carbonic anhydrase-chlorothiazide. The DMSO content in the nanoESI buffer was increased systematically from 0.5 to 8%. For all three model systems, it was shown that the binding affinity decreases upon addition of DMSO. Even 0.5-1% DMSO alters the KD values, in particular for the tight binding system carbonic anhydrase-chlorothiazide. The determined dissociation constant (KD) is up to 10 times higher than for a DMSO-free sample in the case of carbonic anhydrase-chlorothiazide binding. For the trypsin-Pefabloc and lysozyme-NAG3 complexes, the dissociation constants are 7 and 3 times larger, respectively, in the presence of DMSO. This work emphasizes the importance of effects of DMSO as a co-solvent for quantification of protein-ligand binding strengths in the early stages of drug discovery.

Sanjeevini: a freely accessible web-server for target directed lead molecule discovery

Background

Computational methods utilizing the structural and functional information help to understand specific molecular recognition events between the target biomolecule and candidate hits and make it possible to design improved lead molecules for the target.

Results

Sanjeevini represents a massive on-going scientific endeavor to provide to the user, a freely accessible state of the art software suite for protein and DNA targeted lead molecule discovery. It builds in several features, including automated detection of active sites, scanning against a million compound library for identifying hit molecules, all atom based docking and scoring and various other utilities to design molecules with desired affinity and specificity against biomolecular targets. Each of the modules is thoroughly validated on a large dataset of protein/DNA drug targets.

Conclusions

The article presents Sanjeevini, a freely accessible user friendly web-server, to aid in drug discovery. It is implemented on a tera flop cluster and made accessible via a web-interface at http://www.scfbio-iitd.res.in/sanjeevini/sanjeevini.jsp. A brief description of various modules, their scientific basis, validation, and how to use the server to develop in silico suggestions of lead molecules is provided.

Quantifying protein-ligand binding constants using electrospray ionization mass spectrometry: a systematic binding affinity study of a series of hydrophobically modified trypsin inhibitors

NanoESI-MS is used for determining binding strengths of trypsin in complex with two different series of five congeneric inhibitors, whose binding affinity in solution depends on the size of the P3 substituent. The ligands of the first series contain a 4-amidinobenzylamide as P1 residue, and form a tight complex with trypsin. The inhibitors of the second series have a 2-aminomethyl-5-chloro-benzylamide as P1 group, and represent a model system for weak binders. The five different inhibitors of each group are based on the same scaffold and differ only in the length of the hydrophobic side chain of their P3 residue, which modulates the interactions in the S3/4 binding pocket of trypsin. The dissociation constants (K(D)) for high affinity ligands investigated by nanoESI-MS ranges from 15 nM to 450 nM and decreases with larger hydrophobic P3 side chains. Collision-induced dissociation (CID) experiments of five trypsin and benzamidine-based complexes show a correlation between trends in K(D) and gas-phase stability. For the second inhibitor series we could show that the effect of imidazole, a small stabilizing additive, can avoid the dissociation of the complex ions and as a result increases the relative abundance of weakly bound complexes. Here the K(D) values ranging from 2.9 to 17.6 μM, some 1-2 orders of magnitude lower than the first series. For both ligand series, the dissociation constants (K(D)) measured via nanoESI-MS were compared with kinetic inhibition constants (K(i)) in solution.

Virtual screening using molecular simulations

Effective virtual screening relies on our ability to make accurate prediction of protein-ligand binding, which remains a great challenge. In this work, utilizing the molecular-mechanics Poisson-Boltzmann (or Generalized Born) surface area approach, we have evaluated the binding affinity of a set of 156 ligands to seven families of proteins, trypsin β, thrombin α, cyclin-dependent kinase (CDK), cAMP-dependent kinase (PKA), urokinase-type plasminogen activator, β-glucosidase A, and coagulation factor Xa. The effect of protein dielectric constant in the implicit-solvent model on the binding free energy calculation is shown to be important. The statistical correlations between the binding energy calculated from the implicit-solvent approach and experimental free energy are in the range of 0.56-0.79 across all the families. This performance is better than that of typical docking programs especially given that the latter is directly trained using known binding data whereas the molecular mechanics is based on general physical parameters. Estimation of entropic contribution remains the barrier to accurate free energy calculation. We show that the traditional rigid rotor harmonic oscillator approximation is unable to improve the binding free energy prediction. Inclusion of conformational restriction seems to be promising but requires further investigation. On the other hand, our preliminary study suggests that implicit-solvent based alchemical perturbation, which offers explicit sampling of configuration entropy, can be a viable approach to significantly improve the prediction of binding free energy. Overall, the molecular mechanics approach has the potential for medium to high-throughput computational drug discovery.

Absolute binding free energy calculations: on the accuracy of computational scoring of protein-ligand interactions

Calculating the absolute binding free energies is a challenging task. Reliable estimates of binding free energies should provide a guide for rational drug design. It should also provide us with deeper understanding of the correlation between protein structure and its function. Further applications may include identifying novel molecular scaffolds and optimizing lead compounds in computer-aided drug design. Available options to evaluate the absolute binding free energies range from the rigorous but expensive free energy perturbation to the microscopic linear response approximation (LRA/beta version) and related approaches including the linear interaction energy (LIE) to the more approximated and considerably faster scaled protein dipoles Langevin dipoles (PDLD/S-LRA version) as well as the less rigorous molecular mechanics Poisson-Boltzmann/surface area (MM/PBSA) and generalized born/surface area (MM/GBSA) to the less accurate scoring functions. There is a need for an assessment of the performance of different approaches in terms of computer time and reliability. We present a comparative study of the LRA/beta, the LIE, the PDLD/S-LRA/beta, and the more widely used MM/PBSA and assess their abilities to estimate the absolute binding energies. The LRA and LIE methods perform reasonably well but require specialized parameterization for the nonelectrostatic term. The PDLD/S-LRA/beta performs effectively without the need of reparameterization. Our assessment of the MM/PBSA is less optimistic. This approach appears to provide erroneous estimates of the absolute binding energies because of its incorrect entropies and the problematic treatment of electrostatic energies. Overall, the PDLD/S-LRA/beta appears to offer an appealing option for the final stages of massive screening approaches.

A comprehensive examination of the contributions to the binding entropy of protein-ligand complexes

One of the most important requirements in computer-aided drug design is the ability to reliably evaluate the binding free energies. However, the process of ligand binding is very complex because of the intricacy of the interrelated processes that are difficult to predict and quantify. In fact, the deeper understanding of the origin of the observed binding free energies requires the ability to decompose these free energies to their contributions from different interactions. Furthermore, it is important to evaluate the relative entropic and enthalpic contributions to the overall free energy. Such an evaluation is useful for assessing temperature effects and exploring specialized options in enzyme design. Unfortunately, calculations of binding entropies have been much more challenging than calculations of binding free energies. This work is probably the first to present microscopic evaluation of all of the relevant components to the binding entropy, namely configurational, polar solvation, and hydrophobic entropies. All of these contributions are evaluated by the restraint release approach. The calculated results shed an interesting light on major compensation effects in both the solvation and hydrophobic effect and, despite some overestimate, can provide very useful insight. This study also helps in analyzing some problems with the widely used molecular mechanics/Poisson-Boltzmann surface area approach.

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.

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.

Investigation of the interactions of quercetin and morin with trypsin

The interactions of quercetin and morin with trypsin were investigated by UV-vis absorption, fluorescence, synchronous fluorescence and three-dimensional fluorescence spectra techniques under physiological pH 7.40. Quercetin and morin effectively quenched the intrinsic fluorescence of trypsin via static quenching. The process of binding quercetin and morin on trypsin was a spontaneous molecular interaction procedure. The binding constants and thermodynamic parameters at two different temperatures, the binding locality and the binding power were obtained. The conformation of trypsin was discussed by synchronous and three-dimensional fluorescence techniques.

Molecular mechanism of interaction between norfloxacin and trypsin studied by molecular spectroscopy and modeling

The molecular mechanism of the binding of norfloxacin (NRF) to trypsin was investigated by fluorescence, synchronous fluorescence and UV-vis absorbance spectroscopy and molecular modeling at physiological conditions. The quenching mechanism and the binding mode were investigated in terms of the association constants and basic thermodynamic parameters. The results of spectroscopic measurements suggested that NRF have a strong ability to quench the intrinsic fluorescence of trypsin through static quenching procedure. Moreover, fluorescence experiments were also performed at different values of pH to elucidate the effect of pH on the binding. The NRF-trypsin complex was stabilized by hydrophobic forces and hydrogen bonding, via tryptophan residue as indicated from the thermodynamic parameters, which was consistent with the results of molecular docking and accessible surface area calculations.

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.

Standard Free Energy of Binding from a One-Dimensional Potential of Mean Force

A practical approach that enables one to calculate the standard free energy of binding from a one-dimensional potential of mean force (PMF) is proposed. Umbrella sampling and the weighted histogram analysis method are used to generate a PMF along the reaction coordinate of binding. At each point, a restraint is applied orthogonal to the reaction coordinate to make possible the determination of the volume sampled by the ligand. The free energy of binding from an arbitrary unbound volume to the restrained bound form is calculated from the ratio of the PMF integrated over the bound region to that of the unbound. Adding the free energy changes from the standard-state volume to the unbound volume and from the restrained to the unrestrained bound state gives the standard free energy of binding. Exploration of the best choice of binding paths is also made. This approach is first demonstrated on a model binding system and then tested on the benzamidine-trypsin system for which reasonable agreement with experiment is found. A comparison is made with other methods to obtain the standard free energy of binding from the PMF.

DNA-small molecule chimera with responsive protein-binding ability

In this communication, we disclose a generalizable strategy for developing agents with regulable protein-binding ability. In particular, a responsive DNA-small molecule chimera (DC) 1 consisting of two synthetic protein-binding arms and a core oligonucleotide (ODN) domain is discussed. DC 1 can be cycled from a bidentate intramolecular quadruplex form to a monodentate duplex structure, via addition of external ODN stimuli. Importantly, these distinct secondary structures of 1 lead to significantly different protein-binding abilities, with the bidentate conformation showing a 20-fold enhancement (with a 0.8 microM dissociation constant, Kd) in trypsin-binding potency.

Method for Stabilizing Protein−Ligand Complexes in Nanoelectrospray Ionization Mass Spectrometry

The interaction between the bovine pancrease trypsin (Tryp) and its competitive inhibitor benzamidine (1), in solution and the gas phase, is investigated using nanoflow electrospray ionization (nanoES) and Fourier transform ion cyclotron resonance mass spectrometry. In a recent study (Clark, S.M.; Konermann L. Anal. Chem. 2004, 76, 7077-7083), it was reported that the (Tryp + 1) complex could not be detected by ES-MS. Here, it is shown that, with gentle sampling conditions, it is possible to detect gaseous protonated ions of the (Tryp + 1) complex with nanoES-MS. However, the relative abundance of the detected (Tryp + 1)n+ ions is lower than expected, based on solution composition, which suggests that dissociation of (Tryp + 1)n+ ions occurs during MS sampling. The dissociation pathways and corresponding Arrhenius parameters for the protonated (Tryp + 1)n+ ions, at n = 7-9, are determined from time-resolved thermal dissociation experiments, implemented with the blackbody infrared radiative dissociation technique. The gaseous (Tryp + 1)n+ ions are found to have short lifetimes, e.g., <0.6 s, at temperatures of >100 degrees C. The use of solution additives, including polyols, carbohydrates, amino acids, and small organic molecules, to stabilize the (Tryp + 1)n+ ions during nanoES-MS analysis is investigated. Notably, the addition of imidazole to the nanoES solution is shown to preserve the (Tryp + 1)n+ ions. The Kassoc value, (1.9 +/- 0.2) x 104 M-1, determined for the (Tryp + 1) complex by the direct ES-MS method is in agreement with values determined by other analytical methods. The stabilizing effect of imidazole in nanoES-MS is further demonstrated for the interaction between carbonic anhydrase II and 5-(dimethylamino)naphthalene-1-sulfonamide. The stabilizing effect of imidazole may be due to enhanced evaporative cooling achieved by the dissociation of molecules of imidazole, bound nonspecifically, from the protein-ligand complex in the ion source.

Oxalone and Lactone Moieties of Podophyllotoxin Exhibit Properties of Both the B and C Rings of Colchicine in Its Binding with Tubulin

Thermodynamics of podophyllotoxin binding to tubulin and its multiple points of attachment with tubulin has been studied in detail using isothermal titration calorimetry. The calorimetric enthalpy of the association of podophyllotoxin with tubulin is negative and occurs with a negative heat capacity change (DeltaC(p) = -2.47 kJ mol(-)(1) K(-)(1)). The binding is unique with a simultaneous participation of both hydrophobic and hydrogen-bonding forces with unfavorable negative entropic contribution at higher temperature, favored with an enthalpy-entropy compensation. Interestingly, the binding of 2-methoxy-5-(2',3',4'-trimethoxyphenyl)tropone (AC, a colchicine analogue without the B ring) with tubulin is enthalpy-favored. However, the podophyllotoxin-tubulin association depending upon the temperature of the reaction has a favorable entropic and enthalpic component, which resembles both B- and C-ring properties of colchicine. On the basis of the crystal structure of the podophyllotoxin-tubulin complex, distance calculations have indicated a possible interaction between threonine 179 of alpha-tubulin and the hydroxy group on the D ring of podophyllotoxin. To confirm the involvement of the oxalone moiety as well as the lactone ring of podophyllotoxin in tubulin binding, analogues of podophyllotoxin are synthesized with methoxy substitution at the 4' position of ring D along with its isomer and another analogue epimerized at ring E. From these results, involvement of oxalone as well as the lactone ring of the drug in a specific orientation inclusive of ring A is indicated for podophyllotoxin-tubulin binding. Therefore, podophyllotoxin, like colchicine, behaves as a bifunctional ligand having properties of both the B and C rings of colchicine by making more than one point of attachment with the protein tubulin.

Thermodynamic evaluation and modeling of proton and water exchange associated with benzamidine and berenil binding to ß-trypsin

Serine-proteases are involved in vital processes in virtually all species. They are important targets for researchers studying the relationships between protein structure and activity, for the rational design of new pharmaceuticals. Trypsin was used as a model to assess a possible differential contribution of hydration water to the binding of two synthetic inhibitors. Thermodynamic parameters for the association of bovine beta-trypsin (homogeneous material, observed 23,294.4 +/- 0.2 Da, theoretical 23,292.5 Da) with the inhibitors benzamidine and berenil at pH 8.0, 25 degrees C and with 25 mM CaCl2, were determined using isothermal titration calorimetry and the osmotic stress method. The association constant for berenil was about 12 times higher compared to the one for benzamidine (binding constants are K = 596,599 +/- 25,057 and 49,513 +/- 2,732 M(-1), respectively; the number of binding sites is the same for both ligands, N = 0.99 +/- 0.05). Apparently the driving force responsible for this large difference of affinity is not due to hydrophobic interactions because the variation in heat capacity (DeltaCp), a characteristic signature of these interactions, was similar in both systems tested (-464.7 +/- 23.9 and -477.1 +/- 86.8 J K(-1) mol(-1) for berenil and benzamidine, respectively). The results also indicated that the enzyme has a net gain of about 21 water molecules regardless of the inhibitor tested. It was shown that the difference in affinity could be due to a larger number of interactions between berenil and the enzyme based on computational modeling. The data support the view that pharmaceuticals derived from benzamidine that enable hydrogen bond formation outside the catalytic binding pocket of beta-trypsin may result in more effective inhibitors.

Biochemical and molecular modeling analysis of the ability of two p-aminobenzamidine-based sorbents to selectively purify serine proteases (fibrinogenases) from snake venoms

Snake venoms contain several trypsin-like enzymes with equivalent physicochemical characteristics and similar inhibition profiles. These are rather difficult to separate by classical purification procedures and therefore constitute a good model for affinity chromatography analysis. Some of these trypsin homologues present fibrinogenase activity, mimicking one or more features of the central mammalian coagulation enzyme, thrombin. It was previously demonstrated that a number of amidine derivatives are able to interact specifically with some of these serine proteases. To understand the enzyme-sorbent interactions we have investigated the ability of two commercially available benzamidine affinity matrices to purify thrombin-like serine proteases (TLSP) with similar biological properties from two snake venoms (Bothrops jararacussu and Lachesis muta rhombeata). Curiously, each sorbent retained a single but distinct TLSP from each venom with high yield. Molecular modeling analysis suggested that hydrophobic interactions within a specific region on the surface of these enzymes could be generated to explain this exquisite specificity. In addition, it was demonstrated that a specific tandem alignment of the two benzamidine sorbents enables the purification of three other enzymes from B. jararacussu venom.

Sulfonamide drugs binding to the colchicine site of tubulin: thermodynamic analysis of the drug-tubulin interactions by isothermal titration calorimetry

The discovery of several sulfonamide drugs paved the way toward the synthesis of 6 (N-[2-[(4-hydroxyphenyl)amino]-3-pyridinyl]-4-methoxybenzenesulfonamide, E7010) and 7 (N-(3-fluoro-4-methoxyphenyl)pentafluorobenzenesulfonamide, T138067), both of which inhibit tubulin polymerization and are under clinical development. A series of diarylsulfonamides containing an indole scaffold was also found to have antimitotic properties, but their mode of interactions with tubulin has remained unidentified so far. In this study, we demonstrate that these sulfonamide drugs bind to the colchicine site of tubulin in a reversible manner. They quenched intrinsic tryptophan fluorescence of tubulin presumably due to drug-induced conformational changes in the protein, but were unable to modulate GTPase activity of tubulin in contrast to colchicine that enhances the same enzymatic activity. Further investigation using isothermal titration calorimetry (ITC) revealed that 5 (N-(5-chloro-7-indolyl)-4-methoxybenzenesulfonamide) afforded a large positive value of heat capacity change (DeltaC(p)() = +264 cal mol(-1) K(-1)) on binding to tubulin, suggesting a substantial conformational transition in the protein along with partial enthalpy-entropy compensation. On the other hand, the 2-chloro regioisomer 2 gave a large negative value of DeltaC(p)() (-589 cal mol(-1) K(-1)) along with complete enthalpy-entropy compensation. This thermodynamic profile was thought to be attributable to a prominent contribution of van der Waals interaction and hydrogen bonding between specific groups in the drug-tubulin complex. These results indicate that a mere alteration in the position of a single substituent chlorine on the indole scaffold has a great influence on the drug-tubulin binding thermodynamics.

Kinetic analysis of zymogen autoactivation in the presence of a reversible inhibitor

Limited proteolysis is a highly specific irreversible process, which can serve to initiate physiological function by converting a precursor protein into a biologically active form. When the activating enzyme and the activated enzyme coincide, the process is an autocatalytic zymogen activation (i.e. reactions in which the zymogens serves as a substrate for the corresponding active enzyme). The activity of proteases is frequently regulated by the binding of specific protease inhibitors. Thus, to understand the biological regulation of proteolysis, one must understand the role of protease inhibitors. In the present study, a detailed kinetic analysis of autocatalytic reaction modulated by a reversible inhibitor is represented. On the basis of the kinetic equation, a novel procedure is developed to evaluate the kinetic parameters of the reaction. As an example of the application of this method, effects of acetamidine, p-amidinobenzamidine and benzamidine on the autoactivation of trypsinogen by trypsin were studied.

Kinetic and equilibrium characterization of the interaction between bovine trypsin and I-ovalbumin

Serpins are a superfamily of structurally linked proteins with interesting functional properties. Most serpins act as proteinase suicide inhibitors and play a key role in a number of physiological processes. Structural flexibility properties make serpins extremely available to conformational transitions, often causing changes in protein function. Ovalbumin is a member of the serpin family that is not able to inhibit serine proteinases in its native form. In contrast, I-ovalbumin, the product of a heating transition, is a potent reversible serine proteinase inhibitor. In this paper, a detailed equilibrium and kinetic characterization of the interaction between the serpin ovalbumin and bovine trypsin, using a resonant mirror technique, is reported. This methodology revealed that the high affinity interaction between the two binding partners is characterized by high kinetic association constants and low kinetic dissociation constants. The modulation exerted by protons in solution, examined taking into account structural motifs characterizing the binding interface between the two macromolecules, suggests an interaction reminiscent of that observed for the binding of other serine proteases to their ligands.

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

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

Determination of Ligand−Protein Dissociation Constants by Electrospray Mass Spectrometry-Based Diffusion Measurements

A novel approach for the quantification of ligand-protein interactions is presented. Electrospray ionization mass spectrometry (ESI-MS) is used to monitor the diffusion behavior of noncovalent ligands in the presence of their protein receptors. These data allow the fraction of free ligand in solution to be determined, such that the corresponding dissociation constants can be calculated. A set of conditions is developed that provides an "allowable range" of concentrations for this type of assay. The method is tested by applying it to two different inhibitor-enzyme systems. The dissociation constants measured for benzamidine-trypsin and for N,N',N' '-triacetylchitotriose-lysozyme are (50 +/- 10) and (6 +/- 1) mM, respectively. Both of these results are in good agreement with previous data from the literature. In contrast to traditional ESI-MS-based methods, the approach used in this work does not rely on the preservation of specific solution-type noncovalent interactions in the gas phase. It is shown that this method allows an accurate determination of dissociation constants, even in cases in which the ion abundance ratio of free to ligand-bound protein in ESI-MS does not reflect the corresponding concentration ratio in solution.

Prevention of interfacial inactivation of enzymes by coating the enzyme surface with dextran-aldehyde

Interactions between soluble enzymes and interfaces of organic solvent drops or gas bubbles have a very negative effect on the operational stability of the soluble enzymes. In this study, the formation of a hydrophilic shell around the enzyme has been attempted using dextran-aldehyde which would prevent the interaction between enzyme and hydrophobic interfaces with minimal modification of the enzyme surface. After optimizing the size of the dextran (that was found to play a critical role), three different enzymes (glucose oxidase, d-amino acid oxidase, and trypsin) have been conjugated with dextran-aldehyde and their stability towards organic-aqueous and air-liquid interfaces has been evaluated. The treatment itself proved to be very low-cost in terms of activity and was highly stabilizing for the three enzymes assayed. The conjugated preparation of the three assayed enzymes remained fully active in the presence of air-liquid interfaces for at least 10h. However, the unmodified enzymes lost more than 50% of activity within the first hour of the experiments except for trypsin which kept 38% activity after 12h while the trypsin dextran-aldehyde conjugate maintained 100% enzyme activity. Similar results were achieved in the presence of stirred organic solvent-aqueous buffer biphasic system, although in this case some activity was lost by the action of the soluble portion of the organic solvent. In fact, this treatment seems to be also effective to improve the resistance to the action of organic solvent.

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

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

Trypsin revisited: crystallography AT (SUB) atomic resolution and quantum chemistry revealing details of catalysis

A series of crystal structures of trypsin, containing either an autoproteolytic cleaved peptide fragment or a covalently bound inhibitor, were determined at atomic and ultra-high resolution and subjected to ab initio quantum chemical calculations and multipole refinement. Quantum chemical calculations reproduced the observed active site crystal structure with severe deviations from standard stereochemistry and indicated the protonation state of the catalytic residues. Multipole refinement directly revealed the charge distribution in the active site and proved the validity of the ab initio calculations. The combined results confirmed the catalytic function of the active site residues and the two water molecules acting as the nucleophile and the proton donor. The crystal structures represent snapshots from the reaction pathway, close to a tetrahedral intermediate. The de-acylation of trypsin then occurs in true SN2 fashion.

Understanding binding affinity: a combined isothermal titration calorimetry/molecular dynamics study of the binding of a series of hydrophobically modified benzamidinium chloride inhibitors to trypsin

The binding of a series of p-alkylbenzamidinium chloride inhibitors to the serine proteinase trypsin over a range of temperatures has been studied using isothermal titration (micro)calorimetry and molecular dynamics simulation techniques. The inhibitors have small structural variations at the para position of the benzamidinium ion. They show small differences in relative binding affinity but large compensating differences in enthalpy and entropy. Binding affinity decreases with increased branching at the first carbon but increases with increasing the length of a linear alkyl substituent, suggesting that steric hindrance and hydrophobic interactions play dominant roles in binding. Structural analysis showed that the backbone of the enzyme was unaffected by the change of the para substituent. In addition, binding does not correlate strongly with octanol/water partition data. To further characterize this system, the change in the heat capacity on binding, the change in solvent-accessible surface area on binding, the effect of inhibitor binding on the hydration of the active site, the pK(a) of His57, and interactions within the catalytic triad have been investigated. Although the changes in inhibitor structure are small, it is demonstrated that simple concepts such as steric hindrance, hydrophobicity, and buried surface area are insufficient to explain the binding data. Other factors, such as access to the binding site and the cost of dehydration of the active site, are of equal or greater importance.