All-atom molecular dynamics simulations of Synaptotagmin-SNARE-complexin complexes bridging a vesicle and a flat lipid bilayer
Abstract
Synaptic vesicles are primed into a state that is ready for fast neurotransmitter release upon Ca2+-binding to Synaptotagmin-1. This state likely includes trans-SNARE complexes between the vesicle and plasma membranes that are bound to Synaptotagmin-1 and complexins. However, the nature of this state and the steps leading to membrane fusion are unclear, in part because of the difficulty of studying this dynamic process experimentally. To shed light into these questions, we performed all-atom molecular dynamics simulations of systems containing trans-SNARE complexes between two flat bilayers or a vesicle and a flat bilayer with or without fragments of Synaptotagmin-1 and/or complexin-1. Our results need to be interpreted with caution because of the limited simulation times and the absence of key components, but suggest mechanistic features that may control release and help visualize potential states of the primed Synaptotagmin-1-SNARE-complexin-1 complex. The simulations suggest that SNAREs alone induce formation of extended membrane-membrane contact interfaces that may fuse slowly, and that the primed state contains macromolecular assemblies of trans-SNARE complexes bound to the Synaptotagmin-1 C2B domain and complexin-1 in a spring-loaded configuration that prevents premature membrane merger and formation of extended interfaces, but keeps the system ready for fast fusion upon Ca2+ influx.
Data availability
Most files corresponding to our molecular dynamics simulations are available in the dryad database (doi:10.5061/dryad.ns1rn8pw6). Because of the very large size of trajectory files, it was not practical to deposit them in this database, but these files are available from the corresponding author upon reasonable request.
-
Data from: All-atom molecular dynamics simulations of Synaptotagmin-SNARE-complexin complexes bridging a vesicle and a flat lipid bilayerDryad Digital Repository, doi:10.5061/dryad.ns1rn8pw6.
Article and author information
Author details
Funding
National Institute of Neurological Disorders and Stroke (R35 NS097333)
- Josep Rizo
Welch Foundation (I-1304)
- Josep Rizo
National Science Foundation (MCB-2111728)
- Wonpil Im
Natural Science Foundation of Shanghai Grant (19ZR1473600)
- Yife Qi
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2022, Rizo et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 3,262
- views
-
- 993
- downloads
-
- 38
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Structural Biology and Molecular Biophysics
Agonists enhance receptor activity by providing net-favorable binding energy to active over resting conformations, with efficiency (η) linking binding energy to gating. Previously, we showed that in nicotinic receptors, η-values are grouped into five structural pairs, correlating efficacy and affinity within each class, uniting binding with allosteric activation (Indurthi and Auerbach, 2023). Here, we use molecular dynamics (MD) simulations to investigate the low-to-high affinity transition (L→H) at the Torpedo α−δ nicotinic acetylcholine receptor neurotransmitter site. Using four agonists spanning three η-classes, the simulations reveal the structural basis of the L→H transition where: the agonist pivots around its cationic center (‘flip’), loop C undergoes staged downward displacement (‘flop’), and a compact, stable high-affinity pocket forms (‘fix’). The η derived from binding energies calculated in silico matched exact values measured experimentally in vitro. Intermediate states of the orthosteric site during receptor activation are apparent only in simulations, but could potentially be observed experimentally via time-resolved structural studies.
-
- Structural Biology and Molecular Biophysics
We designed novel pre-drug compounds that transform into an active form that covalently modifies particular His residue in the active site, a difficult task to achieve, and applied to carbonic anhydrase (CAIX), a transmembrane protein, highly overexpressed in hypoxic solid tumors, important for cancer cell survival and proliferation because it acidifies tumor microenvironment helping invasion and metastases processes. The designed compounds have several functionalities: (1) primary sulfonamide group recognizing carbonic anhydrases (CA), (2) high-affinity moieties specifically recognizing CAIX among all CA isozymes, and (3) forming a covalent bond with the His64 residue. Such targeted covalent compounds possess both high initial affinity and selectivity for the disease target protein followed by complete irreversible inactivation of the protein via covalent modification. Our designed prodrug candidates bearing moderately active pre-vinylsulfone esters or weakly active carbamates optimized for mild covalent modification activity to avoid toxic non-specific modifications and selectively target CAIX. The lead inhibitors reached 2 pM affinity, the highest among known CAIX inhibitors. The strategy could be used for any disease drug target protein bearing a His residue in the vicinity of the active site.