Conformational Entropy Change Research Articles

Electrochemical Redox Refrigeration

The high conformational entropy change of the Fe(CN)63−/4− redox reaction can be used as the basis for a compact electrochemical refrigerator. This device is comparable to a liquid version of a Peltier cooler, with two distinct advantages: (1) the entropy change per carrier (1.5 mV/K) of the electrochemical refrigerant is more than 5 times larger than that of state-of-the-art solid thermoelectric materials; and (2) the liquid electrolyte can be advected continuously away from the cooling junction, so that Joule heating in the bulk element does not diminish the delivered cooling effect. In this work, we use infrared microscopy to visualize the thermal aspects of Fe(CN)63−/4− redox, and compare the estimated cooling to calculated values with and without electrolyte flow. While the temperature differences achieved in a single cell are small (~50 mK) and not enhanced by electrolyte flow, the cooling power density (~0.5 W/cm3) is high when normalized to the small electrode volume. Non-dimensional figures of merit are proposed to identify electrochemical redox species for maximizing the cooling effect.

The State of Proteins not with Standing, Translational Velocity is Vital for their Function

Background: Targeting of macromolecules by means of translational diffusion for therapeutic reason is generally of interest. But translational velocity that enables directional delivery of any molecule on target is given less attention.
 Objectives: The objectives of this research are to: 1) determine the value of the cohesion factor affecting solutes, 2) determine the translational velocity of porcine pancreatic (PPAA) - and human salivary (HSAA) – alpha amylase and ions, 3) rederive the effective kinetic energy (K.E.) of solutes, 4) determine the thermodynamic parameters for a folded to an unfolded transition and 5) give reasons why the velocity of solution components is generally very important.
 Methods: A theoretical research and experimentation using Bernfeld method.
 Results and Discussion: The K.E. of solution components as re-derived is « 3kBT / 2 (where kB and T are Boltzmann constant and Kelvin temperature respectively.). The velocities of hydrolysis of the substrate with the sucrose-treated PPAA were generally higher than those of HSAA. The values of conformational entropy change (∆Sconf) for PPAA were generally higher than those of HSAA. Expectedly the values of ∆Sconf were positive.
 Conclusion: In conclusion, the square root of the cohesion factor is larger than 22.4 exp (+3) 310.15 K / 273.15 K / 18, accounting for the translational velocity (u) in solution being « gas phase velocity (U). The translational diffusion D and u remain respectively, a function of the hydrodynamic radius of the solutes in particular and the magnitude of D; unfolding of proteins decreases the values of the parameters. Overall, unfolding is entropy driven. Without the mobility of solution components at desired velocity directional delivery of small molecules to site of need such as intrinsically disordered proteins etc may remain impossible.

New Insight on the Interface between Polythiophene and Semiconductors via Molecular Dynamics Simulations.

Polythiophene is considered as an effective dry adhesive and is promising to be a conductive adhesive due to its excellent properties. Here, we used steered molecular dynamics to investigate the interfacial strength between polythiophene and various semiconductors with similar structures including silicon, silicon carbide, and diamond. Energy decomposition was done to have a detailed insight into the adhesive mechanism. Particularly, we laid stress on the entropy difference of the polythiophene chain in different systems. Van der Waals interaction and electrostatic interaction both positively contributed to the adhesion between polythiophene and semiconductors, while the entropy change of polythiophene, including vibrational entropy change and conformational entropy change, weakened the adhesion to some extent. Our results indicated that the combined effect of these three factors made the adhesion between polythiophene and silicon carbide the strongest among the systems we studied. Additionally, it was found that such adhesion was scarcely influenced by temperature. This simple polythiophene-semiconductor interfacial study can help optimize the choice of the semiconductor when applying the polythiophene adhesive.

Role of Conformational Entropy in Molecular Recognition by TAZ1 of CBP.

The globular transcriptional adapter zinc binding 1 (TAZ1) domain of CREB binding protein participates in protein-protein interactions that are involved in transcriptional regulation. TAZ1 binds numerous targets, of which many are intrinsically disordered proteins that undergo a disorder-to-order transition to various degrees. One such target is the disordered transactivation domain of transcription factor RelA (TAD-RelA), which with its interaction with TAZ1 is involved in transcriptional regulation of genes in NF-κB signaling. We have here performed nuclear magnetic resonance backbone and side-chain relaxation studies to investigate the influence of RelA-TA2 (residues 425-490 in TAD-RelA) binding on the subnanosecond internal motions of TAZ1. We find a considerable dynamic response on both the backbone and side-chain levels, which corresponds to a conformational entropy change that contributes significantly to the binding energetics. We further show that the microscopic origins of the dynamic response of TAZ1 vary depending on the target. This study demonstrates that folded protein domains that are able to interact with various targets are not dynamically passive but can have a significant role in the motional response upon target association.

Structural and Thermodynamic Signatures of Ligand Binding to the Enigmatic Chitinase D of Serratia proteamaculans.

The Gram-negative bacteria Serratia marcescens and Serratia proteamaculans have efficient chitinolytic machineries that degrade chitin into N-acetylglucosamine (GlcNAc), which is used as a carbon and energy source. The enzymatic degradation of chitin in these bacteria occurs through the synergistic action of glycoside hydrolases (GHs) that have complementary activities; an endo-acting GH (ChiC) making random scissions on the polysaccharide chains and two exo-acting GHs mainly targeting single reducing (ChiA) and nonreducing (ChiB) chain ends. Both bacteria produce low amounts of a fourth GH18 (ChiD) with an unclear role in chitin degradation. Here, we have determined the thermodynamic signatures for binding of (GlcNAc)6 and the inhibitor allosamidin to SpChiD as well as the crystal structure of SpChiD in complex with allosamidin. The binding free energies for the two ligands are similar (Δ Gr° = -8.9 ± 0.1 and -8.4 ± 0.1 kcal/mol, respectively) with clear enthalpic penalties (Δ Hr° = 3.2 ± 0.1 and 1.8 ± 0.1 kcal/mol, respectively). Binding of (GlcNAc)6 is dominated by solvation entropy change (- TΔ Ssolv° = -17.4 ± 0.4 kcal/mol) and the conformational entropy change dominates for allosamidin binding (- TΔ Sconf° = -9.0 ± 0.2 kcal/mol). These signatures as well as the interactions with allosamidin are very similar to those of SmChiB suggesting that both enzymes are nonreducing end-specific.

Resistance mechanism of the oncogenic β3-αC deletion mutation in BRAF kinase to dabrafenib and vemurafenib revealed by molecular dynamics simulations and binding free energy calculations.

BRAF kinase is an essential target for anti-cancer drug development. Emergence of the β3-αC loop deletion mutation (ΔNVTAP) in BRAF kinase frequently occurred in human cancers seriously compromises the therapeutic efficacy of some BRAF kinase inhibitors, such as dabrafenib and vemurafenib. However, the mechanism of this resistance is still not well understood. In this study, the influence of the β3-αC deletion mutation on the binding profiles of three BRAF kinase inhibitors (AZ628, dabrafenib, and vemurafenib) with BRAFV600E or BRAFΔNVTAP was explored by conventional molecular dynamics (MD) simulations and binding free energy calculations. The simulation results indicated that the β3-αC deletion mutation enhances the flexibility of the αC helix and alters their conformations, which amplify the conformational entropy change (-TΔS) and weaken the interactions between the inhibitors and BRAF. The further per-residue binding free energy decomposition analysis revealed that the ΔNVTAP mutation changed the contributions of a few key residues to the bindings of dabrafenib or vemurafenib, such as L57, L66, W83, C84, F135, G145, and F147, but did not have obvious impact on the contributions of these residues to AZ628. Our results provide valuable clues to understand the mechanisms of drug resistance conferred by the β3-αC deletion mutation.

Scaling of Rates of Vibrational Energy Transfer in Proteins with Equilibrium Dynamics and Entropy.

Theoretical arguments and results of molecular dynamics (MD) simulations of myoglobin at 300 K are presented to relate rates of vibrational energy transfer across nonbonded contacts interacting via short-range potentials to dynamics of the contact. Both theory and the results of the simulations support a scaling relation between the energy transfer rate and the inverse of the variance in the distance between hydrogen-bonded contacts. The results of the MD simulations do not support such a relation for longer-range charged contacts. Instead, the energy transfer rate is found to scale as a power law in the distance between charged groups. The scaling between rates of vibrational energy transfer across nonbonded contacts interacting via short-range potentials and conformational dynamics suggests a relation between vibrational energy transfer rates and entropy associated with the dynamics of interacting residues. The use of time-resolved vibrational spectroscopy to determine change in conformational entropy with change in protein functional state is discussed, and an expression quantifying the connection is provided.

Tuning Free Energy by Backbone Conformational Entropy: A Strategy from Disordered Proteins

Tuning Free Energy by Backbone Conformational Entropy: A Strategy from Disordered Proteins

Structural and thermodynamic insights into β-1,2-glucooligosaccharide capture by a solute-binding protein in Listeria innocua

β-1,2-Glucans are bacterial carbohydrates that exist in cyclic or linear forms and play an important role in infections and symbioses involving Gram-negative bacteria. Although several β-1,2-glucan-associated enzymes have been characterized, little is known about how β-1,2-glucan and its shorter oligosaccharides (Sop n s) are captured and imported into the bacterial cell. Here, we report the biochemical and structural characteristics of the Sop n -binding protein (SO-BP, Lin1841) associated with the ATP-binding cassette (ABC) transporter from the Gram-positive bacterium Listeria innocua Calorimetric analysis revealed that SO-BP specifically binds to Sop n s with a degree of polymerization of 3 or more, with Kd values in the micromolar range. The crystal structures of SO-BP in an unliganded open form and in closed complexes with tri-, tetra-, and pentaoligosaccharides (Sop3-5) were determined to a maximum resolution of 1.6 Å. The binding site displayed shape complementarity to Sop n , which adopted a zigzag conformation. We noted that water-mediated hydrogen bonds and stacking interactions play a pivotal role in the recognition of Sop3-5 by SO-BP, consistent with its binding thermodynamics. Computational free-energy calculations and a mutational analysis confirmed that interactions with the third glucose moiety of Sop n s are significantly responsible for ligand binding. A reduction in unfavorable changes in binding entropy that were in proportion to the lengths of the Sop n s was explained by conformational entropy changes. Phylogenetic and sequence analyses indicated that SO-BP ABC transporter homologs, glycoside hydrolases, and other related proteins are co-localized in the genomes of several bacteria. This study may improve our understanding of bacterial β-1,2-glucan metabolism and promote the discovery of unidentified β-1,2-glucan-associated proteins.

Impact of a non-synonymous Q281R polymorphism on structure of human Lipoprotein-Associated Phospholipase A2 (Lp-PLA2 ).

Non-synonymous single nucleotide polymorphisms (nsSNPs) are genetic variations at single base resulting in an amino acid change which have been associated with various complex human diseases. The human Lipoprotein-associated phospholipase A2 (Lp-PLA2 ) gene harbours a rare Q281R polymorphism which was previously reported to cause loss of enzymatic function. Lp-PLA2 is an important enzyme which catalyzes the hydrolysis of polar phospholipids releasing pro-atherogenic and pro-inflammatory mediators involved in the pathogenesis of atherosclerosis. Our current study is aimed at elucidating the structural and functional consequences of Q281R polymorphism on Lp-PLA2 . The Q281R mutation is classified as deleterious and causes protein instability as deduced from evolutionary, folding free energy changes and Support vector machine (SVM)-based methods. A Q281R mutant structure was deciphered using homology modelling approach and was validated using phi and psi dihedral angles distribution, ERRAT, Verify_3D scores, Protein Structure Analysis (ProSA) energ,y and Z-score. A decreased hydrophobic interactions and weaker substrate binding affinity was observed in the mutant compared to the wild- type (WT) using molecular docking. Further, the mutant displayed enhanced structural flexibility particularly in the low density lipoprotein (LDL) binding domain, decreased solvent accessibility of catalytic residues-Phe274 and Ser273 and increased Cɑ distance between Phe274 and Leu153 and large conformational entropy change as inferred from all-atom molecular dynamics (MD) simulation and essential dynamics (ED) studies. Our results corroborate well with previous experimental studies and thus these aberrations in the Q281R mutant structure may help explain the molecular basis of loss of enzyme activity.

Dynamics, Conformational Entropy, and Frustration in Protein-Protein Interactions Involving an Intrinsically Disordered Protein Domain.

Intrinsically disordered proteins (IDPs) are abundant in the eukaryotic proteome. However, little is known about the role of subnanosecond dynamics and the conformational entropy that it represents in protein-protein interactions involving IDPs. Using nuclear magnetic resonance side chain and backbone relaxation, stopped-flow kinetics, isothermal titration calorimetry, and computational studies, we have characterized the interaction between the globular TAZ1 domain of the CREB binding protein and the intrinsically disordered transactivation domain of STAT2 (TAD-STAT2). We show that the TAZ1/TAD-STAT2 complex retains considerable subnanosecond motions, with TAD-STAT2 undergoing only a partial disorder-to-order transition. We report here the first experimental determination of the conformational entropy change for both binding partners in an IDP binding interaction and find that the total change even exceeds in magnitude the binding enthalpy and is comparable to the contribution from the hydrophobic effect, demonstrating its importance in the binding energetics. Furthermore, we show that the conformational entropy change for TAZ1 is also instrumental in maintaining a biologically meaningful binding affinity. Strikingly, a spatial clustering of very high amplitude motions and a cluster of more rigid sites in the complex exist, which through computational studies we found to overlap with regions that experience energetic frustration and are less frustrated, respectively. Thus, the residual dynamics in the bound state could be necessary for faster dissociation, which is important for proteins that interact with multiple binding partners.

Understanding the resistance mechanism of penicillin binding protein 1a mutant against cefotaxime using molecular dynamic simulation

Antibiotic resistance is a threatening challenge for global health, as the expansion of resistance to current antibiotics has made serious therapeutic problems. Genome mutations are key evolutionary mechanisms conferring antibiotic resistance in bacterial pathogens. For example, penicillin and cephalosporins resistance is mostly mediated by mutations in penicillin binding proteins to change the affinity of the drug. Accordingly, threonine point mutations were reported to develop antibiotic resistance in various bacterial infections including pneumococcal infections. In this study, conventional molecular dynamics simulations, umbrella sampling simulations and MM/GBSA free energy calculations were applied to figure out how the Threonine to Alanine mutation (T to A) at STMK motif affects the binding of cefotaxime to Penicillin Binding Protein 1a and to reveal the resistance mechanism induced by the T to A mutation. The results obtained from the computational methods demonstrate that the T to A mutation increases the flexibility of the binding pocket and changes its conformation, which leads to increased conformational entropy change (−TΔS) and attenuates the bonds between the ligand and the receptor. In brief, our findings indicate that both of the alterations of the conformational enthalpy and entropy contribute to the T to A-induced resistance in the binding of cefotaxime into penicillin binding protein 1a.

Conformational Entropy of FK506 Binding to FKBP12 Determined by Nuclear Magnetic Resonance Relaxation and Molecular Dynamics Simulations.

FKBP12 (FK506 binding protein 12 kDa) is an important drug target. Nuclear magnetic resonance (NMR) order parameters, describing amplitudes of motion on the pico- to nanosecond time scale, can provide estimates of changes in conformational entropy upon ligand binding. Here we report backbone and methyl-axis order parameters of the apo and FK506-bound forms of FKBP12, based on 15N and 2H NMR relaxation. Binding of FK506 to FKBP12 results in localized changes in order parameters, notably for the backbone of residues E54 and I56 and the side chains of I56, I90, and I91, all positioned in the binding site. The order parameters increase slightly upon FK506 binding, indicating an unfavorable entropic contribution to binding of TΔ S = -18 ± 2 kJ/mol at 293 K. Molecular dynamics simulations indicate a change in conformational entropy, associated with all dihedral angles, of TΔ S = -26 ± 9 kJ/mol. Both these values are significant compared to the total entropy of binding determined by isothermal titration calorimetry and referenced to a reactant concentration of 1 mM ( TΔ S = -29 ± 1 kJ/mol). Our results reveal subtle differences in the response to ligand binding compared to that of the previously studied rapamycin-FKBP12 complex, despite the high degree of structural homology between the two complexes and their nearly identical ligand-FKBP12 interactions. These results highlight the delicate dependence of protein dynamics on drug interactions, which goes beyond the view provided by static structures, and reinforce the notion that protein conformational entropy can make important contributions to the free energy of ligand binding.

Conformational Entropy from Slowly Relaxing Local Structure Analysis of 15N-H Relaxation in Proteins: Application to Pheromone Binding to MUP-I in the 283-308 K Temperature Range.

The slowly relaxing local structure (SRLS) approach is applied to 15N-H relaxation from the major urinary protein I (MUP-I), and its complex with pheromone 2-sec-butyl-4,5-dihydrothiazol. The objective is to elucidate dynamics, and binding-induced changes in conformational entropy. Experimental data acquired previously in the 283-308 K temperature range are used. The N-H bond is found to reorient globally with correlation time, τ1,0, and locally with correlation time, τ2,0, where τ1,0 ≫ τ2,0. The local motion is restricted by the potential u = -c02D002, where D002 is the Wigner rotation matrix element for L = 2, K = 0, and c02 evaluates the strength of the potential. u yields straightforwardly the order parameter, ⟨D002⟩, and the conformational entropy, Sk, both given by Peq = exp(-u). The deviation of the local ordering/local diffusion axis from the N-H bond, given by the angle β, is also determined. We find that c02 ≅ 18 ± 4 and τ2,0 = 0-170 ps for ligand-free MUP-I, whereas c02 ≅ 15 ± 4 and τ2,0 = 20-270 ps for ligand-bound MUP-I. β is in the 0-10° range. c02 and τ2,0 decrease, whereas β increases, when the temperature is increased from 283 to 308 K. Thus, SRLS provides physically well-defined structure-related (c02 and ⟨D002⟩), motion-related (τ2,0), geometry-related (β), and binding-related (Sk) local parameters, and their temperature-dependences. Intriguingly, upon pheromone binding the conformational entropy of MUP-I decreases at high temperature and increases at low temperature. The very same experimental data were analyzed previously with the model-free (MF) method which yielded "global" (in this context, "relating to the entire 283-308 K range") amplitude (S2) and rate (τe) of the local motion, and a phenomenological exchange term (Rex). S2 is found to decrease (implying implicitly "global" increase in Sk) upon pheromone binding.

How Does the L884P Mutation Confer Resistance to Type-II Inhibitors of JAK2 Kinase: A Comprehensive Molecular Modeling Study

Janus kinase 2 (JAK2) has been regarded as an essential target for the treatment of myeloproliferative neoplasms (MPNs). BBT594 and CHZ868, Type-II inhibitors of JAK2, illustrate satisfactory efficacy in preclinical MPNs and acute lymphoblastic leukemia (ALL) models. However, the L884P mutation of JAK2 abrogates the suppressive effects of BBT594 and CHZ868. In this study, conventional molecular dynamics (MD) simulations, umbrella sampling (US) simulations and MM/GBSA free energy calculations were employed to explore how the L884P mutation affects the binding of BBT594 and CHZ868 to JAK2 and uncover the resistance mechanism induced by the L884P mutation. The results provided by the US and MD simulations illustrate that the L884P mutation enhances the flexibility of the allosteric pocket and alters their conformations, which amplify the conformational entropy change (−TΔS) and weaken the interactions between the inhibitors and target. Additionally, the structural analyses of BBT594 and CHZ868 in complex with the WT JAK2 illustrate that the drug tail with strong electronegativity and small size located in the allosteric pocket of JAK2 may enhance anti-resistance capability. In summary, our results highlight that both of the changes of the conformational entropies and enthalpies contribute to the L884P-induced resistance in the binding of two Type-II inhibitors into JAK2 kinase.

Entropy in molecular recognition by proteins

Molecular recognition by proteins is fundamental to molecular biology. Dissection of the thermodynamic energy terms governing protein-ligand interactions has proven difficult, with determination of entropic contributions being particularly elusive. NMR relaxation measurements have suggested that changes in protein conformational entropy can be quantitatively obtained through a dynamical proxy, but the generality of this relationship has not been shown. Twenty-eight protein-ligand complexes are used to show a quantitative relationship between measures of fast side-chain motion and the underlying conformational entropy. We find that the contribution of conformational entropy can range from favorable to unfavorable, which demonstrates the potential of this thermodynamic variable to modulate protein-ligand interactions. For about one-quarter of these complexes, the absence of conformational entropy would render the resulting affinity biologically meaningless. The dynamical proxy for conformational entropy or "entropy meter" also allows for refinement of the contributions of solvent entropy and the loss in rotational-translational entropy accompanying formation of high-affinity complexes. Furthermore, structure-based application of the approach can also provide insight into long-lived specific water-protein interactions that escape the generic treatments of solvent entropy based simply on changes in accessible surface area. These results provide a comprehensive and unified view of the general role of entropy in high-affinity molecular recognition by proteins.

15N-H-Related Conformational Entropy Changes Entailed By Plexin-B1 RBD Dimerization: Combined Molecular Dynamics/NMR Relaxation Approach.

We report on a new method for determining function-related conformational entropy changes in proteins. Plexin-B1 RBD dimerization serves as example, and internally mobile N-H bonds serve as probes. Sk (entropy in units of kBT) is given by -∫(PeqlnPeq)dΩ, where Peq = exp(-u) is the probability density for probe orientation, and u the local potential. Previous slowly relaxing local structure (SRLS) analyses of 15N-H relaxation in proteins determined linear combinations of D002(Ω) and (D022(Ω) + D0-22(Ω)) (D0KL(Ω) represents a Wigner rotation matrix element in uniaxial local medium) as "best-fit" form of u. SRLS also determined the "best-fit" orientation of the related ordering tensor. On the basis of this information the coefficients (in the linear combination) of the terms specified above are determined with molecular dynamics (MD) simulations. With the explicit expression for u thus in hand, Sk is calculated. We find that in general Sk decreases, i.e., the local order increases, upon plexin-B1 RBD dimerization. The largest decrease in Sk occurs in the helices α1 and α2, followed by the α2/β6 turn. Only the relatively small peripheral β2 strand, β2/α1 turn, and L3 loop become more disordered. That α-helices dominate ΔSk = Sk(dimer) - Sk(monomer), a few peripheral outliers partly counterbalance the overall decrease in Sk, and the probability density function, Peq, has rhombic symmetry given that the underlying potential function, u, has rhombic symmetry, are interesting features. We also derive S2 (the proxy of u in the simple "model-free (MF)" limit of SRLS) with MD. Its conversion into a potential requires assumptions and adopting a simple axial form of u. Ensuing ΔSk(MF) profiles are u-dependent and differ from ΔSk(SRLS). A method that provides consistent, general, and accurate Sk, atomistic/mesoscopic in nature, has been developed. Its ability to provide new insights in protein research has been illustrated.

Thermodynamics of tunnel formation upon substrate binding in a processive glycoside hydrolase

Glycoside hydrolases (GHs) catalyze the hydrolysis of glycosidic bonds and are key enzymes in carbohydrate metabolism. Efficient degradation of recalcitrant polysaccharides such as chitin and cellulose is accomplished due to synergistic enzyme cocktails consisting of accessory enzymes and mixtures of GHs with different modes of action and active site topologies. The substrate binding sites of chitinases and cellulases often have surface exposed aromatic amino acids and a tunnel or cleft topology. The active site of the exo-processive chitinase B (ChiB) from Serratia marcescens is partially closed, creating a tunnel-like catalytic cleft. To gain insight in the fundamental principles of substrate binding in this enzyme, we have studied the contribution of five key residues involved in substrate binding and tunnel formation to the thermodynamics of substrate binding. Mutation of Trp97, Phe190, Trp220 and Glu221, which are all part of the tunnel walls, resulted in significant less favorable conformational entropy change (ΔS°conf) upon binding (–TΔΔS°conf = ∼5 kcal/mol). This suggest that these residues are important for the structural rigidity and pre-shaping of the tunnel prior to binding. Mutation of Asp316, which, by forming a hydrogen bond to Trp97 is crucial in the active-site tunnel roof, resulted in a more favorable ΔS°conf relative to the wild type (–TΔΔS°conf = −2.2 kcal/mol). This shows that closing the tunnel-roof comes with an entropy cost, as previously suggested based on the crystal structures of GHs with tunnel topologies in complex with their substrates.

Entropy in Molecular Recognition by Proteins

At a fundamental level, biological processes are most often controlled using molecular recognition by proteins. Protein-ligand interactions impact critical events ranging from the catalytic action of enzymes, the assembly of macromolecular structures, complex signaling and allostery, transport phenomena, force generation and so on. The physical origin of high affinity interactions involving proteins continues to be the subject of intense investigation. Conformational entropy represents perhaps the last piece of the thermodynamic puzzle that governs protein structure, stability, dynamics and function. The presence and importance of internal conformational entropy in proteins has been debated for decades but has resisted experimental quantification. Over the past few years we have introduced, developed and validated an NMR-based approach that uses a dynamical proxy to determine changes in conformational entropy. This new approach, which we term the NMR “entropy meter,” requires few assumptions, is empirically calibrated and is apparently robust and universal. Using this “entropy meter,” it can now be quantitatively shown that proteins retain considerable conformational entropy in their native functional states and that this conformational entropy can play a decisive role in the thermodynamics of molecular recognition by proteins. Recent results show that changes conformational entropy of a protein upon binding a high affinity ligand is highly system specific and can vary from strongly inhibiting to even strongly promoting binding and everything in between. Thus one cannot possibly understand comprehensively how proteins work without knowledge of the breadth and underlying principles of the role of conformational entropy in protein function. This approach also allows for the refinement of empirical coefficients that relate changes in accessible surface area to changes in the entropy of water and the determination of the loss of rotational-translational entropy in high affinity protein complexes. Supported by the NIH and the Mathers Foundation.

Role of Conformational Entropy in Extremely High Affinity Protein Interactions

Interactions of extreme affinity (Kd ∼ fM) underlie many biochemical processes necessary to life. The physical determinants of such large binding energies are not well understood. Specific interactions at the interface (ΔHbinding) and the release of solvating water (TΔSsolvation) are usually assumed to dominate the binding energetics. The role of conformational entropy (TΔSconf) in determining binding affinity has remained elusive, in part due to the difficulties in measuring such changes in entropy experimentally. Recent developments in the Wand laboratory have bridged this gap by using solution NMR measurements of dynamics to empirically calibrate a “conformational entropy meter.” It has enabled quantitative measurements of the change in conformational entropy in protein-ligand binding. The toxin-antitoxin system studied here, barnase-barstar, forms a complex with fM affinity (ΔGbinding ∼ 19 kcal/mol) without undergoing any major structural changes and retaining a hydrated interface. To explore the role of conformational entropy, the fast (ps-ns time-scale) motions of backbone and side chains of the two proteins were measured in both the free (unbound) and the complexed (bound) states using NMR spectroscopy. Furthermore, hydration dynamics were measured in water and in the confined space of a reverse micelle. The dynamic response observed leads to an unfavorable change in TΔSconf, with a more rigid, still hydrated interface. This comprehensive study of both protein and “ligand” (in this case, another protein) and the measured site-specific changes in dynamics and hydration sheds light on the thermodynamic contributions that enable fM affinities. Supported by grants from the NIH, The Mathers Foundation and NSF.