Comparison of binding energies of SrcSH2-phosphotyrosyl peptides with structure-based prediction using surface area based empirical parameterization

Heat-dependent opening of TRPV1 in the presence of capsaicin

Transient receptor potential vanilloid member 1 (TRPV1) is a Ca²⁺-permeable cation channel that serves as the primary heat and capsaicin sensor in humans. Using cryo-EM, we have determined the structures of apo and capsaicin-bound full-length rat TRPV1 reconstituted into lipid nanodiscs over a range of temperatures. This has allowed us to visualize the noxious heat-induced opening of TRPV1 in the presence of capsaicin. Notably, noxious heat-dependent TRPV1 opening comprises stepwise conformational transitions. Global conformational changes across multiple subdomains of TRPV1 are followed by the rearrangement of the outer pore, leading to gate opening. Solvent-accessible surface area analyses and functional studies suggest that a subset of residues form an interaction network that is directly involved in heat sensing. Our study provides a glimpse of the molecular principles underlying noxious physical and chemical stimuli sensing by TRPV1, which can be extended to other thermal sensing ion channels.

Correlating structure and energetics in protein-ligand interactions: paradigms and paradoxes

Predicting protein-binding affinities of small molecules, even closely related ones, is a formidable challenge in biomolecular recognition and medicinal chemistry. A thermodynamic approach to optimizing affinity in protein-ligand interactions requires knowledge and understanding of how altering the structure of a small molecule will be manifested in protein-binding enthalpy and entropy changes; however, there is a relative paucity of such detailed information. In this review, we examine two strategies commonly used to increase ligand potency. The first of these involves introducing a cyclic constraint to preorganize a small molecule in its biologically active conformation, and the second entails adding nonpolar groups to a molecule to increase the amount of hydrophobic surface that is buried upon binding. Both of these approaches are motivated by paradigms suggesting that protein-binding entropy changes should become more favorable, but paradoxes can emerge that defy conventional wisdom.

Semi-supervised prediction of SH2-peptide interactions from imbalanced high-throughput data

Src homology 2 (SH2) domains are the largest family of the peptide-recognition modules (PRMs) that bind to phosphotyrosine containing peptides. Knowledge about binding partners of SH2-domains is key for a deeper understanding of different cellular processes. Given the high binding specificity of SH2, in-silico ligand peptide prediction is of great interest. Currently however, only a few approaches have been published for the prediction of SH2-peptide interactions. Their main shortcomings range from limited coverage, to restrictive modeling assumptions (they are mainly based on position specific scoring matrices and do not take into consideration complex amino acids inter-dependencies) and high computational complexity. We propose a simple yet effective machine learning approach for a large set of known human SH2 domains. We used comprehensive data from micro-array and peptide-array experiments on 51 human SH2 domains. In order to deal with the high data imbalance problem and the high signal-to-noise ration, we casted the problem in a semi-supervised setting. We report competitive predictive performance w.r.t. state-of-the-art. Specifically we obtain 0.83 AUC ROC and 0.93 AUC PR in comparison to 0.71 AUC ROC and 0.87 AUC PR previously achieved by the position specific scoring matrices (PSSMs) based SMALI approach. Our work provides three main contributions. First, we showed that better models can be obtained when the information on the non-interacting peptides (negative examples) is also used. Second, we improve performance when considering high order correlations between the ligand positions employing regularization techniques to effectively avoid overfitting issues. Third, we developed an approach to tackle the data imbalance problem using a semi-supervised strategy. Finally, we performed a genome-wide prediction of human SH2-peptide binding, uncovering several findings of biological relevance. We make our models and genome-wide predictions, for all the 51 SH2-domains, freely available to the scientific community under the following URLs: http://www.bioinf.uni-freiburg.de/Software/SH2PepInt/SH2PepInt.tar.gz and http://www.bioinf.uni-freiburg.de/Software/SH2PepInt/Genome-wide-predictions.tar.gz, respectively.

Using genome-wide measurements for computational prediction of SH2-peptide interactions

Peptide-recognition modules (PRMs) are used throughout biology to mediate protein–protein interactions, and many PRMs are members of large protein domain families. Recent genome-wide measurements describe networks of peptide–PRM interactions. In these networks, very similar PRMs recognize distinct sets of peptides, raising the question of how peptide-recognition specificity is achieved using similar protein domains. The analysis of individual protein complex structures often gives answers that are not easily applicable to other members of the same PRM family. Bioinformatics-based approaches, one the other hand, may be difficult to interpret physically. Here we integrate structural information with a large, quantitative data set of SH2 domain–peptide interactions to study the physical origin of domain–peptide specificity. We develop an energy model, inspired by protein folding, based on interactions between the amino-acid positions in the domain and peptide. We use this model to successfully predict which SH2 domains and peptides interact and uncover the positions in each that are important for specificity. The energy model is general enough that it can be applied to other members of the SH2 family or to new peptides, and the cross-validation results suggest that these energy calculations will be useful for predicting binding interactions. It can also be adapted to study other PRM families, predict optimal peptides for a given SH2 domain, or study other biological interactions, e.g. protein–DNA interactions.

Fast predictions of thermodynamics and kinetics of protein-protein recognition from structures: from molecular design to systems biology

The increasing call for an overall picture of the interactions between the components of a biological system that give rise to the observed function is often summarized by the expression systems biology. Both the interpretative and predictive capabilities of holistic models of biochemical systems, however, depend to a large extent on the level of physico-chemical knowledge of the individual molecular interactions making up the network. This review is focused on the structure-based quantitative characterization of protein-protein interactions, ubiquitous in any biochemical pathway. Recently developed, fast and effective computational methods are reviewed, which allow the assessment of kinetic and thermodynamic features of the association-dissociation processes of protein complexes, both in water soluble and membrane environments. The performance and the accuracy of fast and semi-empirical structure-based methods have reached comparable levels with respect to the classical and more elegant molecular simulations. Nevertheless, the broad accessibility and lower computational cost provide the former methods with the advantageous possibility to perform systems-level analyses including extensive in silico mutagenesis screenings and large-scale structural predictions of multiprotein complexes.

Thermodynamics of T-cell receptor-peptide/MHC interactions: progress and opportunities

alphabeta T-cell receptors (TCRs) recognize peptide antigens presented by class I or class II major histocompatibility complex molecules (pMHC). Here we review the use of thermodynamic measurements in the study of TCR-pMHC interactions, with attention to the diversity in binding thermodynamics and how this is related to the variation in TCR-pMHC interfaces. We show that there is no enthalpic or entropic signature for TCR binding; rather, enthalpy and entropy changes vary in a compensatory manner that reflects a narrow free energy window for the interactions that have been characterized. Binding enthalpy and entropy changes do not correlate with structural features such as buried surface area or the number of hydrogen bonds within TCR-pMHC interfaces, possibly reflecting the myriad of contributors to binding thermodynamics, but likely also reflecting a reliance on van't Hoff over calorimetric measurements and the unaccounted influence of equilibria linked to binding. TCR-pMHC binding heat capacity changes likewise vary considerably. In some cases, the heat capacity changes are consistent with conformational differences between bound and free receptors, but there is little data indicating these conformational differences represent the need to organize disordered CDR loops. In this regard, we discuss how thermodynamics may provide additional insight into conformational changes occurring upon TCR binding. Finally, we highlight opportunities for the further use of thermodynamic measurements in the study of TCR-pMHC interactions, not only for understanding TCR binding in general, but also for understanding specifics of individual interactions and the engineering of TCRs with desired molecular recognition properties.

Genome wide analysis of pathogenic SH2 domain mutations

The authors have made a genome-wide analysis of mutations in Src homology 2 (SH2) domains associated with human disease. Disease-causing mutations have been detected in the SH2 domains of cytoplasmic signaling proteins Bruton tyrosine kinase (BTK), SH2D1A, Ras GTPase activating protein (RasGAP), ZAP-70, SHP-2, STAT1, STAT5B, and the p85alpha subunit of the PIP3. Mutations in the BTK, SH2D1A, ZAP70, STAT1, and STAT5B genes have been shown to cause diverse immunodeficiencies, whereas the mutations in RASA1 and PIK3R1 genes lead to basal carcinoma and diabetes, respectively. PTPN11 mutations cause Noonan sydrome and different types of cancer, depending mainly on whether the mutation is inherited or sporadic. We collected and analyzed all known pathogenic mutations affecting human SH2 domains by bioinformatics methods. Among the investigated protein properties are sequence conservation and covariance, structural stability, side chain rotamers, packing effects, surface electrostatics, hydrogen bond formation, accessible surface area, salt bridges, and residue contacts. The majority of the mutations affect positions essential for phosphotyrosine ligand binding and specificity. The structural basis of the SH2 domain diseases was elucidated based on the bioinformatic analysis.

Genome-wide prediction of SH2 domain targets using structural information and the FoldX algorithm

Current experiments likely cover only a fraction of all protein-protein interactions. Here, we developed a method to predict SH2-mediated protein-protein interactions using the structure of SH2-phosphopeptide complexes and the FoldX algorithm. We show that our approach performs similarly to experimentally derived consensus sequences and substitution matrices at predicting known in vitro and in vivo targets of SH2 domains. We use our method to provide a set of high-confidence interactions for human SH2 domains with known structure filtered on secondary structure and phosphorylation state. We validated the predictions using literature-derived SH2 interactions and a probabilistic score obtained from a naive Bayes integration of information on coexpression, conservation of the interaction in other species, shared interaction partners, and functions. We show how our predictions lead to a new hypothesis for the role of SH2 domains in signaling.

Understanding the functional role of every protein in the cell is a long-standing goal of cellular biology. An important step in this direction is to discover how and when proteins interact inside the cell to accomplish their tasks. Many of the cellular functions depend on reversible protein modifications like phosphorylation. To sense these modifications, cells have protein domains capable of binding phosphorylated proteins such as the SH2 domain. In this work, we show that it is possible to use the three-dimensional structure of protein domains to predict its binding preferences. Using a computational tool called FoldX, we have predicted the binding specificity of several human SH2 domains. These predictions, based on the computational analysis of the 3-D structure, were shown to be of similar accuracy as those obtained from experimental binding assays. We show here that it is also possible to understand how a mutation changes the binding preference of protein binding domains, opening the way for better understanding of some disease causing mutations. The combination of this novel computational approach with other sources of information allowed us to provide a set of high-confidence novel interactions for the proteins here studied.

Prediction of Solvation Sites at the Interface of Src SH2 Domain Complexes Using Molecular Dynamics Simulations

Src Homology 2 (SH2) domains are approximately 100 amino acid domains that mediate recognition of tyrosine-phosphorylated sites by signalling proteins. Structures of SH2 domains with bound ligands indicate a potentially important role of water in influencing the binding thermodynamics. In this study, we used molecular dynamics (MD) simulation methods to evaluate solvation sites at the binding interface of the Src SH2 domain. We designed a software, WaRP (Water Residency Potential), to compute the positions of hydration sites from coordinates data of MD simulations and studied the impact of the computed positions on the prediction of the thermodynamics of Src SH2 domain binding to phosphorylated peptides using a method based on accessible surface area buried upon association. Two dually phosphorylated ligands and one monophosphorylated ligand were studied. We showed that the software predicted between 70% and 85% of the crystallographic water molecules depending on complexes. Comparison of the predicted water structures of both the bound and unbound binding partners led to a thorough evaluation of water behaviour during the binding reaction. We also showed that the predicted water structures of all ligand-SH2 domain structures investigated may be used to derive the entropy change provided that the heat capacity change is known. This study is the first to examine the dynamics of the water structure around the SH2 domain binding interface and contributes to our understanding of binding thermodynamics in SH2 domains.

Thermodynamic and Structural Basis for Transition-State Stabilization in Antibody-Catalyzed Hydrolysis

Catalytic antibodies 6D9 and 9C10, which were induced by immunization with a haptenic transition-state analog (TSA), catalyze the hydrolysis of a nonbioactive chloramphenicol monoester derivative to generate a bioactive chloramphenicol. These antibodies stabilize the transition state to catalyze the hydrolysis reaction, strictly according to the theoretical relationship: for 6D9, k(cat)/k(uncat)=895 and K(S)/K(TSA)=900, and for 9C10, k(cat)/k(uncat)=56 and K(S)/K(TSA)=60. To elucidate the molecular basis of the antibody-catalyzed reaction, the crystal structure of 6D9 was determined, and the binding thermodynamics of 6D9 and 9C10 with both the substrate and the TSA were analyzed using isothermal titration calorimetry. The crystal structure of the unliganded 6D9 Fab was determined at 2.25 A resolution and compared with that of the TSA-liganded 6D9 Fab reported previously, showing that the TSA is bound into the hydrophobic pocket of the antigen-combining site in an "induced fit" manner, especially at the L1 and H3 CDR loops. Thermodynamic analyses showed that 6D9 binds the substrate of the TSA with a positive DeltaS, differing from general thermodynamic characteristics of antigen-antibody interactions. This positive DeltaS could be due to the hydrophobic interactions between 6D9 and the substrate or the TSA mediated by Trp H100i. The difference in DeltaG between substrate and TSA-binding to 6D9 was larger than that to 9C10, which is in good correlation with the larger k(cat) value of 6D9. Interestingly, the DeltaDeltaG was mainly because of the DeltaDeltaH. The correlation between k(cat) and DeltaDeltaH is suggestive of "enthalpic strain" leading to destabilization of antibody-substrate complexes. Together with X-ray structural analyses, the thermodynamic analyses suggest that upon binding the substrate, the antibody alters the conformation of the ester moiety in the substrate from the planar Z form to a thermodynamically unstable twisted conformation, followed by conversion into the transition state. Enthalpic strain also contributes to the transition-state stabilization by destabilizing the ground state, and its degree is much larger for the more efficient catalytic antibody, 6D9.

The high-resolution NMR structure of the R21A Spc-SH3:P41 complex: understanding the determinants of binding affinity by comparison with Abl-SH3

Background

SH3 domains are small protein modules of 60–85 amino acids that bind to short proline-rich sequences with moderate-to-low affinity and specificity. Interactions with SH3 domains play a crucial role in regulation of many cellular processes (some are related to cancer and AIDS) and have thus been interesting targets in drug design. The decapeptide APSYSPPPPP (p41) binds with relatively high affinity to the SH3 domain of the Abl tyrosine kinase (Abl-SH3), while it has a 100 times lower affinity for the α-spectrin SH3 domain (Spc-SH3).

Results

Here we present the high-resolution structure of the complex between the R21A mutant of Spc-SH3 and p41 derived from NMR data. Thermodynamic parameters of binding of p41 to both WT and R21A Spc-SH3 were measured by a combination of isothermal titration and differential scanning calorimetry. Mutation of arginine 21 to alanine in Spc-SH3 increases 3- to 4-fold the binding affinity for p41 due to elimination at the binding-site interface of the steric clash produced by the longer arginine side chain. Amide hydrogen-deuterium experiments on the free and p41-bound R21A Spc-SH3 domain indicate that binding elicits a strong reduction in the conformational flexibility of the domain. Despite the great differences in the thermodynamic magnitudes of binding, the structure of the R21A Spc-SH3:P41 complex is remarkably similar to that of the Abl-SH3:P41 complex, with only few differences in protein-ligand contacts at the specificity pocket. Using empirical methods for the prediction of binding energetics based on solvent-accessible surface area calculations, the differences in experimental energetics of binding between the two complexes could not be properly explained only on the basis of the structural differences observed between the complexes. We suggest that the experimental differences in binding energetics can be at least partially ascribed to the absence in the R21A Spc-SH3:P41 complex of several buried water molecules, which have been proposed previously to contribute largely to the highly negative enthalpy and entropy of binding in the Abl-SH3:P41 complex.

Conclusion

Based on a deep structural and thermodynamic analysis of a low and high affinity complex of two different SH3 domains with the same ligand p41, we underline the importance of taking into account in any effective strategy of rational design of ligands, factors different from the direct protein-ligand interactions, such as the mediation of interactions by water molecules or the existence of cooperative conformational effects induced by binding.

The role of water in computational and experimental derivation of binding thermodynamics in SH2 domains

We have studied the role of bound interface water molecules on the prediction of the thermodynamics of SH2 domain binding to tyrosyl phosphopeptides using a method based on accessible surface area buried upon association. We studied three phosphopeptide ligands, which have been shown by Lubman and Waksman (J Mol Biol;328:655, 2003) and Davidson et al. (JACS;124:205, 2002) to have similar binding free energies but very different thermodynamic signatures. The thermodynamic model is semiempirical and applies to the crystal structure of the SH2 domain-bound forms. We explored all possible combinations of bound interfacial waters. We show that the model does not predict the binding thermodynamics of either ligand. However, we identified the empirical formula describing the heat capacity change as the source of the problem. Indeed, systematic exploration of heat capacity change values between 0 and -300 cal/mol deg results in a sharp distribution of the number of ligand/SH2/water-subset structures that provide binding thermodynamics similar to experimental values. The heat capacity change values at which the distributions peak are different for each peptide. This prompted us to experimentally determine the heat capacity change for each of the peptides and we found them to coincide with the values of the peaks. The implications of such findings are discussed.

Grid-based steered thermodynamic integration accelerates the calculation of binding free energies

The calculation of binding free energies is important in many condensed matter problems. Although formally exact computational methods have the potential to complement, add to, and even compete with experimental approaches, they are difficult to use and extremely time consuming. We describe a Grid-based approach for the calculation of relative binding free energies, which we call Steered Thermodynamic Integration calculations using Molecular Dynamics (STIMD), and its application to Src homology 2 (SH2) protein cell signalling domains. We show that the time taken to compute free energy differences using thermodynamic integration can be significantly reduced: potentially from weeks or months to days of wall-clock time. To be able to perform such accelerated calculations requires the ability to both run concurrently and control in realtime several parallel simulations on a computational Grid. We describe how the RealityGrid computational steering system, in conjunction with a scalable classical MD code, can be used to dramatically reduce the time to achieve a result. This is necessary to improve the adoption of this technique and further allows more detailed investigations into the accuracy and precision of thermodynamic integration. Initial results for the Src SH2 system are presented and compared to a reported experimental value. Finally, we discuss the significance of our approach.

Surface Plasmon Resonance Thermodynamic and Kinetic Analysis as a Strategic Tool in Drug Design. Distinct Ways for Phosphopeptides to Plug into Src- and Grb2 SH2 Domains

Thermodynamic and kinetic studies of biomolecular interactions give insight into specificity of molecular recognition processes and advance rational drug design. Binding of phosphotyrosine (pY)-containing peptides to Src- and Grb2-SH2 domains was investigated using a surface plasmon resonance (SPR)-based method. This SPR assay yielded thermodynamic binding constants in solution, and the kinetic information contained in the SPR signal allowed kinetic analysis, which demonstrated distinct ways for pY ligands to interact with the SH2 domains. The results for binding to Src SH2 were consistent with sequestration of water molecules in the interface of the pYEEI peptide/Src SH2 complex. The results for a pYVNV peptide binding to Grb2 SH2 suggested a conformational change for Grb2 SH2 upon binding, which is not observed for Src SH2. Binding of a cyclic construct, allowing the pYVNV sequence in the bound conformation, did not have the expected entropy advantage. The results suggest an alternative binding mode for this construct, with the hydrophobic ring-closing part interacting with the protein. In all cases, except for full-length Grb2 protein, the affinity for the immobilized peptide at the SPR sensor and in solution was identical. This study demonstrates that SPR thermodynamic and kinetic analysis is a useful strategic tool in drug design.

Two different T cell receptors use different thermodynamic strategies to recognize the same peptide/MHC ligand

A6 and B7 are two alphabeta T cell receptors (TCRs) that recognize the Tax peptide presented by the class I major histocompatibility molecule HLA-A2 (Tax/HLA-A2). Despite the fact that the two TCRs have different CDR loops and use different amino acid residues to contact their ligand, both receptors bind ligand with similar diagonal orientations. Here we show that they also bind with very similar binding affinities and kinetics (the DeltaDeltaG degrees for binding is approximately 0.3kcal/mol at 25 degrees C). The two receptors respond similarly to alterations in the MHC molecule, yet differ dramatically in their responses to ionic strength and temperature. The different responses to temperature indicate markedly different binding thermodynamics, which are not predictable from the surface area buried in the interfaces. A6 and B7 thus represent two TCRs that are both compatible with Tax/HLA-A2, although compatibility has been achieved through the use of different thermodynamic strategies. Finally, neither A6 nor B7 are predicted to undergo large conformational adaptations upon binding, distinguishing them from a number of other TCRs whose structure, thermodynamics, and kinetics have been characterized.

Calculation of affinities of peptides for proteins

Several methodologies were employed to calculate the Gibbs standard free energy of binding for a collection of protein-ligand complexes, where the ligand is a peptide and the protein is representative for various protein families. Almost 40 protein-ligand complexes were employed for a continuum approach, which considers the protein and the peptide at the atomic level, but includes solvent as a polarizable continuum. Five protein-ligand complexes were employed for an all-atom approach that relies on a combination of the double decoupling method with thermodynamic integration and molecular dynamics. These affinities were also computed by means of the linear interaction energy method. Although it generally proved rather difficult to predict the absolute free energies correctly, for some protein families the experimental ranking order was correctly reproduced by the continuum and all-atom approach. Considerable attention has also been given to correctly analyze the affinities of charged peptides, where it is required to judge the effect of one or more ions that are being decoupled in an all-atom approach to preserve electroneutrality. The various methods are further judged upon their merits.

Thermodynamic analysis of degenerate recognition by the NKG2D immunoreceptor: not induced fit but rigid adaptation

The homodimeric immunoreceptor NKG2D drives the activation of effector cells following engagement of diverse, conditionally expressed MHC class I-like protein ligands. NKG2D recognition is highly degenerate in that a single surface on receptor monomers binds pairs of distinct surfaces on each structurally divergent ligand, simultaneously accommodating multiple nonconservative ligand allelic or isoform substitutions. In contrast to TCR-pMHC and other NK receptor-ligand interactions, thermodynamic and kinetic analyses of four NKG2D-ligand pairs (MIC-A*001, MIC-B*005, ULBP1, and RAE-1beta) reported here show that the relative enthalpic and entropic terms, heat capacity, association rates, and activation energy barriers are comparable to typical, rigid protein-protein interactions. Rather than "induced-fit" binding, NKG2D degeneracy is achieved using distinct interaction mechanisms at each rigid interface.

A survey of the year 2003 literature on applications of isothermal titration calorimetry

Over the last decade isothermal titration calorimetry (ITC) has developed from a specialist method which was largely restricted in its use to dedicated experts, to a major, commercially available tool in the arsenal directed at understanding molecular interactions. The number of those proficient in this field has multiplied dramatically, as has the range of experiments to which this method has been applied. This has led to an overwhelming amount of new data and novel applications to be assessed. With the increasing number of publications in this field comes a need to highlight works of interest and impact. In this overview of the literature we have attempted to draw attention to papers and issues for which both the experienced calorimetrist and the interested dilettante hopefully will share our enthusiasm.

Specificity is complex and time consuming: mutual exclusivity in tyrosine kinase-mediated signaling

Most fundamental cellular processes are transduced through tyrosine kinase (TK)-mediated pathways. For transduction without corruption, the protein-protein interactions involved have to be mutually exclusive. Many of these proteins bind via homologous domains whose binding characteristics suggest that their innate specificity is not sufficiently high to account for the integrity of signal transduction. Stimulation of TK-mediated signals is often accompanied by recruitment of a precise, multimolecular protein complex that is itself capable of imposing specificity. Furthermore, this complex provides protection against phosphatase activity, controlling the longevity of the active signaling complex, and thus influencing outcomes in subsequent downstream events.

Structural parameterization of the binding enthalpy of small ligands

A major goal in ligand and drug design is the optimization of the binding affinity of selected lead molecules. However, the binding affinity is defined by the free energy of binding, which, in turn, is determined by the enthalpy and entropy changes. Because the binding enthalpy is the term that predominantly reflects the strength of the interactions of the ligand with its target relative to those with the solvent, it is desirable to develop ways of predicting enthalpy changes from structural considerations. The application of structure/enthalpy correlations derived from protein stability data has yielded inconsistent results when applied to small ligands of pharmaceutical interest (MW < 800). Here we present a first attempt at an empirical parameterization of the binding enthalpy for small ligands in terms of structural information. We find that at least three terms need to be considered: (1) the intrinsic enthalpy change that reflects the nature of the interactions between ligand, target, and solvent; (2) the enthalpy associated with any possible conformational change in the protein or ligand upon binding; and, (3) the enthalpy associated with protonation/deprotonation events, if present. As in the case of protein stability, the intrinsic binding enthalpy scales with changes in solvent accessible surface areas. However, an accurate estimation of the intrinsic binding enthalpy requires explicit consideration of long-lived water molecules at the binding interface. The best statistical structure/enthalpy correlation is obtained when buried water molecules within 5-7 A of the ligand are included in the calculations. For all seven protein systems considered (HIV-1 protease, dihydrodipicolinate reductase, Rnase T1, streptavidin, pp60c-Src SH2 domain, Hsp90 molecular chaperone, and bovine beta-trypsin) the binding enthalpy of 25 small molecular weight peptide and nonpeptide ligands can be accounted for with a standard error of +/-0.3 kcal x mol(-1).

Calorimetric and structural studies of 1,2,3-trisubstituted cyclopropanes as conformationally constrained peptide inhibitors of Src SH2 domain binding

Isothermal titration calorimetry and X-ray crystallography have been used to determine the structural and thermodynamic consequences associated with constraining the pTyr residue of the pYEEI ligand for the Src Homology 2 domain of the Src kinase (Src SH2 domain). The conformationally constrained peptide mimics that were used are cyclopropane-derived isosteres whereby a cyclopropane ring substitutes to the N-Calpha-Cbeta atoms of the phosphotyrosine. Comparison of the thermodynamic data for the binding of the conformationally constrained peptide mimics relative to their equivalent flexible analogues as well as a native tetrapeptide revealed an entropic advantage of 5-9 cal mol(-1) K(-1) for the binding of the conformationally constrained ligands. However, an unexpected drop in enthalpy for the binding of the conformationally constrained ligands relative to their flexible analogues was also observed. To evaluate whether these differences reflected conformational variations in peptide binding modes, we have determined the crystal structure of a complex of the Src SH2 domain bound to one of the conformationally constrained peptide mimics. Comparison of this new structure with that of the Src SH2 domain bound to a natural 11-mer peptide (Waksman et al. Cell 1993, 72, 779-790) revealed only very small differences. Hence, cyclopropane-derived peptides are excellent mimics of the bound state of their flexible analogues. However, a rigorous analysis of the structures and of the surface areas at the binding interface, and subsequent computational derivation of the energetic binding parameters, failed to predict the observed differences between the binding thermodynamics of the rigidified and flexible ligands, suggesting that the drop in enthalpy observed with the conformationally constrained peptide mimic arises from sources other than changes in buried surface areas, though the exact origin of the differences remains unclear.

Hierarchy of simulation models in predicting structure and energetics of the Src SH2 domain binding to tyrosyl phosphopeptides

Structure and energetics of the Src Src Homology 2 (SH2) domain binding with the recognition phosphopeptide pYEEI and its mutants are studied by a hierarchical computational approach. The proposed structure prediction strategy includes equilibrium sampling of the peptide conformational space by simulated tempering dynamics with the simplified, knowledge-based energy function, followed by structural clustering of the resulting conformations and binding free energy evaluation of a single representative from each cluster, a cluster center. This protocol is robust in rapid screening of low-energy conformations and recovers the crystal structure of the pYEEI peptide. Thermodynamics of the peptide-SH2 domain binding is analyzed by computing the average energy contributions over conformations from the clusters, structurally similar to the predicted peptide bound structure. Using this approach, the binding thermodynamics for a panel of studied peptides is predicted in a better agreement with the experiment than previously suggested models. However, the overall correlation between computed and experimental binding affinity remains rather modest. The results of this study show that small differences in binding free energies between the Ala and Gly mutants of the pYEEI peptide are considerably more difficult to predict than the structure of the bound peptides, indicating that accurate computational prediction of binding affinities still remains a major methodological and technical challenge.

Heat does not come in different colours: entropy-enthalpy compensation, free energy windows, quantum confinement, pressure perturbation calorimetry, solvation and the multiple causes of heat capacity effects in biomolecular interactions

Modern techniques in microcalorimetry allow us to measure directly the heat changes and associated thermodynamics for biomolecular processes in aqueous solution at reasonable concentrations. All these processes involve changes in solvation/hydration, and it is natural to assume that the heats for these processes should reflect, in some way, such changes in solvation. However, the interpretation of data is still somewhat ambiguous, since different non-covalent interactions may have similar thermodynamic signatures, and analysis is frustrated by large entropy-enthalpy compensation effects. Changes in heat capacity (Delta C(p)) have been related to changes in hydrophobic hydration and non-polar accessible surface areas, but more recent empirical and theoretical work has shown how this need not always be the case. Entropy-enthalpy compensation is a natural consequence of finite Delta C(p) values and, more generally, can arise as a result of quantum confinement effects, multiple weak interactions, and limited free energy windows, giving rise to thermodynamic homeostasis that may be of evolutionary and functional advantage. The new technique of pressure perturbation calorimetry (PPC) has enormous potential here as a means of probing solvation-related volumetric changes in biomolecules at modest pressures, as illustrated with preliminary data for a simple protein-inhibitor complex.

Heat capacity changes upon burial of polar and nonpolar groups in proteins

In this paper we address the question of whether the burial of polar and nonpolar groups in the protein locale is indeed accompanied by the heat capacity changes, DeltaC(p), that have an opposite sign, negative for nonpolar groups and positive for polar groups. To accomplish this, we introduced amino acid substitutions at four fully buried positions of the ubiquitin molecule (Val5, Val17, Leu67, and Gln41). We substituted Val at positions 5 and 17 and Leu at position 67 with a polar residue, Asn. As a control, Ala was introduced at the same three positions. We also replaced the buried polar Gln41 with Val and Leu, nonpolar residues that have similar size and shape as Gln. As a control, Asn was introduced at Gln41 as well. The effects of these amino acid substitutions on the stability, and in particular, on the heat capacity change upon unfolding were measured using differential scanning calorimetry. The effect of the amino acid substitutions on the structure was also evaluated by comparing the (1)H-(15)N HSQC spectra of the ubiquitin variants. It was found that the Ala substitutions did not have a considerable effect on the heat capacity change upon unfolding. However, the substitutions of aliphatic side chains (Val or Leu) with a polar residue (Asn) lead to a significant (> 30%) decrease in the heat capacity change upon unfolding. The decrease in heat capacity changes does not appear to be the result of significant structural perturbations as seen from the HSQC spectra of the variants. The substitution of a buried polar residue (Gln41) to a nonpolar residue (Leu or Val) leads to a significant (> 25%) increase in heat capacity change upon unfolding. These results indicate that indeed the heat capacity change of burial of polar and nonpolar groups has an opposite sign. However, the observed changes in DeltaC(p) are several times larger than those predicted, based on the changes in water accessible surface area upon substitution.

Inhibitors to the Src SH2 domain: a lesson in structure--thermodynamic correlation in drug design

Src homology 2 (SH2) domains play a key role in many tyrosine kinase-mediated intracellular signal transduction pathways. Aberrancies in the interaction of these domains can lead to a range of disease states. As a result, the pharmaceutical industry has made a large temporal and financial investment in the development of specific inhibitors to these domains. Focusing on the interactions of the SH2 domain from the protein Src, we report how the correlation of structural and thermodynamic data allows an assessment of the process of drug design. The binding site of the protein includes two pockets; one interacts with phosphotyrosine groups on cognate ligands, and the other accommodates an aliphatic hydrophobic side chain. The interaction with cognate ligands is also mediated by a network of water molecules. Thermodynamic data from isothermal titration calorimetric studies suggest that modification of the interactions in the SH2 binding site has been largely unsuccessful in producing high-affinity inhibitors. Furthermore, it appears that compounds that disrupt the interfacial water pay the price for the loss of the contribution to the free energy from a network of hydrogen bonds.