Associations between pre-operative cholesterol levels with long-term survival after colorectal cancer surgery: a nationwide propensity score-matched cohort study

PURPOSE

Altered lipid metabolism frequently occurs in patients with solid cancers and dyslipidemia has been associated with poorer outcomes in patients with colorectal cancer. This study sought to investigate whether cholesterol levels are associated with clinical outcomes and can serve as survival predictors.

METHODS

We conducted a retrospective cohort study with Danish patients diagnosed with colorectal cancer who had surgery with curative intent for UICC stages I to III between 2015 and 2020. Using propensity score adjustment, we matched patients in a 1:1 ratio to examine the impact of total cholesterol (TC) > 4 mmol/L vs. ≤ 4 mmol/L within 365 days prior to surgery on overall survival (OS) and disease-free survival (DFS).

RESULTS

A total of 3443 patients were included in the study. Median follow-up time was 3.8 years. Following propensity score matching, 1572 patients were included in the main analysis. There was no statistically significant difference in OS or DFS between patients with TC > 4 mmol/L compared with TC ≤ 4 mmol/L (HR: 0.82, 95% CI, 0.65-1.03, HR: 0.87, 95% CI, 0.68-1.12, respectively.). A subgroup analysis investigating TC > 4 mmol/L as well as low-density lipoprotein (LDL) > 3 mmol/L found a significant correlation with OS (HR: 0.74, 95% CI, 0.54-0.99).

CONCLUSION

TC levels alone were not associated with OS or DFS in patients with colorectal cancer. Interestingly, higher TC and LDL levels were linked to better overall survival, suggesting the need for further exploration of cholesterol's role in colorectal cancer.

TRIAL REGISTRATION

Not applicable.

PIP2 inhibits pore opening of the cyclic nucleotide-gated channel SthK

The signaling lipid phosphatidylinositol-4,5-bisphosphate (PIP2) regulates many ion channels. It inhibits eukaryotic cyclic nucleotide-gated (CNG) channels while activating their relatives, the hyperpolarization-activated and cyclic nucleotide-modulated (HCN) channels. The prokaryotic SthK channel from Spirochaeta thermophila shares features with CNG and HCN channels and is an established model for this channel family. Here, we show SthK activity is inhibited by PIP2. A cryo-EM structure of SthK in nanodiscs reveals a PIP2-fitting density coordinated by arginine and lysine residues from the S4 helix and the C-linker, located between voltage-sensing and pore domains of adjacent subunits. Mutation of two arginine residues weakens PIP2 inhibition with the double mutant displaying insensitivity to PIP2. We propose that PIP2 inhibits SthK by gluing S4 and S6 together, stabilizing a resting channel conformation. The PIP2 binding site is partially conserved in CNG channels suggesting the possibility of a similar inhibition mechanism in the eukaryotic homologs.

Oleic acid released by sensory neurons inhibits TRPV1-mediated thermal hypersensitivity via GPR40

Noxious stimuli activate nociceptive sensory neurons, causing action potential firing and the release of diverse signaling molecules. Several peptides have already been identified to be released by sensory neurons and shown to modulate inflammatory responses and inflammatory pain. However, it is still unclear whether lipid mediators can be released upon sensory neuron activation to modulate intercellular communication. Here, we analyzed the lipid secretome of capsaicin-stimulated nociceptive neurons with LC-HRMS, revealing that oleic acid is strongly released from sensory neurons by capsaicin. We further demonstrated that oleic acid inhibits capsaicin-induced calcium transients in sensory neurons and reverses bradykinin-induced TRPV1 sensitization by a calcineurin (CaN) and GPR40 (FFAR1) dependent pathway. Additionally, oleic acid alleviated zymosan-mediated thermal hypersensitivity via the GPR40, suggesting that the capsaicin-mediated oleic acid release from sensory neurons acts as a protective and feedback mechanism, preventing sensory neurons from nociceptive overstimulation via the GPR40/CaN/TRPV1-axis.

Exploring new avenues of health protection: plant-derived nanovesicles reshape microbial communities

Symbiotic microbial communities are crucial for human health, and dysbiosis is associated with various diseases. Plant-derived nanovesicles (PDNVs) have a lipid bilayer structure and contain lipids, metabolites, proteins, and RNA. They offer unique advantages in regulating microbial community homeostasis and treating diseases related to dysbiosis compared to traditional drugs. On the one hand, lipids on PDNVs serve as the primary substances that mediate specific recognition and uptake by bacteria. On the other hand, due to the multifactorial nature of PDNVs, they have the potential to enhance growth and survival of beneficial bacterial while simultaneously reducing the pathogenicity of harmful bacteria. In addition, PDNVs have the capacity to modulate bacterial metabolism, thus facilitating the establishment of a harmonious microbial equilibrium and promoting stability within the microbiota. These remarkable attributes make PDNVs a promising therapeutic approach for various conditions, including periodontitis, inflammatory bowel disease, and skin infection diseases. However, challenges such as consistency, isolation methods, and storage need to be addressed before clinical application. This review aims to explore the value of PDNVs in regulating microbial community homeostasis and provide recommendations for their use as novel therapeutic agents for health protection.

Exploring the effects of epigallocatechin gallate on lipid metabolism in the rat steatotic liver during normothermic machine perfusion: Insights from lipidomics and RNA sequencing

BACKGROUND

Hepatic steatosis is the leading cause of discarded liver grafts. Defatting steatotic liver grafts using drug combinations during ex vivo normothermic machine perfusion (NMP) has been reported. However, the effectiveness of NMP in reducing fat content using epigallocatechin gallate (EGCG) as a single defatting agent and its effect on lipid metabolism are poorly investigated.

METHODS

In this study, an NMP system was set up to perfuse a steatotic liver from a rat model with 10 mM EGCG. Livers without EGCG served as NMP controls, whereas static cold-preserved livers in the University of Wisconsin medium were used as static cold storage controls. Liver enzyme, reactive oxygen species (ROS), histology, and lipid content assessments were conducted post-perfusion, complemented by lipidomics, RNA sequencing, and western blotting to determine the lipid metabolism changes.

RESULTS

EGCG during NMP reduced hepatocellular injury markers and defatted steatotic liver grafts. Additionally, we observed a significant increase in triglyceride (TG) content in the perfusate post-NMP in the NMP + EGCG group, suggesting TG output from the liver. Furthermore, lipidomics analysis revealed that EGCG primarily affected metabolites involved in glycerophospholipid (GP) and glycerolipid (GL) metabolism. Further, the RNA sequencing indicated the modulation of these metabolic pathways via ECGC, which was associated with the downregulated Lpin1 and Gpat3 expression.

CONCLUSIONS

EGCG defats steatotic livers as a single defatting agent during NMP by promoting GL and GP metabolism via decreasing Lpin1 and Agpat9 levels.

FFAR4 activation inhibits lung adenocarcinoma via blocking respiratory chain complex assembly associated mitochondrial metabolism

Despite notable advancements in the investigation and management of lung adenocarcinoma (LUAD), the mortality rate for individuals afflicted with LUAD remains elevated, and attaining an accurate prognosis is challenging. LUAD exhibits intricate genetic and environmental components, and it is plausible that free fatty acid receptors (FFARs) may bridge the genetic and dietary aspects. The objective of this study is to ascertain whether a correlation exists between FFAR4, which functions as the primary receptor for dietary fatty acids, and various characteristics of LUAD, while also delving into the potential underlying mechanism. The findings of this study indicate a decrease in FFAR4 expression in LUAD, with a positive correlation (P < 0.01) between FFAR4 levels and overall patient survival (OS). Receiver operating characteristic (ROC) curve analysis demonstrated a significant diagnostic value [area under the curve (AUC) of 0.933] associated with FFAR4 expression. Functional investigations revealed that the FFAR4-specific agonist (TUG891) effectively suppressed cell proliferation and induced cell cycle arrest. Furthermore, FFAR4 activation resulted in significant metabolic shifts, including a decrease in oxygen consumption rate (OCR) and an increase in extracellular acidification rate (ECAR) in A549 cells. In detail, the activation of FFAR4 has been observed to impact the assembly process of the mitochondrial respiratory chain complex and the malate-aspartate shuttle process, resulting in a decrease in the transition of NAD+ to NADH and the inhibition of LUAD. These discoveries reveal a previously unrecognized function of FFAR4 in the negative regulation of mitochondrial metabolism and the inhibition of LUAD, indicating its potential as a promising therapeutic target for the treatment and diagnosis of LUAD.

Dynamic lipid interactions in the plasma membrane Na+,K+-ATPase

The function of ion-transporting Na+,K+-ATPases depends on the surrounding lipid environment in biological membranes. Two established lipid-interaction sites A and B within the transmembrane domain have been observed to induce protein activation and stabilization, respectively. In addition, lipid-mediated inhibition has been assigned to a site C, but with the exact location not experimentally confirmed. Also, possible effects on lipid interactions by disease mutants dwelling in the membrane-protein interface remain relatively uncharacterized. We simulated human Na+,K+-ATPase α1β1FXYD homology models in E1 and E2 states in an asymmetric, multicomponent plasma membrane to determine both wild-type and disease mutant lipid-protein interactions. The simulated wild-type lipid interactions at the established sites A and B were in agreement with experimental results thereby confirming the membrane-protein model system. The less well-characterized, proposed inhibitory site C was dominated by lipids lacking inhibitory properties. Instead, two sites hosting inhibitory lipids were identified at the extracellular side and also a cytoplasmic CHL-binding site that provide putative alternative locations of Na+,K+-ATPase inhibition. Three disease mutations, Leu302Arg, Glu840Arg and Met859Arg resided in the lipid-protein interface and caused drastic changes in the lipid interactions. The simulation results show that lipid interactions to the human Na+,K+-ATPase α1β1FXYD protein in the plasma membrane are highly state-dependent and can be disturbed by disease mutations located in the lipid interface, which can open up for new venues to understand genetic disorders.

Structural Studies of Monounsaturated and ω-3 Polyunsaturated Free Fatty Acids in Solution with the Combined Use οf NMR and DFT Calculations-Comparison with the Liquid State

Molecular structures, in chloroform and DMSO solution, of the free fatty acids (FFAs) caproleic acid, oleic acid, α-linolenic acid, eicosapentanoic acid (EPA) and docosahexaenoic acid (DHA) are reported with the combined use of NMR and DFT calculations. Variable temperature and concentration chemical shifts of the COOH protons, transient 1D NOE experiments and DFT calculations demonstrate the major contribution of low molecular weight aggregates of dimerized fatty acids through intermolecular hydrogen bond interactions of the carboxylic groups, with parallel and antiparallel interdigitated structures even at the low concentration of 20 mM in CDCl3. For the dimeric DHA, a structural model of an intermolecular hydrogen bond through carboxylic groups and an intermolecular hydrogen bond between the carboxylic group of one molecule and the ω-3 double bond of a second molecule is shown to play a role. In DMSO-d6 solution, NMR and DFT studies show that the carboxylic groups form strong intermolecular hydrogen bond interactions with a single discrete solvation molecule of DMSO. These solvation species form parallel and antiparallel interdigitated structures of low molecular weight, as in chloroform solution. This structural motif, therefore, is an intrinsic property of the FFAs, which is not strongly affected by the length and degree of unsaturation of the chain and the hydrogen bond ability of the solvent.

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.

Molecular Mechanisms Underlying Caveolin-1 Mediated Membrane Curvature

Caveolin-1 is one of the main protein components of caveolae that acts as a mechanosensor at the cell membrane. The interactions of caveolin-1 with membranes have been shown to lead to complex effects such as curvature and the clustering of specific lipids. Here, we review the emerging concepts on the molecular interactions of caveolin-1, with a focus on insights from coarse-grain molecular dynamics simulations. Consensus structural models of caveolin-1 report a helix-turn-helix core motif with flanking domains of higher disorder that could be membrane composition dependent. Caveolin-1 appears to be mainly surface-bound and does not embed very deep in the membrane to which it is bound. The most interesting aspect of caveolin-1 membrane binding is the interplay of cholesterol clustering and membrane curvature. Although cholesterol has been reported to cluster in the vicinity of caveolin-1 by several approaches, simulations show that the clustering is maximal in membrane leaflet opposing the surface-bound caveolin-1. The intrinsic negative curvature of cholesterol appears to stabilize the negative curvature in the opposing leaflet. In fact, the simulations show that blocking cholesterol clustering (through artificial position restraints) blocks membrane curvature, and vice versa. Concomitant with cholesterol clustering is sphingomyelin clustering, again in the opposing leaflet, but in a concentration-dependent manner. The differential stress due to caveolin-1 binding and the inherent asymmetry of the membrane leaflets could be the determinant for membrane curvature and needs to be further probed. The review is an important step to reconcile the molecular level details emerging from simulations with the mesoscopic details provided by state of the art experimental approaches.

The Role of Lipids in Allosteric Modulation of Dopamine D2 Receptor-In Silico Study

The dopamine D₂ receptor, belonging to the class A G protein-coupled receptors (GPCRs), is an important drug target for several diseases, including schizophrenia and Parkinson’s disease. The D₂ receptor can be activated by the natural neurotransmitter dopamine or by synthetic ligands, which in both cases leads to the receptor coupling with a G protein. In addition to receptor modulation by orthosteric or allosteric ligands, it has been shown that lipids may affect the behaviour of membrane proteins. We constructed a model of a D₂ receptor with a long intracellular loop (ICL3) coupled with G_iα1 or G_iα2 proteins, embedded in a complex asymmetric membrane, and simulated it in complex with positive, negative or neutral allosteric ligands. In this study, we focused on the influence of ligand binding and G protein coupling on the membrane–receptor interactions. We show that there is a noticeable interplay between the cell membrane, G proteins, D₂ receptor and its modulators.

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.

Obesity epidemic through the prism of evolutionary processes

Currently, obesity has become one of the most serious public health problems. It takes millions of lives worldwide every year due to its association with numerous diseases and leads to significant social and economic losses. It is generally accepted that obesity is the result of the interaction of genes and environment, and the predisposition to it lies in our evolutionary past. This review discusses the role of adipose tissue in human evolution, the factors specifying a person’s predisposition to obesity, the main hypotheses for obesity origin, and potential prevention and treatment strategies arising from them. The evolutionary significance of visceral adipose tissue and some ethnic and sex characteristics associated with its distribution are also considered.

Lipids and Phosphorylation Conjointly Modulate Complex Formation of β2-Adrenergic Receptor and β-arrestin2

G protein-coupled receptors (GPCRs) are the largest class of human membrane proteins that bind extracellular ligands at their orthosteric binding pocket to transmit signals to the cell interior. Ligand binding evokes conformational changes in GPCRs that trigger the binding of intracellular interaction partners (G proteins, G protein kinases, and arrestins), which initiate diverse cellular responses. It has become increasingly evident that the preference of a GPCR for a certain intracellular interaction partner is modulated by a diverse range of factors, e.g., ligands or lipids embedding the transmembrane receptor. Here, by means of molecular dynamics simulations of the β₂-adrenergic receptor and β-arrestin2, we study how membrane lipids and receptor phosphorylation regulate GPCR-arrestin complex conformation and dynamics. We find that phosphorylation drives the receptor’s intracellular loop 3 (ICL3) away from a native negatively charged membrane surface to interact with arrestin. If the receptor is embedded in a neutral membrane, the phosphorylated ICL3 attaches to the membrane surface, which widely opens the receptor core. This opening, which is similar to the opening in the G protein-bound state, weakens the binding of arrestin. The loss of binding specificity is manifested by shallower arrestin insertion into the receptor core and higher dynamics of the receptor-arrestin complex. Our results show that receptor phosphorylation and the local membrane composition cooperatively fine-tune GPCR-mediated signal transduction. Moreover, the results suggest that deeper understanding of complex GPCR regulation mechanisms is necessary to discover novel pathways of pharmacological intervention.

The Role of Plasma Membrane Viscosity in the Response and Resistance of Cancer Cells to Oxaliplatin

Simple Summary

Understanding the role of the plasma membrane in the responses of cancer cells to chemotherapy is important because the cell membrane is directly involved in drug transport and the regulation of numerous biological processes. However, the role of the plasma membrane in cell resistance to platinum drugs like oxaliplatin is not fully understood. In this study we identified the changes to plasma membrane viscosity and lipid composition induced by oxaliplatin in responsive, cultured cancer cells and in mouse tumors. It was also found that the acquisition of chemoresistance is accompanied by modification of membrane lipids in ways that preserve the viscous properties unchanged upon further treatment. Therefore, new therapeutic approaches could be developed to reverse chemoresistance based on membrane lipid modifications and the de-stabilisation of membrane viscosity.

Abstract

Maintenance of the biophysical properties of membranes is essential for cell survival upon external perturbations. However, the links between a fluid membrane state and the drug resistance of cancer cells remain elusive. Here, we investigated the role of membrane viscosity and lipid composition in the responses of cancer cells to oxaliplatin and the development of chemoresistance. Plasma membrane viscosity was monitored in live colorectal cancer cells and tumor xenografts using two-photon excited fluorescence lifetime imaging microscopy (FLIM) using the fluorescent molecular rotor BODIPY 2. The lipid profile was analyzed using time-of-flight secondary ion mass spectrometry (ToF-SIMS). It was found that the plasma membrane viscosity increased upon oxaliplatin treatment, both in vitro and in vivo, and that this correlated with lower phosphatidylcholine and higher cholesterol content. The emergence of resistance to oxaliplatin was accompanied by homeostatic adaptation of the membrane lipidome, and the recovery of lower viscosity. These results suggest that maintaining a constant plasma membrane viscosity via remodeling of the lipid profile is crucial for drug resistance in cancer.

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

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

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.

Evaluating Coarse-Grained MARTINI Force-Fields for Capturing the Ripple Phase of Lipid Membranes

Phospholipids, which are an integral component of cell membranes, exhibit a rich variety of lamellar phases modulated by temperature and composition. Molecular dynamics (MD) simulations have greatly enhanced our understanding of phospholipid membranes by capturing experimentally observed phases and phase transitions at molecular resolution. However, the ripple (P_β') membrane phase, observed as an intermediate phase below the main gel-to-liquid crystalline transition with some lipids, has been challenging to capture with MD simulations, both at all-atom and coarse-grained (CG) resolutions. Here, with an aggregate ∼2.5 μs all-atom and ∼122 μs CGMD simulations, we systematically assess the ability of six CG MARTINI 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid and water force-field (FF) variants, parametrized to capture the DPPC gel and fluid phases, for their ability to capture the P_β' phase, and compared observations with those from an all-atom FF. Upon cooling from the fluid phase to below the phase transition temperature with smaller (380-lipid) and larger (>2200-lipid) MARTINI and all-atom (CHARMM36 FF) DPPC lipid bilayers, we observed that smaller bilayers with both all-atom and MARTINI FFs sampled interdigitated P_β' and ripple-like states, respectively. However, while all-atom simulations of the larger DPPC membranes exhibited the formation of the P_β' phase, MARTINI membranes did not sample interdigitated ripple-like states at larger system sizes. We then demonstrated that the ripple-like states in smaller MARTINI membranes were kinetically trapped structures caused by finite size effects rather than being representative of true P_β' phases. We showed that a MARTINI FF variant that could capture the tilted L_β' gel phase, a prerequisite for stabilizing the P_β' phase, was unable to capture the rippled phase upon cooling. Our study reveals that the current MARTINI FFs (including MARTINI3) may require specific reparametrization of the interaction potentials to stabilize lipid interdigitation, a characteristic of the ripple phase.

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

Polyunsaturated Fatty Acid Modulates Membrane-Bound Monomeric α-Synuclein by Modulating Membrane Microenvironment through Preferential Interactions

There is ample evidence that both native functions and pathogenic aggregation of α-synuclein are intimately dependent on lipid interactions and fatty acid type; the regulatory mechanism however remains unclear. In the present work, using extensive atomistic molecular dynamics simulations and enhanced-sampling, we have focused on exploring the mechanism of fatty acid dependent regulation of monomeric α-Syn₁₀₀ in a native synaptic vesicle-like membrane. Our results show that α-Syn₁₀₀ spontaneously binds to the membrane through its N-terminal region (residues 1-34), where the depth of membrane insertion, the structure, and orientation of the membrane-bound α-Syn₁₀₀ and its impact on membrane structure are modulated by docosahexaenoic acid (DHA). DHA is a polyunsaturated fatty acid abundantly found in the brain and known to promote the oligomerization of α-synuclein. We found that DHA exhibits marked propensity to interact with monomeric α-Syn₁₀₀ and modulates the microenvironment of the protein by preferentially sorting DHA-containing phospholipids, depleting other phospholipids and cholesterol as well as increasing the proportion of anionic to neutral lipids in the immediate vicinity of the protein. Owing to the unique conformational flexibility, DHA chains form more lipid-packing defects in the membrane and efficiently coat the membrane-embedded surface of the protein, compared to the saturated and monounsaturated fatty acids. DHA thus makes the bilayer more amiable to protein adsorption and less prone to α-synuclein-induced perturbation associated with cytotoxicity. Indeed, in the absence of DHA, we observed significant thinning of the local bilayer membrane induced by α-Syn₁₀₀. Though α-Syn₁₀₀ is predominantly α-helical in membranes studied here, in the presence of DHA we observe formation of β-sheet/β-strands in the C-terminal region (residues 35-100) of α-Syn₁₀₀, which is extended out from the membrane surface. Notably, DHA induces β structure in the NAC domain of α-Syn₁₀₀ and promotes extended conformations as well as large solvent exposure of this hydrophobic domain, properties that are known to facilitate self-assembly of α-synuclein. To the best of our knowledge, this study for the first time provides the atomistic insights into DHA-induced regulatory mechanism of monomeric α-synuclein, having implications in protein structure and its physiological/pathological functions.

Ion channels as lipid sensors: from structures to mechanisms

Ion channels play critical roles in cellular function by facilitating the flow of ions across the membrane in response to chemical or mechanical stimuli. Ion channels operate in a lipid bilayer, which can modulate or define their function. Recent technical advancements have led to the solution of numerous ion channel structures solubilized in detergent and/or reconstituted into lipid bilayers, thus providing unprecedented insight into the mechanisms underlying ion channel-lipid interactions. Here, we describe how ion channel structures have evolved to respond to both lipid modulators and lipid activators to control the electrical activities of cells, highlighting diverse mechanisms and common themes.

Lipidome-wide characterization of phosphatidylinositols and phosphatidylglycerols on CC location level

Phosphatidylglycerol (PG) and phosphatidylinositol (PI) are two essential classes of glycerophospholipids (GPs), playing versatile roles such as signalling messengers and lipid-protein interaction ligands in cell. Although a majority of PG and PI molecular species contain unsaturated fatty acyl chain(s), conventional tandem mass spectrometry (MS/MS) methods cannot discern isomers different in carbon-carbon double bond (CC) locations. In this work, we paired phosphate methylation with acetone Paternò-Büchi (PB) reaction, aiming to provide a solution for sensitive and structurally informative analysis of these two important classes of GPs down to the location of CC. A liquid chromatography-tandem mass spectrometry (LC-MS/MS) workflow was established. Offline methylated PG or PI mixtures were subjected to hydrophilic interaction chromatographic separation, online acetone PB reaction, and MS/MS via collision-induced dissociation (CID) for CC location determination in positive ion mode. This method was sensitive, offering limit of identification at 5 nM for both PG and PI standards down to CC locations. On molecular species level, 49 PI and 31 PG were identified from bovine liver, while 61 PIs were identified from human plasma. This workflow also enabled ratiometric comparisons of CC location isomers (C18:1 Δ9 vs. Δ11) of a series of PIs from type 2 diabetes (T2D) plasma to that of normal plasma samples. PI 16:0_18:1 and PI 18:0_18:1 were found to exhibit significant changes in CC isomeric ratios between T2D and normal plasma samples. The above results demonstrate that the developed LC-PB-MS/MS workflow is applicable to different classes of lipids and compatible with other established lipid derivatization methods to achieve comprehensive lipid analysis.

Measuring Remodeling of the Lipid Environment Surrounding Membrane Proteins with Lipid Exchange and Native Mass Spectrometry

Due to their crucial biochemical roles, membrane proteins are important drug targets. Although it is clear that lipids can influence membrane protein function, the chemistry of lipid binding remains difficult to study because protein-lipid interactions are polydisperse, competitive, and transient. Furthermore, detergents, which are often used to solubilize membrane proteins in micelles, may disrupt lipid interactions that occur in bilayers. Here, we present two new approaches to quantify protein-lipid interactions in bilayers and understand how membrane proteins remodel their surrounding lipid environment. First, we used mass spectrometry (MS) to measure the exchange of lipids between lipoprotein nanodiscs with and without an embedded membrane protein. Shifts in the lipid distribution toward the membrane protein nanodiscs revealed lipid binding, and titrations allowed measurement of the optimal lipid composition for the membrane protein. Second, we used native or nondenaturing MS to ionize membrane protein nanodiscs with heterogeneous lipids. Ejecting the membrane protein complex with bound lipids in the mass spectrometer revealed enrichment of specific lipids around the membrane protein. Both new approaches showed that the E. coli ammonium transporter AmtB prefers phosphatidylglycerol lipids overall but has a minor affinity for phosphatidylcholine lipids.

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.

G-Protein-Coupled Receptor-Membrane Interactions Depend on the Receptor Activation State

G-protein-coupled receptors (GPCRs) are the largest family of human membrane proteins and serve as primary targets of approximately one-third of currently marketed drugs. In particular, adenosine A1 receptor (A1 AR) is an important therapeutic target for treating cardiac ischemia-reperfusion injuries, neuropathic pain, and renal diseases. As a prototypical GPCR, the A1 AR is located within a phospholipid membrane bilayer and transmits cellular signals by changing between different conformational states. It is important to elucidate the lipid-protein interactions in order to understand the functional mechanism of GPCRs. Here, all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method were performed on both the inactive (antagonist bound) and active (agonist and G-protein bound) A1 AR, which was embedded in a 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) lipid bilayer. In the GaMD simulations, the membrane lipids played a key role in stabilizing different conformational states of the A1 AR. Our simulations further identified important regions of the receptor that interacted distinctly with the lipids in highly correlated manner. Activation of the A1 AR led to differential dynamics in the upper and lower leaflets of the lipid bilayer. In summary, GaMD enhanced simulations have revealed strongly coupled dynamics of the GPCR and lipids that depend on the receptor activation state. © 2019 Wiley Periodicals, Inc.