We use the well-known structural and functional properties of the gramicidin A channel to test the appropriateness of force fields commonly used in molecular dynamics (MD) simulations of ion channels. For this purpose, the high-resolution structure of the gramicidin A dimer is embedded in a dimyristoylphosphatidylcholine bilayer, and the potential of mean force of a K(+) ion is calculated along the channel axis using the umbrella sampling method. Calculations are performed using two of the most common force fields in MD simulations: CHARMM and GROMACS. Both force fields lead to large central barriers for K(+) ion permeation, that are substantially higher than those deduced from the physiological data by inverse methods. In long MD simulations lasting over 60 ns, several ions are observed to enter the binding site but none of them crossed the channel despite the presence of a large driving field. The present results, taken together with many earlier studies, highlights the shortcomings of the standard force fields used in MD simulations of ion channels and calls for construction of more appropriate force fields for this purpose.
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.
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.
Articles you may be interested inTime-dependent four-component relativistic density-functional theory for excitation energies. II. The exchangecorrelation kernel J. Chem. Phys. 123, 054102 (2005); 10.1063/1.1940609 Scalar relativistic all-electron density functional calculations on periodic systems J. Chem. Phys. 122, 084108 (2005); 10.1063/1.1851973 Response to "Comment on 'Four-component relativistic density functional calculations of heavy diatomic molecules'" [Comment on "Four-component relativistic density functional calculations of heavy diatomic molecules" [J. Chem.Molecular density functional calculations in the regular relativistic approximation: Method, application to coinage metal diatomics, hydrides, fluorides and chlorides, and comparison with first-order relativistic calculationsWe perform accurate four-component calculations for heavy closed-shell diatomic molecules in the framework of relativistic density functional theory using local and gradient corrected density functional schemes. As examples we have chosen Cu 2 , Ag 2 , Au 2 , Tl 2 , Pb 2 , Bi 2 , and Pt 2 . The potential energy curves show the quality, and the discrepancies of the density functionals unscreened from any approximation of the relativistic effects.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.