A lever hypothesis for Synaptotagmin-1 action in neurotransmitter release
Data files
Aug 10, 2024 version files 31.67 GB
-
dryadMDSyt1action1.tar.gz
-
dryadMDSyt1action1b.tar
-
README.md
Abstract
Neurotransmitter release is triggered in microseconds by Ca2+-binding to the Synaptotagmin-1 C2 domains and by SNARE complexes that form four-helix bundles between synaptic vesicles and plasma membranes, but the coupling mechanism between Ca2+-sensing and membrane fusion is unknown. Release requires extension of SNARE helices into juxtamembrane linkers that precede transmembrane regions (linker zippering) and binding of the Synaptotagmin-1 C2B domain to SNARE complexes through a ‘primary interface’ comprising two regions (I and II). The Synaptotagmin-1 Ca2+-binding loops were believed to accelerate membrane fusion by inducing membrane curvature, perturbing lipid bilayers or helping bridge the membranes, but SNARE complex binding orients the Ca2+-binding loops away from the fusion site, hindering these putative activities. Molecular dynamics simulations now suggest that Synaptotagmin-1 C2 domains near the site of fusion hinder SNARE action, providing an explanation for this paradox and arguing against previous models of Sytnaptotagmin-1 action. NMR experiments reveal that binding of C2B domain arginines to SNARE acidic residues at region II remains after disruption of region I. These results and fluorescence resonance energy transfer assays, together with previous data, suggest that Ca2+ causes reorientation of the C2B domain on the membrane and dissociation from the SNAREs at region I but not region II. Based on these results and electrophysiological data described in the accompanying paper, we propose that upon Ca2+ binding Synaptotagmin-1 reorients and acts as a lever that pulls the SNARE complex, facilitating linker zippering or other SNARE structural changes that lead to fast membrane fusion.synaptotagmin, SNAREs, neurotransmitter release, membrane fusion, molecular dynamics simulation, synaptic vesicle fusion
README: A lever hypothesis for Synaptotagmin-1 action in neurotransmitter release
https://doi.org/10.5061/dryad.mgqnk9978
This deposition contains files corresponding to the molecular dynamics (MD) simulations, bimane fluorescence quenching assays, NMR spectra and fluorescence resonance energy transfer experiments described in the reference paper.
Description of the data and file structure
The deposition has two tar files:
- dryadMDSyt1action1.tar.gz contains the MD simulation data
- dryadMDSyt1action1b.tar contains the bimane, NMR and FRET data
To extract the files, use the following commands in Linux:
gunzip dryadMDSyt1action1.tar.gz
tar xvf dryadMDSyt1action1.tar
tar xvf dryadMDSyt1action1b.tar
Content of dryadMDSyt1action1.tar.gz
These data sets contain files associated with all-atom molecular dynamics simulations of the neurotransmitter release machinery between a flat bilayer and a vesicle. The simulations were performed with gromacs. Each directory corresponds to one simulation as described briefly below and in more detail in the reference associated with this deposition (see the results part, figure legends and Table S1 of the reference).
cac2absc: four trans-SNARE complexes zippered to distinct extents at the C-terminus, bridging a vesicle and a flat bilayer, and including four synaptotagmin-1 C2AB molecules bound to five Ca2+ ions each and interspersed between the SNARE complexes to investigate whether the synaptotagmin-1 C2 domains spontaneously insert into the bilayers, perturb or bridge the bilayers, and perhaps even initiate membrane fusion
fusiong: similar systems as cac2absc but with all SNARE four-helix bundles almost fully assembled and with the C2B domain Ca2+-binding loops oriented to the center of the membrane-membrane interface, to further explore the possibility that the C2 domains might cooperate with the SNAREs and increase the probability of observing fusion.
nosytfusion: control system analogous to fusiong but without C2AB molecules to examine the possibility that the C2AB molecules might actually hinder SNARE action in this system.
sytfusion2: similar system as fusiong but with two C2AB molecules instead of four; the C2 domains were located between the two membranes with the Ca2+-binding loops oriented toward the center of the membrane-membrane interface such that one of the loops can bind to the vesicle and the other to the flat bilayer, to explore the possibility that Ca2+-binding might favor movement of lipids toward the Ca2+-binding sites to destabilize the bilayers and initiate fusion.
sytfusion3: similar system as fusiong but with the four trans-SNARE complexes closer to the center of the membrane-membrane interface and with four C2AB molecules in positions in which the two C2 domains were poised to bridge the two membranes, to investigate whether , the C2 domains might act as wedges that prevent the membranes from coming closer while the SNARE complexes pull the membranes together in the center, resulting in torque forces that help to bend the membranes to initiate fusion.
s1action2: four almost fully assembled SNAREs complexes bridging a vesicle and a flat bilayer, and bound to one C2AB molecule each through the primary interface to investigate the interactions between C2B domain arginines and SNARE acidic residues at this interface.
For each simulation, filenames have a common beginning (e.g. cac2absc for the simulation in the cac2absc directory) followed by a few letters that indicate the nature or order of the simulation: nvt for temperature equilibration runs; par for pressure equilibration runs and mdx for production simulations, where x is a number that denotes the order of a particular run within the chain of concatenated simulations. The filename extensions follow the common gromacs nomenclature. Here is a list of the types of files included in the deposition.
- Filenamea.sh or Filenameb.sh: text files with the SLURM instructions to run the simulation. These files show which files were used to run the particular simulation.\
- Filename.top: topology file that includes calls to the force field and to other topology files and/or restraint files\
- Filename.itp: topology or restraint file\
- Filename.ndx: index file\
- Filename.mdp: file containing the molecular dynamics parameters\
- Filename.log: log file for the simulation\
- Filename.gro: file with the coordinates and velocities at the end of the simulation\
- Filename.cpt: checkpoint file containing the final coordinates and velocities with high precision to continue the simulation\
- Filaname.edr: energy file
- Filename.pdb: coordinate files used to generate the initial system configuration.
In addition, each simulation directory contains force field, restraint and/or topology directories. The filename.xtc trajectory files that contain coordinates and velocities of frames saved during the simulations are not included because of their large size, which would render the deposition very expensive. These files are available from Josep Rizo upon reasonable request. This is the contact information:
Name: Josep Rizo\
Institution: Department of Biophysics, UT Southwestern Medical Center Dallas\
Address: 6001 Forest Park Road, Dallas, TX 75390-8816\
Email: Jose.Rizo-Rey@UTSouhtwestern.edu
Content of dryadMDSyt1action1b.tar
This tar file includes the following files:
Fig3Bimanedata.xlsx: Excel file with the bimane fluorescence data of Fig. 3 of the reference paper, which shows that Ca2+ does not induce dissociation of Synaptotagmin1 C2AB from membrane-anchored SNARE complex. The graph in the figure shows fluorescence spectra of nanodisc-anchored SNARE complex-59-bimane alone, or in the presence of C2AB-T285W and Mg-ATP plus 1 mM EGTA or 1 mM Ca2+. Mg-ATP was included in the samples to hinder non-specific C2AB-SNARE complex interactions.
Fig6FRETdata.xlsx: Excel file with the FRET data of Fig. 6 of the reference paper, including the replicates used for the statistics of Fig. 6e. These data show that Ca2+ induces reorientation of the Syt1 C2B domain with respect to membrane-anchored SNARE complex. The data corresponding to panels (B-D) and their replicates are fluorescence spectra of liposome-anchored C2AB-SNARE-complexin-1 complex with stabilized Cy3 and Cy5 probes attached to the following residues of SNAP-25 and the C2B domain, respectively: 27-412 (B), 34-412 (C) and 34-391 (D). The spectra were obtained in the presence of Mg-ATP plus 1 mM EGTA , 1 mM Ca2+ or 1% b-OG plus NSF and aSNAP as a control to disrupt the Syt1-SNARE complex. Mg-ATP was included in the samples to hinder non-specific C2AB-SNARE complex interactions. The data corresponding to panel (E) show the quantification of the FRET difference between the spectra obtained in 1 mM Ca2+ and 1 mM EGTA normalized with the FRET observed in the presence of 1 mM EGTA. Values from individual experiments as well as mean values and standard deviations are included.
Fig4andS9toS14 and FigS6toS8andS15: Directories including the NMR data sets corresponding to Figs. 4 and S6-S15 of the reference paper. Within each of these directories there is one directory for each titration included in the corresponding figure and within each titration directory there is one directory for each spectrum acquired during the titration. These directories have extension .fid and their name specifies the synaptotagmin-1 C2B domain version uses for the experiments [wild type (WT) or mutant]. Each of these directories includes files that contain the raw NMR data in Agilent format (fid), the acquisition parameters (procpar), a log file (log), a text file with a description of experimental conditions (text), macros to convert from Agilent format to NMRPipe format (fid.com) and to process the data with NMRPipe (jrtrgood2c.com), an intermediate processing file (tr.ft), the processed NMR spectrum in NMRPipe format (tr.ft2) and the same spectrum converted to NMR view format (same name as the directory but with extension .nv). For each titration, there is also a directory called nvs that contains all the .nv files. To increase the signal to noise ratio, some of the experiments were acquired multiple times and added. For added experiments, an additional directory was created (directory name includes _add) and those nv files were used for further analysis. The spectra included in the addition are indicated by the numbers at the end of the filename. Processed files with extension .nv or .ft2 are ready to be viewed in NMRViewJ or NMRDraw, respectively.
In case of any questions, please do not hesitate to contact Josep Rizo.
Methods
MD simulations. All-atom MD simulations were performed using Gromacs with the CHARMM36 force field. Solvation and ion addition for system setup were performed at the BioHPC supercomputing facility of UT Southwestern. Minimizations, equilibration steps and production molecular dynamics (MD) simulations were carried out on Frontera at the Texas Advanced Computing Center (TACC). Pymol (Schrödinger, LLC) was used for system design, manual manipulation and system visualization. The methodology used to set up the systems and run MD simulations was analogous to that described previously (1, 2). Systems were energy minimized using double precision, whereas the default mixed precision was used in all MD simulations. The systems were heated to the desired temperature running a 1 ns simulation in the NVT ensemble with 1 fs steps, and then equilibrated to 1 atm for 1 ns in the NPT ensemble with isotropic Parrinello-Rahman pressure coupling and 2 fs steps. NPT production MD simulations were performed for the times indicated in Table S1 using 2 fs steps, isotropic Parrinello-Rahman pressure coupling and a 1.1 nm cutoff for non-bonding interactions. Three different groups of atoms were used for Nose-Hoover temperature coupling: i) protein atoms; ii) lipid atoms; and iii) water and KCL ions. Periodic boundary conditions were imposed with Particle Mesh Ewald (PME) summation for long-range electrostatics.
NMR spectroscopy. All NMR spectra were acquired at 25 °C on Agilent DD2 spectrometers equipped with cold probes operating at 600 or 800 MHz. Titrations of WT and mutant 2H,15N-labeled Syt1 C2B domain specifically 13CH3-labeled at the Ile d1 and Met methyl groups (referred to as 15N-C2B for simplicity) with SNARE complex four-helix bundle bound to a fragment spanning residues 26-83 of complexin-1 (referred to as CpxSC) were performed as described in (3). Specific 13CH3 labeling was performed for acquisition of 1H-13C heteronuclear multiple quantum (HMQC) spectra with these samples, although no such spectra are described here. Because the R322E/K325E mutation increases the stability of the C2B domain and slows down H/D exchange, some amide groups in the b-strands are not fully protonated after expression in D2O and purification in buffers containing H2O, resulting in signal loss for the corresponding peaks. All newly purified C2B mutants bearing the R322E/K325E mutation were incubated at RT for one week or at 37 °C for 15 hours to facilitate full exchange, but some amide groups were still not fully exchanged after this procedure. The titrations were performed in NMR buffer containing 20 mM HEPES (pH 7.4), 100 mM KCl, 1 mM EDTA, 1 mM TCEP, 10% D2O and protease inhibitor cocktail [which contained 1 mM Antipain Dihydrochloride (Thermo Fischer Scientific: 50488492), 20 mM Leupeptin (Gold Bio: L01025) and 0.8 mM Aprotinin (Gold Bio: A655100)]. A 1H-15N TROSY-HSQC spectrum was acquired first for 32 μM isolated 15N-C2B domain and additional 1H-15N TROSY-HSQC spectra were acquired after adding increasing concentrations of CpxSC to the sample, resulting in gradual dilution of the 15N-C2B domain. The protein concentrations of each titration step for each mutant are indicated in the legends of Fig. 4 and S15. Soluble SNARE complex was assembled as described in (3), concentrated at room temperature to a concentration above 250 μM using a 30 kDa centrifugation filter (Amicon) and exchanged into NMR buffer using Zeba Spin Desalting Columns, 7K MWCO, 10 mL (Thermo Fisher). Complexin-1 (26–83) was also concentrated above 250 μM using a 3 kDa centrifugation filter (Amicon) and exchanged into NMR buffer. SNARE-Complexin-1 (26–83) complex was preassembled with 20% excess Complexin-1 (26–83) before mixing with 15N-C2B domain. Total acquisition times ranged from 3.5 to 87.5 hr, depending on the sensitivity of the spectra. All NMR data were processed with NMRPipe and analyzed with NMRViewJ.
Bimane fluorescence quenching assay. All fluorescence emission scans were collected on a PTI Quantamaster 400 spectrofluorometer (T-format) at room temperature with slits set to 1.25 mm. For tryptophan-induced bimane fluorescence quenching assays, we used SNARE complexes anchored on nanodiscs as described (3) and labeled with bimane at position R59C of SNAP-25_N. The lipid composition of nanodiscs was 84% POPC, 15% DOPS, 1% PIP2. The experiments were performed in 25 mM HEPES pH 7.4 100 mM KCl 0.1 mM TCEP 2.5 mM MgCl2, 2 mM ATP 1 mM EGTA containing 1.5 μM BSA to prevent sample binding to the cuvette. Fluorescence emission spectra (excitation at 380 nm) were acquired for samples containing 1 mM SC-59-bimane-nanodiscs alone or in the presence of 4 mM C2A-T285W without or with 2.0 mM CaCl2 (1.0 mM free Ca2+).
FRET assays. All fluorescence emission scans were collected on a PTI Quantamaster 400 spectrofluorometer (T-format) at room temperature with excitation at 550 nm and slits set to 1.25 mm. Each sample contained 0.125 μM C2AB-SNARE-complex-liposomes in 25 mM HEPES pH 7.4 100 mM KCl buffer, 1 mM EGTA, 2.5 mM MgCl2, 2 mM ATP, 1.5 μM BSA, and 0.15 μM full-length Complexin-1. Mg-ATP and complexin-1 were included to hinder non-specific interactions. To examine the reproducibility of Ca2+-induced changes in FRET, spectra were acquired for separate samples prepared under identical conditions before and after addition of 2.0 mM CaCl2 (1.0 mM free Ca2+). Because the Syt1 C2AB fragment and SNARE complex were fused covalently via the 37 aa linker and the SNARE complex is resistant to SDS, control spectra to measure the maximum donor fluorescence observable without FRET were acquired after incubation for 5 minutes at 37 oC with 0.4 µM NSF and 2 µM aSNAP (to disassemble the SNARE complex) in 25 mM HEPES pH 7.4 100 mM KCl buffer containing 1 mM EGTA, 2.5 mM MgCl2, 2 mM ATP, 1.5 μM BSA and 1% BOG. The presence of detergent proved to be necessary to recover maximum signal because a subset of the C2AB-SNARE-complexes get reconstituted inside the liposomes and hence are inaccessible to the disassembly machinery.
References
(1) J. Rizo, L. Sari, Y. Qi, W. Im, M. M. Lin, All-atom molecular dynamics simulations of Synaptotagmin-SNARE-complexin complexes bridging a vesicle and a flat lipid bilayer. Elife 11, e76356 (2022).
(2) J. Rizo, L. Sari, K. Jaczynska, C. Rosenmund, M. M. Lin, Molecular mechanism underlying SNARE-mediated membrane fusion enlightened by all-atom molecular dynamics simulations. Proc Natl Acad Sci U S A 121, e2321447121 (2024).
(3) R. Voleti, K. Jaczynska, J. Rizo, Ca(2+)-dependent release of Synaptotagmin-1 from the SNARE complex on phosphatidylinositol 4,5-bisphosphate-containing membranes. Elife 9, e57154 (2020).