Interactions of phosphatase and tensin homologue (PTEN) proteins with phosphatidylinositol phosphates: insights from molecular dynamics simulations of PTEN and voltage sensitive phosphatase

Regulatory mechanisms triggered by enzyme interactions with lipid membrane surfaces

Recruitment of enzymes to intracellular membranes often modulates their catalytic activity, which can be important in cell signaling and membrane trafficking. Thus, re-localization is not only important for these enzymes to gain access to their substrates, but membrane interactions often allosterically regulate enzyme function by inducing conformational changes across different time and amplitude scales. Recent structural, biophysical and computational studies have revealed how key enzymes interact with lipid membrane surfaces, and how this membrane binding regulates protein structure and function. This review summarizes the recent progress in understanding regulatory mechanisms involved in enzyme-membrane interactions.

Molecular dynamic simulation and functional analysis of pathogenic PTEN mutations in glioblastoma

PTEN, a dual-phosphatase and scaffold protein, is one of the most commonly mutated tumour suppressor gene across various cancer types in human. The aim of this study therefore was to investigate the stability, structural and functional effects, and pathogenicity of 12 missense PTEN mutations (R15S, E18G, G36R, N49I, Y68H, I101T, C105F, D109N, V133I, C136Y, R173C and N276S) found by next generation sequencing of the PTEN gene in tissue samples obtained from glioblastoma patients. Computational tools and molecular dynamic simulation programs were used to identify the deleterious effects of these mutations. Furthermore, PTEN mRNA and protein expression levels were evaluated by qRT-PCR, Western Blot, and immunohistochemistry staining methods. Various computational tools predicted strong deleterious effects for the G36R, C105F, C136Y and N276S mutations. Molecular dynamic simulation revealed a significant decrease in protein stability for the Y68H and N276S mutations when compared with the wild type protein; whereas, C105F, D109N, V133I and R173C showed partial stability reduction. Significant residual fluctuations were observed in the R15S, N49I and C136Y mutations and radius of gyration graphs revealed the most compact structure for D109N and least for C136Y. In summary, our study is the first one to show the presence of PTEN E18G, N49I, D109N and N276S mutations in glioblastoma patients; where, D109N is neutral and N276S is a damaging and disease-associated mutation.Communicated by Ramaswamy H. Sarma.

The orchestrated signaling by PI3Kα and PTEN at the membrane interface

The oncogene PI3Kα and the tumor suppressor PTEN represent two antagonistic enzymatic activities that regulate the interconversion of the phosphoinositide lipids PI(4,5)P2 and PI(3,4,5)P3 in membranes. As such, they are defining components of phosphoinositide-based cellular signaling and membrane trafficking pathways that regulate cell survival, growth, and proliferation, and are often deregulated in cancer. In this review, we highlight aspects of PI3Kα and PTEN interplay at the intersection of signaling and membrane trafficking. We also discuss the mechanisms of PI3Kα- and PTEN- membrane interaction and catalytic activation, which are fundamental for our understanding of the structural and allosteric implications on signaling at the membrane interface and may aid current efforts in pharmacological targeting of these proteins.

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

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

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

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

•

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

•

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

•

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

•

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

The mechanism of full activation of tumor suppressor PTEN at the phosphoinositide-enriched membrane

Tumor suppressor PTEN, the second most highly mutated protein in cancer, dephosphorylates signaling lipid PIP₃ produced by PI3Ks. Excess PIP₃ promotes cell proliferation. The mechanism at the membrane of this pivotal phosphatase is unknown hindering drug discovery. Exploiting explicit solvent simulations, we tracked full-length PTEN trafficking from the cytosol to the membrane. We observed its interaction with membranes composed of zwitterionic phosphatidylcholine, anionic phosphatidylserine, and phosphoinositides, including signaling lipids PIP₂ and PIP₃. We tracked its moving away from the zwitterionic and getting absorbed onto anionic membrane that harbors PIP₃. We followed it localizing on microdomains enriched in signaling lipids, as PI3K does, and observed PIP₃ allosterically unfolding the N-terminal PIP₂ binding domain, positioning it favorably for the polybasic motif interaction with PIP₂. Finally, we determined PTEN catalytic action at the membrane, all in line with experimental observations, deciphering the mechanisms of how PTEN anchors to the membrane and restrains cancer.

•

PTEN localizes on membrane microdomains enriched in phosphoinositides, as PI3K does

•

Full PTEN activation requires both signaling lipids, PIP₂ and PIP₃

•

Strong salt bridge interactions sustain stable PTEN membrane localization

•

Substrate-induced P loop conformational change implicates PTEN catalytic activity

Advances in understanding and in multi-disciplinary methodology used to assess lipid regulation of signalling cascades from the cancer cell plasma membrane

The lipid bilayer is a functional component of cells, forming a stable platform for the initiation of key biological processes, including cell signalling. There are distinct changes in the lipid composition of cell membranes during oncogenic transformation resulting in aberrant activation and inactivation of signalling transduction pathways. Studying the role of the cell membrane in cell signalling is challenging, since techniques are often limited to by timescale, resolution, sensitivity, and averaging. To overcome these limitations, combining 'computational', 'wet-lab' and 'semi-dry' approaches offers the best opportunity to resolving complex biological processes involved in membrane organisation. In this review, we highlight analytical tools that have been applied for the study of cell signalling initiation from the cancer cell membranes through computational microscopy, biological assays, and membrane biophysics. The cancer therapeutic potential of extracellular membrane-modulating agents, such as cholesterol-reducing agents is also discussed, as is the need for future collaborative inter-disciplinary research for studying the role of the cell membrane and its components in cancer therapy.

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

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

Molecular dynamic study on PTEN frameshift mutations in breast cancer provide c2 domain as a potential biomarker

PTEN is a tumour suppressor gene known for regulating apoptosis, cell growth, and many other pathways. It is one of the most frequently mutated genes comprising the phosphatase domain (PD) and C terminal domain (C2). Direct therapeutic methods are not applicable for targeting PTEN because once gets mutated, it needs restoration. For mutant detection and restoration using PTEN mRNA there is a need to explore various mutations taking place in PTEN, identify their particular domains, and study their interactions within the cellular system. Here, we have tried to highlight a few such regions in the mutated PTEN of breast cancer patients. In this study, we have selected the top-most-occurring PTEN mutation in breast cancer and compared them to determine the specific properties of each mutation and its effect on functionality. Molecular dynamic simulation for 50 ns was performed on five structures to compare the structural behaviour of mutated PTEN in the system. Our finding suggests that frameshift mutations are more damaging and affect the c2 domain. Frameshift mutant fs_ACTT is the highest occurring as well as the most damaging mutation in all the compared structures. Docking study shows that substitution mutations D92H and R130Q causes loss of binding ability towards PIP2 in normal PTEN, interfering the dephosphorylation process. Overall, the C2 domain is more frequently mutated, and the amino acid residues in the C2 domain show more fluctuations compared to the other regions. Our study can provide the basis for selecting frequently mutated C2 domain as a potential therapeutic marker.Communicated by Ramaswamy H. Sarma.

Structural Mechanisms of PTEN Regulation

The tumor suppressor phosphatase and tensin homolog on chromosome 10 (PTEN) is a tightly regulated enzyme responsible for dephosphorylating the progrowth lipid messenger molecule phosphatidylinositol 3,4,5-trisphosphate (PIP₃) on the plasma membrane. The carboxy-terminal tail (CTT) of PTEN is key for regulation of the enzyme. When phosphorylated, the unstructured CTT interacts with the phosphatase-C2 superdomain to inactivate the enzyme by preventing membrane association. PTEN mutations associated with cancer also inactivate the enzyme. Alternate translation-initiation sites generate extended isoforms of PTEN, such as PTEN-L that has multiple roles in cells. The extended amino-terminal region bears a signal sequence and a polyarginine sequence to facilitate exit from and entry into cells, respectively, and a membrane-binding helix that activates the enzyme. This amino-terminal region also facilitates mitochondrial and nucleolar localization. This review explores PTEN structure and its impact on localization and regulation.

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.

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

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

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

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

Voltage-Sensing Phosphatases: Biophysics, Physiology, and Molecular Engineering

Voltage-sensing phosphatase (VSP) contains a voltage sensor domain (VSD) similar to that in voltage-gated ion channels, and a phosphoinositide phosphatase region similar to phosphatase and tensin homolog deleted on chromosome 10 (PTEN). The VSP gene is conserved from unicellular organisms to higher vertebrates. Membrane depolarization induces electrical driven conformational rearrangement in the VSD, which is translated into catalytic enzyme activity. Biophysical and structural characterization has revealed details of the mechanisms underlying the molecular functions of VSP. Coupling between the VSD and the enzyme is tight, such that enzyme activity is tuned in a graded fashion to the membrane voltage. Upon VSP activation, multiple species of phosphoinositides are simultaneously altered, and the profile of enzyme activity depends on the history of the membrane potential. VSPs have been the obvious candidate link between membrane potential and phosphoinositide regulation. However, patterns of voltage change regulating VSP in native cells remain largely unknown. This review addresses the current understanding of the biophysical biochemical properties of VSP and provides new insight into the proposed functions of VSP.

The hydrophobic nature of a novel membrane interface regulates the enzyme activity of a voltage-sensing phosphatase

Voltage-sensing phosphatases (VSP) contain a voltage sensor domain (VSD) similar to that of voltage-gated ion channels but lack a pore-gate domain. A VSD in a VSP regulates the cytoplasmic catalytic region (CCR). However, the mechanisms by which the VSD couples to the CCR remain elusive. Here we report a membrane interface (named ‘the hydrophobic spine’), which is essential for the coupling of the VSD and CCR. Our molecular dynamics simulations suggest that the hydrophobic spine of Ciona intestinalis VSP (Ci-VSP) provides a hinge-like motion for the CCR through the loose membrane association of the phosphatase domain. Electrophysiological experiments indicate that the voltage-dependent phosphatase activity of Ci-VSP depends on the hydrophobicity and presence of an aromatic ring in the hydrophobic spine. Analysis of conformational changes in the VSD and CCR suggests that the VSP has two states with distinct enzyme activities and that the second transition depends on the hydrophobic spine.

Interaction of monomeric Ebola VP40 protein with a plasma membrane: A coarse-grained molecular dynamics (CGMD) simulation study

Ebola virus is a lipid-enveloped filamentous virus that affects human and non-human primates and consists of several types of protein: nucleoprotein, VP30, VP35, L protein, VP40, VP24, and transmembrane glycoprotein. Among the Ebola virus proteins, its matrix protein VP40 is abundantly expressed during infection and plays a number of critical roles in oligomerization, budding and egress from the host cell. VP40 exists predominantly as a monomer at the inner leaflet of the plasma membrane, and has been suggested to interact with negatively charged lipids such as phosphatidylinositol 4,5-bisphosphate (PIP₂) and phosphatidylserine (PS) via its cationic patch. The hydrophobic loop at the C-terminal domain has also been shown to be important in the interaction between the VP40 and the membrane. However, details of the molecular mechanisms underpinning their interactions are not fully understood. This study aimed at investigating the effects of mutation in the cationic patch and hydrophobic loop on the interaction between the VP40 monomer and the plasma membrane using coarse-grained molecular dynamics simulation (CGMD). Our simulations revealed that the interaction between VP40 and the plasma membrane is mediated by the cationic patch residues. This led to the clustering of PIP₂ around the protein in the inner leaflet as a result of interactions between some cationic residues including R52, K127, K221, K224, K225, K256, K270, K274, K275 and K279 and PIP₂ lipids via electrostatic interactions. Mutation of the cationic patch or hydrophobic loop amino acids caused the protein to bind at the inner leaflet of the plasma membrane in a different orientation, where no significant clustering of PIP₂ was observed around the mutated protein. This study provides basic understanding of the interaction of the VP40 monomer and its mutants with the plasma membrane.

Beyond voltage-gated ion channels: Voltage-operated membrane proteins and cellular processes

Voltage-gated ion channels were believed to be the only voltage-sensitive proteins in excitable (and some non-excitable) cells for a long time. Emerging evidence indicates that the voltage-operated model is shared by some other transmembrane proteins expressed in both excitable and non-excitable cells. In this review, we summarize current knowledge about voltage-operated proteins, which are not classic voltage-gated ion channels as well as the voltage-dependent processes in cells for which single voltage-sensitive proteins have yet to be identified. Particularly, we will focus on the following. (1) Voltage-sensitive phosphoinositide phosphatases (VSP) with four transmembrane segments homologous to the voltage sensor domain (VSD) of voltage-gated ion channels; VSPs are the first family of proteins, other than the voltage-gated ion channels, for which there is sufficient evidence for the existence of the VSD domain; (2) Voltage-gated proton channels comprising of a single voltage-sensing domain and lacking an identified pore domain; (3) G protein coupled receptors (GPCRs) that mediate the depolarization-evoked potentiation of Ca²⁺ mobilization; (4) Plasma membrane (PM) depolarization-induced but Ca²⁺ -independent exocytosis in neurons. (5) Voltage-dependent metabolism of phosphatidylinositol 4,5-bisphosphate (PtdIns[4,5]P₂ , PIP₂ ) in the PM. These recent discoveries expand our understanding of voltage-operated processes within cellular membranes.

Computational and theoretical approaches for studies of a lipid recognition protein on biological membranes

Many cellular functions, including cell signaling and related events, are regulated by the association of peripheral membrane proteins (PMPs) with biological membranes containing anionic lipids, e.g., phosphatidylinositol phosphate (PIP). This association is often mediated by lipid recognition modules present in many PMPs. Here, I summarize computational and theoretical approaches to investigate the molecular details of the interactions and dynamics of a lipid recognition module, the pleckstrin homology (PH) domain, on biological membranes. Multiscale molecular dynamics simulations using combinations of atomistic and coarse-grained models yielded results comparable to those of actual experiments and could be used to elucidate the molecular mechanisms of the formation of protein/lipid complexes on membrane surfaces, which are often difficult to obtain using experimental techniques. Simulations revealed some modes of membrane localization and interactions of PH domains with membranes in addition to the canonical binding mode. In the last part of this review, I address the dynamics of PH domains on the membrane surface. Local PIP clusters formed around the proteins exhibit anomalous fluctuations. This dynamic change in protein-lipid interactions cause temporally fluctuating diffusivity of proteins, i.e., the short-term diffusivity of the bound protein changes substantially with time, and may in turn contribute to the formation/dissolution of protein complexes in membranes.

Domain-to-domain coupling in voltage-sensing phosphatase

Voltage-sensing phosphatase (VSP) consists of a transmembrane voltage sensor and a cytoplasmic enzyme region. The enzyme region contains the phosphatase and C2 domains, is structurally similar to the tumor suppressor phosphatase PTEN, and catalyzes the dephosphorylation of phosphoinositides. The transmembrane voltage sensor is connected to the phosphatase through a short linker region, and phosphatase activity is induced upon membrane depolarization. Although the detailed molecular characteristics of the voltage sensor domain and the enzyme region have been revealed, little is known how these two regions are coupled. In addition, it is important to know whether mechanism for coupling between the voltage sensor domain and downstream effector function is shared among other voltage sensor domain-containing proteins. Recent studies in which specific amino acid sites were genetically labeled using a fluorescent unnatural amino acid have enabled detection of the local structural changes in the cytoplasmic region of Ciona intestinalis VSP that occur with a change in membrane potential. The results of those studies provide novel insight into how the enzyme activity of the cytoplasmic region of VSP is regulated by the voltage sensor domain.

Identification of a conserved 8 aa insert in the PIP5K protein in the Saccharomycetaceae family of fungi and the molecular dynamics simulations and structural analysis to investigate its potential functional role

Homologs of the phosphatidylinositol-4-phosphate-5-kinase (PIP5K), which controls a multitude of essential cellular functions, contain a 8 aa insert in a conserved region that is specific for the Saccharomycetaceae family of fungi. Using structures of human PIP4K proteins as templates, structural models were generated of the Saccharomyces cerevisiae and human PIP5K proteins. In the modeled S. cerevisiae PIP5K, the 8 aa insert forms a surface exposed loop, present on the same face of the protein as the activation loop of the kinase domain. Electrostatic potential analysis indicates that the residues from 8 aa conserved loop form a highly positively charged surface patch, which through electrostatic interaction with the anionic portions of phospholipid head groups, is expected to play a role in the membrane interaction of the yeast PIP5K. To unravel this prediction, molecular dynamics (MD) simulations were carried out to examine the binding interaction of PIP5K, either containing or lacking the conserved signature insert, with two different membrane lipid bilayers. The results from MD studies provide insights concerning the mechanistic of interaction of PIP5K with lipid bilayer, and support the contention that the identified 8 aa conserved insert in fungal PIP5K plays an important role in the binding of this protein with membrane surface. Proteins 2017; 85:1454-1467. © 2017 Wiley Periodicals, Inc.

Molecular dynamics simulations of membrane proteins and their interactions: from nanoscale to mesoscale

Molecular dynamics simulations provide a computational tool to probe membrane proteins and systems at length scales ranging from nanometers to close to a micrometer, and on microsecond timescales. All atom and coarse-grained simulations may be used to explore in detail the interactions of membrane proteins and specific lipids, yielding predictions of lipid binding sites in good agreement with available structural data. Building on the success of protein-lipid interaction simulations, larger scale simulations reveal crowding and clustering of proteins, resulting in slow and anomalous diffusional dynamics, within realistic models of cell membranes. Current methods allow near atomic resolution simulations of small membrane organelles, and of enveloped viruses to be performed, revealing key aspects of their structure and functionally important dynamics.

The Integrin Receptor in Biologically Relevant Bilayers: Insights from Molecular Dynamics Simulations

Integrins are heterodimeric (αβ) cell surface receptors that are potential therapeutic targets for a number of diseases. Despite the existence of structural data for all parts of integrins, the structure of the complete integrin receptor is still not available. We have used available structural data to construct a model of the complete integrin receptor in complex with talin F2-F3 domain. It has been shown that the interactions of integrins with their lipid environment are crucial for their function but details of the integrin/lipid interactions remain elusive. In this study an integrin/talin complex was inserted in biologically relevant bilayers that resemble the cell plasma membrane containing zwitterionic and charged phospholipids, cholesterol and sphingolipids to study the dynamics of the integrin receptor and its effect on bilayer structure and dynamics. The results of this study demonstrate the dynamic nature of the integrin receptor and suggest that the presence of the integrin receptor alters the lipid organization between the two leaflets of the bilayer. In particular, our results suggest elevated density of cholesterol and of phosphatidylserine lipids around the integrin/talin complex and a slowing down of lipids in an annulus of ~30 Å around the protein due to interactions between the lipids and the integrin/talin F2-F3 complex. This may in part regulate the interactions of integrins with other related proteins or integrin clustering thus facilitating signal transduction across cell membranes.

Multiscale Simulations Suggest a Mechanism for the Association of the Dok7 PH Domain with PIP-Containing Membranes

Dok7 is a peripheral membrane protein that is associated with the MuSK receptor tyrosine kinase. Formation of the Dok7/MuSK/membrane complex is required for the activation of MuSK. This is a key step in the complex exchange of signals between neuron and muscle, which lead to neuromuscular junction formation, dysfunction of which is associated with congenital myasthenic syndromes. The Dok7 structure consists of a Pleckstrin Homology (PH) domain and a Phosphotyrosine Binding (PTB) domain. The mechanism of the Dok7 association with the membrane remains largely unknown. Using multi-scale molecular dynamics simulations we have explored the formation of the Dok7 PH/membrane complex. Our simulations indicate that the PH domain of Dok7 associates with membranes containing phosphatidylinositol phosphates (PIPs) via interactions of the β1/β2, β3/β4, and β5/β6 loops, which together form a positively charged surface on the PH domain and interact with the negatively charged headgroups of PIP molecules. The initial encounter of the Dok7 PH domain is followed by formation of additional interactions with the lipid bilayer, and especially with PIP molecules, which stabilizes the Dok7 PH/membrane complex. We have quantified the binding of the PH domain to the model bilayers by calculating a density landscape for protein/membrane interactions. Detailed analysis of the PH/PIP interactions reveal both a canonical and an atypical site to be occupied by the anionic lipid. PH domain binding leads to local clustering of PIP molecules in the bilayer. Association of the Dok7 PH domain with PIP lipids is therefore seen as a key step in localization of Dok7 to the membrane and formation of a complex with MuSK.

Neuromuscular junction formation and maintenance is an essential biological process, the disruption of which leads to congenital myasthenic syndromes and premature death. Dok7 is a key member in formation, maintenance and signaling in neuromuscular junctions. Dok7 is a peripheral membrane protein that is necessary for full activation of the receptor tyrosine kinase MuSK, a receptor tyrosine kinase residing in the postsynaptic membrane. The structure of Dok7 consists of both a PH domain and a PTB domain. The interaction of Dok7 with cell membranes is not well understood. Here, we use molecular simulations to show that the Dok7 PH domain preferentially binds to PIP lipid molecules when associating with a membrane. Dok7 interacts with the bilayer in both a canonical binding mode and an alternative binding mode. Our simulations also demonstrate the presence of both a canonical and an atypical binding site for PIPs in agreement with recent crystallographic studies of the ASAP1 PH domain. Analysis of density landscapes for the interaction of the Dok7 PH domain with PIP-containing lipid bilayers is also able to identify both canonical and alternative binding modes. Our improved understanding of how Dok7 interacts with a membrane is key to examining Dok7/MuSK signaling.

Interactions of Pleckstrin Homology Domains with Membranes: Adding Back the Bilayer via High-Throughput Molecular Dynamics

A molecular simulation pipeline for determining the mode of interaction of pleckstrin homology (PH) domains with phosphatidylinositol phosphate (PIP)-containing lipid bilayers is presented. We evaluate our methodology for the GRP1 PH domain via comparison with structural and biophysical data. Coarse-grained simulations yield a 2D density landscape for PH/membrane interactions alongside residue contact profiles. Predictions of the membrane localization and interactions of 13 PH domains reveal canonical, non-canonical, and dual PIP-binding sites on the proteins. Thus, the PH domains associate with the PIP molecules in the membrane via a highly positively charged loop. Some PH domains exhibit modes of interaction with PIP-containing membranes additional to this canonical binding mode. All 13 PH domains cause a degree of local clustering of PIP molecules upon binding to the membrane. This provides a global picture of PH domain interactions with membranes. The high-throughput approach could be extended to other families of peripheral membrane proteins.

Voltage-dependent motion of the catalytic region of voltage-sensing phosphatase monitored by a fluorescent amino acid

The cytoplasmic region of voltage-sensing phosphatase (VSP) derives the voltage dependence of its catalytic activity from coupling to a voltage sensor homologous to that of voltage-gated ion channels. To assess the conformational changes in the cytoplasmic region upon activation of the voltage sensor, we genetically incorporated a fluorescent unnatural amino acid, 3-(6-acetylnaphthalen-2-ylamino)-2-aminopropanoic acid (Anap), into the catalytic region of Ciona intestinalis VSP (Ci-VSP). Measurements of Anap fluorescence under voltage clamp in Xenopus oocytes revealed that the catalytic region assumes distinct conformations dependent on the degree of voltage-sensor activation. FRET analysis showed that the catalytic region remains situated beneath the plasma membrane, irrespective of the voltage level. Moreover, Anap fluorescence from a membrane-facing loop in the C2 domain showed a pattern reflecting substrate turnover. These results indicate that the voltage sensor regulates Ci-VSP catalytic activity by causing conformational changes in the entire catalytic region, without changing their distance from the plasma membrane.

Role of charged lipids in membrane structures - Insight given by simulations

Lipids and proteins are the main components of cell membranes. It is becoming increasingly clear that lipids, in addition to providing an environment for proteins to work in, are in many cases also able to modulate the structure and function of those proteins. Particularly charged lipids such as phosphatidylinositols and phosphatidylserines are involved in several examples of such effects. Molecular dynamics simulations have proved an invaluable tool in exploring these aspects. This so-called computational microscope can provide both complementing explanations for the experimental results and guide experiments to fruitful directions. In this paper, we review studies that have utilized molecular dynamics simulations to unravel the roles of charged lipids in membrane structures. We focus on lipids as active constituents of the membranes, affecting both general membrane properties as well as non-lipid membrane components, mainly proteins. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Allosteric substrate switching in a voltage-sensing lipid phosphatase

Allostery provides a critical control over enzyme activity, biasing the catalytic site between inactive and active states. We find the Ciona intestinalis voltage-sensing phosphatase (Ci-VSP), which modifies phosphoinositide signaling lipids (PIPs), to have not one but two sequential active states with distinct substrate specificities, whose occupancy is allosterically controlled by sequential conformations of the voltage sensing domain (VSD). Using fast FRET reporters of PIPs to monitor enzyme activity and voltage clamp fluorometry to monitor conformational changes in the VSD, we find that Ci-VSP switches from inactive to a PIP₃-preferring active state when the VSD undergoes an initial voltage sensing motion and then into a second PIP₂-preferring active state when the VSD activates fully. This novel 2-step allosteric control over a dual specificity enzyme enables voltage to shape PIP concentrations in time, and provides a mechanism for the complex modulation of PIP-regulated ion channels, transporters, cell motility and endo/exocytosis.

Structures of the EphA2 Receptor at the Membrane: Role of Lipid Interactions

Ephs are transmembrane receptors that mediate cell-cell signaling. The N-terminal ectodomain binds ligands and enables receptor clustering, which activates the intracellular kinase. Relatively little is known about the function of the membrane-proximal fibronectin domain 2 (FN2) of the ectodomain. Multiscale molecular dynamics simulations reveal that FN2 interacts with lipid bilayers via a site comprising K441, R443, R465, Q462, S464, S491, W467, F490, and P459-461. FN2 preferentially binds anionic lipids, a preference that is reduced in the mutant K441E + R443E. We confirm these results by measuring the binding of wild-type and mutant FN2 domains to lipid vesicles. In simulations of the complete EphA2 ectodomain plus the transmembrane region, we show that FN2 anchors the otherwise flexible ectodomain at the surface of the bilayer. Altogether, our data suggest that FN2 serves a dual function of interacting with anionic lipids and constraining the structure of the EphA2 ectodomain to adopt membrane-proximal configurations.

Biomembranes in atomistic and coarse-grained simulations

The architecture of biological membranes is tightly coupled to the localization, organization, and function of membrane proteins. The organelle-specific distribution of lipids allows for the formation of functional microdomains (also called rafts) that facilitate the segregation and aggregation of membrane proteins and thus shape their function. Molecular dynamics simulations enable to directly access the formation, structure, and dynamics of membrane microdomains at the molecular scale and the specific interactions among lipids and proteins on timescales from picoseconds to microseconds. This review focuses on the latest developments of biomembrane force fields for both atomistic and coarse-grained molecular dynamics (MD) simulations, and the different levels of coarsening of biomolecular structures. It also briefly introduces scale-bridging methods applicable to biomembrane studies, and highlights selected recent applications.

Alchembed: A Computational Method for Incorporating Multiple Proteins into Complex Lipid Geometries

A necessary step prior to starting any membrane protein computer simulation is the creation of a well-packed configuration of protein(s) and lipids. Here, we demonstrate a method, alchembed, that can simultaneously and rapidly embed multiple proteins into arrangements of lipids described using either atomistic or coarse-grained force fields. During a short simulation, the interactions between the protein(s) and lipids are gradually switched on using a soft-core van der Waals potential. We validate the method on a range of membrane proteins and determine the optimal soft-core parameters required to insert membrane proteins. Since all of the major biomolecular codes include soft-core van der Waals potentials, no additional code is required to apply this method. A tutorial is included in the Supporting Information.

Voltage-sensing phosphatase modulation by a C2 domain

The voltage-sensing phosphatase (VSP) is the first example of an enzyme controlled by changes in membrane potential. VSP has four distinct regions: the transmembrane voltage-sensing domain (VSD), the inter-domain linker, the cytosolic catalytic domain, and the C2 domain. The VSD transmits the changes in membrane potential through the inter-domain linker activating the catalytic domain which then dephosphorylates phosphatidylinositol phosphate (PIP) lipids. The role of the C2, however, has not been established. In this study, we explore two possible roles for the C2: catalysis and membrane-binding. The Ci-VSP crystal structures show that the C2 residue Y522 lines the active site suggesting a contribution to catalysis. When we mutated Y522 to phenylalanine, we found a shift in the voltage dependence of activity. This suggests hydrogen bonding as a mechanism of action. Going one step further, when we deleted the entire C2 domain, we found voltage-dependent enzyme activity was no longer detectable. This result clearly indicates the entire C2 is necessary for catalysis as well as for modulating activity. As C2s are known membrane-binding domains, we tested whether the VSP C2 interacts with the membrane. We probed a cluster of four positively charged residues lining the top of the C2 and suggested by previous studies to interact with phosphatidylinositol 4,5-bisphosphate [PI(4,5)P₂] (Kalli et al., 2014). Neutralizing those positive charges significantly shifted the voltage dependence of activity to higher voltages. We tested membrane binding by depleting PI(4,5)P₂ from the membrane using the 5HT2C receptor and found that the VSD motions as measured by voltage clamp fluorometry (VCF) were not changed. These results suggest that if the C2 domain interacts with the membrane to influence VSP function it may not occur exclusively through PI(4,5)P₂. Together, this data advances our understanding of the VSP C2 by demonstrating a necessary and critical role for the C2 domain in VSP function.

Spontaneous adsorption of coiled-coil model peptides K and E to a mixed lipid bilayer

A molecular description of the lipid-protein interactions underlying the adsorption of proteins to membranes is crucial for understanding, for example, the specificity of adsorption or the binding strength of a protein to a bilayer, or for characterizing protein-induced changes of membrane properties. In this paper, we extend an automated in silico assay (DAFT) for binding studies and apply it to characterize the adsorption of the model fusion peptides E and K to a mixed phospholipid/cholesterol membrane using coarse-grained molecular dynamics simulations. In addition, we couple the coarse-grained protocol to reverse transformation to atomistic resolution, thereby allowing to study molecular interactions with high detail. The experimentally observed differential binding of the peptides E and K to membranes, as well as the increased binding affinity of helical over unstructered peptides, could be well reproduced using the polarizable Martini coarse-grained (CG) force field. Binding to neutral membranes is shown to be dominated by initial binding of the positively charged N-terminus to the phospholipid headgroup region, followed by membrane surface-aligned insertion of the peptide at the interface between the hydrophobic core of the membrane and its polar headgroup region. Both coarse-grained and atomistic simulations confirm a before hypothesized snorkeling of lysine side chains for the membrane-bound state of the peptide K. Cholesterol was found to be enriched in peptide vicinity, which is probably of importance for the mechanism of membrane fusion. The applied sequential multiscale method, using coarse-grained simulations for the slow adsorption process of peptides to membranes followed by backward transformation to atomistic detail and subsequent atomistic simulations of the preformed peptide-lipid complexes, is shown to be a versatile approach to study the interactions of peptides or proteins with biomembranes.

Unreported intrinsic disorder in proteins: Building connections to the literature on IDPs

This review opens a new series entitled “Unreported intrinsic disorder in proteins.” The goal of this series is to bring attention of researchers to an interesting phenomenon of missed (or overlooked, or ignored, or unreported) disorder. This series serves as a companion to “Digested Disorder” which provides a quarterly review of papers on intrinsically disordered proteins (IDPs) found by standard literature searches. The need for this alternative series results from the observation that there are numerous publications that describe IDPs (or hybrid proteins with ordered and disordered regions) yet fail to recognize many of the key discoveries and publications in the IDP field. By ignoring the body of work on IDPs, such publications often fail to relate their findings to prior discoveries or fail to explore the obvious implications of their work. Thus, the goal of this series is not only to review these very interesting and important papers, but also to point out how each paper relates to the IDP field and show how common tools in the IDP field can readily take the findings in new directions or provide a broader context for the reported findings.

Phospholipid-binding sites of phosphatase and tensin homolog (PTEN): exploring the mechanism of phosphatidylinositol 4,5-bisphosphate activation

The lipid phosphatase activity of the tumor suppressor phosphatase and tensin homolog (PTEN) is enhanced by the presence of its biological product, phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2). This enhancement is suggested to occur via the product binding to the N-terminal region of the protein. PTEN effects on short-chain phosphoinositide (31)P linewidths and on the full field dependence of the spin-lattice relaxation rate (measured by high resolution field cycling (31)P NMR using spin-labeled protein) are combined with enzyme kinetics with the same short-chain phospholipids to characterize where PI(4,5)P2 binds on the protein. The results are used to model a discrete site for a PI(4,5)P2 molecule close to, but distinct from, the active site of PTEN. This PI(4,5)P2 site uses Arg-47 and Lys-13 as phosphate ligands, explaining why PTEN R47G and K13E can no longer be activated by that phosphoinositide. Placing a PI(4,5)P2 near the substrate site allows for proper orientation of the enzyme on interfaces and should facilitate processive catalysis.

Interactions of peripheral proteins with model membranes as viewed by molecular dynamics simulations

Many cellular signalling and related events are triggered by the association of peripheral proteins with anionic lipids in the cell membrane (e.g. phosphatidylinositol phosphates or PIPs). This association frequently occurs via lipid-binding modules, e.g. pleckstrin homology (PH), C2 and four-point-one, ezrin, radixin, moesin (FERM) domains, present in peripheral and cytosolic proteins. Multiscale simulation approaches that combine coarse-grained and atomistic MD simulations may now be applied with confidence to investigate the molecular mechanisms of the association of peripheral proteins with model bilayers. Comparisons with experimental data indicate that such simulations can predict specific peripheral protein–lipid interactions. We discuss the application of multiscale MD simulation and related approaches to investigate the association of peripheral proteins which contain PH, C2 or FERM-binding modules with lipid bilayers of differing phospholipid composition, including bilayers containing multiple PIP molecules.

Sidekick for Membrane Simulations: Automated Ensemble Molecular Dynamics Simulations of Transmembrane Helices

The interactions of transmembrane (TM) α-helices with the phospholipid membrane and with one another are central to understanding the structure and stability of integral membrane proteins. These interactions may be analyzed via coarse grained molecular dynamics (CGMD) simulations. To obtain statistically meaningful analysis of TM helix interactions, large (N ca. 100) ensembles of CGMD simulations are needed. To facilitate the running and analysis of such ensembles of simulations, we have developed Sidekick, an automated pipeline software for performing high throughput CGMD simulations of α-helical peptides in lipid bilayer membranes. Through an end-to-end approach, which takes as input a helix sequence and outputs analytical metrics derived from CGMD simulations, we are able to predict the orientation and likelihood of insertion into a lipid bilayer of a given helix of a family of helix sequences. We illustrate this software via analyses of insertion into a membrane of short hydrophobic TM helices containing a single cationic arginine residue positioned at different positions along the length of the helix. From analyses of these ensembles of simulations, we estimate apparent energy barriers to insertion which are comparable to experimentally determined values. In a second application, we use CGMD simulations to examine the self-assembly of dimers of TM helices from the ErbB1 receptor tyrosine kinase and analyze the numbers of simulation repeats necessary to obtain convergence of simple descriptors of the mode of packing of the two helices within a dimer. Our approach offers a proof-of-principle platform for the further employment of automation in large ensemble CGMD simulations of membrane proteins.