Hydrogen bonding heterogeneity correlates with protein folding transition state passage time as revealed by data sonification.

Protein-protein and protein-water hydrogen bonding interactions play essential roles in the way a protein passes through the transition state during folding or unfolding, but the large number of these interactions in molecular dynamics (MD) simulations makes them difficult to analyze. Here, we introduce a state space representation and associated "rarity" measure to identify and quantify transition state passage (transit) events. Applying this representation to a long MD simulation trajectory that captured multiple folding and unfolding events of the GTT WW domain, a small protein often used as a model for the folding process, we identified three transition categories: Highway (faster), Meander (slower), and Ambiguous (intermediate). We developed data sonification and visualization tools to analyze hydrogen bond dynamics before, during, and after these transition events. By means of these tools, we were able to identify characteristic hydrogen bonding patterns associated with "Highway" versus "Meander" versus "Ambiguous" transitions and to design algorithms that can identify these same folding pathways and critical protein-water interactions directly from the data. Highly cooperative hydrogen bonding can either slow down or speed up transit. Furthermore, an analysis of protein-water hydrogen bond dynamics at the surface of WW domain shows an increase in hydrogen bond lifetime from folded to unfolded conformations with Ambiguous transitions as an outlier. In summary, hydrogen bond dynamics provide a direct window into the heterogeneity of transits, which can vary widely in duration (by a factor of 10) due to a complex energy landscape.

Elucidation of WW domain ligand binding specificities in the Hippo pathway reveals STXBP4 as YAP inhibitor.

The Hippo pathway, which plays a critical role in organ size control and cancer, features numerous WW domain-based protein-protein interactions. However, ~100 WW domains and 2,000 PY motif-containing peptide ligands are found in the human proteome, raising a "WW-PY" binding specificity issue in the Hippo pathway. In this study, we have established the WW domain binding specificity for Hippo pathway components and uncovered a unique amino acid sequence required for it. By using this criterion, we have identified a WW domain-containing protein, STXBP4, as a negative regulator of YAP. Mechanistically, STXBP4 assembles a protein complex comprising α-catenin and a group of Hippo PY motif-containing components/regulators to inhibit YAP, a process that is regulated by actin cytoskeleton tension. Interestingly, STXBP4 is a potential tumor suppressor for human kidney cancer, whose downregulation is correlated with YAP activation in clear cell renal cell carcinoma. Taken together, our study not only elucidates the WW domain binding specificity for the Hippo pathway, but also reveals STXBP4 as a player in actin cytoskeleton tension-mediated Hippo pathway regulation.

Relationship of Thermostability and Binding Affinity in Metal-binding WW-Domain Minireceptors.

The design of metallo-miniproteins advances our understanding of the structural and functional roles of metals in proteins. We recently designed a metal-binding WW domain, WW-CA-Nle, which displays three histidine residues on its surface for coordination of divalent metals Ni(II), Zn(II) and Cu(II). However, WW-CA-Nle is a molten globule in the apo state and thus showed only moderate binding affinities with Kd values in the µM regime. In this report, we hypothesize that improved thermal stability of the apo state of the metal binding WW-domain scaffold should lead to improved preorganization of the metal-binding site and consequently to higher metal-binding affinities. By redesigning WW-CA-Nle, we obtained WW-CA variants, WW-CA-min and WW-CA-ANG, which were fully folded in the apo states and displayed moderate to excellent thermostabilities in the apo and holo states. We were able to show that the improved thermal stabilities led to improved metal binding, which was reflected in Kd values that were at least one order of magnitude lower compared to WW-CA-Nle. EPR spectroscopy and ITC measurements revealed a better defined and predisposed metal binding site in WW-CA-ANG.

NMR Dynamic View of the Destabilization of WW4 Domain by Chaotropic GdmCl and NaSCN.

GdmCl and NaSCN are two strong chaotropic salts commonly used in protein folding and stability studies, but their microscopic mechanisms remain enigmatic. Here, by CD and NMR, we investigated their effects on conformations, stability, binding and backbone dynamics on ps-ns and µs-ms time scales of a 39-residue but well-folded WW4 domain at salt concentrations ≤200 mM. Up to 200 mM, both denaturants did not alter the tertiary packing of WW4, but GdmCl exerted more severe destabilization than NaSCN. Intriguingly, GdmCl had only weak binding to amide protons, while NaSCN showed extensive binding to both hydrophobic side chains and amide protons. Neither denaturant significantly affected the overall ps-ns backbone dynamics, but they distinctively altered µs-ms backbone dynamics. This study unveils that GdmCl and NaSCN destabilize a protein before the global unfolding occurs with differential binding properties and µs-ms backbone dynamics, implying the absence of a simple correlation between thermodynamic stability and backbone dynamics of WW4 at both ps-ns and µs-ms time scales.

Structural insights into the role of the WW2 domain on tandem WW-PPxY motif interactions of oxidoreductase WWOX.

Class I WW domains are present in many proteins of various functions and mediate protein interactions by binding to short linear PPxY motifs. Tandem WW domains often bind peptides with multiple PPxY motifs, but the interplay of WW-peptide interactions is not always intuitive. The WW domain-containing oxidoreductase (WWOX) harbors two WW domains: an unstable WW1 capable of PPxY binding and stable WW2 that cannot bind PPxY. The WW2 domain has been suggested to act as a WW1 domain chaperone, but the underlying mechanism of its chaperone activity remains to be revealed. Here, we combined NMR, isothermal calorimetry, and structural modeling to elucidate the roles of both WW domains in WWOX binding to its PPxY-containing substrate ErbB4. Using NMR, we identified an interaction surface between these two domains that supports a WWOX conformation compatible with peptide substrate binding. Isothermal calorimetry and NMR measurements also indicated that while binding affinity to a single PPxY motif is marginally increased in the presence of WW2, affinity to a dual-motif peptide increases 10-fold. Furthermore, we found WW2 can directly bind double-motif peptides using its canonical binding site. Finally, differential binding of peptides in mutagenesis experiments was consistent with a parallel N- to C-terminal PPxY tandem motif orientation in binding to the WW1-WW2 tandem domain, validating structural models of the interaction. Taken together, our results reveal the complex nature of tandem WW-domain organization and substrate binding, highlighting the contribution of WWOX WW2 to both protein stability and target binding.

Nedd4-2 binding to 14-3-3 modulates the accessibility of its catalytic site and WW domains.

Neural precursor cells expressed developmentally downregulated protein 4-2 (Nedd4-2), a homologous to the E6-AP carboxyl terminus (HECT) ubiquitin ligase, triggers the endocytosis and degradation of its downstream target molecules by regulating signal transduction through interactions with other targets, including 14-3-3 proteins. In our previous study, we found that 14-3-3 binding induces a structural rearrangement of Nedd4-2 by inhibiting interactions between its structured domains. Here, we used time-resolved fluorescence intensity and anisotropy decay measurements, together with fluorescence quenching and mass spectrometry, to further characterize interactions between Nedd4-2 and 14-3-3 proteins. The results showed that 14-3-3 binding affects the emission properties of AEDANS-labeled WW3, WW4, and, to a lesser extent, WW2 domains, and reduces their mobility, but not those of the WW1 domain, which remains mobile. In contrast, 14-3-3 binding has the opposite effect on the active site of the HECT domain, which is more solvent exposed and mobile in the complexed form than in the apo form of Nedd4-2. Overall, our results suggest that steric hindrance of the WW3 and WW4 domains combined with conformational changes in the catalytic domain may account for the 14-3-3 binding-mediated regulation of Nedd4-2.

Proton TOCSY NMR relaxation rates quantitate protein side chain mobility in the Pin1 WW domain.

Protein side chain dynamics play a vital role in many biological processes, but differentiating mobile from rigid side chains remains a technical challenge in structural biology. Solution NMR spectroscopy is ideally suited for this but suffers from limited signal-to-noise, signal overlap, and a need for fractional 13C or 2H labeling. Here we introduce a simple strategy measuring initial 1H relaxation rates during a 1H TOCSY sequence like DIPSI-2, which can be appended to the beginning of any multi-dimensional NMR sequence that begins on 1H. The TOCSY RF field compels all 1H atoms to behave similarly under the influence of strong coupling and rotating frame cross-relaxation, so that differences in relaxation rates are due primarily to side chain mobility. We apply the scheme to a thermostable mutant Pin1 WW domain and demonstrate that the observed 1H relaxation rates correlate well with two independent NMR measures of side-chain dynamics, cross-correlated 13C relaxation rates in 13CßH2 methylene groups and maximum observable 3J couplings sensitive to the χ1 side chain dihedral angle (3JHα,Hß, 3JN,Hß, and 3JCO,Hß). The most restricted side chains belong to Trp26 and Asn40, which are closely packed to constitute the folding center of the WW domain. None of the other conserved aromatic residues is as immobile as the first tryptophan side chain of the WW domain. The proposed 1H relaxation methodology should make it relatively easy to measure side chain dynamics on uniformly 15N- or 13C-labeled proteins, so long as chemical shift assignments are obtainable.

Target Recognition in Tandem WW Domains: Complex Structures for Parallel and Antiparallel Ligand Orientation in h-FBP21 Tandem WW.

Protein-protein interactions often rely on specialized recognition domains, such as WW domains, which bind to specific proline-rich sequences. The specificity of these protein-protein interactions can be increased by tandem repeats, i.e., two WW domains connected by a linker. With a flexible linker, the WW domains can move freely with respect to each other. Additionally, the tandem WW domains can bind in two different orientations to their target sequences. This makes the elucidation of complex structures of tandem WW domains extremely challenging. Here, we identify and characterize two complex structures of the tandem WW domain of human formin-binding protein 21 and a peptide sequence from its natural binding partner, the core-splicing protein SmB/B'. The two structures differ in the ligand orientation and, consequently, also in the relative orientation of the two WW domains. We analyze and probe the interactions in the complexes by molecular simulations and NMR experiments. The workflow to identify the complex structures uses molecular simulations, density-based clustering, and peptide docking. It is designed to systematically generate possible complex structures for repeats of recognition domains. These structures will help us to understand the synergistic and multivalency effects that generate the astonishing versatility and specificity of protein-protein interactions.

Ultrafast folding kinetics of WW domains reveal how the amino acid sequence determines the speed limit to protein folding.

Protein (un)folding rates depend on the free-energy barrier separating the native and unfolded states and a prefactor term, which sets the timescale for crossing such barrier or folding speed limit. Because extricating these two factors is usually unfeasible, it has been common to assume a constant prefactor and assign all rate variability to the barrier. However, theory and simulations postulate a protein-specific prefactor that contains key mechanistic information. Here, we exploit the special properties of fast-folding proteins to experimentally resolve the folding rate prefactor and investigate how much it varies among structural homologs. We measure the ultrafast (un)folding kinetics of five natural WW domains using nanosecond laser-induced temperature jumps. All five WW domains fold in microseconds, but with a 10-fold difference between fastest and slowest. Interestingly, they all produce biphasic kinetics in which the slower phase corresponds to reequilibration over the small barrier (<3 RT) and the faster phase to the downhill relaxation of the minor population residing at the barrier top [transition state ensemble (TSE)]. The fast rate recapitulates the 10-fold range, demonstrating that the folding speed limit of even the simplest all-ß fold strongly depends on the amino acid sequence. Given this fold's simplicity, the most plausible source for such prefactor differences is the presence of nonnative interactions that stabilize the TSE but need to break up before folding resumes. Our results confirm long-standing theoretical predictions and bring into focus the rate prefactor as an essential element for understanding the mechanisms of folding.

PEGylation Increases the Strength of a Nearby NH-π Hydrogen Bond in the WW Domain.

Here we show that an NH-π interaction between a highly conserved Asn and a nearby Trp stabilizes the WW domain of the human protein Pin1. The strength of this NH-π interaction depends on the structure of the arene, with NH-π interactions involving Trp or naphthylalanine being substantially more stabilizing than those involving Tyr or Phe. Calculations suggest arene size and polarizability are key structural determinants of NH-π interaction strength. Methylation or PEGylation of the Asn side-chain amide nitrogen each strengthens the associated NH-π interaction, though likely for different reasons. We hypothesize that methylation introduces steric clashes that destabilize conformations in which the NH-π interaction is not possible, whereas PEGylation strengthens the NH-π interaction via localized desolvation of the protein surface.

Coupled intra- and interdomain dynamics support domain cross-talk in Pin1.

The functional mechanisms of multidomain proteins often exploit interdomain interactions, or "cross-talk." An example is human Pin1, an essential mitotic regulator consisting of a Trp-Trp (WW) domain flexibly tethered to a peptidyl-prolyl isomerase (PPIase) domain, resulting in interdomain interactions important for Pin1 function. Substrate binding to the WW domain alters its transient contacts with the PPIase domain via means that are only partially understood. Accordingly, we have investigated Pin1 interdomain interactions using NMR paramagnetic relaxation enhancement (PRE) and molecular dynamics (MD) simulations. The PREs show that apo-Pin1 samples interdomain contacts beyond the range suggested by previous structural studies. They further show that substrate binding to the WW domain simultaneously alters interdomain separation and the internal conformation of the WW domain. A 4.5-µs all-atom MD simulation of apo-Pin1 suggests that the fluctuations of interdomain distances are correlated with fluctuations of WW domain interresidue contacts involved in substrate binding. Thus, the interdomain/WW domain conformations sampled by apo-Pin1 may already include a range of conformations appropriate for binding Pin1's numerous substrates. The proposed coupling between intra-/interdomain conformational fluctuations is a consequence of the dynamic modular architecture of Pin1. Such modular architecture is common among cell-cycle proteins; thus, the WW-PPIase domain cross-talk mechanisms of Pin1 may be relevant for their mechanisms as well.

Viruses go modular.

The WW domain is a modular protein structure that recognizes the proline-rich Pro-Pro-x-Tyr (PPxY) motif contained in specific target proteins. The compact modular nature of the WW domain makes it ideal for mediating interactions between proteins in complex networks and signaling pathways of the cell (e.g. the Hippo pathway). As a result, WW domains play key roles in a plethora of both normal and disease processes. Intriguingly, RNA and DNA viruses have evolved strategies to hijack cellular WW domain-containing proteins and thereby exploit the modular functions of these host proteins for various steps of the virus life cycle, including entry, replication, and egress. In this review, we summarize key findings in this rapidly expanding field, in which new virus-host interactions continue to be identified. Further unraveling of the molecular aspects of these crucial virus-host interactions will continue to enhance our fundamental understanding of the biology and pathogenesis of these viruses. We anticipate that additional insights into these interactions will help support strategies to develop a new class of small-molecule inhibitors of viral PPxY-host WW-domain interactions that could be used as antiviral therapeutics.

Requirement for expression of WW domain containing transcription regulator 1 in bovine trophectoderm development.

WW domain-containing transcription regulator 1 (WWTR1) is one of the primary effectors in the Hippo pathway, which plays essential roles in cell differentiation into trophectoderm (TE) and inner cell mass cell lineages at the blastocyst stage. However, little is known about the roles of WWTR1 in preimplantation development. The present study aimed to explore the significance of WWTR1 expression in preimplantation development using an mRNA knockdown (KD) system in bovine embryos. We first quantitated WWTR1 expression at protein and mRNA levels from fertilization to blastocyst stage. WWTR1 proteins gradually shifted from extranuclear localization during the 16-cell stage to nuclear localization by morula stage. WWTR1 mRNA expression was also transiently upregulated at the 16-cell stage. WWTR1 KD efficiently repressed WWTR1 expression at protein and mRNA levels. The WWTR1 KD embryos developed to the blastocyst stage at rates equivalent to those of controls, but TE cell numbers were significantly decreased. Representative TE-expressed genes, including CDX2 and IFNT were also significantly decreased in WWTR1 KD blastocysts. These results provide the first demonstration that WWTR1 expression is responsible for normal TE cell development in preimplantation embryos.

Ifu5, a WW domain-containing protein interacts with Efg1 to achieve coordination of normoxic and hypoxic functions to influence pathogenicity traits in Candida albicans.

Hypoxic adaptation pathways, essential for Candida albicans pathogenesis, are tied to its transition from a commensal to a pathogen. Herein, we identify a WW domain-containing protein, Ifu5, as a determinant of hypoxic adaptation that also impacts normoxic responses in this fungus. Ifu5 activity supports glycosylation homeostasis via the Cek1 mitogen-activated protein kinase-dependent up-regulation of PMT1, under normoxia. Transcriptome analysis of ifu5Δ/Δ under normoxia shows a significant up-regulation of the hypoxic regulator EFG1 and EFG1-dependent genes. We demonstrate physical interaction between Ifu5 by virtue of its WW domain and Efg1 that represses EFG1 expression under normoxia. This interaction is lost under hypoxic growth conditions, relieving EFG1 repression. Hypoxic adaptation processes such as filamentation and biofilm formation are affected in ifu5Δ/Δ cells revealing the role of Ifu5 in hypoxic signalling and modulating pathogenicity traits of C. albicans under varied oxygen conditions. Additionally, the WW domain of Ifu5 facilitates its role in hypoxic adaptation, revealing the importance of this domain in providing a platform to integrate various cellular processes. These data forge a relationship between Efg1 and Ifu5 that fosters the role of Ifu5 in hypoxic adaptation thus illuminating novel strategies to undermine the growth of C. albicans.

Markov state models from hierarchical density-based assignment.

Markov state models (MSMs) have become one of the preferred methods for the analysis and interpretation of molecular dynamics (MD) simulations of conformational transitions in biopolymers. While there is great variation in terms of implementation, a well-defined workflow involving multiple steps is often adopted. Typically, molecular coordinates are first subjected to dimensionality reduction and then clustered into small "microstates," which are subsequently lumped into "macrostates" using the information from the slowest eigenmodes. However, the microstate dynamics is often non-Markovian, and long lag times are required to converge the relevant slow dynamics in the MSM. Here, we propose a variation on this typical workflow, taking advantage of hierarchical density-based clustering. When applied to simulation data, this type of clustering separates high population regions of conformational space from others that are rarely visited. In this way, density-based clustering naturally implements assignment of the data based on transitions between metastable states, resulting in a core-set MSM. As a result, the state definition becomes more consistent with the assumption of Markovianity, and the timescales of the slow dynamics of the system are recovered more effectively. We present results of this simplified workflow for a model potential and MD simulations of the alanine dipeptide and the FiP35 WW domain.

Intrinsic disorder and amino acid specificity modulate binding of the WW2 domain in kidney and brain protein (KIBRA) to synaptopodin.

The second WW domain (WW2) of the kidney and brain scaffolding protein, KIBRA, has an isoleucine (Ile-81) rather than a second conserved tryptophan and is primarily unstructured. However, it adopts the canonical triple-stranded antiparallel ß-sheet structure of WW domains when bound to a two-PPXY motif peptide of the synaptic protein Dendrin. Here, using a series of biophysical experiments, we demonstrate that the WW2 domain remains largely disordered when bound to a 69-residue two-PPXY motif polypeptide of the synaptic and podocyte protein synaptopodin (SYNPO). Isothermal titration calorimetry and CD experiments revealed that the interactions of the disordered WW2 domain with SYNPO are significantly weaker than SYNPO's interactions with the well-folded WW1 domain and that an I81W substitution in the WW2 domain neither enhances binding affinity nor induces substantial WW2 domain folding. In the tandem polypeptide, the two WW domains synergized, enhancing the overall binding affinity with the I81W variant tandem polypeptide 2-fold compared with the WT polypeptide. Solution NMR results showed that SYNPO binding induces small but definite chemical shift perturbations in the WW2 domain, confirming the disordered state of the WW2 domain in this complex. These analyses also disclosed that SYNPO binds the tandem WW domain polypeptide in an antiparallel manner, that is, the WW1 domain binds the second PPXY motif of SYNPO. We propose a binding model consisting of a bipartite interaction mode in which the largely disordered WW2 forms a "fuzzy" complex with SYNPO. This binding mode may be important for specific cellular functions.

Growth arrest specific protein 7 inhibits tau fibrillogenesis.

Here we show that Gas7 inhibits phosphorylated tau fibrillogenesis by binding to phosphorylated tau at its non-WW domain, presumably F-BAR domain. We revealed that Gas7 binds to the third repeat domain of tau, the core element of tau oligomerization and the C-terminal domain of tau and alters the conformation not to form fibrils. These results suggest that Gas7 may serve to protect against Alzheimer's disease and other tauopathies by preventing tau fibrillogenesis.

PEGylation near a Patch of Nonpolar Surface Residues Increases the Conformational Stability of the WW Domain.

Many proteins have one or more surface-exposed patches of nonpolar residues; our observations here suggest that PEGylation near such locations might be a useful strategy for increasing protein conformational stability. Specifically, we show that conjugating a PEG-azide to a propargyloxyphenylalanine via the copper(I)-catalyzed azide-alkyne cycloaddition can increase the conformational stability of the WW domain due to a favorable synergistic effect that depends on the hydrophobicity of a nearby patch of nonpolar surface residues.

Evaluation of the relationships of the WBP1L gene with schizophrenia and the general psychopathology scale based on a case-control study.

WBP1L is a target of microRNA 137 (miR-137) and has been considered a candidate gene for schizophrenia (SCZ). To investigate the relationships between WBP1L and SCZ and its related symptom scales, a total of 5,993 Chinese Han subjects, including 2,128 SCZ patients and 3,865 controls, were enrolled. In addition, an independent sample set for replication study including 1,052 SCZ patients and 2,124 controls were also recruited. Thirty-two tag single nucleotide polymorphisms (SNPs) located within gene region of WBP1L were selected for genotyping and analyzing. The expression quantitative trait loci (eQTL) effects for the targeted SNPs were investigated with gene expression data from multiple human tissues. Rs4147157 (OR = 0.84, p = 1.51 × 10-5 ) and rs284854 (OR = 1.14, p = 7.00 × 10-4 ) were significantly associated with SCZ disease status and these association signals were replicated in our replication sample. A significant association was identified between rs4147157 and the general (ß = -.66, p = .001) and total (ß = -.8, p = .0042) scores of positive and negative syndrome scale scores in SCZ patients. Both SNPs were significant eQTL for genes around WBP1L in human brain tissues including ARL3 and AS3MT. To conclude, SNPs rs4147157 and rs284854 were associated with SCZ in the Chinese Han population. Additionally, rs4147157 was significantly associated with specific symptom features of SCZ.

Substrate clustering potently regulates the activity of WW-HECT domain-containing ubiquitin ligases.

The Nedd4 family of HECT domain-containing E3 ligases ubiquitinate many transcription factors and signaling proteins, and their activity is tightly regulated. Normally, intramolecular interactions curb the catalytic activity of the HECT domain, but these can be broken by the binding of PY motifs, found on substrate molecules and adaptors, to the WW domains characteristic of this E3 ligase family. This raises the prospect of substrates automatically activating the ligases, frustrating the purpose of ligase regulation. Here we show that soluble protein substrates and adaptors such as α arrestins, even with multiple PY elements, cannot activate ligase activity efficiently. However, we found that polymerization or membrane tethering of these substrates dramatically increases the ligase activity both in vivo and in vitro Aggregation of luciferase-containing substrates upon heat shock had a similar effect and could also expose cryptic PY elements in the substrates. We inferred that ligase activation critically requires a substantial array of clustered PY motifs and that the formation of such arrays on membranes or in polymeric aggregates may be an essential step in this mode of ligase regulation. We conclude that recruitment of α arrestins to membrane receptors and aggregation of unstable proteins after heat shock may be physiologically relevant mechanisms for triggering ubiquitination by Nedd4 family HECT domain-containing E3 ligases.