Molecular modeling and structure-activity relationship studies were performed to propose a model for binding of the neurotransmitter serotonin (5-HT) to the human serotonin transporter (hSERT). Homology models were constructed using the crystal structure of a bacterial homologue, the leucine transporter from Aquifex aeolicus, as the template and three slightly different sequence alignments. Induced fit docking of 5-HT into hSERT homology models resulted in two different binding modes. Both show a salt bridge between Asp98 and the charged primary amine of 5-HT, and both have the 5-HT C6 position of the indole ring pointing toward Ala173. The difference between the two orientations of 5-HT is an enantiofacial discrimination of the indole ring, resulting in the 5-hydroxyl group of 5-HT being vicinal to either Ser438/Thr439 or Ala169/Ile172/Ala173. To assess the binding experimentally, binding affinities for 5-HT and 17 analogues toward wild type and 13 single point mutants of hSERT were measured using an approach termed paired mutant-ligand analogue complementation (PaMLAC). The proposed ligand-protein interaction was systematically examined by disrupting it through site-directed mutagenesis and re-establishing another interaction via a ligand analogue matching the mutated residue, thereby minimizing the risk of identifying indirect effects. The interactions between Asp98 and the primary amine of 5-HT and the interaction between the C6-position of 5-HT and hSERT position 173 was confirmed using PaMLAC. The measured binding affinities of various mutants and 5-HT analogues allowed for a distinction between the two proposed binding modes of 5-HT and biochemically support the model for 5-HT binding in hSERT where the 5-hydroxyl group is in close proximity to Thr439.

The plasma membrane, which encapsulates human cells, is composed of a complex mixture of lipids and embedded proteins. Emerging knowledge points towards the lipids as having a regulating role in protein function. Furthermore, insight from protein crystallography has revealed several different types of lipids intimately bound to membrane proteins and peptides, hereby possibly pointing to a site of action for the observed regulation. Cholesterol is among the lipid membrane constituents most often observed to be co-crystallized with membrane proteins, and the cholesterol levels in cell membranes have been found to play an essential role in health and disease. Remarkably little is known about the mechanism of lipid regulation of membrane protein function in health as well as in disease. Herein, we review molecular dynamics simulation studies aimed at investigating the effect of cholesterol on membrane protein and peptide properties. This article is part of a Special Issue entitled: Lipid-protein interactions.

We present a simulation study where different resolutions, namely coarse-grained (CG) and all-atom (AA) molecular dynamics simulations, are used sequentially to combine the long timescale reachable by CG simulations with the high resolution of AA simulations, to describe the complete processes of peptide aggregation and pore formation by alamethicin peptides in a hydrated lipid bilayer. In the 1-micros CG simulations the peptides spontaneously aggregate in the lipid bilayer and exhibit occasional transitions between the membrane-spanning and the surface-bound configurations. One of the CG systems at t = 1 micros is reverted to an AA representation and subjected to AA simulation for 50 ns, during which water molecules penetrate the lipid bilayer through interactions with the peptide aggregates, and the membrane starts leaking water. During the AA simulation significant deviations from the alpha-helical structure of the peptides are observed, however, the size and arrangement of the clusters are not affected within the studied time frame. Solid-state NMR experiments designed to match closely the setup used in the molecular dynamics simulations provide strong support for our finding that alamethicin peptides adopt a diverse set of configurations in a lipid bilayer, which is in sharp contrast to the prevailing view of alamethicin oligomers formed by perfectly aligned helical alamethicin peptides in a lipid bilayer.

The serotonergic system plays an important role in many psychiatric disorders. Its role in depression is well established (1). The majority of antidepressants, including TCAs, 6 cause increased synaptic serotonin (5-HT) levels via blockade of 5-HT reuptake into the presynaptic neuron (2-4) by competitive inhibition of hSERT. TCAs have been in clinical use since the 1950s, with imipramine being the first and most prominent compound (5). In severely depressed hospitalized patients, TCAs appear to be more efficacious than selective serotonin reuptake inhibitors (6). TCAs remain in widespread clinical use, especially for treatment-resistant depression (7).hSERT belongs to the neurotransmitter sodium symporter family (2, 8). These transporters utilize the electrochemical gradient of sodium and chloride ions to accumulate 5-HT against its own gradient (9 -11). No experimentally solved structures of the monoamine transporters exist, including hSERT and the dopamine and norepinephrine transporters. However, the structure of LeuT, a bacterial homolog of the neurotransmitter sodium symporters, in a substrate-occluded conformation, was reported in 2005 (12). Two sodium ions (12) and a chloride ion bind near the central substrate site (13-14) structurally and functionally coupling sodium and chloride binding to substrate binding. Recently, different transport mechanisms have been suggested (15,16).Subsequently, a low affinity noncompetitive binding site for TCAs in the extracellular vestibule of the LeuT 11 Å above the central binding site was identified (17,18). The relevance of the LeuT vestibular site for TCA binding to the physiologically relevant target, hSERT, is a matter of debate. This study identifies the central binding site, not the putative vestibular site, as relevant for TCA binding to hSERT and furthermore describes and validates the orientation of TCAs within this site.In this paper, we present induced fit docking studies of imipramine and selected analogs in the previously described homology models of hSERT (19). We present binding affinity studies of 10 imipramine analogs (Fig. 1 Copenhagen Ø, Denmark. 4 To whom correspondence may be addressed. E-mail: birgit@chem.au.dk. 5 To whom correspondence may be addressed. E-mail: owiborg@post.tele.dk. 6 The abbreviations used are: TCA, tricyclic antidepressant; 5-HT, serotonin; hSERT, human SERT; WT, wild type; PaMLAC, paired mutation ligand analog complementation; MD, molecular dynamics; r.m.s., root mean square.

We study the phase behavior of saturated lipids as a function of temperature and tail length for two coarse-grained models: the soft-repulsive model typically employed with dissipative particle dynamics (DPD) and the MARTINI model. We characterize the simulated transitions through changes in structural properties, and we introduce a computational method to monitor changes in enthalpy, as is done experimentally with differential scanning calorimetry. The lipid system experimentally presents four different bilayer phasessubgel, gel, ripple, and fluidand the DPD model describes all of these phases structurally while MARTINI describes a single orderÀdisorder transition between the gel and the fluid phases. Given both models' varying degrees of success in displaying accurate structural and thermodynamic signatures, there is an overall satisfying extent of agreement for the coarse-grained models. We also study the lipid dynamics displayed by these models for the various phases, discussing this dynamics with relation to fidelity to experiment and computational efficiency.