Key steps in unconventional secretion of fibroblast growth factor 2 reconstituted with purified components

  1. Julia P Steringer
  2. Sascha Lange
  3. Sabína Čujová
  4. Radek Šachl
  5. Chetan Poojari
  6. Fabio Lolicato
  7. Oliver Beutel
  8. Hans-Michael Müller
  9. Sebastian Unger
  10. Ünal Coskun
  11. Alf Honigmann
  12. Ilpo Vattulainen
  13. Martin Hof
  14. Christian Freund
  15. Walter Nickel  Is a corresponding author
  1. Heidelberg University Biochemistry Center, Germany
  2. Freie Universität Berlin, Germany
  3. Academy of Sciences of the Czech Republic, Czech Republic
  4. University of Helsinki, Finland
  5. Tampere University of Technology, Finland
  6. Max Planck Institute of Molecular Cell Biology and Genetics, Germany
  7. Paul Langerhans Institute Dresden, Germany
  8. University of Southern Denmark, United Kingdom

Abstract

FGF2 is secreted from cells by an unconventional secretory pathway. This process is mediated by direct translocation across the plasma membrane. Here, we define the minimal molecular machinery required for FGF2 membrane translocation in a fully reconstituted inside-out vesicle system. FGF2 membrane translocation is thermodynamically driven by PI(4,5)P2-induced membrane insertion of FGF2 oligomers. The latter serve as dynamic translocation intermediates of FGF2 with a subunit number in the range of 8-12 FGF2 molecules. Vectorial translocation of FGF2 across the membrane is governed by sequential and mutually exclusive interactions with PI(4,5)P2 and heparan sulfates on opposing sides of the membrane. Based on atomistic molecular dynamics simulations, we propose a mechanism that drives PI(4,5)P2 dependent oligomerization of FGF2. Our combined findings establish a novel type of self-sustained protein translocation across membranes revealing the molecular basis of the unconventional secretory pathway of FGF2.

https://doi.org/10.7554/eLife.28985.001

Introduction

Fibroblast Growth Factors (FGFs) form a family of more than 20 potent mitogens that stimulate the growth of a wide range of cells including fibroblasts and endothelial cells (Bikfalvi et al., 1997). They are of critical importance for physiological processes such as embryonic development, tissue regeneration, wound repair and hematopoiesis (Bikfalvi et al., 1997). FGF2 is the prototype member of this family that, beyond the functions of FGFs in normal cell growth and differentiation, plays critical roles under pathophysiological conditions (Akl et al., 2016). This is particularly evident in the context of cancer with FGF2 being a major mediator of tumor-induced angiogenesis (Presta et al., 2005). In addition, FGF2 acts as a survival factor that inhibits tumor cell apoptosis by an autocrine secretion-signaling loop (Noh et al., 2014; Pardo et al., 2006). This process is believed to represent a frequent cause of tumor cell resistance against conventional anti-cancer therapies.

As opposed to other FGF family members, FGF2 lacks a signal peptide and is transported into the extracellular space by an ER/Golgi-independent mechanism (La Venuta et al., 2015; Nickel, 2011; Nickel and Rabouille, 2009, 2008). The term unconventional secretion has been used for four different pathways of protein transport towards the plasma membrane and the extracellular space (Rabouille et al., 2012). Two of these pathways concern soluble proteins derived from the cytoplasm that are transported into the extracellular space by either direct protein translocation across the plasma membrane (type I unconventional secretion) or vesicular mechanisms involving endosomal compartments (type III unconventional secretion) (Rabouille et al., 2012; Zhang and Schekman, 2013; Piccioli and Rubartelli, 2013). Based on compelling evidence from both cell-based FGF2 secretion experiments and biochemical in vitro studies, unconventional secretion of FGF2 was shown to follow a type I mechanism that is based on direct protein translocation across the plasma membrane (La Venuta et al., 2015; Zhang and Schekman, 2013; Malhotra, 2013; Nickel, 2005; Schäfer et al., 2004). This process depends on sequential interactions of FGF2 with the phosphoinositide PI(4,5)P2 at the inner leaflet and heparan sulfate proteoglycans at the outer leaflet of the plasma membrane (Nickel, 2011; Nickel and Rabouille, 2009; Nickel, 2007; Temmerman et al., 2008; Zehe et al., 2006). Recruitment by PI(4,5)P2 triggers oligomerization of FGF2 at the inner leaflet resulting in membrane insertion, a key step in FGF2 membrane translocation (La Venuta et al., 2015; Steringer et al., 2012, 2015). The membrane inserted state of FGF2 oligomers has been proposed to be linked to a membrane pore with a toroidal architecture (Steringer et al., 2012). This view is supported by the observation that membrane lipids undergo transbilayer diffusion (Steringer et al., 2012). It appears possible that PI(4,5)P2 molecules themselves redistribute between monolayers when FGF2 oligomers insert into the membrane as it has been reported in other systems (Bucki et al., 2000). Membrane insertion of FGF2 oligomers is stimulated by phosphorylation of tyrosine 81 in FGF2 mediated by the non-receptor tyrosine kinase Tec, a trans-acting factor that is associated with the inner leaflet through an interaction of its PH domain with PI(3,4,5)P3 (Steringer et al., 2012; Ebert et al., 2010; La Venuta et al., 2016). Recently, the integral membrane protein ATP1A1 has been identified as another trans-acting factor in unconventional secretion of FGF2 (La Venuta et al., 2015; Zacherl et al., 2015) FGF2 was shown to directly interact with the cytoplasmic domain of ATP1A1 (La Venuta et al., 2015; Zacherl et al., 2015), however, the precise function of ATP1A1 in FGF2 secretion is unknown. In addition to the cis elements that are important for FGF2 binding to PI(4,5)P2 (K127/R128/K133) and heparan sulfates (K133) as well as Y81 being the target of Tec kinase, membrane insertion of FGF2 oligomers depends on two cysteine residues (C77/C95) on the molecular surface of FGF2 (La Venuta et al., 2015; Müller et al., 2015). Intriguingly, despite a high degree of overall amino acid conservation among FGFs, C77 and C95 are absent from all members of the FGF family carrying signal peptides. C77/C95 have been shown to form intermolecular disulfide bridges, a process that drives efficient oligomerization of FGF2 in a PI(4,5)P2-dependent manner (Müller et al., 2015). Following membrane insertion, cell surface heparan sulfate proteoglycans are required to complete membrane translocation by disassembling FGF2 oligomers at the outer leaflet and trapping of FGF2 in the extracellular space (Nickel, 2007; Zehe et al., 2006; Seelenmeyer et al., 2005). The interaction between FGF2 and cell surface heparan sulfates is mediated by basic residues in the C-terminal part of FGF2 with K133 being an essential component of this binding motif (Nickel and Rabouille, 2009, 2008; Temmerman et al., 2008). Thus, four cis-elements in FGF2 [K127/R128/K133 forming the PI(4,5)P2 binding pocket and K133 involved in FGF2 binding to both PI(4,5)P2 and heparan sulfates, Y81 being the target of Tec kinase and C77/C95 promoting FGF2 oligomerization] and four trans-acting factors [PI(4,5)P2, ATP1A1, Tec kinase and heparan sulfate proteoglycans] have been identified to participate in unconventional secretion of FGF2 from cells. All known trans-acting factors are physically associated with the plasma membrane, the subcellular site of membrane translocation during unconventional secretion of FGF2 (La Venuta et al., 2015).

The molecular mechanism by which FGF2 is secreted from cells might be similar to other proteins secreted by unconventional means. This includes HIV-Tat that has been demonstrated to be secreted from infected T cells in a PI(4,5)P2 dependent manner (Debaisieux et al., 2012; Rayne et al., 2010a, 2010b). Consistently, HIV-Tat was recently found to bind to PI(4,5)P2 concomitant with the formation of membrane pores (Zeitler et al., 2015). Interestingly, a recent study suggested a role for membrane pores in unconventional secretion of Interleukin 1β from macrophages (Martín-Sánchez et al., 2016). While Interleukin 1β on its own cannot form pores in the plasma membrane, its release from macrophages might be coupled to another factor that becomes activated in an inflammasome-dependent manner, gasdermin (Ding et al., 2016). Since gasdermin is capable of forming membrane pores in a PI(4,5)P2 dependent manner, it appears possible that at least one pathway of Interleukin 1β secretion exists that is based upon membrane pores that are formed by gasdermin oligomers upon activation of inflammasomes. Therefore, the molecular mechanisms of unconventional secretion of FGF2, HIV-Tat and even Interleukin 1β might be related in mechanistic terms.

In the current study, we tested key predictions of the FGF2 membrane translocation model as well as investigated the structure-function relationship of membrane-inserted FGF2 oligomers as intermediates in this process. Furthermore, we used atomistic molecular dynamics simulations to shed light on the initial steps of PI(4,5)P2 triggered FGF2 oligomerization and membrane insertion. Using both biochemical and structural approaches, we demonstrate binding to FGF2 of PI(4,5)P2 versus heparin to be mutually exclusive, a key aspect of the vectorial FGF2 membrane translocation model with heparan sulfates forming a molecular trap for FGF2 on cell surfaces. Furthermore, using giant unilamellar vesicles (GUVs) as a model system for the plasma membrane, we define the minimal molecular machinery required for FGF2 membrane translocation in a fully reconstituted inside-out system. These studies revealed only two trans-acting factors to be essential for FGF2 membrane translocation, PI(4,5)P2 on GUV surfaces and long-chain heparins in the lumen of GUVs, the latter being used as mimetics of cell surface heparan sulfate proteoglycans. Consistently, the corresponding cis-elements required for FGF2 binding to PI(4,5)P2 and heparan sulfates as well as the two cysteine residues required for oligomerization and membrane insertion of FGF2 were found essential for FGF2 membrane translocation. In addition, using various kinds of single molecule techniques, we demonstrate membrane-inserted forms of FGF2 to represent highly dynamic oligomers with a subunit number in the range of 8–12 FGF2 molecules. These studies were complemented by molecular dynamics simulations gaining the first insights into the initial steps of PI(4,5)P2 induced oligomerization of FGF2. Most importantly, through simultaneous interactions of FGF2 monomers with multiple PI(4,5)P2 molecules, a high affinity orientation of FGF2 was identified that favors FGF2 dimerization through a C95-C95 disulfide bridge along with additional electrostatic interactions within the dimerization interface. Our combined findings establish both the molecular basis and the minimal molecular machinery required for unconventional secretion of FGF2 from cells that consists of a surprisingly simple set of factors. Our findings further demonstrate the core mechanism of FGF2 membrane translocation to be thermodynamically driven by PI(4,5)P2 dependent FGF2 oligomerization and membrane insertion without a requirement for additional energy sources such as ATP. The combination of PI(4,5)P2 dependent FGF2 oligomerization, membrane insertion and heparan sulfate mediated trapping establishes a new type of self-sustained protein translocation across membranes that explains the molecular basis of how FGF2 is secreted from cells.

Results

Binding to FGF2 of PI(4,5)P2 versus heparin is mutually exclusive

The current model describing the unconventional secretory pathway of FGF2 is based upon direct protein translocation across the plasma membrane. As a prerequisite for directional transport of FGF2 into the extracellular space, a key prediction of this model are sequential and mutually exclusive interactions of FGF2 with PI(4,5)P2 at the inner leaflet and heparan sulfates at the outer leaflet. (La Venuta et al., 2015; Nickel, 2011; Steringer et al., 2015). Thus, employing both structural and biochemical methods, we tested whether interactions of FGF2 with PI(4,5)P2 and heparin (mimicking cell surface heparan sulfates) are mutually exclusive. First, we conducted NMR spectroscopy to compare the binding epitopes in FGF2 for IP3 [the headgroup of PI(4,5)P2] and a defined heparin disaccharide suitable for NMR measurements (Figure 1). IP3 and the heparin disaccharide were titrated with a variant form of FGF2 that is incapable of forming oligomers (FGF2-C77/95S [Müller et al., 2015]; 80 µM; 15N-isotope-labeled). As shown in Figure 1A, B and C, addition of 900 µM IP3 to 15N-labeled FGF2-C77/95S (green contour lines) leads to large chemical shift changes for certain resonances as exemplified by G36, G135 and I145. Titration curves for individual residues within the binding epitope are depicted in Figure 1E. The mean KD value derived from eight resonances was determined as 16.2 ± 4.8 µM. This value is within a range similar to previous data using isothermal titration calorimetry with a KD value in the low micromolar range (Temmerman et al., 2008; Temmerman and Nickel, 2009). In comparison, addition of the small disaccharide variant of heparin led to relatively small changes of an overlapping set of resonances, indicating a weak interaction. For example, as shown in Figure 1A, B and C (purple contour lines), the resonances G36, G135 and I145 are affected by heparin, however, to a significantly smaller extent than the chemical shifts caused by IP3 (green contour lines). Furthermore, titrating IP3 to FGF2 with a pre-bound heparin disaccharide resulted in similar chemical shift changes as for the titration of IP3 to FGF2 alone. This strongly indicates that IP3 is outcompeting the heparin disaccharide under these conditions. Consistently, titrating the heparin disaccharide to FGF2 with a pre-bound IP3 does not outcompete the IP3 binding. Figure 1D shows that the signal of Gly36 of FGF2 is shifting in opposite directions upon IP3 or heparin disaccharide binding, respectively. Addition of the heparin disaccharide to FGF2 with pre-bound IP3 (red contour lines) does not shift the Gly36 signal from ‘IP3’ bound (green contour lines) to ‘heparin bound’ (purple contour lines). When replacing the heparin disaccharide by equimolar amounts of a mixture of long-chain heparins in NMR experiments we obtained a significant loss of signal intensity, indicating slow exchange and strong binding of long-chain heparins to FGF2 (data not shown). Based on the assignment of FGF2 (BMRB code 18995) the binding epitope could be identified as illustrated in Figure 1F. As anticipated, the two epitopes overlap and map to a region of the protein surface that is highly positively charged. These findings are in agreement with the known crystal structure of FGF2 bound to heparin (Faham et al., 1996).

Structural analysis of the FGF2 binding epitopes for IP3 [head group of PI(4,5)P2] and heparin employing NMR spectroscopy.

Enlarged regions of exemplary resonance peaks from a two-dimensional 15N-1H correlation spectrum that are shifted upon addition of IP3 and Heparin or both are shown on panels (A to –D). In panel (E), titration curves for individual 15N-1H resonances are shown. Binding curves were fitted according to a simple two-state binding model and a KD was derived from the mean of the individual titration curves. HSQC for all NMR titration experiments are given in Figure 1—source data 1. In panel F (top), a cartoon of the FGF2 NMR structure is shown (PDB: 1BLD) with the side-chain of residues most significantly shifted upon IP3 binding (blue), heparin binding (yellow) or affected by both binding partners (magenta). In addition, surface mapping of residues shown in panel E using the known structure (PDB: 1BLD) of FGF2 is illustrated in the bottom part of panel (F). Both IP3 and heparin binding epitopes map to the same positively charged region of FGF2 highlighted in blue.

https://doi.org/10.7554/eLife.28985.002
Figure 1—source data 1

Data for Figure 1, panels A-E.

This source file provides HSQC for all NMR titration experiments, of which examples are shown in Figure 1.

https://doi.org/10.7554/eLife.28985.003

The NMR data shown in Figure 1 demonstrate a substantial overlap of the FGF2 binding epitopes for PI(4,5)P2 and heparin and provide initial evidence for these interactions to be mutually exclusive. To challenge these findings using an independent method we conducted biochemical competition experiments. FGF2 was bound to liposomal membranes with a plasma-membrane-like lipid composition including PI(4,5)P2 (Figure 2). Liposomes with FGF2 bound to PI(4,5)P2 were treated with increasing concentrations of either long-chain heparin molecules (consisting of a mixture with various numbers of disaccharide units with high affinity towards FGF2 similar to heparan sulfates on cell surfaces; Figure 2A and B) or a defined low affinity heparin disaccharide (Figure 2C and D). For each condition, bound and unbound FGF2 was analyzed and quantified (Figure 2B and D). These experiments revealed that long-chain heparin molecules directly compete with PI(4,5)P2 for binding to FGF2 with a half-maximal effect at a concentration of ≈5 µM (Figure 2A and B). By contrast, a defined heparin disaccharide was incapable of competing with PI(4,5)P2 for binding to FGF2 (Figure 2C and D). These findings are consistent with the NMR experiments shown in Figure 1D where the low affinity heparin disaccharide fails to replace pre-bound IP3 and provide direct proof for mutually exclusive interactions of FGF2 with PI(4,5)P2 versus heparin. In a cellular context, similar to long-chain heparin molecules used in biochemical reconstitution experiments throughout this study, we conclude that high affinity heparan sulfates on cell surfaces [KD ≈ 100 nM; (Faham et al., 1996)] outcompete PI(4,5)P2 [KD ≈ 5–15 µM; (Temmerman et al., 2008; Temmerman and Nickel, 2009) and Figure 1B] with regard to binding to FGF2. Thus, our data provide direct proof for a central prediction of the FGF2 membrane translocation hypothesis with sequential and competing interactions of FGF2 with PI(4,5)P2 at the inner leaflet and cell surface heparan sulfates at the outer leaflet of the plasma membrane (La Venuta et al., 2015).

Binding of FGF2 to PI(4,5)P2 and heparin is mutually exclusive.

Biochemical analysis employing plasma-membrane-like liposomes either lacking (PM) or containing 2 mol% PI(4,5)P2 (PM +PIP2). Large unilamellar vesicles (LUVs) with bound His-FGF2-Y81pCMF-WT were incubated with increasing concentrations of either long-chain heparins or a defined heparin disaccharide as indicated. After 1 hr of incubation liposome-associated material (bound) was separated from supernatants (unbound). 50% of bound, 13.5% unbound and 14.8% of input material (I) were analyzed by SDS-PAGE. Coomassie-derived signals were quantified and normalized to controls. The fraction of FGF2-Y81pCMF-WT bound to PM-like liposomes containing PI(4,5)P2 in the absence of heparin was set to 100% (red bars in Figure 2B and D). The unbound fraction of FGF2-Y81pCMF-WT was normalized using PM-like liposomes lacking PI(4,5)P2 (blue bars in Figure 2B and D). Mean values with standard deviations are shown (n = 3). Raw and normalized data of individual experiments as well as calculations of mean values with standard deviations are shown in Figure 2—source data 1.

https://doi.org/10.7554/eLife.28985.004
Figure 2—source data 1

Data for Figure 2, panels B and D.

This source file provides raw and normalized data, and calculations of mean values with standard deviations that are plotted in Figure 2, panels B and D.

https://doi.org/10.7554/eLife.28985.005

Biochemical reconstitution of FGF2 membrane translocation with purified components

Based upon our findings demonstrating direct competition between PI(4,5)P2 and long-chain heparins for binding to FGF2 (Figures 1 and 2), we aimed at reconstituting the minimal machinery of FGF2 membrane translocation with purified components (Figure 3). GUVs with a plasma-membrane-like lipid composition were generated that both expose PI(4,5)P2 on their surfaces and contain long-chain heparins in their lumen. These GUVs also contained rhodamine-labelled phosphatidylethanolamine to allow for imaging the lipid bilayer. A phosphomimetic form of wild-type FGF2 [Y81pCMF; (La Venuta et al., 2015; Steringer et al., 2012; Ebert et al., 2010; La Venuta et al., 2016)] was used as a GFP fusion protein (FGF2-Y81pCMF-GFP) to monitor FGF2 membrane translocation into the lumen of GUVs. In addition, a small fluorescent tracer (Alexa647) was used to detect the formation of membrane pores (La Venuta et al., 2015; Steringer et al., 2012). In Figure 3, two experimental conditions were compared analyzing GUVs that either contained (Figure 3A) or lacked (Figure 3B) luminal long-chain heparin. In both cases, FGF2-Y81pCMF-GFP bound efficiently to the membrane surface due to the presence of PI(4,5)P2. By contrast, only in case luminal heparin was included, FGF2-Y81pCMF-GFP membrane translocation into the lumen of GUVs was observed. This phenomenon was quantified based upon measuring GFP fluorescence intensity in the lumen, at the membrane and in the vesicle exterior as described in ‘Methods’. In the given example, the ratio between luminal and external fluorescence was 4.93 in the presence of luminal heparin (Note increased fluorescence intensity on the luminal side in sub-panel e of panel A) and 1.33 in the absence of luminal heparin (Note equal fluorescence intensities on the luminal and external side in sub-panel e of panel B). In panels C and D of Figure 3, cross sections of 3D reconstructions are shown visualizing the dependence of FGF2-GFP membrane translocation on the presence of luminal heparin with a spatial view into the interior of GUVs. Based on a threshold of 1.6, a statistical analysis of 20 to 120 GUVs from at least three independent experiments was conducted providing the percentages of GUVs with FGF2-Y81pCMF-GFP being enriched in their lumen for each experimental condition. This analysis revealed that FGF2-Y81pCMF-GFP membrane translocation occurred in about 23% of the GUVs with luminal heparin (Figure 9). By contrast, less than 2% of GUVs without luminal heparin contained FGF2-Y81pCMF-GFP in their lumen (Figure 9). Simultaneously, the small fluorescent tracer Alexa647 was used to monitor membrane pore formation as indicated by equal fluorescence intensities in the lumen and the exterior of GUVs. Under the two conditions shown in Figure 3, membrane pore formation did not depend on the presence of heparin in the lumen of GUVs. As shown in Figure 9, membrane pore formation occurred in about 33% to 41% of all GUVs, both in the absence and the presence of luminal heparin. These results are consistent with previous findings demonstrating that membrane insertion of FGF2 oligomers concomitant with pore formation does not depend on luminal heparin (Steringer et al., 2012). GUVs containing luminal FGF2-Y81pCMF-GFP in the absence of membrane pore formation were undetectable.

Reconstitution of FGF2 membrane translocation with purified components.

Giant unilamellar vesicles with a plasma membrane-like lipid composition containing PI(4,5)P2 were prepared in the presence (panel A and C) or absence (panel B and D) of long-chain heparins. Rhodamine-PE was incorporated into the lipid bilayer during GUV preparation as membrane marker. After removal of excess heparin by low speed centrifugation, GUVs were incubated with FGF2-Y81pCMF-GFP (200 nM) and a small fluorescent tracer (Alexa647). Following 180 min of incubation luminal penetration of GUVs by FGF2-Y81pCMF-GFP and small tracer molecules was analyzed by confocal microscopy (scale bar = 10 µm). GUVs were analyzed in all three channels using the plugin ‘Radial profile’ of the ImageJ software as explained under ‘Materials and methods’. Profile plots of normalized integrated intensities around concentric circles as a function of distance from the center of the GUV are given in relative fluorescence units (RFU). FGF2 membrane translocation is indicated by increased GFP fluorescence intensity in the lumen of GUVs compared to the exterior as exemplified in sub-panel e of panel A. The dependence of FGF2-Y81pCMF-GFP membrane translocation on luminal heparin is further documented by 3D reconstruction images (panels C and D) providing a spatial view into the GUV interior.

https://doi.org/10.7554/eLife.28985.006

To challenge the results shown in Figure 3, we conducted experiments where high-affinity long-chain heparins were directly compared with the low-affinity heparin disaccharide introduced in Figures 1 and 2. While FGF2-Y81pCMF-GFP membrane binding, membrane pore formation (Alexa647 intensity ratio lumen/exterior = 0.96) and FGF2 membrane translocation (GFP intensity ratio lumen/exterior = 6.44) was fully functional for GUVs containing luminal long-chain heparins (Figure 4A), FGF2-Y81pCMF-GFP membrane translocation was impaired for GUVs containing the heparin disaccharide (Figure 4B; GFP intensity ratio lumen/exterior = 1.07). However, as expected, membrane binding and pore formation was normal under these conditions (Figure 4B; Alexa647 intensity ratio lumen/exterior = 0.95). As depicted in Figure 9, a statistical analysis revealed that substitution of luminal long-chain heparins by a heparin disaccharide causes a drop from about 23% to only 5% of GUVs positive for FGF2-Y81pCMF-GFP membrane translocation. Our combined findings from Figures 3 and 4 establish FGF2 membrane translocation in a minimal system employing GUVs with a plasma-membrane-like lipid composition including PI(4,5)P2 and luminal long-chain heparins mimicking heparan sulfates from cell surfaces.

FGF2 membrane translocation depends on luminal long-chain heparins that cannot be substituted by low affinity heparin disaccharides.

Giant unilamellar vesicles with a plasma membrane-like lipid composition containing PI(4,5)P2 were prepared that contain either long-chain heparins (panel A) or a defined heparin disaccharide (panel B). Incubation with FGF2-Y81pCMF-GFP and data analysis were conducted as described in the legend to Figure 3 and under ‘Materials and methods’.

https://doi.org/10.7554/eLife.28985.007

Dependence of FGF2 membrane translocation on PI(4,5)P2

To analyze a requirement for PI(4,5)P2 on FGF2-Y81pCMF-GFP induced membrane pore formation and membrane translocation (Figure 5), four types of GUVs were used all of which contained luminal long-chain heparins but differed with regard to lipid composition and the presence of PI(4,5)P2 on their membrane surfaces. Similar to the experiments shown in Figures 3 and 4, FGF2-Y81pCMF-GFP membrane binding, membrane pore formation (Alexa647 intensity ratio lumen/exterior = 0.97) and FGF2 membrane translocation (GFP intensity ratio lumen/exterior = 3.27) could be observed when GUVs contained PI(4,5)P2 on their membrane surface (Figure 5A). By contrast, when cholesterol was omitted from the plasma membrane like lipid composition, PI(4,5)P2 dependent membrane recruitment was significantly impaired (Figure 5B). This observation is consistent with earlier findings suggesting an impact of cholesterol dependent microdomain formation on the efficiency of FGF2 binding to PI(4,5)P2 (17,33). In turn, under these conditions, FGF2-Y81pCMF-GFP membrane pore formation (Figure 5B; Alexa647 intensity ratio lumen/exterior = 0.05) and FGF2 membrane translocation (Figure 5B; GFP intensity ratio lumen/exterior = 0.84) could not be observed. These findings were corroborated by a statistical analysis depicted in Figure 9 revealing that only about 6% of GUVs were positive for FGF2-Y81pCMF-GFP membrane translocation when cholesterol was removed from a plasma membrane like lipid composition that included PI(4,5)P2.

PI(4,5)P2 is required for both membrane pore formation and FGF2 membrane translocation.

Giant unilamellar vesicles were prepared with a plasma membrane-like lipid composition either containing PI(4,5)P2 (panel A), containing PI(4,5)P2 but lacking cholesterol (panel B), lacking PI(4,5)P2 (panel C) or containing a Ni-NTA lipid substituting PI(4,5)P2 (panel D). All four types of GUVs contained luminal long-chain heparins. Incubation with FGF2-Y81pCMF-GFP and data analysis were conducted as described in the legend to Figure 3 and under ‘Materials and methods’. Note increased GFP fluorescence in the lumen of GUVs as exemplified in sub-panel e of panel A indicating FGF2-Y81pCMF-GFP membrane translocation.

https://doi.org/10.7554/eLife.28985.008

In the absence of PI(4,5)P2 on the surface of GUVs with luminal long-chain heparins, a complete failure of FGF2-Y81pCMF-GFP membrane recruitment was observed (Figure 5C). This in turn resulted in a lack of membrane pore formation (Alexa647 intensity ratio lumen/exterior = 0.05) and FGF2-Y81pCMF-GFP translocation into the lumen of GUVs (GFP intensity ratio lumen/exterior = 0.7). As shown in Figure 9, in the absence of PI(4,5)P2, less than 2% of GUVs were characterized by membrane pores and FGF2 membrane translocation was undetectable (Figure 9). The fourth type of GUVs in this set of experiments was characterized by substitution of PI(4,5)P2 with a Ni-NTA lipid (Figure 5D). This experimental condition has previously been demonstrated to allow for efficient recruitment of His-tagged FGF2 to membrane surfaces. However, under these conditions, a substantial reduction of membrane pore formation was observed (Steringer et al., 2012). This phenotype is not due to a general block in FGF2 oligomerization but is likely to be related to a requirement for PI(4,5)P2 to facilitate FGF2-dependent formation of a toroidal membrane pore that is characterized by strong membrane curvature (La Venuta et al., 2015). Indeed, in the example given in Figure 5D, the tracer intensity ratio (lumen/exterior) was found to be 0.05 demonstrating pore formation to depend on PI(4,5)P2, even under conditions where membrane recruitment of FGF2-Y81pCMF-GFP is mediated by other means. Consistently, in the statistical analysis shown in Figure 9, only about 10% of GUVs containing the Ni-NTA lipid to recruit His-tagged FGF2-GFP were found to contain membrane pores. This phenotype was accompanied by an almost complete failure of His-tagged FGF2-Y81pCMF-GFP to translocate into the lumen of GUVs when PI(4,5)P2 was substituted by a Ni-NTA lipid (GFP intensity ratio lumen/exterior = 1.15 in Figure 5D). As shown in Figure 9, this was confirmed by a statistical analysis with less than 1% of Ni-NTA GUVs containing luminal FGF2-Y81pCMF-GFP.

Our combined findings from Figures 3, 4 and 5 along with their quantification given in Figure 9 define the minimal molecular machinery required for FGF2 membrane translocation. Two trans-acting factors that previously have been demonstrated to be required for FGF2 secretion from cells (La Venuta et al., 2015; Temmerman et al., 2008; Zehe et al., 2006) are essential for this process, PI(4,5)P2 and long-chain heparins, the latter mimicking cell surface heparan sulfates.

Dependence of FGF2 membrane translocation on cis-elements in FGF2

Following the analysis of trans-acting factors for FGF2 membrane translocation in a minimal system, we aimed at studying cis-elements in FGF2 known to be required for FGF2 secretion from cells (La Venuta et al., 2015). We used three variant forms of purified FGF2-Y81pCMF-GFP that either have a defect in binding to PI(4,5)P2 [K127Q/R128Q; (Nickel, 2011; Temmerman et al., 2008; Temmerman and Nickel, 2009)], a defect in binding to both PI(4,5)P2 and long-chain heparins [K127Q/R128Q/K133Q; (Nickel, 2011; Temmerman et al., 2008; Temmerman and Nickel, 2009)] or a defect in oligomerization and membrane pore formation [C77A/C95A; (La Venuta et al., 2015; Müller et al., 2015)]. The various forms of FGF2-Y81pCMF-GFP were characterized with regard to binding to both long-chain heparins immobilized on beads (Figure 6A and B) and PI(4,5)P2 as part of liposomes using a flow cytometry setup [Figure 6C; (Temmerman et al., 2008; Temmerman and Nickel, 2009)]. In addition, the processed form of interleukin 1β fused to GFP (IL-1β-GFP) and GFP alone were taken along as control proteins. Binding of the various proteins to long-chain heparins and PI(4,5)P2 was compared to the wild-type form of FGF2-Y81pCMF-GFP. As shown in Figure 6A and B, substitution of C77 and C95 by alanines (FGF2-Y81pCMF-C77/95A-GFP) did not affect binding to long-chain heparins. By comparison, FGF2-Y81-pCMF-C77/95A-GFP showed reduced binding to PI(4,5)P2 (Figure 6C). This is caused by a reduction in binding avidity due to the inability of this FGF2 variant form to oligomerize. By contrast, the K127Q/R128Q/K133Q form of FGF2-Y81pCMF-GFP did neither bind to long-chain heparins (Figure 6A and B) nor to PI(4,5)P2 (Figure 6C). A differential phenotype was observed for the K127Q/R128Q of FGF2-Y81pCMF-GFP which did bind to long-chain heparins (Figure 6A and B) but failed to interact with PI(4,5)P2 (Figure 6C). Finally, neither IL-1β-GFP nor GFP alone did bind to long-chain heparins or PI(4,5)P2 (Figure 6A, B and C).

Analysis of FGF2 variant forms with differential defects in binding to PI(4,5)P2 and heparin.

The variant forms of FGF2-Y81pCMF-GFP indicated were tested for binding to heparin beads (panels A and B) and PI(4,5)P2 (panel C). Heparin sepharose beads were incubated with the FGF2 fusion proteins indicated. Bound and unbound material was separated by centrifugation. Bound proteins were eluted with SDS sample buffer and analysed by SDS-PAGE and Coomassie staining [5% input (I), 20% bound (B) and 5% unbound (UB)]. Signals were quantified using a Li-COR Odyssey infrared imaging system. Mean values with standard deviations of three independent experiments are shown (panel B). Raw and normalized data of individual experiments as well as calculations of mean values with standard deviations are shown in Figure 6—source data 1. Binding of the FGF2 fusion proteins to PI(4,5)P2 contained in plasma membrane-like liposomes was assessed using a flow-cytometry assay (Temmerman et al., 2008; Temmerman and Nickel, 2009) (panel C). Data were normalized by defining binding of FGF2-Y81pCMF-GFP to PI(4,5)P2 as 100% binding efficiency. Mean values with standard deviations are shown (n = 4). Consider Figure 6—source data 1 for more details.

https://doi.org/10.7554/eLife.28985.009
Figure 6—source data 1

Data for Figure 6, panels B and C.

This source file provides raw and normalized data of SDS-PAGE and Coomassie analysis (Li-COR Odyssey infrared imaging system) with calculations of mean values with standard deviations that are plotted in Figure 6B. Furthermore, corrected fluorescent units of individual experiments (each experiment represents the mean of a duplicate) plus calculations of mean values with standard deviations of FACS-based experiments that are shown in Figure 6C are given.

https://doi.org/10.7554/eLife.28985.010

The FGF2 variant forms and control proteins described above were tested for activity in membrane pore formation and membrane translocation using GUVs containing both PI(4,5)P2 and long-chain heparins in their lumen (Figures 7 and 8). The wild-type form of FGF2-Y81pCMF-GFP was used as a reference in both assays (Figures 7A and 8A; GFP intensity ratio lumen/exterior = 7.14 and 6.3, respectively). As expected, neither the K127Q/R128Q/K133Q form nor the K127Q/R128Q form of FGF2-Y81pCMF-GFP were capable of binding to PI(4,5)P2-containing GUVs. Consistently, neither the K127Q/R128Q/K133Q (GFP intensity ratio lumen/exterior = 0.61; Figure 7B) nor the K127Q/R128Q form (GFP intensity ratio lumen/exterior = 0.6; Figure 7C) of FGF2-Y81pCMF-GFP were capable of translocating into the lumen of GUVs containing both PI(4,5)P2 and long-chain heparins. This was confirmed by the statistical analysis shown in Figure 9 with membrane translocation being undetectable for K127Q/R128Q/K133Q and below 2% for K127Q/R128Q, respectively. As shown in the example given in Figure 7D, the C77/95A form of FGF2-Y81pCMF-GFP was characterized by a defect in membrane pore formation (Alexa647 intensity ratio lumen/exterior = 0.05) using GUVs containing both PI(4,5)P2 and long-chain heparins. Similarly, FGF2-Y81pCMF-C77/95A-GFP translocation into the lumen was not observed (GFP intensity ratio lumen/exterior = 0.97; Figure 7D). As shown in Figure 9, about 10% of GUVs incubated with FGF2-Y81pCMF-C77/95A-GFP contained membrane pores and luminal FGF2-Y81pCMF-C77/95A-GFP, a significant reduction compared to FGF2-Y81pCMF-GFP. Finally, we tested two control proteins, IL-1β-GFP (Figure 8B) and GFP alone (Figure 8C). Both of them did not show any activity with regard to binding to PI(4,5)P2, membrane pore formation (Alexa647 intensity ratio lumen/exterior = 0.03 and 0.05, respectively) and membrane translocation into the lumen of GUVs (GFP intensity ratio lumen/exterior = 0.56 and 0.48, respectively). These findings were confirmed by the statistical analysis shown in Figure 9. In conclusion, in addition to the trans-acting factors PI(4,5)P2 and long-chain heparins as mimetics of cell surface heparan sulfates, FGF2 membrane translocation depends on cis-elements that allow for binding to PI(4,5)P2 (K127/R128), binding to long-chain heparins (K133) and are required for oligomerization and membrane pore formation (C77/C95).

FGF2 membrane translocation depends on cis-elements mediating binding to PI(4,5)P2 and heparin as well as driving FGF2 oligomerization and membrane pore formation.

Giant unilamellar vesicles with a plasma membrane-like lipid composition containing both PI(4,5)P2 and luminal long-chain heparins were prepared as described in the legend to Figure 3 and under ‘Materials and methods’. GUVs were incubated with variant forms of FGF2-Y81pCMF-GFP as indicated. These included the wild-type form (panel A), the K127Q/R128Q/K133Q form deficient in binding to PI(4,5)P2 and heparin (panel B), the K127Q/R128Q form deficient in binding to PI(4,5)P2 (panel C) and the C77A/C95A form deficient in oligomerization and membrane pore formation (panel D). Incubation conditions and data analysis were conducted as described in the legend to Figure 3 and under ‘Materials and methods’. Note increased GFP fluorescence in the lumen of GUVs as exemplified in sub-panel e of panel A indicating membrane translocation of the wild-type form of FGF2-Y81pCMF-GFP.

https://doi.org/10.7554/eLife.28985.011
Interleukin 1β, a structural homologue of FGF2, is incapable of binding to PI(4,5)P2, membrane pore formation and membrane translocation.

Giant unilamellar vesicles with a plasma membrane-like lipid composition containing both PI(4,5)P2 and luminal long-chain heparins were prepared as described in the legend to Figure 3 and under ‘Materials and methods’. Incubation conditions using the wild-type form of FGF2-Y81pCMF-GFP (panel A), the mature form of Interleukin 1β-GFP (panel B) and GFP as control protein (panel C) as well as data analysis were conducted as described in the legend to Figure 3 and under ‘Materials and methods’. Note increased GFP fluorescence in the lumen of GUVs as exemplified in sub-panel e of panel A indicating membrane translocation of the wild-type form of FGF2-Y81pCMF-GFP.

https://doi.org/10.7554/eLife.28985.012
Quantification and statistical analysis of FGF2-Y81pCMF membrane translocation and its dependence on both cis-elements and trans-acting factors known to be required for FGF2 secretion from cells.

A quantitative analysis of membrane translocation and pore formation by the various proteins indicated was conducted based upon the experiments shown in Figures 3, 4, 5, 7 and 8. Various types of GUVs with a plasma membrane-like composition were used that differed with regard to the presence of the components indicated. For all conditions, data were derived from at least three independent experiments each of which involved the analysis of 20–120 GUVs per experimental condition. Gray bars indicate the percentage of GUVs with membrane pores with a ratio of Alexa647 tracer fluorescence in the lumen versus the exterior of ≥0.6. Green bars indicate the percentage of GUVs where membrane translocation of GFP-tagged proteins had occurred with a ratio of GFP fluorescence in the lumen versus the exterior of ≥1.6 being used as a threshold value. Standard deviations are shown (n ≥ 15 for experiments shown in Figure 3 and n ≥ 3 for all other conditions shown in Figures 4, 5, 7 and 8). Detailed information on each individual experiment is provided in Figure 9—source data 1.

https://doi.org/10.7554/eLife.28985.013
Figure 9—source data 1

Data for Figure 9.

This source file provides additional information on each sample including number of experiments performed, number of GUVs analyzed in each single experiment, original counts and percentage of GUVs positive for the small tracer or FGF2 translocation, respectively, as well as calculations of mean values with standard deviations plotted in Figure 9.

https://doi.org/10.7554/eLife.28985.014

Oligomeric state of membrane inserted FGF2 translocation intermediates in supported lipid bilayers

Beyond the reconstitution of FGF2 membrane translocation with purified components, we aimed at insight into the structure function relationship of membrane-inserted FGF2 oligomers, the key intermediates in unconventional secretion of FGF2 from cells. In particular, we analyzed the subunit number of membrane associated FGF2 oligomers. In a first approach, we used supported lipid bilayers (SLB) containing PI(4,5)P2 to determine the oligomeric state of membrane inserted FGF2 complexes. To obtain high contrast single molecule detection we used FGF2 variant forms carrying a Halo-tag. This allowed for protein labeling with a bright and photo-stable fluorophore (Abberior StarRed) with a degree of labeling of 0.9 per FGF2 monomer. A three-step protocol was used to analyze the oligomeric state distribution of FGF2 translocation clusters (Figure 10; for details see ‘Materials and methods’). Two variant forms of FGF2, FGF2-Y81pCMF-HALO-StarRed and FGF2-Y81pCMF-C77/95A-HALO-StarRed were loaded onto SLBs containing 2 mol% PI(4,5)P2 at a final concentration of 100 nM. While FGF2-Y81pCMF-HALO-StarRed efficiently bound to SLBs, (Figure 10D), binding of FGF2-Y81pCMF-C77/95A-HALO-StarRed was significantly reduced (Figure 10J). These findings are consistent with the data shown in Figures 6 and 7 and reflect the inability of FGF2-Y81pCMF-C77/95A-HALO-StarRed to oligomerize (see above). In a second step, following 10 min of incubation, SLBs were washed with a buffer containing 150 mM NaCl to remove unbound proteins (Figure 10B,E,H and K). At this stage, fluorescence recovery after photo-bleaching (FRAP) was used to determine the mobility of membrane bound FGF2 species (FRAP data available in Figure 11—source data 1). This analysis revealed the majority (98%) of FGF2-Y81pCMF-HALO-StarRed to be highly mobile with a diffusion constant of D = 0.4 µm2/s. These findings indicate that after a short incubation with SLBs of only 10 min, most of FGF2-Y81pCMF-HALO-StarRed did not oligomerize into higher order structures which were integrated in the membrane. In a third step, the highly mobile fraction of FGF2-Y81pCMF-HALO-StarRed monomers bound to PI(4,5)P2 was removed by high salt treatment (500 mM NaCl) (Figure 10C,F,I and L). The remaining population of the protein was largely immobile suggesting membrane insertion of oligomers with contacts to the glass support of SLBs. In contrast, we did not observe significant immobilization of FGF2-Y81pCMF-C77/95A-HALO-StarRed which is due to a failure of oligomerization and membrane insertion (Müller et al., 2015).

Binding and membrane insertion of FGF2-Halo-StarRed fusion proteins into supported lipid bilayers containing PI(4,5)P2.

FGF2-Y81pCMF-Halo-StarRed (Panels A–F) and FGF2-Y81pCMF-C77/95A-Halo-StarRed (panels G–L) were added at a final concentration of 100 nM to supported lipid bilayers (SLBs) containing 68 mol% POPC, 30 mol% cholesterol and 2 mol% PI(4,5)P2 plus trace amounts of DPPE-OregonGreen to image the bilayer. FGF2-Y81pCMF-Halo-StarRed and FGF2-Y81pCMF-C77/95A-Halo-StarRed were bound to SLBs (panels A, D, G and J) followed by a 150 mM NaCl washing procedure (panels B, E, H and K). In a final step, a 500 mM salt wash was applied to remove FGF2 monomers (panels C, F, I and L). FGF2-Y81pCMF-Halo-StarRed and FGF2-Y81pCMF-C77/95A-Halo-StarRed bound to SLBs were imaged as explained in ‘Materials and methods’.

https://doi.org/10.7554/eLife.28985.015

In order to determine the oligomeric state of the immobilized FGF2-Y81pCMF-HALO-StarRed clusters, 2D confocal images were recorded with a long pixel dwell-time of 1 ms. This resulted in high signal to background ratio images without significant photo-bleaching (Figure 11A). Single molecule localization was used to automatically determine the brightness of each cluster via 2D Gaussian fitting. The brightness of a single HALO-StarRed was determined by following the same imaging and fitting procedure with immobilized mEGFP-HALO-StarRed (data not shown). Finally, the number of monomers in each immobilized FGF2 cluster was estimated by normalizing the cluster brightness to the brightness of a single HALO-StarRed (Figure 11A). A Gaussian mixture analysis in MATLAB was used to estimate the number of sub-populations in the whole distribution of more than 1000 clusters from six independent samples. The fit yielded four components with 3, 6, 11 and 17 monomers per cluster, respectively, with the majority of the membrane inserted population after 10 min of incubation being in the trimeric and hexameric state (Figure 11B).

Single molecule imaging and brightness analysis of FGF2-Y81pCMF-Halo-StarRed to determine the oligomeric state of membrane inserted FGF2 clusters.

(A) Immobile FGF2-Y81pCMF-Halo-StarRed clusters associated with SLBs following a high salt wash (Figure 10F) were imaged by confocal microscopy. The brightness of individual clusters was determined by fitting a 2D Gaussian to each diffraction limited spot using a single molecule tracking software in MATLAB. The number of monomers in each cluster (right part of image) was estimated by normalizing the brightness of each cluster to the brightness of monomeric HALO-StarRed). Original peak intensities and cluster analyses are available in Figure 11—source data 1. (B) Gaussian mixture analysis of the oligomeric state of membrane inserted FGF2-Y81pCMF-Halo-StarRed (>1000 clusters from six independent experiments). The distribution of monomers per cluster was complex. A Gaussian mixture analysis found 4 components with 3, 6, 11 and 17 monomers per cluster.

https://doi.org/10.7554/eLife.28985.016
Figure 11—source data 1

Data for Figure 11.

This source file provides peak intensities that were used to determine the number of monomers per cluster and the brightness of the StarRed monomer control (see ‘Material and methods’ for detail). In addition, original data of FRAP experiments mentioned in the text to Figure 11 are shown.

https://doi.org/10.7554/eLife.28985.017

Functional correlation of the oligomeric state of membrane associated FGF2 translocation intermediates and membrane pore formation

Beyond our studies using SLBs (see above) we aimed at using an experimental system that provides direct insight into the structure-function relationship of membrane inserted FGF2 translocation intermediates. Therefore, we conducted experiments to directly correlate oligomerization of FGF2-Y81pCMF-GFP with membrane pore formation in GUVs with a plasma membrane like lipid composition. Using z-scan fluorescence correlation spectroscopy [FCS; (Benda et al., 2003)], we combined a brightness and diffusion analysis of individual FGF2 oligomers with the membrane pore assay introduced in Figure 3 (Steringer et al., 2012). The average oligomeric state (Figure 12A, C and E) and diffusion constants (Figure 12B, D and F) of membrane associated FGF2-Y81pCMF-GFP of individual GUVs were plotted as a function of protein concentration on membrane surfaces [c(FGF2-Y81pCMF-GFP)]. In addition, each GUV was classified regarding membrane pore formation based upon luminal penetration of a small fluorescent tracer. GUVs without membrane pores were characterized by a low average surface protein concentration of c(FGF2-Y81pCMF-GFP)=0.17 ± 0.27 nmol/m2 (Figure 12A, open circles). By contrast, GUVs containing membrane pores were found to have a four-fold higher average protein surface concentration of c(FGF2-Y81pCMF-GFP)=0.71 ± 0.61 nmol/m2 (Figure 12C; closed circles). The average oligomeric state measured on GUVs containing membrane pores was 9.23 ± 2.96 for c(FGF2-Y81pCMF-GFP) larger than 0.4 nmol/m2 (57% of GUVs) and 5.24 ± 3.5 for c(FGF2-Y81pCMF-GFP) smaller than 0.4 nmol/m2 (43% of GUVs). By contrast, 90% of the GUVs without membrane pores were characterized by an average oligomeric state of 5.22 ± 2.86 with c(FGF2-Y81pCMF-GFP)≤0.4 nmol/m2. These experiments revealed a clear correlation between c(FGF2-Y81pCMF-GFP), the formation of higher oligomers and the probability of membrane pore formation.

Functional correlation of the oligomeric state of membrane associated FGF2 translocation intermediates and membrane pore formation.

Giant unilamellar vesicles with a plasma membrane-like lipid composition containing PI(4,5)P2 and the membrane tracer DOPE-Atto633 were prepared as described in ‘Materials and methods’. After pre-incubation with either the wild-type form of FGF2-Y81pCMF-GFP (panels A-D, black circles) or FGF2-Y81pCMF-C77/95A-GFP (panels E and F; red circles) for at least 30 min, z-scan FCS measurements using 515/50 nm (FGF2-GFP) and 697/58 nm (DOPE-Atto633) emission channels were conducted on single GUVs. The small free tracer AlexaFlour532 was added to the buffer in order to visualize FGF2 membrane pore formation. Accordingly, GUVs were classified into two groups with (panels C and D, filled circles) and without membrane pores (panels A, B, E, and F; empty circles). Z-scan measurements and analyses are described in detail under ‘Materials and methods’. Average oligomeric state values (panels A, C, E) and diffusion constants (panels B, D, F) were plotted as a function of protein surface concentration. A total of 60 individual GUVs incubated with FGF2-Y81pCMF-WT-GFP (panels A-–D) and 9 GUVs incubated with FGF2-Y81pCMF-C77/95A (panels E and F) were analyzed. Additional data of Z-scan FCS for each individual GUV, monomer control reference measurements as well as calculations of mean values with standard deviations are provided in Figure 12—source data 1.

https://doi.org/10.7554/eLife.28985.018
Figure 12—source data 1

Data for Figure 12.

This source file provides extensive information of each individual GUV analyzed by Z-scan FCS including pore formation state, particle number, brightness of a cluster, oligomeric state, two independent determinations of surface protein concentration, diffusion times, and diffusion coefficients. Furthermore, reference measurements of monomer units (see ‘Material and methods’ for details) as well as calculations of mean values with standard deviations that are plotted in Figure 12 are given.

https://doi.org/10.7554/eLife.28985.019

In order to verify the determination of average oligomeric state values derived from brightness analyses, diffusion measurements were conducted for FGF2-Y81pCMF-GFP as shown in Figure 12B and D. At low surface protein concentrations [c(FGF2-Y81pCMF-GFP)≤0.4], FGF2-Y81pCMF-GFP clusters moved slightly slower in GUVs with membrane pores (D = 2.5 ± 0.8 µm2/s) than in GUVs without membrane pores (D = 2.8 ± 0.5 µm2/s). At high surface protein concentrations [c(FGF2-Y81pCMF-GFP)>0.4], the diffusion of FGF2-Y81pCMF-GFP membrane pore forming FGF2 oligomers was strongly decreased (D = 1.7 ± 0.4 µm2/s). These data are consistent with the determination of oligomeric states shown in Figure 12A and C.

Finally, the methodology we used to determine oligomeric states and diffusion constants of FGF2-Y81pCMF-GFP on membrane surfaces was further verified using an FGF2 mutant that binds to PI(4,5)P2 containing membranes, however, is incapable of oligomerizing and forming membrane pores [FGF2-Y81pCMF-C77/95A-GFP; Figures 6 and 7; (Müller et al., 2015)]. Indeed, the average oligomeric state of FGF2-Y81pCMF-C77/95A-GFP on membrane surfaces was determined to be 1.17 ± 0.5 along with a high diffusion constant of D = 4.8 ± 1.0 µm2/s (Figure 12E and F). Consistently, GUVs containing membrane pores following incubation with FGF2-Y81pCMF-C77/95A-GFP were undetectable. These results demonstrate that this FGF2 variant form remained a monomer on the membrane surface of GUVs, which is consistent with previous studies and validates our experimental setup to determine the oligomeric state of FGF2-Y81pCMF-GFP on membrane surfaces. Thus, based on the data shown in Figure 12, FGF2-Y81pCMF-GFP oligomers in the membrane of GUVs are characterized by a range of about 8 to 12 subunits of FGF2 under conditions where membrane pore formation can be observed.

Simultaneous interactions of FGF2 monomers with several PI(4,5)P2 molecules stabilizes an FGF2 orientation that triggers dimerization

To complement experiments correlating the oligomeric state of FGF2 assemblies with membrane insertion and pore formation (Figures 10, 11 and 12), we conducted atomistic molecular dynamics simulations to gain insight into the initial molecular events that trigger FGF2 oligomerization in a PI(4,5)P2 dependent manner (Figures 13 and 14). A starting point of this approach were biochemical experiments demonstrating K127, R128 and K133 to be part of a binding pocket that recruits the headgroup of PI(4,5)P2, IP3 (17,33) (Figure 6). However, it is possible that additional residues play a role and the binding stoichiometry between FGF2 and PI(4,5)P2 has not been determined. Furthermore, beyond the formation of intermolecular disulfide bridges involving C77 and C95 (Müller et al., 2015), there is so far no information on protein-protein interfaces that form during PI(4,5)P2 induced oligomerization of FGF2. Thus, we carried out a series of atomistic molecular dynamics simulations (Figure 13 and Video 1; system M1 with five repeats) placing FGF2 monomers with a distance of 1.5 nm above the membrane surface in different orientations (Figure 13 and Figure 13—figure supplement 1). The simulation data showed that FGF2 readily makes contacts with the membrane surface and undergoes changes in orientation as it binds to the membrane. In line with previous biochemical and structural studies (La Venuta et al., 2015; Temmerman et al., 2008; Steringer et al., 2012; Müller et al., 2015), we observed spontaneous binding of FGF2 to PI(4,5)P2 through the key binding site residues (K127, R128, K133) (Video 1). However, as discussed below, a number of other residues were also found to contribute to the binding process depending on the number of PI(4,5)P2 molecules interacting with FGF2 in its vicinity.

Figure 13 with 8 supplements see all
FGF2 orientation on the membrane surface.

(A) High-affinity orientation of FGF2, showing all the known PI(4,5)P2-binding site residues (K127, R128, K133) as well as the previously undetermined binding site residues (K34, K137, K143). (B) Low-affinity orientation of FGF2 in which the binding site residues lose contact with PI(4,5)P2 and point away from the PI(4,5)P2 head groups. FGF2 is rendered as green cartoon, and its C95 and C77 residues are shown as van der Waals (vdW) spheres and highlighted by text in the figure. The key binding pocket residues (K127, R128, K133) are shown as blue vdW spheres, and the additional binding site residues (K34, K137, K143) are shown as purple vdW spheres. Lipids are colored as gray vdW spheres (POPC phosphate atoms), red vdW spheres [PI(4,5)P2 bisphosphates], cyan vdW spheres [inositol ring in PI(4,5)P2], orange vdW spheres [phosphate linking the fatty acid chains and the inositol ring in PI(4,5)P2], and white vdW spheres [fatty acid chains in PI(4,5)P2]. Water molecules and ions are not shown for clarity.

https://doi.org/10.7554/eLife.28985.020
Figure 14 with 3 supplements see all
FGF2 trimer configurations.

Snapshots representing the most populated structures in FGF2 trimer simulations (system T). (Panel A) depicts the top view of the FGF2 trimer aggregate with C95 – C95 and C77 – C77 interfaces labeled, where C95 and C77 are colored as van der Waals spheres. The trimer is split into two dimer interfaces shown in (panels B and C). (B) The interface residues involved in C95 – C95 disulfide-linked dimers. (C) The interface residues involved in C77 – C77 disulfide-linked dimers. The interface residues are depicted in stick representation, where negatively charged residues (D, E) are colored as red and positively charged residues (K, R) as blue. The PI(4,5)P2 binding pocket residues (K127, R128, K133) are rendered as blue van der Waals spheres. For clarity, POPC, PI(4,5)P2, cholesterol, water, and ions are not shown.

https://doi.org/10.7554/eLife.28985.029
Video 1
FGF2 binding simultaneously to several PI(4,5)P2 molecules (system M1).

FGF2 binds to the membrane surface in the high-affinity orientation and interacts simultaneously with multiple PI(4,5)P2, including the key residues in the main binding pocket (K127, R128, K133) and the residues in its vicinity (K34, K137, K143). The final frame of the video at 913 ns matches the snapshot in Figure 13A taken though from a slightly different perspective. The color coding used for the protein and PI(4,5)P2 is consistent with Figure 13. In the membrane, cholesterols are shown in light purple, the phosphorous atoms of POPC molecules are depicted as a transparent surface, and for clarity’s sake water and ions are not shown.

https://doi.org/10.7554/eLife.28985.033

We further identified two different orientations of FGF2 on the membrane surface (Figure 13). The high-affinity orientation (Figure 13A and Video 1) observed in 3 out of 5 repeats is characterized by strong binding to the membrane surface. In this case, FGF2 orients in a manner where C95 is exposed in a way facilitating the formation of a C95–C95 disulfide bridge with a second monomer. In this orientation, C77 is not available to make a contact with a second FGF2 monomer. In this high-affinity orientation, all the three known binding site residues of FGF2 are involved in the interaction with PI(4,5)P2 (Figure 13, Figure 13—figure supplements 2 and 3). The average binding energy for PI(4,5)P2 is highest for K127 (–315.59 kJ/mol), followed by K133 (–301.25 kJ/mol) and R128 (–208.78 kJ/mol). However, the interaction strengths were observed to depend on the number of PI(4,5)P2 head groups interacting with the binding pocket region. In addition to the three key binding residues, a number of so far unidentified residues contributed to membrane binding: K34 (average interaction strength –113.65 kJ/mol), K137 (–244.72 kJ/mol) and K143 (–196.29 kJ/mol) (Figure 13 and Figure 13—figure supplement 4). The number of head groups of PI(4,5)P2 that were found bound to FGF2 varied between one and four with one to two head groups binding to the key residues in the primary binding site (K127, R128, K133) and one to two head groups bound to residues forming the second binding site (K34, K137, K143).

In the second low-affinity orientation (Figure 13B) that was observed in 2 out of 5 simulation repeats, the N- and C- terminal ends of FGF2 are close to the membrane surface, and both C95 and C77 are available to form disulfide bridges with other FGF2 monomers (Figure 13B). However, in this scenario, not all of the three key binding site residues interact with PI(4,5)P2 as they point away from the PI(4,5)P2 head groups. This is also evident from the average binding energies (K127 = –0.03 kJ/mol; R128 = –237.10 kJ/mol and K133 = –0.08 kJ/mol) (Figure 13, Figure 13—figure supplements 5 and 6), which are considerably weaker compared to the high-affinity orientation (Figure 13 and Figure 13—figure supplement 7). Therefore, the weak interaction observed in the low-affinity orientation is likely to represent a transient binding intermediate at the membrane surface that is eventually stabilized when several PI(4,5)P2 molecules are bound to FGF2.

Based on the above, the stability of FGF2-membrane binding is based partially on the FGF2 orientation, but also the number of PI(4,5)P2 molecules that simultaneously bind to a FGF2 monomer is critical for stable membrane interactions of FGF2. This view is supported by additional simulations (system M2), where FGF2 was allowed to interact with only a single PI(4,5)P2 molecule (see Figure 13, Video 2, Figure 13—figure supplements 1B and 8). As expected, FGF2 was observed to bind to PI(4,5)P2 with its key binding site residues, however during the course of the simulation the orientation of FGF2 fluctuated between high-affinity and low-affinity orientations. None of the two orientations was exceptionally stable as the inositol biphosphate ring (central in FGF2 binding) was too flexible to stabilize the structure of the complex. Thus, the simulations indicate that the surface area of FGF2 involved in binding to PI(4,5)P2 is larger than previously assumed and renders interactions with more than one PI(4,5)P2 molecule possible. This leads to strengthening the stability of FGF2 membrane binding suggesting that the high-affinity orientation is the most stable one.

Video 2
FGF2 binding to a single PI(4,5)P2 (system M2).

FGF2 adhering to the membrane surface through binding to a single PI(4,5)P2. In the binding process, FGF2 first interacts with the residues close to the main binding site (K34, K137, K143, shown in purple van der Waals spheres), and then, at about 500 ns, rotates itself and interacts with the residues in the binding pocket (K127, R128, K133, shown in blue van der Waals spheres). The color coding used for the protein and PI(4,5)P2 is consistent with Figure 13. In the membrane, cholesterols are shown in light purple, the phosphorous atoms of POPC molecules are depicted as a transparent surface, and for clarity’s sake water and ions are not shown.

Legends for source data

https://doi.org/10.7554/eLife.28985.034

A model for PI(4,5)P2 dependent FGF2 oligomerization based on C95-C95 disulfide linked dimers and C77-C77 disulfide bridges involved in the formation of higher FGF2 oligomers

The critical and most likely rate-limiting step of FGF2 oligomerization is PI(4,5)P2 dependent dimerization. To identify possible dimerization interfaces along with the residues being involved, we conducted atomistic simulations (Figure 14). In previous FGF2 secretion assays and biochemical reconstitution experiments, FGF2 variant forms lacking both C95 and C77 were found defective in PI(4,5)P2 dependent oligomerization, membrane pore formation and secretion from cells (La Venuta et al., 2015; Müller et al., 2015). When single cysteine substitutions were analyzed, the C95A variant form was characterized by a more severe defect in all of these assays compared to a C77A variant form of FGF2 (Müller, Wegehingel, Steringer and Nickel, unpublished results). To identify the residues involved in dimerization, we carried out two simulations of FGF2 trimers (system T) since they correspond to the minimal aggregation unit we identified in single particle brightness analyses using supported lipid bilayers (Figures 10 and 11). To this end, FGF2 monomers were placed in different orientations at a distance of 0.5 nm above the PI(4,5)P2 head groups. The monomers were arranged such that two of them faced each other at a distance of 0.7 nm between their C95 residues. The third monomer resided 1.5 nm away from the other two monomers (Figure 14 and Figure 14—figure supplement 1). In the first simulation, the two monomers readily oriented themselves into the high-affinity orientation and then dimerized. The dimer was observed to remain stable for the rest of the simulation. A detailed analysis revealed that there were four ion pairs formed across the interface (two pairs of E86 – K118 and E99 – K85 each; Figure 14B) that stabilized the dimer. At the same time, the third monomer diffused towards the dimer structure, however its orientation was closer to the low-affinity orientation than the high-affinity counterpart. Given this, it was quite expected that it bound to the dimer through an orientation where its C77 residue faced C77 in one of the FGF2 monomers in the dimer structure (Figure 14C). The binding across this C77 – C77 interface was based on the R47 – D49, D45 – R41, and R80 – E53 ion pairs that rendered the formation of this complex possible. Moving on, the second simulation of system T also highlighted that FGF2 monomers tend to aggregate, thus supporting the view that there is a force driving FGF2 monomers to form oligomeric structures. However, the time scale needed for orientation changes through rotational diffusion turned out to be long, thus no stable interface for dimerization was found during this type of simulation.

Biochemical and structural experiments have shown that disulfide bridges from C95 – C95 and C77 – C77 residue pairs are present in membrane inserted FGF2 oligomers ([La Venuta et al., 2015; Müller et al., 2015]; Müller and Nickel, unpublished results). Consistently, in our simulations, we observed formation of FGF2 dimers characterized by the C95 – C95 and C77 – C77 pairs at the dimerization interface along with electrostatic interactions by additional residues in the dimer interface. Consequently, we used the dimer structure based on high-affinity FGF2 orientations with a stable C95 – C95 disulfide bridge in further simulations, where this link was established as a covalent bond to create a model of a full dimer (system D). The high-affinity orientation seen in simulations of system T was maintained during the simulations of system D. Importantly, with the dimer structure now fully stable due to the C95 – C95 covalent bond, the simulation of the dimer system D revealed several ion pairs (R80-D98; E86-K118 and R89-E99) (Figure 14, Figure 14—figure supplements 2 and 3). The extensive ion pairing formed at the interface as part of the high-affinity orientation suggests that the formation of FGF2 dimers containing the C95 – C95 disulfide bridge represents the initial step of PI(4,5)P2 dependent FGF2 oligomerization and membrane pore formation.

Discussion

The current study is the first of its kind in which an unconventional mechanism of protein secretion has been reconstituted from purified components establishing the molecular mechanism by which FGF2 is secreted from cells. We define the minimal machinery required for FGF2 membrane translocation and provide novel insights into the structure function relationship of membrane inserted FGF2 oligomers, the key intermediates of this process. Finally, using atomistic molecular dynamics simulations, this study puts forward a mechanism by which FGF2 monomers assemble into dimers and trimers on membrane surfaces in a PI(4,5)P2 dependent manner, the rate limiting step for the formation of higher oligomers that form membrane pores.

The first part of this work provides direct proof for two critical predictions of a previously proposed model describing the molecular mechanism of unconventional secretion of FGF2 from cells (La Venuta et al., 2015; Nickel, 2011; Nickel and Rabouille, 2009, 2008). The first predicted binding of FGF2 to PI(4,5)P2 versus cell surface heparan sulfates to be mutually exclusive to ensure directional translocation of FGF2 across the plasma membrane based on sequential interactions of FGF2 with PI(4,5)P2 at the inner leaflet and heparan sulfates at the outer leaflet. Here, using NMR spectroscopy we demonstrate an overlap of the binding sites on the molecular surface of FGF2 for IP3 [the headgroup of PI(4,5)P2] and a defined low affinity heparin disaccharide. In addition, we provide direct biochemical proof that long-chain heparins (used as mimetics of cell surface heparan sulfates) directly compete with membrane incorporated PI(4,5)P2 for the interaction with FGF2. Based upon the high affinity of FGF2 towards heparan sulfates and long-chain heparins [KD ≈ 100 nM (Faham et al., 1996)] and the lower affinity towards PI(4,5)P2 in the low micromolar range (Temmerman et al., 2008; Temmerman and Nickel, 2009) (Figure 1), our findings provide a direct explanation for vectorial translocation of FGF2 from the cytoplasm to the cell surface.

The second prediction of the FGF2 secretion model proposed the minimal machinery required for FGF2 membrane translocation to be composed of two trans-acting factors, PI(4,5)P2 and cell surface heparan sulfates. This in turn suggests a direct requirement for cis elements in FGF2 that mediate binding to PI(4,5)P2 (K127/R128) and heparan sulfates (K133). In addition, the core mechanism involves FGF2 oligomerization and membrane pore formation, a process that depends on two cysteine residues on the molecular surface of FGF2 (C77/C95). Using GUVs with a plasma membrane like lipid composition along with a FGF2-GFP fusion protein, we demonstrate FGF2 membrane translocation into the lumen of GUVs to depend on the presence of PI(4,5)P2 on membrane surfaces and the presence of long-chain heparins in the lumen of GUVs. Beyond these trans-acting factors, we demonstrate cis-elements in FGF2 required for binding to PI(4,5)P2 and heparin as well as for FGF2 oligomerization and membrane pore formation to be essential for FGF2 membrane translocation in a fully reconstituted system. By contrast, the trans-acting factor ATP1A1, even though required for FGF2 secretion from cells (La Venuta et al., 2015; Zacherl et al., 2015), apparently is dispensable for FGF2 membrane translocation under the conditions used in this study. This indicates that ATP1A1 has regulatory functions in FGF2 secretion from cells rather than belonging to the core machinery mediating physical translocation of FGF2 across the plasma membrane.

In the second part of this study, we used single molecule techniques to gain insight into the structure function relationship of membrane inserted FGF2 oligomers as key intermediates in membrane translocation. First, we used supported lipid bilayers (SLBs) and single molecule confocal microscopy to determine the oligomeric size distribution of membrane-inserted FGF2 species. Following short incubation times of 10 min during which the overall integrity of SLBs was fully maintained, three major oligomeric species could be detected represented by trimers, hexamers and higher oligomers in the range of about 10 to 12 subunits. These findings were corroborated by an independent single molecule system in which the average oligomeric state of FGF2 could be correlated with membrane insertion and pore formation in GUVs. Here, even after incubation times of several hours at high concentrations of FGF2, the average oligomeric state of FGF2 in GUVs leveled out in the range of 8 to 12 subunits. On the one hand, these experiments demonstrate membrane-inserted FGF2 oligomers to represent highly dynamic structures that are not characterized by a uniform number of subunits. On the other hand, our findings demonstrate that PI(4,5)P2 dependent membrane recruitment does not result in the random formation of undefined FGF2 aggregates. Therefore, as discussed above, we propose that membrane inserted FGF2 oligomers serve as dynamic translocation intermediates by controlled assembly through PI(4,5)P2 dependent FGF2 oligomerization at the inner leaflet and controlled disassembly at the outer leaflet mediated by heparan sulfates. Previous evidence suggests that membrane inserted FGF2 oligomers are accommodated within a lipidic membrane pore with a toroidal architecture (La Venuta et al., 2015; Steringer et al., 2012). We propose that these translocation intermediates represent dynamic structures to which FGF2 subunits are constantly added at the cytoplasmic leaflet while FGF2 subunits are continuously removed at the extracellular side of the plasma membrane resulting in FGF2 translocation to the cell surface. It is currently unknown whether these assembly/disassembly units are monomers or disulfide bridged dimers of FGF2. Furthermore, it will be a challenge for future studies to elucidate the three-dimensional architecture of membrane inserted FGF2 oligomers to fully understand how these complexes are functioning as translocation intermediates in unconventional secretion of FGF2.

The third part of this study was concerned with the analysis of the initial events of FGF2 membrane translocation with a particular focus on the molecular mechanism of FGF2 interactions with PI(4,5)P2 and the formation of dimers and trimers that initiate oligomerization of FGF2 in the context of PI(4,5)P2 containing membranes. Atomistic simulations revealed the key binding site residues that interact with PI(4,5)P2, which were found consistent with previous mutational analyses of FGF2 (6,17,19,25,33). Interestingly, the simulations predicted that there are additional previously unknown residues involved in membrane binding when several PI(4,5)P2 head groups bind to FGF2 simultaneously. Therefore, the strength of FGF2 membrane interactions obviously depends on the number of PI(4,5)P2 molecules that simultaenously bind to FGF2. In addition, the simulations provided compelling evidence that the strength of FGF2 membrane interactions also depends on the orientation of FGF2 relative to the membrane. The identified high-affinity orientation of FGF2 indeed favors multiple interactions of FGF2 with several PI(4,5)P2 head groups and also leads to the formation of stable FGF2 dimers that are characterized by C95 – C95 disulfide bridges and additional ion pair interactions between residues within the dimerization interface. By contrast, the low-affinity orientation of FGF2 tilts the protein in a way that impairs simultaneous interactions of FGF2 with several PI(4,5)P2 molecules and also prevents the formation of the C95 – C95 bridge. Given the critical role of C95 in FGF2 oligomerization, membrane pore formation, and secretion from cells, the simulations stress that FGF2 dimerization depends on membrane bound FGF2 monomers in the high-affinity orientation, whose stability is provided by multiple interactions with several PI(4,5)P2 molecules. Therefore, the combined findings from these simulations predict that PI(4,5)P2 dependent oligomerization of FGF2 is initiated through dimerization that is driven by the formation of C95 – C95 disulfide bridges. These FGF2 dimers appear to assemble into higher FGF2 oligomers driven by C77 – C77 disulfide links.

In conclusion, this study defines the minimal machinery required for FGF2 membrane translocation providing direct proof for a previously suggested model of unconventional secretion that was derived from cell-based data. This process depends on sequential and mutually exclusive interactions of FGF2 with PI(4,5)P2 and heparan sulfates and is thermodynamically driven by FGF2 oligomerization and membrane insertion. These translocation intermediates are dynamic structures with a subunit number in the range of 8 to 12 FGF2 molecules. Therefore, depending only on two essential trans-acting factors, PI(4,5)P2 and heparan sulfates, unconventional secretion of FGF2 is based upon a novel type of protein translocation across membranes with the cargo protein forming its own translocation intermediate by oligomerization and membrane insertion.

Materials and methods

Protein expression and purification

Request a detailed protocol

His-tagged variants of FGF2 (pQE30), FGF2-GFP and FGF2-Halo (both pET15b) were expressed in E. coli strains W3110Z1 or BL21 Star (DE3), respectively. For incorporation of the unnatural amino acid p-carboxylmethylphenylalanine (pCMF; custom synthesis by ENAMINE Ltd., Kiev, Ukraine), codon 81 (tyrosine) was replaced by an amber stop codon. Transformation of a strain carrying the pEVOL-pCMF plasmid resulted in expression of recombinant FGF2-Y81pCMF (Young et al., 2010). All proteins were purified in three steps via Ni-NTA affinity chromatography, heparin chromatography (except K127Q/R128Q and K127Q/R128Q/K133Q FGF2 variant forms) and size exclusion chromatography using a Superdex 75 column. In case of FGF2-Halo fusion proteins, desalting was performed using Nap-5 columns (GE Healthcare, Chicago, IL).

NMR spectroscopy

Request a detailed protocol

Recombinant 15N-labeled His-tagged FGF2-C77/95S (pQE30) was expressed in Escherichia coli W3110Z1 cells using M9 minimal medium with 15NH4Cl as the sole nitrogen source. Purification was performed as described above. Purified FGF2-C77/95S was diluted in 25 mM HEPES buffer (pH 7.4), containing 150 mM KCl and 10% D2O, to a final concentration of 160 µM. For the IP3 and heparin disaccharide titration experiments, 500 µl of 80 µM FGF2-C77/95S were titrated with defined volumes of 80 µM FGF2-C77/95S and either 900 µM IP3 (Sigma 74148, Sigma Aldrich, St. Louis, MO) or heparin disaccharide (Sigma H9267). The endpoint of the IP3 titration was further titrated with defined volumes of 80 µM FGF2-C77/95S and 900 µM heparin disaccharide and vice versa. NMR spectra were recorded on a Bruker AV 700 MHz NMR spectrometer equipped with a 5 mm triple resonance cryo-probe at 300 K. For each titration step 2D 15N1H-HSQC spectra were recorded with 1024 points in the 1H dimension and 96 points in the 15N dimension and averaged over eight transients. Spectra were processed with TopSpin (Bruker, Billerica, MA) using CcpNmr Analysis software (Vranken et al., 2005). Signals were assigned using previously published data (accession code: 1BLA; (Moy et al., 1996, 1995)) Chemical shift differences were calculated using the equation Δδ=(δH)2+(0,15δN)2.

Binding of FGF2 to PI(4,5)P2-containing liposomes in the presence and absence of heparin

Request a detailed protocol

Large unilamellar vesicles (LUVs) with a plasma-membrane-like lipid composition either lacking (PM) or containing 2 mol% PI(4,5)P2 (PM +PIP2) were prepared in buffer A (25 mM HEPES, pH 7.4, 150 mM KCl) supplemented with 10% (w/v) sucrose as described previously (Steringer et al., 2012). After blocking LUVs with 3% (w/v) fatty-acid free BSA in buffer A for 1 hr at 25°C, membranes were washed with buffer A and collected by sedimentation (15000x g; 20°C, 10 min). LUVs were resuspended in buffer A containing 2.5 µM His-FGF2-Y81pCMF-WT. After 1 hr of incubation at 25°C, either a mixture of long-chain heparins (Sigma H3149) or a defined heparin disaccharide (Sigma H9267) were added. In case of long-chain heparins, molar concentrations of FGF2 binding sites were defined by heparin disaccharide units with four sodium ions bound (MW 685). Following incubation at 25°C, unbound His-FGF2-Y81pCMF-WT and LUVs with bound His-FGF2-Y81pCMF-WT were separated by centrifugation (15000x g; 20°C, 10 min). While the supernatant (unbound material) was mixed with SDS-sample buffer, the sediment (bound material) was washed with buffer A followed by sedimentation of liposomes. The final pellet was resuspended in SDS-sample buffer. Samples (50% of bound and 14% of unbound material as well as 15% of input material) were analysed on 4–12% Bis-Tris SDS-PAGE in MES buffer (NuPAGE, Thermo Fisher Scientific, Waltham, MA). Proteins were stained with Coomassie Instant Blue (Expedeon, UK) followed by quantification and data normalization using the LI-COR Odyssey infrared imaging platform and Image Studio Lite Software (Version 5.2.5) (LICOR Biosciences, Lincoln, NE). FGF2 binding efficiency to PM-like liposomes containing PI(4,5)P2 in the absence of heparin was set to 100% (bound material). Background binding was defined by PM-like liposomes lacking PI(4,5)P2 (100% unbound material). Mean values with standard deviations (SD) are shown (n = 3).

Characterization of FGF2 variant forms with regard to binding to heparin and PI(4,5)P2

Request a detailed protocol

Heparin Sepharose 6 Fast Flow beads (20 µl; GE Healthcare) were equilibrated in 25 mM HEPES pH 7.4, 150 mM NaCl (buffer B). Heparin beads were incubated with the FGF2 variant forms indicated (2.5 µM protein concentration in 200 µl of buffer B) for 1 hr at room temperature on a rotating wheel. Afterwards, beads were pelleted (500 g; 25°C, 3 min) and the corresponding supernatants were treated with SDS sample buffer (unbound material). Heparin beads with bound proteins were washed three times with buffer B. Bound proteins were eluted with SDS-sample buffer. Input (5%), bound (20%), and unbound material (5%) were analyzed using 4–12% Bis-Tris acrylamide gels with MES buffer (NuPAGE, Thermo Fisher Scientific, Waltham, MA). Coomassie Instant Blue (Expedeon) stained bands were quantified on the LI-COR infrared imaging platform using Image Studio Lite Software (version 5.2.5). Mean values with standard deviations (SD) calculated from three independent experiments are shown.

FGF2 binding to PI(4,5)P2 was quantified using flow cytometry as described previously (Temmerman et al., 2008; Temmerman and Nickel, 2009). PM-like LUVs (containing 1 mol% Rhodamine-PE) with or without 2 mol% PI(4,5)P2 were blocked with 3% (w/v) fatty-acid free BSA in buffer A (25 mM HEPES, pH 7.4, 150 mM KCl) for 1 hr at 25°C. LUVs were sedimented, washed with buffer A and collected by centrifugation (16000x g; 25°C; 10 min). Following resuspension, LUVs were incubated with the various recombinant FGF2-GFP fusion proteins indicated at a final concentration of 2 µM in buffer A. Bound and unbound proteins were separated by sedimentation of LUVs followed by extensive washing in buffer A. The final liposome pellet was resuspended in 300 µl buffer A and analysed by flow cytometry using a FACS Calibur system (BD Biosciences, San Jose, CA). In brief, liposomes were gated by size and rhodamine fluorescence based upon the Rhodamine-PE membrane tracer. Simultaneously, GFP signals were measured to quantify binding of the various FGF2-GFP variant forms. To quantify protein-lipid interactions, raw data (fluorescent units) were corrected for liposome tethering using a shape index as reported previously (Temmerman et al., 2008; Temmerman and Nickel, 2009). Data were normalized defining GFP signals obtained with PM + PIP2 liposomes and FGF2-Y81pCMF-GFP-WT as 100% binding efficiency. Standard deviations are shown (n = 4).

Preparation of giant unilamellar vesicles (GUVs)

Request a detailed protocol

GUVs with a plasma membrane like lipid composition consisting of 30 mol% cholesterol (Chol), 15 mol% sphingomyelin (SM), 34 mol% phosphatidylcholine (PC), 10 mol% phosphatidylethanolamine (PE), 5 mol% phosphatidylserine (PS), 5 mol% phosphatidylinositol (PI) and 1 mol% Biotinyl-PE (Avanti Polar Lipids, Alabaster, AL) were generated based on electro-swelling using platinum electrodes (García-Sáez et al., 2009). Where indicated, GUVs were supplemented with either PI(4,5)P2 or a Ni-NTA lipid at 2 mol % at the expense of PC. In some experiments, GUVs were used lacking cholesterol which was supplemented with PC. For visualization either 0.05 mol% rhodamine B-labelled PE for FGF2 membrane translocation assays or 0.001 mol% DOPE-Atto 633 for z-scan FCS were added. The dried lipid film was hydrated with a 300 mM sucrose solution (300 mOsmol/kg). Where indicated, either long-chain heparins (50 µM; based on disaccharide units s.a.) or a defined heparin disaccharide (Sigma H9267) were included to mimic heparan sulfates in the lumen of GUVs. Osmolality was determined for all three conditions using a Wescor Vapro 5600 instrument and found not to be affected by the addition of long-chain heparins or the heparin disaccharide. Swelling was conducted at 45°C [10 Hz, 1.5 V for 50 min (without heparin) or 70 min (with heparin), 2 Hz, 1.5 V for 25 min]. In order to remove excess amounts of heparin, GUVs were gently washed with buffer B (25 mM HEPES pH7.4, 150 mM NaCl, 310 mOsmol/kg) and collected via centrifugation (1200x g; 25°C; 5 min). The loose GUV pellet was carefully resuspended in a small volume of buffer B and diluted again in 11.5 ml buffer B followed by centrifugation (1200x g; 25°C; 5 min). The supernatant was removed while the loose GUV pellet was carefully resuspended. Imaging chambers (LabTek) were incubated sequentially with 0.1 mg/ml Biotin-BSA (Sigma A8549) and 0.1 mg/ml Neutravidin (Thermo Fisher Scientific A2666) in buffer B. Luminal incorporation of heparin into GUVs was monitored by confocal microscopy in control experiments using a fluorescent derivative (Molecular probes; H7482). Likewise, the presence of PI(4,5)P2 in the bilayers of GUVs was analyzed using a recombinant fusion protein of GFP with the PH domain of PLCδ1, a canonical PI(4,5)P2 marker (Lemmon, 2003; Milosevic et al., 2005).

Imaging and quantification of FGF2 membrane translocation using GUVs

Request a detailed protocol

For FGF2 membrane translocation assays, GUVs were incubated for 3 hr with a small fluorescent tracer (Alexa647) and the FGF2-GFP fusion proteins indicated at a final concentration of 200 nM. Confocal images were recorded at room temperature in multitrack mode using Zeiss LSM510 and LSM780 confocal fluorescence microscopes (Carl Zeiss AG, Oberkochen, Germany) along with a plan apochromat 63 x/1.4 oil immersion objective. Pinholes of the tracks were optimized to 1.2 µm. LSM510: In order to measure (i) GFP-, (ii) Rhodamine-PE-, and (iii) Alexa647-derived signals, samples were excited with (i) an argon laser (488 nm), (ii) a He-Ne-laser (561 nm), or (iii) a He-Ne laser (633 nm) and light was detected after (i) a band pass (BP) filter (505–530 nm), (ii) a BP filter (560–615 nm), or (iii) a long pass (LP) filter (>650 nm). Images were recorded in 8-bit grayscale. LSM780: Samples were excited with (i) an argon laser (488 nm) for GFP, (ii) a He-Ne-laser (561 nm) for Rhodamine-PE, or (iii) a He-Ne laser (633 nm) for Alexa647 and light was split by a spectral beam guide to (i) 498–550 nm, (ii) 585–673, or (iii) 654–759 nm. Images were recorded in 12-bit grayscale and are shown pseudo-coloured in (i) green (Track2-ChS1), (ii) red (Track1-ChS2) and (iii) gray (Track2-Ch2). A radial profile analysis was conducted using ImageJ software (http://rsbweb.nih.gov/ij/) using the plugin ‘radial profile’ (http://rsbweb.nih.gov/ij/plugins/radialprofile.html) to detect intensity differences in the lumen, at the membrane, and in the surroundings of GUVs. Briefly, a circle was drawn around a GUV in a way that the center of the circle matched the center of the GUV in the confocal plane being assessed. From this center point the intensity at any given distance along the radius of the circle is measured and processed resulting in a profile plot of normalized integrated intensities around concentric circles as a function of distance from the center. Luminal fluorescence of individual GUVs was measured and normalized to fluorescence intensity of the surrounding buffer. Per experimental condition, 20 to 120 individual GUVs were analysed as indicated in the corresponding Fig. legends. To allow for a statistical analysis of membrane pore formation and FGF2-GFP membrane translocation across the population of GUVs, thresholds were defined to classify individual GUVs. When the ratio of inside to outside fluorescence of the small Alexa647 tracer was ≥0.6, GUVs were classified as vesicles containing membrane pores. Similarly, when the inside to outside ratio of GFP fluorescence was ≥1.6, the corresponding GUVs were classified as vesicles where FGF2-GFP membrane translocation into the lumen had occurred.

3D reconstruction of GUVs

Request a detailed protocol

3D reconstruction of GUVs (panels C and D of Figure 3) was based upon ImageJ software using the ‘3D Viewer’ plugin (https://imagej.nih.gov/ij/plugins/3d-viewer/index.html). 3D images with a resolution of 512 × 512 pixels (0.07 µm pixel size) were made from z-stacks separated by a distance of 0.37 µm. All three channels – green (FGF2-Y81pCMF-GFP), red (Rhodamine-PE) and magenta (Alexa647) were merged into a single 3D image. Cross sections of GUVs are shown to visualize FGF2-Y81pCMF-GFP membrane translocation into the lumen.

Analysis of FGF2 oligomerization on supported lipid bilayers

Request a detailed protocol

Protein labeling was conducted by incubating FGF2-Y81pCMF-Halo and FGF2-Y81pCMF-C77/95A-Halo, respectively, (20 µM each) with 40 µM StarRed-HTL in 150 mM NaCl, 20 mM HEPES (pH 7.4) for 1 hr at room temperature. Unbound dye was removed by size exclusion chromatography using a Nap-5 column (GE Healthcare). Degree of labelling (average number of dyes per FGF2 monomer) was determined to be DOL = 0.9 for both FGF2-Y81pCMF-Halo and FGF2-Y81pCMF-C77/95A-Halo, by measuring absorption spectra of the labelled protein and calculating the DOL according to:

DOL= A638*εDye(A280A280*CF)*εProt 

with A being the local absorption maxima, εDye=120000, CF=0.32, εProt=79870. Supported lipid bilayers (SLBs) were prepared on 18 mm round cover glasses. Glass surfaces were cleaned in 4% Hellmanex II in a bath sonicator for 10 min, rinsed and dried under nitrogen. Right before usage, cover glasses were plasma cleaned for 1 min. SLBs were formed by spin-coating at 3000 rpm for 1 min using 40 µl of lipid mixture [1.2 mg/ml; 68 mol% POPC, 30 mol% Cholesterol, 2 mol% PI(4,5)P2 and 0.1 mol% DPPE-OregonGreen in methanol/chloroform (1:1)]. Residual solvent was evaporated for 1 hr under vacuum. Lipids were hydrated with a buffer containing 150 mM NaCl and 20 mM HEPES (pH 7.4). This procedure resulted in supported bilayers without defects. Occasionally, patches of double membranes were observed and avoided when studying FGF2 binding and oligomerization. FGF2-Y81pCMF-HALO-StarRed and FGF2-Y81pCMF-C77/95A-HALO-StarRed, respectively, were added to SLBs (final concentration = 100 nM) and incubated for 10 min. Membranes were washed in two steps with 150 mM NaCl buffer to remove unbound proteins and 500 mM NaCl to remove non-oligomerized FGF2 from membranes.

FGF2-HALO-StarRed fusion proteins associated with SLBs were imaged using a confocal laser scanning microscope with a 775 nm super-resolution STED module (Abberior Instruments, Göttingen, Germany). Fluorescence recovery after photo bleaching was performed by bleaching a 2 × 2 µm square with a 405 nm laser followed by measuring the diffusive recovery of FGF2 fusion proteins over time. Diffusion constant and mobile fraction of FGF2 was estimated by fitting the recovery curves to a single exponential in MATLAB.

Single molecule imaging was done using 10 µW laser power (640 nm diode, 40 mMHz) in the backfocal plane of an 100x Oil NA 1.4 (Olympus) objective. Pixel size was 80 nm and pixel dwell-time 1 ms. To reduce background, we removed the first 0.5 ns of the detected photons by time-gating. This effectively reduced the background from scattered laser light. The resulting images were processed in MATLAB to extract the peak brightness per cluster using a 2D Gaussian fitting procedure. The number of monomers in each immobilized FGF2 cluster was estimated by normalizing the cluster brightness to the brightness of a single HALO-StarRed. To determine the brightness of a single HALO-StarRed, we expressed, purified and labelled mEGFP-HALO with StarRed-HTL. mEGFP-HALO-StarRed was then passively immobilized on clean cover glass such that single diffraction limited spots could be clearly separated. Imaging and single molecule brightness determination was identical with the procedure used for FGF2-HALO-StarRed fusion proteins. A Gaussian mixture analysis in MATLAB was conducted to estimate the number of sub-populations in the distribution of single molecule brightness and oligomeric number.

Determination of average oligomeric states of GUV associated FGF2-Y81pCMF and FGF2-Y81pCMF-C77/95A-GFP by z-scan FCS

Request a detailed protocol

Measurements were performed on a self-constructed confocal microscope consisting of an inverted confocal microscope body IX71 (Olympus, Hamburg, Germany) and pulsed diode lasers (LDH-P-C-470, 470 nm, PicoTA, 532 nm, and LDH-D-C-635, 635 nm, PicoQuant, Berlin, Germany). The lasers were pulsing alternately at the repetition rate of 12.5 MHz to avoid artifacts caused by signal bleed-through. The laser light was coupled to a polarization maintaining single mode optical fiber and re-collimated at the output with an air space objective (UPLSAPO 4X, Olympus). The light was up-reflected to a water immersion objective (UPLSAPO 60x, Olympus) with a 465/533/635 dichroic mirror. The signal was split between two single photon avalanche diodes using 515/50 (FGF2-GFP) and 697/58 (DOPE-Atto 633) for FCS measurements or 515/50 (GFP) and 595/50 (small tracer AlexaFluor532) band pass filters (Chroma Rockingham, VT) to analyze membrane pore formation in GUVs selected for a FCS measurement. The laser intensity at the back aperture of the objective was kept below 10 µW for each laser line.

The determination of oligomeric states was based on a comparison of the brightness of a cluster ϕcluster to that one of a monomer ϕmonomer. GUVs were selected and classified with regard to the presence (small tracer AlexaFlour532 equilibrated between lumen and exterior) or absence (no tracer in the vesicle interior) of membrane pores. In the next step, a membrane was placed into the waist of 470 and 635 nm lasers, moved 1.5 µm below the waist and scanned vertically along the z axis in 20 steps spaced 150 nm from each other. FCS measurements (60 s each) using 515/50 nm (FGF2-GFP) and 697/58 nm (DOPE-Atto 633) emission channels were performed at each step. Auto-correlation (AC) curves corresponding to each of the 20 steps were fitted by a model assuming 2D diffusion within the membrane (bound FGF2-GFP/DOPE-Atto-633), 3D diffusion in solution (free FGF2-GFP) and transition of a dye to the triplet state (Widengren et al., 1994):

(1) G(τ)=1+(1PN2D11+(τ/τ2D)+1PN3D11+(τ/τ3D)1+SP(τ/τ3D))1T+T exp(τ/τT)1T

In Equation (1), τ is the lag-time, PNx the particle number, τiD the dye diffusion time, T the fraction of the dye in the triplet state and τT the lifetime of the triplet state. While auto correlation curves belonging to DOPE-Atto 633 were used to judge on quality of the lipid bilayer, those belonging to FGF2-GFP fusion proteins were used for the brightness analysis. Suitable properties of the membrane were indicated by free diffusion of DOPE-Atto 633 in the bilayer and a reasonable diffusion coefficient obtained for this probe. When plotting PN2D(FGF2-GFP) against the vertical position z, a parabolic dependence is obtained due to the Gaussian shape of the excitation beam profile (Benda et al., 2003). The minimum value of PN2D(z, FGF2-GFP), min[PN2D(z, FGF2-GFP)], obtained from such a plot and the corresponding average intensity in counts per second were used for further analysis. Position z corresponds to the middle of the membrane and allows for direct comparison of brightness values obtained for different measurements. Moreover, contribution from the bulk to the overall fluorescence signal is negligible at this position (1/PN3D0). The brightness of a cluster was calculated as ϕcluster=I/min[ PN2D(z,HisFGF2Y81pCMFGFP) ]. The brightness of a monomer ϕmonomer was obtained in a similar manner as ϕcluster, however, the probability that two labelled FGF2 molecules meet randomly in a cluster must be negligible. This has been achieved by (i) using the recombinant protein HisFGF2-Y81pCMF-C77/95A-GFP which binds to DGS-NTA containing bilayer as a dimer at maximum and (ii) diluting HisFGF2-Y81pCMF-C77/95A-GFP by the unlabeled variant HisFGF2-Y81pCMF-C77/95A. The average oligomeric state is calculated as

(2) oligomeric state=ϕclusterϕmonomer

The diffusion coefficient was obtained by plotting τ2D against z (z-scan FCS approach). In analogy to PN2D(z), a parabolic dependence was obtained. D was determined by using the following relationship (Benda et al., 2003):

(3) D=w024min[ τ2D(z) ]

with w0 being the radius of the beam waist.

The protein surface concentration [c(FGF2-GFP)] was determined in two different ways. The first method makes use of the readouts that are obtained during the process of determining the oligomeric state of FGF2:

(4) c(FGF2)=oligomeric state*PN2D(FGF2)πw02

The second method is based on the fact that fluorescence intensity is directly proportional to concentration. Therefore, by having a calibration solution with a given constant concentration of a dye, the actual concentration of FGF2-GFP at the membrane can be determined by the analysis of GUV fluorescence images. Both methods yielded similar results. The data shown in Figure 12 were derived using the first method.

Atomistic molecular dynamics simulations analyzing PI(4,5)P2 dependent oligomerization of FGF2

Request a detailed protocol

For atomistic molecular dynamics (MD) simulations, we used the truncated FGF2 monomer structure [PDB id: 1BFF; (Kastrup et al., 1997)] from residue 26 to 154, without the flexible N-terminus. Cell based as well as biochemical experiments revealed that N-terminal truncations did not show defects in secretion nor oligomerization (André Engling, Julia Steringer, Hans-Michael Müller, Sebastian Unger and Walter Nickel, unpublished results). Here, we use the same residue numbering as in the PDB file (Kastrup et al., 1997), which is consistent with the numbering used in biochemical experiments. In preparation of the FGF2 structure for MD simulations, we processed the PDB structure with CHARMM-GUI (Lee et al., 2016), where we modeled all lysine (K) and arginine (R) residues as positively charged, glutamic acid (E) and aspartate (D) as negatively charged, histidine (H) and glutamine (Q) as uncharged, and tyrosine 81 was phosphorylated.

For the lipid bilayer models, we used a POPC/Cholesterol/PI(4,5)P2 mixture (details below). PI(4,5)P2 had 18:0 and 20:4 chains. The double bonds in the arachidonic fatty acid were located between the carbons C5 = C6, C8 = C9, C11 = C12, and C14 = C15. The overall charge in PI(4,5)P2 was –4e. The bilayer systems prepared using CHARMM-GUI (Lee et al., 2016) were originally symmetric with respect to the transmembrane distribution of POPC and Cholesterol. PI(4,5)P2 molecules were then inserted to the cytosolic leaflet of a membrane by removing a sufficient number of POPC molecules from the cytosolic leaflet and replacing the vacant space by PI(4,5)P2, thus resulting in systems M1 and M2 (Table 1). System M1 had a POPC/Cholesterol (70/30) mixture in the extracellular leaflet and a POPC/Cholesterol/PI(4,5)P2 (65/29.5/5.5) composition in the cytosolic monolayer. System M2 was designed to assess FGF2 binding in the dilute PI(4,5)P2 limit, having only one PI(4,5)P2 on the cytosolic side (Table 1).

Table 1

Descriptions of model systems and their simulation details.

https://doi.org/10.7554/eLife.28985.035
SystemNo. of FGF2 moleculesNo. of POPC moleculesNo. of CHOL moleculesNo. of PI(4,5)P2 moleculesNo. of water moleculesDistance from membrane surface (nm)No. of repeatsSimulation time (ns)
M1134815214347831.551000 × 5
M21177761156841.511000
T334815214315750.521000 × 2
D23481521429770*1130
  1. *The initial FGF2 dimer orientation and position on the membrane surface were the same as those in the final frame of the system T (repeat 2).

In setting up the systems M1 (for five simulation repeats) with a single FGF2 monomer (Table 1), we placed the protein 1.5 nm above the membrane surface in five different orientations such that the key residues (K127, R128, K133) in the binding pocket for PI(4,5)P2 either faced the PI(4,5)P2 head groups or pointed away from them. For system M2, we placed a single FGF2 protein 1.5 nm above the membrane surface with the protein’s binding pocket residues facing the PI(4,5)P2 head group.

The FGF2 trimer systems (T; see Table 1) simulated through two repeats were prepared by placing three FGF2 monomers 0.5 nm above the membrane surface with the binding pocket residues (K127, R128, K133) facing the PI(4,5)P2 head groups, yet each of the monomers was positioned in a slightly different orientation with respect to the membrane surface. Here, the monomers were arranged such that two monomers faced each other with a C95 – C95 distance of 0.7 nm, and the third monomer was 1.5 nm away from the two other monomers.

For the FGF2 dimer system (D; see Table 1), the starting dimer structure was taken from the final frame of the trimer simulation (repeat 2) that demonstrated a stable dimer characterized by a C95 – C95 close-contact bridge with an upright orientation. We processed the dimer by creating a disulfide linkage between the two C95 residues by removing the hydrogen atom from the thiol groups. The C95 – C95 disulfide-linked dimer structure was energy minimized in vacuum with the membrane, thus not altering either the FGF2 dimer orientation or its interaction with the lipids.

The systems were neutralized by an appropriate number of potassium atoms. KCl salt at a concentration of 150 mM was added to mimic experimental conditions. Protein, lipids, and salt ions were described using the CHARMM36 force field (Best et al., 2012; Klauda et al., 2010). For water, we used the TIP3 model (Jorgensen et al., 1983).

All the systems were subjected to energy minimization using the steepest descent algorithm. After minimization, we ran 6 steps of equilibration runs where we gradually reduced the force constant applied to restrain the positions of the protein and the lipids. In the first stage of equilibration under NVT conditions, the Berendsen thermostat (Berendsen et al., 1984) was used to regulate the temperature at 298 K with a time constant of 1.0 ps. In the second stage of equilibration under NpT conditions, we used the Berendsen algorithm for controlling the simulation temperature and pressure along with a semi-isotropic pressure coupling scheme and a time constant set to 5.0 ps. In equilibration under NpT conditions, the reference pressure was set to 1.0 bar and the isothermal compressibility to a value of 4.5 × 10–5 bar–1. For neighbor searching, we used the Verlet scheme with a cut-off distance of 1.2 nm for short-range neighbor search. The electrostatic interactions were calculated using the PME method (Darden et al., 1993; Essmann et al., 1995). The cut-off length of 1.2 nm was used for both electrostatic (real space component) and van der Waals interactions. Hydrogen bonds were constrained using the LINCS algorithm (Hess et al., 1997) and periodic boundary conditions were applied in all directions. For the production MD runs, we removed all the restraints applied to the proteins and the lipids and used the Nose-Hoover thermostat (Nosé, 1984; Hoover, 1985) and Parrinello-Rahman barostat (Parrinello and Rahman, 1981) instead of the Berendsen’s algorithm. The rest of the input parameters for production MD simulations were the same as those used under NpT conditions. The simulations were carried out using an integration time step of 2 fs with coordinates saved every 100 ps.

All simulations were performed using the GROMACS 5.1 (Abraham et al., 2015) simulation package and the analyses were done for the last 200 ns of the simulation trajectories using GROMACS tools and in-house scripts, except for system D, where the analysis was carried out over the last 30 ns. The images were rendered using VMD (Humphrey et al., 1996). The rendering was done for the starting structures and for the most populated structures obtained from clustering analysis based on RMSD values with a cut-off of 0.1 nm (Daura et al., 1999).

References

Decision letter

  1. Christopher G Burd
    Reviewing Editor; Yale School of Medicine, United States

In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.

[Editors’ note: a previous version of this study was rejected after peer review, but the authors submitted for reconsideration. The first decision letter after peer review is shown below.]

Thank you for submitting your work entitled "Minimal Machinery Required for Unconventional Secretion of Fibroblast Growth Factor 2" for consideration by eLife. Your article has been reviewed by three peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by a Senior Editor. The reviewers have opted to remain anonymous.

Our decision has been reached after consultation between the reviewers. Based on these discussions and the individual reviews below, we regret to inform you that your work will not be considered further for publication in eLife.

As you will see from the reviewer's comments, the biochemical reconstitution of FGF2 translocation is considered a biochemical tour de force. Discussion amongst the reviewers, however, resulted in a unanimous view that, at present, the manuscript largely confirms the molecular requirements for FGF2 translocation that were identified previously. The reviewers and I encourage you to continue this line of experimentation, and to this end, they make several suggestions that may help you to provide new information regarding the mechanism of FGF2 translocation. Specific questions that the reviewers felt were especially important to address and expand on regard the nature of the FGF2 translocation intermediate and PIP2 dynamics. Specifically, does FGF2 form a structural pore (versus simply rupturing the membrane) and, related, does PIP2 translocate with FGF2, or remain on the 'donor' leaflet? Though the reviewers concluded that suitable revisions to the manuscript would likely substantial effort, they expressed enthusiasm for considering a manuscript that does contain additional information.

Reviewer #1:

This is a superbly conducted study that investigates the translocation of FGF2 across a synthetic membrane bilayer using a GUV-based experimental system. The authors find that binding of FGF2-GFP to PIP2 is required to initiate translocation and that lumenal heparin is required to accumulate FGF2-GFP within the lumen of the vesicle. Long chain heparin is shown to effectively compete for binding to PIP2, providing an explanation for the vectoral transport (i.e., secretion) previously observed with cells. The results generally support the conclusions of the manuscript, though I think that the authors need to address a couple of points to raise the significance of the manuscript.

1) I commend the authors for a rigorous and convincing reconstitution of FGF translocation, though the manuscript essentially confirms requirements defined by published work using cell-based approaches without revealing anything fundamentally new about the mechanism of translocation. Two suggestions for enhancing the significance of the work are listed below.

2) The conclusion that FGF2-GFP forms a pore, revealed by diffusion of Alexa 647 across the GUV membrane, has important implications regarding the mechanism of translocation. The formation of a toroid intermediate was alluded to (subsection “Biochemical reconstitution of FGF2 membrane translocation with purified components”), but the results don't distinguish the formation of a bona fide pore, which in my mind is constituted by a structure of defined stoichiometry and dimensions, versus a less definite breach of membrane integrity (e.g., due to crowding on the membrane) that results in leakage of the dye across the membrane. I wonder if a systematic test of soluble tracers of varying size might reveal a size limit that would lend support to the conclusion that a pore is an intermediate structure formed during translocation.

3) The conclusion that ATP1A1 plays no core role in FGF2 translocation is not satisfying because it raises the question, does the observed translocation truly reconstitute the physiological mechanism? I realize that this is technically challenging to address, but defining the contribution of ATP1A1, be it a core role or a regulatory function, would, (a) provide new information, and/or (b) might bolster the conclusion that it is not required for translocation.

Reviewer #2:

This is an interesting report that demonstrates a possible mechanism or at least a minimal system allowing translocation of FGF2 across a lipid bilayer based on FGF2 first binding and clustering on one side of a PIP-containing membrane and then flipping over to be captured by heparin on the other side. The in vitro results demonstrating feasibility are generally clear, and the mutational analysis, especially of mutants with perturbed binding to PIP2 but not heparin are consistent with the proposed cellular mechanism. Several issues should be addressed to strength the claim that this effect is what happens in the cell.

What is missing is discussion of how FGF2 traverses the hydrophobic part of the membrane or if it always travels with PIP2 and delvers it to the other leaflet. There is extensive literature about PIP2 flipping or enhancing flipping of other lipids by e.g. Sulpice and colleagues in platelets and other systems that seems highly relevant and should be cited and discussed in this context. E.g. Bucki R, Giraud F, Sulpice (2000) Phosphatidylinositol 4,5-bisphosphate domain inducers promote phospholipid transverse redistribution in biological membranes. Biochemistry 39: 5838-5844.

Is the only binding of FGF2 to PIP2 via the headgroup? If not, the competition between IP3 and heparin disaccharide lacks the full binding energy and it's not clear how to distinguish binding from electrostatic competition.

How specific is the effect of heparin? Is there a similar effect of hyaluronic acid, which might be more commonly bound on the outer cell surface?

How much of the interaction in vitro is due to electrostatics? How dependent is the binding affinity on ionic strength or divalent cations?

A control is needed to determine if heparin in the vesicle lumen changes the osmotic pressure inside and consequently the membrane tension.

A most interesting test would be between two vesicle types that have the same amount of PIP2 with heparin in the lumen but with different cholesterol content to shift from fully mixed to demixed bilayers, since this process is implicated in trafficking. The other controls are required technical controls but have no direct biological relevance.

Reviewer #3:

In their study, the authors reconstitute membrane translocation of FGF2 in vitro, confirming roles for PI(4,5)P2, heparan sulfates, and FGF2 oligomerization in the process. Other factors that promote FGF2 plasma membrane translocation in vivo were not required, suggesting they function only in a regulatory manner. Overall, the study is beautifully performed, and the data presented are rigorous and quantitative. However, the major concern is that the work is largely confirmatory in nature. Please see below.

1) Sequential interactions with PI(4,5)P2 and heparan sulfate was shown previously to be critical in FGF2 translocation across the plasma membrane. Also, the role for PI(4,5)P2 in promoting FGF2 oligomerization and membrane pore formation was also defined previously. The work done here solidifies/confirms our understanding of FGF2 translocation, but does not really extend it. In light of this, can the authors provide any evidence as to the mechanism by which ATP1A1 regulates FGF2 translocation, even if it is not absolutely required?

2) To extend the study, could the author investigate whether the addition of appropriately positioned cysteine residues to other members of the FGF family members, which have signal peptides, enables them to be transported across membranes in vitro?

[Editors’ note: what now follows is the decision letter after the authors resubmitted for further consideration.]

Thank you for resubmitting your work entitled "Key Steps in Unconventional Secretion of Fibroblast Growth Factor 2 Reconstituted with Purified Components" for further consideration at eLife. Your revised article has been favorably evaluated by Vivek Malhotra (Senior editor) and three reviewers, one of whom is a member of our Board of Reviewing Editors.

The manuscript has been improved and the reviewers recommend publication. You will see from the reviews below, which include the opinion of a new reviewer, that there are several points that should be addressed with additional revisions. These require only modification to the text suggested by reviewer 1 (point #3), reviewer # 2 (figure and the legend) and 3.

Reviewer #1:

The major issue raised during review of the original manuscript regarded the confirmatory nature of the manuscript, which confirmed specific requirements that were identified previously using cell based assays. In the revised manuscript, the authors now include several new findings that shed some light onto the mechanism of FGF2 translocation across a membrane. I find the most compelling of these to be the following:

1) Calibrated fluorescence measurements of labeled FGF2 was used to quantify the numbers of FGF molecules present in membrane-associated clusters. This provides some insight into possible mechanism of formation of an initial translocation intermediate.

2) Molecular dynamics simulations of FGF2 binding to PIP2, and FGF2 oligomerization, revealed possible new binding modes. Although these were not validated with 'wet' experiments, they provide some evidence to support the authors' speculation that FGF2 assembles into a larger pore-forming complex that mediates translocation.

3) The new data lead the authors to suggest a toroid translocation intermediate. I have a difficult time appreciating the manner in which FGF subunits might form a toroid intermediate, and especially how such a disulfide-bonded structure would be released from the cell, but the new data provide a foundation for this speculation. As this is not addressed in a substantial manner in the Discussion, the authors may wish to describe their thinking in a bit more detail.

Overall, my major concerns raised with the original submission have been substantially addressed.

Reviewer #2:

Fibroblast Growth Factors (FGFs) participate in a cell signaling pathways that are required for essential biological functions such as embryonic development, hematopoiesis, and wound repair. The focus of this manuscript is FGF2 which functions in normal cell growth and when mis-regulated plays a role in tumor-induced angiogenesis that is thought to be associated with tumor cell resistance to treatment therapies. FGF2 lacks a signal peptide and is transported into extracellular space by an ER/Golgi-independent mechanism. Results from cell-based experiments and in vitro studies, translocation of FGF2 is proposed to happen by direct protein translocation across the plasma membrane. This process depends on FGF2 oligomerization and sequential interaction with PI(4,5)P2 at the inner leaflet and heparan sulfate proteoglycans at the outer surface.

In this report the authors use NMR, in vitro translocation studies, and molecular dynamic simulations to study the mechanism of FGF2 translocation across membranes. Using NMR titration and binding assays the authors show that FGF2 binding to PI(4,5)P2 and long chain heparins are mutually exclusive and that FGF2 has higher binding affinity for heparan sulfates and long-chain heparins. Using GUVs the authors confirm that PI(4,5)P2 and heparan sulfates are sufficient to drive some FGF2 translocation across the membrane). They also confirm that FGF2 must oligomerize for translocation. Using supported lipid bilayers and single molecule techniques FGF2, the authors find that FGF2 oligomers are heterogeneous containing a range of ~8-12 subunits. Molecular dynamic simulations of early steps of FGF2 binding to PI(4,5)P2 confirmed that residues identified in previous mutational analysis studies are important for PI(4,5)P2 binding and that there may be previously unknown residues involved in PI(4,5)P2 binding when several PI(4,5)P2 head groups bind FGF2 at the same time.

The novelty of this manuscript is the establishment of an in vitro system that has allowed the authors to directly test the minimal components FGF2 translocation. However, the results appear to confirm previous finding about FGF2 translocation that were discovered using other approaches. While solid and important work, it may not be of high enough impact for publication in eLife. The paper would be significantly strengthened if the authors were able to show that their work has led to novel, or previously unknown mechanism of FGF2 translocation.

1) It would be helpful to see the entire HSQCs for the binding titrations. Are the residues that bind PI(4,5)P2 and Heparan the only ones that shift?

2) How did the authors confirm that the GUVs contained PI(4,5)P2 on their surfaces and contain long-chain heparins in their lumen?

3) How does the concentration of FGF2 used in the in vitro experiments compare to the amount of FGF2 present in cells. If the amount of FGF2 added to the GUVs is lowered does translocation still occur?

4) The last part of the paper present molecular dynamics studies that suggests that there are previously uncharacterized residues on FGF2 that are important for translocation. Were these tested in the in vitro system and in cells to confirm their importance?

Reviewer #3:

The revision of this manuscript makes important changes and adds significant new data involving the roles of cholesterol, the lack of effect of heparin disaccharide in translocation, and informative new single molecule studies. As a result this already strong report is further improved and will be of interest to a range of readers. Publication in its present form would be appropriate.

There are however, three points that might be considered. Not that they need to be fully resolved before acceptance of the work, but the issues might be worth taking into account for interpretation.

First, the inability of heparin disaccharides to mimic the effect of heparin might be interpreted differently. The result is analogous to a PIP2-binding protein that does not bind IP3 well. If the competition between PIP2 and heparin were analogous to an enzyme inhibitor and a substrate, a charge-matched heparin construct should compete for PIP2 if the concentration is enough. The need for a heparin polymer suggests that polyelectrolyte effects in addition to or instead of lock-key binding is relevant here. That was the reason to ask if HA can mimic heparin if it's packed in the vesicle lumen at high concentration. Heparin has nominally twice the charge density of HA, but maybe high concentration of HA would work. Without more data about e.g. the hydration state of the lipid and heparin bound to the FGF, it seems hard to distinguish traditional "specific" binding from equally specific charge density and polyelectrolyte effects.

Second, the experiments here are all done with only monovalent cations in the buffers. But in vivo, heparin is in a medium with mM Mg2+ and perhaps more important mM Ca2+. Most likely addition of extracellular concentration of Ca2+ or even intracellular Mg2+ makes the experiments intractable, but at least some mention of divalents and what they are known from literature to do to PIP2 and heparin would be useful.

Third, in principle some of the issues can be addressed by MD simulation, but the system here is too briefly described and possibly too limited to get into the atomic detail needed to get an accurate picture of the charge interactions. For example, when PIP2 or heparin bind FGF, are waters lost or retained at the interface of the negative oxygen of the ligand and the positive nitrogen of the protein? The fact that the net charge of PIP2 is –4 suggests that there are limits in the MD. Recent work from Radhakrishnan's group show that a fuller treatment of PIP2 under realistic conditions reduces the net charge closer to –3, which alters the way the lipid orients and binds ligands.

As mentioned above, none of these issues are probably easily resolvable and they need not delay this study. But perhaps addressing the issues would strengthen the interpretation.

https://doi.org/10.7554/eLife.28985.037

Author response

[Editors’ note: the author responses to the first round of peer review follow.]

As you will see from the reviewer's comments, the biochemical reconstitution of FGF2 translocation is considered a biochemical tour de force. Discussion amongst the reviewers, however, resulted in a unanimous view that, at present, the manuscript largely confirms the molecular requirements for FGF2 translocation that were identified previously. The reviewers and I encourage you to continue this line of experimentation, and to this end, they make several suggestions that may help you to provide new information regarding the mechanism of FGF2 translocation. Specific questions that the reviewers felt were especially important to address and expand on regard the nature of the FGF2 translocation intermediate and PIP2 dynamics. Specifically, does FGF2 form a structural pore (versus simply rupturing the membrane) and, related, does PIP2 translocate with FGF2, or remain on the 'donor' leaflet? Though the reviewers concluded that suitable revisions to the manuscript would likely substantial effort, they expressed enthusiasm for considering a manuscript that does contain additional information.

We have submitted a revised version of the manuscript that contains both additional technical controls and two entirely new sections providing novel insights into the molecular mechanism of FGF2 membrane translocation. Beyond the reconstitution of FGF2 membrane translocation with purified components and the demonstration of mutually exclusive interactions of FGF2 with PI(4,5)P2 versus heparin, we have now added a series of single molecule studies that functionally correlate the oligomeric state of PI(4,5)P2 induced FGF2 oligomers with membrane insertion and pore formation. Furthermore, the manuscript now contains a series of atomistic molecular dynamics simulations that provide novel insights into the binding mechanism of FGF2 to PI(4,5)P2 as well as into the molecular mechanism of the initial steps of FGF2 oligomerization triggered by PI(4,5)P2. With these modifications, the revised manuscript addresses the three major aspects that were raised by the reviewers:

1) We have added additional controls to the core part of the original manuscript. This concerns experiments addressing questions of specificity when different forms of heparin are used as luminal traps for FGF2 in GUVs (new Figures 4 and 9) as well as a potential role of cholesterol in PI(4,5)P2 dependent FGF2 membrane translocation (new Figures 5 and 9).

2) We have added an entirely new section providing a functional single molecule

analysis of membrane-inserted FGF2 oligomers (new Figs. 10, 11 and 12). Following the suggestions of the reviewers, these studies provide the first insights into the structure-function relationship of membrane-inserted FGF2 oligomers. Our findings demonstrate PI(4,5)P2 dependent FGF2 oligomers to represent dynamic structures in the range of 10 to 12 subunits that form membrane-inserted translocation intermediates during unconventional secretion of FGF2. These findings are discussed in the context of a toroidal membrane pore and the assembly/disassembly hypothesis of FGF2 membrane translocation we have put forward previously [reviewed in (1)].

3) We added a new section based upon atomistic molecular dynamics simulations providing clues about the molecular mechanism of FGF2 binding to PI(4,5)P2 and the initial steps of PI(4,5)P2 induced oligomerization of FGF2 (new Figure 13 and 14, Figure 13—figure supplements 1-8 and Figure 14—figure supplements 1-3 and Video 1 and 2). Key aspects of this part are (i) the identification of additional binding sites in FGF2 resulting in multiple interactions of one FGF2 molecule with several PI(4,5)P2 molecules and (ii) the identification of a high affinity dimer interface in FGF2 highlighting a disulfide bridge connecting two cysteine 95 residues. These results are discussed in the context of the FGF2 membrane translocation model we are putting forward with the current manuscript.

Reviewer #1:

1.) Oligomeric state of membrane-inserted FGF2 forms and size limit of membrane

pores generated by PI(4,5)P2 induced FGF2 oligomers

In previous studies, we have conducted experiments analyzing a potential size cutoff of membrane pores generated by PI(4,5)P2 dependent FGF2 oligomerization and membrane insertion (2). While we found an Alexa488 dye (1 kDa) to physically traverse the membrane, the same dye coupled to either a dextran molecule (10 kDa) or cytochrome c (12 kDa) did not. Furthermore, we have previously shown that, upon membrane insertion of FGF2 oligomers, transbilayer diffusion of membrane lipids occurs providing direct evidence for the formation of a toroidal pore (2). These findings suggest that FGF2 oligomers insert into membranes through electrostatic interactions with PI(4,5)P2 thereby opening up a lipidic pore with a toroidal structure. In the revised manuscript, we have now included a series of single molecule measurements using both GUVs and supported lipid bilayers functionally correlating the oligomeric state of FGF2 with membrane pore formation (new Figures 10, 11 and 12). We find membrane-inserted FGF2 oligomers to represent dynamic structures with pore forming species characterized by about 10 to 12 FGF2 subunits. As discussed in the revised manuscript, these findings are consistent with a toroidal architecture and point to a dynamic structure that is continuously assembled at the inner leaflet and disassembled at the outer leaflet resulting in net transport of FGF2 molecules across the membrane.

2) The role of ATP1A1 in FGF2 membrane translocation

The fact that FGF2 membrane translocation can be observed in the absence of ATP1A1 does not imply that our reconstitution system is unrelated to the physiological mechanism. The most simple explanation would be that ATP1A1 is required in cells as a first physical contact to accumulate FGF2 at the inner leaflet of the plasma membrane at sites that provide proximity to Tec kinase and/or PI(4,5)P2. This requirement might be simply bypassed in the reconstitution system when the essential components are present in a purified form. However, provided ATP1A1 could be included in a proper way into our reconstitution system, it is possible that the efficiency of FGF2 membrane translocation would increase. Therefore, while ATP1A1 clearly is not essential for FGF2 membrane translocation to occur, we are indeed investigating ways of reconstituting ATP1A1 as part of our FGF2 membrane translocation assay. However, as noted by reviewer #1, reconstituting multi-subunit integral membrane proteins in a functional form is technically challenging and beyond the focus of the current manuscript.

Reviewer #2:

[…] 1) How does FGF2 traverse the hydrophobic part of the membrane?

The results from various studies point at a toroidal architecture of FGF2 induced membrane pores where the FGF2 oligomer is located in the center with each monomer bound to PI(4,5)P2 [reviewed in (1)]. The combined findings led to the assembly/disassembly model of FGF2 membrane translocation where FGF2 subunits are added to membrane-inserted FGF2 oligomers at the inner leaflet followed by removal of FGF2 subunits at the outer leaflet. The latter step is mediated by heparan sulfates that outcompete PI(4,5)P2 with regard to binding to FGF2. This process results in directional transport across the plasma membrane along with FGF2 trapping in the extracellular space where FGF2 remains bound to heparan sulfates on cell surfaces. As also explained in our response to reviewer #1, this view is supported by several lines of evidence. First, a toroidal architecture of FGF2 membrane pores is indicated by transbilayer diffusion of membrane lipids observed following membrane insertion of FGF2 oligomers (2). These results were corroborated by findings demonstrating that FGF2 oligomers can be extracted from membranes using carbonate, a direct indication of electrostatic rather than hydrophobic interactions with the membrane (3). Finally, the assembly/disassembly model as part of a toroidal pore is strongly supported by the current study demonstrating mutually exclusive interactions of FGF2 with PI(4,5)P2 and heparin, the latter mimicking cellular heparan sulfates outcompeting PI(4,5)P2 and thereby trapping FGF2 at the extracellular leaflet. In conclusion, all data from various studies point to a mechanism of FGF2 membrane translocation that is characterized by a lack of physical contacts of FGF2 with the hydrophobic core of the membrane. With regard to transbilayer diffusion of membrane lipids (2), it is indeed likely that this phenomenon includes PI(4,5)P2 itself as part of the assembly/disassembly mechanism of FGF2 membrane translocation. As suggested by reviewer #2, in the revised manuscript these aspects are now discussed in more detail along with the relevant literature on domain induced PIP2 flipping in cellular membranes.

2) Does FGF2 exclusively bind to the headgroup of PI(4,5)P2 ?

There is no evidence that FGF2 interacts with the hydrophobic parts of PI(4,5)P2. Using atomistic molecular dynamics simulations, we have further tested this possibility and, consistent with previous biochemical and structural data, found a high affinity binding site for the headgroup of PI(4,5)P2. (new Figure 13 and 14, Figure 13—figure supplements 1-8 and Figure 14—figure supplements 1-3 and Video 1 and 2). Intriguingly, these studies revealed additional binding sites in FGF2 for PI(4,5)P2, suggesting that FGF2 accumulates PI(4,5)P2 molecules at the sites of oligomerization and membrane pore formation. These findings are of extraordinary relevance for the formation of membrane pores with a toroidal architecture as we have shown previously that PI(4,5)P2 is required for FGF2 membrane insertion based upon its ability to stabilize the high membrane curvature that characterizes toroidal pores (2). The molecular dynamics studies included into the revised manuscript further provide clues on the molecular mechanism of the initial steps of PI(4,5)P2 induced FGF2 oligomerization which will be highly relevant for future studies analyzing the fine structure of FGF2 oligomers as part of toroidal pores.

3) How specific is the effect of heparin? Could it be replaced by for example hyarulonic acid?

FGF family members are well characterized with regard to their ligand binding specificities. FGF2 has been shown to interact exclusively with sulphated glycosaminoglycans as contained in heparan sulfate proteoglycans. As opposed to heparan sulfates and other sulfated proteoglycans, cellular forms of hyaluronic acid do not contain sulfate groups. In any case, reviewer #2 inspired us to systematically compare long-chain heparins with a defined heparin disaccharide (as used in Figures 1 and 2) in FGF2 membrane translocation assays. As shown in the new Figure 4 and quantified in the new Figure 9, replacing long-chain heparins by a heparin disaccharide does not have any impact on membrane pore formation, however, FGF2 membrane translocation is almost fully blocked. These findings are in line with the NMR data and biochemical experiments shown in Figures 1 and 2 demonstrating long-chain heparins to be required to outcompete PI(4,5)P2 with regard to binding to FGF2. Our combined results therefore corroborate the proposed model for the molecular mechanism of directional membrane translocation of FGF2 based on sequential interactions with PI(4,5)P2 and heparan sufates on opposing sides of the membrane.

4) Does heparin change the osmotic pressure and, in turn, membrane tension?

As explained in “Materials and methods”, the generation procedures of all types and batches of GUVs were controlled by osmolality measurements and demonstrated not to be affected by adding either long-chain or disaccharide heparins. Therefore, artifacts related to differences in osmotic pressure and membrane tension can be excluded.

5) Does cholesterol play a role in FGF2 membrane translocation?

As suggested by reviewer #2, we have conducted reconstitution experiments with heparin-containing GUVs that either contained or lacked cholesterol. As shown in the new Figure 5 and quantified in the new Figure 9, omission of cholesterol causes a substantial decrease of both membrane pore formation and translocation of FGF2. The experiments shown in Figure 5 indicate that, in the absence of cholesterol, binding of FGF2 to PI(4,5)P2 is strongly impaired. This is consistent with previous studies from our laboratory demonstrating strong binding of FGF2 to PI(4,5)P2 in liposomes with a plasma membrane like lipid composition whereas binding was found poor when PI(4,5)P2 was reconstituted in pure PC liposomes (4,5). Our combined findings suggest that components of the FGF2 secretion apparatus such as PI(4,5)P2 are organized in cholesterol dependent microdomains that are required for efficient recruitment, oligomerization and membrane translocation of FGF2. The new results are shown in Figures 5 and 9 in the revised manuscript and are discussed accordingly.

Reviewer #3:

1) What could be the role of ATP1A1 as a novel aspect of FGF2 membrane translocation?

In previous studies, we have demonstrated a direct interaction of FGF2 with the cytoplasmic domain of ATP1A1 (6). As explained in our response to reviewer #1, the function of ATP1A1 in FGF2 secretion from cells might be a first physical contact of FGF2 with the inner leaflet of plasma membranes. One possible explanation for a requirement of ATP1A1 in FGF2 secretion from cells might be a need to accumulate FGF2 in plasma membrane microdomains mediating proximity to additional components such as Tec kinase and PI(4,5)P2. In our reconstitution system, this function might be bypassed when the essential components are present in a purified form. While we do run projects analyzing the function of ATP1A1 both in a cellular context and in reconstitution systems, these studies are clearly beyond the focus of the current manuscript.

2) Do cysteine residues transplanted to FGF family members with signal peptides promote membrane translocation similar to FGF2?

We have done such experiments in the context of cells (7). As a model FGF protein carrying a signal peptide, FGF4 has been used as it is a close structural relative of FGF2. Following removal of the signal peptide of FGF4, secretion became undetectable, however, after transplanting the critical cis-elements including Cys77 and Cys95, the corresponding FGF4/2 hybrid protein was indeed found rerouted into the unconventional secretory pathway of FGF2 (7). Accordingly, a double cysteine mutant of FGF2 (C77/95A) was used in the current study and demonstrated to be strongly impaired in membrane translocation (Figures 7 and 9).

References:

1. La Venuta, G., Zeitler, M., Steringer, J. P., Müller, H. M., and Nickel, W. (2015) The Startling Properties of Fibroblast Growth Factor 2: How to Exit Mammalian Cells without a Signal Peptide at Hand. J Biol Chem 290, 27015-27020.

2. Steringer, J. P., Bleicken, S., Andreas, H., Zacherl, S., Laussmann, M., Temmerman, K., Contreras, F. X., Bharat, T. A., Lechner, J., Müller, H. M., Briggs, J. A., Garcia-Saez, A. J., and Nickel, W. (2012) Phosphatidylinositol 4,5-Bisphosphate (PI(4,5)P2)-dependent Oligomerization of Fibroblast Growth Factor 2 (FGF2) Triggers the Formation of a Lipidic Membrane Pore Implicated in Unconventional Secretion. J Biol Chem 287, 27659-27669.

3. Zeitler, M., Steringer, J. P., Müller, H. M., Mayer, M. P., and Nickel, W. (2015) HIV-Tat Protein Forms Phosphoinositide-dependent Membrane Pores Implicated in Unconventional Protein Secretion. J Biol Chem 290, 21976-21984.

4. Temmerman, K., Ebert, A. D., Müller, H. M., Sinning, I., Tews, I., and Nickel, W. (2008) A direct role for phosphatidylinositol-4,5-bisphosphate in unconventional secretion of fibroblast growth factor 2. Traffic 9, 1204-1217.

5. Temmerman, K., and Nickel, W. (2009) A novel flow cytometric assay to quantify interactions between proteins and membrane lipids. J Lipid Res.

6. Zacherl, S., La Venuta, G., Müller, H. M., Wegehingel, S., Dimou, E., Sehr, P., Lewis, J. D., Erfle, H., Pepperkok, R., and Nickel, W. (2015) A direct role for ATP1A1 in unconventional secretion of fibroblast growth factor 2. J Biol Chem 290, 3654-3665.

7. Müller, H. M., Steringer, J. P., Wegehingel, S., Bleicken, S., Munster, M., Dimou, E., Unger, S., Weidmann, G., Andreas, H., Garcia-Saez, A. J., Wild, K., Sinning, I., and Nickel, W. (2015) Formation of Disulfide Bridges Drives Oligomerization, Membrane Pore Formation and Translocation of Fibroblast Growth Factor 2 to Cell Surfaces. J Biol Chem 290, 8925-8937.

[Editors' note: the author responses to the re-review follow.]

[…] Reviewer #1:

The major issue raised during review of the original manuscript regarded the confirmatory nature of the manuscript, which confirmed specific requirements that were identified previously using cell based assays. In the revised manuscript, the authors now include several new findings that shed some light onto the mechanism of FGF2 translocation across a membrane. I find the most compelling of these to be the following:

1) Toroidal translocation intermediate

A lipidic membrane pore with a toroidal architecture containing the FGF2 oligomer in its center is supported by several independent lines of evidence from previous publications and data from our current work. However, the precise molecular mechanism by which membrane-inserted forms of FGF2 oligomers act as translocation intermediates is not fully understood. We have put forward the hypothesis that membrane-inserted FGF2 oligomers are dynamic structures to which FGF2 subunits are constantly added at the inner leaflet concomitant with the removal of subunits at the outer leaflet. In this way, at steady state, membrane-inserted forms of FGF2 have a dynamic oligomeric state in the range of 8 to 12 FGF2 subunits that are continuously remodeled by adding and removing subunits on opposing sides of the membrane. It will be a challenge for future studies to find out whether these assembly/disassembly units are monomers or disulfide-bridged dimers of FGF2, the latter being functional in FGF2 signaling as part of ternary complexes with heparan sulfates and FGF high affinity receptors. Furthermore, a big goal in our lab will be to elucidate the three-dimensional architecture of membrane inserted FGF2 oligomers to fully understand how these complexes are functioning as translocation intermediates in unconventional secretion of FGF2.

As requested by reviewer 1, we have modified the discussion to explain these aspects in more detail.

2) Modifications to Fig. 10 and the corresponding legend

We have streamlined Figure 10 with a simplified labelling of its subpanels that makes this figure more self-explanatory. Furthermore, the lipid tracer is now shown for both variant forms of FGF2 which helps to better distinguish the two experimental conditions being compared. The legend has been adapted accordingly.

Reviewer #2:

1) HSQC data from NMR experiments

The full set of HSQC data from the experiments shown in Figure 1 has been added to the source file with all raw data available to the reader.

2) Measurements of PI(4,5)P2 and heparin incorporation into GUVs

Luminal incorporation of heparin into GUVs was monitored by confocal microscopy in control experiments using a fluorescent derivative. Likewise, the presence of PI(4,5)P2 in the bilayers of GUVs was analyzed using a recombinant fusion protein consisting of GFP and the PH domain of PLCδ1, a canonical PI(4,5)P2 marker. The ‘Materials and methods’ section has been modified accordingly along with the corresponding references.

3) How does the concentration of FGF2 in cells compare with FGF2 concentrations usedin reconstitution experiments?

Though the precise intracellular concentrations of FGF2 in μM has not been reported in the literature, tumor cells massively overexpress FGF2 under pathophysiological conditions which drives an autocrine secretion/signaling loop of FGF2. In our reconstitution experiments, FGF2 is used at low concentrations between 50 and 200 nM (to allow for single molecule measurements) that are likely to be way below the intracellular concentration of FGF2 in tumor cells.

4) Analysis of additional residues in FGF2 with affinity to PI(4,5)P2 as identified by MD simulations

The residues we have identified previously in biochemical experiments (K127, R128 and K133) form a high affinity binding pocket for PI(4,5)P2 [Temmerman et al., (2008) Traffic 9:1204-1217]. Consistent with both biochemical experiments and MD simulations, these residues are essential for interactions of FGF2 with PI(4,5)P2 containing membranes and are required for FGF2 secretion from cells. The identification of additional residues in FGF2 mediating PI(4,5)P2 clustering as identified in MD simulations is a novel finding of the current work with important implications for future studies. However, their detailed analysis in reconstitution experiments and cell-based assays is beyond the scope of the current manuscript.

Reviewer #3:

1) Electrostatic aspects of mutually exclusive interactions of FGF2 with various forms of heparin and PI(4,5)P2

The electrostatic determinants of FGF2 binding to PI(4,5)P2, heparin isaccharides and long-chain heparins are an important aspect with regard to the mechanism of FGF2 membrane translocation we are proposing. However, it is clear from the literature that all of these interactions are not merely based upon the number of electrostatic pairs. Rather, both the orientation of sugar subunits in the heparin backbone and the positioning of sulfate moieties along the ring structure are critical determinants for FGF2 binding. There are in fact examples for different kinds of heparin molecules that are identical in net charge but differ in the positioning of the sulfate groups resulting in strong binding versus a complete lack of FGF2 binding [Li et al., (2014) ACS Chem Biol 9:1712-1717) and papers therein]. Likewise, we have previously demonstrated that FGF2 binding to PI(4,5)P2 is much stronger compared to PI(3,4,5)P3 [Temmerman et al., (2008) Traffic 9:1204-1217], pointing at a defined binding pocket for PI(4,5)P2 in FGF2 as demonstrated by our NMR experiments. Therefore, hyaluronic acid that lacks sulfate groups even at high concentrations is very unlikely to outcompete PI(4,5)P2 with regard to binding to FGF2. As outlined in our manuscript, there are two defined and overlapping binding epitopes in FGF2 for heparin and PI(4,5)P2 resulting in mutually exclusive interactions, a key aspect of the translocation model we are proposing.

2) Potential effects of divalent cations on FGF2-heparin interactions

As indicated by reviewer 3, the use of divalent calcium and magnesium ions in biochemical reconstitution experiments employing GUVs is technically challenging. However, in cell-based FGF2 secretion experiments, we found that both calcium (1.8mM) and magnesium ions (0.8 mM) in the extracellular medium does not affect the ability of heparan sulfates to trigger FGF2 membrane translocation and to retain FGF2 on cell surfaces (Zehe et al., 2006, Proc. Natl. Acad. Sci. U.S.A. 103:15479-15484). These observations in intact cells demonstrate that the mechanism we are proposing is functional even in the presence of millimolar concentrations of divalent cations.

3) MD simulations of the interaction between FGF2 and PI(4,5)P2.

Reviewer 3 commented on the net charge of PI(4,5)P2 and referred to recent work by the Radhakrishnan group [e.g., J Phys Chem B 117, 8322 (2013)]. Their quantum mechanical calculations indicated a net charge of PI(4,5)P2 to be –4 around pH 7. This is fully consistent with our work. However, deviations from the net charge of –4 were suggested in the presence of divalent ions. While the effects of divalent salt and its influence on the net charge of PI(4,5)P2 could be explored in MD simulations, it would require a computational resource that is larger by a factor of 5–10 compared to the resources used for the present work. In addition, the nature of the results is not expected to change as both –3 and –4 are exceptionally large charges for a lipid head group and will be predominate over any other lipid in terms of electrostatic interactions.

Another aspect raised by reviewer 3 concerned the role of water molecules in the course of FGF2 binding to PI(4,5)P2.To address this question, we analyzed our simulation data carefully by computing radial distribution functions, solvent accessible surface areas, number of water molecules around the PI(4,5)P2 binding pocket in FGF2 and number of hydrogen bonds between FGF2 binding pocket residues and water molecules. The best overall picture is given by the number of water molecules around the binding pocket region (residues K127, R128, K133). As shown in Author response image 1, the number of water molecules is about 95 when FGF2 is not in contact with PI(4,5)P2:

Author response image 1
Number of water molecules in the vicinity of the PI(4,5)P2 binding pocket in FGF2 in the course of time.

The PI(4,5)P2 head group binds to the binding pocket at about 500 ns. Here, a water molecule is considered to be in the vicinity of the binding pocket if its distance from the main residues (K127, R128, K133) is 0.6 nm or less.

https://doi.org/10.7554/eLife.28985.036

When PI(4,5)P2binds to FGF2 by occupying the binding pocket at about 500 ns, the number of water molecules around the binding pocket region decreases by ~25% to about 70. This is quite expected, since the PI(4,5)P2 head group blocks some of the water molecules away from the binding pocket region. Nonetheless, the binding pocket remains very hydrated, i.e. the hydration level does not change in any essential manner.

https://doi.org/10.7554/eLife.28985.038

Article and author information

Author details

  1. Julia P Steringer

    Heidelberg University Biochemistry Center, Heidelberg, Germany
    Contribution
    JPS, Conceptualization, Data curation, Investigation
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9418-2762
  2. Sascha Lange

    Institut für Chemie und Biochemie, Freie Universität Berlin, Berlin, Germany
    Present address
    Department of Molecular Biophysics, Leibniz-Forschungsinstitut fürMolekulare Pharmakologie, Berlin, Germany
    Contribution
    SL, Conceptualization, Data curation, Investigation, Methodology
    Contributed equally with
    Sabína Čujová, Radek Šachl, Chetan Poojari, Fabio Lolicato, Oliver Beutel, Ünal Coskun, Alf Honigmann, Ilpo Vattulainen, Martin Hof and Christian Freund
    Competing interests
    The authors declare that no competing interests exist.
  3. Sabína Čujová

    J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Prague, Czech Republic
    Contribution
    SČ, Data curation, Investigation, Methodology
    Contributed equally with
    Sascha Lange, Radek Šachl, Chetan Poojari, Fabio Lolicato, Oliver Beutel, Ünal Coskun, Alf Honigmann, Ilpo Vattulainen, Martin Hof and Christian Freund
    Competing interests
    The authors declare that no competing interests exist.
  4. Radek Šachl

    J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Prague, Czech Republic
    Contribution
    RŠ, Conceptualization, Data curation, Formal analysis, Investigation, Methodology
    Contributed equally with
    Sascha Lange, Sabína Čujová, Chetan Poojari, Fabio Lolicato, Oliver Beutel, Ünal Coskun, Alf Honigmann, Ilpo Vattulainen, Martin Hof and Christian Freund
    Competing interests
    The authors declare that no competing interests exist.
  5. Chetan Poojari

    1. Department of Physics, University of Helsinki, Helsinki, Finland
    2. Department of Physics, Tampere University of Technology, Tampere, Finland
    Contribution
    CP, Data curation, Investigation, Methodology
    Contributed equally with
    Sascha Lange, Sabína Čujová, Radek Šachl, Fabio Lolicato, Oliver Beutel, Ünal Coskun, Alf Honigmann, Ilpo Vattulainen, Martin Hof and Christian Freund
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-6575-221X
  6. Fabio Lolicato

    1. Department of Physics, University of Helsinki, Helsinki, Finland
    2. Department of Physics, Tampere University of Technology, Tampere, Finland
    Contribution
    FL, Data curation, Formal analysis, Investigation, Methodology
    Contributed equally with
    Sascha Lange, Sabína Čujová, Radek Šachl, Chetan Poojari, Oliver Beutel, Ünal Coskun, Alf Honigmann, Ilpo Vattulainen, Martin Hof and Christian Freund
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-7537-0549
  7. Oliver Beutel

    Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany
    Contribution
    OB, Data curation, Methodology
    Contributed equally with
    Sascha Lange, Sabína Čujová, Radek Šachl, Chetan Poojari, Fabio Lolicato, Ünal Coskun, Alf Honigmann, Ilpo Vattulainen, Martin Hof and Christian Freund
    Competing interests
    The authors declare that no competing interests exist.
  8. Hans-Michael Müller

    Heidelberg University Biochemistry Center, Heidelberg, Germany
    Contribution
    H-MM, Data curation, Methodology, Writing—review and editing
    Competing interests
    The authors declare that no competing interests exist.
  9. Sebastian Unger

    Heidelberg University Biochemistry Center, Heidelberg, Germany
    Contribution
    SU, Data curation, Methodology
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2525-8502
  10. Ünal Coskun

    Paul Langerhans Institute Dresden, Dresden, Germany
    Contribution
    ÜC, Conceptualization, Resources, Writing—review and editing
    Contributed equally with
    Sascha Lange, Sabína Čujová, Radek Šachl, Chetan Poojari, Fabio Lolicato, Oliver Beutel, Alf Honigmann, Ilpo Vattulainen, Martin Hof and Christian Freund
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4375-3144
  11. Alf Honigmann

    Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany
    Contribution
    AH, Conceptualization, Data curation, Methodology, Writing—review and editing
    Contributed equally with
    Sascha Lange, Sabína Čujová, Radek Šachl, Chetan Poojari, Fabio Lolicato, Oliver Beutel, Ünal Coskun, Ilpo Vattulainen, Martin Hof and Christian Freund
    Competing interests
    The authors declare that no competing interests exist.
  12. Ilpo Vattulainen

    1. Department of Physics, University of Helsinki, Helsinki, Finland
    2. Department of Physics, Tampere University of Technology, Tampere, Finland
    3. MEMPHYS – Center for Biomembrane Physics, University of Southern Denmark, Denmark, United Kingdom
    Contribution
    IV, Conceptualization, Resources, Data curation, Methodology, Writing—review and editing
    Contributed equally with
    Sascha Lange, Sabína Čujová, Radek Šachl, Chetan Poojari, Fabio Lolicato, Oliver Beutel, Ünal Coskun, Alf Honigmann, Martin Hof and Christian Freund
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-7408-3214
  13. Martin Hof

    J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Prague, Czech Republic
    Contribution
    MH, Conceptualization, Resources, Methodology, Writing—review and editing
    Contributed equally with
    Sascha Lange, Sabína Čujová, Radek Šachl, Chetan Poojari, Fabio Lolicato, Oliver Beutel, Ünal Coskun, Alf Honigmann, Ilpo Vattulainen and Christian Freund
    Competing interests
    The authors declare that no competing interests exist.
  14. Christian Freund

    Institut für Chemie und Biochemie, Freie Universität Berlin, Berlin, Germany
    Contribution
    CF, Conceptualization, Data curation, Writing—review and editing
    Contributed equally with
    Sascha Lange, Sabína Čujová, Radek Šachl, Chetan Poojari, Fabio Lolicato, Oliver Beutel, Ünal Coskun, Alf Honigmann, Ilpo Vattulainen and Martin Hof
    Competing interests
    The authors declare that no competing interests exist.
  15. Walter Nickel

    Heidelberg University Biochemistry Center, Heidelberg, Germany
    Contribution
    WN, Conceptualization, Formal analysis, Supervision, Funding acquisition, Writing—original draft, Project administration, Writing—review and editing
    For correspondence
    walter.nickel@bzh.uni-heidelberg.de
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6496-8286

Funding

Czech Science Foundation (14-03141)

  • Martin Hof
  • Sabína Čujová
  • Radek Šachl

European Research Council (Advanced Grant 290974 CROWDED-PRO-LIPIDS)

  • Ilpo Vattulainen
  • Chetan Poojari
  • Fabio Lolicato

Academy of Finland (Center of Excellence project 272130)

  • Chetan Poojari
  • Fabio Lolicato
  • Ilpo Vattulainen

Sigrid Juselius Foundation

  • Ilpo Vattulainen
  • Chetan Poojari
  • Fabio Lolicato

Academy of Finland (Center of Excellence project 307415)

  • Chetan Poojari
  • Fabio Lolicato
  • Ilpo Vattulainen

Deutsche Forschungsgemeinschaft (SFB 854)

  • Christian Freund

Deutsche Forschungsgemeinschaft (SFB 958)

  • Christian Freund

Deutsche Forschungsgemeinschaft (SFB/TRR 186)

  • Walter Nickel
  • Christian Freund

Deutsche Forschungsgemeinschaft (DFG Ni 423/6-1)

  • Walter Nickel

Deutsche Forschungsgemeinschaft (SFB/TRR 83)

  • Walter Nickel

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

The Abberior StarRed-HALO-Tag ligand was generously provided by Vladimir Belov (MPI-BPC, Göttingen, Germany). WN was supported by the German Research Foundation (DFG SFB/TRR 186, DFG SFB/TRR 83 and DFG Ni 423/6–1) and the DFG Cluster of Excellence at the University of Heidelberg. MH, SC and RS were supported by grant 14–03141 provided by the Czech Science Foundation and Praemium Academie award given by the Czech Academy of Sciences. CP, FL and IV were supported by the European Research Council (Advanced Grant 290974 CROWDED-PRO-LIPIDS), the Academy of Finland (Center of Excellence projects 272130 and 307415) and the Sigrid Juselius Foundation. CF was supported by the German Research Foundation as part of DFG SFB 854, DFG SFB 958 and DFG SFB/TRR 186. The authors wish to acknowledge CSC IT Center for Science (Espoo, Finland) and the Center for Information Services and High Performance Computing (ZIH (Technische Universität Dresden, Germany) for computational resources.

Reviewing Editor

  1. Christopher G Burd, Yale School of Medicine, United States

Publication history

  1. Received: May 25, 2017
  2. Accepted: July 14, 2017
  3. Accepted Manuscript published: July 19, 2017 (version 1)
  4. Accepted Manuscript updated: July 27, 2017 (version 2)
  5. Version of Record published: September 15, 2017 (version 3)
  6. Version of Record updated: October 5, 2017 (version 4)

Copyright

© 2017, Steringer et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 3,652
    Page views
  • 618
    Downloads
  • 49
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Julia P Steringer
  2. Sascha Lange
  3. Sabína Čujová
  4. Radek Šachl
  5. Chetan Poojari
  6. Fabio Lolicato
  7. Oliver Beutel
  8. Hans-Michael Müller
  9. Sebastian Unger
  10. Ünal Coskun
  11. Alf Honigmann
  12. Ilpo Vattulainen
  13. Martin Hof
  14. Christian Freund
  15. Walter Nickel
(2017)
Key steps in unconventional secretion of fibroblast growth factor 2 reconstituted with purified components
eLife 6:e28985.
https://doi.org/10.7554/eLife.28985

Further reading

    1. Biochemistry and Chemical Biology
    2. Structural Biology and Molecular Biophysics
    Dario Segura-Peña, Oda Hovet ... Nikolina Sekulic
    Research Article Updated

    Aurora B, together with IN-box, the C-terminal part of INCENP, forms an enzymatic complex that ensures faithful cell division. The [Aurora B/IN-box] complex is activated by autophosphorylation in the Aurora B activation loop and in IN-box, but it is not clear how these phosphorylations activate the enzyme. We used a combination of experimental and computational studies to investigate the effects of phosphorylation on the molecular dynamics and structure of [Aurora B/IN-box]. In addition, we generated partially phosphorylated intermediates to analyze the contribution of each phosphorylation independently. We found that the dynamics of Aurora and IN-box are interconnected, and IN-box plays both positive and negative regulatory roles depending on the phosphorylation status of the enzyme complex. Phosphorylation in the activation loop of Aurora B occurs intramolecularly and prepares the enzyme complex for activation, but two phosphorylated sites are synergistically responsible for full enzyme activity.

    1. Biochemistry and Chemical Biology
    2. Structural Biology and Molecular Biophysics
    Maura Greiser, Mariusz Karbowski ... Liron Boyman
    Research Article

    Mitochondrial ATP production in cardiac ventricular myocytes must be continually adjusted to rapidly replenish the ATP consumed by the working heart. Two systems are known to be critical in this regulation: mitochondrial matrix Ca2+ ([Ca2+]m) and blood flow that is tuned by local ventricular myocyte metabolic signaling. However, these two regulatory systems do not fully account for the physiological range of ATP consumption observed. We report here on the identity, location, and signaling cascade of a third regulatory system -- CO2/bicarbonate. CO2 is generated in the mitochondrial matrix as a metabolic waste product of the oxidation of nutrients that powers ATP production. It is a lipid soluble gas that rapidly permeates the inner mitochondrial membrane (IMM) and produces bicarbonate (HCO3-) in a reaction accelerated by carbonic anhydrase (CA). The bicarbonate level is tracked physiologically by a bicarbonate-activated adenylyl cyclase, soluble adenylyl cyclase (sAC). Using structural Airyscan super-resolution imaging and functional measurements we find that sAC is primarily inside the mitochondria of ventricular myocytes where it generates cAMP when activated by HCO3-. Our data strongly suggest that ATP production in these mitochondria is regulated by this cAMP signaling cascade operating within the inter-membrane space (IMS) by activating local EPAC1 (Exchange Protein directly Activated by cAMP) which turns on Rap1 (Ras-related protein 1). Thus, mitochondrial ATP production is shown to be increased by bicarbonate-triggered sAC signaling through Rap1. Additional evidence is presented indicating that the cAMP signaling itself does not occur directly in the matrix. We also show that this third signaling process involving bicarbonate and sAC activates the cardiac mitochondrial ATP production machinery by working independently of, yet in conjunction with, [Ca2+]m-dependent ATP production to meet the energy needs of cellular activity in both health and disease. We propose that the bicarbonate and calcium signaling arms function in a resonant or complementary manner to match mitochondrial ATP production to the full range of energy consumption in cardiac ventricular myocytes in health and disease.