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.

A comparison of the location in membranes of curcumin and curcumin-derived bivalent compounds with potential neuroprotective capacity for Alzheimer's disease.

Curcumin and two bivalent compounds, namely 17MD and 21MO, both obtained by conjugation of curcumin with a steroid molecule that acts as a membrane anchor, were comparatively studied. When incorporated into 1,2-dipalmitoyl-sn-glycero-3-phosphocholine the compounds showed a very limited solubility in the model membranes. Curcumin and the two bivalent compounds were also incorporated in membranes of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and quenching the fluorescence of pure curcumin or of the curcumin moiety in the bivalent compounds by acrylamide it was seen that curcumin was accessible to this water soluble quencher but the molecule was somehow located in a hydrophobic environment. This was confirmed by quenching with doxyl-phosphatidylcholines, indicating that the curcumin moieties of 17MD and 21MO were in a more polar environment than pure curcumin itself. 1H NOESY MAS-NMR analysis supports this notion by showing that the orientation of curcumin was parallel to the plane of the membrane surface close to C2 and C3 of the fatty acyl chains, while the curcumin moiety of 17MD and 21MO positioned close to the polar part of the membrane with the steroid moiety in the centre of the membrane. Molecular dynamics studies were in close agreement with the experimental results with respect to the likely proximity of the protons studied by NMR and show that 17MD and 21MO have a clear tendency to aggregate in a fluid membrane. The anchorage of the bivalent compounds to the membrane leaving the curcumin moiety near the polar part may be very important to facilitate the bioactivity of the curcumin moiety when used as anti-Alzheimer drugs.

Non-conservation of folding rates in the thioredoxin family reveals degradation of ancestral unassisted-folding.

Evolution involves not only adaptation, but also the degradation of superfluous features. Many examples of degradation at the morphological level are known (vestigial organs, for instance). However, the impact of degradation on molecular evolution has been rarely addressed. Thioredoxins serve as general oxidoreductases in all cells. Here, we report extensive mutational analyses on the folding of modern and resurrected ancestral bacterial thioredoxins. Contrary to claims from recent literature, in vitro folding rates in the thioredoxin family are not evolutionarily conserved, but span at least a ∼100-fold range. Furthermore, modern thioredoxin folding is often substantially slower than ancestral thioredoxin folding. Unassisted folding, as probed in vitro, thus emerges as an ancestral vestigial feature that underwent degradation, plausibly upon the evolutionary emergence of efficient cellular folding assistance. More generally, our results provide evidence that degradation of ancestral features shapes, not only morphological evolution, but also the evolution of individual proteins.

De novo active sites for resurrected Precambrian enzymes.

Protein engineering studies often suggest the emergence of completely new enzyme functionalities to be highly improbable. However, enzymes likely catalysed many different reactions already in the last universal common ancestor. Mechanisms for the emergence of completely new active sites must therefore either plausibly exist or at least have existed at the primordial protein stage. Here, we use resurrected Precambrian proteins as scaffolds for protein engineering and demonstrate that a new active site can be generated through a single hydrophobic-to-ionizable amino acid replacement that generates a partially buried group with perturbed physico-chemical properties. We provide experimental and computational evidence that conformational flexibility can assist the emergence and subsequent evolution of new active sites by improving substrate and transition-state binding, through the sampling of many potentially productive conformations. Our results suggest a mechanism for the emergence of primordial enzymes and highlight the potential of ancestral reconstruction as a tool for protein engineering.

Mechanism of ß-actin mRNA Recognition by ZBP1.

Zipcode binding protein 1 (ZBP1) is an oncofetal RNA-binding protein that mediates the transport and local translation of ß-actin mRNA by the KH3-KH4 di-domain, which is essential for neuronal development. The high-resolution structures of KH3-KH4 with their respective target sequences show that KH4 recognizes a non-canonical GGA sequence via an enlarged and dynamic hydrophobic groove, whereas KH3 binding to a core CA sequence occurs with low specificity. A data-informed kinetic simulation of the two-step binding reaction reveals that the overall reaction is driven by the second binding event and that the moderate affinities of the individual interactions favor RNA looping. Furthermore, the concentration of ZBP1, but not of the target RNA, modulates the interaction, which explains the functional significance of enhanced ZBP1 expression during embryonic development.

Post-translational modifications modulate ligand recognition by the third PDZ domain of the MAGUK protein PSD-95.

The relative promiscuity of hub proteins such as postsynaptic density protein-95 (PSD-95) can be achieved by alternative splicing, allosteric regulation, and post-translational modifications, the latter of which is the most efficient method of accelerating cellular responses to environmental changes in vivo. Here, a mutational approach was used to determine the impact of phosphorylation and succinimidation post-translational modifications on the binding affinity of the postsynaptic density protein-95/discs large/zonula occludens-1 (PDZ3) domain of PSD-95. Molecular dynamics simulations revealed that the binding affinity of this domain is influenced by an interplay between salt-bridges linking the α3 helix, the ß2-ß3 loop and the positively charged Lys residues in its high-affinity hexapeptide ligand KKETAV. The α3 helix is an extra structural element that is not present in other PDZ domains, which links PDZ3 with the following SH3 domain in the PSD-95 protein. This regulatory mechanism was confirmed experimentally via thermodynamic and NMR chemical shift perturbation analyses, discarding intra-domain long-range effects. Taken together, the results presented here reveal the molecular basis of the regulatory role of the α3 extra-element and the effects of post-translational modifications of PDZ3 on its binding affinity, both energetically and dynamically.

The interconversion between a flexible ß-sheet and a fibril ß-arrangement constitutes the main conformational event during misfolding of PSD95-PDZ3 domain.

The temperature-induced misfolding pathway of PDZ3, the third PDZ domain of the PSD95 neuronal protein, is populated by a trimeric ß-sheet-rich intermediate state that leads to a stepwise and reversible formation of supramacromolecular structures. Using FTIR, we have found that misfolding of this pathway is not due to different ensembles of a variety of precursors, but comes mainly from the interconversion of a flexible ß-sheet of the domain to wormlike fibrils. The appearance of the wormlike fibril FTIR component is also accompanied by a slight decrease of the band that corresponds to loops in the native state, whereas the rest of the regular elements of secondary structure are fairly well maintained upon misfolding. Transmission electron microscope micrographs have confirmed the presence of wormlike fibrils upon heating at 60°C, where the trimeric intermediate is maximally populated. Toxicity assays in the human neuroblastoma cell line SH-SY5Y show that cytotoxicity increases as the aggregation pathway proceeds. NMR analysis of chemical shifts as a function of temperature has revealed, as one of the main conformational aspects of such an interconversion at the residue level, that the ß-sheet arrangement around strand ß3 promotes the change that drives misfolding of the PDZ3 domain.

KH domains with impaired nucleic acid binding as a tool for functional analysis.

In eukaryotes, RNA-binding proteins that contain multiple K homology (KH) domains play a key role in coordinating the different steps of RNA synthesis, metabolism and localization. Understanding how the different KH modules participate in the recognition of the RNA targets is necessary to dissect the way these proteins operate. We have designed a KH mutant with impaired RNA-binding capability for general use in exploring the role of individual KH domains in the combinatorial functional recognition of RNA targets. A double mutation in the hallmark GxxG loop (GxxG-to-GDDG) impairs nucleic acid binding without compromising the stability of the domain. We analysed the impact of the GDDG mutations in individual KH domains on the functional properties of KSRP as a prototype of multiple KH domain-containing proteins. We show how the GDDG mutant can be used to directly link biophysical information on the sequence specificity of the different KH domains of KSRP and their role in mRNA recognition and decay. This work defines a general molecular biology tool for the investigation of the function of individual KH domains in nucleic acid binding proteins.

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.

Analysis of the thermodynamics of binding of an SH3 domain to proline-rich peptides using a chimeric fusion protein.

A complete understanding of the thermodynamic determinants of binding between SH3 domains and proline-rich peptides is crucial to the development of rational strategies for designing ligands for these important domains. Recently we engineered a single-chain chimeric protein by fusing the alpha-spectrin Src homology region 3 (SH3) domain to the decapeptide APSYSPPPPP (p41). This chimera mimics the structural and energetic features of the interaction between SH3 domains and proline-rich peptides. Here we show that analysing the unfolding thermodynamics of single-point mutants of this chimeric fusion protein constitutes a very useful approach to deciphering the thermodynamics of SH3-ligand interactions. To this end, we investigated the contribution of each proline residue of the ligand sequence to the SH3-peptide interaction by producing six single Pro-Ala mutants of the chimeric protein and analysing their unfolding thermodynamics by differential scanning calorimetry (DSC). Structural analyses of the mutant chimeras by circular dichroism, fluorescence and NMR together with NMR-relaxation measurements indicate conformational flexibility at the binding interface, which is strongly affected by the different Pro-Ala mutations. An analysis of the DSC thermograms on the basis of a three-state unfolding model has allowed us to distinguish and separate the thermodynamic magnitudes of the interaction at the binding interface. The model assumes equilibrium between the "unbound" and "bound" states at the SH3-peptide binding interface. The resulting thermodynamic magnitudes classify the different proline residues according to their importance in the interaction as P2 approximately P7 approximately P10>P9 approximately P6>P8, which agrees well with Lim's model for the interaction between SH3 domains and proline-rich peptides. In addition, the thermodynamic signature of the interaction is the same as that usually found for this type of binding, with a strong enthalpy-entropy compensation for all the mutants. This compensation appears to derive from an increase in conformational flexibility concomitant to the weakening of the interactions at the binding interface. We conclude that our approach, based on DSC and site-directed mutagenesis analysis of chimeric fusion proteins, may serve as a suitable tool to analyse the energetics of weak biomolecular interactions such as those involving SH3 domains.

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.

The high-resolution NMR structure of a single-chain chimeric protein mimicking a SH3-peptide complex.

Here we present the high-resolution NMR structure of a chimera (SPCp41) between alpha-spectrin SH3 domain and the decapeptide p41. The tertiary structure mimics perfectly the interactions typically found in SH3-peptide complexes and is remarkably similar to that of the complex between the separate Spc-SH3 domain and ligand p41. Relaxation data confirm the tight binding between the ligand and SH3 part of the chimera. This chimera will serve as a tool for a deeper understanding of the relationship between structure and thermodynamics of binding using a combination of NMR, stability and site-directed mutagenesis studies, which can lead to an effective strategy for ligand design.