A Thermodynamic Analysis of the Binding Specificity between Four Human PDZ Domains and Eight Host, Viral and Designed Ligands.

PDZ domains are binding modules mostly involved in cell signaling and cell-cell junctions. These domains are able to recognize a wide variety of natural targets and, among the PDZ partners, viruses have been discovered to interact with their host via a PDZ domain. With such an array of relevant and diverse interactions, PDZ binding specificity has been thoroughly studied and a traditional classification has grouped PDZ domains in three major specificity classes. In this work, we have selected four human PDZ domains covering the three canonical specificity-class binding mode and a set of their corresponding binders, including host/natural, viral and designed PDZ motifs. Through calorimetric techniques, we have covered the entire cross interactions between the selected PDZ domains and partners. The results indicate a rather basic specificity in each PDZ domain, with two of the domains that bind their cognate and some non-cognate ligands and the two other domains that basically bind their cognate partners. On the other hand, the host partners mostly bind their corresponding PDZ domain and, interestingly, the viral ligands are able to bind most of the studied PDZ domains, even those not previously described. Some viruses may have evolved to use of the ability of the PDZ fold to bind multiple targets, with resulting affinities for the virus-host interactions that are, in some cases, higher than for host-host interactions.

Binding site plasticity in viral PPxY Late domain recognition by the third WW domain of human NEDD4.

The recognition of PPxY viral Late domains by the third WW domain of the HECT-E3 ubiquitin ligase NEDD4 (hNEDD4-WW3) is essential for the completion of the budding process of numerous enveloped viruses, including Ebola, Marburg, HTLV1 or Rabies. hNEDD4-WW3 has been validated as a promising target for the development of novel host-oriented broad spectrum antivirals. Nonetheless, finding inhibitors with good properties as therapeutic agents remains a challenge since the key determinants of binding affinity and specificity are still poorly understood. We present here a detailed structural and thermodynamic study of the interactions of hNEDD4-WW3 with viral Late domains combining isothermal titration calorimetry, NMR structural determination and molecular dynamics simulations. Structural and energetic differences in Late domain recognition reveal a highly plastic hNEDD4-WW3 binding site that can accommodate PPxY-containing ligands with varying orientations. These orientations are mostly determined by specific conformations adopted by residues I859 and T866. Our results suggest a conformational selection mechanism, extensive to other WW domains, and highlight the functional relevance of hNEDD4-WW3 domain conformational flexibility at the binding interface, which emerges as a key element to consider in the search for potent and selective inhibitors of therapeutic interest.

Protein Folding Cooperativity and Thermodynamic Barriers of the Simplest ß-Sheet Fold: A Survey of WW Domains.

Theory and experiments have shown that microsecond folding proteins exhibit characteristic thermodynamic properties that reflect the limited cooperativity of folding over marginal barriers (downhill folding). Those studies have mostly focused on proteins with large α-helical contents and small size, which tend to be the fastest folders. A key open question is whether such properties are also present in the fastest all-ß proteins. We address this issue by investigating the unfolding thermodynamics of a collection of WW domains as representatives of the simplest ß-sheet fold. WW domains are small microsecond folders, although they do not fold as fast as their α-helical counterparts. In previous work on the NEDD4-WW4 domain, we reported deviations from two-state thermodynamics that were less apparent and thus suggestive of an incipient downhill scenario. Here we investigate the unfolding thermodynamics of four other WW domains (NEDD4-WW3, YAP65-WW1(L30K), FBP11-WW1, and FBP11-WW2) by performing all of the thermodynamic tests for downhill folding that have been previously developed on α-helical proteins. This set of five WW domains shares low sequence identity and include examples from two specificity classes, thus providing a comprehensive survey. Thermodynamic analysis of the four new WW domains consistently reveals all of the properties of downhill folding equilibria, which are in all cases more marked than what we found before in NEDD4-WW4. Our results show that fast-folding all-ß proteins do share limited cooperativity and gradual unfolding thermodynamics with fast α-helical proteins and suggest that the free energy barrier to folding of natural proteins is mostly determined by size and fold topology and much less by the specific amino acid sequence.

Approaching the thermodynamic view of protein folding through the reproduction of Anfinsen's experiment by undergraduate physical biochemistry students.

In 1972 Christian B. Anfinsen received the Nobel Prize in Chemistry for "…his work on ribonuclease, especially concerning the connection between the amino acid sequence and the biologically active conformation." The understanding of this principle is crucial for physical biochemistry students, since protein folding studies, bio-computing sciences and protein design approaches are founded on such a well-demonstrated connection. Herein, we describe a detailed and easy-to-follow experiment to reproduce the most relevant assays carried out at Anfinsen's laboratory in the 60s. This experiment provides students with a platform to interpret by themselves the structural and kinetic experiments conceived to understand the protein folding problem. In addition, this three-day experiment brings students a nice opportunity for protein manipulation as well as for the setting up of spectroscopic and chromatographic techniques. © 2018 by The International Union of Biochemistry and Molecular Biology, 46(3):262-269, 2018.

WW domains of the yes-kinase-associated-protein (YAP) transcriptional regulator behave as independent units with different binding preferences for PPxY motif-containing ligands.

YAP is a WW domain-containing effector of the Hippo tumor suppressor pathway, and the object of heightened interest as a potent oncogene and stemness factor. YAP has two major isoforms that differ in the number of WW domains they harbor. Elucidating the degree of co-operation between these WW domains is important for a full understanding of the molecular function of YAP. We present here a detailed biophysical study of the structural stability and binding properties of the two YAP WW domains aimed at investigating the relationship between both domains in terms of structural stability and partner recognition. We have carried out a calorimetric study of the structural stability of the two YAP WW domains, both isolated and in a tandem configuration, and their interaction with a set of functionally relevant ligands derived from PTCH1 and LATS kinases. We find that the two YAP WW domains behave as independent units with different binding preferences, suggesting that the presence of the second WW domain might contribute to modulate target recognition between the two YAP isoforms. Analysis of structural models and phage-display studies indicate that electrostatic interactions play a critical role in binding specificity. Together, these results are relevant to understand of YAP function and open the door to the design of highly specific ligands of interest to delineate the functional role of each WW domain in YAP signaling.

Single-chain protein mimetics of the N-terminal heptad-repeat region of gp41 with potential as anti-HIV-1 drugs.

During HIV-1 fusion to the host cell membrane, the N-terminal heptad repeat (NHR) and the C-terminal heptad repeat (CHR) of the envelope subunit gp41 become transiently exposed and accessible to fusion inhibitors or Abs. In this process, the NHR region adopts a trimeric coiled-coil conformation that can be a target for therapeutic intervention. Here, we present an approach to rationally design single-chain protein constructs that mimic the NHR coiled-coil surface. The proteins were built by connecting with short loops two parallel NHR helices and an antiparallel one with the inverse sequence followed by engineering of stabilizing interactions. The constructs were expressed in Escherichia coli, purified with high yield, and folded as highly stable helical coiled coils. The crystal structure of one of the constructs confirmed the predicted fold and its ability to accurately mimic an exposed gp41 NHR surface. These single-chain proteins bound to synthetic CHR peptides with very high affinity, and furthermore, they showed broad inhibitory activity of HIV-1 fusion on various pseudoviruses and primary isolates.

The impact of extra-domain structures and post-translational modifications in the folding/misfolding behaviour of the third PDZ domain of MAGUK neuronal protein PSD-95.

The modulation of binding affinities and specificities by post-translational modifications located out from the binding pocket of the third PDZ domain of PSD-95 (PDZ3) has been reported recently. It is achieved through an intra-domain electrostatic network involving some charged residues in the ß2-ß3 loop (were a succinimide modification occurs), the α3 helix (an extra-structural element that links the PDZ3 domain with the following SH3 domain in PSD-95, and contains the phosphorylation target Tyr397), and the ligand peptide. Here, we have investigated the main structural and thermodynamic aspects that these structural elements and their related post-translational modifications display in the folding/misfolding pathway of PDZ3 by means of site-directed mutagenesis combined with calorimetry and spectroscopy. We have found that, although all the assayed mutations generate proteins more prone to aggregation than the wild-type PDZ3, those directly affecting the α3 helix, like the E401R substitution or the truncation of the whole α3 helix, increase the population of the DSC-detected intermediate state and the misfolding kinetics, by organizing the supramacromolecular structures at the expense of the two ß-sheets present in the PDZ3 fold. However, those mutations affecting the ß2-ß3 loop, included into the prone-to-aggregation region composed by a single ß-sheet comprising ß2 to ß4 chains, stabilize the trimeric intermediate previously shown in the wild-type PDZ3 and slow-down aggregation, also making it partly reversible. These results strongly suggest that the α3 helix protects to some extent the PDZ3 domain core from misfolding. This might well constitute the first example where an extra-element, intended to link the PDZ3 domain to the following SH3 in PSD-95 and in other members of the MAGUK family, not only regulates the binding abilities of this domain but it also protects PDZ3 from misfolding and aggregation. The influence of the post-translational modifications in this regulatory mechanism is also discussed.

A thermodynamic study of the third PDZ domain of MAGUK neuronal protein PSD-95 reveals a complex three-state folding behavior.

The relevance of the C-terminal α helix of the PDZ3 domain of PSD95 in its unfolding process has been explored by achieving the thermodynamic characterization of a construct where the sequence of the nine residues corresponding to such motif has been deleted. Calorimetric traces at neutral pH require the application of a three-state model displaying three different equilibrium processes in which the intermediate state self-associates upon heating, being stable and populated in a wide temperature range. Temperature scans followed by circular dichroism, Fourier transform infrared spectroscopy and dynamic light scattering support the presence of such oligomeric-partially folded species. This study reveals that the deletion of the α3-helix sequence results in a more complex description of the domain unfolding.

Thermodynamic impact of embedded water molecules in the unfolding of human CD2BP2-GYF domain.

GYF domains are small polyproline-recognition modules adopting a structural arrangement consisting of a single α-helix packed against a small ß-sheet. Although most families of proline-rich recognition modules have been extensively characterized in terms of function, structure, or conformational flexibility, little is known about GYF domain functionality and folding. We have undertaken the thermodynamic characterization of the unfolding of CD2BP2-GYF domain by combining differential scanning calorimetry and circular dichroism under different pH conditions. The experimental data can be well-described in terms of a two-state equilibrium, although an unusually high heat capacity of the native state reflects a considerable conformational flexibility and dynamics of CD2BP2-GYF domain. In addition, the normalized thermodynamic parameters of unfolding (enthalpy, entropy and heat capacity) are roughly a factor of two greater than expected. In contrast, stability curves reveal an ordinary unfolding behavior of CD2BP2-GYF domain in terms of Gibbs energies, incurring thus unusually strong enthalpy-entropy compensation. This phenomenon, previously described as "thermodynamic homeostasis", has been associated in different examples to the contribution of occluded water (solvent) molecules into the protein structure. By means of CASTp server, we have found seven cavities/pockets scattered throughout of the CD2BP2-GYF structure, each able to harbor at least one water molecule. This structural feature provides rationalization for the atypical enthalpy values observed for CD2BP2-GYF because each water molecule is able to organize an extra amount of hydrogen bonds in the native state. In addition, these bound waters increase the vibrational entropy of the protein, which could also be responsible for an increase in protein flexibility and may thus fully explain the homeostatic behavior experimentally observed.

A comparative analysis of the folding and misfolding pathways of the third PDZ domain of PSD95 investigated under different pH conditions.

Equilibrium unfolding at neutral pH of the third PDZ domain of PSD95 is well described by the presence of a partly unfolded intermediate that presents association phenomena. After some days' incubation annular and fibrillar structures form from the oligomers. At pH values below 3, however, differential scanning calorimetry shows that PDZ3 seems to unfold under a two-state scheme. Kinetic measurements followed by dynamic light scattering, ThT and ANS fluorescence reveal that the misfolding pathway still exists despite the absence of any populated intermediates and shows an irreversible assembling of the supramacromolecular structures as well as an appreciable lag-phase, contrary to what is found in similar experiments at neutral pH. Moreover, as shown by transmission-electron-microscopy images, the annular structures seen at neutral pH completely disappear from incubated solutions. According to the structural information, this titration behavior appears to be the consequence of a conformational equilibrium that depends on the protonation of some Glu residues located at the C-terminal α3 helix and at the hairpin formed by strands ß2 and ß3. Our calculations suggest that the enthalpic contribution of these interactions may well be as much as 40kJ·mol(-1). The possible regulatory role of this equilibrium upon PDZ3 functionality and amyloid formation is briefly discussed.

A thermodynamic characterization of the interaction of 8-anilino-1-naphthalenesulfonic acid with native globular proteins: the effect of the ligand dimerization in the analysis of the binding isotherms.

8-Anilino-1-naphthalenesulfonic acid (ANS) is a popular fluorescence probe, broadly used for the analysis of proteins, but the nature of its interaction with proteins and the high increase in the fluorescence intensity that takes place upon such process are still unclear. In the last few years, isothermal titration calorimetry has been used to characterize the nature of the interaction of this dye with proteins. The analysis of the binding isotherms of these studies has not considered the dimerization equilibrium of ANS, which is pH dependent, and it can result in serious errors in the data analysis. In the present work we have developed a suitable data analysis by which this process is taken into account. To study the binding of the dye to proteins at different pH values, we have used the Abl-SH3 domain. Our results suggest that at pH 3 and 5, where the dimerization of the ANS is important, electrostatic interactions are significant for the binding of ANS to the Abl-SH3 domain. However, at pH 7, ANS behaves mostly as monomer and the interaction with the protein is mainly hydrophobic. The pH dependent behavior of the ANS binding to proteins can be explained in terms of ionization states of both, the protein and the ANS.

High-resolution structure of an alpha-spectrin SH3-domain mutant with a redesigned hydrophobic core.

The alpha-spectrin SH3 domain (Spc-SH3) is a small modular domain which has been broadly used as a model protein in folding studies and these studies have sometimes been supported by structural information obtained from the coordinates of Spc-SH3 mutants. The structure of B5/D48G, a multiple mutant designed to improve the hydrophobic core and as a consequence the protein stability, has been solved at 1 A resolution. The crystals belonged to the orthorhombic space group P2(1)2(1)2(1), with unit-cell parameters a=24.79, b=37.23, c=62.95 A. This mutant also bears a D48G substitution in the distal loop and this mutation has also been reported to increase the stability of the protein by itself. The structure of the B5/D48G mutant shows a highly packed hydrophobic core and a more ordered distal loop compared with previous Spc-SH3 structures.

An oligomeric equilibrium intermediate as the precursory nucleus of globular and fibrillar supramacromolecular assemblies in a PDZ domain.

The equilibrium unfolding at neutral pH of the third PDZ domain of PSD95, as followed by DSC, is characterized by the presence of an equilibrium intermediate with clear signs of oligomerization. DLS and SEC measurements indicate that at 60-70 degrees C small oligomers populate, showing a typical beta-sheet far-UV CD spectrum. These intermediate species lead to the formation of rodlike particulates of approximately 12 nm, which remain in solution after 2 weeks incubation and grow until they adopt annular/spherical shapes of approximately 50 nm and protofibrils, which are subsequently fully transformed into fibrils. The fibrils can also disaggregate after the addition of 1:1 buffer dilution followed by cooling to room temperature, thus returning to the initial monomeric state. Growth kinetics, as shown by ThT and ANS fluorescence, show that the organization of the different supramacromolecular structures comes from a common nucleation unit, the small oligomers, which organize themselves before reaching the incubation temperature of 60 degrees C. Our experiments point toward the existence of a well-defined reversible, stepwise, and downhill organization of the processes involved in the association-dissociation of the intermediate. We estimate the enthalpy change accompanying the association-dissociation equilibria to be 130 kJ x mol(-1). Furthermore, the coalescence under essentially reversible conditions of different kinds of supramacromolecular assemblies renders this protein system highly interesting for biophysical studies aimed at our further understanding of amyloid pathological conditions.

Novel conformational aspects of the third PDZ domain of the neuronal post-synaptic density-95 protein revealed from two 1.4A X-ray structures.

The crystal structure of the third PDZ domain of the neuronal post-synaptic density-95 protein (PSD95-PDZ3, residues 302-402) has been solved at 1.4 and 1.35A from two different crystal forms. These structures lack the cloning artefact present in the carboxyl terminal sequence of the former crystallographic structures and they belong to the space groups P4(3) and P1. The new PDZ structures are identical between the two crystal forms and among the four chains of the P1 crystal form. When we compare the new structures with the previous ones, some important conformational differences in the C-terminal alpha-helix and in the loop connecting beta2 and beta3 strands have been found. Additionally, the high resolution of the new structures has allowed us to indentify a succinimide residue at the position corresponding to Asp332 in the beta2-beta3 loop, which may contribute to the alternate conformation of this loop, and at the same time, to the interaction between residues from this loop and the C-terminal alpha-helix. Thus, these features would have implications in the recently proposed allosteric role of this third alpha-helix in the binding of the carboxyl terminal fragments to the PSD95-PDZ3.

Thermodynamic characterization of the folding equilibrium of the human Nedd4-WW4 domain: at the frontiers of cooperative folding.

WW domains are the smallest naturally independent beta-sheet protein structures available to date and constitute attractive model systems for investigating the determinants of beta-sheet folding and stability. Nonetheless, their small size and low cooperativity pose a difficult challenge for a quantitative analysis of the folding equilibrium. We describe here a comprehensive thermodynamic characterization of the conformational equilibrium of the fourth WW domain from the human ubiquitin ligase Nedd4 (hNedd4-WW4) using a combination of calorimetric and spectroscopic techniques with several denaturing agents (temperature, pH, and chemical denaturants). Our results reveal that even though the experimental data can be described in terms of a two-state equilibrium, spectral data together with anomalous values for some thermodynamic parameters (a strikingly low temperature of maximum stability, a higher than expected native-state heat capacity, and a small specific enthalpy of unfolding) could be indicative of more complex types of equilibria, such as one-state downhill folding or alternative native conformations. Moreover, double-perturbation experiments reveal some features that, in spite of the apparent linear correlation between the thermodynamic parameters, seem to be indicative of a complex conformational equilibrium in the presence of urea. In summary, the data presented here point toward the existence of a low-energy barrier between the different macrostates of hNedd4-WW4, placing it at the frontier of cooperative folding.

Evaluation of folding co-operativity of a chimeric protein based on the molecular recognition between polyproline ligands and SH3 domains.

In previous work, we designed a chimeric protein, named SPCp41, to evaluate the thermodynamics of the interaction between SH3 domains and proline-rich ligands by combining thermal unfolding measurements and mutagenesis. Here, we have investigated the energetic integrity of the chain extension corresponding to the ligand sequence into the native structure, since the opposite will produce changes in the folding mechanism of the SH3 domain that may give rise to undesirable contributions to the thermodynamic parameters. We have analysed the folding-unfolding kinetics under standard conditions (50 mM phosphate pH 7). Kinetic evolutions are well described by a bi-exponential where, on top of the main kinetic phase, a low-populated slower phase appears as a consequence of cis-trans isomerisation of Pro39, as demonstrated by the influence of prolyl isomerases and by mutational analysis. There is also a burst phase possibly due to a productive formation of some helical ensembles. The main evolution, accounting for the true folding kinetics of SPCp41, can be considered as a two-state process, where the folding transition state produces essentially the same picture shown by the circular permutant S19-P20s (the 'nucleus' of the design) and the ligand will dock at the latter stages of the two-state process. Thus, all conclusions argue in favour of the effectiveness of SPCp41 to study energetic, dynamic and structural aspects of SH3-ligand interactions.

An error analysis for two-state protein-folding kinetic parameters and phi-values: progress toward precision by exploring pH dependencies on Leffler plots.

The interpretation of phi-values has led to an understanding of the folding transition state ensemble of a variety of proteins. Although the main guidelines and equations for calculating phi are well established, there remains some controversy about the quality of the numerical values obtained. By analyzing a complete set of results from kinetic experiments with the SH3 domain of alphaspectrin (Spc-SH3) and applying classical error methods and error-propagation formulas, we evaluated the uncertainties involved in two-state-folding kinetic experimental parameters and the corresponding calculated phi-values. We show that kinetic constants in water and m values can be properly estimated from a judicious weighting of fitting errors and describe some procedures to calculate the errors in Gibbs energies and phi-values from a traditional two-point Leffler analysis. Furthermore, on the basis of general assumptions made with the protein engineering method, we show how to generate multipoint Leffler plots via the analysis of pH dependencies of kinetic parameters. We calculated the definitive phi-values for a collection of single mutations previously designed to characterize the folding transition state of the alphaspectrin SH3 domain. The effectiveness of the pH-scanning procedure is also discussed in the context of error analysis. Judging from the magnitudes of the error bars obtained from two-point and multipoint Leffler plots, we conclude that the precision obtained for phi-values should be approximately 25%, a reasonable limit that takes into account the propagation of experimental errors.

Sulfate-induced effects in the on-pathway intermediate of the bacterial immunity protein Im7.

Intermediates have now been identified in the folding of a number of small, single-domain proteins. Here we describe experiments to determine the effect of Na(2)SO(4) on the properties of the on-pathway intermediate formed early during the folding of the four-helical protein, Im7. This intermediate, studied previously in 0.4 M Na(2)SO(4), contains three of the four native helices and is fascinating in that several residues in helices I, II, and IV make non-native interactions that stabilize this state. Whether these contacts form as a consequence of the presence of Na(2)SO(4), however, remained unresolved. Using kinetic analysis of the effect of Na(2)SO(4) on the unfolding and refolding kinetics of Im7*, combined with detailed analysis of the resulting chevron plots, we show that decreasing the concentration of Na(2)SO(4) from 0.4 to 0 M destabilizes the intermediate and rate-limiting transition (TS2) states by 7 and 10 kJ mol(-)(1), respectively, and has little effect on the relative compactness of these states compared with that of the unfolded ensemble (beta(I) approximately 0.8, beta(TS2) approximately 0.9 in 0 to 0.4 M Na(2)SO(4)). Analysis of 10 variants of the protein in 0.2 M Na(2)SO(4) using Phi-values showed that the structural properties of the intermediate and TS2 are not altered significantly by the concentration of the kosmotrope. The data demonstrate that the rapid formation of a compact intermediate stabilized by non-native interactions during Im7* folding is not induced by high concentrations of the stabilizing salt, but is a generic feature of the folding of this protein.

A miniprotein scaffold used to assemble the polyproline II binding epitope recognized by SH3 domains.

SH3 domains are molecular-recognition modules that function by interacting with proteins containing sequences in polyproline II (PPII) conformation. The main limitation in designing short-ligand peptides to interact with these domains is the preservation of this helical arrangement, for which a high content of proline is needed. We have overcome this limitation by using a protein scaffold provided by the avian pancreatic polypeptide (APP), a natural hormone of 36 amino acid residues. The APP protein contains a PPII stretch packed against an alpha-helix. We have designed a structure in which some residues of the APP PPII helix are replaced by a sequence motif, named RP1, which interacts with the SH3 domain of the Abelson tyrosine kinase (Abl-SH3). This design, which we call APP-RP1, is folded and, as shown by circular dichroism, has a structural content similar to that of natural APP (APP-WT). The stability of both miniproteins has been compared by unfolding experiments; the designed APP-RP1 is almost 20 deg. C more stable than the wild-type and has a higher Gibbs energy function. This increase in stability has an entropic origin. Isothermal titration calorimetry and fluorescence spectroscopy show that the thermodynamics of the binding of the APP-RP1 molecule to Abl-SH3 is comparable to that of the shorter RP1 peptide. Furthermore, the mutation by Tyr of two proline residues in APP-RP1, which are essential for the binding of some linear peptides to Abl-SH3, demonstrates the effectiveness of the scaffold in enhancing the variability in the design of high-affinity and high-specificity ligands for any SH3 domain. The application of this strategy may help in the design of ligands for other polyproline-recognition domains such as WW, PX or EVH1, and even for the in vivo application of these miniproteins.

Thermodynamic dissection of the binding energetics of proline-rich peptides to the Abl-SH3 domain: implications for rational ligand design.

The inhibition of the interactions between SH3 domains and their targets is emerging as a promising therapeutic strategy. To date, rational design of potent ligands for these domains has been hindered by the lack of understanding of the origins of the binding energy. We present here a complete thermodynamic analysis of the binding energetics of the p41 proline-rich decapeptide (APSYSPPPPP) to the SH3 domain of the c-Abl oncogene. Isothermal titration calorimetry experiments have revealed a thermodynamic signature for this interaction (very favourable enthalpic contributions opposed by an unfavourable binding entropy) inconsistent with the highly hydrophobic nature of the p41 ligand and the Abl-SH3 binding site. Our structural and thermodynamic analyses have led us to the conclusion, having once ruled out any possible ionization events or conformational changes coupled to the association, that the establishment of a complex hydrogen-bond network mediated by water molecules buried at the binding interface is responsible for the observed thermodynamic behaviour. The origin of the binding energetics for proline-rich ligands to the Abl-SH3 domain is further investigated by a comparative calorimetric analysis of a set of p41-related ligands. The striking effects upon the enthalpic and entropic contributions provoked by conservative substitutions at solvent-exposed positions in the ligand confirm the complexity of the interaction. The implications of these results for rational ligand design are discussed.