Interactions between neuronal fusion proteins explored by molecular dynamics

From complex data to clear insights: visualizing molecular dynamics trajectories

Advances in simulations, combined with technological developments in high-performance computing, have made it possible to produce a physically accurate dynamic representation of complex biological systems involving millions to billions of atoms over increasingly long simulation times. The analysis of these computed simulations is crucial, involving the interpretation of structural and dynamic data to gain insights into the underlying biological processes. However, this analysis becomes increasingly challenging due to the complexity of the generated systems with a large number of individual runs, ranging from hundreds to thousands of trajectories. This massive increase in raw simulation data creates additional processing and visualization challenges. Effective visualization techniques play a vital role in facilitating the analysis and interpretation of molecular dynamics simulations. In this paper, we focus mainly on the techniques and tools that can be used for visualization of molecular dynamics simulations, among which we highlight the few approaches used specifically for this purpose, discussing their advantages and limitations, and addressing the future challenges of molecular dynamics visualization.

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.

Computational methods to study enveloped viral entry

The protein-membrane interactions that mediate viral infection occur via loosely ordered, transient assemblies, creating challenges for high-resolution structure determination. Computational methods and in particular molecular dynamics simulation have thus become important adjuncts for integrating experimental data, developing mechanistic models, and suggesting testable hypotheses regarding viral function. However, the large molecular scales of virus-host interaction also create challenges for detailed molecular simulation. For this reason, continuum membrane models have played a large historical role, although they have become less favored for high-resolution models of protein assemblies and lipid organization. Here, we review recent progress in the field, with an emphasis on the insight that has been gained using a mixture of coarse-grained and atomic-resolution molecular dynamics simulations. Based on successes and challenges to date, we suggest a multiresolution strategy that should yield the best mixture of computational efficiency and physical fidelity. This strategy may facilitate further simulations of viral entry by a broader range of viruses, helping illuminate the diversity of viral entry strategies and the essential common elements that can be targeted for antiviral therapies.

Simulations Reveal Multiple Intermediates in the Unzipping Mechanism of Neuronal SNARE Complex

The assembling of the soluble N-ethylmaleimide-sensitive factor attachment protein receptor protein complex is a fundamental step in neuronal exocytosis, and it has been extensively studied in the last two decades. Yet, many details of this process remain inaccessible with the current experimental space and time resolution. Here, we study the zipping mechanism of the soluble N-ethylmaleimide-sensitive factor attachment protein receptor complex computationally by using a coarse-grained model. We explore the different pathways available and analyze their dependence on the computational model employed. We reveal and characterize multiple intermediate states, in agreement with previous experimental findings. We use our model to analyze the influence of single-residue mutations on the thermodynamics of the folding process.

Molecular Dynamics Simulations of the SNARE Complex

Molecular dynamics (MD) simulations enable in silico investigations of the dynamic behavior of proteins and protein complexes. Here, we describe MD simulations of the SNARE complex and its interactions with the neuronal protein complexin. Complexin is an effector of neuronal secretion that inhibits spontaneous fusion and is thought to clamp the fusion process via the interactions with the SNARE complex. We describe MD simulations of the SNARE complex alone and bound to complexin. The MD simulations under external forces imitating the repulsion between lipid bilayers enabled us to investigate unraveling and assembly of the SNARE complex.

A membrane-inserted structural model of the yeast mitofusin Fzo1

Mitofusins are large transmembrane GTPases of the dynamin-related protein family, and are required for the tethering and fusion of mitochondrial outer membranes. Their full-length structures remain unknown, which is a limiting factor in the study of outer membrane fusion. We investigated the structure and dynamics of the yeast mitofusin Fzo1 through a hybrid computational and experimental approach, combining molecular modelling and all-atom molecular dynamics simulations in a lipid bilayer with site-directed mutagenesis and in vivo functional assays. The predicted architecture of Fzo1 improves upon the current domain annotation, with a precise description of the helical spans linked by flexible hinges, which are likely of functional significance. In vivo site-directed mutagenesis validates salient aspects of this model, notably, the long-distance contacts and residues participating in hinges. GDP is predicted to interact with Fzo1 through the G1 and G4 motifs of the GTPase domain. The model reveals structural determinants critical for protein function, including regions that may be involved in GTPase domain-dependent rearrangements.

The influence of cell membrane and SNAP25 linker loop on the dynamics and unzipping of SNARE complex

The soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) complex is composed of three neuronal proteins VAMP2, Syntaxin and SNAP25, which plays a core role during the process of membrane fusion. The zipping assembly of the SNARE complex releases energies and drives the vesicle and cell membrane into close proximity. In this study, we use all-atom molecular dynamics simulations to probe the dynamics of SNARE and its unzipping process in the context of membrane at the atomistic details. Our results indicated that the NTD of SNARE core domain is relatively more stable than CTD, which is in agreement with previous experiments. More importantly, possible interactions between the linker loop (LL) region of SNAP25 and VAMP2 are observed, suggests that the LL region may facilitate VAMP2 binding and SNARE initiation. The forced unzipping of SNARE in the presence of membrane and LL of SNAP25 reveals the possible pathway for energy generation of SNARE zipping, provides information to understand how force may regulate the cooperativity between the membrane and the SNARE complex.

Interaction of the Complexin Accessory Helix with Synaptobrevin Regulates Spontaneous Fusion

Neuronal transmitters are released from nerve terminals via the fusion of synaptic vesicles with the plasma membrane. Vesicles attach to membranes via a specialized protein machinery composed of membrane-attached (t-SNARE) and vesicle-attached (v-SNARE) proteins that zipper together to form a coiled-coil SNARE bundle that brings the two fusing membranes into close proximity. Neurotransmitter release may occur either in response to an action potential or through spontaneous fusion. A cytosolic protein, Complexin (Cpx), binds the SNARE complex and restricts spontaneous exocytosis by acting as a fusion clamp. We previously proposed a model in which the interaction between Cpx and the v-SNARE serves as a spring to prevent premature zippering of the SNARE complex, thereby reducing the likelihood of fusion. To test this model, we combined molecular-dynamics (MD) simulations and site-directed mutagenesis of Cpx and SNAREs in Drosophila. MD simulations of the Drosophila Cpx-SNARE complex demonstrated that Cpx's interaction with the v-SNARE promotes unraveling of the v-SNARE off the core SNARE bundle. We investigated clamping properties in the syx3-69 paralytic mutant, which has a single-point mutation in the t-SNARE and displays enhanced spontaneous release. MD simulations demonstrated an altered interaction of Cpx with the SNARE bundle that hindered v-SNARE unraveling by Cpx, thus compromising clamping. We used our model to predict mutations that should enhance the ability of Cpx to prevent full assembly of the SNARE complex. MD simulations predicted that a weakened interaction between the Cpx accessory helix and the v-SNARE would enhance Cpx flexibility and thus promote separation of SNAREs, reducing spontaneous fusion. We generated transgenic Drosophila with mutations in Cpx and the v-SNARE that disrupted a salt bridge between these two proteins. As predicted, both lines demonstrated a selective inhibition in spontaneous release, suggesting that Cpx acts as a fusion clamp that restricts full SNARE zippering.

Coarse-Grained Model of SNARE-Mediated Docking

Synaptic transmission requires that vesicles filled with neurotransmitter molecules be docked to the plasma membrane by the SNARE protein complex. The SNARE complex applies attractive forces to overcome the long-range repulsion between the vesicle and membrane. To understand how the balance between the attractive and repulsive forces defines the equilibrium docked state we have developed a model that combines the mechanics of vesicle/membrane deformation with an apparently new coarse-grained model of the SNARE complex. The coarse-grained model of the SNARE complex is calibrated by comparison with all-atom molecular dynamics simulations as well as by force measurements in laser tweezer experiments. The model for vesicle/membrane interactions includes the forces produced by membrane deformation and hydration or electrostatic repulsion. Combining these two parts, the coarse-grained model of the SNARE complex with membrane mechanics, we study how the equilibrium docked state varies with the number of SNARE complexes. We find that a single SNARE complex is able to bring a typical synaptic vesicle to within a distance of ∼ 3 nm from the membrane. Further addition of SNARE complexes shortens this distance, but an overdocked state of >4-6 SNAREs actually increases the equilibrium distance.

Role of Internal Water on Protein Thermal Stability: The Case of Homologous G Domains

In this work, we address the question of whether the enhanced stability of thermophilic proteins has a direct connection with internal hydration. Our model systems are two homologous G domains of different stability: the mesophilic G domain of the elongation factor thermal unstable protein from E. coli and the hyperthermophilic G domain of the EF-1α protein from S. solfataricus. Using molecular dynamics simulation at the microsecond time scale, we show that both proteins host water molecules in internal cavities and that these molecules exchange with the external solution in the nanosecond time scale. The hydration free energy of these sites evaluated via extensive calculations is found to be favorable for both systems, with the hyperthermophilic protein offering a slightly more favorable environment to host water molecules. We estimate that, under ambient conditions, the free energy gain due to internal hydration is about 1.3 kcal/mol in favor of the hyperthermophilic variant. However, we also find that, at the high working temperature of the hyperthermophile, the cavities are rather dehydrated, meaning that under extreme conditions other molecular factors secure the stability of the protein. Interestingly, we detect a clear correlation between the hydration of internal cavities and the protein conformational landscape. The emerging picture is that internal hydration is an effective observable to probe the conformational landscape of proteins. In the specific context of our investigation, the analysis confirms that the hyperthermophilic G domain is characterized by multiple states and it has a more flexible structure than its mesophilic homologue.

Interaction of the complexin accessory helix with the C-terminus of the SNARE complex: molecular-dynamics model of the fusion clamp

SNARE complexes form between the synaptic vesicle protein synaptobrevin and the plasma membrane proteins syntaxin and SNAP25 to drive membrane fusion. A cytosolic protein, complexin (Cpx), binds to the SNARE bundle, and its accessory helix (AH) functions to clamp synaptic vesicle fusion. We performed molecular-dynamics simulations of the SNARE/Cpx complex and discovered that at equilibrium the Cpx AH forms tight links with both synaptobrevin and SNAP25. To simulate the effect of electrostatic repulsion between vesicle and membrane on the SNARE complex, we calculated the electrostatic force and performed simulations with an external force applied to synaptobrevin. We found that the partially unzipped state of the SNARE bundle can be stabilized by interactions with the Cpx AH, suggesting a simple mechanistic explanation for the role of Cpx in fusion clamping. To test this model, we performed experimental and computational characterizations of the syx(3-69)Drosophila mutant, which has a point mutation in syntaxin that causes increased spontaneous fusion. We found that this mutation disrupts the interaction of the Cpx AH with synaptobrevin, partially imitating the cpx null phenotype. Our results support a model in which the Cpx AH clamps fusion by binding to the synaptobrevin C-terminus, thus preventing full SNARE zippering.

How SNARE molecules mediate membrane fusion: recent insights from molecular simulations

SNARE molecules are the core constituents of the protein machinery that facilitate fusion of synaptic vesicles with the presynaptic plasma membrane, resulting in the release of neurotransmitter. On a molecular level, SNARE complexes seem to play a quite versatile and involved role during all stages of fusion. In addition to merely triggering fusion by forcing the opposing membranes into close proximity, SNARE complexes are now seen to also overcome subsequent fusion barriers and to actively guide the fusion reaction up to the expansion of the fusion pore. Here, we review recent advances in the understanding of SNARE-mediated membrane fusion by molecular simulations.

GPU-accelerated atom and dynamic bond visualization using hyperballs: a unified algorithm for balls, sticks, and hyperboloids

Ray casting on graphics processing units (GPUs) opens new possibilities for molecular visualization. We describe the implementation and calculation of diverse molecular representations such as licorice, ball-and-stick, space-filling van der Waals spheres, and approximated solvent-accessible surfaces using GPUs. We introduce HyperBalls, an improved ball-and-stick representation replacing tubes, linking the atom spheres by hyperboloids that can smoothly connect them. This type of depiction is particularly useful to represent dynamic phenomena, such as the evolution of noncovalent bonds. It is furthermore well suited to represent coarse-grained models and spring networks. All these representations can be defined by a single general algebraic equation that is adapted for the ray-casting technique and is well suited for execution on the GPU. Using GPU capabilities, this implementation can routinely, accurately, and interactively render molecules ranging from a few atoms up to huge macromolecular assemblies with more than 500,000 particles. In simple cases, based only on spheres, we have been able to display up to two million atoms smoothly.

A unique and potent protein binding nature of liposome containing polyethylenimine and polyethylene glycol: a nondisplaceable property

Most of the currently available targeting vectors are produced via the linkage of targeting molecules. However, the coupling process is complicated, and the covalent linkage may attenuate the activity of certain targeting molecules. In this study, we have developed a cationic liposome complexed with polyethylenimine and polyethylene glycol polymers (LPPC) that can capture various proteins without covalent conjugation. Characterizations of prepared LPPC revealed that the maximal-binding capacity was about 170 µg of bovine serum albumin to 40 µg of sphere-shaped LPPC (180 nm). The proteins were essentially located at or near the surface when analyzed by atomic force or transmission electron microscopy. We demonstrate that polyethylenimine was an essential component to bind the proteins. Upon the saturation of captured proteins, a given protein could not be displaced by other additional proteins and still retained its biological activity. Using a variety of functional proteins, we show some typical examples of the utility of incorporated beta-glucuronidase and antibodies onto the LPPC. The beta-glucuronidase can be used for the study of antigen-antibody interactions, whereas in studies with the antibody complex, we used anti-CD3 as an agonist to stimulate the proliferation of peripheral blood mononuclear cells via a receptor-mediated mechanism and anti-VEGFR for cell staining. In conclusion, the prepared LPPC can provide a platform to capture biologically and biochemically functional proteins on its surface for various applications, such as cell signaling, cell profiling, noncovalent enzyme-linked immunoassays, and others not mentioned.

Complex molecular assemblies at hand via interactive simulations

Studying complex molecular assemblies interactively is becoming an increasingly appealing approach to molecular modeling. Here we focus on interactive molecular dynamics (IMD) as a textbook example for interactive simulation methods. Such simulations can be useful in exploring and generating hypotheses about the structural and mechanical aspects of biomolecular interactions. For the first time, we carry out low-resolution coarse-grain IMD simulations. Such simplified modeling methods currently appear to be more suitable for interactive experiments and represent a well-balanced compromise between an important gain in computational speed versus a moderate loss in modeling accuracy compared to higher resolution all-atom simulations. This is particularly useful for initial exploration and hypothesis development for rare molecular interaction events. We evaluate which applications are currently feasible using molecular assemblies from 1900 to over 300,000 particles. Three biochemical systems are discussed: the guanylate kinase (GK) enzyme, the outer membrane protease T and the soluble N-ethylmaleimide-sensitive factor attachment protein receptors complex involved in membrane fusion. We induce large conformational changes, carry out interactive docking experiments, probe lipid-protein interactions and are able to sense the mechanical properties of a molecular model. Furthermore, such interactive simulations facilitate exploration of modeling parameters for method improvement. For the purpose of these simulations, we have developed a freely available software library called MDDriver. It uses the IMD protocol from NAMD and facilitates the implementation and application of interactive simulations. With MDDriver it becomes very easy to render any particle-based molecular simulation engine interactive. Here we use its implementation in the Gromacs software as an example.

Is assembly of the SNARE complex enough to fuel membrane fusion?

The three key players in the exocytotic release of neurotransmitters from synaptic vesicles are the SNARE (soluble N-ethylmaleimide-sensitive factor attachment protein receptor) proteins synaptobrevin 2, syntaxin 1a, and SNAP-25. Their assembly into a tight four-helix bundle complex is thought to pull the two membranes into close proximity. It is debated, however, whether the energy generated suffices for membrane fusion. Here, we have determined the thermodynamic properties of the individual SNARE assembly steps by isothermal titration calorimetry. We found extremely large favorable enthalpy changes counterbalanced by positive entropy changes, reflecting the major conformational changes upon assembly. To circumvent the fact that ternary complex formation is essentially irreversible, we used a stabilized syntaxin-SNAP-25 heterodimer to study synaptobrevin binding. This strategy revealed that the N-terminal synaptobrevin coil binds reversibly with nanomolar affinity. This suggests that individual, membrane-bridging SNARE complexes can provide much less pulling force than previously claimed.