Mechanism and Energetics of Ligand Release in the Aspartate Transporter GltPh

Observing spontaneous, accelerated substrate binding in molecular dynamics simulations of glutamate transporters

Glutamate transporters are essential for removing the neurotransmitter glutamate from the synaptic cleft. Glutamate transport across the membrane is associated with elevator-like structural changes of the transport domain. These structural changes require initial binding of the organic substrate to the transporter. Studying the binding pathway of ligands to their protein binding sites using molecular dynamics (MD) simulations requires micro-second level simulation times. Here, we used three methods to accelerate aspartate binding to the glutamate transporter homologue Gltph and to investigate the binding pathway. 1) Two methods using user-defined forces to prevent the substrate from diffusing too far from the binding site. 2) Conventional MD simulations using very high substrate concentrations in the 0.1 M range. The final, substrate bound states from these methods are comparable to the binding pose observed in crystallographic studies, although they show more flexibility in the side chain carboxylate function. We also captured an intermediate on the binding pathway, where conserved residues D390 and D394 stabilize the aspartate molecule. Finally, we investigated glutamate binding to the mammalian glutamate transporter, excitatory amino acid transporter 1 (EAAT1), for which a crystal structure is known, but not in the glutamate-bound state. Overall, the results obtained in this study reveal new insights into the pathway of substrate binding to glutamate transporters, highlighting intermediates on the binding pathway and flexible conformational states of the side chain, which most likely become locked in once the hairpin loop 2 closes to occlude the substrate.

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.

Millisecond dynamics of an unlabeled amino acid transporter

Excitatory amino acid transporters (EAATs) are important in many physiological processes and crucial for the removal of excitatory amino acids from the synaptic cleft. Here, we develop and apply high-speed atomic force microscopy line-scanning (HS-AFM-LS) combined with automated state assignment and transition analysis for the determination of transport dynamics of unlabeled membrane-reconstituted Glt_Ph, a prokaryotic EAAT homologue, with millisecond temporal resolution. We find that Glt_Ph transporters can operate much faster than previously reported, with state dwell-times in the 50 ms range, and report the kinetics of an intermediate transport state with height between the outward- and inward-facing states. Transport domains stochastically probe transmembrane motion, and reversible unsuccessful excursions to the intermediate state occur. The presented approach and analysis methodology are generally applicable to study transporter kinetics at system-relevant temporal resolution.

Excitatory amino acid transporters (EAATs) are crucial for the removal of excitatory amino acids from the synaptic cleft. Here authors combined high-speed atomic force microscopy line-scanning with automated state assignment for the determination of transport dynamics of Glt_Ph, a prokaryotic EAAT homologue, with millisecond temporal resolution.

Glutamate transporters: a broad review of the most recent archaeal and human structures

Glutamate transporters play important roles in bacteria, archaea and eukaryotes. Their function in the mammalian central nervous system is essential for preventing excitotoxicity, and their dysregulation is implicated in many diseases, such as epilepsy and Alzheimer's. Elucidating their transport mechanism would further the understanding of these transporters and promote drug design as they provide compelling targets for understanding the pathophysiology of diseases and may have a direct role in the treatment of conditions involving glutamate excitotoxicity. This review outlines the insights into the transport cycle, uncoupled chloride conductance and modulation, as well as identifying areas that require further investigation.

Free-Energy Simulations Resolve the Low-Affinity Na+-High-Affinity Asp Binding Paradox in GltPh

Glutamate transporters clear up excess extracellular glutamate by cotransporting three Na⁺ and one H⁺ with the countertransport of one K⁺. The archaeal homologs are selective to aspartate and only cotransport three Na⁺. The crystal structures of Glt_Ph from archaea have been used in computational studies to understand the transport mechanism. Although some progress has been made with regard to the ligand-binding sites, a consistent picture of transport still eludes us. A major concern is the discrepancy between the computed binding free energies, which predict high-affinity Na⁺-low-affinity aspartate binding, and the experimental results in which the opposite is observed. Here, we show that the binding of the first two Na⁺ ions involves an intermediate state near the Na1 site, where two Na⁺ ions coexist and couple to aspartate with similar strengths, boosting its affinity. Binding free energies for Na⁺ and aspartate obtained using this intermediate state are in good agreement with the experimental values. Thus, the paradox in binding affinities arises from the assumption that the ligands bind to the sites observed in the crystal structure following the order dictated by their binding free energies with no intermediate states. In fact, the presence of an intermediate state eliminates such a correlation between the binding free energies and the binding order. The intermediate state also facilitates transition of the first Na⁺ ion to its final binding site via a knock-on mechanism, which induces substantial conformational changes in the protein consistent with experimental observations.

Dissociation Process of a MDM2/p53 Complex Investigated by Parallel Cascade Selection Molecular Dynamics and the Markov State Model

Recently, we efficiently generated dissociation pathways of a protein-ligand complex without applying force bias with parallel cascade selection molecular dynamics (PaCS-MD) and showed that PaCS-MD in combination with the Markov state model (MSM) yielded a binding free energy comparable to experimental values. In this work, we applied the same procedure to a complex of MDM2 protein and the transactivation domain of p53 protein (TAD-p53), the latter of which is known to be very flexible in the unbound state. Using 30 independent MD simulations in PaCS-MD, we successfully generated 25 dissociation pathways of the complex, which showed complete or partial unfolding of the helical region of TAD-p53 during the dissociation process within an average simulation time of 154.8 ± 46.4 ns. The standard binding free energy obtained in combination with one-dimensional-, three-dimensional (3D)- or C_α-MSM was in good agreement with those determined experimentally. Using 3D-MSM based on the center of mass position of TAD-p53 relative to MDM2, the dissociation rate constant was calculated, which was comparable to those measured experimentally. C_α-MSM based on all C_α coordinates of TAD-p53 reproduced the experimentally measured standard binding free energy, and dissociation and association rate constants. We conclude that the combination of PaCS-MD and MSM offers an efficient computational procedure to calculate binding free energies and kinetic rates.

Kinetic mechanism of coupled binding in sodium-aspartate symporter GltPh

Many secondary active membrane transporters pump substrates against concentration gradients by coupling their uptake to symport of sodium ions. Symport requires the substrate and ions to be always transported together. Cooperative binding of the solutes is a key mechanism contributing to coupled transport in the sodium and aspartate symporter from Pyrococcus horikoshii Glt_Ph. Here, we describe the kinetic mechanism of coupled binding for Glt_Ph in the inward facing state. The first of the three coupled sodium ions, binds weakly and slowly, enabling the protein to accept the rest of the ions and the substrate. The last ion binds tightly, but is in rapid equilibrium with solution. Its release is required for the complex disassembly. Thus, the first ion serves to ‘open the door’ for the substrate, the last ion ‘locks the door’ once the substrate is in, and one ion contributes to both events.

Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations

The expansion of computational power, better parameterization of force fields, and the development of novel algorithms to enhance the sampling of the free energy landscapes of proteins have allowed molecular dynamics (MD) simulations to become an indispensable tool to understand the function of biomolecules. The temporal and spatial resolution of MD simulations allows for the study of a vast number of processes of interest. Here, we review the computational efforts to uncover the conformational free energy landscapes of a subset of membrane proteins: ion channels, transporters and G-protein coupled receptors. We focus on the various enhanced sampling techniques used to study these questions, how the conclusions come together to build a coherent picture, and the relationship between simulation outcomes and experimental observables.

Elucidation of the Role of a Conserved Methionine in Glutamate Transporters and Its Implication for Force Fields

Glutamate transport through the excitatory amino acid transporters is coupled to the cotransport of three Na⁺ ions, the binding sites (Na1-Na3) of which are conserved from archaea to mammalians. Molecular dynamics (MD) simulations reproduce the Na1 and Na3 binding sites observed in the crystal structures but fail in the case of Na2. A distinguishing feature of the Na2 site is that an S atom from a conserved methionine residue is in the coordination shell of Na⁺. We perform MD simulations on the recent Glt_Tk structure and show that the problem with the Na2 site arises from using an inadequate partial charge for S. When methionine is appropriately parametrized, both the position and the binding free energy of Na⁺ at the Na2 site can be reproduced in good agreement with the experimental data. Other properties of methionine, such as its dipole moment and the solvation energy of its side chain analog, also benefit from this reparametrization. Thus, the Na2 site in glutamate transporters provides a good opportunity for a proper parametrization of methionine in MD force fields.

What Can and Cannot Be Learned from Molecular Dynamics Simulations of Bacterial Proton-Coupled Oligopeptide Transporter GkPOT?

We have performed an extensive set of all-atom molecular dynamics (MD) simulations of a bacterial proton-coupled oligopeptide transporter (POT) in an explicit membrane environment. We have characterized both the local and global conformational dynamics of the transporter upon the proton and/or substrate binding, within a statistical framework. Our results reveal a clearly distinct behavior for local conformational dynamics in the absence and presence of the proton at the putative proton binding residue E310. Particularly, we find that the substrate binding conformation is drastically different in the two conditions, where the substrate binds to the protein in a lateral/vertical manner, in the presence/absence of the proton. We do not observe any statistically significant distinctive behavior in terms of the global conformational changes in different simulation conditions, within the time scales of our simulations. Our extensive simulations and analyses call into question the implicit assumption of many MD studies that local conformational changes observed in short simulations could provide clues to the global conformational changes that occur on much longer time scales. The linear regression analysis of quantities associated with the global conformational fluctuations, however, provides an indication of a mechanism involving the concerted motion of the transmembrane helices, consistent with the rocker-switch mechanism.

Characterization of the Inward- and Outward-Facing Substrate Binding Sites of the Prokaryotic Aspartate Transporter, GltPh

Crystal structures of the prokaryotic aspartate transporter, Glt_Ph, have provided important insights into the mechanism of amino acid transport by Glt_Ph and related eukaryotic members of the glutamate transporter family (SLC1A family). Identification of inhibitors of Glt_Ph can provide valuable tools for understanding the molecular basis for substrate and inhibitor specificity and selectivity of SLC1A members, but at present, few inhibitors of Glt_Ph have been identified. We have screened a collection of commercially available aspartate analogues and identified new transportable and nontransportable Glt_Ph inhibitors. We have explored the inhibition profile of Glt_Ph by utilizing a thiol modification assay that isolates sided populations of the transporters reconstituted in liposomes to determine if any aspartate analogues display a preference for either the inwardly or outwardly directed binding sites. Here, we have characterized several new inhibitors of Glt_Ph and identified three β-carbon-substituted molecules that display a strong preference for the outwardly directed binding site of Glt_Ph.

SLC Transporters: Structure, Function, and Drug Discovery

The human Solute Carrier (SLC) transporters are important targets for drug development. Structure-based drug discovery for SLC transporters requires the description of their structure, dynamics, and mechanism of interaction with small molecule ligands and ions. The recent determination of atomic structures of human SLC transporters and their homologs, combined with improved computational power and prediction methods have led to an increased applicability of structure-based drug design methods for human SLC members. In this review, we provide an overview of the SLC transporters' structures and transport mechanisms. We then describe computational techniques, such as homology modeling and virtual screening that are emerging as key tools to discover chemical probes for human SLC members. We illustrate the utility of these methods by presenting case studies in which rational integration of computation and experiment was used to characterize SLC members that transport key nutrients and metabolites, including the amino acid transporters LAT-1 and ASCT2, the SLC13 family of citric acid cycle intermediate transporters, and the glucose transporter GLUT1. We conclude with a brief discussion about future directions in structure-based drug discovery for the human SLC superfamily, one of the most structurally and functionally diverse protein families in human.

Allosteric Mechanisms of Molecular Machines at the Membrane: Transport by Sodium-Coupled Symporters

Solute transport across cell membranes is ubiquitous in biology as an essential physiological process. Secondary active transporters couple the unfavorable process of solute transport against its concentration gradient to the energetically favorable transport of one or several ions. The study of such transporters over several decades indicates that their function involves complex allosteric mechanisms that are progressively being revealed in atomistic detail. We focus on two well-characterized sodium-coupled symporters: the bacterial amino acid transporter LeuT, which is the prototype for the "gated pore" mechanism in the mammalian synaptic monoamine transporters, and the archaeal GltPh, which is the prototype for the "elevator" mechanism in the mammalian excitatory amino acid transporters. We present the evidence for the role of allostery in the context of a quantitative formalism that can reconcile biochemical and biophysical data and thereby connects directly to recent insights into the molecular structure and dynamics of these proteins. We demonstrate that, while the structures and mechanisms of these transporters are very different, the available data suggest a common role of specific models of allostery in their functions. We argue that such allosteric mechanisms appear essential not only for sodium-coupled symport in general but also for the function of other types of molecular machines in the membrane.

Molecular physiology of EAAT anion channels

Glutamate is the major excitatory neurotransmitter in the mammalian central nervous system. After release from presynaptic nerve terminals, glutamate is quickly removed from the synaptic cleft by a family of five glutamate transporters, the so-called excitatory amino acid transporters (EAAT1-5). EAATs are prototypic members of the growing number of dual-function transport proteins: they are not only glutamate transporters, but also anion channels. Whereas the mechanisms underlying secondary active glutamate transport are well understood at the functional and at the structural level, mechanisms and cellular roles of EAAT anion conduction have remained elusive for many years. Recently, molecular dynamics simulations combined with simulation-guided mutagenesis and experimental analysis identified a novel anion-conducting conformation, which accounts for all experimental data on EAAT anion currents reported so far. We here review recent findings on how EAATs accommodate a transporter and a channel in one single protein.

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.

Refinement of the Central Steps of Substrate Transport by the Aspartate Transporter GltPh: Elucidating the Role of the Na2 Sodium Binding Site

Glutamate homeostasis in the brain is maintained by glutamate transporter mediated accumulation. Impaired transport is associated with several neurological disorders, including stroke and amyotrophic lateral sclerosis. Crystal structures of the homolog transporter GltPh from Pyrococcus horikoshii revealed large structural changes. Substrate uptake at the atomic level and the mechanism of ion gradient conversion into directional transport remained enigmatic. We observed in repeated simulations that two local structural changes regulated transport. The first change led to formation of the transient Na2 sodium binding site, triggered by side chain rotation of T308. The second change destabilized cytoplasmic ionic interactions. We found that sodium binding to the transiently formed Na2 site energized substrate uptake through reshaping of the energy hypersurface. Uptake experiments in reconstituted proteoliposomes confirmed the proposed mechanism. We reproduced the results in the human glutamate transporter EAAT3 indicating a conserved mechanics from archaea to humans.

We used the archaeal homolog GltPh of the human glutamate transporters to refine our understanding how large scale conformational changes are translated into substrate translocation. We identified the structural changes that accompany substrate transport and convert the energy stored in the ion gradient into a directional transport. Insights into the mechanics of these transporters are likely to increase our understanding of how they contribute to excitotoxicity and to develop drugs, which preclude the underlying accumulation of glutamate in the synaptic cleft.

Computational approaches for designing potent and selective analogs of peptide toxins as novel therapeutics

Peptide toxins provide valuable therapeutic leads for many diseases. As they bind to their targets with high affinity, potency is usually ensured. However, toxins also bind to off-target receptors, causing potential side effects. Thus, a major challenge in generating drugs from peptide toxins is ensuring their specificity for their intended targets. Computational methods can play an important role in solving such design problems through construction of accurate models of receptor-toxin complexes and calculation of binding free energies. Here we review the computational methods used for this purpose and their application to toxins targeting ion channels. We describe ShK and HsTX1 toxins, high-affinity blockers of the voltage-gated potassium channel Kv1.3, which could be developed as therapeutic agents for autoimmune diseases.

Calculation of free energy changes due to mutations from alchemical free energy simulations

How a mutation affects the binding free energy of a ligand is a fundamental problem in molecular biology/biochemistry with many applications in pharmacology and biotechnology, e.g. design of drugs and enzymes. Free energy change due to a mutation can be determined most accurately by performing alchemical free energy calculations in molecular dynamics (MD) simulations. Here we discuss the necessary conditions for success of free energy calculations using toxin peptides that bind to ion channels as examples. We show that preservation of the binding mode is an essential requirement but this condition is not always satisfied, especially when the mutation involves a charged residue. Otherwise problems with accuracy of results encountered in mutation of charged residues can be overcome by performing the mutation on the ligand in the binding site and bulk simultaneously and in the same system. The proposed method will be useful in improving the affinity and selectivity profiles of drug leads and enzymes via computational design and protein engineering.

Computational characterization of structural dynamics underlying function in active membrane transporters

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

Molecular dynamics simulations elucidate the mechanism of proton transport in the glutamate transporter EAAT3

The uptake of glutamate in nerve synapses is carried out by the excitatory amino acid transporters (EAATs), involving the cotransport of a proton and three Na(+) ions and the countertransport of a K(+) ion. In this study, we use an EAAT3 homology model to calculate the pKa of several titratable residues around the glutamate binding site to locate the proton carrier site involved in the translocation of the substrate. After identifying E374 as the main candidate for carrying the proton, we calculate the protonation state of this residue in different conformations of EAAT3 and with different ligands bound. We find that E374 is protonated in the fully bound state, but removing the Na2 ion and the substrate reduces the pKa of this residue and favors the release of the proton to solution. Removing the remaining Na(+) ions again favors the protonation of E374 in both the outward- and inward-facing states, hence the proton is not released in the empty transporter. By calculating the pKa of E374 with a K(+) ion bound in three possible sites, we show that binding of the K(+) ion is necessary for the release of the proton in the inward-facing state. This suggests a mechanism in which a K(+) ion replaces one of the ligands bound to the transporter, which may explain the faster transport rates of the EAATs compared to its archaeal homologs.

A consistent force field parameter set for zwitterionic amino acid residues

Isolated amino acids play an important role in biochemistry and are therefore an interesting object of study. Atomistic molecular dynamics (MD) simulations can provide a high-resolution picture of the dynamic features of these species, especially in their biological environment. Unfortunately, most standard force field packages lack libraries for isolated amino acids in their zwitterionic form. Although several studies have used ad-hoc parameterizations for single amino acids, a consistent force-field parameter set for these molecules is still missing. Here, we present such a parameter library derived from the widely used parm99SB set from the AMBER program package. The parameter derivation for all 20 proteinogenic amino acids transparently followed established procedures with histidine treated in three different protonation states. All amino acids were subjected to MD simulations in four different forms for comparison: zwitterionic, N-teminally capped with acetyl, C-terminally capped with N-methyl, and capped at both termini. Simulation results show similarities between the different forms. Five zwitterionic amino acids-arginine, glutamate, glycine, phenylalanine, leucine-were simulated in a protein environment. Proteins and ligands generally retained their initial structure. The new parameter set will thus facilitate future atomistic simulations of these species.

Coupled ion binding and structural transitions along the transport cycle of glutamate transporters

Membrane transporters that clear the neurotransmitter glutamate from synapses are driven by symport of sodium ions and counter-transport of a potassium ion. Previous crystal structures of a homologous archaeal sodium and aspartate symporter showed that a dedicated transport domain carries the substrate and ions across the membrane. Here, we report new crystal structures of this homologue in ligand-free and ions-only bound outward- and inward-facing conformations. We show that after ligand release, the apo transport domain adopts a compact and occluded conformation that can traverse the membrane, completing the transport cycle. Sodium binding primes the transport domain to accept its substrate and triggers extracellular gate opening, which prevents inward domain translocation until substrate binding takes place. Furthermore, we describe a new cation-binding site ideally suited to bind a counter-transported ion. We suggest that potassium binding at this site stabilizes the translocation-competent conformation of the unloaded transport domain in mammalian homologues.

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

Molecules of glutamate can carry messages between cells in the brain, and these signals are essential for thought and memory. Glutamate molecules can also act as signals to build new connections between brain cells and to prune away unnecessary ones. However, too much glutamate outside of the cells kills the brain tissue and can lead to devastating brain diseases.

In a healthy brain, special pumps called glutamate transporters move these molecules back into the brain cells, where they can be stored safely. However, when brain cells are damaged—by, for example, a stroke or an injury,—the glutamate stored inside spills out, killing the surrounding cells. This leads to a cascade of dying cells and leaking glutamate, which causes even more damage and slows the recovery.

Glutamate transporters ensure that there are more glutamate molecules inside cells than outside. However, it requires energy to maintain this gradient in the concentration of glutamate molecules. The transporters get this energy by moving three sodium ions into the cell with each glutamate molecule, and moving one potassium ion out of the cell. However, it is not clear how these transporters ensure that they move the glutamate molecules and the sodium ions at the same time.

Now, Verdon, Oh et al. have uncovered the 3D structure of a glutamate transporter homologue at each step of the transport process. These structures reveal that, on the outside of the cell membrane, sodium ions attach to the so-called ‘transporter domain’ and make it better able to bind glutamate. The transporter domain then carries the sodium ions and glutamate through the cell membrane and releases them into the cell. Verdon, Oh et al. suggest that a potassium ion then binds to the empty transport domain, stabilizing it into a more compact shape that easily makes the return trip to the outside of the cell.

Most experiments on glutamate transporters, including the work of Verdon, Oh et al., are carried out on model proteins taken from bacteria. An important challenge for the future will be to obtain structural information on human glutamate transporters, as these could be therapeutic targets for the treatment of various neurological conditions.

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

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.

Computation of standard binding free energies of polar and charged ligands to the glutamate receptor GluA2

Accurate calculation of the binding affinity of small molecules to proteins has the potential to become an important tool in rational drug design. In this study, we use the free energy perturbation (FEP) method with restraints to calculate the standard binding free energy of five ligands (ACPA, AMPA, CNQX, DNQX, and glutamate) to the glutamate receptor GluA2, which plays an essential role in synaptic transmission. To deal with the convergence problem in FEP calculations with charged ligands, we use a protocol where the ligand is coupled in the binding site while it is decoupled in bulk solution simultaneously. The contributions from the conformational, rotational, and translational entropies to the standard binding free energy are determined by applying/releasing respective restraints to the ligand in bulk/binding site. We also employ the confine-and-release approach, which helps to resolve convergence problems in FEP calculations. Our results are in good agreement with the experimental values for all five ligands, including the charged ones which are often problematic in FEP calculations. We also analyze the different contributions to the binding free energy of each ligand to GluA2 and discuss the nature of these interactions.