The extracellular Ca2+ concentration changes locally under certain physiological and pathological conditions. Such variations affect the function of ion channels of the nervous system and consequently also neuronal signalling. We investigated here the mechanisms by which Ca2+ controls the activity of acid-sensing ion channel (ASIC) 3. ASICs are neuronal, H+-gated Na+ channels involved in several physiological and pathological processes, including the expression of fear, learning, pain sensation and neurodegeneration after ischaemic stroke. It was previously shown that Ca2+ negatively modulates the ASIC pH dependence. While protons are default activators of ASIC3, this channel can also be activated at pH7.4 by the removal of the extracellular Ca2+. Two previous studies concluded that low pH opens ASIC3 by displacing Ca2+ ions that block the channel pore at physiological pH. We show here that an acidic residue, distant from the pore, together with pore residues, controls the modulation of ASIC3 by Ca2+. Our study identifies a new regulatory site in ASIC3 and demonstrates that ASIC3 activation involves an allosteric mechanism together with Ca2+ unbinding from the channel pore. We provide a molecular analysis of a regulatory mechanism found in many ion channels.

The CLC family comprises H+-coupled exchangers and Cl- channels, and mutations causing their dysfunction lead to genetic disorders. The CLC exchangers, unlike canonical 'ping-pong' antiporters, simultaneously bind and translocate substrates through partially congruent pathways. How ions of opposite charge bypass each other while moving through a shared pathway remains unknown. Here, we use MD simulations, biochemical and electrophysiological measurements to identify two conserved phenylalanine residues that form an aromatic pathway whose dynamic rearrangements enable H+ movement outside the Cl- pore. These residues are important for H+ transport and voltage-dependent gating in the CLC exchangers. The aromatic pathway residues are evolutionarily conserved in CLC channels where their electrostatic properties and conformational flexibility determine gating. We propose that Cl- and H+ move through physically distinct and evolutionarily conserved routes through the CLC channels and transporters and suggest a unifying mechanism that describes the gating mechanism of both CLC subtypes.

Voltage-gated potassium channels (Kv) allow ion permeation upon changes of the membrane electrostatic potential (Vm). Each subunit of these tetrameric channels is composed of six transmembrane helices, of which the anti-parallel helix bundle S1-S4 constitutes the voltage-sensor domain (VSD) and S5-S6 forms the pore domain. Here, using 82 molecular dynamics (MD) simulations involving 266 replicated VSDs, we report novel responses of the archaebacterial potassium channel KvAP to membrane polarization. We show that the S4 α-helix, which is straight in the experimental crystal structure solved under depolarized conditions (Vm ∼ 0), breaks into two segments when the cell membrane is hyperpolarized (Vm << 0), and reversibly forms a single straight helix following depolarization (Vm = 0). The outermost segment of S4 translates along the normal to the membrane, bringing new perspective to previously paradoxical accessibility experiments that were initially thought to imply the displacement of the whole VSD across the membrane. The novel model is applied through steered and unbiased MD simulations to the recently solved whole structure of KvAP. The simulations show that the resting state involves a re-orientation of the S5 α-helix by ∼ 5–6 degrees in respect to the normal of the bilayer, which could result in the constriction and closure of the selectivity filter. Our findings support the idea that the breakage of S4 under (hyper)polarization is a general feature of Kv channels with a non-swapped topology.

Voltage-gated ion channels allow ion permeation upon changes of the membrane electrostatic potential (Vm). Each subunit of these tetrameric channels is composed of six transmembrane helices, of which the anti-parallel helix bundle S1-S4 constitutes the voltage-sensor domain (VSD) and S5-S6 forms the pore domain. Here, using molecular dynamics (MD) simulations, we report novel responses of the archaebacterial potassium channel KvAP to cell polarization. We show that the S4 helix, which is straight in the experimental crystal structure solved under depolarized conditions (Vm ∼ 0), breaks into two segments when the cell is polarized (Vm << 0), and reversibly forms a single straight helix following depolarization of the cell (Vm =0). The outermost segment of S4 translates along the normal to the membrane, bringing new perspective to previously paradoxical accessibility experiments that were initially thought to imply the displacement of the whole VSD across the membrane. Our simulations of KvAP reveal that the breaking of S4 under polarization is not a feature unique to hyperpolarization activated channel, as might be suggested by recent cryo-EM structures and MD simulations of the HCN channel.

In the arbuscular mycorrhizal (AM) symbiosis, plants satisfy part of their nitrogen (N) requirement through the AM pathway. In sorghum, the ammonium transporters (AMT) AMT3;1, and to a lesser extent AMT4, are induced in cells containing developing arbuscules. Here, we have characterized orthologs of AMT3;1 and AMT4 in four other grasses in addition to sorghum. AMT3;1 and AMT4 orthologous genes are induced in AM roots, suggesting that in the common ancestor of these five plant species, both AMT3;1 and AMT4 were already present and upregulated upon AM colonization. An artificial microRNA approach was successfully used to downregulate either AMT3;1 or AMT4 in rice. Mycorrhizal root colonization and hyphal length density of knockdown plants were not affected at that time, indicating that the manipulation did not modify the establishment of the AM symbiosis and the interaction between both partners. However, expression of the fungal phosphate transporter FmPT was significantly reduced in knockdown plants, indicating a reduction of the nutrient fluxes from the AM fungus to the plant. The AMT3;1 knockdown plants (but not the AMT4 knockdown plants) were significantly less stimulated in growth by AM fungal colonization, and uptake of both 15N and 33P from the AM fungal network was reduced. This confirms that N and phosphorus nutrition through the mycorrhizal pathway are closely linked. But most importantly, it indicates that AMT3;1 is the prime plant transporter involved in the mycorrhizal ammonium transfer and that its function during uptake of N cannot be performed by AMT4.

Potassium channels constitute a super-family of membrane proteins playing a key role in the function of neurons. Of particular importance are the voltage-gated potassium channels (Kv), of which there are around 40 sub-types. These channels have been studied extensively for several decades and great advances in the understanding of the mechanisms governing conduction and gating at the atomic scale have been made since the first crystal structure of a potassium channel was solved almost 20 years ago. Nevertheless, as molecular dynamics (MD) simulations have suggested that the conformational states captured in the crystal structures of voltage-gated potassium channels sustained permeation rates much lower compared to experimental measurements, the details of the ion conduction mechanism and whether the crystal structures truly represent the conducting state are still debated. We are using MD simulations to study the impact of membrane composition and voltage on the conformation of a Kv channel and its ion conductivity. Our results notably suggest that membrane thickness affects the function of the channel. These simulations illustrate how the microenvironment can impact the function of ion channels.

We used targeted molecular dynamics, informed by experimentally determined inter-atomic distances defining the pore region of open and closed states of the KvAP voltage-gated potassium channel, to generate a gating pathway of the pore domain in the absence of the voltage-sensing domains. We then performed umbrella sampling simulations along this pathway to calculate a potential of mean force that describes the free energy landscape connecting the closed and open conformations of the pore domain. The resulting energetic landscape displays three minima, corresponding to stable open, closed, and intermediate conformations with roughly similar stabilities. We found that the extent of hydration of the interior of the pore domain could influence the free energy landscape for pore opening/closing. © 2017 Wiley Periodicals, Inc.

ATP-binding cassette (ABC) transporters use the energy of ATP binding and hydrolysis to move a large variety of compounds across biological membranes. P-glycoprotein, involved in multidrug resistance, is the most investigated eukaryotic family member. Although a large number of biochemical and structural approaches have provided important information, the conformational dynamics underlying the coupling between ATP binding/hydrolysis and allocrite transport remains elusive. To tackle this issue, we performed molecular dynamic simulations for different nucleotide occupancy states of Sav1866, a prokaryotic P-glycoprotein homologue. The simulations reveal an outward-closed conformation of the transmembrane domain that is stabilized by the binding of two ATP molecules. The hydrolysis of a single ATP leads the X-loop, a key motif of the ATP binding cassette, to interfere with the transmembrane domain and favor its outward-open conformation. Our findings provide a structural basis for the unidirectionality of transport in ABC exporters and suggest a ratio of one ATP hydrolyzed per transport cycle.

Acid-sensing ion channels (ASICs) are proton-activated Na+ channels expressed in the nervous system, where they are involved in learning, fear behaviors, neurodegeneration, and pain sensation. In this work, we study the role in pH sensing of two regions of the ectodomain enriched in acidic residues: the acidic pocket, which faces the outside of the protein and is the binding site of several animal toxins, and the palm, a central channel domain. Using voltage clamp fluorometry, we find that the acidic pocket undergoes conformational changes during both activation and desensitization. Concurrently, we find that, although proton sensing in the acidic pocket is not required for channel function, it does contribute to both activation and desensitization. Furthermore, protonation-mimicking mutations of acidic residues in the palm induce a dramatic acceleration of desensitization followed by the appearance of a sustained current. In summary, this work describes the roles of potential pH sensors in two extracellular domains, and it proposes a model of acidification-induced conformational changes occurring in the acidic pocket of ASIC1a.

X-ray crystallography and mutagenesis data of the prokaryotic Cl-/H+ antiporter ClC-ec1 revealed the pathway for chloride ions and suggested that a conserved glutamate, E148, is determinant for both chloride and proton transport. However, the mechanisms regulating the coupled transport of both substrates remain poorly understood. Using molecular dynamic free energy simulations of ClC-ec1, we found that a conformational change of a conserved phenylalanine residue lining the pore, F357, is a novel and essential component of the transport cycle. Our simulations reveal that, in the presence of an ion bound to the central site of the Cl- permeation pathway, the side chains of F357 and E148 undergo a concerted conformational rearrangement that allows the aromatic ring of the phenylalanine to coordinate the protonated E148 carboxylate group through a π;-dipole interaction. The reorientation of the F357 aromatic ring has two major functional consequences: 1) it favors the binding of a second Cl- to the pore, at the external site; 2) it allows for the controlled diffusion of the proton attached to the E148 side chain along a path distinct from the Cl- pathway. Interestingly, in our molecular dynamic simulations, the introduction of a known inhibitory crosslink between residues A399 and A432 prevents rearrangement of the side chain of F357. In further support for the proposed role of F357 in the CLC exchange cycle, we show using an uptake assay and ITC measurements that the F357A mutation drastically impairs transport of CLC-ec1 with minimal effects on the Cl-/H+ stoichiometry. Our results shed mechanistic light into the molecular basis of coupled Cl-/H+ exchange in the CLCs by identifying a novel conformational rearrangement essential for transport of both Cl- and H+.

Potassium channels are opened by ligands and/or membrane potential. In voltage-gated K+ channels and the prokaryotic KcsA channel, conduction is believed to result from opening of an intracellular constriction that prevents ion entry into the pore. On the other hand, numerous ligand-gated K+ channels lack such gate, suggesting that they may be activated by a change within the selectivity filter, a narrow region at the extracellular side of the pore. Using molecular dynamics simulations and electrophysiology measurements, we show that ligand-induced conformational changes in the KcsA channel removes steric restraints at the selectivity filter, thus resulting in structural fluctuations, reduced K+ affinity, and increased ion permeation. Such activation of the selectivity filter may be a universal gating mechanism within K+ channels. The occlusion of the pore at the level of the intracellular gate appears to be secondary.

Flexible polypeptides such as unfolded proteins may access an astronomical number of conformations. The most advanced simulations of such states usually comprise tens of thousands of individual structures. In principle, a comparison of parameters predicted from such ensembles to experimental data provides a measure of their quality. In practice, analyses that go beyond the comparison of unbiased average data have been impossible to carry out on the entirety of such very large ensembles and have, therefore, been restricted to much smaller subensembles and/or nondeterministic algorithms. Here, we show that such very large ensembles, on the order of 10(4) to 10(5) conformations, can be analyzed in full by a maximum entropy fit to experimental average data. Maximizing the entropy of the population weights of individual conformations under experimental χ(2) constraints is a convex optimization problem, which can be solved in a very efficient and robust manner to a unique global solution even for very large ensembles. Since the population weights can be determined reliably, the reweighted full ensemble presents the best model of the combined information from simulation and experiment. Furthermore, since the reduction of entropy due to the experimental constraints is well-defined, its value provides a robust measure of the information content of the experimental data relative to the simulated ensemble and an indication for the density of the sampling of conformational space. The method is applied to the reweighting of a 35,000 frame molecular dynamics trajectory of the nonapeptide EGAAWAASS by extensive NMR (3)J coupling and RDC data. The analysis shows that RDCs provide significantly more information than (3)J couplings and that a discontinuity in the RDC pattern at the central tryptophan is caused by a cluster of helical conformations. Reweighting factors are moderate and consistent with errors in MD force fields of less than 3kT. The required reweighting is larger for an ensemble derived from a statistical coil model, consistent with its coarser nature. We call the method COPER, for convex optimization for ensemble reweighting. Similar advantages of large-scale efficiency and robustness can be obtained for other ensemble analysis methods with convex targets and constraints, such as constrained χ(2) minimization and the maximum occurrence method.

Ligand–receptor interactions that are reinforced by mechanical stress, so-called catch-bonds, play a major role in cell–cell adhesion. They critically contribute to widespread urinary tract infections by pathogenic Escherichia coli strains. These pathogens attach to host epithelia via the adhesin FimH, a two-domain protein at the tip of type I pili recognizing terminal mannoses on epithelial glycoproteins. Here we establish peptide-complemented FimH as a model system for fimbrial FimH function. We reveal a three-state mechanism of FimH catch-bond formation based on crystal structures of all states, kinetic analysis of ligand interaction and molecular dynamics simulations. In the absence of tensile force, the FimH pilin domain allosterically accelerates spontaneous ligand dissociation from the FimH lectin domain by 100,000-fold, resulting in weak affinity. Separation of the FimH domains under stress abolishes allosteric interplay and increases the affinity of the lectin domain. Cell tracking demonstrates that rapid ligand dissociation from FimH supports motility of piliated E. coli on mannosylated surfaces in the absence of shear force.

The ATP binding cassette (ABC) transporters ABCG2 and ABCB1 perform ATP hydrolysis-dependent efflux of structurally highly diverse compounds, collectively called allocrites. Whereas much is known about allocrite-ABCB1 interactions, the chemical nature and strength of ABCG2-allocrite interactions have not yet been assessed. We quantified and characterized interactions of allocrite with ABCG2 and ABCB1 using a set of 39 diverse compounds. We also investigated potential allocrite binding sites based on available transporter structures and structural models. We demonstrate that ABCG2 binds its allocrites from the lipid membrane, despite their hydrophilicity. Hence, binding of allocrite to both transporters is a two-step process, starting with a lipid-water partitioning step, driven mainly by hydrophobic interactions, followed by a transporter binding step in the lipid membrane. We show that binding of allocrite to both transporters increases with the number of hydrogen bond acceptors in allocrites. Scrutinizing the transporter translocation pathways revealed ample hydrogen bond donors for allocrite binding. Importantly, the hydrogen bond donor strength is, on average, higher in ABCG2 than in ABCB1, which explains the higher measured affinity of allocrite for ABCG2. π-π stacking and π-cation interactions play additional roles in binding of allocrite to ABCG2 and ABCB1. With this analysis, we demonstrate that these membrane-mediated weak electrostatic interactions between transporters and allocrites allow for transporter promiscuity toward allocrites. The different sensitivities of the transporters to allocrites' charge and amphiphilicity provide transporter specificity. In addition, we show that the different hydrogen bond donor strengths in the two transporters allow for affinity tuning.

The relation between the sequence of a protein and its three-dimensional structure remains largely unknown. A lasting dream is to elucidate the side-chain-dependent driving forces that govern the folding process. Different structural data suggest that aromatic amino acids play a particular role in the stabilization of protein structures. To better understand the underlying mechanism, we studied peptides of the sequence EGAAXAASS (X = Gly, Ile, Tyr, Trp) through comparison of molecular dynamics (MD) trajectories and NMR residual dipolar coupling (RDC) measurements. The RDC data for aromatic substitutions provide evidence for a kink in the peptide backbone. Analysis of the MD simulations shows that the formation of internal hydrogen bonds underlying a helical turn is key to reproduce the experimental RDC values. The simulations further reveal that the driving force leading to such helical-turn conformations arises from the lack of hydration of the peptide chain on either side of the bulky aromatic side chain, which can potentially act as a nucleation point initiating the folding process.

Oxidative folding in the endoplasmic reticulum (ER) involves ER oxidoreductin 1 (Ero1)-mediated disulfide formation in protein disulfide isomerase (PDI). In this process, Ero1 consumes oxygen (O2) and releases hydrogen peroxide (H2O2), but none of the published Ero1 crystal structures reveal any potential pathway for entry and exit of these reactants. We report that additional mutation of the Cys208–Cys241 disulfide in hyperactive Ero1α (Ero1α-C104A/C131A) potentiates H2O2 production, ER oxidation, and cell toxicity. This disulfide clamps two helices that seal the flavin cofactor where O2 is reduced to H2O2. Through its carboxyterminal active site, PDI unlocks this seal by forming a Cys208/Cys241-dependent mixed-disulfide complex with Ero1α. The H2O2-detoxifying glutathione peroxidase 8 also binds to the Cys208/Cys241 loop region. Supported by O2 diffusion simulations, these data describe the first enzymatically controlled O2 access into a flavoprotein active site, provide molecular-level understanding of Ero1α regulation and H2O2 production/detoxification, and establish the deleterious consequences of constitutive Ero1 activity.

In human cells, membrane proteins of the rhesus (Rh) family excrete ammonium and play a role in pH regulation. Based on high-resolution structures, Rh proteins are generally understood to act as NH3 channels. Given that cell membranes are permeable to gases like NH3, the role of such proteins remains a paradox. Using molecular and quantum mechanical calculations, we show that a crystallographically identified site in the RhCG pore actually recruits NH4+, which is found in higher concentration and binds with higher affinity than NH3, increasing the efficiency of the transport mechanism. A proton is transferred from NH4+ to a signature histidine (the only moiety thermodynamically likely to accept a proton) followed by the diffusion of NH3 down the pore. The excess proton is circulated back to the extracellular vestibule through a hydrogen bond network, which involves a highly conserved and functionally important aspartic acid, resulting in the net transport of NH3.

Ion efflux through potassium channels is regulated by a main gate controlled by ligands and/or the transmembrane voltage. The opening of the main gate let K ions flow down their gradient, diffusing through the selectivity filter of the channel, which can also contribute to the regulation of the ionic current by changing its conformation. Functional and structural data have brought evidence of interaction between the main gate and the selectivity filter of potassium channels. It is understood that opening the main gate favors the inactivated state of the filter, whereas closing the main gate stabilize the conducting state of the filter. We have performed molecular dynamics simulations and potential of mean force calculations showing that the selectivity filter of the KcsA channel in its canonical conducting state, as captured by crystallography, does not allow the diffusion of ions at a rate compatible with experimental data. However, our calculations reveal that the opening of the main gate favors a transient conformational state of the selectivity filter that allows transport rate near the ion diffusion limit. The selectivity filter is thus expected to go from a pre-conducting state with high ion binding affinity to a conducting state with lower ion binding affinity, eventually transiting to the inactivated state. The detailed balance of forces involved in the fine regulation of ion permeation is well described by the latest generation of classical force-fields, which can reveal structural features of ion channels with an unprecedented level of precision.