Repeat-swap homology modeling of secondary active transporters: updated protocol and prediction of elevator-type mechanisms

Conformational coupling of the sialic acid TRAP transporter HiSiaQM with its substrate binding protein HiSiaP

The tripartite ATP-independent periplasmic (TRAP) transporters use an extra cytoplasmic substrate binding protein (SBP) to transport a wide variety of substrates in bacteria and archaea. The SBP can adopt an open- or closed state depending on the presence of substrate. The two transmembrane domains of TRAP transporters form a monomeric elevator whose function is strictly dependent on the presence of a sodium ion gradient. Insights from experimental structures, structural predictions and molecular modeling have suggested a conformational coupling between the membrane elevator and the substrate binding protein. Here, we use a disulfide engineering approach to lock the TRAP transporter HiSiaPQM from Haemophilus influenzae in different conformational states. The SBP, HiSiaP, is locked in its substrate-bound form and the transmembrane elevator, HiSiaQM, is locked in either its assumed inward- or outward-facing states. We characterize the disulfide-locked constructs and use single-molecule total internal reflection fluorescence (TIRF) microscopy to study their interactions. Our experiments demonstrate that the SBP and the transmembrane elevator are indeed conformationally coupled, meaning that the open and closed state of the SBP recognize specific conformational states of the transporter and vice versa.

High-resolution structure of a mercury cross-linked ZIP metal transporter reveals delicate motions and metal relay for regulated zinc transport

Zrt-/Irt-like protein (ZIP) divalent metal transporters play a central role in maintaining trace element homeostasis. The prototypical ZIP from Bordetella bronchiseptica (BbZIP) is an elevator-type transporter, but the dynamic motions and detailed transport mechanism remain to be elucidated. Here, we report a high-resolution crystal structure of a mercury-crosslinked BbZIP variant at 1.95 Å, revealing an upward rotation of the transport domain in the new inward-facing conformation and a water-filled metal release channel that is divided into two parallel pathways by the previously disordered cytoplasmic loop. Mutagenesis and transport assays indicated that the newly identified high-affinity metal binding site in the primary pathway acts as a "metal sink" to reduce the transport rate. The discovery of a hinge motion around an extracellular axis allowed us to propose a sequential hinge-elevator-hinge movement of the transport domain to achieve alternating access. These findings provide key insights into the transport mechanisms and activity regulation.

Structural insights into the elevator-type transport mechanism of a bacterial ZIP metal transporter

The Zrt-/Irt-like protein (ZIP) family consists of ubiquitously expressed divalent metal transporters critically involved in maintaining systemic and cellular homeostasis of zinc, iron, and manganese. Here, we present a study on a prokaryotic ZIP from Bordetella bronchiseptica (BbZIP) by combining structural biology, evolutionary covariance, computational modeling, and a variety of biochemical assays to tackle the issue of the transport mechanism which has not been established for the ZIP family. The apo state structure in an inward-facing conformation revealed a disassembled transport site, altered inter-helical interactions, and importantly, a rigid body movement of a 4-transmembrane helix (TM) bundle relative to the other TMs. The computationally generated and biochemically validated outward-facing conformation model revealed a slide of the 4-TM bundle, which carries the transport site(s), by approximately 8 Å toward the extracellular side against the static TMs which mediate dimerization. These findings allow us to conclude that BbZIP is an elevator-type transporter.

Structure-based virtual screening and molecular dynamics of potential inhibitors targeting sodium-bile acid co-transporter of carcinogenic liver fluke Clonorchis sinensis

Background

Clonorchis sinensis requires bile acid transporters as this fluke inhabits bile juice-filled biliary ducts, which provide an extreme environment. Clonorchis sinensis sodium-bile acid co-transporter (CsSBAT) is indispensable for the fluke’s survival in the final host, as it circulates taurocholate and prevents bile toxicity in the fluke; hence, it is recognized as a useful drug target.

Methodology and principal findings

In the present study, using structure-based virtual screening approach, we presented inhibitor candidates targeting a bile acid-binding pocket of CsSBAT. CsSBAT models were built using tertiary structure modeling based on a bile acid transporter template (PDB ID: 3zuy and 4n7x) and were applied into AutoDock Vina for competitive docking simulation. First, potential compounds were identified from PubChem (holding more than 100,000 compounds) by applying three criteria: i) interacting more favorably with CsSBAT than with a human homolog, ii) intimate interaction to the inward- and outward-facing conformational states, iii) binding with CsSBAT preferably to natural bile acids. Second, two compounds were identified following the Lipinski’s rule of five. Third, other two compounds of molecular weight higher than 500 Da (Mr > 500 Da) were presumed to efficiently block the transporter via a feasible rational screening strategy. Of these candidates, compound 9806452 exhibited the least hepatotoxicity that may enhance drug-likeness properties.

Conclusions

It is proposed that compound 9806452 act as a potential inhibitor toward CsSBAT and further studies are warranted for drug development process against clonorchiasis.

Structural and mechanistic analysis of a tripartite ATP-independent periplasmic TRAP transporter

Tripartite ATP-independent periplasmic (TRAP) transporters are found widely in bacteria and archaea and consist of three structural domains, a soluble substrate-binding protein (P-domain), and two transmembrane domains (Q- and M-domains). HiSiaPQM and its homologs are TRAP transporters for sialic acid and are essential for host colonization by pathogenic bacteria. Here, we reconstitute HiSiaQM into lipid nanodiscs and use cryo-EM to reveal the structure of a TRAP transporter. It is composed of 16 transmembrane helices that are unexpectedly structurally related to multimeric elevator-type transporters. The idiosyncratic Q-domain of TRAP transporters enables the formation of a monomeric elevator architecture. A model of the tripartite PQM complex is experimentally validated and reveals the coupling of the substrate-binding protein to the transporter domains. We use single-molecule total internal reflection fluorescence (TIRF) microscopy in solid-supported lipid bilayers and surface plasmon resonance to study the formation of the tripartite complex and to investigate the impact of interface mutants. Furthermore, we characterize high-affinity single variable domains on heavy chain (VHH) antibodies that bind to the periplasmic side of HiSiaQM and inhibit sialic acid uptake, providing insight into how TRAP transporter function might be inhibited in vivo.

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.

Structural Insights into the Substrate Transport Mechanisms in GTR Transporters through Ensemble Docking

Glucosinolate transporters (GTRs) are part of the nitrate/peptide transporter (NPF) family, members of which also transport specialized secondary metabolites as substrates. Glucosinolates are defense compounds derived from amino acids. We selected 4-methylthiobutyl (4MTB) and indol-3-ylmethyl (I3M) glucosinolates to study how GTR1 from Arabidopsis thaliana transports these substrates in computational simulation approaches. The designed pipeline reported here includes massive docking of 4MTB and I3M in an ensemble of GTR1 conformations (in both inward and outward conformations) extracted from molecular dynamics simulations, followed by clustered and substrate–protein interactions profiling. The identified key residues were mutated, and their role in substrate transport was tested. We were able to identify key residues that integrate a major binding site of these substrates, which is critical for transport activity. In silico approaches employed here represent a breakthrough in the plant transportomics field, as the identification of key residues usually takes a long time if performed from a purely wet-lab experimental perspective. The inclusion of structural bioinformatics in the analyses of plant transporters significantly speeds up the knowledge-gaining process and optimizes valuable time and resources.

Allosteric and transport modulation of human concentrative nucleoside transporter 3 at the atomic scale

Nucleosides are important precursors of nucleotide synthesis in cells, and nucleoside transporters play an important role in many physiological processes by mediating transmembrane transport and absorption. During nucleoside transport, such proteins undergo a significant conformational transition between the outward- and inward-facing states, which leads to alternating access of the substrate-binding site to either side of the membrane. In this work, a variety of molecular simulation methods have been applied to comparatively investigate the motion modes of human concentrative nucleoside transporter 3 (hCNT3) in three states, as well as global and local cavity conformational changes; and finally, a possible elevator-like transport mechanism consistent with experimental data was proposed. The results of the Gaussian network model (GNM) and anisotropic network model (ANM) show that hCNT3 as a whole tends to contract inwards and shift towards a membrane inside, exhibiting an allosteric process that is more energetically favorable than the rigid conversion. To reveal the complete allosteric process of hCNT3 in detail, a series of intermediate conformations were obtained by an adaptive anisotropic network model (aANM). One of the simulated intermediate states is similar to that of a crystal structure, which indicates that the allosteric process is reliable; the state with lower energy is slightly inclined to the inward-facing structure rather than the expected intermediate crystal structure. The final HOLE analysis showed that except for the outward-facing state, the transport channels were gradually enlarged, which was conductive to the directional transport of nucleosides. Our work provides a theoretical basis for the multistep elevator-like transportation mechanism of nucleosides, which helps to further understand the dynamic recognition between nucleoside substrates and hCNT3 as well as the design of nucleoside anticancer drugs.

Thermostability-based binding assays reveal complex interplay of cation, substrate and lipid binding in the bacterial DASS transporter, VcINDY

The divalent anion sodium symporter (DASS) family of transporters (SLC13 family in humans) are key regulators of metabolic homeostasis, disruption of which results in protection from diabetes and obesity, and inhibition of liver cancer cell proliferation. Thus, DASS transporter inhibitors are attractive targets in the treatment of chronic, age-related metabolic diseases. The characterisation of several DASS transporters has revealed variation in the substrate selectivity and flexibility in the coupling ion used to power transport. Here, using the model DASS co-transporter, VcINDY from Vibrio cholerae, we have examined the interplay of the three major interactions that occur during transport: the coupling ion, the substrate, and the lipid environment. Using a series of high-throughput thermostability-based interaction assays, we have shown that substrate binding is Na⁺-dependent; a requirement that is orchestrated through a combination of electrostatic attraction and Na⁺-induced priming of the binding site architecture. We have identified novel DASS ligands and revealed that ligand binding is dominated by the requirement of two carboxylate groups in the ligand that are precisely distanced to satisfy carboxylate interaction regions of the substrate-binding site. We have also identified a complex relationship between substrate and lipid interactions, which suggests a dynamic, regulatory role for lipids in VcINDY's transport cycle.

Cell Physiology and Molecular Mechanism of Anion Transport by Erythrocyte Band 3/AE1

The major transmembrane protein of the red blood cell, known as band 3, AE1, and SLC4A1, has two main functions: 1) catalysis of Cl^-/HCO₃^- exchange, one of the steps in CO₂ excretion; 2) anchoring the membrane skeleton. This review summarizes the 150 year history of research on red cell anion transport and band 3 as an experimental system for studying membrane protein structure and ion transport mechanisms. Important early findings were that red cell Cl^- transport is a tightly coupled 1:1 exchange and band 3 is labeled by stilbenesulfonate derivatives that inhibit anion transport. Biochemical studies showed that the protein is dimeric or tetrameric (paired dimers) and that there is one stilbenedisulfonate binding site per subunit of the dimer. Transport kinetics and inhibitor characteristics supported the idea that the transporter acts by an alternating access mechanism with intrinsic asymmetry. The sequence of band 3 cDNA provided a framework for detailed study of protein topology and amino acid residues important for transport. The identification of genetic variants produced insights into the roles of band 3 in red cell abnormalities and distal renal tubular acidosis. The publication of the membrane domain crystal structure made it possible to propose concrete molecular models of transport. Future research directions include improving our understanding of the transport mechanism at the molecular level and of the integrative relationships among band 3, hemoglobin, carbonic anhydrase, and gradients (both transmembrane and subcellular) of HCO₃^-, Cl^-, O₂, CO₂, pH, and NO metabolites during pulmonary and systemic capillary gas exchange.

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.

Toward unzipping the ZIP metal transporters: structure, evolution, and implications on drug discovery against cancer

The Zrt-/Irt-like protein (ZIP) family consists of divalent metal transporters, ubiquitous in all kingdoms of life. Since the discovery of the first ZIPs in the 1990s, the ZIP family has been expanding to contain tens of thousands of members playing key roles in uptake and homeostasis of life-essential trace elements, primarily zinc, iron and manganese. Some family members are also responsible for toxic metal (particularly cadmium) absorption and distribution. Their central roles in trace element biology, and implications in many human diseases, including cancers, have elicited interest across multiple disciplines for potential applications in biomedicine, agriculture and environmental protection. In this review and perspective, selected areas under rapid progress in the last several years, including structural biology, evolution, and drug discovery against cancers, are summarised and commented. Future research to address the most prominent issues associated with transport and regulation mechanisms are also discussed.

The mannose phosphotransferase system (Man-PTS) - Mannose transporter and receptor for bacteriocins and bacteriophages

Mannose transporters constitute a superfamily (Man-PTS) of the Phosphoenolpyruvate Carbohydrate Phosphotransferase System (PTS). The membrane complexes are homotrimers of protomers consisting of two subunits, IIC and IID. The two subunits without recognizable sequence similarity assume the same fold, and in the protomer are structurally related by a two fold pseudosymmetry axis parallel to membrane-plane (Liu et al. (2019) Cell Research 29 680). Two reentrant loops and two transmembrane helices of each subunit together form the N-terminal transport domain. Two three-helix bundles, one of each subunit, form the scaffold domain. The protomer is stabilized by a helix swap between these bundles. The two C-terminal helices of IIC mediate the interprotomer contacts. PTS occur in bacteria and archaea but not in eukaryotes. Man-PTS are abundant in Gram-positive bacteria living on carbohydrate rich mucosal surfaces. A subgroup of IICIID complexes serve as receptors for class IIa bacteriocins and as channel for the penetration of bacteriophage lambda DNA across the inner membrane. Some Man-PTS are associated with host-pathogen and -symbiont processes.

Structural basis for the reaction cycle of DASS dicarboxylate transporters

Citrate, α-ketoglutarate and succinate are TCA cycle intermediates that also play essential roles in metabolic signaling and cellular regulation. These di- and tricarboxylates are imported into the cell by the divalent anion sodium symporter (DASS) family of plasma membrane transporters, which contains both cotransporters and exchangers. While DASS proteins transport substrates via an elevator mechanism, to date structures are only available for a single DASS cotransporter protein in a substrate-bound, inward-facing state. We report multiple cryo-EM and X-ray structures in four different states, including three hitherto unseen states, along with molecular dynamics simulations, of both a cotransporter and an exchanger. Comparison of these outward- and inward-facing structures reveal how the transport domain translates and rotates within the framework of the scaffold domain through the transport cycle. Additionally, we propose that DASS transporters ensure substrate coupling by a charge-compensation mechanism, and by structural changes upon substrate release.

Functional (un)cooperativity in elevator transport proteins

The activity of enzymes is subject to regulation at multiple levels. Cooperativity, the interconnected behavior of active sites within a protein complex, directly affects protein activity. Cooperativity is a mode of regulation that requires neither extrinsic factors nor protein modifications. Instead, it allows enzymes themselves to modulate reaction rates. Cooperativity is an important regulatory mechanism in soluble proteins, but also examples of cooperative membrane proteins have been described. In this review, we summarize the current knowledge on interprotomer cooperativity in elevator-type proteins, a class of membrane transporters characterized by large rigid-body movements perpendicular to the membrane, and highlight well-studied examples and experimental approaches.

Biosynthesis and Export of Bacterial Glycolipids

Complex carbohydrates are essential for many biological processes, from protein quality control to cell recognition, energy storage, and cell wall formation. Many of these processes are performed in topologically extracellular compartments or on the cell surface; hence, diverse secretion systems evolved to transport the hydrophilic molecules to their sites of action. Polyprenyl lipids serve as ubiquitous anchors and facilitators of these transport processes. Here, we summarize and compare bacterial biosynthesis pathways relying on the recognition and transport of lipid-linked complex carbohydrates. In particular, we compare transporters implicated in O antigen and capsular polysaccharide biosyntheses with those facilitating teichoic acid and N-linked glycan transport. Further, we discuss recent insights into the generation, recognition, and recycling of polyprenyl lipids.

Protein Disulfide Exchange by the Intramembrane Enzymes DsbB, DsbD, and CcdA

The formation of disulfide bonds in proteins is an essential process in both prokaryotes and eukaryotes. In gram-negative bacteria including Escherichia coli, the proteins DsbA and DsbB mediate the formation of disulfide bonds in the periplasm. DsbA acts as the periplasmic oxidant of periplasmic substrate proteins. DsbA is reoxidized by transfer of reducing equivalents to the 4 TM helix membrane protein DsbB, which transfers reducing equivalents to ubiquinone or menaquinone. Multiple structural studies of DsbB have provided detailed structural information on intermediates in the process of DsbB catalyzed oxidation of DsbA. These structures and the insights gained are described. In proteins with more than one pair of Cys residues, there is the potential for formation of non-native disulfide bonds, making it necessary for the cell to have a mechanism for the isomerization of such non-native disulfide bonds. In E. coli, this is mediated by the proteins DsbC and DsbD. DsbC reduces mis-formed disulfide bonds. The eight-TM-helix protein DsbD reduces DsbC and is itself reduced by cytoplasmic thioredoxin. DsbD also contributes reducing equivalents for the reduction of cytochrome c to facilitate heme attachment. The DsbD functional homolog CcdA is a six-TM-helix membrane protein that provides reducing equivalents for the reduction of cytochrome c. A recent structure determination of CcdA has provided critical insights into how reducing equivalents are transferred across the membrane that likely also provides understanding how this is achieved by DsbD as well. This structure and the insights gained are described.

Exploring Conformational Transitions and Free-Energy Profiles of Proton-Coupled Oligopeptide Transporters

Proteins involved in peptide uptake and transport belong to the proton-coupled oligopeptide transporter (POT) family. Crystal structures of POT family members reveal a common fold consisting of two domains of six transmembrane α helices that come together to form a "V" shaped transporter with a central substrate binding site. Proton-coupled oligopeptide transporters operate through an alternate access mechanism, where the membrane transporter undergoes global conformational changes, alternating between inward-facing (IF), outward-facing (OF), and occluded (OC) states. Conformational transitions are promoted by proton and ligand binding; however, due to the absence of crystallographic models of the outward-open state, the role of H⁺ and ligands is still not fully understood. To provide a comprehensive picture of the POT conformational equilibrium, conventional and enhanced sampling molecular dynamics simulations of PepT_st in the presence or absence of ligand and protonation were performed. Free-energy profiles of the conformational variability of PepT_st were obtained from microseconds of adaptive biasing force (ABF) simulations. Our results reveal that both proton and ligand significantly change the conformational free-energy landscape. In the absence of ligand and protonation, only transitions involving IF and OC states are allowed. After protonation of the residue Glu300, the wider free-energy well for Glu300 protonated PepT_st indicates a greater conformational variability relative to the apo system, and OF conformations became accessible. For the Glu300 protonated Holo-PepT_st, the presence of a second free-energy minimum suggests that OF conformations are not only accessible, but also stable. The differences in the free-energy profiles demonstrate that transitions toward outward-facing conformation occur only after protonation, which is likely the first step in the mechanism of peptide transport. Our extensive ABF simulations provide a fully atomic description of all states of the transport process, offering a model for the alternating access mechanism and how protonation and ligand control the conformational changes.

Structural basis for the delivery of activated sialic acid into Golgi for sialyation

The decoration of secretory glycoproteins and glycolipids with sialic acid is critical to many physiological and pathological processes. Sialyation is dependent on a continuous supply of sialic acid into Golgi organelles in the form of CMP-sialic acid. Translocation of CMP-sialic acid into Golgi is carried out by the CMP-sialic acid transporter (CST). Mutations in human CST are linked to glycosylation disorders, and CST is important for glycopathway engineering, as it is critical for sialyation efficiency of therapeutic glycoproteins. The mechanism of how CMP-sialic acid is recognized and translocated across Golgi membranes in exchange for CMP is poorly understood. Here we have determined the crystal structure of a Zea mays CST in complex with CMP. We conclude that the specificity of CST for CMP-sialic acid is established by the recognition of the nucleotide CMP to such an extent that they are mechanistically capable of both passive and coupled antiporter activity.

Multiscale Modeling and Simulation Approaches to Lipid-Protein Interactions

Lipid membranes play a crucial role in living systems by compartmentalizing biological processes and forming a barrier between these processes and the environment. Naturally, a large apparatus of biomolecules is responsible for construction, maintenance, transport, and degradation of these lipid barriers. Additional classes of biomolecules are tasked with transport of specific substances or transduction of signals from the environment across lipid membranes. In this article, we intend to describe a set of techniques that enable one to build accurate models of lipid systems and their associated proteins, and to simulate their dynamics over a variety of time and length scales. We discuss the methods and challenges that allow us to derive structural, mechanistic, and thermodynamic information from these models. We also show how these models have recently been applied in research to study some of the most complex lipid-protein systems to date, including bacterial and viral envelopes, neuronal membranes, and mammalian signaling systems.

Free Energy Landscape of the Complete Transport Cycle in a Key Bacterial Transporter

PepT_So is a proton-coupled bacterial symporter, from the major facilitator superfamily (MFS), which transports di-/tripeptide molecules. The recently obtained crystal structure of PepT_So provides an unprecedented opportunity to gain an understanding of functional insights of the substrate transport mechanism. Binding of the proton and peptide molecule induces conformational changes into occluded (OC) and outward-facing (OF) states, which we are able to characterize using molecular dynamics (MD) simulations. The structural knowledge of the OC and OF state is important to fully understand the major energy barrier associated with the transport cycle. In order to gain functional insight into the interstate dynamics, we performed extensive all atom MD simulations. The Markov state model was constructed to identify the free energy barriers between the states, and kinetic information on intermediate pathways was obtained using the transition pathway theory (TPT). TPT shows that the OF state is obtained by the movement of TM1 and TM7 at the extracellular side approximately 12-16 Å away from each other, and the inward movement of TM4 and TM10 at the intracellular halves to 3-4 Å characterizes the OC state. Helix distance distributions obtained from MD simulations were compared with experimental double electron-electron resonance spectroscopy and were found to be in excellent agreement with previous studies. We also predicted the optimal positions for placement of methane thiosulfonate spin label probes to capture the slowest protein dynamics. Our finding sheds light on the conformational cycle of this key membrane transporter and the functional relationships between the multiple intermediate states.

The molecular mechanism of SLC34 proteins: insights from two decades of transport assays and structure-function studies

The expression cloning some 25 years ago of the first member of SLC34 solute carrier family, the renal sodium-coupled inorganic phosphate cotransporter (NaPi-IIa) from rat and human tissue, heralded a new era of research into renal phosphate handling by focussing on the carrier proteins that mediate phosphate transport. The cloning of NaPi-IIa was followed by that of the intestinal NaPi-IIb and renal NaPi-IIc isoforms. These three proteins constitute the main secondary-active Na⁺-driven pathways for apical entry of inorganic phosphate (P_i) across renal and intestinal epithelial, as well as other epithelial-like organs. The key role these proteins play in mammalian P_i homeostasis was revealed in the intervening decades by numerous in vitro and animal studies, including the development of knockout animals for each gene and the detection of naturally occurring mutations that can lead to P_i-handling dysfunction in humans. In addition to characterising their physiological regulation, research has also focused on understanding the underlying transport mechanism and identifying structure-function relationships. Over the past two decades, this research effort has used real-time electrophysiological and fluorometric assays together with novel computational biology strategies to develop a detailed, but still incomplete, understanding of the transport mechanism of SLC34 proteins at the molecular level. This review will focus on how our present understanding of their molecular mechanism has evolved in this period by highlighting the key experimental findings.

Electrostatic lock in the transport cycle of the multidrug resistance transporter EmrE

EmrE is a small, homodimeric membrane transporter that exploits the established electrochemical proton gradient across the Escherichia coli inner membrane to export toxic polyaromatic cations, prototypical of the wider small-multidrug resistance transporter family. While prior studies have established many fundamental aspects of the specificity and rate of substrate transport in EmrE, low resolution of available structures has hampered identification of the transport coupling mechanism. Here we present a complete, refined atomic structure of EmrE optimized against available cryo-electron microscopy (cryo-EM) data to delineate the critical interactions by which EmrE regulates its conformation during the transport process. With the model, we conduct molecular dynamics simulations of the transporter in explicit membranes to probe EmrE dynamics under different substrate loading and conformational states, representing different intermediates in the transport cycle. The refined model is stable under extended simulation. The water dynamics in simulation indicate that the hydrogen-bonding networks around a pair of solvent-exposed glutamate residues (E14) depend on the loading state of EmrE. One specific hydrogen bond from a tyrosine (Y60) on one monomer to a glutamate (E14) on the opposite monomer is especially critical, as it locks the protein conformation when the glutamate is deprotonated. The hydrogen bond provided by Y60 lowers the [Formula: see text] of one glutamate relative to the other, suggesting both glutamates should be protonated for the hydrogen bond to break and a substrate-free transition to take place. These findings establish the molecular mechanism for the coupling between proton transfer reactions and protein conformation in this proton-coupled secondary transporter.

Conformational transitions of the sodium-dependent sugar transporter, vSGLT

Sodium-dependent transporters couple the flow of Na+ ions down their electrochemical potential gradient to the uphill transport of various ligands. Many of these transporters share a common core structure composed of a five-helix inverted repeat and deliver their cargo utilizing an alternating-access mechanism. A detailed characterization of inward-facing conformations of the Na+-dependent sugar transporter from Vibrio parahaemolyticus (vSGLT) has previously been reported, but structural details on additional conformations and on how Na+ and ligand influence the equilibrium between other states remains unknown. Here, double electron-electron resonance spectroscopy, structural modeling, and molecular dynamics are utilized to deduce ligand-dependent equilibria shifts of vSGLT in micelles. In the absence and presence of saturating amounts of Na+, vSGLT favors an inward-facing conformation. Upon binding both Na+ and sugar, the equilibrium shifts toward either an outward-facing or occluded conformation. While Na+ alone does not stabilize the outward-facing state, gating charge calculations together with a kinetic model of transport suggest that the resting negative membrane potential of the cell, absent in detergent-solubilized samples, may stabilize vSGLT in an outward-open conformation where it is poised for binding external sugars. In total, these findings provide insights into ligand-induced conformational selection and delineate the transport cycle of vSGLT.

Solution structure and elevator mechanism of the membrane electron transporter CcdA

Membrane oxidoredutase CcdA plays a central role in supplying reducing equivalents from the bacterial cytoplasm to the envelope. It transports electrons across the membrane using a single pair of cysteines by a mechanism which has not been elucidated. Here we report an NMR structure of the Thermus thermophilus CcdA (TtCcdA) in an oxidized and outward-facing state. CcdA consists of two inverted structural repeats of three transmembrane helices (2 × 3-TM). We computationally modeled and experimentally validated an inward-facing state, which suggests that CcdA uses an elevator-type movement to shuttle the reactive cysteines across the membrane. CcdA belongs to the LysE superfamily. Its structure may be relevant to other LysE clan transporters. Structure comparisons of CcdA, semiSWEET, Pnu, and major facilitator superfamily (MFS) transporters provide insights about membrane transporter architecture and mechanism.

Asymmetry of inverted-topology repeats in the AE1 anion exchanger suggests an elevator-like mechanism

Anion exchanger 1 catalyzes the transmembrane antiport of chloride and bicarbonate ions through a mechanism that has remained unclear. By modeling its inward-facing state and comparing it with the known outward-facing form, Ficici et al. hypothesize that this transporter features an elevator-like mechanism.

The membrane transporter anion exchanger 1 (AE1), or band 3, is a key component in the processes of carbon-dioxide transport in the blood and urinary acidification in the renal collecting duct. In both erythrocytes and the basolateral membrane of the collecting-duct α-intercalated cells, the role of AE1 is to catalyze a one-for-one exchange of chloride for bicarbonate. After decades of biochemical and functional studies, the structure of the transmembrane region of AE1, which catalyzes the anion-exchange reaction, has finally been determined. Each protomer of the AE1 dimer comprises two repeats with inverted transmembrane topologies, but the structures of these repeats differ. This asymmetry causes the putative substrate-binding site to be exposed only to the extracellular space, consistent with the expectation that anion exchange occurs via an alternating-access mechanism. Here, we hypothesize that the unknown, inward-facing conformation results from inversion of this asymmetry, and we propose a model of this state constructed using repeat-swap homology modeling. By comparing this inward-facing model with the outward-facing experimental structure, we predict that the mechanism of AE1 involves an elevator-like motion of the substrate-binding domain relative to the nearly stationary dimerization domain and to the membrane plane. This hypothesis is in qualitative agreement with a wide range of biochemical and functional data, which we review in detail, and suggests new avenues of experimentation.

Mapping Functionally Important Residues in the Na+/Dicarboxylate Cotransporter, NaDC1

Transporters from the SLC13 family couple the transport of two to four Na⁺ ions with a di- or tricarboxylate, such as succinate or citrate. We have previously modeled mammalian members of the SLC13 family, including the Na⁺/dicarboxylate cotransporter NaDC1 (SLC13A2), based on a structure of the bacterial homologue VcINDY in an inward-facing conformation with one sodium ion bound at the Na1 site. In the study presented here, we modeled the outward-facing conformation of rabbit and human NaDC1 (rbNaDC1 and hNaDC1, respectively) using an outward-facing model of VcINDY as a template and identified residues in or near the putative Na2 and Na3 cation binding sites. Guided by the structural models in both conformations, we performed site-directed mutagenesis in rbNaDC1 for residues proposed to be in the Na⁺ or substrate binding sites. Cysteine substitution of T474 in the predicted Na2 binding site results in an inactive protein. The M539C mutant has a low apparent affinity for both sodium and lithium cations, suggesting that M539 may form part of the putative Na3 binding site. The Y432C and T86C mutants have increased K_m values for succinate, supporting their proposed location in the outward-facing substrate binding site. In addition, cysteine labeling by MTSEA-biotin shows that Y432C is accessible from the outside of the cell, and the accessibility changes in the presence or absence of Na⁺. The results of this study improve our understanding of substrate and ion recognition in the mammalian members of the SLC13 family and provide a framework for developing conformationally specific inhibitors against these transporters.

Rosetta Broker for membrane protein structure prediction: concentrative nucleoside transporter 3 and corticotropin-releasing factor receptor 1 test cases

Background

Membrane proteins are difficult targets for structure prediction due to the limited structural data deposited in Protein Data Bank. Most computational methods for membrane protein structure prediction are based on the comparative modeling. There are only few de novo methods targeting that distinct protein family. In this work an example of such de novo method was used to structurally and functionally characterize two representatives of distinct membrane proteins families of solute carrier transporters and G protein-coupled receptors. The well-known Rosetta program and one of its protocols named Broker was used in two test cases. The first case was de novo structure prediction of three N-terminal transmembrane helices of the human concentrative nucleoside transporter 3 (hCNT3) homotrimer belonging to the solute carrier 28 family of transporters (SLC28). The second case concerned the large scale refinement of transmembrane helices of a homology model of the corticotropin-releasing factor receptor 1 (CRFR1) belonging to the G protein-coupled receptors family.

Results

The inward-facing model of the hCNT3 homotrimer was used to propose the functional impact of its single nucleotide polymorphisms. Additionally, the 100 ns molecular dynamics simulation of the unliganded hCNT3 model confirmed its validity and revealed mobility of the selected binding site and homotrimer interface residues. The large scale refinement of transmembrane helices of the CRFR1 homology model resulted in the significant improvement of its accuracy with respect to the crystal structure of CRFR1, especially in the binding site area. Consequently, the antagonist CP-376395 could be docked with Autodock VINA to the CRFR1 model without any steric clashes.

Conclusions

The presented work demonstrated that Rosetta Broker can be a versatile tool for solving various issues referring to protein biology. Two distinct examples of de novo membrane protein structure prediction presented here provided important insights into three major areas of protein biology. Namely, the dynamics of the inward-facing hCNT3 homotrimer system, the structural changes of the CRFR1 receptor upon the antagonist binding and finally, the role of single nucleotide polymorphisms in both, hCNT3 and CRFR1 proteins, were investigated.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-017-0078-8) contains supplementary material, which is available to authorized users.

Pinning Down the Mechanism of Transport: Probing the Structure and Function of Transporters Using Cysteine Cross-Linking and Site-Specific Labeling

Transporters are crucial in a number of cellular functions, including nutrient uptake, cell signaling, and toxin removal. As such, transporters are important drug targets and their malfunction is related to several disease states. Treating transporter-related diseases and developing pharmaceuticals targeting transporters require an understanding of their mechanism. Achieving a detailed understanding of transporter mechanism depends on an integrative approach involving structural and computational approaches as well as biochemical and biophysical methodologies. Many of the elements of this toolkit exploit the unique and useful chemistry of the amino acid cysteine. Cysteine offers researchers a specific molecular handle with which to precisely modify the protein, which enables the introduction of biophysical probes to assess ligand binding and the conformational ensemble of the transporter, to topologically map transporters and validate structural models, and to assess essential conformational changes. Here, we summarize several uses for cysteine-based labeling and cross-linking in the pursuit of understanding transporter mechanism, the common cysteine-reactive reagents used to probe transporter mechanism, and strategies that can be used to confirm cysteine cross-link formation. In addition, we provide methodological considerations for each approach and a detailed procedure for the cross-linking of introduced cysteines, and a simple screening method to assess cross-link formation.

Cation Interactions and Membrane Potential Induce Conformational Changes in NaPi-IIb

Voltage-dependence of Na(+)-coupled phosphate cotransporters of the SLC34 family arises from displacement of charges intrinsic to the protein and the binding/release of one Na(+) ion in response to changes in the transmembrane electric field. Candidate coordination residues for the cation at the Na1 site were previously predicted by structural modeling using the x-ray structure of dicarboxylate transporter VcINDY as template and confirmed by functional studies. Mutations at Na1 resulted in altered steady-state and presteady-state characteristics that should be mirrored in the conformational changes induced by membrane potential changes. To test this hypothesis by functional analysis, double mutants of the flounder SLC34A2 protein were constructed that contain one of the Na1-site perturbing mutations together with a substituted cysteine for fluorophore labeling, as expressed in Xenopus oocytes. The locations of the mutations were mapped onto a homology model of the flounder protein. The effects of the mutagenesis were characterized by steady-state, presteady-state, and fluorometric assays. Changes in fluorescence intensity (ΔF) in response to membrane potential steps were resolved at three previously identified positions. These fluorescence data corroborated the altered presteady-state kinetics upon perturbation of Na1, and furthermore indicated concomitant changes in the microenvironment of the respective fluorophores, as evidenced by changes in the voltage dependence and time course of ΔF. Moreover, iodide quenching experiments indicated that the aqueous nature of the fluorophore microenvironment depended on the membrane potential. These findings provide compelling evidence that membrane potential and cation interactions induce significant large-scale structural rearrangements of the protein.

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.

Substrate transport and anion permeation proceed through distinct pathways in glutamate transporters

Advances in structure-function analyses and computational biology have enabled a deeper understanding of how excitatory amino acid transporters (EAATs) mediate chloride permeation and substrate transport. However, the mechanism of structural coupling between these functions remains to be established. Using a combination of molecular modeling, substituted cysteine accessibility, electrophysiology and glutamate uptake assays, we identified a chloride-channeling conformer, iChS, transiently accessible as EAAT1 reconfigures from substrate/ion-loaded into a substrate-releasing conformer. Opening of the anion permeation path in this iChS is controlled by the elevator-like movement of the substrate-binding core, along with its wall that simultaneously lines the anion permeation path (global); and repacking of a cluster of hydrophobic residues near the extracellular vestibule (local). Moreover, our results demonstrate that stabilization of iChS by chemical modifications favors anion channeling at the expense of substrate transport, suggesting a mutually exclusive regulation mediated by the movement of the flexible wall lining the two regions.

Visualizing multistep elevator-like transitions of a nucleoside transporter

Membrane transporters move substrates across the membrane by alternating access of their binding sites between the opposite sides of the membrane. An emerging model of this process is the elevator mechanism, in which a substrate-binding transport domain moves a large distance across the membrane. This mechanism has been characterized by a transition between two states, but the conformational path that leads to the transition is not yet known, largely because the available structural information has been limited to the two end states. Here we present crystal structures of the inward-facing, intermediate, and outward-facing states of a concentrative nucleoside transporter from Neisseria wadsworthii. Notably, we determined the structures of multiple intermediate conformations, in which the transport domain is captured halfway through its elevator motion. Our structures present a trajectory of the conformational transition in the elevator model, revealing multiple intermediate steps and state-dependent conformational changes within the transport domain that are associated with the elevator-like motion.

Substituted cysteine accessibility method (SCAM) analysis of the transport domain of human concentrative nucleoside transporter 3 (hCNT3) and other family members reveals features of structural and functional importance

The human SLC28 family of concentrative nucleoside transporter (CNT) proteins has three members: hCNT1, hCNT2, and hCNT3. Na⁺-coupled hCNT1 and hCNT2 transport pyrimidine and purine nucleosides, respectively, whereas hCNT3 transports both pyrimidine and purine nucleosides utilizing Na⁺ and/or H⁺ electrochemical gradients. Escherichia coli CNT family member NupC resembles hCNT1 in permeant selectivity but is H⁺-coupled. Using heterologous expression in Xenopus oocytes and the engineered cysteine-less hCNT3 protein hCNT3(C-), substituted cysteine accessibility method analysis with the membrane-impermeant thiol reactive reagent p-chloromercuribenzene sulfonate was performed on the transport domain (interfacial helix 2, hairpin 1, putative transmembrane domain (TM) 7, and TM8), as well as TM9 of the scaffold domain of the protein. This systematic scan of the entire C-terminal half of hCNT3(C-) together with parallel studies of the transport domain of wild-type hCNT1 and the corresponding TMs of cysteine-less NupC(C-) yielded results that validate the newly developed structural homology model of CNT membrane architecture for human CNTs, revealed extended conformationally mobile regions within transport-domain TMs, identified pore-lining residues of functional importance, and provided evidence of an emerging novel elevator-type mechanism of transporter function.

Microscopic Characterization of Membrane Transporter Function by In Silico Modeling and Simulation

Membrane transporters mediate one of the most fundamental processes in biology. They are the main gatekeepers controlling active traffic of materials in a highly selective and regulated manner between different cellular compartments demarcated by biological membranes. At the heart of the mechanism of membrane transporters lie protein conformational changes of diverse forms and magnitudes, which closely mediate critical aspects of the transport process, most importantly the coordinated motions of remotely located gating elements and their tight coupling to chemical processes such as binding, unbinding and translocation of transported substrate and cotransported ions, ATP binding and hydrolysis, and other molecular events fueling uphill transport of the cargo. An increasing number of functional studies have established the active participation of lipids and other components of biological membranes in the function of transporters and other membrane proteins, often acting as major signaling and regulating elements. Understanding the mechanistic details of these molecular processes require methods that offer high spatial and temporal resolutions. Computational modeling and simulations technologies empowered by advanced sampling and free energy calculations have reached a sufficiently mature state to become an indispensable component of mechanistic studies of membrane transporters in their natural environment of the membrane. In this article, we provide an overview of a number of major computational protocols and techniques commonly used in membrane transporter modeling and simulation studies. The article also includes practical hints on effective use of these methods, critical perspectives on their strengths and weak points, and examples of their successful applications to membrane transporters, selected from the research performed in our own laboratory.

The bacterial dicarboxylate transporter VcINDY uses a two-domain elevator-type mechanism

Secondary transporters use alternating access mechanisms to couple uphill substrate movement to downhill ion flux. Most known transporters utilize a “rocking bundle” motion, where the protein moves around an immobile substrate binding site. However, the glutamate transporter homolog, Glt_Ph, translocates its substrate binding site vertically across the membrane, an “elevator” mechanism. Here, we used the “repeat swap” approach to computationally predict the outward-facing state of the Na⁺/succinate transporter VcINDY, from Vibrio cholerae. Our model predicts a substantial “elevator”-like movement of vcINDY’s substrate binding site, with a vertical translation of ~15 Å and a rotation of ~43°; multiple disulfide crosslinks which completely inhibit transport provide experimental confirmation and demonstrate that such movement is essential. In contrast, crosslinks across the VcINDY dimer interface preserve transport, revealing an absence of large scale coupling between protomers.