Enthalpy-entropy compensation and cooperativity as thermodynamic epiphenomena of structural flexibility in ligand-receptor interactions

Extreme enthalpy‒entropy compensation in the dimerization of small solutes in aqueous solution

This communication summarizes findings from the earliest encounters with extreme enthalpy‒entropy compensation, a phenomenon first detected in the 1950s by a reappraisal of isopiestic and calorimetric measurements on aqueous urea solutions in terms of solute self-association. Because concurrent studies of carboxylic acid association were confined to measurement of the equilibrium constant by conductance, IR spectrophotometry or potentiometric titration measurements, temperature-independence of the dimerization constant was mistakenly taken to signify a value of zero for Δ H o instead of (Δ H o  ‒ TΔ S o ). In those studies of small-solute self-association the extreme enthalpy‒entropy compensation was reflecting the action of water as a reactant whose hydroxyl groups were competing for the solute carbonyl involved in self-association. Such action gives rise to a positive temperature dependence of Δ H o that could well be operating in concert with that responsible for the commonly observed negative dependence for protein‒ligand interactions exhibiting extreme enthalpy‒entropy compensation, where the solvent contribution to the energetics reflects changes in the extent of ordered water structure in hydrophobic environments.

Assembly path dependence of telomeric DNA compaction by TRF1, TIN2, and SA1

Telomeres, complexes of DNA and proteins, protect ends of linear chromosomes. In humans, the two shelterin proteins TRF1 and TIN2, along with cohesin subunit SA1, were proposed to mediate telomere cohesion. Although the ability of the TRF1-TIN2 and TRF1-SA1 systems to compact telomeric DNA by DNA-DNA bridging has been reported, the function of the full ternary TRF1-TIN2-SA1 system has not been explored in detail. Here, we quantify the compaction of nanochannel-stretched DNA by the ternary system, as well as its constituents, and obtain estimates of the relative impact of its constituents and their interactions. We find that TRF1, TIN2, and SA1 work synergistically to cause a compaction of the DNA substrate, and that maximal compaction occurs if all three proteins are present. By altering the sequence with which DNA substrates are exposed to proteins, we establish that compaction by TRF1 and TIN2 can proceed through binding of TRF1 to DNA, followed by compaction as TIN2 recognizes the previously bound TRF1. We further establish that SA1 alone can also lead to a compaction, and that compaction in a combined system of all three proteins can be understood as an additive effect of TRF1-TIN2 and SA1-based compaction. Atomic force microscopy of intermolecular aggregation confirms that a combination of TRF1, TIN2, and SA1 together drive strong intermolecular aggregation as it would be required during chromosome cohesion.

High-resolution structures with bound Mn2+ and Cd2+ map the metal import pathway in an Nramp transporter

Transporters of the Nramp (Natural resistance-associated macrophage protein) family import divalent transition metal ions into cells of most organisms. By supporting metal homeostasis, Nramps prevent diseases and disorders related to metal insufficiency or overload. Previous studies revealed that Nramps take on a LeuT fold and identified the metal-binding site. We present high-resolution structures of Deinococcus radiodurans (Dra)Nramp in three stable conformations of the transport cycle revealing that global conformational changes are supported by distinct coordination geometries of its physiological substrate, Mn2+, across conformations, and by conserved networks of polar residues lining the inner and outer gates. In addition, a high-resolution Cd2+-bound structure highlights differences in how Cd2+ and Mn2+ are coordinated by DraNramp. Complementary metal binding studies using isothermal titration calorimetry with a series of mutated DraNramp proteins indicate that the thermodynamic landscape for binding and transporting physiological metals like Mn2+ is different and more robust to perturbation than for transporting the toxic Cd2+ metal. Overall, the affinity measurements and high-resolution structural information on metal substrate binding provide a foundation for understanding the substrate selectivity of essential metal ion transporters like Nramps.

Thermodynamic analysis of Zα domain-nucleic acid interactions

DNA/RNA molecules adopting the left-handed conformation (Z-form) have been attributed with immunogenic properties. However, their biological role and importance has been a topic of debate for many years. The discovery of Z-DNA/RNA binding domains (Zα domains) in varied proteins that are involved in the innate immune response, such as the interferon inducible form of the RNA editing enzyme ADAR1 (p150), Z-DNA binding protein 1 (ZBP1), the fish kinase PKZ and the poxvirus inhibitor of interferon response E3L, indicates important roles of Z-DNA/RNA in immunity and self/non-self-discrimination. Such Zα domain-containing proteins recognise left-handed Z-DNA/RNA in a conformation-specific manner. Recent studies have implicated these domains in virus recognition. Given these important emerging roles for the Zα domains, it is pivotal to understand the mechanism of recognition of the Z-DNA/Z-RNA by these domains. To this end, we assessed the binding thermodynamics of Zα domain from ORF112 and ADAR1 on T(CG)3 and T(CG)6 oligonucleotides which have high propensity to adopt the Z-conformation. Our study highlights important differences in the mode of oligonucleotide binding by the two Zα domains originating from different proteins. Site-directed mutagenesis was employed together with isothermal titration calorimetry to tease apart finer details of the binding thermodynamics. Our work advances the understanding on binding thermodynamics of Zα domains to their cognate nucleic acid substrates and paves the ground for future efforts to gain a complete appreciation of this process.

Binding Specificity of ASHH2 CW Domain Toward H3K4me1 Ligand Is Coupled to Its Structural Stability Through Its α1-Helix

The CW domain binds to histone tail modifications found in different protein families involved in epigenetic regulation and chromatin remodeling. CW domains recognize the methylation state of the fourth lysine on histone 3 and could, therefore, be viewed as a reader of epigenetic information. The specificity toward different methylation states such as me1, me2, or me3 depends on the particular CW subtype. For example, the CW domain of ASHH2 methyltransferase binds preferentially to H3K4me1, and MORC3 binds to both H3K4me2 and me3 modifications, while ZCWPW1 is more specific to H3K4me3. The structural basis for these preferential bindings is not well understood, and recent research suggests that a more complete picture will emerge if dynamical and energetic assessments are included in the analysis of interactions. This study uses fold assessment by NMR in combination with mutagenesis, ITC affinity measurements, and thermal denaturation studies to investigate possible couplings between ASHH2 CW selectivity toward H3K4me1 and the stabilization of the domain and loops implicated in binding. The key elements of the binding site—the two tryptophans and the α1-helix form and maintain the binding pocket— were perturbed by mutagenesis and investigated. Results show that the α1-helix maintains the overall stability of the fold via the I915 and L919 residues and that the correct binding consolidates the loops designated as η1 and η3, as well as the C-terminal. This consolidation is incomplete for H3K4me3 binding to CW, which experiences a decrease in overall thermal stability on binding. Loop mutations not directly involved in the binding site, nonetheless, affect the equilibrium positions of the key residues.

An Unrevealed Molecular Function of Corannulene Buckybowl Glycoconjugates in Selective Tumor Annihilation by Targeting the Cancer-Specific Warburg Effect

The biomedical application of corannulene π-bowls is historically limited by low solubility and bioavailability despite the potential in their unique electronic properties for new functional materials. Herein, the unexpected role and molecular mechanism of Corranulene π-bowls are uncovered in biomedical applications as an effective anticancer agent for Warburg effect mediated selective tumor targeting. The corannulene triazolyl monosaccharides Cor-sugars exhibit highly potent cytotoxicity against human cancer cells and effectively inhibit xenograft growth of hyperglycolytic tumors. Particularly, the galactose-conjugated Cor-gal exhibits superior in vivo anticancer efficacy in A549 tumor models with outstanding safety profile compared to doxorubicin. Moreover, the combined treatment of Cor-gal with immune checkpoint inhibitor results in an effective synergy in treating H460 human lung carcinoma. An uptake mechanism study reveals that Cor-sugars exploit tumor-specific glucose transporter glucose transporter 1 (GLUT1) for targeted cell delivery and intra-tumoral accumulation through the cancer-specific Warburg effect. Their significant anticancer activity is attributed to multiphasic DNA-binding and cell cycle alteration effects. This study uncovers new molecular properties of corannulene buckybowl and enabling their potential new applications in biomedical engineering.

Entropy Effects in Intermolecular Associations of Crown-Ethers and Cyclodextrins with Amino Acids in Aqueous and in Non-Aqueous Media

The analysis of the ratios of entropy and enthalpy characteristics and their contributions to the change in the Gibbs energy of intermolecular interactions of crown ethers and cyclodextrins with amino acids is carried out. Two different types of macrocycles were chosen for examination: crown ethers with a hydrophilic interior and cyclodextrins with a hydrophobic inner cavity and a hydrophilic exterior. The thermodynamics of complex formation of crown ethers and cyclodextrins with amino acids in water and aqueous-organic solvents of variable composition was examined. The contributions of the entropy solvation of complexes of 18-crown-6 with glycine, alanine, phenylalanine to the change in the entropy of complexation in water-ethanol and water-dimethyl sulfoxide solvents was calculated and analyzed. It was found that the ratios of the entropy and enthalpy solvation of the reagents for these systems have similar trends when moving from water to aqueous-organic mixtures. The relationship between the thermodynamic characteristics and structural features of the complexation processes between cyclodextrins and amino acids has been established. The thermodynamic enthalpy–entropy compensation effect was revealed, and its features for complexation of cyclodextrins and 18-crown-6 were considered. It was concluded that, based on the thermodynamic parameters of molecular complexation, one could judge the mode of the formation of complexes, the main driving forces of the interactions, and the degree of desolvation.

Structure, Function, and Thermodynamics of Lactate Dehydrogenases from Humans and the Malaria Parasite P. falciparum

Temperature adaptation is ubiquitous among all living organisms, yet the molecular basis for this process remains poorly understood. It can be assumed that for parasite-host systems, the same enzymes found in both organisms respond to the same selection factor (human body temperature) with similar structural changes. Herein, we report the existence of a reversible temperature-dependent structural transition for the glycolytic enzyme lactate dehydrogenase (LDH) from the malaria parasite Plasmodium falciparum (pfLDH) and human heart (hhLDH) occurring in the temperature range of human fever. This transition is observed for LDHs from psychrophiles, mesophiles, and moderate thermophiles in their operating temperature range. Thermodynamic analysis reveals unique thermodynamic signatures of the LDH-substrate complexes defining a specific temperature range to which human LDH is adapted and parasite LDH is not, despite their common mesophilic nature. The results of spectroscopic analysis combined with the available crystallographic data reveal the existence of an active center within pfLDH that imparts psychrophilic structural properties to the enzyme. This center consists of two pockets, one formed by the five amino acids (5AA insert) within the substrate specificity loop and the other by the active site, that mutually regulate one another in response to temperature and induce structural and functional changes in the Michaelis complex. Our findings pave the way toward a new strategy for malaria treatments and drug design using therapeutic agents that inactivate malarial LDH selectively at a specific temperature range of the cyclic malaria paroxysm.

MHCBI: a pipeline for calculating peptide-MHC binding energy using semi-empirical quantum mechanical methods with explicit/implicit solvent models

Experimentally estimating peptide-major histocompatibility complex (pMHC) binding affinity has been quite challenging due to the many receptors and the many potential ligands implicated in it. We have thus proposed a straightforward computational methodology considering the different mechanisms involved in pMHC binding to facilitate studying such receptor-ligand interactions. We have developed a pipeline using semi-empirical quantum mechanical methods for calculating pMHC class I and II molecules' binding energy (BE). This pipeline can systematize the methodology for calculating pMHC system BE, enabling the rational design of T-cell epitopes to be used as pharmaceuticals and vaccines.

Role of potassium levels in pkBADH heterogeneity of NAD+ binding site

Betaine aldehyde dehydrogenase (BADH) catalyzes the oxidation of betaine aldehyde to glycine betaine using NAD⁺ as a coenzyme. Studies in porcine kidney BADH (pkBADH) suggested that the enzyme exhibits heterogeneity of active sites and undergoes potassium-induced conformational changes. This study aimed to analyze if potassium concentration plays a role in the heterogeneity of pkBADH active sites through changes in NAD⁺ affinity constants, in its secondary structure content and stability. The enzyme was titrated with NAD⁺ 1 mM at fixed-variable KCl concentration, and the interaction measured by Isothermal Titration Calorimetry (ITC) and Circular Dichroism (CD). ITC data showed that K⁺ increased the first active site affinity in a manner dependent on its concentration; K_D values to the first site were 14.4, 13.1, and 10.4 μM, at 25, 50, and 75 mM KCl. ΔG values showed that the coenzyme binding is a spontaneous reaction without changes between active sites or depending on KCl concentration. ΔH and TΔS_b values showed that NAD⁺ binding to the active site is an endothermic process and is carried out at the expense of changes in entropy. α-Helix content increased as KCl increased, enzyme (T_m)_app values were 2.6 °C and 3.3 °C higher at 20 mM and 200 mM K⁺. PkBADH molecular model showed three different interaction K⁺ sites. Results suggested K⁺ can interact with pkBADH and cause changes in the secondary structure, it provokes changes in the enzyme affinity by the coenzyme, and in the thermostability.

Thermophoresis of biological and biocompatible compounds in aqueous solution

With rising popularity of microscale thermophoresis for the characterisation of protein-ligand binding reactions and possible applications in microfluidic devices, there is a growing interest in considering thermodiffusion in the context of life sciences. But although the understanding of thermodiffusion in non-polar mixtures has grown rapidly in recent years, predictions for associated mixtures like aqueous solutions remain challenging. This review aims to give an overview of the literature on thermodiffusion in aqueous systems, show the difficulties in theoretical description that arise from the non-ideal behaviour of water-mixtures, and highlight the relevance of thermodiffusion in a biological context. We find that the thermodiffusion in aqueous systems is dominated by contributions from heat of transfer, hydrogen bond interactions and charge effects. However, the separation of these effects is often difficult, especially in case of biological systems where a systematic exclusion of contributions may not be feasible.

Lattice theory for binding of linear polymers to a solid substrate from polymer melts. II. Influence of van der Waals interactions and chain semiflexibility on molecular binding and adsorption

The polymeric Langmuir theory, developed in Paper I [J. Dudowicz et al., J. Chem. Phys. 151, 124706 (2019)], is employed to investigate the influence of van der Waals interactions and chain rigidity on the thermodynamics of the reversible molecular binding to interfaces in one component polymer fluids (polymer melts). Both van der Waals interactions and chain stiffness are found to influence the temperature variation of the surface coverage Θ, the binding transition itself, and the cooperativity of molecular binding. Re-entrancy of the surface coverage Θ(T) is found to arise when the intermolecular interactions are sufficiently attractive to cause a liquid-vapor like phase separation in the interfacial region, a phenomenon that can occur in the binding of both small molecules and polymer chains to surfaces.

The Enthalpy-entropy Compensation Phenomenon. Limitations for the Use of Some Basic Thermodynamic Equations

The thermodynamic analyses of proteins, protein-ligands and protein-nucleic acid complexes involves the entropy-enthalpy (S-H) compensation phenomenon. We have examined the question whether the observed compensation is artificial or reflects anything more than the well-known laws of statistical thermodynamics (so-called extra-thermodynamic compensation). We have shown that enthalpy- entropy compensation (EEC) is mainly the trivial consequence of the basic thermodynamic laws and there are no experimental evidences for existence of the extra-thermodynamic compensation. In most cases EEC obtained in the experiments through the plot enthalpies (ΔH) and entropies (TΔS) versus one another is meaningless due to the large correlated errors in ΔH and TΔS, unless special measures are taken to minimize, quantify and propagate these errors. Van't Hoff equation can be used for entropy calculation in limited cases when enthalpy is measured in independent experiments. Eyring equation cannot be used for calculation of entropy in any case and should be excluded from scientific use. Both equation, Van't Hoff and Eyring cannot be used for simultaneous calculation of the enthalpy and entropy values using one set of data. All the data obtained in this way should be recognized as erroneous.

Theoretical explanation for the pharmaceutical incompatibility through the cooperativity effect of the drug-drug intermolecular interactions in the phenobarbital∙∙∙paracetamol∙∙∙H2O complex

In order to reveal the essence of the pharmaceutical incompatibility, the cooperativity effects of the drug-drug intermolecular π∙∙∙π and H∙∙∙O H-bonding interactions involving hydration were evaluated in the phenobarbital∙∙∙paracetamol∙∙∙H₂O complex at the M06-2X/6-311++G** and MP2/6-311++G** levels. The thermodynamic cooperativity effects were also investigated by the statistical thermodynamic method. The results show that the π∙∙∙π stacking ternary complexes with the moderate anti-cooperativity effects are dominant in controling the aggregation process of phenobarbital, paracetamol, and H₂O, as is confirmed by the atoms-in-molecules (AIM) and reduced density gradient (RDG) analyses. Therefore, it can be inferred that the anti-cooperativity effect plays an important role in forming the pharmaceutical incompatibility, and thus a deduction on the formation process of the pharmaceutical incompatibility between phenobarbital and paracetamol, with the hydration effect, is given. Several valuable models that relate the features of molecular surface electrostatic potentials or their statistical parameters, such as the surface areas, average values ([Formula: see text]), variances ([Formula: see text], [Formula: see text] and [Formula: see text]), and product of [Formula: see text] and electrostatic balance parameter (ν) ([Formula: see text]ν), to the values of the cooperativity effects were predicted. The formation of the pharmaceutical incompatibility is a thermodynamic cooperativity process driven by the enthalpy change. Graphical abstract Anti-cooperativity effect plays an important role in forming the pharmaceutical incompatibility.

A Rational Study of the Origin and Generality of Anti-Enthalpy–Entropy Compensation (AEEC) Phenomenon

In addition to enthalpy–entropy compensation (EEC), anti-enthalpy–entropy compensation (AEEC) phenomenon is also found in literature. The reports on the latter are limited, and analyses and justifications are so far unclear. Herein we present demonstration on the nature and possibility of the AEEC phenomenon. Although literature reports so far have mostly shown linear AEEC, we have found both linear and non-linear dependences. The non-linearity we consider arises from large “free energy window” (FEW) like EEC, recently presented and discussed. Attempts have been made to rationalize the observations in terms of solvation–desolvation phenomenon of the involved processes anticipated in previous studies. We have found that enthalpy and entropy of formation of gaseous and solid materials may also exhibit AEEC, it can be thus classified as a global phenomenon. In addition to AEEC, the phenomenon of no-enthalpy–entropy compensation (NEEC) is also reported. Thus, the phenomena EEC, AEEC and NEEC are thermodynamic puzzles that require close attention and analysis. With the help of wide range of physical chemical processes, an elaborate general understanding of the linear and non-linear AEEC phenomena has been attempted. The experimental data herein used are collected from literature reports, new measurements were not done. A variety of examples have supported possible generality of AEEC.

Molecular rigidity and enthalpy-entropy compensation in DNA melting

Enthalpy-entropy compensation (EEC) is observed in diverse molecular binding processes of importance to living systems and manufacturing applications, but this widely occurring phenomenon is not sufficiently understood from a molecular physics standpoint. To gain insight into this fundamental problem, we focus on the melting of double-stranded DNA (dsDNA) since measurements exhibiting EEC are extensive for nucleic acid complexes and existing coarse-grained models of DNA allow us to explore the influence of changes in molecular parameters on the energetic parameters by using molecular dynamics simulations. Previous experimental and computational studies have indicated a correlation between EEC and changes in molecular rigidity in certain binding-unbinding processes, and, correspondingly, we estimate measures of DNA molecular rigidity under a wide range of conditions, along with resultant changes in the enthalpy and entropy of binding. In particular, we consider variations in dsDNA rigidity that arise from changes of intrinsic molecular rigidity such as varying the associative interaction strength between the DNA bases, the length of the DNA chains, and the bending stiffness of the individual DNA chains. We also consider extrinsic changes of molecular rigidity arising from the addition of polymer additives and geometrical confinement of DNA between parallel plates. All our computations confirm EEC and indicate that this phenomenon is indeed highly correlated with changes in molecular rigidity. However, two distinct patterns relating to how DNA rigidity influences the entropy of association emerge from our analysis. Increasing the intrinsic DNA rigidity increases the entropy of binding, but increases in molecular rigidity from external constraints decreases the entropy of binding. EEC arises in numerous synthetic and biological binding processes and we suggest that changes in molecular rigidity might provide a common origin of this ubiquitous phenomenon in the mutual binding and unbinding of complex molecules.

Cooperativity effect involving drug-DNA/RNA intermolecular interaction: A B3LYP-D3 and MP2 theoretical investigation on ketoprofen⋯cytosine⋯H2O system

In order to examine the origin of the drug action and design new DNA/RNA-targeted drugs, the cooperativity effect involving drug-DNA/RNA intermolecular interaction in ketoprofen⋯cytosine⋯H₂O ternary system were investigated by the B3LYP, B3LYP-D3, and MP2 methods with the 6-311++G(2d,p) basis set. The thermodynamic cooperativity was also evaluated at 310.15 K. The N-H⋯O, O-H⋯O, O-H⋯N, C-H⋯N, and C-H⋯O H bonds coexist in ternary complexes. The intermolecular interactions obtained by B3LYP-D3 are close to those calculated by MP2. The steric effects and van der Waals interactions have little influence on the cooperativity effects. The anti-cooperativity effect in ket⋯cyt⋯H₂O is far more notable than the cooperativity effect, and the stability of the cyclic structure with anti-cooperativity effect is higher than that of the linear structure with cooperativity effect, as is confirmed by the AIM (atoms in molecules) and RDG (reduced density gradient) analysis. Thus, it can be inferred that, in the presence of H₂O, the anti-cooperativity effect plays a dominant role in the drug-DNA/RNA interaction, and the nature of the hydration in the binding of drugs to DNA/RNA bases is the H-bonding anti-cooperativity effect. Furthermore, the drug always links simultaneously with DNA/RNA base and H₂O, and only in this way can the biological activity of drugs play a role. In most cases, the enthalpy change is the major factor driving the cooperativity, as is different from most of biomacromolecule complexes.

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.

Enhanced potency of bivalent small molecule gp41 inhibitors

Low molecular weight peptidomimetic inhibitors with hydrophobic pocket binding properties and moderate fusion inhibitory activity against HIV-1 gp41-mediated cell fusion were elaborated by increasing the available surface area for interacting with the heptad repeat-1 (HR1) coiled coil on gp41. Two types of modifications were tested: 1) increasing the overall hydrophobicity of the molecules with an extension that could interact in the HR1 groove, and 2) forming symmetrical dimers with two peptidomimetic motifs that could potentially interact simultaneously in two hydrophobic pockets on the HR1 trimer. The latter approach was more successful, yielding 40-60times improved potency against HIV fusion over the monomers. Biophysical characterization, including equilibrium binding studies by fluorescence and kinetic analysis by Surface Plasmon Resonance, revealed that inhibitor potency was better correlated to off-rates than to binding affinity. Binding and kinetic data could be fit to a model of bidentate interaction of dimers with the HR1 trimer as an explanation for the slow off-rate, albeit with minimal cooperativity due to the highly flexible ligand structures. The strong cooperativity observed in fusion inhibitory activity of the dimers implied accentuated potency due to the transient nature of the targeted intermediate. Optimization of monomer, dimer or higher order structures has the potential to lead to highly potent non-peptide fusion inhibitors by targeting multiple hydrophobic pockets.

Configurational space discretization and free energy calculation in complex molecular systems

We sought to design a free energy calculation scheme with the hope of saving cost for generating dynamical information that is inherent in trajectories. We demonstrated that snapshots in a converged trajectory set are associated with implicit conformers that have invariant statistical weight distribution (ISWD). Since infinite number of sets of implicit conformers with ISWD may be created through independent converged trajectory sets, we hypothesized that explicit conformers with ISWD may be constructed for complex molecular systems through systematic increase of conformer fineness, and tested the hypothesis in lipid molecule palmitoyloleoylphosphatidylcholine (POPC). Furthermore, when explicit conformers with ISWD were utilized as basic states to define conformational entropy, change of which between two given macrostates was found to be equivalent to change of free energy except a mere difference of a negative temperature factor, and change of enthalpy essentially cancels corresponding change of average intra-conformer entropy. By implicitly taking advantage of entropy enthalpy compensation and forgoing all dynamical information, constructing explicit conformers with ISWD and counting thermally accessible number of which for interested end macrostates is likely to be an efficient and reliable alternative end point free energy calculation strategy.

Introduction of a tailor made anion receptor into the side chain of small peptides allows fine-tuning the thermodynamic signature of peptide–DNA binding

Replacing lysine by a tailor made anion binding motif, switches the DNA binding of small peptides from entropy to enthalpy control.

Kinetic binding assays for the analysis of protein-ligand interactions

The importance of binding kinetics in terms of residence time and on-rate in drug discovery has been broadly accepted in the past few years. Furthermore, evidence has accumulated that the optimal binding mechanism of a drug to its target molecule is related to physiological efficacy as well as selectivity and thus drug safety. Homogeneous fluorescence-based binding assays have been shown to enable high throughput kinetics requiring only small amounts of protein and can be developed to elucidate even complex mechanisms of molecular recognition. A generalized approach is proposed that combines high quality kinetic and equilibrium data in an Integrated Global Fit analysis yielding the most probable binding mechanism.

The Thermodynamic Mechanism of Peptide-MHC Class II Complex Formation Is a Determinant of Susceptibility to HLA-DM

Peptides bind MHC class II molecules through a thermodynamically nonadditive process consequent to the flexibility of the reactants. Currently, how the specific outcome of this binding process affects the ensuing epitope selection needs resolution. Calorimetric assessment of binding thermodynamics for hemagglutinin 306-319 peptide variants to the human MHC class II HLA-DR1 (DR1) and a mutant DR1 reveals that peptide/DR1 complexes can be formed with different enthalpic and entropic contributions. Complexes formed with a smaller entropic penalty feature circular dichroism spectra consistent with a non-compact form, and molecular dynamics simulation shows a more flexible structure. The opposite binding mode, compact and less flexible, is associated with greater entropic penalty. These structural variations are associated with rearrangements of residues known to be involved in HLA-DR (DM) binding, affinity of DM for the complex, and complex susceptibility to DM-mediated peptide exchange. Thus, the thermodynamic mechanism of peptide binding to DR1 correlates with the structural rigidity of the complex, and DM mediates peptide exchange by "sensing" flexible complexes in which the aforementioned residues are rearranged at a higher frequency than in more rigid ones.

Quasithermodynamic contributions to the fluctuations of a protein nanopore

Proteins undergo thermally activated conformational fluctuations among two or more substates, but a quantitative inquiry on their kinetics is persistently challenged by numerous factors, including the complexity and dynamics of various interactions, along with the inability to detect functional substates within a resolvable time scale. Here, we analyzed in detail the current fluctuations of a monomeric β-barrel protein nanopore of known high-resolution X-ray crystal structure. We demonstrated that targeted perturbations of the protein nanopore system, in the form of loop-deletion mutagenesis, accompanying alterations of electrostatic interactions between long extracellular loops, produced modest changes of the differential activation free energies calculated at 25 °C, ΔΔG(⧧), in the range near the thermal energy but substantial and correlated modifications of the differential activation enthalpies, ΔΔH(⧧), and entropies, ΔΔS(⧧). This finding indicates that the local conformational reorganizations of the packing and flexibility of the fluctuating loops lining the central constriction of this protein nanopore were supplemented by changes in the single-channel kinetics. These changes were reflected in the enthalpy-entropy reconversions of the interactions between the loop partners with a compensating temperature, TC, of ∼300 K, and an activation free energy constant of ∼41 kJ/mol. We also determined that temperature has a much greater effect on the energetics of the equilibrium gating fluctuations of a protein nanopore than other environmental parameters, such as the ionic strength of the aqueous phase as well as the applied transmembrane potential, likely due to ample changes in the solvation activation enthalpies. There is no fundamental limitation for applying this approach to other complex, multistate membrane protein systems. Therefore, this methodology has major implications in the area of membrane protein design and dynamics, primarily by revealing a better quantitative assessment on the equilibrium transitions among multiple well-defined and functionally distinct substates of protein channels and pores.

Dynamical insights into the differential characteristics of Mycobacterium avium subsp. paratuberculosis peptide binding to HLA-DRB1 proteins associated with multiple sclerosis

Differential properties of MAP binding to HLA proteins in Sardinian MS patients.

Thermodynamics of binding of structurally similar ligands to histone deacetylase 8 sheds light on challenges in the rational design of potent and isozyme-selective inhibitors of the enzyme

Among the different histone deacetylase (HDAC) isozymes, HDAC8 is the most highly malleable enzyme, and it exhibits the potential to accommodate structurally diverse ligands (albeit with moderate binding affinities) in its active site pocket. To probe the molecular basis of this feature, we performed detailed thermodynamic studies of the binding of structurally similar ligands, which differed with respect to the “cap”, “linker”, and “metal-binding” regions of the suberoylanilide hydroxamic acid (SAHA) pharmacophore, to HDAC8. The experimental data revealed that although the enthalpic (ΔH°) and entropic (ΔS°) changes for the binding of individual SAHA analogues to HDAC8 were substantially different, their binding free energies (ΔG°) were markedly similar, conforming to a strong enthalpy–entropy compensation effect. This effect was further observed in the temperature-dependent thermodynamics of binding of all SAHA analogues to the enzyme. Notably, in contrast to other metalloenzymes, our isothermal titration calorimetry experiments (performed in different buffers of varying ionization enthalpies) suggest that depending on the ligand, its zinc-binding group may or may not be deprotonated upon the binding to HDAC8. Furthermore, the heat capacity changes (ΔC_p°) associated with the ligand binding to HDAC8 markedly differed from one SAHA analogue to the other, and such features could primarily be rationalized in light of the dynamic flexibility in the enzyme structure in conjunction with the reorganization of the active site resident water molecules. Arguments are presented that although the binding thermodynamic features described above would facilitate identification of weak to moderately tight-binding HDAC8 inhibitors (by a high-throughput and/or virtual screening of libraries of small molecules), they would pose major challenges for the structure-based rational design of highly potent and isozyme-selective inhibitors of human HDAC8.

A B3LYP and MP2(full) theoretical investigation into cooperativity effects, aromaticity and thermodynamic properties in the Na(+)⋯benzonitrile⋯H2O ternary complex

The cooperativity effects between H-bonding and Na(+)⋯π or Na(+)⋯σ interactions in Na(+)⋯benzonitrile⋯H2O complexes were investigated using the B3LYP and MP2(full) methods with 6-311++G(2d,p) and aug-cc-pVTZ basis sets. The thermodynamic cooperativity and the influence of this cooperativity on aromaticity was evaluated by nucleus-independent chemical shifts (NICS). The results showed that the influence of the Na(+)⋯σ or Na(+)⋯π interaction on the hydrogen bond is more pronounced than that of the latter on the former. The cooperativity effect appeared in the Na(+)⋯σ interaction complex while the anti-cooperativity effect tended to be in the Na(+)⋯π system. The change in enthalpy is the major factor driving cooperativity. Thermodynamic cooperativity is not in accordance with the cooperativity effect evaluated by the change of interaction energy. The ring aromaticity of is weakened while the bond dissociation energy (BDE) of the C-CN bond increases upon ternary complex formation. The cooperativity effect (E coop) correlates with R c (NICS(1)ternary/NICS(1)binary) and ΔΔδ (Δδ ternary - Δδ binary) involving the ring and C ≡ N bond, as well as R BDE(C-CN) [BDE(C-CN)ternary/BDE(C-CN)binary], respectively. AIM (atoms in molecules) analysis confirms the existence of cooperativity.

Antigenic peptide molecular recognition by the DRB1–DQB1 haplotype modulates multiple sclerosis susceptibility

DRB1–DQB1 binding affinities in peptide molecular recognition process. (A) In protective haplotype DRB1 allele displays a preferential affinity for MBP peptide, while (B) in predisposing haplotype DQB1 allele displays a preferential affinity for EBNA1 peptide.

Thermodynamics of Peptide-MHC Class II Interactions: Not all Complexes are Created Equal

The adaptive immune response begins when CD4+ T cells recognize antigenic peptides bound to class II molecules of the Major Histocompatibility Complex (MHCII). The interaction between peptides and MHCII has been historically interpreted as a rigid docking event. However, this model has been challenged by the evidence that conformational flexibility plays an important role in peptide-MHCII complex formation. Thermodynamic analysis of the binding reaction suggests a model of complexation in which the physical-chemical nature of the peptide determines the variability in flexibility of the substates in the peptide-MHC conformational ensemble. This review discusses our understanding of the correlation between thermodynamics of peptide binding and structural features of the resulting complex as well as their impact on HLA-DM activity and on our ability to predict MHCII-restricted epitopes.

For many but not for all: how the conformational flexibility of the peptide/MHCII complex shapes epitope selection

The adaptive immune response starts when CD4+ T cells recognize peptide antigens presented by class II molecules of the Major Histocompatibility Complex (MHCII). Two outstanding features of MHCII molecules are their polymorphism and the ability of each allele to bind a large panoply of peptides. The ability of each MHCII molecule to interact with a limited, though broad, range of amino acid sequences, or "permissive specificity" of binding, is the result of structural flexibility. This flexibility has been identified through biochemical and biophysical studies, and molecular dynamic simulations have modeled the conformational rearrangements that the peptide and the MHCII undergo during interaction. Moreover, there is evidence that the structural flexibility of the peptide/MHCII complex correlates with the activity of the "peptide-editing" molecule DM. In light of the impact that these recent findings have on our ability to predict MHCII epitopes, a review of the structural and thermodynamic determinants of peptide binding to MHCII is proposed.

Entropy-enthalpy compensation: role and ramifications in biomolecular ligand recognition and design

Recent calorimetric studies of interactions between small molecules and biomolecular targets have generated renewed interest in the phenomenon of entropy-enthalpy compensation. In these studies, entropic and enthalpic contributions to binding are observed to vary substantially and in an opposing manner as the ligand or protein is modified, whereas the binding free energy varies little. In severe examples, engineered enthalpic gains can lead to completely compensating entropic penalties, frustrating ligand design. Here, we examine the evidence for compensation, as well as its potential origins, prevalence, severity, and ramifications for ligand engineering. We find the evidence for severe compensation to be weak in light of the large magnitude of and correlation between errors in experimental measurements of entropic and enthalpic contributions to binding, though a limited form of compensation may be common. Given the difficulty of predicting or measuring entropic and enthalpic changes to useful precision, or using this information in design, we recommend ligand engineering efforts instead focus on computational and experimental methodologies to directly assess changes in binding free energy.

Entropy-enthalpy Compensation of Biomolecular Systems in Aqueous Phase: a Dry Perspective

We survey thermodynamic measurements on processes involving biological macromolecules in aqueous solution, which illustrate well the ubiquitous phenomenon of entropy-enthalpy compensation. The processes include protein folding/unfolding and ligand binding/unbinding, with compensation temperatures varying by about 50 K around an average near 293 K. We show that incorporating both near-exact entropy-enthalpy compensation (due to solvent relaxation) and multi-excitation entropy (from vibrational quanta) leads to a compensation temperature in water of about 230 K. We illustrate a general procedure for subtracting solvent and environment-related terms to determine the bare Gibbs free energy changes of chemical processes.