All-atom molecular dynamics simulations of Synaptotagmin-SNARE-complexin complexes bridging a vesicle and a flat lipid bilayer

Non-ionotropic voltage-gated calcium channel signaling

Voltage-gated calcium channels (VGCCs) are the major conduits for calcium ions (Ca2+) within excitable cells. Recent studies have highlighted the non-ionotropic functionality of VGCCs, revealing their capacity to activate intracellular pathways independently of ion flow. This non-ionotropic signaling mode plays a pivotal role in excitation-coupling processes, including gene transcription through excitation-transcription (ET), synaptic transmission via excitation-secretion (ES), and cardiac contraction through excitation-contraction (EC). However, it is noteworthy that these excitation-coupling processes require extracellular calcium (Ca2+) and Ca2+ occupancy of the channel ion pore. Analogous to the "non-canonical" characterization of the non-ionotropic signaling exhibited by the N-methyl-D-aspartate receptor (NMDA), which requires extracellular Ca2+ without the influx of ions, VGCC activation requires depolarization-triggered conformational change(s) concomitant with Ca2+ binding to the open channel. Here, we discuss the contributions of VGCCs to ES, ET, and EC coupling as Ca2+ binding macromolecules that transduces external stimuli to intracellular input prior to elevating intracellular Ca2+. We emphasize the recognition of calcium ion occupancy within the open ion-pore and its contribution to the excitation coupling processes that precede the influx of calcium. The non-ionotropic activation of VGCCs, triggered by the upstroke of an action potential, provides a conceptual framework to elucidate the mechanistic aspects underlying the microseconds nature of synaptic transmission, cardiac contractility, and the rapid induction of first-wave genes.

Making the cut: Multiscale simulation of membrane remodeling

Biological membranes are dynamic heterogeneous materials, and their shape and organization are tightly coupled to the properties of the proteins in and around them. However, the length scales of lipid and protein dynamics are far below the size of membrane-bound organelles, much less an entire cell. Therefore, multiscale modeling approaches are often necessary to build a comprehensive picture of the interplay of these factors, and have provided critical insights into our understanding of membrane dynamics. Here, we review computational methods for studying membrane remodeling, as well as passive and active examples of protein-driven membrane remodeling. As the field advances towards the modeling of key aspects of organelles and whole cells - an increasingly accessible regime of study - we summarize here recent successes and offer comments on future trends.

Biomolecular simulations at the exascale: From drug design to organelles and beyond

The rapid advancement in computational power available for research offers to bring not only quantitative improvements, but also qualitative changes in the field of biomolecular simulation. Here, we review the state of biomolecular dynamics simulations at the threshold to exascale resources becoming available. Both developments in parallel and distributed computing will be discussed, providing a perspective on the state of the art of both. A main focus will be on obtaining binding and conformational free energies, with an outlook to macromolecular complexes and (sub)cellular assemblies.

Neurotransmitter release is triggered by a calcium-induced rearrangement in the Synaptotagmin-1/SNARE complex primary interface

The Ca2+ sensor synaptotagmin-1 triggers neurotransmitter release together with the neuronal SNARE complex formed by syntaxin-1, SNAP25 and synaptobrevin. Moreover, synaptotagmin-1 increases synaptic vesicle priming and impairs spontaneous vesicle release. The synaptotagmin-1 C2B domain binds to the SNARE complex through a primary interface via two regions (I and II), but how exactly this interface mediates distinct functions of synaptotagmin-1, and the mechanism underlying Ca2+-triggering of release is unknown. Using mutagenesis and electrophysiological experiments, we show that region II is functionally and spatially subdivided: binding of C2B domain arginines to SNAP-25 acidic residues at one face of region II is crucial for Ca2+-evoked release but not for vesicle priming or clamping of spontaneous release, whereas other SNAP-25 and syntaxin-1 acidic residues at the other face mediate priming and clamping of spontaneous release but not evoked release. Mutations that disrupt region I impair the priming and clamping functions of synaptotagmin-1 while, strikingly, mutations that enhance binding through this region increase vesicle priming and clamping of spontaneous release, but strongly inhibit evoked release and vesicle fusogenicity. These results support previous findings that the primary interface mediates the functions of synaptotagmin-1 in vesicle priming and clamping of spontaneous release, and, importantly, show that Ca2+-triggering of release requires a rearrangement of the primary interface involving dissociation of region I, while region II remains bound. Together with modeling and biophysical studies presented in the accompanying paper, our data suggest a model whereby this rearrangement pulls the SNARE complex to facilitate fast synaptic vesicle fusion.

A lever hypothesis for Synaptotagmin-1 action in neurotransmitter release

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 molecular modeling, we propose that Synaptotagmin-1 acts as a lever that pulls the SNARE complex when Ca2+ causes reorientation of the C2B domain, facilitating linker zippering and fast membrane fusion. This hypothesis is supported by the electrophysiological data described in the accompanying paper.

Coarse-grained modeling of annexin A2-induced microdomain formation on a vesicle

Annexin A2 (A2)-induced microdomain formation is a key step in biological processes such as Ca2+-mediated exocytosis in neuroendocrine cells. In this work, a total of 15 coarse-grained molecular dynamics simulations were performed on vesicle models having a diameter of approximately 250 Å for 15 μs each using the Martini2 force field. Five simulations were performed in the presence of 10 A2, 5 in the presence of A2 but absence of PIP2, and 5 simulations in the absence of A2 but presence of PIP2. Consistent results were generated among the simulations. A2-induced PIP2 microdomain formation was observed and shown to occur in three phases: A2-vesicle association, localized A2-induced PIP2 clustering, and A2 aggregation driving PIP2 microdomain formation. The relationship between A2 aggregation and PIP2 microdomain formation was quantitatively described using a novel method which calculated the variance among protein and lipid positions via the Fréchet mean. A large reduction in PIP2 variance was observed in the presence of A2 but not in its absence. This reduction in PIP2 variance was proportional to the reduction observed in A2 variance and demonstrates that the observed PIP2 microdomain formation is dependent upon A2 aggregation. The three-phase model of A2-induced microdomain formation generated in this work will serve as a valuable guide for further experimental studies and the development of novel A2 inhibitors. No microdomain formation was observed in the absence of A2 and minimal A2-membrane interaction was observed in the absence of PIP2.

Vesicle and reaction-diffusion hybrid modeling with STEPS

Vesicles carry out many essential functions within cells through the processes of endocytosis, exocytosis, and passive and active transport. This includes transporting and delivering molecules between different parts of the cell, and storing and releasing neurotransmitters in neurons. To date, computational simulation of these key biological players has been rather limited and has not advanced at the same pace as other aspects of cell modeling, restricting the realism of computational models. We describe a general vesicle modeling tool that has been designed for wide application to a variety of cell models, implemented within our software STochastic Engine for Pathway Simulation (STEPS), a stochastic reaction-diffusion simulator that supports realistic reconstructions of cell tissue in tetrahedral meshes. The implementation is validated in an extensive test suite, parallel performance is demonstrated in a realistic synaptic bouton model, and example models are visualized in a Blender extension module.

Dynamic Formation of the Protein-Lipid Pre-fusion Complex

Synaptic vesicles (SVs) fuse with the presynaptic membrane (PM) to release neuronal transmitters. The SV protein Synaptotagmin 1 (Syt1) serves as a Ca2+ sensor for evoked fusion. Syt1 is thought to trigger fusion by penetrating into PM upon Ca2+ binding, however the mechanistic detail of this process is still debated. Syt1 interacts with the SNARE complex, a coiled-coil four-helical bundle that enables the SV-PM attachment. The SNARE-associated protein Complexin (Cpx) promotes the Ca2+-dependent fusion, possibly interacting with Syt1. We employed all-atom molecular dynamics (MD) to investigate the formation of the Syt1-SNARE-Cpx complex interacting with the lipid bilayers of PM and SV. Our simulations demonstrated that the PM-Syt1-SNARE-Cpx complex can transition to a "dead-end" state, wherein Syt1 attaches tightly to PM but does not immerse into it, as opposed to a pre-fusion state, which has the tips of the Ca2+-bound C2 domains of Syt1 inserted into PM. Our simulations unraveled the sequence of Syt1 conformational transitions, including the simultaneous Syt1 docking to the SNARE-Cpx bundle and PM, followed by the Ca2+ chelation and the penetration of the tips of Syt1 domains into PM, leading to the pre-fusion state of the protein-lipid complex. Importantly, we found that the direct Syt1-Cpx interactions are required to promote these transitions. Thus, we developed the all-atom dynamic model of the conformational transitions that lead to the formation of the pre-fusion PM-Syt1-SNARE-Cpx complex. Our simulations also revealed an alternative "dead-end" state of the protein-lipid complex that can be formed if this pathway is disrupted.

Molecular mechanism underlying SNARE-mediated membrane fusion enlightened by all-atom molecular dynamics simulations

The SNAP receptor (SNARE) proteins syntaxin-1, SNAP-25, and synaptobrevin mediate neurotransmitter release by forming tight SNARE complexes that fuse synaptic vesicles with the plasma membranes in microseconds. Membrane fusion is generally explained by the action of proteins on macroscopic membrane properties such as curvature, elastic modulus, and tension, and a widespread model envisions that the SNARE motifs, juxtamembrane linkers, and C-terminal transmembrane regions of synaptobrevin and syntaxin-1 form continuous helices that act mechanically as semirigid rods, squeezing the membranes together as they assemble ("zipper") from the N to the C termini. However, the mechanism underlying fast SNARE-induced membrane fusion remains unknown. We have used all-atom molecular dynamics simulations to investigate this mechanism. Our results need to be interpreted with caution because of the limited number and length of the simulations, but they suggest a model of membrane fusion that has a natural physicochemical basis, emphasizes local molecular events over general membrane properties, and explains extensive experimental data. In this model, the central event that initiates fast (microsecond scale) membrane fusion occurs when the SNARE helices zipper into the juxtamembrane linkers which, together with the adjacent transmembrane regions, promote encounters of acyl chains from both bilayers at the polar interface. The resulting hydrophobic nucleus rapidly expands into stalk-like structures that gradually progress to form a fusion pore, aided by the SNARE transmembrane regions and without clearly discernible intermediates. The propensity of polyunsaturated lipids to participate in encounters that initiate fusion suggests that these lipids may be important for the high speed of neurotransmitter release.

Control of Munc13-1 Activity by Autoinhibitory Interactions Involving the Variable N-terminal Region

Regulation of neurotransmitter release during presynaptic plasticity underlies varied forms of information processing in the brain. Munc13s play essential roles in release via their conserved C-terminal region, which contains a MUN domain involved in SNARE complex assembly, and controls multiple presynaptic plasticity processes. Munc13s also have a variable N-terminal region, which in Munc13-1 includes a calmodulin binding (CaMb) domain involved in short-term plasticity and a C2A domain that forms an inhibitory homodimer. The C2A domain is activated by forming a heterodimer with the zinc-finger domain of αRIMs, providing a link to αRIM-dependent short- and long-term plasticity. However, it is unknown how the functions of the N- and C-terminal regions are integrated, in part because of the difficulty of purifying Munc13-1 fragments containing both regions. We describe for the first time the purification of a Munc13-1 fragment spanning its entire sequence except for a flexible region between the C2A and CaMb domains. We show that this fragment is much less active than the Munc13-1 C-terminal region in liposome fusion assays and that its activity is strongly enhanced by the RIM2α zinc-finger domain together with calmodulin. NMR experiments show that the C2A and CaMb domains bind to the MUN domain and that these interactions are relieved by the RIM2α ZF domain and calmodulin, respectively. These results suggest a model whereby Munc13-1 activity in promoting SNARE complex assembly and neurotransmitter release are inhibited by interactions of the C2A and CaMb domains with the MUN domain that are relieved by αRIMs and calmodulin.

Studying the Effects of Dissolved Noble Gases and High Hydrostatic Pressure on the Spherical DOPC Bilayer Using Molecular Dynamic Simulations

Fine-grained molecular dynamics simulations have been conducted to depict lipid objects enclosed in water and interacting with a series of noble gases dissolved in the medium. The simple point-charge (SPC) water system, featuring a boundary composed of 1,2-Dioleoyl-sn-glycero-3-phosphocholine (DOPC) molecules, maintained stability throughout the simulation under standard conditions. This allowed for the accurate modeling of the effects of hydrostatic pressure at an ambient pressure of 25 bar. The chosen pressure references the 240 m depth of seawater: the horizon frequently used by commercial divers, who comprise the primary patient population of the neurological complication of inert gas narcosis and the consequences of high-pressure neurological syndrome. To quantify and validate the neurological effects of noble gases and discriminate them from high hydrostatic pressure, we reduced the dissolved gas molar concentration to 1.5%, three times smaller than what we previously tested for the planar bilayer (3.5%). The nucleation and growth of xenon, argon and neon nanobubbles proved consistent with the data from the planar bilayer simulations. On the other hand, hyperbaric helium induces only a residual distorting effect on the liposome, with no significant condensed gas fraction observed within the hydrophobic core. The bubbles were distributed over a large volume-both in the bulk solvent and in the lipid phase-thereby causing substantial membrane distortion. This finding serves as evidence of the validity of the multisite distortion hypothesis for the neurological effect of inert gases at high pressure.

Complexin has a dual synaptic function as checkpoint protein in vesicle priming and as a promoter of vesicle fusion

The presynaptic SNARE-complex regulator complexin (Cplx) enhances the fusogenicity of primed synaptic vesicles (SVs). Consequently, Cplx deletion impairs action potential-evoked transmitter release. Conversely, though, Cplx loss enhances spontaneous and delayed asynchronous release at certain synapse types. Using electrophysiology and kinetic modeling, we show that such seemingly contradictory transmitter release phenotypes seen upon Cplx deletion can be explained by an additional of Cplx in the control of SV priming, where its ablation facilitates the generation of a "faulty" SV fusion apparatus. Supporting this notion, a sequential two-step priming scheme, featuring reduced vesicle fusogenicity and increased transition rates into the faulty primed state, reproduces all aberrations of transmitter release modes and short-term synaptic plasticity seen upon Cplx loss. Accordingly, we propose a dual presynaptic function for the SNARE-complex interactor Cplx, one as a "checkpoint" protein that guarantees the proper assembly of the fusion machinery during vesicle priming, and one in boosting vesicle fusogenicity.

Atomistic Insight into the Lipid Nanodomains of Synaptic Vesicles

Membrane curvature, once regarded as a passive consequence of membrane composition and cellular architecture, has been shown to actively modulate various properties of the cellular membrane. These changes could also lead to segregation of the constituents of the membrane, generating nanodomains with precise biological properties. Proteins often linked with neurodegeneration (e.g., tau, alpha-synuclein) exhibit an unintuitive affinity for synaptic vesicles in neurons, which are reported to lack distinct, ordered nanodomains based on their composition. In this study, all-atom molecular dynamics simulations are used to study a full-scale synaptic vesicle of realistic Gaussian curvature and its effect on the membrane dynamics and lipid nanodomain organization. Compelling indicators of nanodomain formation, from the perspective of composition, surface areas per lipid, order parameter, and domain lifetime, are identified in the vesicle membrane, which are absent in a flat bilayer of the same lipid composition. Therefore, our study supports the idea that curvature may induce phase separation in an otherwise fluid, disordered membrane.

Control of Munc13-1 Activity by Autoinhibitory Interactions Involving the Variable N-terminal Region

Regulation of neurotransmitter release during presynaptic plasticity underlies varied forms of information processing in the brain. Munc13s play essential roles in release via their conserved C-terminal region, which contains a MUN domain involved SNARE complex assembly, and control multiple presynaptic plasticity processes. Munc13s also have a variable N-terminal region, which in Munc13-1 includes a calmodulin binding (CaMb) domain involved in short-term plasticity and a C2A domain that forms an inhibitory homodimer. The C2A domain is activated by forming a heterodimer with the zinc-finger domain of αRIMs, providing a link to αRIM-dependent short- and long-term plasticity. However, it is unknown how the functions of the N- and C-terminal regions are integrated, in part because of the difficulty of purifying Munc13-1 fragments containing both regions. We describe for the first time the purification of a Munc13-1 fragment spanning its entire sequence except for a flexible region between the C2A and CaMb domains. We show that this fragment is much less active than the Munc13-1 C-terminal region in liposome fusion assays and that its activity is strongly enhanced by the RIM2α zinc-finger domain together with calmodulin. NMR experiments show that the C2A and CaMb domains bind to the MUN domain and that these interactions are relieved by the RIM2α ZF domain and calmodulin, respectively. These results suggest a model whereby Munc13-1 activity in promoting SNARE complex assembly and neurotransmitter release are inhibited by interactions of the C2A and CaMb domains with the MUN domain that are relieved by αRIMs and calmodulin.

Quantum Computing and the Future of Neurodegeneration and Mental Health Research

Quantum computing and supercomputing are two distinct approaches that can be used to solve complex computational problems [...].

Extension of the iSoLF implicit-solvent coarse-grained model for multicomponent lipid bilayers

iSoLF is a coarse-grained (CG) model for lipid molecules with the implicit-solvent approximation used in molecular dynamics (MD) simulations of biological membranes. Using the original iSoLF (iSoLFv1), MD simulations of lipid bilayers consisting of either POPC or DPPC and these bilayers, including membrane proteins, can be performed. Here, we improve the original model, explicitly treating the electrostatic interactions between different lipid molecules and adding CG particle types. As a result, the available lipid types increase to 30. To parameterize the potential functions of the new model, we performed all-atom MD simulations of each lipid at three different temperatures using the CHARMM36 force field and the modified TIP3P model. Then, we parameterized both the bonded and non-bonded interactions to fit the area per lipid and the membrane thickness of each lipid bilayer by using the multistate Boltzmann Inversion method. The final model reproduces the area per lipid and the membrane thickness of each lipid bilayer at the three temperatures. We also examined the applicability of the new model, iSoLFv2, to simulate the phase behaviors of mixtures of DOPC and DPPC at different concentrations. The simulation results with iSoLFv2 are consistent with those using Dry Martini and Martini 3, although iSoLFv2 requires much fewer computations. iSoLFv2 has been implemented in the GENESIS MD software and is publicly available.

The machinery of vesicle fusion

The compartmentalization of eukaryotic cells is reliant on the fidelity of vesicle-mediated intracellular transport. Vesicles deliver their cargo via membrane fusion, a process requiring membrane tethers, Sec1/Munc18 (SM) proteins, and SNAREs. These components function in concert to ensure that membrane fusion is efficient and accurate, but the mechanisms underlying their cooperative action are still in many respects mysterious. In this brief review, we highlight recent progress toward a more integrative understanding of the vesicle fusion machinery. We focus particular attention on cryo-electron microscopy structures of intact multisubunit tethers in complex with SNAREs or SM proteins, as well as a structure of an SM protein bound to multiple SNAREs. The insights gained from this work emphasize the advantages of studying the fusion machinery intact and in context.

SNARE Protein Snc1 Is Essential for Vesicle Trafficking, Membrane Fusion and Protein Secretion in Fungi

Fungi are an important group of microorganisms that play crucial roles in a variety of ecological and biotechnological processes. Fungi depend on intracellular protein trafficking, which involves moving proteins from their site of synthesis to the final destination within or outside the cell. The soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNARE) proteins are vital components of vesicle trafficking and membrane fusion, ultimately leading to the release of cargos to the target destination. The v-SNARE (vesicle-associated SNARE) Snc1 is responsible for anterograde and retrograde vesicle trafficking between the plasma membrane (PM) and Golgi. It allows for the fusion of exocytic vesicles to the PM and the subsequent recycling of Golgi-localized proteins back to the Golgi via three distinct and parallel recycling pathways. This recycling process requires several components, including a phospholipid flippase (Drs2-Cdc50), an F-box protein (Rcy1), a sorting nexin (Snx4-Atg20), a retromer submit, and the COPI coat complex. Snc1 interacts with exocytic SNAREs (Sso1/2, Sec9) and the exocytic complex to complete the process of exocytosis. It also interacts with endocytic SNAREs (Tlg1 and Tlg2) during endocytic trafficking. Snc1 has been extensively investigated in fungi and has been found to play crucial roles in various aspects of intracellular protein trafficking. When Snc1 is overexpressed alone or in combination with some key secretory components, it results in enhanced protein production. This article will cover the role of Snc1 in the anterograde and retrograde trafficking of fungi and its interactions with other proteins for efficient cellular transportation.

Intrinsic Disorder in α-Synuclein Regulates the Exocytotic Fusion Pore Transition

Today, it is widely accepted that intrinsic disorder is strongly related to the cell cycle, during mitosis, differentiation, and apoptosis. Of particular interest are hybrid proteins possessing both structured and unstructured domains that are critical in human health and disease, such as α-synuclein. In this work, we describe how α-synuclein interacts with the nascent fusion pore as it evolves toward expansion. We unveil the key role played by its intrinsically disordered region as a thermodynamic regulator of the nucleation-expansion energy barrier. By analyzing a truncated variant of α-synuclein that lacks the disordered region, we find that the landscape of protein interactions with PIP₂ and POPS lipids is highly altered, ultimately increasing the energy cost for the fusion pore to transit from nucleation to expansion. We conclude that the intrinsically disordered region in full-length α-synuclein recognizes and allocates pivotal protein:lipid interactions during membrane remodeling in the first stages of the fusion pore.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

Forty years of the adrenal chromaffin cell through ISCCB meetings around the world

This historical review focuses on the evolution of the knowledge accumulated during the last two centuries on the biology of the adrenal medulla gland and its chromaffin cells (CCs). The review emerged in the context of a series of meetings that started on the Spanish island of Ibiza in 1982 with the name of the International Symposium on Chromaffin Cell Biology (ISCCB). Hence, the review is divided into two periods namely, before 1982 and from this year to 2022, when the 21st ISCCB meeting was just held in Hamburg, Germany. The first historical period extends back to 1852 when Albert Kölliker first described the fine structure and function of the adrenal medulla. Subsequently, the adrenal staining with chromate salts identified the CCs; this was followed by the establishment of the embryological origin of the adrenal medulla, and the identification of adrenaline-storing vesicles. By the end of the nineteenth century, the basic morphology, histochemistry, and embryology of the adrenal gland were known. The twentieth century began with breakthrough findings namely, the experiment of Elliott suggesting that adrenaline was the sympathetic neurotransmitter, the isolation of pure adrenaline, and the deciphering of its molecular structure and chemical synthesis in the laboratory. In the 1950s, Blaschko isolated the catecholamine-storing vesicles from adrenal medullary extracts. This switched the interest in CCs as models of sympathetic neurons with an explosion of studies concerning their functions, i.e., uptake of catecholamines by chromaffin vesicles through a specific coupled transport system; the identification of several vesicle components in addition to catecholamines including chromogranins, ATP, opioids, and other neuropeptides; the calcium-dependence of the release of catecholamines; the underlying mechanism of exocytosis of this release, as indicated by the co-release of proteins; the cross-talk between the adrenal cortex and the medulla; and the emission of neurite-like processes by CCs in culture, among other numerous findings. The 1980s began with the introduction of new high-resolution techniques such as patch-clamp, calcium probes, marine toxins-targeting ion channels and receptors, confocal microscopy, or amperometry. In this frame of technological advances at the Ibiza ISCCB meeting in 1982, 11 senior researchers in the field predicted a notable increase in our knowledge in the field of CCs and the adrenal medulla; this cumulative knowledge that occurred in the last 40 years of history of the CC is succinctly described in the second part of this historical review. It deals with cell excitability, ion channel currents, the exocytotic fusion pore, the handling of calcium ions by CCs, the kinetics of exocytosis and endocytosis, the exocytotic machinery, and the life cycle of secretory vesicles. These concepts together with studies on the dynamics of membrane fusion with super-resolution imaging techniques at the single-protein level were extensively reviewed by top scientists in the field at the 21st ISCCB meeting in Hamburg in the summer of 2022; this frontier topic is also briefly reviewed here. Many of the concepts arising from those studies contributed to our present understanding of synaptic transmission. This has been studied in physiological or pathophysiological conditions, in CCs from animal disease models. In conclusion, the lessons we have learned from CC biology as a peripheral model for brain and brain disease pertain more than ever to cutting-edge research in neurobiology. In the 22nd ISCCB meeting in Israel in 2024 that Uri Asheri is organizing, we will have the opportunity of seeing the progress of the questions posed in Ibiza, and on other questions that undoubtedly will arise.

Molecular Dynamics Simulations of the Proteins Regulating Synaptic Vesicle Fusion

Neuronal transmitters are packaged in synaptic vesicles (SVs) and released by the fusion of SVs with the presynaptic membrane (PM). An inflow of Ca²⁺ into the nerve terminal triggers fusion, and the SV-associated protein Synaptotagmin 1 (Syt1) serves as a Ca²⁺ sensor. In preparation for fusion, SVs become attached to the PM by the SNARE protein complex, a coiled-coil bundle that exerts the force overcoming SV-PM repulsion. A cytosolic protein Complexin (Cpx) attaches to the SNARE complex and differentially regulates the evoked and spontaneous release components. It is still debated how the dynamic interactions of Syt1, SNARE proteins and Cpx lead to fusion. This problem is confounded by heterogeneity in the conformational states of the prefusion protein-lipid complex and by the lack of tools to experimentally monitor the rapid conformational transitions of the complex, which occur at a sub-millisecond scale. However, these complications can be overcome employing molecular dynamics (MDs), a computational approach that enables simulating interactions and conformational transitions of proteins and lipids. This review discusses the use of molecular dynamics for the investigation of the pre-fusion protein-lipid complex. We discuss the dynamics of the SNARE complex between lipid bilayers, as well as the interactions of Syt1 with lipids and SNARE proteins, and Cpx regulating the assembly of the SNARE complex.

Cofilin-Membrane Interactions: Electrostatic Effects in Phosphoinositide Lipid Binding

The actin cytoskeleton interacts with the cell membrane primarily through the indirect interactions of actin-binding proteins such as cofilin-1. The molecular mechanisms underlying the specific interactions of cofilin-1 with membrane lipids are still unclear. Here, we performed coarse-grain molecular dynamics simulations of cofilin-1 with complex lipid bilayers to analyze the specificity of protein-lipid interactions. We observed the maximal interactions with phosphoinositide (PIP) lipids, especially PIP₂ and PIP₃ lipids. A good match was observed between the residues predicted to interact and previous experimental studies. The clustering of PIP lipids around the membrane bound protein leads to an overall lipid demixing and gives rise to persistent membrane curvature. Further, through a series of control simulations, we observe that both electrostatics and geometry are critical for specificity of lipid binding. Our current study is a step towards understanding the physico-chemical basis of cofilin-PIP lipid interactions.

The Role of α-Synuclein in SNARE-mediated Synaptic Vesicle Fusion

Neuronal communication depends on exquisitely regulated membrane fusion between synaptic vesicles and presynaptic neurons, which results in neurotransmitter release in precisely timed patterns. Presynaptic dysfunctions are known to occur prior to the onset of neurodegenerative diseases, including Parkinson's disease. Synaptic accumulation of α-synuclein (α-Syn) oligomers has been implicated in the pathway leading to such outcomes. α-Syn oligomers exert aberrant effects on presynaptic fusion machinery through their interactions with synaptic vesicles and proteins. Here, we summarize in vitro bulk and single-vesicle assays for investigating the functions of α-Syn monomers and oligomers in synaptic vesicle fusion and then discuss the current understanding of the roles of α-Syn monomers and oligomers in synaptic vesicle fusion. Finally, we suggest a new therapeutic avenue specifically targeting the mechanisms of α-Syn oligomer toxicity rather than the oligomer itself.

Computational Advances of Protein/Neurotransmitter-membrane Interactions Involved in Vesicle Fusion and Neurotransmitter Release

Vesicle fusion is of crucial importance to neuronal communication at neuron terminals. The exquisite but complex fusion machinery for neurotransmitter release is tightly controlled and regulated by protein/neurotransmitter-membrane interactions. Computational 'microscopies', in particular molecular dynamics simulations and related techniques, have provided notable insight into the physiological process over the past decades, and have made enormous contributions to fields such as neurology, pharmacology and pathophysiology. Here we review the computational advances of protein/neurotransmitter-membrane interactions related to presynaptic vesicle-membrane fusion and neurotransmitter release, and outline the in silico challenges ahead for understanding this important physiological process.

The Core Complex of the Ca2+-Triggered Presynaptic Fusion Machinery

Synaptic neurotransmitter release is mediated by an orchestra of presynaptic proteins that precisely control and trigger fusion between synaptic vesicles and the neuron terminal at the active zone upon the arrival of an action potential. Critical to this process are the neuronal SNAREs (Soluble N-ethylmaleimide sensitive factor Attachment protein REceptor), the Ca2+-sensor synaptotagmin, the activator/regulator complexin, and other factors. Here, we review the interactions between the SNARE complex and synaptotagmin, with focus on the so-called primary interface between synaptotagmin and the SNARE complex that has been validated in terms of its physiological relevance. We discuss several other but less validated interfaces as well, including the so-called tripartite interface, and we discuss the pros and cons for these possible alternative interfaces. We also present new molecular dynamics simulations of the tripartite interface and new data of an inhibitor of the primary interface in a reconstituted system of synaptic vesicle fusion.

Analysis of tripartite Synaptotagmin-1-SNARE-complexin-1 complexes in solution

Characterizing interactions of Synaptotagmin-1 with the SNARE complex is crucial to understand the mechanism of neurotransmitter release. X-ray crystallography revealed how the Synaptotagmin-1 C2 B domain binds to the SNARE complex through a so-called primary interface and to a complexin-1-SNARE complex through a so-called tripartite interface. Mutagenesis and electrophysiology supported the functional relevance of both interfaces, and extensive additional data validated the primary interface. However, ITC evidence suggesting that binding via the tripartite interface occurs in solution was called into question by subsequent NMR data. Here, we describe joint efforts to address this apparent contradiction. Using the same ITC approach with the same C2 B domain mutant used previously (C2 BKA-Q ) but including ion exchange chromatography to purify it, which is crucial to remove polyacidic contaminants, we were unable to observe the substantial endothermic ITC signal that was previously attributed to binding of this mutant to the complexin-1-SNARE complex through the tripartite interface. We were also unable to detect substantial populations of the tripartite interface in NMR analyses of the ITC samples or in measurements of paramagnetic relaxation effects, despite the high sensitivity of this method to detect weak protein complexes. However, these experiments do not rule out the possibility of very low affinity (KD  > 1 mm) binding through this interface. These results emphasize the need to develop methods to characterize the structure of synaptotagmin-1-SNARE complexes between two membranes and to perform further structure-function analyses to establish the physiological relevance of the tripartite interface.

SNARE Proteins in Synaptic Vesicle Fusion

Neurotransmitters are stored in small membrane-bound vesicles at synapses; a subset of synaptic vesicles is docked at release sites. Fusion of docked vesicles with the plasma membrane releases neurotransmitters. Membrane fusion at synapses, as well as all trafficking steps of the secretory pathway, is mediated by SNARE proteins. The SNAREs are the minimal fusion machinery. They zipper from N-termini to membrane-anchored C-termini to form a 4-helix bundle that forces the apposed membranes to fuse. At synapses, the SNAREs comprise a single helix from syntaxin and synaptobrevin; SNAP-25 contributes the other two helices to complete the bundle. Unc13 mediates synaptic vesicle docking and converts syntaxin into the permissive "open" configuration. The SM protein, Unc18, is required to initiate and proofread SNARE assembly. The SNAREs are then held in a half-zippered state by synaptotagmin and complexin. Calcium removes the synaptotagmin and complexin block, and the SNAREs drive vesicle fusion. After fusion, NSF and alpha-SNAP unwind the SNAREs and thereby recharge the system for further rounds of fusion. In this chapter, we will describe the discovery of the SNAREs, their relevant structural features, models for their function, and the central role of Unc18. In addition, we will touch upon the regulation of SNARE complex formation by Unc13, complexin, and synaptotagmin.

On the difficulties of characterizing weak protein interactions that are critical for neurotransmitter release

The mechanism of neurotransmitter release has been extensively characterized, showing that vesicle fusion is mediated by the SNARE complex formed by syntaxin-1, SNAP-25 and synaptobrevin. This complex is disassembled by N-ethylmaleimide sensitive factor (NSF) and SNAPs to recycle the SNAREs, whereas Munc18-1 and Munc13s organize SNARE complex assembly in an NSF-SNAP-resistant manner. Synaptotagmin-1 acts as the Ca²⁺ sensor that triggers exocytosis in a tight interplay with the SNAREs and complexins. Here, we review technical aspects associated with investigation of protein interactions underlying these steps, which is hindered because the release machinery is assembled between two membranes and is highly dynamic. Moreover, weak interactions, which are difficult to characterize, play key roles in neurotransmitter release, for instance by lowering energy barriers that need to be overcome in this highly regulated process. We illustrate the crucial role that structural biology has played in uncovering mechanisms underlying neurotransmitter release, but also discuss the importance of considering the limitations of the techniques used, including lessons learned from research in our lab and others. In particular, we emphasize: (a) the promiscuity of some protein sequences, including membrane-binding regions that can mediate irrelevant interactions with proteins in the absence of their native targets; (b) the need to ensure that weak interactions observed in crystal structures are biologically relevant; and (c) the limitations of isothermal titration calorimetry to analyze weak interactions. Finally, we stress that even studies that required re-interpretation often helped to move the field forward by improving our understanding of the system and providing testable hypotheses.

Energetics, kinetics, and pathways of SNARE assembly in membrane fusion

Fusion of transmitter-containing vesicles with plasma membranes at the synaptic and neuromuscular junctions mediates neurotransmission and muscle contractions, respectively, thereby underlying all thoughts and actions. The fusion process is driven by the coupled folding and assembly of three synaptic SNARE proteins--syntaxin-1 and SNAP-25 on the target plasma membrane (t-SNAREs) and VAMP2 on the vesicular membrane (v-SNARE) into a four-helix bundle. Their assembly is chaperoned by Munc18-1 and many other proteins to achieve the speed and accuracy required for neurotransmission. However, the physiological pathway of SNARE assembly and its coupling to membrane fusion remains unclear. Here, we review recent progress in understanding SNARE assembly and membrane fusion, with a focus on results obtained by single-molecule manipulation approaches and electric recordings of single fusion pores. We describe two pathways of synaptic SNARE assembly, their associated intermediates, energetics, and kinetics. Assembly of the three SNAREs in vitro begins with the formation of a t-SNARE binary complex, on which VAMP2 folds in a stepwise zipper-like fashion. Munc18-1 significantly alters the SNARE assembly pathway: syntaxin-1 and VAMP2 first bind on the surface of Munc18-1 to form a template complex, with which SNAP-25 associates to conclude SNARE assembly and displace Munc18-1. During membrane fusion, multiple trans-SNARE complexes cooperate to open a dynamic fusion pore in a manner dependent upon their copy number and zippering states. Together, these results demonstrate that stepwise and cooperative SNARE assembly drive stagewise membrane fusion.