Backbone Entropy Research Articles

Saddle-Shaped Third Component with Out-of-Plane Electrostatic Dipole for Realizing High-Performance Photovoltaic Donor Terpolymers.

Terpolymer fabrication is an effective methodology for molecular engineering and generating high-performance organic photovoltaic materials to construct highly efficient polymer solar cells. Modification of the polymer PM6 by incorporating a third component resulting in the formation of a ternary copolymer is reported to outperform PM6 in achieving enhanced device performances. However, one of the major challenges in constructing high-performance terpolymers is to counter the molecular disorder caused by the backbone entropy induced by the third moiety. In this work, double B←N bridged bipyridine (BNBP) is used as the third component, which possesses astrong out-of-plane electrostatic dipole owing to the saddle-shaped B←N fused ring structure. The out-of-plane dipole moment introduced in the modified PM6 terpolymer can be used as a means for tuning and optimizing the nanostructures of the blended films. The prepared PM6-BNBP-4 blend polymer with 4% of the benzodithiophene dione monomers replaced by BNBP results in excellent power conversion efficiency of 19.13%. This work demonstrates that the out-of-plane electrostatic dipole moment in saddle-shaped molecules is valuable for achieving high-performance organic photovoltaic donor materials.

Thermodynamics of Conformational Transitions in a Disordered Protein Backbone Model

Conformational entropy is expected to contribute significantly to the thermodynamics of structural transitions in intrinsically disordered proteins or regions in response to protein/ligand binding, posttranslational modifications, and environmental changes. We calculated the backbone (dihedral) conformational entropy of oligoglycine (GlyN), a protein backbone mimic and model intrinsically disordered region, as a function of chain length (N=3, 4, 5, 10, and 15) from simulations using three different approaches. The backbone conformational entropy scales linearly with chain length with a slope consistent with the entropy of folding of well-structured proteins. The entropic contributions of second-order dihedral correlations are predominantly through intraresidue ϕ-ψ pairs, suggesting that oligoglycine may be thermodynamically modeled as a system of independent glycine residues. We find the backbone conformational entropy to be largely independent of global structural parameters, like the end-to-end distance and radius of gyration. We introduce a framework referred to herein as “ensemble confinement” to estimate the loss (gain) of conformational free energy and its entropic component when individual residues are constrained to (released from) particular regions of the ϕ-ψ map. Quantitatively, we show that our protein backbone model resists ordering/folding with a significant, unfavorable ensemble confinement free energy because of the loss of a substantial portion of the absolute backbone entropy. Proteins can couple this free-energy reservoir to distal binding events as a regulatory mechanism to promote or suppress binding.

Peptide Sequence and Solvent as Levers to Control Disulfide Connectivity in Multiple Cysteine Containing Venom Toxins.

Judicious choice of solvent, temperature, and strategic mutations along a peptide backbone can minimize formation of non-native disulfide bond isoforms in chemical synthesis of multiple cysteine containing venom toxins. By exploiting these controls, one can drive the population distribution in favor of a particular isoform. Some chosen ionic liquids (ILs), like 1-ethyl-3-methyl-imidazolium acetate, [Im21][OAc], have proven efficient in favoring the native globular isoform in some conotoxins. To comprehend such a preference, we report an explicit solvent replica exchange molecular dynamics (REMD) study of two conotoxins, AuIB and GI, solvated in either neat water or ∼50% (v/v) mixture of water-[Im21][OAc]. Our simulations indicate that compared to neat water, the probability of obtaining native globular isoform of AuIB significantly increases in a water-IL mixture at 305 K. Strikingly, and aligned with experimental observations, peptide GI does not favor the native connectivity in the water-IL mixture. In presence of IL, strong solvent mediated fluctuations of the GI backbone are observed in our simulations. Uneven ion accumulation along the backbone owing to strong H-bonding interactions of some GI residues with IL ions, especially the anion OAc-, restricts conformational freedom of the peptide. Estimation of backbone entropy and Helmholtz free energy corroborates the lack of conformational freedom in GI as compared to AuIB, especially in the presence of IL. In line with prior experiments, simulations of GI mutants indicate that one could possibly force a given pair of Cys residues to come closer by strategically mutating GI residues with glycine and/or alanine, resulting in the breakage/formation of helix-like motifs.

Local and global anatomy of antibody-protein antigen recognition.

Deciphering antibody-protein antigen recognition is of fundamental and practical significance. We constructed an antibody structural dataset, partitioned it into human and murine subgroups, and compared it with nonantibody protein-protein complexes. We investigated the physicochemical properties of regions on and away from the antibody-antigen interfaces, including net charge, overall antibody charge distributions, and their potential role in antigen interaction. We observed that amino acid preference in antibody-protein antigen recognition is entropy driven, with residues having low side-chain entropy appearing to compensate for the high backbone entropy in interaction with protein antigens. Antibodies prefer charged and polar antigen residues and bridging water molecules. They also prefer positive net charge, presumably to promote interaction with negatively charged protein antigens, which are common in proteomes. Antibody-antigen interfaces have large percentages of Tyr, Ser, and Asp, but little Lys. Electrostatic and hydrophobic interactions in the Ag binding sites might be coupled with Fab domains through organized charge and residue distributions away from the binding interfaces. Here we describe some features of antibody-antigen interfaces and of Fab domains as compared with nonantibody protein-protein interactions. The distributions of interface residues in human and murine antibodies do not differ significantly. Overall, our results provide not only a local but also a global anatomy of antibody structures.

Single-stranded nucleic acid elasticity arises from internal electrostatic tension

Understanding of the conformational ensemble of flexible polyelectrolytes, such as single-stranded nucleic acids (ssNAs), is complicated by the interplay of chain backbone entropy and salt-dependent electrostatic repulsions. Molecular elasticity measurements are sensitive probes of the statistical conformation of polymers and have elucidated ssNA conformation at low force, where electrostatic repulsion leads to a strong excluded volume effect, and at high force, where details of the backbone structure become important. Here, we report measurements of ssDNA and ssRNA elasticity in the intermediate-force regime, corresponding to 5- to 100-pN forces and 50-85% extension. These data are explained by a modified wormlike chain model incorporating an internal electrostatic tension. Fits to the elastic data show that the internal tension decreases with salt, from [Formula: see text]5 pN under 5 mM ionic strength to near zero at 1 M. This decrease is quantitatively described by an analytical model of electrostatic screening that ascribes to the polymer an effective charge density that is independent of force and salt. Our results thus connect microscopic chain physics to elasticity and structure at intermediate scales and provide a framework for understanding flexible polyelectrolyte elasticity across a broad range of relative extensions.

The Dynameomics Entropy Dictionary: A Large-Scale Assessment of Conformational Entropy across Protein Fold Space.

Molecular dynamics (MD) simulations contain considerable information with regard to the motions and fluctuations of a protein, the magnitude of which can be used to estimate conformational entropy. Here we survey conformational entropy across protein fold space using the Dynameomics database, which represents the largest existing data set of protein MD simulations for representatives of essentially all known protein folds. We provide an overview of MD-derived entropies accounting for all possible degrees of dihedral freedom on an unprecedented scale. Although different side chains might be expected to impose varying restrictions on the conformational space that the backbone can sample, we found that the backbone entropy and side chain size are not strictly coupled. An outcome of these analyses is the Dynameomics Entropy Dictionary, the contents of which have been compared with entropies derived by other theoretical approaches and experiment. As might be expected, the conformational entropies scale linearly with the number of residues, demonstrating that conformational entropy is an extensive property of proteins. The calculated conformational entropies of folding agree well with previous estimates. Detailed analysis of specific cases identifies deviations in conformational entropy from the average values that highlight how conformational entropy varies with sequence, secondary structure, and tertiary fold. Notably, α-helices have lower entropy on average than do β-sheets, and both are lower than coil regions.

Residual Conformational Entropies on the Sugar–Phosphate Backbone of Nucleic Acids: An Analysis of the Nucleosome Core DNA and the Ribosome

A nucleic acid folds according to its free energy, but persistent residual conformational fluctuations remain along its sugar-phosphate backbone even after secondary and tertiary structures have been assembled, and these residual conformational entropies provide a rigorous lower bound for the folding free energy. We extend a recently reported algorithm to calculate the residual backbone entropy along a RNA or DNA given configuration of its bases and apply it to the crystallographic structures of the 23S ribosomal subunit and DNAs in the nucleosome core particle. In the 23S rRNAs, higher entropic strains are concentrated in helices and certain tertiary interaction platforms while residues with high surface accessibility and those not involved in base pairing generally have lower strains. Upon folding, residual backbone entropy in the 23S subunit accounts for an average free energy penalty of +0.47 (kcal/mol)/nt (nt = nucleotide) at 310 K. In nucleosomal DNAs, backbone entropies show periodic oscillations with sequence position correlating with the superhelical twist and shifts in the base-pair-step geometries, and nucleosome positioning on the bound DNA exerts strong influence over where entropic strains are located. In contrast to rRNAs, residual backbone entropies account for a free energy penalty of only +0.09 (kcal/mol)/nt in duplex relative to single-stranded DNAs.

On the relationship between NMR-derived amide order parameters and protein backbone entropy changes.

Molecular dynamics simulations are used to analyze the relationship between NMR-derived squared generalized order parameters of amide NH groups and backbone entropy. Amide order parameters (O(2) NH ) are largely determined by the secondary structure and average values appear unrelated to the overall flexibility of the protein. However, analysis of the more flexible subset (O(2) NH < 0.8) shows that these report both on the local flexibility of the protein and on a different component of the conformational entropy than that reported by the side chain methyl axis order parameters, O(2) axis . A calibration curve for backbone entropy vs. O(2) NH is developed, which accounts for both correlations between amide group motions of different residues, and correlations between backbone and side chain motions. This calibration curve can be used with experimental values of O(2) NH changes obtained by NMR relaxation measurements to extract backbone entropy changes, for example, upon ligand binding. In conjunction with our previous calibration for side chain entropy derived from measured O(2) axis values this provides a prescription for determination of the total protein conformational entropy changes from NMR relaxation measurements.

Conformational Entropy of the RNA Phosphate Backbone and Its Contribution to the Folding Free Energy

While major contributors to the free energy of RNA tertiary structures such as basepairing, base-stacking, and charge and counterion interactions have been studied extensively, little is known about the intrinsic free energy of the backbone. To assess the magnitude of the entropic strains along the phosphate backbone and their impact on the folding free energy, we have formulated a mathematical treatment for computing the volume of the main-chain torsion-angle conformation space between every pair of nucleobases along any sequence to compute the corresponding backbone entropy. To validate this method, we have compared the computed conformational entropies against a statistical free energy analysis of structures in the crystallographic database from several-thousand backbone conformations between nearest-neighbor nucleobases as well as against extensive computer simulations. Using this calculation, we analyzed the backbone entropy of several ribozymes and riboswitches and found that their entropic strains are highly localized along their sequences. The total entropic penalty due to steric congestions in the backbone for the native fold can be as high as 2.4 cal/K/mol per nucleotide for these medium and large RNAs, producing a contribution to the overall free energy of up to 0.72 kcal/mol per nucleotide. For these RNAs, we found that low-entropy high-strain residues are predominantly located at loops with tight turns and at tertiary interaction platforms with unusual structural motifs.

Conformational Entropy of the RNA Phosphate-Sugar Backbone

While major contributors to the free energy of RNA tertiary structures such as base pairing, base stacking and charge and counterion interactions have been studied extensively, little is known about the intrinsic free energy of the backbone. To assess the magnitude of the entropic strains along the phosphate backbone and their impact on the folding free energy, we have formulated a mathematical treatment for computing the volume of the mainchain torsion angle conformation space between every pair of nucleobases along any sequence and the corresponding backbone entropy. We compare the computed conformational entropies against a statistical free energy analysis of structures in the crystallographic database from several thousand backbone conformations between nearest-neighbor nucleobases as well as against computer simulations. Using this calculation we have analyzed the backbone entropy of several ribozymes and riboswitches and found that their entropic strains are highly localized. This suggests the folding and stability of the RNA structures are critically dependent on these local entropic strain domains and leave the rest of the sequence relatively flexible. The total entropic penalty in the backbone for the native fold can be as high as 0.7 cal/K/mol per nucleotide for these medium and large RNAs, producing a contribution to the overall free energy of up to 40 kcal/mol for a 200-nucleotide structure. We also look at the correlation of nucleotide conformations along several sample loop sequences to determine the effect of nucleotides beyond the nearest-neighbor.

Backbone Flexibility of CDR3 and Immune Recognition of Antigens

Conformational entropy is an important component of protein–protein interactions; however, there is no reliable method for computing this parameter. We have developed a statistical measure of residual backbone entropy in folded proteins by using the ϕ–ψ distributions of the 20 amino acids in common secondary structures. The backbone entropy patterns of amino acids within helix, sheet or coil form clusters that recapitulate the branching and hydrogen bonding properties of the side chains in the secondary structure type. The same types of residues in coil and sheet have identical backbone entropies, while helix residues have much smaller conformational entropies. We estimated the backbone entropy change for immunoglobulin complementarity-determining regions (CDRs) from the crystal structures of 34 low-affinity T-cell receptors and 40 high-affinity Fabs as a result of the formation of protein complexes. Surprisingly, we discovered that the computed backbone entropy loss of only the CDR3, but not all CDRs, correlated significantly with the kinetic and affinity constants of the 74 selected complexes. Consequently, we propose a simple algorithm to introduce proline mutations that restrict the conformational flexibility of CDRs and enhance the kinetics and affinity of immunoglobulin interactions. Combining the proline mutations with rationally designed mutants from a previous study led to 2400-fold increase in the affinity of the A6 T-cell receptor for Tax-HLAA2. However, this mutational scheme failed to induce significant binding changes in the already-high-affinity C225–Fab/huEGFR interface. Our results will serve as a roadmap to formulate more effective target functions to design immune complexes with improved biological functions.

Detection and Characterisation of Mutations Responsible for Allele-Specific Protein Thermostabilities at the Mn-Superoxide Dismutase Gene in the Deep-Sea Hydrothermal Vent Polychaete Alvinella pompejana

Alvinella pompejana (Polychaeta, Alvinellidae) is one of the most thermotolerant marine eukaryotes known to date. It inhabits chimney walls of deep-sea hydrothermal vents along the East Pacific Rise (EPR) and is exposed to various challenging conditions (e.g. high temperature, hypoxia and the presence of sulphides, heavy metals and radiations), which increase the production of dangerous reactive oxygen species (ROS). Two different allelic forms of a manganese-superoxide dismutase involved in ROS detoxification, ApMnSOD1 and ApMnSOD2, and differing only by two substitutions (M110L and A138G) were identified in an A. pompejana cDNA library. RFLP screening of 60 individuals from different localities along the EPR showed that ApMnSOD2 was rare (2 %) and only found in the heterozygous state. Dynamic light scattering measurements and residual enzymatic activity experiments showed that the most frequent form (ApMnSOD1) was the most resistant to temperature. Their half-lives were similarly long at 65 °C (>110 min) but exhibited a twofold difference at 80 °C (20.8 vs 9.8 min). Those properties are likely to be explained by the occurrence of an additional sulphur-containing hydrogen bond involving the M110 residue and the effect of the A138 residue on the backbone entropy. Our results confirm the thermophily of A. pompejana and suggest that this locus is a good model to study how the extreme thermal heterogeneity of the vent conditions may help to maintain old rare variants in those populations.

Context and Force Field Dependence of the Loss of Protein Backbone Entropy upon Folding Using Realistic Denatured and Native State Ensembles

The loss of conformational entropy is the largest unfavorable quantity affecting a protein's stability. We calculate the reduction in the number of backbone conformations upon folding using the distribution of backbone dihedral angles (ϕ,ψ) obtained from an experimentally validated denatured state model, along with all-atom simulations for both the denatured and native states. The average loss of entropy per residue is TΔS(BB)(U-N) = 0.7, 0.9, or 1.1 kcal·mol(-1) at T = 298 K, depending on the force field used, with a 0.6 kcal·mol(-1) dispersion across the sequence. The average equates to a decrease of a factor of 3-7 in the number of conformations available per residue (f = Ω(Denatured)/Ω(Native)) or to a total of f(tot) = 3(n)-7(n) for an n residue protein. Our value is smaller than most previous estimates where f = 7-20, that is, our computed TΔS(BB)(U-N) is smaller by 10-100 kcal mol(-1) for n = 100. The differences emerge from our use of realistic native and denatured state ensembles as well as from the inclusion of accurate local sequence preferences, neighbor effects, and correlated motions (vibrations), in contrast to some previous studies that invoke gross assumptions about the entropy in either or both states. We find that the loss of entropy primarily depends on the local environment and less on properties of the native state, with the exception of α-helical residues in some force fields.

Calmodulin Conformational Binding Entropy Is Driven by Transient Salt Bridges

Calmodulin (CaM) is a calcium-binding protein that mediates signal transduction through an ability to differentially bind to highly variable binding sequences in target proteins. Experimental measurements show a linear correlation between conformational entropy and overall binding entropy, suggesting that a predictive ability to calculate conformational entropy is critical to the identification of binding partners and their relative affinities. To identify how binding affects low- and high-frequency motions, and their relationships to conformational entropy, we have employed fully atomistic molecular dynamics simulations with explicit solvent. The calculated quasiharmonic entropies of CaM bound to the targets relative to unbound CaM are in agreement with experimentally derived conformational entropies, where extrapolation to infinite simulation time showed that at least 92% of the backbone entropy is captured over the course of the 100 ns simulation. Enhanced side chain entropy results from conformational disruption of the protein backbone within loop regions between sub-domain elements, acting to facilitate the formation of transient salt bridges between Lys148 at the C-terminus and GLU sidechains both in CaM helix A and GLU side-chains present in selected binding sequences (i.e., endothelial and neuronal nitric oxide synthases). Comparison of Cα root-mean-squared fluctuations and MET sidechain dihedral distributions of CaM bound to each of the targets revealed that MET sidechains have more torsional freedom when the backbone is more flexible, showing that MET sidechain rotation mirrors the overall conformational binding entropy. Taken together, these results help to illuminate the interplay between electrostatic, hydrophobic, backbone and sidechain properties in the ability of CaM to recognize and discriminate against targets by tuning its conformational entropy.

Conformational entropy changes upon lactose binding to the carbohydrate recognition domain of galectin-3

The conformational entropy of proteins can make significant contributions to the free energy of ligand binding. NMR spin relaxation enables site-specific investigation of conformational entropy, via order parameters that parameterize local reorientational fluctuations of rank-2 tensors. Here we have probed the conformational entropy of lactose binding to the carbohydrate recognition domain of galectin-3 (Gal3), a protein that plays an important role in cell growth, cell differentiation, cell cycle regulation, and apoptosis, making it a potential target for therapeutic intervention in inflammation and cancer. We used (15)N spin relaxation experiments and molecular dynamics simulations to monitor the backbone amides and secondary amines of the tryptophan and arginine side chains in the ligand-free and lactose-bound states of Gal3. Overall, we observe good agreement between the experimental and computed order parameters of the ligand-free and lactose-bound states. Thus, the (15)N spin relaxation data indicate that the molecular dynamics simulations provide reliable information on the conformational entropy of the binding process. The molecular dynamics simulations reveal a correlation between the simulated order parameters and residue-specific backbone entropy, re-emphasizing that order parameters provide useful estimates of local conformational entropy. The present results show that the protein backbone exhibits an increase in conformational entropy upon binding lactose, without any accompanying structural changes.

α‐helix stabilization by alanine relative to glycine: Roles of polar and apolar solvent exposures and of backbone entropy

The energetics of alpha-helix formation are fairly well understood and the helix content of a given amino acid sequence can be calculated with reasonable accuracy from helix-coil transition theories that assign to the different residues specific effects on helix stability. In internal helical positions, alanine is regarded as the most stabilizing residue, whereas glycine, after proline, is the more destabilizing. The difference in stabilization afforded by alanine and glycine has been explained by invoking various physical reasons, including the hydrophobic effect and the entropy of folding. Herein, the contribution of these two effects and that of hydrophilic area burial is evaluated by analyzing Ala and Gly mutants implemented in three helices of apoflavodoxin. These data, combined with available data for similar mutations in other proteins (22 Ala/Gly mutations in alpha-helices have been considered), allow estimation of the difference in backbone entropy between alanine and glycine and evaluation of its contribution and that of apolar and polar area burial to the helical stabilization typically associated to Gly-->Ala substitutions. Alanine consistently stabilizes the helical conformation relative to glycine because it buries more apolar area upon folding and because its backbone entropy is lower. However, the relative contribution of polar area burial (which is shown to be destabilizing) and of backbone entropy critically depends on the approximation used to model the structure of the denatured state. In this respect, the excised-peptide model of the unfolded state, proposed by Creamer and coworkers (1995), predicts a major contribution of polar area burial, which is in good agreement with recent quantitations of the relative enthalpic contribution of Ala and Gly residues to alpha-helix formation.

Revisiting the Ramachandran plot: hard-sphere repulsion, electrostatics, and H-bonding in the alpha-helix.

What determines the shape of the allowed regions in the Ramachandran plot? Although Ramachandran explained these regions in terms of 1-4 hard-sphere repulsions, there are discrepancies with the data where, in particular, the alphaR, alphaL, and beta-strand regions are diagonal. The alphaR-region also varies along the alpha-helix where it is constrained at the center and the amino terminus but diffuse at the carboxyl terminus. By analyzing a high-resolution database of protein structures, we find that certain 1-4 hard-sphere repulsions in the standard steric map of Ramachandran do not affect the statistical distributions. By ignoring these steric clashes (NH(i+1) and O(i-1)C), we identify a revised set of steric clashes (CbetaO, O(i-1)N(i+1), CbetaN(i+1), O(i-1)Cbeta, and O(i-1)O) that produce a better match with the data. We also find that the strictly forbidden region in the Ramachandran plot is excluded by multiple steric clashes, whereas the outlier region is excluded by only one significant steric clash. However, steric clashes alone do not account for the diagonal regions. Using electrostatics to analyze the conformational dependence of specific interatomic interactions, we find that the diagonal shape of the alphaR and alphaL-regions also depends on the optimization of the NH(i+1) and O(i-1)C interactions, and the diagonal beta-strand region is due to the alignment of the CO and NH dipoles. Finally, we reproduce the variation of the Ramachandran plot along the alpha-helix in a simple model that uses only H-bonding constraints. This allows us to rationalize the difference between the amino terminus and the carboxyl terminus of the alpha-helix in terms of backbone entropy.

Investigations into Sequence and Conformational Dependence of Backbone Entropy, Inter-basin Dynamics and the Flory Isolated-pair Hypothesis for Peptides

The populations and transitions between Ramachandran basins are studied for combinations of the standard 20 amino acids in monomers, dimers and trimers using an implicit solvent Langevin dynamics algorithm and employing seven commonly used force-fields. Both the basin populations and inter-conversion rates are influenced by the nearest neighbor's conformation and identity, contrary to the Flory isolated-pair hypothesis. This conclusion is robust to the choice of force-field, even though the use of different force-fields produces large variations in the populations and inter-conversion rates between the dominant helical, extended β, and polyproline II basins. The computed variation of conformational and dynamical properties with different force-fields exceeds the difference between explicit and implicit solvent calculations using the same force-field. For all force-fields, the inter-basin transitions exhibit a directional dependence, with most transitions going through extended β conformation, even when it is the least populated basin. The implications of these results are discussed in the context of estimates for the backbone entropy of single residues, and for the ability of all-atom simulations to reproduce experimental protein folding data.

Dynamics and entropy of a calmodulin-peptide complex studied by NMR and molecular dynamics.

All-atom, explicit water molecular dynamics simulations of calcium-loaded calmodulin complexed with a peptide corresponding to the smooth muscle myosin light chain kinase target were carried out at 295 and 346 K. Amide and side chain methyl angular generalized order parameters were calculated and analyzed in the context of the protein's structure and dynamics. The agreement between amide order parameters measured by NMR and those from the simulations was found to be good, especially at the higher temperature, indicating both better convergence for the latter and excellent transferrability of the CHARMM parameters to the higher temperature. Subtle dynamical features such as helix fraying were reproduced. A large range of order parameters for the nine calmodulin methionines was observed in the NMR, and reproduced quite well in the simulations. The major determinant of the methionine order parameter was found to be the proximity to side chains of aromatic residues. An upper bound estimate of the difference in backbone entropy between loop and helical regions was extracted from the order parameters using a model of motion in an effective potential. Although loop regions are more flexible than helical regions, it was found that the entropy loss per residue upon folding was only approximately 20% less for loops than for helices. Pairwise correlated motions, which could significantly lower entropy estimates obtained from order parameter analysis alone, were found to be largely absent.

The role of backbone motions in ligand binding to the c-Src SH3 domain

The Src homology 3 (SH3) domain of pp60 c-src (Src) plays dual roles in signal transduction, through stabilizing the repressed form of the Src kinase and through mediating the formation of activated signaling complexes. Transition of the Src SH3 domain between a variety of binding partners during progression through the cell cycle requires adjustment of a delicate free energy balance. Although numerous structural and functional studies of SH3 have provided an in-depth understanding of structural determinants for binding, the origins of binding energy in SH3-ligand interactions are not fully understood. Considering only the protein-ligand interface, the observed favorable change in standard enthalpy (Δ H = −9.1 kcal/mol) and unfavorable change in standard entropy ( TΔ S = −2.7 kcal/mol) upon binding the proline-rich ligand RLP2 (RALPPLPRY) are inconsistent with the predominantly hydrophobic interaction surface. To investigate possible origins of ligand binding energy, backbone dynamics of free and RLP2-bound SH3 were performed via 15N NMR relaxation and hydrogen-deuterium (H/ 2H) exchange measurements. On the ps-ns time scale, assuming uncorrelated motions, ligand binding results in a significant reduction in backbone entropy (−1.5(±0.6) kcal/mol). Binding also suppresses motions on the μs-ms time scale, which may additionally contribute to an unfavorable change in entropy. A large increase in protection from H/ 2H exchange is observed upon ligand binding, providing evidence for entropy loss due to motions on longer time scales, and supporting the notion that stabilization of pre-existing conformations within a native state ensemble is a fundamental paradigm for ligand binding. Observed changes in motion on all three time scales occur at locations both near and remote from the protein-ligand interface. The propagation of ligand binding interactions across the SH3 domain has potential consequences in target selection through altering both free energy and geometry in intact Src, and suggests that looking beyond the protein-ligand interface is essential in understanding ligand binding energetics.