LeuT conformational sampling utilizing accelerated molecular dynamics and principal component analysis

Allosteric modulation of serotonin and dopamine transporters: New insights from computations and experiments

Human monoamine transporters (MATs) are critical to regulating monoaminergic neurotransmission by translocating their substrates from the synaptic space back into the presynaptic neurons. As such, their primary substrate binding site S1 has been targeted by a wide range of compounds for treating neuropsychiatric and neurodegenerative disorders including depression, ADHD, neuropathic pain, and anxiety disorders. We present here a comparative study of the structural dynamics and ligand-binding properties of two MATs, dopamine transporter (DAT) and serotonin transporter (SERT), with focus on the allosteric modulation of their transport function by drugs or substrates that consistently bind a secondary site S2, proposed to serve as an allosteric site. Our systematic analysis of the conformational space and dynamics of a dataset of 50 structures resolved for DAT and SERT in the presence of one or more ligands/drugs reveals the specific residues playing a consistent role in coordinating the small molecules bound to subsites S2-I and S2-II within S2, such as R476 and Y481 in dDAT and E494, P561, and F556 in hSERT. Further analysis reveals how DAT and SERT differ in their two principal modes of structural changes, PC1 and PC2. Notably, PC1 underlies the transition between outward- and inward-facing states of the transporters as well as their gating; whereas PC2 supports the rearrangements of TM helices near the S2 site. Finally, the examination of cross-correlations between structural elements lining the respective sites S1 and S2 point to the crucial role of coupled motions between TM6a and TM10. In particular, we note the involvement of hSERT residues F335 and G338, and E493-E494-T497 belonging to these two respective helices, in establishing the allosteric communication between S1 and S2. These results help understand the molecular basis of the action of drugs that bind to the S2 site of DAT or SERT. They also provide a basis for designing allosteric modulators that may provide better control of specific interactions and cellular pathways, rather than indiscriminately inhibiting the transporter by targeting its orthosteric site.

Occlusion of the human serotonin transporter is mediated by serotonin-induced conformational changes in the bundle domain

The human serotonin transporter (hSERT) terminates neurotransmission by removing serotonin (5HT) from the synaptic cleft, an essential process for proper functioning of serotonergic neurons. Structures of the hSERT have revealed its molecular architecture in four conformations, including the outward-open and occluded states, and show the transporter’s engagement with co-transported ions and the binding mode of inhibitors. In this study, we investigated the molecular mechanism by which the hSERT occludes and sequesters the substrate 5HT. This first step of substrate uptake into cells is a structural change consisting of the transition from the outward-open to the occluded state. Inhibitors such as the antidepressants citalopram, fluoxetine, and sertraline inhibit this step of the transport cycle. Using molecular dynamics simulations, we reached a fully occluded state, in which the transporter-bound 5HT becomes fully shielded from both sides of the membrane by two closed hydrophobic gates. Analysis of 5HT-triggered occlusion showed that bound 5HT serves as an essential trigger for transporter occlusion. Moreover, simulations revealed a complex sequence of steps and showed that movements of bundle domain helices are only partially correlated. 5HT-triggered occlusion is initially dominated by movements of transmembrane helix 1b, while in the final step, only transmembrane helix 6a moves and relaxes an intermediate change in its secondary structure.

Dynamics and conformational propensities of staphylococcal CntA

Enzymes use transition metals as co-factors for catalytic roles in biological processes. Notably, manganese, iron, cobalt, nickel, copper and zinc are abundantly used. Staphylococcus aureus, a commensal bacterium asymptomatically, lies on the human body causing variety of infections. S. aureus is equipped by advanced virulence-regulatory circuits of metal acquisition like Cnt that acquires metals at infection sites by utilizing a nicotianamine-like metallophore staphylopine. Despite significant growth in structural studies, how CntA of Cnt system transmits conformational signal upon staphylopine recognition remains elusive. Here, we analyzed the structural changes adopted by CntA during close-to-open transition by computational approaches. CntA uses a bi-domain architectural form of domain II which performed 37° rigid body rotation and 1.1 Å translation assisted by inter-domain hinge cluster residues. Important clustered communities were found regulating the conformational changes in CntA where communities 4 and 5 are found crucial. Besides open and close states, the fluctuating regions sampled two additional intermediate states which were considered close or open previously. CntA prefers fluctuating the non-conserved regions rather than conserved where domain II turned out to be rigid and maintains a stable fold. Overall, the CntA system is a potential target for structural biologist to hamper such conformational behaviors at family level.Communicated by Ramaswamy H. Sarma.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Mode-Selective Raman Imaging of Dopamine-Human Dopamine Transporter Interaction in Live Cells

Dopamine (DA) is the catecholamine neurotransmitter which interacts with dopamine receptors (DARs) to generate dopaminergic signals in the nervous system. Dopamine transporter (DAT) interacts with DA to maintain DA's homeostasis in synaptic and perisynaptic space. DAT and DARs have great importance in the central nervous system (CNS) because they are associated with the targeted binding of drugs. Interactions of DA, its analogue with DARs, or DAT have been studied extensively to understand the mechanism of the dopaminergic signaling process and several neurodegenerative diseases, including schizophrenia, Parkinson's diseases, addiction, attention deficit hyperactivity disorder, and bipolar disorder. However, there is still a lack of a risk-free, label-free, and minimally invasive imaging approach to probe the interaction between DA and DAT or DARs. Here, we probed the DA, human dopamine transporter (hDAT), and DA-hDAT interactions in live cells using combined approach of two-photon excited (2PE) fluorescence imaging and mode-selective Raman measurement. We utilized the signature Raman peak at 1287 cm^-1 to probe the location of DA and 807 and 1076 cm^-1 to probe the DA-hDAT interaction in live cells. We found that the combined approach of mode-selective Raman imaging, 2PE fluorescence imaging, and computational methods is successful to probe and confirm the DA-hDAT interactions in living cells. The probing of the interactions of DARs or DAT with DA or other targeting drugs is crucial for the diagnosis and cure of several neurodegenerative diseases. Also, this analytical approach could be extended to probe other types of protein-ligand interactions.

Raman Spectroscopic Analysis of Signaling Molecules-Dopamine Receptors Interactions in Living Cells

The selective interaction of signaling compounds including neurotransmitters and drugs with the dopamine receptors (DARs) is extremely important for the treatment of neurodegenerative diseases. Here, we report a method to probe the selective interactions of signaling compounds with D1 and D2 DARs in living cells using the combined approach of theoretical calculation and surface-enhanced Raman spectroscopy (SERS). When signaling compounds such as DA, amphetamine, methamphetamine, and methylenedioxypyrovalerone interact with D1 dopamine receptors (DRD1), the intracellular cyclic adenosine monophosphate (cAMP) level is increased. However, the intracellular level of cAMP is decreased when D2 dopamine receptors (DRD2) interact with the abovementioned signaling compounds. In our experiments, we have internalized the silica-coated silver nanoparticles (AgNP@SiO₂) in living cells to adsorb biologically generated cAMP which was probed by using SERS. Besides adsorptions of cAMP, AgNP@SiO₂ has a crucial role for the enhancement of Raman cross section of the samples. We observed the characteristic SERS peaks of cAMP when DRD1-overexpressed cells interact with the signaling compounds; these peaks were not observed for other cells including DRD2-overexpressed and DRD1-DRD2-coexpressed cells. Our experimental approach is successful to probe the intracellular cAMP and characterize the selectivity of signaling compounds to different types of DARs. Furthermore, our experimental approach is highly capable for in vivo studies because it can probe intracellular cAMP using a low input power of incident laser without significant cell damage. Our experimental results and density functional theory calculations showed that 780 and 1503 cm^-1 are signature Raman peaks of cAMP. The SERS peak at 780 cm^-1 is associated with C-O, C-C, and C-N stretching and symmetric and asymmetric bending of two O-H bonds of cAMP, whereas the SERS peak at 1503 cm^-1 is contributed by the O₉-H₃ bending mode.

Mapping the Dynamic Functions and Structural Features of AcrB Efflux Pump Transporter Using Accelerated Molecular Dynamics Simulations

Multidrug efflux pumps confer resistance to their bacterial hosts by pumping out a diverse range of compounds, including most antibiotics. Being more familiar with the details of functional dynamics and conformations of these types of pumps could help in discovering approaches to stop them functioning properly. Computational approaches, particularly conventional molecular dynamics simulations followed by diverse post simulation analysis, are powerful methods that help researchers by opening a new window to study phenomena that are not detectable in as much detail in vitro or in vivo as they are in silico. In this study, accelerated molecular dynamics simulations were applied to study the dynamics of AcrB efflux pump transporters in interaction with PAβN and tetracycline as an inhibitor and a substrate, respectively, to compare the differences in the dynamics and consequently the mechanism of action of the pump. The different dynamics for PAβN -bound form of AcrB compared to the TET-bound form is likely to affect the rotating mechanism typically observed for AcrB transporter. This shows the dynamics of the active AcrB transporter is different in a substrate-bound state compared to an inhibitor-bound state. This advances our knowledge and helps to unravel the mechanism of tripartite efflux pumps.

Peptidic Macrocycles - Conformational Sampling and Thermodynamic Characterization

Macrocycles are of considerable interest as highly specific drug candidates, yet they challenge standard conformer generators with their large number of rotatable bonds and conformational restrictions. Here, we present a molecular dynamics-based routine that bypasses current limitations in conformational sampling and extensively profiles the free energy landscape of peptidic macrocycles in solution. We perform accelerated molecular dynamics simulations to capture a diverse conformational ensemble. By applying an energetic cutoff, followed by geometric clustering, we demonstrate the striking robustness and efficiency of the approach in identifying highly populated conformational states of cyclic peptides. The resulting structural and thermodynamic information is benchmarked against interproton distances from NMR experiments and conformational states identified by X-ray crystallography. Using three different model systems of varying size and flexibility, we show that the method reliably reproduces experimentally determined structural ensembles and is capable of identifying key conformational states that include the bioactive conformation. Thus, the described approach is a robust method to generate conformations of peptidic macrocycles and holds promise for structure-based drug design.

Raman Spectroscopic Signature Markers of Dopamine-Human Dopamine Transporter Interaction in Living Cells

Dopamine (DA) controls many psychological and behavioral activities in the central nervous system (CNS) through interactions with the human dopamine transporter (hDAT) and dopamine receptors. The roles of DA in the function of the CNS are affected by the targeted binding of drugs to hDAT; thus, hDAT plays a critical role in neurophysiology and neuropathophysiology. An effective experimental method is necessary to study the DA-hDAT interaction and effects of variety of drugs like psychostimulants and antidepressants that are dependent on this interaction. In searching for obtaining and identifying the Raman spectral signatures, we have used surface enhanced Raman scattering (SERS) spectroscopy to record SERS spectra from DA, human embryonic kidney 293 cells (HEK293), hDAT-HEK293, DA-HEK293, and DA-hDAT-HEK293. We have demonstrated a specific 2D-distribution SERS spectral analytical approach to analyze DA-hDAT interaction. Our study shows that the Raman modes at 807, 839, 1076, 1090, 1538, and 1665 cm^-1 are related to DA-hDAT interaction, where Raman shifts at 807 and 1076 cm^-1 are the signature markers for the bound state of DA to probe DA-hDAT interaction. On the basis of density function theory (DFT) calculation, Raman shift of the bound state of DA at 807 cm^-1 is related to combination of bending modes α(C3-O10-H21), α(C2-O11-H22), α(C7-C8-H18), α(C6-C4-H13), α(C7-C8-H19), and α(C7-C8-N9), and Raman shift at 1076 cm^-1 is related to combination of bending modes α(H19-N9-C8), γ(N9-H19), γ(C8-H19), γ(N9-H20), γ(C8-H18), and α(C7-C8-H18). These findings demonstrate that protein-ligand interactions can be confirmed by probing change in Raman shift of ligand molecules, which could be crucial to understanding molecular interactions between neurotransmitters and their receptors or transporters.

Allosteric modulation of human dopamine transporter activity under conditions promoting its dimerization

The human dopamine (DA) transporter (hDAT) is a key regulator of neurotransmission and a target for antidepressants and addictive drugs. Despite the recent resolution of dDAT structures from Drosophila melanogaster, complete understanding of its mechanism of function and even information on its biological assembly is lacking. The resolved dDAT structures are monomeric, but growing evidence suggests that hDAT might function as a multimer, and its oligomerization may be relevant to addictive drug effects. Here, using structure-based computations, we examined the possible mechanisms of hDAT dimerization and its dynamics in a lipid bilayer. Using a combination of site-directed mutagenesis, DA-uptake, and cross-linking experiments that exploited the capacity of Cys-306 to form intermonomeric disulfide bridges in the presence of an oxidizing agent, we tested the effects of mutations at transmembrane segment (TM) 6 and 12 helices in HEK293 cells. The most probable structural model for hDAT dimer suggested by computations and experiments differed from the dimeric structure resolved for the bacterial homolog, LeuT, presumably because of a kink at TM12 preventing favorable monomer packing. Instead, TM2, TM6, and TM11 line the dimer interface. Molecular dynamics simulations of the dimeric hDAT indicated that the two subunits tend to undergo cooperative structural changes, both on local (extracellular gate opening/closure) and global (transition between outward-facing and inward-facing states) scales. These observations suggest that hDAT transport properties may be allosterically modulated under conditions promoting dimerization. Our study provides critical insights into approaches for examining the oligomerization of neurotransmitter transporters and sheds light on their drug modulation.

Ligand Binding in the Extracellular Vestibule of the Neurotransmitter Transporter Homologue LeuT

The human monoamine transporters (MATs) facilitate the reuptake of monoamine neurotransmitters from the synaptic cleft. MATs are linked to a number of neurological diseases and are the targets of both therapeutic and illicit drugs. Until recently, no high-resolution structures of the human MATs existed, and therefore, studies of this transporter family have relied on investigations of the homologues bacterial transporters such as the leucine transporter LeuT, which has been crystallized in several conformational states. A two-substrate transport mechanism has been suggested for this transporter family, which entails that high-affinity binding of a second substrate in an extracellular site is necessary for the substrate in the central binding site to be transported. Compelling evidence for this mechanism has been presented, however, a number of equally compelling accounts suggest that the transporters function through a mechanism involving only a single substrate and a single high-affinity site. To shed light on this apparent contradiction, we have performed extensive molecular dynamics simulations of LeuT in the outward-occluded conformation with either one or two substrates bound to the transporter. We have also calculated the substrate binding affinity in each of the two proposed binding sites through rigorous free energy simulations. Results show that substrate binding is unstable in the extracellular vestibule and the substrate binding affinity within the suggested extracellular site is very low (0.2 and 3.3 M for the two dominant binding modes) compared to the central substrate binding site (14 nM). This suggests that for LeuT in the outward-occluded conformation only a single high-affinity substrate binding site exists.

Effect of Dimerization on the Dynamics of Neurotransmitter:Sodium Symporters

Dimerization is a common feature among the members of the neurotransmitter:sodium symporter (NSS) family of membrane proteins. Yet, the effect of dimerization on the mechanism of action of NSS members is not fully understood. In this study, we examined the collective dynamics of two members of the family, leucine transporter (LeuT) and dopamine transporter (DAT), to assess the significance of dimerization in modulating the functional motions of the monomers. We used to this aim the anisotropic network model (ANM), an efficient and robust method for modeling the intrinsic motions of proteins and their complexes. Transporters belonging to the NSS family are known to alternate between outward-facing (OF) and inward-facing (IF) states, which enables the uptake and release of their substrate (neurotransmitter) respectively, as the substrate is transported from the exterior to the interior of the cell. In both LeuT and DAT, dimerization is found to alter the collective motions intrinsically accessible to the individual monomers in favor of the functional transitions (OF ↔ IF), suggesting that dimerization may play a role in facilitating transport.

The Environment Shapes the Inner Vestibule of LeuT

Human neurotransmitter transporters are found in the nervous system terminating synaptic signals by rapid removal of neurotransmitter molecules from the synaptic cleft. The homologous transporter LeuT, found in Aquifex aeolicus, was crystallized in different conformations. Here, we investigated the inward-open state of LeuT. We compared LeuT in membranes and micelles using molecular dynamics simulations and lanthanide-based resonance energy transfer (LRET). Simulations of micelle-solubilized LeuT revealed a stable and widely open inward-facing conformation. However, this conformation was unstable in a membrane environment. The helix dipole and the charged amino acid of the first transmembrane helix (TM1A) partitioned out of the hydrophobic membrane core. Free energy calculations showed that movement of TM1A by 0.30 nm was driven by a free energy difference of ~15 kJ/mol. Distance measurements by LRET showed TM1A movements, consistent with the simulations, confirming a substantially different inward-open conformation in lipid bilayer from that inferred from the crystal structure.

Crystal structures of the bacterial small amino acid transporter LeuT provided structural evidence for the alternating access model. Thereby, these structures shaped our understanding of the mechanisms underlying substrate translocation by neurotransmitter transporters. However, it has been questioned, if the crystallized inward-open conformation of LeuT can exist in the membrane environment. Here we show that, while stable in detergent micelles, the inward-open conformation of LeuT is of high energy and undergoes structural readjustments. We use a multi-faceted approach including molecular dynamics simulations, scintillation proximity assays, free energy calculations and apply for the first time lanthanide resonance energy transfer measurements to verify the in silico predictions. In silico and in vitro approaches using the same conditions allowed us to combine the macroscopic experimental data with microscopic all atom results from simulations to identify the underlying driving forces: partitioning of charged and polar groups from the hydrophobic membrane interior to the hydrophilic environment. We propose that the inward-facing state shows a much smaller movement of TM1A, but large enough to create an access path to the S1 substrate binding site from the vestibule.

Localization of Millisecond Dynamics: Dihedral Entropy from Accelerated MD

Here, we demonstrate a method to capture local dynamics on a time scale 3 orders of magnitude beyond state-of-the-art simulation approaches. We apply accelerated molecular dynamics simulations for conformational sampling and extract reweighted backbone dihedral distributions. Local dynamics are characterized by torsional probabilities, resulting in residue-wise dihedral entropies. Our approach is successfully validated for three different protein systems of increasing size: alanine dipeptide, bovine pancreatic trypsin inhibitor (BPTI), and the major birch pollen allergen Bet v 1a. We demonstrate excellent agreement of flexibility profiles with both large-scale computer simulations and NMR experiments. Thus, our method provides efficient access to local protein dynamics on extended time scales of high biological relevance.

Structure and Dynamics Study of LeuT Using the Markov State Model and Perturbation Response Scanning Reveals Distinct Ion Induced Conformational States

The bacterial leucine transporter (LeuT), a close homologue of the eukaryote monoamine transporters (MATs), currently serves as a powerful template for computer simulations of MATs. Transport of the amino acid leucine through the membrane is made possible by the sodium electrochemical potential. Recent reports indicate that the substrate transport mechanism is based on structural changes such as hinge movements of key transmembrane domains. In order to further investigate the role of sodium ions in the uptake of leucine, here we present a Markov state model analysis of atomistic simulations of lipid embedded LeuT in different environments, generated by varying the presence of binding pocket sodium ions and substrate. Six metastable conformations are found, and structural differences between them along with transition probabilities are determined. We complete the analysis with the implementation of perturbation response scanning on our system, determining the most sensitive and influential regions of LeuT, in each environment. Our results show that the occupation of sites Na1 and Na2, along with the presence of the substrate, selectively influences the geometry of LeuT. In particular, the occupation of each site Na1/Na2 has strong effects (in terms of changes in influence and/or sensitivity, as compared to the case without ions) in specific regions of LeuT, and the effects are different for simultaneous occupation. Our results strengthen the rationale and provide a conformational mechanism for a putative transport mechanism in which Na2 is necessary, but may not be sufficient, to initiate and stabilize extracellular substrate access to the binding pocket.

Properties of an inward-facing state of LeuT: conformational stability and substrate release

The leucine transporter (LeuT) is a bacterial homolog of the human monoamine transporters, which are important pharmaceutical targets. There are no high-resolution structures of the human transporters available; however, LeuT has been crystallized in several different conformational states. Recently, an inward-facing conformation of LeuT was solved revealing an unexpectedly large movement of transmembrane helix 1a (TM1a). We have performed molecular dynamics simulations of the mutated and wild-type transporter, with and without the cocrystallized Fab antibody fragment, to investigate the properties of this inward-facing conformation in relation to transport by LeuT within the membrane environment. In all of the simulations, local conformational changes with respect to the crystal structure are consistently observed, especially in TM1a. Umbrella sampling revealed a soft potential for TM1a tilting. Furthermore, simulations of inward-facing LeuT with Na(+) ions and substrate bound suggest that one of the Na(+) ion binding sites is fully disrupted. Release of alanine and the second Na(+) ion is also observed, giving insight into the final stage of the translocation process in atomistic detail.

Molecular Mechanism of Dopamine Transport by Human Dopamine Transporter

Dopamine transporters (DATs) control neurotransmitter dopamine (DA) homeostasis by reuptake of excess DA, assisted by sodium and chloride ions. The recent resolution of DAT structure (dDAT) from Drosophila permits us for the first time to directly view the sequence of events involved in DA reuptake in human DAT (hDAT) using homology modeling and full-atomic microseconds accelerated simulations. Major observations are spontaneous closure of extracellular gates prompted by DA binding; stabilization of a holo-occluded intermediate; disruption of N82-N353 hydrogen bond and exposure to intracellular (IC) water triggered by Na2 dislocation; redistribution of a network of salt bridges at the IC surface in the inward-facing state; concerted tilting of IC-exposed helices to enable the release of Na(+) and Cl(-) ions; and DA release after protonation of D79. The observed time-resolved interactions confirm the conserved dynamics of LeuT-fold family, while providing insights into the mechanistic role of specific residues in hDAT.

Insights into the Modulation of Dopamine Transporter Function by Amphetamine, Orphenadrine, and Cocaine Binding

Human dopamine (DA) transporter (hDAT) regulates dopaminergic signaling in the central nervous system by maintaining the synaptic concentration of DA at physiological levels, upon reuptake of DA into presynaptic terminals. DA translocation involves the co-transport of two sodium ions and the channeling of a chloride ion, and it is achieved via alternating access between outward-facing (OF) and inward-facing states of DAT. hDAT is a target for addictive drugs, such as cocaine, amphetamine (AMPH), and therapeutic antidepressants. Our recent quantitative systems pharmacology study suggested that orphenadrine (ORPH), an anticholinergic agent and anti-Parkinson drug, might be repurposable as a DAT drug. Previous studies have shown that DAT-substrates like AMPH or -blockers like cocaine modulate the function of DAT in different ways. However, the molecular mechanisms of modulation remained elusive due to the lack of structural data on DAT. The newly resolved DAT structure from Drosophila melanogaster opens the way to a deeper understanding of the mechanism and time evolution of DAT–drug/ligand interactions. Using a combination of homology modeling, docking analysis, molecular dynamics simulations, and molecular biology experiments, we performed a comparative study of the binding properties of DA, AMPH, ORPH, and cocaine and their modulation of hDAT function. Simulations demonstrate that binding DA or AMPH drives a structural transition toward a functional form predisposed to translocate the ligand. In contrast, ORPH appears to inhibit DAT function by arresting it in the OF open conformation. The analysis shows that cocaine and ORPH competitively bind DAT, with the binding pose and affinity dependent on the conformational state of DAT. Further assays show that the effect of ORPH on DAT uptake and endocytosis is comparable to that of cocaine.

Accelerated molecular dynamics and protein conformational change: a theoretical and practical guide using a membrane embedded model neurotransmitter transporter

Molecular dynamics simulation provides a powerful and accurate method to model protein conformational change, yet timescale limitations often prevent direct assessment of the kinetic properties of interest. A large number of molecular dynamic steps are necessary for rare events to occur, which allow a system to overcome energy barriers and conformationally transition from one potential energy minimum to another. For many proteins, the energy landscape is further complicated by a multitude of potential energy wells, each separated by high free-energy barriers and each potentially representative of a functionally important protein conformation. To overcome these obstacles, accelerated molecular dynamics utilizes a robust bias potential function to simulate the transition between different potential energy minima. This straightforward approach more efficiently samples conformational space in comparison to classical molecular dynamics simulation, does not require advanced knowledge of the potential energy landscape and converges to the proper canonical distribution. Here, we review the theory behind accelerated molecular dynamics and discuss the approach in the context of modeling protein conformational change. As a practical example, we provide a detailed, step-by-step explanation of how to perform an accelerated molecular dynamics simulation using a model neurotransmitter transporter embedded in a lipid cell membrane. Changes in protein conformation of relevance to the substrate transport cycle are then examined using principle component analysis.

Complete mapping of substrate translocation highlights the role of LeuT N-terminal segment in regulating transport cycle

Neurotransmitter: sodium symporters (NSSs) regulate neuronal signal transmission by clearing excess neurotransmitters from the synapse, assisted by the co-transport of sodium ions. Extensive structural data have been collected in recent years for several members of the NSS family, which opened the way to structure-based studies for a mechanistic understanding of substrate transport. Leucine transporter (LeuT), a bacterial orthologue, has been broadly adopted as a prototype in these studies. This goal has been elusive, however, due to the complex interplay of global and local events as well as missing structural data on LeuT N-terminal segment. We provide here for the first time a comprehensive description of the molecular events leading to substrate/Na+ release to the postsynaptic cell, including the structure and dynamics of the N-terminal segment using a combination of molecular simulations. Substrate and Na+-release follows an influx of water molecules into the substrate/Na+-binding pocket accompanied by concerted rearrangements of transmembrane helices. A redistribution of salt bridges and cation-π interactions at the N-terminal segment prompts substrate release. Significantly, substrate release is followed by the closure of the intracellular gate and a global reconfiguration back to outward-facing state to resume the transport cycle. Two minimally hydrated intermediates, not structurally resolved to date, are identified: one, substrate-bound, stabilized during the passage from outward- to inward-facing state (holo-occluded), and another, substrate-free, along the reverse transition (apo-occluded).

Structural dynamics of the monoamine transporter homolog LeuT from accelerated conformational sampling and channel analysis

The bacterial leucine transporter LeuT retains significant secondary structure similarities to the human monoamine transporters (MAT) such as the dopamine and serotonin reuptake proteins. The primary method of computational study of the MATs has been through the use of the crystallized LeuT structure. Different conformations of LeuT can give insight into mechanistic details of the MAT family. A conformational sampling performed through accelerated molecular dynamics simulations testing different combinations of the leucine substrate and bound sodium ions revealed seven distinct conformational clusters. Further analysis has been performed to target salt-bridge residues R30-D404, Y108-F253, and R5-D369 and transmembrane domains on both the seven isolated structures and the total trajectories. In addition, solvent accessibility of LeuT and its substrate binding pockets has been analyzed using a program for calculating channel radii. Occupation of the Na2 site stabilizes the outward conformation and should bind to the open outward conformation before the leucine and Na1 sodium while two possible pathways were found to be available for intracellular transport.

Coupled global and local changes direct substrate translocation by neurotransmitter-sodium symporter ortholog LeuT

Significant advances have been made in recent years in characterizing neurotransmitter:sodium symporter (NSS) family structure and function. Yet, many time-resolved events and intermediates that control the various stages of transport cycle remain to be elucidated. Whether NSSs harbor one or two sites for binding their substrates (neurotransmitters or amino acids), and what the role of the secondary site S2 is, if any, are still unresolved. Using molecular modeling and simulations for LeuT, a bacterial NSS, we present a comprehensive account of substrate-binding and -stabilization events, and subsequently triggered interactions leading to substrate (alanine) release. LeuT instantaneous conformation as it reconfigures from substrate-receiving (outward-facing) to -releasing (inward-facing) state appears to be a determinant of its affinity to bind substrate at site S2. In the outward-facing state, S1 robustly binds alanine and regulates subsequent redistribution of interactions to trigger extracellular gate closure; whereas S2 is only a transient binding site. The substrate-binding affinity at S2 increases in an intermediate close to inward-facing state. LeuT harbors the two substrate-binding sites, and small displacements of second substrate near S2 are observed to induce concerted small translocations in the substrate bound to primary site S1, although complete release requires collective structural rearrangements that fully expose the intracellular vestibule to the cytoplasm.

Insights from molecular dynamics: the binding site of cocaine in the dopamine transporter and permeation pathways of substrates in the leucine and dopamine transporters

The dopamine transporter (DAT) facilitates the regulation of synaptic neurotransmitter levels. As a target for therapeutic and illicit psycho-stimulant drugs like antidepressants and cocaine, DAT has been studied intensively. Despite a wealth of mutational and physiological data regarding DAT, the structure remains unsolved and details of the transport mechanism, binding sites and conformational changes remain debated. A bacterial homolog of DAT, the leucine transporter (LeuT(Aa)) has been used as a template and framework for modeling and understanding DAT. Free energy profiles obtained from Multi-Configuration Thermodynamic Integration simulations allowed us to correctly identify the primary and secondary binding pockets of LeuT(Aa). A comparison of free energy profiles for dopamine and cocaine in DAT suggests that the binding site of cocaine is located in a secondary pocket, not the primary substrate site. Two recurring primary pathways for intracellular substrate release from the primary pocket are identified in both transporters using the Random Acceleration Molecular Dynamics method. One pathway appears to follow transmembranes (TMs) 1a and 6b while the other pathway follows along TMs 6b and 8. Interestingly, we observe that a single sodium ion is co-transported with leucine during both simulation types.