Dilation of fusion pores by crowding of SNARE proteins
Abstract
Hormones and neurotransmitters are released through fluctuating exocytotic fusion pores that can flicker open and shut multiple times. Cargo release and vesicle recycling depend on the fate of the pore, which may reseal or dilate irreversibly. Pore nucleation requires zippering between vesicle-associated v-SNAREs and target membrane t-SNAREs, but the mechanisms governing the subsequent pore dilation are not understood. Here, we probed the dilation of single fusion pores using v-SNARE-reconstituted ~23-nm-diameter discoidal nanolipoprotein particles (vNLPs) as fusion partners with cells ectopically expressing cognate, 'flipped' t-SNAREs. Pore nucleation required a minimum of two v-SNAREs per NLP face, and further increases in v-SNARE copy numbers did not affect nucleation rate. By contrast, the probability of pore dilation increased with increasing v-SNARE copies and was far from saturating at 15 v-SNARE copies per face, the NLP capacity. Our experimental and computational results suggest that SNARE availability may be pivotal in determining whether neurotransmitters or hormones are released through a transient ('kiss and run') or an irreversibly dilating pore (full fusion).
https://doi.org/10.7554/eLife.22964.001Introduction
Fusion pores are nanoscale connections between membrane-enclosed compartments that are key intermediates during membrane fusion reactions such as the exocytotic release of neurotransmitters and hormones (Lindau et al., 2003). Following nucleation by specialized proteins (Chernomordik and Kozlov, 2008), fusion pores flicker repeatedly and then dilate or reseal during the release of hormones (Lindau et al., 2003) or neurotransmitters (Alabi and Tsien, 2013; He and Wu, 2007; Staal et al., 2004), or during fusion mediated by viral proteins (Cohen and Melikyan, 2004). The mechanisms that govern these behaviors are poorly understood, despite the availability of sensitive electrical and electrochemical methods to detect single fusion pores during protein-free fusion (Chanturiya et al., 1997; Mellander et al., 2014), viral-protein-induced fusion (Cohen and Melikyan, 2004), and exocytosis (Lindau, 2012). Even the very nature of the fusion pore intermediate (whether lipid- or protein-lined) is debated (Bao et al., 2016).
During exocytotic neurotransmitter or hormone release, a fusion pore opens as vesicle-associated soluble N-ethylmaleimide-sensitive factor attachment protein receptors (v-SNAREs) pair with cognate t-SNAREs on the target plasma membrane (Südhof and Rothman, 2009). This is a tightly regulated process that requires the coordinated actions of several proteins, including Munc18, Munc13, and others (Rizo and Xu, 2015). Complex formation between the v- and t-SNAREs is likely to start at the membrane distal N-termini and may proceed in stages toward the membrane-proximal regions (Gao et al., 2012). Assembly of the SNARE domains results in a four-helix bundle (SNAREpin) that brings the two bilayers into close proximity, but assembly is thought to be halted at some stage to poise vesicles for fast release. Calcium influx in response to depolarization is thought to lead to further SNARE assembly that promotes pore nucleation. This last step — coupling calcium entry to fusion — also requires Synaptotagmin and Complexin, which may actively contribute to pore opening. The initial fusion pore is a metastable structure that may reseal without ever dilating beyond ~1–2 nm in size. This results in transient ‘kiss and run’ exocytosis, a well-established mode of fusion for hormone-secreting cells (Alabi and Tsien, 2013; Hanna et al., 2009; Fulop et al., 2005). By contrast, whether transient fusion is a relevant mode of release for synaptic vesicle fusion has been a subject of debate (Alabi and Tsien, 2013; He and Wu, 2007; Staal et al., 2004; Pawlu et al., 2004; Chapochnikov et al., 2014), mainly because of technical challenges in probing fusion pores directly during synaptic release. For both neuronal and endocrine release, little is known about the molecular mechanisms that govern pore dilation and that set the balance between transient and full fusion (Alabi and Tsien, 2013), in large part due to a lack of biochemically defined assays that are sensitive to single-pores.
Here, using a recently developed nanodisc-cell fusion system (Wu et al., 2016), we found that the presence of just a few SNARE complexes can nucleate a pore, but that reliable pore dilation necessitates many more.
Results
Fusion between v-SNARE reconstituted nanolipoprotein particles and flipped t-SNARE cells
We used 21–27 nm diameter nanolipoprotein particles (NLPs) (Bello et al., 2016) to determine whether SNAREs alone can catalyze pore dilation. By contrast, most previous studies employed much smaller, 6–18-nm diameter nanodiscs (NDs) (Bao et al., 2016; Wu et al., 2016; Shi et al., 2012), whose dimensions restricted pore diameters to nm (Wu et al., 2016) and SNARE copy numbers to , precluding studies of pore dilation (Bello et al., 2016). We incorporated v-SNAREs into NLPs that were stabilized by a recombinant apolipoprotein E variant consisting of the N-terminal 22-kDa fragment (ApoE422k). We varied the lipid-to-ApoE422k ratio to control the NLP size and the VAMP2-to-ApoE422k ratio to tune the number of v-SNARE copies per NLP (Bello et al., 2016) (Figure 1). We confirmed that vNLPs fused with liposomes that were reconstituted with t-SNAREs in a SNARE-dependent manner using a previously described bulk assay that monitors calcium release through pores connecting v-SNARE nanodiscs with t-SNARE liposomes (Bello et al., 2016; Shi et al., 2012) (Figure 1—figure supplement 1a). Although NLP pores could in principle grow to nm in diameter (Figure 4b), much larger than the nm allowed by the membrane scaffold protein (MSP) based small ND geometry (Wu et al., 2016; Shi et al., 2012), bulk calcium release rates were comparable between vNLP and vMSP NDs loaded with similar v-SNARE copy numbers (Figure 1—figure supplement 1b), confirming that the bulk assay is largely insensitive to pore properties under these conditions (Bello et al., 2016).
We then confirmed lipid mixing between the membranes of vNLPs and flipped t-SNARE cells (tCells) using a previously described protocol (Wu et al., 2016) (Figure 2). NLPs co-labeled with one mole % each of DiI (donor) and DiD (acceptor) were incubated with tCells for 30 min at 4°C, a temperature that allows docking but not fusion. Cells were then rinsed to remove free NLPs and mounted onto a confocal microscope stage held at 37°C to initiate fusion and imaging of DiI and (directly excited) DiD fluorescence. At the concentrations used, the DiI fluorescence is quenched by DiD when the dyes are initially in the NLP membrane. Upon fusion, the dyes become diluted in the plasma membrane and the DiI fluorescence increases. Directly excitated DiD fluorescence provides a measure of the amount of docked NLPs. The ratio of the DiI to DiD intensity normalizes the lipid-mixing signal to the amount of docked NLPs. Normalized lipid-mixing signals increased when tCells were incubated with vNLPs carrying eight v-SNAREs total (vNLP8), but not with empty NLPs (eNLP) or with NLPs loaded with eight copies of a v-SNARE construct, VAMP2-4X, carrying mutations in the C-terminal hydrophobic layers (L70D, A74R, A81D, and L84D) (Figure 2b). These mutations prevent zippering of the C-terminal half of the SNARE domains (Wu et al., 2016; Krishnakumar et al., 2011, 2013), a perturbation that does not affect docking (Figure 2c) but prevents fusion (Figure 2b,d) (Wu et al., 2016; Krishnakumar et al., 2013). A variation of the assay that avoided the 4°C incubation but prevented live imaging in the presence of the labeled NLPs confirmed these results (Figure 2d). We estimate that 4–5% of the docked vNLPs undergo fusion with the flipped t-SNARE cells over the course of ~20 min (Figure 2—figure supplement 1). In comparison, fusion between v-SNARE NLPs and t-SNARE liposomes yields a similar extent of lipid mixing over the same period (Bello et al., 2016).
Lipid mixing could result from the merging of only the proximal lipid bilayer leaflets of the vNLPs and the tCells. To test whether full fusion occurred, we loaded the cells with Fluo-4, a fluorescent calcium probe, and monitored calcium signals. If full fusion occurred, then calcium influx through the fusion pores connecting vNLP and tCell membranes should increase cytosolic Fluo-4 signals (Wu et al., 2016). This was indeed the case for vNLP8 samples, but not for empty NLPs or for NLPs loaded with VAMP2-4X (Figure 3). Using this calcium-influx assay, we also assessed whether pores eventually resealed by washing away the free vNLPs after 5 min of incubation. Cellular calcium levels returned to baseline within a few minutes, suggesting that pores eventually resealed (Figure 3—figure supplement 1).
Dynamics of single fusion pores
Next, we probed single pores connecting vNLPs to tCell membranes (Wu et al., 2016) (Figure 4). We voltage-clamped a tCell, in the cell-attached configuration, that was ectopically expressing ‘flipped’ neuronal/exocytotic t-SNAREs syntaxin1 and SNAP25 (Hu et al., 2003) (Figure 4). NLPs that were reconstituted with eight copies of the complementary neuronal v-SNARE VAMP2/synaptobrevin (vNLPs) were included in the pipette solution (100 nM vNLPs, 120 µM lipid); they diffused to the pipette tip and fused with the patch. Because an NLP is not a closed structure like a vesicle, its fusion with the voltage-clamped membrane patch establishes a direct conduction pathway between the cytosol and the pipette solution, leading to direct-currents whose magnitude reflects pore geometry (Wu et al., 2016).
Currents appeared in bursts with an average frequency of ~0.2 bursts per min, or ~2.5 per patch (Figure 4c). The very low burst frequency, together with small unitary conductances (see below), strongly suggest that each burst represents currents passing through a single pore (Wu et al., 2016). Currents fluctuated and returned to baseline multiple times, as if the fusion pore fluctuated in size and opened and closed repeatedly, i.e. flickered (Figure 4d and Materials and methods). Nucleation was blocked when the cytoplasmic domain of VAMP2 (CDV) was included in the pipette solution, vNLPs were treated with the tetanus neurotoxin light chain (TeNT), or when empty NLPs (eNLPs) were used (Figure 4c). CDV competes with the NLP v-SNAREs for binding to the flipped t-SNAREs on the patch surface, and TeNT cleaves VAMP2. When we used VAMP2-4X, pore nucleation rate was not significantly different than for any of the other negative controls (Figure 4c). Because this construct allows efficient docking (Figure 2c) but is fusion incompetent (Figure 2b,d) (Wu et al., 2016; Krishnakumar et al., 2013), this result indicates that ApoE does not induce pores even when kept in close proximity to the target membrane. Collectively, these observations indicate that, similar to their smaller ND courterparts (Wu et al., 2016), NLPs fuse with liposomes or cell membranes in a strictly SNARE-dependent manner.
Combining data from 64 current bursts, we obtained distributions for vNLP8-tCell fusion pores as shown in Figure 4e–f and Figure 4—figure supplement 1. The number of pore flickers and burst durations were well described by geometric and exponential distributions, respectively, with flickers per burst and s (mean ± S.E.M.), as would be expected for discrete transitions between open, transiently blocked, and closed states (Sakmann and Neher, 2009) (Figure 4—figure supplement 1).
Conductances in the open-state and corresponding radii were broadly distributed (Materials and methods and Figure 4e,f), with mean pS (S.E.M. = pS), and nm (S.E.M nm), respectively. Surprisingly, these values were significantly less than the maximum possible value predicted from NLP dimensions (Figure 4b). This suggested a substantial inherent resistance to pore expansion, independent of the constraints imposed by the NLP dimensions. To quantify the resistance, we computed the apparent pore free energy from the distribution of pore radii, . This suggested that ~2 kT energy was required for every 1 nm increase in pore radius above the most likely value nm (Figure 4f,g).
A few SNARE complexes are sufficient to create a fusion pore, but many more are needed to dilate it
We then varied the number of SNAREs, and found that just a few SNARE complexes are sufficient to create a fusion pore, but many more are needed to dilate it. We repeated the measurements shown in Figure 4 using NLPs loaded with total v-SNARE copy numbers ranging from 1 (vNLP1) to ~30 (vNLP30, ~15 copies per face) (Figure 5). Pore nucleation required at least two v-SNAREs per NLP face and maximal nucleation rates were reached at around the same value (Figure 5a). By contrast, when four v-SNAREs per NLP face were present, the pore conductances (Figure 5b) and radii (Figure 5—figure supplement 1b) were significantly larger than the SNARE-free values and increased dramatically as the copy number per NLP face reached 15. Conductance fluctuations about the mean increased even more sharply (Figure 5—figure supplement 1a), while burst lifetimes and pore open probability showed a more gradual increase (Figure 5—figure supplement 1c,d). Thus, different numbers of SNARE complexes cooperate at the distinct stages of fusion pore nucleation and pore dilation.
-
Figure 5—source data 1
- https://doi.org/10.7554/eLife.22964.013
Is the increase in the mean pore conductance as the SNARE copy numbers are increased (Figure 5b) due to the appearance of multiple small pores per NLP or due to an increase in the mean size of a single pore? The latter is much more likely, for the following reasons. First, a probe that cannot pass through small pores becomes permeant to pores when large copy number vNLPs are employed (Figure 5—figure supplement 4). If multiple small pores were present when vNLP30 are used, then the probe should be equally impermeant. The probe employed was N-methyl-D-glucamine (NMDG+), a large ion of nm in size without its hydration shell (Melikov et al., 2001), which replaced sodium in the pipette solution. Conductance was low when ~15 nm MSP nanodiscs with eight copies of v-SNAREs (vMSP8) were used (Wu et al., 2016), but not affected when ~23 nm NLPs bearing 30 v-SNAREs were employed (vNLP30). These results are consistent with those of Bello et al. (2016), who showed that progressively larger cargo could be released from t-SNARE liposomes during fusion with vNLPs as the v-SNARE copies per NLP was increased. Second, conductance of small pores in a single NLP would be additive, giving total conductance equal to , where is the mean open-pore conductance of a small pore. Doubling the SNARE copies would presumably at most double , and by consequence, total conductance. The fact that we find faster than linear increase in mean pore conductance (Figure 5b) is consistent with each NLP bearing a single pore whose size increases with increasing SNARE copies. Third, if multiple small pores occurred per NLP, this should be evident in the distribution of point-by-point conductance values, with peaks at , where . Instead, for the distribution of mean for vNLP30 we find a peak at ~300 pS, and a broad peak at ~3–14 nS (Figure 6b). If the typical small pore has 300 pS conductance, then to have ~6 nS (typical large conductance), there would have to be ~20 small pores per NLP. It is hard to imagine that this many pores could coexist in this small area. Finally, unless the multiple pores occurred simultaneously, we would also find that the fusion rate increases with copy number. Instead, the rate saturates at around two copies (Figure 5a). In conclusion, although we cannot rule out that, very occasionally, a small number of pores may simultaneously appear in a single NLP, all the evidence suggests that this cannot be very common.
Previous reports suggested that pore nucleation is promoted by the assembly of the v- and t-SNARE transmembrane domains (TMDs) (Wu et al., 2016; Shi et al., 2012; Fdez et al., 2010). To test whether pore dilation also required the TMDs, we replaced the v-SNARE TMDs with lipid anchors (Shi et al., 2012). We used long-chain anchors that span both leaflets of the bilayer because previous work suggested that lipid anchors spanning a single leaflet are not efficient in inducing full fusion (McNew et al., 2000; Chang et al., 2016). Lipid-anchoring VAMP2 into NLPs significantly reduced pore nucleation frequency and increased the mean burst duration (Figure 5—figure supplement 2a,b,e), consistent with previous work using smaller MSP NDs (Wu et al., 2016) and with reduced overall fusion efficiency reported for lipid-anchored VAMP2 (Shi et al., 2012; Chang et al., 2016). Importantly however, fusion pores induced by lipid-anchored v-SNAREs displayed the same trends as their intact counterparts: as VAMP2-C45 copy number increased, so did mean conductance, fluctuations, burst lifetimes, and pore radii, but the pore open probability during a burst varied little (Figure 5—figure supplement 2). These results suggest that specific interactions between v-SNARE and t-SNARE TMDs are not critical for cooperative pore dilation by SNARE proteins.
The target membrane during exocytosis (where the t-SNAREs reside) is the inner leaflet of the plasma membrane, which is rich in acidic phospholipids. By contrast, in our vNLP-tCell fusion assay, the target membrane is the outer leaflet of the plasma membrane, which is largely devoid of negatively charged lipids. In general, a limitation of our system is that the lipid composition of the outer leaflet of the t-SNARE-presenting cell differs substantially from that of the plasma membrane inner leaflet, and lipid composition can play a key role in fusion. To test whether the target membrane composition affected fusion, we swapped the locations of the v- and t-SNAREs and fused flipped v-SNARE cells with t-SNARE NLPs. This allowed us to have a better mimic of the inner plasma membrane leaflet composition on the target membrane (now the tNLP membrane). This swap resulted in similar fusion rates and pore properties for two different SNARE copy numbers per NLP (Figure 5—figure supplement 3), suggesting that fusion mediated by SNAREs alone may not be very sensitive to target membrane composition within a certain range (Stratton et al., 2016).
Interestingly, the increase in mean pore conductance as v-SNARE copy numbers are increased does not occur homogeneously across all pores. For vNLPs bearing four or more total copies of v-SNAREs, we found two types of fusion pores. Most had small mean conductance nS, but with increasing SNARE load, an increasing fraction of pores had much larger conductances of a few nS (Figure 6a). By contrast, for NLPs that contained no copies or just one copy of VAMP2, the pores that occasionally occurred all had small mean conductance nS. The distribution of mean conductances for individual pores revealed two components for vNLP30 (Figure 6b), with a sharp peak at ~300 pS (71% of total mass) and a much broader population centered at 7.21 nS, separated at ~1 nS. Conductance fluctuations increased sharply for larger pores with nS (Figure 6—figure supplement 1), indicating a change in behavior above this threshold. Thus, multiple criteria indicated ~1 nS as a cut-off that separated small and large pores. We applied this cut-off to all NLPs tested, and clustered pore conductances in each NLP group into low ( nS) and high conductance ( nS) states (Figure 6c). The occurrence of high conductance pores increased with increasing SNARE copy number, suggesting that dilation of the pores to >1 nS (corresponding to nm) is facilitated by SNARE crowding.
We defined the pore dilation probability, , as the fraction of pores in the high-conductance state for a given SNARE copy number (Figure 6c), and plotted as a function of v-SNARE copies per face (Figure 6d). Even at the maximum SNARE load of ~15 copies per face, was ~0.30, far from saturating. To estimate how many SNARE complexes would be required to reach saturation, we assumed that the ratio between the probabilities of high and low conductance states, is equal to a Boltzmann factor , with the difference between the energy levels of the two conductance states, and thermal energy. Making the simplest assumption that , where is the number of SNARE complexes involved and is the copy number that would make , we found 19.3 (Figure 6d). Thus, pore dilation with 0 would require complexes (Figure 6d, open black circle), about 10-fold more than is required for nucleation (Figure 5a).
SNARE crowding generates entropic forces that drive pore expansion
The characteristic fusion pore free-energy function progressively softened as the SNARE copy number increased (Figure 7a). The minimum at nm did not shift, but for larger pore radii, the slope decreased and the profile broadened. These free-energy profiles quantify how fusion-pore dilation is driven by SNARE proteins. For example, an energy ~6 kT is required to expand the one-SNARE fusion pore from its preferred radius of ~0.5 nm to a 3-fold larger pore, showing that such an expansion is unlikely to occur spontaneously. On the other hand, with four SNAREs per face, the same expansion requires only ~3 kT, and only ~2 kT with 15 SNAREs (Figure 7a), bringing the expansion within reach of spontaneous fluctuations. The broad and shallow profile suggests that a fully developed pore would then result, if the NLP were replaced by a vesicle as in the physiological setting (Figure 7a). These effects can equivalently be phrased in terms of force: ~24 pN opposes pore expansion with one SNARE, but this is lowered to ~8 pN and ~5 pN by 4 and 15 SNAREs per face, respectively.
To help to elucidate the molecular mechanism underlying SNARE-mediated pore dilation, we developed a coarse-grained mathematical model that assumed that the bilayer-SNARE system is equilibrated, consistent with the long-lived current bursts, so that is then the true thermodynamic free energy (see Materials and methods and Table 1 for model parameters). In our model, for a given pore size, the free energy represents an average over possible fusion pore heights, pore shapes and SNARE complex configurations. SNAREs can be fully zippered at the fusion pore waist, or they may unzipper and roam (Figure 7—figure supplement 1).
The model free energies reproduced experimental profiles with physiologically realistic parameters (Figure 7b). Protein-free pores resisted expansion because a bigger pore has greater area: a 1 nm increase in pore radius required ~3.0 kT work against membrane tension and increased membrane bending energy by ~2.4 kT. Thus, a net force ~22 pN resists pore expansion, close to the experimental value of ~24 pN (Figure 7a). When SNAREs were present, zippering of SNARE linker domains (Gao et al., 2012) and TMDs drove several SNAREs to fully assemble at the fusion pore waist, where crowding generated an entropic pore expansion force (Figure 7c). Bigger pores were associated with more zippered SNAREs at the waist (Figure 7—figure supplement 1). With 15 SNAREs per face, the entropic expansion force of 17 pN was within 5 pN of the 22 pN resistance.
Consistent with our experiments using lipid-anchored v-SNAREs, when we ran the model with lowered excluded volume among zippered SNAREs to mimic the lipid anchor, the best fit total zippering energy was reduced by (see Materials and methods). This suggests that the driving force for SNARE zippering that underlies pore expansion does not rely on putative v- and t-SNARE TMD interactions, but has a significant contribution from non-specific interactions that favor the alignment of membrane anchors.
Discussion
In summary, we find that a few SNAREs can nucleate a fusion pore, consistent with previous findings (Bao et al., 2016; Shi et al., 2012; Sinha et al., 2011; Mohrmann et al., 2010; van den Bogaart et al., 2010; Karatekin et al., 2010), but the pore is highly unlikely to expand significantly without additional forces. Pore dilation is resisted by intrinsic bilayer properties (Chanturiya et al., 1997; Chizmadzhev et al., 1995; Jackson, 2009), but promoted by the action of many SNAREs that cooperatively exert expansion forces of entropic origin (Figure 7c).
In our study, pores fluctuated in size, and closed and opened (flickered) multiple times before resealing, as do exocytotic fusion pores that have been recorded from neuroendocrine cells or neurons (Lindau et al., 2003; Staal et al., 2004; Fulop et al., 2005; He et al., 2006; Klyachko and Jackson, 2002). Each such burst lasted several seconds on average. Confinement of the pore radius by the NLP scaffold to 7–8 nm probably contributed to this long lifetime. However, we suggest the lifetime also reflects the natural tendency of small pores to resist dilation, for the following reasons. First, the most likely pore radius, ~0.5 nm (Figures 4f and 7a), is much smaller than the maximum allowed by the NLP geometry, ~7–8 nm. Second, fusion pores connecting protein-free bilayers flicker for seconds, and do not dilate unless increased membrane tension is applied (Chanturiya et al., 1997). Third, pores lasting of the order of a second or longer have been documented during exocytosis using capacitance recordings made, for example, in beta cells secreting insulin (Hanna et al., 2009; MacDonald et al., 2006) or during synaptic vesicle fusion (He et al., 2006). Amperometry often reports shorter pore lifetimes when compared to capacitance measurements (Chang et al., 2015); it may underestimate pore lifetimes because no signal can be observed once all cargo is released. It is also possible that a pore can reseal after partial dilation giving rise to an amperometric spike, leading to an underestimation of pore lifetime based on the pre-spike foot feature alone (Mellander et al., 2012). Fourth, a recent FRET-based study suggested the existence of long-lived, narrow fusion pores during neuronal SNARE-driven fusion between surface-tethered liposomes, dilation of which was promoted by Synaptotagmin-Ca2+ and Complexin (Lai et al., 2013). Finally, various theoretical models suggested that small pores are metastable (Chizmadzhev et al., 1995; Jackson, 2009; Nanavati et al., 1992).
We measured a low fusion efficiency in our assay (4–5% of docked NLPs undergo lipid mixing within ~20 min (Figure 2—figure supplement 1), possibly due to the absence of factors known to be essential for exocytosis in our assay. Indeed, the low fusion efficiency in SNARE-only reconstitutions that lack other factors that are required in vivo for exocytosis, such as Munc13, is well documented (Bao et al., 2016; Bello et al., 2016; Weber et al., 1998; Liu et al., 2016; Hernandez et al., 2014; Diao et al., 2013). Another possibility is that the NLPs that actually fuse may be biased toward higher copy numbers of proteins, if higher copy number corresponded to higher fusion rates. While we cannot categorically exclude this possibility, we think it unlikely. First, the fusion rates that we measured were statistically indistinguishable for copy numbers four (vNLP4, vNLP8, vNLP15 and vNLP30, Figure 5a). Thus, there is no evidence of a bias due to differential fusion rates. Second, even if there were such a bias, our results would still correctly report the general trend of pore properties versus copy number, as can be seen from the following argument. Assume a Poisson distribution for the copy number up to a maximum of 30, the maximum attainable value in our experiments (presumably a packing constraint). Then, for large mean copy numbers, this distribution has a small width, so that even if the ~5% fused fraction corresponds to the tail of this distribution, the copy numbers involved will not be much greater than the mean value. Thus, the typical copy number of the NLPs whose pore properties are measured would still be an increasing function of the mean value.
A wide range of SNARE copy number requirements for fusion have been reported in the literature (Bao et al., 2016; Shi et al., 2012; Sinha et al., 2011; Mohrmann et al., 2010; van den Bogaart et al., 2010; Karatekin et al., 2010; Hernandez et al., 2014; Domanska et al., 2010; Montecucco et al., 2005), depending on the system studied and the read-out used for fusion. Most studies concluded that only a few copies of neuronal SNAREs are sufficient for calcium-triggered exocytosis and fusion of small liposomes (Bao et al., 2016; Shi et al., 2012; Sinha et al., 2011; Mohrmann et al., 2010; van den Bogaart et al., 2010; Karatekin et al., 2010; Domanska et al., 2010). Despite this, the average synaptic vesicle carries 70 v-SNARE copies (Takamori et al., 2006) and at least as many t-SNAREs are clustered at plasma membrane docking and fusion sites in neuroendocrine cells (Knowles et al., 2010). Our results provide a rationalization for this situation, as they suggest that reliable pore dilation may require the engagement of many SNARE complexes. As the demands for SNARE cooperativity may be different at different stages of the fusion reaction, interpretation of copy number requirements should be made with caution. Methods that rely on lipid mixing or on the exchange of small ions through pores (e.g. capacitance or pH sensing) may measure the requirement for the opening of small fusion pores, which may differ substantially from the requirements for pore dilation.
The action of SNAREs is highly regulated by other proteins during neurotransmitter or hormone release (Südhof and Rothman, 2009). In addition to manipulations of SNAREs (Han et al., 2004; Fang et al., 2008; Kesavan et al., 2007), mutations in Munc18 (Jorgacevski et al., 2011), Synaptotagmin (Wang et al., 2001, 2003a, 2003b), and complexin (Dhara et al., 2014) affect fusion pore properties, linking these proteins to pore dynamics. Thus, one must be cautious when extending our SNARE-only results to physiological neurotransmitter and hormone release. Nevertheless, the concept of the promotion of pore dilation by protein crowding is a general principle that may hold qualitatively in the presence of additional components of the physiological fusion machinery. Indeed, a previous study indicated that the availability of SNAREs affects neurotransmitter release kinetics in neurons (Acuna et al., 2014), while another suggested release occurred faster at sites with more t-SNAREs (Zhao et al., 2013). Thus, we tentatively suggest that some proteins may exert their exocytotic regulatory function by organizing SNARE complexes around the fusion site and thereby controlling the number that participate, or by sequestering SNAREs to limit that number. Given the steep dilation probability curve (Figure 6d), our results suggest a high sensitivity in the balance between transient versus full fusion.
Materials and methods
Stable flipped SNARE and wild-type HeLa cell culture
Request a detailed protocolHeLa cell lines stably co-expressing flipped v-SNAREs (flipped VAMP2 and cytosolic DsRed2-nes, ‘vCells’) and t-SNAREs (flipped Syntaxin-1 and flipped SNAP-25 and the nuclear fluorescent marker CFP-nls, ‘tCells’) were generated in the Rothman laboratory as described (Giraudo et al., 2006). The cells were a generous gift from the Rothman laboratory. The cells were maintained in DMEM (4500 mg/L glucose, L-glutamine, sodium pyruvate and sodium bicarbonate) and 10% (v/v) fetal calf serum at 37°C. A new aliquot of cryopreserved cells was thawed after at most three weeks of cell culture and cultured at least five days before data acquisition.
The flipped t-SNARE HeLa cells were tested by PCR Mycoplasma Test Kit I/C (cat. No. PK-CA91-1048, Promo Kine, Heidelberg, Germany), which showed contamination (Figure 4—figure supplement 3A). Because our assay fuses discs to the surface of these cells, no effect of mycoplasma contamination is expected. Indeed the controls in Figures 2–4 indicate that fusion is SNARE-driven. Nevertheless, we tested any possible impact of mycoplasma contamination on our results by repeating some of our single-pore measurements with cells treated with an antimycoplasma reagent (Plasmocin, cat. code ant-mpt, InvivoGen, California, USA). Fusion rates and pore properties were indistinguishable when untreated or treated cells were used (Figure 4—figure supplement 3, B–E), suggesting that mycoplasma contamination does not affect fusion with NLPs.
Plasmids, protein expression and purification
Request a detailed protocolExpression and purification of the t-SNARE complex used in vND-tSUV fusion experiments is described in Parlati et al. (1999). The cytoplasmic domain of VAMP2 (CDV) was produced using a method that is similar to an earlier protocol (Weber et al., 1998), except that a SUMO vector was used. We followed the methods of Giraudo et al. (2006) for expression and purification of the tetanus neurotoxin light chain, TeNT. VAMP2 proteins were expressed and purified as described earlier by Shi et al. (2012). To produce full-length WT VAMP2, we used the plasmid pET-SUMO-VAMP2 (Shi et al., 2012). To produce lipid-anchored VAMP2, we followed the methods of McNew et al. (2000) and Shi et al. (2012). We first used a previously described construct (McNew et al., 2000) to produce VAMP95Cys containing the entire cytoplasmic domain of VAMP2 (residues 1–95) with a C-terminal cysteine residue. We then coupled this protein to maleimidopropionic acid solanesyl ester (maleimide-C45), produced as previously described (Shi et al., 2012). For producing MSP NDs, we used the vector pET28-MSP1E3D1 (Addgene, Cambridge, MA) to express and purify MSP1E3D1 as described previously (Ritchie et al., 2009), except that we cleaved the MSP proteins directly off the column by TEV protease overnight at 4 °C (Wu et al., 2016).
Plasmid pET32a-Trx-His6X-ApoE422K, which we used to express the N-terminal 22 kDa fragment of apolipoprotein E4 (residues 1–199), was kindly provided by Dr Nicholas Fischer, Lawrence Livermore National Laboratory, CA (Morrow et al., 1999; Blanchette et al., 2008). ApoE422K was expressed and purified as previously described (Morrow et al., 1999) with the following modifications. The His6-ApoE422K was cleaved off the Ni-NTA beads (Qiagen, Germantown, MD) using 100U of Thrombin at 4°C overnight. The protein was eluted in 25 mM HEPES, 140 KCl, pH 7.4 buffer containing 1% octylglucoside (OG), and was functional for up to 4 weeks when stored at 4°C. Protein concentrations were determined using the Bradford assay (Bio-Rad, Hercules, CA) with bovine serum albumin as standard.
Characterization of nanolipoprotein particles (NLP)
Request a detailed protocolDetails are provided by Bello et al. (2016). Briefly, nanolipoprotein particles containing VAMP2 (vNLP) were produced using a modified version of the established protocol to generate SNARE-nanodiscs (Shi et al., 2012, 2013). A palmitoyl-2-oleoylphosphatidylcholine (POPC): 1,2-dioleoyl phosphatidylserine (DOPS) = 85:15 molar ratio lipid mixture (Avanti Polar Lipids, Alabaster, AL) was dried under nitrogen flow, followed by vacuum for 1 hr. The lipid film was re-suspended in 25 mM HEPES, pH 7.4, 140 mM KCl, buffer with 1% OG supplemented by the desired amount of VAMP2. The mixture was vortexed at room temperature (RT) for 1 hr followed by the addition of ApoE422K and vortexed for another hour at RT. The ApoE422K:VAMP2:lipid ratio was varied to tune the v-SNARE copy number per NLP as 1:0.2:180 (1 VAMP2, ‘vNLP1’); 1:1:180 (four copies, ‘vNLP4’); 1:2:180 (eight copies, ‘vNLP8’); 1:4:180 (15 copies, ‘vNLP15’) and 1:8:180 (30 copies, ‘vNLP30’). NLPs containing 1, 4 and 15 copies of VAMP-C45 were obtained using a similar approach. Excess detergent was removed using SM-2 bio-beads (Bio-Rad) overnight at 4°C with constant mixing. The assembled v-NLPs were separated from free proteins and lipids via gel filtration on a Superose six column (Figure 1a). Samples were concentrated using Amicon Ultra (50 KDa cutoff) centrifugal filter units, and analyzed by SDS-PAGE with Coomassie staining (Figure 1b). The number of VAMP2 copies per disc was determined by the VAMP2-to-ApoE ratio by densitometry using ImageJ (NIH). The number of ApoE copies per disc was estimated (Bello et al., 2016) using the calibration of disc size vs the number of ApoE copies previously reported (Blanchette et al., 2008). The size distribution of the v-NLPs was determined for every batch of production using transmission electron microscopy. To do this, the NLP discs were diluted (1:50), mounted onto carbon-coated 400 mesh copper electron microscopy grids, negatively stained with 2% uranyl acetate (w/v) solution, and subsequently examined in an FEI Tecnai-12 electron microscope operated at 120 kV. Micrographs of the specimen were taken on a Gatan Ultrascan4000 CCD camera at a magnification of 42,000. Typical micrographs and a size distribution are shown in Figure 1c,d. The size of the NLP discs with 1:180 ApoE422K: lipid ratio was typically 24 ± 2 nm (100–200 NLP discs were analysed for every production batch). Representative size distributions are shown as box plots for the conditions tested in Figure 1e. At least three independent batches of NLPs were used per condition. tNLPs were produced in a similar fashion, using a t-SNARE:ApoE:lipid ratio of 0.8:1:180 and 3:1:180 for tNLP4 (four copies of t-SNARE complex Stx/SN25 total per NLP) and tNLP 15 (15 total copies of t-SNAREs per NLP) samples, respectively.
Bulk fusion of NLPs with t-SNARE liposomes
Request a detailed protocolWe used a previously established assay (Bello et al., 2016; Shi et al., 2012) to monitor the release of calcium from t-SNARE-reconstituted small unilamellar vesicles (t-SUVs) as they fused with discs loaded with v-SNAREs. 40 μl t-SUVs entrapping 50 mM calcium were mixed with 5 μl vMSP NDs (prepared as described previously [Shi et al., 2012]) or with 10 μl of vNLP-discs in a buffer containing 2 μM of the calcium-sensitive dye mag-fluo-4 (Invitrogen, Carlsbad, CA). The mixture was loaded into a 96-well plate, and the mag-fluo-4 fluorescence (λex= 480 nm, λem=520 nm, 515 cutoff) was recorded by a SpectraMax M5 plate reader (Molecular Devices, Sunnyvale, CA). After 60 min, 15 μl of 5% dodecylmaltoside was added and the mixture was incubated for an additional 20 min to release all remaining entrapped calcium and thus to establish the maximum mag-fluo4 signal. Fusion is reported as percent of maximum fluorescence in Figure 1—figure supplement 1.
Single-cell lipid mixing and calcium-influx assays
Request a detailed protocolThese assays were carried out as described in Wu et al. (2016).Briefly, for lipid mixing, tCells were plated in 35 mm poly-D-lysine-coated glass bottom dishes (MatTek Corporation, MA, USA) and vNLP8s were prepared as described above, except that one mole % each of 1,1'-dioctadecyl-3,3,3',3'-tetramethylindocarbocyanine perchlorate (DiI, cat. no. D282, Molecular Probes, Eugene, OR) and 1,1'-dioctadecyl-3,3,3',3'-tetramethylindodicarbocyanine perchlorate (DiD, cat. no. D307, Molecular Probes) fluorescent lipid labels were included in the lipid composition. For each reaction, 15 µl of vNLP8 was added onto tCells (final NLP lipid concentration was ~54 µM) and incubated for 30 min at 4°C, a temperature at which SNARE complexes assemble but cannot drive fusion (Weber et al., 1998). Fusion was started by raising the temperature to 37°C and monitored by the dequenching of the DiI fluorescence using confocal microscopy. As controls, empty nanodiscs (eNLP) or NLPs bearing VAMP2-4X (harboring the mutations L70D, A74R, A81D and L84D) were used (Figure 2b). To estimate the extent of lipid mixing, at the end of some experiments, DiD fluorescence was completely bleached using direct excitation at 647 nm with 100% laser power. This resulted in the maximum possible donor (DiI) intensity, . We then rescaled the donor fluorescence values to obtain the fraction of maximum DiI fluorescence: , where is the minimum at the beginning of acquision.
To assess lipid mixing, we also used an alternative protocol that avoided the cold incubation step; this protocol could not be used for time-course measurements because of the high background resulting from excess NLPs. The same amount of vNLP8 as above was added to tCells at 37°C. After 15 min incubation, excess NLPs were washed, and DiI, DiD and CFP fluorescence levels were acquired using confocal microscopy (Figure 2c,d).
To measure the influx of calcium through fusion pores, tCells were loaded with 5 μM of Fluo-4 AM (Life Technologies, Carlsbad, CA), a cell-permeant calcium-sensitive fluorescent dye, as previously described (Wu et al., 2016). After washing to remove dye that was not taken in by cells, vNLP (15 µl) were added to tCells at 37°C and the influx of calcium was tracked by imaging of Fluo-4 fluorescence using a confocal microscope (Wu et al., 2016).
Electrophysiology
Request a detailed protocolDetails are given in Wu et al. (2016). Briefly, flipped t-SNARE HeLa cells (tCells) were cultured in 3 cm dishes. For recordings, a dish was placed in a temperature-controlled holder (TC-202A by Harvard Apparatus (Holliston, MA), or Thermo Plate by Tokai Hit (Shizuoka-ken, Japan)) set at 37°C. Cells were visualized with an inverted microscope (Olympus IX71, Olympus Corp., Waltham, MA) using an Andor DU-885K EMCCD camera controlled by Solis software (Andor, South Windsor, CT). Recording pipettes (borosilicate glass, BF 150-86-10, Sutter Instruments, Novato, CA) were pulled using a model P-1000 pipette puller (Sutter Instruments) and polished using a micro-forge (MF-830, Narishige, Tokyo, Japan). Pipette resistances were 5–10 MΩ in NaCl-based solution. The bathing medium contained: 125 mM NaCl, 4 mM KCl, 2 mM CaCl2, 1 mM MgCl2, and 10 mM HEPES, (pH adjusted to 7.2 using NaOH) for the cell-attached recordings. 10 mM glucose was added to the medium before use. All voltage- and current-clamp recordings were made using a HEKA EPC10 Double USB amplifier (HEKA Elektronik Dr. Schulze GmbH, Lambrecht/Pfalz, Germany), controlled by Patchmaster software (HEKA). Currents were digitized at 20 kHz and filtered at 3 kHz.
To measure SNARE-mediated single fusion pore currents in the cell-attached mode (Yang and Sigworth, 1998), electrodes were filled with the pipette solution composed of 125 mM NaCl, 4 mM KCl, 10mM HEPES, 13 mM or 26 mM tetraethylammonium-Cl (TEA-Cl, K+-channel antagonist), adjusted to pH 7.2 using NaOH. This solution had resistivity of Ohm.m, measured using a conductivity cell (DuraProbe, Orion Versa Star, Thermo Scientific). For experiments designed to test the presence of multiple pores connecting large copy number vNLPs to tCells (Figure 5—figure supplement 4), electrodes were filled with a solution containing: 129 mM N-methyl-d-glucamine (NMDG), 10 mM HEPES, 26 mM TEA-Cl, pH adjusted to 7.2 using HCl, resistivity 0.88 Ohm.m, 305 mOsm. The pipette tip was initially filled with 1 μl of NLP-free buffer and back-filled with vNLPs suspended in the same buffer (final [vNLP] = 100 nM, 120 μM lipids). This allowed the establishment of a tight seal ( GOhm) with high success rate, as well as the recording of a stable baseline before the vNLPs diffused to the membrane patch and started fusing with it 2–18 min later. Such a back-filling strategy is typically used in perforated patch measurements (Sakmann and Neher, 2009). All cell-attached recordings were performed using a holding potential of mV relative to bath. With a cell resting membrane potential of mV (mean ± S.D., n = 36), this provided 16 mV driving force across the patch membrane.
Analysis of fusion pore data
Request a detailed protocolThe analysis of fusion pores is described in detail in Wu et al. (2016).Briefly, we developed an interactive graphical user interface in Matlab to help to identify, crop and process single fusion pore currents. Traces were exported from Patchmaster (HEKA Electronik) to Matlab (Mathworks) and low-pass filtered (280 Hz cutoff); frequencies that were the result of line voltage were removed using notch filtering. Zero phase shift digital filtering algorithms (Matlab Signal Processing Toolbox function filtfilt) were employed to prevent signal distortion. Filtered traces were averaged in blocks of 80 points (125 Hz final bandwidth) to achieve rms baseline noise pA. Currents for which pA for at least 250 ms were accepted as fusion pore current bursts. During a burst, rapidly fluctuating currents often returned to baseline multiple times, i.e. pores flickered. To quantify pore flickering, we defined currents pA and lasting ms (15 points) as open pores and currents not meeting these criteria as closed. For a given burst, the number of open periods was equal to the number of flickers, . To estimate the fusion rate for each recording (i.e. the rate at which current bursts appeared), we counted the number of current bursts that fit the set criteria (current amplitude >2 pA for at least 250 ms) and divided this number by the duration of the recording. Rates from different records (patches) were averaged for each condition. We also refer to this rate as the pore nucleation rate. Periods during which the baseline was not stable were excluded from this analysis. Many recordings ended with what seemed to be currents from overlapping fusion pores. Such end-of-record currents were also excluded because they could also be attributed to a loose seal. Thus, the fusion rates that we report may underestimate the true rates, especially for conditions in which fusion activity was high. For distributions of conductances and radii, we used pore open-state values, denoted by the subscript ‘po’. For Figure 4e,f, we first computed the probability density functions (PDFs) for individual pores using a fixed bin width for all, then averaged these to give equal weight to all pores. All distribution fitting was performed using Matlab Statistics Toolbox functions fitdist or mle, using maximum likelihood estimation. Open-pore conductance values were used point-by-point to estimate the open-pore radii, by approximating the pore as a cylinder and using the expression (Hille, 2001) , where is the resistivity of the solution, nm is the length of the cylinder, and is the open-pore conductance. For assessing statistical significance when comparing sample means, we used the two-sample t-test or the nonparametric two-sample Kolmogorov-Smirnov test (ttest2 or kstest2, Matlab Statistics Toolbox), as indicated in the figure legends. We considered each single-pore measurement to be a biological replicate.
For clustering average single-pore conductances for vNLP30 measurements, we used a two-component Gaussian mixture model (Figure 6b) that indicated a boundary between the two components at ~1 nS. Applying this cutoff to all vNLP samples, we produced the boxplot in Figure 6c, where the central red line on each box marks the median, the edges of the box are the 25th and 75th percentiles, and the whiskers extend from q3 +1.5(q3 – q1) to q1 – 1.5(q3 – q1), where q1 and q3 are the 25th and 75th percentiles, respectively. For a given v-SNARE copy number per NLP face, , we defined the probability of achieving a high-conductance pore ( nS) as the fraction of high conductance pores observed for that copy number. For example, only 3 out of 64 pores were large conductance for vNLP8, which had four copies per face (Figure 6c), hence . We plotted as a function of v-SNARE copy number per NLP face in Figure 6d. In Figure 7a, to estimate the energy profiles of fusion pores for a given v-SNARE copy number, we first calculated the probability density function for open-pore radii as in Figure 4f. The probability that the radius is between and is the density at that bin , where is the bin width. We estimated the energy of a pore with radius by , where is an arbitrary constant.
Mathematical model of the fusion pore between a nanodisc and planar membrane in the presence of SNAREs
Membrane free energy
Request a detailed protocolWe modeled the fusion pore as having a toroidal shape formed between a nanolipoprotein particle (NLP) modeled as a planar bilayer of diameter and the tCell membrane modeled as an infinite planar bilayer, both of which are at a constant membrane tension (Figure 7c). This toroidal assumption is similar to that in previous theoretical studies that assumed a toroidal shape of the fusion pore (Chizmadzhev et al., 2000; Jackson, 2010; Kozlov et al., 1989). Experimental studies also observed an hourglass-shaped fusion pore that could be considered approximately toroidal (Curran et al., 1993; Haluska et al., 2006).
The fusion pore is parametrized by the toroidal shape parameters: the radius of the toroid , which corresponds to the fusion pore radius , and the separation of the membranes at the edge of the NLP . The fusion pore is completely toroidal at small pore sizes. The free energy of the fusion pore is calculated using the Helfrich energy form, as was used in previous studies (Chizmadzhev et al., 2000; Jackson, 2010; Kozlov et al., 1989)
Here, the energy due to bending is given by
where and C are the bending modulus and the mean curvature of the membrane, respectively. The energy expended to add membrane area due to membrane tension is the second term, where is the change in total membrane area due to pore formation, given by . Here, is the area of the fusion pore. is the area of both rims of the fusion pore, which is the area that has to be removed from the infinite tCell membrane and the NLP membrane to form the pore. We evaluated all integrals and all areas over the midplane of the membrane forming the pore to give:
where and . Values of in previous studies range from 10–50 (Cohen and Melikyan, 2004; Brochard and Lennon, 1975; Khelashvili et al., 2013; Marsh, 2006). A commonly used value is (Jackson, 2009), which we used here (Table 1). γ was obtained as a best-fit parameter (Table 1).
Free energy due to twisting of the ApoE proteins that line the boundary of the ND
Request a detailed protocolOwing to the finite size of the NLP, toroidal states are not possible for large pores with and that constitute a sizeable fraction of the NLP diameter, . These shapes are partially toroidal and come into existence when , where is the membrane thickness. The ApoE proteins that line the NLP boundary need to be rotated through an angle to form these shapes. An example of one such shape is the right-hand side of Figure 7c.
We assumed these proteins exert a constant torque per unit length of the NLP boundary to resist this rotation. The ApoE proteins exert no torque in the completely toroidal states as vanishes for these states. The free energy of these proteins is
Thus, the ApoE proteins resist pore expansion as increases with the size of the pore. We obtained as a fitting parameter and was measured in this study (values in Table 1). For further details about the partially toroidal states, please see the penultimate subheading ‘Description of partially toroidal states and calculation of membrane free energy’.
Free energy contribution from short-ranged steric-hydration forces
Request a detailed protocolSteric-hydration repulsion between membranes is prominent at small membrane separation. Experimentally measured steric-hydration pressures between planar membranes of separation are of the form (Rand and Parsegian, 1989). Values of and have been measured before for several membrane compositions; is within 0.1–0.3 nm (Rand and Parsegian, 1989). As the pore sizes over which these effects are appreciable () are very small compared with the NLP diameter , only toroidal states are considered for this calculation. We obtained from previous studies and as a best-fit parameter (Table 1).
The steric-hydration forces act in two orthogonal directions on the membranes comprising the fusion pore: to increase the pore radius () and the separation between membranes (). The sum of the work done by these two forces gives the free energy of the steric-hydration interaction
where is the effective pore height, i.e. the height of the section of the pore that contributes substantially to the steric-hydration interaction. To obtain the free energy contribution of the steric-hydration forces, we now calculate the work done by these forces to assemble the fusion pore. We first observe that the work done () to bring two patches of membranes of area to a separation from a large distance apart is
The first term in Equation 7 is the work done to separate the planar part of the membranes to a distance . As the pore area where these forces are relevant is very small compared with the NLP area , we set the area of the planar region in Equation 8 to obtain the first term. The second term in Equation 7 is the work done to separate the membranes to form a pore of diameter . To calculate this, we can imagine that the pore is a cylinder of diameter and height , where , since the change in pore diameter over the height is negligible. Here, is the thickness of the bilayer and is the height over which the cross-sectional diameter of the toroidal shape increases from to .
To obtain , we consider the inner surface of the toroidal pore. This shape is formed by revolving the semicircle given by where in the XZ plane about the Z axis. Here, is the radius of the pore and is the radius of the semicircle, which is also equal to half of the maximum separation between the heads of the monolayers that line the inner surface of the pore, as can be seen from Figure 7c. is that value of at which the cross-sectional radius of the pore increases to . Thus, is obtained by solving the equation of the semicircle to first order in where .
To calculate the steric-hydration contribution from a pore of size and height , one need only consider the pore over that height at which the cross-sectional diameter of the toroidal shape increases from to , since λ is the range of the steric-hydration force. Thus, the area of the cylinder is , and using Equation 8, the work done to set up the pore is , giving the second term of Equation 7.
Free energy contribution from SNARE proteins
Request a detailed protocolIn this section, we calculate the free energy due to the SNARE proteins. We fix the total number of v-SNAREs. Out of these v-SNAREs, are fully zippered and are partially zippered. Only the TMDs and the linker regions of these partially zippered SNAREs are unzippered. For each , we allowed and to vary from 0 to to obtain an equilibrium distribution for both. We calculated the free energy of the SNAREs for , which correspond to the total number of v-SNAREs per NLP face used experimentally (Figure 7a). To match with experiment, we used the assumption that only half of the total number of SNAREs present in the NLP would be present on the side of the NLP that faces the tCell, and that all of these SNAREs would be available to participate in fusion pore expansion.
We assumed that the fully zippered SNAREs form a ring at the waist of the fusion pore (Figure 7c). Their free energy is
The first term is the positional entropy of the zippered TMDs whose diameter is , which we measured by fitting a cylinder to the measured crystal structure (Stein et al., 2009). The second term is the energy released when a partially zippered SNAREpin completes its zippering. The third term is the orientational entropy associated with the zippered SNAREs. We assume that these are very stiff cylindrical rods. Due to their high stiffness, these rods can only explore a small solid angle (Table 1). We calculated this angle based on an atomistic molecular dynamics study of the t-SNARE TMD that shows that these domains explore angles of ~ around their equilibrium position in a bilayer (Knecht and Grubmüller, 2003). The SNAREs fluctuate about their equilibrium orientation, which we assume is the local normal to the membrane.
The free energy of the partially unzippered SNAREs is
The first term is the positional entropy of the TMDs. These partially zippered SNAREs are in a Y-shape with both unzippered TMD domains on the same side of the pore, either on the vNLP or the tCell membrane. The linker domains are also unzippered and this imparts flexibility to these SNAREs (Jahn and Scheller, 2006). This orientational freedom is given by the second term. These SNAREs can adopt all orientations in which they do not intersect with the membrane; this corresponds to a solid angle of steradian. We restricted these SNAREs to a circle of radius equal to that of the NLP, as this considerable orientational freedom is only available when the SNAREpin body is away from the fusion pore lumen.
Calculation of the total free energy as a function of pore size and number of SNAREs
Request a detailed protocolThe probability that a fusion pore accesses a radius in the presence of SNAREs is proportional to in equilibrium, as this is the Boltzmann distribution where is the total free energy. We assumed that the bilayer-SNARE system is equilibrated as the current bursts measured experimentally are long-lived. To calculate this free energy, we summed the Boltzmann factor of all states that comprise such a system:
Here, is the total free energy of one fusion pore state with zippered SNAREs and of membrane separation , given by
We performed the integration and the sum over all states in Equation 11 numerically in MATLAB.
Derivation of best-fit model parameters by fitting model-predicted free energy to experiment
Request a detailed protocolWe performed a numerical calculation using Equation 11 to obtain free -energy curves as a function of SNARE copy numbers and pore size. For the membrane parameters, it is best to fit to data from membranes with no SNAREs. We fit the fusion pore free energy predicted by the model with no SNAREs (setting in Equation 11) to the experimentally measured curve for one SNARE, assuming that such a pore behaves similarly to a protein-free pore. We first fit the model-predicted pore size at the minimum in the free energy to experiment by using the steric-hydration force scale as a best-fit parameter. We obtained the bending modulus from Jackson (2009) and fit the slope following the minimum in the region using the membrane tension as a best-fit parameter, as γ largely determines the slope beyond this minimum.
Using these parameters, we calculated the free energy versus pore radius in the presence of SNAREpins (Figure 7b). As SNAREpins are introduced, the model predicts that the minimum barely shifts, while the slope beyond the minimum decreases with increasing numbers of SNAREs by an amount depending on the zippering energy parameter . We selected a typical experimental curve (vNLP30) and fit the slope of the free energy from simulation to that measured from experiment for and SNAREs (as vNLP30 corresponds to 15 SNAREs per face) using as a fitting parameter, and obtained . This is higher than the zippering energy of the linker domains alone, which was measured to be ~5 kT by Gao et al. (2012).
Further increases in pore size cause increases in free energy, as the fusion pore shapes are partially toroidal and the twisting torque from the ApoE proteins at the NLP boundary resists further expansion. Thus, we fit the slope of the free energy curve in this region () measured from simulation to that measured experimentally for SNAREs to obtain the torque per unit length pN as a fitting parameter.
To understand how lipid-anchored VAMP2 affects pore dilation, we reduced the size of the zippered TMDs by 50% as the zippered SNARE complex lacks the v-SNARE TMD. This is the maximum possible reduction in the excluded volume of the zippered SNAREs. We then varied the zippering energy to ensure that the model-predicted mean pore size was invariant with respect to this reduction, consistent with the invariance in mean pore size observed in experiments (Figure 5—figure supplement 2f). The best-fit zippering energy is 0.43 kT lower than the best-fit value obtained when both TMDs are present. This is an upper bound on the reduction in zippering energy given that we used the largest possible reduction in the excluded volume.
Description of partially toroidal states and calculation of membrane free energy
Request a detailed protocolDue to the finite size of the NLP, toroidal states are not possible for large pores. We instead assumed that these shapes are partially toroidal. These come into existence when , where is the membrane thickness. We set the shape of these states as follows. We constructed a toroidal pore with the shape parameters , and truncated the top half of the toroid at the plane where the cross-sectional diameter of the toroid is equal to the NLP diameter (Figure 7c, right panel). In these partially toroidal states, the ApoE proteins at the edge of the NLP are rotated through an angle compared with the fully toroidal states (Equation 6).
Evaluating the integral from Equation 2
where , and is the bending energy in a fully toroidal state, given by Equation 3. The change in total membrane area due to pore formation in these partially toroidal states is given by
In the partially toroidal states, does not correspond to the membrane separation at the edge of the pore, although is still the pore radius. Thus, we indicate the membrane separation at the NLP boundary by given by
References
-
Perspectives on kiss-and-run: role in exocytosis, endocytosis, and neurotransmissionAnnual Review of Physiology 75:393–422.https://doi.org/10.1146/annurev-physiol-020911-153305
-
Exocytotic fusion pores are composed of both lipids and proteinsNature Structural and Molecular Biology 23:67–73.https://doi.org/10.1038/nsmb.3141
-
Quantifying size distributions of nanolipoprotein particles with single-particle analysis and molecular dynamic simulationsThe Journal of Lipid Research 49:1420–1430.https://doi.org/10.1194/jlr.M700586-JLR200
-
Frequency spectrum of the flicker phenomenon in erythrocytesJournal De Physique 36:1035–1047.https://doi.org/10.1051/jphys:0197500360110103500
-
Lipid-anchored Synaptobrevin provides little or no support for exocytosis or Liposome fusionJournal of Biological Chemistry 291:2848–2857.https://doi.org/10.1074/jbc.M115.701169
-
A structural role for the synaptobrevin 2 transmembrane domain in dense-core vesicle fusion poresJournal of Neuroscience 35:5772–5780.https://doi.org/10.1523/JNEUROSCI.3983-14.2015
-
Mechanics of membrane fusionNature Structural and Molecular Biology 15:675–683.https://doi.org/10.1038/nsmb.1455
-
Membrane mechanics can account for fusion pore dilation in stagesBiophysical Journal 69:2489–2500.https://doi.org/10.1016/S0006-3495(95)80119-0
-
Dynamics of fusion pores connecting membranes of different tensionsBiophysical Journal 78:2241–2256.https://doi.org/10.1016/S0006-3495(00)76771-3
-
Exocytotic fusion pores exhibit semi-stable statesThe Journal of Membrane Biology 133:61–75.https://doi.org/10.1007/BF00231878
-
Complexin synchronizes primed vesicle exocytosis and regulates fusion pore dynamicsThe Journal of Cell Biology 204:1123–1140.https://doi.org/10.1083/jcb.201311085
-
Single vesicle millisecond fusion kinetics reveals number of SNARE complexes optimal for fast SNARE-mediated membrane fusionJournal of Biological Chemistry 285:32158–32166.https://doi.org/10.1074/jbc.A109.047381
-
Transmembrane-domain determinants for SNARE-mediated membrane fusionJournal of Cell Science 123:2473–2480.https://doi.org/10.1242/jcs.061325
-
Activity-dependent differential transmitter release in mouse adrenal chromaffin cellsJournal of Neuroscience 25:7324–7332.https://doi.org/10.1523/JNEUROSCI.2042-05.2005
-
Kiss-and-run exocytosis and fusion pores of secretory vesicles in human beta-cellsPflügers Archiv - European Journal of Physiology 457:1343–1350.https://doi.org/10.1007/s00424-008-0588-0
-
The debate on the kiss-and-run fusion at synapsesTrends in Neurosciences 30:447–455.https://doi.org/10.1016/j.tins.2007.06.012
-
Minimum membrane bending energies of fusion poresJournal of Membrane Biology 231:101–115.https://doi.org/10.1007/s00232-009-9209-x
-
SNARE complex zipping as a driving force in the dilation of proteinaceous fusion poresJournal of Membrane Biology 235:89–100.https://doi.org/10.1007/s00232-010-9258-1
-
SNAREs--engines for membrane fusionNature Reviews Molecular Cell Biology 7:631–643.https://doi.org/10.1038/nrm2002
-
Munc18-1 tuning of vesicle merger and fusion pore propertiesJournal of Neuroscience 31:9055–9066.https://doi.org/10.1523/JNEUROSCI.0185-11.2011
-
Calculating the bending Modulus for multicomponent lipid membranes in different thermodynamic phasesJournal of Chemical Theory and Computation 9:3866–3871.https://doi.org/10.1021/ct400492e
-
Stalk mechanism of vesicle fusion. intermixing of aqueous contentsEuropean Biophysics Journal : EBJ 17:121–129.https://doi.org/10.1007/BF00254765
-
Conformational dynamics of calcium-triggered activation of fusion by synaptotagminBiophysical Journal 105:2507–2516.https://doi.org/10.1016/j.bpj.2013.10.029
-
A conformational switch in complexin is required for synaptotagmin to trigger synaptic fusionNature Structural and Molecular Biology 18:934–940.https://doi.org/10.1038/nsmb.2103
-
The fusion poreBiochimica Et Biophysica Acta (BBA) - Molecular Cell Research 1641:167–173.https://doi.org/10.1016/S0167-4889(03)00085-5
-
High resolution electrophysiological techniques for the study of calcium-activated exocytosisBiochimica Et Biophysica Acta (BBA) - General Subjects 1820:1234–1242.https://doi.org/10.1016/j.bbagen.2011.12.011
-
Elastic curvature constants of lipid monolayers and bilayersChemistry and Physics of Lipids 144:146–159.https://doi.org/10.1016/j.chemphyslip.2006.08.004
-
Two modes of exocytosis in an artificial cellScientific Reports 4:3847.https://doi.org/10.1038/srep03847
-
SNARE complexes and neuroexocytosis: how many, how close?Trends in Biochemical Sciences 30:367–372.https://doi.org/10.1016/j.tibs.2005.05.002
-
Functional characterization of apolipoprotein E isoforms overexpressed in Escherichia coliProtein Expression and Purification 16:224–230.https://doi.org/10.1006/prep.1999.1069
-
The exocytotic fusion pore modeled as a lipidic poreBiophysical Journal 63:1118–1132.https://doi.org/10.1016/S0006-3495(92)81679-X
-
Hydration forces between phospholipid bilayersBiochimica Et Biophysica Acta (BBA) - Reviews on Biomembranes 988:351–376.https://doi.org/10.1016/0304-4157(89)90010-5
-
Chapter 11 - Reconstitution of membrane proteins in phospholipid bilayer nanodiscsMethods in Enzymology 464:211–231.https://doi.org/10.1016/S0076-6879(09)64011-8
-
The Synaptic Vesicle Release MachineryAnnual Review of Biophysics 44:339–367.https://doi.org/10.1146/annurev-biophys-060414-034057
-
Dopamine neurons release transmitter via a flickering fusion poreNature Neuroscience 7:341–346.https://doi.org/10.1038/nn1205
-
Cholesterol increases the openness of SNARE-Mediated flickering fusion poresBiophysical Journal 110:1538–1550.https://doi.org/10.1016/j.bpj.2016.02.019
-
One SNARE complex is sufficient for membrane fusionNature Structural and Molecular Biology 17:358–364.https://doi.org/10.1038/nsmb.1748
-
Synaptotagmin IV induces long-duration kiss-and-run exocytosis through small fusion poresBiophysical Journal 84:209a–225a.
-
Single-channel properties of IKs potassium channelsThe Journal of General Physiology 112:665–678.https://doi.org/10.1085/jgp.112.6.665
Article and author information
Author details
Funding
National Institute of General Medical Sciences (R01GM108954)
- Erdem Karatekin
Kavli Foundation (Neuroscience Scholar Award)
- Erdem Karatekin
Deutsche Forschungsgemeinschaft (VE760/1-1)
- Wensi Vennekate
National Institute of Diabetes and Digestive and Kidney Diseases (R01DK027044)
- Shyam S Krishnakumar
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank F Sigworth and D Zenisek for expert advice and discussions, and J E Rothman, M Caplan, and members of the Karatekin and Rothman laboratories for helpful discussions. We are grateful to James E Rothman, Frederick Sigworth, David Zenisek, Yongli Zhang, and Thomas Melia for critically reading the manuscript. Jing Wang kindly provided C45 lipids. This work was supported by NIH grant R01GM108954 and a Kavli Neuroscience Scholar Award to EK. SMA, OB, and SK are members of the Rothman laboratory, supported by NIH R01DK027044. WV was supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation), fellowship VE760/1-1.
Copyright
© 2017, Wu 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
-
- 2,124
- views
-
- 554
- downloads
-
- 58
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
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)
Further reading
-
- Structural Biology and Molecular Biophysics
In eukaryotes, RNAs transcribed by RNA Pol II are modified at the 5′ end with a 7-methylguanosine (m7G) cap, which is recognized by the nuclear cap binding complex (CBC). The CBC plays multiple important roles in mRNA metabolism, including transcription, splicing, polyadenylation, and export. It promotes mRNA export through direct interaction with a key mRNA export factor, ALYREF, which in turn links the TRanscription and EXport (TREX) complex to the 5′ end of mRNA. However, the molecular mechanism for CBC-mediated recruitment of the mRNA export machinery is not well understood. Here, we present the first structure of the CBC in complex with an mRNA export factor, ALYREF. The cryo-EM structure of CBC-ALYREF reveals that the RRM domain of ALYREF makes direct contact with both the NCBP1 and NCBP2 subunits of the CBC. Comparing CBC-ALYREF with other cellular complexes containing CBC and/or ALYREF components provides insights into the coordinated events during mRNA transcription, splicing, and export.
-
- Structural Biology and Molecular Biophysics
Experimental detection of residues critical for protein–protein interactions (PPI) is a time-consuming, costly, and labor-intensive process. Hence, high-throughput PPI-hot spot prediction methods have been developed, but they have been validated using relatively small datasets, which may compromise their predictive reliability. Here, we introduce PPI-hotspotID, a novel method for identifying PPI-hot spots using the free protein structure, and validated it on the largest collection of experimentally confirmed PPI-hot spots to date. We explored the possibility of detecting PPI-hot spots using (i) FTMap in the PPI mode, which identifies hot spots on protein–protein interfaces from the free protein structure, and (ii) the interface residues predicted by AlphaFold-Multimer. PPI-hotspotID yielded better performance than FTMap and SPOTONE, a webserver for predicting PPI-hot spots given the protein sequence. When combined with the AlphaFold-Multimer-predicted interface residues, PPI-hotspotID yielded better performance than either method alone. Furthermore, we experimentally verified several PPI-hotspotID-predicted PPI-hot spots of eukaryotic elongation factor 2. Notably, PPI-hotspotID can reveal PPI-hot spots not obvious from complex structures, including those in indirect contact with binding partners. PPI-hotspotID serves as a valuable tool for understanding PPI mechanisms and aiding drug design. It is available as a web server (https://ppihotspotid.limlab.dnsalias.org/) and open-source code (https://github.com/wrigjz/ppihotspotid/).