Heidi Koldsø

Lipid clustering correlates with membrane curvature as revealed by molecular simulations of complex lipid bilayers.

Cell membranes are complex multicomponent systems, which are highly heterogeneous in the lipid distribution and composition. To date, most molecular simulations have focussed on relatively simple lipid compositions, helping to inform our understanding of in vitro experimental studies. Here we describe on simulations of complex asymmetric plasma membrane model, which contains seven different lipids species including the glycolipid GM3 in the outer leaflet and the anionic lipid, phosphatidylinositol 4,5-bisphophate (PIP2), in the inner leaflet. Plasma membrane models consisting of 1500 lipids and resembling the in vivo composition were constructed and simulations were run for 5 µs. In these simulations the most striking feature was the formation of nano-clusters of GM3 within the outer leaflet. In simulations of protein interactions within a plasma membrane model, GM3, PIP2, and cholesterol all formed favorable interactions with the model α-helical protein. A larger scale simulation of a model plasma membrane containing 6000 lipid molecules revealed correlations between curvature of the bilayer surface and clustering of lipid molecules. In particular, the concave (when viewed from the extracellular side) regions of the bilayer surface were locally enriched in GM3. In summary, these simulations explore the nanoscale dynamics of model bilayers which mimic the in vivo lipid composition of mammalian plasma membranes, revealing emergent nanoscale membrane organization which may be coupled both to fluctuations in local membrane geometry and to interactions with proteins.

Binding and Orientation of Tricyclic Antidepressants within the Central Substrate Site of the Human Serotonin Transporter

Tricyclic antidepressants (TCAs) have been used for decades, but their orientation within and molecular interactions with their primary target is yet unsettled. The recent finding of a TCA binding site in the extracellular vestibule of the bacterial leucine transporter 11 Å above the central site has prompted debate about whether this vestibular site in the bacterial transporter is applicable to binding of antidepressants to their relevant physiological target, the human serotonin transporter (hSERT). We present an experimentally validated structural model of imipramine and analogous TCAs in the central substrate binding site of hSERT. Two possible binding modes were observed from induced fit docking calculations. We experimentally validated a single binding mode by combining mutagenesis of hSERT with uptake inhibition studies of different TCA analogs according to the paired mutation ligand analog complementation paradigm. Using this experimental method, we identify a salt bridge between the tertiary aliphatic amine and Asp98. Furthermore, the 7-position of the imipramine ring is found vicinal to Phe335, and the pocket lined by Ala173 and Thr439 is utilized by 3-substituents. These protein-ligand contact points unambiguously orient the TCA within the central binding site and reveal differences between substrate binding and inhibitor binding, giving important clues to the inhibition mechanism. Consonant with the well established competitive inhibition of uptake by TCAs, the resulting binding site for TCAs in hSERT is fully overlapping with the serotonin binding site in hSERT and dissimilar to the low affinity noncompetitive TCA site reported in the leucine transporter (LeuT).

Organization and Dynamics of Receptor Proteins in a Plasma Membrane.

The interactions of membrane proteins are influenced by their lipid environment, with key lipid species able to regulate membrane protein function. Advances in high-resolution microscopy can reveal the organization and dynamics of proteins and lipids within living cells at resolutions <200 nm. Parallel advances in molecular simulations provide near-atomic-resolution models of the dynamics of the organization of membranes of in vivo-like complexity. We explore the dynamics of proteins and lipids in crowded and complex plasma membrane models, thereby closing the gap in length and complexity between computations and experiments. Our simulations provide insights into the mutual interplay between lipids and proteins in determining mesoscale (20-100 nm) fluctuations of the bilayer, and in enabling oligomerization and clustering of membrane proteins.

Unbiased Simulations Reveal the Inward-Facing Conformation of the Human Serotonin Transporter and Na+ Ion Release

Monoamine transporters are responsible for termination of synaptic signaling and are involved in depression, control of appetite, and anxiety amongst other neurological processes. Despite extensive efforts, the structures of the monoamine transporters and the transport mechanism of ions and substrates are still largely unknown. Structural knowledge of the human serotonin transporter (hSERT) is much awaited for understanding the mechanistic details of substrate translocation and binding of antidepressants and drugs of abuse. The publication of the crystal structure of the homologous leucine transporter has resulted in homology models of the monoamine transporters. Here we present extended molecular dynamics simulations of an experimentally supported homology model of hSERT with and without the natural substrate yielding a total of more than 1.5 µs of simulation of the protein dimer. The simulations reveal a transition of hSERT from an outward-facing occluded conformation to an inward-facing conformation in a one-substrate-bound state. Simulations with a second substrate in the proposed symport effector site did not lead to conformational changes associated with translocation. The central substrate binding site becomes fully exposed to the cytoplasm leaving both the Na+-ion in the Na2-site and the substrate in direct contact with the cytoplasm through water interactions. The simulations reveal how sodium is released and show indications of early events of substrate transport. The notion that ion dissociation from the Na2-site drives translocation is supported by experimental studies of a Na2-site mutant. Transmembrane helices (TMs) 1 and 6 are identified as the helices involved in the largest movements during transport.

Molecular basis for selective serotonin reuptake inhibition by the antidepressant agent fluoxetine (Prozac).

Inhibitors of the serotonin transporter (SERT) are widely used antidepressant agents, but the structural mechanism for inhibitory activity and selectivity over the closely related norepinephrine transporter (NET) is not well understood. Here we use a combination of chemical, biological, and computational methods to decipher the molecular basis for high-affinity recognition in SERT and selectivity over NET for the prototypical antidepressant drug fluoxetine (Prozac; Eli Lilly, Indianapolis, IN). We show that fluoxetine binds within the central substrate site of human SERT, in agreement with recent X-ray crystal structures of LeuBAT, an engineered monoamine-like version of the bacterial amino acid transporter LeuT. However, the binding orientation of fluoxetine is reversed in our experimentally supported model compared with the LeuBAT structures, emphasizing the need for careful experimental verification when extrapolating findings from crystal structures of bacterial transporters to human relatives. We find that the selectivity of fluoxetine and nisoxetine, a NET selective structural congener of fluoxetine, is controlled by residues in different regions of the transporters, indicating a complex mechanism for selective recognition of structurally similar compounds in SERT and NET. Our findings add important new information on the molecular basis for SERT/NET selectivity of antidepressants, and provide the first assessment of the potential of LeuBAT as a model system for antidepressant binding in human transporters, which is essential for future structure-based drug development of antidepressant drugs with fine-tuned transporter selectivity.

Ligand Induced Conformational Changes of the Human Serotonin Transporter Revealed by Molecular Dynamics Simulations

The competitive inhibitor cocaine and the non-competitive inhibitor ibogaine induce different conformational states of the human serotonin transporter. It has been shown from accessibility experiments that cocaine mainly induces an outward-facing conformation, while the non-competitive inhibitor ibogaine, and its active metabolite noribogaine, have been proposed to induce an inward-facing conformation of the human serotonin transporter similar to what has been observed for the endogenous substrate, serotonin. The ligand induced conformational changes within the human serotonin transporter caused by these three different types of ligands, substrate, non-competitive and competitive inhibitors, are studied from multiple atomistic molecular dynamics simulations initiated from a homology model of the human serotonin transporter. The results reveal that diverse conformations of the human serotonin transporter are captured from the molecular dynamics simulations depending on the type of the ligand bound. The inward-facing conformation of the human serotonin transporter is reached with noribogaine bound, and this state resembles a previously identified inward-facing conformation of the human serotonin transporter obtained from molecular dynamics simulation with bound substrate, but also a recently published inward-facing conformation of a bacterial homolog, the leucine transporter from Aquifex Aoelicus. The differences observed in ligand induced behavior are found to originate from different interaction patterns between the ligands and the protein. Such atomic-level understanding of how an inhibitor can dictate the conformational response of a transporter by ligand binding may be of great importance for future drug design.

The juxtamembrane regions of human receptor tyrosine kinases exhibit conserved interaction sites with anionic lipids

Receptor tyrosine kinases (RTKs) play a critical role in diverse cellular processes and their activity is regulated by lipids in the surrounding membrane, including PIP2 (phosphatidylinositol-4,5-bisphosphate) in the inner leaflet, and GM3 (monosialodihexosylganglioside) in the outer leaflet. However, the precise details of the interactions at the molecular level remain to be fully characterised. Using a multiscale molecular dynamics simulation approach, we comprehensively characterise anionic lipid interactions with all 58 known human RTKs. Our results demonstrate that the juxtamembrane (JM) regions of RTKs are critical for inducing clustering of anionic lipids, including PIP2, both in simple asymmetric bilayers, and in more complex mixed membranes. Clustering is predominantly driven by interactions between a conserved cluster of basic residues within the first five positions of the JM region, and negatively charged lipid headgroups. This highlights a conserved interaction pattern shared across the human RTK family. In particular predominantly the N-terminal residues of the JM region are involved in the interactions with PIP2, whilst residues within the distal JM region exhibit comparatively less lipid specificity. Our results suggest that JM–lipid interactions play a key role in RTK structure and function, and more generally in the nanoscale organisation of receptor-containing cell membranes.

Comparative modeling of the human monoamine transporters: similarities in substrate binding.

The amino acid compositions of the substrate binding pockets of the three human monoamine transporters are compared as is the orientation of the endogenous substrates, serotonin, dopamine, and norepinephrine, bound in these. Through a combination of homology modeling, induced fit dockings, molecular dynamics simulations, and uptake experiments in mutant transporters, we propose a common binding mode for the three substrates. The longitudinal axis of the substrates is similarly oriented with these, forming an ionic interaction between the ammonium group and a highly conserved aspartate, Asp98 (serotonin transporter, hSERT), Asp79 (dopamine transporter, hDAT), and Asp75 (norepinephrine transporter, hNET). The 6-position of serotonin and the para-hydroxyl groups of dopamine and norepinephrine were found to face Ala173 in hSERT, Gly153 in hDAT, and Gly149 in hNET. Three rotations of the substrates around the longitudinal axis were identified. In each mode, an aromatic hydroxyl group of the substrates occupied equivalent volumes of the three binding pockets, where small changes in amino acid composition explains the differences in selectivity. Uptake experiments support that the 5-hydroxyl group of serotonin and the meta-hydroxyl group norepinephrine and dopamine are placed in the hydrophilic pocket around Ala173, Ser438, and Thr439 in hSERT corresponding to Gly149, Ser419, Ser420 in hNET and Gly153 Ser422 and Ala423 in hDAT. Furthermore, hDAT was found to possess an additional hydrophilic pocket around Ser149 to accommodate the para-hydroxyl group. Understanding these subtle differences between the binding site compositions of the three transporters is imperative for understanding the substrate selectivity, which could eventually aid in developing future selective medicines.

Interactions of the EGFR juxtamembrane domain with PIP2-containing lipid bilayers: Insights from multiscale molecular dynamics simulations

Background The epidermal growth factor receptor (EGFR) is the best characterised member of the receptor tyrosine kinases, which play an important role in signalling across mammalian cell membranes. The EGFR juxtamembrane (JM) domain is involved in the mechanism of activation of the receptor, interacting with the anionic lipid phosphatidylinositol 4,5-bisphosphate (PIP2) in the intracellular leaflet of the cell membrane. Methods Multiscale MD simulations were used to characterize PIP2–JM interactions. Simulations of the transmembrane helix plus JM region (TM–JM) dimer (PDB:2M20) in both PIP2-containing and PIP2-depleted lipid bilayer membranes revealed the interactions of the JM with PIP2 and other lipids. Results PIP2 forms strong interactions with the basic residues in the R645–R647 motif of the JM domain resulting in clustering of PIP2 around the protein. This association of PIP2 and the JM domain aids stabilization of JM-A dimer away from the membrane. Mutation (R645N/R646N/R647N) or PIP2-depletion results in deformation of the JM-A dimer and changes in JM–membrane interactions. Conclusions These simulations support the proposal that the positively charged residues at the start of the JM-A domain stabilize the JM-A helices in an orientation away from the membrane surface through binding to PIP2, thus promoting a conformation corresponding to an asymmetric (i.e. activated) kinase. General significance This study indicates that MD simulations may be used to characterise JM/lipid interactions, thus helping to define their role in the mechanisms of receptor tyrosine kinases. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Local Lipid Reorganization by a Transmembrane Protein Domain.

Membrane proteins interact with their lipid bilayer environment via both a transmembrane helix and juxtamembrane regions. The effect of juxtamembrane regions and membrane lipid composition on these interactions has been explored by multiscale molecular dynamics simulations. The consequences of anionic lipids within the inner leaflet of a membrane were studied in combination with membrane spanning protein models differing in their juxtamembrane domains. The simulations reveal sensitivity of the protein-lipid interactions to membrane lipid composition and charged amino acid side chains. Basic residues on the intracellular side of the protein facilitated interactions with anionic lipids. Protein systems without basic residues do not show selectivity for anionic compared with zwitterionic lipids. This reveals the sensitivity to the composition of both the membrane and the protein system when studying membrane-embedded proteins. The results presented here illustrate how even a simple transmembrane domain is able to induce lipid reorganization in a mixed asymmetric bilayer.