Germano Heinzelmann

Position of the Third Na+ Site in the Aspartate Transporter GltPh and the Human Glutamate Transporter, EAAT1

Glutamate transport via the human excitatory amino acid transporters is coupled to the co-transport of three Na+ ions, one H+ and the counter-transport of one K+ ion. Transport by an archaeal homologue of the human glutamate transporters, GltPh, whose three dimensional structure is known is also coupled to three Na+ ions but only two Na+ ion binding sites have been observed in the crystal structure of GltPh. In order to fully utilize the GltPh structure in functional studies of the human glutamate transporters, it is essential to understand the transport mechanism of GltPh and accurately determine the number and location of Na+ ions coupled to transport. Several sites have been proposed for the binding of a third Na+ ion from electrostatic calculations and molecular dynamics simulations. In this study, we have performed detailed free energy simulations for GltPh and reveal a new site for the third Na+ ion involving the side chains of Threonine 92, Serine 93, Asparagine 310, Aspartate 312, and the backbone of Tyrosine 89. We have also studied the transport properties of alanine mutants of the coordinating residues Threonine 92 and Serine 93 in GltPh, and the corresponding residues in a human glutamate transporter, EAAT1. The mutant transporters have reduced affinity for Na+ compared to their wild type counterparts. These results confirm that Threonine 92 and Serine 93 are involved in the coordination of the third Na+ ion in GltPh and EAAT1.

A Potent and Selective Peptide Blocker of the Kv1.3 Channel: Prediction from Free-Energy Simulations and Experimental Confirmation

The voltage-gated potassium channel Kv1.3 is a well-established target for treatment of autoimmune diseases. ShK peptide from a sea anemone is one of the most potent blockers of Kv1.3 but its application as a therapeutic agent for autoimmune diseases is limited by its lack of selectivity against other Kv channels, in particular Kv1.1. Accurate models of Kv1.x-ShK complexes suggest that specific charge mutations on ShK could considerably enhance its specificity for Kv1.3. Here we evaluate the K18A mutation on ShK, and calculate the change in binding free energy associated with this mutation using the path-independent free energy perturbation and thermodynamic integration methods, with a novel implementation that avoids convergence problems. To check the accuracy of the results, the binding free energy differences were also determined from path-dependent potential of mean force calculations. The two methods yield consistent results for the K18A mutation in ShK and predict a 2 kcal/mol gain in Kv1.3/Kv1.1 selectivity free energy relative to wild-type peptide. Functional assays confirm the predicted selectivity gain for ShK[K18A] and suggest that it will be a valuable lead in the development of therapeutics for autoimmune diseases.

Free Energy Simulations of Ligand Binding to the Aspartate Transporter GltPh

Glutamate/Aspartate transporters cotransport three Na(+) and one H(+) ions with the substrate and countertransport one K(+) ion. The binding sites for the substrate and two Na(+) ions have been observed in the crystal structure of the archeal homolog Glt(Ph), while the binding site for the third Na(+) ion has been proposed from computational studies and confirmed by experiments. Here we perform detailed free energy simulations of Glt(Ph), giving a comprehensive characterization of the substrate and ion binding sites, and calculating their binding free energies in various configurations. Our results show unequivocally that the substrate binds after the binding of two Na(+) ions. They also shed light into Asp/Glu selectivity of Glt(Ph), which is not observed in eukaryotic glutamate transporters.

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.

Na+ Interactions with the Neutral Amino Acid Transporter ASCT1

Background: Na+ interactions with the alanine serine cysteine transporters (ASCT) have not been well established. Results: Identified residues involved in coordinating three Na+ ions. Conclusion: ASCT1-mediated transport requires Na+ bound to only two of three Na+ sites; activation of the anion conductance requires only one Na+ bound. Significance: ASCT1 has three Na+ sites, however, full occupancy is not required for function. The alanine, serine, cysteine transporters (ASCTs) belong to the solute carrier family 1A (SLC1A), which also includes the excitatory amino acid transporters (EAATs) and the prokaryotic aspartate transporter GltPh. Acidic amino acid transport by the EAATs is coupled to the co-transport of three Na+ ions and one proton, and the counter-transport of one K+ ion. In contrast, neutral amino acid exchange by the ASCTs does not require protons or the counter-transport of K+ ions and the number of Na+ ions required is not well established. One property common to SLC1A family members is a substrate-activated anion conductance. We have investigated the number and location of Na+ ions required by ASCT1 by mutating residues in ASCT1 that correspond to residues in the EAATs and GltPh that are involved in Na+ binding. Mutations to all three proposed Na+ sites influence the binding of substrate and/or Na+, or the rate of substrate exchange. A G422S mutation near the Na2 site reduced Na+ affinity, without affecting the rate of exchange. D467T and D467A mutations in the Na1 site reduce Na+ and substrate affinity and also the rate of substrate exchange. T124A and D380A mutations in the Na3 site selectively reduce the affinity for Na+ and the rate of substrate exchange without affecting substrate affinity. In many of the mutants that reduce the rate of substrate transport the amplitudes of the substrate-activated anion conductances are not substantially affected indicating altered ion dependence for channel activation compared with substrate exchange.

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.

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.

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.

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.