Increased expression of BCH-motif-containing molecule at the C-terminal region 1 (BMCC1) correlates with a favourable prognosis in neuroblastoma, but the underlying mechanism remains unknown. We here isolated BNIPXL (BNIP2 Extra Long) as a single contig of the extended, in-vitro-assembled BMCC1. Here, we show that in addition to homophilic interactions, the BNIP2 and Cdc42GAP homology (BCH) domain of BNIPXL interacts with specific conformers of RhoA and also mediates association with the catalytic DH-PH domains of Lbc, a RhoA-specific guanine nucleotide exchange factor (RhoGEF). BNIPXL does not recognize the constitutive active G14V and Q63L mutants of RhoA but targets the fast-cycling F30L and the dominant-negative T19N mutants. A second region at the N-terminus of BNIPXL also targets the proline-rich region of Lbc. Whereas overexpression of BNIPXL reduces active RhoA levels, knockdown of BNIPXL expression has the reverse effect. Consequently, BNIPXL inhibits Lbc-induced oncogenic transformation. Interestingly, BNIPXL can also interact with RhoC, but not with RhoB. Given the importance of RhoA and RhoGEF signaling in tumorigenesis, BNIPXL could suppress cellular transformation by preventing sustained Rho activation in concert with restricting RhoA and Lbc binding via its BCH domain. This could provide a general mechanism for regulating RhoGEFs and their target GTPases.
Introduction
Rho proteins are part of the extensive Ras superfamily of small GTPases that cycle between inactive GDP-bound and active GTP-bound states. Despite their structural and biochemical similarities, they target different binding partners and regulate cell morphology, motility, vesicular transport, membrane trafficking, lipid signaling, cell cycle progression and gene transcription (Etienne-Manneville and Hall, 2002; Jaffe and Hall, 2005). The functions of Rho, Rac and Cdc42 have been extensively characterized (Hall, 1998). Similarly to all GTPases, they are tightly regulated by three distinct classes of protein, namely guanine-nucleotide exchange factors (GEFs), GTPase-activating proteins (GAPs) and guanine-nucleotide dissociation inhibitors (GDIs) (Bos et al., 2007). GEFs promote the exchange of GDP for GTP, GAPs stimulate intrinsic GTPase activity and GDIs inhibit nucleotide dissociation and GTP hydrolysis by sequestering Rho in inactive complexes. These regulators modulate Rho activities in response to spatio-temporal cues and may utilise their multiple domains to coordinate crosstalk between different signaling pathways (Rossman et al., 2005; Tcherkezian and Lamarche-Vane, 2007). Dysfunctional regulation of Rho signaling unsurprisingly leads to cancer (Sahai and Marshall, 2002; Ridley, 2004) and neurological disorders (Boettner and van Aelst, 2002). In addition to regulation by lipid metabolites (Rossman et al., 2005) and intrinsic phosphorylation (Moon and Zheng, 2003), it remains to be explored how the activity of these regulators on their cognate GTPases are regulated by additional cellular protein(s).
We previously identified BNIP2 and Cdc42GAP homology (BCH) domain as a conserved multi-functional protein domain, which, in addition to its ability to form homophilic and heterophilic interactions with itself or other homologous BCH domains (Low et al., 1999; Low et al., 2000a; Shang et al., 2003), has emerged as a new class of regulatory GTPase-binding domain. For example, the BNIP2 BCH domain activates Cdc42 to elicit cell protrusions (Zhou et al., 2005) whereas the BNIP-Sα BCH domain sequesters Cdc42GAP/p50RhoGAP via heterophilic BCH interactions while targeting RhoA directly for further activation, leading to apoptosis (Zhou et al., 2002; Zhou et al., 2006). Furthermore, the BCH domain of BPGAP1, a close homolog of Cdc42GAP/p50RhoGAP, also promotes Cdc42-dependent short pseudopodia, which are necessary for cell migration (Shang et al., 2003; Lua and Low, 2004) and Ras/MAPK activation (Lua and Low, 2005). More recently, we demonstrated that the BCH domain of the brain-specific BNIP-H/Caytaxin is involved in glutamate production and glutaminase trafficking (Buschdorf et al., 2006). These findings suggest that the BCH domain exerts an additional regulatory control over the cellular activity of their cognate small GTPases or other specific proteins. Identification and systematic studies of other BCH-domain-containing proteins will therefore provide further insights into the mechanistic functions of this emerging class of protein module.
We isolated the cDNA for BNIPXL (BNIP2 extra long) encoding a protein distinct from BNIP2 (Accession number: NP_004321). The extended isoform, BCH-motif-containing molecule at the C-terminal region 1 (BMCC1), was recently identified as a prognostic marker for childhood neuroblastomas (Machida et al., 2006). BMCC1 enhances neuronal apoptosis upon NGF depletion through as yet unknown mechanisms. Here, we show that BNIPXL could regulate actin cytoskeletal reorganization and cellular transformation. Its BCH domain specifically targets RhoA and RhoC, but not RhoB, whereas separate pools could exist to sequester their RhoGEF, Lbc. Together, such interactions suppress primary foci formation in NIH3T3 fibroblasts. This confirms that the BCH domain acts as a regulatory GTPase-binding domain by targeting GTPases and their regulators, and could represent a general mechanism by which BNIPXL functions as an antagonist of cellular transformation.
Results
BNIPXL as a novel BCH-domain-containing protein
To identify novel BCH-domain-containing proteins, bioinformatics analyses by PSI-BLAST were conducted using full-length Homo sapiens BNIP2 and its BCH domain. Based on the earliest entry KIAA0367 (Nagase et al., 1997), we isolated BNIPXL cDNAs that encode a 769 amino acid (BNIPXLα; accession number, AY439213) and a 732 amino acid protein (BNIPXLβ; accession number, AY439214). The shorter variant results from alternative splicing of exons 11 and 12, which introduces an in-frame stop codon (supplementary material Fig. S1). Multiple sequence analyses of BNIPXL with four other BCH-containing proteins that we previously identified – BNIP2, BNIP-Sα, BNIP-H/Caytaxin and BPGAP1 – indicates the closest homology to BNIP2 across the entire protein (58% identity, 74% similarity) (Fig. 1A), whereas its BCH domain is most similar to that of BNIP-H/Caytaxin (76% identity, 90% similarity) (Fig. 1B and supplementary material Fig. S2). To prevent crosshybridization arising from their significant similarities to BNIP2 and other BCH-domain-containing homologs, and to effectively distinguish BNIPXLα expression from BNIPXLβ, semi-quantitative RT-PCR using specific primers spanning the alternative splice sites demonstrate that both transcripts are expressed in all tissues and cell lines examined except in kidney HEK293T cells, where only the BNIPXLβ isoform was detectable (supplementary material Fig. S3). Recently, Machida et al. (Machida et al., 2006) assembled several cDNA fragments, which together encode a much-extended form, BMCC1. However, in our cell types and culture conditions, we were unable to isolate any sequences upstream of BNIPXL in a single contig or detect the endogenous 350 kDa protein expected of BMCC1 in HEK293T cells, as previously described (Machida et al., 2006). A further MEGABLAST (optimized for highly similar sequences) search using the BMCC1 coding sequence as query against the human expressed sequence tag (EST) database revealed that a 1.5 kb region upstream of BNIPXL was not represented by any ESTs (data not shown). These results imply that BNIPXL could represent another isoform of the BMCC gene family that is differentially expressed. Here, we set out to explore the mechanistic roles of the BMCC gene family by directly investigating how the prototypical BNIPXL engages its cellular protein partners to regulate cell growth.
BNIPXL BCH domain associates with RhoA and RhoC, but not RhoB, and elicits cell morphological changes
Our previous work indicates that the BNIP2 BCH domain activates Cdc42 for cell protrusions (Zhou et al., 2005) whereas BNIP-Sα targets RhoA and Cdc42GAP/p50RhoGAP, leading to concerted RhoA activation and apoptosis (Zhou et al., 2002; Zhou et al., 2006). Phylogenetic analyses indicate that the BNIPXL BCH domain is highly homologous to that of BNIP2 and BNIP-Sα (supplementary material Fig. S2), raising the possibility that BNIPXL could also associate with Rho subfamily GTPases and regulate their biological functions. To test this hypothesis, cells were transfected with Flag-BNIPXL (Fig. 2A) and HA-RhoA or HA-Cdc42 and immunoprecipitated with anti-Flag beads. Bound proteins were resolved for anti-HA western analyses. The results confirmed that both BNIPXL isoforms could indeed interact with RhoA but not with Cdc42 (Fig. 2B). RhoA, RhoB and RhoC are three highly related (sharing 85% identity), but play distinct roles in cells (Wheeler and Ridley, 2004). RhoA participates in oncogenic transformation, whereas RhoC promotes tumor metastasis and cell migration (Wang et al., 2003; Wu et al., 2004; Ridley, 2004; Hakem et al., 2005). By contrast, RhoB is predominantly involved in endosome trafficking (Fernandez-Borja et al., 2005). To examine the functional role(s) of BNIPXL, we determined the specificity of its BCH domain toward these Rho proteins. Coimmunoprecipitation studies with HA-Rho isoforms showed that the CBCH (contains BCH) fragment of BNIPXLβ is highly specific to RhoA and RhoC (Fig. 2C). No binding was observed with other regions of BNIPXL (Fig. 3). Finally, we show that the CBCH fragment but not the N-terminus WBCH (without BCH) or the shorter NBCH (no BCH) fragments of BNIPXL could readily bind to full-length BNIPXLβ (Fig. 2D). This observation confirms that the BNIPXL BCH domain, like all other BCH domains, forms homophilic complexes, adding to the versatility of this protein-protein interaction module. Collectively, these results support the notion that the BNIPXL BCH domain is a novel Rho-targeting domain and is likely to be involved in Rho-mediated cellular processes. Consistently, HeLa cells expressing full-length BNIPXLβ and the CBCH fragment induced short membrane protrusions to which they are highly localized – in addition to being perinuclear and punctate in the cytosol (Fig. 2E, arrows). By contrast, the WBCH fragment presented a uniform cytosolic distribution with no distinct change in cell morphology. Since BNIPXL targets RhoA, such subtle morphological changes imply the possible involvement of RhoA downstream of BNIPXL.
BNIPXL binds to specific RhoA conformers
We further investigated the molecular basis and functional significance of BNIPXL-Rho interaction by focusing on RhoA and BNIPXLβ (since similar results were observed for BNIPXLα). To confirm that BNIPXLβ-RhoA interaction was indeed direct, in vitro transcribed and translated Flag-BNIPXL was incubated with equal amounts of bacterially produced GST or GST-RhoA after guanine nucleotide depletion. The results indicate that the CBCH but not the WBCH fragment forms direct and enhanced interactions with RhoA, whereas full-length BNIPXL interactions with RhoA could be detected only after prolonged film exposure (Fig. 3A). However, no signals were detected for the negative control and WBCH fragment even after prolonged exposure, suggesting that the BNIPXL N-terminus could exert an autoinhibitory effect, restricting BNIPXL-RhoA interaction. Next, we examined whether BNIPXL could interact with different RhoA conformers. GST-RhoA was unloaded or loaded with GDP or GTPγS (a non-hydrolyzable GTP analogue), and incubated with lysates expressing BNIPXL WBCH, CBCH fragments or p160ROCK, a RhoA effector (Ishizaki et al., 1996) used as a positive control for GTP-RhoA loading (Fig. 3B). Unlike p160ROCK binding, which is strictly specific to active GTP-Rho, the BNIPXL CBCH fragment has no distinct preference for the nucleotide-free or GTPγS-RhoA, although less so for GDP-RhoA. The BNIPXL WBCH fragment was not precipitated by any of the RhoA conformers. These results suggest that the BCH domain could be sensitive towards specific RhoA conformations in vivo, requiring specific regions or motifs. To probe this further, BNIPXL interaction with RhoA constitutive active and dominant negative mutants was examined in vivo. The constitutive active G14V and Q63L mutants or the `fast-cycling' F30L mutant cause a net increase in active RhoA. G14V and Q63L are locked in an active conformation mimicking GTP-RhoA in vitro, because residues for GAP recognition in the P-loop (residues 12-19) and switch II (residues 62-68) regions have been mutated (Dvorsky and Ahmadian, 2004; Bos et al., 2007). By contrast, F30L is spontaneously active and demonstrates enhanced GDP to GTP exchange, while retaining regular GTP hydrolysis rates catalyzed by GAPs (Lin et al., 1999). The dominant negative T19N mutant on the other hand, sequesters and inhibits multiple endogenous RhoGEFs, preventing RhoA activation leading to elevated GDP-RhoA levels in vivo. Coimmunoprecipitation of cells transfected with Flag-BNIPXL, CBCH or WBCH fragments and HA-RhoA wild-type or mutants show that the full-length BNIPXL and the BNIPXL CBCH fragment associate weakly with wild-type RhoA, and this was greatly enhanced with the F30L or T19N mutants (Fig. 3C). In strong contrast to the GDP/GTP loading experiments, the G14V or Q63L mutants did not interact with BNIPXL or CBCH, despite being locked in the GTP-bound conformation. Although such preference is not observed in vitro, these results are likely to reflect the greater affinity between BNIPXL and GDP-RhoA in vivo. As a negative control, the WBCH fragment fails to associate with any RhoA conformers. This is consistent with a survey of RhoA structural motifs indicating the importance of the P-loop region for effector binding and GAP recognition (Dvorsky and Ahmadian, 2004). Interestingly, the same regions are important for GEF recognition and nucleotide dissociation (Bos et al., 2007). Taken together, these results demonstrate that the BNIPXL BCH domain can detect RhoA conformational changes, most likely by targeting its P-loop and switch II regions.
Interaction of BNIPXL with RhoA requires an intact BCH domain
We next examined whether the BNIPXL BCH domain contains a unique motif(s) for binding RhoA. Our previous report indicates that the BNIP-Sα BCH domain contains a RhoA-binding domain (RBD) (residues 148-177) (Zhou et al., 2006). Multiple sequence alignments of class I Rho-binding motifs found in several RhoA effectors (Fujisawa et al., 1998) and the BNIP-Sα RBD motif suggests that the proximal region of the BNIPXL BCH domain (residues 615-644) could be important for RhoA binding (supplementary material Fig. S4). A series of three BNIPXL internal deletions (designated Δ0, Δ1 and Δ2) flanking the putative RhoA-binding motif was generated after considering its secondary structure (Fig. 4A, panel i). To minimize disruption of important folds, residues outside the putative α-helices were chosen to define the deletions. Coimmunoprecipitation of lysates expressing these deletions and wild-type RhoA demonstrates that any deletion spanning these regions disables BNIPXL-RhoA binding (Fig. 4A, panel ii). To ensure that the structural integrity and native cellular distribution of these deletions were not altered, they were tested for their ability to form homophilic interactions with wild-type BNIPXL (Fig. 4B). Similarly, these mutants retained a cytosolic distribution when transiently expressed in HeLa cells (supplementary material Fig. S5). These results therefore suggest that multiple motifs within the BCH domain are essential for RhoA recognition.
BNIPXL suppresses RhoA activity and stress fiber formation in vivo
The observation that BNIPXL binds specifically to distinct RhoA mutants, which have altered levels of steady-state GTP-RhoA implies that BNIPXL could directly or indirectly affect RhoA activity. To examine this, RBD pull-down assays were performed to determine active RhoA levels in HEK293T cells expressing the wild-type, Q63L or F30L RhoA mutants in the absence or presence of wild-type BNIPXL or its deletions. To show the efficacy of this assay (Fig. 4C, panel i), GST-RBD greatly precipitated both active Q63L and F30L compared with wild-type RhoA and did not bind T19N. The levels of active wild-type RhoA were greatly reduced in the presence of BNIPXL. Strikingly, the non-interactive mutants (Δ0, Δ1 and Δ2) could also reduce active RhoA levels, but only by 50% (Fig. 4C, panel ii). By contrast, the active Q63L level was not affected, because BNIPXL could not bind to this specific form (Fig. 3C). Consistently, wild-type BNIPXL reduced the overall level of active F30L, albeit less effectively than wild-type RhoA. This is probably because the fast-cycling mutant can counteract the inhibition with its greater intrinsic rate of GDP-GTP exchange. By comparison, the noninteractive mutants of BNIPXL failed to reduce the active F30L levels.
If BNIPXL can indeed reduce active RhoA levels, then suppression of BNIPXL expression should increase Rho-GTP loading. To verify this, short hairpin RNA (shRNA) was used to transiently knockdown BNIPXL expression in 293T cells. Of the six shRNAs used, the sh470 construct demonstrated effective knockdown of HA-BNIPXL expression without affecting endogenous RhoA levels (Fig. 4D, panel i). RBD pull-down assay was performed using lysates cotransfected with wild-type RhoA and sh470 or with two other non-targeting shRNA sequences (sh484 and sh47C). In strong contrast to the ability of BNIPXL to reduce GTP-RhoA levels, transient expression of sh470 but not sh484 and sh47C resulted in an increase in active RhoA levels. Collectively, these results suggest that BNIPXL acts as an upstream and negative regulator of RhoA.
RhoA activation promotes actin stress fibers in many cell types. The effect of BNIPXL and its deletion mutants on stress fiber formation was examined in serum-stimulated NIH3T3 fibroblasts. The extent of transfected cells exhibiting stress fibers was scored and indicated that 50% of cells expressing full-length BNIPXL had reduced phalloidin staining, consistent with its inhibition of RhoA activity. Representative fields of cells are shown in Fig. 4E, panel i. However, similar to the partial inhibition on RhoA activity in the RBD pull-down assays, we observed that a small percentage (25-30%; P<0.01) of cells expressing the deletion mutants that failed to interact with RhoA still exhibited some degree of inhibition on stress fiber formation at a lower potency (Fig. 4E, panel ii). Taken together, these results imply that BNIPXL could exert inhibition on RhoA activity directly by binding through its BCH domain and also via an additional mechanism that could be independent of such an interaction.
BNIPXL targets endogenous proto-Lbc, a RhoA-specific GEF, and suppresses its transforming activity
Since BNIPXL binds RhoA to reduce overall intracellular RhoA activity, it remains unclear whether such binding would only sequester RhoA, thus preventing activation by its upstream activators (such as Rho GEFs) or sequester both RhoA and the immediate RhoGEF from interacting with one another. Several indirect lines of evidence provide clues for the likelihood that BNIPXL could at least prevent RhoA binding to RhoGEF, while at the same time targeting a RhoGEF independently of RhoA binding. First, the BNIPXL BCH domain has a unique binding profile not commonly seen in any RhoGEFs or RhoGAPs. The ability of BNIPXL to recognize unbound and GTP-bound RhoA could indicate the possible sequestration of free RhoA from guanine-nucleotide exchange or for promoting subsequent GTP hydrolysis. Second, it associates strongly with RhoA-T19N, a mutant that serves to target RhoGEFs in a dominant-negative manner. Thus, BNIPXL could potentially compete with RhoGEF binding to RhoA. Third, although BNIPXL reduced GTP-RhoA levels and stress fiber formation, its deletion mutants that failed to interact with RhoA could still partially exert such inhibitory effects, indicative of a mechanism that is independent of direct RhoA binding. To address these possibilities, we adopted a candidate approach and examined whether BNIPXL could interact with the RhoA-specific GEF, proto-Lbc. The Lbc (lymphoid blast crisis) oncogene (Toksoz and Williams, 1994), Lfc (first cousin of Lbc) (Whitehead et al., 1995), and the mouse homologue of p115RhoGEF, Lsc (second cousin of Lbc) (Whitehead et al., 1996) are a subgroup of the Dbl (diffuse B-cell lymphoma) family GEFs that exhibit specificity toward Rho (Hart et al., 1994; Zheng et al., 1995; Glaven et al., 1996). Interestingly, our results show that BNIPXL and BMCC1 could associate with proto-Lbc as observed in coimmunoprecipitation studies (supplementary material Fig. S6). Interestingly, BNIPXL and BNIP-Sα are the only BNIP2 family members that target proto-Lbc (Fig. 5A). To determine whether BNIPXL interacts endogenously with proto-Lbc, we immunoprecipitated Flag-BNIPXL in 293T cells and assayed for the presence of endogenous proto-Lbc in the immunoprecipitates. We demonstrated that BNIPXL could indeed form a complex with endogenous proto-Lbc under normal growth conditions, albeit weakly (Fig. 5B).
To identify the binding regions on proto-Lbc, a series of proto-Lbc truncations were generated, as previously described (Sterpetti et al., 1999) (Fig. 6A), and used in coimmunoprecipitations with BNIPXL (Fig. 6B-E). Full-length BNIPXL bound to all the truncation mutants, although relatively weak interactions were detected with PS. Interestingly, different binding profiles were observed with the CBCH, WBCH and NCBH fragments. CBCH failed to interact with PS, indicating a requirement for the catalytic DH-PH domains, whereas the absence of the proline-rich motif in the PP, α-HEL and ΔC fragments resulted in a loss of binding to WBCH and NBCH. Furthermore, NBCH demonstrated weak interactions with CT when compared with the full-length and PS constructs. These results suggest that the 45 amino acid overlap present in WBCH, but not NBCH, is not essential for interaction with the proline-rich motif, but may instead, stabilize association of BNIPXL N-terminus with proto-Lbc. This is consistent with the observation that the CBCH domain, which contains the 45 amino acid overlap does not bind to PS. Collectively, our results demonstrate that full-length BNIPXL targets proto-Lbc at its catalytic DH-PH domains and the proline-rich motif (Fig. 6F).
Although both RhoA and RhoC demonstrate an ability to induce stress fibers, RhoA participates in oncogenic transformation differently to RhoC, which promotes cell migration and tumor metastasis (Aspenström et al., 2004). Oncogenic transformation is largely attributed to deregulated RhoGEF signaling (Rossman et al., 2005). The loss of C-terminal regulatory sequences in proto-Lbc results in an increasingly cytosolic distribution and oncogenicity, even though the catalytic activities of proto- and onco-Lbc in vivo remain comparable (Sterpetti et al., 1999). To determine the effects of BNIPXL on the transforming activity of onco-Lbc in NIH3T3 fibroblasts, transient transformation assays were performed using DsRed-onco-Lbc and GFP-BNIPXL (Fig. 6G). The CBCH domain abrogates onco-Lbc transformation whereas full-length BNIPXL exhibits moderate attenuation (25%) – a profile very similar to the enhanced binding of CBCH to RhoA when compared with full-length (Fig. 3A). This suggests that the CBCH domain alone is a potent inhibitor of onco-Lbc transformation. Interestingly, the WBCH domain also downregulates the transforming activity (55%). It is possible that the WBCH domain may function in parallel pathways to suppress Lbc via its binding to the proline-rich motif of Lbc, but independently of RhoA downregulation. Taken together, our results suggest that BNIPXL downregulates RhoA activity and transformation through sequestration of RhoA and Lbc, which together restrict sustained RhoA activation.
Discussion
Although Rho-GTPases serve as universal molecular switches for the cytoskeletal network and are tightly controlled by key regulators such as GAPs, GEFs and GDIs, it remains relatively unknown how the activity of these regulators on their cognate GTPases are regulated by additional cellular protein(s). In the current study, we have identified BNIPXL as a novel regulator for both RhoA and its specific RhoGEF, proto-Lbc, leading to RhoA inactivation and suppression of cellular transformation. As BNIPXL represents a single contig of BMCC1 that is a favorable prognostic marker for childhood neuroblastomas, this finding could also provide a general model by which BNIPXL or BMCC1 could function as an antagonist of RhoA- and Lbc-induced oncogenicity.
BNIPXL inhibits RhoA-Lbc interaction by separately binding to their different motifs
One unique feature of BNIPXL is its capacity to interact with Rho and RhoGEF, which both require at least the functional BCH domain. Our competition studies demonstrate that such bimodal interactions essentially deprive RhoA and Lbc from interacting with one another, thus effectively shutting off Lbc signaling to RhoA (data not shown). These distinctive interactions are mediated by the BCH domain recognizing RhoA conformational changes or through targeting the catalytic DH-PH domains of Lbc. An additional region(s) at the N-terminus of BNIPXL is also involved in binding to the proline-rich motif of Lbc. These interactions together provide a general mechanism to regulate Lbc-Rho pathways specifically (Fig. 7) and GTPase signaling in general. The motifs responsible for BNIPXL binding to RhoA are likely to be multiple ones, because short internal deletions have all led to complete loss of RhoA binding. Furthermore, the BCH domain exerted greater potency in suppressing Lbc transformation when compared with the full-length protein, consistent with its enhanced binding to RhoA.
BNIPXL is specific towards subtypes and conformational states of Rho
The unique binding profile of BNIPXL to specific RhoA mutants supports our model, which proposes that the BCH domain sequesters both RhoA and Lbc to limit their activation. The ability of BNIPXL to distinguish different RhoA conformations in vivo suggests that BNIPXL is likely to sequester and retain GDP-RhoA in a complex. Future work will determine whether Gly14 (within the P-loop) and Glu63 (within Switch II) could indeed constitute key binding residues for BNIPXL interaction. In addition, we cannot exclude the possibility that BNIPXL could directly promote intrinsic GTPase activity. Alternatively, it may do so through the recruitment of inactivators, such as the BCH-containing Cdc42GAP/p50RhoGAP or BPGAP1, through heterophilic BCH interactions. Similarly, because Lbc distribution contributes to sustained RhoA activation and oncogenic transformation (Sterpetti et al., 1999), BNIPXL could reduce ectopic GEF activation by affecting Lbc localization.
Although the role of Rho-mediated cytoskeletal rearrangements in cellular transformation remains unclear, early studies suggest that Rho promotes two independent pathways, leading to stress fiber formation and oncogenic Ras transformation (Qiu et al., 1995). Activation of p160ROCK, which is necessary for stress fiber formation (Amano et al., 1997), is required (Sahai et al., 1999), but not sufficient for Rho transformation. Rho contributes to transformation through cell growth regulation and survival (Jaffe and Hall, 2002). More significantly, cells transformed by RhoGEFs, including Lbc, retain normal stress fibers (Zheng et al., 1995). Collectively, these results suggest that BNIPXL targets distinct RhoA pathways, including stress fiber assembly and transformation, and is consistent with its function as an antagonist for cell growth in fibroblasts as reported here, or as a pro-apoptotic factor in neuronal cells (Machida et al., 2006). The identification of a specific BNIPXL mutant that uncouples GEF interactions while retaining RhoA binding should address the specific contributions of BNIPXL in linking actin cytoskeleton rearrangements and oncogenic transformation.
Intriguingly, the BNIPXL BCH domain specifically targets RhoA and RhoC, but not RhoB. Despite having similar subcellular distribution and the ability to induce stress fibers (Aspenström et al., 2004; Adamson et al., 1992), current literature suggests that RhoA participates in cellular transformation whereas RhoC promotes tumor metastasis and cell migration (Wang et al., 2003; Wu et al., 2004; Ridley, 2004; Hakem et al., 2005). Interestingly, Lbc demonstrates comparable GEF activity toward RhoA, RhoB and RhoC (Zheng et al., 1995), but their distinct functions remain unknown. Similarly, the role of RhoC in oncogenic Lbc transformation is not well defined and warrants further investigation, especially in light of its specificity in BNIPXL binding. It is interesting to note that BNIPXL overexpression is not associated with profound cell rounding despite Rho inactivation, suggesting that indirect activation of other GTPases such as Rac and Cdc42 lead to the subtle short protrusions observed in the current study. Such multitude of GTPases crosstalking downstream of BNIPXL awaits further investigation.
The BCH domain as a regulatory domain for small GTPases and their regulators
Our current results indicate that a number of motifs within the BCH domain mediate RhoA interaction, in contrast to the involvement of a single GTPase-binding motif within the BCH domains of BNIP2 and BNIP-Sα for targeting Cdc42 and RhoA, respectively. Significantly, only BNIPXL and BNIP-Sα interact with proto-Lbc amongst the four closely related BNIP2 family proteins. Our group has previously shown that BNIP-Sα can displace Cdc42GAP/p50-RhoGAP from RhoA to elicit RhoA activation and pro-apoptotic events. This implies that unlike BNIPXL, which acts by sequestering Lbc and RhoA separately, leading to RhoA inactivation, BNIP-Sα could instead provide a scaffold for Lbc to activate RhoA. This suggests that homologous BCH domains target specific GTPases and their regulators through a common sequestration mechanism to elicit different cellular outcomes. Collectively, our current findings provide further compelling evidence that the BCH domain is an emerging class of regulatory GTPase-binding domain.
DH-PH domains: a common target for specific BCH or Sec14p-like domains?
In addition to Lbc, BNIPXL also targets p115RhoGEF (supplementary material Fig. S7), another RhoA-specific GEF (Hart et al., 1996). By engaging many RhoGEFs, BNIPXL could exert varying levels of inhibition at different RhoA activation nodes. This leads to our hypothesis that BNIPXL may also target p115RhoGEF via the DH-PH domains. If so, the DH-PH domains could become a new target for selective members of the BNIP2 family, in addition to their small GTPase targets. Intriguingly, the BNIPXL BCH domain shares 17% identity (40% similarity) to the yeast Sec14p lipid-binding domain, and our extensive phylogenetic analyses indicate their likely divergence from a common ancestral domain, with many previously annotated or uncharacterized `Sec14p-like' domains in eukaryotic proteins presenting similar features to the BNIP2 family. These include their ability to form homo- and heterophilic interactions and a selectivity for small GTPases (Zhou and Low, manuscript in preparation). In this regard, recent reports demonstrate that resident Sec14p-like domains target Dbl, Ost (Ueda et al., 2004) and Dbs to Golgi membranes while mediating intramolecular interactions with the PH domain to negatively regulate Dbs transforming activity (Kostenko et al., 2005). These findings provide insight into the regulatory role of the Sec14p-like domain for GEF regulation through cis-inhibitory interactions. This mechanism is consistent with the role of the BNIPXL BCH domain in sequestering Lbc, and possibly p115RhoGEF via their DH-PH domains by interacting in trans. This is in line with previous studies demonstrating that the Lbc DH domain is essential for stress fiber formation but insufficient for efficient transformation, requiring the PH domain for targeting to punctate structures along the stress fibers (Olson et al., 1997). An alignment of the BCH domains of BNIPXL and BNIP2 with the Sec14p-like domains from Dbl, Ost, Dbs and the Sec14p lipid-binding domain indicates that they do indeed exhibit significant homology in several regions (supplementary material Fig. S8). This could pinpoint conserved motifs that execute similar regulatory functions. Therefore, it remains an exciting prospect to investigate whether other BCH or Sec14p-like domains have the capacity to target the tandem DH-PH moiety in selective members of the Dbl family GEFs.
Taken together, our findings provide a new line of evidence that BNIPXL could suppress RhoA signaling through concerted mechanisms that target both Rho and its upstream activators, all via unique motifs in its BCH domain, leading ultimately to suppression of cellular transformation. Therefore, the BCH and Sec14p-like domain families represent a novel class of regulators for small GTPases and their cognate regulators, adding greater plasticity to cellular responses mediated by these molecular switches.
Materials and Methods
Antibodies and reagents
Anti-Lbc antiserum was provided by D. Toksoz (Tufts University School of Medicine, Boston, MA). Anti-Flag, TRITC-conjugated phalloidin, goat anti-mouse and anti-rabbit horseradish peroxidase-conjugated secondary antibodies were from Sigma. Anti-c-Myc 9E10, polyclonal anti-RhoA and anti-Cdc42 antibodies were from Santa Cruz Biotechnology. Polyclonal anti-HA was from Zymed Laboratories. FITC-conjugated goat anti-mouse IgG was from Jackson ImmunoResearch Laboratories and glutathione-Sepharose beads were from Amersham Biosciences.
Expression constructs
The pGEX-RhoA, pGEX-Cdc42, pRK5-Myc-p115RhoGEF, Myc-Lbc and Myc-LbcΔC constructs were provided by A. Hall (Memorial Sloan-Kettering Cancer Center, New York, NY). pCAG-Myc-p160ROCK was provided by S. Narumiya (Kyoto University, Japan) while pGEX-Rhotekin-RBD was courtesy of S. Schoenwaelder (Monash University, Australia). pCAGGS-BMCC1-Flag was provided by A. Nakagawara (Chiba Cancer Center Research Institute, Japan). pXJ40-HA-RhoA, RhoB and RhoC was provided by Y. T. Zhou (National University of Singapore, Singapore). BNIPXL truncations were generated by incorporating XhoI (forward) and NotI (reverse) sites and BNIPXL deletion mutants were generated by site-directed mutagenesis incorporating NheI and cloned into Flag-, HA- and GFP-pXJ40 mammalian expression vectors provided by E. Manser (Institute of Molecular and Cell Biology, Singapore). pRK5-myc-Lbc and LbcΔC were subcloned into Flag- and DsRed-pXJ40. Proto-Lbc truncation mutants incorporating BamHI was generated as previously described (Sterpetti et al., 1999). Constructs were sequenced to confirm sequence fidelity.
Bioinformatics
Novel BCH-domain-containing proteins were identified using the PSI-BLAST program (Altschul et al., 1997) against the non-redundant and EST databases, as described previously (Low et al., 2000b). Multiple sequence alignments were generated by ClustalW (http://www.ebi.ac.uk/clustalw/) and formatted by BOXSHADE (http://bioweb.pasteur.fr/seqanal/interfaces/boxshade.html#fmtseq). Genome analysis was generated using the Human Genome BLAST (http://www.ncbi.nlm.nih.gov/genome/seq/BlastGen/BlastGen.cgi?taxid=9606).
Cloning and semi-quantitative RT-PCR
BNIPXL cDNA was amplified from human adult brain and kidney libraries (Clontech) using Expand Long Template PCR System (Roche Molecular Biochemicals), with 94°C for 4 minutes, 55°C for 30 seconds and 68°C for 4 minutes, followed by 30 cycles of 94°C for 40 seconds, 55°C for 30 seconds, 68°C for 4 minutes and a final extension of 68°C for 10 minutes. Primers used were full-length forward, 5′-ATGGCTTTGTTTGATGGTGATCC-3′ and reverse, 5′-CTTCTTCCAGCATGGCCAACTAAGGC-3′. For RT-PCR, 3 μg total RNA [RNeasy Mini Kit (Qiagen)] was used for oligo(dT)-primed first-strand cDNA synthesis with AMV reverse transcriptase (Promega) and further amplified with DyNAzyme (Finnzymes). Primers used were BCH forward, 5′-ATCATTGTGTTTGCCGCCTG-3′ and full-length reverse. Expression was normalized against GAPDH and verified in at least two independent experiments. The H. sapiens BNIPXLα and BNIPXLβ sequences are deposited under accession numbers AY439213 and AY439214, respectively.
Cell culture and transfection
HEK293T, MCF-7 and SK-BR3 were maintained in RPMI-1640 (Hyclone) with 10% FCS (Invitrogen), 2 mM L-glutamine, 100 U/ml penicillin and 100 μg/ml streptomycin (Hyclone). HeLa and MDA-MB-231 cells were maintained in high-glucose DMEM with 0.1 mM NEAA and 1 mM sodium pyruvate (Invitrogen). NIH3T3 cells were maintained in high-glucose DMEM with 10% bovine serum (Invitrogen) and 1.5 g/l sodium bicarbonate. All cells were grown at 37°C, 5% CO2. HEK293T and HeLa cells were transfected using FuGENE 6 (Roche Molecular Biochemicals) whereas NIH3T3 cells were transfected with Lipofectamine 2000 (Invitrogen).
Immunofluorescence
HeLa cells were processed 24 hours post transfection as previously described (Lua and Low, 2004). Confocal images were captured on Olympus FluoView FV500 microscope with a 60× oil-immersion lens. NIH3T3 cells were starved 24 hours post-transfection for 16 hours and stimulated with growth medium containing 0.5% serum for 10 minutes to induce stress fibers, according to Ridley and Hall (Ridley and Hall, 1992). Cells were fixed, stained and examined by indirect immnofluorescence. Fluorescence images were captured on Zeiss LSM510 META microscope using a 60× water-immersion lens. The percentage of transfected cells exhibiting stress fibers was evaluated for quantitative analysis of stress fiber assembly. At least 100 transfected cells per construct were evaluated for stress fiber assembly. Data are means ± s.d. from three independent experiments. Statistical significance was analyzed using ANOVA and the Student-Newman-Keuls multiple range test (Statgraphics).
Coimmunoprecipitation and western blotting
HEK293T cells were lysed (50mM HEPES, pH 7.4,150mM sodium chloride, 1.5mM magnesium chloride, 5mM EGTA, 10% glycerol, 1% Triton X-100, 5mM sodium orthovanadate, and a mixture of protease inhibitors) and directly analyzed as whole-cell lysates or aliquots used for affinity precipitation with anti-Flag M2 Affinity Gel (Sigma), as previously described (Low et al., 2000a). Bound proteins were eluted and separated by SDS-PAGE for western analyses. For partial endogenous immunoprecipitation, cells were lysed in NP-40 lysis buffer (50 mM Tris-HCl pH 8.0, 150 mM sodium chloride, 1% Nonidet P-40 and a mixture of protease inhibitors). 5 mg total protein was subjected to immunoprecipitation with anti-Flag M2 Affinity Gel overnight at 4°C and processed as described (Park et al., 2002).
Direct protein-binding assays
GST fusion proteins were purified from E. coli BL21(DES) (Shang et al., 2003). In-vitro-transcribed and translated Flag-tagged proteins were synthesized [TNT Quick Coupled Transcription/Translation System (Promega)], and aliquots incubated with GST-RhoA (10 μg) overnight at 4°C. Samples were diluted with 100 μl cell lysis buffer and incubated for 1 hour at 4°C. Beads were washed three times with cell lysis buffer and bound proteins separated by SDS-PAGE for western analyses.
In vitro GTP and GDP loading assays
GST-RhoA (10 μg) were preloaded with 10 mM GDP or GTPγS (Sigma) in binding buffer (25 mM NaCl, 20 mM Tris-HCl pH 7.5, 0.1 mM DTT) and 50 mM EDTA at 30°C for 10 minutes. The reaction was stopped with MgCl2 to a final concentration of 20 mM (Ren and Schwartz, 2000). The beads were incubated with lysates at 4°C in GAP lysis buffer (50 mM HEPES pH 7.4, 150 mM NaCl, 1.5 mM MgCl2, 5 mM EGTA, 10% glycerol, 1% Triton X-100, 5 mM sodium orthovanadate and a mixture of protease inhibitors) and washed with immunoprecipitation wash buffer (20 mM HEPES pH 7.4, 150 mM NaCl, 10 mM MgCl2, 10% glycerol, 0.1% Triton X-100). Bound proteins were separated by SDS-PAGE for western analyses.
RhoA activation assays
HEK293T cells were lysed in lysis buffer (50 mM Tris-HCl pH 7.4, 2 mM magnesium chloride, 1% Triton X-100, 10% glycerol, 100 mM sodium chloride, supplemented with 1 mM DTT and a mixture of protease inhibitors) and 300 μg lysates was incubated with GST-Rhotekin-RBD for 30 minutes at 4°C (Besson et al., 2004). The beads were washed three times with lysis buffer and bound proteins separated by SDS-PAGE for western blot analyses.
Short-hairpin RNA (shRNA) knockdown
Six shRNA oligonucleotides against human BNIPXL were designed using the siRNA target finder (Ambion), synthesized (Operon) and cloned into the pSilencer 2.1-U6 hygro vector (Ambion). The shRNA sequences are sh473, 5′-GATCCCATGCTTTCACCACAGCCTGTTCAAGAGACAGGCTGTGGTGAAAGCATTTTTTTGGAAA-3′; sh477, 5′-GATCCCGTGACACTCATCCTCGGAGATTCAAGAGATCTCCGAGGATGAGTGTCATTTTTTGGAAA-3′; sh479, 5′-GATCCCACAATAGAGGAAGCTGGGTTTCAAGAGAACCCAGCTTCCTCTATTGTTTTTTTGGAAA-3′; sh47C, 5′-GATCCCGAGCAGCGCATTGACATGATTCAAGAGATCATGTCAATGCGCTGCTCTTTTTTGGAAA-3′; sh470, 5′-GATCCGCGCATTGACATGAAGGTCATCTCGAGATGACCTTCATGTCAATGCGCTTTTTTGGAAA-3′; sh484, 5′-GATCCGCCAACAAAGATTCTGGCCAACTCGAGTTGGCCAGAATCTTTGTTGGCTTTTTTGGAAA-3′.Constructs were sequenced to confirm sequence fidelity and knockdown efficiency was determined by cotransfection with GFP- or HA-BNIPXL plasmid at a 10:1 ratio. Lysates were analyzed 72 hours post transfection.
Transformation assays
Primary focus formation assays were performed in NIH3T3 cells using calcium phosphate precipitation (Solski et al., 2000). NIH3T3 cells seeded at 2.5×105 cells/60 mm dish were cotransfected with DsRed-LbcΔC and GFP-BNIPXL, GFP-WBCH or GFP-CBCH and grown in complete medium. Cells were fixed (10% acetic acid, 10% methanol), stained with 0.5% Crystal Violet and foci formation was scored 21 days post transfection. Each experiment was performed at least three times. Data were analyzed for statistical significance using ANOVA and the Student-Newman-Keuls multiple range test (Statgraphics).
Acknowledgements
We thank A. Hall (Memorial Sloan-Kettering Cancer Center, USA), S. Narumiya (Kyoto University, Japan), S. Schoenwaelder (Monash University, Australia), D. Toksoz (Tufts University School of Medicine, USA) and A. Nakagawara (Chiba Cancer Center Research Institute, Japan) for generous gifts of plasmid constructs and antisera. We thank the National University Medical Institutes for the use of confocal microscopy units, Chew Li Li for helpful discussions and Aarthi Ravichandran for critical review of our manuscript. This work is supported in part by a grant from Biomedical Research Council of Singapore and the BMRC Young Investigator Award to B.C.L.