Modes of Interaction of Pleckstrin Homology Domains with Membranes: Toward a Computational Biochemistry of Membrane Recognition

GōMartini 3: From large conformational changes in proteins to environmental bias corrections

Coarse-grained modeling has become an important tool to supplement experimental measurements, allowing access to spatio-temporal scales beyond all-atom based approaches. The GōMartini model combines structure- and physics-based coarse-grained approaches, balancing computational efficiency and accurate representation of protein dynamics with the capabilities of studying proteins in different biological environments. This paper introduces an enhanced GōMartini model, which combines a virtual-site implementation of Gō models with Martini 3. The implementation has been extensively tested by the community since the release of the reparametrized version of Martini. This work demonstrates the capabilities of the model in diverse case studies, ranging from protein-membrane binding to protein-ligand interactions and AFM force profile calculations. The model is also versatile, as it can address recent inaccuracies reported in the Martini protein model. Lastly, the paper discusses the advantages, limitations, and future perspectives of the Martini 3 protein model and its combination with Gō models.

The Influence of Phosphoinositide Lipids in the Molecular Biology of Membrane Proteins: Recent Insights from Simulations

The phosphoinositide family of membrane lipids play diverse and critical roles in eukaryotic molecular biology. Much of this biological activity derives from interactions of phosphoinositide lipids with integral and peripheral membrane proteins, leading to modulation of protein structure, function, and cellular distribution. Since the discovery of phosphoinositides in the 1940s, combined molecular biology, biophysical, and structural approaches have made enormous progress in untangling this vast and diverse cellular network of interactions. More recently, in silico approaches such as molecular dynamics simulations have proven to be an asset in prospectively identifying, characterising, explaining the structural basis of these interactions, and in the best cases providing atomic level testable hypotheses on how such interactions control the function of a given membrane protein. This review details a number of recent seminal discoveries in phosphoinositide biology, enabled by advanced biomolecular simulation, and its integration with molecular biology, biophysical, and structural biology approaches. The results of the simulation studies agree well with experimental work, and in a number of notable cases have arrived at the key conclusion several years in advance of the experimental structures. SUMMARY: Hedger and Yen review developments in simulations of phosphoinositides and membrane proteins.

Lipid transport protein ORP2A promotes glucose signaling by facilitating RGS1 degradation

Heterotrimeric GTP-binding proteins (G proteins) are a group of regulators essential for signal transmission into cells. Regulator of G protein signaling 1 (AtRGS1) possesses intrinsic GTPase-accelerating protein (GAP) activity and could suppress G protein and glucose signal transduction in Arabidopsis (Arabidopsis thaliana). However, how AtRGS1 activity is regulated is poorly understood. Here, we identified a knockout mutant of oxysterol binding protein-related protein 2A, orp2a-1, which exhibits similar phenotypes to the arabidopsis g-protein beta 1-2 (agb1-2) mutant. Transgenic lines overexpressing ORP2A displayed short hypocotyls, a hypersensitive response to sugar, and lower intracellular AtRGS1 levels than the control. Consistently, ORP2A interacted with AtRGS1 in vitro and in vivo. Tissue-specific expression of 2 ORP2A alternative splicing isoforms implied functions in controlling organ size and shape. Bioinformatic data and phenotypes of orp2a-1, agb1-2, and the orp2a-1 agb1-2 double mutant revealed the genetic interactions between ORP2A and Gβ in the regulation of G protein signaling and sugar response. Both alternative protein isoforms of ORP2A localized in the endoplasmic reticulum (ER), plasma membrane (PM), and ER-PM contact sites and interacted with vesicle-associated membrane protein-associated protein 27-1 (VAP27-1) in vivo and in vitro through their two phenylalanines in an acidic track-like motif. ORP2A also displayed differential phosphatidyl phosphoinositide binding activity mediated by the pleckstrin homology domain in vitro. Taken together, the Arabidopsis membrane protein ORP2A interacts with AtRGS1 and VAP27-1 to positively regulate G protein and sugar signaling by facilitating AtRGS1 degradation.

Intrinsically disordered region of talin's FERM domain functions as an initial PIP2 recognition site

Focal adhesions (FAs) mediate the interaction of the cytoskeleton with the extracellular matrix in a highly dynamic fashion. Talin is a central regulator, adaptor protein, and mechano-sensor of FA complexes. For recruitment and firm attachment at FAs, talin's N-terminal FERM domain binds to phosphatidylinositol 4,5-bisphosphate (PIP2)-enriched membranes. A newly published autoinhibitory structure of talin-1, where the known PIP2 interaction sites are covered up, lead us to hypothesize that a hitherto less examined loop insertion of the FERM domain acts as an additional and initial site of contact. We evaluated direct interactions of talin-1 with a PIP2 membrane by means of atomistic molecular dynamics simulations. We show that this unstructured, 33-residue-long loop strongly interacts with PIP2 and can facilitate further membrane contacts, including the canonical PIP2 interactions, by serving as a flexible membrane anchor. Under force as present at FAs, the extensible FERM loop ensures talin maintains membrane contacts when pulled away from the membrane by up to 7 nm. We identify key basic residues of the anchor mediating the highly dynamic talin-membrane interaction. Our results put forward an intrinsically disordered loop as a key and highly adaptable PIP2 recognition site of talin and potentially other PIP2-binding mechano-proteins.

Affinity enhancement of CR3022 binding to RBD; in silico site directed mutagenesis using molecular dynamics simulation approaches

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a disease which caused by a novel beta coronavirus. Structural and non-structural proteins are expressed by the virus gene fragments. The RBD of the S1 protein of the virus has the ability to interact with potent antibodies including CR3022, which was characterized to target the S protein of the virus which can efficiently neutralize the SARS-CoV in vitro and in vivo. In current study, we aimed to design CR3022 based antibody with high affinity compared with wild-type CR3022 using MD simulation method. Two variants were designed based on the amino acid binding conformation and the free binding energy of the critical amino acids which involved in CR3022-RBD interactions were evaluated. In this study three complexes were evaluated; CR3022-RBD, V1-RBD and V2-RBD using molecular dynamics simulations carried out for 100 ns in each case. Then, all the complexes were simulated for 100 ns. In the next step, to calculate the free binding affinity of the wild CR3022 and mutant antibody (V1 and V2) with RBD, the PMF method was performed. The RMSD profile demonstrated that all three complexes were equilibrated after 85 ns. Furthermore, the free binding energy results indicated that the V2-RBD complex has the higher binding affinity than V1-RBD and CR3022-RBD complexes. It should be noted that in above variants, the electrostatic energy and the number of H-bonds between the antibody and RBD increased. Thus, it is suggested that both designed antibodies could be considered as appropriate candidates for covid-19 disease treatment.Communicated by Ramaswamy H. Sarma.

Standard Binding Free Energy and Membrane Desorption Mechanism for a Phospholipase C

Peripheral membrane proteins (PMPs) bind temporarily to cellular membranes and play important roles in signaling, lipid metabolism, and membrane trafficking. Obtaining accurate membrane-PMP affinities using experimental techniques is more challenging than for protein-ligand affinities in an aqueous solution. At the theoretical level, calculation of the standard protein-membrane binding free energy using molecular dynamics simulations remains a daunting challenge owing to the size of the biological objects at play, the slow lipid diffusion, and the large variation in configurational entropy that accompanies the binding process. To overcome these challenges, we used a computational framework relying on a series of potential-of-mean-force (PMF) calculations including a set of geometrical restraints on collective variables. This methodology allowed us to determine the standard binding free energy of a PMP to a phospholipid bilayer using an all-atom force field. Bacillus thuringiensis phosphatidylinositol-specific phospholipase C (BtPI-PLC) was chosen due to its importance as a virulence factor and owing to the host of experimental affinity data available. We computed a standard binding free energy of -8.2 ± 1.4 kcal/mol in reasonable agreement with the reported experimental values (-6.6 ± 0.2 kcal/mol). In light of the 2.3-μs separation PMF calculation, we investigated the mechanism whereby BtPI-PLC disengages from interactions with the lipid bilayer during separation. We describe how a short amphipathic helix engages in transitory interactions to ease the passage of its hydrophobes through the interfacial region upon desorption from the bilayer.

Molecular dynamics simulations reveal membrane lipid interactions of the full-length lymphocyte specific kinase (Lck)

The membrane-bound lymphocyte-specific protein-tyrosine kinase (Lck) triggers T cell antigen receptor signalling to initiate adaptive immune responses. Despite many structure-function studies, the mode of action of Lck and the potential role of plasma membrane lipids in regulating Lck's activity remains elusive. Advances in molecular dynamics simulations of membrane proteins in complex lipid bilayers have opened a new perspective in gathering such information. Here, we have modelled the full-length Lck open and closed conformations  using data available from different crystalographic studies and simulated its interaction with the inner leaflet of the T cell plasma membrane. In both conformations, we found that the unstructured unique domain and the structured domains including the kinase interacted with the membrane with a preference for PIP lipids. Interestingly, our simulations suggest that the Lck-SH2 domain interacts with lipids differently in the open and closed Lck conformations, demonstrating that lipid interaction can potentially regulate Lck's conformation and in turn modulate T cell signalling. Additionally, the Lck-SH2 and kinase domain residues that significantly contacted PIP lipids are found to be conserved among the Src family of kinases, thereby potentially representing similar PIP interactions within the family.

Dissecting peripheral protein-membrane interfaces

Peripheral membrane proteins (PMPs) include a wide variety of proteins that have in common to bind transiently to the chemically complex interfacial region of membranes through their interfacial binding site (IBS). In contrast to protein-protein or protein-DNA/RNA interfaces, peripheral protein-membrane interfaces are poorly characterized. We collected a dataset of PMP domains representative of the variety of PMP functions: membrane-targeting domains (Annexin, C1, C2, discoidin C2, PH, PX), enzymes (PLA, PLC/D) and lipid-transfer proteins (START). The dataset contains 1328 experimental structures and 1194 AphaFold models. We mapped the amino acid composition and structural patterns of the IBS of each protein in this dataset, and evaluated which were more likely to be found at the IBS compared to the rest of the domains' accessible surface. In agreement with earlier work we find that about two thirds of the PMPs in the dataset have protruding hydrophobes (Leu, Ile, Phe, Tyr, Trp and Met) at their IBS. The three aromatic amino acids Trp, Tyr and Phe are a hallmark of PMPs IBS regardless of whether they protrude on loops or not. This is also the case for lysines but not arginines suggesting that, unlike for Arg-rich membrane-active peptides, the less membrane-disruptive lysine is preferred in PMPs. Another striking observation was the over-representation of glycines at the IBS of PMPs compared to the rest of their surface, possibly procuring IBS loops a much-needed flexibility to insert in-between membrane lipids. The analysis of the 9 superfamilies revealed amino acid distribution patterns in agreement with their known functions and membrane-binding mechanisms. Besides revealing novel amino acids patterns at protein-membrane interfaces, our work contributes a new PMP dataset and an analysis pipeline that can be further built upon for future studies of PMPs properties, or for developing PMPs prediction tools using for example, machine learning approaches.

Two cooperative binding sites sensitize PI(4,5)P2 recognition by the tubby domain

Phosphoinositides (PIs) are lipid signaling molecules that operate by recruiting proteins to cellular membranes via PI recognition domains. The dominant PI of the plasma membrane is phosphatidylinositol 4,5-bisphosphate [PI(4,5)P₂]. One of only two PI(4,5)P₂ recognition domains characterized in detail is the tubby domain. It is essential for targeting proteins into cilia involving reversible membrane association. However, the PI(4,5)P₂ binding properties of tubby domains have remained enigmatic. Here, we used coarse-grained molecular dynamics simulations to explore PI(4,5)P₂ binding by the prototypic tubby domain. The comparatively low PI(4,5)P₂ affinity of the previously described canonical binding site is underpinned in a cooperative manner by a previously unknown, adjacent second binding site. Mutations in the previously unknown site impaired PI(4,5)P₂-dependent plasma membrane localization in living cells and PI(4,5)P₂ interaction in silico, emphasizing its importance for PI(4,5)P₂ affinity. The two-ligand binding mode may serve to sharpen the membrane association-dissociation cycle of tubby-like proteins that underlies delivery of ciliary cargo.

Systematic simulation of the interactions of pleckstrin homology domains with membranes

Pleckstrin homology (PH) domains can recruit proteins to membranes by recognition of phosphatidylinositol phosphate (PIP) lipids. Several family members are linked to diseases including cancer. We report the systematic simulation of the interactions of 100 mammalian PH domains with PIP-containing membranes. The observed PIP interaction hotspots recapitulate crystallographic binding sites and reveal a number of insights: (i) The β1 and β2 strands and their connecting loop constitute the primary PIP interaction site but are typically supplemented by interactions at the β3-β4 and β5-β6 loops; (ii) we reveal exceptional cases such as the Exoc8 PH domain; (iii) PH domains adopt different membrane-bound orientations and induce clustering of anionic lipids; and (iv) beyond family-level insights, our dataset sheds new light on individual PH domains, e.g., by providing molecular detail of secondary PIP binding sites. This work provides a global view of PH domain/membrane association involving multivalent association with anionic lipids.

Characterization and computational simulation of human Syx, a RhoGEF implicated in glioblastoma

Structural discovery of guanine nucleotide exchange factor (GEF) protein complexes is likely to become increasingly relevant with the development of new therapeutics targeting small GTPases and development of new classes of small molecules that inhibit protein‐protein interactions. Syx (also known as PLEKHG5 in humans) is a RhoA GEF implicated in the pathology of glioblastoma (GBM). Here we investigated protein expression and purification of ten different human Syx constructs and performed biophysical characterizations and computational studies that provide insights into why expression of this protein was previously intractable. We show that human Syx can be expressed and isolated and Syx is folded as observed by circular dichroism (CD) spectroscopy and actively binds to RhoA as determined by co‐elution during size exclusion chromatography (SEC). This characterization may provide critical insights into the expression and purification of other recalcitrant members of the large class of oncogenic—Diffuse B‐cell lymphoma (Dbl) homology GEF proteins. In addition, we performed detailed homology modeling and molecular dynamics simulations on the surface of a physiologically realistic membrane. These simulations reveal novel insights into GEF activity and allosteric modulation by the plekstrin homology (PH) domain. These newly revealed interactions between the GEF PH domain and the membrane embedded region of RhoA support previously unexplained experimental findings regarding the allosteric effects of the PH domain from numerous activity studies of Dbl homology GEF proteins. This work establishes new hypotheses for structural interactivity and allosteric signal modulation in Dbl homology RhoGEFs.

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.

Computational Approaches to Investigate and Design Lipid-binding Domains for Membrane Biosensing

Association of proteins with cellular membranes is critical for signaling and membrane trafficking processes. Many peripheral lipid-binding domains have been identified in the last few decades and have been investigated for their specific lipid-sensing properties using traditional in vivo and in vitro studies. However, several knowledge gaps remain owing to intrinsic limitations of these methodologies. Thus, novel approaches are necessary to further our understanding in lipid-protein biology. This review briefly discusses lipid-binding domains that act as specific lipid biosensors and provides a broad perspective on the computational approaches such as molecular dynamics (MD) simulations and machine learning (ML)-based techniques that can be used to study protein-membrane interactions. We also highlight the need for de novo design of proteins that elicit specific lipid-binding properties.

Simulations of Kindlin-2 PIP binding domains reveal protonation-dependent membrane binding modes

Kindlin-2, a member of the Kindlin family of peripheral membrane proteins, is important for integrin activation and stabilization of epidermal growth factor receptor. It associates with the cytoplasmic face of the plasma membrane via dedicated phosphatidylinositol phosphate binding domains located in the N-terminal F0 and Pleckstrin Homology domains. These domains have binding affinity for phosphatidylinositol 4,5-bisphosphate and, to a greater degree, phosphatidylinositol 3,4,5-trisphosphate. The biological significance of the differential binding of these phosphatidylinositol phosphates to Kindlin-2 and the mechanism by which they activate Kindlin-2 are not well understood. Recently, ssNMR identified the predominant protonation states of phosphatidylinositol 4,5-bisphosphate and phosphatidylinositol 3,4,5-trisphosphate near physiological pH in the presence of anionic lipids. Here, we perform atomistic simulation of the bound state of the Pleckstrin Homology and F0 domains of Kindlin-2 at membranes containing phosphatidylinositol 4,5-bisphosphate/phosphatidylinositol 3,4,5-trisphosphate with differing protonation states. This computational approach demonstrates that these two phosphatidylinositol phosphates differently modulate Kindlin-2 subdomain binding in a protonation-state-dependent manner. We speculate these variations in binding mode provide a mechanism for intracellular pH and Ca²⁺ influx to control the membrane binding behavior and activity of Kindlin-2.

Membrane-binding mechanism of the EEA1 FYVE domain revealed by multi-scale molecular dynamics simulations

Early Endosomal Antigen 1 (EEA1) is a key protein in endosomal trafficking and is implicated in both autoimmune and neurological diseases. The C-terminal FYVE domain of EEA1 binds endosomal membranes, which contain phosphatidylinositol-3-phosphate (PI(3)P). Although it is known that FYVE binds PI(3)P specifically, it has not previously been described of how FYVE attaches and binds to endosomal membranes. In this study, we employed both coarse-grained (CG) and atomistic (AT) molecular dynamics (MD) simulations to determine how FYVE binds to PI(3)P-containing membranes. CG-MD showed that the dominant membrane binding mode resembles the crystal structure of EEA1 FYVE domain in complex with inositol-1,3-diphospate (PDB ID 1JOC). FYVE, which is a homodimer, binds the membrane via a hinge mechanism, where the C-terminus of one monomer first attaches to the membrane, followed by the C-terminus of the other monomer. The estimated total binding energy is ~70 kJ/mol, of which 50-60 kJ/mol stems from specific PI(3)P-interactions. By AT-MD, we could partition the binding mode into two types: (i) adhesion by electrostatic FYVE-PI(3)P interaction, and (ii) insertion of amphipathic loops. The AT simulations also demonstrated flexibility within the FYVE homodimer between the C-terminal heads and coiled-coil stem. This leads to a dynamic model whereby the 200 nm long coiled coil attached to the FYVE domain dimer can amplify local hinge-bending motions such that the Rab5-binding domain at the other end of the coiled coil can explore an area of 0.1 μm2 in the search for a second endosome with which to interact.

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

Flexible pivoting of dynamin pleckstrin homology domain catalyzes fission: insights into molecular degrees of freedom

The neuronal dynamin1 functions in the release of synaptic vesicles by orchestrating the process of GTPase-dependent membrane fission. Dynamin1 associates with the plasma membrane–localized phosphatidylinositol-4,5-bisphosphate (PIP₂) through the centrally located pleckstrin homology domain (PHD). The PHD is dispensable as fission (in model membranes) can be managed, even when the PHD-PIP₂ interaction is replaced by a generic polyhistidine- or polylysine-lipid interaction. However, the absence of the PHD renders a dramatic dampening of the rate of fission. These observations suggest that the PHD-PIP₂–containing membrane interaction could have evolved to expedite fission to fulfill the requirement of rapid kinetics of synaptic vesicle recycling. Here, we use a suite of multiscale modeling approaches to explore PHD–membrane interactions. Our results reveal that 1) the binding of PHD to PIP₂-containing membranes modulates the lipids toward fission-favoring conformations and softens the membrane, and 2) PHD associates with membrane in multiple orientations using variable loops as pivots. We identify a new loop (VL4), which acts as an auxiliary pivot and modulates the orientation flexibility of PHD on the membrane—a mechanism that we believe may be important for high-fidelity dynamin collar assembly. Together, these insights provide a molecular-level understanding of the catalytic role of PHD in dynamin-mediated membrane fission.

Attracted to membranes: lipid-binding domains in plants

Membranes are essential for cells and organelles to function. As membranes are impermeable to most polar and charged molecules, they provide electrochemical energy to transport molecules across and create compartmentalized microenvironments for specific enzymatic and cellular processes. Membranes are also responsible for guided transport of cargoes between organelles and during endo- and exocytosis. In addition, membranes play key roles in cell signaling by hosting receptors and signal transducers and as substrates and products of lipid second messengers. Anionic lipids and their specific interaction with target proteins play an essential role in these processes, which are facilitated by specific lipid-binding domains. Protein crystallography, lipid-binding studies, subcellular localization analyses, and computer modeling have greatly advanced our knowledge over the years of how these domains achieve precision binding and what their function is in signaling and membrane trafficking, as well as in plant development and stress acclimation.

Resistance to Turnip Mosaic Virus in the Family Brassicaceae

Resistance to diseases caused by turnip mosaic virus (TuMV) in crop species of the family Brassicaceae has been studied extensively, especially in members of the genus Brassica. The variation in response observed on resistant and susceptible plants inoculated with different isolates of TuMV is due to a combination of the variation in the plant resistome and the variation in the virus genome. Here, we review the breadth of this variation, both at the level of variation in TuMV sequences, with one eye towards the phylogeny and evolution of the virus, and another eye towards the nature of the various responses observed in susceptible vs. different types of resistance responses. The analyses of the viral genomes allowed comparisons of pathotyped viruses on particular indicator hosts to produce clusters of host types, while the inclusion of phylogeny data and geographic location allowed the formation of the host/geographic cluster groups, the derivation of both of which are presented here. Various studies on resistance determination in particular brassica crops sometimes led to further genetic studies, in many cases to include the mapping of genes, and in some cases to the actual identification of the genes. In addition to summarizing the results from such studies done in brassica crops, as well as in radish and Arabidopsis (the latter as a potential source of candidate genes for brassica and radish), we also summarize work done using nonconventional approaches to obtaining resistance to TuMV.

Enfuvirtide, an HIV-1 fusion inhibitor peptide, can act as a potent SARS-CoV-2 fusion inhibitor: an in silico drug repurposing study

Regarding the urgency of therapeutic measures for coronavirus disease 2019 (COVID-19) pandemic, the use of available drugs with FDA approval is preferred because of the less time and cost required for their development. In silico drug repurposing is an accurate way to speed up the screening of the existing FDA-approved drugs to find a therapeutic option for COVID-19. The similarity in SARS-CoV-2 and HIV-1 fusion mechanism to host cells can be a key point for Inhibit SARS-CoV-2 entry into host cells by HIV fusion inhibitors. Accordingly, in this study, an HIV-1 fusion inhibitor called Enfuvirtide (Enf) was selected. The affinity and essential residues involving in the Enf binding to the S2 protein of SARS-CoV-2, HIV-1 gp41 protein and angiotensin-converting enzyme 2 (ACE-2) as a negative control, was evaluated using molecular docking. Eventually, Enf-S2 and Enf-gp41 protein complexes were simulated by molecular dynamics (MD) in terms of binding affinity and stability. Based on the most important criteria such as docking score, cluster size, energy and dissociation constant, the strongest interaction was observed between Enf with the S2 protein. In addition, MD results confirmed that Enf-S2 protein interaction was remarkably stable and caused the S2 protein residues to undergo the fewest fluctuations. In conclusion, it can be stated that Enf can act as a strong SARS-CoV-2 fusion inhibitor and demonstrates the potential to enter the clinical trial phase of COVID-19.

Communicated by Ramaswamy H. Sarma

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.

Evaluating inositol phospholipid interactions with inward rectifier potassium channels and characterising their role in disease

Membrane proteins are frequently modulated by specific protein-lipid interactions. The activation of human inward rectifying potassium (hKir) channels by phosphoinositides (PI) has been well characterised. Here, we apply a coarse-grained molecular dynamics free-energy perturbation (CG-FEP) protocol to capture the energetics of binding of PI lipids to hKir channels. By using either a single- or multi-step approach, we establish a consistent value for the binding of PIP2 to hKir channels, relative to the binding of the bulk phosphatidylcholine phospholipid. Furthermore, by perturbing amino acid side chains on hKir6.2, we show that the neonatal diabetes mutation E179K increases PIP2 affinity, while the congenital hyperinsulinism mutation K67N results in a reduced affinity. We show good agreement with electrophysiological data where E179K exhibits a reduction in neomycin sensitivity, implying that PIP2 binds more tightly E179K channels. This illustrates the application of CG-FEP to compare affinities between lipid species, and for annotating amino acid residues.

Regulation of V-ATPase Activity and Organelle pH by Phosphatidylinositol Phosphate Lipids

Luminal pH and the distinctive distribution of phosphatidylinositol phosphate (PIP) lipids are central identifying features of organelles in all eukaryotic cells that are also critical for organelle function. V-ATPases are conserved proton pumps that populate and acidify multiple organelles of the secretory and the endocytic pathway. Complete loss of V-ATPase activity causes embryonic lethality in higher animals and conditional lethality in yeast, while partial loss of V-ATPase function is associated with multiple disease states. On the other hand, many cancer cells increase their virulence by upregulating V-ATPase expression and activity. The pH of individual organelles is tightly controlled and essential for function, but the mechanisms for compartment-specific pH regulation are not completely understood. There is substantial evidence indicating that the PIP content of membranes influences organelle pH. We present recent evidence that PIPs interact directly with subunit isoforms of the V-ATPase to dictate localization of V-ATPase subpopulations and participate in their regulation. In yeast cells, which have only one set of organelle-specific V-ATPase subunit isoforms, the Golgi-enriched lipid PI(4)P binds to the cytosolic domain of the Golgi-enriched a-subunit isoform Stv1, and loss of PI(4)P binding results in mislocalization of Stv1-containing V-ATPases from the Golgi to the vacuole/lysosome. In contrast, levels of the vacuole/lysosome-enriched signaling lipid PI(3,5)P₂ affect assembly and activity of V-ATPases containing the Vph1 a-subunit isoform. Mutations in the Vph1 isoform that disrupt the lipid interaction increase sensitivity to stress. These studies have decoded “zip codes” for PIP lipids in the cytosolic N-terminal domain of the a-subunit isoforms of the yeast V-ATPase, and similar interactions between PIP lipids and the V-ATPase subunit isoforms are emerging in higher eukaryotes. In addition to direct effects on the V-ATPase, PIP lipids are also likely to affect organelle pH indirectly, through interactions with other membrane transporters. We discuss direct and indirect effects of PIP lipids on organelle pH, and the functional consequences of the interplay between PIP lipid content and organelle pH.

Capturing Choline-Aromatics Cation-π Interactions in the MARTINI Force Field

Cation−π interactions play an important role in biomolecular recognition, including interactions between membrane phosphatidylcholine lipids and aromatic amino acids of peripheral proteins. While molecular mechanics coarse grain (CG) force fields are particularly well suited to simulate membrane proteins in general, they are not parameterized to explicitly reproduce cation−π interactions. We here propose a modification of the polarizable MARTINI coarse grain (CG) model enabling it to model membrane binding events of peripheral proteins whose aromatic amino acid interactions with choline headgroups are crucial for their membrane binding. For this purpose, we first collected and curated a dataset of eight peripheral proteins from different families. We find that the MARTINI CG model expectedly underestimates aromatics–choline interactions and is unable to reproduce membrane binding of the peripheral proteins in our dataset. Adjustments of the relevant interactions in the polarizable MARTINI force field yield significant improvements in the observed binding events. The orientation of each membrane-bound protein is comparable to reference data from all-atom simulations and experimental binding data. We also use negative controls to ensure that choline–aromatics interactions are not overestimated. We finally check that membrane properties, transmembrane proteins, and membrane translocation potential of mean force (PMF) of aromatic amino acid side-chain analogues are not affected by the new parameter set. This new version “MARTINI 2.3P” is a significant improvement over its predecessors and is suitable for modeling membrane proteins including peripheral membrane binding of peptides and proteins.

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.

Multiple lipid binding sites determine the affinity of PH domains for phosphoinositide-containing membranes

Association of peripheral proteins with lipid bilayers regulates membrane signaling and dynamics. Pleckstrin homology (PH) domains bind to phosphatidylinositol phosphate (PIP) molecules in membranes. The effects of local PIP enrichment on the interaction of PH domains with membranes is unclear. Molecular dynamics simulations allow estimation of the binding energy of GRP1 PH domain to PIP3-containing membranes. The free energy of interaction of the PH domain with more than two PIP3 molecules is comparable to experimental values, suggesting that PH domain binding involves local clustering of PIP molecules within membranes. We describe a mechanism of PH binding proceeding via an encounter state to two bound states which differ in the orientation of the protein relative to the membrane, these orientations depending on the local PIP concentration. These results suggest that nanoscale clustering of PIP molecules can control the strength and orientation of PH domain interaction in a concentration-dependent manner.

In silico assessment of human Calprotectin subunits (S100A8/A9) in presence of sodium and calcium ions using Molecular Dynamics simulation approach

Calprotectin is a heterodimeric protein complex which consists of two subunits including S100A8 and S100A9. This protein has a major role in different inflammatory disease and various types of cancers. In current study we aimed to evaluate the structural and thermodynamic changes of the subunits and the complex in presence of sodium and calcium ions using molecular dynamics (MD) simulation. Therefore, the residue interaction network (RIN) was visualized in Cytoscape program. In next step, to measure the binding free energy, the potential of mean force (PMF) method was performed. Finally, the molecular mechanics Poisson-Boltzmann surface area (MMPBSA) method was applied as an effective tool to calculate the molecular model affinities. The MD simulation results of the subunits represented their structural changes in presence of Ca2+. Moreover, the RIN and Hydrogen bond analysis demonstrated that cluster interactions between Calprotectin subunits in presence of Ca2+ were greater in comparison with Na+. Our findings indicated that the binding free energy of the subunits in presence of Ca2+ was significantly greater than Na+. The results revealed that Ca2+ has the ability to induce structural changes in subunits in comparison with Na+ which lead to create stronger interactions between. Hence, studying the physical characteristics of the human proteins could be considered as a powerful tool in theranostics and drug design purposes.

Membrane Recognition and Binding by the Phosphatidylinositol Phosphate Kinase PIP5K1A: A Multiscale Simulation Study

Phosphatidylinositol phosphates (PIPs) are lipid signaling molecules that play key roles in many cellular processes. PIP5K1A kinase catalyzes phosphorylation of PI4P to form PIP₂, which in turn interacts with membrane and membrane-associated proteins. We explore the mechanism of membrane binding by the PIP5K1A kinase using a multiscale molecular dynamics approach. Coarse-grained simulations show binding of monomeric PIP5K1A to a model cell membrane containing PI4P. PIP5K1A did not bind to zwitterionic or anionic membranes lacking PIP molecules. Initial encounter of kinase and bilayer was followed by reorientation to enable productive binding to the PI4P-containing membrane. The simulations suggest that unstructured regions may be important for the preferred orientation for membrane binding. Atomistic simulations indicated that the dimeric kinase could not bind to the membrane via both active sites at the same time, suggesting a conformational change in the protein and/or bilayer distortion may be needed for dual-site binding to occur.

•

PIP5K1A kinase interacts with PIP-containing membranes via its activation loop

•

PIP5K1A does not bind to zwitterionic or anionic membranes lacking PIP molecules

•

Initial encounter of protein and bilayer is followed by reorientation and binding

•

Dimeric PIP5K1A binds with membrane contacts via only one catalytic site at a time

Membrane Curvature Sensing by Amphipathic Helices: Insights from Implicit Membrane Modeling

Sensing and generation of lipid membrane curvature, mediated by the binding of specific proteins onto the membrane surface, play crucial roles in cell biology. A number of mechanisms have been proposed, but the molecular understanding of these processes is incomplete. All-atom molecular dynamics simulations have offered valuable insights but are extremely demanding computationally. Implicit membrane simulations could provide a viable alternative, but current models apply only to planar membranes. In this work, the implicit membrane model 1 is extended to spherical and tubular membranes. The geometric change from planar to curved shapes is straightforward but insufficient for capturing the full curvature effect, which includes changes in lipid packing. Here, these packing effects are taken into account via the lateral pressure profile. The extended implicit membrane model 1 is tested on the wild-types and mutants of the antimicrobial peptide magainin, the ALPS motif of arfgap1, α-synuclein, and an ENTH domain. In these systems, the model is in qualitative agreement with experiments. We confirm that favorable electrostatic interactions tend to weaken curvature sensitivity in the presence of strong hydrophobic interactions but may actually have a positive effect when those are weak. We also find that binding to vesicles is more favorable than binding to tubes of the same diameter and that the long helix of α-synuclein tends to orient along the axis of tubes, whereas shorter helices tend to orient perpendicular to it. Adoption of a specific orientation could provide a mechanism for coupling protein oligomerization to tubule formation.

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.

Genome-wide association study reveals new loci involved in Arabidopsis thaliana and Turnip mosaic virus (TuMV) interactions in the field

The genetic architecture of plant response to viruses has often been studied in model nonnatural pathosystems under controlled conditions. There is an urgent need to elucidate the genetic architecture of the response to viruses in a natural setting. A field experiment was performed in each of two years. In total, 317 Arabidopsis thaliana accessions were inoculated with its natural Turnip mosaic virus (TuMV). The accessions were phenotyped for viral accumulation, frequency of infected plants, stem length and symptoms. Genome-wide association mapping was performed. Arabidopsis thaliana exhibits extensive natural variation in its response to TuMV in the field. The underlying genetic architecture reveals a more quantitative picture than in controlled conditions. Ten genomic regions were consistently identified across the two years. RTM3 (Restricted TEV Movement 3) is a major candidate for the response to TuMV in the field. New candidate genes include Dead box helicase 1, a Tim Barrel domain protein and the eukaryotic translation initiation factor eIF3b. To our knowledge, this study is the first to report the genetic architecture of quantitative response of A. thaliana to a naturally occurring virus in a field environment, thereby highlighting relevant candidate genes involved in plant virus interactions in nature.

Flexible lipid nanomaterials studied by NMR spectroscopy

Advances in solid-state nuclear magnetic resonance spectroscopy inform the emergence of material properties from atomistic-level interactions in membrane lipid nanostructures.

Regulation of actin assembly by PI(4,5)P2 and other inositol phospholipids: An update on possible mechanisms

Actin cytoskeleton dynamics depend on a tight regulation of actin filament formation from an intracellular pool of monomers, followed by their linkage to each other or to cell membranes, followed by their depolymerization into a fresh pool of actin monomers. The ubiquitous requirement for continuous actin remodeling that is necessary for many cellular functions is orchestrated in large part by actin binding proteins whose affinity for actin is altered by inositol phospholipids, most prominently PI(4,5)P2 (phosphatidylinositol 4,5-bisphosphate). The kinetics of PI(4,5)P2 synthesis and hydrolysis, its lateral distribution within the lipid bilayer, and coincident detection of PI(4,5)P2 and another signal, all play a role in determining when and where a particular PI(4,5)P2-regulated protein is inactivated or activated to exert its effect on the actin cytoskeleton. This review summarizes a range of models that have been developed to explain how PI(4,5)P2 might function in the complex chemical and structural environment of the cell based on a combination of experiment and computational studies.

Interaction of the human erythrocyte Band 3 anion exchanger 1 (AE1, SLC4A1) with lipids and glycophorin A: Molecular organization of the Wright (Wr) blood group antigen

The Band 3 (AE1, SLC4A1) membrane protein is found in red blood cells and in kidney where it functions as an electro-neutral chloride/bicarbonate exchanger. In this study, we have used molecular dynamics simulations to provide the first realistic model of the dimeric membrane domain of human Band 3 in an asymmetric lipid bilayer containing a full complement of phospholipids, including phosphatidylinositol 4,5–bisphosphate (PIP₂) and cholesterol, and its partner membrane protein Glycophorin A (GPA). The simulations show that the annular layer in the inner leaflet surrounding Band 3 was enriched in phosphatidylserine and PIP₂ molecules. Cholesterol was also enriched around Band 3 but also at the dimer interface. The interaction of these lipids with specific sites on Band 3 may play a role in the folding and function of this anion transport membrane protein. GPA associates with Band 3 to form the Wright (Wr) blood group antigen, an interaction that involves an ionic bond between Glu658 in Band 3 and Arg61 in GPA. We were able to recreate this complex by performing simulations to allow the dimeric transmembrane portion of GPA to interact with Band 3 in a model membrane. Large-scale simulations showed that the GPA dimer can bridge Band 3 dimers resulting in the dynamic formation of long strands of alternating Band 3 and GPA dimers.

Human Band 3 (AE1, SLC4A1), an abundant 911 amino acid glycoprotein, catalyzes the exchange of bicarbonate and chloride across the red blood cell membrane, a process necessary for efficient respiration. Malfunction of Band 3 leads to inherited diseases such as Southeast Asian Ovalocytosis, hereditary spherocytosis and distal renal tubular acidosis. Despite much available structural and functional data about Band 3, key questions about the conformational changes associated with transport and the molecular details of its interaction with lipids and other proteins remain unanswered. In this study, we have used computer simulations to investigate the dynamics of Band 3 in lipid bilayers that resemble the red blood cell plasma membrane. Our results suggest that negatively charged phospholipids and cholesterol interact strongly with Band 3 forming an annulus around the protein. Glycophorin A (GPA) interacts with Band 3 to form the Wright (Wr) blood group antigen. We were able to recreate this complex and show that GPA promotes the clustering of Band 3 in red blood cell membranes. Understanding the molecular details of the interaction of Band 3 with GPA has provided new insights into the nature of the Wright blood group antigen.

Interactions of the EphA2 Kinase Domain with PIPs in Membranes: Implications for Receptor Function

EphA2 is a member of the receptor tyrosine kinase family. Interactions of the cytoplasmic region of EphA2 with the cell membrane are functionally important and yet remain incompletely characterized. Molecular dynamics simulations combined with biochemical studies reveal the interactions of the transmembrane, juxtamembrane (JM), and kinase domains with the membrane. We describe how the kinase domain is oriented relative to the membrane and how the JM region can modulate this interaction. We highlight the role of phosphatidylinositol phosphates (PIPs) in mediating the interaction of the kinase domain with the membrane and, conversely, how positively charged patches at the kinase surface and in the JM region induce the formation of nanoclusters of PIP molecules in the membrane. Integration of these results with those from previous studies enable computational reconstitution of a near complete EphA2 receptor within a membrane, suggesting a role for receptor-lipid interactions in modulation of EphA2.

Sensitivity of peripheral membrane proteins to the membrane context: A case study of phosphatidylserine and the TIM proteins

There is a diverse class of peripheral membrane-binding proteins that specifically bind phosphatidylserine (PS), a lipid that signals apoptosis or cell fusion depending on the membrane context of its presentation. PS-receptors are specialized for particular PS-presenting pathways, indicating that they might be sensitive to the membrane context. In this review, we describe a combination of thermodynamic, structural, and computational techniques that can be used to investigate the mechanisms underlying this sensitivity. As an example, we focus on three PS-receptors of the T-cell Immunoglobulin and Mucin containing (TIM) protein family, which we have previously shown to differ in their sensitivity to PS surface density.