Membrane Docking of the Synaptotagmin 7 C2A Domain: Computation Reveals Interplay between Electrostatic and Hydrophobic Contributions

Ca2+-dependent lipid preferences shape synaptotagmin-1 C2A and C2B dynamics: Insights from experiments and simulations

Signal transmission between neurons requires exocytosis of neurotransmitters from the lumen of synaptic vesicles into the synaptic cleft. Following an influx of Ca2+, this process is facilitated by the Ca2+ sensor synaptotagmin-1. The underlying mechanisms involve electrostatic and hydrophobic interactions tuning the lipid preferences of the two C2 domains of synaptotagmin-1; however, the details are still controversially discussed. We, therefore, follow a multidisciplinary approach and characterize lipid and membrane binding of the isolated C2A and C2B domains. We first target interactions with individual lipid species, and then study interactions with model membranes of liposomes. Finally, we perform molecular dynamics simulations to unravel differences in membrane binding. We found that both C2 domains, as a response to Ca2+, insert into the lipid membrane; however, C2A adopts a more perpendicular orientation while C2B remains parallel. These findings allow us to propose a mechanism for synaptotagmin-1 during membrane fusion.

A conserved electrostatic membrane-binding surface in synaptotagmin-like proteins revealed using molecular phylogenetic analysis and homology modeling

Protein structure prediction has emerged as a core technology for understanding biomolecules and their interactions. Here, we combine homology-based structure prediction with molecular phylogenetic analysis to study the evolution of electrostatic membrane binding among the vertebrate synaptotagmin-like protein (Slp) family. Slp family proteins play key roles in the membrane trafficking of large dense-core secretory vesicles. Our previous experimental and computational study found that the C2A domain of Slp-4 (also called granuphilin) binds with high affinity to anionic phospholipids in the cytoplasmic leaflet of the plasma membrane through a large positively charged protein surface centered on a cluster of phosphoinositide-binding lysine residues. Because this surface contributes greatly to Slp-4 C2A domain membrane binding, we hypothesized that the net charge on the surface might be evolutionarily conserved. To test this hypothesis, the known C2A sequences of Slp-4 among vertebrates were organized by class (from mammalia to pisces) using molecular phylogenetic analysis. Consensus sequences for each class were then identified and used to generate homology structures, from which Poisson-Boltzmann electrostatic potentials were calculated. For comparison, homology structures and electrostatic potentials were also calculated for the five human Slp protein family members. The results demonstrate that the charge on the membrane-binding surface is highly conserved throughout the evolution of Slp-4, and more highly conserved than many individual residues among the human Slp family paralogs. Such molecular phylogenetic-driven computational analysis can help to describe the evolution of electrostatic interactions between proteins and membranes which are crucial for their function.

A Conserved Electrostatic Membrane-Binding Surface in Synaptotagmin-Like Proteins Revealed Using Molecular Phylogenetic Analysis and Homology Modeling

Protein structure prediction has emerged as a core technology for understanding biomolecules and their interactions. Here, we combine homology-based structure prediction with molecular phylogenetic analysis to study the evolution of electrostatic membrane binding among vertebrate synaptotagmin-like proteins (Slps). Slp family proteins play key roles in the membrane trafficking of large dense-core secretory vesicles. Our previous experimental and computational study found that the C2A domain of Slp-4 (also called granuphilin) binds with high affinity to anionic phospholipids in the cytoplasmic leaflet of the plasma membrane through a large positively charged protein surface centered on a cluster of phosphoinositide-binding lysine residues. Because this surface contributes greatly to Slp-4 C2A domain membrane binding, we hypothesized that the net charge on the surface might be evolutionarily conserved. To test this hypothesis, the known C2A sequences of Slp-4 among vertebrates were organized by class (from mammalia to pisces) using molecular phylogenetic analysis. Consensus sequences for each class were then identified and used to generate homology structures, from which Poisson-Boltzmann electrostatic potentials were calculated. For comparison, homology structures and electrostatic potentials were also calculated for the five human Slp protein family members. The results demonstrate that the charge on the membrane-binding surface is highly conserved throughout the evolution of Slp-4, and more highly conserved than many individual residues among the human Slp family paralogs. Such molecular phylogenetic-driven computational analysis can help to describe the evolution of electrostatic interactions between proteins and membranes which are crucial for their function.

Specific interactions of peripheral membrane proteins with lipids: what can molecular simulations show us?

Peripheral membrane proteins (PMPs) can reversibly and specifically bind to biological membranes to carry out functions such as cell signalling, enzymatic activity, or membrane remodelling. Structures of these proteins and of their lipid-binding domains are typically solved in a soluble form, sometimes with a lipid or lipid headgroup at the binding site. To provide a detailed molecular view of PMP interactions with the membrane, computational methods such as molecular dynamics (MD) simulations can be applied. Here, we outline recent attempts to characterise these binding interactions, focusing on both intracellular proteins, such as phosphatidylinositol phosphate (PIP)-binding domains, and extracellular proteins such as glycolipid-binding bacterial exotoxins. We compare methods used to identify and analyse lipid-binding sites from simulation data and highlight recent work characterising the energetics of these interactions using free energy calculations. We describe how improvements in methodologies and computing power will help MD simulations to continue to contribute to this field in the future.

Binding of Ca2+-independent C2 domains to lipid membranes: A multi-scale molecular dynamics study

C2 domains facilitate protein interactions with lipid bilayers in either a Ca²⁺-dependent or -independent manner. We used molecular dynamics (MD) simulations to explore six Ca²⁺-independent C2 domains, from KIBRA, PI3KC2α, RIM2, PTEN, SHIP2, and Smurf2. In coarse-grained MD simulations these C2 domains formed transient interactions with zwitterionic bilayers, compared with longer-lived interactions with anionic bilayers containing phosphatidylinositol bisphosphate (PIP₂). Type I C2 domains bound non-canonically via the front, back, or side of the β sandwich, whereas type II C2 domains bound canonically, via the top loops. C2 domains interacted strongly with membranes containing PIP₂, causing bound anionic lipids to cluster around the protein. Binding modes were refined via atomistic simulations. For PTEN and SHIP2, CG simulations of their phosphatase plus C2 domains with PIP₂-containing bilayers were also performed, and the roles of the two domains in membrane localization compared. These studies establish a simulation protocol for membrane-recognition proteins.

•

Binding of Ca²⁺-independent C2 domains to membranes was explored by MD simulation

•

C2 domains from KIBRA, PI3KC2α, RIM2, PTEN, SHIP2, and Smurf2 were compared

•

C2 domains formed longer-lived interactions with lipid bilayers containing PIP₂

•

For PTEN and SHIP2, simulations of their phosphatase plus C2 domains were performed

Calcium-dependent docking of synaptic vesicles

The concentration of calcium ions in presynaptic terminals regulates transmitter release, but underlying mechanisms have remained unclear. Here we review recent studies that shed new light on this issue. Fast-freezing electron microscopy and total internal reflection fluorescence microscopy studies reveal complex calcium-dependent vesicle movements including docking on a millisecond time scale. Recordings from so-called 'simple synapses' indicate that calcium not only triggers exocytosis, but also modifies synaptic strength by controlling a final, rapid vesicle maturation step before release. Molecular studies identify several calcium-sensitive domains on Munc13 and on synaptotagmin-1 that are likely involved in bringing the vesicular and plasma membranes closer together in response to calcium elevation. Together, these results suggest that calcium-dependent vesicle docking occurs in a wide range of time domains and plays a crucial role in several phenomena including synaptic facilitation, post-tetanic potentiation, and neuromodulator-induced potentiation.

Multivalent lipid targeting by the calcium-independent C2A domain of synaptotagmin-like protein 4/granuphilin

Synaptotagmin-like protein 4 (Slp-4), also known as granuphilin, is a Rab effector responsible for docking secretory vesicles to the plasma membrane before exocytosis. Slp-4 binds vesicular Rab proteins via an N-terminal Slp homology domain, interacts with plasma membrane SNARE complex proteins via a central linker region, and contains tandem C-terminal C2 domains (C2A and C2B) with affinity for phosphatidylinositol-(4,5)-bisphosphate (PIP₂). The Slp-4 C2A domain binds with low nanomolar apparent affinity to PIP₂ in lipid vesicles that also contain background anionic lipids such as phosphatidylserine but much weaker when either the background anionic lipids or PIP₂ is removed. Through computational and experimental approaches, we show that this high-affinity membrane binding arises from concerted interaction at multiple sites on the C2A domain. In addition to a conserved PIP₂-selective lysine cluster, a larger cationic surface surrounding the cluster contributes substantially to the affinity for physiologically relevant lipid compositions. Although the K398A mutation in the lysine cluster blocks PIP₂ binding, this mutated protein domain retains the ability to bind physiological membranes in both a liposome-binding assay and MIN6 cells. Molecular dynamics simulations indicate several conformationally flexible loops that contribute to the nonspecific cationic surface. We also identify and characterize a covalently modified variant that arises through reactivity of the PIP₂-binding lysine cluster with endogenous bacterial compounds and binds weakly to membranes. Overall, multivalent lipid binding by the Slp-4 C2A domain provides selective recognition and high-affinity docking of large dense core secretory vesicles to the plasma membrane.

Otoferlin C2F Domain-Induced Changes in Membrane Structure Observed by Sum Frequency Generation

Proteins that contain C2 domains are involved in a variety of biological processes, including encoding of sound, cell signaling, and cell membrane repair. Of particular importance is the interface activity of the C-terminal C2F domain of otoferlin due to the pathological mutations known to significantly disrupt the protein's lipid membrane interface binding activity, resulting in hearing loss. Therefore, there is a critical need to define the geometry and positions of functionally important sites and structures at the otoferlin-lipid membrane interface. Here, we describe the first in situ probe of the protein orientation of otoferlin's C2F domain interacting with a cell membrane surface. To identify this protein's orientation at the lipid interface, we applied sum frequency generation (SFG) vibrational spectroscopy and coupled it with simulated SFG spectra to observe and quantify the otoferlin C2F domain interacting with model lipid membranes. A model cell membrane was built with equal amounts of phosphatidylserine and phosphatidylcholine. SFG measurements of the lipids that make up the model membrane indicate a 62% increase in amplitude from the SFG signal near 2075 cm^-1 upon protein interaction, suggesting domain-induced changes in the orientation of the lipids and possible membrane curvature. This increase is related to lipid ordering caused by the docking interaction of the otoferlin C2F domain. SFG spectra taken from the amide-I region contain features near 1630 and 1670 cm^-1 related to the C2F domains beta-sandwich secondary structure, thus indicating that the domain binds in a specific orientation. By mapping the simulated SFG spectra to the experimentally collected SFG spectra, we found the C2F domain of otoferlin orients 22° normal to the lipid surface. This information allows us to map what portion of the domain directly interacts with the lipid membrane.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Membrane-Binding Cooperativity and Coinsertion by C2AB Tandem Domains of Synaptotagmins 1 and 7

Synaptotagmin-1 (Syt-1) and synaptotagmin-7 (Syt-7) contain analogous tandem C2 domains, C2A and C2B, which together sense Ca²⁺ to bind membranes and promote the stabilization of exocytotic fusion pores. Syt-1 triggers fast release of neurotransmitters, whereas Syt-7 functions in processes that involve lower Ca²⁺ concentrations such as hormone secretion. Syt-1 C2 domains are reported to bind membranes cooperatively, based on the observation that they penetrate farther into membranes as the C2AB tandem than as individual C2 domains. In contrast, we previously suggested that the two C2 domains of Syt-7 bind membranes independently, based in part on measurements of their liposome dissociation kinetics. Here, we investigated C2A-C2B interdomain cooperativity with Syt-1 and Syt-7 using directly comparable measurements. Equilibrium Ca²⁺ titrations demonstrate that the Syt-7 C2AB tandem binds liposomes lacking phosphatidylinositol-4,5-bisphosphate (PIP₂) with greater Ca²⁺ sensitivity than either of its individual domains and binds to membranes containing PIP₂ even in the absence of Ca²⁺. Stopped-flow kinetic measurements show differences in cooperativity between Syt-1 and Syt-7: Syt-1 C2AB dissociates from PIP₂-free liposomes much more slowly than either of its individual C2 domains, indicating cooperativity, whereas the major population of Syt-7 C2AB has a dissociation rate comparable to its C2A domain, suggesting a lack of cooperativity. A minor subpopulation of Syt-7 C2AB dissociates at a slower rate, which could be due to a small cooperative component and/or liposome clustering. Measurements using an environment-sensitive fluorescent probe indicate that the Syt-7 C2B domain inserts deeply into membranes as part of the C2AB tandem, similar to the coinsertion previously reported for Syt-1. Overall, coinsertion of C2A and C2B domains is coupled to cooperative energetic effects in Syt-1 to a much greater extent than in Syt-7. The difference can be understood in terms of the relative contributions of C2A and C2B domains toward membrane binding in the two proteins.

The high-affinity calcium sensor synaptotagmin-7 serves multiple roles in regulated exocytosis

MacDougall et al. review the structure and function of the calcium sensor synaptotagmin-7 in exocytosis.

Synaptotagmin (Syt) proteins comprise a 17-member family, many of which trigger exocytosis in response to calcium. Historically, most studies have focused on the isoform Syt-1, which serves as the primary calcium sensor in synchronous neurotransmitter release. Recently, Syt-7 has become a topic of broad interest because of its extreme calcium sensitivity and diversity of roles in a wide range of cell types. Here, we review the known and emerging roles of Syt-7 in various contexts and stress the importance of its actions. Unique functions of Syt-7 are discussed in light of recent imaging, electrophysiological, and computational studies. Particular emphasis is placed on Syt-7–dependent regulation of synaptic transmission and neuroendocrine cell secretion. Finally, based on biochemical and structural data, we propose a mechanism to link Syt-7’s role in membrane fusion with its role in subsequent fusion pore expansion via strong calcium-dependent phospholipid binding.

A Polarizable Atomic Multipole-Based Force Field for Molecular Dynamics Simulations of Anionic Lipids

In all of the classical force fields, electrostatic interaction is simply treated and explicit electronic polarizability is neglected. The condensed-phase polarization, relative to the gas-phase charge distributions, is commonly accounted for in an average way by increasing the atomic charges, which remain fixed throughout simulations. Based on the lipid polarizable force field DMPC and following the same framework as Atomic Multipole Optimized Energetics for BiomoleculAr (AMOEBA) simulation, the present effort expands the force field to new anionic lipid models, in which the new lipids contain DMPG and POPS. The parameters are compatible with the AMOEBA force field, which includes water, ions, proteins, etc. The charge distribution of each atom is represented by the permanent atomic monopole, dipole and quadrupole moments, which are derived from the ab initio gas phase calculations. Many-body polarization including the inter- and intramolecular polarization is modeled in a consistent manner with distributed atomic polarizabilities. Molecular dynamics simulations of the two aqueous DMPG and POPS membrane bilayer systems, consisting of 72 lipids with water molecules, were then carried out to validate the force field parameters. Membrane width, area per lipid, volume per lipid, deuterium order parameters, electron density profile, electrostatic potential difference between the center of the bilayer and water are all calculated, and compared with limited experimental data.

Lipid-Coated Gold Nanoparticles and FRET Allow Sensitive Monitoring of Liposome Clustering Mediated by the Synaptotagmin-7 C2A Domain

Synaptotagmin (Syt) family proteins contain tandem C2 domains, C2A and C2B, which insert into anionic membranes in response to increased cytosolic Ca²⁺ concentration and facilitate exocytosis in neuronal and endocrine cells. The C2A domain from Syt7 binds lipid membranes much more tightly than the corresponding domain from Syt1, but the implications of this difference for protein function are not yet clear. In particular, the ability of the isolated Syt7 C2A domain to initiate membrane apposition and/or aggregation has been previously unexplored. Here, we demonstrate that Syt7 C2A induces apposition and aggregation of liposomes using Förster resonance energy transfer (FRET) assays, dynamic light scattering, and spectroscopic techniques involving lipid-coated gold nanoparticles (LCAuNPs). Protein-membrane binding, membrane apposition, and macroscopic aggregation are three separate phenomena with distinct Ca²⁺ requirements: the threshold Ca²⁺ concentration for membrane binding is lowest, followed by apposition and aggregation. However, aggregation is highly sensitive to protein concentration and can occur even at submicromolar Syt7 C2A; thus, highly sensitive assays are needed for measuring apposition without complications arising from aggregation. Notably, the localized surface plasmon resonance of the LCAuNP is sensitive to ≤10 nM Syt7 C2A concentrations. Furthermore, when the LCAuNPs were added into a FRET-based liposome apposition assay, the resultant energy transfer increased; possible explanations are discussed. Overall, LCAuNP-based methods allow for highly sensitive detection of protein-induced membrane apposition under conditions that miminize large-scale aggregation.

Molecular determinants of the N-terminal acetyltransferase Naa60 anchoring to the Golgi membrane

Nα-Acetyltransferase 60 (Naa60 or NatF) was recently identified as an unconventional N-terminal acetyltransferase (NAT) because it localizes to organelles, in particular the Golgi apparatus, and has a preference for acetylating N termini of the transmembrane proteins. This knowledge challenged the prevailing view of N-terminal acetylation as a co-translational ribosome-associated process and suggested a new mechanistic functioning for the enzymes responsible for this increasingly recognized protein modification. Crystallography studies on Naa60 were unable to resolve the C-terminal tail of Naa60, which is responsible for the organellar localization. Here, we combined modeling, in vitro assays, and cellular localization studies to investigate the secondary structure and membrane interacting capacity of Naa60. The results show that Naa60 is a peripheral membrane protein. Two amphipathic helices within the Naa60 C terminus bind the membrane directly in a parallel position relative to the lipid bilayer via hydrophobic and electrostatic interactions. A peptide corresponding to the C terminus was unstructured in solution and only folded into an α-helical conformation in the presence of liposomes. Computational modeling and cellular mutational analysis revealed the hydrophobic face of two α-helices to be critical for membranous localization. Furthermore, we found a strong and specific binding preference of Naa60 toward membranes containing the phosphatidylinositol PI(4)P, thus possibly explaining the primary residency of Naa60 at the PI(4)P-rich Golgi. In conclusion, we have defined the mode of cytosolic Naa60 anchoring to the Golgi apparatus, most likely occurring post-translationally and specifically facilitating post-translational N-terminal acetylation of many transmembrane proteins.

Differential Membrane Binding Mechanics of Synaptotagmin Isoforms Observed in Atomic Detail

Synaptotagmin (Syt) is a membrane-associated protein involved in vesicle fusion through the SNARE complex that is found throughout the human body in 17 different isoforms. These isoforms have two membrane-binding C2 domains, which sense Ca²⁺ and thereby promote anionic membrane binding and lead to vesicle fusion. Through molecular dynamics simulations using the highly mobile membrane mimetic acclerated bilayer model, we have investigated how small protein sequence changes in the Ca²⁺-binding loops of the C2 domains may give rise to the experimentally determined difference in binding kinetics between Syt-1 and Syt-7 isoforms. Syt-7 C2 domains are found to form more close contacts with anionic phospholipid headgroups, particularly in loop 1, where an additional positive charge in Syt-7 draws the loop closer to the membrane and causes the anchoring residue F167 to insert deeper into the bilayer than the corresponding methionine in Syt-1 (M173). By performing additional replica exchange umbrella sampling calculations, we demonstrate that these additional contacts increase the energetic cost of unbinding the Syt-7 C2 domains from the bilayer, causing them to unbind more slowly than their counterparts in Syt-1.

The Effects of Ca2+ Concentration and E200K Mutation on the Aggregation Propensity of PrPC: A Computational Study

The propensity of cellular prion protein to aggregation is reputed essential for the initiation of the amyloid cascade that ultimately lead to the accumulation of neurotoxic aggregates. In this paper, we extended and applied an already reported computational workflow [Proteins 2015; 83: 1751–1765] to elucidate in details the aggregation propensity of PrP protein systems including wild type, wild type treated at different [Ca²⁺] and E200K mutant. The application of the computational procedure to two segments of PrP^C, i.e. 125–228 and 120–231, allowed to emphasize how the inclusion of complete C-terminus and last portion (120–126) of the neurotoxic segment 106–126 may be crucial to unveil significant and unexpected interaction properties. Indeed, the anchoring of N-terminus on H2 domain detected in the wild type resulted to be disrupted upon either E200K mutation or Ca²⁺ binding, and to unbury hydrophobic spots on the PrP^C surface. A peculiar dinuclear Ca²⁺ binding motif formed by the C-terminus and the S2-H2 loop was detected for [Ca²⁺] > 5 mM and showed similarities with binding motifs retraced in other protein systems, thus suggesting a possible functional meaning for its formation. Therefore, we potentiated the computational procedure by including a tool that clusterize the minima of molecular interaction fields of a proteinand delimit the regions of space with higher hydrophobic or higher hydrophilic character, hence, more likely involved in the self-assembly process. Plausible models for the self-assembly of either the E200K mutated or Ca²⁺-bound PrP^C were sketched and discussed. The present investigation provides for structure-based information and new prompts that may represent a starting point for future experimental or computational works on the PrP^C aggregation.