Competition among Ca2+, Mg2+, and Na+ for model ion channel selectivity filters: determinants of ion selectivity

Comediation of voltage gating and ion charge in MXene membrane for controllable and selective monovalent cation separation

Artificial ion channels with controllable mono/monovalent cation separation fulfill important roles in biomedicine, ion separation, and energy conversion. However, it remains a daunting challenge to develop an artificial ion channel similar to biological ion channels due to ion-ion competitive transport and lack of ion-gating ability of channels. Here, we report a conductive MXene membrane with polydopamine-confined angstrom-scale channels and propose a voltage gating and ion charge comediation strategy to concurrently achieve gated and selective mono/monovalent cation separation. The membrane shows a highly switchable "on-off" ratio of ∼9.9 for K+ transport and an excellent K+/Li+ selectivity of 40.9, outperforming the ion selectivity of reported membranes with electrical gating (typically 1.5 to 6). Theoretical simulations reveal that the introduced high-charge cations such as Mg2+ enable the preferential distribution of target K+ over competing Li+ at the channel entrance, and the surface potential reduces the ionic transport energy barrier for allowing K+ to pass quickly through the channel.

Intracellular effects of lithium in aging neurons

Lithium therapy received approval during the 1970s, and it has been used for its antidepressant, antimanic, and anti-suicidal effects for acute and long-term prophylaxis and treatment of bipolar disorder (BPD). These properties have been well established; however, the molecular and cellular mechanisms remain controversial. In the past few years, many studies demonstrated that at the cellular level, lithium acts as a regulator of neurogenesis, aging, and Ca2+ homeostasis. At the molecular level, lithium modulates aging by inhibiting glycogen synthase kinase-3β (GSK-3β), and the phosphatidylinositol (PI) cycle; latter, lithium specifically inhibits inositol production, acting as a non-competitive inhibitor of inositol monophosphatase (IMPase). Mitochondria and peroxisome proliferator-activated receptor γ coactivator-1α (PGC-1α) have been related to lithium activity, and its regulation is mediated by GSK-3β degradation and inhibition. Lithium also impacts Ca2+ homeostasis in the mitochondria modulating the function of the lithium-permeable mitochondrial Na+-Ca2+exchanger (NCLX), affecting Ca2+ efflux from the mitochondrial matrix to the endoplasmic reticulum (ER). A close relationship between the protease Omi, GSK-3β, and PGC-1α has also been established. The purpose of this review is to summarize some of the intracellular mechanisms related to lithium activity and how, through them, neuronal aging could be controlled.

Binding of carboxylate and water to monovalent cations

The interactions of carboxylate anions with water and cations are important for a wide variety of systems, both biological and synthetic. To gain insight on properties of the local complexes, we apply density functional theory, to treat the complex electrostatic interactions, and investigate mixtures with varied numbers of carboxylate anions (acetate) and waters binding to monovalent cations, Li+, Na+ and K+. The optimal structure with overall lowest free energy contains two acetates and two waters such that the cation is four-fold coordinated, similar to structures found earlier for pure water or pure carboxylate ligands. More generally, the complexes with two acetates have the lowest free energy. In transitioning from the overall optimal state, exchanging an acetate for water has a lower free energy barrier than exchanging water for an acetate. In most cases, the carboxylates are monodentate and in the first solvation shell. As water is added to the system, hydrogen bonding between waters and carboxylate O atoms further stabilizes monodentate structures. These structures, which have strong electrostatic interactions that involve hydrogen bonds of varying strength, are significantly polarized, with ChelpG partial charges that vary substantially as the bonding geometry varies. Overall, these results emphasize the increasing importance of water as a component of binding sites as the number of ligands increases, thus affecting the preferential solvation of specific metal ions and clarifying Hofmeister effects. Finally, structural analysis correlated with free energy analysis supports the idea that binding to more than the preferred number of carboxylates under architectural constraints are a key to ion transport.

Protein Ca2+-Sites Prone to Sr2+ Substitution: Implications for Strontium Therapy

Strontium (Sr), an alkali metal with properties similar to calcium, in the form of soluble salts is used to treat osteoporosis. Despite the information accumulated on the role of Sr²⁺ as a Ca²⁺ mimetic in biology and medicine, there is no systematic study of how the outcome of the competition between the two dications depends on the physicochemical properties of (i) the metal ions, (ii) the first- and second-shell ligands, and (iii) the protein matrix. Specifically, the key features of a Ca²⁺-binding protein that enable Sr²⁺ to displace Ca²⁺ remain unclear. To address this, we studied the competition between Ca²⁺ and Sr²⁺ in protein Ca²⁺-binding sites using density functional theory combined with the polarizable continuum model. Our findings indicate that Ca²⁺-sites with multiple strong charge-donating protein ligands, including one or more bidentately bound Asp^-/Glu^- that are relatively buried and rigid are protected against Sr²⁺ attack. On the other hand, Ca²⁺-sites crowded with multiple protein ligands may be prone to Sr²⁺ displacement if they are solvent-exposed and flexible enough so that an extra backbone ligand from the outer shell can bind to Sr²⁺. In addition, solvent-exposed Ca²⁺ sites with only a few weak charge-donating ligands that can rearrange to fit the strontium's coordination requirements are susceptible to Sr²⁺ displacement. We provide the physical basis of these results and discuss potential novel protein targets of therapeutic Sr²⁺.

Lanthanides as Calcium Mimetic Species in Calcium-Signaling/Buffering Proteins: The Effect of Lanthanide Type on the Ca2+/Ln3+ Competition

Lanthanides, the 14 4f-block elements plus Lanthanum, have been extensively used to study the structure and biochemical properties of metalloproteins. The characteristics of lanthanides within the lanthanide series are similar, but not identical. The present research offers a systematic investigation of the ability of the entire Ln3+ series to substitute for Ca2+ in biological systems. A well-calibrated DFT/PCM protocol is employed in studying the factors that control the metal selectivity in biological systems by modeling typical calcium signaling/buffering binding sites and elucidating the thermodynamic outcome of the competition between the “alien” La3+/Ln3+ and “native” Ca2+, and La3+ − Ln3+ within the lanthanide series. The calculations performed reveal that the major determinant of the Ca2+/Ln3+ selectivity in calcium proteins is the net charge of the calcium binding pocket; the more negative the charge, the higher the competitiveness of the trivalent Ln3+ with respect to its Ca2+ contender. Solvent exposure of the binding site also influences the process; buried active centers with net charge of −4 or −3 are characterized by higher Ln3+ over Ca2+ selectivity, whereas it is the opposite for sites with overall charge of −1. Within the series, the competition between La3+ and its fellow lanthanides is determined by the balance between two competing effects: electronic (favoring heavier lanthanides) and solvation (generally favoring the lighter lanthanides).

How Theoretical Evaluations Can Generate Guidelines for Designing/Engineering Metalloproteins with Desired Metal Affinity and Selectivity

Almost half of all known proteins contain metal co-factors. Crucial for the flawless performance of a metalloprotein is the selection with high fidelity of the cognate metal cation from the surrounding biological fluids. Therefore, elucidating the factors controlling the metal binding and selectivity in metalloproteins is of particular significance. The knowledge thus acquired not only contributes to better understanding of the intimate mechanism of these events but, also, significantly enriches the researcher's toolbox that could be used in designing/engineering novel metalloprotein structures with pre-programmed properties. A powerful tool in aid of deciphering the physical principles behind the processes of metal recognition and selectivity is theoretical modeling of metal-containing biological structures. This review summarizes recent findings in the field with an emphasis on elucidating the major factors governing these processes. The results from theoretical evaluations are discussed. It is the hope that the physical principles evaluated can serve as guidelines in designing/engineering of novel metalloproteins of interest to both science and industry.

Factors allowing small monovalent Li+ to displace Ca2+ in proteins

Because Li⁺ and Ca²⁺ differ in both charge and size, the possibility that monovalent Li⁺ could dislodge the bulkier, divalent Ca²⁺ in Ca²⁺ proteins had not been considered. However, our recent density functional theory/continuum dielectric calculations predicted that Li⁺ could displace the native Ca²⁺ from the C2 domain of cytosolic PKCα/γ. This would reduce electrostatic interactions between the Li⁺-bound C2 domain and the membrane, consistent with experimental studies showing that Li⁺ can inhibit the translocation of cytoplasmic PKC to membranes. Besides the trinuclear Ca²⁺-site in the PKCα/γ C2 domain, it is not known whether other Ca²⁺-sites in human proteins may be susceptible to Li⁺ substitution. Furthermore, it is unclear what factors determine the outcome of the competition between divalent Ca²⁺ and monovalent Li⁺. Here we show that the net charge of residues in the first and second coordination shell is a key determinant of the selectivity for divalent Ca²⁺ over monovalent Li⁺ in proteins: neutral/anionic Ca²⁺-carboxylate sites are protected against Li⁺ attack. They are further protected by outer-shell Asp^-/Glu^- and the protein matrix rigidifying the Ca²⁺-site or limiting water entry. In contrast, buried, cationic Ca²⁺-sites surrounded by Arg⁺/Lys⁺, which are found in the C2 domains of PKCα/γ, as well as certain synaptotagmins, are prone to Li⁺ attack.

Theoretical Insight into the Phosphate-Targeted Silver's Antibacterial Action: Differentiation between Gram (+) and Gram (-) Bacteria

Although silver is one of the first metals finding broad applications in everyday life, specific key points of the intimate mechanism of its bacteriostatic/bactericidal activity lack explanation. It is widely accepted that the antimicrobial potential of the silver cation depends on the composition and thickness of the bacterial external envelope: the outer membrane in Gram-negative bacteria is more prone to Ag⁺ attack than the cell wall in Gram-positive bacteria. The major cellular components able to interact strongly with Ag⁺ (teichoic acids, phospholipids, and lipopolysaccharides) contain mono/diesterified phosphate moieties. By applying a reliable DFT/SMD methodology, we modeled the reactions between the aforementioned constituents in typical Gram-positive and Gram-negative bacteria and hydrated Ag⁺ species, thus disclosing the factors that govern the process of metal-model ligand complexation. The conducted research indicates thermodynamically possible reactions in all cases but still a greater preference of the Ag⁺ toward the constituents in Gram-negative bacteria in comparison with their counterparts in Gram-positive bacteria. The observed tendencies shed light on the specific interactions of the silver cation with the modeled phosphate-containing units at the atomic level.

Molecular insights into the interaction of angiotensin I-converting enzyme (ACE) inhibitors and HEXXH motif

Nutraceuticals and functional foods garner a lot of attention as potential alternative therapies for treatment of (pre)hypertension. Food-derived proteins release large variety of bioactive peptides which are similar in structure to peptide sequences acting in the organism and therefore can modulate their physiological functions. Val-Pro-Pro (VPP) is a milk-derived tripeptide with assumed mild inhibitory activity against angiotensin-converting enzyme (ACE). Computational (DFT) methods are applied on simplified models of Zn2+-HEXXH binding motif without/with bound inhibitors in order to assess the ability of two pharmaceutical drugs (Captopril and Lisinopril) and Val-Pro-Pro to coordinate with Zn2+-HEXXH binding motif of ACE. Both drugs have significant affinity towards the active site, while the Val-Pro-Pro tripeptide has weaker affinity. The obtained results shed light on the thermodynamic aspects of the inhibitors coordination to the Zn2+-HEXXH binding motif of ACE.

Metal Affinity/Selectivity of Monophosphate-Containing Signaling/Lipid Molecules

Monophosphate, an essential component of nucleic acids, as well as cell membranes and signaling molecules, is often bound to metal cations. Despite the biological importance of monophosphate-containing cell-signaling or lipid molecules, their propensity to bind the two most abundant cellular dications, Mg²⁺ and Ca²⁺, in a particular mode (inner/outer shell, mono/bidentate) is not well understood. Whether they prefer binding to Mg²⁺ than to Ca²⁺ and if they can outcompete the carboxylates of excitatory Asp/Glu and inhibitory gamma-aminobutyric acid (GABA) neurotransmitters in binding to Mg²⁺/Ca²⁺ remain unclear. To address these questions, we modeled cyclic adenosine/guanosine monophosphate (cAMP/cGMP), nucleoside 2',3'-cyclic phosphate, phosphatidylinositol (PI), phosphatidylserine (PS), and phosphatidylethanolamine (PEA) and determined their most stable metal-binding modes, including those of Asp/Glu and GABA, as well as their selectivity for Mg²⁺/Ca²⁺ using density functional theory combined with the polarizable continuum model. The results obtained, which are consistent with the available experimental findings, reveal that the structurally and functionally diverse monophosphate-containing ligands studied prefer monodentate coordination of Mg²⁺ because of the greater strain encountered upon bidentate coordination, whereas the larger Ca²⁺ imposes less strain upon bidentate binding and has reduced/no preference for monodentate coordination. We further show that in a low-dielectric environment, negatively charged monophosphate-containing ligands favor the better charge-accepting dication, that is, Mg²⁺ rather than Ca²⁺. By promoting Mg²⁺ over Ca²⁺ binding, signaling monophosphates (cAMP/cGMP) do not entrap cellular Ca²⁺ and interfere with signal transduction processes employing Ca²⁺ as a second messenger. In regions with high glutamate cytoplasmic concentration, glutamate may sequester Mg²⁺ bound to isolated five-/six-membered ring phosphates, PI, or neutral PEA, but not anionic phospholipids constituting the inner leaflet of the cell membrane.

Why Cellular Di/Triphosphates Preferably Bind Mg2+ and Not Ca2

Di/triphosphates perform a multitude of essential tasks, being important components of many vital organic cofactors such as adenosine/guanosine di/triphosphate (ADP/GDP, ATP/GTP), flavin adenine dinucleotide, and nicotinamide adenine dinucleotide and its phosphate derivative. They are generally bound to cations inside cells, in particular Mg²⁺ in the case of ATP/GTP. Yet how their metal-binding modes depend on the number, charge, and solvent exposure of the polyphosphate group and how Mg²⁺and Ca²⁺ dications that coexist in cellular fluids compete for di/triphosphates in biological systems remain elusive. Using density functional theory calculations combined with a polarizable continuum model, we have determined the relative free energies and stabilities of the different binding modes of di- and triphosphate groups to Mg²⁺ and Ca²⁺. We show that the thermodynamic outcome of the competition between Mg²⁺ and Ca²⁺ for cellular di/triphosphates depends mainly on the oligomericity/charge and metal-binding mode of the phosphate ligand as well as the solvent exposure of the binding site. Increasing the charge and thus denticity of the phosphate ligand from bi- to tridentate in a buried binding pocket enhances the affinity of the host system for the stronger charge acceptor, Mg²⁺. The cellular di/triphosphates's intrinsic properties and the protein matrix allowing them to bind a dication bi/tridentately, along with the higher cytosolic concentration of Mg²⁺ compared to Ca²⁺, enables Mg²⁺ to outcompete Ca²⁺ in binding to these highly charged anions. This suggests an explanation for why nature has chosen Mg²⁺ but not Ca²⁺ to perform most of the essential tasks associated with biological triphosphates.

Free and Bound Therapeutic Lithium in Brain Signaling

Lithium, a first-line therapy for bipolar disorder, is effective in preventing suicide and new depressive/manic episodes. Yet, how this beguilingly simple monocation with only two electrons could yield such profound therapeutic effects remains unclear. An in-depth understanding of lithium's mechanisms of actions would help one to develop better treatments limiting its adverse side effects and repurpose lithium for treating traumatic brain injury and chronic neurodegenerative diseases. In this Account, we begin with a comparison of the physicochemical properties of Li⁺ and its key native rivals, Na⁺ and Mg²⁺, to provide physical grounds for their competition in protein binding sites. Next, we review the abnormal signaling pathways and proteins found in bipolar patients, who generally have abnormally high intracellular Na⁺ and Ca²⁺ concentrations, high G-protein levels, and hyperactive phosphatidylinositol signaling and glycogen synthase kinase-3β (GSK3β) activity. We briefly summarize experimental findings on how lithium, at therapeutic doses, modulates these abnormal signaling pathways and proteins. Following this survey, we address the following aspects of lithium's therapeutic actions: (1) Can Li⁺ displace Na⁺ from the allosteric Na⁺-binding sites in neurotransmitter transporters and G-protein coupled receptors (GPCRs); if so, how would this affect the host protein's function? (2) Why are certain Mg²⁺-dependent enzymes targeted by Li⁺? (3) How does Li⁺ binding to Mg²⁺-bound ATP/GTP (denoted as NTP) in solution affect the cofactor's conformation and subsequent recognition by the host protein? (4) How do NTP-Mg-Li complexes modulate the properties of the respective cellular receptors and signal-transducing proteins? We show that Li⁺ may displace Na⁺ from allosteric Na⁺-binding sites in certain GPCRs and stabilize inactive conformations, preventing these receptors from relaying signal to the respective G-proteins. It may also displace Mg²⁺ in enzymes containing highly cationic Mg²⁺-binding sites such as GSK3β, but not in enzymes containing Mg²⁺-binding sites with low or zero charge. We further show that Li⁺ binding to Mg²⁺-NTP in water does not alter the NTP conformation, which is locked by all three phosphates binding to Mg²⁺. However, bound lithium in the form of [NTP-Mg-Li]^2- dianions can activate or inhibit the host protein depending on the NTP-binding pocket's shape, which determines the metal-binding mode: The ATP-binding pocket's shape in the P2X receptor is complementary to the native ATP-Mg solution conformation and nicely fits [ATP-Mg-Li]^2-. However, since the ATP βγ phosphates bind Li⁺, bimetallic [ATP-Mg-Li]^2- may be more resistant to hydrolysis than the native cofactor, enabling ATP to reside longer in the binding site and elicit a prolonged P2X response. In contrast, the elongated GTP-binding pockets in G-proteins allow only two GTP phosphates to bind Mg²⁺, so the GTP conformation is no longer "triply-locked". Consequently, Li⁺ binding to GTP-Mg can significantly alter the native cofactor's structure, lowering the activated G-protein level, thus attenuating hyperactive G-protein-mediated signaling in bipolar patients. In summary, we have presented a larger "connected" picture of lithium's diverse effects based on its competition as a free monocation with native cations or as a phosphate-bound polyanionic complex modulating the host protein function.

Novel Insights into Gallium's Mechanism of Therapeutic Action: A DFT/PCM Study of the Interaction between Ga3+ and Ribonucleotide Reductase Substrates

The broadly accepted mechanism of gallium's therapeutic action postulates the inactivation of the upregulated/hyperactive enzyme ribonucleotide reductase (RNR) in cancer cells by substituting the redox-active iron by redox-silent gallium in the enzyme active site. Recently, another hypothesis for the Ga³⁺ curative effect has been put forward: the metal cation can deactivate the enzyme by entrapping its substrates (nucleotide diphosphates; NDPs) into Ga³⁺-NDP complexes, lowering the free substrate levels in the cell. Several questions arise: Does gallium readily form complexes with NDPs? What are the preferable modes of metal binding to NDPs? Does, and if so, to what extent, the metal binding alter the native conformation of the substrate, thus influencing the process of substrate-enzyme recognition? Here, by employing density functional theory (DFT)/polarizable continuum model (PCM) calculations, we attempt to answer these questions. The results, which are in line with the available experimental data, lay support to the recent hypothesis about the curative effect of gallium, revealing that, by engaging the free NDPs in forming metal complexes, on the one side, and producing metal constructs that are not/poorly recognizable by the host enzyme, on the other side, gallium deprives RNR from its substrates, thus reducing the enzyme activity in malignant cells.

Electric field influence on the helical structure of peptides: insights from DFT/PCM computations

The switching of the electric field with a particular directionality could be used for the healing of misfolded proteins.

Extracellular and Intracellular Biomineralization Induced by Bacillus licheniformis DB1-9 at Different Mg/Ca Molar Ratios

Biomineralization has become a research hotspot and attracted widespread attention in the field of carbonate sedimentology. In this study, precipitation of carbonate minerals was induced by Bacillus licheniformis DB1-9 bacteria, (identity confirmed with its phylogenetic tree), to further explore the biomineralization mechanisms. During experiments, lasting up to 24 days with varying Mg/Ca molar ratios and regular monitoring of conditions, ammonia and carbonic anhydrase are released by the bacteria, resulting in a pH increase. Carbonic anhydrase could have promoted carbon dioxide hydration to produce bicarbonate and carbonate ions, and so promoted supersaturation to facilitate the precipitation of carbonate minerals. These include rhombohedral, dumbbell-shaped, and elongated calcite crystals; aragonite appears in the form of mineral aggregates. In addition, spheroidal and fusiform minerals are precipitated. FTIR results show there are organic functional groups, such as C–O–C and C=O, as well as the characteristic peaks of calcite and aragonite; these indicate that there is a close relationship between the bacteria and the minerals. Ultrathin slices of the bacteria analyzed by HRTEM, SAED, EDS, and STEM show that precipitate within the extracellular polymeric substances (EPS) has a poor crystal structure, and intracellular granular areas have no crystal structure. Fluorescence intensity and STEM results show that calcium ions can be transported from the outside to the inside of the cells. This study provides further insights to our understanding of biomineralization mechanisms induced by microorganisms.

Ion Permeation Mechanism in Epithelial Calcium Channel TRVP6

Calcium is the most abundant metal in the human body that plays vital roles as a cellular electrolyte as well as the smallest and most frequently used signaling molecule. Calcium uptake in epithelial tissues is mediated by tetrameric calcium-selective transient receptor potential (TRP) channels TRPV6 that are implicated in a variety of human diseases, including numerous forms of cancer. We used TRPV6 crystal structures as templates for molecular dynamics simulations to identify ion binding sites and to study the permeation mechanism of calcium and other ions through TRPV6 channels. We found that at low Ca²⁺ concentrations, a single calcium ion binds at the selectivity filter narrow constriction formed by aspartates D541 and allows Na⁺ permeation. In the presence of ions, no water binds to or crosses the pore constriction. At high Ca²⁺ concentrations, calcium permeates the pore according to the knock-off mechanism that includes formation of a short-lived transition state with three calcium ions bound near D541. For Ba²⁺, the transition state lives longer and the knock-off permeation occurs slower. Gd³⁺ binds at D541 tightly, blocks the channel and prevents Na⁺ from permeating the pore. Our results provide structural foundations for understanding permeation and block in tetrameric calcium-selective ion channels.

Competition between Li+ and Na+ in sodium transporters and receptors: Which Na+-Binding sites are "therapeutic" Li+ targets?

Li⁺ (turquoise), the better charge acceptor, can displace Na⁺ (purple) bound by only one or two aa residues in buried sites. Thus, Li⁺ can displace Na⁺ bound by Asp^– and Ser in the A_2AAR/β₁AR receptor and enhance the metal site's stability, thus prohibiting structural distortions induced by agonist binding, leading to lower cytosolic levels of activated G-proteins, which are hyperactive in bipolar disorder patients.

Sodium (Na⁺) acts as an indispensable allosteric regulator of the activities of biologically important neurotransmitter transporters and G-protein coupled receptors (GPCRs), which comprise well-known drug targets for psychiatric disorders and addictive behavior. How selective these allosteric Na⁺-binding sites are for the cognate cation over abiogenic Li⁺, a first-line drug to treat bipolar disorder, is unclear. Here, we reveal how properties of the host protein and its binding cavity affect the outcome of the competition between Li⁺ and Na⁺ for allosteric binding sites in sodium transporters and receptors. We show that rigid Na⁺-sites that are crowded with multiple protein ligands are well-protected against Li⁺ attack, but their flexible counterparts or buried Na⁺-sites containing only one or two protein ligands are vulnerable to Li⁺ substitution. These findings suggest a novel possible mode of Li⁺ therapeutic action: By displacing Na⁺ bound by ≤2 protein ligands in buried GPCR sites and stabilizing the receptor's inactive state, Li⁺ could prohibit conformational changes to an active state, leading to lower cytosolic levels of activated guanine nucleotide-binding proteins, which are hyperactive/overexpressed in bipolar disorder patients.

How an electric field can modulate the metal ion selectivity of protein binding sites: insights from DFT/PCM calculations

An electric field (internal or external) is a potent force that can modulate the metal selectivity of a protein binding site.

Lysine and the Na+/K+ Selectivity in Mammalian Voltage-Gated Sodium Channels

Voltage-gated sodium (Na_v) channels are critical in the generation and transmission of neuronal signals in mammals. The crystal structures of several prokaryotic Na_v channels determined in recent years inspire the mechanistic studies on their selection upon the permeable cations (especially between Na⁺ and K⁺ ions), a property that is proposed to be mainly determined by residues in the selectivity filter. However, the mechanism of cation selection in mammalian Na_v channels lacks direct explanation at atomic level due to the difference in amino acid sequences between mammalian and prokaryotic Na_v homologues, especially at the constriction site where the DEKA motif has been identified to determine the Na⁺/K⁺ selectivity in mammalian Na_v channels but is completely absent in the prokaryotic counterparts. Among the DEKA residues, Lys is of the most importance since its mutation to Arg abolishes the Na⁺/K⁺ selectivity. In this work, we modeled the pore domain of mammalian Na_v channels by mutating the four residues at the constriction site of a prokaryotic Na_v channel (Na_vRh) to DEKA, and then mechanistically investigated the contribution of Lys in cation selection using molecular dynamics simulations. The DERA mutant was generated as a comparison to understand the loss of ion selectivity caused by the K-to-R mutation. Simulations and free energy calculations on the mutants indicate that Lys facilitates Na⁺/K⁺ selection by electrostatically repelling the cation to a highly Na⁺-selective location sandwiched by the carboxylate groups of Asp and Glu at the constriction site. In contrast, the electrostatic repulsion is substantially weakened when Lys is mutated to Arg, because of two intrinsic properties of the Arg side chain: the planar geometric design and the sparse charge distribution of the guanidine group.

Gallium as a Therapeutic Agent: A Thermodynamic Evaluation of the Competition between Ga(3+) and Fe(3+) Ions in Metalloproteins

Gallium has been employed (in the form of soluble salts) to fight various forms of cancer, infectious, and inflammatory diseases. The rationale behind this lies in the ability of Ga(3+) cation to mimic closely in appearance the native ferric ion, Fe(3+), thus interfering with the biological processes requiring ferric cofactors. However, Ga(3+) ion cannot participate in redox reactions and, when substituting for the "native" Fe(3+) ion in the enzyme active site, renders it inactive. Although a significant body of information on the Ga(3+)-Fe(3+) competition in biological systems has been accumulated, the intimate mechanism of the process is still not well understood and several questions remain: What are the basic physical principles governing the competition between the two trivalent cations in proteins? What type of metal centers are the most likely targets for gallium therapy? To what extent are the Fe(3+)-binding sites in the key enzyme ribonucleotide reductase vulnerable to Ga(3+) substitution? Here, we address these questions by studying the competition between Ga(3+) and Fe(3+) ions in model metal binding sites of various compositions and charge states. The results obtained are in line with available experimental data and shed light on the intimate mechanism of the Ga(3+)/Fe(3+) selectivity in various model metal binding sites and biological systems such as serum transferrin and ribonucleotide reductase.

Factors controlling the selectivity for Na+over Mg2+in sodium transporters and enzymes

The paper discloses the key factors and physical bases that render a given binding site either Mg²⁺or Na⁺-selective.

Factors Governing the Bridging Water Protonation State in Polynuclear Mg(2+) Proteins

An aqua ligand bridges metal cations in a wide variety of enzymes, many of which are drug targets for various diseases. However, the factors affecting its protonation state and thus biological roles remain elusive. By computing the free energy for replacing the bridging H2O by OH(-) in various model Mg(2+) sites, we have evaluated how the nature of an aqua bridge depends on the site's net charge (i.e., the number of charged ligands in the first and second shell and the number of metal cations), the site's solvent exposure, the ligand's charge-donating ability, the bridging oxygen's hydrogen-bonding interactions, intramolecular proton transfer from the bridging H2O to a nearby carboxylate, and the metal coordination number. The results reveal the key factors dictating the protonation state of bridging H2O and provide guidelines in predicting whether H2O or OH(-) bridges two Mg(2+) in polynuclear sites. This helps to elucidate the nucleophile in the enzyme-catalyzed reaction and the net charge of the metal complex (metal cation and first-shell ligands), which plays a critical role in binding.

Selectivity Mechanism of the Voltage-gated Proton Channel, HV1

Voltage-gated proton channels, HV1, trigger bioluminescence in dinoflagellates, enable calcification in coccolithophores, and play multifarious roles in human health. Because the proton concentration is minuscule, exquisite selectivity for protons over other ions is critical to HV1 function. The selectivity of the open HV1 channel requires an aspartate near an arginine in the selectivity filter (SF), a narrow region that dictates proton selectivity, but the mechanism of proton selectivity is unknown. Here we use a reduced quantum model to elucidate how the Asp-Arg SF selects protons but excludes other ions. Attached to a ring scaffold, the Asp and Arg side chains formed bidentate hydrogen bonds that occlude the pore. Introducing H3O(+) protonated the SF, breaking the Asp-Arg linkage and opening the conduction pathway, whereas Na(+) or Cl(-) was trapped by the SF residue of opposite charge, leaving the linkage intact, thus preventing permeation. An Asp-Lys SF behaved like the Asp-Arg one and was experimentally verified to be proton-selective, as predicted. Hence, interacting acidic and basic residues form favorable AspH(0)-H2O(0)-Arg(+) interactions with hydronium but unfavorable Asp(-)-X(-)/X(+)-Arg(+) interactions with anions/cations. This proposed mechanism may apply to other proton-selective molecules engaged in bioenergetics, homeostasis, and signaling.

Ion selectivity in the selectivity filters of acid-sensing ion channels

Sodium-selective acid sensing ion channels (ASICs), which belong to the epithelial sodium channel (ENaC) superfamily, are key players in many physiological processes (e.g. nociception, mechanosensation, cognition, and memory) and are potential therapeutic targets. Central to the ASIC's function is its ability to discriminate Na⁺ among cations, which is largely determined by its selectivity filter, the narrowest part of an open pore. However, it is unclear how the ASIC discriminates Na⁺ from rival cations such as K⁺ and Ca²⁺ and why its Na⁺/K⁺ selectivity is an order of magnitude lower than that of the ENaC. Here, we show that a well-tuned balance between electrostatic and solvation effects controls ion selectivity in the ASIC1a SF. The large, water-filled ASIC1a pore is selective for Na⁺ over K⁺ because its backbone ligands form more hydrogen-bond contacts and stronger electrostatic interactions with hydrated Na⁺ compared to hydrated K⁺. It is selective for Na⁺ over divalent Ca²⁺ due to its relatively high-dielectric environment, which favors solvated rather than filter-bound Ca²⁺. However, higher Na⁺-selectivity could be achieved in a narrow, rigid pore lined by three weak metal-ligating groups, as in the case of ENaC, which provides optimal fit and interactions for Na⁺ but not for non-native ions.

Ion selectivity strategies of sodium channel selectivity filters

CONSPECTUS: Sodium ion channels selectively transport Na(+) cations across the cell membrane. These integral parts of the cell machinery are implicated in regulating the cardiac, skeletal and smooth muscle contraction, nerve impulses, salt and water homeostasis, as well as pain and taste perception. Their malfunction often results in various channelopathies of the heart, brain, skeletal muscles, and lung; thus, sodium channels are key drug targets for various disorders including cardiac arrhythmias, heart attack, stroke, migraine, epilepsy, pain, cancer, and autoimmune disorders. The ability of sodium channels to discriminate the native Na(+) among other competing ions in the surrounding fluids is crucial for proper cellular functions. The selectivity filter (SF), the narrowest part of the channel's open pore, lined with amino acid residues that specifically interact with the permeating ion, plays a major role in determining Na(+) selectivity. Different sodium channels have different SFs, which vary in the symmetry, number, charge, arrangement, and chemical type of the metal-ligating groups and pore size: epithelial/degenerin/acid-sensing ion channels have generally trimeric SFs lined with three conserved neutral serines and/or backbone carbonyls; eukaryotic sodium channels have EKEE, EEKE, DKEA, and DEKA SFs with an invariant positively charged lysine from the second or third domain; and bacterial voltage-gated sodium (Nav) channels exhibit symmetrical EEEE SFs, reminiscent of eukaryotic voltage-gated calcium channels. How do these different sodium channel SFs achieve high selectivity for Na(+) over its key rivals, K(+) and Ca(2+)? What factors govern the metal competition in these SFs and which of these factors are exploited to achieve Na(+) selectivity in the different sodium channel SFs? The free energies for replacing K(+) or Ca(2+) bound inside different model SFs with Na(+), evaluated by a combination of density functional theory and continuum dielectric calculations, have shed light on these questions. The SFs of epithelial and eukaryotic Nav channels select Na(+) by providing an optimal number and ligating strength of metal ligands as well as a rigid pore whose size fits the cognate Na(+) ideally. On the other hand, the SFs of bacterial Nav channels select Na(+), as the protein matrix attenuates ion-protein interactions relative to ion-solvent interactions by enlarging the pore and allowing water to enter, so the ion interacts indirectly with the conserved glutamates via bridging water molecules. This shows how these various SFs have adapted to the specific physicochemical properties of the native ion, using different strategies to select Na(+) among its contenders.

Molecular strategies to achieve selective conductance in NaK channel variants

A recent crystallization of several ion channels has provided strong impetus for efforts aimed at understanding the different strategies employed by nature for selective ion transport. In this work, we used two variants of the selectivity filter of NaK channel to explore molecular mechanisms that give rise to K(+)-selectivity. We computed one-dimensional (1D) and two-dimensional (2D) potentials of mean force (PMFs) for ion permeation across the channel. The results indicate that the energies for Na(+) and K(+) permeation across the selectivity filter display significant differences in positions of the binding sites and barriers. One characteristic signature of a K(+)-selective channel is the apparent preservation of the site analogous to that of S2 in KcsA. The S2-bound ion can be almost ideally dehydrated and coordinated by 6 to 8 carbonyls. In a striking contrast, the PMFs controlling transport of ions in a nonselective variant show almost identical profiles for either K(+) or Na(+) and significant involvement of water molecules in ion coordination across the entire selectivity filter. An analysis of differences in 1D PMFs for Na(+) and K(+) suggests that coordination number alone is an insufficient predictor of site selectivity, while chemical composition (ratio of carbonyls and water molecules) correlates well with preference for K(+). Multi-ion effects such as dependence of the barriers and wells for permeant ion on the type of copermeant ion were found to play a significant role in the selectivity signature of the channel as well.

The open gate of the K(V)1.2 channel: quantum calculations show the key role of hydration

The open gate of the Kv1.2 voltage-gated potassium channel can just hold a hydrated K(+) ion. Quantum calculations starting from the x-ray coordinates of the channel confirm this, showing little change from the x-ray coordinates for the protein. Water molecules not in the x-ray coordinates, and the ion itself, are placed by the calculation. The water molecules, including their orientation and hydrogen bonding, with and without an ion, are critical for the path of the ion, from the solution to the gate. A sequence of steps is postulated in which the potential experienced by the ion in the pore is influenced by the position of the ion. The gate structure, with and without the ion, has been optimized. The charges on the atoms and bond lengths have been calculated using natural bond orbital calculations, giving K(+) ~0.77 charges, rather than 1.0. The PVPV hinge sequence has been mutated in silico to PVVV (P407V in the 2A79 numbering). The water structure around the ion becomes discontinuous, separated into two sections, above and below the ion. PVPV conservation closely relates to maintaining the water structure. Finally, these results have implications concerning gating.

Molecular mechanisms for generating transmembrane proton gradients

Membrane proteins use the energy of light or high energy substrates to build a transmembrane proton gradient through a series of reactions leading to proton release into the lower pH compartment (P-side) and proton uptake from the higher pH compartment (N-side). This review considers how the proton affinity of the substrates, cofactors and amino acids are modified in four proteins to drive proton transfers. Bacterial reaction centers (RCs) and photosystem II (PSII) carry out redox chemistry with the species to be oxidized on the P-side while reduction occurs on the N-side of the membrane. Terminal redox cofactors are used which have pKas that are strongly dependent on their redox state, so that protons are lost on oxidation and gained on reduction. Bacteriorhodopsin is a true proton pump. Light activation triggers trans to cis isomerization of a bound retinal. Strong electrostatic interactions within clusters of amino acids are modified by the conformational changes initiated by retinal motion leading to changes in proton affinity, driving transmembrane proton transfer. Cytochrome c oxidase (CcO) catalyzes the reduction of O2 to water. The protons needed for chemistry are bound from the N-side. The reduction chemistry also drives proton pumping from N- to P-side. Overall, in CcO the uptake of 4 electrons to reduce O2 transports 8 charges across the membrane, with each reduction fully coupled to removal of two protons from the N-side, the delivery of one for chemistry and transport of the other to the P-side.

The mechanism of Na⁺/K⁺ selectivity in mammalian voltage-gated sodium channels based on molecular dynamics simulation

Voltage-gated sodium (Nav) channels and their Na⁺/K⁺ selectivity are of great importance in the mammalian neuronal signaling. According to mutational analysis, the Na⁺/K⁺ selectivity in mammalian Nav channels is mainly determined by the Lys and Asp/Glu residues located at the constriction site within the selectivity filter. Despite successful molecular dynamics simulations conducted on the prokaryotic Nav channels, the lack of Lys at the constriction site of prokaryotic Nav channels limits how much can be learned about the Na⁺/K⁺ selectivity in mammalian Nav channels. In this work, we modeled the mammalian Nav channel by mutating the key residues at the constriction site in a prokaryotic Nav channel (NavRh) to its mammalian counterpart. By simulating the mutant structure, we found that the Na⁺ preference in mammalian Nav channels is collaboratively achieved by the deselection from Lys and the selection from Asp/Glu within the constriction site.

Distinct interactions of Na+ and Ca2+ ions with the selectivity filter of the bacterial sodium channel Na(V)Ab

Rapid and selective ion transport is essential for the generation and regulation of electrical signaling pathways in living organisms. In this study, we use molecular dynamics simulations and free energy calculations to investigate how the bacterial sodium channel Na(V)Ab (Arcobacter butzleri) differentiates between Na(+) and Ca(2+) ions. Multiple nanosecond molecular dynamics simulations revealed distinct binding patterns for these two cations in the selectivity filter and suggested a high affinity calcium binding site formed by backbone atoms of residues Leu-176 and Thr-175 (S(CEN)) in the sodium channel selectivity filter.