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 Glt(ph) 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.