Ligand induced conformational states of the 5-HT(1A) receptor

Comparative analyses of site-directed mutagenesis of human melatonin MTNR1A and MTNR1B receptors using a yeast fluorescent biosensor

Melatonin is an indoleamine neurohormone made by the pineal gland. Its receptors, MTNR1A and MTNR1B, are members of the G-protein-coupled receptor (GPCR) family and are involved in sleep, circadian rhythm, and mood disorders, and in the inhibition of cancer growth. These receptors, therefore, represent significant molecular targets for insomnia, circadian sleep disorders, and cancer. The yeast Saccharomyces cerevisiae is an attractive host for assaying agonistic activity for human GPCR. We previously constructed a GPCR-based biosensor employing a high-sensitivity yeast strain that incorporated both a chimeric yeast-human Gα protein and a bright fluorescent reporter gene (ZsGreen). Similar approaches have been used for simple and convenient measurements of various GPCR activities. In the current study, we constructed a fluorescence-based yeast biosensor for monitoring the signaling activation of human melatonin receptors. We used this system to analyze point mutations, including previously unreported mutations of the consensus sequences of MTNR1A and MTNR1B melatonin receptors and compared their effects. Most mutations in the consensus sequences significantly affected the signaling capacities of both receptors, but several mutations showed differences between these subtype receptors. Thus, this yeast biosensor holds promise for revealing the functions of melatonin receptors.

Binding Interactions of Ergotamine and Dihydroergotamine to 5-Hydroxytryptamine Receptor 1B (5-HT1b) Using Molecular Dynamics Simulations and Dynamic Network Analysis

Ergotamine (ERG) and dihydroergotamine (DHE), common migraine drugs, have small structural differences but lead to clinically important distinctions in their pharmacological profiles. For example, DHE is less potent than ERG by about 10-fold at the 5-hydroxytrptamine receptor 1B (5-HT_1B). Although the high-resolution crystal structures of the 5-HT_1B receptor with both ligands have been solved, the high similarity between these two complex structures does not sufficiently explain their activity differences and the activation mechanism of the receptor. Hence, an examination of the dynamic motion of both drugs with the receptor is required. In this study, we ran a total of 6.0 μs molecular dynamics simulations on each system. Our simulation data show the subtle variations between the two systems in terms of the ligand-receptor interactions and receptor secondary structures. More importantly, the ligand and protein root-mean-square fluctuations (RMSFs) for the two systems were distinct, with ERG having a trend of lower RMSF values, indicating it to be bound tighter to 5-HT_1B with less fluctuations. The molecular mechanism-general born surface area (MM-GBSA) binding energies illustrate this further, proving ERG has an overall stronger MM-GBSA binding energy. Analysis of several different microswitches has shown that the 5-HT_1B-ERG complex is in a more active conformation state than 5-HT_1B-DHE, which is further supported by the dynamic network model, with reference to mutagenesis data with the critical nodes and the first three low-energy modes from the normal mode analysis. We also identify Trp327^6.48 and Phe331^6.52 as key residues involved in the active state 5-HT_1B for both ligands. Using the detailed dynamic information from our analysis, we made predictions for possible modifications to DHE and ERG that yielded five derivatives that might have more favorable binding energies and reduced structural fluctuations.

Computational and biological evidences on the serotonergic involvement of SeTACN antidepressant-like effect in mice

A series of phenylselanyl-1H-1,2,3-triazole-4-carbonitriles with different substituents were screened for their binding affinity with serotonin transporter (SERT) and dopamine transporter (DAT) by docking molecular. 5-(4methoxyphenyl)-1-(2-(phenylselanyl)phenyl)-1H-1,2,3-triazole-4-carbonitrile (SeTACN) exhibited the best conformation with SERT even higher than fluoxetine and serotonin, suggesting a competitive inhibition. SeTACN demonstrated additional affinity to other serotonergic receptors involved in antidepressant effects: 5HT_1a, 5HT_2a and 5HT₃. In another set of experiments, SeTACN led to significant reductions in the immobility time of mice submitted to forced swimming test (FST) in the dose range of 0.1- 20mg/kg, suggesting an antidepressant-like effect. The possible mechanism of action was investigated using serotonergic and dopaminergic antagonists. The antidepressant-like effect of SeTACN (0.1mg/kg i.g.) was prevented by the pretreatment with WAY100635 (a selective 5HT_1a antagonist), ketanserin (a 5HT2_a/c antagonist) and ondansetron (a selective 5ht₃ antagonist), PCPA (an inhibitor of serotonin synthesis) but not with SCH23390 (dopaminergic D₁ antagonist) and sulpiride (D₂ antagonist). Sub-effective dose of fluoxetine was able to potentiate the effects of a sub-effective dose of SeTACN in FST. None of the treatments affected locomotor activity in open field test (OFT). These results together, suggest that the SeTACN antidepressant-like effect is mediate, at least in parts, by serotonergic system.

From Homology Models to a Set of Predictive Binding Pockets-a 5-HT1A Receptor Case Study

Despite its remarkable importance in the arena of drug design, serotonin 1A receptor (5-HT1A) has been elusive to the X-ray crystallography community. This lack of direct structural information not only hampers our knowledge regarding the binding modes of many popular ligands (including the endogenous neurotransmitter-serotonin), but also limits the search for more potent compounds. In this paper we shed new light on the 3D pharmacological properties of the 5-HT1A receptor by using a ligand-guided approach (ALiBERO) grounded in the Internal Coordinate Mechanics (ICM) docking platform. Starting from a homology template and set of known actives, the method introduces receptor flexibility via Normal Mode Analysis and Monte Carlo sampling, to generate a subset of pockets that display enriched discrimination of actives from inactives in retrospective docking. Here, we thoroughly investigated the repercussions of using different protein templates and the effect of compound selection on screening performance. Finally, the best resulting protein models were applied prospectively in a large virtual screening campaign, in which two new active compounds were identified that were chemically distinct from those described in the literature.

Synthesis and structure evaluation of new complex butylarylpiperazin-1-yl derivatives

A series of arylpiperazine derivatives of 1,16-diphenyl-19-azahexacyclo-[14.5.1.0^2,15.0^3,8.0^9,14.0^17,21]docosa-2,3,5,7,8,9,11,13,14-nonaene-18,20,22-trione and 4,10-diphenyl-1H,2H,3H,5H-indeno[1,2-f]isoindole-1,3,5-trione was synthesized. The pharmacological profile of compound 4 at the 5-HT_1A receptor was measured by binding assay. The title compounds were tested in cell-based assay against the human immunodeficiency virus type-1. The X-ray crystallographic studies of derivatives 2, 6, 7, 11, 19, and 20 were presented.

Linear and Nonlinear Support Vector Machine for the Classification of Human 5-HT1A Ligand Functionality

Upon binding to a receptor, agonists and antagonists can induce distinct biological functions and thus lead to significantly different pharmacological responses. Thus, in silico prediction or in vitro characterization of ligand agonistic or antagonistic functionalities is an important step toward identifying specific pharmacological therapeutics. In this study, we investigated the molecular properties of agonists and antagonists of human 5-hydroxytryptamine receptor subtype 1A (5-HT1A ). Subsequently, intrinsic functions of these ligands (agonists/antagonists) were modeled by support vector machine (SVM), using five 2D molecular fingerprints and the 3D Topomer distance. Five kernel functions, including linear, polynomial, RBF, Tanimoto and a novel Topomer kernel based on Topomer 3D similarity were used to develop linear and nonlinear classifiers. These classifiers were validated through cross-validation, yielding a classification accuracy ranging from 80.4 % to 92.3 %. The performance of different kernels and fingerprints was analyzed and discussed. Linear and nonlinear models were further interpreted through the illustration of underlying classification mechanism. The computation protocol has been automated and demonstrated through our online service. This study expands the scope and applicability of similarity-based methods in cheminformatics, which are typically used for the identification of active molecules against a target protein. Our findings provide a good starting point for further systematic classifications of other GPCR ligands and for the data mining of large chemical libraries.

Simplified modeling approach suggests structural mechanisms for constitutive activation of the C5a receptor

Molecular modeling of conformational changes occurring in the transmembrane region of the complement factor 5a receptor (C5aR) during receptor activation was performed by comparing two constitutively active mutants (CAMs) of C5aR, NQ (I124N/L127Q), and F251A, to those of the wild-type C5aR and NQ-N296A (I124N/L127Q/N296A), which have the wild-type phenotype. Modeling involved comprehensive sampling of various rotations of TM helices aligned to the crystal template of the dark-adapted rhodopsin along their long axes. By assuming that the relative energies of the spontaneously activated states of CAMs should be lower or at least comparable to energies characteristic for the ground states, we selected the plausible models for the conformational states associated with constitutive activation in C5aR. The modeling revealed that the hydrogen bonds between the side chains of D82-N119, S85-N119, and S131-C221 characteristic for the ground state were replaced by the hydrogen bonds D82-N296, N296-Y300, and S131-R134, respectively, in the activated states. Also, conformational transitions that occurred upon activation were hindered by contacts between the side chains of L127 and F251. The results rationalize the available data of mutagenesis in C5aR and offer the first specific molecular mechanism for the loss of constitutive activity in NQ-N296A. Our results also contributed to understanding the general structural mechanisms of activation in G-protein-coupled receptors lacking the "ionic lock", R(3.50) and E/D(6.30). Importantly, these results were obtained by modeling approaches that deliberately simplify many elements in order to explore potential conformations of GPCRs involving large-scale molecular movements.

Synthesis, 3D-QSAR, and structural modeling of benzolactam derivatives with binding affinity for the D(2) and D(3) receptors

A series of 37 benzolactam derivatives were synthesized, and their respective affinities for the dopamine D(2) and D(3) receptors evaluated. The relationships between structures and binding affinities were investigated using both ligand-based (3D-QSAR) and receptor-based methods. The results revealed the importance of diverse structural features in explaining the differences in the observed affinities, such as the location of the benzolactam carbonyl oxygen, or the overall length of the compounds. The optimal values for such ligand properties are slightly different for the D(2) and D(3) receptors, even though the binding sites present a very high degree of homology. We explain these differences by the presence of a hydrogen bond network in the D(2) receptor which is absent in the D(3) receptor and limits the dimensions of the binding pocket, causing residues in helix 7 to become less accessible. The implications of these results for the design of more potent and selective benzolactam derivatives are presented and discussed.

Phe369(7.38) at human 5-HT(7) receptors confers interspecies selectivity to antagonists and partial agonists

BACKGROUND AND PURPOSE

Human and rat 5-HT(7) receptors were studied with a particular emphasis on the molecular interactions involved in ligand binding, searching for an explanation to the interspecies selectivity observed for a set of compounds. We performed affinity studies, molecular modelling and site-directed mutagenesis, with special focus on residue Phe(7.38) of the human 5-HT(7) receptor [Cys(7.38) in rat].

EXPERIMENTAL APPROACH

Competition binding studies were performed for seven 5-HT(7) receptor ligands at three different 5-HT(7) receptors. The functional behaviour was evaluated by measuring 5-carboxytryptamine-stimulated cAMP production. Computational simulations were carried out to explore the structural bases in ligand binding observed for these compounds.

KEY RESULTS

Competition experiments showed a remarkable selectivity for the human receptor when compared with the rat receptor. These results indicate that mutating Cys to Phe at position 7.38 profoundly affects the binding affinities at the 5-HT(7) receptor. Computational simulations provide a structural interpretation for this key finding. Pharmacological characterization of compounds mr25020, mr25040 and mr25053 revealed a competitive antagonistic behaviour. Compounds mr22423, mr22433, mr23284 and mr25052 behaved as partial agonists.

CONCLUSIONS AND IMPLICATIONS

We propose that the interspecies difference in binding affinities observed for the compounds at human and rat 5-HT(7) receptors is due to the nature of the residue at position 7.38. Our molecular modelling simulations suggest that Phe(7.38) in the human receptor is integrated in the hydrophobic pocket in the central part of the binding site [Phe(6.51)-Phe(6.52)] and allows a tighter binding of the ligands when compared with the rat receptor.

Stereoselectivity of 8-OH-DPAT toward the serotonin 5-HT1A receptor: Biochemical and molecular modeling study

The great majority of pharmacological investigations of 5-HT1A receptors' reactivity has been performed using racemic 8-OH-DPAT, therefore the biochemical as well as behavioral profiles of both 8-OH-DPAT enantiomers are not circumstantiated. In the biochemical study capability of racemic 8-OH-DPAT (0.05, 0.1 mg/kg s.c.) and its counterparts R-8-OH-DPAT (0.05, 0.1 mg/kg s.c.) and S-8-OH-DPAT (0.05, 0.1 mg/kg s.c.) to influence 5-HT synthesis rate in rats' prefrontal cortex, hypothalamus, hippocampus and brainstem was evaluated by HPLC/ED technique. Biochemical results are supported by the exhaustive computational study of possible differences between R- and S-enantiomer toward the 5-HT1A receptor. A reliable 3D model of the rat 5-HT1A receptor was constructed from the amino acid sequence using the crystal structure of bovine rhodopsin as a structural template. The structure of the receptor model was validated through docking studies and molecular dynamics simulations that gave results consistent with experimental data. Docking studies and the dynamics of ligand-receptor complexes emphasized different profiles of both enantiomers at the molecular level. The results of both biochemical and computational studies confirmed that R-enantiomer in contrast to S-8-OH-DPAT acts as full and potent agonist, whilst racemic form may display similar pharmacological profile to R-8-OH-DPAT.

From genomics to drug targets

The era of molecular biology and cloning brought new knowledge about the structure and function of drug receptors, and demonstrated that the term 'receptor' must be distinguished from other molecular drug targets such as enzymes, transporters and ion channels. Analysis of the targets of all current therapeutic drugs has shown that more than 95% of these are proteins. The DNA sequencing of the entire human genome has led to identification of many previously unknown proteins that may represent potential drug targets. In order to understand fully the functional mechanisms of a protein, it is crucial to know its three-dimensional molecular structure. This may be determined experimentally by x-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy or electron microscopy, and computationally by structural bioinformatics and molecular modelling. The molecular targets of nearly all current psychotropic drugs are membrane proteins. These have proven extremely difficult to purify and crystallize due to their amphipathic surface, with a hydrophobic area in contact with membrane phospholipids and polar surface areas in contact with the aqueous phases on both sides of the membrane. We have used molecular modelling methods, based on crystal structures of related proteins, to model various neurotransmitter receptors and transporters. The receptor and transporter models have been used to study their structural properties, functional mechanisms and the molecular mechanisms of action of psychotropic drugs. Our results demonstrate the large structural flexibility of transporter and receptor proteins, with substantial movements and conformational changes taking place during substrate translocation in transporters, and by agonist induced receptor stimulation.

3D modeling of the activated states of constitutively active mutants of rhodopsin

The activated (R*) states in constitutively active mutants (CAMs) of G-protein-coupled receptors (GPCRs) are presumably characterized by lower energies than the resting (R) states. If specific configurations of TM helices differing by rotations along the long transmembrane axes possess energies lower than that in the R state for pronounced CAMs, but not for non-CAMs, these particular configurations of TM helices are candidate 3D models for the R* state. The hypothesis was studied in the case of rhodopsin, the only GPCR for which experimentally determined 3D models of the R and R* states are currently available. Indeed, relative energies of the R* state were significantly lower than that of the R state for the rhodopsin mutants G90D/M257Y and E113Q/M257Y (strong CAMs), but not for G90D, E113Q, and M257Y (not CAMs). Next, the developed build-up procedure successfully identified few similar configurations of the TM helical bundle of G90D/M257Y and E113Q/M257Y as possible candidates for the 3D model of the R* state of rhodopsin, all of them being in good agreement with the model suggested by experiment. Since constitutively active mutants are known for many of GPCRs belonging to the large rhodopsin-like family, this approach provides a way for predicting possible 3D structures corresponding to the activated states of the TM regions of many GPCRs for which CAMs have been identified.

Interaction of arylpiperazine ligands with the hydrophobic part of the 5-HT1A receptor binding site

A flexible docking of a series of arylpiperazine derivatives with structurally different aryl part to the binding site of a model of human 5-HT1A receptor was exercised. The influence of structure and hydrophobic properties of aryl moiety on binding affinities was discussed and a model for ligand binding in the hydrophobic part of the binding site was proposed.

CoMFA methodology in structure-activity analysis of hexahydro- and octahydropyrido[1,2-c]pyrimidine derivatives based on affinity towards 5-HT1A, 5-HT2A and alpha1-adrenergic receptors

Structural features of the pyrido[1,2-c]pyrimidine derivatives with arylpiperazine moiety and their affinities towards 5-HT1A, 5-HT2A and alpha1-adrenergic receptors were analyzed using the CoMFA procedure. On the basis of 3D-QSAR models for the 5-HT2A and alpha1-adrenergic receptors, four compounds with expected better affinity/selectivity were proposed and synthesized. The affinities obtained confirm experimentally the usefulness of CoMFA models. Our results suggest that active conformations adopted by the studied molecules when interacting with the receptors are neutral instead of the protonated ones.

Influence of N-1 substituent properties on binding affinities of arylpiperazines to the binding site of 5-HT1A receptor

Serotonin receptors (5-HTRs), especially the 5-HT1A subtype, have been the subject of intensive research for the past decade, due to their function in human physiology. Several structurally different classes of ligands are known to bind to the 5-HT1A receptor, but arylpiperazine derivatives are among the most important ligands. In the work, docking analyses were used to explain the binding affinities of a series of ligands with different N-1 substituent. All ligands had in common the arylpiperazine structure, while the N-1 subsistent was modified to investigate the influence of ligand structure on its binding affinity. The shape and size, as well as the rigidity of the subsistents were altered to investigate the possible effects on the formation of the receptor - ligand complex.

Molecular mechanisms of constitutive activity: mutations at position 111 of the angiotensin AT1 receptor

A possible molecular mechanism for the constitutive activity of mutants of the angiotensin type 1 receptor (AT1) at position 111 was suggested by molecular modeling. This involves a cascade of conformational changes in spatial positions of side chains along transmembrane helix (TM3) from L112 to Y113 to F117, which in turn, results in conformational changes in TM4 (residues I152 and M155) leading to the movement of TM4 as a whole. The mechanism is consistent with the available data of site-directed mutagenesis, as well as with correct predictions of constitutive activity of mutants L112F and L112C. It was also predicted that the double mutant N111G/L112A might possess basal constitutive activity comparable with that of the N111G mutant, whereas the double mutants N111G/Y113A, N111G/F117A, and N111G/I152A would have lower levels of basal activity. Experimental studies of the above double mutants showed significant constitutive activity of N111G/L112A and N111G/F117A. The basal activity of N111G/I152A was higher than expected, and that of N111G/Y113A was not determined due to poor expression of the mutant. The proposed mechanism of constitutive activity of the AT(1) receptor reveals a novel nonsimplistic view on the general problem of constitutive activity, and clearly demonstrates the inherent complexity of the process of G protein-coupled receptor (GPCR) activation.

Atypical and typical antipsychotic drug interactions with the dopamine D2 receptor

A model of the dopamine D2 receptor was used to study the receptor interactions of dopamine, the typical antipsychotics haloperidol and loxapine, and the atypical antipsychotics clozapine and melperone. The atypical antipsychotics interacted with the halogen atom of the ring system in the direction of the transmembrane helices (TMHs) 2, 3 and 7, while the typical had the corresponding halogen atom in the direction of TMH5. Molecular dynamics simulations indicated that the average helical displacement upon binding increased in the order: typical < atypical < dopamine. Upon binding, the atypical induced larger displacements into TMH5 than did the typical. The typical had stronger non-bonded interactions with the receptor than had the atypical, which is in agreement with the experimental observation that the atypical antipsychotic drugs dissociate faster from the receptor than the typical antipsychotic drugs.

Molecular mechanisms of ligand binding, signaling, and regulation within the superfamily of G-protein-coupled receptors: molecular modeling and mutagenesis approaches to receptor structure and function

The superfamily of G-protein-coupled receptors (GPCRs) could be subclassified into 7 families (A, B, large N-terminal family B-7 transmembrane helix, C, Frizzled/Smoothened, taste 2, and vomeronasal 1 receptors) among mammalian species. Cloning and functional studies of GPCRs have revealed that the superfamily of GPCRs comprises receptors for chemically diverse native ligands including (1) endogenous compounds like amines, peptides, and Wnt proteins (i.e., secreted proteins activating Frizzled receptors); (2) endogenous cell surface adhesion molecules; and (3) photons and exogenous compounds like odorants. The combined use of site-directed mutagenesis and molecular modeling approaches have provided detailed insight into molecular mechanisms of ligand binding, receptor folding, receptor activation, G-protein coupling, and regulation of GPCRs. The vast majority of family A, B, C, vomeronasal 1, and taste 2 receptors are able to transduce signals into cells through G-protein coupling. However, G-protein-independent signaling mechanisms have also been reported for many GPCRs. Specific interaction motifs in the intracellular parts of these receptors allow them to interact with scaffold proteins. Protein engineering techniques have provided information on molecular mechanisms of GPCR-accessory protein, GPCR-GPCR, and GPCR-scaffold protein interactions. Site-directed mutagenesis and molecular dynamics simulations have revealed that the inactive state conformations are stabilized by specific interhelical and intrahelical salt bridge interactions and hydrophobic-type interactions. Constitutively activating mutations or agonist binding disrupts such constraining interactions leading to receptor conformations that associates with and activate G-proteins.

Cholesterol modulates ligand binding and G-protein coupling to serotonin(1A) receptors from bovine hippocampus

The serotonin(1A) (5-HT(1A)) receptor is an important member of the superfamily of seven-transmembrane domain G-protein-coupled receptors. We have examined the modulatory role of cholesterol on the ligand binding activity and G-protein coupling of the bovine hippocampal 5-HT(1A) receptor by depleting cholesterol from native membranes using methyl-beta-cyclodextrin (MbetaCD). Removal of cholesterol from bovine hippocampal membranes using varying concentrations of MbetaCD results in a concentration-dependent reduction in specific binding of the agonist 8-OH-DPAT to 5-HT(1A) receptors. This is accompanied by alterations in binding affinity and sites obtained from analysis of binding data. Importantly, cholesterol depletion affected G-protein-coupling of the receptor as monitored by the GTP-gamma-S assay. The concomitant changes in membrane order were reported by changes in fluorescence polarization of membrane probes such as DPH and TMA-DPH, which are incorporated at different locations (depths) in the membrane. Replenishment of membranes with cholesterol led to recovery of ligand binding activity as well as membrane order to a considerable extent. Our results provide evidence, for the first time, that cholesterol is necessary for ligand binding and G-protein coupling of this important neurotransmitter receptor. These results could have significant implications in understanding the influence of the membrane lipid environment on the activity and signal transduction of other G-protein-coupled transmembrane receptors.

Molecular model of the neural dopamine transporter

The dopamine transporter (DAT) regulates the action of dopamine by reuptake of the neurotransmitter into presynaptic neurons, and is the main molecular target of amphetamines and cocaine. DAT and the Na+/H+ antiporter (NhaA) are secondary transporter proteins that carry small molecules across a cell membrane against a concentration gradient, using ion gradients as energy source. A 3-dimensional projection map of the E. coli NhaA has confirmed a topology of 12 membrane spanning domains, and was previously used to construct a 3-dimensional NhaA model with 12 trans-membrane alpha-helices (TMHs). The NhaA model, and site directed mutagenesis data on DAT, were used to construct a detailed 3-dimensional DAT model using interactive molecular graphics and empiric force field calculations. The model proposes a dopamine transport mechanism involving TMHs 1, 3, 4, 5, 7 and 11. Asp79, Tyr252 and Tyr274 were the primary cocaine binding residues. Binding of cocaine or its analogue, (-)-2beta-carbomethoxy-3beta-(4-fluorophenyl)tropane (CFT), seemed to lock the transporter in an inactive state, and thus inhibit dopamine transport. The present model may be used to design further experimental studies of the molecular structure and mechanisms of DAT and other secondary transporter proteins.

Intrinsic activity and comparative molecular dynamics of buspirone analogues at the 5-HT1A receptors

In CNS, the 5-hydroxytryptamine(1A) (5-HT(1A)) receptors exist in two different populations with different behavioural and physiological effects: (1) somatodendritic autoreceptors located pre-synaptically of 5-HT containing neurons and (2) receptors located post-synaptic to 5-HT containing neurons. Clinical studies have shown that 5-HT(1A) partial agonists have anxiolytic properties, while antagonists of pre-synaptical autoreceptors shorten the onset time of selective serotonin reuptake inhibitors (SSRIs). In the present study, the pre- and post-synaptic activity of structural analogues of buspirone was evaluated in animal models. A three dimensional model of the 5-HT(1A) receptor was used to study their interaction modes and helical displacements upon receptor binding. The predicted receptor-ligand interactions indicated similarities in the receptor binding modes for all buspirone analogues, and no clear relationship between receptor contact residues and activity at pre- and post-synaptic receptors. Comparative molecular dynamics (MD) simulations for 650ps indicated that pre-synaptic antagonistic behaviour is connected to large displacements of transmembrane helix (TMH) 7 upon binding, while pre-synaptic agonistic behaviour is connected to large displacements of TMH2 and small displacements of TMH7. Post-synaptic partial agonist behaviour is connected to large displacements of TMH4 and TMH5 upon binding, while post-synaptic antagonists only slightly displace these helices.

A novel ergot alkaloid as a 5-HT(1A) inhibitor produced by Dicyma sp

In the course of a search for small-molecule inhibitors of 5-hydroxytryptamine receptors we have identified a novel ergoline derivative (1) from the fungal culture of Dicyma sp. This compound has a pK(i) of 10.2 versus the 5-hydroxytryptamine(1A) receptor subtype. The structure was elucidated by extensive NMR spectroscopy and mass spectrometry.

Molecular mechanism of citalopram and cocaine interactions with neurotransmitter transporters

The selective serotonin reuptake inhibitors (SSRIs) and cocaine bind to the neural serotonin (5-HT) transporter (SERT) and thus inhibit presynaptic reuptake of 5-HT and elevate its concentration in the synaptic cleft. Cocaine also binds to the dopamine transporter (DAT) and to the noradrenaline transporter (NET) and inhibits presynaptic reuptake of dopamine and noradrenaline. SERT, DAT, and NET belong to the sodium/neurotransmitter symporter family, which is predicted to have a molecular structure with 12 transmembrane alpha-helices (TMHs) and intracellular amino- and carboxy terminals. We used an electron density projection map of the Escherichia coli Na+/H+ anti-porter, and site-directed mutagenesis data on DAT and SERT to construct 3-dimensional molecular models of SERT, DAT and NET. These models were used to simulate the molecular interaction mechanisms of the SSRI, S-citalopram, its less potent enantiomer, R-citalopram and of cocaine with the transporters. In the SERT model, a single amino acid (Tyr95) in TMH1 determined the transporter selectivity of S-citalopram for SERT over DAT and NET. A dipole-dipole interaction was formed between the hydroxy group of Tyr95 in SERT and the nitril group of S-citalopram, but could not be formed by S-citalopram in DAT and NET where the corresponding amino acid is a phenylalanine. The lower binding affinity of R-citalopram may be due to sterical hindrance at the binding site. The tropane ring of cocaine interacted with Tyr95 in SERT and with the corresponding phenylalanines in NET and DAT. This may explain why cocaine, but not S-citalopram, has high binding affinity to all three transporters.

Molecular Dynamics Simulations of the Ligand-Induced Chemical Information Transfer in the 5-HT1A Receptor

Comparative molecular dynamics simulations of the 5-HT(1A) receptor in its empty as well as agonist- (i.e. active) and antagonist-bound (i.e. nonactive) forms have been carried out. The agonists 5-HT and (R)-8-OH-DPAT as well as the antagonist WAY100635 have been employed. The results of this study strengthen the hypothesis that the receptor portions close to the E/DRY/W motif, with prominence to the cytosolic extensions of helices 3 and 6, are particularly susceptible to undergo structural modification in response to agonist binding. Despite the differences in the structural/dynamics behavior of the two agonists when docked into the 5-HT(1A) receptor, they both exert a destabilization of the intrahelical and interhelical interactions found in the empty and antagonist-bound receptor forms between the arginine of the E/DRY sequence and both D133(3.49) and E340(6.30). For both agonists, the chemical information transfer from the extracellular to the cytosolic domains is mediated by a cluster of aromatic amino acids in helix 6, following the ligand interaction with selected amino acids in the extracellular half of the receptor, such as D116(3.32), S199(5.42), Y195(5.38), and F361(6.51). A significant reduction in the bend at P360(6.50), as compared to the empty and the antagonist-bound receptor forms, is one of the features of the agonist-bound forms that is related to the breakage of the interhelical salt bridge between the E/DRY arginine and E340(6.30). Another structural feature, shared by the agonist-bound receptor forms and not by the empty and antagonist-bound forms, is the detachment of helices 2 and 4, as marked by the movement of W161(4.50) away from helix 2, toward the membrane space.

Synthesis and molecular modeling of new 1-aryl-3-[4-arylpiperazin-1-yl]-1-propane derivatives with high affinity at the serotonin transporter and at 5-HT(1A) receptors

It has been proposed that 5-HT(1A) receptor antagonists augment the antidepressant efficacy of selective serotonin (5-HT) reuptake inhibitors. In a search toward new and efficient antidepressants, 1-(aryl)-3-[4-arylpiperazin-1-yl]-1-propane molecular hybrids were designed, synthesized, and evaluated for 5-HT reuptake inhibition and 5-HT(1A) receptor affinity. The design was based in coupling structural moieties related to inhibition of serotonin reuptake, such as benzo[b]thiophene derivatives to arylpiperazines, typical 5-HT(1A) receptor ligands. In binding studies, several compounds showed affinity at the 5-HT transporter and at 5-HT(1A) receptors. Molecular modeling studies predicted the pharmacophore elements required for high affinity binding and the features that enable to discriminate between agonist, partial agonist, or antagonist action at 5-HT(1A) receptors and 5-HT transporter inhibition. Solvent interactions in desolvation prior to the binding step along with enthalpy and enthropy compensations might be responsible to explain agonist, partial agonist, and antagonist character. Hydrogen-bonding capability seems to be important to break hydrogen interhelical hydrogen bonds or alternatively to form other bonds upon ligand binding. Partial agonists and antagonists are unable to do this as the full agonist, which interacts closely by long-range forces or directly. The compounds showing the higher affinity at both the 5-HT transporter (K(i) < 50 nM) and the 5-HT(1A) receptors (K(i) < 20 nM) were further explored for their ability to stimulate [(35)S]GTPgammaS binding or to antagonize 8-hydroxy-2-di-n-propylamino-tetralin (8-OH-DPAT)-stimulated [(35)]GTPgammaS binding to rat hippocampal membranes, an index of agonist/antagonist action at 5-HT(1A) receptors, respectively. Compound 8g exhibited agonist activity (EC(50) = 30 nM) in this assay, whereas compounds 7g and 8h,i behaved as weak partial agonists and 7h-j and 8j,l antagonized the R(+)-8-OH-DPAT-stimulated GTPgammaS binding. Functional characterization was performed by measuring the antagonism to 8-OH-DPAT-induced hypothermia in mice.

Bioinformatics and receptor mechanisms of psychotropic drugs

One important aspect in biotechnology is gene discovery and target validation for drug discovery. Information from the human genome (HUGO) project may be used to deduce the amino acid sequence of all proteins produced in the human body. However, knowing the amino acid sequence of a protein is not the same as knowing its function. Identification of novel molecular targets for discovery of new, safer and more efficient therapeutic drugs from the human genome sequences requires multidisciplinary research efforts, including proteomics, structural biology and bioinformatics. In addition to possible effects on gene expression, most of the currently used therapeutic drugs either have enzymes or membrane proteins as their molecular targets of action. These membrane proteins include transporters of small molecules across cell membranes, ion channels, or receptors that convey signals from one side of a membrane to the other. Our research group as well as others have used computational techniques, along with biotechnology, molecular biology and other experimental techniques, to construct detailed 3-dimensional models of transporter proteins and G-protein coupled receptors (GPCRs), which are the molecular targets of action of psychotropic drugs. The models have been used to simulate the molecular dynamics and study the ligand binding and signal transduction mechanisms of these receptors. The use of bioinformatics, as exemplified in our modelling of GPCRs, is only one of the key factors for success in post-genomic research for new targets for therapeutic drugs.

Interaction of 1,2,4-substituted piperazines, new serotonin receptor ligands, with 5-HT1A and 5-HT2A receptors

In the present paper, we describe affinities to 5-HT1A and 5-HT2A receptors of several new 1,2,4-trisubstituted piperazine derivatives. The affinities were compared with those described earlier for 1,4-disubstituted piperazines and the influence of the third (methyl) substituent on the affinity to both receptors is discussed. The difference between two- and three-substituted derivatives was rationalised in terms of molecular modelling of the respective ligand-receptor complexes. Additionally, the functional activity of some 1,2,4-trisubstituted piperazines for 5-HT1A receptor was examined in behavioural and biochemical models. The obtained results have shown that some trisubstituted compounds exhibited a higher affinity to 5-HT2A receptors than their respective disubstituted analogues (with the affinity to 5-HT1A receptors remaining the same or somewhat improving). The molecular dynamics simulations suggested that the presence of the third substituent in the piperazine ring of those compounds may induce stabilising effect on the ligand-receptor complexes. The results of the in vivo studies have shown that some of the examined trisubstituted piperazines (10-13, 16, 17) exhibited properties of postsynaptic 5-HT1A partial agonists. Moreover, compounds 13 and 16 exhibited features of 5-HT1A presynaptic agonists in in vitro test, and compound 16 also in in vivo tests.

Bioinformatics: from genome to drug targets

The complete sequence determination of the human genome marks the start of a new era in biological science, with focus shifting from sequencing to functional mechanisms of gene products. In addition to effects on gene expression, most of the currently used therapeutic drugs either have enzymes or membrane proteins as their molecular targets of action. These membrane proteins include ion channels and transporters of small molecules, and receptors that convey signals from one side of a membrane to the other. Membrane proteins are thus involved in a variety of cellular processes and have a large potential as targets for new drug discovery. However, detailed structural information is still lacking for the majority of membrane proteins since their association with membrane constituents make NMR (nuclear magnetic resonance) spectroscopic and X-ray diffraction determinations difficult. Molecular modelling by biocomputing is a methodological alternative for structural studies of membrane proteins, but has to be based on experimental structural information in addition to computational techniques. A combination of bioinformatics and experimental techniques was used to model membrane proteins from two different classes, secondary transporters of the sodium:neurotransmitter symporter family (SNF transporters), and G-protein coupled receptors (GPCRs). The protein models were used to examine ligand-protein interactions and signalling/transport mechanisms, and to design experimental site-directed mutagenesis studies. Such studies have provided new insight into the detailed molecular mechanisms of two important classes of membrane proteins, which may be of value in the discovery and development of new pharmaceuticals.

Molecular dynamics of 5-HT1A and 5-HT2A serotonin receptors with methylated buspirone analogues

In the present study experimentally determined ligand selectivity of three methylated buspirone analogues (denoted as MM2, MM5 and P55) towards 5-HT1A and 5-HT2A serotonin receptors was theoretically investigated on a molecular level. The relationships between the ligand structure and 5-HT1A and 5-HT2A receptor affinities were studied and the results were found to be in agreement with the available site-directed mutagenesis and binding affinity data. Molecular dynamics (MD) simulations of ligand-receptor complexes were performed for each investigated analogue, docked twice into the central cavity of 5-HT1A/5-HT2A, each time in a different orientation. Present results were compared with our previous theoretical results, obtained for buspirone and its non-methylated analogues. It was found that due to the presence of the methyl group in the piperazine ring the ligand position alters and the structure of the ligand-receptor complex is modified. Further, the positions of derivatives with pyrimidinyl aromatic moiety and quinolinyl moiety are significantly different at the 5-HT2A receptor. Thus, methylation of such derivatives alters the 3D structures of ligand-receptor complexes in different ways. The ligand-induced changes of the receptor structures were also analysed. The obtained results suggest, that helical domains of both receptors have different dynamical behaviour. Moreover, both location and topography of putative binding sites for buspirone analogues are different at 5-HT1A and 5-HT2A receptors.