PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “molecular docking”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4Linked to original sources

3D QSAR studies of AChE inhibitors based on molecular docking scores and CoMFA.

Three-dimensional quantitative structure-activity relationship (3D QSAR) studies were performed on acetylcholinesterase (AChE) inhibitors, based on molecular docking scores obtained by using FlexX and FlexiDock and comparative molecular field analysis (CoMFA). The docking scores were used as molecular descriptors along with the steric and electrostatic field values of CoMFA, for partial least square (PLS) analysis. The high leave one out (LOO) cross-validated correlation coefficient (q(2)=0.714) reveals that the model is a useful tool for the prediction of test set as well as newly designed structures against AChE activity. The superimposed CoMFA models on the receptor site of AChE are guiding the design of potential inhibitory structures directed against AChE activity.

Acetylcholinesterase↗

Comparative binding energy analysis of haloalkane dehalogenase substrates: modelling of enzyme-substrate complexes by molecular docking and quantum mechanical calculations.

We evaluate the applicability of automated molecular docking techniques and quantum mechanical calculations to the construction of a set of structures of enzyme-substrate complexes for use in Comparative binding energy (COMBINE) analysis to obtain 3D structure-activity relationships. The data set studied consists of the complexes of eighteen substrates docked within the active site of haloalkane dehalogenase (DhlA) from Xanthobacter autotrophicus GJ10. The results of the COMBINE analysis are compared with previously reported data obtained for the same dataset from modelled complexes that were based on an experimentally determined structure of the DhlA-dichloroethane complex. The quality of fit and the internal predictive power of the two COMBINE models are comparable, but better external predictions are obtained with the new approach. Both models show a similar composition of the principal components. Small differences in the relative contributions that are assigned to important residues for explaining binding affinity differences can be directly linked to structural differences in the modelled enzyme-substrate complexes: (i) rotation of all substrates in the active site about their longitudinal axis, (ii) repositioning of the ring of epihalohydrines and the halogen substituents of 1,2-dihalopropanes, and (iii) altered conformation of the long-chain molecules (halobutanes and halohexanes). For external validation, both a novel substrate not included in the training series and two different mutant proteins were used. The results obtained can be useful in the future to guide the rational engineering of substrate specificity in DhlA and other related enzymes.

Binding Sites↗

Functional characterization of peanut serine/threonine/tyrosine protein kinase: molecular docking and inhibition kinetics with tyrosine kinase inhibitors.

Serine/threonine/tyrosine (STY) protein kinase from peanut is developmentally regulated and is induced by abiotic stresses. In addition, STY protein kinase activity is regulated by tyrosine phosphorylation. Kinetic mechanism of plant dual specificity protein kinases is not studied so far. Recombinant STY protein kinase occurs as a monomer in solution as shown by gel filtration chromatography. The relative phosphorylation rate of kinase against increasing enzyme concentrations follows a first-order kinetics indicating an intramolecular phosphorylation mechanism. Moreover, the active recombinant STY protein kinase could not transphosphorylate a kinase-deficient mutant of STY protein kinase. Molecular docking studies revealed that the tyrosine kinase inhibitors bind the protein kinase at the same region as ATP. STY protein kinase activity was inhibited by the tyrosine kinase inhibitors, and the inhibitor potency series against the recombinant STY protein kinase was tyrphostin > genistein > staurosporine. The inhibition constant (K(i)), and the IC(50) value of STY protein kinase for tyrosine kinase inhibitors with ATP and histone are discussed. All the inhibitors competed with ATP. Genistein was an uncompetitive inhibitor with histone, whereas staurosporine and tyrphostin were linear mixed type noncompetitive inhibitors with histone. Molecular docking and kinetic analysis revealed that Y148F mutant of the "ATP-binding loop" and Y297F mutant of the "activation loop" showed a dramatic increase in K(i) values for genistein and tyrphostin with respect to wild-type STY protein kinase. Data presented here provide the direct evidence on the mechanism of inhibition of plant protein kinases by tyrosine kinase inhibitors. This study also suggests that tyrosine kinase inhibitors may be useful in unraveling the plant tyrosine phosphorylation signaling cascades.

Adenosine Triphosphate↗

Influence of stereoisomer of dispiro-1,2,4,5-tetraoxanes on their binding mode with heme and on antimalarial activity: molecular docking studies.

Based on the fact that different isomers may exhibit substantial distinct activities, quantum chemical calculations and automated molecular docking simulations were carried out for 13 dispiro-1,2,4,5-tetraoxane compounds, which experimentally exist as a mixture of several isomers, to elucidate the most probable isomer(s) responsible for their antimalarial activity. The results indicate significant effects of stereoisomer on the binding mode and the activity. Moreover, the antimalarial potency of each compound can be described by the docking results. Compounds 1, 2, 4, 5, 7, and 9 have the most probable isomers coordinate suitably with heme iron and hence they have high activities while the most probable isomer in compounds 3 and 8 could not bind appropriately to heme yielding only moderate activities. On the other hand, the steric hindrance in compounds 11-13 prevents an approach of heme iron to peroxide bonds resulting in a devoid of antimalarial activity. However, compounds 6 and 10 with isopropyl substituents exhibit a different docking character, which is possibly caused by a limitation in molecular flexibility of the available docking technique. Our results can be used as a guideline for stereochemical control in synthesis process to improve drug's potency.

Animals↗

Benchmarking sets for molecular docking.

Ligand enrichment among top-ranking hits is a key metric of molecular docking. To avoid bias, decoys should resemble ligands physically, so that enrichment is not simply a separation of gross features, yet be chemically distinct from them, so that they are unlikely to be binders. We have assembled a directory of useful decoys (DUD), with 2950 ligands for 40 different targets. Every ligand has 36 decoy molecules that are physically similar but topologically distinct, leading to a database of 98,266 compounds. For most targets, enrichment was at least half a log better with uncorrected databases such as the MDDR than with DUD, evidence of bias in the former. These calculations also allowed 40x40 cross-docking, where the enrichments of each ligand set could be compared for all 40 targets, enabling a specificity metric for the docking screens. DUD is freely available online as a benchmarking set for docking at http://blaster.docking.org/dud/.

Adenosine Deaminase↗

A theoretical study to investigate D2DAR/D4DAR selectivity: receptor modeling and molecular docking of dopaminergic ligands.

Molecular modeling methods have been applied to construct three-dimensional models for dopaminergic ligand complexes with D2 and D4 receptor subtypes (D2DAR and D4DAR), using the bovine rhodopsin crystal structure as a template for the modeling study. Different dopaminergic ligands, in particular the N-n-propyl-substituted 3-aryl- and 3-cyclohexylpiperidines, were docked into the D2DAR and the D4DAR, to evaluate the agreement between theoretical and experimental results as regards their D2/D4 selectivity. The different position of an aromatic region in the two receptors might explain the structural basis of this biological property.

Amino Acid Sequence↗

Mechanism exploration of divergent partial denitrification performance under tetracycline stress: Insights from functional gene, electron transport and molecular docking.

Nitrates and antibiotics like tetracycline (TC) coexist in wastewater and inhibit nitrite (NO₂--N) accumulation during partial denitrification (PD), restricting anammox coupling. A moving bed biofilm reactor (PD-MBBR) and a sequencing batch reactor (PD-SBR) were compared under TC stress (0-8 mg/L). The PD-MBBR proved more robust, sustaining a high nitrate transformation ratio (NTR) of 95.11% and ∼53% TC removal. Metagenomic sequencing, quantitative polymerase chain reaction (qPCR), and molecular docking revealed this tolerance stemmed from physical shielding and metabolic compensation. Carrier-attached growth promoted extracellular polymeric substances (EPS) overproduction, forming a dense barrier preventing TC from binding to key denitrifying enzymes. The biofilm maintained stable nitrate reductase (NAR) activity via high narG and napA gene abundances, while nitrite reductase (NIR) was inhibited, ensuring efficient NO₂--N accumulation. This was supported by hyperactivated electron transport chain components, with complex III relative abundance increasing 15.08% and peak enzymatic activity reaching 149.02%. While IntI1-mediated horizontal gene transfer fortified community defense, concentrated antibiotic resistance genes (ARGs) within the biofilm pose a secondary dissemination risk. Thus, PD-MBBR provides an efficient pretreatment strategy for anammox, though downstream ARGs management is warranted.

Denitrification↗

Lead discovery using molecular docking.

As the structures of more and more proteins and nucleic acids become available, molecular docking is increasingly considered for lead discovery. Recent studies consider the hit-rate enhancement of docking screens and the accuracy of docking structure predictions. As more structures are determined experimentally, docking against homology-modeled targets also becomes possible for more proteins. With more docking studies being undertaken, the 'drug-likeness' and specificity of docking hits is also being examined.

Animals↗

Single-coordinate-driving method for molecular docking: application to modeling of guest inclusion in cyclodextrin.

An extension of the computer program CICADA has been developed that allows us to use the single-coordinate-driving (SCD) method for flexible molecular docking. The docking procedure is composed of three independent space rotations, three independent translations, and the torsions selected by the user. One of the coordinates is driven; the other coordinates are relaxed. This procedure follows low-energy wells on the potential energy surface of the entire system. The program allows us to dock more than one ligand molecule to the receptor. We ran two test examples, docking N,N-dimethylformamide into alpha-cyclodextrin and R-phenoxypropionic acid into beta-cyclodextrin. The test examples showed that the SCD approach is able to overcome high-energy barriers and to cover the entire box within which the search is performed. The limitations of molecular dynamics docking in comparison with our approach also are discussed. The philosophy of the newly developed approach is not only to find the best dock for the receptor-ligand(s) system, but also to describe all the important binding modes and provide a good starting point for studying the dynamics within the cavity during the docking process.

Carbohydrate Conformation↗

GEMDOCK: a generic evolutionary method for molecular docking.

We have developed an evolutionary approach for flexible ligand docking. This approval, GEMDOCK, uses a Generic Evolutionary Method for molecular DOCKing and an empirical scoring function. The former combines both discrete and continuous global search strategies with local search strategies to speed up convergence, whereas the latter results in rapid recognition of potential ligands. GEMDOCK was tested on a diverse data set of 100 protein-ligand complexes from the Protein Data Bank. In 79% of these complexes, the docked lowest energy ligand structures had root-mean-square derivations (RMSDs) below 2.0 A with respect to the corresponding crystal structures. The success rate increased to 85% if the structure water molecules were retained. We evaluated GEMDOCK on two cross-docking experiments in which each ligand of a protein ensemble was docked into each protein of the ensemble. Seventy-six percent of the docked structures had RMSDs below 2.0 A when the ligands were docked into foreign structures. We analyzed and validated GEMDOCK with respect to various search spaces and scoring functions, and found that if the scoring function was perfect, then the predicted accuracy was also essentially perfect. This study suggests that GEMDOCK is a useful tool for molecular recognition and may be used to systematically evaluate and thus improve scoring functions.

Algorithms↗

Matching chemistry and shape in molecular docking.

We have added a chemical filter to the ligand placement algorithm of the molecular docking program DOCK. DOCK places ligands in receptors using local shape features. Here we label these shape features by chemical type and insist on complementary matches. We find fewer physically unrealistic complexes without reducing the number of complexes resembling the known ligand-receptor configurations. Approximately 10-fold fewer complexes are calculated and the new algorithm is correspondingly 10-fold faster than the previous shape-only matching. We tested the new algorithm's ability to reproduce three known ligand-receptor complexes: methotrexate in dihydrofolate reductase, deoxyuridine monophosphate in thymidylate synthase and pancreatic trypsin inhibitor in trypsin. The program found configurations within 1 A of the crystallographic mode, with fewer non-native solutions compared with shape-only matching. We also tested the program's ability to retrieve known inhibitors of thymidylate synthase and dihydrofolate reductase by screening molecular databases against the enzyme structures. Both algorithms retrieved many known inhibitors preferentially to other compounds in the database. The chemical matching algorithm generally ranks known inhibitors better than does matching based on shape alone.

Algorithms↗

Orientational sampling and rigid-body minimization in molecular docking revisited: on-the-fly optimization and degeneracy removal.

Strategies for computational association of molecular components entail a compromise between configurational exploration and accurate evaluation. Following the work of Meng et al. [Proteins, 17 (1993) 266], we investigate issues related to sampling and optimization in molecular docking within the context of the DOCK program. An extensive analysis of diverse sampling conditions for six receptor-ligand complexes has enabled us to evaluate the tractability and utility of on-the-fly force-field score minimization, as well as the method for configurational exploration. We find that the sampling scheme in DOCK is extremely robust in its ability to produce configurations near to those experimentally observed. Furthermore, despite the heavy resource demands of refinement, the incorporation of a rigid-body, grid-based simplex minimizer directly into the docking process results in a docking strategy that is more efficient at retrieving experimentally observed configurations than docking in the absence of optimization. We investigate the capacity for further performance enhancement by implementing a degeneracy checking protocol aimed at circumventing redundant optimizations of geometrically similar orientations. Finally, we present methods that assist in the selection of sampling levels appropriate to desired result quality and available computational resources.

Algorithms↗

A structure-activity relationship study of catechol-O-methyltransferase inhibitors combining molecular docking and 3D QSAR methods.

A panel of 92 catechol-O-methyltransferase (COMT) inhibitors was used to examine the molecular interactions affecting their biological activity. COMT inhibitors are used as therapeutic agents in the treatment of Parkinson's disease, but there are limitations in the currently marketed compounds due to adverse side effects. This study combined molecular docking methods with three-dimensional structure-activity relationships (3D QSAR) to analyse possible interactions between COMT and its inhibitors, and to incite the design of new inhibitors. Comparative molecular field analysis (CoMFA) and GRID/GOLPE models were made by using bioactive conformations from docking experiments, which yielded q2 values of 0.594 and 0.636, respectively. The docking results, the COMT X-ray structure, and the 3D QSAR models are in agreement with each other. The models suggest that an interaction between the inhibitor's catechol oxygens and the Mg2+ ion in the COMT active site is important. Both hydrogen bonding with Lys144, Asn170 and Glu199, and hydrophobic contacts with Trp38, Pro174 and Leu198 influence inhibitor binding. Docking suggests that a large R1 substituent of the catechol ring can form hydrophobic contacts with side chains of Val173, Leu198, Met201 and Val203 on the COMT surface. Our models propose that increasing steric volume of e.g. the diethylamine tail of entacapone is favourable for COMT inhibitory activity.

Animals↗

New approach to molecular docking and its application to virtual screening of chemical databases

This paper describes the validation of a molecular docking method and its application to virtual database screening. The code flexibly docks ligand molecules into rigid receptor structures using a tabu search methodology driven by an empirically derived function for estimating the binding affinity of a protein-ligand complex. The docking method has been tested on 70 ligand-receptor complexes for which the experimental binding affinity and binding geometry are known. The lowest energy geometry produced by the docking protocol is within 2.0 A root mean square of the experimental binding mode for 79% of the complexes. The method has been applied to the problem of virtual database screening to identify known ligands for thrombin, factor Xa, and the estrogen receptor. A database of 10,000 randomly chosen "druglike" molecules has been docked into the three receptor structures. In each case known receptor ligands were included in the study. The results showed good separation between the predicted binding affinities of the known ligand set and the database subset.

Journal Article↗

Domesticated Argania spinosa in Eastern Morocco: HPLC-DAD/GC-MS Chemical Profiling, Antioxidant and Antidiabetic Activities, and Network Pharmacology-Guided Molecular Docking.

The argan tree (Argania spinosa) is an endemic Moroccan species known for its primary product, argan oil, which possesses exceptional nutritional and medicinal properties. The current study aimed to evaluate and compare the antidiabetic and antioxidant activities of argan oil obtained from the introduced and native argan tree in eastern Morocco, to analyze its chemical composition using HPLC-DAD and GC-MS, and to investigate the molecular mechanisms behind the obtained pharmacological activities through an in silico pharmacological networking and molecular docking study. The results revealed that argan oil from all three regions of Morocco (Oujda, Agadir, and Chouihya) is rich in oleic and linoleic acids as major constituents, along with the presence of significant tocopherols. Regarding the antioxidant assays, including DPPH radical scavenging and iron-reducing power tests, argan oil from Oujda exhibited the highest activity, with the lowest IC50 values of 15.25 ± 0.022 mg/mL and 28.5 ± 1.7 mg/mL, respectively. Concerning the antidiabetic activity, we found that oil from Chaouihya showed the strongest α-amylase inhibition, while Oujda oil had the highest antiglycation activity, indicating that even introduced argan trees retain potent bioactivity. The results of the in silico investigation suggested that tocopherols may contribute to the antioxidant and antidiabetic potential of argan oil, showing predicted antioxidant activity (Pa = 0.843-0.967) and favorable binding affinities toward iNOS (ΔG = -9.3 kcal mol-1) and α-glucosidase (ΔG = -8.2 kcal mol-1). The identified fatty acids also showed predicted insulin-promoting activity (Pa = 0.59-0.75) and moderate enzyme-binding potential. Pharmacological network analysis identified 51 shared genes associated with antioxidant, antidiabetic, and argan-related targets, with enrichment of the AGE-RAGE signaling pathway. These computational findings provide possible molecular associations that may help explain the observed biological activities, although they remain predictive and require experimental validation. Overall, the in silico analysis suggests that tocopherols could be among the contributors to the multi-target profile of Argania spinosa oil, while fatty acids may provide complementary effects related to glycemic regulation.

Sapotaceae↗

Probing molecular docking in a charged model binding site.

A model binding site was used to investigate charge-charge interactions in molecular docking. This simple site, a small (180A(3)) engineered cavity in cyctochrome c peroxidase (CCP), is negatively charged and completely buried from solvent, allowing us to explore the balance between electrostatic energy and ligand desolvation energy in a system where many of the common approximations in docking do not apply. A database with about 5300 molecules was docked into this cavity. Retrospective testing with known ligands and decoys showed that overall the balance between electrostatic interaction and desolvation energy was captured. More interesting were prospective docking scre"ens that looked for novel ligands, especially those that might reveal problems with the docking and energy methods. Based on screens of the 5300 compound database, both high-scoring and low-scoring molecules were acquired and tested for binding. Out of 16 new, high-scoring compounds tested, 15 were observed to bind. All of these were small heterocyclic cations. Binding constants were measured for a few of these, they ranged between 20microM and 60microM. Crystal structures were determined for ten of these ligands in complex with the protein. The observed ligand geometry corresponded closely to that predicted by docking. Several low-scoring alkyl amino cations were also tested and found to bind. The low docking score of these molecules owed to the relatively high charge density of the charged amino group and the corresponding high desolvation penalty. When the complex structures of those ligands were determined, a bound water molecule was observed interacting with the amino group and a backbone carbonyl group of the cavity. This water molecule mitigates the desolvation penalty and improves the interaction energy relative to that of the "naked" site used in the docking screen. Finally, six low-scoring neutral molecules were also tested, with a view to looking for false negative predictions. Whereas most of these did not bind, two did (phenol and 3-fluorocatechol). Crystal structures for these two ligands in complex with the cavity site suggest reasons for their binding. That these neutral molecules do, in fact bind, contradicts previous results in this site and, along with the alkyl amines, provides instructive false negatives that help identify weaknesses in our scoring functions. Several improvements of these are considered.

Amidines↗

Molecular docking reveals a novel binding site model for fentanyl at the mu-opioid receptor.

The ligand binding modes of a series of fentanyl derivatives are examined using a combination of conformational analysis and molecular docking to the mu-opioid receptor. Condensed-phase molecular dynamics simulations are applied to evaluate potential relationships between ligand conformation and fentanyl substitution and to generate probable "bioactive" structures for the ligand series. Automated docking of the largely populated solution conformers identified a common binding site orientation that places the N-phenethyl group of fentanyl deep in a crevice between transmembrane (TM) helices II and III while the N-phenylpropanamide group projected toward a pocket formed by TM-III, -VI, and -VII domains. An analysis of the binding modes indicates the most potent fentanyl derivatives adopt an extended conformation both in solution and in the bound state, suggesting binding affinity may depend on the conformational preferences of the ligands. The results are consistent with ligand binding data derived from chimeric and mutant receptor studies as well as structure-activity relationship data reported on a wide range of fentanyl analogues. The binding site model is also compared to that of N-phenethylnormorphine. An overlay of the bound conformation of the opiate and cis-3-methylfentanyl shows the N-phenethyl groups occupy equivalent binding domains in the receptor. While the cationic amines of both ligand classes were found docked to an established anchor site (D149 in TM-III), no overlap was observed between the N-phenylpropanamide group and the remaining components of the opiate scaffold. The unique binding mode(s) proposed for the fentanyl series may, in part, explain the difficulties encountered in defining models of recognition at the mu-receptor and suggest opioid receptors may display multiple binding epitopes. Furthermore, the results provide new insight to the design of experiments aimed at understanding the structural basis to the differential selectivities of ligands at the mu-, delta-, and kappa-opioid receptors.

Amino Acid Sequence↗

Precision UV-B irradiation for flavonoid biofortification in indoor-cultivated Morus nigra: Integrated multi-omics and molecular docking insights.

Precision application of UV irradiation represents an effective strategy for improving the quality of functional food crops under controlled environmental conditions. Morus nigra serves as a significant functional food resource in Xinjiang, with its leaves being rich in diverse bioactive compounds with nutritional and health-promoting properties. In this study, a low-dose UV-B treatment regimen was developed to enhance flavonoids in indoor-cultivated M. nigra without growth penalty. Multi-omics revealed a hormone shift (suppressed auxin vs. activated JA signaling). This triggered transcriptional reprogramming of PAL/4CL and CHS/CYP75B1, redirecting carbon flux towards flavonoid biosynthesis. Transient overexpression of MYB, NAC, and TIFY variants validated this regulatory network. The UV-B-induced NAC upregulated key genes and flavonoids, while different TIFY members showed diverse regulatory effects. Molecular docking predicted that the induced flavonoids had hypoglycemic, antioxidant, and anti-browning potential. This study indicates targeted UV-B as a green biofortification strategy for high-value crops in controlled environments.

Flavonoid biosynthesis↗