|Home | About | Journals | Submit | Contact Us | Français|
BECN1/Beclin 1 has a critical role in the early stages of autophagosome formation. Recently, structures of its central and C-terminal domains were reported, however, little structural information is available on the N-terminal domain, comprising a third of the protein. This lack of structural information largely stems from the inability to produce this region in a purified form. Here, we describe the expression and purification of the N-terminal domain of BECN1 (residues 1 to 150) and detailed biophysical characterization, including NMR spectroscopy. Combined, our studies demonstrated at the atomic level that the BECN1 N-terminal domain is intrinsically disordered, and apart from the BH3 subdomain, remains disordered following interaction with a binding partner, BCL2L1/BCL-XL. In addition, the BH3 domain α-helix induced upon interaction with BCL2L1 reverts to a disordered state when the complex is dissociated by exposure to a competitive inhibitor. No significant interactions between N- and C-terminal domains were detected.
Macroautophagy (hereafter referred to as autophagy) is a highly conserved catabolic process involving the formation of autophagosomes.1,2 These vesicles encapsulate cellular components (such as macromolecules and organelles) and pathogens, which are delivered to lysosomes for degradation. BECN1, a haploinsufficient tumor suppressor, is the first identified mammalian autophagy gene.3-5 It plays a central role in the autophagy initiation process involving nucleation of the autophagic vesicle through interactions with other proteins including PIK3R4/Vps15, ATG14/Barkor, and the class III phosphatidylinositol 3-kinase (PtdIns3K; PIK3C3/Vps34)6 which all assemble into a PtdIns3K complex.
Structurally, BECN1 is a 450 amino-acid protein with 3 well-characterized domains (Fig. 1A). Residues 112 to 123 comprise a BH3 domain, which is a short conserved sequence motif important for mediating interactions between members of the BCL2 family in the intrinsic apoptotic pathway. BECN1 was actually discovered as a BCL2-interacting protein3,7 and interactions between the BECN1(BH3) domain and BCL2 and BCL2L1 are well-characterized both biochemically and structurally,8–16 though their physiological relevance has been recently challenged.17,18 A potential helical domain has been proposed (based on unpublished data) between residues 141 to 170,8 which is followed by a coiled-coil region (residues 174 to 266) that mediates BECN1 homodimerization as well as heterodimerization with UVRAG/Vps38 and ATG14.16 A crystal structure of this domain shows it forms an antiparallel dimer characterized by an imperfect interface.19 The C-terminal region of BECN1 after the coiled-coiled domain is called the β-α repeat autophagic-specific (BARA) domain,20,21 which contains a highly conserved region known as the evolutionarily conserved domain (ECD) (residues 244 to 337). The BARA domain is critical for BECN1 activity as it mediates interactions with membrane lipids. A crystal structure of the BARA domain of yeast and human BECN1 orthologs has been solved and showed it to represent a novel fold and a new class of membrane-binding domains.20,22
A low resolution (28 Å) structure, determined by single particle electron microscopy, of one of the entire human PtdIns3K complexes (there are at least 3 such complexes that contain common as well as unique components) has been reported.23 This showed BECN1 has an extended structure that makes contact with ATG14 and PIK3R4. The model also suggested the BECN1 and ATG14 coiled-coil domains are parallel unlike the antiparallel arrangement in homodimeric BECN1.19 Although the BECN1 N terminus was mapped to the pivot point of the V-shaped PtdIns3K complex in the EM structure, no density was visible for this region. Very recently, a crystal structure of the yeast PtdIns3K complex II vacuolar sorting complex (PtdIns3K complex I functions in macroautophagy), that also involves the BECN1 ortholog (Vps30/Atg6), has been determined.24 This complex has a similar overall architecture to the mammalian PtdIns3K autophagy complex, although Atg14 is replaced with Vps38/UVRAG. The Vps30/Atg6 structure shows 2 adjacent coiled-coil domains (a short CC1 and longer CC2) with several contacts between both of these with PIK3R4 and PIK3C3/Vps34. Interestingly, while the Vps30/Atg6 N-terminal domain was mostly undefined in this structure, there was some evidence that this region could potentially contain 3 disconnected helices that interact with Vps38/UVRAG.
Prior to this crystal structure of Vps30/Atg6 in the context of the PtdIns3K complex II, there was no detailed structural information on the entire N-terminal region of BECN1 beyond studies on the ~20 residue BH3 domain, which is unstructured but becomes helical upon engagement of BCL2 prosurvival proteins.8,10-13 As this region comprises approximately one-third of the protein (150 amino acid residues), it is of fundamental interest to characterize it structurally. Recently, extensive primary sequence analysis using multiple algorithms to search for intrinsically disordered regions in autophagy regulators indicated that a large proportion of the BECN1 N terminus (residues 42 to 115) could be intrinsically disordered (consistent with the absence of density in the EM and crystallography studies), although subsequent analyses using NMR spectroscopy could only focus on the BH3 region due to difficulties associated with producing the entire domain.8 An earlier study using NMR to examine a construct that encompasses residues 101 to 154 concluded it was structured,25 although the fragment they purified and analyzed was associated with a “solubility tag” (the 8-kDa streptococcal protein G B1 domain) which itself is a fully folded protein.26,27 Hence, the interpretation of these data must be carefully considered. Intrinsically disordered regions are functionally very important in proteins because they are often the sites of protein-protein interactions and/or post-translational modifications.28
In this paper we report the expression, purification and structural characterization of the entire N-terminal domain of BECN1. We show that it is essentially entirely disordered, and with the exception of the BH3 subdomain, remains unstructured following complex formation with the BCL2 family member BCL2L1. Our results also showed, for the first time, the helical conformation adopted by the BH3 domain of BECN1 upon binding BCL2L1 is reversed to an intrinsically disordered state in solution when the complex is dissociated.
Initial attempts to produce the BECN1 N-terminal domain (residues 1 to 150) were hindered by its sensitivity to bacterial proteases. This was largely overcome by incorporation of a cocktail of protease inhibitors in the purification protocol, and ensuring all procedures were performed at low temperature where possible. As a result of these changes, we were generally able to purify 2 to 4 mg of the various BECN1 constructs per liter of culture media (both minimal and complex). Upon gel-filtration chromatography, the major protein-containing peak eluted at an earlier retention time than expected for a protein of the calculated molecular mass based on the amino acid sequence (Fig. 1B). This could indicate that the protein eluted as some type of oligomer (e.g., a dimer or trimer based on the molecular weight standards). Alternatively, aberrant behavior on gel-filtration chromatography can indicate that the protein is intrinsically disordered, as these proteins have a large hydrodynamic radius compared to globular proteins.29 To ascertain the exact oligomeric state of the BECN1(1–150)[4CS]Y protein (see Fig. 1A and Materials and Methods for details of this and other constructs) we examined its sedimentation velocity by analytical ultracentrifugation. These data demonstrated that the protein exists as a 1.2S monomer (Fig. 1C), hence the behavior on gel-filtration is indicative of the protein being intrinsically disordered rather than oligomeric. The highly elongated shape of the molecule, as determined by its high frictional ratio in the sedimentation velocity experiment (f/f0: 2.2), is also a characteristic of an intrinsically disordered protein.30 Characterization by Coomassie-stained SDS-PAGE indicated the protein in the major gel filtration peak fractions was highly pure (>90%) with minimal proteolysis (Fig. 1B, inset). The slightly anomalous migration on SDS-PAGE is likely due to the high positive charge content of the molecule, another characteristic of an intrinsically disordered protein.31
It is important to note that the purified protein tended to form visible aggregates during attempts to concentrate it following gel-filtration chromatography, therefore subsequent NMR and other analyses were generally performed on the unmanipulated peak fraction (which was typically at ~150 μM) obtained in this purification step. Mutation of the 4 cysteine residues to serine reduced the amount of aggregated material we observed in the gel-filtration step (present in the minor peaks that eluted between 50 to 80 mls) and improved long-term storage without significantly influencing the NMR spectrum (data not shown).
To check the “activity” of the purified BECN1 N-terminal domain (BECN1(1–150)[4CS]Y) we determined its relative affinity for BCL2L1 (Fig. 1D). In solution competition assays, the IC50 (2.0 μM) value of BECN1(1–150)[4CS]Y was similar to that of a synthetic peptide whose sequence overlaps with the BECN1(BH3) domain (IC50 2.9 μM) and consistent with previously published binding data on BECN1(BH3) domain peptides.10,12,16 The similarity in the values between the BECN1(BH3) peptide and the entire N-terminal domain suggests the protein produced is fully functional and that regions outside of the BH3 domain are not involved in interactions with BCL2 family members.
Although various characteristics of the BECN1 N terminus during purification suggested it was intrinsically disordered (sensitivity to proteases, anomalous behavior on gel-filtration and SDS-PAGE), we next performed circular dichroism (CD) analysis to establish whether there was any significant secondary structure present. The CD spectrum showed a single minimum at approximately 200 nm (Fig. 2A) that is indicative of random or disordered secondary structure. Calculation of α-helix content from these data indicated none was present, demonstrating that the BH3 domain in its native context (rather than in short peptides used in previous studies) must only assume the helical conformation observed in complexes with BCL2L111,12 upon binding to the protein.
As all of the data described above provided strong evidence for the BECN1 N terminus being intrinsically disordered, we next performed NMR studies to gain residue-specific structural information. The 1H-15N HSQC spectrum of BECN1(1–150)[4CS]Y (Fig. 2B) displayed very narrow chemical shift dispersion (between 7.7 and 8.7 ppm) on the 1H dimension that is a characteristic feature of intrinsically disordered proteins.32,33 Analyses of triple resonance 3-dimensional spectra of BECN1(1–150)[4CS]Y together with spectra of BECN1(1-76)[2CS]Y (Fig. 1A), which helped to resolve a number of ambiguous resonances due to narrow spectral dispersion of BECN1(1–150)[4CS]Y, led to complete assignments of backbone 13Cα, 13Cβ, 13C, 15N, 1HN, and 1Hα resonances for all non-proline residues, and 13C (13Cα, 13Cβ and 13C) resonances for the prolines except P84. These assignments have been deposited in BioMagResBank (http://www.bmrb.wisc.edu, Accession code: 25384). Experimentally observed backbone chemical shifts of BECN1(1–150)[4CS]Y were then used to evaluate its secondary structure with the improved version of Chemical Shift Index software (CSI 2.0)34 which makes use of all 6 backbone chemical shifts (13Cα, 13Cβ, 13C, 15N, 1HN, and 1Hα), along with sequence-derived features. This analysis confirmed that BECN1 is disordered throughout its entire N terminus (1 to 150) except for a potential half-helical turn involving residues A44 and P45 (Fig. 2C). Furthermore, the nearly identical chemical shifts observed for resonances arising from BECN1(1–76)[2CS]Y and BECN1(1–150)[4CS]Y (apart from those arising from the C-terminal truncation site of BECN1(1–76)[2CS]Y) indicates that it is unlikely there are any direct contacts between BECN1 residues 1 to 76 and BECN1 residues 77 to 150.
Although there was some evidence of electron density associated with the N-terminal domain of the yeast BECN1 ortholog, Vps30/Atg6, in the crystal structure of the yeast PtdIns3K complex II, structural studies on the mammalian BECN1 N terminus have been restricted to the BH3 domain in complex with BCL2L1. However, as we had obtained a fully assigned 1H-15N HSQC spectrum of the N-terminal 150 residues, we were now able to assess whether there were any global structural changes along the whole N-terminal domain upon interaction with BCL2L1. When BCL2L1ΔC25 was titrated into the BECN1(1–150)[4CS]Y sample, only resonances located within or adjacent to the BH3 domain (residues 98 to 132) were significantly perturbed (Fig. 3A to H) with only a small number of new resonances detected in the BCL2L1ΔC25-bound form (highlighted in cyan circles in Fig. 3A). This indicates the exchange between the free and bound conformations of BECN1(1–150)[4CS]Y is intermediate (in the millisecond range) on the NMR chemical shift timescale. Indeed, protein-protein interactions involving intrinsically disordered proteins in solution, as studied by NMR, are frequently hampered by excessive line broadening due to the exchange between free and receptor-bound conformations occurring on this timescale which precludes direct detection of the bound-state resonances.35 Interestingly, the large majority of resonances of BECN1(1–150)[4CS]Y were completely unaffected (both peak intensity and linewidth) by the binding of its BH3 domain with BCL2L1ΔC25 (except for a small loss in intensity in residues adjacent to Q40, Fig. 3B), indicating that nearly all residues outside of the BH3 domain are tumbling as freely when bound to BCL2L1 as in the unbound form, and do not adopt any secondary or tertiary structure elements that are not present in the free form.
Although the conformation of BECN1(1–150)[4CS]Y in complex with BCL2L1 under the current conditions cannot be quantified due to the absence of resonances in the BH3 region, our NMR data combined with published structures of BCL2L1-BECN1(BH3) complexes,11,12 enable us to conclude that residues between 98 to 132 that engage BCL2L1 likely adopt a helical conformation, while the rest of BECN1(1–150)[4CS]Y remains unstructured. To establish whether this secondary structure transition is reversible, we next acquired 1H-15N HSQC spectra following titration of a BCL2L11(BH3)/BIM(BH3) synthetic peptide into the BCL2L1ΔC25-BECN1(1–150)[4CS]Y complex. BCL2L11(BH3) binds BCL2L1 with ~1000-fold higher affinity than BECN1(BH3), hence would be expected to outcompete the associated BECN1 N-terminal domain. Indeed, after addition of the BCL2L11(BH3) peptide, all of the resonances of BECN1(1–150)[4CS]Y residues that were disturbed by binding to BCL2L1ΔC25 (i.e., those within or adjacent to the BH3 domain) began to re-emerge at their original locations (Fig. 3I to N). This is most obvious by examining just the characteristic glycine region of the 1H-15N HSQC spectra (Fig. 3I-M). Moreover, when BCL2L11(BH3) was in molar excess to BCL2L1, the 1H-15N HSQC spectrum of BECN1(1–150)[4CS]Y completely reverted to that seen in the absence of BCL2L1 suggesting that any structural change that occurred was fully reversible (Fig. 3M, N).
Our fully assigned spectrum for the BECN1 N terminus also enabled us to assess whether part(s) of it engages in any significant interactions with other regions on the BECN1 protein, at least when present in cis. Accordingly, a 1H-15N HSQC spectrum for 15N-labeled BECN1(1–150)[4CS]Y was recorded in the presence of an equimolar concentration of purified BECN1(BARA/ECD) domain (residues 248 to 450) produced recombinantly in E. coli (Fig. 4A). When compared to a spectrum of BECN1(1–150)[4CS]Y recorded in the presence of an equal volume of buffer (to account for dilution effects), we observed only very minor changes in peak intensities along the N terminus of BECN1 (e.g., residues 96 to 101, 115, 116, 121, 122, 128, 145) (Fig. 4B, C), unlike the very large changes observed for residues within or adjacent to the BH3 domain upon titration of BCL2L1 (Fig. 3B). Hence, while these data suggest there is some potential for regions within the BECN1 N-terminal region to interact with its C-terminal domain, these interactions are likely very weak and/or transient, and probably not critical for the overall structure of full-length BECN1.
Several reports have demonstrated that deletion of a large portion (residues 88 to 144) of the BECN1 N-terminal domain or indeed the entire N terminus of yeast Vps30/Atg6, does not affect its autophagic activity.7,20,36 Similar deletion studies have also shown the N-terminal domain does not have autophagic activity on its own.37 However, as it contains a BH3 domain, which is the critical functional domain in proapoptotic BH3-only proteins, it was of interest to establish whether the BECN1(BH3) sequence possesses any apoptotic activity. As such we generated a chimeric construct in which the BH3 domain of the proapoptotic BH3-only protein BCL2L11 (BCL2L11Sisoform) was replaced with the corresponding BECN1(BH3) sequence, or those from other proapoptotic BH3-only proteins (BAD/BCL2L8 and PMAIP1/NOXA) as specificity controls. This is a well-established approach for comparing the killing activity of different BH3 domains in vitro and in vivo.38–40 Killing activity was assessed by long-term (7 d) colony formation assays in mouse embryonic fibroblasts (MEFs) of different genotypes. MEFs are dependent on BCL2L1 and MCL1/BCL2L3 for their survival,39 hence both of these proteins must be neutralized for cell-killing activity. BCL2L11S, which can neutralize all BCL2 prosurvival proteins,40 killed all cell lines except those deficient in the critical mediators of apoptosis, BAX and BAK1, indicating its mechanism-based activity. In contrast, the BCL2L11S-BAD(BH3) chimera which targets BCL2L1, BCL2 and BCL2L2/BCLW only killed mcl1 knockout MEFs, while the BCL2L11S-PMAIP1(BH3) chimera (which only targets MCL1) only killed bcl2l1/bclx knockout MEFs, consistent with the specificity profiles of both of these constructs.40 Unlike all of these other constructs, no killing activity was observed with the BCL2L11S-BECN1(BH3) chimera in any cell line tested suggesting its affinity for BCL2 proteins, especially MCL1 and BCL2L1, is too low to functionally neutralize them in cells.
In this paper we provide the first report of the expression and purification of the BECN1 N-terminal domain (residues 1 to 150). This was the critical first step toward the structural evaluation of this previously mostly uncharacterized region that comprises a third of the entire BECN1 protein. Our results from CD and NMR analysis showed that this entire region (residues 1 to 150) is intrinsically disordered, lacking any secondary or tertiary structure. This conclusion was supported by a range of other findings including its susceptibility to proteolysis, aberrant behavior on SDS-PAGE and gel-filtration chromatography, and the high frictional ratio (f/f0) upon analytical centrifugation. This finding is consistent with recent sequence analysis studies that used various algorithms to come to a similar conclusion, though the “consensus” disordered region derived from the results of 4 different prediction programs suggested only residues 42 to 115 were unstructured. It is also consistent with the data on the recent EM structure reported for the PtdIns3K complex where no density was apparent for the BECN1 N terminus (indicating it is disordered),23 and provides an explanation for the difficulties in obtaining crystals of full-length BECN1 in the absence of stabilizing binding partners.
The N-terminal domain of Vps30/Atg6 in the crystal structure of the yeast PtdIns3K vacuolar sorting complex (complex II) was also not well defined, consistent with it being disordered, though there was some electron density that was assigned to this region which appears to form 3 disconnected helices involved in interactions with Vps38 (UVRAG).24 Unfortunately, the identity of these potentially structured residues was not discernable due to the low resolution of the structure. Nevertheless, if these residues can be definitively assigned to the N terminus, then this result is interesting as it suggests that regions outside of the BH3 domain can also undergo disordered-to-ordered structural transitions upon interactions with other proteins. Such structural transitions are often observed in intrinsically disordered proteins or domains.41 It should also be noted, however, that the yeast N-terminal domain has a long insertion (~55 residues) relative to the corresponding domain in mammalian BECN1 so there could be some difference(s) in how this region behaves in different species. It is also unclear from the available data whether the same structural transition occurs in the context of the PtdIns3K autophagy complex.
We were not able to detect any significant interactions between the BECN1 N terminus and C-terminal BARA/ECD domains when examined in cis (Fig. 4), suggesting these domains do not interact in the full-length BECN1 protein. This agrees with the PtdIns3K complex EM structure and PtdIns3K complex II crystal structure in which BECN1/Vps30 appears highly elongated with the N and C termini being located distal to each other, separated by the coiled-coil domains. Similarly, our NMR studies on the truncated BECN1 residues 1 to 76 N-terminal fragment indicated there are no significant intradomain interactions between the N- and C-terminal halves of the N-terminal domain. More importantly, our NMR analysis showed that upon binding to BCL2L1, the whole of the N-terminal region except the BH3 domain remains unstructured, suggesting that this interaction does not “nucleate” additional structural transitions beyond the formation of the α-helical BH3 domain.8,10–12 In addition, we showed that the BH3 domain reassumes its disordered conformation once it dissociates from its binding partner. Although the physiological relevance of interactions between BECN1 and mammalian prosurvival members of the BCL2 family has recently been questioned,17,18 interactions between BECN1 and viral BCL2 orthologs have also been shown to be mediated by the BECN1(BH3) domain.3,7,13,42 Hence, such interactions could be important for virus survival through inhibition of host cell autophagy following infection. As such, structural studies on BECN1 following such interactions, such as described here, will be important for fully understanding BECN1 biology.
In the past 2 decades the extent and importance of structural disorder in proteins has become increasingly apparent and is now recognized as being widespread and a key factor driving protein-protein interactions.31,43-46 In the BCL2 family, the majority of the BH3-only subclass of proteins, critical for triggering the apoptotic cascade, are intrinsically unstructured with only the BH3 domain becoming α-helical after binding the prosurvival family members.47,48 In this context, our data suggests the BECN1 N terminus fits the BH3-only protein structure archetype perfectly. However, in contrast to all BCL2 family BH3-only proteins, the affinity of the BECN1(BH3) domain for prosurvival protein targets is relatively weak, providing an explanation for the lack of cell-killing we showed here in mammalian cell survival assays (Fig. 5). Other studies using overexpressed BECN1 have also shown it to be lacking apoptosis-inducing activity despite possessing a BH3 domain.15 This makes sense because, although a BH3 domain with high affinity for BCL2 prosurvival proteins would better enable autophagy-inhibitory interactions, it would also be more capable of inducing cell death, which would not be desirable. Interestingly the affinity of the BECN1(BH3) domain is significantly higher for some viral BCL2 proteins compared to its mammalian orthologs, suggesting a critical role for the inhibition of autophagy in some viral infections.14,42
Greater than 50% of autophagy effector proteins are predicted to possess an intrinsically disordered region.8 The role of these unstructured regions in autophagy however is still poorly understood, including the function of the N-terminal domain of BECN1 (beyond the BH3 domain). However, intrinsically disordered regions in proteins can be functionally very important as their lack of structure often enables them to engage in protein-protein interactions with multiple binding partners, including those that can create post-translational modifications such as kinases and ubiquitin ligases.28,49 Indeed, the BECN1 N terminus contains a ubiquitination site50 and numerous phosphorylation sites51–60 (Fig. 6) including several within the BH3 domain that can both positively and negatively regulate autophagy. Hence, although the BECN1 N-terminal domain might be dispensable for its autophagic activity,7,20,36 our finding that the BECN1 N-terminal domain is intrinsically disordered is entirely consistent with its critical role as a regulatory hub that influences autophagy through phosphorylation, ubiquitination and other protein:protein interactions.
The human BECN1 DNA encoding residues 1 to 150 was codon-optimized for expression in E. coli and cloned into the pET DUET-1 vector (Site I) (Novagen, 71146–3) at the BamHI/NotI sites, which fuses it to a hexahistdidine purification tag followed by a short linker sequence (MGSSHHHHHHSQDP). This served as a base construct for a subsequent construct in which all cysteines (i.e., residues 18, 21, 137, 140) were mutated to serine plus a C-terminal tyrosine was added to facilitate quantification by UV spectroscopy (called “BECN1(1–150)[4CS]Y”). An additional truncated version of this construct (residues 1 to 76, BECN1(1–76)[2CS]Y), also with a C-terminal tyrosine, was made to facilitate NMR spectrum assignment. The BECN1(BARA/ECD) domain (encoding residues 248 to 450) construct was cloned identically into the pET DUET-1 vector.
To produce BCL2L1 we used a C-terminally truncated construct (BCL2L1ΔC25), which removes the hydrophobic C-terminal transmembrane region that causes the protein to aggregate. This construct was expressed as described previously.40 BECN1 constructs were transformed and expressed in the BL21 DE3 E. coli strain. Cultures were grown at 37°C in Super Broth or minimal medium61 for 15N or 15N/13C labeling to OD600nm 0.5–1.0 before expression was induced by addition of IPTG to a final concentration of 500 µM. Expression was induced for 3 h at 37°C. Cell pellets were stored at −80°C until proteins were purified.
BCL2L1ΔC25 was purified as described previously.40 All BECN1 protein constructs including N-terminal and BARA/ECD domains were histidine-tagged and purified using the same purification steps which were all conducted at 4°C with prechilled solutions. The cell pellet was resuspended in Tris-buffered saline (TBS; 20 mM Tris, pH 8.0, 150 mM NaCl) using ~40 mL for the pellet from a 2 L culture, containing 250 units of DNase 1 (Roche, 04536282001) and 3 tablets of protease inhibitor cocktail (Roche, 05056489001). Cells were lysed by 2 passages through a pressure cell homogenizer (Stansted, model SPCH-10). After centrifugation for 30 min at 48,000g the supernatant fraction was filtered through a 0.2-μm filter prior to loading samples onto a 1 mL HiTrap Chelating HP column (GE Healthcare, 17–0408–01) precharged with nickel sulfate solution. The column was then washed with 20 column volumes TBS containing 20 mM imidazole and the bound protein eluted with 5 column volumes TBS containing 250 mM imidazole. The eluted protein was then injected directly onto a Superdex S75 16/60 column equilibrated with 50 mM phosphate buffer, pH 6.8 containing 50 mM NaCl, at 1 mL/min. Peak fractions (based on absorbance 280 nm for tyrosine-containing constructs or 214 nm for native construct) were analyzed on NuPAGE 4–12% Bis-Tris gels (Invitrogen, NPO335BOX) stained with Coomassie Brilliant Blue R-250. The concentration of the eluted peak protein fractions was determined using the Beer-Lambert equation from the absorption at 280 nm determined using a Nano drop spectrophotometer. The purest fraction(s) were then used directly for subsequent analysis without further manipulation. For long-term storage of NMR samples, NaN3 was added to a final concentration of 0.02 % (w/v).
Peptides were synthesized by Mimotopes (Clayton, Victoria, Australia) and purified by reversed-phase HPLC to >90 % purity. The sequences were: BECN1(BH3) - GEASDGGTMENLSRRLKVTGDLFDIMSGQTDVDH; BCL2L11(BH3)/BIM(BH3) –DMRPEIWIAQELRRIGDEFNAYYARR, used previously for NMR studies.62
Solution competition assays were performed using a Biacore 3000 instrument (GE Healthcare, Little Chalfont, UK) as described previously.63 Briefly, recombinant BCL2L1ΔC25 (5 nM) was incubated with the BECN1(BH3) synthetic peptide or BECN1(1–150)[4CS]Y for at least 2 h prior to injection onto a CM5 sensor chip on which either wild-type BCL2L11(BH3) peptide or an inert BCL2L11(BH3) mutant peptide (BCL2L114E) was immobilized. Specific binding of BCL2L1 to the surface in the presence and absence of BECN1 competitors was quantified by subtracting the signal obtained on the BCL2L11(BH3)-mutant channel from that obtained on the wild-type-BCL2L11(BH3) channel. The ability of the BECN1(BH3) synthetic peptide or BECN1(1–150)[4CS]Y to inhibit BCL2L1 binding to immobilized BCL2L11(BH3) was expressed as IC50, calculated by nonlinear curve fitting of the data by using GraphPad Prism (GraphPad Software)
Sedimentation velocity experiments were conducted in a Beckman model XL-A analytical ultracentrifuge at a temperature of 20°C. BECN1(1–150)[4CS]Y (50 µM in 50 mM sodium phosphate, 50 mM sodium chloride, pH 6.8) was loaded into a conventional double sector quartz cell and mounted in a Beckman 4-hole An-60 Ti rotor (Beckman Coulter Indianapolis, IN, USA). Both the sample (380 µL) and buffer solution (400 µL) were centrifuged at a rotor speed of 116,297 x g, and the data were collected at a single wavelength (236 nm) in continuous mode, using a time interval of 0 s and a step-size of 0.003 cm without averaging. Solvent density (1.0061 g/mL at 20°C) and viscosity (1.0287 cp), and an estimate of the partial specific volume (0.7153 mL/g) were computed using the amino-acid composition of the BECN1(1–150[4CS]Y peptide and the program SEDNTERP.64 Sedimentation velocity data at multiple time points were fitted to continuous size-distribution models65–67 using the program SEDFIT, which is available at www.analyticalultracentrifugation.com.
Circular dichroism spectra were recorded at 20°C using an AVIV Model 420 CD spectrometer (Aviv Biomedical, NJ, USA). Wavelength scans were performed between 195 and 250 nm in aqueous buffer (50 mM sodium phosphate, 50 mM NaCl, pH 6.8) using 0.15 mg/mL of protein in a 1 mm quartz cuvette. The resulting raw data were converted to mean residue elipticity ([θ]) and the α-helical content estimated from ([θ]222nm) as described previously.65
15N-labeled (~0.18 mM) and 13C/15N labeled (~0.12 mM) BECN1(1–150)[4CS]Y were prepared in 50 mM sodium phosphate containing 50 mM sodium chloride and 0.02 % (w/v) sodium azide at pH 6.8. A 15N-labeled N-terminal half of BECN1(1–150)[4CS]Y, BECN1–76[2CS]Y (~0.30 mM), was also prepared in the same buffer solution for the present study. Two-dimensional 1H-15N HSQC spectra of BECN1(1–150)[4CS]Y over the temperature range of 5 to 30°C were first acquired and evaluated. Spectra at 15°C provided a good compromise between molecular rotational tumbling and relatively rapid backbone amide exchange of BECN1(1–150)[4CS]Y. All spectra used for backbone sequential assignments were subsequently recorded at 15°C on a Bruker Avance800 spectrometer (Bruker BioSpin GmbH, Rheinstetten, Germany) equipped with a cryoprobe. Backbone assignments were achieved from standard triple resonance HNCA, HN(CO)CA, HNCACB, CBCA(CO)NH, HNCO and HN(CA)CO spectra68 whereas backbone Hα chemical shifts were obtained from HNHA and 15N-HSQC-TOCSY experiments. NMR data were processed in Topspin (Version 3.2, Bruker) with the 1H chemical shifts referenced indirectly to DSS at 0 ppm via the H2O signal (4.885 ppm at 15°C), and the 13C and 15N chemical shifts were referenced indirectly using absolute frequency ratios.69 Spectra analyses were carried out in XEASY Version 1.3.70 Experimentally determined backbone chemical shifts (13Cα, 13Cβ, 13C, 15N, 1HN, and 1Hα) of BECN1(1–150)[4CS]Y were then used for the identification of secondary structure using the program CSI 2.0.34
To validate the interaction between BECN1(1–150)[4CS]Y and BCL2L1ΔC25, BCL2L1ΔC25 (0.46 mM) was titrated into the 13C/15N BECN1(1–150)[4CS]Y solution (0.145 mM) with the course of interaction monitored by 1H-15N HSQC spectra. After a molar excess of BCL2L1ΔC25 in solution was reached (i.e., BECN1(1–150)[4CS]Y became fully bound), the BCL2L11(BH3) peptide (2 mM) was subsequently titrated into the solution containing the (BECN1(1–150)[4CS]Y)-BCL2L1ΔC25 complex until a molar excess of BCL2L11(BH3) (versus BCL2L1ΔC25) was reached (no bound BECN1(1–150)[4CS]Y present in solution). For studies looking at the interactions between the N-terminal and BARA/ECD domains, 15N BECN1(1–150)[4CS]Y solution (0.150 mM) was mixed in a 1:1 molar ratio with unlabelled BARA/ECD domain or equivalent volume of buffer and the interaction monitored by the 1H-15N HSQC spectra.
Long-term clonogenic assays were performed exactly as described recently.71 Briefly MEFs (wild-type, bax−/− bak1−/−, bcl2l1−/−, mcl1−/−) were infected with retroviruses to express BCL2L11S, BCL2L11S-BECN1(BH3), BCL2L11S-BAD(BH3) or BCL2L11S-PMAIP1(BH3) where their expression was linked to that of green fluorescent protein via an internal ribosome entry site. Infected (green fluorescent protein+) cells were sorted and plated, then colonies counted 7 d later. The BCL2L11S-BECN1(BH3) construct has the BCL2L11(BH3) domain (DMRPEIWIAQELRRIGDEFNAYYARR) replaced with the analogous sequence (DGGTMENLSRRLKVTGDLFDIMSGQT) in BECN1. BCL2L11S-BAD(BH3) and BCL2L11S-PMAIP1(BH3) have been described previously.40
No potential conflicts of interest were disclosed.
We thank Dr Jeff Babon (WEHI) for helpful discussions.
This work was supported by grants and fellowships from the NHMRC of Australia (Project Grant 1049949 and Career Development Fellowship 1024620 to E.F.L). Infrastructure support from NHMRC IRIISS grant #361646 and the Victorian State Government OIS grant is gratefully acknowledged.