Search

Ferroptosis is an iron dependent form of cell death, that is triggered by the discoordination of iron, lipids, and thiols. Its unique signature that distinguishes it from other forms of cell death is the formation and accumulation of lipid hydroperoxides, particularly oxidized forms of polyunsaturated phosphatidylethanolamines (PEs), which drives cell death. These readily undergo iron-catalyzed secondary free radical reactions leading to truncated products which retain the signature PE headgroup and which can readily react with nucleophilic moieties in proteins via their truncated electrophilic acyl chains. Using a redox lipidomics approach, we have identified oxidatively-truncated PE species (trPEox) in enzymatic and non-enzymatic model systems. Further, using a model peptide we demonstrate adduct formation with Cys as the preferred nucleophilic residue and PE(26:2) +2 oxygens, as one of the most reactive truncated PE-electrophiles produced. In cells stimulated to undergo ferroptosis we identified PE-truncated species with sn-2 truncations ranging from 5 to 9 carbons. Taking advantage of the free PE headgroup, we have developed a new technology using the lantibiotic duramycin, to enrich and identify the PE-lipoxidated proteins. Our results indicate that several dozens of proteins for each cell type, are PE-lipoxidated in HT-22, MLE, and H9c2 cells and M2 macrophages after they were induced to undergo ferroptosis. Pretreatment of cells with the strong nucleophile, 2-mercaptoethanol, prevented the formation of PE-lipoxidated proteins and blocked ferroptotic death. Finally, our docking simulations showed that the truncated PE species bound at least as good to several of the lantibiotic-identified proteins, as compared to the non-truncated parent molecule, stearoyl-arachidonoyl PE (SAPE), indicating that these oxidatively-truncated species favor/promote the formation of PEox-protein adducts. The identification of PEox-protein adducts during ferroptosis suggests that they are participants in the ferroptotic process preventable by 2-mercaptoethanol and may contribute to a point of no return in the ferroptotic death process., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

Ferroptosis and necroptosis are two pro-inflammatory cell death programs contributing to major pathologies and their inhibition has gained attention to treat a wide range of disease states. Necroptosis relies on activation of RIP1 and RIP3 kinases. Ferroptosis is triggered by oxidation of polyunsaturated phosphatidylethanolamines (PUFA-PE) by complexes of 15-Lipoxygenase (15LOX) with phosphatidylethanolamine-binding protein 1 (PEBP1). The latter, also known as RAF kinase inhibitory protein, displays promiscuity towards multiple proteins. In this study we show that RIP3 K51A kinase inactive mice have increased ferroptotic burden and worse outcome after irradiation and brain trauma rescued by anti-ferroptotic compounds Liproxstatin-1 and Ferrostatin 16-86. Given structural homology between RAF and RIP3, we hypothesized that PEBP1 acts as a necroptosis-to-ferroptosis switch interacting with either RIP3 or 15LOX. Using genetic, biochemical, redox lipidomics and computational approaches, we uncovered that PEBP1 complexes with RIP3 and inhibits necroptosis. Elevated expression combined with higher affinity enables 15LOX to pilfer PEBP1 from RIP3, thereby promoting PUFA-PE oxidation and ferroptosis which sensitizes Rip3 K51A/K51A kinase-deficient mice to total body irradiation and brain trauma. This newly unearthed PEBP1/15LOX-driven mechanism, along with previously established switch between necroptosis and apoptosis, can serve multiple and diverse cell death regulatory functions across various human disease states., (Copyright © 2022. Published by Elsevier B.V.)

We recently discovered an anti-ferroptotic mechanism inherent to M1 macrophages whereby high levels of NO ● suppressed ferroptosis via inhibition of hydroperoxy-eicosatetraenoyl-phosphatidylethanolamine (HpETE-PE) production by 15-lipoxygenase (15LOX) complexed with PE-binding protein 1 (PEBP1). However, the mechanism of NO ● interference with 15LOX/PEBP1 activity remained unclear. Here, we use a biochemical model of recombinant 15LOX-2 complexed with PEBP1, LC-MS redox lipidomics, and structure-based modeling and simulations to uncover the mechanism through which NO ● suppresses ETE-PE oxidation. Our study reveals that O 2 and NO ● use the same entry pores and channels connecting to 15LOX-2 catalytic site, resulting in a competition for the catalytic site. We identified residues that direct O 2 and NO ● to the catalytic site, as well as those stabilizing the esterified ETE-PE phospholipid tail. The functional significance of these residues is supported by in silico saturation mutagenesis. We detected nitrosylated PE species in a biochemical system consisting of 15LOX-2/PEBP1 and NO ● donor and in RAW264.7 M2 macrophages treated with ferroptosis-inducer RSL3 in the presence of NO ● , in further support of the ability of NO ● to diffuse to, and react at, the 15LOX-2 catalytic site. The results provide first insights into the molecular mechanism of repression of the ferroptotic Hp-ETE-PE production by NO ● .

Ferroptosis, triggered by discoordination of iron, thiols and lipids, leads to the accumulation of 15-hydroperoxy (Hp)-arachidonoyl-phosphatidylethanolamine (15-HpETE-PE), generated by complexes of 15-lipoxygenase (15-LOX) and a scaffold protein, phosphatidylethanolamine (PE)-binding protein (PEBP)1. As the Ca 2+ -independent phospholipase A 2 β (iPLA 2 β, PLA2G6 or PNPLA9 gene) can preferentially hydrolyze peroxidized phospholipids, it may eliminate the ferroptotic 15-HpETE-PE death signal. Here, we demonstrate that by hydrolyzing 15-HpETE-PE, iPLA 2 β averts ferroptosis, whereas its genetic or pharmacological inactivation sensitizes cells to ferroptosis. Given that PLA2G6 mutations relate to neurodegeneration, we examined fibroblasts from a patient with a Parkinson's disease (PD)-associated mutation (fPD R747W ) and found selectively decreased 15-HpETE-PE-hydrolyzing activity, 15-HpETE-PE accumulation and elevated sensitivity to ferroptosis. CRISPR-Cas9-engineered Pnpla9 R748W/R748W mice exhibited progressive parkinsonian motor deficits and 15-HpETE-PE accumulation. Elevated 15-HpETE-PE levels were also detected in midbrains of rotenone-infused parkinsonian rats and α-synuclein-mutant Snca A53T mice, with decreased iPLA 2 β expression and a PD-relevant phenotype. Thus, iPLA 2 β is a new ferroptosis regulator, and its mutations may be implicated in PD pathogenesis.

A huge diversification of phospholipids, forming the aqueous interfaces of all biomembranes, cannot be accommodated within a simple concept of their role as membrane building blocks. Indeed, a number of signaling functions of (phospho)lipid molecules has been discovered. Among these signaling lipids, a particular group of oxygenated polyunsaturated fatty acids (PUFA), so called lipid mediators, has been thoroughly investigated over several decades. This group includes oxygenated octadecanoids, eicosanoids, and docosanoids and includes several hundreds of individual species. Oxygenation of PUFA can occur when they are esterified into major classes of phospholipids. Initially, these events have been associated with non-specific oxidative injury of biomembranes. An alternative concept is that these post-synthetically oxidatively modified phospholipids and their adducts with proteins are a part of a redox epiphospholipidome that represents a rich and versatile language for intra- and inter-cellular communications. The redox epiphospholipidome may include hundreds of thousands of individual molecular species acting as meaningful biological signals. This review describes the signaling role of oxygenated phospholipids in programs of regulated cell death. Although phospholipid peroxidation has been associated with almost all known cell death programs, we chose to discuss enzymatic pathways activated during apoptosis and ferroptosis and leading to peroxidation of two phospholipid classes, cardiolipins (CLs) and phosphatidylethanolamines (PEs). This is based on the available LC-MS identification and quantitative information on the respective peroxidation products of CLs and PEs. We focused on molecular mechanisms through which two proteins, a mitochondrial hemoprotein cytochrome c (cyt c ), and non-heme Fe lipoxygenase (LOX), change their catalytic properties to fulfill new functions of generating oxygenated CL and PE species. Given the high selectivity and specificity of CL and PE peroxidation we argue that enzymatic reactions catalyzed by cyt c /CL complexes and 15-lipoxygenase/phosphatidylethanolamine binding protein 1 (15LOX/PEBP1) complexes dominate, at least during the initiation stage of peroxidation, in apoptosis and ferroptosis. We contrast cell-autonomous nature of CLox signaling in apoptosis correlating with its anti-inflammatory functions vs. non-cell-autonomous ferroptotic signaling facilitating pro-inflammatory (necro-inflammatory) responses. Finally, we propose that small molecule mechanism-based regulators of enzymatic phospholipid peroxidation may lead to highly specific anti-apoptotic and anti-ferroptotic therapeutic modalities., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Kagan, Tyurina, Vlasova, Kapralov, Amoscato, Anthonymuthu, Tyurin, Shrivastava, Cinemre, Lamade, Epperly, Greenberger, Beezhold, Mallampalli, Srivastava, Bayir and Shvedova.)

Hydroperoxy-eicosatetraenoyl-phosphatidylethanolamine (HpETE-PE) is a ferroptotic cell death signal. HpETE-PE is produced by the 15-Lipoxygenase (15LOX)/Phosphatidylethanolamine Binding Protein-1 (PEBP1) complex or via an Fe-catalyzed non-enzymatic radical reaction. Ferrostatin-1 (Fer-1), a common ferroptosis inhibitor, is a lipophilic radical scavenger but a poor 15LOX inhibitor arguing against 15LOX having a role in ferroptosis. In the current work, we demonstrate that Fer-1 does not affect 15LOX alone, however, it effectively inhibits HpETE-PE production by the 15LOX/PEBP1 complex. Computational molecular modeling shows that Fer-1 binds to the 15LOX/PEBP1 complex at three sites and could disrupt the catalytically required allosteric motions of the 15LOX/PEBP1 complex. Using nine ferroptosis cell/tissue models, we show that HpETE-PE is produced by the 15LOX/PEBP1 complex and resolve the long-existing Fer-1 anti-ferroptotic paradox., (Copyright © 2020 The Authors. Published by Elsevier B.V. All rights reserved.)

Cardiotoxin CTII from Naja oxiana cobra venom translocates to the intermembrane space (IMS) of mitochondria to disrupt the structure and function of the inner mitochondrial membrane. At low concentrations, CTII facilitates ATP-synthase activity, presumably via the formation of non-bilayer, immobilized phospholipids that are critical in modulating ATP-synthase activity. In this study, we investigated the effects of another cardiotoxin CTI from Naja oxiana cobra venom on the structure of mitochondrial membranes and on mitochondrial-derived ATP synthesis. By employing robust biophysical methods including 31 P-NMR and 1 H-NMR spectroscopy, we analyzed the effects of CTI and CTII on phospholipid packing and dynamics in model phosphatidylcholine (PC) membranes enriched with 2.5 and 5.0 mol% of cardiolipin (CL), a phospholipid composition that mimics that in the outer mitochondrial membrane (OMM). These experiments revealed that CTII converted a higher percentage of bilayer phospholipids to a non-bilayer and immobilized state and both cardiotoxins utilized CL and PC molecules to form non-bilayer structures. Furthermore, in order to gain further understanding on how cardiotoxins bind to mitochondrial membranes, we employed molecular dynamics (MD) and molecular docking simulations to investigate the molecular mechanisms by which CTII and CTI interactively bind with an in silico phospholipid membrane that models the composition similar to the OMM. In brief, MD studies suggest that CTII utilized the N-terminal region to embed the phospholipid bilayer more avidly in a horizontal orientation with respect to the lipid bilayer and thereby penetrate at a faster rate compared with CTI. Molecular dynamics along with the Autodock studies identified critical amino acid residues on the molecular surfaces of CTII and CTI that facilitated the long-range and short-range interactions of cardiotoxins with CL and PC. Based on our compiled data and our published findings, we provide a conceptual model that explains a molecular mechanism by which snake venom cardiotoxins, including CTI and CTII, interact with mitochondrial membranes to alter the mitochondrial membrane structure to either upregulate ATP-synthase activity or disrupt mitochondrial function., Competing Interests: The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Temporally harmonized elimination of damaged or unnecessary organelles and cells is a prerequisite of health. Under Type 2 inflammatory conditions, human airway epithelial cells (HAECs) generate proferroptotic hydroperoxy-arachidonoyl-phosphatidylethanolamines (HpETE-PEs) as proximate death signals. Production of 15-HpETE-PE depends on activation of 15-lipoxygenase-1 (15LO1) in complex with PE-binding protein-1 (PEBP1). We hypothesized that cellular membrane damage induced by these proferroptotic phospholipids triggers compensatory prosurvival pathways, and in particular autophagic pathways, to prevent cell elimination through programmed death. We discovered that PEBP1 is pivotal to driving dynamic interactions with both proferroptotic 15LO1 and the autophagic protein microtubule-associated light chain-3 (LC3). Further, the 15LO1-PEBP1-generated ferroptotic phospholipid, 15-HpETE-PE, promoted LC3-I lipidation to stimulate autophagy. This concurrent activation of autophagy protects cells from ferroptotic death and release of mitochondrial DNA. Similar findings are observed in Type 2 Hi asthma, where high levels of both 15LO1-PEBP1 and LC3-II are seen in HAECs, in association with low bronchoalveolar lavage fluid mitochondrial DNA and more severe disease. The concomitant activation of ferroptosis and autophagy by 15LO1-PEBP1 complexes and their hydroperoxy-phospholipids reveals a pathobiologic pathway relevant to asthma and amenable to therapeutic targeting., Competing Interests: The authors declare no competing interest.

Ferroptotic death is the penalty for losing control over three processes-iron metabolism, lipid peroxidation and thiol regulation-that are common in the pro-inflammatory environment where professional phagocytes fulfill their functions and yet survive. We hypothesized that redox reprogramming of 15-lipoxygenase (15-LOX) during the generation of pro-ferroptotic signal 15-hydroperoxy-eicosa-tetra-enoyl-phosphatidylethanolamine (15-HpETE-PE) modulates ferroptotic endurance. Here, we have discovered that inducible nitric oxide synthase (iNOS)/NO • -enrichment of activated M1 (but not alternatively activated M2) macrophages/microglia modulates susceptibility to ferroptosis. Genetic or pharmacologic depletion/inactivation of iNOS confers sensitivity on M1 cells, whereas NO • donors empower resistance of M2 cells to ferroptosis. In vivo, M1 phagocytes, in comparison to M2 phagocytes, exert higher resistance to pharmacologically induced ferroptosis. This resistance is diminished in iNOS-deficient cells in the pro-inflammatory conditions of brain trauma or the tumour microenvironment. The nitroxygenation of eicosatetraenoyl (ETE)-PE intermediates and oxidatively truncated species by NO • donors and/or suppression of NO • production by iNOS inhibitors represent a novel redox mechanism of regulation of ferroptosis in pro-inflammatory conditions.

Cobra venom cardiotoxins (CVCs) can translocate to mitochondria to promote apoptosis by eliciting mitochondrial dysfunction. However, the molecular mechanism(s) by which CVCs are selectively targeted to the mitochondrion to disrupt mitochondrial function remains to be elucidated. By studying cardiotoxin from Naja mossambica mossambica cobra (cardiotoxin VII4), a basic three-fingered S-type cardiotoxin, we hypothesized that cardiotoxin VII4 binds to cardiolipin (CL) in mitochondria to alter mitochondrial structure/function and promote neurotoxicity. By performing confocal analysis, we observed that red-fluorescently tagged cardiotoxin rapidly translocates to mitochondria in mouse primary cortical neurons and in human SH-SY5Y neuroblastoma cells to promote aberrant mitochondrial fragmentation, a decline in oxidative phosphorylation, and decreased energy production. In addition, by employing electron paramagnetic resonance (EPR) and protein nuclear magnetic resonance (¹H-NMR) spectroscopy and phosphorescence quenching of erythrosine in model membranes, our compiled biophysical data show that cardiotoxin VII4 binds to anionic CL, but not to zwitterionic phosphatidylcholine (PC), to increase the permeability and formation of non-bilayer structures in CL-enriched membranes that biochemically mimic the outer and inner mitochondrial membranes. Finally, molecular dynamics simulations and in silico docking studies identified CL binding sites in cardiotoxin VII4 and revealed a molecular mechanism by which cardiotoxin VII4 interacts with CL and PC to bind and penetrate mitochondrial membranes.

Ferroptosis is a form of programmed cell death that is pathogenic to several acute and chronic diseases and executed via oxygenation of polyunsaturated phosphatidylethanolamines (PE) by 15-lipoxygenases (15-LO) that normally use free polyunsaturated fatty acids as substrates. Mechanisms of the altered 15-LO substrate specificity are enigmatic. We sought a common ferroptosis regulator for 15LO. We discovered that PEBP1, a scaffold protein inhibitor of protein kinase cascades, complexes with two 15LO isoforms, 15LO1 and 15LO2, and changes their substrate competence to generate hydroperoxy-PE. Inadequate reduction of hydroperoxy-PE due to insufficiency or dysfunction of a selenoperoxidase, GPX4, leads to ferroptosis. We demonstrated the importance of PEBP1-dependent regulatory mechanisms of ferroptotic death in airway epithelial cells in asthma, kidney epithelial cells in renal failure, and cortical and hippocampal neurons in brain trauma. As master regulators of ferroptotic cell death with profound implications for human disease, PEBP1/15LO complexes represent a new target for drug discovery., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Sterile liver inflammation, such as liver ischemia-reperfusion, hemorrhagic shock after trauma, and drug-induced liver injury, is initiated and regulated by endogenous mediators including DNA and reactive oxygen species. Here, we identify a mechanism for redox-mediated regulation of absent in melanoma 2 (AIM2) inflammasome activation in hepatocytes after redox stress in mice, which occurs through interaction with cytosolic high mobility group box 1 (HMGB1). We show that in liver during hemorrhagic shock in mice and in hepatocytes after hypoxia with reoxygenation, cytosolic HMGB1 associates with AIM2 and is required for activation of caspase-1 in response to cytosolic DNA. Activation of caspase-1 through AIM2 leads to subsequent hepatoprotective responses such as autophagy. HMGB1 binds to AIM2 at a non-DNA-binding site on the hematopoietic interferon-inducible nuclear antigen domain of AIM2 to facilitate inflammasome and caspase-1 activation in hepatocytes. Furthermore, binding of HMGB1 to AIM2 is stronger with fully reduced all-thiol HMGB1 than with partially oxidized disulfide-HMGB1, and binding strength corresponds to caspase-1 activation. These data suggest that HMGB1 redox status regulates AIM2 inflammasome activation., Conclusion: These findings suggest a novel and important mechanism for regulation of AIM2 inflammasome activation in hepatocytes during redox stress and may suggest broader implications for how this and other inflammasomes are activated and how their activation is regulated during cell stress, as well as the mechanisms of inflammasome regulation in nonimmune cell types. (Hepatology 2017;65:253-268)., (© 2016 by the American Association for the Study of Liver Diseases.)

Cobra venom cytotoxins are basic three-fingered, amphipathic, non-enzymatic proteins that constitute a major fraction of cobra venom. While cytotoxins cause mitochondrial dysfunction in different cell types, the mechanisms by which cytotoxins bind to mitochondria remain unknown. We analyzed the abilities of CTI and CTII, S-type and P-type cytotoxins from Naja naja oxiana respectively, to associate with isolated mitochondrial fractions or with model membranes that simulate the mitochondrial lipid environment by using a myriad of biophysical techniques. Phosphorus-31 nuclear magnetic resonance (31P-NMR) spectroscopy data suggest that both cytotoxins bind to isolated mitochondrial fractions and promote the formation of aberrant non-bilayer structures. We then hypothesized that CTI and CTII bind to cardiolipin (CL) to disrupt mitochondrial membranes. Collectively, 31P-NMR, electron paramagnetic resonance (EPR), proton NMR (1H-NMR), deuterium NMR (2H-NMR) spectroscopy, differential scanning calorimetry, and erythrosine phosphorescence assays suggest that CTI and CTII bind to CL to generate non-bilayer structures and promote the permeabilization, dehydration and fusion of large unilamellar phosphatidylcholine (PC) liposomes enriched with CL. On the other hand, CTII but not CTI caused biophysical alterations of large unilamellar PC liposomes enriched with phosphatidylserine (PS). Mechanistically, single molecule docking simulations identified putative CL, PS and PC binding sites in CTI and CTII. While the predicted binding sites for PS and PC share a high number of interactive amino acid residues in CTI and CTII, the CL biding sites in CTII and CTI are more divergent as it contains additional interactive amino acid residues. Overall, our data suggest that cytotoxins physically associate with mitochondrial membranes by binding to CL to disrupt mitochondrial structural integrity.

Recognition of injured mitochondria for degradation by macroautophagy is essential for cellular health, but the mechanisms remain poorly understood. Cardiolipin is an inner mitochondrial membrane phospholipid. We found that rotenone, staurosporine, 6-hydroxydopamine and other pro-mitophagy stimuli caused externalization of cardiolipin to the mitochondrial surface in primary cortical neurons and SH-SY5Y cells. RNAi knockdown of cardiolipin synthase or of phospholipid scramblase-3, which transports cardiolipin to the outer mitochondrial membrane, decreased the delivery of mitochondria to autophagosomes. Furthermore, we found that the autophagy protein microtubule-associated-protein-1 light chain 3 (LC3), which mediates both autophagosome formation and cargo recognition, contains cardiolipin-binding sites important for the engulfment of mitochondria by the autophagic system. Mutation of LC3 residues predicted as cardiolipin-interaction sites by computational modelling inhibited its participation in mitophagy. These data indicate that redistribution of cardiolipin serves as an 'eat-me' signal for the elimination of damaged mitochondria from neuronal cells.

Ionotropic glutamate receptors (iGluRs) harbor two extracellular domains: the membrane-proximal ligand-binding domain (LBD) and the distal N-terminal domain (NTD). These are involved in signal sensing: the LBD binds L-glutamate, which activates the receptor channel. Ligand binding to the NTD modulates channel function in the NMDA receptor subfamily of iGluRs, which has not been observed for the AMPAR subfamily to date. Structural data suggest that AMPAR NTDs are packed into tight dimers and have lost their signaling potential. Here, we assess NTD dynamics from both subfamilies, using a variety of computational tools. We describe the conformational motions that underly NMDAR NTD allosteric signaling. Unexpectedly, AMPAR NTDs are capable of undergoing similar dynamics; although dimerization imposes restrictions, the two subfamilies sample similar, interconvertible conformational subspaces. Finally, we solve the crystal structure of AMPAR GluA4 NTD, and combined with molecular dynamics simulations, we characterize regions pivotal for an as-yet-unexplored dynamic spectrum of AMPAR NTDs., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

Binding of the viral spike drives cell entry and infection by HIV-1 to the cellular CD4 and chemokine receptors with associated conformational change of the viral glycoprotein envelope, gp120. Crystal structures of the CD4-gp120-antibody ternary complex reveal a large internal gp120 cavity formed by three domains-the inner domain, outer domain, and bridging sheet domain-and are capped by CD4 residue Phe43. Several structures of gp120 envelope in complex with various antibodies indicated that the bridging sheet adopts varied conformations. Here, we examine bridging sheet dynamics using a crystal structure of gp120 bound to the F105 antibody exhibiting an open bridging sheet conformation and with an added V3 loop. The two strands of the bridging sheet β2/β3 and β20/β21 are dissociated from each other and are directed away from the inner and outer domains. Analysis of molecular dynamics (MD) trajectories indicates that the β2/β3 and β20/β21 strands rapidly rearrange to interact with the V3 loop and the inner and outer domains, respectively. Residue N425 on β20 leads the conformational rearrangement of the β20/β21 strands by interacting with W112 on the inner domain and F382 on the outer domain. An accompanying shift is observed in the inner domain as helix α1 exhibits a loss in helicity and pivots away from helix α5. The two simulations provide a framework for understanding the conformational diversity of the bridging sheet and the propensity of the β20/β21 strand to refold between the inner and outer domains of gp120, in the absence of a bound ligand.

Glutamate transporters clear synaptically released glutamate to maintain precise communication between neurons and limit glutamate neurotoxicity. Although much progress has been made on the topology, structure, and function of these carriers, few studies have addressed large-scale structural motions collectively associated with substrate transport. Here we show that a series of single cysteine substitutions in the helical hairpin HP2 of excitatory amino acid transporter 1 form intersubunit disulfide cross-links within the trimer. After cross-linking, substrate uptake, but not substrate-activated anion conductance, is completely inhibited in these mutants. These disulfide bridges link residue pairs > 40 Å apart in the outward-facing crystal structure, and can be explained by concerted subunit movements predicted by the anisotropic network model (ANM). The existence of these global motions is further supported by the observation that single cysteine substitutions at the extracellular part of the transmembrane domain 8 can also be cross-linked by copper phenanthroline as predicted by the ANM. Interestingly, the transport domain in the un-cross-linked subunit of the trimer assumes an inward-facing orientation, suggesting that individual subunits potentially undergo separate transitions between outward- and inward-facing forms, rather than an all-or-none transition of the three subunits, a mechanism also supported by ANM-predicted intrinsic dynamics. These results shed light on how large collective motions contribute to the functional dynamics of glutamate transporters.

Glutamate transporters regulate excitatory amino acid neurotransmission across neuronal and glial cell membranes by coupling the translocation of their substrate (aspartate or glutamate) into the intracellular (IC) medium to the energetically favorable transport of sodium ions or other cations. The first crystallographically resolved structure of this family, the archaeal aspartate transporter, Glt(Ph), has served as a structural paradigm for elucidating the mechanism of substrate translocation by these transporters. Two helical hairpins, HP2 and HP1, at the core domains of the three subunits that form this membrane protein have been proposed to act as the respective extracellular and IC gates for substrate intake and release. Molecular dynamics simulations using the outward-facing structure have confirmed that the HP2 loop acts as an EC gate. The mechanism of substrate release at atomic scale, however, remained unknown due to the lack of structural data until the recent determination of the inward-facing structure of Glt(Ph). In the present study, we use this recently resolved structure to simulate the release of substrate to the cytoplasm and the roles of HP1 and HP2 in this process. The highly flexible HP2 loop is observed to serve as an activator (or initiator) prompting the release of a gatekeeper Na(+) to the cytoplasm and promoting the influx of water molecules from the cytoplasm, which effectively disrupt substrate-protein interactions and drive the dislodging of the substrate from its binding site. The completion of substrate release and exit, however, entails the opening of the highly stable HP1 loop as well. Overall, the unique conformational flexibility of the HP2 loop, the dissociation of a Na(+), the hydration of binding pocket, and final yielding of the HP1 loop 3-Ser motif emerge as the successive events controlling the release of the bound substrate to the cell interior by glutamate transporters.

Glutamate transporters are membrane proteins found in neurons and glial cells, which play a critical role in regulating cell signaling by clearing glutamate released from synapses. Although extensive biochemical and structural studies have shed light onto different aspects of glutamate transport, the time-resolved molecular mechanism of substrate (glutamate or aspartate) translocation, that is, the sequence of events occurring at the atomic level after substrate binding and before its release intracellularly remain to be elucidated. We identify an energetically preferred permeation pathway of approximately 23 A between the helix HP1b on the hairpin HP1 and the transmembrane helices TM7 and TM8, using the high resolution structure of the transporter from Pyrococcus horikoshii (Glt(Ph)) in steered molecular dynamics simulations. Detailed potential of mean force calculations along the putative pathway reveal 2 energy barriers encountered by the substrate (aspartate) before it reaches the exit. The first barrier is surmounted with the assistance of 2 conserved residues (S278 and N401) and a sodium ion (Na2); and the second, by the electrostatic interactions with D405 and another sodium ion (Na1). The observed critical interactions and mediating role of conserved residues in the core domain, the accompanying conformational changes (in both substrate and transporter) that relieve local strains, and the unique coupling of aspartate transport to Na(+) dislocation provide insights into methods for modulating substrate transport.

Glutamate transporters, also referred to as excitatory amino acid transporters (EAATs), are membrane proteins that regulate glutamatergic signal transmission by clearing excess glutamate after its release at synapses. A structure-based understanding of their molecular mechanisms of function has been elusive until the recent determination of the x-ray structure of an archaeal transporter, Glt(Ph). Glt(Ph) exists as a trimer, with each subunit containing a core region that mediates substrate translocation. In the present study a series of molecular dynamics simulations have been conducted and analyzed in light of new experimental data on substrate binding properties of EAATs. The simulations provide for the first time a full atomic description of the time-resolved events that drive the recognition and binding of substrate. The core region of each subunit exhibits an intrinsic tendency to open the helical hairpin HP2 loop, the extracellular gate, within tens of nanoseconds exposing conserved polar residues that serve as attractors for substrate binding. The NMDGT motif on the partially unwound part of the transmembrane helix TM7 and the residues Asp-390 and Asp-394 on TM8 are also distinguished by their important role in substrate binding and close interaction with mediating water molecules and/or sodium ions. The simulations reveal a Na+ binding site comprised in part of Leu-303 on TM7 and Asp-405 on TM8 and support a role for sodium ions in stabilizing substrate-bound conformers. The functional importance of Leu-303 or its counterpart Leu-391 in human EAAT1 (hEAAT1) is confirmed by site-directed mutagenesis and Na+ dependence assays conducted with hEAAT1 mutants L391C and L391A.

A fundamental question associated with the function of ion channels is the conformational changes that allow for reversibly opening/occluding the pore through which the cations permeate. The recently elucidated crystal structures of potassium channels reveal similar structural motifs at their pore-forming regions, suggesting that they share a common gating mechanism. The validity of this hypothesis is explored by analyzing the collective dynamics of five known K(+) channel structures. Normal-mode analysis using the Gaussian network model strikingly reveals that all five structures display the same intrinsic motions at their pore-forming region despite the differences in their sequences, structures, and activation mechanisms. Superposition of the most cooperative mode profiles shows that the identified common mechanism is a global corkscrew-like counterrotation of the extracellular and cytoplasmic (CP) regions, leading to the opening of the CP end of the pore. A second cooperative mode shared by all five K(+) channels is the extension of the extracellular and/or CP ends via alternating anticorrelated fluctuations of pairs of diagonally opposite monomers. Residues acting as hinges/anchors in both modes are highly conserved across the members of the family of K(+) channel proteins, consistent with their presently disclosed critical mechanical role in pore gating.

In the preceding, accompanying article, we present models of the structure and voltage-dependent gating mechanism of the KvAP bacterial K+ channel that are based on three types of evidence: crystal structures of portions of the KvAP protein, theoretical modeling criteria for membrane proteins, and biophysical studies of the properties of native and mutated voltage-gated channels. Most of the latter experiments were performed on the Shaker K+ channel. Some of these data are difficult to relate directly to models of the KvAP channel's structure due to differences in the Shaker and KvAP sequences. We have dealt with this problem by developing new models of the structure and gating mechanism of the transmembrane and extracellular portions of the Shaker channel. These models are consistent with almost all of the biophysical data. In contrast, much of the experimental data are incompatible with the "paddle" model of gating that was proposed when the KvAP crystal structures were first published. The general folding pattern and gating mechanisms of our current models are similar to some of our earlier models of the Shaker channel., (Copyright 2004 Biophysical Society)

Having inspected the crystal structure of the complete KvAP channel protein, we suspect that the voltage-sensing domain is too distorted to provide reliable information about its native tertiary structure or its interactions with the central pore-forming domain. On the other hand, a second crystal structure of the isolated voltage-sensing domain may well correspond to a native open conformation. We also observe that the paddle model of gating developed from these two structures is inconsistent with many experimental results, and suspect it to be energetically unrealistic. Here we show that the isolated voltage-sensing domain crystal structure can be docked onto the pore domain portion of the full-length KvAP crystal structure in an energetically favorable way to create a model of the open conformation. Using this as a starting point, we have developed rather conventional models of resting and transition conformations based on the helical screw mechanism for the transition from the open to the resting conformation. Our models are consistent with both theoretical considerations and experimental results., (Copyright 2004 Biophysical Society)

GluR0 is a prokaryotic homologue of mammalian glutamate receptors that forms glutamate-activated, potassium-selective ion channels. The topology of its transmembrane (TM) domain is similar to that of simple potassium channels such as KcsA. Two plausible alignments of the sequence of the TM domain of GluR0 with KcsA are possible, differing in the region of the P helix. We have constructed homology models based on both alignments and evaluated them using 6 ns duration molecular dynamics simulations in a membrane-mimetic environment. One model, in which an insertion in GluR0 relative to KcsA is located in the loop between the M1 and P helices, is preferred on the basis of lower structural drift and maintenance of the P helix conformation during simulation. This model also exhibits inter-subunit salt bridges that help to stabilise the TM domain tetramer. During the simulation, concerted K(+) ion-water movement along the selectivity filter is observed, as is the case in simulations of KcsA. K(+) ion exit from the central cavity is associated with opening of the hydrophobic gate formed by the C-termini of the M2 helices. In the intact receptor the opening of this gate will be controlled by interactions with the extramembranous ligand-binding domains.

Potassium channels have been studied intensively in terms of the relationship between molecular structure and physiological function. They provide an opportunity to integrate structural and computational studies in order to arrive at an atomic resolution description of mechanism. We review recent progress in K channel structural studies, focussing on the bacterial channel KcsA. Structural studies can be extended via use of computational (i.e. molecular simulation) approaches in order to provide a perspective on aspects of channel function such as permeation, selectivity, block and gating. Results from molecular dynamics simulations are shown to be in good agreement with recent structural studies of KcsA in terms of the interactions of K(+) ions with binding sites within the selectivity filter of the channel, and in revealing the importance of filter flexibility in channel function. We discuss how the KcsA structure may be used as a template for developing structural models of other families of K channels. Progress in this area is explored via two examples: inward rectifier (Kir) and voltage-gated (Kv) potassium channels. A brief account of structural studies of ancillary domains and subunits of K channels is provided.

Prolines in transmembrane (TM) alpha-helices are believed to play an important structural and/or functional role in membrane proteins. At a structural level a proline residue distorts alpha-helical structure due to the loss of at least one stabilizing backbone hydrogen bond, and introduces flexibility in the helix that may result in substantial kink and swivel motions about the effective "hinge." At a functional level, for example in Kv channels, it is believed that proline-induced molecular hinges may have a direct role in gating, i.e., the conformational change linked to opening/closing the channel to movement of ions. In this article we study the conformational dynamics of the S6 TM helix from of the Kv channel Shaker, which possesses the motif PVP--a motif that is conserved in Kv channels. We perform multiple molecular dynamics simulations of single S6 helices in a membrane-mimetic environment in order to effectively map the kink-swivel conformational space of the protein, exploiting the ability of multiple simulations to achieve greater sampling. We show that the presence of proline locally perturbs the helix, disrupting local dihedral angles and producing local twist and unwinding in the region of the hinge--an effect that is relaxed with distance from the PVP motif. We furthermore show that motions about the hinge are highly anisotropic, reflecting a preferred region of kink-swivel conformation space that may have implications for the gating process., (Copyright 2002 Wiley Periodicals, Inc. Biopolymers 64: 303-313, 2002)

Molecular dynamics simulations of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal significant differences in interactions of the selectivity filter with K(+) compared with Na(+) ions. K(+) ions and water molecules within the filter undergo concerted single-file motion in which they translocate between adjacent sites within the filter on a nanosecond timescale. In contrast, Na(+) ions remain bound to sites within the filter and do not exhibit translocation on a nanosecond timescale. Furthermore, entry of a K(+) ion into the filter from the extracellular mouth is observed, whereas this does not occur for a Na(+) ion. Whereas K(+) ions prefer to sit within a cage of eight oxygen atoms of the filter, Na(+) ions prefer to interact with a ring of four oxygen atoms plus two water molecules. These differences in interactions in the selectivity filter may contribute to the selectivity of KcsA for K(+) ions (in addition to the differences in dehydration energy between K(+) and Na(+)) and the block of KcsA by internal Na(+) ions. In our simulations the selectivity filter exhibits significant flexibility in response to changes in ion/protein interactions, with a somewhat greater distortion induced by Na(+) than by K(+) ions.

The gating mechanism of a bacterial potassium channel, KcsA, has been investigated via multi-nanosecond molecular dynamic simulations of the channel molecules embedded in a fully solvated palmitoyloleoylphosphatidylcholine bilayer. Four events are seen in which a cation (K(+) or, in one case, Na(+)) initially present in the central cavity exits through the intracellular mouth (the presumed gate) of the channel. Whilst in the cavity a cation interacts with the sidechain T107 O gamma atom of one of the subunits prior to its exit from the channel. Secondary structure analysis as a function of time reveals a break in the helicity of one of the M2 helices. This break is expected to lend flexibility to the helices, enabling them to "open" (minimum pore radius >0.13 nm) and "close" (minimum pore radius <0.13 nm) the channel. Fluctuations in the pore radius at the intracellular gate region are of the order of 0.05 nm, with an average radius in the region of the gate of ca. 0.1 nm. However, around the time of exit of a cation, the pore widens to about 0.15 nm. The distances between the C alpha atoms of the inner helices M2 reveal a coupled increase and decrease between the opposite pair of helices at about the time of exit of the ion. This suggests a breathing motion of the M2 helices that may form the basis for a gating mechanism.

Our understanding of ion permeation through K(+) channels, and by extension through other channels, is advancing rapidly. New structural studies, together with computer simulations, have provided profound insights into ion conduction mechanisms.

A number of ion channels contain transmembrane (TM) alpha-helices that contain proline-induced molecular hinges. These TM helices include the channel-forming peptide alamethicin (Alm), the S6 helix from voltage-gated potassium (Kv) channels, and the D5 helix from voltage-gated chloride (CLC) channels. For both Alm and KvS6, experimental data implicate hinge-bending motions of the helix in an aspect of channel gating. We have compared the hinge-bending motions of these TM helices in bilayer-like environments by multi-nanosecond MD simulations in an attempt to describe motions of these helices that may underlie possible modes of channel gating. Alm is an alpha-helical channel-forming peptide, which contains a central kink associated with a Gly-x-x-Pro motif in its sequence. Simulations of Alm in a TM orientation for 10 ns in an octane slab indicate that the Gly-x-x-Pro motif acts as a molecular hinge. The S6 helix from Shaker Kv channels contains a Pro-Val-Pro motif. Modeling studies and recent experimental data suggest that the KvS6 helix may be kinked in the vicinity of this motif. Simulations (10 ns) of an isolated KvS6 helix in an octane slab and in a POPC bilayer reveal hinge-bending motions. A pattern-matching approach was used to search for possible hinge-bending motifs in the TM helices of other ion channel proteins. This uncovered a conserved Gly-x-Pro motif in TM helix D5 of CLC channels. MD simulations of a model of hCLC1-D5 spanning an octane slab suggest that this channel also contains a TM helix that undergoes hinge-bending motion. In conclusion, our simulations suggest a model in which hinge-bending motions of TM helices may play a functional role in the gating mechanisms of several different families of ion channels., (Copyright 2001 Wiley-Liss, Inc.)

KcsA is a bacterial K+ channel that is gated by pH. Continuum dielectric calculations on the crystal structure of the channel protein embedded in a low dielectric slab suggest that side chains E71 and D80 of each subunit, which lie adjacent to the selectivity filter region of the channel, form a proton-sharing pair in which E71 is neutral (protonated) and D80 is negatively charged at pH 7. When K+ ions are introduced into the system at their crystallographic positions the pattern of proton sharing is altered. The largest perturbation is for a K+ ion at site S3, i.e., interacting with the carbonyls of T75 and V76. The presence of multiple K+ ions in the filter increases the probability of E71 being ionized and of D80 remaining neutral (i.e., protonated). The ionization states of the protein side chains influence the potential energy profile experienced by a K+ ion as it is translated along the pore axis. In particular, the ionization state of the E71-D80 proton-sharing pair modulates the shape of the potential profile in the vicinity of the selectivity filter. Such reciprocal effects of ion occupancy on side-chain ionization states, and of side-chain ionization states on ion potential energy profiles will complicate molecular dynamics simulations and related studies designed to calculate ion permeation energetics.

Ion channels mediate electrical excitability in neurons and muscle. Three-dimensional structures for model peptide channels and for a potassium (K+) channel have been combined with computer simulations to permit rigorous exploration of structure-function relations of channels. Water molecules and ions within transbilayer pores tend to diffuse more slowly than in bulk solutions. In the narrow selectivity filter of the bacterial K+ channel (i.e. the region of the channel that discriminates between different species of ions) a column of water molecules and K+ ions moves in a concerted fashion. By combining atomistic simulations (in which all atoms of the channel molecule, water and ions are treated explicitly) with continuum methods (in which the description of the channel system is considerably simplified) it is possible to simulate some of the physiological properties of channels.

A homology model has been generated for the pore-forming domain of Kir6.2, a component of an ATP-sensitive K channel, based on the x-ray structure of the bacterial channel KcsA. Analysis of the lipid-exposed and pore-lining surfaces of the model reveals them to be compatible with the known features of membrane proteins and Kir channels, respectively. The Kir6.2 homology model was used as the starting point for nanosecond-duration molecular dynamics simulations in a solvated phospholipid bilayer. The overall drift from the model structure was comparable to that seen for KcsA in previous similar simulations. Preliminary analysis of the interactions of the Kir6.2 channel model with K(+) ions and water molecules during these simulations suggests that concerted single-file motion of K(+) ions and water through the selectivity filter occurs. This is similar to such motion observed in simulations of KcsA. This suggests that a single-filing mechanism is conserved between different K channel structures and may be robust to changes in simulation details. Comparison of Kir6.2 and KcsA suggests some degree of flexibility in the filter, thus complicating models of ion selectivity based upon a rigid filter.

Potassium channels enable K(+) ions to move passively across biological membranes. Multiple nanosecond-duration molecular dynamics simulations (total simulation time 5 ns) of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal motions of ions, water, and protein. Comparison of simulations with and without K(+) ions indicate that the absence of ions destabilizes the structure of the selectivity filter. Within the selectivity filter, K(+) ions interact with the backbone (carbonyl) oxygens, and with the side-chain oxygen of T75. Concerted single-file motions of water molecules and K(+) ions within the selectivity filter of the channel occur on a 100-ps time scale. In a simulation with three K(+) ions (initially two in the filter and one in the cavity), the ion within the central cavity leaves the channel via its intracellular mouth after approximately 900 ps; within the cavity this ion interacts with the Ogamma atoms of two T107 side chains, revealing a favorable site within the otherwise hydrophobically lined cavity. Exit of this ion from the channel is enabled by a transient increase in the diameter of the intracellular mouth. Such "breathing" motions may form the molecular basis of channel gating.

Isolated pore-lining helices derived from three types of K-channel have been analyzed in terms of their structural and dynamic features in nanosecond molecular dynamics (MD) simulations while spanning a lipid bilayer. The helices were 1) M1 and M2 from the bacterial channel KcsA (Streptomyces lividans), 2) S5 and S6 from the voltage-gated (Kv) channel Shaker (Drosophila melanogaster), and 3) M1 and M2 from the inward rectifier channel Kir6.2 (human). In the case of the Kv and Kir channels, for which x-ray structures are not known, both short and long models of each helix were considered. Each helix was incorporated into a lipid bilayer containing 127 palmitoyloleoylphosphatidylcholine molecules, which was solvated with approximately 4000 water molecules, yielding approximately 20, 000 atoms in each system. Nanosecond MD simulations were used to aid the definition of optimal lengths for the helix models from Kv and Kir. Thus the study corresponds to a total simulation time of 10 ns. The inner pore-lining helices (M2 in KcsA and Kir, S6 in Shaker) appear to be slightly more flexible than the outer pore-lining helices. In particular, the Pro-Val-Pro motif of S6 results in flexibility about a molecular hinge, as was suggested by previous in vacuo simulations (, Biopolymers. 39:503-515). Such flexibility may be related to gating in the corresponding intact channel protein molecules. Analysis of H-bonds revealed interactions with both water and lipid molecules in the water/bilayer interfacial region. Such H-bonding interactions may lock the helices in place in the bilayer during the folding of the channel protein (as is implicit in the two-stage model of membrane protein folding). Aromatic residues at the extremities of the helices underwent complex motions on both short (<10 ps) and long (>100 ps) time scales.