The application of thermodynamic methods in drug design

The optimization of lead compounds as viable drug candidates involves the optimization of their binding affinity towards the selected target. The binding affinity, K_a, is determined by the Gibbs energy of binding, ΔG, which in turn is determined by the enthalpy, ΔH, and entropy, ΔS, changes (ΔG=ΔH−TΔS). In principle, many combinations of ΔH and ΔS values can give rise to the same ΔG value and, therefore, elicit the same binding affinity. However, enthalpically dominated ligands do not behave the same as entropically dominated ligands. Current paradigms in drug design usually generate highly hydrophobic and conformationally constrained ligands. The thermodynamic signature of these ligands is an entropically dominated binding affinity often accompanied by an unfavorable binding enthalpy. Conformationally constrained ligands cannot easily adapt to changes in the geometry of the binding site, being therefore highly susceptible to drug resistance mutations or naturally occurring genetic polymorphisms. The design of ligands with the capability to adapt to a changing target requires the introduction of certain elements of flexibility or the relaxation of some conformational constraints. Since these compounds pay a larger conformational entropy penalty upon binding, the optimization of their binding affinity requires the presence of a favorable binding enthalpy. In this paper, experimental and computational strategies aimed at identifying and optimizing enthalpic ligands will be discussed and applied to the case of HIV-1 protease inhibitors. It is shown that a thermodynamic guide to drug design permits the identification of drug candidates with a lower susceptibility to target mutations causing drug resistance.     

Thermodynamic analysis of a hydrophobic binding site: probing the PDZ domain with nonproteinogenic peptide ligands

[structure: see text] Isothermal titration calorimetry (ITC) is used to study the thermodynamic consequences of systematically modifying the hydrophobic character of a single residue in a series of protein-binding ligands. By substituting standard and nonproteinogenic aliphatic amino acids for the C-terminal valine of the hexapeptide KKETEV, binding to the third PDZ domain (PDZ3) of the PSD-95 protein is characterized by distinct changes in the Gibbs free energy (DeltaG), enthalpy (DeltaH), and entropy (TDeltaS) parameters. One notable observation is that peptide binding affinity can be improved with a nonstandard residue.     

Concepts in receptor optimization: targeting the RGD peptide

Synthetic receptors have a wide range of potential applications, but it has been difficult to design low molecular weight receptors that bind ligands with high, "proteinlike" affinities. This study uses novel computational methods to understand why it is hard to design a high-affinity receptor and to explore the limits of affinity, with the bioactive peptide RGD as a model ligand. The M2 modeling method is found to yield excellent agreement with experiment for a known RGD receptor and then is used to analyze a series of receptors generated in silico with a de novo design algorithm. Forces driving binding are found to be systematically opposed by proportionate repulsions due to desolvation and entropy. In particular, strong correlations are found between Coulombic attractions and the electrostatic desolvation penalty and between the mean energy change on binding and the cost in configurational entropy. These correlations help explain why it is hard to achieve high affinity. The change in surface area upon binding is found to correlate poorly with affinity within this series. Measures of receptor efficiency are formulated that summarize how effectively a receptor uses surface area, total energy, and Coulombic energy to achieve affinity. Analysis of the computed efficiencies suggests that a low molecular weight receptor can achieve proteinlike affinity. It is also found that macrocyclization of a receptor can, unexpectedly, increase the entropy cost of binding because the macrocyclic structure further restricts ligand motion.     

EF-Tu binding peptides identified,dissected, and affinity optimized by phage display

The highly abundant GTP binding protein elongation factor Tu (EF-Tu) fulfills multiple roles in bacterial protein biosynthesis. Phage-displayed peptides with high affinity for EF-Tu were selected from a library of approximately 4.7 x 10(11) different peptides. The lack of sequence homology among the identified EF-Tu ligands demonstrates promiscuous peptide binding by EF-Tu. Homolog shotgun scanning of an EF-Tu ligand was used to dissect peptide molecular recognition by EF-Tu. All homolog shotgun scanning selectants bound to EF-Tu with higher affinity than the starting ligand. Thus, homolog shotgun scanning can simultaneously optimize binding affinity and rapidly provide detailed structure activity relationships for multiple side chains of a polypeptide ligand. The reported peptide ligands do not compete for binding to EF-Tu with various antibiotic EF-Tu inhibitors, and could identify an EF-Tu peptide binding site distinct from the antibiotic inhibitory sites.     

Capillary electrophoretic methods applied to the investigation of peptide complexes

This article gives an overview of the applications of capillary electrophoretic methods to investigate the non‐covalent interactions of peptides (peptide complexes) with variable middle‐ and high‐molecular‐mass receptors (ligands) as well as with small ions and molecules in the period 2007–2014. Different modes of capillary electrophoretic methods, such as mobility shift (vacancy) affinity capillary electrophoresis, multiple injection affinity capillary electrophoresis, partial filling affinity capillary electrophoresis, Hummel–Dryer method, vacancy peak method and (continuous) frontal analysis capillary electrophoresis, are briefly described and their applicability to determination of binding constants of peptide complexes is discussed. In addition, the detailed experimental conditions of individual applications and the values of binding constants of the particular peptide complexes are presented.     

Adsorption of antimicrobial indolicidin-derived peptides on hydrophobic surfaces

The hydrophobic interaction between antimicrobial peptides and membrane hydrophobic cores is usually related to their cytotoxicity. In this study, the adsorption mechanism of five plasma membrane-associated peptides, indolicidin (IL) and its four derivatives, with hydrophobic ligands was investigated to understand the relationship between peptide hydrophobicity and bioactivity. The hydrophobic adsorption mechanisms of IL and its derivatives were interpreted thermodynamically and kinetically by reversed-phase chromatography (RPC) analysis and surface plasmon resonance (SPR) measurement, respectively. IL and its derivatives possess a similar random coil structure in both aqueous and organic solvents. Thermodynamic analysis showed that the binding enthalpy of peptides with higher electropositivity was lower than those with lower electropositivity and exhibited unfavorable binding entropy. Higher electropositivity peptides adsorbed to the hydrophobic surface arising from the less bound solvent on the peptide surface. A comparison with the kinetic analysis showed that IL and its derivatives adopt a two-state binding model (i.e., adsorption onto and self-association on the hydrophobic acyl chain) to associate with the hydrophobic surface, and the binding affinity of peptide self-association correlates well with peptide hemolysis. Consequently, this study provided a novel concept for understanding the action of plasma membrane-associated peptides.     

Tetra- and Hexavalent Siglec-8 Ligands Modulate Immune Cell Activation

Carbohydrate-binding proteins are generally characterized by poor affinities for their natural glycan ligands, predominantly due to the shallow and solvent-exposed binding sites. To overcome this drawback, nature has exploited multivalency to strengthen the binding by establishing multiple interactions simultaneously. The development of oligovalent structures frequently proved to be successful, not only for proteins with multiple binding sites, but also for proteins that possess a single recognition domain. Herein we present the syntheses of a number of oligovalent ligands for Siglec-8, a monomeric I-type lectin found on eosinophils and mast cells, alongside the thermodynamic characterization of their binding. While the enthalpic contribution of each binding epitope was within a narrow range to that of the monomeric ligand, the entropy penalty increased steadily with growing valency. Additionally, we observed a successful agonistic binding of the tetra- and hexavalent and, to an even larger extent, multivalent ligands to Siglec-8 on immune cells and modulation of immune cell activation. Thus, triggering a biological effect is not restricted to multivalent ligands but could be induced by low oligovalent ligands as well, whereas a monovalent ligand, despite binding with similar affinity, showed an antagonistic effect.     

Minimizing the Entropy Penalty for Ligand Binding: Lessons from the Molecular Recognition of the Histo Blood‐Group Antigens by Human Galectin‐3

Ligand conformational entropy plays an important role in carbohydrate recognition events. Glycans are characterized by intrinsic flexibility around the glycosidic linkages, thus in most cases, loss of conformational entropy of the sugar upon complex formation strongly affects the entropy of the binding process. By employing a multidisciplinary approach combining structural, conformational, binding energy, and kinetic information, we investigated the role of conformational entropy in the recognition of the histo blood‐group antigens A and B by human galectin‐3, a lectin of biomedical interest. We show that these rigid natural antigens are pre‐organized ligands for hGal‐3, and that restriction of the conformational flexibility by the branched fucose (Fuc) residue modulates the thermodynamics and kinetics of the binding process. These results highlight the importance of glycan flexibility and provide inspiration for the design of high‐affinity ligands as antagonists for lectins.     

Selection of Peptide Ligands Specific for Baculovirus DNA‐Binding Protein from the Flitrx™ Random Peptide Display Library

Abstract

Autographa californica nucleopolyhedrovirus (AcNPV) is a baculovirus that is widely employed as a vector for the expression of foreign genes and pest control. Although baculoviruses, including AcNPV, efficiently replicate in the nuclei of arthropod cells, the dynamics and mechanism of DNA replication within the infected cell are still poorly understood. It has been found that the DNA‐binding protein (DBP) is an early gene product and appears to be crucial for viral DNA replication.

Presented here is the selection of peptide ligands that specifically bind to DBP for AcNPV from the FliTrx? random peptide display library; this entails the amplification, cloning of the DNA‐binding protein (DBP) gene from AcNPV and the construction of the expression plasmid for DBP, and the expression and purification of the recombinant His.Tag AcNPV DBP that was used as a target molecule for the selection of the peptide ligands specific for AcNPV DBP. The affinity and efficiency of such peptide ligands were then measured by ELISA procedures.

The beneficial aspect of this research is the monospecificity quality of the peptide ligands specific for AcNPV DBP. They could be used for the study of the dynamics of the viral genome and its replication within the infected cell, for the development of a quantitative method for the determination of the presence of baculovirus in various samples, for the development of a peptide ligandbased assay for the determination of baculovirus titers; or they could be immobilized on a chromatographic support for an improved affinity purification of AcNPV DBP.     

Side-chain conformational space analysis (SCSA): A multi conformation-based QSAR approach for modeling and prediction of protein–peptide binding affinities

In this article, the concept of multi conformation-based quantitative structure-activity relationship (MCB-QSAR) is proposed, and based upon that, we describe a new approach called the side-chain conformational space analysis (SCSA) to model and predict protein-peptide binding affinities. In SCSA, multi-conformations (rather than traditional single-conformation) have received much attention, and the statistical average information on multi-conformations of side chains is determined using self-consistent mean field theory based upon side chain rotamer library. Thereby, enthalpy contributions (including electrostatic, steric, hydrophobic interaction and hydrogen bond) and conformational entropy effects to the binding are investigated in terms of occurrence probability of residue rotamers. Then, SCSA was applied into the dataset of 419 HLA-A 0201 binding peptides, and nonbonding contributions of each position in peptide ligands are well determined. For the peptides, the hydrogen bond and electrostatic interactions of the two ends are essential to the binding specificity, van der Waals and hydrophobic interactions of all the positions ensure strong binding affinity, and the loss of conformational entropy at anchor positions partially counteracts other favorable nonbonding effects.     

Additivity and the physical basis of multivalency effects: a thermodynamic investigation of the calcium EDTA interaction

To better understand the origin of multivalency effects in ligand binding, the binding of a series of mono-, bi-, tri- and tetravalent carboxylate ligands to Ca(II) was examined by isothermal titration calorimetry (ITC). The data are inconsistent with an entropic origin of enhanced affinity, but rather show that at least in this instance the multivalency effect is enthalpic in origin. Analysis of binding data using the Jencks model shows the addition of incremental carboxylate "ligands" produces an unfavorable interaction entropy that is more than offset by a strongly favorable interaction enthalpy. The most likely source of this interaction enthalpy is the relief of repulsive Coulombic interactions in the unbound state. The conformational entropy penalty arising from the restriction of flexible dihedrals is negligible, within experimental error. On the other hand, an enthalpic contribution from linker restriction contributes strongly to the overall thermodynamics of ligand binding. Together, these data suggest that enthalpic effects dominate ligand binding, and design strategies should seek to optimize these interactions. The incorporation of unfavorable interactions in the unbound ligand that are relieved during binding provides an important mechanism by which to enhance ligand affinities.     

Epitope affinity chromatography and biophysical studies of monoclonal antibodies and recombinant antibody fragments

Peptide epitope affinity chromatography is a powerful technique for the purification of antibodies. This study aims to demonstrate the versatility of the technique and to show how biophysical techniques such as circular dichroism (CD) and fluorescence quenching (FQ) can aid the rational design of affinity ligands and characterization of antibody-based reagents. The performance of a number of peptide ligands for the purification of a range of different antibodies and recombinant fragments is investigated by automated fast-protein liquid chromatography. Purified products are analyzed for purity by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. They are then radiolabelled and the immunoreactivity is determined. Ligands are analyzed for secondary structural characteristics by CD and for binding affinity by FQ. Finally, a study is performed to investigate the thermal stability of a recombinant antibody fragment by CD analysis. It is found that simple ligand modifications such as the introduction of a C-terminal cysteine residue can improve purification performance. The FQ studies show that the modified peptide has a higher affinity for antibody. The CD analysis shows that it has a tendency to dimerize because of the formation of disulfide bonds. The versatility of epitope affinity is demonstrated through the purification of a recombinant diabody (dbFv) and by the use of a separate peptide matrix for the purification of an unrelated antibody. All studies result in antibody preparations of high purity and immunoreactivity. The CD analysis of the dbFv shows that it is denatured at 37 degrees C and is therefore unsuitable as a targeting reagent for use in humans in its present form. It is concluded that epitope affinity chromatography coupled with biophysical analyses plays an important role in the production and characterization of antibody-based reagents for targeted diagnosis and therapy of human diseases.     

Development and screening of epoxy‐spacer‐phage cryogels for affinity chromatography: Enhancing the binding capacity

Macroporous epoxy cryogels can be used as an alternative for classical matrices in affinity chromatography. Due to the structural properties of cryogels, with pores of up to 100 μm, crude samples can be processed at high speed without previous manipulations such as clarification or centrifugation. Also, we previously used a peptide‐expressing M13 bacteriophage as an affinity ligand. These ligands show high specificity toward the target to be purified. Combination of both, leads to a relative cost‐effective one‐step chromatographic set‐up delivering a high purity sample (>95%), however, so far with limited capacity. To increase the binding capacity of the affinity columns, we now inserted spacers between the chromatographic matrix and the phage ligand. Both linear spacers, di‐amino‐alkanes (C₂–C₁₀), and branched polyethyleneimine spacers with different molecular weights (800 Da–10 kDa) were analyzed. Two types of peptide expressing phage ligands, a linear 15‐mer and a cyclic 6‐mer, were used for screening. Up to a tenfold increase in binding capacity was observed depending on the combination of phage ligand and spacer type.     

Impact of planarity of unfused aromatic molecules on G-quadruplex binding: learning from isaindigotone derivatives

G-quadruplex structures are a new class of attractive targets for DNA-interactive anticancer agents. The primary building block of this structure is the G-quartet, which is composed of four coplanar guanines and serves as the major binding site for small molecules. NMR studies and molecular dynamics simulations have suggested that the planarity of G-quartet surface has been highly dynamic in solution. To better investigate how the planarity of unfused aromatic ligand impacts on its quadruplex binding properties, a variety of planarity controllable isaindigotone derivatives were designed and synthesized. The interaction of G-quadruplex DNA with these designed ligands was systematically explored using a series of biophysical studies. The FRET-melting, SPR, and CD spectroscopy results showed that reducing the planarity of their unfused aromatic core resulted in their decreased binding affinity and stabilization ability for G-quadruplex. NMR studies also suggested that these compounds could stack on the G-quartet surface. Such results are in parallel with subsequent molecular modeling studies. A detailed binding energy analysis indicated that van der Waals energy (ΔE(vdw)) and entropy (TΔS) are responsible for their decreased quadruplex binding and stabilization effect. All these results provided insight information about how quadruplex recognition could be controlled by adjusting the planarity of ligands, which shed light on further development of unfused aromatic molecules as optimal G-quadruplex binding ligands.     

Conformational stabilization of an engineered binding protein

We analyzed the thermodynamic basis for improvement of a binding protein by disulfide engineering. The Z(SPA)(-)(1) affibody binds to its Z domain binding partner with a dissociation constant K(d) = 1.6 microM, and previous analyses suggested that the moderate affinity is due to the conformational heterogeneity of free Z(SPA)(-)(1) rather than to a suboptimal binding interface. Studies of five stabilized Z(SPA)(-)(1) double cystein mutants show that it is possible to improve the affinity by an order of magnitude to K(d) = 130 nM, which is close to the range (20 to 70 nM) observed with natural Z domain binders, without altering the protein-protein interface obtained by phage display. Analysis of the binding thermodynamics reveals a balance between conformational entropy and desolvation entropy: the expected and favorable reduction of conformational entropy in the best-binding Z(SPA)(-)(1) mutant is completely compensated by an unfavorable loss of desolvation entropy. This is consistent with a restriction of possible conformations in the disulfide-containing mutant and a reduction of average water-exposed nonpolar surface area in the free state, resulting in a smaller conformational entropy penalty, but also a smaller change in surface area, for binding of mutant compared to wild-type Z(SPA)(-)(1). Instead, higher Z domain binding affinity in a group of eight Z(SPA)(-)(1) variants correlates with more favorable binding enthalpy and enthalpy-entropy compensation. These results suggest that protein-protein binding affinity can be improved by stabilizing conformations in which enthalpic effects can be fully explored.     

Microcalorimetric studies of the mechanism of interaction between designed peptides and hydrophobic adsorbents

To understand the mechanism of interaction between peptides and peptides with hydrophobic ligands, the oligomers (GWG, GWWG, GWWWG) were designed and synthesized to study adsorption behavior with octyl sepharose and CM-octyl sepharose. By batch equilibrium binding analysis and dilution heat of peptide solution measurement, the binding isotherm and adsorption enthalpy were obtained and the binding thermodynamics parameters were calculated and analyzed. In the isotherm analysis, we reveled that the affinity of GWG for both adsorbents is stronger than that of GWWG and GWWWG. The results demonstrate that the cation-pi interaction between the peptides and the buffer molecules is significant for solutions of peptides with tryptophan residues, and the solvation is competitive with the hydrophobic interaction between the peptides and the hydrophobic ligands. From the dilution heat measurements, we observed an endothermic dilution heat for GWG and exothermic for GWWG and GWWWG. All these results indicate that the increased tryptophan chain length can promote the solvation behavior of the peptides by the peptide-buffer interaction in this buffer system. Comparing the types of ligands reveals that the binding affinities of each peptide for the two adsorbents are similar. However, the mechanism of adsorption for peptides with hydrophobic ligands might be quite different with respect to the binding enthalpy between peptides and adsorbents. The adsorption of the peptides on octyl sepharose is an entropy-driven process for all the peptides. In contrast, the adsorption of CM-octyl sepharose with GWG and GWWG is an enthalpy-driven process, whereas that with GWWWG is entropy-driven. These findings indicate that the amount of tryptophan controls the characteristics of the peptides and the interaction mechanism in the binding procedure. This study of the adsorption mechanism of the designed peptide could provide fundamental information for peptide purification and amino acid residue behavior in peptide drug design.     

Multiple-injection affinity capillary electrophoresis to examine binding constants between glycopeptide antibiotics and peptides

Multiple-injection affinity capillary electrophoresis (MIACE) was used to determine binding constants (K(b)) between vancomycin, ristocetin, and teicoplanin from Streptomyces orientalis, Nocardia lurida, and Actinoplanes teichomyceticus, respectively, and fluorenylmethoxycarbonyl (Fmoc)-(Gly, Ala, Val, and Phe)-D-Ala-D-Ala peptides. In this technique, separate plugs of sample containing non-interacting standards, peptide one, buffer, and peptide two, were injected into the capillary column and electrophoresed. Peptides migrate through the column at similar electrophoretic mobilities but remain as distinct zones due to the buffer plug between peptides. The electrophoresis is then carried out in an increasing concentration of antibiotic in the running buffer. Continued electrophoresis results in a shift in the migration time of the peptides upon binding to the antibiotic. Analysis of the change in the relative migration time ratio (RMTR) of the resultant complexes relative to the non-interacting standards, as a function of the concentration of antibiotic yields a value for K(b). MIACE is a versatile technique that can be used to measure affinity constants between ligands of similar relative molecular mass and charge without the need of separate binding experiments. The findings described, herein, demonstrate the advantages of using MIACE to estimate binding parameters between ligands and receptors.     

Synthesis and binding studies of Alzheimer ligands on solid support

Aminopyrazole derivatives constitute the first class of nonpeptidic rationally designed beta-sheet ligands. Here we describe a double solid-phase protocol for both synthesis and affinity testing. The presented solid-phase synthesis of four types of hybrid compounds relies on the Fmoc strategy and circumvents subsequent HPLC purification by precipitating the final product from organic solution in pure form. Hexa- and octapeptide pendants with internal di- and tetrapeptide bridges are now amenable in high yields to combinatorial synthesis of compound libraries for high-throughput screening purposes. Solid-phase peptide synthesis (SPPS) on an acid-resistant PAM allows us, after PMB deprotection, to subject the free aminopyrazole binding sites in an immobilized state to on-bead assays with fluorescence-labeled peptides. From the fluorescence emission intensity decrease, individual binding constants can be calculated via reference curves by simple application of the law of mass action. Gratifyingly, host/guest complexation can be monitored quantitatively even for those ligands, which are almost insoluble in water.     

Using physics-based pose predictions and free energy perturbation calculations to predict binding poses and relative binding affinities for FXR ligands in the D3R Grand Challenge 2

Computer-aided drug design has become an integral part of drug discovery and development in the pharmaceutical and biotechnology industry, and is nowadays extensively used in the lead identification and lead optimization phases. The drug design data resource (D3R) organizes challenges against blinded experimental data to prospectively test computational methodologies as an opportunity for improved methods and algorithms to emerge. We participated in Grand Challenge 2 to predict the crystallographic poses of 36 Farnesoid X Receptor (FXR)-bound ligands and the relative binding affinities for two designated subsets of 18 and 15 FXR-bound ligands. Here, we present our methodology for pose and affinity predictions and its evaluation after the release of the experimental data. For predicting the crystallographic poses, we used docking and physics-based pose prediction methods guided by the binding poses of native ligands. For FXR ligands with known chemotypes in the PDB, we accurately predicted their binding modes, while for those with unknown chemotypes the predictions were more challenging. Our group ranked #1st (based on the median RMSD) out of 46 groups, which submitted complete entries for the binding pose prediction challenge. For the relative binding affinity prediction challenge, we performed free energy perturbation (FEP) calculations coupled with molecular dynamics (MD) simulations. FEP/MD calculations displayed a high success rate in identifying compounds with better or worse binding affinity than the reference (parent) compound. Our studies suggest that when ligands with chemical precedent are available in the literature, binding pose predictions using docking and physics-based methods are reliable; however, predictions are challenging for ligands with completely unknown chemotypes. We also show that FEP/MD calculations hold predictive value and can nowadays be used in a high throughput mode in a lead optimization project provided that crystal structures of sufficiently high quality are available.     

Multimerization of an Apoptogenic TRAIL‐Mimicking Peptide by Using Adamantane‐Based Dendrons

We have developed a straightforward strategy to multimerize an apoptogenic peptide that mimics the natural tumor necrosis factor‐related apoptosis‐inducing ligand (TRAIL) by using adamantane‐based dendrons as multivalent scaffolds. The selective binding affinity of the ligands to TRAIL receptor 2 (TR2) was studied by surface plasmon resonance, thus demonstrating that the trimeric and hexameric forms of the peptide exert an increased affinity of about 1500‐ and 20 000‐fold, respectively, relative to the monomer. Moreover, only the trimeric and hexameric ligands were able to induce cell death in TR2 expressing cells (BJAB), thus confirming that a multivalent form of the peptide is necessary to trigger a substantial TR2‐dependent apoptotic response in vitro. These results provide interesting insight into the multivalency effect on biological ligand/receptor interactions for future therapeutic applications.