Differential Membrane Binding Mechanics of Synaptotagmin Isoforms Observed in Atomic Detail

Improved Highly Mobile Membrane Mimetic Model for Investigating Protein-Cholesterol Interactions

Cholesterol (CHL) plays an integral role in modulating the function and activity of various mammalian membrane proteins. Due to the slow dynamics of lipids, conventional computational studies of protein-CHL interactions rely on either long-time scale atomistic simulations or coarse-grained approximations to sample the process. A highly mobile membrane mimetic (HMMM) has been developed to enhance lipid diffusion and thus used to facilitate the investigation of lipid interactions with peripheral membrane proteins and, with customized in silico solvents to replace phospholipid tails, with integral membrane proteins. Here, we report an updated HMMM model that is able to include CHL, a nonphospholipid component of the membrane, henceforth called HMMM-CHL. To this end, we had to optimize the effect of the customized solvents on CHL behavior in the membrane. Furthermore, the new solvent is compatible with simulations using force-based switching protocols. In the HMMM-CHL, both improved CHL dynamics and accelerated lipid diffusion are integrated. To test the updated model, we have applied it to the characterization of protein-CHL interactions in two membrane protein systems, the human β2-adrenergic receptor (β2AR) and the mitochondrial voltage-dependent anion channel 1 (VDAC-1). Our HMMM-CHL simulations successfully identified CHL binding sites and captured detailed CHL interactions in excellent consistency with experimental data as well as other simulation results, indicating the utility of the improved model in applications where an enhanced sampling of protein-CHL interactions is desired.

Atomistic Characterization of Beta-2-Glycoprotein I Domain V Interaction with Anionic Membranes

Background

Interaction of beta-2-glycoprotein I ( β 2 GPI) with anionic membranes is crucial in antiphospholipid syndrome (APS), implicating the role of it's membrane bind-ing domain, Domain V (DV). The mechanism of DV binding to anionic lipids is not fully understood.

Objectives

This study aims to elucidate the mechanism by which DV of β 2 GPI binds to anionic membranes.

Methods

We utilized molecular dynamics (MD) simulations to investigate the struc-tural basis of anionic lipid recognition by DV. To corroborate the membrane-binding mode identified in the HMMM simulations, we conducted additional simulations using a full mem-brane model.

Results

The study identified critical regions in DV, namely the lysine-rich loop and the hydrophobic loop, essential for membrane association via electrostatic and hydrophobic interactions, respectively. A novel lysine pair contributing to membrane binding was also discovered, providing new insights into β 2 GPI's membrane interaction. Simulations revealed two distinct binding modes of DV to the membrane, with mode 1 characterized by the insertion of the hydrophobic loop into the lipid bilayer, suggesting a dominant mechanism for membrane association. This interaction is pivotal for the pathogenesis of APS, as it facilitates the recognition of β 2 GPI by antiphospholipid antibodies.

Conclusion

The study advances our understanding of the molecular interactions be-tween β 2 GPI's DV and anionic membranes, crucial for APS pathogenesis. It highlights the importance of specific regions in DV for membrane binding and reveals a predominant bind-ing mode. These findings have significant implications for APS diagnostics and therapeutics, offering a deeper insight into the molecular basis of the syndrome.

Synaptotagmin 7 C2 domains induce membrane curvature stress via electrostatic interactions and the wedge mechanism

Synaptotagmin 7 (Syt-7) is part of the synaptotagmin protein family that regulates exocytotic lipid membrane fusion. Among the family, Syt-7 stands out by its membrane binding strength and stabilization of long-lived membrane fusion pores. Given that Syt-7 vesicles form long-lived fusion pores, we hypothesize that its interactions with the membrane stabilize the specific curvatures, thicknesses, and lipid compositions that support a metastable fusion pore. Using all-atom molecular dynamics simulations and FRET-based assays of Syt-7's membrane-binding C2 domains (C2A and C2B), we found that Syt-7 C2 domains sequester anionic lipids, are sensitive to cholesterol, thin membranes, and generate lipid membrane curvature by two competing, but related mechanisms. First, Syt-7 forms strong electrostatic contacts with the membrane, generating negative curvature stress. Second, Syt-7's calcium binding loops embed in the membrane surface, acting as a wedge to thin the membrane and induce positive curvature stress. These curvature mechanisms are linked by the protein insertion depth as well as the resulting protein tilt. Simplified quantitative models of the curvature-generating mechanisms link simulation observables to their membrane-reshaping effectiveness.

Characterization of promoter elements of isoprene-responsive genes and the ability of isoprene to bind START domain transcription factors

Isoprene has recently been proposed to be a signaling molecule that can enhance tolerance of both biotic and abiotic stress. Not all plants make isoprene, but all plants tested to date respond to isoprene. We hypothesized that isoprene interacts with existing signaling pathways rather than requiring novel mechanisms for its effect on plants. We analyzed the cis-regulatory elements (CREs) in promoters of isoprene-responsive genes and the corresponding transcription factors binding these promoter elements to obtain clues about the transcription factors and other proteins involved in isoprene signaling. Promoter regions of isoprene-responsive genes were characterized using the Arabidopsis cis-regulatory element database. CREs bind ARR1, Dof, DPBF, bHLH112, GATA factors, GT-1, MYB, and WRKY transcription factors, and light-responsive elements were overrepresented in promoters of isoprene-responsive genes; CBF-, HSF-, WUS-binding motifs were underrepresented. Transcription factors corresponding to CREs overrepresented in promoters of isoprene-responsive genes were mainly those important for stress responses: drought-, salt/osmotic-, oxidative-, herbivory/wounding and pathogen-stress. More than half of the isoprene-responsive genes contained at least one binding site for TFs of the class IV (homeodomain leucine zipper) HD-ZIP family, such as GL2, ATML1, PDF2, HDG11, ATHB17. While the HD-zipper-loop-zipper (ZLZ) domain binds to the L1 box of the promoter region, a special domain called the steroidogenic acute regulatory protein-related lipid transfer, or START domain, can bind ligands such as fatty acids (e.g., linolenic and linoleic acid). We tested whether isoprene might bind in such a START domain. Molecular simulations and modeling to test interactions between isoprene and a class IV HD-ZIP family START-domain-containing protein were carried out. Without membrane penetration by the HDG11 START domain, isoprene within the lipid bilayer was inaccessible to this domain, preventing protein interactions with membrane bound isoprene. The cross-talk between isoprene-mediated signaling and other growth regulator and stress signaling pathways, in terms of common CREs and transcription factors could enhance the stability of the isoprene emission trait when it evolves in a plant but so far it has not been possible to say what how isoprene is sensed to initiate signaling responses.

Plant Terpenoid Permeability through Biological Membranes Explored via Molecular Simulations

Plants synthesize small molecule diterpenes composed of 20 carbons from precursor isopentenyl diphosphate and dimethylallyl disphosphate, manufacturing diverse compounds used for defense, signaling, and other functions. Industrially, diterpenes are used as natural aromas and flavoring, as pharmaceuticals, and as natural insecticides or repellents. Despite diterpene ubiquity in plant systems, it remains unknown how plants control diterpene localization and transport. For many other small molecules, plant cells maintain transport proteins that control compound compartmentalization. However, for most diterpene compounds, specific transport proteins have not been identified, and so it has been hypothesized that diterpenes may cross biological membranes passively. Through molecular simulation, we study membrane transport for three complex diterpenes from among the many made by members of the Lamiaceae family to determine their permeability coefficient across plasma membrane models. To facilitate accurate simulation, the intermolecular interactions for leubethanol, abietic acid, and sclareol were parametrized through the standard CHARMM methodology for incorporation into molecular simulations. To evaluate the effect of membrane composition on permeability, we simulate the three diterpenes in two membrane models derived from sorghum and yeast lipidomics data. We track permeation events within our unbiased simulations, and compare implied permeation coefficients with those calculated from Replica Exchange Umbrella Sampling calculations using the inhomogeneous solubility diffusion model. The diterpenes are observed to permeate freely through these membranes, indicating that a transport protein may not be needed to export these small molecules from plant cells. Moreover, the permeability is observed to be greater for plant-like membrane compositions when compared against animal-like membrane models. Increased permeability for diterpene molecules in plant membranes suggest that plants have tailored their membranes to facilitate low-energy transport processes for signaling molecules.

Molecular View into Preferential Binding of the Factor VII Gla Domain to Phosphatidic Acid

Factor VII (FVII) is a serine protease with a key role in initiating the coagulation cascade. It is part of a family of vitamin K-dependent clotting proteins, which require vitamin K for formation of their specialized membrane-binding domains (Gla domains). Membrane binding of the FVII Gla domain is critical to the activity of FVII, mediating the formation of its complex with other clotting factors. While Gla domains among coagulation factors are highly conserved in terms of amino acid sequence and structure, they demonstrate differential binding specificity toward anionic lipids. Although most Gla domain-containing clotting proteins display a strong preference for phosphatidylserine (PS), it has been demonstrated that FVII and protein C instead bind preferentially to phosphatidic acid (PA). We have developed the first model of the FVII Gla domain bound to PA lipids in membranes containing PA, the highly mobile membrane mimetic model, which accelerates slow diffusion of lipids in molecular dynamics simulations and therefore facilitates the membrane binding process and enhances sampling of lipid interactions. Simulations were performed using atomic level molecular dynamics, requiring a fixed charge to all atoms. The overall charge assigned to each PA lipid for this study was -1. We also developed an additional model of the FVII Gla domain bound to a 1:1 PS/PC membrane and compared the modes of binding of PS and PA lipids to FVII, allowing us to identify potential PA-specific binding sites.

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

Binding mode of SARS-CoV-2 fusion peptide to human cellular membrane

Infection of human cells by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) relies on its binding to a specific receptor and subsequent fusion of the viral and host cell membranes. The fusion peptide (FP), a short peptide segment in the spike protein, plays a central role in the initial penetration of the virus into the host cell membrane, followed by the fusion of the two membranes. Here, we use an array of molecular dynamics simulations that take advantage of the highly mobile membrane mimetic model to investigate the interaction of the SARS-CoV2 FP with a lipid bilayer representing mammalian cellular membranes at an atomic level and to characterize the membrane-bound form of the peptide. Six independent systems were generated by changing the initial positioning and orientation of the FP with respect to the membrane, and each system was simulated in five independent replicas, each for 300 ns. In 73% of the simulations, the FP reaches a stable, membrane-bound configuration, in which the peptide deeply penetrated into the membrane. Clustering of the results reveals three major membrane-binding modes (binding modes 1-3), in which binding mode 1 populates over half of the data points. Taking into account the sequence conservation among the viral FPs and the results of mutagenesis studies establishing the role of specific residues in the helical portion of the FP in membrane association, the significant depth of penetration of the whole peptide, and the dense population of the respective cluster, we propose that the most deeply inserted membrane-bound form (binding mode 1) represents more closely the biologically relevant form. Analysis of FP-lipid interactions shows the involvement of specific residues, previously described as the "fusion-active core residues," in membrane binding. Taken together, the results shed light on a key step involved in SARS-CoV2 infection, with potential implications in designing novel inhibitors.

Direct Evidence That Mutations within Dysferlin's C2A Domain Inhibit Lipid Clustering

Mechanical stress on sarcolemma can create small tears in the muscle cell membrane. Within the sarcolemma resides the multidomain dysferlin protein. Mutations in this protein render it unable to repair the sarcolemma and have been linked to muscular dystrophy. A key step in dysferlin-regulated repair is the binding of the C2A domain to the lipid membrane upon increased intracellular calcium. Mutations mapped to this domain cause loss of binding ability of the C2A domain. There is a crucial need to understand the geometry of dysferlin C2A at a membrane interface as well as cell membrane lipid reorientation when compared to that of a mutant. Here, we describe a comparison between the wild-type dysferlin C2A and a mutation to the conserved aspartic acids in the domain binding loops. To identify both the geometry and the cell membrane lipid reorientation, we applied sum frequency generation (SFG) vibrational spectroscopy and coupled it with simulated SFG spectra to observe and quantify the interaction with a model cell membrane composed of phosphotidylserine and phosphotidylcholine. Observed changes in surface pressure demonstrate that calcium-bridged electrostatic interactions govern the initial interaction of the C2A domains docking with a lipid membrane. SFG spectra taken from the amide-I region for the wild type and variant contain features near 1642, 1663, and 1675 cm^-1 related to the C2A domain β-sandwich secondary structure, indicating that the domain binds in a specific orientation. Mapping simulated SFG spectra to the experimentally collected spectra indicated that both wild-type and variant domains have nearly the same orientation to the membrane surface. However, examining the ordering of the lipids that make up a model membrane using SFG, we find that the wild type clusters the lipids as seen by the increase in the ratio of the CD₃ and CD₂ symmetric intensities by 170% for the wild type and by 120% for the variant. This study highlights the capabilities of SFG to probe with great detail biological mutations in proteins at cell membrane interfaces.

Membrane Association and Functional Mechanism of Synaptotagmin-1 in Triggering Vesicle Fusion

Upon Ca²⁺ influx, synaptic vesicles fuse with the presynaptic plasma membrane (PM) to release neurotransmitters. Membrane fusion is triggered by synaptotagmin-1, a transmembrane protein in the vesicle membrane (VM), but the mechanism is under debate. Synaptotagmin-1 contains a single transmembrane helix (TM) and two tandem C2 domains (C2A and C2B). This study aimed to use molecular dynamics simulations to elucidate how Ca²⁺-bound synaptotagmin-1, by simultaneously associating with VM and PM, brings them together for fusion. Although C2A stably associates with VM via two Ca²⁺-binding loops, C2B has a propensity to partially dissociate. Importantly, an acidic motif in the TM-C2A linker competes with VM for interacting with C2B, thereby flipping its orientation to face PM. Subsequently, C2B readily associates with PM via a polybasic cluster and a Ca²⁺-binding loop. The resulting mechanistic model for the triggering of membrane fusion by synaptotagmin-1 reconciles many experimental observations.

Membrane Interactions of Cy3 and Cy5 Fluorophores and Their Effects on Membrane-Protein Dynamics

Organic fluorophores, such as Cy3 and Cy5, have been widely used as chemical labels to probe the structure and dynamics of membrane proteins. Although a number of previous studies have reported on the possibility of some of the water-soluble fluorophores to interact with lipid bilayers, detailed fluorophore-lipid interactions and, more importantly, the potential effect of such interactions on the natural dynamics of the labeled membrane proteins have not been well studied. We have performed a large set of all-atom molecular dynamics simulations employing the highly mobile membrane mimetic model to describe spontaneous partitioning of the fluorophores into lipid bilayers with different lipid compositions. Spontaneous membrane partitioning of Cy3 and Cy5 fluorophores captured in these simulations proceeds in two steps. Electrostatic interaction between the fluorophores and the lipid headgroups facilitates the initial, fast membrane association of the fluorophores, followed by slow insertion of hydrophobic moieties into the lipid bilayer core. After the conversion of the resulting membrane-bound systems to full-membrane representations, biased-exchange umbrella sampling simulations are performed for free energy calculations, revealing a higher energy barrier for partitioning into negatively charged (phosphatidylserine or phosphatidylcholine) membranes than purely zwitterionic (phosphatidylcholine or phosphatidylethanolamine) ones. Furthermore, the potential effect of fluorophore-lipid interactions on membrane proteins has been examined by covalently linking Cy5 to single- and multipass transmembrane helical proteins. Equilibrium simulations show strong position-dependent effects of Cy5-tagging on the structure and natural dynamics of membrane proteins. Interactions between the tagged protein and Cy5 were also observed. Our results suggest that fluorophore-lipid interactions can affect the structure and dynamics of membrane proteins to various extents, especially in systems with higher structural flexibility.

The role of the C2A domain of synaptotagmin 1 in asynchronous neurotransmitter release

Following nerve stimulation, there are two distinct phases of Ca²⁺-dependent neurotransmitter release: a fast, synchronous release phase, and a prolonged, asynchronous release phase. Each of these phases is tightly regulated and mediated by distinct mechanisms. Synaptotagmin 1 is the major Ca²⁺ sensor that triggers fast, synchronous neurotransmitter release upon Ca²⁺ binding by its C₂A and C₂B domains. It has also been implicated in the inhibition of asynchronous neurotransmitter release, as blocking Ca²⁺ binding by the C₂A domain of synaptotagmin 1 results in increased asynchronous release. However, the mutation used to block Ca²⁺ binding in the previous experiments (aspartate to asparagine mutations, syt^D-N) had the unintended side effect of mimicking Ca²⁺ binding, raising the possibility that the increase in asynchronous release was directly caused by ostensibly constitutive Ca²⁺ binding. Thus, rather than modulating an asynchronous sensor, syt^D-N may be mimicking one. To directly test the C₂A inhibition hypothesis, we utilized an alternate C₂A mutation that we designed to block Ca²⁺ binding without mimicking it (an aspartate to glutamate mutation, syt^D-E). Analysis of both the original syt^D-N mutation and our alternate syt^D-E mutation at the Drosophila neuromuscular junction showed differential effects on asynchronous release, as well as on synchronous release and the frequency of spontaneous release. Importantly, we found that asynchronous release is not increased in the syt^D-E mutant. Thus, our work provides new mechanistic insight into synaptotagmin 1 function during Ca²⁺-evoked synaptic transmission and demonstrates that Ca²⁺ binding by the C₂A domain of synaptotagmin 1 does not inhibit asynchronous neurotransmitter release in vivo.

Microscopic Characterization of GRP1 PH Domain Interaction with Anionic Membranes

The pleckstrin homology (PH) domain of general receptor for phosphoionositides 1 (GRP1-PHD) binds specifically to phosphatidylinositol (3,4,5)-triphosphate (PIP3 ), and acts as a second messenger. Using an extensive array of molecular dynamics (MD) simulations employing highly mobile membrane mimetic (HMMM) model as well as complementary full membrane simulations, we capture differentiable binding and dynamics of GRP1-PHD in the presence of membranes containing PC, PS, and PIP3 lipids in varying compositions. While GRP1-PHD forms only transient interactions with pure PC membranes, incorporation of anionic lipids resulted in stable membrane-bound configurations. We report the first observation of two distinct PIP3 binding modes on GRP1-PHD, involving PIP3 interactions at a "canonical" and at an "alternate" site, suggesting the possibility of simultaneous binding of multiple anionic lipids. The full membrane simulations confirmed the stability of the membrane bound pose of GRP1-PHD as captured from our HMMM membrane binding simulations. By performing additional steered membrane unbinding simulations and calculating nonequilibrium work associated with the process, as well as metadynamics simulations, on the protein bound to full membranes, allowing for more quantitative examination of the binding strength of the GRP1-PHD to the membrane, we demonstrate that along with the bound PIP3 , surrounding anionic PS lipids increase the energetic cost of unbinding of GRP1-PHD from the canonical mode, causing them to dissociate more slowly than the alternate mode. Our results demonstrate that concurrent binding of multiple anionic lipids by GRP1-PHD contributes to its membrane affinity, which in turn control its signaling activity. © 2019 Wiley Periodicals, Inc.

The C2A domain of synaptotagmin is an essential component of the calcium sensor for synaptic transmission

The synaptic vesicle protein, synaptotagmin, is the principle Ca²⁺ sensor for synaptic transmission. Ca²⁺ influx into active nerve terminals is translated into neurotransmitter release by Ca²⁺ binding to synaptotagmin’s tandem C2 domains, triggering the fast, synchronous fusion of multiple synaptic vesicles. Two hydrophobic residues, shown to mediate Ca²⁺-dependent membrane insertion of these C2 domains, are required for this process. Previous research suggested that one of its tandem C2 domains (C₂B) is critical for fusion, while the other domain (C₂A) plays only a facilitatory role. However, the function of the two hydrophobic residues in C₂A have not been adequately tested in vivo. Here we show that these two hydrophobic residues are absolutely required for synaptotagmin to trigger vesicle fusion. Using in vivo electrophysiological recording at the Drosophila larval neuromuscular junction, we found that mutation of these two key C₂A hydrophobic residues almost completely abolished neurotransmitter release. Significantly, mutation of both hydrophobic residues resulted in more severe deficits than those seen in synaptotagmin null mutants. Thus, we report the most severe phenotype of a C₂A mutation to date, demonstrating that the C₂A domain is absolutely essential for synaptotagmin’s function as the electrostatic switch.

The postulated role of synaptotagmin’s C₂A domain in triggering neurotransmitter release has fluctuated wildly over the years. Early biochemical experiments suggested that the C₂A domain was essential, while the C₂B domain was superfluous. Then, functional experiments measuring neurotransmitter release in vivo following disruptions in Ca²⁺ binding suggested that C₂B was essential, while C₂A was superfluous. Subsequently, the use of more refined mutations to disrupt Ca²⁺ binding indicated that C₂A played a facilitatory role. Here we show two hydrophobic residues of the C₂A domain are absolutely required for synaptotagmin-triggered neurotransmitter release. Thus, after over twenty years of research, we now demonstrate that the C₂A domain of synaptotagmin is an essential component of the Ca²⁺ sensor for triggering synaptic transmission in vivo.

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.

A network of phosphatidylinositol 4,5-bisphosphate binding sites regulates gating of the Ca2+-activated Cl- channel ANO1 (TMEM16A)

ANO1 (TMEM16A) is a Ca²⁺-activated Cl^- channel that regulates diverse cellular functions including fluid secretion, neuronal excitability, and smooth muscle contraction. ANO1 is activated by elevation of cytosolic Ca²⁺ and modulated by phosphatidylinositol 4,5-bisphosphate [PI(4,5)P₂]. Here, we describe a closely concerted experimental and computational study, including electrophysiology, mutagenesis, functional assays, and extended sampling of lipid-protein interactions with molecular dynamics (MD) to characterize PI(4,5)P₂ binding modes and sites on ANO1. ANO1 currents in excised inside-out patches activated by 270 nM Ca²⁺ at +100 mV are increased by exogenous PI(4,5)P₂ with an EC₅₀ = 1.24 µM. The effect of PI(4,5)P₂ is dependent on membrane voltage and Ca²⁺ and is explained by a stabilization of the ANO1 Ca²⁺-bound open state. Unbiased atomistic MD simulations with 1.4 mol% PI(4,5)P₂ in a phosphatidylcholine bilayer identified 8 binding sites with significant probability of binding PI(4,5)P₂ Three of these sites captured 85% of all ANO1-PI(4,5)P₂ interactions. Mutagenesis of basic amino acids near the membrane-cytosol interface found 3 regions of ANO1 critical for PI(4,5)P₂ regulation that correspond to the same 3 sites identified by MD. PI(4,5)P₂ is stabilized by hydrogen bonding between amino acid side chains and phosphate/hydroxyl groups on PI(4,5)P₂ Binding of PI(4,5)P₂ alters the position of the cytoplasmic extension of TM6, which plays a crucial role in ANO1 channel gating, and increases the accessibility of the inner vestibule to Cl^- ions. We propose a model consisting of a network of 3 PI(4,5)P₂ binding sites at the cytoplasmic face of the membrane allosterically regulating ANO1 channel gating.

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.

Essay on Biomembrane Structure

Of all the macromolecular assemblies of life, the least understood is the biomembrane. This is especially true in regard to its atomic structure. Ideas on biomembranes, developed in the last 200 years, culminated in the fluid mosaic model of the membrane. In this essay, I provide a historical outline of how we arrived at our current understanding of biomembranes and the models we use to describe them. A selection of direct experimental findings on the nano-scale structure of biomembranes is taken up to discuss their physical nature, and special emphasis is put on the surprising insights that arise from atomic scale descriptions.

Phosphatidic acid induces conformational changes in Sec18 protomers that prevent SNARE priming

Eukaryotic cell homeostasis requires transfer of cellular components among organelles and relies on membrane fusion catalyzed by SNARE proteins. Inactive SNARE bundles are reactivated by hexameric N-ethylmaleimide-sensitive factor, vesicle-fusing ATPase (Sec18/NSF)-driven disassembly that enables a new round of membrane fusion. We previously found that phosphatidic acid (PA) binds Sec18 and thereby sequesters it from SNAREs and that PA dephosphorylation dissociates Sec18 from the membrane, allowing it to engage SNARE complexes. We now report that PA also induces conformational changes in Sec18 protomers and that hexameric Sec18 cannot bind PA membranes. Molecular dynamics (MD) analyses revealed that the D1 and D2 domains of Sec18 contain PA-binding sites and that the residues needed for PA binding are masked in hexameric Sec18. Importantly, these simulations also disclosed that a major conformational change occurs in the linker region between the D1 and D2 domains, which is distinct from the conformational changes that occur in hexameric Sec18 during SNARE priming. Together, these findings indicate that PA regulates Sec18 function by altering its architecture and stabilizing membrane-bound Sec18 protomers.

Activity and Thermostability of GH5 Endoglucanase Chimeras from Mesophilic and Thermophilic Parents

Cellulases from glycoside hydrolase family 5 (GH5) are key endoglucanase enzymes in the degradation of diverse polysaccharide substrates and are used in industrial enzyme cocktails to break down biomass. The GH5 family shares a canonical (βα)₈-barrel structure, where each (βα) module is essential for the enzyme's stability and activity. Despite their shared topology, the thermostability of GH5 endoglucanase enzymes can vary significantly, and highly thermostable variants are often sought for industrial applications. Based on the previously characterized thermophilic GH5 endoglucanase Egl5A from Talaromyces emersonii (TeEgl5A), which has an optimal temperature of 90°C, we created 10 hybrid enzymes with elements of the mesophilic endoglucanase Cel5 from Stegonsporium opalus (SoCel5) to determine which elements are responsible for enhanced thermostability. Five of the expressed hybrid enzymes exhibit enzyme activity. Two of these hybrids exhibited pronounced increases in the temperature optimum (10 and 20°C), the temperature at which the protein lost 50% of its activity (T ₅₀) (15 and 19°C), and the melting temperature (T_m ) (16.5 and 22.9°C) and extended half-lives (t _1/2) (∼240- and 650-fold at 55°C) relative to the values for the mesophilic parent enzyme and demonstrated improved catalytic efficiency on selected substrates. The successful hybridization strategies were validated experimentally in another GH5 endoglucanase, Cel5 from Aspergillus niger (AnCel5), which demonstrated a similar increase in thermostability. Based on molecular dynamics (MD) simulations of both the SoCel5 and TeEgl5A parent enzymes and their hybrids, we hypothesize that improved hydrophobic packing of the interface between α₂ and α₃ is the primary mechanism by which the hybrid enzymes increase their thermostability relative to that of the mesophilic parent SoCel5.IMPORTANCE Thermal stability is an essential property of enzymes in many industrial biotechnological applications, as high temperatures improve bioreactor throughput. Many protein engineering approaches, such as rational design and directed evolution, have been employed to improve the thermal properties of mesophilic enzymes. Structure-based recombination has also been used to fuse TIM barrel fragments, and even fragments from unrelated folds, to generate new structures. However, little research has been done on GH5 endoglucanases. In this study, two GH5 endoglucanases exhibiting TIM barrel structure, SoCel5 and TeEgl5A, with different thermal properties, were hybridized to study the roles of different (βα) motifs. This work illustrates the role that structure-guided recombination can play in helping to identify sequence function relationships within GH5 enzymes by supplementing natural diversity with synthetic diversity.

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.

Microscopic view of lipids and their diverse biological functions

Biological membranes and their diverse lipid constituents play key roles in a broad spectrum of cellular and physiological processes. Characterization of membrane-associated phenomena at a microscopic level is therefore essential to our fundamental understanding of such processes. Due to the semi-fluid and dynamic nature of lipid bilayers, and their complex compositions, detailed characterization of biological membranes at an atomic scale has been refractory to experimental approaches. Computational modeling and simulation offer a highly complementary toolset with sufficient spatial and temporal resolutions to fill this gap. Here, we review recent molecular dynamics studies focusing on the diversity of lipid composition of biological membranes, or aiming at the characterization of lipid-protein interaction, with the overall goal of dissecting how lipids impact biological roles of the cellular membranes.

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.

Folding a viral peptide in different membrane environments: pathway and sampling analyses

Flock House virus (FHV) is a well-characterized model system to study infection mechanisms in non-enveloped viruses. A key stage of the infection cycle is the disruption of the endosomal membrane by a component of the FHV capsid, the membrane active γ peptide. In this study, we perform all-atom molecular dynamics simulations of the 21 N-terminal residues of the γ peptide interacting with membranes of differing compositions. We carry out umbrella sampling calculations to study the folding of the peptide to a helical state in homogenous and heterogeneous membranes consisting of neutral and anionic lipids. From the trajectory data, we evaluate folding energetics and dissect the mechanism of folding in the different membrane environments. We conclude the study by analyzing the extent of configurational sampling by performing time-lagged independent component analysis.

PyContact: Rapid, Customizable, and Visual Analysis of Noncovalent Interactions in MD Simulations

Molecular dynamics (MD) simulations have become ubiquitous in all areas of life sciences. The size and model complexity of MD simulations are rapidly growing along with increasing computing power and improved algorithms. This growth has led to the production of a large amount of simulation data that need to be filtered for relevant information to address specific biomedical and biochemical questions. One of the most relevant molecular properties that can be investigated by all-atom MD simulations is the time-dependent evolution of the complex noncovalent interaction networks governing such fundamental aspects as molecular recognition, binding strength, and mechanical and structural stability. Extracting, evaluating, and visualizing noncovalent interactions is a key task in the daily work of structural biologists. We have developed PyContact, an easy-to-use, highly flexible, and intuitive graphical user interface-based application, designed to provide a toolkit to investigate biomolecular interactions in MD trajectories. PyContact is designed to facilitate this task by enabling identification of relevant noncovalent interactions in a comprehensible manner. The implementation of PyContact as a standalone application enables rapid analysis and data visualization without any additional programming requirements, and also preserves full in-program customization and extension capabilities for advanced users. The statistical analysis representation is interactively combined with full mapping of the results on the molecular system through the synergistic connection between PyContact and VMD. We showcase the capabilities and scientific significance of PyContact by analyzing and visualizing in great detail the noncovalent interactions underlying the ion permeation pathway of the human P2X3 receptor. As a second application, we examine the protein-protein interaction network of the mechanically ultrastable cohesin-dockering complex.

Molecular details of spontaneous insertion and interaction of HCV non-structure 3 protease protein domain with PIP2-containing membrane

Hepatitis C virus (HCV), known as the leading cause of liver cirrhosis, viral hepatitis, and hepatocellular carcinoma, has been affecting more than 150 million people globally. The HCV non-structure 3 (NS3) protease protein domain plays a key role in HCV replication and pathogenesis; and is currently a primary target for HCV antiviral therapy. Through unbiased molecular dynamics simulations which take advantage of the novel highly mobile membrane mimetic model, we constructed the membrane-bound state of the protein domain at the atomic level. Our results indicated that protease domain of HCV NS3 protein can spontaneously bind and penetrate to an endoplasmic reticulum complex membrane containing phosphatidylinositol 4,5-bisphosphate (PIP2). An amphipathic helix α₀ and loop S1 show their anchoring role to keep the protein on the membrane surface. Proper orientation of the protein domain at membrane surface was identified through measuring tilt angles of two specific vectors, wherein residue R161 plays a crucial role in its final orientation. Remarkably, PIP2 molecules were observed to bind to three main sites of the protease domain via specific electrostatic contacts and hydrogen bonds. PIP2-interaction determines the protein orientation at the membrane while both hydrophobic interplay and PIP2-interaction can stabilize the NS3 - membrane complex. Simulated results provide us with a detailed characterization of insertion, orientation and PIP2-interaction of NS3 protease domain at membrane environment, thus enhancing our understanding of structural functions and mechanism for the association of HCV non-structure 3 protein with respect to ER membranes.

Coupling X-Ray Reflectivity and In Silico Binding to Yield Dynamics of Membrane Recognition by Tim1

The dynamic nature of lipid membranes presents significant challenges with respect to understanding the molecular basis of protein/membrane interactions. Consequently, there is relatively little known about the structural mechanisms by which membrane-binding proteins might distinguish subtle variations in lipid membrane composition and/or structure. We have previously developed a multidisciplinary approach that combines molecular dynamics simulation with interfacial x-ray scattering experiments to produce an atomistic model for phosphatidylserine recognition by the immune receptor Tim4. However, this approach requires a previously determined protein crystal structure in a membrane-bound conformation. Tim1, a Tim4 homolog with distinct differences in both immunological function and sensitivity to membrane composition, was crystalized in a closed-loop conformation that is unlikely to support membrane binding. Here we have used a previously described highly mobile membrane mimetic membrane in combination with a conventional lipid bilayer model to generate a membrane-bound configuration of Tim1 in silico. This refined structure provided a significantly improved fit of experimental x-ray reflectivity data. Moreover, the coupling of the x-ray reflectivity analysis with both highly mobile membrane mimetic membranes and conventional lipid bilayer molecular dynamics simulations yielded a dynamic model of phosphatidylserine membrane recognition by Tim1 with atomic-level detail. In addition to providing, to our knowledge, new insights into the molecular mechanisms that distinguish the various Tim receptors, these results demonstrate that in silico membrane-binding simulations can remove the requirement that the existing crystal structure be in the membrane-bound conformation for effective x-ray reflectivity analysis. Consequently, this refined methodology has the potential for much broader applicability with respect to defining the atomistic details of membrane-binding proteins.

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.

Extension of the Highly Mobile Membrane Mimetic to Transmembrane Systems through Customized in Silico Solvents

The mechanics of the protein-lipid interactions of transmembrane proteins are difficult to capture with conventional atomic molecular dynamics, due to the slow lateral diffusion of lipids restricting sampling to states near the initial membrane configuration. The highly mobile membrane mimetic (HMMM) model accelerates lipid dynamics by modeling the acyl tails nearest to the membrane center as a fluid organic solvent while maintaining an atomic description of the lipid headgroups and short acyl tails. The HMMM has been applied to many peripheral protein systems; however, the organic solvent used to date caused deformations in transmembrane proteins by intercalating into the protein and disrupting interactions between individual side chains. We ameliorate the effect of the solvent on transmembrane protein structure through the development of two new in silico Lennard-Jones solvents. The parameters for the new solvents were determined through an extensive parameter search in order to match the bulk properties of alkanes in a highly simplified model. Using these new solvents, we substantially improve the insertion free energy profiles of 10 protein side chain analogues across the entire bilayer. In addition, we reduce the intercalation of solvent into transmembrane systems, resulting in native-like transmembrane protein structures from five different topological classes within a HMMM bilayer. The parametrization of the solvents, in addition to their computed physical properties, is discussed. By combining high lipid lateral diffusion with intact transmembrane proteins, we foresee the developed solvents being useful to efficiently identify membrane composition inhomogeneities and lipid binding caused by the presence of membrane proteins.