Showing papers on "Docking (molecular) published in 2012"

Multidrug binding properties of the AcrB efflux pump characterized by molecular dynamics simulations

Abstract: Multidrug resistance in Gram-negative bacteria, to which multidrug efflux pumps such as the AcrB transporter makes a major contribution, is becoming a major public health problem Unfortunately only a few compounds have been cocrystallized with AcrB, and thus computational approaches are essential in elucidating the interaction between diverse ligands and the pump protein We used molecular dynamics simulation to examine the binding of nine substrates, two inhibitors, and two nonsubstrates to the distal binding pocket of AcrB, identified earlier by X-ray crystallography This approach gave us more realistic views of the binding than the previously used docking approach, as the explicit water molecules contributed to the process and the flexible binding site was often seen to undergo large structural changes We analyzed the interaction in detail in terms of the binding energy, hydrophobic surface-matching, and the residues involved in the process We found that all substrates tested bound to the pocket, whereas the binding to this site was not preferred for the nonsubstrates Interestingly, both inhibitors [Phe-Arg-β-naphthylamide and 1-(1-naphtylmethyl)-piperazine] tended to move out of the pocket at least partially, getting into contact with a glycine-rich loop that separates the distal pocket from the more proximal region of the protein and is thought to control the access of substrates to the distal pocket

Templates are available to model nearly all complexes of structurally characterized proteins

Abstract: Traditional approaches to protein–protein docking sample the binding modes with no regard to similar experimentally determined structures (templates) of protein–protein complexes. Emerging template-based docking approaches utilize such similar complexes to determine the docking predictions. The docking problem assumes the knowledge of the participating proteins’ structures. Thus, it provides the possibility of aligning the structures of the proteins and the template complexes. The progress in the development of template-based docking and the vast experience in template-based modeling of individual proteins show that, generally, such approaches are more reliable than the free modeling. The key aspect of this modeling paradigm is the availability of the templates. The current common perception is that due to the difficulties in experimental structure determination of protein–protein complexes, the pool of docking templates is insignificant, and thus a broad application of template-based docking is possible only at some future time. The results of our large scale, systematic study show that, surprisingly, in spite of the limited number of protein–protein complexes in the Protein Data Bank, docking templates can be found for complexes representing almost all the known protein–protein interactions, provided the components themselves have a known structure or can be homology-built. About one-third of the templates are of good quality when they are compared to experimental structures in test sets extracted from the Protein Data Bank and would be useful starting points in modeling the complexes. This finding dramatically expands our ability to model protein interactions, and has far-reaching implications for the protein docking field in general.

Specificity of Linear Motifs that Bind to a Common Mitogen-Activated Protein Kinase Docking Groove.

Abstract: Mitogen-activated protein kinases (MAPKs) have a docking groove that interacts with linear "docking" motifs in binding partners. To determine the structural basis of binding specificity between MAPKs and docking motifs, we quantitatively analyzed the ability of 15 docking motifs from diverse MAPK partners to bind to c-Jun amino-terminal kinase 1 (JNK1), p38α, and extracellular signal-regulated kinase 2 (ERK2). Classical docking motifs mediated highly specific binding only to JNK1, and only those motifs with a sequence pattern distinct from the classical MAPK binding docking motif consensus differentiated between the topographically similar docking grooves of ERK and p38α. Crystal structures of four complexes of MAPKs with docking peptides, representing JNK-specific, ERK-specific, or ERK- and p38-selective binding modes, revealed that the regions located between consensus positions in the docking motifs showed conformational diversity. Although the consensus positions in the docking motifs served as anchor points that bound to common MAPK surface features and mostly contributed to docking in a nondiscriminatory fashion, the conformation of the intervening region between the anchor points mostly determined specificity. We designed peptides with tailored MAPK binding profiles by rationally changing the length and amino acid composition of intervening regions located between anchor points. These results suggest a coherent structural model for MAPK docking specificity that reveals how short linear motifs binding to a common kinase docking groove can mediate diverse interaction patterns and contribute to correct MAPK partner selection in signaling networks.

idTarget: a web server for identifying protein targets of small chemical molecules with robust scoring functions and a divide-and-conquer docking approach

Abstract: Identification of possible protein targets of small chemical molecules is an important step for unravelling their underlying causes of actions at the molecular level. To this end, we construct a web server, idTarget, which can predict possible binding targets of a small chemical molecule via a divide-and-conquer docking approach, in combination with our recently developed scoring functions based on robust regression analysis and quantum chemical charge models. Affinity profiles of the protein targets are used to provide the confidence levels of prediction. The divide-and-conquer docking approach uses adaptively constructed small overlapping grids to constrain the searching space, thereby achieving better docking efficiency. Unlike previous approaches that screen against a specific class of targets or a limited number of targets, idTarget screen against nearly all protein structures deposited in the Protein Data Bank (PDB). We show that idTarget is able to reproduce known off-targets of drugs or drug-like compounds, and the suggested new targets could be prioritized for further investigation. idTarget is freely available as a web-based server at http://idtarget.rcas.sinica.edu.tw.

Structure-based ligand discovery for the protein-protein interface of chemokine receptor CXCR4.

Abstract: G-protein–coupled receptors (GPCRs) are key signaling molecules and are intensely studied. Whereas GPCRs recognizing small-molecules have been successfully targeted for drug discovery, protein-recognizing GPCRs, such as the chemokine receptors, claim few drugs or even useful small molecule reagents. This reflects both the difficulties that attend protein–protein interface inhibitor discovery, and the lack of structures for these targets. Imminent structure determination of chemokine receptor CXCR4 motivated docking screens for new ligands against a homology model and subsequently the crystal structure. More than 3 million molecules were docked against the model and then against the crystal structure; 24 and 23 high-scoring compounds from the respective screens were tested experimentally. Docking against the model yielded only one antagonist, which resembled known ligands and lacked specificity, whereas the crystal structure docking yielded four that were dissimilar to previously known scaffolds and apparently specific. Intriguingly, several were potent and relatively small, with IC50 values as low as 306 nM, ligand efficiencies as high as 0.36, and with efficacy in cellular chemotaxis. The potency and efficiency of these molecules has few precedents among protein–protein interface inhibitors, and supports structure-based efforts to discover leads for chemokine GPCRs.

Characterization of the Modes of Binding between Human Sweet Taste Receptor and Low-Molecular-Weight Sweet Compounds

Abstract: One of the most distinctive features of human sweet taste perception is its broad tuning to chemically diverse compounds ranging from low-molecular-weight sweeteners to sweet-tasting proteins. Many reports suggest that the human sweet taste receptor (hT1R2–hT1R3), a heteromeric complex composed of T1R2 and T1R3 subunits belonging to the class C G protein–coupled receptor family, has multiple binding sites for these sweeteners. However, it remains unclear how the same receptor recognizes such diverse structures. Here we aim to characterize the modes of binding between hT1R2–hT1R3 and low-molecular-weight sweet compounds by functional analysis of a series of site-directed mutants and by molecular modeling–based docking simulation at the binding pocket formed on the large extracellular amino-terminal domain (ATD) of hT1R2. We successfully determined the amino acid residues responsible for binding to sweeteners in the cleft of hT1R2 ATD. Our results suggest that individual ligands have sets of specific residues for binding in correspondence with the chemical structures and other residues responsible for interacting with multiple ligands.

Targeting Aquaporin Function: Potent Inhibition of Aquaglyceroporin-3 by a Gold-Based Compound

Abstract: Aquaporins (AQPs) are membrane channels that conduct water and small solutes such as glycerol and are involved in many physiological functions. Aquaporin-based modulator drugs are predicted to be of broad potential utility in the treatment of several diseases. Until today few AQP inhibitors have been described as suitable candidates for clinical development. Here we report on the potent inhibition of AQP3 channels by gold(III) complexes screened on human red blood cells (hRBC) and AQP3-transfected PC12 cells by a stopped-flow method. Among the various metal compounds tested, Auphen is the most active on AQP3 (IC50 = 0.8±0.08 µM in hRBC). Interestingly, the compound poorly affects the water permeability of AQP1. The mechanism of gold inhibition is related to the ability of Au(III) to interact with sulphydryls groups of proteins such as the thiolates of cysteine residues. Additional DFT and modeling studies on possible gold compound/AQP adducts provide a tentative description of the system at a molecular level. The mapping of the periplasmic surface of an homology model of human AQP3 evidenced the thiol group of Cys40 as a likely candidate for binding to gold(III) complexes. Moreover, the investigation of non-covalent binding of Au complexes by docking approaches revealed their preferential binding to AQP3 with respect to AQP1. The high selectivity and low concentration dependent inhibitory effect of Auphen (in the nanomolar range) together with its high water solubility makes the compound a suitable drug lead for future in vivo studies. These results may present novel metal-based scaffolds for AQP drug development.

Characterization of the interaction of interleukin-8 with hyaluronan, chondroitin sulfate, dermatan sulfate and their sulfated derivatives by spectroscopy and molecular modeling.

Abstract: The interactions between glycosaminoglycans (GAGs), important components of the extracellular matrix, and proteins such as growth factors and chemokines play critical roles in cellular regulation processes. Therefore, the design of GAG derivatives for the development of innovative materials with bio-like properties in terms of their interaction with regulatory proteins is of great interest for tissue engineering and regenerative medicine. Previous work on the chemokine interleukin-8 (IL-8) has focused on its interaction with heparin and heparan sulfate, which regulate chemokine function. However, the extracellular matrix contains other GAGs, such as hyaluronic acid (HA), dermatan sulfate (DS) and chondroitin sulfate (CS), which have so far not been characterized in terms of their distinct molecular recognition properties towards IL-8 in relation to their length and sulfation patterns. NMR and molecular modeling have been in great part the methods of choice to study the structural and recognition properties of GAGs and their protein complexes. However, separately these methods have challenges to cope with the high degree of similarity and flexibility that GAGs exhibit. In this work, we combine fluorescence spectroscopy, NMR experiments, docking and molecular dynamics simulations to study the configurational and recognition properties of IL-8 towards a series of HA and CS derivatives and DS. We analyze the effects of GAG length and sulfation patterns in binding strength and specificity, and the influence of GAG binding on IL-8 dimer formation. Our results highlight the importance of combining experimental and theoretical approaches to obtain a better understanding of the molecular recognition properties of GAG-protein systems.

Sampling protein motion and solvent effect during ligand binding

Abstract: An exhaustive description of the molecular recognition mechanism between a ligand and its biological target is of great value because it provides the opportunity for an exogenous control of the related process. Very often this aim can be pursued using high resolution structures of the complex in combination with inexpensive computational protocols such as docking algorithms. Unfortunately, in many other cases a number of factors, like protein flexibility or solvent effects, increase the degree of complexity of ligand/protein interaction and these standard techniques are no longer sufficient to describe the binding event. We have experienced and tested these limits in the present study in which we have developed and revealed the mechanism of binding of a new series of potent inhibitors of Adenosine Deaminase. We have first performed a large number of docking calculations, which unfortunately failed to yield reliable results due to the dynamical character of the enzyme and the complex role of the solvent. Thus, we have stepped up the computational strategy using a protocol based on metadynamics. Our approach has allowed dealing with protein motion and solvation during ligand binding and finally identifying the lowest energy binding modes of the most potent compound of the series, 4-decyl-pyrazolo[1,5-a]pyrimidin-7-one.

Inhibition Mechanism and Model of an Angiotensin I-Converting Enzyme (ACE)-Inhibitory Hexapeptide from Yeast (Saccharomyces cerevisiae)

Abstract: Angiotensin I-converting enzyme (ACE) has an important function in blood pressure regulation. ACE-inhibitory peptides can lower blood pressure by inhibiting ACE activity. Based on the sequence of an ACE-inhibitory hexapeptide (TPTQQS) purified from yeast, enzyme kinetics experiments, isothermal titration calorimetry (ITC), and a docking simulation were performed. The hexapeptide was found to inhibit ACE in a non-competitive manner, as supported by the structural model. The hexapeptide bound to ACE via interactions of the N-terminal Thr1, Thr3, and Gln4 residues with the residues on the lid structure of ACE, and the C-terminal Ser6 attracted the zinc ion, which is vital for ACE catalysis. The displacement of the zinc ion from the active site resulted in the inhibition of ACE activity. The structural model based on the docking simulation was supported by experiments in which the peptide was modified. This study provides a new inhibitory mechanism of ACE by a peptide which broads our knowledge for drug designing against enzyme targets.

Structure–Activity Relationships, Ligand Efficiency, and Lipophilic Efficiency Profiles of Benzophenone-Type Inhibitors of the Multidrug Transporter P-Glycoprotein

Abstract: The drug efflux pump P-glycoprotein (P-gp) has been shown to promote multidrug resistance (MDR) in tumors as well as to influence ADME properties of drug candidates. Here we synthesized and tested a series of benzophenone derivatives structurally analogous to propafenone-type inhibitors of P-gp. Some of the compounds showed ligand efficiency and lipophilic efficiency (LipE) values in the range of compounds which entered clinical trials as MDR modulators. Interestingly, although lipophilicity plays a dominant role for P-gp inhibitors, all compounds investigated showed LipE values below the threshold for promising drug candidates. Docking studies of selected analogues into a homology model of P-glycoprotein suggest that benzophenones show an interaction pattern similar to that previously identified for propafenone-type inhibitors.

Homology modeling, molecular dynamics, e-pharmacophore mapping and docking study of Chikungunya virus nsP2 protease

Abstract: To date, no suitable vaccine or specific antiviral drug is available to treat Chikungunya viral (CHIKV) fever. Hence, it is essential to identify drug candidates that could potentially impede CHIKV infection. Here, we present the development of a homology model of nsP2 protein based on the crystal structure of the nsP2 protein of Venezuelan equine encephalitis virus (VEEV). The protein modeled was optimized using molecular dynamics simulation; the junction peptides of a nonstructural protein complex were then docked in order to investigate the possible protein-protein interactions between nsP2 and the proteins cleaved by nsP2. The modeling studies conducted shed light on the binding modes, and the critical interactions with the peptides provide insight into the chemical features needed to inhibit the CHIK virus infection. Energy-optimized pharmacophore mapping was performed using the junction peptides. Based on the results, we propose the pharmacophore features that must be present in an inhibitor of nsP2 protease. The resulting pharmacophore model contained an aromatic ring, a hydrophobic and three hydrogen-bond donor sites. Using these pharmacophore features, we screened a large public library of compounds (Asinex, Maybridge, TOSLab, Binding Database) to find a potential ligand that could inhibit the nsP2 protein. The compounds that yielded a fitness score of more than 1.0 were further subjected to Glide HTVS and Glide XP. Here, we report the best four compounds based on their docking scores; these compounds have IDs of 27943, 21362, ASN 01107557 and ASN 01541696. We propose that these compounds could bind to the active site of nsP2 protease and inhibit this enzyme. Furthermore, the backbone structural scaffolds of these four lead compounds could serve as building blocks when designing drug-like molecules for the treatment of Chikungunya viral fever.

Catalytic transitions in the human MDR1 P-glycoprotein drug binding sites.

Abstract: Multidrug resistance proteins that belong to the ATP-binding cassette family like the human P-glycoprotein (ABCB1 or Pgp) are responsible for many failed cancer and antiviral chemotherapies because these membrane transporters remove the chemotherapeutics from the targeted cells. Understanding the details of the catalytic mechanism of Pgp is therefore critical to the development of inhibitors that might overcome these resistances. In this work, targeted molecular dynamics techniques were used to elucidate catalytically relevant structures of Pgp. Crystal structures of homologues in four different conformations were used as intermediate targets in the dynamics simulations. Transitions from conformations that were wide open to the cytoplasm to transition state conformations that were wide open to the extracellular space were studied. Twenty-six nonredundant transitional protein structures were identified from these targeted molecular dynamics simulations using evolutionary structure analyses. Coupled movement of nucleotide binding domains (NBDs) and transmembrane domains (TMDs) that form the drug binding cavities were observed. Pronounced twisting of the NBDs as they approached each other as well as the quantification of a dramatic opening of the TMDs to the extracellular space as the ATP hydrolysis transition state was reached were observed. Docking interactions of 21 known transport ligands or inhibitors were analyzed with each of the 26 transitional structures. Many of the docking results obtained here were validated by previously published biochemical determinations. As the ATP hydrolysis transition state was approached, drug docking in the extracellular half of the transmembrane domains seemed to be destabilized as transport ligand exit gates opened to the extracellular space.

Curcumin Binds to Aβ1–40 Peptides and Fibrils Stronger Than Ibuprofen and Naproxen

Abstract: Binding of curcumin, naproxen, and ibuprofen to Aβ1–40 peptide and its fibrils is studied by docking method and all-atom molecular dynamics simulations. The Gromos96 43a1 force field and simple point charge model of water have been used for molecular dynamics simulations. It is shown that if the receptor is a monomer then naproxen and ibuprofen are bound to the same place that is different from the binding position of curcumin. However all of three ligands have the same binding pocket in fibrillar structures. The binding mechanism is studied in detail showing that the van der Waals interaction between ligand and receptor dominates over the electrostatic interaction. The binding free energies obtained by the molecular mechanic-Poisson–Boltzmann surface area method indicate that curcumin displays higher binding affinity than nonsteroidal anti-inflammatory drugs. Our results are in good agreement with the experiments.

Current Assessment of Docking into GPCR Crystal Structures and Homology Models: Successes, Challenges, and Guidelines

Abstract: The growing availability of novel structures for several G protein-coupled receptors (GPCRs) has provided new opportunities for structure-based drug design of ligands against this important class of targets. Here, we report a systematic analysis of the accuracy of docking small molecules into GPCR structures and homology models using both rigid receptor (Glide SP and Glide XP) and flexible receptor (Induced Fit Docking; IFD) methods. The ability to dock ligands into different structures of the same target (cross-docking) is evaluated for both agonist and inverse agonist structures of the A2A receptor and the β1- and β2-adrenergic receptors. In addition, we have produced homology models for the β1-adrenergic, β2-adrenergic, D3 dopamine, H1 histamine, M2 muscarine, M3 muscarine, A2A adenosine, S1P1, κ-opioid, and C-X-C chemokine 4 receptors using multiple templates and investigated the ability of docking to predict the binding mode of ligands in these models. Clear correlations are observed between the dock...

Protein docking prediction using predicted protein-protein interface

Abstract: Background Many important cellular processes are carried out by protein complexes. To provide physical pictures of interacting proteins, many computational protein-protein prediction methods have been developed in the past. However, it is still difficult to identify the correct docking complex structure within top ranks among alternative conformations.

New insights for drug design from the X-ray crystallographic structures of G-protein-coupled receptors.

Abstract: Methodological advances in X-ray crystallography have made possible the recent solution of X-ray structures of pharmaceutically important G protein-coupled receptors (GPCRs), including receptors for biogenic amines, peptides, a nucleoside, and a sphingolipid. These high-resolution structures have greatly increased our understanding of ligand recognition and receptor activation. Conformational changes associated with activation common to several receptors entail outward movements of the intracellular side of transmembrane helix 6 (TM6) and movements of TM5 toward TM6. Movements associated with specific agonists or receptors have also been described [e.g., extracellular loop (EL) 3 in the A2A adenosine receptor]. The binding sites of different receptors partly overlap but differ significantly in ligand orientation, depth, and breadth of contact areas in TM regions and the involvement of the ELs. A current challenge is how to use this structural information for the rational design of novel potent and selective ligands. For example, new chemotypes were discovered as antagonists of various GPCRs by subjecting chemical libraries to in silico docking in the X-ray structures. The vast majority of GPCR structures and their ligand complexes are still unsolved, and no structures are known outside of family A GPCRs. Molecular modeling, informed by supporting information from site-directed mutagenesis and structure-activity relationships, has been validated as a useful tool to extend structural insights to related GPCRs and to analyze docking of other ligands in already crystallized GPCRs.

Multi-LZerD: multiple protein docking for asymmetric complexes.

Abstract: The tertiary structures of protein complexes provide a crucial insight about the molecular mechanisms that regulate their functions and assembly. However, solving protein complex structures by experimental methods is often more difficult than single protein structures. Here, we have developed a novel computational multiple protein docking algorithm, Multi-LZerD, that builds models of multimeric complexes by effectively reusing pairwise docking predictions of component proteins. A genetic algorithm is applied to explore the conformational space followed by a structure refinement procedure. Benchmark on eleven hetero-multimeric complexes resulted in near-native conformations for all but one of them (a root mean square deviation smaller than 2.5A). We also show that our method copes with unbound docking cases well, outperforming the methodology that can be directly compared with our approach. Multi-LZerD was able to predict near-native structures for multimeric complexes of various topologies.

Mapping of discontinuous conformational epitopes by amide hydrogen/deuterium exchange mass spectrometry and computational docking

Abstract: Understanding antigen-antibody interactions at the sub-molecular level is of particular interest for scientific, regulatory, and intellectual property reasons, especially with increasing demand for monoclonal antibody therapeutic agents. Although various techniques are available for the determination of an epitope, there is no widely applicable, high-resolution, and reliable method available. Here, a combination approach using amide hydrogen/deuterium exchange coupled with proteolysis and mass spectrometry (HDX-MS) and computational docking was applied to investigate antigen-antibody interactions. HDX-MS is a widely applicable, medium-resolution, medium-throughput technology that can be applied to epitope identification. First, the epitopes of cytochrome c-E8, IL-13-CNTO607, and IL-17A-CAT-2200 interactions identified using the HDX-MS method were compared with those identified by X-ray co-crystal structures. The identified epitopes are in good agreement with those identified using high-resolution X-ray crystallography. Second, the HDX-MS data were used as constraints for computational docking. More specifically, the non-epitope residues of an antigen identified using HDX-MS were designated as binding ineligible during computational docking. This approach, termed HDX-DOCK, gave more tightly clustered docking poses than stand-alone docking for all antigen-antibody interactions examined and improved docking results significantly for the cytochrome c-E8 interaction.

An Integrated Study of Tyrosinase Inhibition by Rutin: Progress using a Computational Simulation

Abstract: Tyrosinase inhibition studies have recently gained the attention of researchers due to their potential application values We simulated docking (binding energies for AutoDock Vina: -91 kcal/mol) and performed a molecular dynamics simulation to verify docking results between tyrosinase and rutin The docking results suggest that rutin mostly interacts with histidine residues located in the active site A 10 ns molecular dynamics simulation showed that one copper ion at the tyrosinase active site was responsible for the interaction with rutin Kinetic analyses showed that rutin-mediated inactivation followed a first-order reaction and mono- and biphasic rate constants occurred with rutin The inhibition was a typical competitive type with K(i) = 110±025 mM Measurements of intrinsic and ANS-binding fluorescences showed that rutin showed a relatively strong binding affinity for tyrosinase and one possible binding site that could be a copper was detected accompanying with a hydrophobic exposure of tyrosinase Cell viability testing with rutin in HaCaT keratinocytes showed that no toxic effects were produced Taken together, rutin has the potential to be a potent anti-pigment agent The strategy of predicting tyrosinase inhibition based on hydroxyl group number and computational simulation may prove useful for the screening of potential tyrosinase inhibitors

Consensus Induced Fit Docking (cIFD): methodology, validation, and application to the discovery of novel Crm1 inhibitors

Abstract: We present the Consensus Induced Fit Docking (cIFD) approach for adapting a protein binding site to accommodate multiple diverse ligands for virtual screening. This novel approach results in a single binding site structure that can bind diverse chemotypes and is thus highly useful for efficient structure-based virtual screening. We first describe the cIFD method and its validation on three targets that were previously shown to be challenging for docking programs (COX-2, estrogen receptor, and HIV reverse transcriptase). We then demonstrate the application of cIFD to the challenging discovery of irreversible Crm1 inhibitors. We report the identification of 33 novel Crm1 inhibitors, which resulted from the testing of 402 purchased compounds selected from a screening set containing 261,680 compounds. This corresponds to a hit rate of 8.2 %. The novel Crm1 inhibitors reveal diverse chemical structures, validating the utility of the cIFD method in a real-world drug discovery project. This approach offers a pragmatic way to implicitly account for protein flexibility without the additional computational costs of ensemble docking or including full protein flexibility during virtual screening.

Homology Modeling of Dopamine D2 and D3 Receptors: Molecular Dynamics Refinement and Docking Evaluation

Abstract: Dopamine (DA) receptors, a class of G-protein coupled receptors (GPCRs), have been targeted for drug development for the treatment of neurological, psychiatric and ocular disorders. The lack of structural information about GPCRs and their ligand complexes has prompted the development of homology models of these proteins aimed at structure-based drug design. Crystal structure of human dopamine D3 (hD3) receptor has been recently solved. Based on the hD3 receptor crystal structure we generated dopamine D2 and D3 receptor models and refined them with molecular dynamics (MD) protocol. Refined structures, obtained from the MD simulations in membrane environment, were subsequently used in molecular docking studies in order to investigate potential sites of interaction. The structure of hD3 and hD2L receptors was differentiated by means of MD simulations and D3 selective ligands were discriminated, in terms of binding energy, by docking calculation. Robust correlation of computed and experimental Ki was obtained for hD3 and hD2L receptor ligands. In conclusion, the present computational approach seems suitable to build and refine structure models of homologous dopamine receptors that may be of value for structure-based drug discovery of selective dopaminergic ligands.