Protonation of key acidic residues is critical for the K⁺-selectivity of the Na/K pump

Specific protonation of acidic residues confers K+ selectivity to the gastric proton pump

The gastric proton pump (H+,K+-ATPase) transports a proton into the stomach lumen for every K+ ion exchanged in the opposite direction. In the lumen-facing state of the pump (E2), the pump selectively binds K+ despite the presence of a 10-fold higher concentration of Na+. The molecular basis for the ion selectivity of the pump is unknown. Using molecular dynamics simulations, free energy calculations, and Na+ and K+-dependent ATPase activity assays, we demonstrate that the K+ selectivity of the pump depends upon the simultaneous protonation of the acidic residues E343 and E795 in the ion-binding site. We also show that when E936 is protonated, the pump becomes Na+ sensitive. The protonation-mimetic mutant E936Q exhibits weak Na+-activated ATPase activity. A 2.5-Å resolution cryo-EM structure of the E936Q mutant in the K+-occluded E2-Pi form shows, however, no significant structural difference compared with wildtype except less-than-ideal coordination of K+ in the mutant. The selectivity toward a specific ion correlates with a more rigid and less fluctuating ion-binding site. Despite being exposed to a pH of 1, the fundamental principle driving the K+ ion selectivity of H+,K+-ATPase is similar to that of Na+,K+-ATPase: the ionization states of the acidic residues in the ion-binding sites determine ion selectivity. Unlike the Na+,K+-ATPase, however, protonation of an ion-binding glutamate residue (E936) confers Na+ sensitivity.

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.

Protein pKa Prediction by Tree-Based Machine Learning

Protonation states of ionizable protein residues modulate many essential biological processes. For correct modeling and understanding of these processes, it is crucial to accurately determine their pKa values. Here, we present four tree-based machine learning models for protein pKa prediction. The four models, Random Forest, Extra Trees, eXtreme Gradient Boosting (XGBoost), and Light Gradient Boosting Machine (LightGBM), were trained on three experimental PDB and pKa datasets, two of which included a notable portion of internal residues. We observed similar performance among the four machine learning algorithms. The best model trained on the largest dataset performs 37% better than the widely used empirical pKa prediction tool PROPKA and 15% better than the published result from the pKa prediction method DelPhiPKa. The overall root-mean-square error (RMSE) for this model is 0.69, with surface and buried RMSE values being 0.56 and 0.78, respectively, considering six residue types (Asp, Glu, His, Lys, Cys, and Tyr), and 0.63 when considering Asp, Glu, His, and Lys only. We provide pKa predictions for proteins in human proteome from the AlphaFold Protein Structure Database and observed that 1% of Asp/Glu/Lys residues have highly shifted pKa values close to the physiological pH.

Prediction of protein pK a with representation learning

The behavior of proteins is closely related to the protonation states of the residues. Therefore, prediction and measurement of pK a are essential to understand the basic functions of proteins. In this work, we develop a new empirical scheme for protein pK a prediction that is based on deep representation learning. It combines machine learning with atomic environment vector (AEV) and learned quantum mechanical representation from ANI-2x neural network potential (J. Chem. Theory Comput. 2020, 16, 4192). The scheme requires only the coordinate information of a protein as the input and separately estimates the pK a for all five titratable amino acid types. The accuracy of the approach was analyzed with both cross-validation and an external test set of proteins. Obtained results were compared with the widely used empirical approach PROPKA. The new empirical model provides accuracy with MAEs below 0.5 for all amino acid types. It surpasses the accuracy of PROPKA and performs significantly better than the null model. Our model is also sensitive to the local conformational changes and molecular interactions.

Enantioselective Antiport in Asymmetric Nanochannels

Enantioselective sensing and separation are major challenges. Nanochannel technologies are energy-saving and efficient for membrane separation. Herein, inspired by biological antiporter proteins, artificial nanochannels with antiporter behavior were fabricated for chiral sensing and separation. Tyrosine enantiomers were incorporated into hourglass-shaped nanochannels via stepwise modifications to fabricating multiligand-modified asymmetric channels. Chiral distinction of naproxen enantiomers was amplified in the l-Tyr/d-Tyr channels, with an enantioselectivity coefficient of 524, which was over 100-fold that of one-ligand-modified nanochannels. Furthermore, transport experiments evidenced the spontaneous antiport of naproxen enantiomers in the l-Tyr/d-Tyr channels. The racemic naproxen sample was separated via the chiral antiport process, with an enantiomeric excess of 71.2%. Further analysis using electro-osmotic flow experiments and finite-element simulations confirmed that the asymmetric modified multiligand was key to achieving separation of the naproxen enantiomers. We expect these multiligand-modified asymmetric nanochannels to provide insight into mimicking biological antiporter systems and offer an approach to energy-efficient and robust enantiomer separation.

Beneficent and Maleficent Effects of Cations on Bufadienolide Binding to Na+,K+-ATPase

Kinetic properties and crystal structures of the Na⁺,K⁺-ATPase in complex with cardiotonic steroids (CTS) revealed significant differences between CTS subfamilies (Laursen et al.). Thus, we found beneficial effects of K⁺ on bufadienolide binding, which strongly contrasted with the well-known antagonism between K⁺ and cardenolides. In order to understand this peculiarity of bufalin interactions, we used docking and molecular dynamics simulations of the complexes involving Na⁺,K⁺-ATPase, bufadienolides (bufalin, cinobufagin), and ions (K⁺, Na⁺, Mg²⁺). The results revealed that bufadienolide binding is affected by (i) electrostatic attraction of the lactone ring by a cation and (ii) the ability of a cation to stabilize and "shape" the site constituted by transmembrane helices of the α-subunit (αM1-6). The latter effect was due to varying coordination patterns involving amino acid residues from helix bundles αM1-4 and αM5-10. Substituents on the steroid core of a bufadienolide add to and modify the cation effects. The above rationale is fully consistent with the ion effects on the kinetics of Na⁺,K⁺-ATPase/bufadienolide interactions.

Corneal Abnormalities Are Novel Clinical Feature in Wolfram Syndrome

PURPOSE

To evaluate corneal morphology among patients with Wolfram syndrome (WFS).

DESIGN

Comparative observational longitudinal case series of WFS patients with a laboratory approach in the WFS1 gene knockout (Wfs1KO) mouse model.

METHODS

A group of 12 patients with biallelic mutations in the WFS1 gene recruited from the whole country and a control group composed of 30 individuals with type 1 diabetes (T1D) were evaluated in a national reference center for monogenic diabetes. All subjects (n = 42) underwent a complete ophthalmic examination, computer videokeratography, and corneal thickness and endothelial measurements. Additionally, WFS patients (n = 9) underwent longitudinal videokeratography and Pentacam evaluation. Corneal characteristics were assessed and compared between both groups. Human and mouse corneas were subjected to immunohistochemistry to detect wolframin expression and microscopic evaluation to study corneal morphology ex vivo.

RESULTS

Clinical and topographic abnormalities similar to keratoconus were observed in 14 eyes (58.3%) of 8 WFS patients (66.7%). Flat keratometry, inferior-superior dioptric asymmetry, skewed radial axis, logarithm of keratoconus percentage index, index of surface variance, index of vertical asymmetry, keratoconus index, central keratoconus index, index of height asymmetry, and index of height decentration differed between WFS and T1D patients. Immunohistochemistry demonstrated wolframin expression in human and mouse corneas. Compared with Wfs1WT mice, Wfs1KO mice also presented corneal abnormalities.

CONCLUSIONS

Patients with WFS present a high prevalence of changes in corneal morphology compatible with the diagnosis of early stages of keratoconus. Observations in a mouse model suggest that a mutation in the WFS1 gene may be responsible for corneal abnormalities similar to keratoconus.

Fully Atomistic Multiscale Approach for pKa Prediction

The ionization state of titratable amino acids strongly affects proteins structure and functioning in a large number of biological processes. It is therefore essential to be able to characterize the pK_a of ionizable groups inside proteins and to understand its microscopic determinants in order to gain insights into many functional properties of proteins. A big effort has been devoted to the development of theoretical approaches for the prediction of deprotonation free energies, yet the accurate theoretical/computational calculation of pK_a values is recognized as a current challenge. A methodology based on a hybrid quantum/classical approach is here proposed for the computation of deprotonation free energies. The method is applied to calculate the pK_a of formic acid, methylammonium, and methanethiol, providing results in good agreement with the corresponding experimental estimates. The pK_a is also calculated for aspartic acid and lysine as single residues in solution and for three aspartic/glutamic acids inside a well-characterized protein: hen egg white lysozyme. While for small molecules the method is able to deal with multiple protonation states of all titratable groups, this becomes computationally very expensive for proteins. The calculated pK_a values for the single amino acids (except for the zwitterionic aspartic acid) and inside the protein display a systematic shift with respect to the experimental values that suggests that the fine balance between hydrophobic and polar interactions might be not accurately reproduced by the usual classical force-fields, thus affecting the computation of deprotonation free energies. The calculated pK_a shifts inside the protein are in good agreement with the corresponding experimental ones (within 1 pK_a unit), well reproducing the pK_a changes due to the protein environment even in the case of large pK_a shifts.

Allosteric gate modulation confers K+ coupling in glutamate transporters

Excitatory amino acid transporters (EAATs) mediate glial and neuronal glutamate uptake to terminate synaptic transmission and to ensure low resting glutamate concentrations. Effective glutamate uptake is achieved by cotransport with 3 Na⁺ and 1 H⁺, in exchange with 1 K⁺. The underlying principles of this complex transport stoichiometry remain poorly understood. We use molecular dynamics simulations and electrophysiological experiments to elucidate how mammalian EAATs harness K⁺ gradients, unlike their K⁺‐independent prokaryotic homologues. Glutamate transport is achieved via elevator‐like translocation of the transport domain. In EAATs, glutamate‐free re‐translocation is prevented by an external gate remaining open until K⁺ binding closes and locks the gate. Prokaryotic Glt_Ph contains the same K⁺‐binding site, but the gate can close without K⁺. Our study provides a comprehensive description of K⁺‐dependent glutamate transport and reveals a hitherto unknown allosteric coupling mechanism that permits adaptions of the transport stoichiometry without affecting ion or substrate binding.

A single K+-binding site in the crystal structure of the gastric proton pump

The gastric proton pump (H⁺,K⁺-ATPase), a P-type ATPase responsible for gastric acidification, mediates electro-neutral exchange of H⁺ and K⁺ coupled with ATP hydrolysis, but with an as yet undetermined transport stoichiometry. Here we show crystal structures at a resolution of 2.5 Å of the pump in the E2-P transition state, in which the counter-transporting cation is occluded. We found a single K⁺ bound to the cation-binding site of the H⁺,K⁺-ATPase, indicating an exchange of 1H⁺/1K⁺ per hydrolysis of one ATP molecule. This fulfills the energy requirement for the generation of a six pH unit gradient across the membrane. The structural basis of K⁺ recognition is resolved and supported by molecular dynamics simulations, establishing how the H⁺,K⁺-ATPase overcomes the energetic challenge to generate an H⁺ gradient of more than a million-fold-one of the highest cation gradients known in mammalian tissue-across the membrane.

Cholesterol depletion inhibits Na+,K+-ATPase activity in a near-native membrane environment

Cholesterol's effects on Na+,K+-ATPase reconstituted in phospholipid vesicles have been extensively studied. However, previous studies have reported both cholesterol-mediated stimulation and inhibition of Na+,K+-ATPase activity. Here, using partial reaction kinetics determined via stopped-flow experiments, we studied cholesterol's effect on Na+,K+-ATPase in a near-native environment in which purified membrane fragments were depleted of cholesterol with methyl-β-cyclodextrin (mβCD). The mβCD-treated Na+,K+-ATPase had significantly reduced overall activity and exhibited decreased observed rate constants for ATP phosphorylation (ENa3+ → E2P, i.e. phosphorylation by ATP and Na+ occlusion from the cytoplasm) and K+ deocclusion with subsequent intracellular Na+ binding (E2K2+ → E1Na3+). However, cholesterol depletion did not affect the observed rate constant for K+ occlusion by phosphorylated Na+,K+-ATPase on the extracellular face and subsequent dephosphorylation (E2P → E2K2+). Thus, partial reactions involving cation binding and release at the protein's intracellular side were most dependent on cholesterol. Fluorescence measurements with the probe eosin indicated that cholesterol depletion stabilizes the unphosphorylated E2 state relative to E1, and the cholesterol depletion-induced slowing of ATP phosphorylation kinetics was consistent with partial conversion of Na+,K+-ATPase into the E2 state, requiring a slow E2 → E1 transition before the phosphorylation. Molecular dynamics simulations of Na+,K+-ATPase in membranes with 40 mol % cholesterol revealed cholesterol interaction sites that differ markedly among protein conformations. They further indicated state-dependent effects on membrane shape, with the E2 state being likely disfavored in cholesterol-rich bilayers relative to the E1P state because of a greater hydrophobic mismatch. In summary, cholesterol extraction from membranes significantly decreases Na+,K+-ATPase steady-state activity.

Rationalising pKa shifts in Bacillus circulans xylanase with computational studies

Bacillus circulans xylanase (BcX), a family 11 glycoside hydrolase, catalyses the hydrolysis of xylose polymers with a net retention of stereochemistry. Glu172 in BcX is believed to act as a general acid by protonating the aglycone during glycosylation, and then as a general base to facilitate the deglycosylation step. The key to the dual role of this general acid/base lies in its protonation states, which depend on its intrinsic pK_a value and the specific environment which it resides within. To fully understand the detailed molecular features in BcX to establish the dual role of Glu172, we present a combined study based on both atomistic simulations and empirical models to calculate pK_a shifts for the general acid/base Glu172 in BcX at different functional states. Its pK_a values and those of nearby residues, obtained based on QM/MM free energy calculations, MCCE and PROPKA, show a good agreement with available experimental data. Additionally, our study provides additional insights into the effects of structural and electrostatic perturbations caused by mutations and chemical modifications, suggesting that the local solvation environment and mutagenesis of the residues adjacent to Glu172 establish its dual role during hydrolysis. The strengths and limitations of various methods for calculating pK_as and pK_a shifts have also been discussed.

The CAPOS mutation in ATP1A3 alters Na/K-ATPase function and results in auditory neuropathy which has implications for management

Cerebellar ataxia, areflexia, pes cavus, optic atrophy and sensorineural hearing impairment (CAPOS) is a rare clinically distinct syndrome caused by a single dominant missense mutation, c.2452G>A, p.Glu818Lys, in ATP1A3, encoding the neuron-specific alpha subunit of the Na+/K+-ATPase α3. Allelic mutations cause the neurological diseases rapid dystonia Parkinsonism and alternating hemiplegia of childhood, disorders which do not encompass hearing or visual impairment. We present detailed clinical phenotypic information in 18 genetically confirmed patients from 11 families (10 previously unreported) from Denmark, Sweden, UK and Germany indicating a specific type of hearing impairment-auditory neuropathy (AN). All patients were clinically suspected of CAPOS and had hearing problems. In this retrospective analysis of audiological data, we show for the first time that cochlear outer hair cell activity was preserved as shown by the presence of otoacoustic emissions and cochlear microphonic potentials, but the auditory brainstem responses were grossly abnormal, likely reflecting neural dyssynchrony. Poor speech perception was observed, especially in noise, which was beyond the hearing level obtained in the pure tone audiograms in several of the patients presented here. Molecular modelling and in vitro electrophysiological studies of the specific CAPOS mutation were performed. Heterologous expression studies of α3 with the p.Glu818Lys mutation affects sodium binding to, and release from, the sodium-specific site in the pump, the third ion-binding site. Molecular dynamics simulations confirm that the structure of the C-terminal region is affected. In conclusion, we demonstrate for the first time evidence for auditory neuropathy in CAPOS syndrome, which may reflect impaired propagation of electrical impulses along the spiral ganglion neurons. This has implications for diagnosis and patient management. Auditory neuropathy is difficult to treat with conventional hearing aids, but preliminary improvement in speech perception in some patients suggests that cochlear implantation may be effective in CAPOS patients.

Distinct pH dependencies of Na+/K+ selectivity at the two faces of Na,K-ATPase

The sodium pump (Na,K-ATPase) in animal cells is vital for actively maintaining ATP hydrolysis-powered Na⁺ and K⁺ electrochemical gradients across the cell membrane. These ion gradients drive co- and countertransport and are critical for establishing the membrane potential. It has been an enigma how Na,K-ATPase discriminates between Na⁺ and K⁺, despite the pumped ion on each side being at a lower concentration than the other ion. Recent crystal structures of analogs of the intermediate conformations E2·Pi·2K⁺ and Na⁺-bound E1∼P·ADP suggest that the dimensions of the respective binding sites in Na,K-ATPase are crucial in determining its selectivity. Here, we found that the selectivity at each membrane face is pH-dependent and that this dependence is unique for each face. Most notable was a strong increase in the specific affinity for K⁺ at the extracellular face (i.e. E2 conformation) as the pH is lowered from 7.5 to 5. We also observed a smaller increase in affinity for K⁺ on the cytoplasmic side (E1 conformation), which reduced the selectivity for Na⁺ Theoretical analysis of the pK_a values of ion-coordinating acidic amino acid residues suggested that the face-specific pH dependences and Na⁺/K⁺ selectivities may arise from the protonation or ionization of key residues. The increase in K⁺ selectivity at low pH on the cytoplasmic face, for instance, appeared to be associated with Asp⁸⁰⁸ protonation. We conclude that changes in the ionization state of coordinating residues in Na,K-ATPase could contribute to altering face-specific ion selectivity.

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.

Molecular simulations and free-energy calculations suggest conformation-dependent anion binding to a cytoplasmic site as a mechanism for Na+/K+-ATPase ion selectivity

Na⁺/K⁺-ATPase transports Na⁺ and K⁺ ions across the cell membrane via an ion-binding site becoming alternatively accessible to the intra- and extracellular milieu by conformational transitions that confer marked changes in ion-binding stoichiometry and selectivity. To probe the mechanism of these changes, we used molecular simulation and free-energy perturbation approaches to identify probable protonation states of Na⁺- and K⁺-coordinating residues in E1P and E2P conformations of Na⁺/K⁺-ATPase. Analysis of these simulations revealed a molecular mechanism responsible for the change in protonation state: the conformation-dependent binding of an anion (a chloride ion in our simulations) to a previously unrecognized cytoplasmic site in the loop between transmembrane helices 8 and 9, which influences the electrostatic potential of the crucial Na⁺-coordinating residue Asp⁹²⁶ This mechanistic model is consistent with experimental observations and provides a molecular-level picture of how E1P to E2P enzyme conformational transitions are coupled to changes in ion-binding stoichiometry and selectivity.

Ion channels and ion selectivity

Specific macromolecular transport systems, ion channels and pumps, provide the pathways to facilitate and control the passage of ions across the lipid membrane. Ion channels provide energetically favourable passage for ions to diffuse rapidly and passively according to their electrochemical potential. Selective ion channels are essential for the excitability of biological membranes: the action potential is a transient phenomenon that reflects the rapid opening and closing of voltage-dependent Na⁺-selective and K⁺-selective channels. One of the most critical functional aspects of K⁺ channels is their ability to remain highly selective for K⁺ over Na⁺ while allowing high-throughput ion conduction at a rate close to the diffusion limit. Permeation through the K⁺ channel selectivity filter is believed to proceed as a 'knockon' mechanism, in which 2-3 K⁺ ions interspersed by water molecules move in a single file. Permeation through the comparatively wider and less selective Na⁺ channels also proceeds via a loosely coupled knockon mechanism, although the ions do not need to be fully dehydrated. While simple structural concepts are often invoked to rationalize the mechanism of ion selectivity, a deeper analysis shows that subtle effects play an important role in these flexible dynamical structures.

Glutamate Water Gates in the Ion Binding Pocket of Na+ Bound Na+, K+-ATPase

The dynamically changing protonation states of the six acidic amino acid residues in the ion binding pocket of the Na+, K+ -ATPase (NKA) during the ion transport cycle are proposed to drive ion binding, release and possibly determine Na+ or K+ selectivity. We use molecular dynamics (MD) and density functional theory (DFT) simulations to determine the protonation scheme of the Na+ bound conformation of NKA. MD simulations of all possible protonation schemes show that the bound Na+ ions are most stably bound when three or four protons reside in the binding sites, and that Glu954 in site III is always protonated. Glutamic acid residues in the three binding sites act as water gates, and their deprotonation triggers water entry to the binding sites. From DFT calculations of Na+ binding energies, we conclude that three protons in the binding site are needed to effectively bind Na+ from water and four are needed to release them in the next step. Protonation of Asp926 in site III will induce Na+ release, and Glu327, Glu954 and Glu779 are all likely to be protonated in the Na+ bound occluded conformation. Our data provides key insights into the role of protons in the Na+ binding and release mechanism of NKA.

Ion Pathways in the Na+/K+-ATPase

Na⁺/K⁺-ATPase (NKA) is an essential cation pump protein responsible for the maintenance of the sodium and potassium gradients across the plasma membrane. Recently published high-resolution structures revealed amino acids forming the cation binding sites (CBS) in the transmembrane domain and variable position of the domains in the cytoplasmic headpiece. Here we report molecular dynamic simulations of the human NKA α1β1 isoform embedded into DOPC bilayer. We have analyzed the NKA conformational changes in the presence of Na⁺- or K⁺-cations in the CBS, for various combinations of the cytoplasmic ligands, and the two major enzyme conformations in the 100 ns runs (more than 2.5 μs of simulations in total). We identified two novel cytoplasmic pathways along the pairs of transmembrane helices TM3/TM7 or TM6/TM9 that allow hydration of the CBS or transport of cations from/to the bulk. These findings can provide a structural explanation for previous mutagenesis studies, where mutation of residues that are distal from the CBS resulted in the alteration of the enzyme affinity to the transported cations or change in the enzyme activity.

The selectivity of the Na(+)/K(+)-pump is controlled by binding site protonation and self-correcting occlusion

The Na⁺/K⁺-pump maintains the physiological K⁺ and Na⁺ electrochemical gradients across the cell membrane. It operates via an 'alternating-access' mechanism, making iterative transitions between inward-facing (E₁) and outward-facing (E₂) conformations. Although the general features of the transport cycle are known, the detailed physicochemical factors governing the binding site selectivity remain mysterious. Free energy molecular dynamics simulations show that the ion binding sites switch their binding specificity in E₁ and E₂. This is accompanied by small structural arrangements and changes in protonation states of the coordinating residues. Additional computations on structural models of the intermediate states along the conformational transition pathway reveal that the free energy barrier toward the occlusion step is considerably increased when the wrong type of ion is loaded into the binding pocket, prohibiting the pump cycle from proceeding forward. This self-correcting mechanism strengthens the overall transport selectivity and protects the stoichiometry of the pump cycle.

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

A protein called the sodium-potassium pump resides in the membrane that surrounds living cells. The role of this protein is to 'pump' sodium and potassium ions across the membrane to help restore their concentration inside and outside of the cell. About 25% of the body's energy is used to keep the pump going, rising to nearly 70% in the brain. Problems that affect the pump have been linked to several disorders, including heart, kidney and metabolic diseases, as well as severe neurological conditions.

The sodium-potassium pump must be able to effectively pick out the correct ions to transport from a mixture of many different types of ions. However, it was not clear how the pump succeeds in doing this efficiently. Rui et al. have now used a computational method called molecular dynamics simulations to model how the sodium-potassium pump transports the desired ions across the cell membrane.

The pump works via a so-called 'alternating-access' mechanism, repeatedly transitioning between inward-facing and outward-facing conformations. In each cycle, it binds three sodium ions from the cell’s interior and exports them to the outside. Then, the pump binds to two potassium ions from outside the cell and imports them inside. Although the bound sodium and potassium ions interact with similar binding sites in the pump, the pump sometimes preferentially binds sodium, and sometimes potassium. The study performed by Rui et al. shows that this preference is driven by how protons (hydrogen ions) bind to the amino acids that make up the binding site.

The simulations also suggest that the pump uses a ‘self-correcting’ mechanism to prevent the pump from transporting the wrong types of ions. When incorrect ions are present at the binding sites, the pump cycle pauses temporarily until the ions detach from the pump. Only when the correct ions are bound will the pump cycle continue again.

In the future, Rui et al. hope to use long time-scale molecular dynamics simulations to show the conformational transition in action. In addition, the 'self-correcting' mechanism will be directly tested by letting the wrong and correct ions compete for the binding sites to see whether the pump will transport only the correct ions.

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

Atomic-level mechanisms for phospholamban regulation of the calcium pump

We performed protein pKa calculations and molecular dynamics (MD) simulations of the calcium pump (sarcoplasmic reticulum Ca(2+)-ATPase (SERCA)) in complex with phospholamban (PLB). X-ray crystallography studies have suggested that PLB locks SERCA in a low-Ca(2+)-affinity E2 state that is incompatible with metal-ion binding, thereby blocking the conversion toward a high-Ca(2+)-affinity E1 state. Estimation of pKa values of the acidic residues in the transport sites indicates that at normal intracellular pH (7.1-7.2), PLB-bound SERCA populates an E1 state that is deprotonated at residues E309 and D800 yet protonated at residue E771. We performed three independent microsecond-long MD simulations to evaluate the structural dynamics of SERCA-PLB in a solution containing 100 mM K(+) and 3 mM Mg(2+). Principal component analysis showed that PLB-bound SERCA lies exclusively along the structural ensemble of the E1 state. We found that the transport sites of PLB-bound SERCA are completely exposed to the cytosol and that K(+) ions bind transiently (≤5 ns) and nonspecifically (nine different positions) to the two transport sites, with a total occupancy time of K(+) in the transport sites of 80%. We propose that PLB binding to SERCA populates a novel (to our knowledge) E1 intermediate, E1⋅H(+)771. This intermediate serves as a kinetic trap that controls headpiece dynamics and depresses the structural transitions necessary for Ca(2+)-dependent activation of SERCA. We conclude that PLB-mediated regulation of SERCA activity in the heart results from biochemical and structural transitions that occur primarily in the E1 state of the pump.

MD Simulations of P-Type ATPases in a Lipid Bilayer System

Molecular dynamics (MD) simulation is a computational method which provides insight on protein dynamics with high resolution in both space and time, in contrast to many experimental techniques. MD simulations can be used as a stand-alone method to study P-type ATPases as well as a complementary method aiding experimental studies. In particular, MD simulations have proved valuable in generating and confirming hypotheses relating to the structure and function of P-type ATPases. In the following, we describe a detailed practical procedure on how to set up and run a MD simulation of a P-type ATPase embedded in a lipid bilayer using software free of use for academics. We emphasize general considerations and problems typically encountered when setting up simulations. While full coverage of all possible procedures is beyond the scope of this chapter, we have chosen to illustrate the MD procedure with the Nanoscale Molecular Dynamics (NAMD) and the Visual Molecular Dynamics (VMD) software suites.

Atomistic Characterization of the First Step of Calcium Pump Activation Associated with Proton Countertransport

The calcium pump [sarcoplasmic reticulum (SR) Ca(2+)-ATPase (SERCA)] transports Ca(2+) from the cytosol to the SR lumen at the expense of ATP hydrolysis and proton countertransport, thus playing a central role in Ca(2+) homeostasis and muscle contractility. Proton countertransport via deprotonation of transport site residue Glu309 is a critical first step in SERCA activation because it accelerates the E2-E1 structural transition. Previous studies have suggested that flipping of Glu309 toward the cytosol constitutes the primary mechanism for Glu309 deprotonation, but no conclusive data to support this hypothesis have been published. Therefore, we performed three independent 1 μs molecular dynamics simulations of the E2 state protonated at transport site residues Glu309, Glu771, and Glu908. Structural analysis and pKa calculations showed that Glu309 deprotonation occurs by an inward-to-outward side-chain transition. We also found that Glu309 deprotonation and proton countertransport occur through transient (~113 ps) water wires connecting Glu309 with the cytosol. Although both mechanisms are operational, we found that transient water wire formation, and not Glu309 flipping, is the primary mechanism for Glu309 deprotonation and translocation of protons to the cytosol. The outward-to-inward transition of protonated Glu309 and the presence of water wires suggest that protons from the cytosol might be passively transported to the lumen via Glu309. However, structural analysis indicates that passive SR proton leakage into the lumen unlikely occurs through Glu309 in the E2 state. These findings provide a time-resolved visualization of the first step in the molecular mechanism of SERCA activation and proton transport across the SR.

Mechanism of potassium ion uptake by the Na(+)/K(+)-ATPase

The Na⁺/K⁺-ATPase restores sodium (Na⁺) and potassium (K⁺) electrochemical gradients dissipated by action potentials and ion-coupled transport processes. As ions are transported, they become transiently trapped between intracellular and extracellular gates. Once the external gate opens, three Na⁺ ions are released, followed by the binding and occlusion of two K⁺ ions. While the mechanisms of Na⁺ release have been well characterized by the study of transient Na⁺ currents, smaller and faster transient currents mediated by external K⁺ have been more difficult to study. Here we show that external K⁺ ions travelling to their binding sites sense only a small fraction of the electric field as they rapidly and simultaneously become occluded. Consistent with these results, molecular dynamics simulations of a pump model show a wide water-filled access channel connecting the binding site to the external solution. These results suggest a mechanism of K⁺ gating different from that of Na⁺ occlusion.

During transport by the Na⁺/K⁺-ATPase, Na⁺ and K⁺ ions become occluded between intra- and extracellular gates. Here Castillo et al. measure transient electrical signals arising from K⁺ occlusion and use molecular simulations to describe a K⁺ gating mechanism fundamentally different to that of Na⁺.

Sodium and proton effects on inward proton transport through Na/K pumps

The Na/K pump hydrolyzes ATP to export three intracellular Na (Nai) as it imports two extracellular K (Ko) across animal plasma membranes. Within the protein, two ion-binding sites (sites I and II) can reciprocally bind Na or K, but a third site (site III) exclusively binds Na in a voltage-dependent fashion. In the absence of Nao and Ko, the pump passively imports protons, generating an inward current (IH). To elucidate the mechanisms of IH, we used voltage-clamp techniques to investigate the [H]o, [Na]o, and voltage dependence of IH in Na/K pumps from ventricular myocytes and in ouabain-resistant pumps expressed in Xenopus oocytes. Lowering pHo revealed that Ho both activates IH (in a voltage-dependent manner) and inhibits it (in a voltage-independent manner) by binding to different sites. Nao effects depend on pHo; at pHo where no Ho inhibition is observed, Nao inhibits IH at all concentrations, but when applied at pHo that inhibits pump-mediated current, low [Na]o activates IH and high [Na]o inhibits it. Our results demonstrate that IH is a property inherent to Na/K pumps, not linked to the oocyte expression environment, explains differences in the characteristics of IH previously reported in the literature, and supports a model in which 1), protons leak through site III; 2), binding of two Na or two protons to sites I and II inhibits proton transport; and 3), pumps with mixed Na/proton occupancy of sites I and II remain permeable to protons.

K+ congeners that do not compromise Na+ activation of the Na+,K+-ATPase: hydration of the ion binding cavity likely controls ion selectivity

The Na(+),K(+)-ATPase is essential for ionic homeostasis in animal cells. The dephosphoenzyme contains Na(+) selective inward facing sites, whereas the phosphoenzyme contains K(+) selective outward facing sites. Under normal physiological conditions, K(+) inhibits cytoplasmic Na(+) activation of the enzyme. Acetamidinium (Acet(+)) and formamidinium (Form(+)) have been shown to permeate the pump through the outward facing sites. Here, we show that these cations, unlike K(+), are unable to enter the inward facing sites in the dephosphorylated enzyme. Consistently, the organic cations exhibited little to no antagonism to cytoplasmic Na(+) activation. Na(+),K(+)-ATPase structures revealed a previously undescribed rotamer transition of the hydroxymethyl side chain of the absolutely conserved Thr(772) of the α-subunit. The side chain contributes its hydroxyl to Na(+) in site I in the E1 form and rotates to contribute its methyl group toward K(+) in the E2 form. Molecular dynamics simulations to the E1·AlF4 (-)·ADP·3Na(+) structure indicated that 1) bound organic cations differentially distorted the ion binding sites, 2) the hydroxymethyl of Thr(772) rotates to stabilize bound Form(+) through water molecules, and 3) the rotamer transition is mediated by water traffic into the ion binding cavity. Accordingly, dehydration induced by osmotic stress enhanced the interaction of the congeners with the outward facing sites and profoundly modified the organization of membrane domains of the α-subunit. These results assign a catalytic role for water in pump function, and shed light on a backbone-independent but a conformation-dependent switch between H-bond and dispersion contact as part of the catalytic mechanism of the Na(+),K(+)-ATPase.

Permeation redux: thermodynamics and kinetics of ion movement through potassium channels

The fundamental biophysics underlying the selective movement of ions through ion channels was launched by George Eisenman in the 1960s, using glass electrodes. This minireview examines the insights from these early studies and the explosive progress made since then.

Route, mechanism, and implications of proton import during Na+/K+ exchange by native Na+/K+-ATPase pumps

The Na⁺/K⁺ pump is a hybrid transporter that can also import protons at physiological K⁺ and Na⁺ concentrations.

A single Na⁺/K⁺-ATPase pumps three Na⁺ outwards and two K⁺ inwards by alternately exposing ion-binding sites to opposite sides of the membrane in a conformational sequence coupled to pump autophosphorylation from ATP and auto-dephosphorylation. The larger flow of Na⁺ than K⁺ generates outward current across the cell membrane. Less well understood is the ability of Na⁺/K⁺ pumps to generate an inward current of protons. Originally noted in pumps deprived of external K⁺ and Na⁺ ions, as inward current at negative membrane potentials that becomes amplified when external pH is lowered, this proton current is generally viewed as an artifact of those unnatural conditions. We demonstrate here that this inward current also flows at physiological K⁺ and Na⁺ concentrations. We show that protons exploit ready reversibility of conformational changes associated with extracellular Na⁺ release from phosphorylated Na⁺/K⁺ pumps. Reversal of a subset of these transitions allows an extracellular proton to bind an acidic side chain and to be subsequently released to the cytoplasm. This back-step of phosphorylated Na⁺/K⁺ pumps that enables proton import is not required for completion of the 3 Na⁺/2 K⁺ transport cycle. However, the back-step occurs readily during Na⁺/K⁺ transport when external K⁺ ion binding and occlusion are delayed, and it occurs more frequently when lowered extracellular pH raises the probability of protonation of the externally accessible carboxylate side chain. The proton route passes through the Na⁺-selective binding site III and is distinct from the principal pathway traversed by the majority of transported Na⁺ and K⁺ ions that passes through binding site II. The inferred occurrence of Na⁺/K⁺ exchange and H⁺ import during the same conformational cycle of a single molecule identifies the Na⁺/K⁺ pump as a hybrid transporter. Whether Na⁺/K⁺ pump–mediated proton inflow may have any physiological or pathophysiological significance remains to be clarified.

ATP binding and aspartate protonation enhance photoinduced electron transfer in plant cryptochrome

Cryptochromes are flavoproteins encountered in most vegetal and animal species. They play a role of blue-light receptors in plants and in invertebrates. The putative resting state of the FAD cofactor in these proteins is its fully oxidized form, FADox. Upon blue-light excitation, the isoalloxazine ring (ISO) may undergo an ultrafast reduction by a nearby tryptophan residue W400. This primary reduction triggers a cascade of electron and proton transfers, ultimately leading to the formation of the FADH° radical. A recent experimental study has shown that the yield of FADH° formation in Arabidopsis cryptochrome can be strongly modulated by ATP binding and by pH, affecting the protonation state of D396 (proton donor to FAD°(-)). Here we provide a detailed molecular analysis of these effects by means of combined classical molecular dynamics simulations and time-dependent density functional theory calculations. When ATP is present and D396 protonated, FAD remains in close contact with W400, thereby enhancing electron transfer (ET) from W400 to ISO*. In contrast, deprotonation of D396 and absence of ATP introduce flexibility to the photoactive site prior to FAD excitation, with the consequence of increased ISO-W400 distance and diminished tunneling rate by almost two orders of magnitude. We show that under these conditions, ET from the adenine moiety of FAD becomes a competitive relaxation pathway. Overall, our data suggest that the observed effects of ATP and pH on the FAD photoreduction find their roots in the earliest stage of the photoreduction process; i.e., ATP binding and the protonation state of D396 determine the preferred pathway of ISO* relaxation.

Inhibition of K+ transport through Na+, K+-ATPase by capsazepine: role of membrane span 10 of the α-subunit in the modulation of ion gating

Capsazepine (CPZ) inhibits Na⁺,K⁺-ATPase-mediated K⁺-dependent ATP hydrolysis with no effect on Na⁺-ATPase activity. In this study we have investigated the functional effects of CPZ on Na⁺,K⁺-ATPase in intact cells. We have also used well established biochemical and biophysical techniques to understand how CPZ modifies the catalytic subunit of Na⁺,K⁺-ATPase. In isolated rat cardiomyocytes, CPZ abolished Na⁺,K⁺-ATPase current in the presence of extracellular K⁺. In contrast, CPZ stimulated pump current in the absence of extracellular K⁺. Similar conclusions were attained using HEK293 cells loaded with the Na⁺ sensitive dye Asante NaTRIUM green. Proteolytic cleavage of pig kidney Na⁺,K⁺-ATPase indicated that CPZ stabilizes ion interaction with the K⁺ sites. The distal part of membrane span 10 (M10) of the α-subunit was exposed to trypsin cleavage in the presence of guanidinum ions, which function as Na⁺ congener at the Na⁺ specific site. This effect of guanidinium was amplified by treatment with CPZ. Fluorescence of the membrane potential sensitive dye, oxonol VI, was measured following addition of substrates to reconstituted inside-out Na⁺,K⁺-ATPase. CPZ increased oxonol VI fluorescence in the absence of K⁺, reflecting increased Na⁺ efflux through the pump. Surprisingly, CPZ induced an ATP-independent increase in fluorescence in the presence of high extravesicular K⁺, likely indicating opening of an intracellular pathway selective for K⁺. As revealed by the recent crystal structure of the E₁.AlF₄^-.ADP.3Na⁺ form of the pig kidney Na⁺,K⁺-ATPase, movements of M5 of the α-subunit, which regulate ion selectivity, are controlled by the C-terminal tail that extends from M10. We propose that movements of M10 and its cytoplasmic extension is affected by CPZ, thereby regulating ion selectivity and transport through the K⁺ sites in Na⁺,K⁺-ATPase.

Microsecond molecular dynamics simulations of Mg²⁺- and K⁺-bound E1 intermediate states of the calcium pump

We have performed microsecond molecular dynamics (MD) simulations to characterize the structural dynamics of cation-bound E1 intermediate states of the calcium pump (sarcoendoplasmic reticulum Ca²⁺-ATPase, SERCA) in atomic detail, including a lipid bilayer with aqueous solution on both sides. X-ray crystallography with 40 mM Mg²⁺ in the absence of Ca²⁺ has shown that SERCA adopts an E1 structure with transmembrane Ca²⁺-binding sites I and II exposed to the cytosol, stabilized by a single Mg²⁺ bound to a hybrid binding site I'. This Mg²⁺-bound E1 intermediate state, designated E1•Mg²⁺, is proposed to constitute a functional SERCA intermediate that catalyzes the transition from E2 to E1•2Ca²⁺ by facilitating H⁺/Ca²⁺ exchange. To test this hypothesis, we performed two independent MD simulations based on the E1•Mg²⁺ crystal structure, starting in the presence or absence of initially-bound Mg²⁺. Both simulations were performed for 1 µs in a solution containing 100 mM K⁺ and 5 mM Mg²⁺ in the absence of Ca²⁺, mimicking muscle cytosol during relaxation. In the presence of initially-bound Mg²⁺, SERCA site I' maintained Mg²⁺ binding during the entire MD trajectory, and the cytosolic headpiece maintained a semi-open structure. In the absence of initially-bound Mg²⁺, two K⁺ ions rapidly bound to sites I and I' and stayed loosely bound during most of the simulation, while the cytosolic headpiece shifted gradually to a more open structure. Thus MD simulations predict that both E1•Mg²⁺ and E•2K+ intermediate states of SERCA are populated in solution in the absence of Ca²⁺, with the more open 2K+-bound state being more abundant at physiological ion concentrations. We propose that the E1•2K⁺ state acts as a functional intermediate that facilitates the E2 to E1•2Ca²⁺ transition through two mechanisms: by pre-organizing transport sites for Ca²⁺ binding, and by partially opening the cytosolic headpiece prior to Ca²⁺ activation of nucleotide binding.

Molecular characterization of the α-subunit of Na⁺/K⁺ ATPase from the euryhaline barnacle Balanus improvisus reveals multiple genes and differential expression of alternative splice variants

The euryhaline bay barnacle Balanus improvisus has one of the broadest salinity tolerances of any barnacle species. It is able to complete its life cycle in salinities close to freshwater (3 PSU) up to fully marine conditions (35 PSU) and is regarded as one of few truly brackish-water species. Na⁺/K⁺ ATPase (NAK) has been shown to be important for osmoregulation when marine organisms are challenged by changing salinities, and we therefore cloned and examined the expression of different NAKs from B. improvisus. We found two main gene variants, NAK1 and NAK2, which were approximately 70% identical at the protein level. The NAK1 mRNA existed in a long and short variant with the encoded proteins differing only by 27 N-terminal amino acids. This N-terminal stretch was coded for by a separate exon, and the two variants of NAK1 mRNAs appeared to be created by alternative splicing. We furthermore showed that the two NAK1 isoforms were differentially expressed in different life stages and in various tissues of adult barnacle, i.e the long isoform was predominant in cyprids and in adult cirri. In barnacle cyprid larvae that were exposed to a combination of different salinities and pCO₂ levels, the expression of the long NAK1 mRNA increased relative to the short in low salinities. We suggest that the alternatively spliced long variant of the Nak1 protein might be of importance for osmoregulation in B. improvisus in low salinity conditions.

Relative Free Energies for Hydration of Monovalent Ions from QM and QM/MM Simulations

Methods directly evaluating the hydration structure and thermodynamics of physiologically relevant cations (Na(+), K(+), Cl(-), etc.) have wide ranging applications in the fields of inorganic, physical, and biological chemistry. All-atom simulations based on accurate potential energy surfaces appear to offer a viable option for assessing the chemistry of ion solvation. Although MD and free energy simulations of ion solvation with classical force fields have proven their usefulness, a number of challenges still remain. One of them is the difficulty of force field benchmarking and validation against structural and thermodynamic data obtained for a condensed phase. Hybrid quantum mechanical/molecular mechanical (QM/MM) models combined with sampling algorithms have the potential to provide an accurate solvation model and to incorporate the effects from the surrounding, which is often missing in gas-phase ab initio computations. Herein, we report the results from QM/MM free energy simulations of Na(+)/K(+) and Cl(-)/Br(-) hydration where we simultaneously characterized the relative thermodynamics of ion solvation and changes in the solvation structure. The Flexible Inner Region Ensemble Separator (FIRES) method was used to impose a spatial separation between QM region and the outer sphere of solvent molecules treated with the CHARMM27 force field. FEP calculations based on QM/MM simulations utilizing the CHARMM/deMon2k interface were performed with different basis set combinations for K(+)/Na(+) and Cl(-)/Br(-) perturbations to establish the dependence of the computed free energies on the basis set level. The dependence of the computed relative free energies on the size of the QM and MM regions is discussed. The current methodology offers an accurate description of structural and thermodynamic aspects of the hydration of alkali and halide ions in neat solvents and can be used to obtain thermodynamic data on ion solvation in condensed phase along with underlying structural properties of the ion-solvent system.

New crystal structures of PII-type ATPases: excitement continues

P-type ATPases are ATP-powered ion pumps, classified into five subfamilies (PI-PV). Of these, PII-type ATPases, including Ca2+-ATPase, Na+,K+-ATPase and gastric H+,K+-ATPase, among others, have been the most intensively studied. Best understood structurally and biochemically is Ca2+-ATPase from sarcoplasmic reticulum of fast twitch skeletal muscle (sarco(endo)plasmic reticulum Ca2+-ATPase 1a, SERCA1a). Since publication of the first crystal structure in 2000, it has continuously been a source of excitement, as crystal structures for new reaction intermediates always show large structural changes. Crystal structures now exist for most of the reaction intermediates, almost covering the entire reaction cycle. This year the crystal structure of a missing link, the E1·Mg2+ state, finally appeared, bringing another surprise: bound sarcolipin (SLN). The current status of two other important PII-type ATPases, Na+,K+-ATPase and H+,K+-ATPase, is also briefly described.

Alternative cycling modes of the Na(+)/K(+)-ATPase in the presence of either Na(+) or Rb(+)

A comprehensive study of the interaction between Na(+) and K(+) with the Na(+)/K(+)-ATPase requires dissecting the incidence of alternative cycling modes on activity measurements in which one or both of these cations are absent. With this aim, we used membrane fragments containing pig-kidney Na(+)/K(+)-ATPase to perform measurements, at 25°C and pH=7.4, of ATPase activity and steady-state levels of (i) intermediates containing occluded Rb(+) at different [Rb(+)] in media lacking Na(+), and (ii) phosphorylated intermediates at different [Na(+)] in media lacking Rb(+). Most relevant results are: (1) Rb(+) can be occluded through an ATPasic cycling mode that takes place in the absence of Na(+) ions, (2) the kinetic behavior of the phosphoenzyme formed by ATP in the absence of Na(+) is different from the one that is formed with Na(+), and (3) binding of Na(+) to transport sites during catalysis is not at random unless rapid equilibrium holds.

Protonation state of a conserved acidic amino acid involved in Na(+) binding to the glutamate transporter EAAC1

Substrate transport by glutamate transporters is coupled to the co-transport of 3 Na(+) ions and counter-transport of 1 K(+) ion. The highly conserved Asp454, which may be negatively charged, is of interest as its side chain may coordinate cations and/or contribute to charge compensation. Mutation to the nonionizable Asn resulted in a transporter that no longer catalyzed forward transport. However, Na(+)/glutamate exchange was still functional, as demonstrated by the presence of transient currents following rapid substrate application and voltage jumps. While the kinetics of Na(+)/glutamate exchange were slowed, the apparent valence (z) of the charge moved in EAAC1 D454N (0.71) was similar to that of EAAC1 WT (0.64). Valences calculated using the Poisson-Boltzmann equation were close to the experimental values for EAAC1 D454N (0.55), and with D454 protonated (0.45). In addition, pK(a) calculations performed for the bacterial homologue GltPh revealed a highly perturbed pK(a) (7.6 to >14) for D405 residue (analogous to D454), consistent with this site being protonated at physiological pH. In contrast to the D454N mutation, substitution to alanine resulted in a transporter that still bound glutamate, but could not translocate it. The results are consistent with molecular dynamics simulations, showing that the alanine but not the asparagine mutation resulted in defective Na(+) coordination. Our results raise the possibility that the protonated state of D454 supports transporter function.

Ion binding sites and their representations by reduced models

The binding of small metal ions to complex macromolecular structures is typically dominated by strong local interactions of the ion with its nearest ligands. Progress in understanding the molecular determinants of ion selectivity can often be achieved by considering simplified reduced models comprised of only the most important ion-coordinating ligands. Although the main ingredients underlying simplified reduced models are intuitively clear, a formal statistical mechanical treatment is nonetheless necessary in order to draw meaningful conclusions about complex macromolecular systems. By construction, reduced models only treat the ion and the nearest coordinating ligands explicitly. The influence of the missing atoms from the protein or the solvent is incorporated indirectly. Quasi-chemical theory offers one example of how to carry out such a separation in the case of ion solvation in bulk liquids, and in several ways, a statistical mechanical formulation of reduced binding site models for macromolecules is expected to follow a similar route. However, there are also important differences when the ion-coordinating moieties are not solvent molecules from a bulk phase but are molecular ligands covalently bonded to a macromolecular structure. Here, a statistical mechanical formulation of reduced binding site models is elaborated to address these issues. The formulation provides a useful framework to construct reduced binding site models, and define the average effect from the surroundings on the ion and the nearest coordinating ligands.

Extended protein/water H-bond networks in photosynthetic water oxidation

Oxidation of water molecules in the photosystem II (PSII) protein complex proceeds at the manganese-calcium complex, which is buried deeply in the lumenal part of PSII. Understanding the PSII function requires knowledge of the intricate coupling between the water-oxidation chemistry and the dynamic proton management by the PSII protein matrix. Here we assess the structural basis for long-distance proton transfer in the interior of PSII and for proton management at its surface. Using the recent high-resolution crystal structure of PSII, we investigate prominent hydrogen-bonded networks of the lumenal side of PSII. This analysis leads to the identification of clusters of polar groups and hydrogen-bonded networks consisting of amino acid residues and water molecules. We suggest that long-distance proton transfer and conformational coupling is facilitated by hydrogen-bonded networks that often involve more than one protein subunit. Proton-storing Asp/Glu dyads, such as the D1-E65/D2-E312 dyad connected to a complex water-wire network, may be particularly important for coupling protonation states to the protein conformation. Clusters of carboxylic amino acids could participate in proton management at the lumenal surface of PSII. We propose that rather than having a classical hydrophobic protein interior, the lumenal side of PSII resembles a complex polyelectrolyte with evolutionary optimized hydrogen-bonding networks. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.

Molecular Mechanisms of K+ Selectivity in Na/K Pump

The sodium–potassium (Na/K) pump plays an essential role in maintaining cell volume and secondary active transport of other solutes by establishing the Na+ and K+ concentration gradients across the plasma membrane of animal cells. The recently determined crystal structures of the Na/K pump to atomic resolution provide a new impetus to investigate molecular determinants governing the binding of Na+ and K+ ions and conformational transitions during the functional cycle. The pump cycle is generally described by the alternating access mechanism, in which the pump toggles between different conformational states, where ions can bind from either the intracellular or the extracellular side. However, important issues concerning the selectivity of the Na/K pump remain to be addressed. In particular, two out of the three binding sites are shared between Na+ and K+ and it is not clear how the protein is able to select K+ over Na+ when it is in the outwardly facing phosphorylated conformation (E2P), and Na+ over K+ when it is in the inwardly facing conformation (E1). In this review article, we will first briefly review the recent advancement in understanding the microscopic mechanism of K+ selectivity in the Na/K pump at the E2·Pi state and then outline the remaining challenges to be addressed about ion selectivity.

Protonation states of important acidic residues in the central Ca²⁺ ion binding sites of the Ca²⁺-ATPase: a molecular modeling study

The P-type ATPases are responsible for the transport of cations across cell membranes. The sarco(endo)plasmic reticulum Ca²⁺-ATPase (SERCA) transports two Ca²⁺ ions from the cytoplasm to the lumen of the sarco(endo)plasmic reticulum and countertransports two or three protons per catalytic cycle. Two binding sites for Ca²⁺ ions have been located via protein crystallography, including four acidic amino acid residues that are essential to the ion coordination. In this study, we present molecular dynamics (MD) simulations examining the protonation states of these amino acid residues in a Ca²⁺-free conformation of SERCA. Such knowledge will be important for an improved understanding of atomistic details of the transport mechanism of protons and Ca²⁺ ions. Eight combinations of the protonation of four central acidic residues, Glu309, Glu771, Asp800, and Glu908, are tested from 10 ns MD simulations with respect to protein stability and ability to maintain a structure similar to the crystal structure. The trajectories for the most prospective combinations of protonation states were elongated to 50 ns and subjected to more detailed analysis, including prediction of pK(a) values of the four acidic residues over the trajectories. From the simulations we find that the combination leaving only Asp800 as charged is most likely. The results are compared to available experimental data and explain the observed destabilization upon full deprotonation, resulting in the entry of cytoplasmic K⁺ ions into the Ca²⁺ binding sites during the simulation in which Ca²⁺ ions are absent. Furthermore, a hypothesis for the exchange of protons from the central binding cavity is proposed.