Na+ interactions with the neutral amino acid transporter ASCT1.

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.

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.

BAT2: an Open-Source Tool for Flexible, Automated, and Low Cost Absolute Binding Free Energy Calculations.

Absolute binding free energy (ABFE) calculations with all-atom molecular dynamics (MD) have the potential to greatly reduce costs in the first stages of drug discovery. Here, we introduce BAT2, the new version of the Binding Affinity Tool (BAT.py), designed to combine full automation of ABFE calculations with high-performance MD simulations, making it a potential tool for virtual screening. We describe and test several changes and new features that were incorporated into the code, such as relative restraints between the protein and the ligand instead of using fixed dummy atoms, support for the OpenMM simulation engine, a merged approach to the application/release of restraints, support for cobinders and proteins with multiple chains, and many others. We also reduced the simulation times for each ABFE calculation, assessing the effect on the expected robustness and accuracy of the calculations.

Investigating the Solvent Effects on Binding Affinity of PAHs-ExBox4+ Complexes: An Alchemical Approach.

Polycyclic aromatic hydrocarbons (PAHs) are polluting agents, produced naturally or artificially, widely dispersed in the environment and potentially carcinogenic and immunotoxic to humans and animals, mainly for marine life. Recently, a tetracationic box-shaped cyclophane (ExBox4+) was synthesized, fully characterized, and revealed to form host-guest complexes with PAHs in acetonitrile, demonstrating the potential ability for it to act as a PAHs scavenger. This work investigates, through Molecular Dynamics (MD) simulations, the binding affinity between different PAHs and ExBox4+ in different solvents: chloroform (nonpolar), acetonitrile (polar protic), and water (polar protic). An alchemical method of simultaneous decoupling-recoupling (SDR) was used and implemented in a newly developed Python program called GHOAT, which fully automates the calculation of binding free energies and invokes the AMBER 2020 simulation package. The results showed that the affinity between ExBox4+ and PAHs in water is much larger than in organic media, with free energies between -5 and -20 kcal/mol, being able to act as a PHAs scavenger with great potential for applications in environmental chemistry such as soil washing. The results also reveal a significant correlation with the experimental available ΔG values. The methodology employed presents itself as an important tool for the in silico determination of binding affinities, not only available for charged cyclophanes but also extensible to several other HG supramolecular systems in condensed media, aiding in the rational design of host-guest systems in a significant way.

Absolute binding free energy calculations improve enrichment of actives in virtual compound screening.

We determined the effectiveness of absolute binding free energy (ABFE) calculations to refine the selection of active compounds in virtual compound screening, a setting where the more commonly used relative binding free energy approach is not readily applicable. To do this, we conducted baseline docking calculations of structurally diverse compounds in the DUD-E database for three targets, BACE1, CDK2 and thrombin, followed by ABFE calculations for compounds with high docking scores. The docking calculations alone achieved solid enrichment of active compounds over decoys. Encouragingly, the ABFE calculations then improved on this baseline. Analysis of the results emphasizes the importance of establishing high quality ligand poses as starting points for ABFE calculations, a nontrivial goal when processing a library of diverse compounds without informative co-crystal structures. Overall, our results suggest that ABFE calculations can play a valuable role in the drug discovery process.

Free energy simulations of ligand binding to the aspartate transporter Glt(Ph).

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.

Automation of absolute protein-ligand binding free energy calculations for docking refinement and compound evaluation.

Absolute binding free energy calculations with explicit solvent molecular simulations can provide estimates of protein-ligand affinities, and thus reduce the time and costs needed to find new drug candidates. However, these calculations can be complex to implement and perform. Here, we introduce the software BAT.py, a Python tool that invokes the AMBER simulation package to automate the calculation of binding free energies for a protein with a series of ligands. The software supports the attach-pull-release (APR) and double decoupling (DD) binding free energy methods, as well as the simultaneous decoupling-recoupling (SDR) method, a variant of double decoupling that avoids numerical artifacts associated with charged ligands. We report encouraging initial test applications of this software both to re-rank docked poses and to estimate overall binding free energies. We also show that it is practical to carry out these calculations cheaply by using graphical processing units in common machines that can be built for this purpose. The combination of automation and low cost positions this procedure to be applied in a relatively high-throughput mode and thus stands to enable new applications in early-stage drug discovery.

Shedding Light on the Hydrolysis Mechanism of cis, trans-[Ru(dmso)4Cl2] Complexes and Their Interactions with DNA-A Computational Perspective.

Using several computational tools such as density functional theory analysis, docking, and MD simulations, we performed a study on cis, trans-[Ru(II)(dmso)4Cl2] complexes, which have therapeutic potential as antimetastatic agents, and their association with DNA. Kohn-Sham energy decomposition analysis reveals that dmso ligands have much smaller interaction energies compared to the chlorido ligands, and their substitution by aquo ligands induces an extra stabilization of the other metal-ligand bonds. Once the complex is hydrolyzed, the aquo ligands have the weakest interactions to the metallic center and therefore are more labile for substitution by a DNA atom. Molecular docking and molecular dynamics were employed to understand the complex preassociation to DNA, pointing to a higher affinity of the hydrolyzed complexes, as well as showing spontaneous binding events during the simulations. Our results are consistent with the experimentally available data that suggest a mechanism in which the complexes are quickly hydrolyzed in solution, before forming cross-links with the DNA molecule. We present a set of methods that could be used to optimize these complexes computationally, aiding in the development of new drugs based on transition metals.

Attach-Pull-Release Calculations of Ligand Binding and Conformational Changes on the First BRD4 Bromodomain.

Bromodomains, protein domains involved in epigenetic regulation, are able to bind small molecules with high affinity. In the present study, we report free energy calculations for the binding of seven ligands to the first BRD4 bromodomain, using the attach-pull-release (APR) method to compute the reversible work of removing the ligands from the binding site and then allowing the protein to relax conformationally. We test three different water models, TIP3P, TIP4PEw, and SPC/E, as well as the GAFF and GAFF2 parameter sets for the ligands. Our simulations show that the apo crystal structure of BRD4 is only metastable, with a structural transition happening in the absence of the ligand typically after 20 ns of simulation. We compute the free energy change for this transition with a separate APR calculation on the free protein and include its contribution to the ligand binding free energies, which generally causes an underestimation of the affinities. By testing different water models and ligand parameters, we are also able to assess their influence in our results and determine which one produces the best agreement with the experimental data. Both free energies associated with the conformational change and ligand binding are affected by the choice of water model, with the two sets of ligand parameters affecting their binding free energies to a lesser degree. Across all six combinations of water model and ligand potential function, the Pearson correlation coefficients between calculated and experimental binding free energies range from 0.55 to 0.83, and the root-mean-square errors range from 1.4-3.2 kcal/mol. The current protocol also yields encouraging preliminary results when used to assess the relative stability of ligand poses generated by docking or other methods, as illustrated for two different ligands. Our method takes advantage of the high performance provided by graphics processing units and can readily be applied to other ligands as well as other protein systems.

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.

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.

Mechanism and energetics of ligand release in the aspartate transporter GltPh.

The bacterial aspartate transporter GltPh cotransports three Na(+) ions with the substrate. The mechanism and energetics of ligand binding have previously been studied using molecular dynamics simulations on the crystal structure of GltPh captured in the outward-facing state. Here we use the recent crystal structure of the inward-facing state of GltPh to study the reverse process of unbinding of ligands. Gating behavior is studied in the presence of different ligands. A detailed characterization of the intracellular gate is given, pointing out the differences from the extracellular gate. We then perform free energy simulations to calculate the binding affinities of all the ligands in different combinations, from which the unbinding order is determined as Na2, (gate opens), Asp, Na1, and Na3. The strong coupling between Asp and Na1 is quantified from several free energy calculations. Na3 has the largest affinity to GltPh, and therefore, its unbinding is proposed as the rate-limiting step in the transport cycle. The release time of Na3, estimated from Kramers' rate theory, is shown to be consistent with the experimental turnover rate of the transporter.

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.

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, Glt(Ph), 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 Glt(Ph). In order to fully utilize the Glt(Ph) structure in functional studies of the human glutamate transporters, it is essential to understand the transport mechanism of Glt(Ph) 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 Glt(Ph) 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 Glt(Ph), 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 Glt(Ph) and EAAT1.