Large collective motions regulate the functional properties of glutamate transporter trimers

The Bor1 elevator transport cycle is subject to autoinhibition and activation

Boron, essential for plant growth, necessitates precise regulation due to its potential toxicity. This regulation is achieved by borate transporters (BORs), which are homologous to the SLC4 family. The Arabidopsis thaliana Bor1 (AtBor1) transporter from clade I undergoes slow regulation through degradation and translational suppression, but its potential for fast regulation via direct activity modulation was unclear. Here, we combine cryo-electron microscopy, mutagenesis, and functional characterization to study AtBor1, revealing high-resolution structures of the dimer in one inactive and three active states. Our findings show that AtBor1 is regulated by two distinct mechanisms: an autoinhibitory domain at the carboxyl terminus obstructs the substrate pathway via conserved salt bridges, and phosphorylation of Thr410 allows interaction with a positively charged pocket at the cytosolic face, essential for borate transport. These results elucidate the molecular basis of AtBor1's activity regulation and highlight its role in fast boron level regulation in plants.

Molecular Dynamic Simulations Reveal that Water-Soluble QTY-Variants of Glutamate Transporters EAA1, EAA2 and EAA3 Retain the Conformational Characteristics of Native Transporters

OBJECTIVE

Glutamate transporters play a crucial role in neurotransmitter homeostasis, but studying their structure and function is challenging due to their membrane-bound nature. This study aims to investigate whether water-soluble QTY-variants of glutamate transporters EAA1, EAA2 and EAA3 retain the conformational characteristics and dynamics of native membrane-bound transporters.

METHODS

Molecular dynamics simulations and comparative genomics were used to analyze the structural dynamics of both native transporters and their QTY-variants. Native transporters were simulated in lipid bilayers, while QTY-variants were simulated in aqueous solution. Lipid distortions, relative solvent accessibilities, and conformational changes were examined. Evolutionary conservation profiles were correlated with structural dynamics. Statistical analyses included multivariate analysis to account for confounding variables.

RESULTS

QTY-variants exhibited similar residue-wise conformational dynamics to their native counterparts, with correlation coefficients of 0.73 and 0.56 for EAA1 and EAA3, respectively (p < 0.001). Hydrophobic interactions of native helices correlated with water interactions of QTY- helices (rs = 0.4753, p < 0.001 for EAA1). QTY-variants underwent conformational changes resembling the outward-to-inward transition of native transporters.

CONCLUSIONS

Water-soluble QTY-variants retain key structural properties of native glutamate transporters and mimic aspects of native lipid interactions, including conformational flexibility. This research provides valuable insights into the conformational changes and molecular mechanisms of glutamate transport, potentially offering a new approach for studying membrane protein dynamics and drug interactions.

Synthesis and Structure-Activity Relationships for Glutamate Transporter Allosteric Modulators

Excitatory amino acid transporters (EAATs) are essential CNS proteins that regulate glutamate levels. Excess glutamate release and alteration in EAAT expression are associated with several CNS disorders. Previously, we identified positive allosteric modulators (PAM) of EAAT2, the main CNS transporter, and have demonstrated their neuroprotective properties in vitro. Herein, we report on the structure-activity relationships (SAR) for the analogs identified from virtual screening and from our medicinal chemistry campaign. This work identified several selective EAAT2 positive allosteric modulators (PAMs) such as compounds 4 (DA-023) and 40 (NA-014) from a library of analogs inspired by GT949, an early generation compound. This series also provides nonselective EAAT PAMs, EAAT inhibitors, and inactive compounds that may be useful for elucidating the mechanism of EAAT allosteric modulation.

Discovery of (R)-N-Benzyl-2-(2,5-dioxopyrrolidin-1-yl)propanamide [], a Novel Orally Bioavailable EAAT2 Modulator with Drug-like Properties and Potent Antiseizure Activity In Vivo

(R)-7 [(R)-AS-1] showed broad-spectrum antiseizure activity across in vivo mouse seizure models: maximal electroshock (MES), 6 Hz (32/44 mA), acute pentylenetetrazol (PTZ), and PTZ-kindling. A remarkable separation between antiseizure activity and CNS-related adverse effects was also observed. In vitro studies with primary glia cultures and COS-7 cells expressing the glutamate transporter EAAT2 showed enhancement of glutamate uptake, revealing a stereoselective positive allosteric modulator (PAM) effect, further supported by molecular docking simulations. (R)-7 [(R)-AS-1] was not active in EAAT1 and EAAT3 assays and did not show significant off-target activity, including interactions with targets reported for marketed antiseizure drugs, indicative of a novel and unprecedented mechanism of action. Both in vivo pharmacokinetic and in vitro absorption, distribution, metabolism, excretion, toxicity (ADME-Tox) profiles confirmed the favorable drug-like potential of the compound. Thus, (R)-7 [(R)-AS-1] may be considered as the first-in-class small-molecule PAM of EAAT2 with potential for further preclinical and clinical development in epilepsy and possibly other CNS disorders.

Characterizing unexpected interactions of a glutamine transporter inhibitor with members of the SLC1A transporter family

The solute carrier 1A family comprises a group of membrane proteins that act as dual-function amino acid transporters and chloride (Cl-) channels and includes the alanine serine cysteine transporters (ASCTs) as well as the excitatory amino acid transporters. ASCT2 is regarded as a promising target for cancer therapy, as it can transport glutamine and other neutral amino acids into cells and is upregulated in a range of solid tumors. The compound L-γ-glutamyl-p-nitroanilide (GPNA) is widely used in studies probing the role of ASCT2 in cancer biology; however, the mechanism by which GPNA inhibits ASCT2 is not entirely clear. Here, we used electrophysiology and radiolabelled flux assays to demonstrate that GPNA activates the Cl- conductance of ASCT2 to the same extent as a transported substrate, whilst not undergoing the full transport cycle. This is a previously unreported phenomenon for inhibitors of the solute carrier 1A family but corroborates a body of literature suggesting that the structural requirements for transport are distinct from those for Cl- channel formation. We also show that in addition to its currently known targets, GPNA inhibits several of the excitatory amino acid transporters. Together, these findings raise questions about the true mechanisms of its anticancer effects.

The twisting elevator mechanism of glutamate transporters reveals the structural basis for the dual transport-channel functions

Glutamate transporters facilitate the removal of this excitatory neurotransmitter from the synapse. Increasing evidence indicates that this process is linked to intrinsic chloride channel activity that is thermodynamically uncoupled from substrate transport. A recent cryo-EM structure of Glt_Ph - an archaeal homolog of the glutamate transporters - in an open channel state has shed light on the structural basis for channel opening formed at the interface of two domains within the transporter which is gated by two clusters of hydrophobic residues. These transporters cycle through several conformational states during the transport process, including the chloride conducting state, which appears to be stabilised by protein-membrane interactions and membrane deformation. Several point mutations that perturb the chloride conductance can have detrimental effects and are linked to the pathogenesis of the neurological disorder, episodic ataxia type 6.

Optimized Elastic Network Models With Direct Characterization of Inter-Residue Cooperativity for Protein Dynamics

The elastic network models (ENMs)are known as representative coarse-grained models to capture essential dynamics of proteins. Due to simple designs of the force constants as a decay with spatial distances of residue pairs in many previous studies, there is still much room for the improvement of ENMs. In this article, we directly computed the force constants with the inverse covariance estimation using a ridge-type operater for the precision matrix estimation (ROPE)on a large-scale set of NMR ensembles. Distance-dependent statistical analyses on the force constants were further comprehensively performed in terms of several paired types of sequence and structural information, including secondary structure, relative solvent accessibility, sequence distance and terminal. Various distinguished distributions of the mean force constants highlight the structural and sequential characteristics coupled with the inter-residue cooperativity beyond the spatial distances. We finally integrated these structural and sequential characteristics to build novel ENM variations using the particle swarm optimization for the parameter estimation. The considerable improvements on the correlation coefficient of the mean-square fluctuation and the mode overlap were achieved by the proposed variations when compared with traditional ENMs. This study opens a novel way to develop more accurate elastic network models for protein dynamics.

Sampling of Protein Conformational Space Using Hybrid Simulations: A Critical Assessment of Recent Methods

Recent years have seen several hybrid simulation methods for exploring the conformational space of proteins and their complexes or assemblies. These methods often combine fast analytical approaches with computationally expensive full atomic molecular dynamics (MD) simulations with the goal of rapidly sampling large and cooperative conformational changes at full atomic resolution. We present here a systematic comparison of the utility and limits of four such hybrid methods that have been introduced in recent years: MD with excited normal modes (MDeNM), collective modes-driven MD (CoMD), and elastic network model (ENM)-based generation, clustering, and relaxation of conformations (ClustENM) as well as its updated version integrated with MD simulations (ClustENMD). We analyzed the predicted conformational spaces using each of these four hybrid methods, applied to four well-studied proteins, triosephosphate isomerase (TIM), 3-phosphoglycerate kinase (PGK), HIV-1 protease (PR) and HIV-1 reverse transcriptase (RT), which provide extensive ensembles of experimental structures for benchmarking and comparing the methods. We show that a rigorous multi-faceted comparison and multiple metrics are necessary to properly assess the differences between conformational ensembles and provide an optimal protocol for achieving good agreement with experimental data. While all four hybrid methods perform well in general, being especially useful as computationally efficient methods that retain atomic resolution, the systematic analysis of the same systems by these four hybrid methods highlights the strengths and limitations of the methods and provides guidance for parameters and protocols to be adopted in future studies.

Kinetic mechanism of Na+-coupled aspartate transport catalyzed by GltTk

It is well-established that the secondary active transporters Glt_Tk and Glt_Ph catalyze coupled uptake of aspartate and three sodium ions, but insight in the kinetic mechanism of transport is fragmentary. Here, we systematically measured aspartate uptake rates in proteoliposomes containing purified Glt_Tk, and derived the rate equation for a mechanism in which two sodium ions bind before and another after aspartate. Re-analysis of existing data on Glt_Ph using this equation allowed for determination of the turnover number (0.14 s⁻¹), without the need for error-prone protein quantification. To overcome the complication that purified transporters may adopt right-side-out or inside-out membrane orientations upon reconstitution, thereby confounding the kinetic analysis, we employed a rapid method using synthetic nanobodies to inactivate one population. Oppositely oriented Glt_Tk proteins showed the same transport kinetics, consistent with the use of an identical gating element on both sides of the membrane. Our work underlines the value of bona fide transport experiments to reveal mechanistic features of Na⁺-aspartate symport that cannot be observed in detergent solution. Combined with previous pre-equilibrium binding studies, a full kinetic mechanism of structurally characterized aspartate transporters of the SLC1A family is now emerging.

Trinco et al. measure aspartate uptake rates in proteoliposomes containing purified prokaryotic Na⁺-coupled aspartate transporter GltTk. To overcome limitation of protein orientation, they use synthetic nanobody that blocks transporters from outside and reveal mechanistic features of Na⁺-aspartate symport that cannot be observed in detergent solution.

Elevator-type mechanisms of membrane transport

Membrane transporters are integral membrane proteins that mediate the passage of solutes across lipid bilayers. These proteins undergo conformational transitions between outward- and inward-facing states, which lead to alternating access of the substrate-binding site to the aqueous environment on either side of the membrane. Dozens of different transporter families have evolved, providing a wide variety of structural solutions to achieve alternating access. A sub-set of structurally diverse transporters operate by mechanisms that are collectively named 'elevator-type'. These transporters have one common characteristic: they contain a distinct protein domain that slides across the membrane as a rigid body, and in doing so it 'drags" the transported substrate along. Analysis of the global conformational changes that take place in membrane transporters using elevator-type mechanisms reveals that elevator-type movements can be achieved in more than one way. Molecular dynamics simulations and experimental data help to understand how lipid bilayer properties may affect elevator movements and vice versa.

The high-energy transition state of the glutamate transporter homologue GltPh

Membrane transporters mediate cellular uptake of nutrients, signaling molecules, and drugs. Their overall mechanisms are often well understood, but the structural features setting their rates are mostly unknown. Earlier single-molecule fluorescence imaging of the archaeal model glutamate transporter homologue GltPh from Pyrococcus horikoshii suggested that the slow conformational transition from the outward- to the inward-facing state, when the bound substrate is translocated from the extracellular to the cytoplasmic side of the membrane, is rate limiting to transport. Here, we provide insight into the structure of the high-energy transition state of GltPh that limits the rate of the substrate translocation process. Using bioinformatics, we identified GltPh gain-of-function mutations in the flexible helical hairpin domain HP2 and applied linear free energy relationship analysis to infer that the transition state structurally resembles the inward-facing conformation. Based on these analyses, we propose an approach to search for allosteric modulators for transporters.

Structural ensemble of a glutamate transporter homologue in lipid nanodisc environment

Glutamate transporters are cation-coupled secondary active membrane transporters that clear the neurotransmitter L-glutamate from the synaptic cleft. These transporters are homotrimers, with each protomer functioning independently by an elevator-type mechanism, in which a mobile transport domain alternates between inward- and outward-oriented states. Using single-particle cryo-EM we have determined five structures of the glutamate transporter homologue Glt_Tk, a Na⁺- L-aspartate symporter, embedded in lipid nanodiscs. Dependent on the substrate concentrations used, the protomers of the trimer adopt a variety of asymmetrical conformations, consistent with the independent movement. Six of the 15 resolved protomers are in a hitherto elusive state of the transport cycle in which the inward-facing transporters are loaded with Na⁺ ions. These structures explain how substrate-leakage is prevented - a strict requirement for coupled transport. The belt protein of the lipid nanodiscs bends around the inward oriented protomers, suggesting that membrane deformations occur during transport.

Shared dynamics of LeuT superfamily members and allosteric differentiation by structural irregularities and multimerization

The LeuT-fold superfamily includes secondary active transporters from different functional families, which share a common tertiary structure, despite having a remarkably low sequence similarity. By identifying the common structural and dynamical features upon principal component analysis of a comprehensive ensemble of 90 experimentally resolved structures and anisotropic network model evaluation of collective motions, we provide a unified point of view for understanding the reasons why this particular fold has been selected by evolution to accomplish such a broad spectrum of functions. The parallel identification of conserved sequence features, localized at specific sites of transmembrane helices, sheds light on the role of broken helices (TM1 and TM6 in LeuT) in promoting ion/substrate binding and allosteric interconversion between the outward- and inward-facing conformations of transporters. Finally, the determination of the dynamics landscape for the structural ensemble provides a promising framework for the classification of transporters based on their dynamics, and the characterization of the collective movements that favour multimerization.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

Engineered Cross-Linking to Study the Pore Architecture of the CRAC Channel

ORAI1 constitutes the pore-forming subunit of the calcium release-activated calcium (CRAC) channel, a prototypical store-operated channel that is essential for the activation of cells of the immune system. Here we describe a convenient yet powerful cross-linking approach to examine the pore architecture of CRAC channels using ORAI1 proteins engineered to contain one or two cysteine residues. The generalizable cross-linking in situ approach can also be readily extended to study other integral membrane proteins expressed in various types of cells.

Gaussian network model can be enhanced by combining solvent accessibility in proteins

Gaussian network model (GNM), regarded as the simplest and most representative coarse-grained model, has been widely adopted to analyze and reveal protein dynamics and functions. Designing a variation of the classical GNM, by defining a new Kirchhoff matrix, is the way to improve the residue flexibility modeling. We combined information arising from local relative solvent accessibility (RSA) between two residues into the Kirchhoff matrix of the parameter-free GNM. The undetermined parameters in the new Kirchhoff matrix were estimated by using particle swarm optimization. The usage of RSA was motivated by the fact that our previous work using RSA based linear regression model resulted out higher prediction quality of the residue flexibility when compared with the classical GNM and the parameter free GNM. Computational experiments, conducted based on one training dataset, two independent datasets and one additional small set derived by molecular dynamics simulations, demonstrated that the average correlation coefficients of the proposed RSA based parameter-free GNM, called RpfGNM, were significantly increased when compared with the parameter-free GNM. Our empirical results indicated that a variation of the classical GNMs by combining other protein structural properties is an attractive way to improve the quality of flexibility modeling.

De novo active sites for resurrected Precambrian enzymes

Protein engineering studies often suggest the emergence of completely new enzyme functionalities to be highly improbable. However, enzymes likely catalysed many different reactions already in the last universal common ancestor. Mechanisms for the emergence of completely new active sites must therefore either plausibly exist or at least have existed at the primordial protein stage. Here, we use resurrected Precambrian proteins as scaffolds for protein engineering and demonstrate that a new active site can be generated through a single hydrophobic-to-ionizable amino acid replacement that generates a partially buried group with perturbed physico-chemical properties. We provide experimental and computational evidence that conformational flexibility can assist the emergence and subsequent evolution of new active sites by improving substrate and transition-state binding, through the sampling of many potentially productive conformations. Our results suggest a mechanism for the emergence of primordial enzymes and highlight the potential of ancestral reconstruction as a tool for protein engineering.

Allosteric modulation of human dopamine transporter activity under conditions promoting its dimerization

The human dopamine (DA) transporter (hDAT) is a key regulator of neurotransmission and a target for antidepressants and addictive drugs. Despite the recent resolution of dDAT structures from Drosophila melanogaster, complete understanding of its mechanism of function and even information on its biological assembly is lacking. The resolved dDAT structures are monomeric, but growing evidence suggests that hDAT might function as a multimer, and its oligomerization may be relevant to addictive drug effects. Here, using structure-based computations, we examined the possible mechanisms of hDAT dimerization and its dynamics in a lipid bilayer. Using a combination of site-directed mutagenesis, DA-uptake, and cross-linking experiments that exploited the capacity of Cys-306 to form intermonomeric disulfide bridges in the presence of an oxidizing agent, we tested the effects of mutations at transmembrane segment (TM) 6 and 12 helices in HEK293 cells. The most probable structural model for hDAT dimer suggested by computations and experiments differed from the dimeric structure resolved for the bacterial homolog, LeuT, presumably because of a kink at TM12 preventing favorable monomer packing. Instead, TM2, TM6, and TM11 line the dimer interface. Molecular dynamics simulations of the dimeric hDAT indicated that the two subunits tend to undergo cooperative structural changes, both on local (extracellular gate opening/closure) and global (transition between outward-facing and inward-facing states) scales. These observations suggest that hDAT transport properties may be allosterically modulated under conditions promoting dimerization. Our study provides critical insights into approaches for examining the oligomerization of neurotransmitter transporters and sheds light on their drug modulation.

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.

Dynamical differences of hemoglobin and the ionotropic glutamate receptor in different states revealed by a new dynamics alignment method

A new algorithm for comparison of protein dynamics is presented. Compared protein structures are superposed and their modes of motions are calculated using the anisotropic network model. The obtained modes are aligned using the dynamic programming algorithm of Needleman and Wunsch, commonly used for sequence alignment. Dynamical comparison of hemoglobin in the T and R2 states reveals that the dynamics of the allosteric effector 2,3-bisphosphoglycerate binding site is different in the two states. These differences can contribute to the selectivity of the effector to the T state. Similar comparison of the ionotropic glutamate receptor in the kainate+(R,R)-2b and ZK bound states reveals that the kainate+(R,R)-2b bound states slow modes describe upward motions of ligand binding domain and the transmembrane domain regions. Such motions may lead to the opening of the receptor. The upper lobes of the LBDs of the ZK bound state have a smaller interface with the amino terminal domains above them and have a better ability to move together. The present study exemplifies the use of dynamics comparison as a tool to study protein function. Proteins 2017; 85:1507-1517. © 2014 Wiley Periodicals, Inc.

Identification of Hot Spots in Protein Structures Using Gaussian Network Model and Gaussian Naive Bayes

Residue fluctuations in protein structures have been shown to be highly associated with various protein functions. Gaussian network model (GNM), a simple representative coarse-grained model, was widely adopted to reveal function-related protein dynamics. We directly utilized the high frequency modes generated by GNM and further performed Gaussian Naive Bayes (GNB) to identify hot spot residues. Two coding schemes about the feature vectors were implemented with varying distance cutoffs for GNM and sliding window sizes for GNB based on tenfold cross validations: one by using only a single high mode and the other by combining multiple modes with the highest frequency. Our proposed methods outperformed the previous work that did not directly utilize the high frequency modes generated by GNM, with regard to overall performance evaluated using F1 measure. Moreover, we found that inclusion of more high frequency modes for a GNB classifier can significantly improve the sensitivity. The present study provided additional valuable insights into the relation between the hot spots and the residue fluctuations.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

Revealing an outward-facing open conformational state in a CLC Cl(-)/H(+) exchange transporter

CLC secondary active transporters exchange Cl^- for H⁺. Crystal structures have suggested that the conformational change from occluded to outward-facing states is unusually simple, involving only the rotation of a conserved glutamate (Glu_ex) upon its protonation. Using ¹⁹F NMR, we show that as [H⁺] is increased to protonate Glu_ex and enrich the outward-facing state, a residue ~20 Å away from Glu_ex, near the subunit interface, moves from buried to solvent-exposed. Consistent with functional relevance of this motion, constriction via inter-subunit cross-linking reduces transport. Molecular dynamics simulations indicate that the cross-link dampens extracellular gate-opening motions. In support of this model, mutations that decrease steric contact between Helix N (part of the extracellular gate) and Helix P (at the subunit interface) remove the inhibitory effect of the cross-link. Together, these results demonstrate the formation of a previously uncharacterized 'outward-facing open' state, and highlight the relevance of global structural changes in CLC function.

DOI:http://dx.doi.org/10.7554/eLife.11189.001

Cells have transporter proteins on their surface to carry molecules in and out of the cell. For example, the CLC family of transporters move two chloride ions in one direction at the same time as moving one hydrogen ion in the opposite direction.

To be able to move these ions in opposite directions, transporters have to cycle through a series of shapes in which the ions can only access alternate sides of the membrane. First, the transporter adopts an 'outward-facing' shape when the ions first bind to the transporter, then it switches into the 'occluded' shape to move the ions through the membrane. Finally, the transporter takes on the 'inward-facing' shape to release the ions on the other side of the membrane. However, structural studies of CLCs suggest that the structures of these proteins do not change much while they are moving ions, which suggests that they might work in a different way.

Khantwal, Abraham et al. have now used techniques called “nuclear magnetic resonance” and "double electron-electron resonance" to investigate how a CLC from a bacterium moves ions. The experiments suggest that when the transporter adopts the outward-facing shape, points on the protein known as Y419 and D417 shift their positions. Chemically linking two regions of the CLC prevented this movement and inhibited the transport of chloride ions across the membrane.

Khantwal, Abraham et al. then used a computer simulation to model how the protein changes shape in more detail. This model predicts that two regions of the transporter undergo major rearrangements resulting in a gate-opening motion that widens a passage to allow the chloride ions to bind to the protein. Khantwal, Abraham et al.’s findings will prompt future studies to reveal the other shapes and how CLCs transition between them.

DOI:http://dx.doi.org/10.7554/eLife.11189.002

Computational Studies of Glutamate Transporters

Glutamate is the major excitatory neurotransmitter in the human brain whose binding to receptors on neurons excites them while excess glutamate are removed from synapses via transporter proteins. Determination of the crystal structures of bacterial aspartate transporters has paved the way for computational investigation of their function and dynamics at the molecular level. Here, we review molecular dynamics and free energy calculation methods used in these computational studies and discuss the recent applications to glutamate transporters. The focus of the review is on the insights gained on the transport mechanism through computational methods, which otherwise is not directly accessible by experimental probes. Recent efforts to model the mammalian glutamate and other amino acid transporters, whose crystal structures have not been solved yet, are included in the review.

Blocking dynamics of the SMR transporter EmrE impairs efflux activity

EmrE is a small multidrug resistance transporter that has been well studied as a model for secondary active transport. Because transport requires the protein to convert between at least two states open to opposite sides of the membrane, it is expected that blocking these conformational transitions will prevent transport activity. We have previously shown that NMR can quantitatively measure the transition between the open-in and open-out states of EmrE in bicelles. Now, we have used the antiparallel EmrE crystal structure to design a cross-link to inhibit this conformational exchange process. We probed the structural, dynamic, and functional effects of this cross-link with NMR and in vivo efflux assays. Our NMR results show that our antiparallel cross-link performs as predicted: dramatically reducing conformational exchange while minimally perturbing the overall structure of EmrE and essentially trapping EmrE in a single state. The same cross-link also impairs ethidium efflux activity by EmrE in Escherichia coli. This confirms the hypothesis that transport can be inhibited simply by blocking conformational transitions in a properly folded transporter. The success of our cross-linker design also provides further evidence that the antiparallel crystal structure provides a good model for functional EmrE.

A Mutation in Transmembrane Domain 7 (TM7) of Excitatory Amino Acid Transporters Disrupts the Substrate-dependent Gating of the Intrinsic Anion Conductance and Drives the Channel into a Constitutively Open State

In the mammalian central nervous system, excitatory amino acid transporters (EAATs) are responsible for the clearance of glutamate after synaptic release. This energetically demanding activity is crucial for precise neuronal communication and for maintaining extracellular glutamate concentrations below neurotoxic levels. In addition to their ability to recapture glutamate from the extracellular space, EAATs exhibit a sodium- and glutamate-gated anion conductance. Here we show that substitution of a conserved positively charged residue (Arg-388, hEAAT1) in transmembrane domain 7 with a negatively charged amino acid eliminates the ability of glutamate to further activate the anion conductance. When expressed in oocytes, R388D or R388E mutants show large anion currents that display no further increase in amplitude after application of saturating concentrations of Na(+) and glutamate. They also show a substantially reduced transport activity. The mutant transporters appear to exist preferentially in a sodium- and glutamate-independent constitutive open channel state that rarely transitions to complete the transport cycle. In addition, the accessibility of cytoplasmic residues to membrane-permeant modifying reagents supports the idea that this substrate-independent open state correlates with an intermediate outward facing conformation of the transporter. Our data provide additional insights into the mechanism by which substrates gate the anion conductance in EAATs and suggest that in EAAT1, Arg-388 is a critical element for the structural coupling between the substrate translocation and the gating mechanisms of the EAAT-associated anion channel.

Structure-Encoded Global Motions and Their Role in Mediating Protein-Substrate Interactions

Recent structure-based computational studies suggest that, in contrast to the classical description of equilibrium fluctuations as wigglings and jigglings, proteins have access to well-defined spectra of collective motions, called intrinsic dynamics, encoded by their structure under native state conditions. In particular, the global modes of motions (at the low frequency end of the spectrum) are shown by multiple studies to be highly robust to minor differences in the structure or to detailed interactions at the atomic level. These modes, encoded by the overall fold, usually define the mechanisms of interactions with substrates. They can be estimated by low-resolution models such as the elastic network models (ENMs) exclusively based on interresidue contact topology. The ability of ENMs to efficiently assess the global motions intrinsically favored by the overall fold as well as the relevance of these predictions to the dominant changes in structure experimentally observed for a given protein in the presence of different substrates suggest that the intrinsic dynamics plays a role in mediating protein-substrate interactions. These observations underscore the functional significance of structure-encoded dynamics, or the importance of the predisposition to favor functional global modes in the evolutionary selection of structures.

Differential regulation of two isoforms of the glial glutamate transporter EAAT2 by DLG1 and CaMKII

The gene for EAAT2, the major astrocytic glutamate transporter, generates two carrier isoforms (EAAT2a and EAAT2b) that vary at their C termini as a consequence of alternative RNA splicing. The EAAT2b cytoplasmic C terminus contains a postsynaptic density-95/Discs large/zona occludens-1 (PDZ) ligand, which is absent in EAAT2a. To understand how the distinct C termini might affect transporter trafficking and surface localization, we generated Madin-Darby canine kidney (MDCK) cells that stably express EGFP-EAAT2a or EGFP-EAAT2b and found robust basolateral membrane expression of the EAAT2b isoform. In contrast, EAAT2a displayed a predominant distribution within intracellular vesicle compartments, constitutively cycling to and from the membrane. Addition of the PDZ ligand to EAAT2a as well as its deletion from EAAT2b confirmed the importance of the motif for cell-surface localization. Using EAAT2 constructs with an extracellular biotin acceptor tag to directly assess surface proteins, we observed significant PDZ ligand-dependent EAAT2b surface expression in cultured astrocytes, consistent with observations in cell lines. Discs large homolog 1 (DLG1; SAP97), a PDZ protein prominent in both astrocytes and MDCK cells, colocalized and coimmunoprecipitated with EAAT2b. shRNA knockdown of DLG1 expression decreased surface EAAT2b in both MDCK cells and cultured astrocytes, suggesting that the DLG scaffolding protein stabilizes EAAT2b at the surface. DLG1 can be phosphorylated by Ca(2+)/calmodulin-dependent protein kinase (CaMKII), resulting in disruption of its PDZ-mediated interaction. In murine astrocytes and acute brain slices, activation of CaMKII decreases EAAT2b surface expression but does not alter the distribution of EAAT2a. These data indicate that the surface expression and function of EAAT2b can be rapidly modulated through the disruption of its interaction with DLG1 by CaMKII activation.

Microseconds simulations reveal a new sodium-binding site and the mechanism of sodium-coupled substrate uptake by LeuT

The bacterial sodium-coupled leucine/alanine transporter LeuT is broadly used as a model system for studying the transport mechanism of neurotransmitters because of its structural and functional homology to mammalian transporters such as serotonin, dopamine, or norepinephrine transporters, and because of the resolution of its structure in different states. Although the binding sites (S1 for substrate, and Na1 and Na2 for two co-transported sodium ions) have been resolved, we still lack a mechanistic understanding of coupled Na(+)- and substrate-binding events. We present here results from extensive (>20 μs) unbiased molecular dynamics simulations generated using the latest computing technology. Simulations show that sodium binds initially the Na1 site, but not Na2, and, consistently, sodium unbinding/escape to the extracellular (EC) region first takes place at Na2, succeeded by Na1. Na2 diffusion back to the EC medium requires prior dissociation of substrate from S1. Significantly, Na(+) binding (and unbinding) consistently involves a transient binding to a newly discovered site, Na1″, near S1, as an intermediate state. A robust sequence of substrate uptake events coupled to sodium bindings and translocations between those sites assisted by hydration emerges from the simulations: (i) bindings of a first Na(+) to Na1″, translocation to Na1, a second Na(+) to vacated Na1″ and then to Na2, and substrate to S1; (ii) rotation of Phe(253) aromatic group to seclude the substrate from the EC region; and (iii) concerted tilting of TM1b and TM6a toward TM3 and TM8 to close the EC vestibule.

Global motions exhibited by proteins in micro- to milliseconds simulations concur with anisotropic network model predictions

The Anton supercomputing technology recently developed for efficient molecular dynamics simulations permits us to examine micro- to milli-second events at full atomic resolution for proteins in explicit water and lipid bilayer. It also permits us to investigate to what extent the collective motions predicted by network models (that have found broad use in molecular biophysics) agree with those exhibited by full-atomic long simulations. The present study focuses on Anton trajectories generated for two systems: the bovine pancreatic trypsin inhibitor, and an archaeal aspartate transporter, GltPh. The former, a thoroughly studied system, helps benchmark the method of comparative analysis, and the latter provides new insights into the mechanism of function of glutamate transporters. The principal modes of motion derived from both simulations closely overlap with those predicted for each system by the anisotropic network model (ANM). Notably, the ANM modes define the collective mechanisms, or the pathways on conformational energy landscape, that underlie the passage between the crystal structure and substates visited in simulations. In particular, the lowest frequency ANM modes facilitate the conversion between the most probable substates, lending support to the view that easy access to functional substates is a robust determinant of evolutionarily selected native contact topology.

Exploring the conformational transitions of biomolecular systems using a simple two-state anisotropic network model

Biomolecular conformational transitions are essential to biological functions. Most experimental methods report on the long-lived functional states of biomolecules, but information about the transition pathways between these stable states is generally scarce. Such transitions involve short-lived conformational states that are difficult to detect experimentally. For this reason, computational methods are needed to produce plausible hypothetical transition pathways that can then be probed experimentally. Here we propose a simple and computationally efficient method, called ANMPathway, for constructing a physically reasonable pathway between two endpoints of a conformational transition. We adopt a coarse-grained representation of the protein and construct a two-state potential by combining two elastic network models (ENMs) representative of the experimental structures resolved for the endpoints. The two-state potential has a cusp hypersurface in the configuration space where the energies from both the ENMs are equal. We first search for the minimum energy structure on the cusp hypersurface and then treat it as the transition state. The continuous pathway is subsequently constructed by following the steepest descent energy minimization trajectories starting from the transition state on each side of the cusp hypersurface. Application to several systems of broad biological interest such as adenylate kinase, ATP-driven calcium pump SERCA, leucine transporter and glutamate transporter shows that ANMPathway yields results in good agreement with those from other similar methods and with data obtained from all-atom molecular dynamics simulations, in support of the utility of this simple and efficient approach. Notably the method provides experimentally testable predictions, including the formation of non-native contacts during the transition which we were able to detect in two of the systems we studied. An open-access web server has been created to deliver ANMPathway results.

Many biomolecules are like tiny molecular machines that need to change their shapes and visit many states to perform their biological functions. For a complete molecular understanding of a biological process, one needs to have information on the relevant stable states of the system in question, as well as the pathways by which the system travels from one state to another. We report here an efficient computational method that uses the knowledge of experimental structures of a pair of stable states in order to construct an energetically favoravle pathway between them. We adopt a simple representation of the molecular system by replacing the atoms with beads connected by springs and constructing an energy function with two minima around the end-states. We searched for the structure with highest energy that the system is most likely to visit during the transition and created two paths starting from this structure and proceeding toward the end-states. The combined result of these two paths is the minimum energy pathway between the two stable states. We apply this method to study important structural changes in one enzyme and three large proteins that transport small molecules and ions across the cell membrane.

Molecular dynamics simulations of the mammalian glutamate transporter EAAT3

Excitatory amino acid transporters (EAATs) are membrane proteins that enable sodium-coupled uptake of glutamate and other amino acids into neurons. Crystal structures of the archaeal homolog GltPh have been recently determined both in the inward- and outward-facing conformations. Here we construct homology models for the mammalian glutamate transporter EAAT3 in both conformations and perform molecular dynamics simulations to investigate its similarities and differences from GltPh. In particular, we study the coordination of the different ligands, the gating mechanism and the location of the proton and potassium binding sites in EAAT3. We show that the protonation of the E374 residue is essential for binding of glutamate to EAAT3, otherwise glutamate becomes unstable in the binding site. The gating mechanism in the inward-facing state of EAAT3 is found to be different from that of GltPh, which is traced to the relocation of an arginine residue from the HP1 segment in GltPh to the TM8 segment in EAAT3. Finally, we perform free energy calculations to locate the potassium binding site in EAAT3, and find a high-affinity site that overlaps with the Na1 and Na3 sites in GltPh.

Single molecule analysis reveals coexistence of stable serotonin transporter monomers and oligomers in the live cell plasma membrane

The human serotonin transporter (hSERT) is responsible for the termination of synaptic serotonergic signaling. Although there is solid evidence that SERT forms oligomeric complexes, the exact stoichiometry of the complexes and the fractions of different coexisting oligomeric states still remain enigmatic. Here we used single molecule fluorescence microscopy to obtain the oligomerization state of the SERT via brightness analysis of single diffraction-limited fluorescent spots. Heterologously expressed SERT was labeled either with the fluorescent inhibitor JHC 1-64 or via fusion to monomeric GFP. We found a variety of oligomerization states of membrane-associated transporters, revealing molecular associations larger than dimers and demonstrating the coexistence of different degrees of oligomerization in a single cell; the data are in agreement with a linear aggregation model. Furthermore, oligomerization was found to be independent of SERT surface density, and oligomers remained stable over several minutes in the live cell plasma membrane. Together, the results indicate kinetic trapping of preformed SERT oligomers at the plasma membrane.

Sequence-based Gaussian network model for protein dynamics

Motivation: Gaussian network model (GNM) is widely adopted to analyze and understand protein dynamics, function and conformational changes. The existing GNM-based approaches require atomic coordinates of the corresponding protein and cannot be used when only the sequence is known.

Results: We report, first of its kind, GNM model that allows modeling using the sequence. Our linear regression-based, parameter-free, sequence-derived GNM (L-pfSeqGNM) uses contact maps predicted from the sequence and models local, in the sequence, contact neighborhoods with the linear regression. Empirical benchmarking shows relatively high correlations between the native and the predicted with L-pfSeqGNM B-factors and between the cross-correlations of residue fluctuations derived from the structure- and the sequence-based GNM models. Our results demonstrate that L-pfSeqGNM is an attractive platform to explore protein dynamics. In contrast to the highly used GNMs that require protein structures that number in thousands, our model can be used to study motions for the millions of the readily available sequences, which finds applications in modeling conformational changes, protein–protein interactions and protein functions.

Contact:zerozhua@126.com

Supplementary information:Supplementary data are available at Bioinformatics online.

Mechanisms of glutamate transport

L-Glutamate is the predominant excitatory neurotransmitter in the mammalian central nervous system and plays important roles in a wide variety of brain functions, but it is also a key player in the pathogenesis of many neurological disorders. The control of glutamate concentrations is critical to the normal functioning of the central nervous system, and in this review we discuss how glutamate transporters regulate glutamate concentrations to maintain dynamic signaling mechanisms between neurons. In 2004, the crystal structure of a prokaryotic homolog of the mammalian glutamate transporter family of proteins was crystallized and its structure determined. This has paved the way for a better understanding of the structural basis for glutamate transporter function. In this review we provide a broad perspective of this field of research, but focus primarily on the more recent studies with a particular emphasis on how our understanding of the structure of glutamate transporters has generated new insights.

SLC1 glutamate transporters

The plasma membrane transporters for the neurotransmitter glutamate belong to the solute carrier 1 family. They are secondary active transporters, taking up glutamate into the cell against a substantial concentration gradient. The driving force for concentrative uptake is provided by the cotransport of Na(+) ions and the countertransport of one K(+) in a step independent of the glutamate translocation step. Due to eletrogenicity of transport, the transmembrane potential can also act as a driving force. Glutamate transporters are expressed in many tissues, but are of particular importance in the brain, where they contribute to the termination of excitatory neurotransmission. Glutamate transporters can also run in reverse, resulting in glutamate release from cells. Due to these important physiological functions, glutamate transporter expression and, therefore, the transport rate, are tightly regulated. This review summarizes recent literature on the functional and biophysical properties, structure-function relationships, regulation, physiological significance, and pharmacology of glutamate transporters. Particular emphasis is on the insight from rapid kinetic and electrophysiological studies, transcriptional regulation of transporter expression, and reverse transport and its importance for pathophysiological glutamate release under ischemic conditions.

The SLC1 high-affinity glutamate and neutral amino acid transporter family

Glutamate transporters play important roles in the termination of excitatory neurotransmission and in providing cells throughout the body with glutamate for metabolic purposes. The high-affinity glutamate transporters EAAC1 (SLC1A1), GLT1 (SLC1A2), GLAST (SLC1A3), EAAT4 (SLC1A6), and EAAT5 (SLC1A7) mediate the cellular uptake of glutamate by the co-transport of three sodium ions (Na(+)) and one proton (H(+)), with the counter-transport of one potassium ion (K(+)). Thereby, they protect the CNS from glutamate-induced neurotoxicity. Loss of function of glutamate transporters has been implicated in the pathogenesis of several diseases, including amyotrophic lateral sclerosis and Alzheimer's disease. In addition, glutamate transporters play a role in glutamate excitotoxicity following an ischemic stroke, due to reversed glutamate transport. Besides glutamate transporters, the SLC1 family encompasses two transporters of neutral amino acids, ASCT1 (SLC1A4) and ASCT2 (SLC1A5). Both transporters facilitate electroneutral exchange of amino acids in neurons and/or cells of the peripheral tissues. Some years ago, a high resolution structure of an archaeal homologue of the SLC1 family was determined, followed by the elucidation of its structure in the presence of the substrate aspartate and the inhibitor d,l-threo-benzyloxy aspartate (d,l-TBOA). Historically, the first few known inhibitors of SLC1 transporters were based on constrained glutamate analogs which were active in the high micromolar range but often also showed off-target activity at glutamate receptors. Further development led to the discovery of l-threo-β-hydroxyaspartate derivatives, some of which effectively inhibited SLC1 transporters at nanomolar concentrations. More recently, small molecule inhibitors have been identified whose structures are not based on amino acids. Activators of SLC1 family members have also been discovered but there are only a few examples known.

Large-scale analysis of the dynamics of enzymes

Protein enzymes enable the cell to execute chemical reactions in short time by accelerating the rate of the reactions in a selective manner. The motions or dynamics of the enzymes are essential for their function. Comparison of the dynamics of a set of 1247 nonhomologous enzymes was performed. For each enzyme, the slowest modes of motion are calculated using the Gaussian network model (GNM) and they are globally aligned. Alignment is done using the dynamic programming algorithm of Needleman and Wunsch, commonly used for sequence alignment. Only 96 pairs of proteins were identified to have three similar GNM slow modes with 63 of them having a similar structure. The most frequent slowest mode of motion describes a two domains anticorrelated motion that characterizes at least 23% of the enzymes. Therefore, dynamics uniqueness cannot be accounted for by the slowest mode itself but rather by the combination of several slow modes. Different quaternary structure packing can restrain the motion of enzyme subunits differently and may serve as another mechanism that increases the dynamics uniqueness.

Induced fit substrate binding to an archeal glutamate transporter homologue

Excitatory amino acid transporters (EAATs) are a class of glutamate transporters that terminate glutamatergic synaptic transmission in the mammalian CNS. GltPh, an archeal EAAT homolog from Pyrococcus horikoshii, is currently the only member with a known 3D structure. Here, we studied the kinetics of substrate binding of a single tryptophan mutant (L130W) GltPh in detergent micelles. At low millimolar [Na(+)], the addition of L-aspartate resulted in complex time courses of W130 fluorescence changes over tens of seconds. With increasing [Na(+)], the kinetics were dominated by a fast component [k(obs,fast); K(D) (Na(+)) = 22 ± 3 mM, n(Hill )= 1.7 ± 0.3] with values of k(obs,fast) rising in a saturable manner to ≈ 500 s(-1) (at 6 °C) with increasing [L-aspartate]. The binding kinetics of L-aspartate differed from the binding kinetics of two alternative substrates: L-cysteine sulfinic acid and d-aspartate. L-cysteine sulfinic acid bound with higher affinity than L-aspartate but involved lower saturating rates, whereas the saturating rates after D-aspartate binding were higher. Thus, after the association of two Na(+) to the empty transporter, GltPh binds amino acids by induced fit. Cross-linking and proteolysis experiments suggest that the induced fit results from the closure of helical hairpin 2. This conformational change is faster for GltPh than for most mammalian homologues, whereas the amino acid association rates are similar. Our data reveal the importance of induced fit for substrate selection in EAATs and illustrate how high-affinity binding and the efficient transport of glutamate can be accomplished simultaneously by this class of transporters.

Intracellular gating in an inward-facing state of aspartate transporter Glt(Ph) is regulated by the movements of the helical hairpin HP2

Sodium-coupled neurotransmitter transporters play a key role in neuronal signaling by clearing excess transmitter from the synapse. Structural data on a trimeric archaeal aspartate transporter, Glt_Ph, have provided valuable insights into structural features of human excitatory amino acid transporters. However, the time-resolved mechanisms of substrate binding and release, as well as that of coupling to sodium co-transport, remain largely unknown for this important family. We present here the results of the most extensive simulations performed to date for Glt_Ph in both outward-facing and inward-facing states by taking advantage of significant advances made in recent years in molecular simulation technology. The generated multiple microsecond trajectories consistently show that the helical hairpin HP2, not HP1, serves as an intracellular gate (in addition to its extracellular gating role). In contrast to previous proposals, HP1 can neither initiate nor accommodate neurotransmitter release without prior opening of HP2 by at least 4.0 Å. Aspartate release invariably follows that of a sodium ion located near the HP2 gate entrance. Asp-394 on TM8 and Arg-276 on HP1 emerge as key residues that promote the reorientation and diffusion of substrate toward the cell interior. These findings underscore the significance of examining structural dynamics, as opposed to static structure(s), to make inferences on the mechanisms of transport and key interactions.

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.

Conformational ensemble of the sodium-coupled aspartate transporter

Sodium and aspartate symporter from Pyrococcus horikoshii, Glt_Ph, is a homologue of the mammalian glutamate transporters, homotrimeric integral membrane proteins controlling the neurotransmitter levels in brain synapses. These transporters function by alternating between outward and inward facing states, in which the substrate binding site is oriented toward the extracellular space and the cytoplasm, respectively. Here we employ double electron-electron resonance (DEER) spectroscopy to probe the structure and the state distribution of the subunits in the trimer within distinct hydrophobic environments of detergent micelles and lipid bilayers. Our experiments reveal a conformational ensemble of protomers sampling the outward and inward facing states with nearly equal probabilities, indicative of comparable energies, and independently of each other. On average, the distributions vary only modestly in detergent and in bilayers, but in several mutants unique conformations are stabilized by the latter.

Conformational heterogeneity of the aspartate transporter Glt(Ph)

Glt(Ph) is a Pyrococcus horikoshii homotrimeric Na(+)-coupled aspartate transporter that belongs to the glutamate transporter family. Each protomer consists of a trimerization domain involved in subunit interaction and a transporting domain with the substrate-binding site. Here, we have studied the conformational changes underlying transport by Glt(Ph) using EPR spectroscopy. The trimerization domains form a rigid scaffold, whereas the transporting domains sample multiple conformations, consistent with large-scale movements during the transport cycle. Binding of substrates changed the occupancies of the different conformational states, but the domains remained heterogeneous. The membrane environment favored conformations different from those observed in detergent micelles, but the transporting domain remained structurally heterogeneous in both environments. We conclude that the transporting domains sample multiple conformational states with substantial occupancy regardless of the presence of substrate and coupling ions, consistent with equilibrium constants close to unity between the observed transporter conformations.

Molecular machines

Molecular machines (MM) are essential components of living cells. They conduct mechanical work, transport materials into and out of cells, assist in processing enzymatic reactions, and more. Their operations are frequently combined with significant conformational transitions. Computational studies of these conformational transitions and their coupling to molecular functions are discussed. It is argued that coarse descriptions of these molecules which are based on mass density and shape provide useful information on directions of action. It is further argued that MM are likely to have well focused and narrow reaction pathways. The proposal for such pathways is supported by evolutionary analyses of homologous machines. Finally, these observations are used to build atomically detailed models of these systems that are making the link from structure to functions (kinetics and thermodynamics). For that purpose enhanced sampling techniques are required.

Cysteine cross-linking defines the extracellular gate for the Leishmania donovani nucleoside transporter 1.1 (LdNT1.1)

Equilibrative nucleoside transporters are a unique family of proteins that enable uptake of nucleosides/nucleobases into a wide range of eukaryotes and internalize a myriad of drugs used in the treatment of cancer, heart disease, AIDs, and parasitic infections. In previous work we generated a structural model for such a transporter, the LdNT1.1 nucleoside permease from the parasitic protozoan Leishmania donovani, using ab initio computation. The model suggested that aromatic residues present in transmembrane helices 1, 2, and 7 interact to form an extracellular gate that closes the permeation pathway in the inward-open conformation. Mutation of residues Phe-48(TM1) and Trp-75(TM2) abrogated transport activity, consistent with such prediction. In this study cysteine mutagenesis and oxidative cross-linking were combined to analyze proximity relationships of helices 1, 2, and 7 in LdNT1.1. Disulfide bond formation between introduced paired cysteines at the interface of such helices (A61C(TM1)/F74C(TM2), A61C(TM1)/G350C(TM7), and F74C(TM2)/G350C(TM7)) was analyzed by transport measurement and gel mobility shifts upon oxidation with Cu (II)-(1,10-phenanthroline)(3). In all cases cross-linking inhibited transport. However, if LdNT1.1 ligands were included during cross-linking, inhibition of transport was reduced, suggesting that ligands moved the three gating helices apart. Moreover, all paired cysteine mutants exhibited a mobility shift upon oxidation, corroborating the formation of a disulfide bond. These data support the notion that helices 1, 2, and 7 constitute the extracellular gate of LdNT1.1, thus further validating the computational model and the previously demonstrated importance of F48(TM1) and Trp-75(TM2) in tethering together helices that are part of the gate.

Biomolecular dynamics: order-disorder transitions and energy landscapes

While the energy landscape theory of protein folding is now a widely accepted view for understanding how relatively weak molecular interactions lead to rapid and cooperative protein folding, such a framework must be extended to describe the large-scale functional motions observed in molecular machines. In this review, we discuss (1) the development of the energy landscape theory of biomolecular folding, (2) recent advances toward establishing a consistent understanding of folding and function and (3) emerging themes in the functional motions of enzymes, biomolecular motors and other biomolecular machines. Recent theoretical, computational and experimental lines of investigation have provided a very dynamic picture of biomolecular motion. In contrast to earlier ideas, where molecular machines were thought to function similarly to macroscopic machines, with rigid components that move along a few degrees of freedom in a deterministic fashion, biomolecular complexes are only marginally stable. Since the stabilizing contribution of each atomic interaction is on the order of the thermal fluctuations in solution, the rigid body description of molecular function must be revisited. An emerging theme is that functional motions encompass order-disorder transitions and structural flexibility provides significant contributions to the free energy. In this review, we describe the biological importance of order-disorder transitions and discuss the statistical-mechanical foundation of theoretical approaches that can characterize such transitions.

Charge compensation mechanism of a Na+-coupled, secondary active glutamate transporter

Forward glutamate transport by the excitatory amino acid carrier EAAC1 is coupled to the inward movement of three Na(+) and one proton and the subsequent outward movement of one K(+) in a separate step. Based on indirect evidence, it was speculated that the cation binding sites bear a negative charge. However, little is known about the electrostatics of the transport process. Valences calculated using the Poisson-Boltzmann equation indicate that negative charge is transferred across the membrane when only one cation is bound. Consistently, transient currents were observed in response to voltage jumps when K(+) was the only cation on both sides of the membrane. Furthermore, rapid extracellular K(+) application to EAAC1 under single turnover conditions (K(+) inside) resulted in outward transient current. We propose a charge compensation mechanism, in which the C-terminal transport domain bears an overall negative charge of -1.23. Charge compensation, together with distribution of charge movement over many steps in the transport cycle, as well as defocusing of the membrane electric field, may be combined strategies used by Na(+)-coupled transporters to avoid prohibitive activation barriers for charge translocation.