Sodium recognition by the Na+/Ca2+ exchanger in the outward-facing conformation

Balancing Group 1 Monoatomic Ion-Polar Compound Interactions in the Polarizable Drude Force Field: Application in Protein and Nucleic Acid Systems

An accurate force field (FF) is the foundation of reliable results from molecular dynamics (MD) simulations. In our recently published work, we developed a protocol to generate atom pair-specific Lennard-Jones (known as NBFIX in CHARMM) and through-space Thole dipole screening (NBTHOLE) parameters in the context of the Drude polarizable FF based on readily accessible quantum mechanical (QM) data to fit condensed phase experimental thermodynamic benchmarks, including the osmotic pressure, diffusion coefficient, ionic conductivity, and solvation free energy, when available. In the present work, the developed protocol is applied to generate NBFIX and NBTHOLE parameters for interactions between monatomic ions (specifically Li+, Na+, K+, Rb+, Cs+, and Cl-) and common functional groups found in proteins and nucleic acids. The parameters generated for each ion-functional group pair were then applied to the corresponding functional groups within proteins or nucleic acids followed by MD simulations to analyze the distribution of ions around these biomolecules. The modified FF successfully addresses the issue of overbinding observed in a previous iteration of the Drude FF. Quantitatively, the model accurately reproduces the effective charge of proteins and demonstrates a level of charge neutralization for a double-helix B-DNA in good agreement with the counterion condensation theory. Additionally, simulations involving ion competition correlate well with experimental results, following the trend Li+ > Na+ ≈ K+ > Rb+. These results validate the refined model for group 1 ion-biomolecule interactions that will facilitate the application of the polarizable Drude FF in systems in which group 1 ions play an important role.

19F-NMR Probing of Ion-Induced Conformational Changes in Detergent-Solubilized and Nanodisc-Reconstituted NCX_Mj

Consecutive interactions of 3Na+ or 1Ca2+ with the Na+/Ca2+ exchanger (NCX) result in an alternative exposure (access) of the cytosolic and extracellular vestibules to opposite sides of the membrane, where ion-induced transitions between the outward-facing (OF) and inward-facing (IF) conformational states drive a transport cycle. Here, we investigate sub-state populations of apo and ion-bound species in the OF and IF states by analyzing detergent-solubilized and nanodisc-reconstituted preparations of NCX_Mj with 19F-NMR. The 19F probe was covalently attached to the cysteine residues at entry locations of the cytosolic and extracellular vestibules. Multiple sub-states of apo and ion-bound species were observed in nanodisc-reconstituted (but not in detergent-solubilized) NCX_Mj, meaning that the lipid-membrane environment preconditions multiple sub-state populations toward the OF/IF swapping. Most importantly, ion-induced sub-state redistributions occur within each major (OF or IF) state, where sub-state interconversions may precondition the OF/IF swapping. In contrast with large changes in population redistributions, the sum of sub-state populations within each inherent state (OF or IF) remains nearly unchanged upon ion addition. The present findings allow the further elucidation of structure-dynamic modules underlying ion-induced conformational changes that determine a functional asymmetry of ion access/translocation at opposite sides of the membrane and ion transport rates concurring physiological demands.

Balancing Group I Monatomic Ion-Polar Compound Interactions for Condensed Phase Simulation in the Polarizable Drude Force Field

Molecular dynamics (MD) simulations are a commonly used method for investigating molecular behavior at the atomic level. Achieving reliable MD simulation results necessitates the use of an accurate force field. In the present work, we present a protocol to enhance the quality of group 1 monatomic ions (specifically Li+, Na+, K+, Rb+, and Cs+) with respect to their interactions with common polar model compounds in biomolecules in condensed phases in the context of the Drude polarizable force field. Instead of adjusting preexisting individual parameters for ions, model compounds, and water, we employ atom-pair specific Lennard-Jones (LJ) (known as NBFIX in CHARMM) and through-space Thole dipole screening (NBTHOLE) terms to fine-tune the balance of ion-model compound, ion-water, and model compound-water interactions. This involved establishing a protocol for the optimization of NBFIX and NBTHOLE parameters targeting the difference between molecular mechanical (MM) and quantum mechanical (QM) potential energy scans (PES). It is shown that targeting PES involving complexes that include multiple model compounds and/or ions as trimers and tetramers yields parameters that produce condensed phase properties in agreement with experimental data. Validation of this protocol involved the reproduction of experimental thermodynamic benchmarks, including solvation free energies of ions in methanol and N-methylacetamide, osmotic pressures, ionic conductivities, and diffusion coefficients within the condensed phase. These results show the importance of including more complex ion-model compound complexes beyond dimers in the QM target data to account for many-body effects during parameter fitting. The presented parameters represent a significant refinement of the Drude polarizable force field, which will lead to improved accuracy for modeling ion-biomolecular interactions.

Structural dynamics of Na+ and Ca2+ interactions with full-size mammalian NCX

Cytosolic Ca2+ and Na+ allosterically regulate Na+/Ca2+ exchanger (NCX) proteins to vary the NCX-mediated Ca2+ entry/exit rates in diverse cell types. To resolve the structure-based dynamic mechanisms underlying the ion-dependent allosteric regulation in mammalian NCXs, we analyze the apo, Ca2+, and Na+-bound species of the brain NCX1.4 variant using hydrogen-deuterium exchange mass spectrometry (HDX-MS) and molecular dynamics (MD) simulations. Ca2+ binding to the cytosolic regulatory domains (CBD1 and CBD2) rigidifies the intracellular regulatory loop (5L6) and promotes its interaction with the membrane domains. Either Na+ or Ca2+ stabilizes the intracellular portions of transmembrane helices TM3, TM4, TM9, TM10, and their connecting loops (3L4 and 9L10), thereby exposing previously unappreciated regulatory sites. Ca2+ or Na+ also rigidifies the palmitoylation domain (TMH2), and neighboring TM1/TM6 bundle, thereby uncovering a structural entity for modulating the ion transport rates. The present analysis provides new structure-dynamic clues underlying the regulatory diversity among tissue-specific NCX variants.

Conformational free-energy landscapes of a Na+/Ca2+ exchanger explain its alternating-access mechanism and functional specificity

Secondary-active transporters catalyze the movement of myriad substances across all cellular membranes, typically against opposing concentration gradients, and without consuming any ATP. To do so, these proteins employ an intriguing structural mechanism evolved to be activated only upon recognition or release of the transported species. We examine this self-regulated mechanism using a homolog of the cardiac Na+/Ca2+ exchanger as a model system. Using advanced computer simulations, we map out the complete functional cycle of this transporter, including unknown conformations that we validate against existing experimental data. Calculated free-energy landscapes reveal why this transporter functions as an antiporter rather than a symporter, why it specifically exchanges Na+ and Ca2+, and why the stoichiometry of this exchange is exactly 3:1. We also rationalize why the protein does not exchange H+ for either Ca2+ or Na+, despite being able to bind H+ and its high similarity with H+/Ca2+ exchangers. Interestingly, the nature of this transporter is not explained by its primary structural states, known as inward- and outward-open conformations; instead, the defining factor is the feasibility of conformational intermediates between those states, wherein access pathways leading to the substrate binding sites become simultaneously occluded from both sides of the membrane. This analysis offers a physically coherent, broadly transferable route to understand the emergence of function from structure among secondary-active membrane transporters.

The effect of ferroptosis-related mitochondrial dysfunction in the development of temporal lobe epilepsy

Temporal lobe epilepsy (TLE) is the most common form of epileptic syndrome. It has been established that due to its complex pathogenesis, a considerable proportion of TLE patients often progress to drug-resistant epilepsy. Ferroptosis has emerged as an important neuronal death mechanism in TLE, which is primarily influenced by lipid accumulation and oxidative stress. In previous studies of ferroptosis, more attention has been focused on the impact of changes in the levels of proteins related to the redox equilibrium and signaling pathways on epileptic seizures. However, it is worth noting that the oxidative-reduction changes in different organelles may have different pathophysiological significance in the process of ferroptosis-related diseases. Mitochondria, as a key organelle involved in ferroptosis, its structural damage and functional impairment can lead to energy metabolism disorders and disruption of the excitatory inhibitory balance, significantly increasing the susceptibility to epileptic seizures. Therefore, secondary mitochondrial dysfunction in the process of ferroptosis could play a crucial role in TLE pathogenesis. This review focuses on ferroptosis and mitochondria, discussing the pathogenic role of ferroptosis-related mitochondrial dysfunction in TLE, thus aiming to provide novel insights and potential implications of ferroptosis-related secondary mitochondrial dysfunction in epileptic seizures and to offer new insights for the precise exploration of ferroptosis-related therapeutic targets for TLE patients.

Neuronal and astrocyte NCX isoform/splice variants: How do they participate in Na+ and Ca2+ signalling?

NCX1, NCX2, and NCX3 gene isoforms and their splice variants are characteristically expressed in different regions of the brain. The tissue-specific splice variants of NCX1-3 isoforms show specific expression profiles in neurons and astrocytes, whereas the relevant NCX isoform/splice variants exhibit diverse allosteric modes of Na+- and Ca2+-dependent regulation. In general, overexpression of NCX1-3 genes leads to neuroprotective effects, whereas their ablation gains the opposite results. At this end, the partial contributions of NCX isoform/splice variants to neuroprotective effects remain unresolved. The glutamate-dependent Na+ entry generates Na+ transients (in response to neuronal cell activities), whereas the Na+-driven Ca2+ entry (through the reverse NCX mode) raises Ca2+ transients. This special mode of signal coupling translates Na+ transients into the Ca2+ signals while being a part of synaptic neurotransmission. This mechanism is of general interest since disease-related conditions (ischemia, metabolic stress, and stroke among many others) trigger Na+ and Ca2+ overload with deadly outcomes of downstream apoptosis and excitotoxicity. The recently discovered mechanisms of NCX allosteric regulation indicate that some NCX variants might play a critical role in the dynamic coupling of Na+-driven Ca2+ entry. In contrast, the others are less important or even could be dangerous under altered conditions (e.g., metabolic stress). This working hypothesis can be tested by applying advanced experimental approaches and highly focused computational simulations. This may allow the development of structure-based blockers/activators that can selectively modulate predefined NCX variants to lessen the life-threatening outcomes of excitotoxicity, ischemia, apoptosis, metabolic deprivation, brain injury, and stroke.

Ion channel selectivity through ion-modulated changes of selectivity filter pKa values

In bacterial voltage-gated sodium channels, the passage of ions through the pore is controlled by a selectivity filter (SF) composed of four glutamate residues. The mechanism of selectivity has been the subject of intense research, with suggested mechanisms based on steric effects, and ion-triggered conformational change. Here, we propose an alternative mechanism based on ion-triggered shifts in pKa values of SF glutamates. We study the NavMs channel for which the open channel structure is available. Our free-energy calculations based on molecular dynamics simulations suggest that pKa values of the four glutamates are higher in solution of K+ ions than in solution of Na+ ions. Higher pKa in the presence of K+ stems primarily from the higher population of dunked conformations of the protonated Glu sidechain, which exhibit a higher pKa shift. Since pKa values are close to the physiological pH, this results in predominant population of the fully deprotonated state of glutamates in Na+ solution, while protonated states are predominantly populated in K+ solution. Through molecular dynamics simulations we calculate that the deprotonated state is the most conductive, the singly protonated state is less conductive, and the doubly protonated state has significantly reduced conductance. Thus, we propose that a significant component of selectivity is achieved through ion-triggered shifts in the protonation state, which favors more conductive states for Na+ ions and less conductive states for K+ ions. This mechanism also suggests a strong pH dependence of selectivity, which has been experimentally observed in structurally similar NaChBac channels.

Ca2+ efflux facilitated by co-transport of inorganic phosphate anion in the H+/Ca2+ antiporter YfkE

Ca²⁺ is an important signaling messenger. In microorganisms, fungi, and plants, H⁺/Ca²⁺ antiporters (CAX) are known to play key roles in the homeostasis of intracellular Ca²⁺ by catalyzing its efflux across the cell membrane. Here, we reveal that the bacterial CAX homolog YfkE transports Ca²⁺ in two distinct modes: a low-flux H⁺/Ca²⁺ exchange mode and a high-flux mode in which Ca²⁺ and phosphate ions are co-transported (1:1) in exchange for H⁺. Coupling with phosphate greatly accelerates the Ca²⁺ efflux activity of YfkE. Our studies reveal that Ca²⁺ and phosphate bind to adjacent sites in a central translocation pathway and lead to mechanistic insights that explain how this CAX alters its conserved alpha-repeat motifs to adopt phosphate as a specific "transport chaperon" for Ca²⁺ translocation. This finding uncovers a co-transport mechanism within the CAX family that indicates this class of proteins contributes to the cellular homeostasis of both Ca²⁺ and phosphate.

The role of the Na+/Ca2+-exchanger (NCX) in cancer-associated fibroblasts

Cancer is characterized by uncontrolled growth, invasion, and metastasis. In addition to solid cancer cells, cancer-associated fibroblasts (CAFs) play important roles in cancer pathophysiology. They arise from "healthy" cells but get manipulated by solid cancer cells to supply them and develop a tumor microenvironment (TME) that protects the cancer cells from the immune defense. A wide variety of cell types can differentiate into CAFs, including fibroblasts, endothelial cells, and epithelial cells. Precise Ca²⁺ regulation is essential for each cell including CAFs. The electrogenic Na⁺/Ca²⁺ exchanger (NCX) is one of the ubiquitously expressed regulatory Ca²⁺ transport proteins that rapidly responds to changes of the intracellular ion concentrations. Its transport function is also influenced by the membrane potential and thereby indirectly by the activity of ion channels. NCX transports Ca²⁺ out of the cell (forward mode) or allows its influx (reverse mode), always in exchange for 3 Na⁺ ions that are moved into the opposite direction. In this review, we discuss the functional roles NCX has in CAFs and how these depend on the properties of the TME. NCX activity modifies migration and leads to a reduced proliferation and apoptosis. The effect of the NCX in fibrosis is still largely unknown.

Structure-Based Function and Regulation of NCX Variants: Updates and Challenges

The plasma-membrane homeostasis Na⁺/Ca²⁺ exchangers (NCXs) mediate Ca²⁺ extrusion/entry to dynamically shape Ca²⁺ signaling/in biological systems ranging from bacteria to humans. The NCX gene orthologs, isoforms, and their splice variants are expressed in a tissue-specific manner and exhibit nearly 10⁴-fold differences in the transport rates and regulatory specificities to match the cell-specific requirements. Selective pharmacological targeting of NCX variants could benefit many clinical applications, although this intervention remains challenging, mainly because a full-size structure of eukaryotic NCX is unavailable. The crystal structure of the archaeal NCX_Mj, in conjunction with biophysical, computational, and functional analyses, provided a breakthrough in resolving the ion transport mechanisms. However, NCX_Mj (whose size is nearly three times smaller than that of mammalian NCXs) cannot serve as a structure-dynamic model for imitating high transport rates and regulatory modules possessed by eukaryotic NCXs. The crystal structures of isolated regulatory domains (obtained from eukaryotic NCXs) and their biophysical analyses by SAXS, NMR, FRET, and HDX-MS approaches revealed structure-based variances of regulatory modules. Despite these achievements, it remains unclear how multi-domain interactions can decode and integrate diverse allosteric signals, thereby yielding distinct regulatory outcomes in a given ortholog/isoform/splice variant. This article summarizes the relevant issues from the perspective of future developments.

Exploring the Li+ transporting mutant of NCX_Mj for assigning ion binding sites of mitochondrial NCLX

The plasma membrane (NCX) and mitochondrial (NCLX) Na⁺/Ca²⁺ exchangers are structurally related proteins, although they operate under strictly different ionic conditions and membrane potentials. In contrast with NCX, NCLX can transport either Li⁺ or Na⁺ in exchange for Ca²⁺. Whereas the crystal structure of the archaeal NCX (NCX_Mj) describes the binding sites for alternative binding of 3Na⁺ or 1Ca²⁺, these features remain elusive for NCLX due to the lack of structural information. To elucidate the ion-binding features of mitochondrial NCLX, we analyzed here the Li⁺-transporting NCLX_Mj mutant, produced by replacing the ion-coordinating residues in the archaeal NCX (NCX_Mj) to match the ion-coordinating residues of human NCLX. The NCLX_Mj-mediated Na⁺/Ca²⁺ or Li⁺/Ca²⁺ exchange rates are insensitive to varying voltage, consistent with an electroneutral ion exchange. Molecular dynamics (MD) simulations revealed that NCLX_Mj contains two novel Li⁺ binding sites with four ion-coordinating residues, derived from the three Na⁺ binding sites of NCX_Mj. The ion-coordination modes, observed in the MD simulations, were further supported by two-dimensional infrared (2D IR) spectroscopy and by testing the mutational effects on the ion fluxes. Collectively, our results revealed a structural basis for Li⁺ binding and electroneutral transport (2Na⁺/Li⁺:1Ca²⁺) by NCLX_Mj, meaning that the NCLX-mediated electroneutral transport may predefine mitochondrial Ca²⁺ and Na⁺ signaling to modulate cellular functions.

Mitochondrial calcium exchange in physiology and disease

The uptake of calcium into and extrusion of calcium from the mitochondrial matrix is a fundamental biological process that has critical effects on cellular metabolism, signaling, and survival. Disruption of mitochondrial calcium (mCa2+) cycling is implicated in numerous acquired diseases such as heart failure, stroke, neurodegeneration, diabetes, and cancer and is genetically linked to several inherited neuromuscular disorders. Understanding the mechanisms responsible for mCa2+ exchange therefore holds great promise for the treatment of these diseases. The past decade has seen the genetic identification of many of the key proteins that mediate mitochondrial calcium uptake and efflux. Here, we present an overview of the phenomenon of mCa2+ transport and a comprehensive examination of the molecular machinery that mediates calcium flux across the inner mitochondrial membrane: the mitochondrial uniporter complex (consisting of MCU, EMRE, MICU1, MICU2, MICU3, MCUB, and MCUR1), NCLX, LETM1, the mitochondrial ryanodine receptor, and the mitochondrial permeability transition pore. We then consider the physiological implications of mCa2+ flux and evaluate how alterations in mCa2+ homeostasis contribute to human disease. This review concludes by highlighting opportunities and challenges for therapeutic intervention in pathologies characterized by aberrant mCa2+ handling and by summarizing critical unanswered questions regarding the biology of mCa2+ flux.

The Cardiac Na+ -Ca2+ Exchanger: From Structure to Function

Ca²⁺ homeostasis is essential for cell function and survival. As such, the cytosolic Ca²⁺ concentration is tightly controlled by a wide number of specialized Ca²⁺ handling proteins. One among them is the Na⁺ -Ca²⁺ exchanger (NCX), a ubiquitous plasma membrane transporter that exploits the electrochemical gradient of Na⁺ to drive Ca²⁺ out of the cell, against its concentration gradient. In this critical role, this secondary transporter guides vital physiological processes such as Ca²⁺ homeostasis, muscle contraction, bone formation, and memory to name a few. Herein, we review the progress made in recent years about the structure of the mammalian NCX and how it relates to function. Particular emphasis will be given to the mammalian cardiac isoform, NCX1.1, due to the extensive studies conducted on this protein. Given the degree of conservation among the eukaryotic exchangers, the information highlighted herein will provide a foundation for our understanding of this transporter family. We will discuss gene structure, alternative splicing, topology, regulatory mechanisms, and NCX's functional role on cardiac physiology. Throughout this article, we will attempt to highlight important milestones in the field and controversial topics where future studies are required. © 2021 American Physiological Society. Compr Physiol 12:1-37, 2021.

Corrections in the CHARMM36 Parametrization of Chloride Interactions with Proteins, Lipids, and Alkali Cations, and Extension to Other Halide Anions

The nonpolarizable CHARMM force field is one of the most widely used energy functions for all-atom biomolecular simulations. Chloride is the only halide ion included in the latest version, CHARMM36m, and is used widely in simulation studies, often as an electrolyte ion but also as the biological substrate of transport proteins and enzymes. Here, we find that existing parameters systematically underestimate the interaction of Cl^- with proteins and lipids. Accordingly, when examined in solution, little to no Cl^-association can be observed with most components of the protein, including backbone, polar side chains and aromatic rings. The strength of the interaction with cationic side chains and with alkali ions is also incongruent with experimental measurements, specifically osmotic coefficients of concentrated solutions. Consistent with these findings, a 4-μs trajectory of the Cl^--specific transport protein CLC-ec1 shows irreversible Cl^- dissociation from the so-called S_cen binding site, even in a 150 mM NaCl buffer. To correct for these deficiencies, we formulate a series of pair-specific Lennard-Jones parameters that override those resulting from the conventional Lorentz-Berthelot combination rules. These parameters, referred to as NBFIX, are systematically calibrated against available experimental data as well as ab initio geometry optimizations and energy evaluations, for a wide set of binary and ternary Cl^- complexes with protein and lipid analogs and alkali cations. Analogously, we also formulate parameter sets for the other three biological halide ions, namely, fluoride, bromide, and iodide. The resulting parameters are used to calculate the potential of mean force defining the interaction of each anion and each of the protein and lipid analogues in bulk water, revealing association free energies in the range of -0.3 to -3.3 kcal/mol, with the F^- complexes being the least stable. The NBFIX corrections also preserve the Cl^- occupancy of CLC-ec1 in a second 4-μs trajectory. We posit that these optimized molecular-mechanics models provide a more realistic foundation for all-atom simulation studies of processes entailing changes in hydration, recognition, or transport of halide anions.

The Archaeal Na+/Ca2+ Exchanger (NCX_Mj) as a Model of Ion Transport for the Superfamily of Ca2+/CA Antiporters

The superfamily of Calcium/Cation (Ca²⁺/CA) antiporters extrude Ca²⁺ from the cytosol or subcellular compartments in exchange with Na⁺, K⁺, H⁺, Li⁺, or Mg²⁺ and thereby provide a key mechanism for Ca²⁺ signaling and ion homeostasis in biological systems ranging from bacteria to humans. The structure-dynamic determinants of ion selectivity and transport rates remain unclear, although this is of primary physiological significance. Despite wide variances in the ion selectivity and transport rates, the Ca²⁺/CA proteins share structural motifs, although it remains unclear how the ion recognition/binding is coupled to the ion translocation events. Here, the archaeal Na⁺/Ca²⁺ exchanger (NCX_Mj) is considered as a structure-based model that can help to resolve the ion transport mechanisms by using X-ray, HDX-MS, ATR-FTIR, and computational approaches in conjunction with functional analyses of mutants. Accumulating data reveal that the local backbone dynamics at ion-coordinating residues is characteristically constrained in apo NCX_Mj, which may predefine the affinity and stability of ion-bound species in the ground and transition states. The 3Na⁺ or 1Ca²⁺ binding to respective sites of NCX_Mj rigidify the backbone dynamics at specific segments, where the ion-dependent compression of the ion-permeating four-helix bundle (TM2, TM3, TM7, and TM8) induces the sliding of the two-helix cluster (TM1/TM6) on the protein surface to switch the OF (outward-facing) and IF (inward-facing) conformations. Taking into account the common structural elements shared by Ca²⁺/CAs, NCX_Mj may serve as a model for studying the structure-dynamic and functional determinants of ion-coupled alternating access, transport catalysis, and ion selectivity in Ca²⁺/CA proteins.

Biallelic ADGRV1 variants are associated with Rolandic epilepsy

OBJECTIVE

Rolandic epilepsy (RE) is among the most common focal epilepsies in childhood. For the majority of patients with RE and atypical RE (ARE), the etiology remains elusive. We thus screened patients with RE/ARE in order to detect disease-causing variants..

METHODS

A trios-based whole-exome sequencing approach was performed in a cohort of 28 patients with RE/ARE. Clinical data and EEGs were reviewed. Variants were validated by Sanger sequencing.

RESULTS

Two compound heterozygous missense variants p.Val272Ile/p.Asn3028Ser and p.Ala3657Val/p.Met4419Val of ADGRV1 were identified in two unrelated familial cases of RE/ARE. All the variants were in the calcium exchanger β domain and were suggested to be damaging by at least one web-based prediction tool. These variants are not present or are present at a very low minor allele frequency in the gnomAD database. Previously, biallelic ADGRV1 variants (p.Gly2756Arg and p.Glu4410Lys) have been observed in RE, consistent with the observation in this study and supporting the association between ADGRV1 variants and RE. Additionally, a de novo mutation, p.Asp668Asn, in GRIN2B was identified in a sporadic case of ARE, and a missense variant, p.Asn1551Ser, in RyR2 was identified in a family with RE with incomplete penetrance. These genes are all calcium homeostasis associated genes, suggesting the potential effect of calcium homeostasis in RE/ARE.

CONCLUSIONS

The results from the present study suggest that the genes ADGRV1, GRIN2B, and RyR2 are associated with RE/ARE. These data link defects in neuronal intracellular calcium homeostasis to RE/ARE pathogenesis implicating that these defects plays an important role in the development of these conditions.

Determinants of ion selectivity in ASIC1a- and ASIC2a-containing acid-sensing ion channels

Trimeric acid-sensing ion channels (ASICs) contribute to neuronal signaling by converting extracellular acidification into excitatory sodium currents. Previous work with homomeric ASIC1a implicates conserved leucine (L7') and consecutive glycine-alanine-serine (GAS belt) residues near the middle, and conserved negatively charged (E18') residues at the bottom of the pore in ion permeation and/or selectivity. However, a conserved mechanism of ion selectivity throughout the ASIC family has not been established. We therefore explored the molecular determinants of ion selectivity in heteromeric ASIC1a/ASIC2a and homomeric ASIC2a channels using site-directed mutagenesis, electrophysiology, and molecular dynamics free energy simulations. Similar to ASIC1a, E18' residues create an energetic preference for sodium ions at the lower end of the pore in ASIC2a-containing channels. However, and in contrast to ASIC1a homomers, ion permeation through ASIC2a-containing channels is not determined by L7' side chains in the upper part of the channel. This may be, in part, due to ASIC2a-specific negatively charged residues (E59 and E62) that lower the energy of ions in the upper pore, thus making the GAS belt more important for selectivity. This is confirmed by experiments showing that the L7'A mutation has no effect in ASIC2a, in contrast to ASIC1a, where it eliminated selectivity. ASIC2a triple mutants eliminating both L7' and upper charges did not lead to large changes in selectivity, suggesting a different role for L7' in ASIC2a compared with ASIC1a channels. In contrast, we observed measurable changes in ion selectivity in ASIC2a-containing channels with GAS belt mutations. Our results suggest that ion conduction and selectivity in the upper part of the ASIC pore may differ between subtypes, whereas the essential role of E18' in ion selectivity is conserved. Furthermore, we demonstrate that heteromeric channels containing mutations in only one of two ASIC subtypes provide a means of functionally testing mutations that render homomeric channels nonfunctional.

On the Role of a Conserved Methionine in the Na+-Coupling Mechanism of a Neurotransmitter Transporter Homolog

Excitatory amino acid transporters (EAAT) play a key role in glutamatergic synaptic communication. Driven by transmembrane cation gradients, these transporters catalyze the reuptake of glutamate from the synaptic cleft once this neurotransmitter has been utilized for signaling. Two decades ago, pioneering studies in the Kanner lab identified a conserved methionine within the transmembrane domain as key for substrate turnover rate and specificity; later structural work, particularly for the prokaryotic homologs Glt_Ph and Glt_Tk, revealed that this methionine is involved in the coordination of one of the three Na⁺ ions that are co-transported with the substrate. Albeit extremely atypical, the existence of this interaction is consistent with biophysical analyses of Glt_Ph showing that mutations of this methionine diminish the binding cooperativity between substrates and Na⁺. It has been unclear, however, whether this intriguing methionine influences the thermodynamics of the transport reaction, i.e., its substrate:ion stoichiometry, or whether it simply fosters a specific kinetics in the binding reaction, which, while influential for the turnover rate, do not fundamentally explain the ion-coupling mechanism of this class of transporters. Here, studies of Glt_Tk using experimental and computational methods independently arrive at the conclusion that the latter hypothesis is the most plausible, and lay the groundwork for future efforts to uncover the underlying mechanism.

Conserved binding site in the N-lobe of prokaryotic MATE transporters suggests a role for Na+ in ion-coupled drug efflux

In both prokaryotes and eukaryotes, multidrug and toxic-compound extrusion (MATE) transporters catalyze the efflux of a broad range of cytotoxic compounds, including human-made antibiotics and anticancer drugs. MATEs are secondary-active antiporters, i.e., their drug-efflux activity is coupled to, and powered by, the uptake of ions down a preexisting transmembrane electrochemical gradient. Key aspects of this mechanism, however, remain to be delineated, such as its ion specificity and stoichiometry. We previously revealed the existence of a Na⁺-binding site in a MATE transporter from Pyroccocus furiosus (PfMATE) and hypothesized that this site might be broadly conserved among prokaryotic MATEs. Here, we evaluate this hypothesis by analyzing VcmN and ClbM, which along with PfMATE are the only three prokaryotic MATEs whose molecular structures have been determined at atomic resolution, i.e. better than 3 Å. Reinterpretation of existing crystallographic data and molecular dynamics simulations indeed reveal an occupied Na⁺-binding site in the N-terminal lobe of both structures, analogous to that identified in PfMATE. We likewise find this site to be strongly selective against K⁺, suggesting it is mechanistically significant. Consistent with these computational results, DEER spectroscopy measurements for multiple doubly-spin-labeled VcmN constructs demonstrate Na⁺-dependent changes in protein conformation. The existence of this binding site in three MATE orthologs implicates Na⁺ in the ion-coupled drug-efflux mechanisms of this class of transporters. These results also imply that observations of H⁺-dependent activity likely stem either from a site elsewhere in the structure, or from H⁺ displacing Na⁺ under certain laboratory conditions, as has been noted for other Na⁺-driven transport systems.

An embedded lipid in the multidrug transporter LmrP suggests a mechanism for polyspecificity

Multidrug efflux pumps present a challenge to the treatment of bacterial infections, making it vitally important to understand their mechanism of action. Here, we investigate the nature of substrate binding within Lactococcus lactis LmrP, a prototypical multidrug transporter of the major facilitator superfamily. We determined the crystal structure of LmrP in a ligand-bound outward-open state and observed an embedded lipid in the binding cavity of LmrP, an observation supported by native mass spectrometry analyses. Molecular dynamics simulations suggest that the anionic lipid stabilizes the observed ligand-bound structure. Mutants engineered to disrupt binding of the embedded lipid display reduced transport of some, but not all, antibiotic substrates. Our results suggest that a lipid within the binding cavity could provide a malleable hydrophobic component that allows adaptation to the presence of different substrates, helping to explain the broad specificity of this protein and possibly other multidrug transporters.

New Molecular-Mechanics Model for Simulations of Hydrogen Fluoride in Chemistry and Biology

Hydrogen fluoride (HF) is the most polar diatomic molecule and one of the simplest molecules capable of hydrogen-bonding. HF deviates from ideality both in the gas phase and in solution and is thus of great interest from a fundamental standpoint. Pure and aqueous HF solutions are broadly used in chemical and industrial processes, despite their high toxicity. HF is a stable species also in some biological conditions, because it does not readily dissociate in water unlike other hydrogen halides; yet, little is known about how HF interacts with biomolecules. Here, we set out to develop a molecular-mechanics model to enable computer simulations of HF in chemical and biological applications. This model is based on a comprehensive high-level ab initio quantum chemical investigation of the structure and energetics of the HF monomer and dimer; (HF)_n clusters, for n = 3-7; various clusters of HF and H₂O; and complexes of HF with analogs of all 20 amino acids and of several commonly occurring lipids, both neutral and ionized. This systematic analysis explains the unique properties of this molecule: for example, that interacting HF molecules favor nonlinear geometries despite being diatomic and that HF is a strong H-bond donor but a poor acceptor. The ab initio data also enables us to calibrate a three-site molecular-mechanics model, with which we investigate the structure and thermodynamic properties of gaseous, liquid, and supercritical HF in a wide range of temperatures and pressures; the solvation structure of HF in water and of H₂O in liquid HF; and the free diffusion of HF across a lipid bilayer, a key process underlying the high cytotoxicity of HF. Despite its inherent simplifications, the model presented significantly improves upon previous efforts to capture the properties of pure and aqueous HF fluids by molecular-mechanics methods and to our knowledge constitutes the first parameter set calibrated for biomolecular simulations.

A glimpse into the molecular mechanism of integral membrane proteins through hydrogen-deuterium exchange mass spectrometry

Integral membrane proteins (IMPs) control countless fundamental biological processes and constitute the majority of drug targets. For this reason, uncovering their molecular mechanism of action has long been an intense field of research. They are, however, notoriously difficult to work with, mainly due to their localization within the heterogeneous of environment of the biological membrane and the instability once extracted from the lipid bilayer. High‐resolution structures have unveiled many mechanistic aspects of IMPs but also revealed that the elucidation of static pictures has limitations. Hydrogen–deuterium exchange coupled to mass spectrometry (HDX‐MS) has recently emerged as a powerful biophysical tool for interrogating the conformational dynamics of proteins and their interactions with ligands. Its versatility has proven particularly useful to reveal mechanistic aspects of challenging classes of proteins such as IMPs. This review recapitulates the accomplishments of HDX‐MS as it has matured into an essential tool for membrane protein structural biologists.

Structure-affinity insights into the Na+ and Ca2+ interactions with multiple sites of a sodium-calcium exchanger

Selective recognition and transport of Na⁺ and Ca²⁺ ions by sodium-calcium exchanger (NCX) proteins is a primary prerequisite for Ca²⁺ signaling and homeostasis. Twelve ion-coordinating residues are highly conserved among NCXs, and distinct NCX orthologs contain two or three carboxylates, while sharing a common ion-exchange stoichiometry (3Na⁺ :1Ca²⁺ ). How these structural differences affect the ion-binding affinity, selectivity, and transport rates remains unclear. Here, the mutational effects of three carboxylates (E54, E213, and D240) were analyzed on the ion-exchange rates in the archaeal NCX from Methanococcus jannaschii and ion-induced structure-affinity changes were monitored by attenuated total reflection-Fourier-transform infrared spectroscopy (ATR-FTIR). The D240N mutation elevated the ion-transport rates by twofold to threefold, meaning that the deprotonation of D240 is not essential for transport catalysis. In contrast, mutating E54 or E213 to A, D, N, or Q dramatically decreased the ion-transport rates. ATR-FTIR revealed high- and low-affinity binding of Na⁺ or Ca²⁺ with E54 and E213, but not with D240. These findings reveal distinct structure-affinity states at specific ion-binding sites in the inward-facing (IF) and outward-facing orientation. Collectively, two multidentate carboxylate counterparts (E54 and E213) play a critical role in determining the ion coordination/transport in prokaryotic and eukaryotic NCXs, whereas the ortholog substitutions in prokaryotes (aspartate) and eukaryotes (asparagine) at the 240 position affect the ion-transport rates differently (k_cat ), probably due to the structural differences in the transition state.

Computational Dissection of Membrane Transport at a Microscopic Level

Membrane transporters are key gatekeeper proteins at cellular membranes that closely control the traffic of materials. Their function relies on structural rearrangements of varying degrees that facilitate substrate translocation across the membrane. Characterizing these functionally important molecular events at a microscopic level is key to our understanding of membrane transport, yet challenging to achieve experimentally. Recent advances in simulation technology and computing power have rendered molecular dynamics (MD) simulation a powerful biophysical tool to investigate a wide range of dynamical events spanning multiple spatial and temporal scales. Here, we review recent studies of diverse membrane transporters using computational methods, with an emphasis on highlighting the technical challenges, key lessons learned, and new opportunities to illuminate transporter structure and function.

Hydrogen-Deuterium Exchange Mass-Spectrometry of Secondary Active Transporters: From Structural Dynamics to Molecular Mechanisms

Membrane transporters allow the selective transport of otherwise poorly permeable solutes across the cell membrane and thus, play a key role in maintaining cellular homeostasis in all kingdoms of life. Importantly, these proteins also serve as important drug targets. Over the last decades, major progress in structural biology methods has elucidated important structure-function relationships in membrane transporters. However, structures obtained using methods such as X-ray crystallography and high-resolution cryogenic electron microscopy merely provide static snapshots of an intrinsically dynamic, multi-step transport process. Therefore, there is a growing need for developing new experimental approaches capable of exploiting the data obtained from the high-resolution snapshots in order to investigate the dynamic features of membrane proteins. Here, we present the basic principles of hydrogen-deuterium exchange mass-spectrometry (HDX-MS) and recent advancements in its use to study membrane transporters. In HDX-MS experiments, minute amounts of a protein sample can be used to investigate its structural dynamics under native conditions, without the need for chemical labelling and with practically no limit on the protein size. Thus, HDX-MS is instrumental for resolving the structure-dynamic landscapes of membrane proteins in their apo (ligand-free) and ligand-bound forms, shedding light on the molecular mechanism underlying the transport process and drug binding.

Multipurpose Na+ ions mediate excitation and cellular homeostasis: Evolution of the concept of Na+ pumps and Na+/Ca2+ exchangers

Ionic signalling is the most ancient form of regulation of cellular functions in response to environmental challenges. Signals, mediated by Na⁺ fluxes and spatio-temporal fluctuations of Na⁺ concentration in cellular organelles and cellular compartments contribute to the most fundamental cellular processes such as membrane excitability and energy production. At the very core of ionic signalling lies the Na⁺-K⁺ ATP-driven pump (or NKA) which creates trans-plasmalemmal ion gradients that sustain ionic fluxes through ion channels and numerous Na⁺-dependent transporters that maintain cellular and tissue homeostasis. Here we present a brief account of the history of research into NKA, Na⁺ -dependent transporters and Na⁺ signalling.

Calcium Signaling and Reactive Oxygen Species in Mitochondria

In heart failure, alterations of Na⁺ and Ca²⁺ handling, energetic deficit, and oxidative stress in cardiac myocytes are important pathophysiological hallmarks. Mitochondria are central to these processes because they are the main source for ATP, but also reactive oxygen species (ROS), and their function is critically controlled by Ca²⁺ During physiological variations of workload, mitochondrial Ca²⁺ uptake is required to match energy supply to demand but also to keep the antioxidative capacity in a reduced state to prevent excessive emission of ROS. Mitochondria take up Ca²⁺ via the mitochondrial Ca²⁺ uniporter, which exists in a multiprotein complex whose molecular components were identified only recently. In heart failure, deterioration of cytosolic Ca²⁺ and Na⁺ handling hampers mitochondrial Ca²⁺ uptake and the ensuing Krebs cycle-induced regeneration of the reduced forms of NADH (nicotinamide adenine dinucleotide) and NADPH (nicotinamide adenine dinucleotide phosphate), giving rise to energetic deficit and oxidative stress. ROS emission from mitochondria can trigger further ROS release from neighboring mitochondria termed ROS-induced ROS release, and cross talk between different ROS sources provides a spatially confined cellular network of redox signaling. Although low levels of ROS may serve physiological roles, higher levels interfere with excitation-contraction coupling, induce maladaptive cardiac remodeling through redox-sensitive kinases, and cell death through mitochondrial permeability transition. Targeting the dysregulated interplay between excitation-contraction coupling and mitochondrial energetics may ameliorate the progression of heart failure.

A K+/Na+ co-binding state: Simultaneous versus competitive binding of K+ and Na+ to glutamate transporters

Plasma membrane-associated glutamate transporters play a key role in signaling by the major excitatory neurotransmitter glutamate. Uphill glutamate uptake into cells is energetically driven by coupling to co-transport of three Na⁺ ions. In exchange, one K⁺ ion is counter-transported. Currently accepted transport mechanisms assume that Na⁺ and K⁺ effects are exclusive, resulting from competition of these cations at the binding level. Here, we used electrophysiological analysis to test the effects of K⁺ and Na⁺ on neuronal glutamate transporter excitatory amino acid carrier 1 (EAAC1; the rat homologue of human excitatory amino acid transporter 3 (EAAT3)). Unexpectedly, extracellular K⁺ application to EAAC1 induced anion current, but only in the presence of Na⁺ This result could be explained with a K⁺/Na⁺ co-binding state in which the two cations simultaneously bind to the transporter. We obtained further evidence for this co-binding state, and its anion conductance, by analyzing transient currents when Na⁺ was exchanged for K⁺ and effects of the [K⁺]/[Na⁺] ratio on glutamate affinity. Interestingly, we observed the K⁺/Na⁺ co-binding state not only in EAAC1 but also in the subtypes EAAT1 and -2, which, unlike EAAC1, conducted anions in response to K⁺ only. We incorporated these experimental findings in a revised transport mechanism, including the K⁺/Na⁺ co-binding state and the ability of K⁺ to activate anion current. Overall, these results suggest that differentiation between Na⁺ and K⁺ does not occur at the binding level but is conferred by coupling of cation binding to conformational changes. These findings have implications also for other exchangers.

Mitochondrial Ca2+ Transport: Mechanisms, Molecular Structures, and Role in Cells

Mitochondria are among the most important cell organelles involved in the regulation of intracellular calcium homeostasis. During the last decade, a number of molecular structures responsible for the mitochondrial calcium transport have been identified including the mitochondrial Ca2+ uniporter (MCU), Na+/Ca2+ exchanger (NCLX), and Ca2+/H+ antiporter (Letm1). The review summarizes the data on the structure, regulation, and physiological role of such structures. The pathophysiological mechanism of Ca2+ transport through the cyclosporine A-sensitive mitochondrial permeability transition pore is discussed. An alternative mechanism for the mitochondrial pore opening, namely, formation of the lipid pore induced by saturated fatty acids, and its role in Ca2+ transport are described in detail.

Key residues controlling bidirectional ion movements in Na+/Ca2+ exchanger

Prokaryotic and eukaryotic Na+/Ca2+ exchangers (NCX) control Ca2+ homeostasis. NCX orthologs exhibit up to 104-fold differences in their turnover rates (kcat), whereas the ratios between the cytosolic (cyt) and extracellular (ext) Km values (Kint = KmCyt/KmExt) are highly asymmetric and alike (Kint ≤ 0.1) among NCXs. The structural determinants controlling a huge divergence in kcat at comparable Kint remain unclear, although 11 (out of 12) ion-coordinating residues are highly conserved among NCXs. The crystal structure of the archaeal NCX (NCX_Mj) was explored for testing the mutational effects of pore-allied and loop residues on kcat and Kint. Among 55 tested residues, 26 mutations affect either kcat or Kint, where two major groups can be distinguished. The first group of mutations (14 residues) affect kcat rather than Kint. The majority of these residues (10 out of 14) are located within the extracellular vestibule near the pore center. The second group of mutations (12 residues) affect Kint rather than kcat, whereas the majority of residues (9 out 12) are randomly dispersed within the extracellular vestibule. In conjunction with computational modeling-simulations and hydrogen-deuterium exchange mass-spectrometry (HDX-MS), the present mutational analysis highlights structural elements that differentially govern the intrinsic asymmetry and transport rates. The key residues, located at specific segments, can affect the characteristic features of local backbone dynamics and thus, the conformational flexibility of ion-transporting helices contributing to critical conformational transitions. The underlying mechanisms might have a physiological relevance for matching the response modes of NCX variants to cell-specific Ca2+ and Na+ signaling.

Molecular Dynamics Simulations of Membrane Proteins: An Overview

Simulations of membrane proteins have been rising in popularity in the past decade. Advancements in technology and force fields made it possible to simulate behavior of membrane proteins. Membrane protein simulations can now be used as supporting evidence for experimental findings, for elucidating protein mechanisms, and validating protein crystal structures. Unrelated to experimental data, these simulations can also serve to investigate larger scale processes like protein sorting, protein-membrane interactions, and more. In this review, the history as well as the state-of-the-art methodologies in membrane protein simulations will be summarized. An emphasis will be put on how to set up the system and on the current models for the different components of the simulation system. An overview of the available tools for membrane protein simulation will be given, and current limitations and prospects will also be discussed.

Broadly conserved Na+-binding site in the N-lobe of prokaryotic multidrug MATE transporters

Multidrug and toxic-compound extrusion (MATE) proteins comprise an important but largely uncharacterized family of secondary-active transporters. In both eukaryotes and prokaryotes, these transporters protect the cell by catalyzing the efflux of a broad range of cytotoxic compounds, including human-made antibiotics and anticancer drugs. MATEs are thus potential pharmacological targets against drug-resistant pathogenic bacteria and tumor cells. The activity of MATEs is powered by transmembrane electrochemical ion gradients, but their molecular mechanism and ion specificity are not understood, in part because high-quality structural information is limited. Here, we use computational methods to study PfMATE, from Pyrococcus furiosus, whose structure is the best resolved to date. Analysis of available crystallographic data and additional molecular dynamics simulations unequivocally reveal an occupied Na⁺-binding site in the N-lobe of this transporter, which had not been previously recognized. We find this site to be selective against K⁺ and broadly conserved among prokaryotic MATEs, including homologs known to be Na⁺-dependent such as NorM-VC, VmrA, and ClbM, for which the location of the Na⁺ site had been debated. We note, however, that the chemical makeup of the proposed Na⁺ site indicates it is weakly specific against H⁺, explaining why MATEs featuring this Na⁺-binding motif may be solely driven by H⁺ in laboratory conditions. We further posit that the concurrent coupling to H⁺ and Na⁺ gradients observed for some Na⁺-driven MATEs owes to a second H⁺-binding site, within the C-lobe. In summary, our study provides insights into the structural basis for the complex ion dependency of MATE transporters.

The prokaryotic Na+/Ca2+ exchanger NCX_Mj transports Na+ and Ca2+ in a 3:1 stoichiometry

Sodium–calcium exchangers contribute to the generation of intracellular Ca²⁺ signals in numerous physiological processes. Shlosman et al. determine the ion stoichiometry of the only sodium–calcium exchanger of known atomic structure, revealing its functional similarity to mammalian exchangers.

Intracellular Ca²⁺ signals control a wide array of cellular processes. These signals require spatial and temporal regulation of the intracellular Ca²⁺ concentration, which is achieved in part by a class of ubiquitous membrane proteins known as sodium–calcium exchangers (NCXs). NCXs are secondary-active antiporters that power the translocation of Ca²⁺ across the cell membrane by coupling it to the flux of Na⁺ in the opposite direction, down an electrochemical gradient. Na⁺ and Ca²⁺ are translocated in separate steps of the antiport cycle, each of which is thought to entail a mechanism whereby ion-binding sites within the protein become alternately exposed to either side of the membrane. The prokaryotic exchanger NCX_Mj, the only member of this family with known structure, has been proposed to be a good functional and structural model of mammalian NCXs; yet our understanding of the functional properties of this protein remains incomplete. Here, we study purified NCX_Mj reconstituted into liposomes under well-controlled experimental conditions and demonstrate that this homologue indeed shares key functional features of the NCX family. Transport assays and reversal-potential measurements enable us to delineate the essential characteristics of this antiporter and establish that its ion-exchange stoichiometry is 3Na⁺:1Ca²⁺. Together with previous studies, this work confirms that NCX_Mj is a valid model system to investigate the mechanism of ion recognition and membrane transport in sodium–calcium exchangers.

Endogenous Hydrogen Sulfide Contributes to Tone Generation in Porcine Lower Esophageal Sphincter Via Na+/Ca2+ Exchanger

BACKGROUND AND AIMS

Hydrogen sulfide (H₂S) is a major physiologic gastrotransmitter. Its role in the regulation of the lower esophageal sphincter (LES) function remains unknown. The present study addresses this question.

METHODS

Isometric contraction was monitored in circular smooth muscle strips of porcine LES. Changes in cytosolic Ca²⁺ concentration ([Ca²⁺]_i) and force were simultaneously monitored in fura-2-loaded strips with front-surface fluorometry. The contribution of endogenous H₂S to LES contractility was investigated by examining the effects of inhibitors of H₂S-generating enzymes, including cystathionine-β-synthase, cystathionine-γ-lyase, and 3-mercaptopyruvate sulfurtransferase, on the LES function.

RESULTS

Porcine LES strips myogenically maintained a tetrodotoxin-resistant basal tone. Application of AOA (cystathionine-β-synthase inhibitor) or L-aspartic acid (L-Asp; 3-mercaptopyruvate sulfurtransferase inhibitor) but not DL-PAG (cystathionine-γ-lyase inhibitor), decreased this basal tone. The relaxant effects of AOA and L-Asp were additive. Maximum relaxation was obtained by combination of 1 mM AOA and 3 mM L-Asp. Immunohistochemical analyses revealed that cystathionine-β-synthase and 3-mercaptopyruvate sulfurtransferase, but not cystathionine-γ-lyase, were expressed in porcine LES. AOA+L-Asp-induced relaxation was accompanied by a decrease in [Ca²⁺]_i and inversely correlated with the extracellular Na⁺ concentration ([Na⁺]_o) (25-137.4 mM), indicating involvement of an Na⁺/Ca²⁺ exchanger. The reduction in the basal [Ca²⁺]_i level by AOA was significantly augmented in the antral smooth muscle sheets of Na⁺/Ca²⁺ exchanger transgenic mice compared with wild-type mice.

CONCLUSIONS

Endogenous H₂S regulates the LES myogenic tone by maintaining the basal [Ca²⁺]_i via Na⁺/Ca²⁺ exchanger. H₂S-generating enzymes may be a potential therapeutic target for esophageal motility disorders, such as achalasia.

Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations

The expansion of computational power, better parameterization of force fields, and the development of novel algorithms to enhance the sampling of the free energy landscapes of proteins have allowed molecular dynamics (MD) simulations to become an indispensable tool to understand the function of biomolecules. The temporal and spatial resolution of MD simulations allows for the study of a vast number of processes of interest. Here, we review the computational efforts to uncover the conformational free energy landscapes of a subset of membrane proteins: ion channels, transporters and G-protein coupled receptors. We focus on the various enhanced sampling techniques used to study these questions, how the conclusions come together to build a coherent picture, and the relationship between simulation outcomes and experimental observables.

Dynamic distinctions in the Na+/Ca2+ exchanger adopting the inward- and outward-facing conformational states

Na⁺/Ca²⁺ exchanger (NCX) proteins operate through the alternating access mechanism, where the ion-binding pocket is exposed in succession either to the extracellular or the intracellular face of the membrane. The archaeal NCX_Mj (Methanococcus jannaschii NCX) system was used to resolve the backbone dynamics in the inward-facing (IF) and outward-facing (OF) states by analyzing purified preparations of apo- and ion-bound forms of NCX_Mj-WT and its mutant, NCX_Mj-5L6-8. First, the exposure of extracellular and cytosolic vestibules to the bulk phase was evaluated as the reactivity of single cysteine mutants to a fluorescent probe, verifying that NCX_Mj-WT and NCX_Mj-5L6-8 preferentially adopt the OF and IF states, respectively. Next, hydrogen-deuterium exchange-mass spectrometry (HDX-MS) was employed to analyze the backbone dynamics profiles in proteins, preferentially adopting the OF (WT) and IF (5L6-8) states either in the presence or absence of ions. Characteristic differences in the backbone dynamics were identified between apo NCX_Mj-WT and NCX_Mj-5L6-8, thereby underscoring specific conformational patterns owned by the OF and IF states. Saturating concentrations of Na⁺ or Ca²⁺ specifically modify HDX patterns, revealing that the ion-bound/occluded states are much more stable (rigid) in the OF than in the IF state. Conformational differences observed in the ion-occluded OF and IF states can account for diversifying the ion-release dynamics and apparent affinity (K_m ) at opposite sides of the membrane, where specific structure-dynamic elements can effectively match the rates of bidirectional ion movements at physiological ion concentrations.

A selectivity filter at the intracellular end of the acid-sensing ion channel pore

Increased extracellular proton concentrations during neurotransmission are converted to excitatory sodium influx by acid-sensing ion channels (ASICs). 10-fold sodium/potassium selectivity in ASICs has long been attributed to a central constriction in the channel pore, but experimental verification is lacking due to the sensitivity of this structure to conventional manipulations. Here, we explored the basis for ion selectivity by incorporating unnatural amino acids into the channel, engineering channel stoichiometry and performing free energy simulations. We observed no preference for sodium at the “GAS belt” in the central constriction. Instead, we identified a band of glutamate and aspartate side chains at the lower end of the pore that enables preferential sodium conduction.

DOI:http://dx.doi.org/10.7554/eLife.24630.001

Structure-based dynamic arrays in regulatory domains of sodium-calcium exchanger (NCX) isoforms

Mammalian Na⁺/Ca²⁺ exchangers, NCX1 and NCX3, generate splice variants, whereas NCX2 does not. The CBD1 and CBD2 domains form a regulatory tandem (CBD12), where Ca²⁺ binding to CBD1 activates and Ca²⁺ binding to CBD2 (bearing the splicing segment) alleviates the Na⁺-induced inactivation. Here, the NCX2-CBD12, NCX3-CBD12-B, and NCX3-CBD12-AC proteins were analyzed by small-angle X-ray scattering (SAXS) and hydrogen-deuterium exchange mass-spectrometry (HDX-MS) to resolve regulatory variances in the NCX2 and NCX3 variants. SAXS revealed the unified model, according to which the Ca²⁺ binding to CBD12 shifts a dynamic equilibrium without generating new conformational states, and where more rigid conformational states become more populated without any global conformational changes. HDX-MS revealed the differential effects of the B and AC exons on the folding stability of apo CBD1 in NCX3-CBD12, where the dynamic differences become less noticeable in the Ca²⁺-bound state. Therefore, the apo forms predefine incremental changes in backbone dynamics upon Ca²⁺ binding. These observations may account for slower inactivation (caused by slower dissociation of occluded Ca²⁺ from CBD12) in the skeletal vs the brain-expressed NCX2 and NCX3 variants. This may have physiological relevance, since NCX must extrude much higher amounts of Ca²⁺ from the skeletal cell than from the neuron.

Identification of residues that control Li+ versus Na+ dependent Ca2+ exchange at the transport site of the mitochondrial NCLX

BACKGROUND

The Na⁺/Ca²⁺/Li⁺ exchanger (NCLX) is a member of the Na⁺/Ca²⁺ exchanger family. NCLX is unique in its capacity to transport both Na⁺ and Li⁺, unlike other members, which are Na⁺ selective. The major aim of this study was twofold, i.e., to identify NCLX residues that confer Li⁺ or Na⁺ selective Ca²⁺ transport and map their putative location on NCLX cation transport site.

METHOD

We combined molecular modeling to map transport site of NCLX with euryarchaeal H⁺/Ca²⁺ exchanger, CAX_Af, and fluorescence analysis to monitor Li⁺ versus Na⁺ dependent mitochondrial Ca²⁺ efflux of transport site mutants of NCLX in permeabilized cells.

RESULT

Mutation of Asn149, Pro152, Asp153, Gly176, Asn467, Ser468, Gly494 and Asn498 partially or strongly abolished mitochondrial Ca²⁺ exchange activity in intact cells. In permeabilized cells, N149A, P152A, D153A, N467Q, S468T and G494S demonstrated normal Li⁺/Ca²⁺ exchange activity but a reduced Na⁺/Ca²⁺ exchange activity. On the other hand, D471A showed dramatically reduced Li⁺/Ca²⁺ exchange, but Na⁺/Ca²⁺ exchange activity was unaffected. Finally, simultaneous mutation of four putative Ca²⁺ binding residues was required to completely abolish both Na⁺/Ca²⁺ and Li⁺/Ca²⁺ exchange activities.

CONCLUSIONS

We identified distinct Na⁺ and Li⁺ selective residues in the NCLX transport site. We propose that functional segregation in Li⁺ and Na⁺ sites reflects the functional properties of NCLX required for Ca²⁺ exchange under the unique membrane potential and ion gradient across the inner mitochondrial membrane.

GENERAL SIGNIFICANCE

The results of this study provide functional insights into the unique Li⁺ and Na⁺ selectivity of the mitochondrial exchanger. This article is part of a Special Issue entitled: ECS Meeting edited by Claus Heizmann, Joachim Krebs and Jacques Haiech.

Mitochondrial Calcium Handling in Physiology and Disease

Calcium (Ca²⁺) accumulation inside mitochondria represents a pleiotropic signal controlling a wide range of cellular functions, including key metabolic pathways and life/death decisions. This phenomenon has been first described in the 1960s, but the identity of the molecules controlling this process remained a mystery until just few years ago, when both mitochondrial Ca²⁺ uptake and release systems were genetically dissected. This finally opened the possibility to develop genetic models to directly test the contribution of mitochondrial Ca²⁺ homeostasis to cellular functions. Here we summarize our current understanding of the molecular machinery that controls mitochondrial Ca²⁺ handling and critically evaluate the physiopathological role of mitochondrial Ca²⁺ signaling, based on recent evidences obtained through in vitro and in vivo models.

Structure-Functional Basis of Ion Transport in Sodium-Calcium Exchanger (NCX) Proteins

The membrane-bound sodium-calcium exchanger (NCX) proteins shape Ca2+ homeostasis in many cell types, thus participating in a wide range of physiological and pathological processes. Determination of the crystal structure of an archaeal NCX (NCX_Mj) paved the way for a thorough and systematic investigation of ion transport mechanisms in NCX proteins. Here, we review the data gathered from the X-ray crystallography, molecular dynamics simulations, hydrogen-deuterium exchange mass-spectrometry (HDX-MS), and ion-flux analyses of mutants. Strikingly, the apo NCX_Mj protein exhibits characteristic patterns in the local backbone dynamics at particular helix segments, thereby possessing characteristic HDX profiles, suggesting structure-dynamic preorganization (geometric arrangements of catalytic residues before the transition state) of conserved α₁ and α₂ repeats at ion-coordinating residues involved in transport activities. Moreover, dynamic preorganization of local structural entities in the apo protein predefines the status of ion-occlusion and transition states, even though Na⁺ or Ca2+ binding modifies the preceding backbone dynamics nearby functionally important residues. Future challenges include resolving the structural-dynamic determinants governing the ion selectivity, functional asymmetry and ion-induced alternating access. Taking into account the structural similarities of NCX_Mj with the other proteins belonging to the Ca2+/cation exchanger superfamily, the recent findings can significantly improve our understanding of ion transport mechanisms in NCX and similar proteins.