Recent advances of multifunctional zwitterionic polymers for biomedical application

Zwitterionic polymers possess equal total positive and negative charges in the repeating units, making them electrically neutral overall. This unique property results in superhydrophilicity, which makes the zwitterionic polymers highly effective in resisting protein adsorption, thus endowing the drug carriers with long blood circulation time, inhibiting thrombus formation on biomedical devices in contact with blood, and ensuring the good sensitivity of sensors in biomedical application. Moreover, zwitterionic polymers have tumor-targeting ability and pH-responsiveness, rendering them ideal candidates for antitumor drug delivery. Additionally, the high ionic conductivity of zwitterionic polymers makes them an important raw material for ionic skin. Zwitterionic polymers exhibit remarkable resistance to bacterial adsorption and growth, proving their suitability in a wide range of biomedical applications such as ophthalmic applications, and wound dressings. In this paper, we provide an in-depth analysis of the different structures and characteristics of zwitterionic polymers and highlight their unique qualities and suitability for biomedical applications. Furthermore, we discuss the limitations and challenges that must be overcome to realize the full potential of zwitterionic polymers and present an optimistic perspective for zwitterionic polymers in the biomedical fields. STATEMENT OF SIGNIFICANCE: Zwitterionic polymers have a series of excellent properties such as super hydrophilicity, anti-protein adsorption, antibacterial ability and good ionic conductivity. However, biomedical applications of multifunctional zwitterionic polymers are still a major field to be explored. This review focuses on the design and application of zwitterionic polymers-based nanosystems for targeted and responsive delivery of antitumor drugs and cancer diagnostic agents. Moreover, the use of zwitterionic polymers in various biomedical applications such as biomedical devices in contact with blood, biosensors, ionic skin, ophthalmic applications and wound dressings is comprehensively described. We discuss current results and future challenges for a better understanding of multifunctional zwitterionic polymers for biomedical applications.

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Glycine Transporter 2: Mechanism and Allosteric Modulation

Neurotransmitter sodium symporters (NSS) are a subfamily of SLC6 transporters responsible for regulating neurotransmitter signalling. They are a major target for psychoactive substances including antidepressants and drugs of abuse, prompting substantial research into their modulation and structure-function dynamics. Recently, a series of allosteric transport inhibitors have been identified, which may reduce side effect profiles, compared to orthosteric inhibitors. Allosteric inhibitors are also likely to provide different clearance kinetics compared to competitive inhibitors and potentially better clinical outcomes. Crystal structures and homology models have identified several allosteric modulatory sites on NSS including the vestibule allosteric site (VAS), lipid allosteric site (LAS) and cholesterol binding site (CHOL1). Whilst the architecture of eukaryotic NSS is generally well conserved there are differences in regions that form the VAS, LAS, and CHOL1. Here, we describe ligand-protein interactions that stabilize binding in each allosteric site and explore how differences between transporters could be exploited to generate NSS specific compounds with an emphasis on GlyT2 modulation.

Structural characterisation reveals insights into substrate recognition by the glutamine transporter ASCT2/SLC1A5

Cancer cells undergo a shift in metabolism where they become reliant on nutrients such as the amino-acid glutamine. Glutamine enters the cell via the alanine/serine/cysteine transporter 2 (ASCT2) that is upregulated in several cancers to maintain an increased supply of this nutrient and are therefore an attractive target in cancer therapeutic development. ASCT2 belongs to the glutamate transporter (SLC1A) family but is the only transporter in this family able to transport glutamine. The structural basis for glutamine selectivity of ASCT2 is unknown. Here, we identify two amino-acid residues in the substrate-binding site that are responsible for conferring glutamine selectivity. We introduce corresponding mutations into a prokaryotic homologue of ASCT2 and solve four crystal structures, which reveal the structural basis for neutral amino acid and inhibitor binding in this family. This structural model of ASCT2 may provide a basis for future development of selective ASCT2 inhibitors to treat glutamine-dependent cancers.

Cancer cells are reliant on nutrients such as glutamine, which enter the cell via the alanine/serine/cysteine transporter 2 (ASCT2). Here, authors use crystallography to show which amino-acid residues in the substrate-binding site are responsible for conferring glutamine selectivity to ASCT2.

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.

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.

Effect of External Electric Field on Substrate Transport of a Secondary Active Transporter

Substrate transport across a membrane accomplished by a secondary active transporter (SAT) is essential to the normal physiological function of living cells. In the present research, a series of all-atom molecular dynamics (MD) simulations under different electric field (EF) strengths was performed to investigate the effect of an external EF on the substrate transport of an SAT. The results show that EF both affects the interaction between substrate and related protein's residues by changing their conformations and tunes the timeline of the transport event, which collectively reduces the height of energy barrier for substrate transport and results in the appearance of two intermediate conformations under the existence of an external EF. Our work spotlights the crucial influence of external EFs on the substrate transport of SATs and could provide a more penetrating understanding of the substrate transport mechanism of SATs.

Evaluations of the Absolute and Relative Free Energies for Antidepressant Binding to the Amino Acid Membrane Transporter LeuT with Free Energy Simulations

The binding of ligands to protein receptors with high affinity and specificity is central to many cellular processes. The quest for the development of computational models capable of accurately evaluating binding affinity remains one of the main goals of modern computational biophysics. In this work, free energy perturbation/molecular dynamics simulations were used to evaluate absolute and relative binding affinity for three different antidepressants to a sodium-dependent membrane transporter, LeuT, a bacterial homologue of human serotonin and dopamine transporters. Dysfunction of these membrane transporters in mammals has been implicated in multiple diseases of the nervous system, including bipolar disorder and depression. Furthermore, these proteins are key targets for antidepressants including fluoxetine (aka Prozac) and tricyclic antidepressants known to block transport activity. In addition to being clinically relevant, this system, where multiple crystal structures are readily available, represents an ideal testing ground for methods used to study the molecular mechanisms of ligand binding to membrane proteins. We discuss possible pitfalls and different levels of approximation required to evaluate binding affinity, such as the dependence of the computed affinities on the strength of constraints and the sensitivity of the computed affinities to the particular partial charges derived from restrained electrostatic potential fitting of quantum mechanics electrostatic potential. Finally, we compare the effects of different constraint schemes on the absolute and relative binding affinities obtained from free energy simulations.

Conformational Changes in Two Inter-Helical Loops of Mhp1 Membrane Transporter

Mhp1 is a bacterial secondary transporter with high-resolution crystal structures available for both the outward- and inward-facing conformations. Through molecular dynamics simulations of the ligand-free Mhp1 as well as analysis of its crystal structures, here we show that two inter-helical loops, respectively located at the extra- and intracellular ends of the "hash motif" in the protein, play important roles in the conformational transition. In the outward- and inward-facing states of the protein, the loops adopt different secondary structures, either wrapped to the end of an alpha-helix, or unwrapped to extended conformations. In equilibrium simulations of 100 ns with Mhp1 in explicit lipids and water, the loop conformations remain largely stable. In targeted molecular dynamics simulations with the protein structure driven from one state to the other, the loops exhibit resistance and only undergo abrupt changes when other parts of the protein already approach the target conformation. Free energy calculations on the isolated loops further confirm that the wrapping/unwrapping transitions are associated with substantial energetic barriers, and consist of multiple sequential steps involving the rotation of certain backbone torsion angles. Furthermore, in simulations with the loops driven from one state to the other, a large part of the protein follows the loops to the target conformation. Taken together, our simulations suggest that changes of the loop secondary structures would be among the slow degrees of freedom in the conformational transition of the entire protein. Incorporation of detailed loop structures into the reaction coordinate, therefore, should improve the convergence and relevance of the resulting conformational free energy.

Simulating the distance distribution between spin-labels attached to proteins

EPR/DEER spectroscopy is playing an increasingly important role in the characterization of the conformational states of proteins. In this study, force field parameters for the bifunctional spin-label (RX) used in EPR/DEER are parametrized and tested with molecular dynamics (MD) simulations. The dihedral angles connecting the C_α atom of the backbone to the nitroxide ring moiety of the RX spin-label attached to i and i + 4 positions in a polyalanine α-helix agree very well with those observed in the X-ray crystallography. Both RX_i,i+4 and RX_i,i+3 are more rigid than the monofunctional spin-label (R1) commonly used in EPR/DEER, while RX_i,i+4 is more rigid and causes less distortion in a protein backbone than RX_i,i+3. Simplified dummy spin-label models with a single effective particle representing the RX_i,i+3 and RX_i,i+4 are also developed and parametrized from the all-atom simulations. MD simulations with dummy spin-labels (MDDS) provide distance distributions that can be directly compared to distance distributions obtained from EPR/DEER to rapidly assess if a hypothetical three-dimensional (3D) structural model is consistent with experiment. The dummy spin-labels can also be used in the restrained-ensemble MD (re-MD) simulations to carry out structural refinement of 3D models. Applications of this methodology to T4 lysozyme, KCNE1, and LeuT are shown to provide important insights about their conformational dynamics.

A consistent force field parameter set for zwitterionic amino acid residues

Isolated amino acids play an important role in biochemistry and are therefore an interesting object of study. Atomistic molecular dynamics (MD) simulations can provide a high-resolution picture of the dynamic features of these species, especially in their biological environment. Unfortunately, most standard force field packages lack libraries for isolated amino acids in their zwitterionic form. Although several studies have used ad-hoc parameterizations for single amino acids, a consistent force-field parameter set for these molecules is still missing. Here, we present such a parameter library derived from the widely used parm99SB set from the AMBER program package. The parameter derivation for all 20 proteinogenic amino acids transparently followed established procedures with histidine treated in three different protonation states. All amino acids were subjected to MD simulations in four different forms for comparison: zwitterionic, N-teminally capped with acetyl, C-terminally capped with N-methyl, and capped at both termini. Simulation results show similarities between the different forms. Five zwitterionic amino acids-arginine, glutamate, glycine, phenylalanine, leucine-were simulated in a protein environment. Proteins and ligands generally retained their initial structure. The new parameter set will thus facilitate future atomistic simulations of these species.

The use of LeuT as a model in elucidating binding sites for substrates and inhibitors in neurotransmitter transporters

BACKGROUND

The mammalian neurotransmitter transporters are complex proteins playing a central role in synaptic transmission between neurons by rapid reuptake of neurotransmitters. The proteins which transport dopamine, noradrenaline and serotonin belong to the Neurotransmitter:Sodium Symporters (NSS). Due to their important role, dysfunctions are associated with several psychiatric and neurological diseases and they also serve as targets for a wide range of therapeutic and illicit drugs. Despite the central physiological and pharmacological importance, direct evidence on structure-function relationships on mammalian NSS proteins has so far been unsuccessful. The crystal structure of the bacterial NSS protein, LeuT, has been a turning point in structural investigations.

SCOPE OF REVIEW

To provide an update on what is known about the binding sites for substrates and inhibitors in the LeuT. The different binding modes and binding sites will be discussed with special emphasis on the possible existence of a second substrate binding site. It is the goal to give an insight into how investigations on ligand binding in LeuT have provided basic knowledge about transporter conformations and translocation mechanism which can pave the road for a deeper understanding of drug binding and function of the mammalian transporters.

MAJOR CONCLUSIONS

The LeuT is a suitable model for the structural investigation of NSS proteins including the possible location of drug binding sites. It is still debated whether the LeuT is a suitable model for the molecular mechanisms behind substrate translocation.

GENERAL SIGNIFICANCE

Structure and functional aspects of NSS proteins are central for understanding synaptic transmission. With the purification and crystallization of LeuT as well as the dopamine transporter from Drosophila melanogaster, the application of biophysical methods such as fluorescence spectroscopy, neutron- or x-ray scattering and NMR for understanding its function becomes increasingly available. This article is part of a Special Issue entitled Structural biochemistry and biophysics of membrane proteins.

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.

Chloride binding site of neurotransmitter sodium symporters

Neurotransmitter:sodium symporters (NSSs) play a critical role in signaling by reuptake of neurotransmitters. Eukaryotic NSSs are chloride-dependent, whereas prokaryotic NSS homologs like LeuT are chloride-independent but contain an acidic residue (Glu290 in LeuT) at a site where eukaryotic NSSs have a serine. The LeuT-E290S mutant displays chloride-dependent activity. We show that, in LeuT-E290S cocrystallized with bromide or chloride, the anion is coordinated by side chain hydroxyls from Tyr47, Ser290, and Thr254 and the side chain amide of Gln250. The bound anion and the nearby sodium ion in the Na1 site organize a connection between their coordinating residues and the extracellular gate of LeuT through a continuous H-bond network. The specific insights from the structures, combined with results from substrate binding studies and molecular dynamics simulations, reveal an anion-dependent occlusion mechanism for NSS and shed light on the functional role of chloride binding.

Visualizing functional motions of membrane transporters with molecular dynamics simulations

Computational modeling and molecular simulation techniques have become an integral part of modern molecular research. Various areas of molecular sciences continue to benefit from, indeed rely on, the unparalleled spatial and temporal resolutions offered by these technologies, to provide a more complete picture of the molecular problems at hand. Because of the continuous development of more efficient algorithms harvesting ever-expanding computational resources, and the emergence of more advanced and novel theories and methodologies, the scope of computational studies has expanded significantly over the past decade, now including much larger molecular systems and far more complex molecular phenomena. Among the various computer modeling techniques, the application of molecular dynamics (MD) simulation and related techniques has particularly drawn attention in biomolecular research, because of the ability of the method to describe the dynamical nature of the molecular systems and thereby to provide a more realistic representation, which is often needed for understanding fundamental molecular properties. The method has proven to be remarkably successful in capturing molecular events and structural transitions highly relevant to the function and/or physicochemical properties of biomolecular systems. Herein, after a brief introduction to the method of MD, we use a number of membrane transport proteins studied in our laboratory as examples to showcase the scope and applicability of the method and its power in characterizing molecular motions of various magnitudes and time scales that are involved in the function of this important class of membrane proteins.

Identification of the substrate recognition and transport pathway in a eukaryotic member of the nucleobase-ascorbate transporter (NAT) family

Using the crystal structure of the uracil transporter UraA of Escherichia coli, we constructed a 3D model of the Aspergillus nidulans uric acid-xanthine/H(+) symporter UapA, which is a prototype member of the Nucleobase-Ascorbate Transporter (NAT) family. The model consists of 14 transmembrane segments (TMSs) divided into a core and a gate domain, the later being distinctly different from that of UraA. By implementing Molecular Mechanics (MM) simulations and quantitative structure-activity relationship (SAR) approaches, we propose a model for the xanthine-UapA complex where the substrate binding site is formed by the polar side chains of residues E356 (TMS8) and Q408 (TMS10) and the backbones of A407 (TMS10) and F155 (TMS3). In addition, our model shows several polar interactions between TMS1-TMS10, TMS1-TMS3, TMS8-TMS10, which seem critical for UapA transport activity. Using extensive docking calculations we identify a cytoplasm-facing substrate trajectory (D360, A363, G411, T416, R417, V463 and A469) connecting the proposed substrate binding site with the cytoplasm, as well as, a possible outward-facing gate leading towards the substrate major binding site. Most importantly, re-evaluation of the plethora of available and analysis of a number of herein constructed UapA mutations strongly supports the UapA structural model. Furthermore, modeling and docking approaches with mammalian NAT homologues provided a molecular rationale on how specificity in this family of carriers might be determined, and further support the importance of selectivity gates acting independently from the major central substrate binding site.

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.

The structural pathway for water permeation through sodium-glucose cotransporters

Although water permeation across cell membranes occurs through several types of membrane proteins, the only permeation mechanism resolved at atomic scale is that through aquaporins. Crystallization of the Vibrio parahaemolyticus sodium-galactose transporter (vSGLT) allows investigation of putative water permeation pathways through both vSGLT and the homologous human Na-glucose cotransporter (hSGLT1) using computational methods. Grand canonical Monte Carlo and molecular dynamics simulations were used to stably insert water molecules in both proteins, showing the presence of a water-filled pathway composed of ∼100 water molecules. This provides a structural basis for passive water permeation that is difficult to reconcile with the water cotransport hypothesis. Potential-of-mean-force calculations of water going through the crystal structure of vSGLT shows a single barrier of 7.7 kCal mol(-1), in agreement with previously published experimental data for cotransporters of the SGLT family. Electrophysiological and volumetric experiments performed on hSGLT1-expressing Xenopus oocytes showed that the passive permeation pathway exists in different conformational states. In particular, experimental conditions that aim to mimic the conformation of the crystal structure displayed passive water permeability. These results provide groundwork for understanding the structural basis of cotransporter water permeability.

Atomistic models of ion and solute transport by the sodium-dependent secondary active transporters

The recent determination of high-resolution crystal structures of several transporters offers unprecedented insights into the structural mechanisms behind secondary transport. These proteins utilize the facilitated diffusion of the ions down their electrochemical gradients to transport the substrate against its concentration gradient. The structural studies revealed striking similarities in the structural organization of ion and solute binding sites and a well-conserved inverted-repeat topology between proteins from several gene families. In this paper we will overview recent atomistic simulations applied to study the mechanisms of selective binding of ion and substrate in LeuT, Glt, vSGLT and hSERT as well as its consequences for the transporter conformational dynamics. This article is part of a Special Issue entitled: Membrane protein structure and function.

SLC6 neurotransmitter transporters: structure, function, and regulation

The neurotransmitter transporters (NTTs) belonging to the solute carrier 6 (SLC6) gene family (also referred to as the neurotransmitter-sodium-symporter family or Na(+)/Cl(-)-dependent transporters) comprise a group of nine sodium- and chloride-dependent plasma membrane transporters for the monoamine neurotransmitters serotonin (5-hydroxytryptamine), dopamine, and norepinephrine, and the amino acid neurotransmitters GABA and glycine. The SLC6 NTTs are widely expressed in the mammalian brain and play an essential role in regulating neurotransmitter signaling and homeostasis by mediating uptake of released neurotransmitters from the extracellular space into neurons and glial cells. The transporters are targets for a wide range of therapeutic drugs used in treatment of psychiatric diseases, including major depression, anxiety disorders, attention deficit hyperactivity disorder and epilepsy. Furthermore, psychostimulants such as cocaine and amphetamines have the SLC6 NTTs as primary targets. Beginning with the determination of a high-resolution structure of a prokaryotic homolog of the mammalian SLC6 transporters in 2005, the understanding of the molecular structure, function, and pharmacology of these proteins has advanced rapidly. Furthermore, intensive efforts have been directed toward understanding the molecular and cellular mechanisms involved in regulation of the activity of this important class of transporters, leading to new methodological developments and important insights. This review provides an update of these advances and their implications for the current understanding of the SLC6 NTTs.

A conserved asparagine residue in transmembrane segment 1 (TM1) of serotonin transporter dictates chloride-coupled neurotransmitter transport

Na(+)- and Cl(-)-dependent uptake of neurotransmitters via transporters of the SLC6 family, including the human serotonin transporter (SLC6A4), is critical for efficient synaptic transmission. Although residues in the human serotonin transporter involved in direct Cl(-) coordination of human serotonin transport have been identified, the role of Cl(-) in the transport mechanism remains unclear. Through a combination of mutagenesis, chemical modification, substrate and charge flux measurements, and molecular modeling studies, we reveal an unexpected role for the highly conserved transmembrane segment 1 residue Asn-101 in coupling Cl(-) binding to concentrative neurotransmitter uptake.

The substrate-driven transition to an inward-facing conformation in the functional mechanism of the dopamine transporter

Background

The dopamine transporter (DAT), a member of the neurotransmitter:Na⁺ symporter (NSS) family, terminates dopaminergic neurotransmission and is a major molecular target for psychostimulants such as cocaine and amphetamine, and for the treatment of attention deficit disorder and depression. The crystal structures of the prokaryotic NSS homolog of DAT, the leucine transporter LeuT, have provided critical structural insights about the occluded and outward-facing conformations visited during the substrate transport, but only limited clues regarding mechanism. To understand the transport mechanism in DAT we have used a homology model based on the LeuT structure in a computational protocol validated previously for LeuT, in which steered molecular dynamics (SMD) simulations guide the substrate along a pathway leading from the extracellular end to the intracellular (cytoplasmic) end.

Methodology/Principal Findings

Key findings are (1) a second substrate binding site in the extracellular vestibule, and (2) models of the conformational states identified as occluded, doubly occupied, and inward-facing. The transition between these states involve a spatially ordered sequence of interactions between the two substrate-binding sites, followed by rearrangements in structural elements located between the primary binding site and the cytoplasmic end. These rearrangements are facilitated by identified conserved hinge regions and a reorganization of interaction networks that had been identified as gates.

Conclusions/Significance

Computational simulations supported by information available from experiments in DAT and other NSS transporters have produced a detailed mechanistic proposal for the dynamic changes associated with substrate transport in DAT. This allosteric mechanism is triggered by the binding of substrate in the S2 site in the presence of the substrate in the S1 site. Specific structural elements involved in this mechanism, and their roles in the conformational transitions illuminated here describe, a specific substrate-driven allosteric mechanism that is directly amenable to experiment as shown previously for LeuT.

Structural refinement of the hERG1 pore and voltage-sensing domains with ROSETTA-membrane and molecular dynamics simulations

The hERG1 gene (Kv11.1) encodes a voltage-gated potassium channel. Mutations in this gene lead to one form of the Long QT Syndrome (LQTS) in humans. Promiscuous binding of drugs to hERG1 is known to alter the structure/function of the channel leading to an acquired form of the LQTS. Expectably, creation and validation of reliable 3D model of the channel have been a key target in molecular cardiology and pharmacology for the last decade. Although many models were built, they all were limited to pore domain. In this work, a full model of the hERG1 channel is developed which includes all transmembrane segments. We tested a template-driven de-novo design with ROSETTA-membrane modeling using side-chain placements optimized by subsequent molecular dynamics (MD) simulations. Although backbone templates for the homology modeled parts of the pore and voltage sensors were based on the available structures of KvAP, Kv1.2 and Kv1.2-Kv2.1 chimera channels, the missing parts are modeled de-novo. The impact of several alignments on the structure of the S4 helix in the voltage-sensing domain was also tested. Herein, final models are evaluated for consistency to the reported structural elements discovered mainly on the basis of mutagenesis and electrophysiology. These structural elements include salt bridges and close contacts in the voltage-sensor domain; and the topology of the extracellular S5-pore linker compared with that established by toxin foot-printing and nuclear magnetic resonance studies. Implications of the refined hERG1 model to binding of blockers and channels activators (potent new ligands for channel activations) are discussed.

Modeling and dynamics of the inward-facing state of a Na+/Cl- dependent neurotransmitter transporter homologue

The leucine transporter (LeuT) has recently commanded exceptional attention due mainly to two distinctions; it provides the only crystal structures available for a protein homologous to the pharmacologically relevant neurotransmitter: sodium symporters (NSS), and, it exhibits a hallmark 5-TM inverted repeat ("LeuT-fold"), a fold recently discovered to also exist in several secondary transporter families, underscoring its general role in transporter function. Constructing the transport cycle of "LeuT-fold" transporters requires detailed structural and dynamic descriptions of the outward-facing (OF) and inward-facing (IF) states, as well as the intermediate states. To this end, we have modeled the structurally unknown IF state of LeuT, based on the known crystal structures of the OF state of LeuT and the IF state of vSGLT, a "LeuT-fold" transporter. The detailed methodology developed for the study combines structure-based alignment, threading, targeted MD and equilibrium MD, and can be applied to other proteins. The resulting IF-state models maintain the secondary structural features of LeuT. Water penetration and solvent accessibility calculations show that TM1, TM3, TM6 and TM8 line the substrate binding/unbinding pathway with TM10 and its pseudosymmetric partner, TM5, participating in the extracellular and intracellular halves of the lumen, respectively. We report conformational hotspots where notable changes in interactions occur between the IF and OF states. We observe Na2 exiting the LeuT-substrate- complex in the IF state, mainly due to TM1 bending. Inducing a transition in only one of the two pseudosymmetric domains, while allowing the second to respond dynamically, is found to be sufficient to induce the formation of the IF state. We also propose that TM2 and TM7 may be facilitators of TM1 and TM6 motion. Thus, this study not only presents a novel modeling methodology applied to obtain the IF state of LeuT, but also describes structural elements involved in a possibly general transport mechanism in transporters adopting the "LeuT-fold".

An efficient synthesis of a rationally designed 1,5 disubstituted imidazole AT(1) angiotensin II receptor antagonist: reorientation of imidazole pharmacophore groups in losartan reserves high receptor affinity and confirms docking studies

A new 1,5 disubstituted imidazole AT(1) Angiotensin II (AII) receptor antagonist related to losartan with reversion of butyl and hydroxymethyl groups at the 2-, 5-positions of the imidazole ring was synthesized and evaluated for its antagonist activity (V8). In vitro results indicated that the reorientation of butyl and hydroxymethyl groups on the imidazole template of losartan retained high binding affinity to the AT(1) receptor concluding that the spacing of the substituents at the 2,5- positions is of primary importance. The docking studies are confirmed by binding assay results which clearly show a comparable binding score of the designed compound V8 with that of the prototype losartan. An efficient, regioselective and cost effective synthesis renders the new compound as an attractive candidate for advanced toxicological evaluation and a drug against hypertension.

Ion/substrate-dependent conformational dynamics of a bacterial homolog of neurotransmitter:sodium symporters

Crystallographic, computational and functional analyses of LeuT have revealed details of the molecular architecture of Na(+)-coupled transporters and the mechanistic nature of ion/substrate coupling, but the conformational changes that support a functional transport cycle have yet to be described fully. We have used site-directed spin labeling and electron paramagnetic resonance (EPR) analysis to capture the dynamics of LeuT in the region of the extracellular vestibule associated with the binding of Na(+) and leucine. The results outline the Na(+)-dependent formation of a dynamic outward-facing intermediate that exposes the primary substrate binding site and the conformational changes that occlude this binding site upon subsequent binding of the leucine substrate. Furthermore, the binding of the transport inhibitors tryptophan, clomipramine and octyl-glucoside is shown to induce structural changes that distinguish the resulting inhibited conformation from the Na(+)/leucine-bound state.

Alchemical free energy simulations for biological complexes: powerful but temperamental

Free energy simulations compare multiple ligand:receptor complexes by "alchemically" transforming one into another, yielding binding free energy differences. Since their introduction in the 1980s, many technical and theoretical obstacles were surmounted, and the method ("MDFE," since molecular dynamics are often used) has matured into a powerful tool. We describe its current status, its effectiveness, and the challenges it faces. MDFE has provided chemical accuracy for many systems but remains expensive, with significant human overhead costs. The bottlenecks have shifted, partly due to increased computer power. To study diverse sets of ligands, force field availability and accuracy can be a major difficulty. Another difficulty is the frequent need to consider multiple states, related to sidechain protonation or buried waters, for example. Sophisticated, automated methods to sample these states are maturing, such as constant pH simulations. Meanwhile, combinations of MDFE and simpler approaches, like continuum dielectric models, can be very effective. As illustrations, we show how, with careful force field parameterization, MDFE accurately predicts binding specificities between complex tetracycline ligands and their targets. We describe substrate binding to the aspartyl-tRNA synthetase enzyme, where many distinct electrostatic states play a role, and a histidine and a Mg(2+) ion act as coupled switches that help enforce a strict preference for the aspartate substrate, relative to several analogs. Overall, MDFE has achieved a predictive status, where novel ligands can be studied and molecular recognition elucidated in depth. It should play an increasing role in the analysis of complex cellular processes and biomolecular engineering.

Molecular mechanism of ion-ion and ion-substrate coupling in the Na+-dependent leucine transporter LeuT

Ion-coupled transport of neurotransmitter molecules by neurotransmitter:sodium symporters (NSS) play an important role in the regulation of neuronal signaling. One of the major events in the transport cycle is ion-substrate coupling and formation of the high-affinity occluded state with bound ions and substrate. Molecular mechanisms of ion-substrate coupling and the corresponding ion-substrate stoichiometry in NSS transporters has yet to be understood. The recent determination of a high-resolution structure for a bacterial homolog of Na(+)/Cl(-)-dependent neurotransmitter transporters, LeuT, offers a unique opportunity to analyze the functional roles of the multi-ion binding sites within the binding pocket. The binding pocket of LeuT contains two metal binding sites. The first ion in site NA1 is directly coupled to the bound substrate (Leu) with the second ion in the neighboring site (NA2) only approximately 7 A away. Extensive, fully atomistic, molecular dynamics, and free energy simulations of LeuT in an explicit lipid bilayer are performed to evaluate substrate-binding affinity as a function of the ion load (single versus double occupancy) and occupancy by specific monovalent cations. It was shown that double ion occupancy of the binding pocket is required to ensure substrate coupling to Na(+) and not to Li(+) or K(+) cations. Furthermore, it was found that presence of the ion in site NA2 is required for structural stability of the binding pocket as well as amplified selectivity for Na(+) in the case of double ion occupancy.

Control of ion selectivity in LeuT: two Na+ binding sites with two different mechanisms

The x-ray structure of LeuT, a bacterial homologue of Na(+)/Cl(-)-dependent neurotransmitter transporters, provides a great opportunity to better understand the molecular basis of monovalent cation selectivity in ion-coupled transporters. LeuT possesses two ion binding sites, NA1 and NA2, which are highly selective for Na(+). Extensive all-atom free-energy molecular dynamics simulations of LeuT embedded in an explicit membrane are performed at different temperatures and various occupancy states of the binding sites to dissect the molecular mechanism of ion selectivity. The results show that the two binding sites display robust selectivity for Na(+) over K(+) or Li(+), the competing ions of most similar radii. Of particular interest, the mechanism primarily responsible for selectivity for each of the two binding sites appears to be different. In NA1, selectivity for Na(+) over K(+) arises predominantly from the strong electrostatic field arising from the negatively charged carboxylate group of the leucine substrate coordinating the ion directly. In NA2, which comprises only neutral ligands, selectivity for Na(+) is enforced by the local structural restraints arising from the hydrogen-bonding network and the covalent connectivity of the polypeptide chain surrounding the ion according to a "snug-fit" mechanism.