Computational investigation of cholesterol binding sites on mitochondrial VDAC

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.

ATP functions as a primary alarmin in allergen-induced type 2 immunity

Environmental allergens that interact with the airway epithelium can activate cellular stress pathways that lead to the release of danger signals known as alarmins. The mechanisms of alarmin release are distinct from damage-associated molecular patterns (DAMPs), which typically escape from cells after loss of plasma membrane integrity. Oxidative stress represents a form of allergen-induced cellular stress that stimulates oxidant-sensing mechanisms coupled to pathways, which facilitate alarmin mobilization and efflux across the plasma membrane. In this review, we highlight examples of alarmin release and discuss their roles in the initiation of type 2 immunity and allergic airway inflammation. In addition, we discuss the concept of alarmin amplification, where "primary" alarmins, which are directly released in response to a specific cellular stress, stimulate additional signaling pathways that lead to secretion of "secondary" alarmins that include proinflammatory cytokines, such as IL-33, as well as genomic and mitochondrial DNA that coordinate or amplify type 2 immunity. Accordingly, allergen-evoked cellular stress can elicit a hierarchy of alarmin signaling responses from the airway epithelium that trigger local innate immune reactions, impact adaptive immunity, and exacerbate diseases including asthma and other chronic inflammatory conditions that affect airway function.

Martini 3 Coarse-Grained Force Field for Cholesterol

Cholesterol plays a crucial role in biomembranes by regulating various properties, such as fluidity, rigidity, permeability, and organization of lipid bilayers. The latest version of the Martini model, Martini 3, offers significant improvements in interaction balance, molecular packing, and inclusion of new bead types and sizes. However, the release of the new model resulted in the need to reparameterize many core molecules, including cholesterol. Here, we describe the development and validation of a Martini 3 cholesterol model, addressing issues related to its bonded setup, shape, volume, and hydrophobicity. The proposed model mitigates some limitations of its Martini 2 predecessor while maintaining or improving the overall behavior.

E as in Enigma: The Mysterious Role of the Voltage-Dependent Anion Channel Glutamate E73

The voltage-dependent anion channel (VDAC) is the main passageway for ions and metabolites over the outer mitochondrial membrane. It was associated with many physiological processes, including apoptosis and modulation of intracellular Ca²⁺ signaling. The protein is formed by a barrel of 19 beta-sheets with an N-terminal helix lining the inner pore. Despite its large diameter, the channel can change its selectivity for ions and metabolites based on its open state to regulate transport into and out of mitochondria. VDAC was shown to be regulated by a variety of cellular factors and molecular partners including proteins, lipids and ions. Although the physiological importance of many of these modulatory effects are well described, the binding sites for molecular partners are still largely unknown. The highly symmetrical and sleek structure of the channel makes predictions of functional moieties difficult. However, one residue repeatedly sticks out when reviewing VDAC literature. A glutamate at position 73 (E73) located on the outside of the channel facing the hydrophobic membrane environment was repeatedly proposed to be involved in channel regulation on multiple levels. Here, we review the distinct hypothesized roles of E73 and summarize the open questions around this mysterious residue.

Mitochondrial dynamics, Leydig cell function, and age-related testosterone deficiency

The mitochondrial translocator protein (18 kDa; TSPO) is a high-affinity cholesterol-binding protein that is an integral component of the cholesterol trafficking scaffold responsible for determining the rate of cholesterol import into the mitochondria for steroid biosynthesis. Previous studies have shown that TSPO declines in aging Leydig cells (LCs) and that its decline is associated with depressed circulating testosterone levels in aging rats. However, TSPO's role in the mechanistic decline in LC function is not fully understood. To address the role of TSPO depletion in LC function, we first examined mitochondrial quality in Tspo knockout mouse tumor MA-10 nG1 LCs compared to wild-type MA-10 cells. Tspo deletion caused a disruption in mitochondrial function and membrane dynamics. Increasing mitochondrial fusion via treatment with the mitochondrial fusion promoter M1 or by optic atrophy 1 (OPA1) overexpression resulted in the restoration of mitochondrial function and mitochondrial morphology as well as in steroid formation in TSPO-depleted nG1 LCs. LCs isolated from aged rats form less testosterone than LCs isolated from young rats. Treatment of aging LCs with M1 improved mitochondrial function and increased androgen formation, suggesting that aging LC dysfunction may stem from compromised mitochondrial dynamics caused by the age-dependent LC TSPO decline. These results, taken together, suggest that maintaining or enhancing mitochondrial fusion may provide therapeutic strategies to maintain or restore testosterone levels with aging.

Regulation of PD-L1 through direct binding of cholesterol to CRAC motifs

Cholesterol, an essential molecule for cell structure, function, and viability, plays crucial roles in the development, progression, and survival of cancer cells. Earlier studies have shown that cholesterol-lowering drugs can inhibit the high expression of programmed-death ligand 1 (PD-L1) that contributes to immunoevasion in cancer cells. However, the regulatory mechanism of cell surface PD-L1 abundance by cholesterol is still controversial. Here, using nuclear magnetic resonance and biochemical techniques, we demonstrated that cholesterol can directly bind to the transmembrane domain of PD-L1 through two cholesterol-recognition amino acid consensus (CRAC) motifs, forming a sandwich-like architecture and stabilizing PD-L1 to prevent downstream degradation. Mutations at key binding residues prohibit PD-L1-cholesterol interactions, decreasing the cellular abundance of PD-L1. Our results reveal a unique regulatory mechanism that controls the stability of PD-L1 in cancer cells, providing an alternative method to overcome PD-L1-mediated immunoevasion in cancers.

Structure and Gating Behavior of the Human Integral Membrane Protein VDAC1 in a Lipid Bilayer

The voltage-dependent anion channel (VDAC), the most abundant protein in the outer mitochondrial membrane, is responsible for the transport of all ions and metabolites into and out of mitochondria. Larger than any of the β-barrel structures determined to date by magic-angle spinning (MAS) NMR, but smaller than the size limit of cryo-electron microscopy (cryo-EM), VDAC1’s 31 kDa size has long been a bottleneck in determining its structure in a near-native lipid bilayer environment. Using a single two-dimensional (2D) crystalline sample of human VDAC1 in lipids, we applied proton-detected fast magic-angle spinning NMR spectroscopy to determine the arrangement of β strands. Combining these data with long-range restraints from a spin-labeled sample, chemical shift-based secondary structure prediction, and previous MAS NMR and atomic force microscopy (AFM) data, we determined the channel’s structure at a 2.2 Å root-mean-square deviation (RMSD). The structure, a 19-stranded β-barrel, with an N-terminal α-helix in the pore is in agreement with previous data in detergent, which was questioned due to the potential for the detergent to perturb the protein’s functional structure. Using a quintuple mutant implementing the channel’s closed state, we found that dynamics are a key element in the protein’s gating behavior, as channel closure leads to the destabilization of not only the C-terminal barrel residues but also the α2 helix. We showed that cholesterol, previously shown to reduce the frequency of channel closure, stabilizes the barrel relative to the N-terminal helix. Furthermore, we observed channel closure through steric blockage by a drug shown to selectively bind to the channel, the Bcl2-antisense oligonucleotide G3139.

Lipids in Pathophysiology and Development of the Membrane Lipid Therapy: New Bioactive Lipids

Membranes are mainly composed of a lipid bilayer and proteins, constituting a checkpoint for the entry and passage of signals and other molecules. Their composition can be modulated by diet, pathophysiological processes, and nutritional/pharmaceutical interventions. In addition to their use as an energy source, lipids have important structural and functional roles, e.g., fatty acyl moieties in phospholipids have distinct impacts on human health depending on their saturation, carbon length, and isometry. These and other membrane lipids have quite specific effects on the lipid bilayer structure, which regulates the interaction with signaling proteins. Alterations to lipids have been associated with important diseases, and, consequently, normalization of these alterations or regulatory interventions that control membrane lipid composition have therapeutic potential. This approach, termed membrane lipid therapy or membrane lipid replacement, has emerged as a novel technology platform for nutraceutical interventions and drug discovery. Several clinical trials and therapeutic products have validated this technology based on the understanding of membrane structure and function. The present review analyzes the molecular basis of this innovative approach, describing how membrane lipid composition and structure affects protein-lipid interactions, cell signaling, disease, and therapy (e.g., fatigue and cardiovascular, neurodegenerative, tumor, infectious diseases).

Lysosome activation in peripheral blood mononuclear cells and prognostic significance of circulating LC3B in COVID-19

Coronavirus disease 2019 (COVID-19) has spread rapidly worldwide, causing significant mortality. There is a mechanistic relationship between intracellular coronavirus replication and deregulated autophagosome–lysosome system. We performed transcriptome analysis of peripheral blood mononuclear cells (PBMCs) from COVID-19 patients and identified the aberrant upregulation of genes in the lysosome pathway. We further determined the capability of two circulating markers, namely microtubule-associated proteins 1A/1B light chain 3B (LC3B) and (p62/SQSTM1) p62, both of which depend on lysosome for degradation, in predicting the emergence of moderate-to-severe disease in COVID-19 patients requiring hospitalization for supplemental oxygen therapy. Logistic regression analyses showed that LC3B was associated with moderate-to-severe COVID-19, independent of age, sex and clinical risk score. A decrease in LC3B concentration <5.5 ng/ml increased the risk of oxygen and ventilatory requirement (adjusted odds ratio: 4.6; 95% CI: 1.1–22.0; P = 0.04). Serum concentrations of p62 in the moderate-to-severe group were significantly lower in patients aged 50 or below. In conclusion, lysosome function is deregulated in PBMCs isolated from COVID-19 patients, and the related biomarker LC3B may serve as a novel tool for stratifying patients with moderate-to-severe COVID-19 from those with asymptomatic or mild disease. COVID-19 patients with a decrease in LC3B concentration <5.5 ng/ml will require early hospital admission for supplemental oxygen therapy and other respiratory support.

The Open State Selectivity of the Bean Seed VDAC Depends on Stigmasterol and Ion Concentration

The voltage-dependent anion channel (VDAC) is the major pathway for metabolites and ions transport through the mitochondrial outer membrane. It can regulate the flow of solutes by switching to a low conductance state correlated with a selectivity reversal, or by a selectivity inversion of its open state. The later one was observed in non-plant VDACs and is poorly characterized. We aim at investigating the selectivity inversion of the open state using plant VDAC purified from Phaseolus coccineus (PcVDAC) to evaluate its physiological role. Our main findings are: (1) The VDAC selectivity inversion of the open state occurs in PcVDAC, (2) Ion concentration and stigmasterol affect the occurrence of the open state selectivity inversion and stigmasterol appears to interact directly with PcVDAC. Interestingly, electrophysiological data concerning the selectivity inversion of the PcVDAC open state suggests that the phenomenon probably does not have a significant physiological effect in vivo.

Warburg's Ghost-Cancer's Self-Sustaining Phenotype: The Aberrant Carbon Flux in Cholesterol-Enriched Tumor Mitochondria via Deregulated Cholesterogenesis

Interpreting connections between the multiple networks of cell metabolism is indispensable for understanding how cells maintain homeostasis or transform into the decontrolled proliferation phenotype of cancer. Situated at a critical metabolic intersection, citrate, derived via glycolysis, serves as either a combustible fuel for aerobic mitochondrial bioenergetics or as a continuously replenished cytosolic carbon source for lipid biosynthesis, an essentially anaerobic process. Therein lies the paradox: under what conditions do cells control the metabolic route by which they process citrate? The Warburg effect exposes essentially the same dilemma—why do cancer cells, despite an abundance of oxygen needed for energy-generating mitochondrial respiration with citrate as fuel, avoid catabolizing mitochondrial citrate and instead rely upon accelerated glycolysis to support their energy requirements? This review details the genesis and consequences of the metabolic paradigm of a “truncated” Krebs/TCA cycle. Abundant data are presented for substrate utilization and membrane cholesterol enrichment in tumors that are consistent with criteria of the Warburg effect. From healthy cellular homeostasis to the uncontrolled proliferation of tumors, metabolic alterations center upon the loss of regulation of the cholesterol biosynthetic pathway. Deregulated tumor cholesterogenesis at the HMGR locus, generating enhanced carbon flux through the cholesterol synthesis pathway, is an absolute prerequisite for DNA synthesis and cell division. Therefore, expedited citrate efflux from cholesterol-enriched tumor mitochondria via the CTP/SLC25A1 citrate transporter is fundamental for sustaining the constant demand for cytosolic citrate that fuels the elevated flow of carbons from acetyl-CoA through the deregulated pathway of cholesterol biosynthesis.

VDAC regulation of mitochondrial calcium flux: From channel biophysics to disease

Voltage-dependent anion channel (VDAC), the most abundant mitochondrial outer membrane protein, is important for a variety of mitochondrial functions including metabolite exchange, calcium transport, and apoptosis. While VDAC's role in shuttling metabolites between the cytosol and mitochondria is well established, there is a growing interest in understanding the mechanisms of its regulation of mitochondrial calcium transport. Here we review the current literature on VDAC's role in calcium signaling, its biophysical properties, physiological function, and pathology focusing on its importance in cardiac diseases. We discuss the specific biophysical properties of the three VDAC isoforms in mammalian cells-VDAC 1, 2, and 3-in relationship to calcium transport and their distinct roles in cell physiology and disease. Highlighting the emerging evidence that cytosolic proteins interact with VDAC and regulate its calcium permeability, we advocate for continued investigation into the VDAC interactome at the contact sites between mitochondria and organelles and its role in mitochondrial calcium transport.

Structure, gating and interactions of the voltage-dependent anion channel

The voltage-dependent anion channel (VDAC) is one of the most highly abundant proteins found in the outer mitochondrial membrane, and was one of the earliest discovered. Here we review progress in understanding VDAC function with a focus on its structure, discussing various models proposed for voltage gating as well as potential drug targets to modulate the channel's function. In addition, we explore the sensitivity of VDAC structure to variations in the membrane environment, comparing DMPC-only, DMPC with cholesterol, and near-native lipid compositions, and use magic-angle spinning NMR spectroscopy to locate cholesterol on the outside of the β-barrel. We find that the VDAC protein structure remains unchanged in different membrane compositions, including conditions with cholesterol.

Understanding the complex pharmacology of cannabidiol: Mounting evidence suggests a common binding site with cholesterol

Cannabidiol is claimed to bind to a large number of protein targets based on in vitro assays. This suggests opportunities for a wide range of therapeutic applications. On the other hand, the existence of phytochemical 'nuisance compounds' suggests some measure of caution - these compounds are capable of altering membrane biophysical properties and changing protein function without directly contacting a binding site. Like cannabidiol, cholesterol alters membrane properties, but it also binds directly to membrane proteins through abundant cholesterol recognition sites. We present the evidence that cannabidiol and cholesterol may bind to the same site on some proteins. As a starting point for further research, we also used blind docking to show that cannabidiol binds to a cholesterol binding site on the CB1 receptor. Elucidation of the mechanism(s) of action of cannabidiol will assist the prioritisation of in vitro hits across targets, improve the success rate of medicinal chemistry campaigns, and ultimately benefit patient populations by focusing resources on programs with the most translational potential.

Lipid composition of the cancer cell membrane

Cancer cell possesses numerous adaptations to resist the immune system response and chemotherapy. One of the most significant properties of the neoplastic cells is the altered lipid metabolism, and consequently, the abnormal cell membrane composition. Like in the case of phosphatidylcholine, these changes result in the modulation of certain enzymes and accumulation of energetic material, which could be used for a higher proliferation rate. The changes are so prominent, that some lipids, such as phosphatidylserines, could even be considered as the cancer biomarkers. Additionally, some changes of biophysical properties of cell membranes lead to the higher resistance to chemotherapy, and finally to the disturbances in signalling pathways. Namely, the increased levels of certain lipids, like for instance phosphatidylserine, lead to the attenuation of the immune system response. Also, changes in lipid saturation prevent the cells from demanding conditions of the microenvironment. Particularly interesting is the significance of cell membrane cholesterol content in the modulation of metastasis. This review paper discusses the roles of each lipid type in cancer physiology. The review combined theoretical data with clinical studies to show novel therapeutic options concerning the modulation of cell membranes in oncology.

Olesoxime, a cholesterol-like neuroprotectant restrains synaptic vesicle exocytosis in the mice motor nerve terminals: Possible role of VDACs

Olesoxime is a cholesterol-like neuroprotective compound that targets to mitochondrial voltage dependent anion channels (VDACs). VDACs were also found in the plasma membrane and highly expressed in the presynaptic compartment. Here, we studied the effects of olesoxime and VDAC inhibitors on neurotransmission in the mouse neuromuscular junction. Electrophysiological analysis revealed that olesoxime suppressed selectively evoked neurotransmitter release in response to a single stimulus and 20 Hz activity. Also olesoxime decreased the rate of FM1-43 dye loss (an indicator of synaptic vesicle exocytosis) at low frequency stimulation and 20 Hz. Furthermore, an increase in extracellular Cl^- enhanced the action of olesoxime on the exocytosis and olesoxime increased intracellular Cl^- levels. The effects of olesoxime on the evoked synaptic vesicle exocytosis and [Cl^-]_i were blocked by membrane-permeable and impermeable VDAC inhibitors. Immunofluorescent labeling pointed on the presence of VDACs on the synaptic membranes. Rotenone-induced mitochondrial dysfunction perturbed the exocytotic release of FM1-43 and cell-permeable VDAC inhibitor (but not olesoxime or impermeable VDAC inhibitor) partially mitigated the rotenone-driven alterations in the FM1-43 unloading and mitochondrial superoxide production. Thus, olesoxime restrains neurotransmission by acting on plasmalemmal VDACs whose activation can limit synaptic vesicle exocytosis probably via increasing anion flux into the nerve terminals.

Targeting the Multiple Physiologic Roles of VDAC With Steroids and Hydrophobic Drugs

There is accumulating evidence that endogenous steroids and non-polar drugs are involved in the regulation of mitochondrial physiology. Many of these hydrophobic compounds interact with the Voltage Dependent Anion Channel (VDAC). This major metabolite channel in the mitochondrial outer membrane (MOM) regulates the exchange of ions and water-soluble metabolites, such as ATP and ADP, across the MOM, thus governing mitochondrial respiration. Proteomics and biochemical approaches together with molecular dynamics simulations have identified an impressively large number of non-polar compounds, including endogenous, able to bind to VDAC. These findings have sparked speculation that both natural steroids and synthetic hydrophobic drugs regulate mitochondrial physiology by directly affecting VDAC ion channel properties and modulating its metabolite permeability. Here we evaluate recent studies investigating the effect of identified VDAC-binding natural steroids and non-polar drugs on VDAC channel functioning. We argue that while many compounds are found to bind to the VDAC protein, they do not necessarily affect its channel functions in vitro. However, they may modify other aspects of VDAC physiology such as interaction with its cytosolic partner proteins or complex formation with other mitochondrial membrane proteins, thus altering mitochondrial function.

Fungal allergen-induced IL-33 secretion involves cholesterol-dependent, VDAC-1-mediated ATP release from the airway epithelium

KEY POINTS

Alternaria aeroallergens induce the release of ATP from human bronchial epithelial (HBE) cells by activating a conductive pathway involving voltage-dependent anion channel-1 (VDAC-1) and by exocytosis of ATP localized within membrane vesicles. Inhibition of VDAC-1 blocked Alternaria-evoked Ca²⁺ uptake across the plasma membrane of HBE cells and interleukin (IL)-33 release into the extracellular media. Reducing cholesterol content with a cholesterol scavenger (β-methylcyclodextrin) or statin compound (simvastatin) blocked ATP and IL-33 release by lowering the expression of VDAC-1 in the plasma membrane. Pretreatment with simvastatin for 24 h also inhibited the increase in tight junction macromolecule permeability that occurs following Alternaria exposure. These results establish a novel role for VDAC-1 as a mechanism underlying ATP release induced by fungal allergens and suggests a possible therapeutic use for cholesterol lowering compounds in reducing Alternaria-stimulated allergic inflammation.

ABSTRACT

Human bronchial epithelial (HBE) cells exposed to allergens derived from the common saprophytic fungus, Alternaria alternata release ATP, which in turn stimulates P2X₇ receptor-mediated Ca²⁺ uptake across the plasma membrane. The subsequent increase in intracellular calcium concentration induces proteolytic processing and secretion of interleukin (IL)-33, a critical cytokine involved in the initiation of allergic airway inflammation. A major objective of the present study was to identify the mechanism responsible for conductive ATP release. The results show that pretreatment of HBE cells with inhibitors of the voltage-dependent anion channel-1 (VDAC-1) or treatment with a VDAC-1 selective blocking antibody or silencing mRNA expression of the channel by RNA interference, inhibit Alternaria-evoked ATP release. Moreover, inhibition of VDAC-1 channel activity or reducing protein expression blocked the secretion of IL-33. Similarly, reducing the cholesterol content of HBE cells with simvastatin or the cholesterol scavenger β-methylcyclodextrin also blocked ATP release and IL-33 secretion by decreasing the level of VDAC-1 expression in the plasma membrane. In addition, simvastatin inhibited the increase in tight junction macromolecule permeability that was previously observed after Alternaria exposure. These results demonstrate a novel function for VDAC-1 as the conductive mechanism responsible for Alternaria-induced ATP release, an essential early step in the processing, mobilization and secretion of IL-33 by the airway epithelium. Furthermore, the simvastatin-evoked reduction of VDAC-1 expression in the plasma membrane, suggests the possibility that cholesterol lowering compounds may be beneficial in alleviating allergic airway inflammation induced by fungal allergens.

Modulation of Transient Receptor Potential C Channel Activity by Cholesterol

Changes of cholesterol level in the plasma membrane of cells have been shown to modulate ion channel function. The proposed mechanisms underlying these modulations include association of cholesterol to a single binding site at a single channel conformation, association to a highly flexible cholesterol binding site adopting multiple poses, and perturbation of lipid rafts. These perturbations have been shown to induce reversible targeting of mammalian transient receptor potential C (TRPC) channels to the cholesterol-rich membrane environment of lipid rafts. Thus, the observed inhibition of TRPC channels by methyl-β-cyclodextrin (MβCD), which induces cholesterol efflux from the plasma membrane, may result from disruption of lipid rafts. This perturbation was also shown to disrupt multimolecular signaling complexes containing TRPC channels. The Drosophila TRP and TRP-like (TRPL) channels belong to the TRPC channel subfamily. When the Drosophila TRPL channel was expressed in S2 or HEK293 cells and perfused with MβCD, the TRPL current was abolished in less than 100 s, fitting well the fast kinetic phase of cholesterol sequestration experiments in cells. It was thus suggested that the fast kinetics of TRPL channel suppression by MβCD arise from disruption of lipid rafts. Accordingly, lipid raft perturbation by cholesterol sequestration could give clues to the function of lipid environment in TRPC channel activity and its mechanism.

Multiple neurosteroid and cholesterol binding sites in voltage-dependent anion channel-1 determined by photo-affinity labeling

Voltage-dependent anion channel-1 (VDAC1) is a mitochondrial porin that is implicated in cellular metabolism and apoptosis, and modulated by numerous small molecules including lipids. VDAC1 binds sterols, including cholesterol and neurosteroids such as allopregnanolone. Biochemical and computational studies suggest that VDAC1 binds multiple cholesterol molecules, but photolabeling studies have identified only a single cholesterol and neurosteroid binding site at E73. To identify all the binding sites of neurosteroids in VDAC1, we apply photo-affinity labeling using two sterol-based photolabeling reagents with complementary photochemistry: 5α-6-AziP which contains an aliphatic diazirine, and KK200 which contains a trifluoromethyl-phenyldiazirine (TPD) group. 5α-6-AziP and KK200 photolabel multiple residues within an E73 pocket confirming the presence of this site and mapping sterol orientation within this pocket. In addition, KK200 photolabels four other sites consistent with the finding that VDAC1 co-purifies with five cholesterol molecules. Both allopregnanolone and cholesterol competitively prevent photolabeling at E73 and three other sites indicating that these are common sterol binding sites shared by both neurosteroids and cholesterol. Binding at the functionally important residue E73 suggests a possible role for sterols in regulating VDAC1 signaling and interaction with partner proteins.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

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

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

Ceramides bind VDAC2 to trigger mitochondrial apoptosis

Ceramides draw wide attention as tumor suppressor lipids that act directly on mitochondria to trigger apoptotic cell death. However, molecular details of the underlying mechanism are largely unknown. Using a photoactivatable ceramide probe, we here identify the voltage-dependent anion channels VDAC1 and VDAC2 as mitochondrial ceramide binding proteins. Coarse-grain molecular dynamics simulations reveal that both channels harbor a ceramide binding site on one side of the barrel wall. This site includes a membrane-buried glutamate that mediates direct contact with the ceramide head group. Substitution or chemical modification of this residue abolishes photolabeling of both channels with the ceramide probe. Unlike VDAC1 removal, loss of VDAC2 or replacing its membrane-facing glutamate with glutamine renders human colon cancer cells largely resistant to ceramide-induced apoptosis. Collectively, our data support a role of VDAC2 as direct effector of ceramide-mediated cell death, providing a molecular framework for how ceramides exert their anti-neoplastic activity.

Recent Applications of Diazirines in Chemical Proteomics

The elucidation of substrate-protein interactions is an important component of the drug development process. Due to the complexity of native cellular environments, elucidating these fundamental biochemical interactions remains challenging. Photoaffinity labeling (PAL) is a versatile technique that can provide insight into ligand-target interactions. By judicious modification of substrates with a photoreactive group, PAL creates a covalent crosslink between a substrate and its biological target following UV-irradiation. Among the commonly employed photoreactive groups, diazirines have emerged as the gold standard. In this Minireview, recent developments in the field of diazirine-based photoaffinity labeling will be discussed, with emphasis being placed on their applications in chemical proteomic studies.

Assessing the role of residue E73 and lipid headgroup charge in VDAC1 voltage gating

The voltage-dependent anion channel (VDAC) is the most abundant protein of the mitochondrial outer membrane (MOM) where it regulates transport of ions and metabolites in and out of the organelle. VDAC function is extensively studied in a lipid bilayer system that allows conductance monitoring of reconstituted channels under applied voltage. The process of switching from a high-conductance state, open to metabolites, to a variety of low-conducting states, which excludes metabolite transport, is termed voltage gating and the mechanism remains poorly understood. Recent studies have implicated the involvement of the membrane-solvated residue E73 in the gating process through β-barrel destabilization. However, there has been no direct experimental evidence of E73 involvement in VDAC1 voltage gating. Here, using electrophysiology measurements, we exclude the involvement of E73 in murine VDAC1 (mVDAC1) voltage gating process. With an established protocol of assessing voltage gating of VDACs reconstituted into planar lipid membranes, we definitively show that mVDAC1 gating properties do not change when E73 is replaced by either a glutamine or an alanine. We further demonstrate that cholesterol has no effect on mVDAC1 gating characteristics, though it was shown that E73 is coordinating residue in the cholesterol binding site. In contrast, we found a pronounced gating effect based on the charge of the phospholipid headgroup, where the positive charge stimulates and negative charge suppresses gating. These findings call for critical evaluation of the existing models of VDAC gating and contribute to our understanding of VDAC's role in control of MOM permeability and regulation of mitochondrial respiration and metabolism.

ERK1/2 and the Bcl-2 Family Proteins Mcl-1, tBid, and Bim Are Involved in Inhibition of Apoptosis During Persistent Chlamydia psittaci Infection

Chlamydia psittaci is an obligate intracellular pathogen that can cause zoonosis. Persistent C. psittaci infection can inhibit apoptosis in host cells, thus extending their survival and enabling them to complete their growth cycle. In this study, the antiapoptotic effects of persistent C. psittaci infection, induced by treatment with IFN-γ, were found to be associated with both the death receptor and the mitochondrial pathways of apoptosis. These effects were mediated by Bcl-2 family members, as evidenced by the decreased expression of proapoptotic proteins, such as tBid and Bim. Simultaneously, the antiapoptotic protein Mcl-1 was upregulated by persistent C. psittaci infection. Increased phosphorylation of ERK1/2 was observed; however, the expression of Bad, unlike that of other proapoptotic proteins, did not seem to be involved in this process. In summary, persistent chlamydial infection exerts antiapoptotic effects through both the death receptor and the mitochondrial pathways, in a process that is regulated by the ERK1/2 and apoptotic proteins of the Bcl-2 family.

Lipid composition and salt concentration as regulatory factors of the anion selectivity of VDAC studied by coarse-grained molecular dynamics simulations

The voltage-dependent anion channel (VDAC) is a mitochondrial outer membrane protein whose fundamental function is to facilitate and regulate the flow of metabolites between the cytosol and the mitochondrial intermembrane space. Using coarse-grained molecular dynamics simulations, we investigated the dependence of VDAC selectivity towards small inorganic anions on two factors: the ionic strength and the lipid composition. In agreement with experimental data we found that VDAC becomes less anion selective with increasing salt concentration due to the screening of a few basic residues that point into the pore lumen. The molecular dynamics simulations provide insight into the regulation mechanism of VDAC selectivity by the composition in the lipid membrane and suggest that the ion distribution is differently modulated by POPE compared to the POPC bilayer. This occurs through the more persistent interactions of acidic residues located at both rims of the β-barrel with head groups of POPE which in turn impact the electrostatic potential and thereby the selectivity of the pore. This mechanism occurs not only in POPE single component membranes but also in a mixed POPE/POPC bilayer by an enrichment of POPE over POPC lipids on the surface of VDAC. Thus we show here that computationally-inexpensive coarse-grained simulations are able to capture, in a semi-quantitative way, essential features of VDAC anion selectivity and could pave the way toward a molecular level understanding of metabolite transport in natural membranes.

Hydrophobic Mismatch Modulates Stability and Plasticity of Human Mitochondrial VDAC2

The human mitochondrial outer membrane protein voltage-dependent anion channel isoform 2 (hVDAC2) is a β-barrel metabolite flux channel that is indispensable for cell survival. It is well established that physical forces imposed on a transmembrane protein by its surrounding lipid environment decide protein structure and stability. Yet, how the mitochondrial membrane and protein-lipid interplay together regulate hVDAC2 stability is unknown. Here, we combine experimental biophysical investigations of protein stability with all-atom molecular dynamics simulations to study the effect of the most abundant mitochondrial phosphocholine (PC) lipids on hVDAC2. We demonstrate experimentally that increasing the PC lipid acyl chain length from diC14:0 to diC18:0-PC has a nonlinear effect on the β-barrel. We show that protein stability is highest in diC16:0-PC, which exhibits a negative mismatch with the hVDAC2 barrel. Our simulations also reveal that structural rigidity of hVDAC2 is highest under optimal negative mismatch provided by diC16:0-PC bilayers. Further, we validate our observations by altering the physical properties of PC membranes indirectly using cholesterol. We propose that VDAC plasticity and stability in the mitochondrial outer membrane are modulated by physical properties of the bilayer.

Molecular Dynamics Simulations of Kir2.2 Interactions with an Ensemble of Cholesterol Molecules

Cholesterol is a major regulator of multiple types of ion channels, but the specific mechanisms and the dynamics of its interactions with the channels are not well understood. Kir2 channels were shown to be sensitive to cholesterol through direct interactions with "cholesterol-sensitive" regions on the channel protein. In this work, we used Martini coarse-grained simulations to analyze the long (μs) timescale dynamics of cholesterol with Kir2.2 channels embedded into a model membrane containing POPC phospholipid with 30 mol% cholesterol. This approach allows us to simulate the dynamic, unbiased migration of cholesterol molecules from the lipid membrane environment to the protein surface of Kir2.2 and explore the favorability of cholesterol interactions at both surface sites and recessed pockets of the channel. We found that the cholesterol environment surrounding Kir channels forms a complex milieu of different short- and long-term interactions, with multiple cholesterol molecules concurrently interacting with the channel. Furthermore, utilizing principles from network theory, we identified four discrete cholesterol-binding sites within the previously identified cholesterol-sensitive region that exist depending on the conformational state of the channel-open or closed. We also discovered that a twofold decrease in the cholesterol level of the membrane, which we found earlier to increase Kir2 activity, results in a site-specific decrease of cholesterol occupancy at these sites in both the open and closed states: cholesterol molecules at the deepest of these discrete sites shows no change in occupancy at different cholesterol levels, whereas the remaining sites showed a marked decrease in occupancy.

Challenges and approaches to understand cholesterol-binding impact on membrane protein function: an NMR view

Experimental evidence for a direct role of lipids in determining the structure, dynamics, and function of membrane proteins leads to the term 'functional lipids'. In particular, the sterol molecule cholesterol modulates the activity of many membrane proteins. The precise nature of cholesterol-binding sites and the consequences of modulation of local membrane micro-viscosity by cholesterol, however, is often unknown. Here, we review the current knowledge of the interaction of cholesterol with transmembrane proteins, with a special focus on structural aspects of the interaction derived from nuclear magnetic resonance approaches. We highlight examples of the importance of cholesterol modulation of membrane protein function, discuss the specificity of cholesterol binding, and review the proposed binding motifs from a molecular perspective. We conclude with a short perspective on what could be future trends in research efforts targeted towards a better understanding of cholesterol/membrane protein interactions.

Lipid-Mediated Modulation of Intracellular Ion Channels and Redox State: Physiopathological Implications

Significance: Ion channels play an important role in the regulation of organelle function within the cell, as proven by increasing evidence pointing to a link between altered function of intracellular ion channels and different pathologies ranging from cancer to neurodegenerative diseases, ischemic damage, and lysosomal storage diseases. Recent Advances: A link between these pathologies and redox state as well as lipid homeostasis and ion channel function is in the focus of current research. Critical Issues: Ion channels are target of modulation by lipids and lipid messengers, although in most cases the mechanistic details have not been clarified yet. Ion channel function importantly impacts production of reactive oxygen species (ROS), especially in the case of mitochondria and lysosomes. ROS, in turn, may modulate the function of intracellular channels triggering thereby a feedback control under physiological conditions. If produced in excess, ROS can be harmful to lipids and may produce oxidized forms of these membrane constituents that ultimately affect ion channel function by triggering a "circulus vitiosus." Future Directions: The present review summarizes our current knowledge about the contribution of intracellular channels to oxidative stress and gives examples of how these channels are modulated by lipids and how this modulation may affect ROS production in ROS-related diseases. Future studies need to address the importance of the regulation of intracellular ion channels and related oxidative stress by lipids in various physiological and pathological contexts. Antioxid. Redox Signal. 28, 949-972.

Decrypting Strong and Weak Single-Walled Carbon Nanotubes Interactions with Mitochondrial Voltage-Dependent Anion Channels Using Molecular Docking and Perturbation Theory

The current molecular docking study provided the Free Energy of Binding (FEB) for the interaction (nanotoxicity) between VDAC mitochondrial channels of three species (VDAC1-Mus musculus, VDAC1-Homo sapiens, VDAC2-Danio rerio) with SWCNT-H, SWCNT-OH, SWCNT-COOH carbon nanotubes. The general results showed that the FEB values were statistically more negative (p < 0.05) in the following order: (SWCNT-VDAC2-Danio rerio) > (SWCNT-VDAC1-Mus musculus) > (SWCNT-VDAC1-Homo sapiens) > (ATP-VDAC). More negative FEB values for SWCNT-COOH and OH were found in VDAC2-Danio rerio when compared with VDAC1-Mus musculus and VDAC1-Homo sapiens (p < 0.05). In addition, a significant correlation (0.66 > r² > 0.97) was observed between n-Hamada index and VDAC nanotoxicity (or FEB) for the zigzag topologies of SWCNT-COOH and SWCNT-OH. Predictive Nanoparticles-Quantitative-Structure Binding-Relationship models (nano-QSBR) for strong and weak SWCNT-VDAC docking interactions were performed using Perturbation Theory, regression and classification models. Thus, 405 SWCNT-VDAC interactions were predicted using a nano-PT-QSBR classifications model with high accuracy, specificity, and sensitivity (73–98%) in training and validation series, and a maximum AUROC value of 0.978. In addition, the best regression model was obtained with Random Forest (R² of 0.833, RMSE of 0.0844), suggesting an excellent potential to predict SWCNT-VDAC channel nanotoxicity. All study data are available at https://doi.org/10.6084/m9.figshare.4802320.v2.

The mitochondrial VDAC of bean seeds recruits phosphatidylethanolamine lipids for its proper functioning

The voltage-dependent anion-selective channel (VDAC) is the main pathway for inorganic ions and metabolites through the mitochondrial outer membrane. Studies recently demonstrated that membrane lipids regulate its function. It remains, however, unclear how this regulation takes place. In this study, we show that phospholipids are key regulators of Phaseolus VDAC function and, furthermore, that the salt concentration modulates this regulation. Both selectivity and voltage dependence of Phaseolus VDAC are very sensitive to a change in the lipid polar head from PC to PE. Interestingly enough, this dependence is observed only at low salt concentration. Furthermore, significant changes in VDAC functional properties also occur with the gradual methylation of the PE group pointing to the role of subtle chemical variations in the lipid head group. The dependence of PcVDAC gating upon the introduction of a small mole fraction of PE in a PC bilayer has prompted us to propose the existence of a specific interaction site for PE on the outer surface of PcVDAC. Eventually, comparative modeling and molecular dynamics simulations suggest a potential mechanism to get insight into the anion selectivity enhancement of PcVDAC observed in PE relative to PC.

Photoaffinity labeling with cholesterol analogues precisely maps a cholesterol-binding site in voltage-dependent anion channel-1

Voltage-dependent anion channel-1 (VDAC1) is a highly regulated β-barrel membrane protein that mediates transport of ions and metabolites between the mitochondria and cytosol of the cell. VDAC1 co-purifies with cholesterol and is functionally regulated by cholesterol, among other endogenous lipids. Molecular modeling studies based on NMR observations have suggested five cholesterol-binding sites in VDAC1, but direct experimental evidence for these sites is lacking. Here, to determine the sites of cholesterol binding, we photolabeled purified mouse VDAC1 (mVDAC1) with photoactivatable cholesterol analogues and analyzed the photolabeled sites with both top-down mass spectrometry (MS), and bottom-up MS paired with a clickable, stable isotope-labeled tag, FLI-tag. Using cholesterol analogues with a diazirine in either the 7 position of the steroid ring (LKM38) or the aliphatic tail (KK174), we mapped a binding pocket in mVDAC1 localized to Thr⁸³ and Glu⁷³, respectively. When Glu⁷³ was mutated to a glutamine, KK174 no longer photolabeled this residue, but instead labeled the nearby Tyr⁶² within this same binding pocket. The combination of analytical strategies employed in this work permits detailed molecular mapping of a cholesterol-binding site in a protein, including an orientation of the sterol within the site. Our work raises the interesting possibility that cholesterol-mediated regulation of VDAC1 may be facilitated through a specific binding site at the functionally important Glu⁷³ residue.

Translocator Protein 18 kDa (TSPO): An Old Protein with New Functions?

Translocator protein 18 kDa (TSPO) was previously known as the peripheral benzodiazepine receptor (PBR) in eukaryotes, where it is mainly localized to the mitochondrial outer membrane. Considerable evidence indicates that it plays regulatory roles in steroidogenesis and apoptosis and is involved in various human diseases, such as metastatic cancer, Alzheimer's and Parkinson's disease, inflammation, and anxiety disorders. Ligands of TSPO are widely used as diagnostic tools and treatment options, despite there being no clear understanding of the function of TSPO. An ortholog in the photosynthetic bacterium Rhodobacter was independently discovered as the tryptophan-rich sensory protein (TspO) and found to play a role in the response to changes in oxygen and light conditions that regulate photosynthesis and respiration. As part of this highly conserved protein family found in all three kingdoms, the rat TSPO is able to rescue the knockout phenotype in Rhodobacter, indicating functional as well as structural conservation. Recently, a major breakthrough in the field was achieved: the determination of atomic-resolution structures of TSPO from different species by several independent groups. This now allows us to reexamine the function of TSPO with a molecular perspective. In this review, we focus on recently determined structures of TSPO and their implications for potential functions of this ubiquitous multifaceted protein. We suggest that TSPO is an ancient bacterial receptor/stress sensor that has developed additional interactions, partners, and roles in its mitochondrial outer membrane environment in eukaryotes.

Lipid interaction sites on channels, transporters and receptors: Recent insights from molecular dynamics simulations

Lipid molecules are able to selectively interact with specific sites on integral membrane proteins, and modulate their structure and function. Identification and characterization of these sites are of importance for our understanding of the molecular basis of membrane protein function and stability, and may facilitate the design of lipid-like drug molecules. Molecular dynamics simulations provide a powerful tool for the identification of these sites, complementing advances in membrane protein structural biology and biophysics. We describe recent notable biomolecular simulation studies which have identified lipid interaction sites on a range of different membrane proteins. The sites identified in these simulation studies agree well with those identified by complementary experimental techniques. This demonstrates the power of the molecular dynamics approach in the prediction and characterization of lipid interaction sites on integral membrane proteins. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Combined effect of G3139 and TSPO ligands on Ca(2+)-induced permeability transition in rat brain mitochondria

Permeability of the mitochondrial outer membrane is determined by the activity of voltage-dependent anion channels (VDAC) which are regulated by many factors and proteins. One of the main partner-regulator of VDAC is the 18 kDa translocator protein (TSPO), whose role in the regulation of membrane permeability is not completely understood. We show that TSPO ligands, 1 μM PPIX and PK11195 at concentrations of 50 μM, accelerate opening of permeability transition pores (mPTP) in Ca(2+)-overloaded rat brain mitochondria (RBM). By contrast, PK11195 at 100 nM and anti-TSPO antibodies suppressed pore opening. Participation of VDAC in these processes was demonstrated by blocking VDAC with G3139, an 18-mer phosphorothioate oligonucleotides, which sensitized mitochondria to Ca(2+)-induced mPTP opening. Despite the inhibitory effect of 100 nM PK11195 and anti-TSPO antibodies alone, their combination with G3139 considerably stimulated the mPTP opening. Thus, 100 nM PK11195 and anti-TSPO antibody can modify permeability of the VDAC channel and mPTP. When VDAC channels are closed and TSPO is blocked, permeability of the VDAC for calcium seems to be the highest, which leads to accelerated pore opening.

Dual mechanism of ion permeation through VDAC revealed with inorganic phosphate ions and phosphate metabolites

In the exchange of metabolites and ions between the mitochondrion and the cytosol, the voltage-dependent anion channel (VDAC) is a key element, as it forms the major transport pathway for these compounds through the mitochondrial outer membrane. Numerous experimental studies have promoted the idea that VDAC acts as a regulator of essential mitochondrial functions. In this study, using a combination of molecular dynamics simulations, free-energy calculations, and electrophysiological measurements, we investigated the transport of ions through VDAC, with a focus on phosphate ions and metabolites. We showed that selectivity of VDAC towards small anions including monovalent phosphates arises from short-lived interactions with positively charged residues scattered throughout the pore. In dramatic contrast, permeation of divalent phosphate ions and phosphate metabolites (AMP and ATP) involves binding sites along a specific translocation pathway. This permeation mechanism offers an explanation for the decrease in VDAC conductance measured in the presence of ATP or AMP at physiological salt concentration. The binding sites occur at similar locations for the divalent phosphate ions, AMP and ATP, and contain identical basic residues. ATP features a marked affinity for a central region of the pore lined by two lysines and one arginine of the N-terminal helix. This cluster of residues together with a few other basic amino acids forms a "charged brush" which facilitates the passage of the anionic metabolites through the pore. All of this reveals that VDAC controls the transport of the inorganic phosphates and phosphate metabolites studied here through two different mechanisms.

The influence of cholesterol on membrane protein structure, function, and dynamics studied by molecular dynamics simulations

The plasma membrane, which encapsulates human cells, is composed of a complex mixture of lipids and embedded proteins. Emerging knowledge points towards the lipids as having a regulating role in protein function. Furthermore, insight from protein crystallography has revealed several different types of lipids intimately bound to membrane proteins and peptides, hereby possibly pointing to a site of action for the observed regulation. Cholesterol is among the lipid membrane constituents most often observed to be co-crystallized with membrane proteins, and the cholesterol levels in cell membranes have been found to play an essential role in health and disease. Remarkably little is known about the mechanism of lipid regulation of membrane protein function in health as well as in disease. Herein, we review molecular dynamics simulation studies aimed at investigating the effect of cholesterol on membrane protein and peptide properties. This article is part of a Special Issue entitled: Lipid-protein interactions.