Insights into Membrane Protein-Lipid Interactions from Free Energy Calculations

Drug solubilization in dog intestinal fluids with and without administration of lipid-based formulations

The use of animal experiments can be minimized with computational models capable of reflecting the simulated environments. One such environment is intestinal fluid and the colloids formed in it. In this study we used molecular dynamics simulations to investigate solubilization patterns for three model drugs (carvedilol, felodipine and probucol) in dog intestinal fluid, a lipid-based formulation, and a mixture of both. We observed morphological transformations that lipids undergo due to the digestion process in the intestinal environment. Further, we evaluated the effect of bile salt concentration and observed the importance of interindividual variability. We applied two methods of estimating solubility enhancement based on the simulated data, of which one was in good qualitative agreement with the experimentally observed solubility enhancement. In addition to the computational simulations, we also measured solubility in i) aspirated dog intestinal fluid samples and ii) simulated canine intestinal fluid in the fasted state, and found there was no statistical difference between the two. Hence, a simplified dissolution medium suitable for in vitro studies provided physiologically relevant data for the systems explored. The computational protocol used in this study, coupled with in vitro studies using simulated intestinal fluids, can serve as a useful prescreening tool in the process of drug delivery strategies development.

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

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

Functional diversity among cardiolipin binding sites on the mitochondrial ADP/ATP carrier

Lipid-protein interactions play a multitude of essential roles in membrane homeostasis. Mitochondrial membranes have a unique lipid-protein environment that ensures bioenergetic efficiency. Cardiolipin (CL), the signature mitochondrial lipid, plays multiple roles in promoting oxidative phosphorylation (OXPHOS). In the inner mitochondrial membrane, the ADP/ATP carrier (AAC in yeast; adenine nucleotide translocator, ANT in mammals) exchanges ADP and ATP, enabling OXPHOS. AAC/ANT contains three tightly bound CLs, and these interactions are evolutionarily conserved. Here, we investigated the role of these buried CLs in AAC/ANT using a combination of biochemical approaches, native mass spectrometry, and molecular dynamics simulations. We introduced negatively charged mutations into each CL-binding site of yeast Aac2 and established experimentally that the mutations disrupted the CL interactions. While all mutations destabilized Aac2 tertiary structure, transport activity was impaired in a binding site-specific manner. Additionally, we determined that a disease-associated missense mutation in one CL-binding site in human ANT1 compromised its structure and transport activity, resulting in OXPHOS defects. Our findings highlight the conserved significance of CL in AAC/ANT structure and function, directly tied to specific lipid-protein interactions.

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

In silico characterization of ligand-receptor interactions for U-47700, N,N-didesmethyl-U-47700, U-50488 at mu- and kappa-opioid receptors

The increasing misuse of novel synthetic opioids (NSOs) represents a serious public health concern. In this regard, U-47700 (trans-3,4-dichloro-N-[2-(dimethylamino)cyclohexyl]-N-methylbenzamide) and related "U-compounds" emerged on recreational drug markets as synthetic substitutes for illicit heroin and constituents of counterfeit pain medications. While the pharmacology of U-compounds has been investigated using in vitro and in vivo methods, there is still a lack of understanding about the details of ligand-receptor interactions at the molecular level. To this end, we have developed a molecular modeling protocol based on docking and molecular dynamics simulations to assess the nature of ligand-receptor interactions for U-47700, N,N-didesmethyl U-47700, and U-50488 at the mu-opioid receptor (MOR) and kappa-opioid receptor (KOR). The evaluation of ligand-receptor and ligand-receptor-membrane interaction energies enabled the identification of subtle conformational shifts in the receptors induced by ligand binding. Interestingly, the removal of two key methyl groups from U-47700, to form N,N-didesmethyl U-47700, caused a loss of hydrogen bond contact with tryptophan (Trp)229, which may underlie the lower interaction energy and reduced MOR affinity for the compound. Taken together, our results are consistent with the reported biological findings for U-compounds and provide a molecular basis for the MOR selectivity of U-47700 and KOR selectivity of U-50488.

The energetics and ion coupling of cholesterol transport through Patched1

Patched1 (PTCH1) is a tumor suppressor protein of the mammalian Hedgehog (HH) signaling pathway, implicated in embryogenesis and tissue homeostasis. PTCH1 inhibits the G protein-coupled receptor Smoothened (SMO) via a debated mechanism involving modulating ciliary cholesterol accessibility. Using extensive molecular dynamics simulations and free energy calculations to evaluate cholesterol transport through PTCH1, we find an energetic barrier of ~15 to 20 kilojoule per mole for cholesterol export. In silico data are coupled to in vivo biochemical assays of PTCH1 mutants to probe coupling between cation binding sites, transmembrane motions, and PTCH1 activity. Using complementary simulations of Dispatched1, we find that transition between "inward-open" and solvent "occluded" states is accompanied by Na+-induced pinching of intracellular helical segments. Thus, our findings illuminate the energetics and ion coupling stoichiometries of PTCH1 transport mechanisms, whereby one to three Na+ or two to three K+ couple to cholesterol export, and provide the first molecular description of transitions between distinct transport states.

Molecular Biophysics of Class A G Protein Coupled Receptors-Lipids Interactome at a Glance-Highlights from the A2A Adenosine Receptor

G protein-coupled receptors (GPCRs) are embedded in phospholipid membrane bilayers with cholesterol representing 34% of the total lipid content in mammalian plasma membranes. Membrane lipids interact with GPCRs structures and modulate their function and drug-stimulated signaling through conformational selection. It has been shown that anionic phospholipids form strong interactions between positively charged residues in the G protein and the TM5-TM6-TM 7 cytoplasmic interface of class A GPCRs stabilizing the signaling GPCR-G complex. Cholesterol with a high content in plasma membranes can be identified in more specific sites in the transmembrane region of GPCRs, such as the Cholesterol Consensus Motif (CCM) and Cholesterol Recognition Amino Acid Consensus (CRAC) motifs and other receptor dependent and receptor state dependent sites. Experimental biophysical methods, atomistic (AA) MD simulations and coarse-grained (CG) molecular dynamics simulations have been applied to investigate these interactions. We emphasized here the impact of phosphatidyl inositol-4,5-bisphosphate (PtdIns(4,5)P₂ or PIP₂), a minor phospholipid component and of cholesterol on the function-related conformational equilibria of the human A_2A adenosine receptor (A_2AR), a representative receptor in class A GPCR. Several GPCRs of class A interacted with PIP₂ and cholesterol and in many cases the mechanism of the modulation of their function remains unknown. This review provides a helpful comprehensive overview for biophysics that enter the field of GPCRs-lipid systems.

Conserved cardiolipin-mitochondrial ADP/ATP carrier interactions assume distinct structural and functional roles that are clinically relevant

The mitochondrial phospholipid cardiolipin (CL) promotes bioenergetics via oxidative phosphorylation (OXPHOS). Three tightly bound CLs are evolutionarily conserved in the ADP/ATP carrier (AAC in yeast; adenine nucleotide translocator, ANT in mammals) which resides in the inner mitochondrial membrane and exchanges ADP and ATP to enable OXPHOS. Here, we investigated the role of these buried CLs in the carrier using yeast Aac2 as a model. We introduced negatively charged mutations into each CL-binding site of Aac2 to disrupt the CL interactions via electrostatic repulsion. While all mutations disturbing the CL-protein interaction destabilized Aac2 monomeric structure, transport activity was impaired in a pocket-specific manner. Finally, we determined that a disease-associated missense mutation in one CL-binding site in ANT1 compromised its structure and transport activity, resulting in OXPHOS defects. Our findings highlight the conserved significance of CL in AAC/ANT structure and function, directly tied to specific lipid-protein interactions.

UNILIPID, a Methodology for Energetically Accurate Prediction of Protein Insertion into Implicit Membranes of Arbitrary Shape

The insertion of proteins into membranes is crucial for understanding their function in many biological processes. In this work, we present UNILIPID, a universal implicit lipid-protein description as a methodology for dealing with implicit membranes. UNILIPID is independent of the scale of representation and can be applied at the level of all atoms, coarse-grained particles down to the level of a single bead per amino acid. We provide example implementations for these scales and demonstrate the versatility of our approach by accurately reflecting the free energy of transfer for each amino acid. In addition to single membranes, we describe the analytical implementation of double membranes and show that UNILIPID is well suited for modeling at multiple scales. We generalize to membranes of arbitrary shape. With UNILIPID, we provide a methodological framework for a simple and general parameterization tuned to reproduce a selected reference hydrophobicity scale. The software we provide along with the methodological description is optimized for specific user features such as real-time response, live visual analysis, and virtual reality experiences.

The Energetics and Ion Coupling of Cholesterol Transport Through Patched1

Patched1 (PTCH1) is the principal tumour suppressor protein of the mammalian Hedgehog (HH) signalling pathway, implicated in embryogenesis and tissue homeostasis. PTCH1 inhibits the Class F G protein-coupled receptor Smoothened (SMO) via a debated mechanism involving modulating accessible cholesterol levels within ciliary membranes. Using extensive molecular dynamics (MD) simulations and free energy calculations to evaluate cholesterol transport through PTCH1, we find an energetic barrier of ~15-20 kJ mol -1 for cholesterol export. In simulations we identify cation binding sites within the PTCH1 transmembrane domain (TMD) which may provide the energetic impetus for cholesterol transport. In silico data are coupled to in vivo biochemical assays of PTCH1 mutants to probe coupling between transmembrane motions and PTCH1 activity. Using complementary simulations of Dispatched1 (DISP1) we find that transition between 'inward-open' and solvent 'occluded' states is accompanied by Na + induced pinching of intracellular helical segments. Thus, our findings illuminate the energetics and ion-coupling stoichiometries of PTCH1 transport mechanisms, whereby 1-3 Na + or 2-3 K + couple to cholesterol export, and provide the first molecular description of transitions between distinct transport states.

Comparative Study of Receptor-, Receptor State-, and Membrane-Dependent Cholesterol Binding Sites in A2A and A1 Adenosine Receptors Using Coarse-Grained Molecular Dynamics Simulations

We used coarse-grained molecular dynamics (CG MD) simulations to study protein-cholesterol interactions for different activation states of the A_2A adenosine receptor (A_2AR) and the A₁ adenosine receptor (A₁R) and predict new cholesterol binding sites indicating amino acid residues with a high residence time in three biologically relevant membranes. Compared to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)-cholesterol and POPC-phosphatidylinositol-bisphosphate (PIP₂)-cholesterol, the plasma mimetic membrane best described the cholesterol binding sites previously detected for the inactive state of A_2AR and revealed the binding sites with long-lasting amino acid residues. We observed that using the plasma mimetic membrane and plotting residues with cholesterol residence time ≥2 μs, our CG MD simulations captured most obviously the cholesterol-protein interactions. For the inactive A_2AR, we identified one more binding site in which cholesterol is bound to residues with a long residence time compared to the previously detected, for the active A₁R, three binding sites, and for the inactive A₁R, two binding sites. We calculated that for the active states, cholesterol binds to residues with a much longer residence time compared to the inactive state for both A_2AR and A₁R. The stability of the identified binding sites to A₁R or A_2AR with CG MD simulations was additionally investigated with potential of mean force calculations using umbrella sampling. We observed that the binding sites with residues to which cholesterol has a long residence time in A_2AR have shallow binding free energy minima compared to the related binding sites in A₁R, suggesting a stronger binding for cholesterol to A₁R. The differences in binding sites in which cholesterol is stabilized and interacts with residues with a long residence time between active and inactive states of A₁R and A_2AR can be important for differences in functional activity and orthosteric agonist or antagonist affinity and can be used for the design of allosteric modulators, which can bind through lipid pathways. We observed a stronger binding for cholesterol to A₁R (i.e., generally higher association rates) compared to A_2AR, which remains to be demonstrated. For the active states, cholesterol binds to residues with much longer residence times compared to the inactive state for both A_2AR and A₁R. Taken together, binding sites of active A₁R may be considered as promising allosteric targets.

Supramolecular organization and dynamics of mannosylated phosphatidylinositol lipids in the mycobacterial plasma membrane

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis (TB), a disease that claims ~1.6 million lives annually. The current treatment regime is long and expensive, and missed doses contribute to drug resistance. Therefore, development of new anti-TB drugs remains one of the highest public health priorities. Mtb has evolved a complex cell envelope that represents a formidable barrier to antibiotics. The Mtb cell envelop consists of four distinct layers enriched for Mtb specific lipids and glycans. Although the outer membrane, comprised of mycolic acid esters, has been extensively studied, less is known about the plasma membrane, which also plays a critical role in impacting antibiotic efficacy. The Mtb plasma membrane has a unique lipid composition, with mannosylated phosphatidylinositol lipids (phosphatidyl-myoinositol mannosides, PIMs) comprising more than 50% of the lipids. However, the role of PIMs in the structure and function of the membrane remains elusive. Here, we used multiscale molecular dynamics (MD) simulations to understand the structure-function relationship of the PIM lipid family and decipher how they self-organize to shape the biophysical properties of mycobacterial plasma membranes. We assess both symmetric and asymmetric assemblies of the Mtb plasma membrane and compare this with residue distributions of Mtb integral membrane protein structures. To further validate the model, we tested known anti-TB drugs and demonstrated that our models agree with experimental results. Thus, our work sheds new light on the organization of the mycobacterial plasma membrane. This paves the way for future studies on antibiotic development and understanding Mtb membrane protein function.

Importance of molecular dynamics equilibrium protocol on protein-lipid interaction near channel pore

Multiscale molecular dynamics simulations using Martini coarse-grained (CG) and all-atom (AA) force fields are commonly used in membrane protein studies. In particular, reverse mapping an equilibrated CG model to an AA model offers an efficient way for preparing large membrane protein systems with complex protein shapes and lipid compositions. Here, we report that this hybrid CG-equilibrium-AA-production protocol may artificially increase lipid density and decrease hydration in ion channel pores walled with transmembrane gaps. To understand the origin of this conundrum, we conducted replicas of CG, AA, and CG reverse-mapped AA simulations of the pore domain of the mechanosensitive Piezo1 channel in a nonconducting conformation. Lipid/water density analysis and free energy calculations reveal that the lack of initial pore hydration allows excessive lipids to enter the upper pore lumen through gaps between pore helices during CG simulation. Due to the mismatch between CG and AA lipid kinetics, these pore lipids remain trapped in the subsequent AA simulations, despite unfavorable binding free energy. We tested several CG equilibrium protocols and found that a protocol restraining the whole lipid produces pore hydration consistent with AA results, thus eliminating this artifact for further studies of lipid gating and protein-lipid interactions.

Cardiolipin, and not monolysocardiolipin, preferentially binds to the interface of complexes III and IV

The mitochondrial electron transport chain comprises a series of protein complexes embedded in the inner mitochondrial membrane that generate a proton motive force via oxidative phosphorylation, ultimately generating ATP. These protein complexes can oligomerize to form larger structures called supercomplexes. Cardiolipin (CL), a conical lipid, unique within eukaryotes to the inner mitochondrial membrane, has proven essential in maintaining the stability and function of supercomplexes. Monolysocardiolipin (MLCL) is a CL variant that accumulates in people with Barth syndrome (BTHS). BTHS is caused by defects in CL biosynthesis and characterised by abnormal mitochondrial bioenergetics and destabilised supercomplexes. However, the mechanisms by which MLCL causes pathogenesis remain unclear. Here, multiscale molecular dynamics characterise the interactions of CL and MLCL with yeast and mammalian mitochondrial supercomplexes containing complex III (CIII) and complex IV (CIV). Coarse-grained simulations reveal that both CL and MLCL bind to sites at the interface between CIII and CIV of the supercomplex. Free energy perturbation calculations show that MLCL interaction is weaker than that of CL and suggest that interaction with CIV drives this difference. Atomistic contact analyses show that, although interaction with CIII is similar for CL and MLCL, CIV makes more contacts with CL than MLCL, demonstrating that CL is a more successful "glue" between the two complexes. Simulations of the human CIII₂CIV supercomplex show that this interface site is maintained between species. Our study suggests that MLCL accumulation in people with BTHS disrupts supercomplex stability by formation of relatively weak interactions at the interface lipid binding site.

Exploring the binding kinetics and behaviors of self-aggregated beta-amyloid oligomers to phase-separated lipid rafts with or without ganglioside-clusters

Lipid binding kinetics and energetics of self-aggregated and disordered beta-amyloid oligomers of various sizes, from solution to lipid raft surfaces, were investigated using MD simulations. Our systems include small (monomers to tetramers) and larger (octamers and dodecamers) oligomers. Our lipid rafts contain saturated and unsaturated phosphatidylcholine (PC), cholesterol, and with or without asymmetrically distributed monosialotetrahexosylganglioside (GM1). All rafts exhibited dynamic and structurally diversified domains including liquid-ordered (Lo), liquid-disordered (Ld), and interfacial Lod domains. For rafts without GM1, all oligomers bound to the Lod domain. For GM1-containing rafts, all small oligomers and most larger oligomers bound specifically to the GM1-clusters embedded in the Lo domain. Lipid-protein binding energies followed an order of GM1 >> unsaturated PC > saturated PC > cholesterol for all rafts. In addition, protein-induced membrane structural disruption increased progressively with the size of the oligomer for the annular lipids surrounding the membrane-bound protein in non-GM1-containing rafts. We propose that the tight binding of beta-amyloid oligomers to the GM1-clusters and the structural perturbation of lipids surrounding the membrane-bound proteins at the Lod domain are early molecular events of the beta-amyloid aggregation process on neuronal membrane surfaces that trigger the onset of Alzheimer's.

Probing effects of the SARS-CoV-2 E protein on membrane curvature and intracellular calcium

SARS-CoV-2 contains four structural proteins in its genome. These proteins aid in the assembly and budding of new virions at the ER-Golgi intermediate compartment (ERGIC). Current fundamental research efforts largely focus on one of these proteins – the spike (S) protein. Since successful antiviral therapies are likely to target multiple viral components, there is considerable interest in understanding the biophysical role of its other structural proteins, in particular structural membrane proteins. Here, we have focused our efforts on the characterization of the full-length envelope (E) protein from SARS-CoV-2, combining experimental and computational approaches. Recombinant expression of the full-length E protein from SARS-CoV-2 reveals that this membrane protein is capable of independent multimerization, possibly as a tetrameric or smaller species. Fluorescence microscopy shows that the protein localizes intracellularly, and coarse-grained MD simulations indicate it causes bending of the surrounding lipid bilayer, corroborating a potential role for the E protein in viral budding. Although we did not find robust electrophysiological evidence of ion-channel activity, cells transfected with the E protein exhibited reduced intracellular Ca²⁺, which may further promote viral replication. However, our atomistic MD simulations revealed that previous NMR structures are relatively unstable, and result in models incapable of ion conduction. Our study highlights the importance of using high-resolution structural data obtained from a full-length protein to gain detailed molecular insights, and eventually permitting virtual drug screening.

The guidance and adhesion protein FLRT2 dimerizes in cis via dual small-X3-small transmembrane motifs

Fibronectin Leucine-rich Repeat Transmembrane (FLRT 1-3) proteins are a family of broadly expressed single-spanning transmembrane receptors that play key roles in development. Their extracellular domains mediate homotypic cell-cell adhesion and heterotypic protein interactions with other receptors to regulate cell adhesion and guidance. These in trans FLRT interactions determine the formation of signaling complexes of varying complexity and function. Whether FLRTs also interact at the surface of the same cell, in cis, remains unknown. Here, molecular dynamics simulations reveal two dimerization motifs in the FLRT2 transmembrane helix. Single particle tracking experiments show that these Small-X₃-Small motifs synergize with a third dimerization motif encoded in the extracellular domain to permit the cis association and co-diffusion patterns of FLRT2 receptors on cells. These results may point to a competitive switching mechanism between in cis and in trans interactions, which suggests that homotypic FLRT interaction mirrors the functionalities of classic adhesion molecules.

On Computing Equilibrium Binding Constants for Protein-Protein Association in Membranes

Protein association in lipid membranes is fundamental to membrane protein function and of great biomedical relevance. All-atom and coarse-grained models have been extensively used to understand the protein-protein interactions in the membrane and to compute equilibrium association constants. However, slow translational and rotational diffusion of protein in membrane presents challenges to the effective sampling of conformations defining the ensembles of free and bound states contributing to the association equilibrium and the free energy of dimerization. We revisit the homodimerization equilibrium of the TM region of glycophorin A. Conformational sampling is performed using umbrella sampling along previously proposed one-dimensional collective variables and compared with sampling over a two-dimensional collective variable space using the MARTINI v2.2 force field. We demonstrate that the one-dimensional collective variables suffer from restricted sampling of the native homodimer conformations leading to a biased free energy landscape. Conversely, simulations along the two-dimensional collective variable effectively characterize the thermodynamically relevant native and non-native interactions contributing to the association equilibrium. These results demonstrate the challenges associated with accurately characterizing binding equilibria when multiple poses contribute to the bound state ensemble.

Transition between conformational states of the TREK-1 K2P channel promoted by interaction with PIP2

Members of the TREK family of two-pore domain (K2P) potassium channels are highly sensitive to regulation by membrane lipids, including phosphatidylinositol-4,5-bisphosphate (PIP₂). Previous studies have demonstrated that PIP₂ increases TREK1 channel activity, however, the mechanistic understanding of the conformational transitions induced by PIP₂ remain unclear. Here, we used coarse-grained molecular dynamics (CG-MD) and atomistic MD simulations to model the PIP₂ binding site on both the up and down state conformations of TREK-1. We also calculated the free energy of PIP₂ binding relative to other anionic phospholipids in both conformational states using potential of mean force (PMF) and free energy perturbation (FEP) calculations. Our results identify state-dependent binding of PIP₂ to sites involving the proximal C-terminus and we show that PIP₂ promotes a conformational transition from a down state towards an intermediate that resembles the up state. These results are consistent with functional data for PIP₂ regulation and together provide evidence for a structural mechanism of TREK-1 channel activation by phosphoinositides.

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.

Interaction With the Lipid Membrane Influences Fentanyl Pharmacology

Overdose deaths from fentanyl have reached epidemic proportions in the USA and are increasing worldwide. Fentanyl is a potent opioid agonist that is less well reversed by naloxone than morphine. Due to fentanyl’s high lipophilicity and elongated structure we hypothesised that its unusual pharmacology may be explained by its interactions with the lipid membrane on route to binding to the μ-opioid receptor (MOPr). Through coarse-grained molecular dynamics simulations, electrophysiological recordings and cell signalling assays, we determined how fentanyl and morphine access the orthosteric pocket of MOPr. Morphine accesses MOPr via the aqueous pathway; first binding to an extracellular vestibule, then diffusing into the orthosteric pocket. In contrast, fentanyl may take a novel route; first partitioning into the membrane, before accessing the orthosteric site by diffusing through a ligand-induced gap between the transmembrane helices. In electrophysiological recordings fentanyl-induced currents returned after washout, suggesting fentanyl deposits in the lipid membrane. However, mutation of residues forming the potential MOPr transmembrane access site did not alter fentanyl’s pharmacological profile in vitro. A high local concentration of fentanyl in the lipid membrane, possibly in combination with a novel lipophilic binding route, may explain the high potency and lower susceptibility of fentanyl to reversal by naloxone.

PyLipID: A Python Package for Analysis of Protein-Lipid Interactions from Molecular Dynamics Simulations

Lipids play important modulatory and structural roles for membrane proteins. Molecular dynamics simulations are frequently used to provide insights into the nature of these protein-lipid interactions. Systematic comparative analysis requires tools that provide algorithms for objective assessment of such interactions. We introduce PyLipID, a Python package for the identification and characterization of specific lipid interactions and binding sites on membrane proteins from molecular dynamics simulations. PyLipID uses a community analysis approach for binding site detection, calculating lipid residence times for both the individual protein residues and the detected binding sites. To assist structural analysis, PyLipID produces representative bound lipid poses from simulation data, using a density-based scoring function. To estimate residue contacts robustly, PyLipID uses a dual-cutoff scheme to differentiate between lipid conformational rearrangements while bound from full dissociation events. In addition to the characterization of protein-lipid interactions, PyLipID is applicable to analysis of the interactions of membrane proteins with other ligands. By combining automated analysis, efficient algorithms, and open-source distribution, PyLipID facilitates the systematic analysis of lipid interactions from large simulation data sets of multiple species of membrane proteins.

Cell-Free Expression to Probe Co-Translational Insertion of an Alpha Helical Membrane Protein

The majority of alpha helical membrane proteins fold co-translationally during their synthesis on the ribosome. In contrast, most mechanistic folding studies address refolding of full-length proteins from artificially induced denatured states that are far removed from the natural co-translational process. Cell-free translation of membrane proteins is emerging as a useful tool to address folding during translation by a ribosome. We summarise the benefits of this approach and show how it can be successfully extended to a membrane protein with a complex topology. The bacterial leucine transporter, LeuT can be synthesised and inserted into lipid membranes using a variety of in vitro transcription translation systems. Unlike major facilitator superfamily transporters, where changes in lipids can optimise the amount of correctly inserted protein, LeuT insertion yields are much less dependent on the lipid composition. The presence of a bacterial translocon either in native membrane extracts or in reconstituted membranes also has little influence on the yield of LeuT incorporated into the lipid membrane, except at high reconstitution concentrations. LeuT is considered a paradigm for neurotransmitter transporters and possesses a knotted structure that is characteristic of this transporter family. This work provides a method in which to probe the formation of a protein as the polypeptide chain is being synthesised on a ribosome and inserting into lipids. We show that in comparison with the simpler major facilitator transporter structures, LeuT inserts less efficiently into membranes when synthesised cell-free, suggesting that more of the protein aggregates, likely as a result of the challenging formation of the knotted topology in the membrane.

Lipid–Protein Interactions in Plasma Membrane Organization and Function

Lipid–protein interactions in cells are involved in various biological processes, including metabolism, trafficking, signaling, host–pathogen interactions, and transmembrane transport. At the plasma membrane, lipid–protein interactions play major roles in membrane organization and function. Several membrane proteins have motifs for specific lipid binding, which modulate protein conformation and consequent function. In addition to such specific lipid–protein interactions, protein function can be regulated by the dynamic, collective behavior of lipids in membranes. Emerging analytical, biochemical, and computational technologies allow us to study the influence of specific lipid–protein interactions, as well as the collective behavior of membranes on protein function. In this article, we review the recent literature on lipid–protein interactions with a specific focus on the current state-of-the-art technologies that enable novel insights into these interactions.

Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Entropy in the Molecular Recognition of Membrane Protein-Lipid Interactions

Understanding the molecular driving forces that underlie membrane protein-lipid interactions requires the characterization of their binding thermodynamics. Here, we employ variable-temperature native mass spectrometry to determine the thermodynamics of lipid binding events to the human G-protein-gated inward rectifier potassium channel, Kir3.2. The channel displays distinct thermodynamic strategies to engage phosphatidylinositol (PI) and phosphorylated forms thereof. The addition of a 4'-phosphate to PI results in an increase in favorable entropy. PI with two or more phosphates exhibits more complex binding, where lipids appear to bind two nonidentical sites on Kir3.2. Remarkably, the interaction of 4,5-bisphosphate PI with Kir3.2 is solely driven by a large, favorable change in entropy. Installment of a 3'-phosphate to PI(4,5)P₂ results in an altered thermodynamic strategy. The acyl chain of the lipid has a marked impact on binding thermodynamics and, in some cases, enthalpy becomes favorable.

CG2AT2: an Enhanced Fragment-Based Approach for Serial Multi-scale Molecular Dynamics Simulations

Coarse-grained molecular dynamics provides a means for simulating the assembly and interactions of macromolecular complexes at a reduced level of representation, thereby allowing both longer timescale and larger sized simulations. Here, we describe an enhanced fragment-based protocol for converting macromolecular complexes from coarse-grained to atomistic resolution, for further refinement and analysis. While the focus is upon systems that comprise an integral membrane protein embedded in a phospholipid bilayer, the technique is also suitable for membrane-anchored and soluble protein/nucleotide complexes. Overall, this provides a method for generating an accurate and well-equilibrated atomic-level description of a macromolecular complex. The approach is evaluated using a diverse test set of 11 system configurations of varying size and complexity. Simulations are assessed in terms of protein stereochemistry, conformational drift, lipid/protein interactions, and lipid dynamics.

Relative Affinities of Protein-Cholesterol Interactions from Equilibrium Molecular Dynamics Simulations

Specific interactions of lipids with membrane proteins contribute to protein stability and function. Multiple lipid interactions surrounding a membrane protein are often identified in molecular dynamics (MD) simulations and are, increasingly, resolved in cryo-electron microscopy (cryo-EM) densities. Determining the relative importance of specific interaction sites is aided by determination of lipid binding affinities using experimental or simulation methods. Here, we develop a method for determining protein-lipid binding affinities from equilibrium coarse-grained MD simulations using binding saturation curves, designed to mimic experimental protocols. We apply this method to directly obtain affinities for cholesterol binding to multiple sites on a range of membrane proteins and compare our results with free energies obtained from density-based equilibrium methods and with potential of mean force calculations, getting good agreement with respect to the ranking of affinities for different sites. Thus, our binding saturation method provides a robust, high-throughput alternative for determining the relative consequence of individual sites seen in, e.g., cryo-EM derived membrane protein structures surrounded by an array of ancillary lipid densities.

Recent developments in multiscale free energy simulations

Physics-based free energy simulations enable the rigorous calculation of properties, such as conformational equilibria, solvation or binding free energies. While historically most applications have occurred at the atomistic level of resolution, a range of advances in the past years make it possible now to reliably cross the temporal, spatial and theory scales for the modeling of complex systems or the efficient prediction of results at the accuracy level of expensive quantum-mechanical calculations. In this mini-review, we discuss recent methodological advances as well as opportunities opened up by the introduction of machine learning approaches, which tackle the diverse challenges across the different scales, improve the accuracy and feasibility, and push the boundaries of multiscale free energy simulations.

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.

From Bench to Biomolecular Simulation: Phospholipid Modulation of Potassium Channels

Potassium (K+) ion channels are crucial in numerous cellular processes as they hyperpolarise a cell through K+ conductance, returning a cell to its resting potential. K+ channel mutations result in multiple clinical complications such as arrhythmia, neonatal diabetes and migraines. Since 1995, the regulation of K+ channels by phospholipids has been heavily studied using a range of interdisciplinary methods such as cellular electrophysiology, structural biology and computational modelling. As a result, K+ channels are model proteins for the analysis of protein-lipid interactions. In this review, we will focus on the roles of lipids in the regulation of K+ channels, and how atomic-level structures, along with experimental techniques and molecular simulations, have helped guide our understanding of the importance of phospholipid interactions.

Recent Developments in Free Energy Calculations for Drug Discovery

The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.

Identification and assessment of cardiolipin interactions with E. coli inner membrane proteins

Integral membrane proteins are localized and/or regulated by lipids present in the surrounding bilayer. While bacteria have relatively simple membranes, there is ample evidence that many bacterial proteins bind to specific lipids, especially the anionic lipid cardiolipin. Here, we apply molecular dynamics simulations to assess lipid binding to 42 different Escherichia coli inner membrane proteins. Our data reveal an asymmetry between the membrane leaflets, with increased anionic lipid binding to the inner leaflet regions of the proteins, particularly for cardiolipin. From our simulations, we identify >700 independent cardiolipin binding sites, allowing us to identify the molecular basis of a prototypical cardiolipin binding site, which we validate against structures of bacterial proteins bound to cardiolipin. This allows us to construct a set of metrics for defining a high-affinity cardiolipin binding site on bacterial membrane proteins, paving the way for a heuristic approach to defining other protein-lipid interactions.

Building intuition for binding free energy calculations: Bound state definition, restraints, and symmetry

The theory behind computation of absolute binding free energies using explicit-solvent molecular simulations is well-established, yet somewhat complex, with counter-intuitive aspects. This leads to frequent frustration, common misconceptions, and sometimes erroneous numerical treatment. To improve this, we present the main practically relevant segments of the theory with constant reference to physical intuition. We pinpoint the role of the implicit or explicit definition of the bound state (or the binding site) to make a robust link between an experimental measurement and a computational result. We clarify the role of symmetry and discuss cases where symmetry number corrections have been misinterpreted. In particular, we argue that symmetry corrections as classically presented are a source of confusion and could be advantageously replaced by restraint free energy contributions. We establish that contrary to a common intuition, partial or missing sampling of some modes of symmetric bound states does not affect the calculated decoupling free energies. Finally, we review these questions and pitfalls in the context of a few common practical situations: binding to a symmetric receptor (equivalent binding sites), binding of a symmetric ligand (equivalent poses), and formation of a symmetric complex, in the case of homodimerization.

Structure and lipid dynamics in the maintenance of lipid asymmetry inner membrane complex of A. baumannii

Multi-resistant bacteria are a major threat in modern medicine. The gram-negative coccobacillus Acinetobacter baumannii currently leads the WHO list of pathogens in critical need for new therapeutic development. The maintenance of lipid asymmetry (MLA) protein complex is one of the core machineries that transport lipids from/to the outer membrane in gram-negative bacteria. It also contributes to broad-range antibiotic resistance in several pathogens, most prominently in A. baumannii. Nonetheless, the molecular details of its role in lipid transport has remained largely elusive. Here, we report the cryo-EM maps of the core MLA complex, MlaBDEF, from the pathogen A. baumannii, in the apo-, ATP- and ADP-bound states, revealing multiple lipid binding sites in the cytosolic and periplasmic side of the complex. Molecular dynamics simulations suggest their potential trajectory across the membrane. Collectively with the recently-reported structures of the E. coli orthologue, this data also allows us to propose a molecular mechanism of lipid transport by the MLA system.

Daniel Mann et al. describe a higher-resolution structure of the maintenance of lipid asymmetry inner membrane complex (MlaBDEF) in the Gram-negative pathogen, Acinetobacter baumannii. With this improved structural map, the authors clarify the secondary structure elements of MlaE helices and report on potential lipid dynamics by the MLA system that could inform the development of future therapeutics against A. baumannii infection.

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

Spontaneous lipid binding to the nicotinic acetylcholine receptor in a native membrane

The nicotinic acetylcholine receptor (nAChR) and other pentameric ligand-gated ion channels are native to neuronal membranes with an unusual lipid composition. While it is well-established that these receptors can be significantly modulated by lipids, the underlying mechanisms have been primarily studied in model membranes with few lipid species. Here, we use coarse-grained molecular dynamics simulation to probe specific binding of lipids in a complex quasi-neuronal membrane. We ran a total of 50 μs of simulations of a single nAChR in a membrane composed of 36 species of lipids. Competition between multiple lipid species produces a complex distribution. We find that overall, cholesterol selects for concave inter-subunit sites and polyunsaturated fatty acids select for convex M4 sites, while monounsaturated and saturated lipids are unenriched in the nAChR boundary. We propose the "density-threshold affinity" as a metric calculated from continuous density distributions, which reduces to a standard affinity in two-state binding. We find that the density-threshold affinity for M4 weakens with chain rigidity, which suggests that flexible chains may help relax packing defects caused by the conical protein shape. For any site, PE headgroups have the strongest affinity of all phospholipid headgroups, but anionic lipids still yield moderately high affinities for the M4 sites as expected. We observe cooperative effects between anionic headgroups and saturated chains at the M4 site in the inner leaflet. We also analyze affinities for individual anionic headgroups. When combined, these insights may reconcile several apparently contradictory experiments on the role of anionic phospholipids in modulating nAChR.

Lipid regulation of hERG1 channel function

The lipid regulation of mammalian ion channel function has emerged as a fundamental mechanism in the control of electrical signalling and transport specificity in various cell types. In this work, we combine molecular dynamics simulations, mutagenesis, and electrophysiology to provide mechanistic insights into how lipophilic molecules (ceramide-sphingolipid probe) alter gating kinetics and K+ currents of hERG1. We show that the sphingolipid probe induced a significant left shift of activation voltage, faster deactivation rates, and current blockade comparable to traditional hERG1 blockers. Microseconds-long MD simulations followed by experimental mutagenesis elucidated ceramide specific binding locations at the interface between the pore and voltage sensing domains. This region constitutes a unique crevice present in mammalian channels with a non-swapped topology. The combined experimental and simulation data provide evidence for ceramide-induced allosteric modulation of the channel by a conformational selection mechanism.

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.

How Do Molecular Dynamics Data Complement Static Structural Data of GPCRs

G protein-coupled receptors (GPCRs) are implicated in nearly every physiological process in the human body and therefore represent an important drug targeting class. Advances in X-ray crystallography and cryo-electron microscopy (cryo-EM) have provided multiple static structures of GPCRs in complex with various signaling partners. However, GPCR functionality is largely determined by their flexibility and ability to transition between distinct structural conformations. Due to this dynamic nature, a static snapshot does not fully explain the complexity of GPCR signal transduction. Molecular dynamics (MD) simulations offer the opportunity to simulate the structural motions of biological processes at atomic resolution. Thus, this technique can incorporate the missing information on protein flexibility into experimentally solved structures. Here, we review the contribution of MD simulations to complement static structural data and to improve our understanding of GPCR physiology and pharmacology, as well as the challenges that still need to be overcome to reach the full potential of this technique.

Monolysocardiolipin (MLCL) interactions with mitochondrial membrane proteins

Monolysocardiolipin (MLCL) is a three-tailed variant of cardiolipin (CL), the signature lipid of mitochondria. MLCL is not normally found in healthy tissue but accumulates in mitochondria of people with Barth syndrome (BTHS), with an overall increase in the MLCL:CL ratio. The reason for MLCL accumulation remains to be fully understood. The effect of MLCL build-up and decreased CL content in causing the characteristics of BTHS are also unclear. In both cases, an understanding of the nature of MLCL interaction with mitochondrial proteins will be key. Recent work has shown that MLCL associates less tightly than CL with proteins in the mitochondrial inner membrane, suggesting that MLCL accumulation is a result of CL degradation, and that the lack of MLCL–protein interactions compromises the stability of the protein-dense mitochondrial inner membrane, leading to a decrease in optimal respiration. There is some data on MLCL–protein interactions for proteins involved in the respiratory chain and in apoptosis, but there remains much to be understood regarding the nature of MLCL–protein interactions. Recent developments in structural, analytical and computational approaches mean that these investigations are now possible. Such an understanding will be key to further insights into how MLCL accumulation impacts mitochondrial membranes. In turn, these insights will help to support the development of therapies for people with BTHS and give a broader understanding of other diseases involving defective CL content.

Defining how multiple lipid species interact with inward rectifier potassium (Kir2) channels

Protein-lipid interactions are a key element of the function of many integral membrane proteins. These potential interactions should be considered alongside the complexity and diversity of membrane lipid composition. Inward rectifier potassium channel (Kir) Kir2.2 has multiple interactions with plasma membrane lipids: Phosphatidylinositol (4, 5)-bisphosphate (PIP₂) activates the channel; a secondary anionic lipid site has been identified, which augments the activation by PIP₂; and cholesterol inhibits the channel. Molecular dynamics simulations are used to characterize in molecular detail the protein-lipid interactions of Kir2.2 in a model of the complex plasma membrane. Kir2.2 has been simulated with multiple, functionally important lipid species. From our simulations we show that PIP₂ interacts most tightly at the crystallographic interaction sites, outcompeting other lipid species at this site. Phosphatidylserine (PS) interacts at the previously identified secondary anionic lipid interaction site, in a PIP₂ concentration-dependent manner. There is interplay between these anionic lipids: PS interactions are diminished when PIP₂ is not present in the membrane, underlining the need to consider multiple lipid species when investigating protein-lipid interactions.

Defining how multiple lipid species interact with inward rectifier potassium (Kir2) channels

Ion channels form pores that allow for the selective transport of ions across cell membranes, generating electrical signals in response to a variety of signals. Inward rectifier potassium (Kir) channels in particular are regulated by direct interactions with the complex mixture of lipids that are present in eukaryotic cell membranes. However, the molecular details of these concurrent lipid interactions with Kir channels are not clear and difficult to access via experimental methods. Here, we simulate the Kir2.2 channel in a complex lipid mixture to explore how anionic phospholipids and cholesterol dynamically organize around the membrane protein. In particular we demonstrate a synergy between binding interactions of different anionic phospholipid species which are known to activate Kir channels.

Protein–lipid interactions are a key element of the function of many integral membrane proteins. These potential interactions should be considered alongside the complexity and diversity of membrane lipid composition. Inward rectifier potassium channel (Kir) Kir2.2 has multiple interactions with plasma membrane lipids: Phosphatidylinositol (4, 5)-bisphosphate (PIP₂) activates the channel; a secondary anionic lipid site has been identified, which augments the activation by PIP₂; and cholesterol inhibits the channel. Molecular dynamics simulations are used to characterize in molecular detail the protein–lipid interactions of Kir2.2 in a model of the complex plasma membrane. Kir2.2 has been simulated with multiple, functionally important lipid species. From our simulations we show that PIP₂ interacts most tightly at the crystallographic interaction sites, outcompeting other lipid species at this site. Phosphatidylserine (PS) interacts at the previously identified secondary anionic lipid interaction site, in a PIP₂ concentration-dependent manner. There is interplay between these anionic lipids: PS interactions are diminished when PIP₂ is not present in the membrane, underlining the need to consider multiple lipid species when investigating protein–lipid interactions.

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.

A Mass-Spectrometry-Based Approach to Distinguish Annular and Specific Lipid Binding to Membrane Proteins

Membrane proteins engage in a variety of contacts with their surrounding lipids, but distinguishing between specifically bound lipids, and non‐specific, annular interactions is a challenging problem. Applying native mass spectrometry to three membrane protein complexes with different lipid‐binding properties, we explore the ability of detergents to compete with lipids bound in different environments. We show that lipids in annular positions on the presenilin homologue protease are subject to constant exchange with detergent. By contrast, detergent‐resistant lipids bound at the dimer interface in the leucine transporter show decreased k_off rates in molecular dynamics simulations. Turning to the lipid flippase MurJ, we find that addition of the natural substrate lipid‐II results in the formation of a 1:1 protein–lipid complex, where the lipid cannot be displaced by detergent from the highly protected active site. In summary, we distinguish annular from non‐annular lipids based on their exchange rates in solution.

Best Practices for Alchemical Free Energy Calculations [Article v1.0]

Alchemical free energy calculations are a useful tool for predicting free energy differences associated with the transfer of molecules from one environment to another. The hallmark of these methods is the use of "bridging" potential energy functions representing alchemical intermediate states that cannot exist as real chemical species. The data collected from these bridging alchemical thermodynamic states allows the efficient computation of transfer free energies (or differences in transfer free energies) with orders of magnitude less simulation time than simulating the transfer process directly. While these methods are highly flexible, care must be taken in avoiding common pitfalls to ensure that computed free energy differences can be robust and reproducible for the chosen force field, and that appropriate corrections are included to permit direct comparison with experimental data. In this paper, we review current best practices for several popular application domains of alchemical free energy calculations performed with equilibrium simulations, in particular relative and absolute small molecule binding free energy calculations to biomolecular targets.

Lipid Interactions of a Ciliary Membrane TRP Channel: Simulation and Structural Studies of Polycystin-2

Polycystin-2 (PC2) is a transient receptor potential (TRP) channel present in ciliary membranes of the kidney. PC2 shares a transmembrane fold with other TRP channels, in addition to an extracellular domain found in TRPP and TRPML channels. Using molecular dynamics (MD) simulations and cryoelectron microscopy we identify and characterize PIP₂ and cholesterol interactions with PC2. PC2 is revealed to have a PIP binding site close to the equivalent vanilloid/lipid binding site in the TRPV1 channel. A 3.0-Å structure reveals a binding site for cholesterol on PC2. Cholesterol interactions with the channel at this site are characterized by MD simulations. The two classes of lipid binding sites are compared with sites observed in other TRPs and in Kv channels. These findings suggest PC2, in common with other ion channels, may be modulated by both PIPs and cholesterol, and position PC2 within an emerging model of the roles of lipids in the regulation and organization of ciliary membranes.

•

Lipid interactions of PC2 channels have been explored by MD simulation and cryo-EM

•

PIP₂ binds to a site corresponding to the vanilloid/lipid binding site of TRPV1

•

Cholesterol binds between the S3 and S4 helices and S6 of the adjacent subunit

•

PC2, in common with other channels, may be modulated by PIPs and cholesterol