A model for hydrophobic protrusions on peripheral membrane proteins

Lock and key: Quest to find the most compatible membrane mimetic for studying membrane proteins in native environment

Membrane proteins play crucial roles in cellular signal transduction, molecule transport, host-pathogen interactions, and metabolic processes. However, mutations, changes in membrane properties, and environmental factors can lead to loss of protein function. This results in impaired ligand binding and misfolded structures that prevent proteins from adopting their native conformation. Many membrane proteins are also therapeutic targets in various diseases, where drugs can either restore or inhibit their specific functions. Understanding membrane protein structure and function is vital for advancing cell biology and physiology. Experimental studies often involve extracting proteins from their native environments and reconstituting them in membrane mimetics like detergents, bicelles, amphipols, nanodiscs, and liposomes. These mimetics replicate aspects of native membranes, aiding in the study of protein behavior outside living cells. Scientists continuously explore new, more native-like membrane mimetics to improve experimental accuracy. This dynamic field involves evaluating the advantages and disadvantages of different mimetics and optimizing the reconstitution process to better mimic natural conditions.

Analysis of protein-protein and protein-membrane interactions by isotope-edited infrared spectroscopy

The objective of this work is to highlight the power of isotope-edited Fourier transform infrared (FTIR) spectroscopy in resolving important problems encountered in biochemistry, biophysics, and biomedical research, focusing on protein-protein and protein membrane interactions that play key roles in practically all life processes. An overview of the effects of isotope substitutions in (bio)molecules on spectral frequencies and intensities is given. Data are presented demonstrating how isotope-labeled proteins and/or lipids can be used to elucidate enzymatic mechanisms, the mode of membrane binding of peripheral proteins, regulation of membrane protein function, protein aggregation, and local and global structural changes in proteins during functional transitions. The use of polarized attenuated total reflection FTIR spectroscopy to identify the spatial orientation and the secondary structure of a membrane-bound interfacial enzyme and the mode of lipid hydrolysis is described. Methods of production of site-directed, segmental, and domain-specific labeling of proteins by the synthetic, semisynthetic, and recombinant strategies, including advanced protein engineering technologies such as nonsense suppression and frameshift quadruplet codons are overviewed.

Assembly of Protein Complexes in and on the Membrane with Predicted Spatial Arrangement Constraints

Membrane proteins play crucial roles in various cellular processes, and their interactions with other proteins in and on the membrane are essential for their proper functioning. While an increasing number of structures of more membrane proteins are being determined, the available structure data is still sparse. To gain insights into the mechanisms of membrane protein complexes, computational docking methods are necessary due to the challenge of experimental determination. Here, we introduce Mem-LZerD, a rigid-body membrane docking algorithm designed to take advantage of modern membrane modeling and protein docking techniques to facilitate the docking of membrane protein complexes. Mem-LZerD is based on the LZerD protein docking algorithm, which has been constantly among the top servers in many rounds of CAPRI protein docking assessment. By employing a combination of geometric hashing, newly constrained by the predicted membrane height and tilt angle, and model scoring accounting for the energy of membrane insertion, we demonstrate the capability of Mem-LZerD to model diverse membrane protein-protein complexes. Mem-LZerD successfully performed unbound docking on 13 of 21 (61.9%) transmembrane complexes in an established benchmark, more than shown by previous approaches. It was additionally tested on new datasets of 44 transmembrane complexes and 92 peripheral membrane protein complexes, of which it successfully modeled 35 (79.5%) and 15 (16.3%) complexes respectively. When non-blind orientations of peripheral targets were included, the number of successes increased to 54 (58.7%). We further demonstrate that Mem-LZerD produces complex models which are suitable for molecular dynamics simulation. Mem-LZerD is made available at https://lzerd.kiharalab.org.

A Perspective on the Prospective Use of AI in Protein Structure Prediction

AlphaFold2 (AF2) and RoseTTaFold (RF) have revolutionized structural biology, serving as highly reliable and effective methods for predicting protein structures. This article explores their impact and limitations, focusing on their integration into experimental pipelines and their application in diverse protein classes, including membrane proteins, intrinsically disordered proteins (IDPs), and oligomers. In experimental pipelines, AF2 models help X-ray crystallography in resolving the phase problem, while complementarity with mass spectrometry and NMR data enhances structure determination and protein flexibility prediction. Predicting the structure of membrane proteins remains challenging for both AF2 and RF due to difficulties in capturing conformational ensembles and interactions with the membrane. Improvements in incorporating membrane-specific features and predicting the structural effect of mutations are crucial. For intrinsically disordered proteins, AF2's confidence score (pLDDT) serves as a competitive disorder predictor, but integrative approaches including molecular dynamics (MD) simulations or hydrophobic cluster analyses are advocated for accurate dynamics representation. AF2 and RF show promising results for oligomeric models, outperforming traditional docking methods, with AlphaFold-Multimer showing improved performance. However, some caveats remain in particular for membrane proteins. Real-life examples demonstrate AF2's predictive capabilities in unknown protein structures, but models should be evaluated for their agreement with experimental data. Furthermore, AF2 models can be used complementarily with MD simulations. In this Perspective, we propose a "wish list" for improving deep-learning-based protein folding prediction models, including using experimental data as constraints and modifying models with binding partners or post-translational modifications. Additionally, a meta-tool for ranking and suggesting composite models is suggested, driving future advancements in this rapidly evolving field.

Martini 3 Force Field Parameters for Protein Lipidation Post-Translational Modifications

Protein lipidations are vital co/post-translational modifications that tether lipid tails to specific protein amino acids, allowing them to anchor to biological membranes, switch their subcellular localization, and modulate association with other proteins. Such lipidations are thus crucial for multiple biological processes including signal transduction, protein trafficking, and membrane localization and are implicated in various diseases as well. Examples of lipid-anchored proteins include the Ras family of proteins that undergo farnesylation; actin and gelsolin that are myristoylated; phospholipase D that is palmitoylated; glycosylphosphatidylinositol-anchored proteins; and others. Here, we develop parameters for cysteine-targeting farnesylation, geranylgeranylation, and palmitoylation, as well as glycine-targeting myristoylation for the latest version of the Martini 3 coarse-grained force field. The parameters are developed using the CHARMM36m all-atom force field parameters as reference. The behavior of the coarse-grained models is consistent with that of the all-atom force field for all lipidations and reproduces key dynamical and structural features of lipid-anchored peptides, such as the solvent-accessible surface area, bilayer penetration depth, and representative conformations of the anchors. The parameters are also validated in simulations of the lipid-anchored peripheral membrane proteins Rheb and Arf1, after comparison with independent all-atom simulations. The parameters, along with mapping schemes for the popular martinize2 tool, are available for download at 10.5281/zenodo.7849262 and also as supporting information.

Effects of biophysical membrane properties on recognition of phosphatidylserine, or phosphatidylinositol 4-phosphate by lipid biosensors LactC2, or P4M

Lipid biosensors are molecular tools used both in vivo and in vitro applications, capable of selectively detecting specific types of lipids in biological membranes. However, despite their extensive use, there is a lack of systematic characterization of their binding properties in various membrane conditions. The purpose of this study was to investigate the impact of membrane properties, such as fluidity and membrane charge, on the sensitivity of two lipid biosensors, LactC2 and P4M, to their target lipids, phosphatidylserine (PS) or phosphatidylinositol 4-phosphate (PI4P), respectively. Dual-color fluorescence cross-correlation spectroscopy, employed in this study, provided a useful technique to investigate interactions of these recombinant fluorescent biosensors with liposomes of varying compositions. The results of the study demonstrate that the binding of the LactC2 biosensor to low levels of PS in the membrane is highly supported by the presence of anionic lipids or membrane fluidity. However, at high PS levels, the presence of anionic lipids does not further enhance binding of LactC2. In contrast, neither membrane charge, nor membrane fluidity significantly affect the binding affinity of P4M to PI4P. These findings provide valuable insights into the role of membrane properties on the binding properties of lipid biosensors.

Assembly of Protein Complexes In and On the Membrane with Predicted Spatial Arrangement Constraints

Membrane proteins play crucial roles in various cellular processes, and their interactions with other proteins in and on the membrane are essential for their proper functioning. While an increasing number of structures of more membrane proteins are being determined, the available structure data is still sparse. To gain insights into the mechanisms of membrane protein complexes, computational docking methods are necessary due to the challenge of experimental determination. Here, we introduce Mem-LZerD, a rigid-body membrane docking algorithm designed to take advantage of modern membrane modeling and protein docking techniques to facilitate the docking of membrane protein complexes. Mem-LZerD is based on the LZerD protein docking algorithm, which has been constantly among the top servers in many rounds of CAPRI protein docking assessment. By employing a combination of geometric hashing, newly constrained by the predicted membrane height and tilt angle, and model scoring accounting for the energy of membrane insertion, we demonstrate the capability of Mem-LZerD to model diverse membrane protein-protein complexes. Mem-LZerD successfully performed unbound docking on 13 of 21 (61.9%) transmembrane complexes in an established benchmark, more than shown by previous approaches. It was additionally tested on new datasets of 44 transmembrane complexes and 92 peripheral membrane protein complexes, of which it successfully modeled 35 (79.5%) and 15 (16.3%) complexes respectively. When non-blind orientations of peripheral targets were included, the number of successes increased to 54 (58.7%). We further demonstrate that Mem-LZerD produces complex models which are suitable for molecular dynamics simulation. Mem-LZerD is made available at https://lzerd.kiharalab.org.

A coarse-grained molecular dynamics investigation on spontaneous binding of Aβ1-40 fibrils with cholesterol-mixed DPPC bilayers

Alzheimer's disease is the most common form of dementia. Its aetiology is characterized by the misfolding and aggregation of amyloid-β (Aβ) peptides into β-sheet-rich Aβ oligomers/fibrils. Although multiple experimental studies have suggested that Aβ oligomers/fibrils interact with the cell membranes and perturb their structures and dynamics, the molecular mechanism of this interaction is still not fully understood. In the present work, we have performed a total of 120 μs-long simulations to investigate the interaction between trimeric or hexameric Aβ1-40 fibrils with either a 100% DPPC bilayer, a 70% DPPC-30% cholesterol bilayer or a 50% DPPC-50% cholesterol bilayer. Our simulation data capture the spontaneous binding of the aqueous Aβ1-40 fibrils with the membranes and show that the central hydrophobic amino acid cluster, the lysine residue adjacent to it and the C-terminal hydrophobic residues are all involved in the process. Moreover, our data show that while the Aβ1-40 fibril does not bind to the 100% DPPC bilayer, its binding affinity for the membrane increases with the amount of cholesterol. Overall, our data suggest that two clusters of hydrophobic residues and one lysine help Aβ1-40 fibrils establish stable interactions with a cholesterol-rich DPPC bilayer. These residues are likely to represent potential target regions for the design of inhibitors, thus opening new avenues in structure-based drug design against Aβ oligomer/fibril-membrane interaction.

A novel calcium channel Cavβ2 splice variant with unique properties predominates in the retina

Cavβ subunits are essential for surface expression of voltage-gated calcium channel complexes and crucially modulate biophysical properties like voltage-dependent inactivation. Here, we describe the discovery and characterization of a novel Cavβ2 variant with distinct features that predominates in the retina. We determined spliced exons in retinal transcripts of the Cacnb2 gene, coding for Cavβ2, by RNA-Seq data analysis and quantitative PCR. We cloned a novel Cavβ2 splice variant from mouse retina, which we are calling β2i, and investigated biophysical properties of calcium currents with this variant in a heterologous expression system as well as its intrinsic membrane interaction when expressed alone. Our data showed that β2i predominated in the retina with expression in photoreceptors and bipolar cells. Furthermore, we observed that the β2i N-terminus exhibited an extraordinary concentration of hydrophobic residues, a distinct feature not seen in canonical variants. The biophysical properties resembled known membrane-associated variants, and β2i exhibited both a strong membrane association and a propensity for clustering, which depended on hydrophobic residues in its N-terminus. We considered available Cavβ structure data to elucidate potential mechanisms underlying the observed characteristics but resolved N-terminus structures were lacking and thus, precluded clear conclusions. With this description of a novel N-terminus variant of Cavβ2, we expand the scope of functional variation through N-terminal splicing with a distinct form of membrane attachment. Further investigation of the molecular mechanisms underlying the features of β2i could provide new angles on the way Cavβ subunits modulate Ca2+ channels at the plasma membrane.

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.

Intrinsically Disordered Proteins: An Overview

Many proteins and protein segments cannot attain a single stable three-dimensional structure under physiological conditions; instead, they adopt multiple interconverting conformational states. Such intrinsically disordered proteins or protein segments are highly abundant across proteomes, and are involved in various effector functions. This review focuses on different aspects of disordered proteins and disordered protein regions, which form the basis of the so-called "Disorder-function paradigm" of proteins. Additionally, various experimental approaches and computational tools used for characterizing disordered regions in proteins are discussed. Finally, the role of disordered proteins in diseases and their utility as potential drug targets are explored.

Structure of AQEE-30 of VGF Neuropeptide in Membrane-Mimicking Environments

AQEE-30 is one of the VGF peptides, which are derived from the VGF polypeptide precursor, and related to various physiological phenomena including neuroprotective effects in Huntington's disease and amyotrophic lateral sclerosis (ALS). Although various functions of AQEE-30 have been reported so far, the structure of this peptide has not been reported yet. In this study, the structure of human AQEE-30 was investigated in hexafluoroisopropanol (HFIP) and dodecyl phosphocholine (DPC) micelle solutions, using circular dichroism (CD) and nuclear magnetic resonance (NMR) spectroscopy. CD results showed that AQEE-30 had a partial helical structure in aqueous buffer, and the helical structure was stabilized in the HFIP and DPC micelle solutions. The 3D structures determined by NMR spectroscopy showed that AQEE-30 adopted mainly α-helical structure in both the HFIP and DPC micelle solutions. The surface of AQEE-30 showed that it was predominantly negatively charged. The residues from 601 to 611 in both the HFIP and DPC micelle solutions showed amphiphilicity with four negatively charged residues, glutamate. The C-terminal consecutive arginine residues formed a partial positively charged surface. These results suggest an α-helical active structure of AQEE-30 in the cell-membrane environment.

DREAMM: a web-based server for drugging protein-membrane interfaces as a novel workflow for targeted drug design

SUMMARY

The allosteric modulation of peripheral membrane proteins (PMPs) by targeting protein-membrane interactions with drug-like molecules represents a new promising therapeutic strategy for proteins currently considered undruggable. However, the accessibility of protein-membrane interfaces by small molecules has been so far unexplored, possibly due to the complexity of the interface, the limited protein-membrane structural information and the lack of computational workflows to study it. Herein, we present a pipeline for drugging protein-membrane interfaces using the DREAMM (Drugging pRotein mEmbrAne Machine learning Method) web server. DREAMM works in the back end with a fast and robust ensemble machine learning algorithm for identifying protein-membrane interfaces of PMPs. Additionally, DREAMM also identifies binding pockets in the vicinity of the predicted membrane-penetrating amino acids in protein conformational ensembles provided by the user or generated within DREAMM.

AVAILABILITY AND IMPLEMENTATION

DREAMM web server is accessible via https://dreamm.ni4os.eu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

PGRS domain structures: Doomed to sail the mycomembrane

The impact of artificial intelligence (AI) in understanding biological processes is potentially immense. Structural elucidation of mycobacterial PE_PGRS is sustenance to unveil the role of these enigmatic proteins. We propose a PGRS “sailing” model as a smart tool to diffuse along the mycomembrane, to expose structural motifs for host interactions, and/or to ship functional protein modules at their C-terminus.

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.

Intrinsically disordered proteins play diverse roles in cell signaling

Signaling pathways allow cells to detect and respond to a wide variety of chemical (e.g. Ca2+ or chemokine proteins) and physical stimuli (e.g., sheer stress, light). Together, these pathways form an extensive communication network that regulates basic cell activities and coordinates the function of multiple cells or tissues. The process of cell signaling imposes many demands on the proteins that comprise these pathways, including the abilities to form active and inactive states, and to engage in multiple protein interactions. Furthermore, successful signaling often requires amplifying the signal, regulating or tuning the response to the signal, combining information sourced from multiple pathways, all while ensuring fidelity of the process. This sensitivity, adaptability, and tunability are possible, in part, due to the inclusion of intrinsically disordered regions in many proteins involved in cell signaling. The goal of this collection is to highlight the many roles of intrinsic disorder in cell signaling. Following an overview of resources that can be used to study intrinsically disordered proteins, this review highlights the critical role of intrinsically disordered proteins for signaling in widely diverse organisms (animals, plants, bacteria, fungi), in every category of cell signaling pathway (autocrine, juxtacrine, intracrine, paracrine, and endocrine) and at each stage (ligand, receptor, transducer, effector, terminator) in the cell signaling process. Thus, a cell signaling pathway cannot be fully described without understanding how intrinsically disordered protein regions contribute to its function. The ubiquitous presence of intrinsic disorder in different stages of diverse cell signaling pathways suggest that more mechanisms by which disorder modulates intra- and inter-cell signals remain to be discovered.

Sequence Analysis and Structural Predictions of Lipid Transfer Bridges in the Repeating Beta Groove (RBG) Superfamily Reveal Past and Present Domain Variations Affecting Form, Function and Interactions of VPS13, ATG2, SHIP164, Hobbit and Tweek

Lipid transfer between organelles requires proteins that shield the hydrophobic portions of lipids as they cross the cytoplasm. In the last decade a new structural form of lipid transfer protein (LTP) has been found: long hydrophobic grooves made of beta-sheet that bridge between organelles at membrane contact sites. Eukaryotes have five families of bridge-like LTPs: VPS13, ATG2, SHIP164, Hobbit and Tweek. These are unified into a single superfamily through their bridges being composed of just one domain, called the repeating beta groove (RBG) domain, which builds into rod shaped multimers with a hydrophobic-lined groove and hydrophilic exterior. Here, sequences and predicted structures of the RBG superfamily were analyzed in depth. Phylogenetics showed that the last eukaryotic common ancestor contained all five RBG proteins, with duplicated VPS13s. The current set of long RBG protein appears to have arisen in even earlier ancestors from shorter forms with 4 RBG domains. The extreme ends of most RBG proteins have amphipathic helices that might be an adaptation for direct or indirect bilayer interaction, although this has yet to be tested. The one exception to this is the C-terminus of SHIP164, which instead has a coiled-coil. Finally, the exterior surfaces of the RBG bridges are shown to have conserved residues along most of their length, indicating sites for partner interactions almost all of which are unknown. These findings can inform future cell biological and biochemical experiments.

Molecular Characterization of an Intrinsically Disordered Chaperone Reveals Net-Charge Regulation in Chaperone Action

Molecular chaperones are diverse biomacromolecules involved in the maintenance of cellular protein homeostasis (proteostasis). Here we demonstrate that in contrast to most chaperones with defined three-dimensional structures, the acid-inducible protein Asr in Escherichia coli is intrinsically disordered and exhibits varied aggregation-preventing or aggregation-promoting activities, acting as a "conditionally active chaperone". Bioinformatics and experimental analyses of Asr showed that it is devoid of hydrophobic patches but rich in positive charges and local polyproline II backbone structures. Asr contributes to the integrity of the bacterial outer membrane under mildly acidic conditions in vivo and possesses chaperone activities toward model clients in vitro. Notably, its chaperone activity is dependent on the net charges of clients: on the one hand, it inhibits the aggregation of clients with similar net charges; on the other hand, it stimulates the aggregation of clients with opposite net charges. Mutational analysis confirmed that positively charged residues in Asr are essential for the varied effects on protein aggregation, suggesting that electrostatic interactions are the major driving forces underlying Asr's proteostasis-related activity. These findings present a unique example of an intrinsically disordered molecular chaperone with distinctive dual functions-as an aggregase or as a chaperone-depending on the net charges of clients.

Intrinsically disordered proteins and membranes: a marriage of convenience for cell signalling?

The structure-function paradigm has guided investigations into the molecules involved in cellular signalling for decades. The peripheries of this paradigm, however, start to unravel when considering the co-operation between proteins and the membrane in signalling processes. Intrinsically disordered regions hold distinct advantages over folded domains in terms of their binding promiscuity, sensitivity to their particular environment and their ease of modulation through post-translational modifications. Low sequence complexity and bias towards charged residues are also favourable for the multivalent electrostatic interactions that occur at the surfaces of lipid bilayers. This review looks at the principles behind the successful marriage between protein disorder and membranes in addition to the role of this partnership in modifying and regulating signalling in cellular processes. The HVR (hypervariable region) of small GTPases is highlighted as a well-studied example of the nuanced role a short intrinsically disordered region can play in the fine-tuning of signalling pathways.

Crystal Structure Reveals the Full Ras-Raf Interface and Advances Mechanistic Understanding of Raf Activation

Ras and Raf-kinase interact through the Ras-binding (RBD) and cysteine-rich domains (CRD) of Raf to signal through the mitogen-activated protein kinase pathway, yet the molecular mechanism leading to Raf activation has remained elusive. We present the 2.8 Å crystal structure of the HRas-CRaf-RBD_CRD complex showing the Ras-Raf interface as a continuous surface on Ras, as seen in the KRas-CRaf-RBD_CRD structure. In molecular dynamics simulations of a Ras dimer model formed through the α4-α5 interface, the CRD is dynamic and located between the two Ras protomers, poised for direct or allosteric modulation of functionally relevant regions of Ras and Raf. We propose a molecular model in which Ras binding is involved in the release of Raf autoinhibition while the Ras-Raf complex dimerizes to promote a platform for signal amplification, with Raf-CRD centrally located to impact regulation and function.

Effect of Methionine Sulfoxide on the Synthesis and Purification of Aggregation-Prone Peptides

A two-step synthesis for methionine-containing hydrophobic and/or aggregation-prone peptides is presented that takes advantage of the reversibility of methionine oxidation. The use of polar methionine sulfoxide as a building block in solid-phase peptide synthesis improves the synthesis quality and yields the crude peptide, with significantly improved solubility compared to the reduced species. This facilitates the otherwise often laborious peptide purification by high-performance liquid chromatography. The subsequent reduction proceeds quantitatively. This approach has been optimised with the methionine-rich Tar-DNA-binding protein 43 (307-347), but is also more generally applicable, as demonstrated by the syntheses of human calcitonin and two aggregation-prone peptides from the human prion protein.

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.

Interfacial Aromatics Mediating Cation-π Interactions with Choline-Containing Lipids Can Contribute as Much to Peripheral Protein Affinity for Membranes as Aromatics Inserted below the Phosphates

Membrane-binding interfaces of peripheral proteins are restricted to a small part of their exposed surface, so the ability to engage in strong selective interactions with membrane lipids at various depths in the interface, both below and above the phosphates, is an advantage. Driven by their hydrophobicity, aromatic amino acids preferentially partition into membrane interfaces, often below the phosphates, yet enthalpically favorable interactions with the lipid headgroups, above the phosphate plane, are likely to further stabilize high interfacial positions. Using free-energy perturbation, we calculate the energetic cost of alanine substitution for 11 interfacial aromatic amino acids from 3 peripheral proteins. We show that the involvement in cation-π interactions with the headgroups (i) increases the ΔΔG_transfer as compared with insertion at the same depth without cation-π stabilization and (ii) can contribute at least as much as deeper insertion below the phosphates, highlighting the multiple roles of aromatics in peripheral membrane protein affinity.