Substrate-bound outward-open state of the betaine transporter BetP provides insights into Na+ coupling

MolDy: molecular dynamics simulation made easy

MOTIVATION

Molecular dynamics (MD) is a computational experiment that is crucial for understanding the structure of biological macro and micro molecules, their folding, and the inter-molecular interactions. Accurate knowledge of these structural features is the cornerstone in drug development and elucidating macromolecules functions. The open-source GROMACS biomolecular MD simulation program is recognized as a reliable and frequently used simulation program for its precision. However, the user requires expertise, and scripting skills to carrying out MD simulations.

RESULTS

We have developed an end-to-end interactive MD simulation application, MolDy for Gromacs. This front-end application provides a customizable user interface integrated with the Python and Perl-based logical backend connecting the Linux shell and Gromacs software. The tool performs analysis and provides the user with simulation trajectories and graphical representations of relevant biophysical parameters. The advantages of MolDy are (i) user-friendly, does not requiring the researcher to have prior knowledge of Linux; (ii) easy installation by a single command; (iii) freely available for academic research; (iv) can run with minimum configuration of operating systems; (v) has valid default prefilled parameters for beginners, and at the same time provides scope for modifications for expert users.

AVAILABILITY AND IMPLEMENTATION

MolDy is available freely as compressed source code files with user manual for installation and operation on GitHub: https://github.com/AIBResearchMolDy/Moldyv01.git and on https://aibresearch.com/innovations.

Identification of the potassium-binding site in serotonin transporter

Clearance of serotonin (5-hydroxytryptamine, 5-HT) from the synaptic cleft after neuronal signaling is mediated by serotonin transporter (SERT), which couples this process to the movement of a Na+ ion down its chemical gradient. After release of 5-HT and Na+ into the cytoplasm, the transporter faces a rate-limiting challenge of resetting its conformation to be primed again for 5-HT and Na+ binding. Early studies of vesicles containing native SERT revealed that K+ gradients can provide an additional driving force, via K+ antiport. Moreover, under appropriate conditions, a H+ ion can replace K+. Intracellular K+ accelerates the resetting step. Structural studies of SERT have identified two binding sites for Na+ ions, but the K+ site remains enigmatic. Here, we show that K+ antiport can drive substrate accumulation into vesicles containing SERT extracted from a heterologous expression system, allowing us to study the residues responsible for K+ binding. To identify candidate binding residues, we examine many cation binding configurations using molecular dynamics simulations, predicting that K+ binds to the so-called Na2 site. Site-directed mutagenesis of residues in this site can eliminate the ability of both K+ and H+ to drive 5-HT accumulation into vesicles and, in patch clamp recordings, prevent the acceleration of turnover rates and the formation of a channel-like state by K+ or H+. In conclusion, the Na2 site plays a pivotal role in orchestrating the sequential binding of Na+ and then K+ (or H+) ions to facilitate 5-HT uptake in SERT.

Conformational transition induced in the aspartate:alanine antiporter by L-Ala binding

An aspartate:alanine antiporter (AspT) from the lactic acid bacterium Tetragenococcus halophilus catalyzes the electrogenic aspartate^1-:alanine⁰ exchange reaction. Our previous kinetic analyses of transport reactions mediated by AspT in reconstituted liposomes suggested that, although the substrate transport reactions are physiologically coupled, the putative binding sites of L-aspartate (-Asp) and L-alanine (-Ala) are independently located on AspT. By using the fluorescent probe Oregon Green maleimide (OGM), which reacts specifically with cysteine, we also found that the presence of L-Asp changes the conformation of AspT. In this study, we conducted an OGM labeling assay in the presence of L-Ala. The labeling efficiency of single cysteine mutants (G62C and P79C) in transmembrane helix 3 of the AspT showed novel patterns depending on the presence of L-Ala or analogs. A concentration-dependent shift of AspT from the conformation in the presence of one substrate to that specific to the substrate added subsequently (L-Ala or L-Asp) was observed. Moreover, size-exclusion-chromatography-based thermostability assays indicated that the thermal stability of AspT in the presence of L-Ala differed from that in the presence of L-Asp. From these results, we concluded that L-Ala binding yields a conformation different from the apo or L-Asp binding conformations.

Hyperosmotic stress allosterically reconfigures betaine binding pocket in BetP

The transporter BetP in C. glutamicum is essential in maintaining bacterial cell viability during hyperosmotic stress and functions by co-transporting betaine and Na⁺ into bacterial cells. Hyperosmotic stress leads to increased intracellular K⁺ concentrations which in turn promotes betaine binding. While structural details of multiple end state conformations of BetP have provided high resolution snapshots, how K⁺ sensing by the C-terminal domain is allosterically relayed to the betaine binding site is not well understood. In this study, we describe conformational dynamics in solution of BetP using amide hydrogen/deuterium exchange mass spectrometry (HDXMS). These reveal how K⁺ alters conformation of the disordered C- and N-terminal domains to allosterically reconfigure transmembrane helices 3,8 and 10 (TM 3, 8, 10) to enhance betaine interactions. A map of the betaine binding site, at near single amino acid resolution, reveals a critical extrahelical H-bond mediated by TM3 with betaine.

Principles of Alternating Access in LeuT-fold Transporters: Commonalities and Divergences

Found in all domains of life, transporters belonging to the LeuT-fold class mediate the import and exchange of hydrophilic and charged compounds such as amino acids, metals, and sugar molecules. Nearly two decades of investigations on the eponymous bacterial transporter LeuT have yielded a library of high-resolution snapshots of its conformational cycle linked by solution-state experimental data obtained from multiple techniques. In parallel, its topology has been observed in symporters and antiporters characterized by a spectrum of substrate specificities and coupled to gradients of distinct ions. Here we review and compare mechanistic models of transport for LeuT, its well-studied homologs, as well as functionally distant members of the fold, emphasizing the commonalities and divergences in alternating access and the corresponding energy landscapes. Our integrated summary illustrates how fold conservation, a hallmark of the LeuT fold, coincides with divergent choreographies of alternating access that nevertheless capitalize on recurrent structural motifs. In addition, it highlights the knowledge gap that hinders the leveraging of the current body of research into detailed mechanisms of transport for this important class of membrane proteins.

Na+ homeostasis in Acinetobacter baumannii is facilitated via the activity of the Mrp antiporter

The human opportunistic pathogen Acinetobacter baumannii is a global threat to healthcare institutions worldwide, since it developed very efficient strategies to evade host defense and to adapt to the different environmental conditions of the host. This worked focused on the importance of Na⁺ homeostasis in A. baumannii with regards to pathobiological aspects. In silico studies revealed a homologue of a multicomponent Na⁺ /H⁺ antiporter system. Inactivation of the Mrp antiporter through deletion of the first gene (mrpA') resulted in a mutant that was sensitive to increasing pH values. Furthermore, the strain was highly sensitive to increasing Na⁺ and Li⁺ concentrations. Increasing Na⁺ sensitivity is thought to be responsible for growth impairment in human fluids. Furthermore, deletion of mrpA' is associated with energetic defects, inhibition of motility and survival under anoxic and dry conditions.

Mechanistic basis of choline import involved in teichoic acids and lipopolysaccharide modification

Phosphocholine molecules decorating bacterial cell wall teichoic acids and outer-membrane lipopolysaccharide have fundamental roles in adhesion to host cells, immune evasion, and persistence. Bacteria carrying the operon that performs phosphocholine decoration synthesize phosphocholine after uptake of the choline precursor by LicB, a conserved transporter among divergent species. Streptococcus pneumoniae is a prominent pathogen where phosphocholine decoration plays a fundamental role in virulence. Here, we present cryo-electron microscopy and crystal structures of S. pneumoniae LicB, revealing distinct conformational states and describing architectural and mechanistic elements essential to choline import. Together with in vitro and in vivo functional characterization, we found that LicB displays proton-coupled import activity and promiscuous selectivity involved in adaptation to choline deprivation conditions, and describe LicB inhibition by synthetic nanobodies (sybodies). Our results provide previously unknown insights into the molecular mechanism of a key transporter involved in bacterial pathogenesis and establish a basis for inhibition of the phosphocholine modification pathway across bacterial phyla.

General principles of secondary active transporter function

Transport of ions and small molecules across the cell membrane against electrochemical gradients is catalyzed by integral membrane proteins that use a source of free energy to drive the energetically uphill flux of the transported substrate. Secondary active transporters couple the spontaneous influx of a "driving" ion such as Na+ or H+ to the flux of the substrate. The thermodynamics of such cyclical non-equilibrium systems are well understood, and recent work has focused on the molecular mechanism of secondary active transport. The fact that these transporters change their conformation between an inward-facing and outward-facing conformation in a cyclical fashion, called the alternating access model, is broadly recognized as the molecular framework in which to describe transporter function. However, only with the advent of high resolution crystal structures and detailed computer simulations, it has become possible to recognize common molecular-level principles between disparate transporter families. Inverted repeat symmetry in secondary active transporters has shed light onto how protein structures can encode a bi-stable two-state system. Based on structural data, three broad classes of alternating access transitions have been described as rocker-switch, rocking-bundle, and elevator mechanisms. More detailed analysis indicates that transporters can be understood as gated pores with at least two coupled gates. These gates are not just a convenient cartoon element to illustrate a putative mechanism but map to distinct parts of the transporter protein. Enumerating all distinct gate states naturally includes occluded states in the alternating access picture and also suggests what kind of protein conformations might be observable. By connecting the possible conformational states and ion/substrate bound states in a kinetic model, a unified picture emerges in which the symporter, antiporter, and uniporter functions are extremes in a continuum of functionality. As usual with biological systems, few principles and rules are absolute and exceptions are discussed as well as how biological complexity may be integrated in quantitative kinetic models that may provide a bridge from the structure to function.

Optimizing the Substrate Uptake Rate of Solute Carriers

The diversity in solute carriers arose from evolutionary pressure. Here, we surmised that the adaptive search for optimizing the rate of substrate translocation was also shaped by the ambient extracellular and intracellular concentrations of substrate and co-substrate(s). We explored possible solutions by employing kinetic models, which were based on analytical expressions of the substrate uptake rate, that is, as a function of the microscopic rate constants used to parameterize the transport cycle. We obtained the defining terms for five reaction schemes with identical transport stoichiometry (i.e., Na+: substrate = 2:1). We then utilized an optimization algorithm to find the set of numeric values for the microscopic rate constants, which provided the largest value for the substrate uptake rate: The same optimized rate was achieved by different sets of numerical values for the microscopic rate constants. An in-depth analysis of these sets provided the following insights: (i) In the presence of a low extracellular substrate concentration, a transporter can only cycle at a high rate, if it has low values for both, the Michaelis-Menten constant (KM) for substrate and the maximal substrate uptake rate (Vmax). (ii) The opposite is true for a transporter operating at high extracellular substrate concentrations. (iii) Random order of substrate and co-substrate binding is superior to sequential order, if a transporter is to maintain a high rate of substrate uptake in the presence of accumulating intracellular substrate. Our kinetic models provide a framework to understand how and why the transport cycles of closely related transporters differ.

Biomolecular Modeling and Simulation: A Prospering Multidisciplinary Field

We reassess progress in the field of biomolecular modeling and simulation, following up on our perspective published in 2011. By reviewing metrics for the field's productivity and providing examples of success, we underscore the productive phase of the field, whose short-term expectations were overestimated and long-term effects underestimated. Such successes include prediction of structures and mechanisms; generation of new insights into biomolecular activity; and thriving collaborations between modeling and experimentation, including experiments driven by modeling. We also discuss the impact of field exercises and web games on the field's progress. Overall, we note tremendous success by the biomolecular modeling community in utilization of computer power; improvement in force fields; and development and application of new algorithms, notably machine learning and artificial intelligence. The combined advances are enhancing the accuracy andscope of modeling and simulation, establishing an exemplary discipline where experiment and theory or simulations are full partners.

Heteromeric Solute Carriers: Function, Structure, Pathology and Pharmacology

Solute carriers form one of three major superfamilies of membrane transporters in humans, and include uniporters, exchangers and symporters. Following several decades of molecular characterisation, multiple solute carriers that form obligatory heteromers with unrelated subunits are emerging as a distinctive principle of membrane transporter assembly. Here we comprehensively review experimentally established heteromeric solute carriers: SLC3-SLC7 amino acid exchangers, SLC16 monocarboxylate/H⁺ symporters and basigin/embigin, SLC4A1 (AE1) and glycophorin A exchanger, SLC51 heteromer Ost α-Ost β uniporter, and SLC6 heteromeric symporters. The review covers the history of the heteromer discovery, transporter physiology, structure, disease associations and pharmacology - all with a focus on the heteromeric assembly. The cellular locations, requirements for complex formation, and the functional role of dimerization are extensively detailed, including analysis of the first complete heteromer structures, the SLC7-SLC3 family transporters LAT1-4F2hc, b^0,+AT-rBAT and the SLC6 family heteromer B⁰AT1-ACE2. We present a systematic analysis of the structural and functional aspects of heteromeric solute carriers and conclude with common principles of their functional roles and structural architecture.

Mitigation of Hyper KCl Stress at 42ºC with Externally Existing Sodium Glutamate to a Halotolerant Brevibacterium sp. JCM 6894

Halotolerant Brevibacterium sp. JCM 6894 grew at 37ºC in the presence of 2.3 M KCl, while the growth was repressed with the same concentration of NaCl. When resting cells, 10^{7.4 ± 0.1} (CFU·mL^-1), prepared from cells grown in the absence of salts at 30ºC, were exposed to 3.3 M NaCl for 36 h at 42ºC, reduction of the number of resting cells was maintained within a 1-log cycle in the presence of proline, betaine, or ectoine (50 mM). In the presence of 3.3 M KCl, the most functional osmoprotectant was sodium glutamate (50 mM), and the value was 10^{7.2 ± 0.1} (CFU·mL^-1) when exposed for 72 h at 42ºC. In the absence of osmoprotectants, the value was reduced to four orders of magnitude in each experimental condition. The number of resting cells, 10^{6.8 ± 0.1} (CFU·mL^-1), prepared from grown cells pre-adapted to 2.3 M KCl at 37ºC, was hardly reduced when exposed to 3.3 M KCl in the presence of sodium glutamate more than 50 mM for 72 h at 42ºC. Those results indicate that the isolate can sense the difference in hyper KCl stress as opposed to hyper NaCl stress, and different kinds of osmoadaptation systems can function to cope with each hyper salt stress.

Sodium-coupled glucose transport, the SLC5 family, and therapeutically relevant inhibitors: from molecular discovery to clinical application

Sodium glucose transporters (SGLTs) belong to the mammalian solute carrier family SLC5. This family includes 12 different members in human that mediate the transport of sugars, vitamins, amino acids, or smaller organic ions such as choline. The SLC5 family belongs to the sodium symporter family (SSS), which encompasses transporters from all kingdoms of life. It furthermore shares similarity to the structural fold of the APC (amino acid-polyamine-organocation) transporter family. Three decades after the first molecular identification of the intestinal Na+-glucose cotransporter SGLT1 by expression cloning, many new discoveries have evolved, from mechanistic analysis to molecular genetics, structural biology, drug discovery, and clinical applications. All of these advances have greatly influenced physiology and medicine. While SGLT1 is essential for fast absorption of glucose and galactose in the intestine, the expression of SGLT2 is largely confined to the early part of the kidney proximal tubules, where it reabsorbs the bulk part of filtered glucose. SGLT2 has been successfully exploited by the pharmaceutical industry to develop effective new drugs for the treatment of diabetic patients. These SGLT2 inhibitors, termed gliflozins, also exhibit favorable nephroprotective effects and likely also cardioprotective effects. In addition, given the recent finding that SGLT2 is also expressed in tumors of pancreas and prostate and in glioblastoma, this opens the door to potential new therapeutic strategies for cancer treatment by specifically targeting SGLT2. Likewise, further discoveries related to the functional association of other SGLTs of the SLC5 family to human pathologies will open the door to potential new therapeutic strategies. We furthermore hope that the herein summarized information about the physiological roles of SGLTs and the therapeutic benefits of the gliflozins will be useful for our readers to better understand the molecular basis of the beneficial effects of these inhibitors, also in the context of the tubuloglomerular feedback (TGF), and the renin-angiotensin system (RAS). The detailed mechanisms underlying the clinical benefits of SGLT2 inhibition by gliflozins still warrant further investigation that may serve as a basis for future drug development.

Changes in the Survivability of Marine Vibrio sp. under Hyper KCl Stress in the Presence of Betaine as Well as with Exposure to 37ºC

Survivability at hyper KCl stress was examined at 30ºC and 37ºC in the presence and absence of an osmoprotectant by using resting cells prepared from marine Vibrio sp. grown at early stationary phase. Survivability was decided by counting colonies. The number of initial cells, 10^{7.1 ± 0.2} (CFU·mL^-1), was reduced to 10^{5.1 ± 0.5} and < 10^1.0 (CFU·mL^-1) at 30ºC and 37ºC, respectively, by the exposure of resting cells, that were prepared from cells grown for 8 h at 0.5 M NaCl at 30ºC, to 1.2 M KCl and 50 mM NaCl for 3 h. Betaine externally existed as a final concentration of 50 mM mitigated hyper KCl stress to the resting cells at 37ºC. The number of surviving cells was maintained 10^{4.9 ± 0.3} (CFU·mL^-1) when resting cells, 10^{6.5 ± 0.1} (CFU·mL^-1), that were prepared from pre-adapted cells to relatively high concentration of KCl in the growth for 10 h at 0.8 M KCl and 50 mM NaCl at 37ºC, were exposed to 1.2 M KCl, 50 mM NaCl, and 50 mM betaine at 37ºC for 3 h. The results indicate that osmoadaptation system(s) in resting cells is temperature sensitive and betaine functions to mitigate hyper KCl stress to the resting cells at 37ºC.

Functional Characterization of SLC Transporters Using Solid Supported Membranes

Here, we present a protocol for the functional characterization of the H⁺-coupled human peptide transporter PepT1 and sufficient notes to transfer the protocol to the Na⁺-coupled sugar transporter SGLT1, the organic cation transporter OCT2, the Na⁺/Ca²⁺ exchanger NCX, and the neuronal glutamate transporter EAAT3.The assay was developed for the commercially available SURFE²R N1 instrument (Nanion Technologies GmbH) which applies solid supported membrane (SSM)-based electrophysiology. This technique is widely used for the functional characterization of membrane transporters with more than 100 different transporters characterized so far. The technique is cost-effective, easy to use, and capable of high-throughput measurements.SSM-based electrophysiology utilizes SSM-coated gold sensors to physically adsorb membrane vesicles containing the protein of interest. A fast solution exchange provides the substrate and activates transport. For the measurement of PepT1 activity, we applied a peptide concentration jump to activate H⁺/peptide symport. Proton influx charges the sensor. A capacitive current is measured reflecting the transport activity of PepT1 . Multiple measurements on the same sensor allow for comparison of transport activity under different conditions. Here, we determine EC₅₀ for PepT1-mediated glycylglycine transport and perform an inhibition experiment using the specific peptide inhibitor Lys[Z(NO₂)]-Val.

Locking Two Rigid-body Bundles in an Outward-Facing Conformation: The Ion-coupling Mechanism in a LeuT-fold Transporter

Secondary active transporters use electrochemical gradient of ions to fuel the "uphill" translocation of the substrate following the alternating-access model. The coupling of ions to conformational dynamics of the protein remains one of the least characterized aspects of the transporter function. We employ extended molecular dynamics (MD) simulations to examine the Na+-binding effects on the structure and dynamics of a LeuT-fold, Na+-coupled secondary transporter (Mhp1) in its major conformational states, i.e., the outward-facing (OF) and inward-facing (IF) states, as well as on the OF ↔ IF state transition. Microsecond-long, unbiased MD simulations illustrate that Na+ stabilizes an OF conformation favorable for substrate association, by binding to a highly conserved site at the interface between the two helical bundles and restraining their relative position and motion. Furthermore, a special-protocol biased simulation for state transition suggests that Na+ binding hinders the OF ↔ IF transition. These synergistic Na+-binding effects allosterically couple the ion and substrate binding sites and modify the kinetics of state transition, collectively increasing the lifetime of an OF conformation with high substrate affinity, thereby facilitating substrate recruitment from a low-concentration environment. Based on the similarity between our findings for Mhp1 and experimental reports on LeuT, we propose that this model may represent a general Na+-coupling mechanism among LeuT-fold transporters.

Kinetic Models of Secondary Active Transporters

Kinetic models have been employed to understand the logic of substrate transport through transporters of the Solute Carrier (SLC) family. All SLC transporters operate according to the alternate access model, which posits that substrate transport occurs in a closed loop of partial reactions (i.e., a transport cycle). Kinetic models can help to find realistic estimates for conformational transitions between individual states of the transport cycle. When constrained by experimental results, kinetic models can faithfully describe the function of a candidate transporter at a pre-steady state. In addition, we show that kinetic models can accurately predict the intra- and extracellular substrate concentrations maintained by the transporter at a steady state, even under the premise of loose coupling between the electrochemical gradient of the driving ion and of the substrate. We define the criteria for the design of a credible kinetic model of the SLC transporter. Parsimony is the guiding principle of kinetic modeling. We argue, however, that the level of acceptable parsimony is limited by the need to account for the substrate gradient established by a secondary active transporter, and for random order binding of co-substrates and substrate. Random order binding has consistently been observed in transporters of the SLC group.

A comparison of the transport kinetics of glycine transporter 1 and glycine transporter 2

Erdem et al. compare the kinetics of the SLC6 family glycine transporters GlyT1 and GlyT2. Though the two transporters are rate-limited by distinct reaction steps, they both display high transport capacity, with the kinetics of GlyT1 sufficient to supply extracellular glycine to the NMDA receptor.

Transporters of the solute carrier 6 (SLC6) family translocate their cognate substrate together with Na⁺ and Cl⁻. Detailed kinetic models exist for the transporters of GABA (GAT1/SLC6A1) and the monoamines dopamine (DAT/SLC6A3) and serotonin (SERT/SLC6A4). Here, we posited that the transport cycle of individual SLC6 transporters reflects the physiological requirements they operate under. We tested this hypothesis by analyzing the transport cycle of glycine transporter 1 (GlyT1/SLC6A9) and glycine transporter 2 (GlyT2/SLC6A5). GlyT2 is the only SLC6 family member known to translocate glycine, Na⁺, and Cl⁻ in a 1:3:1 stoichiometry. We analyzed partial reactions in real time by electrophysiological recordings. Contrary to monoamine transporters, both GlyTs were found to have a high transport capacity driven by rapid return of the empty transporter after release of Cl⁻ on the intracellular side. Rapid cycling of both GlyTs was further supported by highly cooperative binding of cosubstrate ions and substrate such that their forward transport mode was maintained even under conditions of elevated intracellular Na⁺ or Cl⁻. The most important differences in the transport cycle of GlyT1 and GlyT2 arose from the kinetics of charge movement and the resulting voltage-dependent rate-limiting reactions: the kinetics of GlyT1 were governed by transition of the substrate-bound transporter from outward- to inward-facing conformations, whereas the kinetics of GlyT2 were governed by Na⁺ binding (or a related conformational change). Kinetic modeling showed that the kinetics of GlyT1 are ideally suited for supplying the extracellular glycine levels required for NMDA receptor activation.

Responses of Microorganisms to Osmotic Stress

The cytoplasm of bacterial cells is a highly crowded cellular compartment that possesses considerable osmotic potential. As a result, and owing to the semipermeable nature of the cytoplasmic membrane and the semielastic properties of the cell wall, osmotically driven water influx will generate turgor, a hydrostatic pressure considered critical for growth and viability. Both increases and decreases in the external osmolarity inevitably trigger water fluxes across the cytoplasmic membrane, thus impinging on the degree of cellular hydration, molecular crowding, magnitude of turgor, and cellular integrity. Here, we assess mechanisms that permit the perception of osmotic stress by bacterial cells and provide an overview of the systems that allow them to genetically and physiologically cope with this ubiquitous environmental cue. We highlight recent developments implicating the secondary messenger c-di-AMP in cellular adjustment to osmotic stress and the role of osmotic forces in the life of bacteria-assembled in biofilms.

Structural and molecular aspects of betaine-GABA transporter 1 (BGT1) and its relation to brain function

ɣ-aminobutyric-acid (GABA) functions as the principal inhibitory neurotransmitter in the central nervous system. Imbalances in GABAergic neurotransmission are involved in the pathophysiology of various neurological diseases such as epilepsy, Alzheimer's disease and stroke. GABA transporters (GATs) facilitate the termination of GABAergic signaling by transporting GABA together with sodium and chloride from the synaptic cleft into presynaptic neurons and surrounding glial cells. Four different GATs have been identified that all belong to the solute carrier 6 (SLC6) transporter family: GAT1-3 (SLC6A1, SLC6A13, SLC6A11) and betaine/GABA transporter 1 (BGT1, SLC6A12). BGT1 has emerged as an interesting target for treating epilepsy due to animal studies that reported anticonvulsant effects for the GAT1/BGT1 selective inhibitor EF1502 and the BGT1 selective inhibitor RPC-425. However, the precise involvement of BGT1 in epilepsy remains elusive because of its controversial expression levels in the brain and the lack of highly selective and potent tool compounds. This review gathers the current structural and functional knowledge on BGT1 with emphasis on brain relevance, discusses all available compounds, and tries to shed light on the molecular determinants driving BGT1 selectivity. This article is part of the issue entitled 'Special Issue on Neurotransmitter Transporters'.

Interrogating the conformational dynamics of BetP transport

Robertson reflects on a new study showing how ensemble-biased metadynamics can be used to interpret BetP protein dynamics data.

Interpretation of spectroscopic data using molecular simulations for the secondary active transporter BetP

Mechanistic understanding of dynamic membrane proteins such as transporters, receptors, and channels requires accurate depictions of conformational ensembles, and the manner in which they interchange as a function of environmental factors including substrates, lipids, and inhibitors. Spectroscopic techniques such as electron spin resonance (ESR) pulsed electron-electron double resonance (PELDOR), also known as double electron-electron resonance (DEER), provide a complement to atomistic structures obtained from x-ray crystallography or cryo-EM, since spectroscopic data reflect an ensemble and can be measured in more native solvents, unperturbed by a crystal lattice. However, attempts to interpret DEER data are frequently stymied by discrepancies with the structural data, which may arise due to differences in conditions, the dynamics of the protein, or the flexibility of the attached paramagnetic spin labels. Recently, molecular simulation techniques such as EBMetaD have been developed that create a conformational ensemble matching an experimental distance distribution while applying the minimal possible bias. Moreover, it has been proposed that the work required during an EBMetaD simulation to match an experimentally determined distribution could be used as a metric with which to assign conformational states to a given measurement. Here, we demonstrate the application of this concept for a sodium-coupled transport protein, BetP. Because the probe, protein, and lipid bilayer are all represented in atomic detail, the different contributions to the work, such as the extent of protein backbone movements, can be separated. This work therefore illustrates how ranking simulations based on EBMetaD can help to bridge the gap between structural and biophysical data and thereby enhance our understanding of membrane protein conformational mechanisms.

The GbsR Family of Transcriptional Regulators: Functional Characterization of the OpuAR Repressor

Accumulation of compatible solutes is a common stress response of microorganisms challenged by high osmolarity; it can be achieved either through synthesis or import. These processes have been intensively studied in Bacillus subtilis, where systems for the production of the compatible solutes proline and glycine betaine have been identified, and in which five transporters for osmostress protectants (Opu) have been characterized. Glycine betaine synthesis relies on the import of choline via the substrate-restricted OpuB system and the promiscuous OpuC transporter and its subsequent oxidation by the GbsAB enzymes. Transcription of the opuB and gbsAB operons is under control of the MarR-type regulator GbsR, which acts as an intracellular choline-responsive repressor. Modeling studies using the X-ray structure of the Mj223 protein from Methanocaldococcus jannaschii as the template suggest that GbsR is a homo-dimer with an N-terminal DNA-reading head and C-terminal dimerization domain; a flexible linker connects these two domains. In the vicinity of the linker region, an aromatic cage is predicted as the inducer-binding site, whose envisioned architecture resembles that present in choline and glycine betaine substrate-binding proteins of ABC transporters. We used bioinformatics to assess the phylogenomics of GbsR-type proteins and found that they are widely distributed among Bacteria and Archaea. Alignments of GbsR proteins and analysis of the genetic context of the corresponding structural genes allowed their assignment into four sub-groups. In one of these sub-groups of GbsR-type proteins, gbsR-type genes are associated either with OpuA-, OpuB-, or OpuC-type osmostress protectants uptake systems. We focus here on GbsR-type proteins, named OpuAR by us, that control the expression of opuA-type gene clusters. Using such a system from the marine bacterium Bacillus infantis, we show that OpuAR acts as a repressor of opuA transcription, where several compatible solutes (e.g., choline, glycine betaine, proline betaine) serve as its inducers. Site-directed mutagenesis studies allowed a rational improvement of the putative inducer-binding site in OpuAR with respect to the affinity of choline and glycine betaine binding. Collectively, our data characterize GbsR-/OpuAR-type proteins as an extended sub-group within the MarR-superfamily of transcriptional regulators and identify a novel type of substrate-inducible import system for osmostress protectants.

Molecular Dynamics Simulation for All

The impact of molecular dynamics (MD) simulations in molecular biology and drug discovery has expanded dramatically in recent years. These simulations capture the behavior of proteins and other biomolecules in full atomic detail and at very fine temporal resolution. Major improvements in simulation speed, accuracy, and accessibility, together with the proliferation of experimental structural data, have increased the appeal of biomolecular simulation to experimentalists-a trend particularly noticeable in, although certainly not limited to, neuroscience. Simulations have proven valuable in deciphering functional mechanisms of proteins and other biomolecules, in uncovering the structural basis for disease, and in the design and optimization of small molecules, peptides, and proteins. Here we describe, in practical terms, the types of information MD simulations can provide and the ways in which they typically motivate further experimental work.

Role of the Extremolytes Ectoine and Hydroxyectoine as Stress Protectants and Nutrients: Genetics, Phylogenomics, Biochemistry, and Structural Analysis

Fluctuations in environmental osmolarity are ubiquitous stress factors in many natural habitats of microorganisms, as they inevitably trigger osmotically instigated fluxes of water across the semi-permeable cytoplasmic membrane. Under hyperosmotic conditions, many microorganisms fend off the detrimental effects of water efflux and the ensuing dehydration of the cytoplasm and drop in turgor through the accumulation of a restricted class of organic osmolytes, the compatible solutes. Ectoine and its derivative 5-hydroxyectoine are prominent members of these compounds and are synthesized widely by members of the Bacteria and a few Archaea and Eukarya in response to high salinity/osmolarity and/or growth temperature extremes. Ectoines have excellent function-preserving properties, attributes that have led to their description as chemical chaperones and fostered the development of an industrial-scale biotechnological production process for their exploitation in biotechnology, skin care, and medicine. We review, here, the current knowledge on the biochemistry of the ectoine/hydroxyectoine biosynthetic enzymes and the available crystal structures of some of them, explore the genetics of the underlying biosynthetic genes and their transcriptional regulation, and present an extensive phylogenomic analysis of the ectoine/hydroxyectoine biosynthetic genes. In addition, we address the biochemistry, phylogenomics, and genetic regulation for the alternative use of ectoines as nutrients.

Guardians in a stressful world: the Opu family of compatible solute transporters from Bacillus subtilis

The development of a semi-permeable cytoplasmic membrane was a key event in the evolution of microbial proto-cells. As a result, changes in the external osmolarity will inevitably trigger water fluxes along the osmotic gradient. The ensuing osmotic stress has consequences for the magnitude of turgor and will negatively impact cell growth and integrity. No microorganism can actively pump water across the cytoplasmic membrane; hence, microorganisms have to actively adjust the osmotic potential of their cytoplasm to scale and direct water fluxes in order to prevent dehydration or rupture. They will accumulate ions and physiologically compliant organic osmolytes, the compatible solutes, when they face hyperosmotic conditions to retain cell water, and they rapidly expel these compounds through the transient opening of mechanosensitive channels to curb water efflux when exposed to hypo-osmotic circumstances. Here, we provide an overview on the salient features of the osmostress response systems of the ubiquitously distributed bacterium Bacillus subtilis with a special emphasis on the transport systems and channels mediating regulation of cellular hydration and turgor under fluctuating osmotic conditions. The uptake of osmostress protectants via the Opu family of transporters, systems of central importance for the management of osmotic stress by B. subtilis, will be particularly highlighted.

Human Concentrative Nucleoside Transporter 3 (hCNT3, SLC28A3) Forms a Cyclic Homotrimer

Many anticancer and antiviral drugs are purine or pyrimidine analogues, which use membrane transporters to cross cellular membranes. Concentrative nucleoside transporters (CNTs) mediate the salvage of nucleosides and the transport of therapeutic nucleoside analogues across plasma membranes by coupling the transport of ligands to the sodium gradient. Of the three members of the human CNT family, CNT3 has the broadest selectivity and the widest expression profile. However, the molecular mechanisms of the transporter, including how it interacts with and translocates structurally diverse nucleosides and nucleoside analogues, are unclear. Recently, the crystal structure of vcCNT showed that the prokaryotic homologue of CNT3 forms a homotrimer. In this study, we successfully expressed and purified the wild type human homologue, hCNT3, demonstrating the homotrimer by size exclusion profiles and glutaraldehyde cross-linking. Further, by creating a series of cysteine mutants at highly conserved positions guided by comparative structure models, we cross-linked hCNT3 protomers in a cell-based assay, thus showing the existence of hCNT3 homotrimers in human cells. The presence and absence of cross-links at specific locations along TM9 informs us of important structural differences between vcCNT and hCNT3. Comparative modeling of the trimerization domain and sequence coevolution analysis both indicate that oligomerization is critical to the stability and function of hCNT3. In particular, trimerization appears to shorten the translocation path for nucleosides across the plasma membrane and may allow modulation of the transport function via allostery.

Core Transmembrane Domain 6 Plays a Pivotal Role in the Transport Cycle of the Sodium/Proline Symporter PutP

Crystal structures of transporters with a LeuT-type structural fold assign core transmembrane domain 6 (TM6') a central role in substrate binding and translocation. Here, the function of TM6' in the sodium/proline symporter PutP, a member of the solute/sodium symporter family, was investigated. A complete scan of TM6' identified eight amino acids as particularly important for PutP function. Of these residues, Tyr-248, His-253, and Arg-257 impact sodium binding, whereas Arg-257 and Ala-260 may participate in interactions leading to closure of the inner gate. Furthermore, the previous suggestion of an involvement of Trp-244, Tyr-248, and Pro-252 in proline binding is further supported. In addition, substitution of Gly-245, Gly-247, and Gly-250 affects the amount of PutP in the membrane. A Cys accessibility analysis suggests an involvement of the inner half of TM6' in the formation of a hydrophilic pathway that is open to the inside in the absence of ligands and closed in the presence of sodium and proline. In conclusion, the results demonstrate that TM6' plays a central role in substrate binding and release on the inner side of the membrane also in PutP and extend the knowledge on functionally relevant amino acids in transporters with a LeuT-type structural fold.

The V-motifs facilitate the substrate capturing step of the PTS elevator mechanism

We propose that the alternative crystal forms of outward open UlaA (which are experimental, not simulated, and contain the substrate in the cavity) can be used to interpret/validate the MD results from MalT (the substrate capture step, which involves the mobile second TMSs of the V-motifs, TMSs 2 and 7). Since the crystal contacts are the same between the two alternative crystal forms of outward open UlaA, the striking biological differences noted, including rearranged hydrogen bonds and salt bridge coordination, are not attributable to crystal packing differences. Using transport assays, we identified G58 and G286 as essential for normal vitamin C transport, but the comparison of alternative crystal forms revealed that these residues to unhinge TMS movements from substrate-binding side chains, rendering the mid-TMS regions of homologous TMSs 2 and 7 relatively immobile. While the TMS that is involved in substrate binding in MalT is part of the homologous bundle that holds the two separate halves of the transport assembly (two proteins) together, an unequal effect of the two knockouts was observed for UlaA where both V-motifs are free from such dimer interface interactions.

Unveiling the Mechanism of Arginine Transport through AdiC with Molecular Dynamics Simulations: The Guiding Role of Aromatic Residues

Commensal and pathogenic enteric bacteria have developed several systems to adapt to proton leakage into the cytoplasm resulting from extreme acidic conditions. One such system involves arginine uptake followed by export of the decarboxylated product agmatine, carried out by the arginine/agmatine antiporter (AdiC), which thus works as a virtual proton pump. Here, using classical and targeted molecular dynamics, we investigated at the atomic level the mechanism of arginine transport through AdiC of E. coli. Overall, our MD simulation data clearly demonstrate that global rearrangements of several transmembrane segments are necessary but not sufficient for achieving transitions between structural states along the arginine translocation pathway. In particular, local structural changes, namely rotameric conversions of two aromatic residues, are needed to regulate access to both the outward- and inward-facing states. Our simulations have also enabled identification of a few residues, overwhelmingly aromatic, which are essential to guiding arginine in the course of its translocation. Most of them belong to gating elements whose coordinated motions contribute to the alternating access mechanism. Their conservation in all known E. coli acid resistance antiporters suggests that the transport mechanisms of these systems share common features. Last but not least, knowledge of the functional properties of AdiC can advance our understanding of the members of the amino acid-carbocation-polyamine superfamily, notably in eukaryotic cells.

Identification of a 3rd Na+ Binding Site of the Glycine Transporter, GlyT2

The Na+/Cl- dependent glycine transporters GlyT1 and GlyT2 regulate synaptic glycine concentrations. Glycine transport by GlyT2 is coupled to the co-transport of three Na+ ions, whereas transport by GlyT1 is coupled to the co-transport of only two Na+ ions. These differences in ion-flux coupling determine their respective concentrating capacities and have a direct bearing on their functional roles in synaptic transmission. The crystal structures of the closely related bacterial Na+-dependent leucine transporter, LeuTAa, and the Drosophila dopamine transporter, dDAT, have allowed prediction of two Na+ binding sites in GlyT2, but the physical location of the third Na+ site in GlyT2 is unknown. A bacterial betaine transporter, BetP, has also been crystallized and shows structural similarity to LeuTAa. Although betaine transport by BetP is coupled to the co-transport of two Na+ ions, the first Na+ site is not conserved between BetP and LeuTAa, the so called Na1' site. We hypothesized that the third Na+ binding site (Na3 site) of GlyT2 corresponds to the BetP Na1' binding site. To identify the Na3 binding site of GlyT2, we performed molecular dynamics (MD) simulations. Surprisingly, a Na+ placed at the location consistent with the Na1' site of BetP spontaneously dissociated from its initial location and bound instead to a novel Na3 site. Using a combination of MD simulations of a comparative model of GlyT2 together with an analysis of the functional properties of wild type and mutant GlyTs we have identified an electrostatically favorable novel third Na+ binding site in GlyT2 formed by Trp263 and Met276 in TM3, Ala481 in TM6 and Glu648 in TM10.

The molecular mechanism of secondary sodium symporters elucidated through the lens of the computational microscope

Transport of molecules across cellular membranes is a key biological process for normal cell function. In this review we describe current state-of-the-art knowledge on molecular mechanism of secondary active transporters obtained by molecular simulations studies.

Lipid-Protein Interactions in the Regulated Betaine Symporter BetP Probed by Infrared Spectroscopy

The Na(+)-coupled betaine symporter BetP senses changes in the membrane state and increasing levels of cytoplasmic K(+) during hyperosmotic stress latter via its C-terminal domain and regulates transport activity according to both stimuli. This intriguing sensing and regulation behavior of BetP was intensively studied in the past. It was shown by several biochemical studies that activation and regulation depends crucially on the lipid composition of the surrounding membrane. In fact, BetP is active and regulated only when negatively charged lipids are present. Recent structural studies have revealed binding of phosphatidylglycerol lipids to functional important parts of BetP, suggesting a functional role of lipid interactions. However, a regulatory role of lipid interactions could only be speculated from the snapshot provided by the crystal structure. Here, we investigate the nature of lipid-protein interactions of BetP reconstituted in closely packed two-dimensional crystals of negatively charged lipids and probed at the molecular level with Fourier transform infrared (FTIR) spectroscopy. The FTIR data indicate that K(+) binding weakens the interaction of BetP especially with the anionic lipid head groups. We suggest a regulation mechanism in which lipid-protein interactions, especially with the C-terminal domain and the functional important gating helices transmembrane helice 3 (TMH3) and TMH12, confine BetP to its down-regulated transport state. As BetP is also activated by changes in the physical state of the membrane, our results point toward a more general mechanism of how active transport can be modified by dynamic lipid-protein interactions.

When two turn into one: evolution of membrane transporters from half modules

The recently increasing number of atomic structures for active transporters has not only revealed strong conservation in the architecture of sequence-unrelated transporter families, but also identified a unifying element called the 'inverted repeat topology,' which is found in nearly all transporter folds to date. Indeed, most membrane transporters consist of two or more domains with similar structure, so-called repeats. It is tempting to speculate that transporters have evolved by duplication of one repeat followed by gene fusion and modification events. An intriguing question is, whether recent genes encoding such a 'half-transporter' still exist as independent folding units. Although it seems likely that the evolution of membrane transport proteins, which harbor internal repeats, is linked to these minimal structural building blocks, their identification in the absence of structural data represents a major challenge, as sequence homology is not an issue. In this review we discuss two protein families, the DedA family and the SWEET family, being potential half-transporters and putative ancestors for two of the most abundant secondary transporter families, the MFS family and the LeuT-fold family.

Coherent energy transport in classical nonlinear oscillators: An analogy with the Josephson effect

By means of a simple theoretical model and numerical simulations, we demonstrate the presence of persistent energy currents in a lattice of classical nonlinear oscillators with uniform temperature and chemical potential. In analogy with the well-known Josephson effect, the currents are proportional to the sine of the phase differences between the oscillators. Our results elucidate general aspects of nonequilibrium thermodynamics and point towards a way to practically control transport phenomena in a large class of systems. We apply the model to describe the phase-controlled spin-wave current in a bilayer nanopillar.

Computational characterization of structural dynamics underlying function in active membrane transporters

Active transport of materials across the cellular membrane is one the most fundamental processes in biology. In order to accomplish this task, membrane transporters rely on a wide range of conformational changes spanning multiple time and size scales. These molecular events govern key functional aspects in membrane transporters, namely, coordinated gating motions underlying the alternating access mode of operation, and coupling of uphill transport of substrate to various sources of energy, for example, transmembrane electrochemical gradients and ATP binding and hydrolysis. Computational techniques such as molecular dynamics simulations and free energy calculations have equipped us with a powerful repertoire of biophysical tools offering unparalleled spatial and temporal resolutions that can effectively complement experimental methodologies, and therefore help fill the gap of knowledge in understanding the molecular basis of function in membrane transporters.

Amphetamines, new psychoactive drugs and the monoamine transporter cycle

In monoaminergic neurons, the vesicular transporters and the plasma membrane transporters operate in a relay. Amphetamine and its congeners target this relay to elicit their actions: most amphetamines are substrates, which pervert the relay to elicit efflux of monoamines into the synaptic cleft. However, some amphetamines act as transporter inhibitors. Both compound classes elicit profound psychostimulant effects, which render them liable to recreational abuse. Currently, a surge of new psychoactive substances occurs on a global scale. Chemists bypass drug bans by ingenuous structural variations, resulting in a rich pharmacology. A credible transport model must account for their distinct mode of action and link this to subtle differences in activity and undesired, potentially deleterious effects.

Identification of a second substrate-binding site in solute-sodium symporters

The structure of the sodium/galactose transporter (vSGLT), a solute-sodium symporter (SSS) from Vibrio parahaemolyticus, shares a common structural fold with LeuT of the neurotransmitter-sodium symporter family. Structural alignments between LeuT and vSGLT reveal that the crystallographically identified galactose-binding site in vSGLT is located in a more extracellular location relative to the central substrate-binding site (S1) in LeuT. Our computational analyses suggest the existence of an additional galactose-binding site in vSGLT that aligns to the S1 site of LeuT. Radiolabeled galactose saturation binding experiments indicate that, like LeuT, vSGLT can simultaneously bind two substrate molecules under equilibrium conditions. Mutating key residues in the individual substrate-binding sites reduced the molar substrate-to-protein binding stoichiometry to ~1. In addition, the related and more experimentally tractable SSS member PutP (the Na(+)/proline transporter) also exhibits a binding stoichiometry of 2. Targeting residues in the proposed sites with mutations results in the reduction of the binding stoichiometry and is accompanied by severely impaired translocation of proline. Our data suggest that substrate transport by SSS members requires both substrate-binding sites, thereby implying that SSSs and neurotransmitter-sodium symporters share common mechanistic elements in substrate transport.