Atomistic models of ion and solute transport by the sodium-dependent secondary active transporters

Procyanidin C1 from Viola odorata L. inhibits Na+,K+-ATPase

Members of the Viola genus play important roles in traditional Asian herbal medicine. This study investigates the ability of Viola odorata L. extracts to inhibit Na⁺,K⁺-ATPase, an essential animal enzyme responsible for membrane potential maintenance. The root extract of V. odorata strongly inhibited Na⁺,K⁺-ATPase, while leaf and seeds extracts were basically inactive. A UHPLC-QTOF-MS/MS metabolomic approach was used to identify the chemical principle of the root extract’s activity, resulting in the detection of 35,292 features. Candidate active compounds were selected by correlating feature area with inhibitory activity in 14 isolated fractions. This yielded a set of 15 candidate compounds, of which 14 were preliminarily identified as procyanidins. Commercially available procyanidins (B1, B2, B3 and C1) were therefore purchased and their ability to inhibit Na⁺,K⁺-ATPase was investigated. Dimeric procyanidins B1, B2 and B3 were found to be inactive, but the trimeric procyanidin C1 strongly inhibited Na⁺,K⁺-ATPase with an IC₅₀ of 4.5 µM. This newly discovered inhibitor was docked into crystal structures mimicking the Na₃E₁∼P·ADP and K₂E₂·P_i states to identify potential interaction sites within Na⁺,K⁺-ATPase. Possible binding mechanisms and the principle responsible for the observed root extract activity are discussed.

Hydrogen-Deuterium Exchange Mass-Spectrometry of Secondary Active Transporters: From Structural Dynamics to Molecular Mechanisms

Membrane transporters allow the selective transport of otherwise poorly permeable solutes across the cell membrane and thus, play a key role in maintaining cellular homeostasis in all kingdoms of life. Importantly, these proteins also serve as important drug targets. Over the last decades, major progress in structural biology methods has elucidated important structure-function relationships in membrane transporters. However, structures obtained using methods such as X-ray crystallography and high-resolution cryogenic electron microscopy merely provide static snapshots of an intrinsically dynamic, multi-step transport process. Therefore, there is a growing need for developing new experimental approaches capable of exploiting the data obtained from the high-resolution snapshots in order to investigate the dynamic features of membrane proteins. Here, we present the basic principles of hydrogen-deuterium exchange mass-spectrometry (HDX-MS) and recent advancements in its use to study membrane transporters. In HDX-MS experiments, minute amounts of a protein sample can be used to investigate its structural dynamics under native conditions, without the need for chemical labelling and with practically no limit on the protein size. Thus, HDX-MS is instrumental for resolving the structure-dynamic landscapes of membrane proteins in their apo (ligand-free) and ligand-bound forms, shedding light on the molecular mechanism underlying the transport process and drug binding.

Quantitative Assessment of the Energetics of Dopamine Translocation by Human Dopamine Transporter

Computational evaluation of the energetics of substrate binding, transport, and release events of neurotransmitter transporters at the molecular level is a challenge, as the structural transitions of these membrane proteins involve coupled global and local changes that span time scales of several orders of magnitude, from nanoseconds to seconds. Here, we provide a quantitative assessment of the energetics of dopamine (DA) translocation through the human DA transporter (hDAT), using a combination of molecular modeling, simulation, and analysis tools. DA-binding and -unbinding events, which generally involve local configurational changes, are evaluated using free-energy perturbation or adaptive biasing force methods. The global transitions between the outward-facing state and the inward-facing state, on the other hand, require a dual-boost accelerated molecular dynamics simulation. We present results on DA-binding/unbinding energetics under different conditions, as well as the conformational energy landscape of hDAT in both DA-bound and -unbound states. The study provides a tractable method of approach for quantitative evaluation of substrate-binding energetics and efficient estimation of conformational energy landscape, in general.

Interactions among triphenyltin degradation, phospholipid synthesis and membrane characteristics of Bacillus thuringiensis in the presence of d-malic acid

Degradation pathway and surface biosorption of triphenyltin (TPT) by effective microbes have been investigated in the past. However, unclear interactions among membrane components and TPT binding and transport are still obstacles to understanding TPT biotransformation. To reveal the mechanism involved, the phospholipid expression, membrane potential, cellular mechanism and molecular dynamics between TPT and fatty acids (FAs) during the TPT degradation process in the presence of d-malic acid (DMA) were studied. The results show that the degradation efficiency of 1 mg L^-1 TPT by Bacillus thuringiensis (1 g L^-1) with 0.5 or 1 mg L^-1 DMA reached values up to approximately 90% due to the promotion of element metabolism and cellular activity, and the depression of FA synthesis induced by DMA. The addition of DMA caused conversion of more linoleic acid into 10-oxo-12(Z)-octadecenoic acid, increased the membrane permeability, and alleviated the decrease in membrane potential, resulting in TPT transport and degradation. Fluorescence analysis reveals that the endospore of B. thuringiensis could act as an indicator for membrane potential and cellular activities. The current findings are advantageous for acceleration of biosorption, transport and removal of pollutants from natural environments.

Glycine Transporters in Glia Cells: Structural Studies

Glycine, besides exerting essential metabolic functions, is an important inhibitory neurotransmitter in caudal areas of the central nervous system and also a positive neuromodulator at excitatory glutamate-mediated synapses. Glial cells provide metabolic support to neurons and modulate synaptic activity. Six transporters belonging to three solute carrier families (SLC6, SLC38, and SLC7) are capable of transporting glycine across the glial plasma membrane. The unique glial glycine-selective transporter GlyT1 (SLC6) is the main regulator of synaptic glycine concentrations, assisted by the neuronal GlyT2. The five additional glycine transporters ATB^0,+, SNAT1, SNAT2, SNAT5, and LAT2 display broad amino acid specificity and have differential contributions to glial glycine transport. Glial glycine transporters are divergent in sequence but share a similar architecture displaying the 5 + 5 inverted fold originally characterized in the leucine transporter LeuT. The availability of protein crystals solved at high resolution for prokaryotic and, more recently, eukaryotic homologues of this superfamily has advanced significantly our understanding of the mechanism of glycine transport.

Properties of an inward-facing state of LeuT: conformational stability and substrate release

The leucine transporter (LeuT) is a bacterial homolog of the human monoamine transporters, which are important pharmaceutical targets. There are no high-resolution structures of the human transporters available; however, LeuT has been crystallized in several different conformational states. Recently, an inward-facing conformation of LeuT was solved revealing an unexpectedly large movement of transmembrane helix 1a (TM1a). We have performed molecular dynamics simulations of the mutated and wild-type transporter, with and without the cocrystallized Fab antibody fragment, to investigate the properties of this inward-facing conformation in relation to transport by LeuT within the membrane environment. In all of the simulations, local conformational changes with respect to the crystal structure are consistently observed, especially in TM1a. Umbrella sampling revealed a soft potential for TM1a tilting. Furthermore, simulations of inward-facing LeuT with Na(+) ions and substrate bound suggest that one of the Na(+) ion binding sites is fully disrupted. Release of alanine and the second Na(+) ion is also observed, giving insight into the final stage of the translocation process in atomistic detail.

Identification of the first sodium binding site of the phosphate cotransporter NaPi-IIa (SLC34A1)

Transporters of the SLC34 family (NaPi-IIa,b,c) catalyze uptake of inorganic phosphate (Pi) in renal and intestinal epithelia. The transport cycle requires three Na(+) ions and one divalent Pi to bind before a conformational change enables translocation, intracellular release of the substrates, and reorientation of the empty carrier. The electrogenic interaction of the first Na(+) ion with NaPi-IIa/b at a postulated Na1 site is accompanied by charge displacement, and Na1 occupancy subsequently facilitates binding of a second Na(+) ion at Na2. The voltage dependence of cotransport and presteady-state charge displacements (in the absence of a complete transport cycle) are directly related to the molecular architecture of the Na1 site. The fact that Li(+) ions substitute for Na(+) at Na1, but not at the other sites (Na2 and Na3), provides an additional tool for investigating Na1 site-specific events. We recently proposed a three-dimensional model of human SLC34a1 (NaPi-IIa) including the binding sites Na2, Na3, and Pi based on the crystal structure of the dicarboxylate transporter VcINDY. Here, we propose nine residues in transmembrane helices (TM2, TM3, and TM5) that potentially contribute to Na1. To verify their roles experimentally, we made single alanine substitutions in the human NaPi-IIa isoform and investigated the kinetic properties of the mutants by voltage clamp and (32)P uptake. Substitutions at five positions in TM2 and one in TM5 resulted in relatively small changes in the substrate apparent affinities, yet at several of these positions, we observed significant hyperpolarizing shifts in the voltage dependence. Importantly, the ability of Li(+) ions to substitute for Na(+) ions was increased compared with the wild-type. Based on these findings, we adjusted the regions containing Na1 and Na3, resulting in a refined NaPi-IIa model in which five positions (T200, Q206, D209, N227, and S447) contribute directly to cation coordination at Na1.

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.

Conformational Changes in Two Inter-Helical Loops of Mhp1 Membrane Transporter

Mhp1 is a bacterial secondary transporter with high-resolution crystal structures available for both the outward- and inward-facing conformations. Through molecular dynamics simulations of the ligand-free Mhp1 as well as analysis of its crystal structures, here we show that two inter-helical loops, respectively located at the extra- and intracellular ends of the "hash motif" in the protein, play important roles in the conformational transition. In the outward- and inward-facing states of the protein, the loops adopt different secondary structures, either wrapped to the end of an alpha-helix, or unwrapped to extended conformations. In equilibrium simulations of 100 ns with Mhp1 in explicit lipids and water, the loop conformations remain largely stable. In targeted molecular dynamics simulations with the protein structure driven from one state to the other, the loops exhibit resistance and only undergo abrupt changes when other parts of the protein already approach the target conformation. Free energy calculations on the isolated loops further confirm that the wrapping/unwrapping transitions are associated with substantial energetic barriers, and consist of multiple sequential steps involving the rotation of certain backbone torsion angles. Furthermore, in simulations with the loops driven from one state to the other, a large part of the protein follows the loops to the target conformation. Taken together, our simulations suggest that changes of the loop secondary structures would be among the slow degrees of freedom in the conformational transition of the entire protein. Incorporation of detailed loop structures into the reaction coordinate, therefore, should improve the convergence and relevance of the resulting conformational free energy.

All-atom molecular dynamics simulations of an artificial sodium channel in a lipid bilayer: the effect of water solvation/desolvation of the sodium ion

All-atom molecular dynamics is used to investigate the transport of Na(+) across a 1,2-dioleoyl-sn-glycero-3-phosphocholine lipid bilayer facilitated by a diazacrown hydraphile. Specifically, the free energy of Na(+) passing through the bilayer is calculated using the adaptive biasing force method to study the free energy associated with the increase in Na(+) transport in the presence of the hydraphile molecule. The results show that water interaction greatly influences Na(+) transport through the lipid bilayer as water is pulled through the bilayer with Na(+) forming a water channel. The hydraphile causes a reduction in the free energy barrier for the transport of Na(+) through the head group part of the lipid bilayer since it complexes the Na(+) reducing the necessity for water to be complexed and, therefore, dragged through with Na(+), an energetically unfavorable process. The free energy associated with Na(+) being desolvated within the bilayer is significantly decreased in the presence of the hydraphile molecule; the hydraphile increases the number of solvation states of Na(+) that can be adopted, and this increase in the number of available configurations provides an entropic explanation for the success of the hydraphile.

Sodium recognition by the Na+/Ca2+ exchanger in the outward-facing conformation

Na(+)/Ca(2+) exchangers (NCXs) are ubiquitous membrane transporters with a key role in Ca(2+) homeostasis and signaling. NCXs mediate the bidirectional translocation of either Na(+) or Ca(2+), and thus can catalyze uphill Ca(2+) transport driven by a Na(+) gradient, or vice versa. In a major breakthrough, a prokaryotic NCX homolog (NCX_Mj) was recently isolated and its crystal structure determined at atomic resolution. The structure revealed an intriguing architecture consisting of two inverted-topology repeats, each comprising five transmembrane helices. These repeats adopt asymmetric conformations, yielding an outward-facing occluded state. The crystal structure also revealed four putative ion-binding sites, but the occupancy and specificity thereof could not be conclusively established. Here, we use molecular-dynamics simulations and free-energy calculations to identify the ion configuration that best corresponds to the crystallographic data and that is also thermodynamically optimal. In this most probable configuration, three Na(+) ions occupy the so-called Sext, SCa, and Sint sites, whereas the Smid site is occupied by one water molecule and one H(+), which protonates an adjacent aspartate side chain (D240). Experimental measurements of Na(+)/Ca(2+) and Ca(2+)/Ca(2+) exchange by wild-type and mutagenized NCX_Mj confirm that transport of both Na(+) and Ca(2+) requires protonation of D240, and that this side chain does not coordinate either ion at Smid. These results imply that the ion exchange stoichiometry of NCX_Mj is 3:1 and that translocation of Na(+) across the membrane is electrogenic, whereas transport of Ca(2+) is not. Altogether, these findings provide the basis for further experimental and computational studies of the conformational mechanism of this exchanger.

Computational and biochemical docking of the irreversible cocaine analog RTI 82 directly demonstrates ligand positioning in the dopamine transporter central substrate-binding site

The dopamine transporter (DAT) functions as a key regulator of dopaminergic neurotransmission via re-uptake of synaptic dopamine (DA). Cocaine binding to DAT blocks this activity and elevates extracellular DA, leading to psychomotor stimulation and addiction, but the mechanisms by which cocaine interacts with DAT and inhibits transport remain incompletely understood. Here, we addressed these questions using computational and biochemical methodologies to localize the binding and adduction sites of the photoactivatable irreversible cocaine analog 3β-(p-chlorophenyl)tropane-2β-carboxylic acid, 4'-azido-3'-iodophenylethyl ester ([(125)I]RTI 82). Comparative modeling and small molecule docking indicated that the tropane pharmacophore of RTI 82 was positioned in the central DA active site with an orientation that juxtaposed the aryliodoazide group for cross-linking to rat DAT Phe-319. This prediction was verified by focused methionine substitution of residues flanking this site followed by cyanogen bromide mapping of the [(125)I]RTI 82-labeled mutants and by the substituted cysteine accessibility method protection analyses. These findings provide positive functional evidence linking tropane pharmacophore interaction with the core substrate-binding site and support a competitive mechanism for transport inhibition. This synergistic application of computational and biochemical methodologies overcomes many uncertainties inherent in other approaches and furnishes a schematic framework for elucidating the ligand-protein interactions of other classes of DA transport inhibitors.

The two Na+ sites in the human serotonin transporter play distinct roles in the ion coupling and electrogenicity of transport

Neurotransmitter transporters of the SLC6 family of proteins, including the human serotonin transporter (hSERT), utilize Na(+), Cl(-), and K(+) gradients to induce conformational changes necessary for substrate translocation. Dysregulation of ion movement through monoamine transporters has been shown to impact neuronal firing potentials and could play a role in pathophysiologies, such as depression and anxiety. Despite multiple crystal structures of prokaryotic and eukaryotic SLC transporters indicating the location of both (or one) conserved Na(+)-binding sites (termed Na1 and Na2), much remains uncertain in regard to the movements and contributions of these cation-binding sites in the transport process. In this study, we utilize the unique properties of a mutation of hSERT at a single, highly conserved asparagine on TM1 (Asn-101) to provide several lines of evidence demonstrating mechanistically distinct roles for Na1 and Na2. Mutations at Asn-101 alter the cation dependence of the transporter, allowing Ca(2+) (but not other cations) to functionally replace Na(+) for driving transport and promoting 5-hydroxytryptamine (5-HT)-dependent conformational changes. Furthermore, in two-electrode voltage clamp studies in Xenopus oocytes, both Ca(2+) and Na(+) illicit 5-HT-induced currents in the Asn-101 mutants and reveal that, although Ca(2+) promotes substrate-induced current, it does not appear to be the charge carrier during 5-HT transport. These findings, in addition to functional evaluation of Na1 and Na2 site mutants, reveal separate roles for Na1 and Na2 and provide insight into initiation of the translocation process as well as a mechanism whereby the reported SERT stoichiometry can be obtained despite the presence of two putative Na(+)-binding sites.

Amino acid secondary transporters: toward a common transport mechanism

Solute carriers (SLC) that transport amino acids are key players in health and diseases in humans. Their prokaryotic relatives are often involved in essential physiological processes in microorganisms, e.g. in homeostasis and acidic/osmotic stress response. High-resolution X-ray structures of the sequence-unrelated amino acid transporters unraveled a striking structural similarity between carriers, which were formerly assigned to different families. The highly conserved fold is characterized by two inverted structural repeats of five transmembrane helices each and indicates common mechanistic transport concepts if not an evolutionary link among a large number of amino acid transporters. Therefore, these transporters are classified now into the structural amino acid-polyamine-organocation superfamily (APCS). The APCS includes among others the mammalian SLC6 transporters and the heterodimeric SLC7/SLC3 transporters. However, it has to be noted that the APCS is not limited entirely to amino acid transporters but contains also transporters for, e.g. amino acid derivatives and sugars. For instance, the betaine-choline-carnitine transporter family of bacterial activity-regulated Na(+)- and H(+)-coupled symporters for glycine betaine and choline is also part of this second largest structural superfamily. The APCS fold provides different possibilities to transport the same amino acid. Arginine can be transported by an H(+)-coupled symport or by antiport mechanism in exchange against agmatine for example. The convergence of the mechanistic concept of transport under comparable physiological conditions allows speculating if structurally unexplored amino acid transporters, e.g. the members of the SLC36 and SLC38 family, belong to the APCS, too. In the kidney, which is an organ that depends critically on the regulated amino acid transport, these different SLC transporters have to work together to account for proper function. Here, we will summarize the basic concepts of Na(+)- and H(+)-coupled amino acid symport and amino acid-product antiport in the light of the respective physiological requirements.

Visualizing functional motions of membrane transporters with molecular dynamics simulations

Computational modeling and molecular simulation techniques have become an integral part of modern molecular research. Various areas of molecular sciences continue to benefit from, indeed rely on, the unparalleled spatial and temporal resolutions offered by these technologies, to provide a more complete picture of the molecular problems at hand. Because of the continuous development of more efficient algorithms harvesting ever-expanding computational resources, and the emergence of more advanced and novel theories and methodologies, the scope of computational studies has expanded significantly over the past decade, now including much larger molecular systems and far more complex molecular phenomena. Among the various computer modeling techniques, the application of molecular dynamics (MD) simulation and related techniques has particularly drawn attention in biomolecular research, because of the ability of the method to describe the dynamical nature of the molecular systems and thereby to provide a more realistic representation, which is often needed for understanding fundamental molecular properties. The method has proven to be remarkably successful in capturing molecular events and structural transitions highly relevant to the function and/or physicochemical properties of biomolecular systems. Herein, after a brief introduction to the method of MD, we use a number of membrane transport proteins studied in our laboratory as examples to showcase the scope and applicability of the method and its power in characterizing molecular motions of various magnitudes and time scales that are involved in the function of this important class of membrane proteins.