A short update on the structure of drug binding sites on neurotransmitter transporters.

BACKGROUND: The dopamine (DAT), noradrenalin (NET) and serotonin (SERT) transporters are molecular targets for different classes of psychotropic drugs. Cocaine and the SSRI (S)-citalopram block neurotransmitter reuptake competitively, but while cocaine is a non-selective reuptake inhibitor, (S)-citalopram is a selective SERT inhibitor. FINDINGS: Here we present comparisons of the binding sites and the electrostatic potential surfaces (EPS) of DAT, NET and SERT homology models based on two different LeuTAa templates; with a substrate (leucine) in an occluded conformation (PDB id 2a65), and with an inhibitor (tryptophan) in an open-to-out conformation (PDB id 3f3a). In the occluded homology models, two conserved aromatic amino acids (tyrosine and phenylalanine) formed a gate between the putative binding pockets, and this contact was interrupted in the open to out conformation. The EPS of DAT and NET were generally negative in the vestibular area, whereas the EPS of the vestibular area of SERT was more neutral. CONCLUSIONS: The findings presented here contribute as an update on the structure of the binding sites of DAT, NET and SERT. The updated models, which have larger ligand binding site areas than models based on other templates, may serve as improved tools for virtual ligand screening.

Incorporating physiological and biochemical mechanisms into pharmacokinetic-pharmacodynamic models: a conceptual framework.

The aim of this conceptual framework paper is to contribute to the further development of the modelling of effects of drugs or toxic agents by an approach which is based on the underlying physiology and pathology of the biological processes. In general, modelling of data has the purpose (1) to describe experimental data, (2a) to reduce the amount of data resulting from an experiment, e.g. a clinical trial and (2b) to obtain the most relevant parameters, (3) to test hypotheses and (4) to make predictions within the boundaries of experimental conditions, e.g. range of doses tested (interpolation) and out of the boundaries of the experimental conditions, e.g. to extrapolate from animal data to the situation in man. Describing the drug/xenobiotic-target interaction and the chain of biological events following the interaction is the first step to build a biologically based model. This is an approach to represent the underlying biological mechanisms in qualitative and also quantitative terms, thus being inherently connected in many aspects to systems biology. As the systems biology models may contain variables in the order of hundreds connected with differential equations, it is obvious that it is in most cases not possible to assign values to the variables resulting from experimental data. Reduction techniques may be used to create a manageable model which, however, captures the biologically meaningful events in qualitative and quantitative terms. Until now, some success has been obtained by applying empirical pharmacokinetic/pharmacodynamic models which describe direct and indirect relationships between the xenobiotic molecule and the effect, including tolerance. Some of the models may have physiological components built in the structure of the model and use parameter estimates from published data. In recent years, some progress toward semi-mechanistic models has been made, examples being chemotherapy-induced myelosuppression and glucose-endogenous insulin-antidiabetic drug interactions. We see a way forward by employing approaches to bridge the gap between systems biology and physiologically based kinetic and dynamic models. To be useful for decision making, the 'bridging' model should have a well founded mechanistic basis, but being reduced to the extent that its parameters can be deduced from experimental data, however capturing the biological/clinical essential details so that meaningful predictions and extrapolations can be made.

Modelling the genesis and treatment of cancer: the potential role of physiologically based pharmacodynamics.

Physiologically based modelling of pharmacodynamics/toxicodynamics requires an a priori knowledge on the underlying mechanisms causing toxicity or causing the disease. In the context of cancer, the objective of the expert meeting was to discuss the molecular understanding of the disease, modelling approaches used so far to describe the process, preclinical models of cancer treatment and to evaluate modelling approaches developed based on improved knowledge. Molecular events in cancerogenesis can be detected using 'omics' technology, a tool applied in experimental carcinogenesis, but also for diagnostics and prognosis. The molecular understanding forms the basis for new drugs, for example targeting protein kinases specifically expressed in cancer. At present, empirical preclinical models of tumour growth are in great use as the development of physiological models is cost and resource intensive. Although a major challenge in PKPD modelling in oncology patients is the complexity of the system, based in part on preclinical models, successful models have been constructed describing the mechanism of action and providing a tool to establish levels of biomarker associated with efficacy and assisting in defining biologically effective dose range selection for first dose in man. To follow the concentration in the tumour compartment enables to link kinetics and dynamics. In order to obtain a reliable model of tumour growth dynamics and drug effects, specific aspects of the modelling of the concentration-effect relationship in cancer treatment that need to be accounted for include: the physiological/circadian rhythms of the cell cycle; the treatment with combinations and the need to optimally choose appropriate combinations of the multiple agents to study; and the schedule dependence of the response in the clinical situation.

Effects of pharmaceuticals and other active chemicals at biological targets: mechanisms, interactions, and integration into PB-PK/PD models.

Biologically active molecules, that is pharmaceuticals and other chemical substances, exert their therapeutic and/or toxic effects by complex interactions with their biological targets/active sites. This review discusses the factors and processes governing the kinetics and effects of active molecules at their cellular targets, the chain of events leading to clinical effects, and the crosstalk between regulatory pathways controlling these processes. Special attention is given to the discussion of effects of a single drug or other chemical on multiple targets, to the interaction of multiple ligands with a single target/receptor and the effects of single ligand-target complexes on multiple signal transduction pathways, and to the control of physiological functions, such as regulation of blood glucose levels, by numerous primary mechanisms acting on different cellular targets. Physiologically based-pharmacokinetic/pharmacodynamic (PB-PK/PD) models are of great value for the design of active principles by the pharmaceutical industry and for the optimization and individualization of patient therapy. Experimental results from in vitro and in vivo studies can be used for building such models. On the other hand, properly designed models and simulation can contribute to a better design of experiments. Much of what is presented in this article applies equally well to drugs and other chemicals. Unless specified otherwise, reference to drugs applies also to other chemicals. This review is based on an expert meeting organized by COST Action B25 held in Eilat, Israel, on 14 - 15 February 2008. The authors have prepared this article to reflect the presentations and discussions at that meeting.

Structure and localisation of drug binding sites on neurotransmitter transporters.

The dopamine (DAT), serotontin (SERT) and noradrenalin (NET) transporters are molecular targets for different classes of psychotropic drugs. The crystal structure of Aquifex aeolicus LeuT(Aa) was used as a template for molecular modeling of DAT, SERT and NET, and two putative drug binding sites (pocket 1 and 2) in each transporter were identified. Cocaine was docked into binding pocket 1 of DAT, corresponding to the leucine binding site in LeuT(Aa), which involved transmembrane helices (TMHs) 1, 3, 6 and 8. Clomipramine was docked into binding pocket 2 of DAT, involving TMHs 1, 3, 6, 10 and 11, and extracellular loops 4 and 6, corresponding to the clomipramine binding site in a crystal structure of a LeuT(Aa)-clomipramine complex. The structures of the proposed cocaine- and tricyclic antidepressant-binding sites may be of particular interest for the design of novel DAT interacting ligands.

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.

Putative drug binding conformations of monoamine transporters.

Structural information about monoamine transporters and their interactions with psychotropic drugs is important for understanding their molecular mechanisms of action and for drug development. The crystal structure of a Major Facilitator Superfamily (MFS) transporter, the lactose permease symporter (lac permease), has provided insight into the three-dimensional structure and mechanisms of secondary transporters. Based on the hypothesis that the 12 transmembrane alpha-helix (TMH) secondary transporters belong to a common folding class, the lac permease structure was used for molecular modeling of the serotonin transporter (SERT), the dopamine transporter (DAT), and the noradrenaline transporter (NET). The molecular modeling methods used included amino acid sequence alignment, homology modeling, and molecular mechanical energy calculations. The lac permease crystal structure has an inward-facing conformation, and construction of outward-facing SERT, DAT, and NET conformations allowing ligand binding was the most challenging step of the modeling procedure. The psychomotor stimulants cocaine and S-amphetamine, and the selective serotonin reuptake inhibitor (SSRI) S-citalopram, were docked into putative binding sites on the transporters to examine their molecular binding mechanisms. In the inward-facing conformation of SERT the translocation pore was closed towards the extracellular side by hydrophobic interactions between the conserved amino acids Phe105, Pro106, Phe117, and Ala372. An unconserved amino acid, Asp499 in TMH10 in NET, may contribute to the low affinity of S-citalopram to NET.

Mechanism-based pharmacokinetic-pharmacodynamic modeling-a new classification of biomarkers.

In recent years, pharmacokinetic/pharmacodynamic (PK/PD) modeling has developed from an empirical descriptive discipline into a mechanistic science that can be applied at all stages of drug development. Mechanism-based PK/PD models differ from empirical descriptive models in that they contain specific expressions to characterize processes on the causal path between drug administration and effect. Mechanism-based PK/PD models have much improved properties for extrapolation and prediction. As such, they constitute a scientific basis for rational drug discovery and development. In this report, a novel classification of biomarkers is proposed. Within the context of mechanism-based PK/PD modeling, a biomarker is defined as a measure that characterizes, in a strictly quantitative manner, a process, which is on the causal path between drug administration and effect. The new classification system distinguishes seven types of biomarkers: type 0, genotype/phenotype determining drug response; type 1, concentration of drug or drug metabolite; type 2, molecular target occupancy; type 3, molecular target activation; type 4, physiological measures; type 5, pathophysiological measures; and type 6, clinical ratings. In this paper, the use of the new biomarker classification is discussed in the context of the application of mechanism-based PK/PD analysis in drug discovery and development.

Molecular modelling of drug targets: the past, the present and the future.

Most currently used therapeutic drugs have an enzyme or a membrane-bound receptor as site of action. The sequencing of the human and other genomes has provided a potential to identify many hitherto unknown proteins that might serve as new drug targets. To achieve this, knowledge about three-dimensional protein structures is crucial for the understanding of their functional mechanisms, and for a rational drug design. Over the last decade atomic resolution crystal structures of soluble proteins have been reported in a rapidly increasing number, but the detailed three-dimensional structures are still unknown for the majority of membrane proteins since their membrane association makes experimental structure determinations complicated. Computerized modelling of protein structures, based on experimentally determined structures of homologue proteins, may be a useful methodological alternative, especially for membrane proteins. In the past, molecular modelling of transporters and G-protein-coupled receptors was based on low-resolution structural data obtained by cryo-electron microscopy. Recent high-resolution crystal structure determinations of a G-protein-coupled receptor, rhodopsin, and several different transporter proteins and ion channels have enabled construction of more accurate receptor and transporter models. For the future, collaborative structural genomics initiatives aim at determining the three-dimensional structure of all known proteins, based on a combination of experimental structure determination and molecular modelling. Development of still more powerful computer hardware and software will enable extensive studies of the protein structure and dynamics of new potential drug targets, but raises a new challenge in the validation and calibration of computerized methods of biosimulations.

Molecular approaches to the identification of biomarkers of exposure and effect--report of an expert meeting organized by COST Action B15. November 28, 2003.

In the past, the term biomarker has been used with several meanings when used in human and environmental toxicology as compared to pharmaceutical development. However, with the advent of molecular approaches and their application in the field of drug development and toxicology, the concept of biomarkers has to be newly defined. In the meeting, the experts found consent in defining the term and described the application of biomarkers in toxicology, drug development and clinical diagnostics. Molecular approaches to biomarker identification and selection lead to a large amount of data. Hence, the statistical analysis is challenging and special statistical problems have to be solved in biomarker characterization, of particular interest are attempts aiming at class discovery and prediction. Reliability and biological relevance are to be demonstrated for biomarkers of exposure and effect which is also true for biomarkers of susceptibility. It is envisaged that the application of biomarkers will expand from current use in pre-clinical toxicology to the risk characterization and risk assessment of chemicals and from early clinical phases of drug development to later phases and even into daily clinical use in diagnostics and disease classification.

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.

Subtype-specific sorting of the ETA endothelin receptor by a novel endocytic recycling signal for G protein-coupled receptors.

We have previously reported that endocytic sorting of ET(A) endothelin receptors to the recycling pathway is dependent on a signal residing in the cytoplasmic carboxyl-terminal region. The aim of the present work was to characterize the carboxyl-terminal recycling motif of the ET(A) receptor. Assay of truncation mutants of the ET(A) receptor with increasing deletions of the carboxyl-terminal tail revealed that amino acids 390 to 406 contained information critical for the ability of the receptor to recycle. This peptide sequence displayed significant sequence similarity to several protein segments confirmed by X-ray crystallography to adopt antiparallel beta-strand structures (beta-finger). One of these segments was the beta-finger motif of neuronal nitric-oxide synthase reported to function as an internal PDZ (postsynaptic density-95/disc-large/zona occludens) domain-binding ligand. Based on these findings, the three-dimensional structure of the recycling motif of ET(A) receptor was predicted to attain a beta-finger conformation acting as an internal PDZ ligand. Site-directed mutagenesis at residues that would be crucial to the structural integrity of the putative beta-finger conformation or PDZ ligand function prevented recycling of the ET(A) receptor. Analysis of more than 300 G protein-coupled receptors (GPCRs) identified 35 different human GPCRs with carboxylterminal sequence patterns that fulfilled the structural criteria of an internal PDZ ligand. Among these are several receptors reported to follow a recycling pathway. In conclusion, recycling of ET(A) receptor is mediated by a motif with the structural characteristics of an internal PDZ ligand. This structural motif may represent a more general principle of endocytic sorting of GPCRs.

Structures and models of transporter proteins.

Transporter proteins in biological membranes may be divided into channels and carriers. Channels function as selective pores that open in response to a chemical or electrophysiological stimulus, allowing movement of a solute down an electrochemical gradient. Active carrier proteins use an energy producing process to translocate a substrate against a concentration gradient. Secondary active transporters use the movement of a solute down a concentration gradient to drive the translocation of another substrate across a membrane. ATP-binding cassette (ABC) transporters couple hydrolysis of adenosine triphosphate (ATP) to the translocation of various substrates across cell membranes. High-resolution three-dimensional structures have now been reported from X-ray crystallographic studies of six different transporters, including two ATP-binding cassette (ABC) transporters. These structures have explained the results from many previous biochemical and biological studies and shed new light on their functional mechanisms. All these transporters have alpha-helical structures of the membrane-spanning domains, as suggested from many previous studies, and some of the helices have irregular shapes with kinks and bends. Together these crystal structures demonstrate the large flexibility of transporter proteins and that substantial movements take place during the substrate translocation process, which to a certain extent may distinguish active carriers from channel proteins. These structures and other low-resolution structures of membrane proteins have served as a basis for construction of three-dimensional protein models that have provided insight into functional mechanisms and molecular structures and enabled formulation of new hypotheses regarding transporter structure and function, which may be experimentally validated.

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.

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.

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.