Munc18-1 catalyzes neuronal SNARE assembly by templating SNARE association

  1. Junyi Jiao
  2. Mengze He
  3. Sarah A Port
  4. Richard W Baker
  5. Yonggang Xu
  6. Hong Qu
  7. Yujian Xiong
  8. Yukun Wang
  9. Huaizhou Jin
  10. Travis J Eisemann
  11. Frederick M Hughson  Is a corresponding author
  12. Yongli Zhang  Is a corresponding author
  1. Yale University School of Medicine, United States
  2. Princeton University, United States
  3. University of California, San Diego, United States
10 figures, 2 videos, 2 tables and 1 additional file

Figures

Two potential pathways for Munc18-1-regulated neuronal SNARE assembly.

(i) Munc18-1 first serves as a syntaxin chaperone and binds syntaxin to inhibit its association with other SNAREs. (ii) Closed syntaxin is opened by Munc13-1, a large multifunctional protein that …

https://doi.org/10.7554/eLife.41771.002
Figure 2 with 7 supplements
Single-molecule manipulation based on optical tweezers revealed a ternary template complex.

(A) Experimental setup and structural model of the template complex. Some key mutations tested in this study are indicated by dots: red (phosphomimetic mutations) or gray (others) for Munc18-1, …

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

MATLAB figure for the FECs shown in Figure 2C.

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

Complete time-dependent instantaneous force, extension, and trap separation obtained in a representative single-molecule experiment in the presence of WT 2 µM Munc18-1.

Data here correspond to Figure 2C (FECs #2–3) and Figure 2—figure supplement 2.

https://doi.org/10.7554/eLife.41771.012
Figure 2—figure supplement 1
Sequences, domains, and crosslinking sites of the SNARE proteins used in this study.

The amino acids in hydrophobic layers (from −7 to +8) and the central ionic layer (0 layer) are colored yellow. The underlined sequences are added to facilitate crosslinking of Qa and R SNAREs, …

https://doi.org/10.7554/eLife.41771.004
Figure 2—figure supplement 2
Time-dependent extension (top panel), force (middle panel), and trap separation (bottom panel) for a typical experiment to test template complex formation (Figure 2—source data 2).

Data here and FECs #2 and #3 in Figure 2C are acquired on the same Qa-R SNARE conjugate, with the same pulling round numbering. Close-up views of different time regions indicated by A-D are shown: (A

https://doi.org/10.7554/eLife.41771.005
Figure 2—figure supplement 3
Extension-time trajectories at two constant mean forces showing the opening-closing transition of the partially closed syntaxin molecule.
https://doi.org/10.7554/eLife.41771.006
Figure 2—figure supplement 4
Force-dependent syntaxin opening probabilities (top panel) and opening and closing rates (bottom panel) obtained by pulling syntaxin from the two N-terminal sites, R198C and I187C (Figure 2—figure supplement 1).

Curves are best model fits to derive the energies and kinetics at zero force associated with the transitions, with solid curves for unfolding and dashed curves for folding. The unfolding energies of …

https://doi.org/10.7554/eLife.41771.007
Figure 2—figure supplement 5
FECs of Qa only (#1) or the Qa-R SNARE conjugate (other FECs) pulled from Site I187C in the presence of 2 μM Munc18-1 without (-) or with (+) 60 nM SNAP-25B.

The wide-type (‘WT’) syntaxin-1 here denotes syntaxin-1A (a.a. 1–265, I187C, C145S; see ‘SNARE protein constructs’), with additional mutations indicated. FECs in each bracket were obtained on the …

https://doi.org/10.7554/eLife.41771.008
Figure 2—figure supplement 6
Extension-time trajectories at two constant mean forces showing the opening-closing transition of the syntaxin molecule pulled from the crosslinking site I187C (Figure 2—figure supplement 1).

The red curves are idealized state transitions derived from hidden-Markov modeling. State 6' represents the fully closed syntaxin (Figure 2B).

https://doi.org/10.7554/eLife.41771.009
Figure 2—figure supplement 7
Extension-time trajectories showing conformational transitions of the template complex transition pulled from Site I187C in the absence (top trace) and presence (bottom) of 60 nM SNAP-25B.

The unfolding energy of the template complex crosslinked at I187C is estimated to be 4.8 ± 0.3 kBT, close to the unfolding energy of 5.2 ± 0.1 kBT of the template complex crosslinked at R198C. …

https://doi.org/10.7554/eLife.41771.010
Figure 3 with 3 supplements
Stability, conformation, and folding kinetics of the template complex.

(A, D, E, I) Extension-time trajectories at constant mean forces with the WT template complex (A) or its variants containing indicated mutations in Munc18-1 (D), VAMP2 (E), or syntaxin (I). The red …

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

MATLAB figure corresponding to Figure 3A with an additional trace at force F = 5.0 pN.

https://doi.org/10.7554/eLife.41771.019
Figure 3—source data 2

MATLAB figure containing expanded traces shown in Figure 3D,E,I,H.

https://doi.org/10.7554/eLife.41771.020
Figure 3—source data 3

MATLAB figure for Figure 3C.

https://doi.org/10.7554/eLife.41771.021
Figure 3—figure supplement 1
FECs obtained in the presence of 2 µM Munc18-1.

Dashed blue ovals mark template complex transitions. Note that the partially closed syntaxin state was abrogated by modifications that are known to destabilize the closed syntaxin, including …

https://doi.org/10.7554/eLife.41771.016
Figure 3—figure supplement 2
Circular Dichroism (CD) spectra show that the mutations tested in our experiments do not significantly alter Munc18-1 folding.

The CD spectra of Munc18-1 mutants that abolished or weakened the template complex are shown, including Munc18-1 F-pocket mutations L247R, T248G, L247A/T248G, disease-related mutations L341P and …

https://doi.org/10.7554/eLife.41771.017
Figure 3—figure supplement 3
Snapshots of the extension-time trajectories at constant mean forces showing sporadic folding of the template complex in the absence of syntaxin NRD.
https://doi.org/10.7554/eLife.41771.018
Stability of the template complex correlates with SNARE-mediated membrane fusion and neurotransmitter release.

(A) Unfolding energy of the WT and mutant template complexes; see also Table 1. Unfolding energy that is less than our detection limit (1.5 kBT) is plotted as zero. The unfolding energy is derived …

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

Data summary table for the results shown in Figure 4.

https://doi.org/10.7554/eLife.41771.023
Spontaneous SNARE assembly in the absence of Munc18-1 is inefficient.

(A) Representative FECs obtained by consecutively pulling and relaxing a single Qa-R SNARE conjugate in the presence of 60 nM SNAP-25B (Figure 5—source data 1). No SNARE assembly is observed. (B) …

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

MATLAB figure corresponding to Figure 5A.

https://doi.org/10.7554/eLife.41771.025
Figure 5—source data 2

MATLAB figure corresponding to Figure 5B (FEC#1–4).

https://doi.org/10.7554/eLife.41771.026
Figure 5—source data 3

MATLAB figure corresponding to Figure 5C.

https://doi.org/10.7554/eLife.41771.027
Comparison of SNARE assembly in the absence and presence of Munc18-1.

Bars indicate probabilities of Munc18-1-independent or spontaneous SNARE assembly (blue), Munc18-1-chaperoned SNARE assembly (red), and SNARE misassembly (black). See also Figure 6—source data 1.

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

Data summary table for the results shown in Figure 6.

Dataset 1. This dataset contains 18 MATLAB figures corresponding to Figures 210, as are listed below. Each plot or curve in a MATLAB figure has its associated data embedded. The data values for all data points in a plot can be obtained using MATLAB command get(plot_handle,’xdata’) or get(plot_handle,’ydata’), where plot_handle is a unique identifier of the plot.

https://doi.org/10.7554/eLife.41771.029
Figure 7 with 2 supplements
Template complex facilitates SNARE assembly.

(A) Representative FECs obtained in the presence of 60 nM SNAP-25B and 2 µM WT Munc18-1 (#1–5 ) or Munc18-1 mutants P335A (#6) or T248G (#7). FECs #1–5 represent consecutive rounds of manipulation …

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

MATLAB figure corresponding to Figure 7A (FECs #1–5).

https://doi.org/10.7554/eLife.41771.033
Figure 7—source data 2

MATLAB figure corresponding to Figure 7B.

https://doi.org/10.7554/eLife.41771.034
Figure 7—figure supplement 1
FECs obtained in the presence of 2 µM Munc18-1 and 60 nM SNAP-25B in the solution.

Red arrows mark events of SNAP-25B binding and SNARE assembly. Note that syntaxin +2 layer mutation I233G/E234G/Y235G significantly weakens CTD zippering.

https://doi.org/10.7554/eLife.41771.031
Figure 7—figure supplement 2
Extension-time trajectories at constant mean forces.

The trajectories exhibit reversible folding and unfolding transitions of the mutant template complexes and irreversible SNAP-25 binding (indicated by red arrows).

https://doi.org/10.7554/eLife.41771.032
Munc18-1 does not significantly accelerate zippering between t- and v-SNAREs.

(A) FECs obtained by pulling single WT SNARE complexes in 2 µM soluble VAMP2 in the absence or presence of 2 µM Munc18-1. Magenta arrows mark binding of the VAMP2 molecules in the solution to the …

https://doi.org/10.7554/eLife.41771.035
Munc18-1 phosphomimetic and disease mutations altered chaperoned SNARE assembly.

(A) FECs for Munc18-1 mutations with 0 nM (#3) or 60 nM (others) SNAP-25B. See also Figure 9—source data 1. (B) Extension-time trajectories at the indicated constant mean forces, some of which (e …

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

MATLAB figure corresponding to Figure 9A (FECs #3–7).

https://doi.org/10.7554/eLife.41771.040
Figure 9—source data 2

MATLAB figure corresponding to Figure 9B.

https://doi.org/10.7554/eLife.41771.041
Figure 10 with 3 supplements
Munc18-3 and Vps33 catalyze SNARE assembly via template complexes.

(A) FECs of the Munc18-3 or Vps33 cognate Qa-R SNARE conjugate in the presence of the indicated protein(s). See also Figure 10—source datas 1 and 2. (B–C) Extension-time trajectories at the …

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

MATLAB figure corresponding to Figure 10A (FECs #1–3).

https://doi.org/10.7554/eLife.41771.046
Figure 10—source data 2

MATLAB figure corresponding to Figure 10A (FECs #4–7).

https://doi.org/10.7554/eLife.41771.047
Figure 10—source data 3

MATLAB figure corresponding to Figure 10B (traces a-e).

https://doi.org/10.7554/eLife.41771.048
Figure 10—source data 4

MATLAB figure corresponding to Figure 10C (traces g-i).

https://doi.org/10.7554/eLife.41771.049
Figure 10—figure supplement 1
FECs obtained by pulling and relaxing a single syntaxin-4-VAMP2 conjugate (#1–4) or Vam3-Nyv1 conjugate (#5–8) in the presence of the indicated protein or proteins.

SNARE CTD transitions and template complex transitions are marked by gray and blue ovals, respectively. Gray arrows indicate SNARE unzipping. Vps33(Δ354–376) is analogous to Munc18-1 Δ324–339 and is …

https://doi.org/10.7554/eLife.41771.043
Figure 10—figure supplement 2
FECs displaying Vps33-catalyzed vacuolar SNARE assembly, marked by red arrows.

Gray arrows indicate SNARE unzipping. Template complex transitions are marked by blue ovals.

https://doi.org/10.7554/eLife.41771.044
Figure 10—figure supplement 3
Probabilities of vacuolar SNARE assembly per relaxation under different conditions.

The insert shows the pulling direction and the region of the Vps33 truncation (yellow).

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

Videos

Video 1
SNARE complex unfolding and subsequent template complex formation as inferred from single-molecule measurements.

The proposed state transitions associated with FEC #2 in Figure 2C or Figure 2—figure supplement 2 are simulated.

https://doi.org/10.7554/eLife.41771.013
Video 2
Template complex facilitates SNAP-25B binding and SNARE assembly.

The extension at a constant mean force of 6.0 pN corresponding to trace b in Figure 7B and its associated state transition are simulated. For simplicity, only the right bead was simulated to move in …

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

Tables

Table 1
Properties of the neuronal template complex.
https://doi.org/10.7554/eLife.41771.014
SNARE
or SM
Mutation
or truncation
Unfolding
energy
(kBT)
Equilibrium
force* (pN)
Folding
rate
(s−1)
Unfolding
rate (s−1)
Partially
closed
syntaxin
Template
formation
SNAP-25
binding
Prob.Prob.N§Prob.N**
WT-5.2 (0.1)5.1 (0.1)1320.70.40.53460.750
Munc18-1L247R1.6 (0.3)2.3 (0.1)--0.30.3990.76
T248G2.9 (0.2)3.1 (0.1)--00.31550.316
L247A/
T248G
<1.5***---00241--
S306D¶¶5.8 (0.1)5.6 (0.1)1840.60.40.91230.953
L307R4.1 (0.2)4.6 (0.1)0.070.431140.5819
S313D¶¶6.1 (0.2)5.7 (0.1)5681.50.411620.870
Δ324–
339††,‡‡
<1.5***--0010500
D326K¶¶6.5 (0.2)5.7 (0.1)4200.60.030.9103127
P335A§§6.0 (0.3)5.9 (0.1)2580.50.020.71550.911
P335L§§4.3 (0.1)4.8 (0.1)170.20.40.32240.836
L341P§§<1.5***--0.060.041760.54
L348R††,‡‡<1.5***--0.020.042220.76
Y473D‡‡4.0 (0.1)4.3 (0.2)--00.13950.524
VAMP2L32G/Q33G3.4 (0.2)3.9 (0.1)310100.40.61700.0633
V39D3.8 (0.4)3.9 (0.2)9020.30.11750.813
M46A5.2 (0.4)5.1 (0.2)1300.70.30.5520.813
E62T††4.1 (0.2)4.8 (0.2)10750.40.51040.423
S61D/
E62T††
3.6 (0.2)4.1 (0.1)0.40.7560.212
Q76A††4.7 (0.2)4.8 (0.1)16620.40.6620.312
F77A‡‡1.5 (0.3)2.3--0.50.11210.56
A81G/A82G5.0 (0.3)4.9 (0.2)1300.80.40.51490.442
Δ85–945.1 (0.2)5.0 (0.1)1200.70.40.5870.729
Syntaxin-1ΔNRD††,‡‡<1.5***---00.081050.212
ΔN-
peptide††,‡‡
3.2 (0.2)4.6 (0.1)4220.030.53280.446
ΔHabc‡‡<1.5***---00.061400.54
L165A/E166A
(LE)¶¶
6.7 (0.2)6.1 (0.1)4060.50.070.7830.926
LE/E76K6.4 (0.2)6.0 (0.2)1230.20.070.9810.730
I202G/I203G3.0 (0.3)3.8 (0.1)240120.40.51770.433
F216A3.7 (0.1)5.1 (0.1)82200.61550.932
I230G/D231G/
R232G†††
3.6 (0.2)4.3 (0.1)--00.51110.47
I233G/E234G/
Y235G†††
3.0 (0.2)4.1 (0.1)--00.61220.730
V237G/E238G/
H239G
5.2 (0.2)4.9 (0.1)1240.70.010.31820.414
T251G/K252G5.2 (0.1)4.9 (0.1)1260.70.50.81970.747
Δ255–2645.4 (0.2)5.1 (0.1)1400.60.50.51340.729
Syntaxin-1L165A/E166A6.6 (0.2)6.2 (0.1)720.10.20.9850.211
Munc18-1D326K¶¶
  1. * Mean of two average forces for the unfolded and folded states when the two states are equally populated (Rebane et al., 2016). The equilibrium force of the template complex generally correlates with its unfolding energy. The number in parentheses is the standard error of the mean.

    † Detected as the syntaxin- and Munc18-1-dependent transition in the force range of 10–15 pN.

  2. Probability per relaxation or pulling measured in the absence of SNAP-25B.

    § Total number of pulling or relaxation FECs acquired, in which transitions of the template complex or syntaxin are scored, including their average equilibrium forces and extension changes.

  3. Probability of SNAP-25B binding and SNARE assembly per relaxation upon formation of the template complex.

    ** Total number of relaxation FECs containing the template complex transition.

  4. †† Mutation that reduces membrane fusion in vitro (Parisotto et al., 2014; Shen et al., 2010; Shen et al., 2007).

    ‡‡ Mutation that impairs exocytosis or neurotransmitter release in vivo (Meijer et al., 2018; Munch et al., 2016; Walter et al., 2010).

  5. §§ Mutation associated with epilepsy (Stamberger et al., 2016).

    ¶¶ Mutation that enhances membrane fusion in vitro or neurotransmitter release in the cell (Genc et al., 2014; Gerber et al., 2008; Lai et al., 2017; Munch et al., 2016; Parisotto et al., 2014; Richmond et al., 2001).

  6. *** Unfolding energy below the detection limit of our method, estimated to be 1.5 kBT, or not available due to no, infrequent, or heterogeneous template complex transition.

    ††† In the observed template complex transition, the template complex frequently dwelled in the unfolded state for an unusually long time (Figure 3—figure supplement 1). Thus, the transition is no longer two-state.

Key resources table
Reagent type
(species)
or resource
DesignationSource or
reference
IdentifiersAdditional
information
Strain, strain
background
(species)
BL21 Gold (DE3)
competent cells
Agilent
echnologies
Cat#230132
Commercial
assay or kit
BirA-500: BirA
biotin-protein
ligase standard
reaction kit
AvidityCat#BirA500
Chemical
compound, drug
dNTP mix (10 mM)InvitrogenCat#18427013
Chemical
compound, drug
2,2'-dithiodipyridine
disulfide (DTDP)
Sigma-AldrichCAS#2127-03-9
Chemical
compound, drug
Protease inhibitor
cocktail tablet,
EDTA free
RocheCat#05892791001
Peptide,
recombinant protein
Catalase from
bovine liver
Sigma-AldrichCAS#9001-05-2
Peptide,
recombinant protein
Glucose Oxidase
from Aspergillus
niger
Sigma-AldrichCAS#9001-37-0
Software,
algorithm
LabVIEW VIs(Gao et al., 2012)instrument control,
data acquisition,
and data analysis
Software,
algorithm
MATLAB scripts(Gao et al., 2012)
(Rebane et al., 2016)
data analysis
Software,
algorithm
GeneiousGeneiousDNA primer design
Software,
algorithm
GraphPad Prism7GraphPad
Software
OtherMicro Bio-
spin six columns
Bio-RADCat#732–6221
OtherNi Sepharose
6 Fast Flow
GE healthcare
Lifesciences
Cat#17531801
OtherAnti-digoxigenin
antibody coated
polystyrene particles
SpherotechCat#DIGP-20–22.1 µm, called
DIG beads
OtherStreptavidin-coated
polystyrene particles
SpherotechCat#SVP-15–51.8 µm
OtherCustomized glass
tubing: bead dispenser
tubes with 100 µm
outer diameter (OD) and
25 µm inner diameter (ID)
King Precision
Glass, Inc
OtherPolyethylene
tubing PE50
Becton
Dickinson
Cat# 22–270835
OtherDual optical
trap setup
(Gao et al., 2012)

Additional files

Download links