Investigation of the Effect of Bilayer Composition on PKCα-C2 Domain Docking Using Molecular Dynamics Simulations

Dynamic Formation of the Protein-Lipid Pre-fusion Complex

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

Insertion Depth Modulates Protein Kinase C-δ-C1b Domain Interactions with Membrane Cholesterol as Revealed by MD Simulations

Protein kinase C delta (PKC-δ) is an important signaling molecule in human cells that has both proapoptotic as well as antiapoptotic functions. These conflicting activities can be modulated by two classes of ligands, phorbol esters and bryostatins. Phorbol esters are known tumor promoters, while bryostatins have anti-cancer properties. This is despite both ligands binding to the C1b domain of PKC-δ (δC1b) with a similar affinity. The molecular mechanism behind this discrepancy in cellular effects remains unknown. Here, we have used molecular dynamics simulations to investigate the structure and intermolecular interactions of these ligands bound to δC1b with heterogeneous membranes. We observed clear interactions between the δC1b-phorbol complex and membrane cholesterol, primarily through the backbone amide of L250 and through the K256 side-chain amine. In contrast, the δC1b-bryostatin complex did not exhibit interactions with cholesterol. Topological maps of the membrane insertion depth of the δC1b-ligand complexes suggest that insertion depth can modulate δC1b interactions with cholesterol. The lack of cholesterol interactions suggests that bryostatin-bound δC1b may not readily translocate to cholesterol-rich domains within the plasma membrane, which could significantly alter the substrate specificity of PKC-δ compared to δC1b-phorbol complexes.

Binding of DEP domain to phospholipid membranes: More than just electrostatics

Over the past decades an extensive effort has been made to provide a more comprehensive understanding of Wnt signaling, yet many regulatory and structural aspects remain elusive. Among these, the ability of Dishevelled (DVL) protein to relocalize at the plasma membrane is a crucial step in the activation of all Wnt pathways. The membrane binding of DVL was suggested to be mediated by the preferential interaction of its C-terminal DEP domain with phosphatidic acid (PA). However, due to the scarcity and fast turnover of PA, we investigated the role on the membrane association of other more abundant phospholipids. The combined results from computational simulations and experimental measurements with various model phospholipid membranes, demonstrate that the membrane binding of DEP/DVL constructs is governed by the concerted action of generic electrostatics and finely-tuned intermolecular interactions with individual lipid species. In particular, while we confirmed the strong preference for PA lipid, we also observed a weak but non-negligible affinity for phosphatidylserine, the most abundant anionic phospholipid in the plasma membrane, and phosphatidylinositol 4,5-bisphosphate. The obtained molecular insight into DEP-membrane interaction helps to elucidate the relation between changes in the local membrane composition and the spatiotemporal localization of DVL and, possibly, other DEP-containing proteins.

Analysis of the MCTP Amino Acid Sequence Reveals the Conservation of Putative Calcium- and Lipid-Binding Pockets Within the C2 Domains In Silico

MCTPs (Multiple C2 Domains and Transmembrane region Proteins) are evolutionarily and structurally related to other C2 proteins, which are central to exocytosis and membrane trafficking; however, their specific function has been little studied. MCTPs are associated with endosomes and the endoplasmic reticulum and possess three C2 domains (C2A-C2C) and two transmembrane regions (TMRs) well conserved in different species. Here, we generated structural models of the MCTP C2 domains of C. elegans and analyzed their putative function by docking, which revealed that these domains possess Ca²⁺- and lipid-binding pockets, suggesting that MCTPs play a significant, calcium-dependent role in membrane physiology.

Multifaceted Analyses of Isolated Mitochondria Establish the Anticancer Drug 2-Hydroxyoleic Acid as an Inhibitor of Substrate Oxidation and an Activator of Complex IV-Dependent State 3 Respiration

The synthetic fatty acid 2-hydroxyoleic acid (2OHOA) has been extensively investigated as a cancer therapy mainly based on its regulation of membrane lipid composition and structure, activating various cell fate pathways. We discovered, additionally, that 2OHOA can uncouple oxidative phosphorylation, but this has never been demonstrated mechanistically. Here, we explored the effect of 2OHOA on mitochondria isolated by ultracentrifugation from U118MG glioblastoma cells. Mitochondria were analyzed by shotgun lipidomics, molecular dynamic simulations, spectrophotometric assays for determining respiratory complex activity, mass spectrometry for assessing beta oxidation and Seahorse technology for bioenergetic profiling. We showed that the main impact of 2OHOA on mitochondrial lipids is their hydroxylation, demonstrated by simulations to decrease co-enzyme Q diffusion in the liquid disordered membranes embedding respiratory complexes. This decreased co-enzyme Q diffusion can explain the inhibition of disjointly measured complexes I–III activity. However, it doesn’t explain how 2OHOA increases complex IV and state 3 respiration in intact mitochondria. This increased respiration probably allows mitochondrial oxidative phosphorylation to maintain ATP production against the 2OHOA-mediated inhibition of glycolytic ATP production. This work correlates 2OHOA function with its modulation of mitochondrial lipid composition, reflecting both 2OHOA anticancer activity and adaptation to it by enhancement of state 3 respiration.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

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

The Accessory Helix of Complexin Stabilizes a Partially Unzippered State of the SNARE Complex and Mediates the Complexin Clamping Function In Vivo

Spontaneous synaptic transmission is regulated by the protein complexin (Cpx). Cpx binds the SNARE complex, a coil-coiled four-helical bundle that mediates the attachment of a synaptic vesicle (SV) to the presynaptic membrane (PM). Cpx is thought to clamp spontaneous fusion events by stabilizing a partially unraveled state of the SNARE bundle; however, the molecular detail of this mechanism is still debated. We combined electrophysiology, molecular modeling, and site-directed mutagenesis in Drosophila to develop and validate the atomic model of the Cpx-mediated clamped state of the SNARE complex. We took advantage of botulinum neurotoxins (BoNTs) B and G, which cleave the SNARE protein synaptobrevin (Syb) at different sites. Monitoring synaptic depression on BoNT loading revealed that the clamped state of the SNARE complex has two or three unraveled helical turns of Syb. Site-directed mutagenesis showed that the Cpx clamping function is predominantly maintained by its accessory helix (AH), while molecular modeling suggested that the Cpx AH interacts with the unraveled C terminus of Syb and the SV lipid bilayer. The developed molecular model was employed to design new Cpx poor-clamp and super-clamp mutations and to tested the predictions in silico employing molecular dynamics simulations. Subsequently, we generated Drosophila lines harboring these mutations and confirmed the poor-clamp and super-clamp phenotypes in vivo. Altogether, these results validate the atomic model of the Cpx-mediated fusion clamp, wherein the Cpx AH inserts between the SNARE bundle and the SV lipid bilayer, simultaneously binding the unraveled C terminus of Syb and preventing full SNARE assembly.

SNARE complex alters the interactions of the Ca2+ sensor synaptotagmin 1 with lipid bilayers

Release of neuronal transmitters from nerve terminals is triggered by the molecular Ca²⁺ sensor synaptotagmin 1 (Syt1). Syt1 is a transmembrane protein attached to the synaptic vesicle (SV), and its cytosolic region comprises two domains, C2A and C2B, which are thought to penetrate into lipid bilayers upon Ca²⁺ binding. Before fusion, SVs become attached to the presynaptic membrane (PM) by the four-helical SNARE complex, which is thought to bind the C2B domain in vivo. To understand how the interactions of Syt1 with lipid bilayers and the SNARE complex trigger fusion, we performed molecular dynamics (MD) simulations at a microsecond scale. We investigated how the isolated C2 modules and the C2AB tandem of Syt1 interact with membranes mimicking either SV or PM. The simulations showed that the C2AB tandem can either bridge SV and PM or insert into PM with its Ca²⁺-bound tips and that the latter configuration is more favorable. Surprisingly, C2 domains did not cooperate in penetrating into PM but instead mutually hindered their insertion into the bilayer. To test whether the interaction of Syt1 with lipid bilayers could be affected by the C2B-SNARE attachment, we performed systematic conformational analysis of the C2AB-SNARE complex. Notably, we found that the C2B-SNARE interface precludes the coupling of C2 domains and promotes their insertion into PM. We performed the MD simulations of the prefusion protein complex positioned between the lipid bilayers mimicking PM and SV, and our results demonstrated in silico that the presence of the Ca²⁺ bound C2AB tandem promotes lipid merging. Altogether, our MD simulations elucidated the role of the Syt1-SNARE interactions in the fusion process and produced the dynamic all-atom model of the prefusion protein-lipid complex.

The Implications for Cells of the Lipid Switches Driven by Protein-Membrane Interactions and the Development of Membrane Lipid Therapy

The cell membrane contains a variety of receptors that interact with signaling molecules. However, agonist-receptor interactions not always activate a signaling cascade. Amphitropic membrane proteins are required for signal propagation upon ligand-induced receptor activation. These proteins localize to the plasma membrane or internal compartments; however, they are only activated by ligand-receptor complexes when both come into physical contact in membranes. These interactions enable signal propagation. Thus, signals may not propagate into the cell if peripheral proteins do not co-localize with receptors even in the presence of messengers. As the translocation of an amphitropic protein greatly depends on the membrane's lipid composition, regulation of the lipid bilayer emerges as a novel therapeutic strategy. Some of the signals controlled by proteins non-permanently bound to membranes produce dramatic changes in the cell's physiology. Indeed, changes in membrane lipids induce translocation of dozens of peripheral signaling proteins from or to the plasma membrane, which controls how cells behave. We called these changes "lipid switches", as they alter the cell's status (e.g., proliferation, differentiation, death, etc.) in response to the modulation of membrane lipids. Indeed, this discovery enables therapeutic interventions that modify the bilayer's lipids, an approach known as membrane-lipid therapy (MLT) or melitherapy.

The energetics of protein-lipid interactions as viewed by molecular simulations

Membranes are formed from a bilayer containing diverse lipid species with which membrane proteins interact. Integral, membrane proteins are embedded in this bilayer, where they interact with lipids from their surroundings, whilst peripheral membrane proteins bind to lipids at the surface of membranes. Lipid interactions can influence the function of membrane proteins, either directly or allosterically. Both experimental (structural) and computational approaches can reveal lipid binding sites on membrane proteins. It is, therefore, important to understand the free energies of these interactions. This affords a more complete view of the engagement of a particular protein with the biological membrane surrounding it. Here, we describe many computational approaches currently in use for this purpose, including recent advances using both free energy and unbiased simulation methods. In particular, we focus on interactions of integral membrane proteins with cholesterol, and with anionic lipids such as phosphatidylinositol 4,5-bis-phosphate and cardiolipin. Peripheral membrane proteins are exemplified via interactions of PH domains with phosphoinositide-containing membranes. We summarise the current state of the field and provide an outlook on likely future directions of investigation.

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.

Lateral distribution of phosphatidylinositol 4,5-bisphosphate in membranes regulates formin- and ARP2/3-mediated actin nucleation

Spatial and temporal control of actin polymerization is fundamental for many cellular processes, including cell migration, division, vesicle trafficking, and response to agonists. Many actin-regulatory proteins interact with phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂) and are either activated or inactivated by local PI(4,5)P₂ concentrations that form transiently at the cytoplasmic face of cell membranes. The molecular mechanisms of these interactions and how the dozens of PI(4,5)P₂-sensitive actin-binding proteins are selectively recruited to membrane PI(4,5)P₂ pools remains undefined. Using a combination of biochemical, imaging, and cell biologic studies, combined with molecular dynamics and analytical theory, we test the hypothesis that the lateral distribution of PI(4,5)P₂ within lipid membranes and native plasma membranes alters the capacity of PI(4,5)P₂ to nucleate actin assembly in brain and neutrophil extracts and show that activities of formins and the Arp2/3 complex respond to PI(4,5)P₂ lateral distribution. Simulations and analytical theory show that cholesterol promotes the cooperative interaction of formins with multiple PI(4,5)P₂ headgroups in the membrane to initiate actin nucleation. Masking PI(4,5)P₂ with neomycin or disrupting PI(4,5)P₂ domains in the plasma membrane by removing cholesterol decreases the ability of these membranes to nucleate actin assembly in cytoplasmic extracts.

Polyphosphoinositide-Binding Domains: Insights from Peripheral Membrane and Lipid-Transfer Proteins

Within eukaryotic cells, biochemical reactions need to be organized on the surface of membrane compartments that use distinct lipid constituents to dynamically modulate the functions of integral proteins or influence the selective recruitment of peripheral membrane effectors. As a result of these complex interactions, a variety of human pathologies can be traced back to improper communication between proteins and membrane surfaces; either due to mutations that directly alter protein structure or as a result of changes in membrane lipid composition. Among the known structural lipids found in cellular membranes, phosphatidylinositol (PtdIns) is unique in that it also serves as the membrane-anchored precursor of low-abundance regulatory lipids, the polyphosphoinositides (PPIn), which have restricted distributions within specific subcellular compartments. The ability of PPIn lipids to function as signaling platforms relies on both non-specific electrostatic interactions and the selective stereospecific recognition of PPIn headgroups by specialized protein folds. In this chapter, we will attempt to summarize the structural diversity of modular PPIn-interacting domains that facilitate the reversible recruitment and conformational regulation of peripheral membrane proteins. Outside of protein folds capable of capturing PPIn headgroups at the membrane interface, recent studies detailing the selective binding and bilayer extraction of PPIn species by unique functional domains within specific families of lipid-transfer proteins will also be highlighted. Overall, this overview will help to outline the fundamental physiochemical mechanisms that facilitate localized interactions between PPIn lipids and the wide-variety of PPIn-binding proteins that are essential for the coordinate regulation of cellular metabolism and membrane dynamics.

Phosphatidylinositol (4, 5)-bisphosphate targets double C2 domain protein B to the plasma membrane

Double C2 domain protein B (DOC2B) is a high-affinity Ca²⁺ sensor that translocates from the cytosol to the plasma membrane (PM) and promotes vesicle priming and fusion. However, the molecular mechanism underlying its translocation and targeting to the PM in living cells is not completely understood. DOC2B interacts in vitro with the PM components phosphatidylserine, phosphatidylinositol (4, 5)-bisphosphate [PI(4, 5)P₂ ] and target SNAREs (t-SNAREs). Here, we show that PI(4, 5)P₂ hydrolysis at the PM of living cells abolishes DOC2B translocation, whereas manipulations of t-SNAREs and other phosphoinositides have no effect. Moreover, we were able to redirect DOC2B to intracellular membranes by synthesizing PI(4, 5)P₂ in those membranes. Molecular dynamics simulations and mutagenesis in the calcium and PI(4, 5)P₂ -binding sites strengthened our findings, demonstrating that both calcium and PI(4, 5)P₂ are required for the DOC2B-PM association and revealing multiple PI(4, 5)P₂ -C2B interactions. In addition, we show that DOC2B translocation to the PM is ATP-independent and occurs in a diffusion-like manner. Our data suggest that the Ca²⁺ -triggered translocation of DOC2B is diffusion-driven and aimed at PI(4, 5)P₂ -containing membranes.