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

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

Dimerization mechanism of an inverted-topology ion channel in membranes.

Many ion channels are multisubunit complexes where oligomerization is an obligatory requirement for function as the binding axis forms the charged permeation pathway. However, the mechanisms of in-membrane assembly of thermodynamically stable channels are largely unknown. Here, we demonstrate a key advance by reporting the dimerization equilibrium reaction of an inverted-topology, homodimeric fluoride channel Fluc in lipid bilayers. While the wild-type channel is a long-lived dimer, we leverage a known mutation, N43S, that weakens Na+ binding in a buried site at the interface, thereby unlocking the complex for reversible association in lipid bilayers. Single-channel recordings show that Na+ binding is required for fluoride conduction while single-molecule microscopy experiments demonstrate that N43S Fluc exists in a dynamic monomer-dimer equilibrium in the membrane, even following removal of Na+. Quantifying the thermodynamic stability while titrating Na+ indicates that dimerization occurs first, providing a membrane-embedded binding site where Na+ binding weakly stabilizes the complex. To understand how these subunits form stable assemblies while presenting charged surfaces to the membrane, we carried out molecular dynamics simulations, which show the formation of a thinned membrane defect around the exposed dimerization interface. In simulations where subunits are permitted to encounter each other while preventing protein contacts, we observe spontaneous and selective association at the native interface, where stability is achieved by mitigation of the membrane defect. These results suggest a model wherein membrane-associated forces drive channel assembly in the native orientation while subsequent factors, such as Na+ binding, result in channel activation.

Mechanism of 4-aminopyridine inhibition of the lysosomal channel TMEM175.

Transmembrane protein 175 (TMEM175) is an evolutionarily distinct lysosomal cation channel whose mutation is associated with the development of Parkinson's disease. Here, we present a cryoelectron microscopy structure and molecular simulations of TMEM175 bound to 4-aminopyridine (4-AP), the only known small-molecule inhibitor of TMEM175 and a broad K+ channel inhibitor, as well as a drug approved by the Food and Drug Administration against multiple sclerosis. The structure shows that 4-AP, whose mode of action had not been previously visualized, binds near the center of the ion conduction pathway, in the open state of the channel. Molecular dynamics simulations reveal that this binding site is near the middle of the transmembrane potential gradient, providing a rationale for the voltage-dependent dissociation of 4-AP from TMEM175. Interestingly, bound 4-AP rapidly switches between three predominant binding poses, stabilized by alternate interaction patterns dictated by the twofold symmetry of the channel. Despite this highly dynamic binding mode, bound 4-AP prevents not only ion permeation but also water flow. Together, these studies provide a framework for the rational design of novel small-molecule inhibitors of TMEM175 that might reveal the role of this channel in human lysosomal physiology both in health and disease.

Bivalent recognition of fatty acyl-CoA by a human integral membrane palmitoyltransferase.

S-acylation, also known as palmitoylation, is the most abundant form of protein lipidation in humans. This reversible posttranslational modification, which targets thousands of proteins, is catalyzed by 23 members of the DHHC family of integral membrane enzymes. DHHC enzymes use fatty acyl-CoA as the ubiquitous fatty acyl donor and become autoacylated at a catalytic cysteine; this intermediate subsequently transfers the fatty acyl group to a cysteine in the target protein. Protein S-acylation intersects with almost all areas of human physiology, and several DHHC enzymes are considered as possible therapeutic targets against diseases such as cancer. These efforts would greatly benefit from a detailed understanding of the molecular basis for this crucial enzymatic reaction. Here, we combine X-ray crystallography with all-atom molecular dynamics simulations to elucidate the structure of the precatalytic complex of human DHHC20 in complex with palmitoyl CoA. The resulting structure reveals that the fatty acyl chain inserts into a hydrophobic pocket within the transmembrane spanning region of the protein, whereas the CoA headgroup is recognized by the cytosolic domain through polar and ionic interactions. Biochemical experiments corroborate the predictions from our structural model. We show, using both computational and experimental analyses, that palmitoyl CoA acts as a bivalent ligand where the interaction of the DHHC enzyme with both the fatty acyl chain and the CoA headgroup is important for catalytic chemistry to proceed. This bivalency explains how, in the presence of high concentrations of free CoA under physiological conditions, DHHC enzymes can efficiently use palmitoyl CoA as a substrate for autoacylation.

Asymmetric drug binding in an ATP-loaded inward-facing state of an ABC transporter.

Substrate efflux by ATP-binding cassette (ABC) transporters, which play a major role in multidrug resistance, entails the ATP-powered interconversion between transporter intermediates. Despite recent progress in structure elucidation, a number of intermediates have yet to be visualized and mechanistically interpreted. Here, we combine cryogenic-electron microscopy (cryo-EM), double electron-electron resonance spectroscopy and molecular dynamics simulations to profile a previously unobserved intermediate of BmrCD, a heterodimeric multidrug ABC exporter from Bacillus subtilis. In our cryo-EM structure, ATP-bound BmrCD adopts an inward-facing architecture featuring two molecules of the substrate Hoechst-33342 in a striking asymmetric head-to-tail arrangement. Deletion of the extracellular domain capping the substrate-binding chamber or mutation of Hoechst-coordinating residues abrogates cooperative stimulation of ATP hydrolysis. Together, our findings support a mechanistic role for symmetry mismatch between the nucleotide binding and the transmembrane domains in the conformational cycle of ABC transporters and is of notable importance for rational design of molecules for targeted ABC transporter inhibition.

MOSAICS: A software suite for analysis of membrane structure and dynamics in simulated trajectories.

Molecular dynamics (MD) simulations have become the predominant computational analysis method in membrane biophysics, as this technique is uniquely suited for investigations of complex molecular systems through the relevant physical principles. Owing to continued improvements in scope and performance, the trajectories generated through this approach contain ever-increasing amounts of information, which must be synthesized and simplified in post-analysis using tools that are not only mechanistically insightful but also computationally efficient and highly scalable. Here, we introduce MOSAICS, a self-contained high-performance suite of C++ software tools designed for advanced analyses of lipid bilayer structure and dynamics from MD trajectories. MOSAICS is to our knowledge the most comprehensive software suite of this kind, enabling analysis of a wide array of morphological and kinetic properties, for both simple and complex membranes, irrespective of system size or resolution. Importantly, MOSAICS is designed to provide spatial distributions of all computed quantities, with built-in masking tools, noise filtering, and statistical significance metrics to facilitate quantitative interpretations of the trajectory data; it is also fully parallelized and can therefore leverage the capabilities of supercomputing facilities. Despite its technical sophistication, MOSAICS is user-friendly and requires minimal computational expertise, making it accessible to researchers of all skill levels. This sofware suite can be freely downloaded at https://github.com/MOSAICS-NIH/.

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

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

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

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

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

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

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

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

TRPV1 pore turret dictates distinct DkTx and capsaicin gating.

Many neurotoxins inflict pain by targeting receptors expressed on nociceptors, such as the polymodal cationic channel TRPV1. The tarantula double-knot toxin (DkTx) is a peptide with an atypical bivalent structure, providing it with the unique capability to lock TRPV1 in its open state and evoke an irreversible channel activation. Here, we describe a distinct gating mechanism of DkTx-evoked TRPV1 activation. Interestingly, DkTx evokes significantly smaller TRPV1 macroscopic currents than capsaicin, with a significantly lower unitary conductance. Accordingly, while capsaicin evokes aversive behaviors in TRPV1-transgenic Caenorhabditis elegans, DkTx fails to evoke such response at physiological concentrations. To determine the structural feature(s) responsible for this phenomenon, we engineered and evaluated a series of mutated toxins and TRPV1 channels. We found that elongating the DkTx linker, which connects its two knots, increases channel conductance compared with currents elicited by the native toxin. Importantly, deletion of the TRPV1 pore turret, a stretch of amino acids protruding out of the channel's outer pore region, is sufficient to produce both full conductance and aversive behaviors in response to DkTx. Interestingly, this deletion decreases the capsaicin-evoked channel activation. Taken together with structure modeling analysis, our results demonstrate that the TRPV1 pore turret restricts DkTx-mediated pore opening, probably through steric hindrance, limiting the current size and mitigating the evoked downstream physiological response. Overall, our findings reveal that DkTx and capsaicin elicit distinct TRPV1 gating mechanisms and subsequent pain responses. Our results also indicate that the TRPV1 pore turret regulates the mechanisms of channel gating and permeation.

Interpretation of HDX Data by Maximum-Entropy Reweighting of Simulated Structural Ensembles.

Hydrogen-deuterium exchange combined with mass spectrometry (HDX-MS) is a widely applied biophysical technique that probes the structure and dynamics of biomolecules without the need for site-directed modifications or bio-orthogonal labels. The mechanistic interpretation of HDX data, however, is often qualitative and subjective, owing to a lack of quantitative methods to rigorously translate observed deuteration levels into atomistic structural information. To help address this problem, we have developed a methodology to generate structural ensembles that faithfully reproduce HDX-MS measurements. In this approach, an ensemble of protein conformations is first generated, typically using molecular dynamics simulations. A maximum-entropy bias is then applied post hoc to the resulting ensemble such that averaged peptide-deuteration levels, as predicted by an empirical model, agree with target values within a given level of uncertainty. We evaluate this approach, referred to as HDX ensemble reweighting (HDXer), for artificial target data reflecting the two major conformational states of a binding protein. We demonstrate that the information provided by HDX-MS experiments and by the model of exchange are sufficient to recover correctly weighted structural ensembles from simulations, even when the relevant conformations are rarely observed. Degrading the information content of the target data-e.g., by reducing sequence coverage, by averaging exchange levels over longer peptide segments, or by incorporating different sources of uncertainty-reduces the structural accuracy of the reweighted ensemble but still allows for useful insights into the distinctive structural features reflected by the target data. Finally, we describe a quantitative metric to rank candidate structural ensembles according to their correspondence with target data and illustrate the use of HDXer to describe changes in the conformational ensemble of the membrane protein LeuT. In summary, HDXer is designed to facilitate objective structural interpretations of HDX-MS data and to inform experimental approaches and further developments of theoretical exchange models.

DHHC20 Palmitoyl-Transferase Reshapes the Membrane to Foster Catalysis.

Cysteine palmitoylation, a form of S-acylation, is a key posttranslational modification in cellular signaling. This type of reversible lipidation occurs in both plasma and organellar membranes, and is catalyzed by a family of integral membrane proteins known as DHHC acyltransferases. The first step in the S-acylation process is the recognition of free acyl coenzyme A (acyl-CoA) from the lipid bilayer. The DHHC enzyme then becomes autoacylated at a site defined by a conserved Asp-His-His-Cys motif. This reaction entails ionization of the catalytic Cys. Intriguingly, in known DHHC structures, this catalytic Cys appears to be exposed to the hydrophobic interior of the lipid membrane, which would be highly unfavorable for a negatively charged nucleophile, thus hindering autoacylation. Here, we use biochemical and computational methods to reconcile these seemingly contradictory facts. First, we experimentally demonstrate that human DHHC20 is active when reconstituted in POPC nanodiscs. Microsecond-long all-atom molecular dynamics simulations are then calculated for human DHHC20 and for different acyl-CoA forms, also in a POPC membrane. Strikingly, we observe that human DHHC20 induces a drastic deformation in the membrane, particularly on the cytoplasmic side, where autoacylation occurs. As a result, the catalytic Cys becomes hydrated and optimally positioned to encounter the cleavage site in acyl-CoA. In summary, we hypothesize that DHHC enzymes locally reshape the membrane to foster a morphology that is specifically adapted for acyl-CoA recognition and autoacylation.

Direct Derivation of Free Energies of Membrane Deformation and Other Solvent Density Variations From Enhanced Sampling Molecular Dynamics.

We report a methodology to calculate the free energy of a shape transformation in a lipid membrane directly from a molecular dynamics simulation. The bilayer need not be homogeneous or symmetric and can be atomically detailed or coarse grained. The method is based on a collective variable that quantifies the similarity between the membrane and a set of predefined density distributions. Enhanced sampling of this "Multi-Map" variable re-shapes the bilayer and permits the derivation of the corresponding potential of mean force. Calculated energies thus reflect the dynamic interplay of atoms and molecules, rather than postulated effects. Evaluation of deformations of different shape, amplitude, and range demonstrates that the macroscopic bending modulus assumed by the Helfrich-Canham model is increasingly unsuitable below the 100-Å scale. In this range of major biological significance, direct free-energy calculations reveal a much greater plasticity. We also quantify the stiffening effect of cholesterol on bilayers of different composition and compare with experiments. Lastly, we illustrate how this approach facilitates analysis of other solvent reorganization processes, such as hydrophobic hydration. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.

The 2.5 Å structure of CD1c in complex with a mycobacterial lipid reveals an open groove ideally suited for diverse antigen presentation.

CD1 molecules function to present lipid-based antigens to T cells. Here we present the crystal structure of CD1c at 2.5 Å resolution, in complex with the pathogenic Mycobacterium tuberculosis antigen mannosyl-ß1-phosphomycoketide (MPM). CD1c accommodated MPM's methylated alkyl chain exclusively in the A' pocket, aided by a unique exit portal underneath the α1 helix. Most striking was an open F' pocket architecture lacking the closed cavity structure of other CD1 molecules, reminiscent of peptide binding grooves of classical major histocompatibility complex molecules. This feature, combined with tryptophan-fluorescence quenching during loading of a dodecameric lipopeptide antigen, provides a compelling model by which both the lipid and peptide moieties of the lipopeptide are involved in CD1c presentation of lipopeptides.

Confidence Analysis of DEER Data and Its Structural Interpretation with Ensemble-Biased Metadynamics.

Given its ability to measure multicomponent distance distributions between electron-spin probes, double electron-electron resonance (DEER) spectroscopy has become a leading technique to assess the structural dynamics of biomolecules. However, methodologies to evaluate the statistical error of these distributions are not standard, often hampering a rigorous interpretation of the experimental results. Distance distributions are often determined from the experimental DEER data through a mathematical method known as Tikhonov regularization, but this approach makes rigorous error estimates difficult. Here, we build upon an alternative, model-based approach in which the distance probability distribution is represented as a sum of Gaussian components, and use propagation of errors to calculate an associated confidence band. Our approach considers all sources of uncertainty, including the experimental noise, the uncertainty in the fitted background signal, and the limited time span of the data collection. The resulting confidence band reveals the most and least reliable features of the probability distribution, thereby informing the structural interpretation of DEER experiments. To facilitate this interpretation, we also generalize the molecular simulation method known as ensemble-biased metadynamics (EBMetaD). This method, originally designed to generate maximal-entropy structural ensembles consistent with one or more probability distributions, now also accounts for the uncertainty in those target distributions exactly as dictated by their confidence bands. After careful benchmarks, we demonstrate the proposed techniques using DEER results from spin-labeled T4 lysozyme.

Membrane plasticity facilitates recognition of the inhibitor oligomycin by the mitochondrial ATP synthase rotor.

Enzymes in the respiratory chain are increasingly seen as potential targets against multi-drug resistance of human pathogens and cancerous cells. However, a detailed understanding of the mechanism and specificity determinants of known inhibitors is still lacking. Oligomycin, for example, has been known to be an inhibitor of the membrane motor of the mitochondrial ATP synthase for over five decades, and yet little is known about its mode of action at the molecular level. In a recent breakthrough, a crystal structure of the S. cerevisiae c-subunit ring with bound oligomycin revealed the inhibitor docked on the outer face of the proton-binding sites, deep into the transmembrane region. However, the structure of the complex was obtained in an organic solvent rather than detergent or a lipid bilayer, and therefore it has been unclear whether this mode of recognition is physiologically relevant. Here, we use molecular dynamics simulations to address this question and gain insights into the mechanism of oligomycin inhibition. Our findings lead us to propose that oligomycin naturally partitions into the lipid/water interface, and that in this environment the inhibitor can indeed bind to any of the c-ring proton-carrying sites that are exposed to the membrane, thereby becoming an integral component of the proton-coordinating network. As the c-ring rotates within the membrane, driven either by downhill proton permeation or ATP hydrolysis, one of the protonated, oligomycin-bound sites eventually reaches the subunit-a interface and halts the rotary mechanism of the enzyme.

On the principle of ion selectivity in Na+/H+-coupled membrane proteins: experimental and theoretical studies of an ATP synthase rotor.

Numerous membrane transporters and enzymes couple their mechanisms to the permeation of Na(+) or H(+), thereby harnessing the energy stored in the form of transmembrane electrochemical potential gradients to sustain their activities. The molecular and environmental factors that control and modulate the ion specificity of most of these systems are, however, poorly understood. Here, we use isothermal titration calorimetry to determine the Na(+)/H(+) selectivity of the ion-driven membrane rotor of an F-type ATP synthase. Consistent with earlier theoretical predictions, we find that this rotor is significantly H(+) selective, although not sufficiently to be functionally coupled to H(+), owing to the large excess of Na(+) in physiological settings. The functional Na(+) specificity of this ATP synthase thus results from two opposing factors, namely its inherent chemical selectivity and the relative availability of the coupling ion. Further theoretical studies of this membrane rotor, and of two others with a much stronger and a slightly weaker H(+) selectivity, indicate that, although the inherent selectivity of their ion-binding sites is largely set by the balance of polar and hydrophobic groups flanking a conserved carboxylic side chain, subtle variations in their structure and conformational dynamics, for a similar chemical makeup, can also have a significant contribution. We propose that the principle of ion selectivity outlined here may provide a rationale for the differentiation of Na(+)- and H(+)-coupled systems in other families of membrane transporters and enzymes.

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

Na(+)/Ca(2+) exchangers (NCXs) are ubiquitous membrane transporters with a key role in Ca(2+) homeostasis and signaling. NCXs mediate the bidirectional translocation of either Na(+) or Ca(2+), and thus can catalyze uphill Ca(2+) transport driven by a Na(+) gradient, or vice versa. In a major breakthrough, a prokaryotic NCX homolog (NCX_Mj) was recently isolated and its crystal structure determined at atomic resolution. The structure revealed an intriguing architecture consisting of two inverted-topology repeats, each comprising five transmembrane helices. These repeats adopt asymmetric conformations, yielding an outward-facing occluded state. The crystal structure also revealed four putative ion-binding sites, but the occupancy and specificity thereof could not be conclusively established. Here, we use molecular-dynamics simulations and free-energy calculations to identify the ion configuration that best corresponds to the crystallographic data and that is also thermodynamically optimal. In this most probable configuration, three Na(+) ions occupy the so-called Sext, SCa, and Sint sites, whereas the Smid site is occupied by one water molecule and one H(+), which protonates an adjacent aspartate side chain (D240). Experimental measurements of Na(+)/Ca(2+) and Ca(2+)/Ca(2+) exchange by wild-type and mutagenized NCX_Mj confirm that transport of both Na(+) and Ca(2+) requires protonation of D240, and that this side chain does not coordinate either ion at Smid. These results imply that the ion exchange stoichiometry of NCX_Mj is 3:1 and that translocation of Na(+) across the membrane is electrogenic, whereas transport of Ca(2+) is not. Altogether, these findings provide the basis for further experimental and computational studies of the conformational mechanism of this exchanger.