Mechanisms of anion conduction by coupled glutamate transporters.

Excitatory amino acid transporters (EAATs) are essential for terminating glutamatergic synaptic transmission. They are not only coupled glutamate/Na(+)/H(+)/K(+) transporters but also function as anion-selective channels. EAAT anion channels regulate neuronal excitability, and gain-of-function mutations in these proteins result in ataxia and epilepsy. We have combined molecular dynamics simulations with fluorescence spectroscopy of the prokaryotic homolog GltPh and patch-clamp recordings of mammalian EAATs to determine how these transporters conduct anions. Whereas outward- and inward-facing GltPh conformations are nonconductive, lateral movement of the glutamate transport domain from intermediate transporter conformations results in formation of an anion-selective conduction pathway. Fluorescence quenching of inserted tryptophan residues indicated the entry of anions into this pathway, and mutations of homologous pore-forming residues had analogous effects on GltPh simulations and EAAT2/EAAT4 measurements of single-channel currents and anion/cation selectivities. These findings provide a mechanistic framework of how neurotransmitter transporters can operate as anion-selective and ligand-gated ion channels.

Accurate prediction of dynamic protein-ligand binding using P-score ranking.

Protein-ligand binding prediction typically relies on docking methodologies and associated scoring functions to propose the binding mode of a ligand in a biological target. Significant challenges are associated with this approach, including the flexibility of the protein-ligand system, solvent-mediated interactions, and associated entropy changes. In addition, scoring functions are only weakly accurate due to the short time required for calculating enthalpic and entropic binding interactions. The workflow described here attempts to address these limitations by combining supervised molecular dynamics with dynamical averaging quantum mechanics fragment molecular orbital. This combination significantly increased the ability to predict the experimental binding structure of protein-ligand complexes independent from the starting position of the ligands or the binding site conformation. We found that the predictive power could be enhanced by combining the residence time and interaction energies as descriptors in a novel scoring function named the P-score. This is illustrated using six different protein-ligand targets as case studies.

Fast prediction of antibiotic permeability through membrane channels using Brownian dynamics.

The efficient permeation across the Gram-negative bacterial membrane is an important step in the overall process of antibacterial action of a molecule and the one that has posed a significant hurdle on the way toward approved antibiotics. Predicting the permeability for a large library of molecules and assessing the effect of different molecular transformations on permeation rates of a given molecule is critical to the development of effective antibiotics. We present a computational approach for obtaining estimates of molecular permeability through a porin channel in a matter of hours using a Brownian dynamics approach. The fast sampling using a temperature acceleration scheme enables the approximate estimation of permeability using the inhomogeneous solubility diffusion model. Although the method is a significant approximation to similar all-atom approaches tested previously, we show that the present approach predicts permeabilities that correlate fairly well with the respective experimental permeation rates from liposome swelling experiments and accumulation rates from antibiotic accumulation assays, and is significantly, i.e., about 14 times, faster compared with a previously reported approach. The possible applications of the scheme in high-throughput screening for fast permeators are discussed.

GPCRmd uncovers the dynamics of the 3D-GPCRome.

G-protein-coupled receptors (GPCRs) are involved in numerous physiological processes and are the most frequent targets of approved drugs. The explosion in the number of new three-dimensional (3D) molecular structures of GPCRs (3D-GPCRome) over the last decade has greatly advanced the mechanistic understanding and drug design opportunities for this protein family. Molecular dynamics (MD) simulations have become a widely established technique for exploring the conformational landscape of proteins at an atomic level. However, the analysis and visualization of MD simulations require efficient storage resources and specialized software. Here we present GPCRmd (http://gpcrmd.org/), an online platform that incorporates web-based visualization capabilities as well as a comprehensive and user-friendly analysis toolbox that allows scientists from different disciplines to visualize, analyze and share GPCR MD data. GPCRmd originates from a community-driven effort to create an open, interactive and standardized database of GPCR MD simulations.

Publisher Correction: GPCRmd uncovers the dynamics of the 3D-GPCRome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Disease related single point mutations alter the global dynamics of a tetratricopeptide (TPR) α-solenoid domain.

Tetratricopeptide repeat (TPR) proteins belong to the class of α-solenoid proteins, in which repetitive units of α-helical hairpin motifs stack to form superhelical, often highly flexible structures. TPR domains occur in a wide variety of proteins, and perform key functional roles including protein folding, protein trafficking, cell cycle control and post-translational modification. Here, we look at the TPR domain of the enzyme O-linked GlcNAc-transferase (OGT), which catalyses O-GlcNAcylation of a broad range of substrate proteins. A number of single-point mutations in the TPR domain of human OGT have been associated with the disease Intellectual Disability (ID). By extended steered and equilibrium atomistic simulations, we show that the OGT-TPR domain acts as an elastic nanospring, and that each of the ID-related local mutations substantially affect the global dynamics of the TPR domain. Since the nanospring character of the OGT-TPR domain is key to its function in binding and releasing OGT substrates, these changes of its biomechanics likely lead to defective substrate interaction. We find that neutral mutations in the human population, selected by analysis of the gnomAD database, do not incur these changes. Our findings may not only help to explain the ID phenotype of the mutants, but also aid the design of TPR proteins with tailored biomechanical properties.

The lipid environment determines the activity of the Escherichia coli ammonium transporter AmtB.

The movement of ammonium across biologic membranes is a fundamental process in all living organisms and is mediated by the ubiquitous ammonium transporter/methylammonium permease/rhesus protein (Amt/Mep/Rh) family of transporters. Recent structural analysis and coupled mass spectrometry studies have shown that the Escherichia coli ammonium transporter AmtB specifically binds 1-palmitoyl-2-oleoyl phosphatidylglycerol (POPG). Upon POPG binding, several residues of AmtB undergo a small conformational change, which stabilizes the protein against unfolding. However, no studies have so far been conducted, to our knowledge, to explore whether POPG binding to AmtB has functional consequences. Here, we used an in vitro experimental assay with purified components, together with molecular dynamics simulations, to characterize the relation between POPG binding and AmtB activity. We show that the AmtB activity is electrogenic. Our results indicate that the activity, at the molecular level, of Amt in archaebacteria and eubacteria may differ. We also show that POPG is an important cofactor for AmtB activity and that, in the absence of POPG, AmtB cannot complete the full translocation cycle. Furthermore, our simulations reveal previously undiscovered POPG binding sites on the intracellular side of the lipid bilayer between the AmtB subunits. Possible molecular mechanisms explaining the functional role of POPG are discussed.-Mirandela, G. D., Tamburrino, G., Hoskisson, P. A., Zachariae, U., Javelle, A. The lipid environment determines the activity of the Escherichia coli ammonium transporter AmtB.

Interfacial self-assembly of a bacterial hydrophobin.

The majority of bacteria in the natural environment live within the confines of a biofilm. The Gram-positive bacterium Bacillus subtilis forms biofilms that exhibit a characteristic wrinkled morphology and a highly hydrophobic surface. A critical component in generating these properties is the protein BslA, which forms a coat across the surface of the sessile community. We recently reported the structure of BslA, and noted the presence of a large surface-exposed hydrophobic patch. Such surface patches are also observed in the class of surface-active proteins known as hydrophobins, and are thought to mediate their interfacial activity. However, although functionally related to the hydrophobins, BslA shares no sequence nor structural similarity, and here we show that the mechanism of action is also distinct. Specifically, our results suggest that the amino acids making up the large, surface-exposed hydrophobic cap in the crystal structure are shielded in aqueous solution by adopting a random coil conformation, enabling the protein to be soluble and monomeric. At an interface, these cap residues refold, inserting the hydrophobic side chains into the air or oil phase and forming a three-stranded ß-sheet. This form then self-assembles into a well-ordered 2D rectangular lattice that stabilizes the interface. By replacing a hydrophobic leucine in the center of the cap with a positively charged lysine, we changed the energetics of adsorption and disrupted the formation of the 2D lattice. This limited structural metamorphosis represents a previously unidentified environmentally responsive mechanism for interfacial stabilization by proteins.

Insights into the function of ion channels by computational electrophysiology simulations.

Ion channels are of universal importance for all cell types and play key roles in cellular physiology and pathology. Increased insight into their functional mechanisms is crucial to enable drug design on this important class of membrane proteins, and to enhance our understanding of some of the fundamental features of cells. This review presents the concepts behind the recently developed simulation protocol Computational Electrophysiology (CompEL), which facilitates the atomistic simulation of ion channels in action. In addition, the review provides guidelines for its application in conjunction with the molecular dynamics software package GROMACS. We first lay out the rationale for designing CompEL as a method that models the driving force for ion permeation through channels the way it is established in cells, i.e., by electrochemical ion gradients across the membrane. This is followed by an outline of its implementation and a description of key settings and parameters helpful to users wishing to set up and conduct such simulations. In recent years, key mechanistic and biophysical insights have been obtained by employing the CompEL protocol to address a wide range of questions on ion channels and permeation. We summarize these recent findings on membrane proteins, which span a spectrum from highly ion-selective, narrow channels to wide diffusion pores. Finally we discuss the future potential of CompEL in light of its limitations and strengths. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Insights into the ion-coupling mechanism in the MATE transporter NorM-VC.

Bacteria have developed a variety of different mechanisms to defend themselves from compounds that are toxic to them, such as antibiotics. One of these defence mechanisms is the expulsion of drugs or other noxious compounds by multidrug efflux pumps. Multidrug and toxic compound extrusion (MATE) transporters are efflux pumps that extrude metabolic waste and a variety of antibiotics out of the cell, using an ion gradient as energy source. They function via an alternating-access mechanism. When ions bind in the outward facing conformation, a large conformational change to the inward facing conformation is induced, from which the ion is released and the extruded chemical compound is bound. NorM proteins, which are usually coupled to a Na+ gradient, are members of the MATE family. However, for NorM-VC from Vibrio cholerae, it has been shown that this MATE transporter is additionally coupled to protons. How H+ and Na+ binding are coupled mechanistically to enable drug antiport is not well understood. In this study, we use molecular dynamics simulations to illuminate the sequence of ion binding events that enable efflux. Understanding this antiport mechanism is important to support the development of novel compounds that specifically inhibit the functional cycle of NorM transporters.

A role for loop G in the ß1 strand in GABAA receptor activation.

KEY POINTS: The role of the ß1 strand in GABAA receptor function is unclear. It lies anti-parallel to the ß2 strand, which is known to participate in receptor activation. Molecular dynamics simulation revealed solvent accessible residues within the ß1 strand of the GABAA ß3 homopentamer that might be amenable to analysis using the substituted Cys accessibility method. Cys substitutions from Asp43 to Thr47 in the GABAA α1 subunit showed that D43C and T47C reduced the apparent potency of GABA. F45C caused a biphasic GABA concentration-response relationship and increased spontaneous gating. Cys43 and Cys47 were accessible to 2-aminoethyl methanethiosulphonate (MTSEA) modification, whereas Cys45 was not. Both GABA and the allosteric agonist propofol reduced MTSEA modification of Cys43 and Cys47. By contrast, modification of Cys64 in the ß2 strand loop D was impeded by GABA but unaffected by propofol. These data reveal movement of ß1 strand loop G residues during agonist activation of the GABAA receptor. ABSTRACT: The GABAA receptor α subunit ß1 strand runs anti-parallel to the ß2 strand, which contains loop D, known to participate in receptor activation and agonist binding. However, a role for the ß1 strand has yet to be established. We used molecular dynamics simulation to quantify the solvent accessible surface area (SASA) of ß1 strand residues in the GABAA ß3 homopentamer structure. Residues in the complementary interface equivalent to those between Asp43 and Thr47 in the α1 subunit have an alternating pattern of high and low SASA consistent with a ß strand structure. We investigated the functional role of these ß1 strand residues in the α1 subunit by individually replacing them with Cys residues. D43C and T47C substitutions reduced the apparent potency of GABA at α1ß2γ2 receptors by 50-fold and eight-fold, respectively, whereas the F45C substitution caused a biphasic GABA concentration-response relationship and increased spontaneous gating. Receptors with D43C or T47C substitutions were sensitive to 2-aminoethyl methanethiosulphonate (MTSEA) modification. However, GABA-evoked currents mediated by α1(F45C)ß2γ2 receptors were unaffected by MTSEA, suggesting that this residue is inaccessible. Both GABA and the allosteric agonist propofol reduced MTSEA modification of α1(D43C)ß2γ2 and α1(T47C)ß2γ2 receptors, indicating movement of the ß1 strand even during allosteric activation. This is in contrast to α1(F64C)ß2γ2 receptors, where only GABA, but not propofol, reduced MTSEA modification. These findings provide the first functional evidence for movement of the ß1 strand during gating of the receptor and identify residues that are critical for maintaining GABAA receptor function.

Crystal structure and functional mechanism of a human antimicrobial membrane channel.

Multicellular organisms fight bacterial and fungal infections by producing peptide-derived broad-spectrum antibiotics. These host-defense peptides compromise the integrity of microbial cell membranes and thus evade pathways by which bacteria develop rapid antibiotic resistance. Although more than 1,700 host-defense peptides have been identified, the structural and mechanistic basis of their action remains speculative. This impedes the desired rational development of these agents into next-generation antibiotics. We present the X-ray crystal structure as well as solid-state NMR spectroscopy, electrophysiology, and MD simulations of human dermcidin in membranes that reveal the antibiotic mechanism of this major human antimicrobial, found to suppress Staphylococcus aureus growth on the epidermal surface. Dermcidin forms an architecture of high-conductance transmembrane channels, composed of zinc-connected trimers of antiparallel helix pairs. Molecular dynamics simulations elucidate the unusual membrane permeation pathway for ions and show adjustment of the pore to various membranes. Our study unravels the comprehensive mechanism for the membrane-disruptive action of this mammalian host-defense peptide at atomistic level. The results may form a foundation for the structure-based design of peptide antibiotics.

Structural basis for cooperativity of CRM1 export complex formation.

In eukaryotes, the nucleocytoplasmic transport of macromolecules is mainly mediated by soluble nuclear transport receptors of the karyopherin-ß superfamily termed importins and exportins. The highly versatile exportin chromosome region maintenance 1 (CRM1) is essential for nuclear depletion of numerous structurally and functionally unrelated protein and ribonucleoprotein cargoes. CRM1 has been shown to adopt a toroidal structure in several functional transport complexes and was thought to maintain this conformation throughout the entire nucleocytoplasmic transport cycle. We solved crystal structures of free CRM1 from the thermophilic eukaryote Chaetomium thermophilum. Surprisingly, unbound CRM1 exhibits an overall extended and pitched superhelical conformation. The two regulatory regions, namely the acidic loop and the C-terminal α-helix, are dramatically repositioned in free CRM1 in comparison with the ternary CRM1-Ran-Snurportin1 export complex. Single-particle EM analysis demonstrates that, in a noncrystalline environment, free CRM1 exists in equilibrium between extended, superhelical and compact, ring-like conformations. Molecular dynamics simulations show that the C-terminal helix plays an important role in regulating the transition from an extended to a compact conformation and reveal how the binding site for nuclear export signals of cargoes is modulated by different CRM1 conformations. Combining these results, we propose a model for the cooperativity of CRM1 export complex assembly involving the long-range allosteric communication between the distant binding sites of GTP-bound Ran and cargo.

The Bacterial Hydrophobin BslA is a Switchable Ellipsoidal Janus Nanocolloid.

BslA is an amphiphilic protein that forms a highly hydrophobic coat around Bacillus subtilis biofilms, shielding the bacterial community from external aqueous solution. It has a unique structure featuring a distinct partition between hydrophilic and hydrophobic surfaces. This surface property is reminiscent of synthesized Janus colloids. By investigating the behavior of BslA variants at water-cyclohexane interfaces through a set of multiscale simulations informed by experimental data, we show that BslA indeed represents a biological example of an ellipsoidal Janus nanoparticle, whose surface interactions are, moreover, readily switchable. BslA contains a local conformational toggle, which controls its global affinity for, and orientation at, water-oil interfaces. This adaptability, together with single-point mutations, enables the fine-tuning of its solvent and interfacial interactions, and suggests that BslA could be a basis for biotechnological applications.

Conserved structure and domain organization among bacterial Slc26 transporters.

The Slc26 proteins are a ubiquitous superfamily of anion transporters conserved from bacteria to humans, among which four have been identified as human disease genes. Our functional knowledge of this protein family has increased but limited structural information is available. These proteins contain a transmembrane (TM) domain and a C-terminal cytoplasmic sulfate transporter and anti-sigma factor (STAS) domain. In a fundamental step towards understanding the structure/function relationships within the family we have used small-angle neutron scattering (SANS) on two distantly related bacterial homologues to show that there is a common, dimeric and structural architecture among Slc26A transporters. Pulsed electron-electron double resonance (PELDOR) spectroscopy supports the dimeric SANS-derived model. Using chimaeric/truncated proteins we have determined the domain organization: the STAS domains project away from the TM core and are essential for protein stability. We use the SANS-generated envelopes to assess a homology model of the TM core.

Crystallographic analysis of Neisseria meningitidis PorB extracellular loops potentially implicated in TLR2 recognition.

Among all Neisseriae species, Neisseria meningitidis and Neisseria gonorrhoeae are the only human pathogens, causative agents of bacterial meningitis and gonorrhoea, respectively. PorB, a pan-Neisseriae trimeric porin that mediates diffusive transport of essential molecules across the bacterial outer membrane, is also known to activate host innate immunity via Toll-like receptor 2 (TLR2)-mediated signaling. The molecular mechanism of PorB binding to TLR2 is not known, but it has been hypothesized that electrostatic interactions contribute to ligand/receptor binding. Strain-specific sequence variability in the surface-exposed loops of PorB which are potentially implicated in TLR2 binding, may explain the difference in TLR2-mediated cell activation in vitro by PorB homologs from the commensal Neisseriae lactamica and the pathogen N. meningitidis. Here, we report a comparative structural analysis of PorB from N. meningitidis serogroup B strain 8765 (63% sequence homology with PorB from N. meningitidis serogroup W135) and a mutant in which amino acid substitutions in the extracellular loop 7 lead to significantly reduced TLR2-dependent activity in vitro. We observe that this mutation both alters the loop conformation and causes dramatic changes of electrostatic surface charge, both of which may affect TLR2 recognition and signaling.

Position of transmembrane helix 6 determines receptor G protein coupling specificity.

G protein coupled receptors (GPCRs) transmit extracellular signals into the cell by binding and activating different intracellular signaling proteins, such as G proteins (Gαßγ, families Gi, Gs, Gq, G12/13) or arrestins. To address the issue of Gs vs Gi coupling specificity, we carried out molecular dynamics simulations of lipid-embedded active ß2-adrenoceptor (ß2AR*) in complex with C-terminal peptides derived from the key interaction site of Gα (GαCT) as surrogate of Gαßγ. We find that GiαCT and GsαCT exploit distinct cytoplasmic receptor conformations that coexist in the uncomplexed ß2AR*. The slim GiαCT stabilizes a ß2AR* conformation, not accessible to the bulkier GsαCT, which requires a larger TM6 outward tilt for binding. Our results suggest that the TM6 conformational heterogeneity regulates the catalytic activity of ß2AR* toward Gi or Gs.

Identification of a cation transport pathway in Neisseria meningitidis PorB.

Neisseria meningitidis is the main causative agent of bacterial meningitis. In its outer membrane, the trimeric Neisserial porin PorB is responsible for the diffusive transport of essential hydrophilic solutes across the bilayer. Previous molecular dynamics simulations based on the recent crystal structure of PorB have suggested the presence of distinct solute translocation pathways through this channel. Although PorB has been electrophysiologically characterized as anion-selective, cation translocation through nucleotide-bound PorB during pathogenesis is thought to be instrumental for host cell death. As a result, we were particularly interested in further characterizing cation transport through the pore. We combined a structural approach with additional computational analysis. Here, we present two crystal structures of PorB at 2.1 and 2.65 Å resolution. The new structures display additional electron densities around the protruding loop 3 (L3) inside the pore. We show that these electron densities can be identified as monovalent cations, in our case Cs(+), which are tightly bound to the inner channel. Molecular dynamics simulations reveal further ion interactions and the free energy landscape for ions inside PorB. Our results suggest that the crystallographically identified locations of Cs(+) form a cation transport pathway inside the pore. This finding suggests how positively charged ions are translocated through PorB when the channel is inserted into mitochondrial membranes during Neisserial infection, a process which is considered to dissipate the mitochondrial transmembrane potential gradient and thereby induce apoptosis.

Functional dynamics in the voltage-dependent anion channel.

The voltage-dependent anion channel (VDAC), located in the outer mitochondrial membrane, acts as a gatekeeper for the entry and exit of mitochondrial metabolites. Here we reveal functional dynamics of isoform one of VDAC (VDAC1) by a combination of solution NMR spectroscopy, Gaussian network model analysis, and molecular dynamics simulation. Micro- to millisecond dynamics are significantly increased for the N-terminal six ß-strands of VDAC1 in micellar solution, in agreement with increased B-factors observed in the same region in the bicellar crystal structure of VDAC1. Molecular dynamics simulations reveal that a charge on the membrane-facing glutamic acid 73 (E73) accounts for the elevation of N-terminal protein dynamics as well as a thinning of the nearby membrane. Mutation or chemical modification of E73 strongly reduces the micro- to millisecond dynamics in solution. Because E73 is necessary for hexokinase-I-induced VDAC channel closure and inhibition of apoptosis, our results imply that micro- to millisecond dynamics in the N-terminal part of the barrel are essential for VDAC interaction and gating.

A cooperative knock-on mechanism underpins Ca2+-selective cation permeation in TRPV channels.

The selective exchange of ions across cellular membranes is a vital biological process. Ca2+-mediated signaling is implicated in a broad array of physiological processes in cells, while elevated intracellular concentrations of Ca2+ are cytotoxic. Due to the significance of this cation, strict Ca2+ concentration gradients are maintained across the plasma and organelle membranes. Therefore, Ca2+ signaling relies on permeation through selective ion channels that control the flux of Ca2+ ions. A key family of Ca2+-permeable membrane channels is the polymodal signal-detecting transient receptor potential (TRP) ion channels. TRP channels are activated by a wide variety of cues including temperature, small molecules, transmembrane voltage, and mechanical stimuli. While most members of this family permeate a broad range of cations non-selectively, TRPV5 and TRPV6 are unique due to their strong Ca2+ selectivity. Here, we address the question of how some members of the TRPV subfamily show a high degree of Ca2+ selectivity while others conduct a wider spectrum of cations. We present results from all-atom molecular dynamics simulations of ion permeation through two Ca2+-selective and two non-selective TRPV channels. Using a new method to quantify permeation cooperativity based on mutual information, we show that Ca2+-selective TRPV channel permeation occurs by a three-binding site knock-on mechanism, whereas a two-binding site knock-on mechanism is observed in non-selective TRPV channels. Each of the ion binding sites involved displayed greater affinity for Ca2+ over Na+. As such, our results suggest that coupling to an extra binding site in the Ca2+-selective TRPV channels underpins their increased selectivity for Ca2+ over Na+ ions. Furthermore, analysis of all available TRPV channel structures shows that the selectivity filter entrance region is wider for the non-selective TRPV channels, slightly destabilizing ion binding at this site, which is likely to underlie mechanistic decoupling.