Biophysical characterization of SARAH domain–mediated multimerization of Hippo pathway complexes in Drosophila

Hippo pathway signaling limits cell growth and proliferation and maintains the stem-cell niche. These cellular events result from the coordinated activity of a core kinase cassette that is regulated, in part, by interactions involving Hippo, Salvador, and dRassF. These interactions are mediated by a conserved coiled-coil domain, termed SARAH, in each of these proteins. SARAH domain–mediated homodimerization of Hippo kinase leads to autophosphorylation and activation. Paradoxically, SARAH domain–mediated heterodimerization between Hippo and Salvador enhances Hippo kinase activity in cells, whereas complex formation with dRassF inhibits it. To better understand the mechanism by which each complex distinctly modulates Hippo kinase and pathway activity, here we biophysically characterized the entire suite of SARAH domain–mediated complexes. We purified the three SARAH domains from Drosophila melanogaster and performed an unbiased pulldown assay to identify all possible interactions, revealing that isolated SARAH domains are sufficient to recapitulate the cellular assemblies and that Hippo is a universal binding partner. Additionally, we found that the Salvador SARAH domain homodimerizes and demonstrate that this interaction is conserved in Salvador's mammalian homolog. Using native MS, we show that each of these complexes is dimeric in solution. We also measured the stability of each SARAH domain complex, finding that despite similarities at both the sequence and structural levels, SARAH domain complexes differ in stability. The identity, stoichiometry, and stability of these interactions characterized here comprehensively reveal the nature of SARAH domain–mediated complex formation and provide mechanistic insights into how SARAH domain–mediated interactions influence Hippo pathway activity.

The Hippo pathway controls tissue growth and differentiation by regulating the rate of cell proliferation and apoptosis (1)(2)(3)(4)(5)(6). Gene expression is silenced upon phosphorylation and cytoplasmic retention of the transcriptional co-factor Yorkie (1,(3)(4)(5)(6)(7)(8). Aberrant pathway activity leads to overgrowth phenotypes in both mammalian and fly models and is associated with tumorigenesis in humans (1,(3)(4)(5)(6)(7)(8)(9)(10)(11). Hippo signal transduction relies on the activity of a core kinase cassette that includes two kinases, Hippo and Warts, and two accessory proteins, Salvador and Mats. These proteins, together with Yorkie, form a core kinase cassette that is responsible for gene regulation. Each component of this core cassette has a mammalian homolog with similar function. The mechanisms regulating the activity of this kinase cassette are also most likely conserved as disruption of the genes encoding either Hippo (hpo) or Mats (mts) in flies results in overgrowth phenotypes that can be rescued by complementation with the mammalian genes encoding the homologous proteins, Mst1/2 or Mob1, respectively (5,7,10).
The activity of the core cassette begins with the activation of Hippo kinase, or Mst1/2 in mammals, which involves phosphorylation of the activation loop (12, 13). Genetic, cellular, and in vitro data suggest that homodimerization mediated by the conserved C-terminal coiled-coil domain, termed SARAH, of either Hippo or Mst1/2 promotes autophosphorylation in trans and kinase activation (6,(14)(15)(16)(17)(18)(19). Activated Hippo then phosphorylates Warts. This phosphorylation event leads to Warts activation, which also involves autophosphorylation and complex formation with Mats (5,(20)(21)(22). Warts kinase can then phosphorylate the transcriptional co-factor Yorkie, leading to the accumulation of phosphorylated Yorkie in the cytoplasm (7,8).
SARAH domains were first identified as a conserved coiledcoil region located at the C termini of three members of the Hippo pathway, Salvador, dRassF, and Hippo, and were named after each (23). Dimeric, trimeric, and tetrameric assemblies of coiled-coils have been proposed, but the preponderance of evidence supports a dimeric assembly with limited evidence for tetrameric assemblies and no evidence for trimeric ones (14,18,(23)(24)(25)(26)(27)(28)(29)(30)(31)(32)(33)(34). Structural studies have each revealed a conserved arrangement of two antiparallel coiled-coils with each mono- cro ARTICLE mer adopting a short 3/10-helix followed by a turn and a long ␣-helix (18,(25)(26)(27)(28)(29)(30)(31)(32). The dimer interface is mediated by interstrand interactions along the length of each ␣-helix and provides a rationale for why each SARAH domain-mediated complex is mutually exclusive. The SARAH domain of Salvador contains an extended N-terminal region that forms an additional ␣-helix that enlarges the dimer interface with Hippo (31,32), and Salvador variants lacking this extension bind less tightly to Hippo SARAH domain (32). Together, these data suggest that Hippo:Salvador complexes are tighter than other SARAH domain-mediated complexes. Analysis of the interactions of SARAH domain complexes have largely focused on the mammalian homologs of Hippo and dRassF with little to no data for Salvador SARAH domain, and reported dissociation constants vary widely, from nanomolar to millimolar (29,30,(35)(36)(37). There has been no systematic study of the entire set of possible SARAH domain-mediated interactions, making it difficult to understand the likelihood of any given complex forming.
Salvador and dRassF bind Hippo via their respective SARAH domains and yet have opposing roles in regulating the activity of Hippo kinase, with Salvador stimulating pathway activity and dRassF inhibiting it (3,24). The mechanisms by which SARAH domain-mediated complex formation influences Hippo kinase activity are still a matter of investigation. The importance of these SARAH domain proteins to Hippo pathway activity is underscored by the observation that flies bearing alleles encoding variants of either Hippo or Salvador that lack SARAH domains display overgrowth phenotypes typical of pathway disruption (4,6). Likewise, combination of mutant alleles of dRassF and Hippo lacking SARAH domains leads to a stronger overgrowth phenotype than either allele separately (24). Each SARAH domain-mediated complex correlates with a distinct phosphorylation state of Hippo. Immunoprecipitation assays revealed that unphosphorylated Hippo bound dRassF and phosphorylated Hippo bound Salvador (24). Evidence suggests two possible routes for the inhibition of the pathway by dRassF. Complex formation between Hippo and dRassF would prevent homodimerization of Hippo and thus block autophosphorylation and thus activation. dRassF could further promote the unphosphorylated state of Hippo by acting as a scaffold to recruit the phosphatase-containing complex dSTRIPAK to Hippo (38). The functions of the mammalian dRassF homologs are more complex, showing either inactivation or activation of Mst1/2, depending on cell type and the level of Hippo pathway activity (16,39). Co-expression of the human homolog of Salvador, hSalvador, with Mst1/2 leads to higher levels of activation loop phosphorylation (5,31,32). The increase in phosphorylation of Mst1/2 can be partly attributed to hSalvador inhibiting the phosphatase activity of the STRIPAK complex (31) but does not address which molecular event promotes the initial phosphorylation of Hippo. Interpretation of existing cellular and genetic data are complicated by an incomplete understanding of how the biophysical characteristics, such as stoichiometry and stability, of each SARAH domain-mediated complex are different. Because these parameters influence the likelihood of complex formation, they also influence the activity of the pathway.
To rationalize the function of SARAH domains in modulating Hippo pathway activity, we undertook a systematic study of the biophysical parameters of complex formation by the entire suite of isolated SARAH domain complexes. We chose to work with the Drosophila SARAH domains because they are a simpler system than their mammalian counterparts; there is only one homolog of each Drosophila SARAH domain compared with multiple isoforms of the human homologs of Hippo and dRassF. We began by purifying each SARAH domain individually to systematically study the entire set of SARAH domain complexes. We report here in vitro binding studies that revealed that a previously unidentified interaction between Salvador SARAH domains that we determine is dimeric, mediates association of the full-length protein in cells, is less stable than other SARAH domain interactions, and is conserved in hSalvador. We also find that the binding patterns of isolated SARAH domains recapitulate those previously observed in cells; the specificity of these interactions, therefore, is derived from the SARAH domains and not from other protein domains or cellular factors. We also report biophysical analysis of SARAH domain interactions and show that each of the interactions is dimeric in solution, but they have varying stabilities. Our results provide biophysical constraints that define the nature of SARAH domain-mediated complex formation and have implications for the molecular mechanisms of these complexes on pathway activity.

The Hippo SARAH domain is a universal binding partner
We previously determined the domain boundaries for the Drosophila melanogaster Salvador SARAH domain (32), which then enabled us to purify each of the three SARAH domains from D. melanogaster, Salvador, Hippo, and dRassF ( Fig. 1).
To determine which complexes could be formed by isolated SARAH domains, we performed an unbiased in vitro pulldown assay using purified proteins. Untagged SARAH domains were Coomassie-stained SDS-polyacrylamide gel of purified SARAH domain proteins following gel-filtration chromatography. From left to right are dRassF, Hippo, and Salvador followed by complexes between Salvador:Hippo and dRassF:Hippo SARAH domains. The migration of each SARAH domain is indicated by a cartoon rectangle with green for dRassF, blue for Hippo, and orange for Salvador.

Biophysical characterization of SARAH domains
incubated with hexahistidine (H6) 2 and thioredoxin-tagged SARAH domains immobilized on nickel-affinity resin; because the SARAH domains have similar mobility on SDS-PAGE, the thioredoxin tag was added to change the migration of one SARAH domain in each pulldown to help with visualization. The resin was collected and washed, and the amount of untagged SARAH domain bound was monitored by Coomassie-stained SDS-PAGE (Fig. 2). Hippo SARAH domain bound each SARAH domain tested, Hippo, Salvador, and dRassF. In contrast, Salvador SARAH domain and dRassF SARAH domain bound only themselves or Hippo SARAH domain. No interactions were detected between Salvador and dRassF SARAH domains, regardless of which protein was the bait or prey in the experiment. To ensure complexes were homogeneous in future assays, SARAH domain heterocomplexes were co-expressed and purified using tandem-affinity purification (Fig. 1).

The Salvador SARAH domain oligomerizes
Because we observed an interaction between Salvador SARAH domains in the pulldown assay (Fig. 2), we asked whether this SARAH domain also mediated homooligomerization of the full-length (FL) protein in cells. Because hSalvador homodimerizes through an atypical WW domain (40,41), binding between differentially tagged FL proteins was used as a positive control. Binding between a FL and a variant containing the isolated SARAH domain was used to detect SARAH domain-mediated homooligomerization. HEK293 cells were transiently co-transfected with plasmids encoding FLAGtagged Salvador FL and variants corresponding to either Salvador FL or the SARAH domain tagged with the HA epitope and monomeric yellow fluorescent protein (mYFP). Complexes were isolated from cell lysates by immunoprecipitation using an anti-HA antibody, and complex formation was monitored by Western blotting (Fig. 3A). FLAG-tagged Salvador FL co-immunoprecipitated with both HA-mYFP Salvador FL and HA-mYFP Salvador SARAH domain. These interactions demonstrate that the SARAH domain of Salvador is sufficient to mediate homooligomerization in cells. To determine whether this interaction is conserved across species, we performed an analogous set of co-immunoprecipitation experiments with hSalvador (Fig. 3B). The HA-mYFP-tagged SARAH domain of hSalvador was able to immunoprecipitate FLAG-tagged hSalvador FL, indicating that homooligomer- 2 The abbreviations used are: H6, hexahistidine; IMAC, immobilized metal affinity; chromatography; FL, full-length; mYFP, monomeric yellow fluorescent protein; K u , equilibrium constant of unfolding; N, native; U, unfolded; ⌬G u , Gibbs free energy of unfolding; ⌬G H2O, Gibbs free energy in water; TCEP, tris(2-carboxyethyl)phosphine hydrochloride; MRE, molar residue ellipticity; R, universal gas constant; T, temperature; Y obs , observed signal; Y n , signal from native protein; F n , fraction native; Y u , signal from unfolded protein; F u , fraction unfolded; a n , y intercept for the native protein; a u , y intercept for the unfolded protein; b n , slope of the baseline for the native region of the sigmoid; b u , slope of the baseline for the native region of the sigmoid; P t , total concentration of protein; HA, hemagglutinin; SUMO, small ubiquitin-like modifier.

Biophysical characterization of SARAH domains
ization mediated by the SARAH domain of Salvador is conserved across species.

SARAH domains are dimeric
We next wanted to determine the stoichiometry of each of the five possible SARAH domain complexes identified (Hippo: Hippo, Salvador:Salvador, dRassF:dRassF, Hippo:Salvador, and Hippo:dRassF) (Fig. 1). We choose to employ native MS as it analyzes the assembly of complexes in physiological buffers without the addition of exogenous tags or chemical cross-linkers that may either disrupt physiologically relevant complexes or stabilize nonphysiologically relevant ones. Each complex was analyzed by electrospray ionization MS under native and denaturing conditions ( Fig. 4 and Fig. S1). For each spectrum, the m/z values of the observed peaks were deconvoluted to the molecular weight, and the identity of the complex was assigned based on those weights ( Table 1).
The highest oligomeric species detected under native conditions for each of the five complexes were dimers-homodimers in the case of Hippo, Salvador, and dRassF SARAH domains and heterodimers for complexes between Hippo and either Salvador or dRassF SARAH domains. Under denaturing conditions, only monomeric species were detected for the Hippo, dRassF, or Hippo:dRassF SARAH domain samples. The charge state distribution of the monomeric peaks was at a lower m/z (higher charge) than those in the native spectra, suggesting that the denaturing buffer reduced the stability of the monomeric SARAH domains. In contrast, peaks corresponding to either homodimers of Salvador SARAH domain or heterodimers of Hippo:Salvador SARAH domains were maintained under denaturing conditions and had the same charge state distribution as in the native conditions, suggesting that homo-and heterodimers containing Salvador SARAH domain were not disrupted by the denaturing buffer. Most likely the hydrophobic

SARAH domain dimers have varying stabilities
To understand the likelihood of each SARAH domainmediated complex forming, we determined the relative stabilities of each using chemically induced equilibrium unfolding. For SARAH domains, complex formation and domain unfolding are coupled; disassociation of homodimers resulted in unfolded monomeric proteins for homodimers of both Mst1 and a mammalian dRassF homolog as well as a heterodimer between the two (29,30,35). Therefore, the fraction of folded protein corresponds to the number of complexes formed, and the equilibrium constant of unfolding (K u ) varies according to the disassociation constant (K d ) (42). In each of the five SARAH domain dimers, the amount of folded protein was determined by measuring the ␣-helical signal by far-UV CD at 222 nm as a function of urea concentration (Fig. 5). For each titration, we ensured that equilibrium conditions were satisfied for each point and that unfolding curves corresponding to denaturing and renaturing experiments were superposable. In each case, the unfolding curve was sigmoidal with a single transition, implying that the unfolding reaction involved a simple twostate transition; SARAH domains went from folded to unfolded without populating any significant intermediate conformations. At low urea concentrations, the unfolding curves for both Salvador and dRassF SARAH domains did reveal a possible second transition. However, in both cases, the data corresponding to a possible intermediate represented a sufficiently small fraction of the system that an intermediate state could not be accurately modeled.
To determine the thermodynamic constants that describe the stability of each SARAH domain dimer, we fit the data to a two-state model of unfolding in which a native folded dimer (N 2 ) transitions to two unfolded monomers (2U), N 2 º 2U. K u can be expressed in terms of the ratio of concentrations of unfolded and folded protein, K u ϭ [U] 2 /[N 2 ], or in terms of free energy of unfolding (⌬G u ), K u ϭ exp(Ϫ⌬G u /RT). ⌬G u can be linearly extrapolated to determine the Gibbs free energy in the absence of urea (⌬G H2O ) for each complex, which has the advantage of allowing for comparison of stabilities between different proteins without being influenced by the denaturant (43). For each of the five SARAH domain complexes, we performed titrations at three different protein concentrations (Fig.   Figure 4. SARAH domains complexes are dimers. Peaks corresponding to either monomeric (left) or dimeric (right) species in the mass spectra for Hippo (A), Salvador (B), dRassF (C), Hippo:Salvador (D), and Hippo:dRassF (E) SARAH domains. The charge states are indicated above selected peaks as are cartoons representing the identity of the species based on the molecular weight assignments with the Hippo SARAH domain represented by blue diamonds, Salvador SARAH by orange, and dRassF by green. Insets correspond to a zoomed-in view of the selected regions of spectra. For A, C, and E, the monomeric and dimeric spectra correspond to the denatured and native spectra, respectively. For B and D, the monomeric and dimeric spectra are from the denatured spectra. The native spectra for B and D are shown in Fig. 1.

Table 1 Predicted and measured molecular weights by native MS of SARAH domain compexes
Biophysical characterization of SARAH domains 5). The midpoint of the transition for each unfolding curve varied with protein concentration, demonstrating that the observed transitions corresponded to the disassociation of the dimers. Therefore, to accurately determine the thermodynamic parameters, data from three titrations performed at different protein concentrations were fit simultaneously using nonlinear least-squares fitting and global analysis ( Table 2 and Fig. 5).
The stabilities of the SARAH domain complexes varied by nearly 10 kcal/mol (

Discussion
Our ability to purify each Drosophila SARAH domain protein allowed us to investigate the biophysical nature of complex formation for the entire suite of SARAH domains and enabled us to directly determine and compare key biophysical parameters among all possible protein complexes (Fig. 1). We started by performing an unbiased pulldown assay using purified SARAH domains to determine which complexes are mediated by isolated SARAH domains. Our results show that the Hippo

Biophysical characterization of SARAH domains
SARAH domain is a universal binding partner, interacting with itself as well as the SARAH domains from dRassF and Salvador. We found that the binding patterns observed in cells can be recapitulated using isolated SARAH domains, implying that specificity arises from the SARAH domains and not additional factors, such as other protein domains or co-factors (Fig. 2). These binding experiments also revealed an unanticipated complex, a homooligomer of the Salvador SARAH domain that we determine also mediates homooligomerization of FL protein in cells in both the fly and human homologs (Figs. 2 and 3). Using native MS, and therefore in solution and free from the addition of any chemical cross-linkers, we determined the stoichiometry of each of the five possible SARAH domain complexes, as identified in our pulldown assay. We show that each complex is dimeric and find no evidence for higher-order complexes ( Fig. 4 and Table 1). The stoichiometries determined here agree with those assigned to homologous mammalian complexes, demonstrating that the organization of SARAH domain complexes is conserved (14,18,(25)(26)(27)(28)(29)(30)(31)(32).
The complexes mediated by the SARAH domain of Hippo, both with itself as well as with the SARAH domains of Salvador and dRassF, influence the activation state of Hippo and, in turn, overall pathway activity. As each possible complex is mutually exclusive, there will be competition between homo-and heterodimer formation (24). To rationalize the effects of each SARAH domain complex on Hippo signal transduction, we need to understand the likelihood of complex formation. We determined the stabilities of each of the five SARAH domain complexes using chemically induced equilibrium unfolding ( Fig. 5 and Table 2). We observed that the midpoints of the transitions in the unfolding curves were dependent on the concentration of the SARAH domains. This dependence confirmed that disassociation of SARAH domain dimers resulted in unfolded monomers and that differences in stabilities directly correlate with differences in affinity (29,30,35,42) (Fig.  5). Because the midpoint of the transitions varied with protein concentration, to accurately determine the thermodynamic parameters associated with unfolding, we needed to globally fit data from a set of unfolding curves collected at multiple protein concentrations. Each of the five possible SARAH domain complexes were fit to a two-state model of unfolding (N 2 º 2U) ( Table 2), reflecting both the single transition of each unfolding curve and the stoichiometries determined by native MS.
Despite having relatively similar structures, the stabilities of the five SARAH domain complexes vary widely, over 10 kcal/ mol. The Hippo:Salvador SARAH domain heterodimer is the most stable (⌬G H2O ϭ 17.7 Ϯ 0.4 kcal/mol) and twice as stable as the weakest complex, the Salvador SARAH domain homo-dimer (⌬G H2O ϭ 7.8 Ϯ 0.4 kcal/mol). The increased stability of the Hippo:Salvador SARAH domain heterodimer is consistent with the complex having a larger buried surface area than other SARAH domain dimers and supports the observation that, by extending the binding interface with Hippo, the N-terminal extension of the Salvador SARAH domain leads to a tighter complex (32). Other factors may also contribute to this elevated stability, such as a partially unfolded intermediate, which was observed for the Mst1 SARAH domain (35). The remaining SARAH domain complexes, Hippo:Hippo, dRassF:dRassF, and Hippo:dRassF, cluster between these two extremes with stabilities that vary only by ϳ1.5 kcal/mol. The stabilities of these three SARAH domain complexes also align with the relative sizes of the buried surface areas observed in the structures of their mammalian homologs (32). Because the rank orders of both the m values and ⌬G H2O of the Drosophila SARAH domain complexes match the rank order of buried surface areas observed in the structures of the mammalian homologs (32), the hierarchy of stabilities of SARAH domain complexes is most likely conserved between the Drosophila and mammalian homologs.
The varying stabilities of SARAH domain complexes provide insight into the likelihood of their formation, which, in turn, will have distinct effects on Hippo kinase activity. The common consensus is that SARAH domain-mediated homodimerization of Hippo leads to autophosphorylation and kinase activation, yet the stability of the Hippo homodimer (⌬G H2O ϭ 11.9 Ϯ 0.2 kcal/mol) is nearly 6 kcal/mol lower than the heterodimer with Salvador (⌬G H2O ϭ 17.7 Ϯ 0.4 kcal/mol). These observations suggest that the presence of Salvador shifts the equilibrium from Hippo homodimers toward Hippo:Salvador heterodimers. Given the difference in stabilities, how would Hippo kinase initially become phosphorylated in the presence of Salvador? Phosphorylation by an upstream kinase, Tao-1, could represent one possible route of Hippo activation. However, inactivation of the TAO kinase genes in mammalian cells does not fully block Hippo pathway activity, indicating that this upstream kinase is not fully responsible for activation of Mst1/2 (45)(46)(47). These results coupled with our analysis of SARAH domain stabilities suggest that there may still be additional cellular events that promote activation of Hippo kinase. The stability of the Hippo homodimer is 1 kcal/mol lower than the heterodimer with dRassF (⌬G H2O ϭ 13.3 Ϯ 0.3 kcal/mol), a negative effector of Hippo activation (24). This difference in stability provides a thermodynamic explanation for the inhibitory effects of dRassF. The presence of dRassF would shift the equilibrium toward dRassF:Hippo heterodimers and away from Hippo homodimers that can undergo autophosphorylation.

Biophysical characterization of SARAH domains
Because the Salvador:Hippo SARAH domain heterodimer is the most stable, perhaps arising from the expanded binding interface mediated by the N-terminal extension of the Salvador SARAH domain (32), it is perplexing that the Salvador SARAH domain homodimer, which contains two N-terminal extensions, is not the most stable complex; instead, the Salvador homodimer is the weakest. To understand this phenomenon, we first compared the structure of the Salvador SARAH domain of D. melanogaster (Fig. 6A) with the human homolog. The organization of the two SARAH domains is highly conserved; they superpose with a root mean square deviation of 1.1 Å over 66 C␣ atoms. Therefore, we hypothesized that the Salvador SARAH domain would likely adopt a similar organization when part of a SARAH domain-mediated homodimer. A conformation other than a coiled-coil has been reported for a monomeric SARAH domain in a larger protein assembly, suggesting that these domains may sample more conformational space than previously appreciated (48). To gain insight into the possible structural organization of a Salvador SARAH domain homodimer, we performed superpositions using the Hippo:Salvador SARAH domain structure as a template. We first superposed the main ␣-helix of Salvador SARAH domain onto the equivalent region of the Hippo-SARAH domain so that intrastrand interactions along the coiled-coil were maintained (Fig. 6B). This superposition resulted in steric clashes between the N-terminal extension of one monomer and C-terminal end of the opposing monomer, rendering this arrangement implausible. We also performed a second superposition in which we aligned the C terminus of Salvador SARAH domain onto the equivalent region of Hippo so that the packing with the N-terminal extension from the other SARAH domain in the dimer would be maintained (Fig. 6C). This superposition resulted in the two ␣-helices becoming splayed apart at the opposite end of the coiled-coil. In this arrangement, only a fraction of the typical interstrand interactions could form. This arrangement is consistent with both the lower stability and smaller m value observed for the Salvador SARAH domain homodimer ( Table 2). Despite this lower stability, the Salvador SARAH domain homodimer is sufficiently tight to be detected in both in vitro and cell-binding experiments (Figs. 2 and 3). We show that this interaction is conserved between flies and mammals (Fig. 3), but the biological implications of Salvador having two points of homodimerization, the SARAH domain and the atypical WW domain, is not yet apparent (40,41).
Early studies established the link between SARAH domainmediated interactions and pathway activity, but how these interactions modulate Hippo kinase activity is still a topic of investigation. Our work provides the biophysical constraints, in the form of specificity, stoichiometry, and stability, for the entire suite of SARAH domain complexes that are necessary to formulate a mechanistic framework for the interplay between SARAH domain complexes. The results presented here provide a mechanistic rationale for how dRassF inhibits Hippo kinase activation. However, these results also raise questions as to the mechanism governing Hippo kinase autophosphorylation in the presence of Salvador and suggest that there may be additional routes to the activation of Hippo kinase. The next step in understanding the regulation of the Hippo kinase cassette will be interpreting how these coiled-coil interactions function within the context of full-length proteins and multicomponent, physiological assemblies.

Protein expression and purification
Nucleotides encoding the SARAH domains of D. melanogaster sav (residues 530 -608), hpo (residues 606 -662), and dRassF (residues 754 -802) were each cloned into a pBAD4 vector derivative downstream from a H6 and SUMO tag (49). Proteins were overexpressed in Escherichia coli and purified as described previously (32). Briefly, proteins were isolated using nickel-charged IMAC, the SUMO tag was removed by incubation with a SUMO-specific protease (SENP) at 4°C overnight,

Biophysical characterization of SARAH domains
and the protein was further purified by anion-exchange and gel-filtration chromatographies.
For protein complexes, genes encoding the Salvador:Hippo SARAH domains and dRassF:Hippo SARAH domains were cloned into pRSF-Duet Vector (EMD-Millipore) as SUMOtagged fusions with an additional SBP tag for Hippo and an additional H6 tag for either dRassF or Salvador. Proteins were isolated by sequential affinity purification and IMAC followed by Streptactin (GE Healthcare) and were further purified as described previously (32). All SARAH domain proteins were flash-frozen in 10 mM Tris, pH 8, 200 mM NaCl, 1 mM Tris(2carboxyethyl)phosphine hydrochloride (TCEP).
For proteins expressed as fusions to thioredoxin, genes encoding the SARAH domains of D. melanogaster sav (residues 530 -608), hpo (residues 606 -662), or dRassF (residues 754 -802) were cloned into a derivative of the pBAD4 vector downstream from H6 and thioredoxin tags (49). T7 Express cells (New England BioLabs) transformed with the plasmids encoding either Hippo or dRassF SARAH domains were grown at 37°C, and protein expression was induced at 0.5 A 600 with 0.5 mM isopropyl ␤-D-1-thiogalactopyranoside at 37°C for 3 h. Rosetta 2 (DE3) (Novagen) cells transformed with a vector encoding Salvador SARAH domain were grown at 37°C, and protein expression was induced at 0.5 A 600 with 0.5 mM isopropyl ␤-D-1-thiogalactopyranoside at 20°C for 16 h. Cells were lysed in 40 mM Tris, pH 8, and 400 mM NaCl, 0.1% Nonidet P-40 and supplemented with protease inhibitor mixture (Sigma). For Salvador SARAH domain, the lysis buffer was further supplemented with 10% glycerol. The thioredoxin fusion proteins were purified by sequential IMAC and anion-and gel-filtration chromatographies. Final proteins were stored in 10 mM Tris, pH 8, 200 mM NaCl, 1 mM TCEP and flash-frozen in liquid nitrogen at 10 mg/ml.

In vitro pulldown experiments
25-l reactions containing a 100 M concentration of the indicated untagged SARAH domains and a 25 M concentration of the indicated H6-thioredoxin-tagged SARAH domain in 10 mM Tris, pH 8.0, 200 mM NaCl, and 1 mM TCEP were incubated for 1 h at room temperature. To promote more robust binding, binding reactions with untagged dRassF SARAH domain were modified to include 500 M untagged dRassF SARAH domain and a reaction temperature of 30°C. To each binding reaction, 10 l of nickel-charged Profinity IMAC resin (Bio-Rad) was added, and the reaction was further incubated for 30 min at room temperature. The resin was collected and washed three times with 500 l of 10 mM Tris, pH 8.0, and 200 mM NaCl. The resin was boiled in SDS-loading buffer. Samples were analyzed by SDS-PAGE stained with Coomassie Brilliant Blue. Gels were imaged using the Odyssey IR imaging system (LI-COR). Band intensities were quantified in ImageJ (50). The amount of untagged SARAH domain bound above background in each complex was calculated by subtracting the amount of untagged protein bound to the H6-thioredoxin control; to normalize between replicates, the difference was divided by the amount of untagged protein bound to H6-thioredoxin for each set. Prism software (GraphPad Soft-ware, La Jolla, CA) was used to plot data and for statistical analysis. Ratios were analyzed by a paired t test.

Co-immunoprecipitation experiments
Nucleotides encoding the FL (residues 1-608) D. melanogaster sav were cloned into pCDNA 3.1 (Invitrogen) downstream from either the FLAG tag or tandem HA and mYFP (51) tags. Additionally, nucleotides corresponding to the SARAH domain (residues 530 -608) from D. melanogaster sav were also cloned downstream of HA and mYFP tags. An equivalent set of plasmids were cloned using the nucleotides encoding either the FL (residues 1-383) or SARAH domain (residues 299 -383) of SAV1 from Homo sapiens. HEK293T cells (ATCC) were cultured in Dulbecco's modified Eagle's medium (Gibco) supplemented with 5% fetal bovine serum (VWR) and 2 mM glutamine at 37°C and 5% CO 2 . Cells were transfected with the indicated plasmids using X-tremeGENE HP DNA transfection reagent (Sigma-Aldrich) according to the manufacturer's protocol. Cells were harvested 24 h following transfection and lysed in ice-cold lysis buffer (20 mM Tris, pH 8, 150 mM NaCl, 1% Nonidet P-40, 10% glycerol) supplemented with 1 mM phenylmethylsulfonyl fluoride, 1 mM NaVO 4 , 5 mM NaF, 2.5 mM Na 2 P 2 O 7 , 0.1 mM ␤-glycerophosphate, 5 mM EDTA, Protease Inhibitor Mixture (Sigma), and Universal Nuclease (Pierce). Normalized, clarified lysates were incubated with Protein G resin and HA antibody (Roche Applied Science, lot 34502100) for 2 h at 4°C. Resin was collected, washed with lysis buffer, and boiled in SDS-loading buffer. Samples were analyzed by Western blotting using primary antibodies that recognized the HA (Roche Applied Science, lot 34502100) or FLAG (Sigma-Aldrich, lot SLCD3990) epitopes followed by either IRDye 800CW goat anti-rat (LI-COR, lot C90122-02) or IRDye 800CW goat anti-mouse (LI-COR, lot C81106-01) secondary antibodies, respectively. Blots were scanned on an Odyssey IR imaging system (LI-COR). Bands were quantified using ImageJ (50). The relative amount of FLAG-tagged protein bound was determined by first correcting for background binding by subtracting the intensity of the FLAG-tagged protein bound in the negative control, lysates lacking HA-mYFP, from each experiment. Relative binding was then determined by calculating the ratio of FLAG-tagged protein in each lane to the amount bound by the positive control, HA-mYFP Salvador FL. Ratios were analyzed by a paired t test. Prism software (GraphPad Software) was used to plot data and for statistical analysis.

Native mass spectrometry
MS analysis was carried out as described previously (52) using a Waters SYNAPT G2-Si instrument. Immediately prior to analysis, purified SARAH domain complexes were buffer-exchanged into 200 mM ammonium acetate, pH 7.6, using a 3K molecular weight cutoff filter (Pall Life Sciences). Before native analysis, SARAH domain samples were bufferexchanged and then diluted 1:100 (v/v) into 200 mM ammonium acetate for a final concentration of each SARAH domain of 0.023 mg/ml for Hippo, 0.04 mg/ml for Salvador, 0.01 mg/ml for dRassF, 0.1 mg/ml for Hippo:Salvador, and 0.012 mg/ml for Hippo:dRassF. Concentrations were chosen to maximize the ratio of signal to noise of the spectra. Sam-

Biophysical characterization of SARAH domains
ples were then infused from in-house prepared gold-coated borosilicate glass capillaries into the electrospray source at ϳ90 nl/min with the following settings: source temperature 30°C, capillary voltage 1.75 kV, argon collision cell flow rate 7 ml/min, transfer collision energy 5 V, trap collision energy 20 V. Comparisons between samples were only made in cases where settings were identical. For denatured analyses, each sample was diluted 1:10 (v/v) into 200 mM ammonium acetate, pH 7.6, and then 1:10 (v/v) into 50% methanol plus 1% formic acid or 50% methanol plus 0.5% formic acid for the Salvador SARAH domain. Denatured samples were analyzed using the same settings as used for the native analysis. Data analysis was performed using the Waters MassLynx version 4.1 software with the Maximum Entropy add-on.

Urea-induced equilibrium denaturation
CD data were collected with an Aviv model 420 spectropolarimeter (Lakewood, NJ) equipped with a temperature control cell holder and an automatic syringe titrator (Hamilton Microlab 500 titrator). Measurements were collected at 25°C in a cuvette with a 1-cm path length at 222 nm and averaged over 30 s. Titrations were performed at the indicated protein concentrations, where protein concentration corresponds to molar amount of total protein chains, and the sample was equilibrated for 300 s between injections. Native samples were equilibrated in 10 mM Tris, pH 8.0, 200 mM NaCl, 1 mM TCEP (Thermo), and denatured samples were equilibrated in the same buffer supplemented with ϳ8 M urea (VWR). Urea was purified using mixed bed resin (Bio-Rad) (53), and the concentration of urea stock solutions was determined by refractive index (43). Three denaturing titrations and one renaturing titration were performed for each sample.

Global fitting of denaturation data using a two-state model
Raw CD data were converted to molar residue ellipticity (MRE) (degrees ϫ cm 2 /dmol) to normalize signal arising from differences in protein concentration, and fits were performed on MRE values. SARAH domain unfolding was fit to a dimeric, two-state folding model in which a native folded dimer (N 2 ) transitions to two unfolded monomers (2U) (29,30,35). This process can be expressed as follows. where R is the Universal gas constant and T is the temperature in K. ⌬G u is linearly dependent on the concentration of urea (43) and can be expressed as follows, ⌬G u ϭ ⌬G H2O ϩ m[urea] (Eq. 4) where ⌬G H2O is the change in the Gibbs free energy of the SARAH domain dimer in water and m is the slope of the urea dependence.
The observed signal for each data point (Y obs ) can be expressed as the sum of the signal (Y n ) arising from the fraction native (F n ) and the signal (Y u ) from the fraction unfolded (F n ).
͑Y obs ϭ ͑1 Ϫ F u ͒Y n ϩ F u ϫ Y u ͒ (Eq. 5) The native and unfolded signals are linearly dependent on the concentration of urea and can be rewritten as follows, Y n ϭ a n ϩ b n [urea] (Eq. 6) Y u ϭ a u ϩ b u [urea] (Eq. 7) where a n and a u are the y intercepts, and b n and b u are the slopes of the baselines corresponding to the native and unfolded regions of the sigmoid, respectively. The total concentration of protein (P t ) is equal to the sum of the concentrations of the native and unfolded chains.  By substituting Equation 9 into Equation 2, F u can also be expressed as follows.
Isolation of F u allows for substitution of Equation 10 into Equation 5, resulting in an expression in which Y obs depends on urea and which can be used to determine ⌬G H2O .
Three data sets, at varying protein concentrations, for each SARAH domain were fitted globally using the nonlinear fitting function in Mathematica (Wolfram Research), constraining ⌬G H2O and m, while treating a n , b n , a u , and b u for each concentration as locally adjustable parameters (54).
To enable easy comparison among the sigmoids and fits from different SARAH domains, the data and fits in the figure were converted to fraction folded using the equation, using Y u and Y n values obtained from the fits. This normalization does result in variations in the baseline regions that were not present in the raw data.

Structural analysis
Superpositions were performed using least-squares fitting in CCP4 (55) and visually inspected in Coot (56). Atomic model representations were generated in PyMOL (Schrodinger, LLC).

Data availability
Raw MS data have been deposited in the Mass Spectrometry Interactive Virtual Environment (MassIVE) with the project identifier MSV000085022. All other data contained within the article are available upon request from the corresponding author, Jennifer Kavran (jkavran@jhu.edu).