The isolated, purified TLDc domains of Ncoa7, Oxr1 and Tldc2 interact with the V-ATPase in GST pull-down assays.
As mentioned in introduction, the interaction between V-ATPase and Ncoa7 was originally identified in immunoprecipitation studies with anti-B1 subunit antibodies, Therefore we initially focused on details of the interaction between Ncoa7 and the V-ATPase B1 subunit. Mouse Ncoa7 is a large 943 amino acid long protein (Fig.1A). To study the molecular details of its interaction with V-ATPase, we first narrowed down its interaction site by dividing it into three approximately equal parts: 1) N-terminal Ncoa7 (2-353), which contains a LysM domain of unknown function, 2) middle Ncoa7 (354-592), without known motifs or domains, and 3) C-terminal Ncoa7 (593-943), which contains the TLDc domain and an upstream poly-E motif. The numbering is based on protein sequence NP_766083 of the Ncoa7 long isoform and boundaries of the domains are predicted using NCBI’s Conserved Domains Search tool (CD-Search) 26. These three regions of mouse Ncoa7, as well as its stand-alone TLDc domain (775-943), were expressed in bacteria with N-terminal GST- and C-terminal 6xHis-tags, purified and used as baits in GST pull-down experiments with mouse kidney protein lysates. Interaction with the V-ATPase B1 subunit was assessed by western blot analysis with anti-B1 subunit antibodies. The C-terminal part of Ncoa7 (593-943) interacted with the B1 subunit, while N-terminal and middle regions did not (Fig.1B). Then, using the same approach, we further narrowed down the interaction site to the isolated TLDc domain of Ncoa7 (775-943) (Fig 1B, last lane).
We then asked if the TLDc domains of Oxr1, Tbc1d24, Tldc1 and Tldc2 also interact with the V-ATPase. To answer this question, we first identified the TLDc domain boundaries in these proteins by performing multiple sequence alignments of mouse Ncoa7, Oxr1, Tbc1d24, Tldc1 and Tldc2. In addition, we examined human Tbc1d24, in which disease-causing mutations are well characterized, and zebrafish Oxr2, whose structure is solved 13 (Fig. S1). Based on this alignment, the amino acids that correspond to P775 found at the start of the conserved TLDc domain of Ncoa7, are P698 in Oxr1, H336 in Tbc1d24, D235 in Tldc1 and P44 in Tldc2 (numbering is based on the longest isoforms of all proteins, Fig. S1). The conserved regions extend to the very end of Ncoa7, Oxr1 and Tldc2, while Tbc1d24 and Tldc1 contain additional non-conserved amino acid sequences at the C-terminus, and Tbc1d24 contains a long non-conserved loop within its TLDc domain (Fig. S1). In our initial experiments, these non-conserved regions were included in the TLDc domain constructs to better resemble the native proteins. Thus, Oxr1 (698-866), Tbc1d24 (336-561), Tldc1 (235-455) and Tldc2 (44-212), representing TLDc domains of the corresponding proteins, were expressed in bacteria, purified and used in GST pull-down experiments with mouse kidney lysates. We found, that the TLDc domains of Oxr1 and Tldc2, but not Tbc1d24 or Tldc1 interacted with the V-ATPase in these experiments (Fig. 1C). Using the same approach, we also found that even the full-length purified Tldc1 construct (2-455) did not interact with the V-ATPase in GST pull-down experiments (Fig. 2B, last lane). The absence of interaction of Tbc1d24 or Tldc1 and the V-ATPase was unexpected, so we performed additional GST pull-down experiments with truncated versions of the TLDc domains of Tbc1d24 and Tldc1 and also co-immunoprecipitation studies.
Truncated versions of Tbc1d24 and Tldc1 TLDc domains, lacking the non-conserved insertion and non-conserved C-terminal extensions, respectively, did not interact with V-ATPase.
As mentioned above, Tbc1d24 contains a long non-conserved loop (447-495 a.a. in mouse Tbc1d24), in addition both Tbc1d24 and Tldc1 (but not Ncoa7, Oxr1 or Tldc2) have non-conserved C-terminal sequences, extending beyond their TLDc domains (557-561 a.a. in mouse Tbc1d24 and 411-455 a.a. in mouse Tldc1, Fig. 2A). The indicated boundaries are based on the alignment between all 5 TLDc-domain containing mouse proteins (Fig.S1). We hypothesized that these additional regions in Tbc1d24 and Tldc1 might play a regulatory role and prevent interaction with the V-ATPase under our experimental conditions. Therefore, truncated versions of Tbc1d24 and Tldc1, referred to as Tbc1d24 (336-446_496-556) and Tldc1 (235-410), were produced and used as baits in GST pull-downs. However, similar to the non-modified TLDc domains, these truncated versions still did not bind to the V-ATPase (Fig. 2B). These results suggest that the non-conserved extensions in the TLDc domains of Tbc1d24 and Tldc1 are not the explanation for the lack of interaction of these two TLDc domains with the V-ATPase.
Tbc1d24 protein, but not Tldc1 and Tldc2, is detected in the kidney; endogenous kidney Tbc1d24 and exogenous Tldc1 and Tldc2 overexpressed in HEK293T cells interact with V-ATPase in co-immunoprecipitation studies.
In parallel with GST pull-down experiments with the isolated purified TLDc domains, we also studied interaction between the V-ATPase and endogenous full-length Tbc1d24, Tldc1 and Tldc2 by co-immunoprecipitation of these proteins from mouse kidney lysates. Mouse kidney lysate was chosen for these studies initially, because kidneys express high levels of V-ATPase and we found previously that Ncoa7 and Oxr1 are expressed and interact with V-ATPase in kidney 9. In preliminary experiments the specificity of commercially available antibodies against Tbc1d24 was confirmed by western blot analysis of lysates of mouse cortical collecting duct M-1 cells transfected with a plasmid expressing mouse Tbc1d24 with or without Tbc1d24-specific siRNAs (Fig 3A). These anti-Tbc1d24 antibodies were able to detect the endogenous Tbc1d24 in the kidney by western blotting (Fig. 3B, input lane) and in kidney intercalated cells in both the cortical (Fig. 3C) and medullary collecting ducts (Fig. 3D) by immunofluorescence. We then performed a co-immunoprecipitation experiment using anti-B1 and anti-B2 antibodies and mouse kidney lysate. Interestingly, Tbc1d24 was detected only in the anti-B1 co-immunoprecipitated material, demonstrating its interaction with V-ATPase and also a specificitypreference towards interaction with the B1 subunit of V-ATPase in the kidney (Fig. 3B).
We also confirmed the specificity of commercial anti-Tldc1 and anti-Tldc2 antibodies, using lysates of HEK293T cells overexpressing mouse Tldc1 (Fig. 4A) or Tldc2 (Fig. 4B). However, these antibodies were not able to detect endogenous levels of Tldc1 (Fig. 4A) or Tldc2 (Fig. 4B) in total or medullary mouse kidney lysates. This could be due to a low level of expression of Tldc1 and Tldc2 in the kidney, or their expression could be limited to a very small number of cells in the kidney. Because proteins can be concentrated during immunoprecipitation and sometimes reach detectable levels, we proceeded with studies using anti-B1 and anti-B2 antibodies for immunoprecipitation, followed by western blotting with anti-Tldc1 and anti-Tldc2 antibodies. However, we were not able to detect Tldc1 or Tldc2 in the co-immunoprecipitated material using kidney lysates (Fig. S2).
Since Tldc1 and Tldc2 were not detected in kidney by western blotting, possibly due to their low expression levels, we overexpressed them in HEK293T cells. We then studied their co-immunoprecipitation with V-ATPase using kidney lysate, mixed with the lysates of HEK293T cells, overexpressing Tldc1 or Tldc2, using anti-B2 and anti-B1 antibodies. In this way, we were able to show that both Tldc1 (Fig. 4C) and Tldc2 (Fig. 4D) interacted with the B1 subunit of V-ATPase. In addition, Tldc1 but not Tldc2 interacted with B2 subunit of V-ATPase (Fig. 4C and 4D), demonstrating that Tldc2 but not Tldc1 interacts specificallypreferentially with the B1 subunit of V-ATPase.
A poly-E rich motif, located upstream of the TLDc domain in Ncoa7 and Tldc2, enhancesexhibit a (nonsignificant) trend to enhance their interaction with the V-ATPase.
Both Ncoa7 and Tldc2, but not Oxr1, Tbc1d24 or Tldc1, contain a distinct poly-E rich motif upstream of their TLDc domains (Fig. S1 and 5A). We hypothesized that this poly-E rich motif could potentially contribute to the interaction of Ncoa7 and Tldc2 with the V-ATPase. To explore this possibility, we expressed and purified the following constructs: 1) Ncoa7 (593-750), lacking both the poly-E rich motif and the TLDc domain, 2) Ncoa7 (593-774), containing the poly-E rich motif but lacking the TLDc domain, 3) Ncoa7 (593-943), containing both the poly-E rich motif and the TLDc domain in addition to a portion of the C-terminus, 4) Ncoa7 (751-943), containing only the poly-E rich motif and the TLDc domain, 5) full-length Ncoa7 (775-943), 6) Tldc2 (44-212), lacking the poly-E motif, and 7) full-length Tldc2 (2-212), which consists almost entirely of the poly-E rich motif and the TLDc domain (Fig. 5A). These constructs were then used as bait in GST pull-down assays with mouse kidney lysates. Both constructs lacking the TLDc domain, Ncoa7 (593-750) and Ncoa7 (593-774), did not interact with the V-ATPase, indicating that the presence of a poly-E rich motif without a TLDc domain is not sufficient to produce significant interaction with the V-ATPase (Fig. 5B). In contrast, the constructs with both the poly-E motif and the TLDc domain, Ncoa7 (751-943), Ncoa7 (593-943) and Tldc2 (2-212), showed a trend toward stronger interactions with the V-ATPase in comparison with the TLDc domain-only constructs, Ncoa7 (775-943) and Tldc2 (44-212) (Fig. 5B and 5C). Thus,However, the differences were not significant at the P<0.05 level, due to high variability in the presence of a TLDc domain, the poly-E rich motif enhances interaction of both Ncoa7 and Tldc2 with the V-ATPase.data (Fig. 5C). Therefore, whilewe conclude that the TLDc domain is necessary and sufficient for interaction with the V-ATPase, and that the poly-E rich motif enhances this interaction, but is domain is not sufficient by itself to interact with the V-ATPase. It may, however, enhance the binding of the TLDc domain, although more work is needed to prove that this is a significant effect.
Alanine scanning mutagenesis of the evolutionarily conserved amino acids in Ncoa7 revealed that the G815A, G845A and G896A point mutants did not interact with V-ATPase, while the S817A, L926A and E938A mutations decreased the interaction.
Next, we asked which particular amino acids within the TLDc domain of Ncoa7 are important for its interaction with V-ATPase. Several amino acids in the TLDc domain of Ncoa7 are highly conserved across many species 10, such as G815, S817, G845, G896, L926 and E938 (numbering is based on protein sequence NP_766083 of the Ncoa7 long isoform). Substitution of each of these amino acids with alanine in the short Ncoa7-B isoform resulted in a less efficient protection against oxidative stress in comparison with wild type isoform, suggesting that these amino acids play an important role in this function of Ncoa7 10. Therefore, we examined whether these mutations also affect the ability of Ncoa7 to interact with V-ATPase, using the GST pull-down approach described above. In addition, a new alanine mutation of the non-conserved G802, predicted to be present in the unstructured loop on the surface of the Ncoa7 TLDc domain (Fig 6A), was used as a control mutation that we hypothesized would not affect Ncoa7 interaction with the V-ATPase.
The residues G802, G815, S817, G845, G896, L926 and E938 in mouse Ncoa7 are homologous to the residues G660, G673, S675, G703, G754, L784 and E796 in the zebrafish Oxr2, whose crystal structure was solved 13 (Fig. 6A). Based on the position of the homologous zebrafish amino acids in the structure, the Ncoa7 residues G802, G815 and S817 are exposed on the surface, while residues G845, G896, L926 and E938 are buried within the protein (Fig. 6A). Mutations of the buried residues are expected to affect the overall TLDc domain structure and stability. In agreement with that, when overexpressed in bacteria the G845A, G896A, L926A and E938A TLDc mutants were mostly insoluble, suggesting that the majority of the mutant protein molecules were not folded properly, and aggregated into the insoluble material (Fig. S3). However, some protein remained soluble and we were able to purify an amount sufficient to perform GST pull-down assays (Fig. S3). In contrast, G802A, G815A and S817A mutations on the protein surface did not affect the overall folding and, when overexpressed in bacteria, were very soluble (Fig. S4). All of these bacterially expressed Ncoa7 TLDc domain mutants were then purified and used as baits in GST pull-down experiments together with the wild type TLDc domain. Interestingly, in the GST pull down experiments, alanine mutations of surface exposed residues G815 and S817, which are located in the same short loop very close to each other on the structure (Fig. 6A), affected interaction with V-ATPase differently: while S817A only partially decreased the interaction, G815A completely abolished it (Fig. 6B and 6C). As expected, the G802A “control” mutation present in another loop further away from G815 and S817 did not decrease interaction with the V-ATPase at all (Fig. 6B and 6C). Alanine mutations of the buried G845, G896, L926 and E938 residues, that affected folding of Ncoa7 TLDc domain, all decreased interaction with V-ATPase, but to various degrees: G845A and G896A were most detrimental and no interaction between them and the V-ATPase was detected at all; interaction with L926A was barely detectable, while interaction with E938A was only partially decreased relative to the wild-type TLDc domain (Fig. 6B and 6C). This may be explained by the fact that in comparison to G845, G896 and L926, E938 is located closer to the C-terminus and to the surface of the protein. Therefore, the E938A mutation may not affect the overall 3D structure and interaction with V-ATPase as much as the three other mutations.
In addition, we performed a computational analysis of the effect of G802A, G815A, S817A, G845A, G896A, L926A and E938A mutations on Ncoa7 TLDc domain stability using PoPMuSiCv3.1 predictive software27, as described in methods. The calculated ΔΔG values predict the effect of mutations on thermodynamic stability of TLDc domain: the higher ΔΔG values indicate a greater destabilizing effect of the mutation on the TLDc domain folding. With the exception of L926A mutation, predicted ΔΔG values of less than 1 correlated with a medium to high degree of interaction between the mutants and V-ATPase, while values above 1 correlated with a loss of interaction. However, despite having the largest predicted ΔΔG value, the L926A mutant still interacted with V-ATPase, although very weakly (Fig. 6, Table 1). Overall, we found that the higher ΔΔG values strongly correlate with the loss of interaction of the mutants with the V-ATPase, further supporting the idea that Ncoa7 TLDc domain interaction with V-ATPase relies on the 3D folding stability of the TLDc domain (Table 1).
Table 1
Correlation between the predicted effect of the mutations on TLDc domain stability and interaction between the mutants and V-ATPase. ΔΔG, the change in the change in Gibbs free energy between the folded and unfolded states (ΔG folding ) and the change in ΔG folding when a point mutation is present, were calculated by PoPMuSiCv3.1. The smaller value of ΔΔG predicts smaller effect on TLDc domain stability. -, interaction undetectable; +++, interaction is similar to non-mutated (wild-type) protein; ++ and +, reduced and strongly reduced interaction.
Mutation: mouse Ncoa7 / zebrafish Oxr2
|
ΔΔG, kcal/mol
|
Predicted relative effect of the mutation on TLDc domain stability
|
Interaction with V-ATPase (in pull-downs)
|
G802A / G660A
|
0.84
|
small (<1)
|
+++
|
G815A / G673A
|
1.97
|
medium (>1, <2)
|
-
|
S817A / S675A
|
0.71
|
small (<1)
|
++
|
G845A / G703A
|
1.44
|
medium (>1, <2)
|
-
|
G896A / G754A
|
1.31
|
medium (>1, <2)
|
-
|
L926A / L784A
|
2.42
|
large (>2)
|
+
|
E938A / E796A
|
0.74
|
small (<1)
|
++
|
In conclusion, alanine mutations of the evolutionarily conserved amino acids G815, S817, G845, G896, L926 and E938 affected interaction with the V-ATPase differently. The G815A, G845A and G896A mutants of the Ncoa7 TLDc domain did not interact with V-ATPase at all and are, therefore, critical residues for this interaction. S817A, L926A and E938A mutations resulted in a decreased interaction with the V-ATPase, while the non-conserved control G802A mutation interacted with V-ATPase to a similar degree as the wild type construct. Alanine mutations of the buried residues G845, G896, L926 and E938 also affected proper folding of TLDc domain and decreased interaction with V-ATPase, suggesting that the overall native three-dimensional structure of the TLDc domain is necessary for its interaction with V-ATPase.