Identification of an allosteric binding site on the human glycine transporter, GlyT2, for bioactive lipid analgesics.

The treatment of chronic pain is poorly managed by current analgesics, and there is a need for new classes of drugs. We recently developed a series of bioactive lipids that inhibit the human glycine transporter GlyT2 (SLC6A5) and provide analgesia in animal models of pain. Here, we have used functional analysis of mutant transporters combined with molecular dynamics simulations of lipid-transporter interactions to understand how these bioactive lipids interact with GlyT2. This study identifies a novel extracellular allosteric modulator site formed by a crevice between transmembrane domains 5, 7, and 8, and extracellular loop 4 of GlyT2. Knowledge of this site could be exploited further in the development of drugs to treat pain, and to identify other allosteric modulators of the SLC6 family of transporters.

Mutation p.R356Q in the Collybistin Phosphoinositide Binding Site Is Associated With Mild Intellectual Disability.

The recruitment of inhibitory GABAA receptors to neuronal synapses requires a complex interplay between receptors, neuroligins, the scaffolding protein gephyrin and the GDP-GTP exchange factor collybistin (CB). Collybistin is regulated by protein-protein interactions at the N-terminal SH3 domain, which can bind neuroligins 2/4 and the GABAAR α2 subunit. Collybistin also harbors a RhoGEF domain which mediates interactions with gephyrin and catalyzes GDP-GTP exchange on Cdc42. Lastly, collybistin has a pleckstrin homology (PH) domain, which binds phosphoinositides, such as phosphatidylinositol 3-phosphate (PI3P/PtdIns3P) and phosphatidylinositol 4-monophosphate (PI4P/PtdIns4P). PI3P located in early/sorting endosomes has recently been shown to regulate the postsynaptic clustering of gephyrin and GABAA receptors and consequently the strength of inhibitory synapses in cultured hippocampal neurons. This process is disrupted by mutations in the collybistin gene (ARHGEF9), which cause X-linked intellectual disability (XLID) by a variety of mechanisms converging on disrupted gephyrin and GABAA receptor clustering at central synapses. Here we report a novel missense mutation (chrX:62875607C>T, p.R356Q) in ARHGEF9 that affects one of the two paired arginine residues in the PH domain that were predicted to be vital for binding phosphoinositides. Functional assays revealed that recombinant collybistin CB3SH3- R356Q was deficient in PI3P binding and was not able to translocate EGFP-gephyrin to submembrane microaggregates in an in vitro clustering assay. Expression of the PI3P-binding mutants CB3SH3- R356Q and CB3SH3- R356N/R357N in cultured hippocampal neurones revealed that the mutant proteins did not accumulate at inhibitory synapses, but instead resulted in a clear decrease in the overall number of synaptic gephyrin clusters compared to controls. Molecular dynamics simulations suggest that the p.R356Q substitution influences PI3P binding by altering the range of structural conformations adopted by collybistin. Taken together, these results suggest that the p.R356Q mutation in ARHGEF9 is the underlying cause of XLID in the probands, disrupting gephyrin clustering at inhibitory GABAergic synapses via loss of collybistin PH domain phosphoinositide binding.

Is protein structure enough? A review of the role of lipids in SLC6 transporter function.

SLC6 neurotransmitter transporters facilitate the Na+- and Cl--dependent uptake of amino acids and amino acid derivatives into cells. Disrupting transport leads to a range of neurological disorders. However, the SLC6 substrate transport mechanism is a topic of ongoing debate. Here, we review the prominent SLC6 substrate transport mechanisms through the lens of molecular dynamics simulations. SLC6 transporters are membrane proteins, yet their transport mechanism(s) have largely been studied without considering the impacts of synaptic lipid composition, or endogenous lipid modulators, on transporter structure and function. In this review, we highlight the importance of studying membrane transporters in an appropriate membrane model, and present opportunities for the community to glean understanding and insight into SLC6 transporter structure and function-in particular transport mechanism(s)-when both membrane lipids and endogenous lipid modulators are considered.

The effects of oxidised phospholipids and cholesterol on the biophysical properties of POPC bilayers.

Oxidation of unsaturated membrane phospholipids by oxidative stress is associated with inflammation, infection, numerous diseases and neurodegenerative disorders. Lipid oxidation is observed in experimental samples when the parent lipid is exposed to oxidative stressors. The effect of phospholipid oxidation on the properties of biological membranes are still being explored, while low concentrations (0.1-2.0 mol%) of oxidised phospholipids are associated with disease states [1]. Previous computational studies have focused on the effect of high concentrations (~50 mol%) of oxidised phospholipids on binary lipid bilayers. This work systematically characterises the effect of lower concentrations (~10 mol%) of two oxidised lipid species, PoxnoPC (1-palmitoyl-2-(9'-oxo-nonanoyl)-sn-glycero-3-phosphocholine) or PazePC (1-palmitoyl-2-azelaoyl-sn-glycero-3-phosphocholine), on POPC/cholesterol and pure POPC bilayers. During µs atomistic simulations in pure POPC bilayers, PoxnoPC and PazePC reoriented their oxidised sn-2 acyl chains towards the solution, and PazePC adopted an extended conformation. The addition of 20 mol% cholesterol not only modulated the fluidity of the bilayers; it also modulated the flexibility of the PoxnoPC oxidised sn-2 tail, reducing bilayer disorder. In contrast, the addition of cholesterol had little effect on bilayers containing PazePC. Our studies show that the effect of oxidised lipids on the biophysical properties of a multicomponent bilayer cannot be intuitively extrapolated from a binary lipid system.

Lipid-Based Inhibitors Act Directly on GlyT2.

The endogenous lipids N-arachidonylglycine and oleoyl-l-carnitine are potential therapeutic leads in the treatment of chronic pain through their inhibition of the glycine transporter GlyT2. However, their mechanism of action is unknown. It has been hypothesized that these "bioactive" lipids either inhibit GlyT2 indirectly, by significantly perturbing the biophysical properties of the membrane, or directly, by binding directly to the transporter (either from a membrane-exposed or solvent-exposed binding site). Here, we used molecular dynamics simulations to study the effects of the lipids anandamide, N-arachidonylglycine, and oleoyl-l-carnitine on (a) the biophysical properties of the bilayer and (b) direct binding interactions with GlyT2. During the simulations, the biophysical properties of the bilayer itself, for example, the area per lipid, bilayer thickness, and order parameters, were not significantly altered by the presence or type of bioactive lipid, regardless of the presence of GlyT2. Our work, together with previous computational and experimental data, suggests that these acyl-inhibitors of GlyT2 inhibit the transporter by directly binding to it. However, these bioactive lipids bound to various parts of GlyT2 and did not prefer a single binding site during 4.5 µs of simulation. We postulate that the binding site is located at the solvent-exposed regions of GlyT2. Understanding the mechanism of action of these and related bioactive lipids is essential in effectively developing high-affinity GlyT2 inhibitors for the treatment of pain.

Probing the Pharmacological Binding Sites of P-Glycoprotein Using Umbrella Sampling Simulations.

The human multidrug transporter P-glycoprotein (P-gp) transports over 200 chemically diverse substrates, influencing their bioavailability and tissue distribution. Pharmacological studies have identified both competitive and noncompetitive P-gp substrates, but neither the precise location of the substrate binding sites, nor the basis of competitive and noncompetitive interactions has been fully characterized. Here, potential of mean force (PMF) calculations are used to identify the transport-competent minimum free energy binding locations of five compounds, Hoechst 33342, Rhodamine 123, paclitaxel, tariquidar, and verapamil to P-gp. Unrestrained molecular dynamics simulations were also performed to confirm the substrates were stable in the energy wells determined using the PMF calculations. All compounds had energy minima within the P-gp transmembrane (TM) pore. For Hoechst 33342 and Rhodamine 123, a second minimum outside the TM pore was also identified. Based on this and previous studies of nicardipine and morphine [ Subramanian et al. J. Chem. Inf. Model. 2015 , 55 , 1202 ], a general scheme that accounts for the observed noncompetitive and competitive substrate interactions with P-gp is proposed.

A potential new, stable state of the E-cadherin strand-swapped dimer in solution.

E-cadherin is a transmembrane glycoprotein that facilitates inter-cellular adhesion in the epithelium. The ectodomain of the native structure is comprised of five repeated immunoglobulin-like domains. All E-cadherin crystal structures show the protein in one of three alternative conformations: a monomer, a strand-swapped trans homodimer and the so-called X-dimer, which is proposed to be a kinetic intermediate to forming the strand-swapped trans homodimer. However, previous studies have indicated that even once the trans strand-swapped dimer is formed, the complex is highly dynamic and the E-cadherin monomers may reorient relative to each other. Here, molecular dynamics simulations have been used to investigate the stability and conformational flexibility of the human E-cadherin trans strand-swapped dimer. In four independent, 100 ns simulations, the dimer moved away from the starting structure and converged to a previously unreported structure, which we call the Y-dimer. The Y-dimer was present for over 90% of the combined simulation time, suggesting that it represents a stable conformation of the E-cadherin dimer in solution. The Y-dimer conformation is stabilised by interactions present in both the trans strand-swapped dimer and X-dimer crystal structures, as well as additional interactions not found in any E-cadherin dimer crystal structures. The Y-dimer represents a previously unreported, stable conformation of the human E-cadherin trans strand-swapped dimer and suggests that the available crystal structures do not fully capture the conformations that the human E-cadherin trans homodimer adopts in solution.

Understanding the accumulation of P-glycoprotein substrates within cells: The effect of cholesterol on membrane partitioning.

The apparent activity of the multidrug transporter P-glycoprotein (P-gp) is enhanced by the presence of cholesterol. Whether this is due to the direct effect of cholesterol on the activity of P-gp, its effect on the local concentration of substrate in the membrane, or its effect on the rate of entry of the drug into the cell, is unknown. In this study, molecular dynamics simulation techniques coupled with potential of mean force calculations have been used to investigate the role of cholesterol in the movement of four P-gp substrates across a POPC bilayer in the presence or absence of 10% cholesterol. The simulations suggest that the presence of cholesterol lowers the free energy associated with entering the middle of the bilayer in a substrate-specific manner. These findings suggest that P-gp substrates may preferentially accumulate in cholesterol-rich regions of the membrane, which may explain its enhanced transport activity.