The dopamine D(2) receptor dimer and its interaction with homobivalent antagonists: homology modeling, docking and molecular dynamics

In order to apply structure-based drug design techniques to G protein-coupled receptor complexes, it is essential to model their 3D structure and to identify regions that are suitable for selective drug binding. For this purpose, we have developed and tested a multi-component protocol to model the inactive conformation of the dopamine D(2) receptor dimer, suitable for interaction with homobivalent antagonists. Our approach was based on protein–protein docking, applying the Rosetta software to obtain populations of dimers as present in membranes with all the main possible interfaces. Consensus scoring based on the values and frequencies of best interfaces regarding four scoring parameters, Rosetta interface score, interface area, free energy of binding and energy of hydrogen bond interactions indicated that the best scored dimer model possesses a TM4–TM5–TM7–TM1 interface, which is in agreement with experimental data. This model was used to study interactions of the previously published dopamine D(2) receptor homobivalent antagonists based on clozapine,1,4-disubstituted aromatic piperidines/piperazines and arylamidoalkyl substituted phenylpiperazine pharmacophores. It was found that the homobivalent antagonists stabilize the receptor-inactive conformation by maintaining the ionic lock interaction, and change the dimer interface by disrupting a set of hydrogen bonds and maintaining water- and ligand-mediated hydrogen bonds in the extracellular and intracellular part of the interface. [Figure: see text] ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s00894-016-3065-2) contains supplementary material, which is available to authorized users.