Computational modeling toward understanding agonist binding on dopamine 3

In Silico Repositioning of Dopamine Modulators with Possible Application to Schizophrenia: Pharmacophore Mapping, Molecular Docking and Molecular Dynamics Analysis

We have performed theoretical calculations with 70 drugs that have been considered in 231 clinical trials as possible candidates to repurpose drugs for schizophrenia based on their interactions with the dopaminergic system. A hypothesis of shared pharmacophore features was formulated to support our calculations. To do so, we have used the crystal structure of the D2-like dopamine receptor in complex with risperidone, eticlopride, and nemonapride. Linagliptin, citalopram, flunarizine, sildenafil, minocycline, and duloxetine were the drugs that best fit with our model. Molecular docking calculations, molecular dynamics outcomes, blood-brain barrier penetration, and human intestinal absorption were studied and compared with the results. From the six drugs selected in the shared pharmacophore features input, flunarizine showed the best docking score with D2, D3, and D4 dopamine receptors and had high stability during molecular dynamics simulations. Flunarizine is a frequently used medication to treat migraines and vertigo. However, its antipsychotic properties have been previously hypothesized, particularly because of its possible ability to block the D2 dopamine receptors.

Receptor Ligands as Helping Hands to L-DOPA in the Treatment of Parkinson's Disease

Levodopa (LD) is the most effective drug in the treatment of Parkinson’s disease (PD). However, although it represents the “gold standard” of PD therapy, LD can cause side effects, including gastrointestinal and cardiovascular symptoms as well as transient elevated liver enzyme levels. Moreover, LD therapy leads to LD-induced dyskinesia (LID), a disabling motor complication that represents a major challenge for the clinical neurologist. Due to the many limitations associated with LD therapeutic use, other dopaminergic and non-dopaminergic drugs are being developed to optimize the treatment response. This review focuses on recent investigations about non-dopaminergic central nervous system (CNS) receptor ligands that have been identified to have therapeutic potential for the treatment of motor and non-motor symptoms of PD. In a different way, such agents may contribute to extending LD response and/or ameliorate LD-induced side effects.

Computational Study and Modified Design of Selective Dopamine D3 Receptor Agonists

Dopamine D3 receptor (D3 R) is considered as a potential target for the treatment of nervous system disorders, such as Parkinson's disease. Current research interests primarily focus on the discovery and design of potent D3 agonists. In this work, we selected 40 D3 R agonists as the research system. Comparative molecular field analysis (CoMFA) of three-dimensional quantitative structure-activity relationship (3D-QSAR), structure-selectivity relationship (3D-QSSR), and molecular docking was performed on D3 receptor agonists to obtain the details at atomic level. The results indicated that both the CoMFA model (r(2) = 0.982, q(2) = 0.503, rpred2 = 0.893, SEE  = 0.057, F = 166.308) for structure-activity and (r(2) = 0.876, q(2) = 0.436, rpred2 = 0.828, F = 52.645) for structure-selectivity have good predictive capabilities. Furthermore, docking studies on three compounds binding to D3 receptor were performed to analyze the binding modes and interactions. The results elucidate that agonists formed hydrogen bond and hydrophobic interactions with key residues. Finally, we designed six molecules under the guidance of 3D-QSAR/QSSR models. The activity and selectivity of designed molecules have been improved, and ADMET properties demonstrate they have low probability of hepatotoxicity (<0.5). These results from 3D-QSAR/QSSR and docking studies have great significance for designing novel dopamine D3 selective agonists in the future.

Identification of key residues involved in the activation and signaling properties of dopamine D3 receptor

The dopamine D3 receptor exhibits agonist-dependent tolerance and slow response termination (SRT) signaling properties that distinguish it from the closely-related D2 receptors. While amino acid residues important for D3 receptor ligand binding have been identified, the residues involved in activation of D3 receptor signaling and induction of signaling properties have not been determined. In this paper, we used cis and trans isomers of a novel D3 receptor agonist, 8-OH-PBZI, and site-directed mutagenesis to identify key residues involved in D3 receptor signaling function. Our results show that trans-8-OH-PBZI, but not cis-8-OH-PBZI, elicit the D3 receptor tolerance and SRT properties. We show that while both agonists require a subset of residues in the orthosteric binding site of D3 receptors for activation of the receptor, the ability of the two isomers to differentially induce tolerance and SRT is mediated by interactions with specific residues in the sixth transmembrane helix and third extracellular loop of the D3 receptor. We also show that unlike cis-8-OH-PBZI, which is a partial agonist at the dopamine D2S receptor and full agonist at dopamine D2L receptor, trans-8-OH-PBZI is a full agonist at both D2S and D2L receptors. The different effect of the two isomers on D3 receptor signaling properties and D2S receptor activation correlated with differential effects of the isomers on agonist-induced mouse locomotor activity. The two isomers of 8-OH-PBZI represent novel pharmacological tools for in silico D3 and D2 receptor homology modeling and for determining the role of D3 receptor tolerance and SRT properties in signaling and behavior.

Molecular modeling study on the allosteric inhibition mechanism of HIV-1 integrase by LEDGF/p75 binding site inhibitors

HIV-1 integrase (IN) is essential for the integration of viral DNA into the host genome and an attractive therapeutic target for developing antiretroviral inhibitors. LEDGINs are a class of allosteric inhibitors targeting LEDGF/p75 binding site of HIV-1 IN. Yet, the detailed binding mode and allosteric inhibition mechanism of LEDGINs to HIV-1 IN is only partially understood, which hinders the structure-based design of more potent anti-HIV agents. A molecular modeling study combining molecular docking, molecular dynamics simulation, and binding free energy calculation were performed to investigate the interaction details of HIV-1 IN catalytic core domain (CCD) with two recently discovered LEDGINs BI-1001 and CX14442, as well as the LEDGF/p75 protein. Simulation results demonstrated the hydrophobic domain of BI-1001 and CX14442 engages one subunit of HIV-1 IN CCD dimer through hydrophobic interactions, and the hydrophilic group forms hydrogen bonds with HIV-1 IN CCD residues from other subunit. CX14442 has a larger tert-butyl group than the methyl of BI-1001, and forms better interactions with the highly hydrophobic binding pocket of HIV-1 IN CCD dimer interface, which can explain the stronger affinity of CX14442 than BI-1001. Analysis of the binding mode of LEDGF/p75 with HIV-1 IN CCD reveals that the LEDGF/p75 integrase binding domain residues Ile365, Asp366, Phe406 and Val408 have significant contributions to the binding of the LEDGF/p75 to HIV1-IN. Remarkably, we found that binding of BI-1001 and CX14442 to HIV-1 IN CCD induced the structural rearrangements of the 140 s loop and oration displacements of the side chains of the three conserved catalytic residues Asp64, Asp116, and Glu152 located at the active site. These results we obtained will be valuable not only for understanding the allosteric inhibition mechanism of LEDGINs but also for the rational design of allosteric inhibitors of HIV-1 IN targeting LEDGF/p75 binding site.

Structural modifications of neuroprotective anti-Parkinsonian (-)-N6-(2-(4-(biphenyl-4-yl)piperazin-1-yl)-ethyl)-N6-propyl-4,5,6,7-tetrahydrobenzo[d]thiazole-2,6-diamine (D-264): an effort toward the improvement of in vivo efficacy of the parent molecule

In our overall goal to develop multifunctional dopamine D2/D3 agonist drugs for the treatment of Parkinson's disease (PD), we previously synthesized potent D3 preferring agonist D-264 (1a), which exhibited neuroprotective properties in two animal models of PD. To enhance the in vivo efficacy of 1a, a structure-activity relationship study was carried out. Competitive binding and [(35)S]GTPγS functional assays identified compound (-)-9b as one of the lead molecules with preferential D3 agonist activity (EC50(GTPγS); D3 = 0.10 nM; D2/D3 (EC50): 159). Compounds (-)-9b and (-)-8b exhibited high in vivo activity in two PD animal models, reserpinized and 6-hydroxydopamine (OHDA)-induced unilateral lesioned rats. On the other hand, 1a failed to show any in vivo activity in these models unless the compound was dissolved in 5-10% beta-hydroxy propyl cyclodextrin solution. Lead compounds exhibited appreciable radical scavenging activity. In vitro experiments with dopaminergic MN9D cells indicated neuroprotection by both 1a and (-)-9b from toxicity of MPP+.

Understanding the Structural Requirements of Hybrid (S)-6-((2-(4-Phenylpiperazin-1-yl)ethyl)(propyl)amino)-5,6,7,8-tetrahydronaphthalen-1-ol and its Analogs as D2/D3 Receptor Ligands: A Three-Dimensional Quantitative Structure-Activity Relationship (3D QSAR) Investigation

To gain insights into the structural requirements for dopamine D2 and D3 agonists in the treatment of Parkinson's disease (PD) and to elucidate the basis of selectivity for D3 over D2 (D2/D3), 3D quantitative structure-activity relationship (3D QSAR) investigations using CoMFA (comparative molecular field analysis) and CoMSIA (comparative molecular similarity indices analysis) methods were performed on a series of 45 structurally related D2 and D3 dopaminergic ligands. Two alignment methods (atom-based and flexible) and two charge calculation methods (Gasteiger-Hückel and AM1) were used in the present study. Overall, D2 affinity and selectivity (D2/D3) models performed better with r²_cv values of 0.71 and 0.63 for CoMFA and 0.71 and 0.79 for CoMSIA, respectively. The corresponding predictive r² values for the CoMFA and CoMSIA models were 0.92 and 0.86 and 0.91 and 0.78, respectively. The CoMFA models generated using flexible alignment outperformed the models with the atom-based alignment in terms of relevant statistics and interpretability of the generated contour maps while CoMSIA models obtained using atom-based alignment showed superiority in terms of internal and external predictive abilities. The presence of carbonyl group (C=O) attached to the piperazine ring and the hydrophobic biphenyl ring were found to be the most important features responsible for the D3 selectivity over D2. This study can be further utilized to design and develop selective and potent dopamine agonists to treat PD.

Insights into the structural biology of G-protein coupled receptors impacts drug design for central nervous system neurodegenerative processes

In the last few years, there have been important new insights into the structural biology of G-protein coupled receptors. It is now known that allosteric binding sites are involved in the affinity and selectivity of ligands for G-protein coupled receptors, and that signaling by these receptors involves both G-protein dependent and independent pathways. The present review outlines the physiological and pharmacological implications of this perspective for the design of new drugs to treat disorders of the central nervous system. Specifically, new possibilities are explored in relation to allosteric and orthosteric binding sites on dopamine receptors for the treatment of Parkinson's disease, and on muscarinic receptors for Alzheimer's disease. Future research can seek to identify ligands that can bind to more than one site on the same receptor, or simultaneously bind to two receptors and form a dimer. For example, the design of bivalent drugs that can reach homo/hetero-dimers of D2 dopamine receptor holds promise as a relevant therapeutic strategy for Parkinson's disease. Regarding the treatment of Alzheimer's disease, the design of dualsteric ligands for mono-oligomeric rinic receptors could increase therapeutic effectiveness by generating potent compounds that could activate more than one signaling pathway.

Progress in the structural prediction of G protein-coupled receptors: D3 receptor in complex with eticlopride

Predicting the three-dimensional structure of ligand-receptor complexes involving G protein-coupled receptors (GPCRs) is still a challenging task in rational drug design. To evaluate the reliability of the GPCR structural prediction, only a couple of community-wide assessments have been carried out. Our participation in the last edition, DOCK2010, involved the blind prediction of the dopaminergic D(3) receptor in complex with the D(2)/D(3) selective antagonist eticlopride for which the crystal structure has been recently released. Here, we describe a methodology that succeeded to produce a correctly predicted eticlopride-D(3) receptor complex out of three submitted models. Ranking the obtained models in the correct order is the main challenge due to subtle structural differences in the complex that are not sufficiently captured by conventional scoring functions. Importantly, our work reveals that a correct ranking is obtained by including a more sophisticated description of conformational ligand energy on binding. All in all, this case study highlights the current progress in modeling GPCR complexes and underlines that in silico modeling can be a valuable complement in GPCR drug discovery.

Recent advances in the search for D3- and D4-selective drugs: probes, models and candidates

Dopamine D(2)-like receptors (including D(2), D(3) and D(4)) belong to the 'rhodopsin-like' family of G protein-coupled receptors (GPCRs), which represent the largest group of targets for bioactive molecules. Due to their high sequence similarity, the design of subtype-selective ligands requires rational and effective strategies. The general formula of 1,4-disubstituted aromatic piperidines and piperazines (1,4-DAPs) was extracted from classical dopaminergic drugs. The biological properties of this compound family are encoded by an aromatic head group that controls intrinsic activity, an amine moiety and a lipophilic appendage. D(3)- and D(4)-selective molecular probes and drug candidates were generated from the general formula of 1,4-DAP. Formal structural rearrangement led to investigational drugs beyond the 1,4-DAP structure. The very recent publication of the X-ray crystal structure of D(3) should facilitate efficient discovery of unprecedented chemotypes. However, the development of D(3)-selective agonists, functionally selective ligands and the exploitation of homo- and heteromers remain challenging.