If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Department of Microbiology and Molecular Genetics, Faculty of Medicine, Institute of Medical Research Israel-Canada, The Hebrew University of Jerusalem, Jerusalem, Israel
Department of Microbiology and Molecular Genetics, Faculty of Medicine, Institute of Medical Research Israel-Canada, The Hebrew University of Jerusalem, Jerusalem, Israel
Department of Microbiology and Molecular Genetics, Faculty of Medicine, Institute of Medical Research Israel-Canada, The Hebrew University of Jerusalem, Jerusalem, Israel
Department of Microbiology and Molecular Genetics, Faculty of Medicine, Institute of Medical Research Israel-Canada, The Hebrew University of Jerusalem, Jerusalem, Israel
Wolfson Centre for Applied Structural Biology, Hebrew University of Jerusalem, Jerusalem, IsraelDepartment of Pharmaceutical Engineering, Azrieli College of Engineering, Jerusalem, Israel
Department of Microbiology and Molecular Genetics, Faculty of Medicine, Institute of Medical Research Israel-Canada, The Hebrew University of Jerusalem, Jerusalem, Israel
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.
WW domains are small 38 to 40 residue modules that adopt a three-stranded antiparallel β-sheet fold. They are named for their two conserved tryptophan (W) residues. The first W helps stabilize the domain, whereas the second W facilitates binding to short linear motifs rich in proline (P) (
). WW domains mediate protein–protein interactions involved in a range of protein functions, from degradation by ubiquitination to nuclear transport. These protein functions are instrumental in determining cell fate by mediating apoptosis and cell proliferation.
WW domains usually occur in tandem, allowing for fine-tuned regulation through a combination of binding events (
). This complexity is enhanced by the presence of multiple proline-rich binding motifs, such as the class I PPxY motif, on partner proteins. Many of the reported interactions between WW domain proteins and PPxY motif partners involve more than one WW domain–PPxY motif pair. For example, the tandem WW domain of Yki, the Drosophila melanogaster homolog of human Yes-associated protein (YAP), simultaneously binds two Tgi (Tondu domain–containing growth inhibitor) PPxY peptide motifs (
). Moreover, the WW tandem domains of proteins, such as formin-binding protein 21 (FBP21), YAP, and TAZ (PDZ-binding motif), bind dual PPxY peptide motifs with higher affinity than a corresponding single PY peptide motif (
). It was also shown that while a minimal stretch that contains the three PPxY motifs of NDFIP2 retains the ability to activate ITCH, mutation of any one of the PPxY elements reduces activity (
Human WW domain–containing oxidoreductase (WWOX) is a tumor suppressor involved in many biological functions, such as apoptosis, DNA damage, inhibition of cell growth, cellular metabolism, and proper neurodevelopment (
). WWOX contains two N-terminal WW domains (WW1 and WW2) separated by a short linker as well as a C-terminal short-chain dehydrogenase/reductase (SDR) domain (Fig. 1) (
). While the WWOX domain architecture is conserved, the two WWOX homologs found in Caenorhabditis elegans contain the SDR domains with a similar substrate-binding site but lack the WW domains altogether, suggesting independent SDR and WW activity (
). In contrast to the SDR domain, much more is known about the WW domain activity. WWOX binds many class I proline-rich motif-containing proteins, including ErbB4 and p73 (
), and competition with other WW domain proteins, such as YAP, for these binding partners has been shown to influence cellular behavior. Thus, WWOX sequesters ErbB4 in the cytosol by preventing its binding to YAP and its subsequent shuttling to the nucleus for the initiation of a transcriptional program leading to cellular proliferation (
). While it would seem that the two WWOX domains are a prototypical case of the dual interactions described previously, the second of these, WW2, is in fact considered incapable of independently binding PPxY peptides with measurable affinity, because of the W85Y mutation in the conserved binding pocket (In fact, the reverse Y85W mutation re-establishes binding for this domain (
)). WWOX therefore seems to bind PPxY partners exclusively through WW1. What then is the role of WW2, and how is WWOX able to compete with WW proteins with two fully functioning WW domains, especially for partner proteins that contain multiple PPxY motifs? It has been suggested that WW2 functions as a chaperone that stabilizes WW1 to improve the latter’s binding affinity to its partners (
), but the details of such a possible mechanism remain to be elucidated.
Figure 1WWOX and its substrate ErbB4. WWOX contains two WW domains (residues 16–50 and 57–91, respectively) and an SDR domain (residues 122–383; domain position annotations according to InterPro (
)). Isolated WW1 is unstable (dashed lines) but binds to PPxY motifs (filled crescent), whereas WW2 is stable (solid lines) but was not reported to bind to PPxY motifs with measurable affinity (dashed crescent (
)). WWOX partner ErbB4 is cleaved upon maturation, releasing the C-terminal cytoplasmic region (residues 676–1308) that contains a kinase domain (positions 718–974) followed by an extended disordered region (gray curve) that contains several recognition motifs, including the three WWOX-binding PPxY motifs investigated here (PY1 1032–1035, PY2 1053–1056, and PY3 1298–1301; black ovals). SDR, short-chain dehydrogenase/reductase; WWOX, WW domain–containing oxidoreductase.
Here, we have studied the structural and functional role of WWOX WW2 in binding WWOX partners, focusing on ErbB4 and its three PPxY motifs in the intracellular domain that is released to the cytoplasm by γ-secretase after receptor stimulation (
) (Fig. 1). Using CD and NMR experiments, we found that the presence of the WW2 domain, as well as the binding of the peptide, induces WW1 structural stabilization. Furthermore, our isothermal titration calorimetry (ITC) and NMR protein–peptide binding assays show increased affinity of a double-PPxY peptide to the WWOX tandem WW domains when compared with a single PPxY peptide binding to the tandem domain or to isolated WW1 binding a double-PPxY. This suggests that WW2 can engage a suitably oriented PPxY motif, despite its missing key tryptophan residue. Also, NMR spectra establish that the double PPxY peptide induces a significantly more structured conformation of WWOX than the single PPxY peptide, once again supporting direct involvement of WW2. By comparing affinities and chemical shift perturbation (CSP) effects induced by native, engineered, and mutated double-PPxY motif peptides, we deduced the binding pose and directionality of two-site ligand binding and the contribution of different molecular determinants to affinity. In aggregate, these data suggest a plausible model for the relative orientation between the two WWOX WW domains and reveal details of the interaction between WWOX WW2 and WW1 domains and substrate single- and double-motif peptides, which are discussed in the functional context of WWOX–substrate interactions.
Results
To elucidate the role played by WW2 in WWOX functionality, we investigated its influence on the stability and binding of WWOX. We compared single WW1 and WW2 domains to the tandem domain WW1–WW2 in terms of structural stability and binding affinity to different peptides derived from ErbB4, a known WWOX substrate (
We used both denaturation experiments and CD spectra to characterize the stability of the isolated and tandem WWOX domains. Monitoring tryptophan fluorescence changes in urea denaturation experiments showed that WW1 is structured only within the context of the WW1–WW2 tandem domain, but not in its isolated form, indicating significant stabilization of WW1 by WW2 (Fig. 2A). Folded WW domains exhibit a CD spectrum dominated by a strong positive peak at 225 to 230 nm contributed by ordered aromatic side chains and a negative peak at >205 nm (
Characterization of the structure and function of W--> F WW domain variants: identification of a natively unfolded protein that folds upon ligand binding.
). The CD spectrum of isolated WW2 indeed conforms to this characteristic spectrum, albeit with a relatively weak positive peak and a negative peak at <205 nm, expected because of the W85Y mutation, whereas isolated WW1 is predominantly unstructured (Fig. 2B). In comparison, two important spectral changes appear in the spectrum of tandem WW1–WW2. The WW1–WW2 positive peak is synergistically stronger than the sum of the WW1 and WW2 contributions, indicating mutual influence of the domains on each other resulting in a structure closer to a typical WW domain (
Characterization of the structure and function of W--> F WW domain variants: identification of a natively unfolded protein that folds upon ligand binding.
). In addition, the negative band in the CD spectrum shifts to the right from 205 nm, indicating that the WW domain adopts a more folded antiparallel β-sheet structure. Together, these changes reinforce the view that WW2 significantly stabilizes WW1, as described before (
Figure 2WWOX domain 1 (WW1) is stabilized by domain 2 (WW2).A, tryptophan fluorescence assay for protein in native conditions (full) and in 6 M urea (empty). Whereas isolated WW1 (circles) exhibits little change, both WW2 and tandem WW1–WW2 domains (triangles and squares, respectively) exhibit significant change indicating a loss of structure. The larger change in the tandem domain as compared with WW2 indicates stabilization of WW1 by WW2. B, CD curves for isolated WW1 (circles), isolated WW2 (triangles), and the tandem WW1–WW2 domains (line). Whereas WW1 is mostly unfolded, the latter two are ordered, showing a characteristic positive peak for WW domains around 225 to 230 nm. The tandem domain also shows a shift of the negative peak away from 205 nm, indicative of an extended β-sheet structure, and is significantly more structured than what would be observed from combining the spectra of the two individual domains (plus signs), highlighting the stabilization effect of WW2 on WW1. WWOX, WW domain–containing oxidoreductase.
A molecular description of the WW2–WW1 interaction
For a more detailed picture of this stabilization event, we investigated its accompanying structural changes using well-established NMR-based methods. As 1H,15N-heteronuclear single quantum coherence (HSQC) spectra of proteins are highly sensitive to the local electronic environment, shifted crosspeaks in this “fingerprint” spectrum (known as CSPs) report on local and global changes in structure and dynamics on a per-residue basis. The WW1–WW2 HSQC showed a striking presence of two subpopulations of crosspeaks, one exhibiting line widths consistent with a small protein (below 10 kDa) and the other suffering from extensive line broadening that was aggravated at higher temperatures (data not shown). By recording spectra for the single domains, we confirmed that the broadened peaks belong to WW1 residues, whereas WW2 affords a spectrum with narrow peaks (Fig. 3, A and B). This broadening and its temperature dependence suggested that WW1 loses signal intensity because of solvent-exchange, reflecting a less-folded and more flexible domain when compared with WW2.
Figure 3NMR HSQC spectra of single and tandem WWOX WW domains reveal an interdomain-stabilizing interaction.A, the 1H,15N-HSQC spectrum of WW1 (red) shows a typical spectral dispersion of an unfolded domain, which changes to defined peaks in the WWOX tandem domain (blue), indicating that WW2 stabilizes WW1 (see Fig. S1 for peak assignments). B, in contrast, structured WW2 (red) gives a well-behaved spectrum with many overlapping peaks with the WWOX tandem domain (blue). Peaks discussed in the text are annotated and highlighted by circles. C, chemical shift perturbations (CSPs) of WW2 residues in the presence of WW1 (tandem domain). Significant CSPs (Δδ > 0.05 ppm) are color coded for their main regions, including the β2/β3 hairpin (green) and the C-terminal region (cyan); peaks that disappear/appear upon addition of WW1 (not fully assigned but clearly affected) are shown by blue (vertically/horizontally) hatched bars; stable peaks in the β1/β2 hairpin region are colored yellow. D, mapping of changes onto the structure of the WW2 domain (Protein Data Bank accession code: 1WMV) using the same color scheme. Note the close vicinity of H78 (in blue) and G62 at the basis of the beta hairpin formed by strands β1 and β2 and other residues. Right, detailed structure of the β2/β3 hairpin. Figures of structures were generated using PyMOL (
An analysis of differences between the two single-domain spectra and that of the WW1–WW2 tandem domain was instructive. WW1 crosspeaks in the tandem domain exhibited a significant global change, specifically an appearance of well-dispersed peaks and a reduction of peaks in the region characteristic of unstructured proteins, suggesting a stabilizing effect of WW2 upon the neighboring domain (Fig. 3A). At the same time, WW2 peaks in two regions exhibited small yet significant intensity changes or cross-peak shifts (Fig. 3, B and C; see the Experimental procedures section for choice of the 0.05 ppm significance threshold). The first includes a cluster of residues on one side of the WW2 β-sheet, involving G62 (preceding β1) and the β2β3 turn (particularly H78, N80, and R82). The second region consists of C-terminal residues R89–L90–A91, whose chemical shifts suggest that this segment packs against the WW2 domain core, rather than adopting a free random-coil conformation, and that this interaction with the domain core is modified upon linking WW2 to WW1. These two effects could be attributed to direct interaction with WW1, or alternatively, to indirect effects propagated through the WW2 β-sheet. In summary, we find that the WW2 “tip,” comprising pre-β1 and β2β3-turn residues, plays an important role, or is at least significantly perturbed upon WW1 stabilization, and that the two WW-domain cores, linker, and C terminus form a contiguous set of interactions that further contribute to the WW1–WW2 structure.
Binding of ErbB4 ligands to WWOX
How does the tandem domain bind its peptide substrates, and what does WW2 contribute to this affinity? Our previous study highlighted the variety of mutual influences between individual WW domains in tandem repeats upon substrate-binding affinity and specificity (
). WWOX stood out because of the inability of its WW2 to bind to canonical substrates (e.g., PPxY motifs) at measurable affinities, because of a tryptophan-to-tyrosine mutation at the second characteristic W-position. To study the details of the suggested chaperone role of WW2 in WWOX binding of PPxY motifs, we compared the affinities of isolated WW1 domain and tandem WW1–WW2 to various ErbB4-derived PPxY-containing peptides (Fig. 1 and Table 1). Affinities were measured using both ITC and NMR, a combination that allowed us to cover a wide range of affinities including often-inaccessible KD values in the 200 to 1000 μM range.
In terms of affinity, ITC results differentiated between the affinity of tandem WW1WW2 to PY3 (KD = 30 μM) and PY1 or PY2 (estimated weak binding, KD > 100 μM) (Table 2), in qualitative agreement with previous studies (
). However, although no significant affinity to the isolated WW2 domain was detected for either peptide (KD >> 200 μM), ITC showed an approximate twofold decrease in PY3-binding affinity when the WW2 domain was removed (KD of 78 versus 30 μM) (Fig. 4, A and B). The 15N,1H–HSQC spectrum shows that binding PY3 induces a stabilization of WW1 in the tandem domain, as demonstrated by a twofold increase in visible crosspeaks, several of which representing β-sheet residues (see Fig. 5 discussed later and Figs. S2 and S5A). In terms of stability, changes in the CD spectrum of WW1WW2 upon addition of PY3 reveal that the unusually unstable WW1 is subject to further ligand-induced stabilization (Fig. 4D), in addition to the stabilization by the WW2 in the tandem domain (Fig. 3A).
Table 2Binding affinities of PPxY peptides to WWOX single and tandem WW domains as determined by ITC (see Experimental procedures section)
Figure 4Substrate-binding affinity is increased in the presence of WW2 and for a double-PPxY peptide.A–C, isothermal titration calorimetry (ITC) curves of (A) WW1, (B) WW1–WW2 binding to single-motif peptide PY3, and (C) of WW1–WW2 binding to double-motif peptide PY3PY3. Calculated affinities are shown. D, CD curves of WWOX WW1–WW2 with no peptide (solid line), when bound to single-motif peptide PY3 (dots), or double-motif peptide PY3PY3 (dashed line). The higher intensity at 225 nm is due to increased order in aromatic side chains and to a polyproline type II (PPII) helical conformation of the peptide induced by the binding. Peptide CD curves were subtracted to highlight WW contributions to the spectra (For curves before subtraction and curves of isolated peptides, see Fig. S3).
Figure 5Binding of dual-PPxY peptides to WWOX.A, details of 15N,1H-HSQC spectra of the unbound tandem domains (gray) and after addition of single and dual-PPxY peptides at 2:1 peptide:WWOX ratio (PY3 [black] and PY3PY3 [purple], respectively). B, specific peaks with distinct perturbation patterns upon addition of various peptides, including PY3, PY3PY3, and also PY1PY3 (blue), PY1PY2 (green), PY3AY3 (pink), and AY3PY3 (red) (see Table 1 for peptide sequences). Peaks that disappear are shown in dotted lines. C, plots summarizing the CSPs observed after addition of 2 M equivalents of the various peptides. CSP of 0.3 designates a peak that was broadened beyond detection. D, corresponding plots for mutant peptides. CSP, chemical shift perturbation; HSQC, heteronuclear single quantum coherence; WWOX, WW domain–containing oxidoreductase.
Assuming that binding of PY3 involves predominantly WW1 (as suggested from Fig. 5C, first panel, shown later), analysis of enthalpy and entropy contributions to binding allows us to quantify the stabilization of WW1 by WW2: The reduced binding enthalpy when adding PY3 to the tandem domain compared with the single domain (as well as the corresponding decrease in loss of conformational entropy) reflects the ordering of WW1 prior to binding. Thus, half of the binding enthalpy (or more) gained in binding of WW1 to PY3 (−11 of −23 kcal/mol) is contributed by WW2-mediated stabilization. Of note, ΔGbinding does not change significantly; both constructs lead to similar overall binding affinity—only 2.5× increase (from a KD of 78 μM to 30 μM). WW2 stabilizes WW1 in a peptide binding–like conformation, as suggested by the similarity of the HSQC spectra of PY3 bound to isolated WW1, and WW1 connected to WW2 in the tandem domain context (Fig. S2).
A significant increase in binding affinity is observed for double binding motifs
While WW domain crosstalk has been extensively discussed (
), less is known about how binding motif multiplicity in the partner affects affinity and specificity. To unravel the effect of motif repeats in the substrate on WW domain binding, we generated several double-motif peptides and measured their binding affinity to WWOX (Tables 1 and 2). In these double-motif ligands, the two PPxY motifs are connected by a natural linker segment in the case of the proximal PY1PY2, and a designed polyglycine linker for the non-native PY1PY3 and PY3PY3 peptides (as well as the PY1PY2S peptide, in which the natural linker was replaced by the same polyglycine linker). Such duplication of motifs significantly increased binding affinity, often by an order of magnitude, even when the isolated motifs (i.e., PY1 or PY2) failed to bind. Thus, PY1PY2 and PY1PY2S bind WW1WW2 with a KD of 14 and 21 μM, respectively (as opposed to >100 μM for PY1 and PY2 alone), PY1PY3 binds with a KD of 8 μM (stronger than PY1 and PY3 alone), and PY3PY3 binds with a KD of 3 μM (Fig. 4C, 30 μM for PY3 alone). The PY3PY3 affinity was the strongest measured in this study. This increase in affinity for the dual motif peptides could be seen in CD experiments as well, as these suggest that peptide-induced tightening of the WW domain was more pronounced for double-motif ligands (Fig. 4D). Although an increase in β-sheet content induced by the peptides could not be sufficiently separated from absorbance of the free peptide molecules (Fig. S3), these results do suggest that tandem PPxY motif binding induces further structural changes beyond those induced by a single PPxY ligand.
To quantify the contribution of an avidity effect brought about by the two proximal ligands, we measured the binding of PY3PY3 to the single WW1 domain and observed only a small change in affinity (less than twofold). Moreover, a size-exclusion chromatography–multiangle light scattering (SEC–MALS) WW1WW2 elution peak had a calculated molecular weight of ca. 10 kDa (Fig. S4), ruling out the possibility that it behaves as a dimer (as occurs, e.g., in the SAV1 protein (
)). These results prove that WW2 plays an active role in peptide binding and that increased binding affinity is only partially due to avidity.
Deconvoluting domain contributions to double binding motifs
NMR offers a dual advantage in characterizing the WW–PPxY interactions. First, it is capable of identifying weak affinities (>100 μM) by following chemical shifts under conditions of fast exchange on the NMR timescale. Second, and more importantly, since NMR monitors local changes at each residue, it provides domain-specific binding information, as opposed to the global view afforded by ITC. This was clearly demonstrated by the ITC-derived affinities for PY1PY2S with the polyglycine linker (21 μM), compared with PY1PY2 with the natural linker (14 μM); there a convolution of two binding contributions actually resulted in a seemingly lower affinity (for PY1PY2S) when compared with single binding site fitting (for PY1PY2). We exploited this feature of the WWOX NMR fingerprint spectra to analyze the relative contributions of WW1 and WW2 to the binding of double-PPxY peptides. Figure 5, A and B shows the effects of peptides binding to WW1–WW2, and Figure 5C summarizes CSPs in the WW1WW2 spectrum for the single-motif PY3 and the four double-motif peptides PY1PY2 (native linker), PY1PY3, PY3PY3, and PY1PY2S (with polyglycine linker), reporting on contacts between WWOX and the various ligands. The HSQC spectrum (Fig. 5A) shows that single-motif PY3 induced significant peak shifts for most assigned WW1 peaks, in particular for beta-strand residues and the W44 indole signal (Fig. 5B), in agreement with the canonical WW–PPxY binding mode. Concomitantly, only a few WW2 peaks, located in the β3-strand, changed significantly upon binding of PY3. Thus, increase in local concentration of PPxY ligand in the PY3PY3 peptide allows the interaction of the W-deficient WW2 domain to the PPxY motif. However, changes induced by the three double-motif peptides draw a more complex picture (Fig. 5C). The PY3 CSP pattern in WW1 is echoed in the spectra of double-motif peptides, suggesting they all employ a similar binding mode. In contrast, the WW2 CSP pattern distinguished between native sequence PY1PY2, reminiscent of PY3, and PY1PY2S (poly G linker), PY1PY3, and PY3PY3 for which large WW2 CSPs suggest a strong interaction with this domain as well. To some extent, this is due to the longer native linker connecting PY1 and PY2, since the interaction grew stronger when a shorter flexible non-native polyglycine linker (similar to the one in the other two double-motif peptides) was used (Fig. S5B). Overall, the size of the average WW2 CSP was consistent with the affinities of the double motifs with the shorter polyglycine linker, PY1PY2S < PY1PY3 < PY3PY3, suggesting a common WW1 interaction mode and a difference in contribution of WW2 to the affinity. For a more quantitative view of domain-specific contributions to affinity, we assumed a two-site binding event with distinct WW1- and WW2-related affinities and followed the concentration-dependent CSP along a titration curve for WW1 and WW2 crosspeaks separately (e.g., the convenient W44 indole and V73 backbone crosspeaks, easily identified in all spectra; Fig. S6). This resulted in domain-specific values for apparent affinities (Table 3 and Experimental procedures section). PY1PY2 exhibits apparent affinity of ∼1000 μM to WW2, whereas PY1PY3, PY3PY3, and PY1PY2S exhibit appreciable apparent WW2 affinities of 80, 40, and 70 μM, respectively. Altogether, this demonstrates the role of the linker, beyond the affinities of the individual peptide motifs.
Table 3Binding affinities of single and dual PPxY peptides to WWOX tandem domain as determined by NMR (see Experimental procedures section)
WWOX
Ligand
KD WW1 (μM)
KD WW2 (μM)
Based on W44
Based on V73
WW1WW2
PY3
10 ± 2
>1000
WW1WW2
PY1
21 ± 2
ND
WW1WW2
PY2
15 ± 4
ND
WW1WW2
PY3PY3
0.6 ± 0.3
40 ± 20
WW1WW2
PY1PY3
6 ± 1
80 ± 20
WW1WW2
PY1PY2
10 ± 2
∼1000
WW1WW2
PY1PY2S
2 ± 1
70 ± 20
See Table 1 for details on WW domains and peptide motifs used.
Interestingly, WW2 CSPs induced by binding to the double-motif peptides are correlated with additional residues, namely WW1 β1 residue T27 and C-terminal residues R89–L90–A91 (Fig. 5B), all located outside the peptide-binding interface. A structural model of the WWOX tandem domain (see later) colocalizes these four residues at the interdomain core. It is thus reasonable that these CSPs reflect the bridging of the two domains by the dual motif peptides, a clear indication that two-site binding is occurring. PY1PY3, PY3PY3, and PY1PY2S form two contact surfaces, whereas PY1PY2 binds mostly with a single motif to WW1. Since PY1PY2 does exhibit significant affinity to WW1WW2, despite this difference, a more complex mode of binding for this peptide may be indicated.
Binding orientation of the double-motif peptide on the tandem domain
To further elucidate the binding mode of the double-motif PY3PY3 peptide to the tandem domain, in particular, the directionality of binding (parallel versus antiparallel), we generated two mutant tandem peptides in which either the N-terminal (termed AY3PY3) or the C-terminal motif (termed PY3AY3) was obviated by PP-to-AA mutations (Tables 1 and 4; mutation of these residues was shown before to abolish binding (
)). Mutation in the first (AY3PY3) or the second (PY3AY3) peptide motif afforded ITC KD values closer to the affinity of the single-motif peptide (PY3, 30 μM) than to that of the double-motif peptide (PY3PY3, 3 μM); AY3PY3 bound slightly stronger (18 μM) than PY3AY3 (27 μM), indicating nonsymmetrical effects of the two mutations on WW1 binding. NMR titration curves afforded poor fits to single-site binding (although KD values in the 10–20 μM range were consistent with ITC results), most likely because of the presence of a residual binding effect of the mutated motif. AY3PY3 did induce slightly larger W44 CSPs (Fig. 5B), again in agreement with the ITC measurements. On the other hand, of the two peptides, PY3AY3 (and not AY3PY3) induced the CSP pattern of interdomain core residues T27/R89–A91 (Figs. 5D and 6, top two spectra) previously established as an indicator of domain bridging and two-site binding. A plausible explanation is found in the WW2 CSP values, particularly for residues T83–T84 that are the largest WW2 change observed in most peptides. While these titrations reflect weak WW2–PY interactions in all cases, in AY3PY3 (but not PY3AY3), the slope of peak displacement increases with increasing peptide concentration, affording an isotherm that cannot be fitted to a single binding event, indicating biphasic behavior (Fig. 6). This can be explained by assuming that single-site WW1–PY3 interactions dominate at low concentrations, whereas at higher concentrations, WW1–AY3 interactions increase, allowing some WW2–PY3 encounters to occur. Conversely, for the PY3AY3 peptide, we observe a normal isotherm corresponding to a single binding event (of very weak affinity), suggesting that the WW1–PY3 interaction is strongest at all concentrations (Fig. 6). Enthalpy change values observed in ITC measurements are in agreement with this hypothesis, as the PY3AY3 and PY3 share indistinguishable ΔH values (12–13 kcal/mol), whereas the AY3PY3 ΔH (15.6 kcal/mol) indicates an additional interaction (Table 4). Using this model to address the question of directionality, we note that PY3AY3, for which the first motif is well anchored to WW1, exhibits an interaction with the interdomain core, whereas AY3PY3, actually the stronger ligand, draws its higher affinity from a combination of interactions but does not appreciably bridge the two domains. For PY3AY3, this supports a mode in which WW1 interacts strongly with the PY3 motif, and WW2 interacts with the AY3 motif (only at high concentrations because of the weakness of the interaction), and in doing so “crosses over” the interdomain core, affording a parallel binding orientation. In contrast to AY3PY3, the alternative PY3AY3 binding pose including a WW2–PY3 interaction (as a single domain) appears to be highly unfavorable.
Table 4Binding affinities of peptides mutated in one of the binding motifs
WWOX
Ligand
KD (μM)
ΔH (kcal/mol)
TΔS (kcal/mol)
WW1WW2
PY3
30 ± 2.3
−12.1 ± 0.5
6.1 ± 0.5
WW1WW2
PY3PY3
3.1 ± 0.3
−24.9 ± 0.9
17.5 ± 0.9
WW1WW2
AY3PY3
18 ± 3.3
−15 ± 1.4
9 ± 1.6
WW1WW2
PY3AY3
27.1∗
−12.9
6.8
KD and ΔH/TΔS are given in micromolar and kilocalorie/mole, respectively. See Table 1 for details on WW domains and peptide motifs used. Values were compiled from n = 3 independent experiments (except for ∗: average of two experiments).
Figure 6Scheme of proposed interaction between WWOX tandem WW domain and double-motif peptides, based on distinct effects of mutation of the first or second PPXY motif on binding affinity (from ITC: values at the top; and NMR titration) and perturbation of WW1 and WW2 domains (from NMR HSQC spectra). Models are based on a parallel binding orientation, see Figure 7. Left, the AY3PY3 peptide with spectral insets showing residues T27, R89, T83, and T84 at 0, 0.5, 1, and 2 mol:mol peptide:WWOX ratios. Right, a similar presentation for the PY3AY3 peptide, PY3AY3 shows no sign of a second binding pose (see text). WW1 and WW2 are shown in light green and orange, respectively, The PY3 (AY3) motifs are depicted in dark green (red). Location of residues T27, R89, L90, and A91 is designated with yellow dots. N and C designated the N and C termini of the double-motif polypeptide. HSQC, heteronuclear single quantum coherence; ITC, isothermal titration calorimetry; WWOX, WW domain–containing oxidoreductase.
One of the most challenging aspects of understanding the interaction between a double WW domain protein and a double binding motif has been providing a structural view that unifies all experimental findings. While most probably the WWOX WW1WW2 tandem domain adopts an ensemble of interchanging conformations, our experiments still suggest that one or a few will strongly dominate this ensemble. Successful crystallization of such assemblies is plagued by relatively flexible domains (as in the case of WW1) and linkers as well as mediocre affinity levels that hamper such efforts. However, in light of the combination of the more global ITC/CD/fluorescence experiments alongside NMR experiments affording local per-residue information, we are now in a position to tackle this structural question using models. In formulating this, we draw upon two important conclusions of the current study: (i) residues T27 (WW1) and R89–L90–A91 (WW2), exhibiting correlated CSPs in all titration experiments, participate in an interdomain core that stabilizes WW1 and (ii) double-motif peptides interact with both WW domains in a parallel orientation, with the N- and C-terminal PPxY motifs binding WW1 and WW2, respectively.
Thanks to significant recent advances in structure prediction, spearheaded by AlphaFold2 (
), we were able to generate a model that agrees well with these constraints (Fig. 7A). This model positions the double-motif peptide PY3PY3 onto the tandem domain in a parallel orientation, in agreement with our experimental results. Mapping of the electrostatic potential on the surface highlights the acidic patches located near the C terminus of the PPxY binding motif (Fig. 7B), explaining the preference in this region for basic residues, as already reported previously (
). The second peptide motif is removed from the WW2 binding site in a corresponding model of PY1PY2S with a synthetic linker, and a model for PY1PY2 connected by a longer natural linker positioned only the first motif into the WW1 binding site, whereas the second motif did not form the canonical interaction between the first proline in the motif and Y85 and moved away.
Figure 7Structural model of the WWOX WW1WW2 tandem domain and its interaction with the tandem PPXY peptide.A, model of WW1WW2 bound to PY3PY3, generated with AlphaFold2 multimer (
). WW1 and WW2 domains are colored in magenta and green, respectively; peptide motifs PPPY are shown in yellow sticks. The parallel orientation adapted by the double-peptide motif when binding to the tandem WW domain pair is apparent. The C-terminal residues of WW2 89RLA91 and WW1 T27 are positioned in the center in this model, explaining the concerted shifts observed for these residues upon binding the double-motif peptide (Fig. 5). In turn, the shifts observed for the beta-sheet edge (namely G62 at the beginning of β1 and H78 and N80 in the β2β3 turn), upon addition of WW1 to WW2, could be explained by a change in the beta sheet, or the strain of the linker that is propagated via G62 to that region. B, electrostatic map of WW1WW2, highlighting the strong negative patches, in particular in the WW2 domain, which is contacted by the positively charged flanking regions RHR in the peptide (in cyan). Positions discussed in the text are labeled. WWOX, WW domain–containing oxidoreductase.
All these features indicate that the tandem domain adopts an open conformation that allows peptide binding, rather than necessitating large rearrangements to free the binding sites. While this suggested model might be the dominant conformation, it is not necessarily the only one (as also indicated by the wide range of different domain–domain orientations that we observe in a larger set of models generated using a number of different approaches, data not shown). The equilibrium between conformations could be shifted to this conformation in particular after binding double-motif peptide PY3PY3, whereas to other conformations upon binding of different peptides, such as PY1PY2.
Discussion
WW2 stabilizes the partially unfolded WW1 domain
As many other small peptide-binding domains, WW domains tend to occur in tandem, which allows for the interaction with several corresponding peptide motifs. The ample information on different tandem motif–containing proteins, such as YAP, WWOX, and NEDD4 has revealed a wide variety of different strategies that such a framework provides for that purpose (
). WWOX is a particularly interesting example as it takes the way of cooperation to its extreme by optimizing two domains for a different aspect of the interaction, namely stability (WWOX WW2) and binding (WWOX WW1). In this study, we have investigated the details of this cooperation using a range of complementary biophysical approaches and modeling. Our first important finding is that while the isolated WW1 domain is unstable, WW2 stabilizes this domain through a network of interdomain contacts. As a consequence, despite the fact that WW2 lacks significant inherent affinity to PY3, the stabilized WWOX tandem domain binds the ErbB4 PY3 substrate with stronger affinity than an isolated WW1 domain (30 versus 78 μM; a ratio similar to previous reports (
Delineating WWOX protein interactome by tandem affinity purification-mass spectrometry: identification of top interactors and key metabolic pathways involved.
). Indeed, we too find that WW2 will not bind a single PPxY ligand independently, and that in tandem, it is outcompeted by WW1. However, in our second important finding, using ITC and NMR experiments, capable of detecting intermediate and weaker interactions, we demonstrate that WW2 does participate in the binding of peptides that contain a dual PPXY motif. Addition of a second PPXY motif consistently increased affinity by up to 10-fold, and this increase cannot be attributed to avidity, since the corresponding affinity increase for the single WW1 domain is less than twofold (78 μM for PY3 versus 50 μM for PY3PY3) (Fig. 4 and Table 2). This increase in binding affinity necessitates active involvement of the WW2 peptide–binding pocket despite its inherent low affinity for peptide substrates (because of the missing canonical tryptophan residue), presumably because of increased local effective concentration (Figure 4, Figure 5, Figure 6). It is true that such effects are seen in other systems, for example, tandem SH2 domains bind double phosphotyrosine peptides with >1000-fold higher affinity than a single domain (
), and the fibronectin (Fn) N-terminal domain contains five type 1 modules and binds short repeat motifs in the SfbI protein of Streptococcus pyogenes (
). However, in most of these and other cases, both domains bind the PPxY motif in isolation. Here, we find a dual role for WW2: A stabilizing influence on the partially unfolded WW1 domain as well as substrate binding in the case of double-PPxY motifs.
A molecular view of the parallel tandem-WW/dual-PPxY assembly
The essence of biological regulation by the interactions of WW domains and their PPxY ligands is the ability of a tandem WW domain to capture polypeptides containing proximal PPxY sequences. We explored the molecular factors governing this binding interaction for WWOX and a series of designed polypeptides. As expected, one factor is the inherent affinity of the PPxY sequence, adhering to previously reported relative binding strengths (PY3 > PY1 > PY2, determined by the presence of a key basic residue following the PPxY motif (
)). However, a second factor is the relative location of the two binding sequences. The connecting linker length is a well-recognized determinant of tandem binding affinities, as in polypeptides with multiple PPxY motifs binding predominantly occurs for proximal motifs (separated by ∼15–20 residues) and only rarely for distal motifs (
). These two factors embody enthalpic and entropic contributions, respectively, to binding. Here, we report a significant difference in the binding mode of the PY3PY3 and PY1PY3 peptides compared to PY1PY2, where WW2 is bound appreciably only by the former, whereas the latter binds predominantly to WW1 (Fig. 5). Replacement of the native linker by a short polyglycine sequence in PY1PY2S (as in PY1PY3 and PY3PY3) also leads to stronger interaction with WW2 (Figs. 5C and S5B; Tables 2 and 3).
Our third important finding is the establishment of a parallel binding mode, by investigating the effect of independently “knocking out” each of the PPxY motifs in the PY3PY3 peptide. The first PPxY interacts with WW1, and the second binds weakly to WW2 (Fig. 6). This binding mode is supported by a structural model of the WWOX tandem WW domain interacting with the PY3PY3 double-motif peptide (Fig. 7). Our model reconciles our experimental findings: (1) It reconfirms parallel binding of the peptide to WWOX, positioning the first motif into the first WW1 domain; (2) it provides an explanation for binding to WW2, where a strong electrostatic acidic patch attracts the basic region that flanks the binding motif. Thus, even without a canonical tryptophan, this domain can bind to a peptide motif brought into proximity by its adjacent second motif bound to WW1. (3) It is consistent with prearrangement in a peptide-binding compatible orientation. This conformation is able to bind the double-motif peptide in parallel orientation, where the two peptide motifs wrap around WWOX, bringing the two domains closer. Thus, residue T27 in WW1 at the domain interface, as well as the C-terminal region of WW2, is only perturbed by peptides binding also to WW2 (Fig. 6).
Of note, this model differs from a model suggested in a previous study, which was generated using molecular dynamics simulations starting from the template structure of FBP21 (
). That model suggested that the WW1 peptide–binding site is occupied by the WW2 domain and consequently, binding of peptide substrates would involve major domain rearrangement to allow access of the peptide to the binding site. Our NMR results do not provide any support for such a rearrangement, and moreover, regions that show significant change in surrounding when comparing the single and tandem WW domains (as reflected by our measured CSPs) are not located at the domain interface.
Beyond the main model presented in Figure 7 however, WWOX is functional as a dynamic ensemble of different conformations, as suggested from our NMR experiments, as well from previous studies (
). While different AlphaFold2 simulations converged on the same overall fold of the tandem WW domain, models that we generated using additional protocols (including Rosetta ab initio folding with GREMLIN constraints (
), data not shown) exhibit a range of additional possible relative orientations between WW1 and WW2, with predominantly similar conformation of the individual WW domains, but a wide variability for the linker and the tails. While only few comply with our experimental results, as an ensemble, they may set the stage for a more detailed investigation of the possible contribution of different conformations to distinct functional contexts. After all, even if some specific conformations dominate, it is often the ensemble of conformations of a protein that will define its function (
This highlights the advantages as well as the limitations of structural models of proteins and protein interactions generated by AlphaFold2: While the model of the WW1WW2 tandem domain bound to its substrate (Fig. 7) provides a structural explanation of our experimental findings, AlphaFold2 is limited in its scope and cannot provide a picture of the dynamics of the system (although attempts to increase diversity by manipulating the input multiple sequence alignment have been reported, e.g., (
)). All its models are similar: the isolated unstable WW1 domain adopts the same predicted conformation as within the tandem context, and the tandem domain models do not change significantly when adding single- or double-motif peptides PY3 and PY3PY3 (data not shown). This is in line with previous reports that AlphaFold2 often generates bound conformations even in the absence of the ligand (e.g., (
)). For similar reasons, in its quest for the most stable conformation, AlphaFold2 cannot be used to assess effects of point mutations, since it relates to point mutations as “local noise” and will by default converge on the same result (e.g., a very similar structure is generated for the complex of a tandem domain that reverts to the canonical tryptophan in WW2, Y85W; and peptides with different affinities cannot necessarily be distinguished).
Evolution of tandem domain binding: WWOX adopts a new WW-binding mode
This study and others have closely examined the structural factors governing binding affinities between tandem WW domains and their single (or dual) PPxY ligands. Comparison of tandem WW domains demonstrates that their interaction with their natural ligands is controlled by a multifactorial array of sequence- and structure-related parameters. This is schematically summarized in Figure 8A. Variants of the WW domain—stable, “unfolded,” and mutated (e.g., loss of tryptophan)—exhibiting different inherent PPxY affinities are connected by linkers differing in length and structural flexibility. We find WWOX to be relatively unstructured in WW1 (similarly to Su(dx)), lacking the canonical tryptophan in WW2 (as in Su(dx) and KIBRA), and connected by a flexible linker (as in FBP21, shorter than the YAP linker, Fig. 8B). A degree of coupling between the two domains is observed, with significant WW2-induced stabilization of WW1, although effects of PPxY binding to WW1 are not propagated to WW2. Despite these partial similarities, binding of dual PPxY by WWOX does not fully resemble any solved structure of a tandem WW domain–peptide interaction (Fig. 8C).
Figure 8Multilayer modulation of WW–ligand interactions and affinities.A, components of modulation including the nature of the WW domain (wildtype, mutant, “unfolded”) and the linker domain (length, flexible/rigid). B, schematic representation of known tandem-WW structures (not necessarily an exhaustive list) using the aforementioned building blocks. Shown are YAP (
). C, representation of WWOX, the subject of this study, using the same building blocks. The presence of WW2 exerts a stabilizing effect on WW1 via interactions with the linker. Similar stabilizing effects occur for other proteins (not shown for clarity). FBP21, formin-binding protein 21; WWOX, WW domain–containing oxidoreductase; YAP, Yes-associated protein.
Additional modes may be considered as well. Just as the Tondu domain–containing growth inhibitor protein bears one PPxY motif far from two “low-affinity” close motifs (
), the PY3 of ErbB4 may be the “binding initiation” site leading to interaction with the PY1 or PY2 site, which would otherwise exhibit low affinity. There are cases where tandem motif binding does not improve binding, as for example, in the YAP–LATS interaction, where the tandem motif binds with affinity similar to that of the single motif (
). This suggests that this kind of binding interplay does not always play a role in the regulation. Finally, post-translational modifications may also modulate affinity (
). The result is a range of seemingly similar tandem interaction modules actually capable of a wide range of affinities, and, commensurately, biological functions. Together with factors present in the ligand, such as variations in the binding motif and the intermotif distance, all these impact the balance of enthalpic and entropic factors that eventually determine the nature of the ensuing interaction between tandem WW domain and PPxY ligands.
We find this consistent from an evolutionary viewpoint—since WW–PPxY interactions (barring the missing tryptophan residue in some domains and occasionally leucines replacing the first proline of the PPxY motif) are relatively similar in structural determinants and consequently in affinity, other factors have evolved to allow a fine-tuning of these interactions in a setting-dependent manner. Here, WWOX provides a case in point, exemplifying that nature can create a less stable unstructured domain for modifying binding affinity (and possibly specificity) and compensate for this destabilization by an adjacent chaperone domain. This is in agreement with a large-scale analysis of the (calculated) stability of single domains versus their corresponding occurrence within a multidomain setting. Domains that lack independent stability were stabilized by favorable interactions with other domains, stabilization of such folds was optimized by evolution, and specific mutations at interdomain hotspots could rescue the overall stability (
), again underlining the importance of noncanonical factors as an evolutionary tool for fine-tuned control of this family of protein–protein interactions.
To summarize, we have mapped in this study the mutual influence of WW1 and WW2 tandem domains and their involvement in the binding of single and double motif–containing peptide substrates in WWOX, using a number of complementary experiments including ITC binding experiments, NMR CSP changes and titrations, structural modeling, and more. We have shown a dramatic increase in binding affinity of double motif–containing peptide substrates and demonstrated that WW2 can play an active and important role in the binding of a second peptide motif, beyond merely increasing WW1 stability. These results join and confirm others showing that WW domains have evolved to utilize several structural and sequence factors to modify the affinity for their canonical ligands. Further studies are needed to determine how this combination of two WW domains, one with reduced binding and one with reduced stability, contributes to defining the binding specificity, and consequently functionality of WWOX within cells.
Experimental procedures
Expression plasmids
PCR products WWOX WW1–WW2 (amino acids 16–91), WWOX ww1 (amino acids 16–50), and WWOX WW2 (amino acids 57–91) were cloned into pETM 30 YAP 171 to 264 plasmid (a kind gift of Dr Maria J. Macias, Institute for Research in Biomedicine, Barcelona) containing an N-terminal HisX6 tag followed by a glutathione-S-transferase and tobacco etch virus protease, instead of the YAP sequence in the NcoI and HindIII restriction sites.
Protein expression and purification
All the WWOX constructs were expressed in Escherichia coli BL21 pLysS cells (Novagen). Cells were grown in 2× YT medium. Induction was done at absorbance of 0.6 to 0.8 at 600 nm with 0.1 mM IPTG, and cells were grown overnight at 20 °C. 15N-labeled proteins for NMR measurements were expressed in M9 minimal medium (
) supplemented with 1 g/l 15NH4Cl. 15N,13C-double labeled samples were also grown with 2.5 g/l 13C6-d-glucose (only WWOX WW1–WW2) to an absorbance of 0.8 at 600 nm, induced with 0.15 mM IPTG and grown overnight at 27 °C. Cells were harvested by centrifugation and stored at −80 °C.
Cell pellets containing expressed WWOX constructs were resuspended in lysis buffer (50 mM NaH2PO4 pH 7.0, 300 mM NaCl, 10 mM imidazole, and 5 mM β-mercaptoethanol), supplemented with 1 mM PMSF and DNase. Cells were disrupted using a microfluidizer (Microfluidics). Lysate was cleared by centrifugation and subjected to 5 ml His-Trap columns (GE Healthcare). The protein was eluted with a linear imidazole gradient of 15 to 250 mM in 30 column volumes. Fractions containing the purified protein were pooled and dialyzed overnight at 4 °C against dialysis buffer (20 mM NaH2PO4 pH 7.0, 150 mM NaCl, and 5 mM β-mercaptoethanol) in the presence of tobacco etch virus protease. Cleaved protein was then subjected to a second round of His-Trap column, and flow-through containing the cleaved protein was collected. The proteins were further purified using 16/600 Superdex 75 pg size-exclusion chromatography columns (GE Healthcare) equilibrated in protein buffer (20 mM Tris pH 7.8 and 150 mM NaCl). NMR samples were prepared in 20 mM NaH2PO4 pH 6.8, 100 mM NaCl, 7% 2H2O, and supplemented with EDTA-free protease inhibitor cocktail (Roche). Final protein concentrations were 0.3 to 0.6 mM. All proteins were concentrated, flash-frozen in liquid nitrogen, and stored at −80 °C.
Peptide synthesis
Peptides were purchased from PHTD peptide and GenScript (HK Limited) with 90 to 95% purity.
CD
CD spectra of 50 μM WWOX WW1–WW2, WWOX WW1, and WWOX WW2 were recorded using a J-810 spectropolarimeter (Jasco) in protein buffer, in a quartz cuvette for far-UV CD spectroscopy. Far-UV CD spectra were collected in a spectral range of 190 to 260 nm. For measuring the spectra of interactions, 50 μM WWOX WW1–WW2 was incubated with 50 μM ErbB4 1291—1305 (PY3) or with 50 μM ErbB4 PY3PY3 in protein buffer. Background scans of the peptides were conducted in buffer alone and subtracted from the protein with peptide scans.
We note that our CD spectra of WWOX WW domain fragments differ from those reported previously by others (
). This is mainly because of the use of a slightly different construct (including only four extra residues at its N terminus) and the performance of the experiments at different buffers and temperatures.
Protein fluorescence
About 5 μM protein in protein buffer was incubated with or without 6 M urea. Protein tryptophans were excited at 295 nm, and emission was measured at 325 to 400 nm in 96-well plates with Cytation3 imaging reader (BioTek).
ITC
ITC measurements were performed on an ITC-200 or a PEAQ-ITC calorimeter (MicroCal, GE Healthcare) at 20 °C. About 1 to 2 mM PY peptides were titrated to 100 to 200 μM of the different WWOX constructs in protein buffer. The data were fitted using the analysis software PEAQ-ITC, to the single-site binding model, and N was set to 1 when the KD was weaker than 30 μM. The integrated peak of the first injection was excluded from the fit because of the large errors in the first step. The average value and the standard deviations of at least three independent experiments are reported. We note that previous studies have reported slightly different absolute binding affinities (
), mostly because of differences in buffer conditions, but the relative affinities are unchanged.
SEC–MALS
SEC–MALS experiments were performed with a pre-equilibrated analytical SEC column (Superdex 200 10/300 GL; GE Healthcare Life Sciences) with protein buffer, as described in the study by Mashahreh et al. (
NMR spectra were recorded on a DRX700 Bruker spectrometer using a cryogenic triple-resonance TCI probehead equipped with z-axis pulsed field gradients. Spectra were measured at 16.4 T and 286 K. 1H,15N-HSQC spectra for sample characterization and optimization were run for 30 to 40 min with acquisition times of 91.4 (70) ms and 1024 (100) complex points in the F2 (F1) dimensions, respectively. Triple-resonance HNCO, HNCA, HN(CO)CACB, and HNCACB spectra, using sensitivity-enhanced echo–antiecho detection, were acquired for uniformly 13C,15N-labeled samples. All triple-resonance and 15N-edited experiments were typically acquired with 40 to 48 complex points and an acquisition time of 20 to 24.1 ms in the 15N dimension, and with 1024 complex points and an acquisition time of 91.8 to 104.4 (145) ms in the observed proton (13C′) dimension. For indirect 13C dimensions, experiments with 13CO(13Cα) evolution were acquired with 60 (44) complex points and an acquisition time of 28.3 (7.8) ms, and experiments with 13Cα/β evolution were acquired with 44 to 52 complex points and an acquisition time of 4.1 to 4.8 ms. Peak assignment was based on these triple resonance spectra; in cases of ambiguity and/or significant peak broadening, the assignment was assisted by data acquired for a WW1 A35V mutant that exhibited reduced exchange broadening (data not shown). Processing and analysis of all spectra were performed using the TopSpin 3.2 package (Bruker BioSpin).
Binding of ErbB4 peptides to the proteins was monitored by repeating the 1H,15N-HSQC spectrum after serial additions of the desired peptide up to 2 to 4 M equivalents taking care to maintain a constant WW-domain concentration. Using these spectra, the affinity could be estimated by plotting the CSP as a function of the number of peptide equivalents added (
) (Fig. S6). In CSP experiments, changes of CSP >0.05 ppm (where CSP = √[ΔH2+(ΔN/5)2]) were considered significant, representing a change beyond 2.5 standard deviations of shifts in unaffected regions.
Generation of structural models of WWOX
Structural models of the tandem domain were generated with DeepMind AlphaFold2 multimer mode (
). The peptide was provided as a separate chain. AlphaFold, version 2.1.0 was used, with default parameters. The models presented in Figure 7 were generated on November 9, 2021 and are available at https://zenodo.org/record/6532331.
We also generated additional models using a range of different approaches, including RoseTTAFold (http://new.robetta.org (
). None of these models provided good agreement with all experimental results (data not shown).
Data availability
Structural models used to generate Figure 7 (WWOX WW1 WW2 tandem domain bound to PY3PY3) as well as corresponding models of the tandem domain bound to PY1PY2 and PY1PY2S are available on zenodo (https://zenodo.org/record/6532331). Raw data files from experiments that form the basis of the results presented in this study are provided upon request by the corresponding authors.
The authors declare that they have no conflicts of interest with the contents of this article.
Acknowledgments
We thank Amjad Farooq and Rami Aqeilan for providing the WWOX plasmids.
Author contributions
S. R.-B., J. F., D. E. S., J. H. C., and O. S.-F. conceptualization; S. R.-B., J. F., K. K.-A., T. T., O. A., J. H. C., and O. S.-F. methodology; T. T. and O. A. software; S. R.-B., J. F., K. K.-A., J. H. C., and O. S.-F. formal analysis; S. R.-B., J. F., and K. K.-A. investigation; S. R.-B., J. H. C., and O. S.-F. writing–original draft; J. F., K. K.-A., and D. E. S. writing—review & editing; S. R.-B., J. F., K. K.-A., T. T., O. A., J. H. C., and O. S.-F. visualization; J. H. C. and O. S.-F. supervision; J. H. C. and O. S.-F. project administration; J. H. C. and O. S.-F. funding acquisition.
Funding and additional information
This work was supported, in whole or in part, by the Israel Science Foundation, founded by the Israel Academy of Science and Humanities (grant number: 717/2017 [to O. S.-F.] and grant number: 964/19 [to J. H. C.]) and by the grant agreement of the European Research Council (grant number: 310873 [to O. S.-F.]).
Characterization of the structure and function of W--> F WW domain variants: identification of a natively unfolded protein that folds upon ligand binding.
Delineating WWOX protein interactome by tandem affinity purification-mass spectrometry: identification of top interactors and key metabolic pathways involved.