Assessing the ability of sequence-based methods to provide functional insight within membrane integral proteins: a case study analyzing the neurotransmitter/Na+ symporter family

Resculpting the binding pocket of APC superfamily LeuT-fold amino acid transporters

Amino acid transporters are essential components of prokaryote and eukaryote cells, possess distinct physiological functions, and differ markedly in substrate specificity. Amino acid transporters can be both drug targets and drug transporters (bioavailability, targeting) with many monogenic disorders resulting from dysfunctional membrane transport. The largest collection of amino acid transporters (including the mammalian SLC6, SLC7, SLC32, SLC36, and SLC38 families), across all kingdoms of life, is within the Amino acid-Polyamine-organoCation (APC) superfamily. The LeuT-fold is a paradigm structure for APC superfamily amino acid transporters and carriers of sugars, neurotransmitters, electrolytes, osmolytes, vitamins, micronutrients, signalling molecules, and organic and fatty acids. Each transporter is specific for a unique sub-set of solutes, specificity being determined by how well a substrate fits into each binding pocket. However, the molecular basis of substrate selectivity remains, by and large, elusive. Using an integrated computational and experimental approach, we demonstrate that a single position within the LeuT-fold can play a crucial role in determining substrate specificity in mammalian and arthropod amino acid transporters within the APC superfamily. Systematic mutation of the amino acid residue occupying the equivalent position to LeuT V104 titrates binding pocket space resulting in dramatic changes in substrate selectivity in exemplar APC amino acid transporters including PAT2 (SLC36A2) and SNAT5 (SLC38A5). Our work demonstrates how a single residue/site within an archetypal structural motif can alter substrate affinity and selectivity within this important superfamily of diverse membrane transporters.

Revised Ion/Substrate Coupling Stoichiometry of GABA Transporters

The purpose of this review is to highlight recent evidence in support of a 3 Na⁺: 1 Cl^-: 1 GABA coupling stoichiometry for plasma membrane GABA transporters (SLC6A1 , SLC6A11 , SLC6A12 , SLC6A13 ) and how the revised stoichiometry impacts our understanding of the contribution of GABA transporters to GABA homeostasis in synaptic and extrasynaptic regions in the brain under physiological and pathophysiological states. Recently, our laboratory probed the GABA transporter stoichiometry by analyzing the results of six independent measurements, which included the shifts in the thermodynamic transporter reversal potential caused by changes in the extracellular Na⁺, Cl^-, and GABA concentrations, as well as the ratio of charge flux to substrate flux for Na⁺, Cl^-, and GABA under voltage-clamp conditions. The shifts in the transporter reversal potential for a tenfold change in the external concentration of Na⁺, Cl^-, and GABA were 84 ± 4, 30 ± 1, and 29 ± 1 mV, respectively. Charge flux to substrate flux ratios were 0.7 ± 0.1 charges/Na⁺, 2.0 ± 0.2 charges/Cl^-, and 2.1 ± 0.1 charges/GABA. We then compared these experimental results with the predictions of 150 different transporter stoichiometry models, which included 1-5 Na⁺, 0-5 Cl^-, and 1-5 GABA per transport cycle. Only the 3 Na⁺: 1 Cl^-: 1 GABA stoichiometry model correctly predicts the results of all six experimental measurements. Using the revised 3 Na⁺: 1 Cl^-: 1 GABA stoichiometry, we propose that the GABA transporters mediate GABA uptake under most physiological conditions. Transporter-mediated GABA release likely takes place under pathophysiological or extreme physiological conditions.

Functional consequences of sulfhydryl modification of the γ-aminobutyric acid transporter 1 at a single solvent-exposed cysteine residue

The aims of this study were to optimize the experimental conditions for labeling extracellularly oriented, solvent-exposed cysteine residues of γ-aminobutyric acid transporter 1 (GAT1) with the membrane-impermeant sulfhydryl reagent [2-(trimethylammonium)ethyl]methanethiosulfonate (MTSET) and to characterize the functional and pharmacological consequences of labeling on transporter steady-state and presteady-state kinetic properties. We expressed human GAT1 in Xenopus laevis oocytes and used radiotracer and electrophysiological methods to assay transporter function before and after sulfhydryl modification with MTSET. In the presence of NaCl, transporter exposure to MTSET (1–2.5 mM for 5–20 min) led to partial inhibition of GAT1-mediated transport, and this loss of function was completely reversed by the reducing reagent dithiothreitol. MTSET treatment had no functional effect on the mutant GAT1 C74A, whereas the membrane-permeant reagents N-ethylmaleimide and tetramethylrhodamine-6-maleimide inhibited GABA transport mediated by GAT1 C74A. Ion replacement experiments indicated that MTSET labeling of GAT1 could be driven to completion when valproate replaced chloride in the labeling buffer, suggesting that valproate induces a GAT1 conformation that significantly increases C74 accessibility to the extracellular fluid. Following partial inhibition by MTSET, there was a proportional reduction in both the presteady-state and steady-state macroscopic signals, and the functional and pharmacological properties of the remaining signals were indistinguishable from those of unlabeled GAT1. Therefore, covalent modification of GAT1 at C74 results in completely nonfunctional as well as electrically silent transporters.

Topology improves phylogenetic motif functional site predictions

Prediction of protein functional sites from sequence-derived data remains an open bioinformatics problem. We have developed a phylogenetic motif (PM) functional site prediction approach that identifies functional sites from alignment fragments that parallel the evolutionary patterns of the family. In our approach, PMs are identified by comparing tree topologies of each alignment fragment to that of the complete phylogeny. Herein, we bypass the phylogenetic reconstruction step and identify PMs directly from distance matrix comparisons. In order to optimize the new algorithm, we consider three different distance matrices and 13 different matrix similarity scores. We assess the performance of the various approaches on a structurally nonredundant data set that includes three types of functional site definitions. Without exception, the predictive power of the original approach outperforms the distance matrix variants. While the distance matrix methods fail to improve upon the original approach, our results are important because they clearly demonstrate that the improved predictive power is based on the topological comparisons. Meaning that phylogenetic trees are a straightforward, yet powerful way to improve functional site prediction accuracy. While complementary studies have shown that topology improves predictions of protein-protein interactions, this report represents the first demonstration that trees improve functional site predictions as well.

GATMD: γ-aminobutyric acid transporter mutagenesis database

Since the cloning of the first γ-aminobutyric acid (GABA) transporter (GAT1; SLC6A1) from rat brain in 1990, more than 50 published studies have provided structure-function information on investigator-designed rat and mouse GAT1 mutants. To date, more than 200 of 599 GAT1 residues have been subjected to mutagenesis experiments by substitution with different amino acids, and the resulting transporter functional properties have significantly advanced our understanding of the mechanism of Na+- and Cl⁻-coupled GABA transport by this important member of the neurotransmitter:sodium symporter family. Moreover, many studies have addressed the functional consequences of amino acid deletion or insertion at various positions along the primary sequence. The enormity of this growing body of structure-function information has prompted us to develop GABA Transporter Mutagenesis Database (GATMD), a web-accessible, relational database of manually annotated biochemical, functional and pharmacological data reported on GAT1-the most intensely studied GABA transporter isoform. As of the last update of GATMD, 52 GAT1 mutagenesis papers have yielded 3360 experimental records, which collectively contain a total of ∼100 000 annotated parameters. Database URL: http://physiology.sci.csupomona.edu/GATMD/

Improving position-specific predictions of protein functional sites using phylogenetic motifs

MOTIVATION

Accurate computational prediction of protein functional sites is critical to maximizing the utility of recent high-throughput sequencing efforts. Among the available approaches, position-specific conservation scores remain among the most popular due to their accuracy and ease of computation. Unfortunately, high false positive rates remain a limiting factor. Using phylogenetic motifs (PMs), we have developed two combined (conservation + PMs) prediction schemes that significantly improve prediction accuracy.

RESULTS

Our first approach, called position-specific MINER (psMINER), rank orders alignment columns by conservation. Subsequently, positions that are also not identified as PMs are excluded from the prediction set. This approach improves prediction accuracy, in a statistically significant way, compared to the underlying conservation scores. Increased accuracy is a general result, meaning improvement is observed over several different conservation scores that span a continuum of complexity. In addition, a hybrid MINER (hMINER) that quantitatively considers both scoring regimes provides further improvement. More importantly, it provides critical insight into the relative importance of phylogeny versus alignment conservation. Both methods outperform other common prediction algorithms that also utilize phylogenetic concepts. Finally, we demonstrate that the presented results are critically sensitive to functional site definition, thus highlighting the need for more complete benchmarks within the prediction community.