Phage display identification of high-affinity ligands for human TSG101-UEV: A structural and thermodynamic study of PTAP recognition

The ubiquitin E2 variant domain of TSG101 (TSG101-UEV) plays a pivotal role in protein sorting and virus budding by recognizing PTAP motifs within ubiquitinated proteins. Disrupting TSG101-UEV/PTAP interactions has emerged as a promising strategy for the development of novel host-oriented antivirals with a broad spectrum of action. 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. Here we present a detailed thermodynamic, structural, and dynamic characterization viral PTAP Late domain recognition by TSG101-UEV, combining isothermal titration calorimetry, X-ray diffraction structural studies, molecular dynamics simulations, and computational analysis of intramolecular communication pathways. Our analysis highlights key contributions from conserved hydrophobic contacts and water-mediated hydrogen bonds at the PTAP binding interface. We have identified additional electrostatic hotspots adjacent to the core motif that modulate affinity. Using competitive phage display screening we have improved affinity by 1-2 orders of magnitude, producing novel peptides with low micromolar affinities that combine critical elements found in the best natural binders. Molecular dynamics simulations revealed that optimized peptides engage new pockets on the UEV domain surface. This study provides a comprehensive view of the molecular forces directing TSG101-UEV recognition of PTAP motifs, revealing that binding is governed by conserved structural elements yet tuneable through targeted optimization. These insights open new venues to design inhibitors targeting TSG101-dependent pathways with potential application as novel broad-spectrum antivirals.

Structural dynamics at the active site of the cancer-associated flavoenzyme NQO1 probed by chemical modification with PMSF

A large conformational heterogeneity of human NAD(P)H:quinone oxidoreductase 1 (NQO1), a flavoprotein associated with various human diseases, has been observed to occur in the catalytic site of the enzyme. Here, we report the X-ray structure of NQO1 with phenylmethylsulfonyl fluoride (PMSF) at 1.6 Å resolution. Activity assays confirmed that, despite being covalently bound to the Tyr128 residue at the catalytic site, PMSF did not abolish NQO1 activity. This may indicate that the PMSF molecule does not reduce the high flexibility of Tyr128, thus allowing NADH and DCPIP substrates to bind to the enzyme. Our results show that targeting Tyr128, a key residue in NQO1 function, with small covalently bound molecules could possibly not be a good drug discovery strategy to inhibit this enzyme.

A Systematic Compilation of Human SH3 Domains: A Versatile Superfamily in Cellular Signaling

SRC homology 3 (SH3) domains are fundamental modules that enable the assembly of protein complexes through physical interactions with a pool of proline-rich/noncanonical motifs from partner proteins. They are widely studied modular building blocks across all five kingdoms of life and viruses, mediating various biological processes. The SH3 domains are also implicated in the development of human diseases, such as cancer, leukemia, osteoporosis, Alzheimer's disease, and various infections. A database search of the human proteome reveals the existence of 298 SH3 domains in 221 SH3 domain-containing proteins (SH3DCPs), ranging from 13 to 720 kilodaltons. A phylogenetic analysis of human SH3DCPs based on their multi-domain architecture seems to be the most practical way to classify them functionally, with regard to various physiological pathways. This review further summarizes the achievements made in the classification of SH3 domain functions, their binding specificity, and their significance for various diseases when exploiting SH3 protein modular interactions as drug targets.

A calorimetric and structural analysis of cooperativity in the thermal unfolding of the PDZ tandem of human Syntenin-1

Syntenin-1 is a multidomain protein containing a central tandem of two PDZ domains flanked by two unnamed domains. Previous structural and biophysical studies show that the two PDZ domains are functional both isolated and in tandem, occurring a gain in their respective binding affinities when joined through its natural short linker. To get insight into the molecular and energetic reasons of such a gain, here, the first thermodynamic characterization of the conformational equilibrium of Syntenin-1 is presented, with special focus on its PDZ domains. These studies include the thermal unfolding of the whole protein, the PDZ-tandem construct and the two isolated PDZ domains using circular dichroism, differential scanning fluorimetry and differential scanning calorimetry. The isolated PDZ domains show low stability (ΔG < 10 kJ·mol-1) and poor cooperativity compared to the PDZ-tandem, which shows higher stability (20-30 kJ·mol-1) and a fully cooperative behaviour, with energetics similar to that previously described for archetypical PDZ domains. The high-resolution structures suggest that this remarkable increase in cooperativity is associated to strong, water-mediated, interactions at the interface between the PDZ domains, associated to nine conserved hydration regions. The low Tm value (45 °C), the anomalously high unfolding enthalpy (>400 kJ·mol-1), and native heat capacity values (above 40 kJ·K-1·mol-1), indicate that these interfacial buried waters play a relevant role in Syntenin-1 folding energetics.

Phage display identification of nanomolar ligands for human NEDD4-WW3: Energetic and dynamic implications for the development of broad-spectrum antivirals

The recognition of PPxY viral Late domains by the third WW domain of the human HECT-E3 ubiquitin ligase NEDD4 (NEDD4-WW3) is essential for the budding of many viruses. Blocking these interactions is a promising strategy to develop broad-spectrum antivirals. As all WW domains, NEDD4-WW3 is a challenging therapeutic target due to the low binding affinity of its natural interactions, its high conformational plasticity, and its complex thermodynamic behavior. In this work, we set out to investigate whether high affinity can be achieved for monovalent ligands binding to the isolated NEDD4-WW3 domain. We show that a competitive phage-display set-up allows for the identification of high-affinity peptides showing inhibitory activity of viral budding. A detailed biophysical study combining calorimetry, nuclear magnetic resonance, and molecular dynamic simulations reveals that the improvement in binding affinity does not arise from the establishment of new interactions with the domain, but is associated to conformational restrictions imposed by a novel C-terminal -LFP motif in the ligand, unprecedented in the PPxY interactome. These results, which highlight the complexity of WW domain interactions, provide valuable insight into the key elements for high binding affinity, of interest to guide virtual screening campaigns for the identification of novel therapeutics targeting NEDD4-WW3 interactions.

Understanding binding affinity and specificity of modular protein domains: A focus in ligand design for the polyproline-binding families

Within the modular protein domains there are five families that recognize proline-rich sequences: SH3, WW, EVH1, GYF and UEV domains. This chapter reviews the main strategies developed for the design of ligands for these families, including peptides, peptidomimetics and drugs. We also describe some studies aimed to understand the molecular reasons responsible for the intrinsic affinity and specificity of these domains.

Rapid Quantification of Protein-Ligand Binding via 19F NMR Lineshape Analysis

Fluorine incorporation is ideally suited to many NMR techniques, and incorporation of fluorine into proteins and fragment libraries for drug discovery has become increasingly common. Here, we use one-dimensional ¹⁹F NMR lineshape analysis to quantify the kinetics and equilibrium thermodynamics for the binding of a fluorine-labeled Src homology 3 (SH3) protein domain to four proline-rich peptides. SH3 domains are one of the largest and most well-characterized families of protein recognition domains and have a multitude of functions in eukaryotic cell signaling. First, we showe that fluorine incorporation into SH3 causes only minor structural changes to both the free and bound states using amide proton temperature coefficients. We then compare the results from lineshape analysis of one-dimensional ¹⁹F spectra to those from two-dimensional ¹H-¹⁵N heteronuclear single quantum coherence spectra. Their agreement demonstrates that one-dimensional ¹⁹F lineshape analysis is a robust, low-cost, and fast alternative to traditional heteronuclear single quantum coherence-based experiments. The data show that binding is diffusion limited and indicate that the transition state is highly similar to the free state. We also measured binding as a function of temperature. At equilibrium, binding is enthalpically driven and arises from a highly positive activation enthalpy for association with small entropic contributions. Our results agree with those from studies using different techniques, providing additional evidence for the utility of ¹⁹F NMR lineshape analysis, and we anticipate that this analysis will be an effective tool for rapidly characterizing the energetics of protein interactions.

Entropy Hotspots for the Binding of Intrinsically Disordered Ligands to a Receptor Domain

Proline-rich motifs (PRMs) are widely used for mediating protein-protein interactions with weak binding affinities. Because they are intrinsically disordered when unbound, conformational entropy plays a significant role for the binding. However, residue-level differences of the entropic contribution in the binding of different ligands remain not well understood. We use all-atom molecular dynamics simulation and the maximal information spanning tree formalism to analyze conformational entropy associated with the binding of two PRMs, one from the Abl kinase and the other from the nonstructural protein 1 of the 1918 Spanish influenza A virus, to the N-terminal SH3 (nSH3) domain of the CrkII protein. Side chains of the stably folded nSH3 experience more entropy change upon ligand binding than the backbone, whereas PRMs involve comparable but heterogeneous entropy changes among the backbone and side chains. In nSH3, two conserved nonpolar residues forming contacts with the PRM experience the largest side-chain entropy loss. In contrast, the C-terminal charged residues of PRMs that form polar contacts with nSH3 experience the greatest side-chain entropy loss, although their "fuzzy" nature is attributable to the backbone that remains relatively flexible. Thus, residues that form high-occupancy contacts between nSH3 and PRM do not reciprocally contribute to entropy loss. Furthermore, certain surface residues of nSH3 distal to the interface with PRMs gain entropy, indicating a nonlocal effect of ligand binding. Comparing between the PRMs from cAbl and nonstructural protein 1, the latter involves a larger side-chain entropy loss and forms more contacts with nSH3. Consistent with experiments, this indicates stronger binding of the viral ligand at the expense of losing the flexibility of side chains, whereas the backbone experiences less entropy loss. The entropy "hotspots" as identified in this study will be important for tuning the binding affinity of various ligands to a receptor.

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.

Site-specific 2D IR spectroscopy: a general approach for the characterization of protein dynamics with high spatial and temporal resolution

The conformational heterogeneity and dynamics of protein side chains contribute to function, but investigating exactly how is hindered by experimental challenges arising from the fast timescales involved and the spatial heterogeneity of protein structures. The potential of two-dimensional infrared (2D IR) spectroscopy for measuring conformational heterogeneity and dynamics with unprecedented spatial and temporal resolution has motivated extensive effort to develop amino acids with functional groups that have frequency-resolved absorptions to serve as probes of their protein microenvironments. We demonstrate the full advantage of the approach by selective incorporation of the probe p-cyanophenylalanine at six distinct sites in a Src homology 3 domain and the application of 2D IR spectroscopy to site-specifically characterize heterogeneity and dynamics and their contribution to cognate ligand binding. The approach revealed a wide range of microenvironments and distinct responses to ligand binding, including at the three adjacent, conserved aromatic residues that form the recognition surface of the protein. Molecular dynamics simulations performed for all the labeled proteins provide insight into the underlying heterogeneity and dynamics. Similar application of 2D IR spectroscopy and site-selective probe incorporation will allow for the characterization of heterogeneity and dynamics of other proteins, how heterogeneity and dynamics are affected by solvation and local structure, and how they might contribute to biological function.

Binding Mechanism of the N-Terminal SH3 Domain of CrkII and Proline-Rich Motifs in cAbl

The N-terminal Src homology 3 (nSH3) domain of a signaling adaptor protein, CT-10 regulator of kinase II (CrkII), recognizes proline-rich motifs (PRMs) of binding partners, such as cAbl kinase. The interaction between CrkII and cAbl kinase is involved in the regulation of cell spreading, microbial pathogenesis, and cancer metastasis. Here, we report the detailed biophysical characterizations of the interactions between the nSH3 domain of CrkII and PRMs in cAbl. We identified that the nSH3 domain of CrkII binds to three PRMs in cAbl with virtually identical affinities. Structural studies, by using x-ray crystallography and NMR spectroscopy, revealed that the binding modes of all three nSH3:PRM complexes are highly similar to each other. Van 't Hoff analysis revealed that nSH3:PRM interaction is associated with favorable enthalpy and unfavorable entropy change. The combination of experimentally determined thermodynamic parameters, structure-based calculations, and (15)N NMR relaxation analysis highlights the energetic contribution of conformational entropy change upon the complex formation, and water molecules structured in the binding interface of the nSH3:PRM complex. Understanding the molecular basis of nSH3:PRM interaction will provide, to our knowledge, new insights for the rational design of small molecules targeting the interaction between CrkII and cAbl.

Grb2 carboxyl-terminal SH3 domain can bivalently associate with two ligands, in an SH3 dependent manner

Src homology domain containing leukocyte protein of 65 kDa (SLP65), the growth factor receptor binding protein 2 (Grb2), and the guanine nucleotide exchange factor for the Rho family GTPases (Vav), self associate in unstimulated B cells as components of the preformed B cell receptor transducer module, in an SH3-dependent manner. The complex enables the B cell to promptly respond to BCR aggregation, resulting in signal amplification. It also facilitates Vav translocation to the membrane rafts, for activation. Here we uncover the molecular mechanism by which the complex may be formed in the B cell. The C-terminal SH3 domain (SH3C) of Grb2 bivalently interacts with the atypical non-PxxP proline rich region of SLP65, and the N-terminal SH3 domain (SH3N) of Vav, both the interactions crucial for the proper functioning of the B cell. Most surprisingly, the two ligands bind the same ligand binding site on the surface of Grb2 SH3C. Addition of SLP65 peptide to the Grb2-Vav complex abrogates the interaction completely, displacing Vav. However, the addition of Vav SH3N to the SLP65-Grb2 binary complex, results in a trimeric complex. Extrapolating these results to the in vivo conditions, Grb2 should bind the SLP65 transducer module first, and then Vav should associate.

The role of water molecules in the binding of class I and II peptides to the SH3 domain of the Fyn tyrosine kinase

Interactions of proline-rich motifs with SH3 domains are present in signal transduction and other important cell processes. Analysis of structural and thermodynamic data suggest a relevant role of water molecules in these protein-protein interactions. To determine whether or not the SH3 domain of the Fyn tyrosine kinase shows the same behaviour, the crystal structures of its complexes with two high-affinity synthetic peptides, VSL12 and APP12, which are class I and II peptides, respectively, have been solved. In the class I complexes two water molecules were found at the binding interface that were not present in the class II complexes. The structures suggest a role of these water molecules in facilitating conformational changes in the SH3 domain to allow the binding of the class I or II peptides. In the third binding pocket these changes modify the cation-π and salt-bridge interactions that determine the affinity of the binding. Comparison of the water molecules involved in the binding of the peptides with previous reported hydration spots suggests a different pattern for the SH3 domains of the Src tyrosine kinase family.

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.

Structure of the c-Src-SH3 domain in complex with a proline-rich motif of NS5A protein from the hepatitis C virus

The non-structural hepatitis C virus proteins NS5A and NS5B form a complex through interaction with the SH2 and SH3 domains of the non-receptor Src tyrosine kinase, which seems essential for viral replication. We have crystallized the complex between the SH3 domain of the c-Src tyrosine kinase and the C-terminal proline rich motif of the NS5A protein (A349PPIPPPRRKR359). Crystals obtained at neutral pH belong to the space group I41, with a single molecule of the SH3/NS5A complex at the asymmetric unit. The NS5A peptide is bound in a reverse orientation (class II) and the comparison of this structure with those of the high affinity synthetic peptides APP12 and VSL12 shows some important differences at the salt bridge that drives the peptide orientation. Further conformational changes in residues placed apart from the binding site also seem to play an important role in the binding orientation of this peptide. Our results show the interaction of the SH3 domain of the c-Src tyrosine kinase with a proline rich motif in the NS5A protein and point to their potential interaction in vivo.

Molecular origin of the binding of WWOX tumor suppressor to ErbB4 receptor tyrosine kinase

The ability of WWOX tumor suppressor to physically associate with the intracellular domain (ICD) of ErbB4 receptor tyrosine kinase is believed to play a central role in downregulating the transcriptional function of the latter. Herein, using various biophysical methods, we show that while the WW1 domain of WWOX binds to PPXY motifs located within the ICD of ErbB4 in a physiologically relevant manner, the WW2 domain does not. Importantly, while the WW1 domain absolutely requires the integrity of the PPXY consensus sequence, nonconsensus residues within and flanking this motif do not appear to be critical for binding. This strongly suggests that the WW1 domain of WWOX is rather promiscuous toward its cellular partners. We also provide evidence that the lack of binding of the WW2 domain of WWOX to PPXY motifs is due to the replacement of a signature tryptophan, lining the hydrophobic ligand binding groove, with tyrosine (Y85). Consistent with this notion, the Y85W substitution within the WW2 domain exquisitely restores its binding to PPXY motifs in a manner akin to the binding of the WW1 domain of WWOX. Of particular significance is the observation that the WW2 domain augments the binding of the WW1 domain to ErbB4, implying that the former serves as a chaperone within the context of the WW1-WW2 tandem module of WWOX in agreement with our findings reported previously. Altogether, our study sheds new light on the molecular basis of an important WW-ligand interaction involved in mediating a plethora of cellular processes.

DFT and docking studies of rhodostreptomycins A and B and their interactions with solvated/nonsolvated Mg²⁺ and Ca²⁺ ions

The interactions of L-aminoglucosidic stereoisomers such as rhodostreptomycins A (Rho A) and B (Rho B) with cations (Mg(2+), Ca(2+), and H(+)) were studied by a quantum mechanical method that utilized DFT with B3LYP/6-311G. Docking studies were also carried out in order to explore the surface recognition properties of L-aminoglucoside with respect to Mg(2+) and Ca(2+) ions under solvated and nonsolvated conditions. Although both of the stereoisomers possess similar physicochemical/antibiotic properties against Helicobacter pylori, the thermochemical values for these complexes showed that its high affinity for Mg(2+) cations caused the hydration of Rho B. According to the results of the calculations, for Rho A-Ca(2+)(H2O)6, ΔH = -72.21 kcal mol(-1); for Rho B-Ca(2+)(H2O)6, ΔH = -72.53 kcal mol(-1); for Rho A-Mg(2+)(H2O)6, ΔH = -72.99 kcal mol(-1) and for Rho B-Mg(2+)(H2O)6, ΔH = -95.00 kcal mol(-1), confirming that Rho B binds most strongly with hydrated Mg(2+), considering the energy associated with this binding process. This result suggests that Rho B forms a more stable complex than its isomer does with magnesium ion. Docking results show that both of these rhodostreptomycin molecules bind to solvated Ca(2+) or Mg(2+) through hydrogen bonding. Finally, Rho B is more stable than Rho A when protonation occurs.

Invited review: Caseins and the casein micelle: their biological functions, structures, and behavior in foods

A typical casein micelle contains thousands of casein molecules, most of which form thermodynamically stable complexes with nanoclusters of amorphous calcium phosphate. Like many other unfolded proteins, caseins have an actual or potential tendency to assemble into toxic amyloid fibrils, particularly at the high concentrations found in milk. Fibrils do not form in milk because an alternative aggregation pathway is followed that results in formation of the casein micelle. As a result of forming micelles, nutritious milk can be secreted and stored without causing either pathological calcification or amyloidosis of the mother's mammary tissue. The ability to sequester nanoclusters of amorphous calcium phosphate in a stable complex is not unique to caseins. It has been demonstrated using a number of noncasein secreted phosphoproteins and may be of general physiological importance in preventing calcification of other biofluids and soft tissues. Thus, competent noncasein phosphoproteins have similar patterns of phosphorylation and the same type of flexible, unfolded conformation as caseins. The ability to suppress amyloid fibril formation by forming an alternative amorphous aggregate is also not unique to caseins and underlies the action of molecular chaperones such as the small heat-shock proteins. The open structure of the protein matrix of casein micelles is fragile and easily perturbed by changes in its environment. Perturbations can cause the polypeptide chains to segregate into regions of greater and lesser density. As a result, the reliable determination of the native structure of casein micelles continues to be extremely challenging. The biological functions of caseins, such as their chaperone activity, are determined by their composition and flexible conformation and by how the casein polypeptide chains interact with each other. These same properties determine how caseins behave in the manufacture of many dairy products and how they can be used as functional ingredients in other foods.

Atomic resolution structures of the c-Src SH3 domain in complex with two high-affinity peptides from classes I and II

The atomic resolution crystal structures of complexes between the SH3 domain of the c-Src tyrosine kinase and two high-affinity peptides belonging to class I and class II have been solved. The crystals of the Thr98Asp and Thr98Glu mutants in complex with the APP12 peptide (APPLPPRNRPRL) belonged to the trigonal space group P3121 and in both cases the asymmetric unit was composed of one molecule of the SH3-APP12 complex. The crystals of the Thr98Glu mutant in complex with the VSL12 peptide (VSLARRPLPLP) belonged to the trigonal space group P3221 and the asymmetric unit was also composed of a single molecule of the SH3-VSL12 complex. All crystals were obtained in the presence of PEG 300 under the same conditions as reported for the intertwined dimeric structure of the c-Src SH3 domain, but the presence of the peptide stabilizes the monomeric form of the domain. These structures allow a detailed analysis of the role of salt bridges, cation-π interactions and hydrogen bonds in the binding of proline-rich motifs to the c-Src SH3 domain. Moreover, these crystallographic structures allow the role of water molecules in the binding of these motifs to the c-Src SH3 domain to be studied for the first time.

Proline: the distribution, frequency, positioning, and common functional roles of proline and polyproline sequences in the human proteome

Proline is an anomalous amino acid. Its nitrogen atom is covalently locked within a ring, thus it is the only proteinogenic amino acid with a constrained phi angle. Sequences of three consecutive prolines can fold into polyproline helices, structures that join alpha helices and beta pleats as architectural motifs in protein configuration. Triproline helices are participants in protein-protein signaling interactions. Longer spans of repeat prolines also occur, containing as many as 27 consecutive proline residues. Little is known about the frequency, positioning, and functional significance of these proline sequences. Therefore we have undertaken a systematic bioinformatics study of proline residues in proteins. We analyzed the distribution and frequency of 687,434 proline residues among 18,666 human proteins, identifying single residues, dimers, trimers, and longer repeats. Proline accounts for 6.3% of the 10,882,808 protein amino acids. Of all proline residues, 4.4% are in trimers or longer spans. We detected patterns that influence function based on proline location, spacing, and concentration. We propose a classification based on proline-rich, polyproline-rich, and proline-poor status. Whereas singlet proline residues are often found in proteins that display recurring architectural patterns, trimers or longer proline sequences tend be associated with the absence of repetitive structural motifs. Spans of 6 or more are associated with DNA/RNA processing, actin, and developmental processes. We also suggest a role for proline in Kruppel-type zinc finger protein control of DNA expression, and in the nucleation and translocation of actin by the formin complex.

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.