SEARCH · PubMed Health

Results for “Molecular Docking Simulation”

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 1,171 records · Page 65Linked to original sources

PubMed2002

On the molecular interaction between lactoferrin and the dye Red HE-3b. A novel approach for docking a charged and highly flexible molecule to protein surfaces.

Lactoferrin (Lf) is a non-heme, iron binding protein present in many physiological fluids of vertebrates where its main role is the microbicidal activity. It has been isolated by different methods, including dye-affinity chromatography. Red HE-3B is one of the most common triazinic dyes applied in protein purification, but scant knowledge is available on structural details and on the energetics of its interaction with proteins. In this work we present a computational approach useful for identifying possible binding sites for Red HE-3B in apo and holo forms of Lfs from human and bovine source. A new geometrical description of Red HE-3B is introduced which greatly simplifies the conformational analysis. This approach proved to be of particular advantage for addressing conformational ensembles of highly flexible molecules. Predictions from this analysis were correlated with experimentally observed dye-binding sites, as mapped by protection from proteolysis in Red HE-3B/Lf complexes. This method could bear relevance for the screening of possible dye-binding sites in proteins whose structure is known and as a potential tool for the design of engineered protein variants which could be purified by dye-affinity chromatography.

Amino Acid Sequence↗

PubMed1995

A three dimensional model for the interaction of MSH with the melanocortin-1 receptor.

A model for the interaction of the melanocortin-1 (MC1) receptor with MSH peptide has been constructed. The model was built by homology modelling using bacteriorhodopsin as template. A cyclic analogue of MSH could be docked into a binding pocket located between transmembrane (TM) domains 2, 3 and 6 of the receptor. The most significant receptor to ligand interactions occur between D117 in TM3 of receptor with histidine in cyclic MSH-peptide, H260 in TM6 with glutamic in peptide and D121 in TM3 with arginine in peptide. The model finds support from earlier mutagenesis data.

Amino Acid Sequence↗

PubMed2005

Three-dimensional structure of Plasmodium falciparum Ca2+ -ATPase(PfATP6) and docking of artemisinin derivatives to PfATP6.

Construction of the 3D structure of PfATP6 by homology modeling and docking simulation of artemisinin derivatives to this protein model are reported. Docking and consequent LUDI scores show good relation with in vitro antimalarial activities. The main binding source of artemisinins to the PfATP6 is hydrophobic interaction and biologically important peroxide bonds were exposed to outside of the binding pocket. This study suggests binding of artemisinin to PfATP6 precedes activation of peroxide bond by Fe(2+) species.

Amino Acid Sequence↗

PubMed2006

Imaging the migration pathways for O2, CO, NO, and Xe inside myoglobin.

Myoglobin (Mb) is perhaps the most studied protein, experimentally and theoretically. Despite the wealth of known details regarding the gas migration processes inside Mb, there exists no fully conclusive picture of these pathways. We address this deficiency by presenting a complete map of all the gas migration pathways inside Mb for small gas ligands (O2, NO, CO, and Xe). To accomplish this, we introduce a computational approach for studying gas migration, which we call implicit ligand sampling. Rather than simulating actual gas migration events, we infer the location of gas migration pathways based on a free-energy perturbation approach applied to simulations of Mb's dynamical fluctuations at equilibrium in the absence of ligand. The method provides complete three-dimensional maps of the potential of mean force of gas ligand placement anywhere inside a protein-solvent system. From such free-energy maps we identify each gas docking site, the pathways between these sites, to the heme and to the external solution. Our maps match previously known features of these pathways in Mb, but also point to the existence of additional exits from the protein matrix in regions that are not easily probed by experiment. We also compare the pathway maps of Mb for different gas ligands and for different animal species.

Carbon Monoxide↗

PubMedDate not supplied

Modeling ligand-receptor interaction for some MHC class II HLA-DR4 peptide mimetic inhibitors using several molecular docking and 3D QSAR techniques.

The ligand-receptor interaction between some peptidomimetic inhibitors and a class II MHC peptide presenting molecule, the HLA-DR4 receptor, was modeled using some three-dimensional (3D) quantitative structure-activity relationship (QSAR) methods such as the Comparative Molecular Field Analysis (CoMFA), Comparative Molecular Similarity Indices Analysis (CoMSIA), and a pharmacophore building method, the Catalyst program. The structures of these peptidomimetic inhibitors were generated theoretically, and the conformations used in the 3D QSAR studies were defined by docking them into the known structure of HLA-DR4 receptor through the GOLD, GLIDE Rigidly, GLIDE Flexible, and Xscore programs. Some of the parameters used in these docking programs were selected by docking an X-ray ligand into the receptor and comparing the root-means-square difference (RMSD) computed between the coordinates of the X-ray and docked structure. However, the goodness of a docking result for docking a series of peptidomimetic inhibitors into the HLA-DR4 receptor was judged by comparing the Spearman's rank correlation coefficient computed between each docking result and the activity data taken from the literature. The best CoMFA and CoMSIA models were constructed using the aligned structures of the best docking result. The CoMSIA was conducted in a stepwise manner to identify some important molecular features that were further employed in a pharmacophore building process by the Catalyst program. It was found that most inhibitors of the training set were accurately predicted by the best pharmacophore model, the Hypo1 hypothesis constructed. The deviation or conflict found between the actual and predicted activities of some inhibitors of both the training and the test sets were also investigated by mapping the Hypo1 hypothesis onto the corresponding structures of the inhibitors.

Binding Sites↗

PubMed2004

Saturation transfer difference NMR and computational modeling of a sialoadhesin-sialyl lactose complex.

The siglecs are a family of I-type lectins binding to sialic acids on the cell surface. Sialoadhesin (siglec-1) is expressed at much higher levels in inflammatory macrophages and specifically binds to alpha-2,3-sialylated N-acetyl lactosamine residues of glycan chains. The terminal disaccharide alpha-D-Neu5Ac-(2-->3)-beta-D-Gal is thought to be the main epitope recognized by sialoadhesin. To understand the basis of this biological recognition reaction we combined NMR experiments with a molecular modeling study. We employed saturation transfer difference (STD) NMR experiments to characterize the binding epitope of alpha-2,3-sialylated lactose, alpha-D-Neu5Ac-(2-->3)-beta-D-Gal-(1-->4)-D-Glc 1 to sialoadhesin at atomic resolution. The experimental results were compared to a computational docking model and to X-ray data of a complex of sialyl lactose and sialoadhesin. The data reveal that sialoadhesin mainly recognizes the N-acetyl neuraminic acid and a small part of the galactose moiety of 1. The crystal structure of a complex of sialoadhesin with sialyl lactose 1 was used as a basis for a modeling study using the FlexiDock algorithm. The model generated was very similar to the original crystal structure. Therefore, the X-ray data were used to predict theoretical STD values utilizing the CORCEMA-STD protocol. The good agreement between experimental and theoretical STD values indicates that a combined modeling/STD NMR approach yields a reliable structural model for the complex of sialoadhesin with alpha-D-Neu5Ac-(2-->3)-beta-D-Gal-(1-->4)-D-Glc 1 in aqueous solution.

Carbohydrate Sequence↗

PubMed2005

Comparison of automated docking programs as virtual screening tools.

The performance of several commercially available docking programs is compared in the context of virtual screening. Five different protein targets are used, each with several known ligands. The simulated screening deck comprised 1000 molecules from a cleansed version of the MDL drug data report and 49 known ligands. For many of the known ligands, crystal structures of the relevant protein-ligand complexes were available. We attempted to run experiments with each docking method that were as similar as possible. For a given docking method, hit rates were improved versus what would be expected for random selection for most protein targets. However, the ability to prioritize known ligands on the basis of docking poses that resemble known crystal structures is both method- and target-dependent.

Binding Sites↗

PubMed2001

Phalloidin synthetic analogues: structural requirements in the interaction with F-actin.

Synthetic derivatives of phalloidin have been investigated in solution by circular dichroism (CD) and NMR spectroscopy. They differ from natural phalloidin (PHD). bicyclo(Ala1-D-Thr2-Cys3-cis-4-hydroxy-Pro4-Ala5-2-mercapto-Trp6-(OH)2Leu7)(S-3 --> 6), in that they are modified at positions 2, 3, and 7. Among these synthetic analogues, structural differences and varying degrees of atropisomerism are found. By comparing the respective molecular models obtained by restrained molecular dynamics (RMD) simulations based on experimental NMR data, structural features that may be responsible for the different biological behavior become apparent. Our results indicate that the structural changes that result from an inversion of chirality of residue 3 lead to a complete loss of toxicity. Conversely, toxicity is less affected by the structural changes that stem from an inversion of chirality of residue 2. Moreover, unlike the other phallotoxins, when the thioether unit bridges to the opposite face of the main peptide ring, in contrast to the situation in other phallotoxins, large structural changes are observed as well as a total loss of activity. Molecular models of the synthetic phalloidin analogues have been used to investigate the necessary structural requirements for the interaction with F-actin. To this end, the F-actin/PHD model of M. Lorenz et al. was employed; docking experiments of our molecular models in the PHD binding site are presented.

Actins↗

PubMedDate not supplied

3-D QSAR studies on histone deacetylase inhibitors. A GOLPE/GRID approach on different series of compounds.

Docking simulation and three-dimensional quantitative structure-activity relationships (3D-QSARs) analyses were conducted on four series of HDAC inhibitors. The studies were performed using the GRID/GOLPE combination using structure-based alignment. Twelve 3-D QSAR models were derived and discussed. Compared to previous studies on similar inhibitors, the present 3-D QSAR investigation proved to be of higher statistical value, displaying for the best global model r2, q2, and cross-validated SDEP values of 0.94, 0.83, and 0.41, respectively. A comparison of the 3-D QSAR maps with the structural features of the binding site showed good correlation. The results of 3D-QSAR and docking studies validated each other and provided insight into the structural requirements for anti-HDAC activity. To our knowledge this is the first 3-D QSAR application on a broad molecular diversity training set of HDACIs.

Enzyme Inhibitors↗

PubMed2001

Ligand-receptor docking with the Mining Minima optimizer.

The optimizer developed for the Mining Minima algorithm, which uses ideas from Genetic Algorithms, the Global Underestimator Method, and Poling, has been adapted for use in ligand-receptor docking. The present study describes the resulting methodology and evaluates its accuracy and speed for 27 test systems. The performance of the new docking algorithm appears to be competitive with that of previously published methods. The energy model, an empirical force field with a distance-dependent dielectric treatment of solvation, is adequate for a number of test cases, although incorrect low-energy conformations begin to compete with the correct conformation for larger sampling volumes and for highly solvent-exposed binding sites that impose little steric constraint on the ligand.

Algorithms↗

PubMed2005

The orientations of cytochrome c in the highly dynamic complex with cytochrome b5 visualized by NMR and docking using HADDOCK.

The interaction of bovine microsomal ferricytochrome b5 with yeast iso-1-ferri and ferrocytochrome c has been investigated using heteronuclear NMR techniques. Chemical-shift perturbations for 1H and 15N nuclei of both cytochromes, arising from the interactions with the unlabeled partner proteins, were used for mapping the interacting surfaces on both proteins. The similarity of the binding shifts observed for oxidized and reduced cytochrome c indicates that the complex formation is not influenced by the oxidation state of the cytochrome c. Protein-protein docking simulations have been performed for the binary cytochrome b5-cytochrome c and ternary (cytochrome b5)-(cytochrome c)2 complexes using a novel HADDOCK approach. The docking procedure, which makes use of the experimental data to drive the docking, identified a range of orientations assumed by the proteins in the complex. It is demonstrated that cytochrome c uses a confined surface patch for interaction with a much more extensive surface area of cytochrome b5. Taken together, the experimental data suggest the presence of a dynamic ensemble of conformations assumed by the proteins in the complex.

Animals↗

PubMed1998

A computational system for modelling flexible protein-protein and protein-DNA docking.

A computational system is described that predicts the structure of protein/protein and protein/DNA complexes starting from unbound coordinate sets. The approach is (i) a global search with rigid-body docking for complexes with shape complementarity and favourable electrostatics; (ii) use of distance constraints from experimental (or predicted) knowledge of critical residues; (iii) use of pair potential to screen docked complexes and (iv) refinement and further screening by protein-side chain optimisation and interfacial energy minimisation. The system has been applied to model ten protein/protein and eight protein-repressor/DNA (steps i to iii only) complexes. In general a few complexes, one of which is close to the true structure, can be generated.

Algorithms↗

PubMed1991

Computer-aided active-site-directed modeling of the herpes simplex virus 1 and human thymidine kinase.

Thymidine kinase (TK), which is induced by Herpes Simplex Virus 1 (HSV1), plays a key role in the antiviral activity of guanine derivatives such as aciclovir (ACV). In contrast, ACV shows only low affinity to the corresponding host cell enzyme. In order to define the differences in substrate binding of the two enzymes on molecular level, models for the three-dimensional (3-D) structures of the active sites of HSV1-TK and human TK were developed. The reconstruction of the active sites of HSV1-TK and human TK were developed. The reconstruction of the active sites started from primary and secondary structure analysis of various kinases. The results were validated to homologous enzymes with known 3-D structures. The models predict that both enzymes consist of a central core beta-sheet structure, connected by loops and alpha-helices very similar to the overall structure of other nucleotide binding enzymes. The phosphate binding site is made up of a highly conserved glycine-rich loop at the N-terminus of the proteins and a conserved region at the C-terminus. The thymidine recognition site was found about 100 amino acids downstream from the phosphate binding loop. The differing substrate specificity of human and HSV1-TK can be explained by amino-acid substitutions in the homologous regions. To achieve a better understanding of the structure of the active site and how the thymidine kinase proteins interact with their substrates, the corresponding complexes of thymidine and dihydroxypropoxyguanine (DHPG) with HSV1 and human TK were built. For the docking of the guanine derivative, the X-ray structure of Elongation Factor Tu (EF-Tu), co-crystallized with guanosine diphosphate, was taken as reference. Fitting of thymidine into the active sites was done with respect to similar interactions found in thymidylate kinase. To complement the analysis of the 3-D structures of the two kinases and the substrate enzyme interactions, site-directed mutagenesis of the thymidine recognition site of HSV1-TK has been undertaken, changing Asp162 in the thymidine recognition site into Asn. First investigations reveal that the enzymatic activity of the mutant protein is destroyed.

Amino Acid Sequence↗

PubMedDate not supplied

Interactions of aryloxyphenoxypropionic acids with sensitive and resistant acetyl-coenzyme a carboxylase by homology modeling and molecular dynamic simulations.

Acetyl-coenzyme A carboxylase (ACCase) has been identified as one of the most important targets of herbicides. In the present study, we constructed homology models of the carboxyl-transferase (CT) domain of ACCase from sensitive and resistant foxtail and used these models as templates to study the molecular mechanism of herbicide resistance and stereochemistry-activity relationships of aryloxyphenoxypropionates (APPs). In the homology modeling structures, the dimer of the CT domain was formed by the side-to-side arrangement of the two monomers, in such a way that the N domain of one molecule is placed next to the C domain of the other. The dimeric association of sensitive foxtail CT was found to differ from that of resistant foxtail CT, and the spatial orientation of two key residues, Leu-695 and Ile-695, in these dimers also differed. The mutation of Ile to Leu may perturb the conformation of the dimeric interface, which may account for the molecular mechanism of herbicide resistance. Further docking analysis indicated that the binding model of high-active compounds is similar to that in the crystal structure of the enzyme-ligand complex. The different spatial orientations of ester groups of the isomers of APPs may explain the stereochemistry-activity relationship. Ser-698 formed a H-bonding interaction with all of the docked ligands, while Tyr-728 formed a pi-pi stacking interaction with some of the APPs. These findings may enhance our understanding of the molecular mechanism of herbicide resistance and stereochemistry-activity relationships, which may provide a new starting point for the identification of more potent inhibitors against both sensitive and resistant ACCase.

Acetyl-CoA Carboxylase↗

PubMed2004

Design and synthesis of Rho kinase inhibitors (I).

Several structurally unrelated scaffolds of the Rho kinase inhibitor were designed using pharmacophore information obtained from the results of a high-throughput screening and structural information from a homology model of Rho kinase. A docking simulation using the ligand-binding pocket of the Rho kinase model helped to comprehensively understand and to predict the structure-activity relationship of the inhibitors. This understanding was useful for developing new Rho kinase inhibitors of higher potency and selectivity. We identified several potent platforms for developing the Rho kinase inhibitors, namely, pyridine, 1H-indazole, isoquinoline, and phthalimide.

Amino Acid Sequence↗

PubMed2000

Synthesis and aldose reductase inhibitory activities of novel N-nitromethylsulfonanilide derivatives.

A novel series of 14 N-nitromethylsulfonanilide derivatives were synthesized and evaluated for their ability to inhibit recombinant aldose reductase. Computational docking simulations provided a good explanation for the observed structure-activity relationships. Kinetic analysis of (2-fluoro-5-methyl-N-methyl)-N-nitromethylsulfonanilide, 11, one of the most potent compounds in this series with an IC50 = 0.35 M, showed uncompetitive inhibition. Subsequent in vitro culture studies of rat lenses with 11 indicated that this series of aldose reductase inhibitors are effective in either preventing or retarding sugar cataract formation associated with diabetes.

Aldehyde Reductase↗

PubMed1997

Modelling of the binding site of the human m1 muscarinic receptor: experimental validation and refinement.

Our model of the human m1 muscarinic receptor has been refined on the basis of the recently published projection map of bovine rhodopsin. The refined model has a slightly different helix arrangement, which reveals the presence of an extra hydrophobic pocket located between helices 3, 4 and 5. The interaction of series of agonists and antagonists with the m1 muscarinic receptor has been studied experimentally by site-directed mutagenesis. In order to account for the observed results, three-dimensional models of m1 ligands docked in the target receptor are proposed. Qualitatively, the obtained models are in good agreement with the experimental observations. Agonists and partial agonists have a relatively small size. They can bind to the same region of the receptor using, however, different anchoring receptor residues. Antagonists are usually larger molecules, filling almost completely the same pocket as agonists. They can usually produce much stronger interactions with aromatic residues. Experimental data combined with molecular modelling studies highlight how subtle and diverse receptor-ligand interactions could be.

Amino Acid Sequence↗

PubMed2006

Allosteric interactions with muscarinic acetylcholine receptors: complex role of the conserved tryptophan M2422Trp in a critical cluster of amino acids for baseline affinity, subtype selectivity, and cooperativity.

In general, the M2 subtype of muscarinic acetylcholine receptors has the highest sensitivity for allosteric modulators and the M5 subtype the lowest. The M2/M5 selectivity of some structurally diverse allosteric agents is known to be completely explained by M2 177Tyr and M2 423Thr in receptors whose orthosteric site is occupied by the conventional ligand N-methylscopolamine (NMS). This study explored the role of the conserved M2 422Trp and the adjacent M2 423Thr in the binding of alkane-bisammonio type modulators, gallamine, and diallylcaracurine V. Experiments were performed with human M2 or M5 receptors or mutants thereof. It was found that M2 422Trp and M2 423Thr independently influenced allosteric agent binding. The presence of M2 423Thr may enhance the affinity of binding, depending on the allosteric agent, either directly or indirectly (by avoiding sterical hindrance through its M5 counterpart 478His). Replacement of M2 422Trp and of the corresponding M5 477Trp by alanine revealed a pronounced contribution of these epitopes to subtype independent baseline affinity in NMS-bound and NMS-free receptors for all agents except diallylcaracurine V. In a few instances, this tryptophan also influenced cooperativity and subtype selectivity. Docking simulations using a three-dimensional M2 receptor model revealed that the aromatic rings of M2 177Tyr and M2 422Trp, in a concerted action, might fix one of the aromatic moieties of alkane-bisammonio compounds between them. Thus, M2 422Trp and the spatially adjacent M2 177Tyr, as well as M2 423Thr, form a cluster of amino acids within the allosteric binding cleft that is pivotal for both M2/M5 subtype selectivity and baseline affinity of allosteric agents.

Allosteric Regulation↗