Search

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. [ABSTRACT FROM AUTHOR]

Transmembrane protein 175 (TMEM175) is an evolutionary 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The atomic mechanisms of isomerization of ATP-Mg2+ in solution are characterized using the recently developed String Method with Optimal Molecular Alignment (SOMA) and molecular-dynamics simulations. Bias-Exchange Metadynamics simulations are first performed to identify the primary conformers of the ATP-Mg2+ complex and their connectivity. SOMA is then used to elucidate the minimum free-energy path (MFEP) for each transition, in a 48-dimensional space. Analysis of the per-atom contributions to the global free-energy profiles reveals that the mechanism of these transitions is controlled by the Mg2+ ion and its coordinating oxygen atoms in the triphosphate moiety, as well as by the ion-hydration shell. Metadynamics simulations in path collective variables based on the MFEP demonstrate these isomerizations proceed across a narrow channel of configurational space, thus validating the premise underlying SOMA. This study provides a roadmap for the examination of conformational changes in biomolecules, based on complementary enhanced-sampling techniques with different strengths. © 2015 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

ATG9, the only transmembrane protein in the core macroautophagy/autophagy machinery, is a key player in the early stages of autophagosome formation. Yet, the lack of a high-resolution structure of ATG9 was a major impediment in understanding its three-dimensional organization and function. We recently solved a high-resolution cryoEM structure of the ubiquitously expressed human ATG9A isoform. The structure revealed that ATG9A is a domain-swapped homotrimer with a unique fold, and has an internal network of branched cavities. In cellulo analyses demonstrated the functional importance of the cavity-lining residues. These cavities could serve as conduits for transport of hydrophilic moieties, such as lipid headgroups, across the bilayer. Finally, structure-guided molecular dynamics predicted that ATG9A has membrane-bending properties, which is consistent with its localization to highly curved membranes. [ABSTRACT FROM AUTHOR]

The proton pore of the F1Fo ATP synthase consists of a ring of c subunits, which rotates, driven by downhill proton diffusion across the membrane. An essential carboxylate side chain in each subunit provides a proton-binding site. In all the structures of c-rings reported to date, these sites are in a closed, ion-locked state. Structures are here presented of the c10 ring from Saccharomyces cerevisiae determined at pH 8.3, 6.1 and 5.5, at resolutions of 2.0 Å, 2.5 Å and 2.0 Å, respectively. The overall structure of this mitochondrial c-ring is similar to known homologs, except that the essential carboxylate, Glu59, adopts an open extended conformation. Molecular dynamics simulations reveal that opening of the essential carboxylate is a consequence of the amphiphilic nature of the crystallization buffer. We propose that this new structure represents the functionally open form of the c subunit, which facilitates proton loading and release. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations based on the replica-exchange framework (REMD) are emerging as a useful tool to characterize the conformational variability that is intrinsic to most chemical and biological systems. In this work, it is shown that a simple extension of the replica-exchange method, known as Hamiltonian REMD, greatly facilitates the characterization of conformational equilibria across large energetic barriers, or in the presence of substantial entropic effects, overcoming some of the difficulties of REMD based on temperature alone. In particular, a comparative assessment of the HREMD and TREMD approaches was made, through computation of the gas-phase free-energy difference between the so-called D2d and S4 states of tetrabutylammonium (TBA), an ionic compound of frequently used in biophysical studies of ion channels. Taking advantage of the greater efficiency of the HREMD scheme, the conformational equilibrium of TBA was characterized in a variety of conditions. Simulation of the gas-phase equilibrium in the 100–300 K range allowed us to compute the entropy difference between these states as well as to describe its temperature dependence. Through HREMD simulations of TBA in a water droplet, the effect of solvation on the conformational equilibrium was determined. Finally, the equilibrium of TBA in the context of a simplified model of the binding cavity of the KcsA potassium channel was simulated, and density maps for D2d and S4 states analogous to those derived from X-ray crystallography were constructed. Overall, this work illustrates the potential of the HREMD approach in the context of computational drug design, ligand-receptor structural prediction and more generally, molecular recognition, where one of the most challenging issues remains to account for conformational flexibility as well for the solvation and entropic effects thereon. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2007 [ABSTRACT FROM AUTHOR]

Abstract: The mechanism of intracellular blockade of the KcsA potassium channel by tetrabutylammonium (TBA) is investigated through functional, structural and computational studies. Using planar-membrane electrophysiological recordings, we characterize the binding kinetics as well as the dependence on the transmembrane voltage and the concentration of the blocker. It is found that the apparent affinity of the complex is significantly greater than that of any of the eukaryotic K+ channels studied previously, and that the off-rate increases with the applied transmembrane voltage. In addition, we report a crystal structure of the KcsA–TBA complex at 2.9 Å resolution, with TBA bound inside the large hydrophobic cavity located at the center of the channel, consistent with the results of previous functional and structural studies. Of particular interest is the observation that the presence of TBA has a negligible effect on the channel structure and on the position of the potassium ions occupying the selectivity filter. Inspection of the electron density corresponding to TBA suggests that the ligand may adopt more than one conformation in the complex, though the moderate resolution of the data precludes a definitive interpretation on the basis of the crystallographic refinement methods alone. To provide a rationale for these observations, we carry out an extensive conformational sampling of an atomic model of TBA bound in the central cavity of KcsA, using the Hamiltonian replica-exchange molecular dynamics simulation method. Comparison of the simulated and experimental density maps indicates that the latter does reflect at least two distinct binding orientations of TBA. The simulations show also that the relative population of these binding modes is dependent on the ion configuration occupying the selectivity filter, thus providing a clue to the nature of the voltage-dependence of the binding kinetics. [Copyright &y& Elsevier]

In this paper we describe a method for setting up an atomistic simulation of a membrane protein in a hydrated lipid bilayer and report the effect of differing electrostatic parameters on the drift in the protein structure during the subsequent simulation. The method aims to generate a suitable cavity in the interior of a lipid bilayer, using the solvent-accessible surface of the protein as a template, during the course of a short steered molecular dynamics simulation of a solvated lipid membrane. This is achieved by a two-stage process: firstly, lipid molecules whose headgroups are inside a cylindrical volume equivalent to that defined by the protein surface are removed; then the protein-lipid interface is optimized by applying repulsive forces perpendicular to the protein surface, and of gradually increased magnitude, to the remaining lipid atoms inside the volume occupied by the protein surface until it is emptied. The protein itself may then be inserted. Using the bacterial membrane proteins KcsA and FhuA as test cases, we show how the method achieves the formation of a suitable cavity in the interior of a dimyristoylphosphatidylcholine lipid bilayer without perturbing the configuration of the non-interfacial regions of the previously equilibrated lipid bilayer, even in cases of membrane proteins with irregular geometrical shapes. In addition, we compare subsequent simulations in which the long-range electrostatic interactions are treated via either a cut-off or a particle-mesh Ewald method. The results show that the drift from the initial structure is less in the latter case, especially for KcsA and for the non-core secondary structural elements (i.e. surface loops) of both proteins. [ABSTRACT FROM AUTHOR]

Bedaquiline (BDQ, Sirturo) has been approved to treat multidrug resistant forms of Mycobacterium tuberculosis. Prior studies suggested that BDQ was a selective inhibitor of the ATP synthase from M. tuberculosis. However, Sirturo treatment leads to an increased risk of cardiac arrhythmias and death, raising the concern that this adverse effect results from inhibition at a secondary site. Here we show that BDQ is a potent inhibitor of the yeast and human mitochondrial ATP synthases. Single-particle cryo-EM reveals that the site of BDQ inhibition partially overlaps with that of the inhibitor oligomycin. Molecular dynamics simulations indicate that the binding mode of BDQ to this site is similar to that previously seen for a mycobacterial enzyme, explaining the observed lack of selectivity. We propose that derivatives of BDQ ought to be made to increase its specificity toward the mycobacterial enzyme and thereby reduce the side effects for patients that are treated with Sirturo. Luo, Zhou et al. show that Bedaquiline (BDQ, Sirturo), approved to treat multi-drug-resistant tuberculosis, inhibits the yeast and human mitochondrial ATP synthases in addition to its intended target, the Mycobacterium tuberculosis ATP synthase. The structure of the mitochondrial ATP synthase bound to BDQ suggests a means to modify this inhibitor to increase its specificity for the M. tuberculosis enzyme, thereby reducing its side effects for patients. [ABSTRACT FROM AUTHOR]

Autophagy is a catabolic process involving capture of cytoplasmic materials into double-membraned autophagosomes that subsequently fuse with lysosomes for degradation of the materials by lysosomal hydrolases. One of the least understood components of the autophagy machinery is the transmembrane protein ATG9. Here, we report a cryoelectron microscopy structure of the human ATG9A isoform at 2.9-Å resolution. The structure reveals a fold with a homotrimeric domain-swapped architecture, multiple membrane spans, and a network of branched cavities, consistent with ATG9A being a membrane transporter. Mutational analyses support a role for the cavities in the function of ATG9A. In addition, structure-guided molecular simulations predict that ATG9A causes membrane bending, explaining the localization of this protein to small vesicles and highly curved edges of growing autophagosomes. • The transmembrane autophagy protein ATG9A is a domain-swapped homotrimer • Each ATG9A protomer comprises four transmembrane domains • The ATG9A homotrimer exhibits an internal network of branched cavities • Molecular dynamics simulations show that ATG9A trimers deform membranes Guardia et al. report a high-resolution cryo-EM structure of human ATG9A, the only transmembrane protein of the core autophagy machinery. The structure shows that ATG9A is a domain-swapped homotrimer with a complex network of internal cavities. Structure-based computational simulations predict that ATG9A has membrane-bending properties. [ABSTRACT FROM AUTHOR]

S-acylation, whereby a fatty acid chain is covalently linked to a cysteine residue by a thioester linkage, is the most prevalent kind of lipid modification of proteins. Thousands of proteins are targets of this post-translational modification, which is catalyzed by a family of eukaryotic integral membrane enzymes known as DHHC protein acyltransferases (DHHC-PATs). Our knowledge of the repertoire of S-acylated proteins has been rapidly expanding owing to development of the chemoproteomic techniques. There has also been an increasing number of reports in the literature documenting the importance of S-acylation in human physiology and disease. Recently, the first atomic structures of two different DHHC-PATs were determined using X-ray crystallography. This review will focus on the insights gained into the molecular mechanism of DHHC-PATs from these structures and highlight representative data from the biochemical literature that they help explain. Unlabelled Image • Protein S-acylation is one of the most abundant forms of protein lipidation. • Members of the DHHC family of integral membrane enzymes catalyze protein S-acylation in eukaryotes. • High-resolution structures of two members of the family revealed a tremendous amount of information about organization and mechanism of DHHC enzymes. • Atomistic molecular dynamics simulation revealed DHHC enzymes deform the membrane to facilitate catalysis. [ABSTRACT FROM AUTHOR]

Adenosine triphosphate (ATP) synthases populate the inner membranes of mitochondria, where they produce the majority of the ATP required by the cell. From yeast to vertebrates, cryoelectron tomograms of these membranes have consistently revealed a very precise organization of these enzymes. Rather than being scattered throughout the membrane, the ATP synthases form dimers, and these dimers are organized into rows that extend for hundreds of nanometers. The rows are only observed in the membrane invaginations known as cristae, specifically along their sharply curved edges. Although the presence of these macromolecular structures has been irrefutably linked to the proper development of cristae morphology, it has been unclear what drives the formation of the rows and why they are specifically localized in the cristae. In this study, we present a quantitative molecular-simulation analysis that strongly suggests that the dimers of ATP synthases organize into rows spontaneously, driven by a long-range attractive force that arises from the relief of the overall elastic strain of the membrane. The strain is caused by the V-like shape of the dimers, unique among membrane protein complexes, which induces a strong deformation in the surrounding membrane. The process of row formation is therefore not a result of direct protein-protein interactions or a specific lipid composition of the membrane. We further hypothesize that, once assembled, the ATP synthase dimer rows prime the inner mitochondrial membrane to develop folds and invaginations by causing macroscopic membrane ridges that ultimately become the edges of cristae. In this way, mitochondrial ATP synthases would contribute to the generation of a morphology that maximizes the surface area of the inner membrane, and thus ATP production. Finally, we outline key experiments that would be required to verify or refute this hypothesis. [ABSTRACT FROM AUTHOR]

The force-from-lipids hypothesis of cellular mechanosensation posits that membrane channels open and close in response to changes in the physical state of the lipid bilayer, induced for example by lateral tension. Here, we investigate the molecular basis for this transduction mechanism by studying the mechanosensitive ion channel MscS from Escherichia coli and its eukaryotic homolog MSL1 from Arabidopsis thaliana. First, we use single-particle cryo-electron microscopy to determine the structure of a novel open conformation of wild-type MscS, stabilized in a thinned lipid nanodisc. Compared with the closed state, the structure shows a reconfiguration of helices TM1, TM2, and TM3a, and widening of the central pore. Based on these structures, we examined how the morphology of the membrane is altered upon gating, using molecular dynamics simulations. The simulations reveal that closed-state MscS causes drastic protrusions in the inner leaflet of the lipid bilayer, both in the absence and presence of lateral tension, and for different lipid compositions. These deformations arise to provide adequate solvation to hydrophobic crevices under the TM1-TM2 hairpin, and clearly reflect a high-energy conformation for the membrane, particularly under tension. Strikingly, these protrusions are largely eradicated upon channel opening. An analogous computational study of open and closed MSL1 recapitulates these findings. The gating equilibrium of MscS channels thus appears to be dictated by opposing conformational preferences, namely those of the lipid membrane and of the protein structure. We propose a membrane deformation model of mechanosensation, which posits that tension shifts the gating equilibrium towards the conductive state not because it alters the mode in which channel and lipids interact, but because it increases the energetic cost of the morphological perturbations in the membrane required by the closed state. [ABSTRACT FROM AUTHOR]

Structures of the human lysosomal K+ channel transmembrane protein 175 (TMEM175) in open and closed states revealed a novel architecture lacking the canonical K+ selectivity filter motif present in previously known K+ channel structures. A hydrophobic constriction composed of four isoleucine residues was resolved in the pore and proposed to serve as the gate in the closed state, and to confer ion selectivity in the open state. Here, we achieve higher-resolution structures of the open and closed states and employ molecular dynamics simulations to analyze the conducting properties of the putative open state, demonstrating that it is permeable to K+ and, to a lesser degree, also Na+. Both cations must dehydrate significantly to penetrate the narrow hydrophobic constriction, but ion flow is assisted by a favorable electrostatic field generated by the protein that spans the length of the pore. The balance of these opposing energetic factors explains why permeation is feasible, and why TMEM175 is selective for K+ over Na+, despite the absence of the canonical selectivity filter. Accordingly, mutagenesis experiments reveal an exquisite sensitivity of the channel to perturbations that mitigate the constriction. Together, these data reveal a novel mechanism for selective permeation of ions by TMEM175 that is unlike that of other K+ channels. [ABSTRACT FROM AUTHOR]

In both prokaryotes and eukaryotes, multidrug and toxiccompound extrusion (MATE) transporters catalyze the efflux of a broad range of cytotoxic compounds, including humanmade 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 Nterminal 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. [ABSTRACT FROM AUTHOR]

Pathological metabolic conditions such as ischemia induce the rupture of the mitochondrial envelope and the release of pro-apoptotic proteins, leading to cell death. At the onset of this process, the inner mitochondrial membrane becomes depolarized and permeable to osmolytes, proposedly due to the opening of a non-selective protein channel of unknown molecular identity. A recent study purports that this channel, referred to as Mitochondrial Permeability Transition Pore (MPTP), is formed within the c-subunit ring of the ATP synthase, upon its dissociation from the catalytic domain of the enzyme. Here, we examine this claim for two c-rings of different lumen width, through calculations of their ion conductance and selectivity based on allatom molecular dynamics simulations. We also quantify the likelihood that the lumen of these c-rings is in a hydrated, potentially conducting state rather than empty or blocked by lipid molecules. These calculations demonstrate that the structure and biophysical properties of a correctly assembled c-ring are inconsistent with those attributed to the MPTP. [ABSTRACT FROM AUTHOR]

Abstract: The functional mechanism of the F1Fo ATP synthase, like many membrane transporters and pumps, entails a conformational cycle that is coupled to the movement of H+ or Na+ ions across its transmembrane domain, down an electrochemical gradient. This coupling is an efficient means of energy transduction and regulation, provided that ion binding to the membrane domain, known as Fo, is appropriately selective. In this study we set out to establish the structural and energetic basis for the ion-binding selectivity of the membrane-embedded Fo rotors of two representative ATP synthases. First, we use a biochemical approach to demonstrate the inherent binding selectivity of these rotors, that is, independently from the rest of the enzyme. We then use atomically detailed computer simulations of wild-type and mutagenized rotors to calculate and rationalize their selectivity, on the basis of the structure, dynamics and coordination chemistry of the binding sites. We conclude that H+ selectivity is most likely a robust property of all Fo rotors, arising from the prominent presence of a conserved carboxylic acid and its intrinsic chemical propensity for protonation, as well as from the structural plasticity of the binding sites. In H+-coupled rotors, the incorporation of hydrophobic side chains to the binding sites enhances this inherent H+ selectivity. Size restriction may also favor H+ over Na+, but increasing size alone does not confer Na+ selectivity. Rather, the degree to which Fo rotors may exhibit Na+ coupling relies on the presence of a sufficient number of suitable coordinating side chains and/or structural water molecules. These ligands accomplish a shift in the relative binding energetics, which under some physiological conditions may be sufficient to provide Na+ dependence. [Copyright &y& Elsevier]

Abstract: The recently reported crystal structures of the membrane-embedded proton-dependent c-ring rotors of a cyanobacterial F1Fo ATP synthase and a chloroplast F1Fo ATP synthase have provided new insights into the mechanism of this essential enzyme. While the overall features of these c-rings are similar, a discrepancy in the structure and hydrogen-bonding interaction network of the H+ sites suggests two distinct binding modes, potentially reflecting a mechanistic differentiation. Importantly, the conformation of the key glutamate side chain to which the proton binds is also altered. To investigate the nature of these differences, we use molecular dynamics simulations of both c-rings embedded in a phospholipid membrane. We observe that the structure of the c15 ring from Spirulina platensis is unequivocally stable within the simulation time. By contrast, the proposed structure of the H+ site in the chloroplast c14 ring changes rapidly and consistently into that reported for the c15 ring, indicating that the latter represents a common binding mode. To assess this hypothesis, we have remodeled the c14 ring by molecular replacement using the published structure factors. The resulting structure provides clear evidence in support of a common binding site conformation and is also considerably improved statistically. These findings, taken together with a sequence analysis of c-subunits in the ATP synthase family, indicate that the so-called proton-locked conformation observed in the c15 ring may be a common characteristic not only of light-driven systems such as chloroplasts and cyanobacteria but also of a selection of other bacterial species. [Copyright &y& Elsevier]

Summary: The membrane-embedded rotors of Na+-dependent F-ATP synthases comprise 11 c-subunits that form a ring, with 11 Na+ binding sites in between adjacent subunits. Following an updated crystallographic analysis of the c-ring from Ilyobacter tartaricus, we report the complete ion-coordination structure of the Na+ sites. In addition to the four residues previously identified, there exists a fifth ligand, namely, a buried structural water molecule. This water is itself coordinated by Thr67, which, sequence analysis reveals, is the only residue involved in binding that distinguishes Na+ synthases from H+-ATP synthases known to date. Molecular dynamics simulations and free-energy calculations of the c-ring in a lipid membrane lend clear support to the notion that this fifth ligand is a water molecule, and illustrate its influence on the selectivity of the binding sites. Given the evolutionary ascendancy of sodium over proton bioenergetics, this structure uncovers an ancient strategy for selective ion coupling in ATP synthases. [Copyright &y& Elsevier]