Browse
Submit Data
Databases
API
Help

Dataset Information

21 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

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

ABSTRACT: 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.

SUBMITTER: Wang J

PROVIDER: S-EPMC8064580 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Binding Analysis Using Accelerated Molecular Dynamics Simulations and Future Perspectives.

Project description:Biomolecular recognition such as binding of small molecules, nucleic acids, peptides and proteins to their target receptors plays key roles in cellular function and has been targeted for therapeutic drug design. Molecular dynamics (MD) is a computational approach to analyze these binding processes at an atomistic level, which provides valuable understandings of the mechanisms of biomolecular recognition. However, the rather slow biomolecular binding events often present challenges for conventional MD (cMD), due to limited simulation timescales (typically over hundreds of nanoseconds to tens of microseconds). In this regard, enhanced sampling methods, particularly accelerated MD (aMD), have proven useful to bridge the gap and enable all-atom simulations of biomolecular binding events. Here, we will review the recent method developments of Gaussian aMD (GaMD), ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD), which have greatly expanded our capabilities to simulate biomolecular binding processes. Spontaneous binding of various biomolecules to their receptors has been successfully simulated by GaMD. Microsecond LiGaMD and Pep-GaMD simulations have captured repetitive binding and dissociation of small-molecule ligands and highly flexible peptides, and thus enabled ligand/peptide binding thermodynamics and kinetics calculations. We will also present relevant application studies in simulations of important drug targets and future perspectives for rational computer-aided drug design.

| S-EPMC8747661 | biostudies-literature

Monoamine transporters: insights from molecular dynamics simulations.

Project description:The human monoamine transporters (MATs) facilitate the reuptake of the neurotransmitters serotonin, dopamine, and norepinephrine from the synaptic cleft. Imbalance in monoaminergic neurotransmission is linked to various diseases including major depression, attention deficit hyperactivity disorder, schizophrenia, and Parkinson's disease. Inhibition of the MATs is thus an important strategy for treatment of such diseases. The MATs are sodium-coupled transport proteins belonging to the neurotransmitter/Na(+) symporter (NSS) family, and the publication of the first high-resolution structure of a NSS family member, the bacterial leucine transporter LeuT, in 2005, proved to be a major stepping stone for understanding this family of transporters. Structural data allows for the use of computational methods to study the MATs, which in turn has led to a number of important discoveries. The process of substrate translocation across the membrane is an intrinsically dynamic process. Molecular dynamics simulations, which can provide atomistic details of molecular motion on ns to ms timescales, are therefore well-suited for studying transport processes. In this review, we outline how molecular dynamics simulations have provided insight into the large scale motions associated with transport of the neurotransmitters, as well as the presence of external and internal gates, the coupling between ion and substrate transport, and differences in the conformational changes induced by substrates and inhibitors.

| S-EPMC4607855 | biostudies-literature

Accelerated molecular dynamics simulations of protein folding.

Project description:Folding of four fast-folding proteins, including chignolin, Trp-cage, villin headpiece and WW domain, was simulated via accelerated molecular dynamics (aMD). In comparison with hundred-of-microsecond timescale conventional molecular dynamics (cMD) simulations performed on the Anton supercomputer, aMD captured complete folding of the four proteins in significantly shorter simulation time. The folded protein conformations were found within 0.2-2.1 Å of the native NMR or X-ray crystal structures. Free energy profiles calculated through improved reweighting of the aMD simulations using cumulant expansion to the second-order are in good agreement with those obtained from cMD simulations. This allows us to identify distinct conformational states (e.g., unfolded and intermediate) other than the native structure and the protein folding energy barriers. Detailed analysis of protein secondary structures and local key residue interactions provided important insights into the protein folding pathways. Furthermore, the selections of force fields and aMD simulation parameters are discussed in detail. Our work shows usefulness and accuracy of aMD in studying protein folding, providing basic references in using aMD in future protein-folding studies.

| S-EPMC4487363 | biostudies-literature

Molecular Determinants of Substrate Specificity in Sodium-coupled Glutamate Transporters.

Project description:Crystal structures of the archaeal homologue GltPh have provided important insights into the molecular mechanism of transport of the excitatory neurotransmitter glutamate. Whereas mammalian glutamate transporters can translocate both glutamate and aspartate, GltPh is only one capable of aspartate transport. Most of the amino acid residues that surround the aspartate substrate in the binding pocket of GltPh are highly conserved. However, in the brain transporters, Thr-352 and Met-362 of the reentrant hairpin loop 2 are replaced by the smaller Ala and Thr, respectively. Therefore, we have studied the effects of T352A and M362T on binding and transport of aspartate and glutamate by GltPh. Substrate-dependent intrinsic fluorescence changes were monitored in transporter constructs containing the L130W mutation. GltPh-L130W/T352A exhibited an ~15-fold higher apparent affinity for l-glutamate than the wild type transporter, and the M362T mutation resulted in an increased affinity of ~40-fold. An even larger increase of the apparent affinity for l-glutamate, around 130-fold higher than that of wild type, was observed with the T352A/M362T double mutant. Radioactive uptake experiments show that GltPh-T352A not only transports aspartate but also l-glutamate. Remarkably, GltPh-M362T exhibited l-aspartate but not l-glutamate transport. The double mutant retained the ability to transport l-glutamate, but its kinetic parameters were very similar to those of GltPh-T352A alone. The differential impact of mutation on binding and transport of glutamate suggests that hairpin loop 2 not only plays a role in the selection of the substrate but also in its translocation.

| S-EPMC4661411 | biostudies-literature

Coupling Accelerated Molecular Dynamics Methods with Thermodynamic Integration Simulations.

Project description:In this work we propose a straightforward and efficient approach to improve accuracy and convergence of free energy simulations in condensed-phase systems. We also introduce a new accelerated Molecular Dynamics (MD) approach in which molecular conformational transitions are accelerated by lowering the energy barriers while the potential surfaces near the minima are left unchanged. All free energy calculations were performed on the propane-to-propane model system. The accuracy of free energy simulations was significantly improved when sampling of internal degrees of freedom of solute was enhanced. However, accurate and converged results were only achieved when the solvent interactions were taken into account in the accelerated MD approaches. The analysis of the distribution of boost potential along the free energy simulations showed that the new accelerated MD approach samples efficiently both low- and high-energy regions of the potential surface. Since this approach also maintains substantial populations in regions near the minima, the statistics are not compromised in the thermodynamic integration calculations, and, as a result, the ensemble average can be recovered.

| S-EPMC2646661 | biostudies-other

Accelerated molecular dynamics simulations of ligand binding to a muscarinic G-protein-coupled receptor.

Project description:Elucidating the detailed process of ligand binding to a receptor is pharmaceutically important for identifying druggable binding sites. With the ability to provide atomistic detail, computational methods are well poised to study these processes. Here, accelerated molecular dynamics (aMD) is proposed to simulate processes of ligand binding to a G-protein-coupled receptor (GPCR), in this case the M3 muscarinic receptor, which is a target for treating many human diseases, including cancer, diabetes and obesity. Long-timescale aMD simulations were performed to observe the binding of three chemically diverse ligand molecules: antagonist tiotropium (TTP), partial agonist arecoline (ARc) and full agonist acetylcholine (ACh). In comparison with earlier microsecond-timescale conventional MD simulations, aMD greatly accelerated the binding of ACh to the receptor orthosteric ligand-binding site and the binding of TTP to an extracellular vestibule. Further aMD simulations also captured binding of ARc to the receptor orthosteric site. Additionally, all three ligands were observed to bind in the extracellular vestibule during their binding pathways, suggesting that it is a metastable binding site. This study demonstrates the applicability of aMD to protein-ligand binding, especially the drug recognition of GPCRs.

| S-EPMC5435230 | biostudies-literature

Improved Modeling of Peptide-Protein Binding Through Global Docking and Accelerated Molecular Dynamics Simulations.

Project description:Peptides mediate up to 40% of known protein-protein interactions in higher eukaryotes and play a key role in cellular signaling, protein trafficking, immunology, and oncology. However, it is challenging to predict peptide-protein binding with conventional computational modeling approaches, due to slow dynamics and high peptide flexibility. Here, we present a prototype of the approach which combines global peptide docking using ClusPro PeptiDock and all-atom enhanced simulations using Gaussian accelerated molecular dynamics (GaMD). For three distinct model peptides, the lowest backbone root-mean-square deviations (RMSDs) of their bound conformations relative to X-ray structures obtained from PeptiDock were 3.3-4.8 Å, being medium quality predictions according to the Critical Assessment of PRediction of Interactions (CAPRI) criteria. GaMD simulations refined the peptide-protein complex structures with significantly reduced peptide backbone RMSDs of 0.6-2.7 Å, yielding two high quality (sub-angstrom) and one medium quality models. Furthermore, the GaMD simulations identified important low-energy conformational states and revealed the mechanism of peptide binding to the target proteins. Therefore, PeptiDock+GaMD is a promising approach for exploring peptide-protein interactions.

| S-EPMC6835073 | biostudies-literature

Trypsinogen activation as observed in accelerated molecular dynamics simulations.

Project description:Serine proteases are involved in many fundamental physiological processes, and control of their activity mainly results from the fact that they are synthetized in an inactive form that becomes active upon cleavage. Three decades ago Martin Karplus's group performed the first molecular dynamics simulations of trypsin, the most studied member of the serine protease family, to address the transition from the zymogen to its active form. Based on the computational power available at the time, only high frequency fluctuations, but not the transition steps, could be observed. By performing accelerated molecular dynamics (aMD) simulations, an interesting approach that increases the configurational sampling of atomistic simulations, we were able to observe the N-terminal tail insertion, a crucial step of the transition mechanism. Our results also support the hypothesis that the hydrophobic effect is the main force guiding the insertion step, although substantial enthalpic contributions are important in the activation mechanism. As the N-terminal tail insertion is a conserved step in the activation of serine proteases, these results afford new perspective on the underlying thermodynamics of the transition from the zymogen to the active enzyme.

| S-EPMC4241106 | biostudies-literature

Enhancing protein adsorption simulations by using accelerated molecular dynamics.

Project description:The atomistic modeling of protein adsorption on surfaces is hampered by the different time scales of the simulation ([Formula: see text][Formula: see text]s) and experiment (up to hours), and the accordingly different 'final' adsorption conformations. We provide evidence that the method of accelerated molecular dynamics is an efficient tool to obtain equilibrated adsorption states. As a model system we study the adsorption of the protein BMP-2 on graphite in an explicit salt water environment. We demonstrate that due to the considerably improved sampling of conformational space, accelerated molecular dynamics allows to observe the complete unfolding and spreading of the protein on the hydrophobic graphite surface. This result is in agreement with the general finding of protein denaturation upon contact with hydrophobic surfaces.

| S-EPMC3670854 | biostudies-literature

Molecular dynamics simulations of the mammalian glutamate transporter EAAT3.

Project description: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.

| S-EPMC3958442 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data