Fluoroalkyl and alkyl chains have similar hydrophobicities in binding to the "hydrophobic wall" of carbonic anhydrase

The hydrophobic effect, the free-energetically favorable association of nonpolar solutes in water, makes a dominant contribution to binding of many systems of ligands and proteins. The objective of this study was to examine the hydrophobic effect in biomolecular recognition using two chemically different but structurally similar hydrophobic groups, aliphatic hydrocarbons and aliphatic fluorocarbons, and to determine whether the hydrophobicity of the two groups could be distinguished by thermodynamic and biostructural analysis. This paper uses isothermal titration calorimetry (ITC) to examine the thermodynamics of binding of benzenesulfonamides substituted in the para position with alkyl and fluoroalkyl chains (H(2)NSO(2)C(6)H(4)-CONHCH(2)(CX(2))(n)CX(3), n = 0-4, X = H, F) to human carbonic anhydrase II (HCA II). Both alkyl and fluoroalkyl substituents contribute favorably to the enthalpy and the entropy of binding; these contributions increase as the length of chain of the hydrophobic substituent increases. Crystallography of the protein-ligand complexes indicates that the benzenesulfonamide groups of all ligands examined bind with similar geometry, that the tail groups associate with the hydrophobic wall of HCA II (which is made up of the side chains of residues Phe131, Val135, Pro202, and Leu204), and that the structure of the protein is indistinguishable for all but one of the complexes (the longest member of the fluoroalkyl series). Analysis of the thermodynamics of binding as a function of structure is compatible with the hypothesis that hydrophobic binding of both alkyl and fluoroalkyl chains to hydrophobic surface of carbonic anhydrase is due primarily to the release of nonoptimally hydrogen-bonded water molecules that hydrate the binding cavity (including the hydrophobic wall) of HCA II and to the release of water molecules that surround the hydrophobic chain of the ligands. This study defines the balance of enthalpic and entropic contributions to the hydrophobic effect in this representative system of protein and ligand: hydrophobic interactions, here, seem to comprise approximately equal contributions from enthalpy (plausibly from strengthening networks of hydrogen bonds among molecules of water) and entropy (from release of water from configurationally restricted positions).     

Insights into Modeling Approaches in Chemistry: Assessing Ligand-Protein Binding Thermodynamics Based on Rigid-Flexible Model Molecules

In the field of chemistry, model compounds find extensive use for investigating complex objects. One prime example of such object is the protein-ligand supramolecular interaction. Prediction the enthalpic and entropic contribution to the free energy associated with this process, as well as the structural and dynamic characteristics of protein-ligand complexes poses considerable challenges. This review exemplifies modeling approaches used to study protein-ligand binding (PLB) thermodynamics by employing pairs of conformationally constrained/flexible model molecules. Strategically designing the model molecules can reduce the number of variables that influence thermodynamic parameters. This enables scientists to gain deeper insights into the enthalpy and entropy of PLB, which is relevant for medicinal chemistry and drug design. The model studies reviewed here demonstrate that rigidifying ligands may induce compensating changes in the enthalpy and entropy of binding. Some “rules of thumb” have started to emerge on how to minimize entropy-enthalpy compensation and design efficient rigidified or flexible ligands.     

Structural and ITC Characterization of Peptide-Protein Binding: Thermodynamic Consequences of Cyclization Constraints,a Case Study on Vascular Endothelial Growth Factor Ligands

Macrocyclization constraints are widely used in the design of protein ligands to stabilize their bioactive conformation and increase their affinities. However, the resulting changes in binding entropy can be puzzling and uncorrelated to affinity gains. Here, the thermodynamic (Isothermal Titration Calorimetry) and structural (X-ray, NMR and CD) analysis of a complete series of lactam-bridged peptide ligands of the vascular endothelial growth factor, and their unconstrained analogs are reported. It is shown that differences in thermodynamics arise mainly from the folding energy of the peptide upon binding. The systematic reduction in conformational entropy penalty due to helix pre-organization can be counterbalanced by an unfavorable vibrational entropy change if the constraints are too rigid. The gain in configurational entropy partially escapes the enthalpy/entropy compensation and leads to an improvement in affinity. The precision of the analytical ITC method makes this study a possible benchmark for constrained peptides optimization.     

Thermodynamic Parameters for the Association of Fluorinated Benzenesulfonamides with Bovine Carbonic Anhydrase II

This paper describes a calorimetric study of the association of a series of seven fluorinated benzenesulfonamide ligands (C₆H_nF_5?nSO₂NH₂) with bovine carbonic anhydrase II (BCA). Quantitative structure–activity relationships between the free energy, enthalpy, and entropy of binding and pK_a and log P of the ligands allowed the evaluation of the thermodynamic parameters in terms of the two independent effects of fluorination on the ligand: its electrostatic potential and its hydrophobicity. The parameters were partitioned to the three different structural interactions between the ligand and BCA: the Zn^II cofactor–sulfonamide bond (≈65 % of the free energy of binding), the hydrogen bonds between the ligand and BCA (≈10 %), and the contacts between the phenyl ring of the ligand and BCA (≈25 %). Calorimetry revealed that all of the ligands studied bind in a 1:1 stoichiometry with BCA; this result was confirmed by ¹⁹F NMR spectroscopy and X‐ray crystallography (for complexes with human carbonic anhydrase II).     

An enthalpic component in cooperativity: the relationship between enthalpy, entropy, and noncovalent structure in weak associations.

Attempts to quantify binding interactions of noncovalent complexes in aqueous solution have been stymied by complications arising from enthalpy-entropy compensation and cooperativity. We have extended work detailing the relationship between noncovalent structure and free energy of binding to include the roles of enthalpy and entropy of association. On the basis of van't Hoff measurements of the dimerization of vancomycin type antibiotics, we demonstrate that positive cooperativity manifests itself in a more favorable enthalpy of association and a partially compensating less favorable entropy of association. Finally, we extend these results to rationalize thermodynamic observations in unrelated systems.     

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.     

Residual ligand entropy in the binding of p-substituted benzenesulfonamide ligands to bovine carbonic anhydrase II

In studies on the thermodynamics of ligand-protein interactions, it is often assumed that the configurational and conformational entropy of the ligand is zero in the bound state (i.e., the ligand is rigidly fixed in the binding pocket). However, there is little direct experimental evidence for this assumption, and in the case of binding of p-substituted benzenesulfonamide inhibitors to bovine carbonic anhydrase II (BCA II), the observed thermodynamic binding signature derived from isothermal titration calorimetry experiments leads indirectly to the conclusion that a considerable degree of residual entropy remains in the bound ligand. Specifically, the entropy of binding increases with glycine chain length n, and strong evidence exists that this thermodynamic signature is not driven by solvent reorganization. By use of heteronuclear (15)N NMR relaxation measurements in a series (n = 1-6) of (15)N-glycine-enriched ligands, we find that the observed thermodynamic binding signature cannot be explained by residual ligand dynamics in the bound state, but rather results from the indirect influence of ligand chain length on protein dynamics.     

A free energy calculation study of the effect of H-->F substitution on binding affinity in ligand-antibody interactions

Changes in binding affinity to catalytic antibody 6D9 of chloramphenicol phosphonate derivatives (CPDs) containing H or F were investigated by performing free energy calculations based on molecular dynamics simulations. We calculated the binding free energy, enthalpy, and entropy changes (DeltaDeltaG, DeltaDeltaH, and -TDeltaDeltaS) attributable to H-->F substitution by comparing results for CPDs containing a trifluoroacetylamino group (CPD-F) or an acetylamino group (CPD-H). The calculated DeltaDeltaG, DeltaDeltaH, and -TDeltaDeltaS values were -2.9, -6.3, and 3.5 kcal mol(-1) and close to experimental values observed for a series of similar ligands, chloramphenicol phosphonates with F and H (-1.4, -3.5, and 2.1 kcal mol(-1)). Therefore, CPD-F binds more strongly to 6D9 than does CPD-H. To clarify the origin of the large difference in DeltaDeltaG, we apportioned the calculated values of DeltaDeltaG and DeltaG for the associated and dissociated states into contributions from various atomic interactions. We found that the H-->F substitution increased the binding affinity mainly by decreasing the hydration free energy and not by increasing favorable interactions with the antibody. The decreased hydration free energy of the ligand was mainly due to unfavorable coulombic interactions between the trifluoroacetylamino group and solvent waters, which increased the free energy of the dissociated state (by about 3.7 kcal mol(-1)). Also, the trifluoroacetylamino group slightly increased the free energy level of the associated state (about 0.8 kcal mol(-1)) because favorable van der Waals interactions compensated for unfavorable coulombic interactions with antibody atoms. In addition, the enthalpy and entropy changes, DeltaDeltaH and -TDeltaDeltaS (computationally -6.3 and 3.5 kcal mol(-1)), originated mainly from a decrease in hydration free energy in the dissociated state. The CPD-F and CPD-H ligands had substantially different structures in the dissociated and complexed states.     

Involvement of water in carbohydrate-protein binding.

The interactions of trimannosides 1 and 2 with Con A were studied to reveal the effects of displacement of well-ordered water molecules on the thermodynamic parameters of protein-ligand complexation. Trisaccharide 2 is a derivative of 1, in which the hydroxyl at C-2 of the central mannose unit is replaced by a hydroxyethyl moiety. Upon binding, this moiety displaces a conserved water molecule present in the Con A binding site. Structural studies by NMR spectroscopy and MD simulations showed that the two compounds have very similar solution conformational properties. MD simulations of the complexes of Con A with 1 and 2 demonstrated that the hydroxyethyl side chain of 2 can establish the same hydrogen bonds in a low energy conformation with the protein binding site as those mediated by the water molecule in the complex of 1 with Con A. Isothermal titration microcalorimetry (ITC) measurements showed that 2 has a more favorable entropy of binding compared to 1. This term, which was expected, arises from the return of the highly ordered water molecule to bulk solution. The favorable entropy term was, however, offset by a relatively large unfavorable enthalpy term. This observation was rationalized by comparing the extent of hydrogen bond and solvation changes during binding. It is proposed that an indirect interaction through a water molecule will provide a larger number of hydrogen bonds in the complex that have higher occupancies than in bulk solution, thereby stabilizing the complex.     

Protein-ligand interactions: thermodynamic effects associated with increasing nonpolar surface area

Thermodynamic parameters were determined for complex formation between the Grb2 SH2 domain and Ac-pTyr-Xaa-Asn derived tripeptides in which the Xaa residue is an α,α-cycloaliphatic amino acid that varies in ring size from three- to seven-membered. Although the six- and seven-membered ring analogs are approximately equipotent, binding affinities of those having three- to six-membered rings increase incrementally with ring size because increasingly more favorable binding enthalpies dominate increasingly less favorable binding entropies, a finding consistent with an enthalpy-driven hydrophobic effect. Crystallographic analysis reveals that the only significant differences in structures of the complexes are in the number of van der Waals contacts between the domain and the methylene groups in the Xaa residues. There is a positive correlation between buried nonpolar surface area and binding free energy and enthalpy, but not with ΔC(p). Displacing a water molecule from a protein-ligand interface is not necessarily reflected in a favorable change in binding entropy. These findings highlight some of the fallibilities associated with commonly held views of relationships of structure and energetics in protein-ligand interactions and have significant implications for ligand design.     

非离子表面活性剂和β-环状糊精包结作用的微量热法研究

用Model-4200型恒温滴定量热系统,研究了非离子表面活性剂(TritonX-114,TritonX-100,TritonX-102,TritonX-405)和β-环状糊精的包结的热力学性质,结果表明:对所有体系,均形成1:1的包结物,ΔHⅲ和ΔSⅲ值均为负值,形成包结物的过程是焓有利而熵不利的过程,而且-ΔHⅲ和ΔSⅲ值随非离子表面活性剂环氧乙烷(EO)链长的增加而减少。另外,非离子表面活性剂环氧乙烷(EO)链越长,包结物稳定常数越小。     

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.     

Computer simulation of the distribution of hexane in a lipid bilayer: spatially resolved free energy, entropy, and enthalpy profiles

The partitioning behavior of small molecules in lipid bilayers is important in a variety of areas including membrane protein folding and pharmacology. However, the inhomogeneous nature of lipid bilayers on a nanometer length scale complicates experimental studies of membrane partitioning. To gain more insight in the partitioning of a small molecule into the lipid bilayer, we have carried out atomistic computer simulations of hexane in a dioleoyl phosphatidylcholine model membrane. We have been able to obtain spatially resolved free energy, entropy, enthalpy, and heat capacity profiles based on umbrella sampling calculations at three different temperatures. In agreement with experiment, hexane partitions preferentially to the center of the bilayer. This process is driven almost entirely by a favorable entropy change, consistent with the hydrophobic effect. In contrast, partitioning to the densest region of the acyl chains is dominated by a favorable enthalpy change with a small entropy change, which is consistent with the "nonclassical" hydrophobic effect or "bilayer" effect. We explain the features of the entropy and enthalpy profiles in terms of density and free volume in the system.     

Relative contributions of desolvation, inter- and intramolecular interactions to binding affinity in protein kinase systems

In several previous studies, we performed sensitivity analysis to gauge the relative importance of different atomic partial charges in determining protein-ligand binding. In this work, we gain further insights by decomposing these results into three contributions: desolvation, intramolecular interactions, and intermolecular interactions, again based on a Poisson continuum electrostatics model. Three protein kinase-inhibitor systems have been analyzed: CDK2-deschloroflavopiridol, PKA-PKI, and LCK-PP2. Although our results point out the importance of specific intermolecular interactions to the binding affinity, they also reveal the remarkable contributions from the solvent-mediated intramolecular interactions in some cases. Thus, it is necessary to look beyond analyzing protein-ligand interactions to understand protein-ligand recognition or to gain insights into designing ligands and proteins. In analyzing the contributions of the three components to the overall binding free energy, the PKA-PKI system with a much larger ligand was found to behave differently from the other two systems with smaller ligands. In the former case, the intermolecular interactions are very favorable, and together with the favorable solvent-mediated intramolecular interactions, they overcome the large desolvation penalties to give a favorable electrostatics contribution to the overall binding affinity. On the other hand, the other two systems with smaller ligands only present modest intermolecular interactions and they are not or are only barely sufficient to overcome the desolvation penalty even with the aid of the favorable intramolecular contributions. As a result, the binding affinity of these two systems do not or only barely benefit from electrostatics contributions.     

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.     

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.     

Binding of methyl orange and its homologs by powdered nylon 66

Powdered Nylon 66 was prepared as a model of amorphous polymers. The resultant powder polyamide was composed of only amorphous regions. The extent of uptake of the acid azo dyes, a homologous series of methyl orange derivatives, by the polymer was measured in an aqueous solution. The first binding constants and the thermodynamic parameters in the course of the binding were evaluated. The thermodynamic behaviors obtained are very similar to those of crosslinked polyvinylpyrrolidone. The favorable free energy of the binding is accompanied by an entropy gain and an exothermic enthalpy change. The shorter the alkyl chain of the dyes, the more negative is the enthalpy change and, hence, the smaller is the entropy change. The thermodynamic data for butyl orange showed that the binding process is athermal and is wholly an entropic effect. The binding of the dyes to the matrix is entropically favorable as a result of the operation of the hydrophobic effect. In addition, an electrostatic force is operative between the sulfonate group on the dyes and the terminal amino groups on the polyamide.     

Strong solute-solute dispersive interactions in a protein-ligand complex

The contributions of solute-solute dispersion interactions to binding thermodynamics have generally been thought to be small, due to the surmised equality between solute-solvent dispersion interactions prior to the interaction versus solute-solute dispersion interactions following the interaction. The thermodynamics of binding of primary alcohols to the major urinary protein (MUP-I) indicate that this general assumption is not justified. The enthalpy of binding becomes more favorable with increasing chain length, whereas the entropy of binding becomes less favorable, both parameters showing a linear dependence. Despite the hydrophobicity of the interacting species, these data show that binding is not dominated by the classical hydrophobic effect, but can be attributed to favorable ligand-protein dispersion interactions.     

Amperometric method for determining the degree of complexation of polyelectrolytes with cationic surfactants

The complexation of sodium polystyrene sulfonate with monovalent cationic surfactants at a microsized liquid/liquid interface has been studied using electrochemistry. The method is based on measurement of surfactant ion transfer across the interface between two immiscible electrolyte solutions (ITIES). The complexation of various cationic surfactants (alkylpyridinium- and trimethylammonium-) with oligosized polystyrene sulfonate was measured. Binding isotherms were used to determine the degree of binding as a function of the surfactant chain length and type of head group. It was found that the hydrophobicity of the surfactant was the predominant factor. The effect of the polyelectrolyte chain length on the binding mechanism was studied using cetylpyridinium chloride as a complexing agent. It was found that binding affinity, as well as cooperativity of the binding process, decreases with decreasing polyelectrolyte chain length. Thermodynamics of surfactant binding was measured using titration microcalorimetry. The thermodynamic data obtained show that the enthalpy of surfactant binding is not dependent on polymer chain length, but an increase in chain length makes the binding process entropically more favorable.     

人血清白蛋白与季铵盐双子表面活性剂的相互作用

在298.15 K下, 应用等温滴定量热法研究了人血清白蛋白(HSA)与两种季铵盐双子表面活性剂[(C_nN)₂Cl₂, n＝12, 14]在缓冲溶液(pH＝7.0)中相互作用的热力学性质. 实验结果表明, HSA对这两种表面活性剂有两类结合位点, 分别为结合时需要吸收热量的强结合位点和可放出热量的弱结合位点. 两种表面活性剂对应的第一类结合——强结合为熵驱动过程, 且该结合位点对应的结合位点数、结合常数和热力学参数差别不大. 至于第二类结合——弱结合位点, 由于 (C₁₄N)₂Cl₂ 疏水链过长, 只有部分进入HSA的疏水空腔内, 因此相应的的结合位点数和放热量减小, 而熵变增加, 为焓和熵共同驱动的反应. 圆二色研究表明(C_nN)₂Cl₂的加入使HSA的二级结构发生变化, 这说明(C_nN)₂Cl₂与HSA的相互作用既包含结合反应也包含(C_nN)₂Cl₂诱导该蛋白部分结构改变的过程.