Improved protein–ligand docking using GOLD

ABSTRACT:We previously reported on the synthesis of substituted phenyl-4-hydroxy-1-piperidyl indole analogues with nanomolar affinity at D2 dopamine receptors, ranging from 10-to 100-fold selective for D2 compared to the D3 dopamine receptor subtype. More recently, we evaluated a panel of aripiprazole analogues, identifying several analogues that also exhibit D2 vs D3 dopamine receptor binding selectivity. These studies further characterize the intrinsic efficacy of the compound with the greatest binding selectivity from each chemical class, 1-((5-methoxy-1H-indol-3-yl)methyl)-4-(4-(methylthio)phenyl)piperidin-4-ol (SV 293) and 7-(4-(4-(2-methoxyphenyl)piperazin-1-yl)butoxy)-3,4-dihydroquinolin-2(1H)-one (SV-III-130s), using an adenylyl cyclase inhibition assay, a G-protein-coupled inward-rectifying potassium (GIRK) channel activation assay, and a cell based phospho-MAPK (pERK1/2) assay. SV 293 was found to be a neutral antagonist at D2 dopamine receptors using all three assays. SV-III-130s is a partial agonist using an adenylyl cyclase inhibition assay but an antagonist in the GIRK and phospho ERK1/2 assays. To define the molecular basis for the binding selectivity, the affinity of these two compounds was evaluated using (a) wild type human D2 and D3 receptors and (b) a panel of chimeric D2/D3 dopamine receptors. Computer-assisted modeling techniques were used to dock these compounds to the human D2 and D3 dopamine receptor subtypes. It is hoped that these studies on D2 receptor selective ligands will be useful in the future design of (a) receptor selective ligands used to define the function of D2-like receptor subtypes, (b) novel pharmacotherapeutic agents, and/or (c) in vitro and in vivo imaging agents. KEYWORDS: Dopamine receptors, binding selectivity, functional selectivity, GPCR structure, D2-like dopamine receptors, GPCR model building, ligand−receptor docking T here are three dopaminergic pathways in the brain: the nigrostriatal pathway, the mesocorticolimbic pathway, and the tuberoinfundibular pathway. These pathways are involved in movement coordination, cognition, emotion, memory, reward, and regulation of prolactin secretion. Alterations in the dopaminergic pathways are thought to be involved in the pathogenesis of neurological, neuropsychiatric, and hormonal disorders. 1−6 Modulation of the dopaminergic pathways is also thought to occur as a consequence of acute or chronic abuse of pyschostimulants. 7,8 Previous studies have defined two types of dopamine receptors, the D1-like (D1 and D5 subtypes) and D2-like (D2, D3, and D4 subtypes) receptors. D1-like receptors are linked to the activation of adenylyl cyclase via coupling to the Gs/Golf class of G proteins. 9 Stimulation of the D2-like receptors results in coupling with the Gi/Go class of G proteins, leading to the inhibition of adenylyl cyclase activity. 10,11 Agonist activation of D2-like receptors can also lead to (a) activation of G-protein-coupled inward rectifying potassium (GIRK) channels, (b) stimulation of mitogenesis, (c) an increase in...

show abstract

“…57 For all the ligands, the piperidine rings were treated as rigid in a chair conformation. All the other rotatable bonds were set as flexible.…”

Section: ■ Methodsmentioning

confidence: 99%

Comparison of the Binding and Functional Properties of Two Structurally Different D2 Dopamine Receptor Subtype Selective Compounds

Luedtke

Murti

Wang

et al. 2012

ACS Chem. Neurosci.

View full text Add to dashboard Cite

show abstract

“…These calculations, which considered the variations in the distances between all lysine residues during the transition from Lf to LfFe 2 , were performed using an additional program written in house in Python and implemented in the UCSF Chimera package [32]. The script provided a 46 x 46 matrix of the entire set of Lys-Lys distances.…”

Section: Fig 5 Lowest-energy Solutions For the Conformation Of Af430mentioning

confidence: 99%

Monitoring lactoferrin iron levels by fluorescence resonance energy transfer: a combined chemical and computational study

et al. 2014

View full text Add to dashboard Cite

Abstract.Three forms of lactoferrin (Lf) that differed in their levels of iron loading (Lf, LfFe, and LfFe 2 ) were simultaneously labeled with the fluorophores AF350 and AF430. All three resulting fluorescent lactoferrins exhibited fluorescence resonance energy transfer (FRET), but they all presented different FRET patterns. Whereas only partial FRET was observed for Lf and LfFe, practically complete FRET was seen for the holo form (LfFe 2 ). For each form of metal-loaded lactoferrin, the AF350-AF430 distance varied depending on the protein conformation, which in turn depended on the level of iron loading. Thus, the FRET patterns of these lactoferrins were found to correlate with their iron loading levels. In order to gain greater insight into the number of fluorophores and the different FRET patterns observed (i.e., their iron levels), a computational analysis was performed. The results highlighted a number of lysines that have the greatest influence on the FRET profile. Moreover, despite the lack of an X-ray structure for any LfFe species, our study also showed that this species presents modified subdomain organization of the N-lobe, which narrows its iron-binding site. Complete domain rearrangement occurs during the LfFe to LfFe 2 transition. This latter study demonstrated that lactoferrin supplies iron to this bacterium, and suggested that this process occurs with no protein internalization.

show abstract

“…[53] 2.3 Molecular docking Selection of dienophiles and generation of initial studies were accomplished by docking simulations, using the GOLD Suite. [54,55,56,57] The docking procedure was conducted according to the following protocol: The built-in genetic algorithm (GA) was used to generate different poses, which were evaluated with the ChemScore scoring function. [58,59] The space of possible poses was limited to relevant ones by inducing hydrogen bond constraints between the dienophile carbonyl oxygen and T40-NH, T40-O γ H and Q106-NH.…”

Section: Enzyme Variantsmentioning

confidence: 99%

Computational design of a lipase for catalysis of the Diels-Alder reaction

et al. 2010

View full text Add to dashboard Cite

Combined molecular docking, molecular dynamics (MD) and density functional theory (DFT) studies have been employed to study catalysis of the Diels-Alder reaction by a modified lipase. Six variants of the versatile enzyme {\em Candida Antarctica} lipase B (CALB) have been rationally engineered {\em in silico} based on the specific characteristics of the pericyclic addition. A kinetic analysis reveals that hydrogen bond stabilization of the transition state and substrate binding are key components of the catalytic process. In the case of substrate binding, which has the greater potential for optimization, both binding strength and positioning of the substrates are important for catalytic efficiency. The binding strength is determined by hydrophobic interactions and can be tuned by careful selection of solvent and substrates. The MD simulations show that substrate positioning is sensitive to cavity shape and size, and can be controlled by a few rational mutations. The well-documented S105A mutation is essential to enable sufficient space in the vicinity of the oxyanion hole. Moreover, bulky residues on the edge of the active site hinders the formation of a sandwich-like near-attack conformer (NAC), and the I189A mutation is needed to obtain enough space above the face of the $\alpha$,$\beta$-double bond on the dienophile. The double mutant S105A/I189A performs quite well for two of three dienophiles. Based on binding constants and NAC energies obtained from MD simulations combined with activation energies from DFT computations, relative catalytic rates ($v_{cat}/v_{uncat}$) of up to $10^3$ are predicted.Response to Reviewers: We are happy with the favorable comments by the reviewers. We have have corrected the manuscript as suggested by reviewer 1 1. A sentence explaining catalytic promiscuity has been added. 2. A paragraph discussing the relationship between the NAC concept and the Reaction Force has been added.We have shortened the manuscript as suggested by reviewer 2. In particular, we have moved information regarding the structural relaxation of protein structure in MD simulations to the supplementary material. Abstract Combined molecular docking, molecular dynamics (MD) and density functional theory (DFT) studies have been employed to study catalysis of the DielsAlder reaction by a modified lipase. Six variants of the versatile enzyme Candida Antarctica lipase B (CALB) have been rationally engineered in silico based on the specific characteristics of the pericyclic addition. A kinetic analysis reveals that hydrogen bond stabilization of the transition state and substrate binding are key components of the catalytic process. In the case of substrate binding, which has the greater potential for optimization, both binding strength and positioning of the substrates are important for catalytic efficiency. The binding strength is determined by hydrophobic interactions and can be tuned by careful selection of solvent and substrates. The MD simulations show that substrate positioning is sensitive to cavity s...

show abstract

Improved protein–ligand docking using GOLD

Cited by 2,521 publications

References 30 publications

Comparison of the Binding and Functional Properties of Two Structurally Different D2 Dopamine Receptor Subtype Selective Compounds

Comparison of the Binding and Functional Properties of Two Structurally Different D2 Dopamine Receptor Subtype Selective Compounds

Monitoring lactoferrin iron levels by fluorescence resonance energy transfer: a combined chemical and computational study

Computational design of a lipase for catalysis of the Diels-Alder reaction

Contact Info

Product

Resources

About