Protein docking
Dotaz
Zobrazit nápovědu
G protein-coupled receptors (GPCRs) are hard to crystallize. However, attempts to predict their structure have boomed as a result of advancements in crystallographic techniques. This trend has allowed computer-aided molecular modeling of GPCRs. We analyzed the performance of four molecular modeling programs in pose evaluation of re-docked antagonists / inverse agonists to 11 original crystal structures of aminergic GPCRs using an induced fit-docking procedure. AutoDock and Glide were used for docking. AutoDock binding energy function, GlideXP, Prime MM-GB/SA, and YASARA binding function were used for pose scoring. Root mean square deviation (RMSD) of the best pose ranged from 0.09 to 1.58 Å, and median RMSD of the top 60 poses ranged from 1.47 to 3.83 Å. However, RMSD of the top pose ranged from 0.13 to 7.33 Å and ranking of the best pose ranged from the 1st to 60th out of 60 poses. Moreover, analysis of ligand-receptor interactions of top poses revealed substantial differences from interactions found in crystallographic structures. Bad ranking of top poses and discrepancies between top docked poses and crystal structures render current simple docking methods unsuitable for predictive modeling of receptor-ligand interactions. Prime MM-GB/SA optimized for 3NY9 by multiple linear regression did not work well at 3NY8 and 3NYA, structures of the same receptor with different ligands. However, 9 of 11 trajectories of molecular dynamics simulations by Desmond of top poses converged with trajectories of crystal structures. Key interactions were properly detected for all structures. This procedure also worked well for cross-docking of tested β2-adrenergic antagonists. Thus, this procedure represents a possible way to predict interactions of antagonists with aminergic GPCRs.
Bringing a new drug to the market is a lengthy, risky and expensive endeavor. Money spent on developing new drugs keeps going up each year, which is disproportional to the number of drugs brought to the market. Therefore, it is important to find ways to reduce costs and improve the effectiveness of drug discovery and development. Underlined by fast-paced developments in algorithms and processing power of modern hardware, computational methods have shown great potential in achieving this goal and molecular docking is an important tool in this toolbox. In this work, we briefly introduce the very basic principles of molecular docking and review some important contemporary challenges and developments in this field.
- Klíčová slova
- skórování, změna vazebné energie,
- MeSH
- algoritmy MeSH
- počítačová simulace MeSH
- racionální návrh léčiv MeSH
- simulace molekulového dockingu * metody trendy MeSH
- termodynamika MeSH
- vazba proteinů MeSH
- Publikační typ
- práce podpořená grantem MeSH
MOTIVATION: Protein tunnels and channels are key transport pathways that allow ligands to pass between proteins' external and internal environments. These functionally important structural features warrant detailed attention. It is difficult to study the ligand binding and unbinding processes experimentally, while molecular dynamics simulations can be time-consuming and computationally demanding. RESULTS: CaverDock is a new software tool for analysing the ligand passage through the biomolecules. The method uses the optimized docking algorithm of AutoDock Vina for ligand placement docking and implements a parallel heuristic algorithm to search the space of possible trajectories. The duration of the simulations takes from minutes to a few hours. Here we describe the implementation of the method and demonstrate CaverDock's usability by: (i) comparison of the results with other available tools, (ii) determination of the robustness with large ensembles of ligands and (iii) the analysis and comparison of the ligand trajectories in engineered tunnels. Thorough testing confirms that CaverDock is applicable for the fast analysis of ligand binding and unbinding in fundamental enzymology and protein engineering. AVAILABILITY AND IMPLEMENTATION: User guide and binaries for Ubuntu are freely available for non-commercial use at https://loschmidt.chemi.muni.cz/caverdock/. The web implementation is available at https://loschmidt.chemi.muni.cz/caverweb/. The source code is available upon request. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
A comparative study of interaction between chicken egg white lysozyme (Lyz) with two hexavalent chromate ions; chromate and dichromate; which are prevalently known for their toxicity, was investigated using different spectroscopic techniques along with a molecular docking study. Both steady-state and time-resolved studies revealed that the addition of chromate/dichromate is responsible for strong quenching of intrinsic fluorescence in Lyz and the quenching is caused by both static and dynamic quenching mechanisms. Different binding and thermodynamic parameters were also calculated at different temperatures from the intrinsic fluorescence of Lyz. The conformational change in Lyz and thermodynamic parameters obtained during the course of interaction with chromate/dichromate were well-supported by the molecular docking results.
Formation of transient complexes of cytochrome P450 (P450) with another protein of the endoplasmic reticulum membrane, cytochrome b5 (cyt b5), dictates the catalytic activities of several P450s. Therefore, we examined formation and binding modes of the complex of human P450 1A2 with cyt b5. Docking of soluble domains of these proteins was performed using an information-driven flexible docking approach implemented in HADDOCK. Stabilities of the five unique binding modes of the P450 1A2-cyt b5 complex yielded by HADDOCK were evaluated using explicit 10 ns molecular dynamics (MD) simulations in aqueous solution. Further, steered MD was used to compare the stability of the individual P450 1A2-cyt b5 binding modes. The best binding mode was characterized by a T-shaped mutual orientation of the porphyrin rings and a 10.7 Å distance between the two redox centers, thus satisfying the condition for a fast electron transfer. Mutagenesis studies and chemical cross-linking, which, in the absence of crystal structures, were previously used to deduce specific P450-cyt b5 interactions, indicated that the negatively charged convex surface of cyt b5 binds to the positively charged concave surface of P450. Our simulations further elaborate structural details of this interface, including nine ion pairs between R95, R100, R138, R362, K442, K455, and K465 side chains of P450 1A2 and E42, E43, E49, D65, D71, and heme propionates of cyt b5. The universal heme-centric system of internal coordinates was proposed to facilitate consistent classification of the orientation of the two porphyrins in any protein complex.
Networks of protein-protein interactions (PPI) constitute either stable or transient complexes in every cell. Most of the cellular complexes keep their function, and therefore stay similar, during evolution. The evolutionary constraints preserve most cellular functions via preservation of protein structures and interactions. The evolutionary conservation information is utilized in template-based approaches, like protein structure modeling or docking. Here we use the combination of the template-free docking method with conservation-based selection of the best docking model using our newly developed COZOID tool.We describe a step-by-step protocol for visual selection of docking models, based on their similarity to the original protein complex structure. Using the COZOID tool, we first analyze contact zones of the original complex structure and select contact amino acids for docking restraints. Then we model and dock the homologous proteins. Finally, we utilize different analytical modes of our COZOID tool to select the docking models most similar to the original complex structure.
We report on the synthesis, activity testing, docking, and quantum mechanical scoring of novel imidazo[1,2-c]pyrimidin-5(6H)-one scaffold for cyclin-dependent kinase 2 (CDK2) inhibition. A series of 26 compounds substituted with aromatic moieties at position 8 has been tested in in vitro enzyme assays and shown to inhibit CDK2. 2D structure-activity relationships have ascertained that small substituents at position 8 (up to the size of naphtyl or methoxyphenyl) generally lead to single-digit micromolar IC50 values, whereas bigger substituents (substituted biphenyls) decreased the compounds' activities. The binding modes of the compounds obtained using Glide docking have exhibited up to 2 hinge-region hydrogen bonds to CDK2 and differed in the orientation of the inhibitor core and the placement of the 8-substituents. Semiempirical quantum mechanics-based scoring identified probable favourable binding modes, which will serve for future structure-based design and synthetic optimization of substituents of the heterocyclic core. In summary, we have identified a novel core for CDK2 inhibition and will explore it further to increase the potencies of the compounds and also monitor selectivities against other protein kinases.
- MeSH
- cyklin-dependentní kinasa 2 antagonisté a inhibitory chemie MeSH
- inhibitory proteinkinas chemická syntéza chemie farmakologie MeSH
- kvantová teorie MeSH
- lidé MeSH
- pyrimidiny chemická syntéza chemie farmakologie MeSH
- racionální návrh léčiv MeSH
- simulace molekulového dockingu MeSH
- vazba proteinů MeSH
- vztahy mezi strukturou a aktivitou MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Primaquine is a traditional antimalarial drug with low parasitic resistance and generally good acceptance at higher doses, which has been used for over 60 years in malaria treatment. However, several limitations related to its hematotoxicity have been reported. It is believed that this toxicity comes from the hydroxylation of the C-5 and C-6 positions of its 8-aminoquinoline ring before binding to the molecular target: the quinone reductase II (NQO2) human protein. In this study we propose primaquine derivatives, with substitution at position C-6 of the 8-aminoquinoline ring, planned to have better binding to NQO2, compared to primaquine, but with a reduced toxicity related to the C-5 position being possible to be oxidized. On this sense the proposed analogues were suggested in order to reduce or inhibit hydroxylation and further oxidation to hemotoxic metabolites. Five C-6 substituted primaquine analogues were selected by de novo design and further submitted to docking and molecular dynamics simulations. Our results suggest that all analogues bind better to NQO2 than primaquine and may become better antimalarials. However, the analogues 3 and 4 are predicted to have a better activity/toxicity balance.
- MeSH
- chinonreduktasy antagonisté a inhibitory chemie MeSH
- inhibitory enzymů chemie MeSH
- katalytická doména MeSH
- lidé MeSH
- primachin analogy a deriváty chemie MeSH
- sekundární struktura proteinů MeSH
- simulace molekulární dynamiky MeSH
- simulace molekulového dockingu MeSH
- termodynamika MeSH
- vazba proteinů MeSH
- vodíková vazba MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Biogenic amines are common biologically active substances extended within the whole animal kingdom where they play vital roles as signal transducer as well as regulator of cell functions. One of these biogenic amines called octopamine (OA) is synthesized from tyramine (TA) by the catalysis of tyramine-β-hydroxylase (TβH) originated in the insect nervous system. Both TA and OA act as neurotransmitters, neurohormones and neuromodulators in the arthropod nervous system. Herein, the inhibitory activity of 1-arylimidazole-2(3H)-thiones (AITs) was tested on cloned Drosophila tyramine-β-hydroxylase (DmTβH) expressed in Bombyx mori strain. Radiolabelled (3)H-TA was used to analyze the activity of AITs exhibited inhibitory effects on DmTβH, whose ID50 values range from 0.02 to 2511nM where DmTβH was inhibited in a dose-dependent manner at pH 7.6 and 25°C during a 30min of incubation. To understand the catalytic role of the TβH, a three dimensional structure of the TβH from Drosophila melanogaster was constructed by homology modeling using the Phyre2 web server with 100% confidence. The modeled three-dimensional structure of TβH was used to perform the docking study with AITs. This may give more insights to precise design of inhibitors for TβH to control insect's population.
- MeSH
- Drosophila melanogaster enzymologie MeSH
- inhibitory enzymů chemie MeSH
- katalytická doména MeSH
- konformace proteinů, beta-řetězec MeSH
- oxygenasy se smíšenou funkcí antagonisté a inhibitory chemie MeSH
- proteiny Drosophily antagonisté a inhibitory chemie MeSH
- sekvence aminokyselin MeSH
- simulace molekulární dynamiky MeSH
- simulace molekulového dockingu MeSH
- vazba proteinů MeSH
- zvířata MeSH
- Check Tag
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
Targeting of epigenetic mechanisms, such as the hydroxymethylation of DNA, has been intensively studied, with respect to the treatment of many serious pathologies, including oncological disorders. Recent studies demonstrated that promising therapeutic strategies could potentially be based on the inhibition of the TET1 protein (ten-eleven translocation methylcytosine dioxygenase 1) by specific iron chelators. Therefore, in the present work, we prepared a series of pyrrolopyrrole derivatives with hydrazide (1) or hydrazone (2-6) iron-binding groups. As a result, we determined that the basic pyrrolo[3,2-b]pyrrole derivative 1 was a strong inhibitor of the TET1 protein (IC50 = 1.33 μM), supported by microscale thermophoresis and molecular docking. Pyrrolo[3,2-b]pyrroles 2-6, bearing substituted 2-hydroxybenzylidene moieties, displayed no significant inhibitory activity. In addition, in vitro studies demonstrated that derivative 1 exhibits potent anticancer activity and an exclusive mitochondrial localization, confirmed by Pearson's correlation coefficient of 0.92.