SNAP25 disease mutations change the energy landscape for synaptic exocytosis due to aberrant SNARE interactions

SNAP25 is one of three neuronal SNAREs driving synaptic vesicle exocytosis. We studied three mutations in SNAP25 that cause epileptic encephalopathy: V48F, and D166Y in the synaptotagmin-1 (Syt1)-binding interface, and I67N, which destabilizes the SNARE complex. All three mutations reduced Syt1-dependent vesicle docking to SNARE-carrying liposomes and Ca2+-stimulated membrane fusion in vitro and when expressed in mouse hippocampal neurons. The V48F and D166Y mutants (with potency D166Y > V48F) led to reduced readily releasable pool (RRP) size, due to increased spontaneous (miniature Excitatory Postsynaptic Current, mEPSC) release and decreased priming rates. These mutations lowered the energy barrier for fusion and increased the release probability, which are gain-of-function features not found in Syt1 knockout (KO) neurons; normalized mEPSC release rates were higher (potency D166Y > V48F) than in the Syt1 KO. These mutations (potency D166Y > V48F) increased spontaneous association to partner SNAREs, resulting in unregulated membrane fusion. In contrast, the I67N mutant decreased mEPSC frequency and evoked EPSC amplitudes due to an increase in the height of the energy barrier for fusion, whereas the RRP size was unaffected. This could be partly compensated by positive charges lowering the energy barrier. Overall, pathogenic mutations in SNAP25 cause complex changes in the energy landscape for priming and fusion.


Introduction
The fusion machinery responsible for chemical synaptic transmission is well known: it consists of the SNARE complex, a ternary complex formed by the proteins VAMP2, syntaxin-1, and SNAP25 (Sutton et al., 1998).This complex is under tight control by upstream partner protein Munc18-1, which acts as a template for SNARE complex formation, and Munc13s, which assist in the transitions required along the pathway of assembly (Rizo, 2022).SNARE complexes assemble in a zipper-like manner; partially assembled SNAREpins bind synaptotagmin-1 (Syt1) and complexin within two separate interaction sites: the primary interface formed between a synaptotagmin-1 molecule and syntaxin-1 and SNAP25 (Zhou et al., 2015), and the tripartite interface formed by syntaxin-1 and VAMP2 with another Syt1 molecule and complexin (Zhou et al., 2017).Upon arrival of an action potential, Ca 2+ binds to the two C2 domains of Syt1, which results in rapid vesicle-plasma membrane fusion and release of neurotransmitter within a fraction of a millisecond (Südhof, 2013).
The strong functional integration and specialization of the neuronal SNARE for speed have rendered the release machinery exquisitely susceptible to insults.De novo mutations in SNAREs and associated proteins lead to complex neurological disease, characterized by drug-resistant epilepsy, intellectual disability, movement disorders, and often autism; a syndrome which has been denoted 'SNAREopathy' (Verhage and Sørensen, 2020).Although rare, these are devastating conditions for the patients and their families.Consequently, there is considerable interest in revealing the molecular/cellular mechanisms for these conditions, which is seen as key to the development of treatment.Diseasecausing mutations in the SNARE-machinery fall into distinct categories according to the nature of the defect (Verhage and Sørensen, 2020).These include different forms of haploinsufficiency, where the mutated protein is either lost altogether or has lost its functionality, and dominant-negative or recessive variants, as well as variants with new or changed protein interactions, referred to as neomorphs.
Most SNAREopathy mutations have been described in STXBP1 (Abramov et al., 2021a;Verhage and Sørensen, 2020;Xian et al., 2022), the gene encoding Munc18-1.These mutations are generally found to cause Munc18-1 hypo-expression, probably due to protein instability (Guiberson et al., 2018;Kovacevic et al., 2018;Martin et al., 2014;Saitsu et al., 2008), and thus the mutations belong in the haploinsufficiency category.Synaptic phenotypes of Munc18-1 hypo-expression have been identified in the Drosophila neuromuscular junction (Wu et al., 1998), in cultured mouse neurons eLife digest Neurons in the brain communicate with one another by passing molecules called neurotransmitters across the synapse connecting them together.Mutations in the machinery that controls neurotransmitter release can lead to epilepsy or developmental delays in early childhood, but how exactly is poorly understood.
Neurotransmitter release is primarily controlled by three proteins that join together to form the SNARE complex, and another protein called synaptotagmin-1.This assembly of proteins primes vesicles containing neurotransmitter molecules to be released from the neuron.When calcium ions bind to synaptotagmin-1, this triggers vesicles in this readily releasable pool to then fuse with the cell membrane and secrete their contents into the small gap between the communicating neurons.
Mutations associated with epilepsy and developmental delays have been found in all components of this release machinery.Here, Kádková, Murach, Østergaard et al. set out to find how three of these mutations, which are found in a protein in the SNARE complex called SNAP25, lead to aberrant neurotransmitter release.
Two of these mutations are located in the interface between the SNARE complex and synaptotagmin-1, while the other is found within the bundle of proteins that make up the SNARE complex.In vitro and ex vivo experiments in mice revealed that the two interface mutations led to defects in vesicle priming, while at the same time bypassing the control by synaptotagmin-1, resulting in vesicles spontaneously fusing with the cell membrane in an unregulated manner.These mutations therefore combine loss-of-function and gain-of-function features.
In contrast, the bundle mutation did not impact the number of vesicles in the releasable pool but reduced spontaneous and calcium ion evoked vesicle fusion.This was due to the mutation destabilizing the SNARE complex, which reduced the amount of energy available for merging vesicles to the membrane.
These findings reveal how SNAP25 mutations can have different effects on synapse activity, and how these defects disrupt the release of neurotransmitters.This experimental framework could be used to study how other synaptic mutations lead to diseases such as epilepsy.Applying this approach to human neurons and live model organisms may lead to the discovery of new therapeutic targets for epilepsy and delayed development.(Toonen et al., 2006) and in brain slices of STXBP1 heterozygous mice (Chen et al., 2020;dos Santos et al., 2023), as well as in human neurons (Patzke et al., 2015), whereas expression of human missense mutations in STXBP1 null mouse neurons led to synaptic phenotypes for some mutations, but not for others (Kovacevic et al., 2018).Attempts to rescue the disease phenotypes have focused on mechanisms for increasing expression levels, for instance chemical chaperones that might prevent Munc18-1 misfolding and degradation (Abramov et al., 2021b;Guiberson et al., 2018).
The situation is different for disease mutations in Syt1 and SNAP25.In these cases, mutation generally does not cause protein instability, but instead changes function in different ways.For Syt1, disease mutations were found to cluster in the C2B domain around the top loops that coordinate Ca 2+ (Baker et al., 2018).Expression of these mutants in neurons demonstrated reduced evoked release, in some cases an increase in spontaneous release, and a dominant-negative phenotype when co-expressed with wildtype (WT) protein (Bradberry et al., 2020).The molecular mechanism was identified as a decrease in Ca 2+ -dependent lipid binding (Bradberry et al., 2020).In SNAP25, disease-causing mutations are found within the SNARE domains (Hamdan et al., 2017;Klöckner et al., 2021;Rohena et al., 2013;Shen et al., 2014).Alten et al., 2021 studied a selection of SNAP25 mutants, and found no changes in expression levels, but changes in both spontaneous and evoked release.Specifically, mutations in the primary Syt1:SNARE interface (Zhou et al., 2015) caused an increase in spontaneous release rates, and a decrease in evoked release amplitudes.Conversely, C-terminal mutations in the so-called 'layer residues', whose side chains point to the center of the SNARE complex and are involved in SNAREpin zippering (Sutton et al., 1998), led to a decrease in both evoked and spontaneous release (Alten et al., 2021).This is expected because the C-terminal end of the SNARE complex is required for both types of release (Weber et al., 2010).Alten et al., 2021 described striking phenotypes for most SNAP25 mutations tested, but the molecular reason for these phenotypes remains incompletely understood, and a few findings were surprising.For instance, Alten et al. reported that the V48F and D166Y mutants supported an unchanged readily releasable pool (RRP) of vesicles.The RRP is the pool of vesicles available for immediate release upon invasion of the presynapse by the action potential, and it is often assessed by applying a hyperosmotic solution (often 0.5 M sucrose), which causes primed vesicles to fuse (Rosenmund and Stevens, 1996).Previously, we concluded that the primary Syt1:SNARE interface is involved in vesicle priming (Schupp et al., 2016), which agrees with the suggestion that Syt1 binds to the SNAREs before Ca 2+ arrival (Zhou et al., 2015).Conversely, Alten et al. reported a smaller RRP for the I67N mutant, but the I67 residue is present in the internal of layer +4, which we expected to affect final SNARE complex zippering causing membrane fusion rather than priming (Sørensen et al., 2006;Weber et al., 2010).These discrepancies might be explained by the fact that Alten et al. used mixed cultures and mainly studied GABAergic transmission.Glutamatergic transmission is hard to investigate in mixed cultures due to reverberating activity, but it can be studied using autaptic neurons, which was our approach (Schupp et al., 2016;Weber et al., 2010).Thus, there is a need to further study these disease mutations in glutamatergic neurons.Another open question is whether the phenotype of the mutants in the primary interface (V48F and D166Y) is explainable solely by the loss of Syt1 binding, or whether other features of these mutations add to, or detract from, the phenotype.
Here, we reexamined three SNAP25 disease mutations using glutamatergic autaptic neurons: I67N, V48F, and D166Y.We confirmed the increase in spontaneous release rate and decrease in evoked release by V48F and D166Y reported previously (Alten et al., 2021).Through in vitro analysis we find that both mutations lower Syt1 association to SNARE-protein liposomes.Additional experiments both in vitro and in cells demonstrate that V48F and especially D166Y represent partial gainof-function mutations that increase association to partner SNAREs and lower the energy barrier for fusion, bypassing Syt1-dependent control.Thus, these mutants do not phenocopy the loss of Syt1, but combine loss of Syt1 binding with a gain-of-function phenotype.At the same time, the mutants act as loss-of-function in upstream reactions, through effects on priming.The I67N is a classical dominant-negative mutation, which increases the energy barrier for fusion, but does not change the size of the RRP if probed by a sufficiently high concentration of sucrose, nor does the I67N change the electrostatics of triggering itself.Thus, for V48F and D166Y, loss-of-function and gain-of-function features combine to change the energy landscape for vesicle priming and fusion in a complex way.These findings have consequences for our understanding of the simultaneous role of the primary SNARE:Syt1 interface in vesicle priming and release clamping.It further demonstrates the challenge faced by finding mechanism-based treatments of these disorders in the presence of multiple effects caused by single-point mutations.

Results
We first investigated two SNAP25 disease-causing mutations within the primary Syt1:SNARE interface (V48F and D166Y; Figure 1A-C).The mutations occurred de novo and were identified in single heterozygous subjects.The V48F mutant was identified in a 15-year-old female with encephalopathy, intellectual disability and generalized epilepsy with seizures started at 5 months of age; Magnetic Resonance Imaging (MRI) was normal except for delayed myelination (Rohena et al., 2013).The D166Y mutants were found in a 23-year-old male with global developmental delay, nocturnal tonicclonic seizures, and moderate intellectual disability; the MRI showed mild diffuse cortical atrophy (Hamdan et al., 2017).We aimed to achieve a detailed understanding of the reason for synaptic dysfunction.

The V48F and D166Y mutants disinhibit spontaneous release and desynchronize evoked release
We constructed lentiviral vectors, which expressed WT or mutant SNAP25b fused N-terminally to EGFP as an expression marker (Delgado-Martínez et al., 2007).In the absence of SNAP25, neuronal viability is compromised (Delgado-Martínez et al., 2007;Peng et al., 2013;Santos et al., 2017;Weber et al., 2010) with low-density autaptic glutamatergic neurons dying within a few days, whereas Snap25 KO neurons growing at higher densities, or in the intact brain, can survive longer, but also eventually degenerate (Bronk et al., 2007;Hoerder-Suabedissen et al., 2019).We therefore first examined the morphology and survival of autaptic hippocampal neurons from Snap25 KO after expressing mutated or WT EGFP-SNAP25b (henceforth denoted 'WT' in rescue experiments; not to be confused with 'Syt1 WT', which refers to wildtype littermates of Syt1 KO mice).
As expected, Snap25 KO autaptic neurons did not survive in the absence of expression of exogeneous SNAP25 (Figure 2G), whereas expression of WT EGFP-SNAP25b restored survival (Delgado-Martínez et al., 2007;Ruiter et al., 2019;Weber et al., 2010).Both mutations (V48F and D166Y) caused rescue of survival; however, the number of neurons per islet was lower than in WT expressing neurons prepared in parallel.This difference was statistically significant for D166Y, but not for V48F (Figure 2G).Since patients harbor one mutated and one WT allele, we co-expressed WT EGFP-SNAP25b with mutated EGFP-SNAP25b, by combining infection with separate viruses encoding WT and mutant protein in a 1:1 ratio, as done previously by others for Syt1 (Bradberry et al., 2020).As a prerequisite for this, we demonstrated by Western blot analysis that both mutant viruses expressed similar amounts of protein as WT viruses (Figure 2A-C).In co-expressing neurons, we added half the volume of WT and mutant virus as compared to the WT condition, thus keeping the total amount of virus constant.We preferred to combine two single-cistronic viruses, rather than to construct bicistronic viruses, since the larger insert would be expected to lower the titer of viruses, which  might compromise survival rescue of Snap25 KO neurons.Using co-expression of WT and mutant SNAP25b, neuronal survival was mildly and significantly reduced in the V48F + WT condition from WT (Figure 2G).Overall, mutations seemed to mildly reduce survivability of neurons.Staining of neurons against VGlut1 (a marker for glutamatergic synapses) and MAP2 (a dendritic marker) revealed no significant difference in the dendritic length or the number of synapses in neurons expressing V48F, and D166Y alone (Figure 2D-F).
Patch-clamp was performed on days 10-14 on autaptic neurons expressing either the WT EGFP-SNAP25b, the mutants, or both the mutant and WT in a 1:1 ratio.Spontaneous miniature events were recorded at a holding voltage of −70 mV.Both primary interface mutations (V48F and D166Y) strongly increased the frequency of mEPSCs (Figure 3A, B, D, E); the mEPSC amplitude was significantly increased in the V48F, and insignificantly (p = 0.14) increased for the D166Y (Figure 3C, F).The mEPSC frequency was at least as high as in Syt1 KO neurons measured in separate experiments (Figure 3G-I).V48F and D166Y co-expressed with WT resulted in mEPSC frequencies close to the arithmetic mean between frequencies in WT and the mutant alone (Figure 3B, E), indicating that the mutations are incompletely dominant, or co-dominant.
Brief depolarization elicits an unclamped action potential in the axon, which makes it possible to study evoked release, which is essentially absent in Snap25 KO neurons (Bronk et al., 2007;Delgado-Martínez et al., 2007).The V48F and D166Y mutants both supported evoked release, but the evoked EPSC (eEPSC) amplitude was significantly reduced in the mutant condition (Figure 4A, B, E, F), whereas when co-expressed with the WT, only the D166Y mutant had a significantly reduced eEPSC (Figure 4B,F).Integration of individual eEPSCs allowed determination of the total charge and the assessment of synchronous and asynchronous release components (Figure 4-figure supplement 1).The eEPSC charge was significantly reduced in the D166Y, but not in the V48F (Figure 4C, G), whereas in both cases the kinetics was significantly shifted in the direction of more asynchronous release, and the fast time constant was prolonged (Figure 4D, H; Figure 4-figure supplement 2).We compared the data obtained from V48F and D166Y SNAP25 to Syt1 KO neurons recorded in separate experiments.Syt1 WT neurons in this set of experiments had larger eEPSC amplitudes, and a larger charge than WT-rescued Snap25 KO neurons (Figure 4I-K), which might be caused by differences between cell cultures, or animal lines (SNAP25 vs. Syt1 lines).Note that our experiments using WT, mutant, and WT + mutant SNAP25 were always carried out in neurons prepared and recorded in parallel from the same Snap25 KO embryos (Materials and methods).In the Syt1 KO, the kinetic change was similar to the V48F and D166Y (Figure 4L), but the total charge was also strongly reduced (Figure 4K), unlike the two mutations that displayed at most a mild reduction (Figure 4C,G).
Overall, D166Y and V48F caused a strong disinhibition of spontaneous release and a desynchronization of evoked release, as demonstrated before (Alten et al., 2021), and consistent with previous mutational studies of the primary interface, both in SNAP25 (Schupp et al., 2016) and Syt1 (Zhou et al., 2015).Accordingly, these phenotypes are similar to Syt1 KO, but the effect of V48F and D166Y on total evoked charge was milder than in the Syt1 KO.The preserved overall charge of evoked release in V48F and mild reduction in D166Y might point to a compensatory gain-of-function aspect to these mutations, in addition to the impaired Syt1 interaction.
V48F and D166Y mutations lower the energy barrier for release and reduce the RRP size The role of Syt1 and Syt1:SNARE interactions in vesicle priming has been controversially discussed.The RRP is often assessed by applying a pulse of hypertonic solution to the neurons, usually 0.5 M sucrose (Rosenmund and Stevens, 1996;Schotten et al., 2015).However, in some experiments this did not lead to a change in RRP size in the absence of Syt1 (Bacaj et al., 2015;Xu et al., 2009), whereas experiments both from our laboratories and others showed a decrease in RRP in the absence of Syt1 (Bouazza-Arostegui et al., 2022;Chang et al., 2018;Courtney et al., 2019;Huson et al., Source data 2. Original files for the Western blot analysis in Figures .
Source data 3. PDF containing Figure 2A and 10A, and original scans of the Western blots with highlighted bands and sample labels.Liu et al., 2009;Ruiter et al., 2019).The reason for this discrepancy is likely partly technical and has to do with how fast sucrose can be applied to the dendritic tree (see also Discussion).Using neurons growing on small autaptic islands makes it possible to apply sucrose to the entire dendritic tree (which is confined to the 30-50 µm island) within tens of milliseconds, which is fast enough to    dissect the RRP with a short (few s) sucrose application.We set out to understand whether the V48F and D166Y mutants changed the size of the RRP.Application of 0.5 M sucrose to neurons with the V48F or D166Y mutation resulted in estimates of the RRP (denoted RRP 0.5 ) that were significantly reduced compared to the WT condition (Figure 5A, The online version of this article includes the following source data for figure 5: Source data 1.Quantitative data. C, F, H).Also the co-expressed condition displayed significantly reduced RRP 0.5 .Application of two different sucrose concentrations is used to probe the size of the energy barrier for release (Basu et al., 2007;Schotten et al., 2015), because the use of a lower sucrose concentration (typically 0.25 M) will only release a fraction of the RRP, and this fraction depends sensitively on the energy barrier for release; the size of the RRP at 0.5 M sucrose (RRP 0.5 ) is used for normalization.Application of 0.25 M sucrose to V48F and D166Y expressing neurons strikingly led to unchanged pool release (RRP 0.25 ) (Figure 5B, G).Consequently, the ratio of pools (RRP 0.25 /RRP 0.5 ) was significantly increased for both the V48F and the D166Y mutations (Figure 5D, I), indicating that these two mutations decrease the apparent energy barrier for fusion (Schotten et al., 2015).The co-expressed conditions were in-between WT and mutant and did not reach statistical significance.Consistent with these results, both mutations on average increased the release probability, calculated as the ratio between the eEPSC charge and the RRP 0.5 charge measured in the same cell.The increase was statistically significant for the D166Y mutation (Figure 5J), but not for V48F (Figure 5E).This is different from the situation in the Syt1 KO, where the RRP 0.25 /RRP 0.5 and the release probability were both significantly decreased (Figure 5K-O), although in this dataset the reduction in RRP 0.5 size by removing Syt1 did not reach statistical significance (p = 0.0548, unpaired t-test).Thus, the increased mEPSC release rate from the Syt1 KO does not correlate with a reduced energy barrier as assayed by sucrose, which was reported before (Huson et al., 2020), and the V48F and D166Y have specific gain-of-function features that lower the energy barrier for vesicle fusion.
To investigate the reasons for the change in RRP size by D166Y and V48F, we considered a onepool model for the RRP (Figure 6A), where the RRP is filled by priming (k 1 ) from an upstream pool and depleted by depriming (k −1 ) or spontaneous fusion (k f ).Dividing the miniature release rates (r mini ) with the size of the RRP yields the spontaneous fusion rate k f .The current plateau during 0.5 M sucrose application (Figure 6B) essentially reports on the priming rate (k 1 ), providing that the fusion rate is sufficiently increased by sucrose (Materials and methods).However, high sucrose concentrations can change the baseline current level, which will cause an error in estimation of k 1 .We therefore corrected the plateau level using a plot of the variance versus mean during the sucrose application (Figure 6C); this plot is linear for the type of noise generated by synaptic transmission (shot noise).Back-extrapolation of a regression line allows a determination of the baseline current level in the presence of sucrose (Figure 6B, C; see Materials and methods).Combining RRP size with the estimates of k 1 and k f allows determining the depriming rate, k −1 .
These calculations showed that k f , the spontaneous fusion rate was strongly increased in the V48F and D166Y, and this increase was much larger for the mutations than for the Syt1 KO (Figure 6D) fulfilling D166Y > V48F > Syt1 KO.Further analysis showed that in both mutations, the forward priming rate, k 1 , and the depriming rate, k −1 , were both decreased -the latter effect was only significant for D166Y (Table 1).Summing up the effects of changes in the three parameters on the RRP size for the V48F, D166Y, and Syt1 KO (Figure 6E), we can conclude that for the V48F and even more for the D166Y, spontaneous release contributes to the reduction in RRP size, whereas this effect is minimal in the Syt1 KO.However, the major reason for the smaller RRP size is a reduction in priming rate, k 1 , which is partly counteracted by the decrease in k −1 (which would increase RRP size).Overall, changes in priming, depriming, and spontaneous fusion rates combine to change RRP size.
Repetitive stimulation to determine the RRPs often results in lower estimates than sucrose application, because action potentials draw on a sub-pool of the RRP (Moulder and Mennerick, 2005), whereas sucrose releases the entire RRP (Rosenmund and Stevens, 1996).To determine the RRP sub-pool that evoked release draws on (denoted RRP ev ) we applied repetitive stimulation (50 APs @ 40 Hz) and used back-extrapolation to determine the RRP (Neher, 2015).Performing this in our standard 2 mM Ca 2+ -containing extracellular solution resulted in overall smaller estimates for RRP ev for V48F and D166Y compared to WT (Figure 7-figure supplement 1); however, the differences were not statistically significant.The back-extrapolation method works best with high release probabilities (Neher, 2015); therefore, we repeated these experiments in the presence of 4 mM extracellular Ca 2+ .Under these conditions, the RRP ev was significantly reduced for both the V48F and D166Y mutation (Figure 7C, G; Figure 7-figure supplement 2).Strikingly, the release probabilities (calculated as the charge of the first eEPSC of the train divided by the RRP ev ) were decreased for both mutations (Figure 7D, H), which correlated with an increased paired-pulse ratio (Figure 7A, E, inserts) (Zucker and Regehr, 2002).Thus, although the release probability calculated by normalizing evoked charge to

Corrected baseline Current plateau
Var/mean @ baseline Var/mean @ sucrose  the sucrose-determined RRP (RRP 0.5 ) was increased (see above), the release probability when normalizing to RRP ev was decreased.This points to a difference in the organization of the RRP ev and RRP 0.5 (see Discussion).The forward priming rate is determined as the slope of the linear fit used for the back-extrapolation; this parameter was reduced in both mutations (Figure 7B, F; Figure 7-figure supplement 2).Thus, the difference in priming rate extends to the RRP ev .
Overall, these data demonstrate a rather complex phenotype of the D166Y and V48F mutations, which combine a lowering of the energy barrier -a gain-of-function feature -with a loss of vesicle priming -a loss-of-function feature.The D166Y and V48F phenotypes can be summarized as (1) desynchronized eEPSCs with at most mildly reduced total charge, (2) lowered energy barrier for fusion, (3) increased release probability when normalized to the sucrose pools, (4) decreased RRP size due to unclamped spontaneous release and lowered forward priming rates, and (5) short-term facilitation.These phenotypes are distinctive from the Syt1 KO, which does not have a preserved charge of the eEPSC, or lowered energy barrier when probed by sucrose, or an increased release probability.Thus, the D166Y and V48F cannot be understood solely in terms of a lack of Syt1 coupling; instead, gain-offunction features are present in the mutants which are absent upon deletion of Syt1.

V48F and D166Y mutants show increased partner SNARE interactions and cause unregulated fusion in vitro
We next tried to identify the biochemical properties of V48F and D166Y, which could support a gainof-function phenotype during exocytosis.These data are displayed in Figures 8 and 9 -data on the I67N mutant were obtained in parallel and will be presented later.To test to which degree the mutants may change the stability of SNARE complexes, full-length t-SNAREs (syntaxin-1 and SNAP25) were incubated with the cytosolic domain of VAMP2 (VAMP2cd) overnight.Cis-SNARE complex stability was tested in the presence of SDS at the indicated temperatures (Figure 8A), and the release of syntaxin-1 as a single-protein band was used as a measure of the complex dissociation.The WT, V48F, and D166Y v-/t-SNARE complexes showed a similar stability with half-maximal dissociation occurring at approximately 71°C.
To understand how the mutations affect fusion in a well-defined reconstituted membrane fusion system, we performed in vitro lipid mixing assays using GUVs containing both t-SNAREs (syntaxin-1 and SNAP25) with 1% PI(4,5)P 2 , and 0.5% Atto488/0.5% Atto550-labeled SUVs containing Syt1 and VAMP2 in the presence or absence of complexin-II (6 µM) (Kedar et al., 2015;Malsam et al., 2012).Fusion was measured at 37°C by Atto488 fluorescence dequenching, which occurs upon lipid mixing with GUVs.Calcium (100 µM free Ca 2+ final) was added after 2 min to the t-SNARE GUV assay.Measurements were continued for another 5 min.SUVs treated with botulinum toxin D, which cleaves VAMP2 and abolishes membrane fusion, served as negative control and the corresponding    Complexin-II clamps spontaneous Ca 2+ -independent membrane fusion in the reconstituted assay (i.e.fusion before addition of Ca 2+ ), via laterally binding the membrane-proximal C-terminal ends of SNAP25 and VAMP2 (Malsam et al., 2020).Remarkably, V48F and D166Y showed impaired clamping by complexin, as apparent by increased fusion before Ca 2+ addition (Figure 8D).The decreased clamping is likely caused by the reduced interaction of V48F and D166Y with Syt1.As a note, the clamping function of Syt1 becomes only obvious in the presence of complexin, whereas the clamping function of complexin depends on the presence of Syt1 (Malsam et al., 2012); thus, the clamping function of the two proteins cannot be separated.After Ca 2+ triggering, WT SNAP25 supported the largest amount of fusion, followed by V48F, whereas D166Y Ca 2+ -dependent fusion was clearly reduced (Figure 8D).Notably, this is the same sequence as found for evoked release in synapses, considering either the amplitude or the charge of the eEPSC (Figure 4B, C, F, G).Thus, the in vitro assay reproduces both the increased spontaneous release and the reduced Ca 2+ -dependent release found in neurons, indicating that these features are present within the minimal set of fusion proteins included in this assay (i.e. the SNAREs, Syt1 and complexin).
There is evidence that SNAP25 might enter the SNARE complex last, after syntaxin-1 and VAMP2 are joined by Munc18-1 (Baker et al., 2015;Jiao et al., 2018;Sitarska et al., 2017), although another view is that a syntaxin/SNAP25 dimer bound to Munc18-1 acts as an intermediary (Jakhanwal et al., 2017).In the former case, mutations changing association of SNAP25 to the SNAREs might change exocytosis efficiency.To test this, we used a syntaxin-1 GUV assay, where the incorporation of soluble SNAP25 into the SNARE complex becomes a rate-limiting step.In this assay, fusion kinetics are much slower when compared with fusion reactions containing GUVs with preassembled t-SNAREs.Accordingly, we allowed 30 min for pre-stimulation fusion to take place, and after addition of Ca 2+ (100 µM free Ca 2+ ), fusion was followed for another 30 min (Figure 9A).The assay was performed in the presence and absence of complexin-II.D166Y and to a lesser degree V48F revealed enhanced stimulation of fusion in comparison to WT before calcium was added, regardless of the presence or absence of complexin (Figure 9B).Comparing data with and without complexin established that complexin barely suppressed spontaneous fusion in reactions containing the V48F mutant.The presence of complexin partially reduced Ca 2+ -independent fusion, but D166Y still stimulated pre-Ca 2+ fusion compared to WT.Overall, these data demonstrate that D166Y and V48F are gain-of-function mutants under conditions where SNAP25 association to the other SNAREs is rate limiting.
To test directly whether V48F and especially D166Y enhance SNARE interactions, co-flotation assays were performed.SUVs containing either syntaxin-1 (Stx-1), or VAMP2, or Syt1, or the Syt1/ VAMP2 combination were incubated with SNAP25 and re-isolated by flotation using a Nycodenz density gradient.An additional reaction, reflecting the fusion assay, contained SNAP25 in combination with both Syt1/VAMP2 SUVs and Syntaxin-1 SUVs (Figure 9D).SNAP25 recruitment for each condition was determined by SDS-PAGE followed by Coomassie Blue and silver staining.Silver stained SNAP25 bands were quantified, and the mutants were plotted relatively to the WT SNAP25 (Figure 9D and Figure 9-figure supplement 1).
SNAP25 WT did not show any binding to Syt1 SUVs (Figure 9-figure supplement 1), although direct interactions with Syt1 would be predicted based on the primary interface, indicating that such interactions are not stable under the employed conditions, which is expected because SNAP25 is unstructured until it binds its SNARE partners (Fasshauer et al., 1997).D166Y showed profoundly increased interactions with SUVs containing either syntaxin-1, or VAMP2, or Syt1/VAMP2, and the combination used in the fusion assay (Figure 9D; Figure 9-figure supplement 1).V48F displayed mildly increased binding to syntaxin-1 and Syt1/VAMP1 SUVs, and a tendency to increased association to VAMP2 SUVs (p = 0.0629, one-sample t-test).These data show that loss of Syt1 interaction upon mutation in the primary interface can be accompanied by a gain-of-function phenotype stimulating  interactions with the other SNARE partners.This association of SNAP25 to the other SNAREs might happen as one of the last steps toward fusion; consequently, when D166Y and V48F join the complex prematurely, it will bypass layers of control and result in uncontrolled fusion.
The I67N mutation supports an intact RRP, but an increased energy barrier for fusion We next addressed the phenotype of the I67N disease mutation, which was found in an 11-year-old female, who suffered from myasthenia, cortical hyperexcitability, ataxia, and intellectual disability, but with normal brain MRI (Shen et al., 2014).The I67 is found within the interaction layer +4 (Fasshauer et al., 1998), which helps in assembly of the C-terminal of the SNARE complex (Gao et al., 2012), and it might therefore have a different synaptic phenotype than V48F and D166Y.
In in vitro experiments, SNARE complexes formed with the I67N mutant displayed a lower stability, with the melting temperature reduced from 71°C to approximately 56°C (Figure 8A), as expected for a mutation that destabilizes the SNARE complex.I67N also showed a strong decrease in SUV docking (Figure 8B, C), which is likely caused by the destabilization of the t-SNARE complex, which indirectly perturbs the Syt1-binding interface(s).In lipid mixing assays, I67N strongly reduced both Ca 2+ -independent and Ca 2+ -dependent fusion, whether the t-SNAREs were preassembled (Figure 8D,  E) or not (Figure 9B, C).The co-flotation assay did not display any binding of I67N to t-SNAREs or Syt1 (Figure 9D).The binding to Syt1/VAMP2 SUVs was reduced compared to SNAP25 WT, but since binding of the WT protein is already very low, the biological significance of this result is unclear.Overall, these data indicate that I67N is inferior in membrane fusion.
Lentiviruses encoding the I67N mutant N-terminally fused to EGFP expressed similar amounts as WT EGFP-SNAP25b (Figure 10A-C).Expression in SNAP25 KO neurons resulted in reduced rescue of survival in neurons expressing I67N alone, whereas neurons co-expressing WT and I67N had intermediate survival, not significantly different from WT (Figure 10D).Staining against MAP2 (dendritic marker) and VGlut1 (synapse marker) showed that the number of synapses on average was reduced in the I67N (Figure 10E, F), and the dendritic length was on average reduced (Figure 10G), but the changes did not reach statistical significance.Patch-clamp measurements demonstrated strongly reduced spontaneous release frequencies (Figure 10H-J) and evoked release amplitude (Figure 10K, L) with the I67N mutation, see also Alten et al., 2021.mEPSCs were absent in most I67N expressing neurons, whereas WT and I67N co-expressing neurons had a very low mEPSC rate, much closer to the I67N than the WT phenotype; similar for eEPSC amplitudes and charges (Figure 10L-M).The I67N therefore is dominant negative for both types of release, in contrast to the incompletely dominant phenotypes of the V48F and D166Y mutants (see above).The fraction of synchronous release was unchanged in WT and I67N co-expressing neuron (Figure 10N); this number could not be estimated for I67N expressing neurons due to the low amount of release.
Reduced spontaneous and evoked release could result from a decrease in priming, or fusion, or both.To distinguish between these possibilities, we turned to sucrose applications.Application of 0.25 M sucrose did not lead to any measurable release in I67N expressing cells, and only minimal release in cells co-expressing WT and I67N (Figure 11B).This indicates that the energy barrier is 0.0001; ***p < 0.001; **p < 0.01; *p < 0.05, analysis of variance (ANOVA) with Dunnett's multiple comparisons test, comparing each mutation to the corresponding wildtype (WT) condition.(D) SNAP25b D166Y showed enhanced interactions with SUVs carrying reconstituted syntaxin-1A (Stx-1), VAMP2, Syt1/VAMP2, or an SUV mixture containing Syntaxin-1A and VAMP2/Syt1 in co-flotation assays, whereas V48F displayed weaker increases in interactions with SUVs containing syntaxin-1A, or Syt1/VAMP2.Shown is mean ± SEM on a logarithmic scale.***p < 0.001, **p < 0.01, *p < 0.05, twotailed one-sample t-test comparing to 1.
The online version of this article includes the following source data and figure supplement(s) for figure 9: Source data 1.Excel file containing quantitative data.increased in amplitude, and therefore the RRP might be underestimated when probed by 0.5 M sucrose (Schotten et al., 2015).Indeed, 0.5 M sucrose displayed reduced release in the I67N and WT + I67N condition (Figure 11C), but this could be due to defects in priming or fusion.To investigate this, we applied a stronger stimulus, 0.75 M sucrose, to these cells (Schotten et al., 2015).Strikingly, the RRP as assessed by 0.75 M sucrose (RRP 0.75 ) was unchanged between WT, I67N, and WT + I67N co-expressed cells (Figure 11E, G).Application of 0.375 M sucrose led to small amounts of release in the I67N, but more in the WT + I67N co-expressed situation (Figure 11F).Forming the ratio RRP 0.375 / RRP 0.75 revealed a statistically significant reduction in I67N compared to WT (Figure 11H).Therefore, the RRP per se appears intact in the I67N (when probed by sufficiently high concentrations of sucrose), but the vesicles face a higher apparent fusion barrier, which makes the RRP appear smaller if assessed by 0.5 M sucrose.The higher apparent fusion barrier explains the lower frequency of mEPSC in the I67N mutation, the lower degree of spontaneous fusion in in vitro assays, as well as the lower amount of Ca 2+ -dependent release in vitro and in the cell.The I67N mutation profoundly affected trains in 2 mM extracellular Ca 2+ (Figure 11I) leading to strong facilitation, which is expected due to the strong phenotype of this mutation, which radically lowers release probability.Even when co-expressed with WT protein, the train facilitated over the first several stimulations, attesting to the strong dominant-negative feature of the I67N mutation (Figure 11I-J).Consequently, the paired-pulse ratio was increased in the I67N and intermediate in the WT + I67N condition (Figure 11K).Back-extrapolation of these trains was not reliable, because the low release probability in the I67N mutation made it impossible to achieve sufficient depletion of the RRP.
The energy barrier for fusion is exquisitely sensitive to the charges on the surface of the SNARE complex, with positive charges decreasing and negative charges increasing the fusion barrier amplitude (Ruiter et al., 2019).To investigate whether the same electrostatic mechanism applies to the I67N, we combined the I67N mutation with a mutation of four amino acids ('4K' = SNAP25 E183K/ S187K/T190K/E194K) in the second SNARE domain of SNAP25, constructing the quintuple mutation ('I67N/4K' = SNAP25 I67N/E183K/S187K/T190K/E194K).The 4K mutation lowers the energy barrier for fusion by increasing the charge of the SNARE complex surface by +6 via charge introduction and charge reversal (Ruiter et al., 2019).The 4K mutation increased the mEPSC release rate compared to WT (Ruiter et al., 2019), whereas in the combined I67N/4K mutation the spontaneous release rate was indistinguishable from WT (Figure 12A-C), showing that increased positive charges rescued the defect of spontaneous release in the I67N mutant.Evoked release was also increased in the 4K mutation compared to WT (Ruiter et al., 2019), but in the I67N/4K mutation, evoked release was still strongly depressed compared to WT (Figure 12D, E).Nevertheless, evoked release was noticeable in the combined mutation, whereas it was almost absent in I67N expressing cells (comp.Figure 10K, L), indicating a positive effect of the 4K mutation, which amounted to an increase by a factor ~5 (eEPSC amplitude, I67N: 0.0475 ± 0.0087; I67N/4K: 0.2668 ± 0.12 nA; Mann-Whitney test, p = 0.035).
We previously created a simple mathematical model that links the release rate to the number of charges added to the SNARE complex (Ruiter et al., 2019).This model includes both spontaneous and evoked release, which are separated by the addition of 35 positive charges in the latter case (Figure 12F, black points are WT spontaneous and evoked release rates; blue line is the model fitted to WT data).Placing the spontaneous release rates for I67N and I67N/4K on this curve (by dividing the spontaneous release rate with RRP size and finding a corresponding charge value using the model) resulted in two points (red) separated by 5.6 charges, which is close to the nominal 6 charges added by the 4K mutation (Figure 12F).Similarly, evoked release in the I67N and the I67N/4K mutant neurons.(H) Example traces of mEPSC release for WT, mutant (I67N), and 1:1 co-expression of WT and SNAP25 mutant.(I) The mini frequency was decreased in both I67N mutant and the WT + I67N combination (I67N: n = 39, 36, 30 for WT, co-expressed and mutant).****p < 0.0001, ***p < 0.001, Kruskal-Wallis with Dunn's multiple comparisons.(J) mEPSC amplitudes were unchanged by the I67N mutation.(K) Example evoked excitatory postsynaptic currents (eEPSC) for WT, mutant (I67N), and co-expressed WT and mutant.(L) eEPSC amplitude was decreased by the I67N mutations (I67N: n = 39, 37, 30 for WT, co-expressed and mutant conditions, respectively).SNAP25b mutations: ****p < 0.0001, **p < 0.01, Brown-Forsythe ANOVA test with Dunnett's multiple comparisons test.(M) Overall evoked charge after a single depolarization (I67N: 24, 10, 0 for WT, co-expressed and mutant conditions, respectively).(N) Fractional contribution of the synchronous release component to the overall charge (I67N: 24, 10, 0 for WT, coexpressed, and mutant conditions, respectively).Source Data containing quantitative data are found in the Source Data files for Figures 2-4.were separated by 5.9 charges (Figure 12F).The fact that these numbers are close to the nominal 6 charges introduced shows that even for the I67N mutation, the same basal electrostatic model still applies, but the deleterious effect of the I67N on evoked release rates is larger than on spontaneous rates and therefore the rescue of evoked release rate by positive charges is insufficient to reach WT values.Moreover, because of the saturating form of the curve (i.e. the model) adding positive charges is an effective way of rescuing spontaneous, but not evoked release.Overall, the I67N disease mutation increases the amplitude of the energy barrier for fusion, and it does so more for evoked than for spontaneous release, but the electrostatic mechanism, which we assume is part of release triggering (Ruiter et al., 2019), appears to be intact.Ruiter et al., 2019) refitted to WT spontaneous and evoked data points (black points).Fitted parameters: rate 0.00029 s −1 at zero (0) charge (Z); fraction f = 0.030; the maximum rate was fixed at 6000 s −1 .WT (black points), I67N, I67N/4K (red points): means of log-transformed data.The charge values (Z, horizontal axis) for I67N and I67N/4K were found by interpolation in the model; the two spontaneous points (I67N, I67N/4K) are separated by 5.6 charges.For evoked release, rates were found by deconvolution and normalizing to RRP 0.5 (Ruiter et al., 2019).The Z-values for evoked release were found by interpolation in the model; the two mutants (I67N, I67N/4K) are separated by 5.9 charges.
The online version of this article includes the following source data for figure 12: Source data 1.Excel file containing quantitative data.

Discussion
We have shown that two SNAP25 mutations (V48F and D166Y) that compromise interaction with Syt1 lead to complex phenotypes characterized by a combination of loss-of-function and gain-of-function features.Thereby, the mutations fall into the 'neomorph' category, where the mutated protein has novel or changed interactions or functions (Verhage and Sørensen, 2020).In contrast, the I67N substitution within the SNARE bundle is a dominant-negative mutation.

SNAP25 disease mutations change protein-protein interactions and the energy landscape of fusion
Both the V48F and the D166Y resulted in a decrease in the amplitude of the apparent energy barrier for fusion, whereas the I67N increased the amplitude of the apparent fusion barrier.A vesicle's release willingness can only be assessed by fusing it; therefore, it is not possible to distinguish between effects on the fusion barrier per se, and effects on the fusion machinery.In recognition of this fact, we here refer to the 'apparent energy barrier'.Using Arrhenius' equation to convert fusion, priming and depriming rates to their respective energy barrier heights (see Materials and methods for the assumptions behind this procedure), we can derive the energy landscape for fusion of the three mutants 13A-D).This shows the multiple changes in the V48F and D166Y, which affect at least two different barriers (priming and fusion, Figure 13A, B), leading to a complex phenotype, whereas for I67N the fusion barrier is primarily (or solely) affected (Figure 13D).
When combining the I67N with the 4K mutation, which introduces six extra positive fixed charges, we could place our data within the framework of our model for electrostatic triggering (Ruiter et al., 2019) and show that the effect of charge per se is approximately the same in the I67N mutant as in the WT.Note that there are endogenous positively charged amino acids toward the C-terminal end of SNAP25 that are important for release rates (Fang et al., 2015).Rescue of spontaneous release was completed by adding six positive charges, which is consistent with the idea that the assembly of the C-terminal end of the SNARE complex, which is compromised by I67N (Rebane et al., 2018), works against the electrostatic energy barrier, which is affected by the SNARE surface charge (Ruiter et al., 2019).In contrast, rescue of evoked release by charges were incomplete, due to the larger effect of I67N on evoked release, combined with the shallow effect of charges on evoked release (Ruiter et al., 2019).The larger susceptibility of evoked release to C-terminal mutation of the SNAREs might be partly due to the higher number of SNARE complexes involved in evoked than in sustained/spontaneous release (Mohrmann et al., 2010).
The effect of the I67N mutation in the energy domain at rest can be calculated from the spontaneous release rate, which was reduced by a factor 22.4 (from 1.31 ± 0.36 to 0.0583 ± 0.0274 Hz).This corresponds to an effect in the energy domain of 3.1 k B T (where k B is Boltzmann's constant; assuming unchanged RRP size, pool normalization is not required).Work with single-molecular optical tweezers showed that the I67N mutation destabilizes the overall SNARE C-terminal and linker domain, which are supposed to deliver the power stroke for membrane fusion, by 14 k B T (Rebane et al., 2018), which is substantially more.The reduction in spontaneous release rate is more comparable to the reduction in the transition rate of folding by the C-terminal and linker domain by a factor of ~10 (Rebane et al., 2018).Since at least three SNARE complexes, possibly more, contribute to vesicle fusion (Bao et al., 2018;Manca et al., 2019;Mohrmann et al., 2010;Shi et al., 2012), folding kinetics correlates better to spontaneous fusion rates than overall SNARE complex stability.
The D166Y and V48F mutations lead to increases in spontaneous release, and more asynchronous eEPSCs, consistent with their localization in the primary SNARE:Syt1 interface (Schupp et al., 2016;Zhou et al., 2015;Zhou et al., 2017), and the demonstrated impaired Syt1 binding (Figure 8B, C), see also Alten et al., 2021.These phenotypes are at first glance similar to Syt1 knockout/knockdown (Bouazza-Arostegui et al., 2022;Chang et al., 2018;Huson et al., 2020;Ruiter et al., 2019).However, when normalized to the RRP size, mEPSC frequencies were much higher in the V48F and the D166Y than in the Syt1 KO (Figure 6D), and dual sucrose applications indicated a decrease in the amplitude of the sucrose-probed apparent energy barrier and increased release probability, features not found in the Syt1 KO (Bouazza-Arostegui et al., 2022;Huson et al., 2020).This indicates a gainof-function feature of these mutations, which fulfills D166Y > V48F; such a feature was identified as an increased interaction of V48F and D166Y with VAMP2-and syntaxin-1-containing SUVs (Figure 9D); the interaction was stronger for D166Y than for V48F.Molecular dynamics (MD) simulations of the SNAP25 mutants did not reveal major structural changes in the SNAP25 backbone compared to WT (Figure 9-figure supplement 2A-C).Nevertheless, the calculation of electrostatic (Coulomb) and van der Waals (Lennard-Jonson, LJ) interactions for residues 48-52 (Figure 9-figure supplement 2D) and residues 162-166 (Figure 9-figure supplement 2E) shows that for D166Y the interaction energy is substantially more negative.This suggests that the interaction between D166Y and nearby residue H162 is stronger than in WT.Overall, this may result in stabilization of a structure consistent with SNARE complex formation and explain why D166Y is a stronger gain-of-function mutation than V48F.However, note that the AlphaFold prediction (Jumper et al., 2021;Mirdita et al., 2022) used as a starting point for these simulations is identical to the structure of SNAP25 in the assembled SNARE complex (Zhou et al., 2015), whereas unassembled SNAP25 is likely less structured or unstructured (Fasshauer et al., 1997).Overall, these disease mutants do not only fail in their interaction with Syt1, they bypass fusion control, resulting in premature SNARE complex assembly.Although the increased spontaneous release rate might cause patient symptoms, Alten et al., 2021 suggested that the resulting postsynaptic depolarization might compensate for the smaller eEPSC amplitude to normalize the overall firing rate.
In an SUV:GUV fusion assay, where folding of SNAP25 onto syntaxin-1/VAMP2 is limiting, D166Y and V48F caused an increase in pre-stimulation fusion rates.This is consistent with recent data showing that folding of SNAP25 onto a template formed by VAMP2 and syntaxin-1 held in place by Munc18-1 might be a late, rate-limiting, step in exocytosis (Jiao et al., 2018).This process is regulated and sped up by Munc13-1 (Kalyana Sundaram et al., 2021;Shu et al., 2020;Wang et al., 2019).The fact that a similar effect on spontaneous fusion was seen in the assay with preformed t-SNAREs in the presence of complexin indicates that there is an assembly step, even in that assay, which can be sped up by the mutations.This aligns with the demonstration by single-molecule FRET of a further assembled (tighter) state of the trans-SNARE complex induced by Ca 2+ -unbound Syt1, which becomes committed for fusion once Ca 2+ binds (Das et al., 2020).
V48F and D166Y change the size of the RRP via effects on priming, depriming and fusion We found that V48F and D166Y cause a decrease in RRP size, whether measured by sucrose application or by train stimulation.The reduced RRP appears to be at variance with the publication by Alten et al., 2021, who used sucrose application to larger mixed cultures.There might be several reasons for this discrepancy.First, Alten et al. used longer duration of sucrose application (~1 min) and applied sucrose to a large mixed culture, which might result in a variable delay such that all synapses are not stimulated simultaneously.This will result in temporally overlapping release of RRP and upstream vesicle pools, which then cannot be distinguished from each other.When grown on 50 µm microislands, sucrose can be applied acutely (within ~0.05 s) to the entire dendritic tree and all synapses using a local perfusion system, which allows distinguishing the RRP from upstream pools.Second, Alten assessed the RRP for GABAergic synaptic transmission in mixed culture, whereas we measured the RRP for glutamatergic neurons in autaptic culture.The RRP has a different substructure in the two types of neurons (Moulder and Mennerick, 2005), which might lead to different findings.Third, it was recently shown that sufficient neuronal maturation is necessary to detect the decrease in RRP upon Syt1 elimination (Bouazza-Arostegui et al., 2022), which is likely modulated by the presence of other neurons (Chang et al., 2018;Liu et al., 2009;Wierda and Sørensen, 2014).This adds an additional layer of complexity since neuronal maturation likely varies between laboratories.
Dissection of the three rates that determine RRP size (priming, depriming, and fusion rates) showed that all three are changed in the V48F and D166Y mutations.Especially for D166Y, but also for V48F, spontaneous release contributed to RRP depletion by triggering premature fusion.Significant RRP depletion has not been expected for moderate increases in mEPSC frequency (Rhee et al., 2005;Ruiter et al., 2019), based on the argument that RRP refilling should be fast enough to counteract depletion.However, to properly make this argument, priming, depriming, and spontaneous release rates must be compared in the steady-state situation.Similarly, we recently showed by cryo-electron tomography that the 4K mutation, which increased the mEPSC frequency to ~30 Hz and had a reduced RRP (Ruiter et al., 2019) caused a loss of synaptic vesicles tethered to the membrane with three tethers (Radecke et al., 2023), which is the structural correlate of the RRP (Fernández-Busnadiego et al., 2010).
The major effect on RRP size is caused by a reduction in forward priming rate by D166Y and V48F.Comparison to the Syt1 KO showed qualitatively similar changes, but of a smaller magnitude, with spontaneous release playing a negligible role for RRP size.Notably, in all three cases a reduction in depriming rate partly counteracted the lowered priming rate (Table 1) -this was only significant for D166Y.Indeed, part of the role of Syt1:SNARE interaction might be catalytic, lowering the energy level of a transition state along the path to priming, which will affect both rates (Walter et al., 2013).This might happen because transient binding to Syt1 might structure SNAP25 and assist in formation of the SNARE complex, whereas SNAP25 mutants might prestructure the protein, bypassing the need for Syt1.In an energy diagram (Figure 13), it becomes clear that vesicle priming and regulation of spontaneous release are interdependent.Stabilization of the RRP state will both increase RRP size at rest and reduce spontaneous release, since lowering the energy level of (i.e.stabilizing) RRP vesicles will increase the size of the energy barrier that the RRP vesicles face (Figure 13).Consistently, downregulation (clamping) of spontaneous release and upregulation of evoked release are often interdependent under conditions where Syt1 expression level, or Syt1 interaction with the SNAREs, are up-or downregulated (Courtney et al., 2021;Courtney et al., 2019;Vevea and Chapman, 2020;Zhou et al., 2015;Zhou et al., 2017).This does not rule out the of mutations that can affect one mode of release more than the other (e.g. the I67N/4K mutation).By inference, assembly of the primary Syt1:SNARE interface (Zhou et al., 2017) is most likely involved in both clamping release and setting up an RRP.In further support of this, the minimal in vitro assay with preassembled t-SNARE dimers displayed a qualitatively similar reduction in calcium-dependent release, with D166Y being more impaired than V48F (Figure 8D).Overall, the energetic contribution of Syt1:SNARE interaction to Ca 2+ -triggered release is a stabilization of upstream steps, and the increase of the fusion barrier downstream of the RRP (Figure 13).The electrostatic nature of the fusion barrier (Ruiter et al., 2019) ensures its rapid dissolution by Ca 2+ , possibly by unbinding of Syt1 from the SNARE complex (Voleti et al., 2020).
For the I67N mutation, the sucrose RRP was also reduced in size when using 0.5 M sucrose, which is consistent with previous observations (Alten et al., 2021).However, when the energy barrier for release is increased, 0.5 M sucrose is insufficient to deplete the RRP (Schotten et al., 2015).Accordingly, 0.75 M sucrose released an RRP of similar size in WT, I67N, and I67N + WT co-expressing cells, showing that the vesicle priming reaction is intact, but the vesicles face a larger fusion barrier (Figure 11).In this mutant, train stimulations in the I67N or co-expressed WT + I67N case resulted in a phenotype quite distinct from the WT, with strong facilitation throughout the train.Thus, the I67N is a strongly dominant-negative mutant, whether considering mEPSC frequency, eEPSC amplitude, or train stimulations, whereas V48F and D166Y are incompletely dominant.This can be explained by the different function of SNAP25 domains during fusion, where V48 and D166 help set up an arrested primed vesicle state by interacting with Syt1 (Schupp et al., 2016;Zhou et al., 2015;Zhou et al., 2017), whereas I67 participates in the final conformational change, the 'stroke' that leads to assembly of the C-terminal end of the complex and the linker domain.Due to their defect in priming, V48F and D166Y might not enter the super-complex (the complex of SNARE complexes driving fusion) as often as WT protein.In contrast, I67N will readily enter the super-complex and compromise its function, which will lead to a dominant-negative phenotype due to the multiple SNARE complexes involved in fusion (Bao et al., 2018;Mohrmann et al., 2010;Shi et al., 2012).
Back-extrapolation of 40 Hz trains (@4 mM Ca 2+ ) also led to the conclusion that V48F and D166Y have a reduced RRP ev , the RRP sub-pool that action potentials draw on, and a reduced forward priming rate.Interestingly, the release probability when normalized to the sucrose pool (RRP 0.5 ) was increased non-significantly for the V48F and significantly for the D166Y mutation, but when considering 40 Hz trains, the release probabilities of both mutations were decreased, consistent with a shift toward facilitation.The latter finding is likely caused by the defective interaction with Syt1, which leads to suboptimal priming and/or defective super-priming, which is an additional priming step after entry of the vesicle into the RRP (Lee et al., 2013;Taschenberger et al., 2016).This will lead to reduction of the first eEPSC of a train and thereby a lowered nominal release probability.However, when comparing the eEPSC charge to the RRP 0.5 , the strong reduction in the RRP 0.5 pool (especially in the D166Y mutation) accounts for the increase in the overall release probability.This can be explained if spontaneous release causes a disproportional depletion of vesicles, which are in the RRP 0.5 , but not in the RRP ev .

Conclusion
SNAP25 disease missense mutations change the function of the protein without compromising its expression, leading to dominant negative or neomorphic mutations.Missense mutations in the primary SNARE:Syt1 interface (V48F and D166Y) result in a complex phenotype characterized by loss-of-function in the priming step and gain-of-function in the fusion step.Missense mutation in the SNARE bundle (I67N) leads to an increased amplitude of the energy barrier for fusion.In addition, disease mutations display inefficient rescue of neuronal survival.Overall, SNAP25 encephalopathy caused by single missense point mutations presents with interdependent functional deficits, which must be overcome for successful treatment.

Animals
Snap25 KO (Washbourne et al., 2002) and Synaptotagmin-1 (Syt1) KO (Geppert et al., 1994) C57BL/6 mice: Heterozygous animals were routinely backcrossed to C57BL/6J to generate new heterozygotes.The strain was kept in the heterozygous condition and timed were used to recover KO embryos by cesarean section at embryonic day 18 (E18).Pregnant females were killed by cervical dislocation; E18 embryos of either sex were collected and killed by decapitation.Permission to keep, breed, and use Snap25 and Syt1 mice was obtained from the Danish Animal Experiments Inspectorate (2018-15-0202-00157) and followed institutional guidelines as overseen by the Institutional Animal Care and Use Committee (IACUC).CD1 outbred mice were used to create astrocytic cultures and mass cultures for Western blotting.Newborns (P0-P2) of either sex were used.Pups were killed by decapitation.

Cell lines
HEK293-FT cells for production of lentiviruses were obtained from the Max-Planck-Institute for biophysical chemistry.The cells were passaged once a week, and they were used between passages 11 and 25 for generation of lentiviral particles.The cells were kept in Dulbecco's modified Eagle medium (DMEM) + Glutamax (Gibco, cat.31966047) supplemented with fetal bovine serum (Gibco, cat.10500064), Pen/Strep (Gibco, cat.15140122), and Geneticin G418 (Gibco, cat.11811064) at 37°C in 5% CO 2 .Negative mycoplasma status was confirmed using a commercial kit (Venor GeM OneStep,Minerva biolabs,.
Glass coverslips were washed overnight in 1 M HCl; for an hour in 1 M NaOH and washed with water before storage in 96% ethanol.Coverslips were first coated by 0.15% agarose and islands were made by stamping the coating mixture (3 parts acetic acid [17 mM], 1 part collagen [4 mg/ml], and 1 part poly-D-lysine [0.5 mg/ml]) onto the glass coverslips using a custom rubber stamp.Glial cells were washed with prewarmed HBSS-HEPES.Trypsin was added and the flasks were incubated at 37°C for 10 min.Cells were triturated and counted in a Bürker chamber before plating onto the glass coverslip with DMEM solution.After 2-5 days, neurons were plated on the islands.
Neurons for autaptic culture were isolated from E18 Snap25 KO mice of either sex.Pups were selected based on the absence of motion after tactile stimulation and bloated neck (Washbourne et al., 2002); the genotype was confirmed by Polymerase Chain Reaction (PCR) in all cases.The cortices were isolated from the brains in the HBSS-HEPES medium.The meninges were removed and the hippocampi were cut from the cortices before being transferred to 0.25% trypsin dissolved in HBSS-HEPES solution.The hippocampi were incubated for 20 min at 37°C, subsequently washed with HBSS-HEPES, triturated and the cell count was determined with a Bürker chamber before plating on the islands (7000-8000 neurons per well).Cells were incubated in the NB medium (Neurobasal with 2% B-27, 1 M HEPES, 0.26% Glutamax, 14.3 mM β-mercaptoethanol, 20,000 IU penicillin, 20 mg streptomycin) and used for measurements between DIV10 and 14.
Hippocampal neurons for high-density cell culture for Western blotting were obtained from P0 to P1 CD1 outbred mice.The dissection, tissue digestion, and cell counting were performed the same way as the neurons for the autaptic culture, the high-density culture (600,000 neurons per well) was then kept in NB-A medium (Neurobasal-A with 2% B-27, 1 M HEPES, 0.26% Glutamax, 20,000 IU penicillin, 20 mg streptomycin) and half the volume of the media was replaced every 2-3 days before harvesting the cells at DIV14.

Constructs for rescue experiments
SNAP25b was N-terminally fused to GFP and cloned into a pLenti construct with a CMV promoter (Delgado-Martínez et al., 2007).Mutations were made using the QuikChange II XL kit (Agilent).
Primers were ordered from TAGC Copenhagen.All mutations were verified by sequencing before virus production.Viruses were prepared previously described using transfection of HEK293FT cells (Naldini et al., 1996).Neurons were infected with lentiviruses on DIV 0-1, 30 µl total virus per well (WT or mutant; 15 µl + 15 µl WT + mutant virus for co-expressed condition).

Immunostaining and confocal microscopy
Autaptic hippocampal neurons were fixed at DIV14 in 2% paraformaldehyde (PFA) in culture medium for 10 min and subsequently in 4% PFA for 10 min.Cells were then washed with phosphate-buffered saline (PBS), permeabilized by 0.5% Triton X-100 for 5 min, and blocked with 4% normal goat serum in 0.1% Triton X-100 (blocking solution) for 30 min.Cells were incubated with primary antibodies diluted in blocking solution (anti-MAP2, 1:500, chicken, ab5392, Abcam; and anti-vGlut1 1:1000, guinea pig, AB5905, Merck Millipore) for 2 hr at room temperature (RT).After washing with PBS, the cells incubated with secondary antibodies in blocking solution for 1 hr at RT in the dark (anti-chicken Alexa 568, 1:1000, A11041, Thermo Fisher Scientific; and anti-guinea pig Alexa 647, 1:1000, A-21450, Thermo Fisher Scientific) and washed again.Coverslips were mounted with FluorSave and imaged on Zeiss CellObserver spinning disc confocal microscope (×40 water immersion objective; NA 1.2) with Zeiss Zen Blue 2012 software.Images were acquired as Z-stack and 9 images per plane to capture the whole island in the field of view.The images were post-processed with Zeiss Zen Black software and neuronal morphology was analyzed using SynD automated image analysis (Schmitz et al., 2011).

Western blotting
Harvesting the high-density hippocampal culture, Bicinchoninic acid (BCA) assay and transferring the protein samples on a polyvinylidene difluoride (PVDF) membrane were performed as described before (Ruiter et al., 2019).Incubation in primary antibodies (a-SNAP25: mouse, 1:10,000, SYSY 111011, Synaptic Systems; a-VCP: mouse, 1:2000, ab11433, Abcam) was performed overnight with 70 rpm shaking at 4°C, followed by washing in Tris-buffered saline with 0.1% Tween 20 Detergent (TBST) and a 1-hr incubation in secondary antibody (goat a-mouse-HRP: 1:10,000, P0447, Dako).After washing, Pierce ECL Western blotting substrate was added and chemiluminescence was visualized with FluorChem E (Proteinsimple).The Western blots were quantified using ImageJ (1.52a) to measure the signal intensity of the protein bands relative to the loading control (VCP).All relative levels of the target protein were normalized to the average relative level of the same target protein in the WT samples.
where E A is the activation energy, and T are the gas constants and the absolute temperature, respectively, and A is an empirical constant that depends on collision rates (Schotten et al., 2015).Solving for the activation energy, we get: Since A is unknown, we cannot use this equation to calculate the absolute values of the transition energies; however, when we compare a mutation to the WT condition, and under the assumption that A is unchanged by mutation, we can calculate the difference in energy level of the transition states: Using this equation sequentially, for the priming rate, the depriming rate and the fusion rate, we can derive the entire energy diagram, under the additional assumption that the energy in the pre-primed state is identical between WT and mutation, and using that at room temperature, RT = 2.479 kJ/mol.
The assumption that the empiric factor, A, is unchanged by mutation is likely to hold for the fusion reaction, which depends on conformational changes in a preformed complex.In contrast, collision rates might be involved in priming; in that case, the effect of the V48F and D166Y mutations on priming, which we here attribute to an increase in the energy level of the priming transition state might reflect a lower collision rate between vesicles and plasma membrane fusion machinery, and/or a lower energy level of the pre-primed state.

Protein expression and purification
In general, the expression vectors, encoding the desired protein constructs were transfected into Escheria coli BL21 (DE3) (Stratagene).At an OD 660 of 0.8, protein expression was induced by the addition of 0.3 mM Isopropyl β-D-1-thiogalactopyranoside (IPTG).Alternatively, proteins were expressed by autoinduction using buffered media containing lactose (Studier, 2005).Cells were harvested by centrifugation (3500 rpm, 15 min, H-12000 rotor, Sorvall) and lysed using the high-pressure pneumatic processor 110 L (Microfluidics).Cell fragments were removed by centrifugation at 60,000 rpm (70Ti rotor, Beckman) for 1 hr and the clarified supernatant was snap frozen in liquid nitrogen.
The purification of full-length VAMP2 was performed as described previously (Kedar et al., 2015) with the following modifications.Cells were grown in ZYM media (Studier, 2005) and protein expression was induced with IPTG for 3 hr at 25°C.The purification and expression of the GST-tagged cytosolic domain of VAMP2 were described previously (Ruiter et al., 2019).Synaptotagmin-1-His6 lacking the luminal domain was purified as described previously (Malsam et al., 2012) with the following modification.After dilution to 50 mM salt, the protein was further purified on a MonoS Sepharose column (GE applying a gradient of 50-500 mM KCl in 25 mM HEPES-KOH (pH 7.4).
His6-syntaxin-1A purification was performed as outlined by Schollmeier et al., 2011 with the following modifications.Briefly, cells were grown in ZYM media (Studier, 2005) and autoinduction was used for the protein expression at 22°C overnight.Syntaxin-1A was eluted from Ni-NTA beads (QIAGEN) by overnight cleavage with Prescission protease (GE Healthcare) at 4°C, removing the His6 tag.After dilution to 80 mM salt, the protein was further purified on a MonoQ Sepharose column (GE Healthcare) applying a gradient of 50-500 mM KCl in 25 mM HEPES-KOH (pH 7.4).
His6-SNAP25 was expressed as depicted for syntaxin-1A and purified via Ni-NTA beads, followed by MonoQ Sepharose column chromatography (Ruiter et al., 2019).Preassembled full-length t-SNARE complexes were expressed and purified as described previously (Weber et al., 1998).His6-Complexin II expression and purification were performed according to Malsam et al., 2012 with the following modifications.His6-CpxII was expressed in BL21-DE3 codon + bacteria for 2 hr at 27°C.
The concentrations of purified proteins were determined by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and Coomassie Blue staining using Bovine Serum Albumin (BSA) as a standard and the Fiji software for quantification.Furthermore, mutant and WT protein concentrations were directly compared on a single gel.
Proteins were reconstituted as described previously (Malsam et al., 2012).For the docking and fusion assays, t-SNARE WT and mutants were reconstituted at a protein to lipid ratio of 1:900.For the syntaxin-1A membrane fusion assay, syntaxin-1A was reconstituted at a protein to lipid ratio of 1:1000.Briefly, 5 µmol dried lipids were dissolved in 0.7 ml reconstitution buffer (25 mM HEPES-KOH, pH 7.4, 550 mM KCl, 1 mM Ethylenediaminetetraacettic acid (EDTA-NaOH), 1 mM dithiothreitol (DTT)) containing final 1.4 wt% octyl-β-D-glucopyranoside (β-OG) and either 6.5 nmol (390 µg) t-SNARE complex or 5 nmol (165 µg) syntaxin-1A.SUVs containing either t-SNARE or syntaxin-1A were formed by rapid β-OG dilution below the critical micelle concentration by adding 1.4 ml reconstitution buffer.For the quantification of lipid recovery, a 20-µl aliquot (GUV input) was removed and stored at −20 °C.The liposome suspension was desalted using a PD10 column (GE Healthcare) equilibrated with desalting buffer 1 (1 mM HEPES-KOH, pH 7.5, 1 wt/vol% glycerol, 10 µM EGTA-KOH, 1 mM DTT) and snap frozen in four aliquots in liquid nitrogen and stored in −80°C.These syntaxin-1A and t-SNARE SUVs were later used to prepare GUVs.The final protein to lipid ratios were determined by SDS-PAGE and Coomassie Blue staining of the proteins and Atto647 fluorescence intensity measurements of the lipids.
To determine the respective inputs, 28 nmol GUV lipids or 5 nmol SUV lipids were treated with 1% SDS/1% DDM (final) and the corresponding Atto647 and Atto488 fluorescence was measured at λ ex = 620/40 nm, λ em = 680/30 nm and λ ex = 485/20 nm, λ em = 528/20 nm, respectively.A sample lacking SUVs was used to determine the GUV recovery (usually 80-95%).GUV recovery of each sample was used to normalize the respective SUV docking.A sample without GUVs was used to determine the absolute background (usually <15%), which was subtracted from all samples.

SNAP25 recruitment assay/SUV flotation assay
SUVs containing 35 pmol syntaxin-1A or 210 pmol VAMP2 and/or 100 pmol Synaptotagmin-1 (66 nmol for Stx1-, V2/Syt1-, or Syt1-SUVs, 40 nmol lipid for V2-SUVs) were mixed with 180 pmol SNAP25 in a final assay volume of 50 µl in fusion buffer (20 mM MOPS-KOH, pH 7.4, 135 mM KCl, 1 mM DTT).After 2 hr on ice, allowing the SNARE complex formation, the samples were diluted five times with fusion buffer and mixed with the equivalent volume of 80% nycodenz solution.After overlaying the sample with 100 µl 35% nycodenz, 25 µl 20% nycodenz solution, and finally with 5 µl fusion buffer in an ultra-clear tube (5 × 41 mm), the liposomes were isolated by centrifugation for 3 hr 40 min at 55,000 rpm at 4°C in a SW 60 rotor (SW 60 Ti, Beckman).20 µl were harvested from the top of the gradient and mixed with 8 µl of 4× Laemmli buffer (final 62.5 mM Tris-HCl, pH 6.8, 10% glycerol, 2% SDS, 50 mM β-mercaptoethanol, 0.1% bromphenol blue).18 µl of this mixture were used to quantify the amount of recruited SNAP25b by SDS-PAGE followed by Coomassie Staining and Silver Staining.Using the Fiji software (ImageJ based) the Coomassie stained band intensity of syntaxin-1A or VAMP2 or Syt1, respectively, were determined and normalized to the respective mean.Subsequently, Silver Staining was used to quantify the band intensities of SNAP25b, and these values were normalized to the intensities of the relative protein (e.g.Syntaxin-1A) based on the Coomassie Staining.From this, the ratios between the SNAP25b mutants and WT were determined.

Silver Staining
Coomassie prestained gels from SDS-PAGE were destained overnight in destain solution (30% methanol, 10% acetic acid).Gels were gently washed 30 min with 10% ethanol and 1 min with 0.02% sodium thiosulfate.After short washing with deionized water, gels were stained (0.03% PFA, 0.002% silver nitrate) for 15 min and again quickly washed with water.Developing solution (0.06% sodium carbonate, 0.018% PFA, 0.0002% sodium thiosulfate) was applied to the gels until protein band intensities were satisfactorily stained.The reactions were stopped by replacing the staining solution with 0.07% acetic acid.Gels were scanned and quantified by using the Fiji (ImageJ based) software.

Atomistic MD simulations
Simulations were carried out at the atomic level using classical MD.We used the CHARMM36m force field for protein, CHARMM TIP3P for water, and the standard CHARMM36 for ions (Huang et al., 2017).GROMACS 2022 software package was used for performing these simulations (Abraham et al., 2015).The 3D structure of SNAP25 WT protein was generated through AlphaFold2, utilizing the first ranked structure obtained using the default settings of ColabFold v1.5.2 (Mirdita et al., 2022).Both WT and the V48F and D166Y mutant variants' topology files and initial 3D structures, inclusive of water and ions, were produced via the CHARMM-GUI web interface (Lee et al., 2016).All systems were hydrated, neutralized with counter ions, and supplemented with 150 mM potassium chloride to replicate experimental conditions.Following the CHARMM-GUI recommended protocol, systems were energy-minimized, equilibrated under NpT conditions, temperature-stabilized at 310 K by the Nose-Hoover thermostat with a 1.0 ps time constant (Evans and Holian, 1985), and maintained a constant 1 atm pressure via the Parrinello-Rahman barostat, setting the time constant at 5.0 ps and an isothermal compressibility at 4.5 × 10 −5 bar −1 (Parrinello and Rahman, 1981).Isotropic pressure coupling was utilized in the simulations.The Verlet scheme determined neighbor searches, updating once every 20 steps (Verlet, 1967).Electrostatics were using the Particle Mesh Ewald method (Darden et al., 1993) with parameters set to a 0.12-nm spacing, a tolerance 10 −5 , and a 1.2-nm cutoff.Periodic boundary conditions were applied in all three dimensions.Simulations ran with a 2-fs timestep until 800 ns was achieved.As shown in Figure 9-figure supplement 2 panel A, the structures used for protein alignment represent the predominant structure from the most populated cluster.Clustering was based on the RMSD value using the GROMACS 'gmx cluster' tool and the Gromos algorithm (Daura et al., 1999), setting an RMSD cutoff for neighbor structures at 0.6 nm (Legrand et al., 2020).The RMSD analysis was conducted using the average structure of the WT as a reference.Interaction energies, including both short-range Coulomb and Lennard-Jones forces, were computed throughout the 800 ns trajectory.The autocorrelation function of the data indicated that correlations diminished substantially after approximately 6 ns.Given this reduced correlation, we adopted a block size of 200 ns to ensure statistical independence between blocks.
All structures generated using AlphaFold2, as well as the initial structures and topology files for atomistic MD simulations, were deposited to ZENODO (https://doi.org/10.5281/zenodo.10051665).

Statistical analysis
Graphs (bar and line) display mean ± standard error of the mean with all points displayed, except when otherwise noted; for electrophysiological experiments n denotes the number of cells recorded and is given in the legends.For in vitro experiments the number of biological replicates was 3 unless stated otherwise in the legends.Statistics were performed using GraphPad Prism 9. Unless otherwise noted, statistical differences between several groups were determined by one-way analysis of variance (ANOVA); post-test was Dunnett's test comparing to the WT condition as a reference, unless otherwise noted.Equal variance of groups was tested by the Brown-Forsythe test; in case of a significant test, the Brown-Forsythe ANOVA test, which does not assume equal variances, was used instead.Kruskal-Wallis test was used in cases, where the data structure contained many identical values (zeros).Pairwise testing was carried out using unpaired t-test, or Welch's t-test in case of significantly different variances as determined by an F-test.The test is mentioned in the figure legend; if no test is mentioned, the difference was not significantly different.Significance was assumed when p < 0.05 and the level of significance is indicated by asterisks: *p < 0.05, **p < 0.01, ***p < 0.001, ****p <0.0001.

Materials availability
New materials generated during this project are available upon request.
The following dataset was generated:

Figure 1 .Figure 2 .
Figure 1.Localization of three pathogenic mutations in SNAP25.(A) Schematic of the neuronal SNARE complex interacting with C2B domain of synaptotagmin-1 (Syt1; not to scale) via the primary interface.Position of the I67N mutation in the first SNARE domain of SNAP25 is depicted by an asterisk.(B) Interaction site of the C2B domain of Syt1 and SNAP25.Syt1 interacts with SNAP25 both electrostatically (regions I and II) and within the hydrophobic patch (HP patch)(Zhou et al., 2015).(C) Position of the disease-linked mutations V48F (hydrophobic patch) and D166Y (region I) in the SNARE complex.

Figure 2
Figure 2 continued on next page Figure2 continued

Figure 6 .
Figure 6.Dissection of the readily releasable pool (RRP) reduction in V48F and D166Y mutations.(A) One-pool model of the RRP.k 1 is the rate of priming (units vesicles/s), k −1 is the rate of depriming (s −1 ), and k f is the rate of fusion (s −1 ).(B) Estimation of the three parameters from the response to 0.5 M sucrose and a measurement of the spontaneous release rate.(C) Variance-mean analysis in 50-ms intervals during the sucrose application allows determination of the corrected baseline by back-extrapolation of a regression line to the variance of the baseline.(D) Normalized mEPSC frequency (k f ) for V48F, D166Y, and Syt1 knockout (KO) (V48F: n = 23, 24 for wildtype [WT] and mutant conditions, respectively; D166Y: n = 19, 19; Syt1: n = 23, 26).Brown-Forsythe analysis of variance (ANOVA) test with Dunnett's multiple comparison test, testing the three mutant conditions against each other.****p < 0.0001, ***p < 0.001.(E) Normalized RRP size for WT and mutant conditions, with indications of the effect of the mutant-induced changes in k 1 , k −1 , and k f on the RRP size.The online version of this article includes the following source data and figure supplement(s) for figure 6: Source data 1.Excel file containing quantitative data.

Figure supplement 1 .
Figure supplement 1.Effect of sucrose stimulation on estimates of k 1 and k −1 .

Figure
Figure supplement 1-source data 1.Excel file containing quantitative data.
the reaction.Shown is mean ± SEM (n = 3).****p < 0.0001; **p < 0.01, ANOVA with Dunnett's multiple comparisons test, comparing each mutation to the corresponding WT condition.The online version of this article includes the following source data for figure8:Source data 1.Excel file containing quantitative data.

Figure supplement 1
Figure supplement 1-source data 1.Original files for the analysis by Coomassie and silver stained gels.

Figure supplement 1
Figure supplement 1-source data 2. PDF containing original pictures of gels with highlighted bands and sample labels.

Figure supplement 2 .
Figure supplement 2. Molecular dynamics simulations of mutants.

Figure supplement 2
Figure supplement 2-source data 1.Excel file with quantitative data.

"Figure 13 .
Figure13.Energy landscapes.The energy landscapes of wildtype (WT) and mutants were calculated as explained in Materials and methods and displayed to scale.Energy landscapes for D166Y (A), V48F (B), and Syt1 knockout (KO) (C) are shown at rest and are characterized by a higher priming barrier ('loss-of-function' phenotype), a destabilized readily releasable pool (RRP), and a lower fusion barrier ('gain-of-function' phenotype).The I67N (D) is characterized by a higher fusion barrier ('loss-of-function' phenotype).The relative increase in the fusion barrier by the I67N mutation is higher during stimulation than at rest.Dotted lines represent energy levels for which less is known.